CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Running, Solving & CFD

zeroGradient pressure and moving mesh

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

LinkBack Thread Tools Display Modes
Old   July 22, 2009, 09:28
Default zeroGradient pressure and moving mesh
New Member
Christophe Kassiotis
Join Date: Mar 2009
Location: Paris
Posts: 17
Rep Power: 8
kassiotis is on a distinguished road
Dear all,

I'm currently trying to solved one of the benchmark fsi problem proposed by Wall and co-authors (see for instance Kütller and Wall (2008) Fixed-point fluid–structure interaction solvers with dynamic relaxation doi:10.1007/s00466-008-0255-5). The problem consist in a modified lid-driven cavity with a moving bottom, with two unconstrained nodes at the top right and left vertical walls. In the first attached picture, the discretized problem by 32x32 elements is represented, as well as the boundary condition (Picture from the publication given above). The time step considered is 0.1s

For the sake of simplicity, let us consider for the moment the fluid problem only. For the velocity (U) boundary condition, I consider imposed velocity (0 0 0) at walls, (1-cos(2×Pi×t/T) 0.0 0.0) at the top and zeroGradient for the unconstrained faces.
As one can see, the pressure (p) boundary conditions are not define, and therefore, I used zeroGradient everywhere. It is well known that then p is defined to a constant, and playing with pRefCell and pRefValue, i managed to get approximatively good pressure forces ( p×dS ) at the bottom of my cavity (fixing the top right cell to 0).

The problem start when I want to move the mesh. I consider for the moment only imposed motion of the bottom (0.0 L(1-cos(2×Pi×t/T)) 0.0) with L = 0.1. Then, if i'm not fixing one of the pressure, the computation quickly crash (at t=0.9s)

Continuity error cannot be removed by adjusting the outflow.
Please check the velocity boundary conditions and/or run potentialFoam to initialise the outflow.
Total flux              : 0.0509994
Specified mass inflow   : 0.000944722
Specified mass outflow  : 1.87515e-18
Adjustable mass outflow : 0
The problem is that when I fixed the pressure to one (or two) of the output, then the pressure oscillations become to big compared to what they supposed to be.

My questions are the following:
Q1 : is it normal that the problem occurs in Finite Volume, and not with their stabilized Finite Element ? I get a message from Küttler that say that the pressure is not fixed
in their computation, and that pose no problem. Maybe I'm wrong, but this does not depend of the discretisation technic used, as the continuum equation are the same ? So maybe something miss here...
Q2 : as the force applied to the body directly depend of p (and not of its gradient) the constant used here to fixe p is of importance, isn't it ?
Q3 : is it possible to solve a problem with icoDyMFoam and zeroGradient pressure every where for the pressure ? Or is my crash a desired behavior ? Maybe is it hard to adjust phi when the domain is moving ? Is their a way to fix this problem ?
Q4 : as someone any advices for the boundary condition ? I have to say that I try many of them (fixedValue, totalPressure at the inlet, one gives results as close as them for the force at the bottom as the one with zeroGradient) ?

Thank you in advance.

Attached Images
File Type: png problemDefinition.png (29.4 KB, 120 views)
File Type: png p1023.png (6.6 KB, 122 views)
kassiotis 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
pressure oscillations in moving mesh cases david OpenFOAM Running, Solving & CFD 5 April 25, 2011 11:16
moving (dynamic) mesh gives fluctuating pressure daniel FLUENT 1 July 31, 2010 14:32
VOF with a moving mesh Jeremie FLUENT 1 November 26, 2008 09:55
CavitatingFoam with a moving mesh mass conservation error idosil OpenFOAM Running, Solving & CFD 3 November 27, 2007 18:27
Gas pressure question Dan Moskal Main CFD Forum 0 October 24, 2002 22:02

All times are GMT -4. The time now is 06:52.