|
[Sponsors] |
transientSimpleFoam - SIMPLE with time derivatives |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
December 14, 2015, 11:21 |
transientSimpleFoam - SIMPLE with time derivatives
|
#1 |
Senior Member
|
The explanation in the top of the source file in this solver says " For non-Newtonian fluid'' which questions me. Was this Solver really built for Non-Newtonian fluids Bruno?
[ Moderator note: Moved from this thread: http://www.cfd-online.com/Forums/ope...implefoam.html ] Last edited by wyldckat; December 28, 2015 at 16:25. Reason: see "Moderator note:" |
|
December 26, 2015, 14:01 |
transientSimpleFoam - SIMPLE with time derivatives
|
#2 |
Senior Member
|
Hello Everyone;
Hope every one is doing great ! Has any one used this solver which is supposed to be SIMPLE but with time derivatives so ( unsteady simple foam). My only concern is it says that it is for non- newtonian flow in the source file. Is this true? Here is the link: https://github.com/wyldckat/transien...ientSimpleFoam I saw a post where Prof Jasak said that transientSimple algorithm can be used for large courant number flows. So I am really keen to use this solver. Furthermore has any one coded SIMPLEC or SIMPLER for transient ? Regards and have a good day!! Shereez |
|
December 26, 2015, 23:13 |
|
#3 |
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,912
Rep Power: 36 |
You can use pimpleFoam instead.
__________________
Alberto Passalacqua GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats (current status: http://albertopassalacqua.com/?p=1541) OpenQBMM - An open-source implementation of quadrature-based moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using. |
|
December 27, 2015, 11:05 |
|
#4 |
Senior Member
|
Dear alberto;
Many thanks for the reply. I have had previous tries trying to run pimpleFoam for cases. The thing is I am running airfoil vortex shedding cases and depending on the strouhal number some of the Time periods are 0.2 seconds, 1 seconds or even 5 seconds. And pimpleFoam stability blows up above any courant number 20 not even menioning about the accuracy of the solution. however the transient Simple Foam I have compiled can handle and converge any courant number. My only concern is that it says in the source file description that it is for non- newtonian flows which I don't understand why. But once again, many thanks for your suggestion. Regards Shereez |
|
December 28, 2015, 04:00 |
|
#5 |
Member
W.T
Join Date: Oct 2012
Posts: 35
Rep Power: 13 |
In transportProperties file yuo can chose any single phase transport model, including nonNewtonian and IMO thats is a reason.
|
|
December 28, 2015, 13:56 |
|
#6 | |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,974
Blog Entries: 45
Rep Power: 128 |
Greetings to all!
@Shereez: Quote:
Code:
singlePhaseTransportModel laminarTransport(U, phi); //autoPtr<incompressible::RASModel> turbulence autoPtr<incompressible::turbulenceModel> turbulence ( //incompressible::RASModel::New(U, phi, laminarTransport) incompressible::turbulenceModel::New(U, phi, laminarTransport) ); How do I know it is the file "transportProperties" that is loaded? Check the file pointed out by this command: Code:
echo $FOAM_SRC/transportModels/incompressible/singlePhaseTransportModel/singlePhaseTransportModel.C Best regards, Bruno
__________________
Last edited by wyldckat; December 28, 2015 at 16:26. Reason: edited the header to greet everyone, not just Shereez |
||
December 28, 2015, 14:31 |
|
#7 | |
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,912
Rep Power: 36 |
Quote:
For what I know transientSimpleFoam is not part of OpenFOAM either...
__________________
Alberto Passalacqua GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats (current status: http://albertopassalacqua.com/?p=1541) OpenQBMM - An open-source implementation of quadrature-based moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using. |
||
December 28, 2015, 16:20 |
|
#8 | |
Senior Member
|
Quote:
But for instance consider a case like this : A vortex shedding case for which Force convergence has a periodic solution which has a Frequency of 4Hz : --> so Time period = 1/4 = 0.25 seconds. Let's assume that this case is for reynolds number = 1 million. Which means to have a courant number of 1 or below I will need a DT = 1e-05 seconds or lower. So one complete wave length of the force curve can be computed using 20,000 time steps. Which is computationally expensive. One of my professors pointed out to me that if I can place 100 or 200 points in one time period ( One wave length) then this should lead to a good estimate of the periodic shedding pattern that we are interested in. So I believe if I can have a DT = 0.002 seconds I should still be able to capture the unsteadiness that I am interested in. Or not Well, I will see. Best Regards Shereez |
||
December 29, 2015, 08:20 |
|
#9 | |
Senior Member
|
Quote:
Bruno; I did not see your reply until now. that's great news. thanks very much Regards shereez |
||
August 24, 2018, 16:54 |
|
#10 |
Senior Member
Guilherme
Join Date: Apr 2017
Posts: 216
Rep Power: 9 |
Hi,
I edited the original code for OpenFOAM 5.0. Is attached. I kindly ask them to verify the code (Bruno?), since I did not test them and feedback would be extremely important. |
|
August 24, 2018, 18:26 |
|
#11 | |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,974
Blog Entries: 45
Rep Power: 128 |
Quick answer: Sorry, but... I already provided the port of this utility for OpenFOAM 5, roughly 30 minutes before your post above, also as an answer to your associated request on that thread:
Quote:
|
||
August 25, 2018, 08:31 |
|
#12 | |
Senior Member
Guilherme
Join Date: Apr 2017
Posts: 216
Rep Power: 9 |
Quote:
About the code, I realized that it's a mix between SIMPLE and PISO codes, correct?! Would it be possible to create a SIMPLE(C) with discretization of time?! Without the PISO loop... |
||
August 25, 2018, 09:14 |
|
#13 | |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,974
Blog Entries: 45
Rep Power: 128 |
Quote:
You will have to study the solvers in more detail to see that for yourself, as well as look for information online about SIMPLEC in general. |
||
August 26, 2018, 12:33 |
|
#14 | |
Senior Member
Guilherme
Join Date: Apr 2017
Posts: 216
Rep Power: 9 |
Quote:
I have a question and maybe you can help me understand,... I was studying both codes (transientSimpleFoam and simpleFoam) and both contain the term ''p.storePrevIter( )''. The code (simpleFoam) compiled in current versions of OF does not contain this term... but I think they are included, as this would misrepresent the model of both. Could you tell where this variable is in the current codes? |
||
September 1, 2018, 15:08 |
|
#15 |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,974
Blog Entries: 45
Rep Power: 128 |
Quick answer:
|
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[Other] refineWallLayer Error | Yuby | OpenFOAM Meshing & Mesh Conversion | 2 | November 11, 2021 12:04 |
decomposePar problem: Cell 0contains face labels out of range | vaina74 | OpenFOAM Pre-Processing | 37 | July 20, 2020 06:38 |
Floating point exception error | lpz_michele | OpenFOAM Running, Solving & CFD | 53 | October 19, 2015 03:50 |
pisoFoam with k-epsilon turb blows up - Some questions | Heroic | OpenFOAM Running, Solving & CFD | 26 | December 17, 2012 04:34 |
calculation diverge after continue to run | zhajingjing | OpenFOAM | 0 | April 28, 2010 05:35 |