[wvn]

Hello all,

I am attempting to validate a mathematical model for a chemical reaction that will eventually involve some nucleation physics. I am implementing the Convection-Diffusion equation multiple times for a number of species. The mathematical model involves a source term for species A being proportional (equal to for now) to the concentrations of species B and C. I have implemented this as follows:

in "AEqn.H":

solve
(
 fvm::ddt(A)
 +fvm::div(phi,A)
 ==fvm::laplacian(D_A,A)
 +(C*B)
);

However, when I compile the solver and run a test case, the results suggest that the "+(C*B)" term is ignored. C and B have their own corresponding equations that need to implement similar terms.

Thanks!

[WildeCat]

Originally Posted by wvn

Hello all,

I am attempting to validate a mathematical model for a chemical reaction that will eventually involve some nucleation physics. I am implementing the Convection-Diffusion equation multiple times for a number of species. The mathematical model involves a source term for species A being proportional (equal to for now) to the concentrations of species B and C. I have implemented this as follows:

in "AEqn.H":

solve
(
 fvm::ddt(A)
 +fvm::div(phi,A)
 ==fvm::laplacian(D_A,A)
 +(C*B)
);

However, when I compile the solver and run a test case, the results suggest that the "+(C*B)" term is ignored. C and B have their own corresponding equations that need to implement similar terms.

Thanks!

You need to provide a better picture of what you did for anyone to be able to help you

[wvn]

Fortunately, this issue has been fixed. Thanks!