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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Western Sector In-Situ Chemical Oxidation Project: Additional Injection Pilot Testing Supplemental Results (U)

Trichloroethylene (TCE) and tetrachloroethylene (PCE) are persistent at elevated concentrations in the Western Sector of the M-Area Hazardous Waste Management Facility (HWMF) groundwater plume. Migration of these chlorinated ethenes as Dense Non-Aqueous Phase Liquids (DNAPLs) occurred from the M-Area Settling Basin along geologic strata resulting in a high concentration groundwater plume in Western Sector. Pump-and-treat has been active for over 30 years at the M-Area HWMF to hydraulically contain the primary sources of the plume in addition to other remediation strategies throughout the contaminated zone, like soil vapor extraction, dynamic underground stripping, groundwater recirculation wells, and bioremediation. In-situ chemical oxidation (ISCO) has been identified as a viable treatment method for removal of these volatile organic compounds (VOCs) due to the complex lithology and high concentrations in this part of the aquifer. In addition to ISCO, a new recovery well (i.e., RWM018) was also installed in Western Sector to contain the high concentration VOC plume. The ISCO projects have utilized the gradient induced by RWM018 to help plan the injection and monitoring locations. The combination of ISCO and the additional pump and treat in Western Sector with the installation of RWM018 will help to control the plume and maximize remediation efforts.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Accumulation of Antifreeze Proteins on Ice Is Determined by Adsorption

Antifreeze proteins (AFPs) enable the survival of organisms in freezing environments by binding to ice crystals and suppressing their growth. The rate at which AFPs accumulate on ice is determined by an interplay between their diffusion from the aqueous solution and their adsorption to the ice crystal. To interrogate the relative importance of these two processes, we combined fluorescence experiments with reaction–diffusion modeling. We found that as AFPs accumulate on ice, their concentration in the aqueous solution does not develop a gradient but remains equal to its bulk concentration throughout our experiments. Furthermore, these findings led us to conclude that AFP accumulation on ice crystals, which are smaller than 100 μm in radius, is not limited by the diffusion of AFPs, but by the kinetics of AFP adsorption.

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

Harnessing deep reinforcement learning to construct time-dependent optimal fields for quantum control dynamics

Here, we present an efficient deep reinforcement learning (DRL) approach to automatically construct time-dependent optimal control fields that enable desired transitions in dynamical chemical systems. Our DRL approach gives impressive performance in constructing optimal control fields, even for cases that are difficult to converge with existing gradient-based approaches. We provide a detailed description of the algorithms and hyperparameters as well as performance metrics for our DRL-based approach. Our results demonstrate that DRL can be employed as an effective artificial intelligence approach to efficiently and autonomously design control fields in quantum dynamical chemical systems.

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

PyFLOSIC: Python-based Fermi–Löwdin orbital self-interaction correction

We present pyflosic, an open-source, general-purpose python implementation of the Fermi–Löwdin orbital self-interaction correction (FLO-SIC), which is based on the python simulation of chemistry framework (pyscf) electronic structure and quantum chemistry code. Thanks to pyscf, pyflosic can be used with any kind of Gaussian-type basis set, various kinds of radial and angular quadrature grids, and all exchange-correlation functionals within the local density approximation, generalized-gradient approximation (GGA), and meta-GGA provided in the libxc and xcfun libraries. A central aspect of FLO-SIC is the Fermi-orbital descriptors, which are used to estimate the self-interaction correction. Importantly, they can be initialized automatically within pyflosic; they can also be optimized within pyflosic with an interface to the atomic simulation environment, a python library that provides a variety of powerful gradient-based algorithms for geometry optimization. Although pyflosic has already facilitated applications of FLO-SIC to chemical studies, it offers an excellent starting point for further developments in FLO-SIC approaches, thanks to its use of a high-level programming language and pronounced modularity.

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

Permutation invariant polynomial neural network based diabatic ansatz for the ( E + A ) × ( e + a ) Jahn–Teller and Pseudo-Jahn–Teller systems

In this work, the permutation invariant polynomial neural network (PIP-NN) approach is employed to construct a quasi-diabatic Hamiltonian for system with non-Abelian symmetries. It provides a flexible and compact NN-based diabatic ansatz from the related approach of Williams, Eisfeld, and co-workers. The example of $H^3_+$ is studied, which is an (E + A) × (e + a) Jahn–Teller and Pseudo-Jahn–Teller system. The PIP-NN diabatic ansatz is based on the symmetric polynomial expansion of Viel and Eisfeld, the coefficients of which are expressed with neural network functions that take permutation-invariant polynomials as input. This PIP-NN-based diabatic ansatz not only preserves the correct symmetry but also provides functional flexibility to accurately reproduce ab initio electronic structure data, thus resulting in excellent fits. Further, the adiabatic energies, energy gradients, and derivative couplings are well reproduced. A good description of the local topology of the conical intersection seam is also achieved. Therefore, this diabatic ansatz completes the PIP-NN based representation of DPEM with correct symmetries and will enable us to diabatize even more complicated systems with complex symmetries.

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

Non-local large-strain FFT-based formulation and its application to interface-dominated plasticity of nano-metallic laminates

Here, this paper presents a novel formulation and its robust numerical implementation of strain-gradient (SG) crystal plasticity within a large-strain (LS) elasto-viscoplastic (EVP) fast Fourier transform (FFT)-based micromechanical model. The resulting non-local SG-LS-EVPFFT formulation is used to model and understand the process of kink band formation during layer-parallel compression of nano-metallic laminates (NMLs). NMLs are layered composites with nanoscale thicknesses, thus requiring consideration of the interaction between dislocations and interfaces within the micromechanical model. The length-scale parameter of the SG model is calibrated by simulating a double pile-up and comparing predictions to analytical solution. This required new expressions for the defect energy, resulting in more accurate double pile-up predictions. The calibrated SG-LS-EVPFFT model is then used to simulate layer-parallel compression of copper-niobium NML. Formation of kink bands is predicted, and the model is used to rationalize the microscopic mechanisms enabling the formation process. It is found that accumulation of dislocations at interfaces leads to activation of layer-parallel slip, which in turn leads to kink band formation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Enriching Cations Using Electric Field Gradients Generated by Bipolar Electrodes in the Absence of Buffer

We discuss that ionic current gradients and corresponding electric field gradients formed at bipolar electrodes are useful for enriching cations up to ~350-fold at specific locations in microfluidic devices. The enrichment mechanism is determined by collecting current data, measuring solution conductivity, and by performing pH indicator experiments. These data reveal the dynamic nature of the enrichment mechanism. Specifically, they show the complex relationship between the ionic current gradients formed at the bipolar electrode and the extent of cation enrichment. Finally, we demonstrate cation enrichment is possible in buffer-free solutions containing up to 50.0 mM KCl. This finding is significant because it broadens the range of solutions in which membrane-free, electrochemical separations using bipolar electrodes can be performed.

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

One-pot ternary sequential reactions for photopatterned gradient multimaterials

Seamless multimaterial construction is a common motif in animal physiology. Such continuous mechanical gradients remain challenging to reproduce in engineered systems, as current resin chemistries typically result in a single fixed set of properties. As an alternative to single-property materials, we introduce a thiol-ene-epoxy-based photothermal reaction scheme that produces multimaterials by altering the polymer microstructure within a single resin. In this system, the photodosage during the first stage of processing dictates the extent of conversion for each subsequent reaction. As a result, our photochemistry can exhibit a diverse range of soft (Young’s modulus, E ~ 400 kPa; elongation, dL/L0 ~ 300%) and stiff (E ~ 1.6 GPa; dL/L 0 ~ 3%) mechanical properties. Furthermore, we pattern photostable and mechanically robust modulus gradients (d[E r, stiff /E r , soft] /dx > 1,000 mm -1 ) that exceed those found in squid beaks and human knee entheses. We demonstrate the ability to build intricate multimaterial architectures including a soft, wearable braille display.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Electronic Structure Calculations of the Fundamental Interactions in Solvent Extraction Desalination

Solvent extraction desalination (SED) is a promising membrane-free desalination technology, but new solvents are still needed that can improve the molecular-level water selectivity, while minimizing solvent contamination in the aqueous phase. Improvements in solvent design for SED require a fundamental understanding of the molecular interactions governing solvation behavior. Herein, quantum chemical calculations are used to delineate different fundamental interactions in water-solvent systems involving 38 different solvents, with a particular focus on the role of the electrostatic potential characteristics. Furthermore, these interactions are rigorously analyzed using the independent gradient model based on Hirshfeld partition (IGMH), Symmetry-Adapted Perturbation Theory (SAPT), as well as the energy decomposition analysis based on force field (EDA-FF) to compare against a standard force field model. The solvation free energies and partition coefficients are also calculated to estimate water extraction performance in a bulk solvent system.

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

Analytically differentiable metrics for phase stability

Here, in this work, a long-established but sparsely documented method of obtaining semi-analytic derivatives of thermodynamic properties with respect to equilibrium conditions is briefly reviewed and rigorously derived. This procedure is then leveraged to construct general forms of derivatives of the residual driving force, a metric for measuring phase stability used in CALPHAD model optimization, with respect to overall system and individual phase compositions. Applied examples – calculating heat capacity in the Al-Fe system, thermodynamic factors in the Nb-V-W system, and residual driving force derivatives in the Ni-Ti system – demonstrate the versatility, accuracy, and extensibility of this method. Using the developed method, residual driving force gradients can be applied directly in CALPHAD model optimizers, as well as in materials design frameworks, to identify regions of phase stability with an efficient, gradient-based approach.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The role of gas flow distributions on CO 2 mineralization within monolithic cemented composites: coupled CFD-factorial design approach

The carbonation kinetics of monolithic cementing composites are strongly affected by gas transport which is, in turn, influenced by microstructural resistances and the presence of liquid water within pore networks. The non-uniform gas flow distribution within the CO 2 mineralization reactor can impart mass transfer resistance in the monolith microstructure, which affects the uptake of CO 2 (“carbonation”) of the cementing composites. This paper demonstrates how the gas spatial distribution (velocity and flow rate; quantified by CFD analysis) and processing conditions (temperature, relative humidity, and flow rate; quantified by factorial design) affect drying and carbonation, and in turn, the engineering properties of a representative ‘monolithic’ carbonate-cemented concrete component (i.e., herein concrete masonry unit: CMUs, also known as concrete block). It is shown that the gas flow distribution affects drying front penetration and results in moisture and carbonation gradients within the monolith. Particularly, variations in drying kinetics caused by non-uniformity of the contacting gas velocity impose gradients in moisture saturation, which results in increasing microstructural resistance to CO 2 transport. The resultant non-uniform carbonate-mineral formation (i.e., carbonate cementation), if not controlled, can produce gradients in mechanical properties and may alter failure patterns upon loading. Finally, these insights inform the optimal design of gas flow distribution systems and processing conditions within CO 2 mineralization reactors for the manufacturing of low-CO 2 concrete components using CO 2 -dilute industrial flue gas streams.

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

Self-consistent theory of the Darrieus–Landau and Rayleigh–Taylor instabilities with self-generated magnetic fields

The Rayleigh–Taylor (RT) and Darrieus–Landau (DL) instabilities are studied in an inertial confinement fusion context within the framework of small critical-to-shell density ratio D R and weak acceleration regime, i.e., large Froude number Fr. The quasi-isobaric analysis in Sanz et al. is completed with the inclusion of non-isobaric and self-generated magnetic-field effects. The analysis is restricted to perturbation wavelengths k~1 larger than the conduction length scale at the ablation front, yet its validity ranges from wavelengths shorter and larger than the conduction layer width (distance between the ablation front and the critical surface). The use of a sharp boundary model leads to a single analytical expression of the dispersion relation encompassing both instabilities. The two new effects come into play by modifying the perturbed mass and momentum fluxes at the ablation front. The momentum flux (perturbed pressure at the spike) is the predominant stabilizing mechanism in the RT instability (overpressure) and the driving mechanism in the DL instability (underpressure). The non-isobaric effects notably modify the scaling laws in the DL limit, leading to an underpressure scaling as ~k –11/15 rather than ~k –2/5 obtained in the quasi-isobaric model. The magnetic fields are generated due to misalignment between pressure and density gradients (Biermann battery effect). They affect the hydrodynamics by bending the heat flux lines. Within the framework of this paper, they enhance ablation, resulting in a stabilizing effect that peaks for perturbation wavelengths comparable to the conduction layer width. Here, the combination of parameters D R Fr 2/3 defines the region of predominance of each instability in the dispersion relation. It is proven that the DL region falls outside of the parameter range in inertial confinement fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Density Matrix Implementation of the Fermi–Löwdin Orbital Self-Interaction Correction Method

The Fermi–Löwdin orbital self-interaction correction (FLOSIC) method effectively provides a transformation from canonical orbitals to localized Fermi–Löwdin orbitals which are used to remove the self-interaction error in the Perdew–Zunger (PZ) framework. This transformation is solely determined by a set of points in space, called Fermi–Löwdin descriptors (FODs), and the occupied canonical orbitals or the density matrix. In this work, we provide a detailed workflow for the implementation of the FLOSIC method for removal of self-interaction error in DFT calculations in an orbital-by-orbital basis that takes advantage of the unitary invariant nature of the FLOSIC method. In this way, it is possible to cast the self-consistent energy minimization at fixed FODs in the same manner than standard Kohn–Sham with one additional term in the Kohn–Sham Hamiltonian that introduces the PZ self-interaction correction. Each energy minimization iteration is divided in two substeps, one for the density matrix and one for the FODs. Expressions for the effective Kohn–Sham matrix and FOD gradients are provided such that its implementation is suitable for most electronic structure codes. Here, we analyze the convergence characteristics of the algorithm and present applications for the evaluation of NMR shielding constants and real-time time-dependent DFT simulations based on the Liouville–von Neumann equation to calculate excitation energies.

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

Post-Disaster Microgrid Formation for Enhanced Distribution System Resilience

This paper proposes a deep reinforcement learning (DRL) based approach for post-disaster critical load restoration in active distribution systems to form microgrids through network reconfiguration to minimize critical load curtailments. Distribution networks are represented as graph networks, and optimal network configurations with microgrids are obtained by searching for the optimal spanning forest. The constraints to the research question being explored are the radial topology and power balance. Unlike existing analytical and population-based approaches, which necessitate the repetition of entire analyses and computation for each outage scenario to find the optimal spanning forest, the proposed approach, once properly trained, can quickly determine the optimal, or near-optimal, spanning forest even when outage scenarios change. When multiple lines fail in the system, the proposed approach forms microgrids with distributed energy resources in active distribution systems to reduce critical load curtailment. The proposed DRL-based model learns the action-value function using the REINFORCE algorithm, which is a model-free reinforcement learning technique based on stochastic policy gradients. A case study was conducted on a 33-node distribution test system, demonstrating the effectiveness of the proposed approach for post-disaster critical load restoration.

active distribution systems↗
DOE OSTI2023

Perfect Crystals: microgravity capillary counterdiffusion crystallization of human manganese superoxide dismutase for neutron crystallography

The NASA mission Perfect Crystals used the microgravity environment on the International Space Station (ISS) to grow crystals of human manganese superoxide dismutase (MnSOD)—an oxidoreductase critical for mitochondrial vitality and human health. The mission’s overarching aim is to perform neutron protein crystallography (NPC) on MnSOD to directly visualize proton positions and derive a chemical understanding of the concerted proton electron transfers performed by the enzyme. Large crystals that are perfect enough to diffract neutrons to sufficient resolution are essential for NPC. This combination, large and perfect, is hard to achieve on Earth due to gravity-induced convective mixing. Capillary counterdiffusion methods were developed that provided a gradient of conditions for crystal growth along with a built-in time delay that prevented premature crystallization before stowage on the ISS. Here, we report a highly successful and versatile crystallization system to grow a plethora of crystals for high-resolution NPC.

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

Structure optimization with stochastic density functional theory

Linear-scaling techniques for Kohn–Sham density functional theory are essential to describe the ground state properties of extended systems. Still, these techniques often rely on the localization of the density matrix or accurate embedding approaches, limiting their applicability. In contrast, stochastic density functional theory (sDFT) achieves linear- and sub-linear scaling by statistically sampling the ground state density without relying on embedding or imposing localization. In return, ground state observables, such as the forces on the nuclei, fluctuate in sDFT, making optimizing the nuclear structure a highly non-trivial problem. In this work, we combine the most recent noise-reduction schemes for sDFT with stochastic optimization algorithms to perform structure optimization within sDFT. We compare the performance of the stochastic gradient descent approach and its variations (stochastic gradient descent with momentum) with stochastic optimization techniques that rely on the Hessian, such as the stochastic Broyden–Fletcher–Goldfarb–Shanno algorithm. In conclusion, we further provide a detailed assessment of the computational efficiency and its dependence on the optimization parameters of each method for determining the ground state structure of bulk silicon with varying supercell dimensions.

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

Pressure Modification or Barrier Issues during Polymer Flooding Enhanced Oil Recovery

Effective oil displacement from a reservoir requires adequate and properly directed pressure gradients in areas of high oil saturation. If the polymer bank is too large or too viscous during a polymer flood, the pressure drops from the injection well to the polymer front may act as a pressure modification or barrier by usurping most of the downstream driving force for oil displacement. Polymer injection pressures must be limited. The maximum allowable injection pressure is commonly constrained by caprock integrity, injection equipment, and/or regulations, even though fractures can be beneficial to polymer injectivity. This paper examines when the pressure-barrier concept limits the size and viscosity of the polymer bank during a polymer flood. Analytical and numerical methods are used to address this issue. We examine the relevance of the pressure modification concept for a wide variety of circumstances, including oil viscosities ranging from 10 cp to 1650 cp, vertical wells versus horizontal wells, single versus multiple layered reservoirs, permeability contrast, and crossflow between layers. We also examine the relation between the pressure-barrier concept and fractures and fracture extension during polymer injection. We demonstrate that in reservoirs with single layers, the pressure-barrier concept only limits the optimum viscosity of the injected polymer if the mobility of the polymer bank is less than the mobility of the displaced oil bank. The same is true for multizoned reservoirs with no crossflow between layers. Thus, for these cases, the optimum polymer viscosity is likely to be dictated by the mobility of the oil bank, unless other factors intervene. For multizoned reservoirs with free crossflow between layers, the situation is different. A compromise must be reached between injected polymer viscosity and the efficiency of oil recovery. This work is particularly relevant to viscous oil reservoirs where polymer viscosities are substantially lower than the oil viscosity.

58 GEOSCIENCES↗
DOE OSTI2026

Historic climate, cosmogenic 10Be, denudation-rate, and geospatial datasets from the Pikes Peak region, Colorado, USA

This data package contains geographic information system (GIS) layers and tabular datasets associated with the study of elevation-dependent denudation rates on Pikes Peak in the Front Range of the Rocky Mountains, Colorado, USA. The package includes GIS layers used to produce the study-area map, including sample locations, sample watershed boundaries, the Pikes Peak batholith, Pleistocene glacier extent, weather station locations, and elevation and hillshade rasters, together with comma-separated value (CSV) tables and matching CSV data dictionaries. These mapped layers provide the geographic framework for interpreting denudation patterns across the Pikes Peak region and for relating sample locations to watershed geometry, bedrock setting, glacial history, and nearby climate stations. The first group of tables reports climate and geospatial context for the study area. These files include station-based temperature and precipitation data used to characterize elevational gradients in mean annual climate and monthly climate seasonality, sample locations, denudation-rate and topographic metrics, fixed frost-cracking model parameters, frost-cracking intensity and precipitation-frequency metrics, and stream-power inversion results. Together, these data provide the basis for evaluating how denudation varies with elevation, climate, and landscape form across sampled catchments on Pikes Peak. The second group of tables reports cosmogenic nuclide and erosion-model results used in the denudation analysis. Included files contain accelerator mass spectrometry (AMS) measurements for in situ-produced cosmogenic beryllium-10 (10Be), including sample identifiers, measured 10Be:9Be ratios, analytical uncertainties, carrier mass, quartz mass, blank corrections, blank-group statistics, and calculated 10Be concentrations and uncertainties. Additional tables summarize stream-power-law inversion results for sampled catchments, including optimized model parameters, predicted erosion rates, residual metrics, channel-pixel counts, and convergence status, as well as regression equations and summary statistics used to evaluate relationships among elevation, climate, frost cracking, precipitation forcing, and denudation rate. The package contains GIS files, comma-separated value files (.csv), Microsoft Excel files (.xlsx), CSV data dictionaries, a file-level metadata table, and a readme text file.

10Be cosmogenic nuclides↗