Commit faaf4fc8 authored by Matthieu Schaller's avatar Matthieu Schaller
Browse files

Updated the Sedov blast test case. The need for a better time step size estimator is now reduced.


Former-commit-id: d132808dfb49add997be399b9790a4444d8f600b
parent c969365b
...@@ -26,11 +26,11 @@ from numpy import * ...@@ -26,11 +26,11 @@ from numpy import *
# Parameters # Parameters
periodic= 1 # 1 For periodic box periodic= 1 # 1 For periodic box
boxSize = 1. boxSize = 10.
L = 101 # Number of particles along one axis L = 101 # Number of particles along one axis
rho = 1. # Density rho = 1. # Density
P = 1.e-5 # Pressure P = 1.e-5 # Pressure
E0= 1.e5 # Energy of the explosion E0= 1.e2 # Energy of the explosion
pert = 0.1 pert = 0.1
gamma = 5./3. # Gas adiabatic index gamma = 5./3. # Gas adiabatic index
fileName = "sedov.hdf5" fileName = "sedov.hdf5"
...@@ -71,7 +71,7 @@ for i in range(L): ...@@ -71,7 +71,7 @@ for i in range(L):
u[index] = internalEnergy u[index] = internalEnergy
ids[index] = index ids[index] = index
if sqrt((x - boxSize/2.)**2 + (y - boxSize/2.)**2 + (z - boxSize/2.)**2) < 2.01 * boxSize/L: if sqrt((x - boxSize/2.)**2 + (y - boxSize/2.)**2 + (z - boxSize/2.)**2) < 2.01 * boxSize/L:
u[index] = u[index] + E0 / 33. u[index] = u[index] + E0 / (33. * mass)
coords[index,0] = x + random.random() * pert * boxSize/(2.*L) coords[index,0] = x + random.random() * pert * boxSize/(2.*L)
coords[index,1] = y + random.random() * pert * boxSize/(2.*L) coords[index,1] = y + random.random() * pert * boxSize/(2.*L)
coords[index,2] = z + random.random() * pert * boxSize/(2.*L) coords[index,2] = z + random.random() * pert * boxSize/(2.*L)
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment