Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical relativity”

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

Dispersion relation for gauge-free electromagnetic drift kinetics

Recently, a new approach to gyrokinetics, invariant under electromagnetic gauge transformations, was developed. The gyrocenter equations of motion are now expressed in terms of the perturbed fields instead of the potentials, in a form suitable for numerical simulations and analytic studies. In this paper, we verify that the long-wavelength limit, i.e., the drift-kinetic limit of the new gyrokinetic theory, is in line with existing work, providing a solid foundation for simulations. Here, we compute the dispersion relation of the new drift-kinetic theory in slab geometry and find agreement with a long-wavelength limit of the full Vlasov–Maxwell model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Plant Defense Proteins: Recent Discoveries and Applications

Proteins play pivotal roles in safeguarding plants against numerous biotic and abiotic stresses. Understanding their biological functions and mechanisms of action is essential for advancing plant biology, agriculture, and biotechnology. This review considers the diversity and potential applications of plant defense proteins including pathogenesis-related (PR) proteins, chitinases, glucanases, protease inhibitors, lectins, and antimicrobial peptides. Recent advances, such as the omics technologies, have enabled the discovery of new plant defense proteins and regulatory networks that govern plant defense responses and unveiled numerous roles of plant defense proteins in stress perception, signal transduction, and immune priming. The molecular affinities and enzymatic activities of plant defense proteins are essential for their defense functions. Applications of plant defense proteins span agriculture, biotechnology, and medicine, including the development of resistant crop varieties, bio-based products, biopharmaceuticals, and functional foods. Future research directions include elucidating the structural bases of defense protein functions, exploring protein interactions with ligands and other proteins, and engineering defense proteins for enhanced efficacy. Overall, this review illuminates the significance of plant defense proteins against biotic stresses in plant biology and biotechnology, emphasizing their potential for sustainable agriculture and environmental management.

G-proteins↗
DOE OSTI2025

Numerical design considerations for vapor transport deposition of metal-halide perovskite thin films

While metal-halide perovskites (MHPs) offer high efficiency and potential application in single junction and tandem solar cells, challenges remain in translating solution-based methods to commercial-scale production. Vapor transport deposition (VTD) offers advantages such as precise control over film composition and purity, as well as the possibility of achieving high thin film deposition rates for scale-up. In this work, a numerical model is developed to simulate the deposition process, analyzing the impact of parameters including source material temperature, carrier gas and dilution gas flow rates, and chamber pressure to optimize deposition rates and film uniformity. Geometrical parameters, such as throw distance, are also found to significantly affect rate and uniformity. Furthermore, the modeling results presented here offer valuable insight for improving the scalability of VTD for the processing of metal-halide perovskites.

Chemical vapor deposition↗
DOE OSTI2021

Meteoroid Stream of Comet C/1961 T1 (Seki) and Its Relation to the December ρ-Virginids and γ-Sagittariids

We answer the question of whether long-period comet C/1961 T1 (Seki) is the parent body of a meteor shower observed in Earth’s atmosphere. We model discrete parts of a meteoroid stream associated with this comet and follow the orbital evolution of its particles via a numerical integration of their orbits. The gravitational perturbations of major planets and the Poynting–Robertson effect, which influences the motion of the particles, are taken into account. We predict two showers which are related to C/1961 T1. One of these showers was identified with the December ρ-Virginids, #502. This identification confirms an earlier suggested relationship between the shower and the comet. The second predicted shower was identified with the γ-Sagittariids, #657. The radiant areas of both showers are symmetrical in respect to the apex of Earth’s heliocentric motion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Dynamic Behavior of Natural Seep Vents: Analysis of Field and Laboratory Observations and Modeling (Final Scientific/Technical Report)

In this project, we have analyzed data collected by the U.S. Department of Energy (DOE), National Energy Technology Laboratory (NETL) in a high pressure water tunnel (HPWT) and data from two research cruises to natural seeps in the Gulf of Mexico to adapt and validate a numerical model to predict the dynamics of natural seeps in the deep oceans. The HPWT data include video observations of the shrinkage rate of individual methane and natural gas bubbles under simulated deep-water conditions. Field data were collected during two cruises by the Gulf Integrated Spill Research (GISR) Consortium led by Texas A&M University and funded by the Gulf of Mexico Research Initiative (GoMRI). These data included in situ observations from a remotely operated vehicle (ROV) of gas bubbles at two natural seep sites in the Gulf and acoustic observations of the natural seep bubble flares in the ocean water column. The acoustic data were from multibeam echosounders, one mounted in a forward-looking orientation on the ROV and another mounted down-looking in the haul of the ship. All of these laboratory and field data were focused on the dynamics of natural gas bubbles at temperatures and pressures favorable for clathrate hydrate formation between the gas and water. Our analyses of this data focused on understanding the mechanisms responsible for gas bubble dissolution within the hydrate stability zone (HSZ) of the oceans. Ice-like hydrate shells may form on the bubble-water interface under these conditions, and it was unknown how this might affect the mass transfer of gas into the ocean. We were able to extract bubble shrinkage rates from the HPWT datasets. Using this data we determined that mass transfer coefficients with and without a hydrate shell match empirical values for bubbles in contaminated systems (so-called dirty bubbles contaminated by naturally occurring surfactants). We also showed that free gas, and not gas hydrate, is the dominant dissolving phase when the hydrate sub-cooling is below 11 degree Celsius (temperature difference between hydrate the hydrate formation temperature and ambient temperator) or the pressure is reducing as bubbles rise through the ocean water column. Using this mass transfer model, our numerical model of bubble dissolution matched the over 200 HPWT experiments with an average error of 10% for predicting the bubble size at the end of an experiment. From field data in the literature, we also observed that gas bubbles dissolve faster when they are initially released, following mass transfer coefficients for so-called clean-bubbles (those not yet contaminated by surfactants). Shortly after release within the HSZ, a hydrate shell forms on the bubble-water interface, and the mass transfer reduces to rates matching those of dirty bubbles. We correlated this transition time from clean to dirty bubble behavior with the initial bubble surface area and the hydrate sub-cooling. With this model for hydrate formation time and using the mass transfer coefficients deduced from the HPWT data, we validated our numerical model for predicting the rise heights of natural seep flares in the oceans. Flare heights are commonly observed in haul-mounted acoustic multibeam data. The numerical model predicts bubbles to rise high in the ocean water column owing to the slower mass transfer rates for dirty bubbles that accompany the majority of their rise time. We found that the numerical model predictions matched the observed flare heights within 5% to 10% accuracy when we compared the rise heights of the largest bubbles released from the seafloor with the bubbles acoustically visible in the multibeam data. Bubbles become acoustically transparent as they shrink to sizes of order 1 mm in diameter for the multibeam frequencies used in the field. The forward-looking multibeam on the ROV also provided data on the lateral spreading of bubbles in natural seep flares. Our analysis of this data showed that spreading follows a diffusion process, with the effective diffusivity correlating with the wobbling length scale of these ellipsoidal bubbles. When we apply this diffusivity in a random displacement model of bubble spreading, our numerical simulations match closely the lateral spread observed by the M3 in the ocean water column. Finally, we compared the seep model predictions for the acoustic properties of these natural seep plumes with that observed by the acoustic instruments in the field. The M3 and EM 302 observations were converted to relative values of target strength using a calibration we obtained in the laboratory for the M3 and using an algorithm from the manufacturer for the EM 302. Comparing the numerical seep model to these data, we obtain good agreement over the whole height of rise of these bubble flares. This further validates the numerical model. Overall, our validated seep model captures the key dynamics of gas bubbles released from natural seeps in the oceans and helps to predict the fate of methane in the water column.

03 NATURAL GAS↗
DOE OSTI2021

Structure-preserving numerical discretizations for domains with boundaries

This SAND report documents Exploratory Express LDRD Project 223790, "Structure-preserving numerical discretizations for domains with boundaries", which developed a method to incorporate consistent treatment of domain boundaries and arbitrary boundary conditions in discrete exterior calculus (DEC) for arbitrary polygonal (2D) and tensor-product structure prism (3D) grids. The new DEC required the development of novel discrete exterior derivatives, boundary operators, wedge products and Hodge stars. This was accomplished through the use of boundary extension and the blending of known 2D operators on the interior with 1D operators on the boundary. The Hodge star was based on the Voronoi Hodge star, and retained the limitation of a triangular circumcentric primal or dual grid along with low-order accuracy. In addition to the new DEC, two related software packages were written: one for the study of DEC operators on arbitrary polygonal and polyhedral grids using both symbolic and numerical approaches and one for a (thermal) shallow water testbed using TRiSK-type numerics. Immediately relevant (already funded, through CANGA) followup work is the development of a high-order, geometrically flexible Hodge star and structure-preserving, high-order, oscillation-limiting transport operators (using WENO) for n-forms on arbitrary 2D and 3D grids. This will provide all of the machinery required for a high-order version of TRiSK with boundaries on arbitrary 2D and tensor-product 3D grids, which is applicable to both the atmospheric (CRM in E3SM-MMF) and oceanic (MPAS-O) components of E3SM.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Can lidars assess wind plant blockage in simple terrain? A WRF-LES study

Wind plant blockage reduces wind velocity upstream of wind plants, reducing the power generated by turbines adjacent to the inflow, and potentially throughout the plant as well. The nature of the mechanism that amplifies blockage as well as the velocity reductions in both the induction zone and potentially deeper into the array are not well understood. Field observations can provide valuable insight into the characteristics of the induction zone and the mechanisms that amplify it. However, the relatively small velocity reductions that have been measured experimentally pose a challenge in quantifying blockage, especially in onshore environments with flow heterogeneities that may be of the same scale as the blockage effect itself. We simulate the flow around the King Plains wind plant in the relatively simple terrain of Oklahoma, the location of the American WAKE experimeNt, to evaluate wind plant blockage in this environment. Using numerical simulations, we find the largest velocity deceleration (0.64 m s−1; 8%) immediately upstream of the wind plant, and 1% velocity deficits 24 rotor diameters upstream of the first turbine row. We also use virtual measurements upstream of the wind plant to analyze the uncertainties and difficulties in measuring blockage using a scanning lidar on shore. Based on our virtual lidar study, the induction zone of land-based wind plants can be incorrectly estimated using observations if the effects of nonuniform terrain on the flow are not carefully considered. Changes in terrain elevation produce local variations in wind speed (as measured by a scanning lidar) that exceed in magnitude the deceleration within the induction zone. We refer to these local changes in wind speed as terrain effects. A methodology to differentiate between terrain effects and blockage in experimental settings is proposed and evaluated herein, highlighting the difficulties and uncertainties associated with measurement and simulation of blockage in even relatively simple onshore environments.

17 WIND ENERGY↗
DOE OSTI2020

Evaluation of Methods for Causal Discovery in Hydrometeorological Systems

Understanding causal relations is of utmost importance in hydrology and climate research for systems identification, prediction, and understanding systems behavior in a changing climate. Traditionally, researchers in hydrometeorology attempted to study causal questions by conducting controlled experiments using numerical models. This approach, however, in most cases of interest provides uncertain results because the models are approximate representation of the natural system. An alternative approach that has recently drawn significant attention in several fields is to infer causal relations from purely observational data. It possesses several traits to its utility particularly in hydrometeorology due to the rapid accumulation of in situ and remotely sensed data records. The first objective of this study is to present a brief description of four causal discovery methods (Granger causality, Transfer Entropy, graph-based algorithms, and Convergent Cross Mapping) with special emphasis on the assumptions on which they are built. Second, using synthetic data generated from a hydrological model, we assess their performance in retrieving causal information taking into account sensitivity to sample size and presence of noise. Last, we use causal analysis to examine and formulate hypotheses on causal drivers of evapotranspiration in a shrubland region during summer and winter seasons. An interpretation of the hypotheses based on canopy seasonal dynamics and evapotranspiration processes is presented. It is hoped that the results presented here can be useful in guiding researchers studying hydrometeorological systems as to which causal method is most appropriate to the characteristics of the system under study.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

The First Tri-Lab Workshop on Formal Verification: Capabilities, Challenges, Research Opportunities, and Exemplars

The First Tri-Lab Workshop on Formal Verification was held in Santa Fe, New Mexico, on December 5th, 2023. This workshop gathered staff from Sandia, Los Alamos, and Lawrence Livermore National Laboratories and NASA’s Jet Propulsion Laboratory. This report summarizes and expands on the presentations given and discussion had at this workshop. In this report, we describe the current capabilities and research needs related to formal methods at the NNSA labs. In particular, we identify medium-term and long-term research gaps in programming languages, formalization efforts of complex systems, embedded systems verification, hardware verification, cybersecurity, formal methods usability, workflows, numerical methods, the use of formal methods for artificial intelligence (and its converse, artificial intelligence for formal methods), and collaboration opportunities and considerations on these topics. We conclude with a small number of exemplar research problems related to these topics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Changes in the Regional Water Cycle and Their Impact on Societies

ABSTRACT Changes in “blue water”, which is the total supply of fresh water available for human extraction over land, are quite closely related to changes in runoff or equivalently precipitation minus evaporation, . This article examines how climate change‐driven recent past and future changes in the regional water cycle relate to blue water availability and changes in human blue water demand. Although at the largest scales theoretical and numerical model predictions are in broad agreement with observations, at continental scales and below models predict large ranges of possible future and runoff especially at the scale of individual river catchments and for shorter timescale subseasonal floods and droughts. Nevertheless, it is expected that the occurrence and severity of floods will increase and that of droughts may increase, possibly compounded by human‐driven non‐climatic changes such as changes in land use, dam water impoundment, irrigation and extraction of groundwater. Contemporary assessments predict that increases in 21st century human water extraction in many highly‐populated regions are unlikely to be sustainable given projections of future . To reduce uncertainty in future predictions, there is an urgent need to improve modeling of atmospheric, land surface and human processes and how these components are coupled. This should be supported by maintaining the observing network and expanding it to improve measurements of land surface, oceanic and atmospheric variables. This includes the development of satellite observations stable over multiple decades and suitable for building reanalysis datasets appropriate for model evaluation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Investigation of inter-subject variation in ultrafine particle deposition across human nasal airways: A study involving children, adults, and the elderly

Occupational and environmental exposure to toxic nanoparticles, driven by the rapid expansion of nanotechnology, raises significant respiratory health concern. Numerous studies have explored airflow and particle dynamics in adult nasal airways, but understanding the impact of age-related anatomical changes in children and the elderly remains limited. This study systematically investigates age-related anatomical variations and associated influence on nasal airflow dynamics and ultrafine particle deposition characteristics. Using Computational Fluid-Particle Dynamics (CFPD) method, simulation was conducted under diverse inhalation conditions spanning a wide age range, including: two children (5 years old), two young adults (in their twenties), and two elderly (over 77 years old). Our results reveal distinctive variations across age groups in anatomical dimensions, which affect distribution of wall shear stress where the elderly and children display unique patterns distinct from the young adults. While total deposition efficiency differs significantly between children and adults, filtration efficiency in the subregion with most deposition, main respiratory, remains consistent. However, inter-subject differences are observed in the vestibular and olfactory regions,emphasizing nuanced impact of age-related anatomical variations. Overall and subregional empirical equations for deposition efficiency were developed by incorporating the combined diffusion parameter, Sc a Δ b , corroborating the use of geometrical characteristic parameters for each specific subject in predicting nasal deposition efficiency across age groups. Our findings contribute to predictive nanoparticle exposure analysis in nasal airways across different age groups, thereby enhancing respiratory healthcare for individuals across the life span.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Phenolic acid-degrading Paraburkholderia prime decomposition in forest soil

Abstract Plant-derived phenolic acids are catabolized by soil microorganisms whose activity may enhance the decomposition of soil organic carbon (SOC). We characterized whether phenolic acid-degrading bacteria enhance SOC mineralization in forest soils when primed with 13C-labeled p-hydroxybenzoic acid (pHB). We further tested whether pHB-induced priming could explain differences in SOC content among mono-specific tree plantations in a 70-year-old common garden experiment. pHB addition primed significant losses of SOC (3–13 µmols C g−1 dry wt soil over 7 days) compared to glucose, which reduced mineralization (-3 to -8 µmols C g−1 dry wt soil over 7 days). The principal degraders of pHB were Paraburkholderia and Caballeronia in all plantations regardless of tree species or soil type, with one predominant phylotype (RP11ASV) enriched 23-fold following peak pHB respiration. We isolated and confirmed the phenolic degrading activity of a strain matching this phylotype (RP11T), which encoded numerous oxidative enzymes, including secretion signal-bearing laccase, Dyp-type peroxidase and aryl-alcohol oxidase. Increased relative abundance of RP11ASV corresponded with higher pHB respiration and expression of pHB monooxygenase (pobA), which was inversely proportional to SOC content among plantations. pobA expression proved a responsive measure of priming activity. We found that stimulating phenolic-acid degrading bacteria can prime decomposition and that this activity, corresponding with differences in tree species, is a potential mechanism in SOC cycling in forests. Overall, this study highlights the ecology and function of Paraburkholderia whose associations with plant roots and capacity to degrade phenolics suggest a role for specialized bacteria in the priming effect.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A resource adequacy assessment of correlated wide-area outages in the power grid

As the power grid is undergoing rapid transformations, numerous questions are emerging about its vulnerability to wide-area extreme events (WAEE), which could influence its operations. Relatively few analyses have been conducted to date regarding the impact of correlated outages during WAEEs on the grid’s ability to balance resources with demand. This study addresses this gap by conducting a resource adequacy analysis for a hurricane-inspired WAEE on a 2035 synthetic power grid system for the United States. A sensitivity analysis was also conducted to characterize the relative impact of weather on unserved energy. Our results indicate that although the magnitude and duration of the shortfalls vary depending on weather conditions, persistent shortfalls are observed in some regions. Initial explorations indicate a strong correlation between transmission-constrained regions and regions with persistent shortfalls. Future work could generate empirically-grounded representations for generator outages as well as conduct causal analyses of these shortfalls to improve understanding of drivers as well as possible mitigation strategies. Continued exploration of extreme weather impacts on the grid is important to develop more robust understanding of the reliability and resilience of our power systems, especially as they undergo rapid transformations.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

A phenomenological wobbling model for isolated pulsars and the braking index

ABSTRACT An isolated pulsar is a rotating neutron star possessing a very high magnetic dipole moment, thus providing a powerful radiating mechanism. These stars loose rotational energy E through various processes, including a plasma wind originating from a highly magnetized magnetosphere and the emission of magnetic dipole radiation (MDR). Such phenomena produce a time decreasing angular velocity Ω(t) of the pulsar that is usually quantified in terms of its braking index. Although these mechanisms are widely acknowledged as the primary drivers of the spin evolution of isolated pulsars, it is plausible that other contributing factors influencing this effect have yet to be comprehensively investigated. Most of young isolated pulsars present a braking index different from that given by the MDR and plasma wind processes. Working in the weak field (Newtonian) limit, we take in this work a step forward in describing the evolution of such a system by allowing the star’s shape to wobble around an ellipsoidal configuration as a backreaction effect produced by its rotational deceleration. It is assumed that an internal damping of the oscillations occurs, thus introducing another form of energy loss in the system, and this phenomenon may be related to the deviation of the braking index from the models based on $\dot{E} \sim -\Omega ^4$ predictions. Numerical calculations suggest that the average braking index for typical isolated pulsars can be thus explained by a simple phenomenological model.

Araujo, E. C. A. (ORCID:0000000213520879)↗
DOE OSTI2025

Connecting GRBs from binary neutron star mergers to nuclear properties of neutron stars

The fate of the binary neutron star (NS) merger remnants hinges sensitively upon the NS equation of state and the threshold mass, M ls , that separates a long-lived from a short-lived NS remnant. The nature of the electromagnetic counterparts is also influenced by the remnant type, particularly in determining whether a gamma-ray burst from a compact binary merger (cbGRB) is of short or long duration. We propose a novel approach to probe Mls by linking it to the estimated observed ratio of long to short cbGRBs. We find that current observations broadly favor a relatively high value for this transition, M ls ≃ 1.3M TOV , for which M TOV ≲ 2.6M ⊙ , consistent with numerical simulations, as also shown here. Furthermore, our results disfavor nuclear physics scenarios that would lead to catastrophic pressure loss at a few times nuclear density and temperatures of tens of MeV, leading to a rapid gravitational collapse of binaries with total mass M ≲ 1.3M TOV . Future individual gravitational wave events with on-axis cbGRBs can further bound Mls.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Large-Scale Optimization with Linear Equality Constraints Using Reduced Compact Representation

For optimization problems with linear equality constraints, we prove that the (1,1) block of the inverse KKT matrix remains unchanged when projected onto the nullspace of the constraint matrix. In this work, we develop reduced compact representations of the limited-memory inverse BFGS Hessian to compute search directions efficiently when the constraint Jacobian is sparse. Orthogonal projections are implemented by a sparse QR factorization or a preconditioned LSQR iteration. In numerical experiments two proposed trust-region algorithms improve in computation times, often significantly, compared to previous implementations of related algorithms and compared to IPOPT.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Acute stress reduces population-level metabolic and proteomic variation

Abstract Background Variation in omics data due to intrinsic biological stochasticity is often viewed as a challenging and undesirable feature of complex systems analyses. In fact, numerous statistical methods are utilized to minimize the variation among biological replicates. Results We demonstrate that the common statistics relative standard deviation (RSD) and coefficient of variation (CV), which are often used for quality control or part of a larger pipeline in omics analyses, can also be used as a metric of a physiological stress response. Using an approach we term Replicate Variation Analysis (RVA), we demonstrate that acute physiological stress leads to feature-wide canalization of CV profiles of metabolomes and proteomes across biological replicates. Canalization is the repression of variation between replicates, which increases phenotypic similarity. Multiple in-house mass spectrometry omics datasets in addition to publicly available data were analyzed to assess changes in CV profiles in plants, animals, and microorganisms. In addition, proteomics data sets were evaluated utilizing RVA to identify functionality of reduced CV proteins. Conclusions RVA provides a foundation for understanding omics level shifts that occur in response to cellular stress. This approach to data analysis helps characterize stress response and recovery, and could be deployed to detect populations under stress, monitor health status, and conduct environmental monitoring.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Impact of the optical parametric amplification phase on laser pulse compression

Optical parametric chirped-pulse amplification (OPCPA) is an effective way to generate ultrashort pulses that has been used extensively for a variety of applications requiring high peak intensities. Precise control and measurement of a system’s spectral and spatial phases are required for Fourier transform–limited pulse compression and diffraction-limited focusing. Phase accumulated during optical parametric amplification (OPA) can degrade the compressibility and focusability of the pulse, reducing peak intensity. We used analytic and numerical analysis of OPA to study the influence of crystal parameters, wavefront of the pump and signal, and their relative optical alignment on the accumulated phase. Here, we show that the accumulated signal phase is largely independent of amplifier saturation and, with significant local wavefront gradients in the signal or pump beam, the quality of the compressed pulses can be degraded. We use first-order expressions for the pump- and signal-angle sensitivity to evaluate an OPCPA system consisting of a highly deuterated potassium dihydrogen phosphate (DKDP) amplifier designed to support bandwidth for 15 fs pulses centered at 920 nm.

47 OTHER INSTRUMENTATION↗