CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM

Transport of active scalars wihtin Operator Splitting Approach

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 4, 2011, 12:02
Default Transport of active scalars wihtin Operator Splitting Approach
  #1
New Member
 
Evren
Join Date: Mar 2010
Posts: 20
Rep Power: 16
pbe_cfd is on a distinguished road
Hi,

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
}

end

2. Step: Transport each scalar field, C[i]
for (i=0; i<n; i++)
{
solve
(
fvm::ddt(C[i])
+ fvm::div(phi, C[i])
- fvm::laplacian(nuT, C[i]),
mesh.solver("class")
);
}

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,
Evren

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

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
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 11:38
PROBLEM: Scalar mass fraction Lorenzo Siemens 1 July 14, 2008 04:50


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