CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

DarcyForchheimer crashes when sudden changes in 'D'

Register Blogs Community New Posts Updated Threads Search

Like Tree2Likes
  • 1 Post By sylvester

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 10, 2023, 22:42
Talking DarcyForchheimer crashes when sudden changes in 'D'
  #1
Senior Member
 
Desh
Join Date: Mar 2021
Location: Sydney
Posts: 118
Rep Power: 5
dasith0001 is on a distinguished road
Hi Foarmers,

I am running a chtmultiRegionFoam and there is a flow channel that I want to close and open time to time. And I am using OF10.

After bit of survey, I thought it would be easiest to use ''explicitPorositySource'' in the fluid region with high resistance values of 'D' (D=1e20) when I want to block the flow and ''D=0'' when I want the flow to go freely.

The D=1e20 stops the flow effectively and D=0 opens it and I can clearly see the changes in the velocity profile when I manually change the D value.

Problem now is the solver becomes so unstable (crashes ) when I suddenly change the 'D' value from 1e20 to 0. Any ideas resolving the issue?

Is there a way to change the D gradually with respect to time ?

Thank you
Dasith
dasith0001 is offline   Reply With Quote

Old   July 12, 2023, 07:36
Default timeActivatedFileUpdate
  #2
Member
 
Roland
Join Date: Mar 2009
Location: Netherlands
Posts: 93
Rep Power: 17
sylvester is on a distinguished road
timeActivatedFileUpdate?
https://github.com/OpenFOAM/OpenFOAM...edFileUpdate.H

(BTW, 1e20 seems a bit excessive to me. A lower value might produce similar results with a more stable simulation.)
dasith0001 likes this.
sylvester is offline   Reply With Quote

Old   July 12, 2023, 21:33
Default
  #3
Senior Member
 
Desh
Join Date: Mar 2021
Location: Sydney
Posts: 118
Rep Power: 5
dasith0001 is on a distinguished road
Quote:
Originally Posted by sylvester View Post
timeActivatedFileUpdate?
https://github.com/OpenFOAM/OpenFOAM...edFileUpdate.H

(BTW, 1e20 seems a bit excessive to me. A lower value might produce similar results with a more stable simulation.)
Hi Ronald,

Thank you very much, works perfectly well now. if anyone interested, following is what I used in the functions (in controlDict )

PHP Code:
    fileUpdate1
    
   
{
        
type                       timeActivatedFileUpdate;
        
functionObjectLibs         ("libutilityFunctionObjects.so");
        
outputControl              timeStep;
        
outputInterval             20;
        
fileToUpdate               "$FOAM_CASE/constant/fluidRegion/fvModels";
        
timeVsFile
        
(
            (
6220    "$FOAM_CASE/constant/fluidRegion/fvModels_1")
            (
6240    "$FOAM_CASE/constant/fluidRegion/fvModels_2")
            (
6260    "$FOAM_CASE/constant/fluidRegion/fvModels_3")
        );
    } 
Cheers
dasith0001 is offline   Reply With Quote

Reply

Tags
darcyforchheimer, explicitporositysource, unstable behaviour


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
Need help how to meshing in sudden contraction and sudden expansion 2D geometry pipe Khai123 Main CFD Forum 2 April 24, 2022 04:15
fvOptions- DarcyForchheimer Modification us7 OpenFOAM Running, Solving & CFD 3 August 6, 2020 04:35
[ICEM] best blocking strategy for sudden expansion pulplus ANSYS Meshing & Geometry 3 August 3, 2020 12:46
LES, Courant Number, Crash, Sudden Alhasan OpenFOAM Running, Solving & CFD 5 November 22, 2019 03:05
flo-efd v11.0.0 crashes YoavF FloEFD, FloWorks & FloTHERM 3 June 21, 2012 13:37


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