SEARCH · Engineering Papers

Results for “3-D”

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.

SourceResearch connectionFrom yearThrough year

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

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Investigation of fuel film formation and soot emissions in a GDI engine during cold-start with Split-injection strategies

This study investigates the impact of fuel-film formation on engine-out soot emissions in a gasoline direct injection (GDI) engine under cold-start conditions. Split-injection strategies were applied by varying the number of injections, injection duration, and total fuel quantity to affect wall wetting and control the average in-cylinder equivalence ratio. A combined experimental and numerical approach was employed to analyze fuel-film deposition, combustion efficiency, and engine-out soot and unburnt hydrocarbons (UHC) emissions. In particular, fuel film distribution estimated by means of non-reacting, 3-D, computational fluid dynamics (CFD) simulations, together with experimentally measured soot data, were used to investigate fuel film formation and its role in soot generation. In the experiments, a skip-fired engine control strategy was applied to mimic the transient nature of engine cold-start operation. The results indicate that, under the same total number of injections, increasing the average in-cylinder equivalence ratio through longer injection durations improves combustion stability, as indicated by the decrease of the coefficient of variation of IMEP n (Net Indicated Mean Effective Pressure) from 6.0 % to 0.6 %. However, this strategy leads to higher soot emissions, which increased by nearly an order of magnitude, primarily due to enhanced wall-film formation. In contrast, increasing the number of injections while maintaining a constant equivalence ratio significantly impacts fuel-film deposition and, consequently, soot emissions, with a fivefold reduction of the measured engine-out soot, decreasing from 3.5 mg to 0.7 mg. A soot-film correlation was developed and achieved a high coefficient of determination (r 2 = 0.95) and was further extended to account for spark timing effects. These findings confirm the effectiveness of split-injection for avoiding wall film formation and soot emissions, and the critical role of fuel film in soot generation, supporting the hypothesis that pool fires play a crucial role in contributing to soot formation under these cold-start conditions. In conclusion, the study also indicates the value of a predictive soot-film correlation for developing cold-start emission control strategies.

Engine Cold-start↗

DOE OSTI2025

Understanding and promoting the reaction kinetics of photothermally driven Diels–Alder reaction with minimized side reactions

The thermoreversible Diels–Alder (DA) reaction involving furan and maleimide precursors has been extensively exploited to develop reversible thermosets for circular economy and additive manufacturing. Recently, applications of photothermal nanoparticles, which can absorb light and generate nanoscale heat, have gained significant attention in driving the Diels–Alder (DA) reaction more efficiently. This approach can be utilized by various applications, for example, solar light-enhanced recycling, self-healing at targeted area, and photothermal 3-D printing. Here, in this work, we address an important fundamental question on the photothermal approach: how does the reaction kinetics of the photothermally driven DA reaction compared to that of the conventional heat driven one? We found that the forward DA reaction kinetics varies dramatically depending on a type of heat sources, i.e., light-induced photothermal heat or conventional heat, despite of a similar bulk temperature. We observed that the light-induced photothermal heat significantly enhances the DA reaction rate in comparison to the bulk heating in an inert environment. On the contrary, the photothermally driven forward DA reaction is significantly slower than that of the heat driven reaction in the presence of oxygen. This is attributed to the photoexcited electrons in photothermal nanoparticles that can generate singlet oxygen ( 1 O 2 ) in the presence of visible light. This singlet oxygen may facilitate an endoperoxide side reaction on the furan group, ultimately, retard the DA reaction from the furan group. These results reveal the effectiveness of photothermal nanoparticles in accelerating DA reaction while underscoring the need to minimize oxygen exposure to prevent unfavorable side reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Magnetosynthesis Effect on the Structure and Ground State of Cu 2+ -Based Antiferromagnets

Synthetic variables can have an outsized influence on the crystal structure and magnetic properties of a material, particularly those of quantum materials. In this work, we investigate the impact of synthesis under a magnetic field (magnetosynthesis) on the crystal structure and magnetic properties of several Cu 2+ (S = 1/2)-based materials with antiferromagnetic interactions and varying levels of magnetic frustration, from simple antiferromagnets to a quantum spin liquid. Here, we develop methods to apply small (0.09–0.37 T) magnetic fields during low-temperature hydrothermal, evaporative, and rehydration syntheses of the simple antiferromagnet CuCl 2 ·2H 2 O, the canted antiferromagnet (Cu,Zn) 3 Cl 4 (OH) 2 ·2H 2 O, the frustrated and canted antiferromagnet atacamite Cu 2 (OH) 3 Cl, and the highly frustrated quantum spin liquid herbertsmithite Cu 3 Zn(OH) 6 Cl 2 . We report the first single-crystal X-ray structural determination of the Cu 3 Cl 4 (OH) 2 ·2H 2 O structure type and probe the stability of this phase both experimentally and computationally. Atacamite Cu 2 (OH) 3 Cl synthesized under a 0.19 T field experiences a 0.15 K (∼3%) decrease in its Néel transition temperature. This result suggests that magnetosynthesis with small applied fields may have a very subtle influence upon the magnetic properties of moderately magnetically frustrated 3 d materials.

36 MATERIALS SCIENCE↗

DOE OSTI2026

In Situ Insights into Enhanced Cooperative Ligand Exchange Kinetics via Solvent-Induced Restacking in a 2D Metal–Organic Framework

Understanding the reaction kinetics at catalytically active sites is crucial for integrating catalytic two-dimensional (2D) materials into industrial processes. This study focuses on in situ observation of ligand exchange kinetics and solvent-assisted structural restacking transition in the 2D paddle wheel-based MOF [Cu 2 (dttc) 2 ] n (DUT-134(Cu), dttc = dithieno[3,2-b:2′,3′- d]thiophene-2,6-dicarboxylate). The ligand exchange process, involving the replacement of dimethylformamide (DMF) with nitriles such as acetonitrile (ACN), pentanenitrile, and heptanenitrile, was investigated using advanced in situ characterization techniques with high temporal resolution, including powder X-ray diffraction and Raman spectroscopy. The larger analytes exhibited reduced exchange rates, consistent with enhanced steric hindrance and greater diffusion constraints. Interestingly, the study revealed that the exchange of DMF with ACN induces a structural transition to higher symmetry within few seconds, a transition from AB to AA stacking mode of the layers, and a widening of the interlayer distance. Crucially, this structural transition dramatically accelerates the solvent exchange process through cooperative effects, offering critical advantages for catalytic applications. Notably, the reverse exchange from ACN to DMF proceeds more slowly and does not reverse the structural changes, but a new phase is formed with preserved AA stacking. By isotope labeling of linker molecules in combination with two complementary theoretical vibrational simulation methods, the precise assignment of Raman bands and the vibrational modes associated with the ligand exchange process could be achieved. These pioneering insights into the dynamic behavior of 2D MOFs, coupled with ligand exchange, establish a highly promising and transformative approach to achieving enhanced tunability and responsiveness in future catalytic applications.

Layers↗

DOE OSTI2025

Terrestrial atmospheric ion implantation occurred in the nearside lunar regolith during the history of Earth’s dynamo

Light volatile elements in lunar regolith are thought to be a mixture of the solar wind and Earth’s atmosphere, the latter sourced in the absence of geomagnetic field. However, the extent to which both the current and primitive geodynamo influence the transport of terres trial ions still remains unclear, and this uncertainty is further complicated by the enigmatic composition and poorly constrained location of the Eoarchean exosphere. Here we use 3-D MHD numerical simulations with contemporary magnetized and Archean unmagnetized atmospheres to investigate how Earth’s intrinsic magnetic field affects this transfer, aiming to constrain how and when the lunar isotopic signature provides a record of Earth’s pale-oatmosphere. We find that atmospheric transfer is efficient only when the Moon is within Earth’s magnetotail. The non-solar contribution to the lunar soil is best explained by implantation during the long history of the geodynamo with present-day solar wind conditions, rather than any short, putatively unmagnetized epoch of early Archean Earth. This further suggests the history of the terrestrial atmosphere, spanning billions of years, could be preserved in buried lunar soils. Our results indicate that the elemental abundances of Apollo samples are very sensitive to Earth’s hydrodynamic exobase altitude, which, at the time of ion implantation, was never smaller than 190 km.

Atmospheric dynamics↗

DOE OSTI2026

Activation of methane by U + studied by guided ion beam tandem mass spectrometry and quantum chemistry

Reaction pathways of all products formed in the U + + CH 4 (CD 4 ) reaction were explored as a function of kinetic energy using guided ion beam tandem mass spectrometry and quantum chemical calculations. UH + , UC + , UCH + , UCH 2 + , and UCH 3 + (and their perdeuterated analogues) are formed in endothermic reactions. In both systems, the UCH 2 + (UCD 2 + ) dehydrogenated product was the dominant product in the low-energy region, whereas the UH + (UD + ) hydride product became predominant at high energies. The kinetic energy behavior of the various products is consistent with a common intermediate of H–U + –CH 3 (D–U + –CD 3 ). Here, the kinetic energy dependence of all product cross sections was modeled to obtain experimental bond dissociation energies at 0 K (in eV): D 0 (U + –H) = 2.42 ± 0.10, D 0 (U + –C) = 3.95 ± 0.12, D 0 (U + –CH) = 4.91 ± 0.09, D 0 (U + –CH 2 ) = 4.11 ± 0.04, and D 0 (U + –CH 3 ) = 2.41 ± 0.09. Quantum chemical calculations using the UCCSD(T) and UB3LYP approaches with the cc-pwCVXZ-PP basis set with MDF-60 pseudopotential for U + and the aug-cc-pCVXZ and aug-cc-pVXZ (X = T, Q) basis set for carbon and hydrogen, respectively, validate the experimental bond dissociation energies and outline the potential energy surface for all reactions observed. In addition, spin–orbit corrections of the bond energies for all products were calculated at a CASSCF-CASPT2-RASSI level.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Data-informed grid refinement to improve traveltime accuracy in the regional seismic traveltime (RSTT) model

The regional seismic traveltime (RSTT) model predicts traveltimes of regional seismic phases accounting for 3-D structure of the crust and the upper mantle on a global scale. Previous versions of the RSTT model have been implemented using nodes separated by ∼1° spacing across the globe. A regional-scale study using regional Pn and Pg traveltimes across Israel and the Middle East demonstrated that data driven, systematic grid refinement reduces traveltime residuals and enhances resolution of smaller tectonic features in regions having dense ray coverage. High density Pn ray coverage in the western US, Europe, Middle East and East Asia can likewise provide the resolution that allows systematic global grid refinement of the RSTT model. In this study, we use a large number of Pn ray paths originating from events located with an epicentral location uncertainty of 25 km (GT25) or better. We conduct targeted grid refinements at 1.0°, 0.5°, 0.25° and 0.125° on a global scale, producing a refined RSTT model that yields a 21.6 per cent reduction in median event location error in Europe and the Middle East, when compared with the original global RSTT model presented in Begnaud et al. The new model also resolves finer tectonic structures in regions with high Pn ray density.

58 GEOSCIENCES↗

DOE OSTI2026

Stellarator Design Exploration Using Symbolic-Regression Neutronics Surrogates

Systems codes require fast, simplified models to rapidly evaluate fusion power plant concepts, but neutronics analyses are often a computational bottleneck. Here, to address this, surrogate models for key neutronics responses have been developed using 3-D neutronics-ready models built with the open-source code ParaStell from a database of stellarator equilibria. Neutronics responses such as tritium breeding ratio (TBR), nuclear heating, and neutron-induced radiation damage displacements per atom (dpa) were simulated using OpenMC. Through sensitivity analysis and symbolic regression (SR), simple power-law formulas were derived connecting these neutronics responses to global stellarator parameters, including fusion power, plasma surface area, and plasma elongation. Validation shows these formulas can predict the simulation results with low error, enabling quick and accurate assessment of neutronics requirements in stellarator design exploration activities with systems codes.

Modeling↗

DOE OSTI2026

Multiphysics Analysis of Li Cooled Divertor Substrate During Loss of Coolant Accident (LOCA)

In the ongoing study of potential designs for liquid metal (LM) plasma-facing components (PFCs), so-called “slow” and “fast” Li flow divertor concepts are under investigation. In the previous studies on design and analysis of the slow Li flow divertor and comparison with the fast Li flow divertor, the magnetohydrodynamics (MHD)/heat transfer effects of the Li flowing inside the substrate as a second coolant were comprehensively investigated under the normal steady-state operation conditions. Here, in the present study, the multiphysics analysis is extended to the unsteady abnormal divertor scenario where the Li layer on top of the substrate does not provide full coverage or even totally disappears for a certain period of time, so that the substrate becomes directly exposed to the incident high plasma heat flux. Such an unwanted event may happen regardless of the concept of the divertor and is worth detailed investigations, typically referred to as a loss of coolant accident (LOCA). To address this situation, a simplified scoping analysis is conducted first in 2-D, and then an integrated 3-D modeling is performed using a time-dependent multiphysics model in COMSOL Multiphysics that integrates LM MHD, heat transfer, and solid mechanics. The main goal is to evaluate conditions under which the major material limits, such as the maximum allowable temperature, stress, and displacement of the substrate, can still be met. It was shown that the maximum time over which the substrate of RAFM steel can retain structural integrity during the LOCA is around 0.2 ∼ 0.3 s. Any divertor concept that utilizes RAFM steel as a substrate material and liquid Li as a second coolant should take such a permitted time into consideration.

divertor↗

DOE OSTI2025

The Sensitivity of Nickel Substitution in the Structural Design of Manganese-Based Layered and Lithiated Spinel Electrodes for Li-Ion Batteries

A concerted effort is being made at Argonne National Laboratory to explore and develop high-performance lithium-manganese-nickel-oxide cathode materials for the lithium battery industry. They are both energy- and cost-competitive relative to cobalt- and nickel-rich systems, such as LiCoO 2 , LiNi 0.8 Co 0.15 Al 0.05 O 2 , and LiNi 0.8 Mn 0.1 Co 0.1 O 2 . The recent discovery of the structural versatility of LiNi 0.5 Mn 0.5 O 2 , has revealed that it can exist in three configurations: lithiated spinel, partially disordered layered, and completely disordered rock salt within a common cubic-close-packed oxygen array. This versatility opens up the possibility of designing structurally integrated electrode components with interlinking 2-D (layered) and 3-D (spinel) channels for Li-ion transport. This paper highlights the sensitivity of Ni substitution and synthesis temperature on the structural and electrochemical properties of LiNi 0.50 Mn 0.50 O 2 , LiNi 0.53 Mn 0.47 O 2 , LiNi 0.47 Mn 0.53 O 2 , and LiNi 0.25 Mn 0.75 O 2 electrodes, as well as a cobalt-containing composition, LiNi 0.333 Mn 0.333 Co 0.333 O 2 , when synthesized at a relatively low temperature of 400 °C.

25 ENERGY STORAGE↗

DOE OSTI2026

HELIX

HELIX — Hybrid Envelope-multiparticle LInac eXplorer: TraceWin-compatible linac simulation with 3-D PIC space charge, matching, a PyQt6 workbench, and a voice-driven AI copilot. Built at Fermilab for PIP-II.

Pathak, Abhishek [Fermi National Accelerator Labor↗

DOE OSTI2026

Improving High Resolution Offshore Wind Resource Assessments and Forecasts using Observations in the MA/RI Lease Areas

The third Wind Forecast Improvement Project (WFIP3) sought to improve understanding of the physical phenomena in the atmosphere and ocean that dictate the structure and variability of wind and thermodynamic fields within the Marine Atmospheric Boundary Layer (MABL). WFIP3 conducted a comprehensive 18-month observational study over the Northeast U.S. outer continental shelf, a high use coastal zone, using a 3-D multiscale sensor array to highly resolve the temporal, vertical, and horizontal structure of the coupled atmospheric and oceanic boundary layers. Multiple land-based study sites adjacent to the coastal ocean observed surface meteorology and vertical profiles of atmospheric properties via passive infrared and microwave radiometers, active lidars and radars, and radiosondes. At sea, an array of surface flux buoys and two vertical profiling lidar buoys observed both atmospheric and oceanic properties, augmented by land-based oceanographic radar systems and routine ship-based surveys. Intensive observations of the MABL over the ocean was done from an air-sea interaction flux tower and extended deployments of a large barge platform. WFIP3 focused on mesoscale and sub-mesoscale flows -- including sea breezes, low-level jets, low-level clouds, and coastal storms -- and the ability of advanced numerical model parameterizations to represent them within fully coupled oceanic and atmospheric modeling systems and foundational weather forecast models. Numerous critical forecasting phenomena were observed, however the project was terminated prior to the completion of the field observational period and the analysis period.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Numerical Investigation of Fluid Flow and Space Charge in Liquid Argon Time Projection Chamber (LArTPC) Detectors

Overview This project focused on developing a high-fidelity numerical framework to simulate the multiphysics environment within Liquid Argon Time Projection Chamber (LArTPC) detectors. The primary objective was to characterize the complex interplay between ion transport, background fluid dynamics, and electric field distortions—a critical factor for the calibration and sensitivity of next-generation High Energy Physics experiments, such as DUNE. Technical Achievements The research successfully yielded a hybrid numerical space-charge solver utilizing a Cell-Centered Finite Volume Method (FVM) for ion transport coupled with a Finite Element Method (FEM) for electric potential. Key accomplishments include: • Verification & Validation: The 3-D solver was rigorously verified against 1-D analytical solutions, demonstrating high numerical accuracy in predicting space-charge-induced field deviations. • Field Distortion Analysis: 3D simulations revealed that space charge effects introduce significant non-uniformities in the electric field. Critically, the research identified that background LAr flow velocities, when comparable to ion drift velocities, markedly exacerbate these distortions. • Technology Transfer: The resulting source code and comprehensive user manuals were successfully transferred to collaborators at Fermilab, providing a portable computational tool for the broader scientific community. Challenges and Future Directions While the space-charge solver achieved all performance metrics, the integrated fluid dynamics modeling encountered convergence challenges stemming from the extreme 200-fold disparity in length scales between the detector's 37 mm inlet pipes and the 8-meter global domain. To address this, the project has identified a clear technical pivot toward Hierarchical Geometric Adaptive Mesh Refinement (HG-AMR). By implementing an h-type refinement strategy with hanging nodes, future iterations of this solver will be capable of resolving localized high-gradient inlet flows without the prohibitive computational costs of regular grids. This advancement, combined with data-driven uncertainty quantification based on MicroBooNE-style calibration, will enable the precise modeling of detector responses in large-scale cryogenic environments where direct measurement remains difficult. Impact The computational tools developed under this award provide a foundation for enhancing the energy resolution and spatial reconstruction of noble liquid detectors. By bridging the gap between theoretical fluid dynamics and experimental field calibration, this work supports the DOE’s mission to advance the frontiers of neutrino physics and dark matter detection.

42 ENGINEERING↗

DOE OSTI2025

Water Column Respiration in the Yakima River Basin is Explained by Temperature, Nutrients and Suspended Solids

Understanding aquatic ecosystem metabolism involves the study of two key processes: carbon fixation via primary production and organic C mineralization as total ecosystem respiration (ERtot). In streams and rivers, ERtot includes respiration in the water column (ERwc) and in the sediments (ERsed). While literature surveys suggest that ERsed is often a dominant contributor to ERtot, recent studies indicate that the relative influence of sediment-associated processes versus water column processes can fluctuate along the river continuum. Still, a comprehensive understanding of the factors contributing to these shifts within basins and across stream orders is needed. Here we contribute to this need by measuring ERwc and collecting water samples across 47 sites in the Yakima River basin, Washington, USA. We found that ERwc rates varied throughout the basin during baseflow conditions, ranging from –7.38 to 0.36 g O2 m?3 d?1, and encompassed the range of ERwc literature values. Additionally, by comparing to ERtot estimates for rivers across the contiguous United States, we suggest that the contribution of ERwc rates to reach-scale ERtot rates across the Yakima River was likely highly variable, but we did not test this directly. We observed that temperature, nutrient concentrations (dissolved organic carbon, total dissolved nitrogen), and total suspended solids explained 41% of ERwc variability across the basin. Our findings highlight the potential relevance of water column processes in aquatic ecosystem metabolism, with the Yakima River basin serving as an environmentally diverse river network representative of the larger Columbia River basin that spans much of the Pacific Northwest region of the United States. Our results are generally congruent with previous work, suggesting that the observed variability and suite of associated environmental factors influencing ERwc are potentially transferable across basins.

Laan, Maggi M.↗

NASA NTRS1991

Performance capabilities of a JPL dual-arm advanced teleoperation system

The system comprises: (1) two PUMA 560 robot arms, each equipped with the latest JPL developed smart hands which contain 3-D force/moment and grasp force sensors; (2) two general purpose force reflecting hand controllers; (3) a NS32016 microprocessors based distributed computing system together with JPL developed universal motor controllers; (4) graphics display of sensor data; (5) capabilities for time delay experiments; and (6) automatic data recording capabilities. Several different types of control modes are implemented on this system using different feedback control techniques. Some of the control modes and the related feedback control techniques are described, and the achievable control performance for tracking position and force trajectories are reported. The interaction between position and force trajectory tracking is illustrated. The best performance is obtained by using a novel, task space error feedback technique.

Szakaly, Z. F.↗

NASA NTRS1991

Preparation of the Ice Certification of the Dornier 328 Regional Airliner By Numerical Simulation and By Ground Test

The Dornier 328, a new regional airliner, is to be qualified to FAR25/JAR25 requirements for operation into known icing conditions. All lifting surfaces are pneumatically deiced. Impingement limits were determined using a validated computational fluid dynamic (CFD) code. Wind tunnel tests with a model of the whole aircraft and the isolated empennage equipped with artificial ice shapes showed a degradation of handling characteristics due to ice accretion. Extensive two and three dimensional CFD calculations evaluated ice accretion at the deflection elevator horn. Icing tunnel tests were performed with a model of the horizontal tail with icing similitude fulfilled. The Dornier Do328 powerplant air induction system has successfully completed a program of icing tests in support of aircraft certification for operation into known icing conditions. Potential flight icing conditions were extensively analyzed using a CFD potential flow code and 3-D particle tracking routine to identify critical flight icing conditions. The tunnel test program verified the analytically predicted critical ice accretion surfaces of the air induction system, the adequacy of the ice protection provisions and demonstrated compliance with the applicable (JAR25) regulatory material.

D Welte↗

DOE OSTI2026

A New Window into the Baryon Cycle at Cosmic Noon with Line Intensity Mapping: Forecasts for auto- and cross-correlations in [CII]-158$μ$m, HI 21 cm, CO$_{J+1\rightarrow J}$, and H$α$ galaxies

Across the peak of cosmic star formation at $z\sim1-2$, inflow, processing, and feedback drive rapid changes in the spatial distribution and chemical composition of baryons in galaxies and surrounding reservoirs; this baryon cycle can be tomographically mapped by line intensity mapping (LIM) of atomic hydrogen, ionized carbon, and carbon monoxide. We present a simulation-based forecasting framework for detecting auto- and cross-power spectra between spectroscopic surveys of four such tracers at $z\sim0.5-1.7$ mapping the same deep field - TIM, EoRSpec/FYST, MeerKAT, & Euclid. We forward-model 3-D distributions for these tracers from magnetohydrodynamic simulations, directly capturing the two-halo, one-halo, and shot statistics without relying on analytical decompositions. We further detail a signal-to-noise formalism, tailored to LIM surveys with highly anisotropic geometries and Fourier-space coverage. We demonstrate that galaxy cross-correlations will be the dominant discovery channel for current-generation surveys. These instruments will detect the auto-spectra for CO and HI 21 cm and the CO $\times$ 21 cm cross-spectrum at modest S/N $\sim 1-10$, while placing upper limits on the [CII]-158$μ$m signals. [CII], CO, and HI LIM will be $\sim3-30\times$ ($0.5-1.5$ dex) more sensitive to cross-correlation with the Euclid survey, however, than their respective auto-correlations, constraining all three models of line emission at high significance (S/N $\sim 10-40$) within this decade. Finally, we formulate a staged instrumental trajectory with planned or reasonable improvements, including the as-proposed SKA-Mid. We forecast advancing the per-$k$-mode sensitivities of each auto-, galaxy-line, and line-line spectrum by several orders of magnitude, enabling new percent- and sub-percent level constraints on cosmology and the redshift evolution of star formation and the baryon cycle.

Agrawal, Shubh [Pennsylvania U., Dept. Math.]↗

DOE OSTI2026

Combined effects of horizontal and vertical resolution on reliable turbulence prediction at tidal energy sites: A systematic study in the Salish Sea, WA

Predicting turbulence characteristics with coastal ocean models is essential for tidal energy converter deployment. While large eddy simulation provides a detailed turbulence representation, computational limitations restrict its use to smaller domains. We systematically evaluate whether well-configured coastal models can provide reliable turbulence prediction through progressive refinement of 3D model representation. We implemented four model configurations (Levels 1–4) using terrain-following coordinates, isolating the impacts of horizontal resolution, vertical resolution, and layer distribution. We validated all configurations against field measurements from the Salish Sea, WA. Results show that tidal current velocity predictions remain unchanged regarding model configuration, but turbulence properties are sensitive to resolution refinement. Increasing vertical resolution alone proved insufficient; even with vertical sigma-levels rising from 11 to 41, significant underprediction persisted until finer horizontal resolution better captured bathymetric variations. The Level 4 configuration, incorporating geometric sigma-levels distribution, achieved turbulence prediction skill scores exceeding 0.90. Turbulence closure comparison revealed Mellor- Yamada 2.5 outperformed k-epsilon in TKE prediction (skill scores 0.84–0.94 versus 0.72–0.81) due to better boundary layer parameterization. This study shows that well-configured coastal models effectively bridge the gap between simplified tools and costly high-fidelity modeling, offering the tidal energy industry practical and cost-effective turbulence data at commercially relevant scales.

Marine Energy↗