Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FORMULATION”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

An implicit-in-time DPG formulation of the 1D1V Vlasov-Poisson equations

Efficient solution of the Vlasov equation, which can be up to six-dimensional, is key to the simulation of many difficult problems in plasma physics. The discontinuous Petrov-Galerkin (DPG) finite element methodology provides a framework for the development of stable (in the sense of Ladyzhenskaya–Babuška–Brezzi conditions) finite element formulations, with built-in mechanisms for adaptivity. While DPG has been studied extensively in the context of steady-state problems and to a lesser extent with space-time discretizations of transient problems, relatively little attention has been paid to time-marching approaches. In the present work, we study a first application of time-marching DPG to the Vlasov equation, using backward Euler for a Vlasov-Poisson discretization. We demonstrate adaptive mesh refinement for two problems: the two-stream instability problem, and a cold diode problem. Furthermore, we believe the present work is novel both in its application of unstructured adaptive mesh refinement (as opposed to block-structured adaptivity, which has been studied previously) in the context of Vlasov-Poisson, as well as in its application of DPG to the Vlasov-Poisson system. We also discuss extensive additions to the Camellia library in support of both the present formulation as well as extensions to higher dimensions, Maxwell equations, and space-time formulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Modelling of transport processes in concrete exposed to elevated temperatures – An alternative formulation for sorption isotherms

There is a significant need to understand, analyse and assess moisture transport in cementitious materials exposed to elevated temperatures in order to confidently predict the behaviour and ultimately the development of damage in safety critical applications such as nuclear reactor vessels, structures exposed to fire and well bore grouts. In view of this need a rigorous and robust formulation to describe water retention curves (sorption isotherms) as a function of temperature based on the evolution of physical parameters is presented. The model formulation is successfully validated against independent sets of experimental data up to temperatures of 80 °C. It is then further validated under isothermal drying conditions and then high temperature conditions through the numerical reproduction of laboratory experiments following implementation in a fully coupled hygro-thermo-mechanical finite element model. The new formulation is found to work well under a variety of conditions in a variety of cementitious material types.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

A reconstructed discontinuous Galerkin method for compressible flows in Lagrangian formulation

In this work, a high-order accurate reconstructed discontinuous Galerkin (rDG) method is developed for solving two-dimensional hydrodynamic problems in cell-centered updated Lagrangian formulation. This method is the Lagrangian limit of the unsplit rDG-ALE formulation, and is obtained by assuming the equality of the grid velocity to the fluid velocity only at cell boundaries. The conservative variables and the Taylor basis defined on the time-dependent moving mesh, provide the piece-wise polynomial expansion in the updated Lagrangian formulation. A multi-directional nodal Riemann solver is implemented for computing the grid velocity at the vertices and the numerical flux at the cell boundaries. A characteristic limiting procedure is extended from the primitive variable version to the conservative variable version, and its performance is compared with the limiter on physical variables. A number of benchmark test cases are conducted to assess the accuracy, robustness, and non-oscillatory property of the DG(P0), DG(P1) and rDG(P1P2) methods. The numerical experiments demonstrate that the developed rDG method is able to attain the designed order of accuracy and the characteristic limiting procedure outperforms the limiter on physical variables in terms of the monotonicity and symmetry preservation for shock problems.

42 ENGINEERING↗
DOE OSTI2015

A thermodynamically consistent discontinuous Galerkin formulation for interface separation

Our paper describes the formulation of an interface damage model, based on the discontinuous Galerkin (DG) method, for the simulation of failure and crack propagation in laminated structures. The DG formulation avoids common difficulties associated with cohesive elements. Specifically, it does not introduce any artificial interfacial compliance and, in explicit dynamic analysis, it leads to a stable time increment size which is unaffected by the presence of stiff massless interfaces. This proposed method is implemented in a finite element setting. Convergence and accuracy are demonstrated in Mode I and mixed-mode delamination in both static and dynamic analyses. Significantly, numerical results obtained using the proposed interface model are found to be independent of the value of the penalty factor that characterizes the DG formulation. By contrast, numerical results obtained using a classical cohesive method are found to be dependent on the cohesive penalty stiffnesses. The proposed approach is shown to yield more accurate predictions pertaining to crack propagation under mixed-mode fracture because of the advantage. Furthermore, in explicit dynamic analysis, the stable time increment size calculated with the proposed method is found to be an order of magnitude larger than the maximum allowable value for classical cohesive elements.

42 ENGINEERING↗
DOE OSTI2023

Piecewise linear approximation with minimum number of linear segments and minimum error: A fast approach to tighten and warm start the hierarchical mixed integer formulation

In several areas of economics and engineering, it is often necessary to fit discrete data points or approximate nonlinear functions with continuous functions. Piecewise linear (PWL) functions are a convenient way to achieve this. PWL functions can be modeled in mathematical problems using only linear and integer variables. Moreover, there is a computational benefit in using PWL functions that have the least possible number of segments. This work proposes a novel hierarchical mixed integer linear programming (MILP) formulation that identifies a continuous PWL approximation with minimum number of linear segments for a given target maximum error. The proposed MILP formulation also identifies the solution with the least maximum error among the solutions with minimum number of segments. Then, this work proposes a fast iterative algorithm that identifies non necessarily continuous PWL approximations by solving O(S log N) linear programming (LP) problems, where N is the number of data points and S is the minimum number of segments in the non necessarily continuous case. This work demonstrates that tight bounds for the MILP problem can be derived from these approximations. Next, a fast algorithm is introduced to transform a non necessarily continuous PWL approximation into a continuous one. Finally, the tight bounds and the continuous PWL approximations are used to tighten and warm start the MILP problem. The tightened formulation is shown in experimental results to be more efficient, especially for large data sets, with a solution time that is up to two orders of magnitude less than the existing literature.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

On the convexity of phase-field fracture formulations: Analytical study and comparison of various degradation functions

Efficient and accurate fracture modeling is of great importance in applications where catastrophic outcomes under extreme scenarios are possible. The phase-field (PF) approach to fracture received significant attention over the past decade, due to its capability to capture complicated fracture patterns (e.g., crack merging and branching). Specifically, crack initiation and propagation are modeled via minimization of the total energy functional, which is regularized with the aid of a phase field. Despite the promising results and modeling capabilities of the PF method in many applications, the solution of fracture problems remains computationally challenging mainly due to the non-convexity of the total energy functional with respect to the combined unknown (phase field and displacement) fields. Understanding the effects of their coupling on convexity is crucial in order to address frequently encountered hurdles in fracture modeling (e.g., inefficient solvers and non-physical crack nucleation). In this paper, we develop convexity criteria for a wide class of PF fracture formulations. For this class of formulations, the second variation of the total energy functional is expressed in terms of Hessian matrices (evaluated at individual material points). Depending on the choice of geometric crack functions and degradation functions, we classify the formulations into three categories and analytically study each one separately. To study the sign of the second variation, we derive inequalities which are satisfied at material points when the Hessian matrix is locally positive semi-definite. These inequalities provide objective criteria for comparing degradation functions. Finally, the applicability of the proposed convexity criteria is demonstrated in the context of a one-dimensional problem, solved using a conventional monolithic solver.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Modeling of intragranular misorientation and grain fragmentation in polycrystalline materials using the viscoplastic self-consistent formulation

In this paper, the recently established methodology to use known algorithmic expressions of the second moments of the stress field in the grains of a polycrystalline aggregate for calculating average fluctuations of lattice rotation rates and the associated average intragranular misorientation distributions using the mean-field viscoplastic self-consistent (VPSC) formulation is extended to solve the coupled problem of considering the effect of intragranular misorientations on stress and rotation rate fluctuations. In turn, these coupled expressions are used to formulate and implement a grain fragmentation (GF) model in VPSC. Case studies, including tension and plane-strain compression of face-centered cubic polycrystals are used to illustrate the capabilities of the new model. GF-VPSC predictions of intragranular misorientation distributions and texture evolution are compared with experiments and full-field numerical simulations, showing good agreement. In particular, the inclusion of misorientation spreads reduced the intensity of the deformed texture and thus improved the texture predictions. Finally and moreover, considering that intragranular misorientations act as driving forces for recrystallization, the new GF-VPSC formulation is shown to enable modeling of microstructure evolution during deformation and recrystallization, in a computationally efficient manner.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Excipient screening by lyophilization provides insights into spray drying formulations for nanoparticle vaccines

Nanoparticles have shown great promise as delivery platforms in the development of tunable and safe vaccines. Nanolipoprotein particles (NLPs), also known as nanodiscs, are discoidal nanoparticles composed of a lipid bilayer stabilized at their periphery by apolipoproteins. Under the right conditions, the NLP self-assembly process is highly customizable in terms of lipids and apolipoprotein constituents, allowing for tunable physical and chemical characteristics. This flexibility allows a wide range of vaccine antigens and adjuvants to be incorporated onto the NLP platform for tailored vaccine design. The stability of NLPs during long term storage is a very important factor in developing a vaccine delivery platform suitable for widespread global use. When stored in a solution for extended periods of time, NLPs dissociate into their corresponding lipids and protein constituents, leading to particle degradation. Proper stabilization of NLPs can often be achieved by lyophilization (i.e. freeze-drying), a method widely used for various applications including pharmaceuticals. This process, however, can be damaging to particles without the presence of lyoprotectants or excipients that help maintain particle stability during lyophilization. Another method used to stabilize vaccines and pharmaceuticals is spray drying, a process that converts liquid formulations into dry powders through controlled heating and airflow. While spray drying is rapid, scalable, and cost-effective, lyophilization is typically a gentler process that better retains biomolecule structure and function. Both processes use excipients for particle stabilization, so lyophilization can be used as a surrogate to test stability of NLPs, to down-select formulations that may withstand the harsher conditions of spray drying. To screen different formulations, NLPs were synthesized and purified to homogeneity by size exclusion chromatography (SEC) and samples were prepared with a wide range of lyoprotectants and/or excipients. To assess the protective effects of excipients on NLPs upon spray drying, both pre- and post-lyophilized samples were analyzed by SEC. To assess the protective effects upon heating (encountered during the spray drying process), NLP samples were incubated at elevated temperatures prior to SEC analysis. The lyoprotectants and excipients evaluated in this study had different efficiencies in protecting NLPs during lyophilization and heating tests. Trehalose, for example, exhibits stabilization on NLPs both upon lyophilization and heating whereas leucine accelerated NLP dissociation. Although some lyoprotectants are effective by themselves, different combinations can decrease the stabilization of NLPs. Shelf-stable vaccines that do not require cold-chain storage are essential for global accessibility and our findings provide fundamental insight into how to advance NLP-based vaccines for these applications.

Serrano, Litzay J [Lawrence Livermore National Lab↗
DOE OSTI2015

An adhesive contact mechanics formulation based on atomistically induced surface traction

In this work, we have developed a novel multiscale computational contact formulation based on the generalized Derjuguin approximation for continua that are characterized by atomistically enriched constitutive relations in order to study macroscopic interaction between arbitrarily shaped deformable continua. The proposed adhesive contact formulation makes use of the microscopic interaction forces between individual particles in the interacting bodies. In particular, the double-layer volume integral describing the contact interaction (energy, force vector, matrix) is converted into a double-layer surface integral through a mathematically consistent approach that employs the divergence theorem and a special partitioning technique. The proposed contact model is formulated in the nonlinear continuum mechanics framework and implemented using the standard finite element method. With no large penalty constant, the stiffness matrix of the system will in general be well-conditioned, which is of great significance for quasi-static analysis. Three numerical examples are presented to illustrate the capability of the proposed method. Results indicate that with the same mesh configuration, the finite element computation based on the surface integral approach is faster and more accurate than the volume integral based approach. In addition, the proposed approach is energy preserving even in a very long dynamic simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Higher-order finite-difference formulation of periodic Orbital-free Density Functional Theory

We present a real-space formulation and higher-order finite-difference implementation of periodic Orbital-free Density Functional Theory (OF-DFT). Specifically, utilizing a local reformulation of the electrostatic and kernel terms, we develop a generalized framework for performing OF-DFT simulations with different variants of the electronic kinetic energy. In particular, we propose a self-consistent field (SCF) type fixed-point method for calculations involving linear-response kinetic energy functionals. In this framework, evaluation of both the electronic ground-state and forces on the nuclei are amenable to computations that scale linearly with the number of atoms. We develop a parallel implementation of this formulation using the finite-difference discretization. We demonstrate that higher-order finite-differences can achieve relatively large convergence rates with respect to mesh-size in both the energies and forces. Additionally, we establish that the fixed-point iteration converges rapidly, and that it can be further accelerated using extrapolation techniques like Anderson's mixing. We validate the accuracy of the results by comparing the energies and forces with plane-wave methods for selected examples, including the vacancy formation energy in Aluminum. Overall, the suitability of the proposed formulation for scalable high performance computing makes it an attractive choice for large-scale OF-DFT calculations consisting of thousands of atoms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

WeakIdent: Weak formulation for identifying differential equation using narrow-fit and trimming

Data-driven identification of differential equations is an interesting but challenging problem, especially when the given data are corrupted by noise. When the governing differential equation is a linear combination of various differential terms, the identification problem can be formulated as solving a linear system, with the feature matrix consisting of linear and nonlinear terms multiplied by a coefficient vector. This product is equal to the time derivative term, and thus generates dynamical behaviors. The goal is to identify the correct terms that form the equation to capture the dynamics of the given data. We propose a general and robust framework to recover differential equations using a weak formulation with two new mechanisms, narrow-fit and trimming, for both ordinary and partial differential equations (ODEs and PDEs). The weak formulation facilitates an efficient and robust way to handle noise, and two new mechanisms, narrow-fit and trimming, improve the coefficient support and value recoveries respectively. For each sparsity level, Subspace Pursuit is utilized to find an initial set of support from the large dictionary. Then, we focus on highly dynamic regions (rows of the feature matrix), and error normalize the feature matrix in the narrow-fit step. The support is further updated via trimming the terms that contribute the least. Finally, the support set of features with the smallest Cross-Validation error is chosen as the result. A comprehensive set of numerical experiments are presented for both systems of ODEs and PDEs with various noise levels. The proposed method gives a robust recovery of the coefficients, and a significant denoising effect which can handle up to 100% noise-to-signal ratio for some equations. We compare the proposed method with several state-of-the-art algorithms for the recovery of differential equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Primal interface debonding formulation for finite strain isotropic plasticity

In this work, a framework is developed for modeling ductile damage of nonlinear materials whose plastic deformation is characterized using rate independent classical plasticity. This method relies on the assumption that the free energy can be decomposed into elastic, plastic and damage parts. A thermodynamically consistent method is derived which satisfies the second law of thermodynamics in the Clausius–Duhem inequality form. The dissipation associated with plasticity takes place in the domain only, while damage dissipation is localized to the interface. The method is developed using Variational Multiscale ideas to obtain definitions of the interface fluxes within a primal formulation analogous to the Discontinuous Galerkin method, which ensures weakly vanishing interface gap prior to reaching a damage initiation criterion. The local nonlinear problem to calculate both plastic deformation gradient and damage variable follows an incremental approach similar to classical plasticity return mapping algorithm. This elastoplastic damage formulation is developed for material undergoing finite strain, and it naturally accommodates a trapezoidal traction separation law (TSL) whose shape can be varied to model either ductile interface behavior or brittle interface behavior. The formulation's performance is assessed through modeling a patch test and a compact tension specimen.

42 ENGINEERING↗
DOE OSTI2013

Aggregate material formulated with MSWI bottom ash and APC fly ash for use as secondary building material

Highlights: ► A concrete formulation was optimized using Bottom Ash and APC ash. ► 10% of APC ash achieves good compromise between economic and performance aspects. ► The crushed concrete was evaluated as secondary building granular material. ► The environmental behavior allows its use as secondary material. ► The abrasion resistance is not good enough for its use as a road sub-base material. - Abstract: The main goal of this paper is to obtain a granular material formulated with Municipal Solid Waste Incineration (MSWI) bottom ash (BA) and air pollution control (APC) fly ash to be used as secondary building material. Previously, an optimum concrete mixture using both MSWI residues as aggregates was formulated. A compromise between the environmental behavior whilst maximizing the reuse of APC fly ash was considered and assessed. Unconfined compressive strength and abrasion resistance values were measured in order to evaluate the mechanical properties. From these results, the granular mixture was not suited for certain applications owing to the high BA/APC fly ash content and low cement percentages used to reduce the costs of the final product. Nevertheless, the leaching test performed showed that the concentrations of all heavy metals were below the limits established by the current Catalan legislation for their reutilization. Therefore, the material studied might be mainly used in embankments, where high mechanical properties are not needed and environmental safety is assured.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2018

Laguerre–Hermite pseudo-spectral velocity formulation of gyrokinetics

First-principles simulations of tokamak turbulence have proven to be of great value in recent decades. We develop a pseudo-spectral velocity formulation of the turbulence equations that smoothly interpolates between the highly efficient but lower resolution three-dimensional (3-D) gyrofluid representation and the conventional but more expensive 5-D gyrokinetic representation. Our formulation is a projection of the nonlinear gyrokinetic equation onto a Laguerre–Hermite velocity-space basis. We discuss issues related to collisions, closures and entropy. While any collision operator can be used in the formulation, we highlight a model operator that has a particularly sparse Laguerre–Hermite representation, while satisfying conservation laws and the H theorem. Free streaming, magnetic drifts and nonlinear phase mixing each give rise to closure problems, which we discuss in relation to the instabilities of interest and to free energy conservation. We show that the model is capable of reproducing gyrokinetic results for linear instabilities and zonal flow dynamics. Thus the final model is appropriate for the study of instabilities, turbulence and transport in a wide range of geometries, including tokamaks and stellarators.

Physics↗
DOE OSTI2023

Flash NanoPrecipitation as an Agrochemical Nanocarrier Formulation Platform: Phloem Uptake and Translocation after Foliar Administration

The increasing severity of pathogenic and environmental stressors that negatively affect plant health has led to interest in developing next-generation agrochemical delivery systems capable of precisely transporting active agents to specific sites within plants. In this work, we adapt Flash NanoPrecipitation (FNP), a scalable nanocarrier (NC) formulation technology used in the pharmaceutical industry, to prepare organic core–shell NCs and study their efficacy as foliar or root delivery vehicles. NCs ranging in diameter from 55 to 200 nm, with surface zeta potentials from -40 to +40 mV, and with seven different shell material properties were prepared and studied. Shell materials included synthetic polymers poly(acrylic acid), poly(ethylene glycol), and poly(2-(dimethylamino)ethyl methacrylate), naturally occurring compounds fish gelatin and soybean lecithin, and semisynthetic hydroxypropyl methylcellulose acetate succinate (HPMCAS). NC cores contained a gadolinium tracer for tracking by mass spectrometry, a fluorescent dye for tracking by confocal microscopy, and model hydrophobic compounds (alpha tocopherol acetate and polystyrene) that could be replaced by agrochemical payloads in subsequent applications. After foliar application onto tomato plants with Silwet L-77 surfactant, internalization efficiencies of up to 85% and NC translocation efficiencies of up to 32% were observed. Significant NC trafficking to the stem and roots suggests a high degree of phloem loading for some of these formulations. Results were corroborated by confocal microscopy and synchrotron X-ray fluorescence mapping. NCs stabilized by cellulosic HPMCAS exhibited the highest degree of translocation, followed by formulations with a significant surface charge. The results from this work indicate that biocompatible materials like HPMCAS are promising agrochemical delivery vehicles in an industrially viable pharmaceutical nanoformulation process (FNP) and shed light on the optimal properties of organic NCs for efficient foliar uptake, translocation, and delivery.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Formulation and Performance Evaluation of Epoxy Sealant Systems for Double-Shell Tank Bottom Refurbishment

The performance of epoxy sealants used in the refurbishment of double-shell tank (DST) systems requires balancing processability, thermomechanical stability, and adhesion to cementitious substrates. This study incorporates Heloxy 8 as a reactive diluent into Westlake 862 epoxy to tailor workability and cured-state properties. Rheological time-sweep analysis demonstrates that increasing the diluent content significantly reduces complex viscosity and extends workability, thereby improving pumpability and flow for large-area applications. However, the targeted 2-hour processing window is not fully achieved. Differential scanning calorimetry (DSC) confirms that all formulations cure at room temperature to glass transition temperatures ( T g ) at least 20 °C above the maximum DST operating temperature (27 °C), thereby ensuring service in the glassy regime. Dynamic mechanical analysis (DMA) reveals formulation-dependent reductions in tan delta and increases in storage modulus, indicating increasingly elastic and mechanically stable networks with diluent incorporation. Pull-off adhesion testing shows that modified formulations (70–90% Westlake epoxy) exhibit significantly higher adhesion strengths than the unmodified system. Grout cohesive failure indicates that interfacial bonding exceeds substrate strength. Collectively, these results demonstrate that controlled reactive diluent incorporation enables optimization of processing behavior, interfacial adhesion, and thermomechanical performance, supporting the suitability of the modified epoxy systems as durable sealant layers for cementitious barrier applications in hazardous waste containment infrastructure.

Differential scanning calorimetry↗
DOE OSTI2018

A Structure Function Model Recovers the Many Formulations for Air-Water Gas Transfer Velocity: Air-Water Gas Transfer

Two ideas regarding the structure of turbulence near a clear air-water interface are used to derive a waterside gas transfer velocity kL for sparingly and slightly soluble gases. The first is that k L is proportional to the turnover velocity described by the vertical velocity structure function D ww (r), where r is separation distance between two points. The second is that the scalar exchange between the air-water interface and the waterside turbulence can be suitably described by a length scale proportional to the Batchelor scale lB=ηSc -1/2 , where Sc is the molecular Schmidt number and η is the Kolmogorov microscale defining the smallest scale of turbulent eddies impacted by fluid viscosity. Using an approximate solution to the von Kármán-Howarth equation predicting D ww (r) in the inertial and viscous regimes, prior formulations for k L are recovered including (i) k L =$\sqrt{2/15}$Sc -1/2 $v$, v K is the Kolmogorov velocity defined by the Reynolds number v K η/ν = 1 and ν is the kinematic viscosity of water; (ii) surface divergence formulations; (iii) k L ∝ SC -1/2 u * , where u * is the waterside friction velocity; (iv) k L ∝ Sc -1/2 $\sqrt{gv/u_*}$ for Keulegan numbers exceeding a threshold needed for long-wave generation, where the proportionality constant varies with wave age, g is the gravitational acceleration; and (v) k L = $\sqrt{2/15}$Sc -1/2 (vgβ o q o ) 1/4 in free convection, where qo is the surface heat flux and β o is the thermal expansion of water. The work demonstrates that the aforementioned k L formulations can be recovered from a single structure function model derived for locally homogeneous and isotropic turbulence.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Feasibility of Formulating Ecosystem Biogeochemical Models From Established Physical Rules

Abstract To improve the predictive capability of ecosystem biogeochemical models (EBMs), we discuss the feasibility of formulating biogeochemical processes using physical rules that have underpinned the many successes in computational physics and chemistry. We argue that the currently popular empirically based approaches, such as multiplicative empirical response functions and the law of the minimum, will not lead to EBM formulations that can be continuously refined to incorporate improved mechanistic understanding and empirical observations of biogeochemical processes. Instead, we propose that EBM parameterizations, as a lossy data compression problem, can be better formulated using established physical rules widely used in computational physics and chemistry, and different biogeochemical processes can be more robustly integrated within a reactive‐transport framework. Through several examples, we demonstrate how mathematical representations derived from physical rules can improve understanding of relevant biogeochemical processes and enable more effective communication between modelers, observationalists, and experimentalists regarding essential questions, such as what measurements are needed to meaningfully inform models and how can models generate new process‐level hypotheses to test in empirical studies. Finally, while empirical models with more parameters are often less robust, physical rules‐based models can be more robust and show lower predictive equifinality, stemming from their enhanced consistency in representations of processes, interactions and spatial scaling.

54 ENVIRONMENTAL SCIENCES↗