From 32633fd2013fd1795942ab935926ba281ecc4bca Mon Sep 17 00:00:00 2001 From: bangerth Date: Thu, 8 May 2008 22:15:43 +0000 Subject: [PATCH] Add input.prm and reference it. git-svn-id: https://svn.dealii.org/trunk@16063 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/examples/step-33/doc/results.dox | 2 +- deal.II/examples/step-33/input.prm | 86 ++++++++++++++++++++++++ 2 files changed, 87 insertions(+), 1 deletion(-) create mode 100644 deal.II/examples/step-33/input.prm diff --git a/deal.II/examples/step-33/doc/results.dox b/deal.II/examples/step-33/doc/results.dox index 586a91170e..6bdcb6c4c4 100644 --- a/deal.II/examples/step-33/doc/results.dox +++ b/deal.II/examples/step-33/doc/results.dox @@ -3,7 +3,7 @@ We run the problem with the mesh slide.inp (this file is in the same directory as the source code for this program) and the following input -deck: +deck (available as input.prm in the same directory): @verbatim # Listing of Parameters # --------------------- diff --git a/deal.II/examples/step-33/input.prm b/deal.II/examples/step-33/input.prm new file mode 100644 index 0000000000..f2f5e6edd2 --- /dev/null +++ b/deal.II/examples/step-33/input.prm @@ -0,0 +1,86 @@ +# Listing of Parameters +# --------------------- + +# The input grid +set mesh = slide.inp + +# Stabilization parameter +set diffusion power = 2.0 + +# Scaled value for gravity. Positive means gravity points down. +set gravity = 1.0 + +# -------------------------------------------------- +# Boundary conditions +# We may specify boundary conditions for up to MAX_BD boundaries. +# Your .inp file should have these boundaries designated. +subsection boundary_1 + set no penetration = true # reflective boundary condition +end + +subsection boundary_2 + # outflow boundary + # set w_2 = pressure + # set w_2 value = 1.5 - y +end + +subsection boundary_3 + set no penetration = true # reflective + # set w_3 = pressure + # set w_3 value = 1.0 +end + +subsection boundary_4 + set no penetration = true #reflective +end + +# -------------------------------------------------- +# Initial Conditions +# We set the initial conditions of the conservative variables. These lines +# are passed to the expression parsing function. You should use x,y,z for +# the coordinate variables. + +subsection initial condition + set w_0 value = 0 + set w_1 value = 0 + set w_2 value = 10*(x<-0.7)*(y> 0.3)*(y< 0.45) + (1-(x<-0.7)*(y> 0.3)*(y< 0.45))*1.0 + set w_3 value = (1.5-(1.0*1.0*y))/0.4 +end + +# -------------------------------------------------- +# Time stepping control +subsection time stepping + set final time = 10.0 # simulation end time + set time step = 0.02 # simulation time step +end + +subsection linear solver + set output = quiet + set method = gmres + set ilut fill = 1.5 + set ilut drop tolerance = 1e-6 + set ilut absolute tolerance = 1e-6 + set ilut relative tolerance = 1.0 +end + +# -------------------------------------------------- +# Output frequency. +# You may wish to set this > time step if you dont want output at every step +subsection output + set step = 0.01 +end + +# -------------------------------------------------- +# Refinement control +subsection refinement + set refinement = shock # none only other option + set shock value = 1.5 + set shock levels = 1 # how many levels of refinement to allow +end + +# -------------------------------------------------- +# Flux parameters +subsection flux + set stab = mesh + #set stab value = 1.0 +end -- 2.39.5