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

Eulerian thermo-mechanical simulations of heterogeneous solid propellants using an approximate projection method

We present a mathematical model that describes thermo-mechanical deformations and thermal gradients on the unsteady burning of a heterogeneous solid propellant. A scaling study shows that the deformations in solid at combustion timescales can be treated as quasi-static. The resulting thermo-mechanical formulation is formulated on a Cartesian grid and makes use of a weak form of Chorin-type projection method to deal with large difference in shear modulus of constituent materials. Further, a one-dimensional verification study is carried out by comparing numerical simulations with those of an analytical model. In addition, convergence studies for a two-dimensional propellant sandwich configuration are presented for the stress, velocity, and reference map components. Finally, simulations are carried out for a two-dimensional random propellant pack and the time-averaged burn rate and strains are reported.

42 ENGINEERING↗
DOE OSTI2022

Systematic Errors Induced by Dipole Component Change in Quadrupole Beam-Based Alignment Measurements for Circular Accelerators

Beam-based alignment (BBA) for quadrupole magnets is now a well-established process for circular accelerators to steer beam orbit through the magnetic centers such that beam orbit is not disturbed when the gradients of quadrupoles are changed. At some facilities, the centers estimated from BBA techniques are also used as a proxy for steering beam through the centers of nearby nonlinear magnets to minimize degradation of nonlinear dynamics caused by off-centered orbit. The random errors associated with BBA are well known, but very little attention has been paid to the systematic error generated by a dipole component change, or equivalently a shift in the magnetic center, when the quadrupole strength is varied. The analytical formulas for this error and its amplification factor (potentially a few orders of magnitude) with respect to the magnetic center motion have been derived and confirmed with simulations. We address the significance of this error as well as the part that cannot be estimated from the beam-based measurement alone. Finally, the BBA measurements obtained at NSLS-II and their systematic errors on the order of hundreds of microns are discussed.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Numerical investigation of active control of tearing mode by magnetic coils and the role of Δ'

Magnetic feedback stabilization of the tearing mode (TM) is numerically investigated, utilizing the MARS-F code (Liu et al 2000 Phys. Plasmas 7 3681) for toroidal tokamak equilibria. With control coil configurations assumed in this study, magnetic feedback partially or fully stabilizes the TM, with either vanishing or finite equilibrium pressure. The best control is achieved by the combination of internal active coils and internal poloidal sensors. The internal and external tearing indices are evaluated for the close-loop system, based on the MARS-F computed mode eigenvalue and eigenfunction, respectively. In the absence of the favorable curvature effect, these two indices are real-valued and quantitatively agree well with each other. For the equilibrium with finite pressure gradient at the mode rational surface, the favorable average curvature effect becomes important and the close-loop tearing index also becomes complex-valued, partly due to interaction of the feedback system with the dissipative wall eddy current response. Isolating the inner layer and outer region response to magnetic feedback, with either proportional or proportional-derivative actions, allows to establish that feedback stabilization of the TM occurs mainly due to modification of the behavior of the external ideal solution, further confirming the analytic result reported in He et al 2021 Phys. Plasmas 28 012504.

Physics↗
DOE OSTI2025

Flux Hypothesis for Odd Transport Phenomena

Onsager's regression hypothesis makes a fundamental connection between macroscopic transport phenomena and the average relaxation of spontaneous microscopic fluctuations. This relaxation, however, is agnostic to odd transport phenomena, in which fluxes run orthogonal to the gradients driving them. To account for odd transport, we generalize the regression hypothesis, postulating that macroscopic linear constitutive laws are, on average, obeyed by microscopic fluctuations, whether they contribute to relaxation or not. From this "flux hypothesis," Green-Kubo and reciprocal relations follow, elucidating the separate roles of broken time-reversal and parity symmetries underlying various odd transport coefficients. As an application, we derive and verify the Green-Kubo relation for odd collective diffusion in chiral active matter, first in an analytically tractable model and subsequently through molecular dynamics simulations of concentrated active spinners.

Hargus, Cory↗
DOE OSTI2020

TeraChem Cloud: A High-Performance Computing Service for Scalable Distributed GPU-Accelerated Electronic Structure Calculations

The encapsulation and commoditization of electronic structure arise naturally as interoperability, and the use of nontraditional compute resources (e.g., new hardware accelerators, cloud computing) remains important for the computational chemistry community. Here, we present TERACHEM CLOUD, a high-performance computing service (HPCS) that offers on-demand electronic structure calculations on both traditional HPC clusters and cloud-based hardware. The framework is designed using off-the-shelf web technologies and containerization to be extremely scalable and portable. Within the HPCS model, users can quickly develop new methods and algorithms in an interactive environment on their laptop while allowing TERACHEM CLOUD to distribute ab initio calculations across all available resources. This approach greatly increases the accessibility of hardware accelerators such as graphics processing units (GPUs) and flexibility for the development of new methods as additional electronic structure packages are integrated into the framework as alternative backends. Cost-performance analysis indicates that traditional nodes are the most cost-effective long-term solution, but commercial cloud providers offer cutting-edge hardware with competitive rates for short-term large-scale calculations. We demonstrate the power of the TERACHEM CLOUD framework by carrying out several showcase calculations, including the generation of 300,000 density functional theory energy and gradient evaluations on medium-sized organic molecules and reproducing 300 fs of nonadiabatic dynamics on the B800-B850 antenna complex in LH2, with the latter demonstration using over 50 Tesla V100 GPUs in a commercial cloud environment in 8 h for approximately $1250.

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

Predictive understanding of the surface tension and velocity of sound in ionic liquids using machine learning

Knowledge of the physical properties of ionic liquids (ILs), such as the surface tension and speed of sound, is important for both industrial and research applications. Unfortunately, technical challenges and costs limit exhaustive experimental screening efforts of ILs for these critical properties. Previous work has demonstrated that the use of quantum-mechanics-based thermochemical property prediction tools, such as the conductor-like screening model for real solvents, when combined with machine learning (ML) approaches, may provide an alternative pathway to guide the rapid screening and design of ILs for desired physiochemical properties. However, the question of which machine-learning approaches are most appropriate remains. In the present study, we examine how different ML architectures, ranging from tree-based approaches to feed-forward artificial neural networks, perform in generating nonlinear multivariate quantitative structure–property relationship models for the prediction of the temperature- and pressure-dependent surface tension of and speed of sound in ILs over a wide range of surface tensions (16.9–76.2 mN/m) and speeds of sound (1009.7–1992 m/s). The ML models are further interrogated using the powerful interpretation method, shapley additive explanations. We find that several different ML models provide high accuracy, according to traditional statistical metrics. The decision tree-based approaches appear to be the most accurate and precise, with extreme gradient-boosting trees and gradient-boosting trees being the best performers. However, our results also indicate that the promise of using machine-learning to gain deep insights into the underlying physics driving structure–property relationships in ILs may still be somewhat premature.

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

Transport Barriers in magnetized plasmas- general theory with dynamical constraints

Abstract A fundamental dynamical constraint—that fluctuation induced charge-weighted particle flux must vanish- can prevent instabilities from accessing the free energy in the strong gradients characteristic of Transport Barriers (TBs). Density gradients, when large enough, lead to a violation of the constraint and hence preclude unstable modes and turbulent transport. This mechanism, then, broadens the class of configurations (in magnetized plasmas) where these high confinement states can be formed and sustained. The need for velocity shear, the conventional agent for TB formation, is obviated. The most important ramifications of the constraint is to permit a charting out of the domains conducive to TB formation and hence to optimally confined fusion worthy states; the detailed investigation is conducted through new analytic methods and extensive gyrokinetic simulations.

Physics↗
DOE OSTI2022

Energy-Efficient Microbial Electrochemical Lignin and Alkaline Hydroxide Recovery from DMR Black Liquor

Black liquor is a complex waste stream generated from the fast-expanding biorefinery industry, and the high management cost has been known as a major development barrier. In this work, we demonstrated a rationally designed 4-chamber microbial electrolysis process that was capable of recovering the value-added lignin, NaOH, and H2 products while removing waste organics from the Deacetylation and Mechanical Refining (DMR) black liquor, which could potentially reduce the disposal cost by up to 40%. The lab-scale MEC achieved 52% organic removal from a highly concentrated feedstock (28,960 ± 212 mg/L COD) and recovered 83.2 mM NaOH under a small voltage application (1.2 V). Lignin recovery was observed to be as high as 13.98 +/- 0.28 g/L, and 0.35 L/L/day H2 was produced with a faradaic efficiency of 93.5%. This work provided a good example on how complex liquid waste stream can be valorized via multi-functional microbial electrochemical technology, and how the traditionally unfavored pH gradient in an electrochemical reactor can be manipulated for lignin precipitation and alkaline recovery at the same time. The versatility and good performance in concurrent contaminants removal, lignin precipitation, H2 production and NaOH recovery demonstrate a good potential for microbial electrolysis in black liquor treatment.

BIOMASS FUELS,INORGANIC, ORGANIC, PHYSICAL, AND AN↗
DOE OSTI2021

Fully variational incremental CASSCF

We report the complete-active-space self-consistent field (CASSCF) method is a canonical electronic structure theory that holds a central place in conceptualizing and practicing first principles simulations. For application to realistic molecules, however, the CASSCF must be approximated to circumvent its exponentially scaling computational costs. Applying the many-body expansion - also known as the method of increments - to CASSCF (iCASSCF) has been shown to produce a polynomially scaling method that retains much of the accuracy of the parent theory and is capable of treating full valence active spaces. Due to an approximation made in the orbital gradient, the orbital parameters of the original iCASSCF formulation could not be variationally optimized, which limited the accuracy of its nuclear gradient. Herein, a variational iCASSCF is introduced and implemented, where all parameters are fully optimized during energy minimization. This method is able to recover electronic correlations from the full valence space in large systems, produce accurate gradients, and optimize stable geometries as well as transition states. Demonstrations on challenging test cases, such as the oxoMn(salen)Cl complex with 84 electrons in 84 orbitals and the automerization of cyclobutadiene, show that the fully variational iCASSCF is a powerful tool for describing challenging molecular chemistries.

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

The effect of curvature induced broken potential vorticity conservation on drift wave turbulences

Abstract The correlation theory of turbulence suppression (Zhang and Mahajan 1993 Phys. Fluids B5 2000) by velocity shear was constructed by invoking the ansatz of potential vorticity conservation (PVC) that holds for relatively simple (slab) models of drift wave turbulence. It is, therefore, surprising that a detailed modern simulation of the H mode pedestal, using the gyrokinetic code GENE (Hatch et al 2018 Plasma Phys. Control. Fusion 60 084003), found ‘striking agreement’ with the predictions of the analytic model. To understand the reasons for this remarkable agreement, an extended theory that contains finite (magnetic) curvature, and which does not conserve potential vorticity, is developed and ‘solved’ by calculating an inhomogeneous Green function reflecting the fact that the new system has a potential vorticity source. It is, then, demonstrated that the effect of the broken PVC is insignificant for the normal operation parameters in tokamaks; the correction due to curvature is at the order of 2 L 0 / R , where L 0 is the scale length of the local gradient and R is the major radius. The excellent agreement between simulation and slab model 1993 theory is therefore not accidental; the latter can be applied with confidence to the tokamak pedestal.

Physics↗
DOE OSTI2022

Unconventional Charge Transport in MgCr 2 O 4 and Implications for Battery Intercalation Hosts

Ion transport in solid-state cathode materials prescribes a fundamental limit to the rates batteries can operate; therefore, an accurate understanding of ion transport is a critical missing piece to enable new battery technologies, such as magnesium batteries. Based on our conventional understanding of lithium-ion materials, MgCr 2 O 4 is a promising magnesium-ion cathode material given its high capacity, high voltage against an Mg anode, and acceptable computed diffusion barriers. Electrochemical examinations of MgCr 2 O 4 , however, reveal significant energetic limitations. Motivated by these disparate observations; herein, we examine long-range ion transport by electrically polarizing dense pellets of MgCr 2 O 4 . Our conventional understanding of ion transport in battery cathode materials, e.g., Nernst-Einstein conduction, cannot explain the measured response since it neglects frictional interactions between mobile species and their nonideal free energies. In this work, we propose an extended theory that incorporates these interactions and reduces to the Nernst-Einstein conduction under dilute conditions. This theory describes the measured response, and we report the first study of long-range ion transport behavior in MgCr 2 O 4 . We conclusively show that the Mg chemical diffusivity is comparable to lithium-ion electrode materials, whereas the total conductivity is rate-limiting. Given these differences, energy storage in MgCr 2 O 4 is limited by particle-scale voltage drops, unlike lithium-ion particles that are limited by concentration gradients. Future materials design efforts should consider the interspecies interactions described in this extended theory, particularly with respect to multivalent-ion systems and their resultant effects on continuum transport properties.

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

Computationally inexpensive part-scale thermal history of additive friction-stir deposition

This study presents an analytical model for steady-state power generation and tool heat loss in additive friction-stir deposition (AFSD), developed to enable part-scale thermal simulation while remaining computationally inexpensive. The model predicts total generated power, yielding 3.7–4.7 kW across deposition temperature setpoints of 400–460 °C for the deposition of AA6061 with a Be-Cu tool. This corresponds to 90–95% of the reported spindle power. Tool heat loss is experimentally determined by calibrating a steady-state energy balance between the generated power, the substrate-deposition thermal gradient, and a temperature dependent tool heat loss term: q tool (T) = a + b (T - 400°C) with a = 2.7 x 10 6 Wm -2 and b = 9.5 x 10 3 Wm -2 K -1 . The calibration indicates that about 69% of the generated heat is conducted into the tool for this configuration, which is much higher than previously reported. The calibrated heat-source is implemented in finite element software (Adamantine) to simulate the transient thermal history of a 100 cm 3 representative build in 8 min on a standard desktop (at 0.635 mm build-height resolution). For the first three layers, the substrate temperatures between simulation and experiment are within 10% mean absolute percentage error. Sensitivity analysis indicates that uncertainties in average deposition temperature and deformation localization (stir-zone geometry, depth, and spatial dependance of strain-rate and flow stress) dominate model variance, motivating additional experimental verification.

Additive Friction-Stir Deposition↗
DOE OSTI2022

Chemical design of self-propelled Janus droplets

Solubilizing, self-propelling droplets have emerged as a rich chemical platform for the exploration of active matter, but isotropic droplets rely on spontaneous symmetry breaking to sustain motion. The introduction of permanent asymmetry, e.g., in the form of a biphasic Janus droplet, has not been explored as a comprehensive design strategy for active droplets, despite the widespread use of Janus structures in motile solid particles. Here, we uncover the chemomechanical framework underlying the self-propulsion of biphasic Janus oil droplets solubilizing in aqueous surfactant. We elucidate how droplet propulsion is influenced by the degree of oil mixing, droplet shape, and oil solubilization rates for a range of oil combinations. In addition, spatiotemporal control over droplet swimming speed and orientation is demonstrated through the application of thermal gradients applied via joule heating and laser illumination. Finally, we also explore the interactions between collections of Janus droplets, including the spontaneous formation of spinning multi-droplet clusters.

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

Prediction of residual stresses in additively manufactured parts using lumped capacitance and classical lamination theory

Several industries are interested in Laser Powder Bed Fusion (L-PBF) Additively Manufactured (AM) metal parts because their designs can be made arbitrarily complex while retaining bulk-type material properties. However, the residual stresses (RS) and distortions caused by the heat gradients inherent to L-PBF processes are detrimental to the structural integrity of the parts and must be taken into consideration during the part design cycle. Predicting the state of stresses in as-built 3D printed parts is a difficult problem that is typically approached with the use of transient thermomechanical Finite Element Models (FEMs). However, the nonlinearities associated with AM processes are difficult to capture in these FEMs without increasing the computational cost of the simulation, limiting their ability to be incorporated into practical design cycles. This work presents a novel analytical framework that combines lumped capacitance nonlinear heat transfer with time dependent classical lamination theory to efficiently and accurately predict RS in as-built L-PBF parts without the need of FEMs. The simulation was compared to Neutron Diffraction (ND) residual strain measurements taken at Oak Ridge National Laboratories (ORNL) as well as Synchrotron X-ray Diffraction (XRD) strain data published by the National Institute of Standards and Technology (NIST). The simulation predictions and the experimental data showed excellent agreement for the in-plane strain directions, and general agreement for the out of plane strain component, highlighting an area where further development can be implemented.

42 ENGINEERING↗
DOE OSTI2020

Anomalously Large Seebeck Coefficient of CuFeS 2 Derives from Large Asymmetry in the Energy Dependence of Carrier Relaxation Time

The Seebeck effect in a material originates from the distribution of asymmetry in the electron transport under a temperature gradient, which has contributions from the energy-dependent electronic density-of-states and carrier mobility. However, because the energy dependence of common electron scattering mechanisms is weak, the mobility-driven Seebeck coefficient has long been ignored in most thermoelectric materials, and the energy asymmetry of the density-of-states has been considered the dominant contribution. In this work, we describe a hopping transport behavior observed in CuFeS 2 , and a large carrier Hall mobility gradient of dμH/dT that creates an unusually large energy-dependent mobility contribution to the Seebeck coefficient. This work offers several ideas regarding the mobility-driven Seebeck effect and its potential utilization in the design of thermoelectric materials.

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

Temperature Dependence of Nuclear Quadrupole Resonance and the Observation of Metal–Ligand Covalency in Actinide Complexes: 35 Cl in Cs 2 UO 2 Cl 4

We report a study of the temperature dependence of 35 Cl nuclear quadrupole resonance (NQR) transition energies and spin–lattice relaxation times (T 1 ) for 235 U-depleted dicesium uranyl tetrachloride (Cs 2 UO 2 Cl 4 ) aimed at elucidating electronic interactions between the uranium center and atoms in the equatorial plane of the UO 2 2+ ion. The transition frequency decreases slowly with temperature below 75 K and with a more rapid linear dependence above this temperature. The spin–lattice relaxation time becomes shorter with temperature, and as temperatures increase, the T 1 decrease becomes nearly quadratic. The observed trends are reproduced by a model that assumes phonon-induced fluctuations of the electric field gradient tensor and partial electron delocalization from Cl to U. The fit of the theoretical model to experimental data allows a Debye temperature of 96 K to be estimated. Finally, the generalization of this approach to investigations of covalency in actinide–ligand bonding is examined.

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

Quest for an Optimal Methane Hydrate Formation in the Pores of Hydrolytically Stable Metal–Organic Frameworks

Porous metal–organic frameworks (MOFs) capable of storing a relatively high amount of dry methane (CH 4 ) in the adsorbed phase are largely explored; however, solid CH 4 storage in confined pores of MOFs in the form of hydrates is yet to be discovered. Here we report a rational approach to form CH 4 hydrates by taking advantage of the optimal pore confinement in relatively narrow cavities of hydrolytically stable MOFs. Unprecedentedly, we were able to isolate methane hydrate (MH) nanocrystals with an sI structure encapsulated inside MOF pores with an optimal cavity dimension. Furthermore, it was found that confined nanocrystals require cavities slightly larger than the unit cell crystal size of MHs (1.2 nm), as exemplified in the experimental case study performed on Cr- soc -MOF-1 vs smaller cavities of Y- shp -MOF-5. Under these conditions, the excess amount of methane stored in the pores of Cr- soc -MOF-1 in the form of MH was found to be ≈50% larger than the corresponding dry adsorbed amount at 10 MPa. More importantly, the pressure gradient driving the CH 4 storage/delivery process could be drastically reduced compared to the conventional CH 4 -adsorbed phase storage on the dry Cr- soc -MOF-1 (≤3 MPa vs 10 MPa).

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

Dimensional interpolation for metallic hydrogen

In this work, we employ a simple and mostly accurate dimensional interpolation formula using dimensional limits D = 1 and D = ∞ to obtain D = 3 ground-state energy of metallic hydrogen. We also present results describing the phase transitions for different symmetries of three-dimensional structure lattices. The interpolation formula not only predicts fairly accurate energies but also predicts a correct functional form of the energy as a function of the lattice parameters. That allows us to calculate different physical quantities such as the bulk modulus, Debye temperature, and critical transition temperature, from the gradient and the curvature of the energy curve as a function of the lattice parameters. These theoretical calculations suggest that metallic hydrogen is a likely candidate for high temperature superconductivity. The dimensional interpolation formula is robust and might be useful to obtain the energies of complex many-body systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗