solver for subsonic compressible turbulent flow in OF 1.7
 

Hi all,

I want to solve a susonic, turbulent compressible flow in OF 1.7. Which solver should i use??? Please help..its argent

Depending upon the case (steady/unsteady), take a look at:

• rhoSimpleFoam (steady)
• rhoPisoFoam / rhoPimpleFoam (unsteady / pseudo-transient)

Best,

Thank you Alberto for the reply.

My doubt it can i use these solvers for the flow which has Mach numbers upto 1..??? Are the pressure based solvers able to handle higher subsonic mach numbers???

You asked about subsonic flows. Not highly supersonic. :-)

There are other solvers able to do that (see rhoCentralFoam and the various "sonic" solvers).

Best,

With all due respect Sir, i said "highly subsonic" means mach number arnd 0.8 or 0.9. Sorry if my statement was ambiguous.

Thanks for your quick replies. One more question, my friend suggested me that "rhoturbFoam" is the best choice for my problem but its not there in OF 1.7.. So what is the difference between "rhoturbForm" and say "rhosonicForm"??

I'm really sorry: I misread. Result of answering late in the evening ;-)

rhoTurbFoam was a RANS solver in OpenFOAM 1.5, and was merged in rhoPisoFoam in OpenFOAM 1.6 and 1.7, with RANS and LES capabilities.

rhoSonicFoam is a density-based inviscid solver.

About your original question, yes, in theory pressure-based solver should work with those values of Mach number (See for example Peric's papers and book on the topic).

Best,

Thank you very much Sir.

Your lucid answer really helped me a lot. This lack of documentation thing really makes it a difficult task for a newbie like me to get going with OF.

You are welcome. If you need further help, please don't hesitate to ask. The forum is a often a quick way to get answers.

Best,

Hello Alberto Sir,

As you have suggested, i am using rhopisoFoam for my problem. :)

Brief description of my problem:

Its cold flow simulation of mixing in a secondary combustor of a rocket motor. In this air and fuel streams are getting mixed together. I am planning to use RNGkepsilon model for this using rhoPisoFoam.

BCs:

Air inlet: Mach 0.3, stgT = 840K, staticT = 825K, Vin = 172m/s, massflow = 5.19kg/s, k and eps can be calculated using formulae

Fuel Inlet : Mach = 1, staticT = 650 >>Vin 244m/s, massflow = 0.2581kg/s, k and eps can be calculated using formulae

Wall: Adiabetic

Outlet: there is a chocked nozzle at the end of combustor which i am not simulating, but from chocked condition the pressure at the inlet of nozzle can be calculated as 7.2 bar.


I saw the tut for rhoPisoFoam. It uses SST kw model..but the "0" directory contains files for alphat, p, t, u, k, eps, omega, R. I got why we need u,p,t,k,omega. but why other files are present? I know what eps is, alphat is probably diffusivity but i dont know what is R.

Also i am confused that, out of the known BC parameters which should be specified???
Also i found out some of the BCs as:

what are these types??


Can you give me any direction???

Just Alberto ;)

R is the Reynolds stress tensor.

alphat is the turbulent thermal diffusivity, and it is used essentially in wall-functions.

The 'type' keyword specifies the type of the BC you are setting. You find a summary in section 5.2 of the user guide, while wall functions are briefly discussed in section 7.2.2.

Best,

Alberto,

Thanks for the reply.:)

What i am confused about is, do i really need to specify the values for omega, mut, alphat and R??? if yes how should i obtain them based on the BCs i have??? :confused:

I know the formula for calculating omega. but not for the others...:(

A thought:
Is it like, the variables we need to depend are based on the turbulence model we use?? like for k-eps we need k, eps...for k-w we need k and omega, for RSM mode we would need k, and R..and alphat in case of heat transfer is needed to be calculated????

You said you want to use the RNG k-epsilon model, as a consequence omega and R are not required. You need omega if you use for example k-omega, and you need R if you use Reynolds stress models.

You should have mut and alphat, which are specified using wall-functions at walls, and are calculated elsewhere. Take a look at this tutorial

/OpenFOAM-1.7.x/tutorials/compressible/rhoPimpleFoam/angledDuct/0/

It is for rhoPimpleFoam, however the only difference is in the iterative procedure. There you can see how alphat and mut are set up.

Correct. However the energy equation is always solved in compressible codes, as consequence you need alphat :-)

Best,

Thanks..that makes the things more clear now.

Actually i have worked on a k-eps research code for long time, and there such parameters as alpha or mut would be calculated internally, so i never had to specify them externally. So it didnt occur to me quickly..:(

But thank you very much for the quick replies.:D

I have one more doubt if i may trouble some more..

I need to have info on mixing, which means air/fuel mixture fraction distribution which obviously not gonna be there in rhoPisoFoam. So how to get this data?? One approach i could think of is to define a scaler which will have value 1 for fuel and 0 for air and thus can give the mixture fraction straightaway. But i think i will need to solve transport equation for that. Is there any better way?? If not how to do this scalar thing???

They are still calculated internally, but you can choose how in OpenFOAM ;-)


Depending on you needs you can add the scalar transport equation, or use one of the solvers which include species transport (reacting/diesel solvers), which are more flexible and can use chemical databases (chemkin).

Best,

Thanks, that would be all for now. i am setting up case for rhopisoFoam now.:o

Hello Alberto, :)

I am facing this error for rhoPisoFoam. Do you know what it means???


Thank you.

Let me grab my crystal ball :D

It means that there is some problem in your setup. Compare with a tutorial case, and start changing it step by step to adapt it to your case.

Best,

Hi,

I somehow managed to make the solver run. But now i am facing one more problem.

The Co number keeps increasing to large numbers. I cannot make the mesh more coarser and the Co problem makes the time step so small that the simulation will take thousand years to complete..:eek: (Thank god i have i7 processor;))

Now i am using the default solver settings of rhoPisoFoam from the tut of OF1.7. I want to know how can i tackle this Co problem without changing mesh and timestep, but modifying say relaxation factors and scheme only??? For this which scheme can be better or which UR factors can be good. (This is part of experience which i dont have much. :( so asking here).

Any suggestions???

PFA my case files.

Hi,

set "adjustTimeStep" to yes and perhaps start with smaller delta T

Regards,
Christian

Hi chris,

I have tried that. It does not help.

hii
 

i m doing work in supersonic mixing ... i need to solve density based solver... i was using fluent ..now i would like to do the work in open foam ......
inlet pressures
primary 10bar and
secondary 2bar
i need to solve a turbulent case
which foam should i use
looking for ur advices
if possible send me ans to my mailid
krishnaraja88@gmail.com