|
[Sponsors] |
How should I implement this scalar transport equation |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
April 14, 2014, 09:38 |
How should I implement this scalar transport equation
|
#1 |
Senior Member
Dongyue Li
Join Date: Jun 2012
Location: Beijing, China
Posts: 838
Rep Power: 17 |
Hi guys,
I wanna solve this equation: I code it like this: Code:
fvScalarMatrix nEqn ( fvm::ddt(j) + fvm::div(phi1,j) == C1*j* epsilon -fvm::Sp(C2*pow(j, 2.0/9.0)* epsilon, j) ); Last edited by sharonyue; April 14, 2014 at 21:55. |
|
April 14, 2014, 21:32 |
|
#2 |
Member
xuhe-openfoam
Join Date: Aug 2013
Location: DaLian,china
Posts: 82
Rep Power: 12 |
hello!
maybe you could give more information , for example: C1>0 or C1<0, epsilon? |
|
April 14, 2014, 21:37 |
|
#3 |
Senior Member
Dongyue Li
Join Date: Jun 2012
Location: Beijing, China
Posts: 838
Rep Power: 17 |
||
April 14, 2014, 22:53 |
|
#4 |
Member
xuhe-openfoam
Join Date: Aug 2013
Location: DaLian,china
Posts: 82
Rep Power: 12 |
I suppose that :
1 "j+jSmall" is "j" 2 there is no limit of the value of "j" then , modify the source term " -C2*pow((j+jSmall),11.0/9.0)* epsilon " to : Code:
- 11.0/9.0*C2*fvm::Sp(epsilon,j) + 2.0/9.0*C2*epsilon*pow(j,11.0/9.0) of cource maybe there is some other reason for wired result, not only the discretization of source term. |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
how to make a scalar transport equation to be conserved? | immortality | OpenFOAM Programming & Development | 3 | May 5, 2013 08:11 |
Transport Equation for a scalar - Compressible Solver | alessio.nz | OpenFOAM Programming & Development | 5 | April 11, 2013 09:04 |
new scalar transport eqn. with new density | Tobi | OpenFOAM Programming & Development | 0 | September 3, 2012 08:18 |
Implement species transport equation | Tobi | OpenFOAM Programming & Development | 0 | June 2, 2012 13:26 |
how solve a scalar equation with Fluent | tomik | FLUENT | 0 | January 5, 2006 06:38 |