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

compressible scalar transport euqation

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 22, 2010, 11:35
Default compressible scalar transport euqation
  #1
Senior Member
 
Join Date: Jan 2010
Location: Stuttgart
Posts: 150
Rep Power: 16
Chrisi1984 is on a distinguished road
Hi all,

I added three scalar transport equations to my compressible solver like this:

Quote:
{
fvScalarMatrix w4Eqn
(
fvm::ddt(rho, w4)
+ fvm::div(phi, w4)
- fvm::laplacian(turbulence->muEff(), w4)
);

w4Eqn.solve();


Info<< "Phase 1 mass fraction = "
<< w4.weightedAverage(mesh.V()).value()
<< " Min(w4) = " << min(w4).value()
<< " Max(w4) = " << max(w4).value()
<< endl;
}
My species are all the same material (air), so I thought I need not change something in the energy equation hEqn. I only want to track air flow through different inlets in my flow.

As long as my temperatures don't differ very much in my flow, I am getting proper results. But when I get regions with higher or lower temperatures, then the sum of my species exceeds 1 in the hot regions, where rho is small. In the colder regions, where rho is high, the sum of my species is too small.

So there might be a problem with the coupling of the transport equations with the energy equation.

I know that in other Solvers like rhoReactingFoam there is applied a multivariant scheme (mvConvection) to couple the energy equation with the species transport. But I thought if my species are all the same this is unnecessary.

Can anyone help me to solve my problem or has anyone any experience with similar problems?

Thanks in advance!


Best regards Chrisi
Chrisi1984 is offline   Reply With Quote

Old   October 18, 2010, 06:17
Default Solved
  #2
Senior Member
 
Join Date: Jan 2010
Location: Stuttgart
Posts: 150
Rep Power: 16
Chrisi1984 is on a distinguished road
Hi,

perhabs someone will once have a similar problem.

Here is the solution for my problem:

I changed the interpolation schemes and than everything was all right:

Quote:
interpolationSchemes
{
default linear;
div(U,p) upwind phi;
div(U,rho) upwind rho;
div(U,h) upwind h;
}
Regards Chrisi
Chrisi1984 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
Question about scalar transport osimonsimon OpenFOAM Running, Solving & CFD 33 August 21, 2020 11:45
dieselFoam problem!! trying to introduce a new heat transfer model vivek070176 OpenFOAM Programming & Development 10 December 24, 2014 00:48
Is it possible to use divergence in the source term of an scalar transport equation? jannnesss CFX 0 January 8, 2010 20:53
Unsteady kE with scalar transport - Outlet-geometry heavely effects scalar transport christian_andersen OpenFOAM 0 June 25, 2009 04:00


All times are GMT -4. The time now is 07:23.