Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Refining”

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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Minimum feature size control in level set topology optimization via density fields

A level set topology optimization approach that uses an auxiliary density field to nucleate holes during the optimization process and achieves minimum feature size control in optimized designs is explored. The level set field determines the solid-void interface and the density field describes the distribution of a fictitious porous material using the solid isotropic material with penalization. These fields are governed by two sets of independent optimization variables which are initially coupled using a penalty for hole nucleation. The strength of the density field penalization and projection is gradually increased during the optimization process to promote a 0-1 density distribution. In addition, a second penalty regulates the evolution of the density field in the void phase. The treatment of the density field combined with the second penalty mitigate the appearance of small design features. The minimum feature size of optimized designs is controlled by the radius of the linear filter applied to the density optimization variables. The structural response is predicted by the extended finite element method, the sensitivities by the adjoint method, and the optimization variables are updated by a gradient-based optimization algorithm. Numerical examples investigate the robustness of this approach with respect to algorithmic parameters and mesh refinement. The results show the applicability of the combined density level set topology optimization approach for both optimal hole nucleation and for minimum feature size control in 2D and 3D. This comes, however, at the cost of a more complex problem formulation and additional computational cost due to an increased number of optimization variables.

42 ENGINEERING↗
DOE OSTI2020

Aggregation and analysis of indication-symptom relationships for drugs approved in the USA

Abstract Purpose Drug indications and disease symptoms often confound adverse event reports in real-world datasets, including electronic health records and reports in the FDA Adverse Event Reporting System (FAERS). A thorough, standardized set of indications and symptoms is needed to identify these confounders in such datasets for drug research and safety assessment. The aim of this study is to create a comprehensive list of drug-indication associations and disease-symptom associations using multiple resources, including existing databases and natural language processing. Methods Drug indications for drugs approved in the USA were extracted from two databases, RxNorm and Side Effect Resource (SIDER). Symptoms for these indications were extracted from MedlinePlus and using natural language processing from PubMed abstracts. Results A total of 1361 unique drugs, 1656 unique indications, and 2201 unique symptoms were extracted from a wide variety of MedDRA System Organ Classes. Text-mining precision was maximized at 0.65 by examining Term Frequency Inverse Document Frequency (TF-IDF) scores of the disease-symptom associations. Conclusion The drug-indication associations and disease-symptom associations collected in this study may be useful in identifying confounders in other datasets, such as safety reports. With further refinement and additional drugs, indications, and symptoms, this dataset may become a quality resource for disease symptoms.

Punyala, Ananth↗
DOE OSTI2025

The Observed and Projected Changes of Global Monsoons: Current Status and Future Perspectives

The global monsoon system, encompassing the Asian-Australian, African, and American monsoons, sustains two-thirds of the world’s population by regulating water resources and agriculture. Monsoon anomalies pose severe risks, including floods and droughts. Recent research associated with the implementation of the Global Monsoons Model Intercomparison Project under the umbrella of CMIP6 has advanced our understanding of its historical variability and driving mechanisms. Observational data reveal a 20th-century shift: increased rainfall pre-1950s, followed by aridification and partial recovery post-1980s, driven by both internal variability (e.g., Atlantic Multidecadal Oscillation) and external forcings (greenhouse gases, aerosols), while ENSO drives interannual variability through ocean-atmosphere interactions. Future projections under greenhouse forcing suggest long-term monsoon intensification, though regional disparities and model uncertainties persist. Models indicate robust trends but struggle to quantify extremes, where thermodynamic effects (warming-induced moisture rise) uniformly boost heavy rainfall, while dynamical shifts (circulation changes) create spatial heterogeneity. Volcanic eruptions and proposed solar radiation modification (SRM) further complicate predictions: tropical eruptions suppress monsoons, whereas high-latitude events alter cross-equatorial flows, highlighting unresolved feedbacks. The emergent constraint approach is booming in terms of correcting future projections and reducing uncertainty with respect to the global monsoons. Critical challenges remain. Model biases and sparse 20th-century observational data hinder accurate attribution. The interplay between natural variability and anthropogenic forcings, along with nonlinear extreme precipitation risks under warming, demands deeper mechanistic insights. Additionally, SRM’s regional impacts and hemispheric monsoon interactions require systematic evaluation. Addressing these gaps necessitates enhanced observational networks, refined climate models, and interdisciplinary efforts to disentangle multiscale drivers, ultimately improving resilience strategies for monsoon-dependent regions.

climate extreme events↗
DOE OSTI2021

High-pressure crystal structure and equation of state of ferromagnesian jeffbenite: implications for stability in the transition zone and uppermost lower mantle

Jeffbenite, ideally Mg 3 Al 2 Si 3 O 12 , has been identified as inclusions in super-deep diamonds originating from depths that exceed 300 km. Although Mg-end member jeffbenite has limited stability at upper-mantle conditions, iron-bearing jeffbenite may have broader P–T stability that extends to the transition zone or uppermost lower mantle, incorporating significant amounts of ferric iron. Using synchrotron-based, single-crystal X-ray diffraction (XRD) and synchrotron Mössbauer spectroscopy (SMS) at pressures up to 29 GPa, we report the crystal structure, compressibility, and likely spin transition of iron in ferromagnesian jeffbenite (Mg 2.32 Al 0.03 Fe 2+ 1.28 Fe 3+ 1.77 Si 2.85 O 12 ). High-pressure structure refinements reveal that Fe 3+ substitution for Si in the T2 site, which shares edges with the M2 octahedron, likely stabilizes jeffbenite at high pressure, because it increases the cation-to-cation distance between these sites. Although ferromagnesian jeffbenite does not undergo a structural phase transition below 30 GPa, SMS hyperfine parameters suggest the onset of an electronic spin transition of iron from high-spin (HS) to low-spin (LS) at around 22 GPa, which may increase its stability at high pressures. Pressure–volume data were fit to a third order Birch–Murnaghan equation of state, resulting in V o = 816.54(9), K To = 181.54(1.39), and K' T0 = 2.76(14). These equation of state parameters are applicable to evaluating the encapsulation pressures of super-deep diamonds. The density and bulk modulus of ferromagnesian jeffbenite are similar to or higher than pyrope–almandine, pyrope–majorite, and skiagite–majorite solid solution garnets, further suggesting that jeffbenite may be an important ferric–iron silicate in the deeper parts of the mantle transition zone and uppermost lower mantle. Furthermore, future studies on the influence of temperature and oxidation state on the stability and equations of state of iron-bearing jeffbennite are still needed to determine what role, if any, jeffbenite plays in transition-zone mineralogy.

58 GEOSCIENCES↗
DOE OSTI2022

Chromosome aberrations among atomic-bomb survivors exposed in utero: updated analysis accounting for revised radiation doses and smoking

A previous study of peripheral blood lymphocyte translocations around age 40 among atomic-bomb survivors exposed in utero revealed no overall association with radiation dose—despite a clear association between translocations and dose among their mothers—but the data suggested an increase at doses below 100 mGy with a definite peak. That analysis of the in utero-exposed survivors did not adjust for their subsequent smoking behavior, an established cause of chromosomal aberrations, or their subsequent exposures to medical irradiation, a potential mediator. In addition, atomic-bomb survivor radiation dose estimates have subsequently been updated and refined. We therefore re-estimated the dose response using the latest DS02R1 dose estimates and adjusting for smoking as well as for city and proximal–distal location at the time of exposure to the atomic bomb. Sex of the survivor, mother’s age around the time of conception, and approximate trimester of gestation at the time of exposure were also considered as explanatory variables and modifiers. Precision of the estimated dose response was slightly lower due to greater variability near zero in the updated dose estimates, but there was little change in evidence of a low-dose increase and still no suggestion of an overall increase across the entire dose range. Adjustment for smoking behavior led to a decline in background number of translocations (the dose–response intercept), but smoking did not interact with dose overall (across the entire dose range). Adjustment for medical irradiation did not alter the association between dose and translocation frequency. Sex, mother’s age, and trimester were not associated with number of translocations, nor did they interact with dose overall. Interactions with dose in the low-dose range could not be evaluated because of numerical instability.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Volcanic electrification: recent advances and future perspectives

Abstract The electrification of volcanic plumes has been described intermittently since at least the time of Pliny the Younger and the 79 AD eruption of Vesuvius. Although sometimes disregarded in the past as secondary effects, recent work suggests that the electrical properties of volcanic plumes reveal intrinsic and otherwise inaccessible parameters of explosive eruptions. An increasing number of volcanic lightning studies across the last decade have shown that electrification is ubiquitous in volcanic plumes. Technological advances in engineering and numerical modelling, paired with close observation of recent eruptions and dedicated laboratory studies (shock-tube and current impulse experiments), show that charge generation and electrical activity are related to the physical, chemical, and dynamic processes underpinning the eruption itself. Refining our understanding of volcanic plume electrification will continue advancing the fundamental understanding of eruptive processes to improve volcano monitoring. Realizing this goal, however, requires an interdisciplinary approach at the intersection of volcanology, atmospheric science, atmospheric electricity, and engineering. Our paper summarizes the rapid and steady progress achieved in recent volcanic lightning research and provides a vision for future developments in this growing field.

58 GEOSCIENCES↗
DOE OSTI2023

Perspectives on Few-Body Cluster Structures in Exotic Nuclei

It is a fascinating phenomenon in nuclear physics that states with a pronounced few-body structure can emerge from the complex dynamics of many nucleons. Such halo or cluster states often appear near the boundaries of nuclear stability. As such, they are an important part of the experimental program beginning at the Facility for Rare Isotope Beams (FRIB). A concerted effort of theory and experiment is necessary both to analyze experiments involving effective few-body states, as well as to constrain and refine theories of the nuclear force in light of new data from these experiments. Here, as a contribution to exactly this effort, this paper compiles a collection of “perspectives” that emerged out of the Topical Program “Few-body cluster structures in exotic nuclei and their role in FRIB experiments” that was held at FRIB in August 2022 and brought together theorists and experimentalists working on this topic.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Ecological stoichiometry as a foundation for omics-enabled biogeochemical models of soil organic matter decomposition

Coupled biogeochemical cycles drive ecosystem ecology by influencing individual-to-community scale behaviors; yet the development of process-based models that accurately capture these dynamics remains elusive. Soil organic matter (SOM) decomposition in particular is influenced by resource stoichiometry that dictates microbial nutrient acquisition (‘ecological stoichiometry’). Despite its basis in biogeochemical modeling, ecological stoichiometry is only implicitly considered in high-resolution microbial investigations and the metabolic models they inform. State-of-science SOM decomposition models in both fields have advanced largely separately, but they agree on a need to move beyond seminal pool-based models. This presents an opportunity and a challenge to maximize the strengths of various models across different scales and environmental contexts. To address this challenge, we contend that ecological stoichiometry provides a framework for merging biogeochemical and microbiological models, as both explicitly consider substrate chemistries that are the basis of ecological stoichiometry as applied to SOM decomposition. We highlight two gaps that limit our understanding of SOM decomposition: (1) understanding how individual microorganisms alter metabolic strategies in response to substrate stoichiometry and (2) translating this knowledge to the scale of biogeochemical models. We suggest iterative information exchange to refine the objectives of high-resolution investigations and to specify limited dynamics for representation in large-scale models, resulting in a new class of omics-enabled biogeochemical models. Assimilating theoretical and modelling frameworks from different scientific domains is the next frontier in SOM decomposition modelling; advancing technologies in the context of stoichiometric theory provides a consistent framework for interpreting molecular data, and further distilling this information into tractable SOM decomposition models.

Carbon cycling↗
DOE OSTI2020

Toward indicators of the performance of US infrastructures under climate change risks

Built infrastructures are increasingly disrupted by climate-related extreme events. Being able to monitor what climate change implies for US infrastructures is of considerable importance to all levels of decision-makers. A capacity to develop cross-cutting, widely applicable indicators for more than a dozen different kinds of infrastructure, however, is severely limited at present. The development of such indicators must be considered an ongoing activity that will require expansion and refinement. A number of recent consensus reports suggest four priorities for indicators that portray the impacts of climate change, climate-related extreme events, and other driving forces on infrastructure. These are changes in the reliability of infrastructure services and the implications for costs; changes in the resilience of infrastructures to climate and other stresses; impacts due to the interdependencies of infrastructures; and ongoing adaptation in infrastructures.

adaptation, indicators, infrastructure↗
DOE OSTI2025

Projections of future climate for U.S. national assessments: past, present, future

Climate assessments consolidate our understanding of possible future climate conditions as represented by climate projections, which are largely based on the output of global climate models. Over the past 30 years, the scientific insights gained from climate projections have been refined through model structural improvements, emerging constraints on climate feedbacks, and increased computational efficiency. Within the same period, the process of assessing and evaluating information from climate projections has become more defined and targeted to inform users. As the size and audience of climate assessments has expanded, the framing, relevancy, and accessibility of projections has become increasingly important. This paper reviews the use of climate projections in national climate assessments (NCA) while highlighting challenges and opportunities that have emerged over time. Reflections and lessons learned address the continuous process to understand the broadening assessment audience and evolving user needs. Insights for future NCA development include (1) identifying benchmarks and standards for evaluating downscaled datasets, (2) expanding efforts to gather research gaps and user needs to inform how climate projections are presented in the assessment (3) providing practitioner guidance on the use, interpretation, and reporting of climate projections and uncertainty to better inform decision-making.

Assessment↗
DOE OSTI2021

Convergence of sum-up rounding schemes for cloaking problems governed by the Helmholtz equation

In this work, we consider the problem of designing a cloak for waves described by the Helmholtz equation from an integer programming point of view. The problem can be modeled as a PDE-constrained optimization problem with integer-valued control inputs that are distributed in the computational domain. A first-discretize-then-optimize approach results in a large-scale mixed-integer nonlinear program that is in general intractable because of the large number of integer variables that arise from the discretization of the domain. Instead, we propose an efficient algorithm that is able to approximate the local infima of the underlying nonconvex infinite-dimensional problem arbitrarily close without the need to solve the discretized finite-dimensional integer programs to optimality. We optimize only the continuous relaxations of the approximations for local minima and then apply the sum-up rounding methodology to obtain integer-valued controls. If the solutions of the discretized continuous relaxations converge to a local minimizer of the continuous relaxation, then the resulting discrete-valued control sequence converges weakly \(^*\) in \(L^\infty\) to the same local minimizer. These approximation properties follow under suitable refinements of the involved discretization grids. Our results use familiar concepts arising from the analytical properties of the underlying PDE and complement previous results, derived from a topology optimization point of view.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Fundamental complement of a gravitating region

Any gravitating region a in any spacetime gives rise to a generalized entanglement wedge, the hologram e(a). Holograms exhibit properties expected of fundamental operator algebras, such as strong subadditivity, nesting, and no-cloning. But the entanglement wedge EW of an AdS boundary region B with commutant $\overline{B}$ satisfies an additional condition, complementarity: EW (B) is the spacelike complement of EW ($\overline{B}$) in the bulk. Here we identify an analogue of the boundary commutant $\overline{B}$ in general spacetimes: given a gravitating region a, its fundamental complement is the smallest wedge that contains all infinite world lines contained in the spacelike complement a' of a. We refine the definition of e(a) by requiring that it be spacelike to $\tilde{a}$. We prove that e(a) is the spacelike complement of e($\tilde{a}$) when the latter is computed in a'. We exhibit many examples of and of e(a) in de Sitter, flat, and cosmological spacetimes. Here, we find that a Big Bang cosmology (spatially closed or not) is trivially reconstructible: the whole universe is the entanglement wedge of any wedge inside it. But de Sitter space is not trivially reconstructible, despite being closed. We recover the AdS/CFT prescription by proving that EW (B) = e(causal wedge of B).

Bousso, Raphael [University of California, Berkele↗
DOE OSTI2025

Robust Implicit Adaptive Low Rank Time-Stepping Methods for Matrix Differential Equations

In this work, we develop implicit rank-adaptive schemes for time-dependent matrix differential equations. The dynamic low rank approximation (DLRA) is a well-known technique to capture the dynamic low rank structure based on Dirac–Frenkel time-dependent variational principle. In recent years, it has attracted a lot of attention due to its wide applicability. Our schemes are inspired by the three-step procedure used in the rank adaptive version of the unconventional robust integrator (the so called BUG integrator) (Ceruti et al. in BIT Numer Math 62(4):1149–1174, 2022) for DLRA. First, a prediction (basis update) step is made computing the approximate column and row spaces at the next time level. Second, a Galerkin evolution step is invoked using an implicit solves for the small core matrix. Finally, a truncation is made according to a prescribed error threshold. Since the DLRA is evolving the differential equation projected on to the tangent space of the low rank manifold, the error estimate of the BUG integrator contains the tangent projection (modeling) error which cannot be easily controlled by mesh refinement. This can cause convergence issue for equations with cross terms. To address this issue, we propose a simple modification, consisting of merging the row and column spaces from the explicit step truncation method together with the BUG spaces in the prediction step. In addition, we propose an adaptive strategy where the BUG spaces are only computed if the residual for the solution obtained from the prediction space by explicit step truncation method, is too large. Here, we prove stability and estimate the local truncation error of the schemes under assumptions. We benchmark the schemes in several tests, such as anisotropic diffusion, solid body rotation and the combination of the two, to show robust convergence properties.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Analysis of Glycolate in Radioactive Waste by Ion Chromatography (IC) and Proton Nuclear Magnetic Resonance (H NMR)

Here, in preparation for implementing the Nitric-Glycolic (NG) acid flowsheet for the Savannah River Site (SRS) Liquid Waste System (LWS), analytical methods for determining glycolate at low concentration, below 20 mg/L in radioactive samples, were developed to support system management and safety. To accommodate the wide range of LWS matrix conditions, two alternative methods were developed, refined, and demonstrated for glycolate analysis in radioactive waste samples: ion chromatography (IC) and a proton nuclear magnetic resonance (H NMR). Investigators validated IC and H NMR methods for glycolate analysis, defined the range of applicability, and demonstrated key supporting analytical protocols. The deployed IC method is applicable in low to moderate ionic strength samples and requires sample pretreatment using a Dionex OnGuard II H + cartridge. The deployed H NMR method is more labor intensive but provides options for a broader range of matrices. Based on the results, high quality glycolate analysis of the Defense Waste Processing Facility (DWPF) condensate in Tank 22 is feasible by IC down to approximately 12 mg/L. Using H NMR, glycolate may be determined to 8 mg/L or lower depending on the run time with the potential for broader applicability of the method to higher ionic strength conditions in other tanks of the SRS LWS.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2023

Development, feasibility, and uncertainty of radioactive 22 Na tracer dilution and gamma spectroscopy for mass determination of molten salt for pyroprocessing spent nuclear fuels

Accurate knowledge of the total mass of molten salts both for pyroprocessing spent nuclear fuels and for molten salt reactors is necessary for safeguards purposes. However, it is challenging to know or measure the total mass of molten salt due to the complicated shapes as well as the compositional—and thus density—changes that occur over time during operation. Here, we investigated radioactive tracer dilution (RTD) as a potential safeguards technique for application to the total mass measurement of LiCl–KCl–UCl 3 salt used for uranium electrorefining, with a focus on the feasibility, uncertainty, and fission product effect of RTD for salt mass measurement. To this end, experimental-scale RTD tests (600–700 g of total mass) were initially carried out by adding 22 Na tracer salt into LiCl–KCl salt containing 0, 5, and 28wt% radioactive salts from an electrorefiner (ER) for refining uranium. Upon the completion of the experimental-scale RTD tests, we performed a scale-up RTD test by directly adding 22 Na tracer salts to the ER, which uses about 12 kg of LiCl–KCl–UCl 3 salt during normal operation. This paper reports the main features of RTD for total mass measurement of molten salt and summarizes the results of the experimental-scale and scale-up RTD tests.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

A flow-through gamma detection system for sub-surface radioxenon tracer measurements in LYNM PE1

A continuous flow-through, real-time radioxenon sensor array has been designed to monitor either 127 Xe or 133 Xe tracer gas as part of a series of underground chemical explosions in the Low-Yield Nuclear Monitoring (LYNM) Physics Experiment 1 (PE1). The sensors are pairs of shielded thallium-doped sodium iodide gamma-ray detectors that surround a gas measurement chamber. Each measurement chamber is connected to a sampling location in the geology surrounding the chemical explosion or in the access tunnel. The first experiment in the series, PE1 A, occurred in 2023 and the system quantified the 127 Xe tracer gas in the gas sampling lines for four weeks following the experiment execution. The array generated a set of radioxenon concentration time series with a 15-minute resolution that will be used to refine subsurface gas transport models.

Gamma-ray coincidence detection↗
DOE OSTI2026

Estimating the Contribution of the Nickel to Protium Loading of Full-Length Getters

The full-length getter (FLG) is a critical component of the tritium producing burnable absorber rod (TPBAR), designed to capture tritium produced within the lithium aluminate pellets. However, due to the chemical similarity between protium (¹H) and tritium (³H), the FLG readily absorbs protium, reducing its capacity to absorb tritium and increasing the risk of tritium permeation into the reactor coolant. This study evaluates protium produced from neutron irradiation of nickel plating on the FLG through 5?Ni(n,p) reactions, quantifies its contribution to the total protium observed, and informs models of tritium and hydrogen transport within TPBARs. During irradiation, neutron capture by 58Ni results in the formation of 5?Ni, which undergoes neutron bombardment to produce 4He via 5?Ni(n,a) reactions and protium via 5?Ni(n,p) reactions. The measured helium content post-irradiation provides insight into the neutron capture processes within the getter. The 4He measured within the getter may be useful in determining hydrogen produced by the nickel plating on FLG because 59Ni also produces protium in a 59Ni(n, p) reaction. The contribution of 1H from the nickel in FLG contributed less than 1% of the measured H2 gas in PIE, ranging from 0.005% to 0.614%. These findings refine current knowledge of the protium-tritium interplay in TPBARs and support the development of improved transport models for tritium and hydrogen, ultimately aiding in the optimization of TPBAR design and reactor operations.

Arbova, Dana L.↗
DOE OSTI2021

Multilevel Graph Partitioning for Three-Dimensional Discrete Fracture Network Flow Simulations

We present a topology-based method for mesh-partitioning in three-dimensional discrete fracture network (DFN) simulations that takes advantage of the intrinsic multi-level nature of a DFN. DFN models are used to simulate flow and transport through low-permeability fractured media in the subsurface by explicitly representing fractures as discrete entities. The governing equations for flow and transport are numerically integrated on computational meshes generated on the interconnected fracture networks. Modern high-fidelity DFN simulations require high-performance computing on multiple processors where performance and scalability depends partially on obtaining a high-quality partition of the mesh to balance work-loads and minimize communication across all processors. The discrete structure of a DFN naturally lends itself to various graph representations, which can be thought of as coarse-scale representations of the computational mesh. Using this concept, we develop two applications of the multilevel graph partitioning algorithm to partition the mesh of a DFN. In the first, we project a partition of the graph based on the DFN topology onto the mesh of the DFN and in the second, this DFN-based projection is used as the initial condition for further partitioning refinement of the mesh. We compare the performance of these methods with standard multi-level graph partitioning using graph-based metrics (cut, imbalance, partitioning time), computational-based metrics (FLOPS, iterations, solver time), and total run time. The DFN-based and the mesh-based partitioning methods are comparable in terms of the graph-based metrics, but the time required to obtain the partition is several orders of magnitude faster using the DFN-based partitions. The computation-based metrics show comparable performance between both methods so, in combination, the DFN-based partitions are several orders of magnitude faster than the mesh-based partition. Furthermore, the method which uses the DFN-partition solution as the initial condition of the mesh partition provided cut and imbalance values that were close to the mesh-based partition but in a fraction of the time. In turn, this hybrid method outperformed both of the other methods in terms of the total run time.

58 GEOSCIENCES↗