CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

dynamic mesh refinement and rhoCentralFoam

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

LinkBack Thread Tools Search this Thread Display Modes
Old   January 16, 2013, 22:14
Default dynamic mesh refinement and rhoCentralFoam
Join Date: Aug 2012
Location: Calgary, Alberta, Canada
Posts: 77
Rep Power: 13
ChrisA is on a distinguished road
I'm trying to implement dynamic mesh refinement in my modified rhoCentralFoam solver to refine my mesh around shocks that form in my simulation. I believe I have set up the code properly (I used the tutorial for dynamic mesh implementation for dieselFoam) but its crashing when I run it. I believe this is due to not specifying the proper fluxes in the dynamicFvMesh dictionary correctFluxes entry. I notice theres a couple threads around here floating around with people who have done exactly what I am attempting to do. Which fluxes are required? Looking at the way rhoCentralFoam solves, most of the fields except phi and rhoU are declared each time the loop runs at the start of the loop and shouldn't need to be corrected if I refine my mesh at the end of the loop. I could however be mistaken as to how the dynamic meshing works.

Edit: So it looks like rho as a refine variable is a bad idea, it keeps tagging every cell for refinement up until the max number of cells. This is obviously not true. The case I'm running should dynamically refine to under 400,000 cells (work in FLUENT for a very similar case has given me roughly that many) and refining based on rho has gotten as high as 3million cells. Starting with a mesh of ~130,000 cells and using p as a refine field solves to completion without refining a single cell. Anyone with some knowledge regarding dynamic meshing and rhoCentralFoam able to give some advice?

Edit2: I think I've figured out a bit of how the mesh dict works, at which point the above comments sound fairly silly. It only refines cells with the field value between the two values indicated in the dict. I thought it was a bit more clever than that and worked off the relative values between neighbouring cells or something of the sort. Which means I need to code something that will establish that and refine based off that. Again, advice from those who have done this in rhoCentralFoam to accurately resolve shocks would be appreciated.

Last edited by ChrisA; January 17, 2013 at 20:07.
ChrisA is offline   Reply With Quote

Old   March 21, 2013, 08:00
Join Date: Nov 2012
Posts: 83
Rep Power: 13
Henning86 is on a distinguished road
did you make any progress with refinement around shocks?

It would be great if you could post the solver with a simple example (e.g. forward step).
Henning86 is offline   Reply With Quote


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
rhoCentralFoam and wallHeatFlux RomW OpenFOAM Post-Processing 8 November 8, 2012 10:22
Mesh refinement in parallel KrisT OpenFOAM Running, Solving & CFD 1 July 4, 2010 02:38

All times are GMT -4. The time now is 19:41.