Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MODEL”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Mixing Cell Model: A One-Dimensional Numerical Model for Assessment of Water Flow and Contaminant Transport in the Unsaturated Zone

This report describes the Mixing Cell Model code, a one dimensional model for water flow and solute transport in the unsaturated zone under steady state or transient flow conditions. The model is based on the principles and assumptions underlying mixing-cell model formulations. The unsaturated zone is discretized into a series of independent mixing cells. Each cell may have unique hydrologic, lithologic, and sorptive properties. Ordinary differential equations describe the material (water and solute) balance within each cell. Water-flow equations are derived from the continuity equation, assuming that unit gradient conditions exist at all times in each cell. Pressure gradients are considered implicitly through model discretization. Unsaturated hydraulic conductivity and moisture contents are determined by the material specific moisture-characteristic curves. Solute-transport processes include explicit treatment of advective processes, first order chain decay, and linear sorption reactions. Dispersion is addressed through implicit and explicit dispersion. Implicit dispersion is an inherent feature of all mixing-cell models and originates from the formulation of the problem in terms of mass balance around fully mixed volume elements. Expressions are provided that relate implicit dispersion to the physical dispersion of the system. Two FORTRAN codes were developed to solve the water flow and solute-transport equations: (1) the Mixing Cell Model for Flow (MCMF) solves transient water-flow problems and (2) the Mixing Cell Model for Transport (MCMT) solves the solute-transport problem. The transient water-flow problem is typically solved first by estimating the water flux through each cell in the model domain as a function of time using the MCMF code. These data are stored in either ASCII or binary files that are later read by the solute transport code (MCMT). Code output includes solute pore water concentrations, water and solute inventories in each cell and at each specified output time, and water and solute fluxes through each cell and specified output time. Computer run times for coupled transient water flow and solute transport were typically several seconds on a 2 GHz Intel Pentium IV desktop computer. The model was benchmarked against analytical solutions and finite element approximations to the partial differential equations (PDE) describing unsaturated flow and transport. Differences between the maximum solute flux estimated by the mixing cell model and the PDE models were typically less than two percent. This revision includes an option for a fixed concentration lower boundary condition for diffusive fluxes for versions 020321 and later.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Validation of Wind Power Plant Modeling Approaches in Complex Terrain: Preprint

The effect of terrain on wind-plant aerodynamics is often neglected due to the complex dynamics that make up terrain-atmosphere-turbine interactions and the computational effort required to adequately resolve the relevant physics. In this work, a comprehensive set of measurements (from the Wind Forecast Improvement Project 2) and two distinct modeling approaches are used to investigate the effects of topography on plant performance. The first approach consists of high-fidelity large-eddy simulations conducted for a one- hour period on a terrain-resolving mesh, and producing unsteady, heterogeneous flow that is sensitive to the local topography. The second approach makes use of a steady-state engineering model that has been modified to account for terrain effects. In situ and remote- sensing measurements upstream of the wind plant are used to prescribe initial and boundary conditions to the models, while wind turbine power measurements are used to assess the ability of the different models to capture the effects of terrain on turbine wake dynamics and therefore power production. The full paper will quantify the degree to which our high-fidelity and engineering models are able to capture realistic wind plant aerodynamics in complex terrain, and through improved understanding of the underlying physics, inform our model development efforts.

17 WIND ENERGY↗
NASA NTRS2023

The Emissions Model Intercomparison Project (Emissions-MIP): Quantifying Model Sensitivity to Emission Characteristics

Anthropogenic emissions of aerosols and precursor compounds are known to significantly affect the energy balance of the Earth–atmosphere system, alter the formation of clouds and precipitation, and have a substantial impact on human health and the environment. Global models are an essential tool for examining the impacts of these emissions. In this study, we examine the sensitivity of model results to the assumed height of SO 2 injection, seasonality of SO 2 and black carbon (BC) particulate emissions, and the assumed fraction of SO 2 emissions that is injected into the atmosphere as particulate phase sulfate (SO 4 ) in 11 climate and chemistry models, including both chemical transport models and the atmospheric component of Earth system models. We find large variation in atmospheric lifetime across models for SO 2 , SO 4 , and BC, with a particularly large relative variation for SO 2 , which indicates that fundamental aspects of atmospheric sulfur chemistry remain uncertain. Of the perturbations examined in this study, the assumed height of SO 2 injection had the largest overall impacts, particularly on global mean net radiative flux (maximum difference of −0.35 W m −2 ), SO 2 lifetime over Northern Hemisphere land (maximum difference of 0.8 d), surface SO 2 concentration (up to 59 % decrease), and surface sulfate concentration (up to 23 % increase). Emitting SO 2 at height consistently increased SO 2 and SO 4 column burdens and shortwave cooling, with varying magnitudes, but had inconsistent effects across models on the sign of the change in implied cloud forcing. The assumed SO 4 emission fraction also had a significant impact on net radiative flux and surface sulfate concentration. Because these properties are not standardized across models this is a source of inter-model diversity typically neglected in model intercomparisons. These results imply a need to ensure that anthropogenic emission injection height and SO 4 emission fraction are accurately and consistently represented in global models.

Emissions Model Intercomparison Project↗
DOE OSTI2025

Investigation of Main Bearing Fatigue Estimate Sensitivity to Synthetic Turbulence Models Using a Novel Drivetrain Model Implemented in OpenFAST

ABSTRACT A coupled medium‐fidelity drivetrain model is developed and implemented in OpenFAST for a 10‐MW land‐based reference turbine. The implementation is verified against a fully coupled multibody wind turbine model, including a detailed drivetrain. The new model can simultaneously and accurately estimate main bearing loads and represent elastic bending of the drivetrain. It has low computational cost and is useful for early design phases, sensitivity analyses and complex systems like wind farms (where computational expense must be expended elsewhere). Here, the model is implemented for a monopile offshore wind turbine and used to investigate the sensitivity of main bearing basic rating life to different synthetic turbulence models. Large‐eddy simulations (LES) targeting stable, neutral, and unstable atmospheric conditions at below‐, near‐ and above‐rated wind speeds are used as a reference. The turbulence models recommended by the International Electrotechnical Commission, the Mann spectral tensor model, and the Kaimal spectral model with exponential coherence are fitted to the LES data. Additionally, a constrained turbulence generator, PyConTurb (short for Python Constrained Turbulence ), based on LES data, is applied in the aero‐hydro‐servo‐elastic simulations. Taking PyConTurb as the baseline, the Kaimal model significantly underestimates fatigue of the downwind main bearing, with between 10% and 40% less damage. The Mann model also underestimates the downwind main bearing fatigue by up to 30%. The upwind main bearing damage is driven by mean loads, and differences between models are less significant, although the trends are similar. Reasons for these discrepancies are investigated and attributed to differences in spatial and temporal variations among the turbulence models.

17 WIND ENERGY↗
DOE OSTI2019

Looking under the hood: A comparison of techno-economic assumptions across national and global integrated assessment models

Integrated assessment models are extensively used in the analysis of climate change mitigation and are informing national decision makers as well as contribute to international assessments such as those of the Intergovernmental Panel on Climate Change. This paper conducts a comprehensive review of techno-economic assumptions in the electricity sector among sixteen different global and national integrated assessment models. Particular focus is given to six major economies in the world: Brazil, China, the EU, India, Japan and the US. The comparison reveals that techno-economic characteristics are quite different across integrated assessment models, both for the base year and future years. It is, however, important to recognize that techno-economic assessments from the literature exhibit an equally large range of parameters as the integrated assessment models reviewed. Beyond numerical differences, the representation of technologies also differs among models which needs to be taken into account when comparing numerical parameters. While desirable from a first principles perspective and often called for, it seems difficult to fully harmonize techno-economic parameters across a broader range of models due to structural differences in the representation of technology. Therefore, in the future making techno-economic parameters available together with of the technology representation as well as the exact definitions of the parameters should become the standard approach as it allows an open discussion of appropriate assumptions.

Krey, Volker↗
DOE OSTI2016

A comparison of simple global kinetic models for coal devolatilization with the CPD model

Simulations of coal combustors and gasifiers generally cannot incorporate the complexities of advanced pyrolysis models, and hence there is interest in evaluating simpler models over ranges of temperature and heating rate that are applicable to the furnace of interest. In this paper, six different simple model forms are compared to predictions made by the Chemical Percolation Devolatilization (CPD) model. The model forms included three modified one-step models, a simple two-step model, and two new modified two-step models. These simple model forms were compared over a wide range of heating rates (5 × 10 3 to 10 6 K/s) at final temperatures up to 1600 K. Comparisons were made of total volatiles yield as a function of temperature, as well as the ultimate volatiles yield. Advantages and disadvantages for each simple model form are discussed. In conclusion, a modified two-step model with distributed activation energies seems to give the best agreement with CPD model predictions (with the fewest tunable parameters).

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2015

Signatures of Planck corrections in a spiralling axion inflation model

The minimal sub-Planckian axion inflation model accounts for a large scalar-to-tensor ratio via a spiralling trajectory in the field space of a complex field Φ. Here we consider how the predictions of the model are modified by Planck scale-suppressed corrections. In the absence of Planck corrections the model is equivalent to a ϕ{sup 4/3} chaotic inflation model. Planck corrections become important when the dimensionless coupling ξ of |Φ|{sup 2} to the topological charge density of the strongly-coupled gauge sector FF{sup ~} satisfies ξ∼1. For values of |Φ| which allow the Planck corrections to be understood via an expansion in powers of |Φ|{sup 2}/M{sub Pl}{sup 2}, we show that their effect is to produce a significant modification of the tensor-to-scalar ratio from its ϕ{sup 4/3} chaotic inflation value without strongly modifying the spectral index. In addition, to leading order in |Φ|{sup 2}/M{sub Pl}{sup 2}, the Planck modifications of n{sub s} and r satisfy a consistency relation, Δn{sub s}=−Δr/16. Observation of these modifications and their correlation would allow the model to be distinguished from a simple ϕ{sup 4/3} chaotic inflation model and would also provide a signature for the influence of leading-order Planck corrections.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

A Process Model for ITCZ Narrowing under Warming Highlights Clear-Sky Water Vapor Feedbacks and Gross Moist Stability Changes in AMIP Models

Tropical areas with mean upward motion—and as such the zonal-mean intertropical convergence zone (ITCZ)—are projected to contract under global warming. To understand this process, a simple model based on dry static energy and moisture equations is introduced for zonally symmetric overturning driven by sea surface temperature (SST). Processes governing ascent area fraction and zonal mean precipitation are examined for insight into Atmospheric Model Intercomparison Project (AMIP) simulations. Bulk parameters governing radiative feedbacks and moist static energy transport in the simple model are estimated from the AMIP ensemble. Uniform warming in the simple model produces ascent area contraction and precipitation intensification—similar to observations and climate models. Contributing effects include stronger water vapor radiative feedbacks, weaker cloud-radiative feedbacks, stronger convection-circulation feedbacks, and greater poleward moisture export. The simple model identifies parameters consequential for the inter-AMIP-model spread; an ensemble generated by perturbing parameters governing shortwave water vapor feedbacks and gross moist stability changes under warming tracks inter-AMIP-model variations with a correlation coefficient ~0.46. Here, the simple model also predicts the multimodel mean changes in tropical ascent area and precipitation with reasonable accuracy. Furthermore, the simple model reproduces relationships among ascent area precipitation, ascent strength, and ascent area fraction observed in AMIP models. A substantial portion of the inter-AMIP-model spread is traced to the spread in how moist static energy and vertical velocity profiles change under warming, which in turn impact the gross moist stability in deep convective regions—highlighting the need for observational constraints on these quantities.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Solar Models in Light of New High Metallicity Measurements from Solar Wind Data

We study the impact of new metallicity measurements, from solar wind data, on the solar model. The “solar modeling problem” refers to the persisting discrepancy between helioseismological observations and predictions of solar models computed implementing state-of-the-art photospheric abundances. We critically reassess the problem, in particular considering the new set of abundances of von Steiger and Zurbuchen, determined through the in situ collection of solar wind samples from polar coronal holes. This new set of abundances indicates a solar metallicity Z{sub ⊙}⩾0.0196±0.0014, significantly higher than the currently established value. The new values hint at an abundance of volatile elements (i.e., C, N, O, Ne) close to previous results of Grevesse and Sauval, whereas the abundance of refractory elements (i.e., Mg, Si, S, Fe) is considerably increased. Using the Linear Solar Model formalism, we determine the variation of helioseismological observables in response to the changes in elemental abundances, in order to explore the consistency of these new measurements with constraints from helioseismology. We find that for observables that are particularly sensitive to the abundance of volatile elements, in particular the radius of the convective zone boundary (CZB) and the sound speed around the radius of CZB, improved agreement over previous models is obtained. Conversely, the high abundance of refractories correlates with a higher core temperature, resulting in an overproduction of neutrinos and a huge increase in the surface helium abundance. We conclude that the “solar modeling problem” remains unsolved.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2018

Chemical Modeling of the Reactivity of Short-Lived Greenhouse Gases: A Model Inter-Comparison Prescribing a Well-Measured, Remote Troposphere

We develop a new protocol for merging in situ measurements with 3-D model simulations of atmospheric chemistry with the goal of integrating over the data to identify the most reactive air parcels in terms of tropospheric production and loss of the greenhouse gases ozone and methane. Presupposing that we can accurately measure atmospheric composition, we examine whether models constrained by such measurements agree on the chemical budgets for ozone and methane. In applying our technique to a synthetic data stream of 14,880 parcels along 180W, we are able to isolate the performance of the photochemical modules operating within their global chemistry-climate and chemistry-transport models, removing the effects of modules controlling tracer transport, emissions, and scavenging. Differences in reactivity across models are driven only by the chemical mechanism and the diurnal cycle of photolysis rates, which are driven in turn by temperature, water vapor, solar zenith angle, clouds, and possibly aerosols and overhead ozone, which are calculated in each model. We evaluate six global models and identify their differences and similarities in simulating the chemistry through a range of innovative diagnostics. All models agree that the more highly reactive parcels dominate the chemistry (e.g., the hottest 10% of parcels control 25-30% of the total reactivities), but do not fully agree on which parcels comprise the top 10%. Distinct differences in specific features occur, including the regions of maximum ozone production and methane loss, as well as in the relationship between photolysis and these reactivities. Unique, possibly aberrant, features are identified for each model, providing a benchmark for photochemical module development. Among the 6 models tested here, 3 are almost indistinguishable based on the inherent variability caused by clouds, and thus we identify 4, effectively distinct, chemical models. Based on this work, we suggest that water vapor differences in model simulations of past and future atmospheres may be a cause of the different evolution of tropospheric O3 and CH4, and lead to different chemistry-climate feedbacks across the models.

greenhouse gases↗
NASA NTRSDate not supplied

Aeroelastic Modeling and CFD Simulation of Wind-Tunnel Scale Aspect Ratio 13.5 Common Research Model

This paper presents an aeroelastic modeling and simulation study of an aspect ratio 13.5 wind-tunnel scale Common Research Model (CRM) with distributed flaps. A vortex-lattice VSPAERO model of the CRM model is developed. A transonic small disturbance/integral boundary layer correction method is implemented in the VSPAERO model to account for the transonic and viscous flow effects. The structural deformation of the CRM model is calculated using a NASTRAN equivalent beam model. The VSPAERO model is coupled to the NASTRAN equivalent beam model to provide a rapid aero-structural analysis. A validation of the VSPAERO aeroelastic model is conducted by comparing the results to FUN3D CFD aeroelastic simulation results. An aerodynamic database is generated using the developed VSPAERO aeroelastic model for the real-time drag optimization and maneuver load alleviation study of the wind-tunnel scale CRM model.

Aeroelastic Modeling↗
DOE OSTI2026

Machine learning-enhanced hybrid modeling approach for better identification of a building thermal network model and improved prediction

The gray-box modeling approach, which uses a semi-physical thermal network model, has been widely used in building prediction applications, such as model predictive control (MPC). However, unmeasured disturbances, such as occupants, lighting, and in/exfiltration loads, make it challenging to apply this approach to practical buildings. In this word, we propose a hybrid modeling approach that integrates the gray-box model with a model for unmeasured disturbance. After reviewing several system identification approaches, we systematically designed the unmeasured disturbance model with a model selection process based on statistical tests to make it robust. We generated data based on the building model calibrated by real operational data and then trained the hybrid model for two different weather conditions. The hybrid model approach demonstrates an RMSE reduction of approximately 0.2–0.9 °C and 0.3–2 °C on 1-day ahead temperature prediction compared to the Conventional approach for mild (Berkeley, CA) and cold (Chicago, IL) climates, respectively. In addition, this approach was applied to experimental data obtained from the laboratory building to be used for the MPC application, showing superior prediction performances.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Model America: Data and Models for every U.S. Building

The 5-year goal of the “Model America” concept was to generate a model of every building in the United States. This data repository delivers on that goal with "Model America v1". Oak Ridge National Laboratory (ORNL) has developed the Automatic Building Energy Modeling (AutoBEM) software suite to process multiple types of data, extract building-specific descriptors, generate building energy models, and simulate them on High Performance Computing (HPC) resources. For more information, see AutoBEM-related publications (bit.ly/AutoBEM). There were 125,715,609 buildings detected in the United States. Of this number, 122,146,671 (97.2%) buildings resulted in a successful generation and simulation of a building energy model. This dataset includes the full 125 million buildings. Future updates may include additional buildings, data improvements, or other algorithmic model enhancements in "Model America v2". This dataset contains OSM and IDF zip files for every U.S. county. Each zip file contains the generated buildings from that county. The .csv input data contains the following data fields: 1. ID - the Unique Building Identifier (UBID), generated using the Pacific Northwest National Laboratory (PNNL) BuildingID framework 2. Centroid - building center location in latitude/longitude (from Footprint2D) 3. Footprint2D - building polygon of 2D footprint (lat1/lon1_lat2/lon2_...) 4. State_abbr - state name 5. Area - estimate of total conditioned floor area (ft2) 6. Area2D - footprint area (ft2) 7. Height - building height (ft) 8. NumFloors - number of floors (above-grade) 9. WWR_surfaces - percent of each facade (pair of points from Footprint2D) covered by fenestration/windows (average 14.5% for residential, 40% for commercial buildings) 10. CZ - ASHRAE Climate Zone designation 11. BuildingType - DOE prototype building designation (IECC=residential) as implemented by OpenStudio-standards 12. Standard - building vintage This data is made free and openly available in hopes of stimulating any simulation-informed use case. Data is provided as-is with no warranties, express or implied, regarding fitness for a particular purpose. We wish to thank our sponsors which include Oak Ridge National Laboratory (ORNL) Laboratory Directed Research and Development (LDRD), U.S. Dept. of Energy’s (DOE) Building Technologies Office (BTO), Office of Electricity (OE), Biological and Environmental Research (BER), and National Nuclear Security Administration (NNSA). Update (September 23, 2025): We corrected the ID field in all state-level.csv input files to ensure one-to-one consistency with the corresponding .osm and .idf output files. The schema and file structure are unchanged; only the values in the ID column were modified. No files were added or removed, and the .zip bundles (containing .osm / .idf) are unchanged. The corrected .csv inputs were re-extracted in March 2025 from the original data generated ~ 2021 (Theta supercomputer runs), and published here to align input IDs with model outputs. Update (September 6, 2026): The Model America dataset was updated to replace the previous building ID field with the Unique Building Identifier (UBID), using the Pacific Northwest National Laboratory (PNNL) BuildingID framework. UBIDs provide standardized, location-based identifiers for individual building footprints and improve interoperability with other building and geospatial datasets. The data files containing the previous building identifiers were updated to include UBIDs. This update standardizes building identification; the underlying Model America building characteristics and energy simulation results were not recomputed as part of this update.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Concurrent Inter-Model Spread of Boreal Winter Westerly Jet Meridional Positions Between the Northern and Southern Hemispheres in CMIP6 Models

Here, this study investigates the inter-model spread of climatological extratropical westerly jets in boreal winter, using the historical simulation of 52 Coupled Model Intercomparison Project phase 6 (CMIP6) models from 1851 to 2014. The results show that there is a substantial spread in the latitude of the upper-tropospheric westerly jet across models, characterised by large inter-model standard deviations to both the poleward and equatorward sides of the jet axis, although the multi-model ensemble mean (MME) performs well in simulating meridional position of westerly jets. Furthermore, we detect the consistency of inter-model jet position spread between the Northern and Southern Hemispheres, based on the inter-model empirical orthogonal function (EOF) decomposition and correlation of regional-averaged zonal winds. Specifically, the models that simulate the westerly jets poleward/equatorward relative to the MME position in one hemisphere also tend to simulate the jets poleward/equatorward in the other hemisphere. Accordingly, we define a global jet spread index to depict the concurrence of jet shift in the two hemispheres. The results of inter-model regression analyses based on this index indicate that the models positioning the jets poleward than the MME tend to simulate a wider Hadley Cell, a poleward-shifted Ferrel Cell in the Southern Hemisphere, enhanced precipitation in the subtropics and suppressed precipitation in the tropics, and warmer sea surface temperatures in the subtropics and mid-latitudes. The present results suggest that improving the simulation of jet positions in climate models requires a comprehensive consideration of thermal states in the tropics and subtropics/mid latitudes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

High-Performance Computing for Earth System Modeling

High-performance computing (HPC) plays an important role during the development of Earth system models. This chapter reviews HPC efforts related to Earth system models, including community Earth system models and energy exascale Earth system models. Specifically, this chapter evaluates computational and software design issues, analyzes several current HPC-related model developments, and provides an outlook for some promising areas within Earth system modeling in the era of exascale computing.

Wang, Dali↗
DOE OSTI2024

Implementation of an extensible property modeling framework in ESPEI with applications to molar volume and elastic stiffness models

Property models are becoming more widely adopted by commercial Calphad databases, but they are not nearly as common in non-commercial or traditional academic Calphad databases. A primary driver is that user-friendly Calphad modeling tools that support property models are not widely available. Here we present new property modeling capabilities that have been implemented in ESPEI (the Extensible, Self-optimizing Phase Equilibrium Infrastructure). These capabilities include both generating property model parameters from data and improvements to the algorithmic selection of the most appropriate model from a series of candidates. Additionally, two illustrative examples are given that use ESPEI to fit different property models. First, we generate molar volume model parameters for Group IV, V, and VI refractory BCC alloys based on the model by Lu et al. (2005). Second, we demonstrate the extensibility of ESPEI’s property modeling capabilities by implementing a custom PyCalphad model for BCC elastic stiffness parameters to generate and compare parameters to the ones assessed by Marker et al. (2018) using the same data. Property models generated by ESPEI can be used in PyCalphad or further optimized with uncertainty quantification using ESPEI.

36 MATERIALS SCIENCE↗
DOE OSTI2019

A two-set order parameters phase-field modeling of crack deflection/penetration in a heterogeneous microstructure

Modeling of crack deflection/penetration or manifested as intergranular/transgranular fracture in polycrystals has long been a challenge in both fracture mechanics and materials science. A phase-field model formulated with two-set order parameters describing the crack field and the microstructure field respectively, is herein established to investigate the competition between crack penetration and deflection at an interface. The advantages of the proposed model are three-fold: (1) the fracture toughness of grain boundaries and/or interfaces, besides single crystal properties (e.g., Young’s modulus, Poisson ratio, and fracture toughness), can be automatically incorporated, (2) cohesive interfaces are treated in one phase-field framework in comparison with the integrated phase-field and cohesive zone model in which cohesive elements have to be manually placed at the interface, and (3) two-set order parameters enable to describe the concurrent crack propagation and microstructure evolution. The proposed model is validated by comparing the simulated critical deflection angles with the results from theoretical analysis of the linear elastic fracture mechanics and the experiments. The model is then applied to a heterogeneous Silicon Nitride polycrystal to determine the effect of interface fracture toughness on the intergranular/transgranular fracture transition and the effective toughness of the materials.

Chen, Hu↗