|July 12, 2007, 13:13||
Compressible inviscid 1d spherical symmetry flow
I am solving 1d Euler spherical symmetry flow.
in other forms, rho_t+(rho u)_r=-alpha/r*(rho u) rhou_t+(rhou^2+p)_r=-alpha/r*(rho u^2) rhoe_t+(u(rhoe+p))_r=-alpha/r*(u(rhoe+p))
To solve this inhomogeneous hyperbolic equation, I used the fractional-step method for two equations
U_t+F_r=0 ! homogeneous conservation laws U_t=-alpha/r*F(U) ! ODE for the geometric source term
3rd RK and 2nd RK are employed for the first and the second equation above. I am thinking r is the distance from the orgin. That is, r is the same as the x in the cartesian system if the origin is located at the left-end.
Unfortunately, I couldn't obtain a desirable results without the oscillation.
I was trying to find any book containing this kind of problem but I couldn't find it.
Is there any misunderstanding about aboves? Have you seen any good reference how to discretize the geometric source term? Currently, I am using DGM method.
Please advice me. Thanks in advance!
|Thread||Thread Starter||Forum||Replies||Last Post|
|3-D Compressible Flow Boundary Conditions||Tyler||FLUENT||4||February 5, 2009 20:58|
|Convergence for inviscid flow over sylinder||Rjakk||Main CFD Forum||2||March 21, 2007 11:53|
|Can 'shock waves' occur in viscous fluid flows?||diaw||Main CFD Forum||104||February 16, 2006 06:44|
|Boundary Layer created by Euler Solvers||Jim||Main CFD Forum||31||November 18, 2001 00:18|
|compressible channel flow..||R.D.Prabhu||Main CFD Forum||0||July 17, 1998 17:23|