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

Version/Machine dependent stability?

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 4, 2019, 07:10
Default Version/Machine dependent stability?
  #1
Member
 
Piotr Ładyński
Join Date: Apr 2017
Posts: 55
Rep Power: 9
piotr.mecht is on a distinguished road
I prepared some simulation setups for buoyantPimpleFoam which worked perfectly fine whith every scenario so far on my home PC, which runs OpenFOAM v1906 Windows10-native on six core 9th gen Intel processor (I used variable number of cores). I moved my simulation files on some offhand PC at my workplace, so I could do other stuff at home. This machine runs Openfoam 7 with Linux Mint on some older 4 core Intel Core 2 Quad Q8200 processor. I had issues with correct docker setup (i'll work on this when i'll find some more time for this), so i switched to foundation version for easy installation and usage.


My simulation quickly goes unstable after first iterations on the second machine with the same files that worked on the first machine. When I reduced time step vale the max Courant number grows even quicker.
I've seen this symptom before, while using icoPolynomial density model with pressure-work included to energy equation (sensibleEnthalpy). I disabled this with:
Code:
dpdt off;
in thermophysicalProperties dictionary and this worked well on my home machine. Is this input unrecognized with foundation version maybe?


What shoud I do to expect the same output?

Update:
I've managed to install Docker on Linux machine later taht day and the simulation in v1906 works like on my home computer now, so it's definitely .org/.com build dependent. Despite this, I'm curious about differences to my input, that would possibly produce the same results. I'm not convinced to use OF commads in separate container as it's harder to navigate between e.g. the text editor and other stuff.
Attached Files
File Type: zip dictionariesOnly.zip (20.2 KB, 0 views)

Last edited by piotr.mecht; November 4, 2019 at 14:07. Reason: Update
piotr.mecht is offline   Reply With Quote

Old   November 5, 2019, 05:35
Default
  #2
Senior Member
 
Joachim Herb
Join Date: Sep 2010
Posts: 650
Rep Power: 22
jherb is on a distinguished road
I see one major difference in the solvers between the Foundation and the ESI (OpenFOAM-plus) development lines:


In the pEqn.H of the ESI version, the following lines have been added (assuming, that you have a boundary condition on at least one patch with a fixed pressure

Code:
thermo.correctRho(psi*p - psip0);
It would be really interesting, what happens, if you remove this line form the solver, recompile it and test it:
https://develop.openfoam.com/Develop...am/pEqn.H#L100
jherb is offline   Reply With Quote

Reply


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
Dependent variables for linearized source terms ASimonsen Fluent UDF and Scheme Programming 0 April 5, 2017 03:20
Time dependent and geometry dependent boundary conditions in ACUSOLVE arsroy Main CFD Forum 1 December 28, 2016 08:32
Question Regarding Stability of Nonlinear Terms H0T_S0UP Main CFD Forum 1 January 3, 2014 04:29
when I use multi block the condition of stability can not be fulfilled mostafa_khan FLOW-3D 2 April 19, 2012 13:17
Solving linear stability eq. K S Chang Main CFD Forum 3 January 19, 2004 17:01


All times are GMT -4. The time now is 01:13.