Mathematical Model UNTRIM
From BAWiki
Contents
Short Description
The numerical method UNTRIM was developed by Prof. Vincenzo Casulli (Trento University, Italy). UNTRIM is a semi-implicit finite difference (-volume) model based on the three-dimensional shallow water equations as well as on the three-dimensional transport equation for salt, heat, dissolved matter and suspended sediments.
UNTRIM is able to work on unstructured orthogonal grids (UOG). The modelling domain is covered by a grid consisting of a set of non-overlapping convex polygons, usually either triangles or quadrilaterals. The grid is said to be an unstructured orthogonal grid if within each polygon a point (hereafter called a center) can be identified in such a way that the segment joining the center of two adjacent polygons and the side shared by the two polygons, have a non-empty intersection and are orthogonal to each other.
Physical Processes
- reynolds-averaged Navier-Stokes equations (RANS)
- local acceleration (inertia)
- advective acceleration
- Coriolis acceleration
- barotropic pressure gradient
- baroclinic pressure gradient
- hydrostatic or non-hydrostatic pressure
- horizontal turbulent viscosity
- vertical turbulent viscosity influenced by density stratification
- bottom friction
- wind friction
- sources and sinks
- horizontal acceleration due to wave effects (by means of radiation stress)
- transport of tracers
- local rate of change of concentration
- advective rate of change of concentration
- optional flux limiter : Minmod, van Leer or Superbee
- horizontal turbulent diffusivity
- vertical turbulent diffusivity influenced by density stratification
- settling of particles, deposition and erosion (for suspended sediments)
- heat-transfer to/from the atmosphere and to/from the bottom
- sources and sinks
- sinks with immediate return inflow at a different location, with optional modification of inflow-temperature as well as -salinity
Computational Results
- water level elevation at the free surface
- current velocity
- tracer concentration (e.g. salinity, temperature, suspended sediments)
- hydrodynamic pressure
Note: when UNTRIM is used in two-dimensional (depth-integrated) mode, results correspond to the depth-averaged values for the above-mentioned quantities.
Publications
- Casulli, V. and Zanolli, P. (1998), A Three-Dimensional Semi-Implicit Algorithm for Environmental Flows on Unstructured Grids, Proc. of Conf. on Num. Methods for Fluid Dynamics, University of Oxford.
- Casulli, V. (1999), A Semi-Implicit Finite Difference Method for Non-Hydrostatic, Free-Surface Flows, International Journal for Numerical Methods in Fluids, 30: 425 - 440.
- Casulli, V. and R.A. Walters (2000), An unstructured grid, three-dimensional model based on the shallow water equations. International Journal for Numerical Methods in Fluids Volume 32, Issue 3, pages 331 - 348.
- Casulli, V. and Zanolli, P. (2002), Semi-Implicit Numerical Modelling of Non-Hydrostatic Free-Surface Flows for Environmental Problems, Mathematical and Computer Modelling, 36: 1131 - 1149.
- Casulli, V. and Zanolli, P. (2004), High Resolution Methods for Multidimensional Advection-Diffusion Problems in Free-Surface Hydrodynamics, Ocean Modelling.
Other Users
Different versions of the computational core UnTRIM developed by Prof. V. Casulli are used at the following organizations:
- Delta Modeling Associates, Inc., San Francisco, California, USA
- Resource Management Associates, Fairfield, California, USA
- Stanford University, Stanford, California, USA
- University of California, Davis, California, USA
- Virginia Institute of Marine Science, Gloucester Point, Maryland, USA
Similar Software
The following models were inspired by some basics UnTRIM concepts developed and published by Prof. V. Casulli:
- D-Flow Flexible Mesh
- ELCIRC, SELFE
- SUNTANS
Validation Document
In addition to a quite general introduction to UNTRIM the validation document contains also more detailed informations related to the following topics:
- physical system,
- model functionality,
- conceptual model,
- algorithmic implementation,
- software implementation,
- validation studies, and
- literature.
A PDF-version of the validation document is freely available for download:
User Interface Description
This document contains a detailed description of all interface functions and routines available to the user. The following topics are dealt with in this document:
- set data (set-interfaces),
- get data (get-interfaces),
- check grid consistency as well as accuracy of iteratively computed results (check-routines),
- external routines called by the computational core (user-interface-routines) which are required to,
- define paths and names of the standard input data files, to
- define (set) the inital state (initial data), to
- set the forcing terms (e.g. along open boundaries) for each time step, and to
- retrieve the computational results.
- tables with short descriptions of all get- and set-interfaces available, and
- example standard input data files.
A PDF-version of the user interface description document is freely available for download:
MPI-Parallelisation
The core of the mathematical model UnTRIM has been parallelized using MPI by Jacek Jankowski. A detailed description is available in the technical report (approx. 1.1 MB) MPI Version Manual. This version is actually applied internally mainly in the department of Hydraulic Engineering in Inland Areas.
BAW-Specific Informations
Grid Generation
An unstructured orthogonal grid for UNTRIM can be prepared using JANET grid generator software, made by SmileConsult. For further informations related to the integration of JANET into BAW's programming environment please visit JANET program description.
Simulation
The mathematical model UNTRIM is fully integrated into BAW's programming environment. More detailed information concerning it's integration can be found visiting UNTRIM program description.
Graphical Presentation of Computed Results
To display UnTRIM results currently several methods are used at BAW. The more important ones are,
- HVIEW2D, for data available throughout the computational domain,
- VVIEW2D and/or LQ2PRO, for data at longitudinal- and/or cross-sections, as well as
- GVIEW2D, for data at specific locations.
Analyses of Computational Results
A great variety of methods for analyses of computational results is available which enables the user to respond to many different questions.
Coupling to Independent Sub-Models
At BAW UNTRIM can be used together with the following independent sub-models:
- spectral wave model k-modell (see unk.dat);
- sedimentological model SediMorph (see sedimorph.dat).
The above mentioned sub-models can be used in direct coupling with the computational core of UNTRIM.
back to Mathematical Models for Coastal Areas and Estuaries