Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CFD Modeling”

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 721 records · Page 40

A Finite Difference informed Random Walk solver for simulating radiation defect evolution in polycrystalline structures with strongly inhomogeneous diffusivity

Diffusivity of species and defects on grain boundaries is usually several orders of magnitude larger than that inside grains. Such strongly inhomogeneous diffusivity requires prohibitively high computational demands for modeling microstructural evolution. Here, this paper presents a highly-efficient numerical solver, combining the Finite Difference method and Random Walk model, designed for accurately modeling strongly inhomogeneous diffusion within polycrystalline structures. The proposed solver, termed Finite Difference informed Random Walk (FDiRW), integrates a customized Finite Difference (cFD) scheme tailored for fast diffusion along thin grain boundaries represented by a single-layer of nodes. Numerical experiments demonstrate that the FDiRW solver achieves an impressive efficiency gain of 1560x compared to traditional Finite Difference methods while maintaining accuracy, making it feasible for personal computer machines to handle diffusional systems with strongly inhomogeneous diffusivity across static polycrystalline microstructures. The model has been successfully applied to simulate radiation defect evolution, showcasing its scalability to engineering scales in both length and time dimensions.

36 MATERIALS SCIENCE↗

Mixing in Low Reynolds Number Reacting Impinging Jets in Crossflow

Previous efforts to model uranyl fluoride formation in an impinging jet gas reactor underpredicted spatial mixing and overpredicted chemical conversion into particulates. The previous fluid dynamics model was based on the solution of the Reynolds Averaged Navier Stokes equations. After simulating fluid dynamics, aerosol dynamics were superimposed onto CFD-simulated gas reactant species concentrations. The current work explores the influence of complex unsteady flow features on the overall flow physics and chemistry for a low Reynolds number, opposed flow, impinging jet gas reactor where there is a low Reynolds number cross flow. The objective of this study was to assess the impact of model formulation on scalar mixing and transport. Here, transient flow simulations were performed using Scale Resolving Simulations. Large-Eddy Simulations with the dynamic Smagorinsky turbulence model were performed along with simulations which directly resolved the flow. Average and root-mean-square (RMS) velocities and species concentrations were computed along with modeled and resolved turbulence kinetic energy (TKE), modeled turbulence dissipation, and modeled turbulent viscosity. Lagrangian flow tracers were also used to quantify species concentrations along path lines emanating from the jet tips. Transient simulation data were compared to results from RANS simulations using the k-ω shear stress transport (SST) model and Reynolds Stress Model (RSM). Transient simulations showed spatial mixing patterns which were more consistent with experimental data and helped elucidate the process of particle formation observed in experiments.

42 ENGINEERING↗

Influence of Blood Rheology and Turbulence Models in the Numerical Simulation of Aneurysms

An aneurysm is a vascular malformation that can be classified according to its location (cerebral, aortic) or shape (saccular, fusiform, and mycotic). Recently, the study of blood flow interaction with aneurysms has gained attention from physicians and engineers. Shear stresses, oscillatory shear index (OSI), gradient oscillatory number (GON), and residence time have been used as variables to describe the hemodynamics as well as the origin and evolution of aneurysms. However, the causes and hemodynamic conditions that promote their growth are still under debate. The present work presents numerical simulations of three types of aneurysms: two aortic and one cerebral. Simulation results showed that the blood rheology is not relevant for aortic aneurysms. However, for the cerebral aneurysm case, blood rheology could play a relevant role in the hemodynamics. The evaluated turbulence models showed equivalent results in both cases. Lastly, a simulation considering the fluid–structure interaction (FSI) showed that this phenomenon is the dominant factor for aneurysm simulation.

59 BASIC BIOLOGICAL SCIENCES↗

The Effect of Air Separations on Fast Pyrolysis Products for Forest Residue Feedstocks

This study investigates the intricate relationship between biomass preprocessing and pyrolysis product yields, employing the air classification technique for the treatment of loblolly pine residues with varying moisture content. A comprehensive exploration of the physicochemical properties of air-classified loblolly pine informs a sophisticated pyrolysis simulation model. Given the complex and multifaceted nature of biomass pyrolysis, operating across diverse temporal and spatial scales, a pyrolysis kinetics-based CFD–DEM simulation method is employed to predict product yields. Results showed that the elevated moisture content amplifies particle adhesiveness, necessitating augmented air velocities for effective separation, thereby influencing the efficiency of the separation process. While carbon and hydrogen contents exhibit relative stability across diverse moisture contents and blower frequencies, the oxygen content undergoes noticeable changes. For example, the oxygen contents were measured as 29.2 and 38.6 wt% in the light fraction of 30% moisture content sample at blower frequencies of 10 and 20 Hz, respectively. An intriguing finding emerges from pyrolysis simulation, indicating that a lower blower frequency in air classification moderately enhances bio-oil yield and significantly improves its quality, particularly in terms of water content. For instance, the water content in the bio-oil was about 1.5% and 10% in the heavy and light fractions, respectively from 10% moisture sample under 15 Hz blower frequency. In summary, a detailed understanding and strategic manipulation of critical material attributes in biomass through efficient fractionation techniques are imperative for advancing fast pyrolysis as a sustainable avenue for renewable energy and chemical production.

09 BIOMASS FUELS↗

Applying a Compact Porous Media Model to Numerically Derive Resistance Coefficients for Lattice Structures

Additive Manufacturing allows for exploring various geometries to achieve specific engineering criteria. Lattices are one geometry with unique properties, including being periodically repeating structures which allow flow through them to be represented as a porous media according to Darcy-Forchheimer equations. These equation’s coefficients are generally experimentally derived, but this work demonstrates the ability to numerically derive them with CFD. Simulations were performed using three-dimensional stead state Reynolds-averaged Navier-Stokes with a k-ω Shear Stress Transport turbulence model using Ansys Fluent. Three lattice geometries were investigated and drag coefficients were derived. The method was validated against externally published data for similar geometries demonstrating strong agreement, and grid convergence for all simulations was calculated with a Grid Convergence Index method. Wall roughness is demonstrated to have a non-negligible impact on results and roughness values are considered for the primary focus Octahedral geometry where both smooth wall and rough wall coefficients were derived. The porosity coefficients for the Octahedral geometry at 1.0 [m/s] were found to be 2.89×10 6 and 2.90×10 6 [1/(Pa*m*s)] for the permeability coefficients, 6.37×10 1 and 5.44×10 1 [m 2 /kg] for the inertial resistance coefficients, and with a max pressure drop of 5116.7 [Pa] and 4429.5 [Pa] for the smooth walls and rough walls, respectively. The derived numerical method enables rapid exploration and optimization of new lattice designs for diverse engineering applications.

42 ENGINEERING↗

High-Fidelity Simulation Aerodynamics Dataset of NACA 0012 and 0021 Airfoils

The repository contains time series data of pressure, viscous and moment forces for the NACA 0012 and 0021 airfoils at angles of attack (AOA), 5, 17, 30, 45, 60 and 90 degrees. The data was generated during the high-fidelity simulations of the airfoils using the open-source CFD code, Nalu-Wind (https://github.com/Exawind/nalu-wind). The research objective was to study the effects of mesh resolution and turbulence model on the three-dimensional deep stall aerodynamics of the airfoils, which has been published in the Journal of Turbulence (https://doi.org/10.1080/14685248.2023.2225141). Several meshes of the O-grid type were generated using the modified Pointwise Glyph script of Carrigan (https://github.com/pointwise/AirfoilMesh) by varying the wall normal and spanwise resolutions. Note that, the geometry was extruded by four chord lengths in the spanwise direction. In the paper, we showed that Improved Delayed Detached Eddy Simulation (IDDES) hybrid RANS/LES turbulence model with a minimum of 24 cells per chord length in the spanwise direction is necessary to correctly predict the loads in the deep stall regime. Under the NACA 0012 airfoil directory, the data for several combinations of wall-normal and spanwise resolutions are presented in each of the AOA subdirectories. The corresponding experimental data for the lift and drag polars can be found in the "exp_data" subdirectory. The directory for the NACA 0021 airfoil consists of CFD force data and airfoil polars for Reynolds numbers 2.7e5 and 2.0e6.

17 WIND ENERGY↗

Simulation and Analysis on Reactor Pressure Vessel (RPV) subjected to Pressurized Thermal Shock (PTS) under SBLOCA scenario by using Cardinal to support the fracture mechanics analyses

The structural components that comprise nuclear reactors and their supporting structures are subjected to harsh operating environments that can challenge their integrity, especially after exposure for extended durations or under accident condition. As one of the most significant components of a Reactor, the Reactor Pressure Vessel (RPV) is exposed to an aggressive environment during the operation time (e.g. more than 40 years). Ageing degradation mechanisms (e.g. thermos-fatigue) could grow initial defects up to a critical size, increasing the susceptibility to failure in the RPV. The conventional methods are mostly based on simple crack and structure geometries. Very limited studies consider the real conditions of the RPV subjected to a thermal shock due to a Loss of Coolant Accident (LOCA). During a LOCA event, the most severe conditions take place when the emergency core cooling (ECC) water is injected inside the cold legs filled initially with hotter water and/or steam. The rapid cooling of the down-comer and the internal RPV surface followed probably by re-pressurization of the RPV causes large temperature gradients and variation of pressure which induces thermal-mechanical stresses. In order to develop the model for integrity assessment of a reactor pressure vessel (RPV) subjected to pressurized thermal shock (PTS), a multi-physics simulation, which includes the thermo-hydraulic, thermo-mechanical and fracture mechanics analyses is necessary. The multi-physics simulations are performed using Cardinal, a wrapping of the GPU-oriented spectral element Computational Fluid Dynamics (CFD) code NekRS and other multi-physics sub-modules within the MOOSE framework. Cardinal now fully supports MOOSE stochastic perturbations of NekRS models with varying boundary conditions, initial conditions, material properties, and any other quantity which is defined by a kernel (such as coefficients in a momentum source model). The implementation is designed in a flexible manner so that scalar values are sent from MOOSE into a user scratch space in NekRS, which can then be applied for any purpose within the NekRS case files (both on the host and device). When modeling PTS, several factors can impact the results significantly. In this report, the impacts of the geometry of the model, Reynolds number and buoyancy effect are investigated. Two geometry, i.e., a simplified model and a realistic RPV model, with both laminar and turbulent flow condition are adopted for the PTS simulation with and without buoyancy effect. The purpose of the investigation is to understand the impact of these factors on the prediction of temperature history of RPV. The accurate prediction on the temperature evolution, which will be exported to Grizzly code for further analyses on the progression of aging mechanisms and their effects on the integrity of RPV structures, is very crucial. Based on the understanding of these factors, a more sophisticated model is built to analysis the PTS under SBLOCA scenario. A literature survey is conducted to pick the SBLOCA scenario for the multi-physics simulation. The analysis helps to explain the form and the transformation of the cold plum when the ECC is activated under SBLOCA. This model can be can be applied to study the PTS effect for different RPV configurations. The results can help to assess structural component degradation for advanced reactors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Hydrogen Equipment Enclosure Risk Reduction through Earlier Detection of Component Failures

Hydrogen component reliability and the hazard associated with failure rates is a critical area of research for the successful implementation and growth of hydrogen technology across the globe. The research team has partnered to quantify system risk reduction through earlier detection of hydrogen component failures. A model of hydrogen dispersion in a hydrogen equipment enclosure has been developed utilizing experimentally quantified hydrogen component leak rates as inputs. This model provides insight into the impact of hydrogen safety sensors and ventilation on the flammable mass within a hydrogen equipment enclosure. This model also demonstrates the change in safety sensor response time due to detector placement under various leak scenarios. The team looks to improve overall hydrogen system safety through an improved understanding of hydrogen component reliability and risk mitigation methods. This collaboration fits under the work program of IEA Hydrogen Task 43 Subtask E Hydrogen System Safety.

CFD↗

SMR safety through HTTF modeling and benchmark efforts for code validation for gas-cooled reactor applications

Accurate modeling and simulation tools for thermal-hydraulics calculations are a key element needed to design and license new advanced reactors including Small Modular Reactors (SMR) and Microreactors. Uncertainties in modeling and simulation can have significant safety and economic implications. The High Temperature Test Facility (HTTF) at Oregon State University (OSU) is a scaled integral effects experiment designed to investigate transient behavior in high-temperature gas-cooled prismatic-block nuclear reactors. High-quality measurement data is available from the HTTF that is suitable for a thermal-hydraulics code validation benchmark for gas-cooled reactor simulations. Here, this paper summarizes individual HTTF modeling efforts to date for tool validation at Idaho National Laboratory (INL), Argonne National Laboratory (ANL), Oregon State University (OSU) and Canadian Nuclear Laboratories (CNL) using system thermal-hydraulics codes, Computational Fluid Dynamics (CFD) codes and system-CFD code couplings. Also, the paper introduces the ongoing OECD Nuclear Energy Agency (NEA) High Temperature Gas Reactor Thermal-Hydraulics (HTGR T/H) benchmark that allows for better comparisons of results between different international modeling teams. The benchmark provides well defined computational problems that include code-to-code comparisons and comparisons to measured data. These problems provide an avenue for quantifying accuracy and identifying sources of uncertainty in thermal-hydraulics calculations, including in measured thermophysical properties, as part of validation for gas-cooled reactor simulation tools.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

Experimental and Computational Heat Transfer Study of sCO 2 Single-Jet Impingement

This study experimentally and computationally investigates the heat transfer capability of supercritical carbon dioxide (sCO 2 ) single jet impingement. The evaluated jet Reynolds number range is between 80,000 and 600,000, with a nondimensional jet-to-target surface spacing of 2.8. CO 2 -impinging jet stagnation conditions were maintained at approximately 200 bar and a temperature of ∼400 °C for most experiments. The goal is to understand how changes in the aforementioned parameters influence heat transfer between the working fluid and the heated surface. Additionally, due to the elevated Reynolds numbers and difference in thermodynamic properties between air and CO 2 , air-derived impingement correlations may not be appropriate for CO 2 impingement; these correlations will be evaluated against experimental sCO 2 impingement data. At the time of this study, no sCO 2 impingement data was available relevant to sCO 2 power cycles. The target surface is a 1.5-in. diameter copper block centered on the 3 mm jet orifice. A mica heating element bolted to the bottom of the copper block provides a uniform heat flux. Thermocouples embedded in the copper block are used to determine the surface temperature. Nusselt numbers obtained from experimental sCO 2 test data are compared to area-averaged Nusselt numbers from air-derived correlations. The comparisons showed that air correlations drastically underpredict the heat transfer when sCO 2 is used as the working fluid. A modified sCO 2 correlation using experimental data at discussed conditions is derived based on an existing air correlation. A CFD study is also performed to further investigate sCO 2 heat transfer characteristics, and assess the numerical model applicability to this problem type.

42 ENGINEERING↗

CFD Study of Countercurrent Flow in Triply Periodic Minimal Surfaces with CO 2 BOL Solvent

Triply periodic minimal surfaces (TPMS) are promising structures for heat exchangers and ultrafiltration due to their high heat and mass transfer efficiency compared to traditional devices. Because of this, the TPMS structures are a potential candidate for building the multifunctional, intensified devices in CO 2 capture tasks. With the Lawrence Livermore National Laboratory’s (LLNL’s) additive manufacturing (AM) capability, these complicated TPMS structures can be printed in the packed column for experimental tests. To optimize the TPMS structure and enable fast prototyping, a computationally-efficient computational fluid dynamics (CFD) framework was established to explore the TPMS structure’s hydrodynamics and mass transfer performance. The established model was able to simulate the countercurrent flows by integrating the Pacific Northwest National Laboratory (PNNL) CO 2 binding organic liquids (CO 2 BOL) solvent into the TPMS structures. Three typical TPMS structures - namely, Gyroid, Schwarz-D, and Schwarz-P - were investigated using the periodic boundary set up in all three directions. Two different unit cell size (1 cm and 2 cm) geometries are tested for each type of TPMS. The solvent flow rate in this study covered a range of [0, 0.1] m/s with an initial liquid film thickness of [0.1, 1.8] mm. These simulations provided a preliminary understanding of the TPMS structure behaviors in countercurrent flow conditions. The study also related the TPMS geometry requirement to the given CO 2 BOL solvent physical property and covered reasonable packed column operation ranges.

36 MATERIALS SCIENCE↗

Hydrodynamic Analysis and Optimization of Aquantis Marine Turbine: Cooperative Research and Development (Final Report)

The primary aim of this proposal is to improve the accurate prediction of hydrodynamic performance and dynamic load responses of the AQ10 floating axial-flow tidal turbine with a tri-cat mooring configuration. The validation of reduced-order modeling approaches with high-fidelity model will be implemented. Additionally, the frequency response domain, Response Amplitude Floating Wind (RAFT) toolbox plus an optimizer expanded for marine hydrokinetic turbines under the Submarine Hydrokinetic And Riverine Kilo-megawatt. Systems (SHARKS) program will be used for designing and exploring different key design parameters (platform dimension, mooring layout and its parameters) of next marine hydrokinetic (MHK) turbine generation.

16 TIDAL AND WAVE POWER↗

Latent heat thermal energy storage performance maps enabling fast & accurate building energy simulations

Thermal energy storage (TES) using phase change materials (PCMs) has gained attention as an effective approach to manage energy demand fluctuations and shift peak building loads. PCM embedded heat exchangers (PCM-HXs) offer high energy storage density and low temperature variation during phase change, being suitable for load-shifting applications. However, this component is typically evaluated using computationally expensive methods, which present significant challenges when the ultimate goal is to assess the performance of PCM-HX integrated thermal energy storage systems in the full building context. In this paper, we present a methodology to generate highly accurate and computationally efficient PCM-HX performance maps which can be easily integrated into building energy simulation tools to analyze the feasibility of space conditioning systems with latent heat PCM-based TES. The performance maps are generated using a computationally efficient PCM-HX simulation tool based on a Generalized Resistance-Capacitance Model (GRCM) which can simulate arbitrary PCM-HXs with high accuracy and significantly less computational effort compared to full CFD simulations. The methodology was verified for a case study considering a 5-ton (~17.5 kW) air-to-water heat pump-thermal energy storage system (HP-TES), which was co-simulated in Modelica for a DOE prototype small-office building in Vienna, Austria, using Spawn of EnergyPlus™. The TES performance maps provided accurate predictions of PCM-HX behavior when used as Modelica component, with deviations within 2-4% while also achieving at least 103 computational time reduction. Leveraging this faster prediction capability, four PCMs with different melting temperatures for cooling (12°C, 16°C) and heating (31°C, 36°C) were assessed to investigate their impact on system performance. This work highlights the importance of robust PCM-HX models for efficient and high-fidelity building-level simulations, presenting new opportunities for advanced control strategy development and parametric analysis of TES configurations in a computationally efficient manner

Modelica Building Simulations↗

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↗

Coupled model for liquid lithium plasma facing components

Numerical analysis provides the design choice and operating window of liquid metal Plasma Facing Components (PFC) concepts. Coupled analysis of boundary plasma together with the surrounding boundary structures is required. Here, to achieve this goal, PPPL is developing a comprehensive multi-physics model for modeling of PFCs in fusion devices. The model includes the fluid-kinetic code SOLPS-ITER and the flow and heat transfer code CFX from ANSYS. SOLPS-ITER was augmented with a liquid metal boundary condition algorithm, allowing direct two-way coupling of the plasma analysis with the two-dimensional analytical slab flow model which includes heat convection in the liquid metal PFC. The target heat flux resulting from this coupled analysis is used as a boundary condition for detailed 3D Computational Fluid Dynamics (CFD) Magneto Hydro Dynamics (MHD) and heat transfer analysis. A new formulation of MHD equations is introduced in the numerical procedure ensuring current conservation of the discretized equations. Results of the 3D analysis are used for final validation of the coupled model. A PFC design where a porous wall is used to stabilize the liquid metal surface, while MHD drive is used to push the liquid metal flow inside the PFC, will be investigated in the regimes where vapor shielding is created for enhanced volumetric plasma heat dissipation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Thermal-Hydrological-Mechanical Modelling of Stockton University Reservoir Cooling System, Fine Scale Stress Test Modelling

Mesh, properties, initial conditions, injection/withdrawal rates for modelling thermal, hydrological, and mechanical effects of fluid injection to and withdrawal from ground for Stockton University reservoir cooling system (aquifer storage cooling system), Galloway, New Jersey, for unscheduled two hour injection at 133 % designed capacity, on fine scale grid, with some results. Second simulation of J.T. Smith, E. Sonnenthal, P. Dobson, P. Nico, and M. Worthington, 2021. Thermal-hydrological-mechanical modeling of Stockton University reservoir cooling system, Proceedings of the 46th Workshop on Geothermal Reservoir Engineering, Stanford University, SGP-TR-218, from which Figures 6-9, pertain.

15 GEOTHERMAL ENERGY↗

Q3 Report for FY25 Theory and Simulation Performance Target: Development of an integrated modeling framework for fusion reactor design and assessment

This report describes the work and activities carried out towards the completion of each of the following milestones in FY25 Q3: 1. Demonstrate workflow for generating self-consistent CESOL plasma profiles + first wall and divertor loading prediction and generate the CAT plasma and neutron loading needed for further engineering analysis. • Benchmark between two first wall heat flux mapping methods, identify importance of various heat flux sources and physics impact of using fully coupled CESOL vs post-analysis evaluation. 2. Generate medium fidelity parametrized CAD. • Generate parametrized CAD components for the CAT example case via either user-defined modules called within the geometry generation or by defeatured/parametrized CAD, including DCLL blanket matched to divertor boundary and magnets. Define materials, labels, and boundary conditions for passing the mesh to CFD tools. 3. Demonstrate multiphysics magnet analysis. • Demonstrate magnet analysis workflow called from the FREDA workflow, and 4. Demonstrate nuclear analysis. • Add model to OpenFOAM and/or other codes possibly including Diablo to account for tritium diffusion in solids.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Understanding flow dynamics in membrane distillation: Effects of reactor design on polarization

Optimization and design of full-scale membrane distillation (MD) systems usually require Sherwood and Nusselt correlations that are developed from lab-scale systems. However, entrance effects in lab-scale systems can significantly impact heat, mass and momentum transfer in the reactor, therefore affect the accuracy of the developed experimental Sherwood and Nusselt correlations. Here, Computational Fluid Dynamics (CFD) simulations using OpenFOAM are performed to understand the effects of right-angled bends and inlet design on flow dynamics, temperature and concentration polarization in MD systems. Simulation results show that the presence of right-angled bends and inlets with sudden expansions lead to the formation of Dean vortices. Dean vortices enhance perpendicular mixing in MD systems and reduce both temperature and concentration polarization. Temperature and concentration polarization coefficients in MD systems with right-angled bends and inlets with sudden expansions vary significantly for the same volumetric flow rate. Our studies show that lab-scale systems with the same volumetric flow rate but different designs lead to significantly different Nusselt and Sherwood correlations. Here, this study demonstrates the importance of CFD-informed design of lab-scale systems to minimize entrance effects and suppress Dean vortices for consistent model development and calibration across multiple scales.

42 ENGINEERING↗