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

Development of an All mach number fractional step method

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   December 16, 2016, 21:48
Default Development of an All mach number fractional step method
  #1
Member
 
Join Date: Jul 2012
Posts: 67
Rep Power: 13
Carlen is on a distinguished road
Dear Foamers,

I have recently come across this paper "A Primitive Variable Central Flux Scheme for All
Mach Number Flows" that triggered me to realise this solver in OpenFOAM.
Here is reference for this paper:

Shi, Ke; Morris, Scott; Jemcov, Aleksandar
A Primitive Variable Central Flux Scheme for All Mach Number Flows (Conference)
53rd AIAA Aerospace Sciences Meeting, 2015, At Kissimmee, FL, 2015.

However, I am stuck at the central flux part:

Code:
surfaceScalarField rho_left
(
     "rho_left",
      fvc::interpolate(rho, pos, "reconstruct(rho)")
);
surfaceScalarField rho_right
(
     "rho_right",
     fvc::interpolate(rho, neg, "reconstruct(rho)")
);

surfaceVectorField U_left
(
     "U_left",
     fvc::interpolate(U, pos, "reconstruct(U)")
);

surfaceVectorField U_right
(
     "U_right",
     fvc::interpolate(U, neg, "reconstruct(U)")
);

surfaceScalarField phil((rho_left*U_left) & mesh.Sf());
surfaceScalarField phir((rho_right*U_right) & mesh.Sf());
surfaceScalarField pRML(runTime.deltaT()*fvc::snGrad(p)*mesh.magSf());
surfaceScalarField phif(0.5*(phil + phir) + pRML);
which corresponds to Eqs 11-17 of the paper. With this implementation, the continuity error is at an order of 1e-02 which is not at all acceptable. I have attached the solver and a forwardStep case for you to help me with this problem.

Many thanks,
Carlen
Attached Files
File Type: gz cpbnsFSM.tar.gz (3.8 KB, 6 views)
File Type: gz forwardStep-cFSM.tar.gz (2.9 KB, 2 views)
Carlen is offline   Reply With Quote

Old   December 17, 2016, 04:49
Default
  #2
Senior Member
 
Santiago Lopez Castano
Join Date: Nov 2012
Posts: 354
Rep Power: 15
Santiago is on a distinguished road
Can you put the paper?

Sent from my GT-I8190L using CFD Online Forum mobile app
Santiago is offline   Reply With Quote

Old   December 17, 2016, 05:35
Default
  #3
Member
 
Join Date: Jul 2012
Posts: 67
Rep Power: 13
Carlen is on a distinguished road
Quote:
Originally Posted by Santiago View Post
Can you put the paper?

Sent from my GT-I8190L using CFD Online Forum mobile app
Here is the link for this paper: http://sci-hub.cc/10.2514/6.2015-1265
Carlen is offline   Reply With Quote

Old   December 17, 2016, 09:49
Default
  #4
Senior Member
 
Santiago Lopez Castano
Join Date: Nov 2012
Posts: 354
Rep Power: 15
Santiago is on a distinguished road
Thanks! I'll check it out this weekend and see how I can give you a hand

Sent from my GT-I8190L using CFD Online Forum mobile app
Santiago is offline   Reply With Quote

Old   December 18, 2016, 20:09
Smile
  #5
Member
 
Join Date: Jul 2012
Posts: 67
Rep Power: 13
Carlen is on a distinguished road
Anyone would like to give me a hand
Carlen 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
Foam::error::PrintStack almir OpenFOAM Running, Solving & CFD 91 December 21, 2022 05:50
[snappyHexMesh] Error snappyhexmesh - Multiple outside loops avinashjagdale OpenFOAM Meshing & Mesh Conversion 53 March 8, 2019 10:42
Local mach number jacobH Main CFD Forum 3 October 1, 2014 11:20
[blockMesh] --> foam fatal error: lillo763 OpenFOAM Meshing & Mesh Conversion 0 March 5, 2014 11:27
Micro Scale Pore, icoFoam gooya_kabir OpenFOAM Running, Solving & CFD 2 November 2, 2013 14:58


All times are GMT -4. The time now is 03:15.