wiki:u/erica/scratch5

Version 3 (modified by Erica Kaminski, 9 years ago) ( diff )

Set the following scales in physics.data:

lscale 3.08 e+18 cm
rscale 2.02 e-16 g/cm3
tempscale 1 K

From that initialized an ambient object with the following density and pressure,

ambient%rho=2.03d-16/rscale
ambient%pressure=(rhoOut/rscale)*(tempOut/tempscale)

Have my grid running from, gxbounds= -.0068 -.0068 0 .0068 .0068

I ran the sim to frame 0 to get the timescale from scales.data.

I used this to get the tsim=tdiff in computational units,

tsim=2.074e+6s/!338402526430554s = 6.13e-9

Next, I calculated the energy injection rate using the luminosity of a FHSC:

L=3.9e+29 s
E(radiated in a diffusion time)=L*tdiff=8e+35 erg

I convert this into an energy density, in computational units, by dividing by pscale,

e(cu)=8e+35erg/Pscale=4.74e+43

From this I can get an energy injection rate (i.e. the amount of energy density injected into the kernel each timestep),

de/dt=e(cu)/tsim(cu)=4.74e+43/6.13e-9=7.73246e+51.

Running this high of an de/dt gave problems, so I decreased this to:

de/dt=e(cu)/tsim(cu)=4.74e+43/6.13e-9=7.73246e+21.

Next came scaling the opacities.

The code takes specific opacity (this is in cm2/g), and then asks whether it should be constant, or multiplied by density (or temperature) to some power. I set it to be a power law scaling.

tpow=0
dpow=1

so that,

opacity=kappa*rho

I scaled the planck and rosseland opacities into computational units by:

kappa_r=0.23*mscale/(lscale)2=144.07
kappa_p=.4*mscale/(lscale)2=250.56

Radiation was set to 'diffusion and energy exchange only', rad limiter was set to 'diffusion limit'.

This all seems right to me, so then here are the results.

Run 1

de/dt=7.73246e+51

The code failed to run past frame zero.

Attachments (51)

Note: See TracWiki for help on using the wiki.