GfsPhysicalParams

From Gerris

Revision as of 04:04, 21 September 2007 GeordieMcBain (Talk | contribs) (markup typo) ← Previous diff		Current revision Popinet (Talk | contribs)
Line 1:		Line 1:
-	Use GfsPhysicalParams to set the gravitational field strength <code>g</code> (a number) and the specific volume (the reciprocal of density) <code>alpha</code> (a scalar field).	+	GfsPhysicalParams controls the following physical parameters
		+
		+	{
		+	L = NUMBER
		+	g = NUMBER
		+	alpha = [ [[GfsFunction]] ]
		+	}
-	For example, from the [http://gfs.sourceforge.net/tests/tests/geo.html Geostrophic adjustment] test case	+	with
-	PhysicalParams { g = 9.4534734306584e-4 }	+	;<code>L</code>: the size of a [[GfsBox]] (default is unity),
		+	;<code>g</code>: for the [[GfsOcean]] and [[GfsRiver]] solvers only, the strength of the gravitational field (default is unity),
		+	;<code>alpha</code>: for the [[GfsSimulation]] solver, the specific volume (i.e the inverse of density) (default is unity).
-	and from the [http://gfs.sourceforge.net/examples/examples/rt.html Rayleigh&ndash;Taylor] example	+	<examples/>
-		+
-	VariableTracerVOF {} T	+
-	InitFraction {} T (0.05*cos (2.*M_PI*x) + y) { ty = 0.5 }	+
-	PhysicalParams { alpha = 1./(T*1.225 + (1. - T)*0.1694) }	+
-		+
-	(here the fluid with characteristic function <code>T</code> has density 1.225 and the other fluid has density 0.1694).	+

---

Current revision

GfsPhysicalParams controls the following physical parameters

{
  L = NUMBER
  g = NUMBER
  alpha = [ GfsFunction ]
}

with

L

the size of a GfsBox (default is unity),

g

for the GfsOcean and GfsRiver solvers only, the strength of the gravitational field (default is unity),

alpha

for the GfsSimulation solver, the specific volume (i.e the inverse of density) (default is unity).

Examples

• Rayleigh-Taylor instability

  PhysicalParams { alpha = 1./(T*1.225 + (1. - T)*0.1694) }

• Collapse of a column of grains

    PhysicalParams { L = LDOMAIN }

    PhysicalParams { alpha = 1./RHO(T) }

• Savart--Plateau--Rayleigh instability of a water column

    PhysicalParams { alpha = 1./RHO(T1) }

• Atomisation of a pulsed liquid jet

    PhysicalParams { alpha = 1./rho(T) }

• Air-water flow around a Series 60 cargo ship

    PhysicalParams { alpha = 1./VAR(T1,RATIO,1.) }

• Lunar tides in Cook Strait, New Zealand

    PhysicalParams { L = 500e3 }

    PhysicalParams { g = 9.81 }

• Dam break on complex topography

    PhysicalParams { L = 8. }

    PhysicalParams { g = 1. }

• Small amplitude solitary wave interacting with a parabolic hump

    PhysicalParams { g = 1 }

• Shock reflection by a circular cylinder

    PhysicalParams { L = 5 g = 9.81 }

• Tsunami runup onto a complex three-dimensional beach

    PhysicalParams { L = 3.402 g = 9.81 }

• The 2004 Indian Ocean tsunami

    PhysicalParams { 
	# length of the domain (m)
	L = LENGTH 
        # gravity is 9.81 m/s^2
	g = 9.81
	# from now on, units have been chosen to be metres and seconds
    }

• "Garden sprinkler effect" in wave model

    PhysicalParams { L = 5000 }

• Cyclone-generated wave field

    PhysicalParams { L = 3328 }

• Rotation of a straight interface

    PhysicalParams { L = 2. }

• Flow created by a cylindrical volume source

    PhysicalParams { L = 2 }

• Potential flow around a sphere

    PhysicalParams { L = 50 }

• Viscous flow past a sphere

    PhysicalParams { L = 50 }

• Boundary layer on a rotating disk

    PhysicalParams { L = 3 }

• Poiseuille flow with multilayer Saint-Venant

    PhysicalParams { L = 1. }

• Momentum conservation for large density ratios

    PhysicalParams { alpha = 1./rho(T1) }

• Coriolis formulation in 3-D

  PhysicalParams { L = 60 }

• Multi-layer Saint-Venant solver

    PhysicalParams { L = RATIO g = 100./RE }

• Wind-driven stratified lake

    PhysicalParams {
	L = RATIO g = 1./(ALPHA*RE)
	alpha = 1./(1. + DRHO)
    }

• Air-Water capillary wave

  PhysicalParams { alpha = 1./RHO(T1) }

• Fluids of different densities

  PhysicalParams { alpha = 1./(T + 0.1*(1. - T)) }

• Pure gravity wave

  PhysicalParams { alpha = 1./(T + 0.1*(1. - T)) }

• Shape oscillation of an inviscid droplet

    PhysicalParams { alpha = 1./RHO(T1) }

• Scalings for Plateau--Rayleigh pinchoff

   PhysicalParams { alpha = 1./(T + 1e-2*(1. - T)) }

• Sessile drop

    PhysicalParams { alpha = 1./(T + 0.01*(1. - T)) }

• Geostrophic adjustment

  PhysicalParams { L = L0 g = G }

• Geostrophic adjustment on a beta-plane

  PhysicalParams { g = 9.4534734306584e-4 }

• Geostrophic adjustment with Saint-Venant

  PhysicalParams { L = L0 g = G }

• Non-linear geostrophic adjustment

  PhysicalParams { L = 1 }

• Coastally-trapped waves

  PhysicalParams { g = 5.87060327757e-3 }

• Coastally-trapped waves with adaptive refinement

  PhysicalParams { g = 5.87060327757e-3 }

• Gravity waves in a realistic ocean basin

    PhysicalParams { g = 19.62 }

• Oscillations in a parabolic container

    PhysicalParams { L = 10000 }

    PhysicalParams { g = G }

• Parabolic container with embedded solid

    PhysicalParams { L = 10000 }

    PhysicalParams { g = G }

• Transcritical flow over a bump

    PhysicalParams { L = 25 g = 9.81 }

• Transcritical flow with multiple layers

    PhysicalParams { L = LENGTH g = G }

• Tsunami runup onto a plane beach

    PhysicalParams { L = 60000 }

    PhysicalParams { g = 9.81 }

• Lake-at-rest balance in an inclined domain with cut cells

    PhysicalParams { L = 10 g = 9.81 }

• Lake-at-rest balance in an inclined domain with bipolar metric

    PhysicalParams { g = 9.81 }

• Terrain reconstruction

    PhysicalParams { L = 8 }

• Circular dam break on a sphere

    PhysicalParams { L = LENGTH }

• Circular dam break on a rotating sphere

    PhysicalParams { L = LENGTH }

• Circular dam break on a ``cubed sphere''

    PhysicalParams { L = 2.*M_PI/4. }

• Advection of a cosine bell around the sphere

  PhysicalParams { L = 2.*M_PI/4. }

• Poisson problem with a pure spherical harmonics solution

  PhysicalParams { L = 2.*M_PI/4. }

• Spherical harmonics with longitude-latitude coordinates

  PhysicalParams { L = M_PI }

• Poisson equation on a sphere with Gaussian forcing

  PhysicalParams { L = 2.*M_PI/4. }

• Gaussian forcing using longitude-latitude coordinates

  PhysicalParams { L = M_PI }

• Creeping Couette flow between eccentric cylinders

  PhysicalParams { L = 2.5 }

• Flow between eccentric cylinders on a stretched grid

  PhysicalParams { L = 2.5 }

• Rossby--Haurwitz wave

    PhysicalParams { L = 2.*M_PI*AR/4. }

• Rossby--Haurwitz wave with a free surface

    PhysicalParams { 
	L = 2.*M_PI*AR/4. 
        # g*H0
	g = G*8e3
    }

• Rossby--Haurwitz wave with Saint-Venant

    PhysicalParams { L = 2.*M_PI*AR/4. g = G }

• Charge relaxation in a planar cross-section

    PhysicalParams { L = 2 }

• Equilibrium of a droplet suspended in an electric field

    PhysicalParams { L = 2 }

• Groundwater flow with piecewise constant permeability

    PhysicalParams { alpha = k }