Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density gradient function”

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 379 records · Page 21

Diffusion-Model-Assisted Supervised Learning of Generative Models for Density Estimation

Here, we present a supervised learning framework of training generative models for density estimation. Generative models, including generative adversarial networks (GANs), normalizing flows, and variational auto-encoders (VAEs), are usually considered as unsupervised learning models, because labeled data are usually unavailable for training. Despite the success of the generative models, there are several issues with the unsupervised training, e.g., requirement of reversible architectures, vanishing gradients, and training instability. To enable supervised learning in generative models, we utilize the score-based diffusion model to generate labeled data. Unlike existing diffusion models that train neural networks to learn the score function, we develop a training-free score estimation method. This approach uses mini-batch-based Monte Carlo estimators to directly approximate the score function at any spatial-temporal location in solving an ordinary differential equation (ODE), corresponding to the reverse-time stochastic differential equation (SDE). This approach can offer both high accuracy and substantial time savings in neural network training. Once the labeled data are generated, we can train a simple, fully connected neural network to learn the generative model in the supervised manner. Compared with existing normalizing flow models, our method does not require the use of reversible neural networks and avoids the computation of the Jacobian matrix. Compared with existing diffusion models, our method does not need to solve the reverse-time SDE to generate new samples. As a result, the sampling efficiency is significantly improved. We demonstrate the performance of our method by applying it to a set of 2D datasets as well as real data from the University of California Irvine (UCI) repository.

97 MATHEMATICS AND COMPUTING↗

A model for simulating random atmospheres as a function of latitude, season, and time

An empirical stochastic computer model was developed with the capability of generating random thermodynamic profiles of the atmosphere below an altitude of 99 km which are characteristic of any given season, latitude, and time of day. Samples of temperature, density, and pressure profiles generated by the model are statistically similar to measured profiles in a data base of over 6000 rocket and high-altitude atmospheric soundings; that is, means and standard deviations of modeled profiles and their vertical gradients are in close agreement with data. Model-generated samples can be used for Monte Carlo simulations of aircraft or spacecraft trajectories to predict or account for the effects on a vehicle's performance of atmospheric variability. Other potential uses for the model are in simulating pollutant dispersion patterns, variations in sound propagation, and other phenomena which are dependent on atmospheric properties, and in developing data-reduction software for satellite monitoring systems.

Campbell, J. W.↗

Low-Earth Orbit Trajectory Optimization in the Presence of Atmospheric Uncertainty

The previous 20 to 25 years have seen a tremendous increase in space exploration, and with that an increase in the level of logistics planning needed to ensure mission success. For spacecraft that are designed to be periodically re-supplied, a key logistics consumable is propellant, as it constitutes the greatest up-mass on re-supply vehicles. A trajectory design strategy is therefore desired that minimizes propellant usage in order to ease the demand for propellant re-supply missions. This thesis develops such a strategy in three stages, and uses the International Space Station (ISS) as its testbed, as no other LEO spacecraft is more challenging from a space logistics standpoint. First, the ISS trajectory planning problem is formulated as a constrained burn optimization problem assuming a deterministic atmosphere. The cost function is total ∆v, with constraints imposed on longitude of ascending viii node (LAN) and semi-major axis (SMA) altitude. Analytic derivatives are constructed for both the cost and constraints, which are necessary given the 6-week to 2-year time frames being considered. A gradient-based optimizer is then utilized to find locally-optimal solutions to real-world ISS trajectory planning problems. Second, atmospheric uncertainty is addressed by constructing a probabilistic model of space weather data using Gaussian Processes (GPs). Bayesian inference is performed using the GP model to generate mean and covariance estimates for space weather predictions, whose pedigree is assessed against test data. The predictions are then mapped into atmospheric density via the analytic Jacchia-Roberts density model, and the effect of space weather uncertainty on orbital lifetime is examined. Third, an ISS burn execution uncertainty model is developed. This model, along with the space weather uncertainty model, are deployed in a linear covariance analysis to ascertain their combined effect on LAN and SMA altitude dispersions. The deterministic constraints from the original problem are re-formulated as stochastic constraints, where now the constraint uncertainty interval is required to fall within specified bounds. An updated optimization framework is constructed using the original ∆v cost function along with the stochastic constraints to solve the trajectory optimization problem under atmospheric uncertainty. Finally, the complete architecture is summarized for deployment in an operational setting.

Trajectory Optimization↗

System noise analysis of the dumbbell tethered satellite for gravity-gradient measurements

An analysis of the dumbbell gravity gradiometer concept for measuring short wavelength variations in the earth's gravity gradient is presented. Variations in the gradient are recorded by measuring tension variations in a vertically stabilized satellite consisting of heavy masses connected by a long wire or rod. Tension noise arises from the excitation of various mechanical oscillations of the system. The principal noise sources that were identified are fluctuations in atmospheric drag heating and drag force resulting from density variations and winds. Approximate analytical expressions are presented for the tension noise as a function of the system design parameters for various possible configurations. Computer simulations using numerical integration were performed to study the tension noise for several sample cases. Three designs consistent with Shuttle launch capabilities are discussed.

Colombo, G.↗

Predicting Morphological Evolution during Coprecipitation of MnCO 3 Battery Cathode Precursors Using Multiscale Simulations Aided by Targeted Synthesis

The performance of lithium-ion batteries is intimately linked to both the structure and the morphology of the cathode material, which in turn is critically linked to the synthesis conditions. However, few studies focus on understanding synthesis, especially during the coprecipitation of metal oxide precursors, a process that largely determines the final morphology of the material. In this paper, we go beyond the typical equilibrium particle shape analysis conducted in the literature and incorporate kinetic aspects of morphology evolution. We perform these studies using controlled synthesis on a well-defined metal salt system (MnCO 3 ) combined with multiscale simulations and high-resolution microscopy. Results show that with increasing metal concentration, the particles transition from rhombohedral to cubic to spherical shapes. Computational analysis using density functional theory (DFT) reveals that rhombohedral shaped particles evolve under equilibrium conditions. Phase field techniques indicate that at higher metal concentrations, fast growth kinetics of the precipitates result in the transition to cubic and, subsequently, spherical shapes, accompanied by a decrease in particle size. This study, while limited to the one metal salt system, provides an approach to shed light on the synthesis process of mixed transition metal salts, gradient materials, and other cathode materials of interest to the battery community.

25 ENERGY STORAGE↗

The Influence of Tungsten on the Chemical Composition of a Temporally Evolving Nanostructure of a Model Ni-Al-Cr Superalloy

The influence of W on the temporal evolution of gamma' precipitation toward equilibrium in a model Ni-Al-Cr alloy is investigated by three-dimensional atom-probe (3DAP) microscopy and transmission electron microscopy (TEM). We report on the alloys Ni-10 Al-8.5 Cr (at.%) and Ni-10 Al-8.5 Cr-2 W (at.%), which were aged isothermally in the gamma + gamma' two-phase field at 1073 K, for times ranging from 0.25 to 264 h. Spheroidal-shaped gamma' precipitates, 5-15 nm diameter, form during quenching from above the solvus temperature in both alloys at a high number density (approx. 10(exp 23/cu m). As gamma' precipitates grow with aging at 1073 K, a transition from spheriodal- to cuboidal-shaped precipitates is observed in both alloys. The elemental partitioning and spatially resolved concentration profiles across the gamma' precipitates are obtained as a function of aging time from three-dimensional atom-by-atom reconstructions. Proximity histogram concentration profiles of the quaternary alloy demonstrate that W concentration gradients exist in gamma' precipitates in the as-quenched and 0.25-h aging states, which disappear after 1 h of aging. The diffusion coefficient of W in gamma' is estimated to be 6.2 x 10(exp -20) sq m/s at 1073 K. The W addition decreases the coarsening rate constant, and leads to stronger partitioning of Al to gamma' and Cr to gamma.

Sudbrack, Chantal K.↗

Tailoring the Physicochemical Properties of Nb Thin Films via Surface Engineering Methods

The modification of surface oxide layers formed on niobium (Nb) thin films via chemical mechanical planarization (CMP) and accelerated neutral atom beam (ANAB) processing provides a promising route toward tailoring their emergent properties and performance when used as superconducting qubits. Here, in this study, we show that CMP- and ANAB-formed Nb oxides are significantly thinner and smoother than the native oxide, as revealed by transmission electron microscopy (TEM) and atomic force microscopy. Scanning TEM and energy-dispersive X-ray spectroscopy along with X-ray photoelectron spectroscopy identified an oxidation gradient within the native and surface-engineered oxides. The topside layer is dominated by Nb 5+ (Nb 2 O 5 ), with various Nb suboxides present closer to the oxide/metal interface. Time-of-flight secondary ion mass spectrometry (ToF-SIMS) depth profiling confirmed the presence of an oxygen content gradient and demonstrated the enhanced resistance of the CMP- and ANAB-formed oxides to oxygen surface exchange and subsequent diffusion via 18 O 2 isotopic labeling experiments. ToF-SIMS also identified an interfacial layer containing trapped hydrogen (H)-containing species at the Nb oxide/metal interface. In situ ToF-SIMS and TEM revealed migration of the H/OH interfacial layer coinciding with decomposition of the surface oxide. Furthermore, our density functional theory calculations indicated that both H from moisture present in ambient air and bulk H in Nb films tend to segregate at the interface. These findings underscore the importance of understanding surface oxidation mechanisms, hydrogen incorporation, and their impact on the designed functionalities of Nb-based devices.

36 MATERIALS SCIENCE↗

Compressible boundary layer with normal pressure gradients: Quasi-similarity equations - Their properties at the wall and at sharp and blunt leading edges.

The failure of most viscous-inviscid interaction methods at strong interactions is attributed to the presence of a normal pressure gradient. A new theory is proposed for supersonic laminar boundary layers that can generate normal pressure gradients. The Navier-Stokes equations are reexamined by an order of magnitude analysis and all first and second order terms are retained. The approximation is found to be dependent not only on the boundary layer thickness but also on the ratio of the dimensionless viscosity and density. The equations are transformed into two quasi-similar, nonlinear, third order, ordinary integro-differential equations for the velocity and pressure as functions of a single transverse variable. The properties of the equations at the boundaries are discussed.

Seginer, A.↗

Solar Wind Helium Abundance as a Function of Speed and Heliographic Latitude: Variation through a Solar Cycle

We present a study of the variation of the relative abundance of helium to hydrogen in the solar wind as a function of solar wind speed and heliographic latitude over the previous solar cycle. The average values of A(sub He), the ratio of helium to hydrogen number densities, are calculated in 25 speed intervals over 27-day Carrington rotations using Faraday Cup observations from the Wind spacecraft between 1995 and 2005. The higher speed and time resolution of this study compared to an earlier work with the Wind observations has led to the discovery of three new aspects of A(sub He), modulation during solar minimum from mid-1995 to mid-1997. First, we find that for solar wind speeds between 350 and 415 km/s, A(sub He), varies with a clear six-month periodicity, with a minimum value at the heliographic equatorial plane and a typical gradient of 0.01 per degree in latitude. For the slow wind this is a 30% effect. We suggest that the latitudinal gradient may be due to an additional dependence of coronal proton flux on coronal field strength or the stability of coronal loops. Second, once the gradient is subtracted, we find that A(sub He), is a remarkably linear function of solar wind speed. Finally, we identify a vanishing speed, at which A(sub He), is zero, is 259 km/s and note that this speed corresponds to the minimum solar wind speed observed at one AU. The vanishing speed may be related to previous theoretical work in which enhancements of coronal helium lead to stagnation of the escaping proton flux. During solar maximum the A(sub He), dependences on speed and latitude disappear, and we interpret this as evidence of two source regions for slow solar wind in the ecliptic plane, one being the solar minimum streamer belt and the other likely being active regions.

Kasper, J. C.↗

Force-Free Identification of Minimum-Energy Pathways and Transition States for Stochastic Electronic Structure Theories

Here, the accurate mapping of potential energy surfaces (PESs) is crucial to our understanding of the numerous physical and chemical processes mediated by atomic rearrangements, such as conformational changes and chemical reactions, and the thermodynamic and kinetic feasibility of these processes. Stochastic electronic structure theories, e.g., Quantum Monte Carlo (QMC) methods, enable highly accurate total energy calculations that in principle can be used to construct the PES. However, their stochastic nature poses a challenge to the computation and use of forces and Hessians, which are typically required in algorithms for minimum-energy pathway (MEP) and transition state (TS) identification, such as the nudged elastic band (NEB) algorithm and its climbing image formulation. Here, we present strategies that utilize the surrogate Hessian line-search method, previously developed for QMC structural optimization, to efficiently identify MEP and TS structures without requiring force calculations at the level of the stochastic electronic structure theory. By modifying the surrogate Hessian algorithm to operate in path-orthogonal subspaces and at saddle points, we show that it is possible to identify MEPs and TSs by using a force-free QMC approach. We demonstrate these strategies via two examples, the inversion of the ammonia (NH 3 ) molecule and the nucleophilic substitution (S N 2) reaction F – + CH 3 F → FCH 3 + F – . We validate our results using Density Functional Theory (DFT)- and Coupled Cluster (CCSD, CCSD(T))-based NEB calculations. We then introduce a hybrid DFT-QMC approach to compute thermodynamic and kinetic quantities, free energy differences, rate constants, and equilibrium constants that incorporates stochastically optimized structures and their energies, and show that this scheme improves upon DFT accuracy. Our methods generalize straightforwardly to other systems and other high-accuracy theories that similarly face challenges computing energy gradients, paving the way for highly accurate PES mapping, transition state determination, and thermodynamic and kinetic calculations at significantly reduced computational expense.

Iyer, Gopal R.↗

Embedded training of neural-network subgrid-scale turbulence models

We report that the weights of a deep neural-network model are optimized in conjunction with the governing flow equations to provide a model for subgrid-scale stresses in a temporally developing plane turbulent jet at Reynolds number Re 0 = 6000 . The objective function for training is first based on the instantaneous filtered velocity fields from a corresponding direct numerical simulation, and the training is by a stochastic gradient descent method, which uses the adjoint Navier-Stokes equations to provide the end-to-end sensitivities of the model weights to the velocity fields. In-sample and out-of-sample testing on multiple dual-jet configurations show that its required mesh density in each coordinate direction for prediction of mean flow, Reynolds stresses, and spectra is half that needed by the dynamic Smagorinsky model for comparable accuracy. The same neural-network model trained directly to match filtered subgrid-scale stresses, without the constraint of being embedded within the flow equations during the training, fails to provide a qualitatively correct prediction. The coupled formulation is generalized to train based only on mean-flow and Reynolds stresses, which are more readily available in experiments. The mean-flow training provides a robust model, which is important, though a somewhat less accurate prediction for the same coarse meshes, as might be anticipated due to the reduced information available for training in this case. The anticipated advantage of the formulation is that the inclusion of resolved physics in the training increases its capacity to extrapolate. This is assessed for the case of passive scalar transport, for which it outperforms established models due to improved mixing predictions.

42 ENGINEERING↗

Interchange Method in Compressible Magnetized Couette Flow: Magnetorotational and Magnetoconvective Instabilities

We obtain the general forms of the axisymmetric stability criteria in a magnetized compressible Couette flow using an energy variational principle, the so-called interchange or Chandrasekhar s met hod, which we applied successfully in the incompressible case. This formulation accounts for the simultaneous presence of gravity, rotation, a toroidal magnetic field, a weak axial magnetic field, entropy gradients, and density gradients in the initial equilibrium state. The power of the method lies in its simplicity which allows us to derive extremely compact and physically clear expressions for the relevant stability criteria despite the inclusion of so many physical effects. In the implementation of the method, all the applicable conservation laws are explicitly taken into account during the variations of a quantity with dimensions of energy which we call the free energy function. As in the incompressible case, the presence of an axial field invalidates the conservation laws of angular momentum and azimuthal magnetic flux and introduces instead isorotation and axial current conservation along field lines. Our results are therefore markedly different depending on whether an axial magnetic field is present, and generalize in two simple expressions all previously known, partial stability criteria for the appearance of magnetorotational instability. Furthermore, the coupling between magnetic tension and buoyancy and its influence to the dynamics of nonhomoentropic magnetized flows becomes quite clear from our results. In the limits of plane-parallel atmospheres and homoentropic flows, our formulation easily recovers the stability criteria for suppression of convective and Parker instabilities, as well as some related special cases studied over 40 years ago by Newcomb and Tserkovnikov via laborious variational techniques.

Christodoulou, Dimitris M.↗

Genomic fingerprints of the world’s soil ecosystems

Despite the explosion of soil metagenomic data, we lack a synthesized understanding of patterns in the distribution and functions of soil microorganisms. These patterns are critical to predictions of soil microbiome responses to climate change and resulting feedbacks that regulate greenhouse gas release from soils. To address this gap, we assay 1,512 manually curated soil metagenomes using complementary annotation databases, read-based taxonomy, and machine learning to extract multidimensional genomic fingerprints of global soil microbiomes. Our objective is to uncover novel biogeographical patterns of soil microbiomes across environmental factors and ecological biomes with high molecular resolution. We reveal shifts in the potential for (i) microbial nutrient acquisition across pH gradients; (ii) stress-, transport-, and redox-based processes across changes in soil bulk density; and (iii) greenhouse gas emissions across biomes. We also use an unsupervised approach to reveal a collection of soils with distinct genomic signatures, characterized by coordinated changes in soil organic carbon, nitrogen, and cation exchange capacity and in bulk density and clay content that may ultimately reflect soil environments with high microbial activity. Genomic fingerprints for these soils highlight the importance of resource scavenging, plant-microbe interactions, fungi, and heterotrophic metabolisms. Across all analyses, we observed phylogenetic coherence in soil microbiomes—more closely related microorganisms tended to move congruently in response to soil factors. Collectively, the genomic fingerprints uncovered here present a basis for global patterns in the microbial mechanisms underlying soil biogeochemistry and help beget tractable microbial reaction networks for incorporation into process-based models of soil carbon and nutrient cycling.

59 BASIC BIOLOGICAL SCIENCES↗

Charge Separation in Monolayer WSe 2 by Strain Engineering: Implications for Strain-Induced Diode Action

Strain-engineering band structure in transition-metal dichalcogenides (TMDC) is a promising avenue toward capabilities in optoelectronics. For example, controlling the flow of optically generated quasiparticles can be achieved by a localized strain field which reduces the bandgap and generates an energy-band gradient that funnels neutral excitons to the strain apex. It would be even more advantageous to mimic a diode’s internal field, where both conduction and valence bands bend in the same direction, to separate electrons and holes. This can be achieved if the strain in the TMDC layer lowers both the conduction band minimum as well as the valence band maximum during strain-induced band narrowing. Here, we have used density functional theory (DFT) calculations of monolayer WSe 2 electronic structure under biaxial strain to show that WSe 2 has this property. In this work, to test the band bending experimentally, we combined localized strain with electrostatic doping to follow photoluminescence from excitons and positive or negative trions. In unstrained WSe 2 , both positive and negative trion emissions dominate over excitons away from charge neutrality. In contrast, for strained areas, negative trions accumulate, while positive trion emission is near zero away from charge neutrality, indicating a lack of holes. Hence, strain bends both conduction and valence bands down, similarly to the band bending in a PN-diode depletion region, providing an opportunity to separate electrons and holes via localized strain.

2D materials↗

An analysis of the Voyager 2 Ultraviolet Spectrometer occultation data at Uranus - Inferring heat sources and model atmospheres

Heat source information is derived here from the Voyager 2 Ultraviolet Spectrometer occultation data of Uranus. Analytic functions for the local heat dependence on altitude are used to obtain a temperature profile by solving the heat equation. The stellar entrance and exit occultation and a solar occultation are used to infer the thermal and density structure of the atmosphere. The least squares fit solution to the solar occultation data gives one source located at 1.8 x 10 exp -5 microbar with a strength of 0.056 +/- 0.01 erg/sq cm/s. Latitudinal temperature gradients are obtained.

Stevens, Michael H.↗

A predictive strain-gradient model with no undetermined constants or length scales

A general meso–scale (GM) crystal plasticity (CP) model was developed that accounts for lower-order (strain hardening) and higher-order (internal stress) effects of geometrically necessary dislocations (GNDs). It is predictive: no arbitrary parameters or length scales were invoked and no ad hoc numerical techniques were employed. It uses general stress field equations for GND content and a novel harmonization technique to enforce consistency of elastic long-range singular defect fields with applied elastic-plastic fields. The model facilitates implementation in commercial finite element programs without requiring special elements, special boundary conditions, or access to element shape functions. GM simulations confirmed, with improved accuracy, previously published predictions of the Hall-Petch effect, Bauschinger effect, and anelasticity. Previously unpredicted phenomena were also predicted: anelasticity and hysteresis for single Ta crystals and strain-hardening stagnation. The internal stresses (higher-order effect) dominate at large length scales, while at small length scales, the GND density hardening (lower-order effect) dominates. As a result, GM predicts that strain heterogeneity and consequent GND internal stresses are important factors in anelasticity.

36 MATERIALS SCIENCE↗

Highly Rectifying Fluidic Diodes Based on Asymmetric Layer-by-Layer Nanofilms on Nanochannel Membranes

Nanochannel-based fluidic diodes display ion selectivity and ion-current rectification (ICR), which may prove important in energy harvesting and biosensors. This paper reports asymmetric functionalization of the outer surface of a flexible nanochannel polymer membrane to create fluidic diodes that give ICR. Layer-by-layer (LbL) adsorption with cross-linking of only the underlying part of the polyelectrolyte nanofilm leads to a porosity step across the film. The combination of a high effective surface-charge density and the porosity step in the film leads to a remarkable maximum ICR factor of ~200 with a pH gradient across the film. Incorporation of pH-sensitive polyelectrolyte components enables the ICR factor to increase an order of magnitude on going from pH 8 to pH 3. Furthermore, the coated membrane shows excellent anion-selectivity. Thus, LbL adsorption with partial cross-linking provides a simple method for creating coated nanochannel membranes that serve as pH-responsive ionic diodes for potential chemical/biosensors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Random Phase Approximation Correlation Energy Using Real-Space Density Functional Perturbation Theory

We present a real-space method for computing the random phase approximation (RPA) correlation energy within Kohn–Sham density functional theory, leveraging the low-rank nature of the frequency-dependent density response operator. In particular, we employ a cubic-scaling formalism based on density functional perturbation theory that circumvents the calculation of the response function matrix, instead relying on the ability to compute its product with a vector through the solution of the associated Sternheimer linear systems. We develop a large-scale parallel implementation of this formalism using the subspace iteration method in conjunction with the spectral quadrature method while employing the Kronecker product-based method for the application of the Coulomb operator and the conjugate orthogonal conjugate gradient method for the solution of the linear systems. We demonstrate convergence with respect to key parameters and verify the method’s accuracy by comparing with plane-wave results. We show that the framework achieves good strong scaling to many thousands of processors, reducing the time to solution for a lithium hydride system with 128 electrons to around 150 s on 4608 processors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗