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

Dynamic Mesh - Strange Pressure behaviour in the corner

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 2, 2024, 06:51
Smile Dynamic Mesh - Strange Pressure behaviour in the corner
  #1
New Member
 
Join Date: Jan 2024
Posts: 8
Rep Power: 2
danylo is on a distinguished road
Hello Foamers!

Hope everyone is doing well.

I am trying to model a phase change of water to ice, which includes an expansion of my system due to the density difference of the phases. The cooling comes from the wall on the right side and the idea is, that expansion is inversely proportional to ice present in the top cell of a column, so frozen regions dont expand any further and my system is bulge shaped in the end. The expansion is realized by applying a calculated Point-Displacement-value to the top boundary. I attached screenshots, where you can see the beginning of this happening.

Now my problem is, that once the first frozen region slows its expansion down, the pressure field shows a weird behaviour in the corner. The cell in the corner has a higher pressure than the system (1e5), while the one right beneath it has a lower one. This trend accelerates as the simulation goes on and leads to a crash. And the pressure field is not the only one showing weird behavior due to the mesh motion (e.g. ice fraction tends to go over 100 % in this case). I guess the reason for this is, that the CellMotion, which OpenFOAM calculates, tends to zero in the top right corner, but the cells beneath it still have a relevant CellMotion value due to the moving neighbour cells (Screenshot of this is attached), but I dont know how I could counteract this.

My question is, what are the internal mechanisms of the OpenFOAM mesh-motion, which could lead to this an what can I do to prevent it from happening? Any help would be appreciated!

Some extra Info on my case:
BCs:
U
Top: pressureInletOutletVelocity (0 0 0)
Rest: noSlip

p
Top: totalPressure 1e5
Rest: zeroGradient

motionSolver:
velocityComponentLaplacian

diffusivity:
quadratic inverse distance to bottom wall

Down below are screenshots of the CellMotion in y direction, Ice Fraction and temperature.

Thank you in advance!
Attached Images
File Type: png CellMotionUy.png (47.3 KB, 9 views)
File Type: png IceFraction.png (20.5 KB, 7 views)
File Type: png Pressure.png (19.2 KB, 5 views)
danylo is offline   Reply With Quote

Old   April 2, 2024, 07:00
Arrow Addition
  #2
New Member
 
Join Date: Jan 2024
Posts: 8
Rep Power: 2
danylo is on a distinguished road
Addition:

I am using a customized version of rhoPimpleDyMFoam on OpenFOAM 5.
danylo is offline   Reply With Quote

Old   April 8, 2024, 04:09
Default Update
  #3
New Member
 
Join Date: Jan 2024
Posts: 8
Rep Power: 2
danylo is on a distinguished road
Update on my question

I improved the code insofar as the equations work, when the Mesh Motion is distributed evenly across the Mesh. To come by this problem of having a CellMotion in the middle of the cell column and none at the top, I now try to work with directional diffusivity, where I disable the diffusivity in x and z direction. So the line in dynamicMeshDict looks like this:

diffusivity: directional (0 1 0);
(I attached the rest of dynamicMeshDict as a screenshot)

In the beginning this works out nicely and I have next to no CellMotion in the Frozen Regions. But after a while the the internal cells start to have some strange Motion, while the top layer has none at all. So it seems like there still is a diffusivity in the x direction even though it should not be there. Some cells even move with a negative velocity which should not be possible at all (Either the whole mesh expands or the whole mesh contracts) and there are opposing movements in the same column. This skews my mesh and leads to a solver crash once again

I attached screenshots of how that looks.

Does anybody have an idea what could cause this?
Attached Images
File Type: png Cells tilted from right to left.png (38.8 KB, 3 views)
File Type: png WeirdCellMotion.png (57.0 KB, 2 views)
File Type: png DynamicMeshDict.png (111.2 KB, 2 views)

Last edited by danylo; April 8, 2024 at 04:17. Reason: added attachments
danylo is offline   Reply With Quote

Reply

Tags
corner, dynamic mesh, expansion, phase change, pressure

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
foam-extend-4.1 release hjasak OpenFOAM Announcements from Other Sources 19 July 16, 2021 05:02
Fail to converge when solving with a fabricated solution zizhou FLUENT 0 March 22, 2021 06:33
total pressure dynamic pressure static pressure roi247 FLUENT 19 February 4, 2020 12:52
[snappyHexMesh] No layers in a small gap bobburnquist OpenFOAM Meshing & Mesh Conversion 6 August 26, 2015 09:38
Hydrostatic pressure in 2-phase flow modeling (CFX4.2) HB &DS CFX 0 January 9, 2000 13:19


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