Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Spreading models”

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

Factors Governing Cloud Growth and Entrainment Rates in Shallow Cumulus and Cumulus Congestus During GoAmazon2014/5

Shallow cumulus and cumulus congestus clouds play an important role in the large-scale tropical circulation by mixing heat and moisture vertically and preconditioning the environment for deeper convection. Different representations of these shallow clouds account for much of the spread in General Circulation Model (GCMs) climate sensitivity, potentially because of how entrainment is represented in GCM parameterizations. This study uses observations from the Department of Energy's Atmospheric Radiation Measurement (ARM) mobile facility deployed at Manacapuru, Brazil, during the Green Ocean Amazon (GoAmazon2014/5) Campaign. Environmental thermodynamic profiles and observations of cloud top height (CTH) are used to constrain an entraining plume model to estimate bulk entrainment rates (ERs). Estimates of CTH are obtained from a combination of vertically pointing W-band ARM cloud radar and 1,290 MHz Radar Wind Profiler observations. A combination of radiosonde, microwave radiometer profiler, and microwave radiometer observations provides new best estimates of the environmental thermodynamic state. We quantify uncertainty in ERs considering uncertainties in estimated CTH, environmental thermodynamic properties, and assumed initial parcel characteristics. We find ERs ranging from 0.16 to 2.8 km -1 with an average of 0.58 ± 0.10 km -1 over a selected population of 469 shallow cumulus and cumulus congestus clouds. Using the retrieved estimates of ER, we evaluate several entrainment closures that are currently used in atmospheric models or have been proposed based on theory or large eddy simulation. Finally, entrainment rates in cumulus clouds are weakly correlated with low-level buoyancy, cloud depth, and cloud size.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Accurate shear estimation with fourth-order moments

ABSTRACT As imaging surveys progress in exploring the large-scale structure of the Universe through the use of weak gravitational lensing, achieving sub-per cent accuracy in estimating shape distortions caused by lensing, or shear, is imperative for precision cosmology. In this paper, we extend the Fourier power function shapelets (FPFS) shear estimator using fourth-order shapelet moments and combine it with the original second-order shear estimator to reduce galaxy shape noise. We calibrate this novel shear estimator analytically to a sub-per cent level-accuracy using the AnaCal framework. This higher order shear estimator is tested with realistic image simulations, and after analytical correction for the detection/selection bias and noise bias, the multiplicative shear bias $|m|$ is below $3\times 10^{-3}$ (99.7 per cent confidence interval) for both isolated and blended galaxies. Once combined with the second-order FPFS shear estimator, the shape noise is reduced by $\sim 35~{{\ \rm per\ cent}}$ for isolated galaxies in simulations with Hyper Suprime-Cam and Vera C. Rubin Observatory Legacy Survey of Space and Time observational conditions. However, for blended galaxies, the effective number density does not significantly improve with the combination of the two estimators. Based on these results, we recommend exploration of how this framework can further reduce the systematic uncertainties in shear due to point spread function leakage and modelling error, and potentially provide improved precision in shear inference in high-resolution space-based images.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Fully Bayesian Analysis With Model Inadequacy Correction For Nuclear Graphite Property Models With Hierarchical Variance Structure

Nuclear-grade graphites are extensively utilized in the core designs of various advanced nuclear reactors. Within the reactor environment, graphite is subjected to prolonged exposure to extreme conditions, including high temperatures, radiation, and potentially molten salt and oxygen. Such exposure can induce several degradation mechanisms in graphite, such as nonuniform volumetric strains caused by irradiation and thermal expansion, leading to stresses that may compromise the performance of graphite components. Assessing component integrity, forecasting component performance over the reactor's lifespan, and developing design standards necessitate robust tools for predicting fracture initiation and propagation in graphite structural components within nuclear reactors. This code enables the Bayesian calibration of properties for nuclear-grade graphites. Using a hierarchical Bayesian approach, multiple experimental data sources are combined to develop Gaussian process models for the properties. Using the Kennedy O'Hagan framework, the uncertainties due inadequacies in the model and the inherent spread in the experimental data are quantified.

Dhulipala, Som Lakshmi NarasimhaLakshmi Narasimha ↗
DOE OSTI2024

Towards reverse mode automatic differentiation of Kokkos-based codes

Derivative computation is a key component of optimization, sensitivity analysis, uncertainty quantification, and the solving of nonlinear problems. Automatic differentiation (AD) is a powerful technique for evaluating such derivatives, and in recent years, has been integrated into programming environments such as Jax, PyTorch, and TensorFlow to support derivative computations needed for training of machine learning models, facilitating wide-spread use of these technologies. The C++ language has become the de facto standard for scientific computing due to numerous factors, yet language complexity has made the wide-spread adoption of AD technologies for C++ difficult, hampering the incorporation of powerful differentiable programming approaches into C++ scientific simulations. This is exacerbated by the increasing emergence of architectures, such as GPUs, with limited memory capabilities and requiring massive thread-level concurrency. C++ AD tools must effectively use these environments to bring novel scientific simulations to next-generation DOE experimental and observational facilities. In this project, we investigated source transformation-based automatic differentiation using LLVM compiler infrastructure to automatically generate portable and efficient gradient computations of Kokkos-based code. We have demonstrated that our proposed strategy is feasible by investigating the usage of a prototype LLVM-based source transformation tool to generate gradients of simple functions made of sequences of simple Kokkos parallel regions. Speedups of up to 500x compared to Sacado were observed on NVIDIA V100 GPU.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

No constraint on long-term tropical land carbon-climate feedback uncertainties from interannual variability

Unraveling drivers of the interannual variability of tropical land carbon cycle is critical for understanding land carbon-climate feedbacks. Here we utilize two generations of factorial model experiments to show that interannual variability of tropical land carbon uptake under both present and future climate is consistently dominated by terrestrial water availability variations in Earth system models. The magnitude of this interannual sensitivity of tropical land carbon uptake to water availability variations under future climate shows a large spread across the latest 16 models (2.3 ± 1.5 PgC/yr/Tt H 2 O), which is constrained to 1.3 ± 0.8 PgC/yr/Tt H 2 O using observations and the emergent constraint methodology. However, the long-term tropical land carbon-climate feedback uncertainties in the latest models can no longer be directly constrained by interannual variability compared with previous models, given that additional important processes are not well reflected in interannual variability but could determine long-term land carbon storage. Our results highlight the limited implication of interannual variability for long-term tropical land carbon-climate feedbacks and help isolate remaining uncertainties with respect to water limitations on tropical land carbon sink in Earth system models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Analysis of picosecond coherent anti-Stokes Raman spectra for gas-phase diagnostics

We present a hybrid frequency- and time-domain solution, applicable to the case of picosecond coherent anti-Stokes Raman scattering (CARS), for gas-phase diagnostics. A solution has been derived based on both physical arguments and four-wave mixing equations for picosecond CARS, with pulse durations that are comparable to the dephasing time scale for gas-phase Raman coherence—a regime where commonly employed solutions for impulsive (femtosecond) or cw (nanosecond) pump/Stokes forcing are not strictly valid. We present the ps-CARS spectrum in the form of incoherent sums of CARS intensity spectra, calculated from the fundamental solution for impulsive pump/Stokes Raman preparation. The solution was examined for temperatures from 1000–3000 K, for four plausible experimental configurations, with laser pulse durations of 50–150 ps, and probe pulse delays from −20 to 240 ps. Approximations based on cw and impulsive pump/Stokes preparation to fit picosecond CARS spectra at atmospheric pressure were examined and the relative thermometric accuracy and computational cost of these approximations were quantified for the case of a zero nonresonant CARS contribution, and a nonresonant susceptibility equal to 10% of the Raman-resonant value at the N 2 bandhead. The nanosecond CARS approximation can result in large fitting errors when the probe pulse time delay is less than the probe pulse duration. Errors as large as 10–20% are observed in the fit temperatures for a zero picosecond probe pulse delay, when the nonresonant background is neglected, largely due to an inability of the time-independent cw model to capture transient frequency spread dephasing effects at the Q -branch bandhead. The inclusion of a nonresonant background results in 40–60% thermometry errors with a nanosecond model at a zero-probe delay. Time-dependent impulsive calculations used for femtosecond CARS better approximate the structure of the N 2 bandhead, reducing temperature fitting errors to 5–10% at a short probe pulse delay. The impulsive approximation results in errors up to 10% at intermediate probe pulse delays, where the coherence of the pump and probe pulses leads to multiple terms in the picosecond CARS solution. Both approximations improve as the probe pulse delay exceeds the probe duration. The nanosecond approximation results in a 2–3% error, while the impulsive model results in differences of less than 1% in some cases. Fits to experimental data obtained using short, ∼60ps pulses at a zero probe time delay and longer 100 ps pulses at a substantial 200 ps delay are presented with accuracies of 1–3% in the fit temperature.

Kearney, Sean P.↗
DOE OSTI2020

Probabilistic Topography Project Update 5/20 [Slides]

Problem/Goal: Earth System Models (ESMs) do not accurately predict fire spread. There are no effects from topography/vegetation heterogeneity below the model’s scale (100-200km). Use physics based models (FIRETEC and QUICFIRE) to produce an archive of correction factors for SPITFIRE (chosen ESM). This project will: Use datasets to create PDFs of slope & aspect, vegetation distribution; Create joint distribution and sample to make new domain; Run FIRETEC & QUICFIRE on created domain, record output (area burned, intensity, etc.); and Use output to create PDF of correction factors for SPITFIRE runs.

58 GEOSCIENCES↗
DOE OSTI2020

Potentially adaptive SARS-CoV-2 mutations discovered with novel spatiotemporal and explainable AI models

Abstract Background A mechanistic understanding of the spread of SARS-CoV-2 and diligent tracking of ongoing mutagenesis are of key importance to plan robust strategies for confining its transmission. Large numbers of available sequences and their dates of transmission provide an unprecedented opportunity to analyze evolutionary adaptation in novel ways. Addition of high-resolution structural information can reveal the functional basis of these processes at the molecular level. Integrated systems biology-directed analyses of these data layers afford valuable insights to build a global understanding of the COVID-19 pandemic. Results Here we identify globally distributed haplotypes from 15,789 SARS-CoV-2 genomes and model their success based on their duration, dispersal, and frequency in the host population. Our models identify mutations that are likely compensatory adaptive changes that allowed for rapid expansion of the virus. Functional predictions from structural analyses indicate that, contrary to previous reports, the Asp 614 Gly mutation in the spike glycoprotein (S) likely reduced transmission and the subsequent Pro 323 Leu mutation in the RNA-dependent RNA polymerase led to the precipitous spread of the virus. Our model also suggests that two mutations in the nsp13 helicase allowed for the adaptation of the virus to the Pacific Northwest of the USA. Finally, our explainable artificial intelligence algorithm identified a mutational hotspot in the sequence of S that also displays a signature of positive selection and may have implications for tissue or cell-specific expression of the virus. Conclusions These results provide valuable insights for the development of drugs and surveillance strategies to combat the current and future pandemics.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Reliable Measures of Spread in High Dimensional Latent Spaces

Understanding geometric properties of the latent spaces of natural language processing models allows the manipulation of these properties for improved performance on downstream tasks. One such property is the amount of data spread in a model’s latent space, or how fully the available latent space is being used. We demonstrate that the commonly used measures of data spread, average cosine similarity and a partition function min/max ratio I (V), do not provide reliable metrics to compare the use of latent space across data distributions. We propose and examine six alternative measures of data spread, all of which improve over these current metrics when applied to seven synthetic data distributions. Of our proposed measures, we recommend one principal component-based measure and one entropy-based measure that provide reliable, relative measures of spread and can be used to compare models of different sizes and dimensionalities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Evaluating Climate Models’ Cloud Feedbacks Against Expert Judgment

Abstract The persistent and growing spread in effective climate sensitivity (ECS) across global climate models necessitates rigorous evaluation of their cloud feedbacks. Here we evaluate several cloud feedback components simulated in 19 climate models against benchmark values determined via an expert synthesis of observational, theoretical, and high‐resolution modeling studies. We find that models with smallest feedback errors relative to these benchmark values generally have moderate total cloud feedbacks (0.4–0.6 W m −2 K −1 ) and ECS (3–4 K). Those with largest errors generally have total cloud feedback and ECS values that are too large or too small. Models tend to achieve large positive total cloud feedbacks by having several cloud feedback components that are systematically biased high rather than by having a single anomalously large component, and vice versa. In general, better simulation of mean‐state cloud properties leads to stronger but not necessarily better cloud feedbacks. The Python code base provided herein could be applied to developmental versions of models to assess cloud feedbacks and cloud errors and place them in the context of other models and of expert judgment in real‐time during model development.

58 GEOSCIENCES↗
DOE OSTI2024

Congruity of genomic and epidemiological data in modelling of local cholera outbreaks

Cholera continues to be a global health threat. Understanding how cholera spreads between locations is fundamental to the rational, evidence-based design of intervention and control efforts. Traditionally, cholera transmission models have used cholera case-count data. More recently, whole-genome sequence data have qualitatively described cholera transmission. Integrating these data streams may provide much more accurate models of cholera spread; however, no systematic analyses have been performed so far to compare traditional case-count models to the phylodynamic models from genomic data for cholera transmission. Here, we use high-fidelity case-count and whole-genome sequencing data from the 1991 to 1998 cholera epidemic in Argentina to directly compare the epidemiological model parameters estimated from these two data sources. We find that phylodynamic methods applied to cholera genomics data provide comparable estimates that are in line with established methods. Our methodology represents a critical step in building a framework for integrating case-count and genomic data sources for cholera epidemiology and other bacterial pathogens.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

IoT-based retrofit information diffusion in future smart communities

Community-scale building retrofits are not merely scaled-up versions of single-building retrofits. They involve complex challenges, such as reconciling individual interests with collective goals and managing the dynamic interplay between buildings through mechanisms like power grids and social connections. Internet of Things (IoT) connectivity holds the potential to leverage these interplays to balance individual and collective interests effectively in smart communities. One critical aspect of this interplay is information diffusion, which shapes how retrofit decisions spread among neighbors, influencing individual choices and ultimately impacting community-level retrofit outcomes. In other words, IoT-based smart devices automatically push tailored retrofit notifications to homeowners, which completely changes the format of information diffusion in the future. To investigate this influence by such information diffusion, the study used CityBES to simulate energy performance for different retrofits and applied an information diffusion model to analyze how decisions spread in a networked community of 192 buildings. The diffusion process was modeled on a weighted, directed network, capturing the dynamics of information flow and decision-making across 16 scenarios. Individual retrofit benefits were evaluated through payback years, while community-level retrofit outcomes were assessed using greenhouse gas (GHG) emission reductions. The results demonstrate that easier information diffusion among neighbors encourages households to prioritize retrofit measures that align with the majority’s optimal choices, even at the expense of individual financial benefits. In this case, such collective prioritization enhanced community-level retrofit performance, increasing GHG emission reductions by up to 29.4 %. However, this improvement came with trade-offs, as the average payback period for households extended by approximately 1.74 years. These findings highlight the potential of IoT-based information diffusion in future smart communities to coordinate individual interests with collective goals, ultimately accelerating community-level building retrofits.

Shu, Lei↗
DOE OSTI2020

Effect of Particle Spreading Dynamics on Powder Bed Quality in Metal Additive Manufacturing

Powder spreading precedes creation of every new layer in powder bed additive manufacturing (AM). The powder spreading process can lead to powder layer defects such as porosity, poor surface roughness and particle segregation. Therefore, the creation of homogeneous layers is the first task for optimal part printing. Discrete element methods (DEM) powder spreading simulations are typically limited to a single layer and/or small number of particles. Therefore, results from such model configurations may not be generalized to multiple layer processes. In this study, a computationally efficient multi-layer powder spreading DEM simulation model is proposed. The model is calibrated experimentally using static Angle of Repose measurements. The adhesion model parameter, cohesive energy density is related to adhesive surface energy and strain energy release rate parameters. The model results show that interaction between particle and the powder spreading rake leads to noticeable variation in packing density, surface roughness, dynamic angle of repose (AOR), particle size distribution, and particle segregation. Finally, the powder model is experimentally validated using a recoater spreading rig to measure the dynamic AOR at spreading speeds consistent with recoating speeds and layer heights used in AM processes.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Effective radiative forcing and adjustments in CMIP6 models

Abstract. The effective radiative forcing, which includes the instantaneous forcing plus adjustments from the atmosphere and surface, has emerged as the key metric of evaluating human and natural influence on the climate. We evaluate effective radiative forcing and adjustments in 17 contemporary climate models that are participating in the Coupled Model Intercomparison Project (CMIP6) and have contributed to the Radiative Forcing Model Intercomparison Project (RFMIP). Present-day (2014) global-mean anthropogenic forcing relative to pre-industrial (1850) levels from climate models stands at 2.00 (±0.23) W m−2, comprised of 1.81 (±0.09) W m−2 from CO2, 1.08 (± 0.21) W m−2 from other well-mixed greenhouse gases, −1.01 (± 0.23) W m−2 from aerosols and −0.09 (±0.13) W m−2 from land use change. Quoted uncertainties are 1 standard deviation across model best estimates, and 90 % confidence in the reported forcings, due to internal variability, is typically within 0.1 W m−2. The majority of the remaining 0.21 W m−2 is likely to be from ozone. In most cases, the largest contributors to the spread in effective radiative forcing (ERF) is from the instantaneous radiative forcing (IRF) and from cloud responses, particularly aerosol–cloud interactions to aerosol forcing. As determined in previous studies, cancellation of tropospheric and surface adjustments means that the stratospherically adjusted radiative forcing is approximately equal to ERF for greenhouse gas forcing but not for aerosols, and consequentially, not for the anthropogenic total. The spread of aerosol forcing ranges from −0.63 to −1.37 W m−2, exhibiting a less negative mean and narrower range compared to 10 CMIP5 models. The spread in 4×CO2 forcing has also narrowed in CMIP6 compared to 13 CMIP5 models. Aerosol forcing is uncorrelated with climate sensitivity. Therefore, there is no evidence to suggest that the increasing spread in climate sensitivity in CMIP6 models, particularly related to high-sensitivity models, is a consequence of a stronger negative present-day aerosol forcing and little evidence that modelling groups are systematically tuning climate sensitivity or aerosol forcing to recreate observed historical warming.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Mid-latitude clouds contribute to Arctic amplification via interactions with other climate feedbacks

Abstract Traditional feedback analyses, which assume that individual climate feedback mechanisms act independently and add linearly, suggest that clouds do not contribute to Arctic amplification. However, feedback locking experiments, in which the cloud feedback is disabled, suggest that clouds, particularly outside of the Arctic, do contribute to Arctic amplification. Here, we reconcile these two perspectives by introducing a framework that quantifies the interactions between radiative feedbacks, radiative forcing, ocean heat uptake, and atmospheric heat transport. We show that including the cloud feedback in a comprehensive climate model can result in Arctic amplification because of interactions with other radiative feedbacks. The surface temperature change associated with including the cloud feedback is amplified in the Arctic by the surface-albedo, Planck, and lapse-rate feedbacks. A moist energy balance model with a locked cloud feedback exhibits similar behavior as the comprehensive climate model with a disabled cloud feedback and further indicates that the mid-latitude cloud feedback contributes to Arctic amplification via feedback interactions. Feedback locking in the moist energy balance model also suggests that the mid-latitude cloud feedback contributes substantially to the intermodel spread in Arctic amplification across comprehensive climate models. These results imply that constraining the mid-latitude cloud feedback will greatly reduce the intermodel spread in Arctic amplification. Furthermore, these results highlight a previously unrecognized non-local pathway for Arctic amplification.

58 GEOSCIENCES↗
DOE OSTI2020

The Influence of Aqueous Dispersivity on Transport of Radiological Material at LANL MDA G - 20450

Aqueous dispersivity is a phenomenon by which a solute moving in groundwater mixes with water that does not contain the solute and becomes diluted. Solutes take paths of varying lengths through pore networks or experience slower velocities due to interactions with pore walls. Aqueous dispersivity leads to breakthrough curves with earlier arrival than expected in purely advective systems and persistent 'tails,' or late-arriving, dilute concentrations of the solute. Dispersivity can be longitudinal (in the direction of flow) or transverse (orthogonal to flow). LANL has disposed of radioactive material at Material Disposal Area (MDA) G since the late 1950's. Most of the disposed materials at MDA G are low-level waste (LLW), but other types of radioactive wastes are disposed or stored at the site. For example, there is a collection of remote-handled transuranic waste (RH-TRU) waste in 33 shafts, which were intended to be stored retrievably. Aqueous-phase radionuclides included in models of contaminant transport at MDA G are subject to dispersive processes in porous media of the vadose zone and saturated zone. To model dispersivity at MDA G probabilistically, we created site-specific distributions of aqueous dispersivity. First, we performed a comprehensive review of peer-reviewed publications, reports from other performance assessments of MDA G, and documentation from sites that were considered similar based on expert opinion, such as Yucca Mountain. Data were included from a wide range of sources, so that distributions developed based on the data reflect the current state of knowledge in properties of the subsurface environment. Data were grouped based on value type (e.g., determined via experiment or used in a different model), rock type, zone type (saturated or unsaturated), and site, and were weighted in the distributions according to their relevance to the site and zone type. We combined Monte Carlo simulation and bootstrapping techniques to produce approximate distributions of mean aqueous dispersivity. Distributions of dispersivity for saturated and unsaturated zone transport are differentiated by data weight and the length scales over which transport occurs. We find that underestimating or disregarding transverse dispersivity may lead to overstating simulated well capture in the saturated zone at MDA G. Neglecting transverse dispersivity is a common assumption due in part to use of one-dimensional transport models that do not consider solute spreading orthogonal to flow. Also, we find that variability in longitudinal dispersivity is associated with substantial travel time alterations in the MDA G unsaturated zone and in simulated contaminant breakthrough times at a downstream observation plane in the MDA G saturated zone. Finally, at large scales (>1000 m), a linear trend between the dispersivity coefficient and length over which transport occurs may overestimate the dispersivity coefficient. Future modeling work with application to solute spreading could include interacting advective and diffusive zones in the highly fractured basalt below MDA G. (authors)

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

Simulations of flaming combustion and flaming-to-smoldering transition in wildland fire spread at flame scale

Our objective in the present study is to provide basic insights into the coupling between external-gas and solid biomass vegetation processes that control the dynamics of flame spread in wildland fire problems. Here we focus on a modeling approach that resolves processes occurring at vegetation and flame scales, i.e., the formation of flammable vapors due to the thermal degradation of the solid biomass, the subsequent combustion in ambient air, the thermal feedback to the biomass through radiative and convective heat transfer, and the possible transition from flaming combustion (taking place outside of the solid biomass) to smoldering combustion (taking place inside the solid biomass). The capability uses a multiphase combustion framework and treats external-gas processes through a Large Eddy Simulation solver and solid biomass processes through a discrete particle model. The discrete particle model adopts a one-dimensional porous medium formulation, includes descriptions of drying, thermal pyrolysis, oxidative pyrolysis, and char oxidation, as well as a description of the external-gas-to-solid-biomass diffusion of oxygen mass; the discrete particle model thereby provides a treatment of in-depth oxidative processes and allows the simulation of smoldering combustion. The modeling capability is applied to the simulation of fire spread across a surrogate biomass vegetation bed corresponding to a discrete array of cylindrical-shaped, vertically-oriented, pine wood sticks, characterized by a monomodal size distribution, in horizontal flat terrain and under wind-aided conditions. The numerical results demonstrate that the model can simulate successful flaming-to-smoldering transition followed by complete biomass consumption.

01 COAL, LIGNITE, AND PEAT↗