CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Abrupt jump at inlet in rhoCentralFoam for fixed pressure & velocity

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

Reply
 
LinkBack Thread Tools Display Modes
Old   November 9, 2016, 10:23
Default Abrupt jump at inlet in rhoCentralFoam for fixed pressure & velocity
  #1
Member
 
Join Date: Sep 2015
Location: Singapore
Posts: 42
Rep Power: 3
usv001 is on a distinguished road
I am trying to simulate compressible pipe flow in rhoCentralFoam. I am providing a fixed pressure and inlet velocity at the inlet but this causes an abrupt jump in properties from the inlet patch to the first cell and another jump from the first cell to the second.

I do not encounter the problem if I specify the outlet pressure (waveTransmissive with fieldInf), let the inlet pressure be waveTransmissive (without any fieldInf) and fix the inlet velocity. However, this is impractical since I would like to see how the flow develops from fixed inlet conditions.

I have a feeling that the problem is due to the manner in which fluxes are calculated in rhoCentralFoam. I could be wrong though.

Does anyone know of any correction to apply for the incoming flux to resolve this issue?

Many thanks.
usv001 is offline   Reply With Quote

Old   November 9, 2016, 13:51
Default
  #2
Senior Member
 
Mahdi Hosseinali
Join Date: Apr 2009
Location: NB, Canada
Posts: 270
Rep Power: 10
anishtain4 is on a distinguished road
Mathematically you can't do that. Read the book by Malalasgar about the boundary conditions.
anishtain4 is offline   Reply With Quote

Old   November 10, 2016, 01:29
Default
  #3
Member
 
Join Date: Sep 2015
Location: Singapore
Posts: 42
Rep Power: 3
usv001 is on a distinguished road
Thank you for the quick reply.

So, am I to understand that, currently, one cannot specify fixedValue inlet conditions for pressure and velocity simultaneously in rhoCentralFoam? And, is my guess, that the central-upwind flux scheme is the reason behind it, correct?

Many thanks!
usv001 is offline   Reply With Quote

Old   November 14, 2016, 14:20
Default
  #4
Senior Member
 
Mahdi Hosseinali
Join Date: Apr 2009
Location: NB, Canada
Posts: 270
Rep Power: 10
anishtain4 is on a distinguished road
You may find a better explanation of this in Hirsch, but this should also give you the idea:
Look up an-introduction-to-computational-fluid-dynamics by versteeg,
Read page 27, 2.6 Classification of physical behaviour.
anishtain4 is offline   Reply With Quote

Old   November 19, 2016, 00:08
Default
  #5
Member
 
Join Date: Sep 2015
Location: Singapore
Posts: 42
Rep Power: 3
usv001 is on a distinguished road
Thanks for the tip!

I have been reading up on hyperbolic conservation laws since I started to use rhoCentralFoam about a year ago. I still haven't quite gotten one aspect of it: the energy equation. I have already posted this question [Energy equation rhoCentralFoam (revisited)] before but none answered. I am repeating it here in the hope that you (or someone else) might direct me in the right direction.

In the inviscid energy equation, why does sigmaDotU, clearly a viscous term, appear? Furthermore, U here is interpolated using Kurganov's flux scheme which was developed for the hyperbolic conservation laws, Euler equations. Even the wave speeds are computed based on Euler equations. My understanding is that the flux scheme is to compute a numerical flux but here they have used it to interpolate just U. I don't quite follow...

Thanks once again.

USV
usv001 is offline   Reply With Quote

Reply

Tags
flux calculation, pressure and velocity, rhocentralfoam

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Wind tunnel Boundary Conditions in Fluent metmet FLUENT 3 October 8, 2016 08:41
Pressure loss Velocity coupling CFXMUFFIN CFX 1 February 6, 2016 05:43
Unphysical jump in pressure at the outlet predicted by rhoCentralFoam sathya123 OpenFOAM Verification & Validation 0 January 23, 2016 09:49
Pressure Inlet VS velocity Inlet difference Mohsin FLUENT 6 January 15, 2015 18:15
Pressure farfield vs Velocity Inlet BC DarrenC FLUENT 1 September 16, 2009 02:23


All times are GMT -4. The time now is 13:33.