
[Sponsors] 
July 22, 2009, 09:28 
zeroGradient pressure and moving mesh

#1 
New Member
Christophe Kassiotis
Join Date: Mar 2009
Location: Paris
Posts: 17
Rep Power: 9 
Dear all,
I'm currently trying to solved one of the benchmark fsi problem proposed by Wall and coauthors (see for instance Kütller and Wall (2008) Fixedpoint fluid–structure interaction solvers with dynamic relaxation doi:10.1007/s0046600802555). The problem consist in a modified liddriven 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, (1cos(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(1cos(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) Code:
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.87515e18 Adjustable mass outflow : 0 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, inletOutlet...no one gives results as close as them for the force at the bottom as the one with zeroGradient) ? Thank you in advance.  ck 

Thread Tools  
Display Modes  


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 