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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Delocalization error poisons the density-functional many-body expansion

The many-body expansion is a fragment-based approach to large-scale quantum chemistry that partitions a single monolithic calculation into manageable subsystems. This technique is increasingly being used as a basis for fitting classical force fields to electronic structure data, especially for water and aqueous ions, and for machine learning. Here, we show that the many-body expansion based on semilocal density functional theory affords wild oscillations and runaway error accumulation for ion–water interactions, typified by F − (H 2 O) N with N ≳ 15. We attribute these oscillations to self-interaction error in the density-functional approximation. The effect is minor or negligible in small water clusters, explaining why it has not been noticed previously, but grows to catastrophic proportion in clusters that are only moderately larger. This behavior can be counteracted with hybrid functionals but only if the fraction of exact exchange is ≳50%, whereas modern meta-generalized gradient approximations including ωB97X-V, SCAN, and SCAN0 are insufficient to eliminate divergent behavior. Other mitigation strategies including counterpoise correction, density correction (i.e., exchange–correlation functionals evaluated atop Hartree–Fock densities), and dielectric continuum boundary conditions do little to curtail the problematic oscillations. In contrast, energy-based screening to cull unimportant subsystems can successfully forestall divergent behavior. These results suggest that extreme caution is warranted when the many-body expansion is combined with density functional theory.

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

A Two-Dimensional Non-Linear Magnetic Equivalent Circuit Model to Facilitate the Preliminary Design of a Normal Conducting Quadrupole

Normal conducting quadrupoles have been used to focus charged particle beams in synchrotrons, beam transfer lines, medical linacs, etc. for a long time. Optimization techniques based on analytical expressions combined with the use of a numerical field analysis tool exist in the literature for the design of an optimal pole tip shape. However, the initial shape and dimensions of the remaining yoke (including the pole itself and yoke base that act as the return path) are usually less well defined. This article discusses a design methodology for a normal conducting quadrupole based on a two-dimensional magnetic equivalent circuit. Here, this approach considers the geometry of the entire magnet and the non-linear behavior of the yoke material, thereby eliminating the initial iterations. The design outcome of this exercise serves as a competent starting point and can then be used to refine the pole tip, pole taper, yoke size, and add other geometrical features to achieve the required field quality, gradient, coil considerations, etc. by employing a finite element analysis tool. An example design study to demonstrate the proposed methodology is presented.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2015

Variable-Domain Displacement Transfer Functions for Converting Surface Strains into Deflections for Structural Deformed Shape Predictions

Variable-Domain Displacement Transfer Functions were formulated for shape predictions of complex wing structures, for which surface strain-sensing stations must be properly distributed to avoid jointed junctures, and must be increased in the high strain gradient region. Each embedded beam (depth-wise cross section of structure along a surface strain-sensing line) was discretized into small variable domains. Thus, the surface strain distribution can be described with a piecewise linear or a piecewise nonlinear function. Through discretization, the embedded beam curvature equation can be piece-wisely integrated to obtain the Variable-Domain Displacement Transfer Functions (for each embedded beam), which are expressed in terms of geometrical parameters of the embedded beam and the surface strains along the strain-sensing line. By inputting the surface strain data into the Displacement Transfer Functions, slopes and deflections along each embedded beam can be calculated for mapping out overall structural deformed shapes. A long tapered cantilever tubular beam was chosen for shape prediction analysis. The input surface strains were analytically generated from finite-element analysis. The shape prediction accuracies of the Variable- Domain Displacement Transfer Functions were then determined in light of the finite-element generated slopes and deflections, and were fofound to be comparable to the accuracies of the constant-domain Displacement Transfer Functions

displacement transfer functions↗
DOE OSTI2020

Method and device to synthesize boron nitride nanotubes and related nanoparticles

Methods and apparatus for producing chemical nanostructures having multiple elements, such as boron and nitride, e.g. boron nitride nanotubes, are disclosed. The method comprises creating a plasma jet, or plume, such as by an arc discharge. The plasma plume is elongated and has a temperature gradient along its length. It extends along its length into a port connector area having ports for introduction of feed materials. The feed materials include the multiple elements, which are introduced separately as fluids or powders at multiple ports along the length of the plasma plume, said ports entering the plasma plume at different temperatures. The method further comprises modifying a temperature at a distal portion of or immediately downstream of said plasma plume; and collecting said chemical nanostructures after said modifying.

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

Impacts of Latent Energy and Snow Fall Speed on a Wintertime Midlatitude Cyclone

This study aimed to understand the impacts of latent energy as well as snow fall speeds on precipitation properties and synoptic-scale storm characteristics of a wintertime midlatitude cyclone. Simulations of a potent winter storm that impacted the Rocky Mountains and northern Great Plains of the United States in February 2017 were performed using the Weather Research and Forecasting model with latent heating or cooling from individual microphysical processes systematically turned off and fall speeds adjusted by ±50%. Results indicated substantial impacts on the microphysical characteristics of the simulated storm to fall speed, cooling from sublimation, and warming from deposition. The impacts of cooling and warming were manifested as differences in accumulated snowfall. Increased (decreased) fall speeds led to smaller (larger) ice crystals and total mass, resulting in offsetting effects on the precipitation flux, and thus minimal impacts on snowfall and large-scale characteristics of the storm. Warming and cooling associated with deposition and sublimation, respectively, impacted the synoptic-scale dynamics, whereby removing warming from deposition resulted in an increased meridional temperature gradient near the jet stream, thus increasing the jet strength and causing it to be more westerly with less curvature aloft. This in turn limited upper-level divergence, creating a weaker surface low and shifting the precipitation shield southward. The opposite occurred with the removal of latent cooling due to sublimation. This study highlights the potential importance of latent energy associated with ice sublimation and deposition and fall speeds in the larger-scale characteristics of winter storms.

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

Assessing MP2 frozen natural orbitals in relativistic correlated electronic structure calculations

The high computational scaling with the basis set size and the number of correlated electrons is a bottleneck limiting applications of coupled cluster algorithms, in particular for calculations based on two- or four-component relativistic Hamiltonians, which often employ uncontracted basis sets. This problem may be alleviated by replacing canonical Hartree–Fock virtual orbitals by natural orbitals (NOs). Here, in this paper, we describe the implementation of a module for generating NOs for correlated wavefunctions and, in particular, second order Møller–Plesset perturbation frozen natural orbitals (MP2FNOs) as a component of our novel implementation of relativistic coupled cluster theory for massively parallel architectures [Pototschnig et al. J. Chem. Theory Comput. 17, 5509, (2021)]. Our implementation can manipulate complex or quaternion density matrices, thus allowing for the generation of both Kramers-restricted and Kramers-unrestricted MP2FNOs. Furthermore, NOs are re-expressed in the parent atomic orbital (AO) basis, allowing for generating coupled cluster singles and doubles NOs in the AO basis for further analysis. By investigating the truncation errors of MP2FNOs for both the correlation energy and molecular properties—electric field gradients at the nuclei, electric dipole and quadrupole moments for hydrogen halides HX (X = F–Ts), and parity-violating energy differences for H 2 Z 2 (Z = O–Se)—we find MP2FNOs accelerate the convergence of the correlation energy in a roughly uniform manner across the Periodic Table. It is possible to obtain reliable estimates for both energies and the molecular properties considered with virtual molecular orbital spaces truncated to about half the size of the full spaces.

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

Microwave Heating-Induced Temperature Gradients in Liquid–Liquid Biphasic Systems

Microwaves (MWs) can enable the electrification and intensification of chemical manufacturing. They have been applied to various unit separations, such as drying, distillation, and extraction, entailing gas–liquid and solid–liquid systems. However, a limited quantitative understanding of MW-heated liquid–liquid biphasic systems related to extraction exists. This work measures the temporal and spatial temperature difference between an aqueous and an organic phase in batch and continuous microfluidic modes. We demonstrate permanent temperature differences between phases over 35 °C and spatiotemporal periodic and quasiperiodic oscillations modulated by the flow patterns. The temperature differences are primarily driven by the faster absorption rate of MW irradiation by the aqueous phase versus the slower heat transfer from the aqueous phase to the organic phase. These are amplified by low specific interfacial area and modifications of the electromagnetic field. We employ a multiphysics model to predict the temperature difference in a batch system. The model is in good agreement with the experiments. We demonstrate a strong effect of input power, dielectric properties of organic solvents, the volume of solvents, and the volume ratio between phases on the temperature difference. A simple analytical model describes the temperature difference and provides design principles. Furthermore, the combined approach offers new insights into the design and optimization of the MW-heated biphasic systems.

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

Topological Guided Detection of Extreme Wind Phenomena: Implications for Wind Energy

Extreme wind phenomena play a crucial role in the efficient operation of wind farms for renewable energy generation. However, existing detection methods are computationally expensive, limited to specific coordinate. In real-world scenarios, understanding the occurrence of these phenomena over a large area is essential. Therefore, there is a significant demand for a fast and accurate approach to forecast such events. In this paper, we propose a novel method for detecting wind phenomena using topological analysis, leveraging the gradient of wind speed or critical points in a topological framework. By extracting topological features from the wind speed profile within a defined region, we employ topological distance to identify extreme wind phenomena. Our results demonstrate the effectiveness of utilizing topological features derived from regional wind speed profiles. We validate our approach using high-resolution simulations with the Weather Research and Forecasting model (WRF) over a month in the US East Coast.

MATHEMATICS AND COMPUTING,WIND ENERGY↗
DOE OSTI2023

Topological Guided Detection of Extreme Wind Phenomena: Implications for Wind Energy: Preprint

Extreme wind phenomena play a crucial role in the efficient operation of wind farms for renewable energy generation. However, existing detection methods are computationally expensive, limited to specific coordinate. In real-world scenarios, understanding the occurrence of these phenomena over a large area is essential. Therefore, there is a significant demand for a fast and accurate approach to forecast such events. In this paper, we propose a novel method for detecting wind phenomena using topological analysis, leveraging the gradient of wind speed or critical points in a topological framework. By extracting topological features from the wind speed profile within a defined region, we employ topological distance to identify extreme wind phenomena. Our results demonstrate the effectiveness of utilizing topological features derived from regional wind speed profiles. We validate our approach using high-resolution simulations with the Weather Research and Forecasting model (WRF) over a month in the US East Coast.

MATHEMATICS AND COMPUTING,WIND ENERGY↗
DOE OSTI2025

Feature learning and generalization in deep networks with orthogonal weights

Fully-connected deep neural networks with weights initialized from independent Gaussian distributions can be tuned to criticality, which prevents the exponential growth or decay of signals propagating through the network. However, such networks still exhibit fluctuations that grow linearly with the depth of the network, which may impair the training of networks with width comparable to depth. We show analytically that rectangular networks with tanh activations and weights initialized from the ensemble of orthogonal matrices have corresponding preactivation fluctuations which are independent of depth, to leading order in inverse width. Moreover, we demonstrate numerically that, at initialization, all correlators involving the neural tangent kernel (NTK) and its descendants at leading order in inverse width—which govern the evolution of observables during training—saturate at a depth of ~20, rather than growing without bound as in the case of Gaussian initializations. We speculate that this structure preserves finite-width feature learning while reducing overall noise, thus improving both generalization and training speed in deep networks with depth comparable to width. We provide some experimental justification by relating empirical measurements of the NTK to the superior performance of deep non-linear orthogonal networks trained under full-batch gradient descent on the MNIST and CIFAR-10 classification tasks.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Investigation of Scanning Droplet Cell Technology for Electrochemical Deposition of Custom Three-Dimensional Alloys

The goal of this project was to assess the feasibility of a new method of electroplating, which we have termed “electroprinting”, for fabricating millimeter to centimeter metal parts with full density, and eventually, with bespoke 3D internal density patterns. Alloys, gradients and varying density is beyond the scope of this work, and efforts were focused on developing the technology to be capable of printing solid parts, beyond the lines and columns previously reported in the literature. Enabling this technology would expand the design space possible for the WPD program, allowing for smooth and complex density gradients in parts rather than discrete density steps between multilayers. We successfully designed and built an electroprinting apparatus, capable of printing copper in customizable 1-D patterns, which can be printed in stacked layers to form 3D parts. We successfully characterized the flat printed patterns, however, have encountered difficulty in characterizing multilayer prints. We have partially addressed the feasibility question, by developing the method for electroprinting 3D parts, however, some questions about the internal porosity and density of these parts still remain.

36 MATERIALS SCIENCE↗
NASA NTRS1996

A Thermostructural Analysis of a Diboride Composite Leading Edge

In an effort to support the design of zirconium diboride composite leading edges for hypersonic vehicles, a finite element model (FEM) of a prototype leading edge was created and finite element analysis (FEA) was employed to assess its thermal and structural response to aerothermal boundary conditions. Unidirectional material properties for the structural components of the leading edge, a continuous fiber reinforced diboride composite, were computed with COSTAR. These properties agree well with those experimentally measured. To verify the analytical approach taken with COSMOS/M, an independent FEA of one of the leading edge assembly components was also done with COSTAR. Good agreement was obtained between the two codes. Both showed that a unidirectional lay-up had the best margin of safety for a simple loading case. Both located the maximum stress in the same region and ply. The magnitudes agreed within 4 percent. Trajectory based aerothermal heating was then applied to the leading edge assembly FEM created with COSMOS/M to determine steady state temperature response, displacement, stresses, and contact forces due to thermal expansion and thermal strains. Results show that the leading edge stagnation line temperature reached 4700 F. The maximum computed failure index for the laminated composite components peaks at 4.2, and is located at the bolt flange in layer 2 of the side bracket. The temperature gradient in the tip causes a compressive stress of 279 ksi along its width and substantial tensile stresses within its depth.

Kowalski, Tom↗
NASA NTRS1998

Numerical Study of Interaction of a Vortical Density Inhomogeneity with Shock and Expansion Waves

We studied the interaction of a vortical density inhomogeneity (VDI) with shock and expansion waves. We call the VDI the region of concentrated vorticity (vortex) with a density different from that of ambiance. Non-parallel directions of the density gradient normal to the VDI surface and the pressure gradient across a shock wave results in an additional vorticity. The roll-up of the initial round VDI towards a non-symmetrical shape is studied numerically. Numerical modeling of this interaction is performed by a 2-D Euler code. The use of an adaptive unstructured numerical grid makes it possible to obtain high accuracy and capture regions of induced vorticity with a moderate overall number of mesh points. For the validation of the code, the computational results are compared with available experimental results and good agreement is obtained. The interaction of the VDI with a propagating shock wave is studied for a range of initial and induced circulations and obtained flow patterns are presented. The splitting of the VDI develops into the formation of a non-symmetrical vortex pair and not in a set of vortices. A method for the analytical computation of an overall induced circulation Gamma(sub 1) as a result of the interaction of a moving VDI with a number of waves is proposed. Simplified, approximated, expressions for Gamma(sub 1) are derived and their accuracy is discussed. The splitting of the VDI passing through the Prandtl-Meyer expansion wave is studied numerically. The obtained VDI patterns are compared to those for the interaction of the VDI with a propagating shock wave for the same values of initial and induced circulations. These patterns have similar shapes for corresponding time moments.

Povitsky, A.↗
NASA NTRS2005

Computation of Sound Propagation by Boundary Element Method

This report documents the development of a Boundary Element Method (BEM) code for the computation of sound propagation in uniform mean flows. The basic formulation and implementation follow the standard BEM methodology; the convective wave equation and the boundary conditions on the surfaces of the bodies in the flow are formulated into an integral equation and the method of collocation is used to discretize this equation into a matrix equation to be solved numerically. New features discussed here include the formulation of the additional terms due to the effects of the mean flow and the treatment of the numerical singularities in the implementation by the method of collocation. The effects of mean flows introduce terms in the integral equation that contain the gradients of the unknown, which is undesirable if the gradients are treated as additional unknowns, greatly increasing the sizes of the matrix equation, or if numerical differentiation is used to approximate the gradients, introducing numerical error in the computation. It is shown that these terms can be reformulated in terms of the unknown itself, making the integral equation very similar to the case without mean flows and simple for numerical implementation. To avoid asymptotic analysis in the treatment of numerical singularities in the method of collocation, as is conventionally done, we perform the surface integrations in the integral equation by using sub-triangles so that the field point never coincide with the evaluation points on the surfaces. This simplifies the formulation and greatly facilitates the implementation. To validate the method and the code, three canonic problems are studied. They are respectively the sound scattering by a sphere, the sound reflection by a plate in uniform mean flows and the sound propagation over a hump of irregular shape in uniform flows. The first two have analytical solutions and the third is solved by the method of Computational Aeroacoustics (CAA), all of which are used to compare the BEM solutions. The comparisons show very good agreements and validate the accuracy of the BEM approach implemented here.

Guo, Yueping↗
DOE OSTI2020

Analysis of chemical stress and the propensity for cracking during the vertical Bridgman growth of BaBrCl:Eu

Computational models are employed to analyze residual chemical stresses arising from compositional variations in europium-doped BaBrCl crystals grown by the vertical Bridgman method. We find that significant chemical stress is produced by radial segregation of Eu in this system. In particular, the distribution of normal stresses is set by the radial concentration gradient, whose changing sign produces surface states in tension or compression. Additionally, crack opening from surface flaws will be promoted or suppressed by tensile or compressive surface stresses, respectively. Thus, crystal growth processing strategies that change the radial dopant concentration gradients are posited to affect the propensity for cracking. For this system, surface stresses are changed from states of tension to compression when the growth rate is increased, thus improving the chances to avoid cracking—a strategy that defies classical wisdom that dictates slower growth to improve outcomes. Similar strategies affecting segregation may prove beneficial to tailor chemical stress fields to reduce cracking in other crystal growth systems.

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

Isotopic Substitution Reveals the Importance of Aluminate Diffusion Dynamics in Gibbsite (Al(OH) 3 ) Crystallization from Alkaline Aqueous Solution

Understanding molecular-scale factors governing the precipitation of aluminum hydroxides, such as gibbsite, under alkaline conditions is important for the formation of laterite deposits, as well as aluminum processing. However, mechanisms enabling tetrahedral aluminate ions to assemble into octahedral sites of the gibbsite lattice remain unclear. Formation of oligomeric complexes has been hypothesized as a critical intermediate step. Here, we report a study of gibbsite solubility in highly alkaline solutions using deuterium substitution to probe equilibrium and kinetic factors that could affect oligomeric intermediate formation, including the reactivity and the diffusivity of aluminate ions. When substituting sodium hydroxide with sodium deuteroxide, solution analysis shows a nearly 40% and 50% increase in gibbsite solubility in 2.4 and 3.3 mol·kg –1 total sodium solutions, respectively. Raman spectroscopy indicated that monomeric aluminate ions are the predominant species in solution irrespective of deuteration. However, both 27 Al and 23 Na nuclear magnetic resonance (NMR) spectroscopy revealed significant differences in chemical shifts in deuterated solutions, and analysis of 1 H, 23 Na, and 27 Al diffusion coefficients with pulsed-field gradient, stimulated echo NMR spectroscopy shows a decrease in ion diffusivity with increasing total deuterium in solution, consistent with the notion of increasing strength in the hydrogen/deuterium bonding network. In conclusion, because the relative change in 1 H diffusion coefficients are commensurate with the difference in apparent solubility constants, there is evidence for an oligomeric intermediate whose steady-state concentration is maintained by the collision frequency of aluminate ions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1988

Barotropic stability of realistic stratospheric jets

The stability of realistic jets is examined in a linearized barotropic model on a sphere. Approximately non-dispersive modes associated with a region of negative basic state absolute vorticity gradient on the poleward side of the jet are examined in detail. As in previous studies, broader jets and those which peak at higher latitudes produce poleward modes that are less dispersive. Jet profiles derived from observational data at 10, 5, and 2 mb for three Southern Hemisphere winter months are used in the model, and the results are compared with quasi-nondispersive features (QNDF) which have been observed in satellite data in the Southern Hemisphere winter stratosphere. Characteristics of the barotropically unstable modes compare remarkably well with those of the observed modes. The barotropic model results for a month in which these features are not observed indicate the presence of equatorward modes at wavenumbers 3 and 4 which grow considerably faster than the quasi-nondispersive poleward modes. The appearance of westward moving modes in the summer hemisphere during June and in analytical profiles with a realistic global structure is also noted.

Manney, Gloria L.↗
NASA NTRS1992

Adjustment of footprint correction for airborne flux mapping over the FIFE site

A meaningful interpretation of airborne flux estimates over nonuniform terrain must consider local advection from nonhomogeneous source distributions (footprint correction). An empirical procedure for footprint correction of airborne flux data obtained in the First ISLSCP Field Experiment (FIFE) 1989 is presented, based on simple analytical solutions of the diffusion equation which have been adjusted to approximate more realistic solutions. Optimized estimates of surface flux 'maps' derived from airborne observations are then constructed for maximum correlation between flux estimates and independently observed surface characteristics. This paper also includes an analysis of resolution for the given procedure, in terms of amplitude recovery of hypothetical square-wave distributions of surface flux density, and discusses its implications for the interpretation of the optimized estimates of FIFE surface flux maps. Results show that corrected spatial sink distributions are obtained by the current empirical procedure at the cost of considerable reduction in amplitude recovery of small-scale variations in surface flux density. The high spatial correlation observed between corrected flux maps and surface characteristics such as greenness or surface temperature excess is ascribed to the approximately correct positioning of the main flux gradients across the site, which can be expected to have been reproduced with amplitude recovery equal to or greater than 60 percent.

Schuepp, P. H.↗