Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Analytic gradients”

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 811 records · Page 45

Thermal Gradient Effects on Local Hotspot Ignition in 1,3,5,7‐Tetranitro‐1,3,5,7‐tetrazocane (HMX)

Understanding hotspot ignition, growth, and criticality, as well as the timescales of each, is crucial for parameterizing mesoscale and continuum‐level models that rely on a statistical understanding of hotspots. However, these models often consider hotspots to have uniform temperatures or for hotspots of a given temperature and size to always behave the same. Therefore, using molecular dynamics simulations, we assess the influence of thermal distribution effects on hotspot local ignition and time to ignition in 1,3,5,7‐tetranitro‐1,3,5,7‐tetrazocane (HMX) nanoscale hotspots. Finally, by assessing hotspots with a gradient driven from an initial two‐temperature core–shell setup, we show that small increases in the core temperature under a constant average temperature can lead to order‐of‐magnitude effects on reaction and local ignition timescales.

36 MATERIALS SCIENCE↗

A molecular dynamics simulation study of thermal conductivity anisotropy in β-octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (β-HMX)

Molecular dynamics (MD) simulations were used to predict the thermal conductivity of β-octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (β-HMX) along directions normal to the (011), (110), and (010) crystal planes. These directions were selected based on the measured morphological importance of the corresponding crystal surfaces. A reverse non-equilibrium MD approach was used wherein a constant heat flux is imposed along a prescribed direction and the resulting steady-state temperature gradient determined. The coefficient of thermal conductivity λ is the quotient of heat flux and temperature gradient (i.e. Fourier's law). Finite-size effects and sensitivity to imposed heat flux were investigated. The results reveal a modest dependence of the conductivity on crystal orientation, significant finite-size effects, and low sensitivity to imposed flux so long as the Fourier's law analysis is limited to the spatial interval in the simulation cell for which the temperature gradient is constant. Infinite-length thermal conductivities were estimated by linear regression of λ -1 (L) versus reciprocal cell length L -1 for each direction. The predicted values are systematically larger, but within a factor of two, than most published experimental determinations, the latter of which were obtained for pressed-powder or composite samples rather than oriented single crystals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Persistent compositions of non-stoichiometric compounds with low bulk diffusivity: A theory and application to Nb3Sn superconductors

Non-stoichiometric compounds may develop a composition gradient when they are formed by reactive diffusional processes. This paper reports an interesting phenomenon that in compounds with low bulk diffusivities, which rely mainly on grain boundary diffusion for their growth, the final bulk compositions may be far from equilibrium, with a very low bulk diffusivity leading to fixed (persistent) compositions – one such example is Nb 3 Sn, a superconductor. We investigated the microchemistry at the reactive interface using atom probe tomography to clarify the diffusion reaction mechanism for this low-bulk-diffusivity case and thus propose a theory for what determines the compound composition, using Nb 3 Sn as an example for concreteness. Using certain approximations, we derive an explicit analytical equation that illustrates what factors determine its composition profile. We compare our model with the known facts of Nb 3 Sn and see good agreement. In particular, this model predicts that internal oxidation may lead to higher Sn contents than conventional, non-oxidized Nb 3 Sn. Our measurements show that this is indeed true, and that the internally-oxidized Nb 3 Sn also has higher upper critical fields, achieving up to 28.2 T at 4.2 K. We discuss the general applicability of this model to non-stoichiometric compounds with low bulk diffusivity, and propose it as a tool to help in the design and processing of such materials for compositional control.

36 MATERIALS SCIENCE↗

Radiation-induced defects in SiO2 as determined with XPS

An analytical method based on oxide/substrate intensity ratios is used to deduce the product of the atomic number density and electron mean free path as a function of depth, for the cases of device-quality gate oxides grown on Si substrates and irradiated with zero to 20 eV electrons during in situ X-ray photoelectron spectroscopy. These structures had been thinned stepwise to 25-60 A by means of a wet chemical depth-profiling procedure. Si(3+) species are formed in the Si/SiO2 interface region, and their observation during their relaxation and annihilation is correlated with strained Si-O-Si bonds. The observation of bond of cleavage and bond strain gradients in these samples is used to extend silica devitrification models for the explanation of fixed oxide charge generation and interface states.

Grunthaner, F. J.↗

An analytical investigation of transient effects on rewetting of heated thin flat plates

The rewetting of a hot surface is a problem of prime importance in the microgravity application of heat pipe technology, where rewetting controls the time before operations can be re-established following depriming of a heat pipe. Rewetting is also important in the nuclear industry (in predicting behavior during loss-of-coolant accidents), as well as in the chemical and petrochemical industries. Recently Chan and Zhang have presented a closed-form solution for the determination of the rewetting speed of a liquid film flowing over a finite (but long) hot plate subject to uniform heating. Unfortunately, their physically unreasonable initial conditions preclude a meaningful analysis of start-up transient behavior. A new nondimensionalization and closed-form solution for an infinitely-long, uniformly-heated plate is presented. Realistic initial conditions (step change in temperature across the wetting front) and boundary conditions (no spatial temperature gradients infinitely far from the wetting front) are employed. The effects of parametric variation on the resulting simpler closed-form solution are presented and compared with the predictions of a 'quasi-steady' model. The time to reach steady-state rewetting is found to be a strong function of the initial dry-region plate temperature. For heated plates it is found that in most cases the effect of the transient response terms cannot be neglected, even for large times.

Platt, J. A.↗

Atomically Local Electric Field Induced Interface Water Reorientation for Alkaline Hydrogen Evolution Reaction

The slow water dissociation process in alkaline electrolyte severely limits the kinetics of HER. The orientation of H 2 O is well known to affect the dissociation process, but H 2 O orientation is hard to control because of its random distribution. Herein, an atomically asymmetric local electric field was designed by IrRu dizygotic single-atom sites (IrRu DSACs) to tune the H 2 O adsorption configuration and orientation, thus optimizing its dissociation process. The electric field intensity of IrRu DSACs is over 4.00×10 10 N/C. The ab initio molecular dynamics simulations combined with in situ Raman spectroscopy analysis on the adsorption behavior of H 2 O show that the M–H bond length (M=active site) is shortened at the interface due to the strong local electric field gradient and the optimized water orientation promotes the dissociation process of interfacial water. This work provides a new way to explore the role of single atomic sites in alkaline hydrogen evolution reaction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Heterogeneity and Hysteresis in the Polymer Collapse of Single Core–Shell Stimuli-Responsive Plasmonic Nanohybrids

Broad application of polymeric stimuli-responsive smart nanohybrids requires understanding the mechanisms governing active control. Ensemble techniques have identified inhomogeneous polymer collapse in microgels that potentially arise from heterogeneous interchain interactions and differences in core size. A single-particle examination would establish the influence of core size and internal polymer network heterogeneity on local interactions that contribute to the observed inhomogeneous polymer collapse dynamics of nanohybrids. Using single-particle dark-field spectroscopy, we investigated the complex polymer collapse profiles of core–shell plasmonic nanohybrids comprising thermoresponsive poly(N-isopropylacrylamide) (pNIPAM)-encapsulated gold nanorods (AuNRs). We report that the polymer collapse behavior was independent of the core size. For thinner polymer shells, we observed hysteresis in the collapse of AuNR@pNIPAMs, likely related to local pNIPAM aggregation due to interchain hydrogen bonding. For thicker polymer shells, we observed a broad polymer collapse distribution that we attributed to a two-step phase transition that arises from a polymer network density gradient. Furthermore, our single-particle approach relates the internal heterogeneity of the polymer network of nanohybrids to the mechanisms underlying heterogeneous phase transitions that traditional, ensemble-averaged approaches are unable to discern.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Solid-state polymer-particle hybrid electrolytes: Structure and electrochemical properties

Solid-state electrolytes (SSEs) are challenged by complex interfacial chemistry and poor ion transport through the interfaces they form with battery electrodes. Here, we investigate a class of SSE composed of micrometer-sized lithium oxide (Li 2 O) particles dispersed in a polymerizable 1,3-dioxolane (DOL) liquid. Ring-opening polymerization (ROP) of the DOL by Lewis acid salts inside a battery cell produces polymer-inorganic hybrid electrolytes with gradient properties on both the particle and battery cell length scales. These electrolytes sustain stable charge-discharge behavior in Li||NCM811 and anode-free Cu||NCM811 electrochemical cells. On the particle length scale, Li 2 O retards ROP, facilitating efficient ion transport in a fluid-like region near the particle surface. On battery cell length scales, gravity-assisted settling creates physical and electrochemical gradients in the hybrid electrolytes. By means of electrochemical and spectroscopic analyses, we find that Li 2 O particles participate in a reversible redox reaction that increases the effective CE in anode-free cells to values approaching 100%, enhancing battery cycle life.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Positron Beam Loading and Acceleration in the Blowout Regime of a Plasma Wakefield Accelerator

Plasma wakefield acceleration in the nonlinear blowout regime has achieved marked milestones in electron beam acceleration, demonstrating high acceleration gradients and energy efficiency while preserving excellent beam quality. However, this regime is deemed unsuitable for achieving positron acceleration of comparable results, which is vital for future compact electron–positron colliders. In this article, we find that an intense positron beam loaded at the back of beam-driven blowout cavity can self-consistently induce the focusing field and flatten the longitudinal wakefield, leading to stable, high-efficiency, and high-quality positron acceleration. This is achieved through the formation of an on-axis electron filament induced by positron beam load, which shapes the plasma wakefield in a distinct way compared to electron beam load in the blowout regime. Via a nonlinear analytic model and numerical simulations, we explain the novel beam loading effects of the interaction between the on-axis filament and the blowout cavity. High-fidelity simulations show that a high-charge positron beam can be accelerated with >20% energy transfer efficiency, ~1% energy spread, and ~1 mm·mrad normalized emittance, while considerably depleting the energy of the drive beam. The concept can also be extended to simultaneous acceleration of electron and positron beams and high transformer ratio positron acceleration as well. This development offers a new route for the application of plasma wakefield acceleration into particle physics.

Science & Technology - Other Topics↗

Retinoid quantification by HPLC/MS(n)

Retinoic acid (RA) mediates most of the biological effects of vitamin A that are essential for vertebrate survival. It acts through binding to receptors that belong to the nuclear receptor transcription factor superfamily (Mangelsdorf et al. 1994). It is also a highly potent vertebrate teratogen. To determine the function and effects of endogenous and exogenous RA, it is important to have a highly specific, sensitive, accurate, and precise analytical procedure. Current analyses of RA and other retinoids are labor intensive, of poor sensitivity, have limited specificity, or require compatibility with RA reporter cell lines (Chen et al. 1995. BIOCHEM: Pharmacol. 50: 1257-1264; Creech Kraft et al. 1994. BIOCHEM: J. 301: 111-119; Lanvers et al. 1996. J. Chromatogr. B Biomed. Appl. 685: 233-240; Maden et al. 1998. DEVELOPMENT: 125: 4133-4144; Wagner et al. 1992. DEVELOPMENT: 116: 55-66). This paper describes an HPLC/mass spectrometry/mass spectrometry product ion scan (HPLC/MS(n)) procedure for the analysis of retinoids that employs atmospheric pressure chemical ionization MS. The retinoids are separated by normal-phase column chromatography with a linear hexane-isopropanol-dioxane gradient. Each retinoid is detected by a unique series of MS(n) functions set at optimal collision-induced dissociation energy (30% to 32%) for all MS(n) steps. The scan events are divided into three segments, based on HPLC elution order, to maximize the mass spectrometer duty cycle. The all-trans, 9-cis, and 13-cis RA isomers are separated, if desired, by an isocratic hexane-dioxane-isopropanol mobile phase. This paper describes an HPLC/MS(n) procedure possessing high sensitivity and specificity for retinoids.

NASA Discipline Developmental Biology↗

Cyclic plasticity and failure of structural components

An analytical approach for low-cycle fatigue prediction is presented. The approach combines a cyclic plasticity model with the finite element method and a damage accumulation criterion for ductile metals. The cyclic plasticity model is based on the concept of the combination of several yield surfaces. The surfaces are related to the material uniaxial stress-strain curve idealized by piecewise linear segments. The damage criterion is based on the Coffin-Manson formulae modified for the mean stress variation effect. It is extended to the multiaxial varying stress-strain field and applied for both the crack initiation and the crack growth processes. The stable slow crack growth rate is approximated by the damage accumulation gradient computed from the cracked finite element models. This procedure requires fatigue testing data of only smooth specimens under constant strain amplitudes. The present approach is illustrated by numerical examples of an aircraft wing stiffened panel subjected to compression, which causes material yielding and residual tension.

Kalev, I.↗

Thermomechanical loading of multilayered cylindrical geometries in thermal cycling from 300 to 1300 K

The principal of multiple material layering is well known as an effective method of reducing heat transfer; however, thermal gradients can impose significant loads and lead to delamination and subsequent component failure. An analysis is developed and experimental data are discussed for the thermomechanical effects of multilayered materials on a heat sink substrate of cylindrical geometry subject to thermal cycling. The geometry is heated in cross-flow by a high-velocity flame and cooled in crossflow by ambient-temperature air from a critical flow orifice. Each layer of material possesses a threshold beyond which small changes in temperature or mechanical loading greatly influence the life in thermal cycling of the layered materials. Comparisons are made between the thermomechanical loads predicted by various numerical codes for the linear case and by a simplified analytic model.

Hendricks, R. C.↗

Study of Magnetic Damping Effect on Convection and Solidification Under G-Jitter Conditions

As shown by NASA resources dedicated to measuring residual gravity (SAMS and OARE systems), g-jitter is a critical issue affecting space experiments on solidification processing of materials. This study aims to provide, through extensive numerical simulations and ground based experiments, an assessment of the use of magnetic fields in combination with microgravity to reduce the g-jitter induced convective flows in space processing systems. We have so far completed asymptotic analyses based on the analytical solutions for g-jitter driven flow and magnetic field damping effects for a simple one-dimensional parallel plate configuration, and developed both 2-D and 3-D numerical models for g-jitter driven flows in simple solidification systems with and without presence of an applied magnetic field. Numerical models have been checked with the analytical solutions and have been applied to simulate the convective flows and mass transfer using both synthetic g-jitter functions and the g-jitter data taken from space flight. Some useful findings have been obtained from the analyses and the modeling results. Some key points may be summarized as follows: (1) the amplitude of the oscillating velocity decreases at a rate inversely proportional to the g-jitter frequency and with an increase in the applied magnetic field; (2) the induced flow approximately oscillates at the same frequency as the affecting g-jitter, but out of a phase angle; (3) the phase angle is a complicated function of geometry, applied magnetic field, temperature gradient and frequency; (4) g-jitter driven flows exhibit a complex fluid flow pattern evolving in time; (5) the damping effect is more effective for low frequency flows; and (6) the applied magnetic field helps to reduce the variation of solutal distribution along the solid-liquid interface. Work in progress includes numerical simulations and ground-based measurements. Both 2-D and 3-D numerical simulations are being continued to obtain further information on g-jitter driven flows and magnetic field effects. A physical model for ground-based measurements is completed and some measurements of the oscillating convection are being taken on the physical model. The comparison of the measurements with numerical simulations is in progress. Additional work planned in the project will also involve extending the 2-D numerical model to include the solidification phenomena with the presence of both g-jitter and magnetic fields.

Li, Ben Q.↗

Three-dimensional analytical model for the spatial variation of the foreshock electron distribution function - Systematics and comparisons with ISEE observations

A model which is consistent with the solar wind and shock surface boundary conditions for the foreshock electron distribution in the absence of wave-particle effects is formulated for an arbitrary location behind the magnetic tangent to the earth's bow shock. Variations of the gyrophase-averaged velocity distribution are compared and contrasted with in situ ISEE observations. It is found that magnetic mirroring of solar wind electrons is the most important process by which nonmonotonic reduced electron distributions in the foreshock are produced. Leakage of particles from the magnetosheath is shown to be relatively unimportant in determining reduced distributions that are nonmonotonic. The two-dimensional distribution function off the magnetic field direction is the crucial contribution in producing reduced distributions which have beams. The time scale for modification of the electron velocity distribution in velocity space can be significantly influenced by steady state spatial gradients in the background imposed by the curved shock geometry.

Fitzenreiter, R. J.↗

Modeling Liquid Water by Climbing up Jacob’s Ladder in Density Functional Theory Facilitated by Using Deep Neural Network Potentials

Within the framework of Kohn–Sham density functional theory (DFT), the ability to provide good predictions of water properties by employing a strongly constrained and appropriately normed (SCAN) functional has been extensively demonstrated in recent years. Here, we further advance the modeling of water by building a more accurate model on the fourth rung of Jacob’s ladder with the hybrid functional, SCAN0. In particular, we carry out both classical and Feynman path-integral molecular dynamics calculations of water with the SCAN0 functional and the isobaric–isothermal ensemble. To generate the equilibrated structure of water, a deep neural network potential is trained from the atomic potential energy surface based on ab initio data obtained from SCAN0 DFT calculations. For the electronic properties of water, a separate deep neural network potential is trained by using the Deep Wannier method based on the maximally localized Wannier functions of the equilibrated trajectory at the SCAN0 level. The structural, dynamic, and electric properties of water were analyzed. The hydrogen-bond structures, density, infrared spectra, diffusion coefficients, and dielectric constants of water, in the electronic ground state, are computed by using a large simulation box and long simulation time. For the properties involving electronic excitations, we apply the GW approximation within many-body perturbation theory to calculate the quasiparticle density of states and bandgap of water. Compared to the SCAN functional, mixing exact exchange mitigates the self-interaction error in the meta-generalized-gradient approximation and further softens liquid water toward the experimental direction. For most of the water properties, the SCAN0 functional shows a systematic improvement over the SCAN functional. However, some important discrepancies remain. The H-bond network predicted by the SCAN0 functional is still slightly overstructured compared to the experimental results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An analytical approach to thermal modeling of Bridgman-type crystal growth. I - One-dimensional analysis

A relatively simple one-dimensional thermal model of the Bridgman growth process has been developed which is applicable to the growth of small diameter samples with conductivities similar to those of metallic alloys. The heat flow in a translating rod is analyzed in a way that is applicable to Biot numbers less than unity. The model accommodates an adiabatic zone, different heat transfer coefficients in the hot and cold zones, and changes in sample material properties associated with phase change. The analysis is applied to several simplified cases. The effect of the rod's motion is studied in a three-zone furnace for a rod sufficiently long that end effects can be neglected; end effects are then investigated for a motionless rod. Finally, the addition of a fourth zone, an independently controlled booster heater between the main heater and the adiabatic zone, is evaluated for its ability to increase the gradient in the sample at the melt interface and to control the position of the interface.

Naumann, R. J.↗

Robust scalable initialization for Bayesian variational inference with multi-modal Laplace approximations

Predictive modeling typically relies on Bayesian model calibration to provide uncertainty quantification. Variational inference utilizing fully independent (“mean-field”) Gaussian distributions are often used as approximate probability density functions. This simplification is attractive since the number of variational parameters grows only linearly with the number of unknown model parameters. However, the resulting diagonal covariance structure and unimodal behavior can be too restrictive to provide useful approximations of intractable Bayesian posteriors that exhibit highly non-Gaussian behavior, including multimodality. High-fidelity surrogate posteriors for these problems can be obtained by considering the family of Gaussian mixtures. Gaussian mixtures are capable of capturing multiple modes and approximating any distribution to an arbitrary degree of accuracy, while maintaining some analytical tractability. Unfortunately, variational inference using Gaussian mixtures with full-covariance structures suffers from a quadratic growth in variational parameters with the number of model parameters. The existence of multiple local minima due to strong nonconvex trends in the loss functions often associated with variational inference present additional complications, These challenges motivate the need for robust initialization procedures to improve the performance and computational scalability of variational inference with mixture models. In this work, we propose a method for constructing an initial Gaussian mixture model approximation that can be used to warm-start the iterative solvers for variational inference. The procedure begins with a global optimization stage in model parameter space. In this step, local gradient-based optimization, globalized through multistart, is used to determine a set of local maxima, which we take to approximate the mixture component centers. Around each mode, a local Gaussian approximation is constructed via the Laplace approximation. Finally, the mixture weights are determined through constrained least squares regression. The robustness and scalability of the proposed methodology is demonstrated through application to an ensemble of synthetic tests using high-dimensional, multimodal probability density functions. Here, the practical aspects of the approach are demonstrated with inversion problems in structural dynamics.

97 MATHEMATICS AND COMPUTING↗

Extended adiabatic blast waves and a model of the soft X-ray background

An analytical approximation is generated which follows the development of an adiabatic spherical blast wave in a homogeneous ambient medium of finite pressure. An analytical approximation is also presented for the electron temperature distribution resulting from coulomb collisional heating. The dynamical, thermal, ionization, and spectral structures are calculated for blast waves of energy E sub 0 = 5 x 10 to the 50th power ergs in a hot low-density interstellar environment. A formula is presented for estimating the luminosity evolution of such explosions. The B and C bands of the soft X-ray background, it is shown, are reproduced by such a model explosion if the ambient density is about .000004 cm, the blast radius is roughly 100 pc, and the solar system is located inside the shocked region. Evolution in a pre-existing cavity with a strong density gradient may, it is suggested, remove both the M band and OVI discrepancies.

Cox, D. P.↗