Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stratified”

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 163 records · Page 9

Histopathologic and proteogenomic heterogeneity reveals features of clear cell renal cell carcinoma aggressiveness

Clear cell renal cell carcinomas (ccRCCs) represent ∼75% of RCC cases and account for most RCC-associated deaths. Inter- and intratumoral heterogeneity (ITH) results in varying prognosis and treatment outcomes. To obtain the most comprehensive profile of ccRCC, we perform integrative histopathologic, proteogenomic, and metabolomic analyses on 305 ccRCC tumor segments and 166 paired adjacent normal tissues from 213 cases. Combining histologic and molecular profiles reveals ITH in 90% of ccRCCs, with 50% demonstrating immune signature heterogeneity. High tumor grade, along with BAP1 mutation, genome instability, increased hypermethylation, and a specific protein glycosylation signature define a high-risk disease subset, where UCHL1 expression displays prognostic value. Single-nuclei RNA sequencing of the adverse sarcomatoid and rhabdoid phenotypes uncover gene signatures and potential insights into tumor evolution. In vitro cell line studies confirm the potential of inhibiting identified phosphoproteome targets. This study molecularly stratifies aggressive histopathologic subtypes that may inform more effective treatment strategies.

CPTAC↗

On the effects of reactant stratification and wall curvature in non-premixed rotating detonation combustors

The optimization of non-premixed rotating detonation combustors (RDCs) requires improved understanding of the coupled effects of reactant stratification, fluid property gradients, and complex shock-wave interactions on the detonation wave structure within annular geometries. In the current work, simultaneous orthogonal views of chemiluminescence and hydroxyl planar laser-induced fluorescence (PLIF) are utilized to establish the existence of a dual-wave system characterized by leading and trailing detonation waves that are closely coupled by the local flow physics. These features are persistent over a wide range of mass flow rates and are consistent with prior observations of non-premixed rotating detonations in annular geometries. The detailed instantaneous time sequences are compared with a 3D reactive unsteady Reynolds averaged Navier-Stokes (URANS) simulation to more clearly elucidate the in-situ combustion dynamics and the sensitivity to reactant inlet conditions. It is found that the dual-wave system results from unburned reactants that survive the leading detonation wave in the injector near field and are consumed within a trailing azimuthal reflected-shock combustion (ARSC) zone. By contrast, the injector far field is characterized by rapid mixing due to a sudden drop to subsonic conditions, and the bifurcated detonation wave structure collapses into a stronger, single-wave detonation front with higher overall pressure ratio as compared with the dual-wave system. While each RDC will have different inflow, mixing, and combustion characteristics, the underlying interactions between the stratified reactants and azimuthal wave dynamics identified through the combination of advanced MHz-rate diagnostics and 3D numerical simulations have important implications for the study of detonation wave stability, mode transition, and combustion efficiency in non-premixed annular RDCs.

3D URANS↗

A comprehensive model to capture electrical discharge and spark channel evolution during spark-ignition processes

Emissions reduction through engine efficiency improvements is a priority for automakers who have turned to unconventional engine operation such as highly dilute, boosted, and stratified charge. Given its importance to flame initiation and sustained turbulent flame propagation, reliable and accurate spark ignition models are necessary to design ignition systems that reduce cyclic indicated mean effective pressure (IMEP) variability and increase engine efficiency in these operation modes. In this paper, secondary electric circuit, short-circuit, blowout, and re-strike sub-models are added to the Lagrangian-Eulerian spark ignition (LESI) model to simulate electrical discharge and spark channel elongation in an inert cross-flow combustion vessel. First, the physics of the underlying sub-models are described and the governing equations discussed with the spark channel voltage expression playing a critical role. Then, the experimental and simulations setups are presented. The results section begins with the derivation of spark channel voltage from experimental results. Then, electrical discharge and spark channel elongation simulations in inert flows are carried out using LESI and compared against experimental results. Finally, the results validate the model's ability to accurately predict spark channel elongation, as well as the occurrence of short-circuits, blowouts, re-strikes, and end of discharge.

42 ENGINEERING↗

Detonation structure in the presence of mixture stratification using reaction-resolved simulations

Many investigations of detonation-based combustors have identified reactant mixture inhomogeneity as having a leading-order impact on wave dynamics and combustion efficiency. To examine this phenomenon in a simplified context, an array of two- and three-dimensional channel detonation simulations are conducted in the present work. The reactant mixture consists of stratified fuel and air, wherein the randomly distributed equivalence ratio field features a characteristic stratification length scale. Detailed chemical kinetics are implemented in an adaptive mesh refinement solution framework where the region near the shock front is resolved with Ο (100) cells per representative ZND induction length. The results show that in comparison to baseline cases with uniform reactant mixtures, reactant stratification has a marked impact on the detonation structure. Increasing the stratification length scale increases the size and irregularity of the detonation cells, yielding larger variations in wave speed. Triple point collisions in fuel-rich regions lead to local wave speeds above the notional mean CJ speed, but wave passage through inert regions causes the local wave speed and strength to diminish. Further, conditional statistics show that increasing the stratification length scale increases the variance in pressure and temperature in the primary reaction zone, as well as the variance in heat release over a range of mixture conditions. In addition to the reactant mixture, the impact of the boundary condition behind the detonation is also investigated. The results show that an inflow boundary condition acts to over-drive the wave, leading to higher peak pressures, smaller detonation cells, and increased reactant consumption. On the other hand, cases with a wall behind the wave exhibit weaker waves with lower peak pressures and heat release rates, as well as greater variance in conditional quantities. Comparisons between complementary two- and three-dimensional simulations show reasonable qualitative agreement in wave structure, speed, and conditional statistics.

42 ENGINEERING↗

Non-monotonic behavior of the lateral diffusivity in an adsorbate as a function of the surface coverage

Diffusion dynamics of adsorbate particles confined on a surface strongly depends on the adsorbate surface coverage. Our molecular dynamics simulation reveals that the average lateral diffusivity of the Lennard-Jones adsorbate on a flat substrate varies non-monotonically with the surface coverage. Given that the adsorbate becomes stratified into distinct layers, we propose a novel approach to analyze the lateral diffusion of the adsorbate particles in different layers, via dividing the total lateral diffusion into consecutive lateral diffusion processes, separated by the layer boundary crossing events. A new definition of the lateral diffusivity in the specific layer has been introduced, and it is revealed that the total lateral diffusivity is equal to the sum of lateral diffusivities in different layers. Furthermore, the non-monotonic behavior of the average lateral diffusivity with the surface coverage is found to result from the difference in behavior of the lateral diffusivity between the first and other layers.

36 MATERIALS SCIENCE↗

Direct numerical simulations of turbulent reacting flows with shock waves and stiff chemistry using many-core/GPU acceleration

Compressible reacting flows may display sharp spatial variation related to shocks, contact discontinuities or reactive zones embedded within relatively smooth regions. The presence of such phenomena emphasizes the relevance of shock-capturing schemes such as the weighted essentially non-oscillatory (WENO) scheme as an essential ingredient of the numerical solver. However, these schemes are complex and have more computational cost than the simple high-order compact or non-compact schemes. In this paper, we present the implementation of a seventh-order, minimally-dissipative mapped WENO (WENO7M) scheme in a newly developed direct numerical simulation (DNS) code called KAUST Adaptive Reactive Flows Solver (KARFS). In order to make efficient use of the computer resources and reduce the solution time, without compromising the resolution requirement, the WENO routines are accelerated via graphics processing unit (GPU) computation. The performance characteristics and scalability of the code are studied using different grid sizes and block decomposition. Furthermore, the performance portability of KARFS is demonstrated on a variety of architectures including NVIDIA Tesla P100 GPUs and NVIDIA Kepler K20X GPUs. In addition, the capability and potential of the newly implemented WENO7M scheme in KARFS to perform DNS of compressible flows is also demonstrated with model problems involving shocks, isotropic turbulence, detonations and flame propagation into a stratified mixture with complex chemical kinetics.

97 MATHEMATICS AND COMPUTING↗

Independent and interactive effects of wet bulb globe temperature and air pollution exposures on suicide mortality

Background: Individual components of the ambient environment, such as temperature and air pollution, exist as part of a complex mixture and have been associated with suicide; however, their interactive effects remain poorly understood. This study examined the independent and interactive effects of wet bulb globe temperature (WBGT), nitrogen dioxide (NO 2 ), and fine particulate matter (PM 2.5 ) on suicide mortality. Methods: We identified 7,551 suicide cases in Utah, USA, from 2000 to 2016 and assigned exposure to daily maximum WBGT (sourced from the European Center for Medium-Range Weather Forecasts) and PM 2.5 and NO 2 concentrations (sourced from a national spatiotemporal ensemble model) using decedent’s residential address at the time of death. A case-crossover design with conditional logistic regression was used to estimate the independent and interactive effects of WBGT max , PM 2.5 , and NO 2 on suicide. For exposure windows, we considered single days preceding suicide (lag 0 to 6) and their averages across preceding days (lag 0–1, 0–3, and 0–6). Analyses were stratified by season. Results: We identified a significant association between WBGT max and suicide across all seasons (odds ratio [OR] = 1.05, 95% confidence interval [CI]: 1.01, 1.10; per 5 °C increase on lag 0–3 days). The associations were stronger in the warm season (March 22 to September 21), with ORs and 95% CIs ranging from 1.08 (1.02, 1.15) to 1.20 (1.10, 1.30) per 5 °C increase depending on the lag periods. We observed synergistic interactions between WBGT max and PM 2.5 and NO 2 in the warm season, associated with higher odds of suicide. The associations of WBGT max with suicide were most pronounced at high NO 2 levels. Conclusions: We found evidence of synergistic interactions between WBGT max and PM 2.5 and NO 2 on suicide in the warm season, emphasizing the need for considering the combined effects of heat stress and air pollution in suicide prevention strategies.

Fine particulate matter↗

How to estimate soil organic carbon stocks of agricultural fields? perspectives using ex-ante evaluation

Estimating soil organic carbon (SOC) stocks of agricultural fields has a range of important applications from development of sustainable management practices to monitoring carbon stocks. There are many estimation strategies with the potential for more reliable estimates of SOC stock and more efficient use of soil sampling and analysis resources, especially by leveraging readily available auxiliary information such as remote sensing. However, concrete guidance for strategy selection is lacking. This study narrows this gap with a comparison of strategies for estimating deep SOC stock (0–60 cm) in a prototypical field. Using high density SOC stock measurements and simulation, we built on past studies by 1) ex-ante evaluating a large number of strategy options, 2) using a Bayesian approach to quantify the uncertainty of the comparison, and 3) considering multiple Bayesian models to assess sensitivity to this modeling choice. We found that, using readily available auxiliary information, both balanced and stratified sampling offer substantial improvements over simple random sampling. The auxiliary information most important for this improvement is a Sentinel-2 SOC index = blue / (green × red), followed by the topographic wetness index. We found that these results are robust to the choice of mapping method, but that there is uncertainty in the magnitude of improvement. Here, we recommend future studies implement this Bayesian approach for simulated ex-ante evaluation of SOC stock estimation strategies across more fields to investigate the generalizability of these findings.

54 ENVIRONMENTAL SCIENCES↗

Using optical fibers to examine thermal mixing of liquid sodium in a pool-type geometry

Sodium fast reactors (SFR) are poised to be a leading candidate for the next generation of commercial nuclear reactor deployment. Companies are commissioning designs for SFR increasing the need for experimental and computational analysis to improve the economics and safety of these reactors. Analysis of SFR behavior during reactor transients directly informs the understanding of reactor safety. Particularly, investigating the thermal transients associated with the reactor loss of flow transient allows for informed design decisions to be made regarding reactor safety. During loss of flow transients the flow rate of the coolant is significantly reduced, but not necessarily stopped. The reduction in flow rate results in coolant temperatures exiting the reactor core that can be drastically different than the coolant in upper pool. When these different temperature fluids interact, thermally stratified layers can form causing cyclic thermal fatigue on the reactor vessel potentially leading to pre-material failure. Past experimental studies have lacked high resolution temperature measurements. This research focuses on demonstrating the reliability of novel techniques for acquiring high resolution temperature measurements in a scaled liquid sodium facility. Tests are conducted that simulate postulated loss of flow transients and provide unprecedented spatial temperature distribution measurements.

42 ENGINEERING↗

CFD unified approach under Eulerian–Lagrangian framework for methanol and gasoline direct injection sprays in evaporative and flash boiling conditions

Innovative synthetic fuels for advanced propulsion systems, such as methanol and ammonia, and synthetic blended fuels (E00, E10, and E30), known for their high volatility, are often injected directly into combustion chambers. It follows that Eulerian–Lagrangian spray models need to accurately capture the spray collapse as a consequence of flash boiling onset and be capable of proficiently handling the preferential evaporation of multi-component fuels in evaporative scenarios. So, we performed the assessment of an Eulerian–Lagrangian CFD code for simulating methanol and E00 gasoline blend sprays in both early and late injection conditions involving flash boiling conditions and preferential evaporation. The adoption of an effervescent breakup model and of a non-equilibrium phase transition model for the discrete phase allows the adoption of a setup that is almost completely free from specific constant tuning, especially for what concerns the breakup model. We validated the simulations using experimental PLV maps of methanol and E00 sprays issued from the ECN Spray M injector. The results highlight a significantly different morphology of the methanol spray compared to the E00 one under late injection conditions. Under stratified combustion, low-volatile fuels are likely to be ignited first, and the flame propagates toward the high-volatile fuels. In conclusion, the spray collapse was also correctly reproduced, inducing the presence of a low-pressure zone and modifying the spray morphology.

E00↗

Integration of Risk Survival Measures Estimated From Pre- and Posttreatment Computed Tomography Scans Improves Stratification of Patients With Early-Stage Non-small Cell Lung Cancer Treated With Stereotactic Body Radiation Therapy

To predict overall survival of patients receiving stereotactic body radiation therapy (SBRT) for early-stage non-small cell lung cancer (ES-NSCLC), we developed a radiomic model that integrates risk of death estimates and changes based on pre- and posttreatment computed tomography (CT) scans. We hypothesize this innovation will improve our ability to stratify patients into various oncologic outcomes with greater accuracy.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Forest aboveground biomass estimation through integration of sentinel-2 and PALSAR-2 time series: assessing models trained on GEDI and field inventory benchmarks

Accurate and spatially explicit forest Aboveground Biomass (AGB) mapping through remote sensing is critical for quantifying terrestrial carbon stocks and informing effective forest management strategies. However, AGB estimation in dense forests with complex terrain remains challenging due to satellite sensor signal saturation problem (saturation issue occurs in high biomass forests), structural complexity, and limited ground truth for calibration. This study presents a novel framework that integrates multi-temporal Sentinel-2 optical imagery, ALOS PALSAR-2 Synthetic Aperture Radar (SAR) data, and topographic variables with explainable Machine Learning to map AGB across mountainous forests within subtropical and temperate oceanic climate zones of Mexico. We evaluate the effects of temporal granularity and sensor synergy by comparing multiple temporal inputs and sensor configurations (Sentinel-2, PALSAR-2, and their fusion), and assess model performance using two reference datasets: NASA GEDI LiDAR-derived biomass and Mexico’s National Forest and Soil Inventory (INFyS). Our results showed that models trained on INFyS consistently outperformed those trained on GEDI, highlighting limitations in GEDI’s reliability in biomass estimates within this study region. Furthermore, the integration of Sentinel-2 and PALSAR-2 provided improved predictions compared to single-sensor models, particularly when combined with temporally explicit yearly statistics. The best-performing model, which was trained on INFyS data, and considered both Sentinel-2 and PALSAR-2 yearly statistics, as well as topographic variables, achieved an R2 of 0.64, RMSE of 51.10 Mg/ha, and relative RMSE (rRMSE) of 58.69%. Explainable ML analysis identified Sentinel-2 spectral indices and topographic features as key predictors, while PALSAR-2 metrics provided complementary information, partially mitigating saturation effects in high-biomass areas. Specifically, integrating both sensors substantially improved AGB estimation in high biomass forest (≥200 Mg/ha), yielding 98% gains over optical-only model, with resulting estimates exceeding GEDI L4B by 29% and ESA-CCI-BIOMASS by 174%. Terrain-stratified analysis indicated close agreement with GEDI in low-slope areas, with increasing divergence as slope steepness increased, while estimates remained consistently higher than ESA-CCI-BIOMASS across all slope classes. The proposed approach advances multi-sensor fusion and temporal feature engineering for AGB mapping using open-access satellite datasets, providing a scalable and reproducible framework for annual biomass monitoring in topographically complex mountainous forests. The resulting 25 m resolution biomass product has the potential to provide spatially detailed information for forest monitoring and may support applications in carbon accounting and forest management.

54 ENVIRONMENTAL SCIENCES↗

Interferometric SAR modelling of near surface data to improve geological model in the Surat Basin, Australia

Here we present a study on geophysical inverse modelling for subsurface structural properties of an unconventional hydrocarbon site that was monitored previously by Interferometric Synthetic Aperture Radar (InSAR) technology for surface deformation. A static three-dimensional geomodel along with extracted property maps replicates the depth of each underlying stratigraphic unit and structural feature with the density of each geological layer. We examine the hypothesis that integration of elastic properties of each formation layer with InSAR observations in a stratified elastic medium will lead to a viscoelastic geophysical inverse problem that can be solved to estimate fractional volume change at the reservoir level. Moreover, we examine synthetic scenarios in which the elastic properties of the formations are perturbed before determining the resulting impact on the rate of surface deformation. The results show that although the slope of underlying formations, their density and depth can define the extent and pattern of a deformation signal, their properties have a marginal impact on volumetric change compared to the dense network of shallow depth Coal Seam Gas (CSG) mining wells. Besides, it is also demonstrated that the inversion of InSAR deformation maps can resolve the uncertainties associated with low-resolution seismic interpretation as well as filling the data gaps within seismic acquisitions. A significant contribution of this investigation to the geological basin modelling involves a) introducing a remote and non-invasive technology such as InSAR to improve geophysical mapping of subsurface structures such as faults in areas with sparse or no reflective seismic information, and b) applying a multi-layer viscoelastic geophysical source model for an unconventional hydrocarbon reservoir such as CSG.

58 GEOSCIENCES↗

Morphological and Phenotypic Characteristics of the Bovine Nasopharyngeal Mucosa and Associated Lymphoid Tissue

The mammalian nasopharynx is an anatomically complex region of the upper respiratory tract that directly communicates with the nasal cavity, laryngopharynx, oesophagus and trachea. The nasopharyngeal mucosa contains moderate quantities of mucosa-associated lymphoid tissue (MALT) that is appropriately located for immunological sampling but also creates vulnerability to pathogens. In recent years, the nasopharynx has been inculpated in the pathogenesis of important diseases of cattle (foot-and-mouth disease) and humans (COVID-19), yet the tissue has never been described in detail in any species. In order to characterize the morphology and cellular composition of the bovine nasopharynx, samples of mucosa were collected from the nasopharynx of five 8–13-month-old steers and examined using light microscopy, immunohistochemistry and multichannel immunofluorescence. Morphologically, the nasopharyngeal epithelium was highly heterogeneous, with a continuum ranging from stratified squamous epithelium to highly attenuated, follicle-associated epithelium (FAE). Distribution of MALT was similarly regionally variable ranging from absent to clusters of multiple lymphoid follicles. Phenotypic characterization demonstrated dense distributions of dendritic cells and T lymphocytes surrounding lymphoid follicles, which comprised mostly B lymphocytes. The FAE overlaying the lymphoid follicles also contained higher numbers of dendritic cells and lymphocytes compared with the adjacent non-lymphoid epithelium, although cytotoxic T cells were notably scarce in the FAE. The bovine nasopharyngeal lymphoid tissue had comparable elements to other MALTs with specific differences that may help to elucidate the pathogenesis of infectious agents that have specific tropism for this tissue.

59 BASIC BIOLOGICAL SCIENCES↗

Manganese in drinking-water reservoirs: a multi-disciplinary review of current issues, biogeochemical controls, and oxygenation-based management

Decreased water quality and increased treatment costs due to excess manganese (Mn) in drinking-water supplies are critical issues globally. To combat on-going and emerging taste and odour issues with Mn and other contaminants (e.g., algal toxins), many utilities are using engineered oxygenation or aeration (EOA) systems to improve water quality in lakes and reservoirs. Resultant shifts in key biogeochemical and physical processes are still poorly understood, often leading to inefficiently managed systems. Paired with knowledge gaps regarding environmental drivers of Mn and the complexity of Mn redox kinetics, Mn problems persist. This review presents the state of current research in areas critical to optimisation of EOA mitigation of Mn, with focus on i) Mn biogeochemical cycling within drinking-water reservoirs; ii) influences of local catchment geology, hydrology, and land use on Mn dynamics; and iii) Mn management using different EOA approaches. The importance of considering the combined implications of these factors for successful Mn management in reservoirs is highlighted by an evaluation of relevant field-based studies; a wide range in EOA performance is observed, from a 97 % decrease in soluble Mn up to a ∼400 % increase in total Mn. Despite the breadth of studies that consider Mn in water-supply systems, there are still several areas of research which warrant further investigation, including: the influence of natural sources and anthropogenic activities on Mn within a given catchment, Mn speciation and transport in stratified and destratified lakes and reservoirs, and optimal site-specific EOA strategies for Mn mitigation.

Aeration↗

Examining the effects of soil entrainment during nuclear cloud rise on fallout predictions using a multiscale atmospheric modeling framework

Current operational models for nuclear cloud rise over land were developed and validated using observations from shallow-buried or surface detonations, where lofted soil quickly mixed with fission products from the detonation. These models poorly predict fallout from elevated detonations near the fallout-free height of burst (FFHOB), where interactions with the ground are limited and the mixing of fission products and lofted soil is incomplete. Fallout-free is a misnomer at this HOB, as fallout was observed in these cases, but was below the levels of concern, especially off-grounds of the nuclear test site. To correctly characterize and model fallout from detonations near the FFHOB, models must be developed which can capture the stratified nature of the particle and activity-size distributions within the cloud. Previously, it was shown that the Weather Research and Forecasting (WRF) model can accurately simulate nuclear cloud rise for airbursts with little to no ground interactions (Arthur et al., 2021). That work is expanded here by (1) using a radiation-hydrodynamics code to improve the fireball initialization in WRF, (2) further developing an aerosol package from WRF-Chem to simulate lofted soil, and (3) combining the WRF cloud rise simulations with the operational models used at the National Atmospheric Release Advisory Center (NARAC) for fallout modeling. Using this combination of codes, the Upshot-Knothole Grable detonation, which was just below the FFHOB, is simulated from seconds after detonation through cloud rise and fallout, and results are compared to historical test data. Here, the results show improved prediction of dose rate and highlight the need to correctly characterize the entrainment of material into the cloud and the subsequent mixing of fission products with entrained material.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Intercomparison of flood inundation models across land use types and hydrological flood stages

Flood Inundation Mapping (FIM) model selection is a key operational decision because accurate, rapid mapping underpins early warning and resource allocation. FIM performance is context-dependent and can vary with hydrograph phase, land-use/land-cover (LULC), and the evaluation benchmark. Intercomparison studies typically assess a single near-peak snapshot against one reference dataset. Here, we provide a context-stratified intercomparison across (i) multiple hydrograph phases, (ii) LULC classes, and (iii) benchmark types, for five FIM approaches spanning a wide range of physical complexity and operational cost (TRITON, LISFLOOD-FP, HEC-RAS 2D, ARC-Curve2Flood, and OWP HAND-FIM). We use the Hurricane Matthew flood (2016) in the Neuse River Basin, North Carolina, USA, as a case study. Using high-resolution remote sensing-derived flood inundation maps, hand-labeled points, and building footprints, we assess model skill across two rising and two falling hydrograph limbs and across major LULC types. Results show that model rankings shift systematically across contexts: LISFLOOD-FP ranks highest in three of four flood phases, while TRITON leads during one rising limb phase; LISFLOOD-FP performs best in vegetated areas, whereas HEC-RAS improves relative performance in agricultural and urban areas; and benchmark choice influences conclusions, with LISFLOOD-FP performing best for flooded-building detection in the late falling limb, while TRITON ranks highest against hand-labeled points. We also report representative wall-clock runtimes for each workflow to provide use-case context for operational feasibility. Together, these results offer transferable guidance for model selection and for designing large-scale, benchmark-aware FIM intercomparison studies.

Nikrou, Parvaneh [University of Alabama]↗

A nonhydrostatic formulation for MPAS-Ocean

The Model for Prediction Across Scales-Ocean (MPAS-Ocean) is an open-source, global ocean model and is one component of a family of climate models within the MPAS framework, including atmosphere, sea-ice, and land-ice models. Here, in this work, a new formulation for the ocean model is presented that solves the nonhydrostatic, incompressible Boussinesq equations on an unstructured, staggered, z-level grid. The introduction of this nonhydrostatic capability is necessary for the resolution of internal wave dynamics and large eddy simulations. Compared to the standard, hydrostatic formulation, a nonhydrostatic pressure solver and a vertical momentum equation are added, where the PETSc (Portable Extensible Toolkit for Scientific Computation) library is used for the inversion of a large sparse system for the nonhydrostatic pressure. Numerical results on a stratified seiche, internal solitary wave, overflow and lock-exchange test cases are presented, and the parallel efficiency of the code is evaluated using up to 1024 processors.

3D Poisson equation↗