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.

SourceResearch connectionFrom yearThrough year

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

At least 397 records · Page 22

All sources NASA NTRS DOE OSTI

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

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↗

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↗

DOE OSTI2024

Sampling in Long-Screened Wells: Issues, Misconceptions, and Solutions

The issues associated with long-screened wells (LSWs) (and open boreholes) at contaminated sites are well documented in the groundwater literature but are still not fully appreciated in practice. As established in seminal and review papers going back over three decades, the interpretation of sampling results from LSWs is challenging in the presence of vertical hydraulic gradients and borehole flow; furthermore, LSWs allow for vertical redistribution of contamination between aquifer layers. Acknowledgment of these issues has led to the development of new technologies and well designs to enable discrete-zone monitoring (DZM), yet LSWs remain common for many reasons, for example, as multipurpose wells, for geophysical logging, and (or) as legacy installations. Despite the literature on LSWs and despite the adoption of DZM at many sites, the use of LSWs persists and the challenges of interpreting sampling results from LSWs remain. In this issue paper, we provide a conceptual overview of the problems posed by LSWs and review existing literature and past work to improve the interpretation of sampling in LSWs. We draw on experience from previous studies at the Hanford Site in eastern WA, USA, and use synthetic examples to illustrate key concepts and challenges for interpretation. A recently published analytical modeling framework is used to develop illustrative synthetic examples and demonstrate a workflow for building scientific intuition to understand issues around interpreting samples from LSWs, which is critical to effective characterization and groundwater remediation at sites with LSWs.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Equatorial Electronic Structure in the Uranyl Ion: Cs 2 UO 2 Cl 4 and Cs 2 UO 2 Br 4

Electric field gradient (EFG) tensors in the equatorial plane of the linear UO 2 2+ ion have been measured by nuclear magnetic resonance (NMR) and nuclear quadrupole resonance (NQR) experiments and computed by relativistic Kohn–Sham methods with and without environment embedding for Cs 2 UO 2 Cl 4 and Cs 2 UO 2 Br 4 . This approach expands the possibilities for probing the electronic structure in uranyl complexes beyond the strongly covalent U–O bonds. The combined analyses find that one of the two largest principal EFG tensor components at the halogen sites points along the U–X bond (X = Cl, Br), and the second is parallel to the UO 2 2+ ion; in Cs 2 UO 2 Cl 4 , the components are nearly equal in magnitude, whereas in Cs 2 UO 2 Br 4 , due to short-range bromide–cesium interactions, the equatorial component is dominant for one pair of Br sites and the axial component is larger for the second pair. Additionally, the directions and relative magnitudes of the field gradient principal axes are found to be sensitive to the σ and π electron donation by the ligands and the model of the environment. Chlorine-35 NQR spectra of 235 U-depleted and 235 U-enriched Cs 2 UO 2 Cl 4 exhibited no uranium-isotope-dependent shift, but the resonance of the depleted sample displayed a 58% broader line width.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

RASER MRI: Magnetic resonance images formed spontaneously exploiting cooperative nonlinear interaction

The spatial resolution of magnetic resonance imaging (MRI) is limited by the width of Lorentzian point spread functions associated with the transverse relaxation rate 1/T 2 * . Here, we show a different contrast mechanism in MRI by establishing RASER (radio-frequency amplification by stimulated emission of radiation) in imaged media. RASER imaging bursts emerge out of noise and without applying radio-frequency pulses when placing spins with sufficient population inversion in a weak magnetic field gradient. Small local differences in initial population inversion density can create stronger image contrast than conventional MRI. This different contrast mechanism is based on the cooperative nonlinear interaction between all slices. On the other hand, the cooperative nonlinear interaction gives rise to imaging artifacts, such as amplitude distortions and side lobes outside of the imaging domain. Contrast mechanism and artifacts are explored experimentally and predicted by simulations on the basis of a proposed RASER MRI theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Self-consistent implementation of locally scaled self-interaction-correction method

Recently proposed local self-interaction correction (LSIC) method is a one-electron self-interaction-correction (SIC) method that uses an iso-orbital indicator to apply the SIC at each point in space by scaling the exchange–correlation and Coulomb energy densities. The LSIC method is exact for the one-electron densities, also recovers the uniform electron gas limit of the uncorrected density functional approximation, and reduces to the well-known Perdew–Zunger SIC (PZSIC) method as a special case. This article presents the self-consistent implementation of the LSIC method using the ratio of Weizsäcker and Kohn–Sham kinetic energy densities as an iso-orbital indicator. The atomic forces as well as the forces on the Fermi-Löwdin orbitals are also implemented for the LSIC energy functional. Results show that LSIC with the simplest local spin density functional predicts atomization energies of the AE6 dataset better than some of the most widely used generalized-gradient-approximation (GGA) functional and barrier heights of the BH6 database better than some of the most widely used hybrid functionals. The LSIC method [a mean absolute error (MAE) of 0.008 Å] predicts bond lengths of a small set of molecules better than the PZSIC-LSDA (MAE 0.042 Å) and LSDA (0.011 Å). Furthermore, this work shows that accurate results can be obtained from the simplest density functional by removing the self-interaction-errors using an appropriately designed SIC method.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Predicting solid state material platforms for quantum technologies

Semiconductor materials provide a compelling platform for quantum technologies (QT). However, identifying promising material hosts among the plethora of candidates is a major challenge. Therefore, we have developed a framework for the automated discovery of semiconductor platforms for QT using material informatics and machine learning methods. Different approaches were implemented to label data for training the supervised machine learning (ML) algorithms logistic regression, decision trees, random forests and gradient boosting. We find that an empirical approach relying exclusively on findings from the literature yields a clear separation between predicted suitable and unsuitable candidates. In contrast to expectations from the literature focusing on band gap and ionic character as important properties for QT compatibility, the ML methods highlight features related to symmetry and crystal structure, including bond length, orientation and radial distribution, as influential when predicting a material as suitable for QT.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Nonlinear Shear Rheology of Entangled Polymer Rings

Steady-state shear viscosity η($\bar{γ}$) of unconcatenated ring polymer melts as a function of the shear rate $\bar{γ}$ is studied by a combination of experiments, simulations, and theory. Experiments using polystyrenes with Z ≈ 5 and Z ≈ 11 entanglements indicate weaker shear thinning for rings compared to linear polymers exhibiting power law scaling of shear viscosity η ~ $\bar{γ}$ –0.56 ± 0.02 , independent of chain length, for Weissenberg numbers up to about 10 2 . Nonequilibrium molecular dynamics simulations using the bead-spring model reveal a similar behavior with η ~ $\bar{γ}$ -0.57 ± 0.08 for 4 ≤ Z ≤ 57. Viscosity decreases with chain length for high $\bar{γ}$. In our experiments, we see the onset of this regime, and in simulations, which we extended to Wi ~ 10 4 , the nonuniversality is fully developed. In addition to a naive scaling theory yielding for the universal regime η ~ $\bar{γ}$ –0.57 , we developed a novel shear slit model explaining many details of observed conformations and dynamics as well as the chain length-dependent behavior of viscosity at large $\bar{γ}$. The signature feature of the model is the presence of two distinct length scales: the size of tension blobs and much larger thickness of a shear slit in which rings are self-consistently confined in the velocity gradient direction and which is dictated by the size of a chain section with relaxation time 1/$\bar{γ}$. These two length scales control the two normal stress differences. In this model, the chain length-dependent onset of nonuniversal behavior is set by tension blobs becoming as small as about one Kuhn segment. This model explains the approximate applicability of the Cox–Merz rule for ring polymers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

CFD simulations of electric motor end ring cooling for improved thermal management

Proper thermal management of an electric motor for vehicle applications extends its operating range. One cooling approach is to impinge Automatic Transmission Fluid (ATF) onto the rotor end ring. Increased ATF coverage correlates to enhanced heat transfer. Computational Fluid Dynamics (CFD) analytical tools provide a mechanism to assess motor thermal management prior to hardware fabrication. The complexity of the fluid flow (e.g., jet atomization, interface tracking, wall impingement) and heat transfer makes these simulations challenging. Computational costs are high when solving these flows on high-speed rotating meshes. Typically, a Volume-of Fluid (VOF) technique (i.e., two-fluid system) is used to resolve ATF dynamics within this rotating framework. Suitable numerical resolution of the relevant physics for thin films under strong inertial forces at high rotor speeds is computationally expensive, further increasing the run times. In this work, a numerical study of rotor-ring cooling by ATF is presented using a patent automated Cartesian cut-cell based method coupled with Automatic Mesh Refinement (AMR). This approach automatically creates the Cartesian mesh on-the-fly and can effectively handle complex rotating geometries by adaptively refining the mesh based on local gradients in the flow field which results in better resolution of the air-ATF interface. A Single non-inertial Reference Frame (SRF) approach is used to account for the rotating geometry and to further improve the overall computational efficiency. Quasi-steady state conditions are targeted in the analysis of the results. Important physics such as ATF jet structure, velocity detail near the air-jet interface, ATF coverage/accumulation on the ring surface, and cooling capacity are presented for a low-resolution Reynolds averaged Navier-Stokes (RANS), high-resolution RANS, and high-resolution Large-Eddy Simulation (LES) models. Computations are scaled over hundreds of cores on a supercomputer to maximize turnaround time. Each numerical approach is shown to capture the general trajectory of the oil jet prior to surface impingement. The high-resolution LES simulation, however, is superior in capturing small scale details and heat transfer between the free jet and surrounding air.

42 ENGINEERING↗

DOE OSTI2025

Introducing GPU Acceleration into the Python-Based Simulations of Chemistry Framework

We introduce the first version of GPU4P Y SCF, a module that provides GPU acceleration of methods in P Y SCF. As a core functionality, this provides a GPU implementation of two-electron repulsion integrals (ERIs) for contracted basis sets comprising up to g functions using the Rys quadrature. As an illustration of how this can accelerate a quantum chemistry workflow, we describe how to use the ERIs efficiently in the integral-direct Hartree–Fock build and nuclear gradient construction. Benchmark calculations show a significant speedup of 2 orders of magnitude with respect to the multithreaded CPU Hartree–Fock code of P Y SCF and the performance comparable to other open-source GPU-accelerated quantum chemical packages, including GAMESS and QUICK, on a single NVIDIA A100 GPU.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗