Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “random effects”

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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Interface models for three-dimensional Rayleigh–Taylor instability

We derive interface models for three-dimensional Rayleigh–Taylor instability (RTI), making use of a novel asymptotic expansion in the non-locality of the fluid flow. These interface models are derived for the purpose of studying universal features associated with RTI such as the Froude number in single-mode RTI, the predicted quadratic growth of the interface amplitude under multi-mode random perturbations, the optimal (viscous) mixing rates induced by the RTI and the self-similarity of horizontally averaged density profiles and the remarkable stabilization of the mixing layer growth rate which arises for the three-fluid two-interface heavy–light–heavy configuration, in which the addition of a third fluid bulk slows the growth of the mixing layer to a linear rate. Our interface models can capture the formation of small-scale structures induced by severe interface roll-up, reproduce experimental data in a number of different regimes and study the effects of multiple interface interactions even as the interface separation distance becomes exceedingly small. Compared with traditional numerical schemes used to study such phenomena, our models provide a computational speed-up of at least two orders of magnitude.

Mechanics↗
DOE OSTI2025

Computing Nonequilibrium Responses with Score-Shifted Stochastic Differential Equations

Using equilibrium fluctuations to understand the response of a physical system to an externally imposed perturbation is the basis for linear response theory, which is widely used to interpret experiments and shed light on microscopic dynamics. For nonequilibrium systems, perturbations cannot be interpreted simply by monitoring fluctuations in a conjugate observable and general response results rely on path ensemble averaging. Furthermore, these techniques do not apply to perturbations that affect the diffusion tensor in a stochastic system. Here, we introduce an “effective” physical process that represents the diffusion perturbed dynamics and enables accurate calculations of responses to a change in the diffusion. Interestingly, the effective dynamics contain an additional drift involving the instantaneous “score” of the system, and we leverage score matching algorithms to carry out nonequilibrium response calculations on systems for which the exact stationary distribution is unknown.

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

Quantum Control Landscapes Beyond the Dipole Approximation: Controllability, Singular Controls, and Resources

We investigate the control landscapes of closed n-level quantum systems beyond the dipole approximation by including a polarizability term in the Hamiltonian. The latter term is quadratic in the control field. Theoretical analysis of singular controls is presented, which are candidates for producing landscape traps. The results for considering the presence of singular controls are compared to their counterparts in the dipole approximation (i.e., without polarizability). A numerical analysis of the existence of traps in control landscapes for generating unitary transformations beyond the dipole approximation is made upon including the polarizability term. An extensive exploration of these control landscapes is achieved by creating many random Hamiltonians which include terms linear and quadratic in a single control field. The discovered singular controls are all found not to be local optima. This result extends a great body of recent work on typical landscapes of quantum systems where the dipole approximation is made. We further investigate the relationship between the magnitude of the polarizability and the fluence of the control resulting from optimization. It is also shown that including a polarizability term in an otherwise uncontrollable dipole coupled system removes traps from the corresponding control landscape by restoring controllability. We numerically assess the effect of a polarizability term on a known example of a particular three-level Λ-system with a second order trap in its control landscape. It is found that the addition of the polarizability removes the trap from the landscape. The general practical control implications of these simulations are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Countercurrent flow characteristics of next generation solvent in novel 3D printed columns for carbon capture

Solvent based absorption for carbon dioxide capture in a packed column is being considered as an efficient technology for the decarbonization process of point source. Significant efforts are underway to improve the design of packings with the cutting edge 3D printing technology for efficient carbon capture. Accordingly, multiphase flow studies were conducted to assess the performance of novel 3D printed columns with various triply periodic minimal surface (TPMS) designs. The effects of solvent properties, liquid and gas loads on the performance of TPMS columns are extensively explored. Hydrodynamics of the potential water-lean solvents (EEMPA) as well as aqueous monoethanolamine (MEA) solvent for carbon capture are evaluated and compared. The interfacial area and liquid holdup increase with increasing liquid loads (q L ) for TPMS columns. Schwarz column consistently shows the highest liquid holdup value. The EEMPA exhibits higher values for the interfacial area (10–15 %) and liquid holdup (~4 %) in comparison to the MEA. TPMS columns exhibit the intermediate value of dry pressure drop between the random and the structure packed columns. Among selected TPMS columns, the gyroid packing shows the lowest pressure drop. The gas load has marginal impact on the interfacial area at lower value while a higher gas load leads to column’s flooding. Prior to the flooding, the interfacial area in TMPS packings rises with increased gas load at a fixed liquid load except Schwarz packing where interfacial area is incentive to the gas load. Additionally, the liquid holdup and wet pressure drop rise as gas load increases in TPMS columns. Overall, Fisher Koch packing is more susceptible to flood as compared to other TPMS packings. Flow regimes: loading and flooding are also delineated in the TPMS packings.

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

Deep Learning for Intelligent Bubble Size Detection in the Spallation Neutron Source Visual Target

The Spallation Neutron Source (SNS) at Oak Ridge National Laboratory (ORNL) will undergo proton power upgrade (PPU), increasing the proton beam power from 1.4 MW to 2.8 MW. From 2.8 MW, 2.0 MW will go to the current First Target Station and the rest will go to the future Second Target Station (STS). The First Target Station uses a liquid mercury target that is contained in a 316L stainless steel vessel. The proton beam is pulsed at 60 Hz, with a pulse of about 0.7μs. When the proton beam hits the target, the intense energy deposition leads to a rapid rise in temperature in the mercury. This temperature rise creates pressure waves that propagate through the mercury and cause cavitation erosion. The power upgrade will cause stronger pressure waves that will further increase damage because of cavitation. Injecting small helium bubbles in the mercury has been an efficient method of mitigating the pressure wave at 1.4 MW. However, at higher power, additional mitigation is necessary. Therefore, the 2 MW target vessel will be equipped with swirl bubblers and an additional gas injection port near the nose to inject more gas in the target. To develop a gas injection strategy and design, flow visualization in water with a transparent prototypical target (“visual target”) was performed. Bubble sizes and their spatial distribution in the flow loop are crucial to understanding the effectiveness of the bubbles in mitigating pressure waves. Bubbles were generated in the visual target under varied conditions of input pressures with helium and air. Images were captured using a high-speed camera at varied frame rates at different positions away from the swirl bubbler and different depths in the flow loop under varying lighting conditions. Initially, methods such as circular Hough transforms were applied after a series of images processing to obtain a general distribution of bubble sizes. Bubbles smaller than 500 μm are preferred to effectively mitigate the effect of pressure waves, which demands an accurate bubble detection and sizing system. Intelligent detection and identification of bubble sizes alleviate misdetection and improves accuracies. Employing neural networks, intelligent detection of bubble sizes and their distribution was developed and provides a robust alternative to traditional techniques. Human intervention was employed to label in-focus and out-of-focus bubbles in the set of training images. An object detection network using a pretrained convolutional neural network was created that extracted the features from the training images. Data augmentation was used to improve network accuracy through a random transformation of the original data.

Rasheed, Fayaz↗
DOE OSTI2023

How can a diverse set of integral and semi-integral measurements inform identification of discrepant nuclear data?

Nuclear data are used for a variety of applications, including criticality safety, reactor performance, and material safeguards. Despite the breadth of use-cases, the effective neutron multiplication factor, keff, of ICSBEP critical assemblies are primarily used for nuclear data validation; these are sensitive to specific energy regions and nuclides and are unable to uniquely constrain nuclear data. As a consequence, general-purpose nuclear data libraries, such as ENDF/B-VIII.0, may have deficiencies that, while not apparent in criticality applications, negatively impact other applications, such as non-destructive analysis of special nuclear material and neutron diagnosed subcritical experiments. Recent work by the Experiments Underpinned by Computational Learning for Improvements in Nuclear Data (EUCLID) project developed a machine learning tool, RAFIEKI, which uses random forests and the SHAP metric to determine which nuclear data contribute most to predicted bias between measured and simulated responses (e.g. keff). This paper contrasts RAFIEKI analysis applied to keff only against RAFIEKI analysis with keff paired with either LLNL pulsed sphere measurements or subcritical benchmarks. Two examples show that a) including pulsed sphere measurements substantially increases 9Be nuclear data importance to bias between 2 and 15 MeV, and b) including subcritical benchmarks has the potential for disentangling compensating errors between 240Pu (n,el) and (n,il) cross-sections between 0.1 and 10 MeV. These results show that RAFIEKI analysis applied to response sets that include, but go beyond, keff can aid nuclear data evaluators in identifying issues in nuclear data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

A Radiative Transfer Model for Patchy Landscapes Based on Stochastic Radiative Transfer Theory

The availability of global high-resolution land cover maps provides promising a priori knowledge for characterizing subpixel heterogeneity and improving predictions of directional reflectance of coarse-resolution pixels. Due to mutual shadowing and sheltering effects between the adjacent forest and cropland patches, the spectral nonlinear mixing of patchy ecotones is significant, especially when the sun illuminates the ecotone from the forest side with high solar zenith angle. The spectral linear mixture (SLM) approach leads to overestimation of the bidirectional reflectance factor (BRF) in the red band in the principal plane (PP), with a maximum absolute error (MAE) of 0.0063 and a maximum relative error (MRE) of 52.5%, and to underestimation in the near-infrared band in PP with an MAE of 0.0940 and an MRE of 14.5%. In a scenario with randomly distributed boundary orientations, the overestimation of SLM increases with the degree of fragmentation and the view zenith angle. We propose a Radiative Transfer model for patchy ECotones (RTEC). which improves R² from 0.61 to 0.94 in the red band of Landsat-8 directional reflectance at the validation site. The RTEC model provides an efficient and analytical approach for directional reflectance predictions over heterogeneous patchy landscapes at coarse resolution and will be used for biophysical parameter retrievals [e.g., the leaf area index (LAI)] in future applications.

Zeng, Yelu↗
DOE OSTI2022

SDSS-IV MaNGA: The MaNGA Dwarf Galaxy Sample Presentation

We present the MaNGA Dwarf galaxy (MaNDala) Value Added Catalog (VAC), from the final release of the Sloan Digital Sky Survey-IV program. MaNDala consists of 136 randomly selected bright dwarf galaxies with M * < 10 9.1 M ⊙ and M g > -18.5, making it the largest integral field spectroscopy homogeneous sample of dwarf galaxies. We release a photometric analysis of the g, r, and z broadband imaging based on the DESI Legacy Imaging Surveys, as well as a spectroscopic analysis based on the Pipe3D SDSS-IV VAC. Our release includes the surface brightness (SB), geometric parameters, and color profiles, Sérsic fits as well as stellar population properties (such as stellar ages, metallicities, and star formation histories), and emission lines' fluxes within the FOV and the effective radii of the galaxies. We find that the majority of the MaNDala galaxies are star-forming late-type galaxies with $\langle$n Sersic,r $\rangle$ ~ 1.6 that are centrals (central/satellite dichotomy). MaNDala covers a large range of SB values (we find 11 candidate ultra-diffuse galaxies and three compact ones), filling the gap between classical dwarfs and low-mass galaxies in the Kormendy Diagram and in the size–mass/luminosity relation, which seems to flatten at 10 8 < M * /M ⊙ < 10 9 with $\langle$R e,r $\rangle$ ~ 2.7 kpc. A large fraction of MaNDala galaxies formed from an early low-metallicity burst of SF, but also from late SF events from more metal-enriched gas: half of the MaNDala galaxies assembled 50% of their mass at $\langle$z$\rangle$ > 2, while the last 20% was at $\langle$z$\rangle$ < 0.3. Finally, a bending of the sSFR-M * relation at M * ~ 10 9 M ⊙ for the main-sequence galaxies seems to be supported by MaNDala.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Realization of Non‐Equilibrium Wurtzite Structure in Heterovalent Ternary MgSiN 2 Film Grown by Reactive Sputtering

The piezoelectric and ferroelectric applications of heterovalent ternary materials are not well explored. Epitaxial MgSiN 2 films are grown at 600 °C on (111)Pt//(001)Al 2 O 3 substrates by the reactive sputtering method using metallic Mg and Si under the N 2 atmosphere. Detailed X-ray diffraction measurements and transmission electron microscopy observations revealed that the epitaxially grown films on the substrates have a hexagonal wurtzite structure with c-axis out-of-plane orientation. The random occupation of this structure by Mg and Si differs from that of the previously reported structure in which these two cations periodically occupy the cationic sites. However, the lattice spacings closely approximate those that are previously reported, irrespective of the ordering, and they are almost comparable with those of (Al 0.8 Sc 0.2 )N. The wide bandgap of >5.0 eV in deposited MgSiN 2 is compatible with that of AlN and suggests durability against the application of strong external electric fields, possibly to induce polarization switching. In addition, MgSiN 2 is shown to have piezoelectric properties with an effective d 33 value of 2.3 pm V −1 for the first time. This work demonstrates the compositional expansion of hexagonal wurtzite to heterovalent ternary nitrides for novel piezoelectric materials, whose ferroelectricity is expected.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Highly Oriented Organic Ferroelectric Films with Single-Crystal-Level Properties from Restrained Crystallization

Device applications often require thin film growth with quasi-two-dimensional morphology and crystallization which are not always compatible. In this work, we exploit the method of low-temperature deposition followed by restrained crystallization (LDRC) in the growth of crystalline organic molecular thin films, which tend to grow three-dimensionally with random orientations. We demonstrate that, for 2-methylbenzimidazole (MBI), a molecular ferroelectric that tends to crystalize in spherulites, the quasi-two-dimensional films can be grown using LDRC with highly oriented polar axes and single-crystal-level ferroelectric properties. The crystallization process was shown to occur during the post-deposition annealing process using the in situ electrical measurements. The limited diffusion, low nucleation density, and low activation energy was found critical for the formation of the plate-shaped quasi-two-dimensional films. These results mark an important step in elucidating the LDRC as an effective and general approach for fabricating films with balanced crystallinity and morphology which are critical for applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

DDStore: Distributed Data Store for Scalable Training of Graph Neural Networks on Large Atomistic Modeling Datasets

Graph neural networks (GNNs) are a class of Deep Learning models used in designing atomistic materials for effective screening of large chemical spaces. To ensure robust prediction, GNN models must be trained on large volumes of atomistic data on leadership class supercomputers. Even with the advent of modern architectures that consist of multiple storage layers that include node-local NVMe devices in addition to device memory for caching large datasets, extreme-scale model training faces I/O challenges at scale.We present DDStore, an in-memory distributed data store designed for GNN training on large-scale graph data. DDStore provides a hierarchical, distributed, data caching technique that combines data chunking, replication, low-latency random access, and high throughput communication. DDStore achieves near-linear scaling for training a GNN model using up to 1000 GPUs on the Summit and Perlmutter supercomputers, and reaches up to a 6.15x reduction in GNN training time compared to state-of-the-art methodologies.

Choi, Jong Youl↗
DOE OSTI2025

Floquet–Bloch manipulation of the Dirac gap in a topological antiferromagnet

Floquet–Bloch manipulation, achieved by driving a material periodically with a laser pulse, is a method that enables the engineering of electronic and magnetic phases in solids by effectively modifying the structure of their electronic bands. However, the application of Floquet–Bloch manipulation in topological magnetic systems, particularly those with inherent disorder, remains largely unexplored. Here, in this study, we realize Floquet–Bloch manipulation of the Dirac surface-state mass of the topological antiferromagnet MnBi 2 Te 4 . Using time- and angle-resolved photoemission spectroscopy, we show that opposite helicities of mid-infrared circularly polarized light result in substantially different Dirac mass gaps in the antiferromagnetic phase, despite the equilibrium Dirac cone being massless. We explain our findings in terms of a Dirac fermion with a random mass. Our results underscore Floquet–Bloch manipulation as a powerful tool for controlling topology, even in the presence of disorder, and for uncovering properties of materials that may elude conventional probes.

42 ENGINEERING↗
DOE OSTI2024

Geographical Insights into Suicide Mortality Through Spatial Machine Learning

Suicide mortality is a leading cause of death in the United States, with an upward trend that emphasizes its significance as a public health issue. Previous research has employed global models like ordinary least squares (OLS) regression and local models such as geographically weighted regression (GWR). While local models are useful for analyzing spatial variations in suicide mortality, they share limitations with traditional global models, particularly about their inability to handle multi-collinearity and non-linear relationships. Machine learning approaches, like random forests (RF), can address some of these limitations but often fail to account for spatial variability. This gap highlights the need for spatial ML models specifically designed to tackle suicide mortality. This research seeks to fill this void by using a geographically weighted random forest model (GWRF) to examine the associations between county-level suicide mortality in the U.S. from 2010 to 2020 and various social and environmental determinants of health. A key aspect of our methodology is disciplined feature selection, which reduces the pool of explanatory variables by about 90%. This refinement enhances the explanatory power of both global (R2 improved from 0.59 to 0.67) and local (R2 improved from 0.64 to 0.67) RF models while reducing their run times. An analysis of the importance scores for these selected features reveals that the drivers of suicide mortality vary by context. Thus, to effectively address regional disparities and inform targeted public health interventions, a holistic approach that incorporates multiple county-level characteristics is essential.

Lebakula, Viswadeep [ORNL] (ORCID:0000000152935914↗
DOE OSTI2020

Survey of Image Denoising Methods for Medical Image Classification

Medical imaging devices, such as X-ray machines, inherently produce images that suffer from visual noise. Our objectives were to (i.) determine the effect of image denoising on a medical image classification task, and (ii.) determine if there exists a correlation between image denoising performance and medical image classification performance. We performed the medical image classification task on chest X-rays using the DenseNet-121 convolutional neural network (CNN) and used the peak signal-to-noise ratio (PSNR) and structural similarity (SSIM) metrics as the image denoising performance measures. We first found that different denoising methods can make a statistically significant difference in classification performance for select labels. We also found that denoising methods affect fine-tuned models more than randomly-initialized models and that fine-tuned models have significantly higher and more uniform performance than randomly-initialized models. Lastly, we found that there is no significant correlation between PSNR and SSIM values and classification performance for our task.

Michael, Peter↗
DOE OSTI2022

Augmented Gaussian random field: Theory and computation

We propose the novel augmented Gaussian random field (AGRF), which is a universal framework incorporating the data of observable and derivatives of any order. Rigorous theory is established. We prove that under certain conditions, the observable and its derivatives of any order are governed by a single Gaussian random field, which is the aforementioned AGRF. As a corollary, the statement "the derivative of a Gaussian process remains a Gaussian process" is validated, since the derivative is represented by a part of the AGRF. Moreover, a computational method corresponding to the universal AGRF framework is constructed. Both noiseless and noisy scenarios are considered. Formulas of the posterior distributions are deduced in a nice closed form. A significant advantage of our computational method is that the universal AGRF framework provides a natural way to incorporate arbitrary order derivatives and deal with missing data. We use four numerical examples to demonstrate the effectiveness of the computational method. The numerical examples are composite function, damped harmonic oscillator, Korteweg-De Vries equation, and Burgers' equation.

Mathematics↗
DOE OSTI2026

Telomere-to-telomere assemblies of chromosome 10 reveal complex adaptive variation of 3-ketoacyl-CoA-synthases in Populus trichocarpa likely driven by Helitrons

The model woody plant Populus trichocarpa displays an atypical alkene-diverse wax cuticle likely driven by copy number variation (CNV) of 3-ketoacyl-CoA synthases ( KCS ), which has been difficult to confirm with short-read assemblies. Long-read sequencing enables the development of telomere-to-telomere resources to detect cryptic variation, including CNVs, which are currently missed. Integrating this information can improve genomic prediction for breeding and provide insights into the evolutionary basis of important traits. Our analysis of 78 long-read haplotypes from chromosome 10 identified more than twice as many KCS genes as previously reported, and numerous intragenic non-synonymous substitutions. Random Forest predictive models highlighted the importance of Potri.010G079500 in producing very long chain alkenes; however, its absence did not predict previously reported alkene-deficient phenotypes. Instead, alkene levels are best predicted by the combinations of KCS copies. Additionally, amino acid substitutions clustered around ligand and donor binding pockets, suggesting they contribute to differing wax cuticle composition. Finally, each KCS gene and copy was linked to a Helitron transposon. A phylogenetic analysis suggests Helitrons are the evolutionary mechanism for generating KCS tandem arrays. Long-read generated telomere-to-telomere assemblies of P. trichocarpa chromosome 10 revealed large-effect loci critical to genetic studies that are unattainable from short-reads. This new resource produced novel insights into genome structure and function, and a novel mechanism for generating tandem gene duplication. Our results highlight that, given current challenges in annotation and assembly, detailed and focused long-read sequences are key to interpreting complex genomic regions that contain tandem copy number variants.

09 BIOMASS FUELS↗
DOE OSTI2024

Ensemble variational Fokker-Planck methods for data assimilation

Particle flow filters solve Bayesian inference problems by smoothly transforming a set of particles into samples from the posterior distribution. Particles move in state space under the flow of an McKean-Vlasov-Itˆo process. This work introduces the Variational Fokker-Planck (VFP) framework for data assimilation, a general approach that includes previously known particle flow filters as special cases. The McKean-Vlasov-Itˆo process that transforms particles is defined via an optimal drift that depends on the selected diffusion term. It is established that the underlying probability density - sampled by the ensemble of particles - converges to the Bayesian posterior probability density. For a finite number of particles the optimal drift contains a regularization term that nudges particles toward becoming independent random variables. Based on this analysis, we derive computationally-feasible approximate regularization approaches that penalize the mutual information between pairs of particles, and avoid particle collapse. Moreover, the diffusion plays a role akin to a particle rejuvenation approach that aims to alleviate particle collapse. The VFP framework is very flexible. Different assumptions on prior and intermediate probability distributions can be used to implement the optimal drift, and localization and covariance shrinkage can be applied to alleviate the curse of dimensionality. A robust implicit-explicit method is discussed for the efficient integration of stiff McKean- Vlasov-Itˆo processes. Here, the effectiveness of the VFP framework is demonstrated on three progressively more challenging test problems, namely the Lorenz ’63, Lorenz ’96 and the quasi-geostrophic equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Constructing a spatiotemporally coherent long-term PM2.5 concentration dataset over China during 1980–2019 using a machine learning approach

The lack of long-term observations and satellite retrievals of health-damaging fine particulate matter in China has demanded the estimates of historical PM 2.5 (particulate matter less than 2.5 µm in diameter) concentrations. This study constructs a gridded near-surface PM 2.5 concentration dataset across China covering 1980–2019 using the space-time random forest model with atmospheric visibility observations and other auxiliary data. The modeled daily PM2.5 concentrations are in excellent agreement with ground measurements during 2015–2019, with a coefficient of determination of 0.95 and mean relative error of 12%. Besides the atmospheric visibility which explains 30% of total importance of variables in the model, both emissions and meteorological conditions are also essential factors affecting PM2.5 predictions. From 1980 to 2014, the model-predicted PM2.5 concentrations increased constantly with the maximum growth rate of 5–10 µg/m3/decade over eastern China. Due to the clean air actions, PM 2.5 concentrations have decreased effectively at a rate over 50 µg/m 3 /decade in the North China Plain and 20–50 µg/m 3 /decade over eastern China during 2014–2019. The newly generated dataset of 1-degree gridded PM 2.5 concentrations for the past 40 years across China provides a useful means for investigating interannual and decadal environmental and climate impacts related to aerosols.

54 ENVIRONMENTAL SCIENCES↗