My colleagues and I are trying to simulate jets emitted inside
.the envelope of a Red Giant Star
The shocks caused by the expansion of the jets create dense
and hot fronts, where the radiation pressure, (a*T^4)/3 and the
,radiation energy density
.e=a*T^4 are dominant over the gas pressure
?Has anyone implemented these terms in the code in the past
Ron