Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “REPLICAS”

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

Nonparametric, data-based kernel interpolation for particle-tracking simulations and kernel density estimation

Traditional interpolation techniques for particle tracking include binning and convolutional formulas that use pre-determined (i.e., closed-form, parameteric) kernels. In many instances, the particles are introduced as point sources in time and space, so the cloud of particles (either in space or time) is a discrete representation of the Green’s function of an underlying PDE. As such, each particle is a sample from the Green’s function; therefore, each particle should be distributed according to the Green’s function. In short, the kernel of a convolutional interpolation of the particle sample “cloud” should be a replica of the cloud itself. This idea gives rise to an iterative method by which the form of the kernel may be discerned in the process of interpolating the Green’s function. When the Green’s function is a density, this method is broadly applicable to interpolating a kernel density estimate based on random data drawn from a single distribution. We formulate and construct the algorithm and demonstrate its ability to perform kernel density estimation of skewed and/or heavy-tailed data including breakthrough curves.

42 ENGINEERING↗
DOE OSTI2015

Excited-state quantum phase transitions in systems with two degrees of freedom: II. Finite-size effects

This article extends our previous analysis Stránský et al. (2014) of Excited-State Quantum Phase Transitions (ESQPTs) in systems of dimension two. We focus on the oscillatory component of the quantum state density in connection with ESQPT structures accompanying a first-order ground-state transition. It is shown that a separable (integrable) system can develop rather strong finite-size precursors of ESQPT expressed as singularities in the oscillatory component of the state density. The singularities originate in effectively 1-dimensional dynamics and in some cases appear in multiple replicas with increasing excitation energy. Using a specific model example, we demonstrate that these precursors are rather resistant to proliferation of chaotic dynamics. - Highlights: • Oscillatory components of state density and spectral flow studied near ESQPTs. • Enhanced finite-size precursors of ESQPT caused by fully/partly separable dynamics. • These precursors appear due to criticality of a subsystem with lower dimension. • Separability-induced finite-size effects disappear in case of fully chaotic dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Non-ergodic delocalized phase in Anderson model on Bethe lattice and regular graph

We develop a novel analytical approach to the problem of single particle localization in infinite dimensional spaces such as Bethe lattice and random regular graph models. The key ingredient of the approach is the notion of the inverted order thermodynamic limit (IOTL) in which the coupling to the environment goes to zero before the system size goes to infinity. Using IOTL and Replica Symmetry Breaking (RSB) formalism we derive analytical expressions for the fractal dimension D{sub 1} that distinguishes between the extended ergodic, D{sub 1}=1, and extended non-ergodic (multifractal), 0<D{sub 1}<1 states on the Bethe lattice and random regular graphs with the branching number K. We also employ RSB formalism to derive the analytical expression lnS{sub typ}{sup −1}=−〈lnS〉∼(W{sub c}−W){sup −1} for the typical imaginary part of self-energy S{sub typ} in the non-ergodic phase close to the Anderson transition in the conventional thermodynamic limit. We prove the existence of an extended non-ergodic phase in a broad range of disorder strength and energy and establish the phase diagrams of the models as a function of disorder and energy. The results of the analytical theory are compared with large-scale population dynamics and with the exact diagonalization of Anderson model on random regular graphs. We discuss the consequences of these results for the many body localization.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Effect of reaction media on hydrogenolysis of polyethylene plastic waste: Polymer-surface interactions in small alkane/polymer blends

The polymer reaction media and its properties can be altered by recycling a fraction of liquid products or adding alkane solvents. Less clear is whether this strategy affects hydrogenolysis. Herein, we investigated the effect of short-chain alkanes C n consisting of n carbons (n=8, 16, and 32) on the upcycling of high-density polyethylene (HDPE) plastic waste to lubricant-range products over Ru/TiO 2 catalysts by multiscale simulations and experiments. First, we trained a force field for polymer/surface interactions on a Ru 22 nanoparticle (NP) supported on TiO 2 . Using replica exchange molecular dynamics simulations, we studied the effect of small hydrocarbons on the adsorption of a surrogate polymer, C 142 , on the catalyst. We found segregation of long chains (C 142 ) at the catalyst surface due to the enthalpy gained by adsorbing more C-C bonds of the long chains, compensating for entropic losses upon adsorption. Short-chain molecules decrease the adsorbed carbons of long chains on the Ru NP due to blocking Ru active sites. Compared to the bulk chains, competitive adsorption results in a broader, heavy-tailed distribution of end-to-end distance of adsorbed chains. Our experiments demonstrated that catalyst activity declines significantly beyond simple dilution due to changes in polymer adsorption, and tuning the reaction media by creating suitable blends impacts hydrogenolysis. Density distributions for a 50:50%wt mixture of PP and PE show that PE chains are segregated at the surface, so they are prone to C-C bond breaking much faster than PP chains. H/D exchange experiments show preferential deuteration of PE, while CH 3 groups of PP remain undeuterated. Furthermore, this may be explained by the preferential sorption of PE over PP, leading to specific distribution in the polymer blend.

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

A digital twin for rapid qualification of 3D printed metallic components

The customized production of complex components by 3D printing has been hailed as a potentially transformative tool in manufacturing with important applications in health care, automotive and aerospace industries. However, after about a quarter of a century of research and development, only a handful of commercial alloys can be printed and the market value of all 3D printed products now amounts to a negligible portion of the manufacturing economy. This difficulty is attributable to a remarkable diversity in structure and properties of the printed components and susceptibility to defects. In addition, the current practice of qualifying components by prolonged trial and error with expensive printing equipment and feed stock material confine the printed products to a niche market where the high product cost and the delay in the qualification are not critical factors. In this report we explain how a digital twin or a digital replica of the printing machine will reduce the number of trial and error tests to obtain desired product attributes and reduce the time required for part qualification to make the printed components cost effective. It is shown that a comprehensive digital twin of 3D printing machine consisting of mechanistic, control and statistical models of 3D printing, machine learning and big data can reduce the volume of trial and error testing, reduce defects and shorten time between the design and production.

36 MATERIALS SCIENCE↗
DOE OSTI2014

The role of histone tails in the nucleosome: A computational study

Histone tails play an important role in gene transcription and expression. We present here a systematic computational study of the role of histone tails in the nucleosome, using replica exchange molecular dynamics simulations with an implicit solvent model and different well-established force fields. We performed simulations for all four histone tails, H4, H3, H2A, and H2B, isolated and with inclusion of the nucleosome. The results confirm predictions of previous theoretical studies for the secondary structure of the isolated tails but show a strong dependence on the force field used. In the presence of the entire nucleosome for all force fields, the secondary structure of the histone tails is destabilized. Specific contacts are found between charged lysine and arginine residues and DNA phosphate groups and other binding sites in the minor and major DNA grooves. Using cluster analysis, we found a single dominant configuration of binding to DNA for the H4 and H2A histone tails, whereas H3 and H2B show multiple binding configurations with an equal probability. The leading stabilizing contribution for those binding configurations is the attractive interaction between the positively charged lysine and arginine residues and the negatively charged phosphate groups, and thus the resulting charge neutralization. Lastly, we present results of molecular dynamics simulations in explicit solvent to confirm our conclusions. Results from both implicit and explicit solvent models show that large portions of the histone tails are not bound to DNA, supporting the complex role of these tails in gene transcription and expression and making them possible candidates for binding sites of transcription factors, enzymes, and other proteins.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Insight into molecular basis and dynamics of full-length CRaf kinase in cellular signaling mechanisms

Raf kinases play key roles in signal transduction in cells for regulating proliferation, differentiation, and survival. Despite decades of research into functions and dynamics of Raf kinases with respect to other cytosolic proteins, understanding Raf kinases is limited by the lack of their full-length structures at the atomic resolution. Here, we present the first model of the full-length CRaf kinase obtained from artificial intelligence/machine learning algorithms with a converging ensemble of structures simulated by large-scale temperature replica exchange simulations. Our model is validated by comparing simulated structures with the latest cryo-EM structure detailing close contacts among three key domains and regions of the CRaf. Our simulations identify potentially new epitopes of intramolecule interactions within the CRaf and reveal a dynamical nature of CRaf kinases, in which the three domains can move back and forth relative to each other for regulatory dynamics. The dynamic conformations are then used in a docking algorithm to shed insight into the paradoxical effect caused by vemurafenib in comparison with a paradox breaker PLX7904. In this study, we propose a model of Raf-heterodimer/KRas-dimer as a signalosome based on the dynamics of the full-length CRaf.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Monte Carlo simulation of order-disorder transition in refractory high entropy alloys: A data-driven approach

High entropy alloys (HEAs) are a series of novel materials that demonstrate many exceptional mechanical properties. To understand the origin of these attractive properties, it is important to investigate the thermodynamics and elucidate the evolution of various chemical phases. In this work, we introduce a machine learning approach to construct the effective Hamiltonian and study the thermodynamics of HEAs through canonical Monte Carlo simulation. The main characteristic of our method is to use pairwise interactions between atoms as the machine learning features and systematically improve the representativeness of the dataset using active learning from the Monte Carlo samples. We find this method produces highly robust and accurate effective Hamiltonians that give less than 0.1 mRy test error for all the three refractory HEAs: MoNbTaW, MoNbTaVW, and MoNbTaTiW. Using replica exchange to speed up the MC simulation, we calculated the specific heats and short-range order parameters in a wide range of temperatures. For all the studied materials, we find there are two major order-disorder transitions occurring respectively at T 1 and T 2 , where T 1 is near room temperature but T 2 is much higher. Further, we demonstrate that the transition at T 1 is caused by W and Nb while the one at T 2 is caused by the other elements. By comparing with experiments, the results demonstrate that the order-disorder transitions have a profound impact on the strength and ductility of HEAs.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Flow visualisation in real-size optical injectors of conventional, additised, and renewable gasoline blends

Research on renewable and alternative fuels is crucial for improving the energy and environmental efficiency of modern gasoline internal combustion engines. To highlight the influence of fuel rheological and thermodynamic properties on phase change and atomisation processes, three types of gasoline blends were tested. More specifically, the campaign comprised a reference gasoline, an ethanol/gasoline blend (10% v/v) representative of renewable fuels, and an additised gasoline sample treated with viscoelasticity-inducing agents. High-speed imaging of the transient two-phase flow field arising in the internal geometry and the near-nozzle spray region of gasoline injectors was performed employing Diffuse Backlight Illumination. The metallic body of a commercial injector was modified to fit transparent tips realising two nozzle layouts, namely a two-hole real size model resembling the Engine Combustion Network spray G injector and an enraged replica with an offset hole. Experiments were conducted at realistic operating conditions comprising an injection pressure of 100 bar and ambient pressures in the range of 0.1–6.0 bar to cover the entire range of chamber pressures prevailing in Gasoline Direct Injection engines. The action of viscoelastic additives was verified to have a suppressive effect on in-nozzle cavitation (6% reduction in cavitation extent) , while also enhancing spray atomisation at flash-boing conditions, in a manner resembling the more volatile gasoline/ethanol blends. Finally, persisting liquid ligaments were found to form after the end of injection for the additised sample, owing to the surfactant nature of the additives.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2022

Reticulated porous lanthanum strontium manganite structures for solar thermochemical hydrogen production

Porous ceramic foams are utilized in various fields and applications due to their advantageous properties such as high porosity, high specific surface area, and controlled permeability. Non-stoichiometric oxide foams are of particular interest for the solar thermochemical production of hydrogen in directly irradiated solar receivers due to their enhanced heat transfer capabilities and structural integrity. Here we report the fabrication of lanthanum strontium manganite replica foams and their solar thermochemical water splitting performance. Highly porous (>80%) foams were produced with controlled structure and tested in a high temperature water splitting furnace to demonstrate stability and overall performance. The temperature profile of the burnout and sintering of the foam structure were altered based on simultaneous differential scanning calorimetry and thermal gravimetry analysis. Two-step sintering was used both to retain micron-sized pores (D50 = 0.4 μm) on the foam's struts and to limit grain growth to a grain size of ~1 μm. The structural characteristics were confirmed by porosimetry and scanning electron microscopy. The foams were repeatedly cycled, producing approximately 200 μmol/g of H 2 per cycle. Water splitting performance was tested for 50 cycles and approximately 35 h at 1400 °C. No statistical difference in the effective hydrogen production per cycle was observed, which is promising for long-term redox cycle stability.

08 HYDROGEN↗
DOE OSTI2016

Modeling of inertial deposition in scaled models of rat and human nasal airways: Towards in vitro regional dosimetry in small animals

Rodents are routinely used in inhalation toxicology tests as human surrogates. However, in vitro dosimetry tests in rodent casts are still scarce due to small rodent airways and in vitro tests to quantify sub-regional dosimetry are still impractical. We hypothesized that for inertial particles whose deposition is dominated by airflow convection (Reynolds number) and particle inertia (Stokes number), the deposition should be similar among airway replicas of different scales if their Reynolds and Stokes numbers are kept the same. In this study, we aimed to (1) numerically test the hypothesis in three airway geometries: a USP induction port, a human nose model, and a Sprague-Dawley rat nose model, and (2) find the range of applicability of this hypothesis. Five variants of the USP and human nose models and three variants of the rat nose model were tested. Inhalation rates and particle sizes were scaled to match the Reynolds number and Stokes numbers. A low-Reynolds-number k–ω model was used to resolve the airflow and a Lagrangian tracking algorithm was used to simulate the particle transport and deposition. Statistical analysis of predicted doses was conducted using ANOVA. For normal inhalation rates and particle dia- meters ranging from 0.5 to 24 mm, the deposition differences between the life-size and scaled models are insignificant for all airway geometries considered (i.e., human nose, USP, and rat nose). Furthermore, the deposition patterns and exit particle profiles also look similar among scaled models. However, deposition rates and patterns start to deviate if inhalation rates are too low, or particle sizes are too large. For the rat nose, the threshold velocity was found to be 0.71 m/s and the threshold Froude number to be 50. Results of this study provide a theoretical foundation for sub-regional in vitro dosimetry tests in small animals and for interpretation of data from inter-species or intra-species with varying body sizes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Molecular simulation of lignin-related aromatic compound permeation through gram-negative bacterial outer membranes

Lignin, an abundant aromatic heteropolymer in secondary plant cell walls, is the single largest source of renewable aromatics in the biosphere. Leveraging this resource for renewable bioproducts through targeted microbial action depends on lignin fragment uptake by microbial hosts and subsequent enzymatic action to obtain the desired product. Recent computational work has emphasized that bacterial inner membranes are permeable to many aromatic compounds expected from lignin depolymerization processes. In this study, we expand on these findings through simulations for 42 lignin-related compounds across a gram-negative bacterial outer membrane model. Unbiased simulation trajectories indicate that spontaneous crossing for the full outer membrane is relatively rare at molecular simulation timescales, primarily due to preferential membrane partitioning and slow diffusion within the lipopolysaccharide layer within the outer membrane. Membrane partitioning and permeability coefficients were determined through replica exchange umbrella sampling simulations to overcome sampling limitations. We find that the glycosylated lipopolysaccharides found in the outer membrane increase the permeation barrier to many lignin-related compounds, particularly the most hydrophobic compounds. However, the effect is relatively modest; at industrially relevant concentrations, uncharged lignin-related compounds will readily diffuse across the outer membrane without the need for specific porins. Together, our results provide insight into the permeability of the bacterial outer membrane for assessing lignin fragment uptake and the future production of renewable bioproducts.

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

Uncertainty quantification methodology for hyperbolic systems with application to blood flow in arteries

We present a Stochastic Finite Volume - ADER (SFV-ADER) methodology for Uncertainty Quantification (UQ) in the general framework of systems of hyperbolic balance laws with uncertainty in parameters. The resulting method is weakly intrusive, meaning that deterministic solvers need minor modifications to include random contributions, and has no theoretical accuracy barrier. An illustration of the second-order version for the viscous Burgers equation with uncertain viscosity coefficient is first given in detail, under different deterministic initial conditions; the attainment of the theoretically expected convergence rate is demonstrated empirically; results and features of our scheme are discussed. We then extend the SFV-ADER method to a non-linear hyperbolic system with source terms that models one-dimensional blood flow in arteries, assuming uncertainties in a set of parameters of the problem. Results are then compared with measurements in a well-defined 1:1 experimental replica of the network of largest arteries in the human systemic circulation. For a 99%-confidence level, severe influence of parameter fluctuations is seen in pressure profiles, while secondary effects on flow rate profiles become visible, particularly in terminal branches. In conclusion, results suggest that the proposed methodology could successfully be applied to other physical problems involving hyperbolic balance laws.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Discrete dislocation dynamics for crystal RVEs. Part 1: Periodic network kinematics

A novel implementation of the dislocation flux boundary condition in discrete dislocation dynamics is presented. The continuity of the individual dislocation loops in a periodic representative crystal volume (RVE) is enforced across the boundary of the RVE with the help of a dual topological description for representing dislocation line kinematics in two equivalent spaces representing the deforming crystal, the RVE and the unbounded crystal spaces. The former describes the motion of the dislocations in the simulated crystal RVE whereas the latter represents the motion of dislocations in an infinite space containing all replicas of the RVE. A mapping between the two spaces forms the basis of the implementation of flux boundary condition. The implementation details are discussed in the context of statistical homogeneity of bulk crystals undergoing macroscopically homogeneous plastic deformation. In this case, the boundary nodes associated with dislocation segments bear no relevance in the motion of the dislocations. Finally, some test cases are presented and discussed to establish the proposed approach.

42 ENGINEERING↗
DOE OSTI2026

Physics-aware adaptive checkpointing with shadow systems for nonlinear PDE simulations

Large-scale simulations of nonlinear partial differential equations (PDEs) that exhibit strongly transient behavior and pattern-forming dynamics produce enormous amounts of data, which, even with modern storage systems, cannot be stored for later curation. Current I/O strategies either write dense time series of snapshots, which is often prohibitive in I/O and storage, or store a few checkpoints that enable restart but incur expensive recomputation cost and provide no control over post-restart error growth, especially when lossy compression is used. Moreover, most, if not all, existing strategies take no account of the actual physical state of the system. Here, we present a simple physics-aware I/O framework in which a low-cost shadow system adaptively triggers lossy checkpoints when the shadow system deviates from the fine-scale simulation. The shadow system can be a coarsened replica of the fine-scale simulation that evolves concurrently. This means that checkpoints are taken based on the physical state of the system: fewer checkpoints are triggered when the system is quiescent while more are taken when the system undergoes a rapid change. This type of behavior is observed in many systems such as Brusselator and FitzHugh–Nagumo. We illustrate that our framework maintains stable restarts, keeps fine-scale restart errors bounded by shadow errors, and reconstructs the time history with significantly lower error and storage than interpolating fixed-interval snapshots, with low-cost shadow replay and modest online synchronization overhead.

Gong, Qian [ORNL] (ORCID:0000000235704142)↗
DOE OSTI2025

Entanglement entropy of a color flux tube in (1+1)D Yang–Mills theory

In recent work Amorosso et al. (2024), we computed a novel flux tube entanglement entropy (FTE) of the color flux tube stretched between a heavy quark-antiquark pair on a Euclidean lattice in (2+1)D Yang–Mills theory. Our numerical results suggested that FTE can be partitioned into an internal color entanglement entropy and a vibrational entropy corresponding to the transverse excitations of a QCD string, with the latter described by a thin string model. Since the color flux tube does not have transverse excitations in (1+1)D Yang–Mills theory, we use this simpler framework to perform an exact analytical computation of the contribution of the internal color degrees of freedom to FTE. For the multipartite partitioning of the color flux tube, we find the remarkable result that FTE depends only on the dimension of the color group representation and the number of times the flux tube crosses the boundary between the traced and untraced spatial regions but not on the string length. Our proof is independent of whether the replica and region boundaries on the lattice are placed on vertices or in plaquette centers.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Quantification of water content by laser induced breakdown spectroscopy on Mars

Laser induced breakdown spectroscopy (LIBS), as performed by the ChemCam instrument, provides a new technique to measure hydrogen at the surface of Mars. Using a laboratory replica of the LIBS instrument onboard the Curiosity rover, different types of hydrated samples (basalts, calcium and magnesium sulfates, opals and apatites) covering a range of targets observed on Mars have been characterized and analyzed in this paper. A number of factors related to laser parameters, atmospheric conditions and differences in targets properties can affect the standoff LIBS signal, and in particular the hydrogen emission peak. Dedicated laboratory tests were run to identify a normalization of the hydrogen signal which could best compensate for these effects and enable the application of the laboratory calibration to Mars data. We check that the hydrogen signal increases linearly with water content; and normalization of the hydrogen emission peak using to oxygen and carbon emission peaks (related to the breakdown of atmospheric carbon dioxide) constitutes a robust approach. Finally, moreover, the calibration curve obtained is relatively independent of the samples types.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗