From faaf4fc8cc05a6edadb07acce27ad6ee0f9e3b16 Mon Sep 17 00:00:00 2001
From: Matthieu Schaller <matthieu.schaller@durham.ac.uk>
Date: Sun, 10 Feb 2013 19:53:54 +0000
Subject: [PATCH] Updated the Sedov blast test case. The need for a better time
step size estimator is now reduced.
Former-commit-id: d132808dfb49add997be399b9790a4444d8f600b
---
examples/SedovBlast/makeIC.py | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/examples/SedovBlast/makeIC.py b/examples/SedovBlast/makeIC.py
index 7765149aa5..6a142334ed 100644
--- a/examples/SedovBlast/makeIC.py
+++ b/examples/SedovBlast/makeIC.py
@@ -26,11 +26,11 @@ from numpy import *
# Parameters
periodic= 1 # 1 For periodic box
-boxSize = 1.
+boxSize = 10.
L = 101 # Number of particles along one axis
rho = 1. # Density
P = 1.e-5 # Pressure
-E0= 1.e5 # Energy of the explosion
+E0= 1.e2 # Energy of the explosion
pert = 0.1
gamma = 5./3. # Gas adiabatic index
fileName = "sedov.hdf5"
@@ -71,7 +71,7 @@ for i in range(L):
u[index] = internalEnergy
ids[index] = index
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,1] = y + random.random() * pert * boxSize/(2.*L)
coords[index,2] = z + random.random() * pert * boxSize/(2.*L)
--
GitLab