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

All Mach number implicit solver with Kurganov-Tadmore scheme - pisoCentralFoam

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

Like Tree20Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 15, 2015, 11:38
Default
  #81
Senior Member
 
mkraposhin's Avatar
 
Matvey Kraposhin
Join Date: Mar 2009
Location: Moscow, Russian Federation
Posts: 343
Rep Power: 18
mkraposhin is on a distinguished road
Quote:
Originally Posted by chriss85 View Post
Just a short question on this topic. Why are you not separating the flux in positive/negative parts in the continuity equation?
Code:
solve
    (
        fvm::ddt(rho) + fvc::div(phi)
    );
Because phi = phiPos + phiNeg
This means, that for continuity equation, formulated with explicit scheme, both of two versions are equals mathematically and numerically:
div(phi) = div(phiPos + phiNeg) = div(phiPos) + div(phiNeg)

Last edited by mkraposhin; November 16, 2015 at 07:17. Reason: grammar
mkraposhin is offline   Reply With Quote

Old   November 16, 2015, 06:40
Default
  #82
Senior Member
 
Join Date: Oct 2013
Posts: 394
Rep Power: 15
chriss85 will become famous soon enough
Thanks, that makes sense. I need to read up on the Kurganov-Tadmor schemes when I find some time.
chriss85 is offline   Reply With Quote

Old   December 15, 2015, 11:55
Default
  #83
New Member
 
Lewis
Join Date: Dec 2015
Posts: 14
Rep Power: 7
levvis is on a distinguished road
Hi Matvey

I'm Lewis, a student from the UK studying Aerospace Engineering.

For my Final Year Project, I have been allocated "Development of an OpenFOAM code for Aerospace Engineering".

I have run some tutorial cases and got to grips with OpenFOAM and ParaView, using Ubuntu and the terminal window.

As your code is a newer hybrid solver, with scope for expansion, I am very interested in using your published code as a foundation. Therefore I will adapt it to my own need, the beauty of OpenFOAM.

I see you've already run cases using a cylinder, which is ideal for me, as my plan with your code is to run it with a cylinder within the geometry. And progress this shape into an aerofoil, for greater use in the Aviation Industry.


What are your opinions on use of your code to adapt and create this application? As I would like to be able to run the simulations at a sub-sonic velocity, and then supersonic.

Many thanks in advance
levvis is offline   Reply With Quote

Old   December 16, 2015, 04:47
Default
  #84
Senior Member
 
mkraposhin's Avatar
 
Matvey Kraposhin
Join Date: Mar 2009
Location: Moscow, Russian Federation
Posts: 343
Rep Power: 18
mkraposhin is on a distinguished road
Hi, Lewis

I'm not sure that i undestand your question correctly.

Because you are talking about
Code:
Development of an OpenFOAM code for Aerospace Engineering
,
or
Code:
 I am very interested in using your published code as a foundation
but you didn't mentioned which physical case you want to consider. So i can't understand for which case you want to extend pisoCentralFoam solver.

As for simulation of external flow over bodies, i can tell you that this solver was test for simple and complex geomtries for Mach humbers from 0.01 to 6 by me and my colleagues.

So, if you need to use it only for gas dynamics, then you don't need to change it. In this case, we need to create a separate thread on forum in "Solving" section for questions about setup of numerical case.

And if you are talking about extending this solver, we can proceed here, but you must supply more information on your physical problem
mkraposhin is offline   Reply With Quote

Old   December 16, 2015, 11:09
Default
  #85
New Member
 
Lewis
Join Date: Dec 2015
Posts: 14
Rep Power: 7
levvis is on a distinguished road
Hi Matvey

My question is, do you think there is an opportunity for me to use this for airflow over an aerofoil, at transonic speeds? Will it converge?

I suppose you have already answered that question, when you said you and your colleagues used it with simple and complex geometries.

I'd like to use this solver, just have to edit some of the code to show I know what I'm doing, and why I did it.


Thanks

Lewis
levvis is offline   Reply With Quote

Old   December 16, 2015, 14:51
Default
  #86
Senior Member
 
mkraposhin's Avatar
 
Matvey Kraposhin
Join Date: Mar 2009
Location: Moscow, Russian Federation
Posts: 343
Rep Power: 18
mkraposhin is on a distinguished road
Quote:
Originally Posted by levvis View Post
Hi Matvey

My question is, do you think there is an opportunity for me to use this for airflow over an aerofoil, at transonic speeds? Will it converge?

I suppose you have already answered that question, when you said you and your colleagues used it with simple and complex geometries.

I'd like to use this solver, just have to edit some of the code to show I know what I'm doing, and why I did it.


Thanks

Lewis
'

Hi, i think, yes, it will converge. We have plans to make it LTS, but not for this year.
Also, we make publication YSC 2015 procedia, so you can download pdf from here

PisoCentralFoam at researchgate

or from elsevier

doi:10.1016/j.procs.2015.11.007

Also, we made application of this scheme for multicomponent reacting flow

It will be nice if you will help us during your work to investigate range of application of this scheme (solver)
mkraposhin is offline   Reply With Quote

Old   December 18, 2015, 06:15
Default
  #87
Senior Member
 
mkraposhin's Avatar
 
Matvey Kraposhin
Join Date: Mar 2009
Location: Moscow, Russian Federation
Posts: 343
Rep Power: 18
mkraposhin is on a distinguished road
Dear colleagues!

We made update of pisoCentralFoam solver for OpenFOAM 3.0.0. In new version LTS time scheme is included and pisoCentralFoam can be used for steady state simulations
hk318i, Alexee and OlegSutyrin like this.
mkraposhin is offline   Reply With Quote

Old   January 14, 2016, 04:17
Default
  #88
New Member
 
sethu
Join Date: Oct 2011
Posts: 7
Rep Power: 11
Sethu is on a distinguished road
Dear Matvey Kraposhin,

Thank you very much for your pisoCentralFoam LTS version. Im wondering whether I can use your solver for application like steady state simulation of NASA rotor 37 is an isolated transonic axial compressor rotor where i need to account for rotational effect either in SRF or MRF (for rotor+stator). Please give me your comments on this.

I had modified rhoSimpleFoam as rhoSRFSimpleFoam for a single rotor case. But, when i executed the solver im getting error as follows,
[Maximum number of iterations exceeded

From function thermo<Thermo, Type>::T(scalar f, scalar T0, scalar (thermo<Thermo, Type>::*F)(const scalar) const, scalar (thermo<Thermo, Type>::*dFdT)(const scalar) const, scalar (thermo<Thermo, Type>::*limit)(const scalar) const) const
in file /home/nrcs-cae/OpenFOAM/OpenFOAM-2.3.x/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 76.

FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::error::abort() at ??:?
#2 Foam::species::thermo<Foam::hConstThermo<Foam:er fectGas<Foam::specie> >, Foam::sensibleInternalEnergy>::TEs(double, double, double) const at ??:?
#3 Foam::hePsiThermo<Foam:siThermo, Foam:ureMixture<Foam::constTransport<Foam::speci es::thermo<Foam::hConstThermo<Foam:erfectGas<Foa m::specie> >, Foam::sensibleInternalEnergy> > > >::calculate() at ??:?
#4 Foam::hePsiThermo<Foam:siThermo, Foam:ureMixture<Foam::constTransport<Foam::speci es::thermo<Foam::hConstThermo<Foam:erfectGas<Foa m::specie> >, Foam::sensibleInternalEnergy> > > >::correct() at ??:?
#5
at ??:?
#6 __libc_start_main in "/lib/i386-linux-gnu/libc.so.6"
#7
at ??:?
Aborted (core dumped)


Any advice on selecting solvers for my steady state (SRF,MRF) turbo machinery compressible cases?


Regards,

Sethu
Sethu is offline   Reply With Quote

Old   January 14, 2016, 10:44
Default
  #89
Senior Member
 
mkraposhin's Avatar
 
Matvey Kraposhin
Join Date: Mar 2009
Location: Moscow, Russian Federation
Posts: 343
Rep Power: 18
mkraposhin is on a distinguished road
Hi, Sethu

I added MRF capabilities to pisoCentralFoam, you can download new solver and tutorials from github.

As for rhoSRFSimpleFoam, you must create new thread, because this thread dedicated to pisoCentralFoam
mkraposhin is offline   Reply With Quote

Old   March 29, 2016, 19:30
Default Species diffusion in reactingCentralFoam
  #90
New Member
 
Bryan Schmidt
Join Date: Sep 2013
Posts: 18
Rep Power: 10
Bryan S is on a distinguished road
Dear Dr. Kraposhin,
I have been using reactingCentralFoam to compute a nonreacting supersonic flow with foreign gas injection in a boundary layer. The mixtures are binary, but the foreign gas often has a molecular weight that is quite different from that of the free stream gas (nitrogen). I have a question about how species diffusion is handled. It seems to me that a unity Schmidt number assumption is made using the mixture-averaged viscosity, such that all species effectively have the same diffusivity. Am I correct in this understanding?

I also see that the equations invoke Fick's Law for diffusion, except I don't see a correction velocity as is typically done when mass fractions are implemented as opposed to mole fractions in Fick's Law. If I'm correct, do you have a sense of what the errors are associated with this both in the species transport and energy equations? Thanks.
Bryan S is offline   Reply With Quote

Old   April 6, 2016, 12:16
Default
  #91
New Member
 
Lewis
Join Date: Dec 2015
Posts: 14
Rep Power: 7
levvis is on a distinguished road
Hi Matvey

I'm struggling to understand this error, all I have done is simply created a new geometry (a double-wedge aerofoil) in GMSH and meshed it. And then substituted this in for the cylinder geometry.

I've changed a few properties too, like rho = 1.225 and velocity to simulate conditions at Sea Level. And upon running the command I get a printstack error:
Quote:
Mean and max Acoustic Courant Numbers = 0.0111359 0.102596
Courant Number mean: 0.0046786 max: 35.2339
Time = 0.00248

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::GaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) at ??:?
#4 Foam::GaussSeidelSmoother::smooth(Foam::Field<doub le>&, Foam::Field<double> const&, unsigned char, int) const at ??:?
#5 Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#6 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:?
#7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
#8 ? at ??:?
#9 ? at ??:?
#10 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#11 ? at ??:?
Floating point exception (core dumped)
i haven't changed anything else, so haven't even delved into the physics of the solver.

Can you, or anyone else, please advise?

Many thanks
levvis is offline   Reply With Quote

Old   April 7, 2016, 04:41
Default
  #92
Senior Member
 
Join Date: Oct 2013
Posts: 394
Rep Power: 15
chriss85 will become famous soon enough
I suggest trying another linear equation solver. Matvey (or someone else from his group) has ported the BiCGStab solver which is more stable and crashes less often. Check out the gits, it's in one of them.
chriss85 is offline   Reply With Quote

Old   April 10, 2016, 06:14
Default
  #93
New Member
 
Lewis
Join Date: Dec 2015
Posts: 14
Rep Power: 7
levvis is on a distinguished road
Quote:
Originally Posted by chriss85 View Post
I suggest trying another linear equation solver. Matvey (or someone else from his group) has ported the BiCGStab solver which is more stable and crashes less often. Check out the gits, it's in one of them.
I shall try this too, but just trying to problem solve the pisoCentralFoam.
It has been edited ever so slightly and is now running without crashing, but the flow is wrong.
It's an MH45 aerofoil at an angle of attack of 10 degrees.

But when running it, the lower velocity just seems to grow out from the aerofoil.
Has anyone seen this before? Can anyone advise? Am I missing files or something?

Picture attached

EDIT: Problem solving files, the time step is very small and thus the end time is small, this is probably the reason for the simulation? Because the simulation is correct in respective to such a small time step?
Am I right in thinking so? I've changed it and have it running but obviously it's a long compute time, even with a higher time step.

Just thinking as the difference between the scale is very small, and looking at solvers like simpleFoam aerofoil case, if converted to a small time step, it shows a similar picture.
Attached Images
File Type: jpg Screenshot from 2016-04-09 21:18:22.jpg (86.5 KB, 50 views)

Last edited by levvis; April 10, 2016 at 10:48.
levvis is offline   Reply With Quote

Old   May 30, 2016, 14:57
Default
  #94
New Member
 
Jim
Join Date: Feb 2014
Location: UK
Posts: 22
Rep Power: 9
jimteb is on a distinguished road
Hi Matvey,

I saw from an earlier post that pisoCentralFoam is not bound to the acoustic Courant number, but rather the flow Courant number.

Does this mean that this solver is unsuitable for cases where you are interested in acoustic waves?

Thanks,

James
jimteb is offline   Reply With Quote

Old   June 3, 2016, 11:03
Default
  #95
Senior Member
 
mkraposhin's Avatar
 
Matvey Kraposhin
Join Date: Mar 2009
Location: Moscow, Russian Federation
Posts: 343
Rep Power: 18
mkraposhin is on a distinguished road
Hi,

We tested this solver for several simple cases (monopole and diploe). It was found, that it can be used to reproduce acoustic waves if you will preserve next conditions:

1) Acoustics courant < 1 (Or better < 0.1)
2) Mesh resolution is > 20 cells per resolved wave
fumiya likes this.
mkraposhin is offline   Reply With Quote

Old   June 21, 2016, 05:34
Default
  #96
New Member
 
Myth
Join Date: Jul 2011
Location: Unitek kingdom
Posts: 17
Rep Power: 12
mg.mithun is on a distinguished road
Send a message via Yahoo to mg.mithun Send a message via Skype™ to mg.mithun
Hi mkraposhin,
Have you tried this solver for cavitating flows?
hk318i likes this.
mg.mithun is offline   Reply With Quote

Old   July 8, 2016, 16:48
Default
  #97
New Member
 
Eduardo Alvarez
Join Date: Dec 2015
Posts: 3
Rep Power: 7
AlvarezEdo is on a distinguished road
Hi Mkraposhin,

I'm modeling a shock tube with multiple sections of different gases. I've been an OpenFOAM front-end user for a while, but I have never done anything complicated like mixing two solvers. I'm using your solver to model the shock with turbulent boundary layers along the shock tube, and is working very well, but now I need to add the different immiscible stages of gases.

How should I go about doing so? Is there a specific solver (compressibleMultiphaseInterFoam, for example) you would recommend mixing into your solver?



Thanks for your help, and great job in developing pisoCentralFoam!
AlvarezEdo 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
DPMFoam - Serious Error --particle-laden flow in simple geometric config benz25 OpenFOAM Running, Solving & CFD 27 December 19, 2017 20:47
foam-extend_3.1 decompose and pyfoam warning shipman OpenFOAM 3 July 24, 2014 08:14
Solver is finishing with huge Mach number Fonzie CFX 1 March 12, 2007 14:15
High Mach number solver error Luke CFX 3 January 31, 2007 22:26
TVD scheme at low Mach number Axel Rohde Main CFD Forum 5 August 6, 1999 02:01


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