Context Navigation

← Previous Post
Next Post →

Adjusting Gravity in 3D Ablative RT module

Gravity value with different bottom heat flux for the current code

Heat flux at the bottom seems too low. Tried randomly bigger flux. Gravity increases first then drops. Seems something very wrong.

	flb=6.0E+21 (calculated)	flb=6.8E+21	flb=7.13E+21
gravity & total mass
Int(P+rho*V2)

Compare with 2D Case

1. Initial Profile

Initial profile are close enough except the momentum along the gravity direction.

Rho&T
py or pz

2. heat diffusion check

Turn hydro off. Compare the 2D (gravity along -y direction) and 3D (gravity along -z direction). The solver seems OK.

2D, flb=0	2D, flb=6E21
3D, flb=0	3D, flb=6E21

3. Pure hydro test

Turn off the heat diffusion off. Compare the density and the momentum along gravity direction -py for 2D and pz for 3D — plots are along the center line.

2D rho	2D, momentum	2D, pressure
3D rho	3D, momentum	3D, pressure

Here's a picture to compare the py (2D) and pz (3D) — both along the center line at time=1.345E-5 (cu)

4. compare top and bottom integral of P+Rho*V2

The integrals (pressure + rho*v2) were calculated with queries of integrate (2D) and weighted variable sum (3D) in visit.

Bottom	Top

Check the derivative of momentum

2D	3D

astrobear