Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “gradient methods”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Heavy metal separations using strongly paramagnetic column packings in a nonhomogeneous magnetic field

The present invention is directed to a method for the separation of an actinide from another metal. The method comprises the following steps: (a) establishing a non-homogeneous magnetic field across a separation column containing a paramagnetic packing material and (b) providing a fluid containing the actinide and the another metal to the separation column wherein the fluid and the paramagnetic packing material are exposed to the non-homogeneous magnetic field. The non-homogeneous magnetic field is produced by a magnet having a first pole for magnetic interaction with a second pole of the magnet wherein the first pole has a different surface area than the second pole. The non-homogeneous magnetic field has a magnetic field gradient of about 500 lines/cm 2 /cm or more. In addition, the present invention is also directed to a method for the separation of one heavy metal from another heavy metal.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A Kaczmarz-inspired approach to accelerate the optimization of neural network wavefunctions

Neural network wavefunctions optimized using the variational Monte Carlo method have been shown to produce highly accurate results for the electronic structure of atoms and small molecules, but the high cost of optimizing such wavefunctions prevents their application to larger systems. We propose the Subsampled Projected-Increment Natural Gradient Descent (SPRING) optimizer to reduce this bottleneck. SPRING combines ideas from the recently introduced minimum-step stochastic reconfiguration optimizer (MinSR) and the classical randomized Kaczmarz method for solving linear least-squares problems. We demonstrate that SPRING outperforms both MinSR and the popular Kronecker-Factored Approximate Curvature method (KFAC) across a number of small atoms and molecules, given that the learning rates of all methods are optimally tuned. For example, on the oxygen atom, SPRING attains chemical accuracy after forty thousand training iterations, whereas both MinSR and KFAC fail to do so even after one hundred thousand iterations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Analytic Nuclear Gradients for Complete Active Space Linearized Pair-Density Functional Theory

Accurately modeling photochemical reactions is difficult due to the presence of conical intersections and locally avoided crossings, as well as the inherently multiconfigurational character of excited states. As such, one needs a multistate method that incorporates state interaction in order to accurately model the potential energy surface at all nuclear coordinates. The recently developed linearized pair-density functional theory (L-PDFT) is a multistate extension of multiconfiguration PDFT, and it has been shown to be a cost-effective post-MCSCF method (as compared to more traditional and expensive multireference many-body perturbation methods or multireference configuration interaction methods) that can accurately model potential energy surfaces in regions of strong nuclear–electronic coupling in addition to accurately predicting Franck–Condon vertical excitations. Here, in this paper, we report the derivation of analytic gradients for L-PDFT and their implementation in the PySCF-forge software, and we illustrate the utility of these gradients for predicting ground- and excited-state equilibrium geometries and adiabatic excitation energies for formaldehyde, s-trans-butadiene, phenol, and cytosine.

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

The ellipticity parameterization for an NFW profile: An overlooked angular structure in strong lens modeling

Galaxy-scale gravitational lenses are often modeled with two-component mass profiles where one component represents the stellar mass and the second is a Navarro Frenk White (NFW) profile representing the dark matter. Outside of the spherical case, the NFW profile is costly to implement, and so it is approximated via two different methods; ellipticity can be introduced via the lensing potential (NFWp) or via the mass by approximating the NFW profile as a sum of analytical profiles (NFWm). While the NFWp method has been the default for lensing applications, it gives a different prescription of the azimuthal structure, which we show introduces ubiquitous gradients in ellipticity and boxiness in the mass distribution rather than having a constant elliptical shape. Because an unmodeled azimuthal structure has been shown to be able to bias lens model results, we explored the degree to which this azimuthal structure that was introduced can affect the model accuracy. We constructed input profiles using composite models using both the NFWp and NFWm methods and fit these mocks with a power-law elliptical mass distribution (PEMD) model with external shear. As a measure of the accuracy of the recovered lensing potential, we calculated the value of the Hubble parameter H 0 one would determine from the lensing fit. We found that the fits to the NFWp input return H 0 values that are systematically biased by about 3% lower than the NFWm counterparts. We explored whether such an effect is attributable to the mass sheet transformation (MST) by using an MST-independent quantity, ξ 2 . We show that, as expected, the NFWm mocks are degenerate with PEMD through an MST. For the NFWp, an additional bias was found beyond the MST due to the azimuthal structure exterior to the Einstein radius. We recommend modelers use an NFWm prescription in the future, such that the azimuthal structure can be introduced explicitly rather than implicitly.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Enhanced Data Efficiency Using Deep Neural Networks and Gaussian Processes for Aerodynamic Design Optimization

Adjoint-based optimization methods are attractive for aerodynamic shape design primarily due to their computational costs being independent of the dimensionality of the input space and their ability to generate high-fidelity gradients that can then be used in a gradient-based optimizer. This makes them very well suited for high-fidelity simulation based aerodynamic shape optimization of highly parametrized geometries such as aircraft wings. However, the development of adjoint-based solvers involve careful mathematical treatment and their implementation require detailed software development. Furthermore, they can become prohibitively expensive when multiple optimization problems are being solved, each requiring multiple restarts to circumvent local optima. In this work, we propose a machine learning enabled, surrogate-based framework that replaces the expensive adjoint solver, without compromising on predicting predictive accuracy. Specifically, we first train a deep neural network (DNN) from training data generated from evaluating the high-fidelity simulation model on a model-agnostic design of experiments on the geometry shape parameters. The optimum shape may then be computed by using a gradient-based optimizer coupled with the trained DNN. Subsequently, we also perform a gradient-free Bayesian optimization, where the trained DNN is used as the prior mean. We observe that the latter framework (DNN-BO) improves upon the DNN-only based optimization strategy for the same computational cost. Overall, this framework predicts the true optimum with very high accuracy, while requiring far fewer high-fidelity function calls compared to the adjoint-based method. Furthermore, we show that multiple optimization problems can be solved with the same machine learning model with high accuracy, to amortize the offline costs associated with constructing our models. Our methodology finds applications in the early stages of aerospace design. (C) 2021 Published by Elsevier Masson SAS.

Renganathan, S. Ashwin↗
DOE OSTI2020

Coupled Time-Lapse Full-Waveform Inversion for Subsurface Flow Problems Using Intrusive Automatic Differentiation

We describe a novel framework for estimating subsurface properties, such as rock permeability and porosity, from time-lapse observed seismic data by coupling full-waveform inversion (FWI), subsurface flow processes, and rock physics models. For the inverse modeling, we handle the back propagation of gradients by an intrusive automatic differentiation strategy that offers three levels of user control: (1) At the wave physics level, we adopted the discrete adjoint method in order to use our existing high-performance FWI code; (2) at the rock physics level, we used built-in automatic differentiation operators from the TensorFlow backend; (3) at the flow physics level, we implemented customized partial differential equation (PDE) operators for the multiphase flow equations. The three-level coupled inversion strategy strikes a good balance between computational efficiency and programming efforts, and when the gradients are chained together, it constitutes a coupled inverse system. Our numerical experiments demonstrate that the three-level coupled inverse problem is superior in terms of accuracy to a traditional decoupled inversion strategy. Additionally, our method is able to simultaneously invert for parameters in empirical relationships such as the rock physics models. Our proposed inverted model can be used for reservoir performance prediction and reservoir management/optimization purposes.

54 ENVIRONMENTAL SCIENCES↗
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 OSTI2023

Efficient Hydraulic-to-electric Energy Conversion for PCM-based Ocean Thermal Gradient Energy System to Power Uncrewed Underwater Vehicles

Since the 2000, there have been numerous investigations into electrical energy storage types that use ocean thermal-to-electricity conversion techniques. Field tests have demonstrated that this technology is feasible for UUVs with long-term missions of over 3,000 km per deployment. Nonetheless, the energy generated by this method is not sufficient to meet the energy requirements of Unscrewed Underwater Vehicles (UUVs). The objective of this study is to design a highly effective ocean thermal gradient energy system to support long-term missions of a UUV. To enhance the performance of the thermal gradient energy system, we have developed a hybrid type (buoyancy regulation and electrical storage). Phase Change Material (PCM)-based energy harvesting system to facilitate efficient UUV operations. A prototype of the hydraulic-to-electric conversion system was implemented, and experimental findings indicate that the hybrid system successfully delivers hydraulic energy for buoyancy change and improved efficiency of hydraulic-to-electric energy conversion, demonstrating its potential as a promising technology for supporting long-term UUV missions.

Deng, Zhiqun↗
DOE OSTI2021

Near total reflection x-ray photoelectron spectroscopy: quantifying chemistry at solid/liquid and solid/solid interfaces

Near total reflection regime has been widely used in x-ray science, specifically in grazing incidence small angle x-ray scattering and in hard x-ray photoelectron spectroscopy (XPS). In this work, we introduce some practical aspects of using near total reflection (NTR) in ambient pressure XPS and apply this technique to study chemical concentration gradients in a substrate/photoresist system. Experimental data are accompanied by x-ray optical and photoemission simulations to quantitatively probe the photoresist and the interface with the depth accuracy of ~1 nm. Together, our calculations and experiments confirm that NTR XPS is a suitable method to extract information from buried interfaces with highest depth-resolution, which can help address open research questions regarding our understanding of concentration profiles, electrical gradients, and charge transfer phenomena at such interfaces. Furthermore, the presented methodology is especially attractive for solid/liquid interface studies, since it provides all the strengths of a Bragg-reflection standing-wave spectroscopy without the need of an artificial multilayer mirror serving as a standing wave generator, thus dramatically simplifying the sample synthesis.

solid/solid interface↗
DOE OSTI2020

Efficient optimization method for finding minimum energy paths of magnetic transitions

Here, efficient algorithms for the calculation of minimum energy paths of magnetic transitions are implemented within the geodesic nudged elastic band (GNEB) approach. While an objective function is not available for GNEB and a traditional line search can, therefore, not be performed, the use of limited memory Broyden–Fletcher–Goldfarb–Shanno (LBFGS) and conjugate gradient algorithms in conjunction with orthogonal spin optimization (OSO) approach is shown to greatly outperform the previously used velocity projection and dissipative Landau–Lifschitz dynamics optimization methods. The implementation makes use of energy weighted springs for the distribution of the discretization points along the path and this is found to improve performance significantly. The various methods are applied to several test problems using a Heisenberg-type Hamiltonian, extended in some cases to include Dzyaloshinskii–Moriya and exchange interactions beyond nearest neighbours. Minimum energy paths are found for magnetization reversals in a nano-island, collapse of skyrmions in two-dimensional layers and annihilation of a chiral bobber near the surface of a three-dimensional magnet. The LBFGS-OSO method is found to outperform the dynamics based approaches by up to a factor of 8 in some cases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Analytic evaluation of energy first derivatives for spin–orbit coupled-cluster singles and doubles augmented with noniterative triples method: General formulation and an implementation for first-order properties

The paper reports a formulation of analytic energy first derivatives for the coupled-cluster singles and doubles augmented with noniterative triples [CCSD(T)] method with spin–orbit coupling included at the orbital level and an implementation for evaluation of first-order properties. Here, the new development extends the applicability of relativistic coupled cluster methods for heavy-element chemistry and spectroscopy. Benchmark calculations of first-order properties including dipole moments and electric-field gradients using the corresponding exact two-component property integrals are presented for heavy-element containing molecules to demonstrate the usefulness of the present analytic scheme.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Review: Pre-Darcy flows in low-permeability porous media

The widely used Darcy’s law specifies a linear relation between the Darcy velocity of fluid flow and the pressure gradient that drives the flow. However, studies have shown that Darcy velocity can exhibit a nonlinear dependence on the pressure gradient in low-permeability porous media such as clay and shale when the pressure gradient is adequately low. This phenomenon is referred to as low-velocity non-Darcian flow or pre-Darcy flow. This paper provides a comprehensive review of the theories, experimental data, and modeling methods for pre-Darcy flow in low-permeability porous media. The review begins by outlining the fundamental mechanisms underlying pre-Darcy flow that regulate the unique characteristics such as nonlinear dependence of the Darcy velocity on the pressure gradient and its relevance to fluid–rock interactions. The review then proceeds to present a thorough compilation of experimental investigations performed in various low-permeability geomaterials including tight sandstones, shales, and clays. Next, empirical and theoretical models and simulation methods that have been developed to fit and interpret experimental data are reviewed. Finally, the review underscores the challenges encountered in conducting and interpreting pre-Darcy flow experiments and suggests future research directions. By analyzing previous experimental investigations, this review aims to offer a valuable resource for researchers and practitioners seeking to enhance their understanding of fluid dynamics in low-permeability geomaterials. This provides insights into the application of pre-Darcy flow in numerous natural and engineered processes such as shale oil and gas recovery, contaminant transport in low-permeability aquifers, and geological disposal of nuclear waste.

58 GEOSCIENCES↗
DOE OSTI2024

The Influence of Residual Stress on Fatigue Crack Growth Rates in Stainless Steel Processed by Different Additive Manufacturing Methods

The properties and microstructure of Type 304L stainless steel produced by two additive manufacturing (AM) methods—directed energy deposition (DED) and powder bed fusion (PBF)—are evaluated and compared. Localized heating and steep temperature gradients of AM processes lead to significant residual stress and distinctive microstructures, which may be process-specific and influence mechanical behavior. Test data show that materials produced by DED and PDF have small differences in tensile strengths but clear differences in residual stress and microstructural features. Measured fatigue crack growth rates (FCGRs) for cracks propagating parallel to and perpendicular to the build directions differ between the two AM materials. To separate the influences of residual stress and microstructure, K-control test procedures with decreasing and constant stress intensity factor ranges are used to measure FCGRs in the near-threshold regime (crack growth rates ≤ 1 × 10 −8 m/cycle). Residual stress is quantified by the residual stress intensity factor, K res , measured by the online crack compliance method. Correcting the FCGR data for differences in K res brings results for specimens of the two AM materials into agreement with each other and with results for wrought specimens, when the latter are corrected for crack closure. Differences in microstructure and tensile strength have an insignificant influence on FCGRs in these tests.

additive manufacturing↗
DOE OSTI2021

A new procedure for implementing the modified inherent strain method with improved accuracy in predicting both residual stress and deformation for laser powder bed fusion

As a metal additive manufacturing (AM) process, laser powder bed fusion (L-PBF) has been widely used to produce parts with complex geometries. The large thermal gradient caused by the fast, intense, and repeated laser scanning induces significant residual deformation and stress to the as-built parts, which increase manufacturing difficulty and geometrical inaccuracy as a result. The modified inherent strain (MIS) method exploiting multiscale process simulations was developed to simulate residual deformation accurately and efficiently. However, the existing procedure of implementing the MIS method is found to give inaccurate residual stress prediction. Here in this work, a new implementation procedure for the MIS method is proposed to improve the simulation accuracy of residual stress without degrading the residual deformation prediction. The new procedure concerns the application of inherent strains to the part-scale layer-by-layer finite element model to obtain residual stress and deformation field. While the existing implementation of the part-scale MIS model involves only mechanical properties at ambient temperature, the new procedure adds one more solution step employing mechanical properties at an elevated temperature determined from the inherent strain extraction step. Both numerical and experimental studies are conducted to validate the proposed new implementation procedure. It shows that by using the new procedure, the MIS-based simulation can predict both residual stress and deformation of as-built L-PBF metal parts with good accuracy.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A variational mimetic finite difference method for elliptic interface problems on non-matching polytopal meshes with geometric interface inconsistencies

A new variational mimetic finite difference method for elliptic interface problems with perfect and imperfect thermal contacts on non-matching polytopal meshes with geometric interface inconsistencies is developed and analyzed theoretically and numerically. The method is defined on multiple non-matching submeshes with gaps and overlaps along their interfaces. The discrete equations are derived from a minimization problem for the augmented Dirichlet functional. For a perfect thermal contact, the functional uses a modified mimetic gradient with extended stencil which couples unknowns from both sides of an interface, as well as penalty terms to enforce weak continuity of temperature across the interface. The method leads to a symmetric positive definite matrix for any scaling of the penalty terms. For an imperfect thermal contact, the Dirichlet functional is supplemented with a quadratic jump term along the interface related to the interface thermal resistance. We prove that the method conserves the total heat flux across each interface. In conclusion, the obtained results are verified with numerical experiments showing convergence in the discrete L 2 and L ∞ norms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Learning-Based Load Control to Support Resilient Networked Microgrid Operations

Microgrids have proven to be an effective option for increasing the resiliency of critical end-use loads during extreme events. Building on past operational experiences, some microgrid operators are examining the potential to network microgrids to further improve resiliency. However, the frequency deviations experienced on isolated microgrids during transient events, such as switching operations, step increases in load, and loss of generation, are significantly larger than those typically seen on bulk transmission systems. The larger frequency deviations can cause a loss of inverter-connected assets, resulting in a loss of power to critical end-use loads. This paper presents a method of mitigating the impact of transient events by engaging end-use loads using Grid-Friendly Appliance TM (GFA) controllers. An online, i.e., real-time, device-level algorithm is presented, which adjusts individual GFA controller frequency set-points based on the operational characteristics of each end-use load, and on the changing grid dynamic characteristics. The presented method improves the dynamic stability of the networked microgrid operations while minimizing the interruptions to end-use loads. The presented work is validated with dynamic simulations using a modified version of the IEEE 123-node test system with three microgrids, using the GridLAB-D TM simulation environment.

Radhakrishnan, Nikitha↗
DOE OSTI2024

Simple strategy for the simulation of axially symmetric large-area metasurfaces

Metalenses are composed of nanostructures for focusing light and have been widely explored in many exciting applications. However, their expanding dimensions pose simulation challenges. We propose a method to simulate metalenses in a timely manner using vectorial wave and ray tracing models. We sample the metalens’s radial phase gradient and locally approximate the phase profile by a linear phase response. Each sampling point is modeled as a binary blazed grating, employing the chosen nanostructure, to build a transfer function set. The metalens transmission or reflection is then obtained by applying the corresponding transfer function to the incoming field on the regions surrounding each sampling point. Fourier optics is used to calculate the scattered fields under arbitrary illumination for the vectorial wave method, and a Monte Carlo algorithm is used in the ray tracing formalism. We validated our method against finite-difference time domain simulations at 632 nm, and we were able to simulate metalenses larger than 3000 wavelengths in diameter on a personal computer.

Martins, Augusto (ORCID:0000000295546481)↗
DOE OSTI2022

The Transference Number

The performance of rechargeable batteries and other electrochemical systems depends on the rate at which the working ion (often a cation) is transported from one electrode to the other. The cation transference number is an important transport parameter that affects this rate. The purpose of this perspective is to distinguish between approximate and rigorous methods used in the literature to measure the transference number. We emphasize the fact that this parameter is dependent on the reference frame used in the analysis; care must be taken when comparing values obtained from different sources to account for differences in reference frames. We present data obtained from a well‐characterized electrolyte. We compare rigorously determined transference numbers in two reference frames with values obtained by an approximate method. We conclude with a qualitative discussion of the relationship between the transference number and salt concentration gradients that are obtained when current is drawn through a battery.

Gao, Kevin W.↗