Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computer science, Physics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Jarzynski Equality for Conditional Stochastic Work

It has been established that the inclusive work for classical, Hamiltonian dynamics is equivalent to the two-time energy measurement paradigm in isolated quantum systems. However, a plethora of other notions of quantum work has emerged, and thus the natural question arises whether any other quantum notion can provide motivation for purely classical considerations. In the present analysis, we propose the conditional stochastic work for classical, Hamiltonian dynamics, which is inspired by the one-time measurement approach. This novel notion is built upon the change of expectation value of the energy conditioned on the initial energy surface. As main results, we obtain a generalized Jarzynski equality and a sharper maximum work theorem, which account for how non-adiabatic the process is. Our findings are illustrated with the parametric harmonic oscillator.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

State Space Truncation with Quantified Errors for Accurate Solutions to Discrete Chemical Master Equation

The discrete chemical master equation (dCME) provides a general framework for studying stochasticity in mesoscopic reaction networks. Since its direct solution rapidly becomes intractable due to the increasing size of the state space, truncation of the state space is necessary for solving most dCMEs. It is therefore important to assess the consequences of state space truncations so errors can be quantified and minimized. Here we describe a novel method for state space truncation. By partitioning a reaction network into multiple molecular equivalence groups (MEGs), we truncate the state space by limiting the total molecular copy numbers in each MEG. We further describe a theoretical framework for analysis of the truncation error in the steady-state probability landscape using reflecting boundaries. By aggregating the state space based on the usage of a MEG and constructing an aggregated Markov process, we show that the truncation error of a MEG can be asymptotically bounded by the probability of states on the reflecting boundary of the MEG. Furthermore, truncating states of an arbitrary MEG will not undermine the estimated error of truncating any other MEGs. We then provide an overall error estimate for networks with multiple MEGs. To rapidly determine the appropriate size of an arbitrary MEG, we also introduce an a priori method to estimate the upper bound of its truncation error. This a priori estimate can be rapidly computed from reaction rates of the network, without the need of costly trial solutions of the dCME. As examples, we show results of applying our methods to the four stochastic networks of (1) the birth and death model, (2) the single gene expression model, (3) the genetic toggle switch model, and (4) the phage lambda bistable epigenetic switch model. We demonstrate how truncation errors and steady-state probability landscapes can be computed using different sizes of the MEG(s) and how the results validate our theories. Overall, the novel state space truncation and error analysis methods developed here can be used to ensure accurate direct solutions to the dCME for a large number of stochastic networks.

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

Computational study of full-scale VHTR lower plenum for turbulent mixing assessment

Next-generation nuclear reactors are poised to efficiently provide reliable power at enhanced safety levels for many years to come. Among the fundamental designs for these reactors is the very-hightemperature reactor (VHTR), which employs helium as the primary coolant and has the potential of reaching elevated temperatures such that chemical processing (e.g., electrolysis for hydrogen production) is viable alongside traditional electricity generation. This and other next-generation reactors also represent a strategic bridge for ultimately transitioning from fossil fuel-dominated energy portfolios to one where renewable options abound. For the VHTR, the ultimate goal of constructing a new plant will be preceded by additional fundamental research aimed at generating trusted models for behavior prediction under normal and accident scenarios. This paper presents a detailed computational fluid dynamics simulation of the lower plenum, where hot coolant from the core mixes together in a turbulent fashion before traveling to power conversion equipment. Because the flow is expected to enter the lower plenum across a wide range of temperatures and velocities, concerns exist when the mixing is incomplete. The potential hot spots on cylindrical support pedestals and the uniformity of the main outlet are the primary metrics of interest in this study, along with flow velocities and temperatures at different locations within the lower plenum. Three locations are probed in detail and suggest a large range exists for approach velocities and temperatures. A large degree of stratification is also seen on the surfaces of the support pedestals, suggesting a facility should accommodate testing for such behavior. The characterizations presented provide the valuable data needed in order to design appropriate experimental testbeds, where scaled modeling can be carried out in a manner meaningful in predicting the full-scale behavior.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Communication-efficient certificate revocation management for Advanced Metering Infrastructure and IoT Integration

Advanced Metering Infrastructure forms a communication network for the collection of power data from smart meters in Smart Grid. As the communication between smart meters could be secured utilizing public-key cryptography, however, public-key cryptography still has certain challenges in terms of certificate revocation and management particularly related distribution and storage overhead of revoked certificates. To address this challenge, in this work, we propose a novel revocation management approach by utilizing cryptographic accumulators which reduces the space requirements for revocation information significantly and thus enables efficient distribution of such information to all smart meters. We implemented the proposed approach on both ns-3 network simulator and a testbed. We demonstrated its superior performance with respect to traditional methods for revocation management.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Machine learning moment closure models for the radiative transfer equation I: Directly learning a gradient based closure

In this paper, we take a data-driven approach and apply machine learning to the moment closure problem for the radiative transfer equation in slab geometry. Instead of learning the unclosed high order moment, we propose to directly learn the gradient of the high order moment using neural networks. This new approach is consistent with the exact closure we derive for the free streaming limit and also provides a natural output normalization. Finally, a variety of benchmark tests, including the variable scattering problem, the Gaussian source problem with both periodic and reflecting boundaries, and the two-material problem, show both good accuracy and generalizability of our machine learning closure model.

Computer science↗
DOE OSTI2013

Structural studies of TiC{sub 1−x}O{sub x} solid solution by Rietveld refinement and first-principles calculations

The lattice parameters, structural stability and electronic structure of titanium oxycarbides (TiC{sub 1−x}O{sub x}, 0≤x≤1) solid solution were investigated by Rietveld refinement and first-principles calculations. Series of TiC{sub 1−x}O{sub x} were precisely synthesized by sintering process under the vacuum. Rietveld refinement results of XRD patterns show the properties of continuous solid solution in TiC{sub 1−x}O{sub x} over the whole composition range. The lattice parameters vary from 0.4324 nm to 0.4194 nm decreasing with increasing oxygen concentration. Results of first-principles calculations reveal that the disorder C/O structure is stable than the order C/O structure. Further investigations of the vacancy in Ti{sub 1−Va}(C{sub 1−x}O{sub x}){sub 1−Va} solid solution present that the structure of vacancy segregated in TiO-part is more stable than the disorder C/O structure, which can be ascribed to the Ti–Ti bond across O-vacancy and the charge redistributed around Ti-vacancy via the analysis of the electron density difference plots and PDOS. - Graphical abstract: XRD of series of titanium oxycarbides (TiC{sub 1−x}O{sub x}, 0≤x≤1) solid solution prepared by adjusting the proportion of TiO in the starting material. Highlights: • Titanium oxycarbides were obtained by sintering TiO and TiC under carefully controlled conditions. • Rietveld refinement results show continuous solid solution with FCC structure in TiC{sub 1−x}O{sub x}. • The disorder C/O structure is stable than the order C/O structure. • Introduction of vacancy segregated in TiO-part is more stable than disorder C/O structure. • Ti–Ti bond across O-vacancy and the charge redistributed around Ti-vacancy enhance structural stability.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Luminescent microporous metal–organic framework with functional Lewis basic sites on the pore surface: Quantifiable evaluation of luminescent sensing mechanisms towards Fe{sup 3+}

A systematic study has been conducted on a novel luminescent metal-organic framework, ([Zn(bpyp)(L-OH)]·DMF·2H{sub 2}O){sub n} (1), to explore its sensing mechanisms to Fe{sup 3+}. Structure analyses show that compound 1 exist pyridine N atoms and -OH groups on the pore surface for specific sensing of metal ions via Lewis acid-base interactions. On this consideration, the quenching mechanisms are studied and the processes are controlled by multiple mechanisms in which dynamic and static mechanisms are calculated, achieving the quantification evaluation of the quenching process. This work not only achieves the quantitative evaluation of the luminescence quenching but also provides certain insights into the quenching process, and the possible mechanisms explored in this work may inspire future research and design of target luminescent metal-organic frameworks (LMOFs) with specific functions. - Graphical abstract: A systematic study has been conducted on a novel luminescent metal-organic framework to explore its sensing mechanisms to Fe{sup 3+}. The quenching mechanisms are studied and the processes are controlled by multiple mechanisms in which dynamic and static mechanisms are calculated, achieving the quantification evaluation of the quenching process. - Highlights: • A novel porous luminescent MOF containing uncoordinated groups in interlayer channels was successfully synthesized. • The compound 1 can exhibit significant luminescent sensitivity to Fe{sup 3+}, which make its good candidate as luminescent sensor. • The corresponding dynamic and static quenching constants are calculated, achieving the quantification evaluation of the quenching process.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Electronic, elastic, and optical properties of monolayer BC{sub 2}N

The structural stability, electronic structure, elasticity, and optical properties of four types of monolayer BC{sub 2}N have been investigated from first principles using calculation based on density functional theory. The results show that the structural stability of BC{sub 2}N increases with the number of C–C and B–N bonds. By calculating the two-dimensional Young's modulus, shear modulus, Poisson's ratio, and shear anisotropic factors in different directions, four structures present various anisotropies and the most stable structure is almost isotropic. For C-type BC{sub 2}N, the values of two-dimensional Young's modulus, shear modulus, and bulk modulus (309, 128, 195 GPa m{sup −1}), are smaller than those of graphene (343, 151, 208) but bigger than those of h-BN (286, 185, 116). Furthermore, the dielectric function, refractive index, reflectivity, absorption coefficient, and energy loss spectrum are also calculated to investigate the mechanism underpinning the optical transitions in BC{sub 2}N, revealing monolayer BC{sub 2}N as a candidate window material. - Graphical abstract: Schematic diagram of BC{sub 2}N under the biaxial tensile strain. Changes in the valence-band top and the conduction-band bottom of BC{sub 2}N with increasing strain.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Invariant Molecular Representations for Heterogeneous Catalysis

Catalyst screening is a critical step in the discovery and development of heterogeneous catalysts, which are vital for a wide range of chemical processes. In recent years, computational catalyst screening, primarily through density functional theory (DFT), has gained significant attention as a method for identifying promising catalysts. However, the computation of adsorption energies for all likely chemical intermediates present in complex surface chemistries is computationally intensive and costly due to the expensive nature of these calculations and the intrinsic idiosyncrasies of the methods or data sets used. This study introduces a novel machine learning (ML) method to learn adsorption energies from multiple DFT functionals by using invariant molecular representations (IMRs). To do this, we first extract molecular fingerprints for the reaction intermediates and later use a Siamese-neural-network-based training strategy to learn invariant molecular representations or the IMR across all available functionals. Our Siamese network-based representations demonstrate superior performance in predicting adsorption energies compared with other molecular representations. Notably, when considering mean absolute values of adsorption energies as 0.43 eV (PBE-D3), 0.46 eV (BEEF-vdW), 0.81 eV (RPBE), and 0.37 eV (scan+rVV10), our IMR method has achieved the lowest mean absolute errors (MAEs) of 0.18 0.10, 0.16, and 0.18 eV, respectively. These results emphasize the superior predictive capacity of our Siamese network-based representations. The empirical findings in this study illuminate the efficacy, robustness, and dependability of our proposed ML paradigm in predicting adsorption energies, specifically for propane dehydrogenation on a platinum catalyst surface.

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

Atomistic Characterization of Gramicidin Channel Formation

We investigated gramicidin A (gA) subunit dimerization in lipid bilayers using microsecond-long replica-exchange umbrella sampling simulations, millisecond-long unbiased molecular dynamics simulations, and machine learning. Our simulations led to a dimer structure that is indistinguishable from the experimentally determined gA channel structures, with the two gA subunits joined by six hydrogen bonds (6HB). The simulations also uncovered two additional dimer structures, with different gA–gA stacking orientations that were stabilized by four or two hydrogen bonds (4HB or 2HB). When examining the temporal evolution of the dimerization, we found that two bilayer-inserted gA subunits can form the 6HB dimer directly, with no discernible intermediate states, as well as through paths that involve the 2HB and 4HB dimers.

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

Generating Protein Structures for Pathway Discovery Using Deep Learning

Resolving the intricate details of biological phenomena at the molecular level is fundamentally limited by both length- and time scales that can be probed experimentally. Molecular dynamics (MD) simulations at various scales are powerful tools frequently employed to offer valuable biological insights beyond experimental resolution. However, while it is relatively simple to observe long-lived, stable configurations of, for example, proteins, at the required spatial resolution, simulating the more interesting rare transitions between such states often takes orders of magnitude longer than what is feasible even on the largest supercomputers available today. One common aspect of this challenge is pathway discovery, where the start and end states of a scientific phenomenon are known or can be approximated, but the mechanistic details in between are unknown. Here, we propose a representation-learning-based solution that uses interpolation and extrapolation in an abstract representation space to synthesize potential transition states, which are automatically validated using MD simulations. The new simulations of the synthesized transition states are subsequently incorporated into the representation learning, leading to an iterative framework for targeted path sampling. Our approach is demonstrated by recovering the transition of a RAS-RAF protein domain (CRD) from membrane-free to interacting with the membrane using coarse-grain MD simulations.

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

Direct Ab Initio Simulation of the Synthesis of BaZrO 3 and the Microstructure Impacts on Proton Transport

Controlling and predicting the processing-structure-performance relationship in functional materials is a grand challenge in materials science, with important implications for a wide range of emerging applications; a high fidelity understanding of the performance impact of microstructures formed under synthesis conditions is required to develop advanced materials, such as solid-state fuel cells and electrolyzers. Using the ceramic BaZrO 3 as a case study, we directly simulate the synthesis and investigate how proton transport is dictated by microstructures. We develop a framework that couples density functional theory (DFT), machine-learning interatomic potential (MLIP) driven molecular dynamics, and grand canonical Monte Carlo to perform large-scale, microstructure-resolved, atomistic simulations of proton transport in experimentally representative polycrystalline structures. Our fully ab initio approach, using a MLIP as a proxy for DFT, allows us to quantify the competition between two distinct diffusion mechanisms: one associated with grain-boundary regions and another within grains. When the impacts of grain boundaries are taken into account, proton transport exhibits substantial deviation from the bulk oxide limit. This addresses long-standing discrepancies between theory and experiments. Our integrated approach provides atomistic insight into microstructure-dependent proton pathways in BaZrO 3 and establishes a general protocol for predicting processing-structure-performance relationships.

organic↗
DOE OSTI2014

Hydroxide Degradation Pathways for Imidazolium Cations. A DFT Study

Imidazolium cations are promising candidates as covalently tetherable cations for application in anion exchange membranes. They have generated specific interest in alkaline membrane fuel cell applications where ammonium-based cations have been the most commonly applied but have been found to be susceptible to hydroxide attack. In the search for high stability cations, a detailed understanding of the degradation pathways and reaction barriers is required. In this work, we investigate imidazolium and benzimidazolium cations in the presence of hydroxide using density functional theory calculations for their potential in alkaline membrane fuel cells. Moreover, the dominant degradation pathway for these cations is predicted to be the nucleophilic addition–elimination pathway at the C-2 atom position on the imidazolium ring. Steric interferences, introduced by substitutions at the C-2, C-4, and C-5 atom positions, were investigated and found to have a significant, positive impact on calculated degradation energy barriers. Benzimidazolium cations, with their larger conjugated systems, are predicted to degrade much faster than their imidazolium counterparts. Our results provide important insight into designing stable cations for anion exchange membranes. Some of the molecules studied have significantly increased degradation energy barriers suggesting that they could possess significantly improved (several orders of magnitude) durability compared to traditional cations and potentially enable new applications.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2017

Simple spatial scaling rules behind complex cities

Although most of wealth and innovation have been the result of human interaction and cooperation, we are not yet able to quantitatively predict the spatial distributions of three main elements of cities: population, roads, and socioeconomic interactions. By a simple model mainly based on spatial attraction and matching growth mechanisms, we reveal that the spatial scaling rules of these three elements are in a consistent framework, which allows us to use any single observation to infer the others. All numerical and theoretical results are consistent with empirical data from ten representative cities. In addition, our model can also provide a general explanation of the origins of the universal super- and sub-linear aggregate scaling laws and accurately predict kilometre-level socioeconomic activity. Our work opens a new avenue for uncovering the evolution of cities in terms of the interplay among urban elements, and it has a broad range of applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Visualizing group II intron dynamics between the first and second steps of splicing

Group II introns are ubiquitous self-splicing ribozymes and retrotransposable elements evolutionarily and chemically related to the eukaryotic spliceosome, with potential applications as gene-editing tools. Recent biochemical and structural data have captured the intron in multiple conformations at different stages of catalysis. Here, we employ enzymatic assays, X-ray crystallography, and molecular simulations to resolve the spatiotemporal location and function of conformational changes occurring between the first and the second step of splicing. We show that the first residue of the highly-conserved catalytic triad is protonated upon 5’-splice-site scission, promoting a reversible structural rearrangement of the active site (toggling). Protonation and active site dynamics induced by the first step of splicing facilitate the progression to the second step. Our insights into the mechanism of group II intron splicing parallels functional data on the spliceosome, thus reinforcing the notion that these evolutionarily-related molecular machines share the same enzymatic strategy.

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

Automated exploitation of the big configuration space of large adsorbates on transition metals reveals chemistry feasibility

Mechanistic understanding of large molecule conversion and the discovery of suitable heterogeneous catalysts have been lagging due to the combinatorial inventory of intermediates and the inability of humans to enumerate all structures. Here, we introduce an automated framework to predict stable configurations on transition metal surfaces and demonstrate its validity for adsorbates with up to 6 carbon and oxygen atoms on 11 metals, enabling the exploration of ~10 8 potential configurations. It combines a graph enumeration platform, force field, multi-fidelity DFT calculations, and first-principles trained machine learning. Clusters in the data reveal groups of catalysts stabilizing different structures and expose selective catalysts for showcase transformations, such as the ethylene epoxidation on Ag and Cu and the lack of C-C scission chemistry on Au. Deviations from the commonly assumed atom valency rule of small adsorbates are also manifested. This library can be leveraged to identify catalysts for converting large molecules computationally.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗