Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “example”

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 1,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Microscale THMC Modeling of Pressure Solution in Salt Rock: Impacts of Geometry and Temperature

Abstract Pressure solution, a mechanism that involves tight coupling between the geometry and thermal-hydro-mechanical-chemical (THMC) processes, plays an important role in diagenesis. In this study, we make the first attempt to conduct microscale THMC modeling to understand and quantify the impacts of geometry and temperature on pressure solution, taking natural salt rock as an example. This modeling capability is achieved by expanding a novel MC code that we developed previously (Hu et al. J Geophys Res: Solid Earth 126:e2021JB023112, 2021) to include temperature effects. We first conduct a simulation of an example that involves a single brine inclusion within a single halite grain and find that the temperature impact is limited for that case. We then extract geometry from an image of a natural salt rock and conduct simulations with different cases: (A) only temperature and no stress, (B) only stress and no temperature, and (C) with both stress and temperature. These different cases result in quite different phenomena. In case A, dissolution and precipitation occur across the entire system due to isolated pore space reaching a localized mass balance between dissolution, precipitation, and diffusion. In case B, intense geometric features (e.g., major asperities, inclusions) in one area undergo stress concentration, thus dominating pressure solution in that area. In case C, pressure solution is spread out at contacting highly stressed geometric features close to the hotter side. We conclude that geometric features dominate stress distribution, thus dominating pressure solution in a natural salt rock that may be affected by the temperature if a sufficient temperature gradient is applied.

58 GEOSCIENCES↗
DOE OSTI2021

A solution framework for linear PDE-constrained mixed-integer problems

Abstract We present a general numerical solution method for control problems with state variables defined by a linear PDE over a finite set of binary or continuous control variables. We show empirically that a naive approach that applies a numerical discretization scheme to the PDEs to derive constraints for a mixed-integer linear program (MILP) leads to systems that are too large to be solved with state-of-the-art solvers for MILPs, especially if we desire an accurate approximation of the state variables. Our framework comprises two techniques to mitigate the rise of computation times with increasing discretization level: First, the linear system is solved for a basis of the control space in a preprocessing step. Second, certain constraints are just imposed on demand via the IBM ILOG CPLEX feature of a lazy constraint callback. These techniques are compared with an approach where the relations obtained by the discretization of the continuous constraints are directly included in the MILP. We demonstrate our approach on two examples: modeling of the spread of wildfire and the mitigation of water contamination. In both examples the computational results demonstrate that the solution time is significantly reduced by our methods. In particular, the dependence of the computation time on the size of the spatial discretization of the PDE is significantly reduced.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A proximal trust-region method for nonsmooth optimization with inexact function and gradient evaluations

Many applications require minimizing the sum of smooth and nonsmooth functions. For example, basis pursuit denoising problems in data science require minimizing a measure of data misfit plus an $\ell^1$-regularizer. Similar problems arise in the optimal control of partial differential equations (PDEs) when sparsity of the control is desired. Here, we develop a novel trust-region method to minimize the sum of a smooth nonconvex function and a nonsmooth convex function. Our method is unique in that it permits and systematically controls the use of inexact objective function and derivative evaluations. When using a quadratic Taylor model for the trust-region subproblem, our algorithm is an inexact, matrix-free proximal Newton-type method that permits indefinite Hessians. We prove global convergence of our method in Hilbert space and demonstrate its efficacy on three examples from data science and PDE-constrained optimization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

McCormick envelopes in mixed-integer PDE-constrained optimization

McCormick envelopes are a standard tool for deriving convex relaxations of optimization problems that involve polynomial terms. Such McCormick relaxations provide lower bounds, for example, in branch-and-bound procedures for mixed-integer nonlinear programs but have not gained much attention in PDE-constrained optimization so far. This lack of attention may be due to the distributed nature of such problems, which on the one hand leads to infinitely many linear constraints (generally state constraints that may be difficult to handle) in addition to the state equation for a pointwise formulation of the McCormick envelopes and renders bound-tightening procedures that successively improve the resulting convex relaxations computationally intractable. We analyze McCormick envelopes for a model problem class that is governed by a semilinear PDE involving a bilinearity and integrality constraints. We approximate the nonlinearity and in turn the McCormick envelopes by averaging the involved terms over the cells of a partition of the computational domain on which the PDE is defined. This yields convex relaxations that underestimate the original problem up to an a priori error estimate that depends on the mesh size of the discretization. These approximate McCormick relaxations can be improved by means of an optimization-based bound-tightening procedure. We show that their minimizers converge to minimizers to a limit problem with a pointwise formulation of the McCormick envelopes when driving the mesh size to zero. We provide a computational example, for which we certify all of our imposed assumptions. The results point to both the potential of the methodology and the gaps in the research that need to be closed. Our methodology provides a framework first for obtaining pointwise underestimators for nonconvexities and second for approximating them with finitely many linear inequalities in an infinite-dimensional setting.

Approximations and Expansions↗
DOE OSTI2024

Evaluating county-level lung cancer incidence from environmental radiation exposure, PM 2.5 , and other exposures with regression and machine learning models

Characterizing the interplay between exposures shaping the human exposome is vital for uncovering the etiology of complex diseases. For example, cancer risk is modified by a range of multifactorial external environmental exposures. Environmental, socioeconomic, and lifestyle factors all shape lung cancer risk. However, epidemiological studies of radon aimed at identifying populations at high risk for lung cancer often fail to consider multiple exposures simultaneously. For example, moderating factors, such as PM 2.5 , may affect the transport of radon progeny to lung tissue. This ecological analysis leveraged a population-level dataset from the National Cancer Institute’s Surveillance, Epidemiology, and End-Results data (2013–17) to simultaneously investigate the effect of multiple sources of low-dose radiation (gross γ activity and indoor radon) and PM 2.5 on lung cancer incidence rates in the USA. County-level factors (environmental, sociodemographic, lifestyle) were controlled for, and Poisson regression and random forest models were used to assess the association between radon exposure and lung and bronchus cancer incidence rates. Tree-based machine learning (ML) method perform better than traditional regression: Poisson regression: 6.29/7.13 (mean absolute percentage error, MAPE), 12.70/12.77 (root mean square error, RMSE); Poisson random forest regression: 1.22/1.16 (MAPE), 8.01/8.15 (RMSE). The effect of PM 2.5 increased with the concentration of environmental radon, thereby confirming findings from previous studies that investigated the possible synergistic effect of radon and PM 2.5 on health outcomes. In summary, the results demonstrated (1) a need to consider multiple environmental exposures when assessing radon exposure’s association with lung cancer risk, thereby highlighting (1) the importance of an exposomics framework and (2) that employing ML models may capture the complex interplay between environmental exposures and health, as in the case of indoor radon exposure and lung cancer incidence.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2024

Opportunities for Earth Observation to Inform Risk Management for Ocean Tipping Points

Abstract As climate change continues, the likelihood of passing critical thresholds or tipping points increases. Hence, there is a need to advance the science for detecting such thresholds. In this paper, we assess the needs and opportunities for Earth Observation (EO, here understood to refer to satellite observations) to inform society in responding to the risks associated with ten potential large-scale ocean tipping elements: Atlantic Meridional Overturning Circulation; Atlantic Subpolar Gyre; Beaufort Gyre; Arctic halocline; Kuroshio Large Meander; deoxygenation; phytoplankton; zooplankton; higher level ecosystems (including fisheries); and marine biodiversity. We review current scientific understanding and identify specific EO and related modelling needs for each of these tipping elements. We draw out some generic points that apply across several of the elements. These common points include the importance of maintaining long-term, consistent time series; the need to combine EO data consistently with in situ data types (including subsurface), for example through data assimilation; and the need to reduce or work with current mismatches in resolution (in both directions) between climate models and EO datasets. Our analysis shows that developing EO, modelling and prediction systems together, with understanding of the strengths and limitations of each, provides many promising paths towards monitoring and early warning systems for tipping, and towards the development of the next generation of climate models.

Wood, Richard A. (ORCID:0000000239609513)↗
DOE OSTI2022

Applications of a dual-column technique in actinide separations

Selective separation of an individual actinide of interest from other actinides and accompanying lanthanides is a challenging task in radiochemistry and radioanalytical chemistry. This paper illustrates a dual column technique (i.e., stacked columns of two appropriate resins) to separate an individual actinide of interest, where the selection of two resins and a common effluent running through the stacked columns are key parameters in design of an appropriate dual column. This paper describes the dual column design and practice with two examples, one for a heavy actinide 249Bk and the other for a light actinide 230U. In the former case, stacked columns of anion exchange (AX) resin and LN resin replaced the traditional CX-AHIB method used at Oak Ridge National Laboratory for 60 years. In the latter case, an elution process with an AX column followed with a stacked column of AX resin and a DGA resin was applied, instead of traditional methods, e.g., PUREX or U-TEVA methods. Application results of the dual column method to the two example actinides are displayed, while the considerations in method design and the required conditions are discussed.

Du, Miting↗
DOE OSTI2020

Uncoupling Electrokinetic Flow Solutions

The continuum-scale electrokinetic porous-media flow and excess charge redistribution equations are uncoupled using eigenvalue decomposition. The uncoupling results in a pair of independent diffusion equations for “intermediate” potentials subject to modified material properties and boundary conditions. The fluid pressure and electrostatic potential are then found by recombining the solutions to the two intermediate uncoupled problems in a matrix-vector multiplication. Expressions for the material properties or source terms in the intermediate uncoupled problem may require extended precision or careful rewriting to avoid numerical cancellation, but the solutions themselves can typically be computed in double precision. The approach works with analytical or gridded numerical solutions and is illustrated through two examples. The solution for flow to a pumping well is manipulated to predict streaming potential and electroosmosis, and a periodic one-dimensional analytical solution is derived and used to predict electroosmosis and streaming potential in a laboratory flow cell subjected to low frequency alternating current and pressure excitation. The examples illustrate the utility of the eigenvalue decoupling approach, repurposing existing analytical solutions or numerical models and leveraging solutions that are simpler to derive for coupled physics.

58 GEOSCIENCES↗
DOE OSTI2026

Consequential improvement acquisition function for efficient multi-fidelity Bayesian optimization

Abstract Surrogate-based Bayesian optimization has been widely applied in design optimization to increase sampling efficiency. However, the cost for each evaluation of the objective function can still be very high when physical experiments or large-scale simulations are involved. Multi-fidelity Bayesian optimization is the new approach to further improve the sampling efficiency by reducing the number of expensive samples at the highest fidelity level and supplementing them with less expensive ones at low-fidelity levels. In this paper, a new consequential improvement (CI) acquisition function is proposed to allow for the simultaneous selection of the solution and the fidelity level in problems with a known hierarchy of fidelity levels. The new CI acquisition function incorporates the consequential effectiveness of objective improvement with the considerations of cost, accuracy, and validity differences between high- and low-fidelity samples in engineering practice. The new method of multi-fidelity Bayesian optimization based on the CI is demonstrated with several analytical and simulation-based design examples. In the simulation-based design optimization example, the results show that the CI acquisition function has a decisive advantage in the sampling efficiency over the other methods of multi-fidelity Bayesian optimization with simultaneous selection. The results indicate that the proposed method is particularly advantageous in solving high-dimensional problems and when large cost ratios between high- and low-fidelity evaluations exist and high-fidelity validation is mandatory. Furthermore, the method robustly avoids the prevalent issue of over sampling at low-fidelity levels.

Aydogdu, Ibrahim [Georgia Institute of Technology,↗
DOE OSTI2020

Heterogeneities at multiple length scales in 2D layered materials: From localized defects and dopants to mesoscopic heterostructures

Two-dimensional (2D) materials hold great promise for applications in optoelectronics, quantum information science, and energy conversion due to their remarkable properties imbued by their physical characteristics. Although heterogeneities in their intrinsic structure are the major challenges limiting their synthesis and predictable properties, they also provide a pathway to controllably tune the properties and broaden the potential of 2D materials. Heterogeneities that can be tailored, including defects, dopants, strain, edges, and layer stackings offer transformative opportunities in heterogeneous 2D materials through the introduction of novel properties for technological applications. This article provides a review of recent progress in studying heterogeneities in 2D materials. The review uses examples from our work to develop a strategy to understand the heterogeneities across multiple length scales to link the effect of heterogeneity at the nanoscale with the macroscale properties of 2D materials. We describe specific types of heterogeneities and explore novel synthesis and processing methods for their controlled production with example of the potential impact and applications enabled by their intriguing properties. Finally, we provide a perspective on how to extend the range of tunable properties through further engineering the heterogeneities in 2D materials.

2D materials↗
DOE OSTI2020

Using yttrium as an indicator to estimate total rare earth element concentration: a case study of anthracite-associated clays from northeastern Pennsylvania

Abstract This study demonstrated using yttrium (Y) as an indicator to estimate the total rare earth element and Y contents (REY) in coal-associated samples and to facilitate selection of samples with high REY assays in a fast and inexpensive manner. More than 10 anthracite-associated samples were collected from each of three Pennsylvanian sites (sites B, J and C) based on Thorium gamma ray logging suggesting high REY content. Several samples from each site were analyzed by ICP-MS to determine the rare earth distribution patterns and to establish the site-specific linear equations of Y and REY. The Y contents of the remaining samples were measured by a portable X-ray fluorescence analyzer, and the REY values were estimated based on the site-specific linear equation developed earlier. R-squared values above 0.70 were obtained for all the estimation equations from all three sites on both a whole sample basis and an ash basis. Previously, ash content has been widely used as an indicator of high REY content. This may not be applicable for a specific site. Site B in this study is an example where ash contents could not be statistically correlated with REY, so using Y for estimation is more applicable. The demonstrated sample screening process is suitable for samples from sites that share more similar distribution patterns (either MREY or LREY or HREY) as well as for samples from sites that share multiple distribution patterns (LREY/MREY/HREY) depending on the desirable accuracy. The demonstrated process lowers the analytical cost from $70 to 80 dollars per sample to $10–15 per sample while significantly reducing the processing time and acid consumption for ICP digestion. This is particularly true when a relatively large sample size is involved, for example, 100 samples from one site analyzed by ICP-MS/OES.

Yang, Xiaojing↗
DOE OSTI2022

Protein–chromophore interactions controlling photoisomerization in red/green cyanobacteriochromes

Abstract Photoreceptors in the phytochrome superfamily use 15,16-photoisomerization of a linear tetrapyrrole (bilin) chromophore to photoconvert between two states with distinct spectral and biochemical properties. Canonical phytochromes include master regulators of plant growth and development in which light signals trigger interconversion between a red-absorbing 15Z dark-adapted state and a metastable, far-red-absorbing 15E photoproduct state. Distantly related cyanobacteriochromes (CBCRs) carry out a diverse range of photoregulatory functions in cyanobacteria and exhibit considerable spectral diversity. One widespread CBCR subfamily typically exhibits a red-absorbing 15Z dark-adapted state similar to that of phytochrome that gives rise to a distinct green-absorbing 15E photoproduct. This red/green CBCR subfamily also includes red-inactive examples that fail to undergo photoconversion, providing an opportunity to study protein–chromophore interactions that either promote photoisomerization or block it. In this work, we identified a conserved lineage of red-inactive CBCRs. This enabled us to identify three substitutions sufficient to block photoisomerization in photoactive red/green CBCRs. The resulting red-inactive variants faithfully replicated the fluorescence and circular dichroism properties of naturally occurring examples. Converse substitutions restored photoconversion in naturally red-inactive CBCRs. This work thus identifies protein–chromophore interactions that control the fate of the excited-state population in red/green cyanobacteriochromes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Accident-Tolerant Fuel

Advanced light water reactor (LWR) fuels have been developed incrementally for more than 60 years before the term Accident Tolerant Fuel (ATF) was created. There are many different categories of ATF concepts, which range from near-term marginal enhancements to fission gas release in fuel or hydrogen pickup in the cladding, to transformational fuel types with new engineered fission product barriers. In this work we characterize these different types of ATF, and give several important examples for each.We conclude that transformational ATF concepts have the greatest potential to enhance the traditional definition of defense-in-depth by enhancing current barriers and/or by providing additional barriers to fission product release. This objective can primarily be achieved by introducing a new fission product barrier that enhances defense-in-depth by design, for example the added barriers to fission product release in some engineered ATF concepts. We recommend future gap analyses and a phenomena identification and ranking to assess the applicability of existing safe acceptable fuel design limits and the standard review plan to specific ATF concepts.

Brown, Nicholas↗
DOE OSTI2023

Relativistic Effects From Coupled-Cluster Theory

We discuss the theory and computational challenges of the relativistic coupled-cluster methods. Example calculations of heavy-atom-containing molecules are then presented to demonstrate the importance of scalar-relativistic, spin-orbit coupling, and electron-correlation effects on molecular properties as well as the applicability and usefulness of relativistic coupled-cluster methods in calculations aiming at high-accuracy results. Here, a unique applicability of the spinor-based relativistic coupled-cluster methods is also highlighted using the calculations of open shell actinide-containing small molecules. As a result, a summary is given together with an outlook into future developments.

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

Differentiating densities on smooth manifolds

Lebesgue integration of derivatives of strongly-oscillatory functions is a recurring challenge in computational science and engineering. Integration by parts is an effective remedy for huge computational costs associated with Monte Carlo integration schemes. In case of Lebesgue integrals over a smooth manifold, however, integration by parts gives rise to a derivative of the density implied by charts describing the domain manifold. This paper focuses on the computation of that derivative, which we call the density gradient function, on general smooth manifolds. We analytically derive formulas for the density gradient and present examples of manifolds determined by popular differential equation-driven systems. Furthermore, we highlight the significance of the density gradient by demonstrating a numerical example of Monte Carlo integration involving oscillatory integrands.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Cold spray technology in nuclear energy applications: A review of recent advances

Cold spray technology is a promising solid-state, powder-based deposition methods for fabrication of coatings, near-net-shape manufacturing, and component repair. In recent years, this technology has been investigated for many applications in the nuclear energy sector. For example, it has been explored for the deposition of coatings of corrosion and oxidation materials on zirconium-alloys fuel claddings in light water reactors (LWR) for achieving more accident-tolerant fuel (ATF) cladding. Other examples of its use for nuclear energy applications include coatings for enhancements in wear resistance and heat transfer, near-net shape manufacturing of next-generation fuel claddings of oxide dispersion-strengthened (ODS) steels, and for mitigation and repair of potential chloride induced stress corrosion cracking in used fuel storage canister systems. This work reviews published literature on the application of cold spray process for nuclear energy applications, while pointing to gaps that have yet to be addressed.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Parametrizations of the spin density matrix

Highlights: • Proposed two methods for parametrizing the SDM of spinning, unstable particles. • All non-negativity conditions fulfilled; adaptability to various constraints. • First method: based on a theorem; applicable in a simple and flexible way. • Second method: particularly suitable for especial conditions. • Applications: two of them interesting in the search for new physics. We propose for the spin density matrix two parametrizations which automatically fulfill the non-negativity conditions, without setting any bound on the parameters. The first one relies on a theorem, that we prove, and it is rather simple and easily adaptable to some specific reactions, where, for example, parity is conserved or angular momentum conservation entails selection rules. Moreover, in the case when the rank is less than the order of the density matrix, we show how to improve the fits to the data, either by implementing previous suggestions, or by elaborating an alternative method, for which we prove a second theorem. Our second parametrization is a variant of previous treatments, it appears suitable for some particular processes. Moreover, we discuss about the possibility of inferring the elements of the density matrix from the differential decay width. Last, we illustrate various examples of current interest, both in strong and weak decays; some of them may be helpful in the investigation of physics beyond the standard model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Modeling and analysis of heat emissions from buildings to ambient air

Heat emissions from buildings is a significant source of anthropogenic heat influencing the urban microclimate; however, they are usually oversimplified in urban climate and microclimate modeling. This study developed a bottom-up physics-based approach to calculate heat emissions from buildings to the ambient air and implemented the approach in EnergyPlus. In this work, a simple result verification was conducted by comparing the EnergyPlus simulated results against the spreadsheet calculations. Simulations covering 16 commercial building types, four climates, and two energy efficiency levels were conducted to understand and evaluate the building heat emissions and their temporal patterns as well as three major components: (1) building envelope (convective heat transfer to ambient air), (2) zones (air exfiltration and exhaust air), and (3) HVAC systems (relief air and heat rejection from condensers or cooling towers). The main findings are: (1) heat emissions are usually higher than the site energy use (about 2.5 times), and their dynamics should be considered; (2) building characteristics and their energy systems lead to differences in heat emission contributions from the three components, and their dynamics, for example, in the warehouse models, the envelope component accounts for 90.4%, while it is 12.7% for the large office models; (3) for most building typologies, the climate has a strong impact on heat emissions, for example, buildings with dominant heat emissions from the zone exhaust air and/or the HVAC reject heat, a general decrease in heat emissions in hotter climates is observed, while envelope-dominated buildings show the opposite; and (4) building technologies that reduce energy use in buildings may perform differently in reducing heat emissions. The developed heat emissions calculation method can be adopted in EnergyPlus and most other building energy modeling programs. It can provide dynamic building heat emissions as an input to urban climate computational fluid dynamics (CFD) models at a higher spatial and temporal resolution than is currently available, to improve the simulation accuracy of the urban microclimate and capture the urban heat island effect and urban overheating.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗