CFD Online Logo CFD Online URL
Home > Forums > CFX

x,y,z,t- Dependent Pressure Initialization Using CEL

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

LinkBack Thread Tools Display Modes
Old   March 26, 2015, 14:51
Default x,y,z,t- Dependent Pressure Initialization Using CEL
Kegan Leckness
Join Date: Mar 2015
Posts: 38
Rep Power: 4
KeganLeckness is on a distinguished road
Hi all,

I'm a graduate student simulating pressure propagation through the ear. I've used an Opening boundary condition at the inlet of the ear canal and
if(t>0.0 [s] && t<0.003 [s], (-6666.6666666666666666666666667 [kPa/s] * t) +20 [kPa], 0.0 [Pa])
to initialize the pressure.

The problem I'm having is due to an artificial wall at the Opening. Although the manuals state that the Opening BC does not add artificial walls like the Outlet or Inlet BCs do, it appears that the Opening condition does add a wall when a pressure like the one above is set at the location of the Opening. I get the negative reflected pressure, as discussed in the description of the artificial walls inherent in Outlet and Inlet conditions. I want to get away from this.

I've tried to implement the beta-feature Acoustic Reflectivity, but this is not applicable. I've also tried to separate a domain for pressure initialization, but this only applies the pressure for t=0 and not t=0 to t=3 ms as desired.

I wish to add spacial constraints to the time constraint listed in the above IF statement so that I may apply this pressure at various locations in the canal's domain. I've read the various guides discussing CCL, CEL, Expressions, Functions, and Routines, but I am having a difficult time understanding and putting it all together.

Would someone be kind enough to elucidate the relationship between the aforementioned topics and tell me what I need to define a pressure depending on x,y,z,&t? If I need to use a library, please explain the purpose of that, too.

I greatly appreciate any assistance,
KeganLeckness is offline   Reply With Quote

Old   March 26, 2015, 17:56
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 13,535
Rep Power: 104
ghorrocks is a jewel in the roughghorrocks is a jewel in the roughghorrocks is a jewel in the rough
There are several options for pressure openings. I was not aware that any of them generated walls to stop reverse flow. So if you are having problems with the boundary condition you are using I would consider changing to a different type. For instance there is opening pressure and direction versus Static pressure and direction, and there is entrainment. Try a few combinations of these.

The problem with all these approaches is that they are going to reflect pressure waves and generate spurious waves. I would investigate the non-reflecting boundary condition beta feature. I would also investigate moving the opening boundary away from the area of interest so the reflected waves do not reach the region of interest (or at least are minimised).
ghorrocks is offline   Reply With Quote


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
question regarding LES of pipe flow - pimpleFoam Dan1788 OpenFOAM Running, Solving & CFD 35 December 14, 2016 05:34
Calculation of the Governing Equations Mihail CFX 7 September 7, 2014 06:27
CEL, CP, Pressure Saleeman1234 CFX 10 August 11, 2014 10:12
time dependent pressure BC on outlet sometimes1900 CFX 4 October 24, 2008 11:14
Hydrostatic pressure in 2-phase flow modeling (long) DS & HB Main CFD Forum 0 January 8, 2000 16:00

All times are GMT -4. The time now is 12:17.