TITAN2D

Last updated
TITAN2D
Developer(s) Geophysical Mass Flow Group
Stable release
2.0.0 / July 21, 2007;15 years ago (2007-07-21)
Operating system Unix-like
Type Geoflow Simulator
License NCSA Open Source License
Website www.gmfg.buffalo.edu github.com/TITAN2D/titan2d

TITAN2D is a geoflow simulation software application, intended for geological researchers. It is distributed as free software.

Contents

Overview

TITAN2D is a free software application developed by the Geophysical Mass Flow Group at the State University of New York (SUNY) at Buffalo. TITAN2D was developed for the purpose of simulating granular flows (primarily geological mass flows such as debris avalanches and landslides) over digital elevation models (DEM)s of natural terrain. The code is designed to help scientists and civil protection authorities assess the risk of, and mitigate, hazards due to dry debris flows and avalanches. TITAN2D combines numerical simulations of a flow with digital elevation data of natural terrain supported through a Geographical Information System (GIS) interface such as GRASS.

TITAN2D is capable of multiprocessor runs. A Message Passing Interface (MPI) Application Programming Interface (API) allows for parallel computing on multiple processors, which effectively increases computational power, decreases computing time, and allows for the use of large data sets.

Adaptive gridding allows for the concentration of computing power on regions of special interest. Mesh refinement captures the complex flow features that occur at the leading edge of a flow, as well as locations where rapid changes in topography induce large mass and momentum fluxes. Mesh unrefinement is applied where solution values are relatively constant or small to further improve computational efficiency.

TITAN2D requires an initial volume and shape estimate for the starting material, a basal friction angle, and an internal friction angle for the simulated granular flow. The direct outputs of the program are dynamic representations of a flow's depth and momentum. Secondary or derived outputs include flow velocity, and such field-observable quantities as run-up height, deposit thickness, and inundation area.

Mathematical Model

The TITAN2D program is based upon a depth-averaged model for an incompressible Coulomb continuum, a “shallow-water” granular flow. The conservation equations for mass and momentum are solved with a Coulomb-type friction term for the interactions between the grains of the media and between the granular material and the basal surface. The resulting hyperbolic system of equations is solved using a parallel, adaptive mesh, Godunov scheme. The basic form of the depth-averaged governing equations appear as follows.

The depth-averaged conservation of mass is:

The depth-averaged x,y momentum balances are:


See also

Related Research Articles

<span class="mw-page-title-main">Navier–Stokes equations</span> Equations describing the motion of viscous fluid substances

The Navier–Stokes equations are partial differential equations which describe the motion of viscous fluid substances, named after French engineer and physicist Claude-Louis Navier and Anglo-Irish physicist and mathematician George Gabriel Stokes. They were developed over several decades of progressively building the theories, from 1822 (Navier) to 1842-1850 (Stokes).

In vector calculus, the divergence theorem, also known as Gauss's theorem or Ostrogradsky's theorem, is a theorem which relates the flux of a vector field through a closed surface to the divergence of the field in the volume enclosed.

In fluid dynamics, turbulence or turbulent flow is fluid motion characterized by chaotic changes in pressure and flow velocity. It is in contrast to a laminar flow, which occurs when a fluid flows in parallel layers, with no disruption between those layers.

In fluid dynamics, Stokes' law is an empirical law for the frictional force – also called drag force – exerted on spherical objects with very small Reynolds numbers in a viscous fluid. It was derived by George Gabriel Stokes in 1851 by solving the Stokes flow limit for small Reynolds numbers of the Navier–Stokes equations.

<span class="mw-page-title-main">Boundary layer</span> Layer of fluid in the immediate vicinity of a bounding surface

In physics and fluid mechanics, a boundary layer is the thin layer of fluid in the immediate vicinity of a bounding surface formed by the fluid flowing along the surface. The fluid's interaction with the wall induces a no-slip boundary condition. The flow velocity then monotonically increases above the surface until it returns to the bulk flow velocity. The thin layer consisting of fluid whose velocity has not yet returned to the bulk flow velocity is called the velocity boundary layer.

A continuity equation or transport equation is an equation that describes the transport of some quantity. It is particularly simple and powerful when applied to a conserved quantity, but it can be generalized to apply to any extensive quantity. Since mass, energy, momentum, electric charge and other natural quantities are conserved under their respective appropriate conditions, a variety of physical phenomena may be described using continuity equations.

<span class="mw-page-title-main">Stellar dynamics</span>

Stellar dynamics is the branch of astrophysics which describes in a statistical way the collective motions of stars subject to their mutual gravity. The essential difference from celestial mechanics is that the number of body

<span class="mw-page-title-main">Open-channel flow</span> Type of liquid flow within a conduit

In fluid mechanics and hydraulics, open-channel flow is a type of liquid flow within a conduit with a free surface, known as a channel. The other type of flow within a conduit is pipe flow. These two types of flow are similar in many ways but differ in one important respect: open-channel flow has a free surface, whereas pipe flow does not.

In fluid dynamics, turbulence kinetic energy (TKE) is the mean kinetic energy per unit mass associated with eddies in turbulent flow. Physically, the turbulence kinetic energy is characterised by measured root-mean-square (RMS) velocity fluctuations. In the Reynolds-averaged Navier Stokes equations, the turbulence kinetic energy can be calculated based on the closure method, i.e. a turbulence model.

<span class="mw-page-title-main">Shallow water equations</span> Set of partial differential equations that describe the flow below a pressure surface in a fluid

The shallow-water equations (SWE) are a set of hyperbolic partial differential equations that describe the flow below a pressure surface in a fluid. The shallow-water equations in unidirectional form are also called Saint-Venant equations, after Adhémar Jean Claude Barré de Saint-Venant.

The turbulent Prandtl number (Prt) is a non-dimensional term defined as the ratio between the momentum eddy diffusivity and the heat transfer eddy diffusivity. It is useful for solving the heat transfer problem of turbulent boundary layer flows. The simplest model for Prt is the Reynolds analogy, which yields a turbulent Prandtl number of 1. From experimental data, Prt has an average value of 0.85, but ranges from 0.7 to 0.9 depending on the Prandtl number of the fluid in question.

Pressure-correction method is a class of methods used in computational fluid dynamics for numerically solving the Navier-Stokes equations normally for incompressible flows.

In physical oceanography and fluid dynamics, the wind stress is the shear stress exerted by the wind on the surface of large bodies of water – such as oceans, seas, estuaries and lakes. Stress is the quantity that describes the magnitude of a force that is causing a deformation of an object. Therefore, stress is defined as the force per unit area and its SI unit is the Pascal. When the deforming force acts parallel to the object's surface, this force is called a shear force and the stress it causes is called a shear stress. When wind is blowing over a water surface, the wind applies a wind force on the water surface. The wind stress is the component of this wind force that is parallel to the surface per unit area. Also, the wind stress can be described as the flux of horizontal momentum applied by the wind on the water surface. The wind stress causes a deformation of the water body whereby wind waves are generated. Also, the wind stress drives ocean currents and is therefore an important driver of the large-scale ocean circulation. The wind stress is affected by the wind speed, the shape of the wind waves and the atmospheric stratification. It is one of the components of the air–sea interaction, with others being the atmospheric pressure on the water surface, as well as the exchange of energy and mass between the water and the atmosphere.

In nonideal fluid dynamics, the Hagen–Poiseuille equation, also known as the Hagen–Poiseuille law, Poiseuille law or Poiseuille equation, is a physical law that gives the pressure drop in an incompressible and Newtonian fluid in laminar flow flowing through a long cylindrical pipe of constant cross section. It can be successfully applied to air flow in lung alveoli, or the flow through a drinking straw or through a hypodermic needle. It was experimentally derived independently by Jean Léonard Marie Poiseuille in 1838 and Gotthilf Heinrich Ludwig Hagen, and published by Poiseuille in 1840–41 and 1846. The theoretical justification of the Poiseuille law was given by George Stokes in 1845.

<span class="mw-page-title-main">Boundary current</span> Ocean current with dynamics determined by the presence of a coastline

Boundary currents are ocean currents with dynamics determined by the presence of a coastline, and fall into two distinct categories: western boundary currents and eastern boundary currents.

In fluid dynamics, Airy wave theory gives a linearised description of the propagation of gravity waves on the surface of a homogeneous fluid layer. The theory assumes that the fluid layer has a uniform mean depth, and that the fluid flow is inviscid, incompressible and irrotational. This theory was first published, in correct form, by George Biddell Airy in the 19th century.

In engineering, physics, and chemistry, the study of transport phenomena concerns the exchange of mass, energy, charge, momentum and angular momentum between observed and studied systems. While it draws from fields as diverse as continuum mechanics and thermodynamics, it places a heavy emphasis on the commonalities between the topics covered. Mass, momentum, and heat transport all share a very similar mathematical framework, and the parallels between them are exploited in the study of transport phenomena to draw deep mathematical connections that often provide very useful tools in the analysis of one field that are directly derived from the others.

<span class="mw-page-title-main">Geophysical fluid dynamics</span> Dynamics of naturally occurring flows

Geophysical fluid dynamics, in its broadest meaning, refers to the fluid dynamics of naturally occurring flows, such as lava flows, oceans, and planetary atmospheres, on Earth and other planets.

<span class="mw-page-title-main">Radiation stress</span> Term in physical oceanography

In fluid dynamics, the radiation stress is the depth-integrated – and thereafter phase-averaged – excess momentum flux caused by the presence of the surface gravity waves, which is exerted on the mean flow. The radiation stresses behave as a second-order tensor.

Nonlinear tides are generated by hydrodynamic distortions of tides. A tidal wave is said to be nonlinear when its shape deviates from a pure sinusoidal wave. In mathematical terms, the wave owes its nonlinearity due to the nonlinear advection and frictional terms in the governing equations. These become more important in shallow-water regions such as in estuaries. Nonlinear tides are studied in the fields of coastal morphodynamics, coastal engineering and physical oceanography. The nonlinearity of tides has important implications for the transport of sediment.

References