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

wavetransmissive boundary for pressure and velocity

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

Like Tree1Likes
  • 1 Post By SanchoPanza

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   December 8, 2014, 09:28
Default wavetransmissive boundary for pressure and velocity
  #1
New Member
 
Ignacio Duran
Join Date: Dec 2014
Posts: 2
Rep Power: 0
SanchoPanza is on a distinguished road
Dear all,

This is my first post here since I'm new to OPENFOAM (though not new to CFD).

My aim is to perform a rather simple compressible simulation using rhoCentralFoam. The geometry is a simple square with zeroGradient boundaries left and right. I've added a rho*g term in the momentum equation to account for gravity. I initialize at t=0 with a pressure gradient to account for this gravity term.

I would like to use wavetransmissive boundary conditions (at the top and at the bottom). Its use is described in a rather simple way in every post I've found:

0/p:
top
{
type waveTransmissive;
field p;
phi phi;
rho rho;
psi thermosi_0;
gamma 1.4;
fieldInf 99884.8;
lInf 100;
value uniform 99884.8;
}

where lInf represents how far the boundary condition is (I'm 100% aware of how the NSCBC system of wavetransmissive works!)

However, OpenFoam requires a boundary condition for U as well! This is surprising, since the NSCBC system is written in a way in which one will prescribe p only (it can be written to prescribe U also, but not both at the same time).

I've tried several boundary conditions with no success:
1. calculated: this will not run (though, by the name it was my first guess)
2. pressureInletVelocity (or pressureInletOutletVelocity), which runs, but pressure oscillates without leaving the domain (this indicates that the boundary condition is reflective)
3. wavetransmissive in U: where pressure/velocity drifts away slowly from the expected value.

My questions are then:

What is the 'partner' boundary condition in velocity when one specifies wavetransmissive in pressure?

can wavetransmissive be used with a gravity term? (I am aware that the NSCBC conditions have to compute gradients normal to the flow... this could be the reason for my problems...


Many thanks for your help!

Sancho
lukasf likes this.
SanchoPanza is offline   Reply With Quote

Old   December 9, 2014, 05:23
Default
  #2
Senior Member
 
Olivier
Join Date: Jun 2009
Location: France, grenoble
Posts: 272
Rep Power: 17
olivierG is on a distinguished road
Hello,

I would say you should use also wavetransmissive for U, but you can try "advective" BC for U instead.
NB: you may should use advective for other field too.

regards,
olivier
olivierG is offline   Reply With Quote

Old   December 9, 2014, 07:37
Default
  #3
New Member
 
Ignacio Duran
Join Date: Dec 2014
Posts: 2
Rep Power: 0
SanchoPanza is on a distinguished road
Thanks for your help!

It seems indeed that the best choice would be to use the wavetransmissive BC in U as well... though this is not obvious, since a single wavetransmissive BC should set a combination of P and U, and this is enough. Having 2 of them is therefore redundant.

I've had a look at the rhoCentralFoam solver, and it seems that the boundary conditions on U and on P are computed separately. This may be the reason why I have to set wavetrasmissive in both fields.

Many thanks!
SanchoPanza is offline   Reply With Quote

Old   September 7, 2020, 11:46
Default what is the correct approch for using wave transmissive boundary condition in OPENFOA
  #4
New Member
 
zein elserfy
Join Date: May 2018
Posts: 25
Rep Power: 7
zeinelserfy is on a distinguished road
do you test your case? is the best way is to apply wave transmissive for both p and u?

Quote:
Originally Posted by SanchoPanza View Post
Thanks for your help!

It seems indeed that the best choice would be to use the wavetransmissive BC in U as well... though this is not obvious, since a single wavetransmissive BC should set a combination of P and U, and this is enough. Having 2 of them is therefore redundant.

I've had a look at the rhoCentralFoam solver, and it seems that the boundary conditions on U and on P are computed separately. This may be the reason why I have to set wavetrasmissive in both fields.

Many thanks!
zeinelserfy is offline   Reply With Quote

Old   April 11, 2021, 06:11
Default
  #5
New Member
 
Join Date: May 2020
Posts: 4
Rep Power: 5
labyrinth01 is on a distinguished road
Quote:
Originally Posted by SanchoPanza View Post
Dear all,

This is my first post here since I'm new to OPENFOAM (though not new to CFD).

My aim is to perform a rather simple compressible simulation using rhoCentralFoam. The geometry is a simple square with zeroGradient boundaries left and right. I've added a rho*g term in the momentum equation to account for gravity. I initialize at t=0 with a pressure gradient to account for this gravity term.

I would like to use wavetransmissive boundary conditions (at the top and at the bottom). Its use is described in a rather simple way in every post I've found:

0/p:
top
{
type waveTransmissive;
field p;
phi phi;
rho rho;
psi thermosi_0;
gamma 1.4;
fieldInf 99884.8;
lInf 100;
value uniform 99884.8;
}

where lInf represents how far the boundary condition is (I'm 100% aware of how the NSCBC system of wavetransmissive works!)

However, OpenFoam requires a boundary condition for U as well! This is surprising, since the NSCBC system is written in a way in which one will prescribe p only (it can be written to prescribe U also, but not both at the same time).

I've tried several boundary conditions with no success:
1. calculated: this will not run (though, by the name it was my first guess)
2. pressureInletVelocity (or pressureInletOutletVelocity), which runs, but pressure oscillates without leaving the domain (this indicates that the boundary condition is reflective)
3. wavetransmissive in U: where pressure/velocity drifts away slowly from the expected value.

My questions are then:

What is the 'partner' boundary condition in velocity when one specifies wavetransmissive in pressure?

can wavetransmissive be used with a gravity term? (I am aware that the NSCBC conditions have to compute gradients normal to the flow... this could be the reason for my problems...


Many thanks for your help!

Sancho
This message is made for the followers since it is a long time after the last reply.

There is one example case in OF (reactingFoam-sandiaD). It uses 'waveTransmissive' for 'p' but 'pressureInletOutletVelocity' for 'u'.
labyrinth01 is offline   Reply With Quote

Old   June 29, 2021, 07:48
Default So it works?
  #6
New Member
 
deewakar
Join Date: Dec 2011
Posts: 18
Rep Power: 14
alligngr8 is on a distinguished road
Quote:
Originally Posted by labyrinth01 View Post
This message is made for the followers since it is a long time after the last reply.

There is one example case in OF (reactingFoam-sandiaD). It uses 'waveTransmissive' for 'p' but 'pressureInletOutletVelocity' for 'u'.
So it implies that waveTransmission for pressure and zero gradient for velocity.
However, I have come across few articles which say this may not be very well implemented in openFoam and the best is to implement NSCBC.

Any leads or comments on that please?
alligngr8 is offline   Reply With Quote

Reply

Tags
boundary conditions, openfoam 2.1.x, wavetransmissive

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
Access to the velocity at boundary cells mathslw OpenFOAM 14 January 10, 2018 00:31
Domain Imbalance HMR CFX 5 October 10, 2016 06:57
GETVAR Error in Multiband Monte Carlo Radiation Simulation with Directional Source silvan CFX 3 June 16, 2014 10:49
help with velocity boundary condition potentialFoam hfc OpenFOAM Running, Solving & CFD 2 May 31, 2012 20:52
RPM in Wind Turbine Pankaj CFX 9 November 23, 2009 05:05


All times are GMT -4. The time now is 05:46.