SEARCH · Search NASA

Results for “explainability”

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 361 records · Page 20

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Iron Oxide Reduction Rate Affects Iron Product Morphology

Magnetite concentrates from different parts of the Mesabi Range, exhibiting varying tendencies to encapsulate, were reduced using thermogravimetric analysis (TGA) and a high-temperature confocal scanning laser microscope (CSLM). Gas-phase mass-transfer conditions were intentionally varied by employing two reactors with different geometries: the TGA produced significantly slower reduction, whereas the CSLM setup yielded reduction rates up to 20 times faster under mass-transfer-controlled conditions. A preliminary mechanism is proposed to explain the transition from a porous iron product to a dense encapsulating layer. Encapsulation is interpreted as the outcome of competition between pore creation, driven by the reduction rate, and pore elimination, driven by surface diffusion. When pore formation is insufficient to counteract pore coarsening, a dense iron layer develops and restricts further reduction.

encapsulation↗

DOE OSTI2024

Physical, socio-psychological, and behavioural determinants of household energy consumption in the UK

Determining which attitudes and behaviours predict household energy consumption can help accelerate the low-carbon energy transition. Conventional approaches in this domain are limited, often relying on survey methods that produce data on individuals’ motivations and self-reported activities without pairing these with actual energy consumption records, which are particularly hard to collect for large, nationally representative samples. This challenge precludes the development of empirical evidence on which attitudes and behaviours influence patterns of energy consumption, thus limiting the extent to which these can inform energy interventions or conservation programs. This study demonstrates a novel methodology for estimating energy consumption in the absence of actual energy records by using a large, publicly available data set of energy consumption in the UK. We develop a predictive model using the Smart Energy Research Laboratory (SERL) data portal (with records from nearly 13,000 UK households) and then use this model to predict energy consumption (both electric and gas) for a sample of 1,000 UK householders for which we separately collect over 200 variables relating to climate change attitudes and practices. Our approach uses a set of over 50 independent variables that are shared between the data sets, allowing us to train a model on the SERL data and use it to analyse the relationship between energy consumption and the opinions, motivations, and daily practices of survey respondents. Results show that electricity consumption is influenced by a broader range of factors compared to gas. Household energy use is best explained by physical dwelling characteristics, socio-demographic variables, and certain behavioural and attitudinal measures. Notably, pro-environmental attitudes, frugality, and conscientiousness correlate with lower energy use, while income and consumerism are linked to higher consumption. We discuss how these findings can inform efforts to decarbonise home energy use in the UK.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2024

Eco-driving Profile Optimization by Dynamic Programming for Battery Electric Vehicles

Although full automation has not yet been achieved, automated vehicles are a valid research area. Not only would automated vehicles provide ultimate driver convenience, but they would maximize energy efficiency by eliminating undesired human driving behaviors and optimally controlling the powertrain. From the perspective of control related to energy saving, speed profile optimization is important for improving system efficiency and satisfying passenger demands. This study employs Dynamic Programming (DP) to solve the constrained optimal problem for travel time, distance, and speed limit by exploring all possible control options. The solutions obtained by DP demonstrate consistent control patterns combining four control modes-acceleration, cruising, coasting, and braking, with cruising or coasting being selective depending on the boundary conditions. Further, this study introduces DP-based simulation results and attempts to provide comprehensive interpretations of the optimal policy by analyzing the essential factors that affect the control problem, including boundary conditions, road load, and powertrain characteristics. Based on these interpretations, the control concepts can be explained as the optimal policy selecting the best control option based on system efficiency and boundary conditions. The results of DP are compared with a human-like driver model to show that the optimal speed profiles can effectively reduce energy consumption.

Autonomous vehicles↗

DOE OSTI2025

Magnet-superconductor hybrid quantum systems: a materials platform for topological superconductivity

Magnet–superconductor hybrid (MSH) systems have recently emerged as one of the most significant developments in condensed matter physics. This has generated, in the last decade, a steadily rising interest in the understanding of their unique properties. They have been proposed as one of the most promising platforms for the establishment of topological superconductivity, which holds high potential for application in future quantum information technologies. Their emergent electronic properties stem from the exchange interaction between the magnetic moments and the superconducting condensate. Given the atomic-level origin of such interaction, it is of paramount importance to investigate new magnet–superconductor hybrids at the atomic scale. In this regard, scanning tunneling microscopy (STM) and spectroscopy are playing a crucial role in the race to unveil the fundamental origin of the unique properties of MSH systems, with the aim to discover new hybrid quantum materials capable of hosting topologically non-trivial unconventional superconducting phases. In particular, the combination of STM studies with tight-binding model calculations have represented, so far, the most successful approach to unveil and explain the emergent electronic properties of MSHs. The scope of this review is to offer a broad perspective on the field of MSHs from an atomic-level investigation point-of-view. The focus is on discussing the link between the magnetic ground state hosted by the hybrid system and the corresponding emergent superconducting phase. This is done for MSHs with both one-dimensional (atomic chains) and two-dimensional (atomic lattices and thin films) magnetic systems proximitized to conventional s-wave superconductors. We present a systematic categorization of the experimentally investigated systems with respect to defined experimentally accessible criteria to verify or falsify the presence of topological superconductivity and Majorana edge modes. The discussion will start with an introduction to the physics of Yu–Shiba–Rusinov bound states at magnetic impurities on superconducting surfaces. This will be used as a base for the discussion of magnetic atomic chains on superconductors, distinguishing between ferromagnetic, antiferromagnetic and non-collinear magnetic ground states. A similar approach will be used for the discussion of magnetic thin film islands on superconductors. Given the vast number of publications on the topic, we limit ourselves to discuss works which are most relevant to the search for topological superconductivity.

Majorana states↗

DOE OSTI2025

Simulated annealing of reduced magnetohydrodynamic systems

Theory of simulated annealing (SA), a method for equilibrium and stability analyses for Hamiltonian systems, is reviewed. The SA explained in this review is based on a double bracket formulation that derives from Hamiltonian structure. In addition to general theoretical aspects, the explicit formulation as well as numerical applications are presented. Both finite and infinite degree-of-freedom systems are treated, in particular, the heavy top, a toy model mimicking low-beta reduced magnetohydrodynamics (MHD) and low- and high-beta reduced MHD. Furthermore, the numerical results successfully demonstrate the usefulness of SA for equilibrium and stability analyses. At the same time, the results raise some future issues that are discussed in the paper.

Poisson Bracket↗

DOE OSTI2026

Water availability modulates maximum canopy heights of low-elevation Amazonian second-growth forests

Tropical second-growth forests of the Amazon sequester large amounts of carbon and are important carbon sinks, contributing substantially to climate change mitigation, biodiversity conservation, and providing crucial ecosystem services. Deforestation due to selective logging and shifting cultivation is expanding second-growth forest areas in tropical forest regions, which if well managed, regenerate rapidly over time. Maximum forest canopy height is an important metric of biomass and carbon accumulation in second-growth forests and is strongly influenced by water availability. The water limitation hypothesis explains the positive influence of water availability on maximum tree heights and has been examined and demonstrated at a small-scale using field data, and at a global scale, with limited accuracy, using remote sensing data in tropical ecosystems. However, this hypothesis concerning maximum canopy height has not been much studied at regional and national scales for tropical second-growth forests. In this study, we leveraged NASA GEDI spaceborne lidar data across the Brazilian Amazon and derived second-growth forest relative height metrics for delineating the influence of water availability, second-growth forest age, and topographic elevation on maximum canopy height. Water availability was found to significantly influence the maximum canopy height of second-growth forest trees, of age range from 30 to 35 years, at elevations less than 500 m and maximum precipitation thresholds of 1500 mm. Our results indicate that changing precipitation patterns or increased drought conditions under different climate change regimes could impact forest structure, plant communities, ecosystem functioning, and carbon sequestration capabilities of tropical second-growth forests in the Amazon.

Carbon sequestration↗

DOE OSTI2024

The many faces of θ' -Al 2 Cu precipitates: Energetics of pristine and solute segregated Al/ θ' semi-coherent interfaces

θ'-Al 2 Cu precipitates in Al-Cu alloys have various distorted octagon shapes, which can be explained by the competition between {100} and {110} type semi-coherent interfaces with the Al matrix. While most prior studies on the semi-coherent Al/θ' interfaces have focused on the {100} orientation, little is known about the {110} interface. We have investigated the energetics of pristine and solute-segregated {110} semi-coherent Al/θ' interfaces with advanced characterization and first-principles studies. We report interfacial, strain, and solute segregation energetics of the {110} Al/θ' semi-coherent interface for 39 elements and compared them with previously reported values of the {100} interface. We discuss the atomic features and atomic local structures to identify similarities and differences between the two types of Al/θ' semi-coherent interfaces. Here, the isotropy in pristine Al/θ' semi-coherent interfacial energy and the anisotropy resulting from solute segregation provide insight into the formation of different types of θ' precipitate “faces” reported in the literature.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Rigorous computation of short-range order unifies its controversial effects in complex concentrated alloys

Direct experimental observations of chemical short-range order (SRO) in complex concentrated alloys (CCAs) have triggered high interest. However, the reported effects of SRO on yield stresses are controversial, and their atomic-scale mechanisms are elusive, which limits our ability to utilize SRO in alloy design. Here we tackle this challenge using an advanced computational approach that rigorously takes into account the critical lattice distortion in CCAs and further verify our theoretical predictions with experiments. We show that the CoCrNi model alloy has a narrow temperature window around 670 °C for SRO formation. This explains why the mechanical effect of SRO is observed in some experiments but not in others. Here, we propose an effective alloy-doping method to control SRO and reveal atomic-bonding types that dominate SRO formation for different alloys. The strategies and insights generally apply to a broad spectrum of alloys, laying the foundation for designing advanced alloys by manipulating their SRO.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Assessment of metadynamic recrystallization in single copper particle impacts by focused ion beam tomography

We study single Cu-on-Cu impacts relevant to cold spray deposition and quantitatively analyze the metadynamic recrystallization (mDRX) that takes place after the impact by virtue of lingering impact adiabatic heat. Unlike prior studies, the current full 3D tomographic analysis of the mDRX volume shows that mDRX is extremely common in such impacts, although it is often missed when examining 2D sections. We also report an unexpected trend: there is a “sour spot” for mDRX at velocities about 20–40 % above the velocity for particle adhesion. This non-monotonic trend is contrary to the expectations based on increasing adiabatic heating with velocity. With a schematic model, we show that the trend can be explained on the basis of heat transfer: cooling of the heat-affected region is limited by transport through the bonded regions at the particle-substrate interface. Thus, bonding has a prominent role in the heat dissipation process and the best bonded particles most rapidly bulk quench, avoiding mDRX. Here, the developed semi-empirical model aligns with the experimental findings and may help inform microstructural evolution during cold spray and post-spray processing.

FIB-SEM tomography↗

DOE OSTI2025

Chemical trends favoring interstitial cluster formation in bcc high-entropy alloys from first-principles calculations

Achieving high strength and ductility is a common goal in structural alloy design. Body-centered cubic high-entropy alloys (HEAs) commonly highlight the conflict between these properties, with stronger alloys being brittle and vice versa. Recent reports suggest interstitial solutes can be used to overcome this trade-off, in some cases providing both strength and ductility enhancements. This effect has been correlated with interstitial cluster formation, although the conditions favoring their formation remain incompletely understood. Using first-principles calculations of solution energies and diffusivities, we provide insights into thermodynamic and kinetic factors favoring interstitial solute clusters. Among C, N and O solutes, O interstitials display most desirable diffusion kinetics. Further, the results highlight the importance of local composition fluctuations in the HEAs to enable the formation of clusters of appreciable size. The results are explained in terms of bonding and distortion trends across solutes and HEA compositions to provide guidelines for alloy design.

Borges, Pedro P P O↗

DOE OSTI2025

Decoding anomalous grain growth at room temperature during pressure-induced phase transformations

Significant grain growth is observed during the high-pressure phase transformations (PTs) at room temperature in various materials. The main focus here is grain growth from a few hundred nanometers to 10 μm within an hour during α → ω PT in Zr. No existing theory explains this phenomenon since without PT, Zr nanocrystals do not grow at room temperature even for up to 10 years. Here, in this study, a multistep mechanism for the grain growth during α → ω PT in Zr is suggested. Phase interfaces (PI) and grain boundaries (GBs) coincide and move together as a combined PI-GBs under the action of the combined thermodynamic driving force. Such a combined motion changes the diffusional grain growth mechanism to the transformational one and the martensitic mechanism of PT to a reconstructive one via an intermediate disordered phase. The primary condition is that the GB energy of the ω phase is smaller than that of the α phase, which promotes the nucleation of ω-Zr and is consistent with the absence of the reverse PT and reduction in the PT pressure with reducing grain size. Several intermediate steps for such motion are suggested and justified kinetically. Nonhydrostatic stresses due to volume reduction in the growing ω grain promote continuous growth of the existing ω grain instead of a new nucleation at other GBs. In situ synchrotron Laue diffraction experiments confirm the main predictions of the theory. The suggested mechanism provides a new insight into synergistic interaction between PTs and microstructure evolution.

anomalous grain growth during phase transformation↗

DOE OSTI2026

Formation of carbon homonuclear bonds in β-SiC under neutron irradiation at various temperatures and neutron doses

To elucidate radiation defect processes in SiC, Raman spectroscopy was systematically applied to high-purity, polycrystalline β-SiC that was neutron irradiated at a range of temperature and dose conditions. The analysis specifically focused on formation of carbon homonuclear bonds by irradiation; these bonds were indicated by D and G bands and amorphous carbon peaks. Intensity of the carbon peaks relative to SiC peaks significantly decreased in the case of high temperature and/or high neutron dose of 500 °C to 29 displacements per atom (dpa) and about 800 °C to 1.38 and 29 dpa. The absence of carbon bond peaks under those conditions was explained by growth of stoichiometric defect clusters, consistent with previous atomistic simulations on SiC defect stability. The lack of Raman bands associated with carbon clusters under high-temperature and high-dose radiation conditions accounts for the resistance of SiC to phase separation under irradiation. The findings further suggest that material compositions and chemical properties that are inherently resistant to chemical disordering under high-dose radiation conditions are indicative of the long-term durability of ceramic compounds in radiation environments.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Understanding the influence of boron in additively manufactured GammaPrint®-700 CoNi-based superalloy

Boron is commonly added to superalloys in small amounts to enhance creep resistance, but can lead to cracking at high concentrations, especially during the additive manufacturing process. Two variants of CoNi-based GammaPrint®-700 superalloy with different B contents (0.08 at% vs 0.16 at%) were printed via laser powder bed fusion (LPBF) with the same printing parameters, with only the high B alloy exhibiting solidification cracking. Atom probe tomography (APT) revealed stronger segregation behaviors in the high B alloy compared to the low B alloy at both the inter-dendritic regions and grain boundaries (GBs). The segregation behavior at inter-dendritic regions was well captured with Scheil simulation and can correlate with the existing cracking susceptibility index (CSI) on cracking tendencies, although high angle GBs are where cracking occurs according to electron backscatter diffraction (EBSD) measurements. Additionally, the extent of GB segregation was compared between the high B and low B alloy. Higher B additions led to significantly more GB B segregation in the high B alloy compared to the low B alloy. Further, for the high B alloy, the cracked region of one GB exhibited higher levels of B compared to the uncracked region of the same GB. However, much higher B contents were also found in two other uncracked GBs in the high B alloy, which demonstrates that higher GB B concentrations are not fully responsible for the cracking. A much larger variance in GB B segregation content was found in the high B alloy compared to the low B alloy. These phenomena were explained with a solidification model with the GB segregation content expressed explicitly by a modified Langmuir-McLean equation. This model linked the GB segregation content with solidification undercooling, which can be used as quantitative cracking criteria for future builds.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Dissolution zone model of the oxide structure in additively manufactured dispersion-strengthened alloys

The structural evolution of oxides in dispersion-strengthened superalloys during laser-powder bed fusion is considered in detail. Alloy chemistry and process parameter effects on oxide structure are assessed through a parameter study on the model alloy Ni-20Cr, doped with varying concentrations of Y 2 O 3 and Al. Small angle neutron scattering measurements of the dispersoid size distribution show the dispersoid size increases with higher laser power, slower scan speed, and increasing Y 2 O 3 and Al content. Complementary electron microscopy measurements reveal reactions between Y 2 O 3 and Al, even in nanoscale dispersoids, and the presence of micron-scale oxide slag inclusions in select specimens. A scaling analysis of mass and momentum transport within the melt pool, presented here, establishes that diffusional structural evolution mechanisms dominate for nanoscale dispersoids, while fluid forces and advection become significant for larger slag inclusions. These findings are developed into a theory of dispersoid structural evolution, integrating quantitative models of diffusional processes – dispersoid dissolution, nucleation, growth, coarsening – with a reduced order model of time-temperature trajectories of fluid parcels within the melt pool. Calculations of the dispersoid size in single-pass melting reveal a zone in the center of the melt track in which the oxide feedstock fully dissolves. Within this zone the final Y 2 O 3 size is independent of feedstock size and determined by nucleation and growth kinetics. If the dissolution zones of adjacent melt tracks overlap sufficiently with each other to dissolve large oxides, formed during printing or present in the powder feedstock, then the dispersoid structure throughout the build volume is homogeneous and matches that from a single pass within the dissolution zone. Gaps between adjacent dissolution zones result in oxide accumulation into larger slag inclusions. Predictions of final dispersoid size and slag formation using this dissolution zone model match the present experimental data and explain process-structure linkages speculated in the open literature.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Hysteresis area at the canopy level during and after a drought event in the Central Amazon

Understanding forest water limitation during droughts within a warming climate is essential for accurate predictions of forest-climate interactions. In hyperdiverse ecosystems like the Amazon forest, the mechanisms shaping hysteresis patterns in transpiration relative to environmental factors are not well understood. From this perspective, we investigated these dynamics by conducting in situ leaf-level measurements throughout and after the 2015 El Niño-Southern Oscillation (ENSO) drought. Our findings indicate a substantial increase in the hysteresis area (H area ) among transpiration (E), vapor pressure deficit (VPD), and stomatal conductance (g s ) at canopy level during the ENSO peak, attributed to both temporal lag and differences in magnitude between g s and VPD peaks. Specifically, the canopy species Pouteria anomala exhibited an increased H area , due to earlier maximum g s rates leading to a greater temporal lag with VPD compared to the post-drought period. Additionally, leaf water potential (ψ L ) and canopy temperature (T canopy ) showed larger H area during the ENSO peak compared to post-drought conditions across all studied species, suggesting that stomatal closure, particularly during the afternoon, acts to minimize water loss and may explain the counterclockwise hysteresis observed between ψ L and T canopy . Here, the pronounced H area during the drought points to a potential imbalance between water supply and demand, underlining the role of stomatal behavior of isohydric species in response to drought.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Refining water and carbon fluxes modeling in terrestrial ecosystems via plant hydraulics integration

Plant hydraulics substantially affects terrestrial water and carbon cycles by modulating water transport and carbon assimilation. Despite improved drought simulations in certain ecosystems through their integration into land surface models (LSMs), the broader application of plant hydraulics in diverse ecosystems and hydroclimates is still underexplored. Here, in this study, we implemented the recently developed Noah-Multiparameterization Land Surface Model (Noah-MP LSM) equipped with a plant hydraulics scheme (Noah-MP-PHS) across 40 FLUXNET sites globally. Employing the Shuffled Complex Evolution-University of Arizona (SCE-UA) auto-calibration algorithm, we optimized key plant hydraulics parameters for these sites spanning eight vegetation types in both arid and humid climates. Noah-MP-PHS significantly improves the simulation of evapotranspiration (ET) and gross primary production (GPP) by better representing atmospheric and soil water stress compared to traditional soil hydraulic schemes (SHSs, such as Noah and CLM). The augmented Noah-MP-PHS models reduce surface flux overestimation and underestimation, exhibiting an average increase of 0.14 and 0.15 in Kling-Gupta Efficiency (KGE) compared to Noah and CLM, respectively. The explicit consideration of plant capacitance in PHS reveals substantial deep-layer and nocturnal root water uptake especially under dry conditions. We employed eXplainable Machine learning (XML) to quantify the model’s relative sensitivity to newly introduced leaf-, stem and root-related parameters in PHS. The sensitivity analysis reveals a rise in root parameter importance and a decline in leaf and stem parameters as conditions shift from humid to arid. These findings indicate that as aridity states vary, the most influential parameters affecting surface fluxes variation may change in parameter calibration for PHS applications. Our findings underscore the importance of incorporating plant hydraulics into LSMs to enhance simulations of terrestrial water and carbon dynamics. These findings are crucial for understanding ecosystem responses to global climate changes and guide the broader application of PHS at larger scales.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Better practices for inferring ecosystem water use strategy from eddy covariance data

Eddy covariance data are critical for inferring ecosystem water use strategies. Yet, such inferences are sensitive to a range of assumptions applied across studies, hindering our understanding of water use strategies within and across eddy covariance sites. A recent analysis across 151 FLUXNET2015 and AmeriFlux-FLUXNET datasets found that poor model performance was the key driver of non-robust inferences of ecosystem water use strategies. Here, we leverage this previous analysis to (i) identify the specific assumptions that improve inference model performance across most sites, (ii) explain the mechanisms behind the performance improvements, and (iii) check whether better performance improves water use inference. We find that the common practice of fitting a model to canopy conductance (G c ) derived from the evapotranspiration (ET) observations, rather than to observed ET itself, artificially amplifies data errors and degrades the model performance. Next, accounting for vegetation dynamics by applying a growing season filter or incorporating satellite LAI data improves performance, but the former practice may remove soil water stress periods. Lastly, using the leaf-to-air vapor pressure deficit (VPD l ) derived from ET observations as a model input may artificially inflate performance. Based on these results, we recommend selecting observed ET (rather than derived G c ) as the response variable, carefully accounting for vegetation dynamics, and avoiding derived VPD l as a model input; these best practices improve model performance by c. 20% and robustness by c. 80% across all eddy covariance sites. Nevertheless, the performance improvements do not always correspond to more robust inference of water use strategies, as model parameter selection and surface energy budget closure corrections still strongly influence the ecosystem water use parameter estimation in a site-specific manner.

AmeriFlux↗

DOE OSTI2024

Numerical simulation of involute-plate research reactor flow behavior using RANS, LES and DNS

This paper investigates the flow behavior of involute-plate research reactors by performing Reynolds-Averaged Navier Stokes simulation (RANS), Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS) of the channel flow between fuel plates. By modeling turbulence with different numerical approaches, this study provides data with three levels of fidelity. For the RANS simulation, three widely used turbulence models, i.e., k-ε, k-ω, Reynolds Stress Turbulence model (RST) are applied by using the commercial CFD code STAR-CCM +. For LES and DNS, the open-source CFD code, Nek5000, is used given its outstanding scalability on High Performance Computer (HPC) and high-order technique. The results from RANS simulations are compared with that from LES and DNS for benchmarking. Both macroscale parameters and turbulence statistics, such as velocity magnitude, lateral velocity and turbulence kinetic energy, are presented and analyzed. The results from RANS simulation achieve good agreement with LES and DNS on velocity and turbulence kinetic energy prediction. The RST turbulence model predicts the most similar flow pattern of lateral velocity as compared to LES and DNS. The Lambda-2 (λ2) criterion with a reasonable threshold is used to demonstrate the instantaneous vortices distribution in the involute channel from both LES and DNS calculation. The DNS simulation captures more detailed turbulence especially near the corner, which explains the discrepancy between LES and DNS results near the corner. The normalized RMS error are defined and calculated to assess the performance of those turbulence models. The RST model captures the anisotropic feature of turbulence, which enable it to outperform other turbulence models for predicting the flow behavior in an involute channel. Although some discrepancies are found between LES and DNS results in the corner, the overall deviations between LES and DNS are found to be small. In conclusion, given that the computational cost of DNS calculation is an order of magnitude higher, using LES data for benchmarking RANS model is a cost-effective approach.

DNS↗