Abrupt jump at inlet in rhoCentralFoam for fixed pressure & velocity

usv001 · November 9, 2016, 09:23

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.

anishtain4 · November 9, 2016, 12:51

Mathematically you can't do that. Read the book by Malalasgar about the boundary conditions.

usv001 · November 10, 2016, 00:29

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!

anishtain4 · November 14, 2016, 13:20

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.

usv001 · November 18, 2016, 23:08

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

November 9, 2016, 09:23	Abrupt jump at inlet in rhoCentralFoam for fixed pressure & velocity	#1
usv001 Senior Member Join Date: Sep 2015 Location: Singapore Posts: 102 Rep Power: 10	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.

November 14, 2016, 13:20		#4
anishtain4 Senior Member Mahdi Hosseinali Join Date: Apr 2009 Location: NB, Canada Posts: 273 Rep Power: 18	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. febriyan91 likes this.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Pressure Inlet VS velocity Inlet difference	Mohsin	FLUENT	9	January 4, 2021 10:34
Wind tunnel Boundary Conditions in Fluent	metmet	FLUENT	6	October 30, 2019 12:23
Pressure loss Velocity coupling	CFXMUFFIN	CFX	1	February 6, 2016 04:43
Unphysical jump in pressure at the outlet predicted by rhoCentralFoam	sathya123	OpenFOAM Verification & Validation	0	January 23, 2016 08:49
Pressure farfield vs Velocity Inlet BC	DarrenC	FLUENT	1	September 16, 2009 02:23

November 9, 2016, 12:51		#2
anishtain4 Senior Member Mahdi Hosseinali Join Date: Apr 2009 Location: NB, Canada Posts: 273 Rep Power: 18	Mathematically you can't do that. Read the book by Malalasgar about the boundary conditions.

November 10, 2016, 00:29		#3
usv001 Senior Member Join Date: Sep 2015 Location: Singapore Posts: 102 Rep Power: 10	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!

November 18, 2016, 23:08		#5
usv001 Senior Member Join Date: Sep 2015 Location: Singapore Posts: 102 Rep Power: 10	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