Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model parameterization”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

An Evaluation of Marine Boundary Layer Cloud Property Simulations in the Community Atmosphere Model Using Satellite Observations: Conventional Subgrid Parameterization versus CLUBB

This paper presents a satellite-observation-based evaluation of the marine boundary layer (MBL) cloud properties from two Community Atmosphere Model, version 5 (CAM5), simulations, one with the standard parameterization schemes (CAM5–Base) and the other with the Cloud Layers Unified by Binormals scheme (CAM5–CLUBB). When comparing the direct model outputs, the authors find that CAM5–CLUBB produces more MBL clouds, a smoother transition from stratocumulus to cumulus, and a tighter correlation between in-cloud water and cloud fraction than CAM5–Base. In the model-to-observation comparison using the COSP satellite simulators, the authors find that both simulations capture the main features and spatial patterns of the observed cloud fraction from MODIS and shortwave cloud radiative forcing (SWCF) from CERES. However, CAM5–CLUBB suffers more than CAM5–Base from a problem that can be best summarized as “undetectable” clouds (i.e., a significant fraction of simulated MBL clouds are thinner than the MODIS detection threshold). This issue leads to a smaller COSP–MODIS cloud fraction and a weaker SWCF in CAM5–CLUBB than the observations and also CAM5–Base in the tropical descending regions. Finally, the authors compare modeled radar reflectivity with CloudSat observations and find that both simulations, especially CAM5–CLUBB, suffer from an excessive drizzle problem. Further analysis reveals that the subgrid precipitation enhancement factors in CAM5–CLUBB are unrealistically large, which makes MBL clouds precipitate too excessively, and in turn results in too many undetectable thin clouds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Phase-field modeling of stored-energy-driven grain growth with intra-granular variation in dislocation density

Abstract We present a phase-field (PF) model to simulate the microstructure evolution occurring in polycrystalline materials with a variation in the intra-granular dislocation density. The model accounts for two mechanisms that lead to the grain boundary migration: the driving force due to capillarity and that due to the stored energy arising from a spatially varying dislocation density. In addition to the order parameters that distinguish regions occupied by different grains, we introduce dislocation density fields that describe spatial variation of the dislocation density. We assume that the dislocation density decays as a function of the distance the grain boundary has migrated. To demonstrate and parameterize the model, we simulate microstructure evolution in two dimensions, for which the initial microstructure is based on real-time experimental data. Additionally, we applied the model to study the effect of a cyclic heat treatment (CHT) on the microstructure evolution. Specifically, we simulated stored-energy-driven grain growth during three thermal cycles, as well as grain growth without stored energy that serves as a baseline for comparison. We showed that the microstructure evolution proceeded much faster when the stored energy was considered. A non-self-similar evolution was observed in this case, while a nearly self-similar evolution was found when the microstructure evolution is driven solely by capillarity. These results suggest a possible mechanism for the initiation of abnormal grain growth during CHT. Finally, we demonstrate an integrated experimental-computational workflow that utilizes the experimental measurements to inform the PF model and its parameterization, which provides a foundation for the development of future simulation tools capable of quantitative prediction of microstructure evolution during non-isothermal heat treatment.

Materials Science↗
DOE OSTI2015

A review on the relativistic effective field theory with parameterized couplings for nuclear matter and neutron stars

Nuclear science has developed many excellent theoretical models for many-body systems in the domain of the baryon-meson strong interaction for the nucleus and nuclear matter at low, medium and high densities. However, a full microscopic understanding of nuclear systems in the extreme density domain of compact stars is still lacking. The aim of this contribution is to shed some light on open questions facing the nuclear many-body problem at the very high density domain. Here we focus our attention on the conceptual issue of naturalness and its role in shaping the baryon-meson phase space dynamics in the description of the equation of state (EoS) of nuclear matter and neutrons stars. In particular, in order to stimulate possible new directions of research, we discuss relevant aspects of a recently developed relativistic effective theory for nuclear matter within Quantum Hadrodynamics (QHD) with genuine many-body forces and derivative natural parametric couplings. Among other topics we discuss in this work the connection of this theory with other known effective QHD models of the literature and its potentiality in describing a new physics for dense matter. The model with parameterized couplings exhausts the whole fundamental baryon octet (n, p, Σ{sup −}, Σ{sup 0}, Σ{sup +}, Λ, Ξ{sup −}, Ξ{sup 0}) and simulates n-order corrections to the minimal Yukawa baryon couplings by considering nonlinear self-couplings of meson fields and meson-meson interaction terms coupled to the baryon fields involving scalar-isoscalar (σ, σ∗), vector-isoscalar (ω, Φ), vector-isovector (ϱ) and scalar-isovector (δ) virtual sectors. Following recent experimental results, we consider in our calculations the extreme case where the Σ{sup −} experiences such a strong repulsion that its influence in the nuclear structure of a neutron star is excluded at all. A few examples of calculations of properties of neutron stars are shown and prospects for the future are discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Effects of random forest modeling decisions on biogeochemical time series predictions

Abstract Random forests (RF) are an increasingly popular machine learning approach used to model biogeochemical processes in the Earth system. While RF models are robust to many assumptions that complicate deterministic models, there are several important parameterization decisions for appropriate use and optimal model fit. We explored the role that parameter decisions, including training/testing data splitting strategies, variable selection, and hyperparameters play on RF goodness‐of‐fit by constructing models using 1296 unique parameter combinations to predict concentrations of nitrate, a key nutrient for biogeochemical cycling in aquatic ecosystems. Models were built on long‐term, publicly available water quality and meteorology time series collected by the National Estuarine Research Reserve monitoring network for two contrasting ecosystems representing freshwater and brackish estuaries. We found that accounting for temporal dependence when splitting data into training and testing subsets was key for avoiding over‐estimation of model predictive power. In addition, variable selection, the ratio of training to testing data, and to a lesser degree, variables per split and number of trees, were significant parameters for optimizing RF goodness‐of‐fit. We also explored how model parameter decisions influenced interpretation of the relative importance of predictors to the model, and model predictor‐dependent variable relationships, with results suggesting that both data structure and model parameterization influence these factors. Because much of the current RF literature is written for the computational and statistical science communities, the primary goal of this study is to provide guidelines for aquatic scientists new to machine learning to apply RF techniques appropriately to aquatic biogeochemical datasets.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Polynomial Chaos–Based Bayesian Inference of K -Profile Parameterization in a General Circulation Model of the Tropical Pacific

The authors present a polynomial chaos (PC)–based Bayesian inference method for quantifying the uncertainties of the K-profile parameterization (KPP) within the MIT general circulation model (MITgcm) of the tropical Pacific. The inference of the uncertain parameters is based on a Markov chain Monte Carlo (MCMC) scheme that utilizes a newly formulated test statistic taking into account the different components representing the structures of turbulent mixing on both daily and seasonal time scales in addition to the data quality, and filters for the effects of parameter perturbations over those as a result of changes in the wind. To avoid the prohibitive computational cost of integrating the MITgcm model at each MCMC iteration, a surrogate model for the test statistic using the PC method is built. Because of the noise in the model predictions, a basis-pursuit-denoising (BPDN) compressed sensing approach is employed to determine the PC coefficients of a representative surrogate model. The PC surrogate is then used to evaluate the test statistic in the MCMC step for sampling the posterior of the uncertain parameters. Results of the posteriors indicate good agreement with the default values for two parameters of the KPP model, namely the critical bulk and gradient Richardson numbers; while the posteriors of the remaining parameters were barely informative.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Molecular thermodynamics for scaling prediction: Case of membrane distillation

Membrane distillation (MD) for water treatment is significantly impaired by the scaling of dissolved minerals. The type and content of minerals generally measured as total dissolved solids (TDS) in hypersaline brines not only reduce the MD flux but also control the scaling behavior on the membrane surface. The scaling-induced pore blockage further reduces water flux and eventually leads to membrane wetting. The scaling problem is even more pronounced in the treatment of produced water (PW) as it contains 3–7 times higher TDS concentrations, compared to seawater. Theoretically, the necessary conditions for a salt to precipitate can be traced from its solubility product constant and activity of the constituents within the solution. Therefore, a comprehensive thermodynamic model is necessary to represent the electrolyte behavior and to predict the precipitation of different salts in a complex solution like PW. We pursued electrolyte Nonrandom Two Liquid Theory (eNRTL), a state-of-the-art electrolyte model, to address the PW fluid phase equilibria. With a fully parameterized eNRTL model, we predicted salt precipitations in two different PW samples and compared the results against the experimental findings. Furthermore, here we forecast the precipitation of salts in response to the change in PW concentration and temperature in the MD operation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Characterizing the Influence of Fracture Density on Network Scale Transport

The topology of natural fracture networks is inherently linked to the structure of the fluid velocity field and transport therein. Here we study the impact of network density on flow and transport behaviors. We stochastically generate fracture networks of varying density and simulate flow and transport with a discrete fracture network model, which fully resolves network topology at the fracture scale. We study conservative solute trajectories with Lagrangian particle tracking and find that as fracture density decreases, solute channelization to large local fractures increases, thereby reducing plume spreading. Furthermore, in sparse networks mean particle travel distance increases and local network features, such as velocity zones where flow is counter to the primary pressure gradient, become increasingly important for transport. As the network density increases, network statistics homogenize and such local features have a reduced impact. We quantify local topological influence on transport behavior with an effective tortuosity parameter, which measures the ratio of total advective distance to linear distance at the fracture scale; large tortuosity values are correlated to slow-velocity regions. Additionally, these large tortuosity, slow-velocity regions delay downstream transport and enhance tailing on particle breakthrough curves. Finally, we predict transport with an upscaled, Bernoulli spatial Markov random walk model and parameterize local topological influences with a novel tortuosity parameter. Furthermore, Bernoulli model predictions improve when sampling from a tortuosity distribution, as opposed to a fixed value as has previously been done, suggesting that local network topological features must be carefully considered in upscaled modeling efforts of fracture network systems.

58 GEOSCIENCES↗
DOE OSTI2025

Dark Energy Survey: Modeling strategy for multiprobe cluster cosmology and validation for the Full Six-year Dataset

We introduce an updated To&Krause2021 model for joint analyses of cluster abundances and large-scale two-point correlations of weak lensing and galaxy and cluster clustering (termed CL+3x2pt analysis) and validate that this model meets the systematic accuracy requirements of analyses with the statistical precision of the final Dark Energy Survey (DES) Year 6 (Y6) dataset. The validation program consists of two distinct approaches, (1) identification of modeling and parameterization choices and impact studies using simulated analyses with each possible model misspecification (2) end-to-end validation using mock catalogs from customized Cardinal simulations that incorporate realistic galaxy populations and DES-Y6-specific galaxy and cluster selection and photometric redshift modeling, which are the key observational systematics. In combination, these validation tests indicate that the model presented here meets the accuracy requirements of DES-Y6 for CL+3x2pt based on a large list of tests for known systematics. In addition, we also validate that the model is sufficient for several other data combinations: the CL+GC subset of this data vector (excluding galaxy--galaxy lensing and cosmic shear two-point statistics) and the CL+3x2pt+BAO+SN (combination of CL+3x2pt with the previously published Y6 DES baryonic acoustic oscillation and Y5 supernovae data).

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2017

Final Technical Report for Award SC0008613

Discovering how aerosol particles, present in the atmosphere in sizes of a few nanometers to hundred micrometers, initiate ice crystal formation represents a great challenge. Atmospheric ice nucleation is important because ice crystals alter the radiative properties of clouds and thus climate, and impact precipitation and thus the hydrological cycle. The difficulty in predicting atmospheric ice formation is attributable at least in part, to the diversity of ice nucleation pathways, the physical and chemical complexity of the ice nucleating particles (INPs), and the relatively small numbers of INPs (compared with all other aerosol particles), sometimes less than one in 100000. These factors in turn makes constraining ice nucleation parameterizations for modeling applications challenging. The majority of airborne particles are known to be organic in nature or contain organic biogenic material. The presence of organic material adds to the complexity of the particles and therefore the predictability of ice nucleation events since the organic species can display different phase states, e.g. liquid or solid, in response to temperature and humidity. The award DE-SC0008613 to PI Prof. Daniel Knopf at Stony Brook University, “Relating the Chemical and Physical Properties of Aerosols to the Water Uptake and Ice Nucleation Potential of Particles Collected During the Carbonaceous Aerosols and Radiative Effects Study (CARES)”, allowed examination of laboratory generated aerosol particles and field-collected particles for their propensity to nucleate ice under typical tropospheric conditions and relate ice nucleation to the physicochemical properties of the particles including their morphology and chemical composition. This in turn allowed for development of ice nucleation parameterizations for implementation in cloud models. The award resulted in 10 peer-reviewed publications and more than 20 seminar and conference presentations. We demonstrated that the rate of immersion freezing, when a particle immersed in an aqueous droplet acts as an ice nucleus (IN), can be predicted by knowledge of the IN particle type present and the droplet’s water activity which is equal to ambient relative humidity. Our water activity based immersion freezing model is successful in predicting freezing data including INPs such as mineral dusts, marine biological material, organic species, and surfactant molecules. Its mathematical simplicity makes it an ideal candidate for implementation in cloud and climate models. Furthermore, we could show that this model can reproduce many past laboratory measurements which were generated using a variety of instruments to study immersion freezing. Lastly, we have demonstrated that this model can also be applied to field collected particles. Currently, this novel physical parameterization of immersion freezing is being implemented in a cloud model. We examined the physicochemical properties and the ice nucleation potential of particles collected during CARES applying a novel experimental method that allows identification of the individual INPs within a large population of particles sampled from an ambient environment. Taking advantage of a variety of micro-spectroscopic techniques, we characterized the composition and morphology of IN and non-IN particles present in the airborne population. We developed a new parameterization for quantifying the mixing state of the entire aerosol populations by introducing a mixing state index. We found that the identified INPs belong to the most common particle-type classes observed in the CARES field samples and as such are not special or rare particles. In other words, the INPs can be shown not to be unique in contrast to the common paradigm of being rare and exceptional. Either there are differences between particles acting as IN and particles not acting as IN which are beyond our current detection limit or nucleation occurs randomly on the surface of any one of these compositionally equivalent particles. These results suggest that total particle surface area of the different particle types present in the aerosol population is also a crucial factor when predicting ice nucleation in an air mass. We also observed that ambient organic aerosol particles can initiate ice nucleation and corroborated these findings using laboratory generated organic particles. These experiments demonstrated that information of the organic phase state is important for predicting the ice nucleation pathway and the ability of an organic particle to participate in atmospheric ice crystal formation. This award resulted in substantial new insights in the processes governing immersion freezing, the role of organic aerosol particles in ice cloud formation, and the importance of the ambient aerosol population for prediction of ice nucleation in an air parcel. These findings have significant implications for modeling and field measurement strategies of atmospheric ice nucleation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Development of Heterogeneous Ice Nucleation Rate Coefficient Parameterizations From Ambient Measurements

Ice nucleating particles (INPs) are a rare subset of particles that can have an outsized impact relative to their prevalence. To simulate INP variability, models require parameterizations for the most important sources of INPs. Most parameterizations in the literature were developed from laboratory experiments that used proxies for ambient particles. Whether these laboratory models effectively represent ambient aerosol is still not well understood and parameterizations developed from field measurements may be more atmospherically relevant. Expanding on previous work in which we used an ice chamber and a single-particle mass spectrometer (SPMS) to characterize INP composition, we calculate heterogeneous ice nucleation rate coefficients for ambient particles. We find good agreement between our dust and previous measurements of airborne dust. We also use a Monte Carlo approach to assess the relative contribution of each measurement to the uncertainty and find that the biggest source is the sampling efficiency of the SPMS.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Dataset for "Development of heterogeneous ice nucleation rate coefficient parameterizations from ambient measurements"

Ice nucleating particles (INPs) are a rare subset of particles that can have an outsized impact relative to their prevalence. To simulate INP variability, models require parameterizations for the most important sources of INPs. Most parameterizations in the literature were developed from laboratory experiments that used proxies for ambient particles. Whether these laboratory models effectively represent ambient aerosol is still not well understood and parameterizations developed from field measurements may be more atmospherically relevant. Expanding on previous work in which we used an ice chamber and an single particle mass spectrometer (SPMS) to characterize INP composition, we use the combined measurements to calculate heterogeneous ice nucleation rate coefficients for ambient particles. We find good agreement between our results for dust and sea spray aerosol with previous results. We also determine the relative contribution of each measurement to the total measurement uncertainty using a Monte Carlo approach and find that the biggest source of uncertainty is the sampling efficiency of the SPMS.

Cornwell, Gavin C↗
DOE OSTI2021

Wind Energy Forecasting with the Weather Research and Forecasting Model

This was a collaborative effort between Lawrence Livermore National Security, LLC as manager and operator of Lawrence Livermore National Laboratory (LLNL) and Siemens Energy, Inc. (Siemens) to develop a wind resource forecasting tool. LLNL was to develop an independent high-resolution mesoscale modeling capability forecasting tool that could be implemented in conjunction with existing wind farm control and monitoring software to provide forecasting of wind resources using local observations of winds and temperature. Research with LLNL’s state-of-the-art large-eddy simulation meteorological prediction model, based on the community WRF model and innovative turbulence parameterizations, would improve that model’s applicability to large wind farms offshore and in complex terrain. The modeling capability would include uncertainty quantification. Finally, the application of the modeling tool and existing global climate change predictions would enable the delineation of the likely effects of climate change on wind resources. Siemens was to provide high time resolution hub-height wind speed and other meteorological data streams, including temperature profiles from wind farms, for LLNL to incorporate into the modeling system, to validate and tune this forecasting model for their locations of interest. These data streams would also be used for longer-term studies of correlations of wind resources to climate oscillations to indicate how long-term climate change trends may affect the available wind resource. Siemens would also provide information and observations of turbine wakes for incorporation into the modeling tool. By implementing state-of-the-art turbulence parameterizations into a simulation model and/or ensembles of simulation models, and by integrating real-time hub height wind speed and other meteorological datastreams from wind farms into that model or ensemble of models, LLNL would develop a forecasting tool that could be implemented by Siemens as an add-on to existing wind farm control and monitoring software to provide owners with useful resource forecasting. The desired outcome was that the accuracy level of the output would be sufficient to substantiate power output commitments. The final deliverable for this work would consist of a document outlining the algorithms and software tools that could be integrated into Siemens Wind Park Supervisor.

17 WIND ENERGY↗
DOE OSTI2015

Evaluating convective parameterization closures using cloud‐resolving model simulation of tropical deep convection

Abstract Closure is an important component of a mass flux‐based convective parameterization scheme, and it determines the amount of convection with the aid of a large‐scale variable (closure variable) that is sensitive to convection. In this study, we have evaluated and quantified the relationship between commonly used closure variables and convection for a range of global climate model (GCM) horizontal resolutions, taking convective precipitation and mass flux at 600 hPa as measures for deep convection. We have used cloud‐resolving model simulation data to create domain averages representing GCM horizontal resolutions of 128 km, 64 km, 32 km, 16 km, 8 km, and 4 km. Lead‐lag correlation analysis shows that except moisture convergence and turbulent kinetic energy, none of the other closure variables evaluated in this study show any relationship with convection for the six subdomain sizes. It is found that the correlation between moisture convergence and convective precipitation is largest when moisture convergence leads convection. This correlation weakens as the subdomain size decreases to 8 km or smaller. Although convective precipitation and mass flux increase with moisture convergence at a given subdomain size, as the subdomain size increases, the rate at which they increase becomes smaller. This suggests that moisture convergence‐based closure should scale down the predicted mass flux for a given moisture convergence as GCM resolution increases.

Suhas, E.↗
DOE OSTI2017

Simultaneous event-specific estimates of transport, loss, and source rates for relativistic outer radiation belt electrons: Event-Specific 1-D Modeling

The most significant unknown regarding relativistic electrons in Earth’s outer Van Allen radiation belt is the relative contribution of loss, transport, and acceleration processes within the inner magnetosphere. Detangling each individual process is critical to improve the understanding of radiation belt dynamics, but determining a single component is challenging due to sparse measurements in diverse spatial and temporal regimes. However, there are currently an unprecedented number of spacecraft taking measurements that sample different regions of the inner magnetosphere. With the increasing number of varied observational platforms, system dynamics can begin to be unraveled. In this work, we employ in-situ measurements during the 13-14 January 2013 enhancement event to isolate transport, loss, and source dynamics in a one dimensional radial diffusion model. We then validate the results by comparing them to Van Allen Probes and THEMIS observations, indicating that the three terms have been accurately and individually quantified for the event. Finally, a direct comparison is performed between the model containing event-specific terms and various models containing terms parameterized by geomagnetic index. Models using a simple 3/Kp loss timescale show deviation from the event specific model of nearly two orders of magnitude within 72 hours of the enhancement event. However, models using alternative loss timescales closely resemble the event specific model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Benchmark Modeling of the Near-Field and Far-Field Wave Effects of Wave Energy Arrays

This project is an industry-led partnership between Columbia Power Technologies and Oregon State University that will perform benchmark laboratory experiments and numerical modeling of the near-field and far-field impacts of wave scattering from an array of wave energy devices. These benchmark experimental observations will help to fill a gaping hole in our present knowledge of the near-field effects of multiple, floating wave energy converters and are a critical requirement for estimating the potential far-field environmental effects of wave energy arrays. The experiments will be performed at the Hinsdale Wave Research Laboratory (Oregon State University) and will utilize an array of newly developed Buoys' that are realistic, lab-scale floating power converters. The array of Buoys will be subjected to realistic, directional wave forcing (1:33 scale) that will approximate the expected conditions (waves and water depths) to be found off the Central Oregon Coast. Experimental observations will include comprehensive in-situ wave and current measurements as well as a suite of novel optical measurements. These new optical capabilities will include imaging of the 3D wave scattering using a binocular stereo camera system, as well as 3D device motion tracking using a newly acquired LED system. These observing systems will capture the 3D motion history of individual Buoys as well as resolve the 3D scattered wave field; thus resolving the constructive and destructive wave interference patterns produced by the array at high resolution. These data combined with the device motion tracking will provide necessary information for array design in order to balance array performance with the mitigation of far-field impacts. As a benchmark data set, these data will be an important resource for testing of models for wave/buoy interactions, buoy performance, and far-field effects on wave and current patterns due to the presence of arrays. Under the proposed project we will initiate high-resolution (fine scale, very near-field) fluid/structure interaction simulations of buoy motions, as well as array-scale, phase-resolving wave scattering simulations. These modeling efforts will utilize state-of-the-art research quality models, which have not yet been brought to bear on this complex problem of large array wave/structure interaction problem.

16 TIDAL AND WAVE POWER↗
DOE OSTI2019

Radiation-induced Swelling and Microcracking in SiC Cladding for LWRs

After the Fukushima accident in 2011, SiC has attracted considerable attention as a potential material for fuel cladding that can offer increased accident tolerance in light water reactors (LWRs) compared to zircaloy cladding. Existing measurements on nuclear grade SiC/SiC composites made of high-purity, stoichiometric SiC have already shown very good radiation tolerance in various harsh environments. However, there are a number of challenges that need to be addressed before SiC composites can be successfully used in LWR technologies. One of them is irradiation-induced swelling, which can lead to development of significant stresses in the matrix, followed by microcracking and subsequent release of fission products. Models of irradiation-induced swelling had been largely empirical and therefore they cannot be used to predict the microstructural dependence of microcracking. Existing models are also limited by the lack of understanding of defects that form due to irradiation in the range of temperatures relevant to fuel cladding in LWRs (<1000°C). Many of the defects in this regime of temperature are too small to be detected with traditional transmission electron microscopy (TEM) techniques. In fact the defects observed in traditional TEM account only for 10-45% of the swelling measured in irradiated SiC. Here, we developed a multi-scale simulation methodology combined with state of the art experimental imaging techniques based on aberration corrected scanning transmission electron microscopy (STEM) to provide a scientific basis for predictive models of radiation swelling and microcracking in SiC. Specifically, we have developed a phase field fracture model that predicts microcrack initiation and growth in polycrystalline SiC. The model accounts for single crystal anisotropy, including anisotropy in the elastic constants and anisotropic surface energy. It also considers the reduced fracture strength of grain boundaries. The model was implemented using the MARMOT mesoscale nuclear materials tool. We have parameterized the model for SiC using molecular dynamics simulation results. The model has been verified using single crystal, bicrystal, and polycrystal examples. Swelling model has also been developed based on cluster dynamics simulations of defect evolution in SiC. We found that point defects and small defect clusters have a significant contribution to swelling. This finding is important because point defects and very small defect clusters are not visible in experiments (e.g., in TEM). Combined modeling and experiments can bring critical insights into how swelling depends on irradiation conditions and the microstructure. In situ nanoindentation tests and ex situ nanoindentation tests were performed on assynthesized 3C-SiC and 3C-SiC irradiated under different radiation conditions. Without radiation, bend contour movement and residual bend contours after indentation demonstrate some room temperature plastic flow and residual plastic deformation. Plastic deformation recovery is also observed on unirradiated samples during unloading. Radiation-induced embrittlement is demonstrated by differences in crack geometry in in situ tests, and is consistent with ex situ test results. Failure by single, straight, cracks indicates decreased fracture toughness due to irradiation in 3C-SiC. Future work includes coupling of the swelling and fracture models and development of a master equation that can be used in the BISON code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2015

Complex functionality with minimal computation: Promise and pitfalls of reduced‐tracer ocean biogeochemistry models

Abstract Earth System Models increasingly include ocean biogeochemistry models in order to predict changes in ocean carbon storage, hypoxia, and biological productivity under climate change. However, state‐of‐the‐art ocean biogeochemical models include many advected tracers, that significantly increase the computational resources required, forcing a trade‐off with spatial resolution. Here, we compare a state‐of‐the art model with 30 prognostic tracers (TOPAZ) with two reduced‐tracer models, one with 6 tracers (BLING), and the other with 3 tracers (miniBLING). The reduced‐tracer models employ parameterized, implicit biological functions, which nonetheless capture many of the most important processes resolved by TOPAZ. All three are embedded in the same coupled climate model. Despite the large difference in tracer number, the absence of tracers for living organic matter is shown to have a minimal impact on the transport of nutrient elements, and the three models produce similar mean annual preindustrial distributions of macronutrients, oxygen, and carbon. Significant differences do exist among the models, in particular the seasonal cycle of biomass and export production, but it does not appear that these are necessary consequences of the reduced tracer number. With increasing CO 2 , changes in dissolved oxygen and anthropogenic carbon uptake are very similar across the different models. Thus, while the reduced‐tracer models do not explicitly resolve the diversity and internal dynamics of marine ecosystems, we demonstrate that such models are applicable to a broad suite of major biogeochemical concerns, including anthropogenic change. These results are very promising for the further development and application of reduced‐tracer biogeochemical models that incorporate “sub‐ecosystem‐scale” parameterizations.

58 GEOSCIENCES↗
DOE OSTI2020

Dynamic projection of anthropogenic emissions in China: methodology and 2015–2050 emission pathways under a range of socio-economic, climate policy, and pollution control scenarios

Abstract. Future trends in air pollution and greenhouse gas (GHG)emissions for China are of great concern to the community. A set of globalscenarios regarding future socio-economic and climate developments, combiningshared socio-economic pathways (SSPs) with climate forcing outcomes asdescribed by the Representative Concentration Pathways (RCPs), was createdby the Intergovernmental Panel on Climate Change (IPCC). Chinese researchers have also developed various emission scenarios by considering detailed local environmental and climate policies. However, a comprehensive scenario set connecting SSP–RCP scenarios with local policies and representing dynamic emission changes under local policies is still missing. In this work, to fill this gap, we developed a dynamic projection model, the Dynamic Projection model for Emissions in China (DPEC), to explore China'sfuture anthropogenic emission pathways. The DPEC is designed tointegrate the energy system model, emission inventory model, dynamicprojection model, and parameterized scheme of Chinese policies. The modelcontains two main modules, an energy-model-driven activity rate projectionmodule and a sector-based emission projection module. The activity rateprojection module provides the standardized and unified future energyscenarios after reorganizing and refining the outputs from the energy systemmodel. Here we use a new China-focused version of the Global ChangeAssessment Model (GCAM-China) to project future energy demand and supply inChina under different SSP–RCP scenarios at the provincial level. Theemission projection module links a bottom-up emission inventory model, theMulti-resolution Emission Inventory for China (MEIC), to GCAM-China andaccurately tracks the evolution of future combustion and production technologiesand control measures under different environmental policies. We developedtechnology-based turnover models for several key emitting sectors (e.g.coal-fired power plants, key industries, and on-road transportationsectors), which can simulate the dynamic changes in the unit/vehicle fleetturnover process by tracking the lifespan of each unit/vehicle on an annualbasis. With the integrated modelling framework, we connected five SSP scenarios(SSP1–5), five RCP scenarios (RCP8.5, 7.0, 6.0, 4.5, and 2.6), and threepollution control scenarios (business as usual, BAU; enhanced controlpolicy, ECP; and best health effect, BHE) to produce six combined emissionscenarios. With those scenarios, we presented a wide range of China's futureemissions to 2050 under different development and policy pathways. We foundthat, with a combination of strong low-carbon policy and air pollutioncontrol policy (i.e. SSP1-26-BHE scenario), emissions of major airpollutants (i.e. SO 2 , NO x , PM 2.5 , and non-methane volatile organic compounds – NMVOCs) in China willbe reduced by 34%–66% in 2030 and 58%–87% in 2050 compared to 2015. End-of-pipe control measures are more effective for reducing air pollutant emissions before 2030, while low-carbon policy will play a more important rolein continuous emission reduction until 2050. In contrast, China's emissionswill remain at a high level until 2050 under a reference scenario without activeactions (i.e. SSP3-70-BAU). Compared to similar scenarios set from theCMIP6 (Coupled Model Intercomparison Project Phase 6), our estimates ofemission ranges are much lower than the estimates from the harmonized CMIP6 emissions dataset in2020–2030, but their emission ranges become similar in the year 2050.

54 ENVIRONMENTAL SCIENCES↗