|
[Sponsors] |
Temperature instability in transonic simulation |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
March 27, 2021, 08:43 |
Temperature instability in transonic simulation
|
#1 |
New Member
Michelangelo
Join Date: Apr 2015
Posts: 17
Rep Power: 11 |
Hi Fomers,
I have a problem with a transonic simulation. I’m using rhoPimpleFoam (with "transonic yes") and I’m progressively reaching the outlet target solution. However when the solution is close to the sonic condition, an instability of temperature field appears as shown in the attached images. Has anyone already found this instability? |
|
March 29, 2021, 04:11 |
|
#2 |
Senior Member
Join Date: Dec 2019
Posts: 215
Rep Power: 7 |
Can you share your fvSchemes and controldict?
I suppose you are using high order schemes with a big timestep, thats why those oscillations might occure. |
|
March 29, 2021, 07:21 |
|
#3 |
New Member
Michelangelo
Join Date: Apr 2015
Posts: 17
Rep Power: 11 |
Hi shock77,
Code:
FoamFile { version 2.0; format ascii; class dictionary; object fvSchemes; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // ddtSchemes { default Euler; } gradSchemes { default Gauss linear; grad(U) cellLimited Gauss linear 1; } divSchemes { default none; div(phi,U) Gauss limitedLinearV 1; turbulence Gauss upwind; energy Gauss limitedLinear 1; div(phi,k) $turbulence; div(phi,omega) $turbulence; div(phi,e) $energy; div(phi,h) $energy; div(phi,K) $energy; div(phi,Ekp) $energy; div(phid,p) Gauss limitedLinear 1; div(phiv,p) Gauss limitedLinear 1; div((phi|interpolate(rho)),p) Gauss linear; div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear; }laplacianSchemes { default Gauss linear corrected; } interpolationSchemes { default linear; } snGradSchemes { default corrected; } wallDist { method meshWave; } and Code:
FoamFile { version 2.0; format ascii; class dictionary; location "system"; object controlDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // application rhoPimpleFoam; startFrom latestTime; startTime 0; stopAt endTime; endTime 0.5; deltaT 1e-5; writeControl adjustable; writeInterval 0.001; purgeWrite 0; writeFormat ascii; writePrecision 10; writeCompression off; timeFormat general; timePrecision 9; runTimeModifiable true; adjustTimeStep yes; maxCo 1; maxDeltaT 1e-3; minDeltaT 1e-6; Yesterday I tested Gauss limitedVanLeer -2.0 3.0; for $energy and it improved the instability but it crashed anyway near the sonic condition. |
|
March 29, 2021, 07:57 |
|
#4 |
Senior Member
Join Date: Dec 2019
Posts: 215
Rep Power: 7 |
First, I think you should try that and see if it helps:
Code:
gradSchemes { default cellLimited Gauss linear 1; } |
|
March 29, 2021, 16:26 |
|
#5 |
New Member
Michelangelo
Join Date: Apr 2015
Posts: 17
Rep Power: 11 |
Thanks,
now I try as you advised me. Updates soon. |
|
April 1, 2021, 16:10 |
|
#6 |
New Member
Michelangelo
Join Date: Apr 2015
Posts: 17
Rep Power: 11 |
I tried with cellLimited but the problem is still there.
I'll try other setups. Thanks anyway for the advice. |
|
April 1, 2021, 16:21 |
|
#7 |
Senior Member
Join Date: Dec 2019
Posts: 215
Rep Power: 7 |
You should try to run your simulation first order first. Than step by step increase the order.
|
|
Tags |
rhopimplefoam, temperature instability, transonic |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[openSmoke] libOpenSMOKE | Tobi | OpenFOAM Community Contributions | 562 | January 25, 2023 09:21 |
Surface Source - Fixed Temperature? | robtheslob | FloEFD, FloWorks & FloTHERM | 18 | May 12, 2017 02:28 |
simulation of particle temperature with own routine | amd18 | CFX | 0 | March 27, 2017 05:41 |
ThermophysicalProperties do not vary with temperature along the simulation | DuarteMagalhaes | OpenFOAM Running, Solving & CFD | 0 | July 10, 2014 06:08 |
how to set up high temperature gas turbine flow simulation? | adam2008 | CFX | 1 | July 22, 2009 18:33 |