BUNDESANSTALT FUR¨ WASSERBAU Karlsruhe • Hamburg • Ilmenau
Technical Report
Mathematical Model UnTRIM
Validation Document
– Version June 2004 (1.0) – BUNDESANSTALT FUR¨ WASSERBAU Karlsruhe • Hamburg • Ilmenau
Technical Report
Mathematical Model UnTRIM
Validation Document
– Version June 2004 (1.0) –
Contributors: Vincenzo Casulli (Trento University, Italy) G¨unther Lang (BAW, Germany)
Version-Date: June 2004
Version-no.: 1.0
Trento ¡ Hamburg — August 2004
¡ ¡ The Federal Waterways Engineering and Research Institute ¡ Wedeler Landstr. 157 D-22559 Hamburg ([+49] (0)40 81908 - 0 The Federal Waterways Engineering and Research Institute (BAW) Mathematical Model UnTRIM Validation Document – Version June 2004 (1.0)
Summary
This document is the validation document for the mathematical model UnTRIM. The docu- ment is organized conforming to Guidelines for documenting the validity of computational mod- elling software [IAHR, 1994]. The subject of this document is the validation of a computational model. The term compu- tational model refers to software whose primary function is to model a certain class of phys- ical systems, and may include pre- and post-processing components and other necessary ancillary programmes. Validation applies primarily to the theoretical foundation and to the computational techniques that form the basis for the numerical and graphical results pro- duced by the software. In the context of this document, validation of the model is viewed as the formulation and substantiation of explicit claims about applicability and accuracy of the computational results. This preface explains the approach that has been adopted in organizing and presenting the information contained in this document.
Standard validation document
This document conforms to a standard system for validation documentation. This system, the Standard Validation Document, has been developed by the hydraulic research indus- try in order to address the need for useful and explicit information about the validity of computational models. Such information is summarized in a validation document, which accompanies the technical reference documentation associated with a computational model. In conforming to the Standard, this validation document meets the following require- ments:
1. It has a prescribed table of contents, based on a framework that allows separate quality issues to be clearly distinguished and described.
2. It includes a comprehensive list of the assumptions and approximations that were made during the design and implementation of the model.
3. It contains claims about the performance of the model, together with statements that point to the available substantiating evidence for these claims.
4. Claims about the model made in this document are substantiable and bounded: they can be tested, justified, or supported by means of physical or computational experi- ments, theoretical analysis, or case studies.
5. Claims are substantiated by evidence contained within this document, or by specific reference to accessible publications.
6. Results of validation studies included or referred to in this document are reproducible. Consequently the contents of this document are consistent with the current version of the software.
Page i The Federal Waterways Engineering and Research Institute (BAW) Mathematical Model UnTRIM Validation Document – Version June 2004 (1.0)
7. This document will be updated as the process of validating the model progresses.
Organization of this document
Chapter 1 contains a short overview of the computational model and introduces the main issues to be addressed by the validation process. The model overview includes information about the purpose of the model, about pre- and post-processing options and other software features, and about reference versions of the software. Validation priorities and approaches are briefly described, and a list of related documents is included. Chapter 2 summarizes the available information about the validity of the computational core of the model. In this chapter, claims are made about the range of applicability of the model and about the accuracy of computational results. Each claim is followed by a brief statement regarding its substantiation. This statement indicates the extent to which the claim has in fact been substantiated and points to the available evidence. Chapter 3 contains such evidence, in the form of brief descriptions of relevant validation studies. Each description includes information about the purpose and approach of the study, and a summary of main results and implications. A glossary and complete list of references are contained in this document too.
A word of caution
This document contains information about the quality of a complex modelling tool. Its pur- pose is to assist the user in assessing the reliability and accuracy of computational results, and to provide guidelines with respect to the applicability and judicious employment of this tool. This document does not, however, provide mathematical proof of the correctness of re- sults for a specific application. The reader is referred to the License Agreement for pertinent legal terms and conditions associated with the use of the software. The contents of this validation document attest to the fact that computational modelling of complex physical systems requires great care and inherently involves a number of uncer- tain factors. In order to obtain useful and accurate results for a particular application, the use of high-quality modelling tools is necessary but not sufficient. Ultimately, the quality of the computational results that can be achieved will depend upon the adequacy of available data as well as a suitable choice of model and modelling parameters.
Electronic standard validation document
This document is also available in electronic form in Portable Document Format (PDF). The electronic version may be read using the ACROBAT READER software which is available for many computer platforms.
Page ii The Federal Waterways Engineering and Research Institute (BAW) Mathematical Model UnTRIM Validation Document – Version June 2004 (1.0)
Acknowledgements
Ralph T. Cheng is responsible for the development of a version of the UnTRIM at the U. S. Geological Survey (USGS). Continuing discussions and exchanging of modeling ideas with regards to UnTRIM validations and applications with Ralph T. Cheng are acknowledged. During the first International UnTRIM Users Workshop in Verona, June 7 – 9, 2004, it was decided to update this document, to conform to the actual version of the software (June 2004). Positive and stimulating feedback from the users to continue working on this type of document is herewidth gratefully acknowledged too.
Page iii The Federal Waterways Engineering and Research Institute (BAW) Mathematical Model UnTRIM Validation Document – Version June 2004 (1.0)
Contents
1 Introduction 1 1.1 Model overview ...... 1 1.1.1 Purpose ...... 1 1.1.2 Properties of the computational model ...... 2 1.1.3 Unstructured orthogonal grid ...... 2 1.1.4 Pre- and post-processing and other software features ...... 2 1.1.5 Version information ...... 4 1.2 Validation priorities and approaches ...... 4 1.3 Related documents ...... 4
2 Model validity 5 2.1 Physical system ...... 5 2.2 Model functionality ...... 6 2.2.1 Applications ...... 6 2.2.2 Processes ...... 11 2.3 Conceptual model ...... 14 2.3.1 Governing equations ...... 14 2.3.2 Assumptions and approximations ...... 17 2.3.3 Claims and substantiations ...... 18 2.4 Algorithmic implementation ...... 21 2.4.1 Unstructured orthogonal grid ...... 21 2.4.2 Assumptions and approximations ...... 24 2.4.3 Claims and substantiations ...... 33 2.5 Software implementation ...... 35 2.5.1 Implementation techniques ...... 35 2.5.2 Claims and substantiations ...... 37
3 Validation studies 40 3.1 Analytical test cases ...... 41 3.1.1 Wave propagation: rectangular basin ...... 41 3.1.2 Free barotropic oscillations: rectangular basin ...... 42 3.1.3 Free barotropic oscillations: circular basin ...... 43 3.1.4 Free oscillations: non-hydrostatic pressure ...... 44 3.1.5 Tidal forcing: flat bottom ...... 46 3.1.6 Tidal forcing: varying bathymetry ...... 47 3.1.7 Wind driven flow: flat bottom ...... 48 3.1.8 Internal seiches ...... 49 3.1.9 Steady density induced flow ...... 50 3.1.10 Wetting and Drying ...... 51 3.1.11 Solitary wave: flat bottom ...... 52 3.1.12 Solitary wave: varying bottom ...... 53
Page iv The Federal Waterways Engineering and Research Institute (BAW) Mathematical Model UnTRIM Validation Document – Version June 2004 (1.0)
3.2 Laboratory test cases ...... 54 3.2.1 Wave flume ...... 54 3.3 Schematic test cases ...... 56 3.3.1 Lock exchange flow: hydrostatic pressure ...... 56 3.3.2 Lock exchange flow: non-hydrostatic pressure ...... 58 3.3.3 Wave pattern in a square basin ...... 59 3.3.4 Short waves in a harbour basin ...... 60 3.3.5 Advection in a curved channel ...... 63 3.4 Examples from real-world applications ...... 66 3.4.1 Hydrostatic and non-hydrostatic flow in Venice Lagoon ...... 66 3.4.2 Tidal flow and salt transport ...... 70 3.4.3 River flow ...... 71 3.4.4 Storm surge ...... 72 3.4.5 Suspended sediment transport ...... 73 3.4.6 Transport of cooling water from a power plant ...... 74
References 75
A Glossary 78
Page v The Federal Waterways Engineering and Research Institute (BAW) Mathematical Model UnTRIM Validation Document – Version June 2004 (1.0)
List of Figures
1 Unstructured orthogonal grid ...... 3 2 Definition of total water depth ...... 15 3 Orthogonal unstructured grid ...... 22 4 Special unstructured orthogonal grids ...... 22 5 Mixed grid with grid refinement ...... 23 6 Software concept ...... 38 7 Free oscillations: analytical, hydrostatic and non-hydrostatic solution . . . . . 45 8 Wave flume: flume geometry ...... 54 9 Wave flume: results at 13.5 m from open boundary ...... 55 10 Wave flume: results at 15.7 m from open boundary ...... 55 11 Wave flume: results at 19 m from open boundary ...... 55 12 Lock exchange flow: hydrostatic computation ...... 57 13 Lock exchange flow: non-hydrostatic computation ...... 58 14 Harbour basin: grid ...... 61 15 Harbour basin: incoming waves ...... 61 16 Harbour basin: diffraction of waves ...... 62 17 Harbour basin: diffraction and reflection of waves ...... 62 18 Advection in U-channel: square grid ...... 64 19 Advection in U-channel: flow aligned grid ...... 64 20 Advection in U-channel: upwind on a square grid ...... 65 21 Advection in U-channel: flux limiter on a square grid ...... 65 22 Advection in U-channel: flux limiter on a flow aligned grid ...... 65 23 Venice Lagoon: unstructured orthogonal grid ...... 68 24 Venice Lagoon: vertical velocity for hydrostatic/non-hydrostatic pressure . . 69
Page vi The Federal Waterways Engineering and Research Institute (BAW) Mathematical Model UnTRIM Validation Document – Version June 2004 (1.0)
List of Tables
1 Venice Lagoon: model performance ...... 66
Page vii The Federal Waterways Engineering and Research Institute (BAW) Mathematical Model UnTRIM Validation Document – Version June 2004 (1.0)
List of Abbreviations
Abbreviation Full Name BAW The Federal Waterways Engineering and Research Institute USGS U. S. Geological Survey r. h. s. right hand side terms of a linear system
Page viii The Federal Waterways Engineering and Research Institute (BAW) Mathematical Model UnTRIM Validation Document – Version June 2004 (1.0)
List of Variables
Name Unit Description C – tracer concentration
Ci k – C for the i-th polygon and k-th layer
CB – prescribed bottom concentration
CB i – CB for i-th polygon
CT – prescribed surface concentration
CT i – CT for i-th polygon 1/2 Cz m /s Chezy coefficient 3
Di k m /s diffusive flux coefficient for i-th polygon and k-th layer 3
D j k m /s diffusive flux coefficient for j-th side and k-th layer
I3 – number of computational prisms £ I3 – I3 at open boundaries (with prescribed water level) Ii – number of computational prisms above polygon i
J3 – number of computational faces £ J3 – J3 at open boundaries (with prescribed water level) Jj – number of computational faces above side j H m total water depth
Hi m H for the i-th polygon
Hj m H for the j-th side
Hmin m minimum allowed water depth H Kh m2/s horizontal eddy diffusivity Kh m2/s Kh for j-th side and k-th layer j k Kv m2/s vertical eddy diffusivity v 2 v
K 1 m /s K at i-th polygon and top/bottom of the k-th layer ¦ i k 2 L m basin length M – index for the surface z-layer
Nc – number of species
Nd – number of source/sink locations
Np – number of polygons/linear tridiagonal systems
Npr – number of red polygons (red-black-sorting) £ Np – Np along open boundaries (with prescribed water level) Ns – number of sides
Nsi – number of internal sides (share two polygons)
Ns f – number of last side with Dirichlet boundary condition Nv – number of vertices
Nz – number of level surfaces Nτ – number of substeps
Page ix The Federal Waterways Engineering and Research Institute (BAW) Mathematical Model UnTRIM Validation Document – Version June 2004 (1.0)
List of Variables (continued)
Name Unit Description 2 Pi m area of the i-th polygon 3
Qi k m /s advective flux coefficient for i-th polygon and k-th layer 3
Q j k m /s advective flux coefficient for j-th side and k-th layer
Si – number of sides for the i-th polygon
S · – sides j for outflow-faces of polygon i, layer k i k S – sides j for inflow-faces of polygon i, layer k i k
νv ∆ ∆
¦ ¦
a j k 1 m abbreviation for t/ z j k 1 ¦ j k 1 c m/s wave celerity 3
di k m /s flux coefficient for polygon i and k-th layer 3
d j k m /s flux coefficient for side j and k-th layer f 1/s Coriolis parameter s fdep – probability of deposition for suspended sediments s fres – probability of resuspension for suspended sediments g m/s2 gravitational acceleration h m bathymetric depth
h j m h for the j-th side
hL m depth for permanently dry land
i – index for the i-th polygon of a grid