Calculations with VisIt

Focusing on mass loss rate at the moment. Testing on Run5 (Aniso) so I can compare to Jonathan's paper. Created all the expressions I believe are required (see attachment:expressions.xml​), with the mass flux defined as density times radial velocity² (properly scaled). When I perform a weighted query on the mass flux over a spherical isosurface (defined in terms of distance from origin), however, I get a value about 3 magnitudes too large (1.14*10¹³). Also tried calculating flux with the mesh flux operator, but velocity times density doesn't seem to be the appropriate quantity. It does make nice pictures, though:

First is at r = 0.1, second is at r = 0.05 (Plot is pseudocolor of operator/Flux/Mesh, SphericalSlice centered at origin). Note that values are negative because the normals of the spherical slice point inwards.

Fixed scaling of expressions. Now summing over an isosurface gives the expected values.

Parameter Space for Planet Sims

Successfully recreated Run 5 (Aniso) from Jonathan's simulations, after meeting and fixing some of the code. Most runs are taking 5-8 hours with 1 level AMR (did check, optimization flag is set to level 3 in makefile). Used momenta to compare relative strength of outflows, made density movies.

In comparison to Aniso, changing lambda performs as expected - with lambda = 5, there is a marked increase in the strength of the wind, and with lambda = 15 there is very little to no wind (size of planet changes as well).

Changing the mass of the planet wasn't quite as clear. Overall, it doesn't appear to have affected the wind very much. Again, the planet size is related to the mass (since density is also specified). It is hard to see a wind on the smallest (.25 M_J) planet, but momentum plots appear to support its existence.

Changing the temperature of the planet also affects the radius (through lambda). A 5x hotter planet appears to have a stronger outflow at the end, although it takes a significant amount of time to gain significant size (perhaps should run this simulation for longer). A planet with .5 T doesn't appear to be very different, qualitatively, from the original run.

Finally (for the moment), changed ambient temperature to lower values. Would also like to explore higher values to see if they have much effect. The 50K, 25K, and 3K runs all have similar maximum momenta, approx. 2x the Aniso run, which seems to suggest a stronger outflow, but the density plots look nearly identical.

What parameters are best to plot?

Experimenting with AstroBEAR

Ran OutflowWind simulations with varying parameters. Thickness and wind velocity don't appear to have any effect - perhaps an effect of rescaling? Velocity speeds up the simulation as a whole. The rest are as follows:

Reference (default parameters, 100 frames, 2 units computational time (conversion?)):

Radius increases size of wind, but otherwise has little effect, at least at small values (expected).

Radius of 2:

Increasing density appears to decrease relative density at the center of the front of the outflow (slightly counter-intuitive, but makes sense).

Density of 5: Density of 10: Density of 20:

Increasing temp. creates a more diffuse cloud at the end (predictably). Temp = 0 appears to cause buggy behavior.

Temp = 0: Temp = 1: Temp = 5:

Also tried a run at 1000 frames, but definitely overkill for time.