CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Error caused by unexpected mesh displacement (OpenFOAM)

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 19, 2017, 13:02
Default Error caused by unexpected mesh displacement (OpenFOAM)
  #1
New Member
 
Join Date: Dec 2010
Posts: 4
Rep Power: 8
pille is on a distinguished road
Hello everyone,

I've got a problem probably caused by the mesh displacement in OpenFOAM 4.1.
It only occurs when the displacement is relatively high and I really don't understand why it moves the points in this way.

The case:
- 2D
- rectangular box with walls on all sides
- a piston in the middle which moves horizontally with a given frequency and amplitude
- the mesh was made by blockMesh
- the solver is "compressibleInterDyMFoam" because there are two phases (water, air)
and compressibility may have an influence at high frequencies
- a simplyfied mesh is shown in the attached picture


Description of the problem:

For small amplitudes everything works fine. But if the amplitude goes above a certain level the simulations aborts.
I tried some variations in the dynamicMeshDict such as changing from displacementLaplacian to displacementSBRStress or changing the diffusivity parameters
but the problem only changed the place place or time of occurrence.
In the attached pictures the problem is shown by a exemplarily case.
2 of the last four time steps are shown.
piston on the left, wall on the right side. colored by pressure

dynamicMeshDict:
Code:
dynamicFvMesh dynamicMotionSolverFvMesh;
motionSolverLibs ( "libfvMotionSolvers.so" );
solver displacementLaplacian;
displacementLaplacianCoeffs
{
      diffusivity quadratic inverseDistance (cylinder walls);
}
The test case is attached too, for anybody who may wants to give it a try.
(you have to execute blockMesh first)


Any help will be appreciated.
Best regards,
Henry
Attached Images
File Type: jpg mesh_time_0.0.jpg (192.1 KB, 5 views)
File Type: png p_displacementLaplacian_inverseDist_1.png (16.3 KB, 6 views)
File Type: png p_displacementLaplacian_inverseDist_4.png (17.0 KB, 7 views)
Attached Files
File Type: gz test_displacement.tar.gz (13.2 KB, 0 views)
pille is offline   Reply With Quote

Reply

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
Exporting mesh for an OpenFOAM chtMultiRegion case akrasemann enGrid 14 October 4, 2015 16:26
No layers in a small gap bobburnquist OpenFOAM Native Meshers: snappyHexMesh and Others 6 August 26, 2015 09:38
snappyhexmesh remove blockmesh geometry philipp1 OpenFOAM Running, Solving & CFD 2 December 12, 2014 11:58
snappyHexMesh won't work - zeros everywhere! sc298 OpenFOAM Native Meshers: snappyHexMesh and Others 2 March 27, 2011 21:11
fluent add additional zones for the mesh file SSL FLUENT 2 January 26, 2008 12:55


All times are GMT -4. The time now is 13:33.