Changes between Version 8 and Version 9 of u/erica/JeansTest


Ignore:
Timestamp:
06/24/13 14:39:54 (11 years ago)
Author:
Erica Kaminski
Comment:

—

Legend:

Unmodified
Added
Removed
Modified

  • u/erica/JeansTest

    v8 v9  
    7575= The code =
    7676
     77The following is for a 1D isothermal sim of the Jeans Instability.
     78
     79In !ProblemModuleInit, I read in ambient density and temperature. These are used to uniformly set a corresponding Ambient Object.
     80
     81[[latex($READ ~ rho$)]]
     82
     83
     84[[latex($Info \% q(i, :, :, 1) = rho$)]]
     85
     86[[latex($Ambient\%density = info \% q(i, :, :, 1) $)]]
     87
     88
     89In ProblemGridInit, I loop over the cells in x (indexed by i) adding the following density perturbation:
     90
     91[[latex($info \% q(i, :, :, 1) = info\%q(i, :, :, 1) + 0.1*info\%q(i, :, :, 1)*Cos(kx)*e^0 $)]]
     92
     93
     94I am not 100% sure, but I think if I wanted to place the analytical function in astrobear, I would then add to !ProblemBeforeStep a loop over cells in i that continues to perturb the density function in time t:
     95
     96[[latex($info \% q(i, :, :, 1) = info\%q(i, :, :, 1, t=0)*e^{\omega t} $)]]
     97
     98Thus we should see the characteristic growth rate of [[latex($\omega$)]].
     99
     100
    77101= Results =
     102
     103Analytic function -
     104
     105Since I am having difficulty saving an initial density function in astrobear for use with !ProblemBeforeStep, I just plotted the analytical function for the initial density perturbation and the perturbation after an e-folding time in Mathematica. The e-folding time in physical units is -  These