CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM

Transport of active scalars wihtin Operator Splitting Approach

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Display Modes
Old   October 4, 2011, 11:02
Default Transport of active scalars wihtin Operator Splitting Approach
New Member
Join Date: Mar 2010
Posts: 20
Rep Power: 9
pbe_cfd is on a distinguished road

I need to transport n (>>1) active scalar within Oper. Split. Approach. I have a code to solve the reaction step so, what I do is first solving the reaction step for cell based values then solving the transport step. A pseudo code is given here

1. Step: Solve reaction for each cell

for (i=0; i<C[0].size(); i++) // C := PtrList<volScalarField>
// Assign values of all the scalars at the i'th cell to Ci
for ( j=0; j<n; j++)
Ci[j]=C[j].internalField()[i]; Ci: is an array of length n, for the values of n Scalar at the i'th cell
solveReaction(Ci) // Ci: is an in/output and updated in the routine solveReaction
for ( j=0; j<n; j++)
C[j].internalField()[i]=Ci[j]=; Updating the values of C at the i'th cell with the updated values of Ci


2. Step: Transport each scalar field, C[i]
for (i=0; i<n; i++)
+ fvm::div(phi, C[i])
- fvm::laplacian(nuT, C[i]),

I implemented this pseudo code by extending scalarTranportFoam. The problem is, although each step works fine independently:
1) Cancel 2. step and update the scalar fields by reaction. Result: OK
2) Cancel 1. step and solve the passive transport problem for n Scalar. Result: OK

When I calculate them together, the results are as if I had not performed the reaction step. What can be the problem??? And, how can it be fixed ???

Thanks in advance,

PS. I can post or can send the piece of code via e-mail to whom are interested in.

PPS. Boundary values are another problem which is hacked by using symmetryPlane BC for the moment.
pbe_cfd is offline   Reply With Quote


Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
How to apply Transport Equation on DPM Scalars? yuhai FLUENT 0 May 19, 2010 10:38
PROBLEM: Scalar mass fraction Lorenzo Siemens 1 July 14, 2008 03:50

All times are GMT -4. The time now is 18:56.