CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Native Meshers: blockMesh

How to mesh this?

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

LinkBack Thread Tools Display Modes
Old   November 17, 2014, 21:19
Default How to mesh this?
Tony Ladd
Join Date: Aug 2013
Posts: 37
Rep Power: 6
tladd is on a distinguished road
I have been running some simulations of a Stefan problem. Basically two domains with a Laplace equation in each domain and a jump condition at the interface to include the (step) change in permeability. We solve a convection diffusion equation for concentration (C) in the inner region and advect the interface according to snGrad(C). The steady-state (more or less) finger shape at Pe = 50 (with the jump in perm = 10) is shown in C.jpg. There are 256^2 zones in this region - too many to show; coarser meshes develop artefacts before the steady state is reached.

Problem is the mesh is not very good; there are many non-orthogonal faces - see corner.jpg and tip.jpg. However, the max angle is still less than 80 deg. I make the mesh using blockMesh with a single block and the boundary defined with polyLine (OF2.3.0). I can see I might do better with foamyQuadMesh but before that I was wondering if anyone could suggest a better strategy for blockMesh.


Attached Images
File Type: jpg C.jpg (12.2 KB, 24 views)
File Type: jpg corner.jpg (59.7 KB, 18 views)
File Type: jpg tip.jpg (70.3 KB, 21 views)
tladd is offline   Reply With Quote

Old   November 18, 2014, 06:50
Senior Member
Paulo Vatavuk
Join Date: Mar 2009
Location: Campinas, Brasil
Posts: 162
Rep Power: 10
vatavuk is on a distinguished road
Hi Tony,
If have time (and patience) you could try this topology with 10 blocks.
Best Regards,

Last edited by vatavuk; November 21, 2014 at 06:51.
vatavuk is offline   Reply With Quote

Old   November 18, 2014, 08:15
Tony Ladd
Join Date: Aug 2013
Posts: 37
Rep Power: 6
tladd is on a distinguished road

Thanks for the suggestion. It is a bit or work and I think I will have to remesh as I go. I will let you you know if I have any success. But its clever how you keep the hex blocks - I did not see that.

tladd 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
Star CCM Overset Mesh Error (Rotating Turbine) thezack Siemens 7 October 12, 2016 11:14
[ICEM] surface mesh merging problem everest ANSYS Meshing & Geometry 44 April 14, 2016 06:41
sliding mesh problem in CFX Saima CFX 45 September 22, 2015 10:53
[ICEM] Problem making structural mesh on a surface froztbear ANSYS Meshing & Geometry 1 November 10, 2011 09:52
snappyHexMesh won't work - zeros everywhere! sc298 OpenFOAM Native Meshers: snappyHexMesh and Others 2 March 27, 2011 21:11

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