|
[Sponsors] |
Wall Temperature in Nozzle with RhoCentralFoam |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
February 10, 2016, 20:16 |
Wall Temperature in Nozzle with RhoCentralFoam
|
#1 |
Member
Hilbert
Join Date: Aug 2015
Location: Australia
Posts: 50
Rep Power: 11 |
Dear Foamers,
I have been working on a nozzle simulation with rhoCentralFoam. I have a CFX solution which I am currently using as a baseline. The core-flow matches up beautifully with the CFX results. The wall pressure is also looks very nice, and compares very well. My current issue is with the wall temperature. I have a large drop in wall temperature in my throat, which recovers after my throat to about 3 degrees lower then my CFX baseline. My temperature B.C at the well is zero gradient, i.e adiabatic wall. If I compute the adiabatic wall temperature for my case it look like the CFX solution is right. The wall shear is also over predicted w.r.t the CFX solution. Which would suggest an over prediction of the velocity wall gradient. Yplus is of order 1, and checkmesh also checks out. My fvSchemes file is currently set to the following: Code:
fluxScheme Tadmor; ddtSchemes { default Euler; } gradSchemes { default leastSquares; grad(U) cellLimited leastSquares 1; grad(rho) cellLimited leastSquares 1; grad(omega) cellLimited leastSquares 1; grad(k) cellLimited leastSquares 1; } divSchemes { default none; div(tauMC) Gauss linear; div(phi,U) Gauss linear; div(phid,p) Gauss linear; div(phi,K) Gauss linear; div(phi,h) Gauss Gamma 1; div(phi,k) Gauss Gamma 1; div(phi,epsilon) Gauss linear; div(phi,R) Gauss Gamma 1; div(phi,omega) Gauss Gamma 1; div((rho*R)) Gauss Gamma 1; div(R) Gauss Gamma 1; div(U) Gauss linear; div((muEff*dev2(T(grad(U))))) Gauss linear; div(yPhi,yWall) Gauss linear ; } laplacianSchemes { default none; laplacian(muEff,U) Gauss Gamma phi 1 corrected; laplacian(alphaEff,e) Gauss Gamma phi 1 corrected; laplacian((rho*DomegaEff),omega) Gauss Gamma phi 1 corrected; laplacian((rho*DkEff),k) Gauss Gamma phi 1 corrected; laplacian(yPsi) Gauss Gamma phi 1 corrected; laplacian(yWall) Gauss Gamma phi 1 corrected; } interpolationSchemes { default linear; reconstruct(rho) Gamma 1; reconstruct(U) GammaV 1; reconstruct(T) Gamma 1; interpolate(rhoU) Gamma pos 1; interpolate(rho) Gamma pos 1; } snGradSchemes { default corrected; } wallDist { method advectionDiffusion; // Optional entry enabling the calculation // of the normal-to-wall field nRequired false; advectionDiffusionCoeffs { method Poisson; epsilon 0.4; tolerance 1e-3; maxIter 20; } } |
|
August 25, 2024, 09:03 |
|
#2 |
New Member
Gaoqiang Yang
Join Date: Apr 2018
Posts: 14
Rep Power: 8 |
Hi, have you solved this problem?
I have met same problems, i.e. pressure and temperature seems correct at the center line of nozzle. However, the wall temperature is not correct, which is increasing along the flow direction (It should be decreasing ) |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[Commercial meshers] Fluent3DMeshToFoam | simvun | OpenFOAM Meshing & Mesh Conversion | 50 | January 19, 2020 16:33 |
Water subcooled boiling | Attesz | CFX | 7 | January 5, 2013 04:32 |
How to adjust wall temperature | nutow | FLUENT | 0 | September 4, 2011 06:37 |
Errors running allwmake in OpenFOAM141dev with WM_COMPILE_OPTION%3ddebug | unoder | OpenFOAM Installation | 11 | January 30, 2008 21:30 |
About the wall temperature | zt_zt | FLUENT | 3 | May 12, 2003 09:09 |