coupled source term in 3 transport equations
I want to revise icoFoam to include three transport equations solving for C1, C2, and C3. the problem is that they are coupled in their source term. I have included "CEqn.H" in my transportFoam (revised icoFoam).
I have a constant depending on C3 (third species concentration) which is used in C1, and C2. There are also source term in C2, and C3 equations depending on C1. I faced some questions regarding this:
1st question: Is it correct to use "fvm:Sp(kpa,C1)" in transport equation of C2? or I can simply use "...+ kpa*C1" so it is accounted explicitly? (in the equations kpa and landJ are constants which are calculated and given in transportProperties, respectively)
2nd question: shall I have a loop to converge (solving equations sequentially)? If I do, I have to calculate the difference between each iteration results. Shall I declare a new variable to save each iteration results so I can have comparison in each iteration.
3rd Question:Is it correct to declare and initialize some constants (e.g. omega, kpa in the code below) in the "CEqn.H"? they are to be calculated after C3 equation is solved.
double kpa; // constant calculated after C3 is calculated
double omega; //constant calculated after C3 is calculated
double C3Crit=2; //a constant
I would appreciate your help. thanks in forward.
|All times are GMT -4. The time now is 17:34.|