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

hack - using porous region to fix local instability

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

Like Tree1Likes
  • 1 Post By ufocfd

LinkBack Thread Tools Search this Thread Display Modes
Old   August 14, 2021, 10:49
Default hack - using porous region to fix local instability
Giles Richardson
Join Date: Jun 2012
Location: Cambs UK
Posts: 84
Rep Power: 12
ufocfd is on a distinguished road
I recently had some dodgy cells which were causing some very spikey convergence and creating some flow instability in the region due to very high unphysical velocity at particular location. If you dont want to fix the mesh (for whatever reason) you can apply some porosity in the area to stabilise the flow. You just need to apply "topoSet - latestTime" (to create a cell zone), and then modify the fvOptions file (to apply some porosity to that cell zone). Its not an ideal approach but quite a useful "hack" if you need a quick fix sometimes. I found it gave an instant and significant improvement of the convergence and flow solution. I wonder if there are any other methods?

// radiator cell set
name radiator_cellSet;
type cellSet;
action new;
source sphereToCell;
centre (0.45 0.014 0);
radius 0.03;
// radiator cell zone
name radiator_cellZone;
type cellZoneSet;
action new;
source setToCellZone;
set radiator_cellSet;

type explicitPorositySource;
active yes;

muName mu;
type DarcyForchheimer;
selectionMode cellZone;
cellZone radiator_cellZone;

d d [0 -2 0 0 0 0 0] (1e10 1e10 1e10); // viscous
f f [0 -1 0 0 0 0 0] (0 0 0); // inertia

type cartesian;
origin (0 0 0);
type axesRotation;
e1 (1 0 0);
e3 (0 1 0);
sylvester likes this.
ufocfd is offline   Reply With Quote


convergence, hack, instability, mesh quality, porosity

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
Problem with chtMultiregionFoam radiation boundary condition baran_foam OpenFOAM Running, Solving & CFD 10 December 17, 2019 18:36
chtMultiRegionFoam- too slow sandymech1 OpenFOAM Running, Solving & CFD 4 November 20, 2017 19:51
HeatSource BC to the whole region in chtMultiRegionHeater xsa OpenFOAM Running, Solving & CFD 3 November 7, 2016 06:07
How to write k and epsilon before the abnormal end xiuying OpenFOAM Running, Solving & CFD 8 August 27, 2013 16:33
Could anybody help me see this error and give help liugx212 OpenFOAM Running, Solving & CFD 3 January 4, 2006 19:07

All times are GMT -4. The time now is 00:16.