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

Total Pressure BC for High-speed compressible flow - erroneous activation conditions

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 5, 2019, 17:14
Default Total Pressure BC for High-speed compressible flow - erroneous activation conditions
  #1
New Member
 
A. Z.
Join Date: Jan 2019
Posts: 1
Rep Power: 0
Zi_Aun is on a distinguished road
Hello,
Trying to derive the totalPressure boundary condition for a specific use, I found that the conditions for the selection between transonic and supersonic regimes for High-speed compressible flow are not consistent with the aerodynamic definition of these regimes. These conditions are associated with \gamma (specific heat ratios), where they should depend on the speed at least (Mach number).
See the script here below from: fields/fvPatchFields/derived/totalPressure/totalPressureFvPatchScalarField.C (OpenFOAM-5.x)

################################################## ##########################
Code:
        else
          {
              // High-speed compressible flow
   
              const fvPatchField<scalar>& psip =
                  patch().lookupPatchField<volScalarField, scalar>(psiName_);
   
              if (gamma_ > 1)
              {
                  scalar gM1ByG = (gamma_ - 1)/gamma_;
   
                  operator==
                  (
                      p0p
                     /pow
                      (
                          (1.0 + 0.5*psip*gM1ByG*(1.0 - pos0(phip))*magSqr(Up)),
                          1.0/gM1ByG
                      )
                  );
              }
              else
              {
                  operator==(p0p/(1.0 + 0.5*psip*(1.0 - pos0(phip))*magSqr(Up)));
              }
          }
################################################## #########################
This may induce errors, any comment or suggestion regarding the way to proceed with this?
Thanks.
Zi Aun.

Last edited by Zi_Aun; August 5, 2019 at 20:16. Reason: Add of the traditional "Hello"
Zi_Aun is offline   Reply With Quote

Reply

Tags
gamma, totalpressure bc

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
Pressure Outlet Targeted Mass Flow Rate LuckyTran FLUENT 1 November 23, 2016 11:40
Wrong flow in ratating domain problem Sanyo CFX 17 August 15, 2015 07:20
Problem with SIMPLEC-like finite volume channel flow boundary conditions ghobold Main CFD Forum 3 June 15, 2015 12:14
Boundary Conditions : Total Pressure or Velocity Gearb0x OpenFOAM Running, Solving & CFD 2 February 28, 2011 22:18
Total pressure in real gas (compressible flow) Bart Prast Main CFD Forum 3 November 14, 2000 11:44


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