CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Native Meshers: blockMesh

How to mesh this?

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

Reply
 
LinkBack Thread Tools Display Modes
Old   November 17, 2014, 21:19
Default How to mesh this?
  #1
Member
 
Tony Ladd
Join Date: Aug 2013
Posts: 34
Rep Power: 3
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.

Thanks

Tony
Attached Images
File Type: jpg C.jpg (12.2 KB, 20 views)
File Type: jpg corner.jpg (59.7 KB, 15 views)
File Type: jpg tip.jpg (70.3 KB, 19 views)
tladd is offline   Reply With Quote

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

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

Old   November 18, 2014, 08:15
Default
  #3
Member
 
Tony Ladd
Join Date: Aug 2013
Posts: 34
Rep Power: 3
tladd is on a distinguished road
Paulo

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.

Tony
tladd 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
sliding mesh problem in CFX Saima CFX 35 May 12, 2014 20:28
Star CCM Overset Mesh Error (Rotating Turbine) thezack CD-adapco 3 December 11, 2013 04:09
[ICEM] surface mesh merging problem everest ANSYS Meshing & Geometry 39 June 5, 2013 19:02
[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 08:17.