|October 6, 2013, 13:06||
Inverting speed with Godunov solver
Join Date: Sep 2013
Posts: 6Rep Power: 5
I have written a simple Godunov upwind scheme and a HLLC solver to study an astrophysical problem. I tested it with the Sod test and it works smoothly.
Then, I change the initial conditions to the astrophysical ones of my interest and still works fine, until something happens. The problem has spherical symmetry and the gas initially moves with increasing radial distance from the star. As a consequence of gravity, the gas slows down and eventually inverts the motion. That is, the velocity is positive at the beginning and negative after some time. At the inverting point, the program crashes. I suppose that I failed to implement correctly the Godunov flux, writing it only in the positive velocity case. Could you give me a reference where I can find a clear description of this issue?
Thank you very much
|numerical methods, numerical scheme|
|Thread||Thread Starter||Forum||Replies||Last Post|
|Creating New Solver: For particle-laden compressible jets||sankarv||OpenFOAM Running, Solving & CFD||17||December 3, 2014 20:41|
|OpenCL linear solver for OpenFoam 1.7 (alpha) will come out very soon||qinmaple||OpenFOAM Announcements from Other Sources||4||August 10, 2012 11:00|
|Creating New Solver: For particle-laden compressible jets||sankarv||OpenFOAM||0||April 4, 2010 18:06|
|Problem with implicit unsteady solver||CCMuser||STAR-CCM+||2||March 3, 2010 12:20|
|Does it metter to increase RAM for solver speed??||raima||Main CFD Forum||1||February 28, 2008 11:47|