Changes between Version 5 and Version 6 of DianaICSolver

Timestamp:

Sep 12, 2011, 2:36:16 PM (14 years ago)

Author:

miroshkin

Comment:

--

DianaICSolver

@@ -1 @@
-= Examples =
-== Models requiring ICSolver ==
+= Description =
+To simulate chemical or biological process in Diana the Initial Value Problem should be solved at first. A guess for initial conditions can be set in model file. [[BR]] Standard {{{ IDACalcIC() }}} method of Sundials package can calculate consistent initial conditions for homogeneous differential algebraic equations systems of the type
 
-The following example shows how ICSolver can be used to find
-consistent initial conditions where just applying the standard
-IDA solver throws an exception:
+{{{
+#!html
+<table border="0" width="50%" align="center">
+<tr>
+  <td align = "center" fontsize ="50%" width="90%">
+  dx<sub>1</sub>/dt = f<sub>1</sub>(x<sub>1</sub>, x<sub>2</sub>, ..., x<sub>n</sub>, t)<br>
+  dx<sub>2</sub>/dt = f<sub>2</sub>(x<sub>1</sub>, x<sub>2</sub>, ..., x<sub>n</sub>, t)<br>
+             ...            <br>
+  dx<sub>n</sub>/dt = f<sub>n</sub>(x<sub>1</sub>, x<sub>2</sub>, ..., x<sub>n</sub>, t)
+  </td>
+  <td align="center">(1)
+  </td>
+</tr>
+</table>
+}}}
 
- * Adsorption column[[BR]]
-   [attachment:AdsColumn.py "script"] [attachment:AdsColumn.tar "script + model + results"]
+But in case of non-homogeneous differential algebraic equations systems of the type
+
+{{{
+#!html
+<table border="0" width="50%" align="center">
+<tr>
+  <td align = "center" fontsize ="50%" width="90%">
+  F<sub>1</sub>(dx<sub>1</sub>/dt, dx<sub>2</sub>/dt, ..., dx<sub>n</sub>/dt) = f<sub>1</sub>(x<sub>1</sub>, x<sub>2</sub>, ..., x<sub>n</sub>, t)<br>
+  F<sub>2</sub>(dx<sub>1</sub>/dt, dx<sub>2</sub>/dt, ..., dx<sub>n</sub>/dt) = f<sub>2</sub>(x<sub>1</sub>, x<sub>2</sub>, ..., x<sub>n</sub>, t)<br>
+             ...            <br>
+  F<sub>n</sub>(dx<sub>1</sub>/dt, dx<sub>2</sub>/dt, ..., dx<sub>n</sub>/dt) = f<sub>n</sub>(x<sub>1</sub>, x<sub>2</sub>, ..., x<sub>n</sub>, t)
+  </td>
+  <td align="center">(2)
+  </td>
+</tr>
+</table>
+}}}
+
+{{{IDACalcIC()}}} generates the exception
+{{{
+   ECapeUnknown in IDASolver (Diana) startIntegrator, IDACalcIC failed
+}}}
+
+Consistent initial conditions for differential algebraic equations system  of type (2) can be found with {{{IDACalcIC()}}} method, but the system should be manually changed to form (1) before.[[BR]]
+'''''ICSolver''''' can calculate initial conditions for systems of the type (2) without any conversions of the system.
+
+= Example =
+
+Adsorption column model contains the equations system of type (2) example[[BR]]
+Python script with calling standard IDACalcIC() method generates the exception[[BR]]
+Python script with using ICSolver and nonzero output verbosity level shows the initial state vector and message
+
+{{{
+   ICSolver: Consistent initial conditions have been found
+}}}
+
+== Additional information ==
+
+{{{ eso.GetAllVariables() }}} returns state vector values of model[[BR]]
+{{{ eso.GetAllResiduals() }}} returns vector of right hand side values of equations system[[BR]]
+With correct initial values of state vector for model equations all residuals should be equal to zeros.
 
 == Additional examples ==
@@ -21 +72 @@
  * MATLAB (with comments)[[BR]]
    [attachment:matlab.tar "files"]
-