Changes between Version 3 and Version 4 of u/erica/TestingOpSplitting


Ignore:
Timestamp:
09/15/13 12:16:51 (11 years ago)
Author:
Erica Kaminski
Comment:

Legend:

Unmodified
Added
Removed
Modified
  • u/erica/TestingOpSplitting

    v3 v4  
    11= Initialization =
    22
    3 To test the operator splitting, I initialized the grid with the same conditions as the Jeans test I did with astrobear (see that wiki page under the self-gravity section of my user page). Namely, I
    4 
    5 computed the Jeans length (in cm) for a gas of density rho = 3.34d-21 g/cc, temperature T = 100 K, and gamma = 1.0001. I then seeded the grid (that was 50x Jeans length long) with the prescribed
    6 
    7 density, velocity, and pressure perturbations as derived before. The density form was straightforward,
     3To test the operator splitting, I initialized the grid with the same conditions as the Jeans test I did with astrobear (see that wiki page under the self-gravity section of my user page). Namely, I computed the Jeans length (in cm) for a gas of density rho = 3.34d-21 g/cc, temperature T = 100 K, and gamma = 1.0001. I then seeded the grid (that was 50x Jeans length long) with the prescribed density, velocity, and pressure perturbations as derived before. The density form was straightforward,
    84
    95
     
    2420[[latex($\Gamma = 4 \pi G \rho_0[1-(\lambda_J/ \lambda)^2]$)]]
    2521
    26 In order to keep the velocity perturbation dimensionless, the factor out in front of the Sin(kx) needed to be scaled to computational units. This involved a factor of lscale for k, rhoscale for
    27 
    28 rho_0 in the growthrate, and a term for scaled gravity in the growthrate.
     22In order to keep the velocity perturbation dimensionless, the factor out in front of the Sin(kx) needed to be scaled to computational units. This involved a factor of lscale for k, rhoscale for rho_0 in the growthrate, and a term for scaled gravity in the growthrate.
    2923
    3024= Scaling gravitational constant =