|February 9, 2006, 02:29||
Hi everybody, I have tried t
Join Date: Mar 2009
Posts: 24Rep Power: 6
I have tried to manage a viscous version of rhopSonicFoam to solve high-Mach number flows.
Since the equations are written in conservative form, the viscous terms have to be formulated in an explicit manner (as far as I know):
+fv::gaussConvectionScheme<vector>(mesh, phiv, rhoUScheme).fvmDiv(phiv, rhoU)
solve(rhoUEqn == - fvc::grad(p)
+ fvc::laplacian(mu, U)
+ mvConvection.fvmDiv(phiv, rhoE)
solve(rhoEEqn == - mvConvection.fvcDiv(phiv, p)
+ fvc::laplacian(mu * Cp / Pr, T)
Question #1: is this the only way to include them in a conservative form of N-S equations?
The computation of a bluff body at Mach 15 is performed easily, and a grid-independent solution is reached. On the other hand, a flat-plate parallel to the flow (Mach 9.9) shows instabilities at the first cells as soon as the grid is refined (gradients are sharper).
I have also tried to increase the nOuterCorrectors and nCorrectors but the solution diverges within the timestep!
Question #2: any idea (or reference) about how to solve the problem?
Question #3: have you references where I can get a description of the multivariateGaussConvectionScheme? I would like to know how it works and why it is used in this solver...
|Thread||Thread Starter||Forum||Replies||Last Post|
|SonicTurbFoam Vs rhopSonicFoam||vishal||OpenFOAM Running, Solving & CFD||0||January 29, 2009 00:58|
|Canbt run rhopSonicFoam from the tutorials for version 13||rinaoldek||OpenFOAM Running, Solving & CFD||0||January 11, 2009 04:06|
|How do about total temperature BC in rhopSonicFoam||yijin_li||OpenFOAM Running, Solving & CFD||0||December 3, 2008 01:34|
|RhopSonicFoam reference||abhishek||OpenFOAM Running, Solving & CFD||0||November 5, 2006 22:43|
|RhopSonicFoam||didomenico||OpenFOAM Running, Solving & CFD||0||October 3, 2005 10:53|