Meeting update

Here is my page on operator splitting, with derived update formulas I am using-

​https://astrobear.pas.rochester.edu/trac/astrobear/wiki/u/erica/PoissonPlusHydro

Here is my page on setting up the Jeans instability for my code, including how I scaled the gravitational constant,

​https://astrobear.pas.rochester.edu/trac/astrobear/wiki/u/erica/TestingOpSplitting

The results section is not filled in, as I am still working out some kinks in the code. There seems to be a bug somewhere that is causing nans…

Along with this, will be running some new BE sphere simulations, and preparing for the MHD colliding flows simulations, as well as editing the paper with Phil's suggestions.

I am not sure why this code is not giving me an oscillation in p init (W3 array),

DO i = 1, cells

x = REAL((REAL(i)-0.5)*dx)

W1(i) = rho0 + delta*Cos(k*x)

W2(i) = -(delta*growthrate*lambda)/(rho0*2*Pi)*Sin(k*x)

W3(i) = p0 + delta*((W1(i)*rhoscale*Kb*100)/mh)/pscale !ideal gas law in cgs for pressure, over pscale

WRITE(85, *) X, W1(i), W2(i), W3(i) !Compute conserved quantities

CS(1,i) = W1(i)

CS(2,i) = W2(i)*W1(i)

CS(3,i) = 0.5*w1(i)*w2(i)*w2(i) + w3(i)/(gamma - 1.0)

END DO

Comments

No comments.