CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   Main CFD Forum (https://www.cfd-online.com/Forums/main/)
-   -   ALE problem (https://www.cfd-online.com/Forums/main/6328-ale-problem.html)

Neil July 7, 2003 05:17

ALE problem
 
I have implemented an ALE formulation into a 2D Euler solver, which Roe's Flux difference splitting and a ghost cell method for enforcing boundary conditions. I have found that when solid boundaries at moved at certain rates the solution fails (due to negative pressure developing in some cells). I would be grateful if anyone could tell me why this might be occurring. Thanks for your help

Neil

Kourosh Hejazi July 7, 2003 06:47

Re: ALE problem
 
Do you update the geometry and re-distribute the cells in every time-step? That is needed for interior cells for each time-step. The other thing that I can think of is when the boudnary is moving very rapidly and the deformation happens in more than one layer of the grids, in which case (for free-surafec for example) is best to keep the changes within one layer first and get the simulation running. This avoids the mesh crossing in sharp gradinets. If you think this is any close to what you think may cause the problem then we can discuss it in more detail.

Neil July 7, 2003 07:14

Re: ALE problem
 
The boundary geometry is updated after each time-step and laplacian smoothing is preformed to redistribute the interior points. At the point of failure of the solution the grid quality is good, so this does not seem to be the problem. Also the deformationis such that the boundary does not cross the first layer

Tamer July 7, 2003 13:17

Re: ALE problem
 
Explain your problem more. Are you solving for a phase change such as solidification, why would the solid boundary move in reality? The solution you gets may be unrealistic because the boundary conditions are unrealistic.

Neil July 8, 2003 03:59

Re: ALE problem
 
I am attempting to model the simple case of aerofoil undergoing forced oscillations of a sinusoidal nature. I found that as the amplitude of the oscilations was increased, but the period kept constant, there was a point at which the solution failed due to negative pressure, as explained earlier

Neil


All times are GMT -4. The time now is 15:32.