Search NASA⌕ Search

SEARCH · Search NASA

Results for “Unstructured Grid”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Algebraic multigrid preconditioning of the Hessian in optimization constrained by a partial differential equation

Summary We construct an algebraic multigrid (AMG) based preconditioner for the reduced Hessian of a linear‐quadratic optimization problem constrained by an elliptic partial differential equation. While the preconditioner generalizes a geometric multigrid preconditioner introduced in earlier works, its construction relies entirely on a standard AMG infrastructure built for solving the forward elliptic equation, thus allowing for it to be implemented using a variety of AMG methods and standard packages. Our analysis establishes a clear connection between the quality of the preconditioner and the AMG method used. The proposed strategy has a broad and robust applicability to problems with unstructured grids, complex geometry, and varying coefficients. The method is implemented using the Hypre package and several numerical examples are presented.

Barker, Andrew T.↗
DOE OSTI2024

ExaWind: Open‐source CFD for hybrid‐RANS/LES geometry‐resolved wind turbine simulations in atmospheric flows

Abstract Predictive high‐fidelity modeling of wind turbines with computational fluid dynamics, wherein turbine geometry is resolved in an atmospheric boundary layer, is important to understanding complex flow accounting for design strategies and operational phenomena such as blade erosion, pitch‐control, stall/vortex‐induced vibrations, and aftermarket add‐ons. The biggest challenge with high‐fidelity modeling is the realization of numerical algorithms that can capture the relevant physics in detail through effective use of high‐performance computing. For modern supercomputers, that means relying on GPUs for acceleration. In this paper, we present ExaWind, a GPU‐enabled open‐source incompressible‐flow hybrid‐computational fluid dynamics framework, comprising the near‐body unstructured grid solver Nalu‐Wind, and the off‐body block‐structured‐grid solver AMR‐Wind, which are coupled using the Topology Independent Overset Grid Assembler. Turbine simulations employ either a pure Reynolds‐averaged Navier–Stokes turbulence model or hybrid turbulence modeling wherein Reynolds‐averaged Navier–Stokes is used for near‐body flow and large eddy simulation is used for off‐body flow. Being two‐way coupled through overset grids, the two solvers enable simulation of flows across a huge range of length scales, for example, 10 orders of magnitude going from O(μm) boundary layers along the blades to O(10 km) across a wind farm. In this paper, we describe the numerical algorithms for geometry‐resolved turbine simulations in atmospheric boundary layers using ExaWind. We present verification studies using canonical flow problems. Validation studies are presented using megawatt‐scale turbines established in literature. Additionally presented are demonstration simulations of a small wind farm under atmospheric inflow with different stability states.

17 WIND ENERGY↗
DOE OSTI2015

Particle tracking approach for transport in three-dimensional discrete fracture networks: Particle tracking in 3-D DFNs

The discrete fracture network (DFN) model is a method to mimic discrete pathways for fluid flow through a fractured low-permeable rock mass, and may be combined with particle tracking simulations to address solute transport. However, experience has shown that it is challenging to obtain accurate transport results in three-dimensional DFNs because of the high computational burden and difficulty in constructing a high-quality unstructured computational mesh on simulated fractures. We present a new particle tracking capability, which is adapted to control volume (Voronoi polygons) flow solutions on unstructured grids (Delaunay triangulations) on three-dimensional DFNs. The locally mass-conserving finite-volume approach eliminates mass balance-related problems during particle tracking. The scalar fluxes calculated for each control volume face by the flow solver are used to reconstruct a Darcy velocity at each control volume centroid. The groundwater velocities can then be continuously interpolated to any point in the domain of interest. The control volumes at fracture intersections are split into four pieces, and the velocity is reconstructed independently on each piece, which results in multiple groundwater velocities at the intersection, one for each fracture on each side of the intersection line. This technique enables detailed particle transport representation through a complex DFN structure. Verified for small DFNs, the new simulation capability enables numerical experiments on advective transport in large DFNs to be performed. As a result, we demonstrate this particle transport approach on a DFN model using parameters similar to those of crystalline rock at a proposed geologic repository for spent nuclear fuel in Forsmark, Sweden.

58 GEOSCIENCES↗
DOE OSTI2020

Analysis of the SBP-SAT Stabilization for Finite Element Methods Part I: Linear Problems

In the hyperbolic community, discontinuous Galerkin (DG) approaches are mainly applied when finite element methods are considered. As the name suggested, the DG framework allows a discontinuity at the element interfaces, which seems for many researchers a favorable property in case of hyperbolic balance laws. On the contrary, continuous Galerkin methods appear to be unsuitable for hyperbolic problems and there exists still the perception that continuous Galerkin methods are notoriously unstable. To remedy this issue, stabilization terms are usually added and various formulations can be found in the literature. However, this perception is not true and the stabilization terms are unnecessary, in general. In this paper, we deal with this problem, but present a different approach. We use the boundary conditions to stabilize the scheme following a procedure that are frequently used in the finite difference community. Here, the main idea is to impose the boundary conditions weakly and specific boundary operators are constructed such that they guarantee stability. This approach has already been used in the discontinuous Galerkin framework, but here we apply it with a continuous Galerkin scheme. No internal dissipation is needed even if unstructured grids are used. Further, we point out that we do not need exact integration, it suffices if the quadrature rule and the norm in the differential operator are the same, such that the summation-by-parts property is fulfilled meaning that a discrete Gauss Theorem is valid. This contradicts the perception in the hyperbolic community that stability issues for pure Galerkin scheme exist. In numerical simulations, we verify our theoretical analysis.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

A modeling study of coastal inundation induced by storm surge, sea-level rise, and subsidence in the Gulf of Mexico

The northern coasts of the Gulf of Mexico are highly vulnerable to the direct threats of climate change, such as hurricane-induced storm surge, and such risks can be potentially exacerbated by land subsidence and global sea level rise. This paper presents an application of a coastal storm surge model to study the coastal inundation process induced by tide and storm surge, and its response to the effects of land subsidence and sea level rise in the northern Gulf coast. An unstructured-grid Finite Volume Coastal Ocean Model was used to simulate tides and hurricane-induced storm surges in the Gulf of Mexico. Simulated distributions of co-amplitude and co-phase of semi-diurnal and diurnal tides are in good agreement with previous modeling studies. The storm surges induced by four historical hurricanes (Rita, Katrina, Ivan and Dolly) were simulated and compared to observed water levels at National Oceanic and Atmospheric Administration tide stations. Effects of coastal subsidence and future global sea level rise on coastal inundation in the Louisiana coast were evaluated using a parameter “change of inundation depth” through sensitivity simulations that were based on a projected future subsidence scenario and 1-m global sea level rise by the end of the century. Model results suggested that hurricane-induced storm surge height and coastal inundation could be exacerbated by future global sea level rise and subsidence, and that responses of storm surge and coastal inundation to the effects of sea level rise and subsidence are highly nonlinear and vary on temporal and spatial scales.

Yang, Zhaoqing↗
DOE OSTI2023

Simulation of compound flooding in Japan using a nationwide model

For this study, a high-resolution, summit-to-sea unstructured-grid model was used to simulate two compound flooding events in different geomorphic settings in Japan: the July 2012 flood in Kumamoto city and July 2018 flood in Okayama City, both caused by a torrential rainfall during seasonal Meiyu front but otherwise exhibiting quite different compound flood characteristics. The model shows good performance in simulating flooding extent; e.g., the Hit Rate for the inundation event in 2018 exceeds 0.9. Sensitivity tests were conducted to determine the roles/significance of levee structures along major rivers in coastal regions. Our results indicate that riverine and oceanic factors as well as levees significantly contribute to the observed flooding extents, and we are able to quantify the contributions from each factor at different locations. The model has been implemented operationally as the first nationwide 3-day compound flooding forecast system for Japan.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

A minimally-dissipative low-Mach number solver for complex reacting flows in Open$\mathrm{FOAM}$

Large eddy simulation (LES) has become the de-facto computational tool for modeling complex reacting flows, especially in gas turbine applications. However, readily usable general-purpose LES codes for complex geometries are typically academic or proprietary/commercial in nature. In this report the objective is to develop and disseminate an open source LES tool for low-Mach number turbulent combustion using the OpenFOAM framework. In particular, a collocated-mesh approach suited for unstructured grid formulation is provided. Unlike other fluid dynamics models, LES accuracy is intricately linked to so-called primary and secondary conservation properties of the numerical discretization schemes. This implies that although the solver only evolves equations for mass, momentum, and energy, the implied discrete equation for kinetic energy (square of velocity) should be minimally-dissipative. Here, a specific spatial and temporal discretization is imposed such that this kinetic energy dissipation is minimized. The method is demonstrated using manufactured solutions approach on regular and skewed meshes for a canonical flow and a lab-scale turbulent flow problem.

42 ENGINEERING↗
DOE OSTI2025

Direct Discontinuous Galerkin methods for the reacting multi-component flow equations

The Direct Discontinuous Galerkin (DDG (Liu and Yan, 2008)) method and a counterpart with Interface Correction (DDGIC (Danis and Yan, 2022)) are extended to compute diffusion terms that arise when solving the compressible multi-component flow equations in thermochemical nonequilibrium. Thermodynamic properties, transport properties, chemical reaction rates, and energy exchange terms are computed using Mutation++ (Scoggins et al., 2020). The DG method is applied on unstructured grids, where the accuracy and convergence rates can be sensitive to the numerical method chosen for parabolic terms. A method for determining the homogeneity tensor of the flow equations required for DDGIC is shown. The convergence properties of the DDG methods are studied and compared to the Interior Penalty (IP) method. A number of numerical experiments are conducted to assess the accuracy and performance of the method. The numerical results and convergence studies indicate that DDG and DDGIC provide accurate solutions and perform well for general flows in thermochemical nonequilibrium.

Diffusion↗
DOE OSTI2015

Understanding the Flushing Capability of Bellingham Bay and Its Implication on Bottom Water Hypoxia

In this study, an unstructured-grid finite-volume coastal ocean model (FVCOM) was used to simulate hydrodynamic circulation and assess the flushing capability in Bellingham Bay, Washington, USA. The model was reasonably calibrated against field observations for water level, velocity and salinity, and was further used to calculate residence time distributions in the study site. The model results suggest that, despite the large tidal ranges (~4 m during spring tide), tidal currents are relatively weak in Bellingham Bay with surface currents generally below 0.5 m/s. The local residence time in Bellingham Bay varies from to near zero to as long as 15 days, depending on the location and river flow condition. In general, Bellingham Bay is a well-flushed coastal embayment affected by freshwater discharge, tides, wind, and density-driven circulation. The basin-wide global residence time ranges from 5-7 days. The model results also provide useful information on possible causes of the emerging summertime hypoxia problem in the north central region of Bellingham Bay. It was concluded that the formation of the bottom hypoxic water should result from the increased consumption rate of oxygen in the bottom oceanic inflow with low dissolved oxygen by organic matters accumulated at the regions characterized with relatively long residence time in summer months.

Bellingham Bay↗
DOE OSTI2019

Modeling Analysis of the Swell and Wind-Sea Climate in the Salish Sea

Accurate model hindcast of wave climate in complex estuarine systems is challenging because high-resolution wave model and wind forcing are required. In this study, a modeling approach using the unstructured-grid Simulating WAves Nearshore and wind product from a high-resolution regional Weather Research and Forecasting hindcast was used to simulate the swell and wind-sea climate in the Salish Sea, a large estuary with many interconnected waterways on the Pacific Northwest coast of North America. The model hindcast was validated with observed data at four wave buoys in the Salish Sea. Spatial distribution and seasonal variations of wave climate in the Salish Sea were analyzed. Of the three major basins in the Salish Sea, the Strait of Juan de Fuca has the largest waves and is dominated by swells propagated from the Pacific Ocean. Significant wave heights in the Strait of Georgia have spatial and seasonal distribution patterns similar to those found in the Strait of Juan de Fuca. Waves in Puget Sound are small and primarily dominated by the wind-sea climate. Strong seasonal variations are observed in the Strait of Juan de Fuca and Strait of Georgia, but there is little seasonality of wave climate in Puget Sound. The high-resolution wave hindcast conducted in this study provides a comprehensive and important data set for better understanding the role of wave climate in coastal processes and natural hazards assessment in the Salish Sea.

wave climate, salish sea, modeling, swell, sea sta↗
DOE OSTI2019

Cross-scale modeling of storm surge, tide, and inundation in Mid-Atlantic Bight and New York City during Hurricane Sandy, 2012

Driven by high-resolution NAM (North America Mesoscale) and ECMWF (European Centre for Medium-Range Weather Forecasts) atmospheric fields, a 3-D unstructured grid SCHISM (Semi-implicit Cross-scale Hydroscience Integrated System Model) was applied to simulate total water level in the Mid-Atlantic Bight during Hurricane Sandy in 2012. The simulated storm surge, tide, significant wave height, and wave peak period results compared favorably with NOAA observations along US East Coast, Long Island Sound, and New York Harbor. The maximum total water level at The Battery in the New York Harbor was accurately simulated with an absolute error of less than 0.08 m (out of 2.9 m peak surge) and a timing difference less than 10 min, which includes the effects on the order of 0.1–0.3 m (5–10%) from the coastal setup induced by the surface waves. The scenario comparisons of (1) “NAM” versus “ECMWF”, (2) “2-D” versus “3-D”, and (3) “with” versus “without” wind wave model were examined. The 3-D barotropic model forced by ECMWF including the effects of wind waves performs the best, attributed to wave-induced radiation stress and reduction of bottom stress when the 3-D version is used. Simultaneously, the ELCIRC-Sub model, using a 5-m topography/bathymetry sub-grid and a regular 200-m resolution finite volume computational model grid, was developed to simulate the street-level inundation in New York City. The momentum and mass fluxes calculated by the coarser base grid model were effectively coupled with the sub-grid so that running a full-blown high-resolution base model is not required. The ELCIRC-Sub model uses MPI (Message Passing Interface) parallel computing to enlarge the coverage of the land surface and an efficient non-linear solver to improve the accuracy of the wetting-and-drying processes. The temporal comparisons of modeled water level with NOAA's tidal stations and USGS' rapid-deployed gauges showed overall performance with an average error of 0.1 m. Particularly, ELCIRC-Sub captured the profile of the highest peak surge (3.9 m), which rose 2 m within 3 hour at Kings Point, NY. The spatial comparisons of the modeled peak water level at 80 surveyed locations showed an average error less than 0.13 m. The modeled maximum inundation extent also matched well with an 80% hit rate against FEMA's Hurricane Sandy maximum flooding extent.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Explicit quantification of residence and flushing times in the Salish Sea using a sub-basin scale shoreline resolving model

The Salish Sea, located in the Pacific Northwest region of North America has complex currents and circulation features distributed over numerous interconnected deep basins with islands. Increased risk of exposure to oil spills and untreated wastewater from maritime emergencies and treatment plant failures have led to a need for quantifying residence and flushing characteristics at a sub-basin scale using a shoreline resolving hydrodynamic model. An unstructured grid model of the Salish Sea was developed using FVCOM with a ≈75-100m shoreline resolution. In addition to 7 tides and 23 salinity and temperature monitoring stations, an extensive currents data set from 135 stations collected over a span of three years was used for skill assessment and validation. Explicit forward computations were then conducted to define and quantify residence and flushing times in various sub-basins of interest using (a) Lagrangian particles and (b) Numerical/virtual dye experiments. The results in most basins show expected seasonal variability with longer flushing time associated with summer lower tides and lower freshwater inflows. However, contrary to expectation, flushing time is significantly longer in wintertime in large fjord-like basins such as Hood Canal (≈138 days), likely due to increased stratification and reduced mixing. The flushing time for the Puget Sound region of the Salish Sea is ≈115 days, while Georgia Basin is 240 days when analyzed as stand-alone basins with zero background concentrations. When examined as part of the flushing of the entire Puget Sound filled with virtual dye, the compounded flushing times for embedded sub-basins can be significantly longer in order of magnitudes and largely dictated by the flushing time of Puget Sound. The computed residence and flushing time scales tabulated over 36 sub-basins provide an improved understanding of water renewal in the system, informing pollution management actions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Characteristics and Variability of the Nearshore Wave Resource on the U.S. West Coast

Characterizing the nearshore wave resource at a regional scale poses a great challenge because high-quality wave data are required. This paper presents a detailed analysis of nearshore wave climate and resource characterization on the U.S. West Coast based on a 32-year, regional wave hindcast. Resource characterization closely followed the International Electrotechnical Commission standards. The wave hindcast data were generated from an unstructured-grid SWAN model with approximately 300 m grid resolution in the nearshore region. Statistics of six IEC wave resource parameters at the 50 m water depth were calculated for four sub-regions along the coast, including Washington, Oregon, and Northern and Southern California. Results indicate that the Washington and Oregon coast has a similar nearshore wave resource. The wave resource in Southern California is significantly less than that of the Washington and Oregon coast. Alongshore and cross-shore variabilities of wave power were assessed. Temporal variability, at seasonal, inter-annual, and decadal scales were analyzed. Overall, strong seasonal variation is observed; high wave energy occurs in the winter months and calm seas in the summer. Finally, extreme environmental conditions, i.e., the 100-year extreme environmental contours, for different locations along the West Coast were investigated.

wave resource characterization, nearshore, west co↗
DOE OSTI2017

On the utility of GPU accelerated high-order methods for unsteady flow simulations: A comparison with industry-standard tools

First- and second-order accurate numerical methods, implemented for CPUs, underpin the majority of industrial CFD solvers. Whilst this technology has proven very successful at solving steady-state problems via a Reynolds Averaged Navier–Stokes approach, its utility for undertaking scale-resolving simulations of unsteady flows is less clear. High-order methods for unstructured grids and GPU accelerators have been proposed as an enabling technology for unsteady scale-resolving simulations of flow over complex geometries. In this study we systematically compare accuracy and cost of the high-order Flux Reconstruction solver PyFR running on GPUs and the industry-standard solver STAR-CCM+ running on CPUs when applied to a range of unsteady flow problems. Specifically, we perform comparisons of accuracy and cost for isentropic vortex advection (EV), decay of the Taylor–Green vortex (TGV), turbulent flow over a circular cylinder, and turbulent flow over an SD7003 aerofoil. We consider two configurations of STAR-CCM+: a second-order configuration, and a third-order configuration, where the latter was recommended by CD-adapco for more effective computation of unsteady flow problems. Results from both PyFR and STAR-CCM+ demonstrate that third-order schemes can be more accurate than second-order schemes for a given cost e.g. going from second- to third-order, the PyFR simulations of the EV and TGV achieve 75× and 3× error reduction respectively for the same or reduced cost, and STAR-CCM+ simulations of the cylinder recovered wake statistics significantly more accurately for only twice the cost. Moreover, advancing to higher-order schemes on GPUs with PyFR was found to offer even further accuracy vs. cost benefits relative to industry-standard tools.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Analytical and variational numerical methods for unstable miscible displacement flows in porous media

The miscible displacement of one fluid by another in a porous medium has received considerable attention in subsurface, environmental and petroleum engineering applications. When a fluid of higher mobility displaces another of lower mobility, unstable patterns – referred to as viscous fingering – may arise. Their physical and mathematical study has been the object of numerous investigations over the past century. In writing this paper the objective is to present a review of these contributions with particular emphasis on variational methods. These algorithms are tailored to real field applications thanks to their advanced features: handling of general complex geometries, robustness in the presence of rough tensor coefficients, low sensitivity to mesh orientation in advection dominated scenarios, and provable convergence with fully unstructured grids.

42 ENGINEERING↗
DOE OSTI2018

A robust and efficient finite volume method for compressible inviscid and viscous two-phase flows [A robust and efficient finite volume method for compressible viscous two-phase flows]

A robust and efficient density-based finite volume method is developed for solving the six-equation single pressure system of two-phase flows at all speeds on hybrid unstructured grids. Unlike conventional approaches where an expensive exact Riemann solver is normally required for computing numerical fluxes at the two-phase interfaces in addition to AUSM-type fluxes for single-phase interfaces in order to maintain stability and robustness in cases involving interactions of strong pressure and void-fraction discontinuities, a volume-fraction coupling term for the AUSM + -up fluxes is introduced in this work to impart the required robustness without the need of the exact Riemann solver. The resulting method is significantly less expensive in regions where otherwise the Riemann solver would be invoked. A transformation from conservative variables to primitive variables is presented and the primitive variables are then solved in the implicit method in order for the current finite volume method to be able to solve, effectively and efficiently, low Mach number flows in traditional multiphase applications, which otherwise is a great challenge for the standard density-based algorithms. Here, a number of benchmark test cases are presented to assess the performance and robustness of the developed finite volume method for both inviscid and viscous two-phase flow problems. The numerical results indicate that the current density-based method provides an attractive and viable alternative to its pressure-based counterpart for compressible two-phase flows at all speeds.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Interface tracking simulations of bubbly flows in PWR relevant geometries

[en] Highlights: • Simulations were performed for turbulent bubbly flows in PWR subchannel geometry. • Liquid turbulence is fully resolved by direct numerical simulation approach. • Bubble behavior is captured using level-set interface tracking method. • Time-averaged single- and two-phase turbulent flow statistical quantities are obtained. - Abstract: The advances in high performance computing (HPC) have allowed direct numerical simulation (DNS) approach coupled with interface tracking methods (ITM) to perform high fidelity simulations of turbulent bubbly flows in various complex geometries. In this work, we have chosen the geometry of the pressurized water reactor (PWR) core subchannel to perform a set of interface tracking simulations (ITS) with fully resolved liquid turbulence. The presented research utilizes a massively parallel finite-element based code, PHASTA, for the subchannel geometry simulations of bubbly flow turbulence. The main objective for this research is to demonstrate the ITS capabilities in gaining new insight into bubble/turbulence interactions and assisting the development of improved closure laws for multiphase computational fluid dynamics (M-CFD). Both single- and two-phase turbulent flows were studied within a single PWR subchannel. The analysis of numerical results includes the mean gas and liquid velocity profiles, void fraction distribution and turbulent kinetic energy profiles. Two sets of flow rates and bubble sizes were used in the simulations. The chosen flow rates corresponded to the Reynolds numbers of 29,079 and 80,775 based on channel hydraulic diameter (D_h) and mean velocity. The finite element unstructured grids utilized for these simulations include 53.8 million and 1.11 billion elements, respectively. This has allowed to fully resolve all the turbulence scales and the deformable interfaces of individual bubbles. For the two-phase flow simulations, a 1% bubble volume fraction was used which resulted in 17 bubbles in the smaller case and 262 bubbles in the larger case. In the larger simulation case the size of the resolved bubbles is 0.65 mm in diameter, and the bulk mesh cell size is about 30 microns. Those large-scale simulations provide new level of details previously unavailable and were enabled by the excellent scaling performance of our two-phase flow solver and access to the state-of-the-art supercomputing resources. The presented simulations used up to 256 thousand processing threads on the IBM BG/Q supercomputer “Mira” (Argonne National Laboratory).

Nuclear Science & Technology↗
DOE OSTI2024

Debunking common myths in coastal circulation modeling

Despite tremendous progress in algorithm development, computational efficiency and transition into operations over the past two decades, coastal modeling still lacks scientific rigor due to proliferation of many ‘gray’ areas related to various modeling choices made by modelers. Here, in this paper, we propose some guiding principles for the modeling community to improve performance, and we also debunk commonly held myths that make the coastal modeling lack rigor. Using our own experience in developing seamless cross-scale unstructured-grid based models for the past two decades, we describe in unprecedented detail the end-to-end modeling process (i.e., from digital elevation models (DEMs) to mesh generation to post analysis), and demonstrate that defensible modeling is within reach for any end user by following three guiding principles: (1) Bathymetry is a first order forcing in coastal domains and thus should be respected in all aspects of modeling; (2) Oceanographic processes are driven across multiple spatial scales and so models should enable appropriate resolution as needed; and (3) Model assessment should focus on physical processes. Through qualitative and quantitative model assessments, we demonstrate the fundamental role played by bathymetry/topography as embedded in DEMs in making the results defensible, which is unfortunately glossed over in many modeling studies. Focusing on process-based assessment simplifies the calibration process. A major conclusion of this work is that model developers and operators should maximize the scientific rigor for in silico oceanography by avoiding some common pitfalls that rely on error compensation at the expense of representation of physical system processes. We present some best practice procedures for defensive and trustworthy numerical modeling.

54 ENVIRONMENTAL SCIENCES↗