Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complex systems theory”

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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Effect of lithium diffusion into Ga 2 O 3 thin films

The integration of lithium based compounds (e.g., Li:NiO/Ga 2 O 3 , LiGa 5 O 8 /Ga 2 O 3 ) in Ga 2 O 3 based pn-heterojunctions raises concerns about interface stability, i.e., Li diffusion effects on Ga 2 O 3 properties. In this work the ex-situ diffusion of Li is investigated in three different Ga 2 O 3 epilayers systems [(001) κ-Ga 2 O 3 and (−201) β-Ga2O3 heteroepitaxy on (001) α-Al 2 O 3 , and (010) β-Ga 2 O 3 homoepitaxy] at relevant temperatures for the synthesis / processing of Li-based epilayers. It is here experimentally demonstrated and quantified the Li diffusion in all the investigated Ga 2 O 3 epilayers systems and Li bulk (D Li,bulk ) and 2D defects (D Li,2D ) diffusion coefficients are provided. In the case of the (010) β-Ga 2 O 3 homoepitaxial layer (nominally free of structural defects), hybrid functional theory calculations foresee a diffusion mechanism mediated by Ga vacancies (VGa). Moreover, in the (010) β-Ga 2 O 3 homo-layer a significant effect on its functional properties (e.g., additional Raman vibrational modes, induced conductivity in an otherwise insulating sample) upon the Li-diffusion process is experimentally highlighted and tentatively related to the passivation of acceptor defects (i.e., formation of V Ga -nLi complexes).

Defects↗
NASA NTRS1975

Frequency stability requirements for two way range rate tracking

Accuracy limitations to two way range rate Doppler tracking due to master (reference) oscillator frequency instabilities are discussed. Theory is developed to treat both the effects of random and nonrandom oscillator instabilities. The nonrandom instabilities treated are drift, environmental effects, and coherent phase modulation. The effects of random instabilities on range rate accuracy are shown to be describable in terms of sigma y (2, T, tau). For the typical noise processes encountered in precision oscillators, range rate error is related to the more familiar sigma y (tau) and script L (f). Three examples are discussed to show how to determine range rate error from given sigma y (tau) or script L (f) curves, and approximations are developed to simplify the treatment of complex systems. An error analysis of range determined from rate data is also given.

Reinhardt, V.↗
NASA NTRS1987

Analysis of a parallelized nonlinear elliptic boundary value problem solver with application to reacting flows

A parallelized finite difference code based on the Newton method for systems of nonlinear elliptic boundary value problems in two dimensions is analyzed in terms of computational complexity and parallel efficiency. An approximate cost function depending on 15 dimensionless parameters is derived for algorithms based on stripwise and boxwise decompositions of the domain and a one-to-one assignment of the strip or box subdomains to processors. The sensitivity of the cost functions to the parameters is explored in regions of parameter space corresponding to model small-order systems with inexpensive function evaluations and also a coupled system of nineteen equations with very expensive function evaluations. The algorithm was implemented on the Intel Hypercube, and some experimental results for the model problems with stripwise decompositions are presented and compared with the theory. In the context of computational combustion problems, multiprocessors of either message-passing or shared-memory type may be employed with stripwise decompositions to realize speedup of O(n), where n is mesh resolution in one direction, for reasonable n.

Keyes, David E.↗
DOE OSTI2023

Theoretical Investigations on the Plasmon-Mediated Dissociation of Small Molecules in the Presence of Silver Atomic Wires

Plasmonic nanoparticles can promote bond activation in adsorbed molecules under relatively benign conditions via excitation of the nanoparticle’s plasmon resonance. As the plasmon resonance often falls within the visible light region, plasmonic nanomaterials are a promising class of catalysts. However, the exact mechanisms through which plasmonic nanoparticles activate the bonds of nearby molecules are still unclear. Herein, we evaluate Ag 8 –X 2 (X = N, H) model systems via real-time time-dependent density functional theory (RT-TDDFT), linear response time-dependent density functional theory (LR-TDDFT), and Ehrenfest dynamics in order to better understand the bond activation processes of N 2 and H 2 facilitated by the presence of the atomic silver wire under excitation at the plasmon resonance energies. We find that dissociation is possible for both small molecules at high electric field strength. Activation of each adsorbate is symmetry- and electric field-dependent, and H 2 activates at lower electric field strengths than N 2 . Furthermore, this work serves as a step toward understanding the complex time-dependent electron and electron–nuclear dynamics between plasmonic nanowires and adsorbed small molecules.

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

Real-Frequency Response Functions at Finite Temperature

Building on previous developments [A. Taheridehkordi, S. H. Curnoe, and J. P. F. LeBlanc, Phys. Rev. B 99, 035120 (2019); A. Taheridehkordi, S. H. Curnoe, and J. P. F. LeBlancPhys. Rev. B101, 125109 (2020); A. Taheridehkordi, S. H. Curnoe, and J. P. F. LeBlancPhys. Rev. B102, 045115 (2020), B. Holm and U. von Barth, Phys. Rev. B 57, 2108 (1998), J. Vičičevic and M. Ferrero, Phys. Rev. B 101, 075113 (2020)], we show that the diagrammatic Monte Carlo technique allows us to compute finite-temperature response functions directly on the real-frequency axis within any field-theoretical formulation of the interacting fermion problem. There are no limitations on the type and nature of the system’s action or whether partial summation and self-consistent treatment of certain diagram classes are used. In particular, by eliminating the need for numerical analytic continuation from a Matsubara representation, our scheme allows us to study spectral densities of arbitrary complexity with controlled accuracy in models with frequency-dependent effective interactions. Finally, for illustrative purposes we consider the problem of the plasmon linewidth in a homogeneous electron gas (jellium).

3-dimensional systems↗
DOE OSTI2025

Derivation of A Representative Elementary Volume (REV) for Upscaled Two-Phase Flow in Porous Media

Relative permeability plays an important role in the upscaling of multiphase flow in porous media from the pore scale to the Darcy scale. The entire concept of relative permeability is contingent on the existence of a representative elementary volume (REV). As we move to smaller samples to measure relative permeability, such as with digital core analysis, the concept of a classical REV has become increasingly unlikely when using the conventional approach to defining a representative volume. The “‘conventional”’ understanding of an REV is that a large enough volume must be considered such that spatial variability averages out. In digital rock methods, such as pore-scale simulations based on micro-computed tomography (CT) images, the domain size is typically 2 to 4 mm. This is approximately the length scale of a single-phase flow REV using the classic REV approach. However, the single-phase perspective does not consider the complex dynamics and fluctuations often observed in multiphase flow systems, even at centimeter-scale experiments and/or simulations. A fundamental question is, therefore, whether the domain size commonly used in digital rock simulations can provide a consistent energy budget such that the concept of relative permeability exists. Based on first principles, relative permeability accounts for the rate of energy dissipated in a stationary process. If the dynamics are fluctuating, the energy dissipated can vary but will average out over a long enough timescale. The key to determining the validity of the relative permeability is the timescale of the measurement, not the spatial scale. The conventional REV theory assumes that spatial, temporal, and ensemble averages are equivalent in an ergodic system, but it does not provide a way to test this assumption. We provide a formal way to identify the timescale where the relative permeability accurately captures energy dissipation as a way to validate relative permeability measurements and quantitatively assess their accuracy. This result will be tested for a practical SCAL test, determining how long a flow experiment needs to be run to accurately characterize the rate of energy dissipation by the flow. The outcome will be a best practice guide for the determination of relative permeability from core-scale experiments and/or digital core simulations that ensure the energy budget is fully accounted for in the relative permeability coefficient.

Mcclure, James [Virginia Tech, Blacksburg]↗
NASA NTRS1990

Modeling the pressure-strain correlation of turbulence: An invariant dynamical systems approach

The modeling of the pressure-strain correlation of turbulence is examined from a basic theoretical standpoint with a view toward developing improved second-order closure models. Invariance considerations along with elementary dynamical systems theory are used in the analysis of the standard hierarchy of closure models. In these commonly used models, the pressure-strain correlation is assumed to be a linear function of the mean velocity gradients with coefficients that depend algebraically on the anisotropy tensor. It is proven that for plane homogeneous turbulent flows the equilibrium structure of this hierarchy of models is encapsulated by a relatively simple model which is only quadratically nonlinear in the anisotropy tensor. This new quadratic model - the SSG model - is shown to outperform the Launder, Reece, and Rodi model (as well as more recent models that have a considerably more complex nonlinear structure) in a variety of homogeneous turbulent flows. Some deficiencies still remain for the description of rotating turbulent shear flows that are intrinsic to this general hierarchy of models and, hence, cannot be overcome by the mere introduction of more complex nonlinearities. It is thus argued that the recent trend of adding substantially more complex nonlinear terms containing the anisotropy tensor may be of questionable value in the modeling of the pressure-strain correlation. Possible alternative approaches are discussed briefly.

Speziale, Charles G.↗
NASA NTRS1991

Modelling the pressure-strain correlation of turbulence - An invariant dynamical systems approach

The modeling of the pressure-strain correlation of turbulence is examined from a basic theoretical standpoint with a view toward developing improved second-order closure models. Invariance considerations along with elementary dynamical systems theory are used in the analysis of the standard hierarchy of closure models. In these commonly used models, the pressure-strain correlation is assumed to be a linear function of the mean velocity gradients with coefficients that depend algebraically on the anisotropy tensor. It is proven that for plane homogeneous turbulent flows the equilibrium structure of this hierarchy of models is encapsulated by a relatively simple model which is only quadratically nonlinear in the anisotropy tensor. This new quadratic model - the SSG model - is shown to outperform the Launder, Reece, and Rodi model (as well as more recent models that have a considerably more complex nonlinear structure) in a variety of homogeneous turbulent flows. Some deficiencies still remain for the description of rotating turbulent shear flows that are intrinsic to this general hierarchy of models and, hence, cannot be overcome by the mere introduction of more complex nonlinearities. It is thus argued that the recent trend of adding substantially more complex nonlinear terms containing the anisotropy tensor may be of questionable value in the modeling of the pressure-strain correlation. Possible alternative approaches are discussed briefly.

Speziale, Charles G.↗
DOE OSTI2021

Spectroscopic factors in dripline nuclei

Single-nucleon knockout reaction studies of the proton-dripline nuclei 9 C and 13 O suggest an appreciable suppression of spectroscopic factors. In this work, we calculate the one-neutron and one-proton spectroscopic factors for the mirror pair 9 C – 9 Li and 13 O using two variants of the continuum shell model: the complex-energy Gamow shell model and the real-energy shell model embedded in the continuum. Furthermore, our results indicate that the continuum effects strongly suppress the spectroscopic factors of well-bound orbits in the dripline systems but have less impact on the spectroscopic factors of weakly bound states.

6 ≤ A ≤ 19↗
DOE OSTI2021

Stochastic density functional theory: Real- and energy-space fragmentation for noise reduction

Stochastic density functional theory (sDFT) is becoming a valuable tool for studying ground-state properties of extended materials. The computational complexity of describing the Kohn-Sham orbitals is replaced by introducing a set of random (stochastic) orbitals leading to linear and often sub-linear scaling of certain ground-state observables at the account of introducing a statistical error. Schemes to reduce the noise are essential, for example, for determining the structure using the forces obtained from sDFT. Recently, we have introduced two embedding schemes to mitigate the statistical fluctuations in the electron density and resultant forces on the nuclei. Both techniques were based on fragmenting the system either in real space or slicing the occupied space into energy windows, allowing for a significant reduction in the statistical fluctuations. For chemical accuracy, further reduction of the noise is required, which could be achieved by increasing the number of stochastic orbitals. However, the convergence is relatively slow as the statistical error scales as 1 / N χ according to the central limit theorem, where N χ is the number of random orbitals. In this paper, we combined the embedding schemes mentioned above and introduced a new approach that builds on overlapped fragments and energy windows. The new approach significantly lowers the noise for ground-state properties, such as the electron density, total energy, and forces on the nuclei, as demonstrated for a G-center in bulk silicon.

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

Thermal adaptation of soil microbial growth traits in response to chronic warming

ABSTRACT Adaptation of soil microbes due to warming from climate change has been observed, but it remains unknown what microbial growth traits are adaptive to warming. We studied bacterial isolates from the Harvard Forest Long-Term Ecological Research site, where field soils have been experimentally heated to 5°C above ambient temperature with unheated controls for 30 years. We hypothesized that Alphaproteobacteria from warmed plots have (i) less temperature-sensitive growth rates; (ii) higher optimum growth temperatures; and (iii) higher maximum growth temperatures compared to isolates from control plots. We made high-throughput measurements of bacterial growth in liquid cultures over time and across temperatures from 22°C to 37°C in 2–3°C increments. We estimated growth rates by fitting Gompertz models to the growth data. Temperature sensitivity of growth rate, optimum growth temperature, and maximum growth temperature were estimated by the Ratkowsky 1983 model and a modified Macromolecular Rate Theory (MMRT) model. To determine evidence of adaptation, we ran phylogenetic generalized least squares tests on isolates from warmed and control soils. Our results showed evidence of adaptation of higher optimum growth temperature of bacterial isolates from heated soils. However, we observed no evidence of adaptation of temperature sensitivity of growth and maximum growth temperature. Our project begins to capture the shape of the temperature response curves, but illustrates that the relationship between growth and temperature is complex and cannot be limited to a single point in the biokinetic range. IMPORTANCE Soils are the largest terrestrial carbon sink and the foundation of our food, fiber, and fuel systems. Healthy soils are carbon sinks, storing more carbon than they release. This reduces the amount of carbon dioxide released into the atmosphere and buffers against climate change. Soil microbes drive biogeochemical cycling and contribute to soil health through organic matter breakdown, plant growth promotion, and nutrient distribution. In this study, we determined how soil microbial growth traits respond to long-term soil warming. We found that bacterial isolates from warmed plots showed evidence of adaptation of optimum growth temperature. This suggests that increased microbial biomass and growth in a warming world could result in greater carbon storage. As temperatures increase, greater microbial activity may help reduce the soil carbon feedback loop. Our results provide insight on how atmospheric carbon cycling and soil health may respond in a warming world.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Efficient Hamiltonian encoding algorithms for extracting quantum control mechanism as interfering pathway amplitudes in the Dyson series

Hamiltonian encoding is a methodology for revealing the mechanism behind the dynamics governing controlled quantum systems. In this paper, following Mitra and Rabitz \cite{abhra_1}, we define mechanism via pathways of eigenstates that describe the evolution of the system, where each pathway is associated with a complex-valued amplitude corresponding to a term in the Dyson series. The evolution of the system is determined by the constructive and destructive interference of these pathway amplitudes. Pathways with similar attributes can be grouped together into pathway classes. The amplitudes of pathway classes are computed by modulating the Hamiltonian matrix elements and decoding the subsequent evolution of the system rather than by direct computation of the individual terms in the Dyson series. The original implementation of Hamiltonian encoding was computationally intensive and became prohibitively expensive in large quantum systems. This paper presents two new encoding algorithms that calculate the amplitudes of pathway classes by using techniques from graph theory and algebraic topology to exploit patterns in the set of allowed transitions, greatly reducing the number of matrix elements that need to be modulated. These new algorithms provide an exponential decrease in both computation time and memory utilization with respect to the Hilbert space dimension of the system. To demonstrate the use of these techniques, they are applied to two illustrative state-to-state transition problems.

Abrams, Erez [Princeton University, Massachusetts ↗
NASA NTRS1987

Space station integrated wall design and penetration damage control

The analysis code BUMPER executes a numerical solution to the problem of calculating the probability of no penetration (PNP) of a spacecraft subject to man-made orbital debris or meteoroid impact. The codes were developed on a DEC VAX 11/780 computer that uses the Virtual Memory System (VMS) operating system, which is written in FORTRAN 77 with no VAX extensions. To help illustrate the steps involved, a single sample analysis is performed. The example used is the space station reference configuration. The finite element model (FEM) of this configuration is relatively complex but demonstrates many BUMPER features. The computer tools and guidelines are described for constructing a FEM for the space station under consideration. The methods used to analyze the sensitivity of PNP to variations in design, are described. Ways are suggested for developing contour plots of the sensitivity study data. Additional BUMPER analysis examples are provided, including FEMs, command inputs, and data outputs. The mathematical theory used as the basis for the code is described, and illustrates the data flow within the analysis.

Coronado, A. R.↗
DOE OSTI2021

Multirate linearly-implicit GARK schemes

Many complex applications require the solution of initial-value problems where some components change fast, while others vary slowly. Multirate schemes apply different step sizes to resolve different components of the system, according to their dynamics, in order to achieve increased computational efficiency. The stiff components of the system, fast or slow, are best discretized with implicit base methods in order to ensure numerical stability. To this end, linearly implicit methods are particularly attractive as they solve only linear systems of equations at each step. This paper develops the Multirate GARK-ROS/ROW (MR-GARK-ROS/ROW) framework for linearly-implicit multirate time integration. The order conditions theory considers both exact and approximative Jacobians. The effectiveness of implicit multirate methods depends on the coupling between the slow and fast computations; an array of efficient coupling strategies and the resulting numerical schemes are analyzed. Multirate infinitesimal step linearly-implicit methods, that allow arbitrarily small micro-steps and offer extreme computational flexibility, are constructed. The new unifying framework includes existing multirate Rosenbrock(-W) methods as particular cases, and opens the possibility to develop new classes of highly effective linearly implicit multirate integrators.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Predicting the oxidation states of Mn ions in the oxygen-evolving complex of photosystem II using supervised and unsupervised machine learning

Abstract Serial Femtosecond Crystallography at the X-ray Free Electron Laser (XFEL) sources enabled the imaging of the catalytic intermediates of the oxygen evolution reaction of Photosystem II (PSII). However, due to the incoherent transition of the S-states, the resolved structures are a convolution from different catalytic states. Here, we train Decision Tree Classifier and K-means clustering models on Mn compounds obtained from the Cambridge Crystallographic Database to predict the S-state of the X-ray, XFEL, and CryoEM structures by predicting the Mn’s oxidation states in the oxygen-evolving complex. The model agrees mostly with the XFEL structures in the dark S 1 state. However, significant discrepancies are observed for the excited XFEL states (S 2 , S 3, and S 0 ) and the dark states of the X-ray and CryoEM structures. Furthermore, there is a mismatch between the predicted S-states within the two monomers of the same dimer, mainly in the excited states. We validated our model against other metalloenzymes, the valence bond model and the Mn spin densities calculated using density functional theory for two of the mismatched predictions of PSII. The model suggests designing a more optimized sample delivery and illumiation systems are crucial to precisely resolve the geometry of the advanced S-states to overcome the noncoherent S-state transition. In addition, significant radiation damage is observed in X-ray and CryoEM structures, particularly at the dangler Mn center (Mn4). Our model represents a valuable tool for investigating the electronic structure of the catalytic metal cluster of PSII to understand the water splitting mechanism.

Plant Sciences↗
DOE OSTI2021

Adiabatic quantum decoherence in many non-interacting subsystems induced by the coupling with a common boson bath

Highlights: • System–environment quantum correlation: a main solid state NMR decoherence channel. • Non-separable system–environment model yields realistic spin decoherence rates. • New open quantum system approach explains irreversible decay of refocused NMR echoes. • Adiabatic quantum decoherence is inherently irreversible and eigen-selective. This work addresses adiabatic quantum decoherence of many-body spin systems coupled with a boson field in the framework of open quantum systems theory. We generalize the traditional spin-boson model by considering a system–environment interaction Hamiltonian that represents a partition of non-interacting subsystems and highlights the collective correlation that appears exclusively due to the coupling with a common environment. Remarkably, this simple, exactly solvable model encompasses relevant aspects of a many-body open quantum system and features the subtle quantum effects that arise when the size scales up to a macroscopic level. We derive an analytical expression for the time dependence of the density matrix elements (in the preferred basis) without assuming coarse-graining. The resulting decoherence function is eigen-selective and is a complex exponential whose exponent has a real part that introduces a decay similar to that in the spin-boson model. On the contrary, the imaginary part depends on the quantum numbers and geometry of the whole partition and does not reflect the system temperature. Motivated by decoherence in solid-state NMR, and in search of realistic numerical estimations, we apply the theoretical results to a partition of dipole-coupled spin pairs in contact with a common phonon bath, using typical parameters of hydrated salts. The proposal allows estimating the decoherence time scale in terms of the system physical constants: sound velocity and eigenvalue distribution width. As a significant novelty, the decoherence function phase depends on the eigenvalue distribution throughout the sample. It plays the leading role, overshadowing the mechanism associated with the bath thermal state. Finally, we apply the formalism to describe decoherence in the “magic echo” NMR reversal experiment. We find that the system–environment correlation explains the origin of irreversibility, and both the decoherence rate value and its dependence on the dipolar frequency, are remarkably similar to the experiment.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Surface vs Homogeneous Organo-Hafnium Catalyst Ion-Pairing and Ligand Effects on Ethylene Homo- and Copolymerizations

Heterogeneous catalysts have long dominated polyethylene and polypropylene production, but understanding their catalysis is challenged by uncertainties in active site structures and percentages. Surface-bound organometallic catalysts are an emerging strategy to combine successful homogeneous catalysts having well-understood structures, relatively high percentages of active sites, and exceptional control of selectivity, with the attractions of surface catalysts, to transition promising homogeneous systems to large-scale heterogeneous ones. Nevertheless, surface-bound olefin polymerization catalysts typically produce ultrahigh M w ’s but with low activity and comonomer selectivity. Here, we report the systematic synthesis and characterization of a series of pyridylamido–Hf complexes and their corresponding surface catalysts chemisorbed on sulfated alumina (AlS) and zirconia (ZrS). Comparative ethylene homo- and 1-octene copolymerizations reveal similar activity and 1-octene selectivity trends in the homogeneous and heterogeneous systems. For the surface pyridylamido–Hf catalyst series, large variations in activity (up to 10×) and 1-octene incorporation (up to 28×) are achieved by ligand and support manipulation. Interestingly, while the homogeneous catalysts exhibit positive comonomer effects in ethylene/1-octene copolymerization, the surface catalysts behave oppositely. Extended X-ray absorption fine structure (EXAFS) reveals significantly elongated Hf···O bond distances vs typical Hf–O covalent bonds (2.06 vs 1.97 Å). Furthermore, density functional theory (DFT) analysis of the heterolytic ion pair separation enthalpies, olefin insertion energetics, and NBO/Bader charges also suggest electrostatic Hf cation–anionic support binding and catalytic patterns, which are modulated by the ion-pairing energetics and ligand architecture.

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

Deeply Virtual Compton Scattering at Hall A, Jefferson Lab

The Standard Model of particle physics de?nes quarks and leptons as the basic building blocks of all matter. The interaction between them are mediated by force carrying gauge bosons. Quantum Chromo Dynamics (QCD), the theory that explains the strong interaction is still not complete enough to derive the physical observables of a Quark-Gluon system from the fundamental degrees of freedom of its constituents. Experimentally observable single particle densities provide important insights into our understanding of the quark-gluon system and hence help fill in the gaps of QCD. Generalized Parton Distributions (GPDs) provide simultaneous information of both spacial and longitudinal momentum distributions of constituents of a quark-gluon system. Deeply Virtual Compton Scattering (DVCS) is understood to be the simplest and cleanest process to access GPDs. Even though the exclusive DVCS is simple to understand, the experimental process however, is complex with the Bethe-Heitler and e + p ? e + ? + N? being in the mix of the electron proton scattering e+p ? e+?+X. Over the years, 3 generations of DVCS experiments have been conducted in the Experimental Hall-A of Thomas Jefferson National Accelerator Facility (JLab). This thesis presents the extraction of DVCS cross section in 9 total kinematic points from the 3rd generation experiment (DVCS3) conducted after the 12 GeV upgrade of the Continuous Electron Beam Accelerator Facility (CEBAF) of JLab during Fall 2014 - Fall 2016.

Hashir Rashad, Mohamed Nuhman↗