|
[Sponsors] |
November 17, 2014, 20:19 |
How to mesh this?
|
#1 |
Member
Tony Ladd
Join Date: Aug 2013
Posts: 48
Rep Power: 12 |
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 |
|
November 18, 2014, 05:50 |
|
#2 |
Senior Member
Paulo Vatavuk
Join Date: Mar 2009
Location: Campinas, Brasil
Posts: 196
Rep Power: 17 |
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 05:51. |
|
November 18, 2014, 07:15 |
|
#3 |
Member
Tony Ladd
Join Date: Aug 2013
Posts: 48
Rep Power: 12 |
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 |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
sliding mesh problem in CFX | Saima | CFX | 46 | September 11, 2021 07:38 |
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 |
[ICEM] Problem making structural mesh on a surface | froztbear | ANSYS Meshing & Geometry | 1 | November 10, 2011 08:52 |
[snappyHexMesh] snappyHexMesh won't work - zeros everywhere! | sc298 | OpenFOAM Meshing & Mesh Conversion | 2 | March 27, 2011 21:11 |