Context Navigation

← Previous Post
Next Post →

OutflowWind: spline soft for Point gravity

This is to try to understand/fix the tight sonic surface found in 3D co-rotating frame with omega=0 as shown in the 3rd part of blog:bliu06112015_2 and the different line plots for the 3D co-rotating and 2.5D as shown in the 4th part of the same blog post… In 2.5D and 3D simulations before the problem module hardcoded PLUMERSOFT and soft radius=1 for the pointgravity.

        SPLINESOFT = 1, & !g ~ r/r^2 for r < r_soft and then goes to 0 as r-> 0 
        PLUMMERSOFT = 2 !g ~ (r^2+r_soft^2)^(-3/2) r

1. Code modification: Splinesoft

Since in the 3D case we use r=1 but in 2.5D we use r=2. So the g in 3D is about twice of the g in 2.5D.

To handle point gravity more properly, we change the code to use SPLINESOFT instead. Here's the new update to the OutflowWind module

        outflow_radius=0 ! outflow radius is 0 always
        outflow_thickness=planet_radius 
        pointgravity_soft=1 !splinesoft always
        pointgravity_r_soft=0.5*planet_radius

2. Testing Results with Splinesoft =

	New Results with SplineSoft (low res)	Old Results with PlummerSoft (high res)
2.5D no stellar wind	;movie 8 zones per radii	;high res movie
2.5D stellar wind	;movie 8 zones per radii	;high res movie
3D Co-rotating	; density 5 to 10 zones per radii	; density with 40 to 80 zones per radii
3D Co-rot and 2D plots compare	;

astrobear

Context Navigation

OutflowWind: spline soft for Point gravity

1. Code modification: Splinesoft

2. Testing Results with Splinesoft =

Comments

Download in other formats: