DateAuthorLogOptionsExposure
2013-12-05Randall BrittenAdded ?SaveDefinitions -> True? to ?Manipulate? blocks so that CDF export works. [files][create_exposure]
2013-12-05Randall BrittenSeem to have successfully resolved all units issues now, and also worked around the quirks in the Mathematica 3D plotting functions. Works quite nicely now. [files][create_exposure]
2013-12-05Randall BrittenAttempting to clean up issues with units and unit conversions by more disciplined use of units suffixes in variable names. Some progress made so far. Still having issues showing trajectory relative to 3D plot. [files][create_exposure]
2013-12-04Randall BrittenWorking nicely now. Shows contour plot, 3D plot and individual time series of each reactant/product (i.e. reaction species), as well as final value. [files][create_exposure]
2013-12-04Randall BrittenCool! Got parametric plot and cursor interaction combined with contour plot. [files][create_exposure]
2013-12-04Randall BrittenGot cursor to determine initial conditions :-) [files][create_exposure]
2013-12-04Randall BrittenAttempts to get parametric plot of trajectory of point in {H, HCO3, CO2} space. Partially working. Haven?t succeeded yet in setting initial conditions based on cursor position. [files][create_exposure]
2013-12-02Randall BrittenExtracted the ODE numerical solve call to make Manipulate block cleaner. Also, showing final values of pH and HCO3. [files][create_exposure]
2013-12-01Randall BrittenAllow time range of ODE solution to be adjusted. [files][create_exposure]
2013-12-01Randall BrittenManipulate now shows ODE results allowing for initial pH and initial CO2 level to be adjusted. [files][create_exposure]
2013-12-01Randall BrittenAdded 3D plot. Also, made ODE numerical solve parameterised by initial conditions, which is work towards being able to pick IV?s from Manipulate UI and see solution interactively. [files][create_exposure]
2013-09-13Randall BrittenInitial value for HCO3 is now calculated from given initial pH and initial CO2 partial pressure, in the same way as in the Mathematica document. The system now remains at equilibrium :-) Also, using kf = K kb, with kb=1 and K calculated from given pK = 6.1, which is apparently the accepted value in the textbooks. Also, now actually using the initial values for CO2 concentration that was calculated from given partial pressure. [files][create_exposure]
2013-09-13Randall BrittenSimplified calculation of initial HCO3. [files][create_exposure]
2013-09-13Randall BrittenFinally, got ODE to match Davenport plot nicely. Using kb = 1, and deriving kf from K derived from textbook pK = 6.1, then deriving initialHCO3. [files][create_exposure]
2013-09-13Randall BrittenUpdated Mathematica to explore what equilibrium the textbook value of pK=6.1 gives compared to using the initial values read off Peter's slides. Change pK = pKbyInitial to pK = pKfromTextbook to see the difference: textbook value does not match initial conditions on slide being at equilibrium, and ODE runs to a pH of around 4.5! [files][create_exposure]
2013-09-12Randall BrittenCorrected conversion between partial pressure and concentration: based on solubility, not ideal gas law since we are modelling the solution. Seems obvious. [files][create_exposure]
2013-09-12Randall BrittenMathematica Davenport plot using contours of pCO2 as a function of pH and [HCO3-]. [files][create_exposure]
2013-09-12Randall BrittenBackout changeset f90bdc4dcf7f06125137b0272f5cb6a7d4b973ad [files][create_exposure]
2013-09-12Randall BrittenMathematica Davenport plot using contours of pCO2 as a function of pH and [HCO3-]. [files][create_exposure]
2013-09-11Randall BrittenTurns out that you have to use the reciprocal as the multiplier! [files][create_exposure]
2013-09-11Randall BrittenAlthough this runs, it is not producing numerically sensible results yet. Current state is partly a debugging state, and so the PhysicalConstants.cellml model isn't being used yet. [files][create_exposure]
2013-09-11Randall BrittenFirst attempt at introducing partial pressure as units for CO2 concentration. Still seem to be some problems with units and scaling, and a minor problem with encapsulation that should be easy to fix. [files][create_exposure]
2013-09-10Randall BrittenIncluding pH in model, setting initial pH to 7.4, which I think is physiologically realistic. Other parameters still seem far from physiological though. [files][create_exposure]
2013-07-15Randall BrittenAdded an H+ ion source (i.e. a proton source). [files][create_exposure]
2013-07-15Randall BrittenAdded a CO2 (Carbon dioxide) source and sync term. [files][create_exposure]
2013-07-15Randall BrittenOriginal versions from Peter Hunter. [files][create_exposure]