Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “unstructured ocean model”

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.

34 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Comparing the Performance of Julia on CPUs versus GPUs and Julia-MPI versus Fortran-MPI: a case study with MPAS-Ocean (Version 7.1)

Abstract. Some programming languages are easy to develop at the cost of slow execution, while others are fast at runtime but much more difficult to write. Julia is a programming language that aims to be the best of both worlds – a development and production language at the same time. To test Julia's utility in scientific high-performance computing (HPC), we built an unstructured-mesh shallow water model in Julia and compared it against an established Fortran-MPI ocean model, the Model for Prediction Across Scales–Ocean (MPAS-Ocean), as well as a Python shallow water code. Three versions of the Julia shallow water code were created: for single-core CPU, graphics processing unit (GPU), and Message Passing Interface (MPI) CPU clusters. Comparing identical simulations revealed that our first version of the Julia model was 13 times faster than Python using NumPy, where both used an unthreaded single-core CPU. Further Julia optimizations, including static typing and removing implicit memory allocations, provided an additional 10–20× speed-up of the single-core CPU Julia model. The GPU-accelerated Julia code was almost identical in terms of performance to the MPI parallelized code on 64 processes, an unexpected result for such different architectures. Parallelized Julia-MPI performance was identical to Fortran-MPI MPAS-Ocean for low processor counts and ranges from 2× faster to 2× slower for higher processor counts. Our experience is that Julia development is fast and convenient for prototyping but that Julia requires further investment and expertise to be competitive with compiled codes. We provide advice on Julia code optimization for HPC systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

The Energy Exascale Earth System Model Version 3: 2. Overview of the Coupled System

The Energy Exascale Earth System Model version 3 (E3SMv3) represents the latest advancement in Earth system modeling developed by the U.S. Department of Energy (DOE). Building upon previous versions, E3SMv3 introduces significant updates across its coupled components to enhance capability and improve fidelity. The atmosphere component incorporates advancements in chemistry, aerosol-cloud interactions, convection, and microphysics. The ocean features a new time-stepping scheme and a higher-resolution unstructured mesh with sub-ice-shelf cavities, while the sea ice model integrates advanced snow and ice physics for more realistic cryospheric simulations. The land model introduces prognostic vegetation dynamics and a new sub-grid topographic treatment of solar radiation. A new tri-grid configuration harmonizes the horizontal grids of the land and river components for improved process coupling. It is enabled by a new non-linear remapping between the atmosphere and land. E3SMv3 underwent extensive testing through a comprehensive simulation campaign, including pre-industrial control, idealized CO 2 experiments, and historical simulations spanning 1850–2024. The model demonstrates significant improvements in simulating the evolution of the historical surface temperature, particularly addressing the “pothole cooling” bias in earlier versions. Reduced aerosol-related forcing contributes to more realistic radiative forcing and better alignment with the observational record. Ocean heat content (OHC) and sea ice trends are also improved as a result.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Multi-decade high-resolution regional hindcasts for wave energy resource characterization in U.S. coastal waters

Long-term, high-resolution, regional wave hindcast datasets were generated using unstructured-grid Simulating WAves Nearshore (SWAN) models for the U.S. coastal waters to support nearshore wave energy development in the U.S. including those bordering U.S. territorial islands. The model domains resolved the entire U.S. exclusive economic zones, with a spatial resolution of approximately 200 m nearshore. The regional SWAN models were driven by the global WAVEWATCH III® model outputs and run for a 42-year period from 1979 to 2020. Extensive model validations were performed using buoy observations and altimeter data. Regional resource characterization was performed based on hindcast data points at 2 km from shore and along the 100 m isobath. Aggregations of wave resource parameters were produced, and spatial and seasonal variations were analyzed for all the regions. Wave resource metrics recommended by international standards, including a 3-hour time series of six resource parameters, hourly frequency- and directionally resolved wave spectra at selected “virtual buoy” locations, and average-annual values of omni-directional wave power, significant wave height, and energy period are publicly disseminated through an Amazon Web Service and a Marine Energy Atlas web application tool to facilitate wave energy research and a wide range of coastal ocean applications.

16 TIDAL AND WAVE POWER↗
DOE OSTI2020

Robust Control of Wave Energy Converters Using Unstructured Uncertainty

In the design of ocean wave energy converters, proper control design is essential to the maximization of the power generation performance for the device. However, in realistic applications, this control design must be undertaken in the presence of model uncertainty. This paper considers the use of robust control theory to optimize the nominal performance for a wave energy converter in stochastic waves, subject to the constraint that the controller be stability-robust to unstructured uncertainties. We formulate the problem as a multi-objective optimal control problem, in which the primary objective is the maximization of power generation for the nominal system, and the competing objective is the norm of the uncertainty input/output channel. This optimal control problem is nonconvex, and we therefore propose an iterative algorithm that can be used to arrive at a local optimal solution. This iterative approach employs the concept of Iterative Convex Overbounding, in the context of the classical Method of Centers. Here, the methodology is demonstrated on a model of a single, buoy- type wave energy converter.

16 TIDAL AND WAVE POWER↗
DOE OSTI2026

High-Resolution Regional Atmosphere–Ocean–Wave Coupled Simulations of Hurricane Henri (2021)

To explore the integrated effects of ocean and ocean surface wave related physical processes on TC simulations, a set of three model simulations is performed. * In experiment 'A', the event is modeled using the Weather Research Forecasting (WRF) model alone with prescribed Sea Surface Temperature (SST) at 6-hour intervals. * In experiment ‘AO,’ WRF is coupled with the Finite Volume Community Ocean Model (FVCOM), enabling variable exchange between atmosphere and ocean, but without considering ocean surface wave-related physical processes. * In experiment ‘AOW’, WRF, FVCOM, and Simulating WAves Nearshore (SWAN) exchange variables with each other every hour through the OASIS3-MCT Coupler to allow direct and indirect atmosphere-ocean-wave interactions. * Observational data are also included in this dataset (Dropsonde, HRD-Radar, NDBC_wave). All simulations are initialized at 18:00 UTC on August 19, 2021, within a domain encompassing the western North Atlantic Ocean. The atmospheric domain features a horizontal resolution of 3 km. The ocean domain, which covers a substantial portion of the WRF ocean domain, employs an unstructured triangular grid with resolutions ranging from 3 km near the coast to 9 km in the open ocean, effectively resolving the complex coastline of the U.S. Northeast Coast. Initial and boundary conditions for the atmosphere model are obtained from the 6-hourly 0.25° NCEP (National Centers for Environmental Prediction) Global Forecast System (GFS; NCEP, 2015) data.

17 WIND ENERGY↗
DOE OSTI2026

High-Resolution Regional Atmosphere–Ocean–Wave Coupled Simulations of Hurricane Henri (2021)

To explore the integrated effects of ocean and ocean surface wave related physical processes on TC simulations, a set of three model simulations is performed. * In experiment 'A', the event is modeled using the Weather Research Forecasting (WRF) model alone with prescribed Sea Surface Temperature (SST) at 6-hour intervals. * In experiment ‘AO,’ WRF is coupled with the Finite Volume Community Ocean Model (FVCOM), enabling variable exchange between atmosphere and ocean, but without considering ocean surface wave-related physical processes. * In experiment ‘AOW’, WRF, FVCOM, and Simulating WAves Nearshore (SWAN) exchange variables with each other every hour through the OASIS3-MCT Coupler to allow direct and indirect atmosphere-ocean-wave interactions. * Observational data are also included in this dataset (Dropsonde, HRD-Radar, NDBC_wave). All simulations are initialized at 18:00 UTC on August 19, 2021, within a domain encompassing the western North Atlantic Ocean. The atmospheric domain features a horizontal resolution of 3 km. The ocean domain, which covers a substantial portion of the WRF ocean domain, employs an unstructured triangular grid with resolutions ranging from 3 km near the coast to 9 km in the open ocean, effectively resolving the complex coastline of the U.S. Northeast Coast. Initial and boundary conditions for the atmosphere model are obtained from the 6-hourly 0.25° NCEP (National Centers for Environmental Prediction) Global Forecast System (GFS; NCEP, 2015) data.

17 WIND ENERGY↗
DOE OSTI2026

High-Resolution Regional Atmosphere–Ocean–Wave Coupled Simulations of Hurricane Henri (2021)

To explore the integrated effects of ocean and ocean surface wave related physical processes on tropical cyclone simulations, a set of three model simulations is performed. * In experiment 'A', the event is modeled using the Weather Research Forecasting (WRF) model alone with prescribed Sea Surface Temperature (SST) at 6-hour intervals. * In experiment ‘AO,’ WRF is coupled with the Finite Volume Community Ocean Model (FVCOM), enabling variable exchange between atmosphere and ocean, but without considering ocean surface wave-related physical processes. * In experiment ‘AOW’, WRF, FVCOM, and Simulating WAves Nearshore (SWAN) exchange variables with each other every hour through the OASIS3-MCT Coupler to allow direct and indirect atmosphere-ocean-wave interactions. * Observational data are also included in this dataset (Dropsonde, HRD-Radar, NDBC_wave). All simulations are initialized at 18:00 UTC on August 19, 2021, within a domain encompassing the western North Atlantic Ocean. The atmospheric domain features a horizontal resolution of 3 km. The ocean domain, which covers a substantial portion of the WRF ocean domain, employs an unstructured triangular grid with resolutions ranging from 3 km near the coast to 9 km in the open ocean, effectively resolving the complex coastline of the U.S. Northeast Coast. Initial and boundary conditions for the atmosphere model are obtained from the 6-hourly 0.25° NCEP (National Centers for Environmental Prediction) Global Forecast System (GFS; NCEP, 2015) data. These CSV files are derived from the NetCDF files in the c0 dataset. Unlike the original format, where geographic coordinates were stored in a separate file, each CSV now embeds the corresponding latitude and longitude values alongside the measured or simulated variables.

17 WIND ENERGY↗
DOE OSTI2023

Topological Relationship–Based Flow Direction Modeling: Mesh–Independent River Networks Representation

River networks are important features in surface hydrology. However, accurately representing river networks in spatially distributed hydrologic and Earth system models is often sensitive to the model's spatial resolution. Specifically, river networks are often misrepresented because of the mismatch between the model's spatial resolution and river network details, resulting in significant uncertainty in the projected flow direction. In this study, we developed a topological relationship-based river network representation method for spatially distributed hydrologic models. This novel method uses (a) graph theory algorithms to simplify real-world vector-based river networks and assist in mesh generation; and (b) a topological relationship-based method to reconstruct conceptual river networks. The main advantages of our method are that (a) it combines the strengths of vector-based and DEM raster-based river network extraction methods; and (b) it is mesh-independent and can be applied to both structured and unstructured meshes. This method paves a path for advanced terrain analysis and hydrologic modeling across different scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A High-resolution Regional Wave Resource Characterization For The U.S. West Coast

Objectives/Scope: Wave resource characterization is a critical step for wave energy converter deployment in the coastal ocean and relies on long-term, high-resolution wave datasets. This study presents a detailed modeling study of the wave resource along the U.S. West Coast (Washington, Oregon, and California), a coastal region that was identified with high wave energy potential in earlier studies. Methods, Procedures, Process: The wave hindcast covers a 32-year period from 1979 to 2010 and is based on a multi-resolution, unstructured-grid SWAN model framework. Model configuration closely follows and meets the requirements recommended by the International Electrotechnical Commission Technical Specification (IEC TS) for wave energy resource assessment and characterization (Class 2 - feasibility study). The model domain covers the entire U.S. Exclusive Economic Zone (EEZ) in the West Coast and has a spatial resolution varying from ~300 m in the nearshore region (20 km from the shoreline) to ~2500 m within the EEZ and ~5000 m at the open boundary, which extends beyond the EEZ. The model was forced by hourly 2-D wave spectra produced by a two-way nested WaveWatch III model, which covers the global ocean domain and the broader U.S. West Coast region domain with spatial resolutions of 0.5 degree and 10 arc-minutes, respectively. Both wave models are forced by hourly, 0.5-degree wind forcing obtained from NCEP’s Climate Forecast System Reanalysis (CFSR) product. Results, Observations, Conclusions: The standard model output for the SWAN model includes 3-hourly output for the six IEC wave resource parameters (e.g., omnidirectional wave power) at each grid point and hourly 2-D spectra at more than 50 NDBC buoys. Extensive model validation was achieved by comparing the six model-predicted IEC parameters with those derived from field observations at representative NDBC buoys. The error statistics indicated the model’s satisfactory performance. Further analyses were conducted to systematically evaluate the temporal and spatial distributions of wave energy potential and wave climate along the U.S. West Coast. Results suggest that Washington and Oregon coasts have similar nearshore wave resource, which is significantly higher than resources in Southern California. Strong seasonal variations are also observed, e.g., high wave energy tends to occur in the winter months. In summary, this study produced the first high-resolution, comprehensive dataset on wave energy distribution along the U.S. West Coast. Novel/Additive Information: The results are being used by the National Renewable Energy Laboratory to update the MHK Atlas, which was originally derived from NOAA’s 4-arc-minute WaveWatch III model output. In addition, the monthly averaged wave energy climatology dataset can be readily shared to support a variety of research and application efforts within the EEZ of the U.S. West Coast.

Wang, Taiping↗
DOE OSTI2024

Total Water Level Prediction at Continental Scale: Coastal Ocean

We demonstrate recent progress made in the simulation of total water level (TWL) at continental scale, using the coastal ocean of US East Coast/Gulf of Mexico coast as an example. A key difference between the continental-scale and small-scale modeling is that the former requires a more accurate vertical datum. Using a geoid-based datum (xGEOID20b), a satellite altimetry product, and a state-of-the-art 3D unstructured-grid model, we significantly improve the accuracy for TWL both near- and off-shore. The average root-mean-square error at all NOAA stations is 14 cm. The non-tidal signals are found to be sensitive to the representation of a large-scale current system near the boundary and extending the domain extent to accommodate this system improves these signals.

Cui, Linlin↗
DOE OSTI2023

Improving the representation of isopycnal mixing in E3SM (Final Report)

The DOE-supported Energy Exascale Earth System Model (E3SM) is a major attempt by the Department of Energy to develop a new Earth System Model using unstructured grids, which allow for the model to concentrate resolution where it is most needed and avoid numerical artifacts experienced by previous Earth System Models (ESMs) using regular grids in which the lateral resolution becomes extremely fine as the resolution shrinks at the poles. This effort required rewriting many of the algorithms previously developed for regular grids. One algorithm that was problematic in the first version of E3SM was the representation of mixing due to turbulent ocean eddies with scales smaller than the model grid. These mesoscale eddies are the primary way by which tracers are stirred along density surfaces in the ocean interior. In particular, the process of isopycnal mixing exchanges fresher, more oxygenated waters from polar regions with saltier, nutrient-rich and oxygen poor waters from tropical regions. Previous work in our group has shown that this mechanism is vitally important for bringing oxygen into poorly ventilated tropical regions (Gnanadesikan et al., 2012, 2013) and plays a significant role in determining the uptake of anthropogenic carbon dioxide (Gnanadesikan et al., 2015). However, in E3SMv1 this process was turned off, as turning it on caused the model to become unstable. Additionally, the rate of isopycnal mixing is directly proportional to a mixing coefficient $A_{Redi}$ whose value varies from less than 400 m 2 /s to 2000 m 2 /s across contemporary climate models. The main thrusts of this proposal therefore were to: 1.) Improve the numerics of isopycnal mixing in E3SM; 2.) Identify ways in which climate models are sensitive to the isopycnal mixing coefficient; and 3.) Develop better representations of the isopycnal mixing coefficient.

54 ENVIRONMENTAL SCIENCES↗
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 OSTI2022

Dissipation and Bathymetric Sensitivities in an Unstructured Mesh Global Tidal Model

Abstract The mechanisms and geographic distribution of global tidal dissipation in barotropic tidal models are examined using a high resolution unstructured mesh finite element model. Mesh resolution varies between 2 and 25 km and is especially focused on inner shelves and steep bathymetric gradients. Tidal response sensitivities to bathymetric changes are examined to put into context response sensitivities to frictional processes. We confirm that the Ronne Ice Shelf dramatically affects Atlantic tides but also find that bathymetry in the Hudson Bay system is a critical control. We follow a sequential frictional parameter optimization process and use TPXO9 data‐assimilated tidal elevations as a reference solution. From simulated velocities and depths, dissipation within the global model is estimated and allows us to pinpoint dissipation at high resolution. Boundary layer dissipation is extremely focused with 1.4% of the ocean accounting for 90% of the total. Internal tide friction is much more distributed with 16.7% of the ocean accounting for 90% of the total. Often highly regional dissipation can impact basin‐scale and even ocean wide tides. Optimized boundary layer friction parameters correlate very well with the physical characteristics of the locality with high friction factors associated with energetic tidal regions, deep ocean island chains, and ice covered areas. Global complex M 2 tide errors are 1.94 cm in deep waters. Total global boundary layer and internal tide dissipation are estimated, respectively, at 1.83 and 1.49 TW. This continues the trend in the literature toward attributing more dissipation to internal tides.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Energy Exascale Earth System Model v2.0.1

First patch release of v2.0.0 Changes since v2.0.0 [Important change] Fix ocean threading bug seen in debug cases on Chrysalis with Intel 20.0.4. Was introduced around time of v2.0.0 tag. Does not change v2.0.0 answers on Chrysalis because those didn't use threading or debugging. [EAM] Add semi-lagrangian tracer transport for theta-l (F90 and C++), add new algorithm for finding tropopause, add DSCREAM to allow v2 and SCREAM settings in same code such as adjust_ps [EAM-MMF] 60L default, allow C++ back end of RRTMGP (EAM too). [EAMxx] add nu-top functionality, fix forcing functor, add ttype9 and dcmip2012 tests 2.1, 2.2, and 3 HOMME: remove obsolete remap algs, option to specify dynamics alg indep of tracer, new sponge layer, add imex tests [ELM] Add topography-based subgrid (topounits), add FATES-ELM Nitro., Phos. and CH4 coupling, add land-use ts for NARRM, add lulc for SSP3 RCP7, Fix nutrient fertilization exp test and carbon isotope flux, Fix xactive lnd dry deposition, add lake water storage option, fix plant hydraulics 2d params, fix carbon budget calc, fix soil nutrient conc. bug, fix mosart dam bug, add test for new ELM, MOSART features, fix bug in O3 dry dep stomatal resistances, fix plant hydraulics restart BFB error, update mkmapdata. [MOSART] fix bug for reading the latitude from an unstructured input file, fix oversat in bubble test. [MPAS-ocean] Add CFC11, CFC12 tracers, add 2D spherical transport tests, fix del4 tracer mixing, add MARBL ocean tracer mixing, modify harmonic analysis options, add GPU port of vmix routines, fix calc of ML-averaged BV freq. [MPAS-seaice] Change extents of initial polar disks for oRRS18to6v3 grid, fix ice BGC with MARBL, update spherical test cases, fix DON coupling, Remove Cf from sea ice constants. [MPAS-landice] add CRYO1850-4xCO2 compset [CIME] add GCP, ANL GCE, Spock, Perlmutter, deprecate config_compilers.xml, fix and clean-up cmake macros, fix slurm bindings, refactor CIME internal testing, cleanup SCORPIO perf data, allow position independent compset naming, [also] update v2 benchmarking suite, extend e3sm_prod with throughput and memory checks

E3SM Project, DOE↗
DOE OSTI2020

High-resolution Hindcasts for U.S. Wave Energy Resource Characterization

The marine and hydrokinetic (MHK) industry is at an early stage of development and has the potential to play a significant role in diversifying the U.S. energy portfolio and reducing the U.S. carbon footprint. Wave energy is the largest among all the U.S. MHK energy resources, which include wave energy, ocean current, tidal-instream, ocean thermal energy conversion, and river-instream. Wave resource characterization is an essential step for regional wave energy assessments, Wave Energy Converter (WEC) project development, site selection and WEC design. The present paper provides an overview of a joint modelling effort by the Pacific Northwest National Laboratory and Sandia National Laboratories on high-resolution wave hindcasts to support the U.S. Department of Energy’s Water Power Technologies Office’s program of wave resource characterization, assessment and classifications in all US coastal regions. Topics covered include the modelling approach, model input requirements, model validation strategies, high performance computing resource requirements, model outputs and data management strategies. Examples of model setup and validation for different regions are provided along with application to development of classification systems, and analysis of regional wave climates. Lessons learned and technical challenges of the long-term, high-resolution regional wave hindcast are discussed.

Yang, Zhaoqing↗
DOE OSTI2025

Characterization of Gulf Stream energy resources: A 30-year high-resolution simulation study

The strength and proximity of the Gulf Stream (GS) along the U.S. coast offers opportunities for marine energy extraction. Leveraging high-resolution (up to 400 m) long-term (30-year) hindcasts based on an unstructured grid model, this study characterizes the energy potential of the GS in two viable locations: the Florida Straits (FS) and Cape Hatteras (CH). Results reveal distinct characteristics in the kinetic energy flux (KEF) and kinetic energy density (KED) for the two regions. The 30-year mean KEF at CH is 27.50 GW, surpassing the 19.74 GW estimated for the FS, which is attributed to the wider cross-section of the GS in CH. Additionally, CH exhibits higher KEF variability due to GS’s larger path shifts and meanders. Estimated KEDs at ADCP site locations reach up to 2,908 W/m 2 in the FS and 1,512 W/m 2 in the CH at a depth of 20 m. High-energy zones are closer to shore in the FS, whereas energy-rich areas at CH are located further offshore. A depth-based KED analysis shows significant decreases across typical energy extraction depths (20 to 100 m), with reductions up to 44%. Finally, this study provides high-resolution resource characterization, forming a foundation for optimized site selection and marine hydrokinetic system design.

16 TIDAL AND WAVE POWER↗