plotExpression(apply(-log10(CarbonDioxideInWater.H3O.solute)), false, "pH", 1);
Please note, that OH- (and CO3^-2) can be neglected from electroneutrality calculation, because of very small concentrations (in physiological pH) anyway.
And if SID>0 then also H3O+ can be also neglected from electroneutrality, because only bicarbonate anions HCO3- (or CO3^-2) are needed there to balance the electroneutrality.
-
The partial pressure of CO2 in gas are input parameter. Outputs are an amount of free disolved CO2 in liquid and an amount of HCO3-.
The partial pressure of CO2 in gas are input parameter. Outputs are an amount of free dissolved CO2 in liquid and an amount of HCO3-.
The titration slope der(pH)/der(SID)=17.5 1/(mol/L) at pH=7.4 and pCO2=40 mmHg.
Molar heat of formation (aqueous):
CO2: -413.5 kJ/mol (gas: -393.5 kJ/mol )
@@ -3291,7 +3291,7 @@ package Chemical "Domain with Molar Concentration and Molar Flow" smooth=Smooth.None)); annotation ( Documentation(info="Henderson-Hasselbalch equation in ideal buffered solution, where pH remains constant.
-The partial pressure of CO2 in gas are input parameter. Outputs are an amount of free disolved CO2 in liquid and an amount of HCO3-.
+The partial pressure of CO2 in gas are input parameter. Outputs are an amount of free dissolved CO2 in liquid and an amount of HCO3-.
", revisions="2014
Marek Matejak, Charles University, Prague, Czech Republic
@@ -3868,7 +3868,7 @@ package Chemical "Domain with Molar Concentration and Molar Flow"TK
base temperature of tabulated coeficients
base temperature of tabulated coefficients
ω
TK
base temperature of tabulated coeficients
base temperature of tabulated coefficients
Ω
This block identifies one specific chemical form of one macromolecule defined by forms of its subunits (one chosen chemical species - called specie).
Only main connector called species is designed for inflow and outflow of macromolecule to/from system. The concentration in this connector is the concentration of its specific specie.
Connectors subunitSpecies[:] represent specific forms of the macromolecule subunit types. If the subnunit type occures n-times in macromolecule, the inflow is n-time greater than the inflow of macromolecule.
-
Initial total concentrations of subunits must be set to be right distribution of total macromolecule concentration. So the ratios between subunit concentrations are the ratios of their occurence in macromolecule. In equilibrium are this proporties fullfiled.
Initial total concentrations of subunits must be set to be right distribution of total macromolecule concentration. So the ratios between subunit concentrations are the ratios of their occurrence in macromolecule. In equilibrium are this proporties fulfilled.
For example: If the macromolecule has four identical independent subunits and each subunit can occur in two form F1 and F2, then the concentration of macromolecule specie composed only from four subunits in form F1 is species.conc=conc*fF1^4.
Where:
conc is totat concentration of macromolecule in system accumulated by species.q,
@@ -5104,7 +5104,7 @@ Connector with one flow signal of type Real. end ChemicalPort_b; partial model OnePort - "Partial transfer of solute beween two ports without its accumulation" + "Partial transfer of solute between two ports without its accumulation" ChemicalPort_b q_out annotation (Placement(transformation(extent={{90,-10},{110,10}}))); diff --git a/Physiolibrary/Hydraulic.mo b/Physiolibrary/Hydraulic.mo index ad5b01c..3e4b82c 100644 --- a/Physiolibrary/Hydraulic.mo +++ b/Physiolibrary/Hydraulic.mo @@ -511,7 +511,7 @@ package Hydraulic "Domain with Pressure and Volumetric Flow" Physiolibrary.Types.RealIO.VolumeFlowRateOutput volumeflowrate annotation (Placement(transformation(extent={{80,-10},{100,10}}), iconTransformation(extent={{80,-10},{100,10}}))); - discrete Time T0 "begining of cardiac cycle"; + discrete Time T0 "beginning of cardiac cycle"; Boolean b(start=false); discrete Time HP "duration of cardiac cycle"; parameter Frequency HR=1.2; @@ -522,7 +522,7 @@ package Hydraulic "Domain with Pressure and Volumetric Flow" equation b = time - pre(T0) >= pre(HP) "true if new cardiac cycle begins"; when {initial(),b} then - T0 = time "set begining of cardiac cycle"; + T0 = time "set beginning of cardiac cycle"; HP = 1/HR "update length of carciac cycle"; TD2 = TD1 + (2/5)*HP "compute end time of systole"; end when; diff --git a/Physiolibrary/Population.mo b/Physiolibrary/Population.mo index f09058f..01f6bcd 100644 --- a/Physiolibrary/Population.mo +++ b/Physiolibrary/Population.mo @@ -98,7 +98,7 @@ package Population annotation ( HideResult=true, Dialog(group="Initialization")); parameter Types.Population NominalPopulation = 1 - "Numerical scale. Default is 1, but for huge amount of cells it should be any number in the apropriate numerical order of typical amount." + "Numerical scale. Default is 1, but for huge amount of cells it should be any number in the appropriate numerical order of typical amount." annotation ( HideResult=true, Dialog(tab="Solver",group="Numerical support of very huge populations")); parameter Types.PopulationChange NominalPopulationChange = 1/(60*60*24) "Numerical scale. Default change is 1 individual per day, but for much faster or much slower chnages should be different." @@ -382,7 +382,7 @@ Connector with one flow signal of type Real. end PopulationPort_b; partial model OnePort - "Partial change of population beween two ports without its accumulation" + "Partial change of population between two ports without its accumulation" PopulationPort_b port_b annotation (Placement(transformation(extent={{90,-10},{110,10}}))); diff --git a/Physiolibrary/Resources/Install/SystemModeler/README.txt b/Physiolibrary/Resources/Install/SystemModeler/README.txt index a464089..5159380 100644 --- a/Physiolibrary/Resources/Install/SystemModeler/README.txt +++ b/Physiolibrary/Resources/Install/SystemModeler/README.txt @@ -1,4 +1,4 @@ -Be carefull, this installation is not tested! It should copy the Physiolibrary files into ProgramFiles/Wolfram/SystemModeler directory. +Be careful, this installation is not tested! It should copy the Physiolibrary files into ProgramFiles/Wolfram/SystemModeler directory. To add library WITHOUT display units to Wolfram SystemModeler 4.0 please run "install.bat". diff --git a/Physiolibrary/SteadyStates.mo b/Physiolibrary/SteadyStates.mo index 7c66e95..92b52b5 100644 --- a/Physiolibrary/SteadyStates.mo +++ b/Physiolibrary/SteadyStates.mo @@ -1316,7 +1316,7 @@ package SteadyStates "Dynamic Simulation / Steady State" if Simulation <> Types.SimulationType.SteadyState then der(state) = change; - elseif not isDependent then /*** this test and equation exclusion could be done automatically, if the solver will be so smart that it removes all this dependend equations from the total equilibrated system. The most probable form of this dependent equation in equilibrium setting is (0 = 0). ***/ + elseif not isDependent then /*** this test and equation exclusion could be done automatically, if the solver will be so smart that it removes all this dependent equations from the total equilibrated system. The most probable form of this dependent equation in equilibrium setting is (0 = 0). ***/ change = 0; end if; @@ -1419,7 +1419,7 @@ package SteadyStates "Dynamic Simulation / Steady State" der(state) = change; else for i in 1:n loop - if not isDependent[n] then /*** this test and equation exclusion could be done automatically, if the solver will be so smart that it removes all this dependend equations from the total equilibrated system. The most probable form of this dependent equation in equilibrium setting is (0 = 0). ***/ + if not isDependent[n] then /*** this test and equation exclusion could be done automatically, if the solver will be so smart that it removes all this dependent equations from the total equilibrated system. The most probable form of this dependent equation in equilibrium setting is (0 = 0). ***/ change[i] = 0; end if; end for; diff --git a/Physiolibrary/Thermal.mo b/Physiolibrary/Thermal.mo index 1c18065..23317ea 100644 --- a/Physiolibrary/Thermal.mo +++ b/Physiolibrary/Thermal.mo @@ -642,7 +642,7 @@ i.e., it defines a fixed temperature as a boundary condition. transformation(extent={{90,-10},{110,10}}), iconTransformation( extent={{90,-10},{110,10}}))); parameter Types.SpecificEnergy dH=0 - "Enthalpy of incomming substance (i.e. enthalpy of solvation)"; // or 2428344 for water vaporization + "Enthalpy of incoming substance (i.e. enthalpy of solvation)"; // or 2428344 for water vaporization parameter Types.SpecificHeatCapacity SpecificHeat=4186.8 "Of inflowing medium"; //default heat capacity of water is 1 kcal/(degC.kg) diff --git a/Physiolibrary/Types.mo b/Physiolibrary/Types.mo index 338bceb..0f1b6c1 100644 --- a/Physiolibrary/Types.mo +++ b/Physiolibrary/Types.mo @@ -34,7 +34,7 @@ package Types "Physiological units with nominals"Blocks
Usefull blocks, that are missing in package Modelica.Blocks (MSL 3.2), cubic interpolation curves, multiplication factors.
Useful blocks, that are missing in package Modelica.Blocks (MSL 3.2), cubic interpolation curves, multiplication factors.
Chemical
Usefull blocks, that are missing in package Modelica.Blocks (MSL 3.2), cubic interpolation curves, multiplication factors.
Useful blocks, that are missing in package Modelica.Blocks (MSL 3.2), cubic interpolation curves, multiplication factors.