Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “generalized functions”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

GPU acceleration of local and semilocal density functional calculations in the SPARC electronic structure code

We present a Graphics Processing Unit (GPU)-accelerated version of the real-space SPARC electronic structure code for performing Kohn–Sham density functional theory calculations within the local density and generalized gradient approximations. In particular, we develop a modular math-kernel based implementation for NVIDIA architectures wherein the computationally expensive operations are carried out on the GPUs, with the remainder of the workload retained on the central processing units (CPUs). Here, using representative bulk and slab examples, we show that relative to CPU-only execution, GPUs enable speedups of up to 6× and 60× in node and core hours, respectively, bringing time to solution down to less than 30 s for a metallic system with over 14 000 electrons and enabling significant reductions in computational resources required for a given wall time.

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

Endothelial, cardiac muscle and skeletal muscle exhibit different viscous and elastic properties as determined by atomic force microscopy

This study evaluated the hypothesis that, due to functional and structural differences, the apparent elastic modulus and viscous behavior of cardiac and skeletal muscle and vascular endothelium would differ. To accurately determine the elastic modulus, the contribution of probe velocity, indentation depth, and the assumed shape of the probe were examined. Hysteresis was observed at high indentation velocities arising from viscous effects. Irreversible deformation was not observed for endothelial cells and hysteresis was negligible below 1 microm/s. For skeletal muscle and cardiac muscle cells, hysteresis was negligible below 0.25 microm/s. Viscous dissipation for endothelial and cardiac muscle cells was higher than for skeletal muscle cells. The calculated elastic modulus was most sensitive to the assumed probe geometry for the first 60 nm of indentation for the three cell types. Modeling the probe as a blunt cone-spherical cap resulted in variation in elastic modulus with indentation depth that was less than that calculated by treating the probe as a conical tip. Substrate contributions were negligible since the elastic modulus reached a steady value for indentations above 60 nm and the probe never indented more than 10% of the cell thickness. Cardiac cells were the stiffest (100.3+/-10.7 kPa), the skeletal muscle cells were intermediate (24.7+/-3.5 kPa), and the endothelial cells were the softest with a range of elastic moduli (1.4+/-0.1 to 6.8+/-0.4 kPa) depending on the location of the cell surface tested. Cardiac and skeletal muscle exhibited nonlinear elastic behavior. These passive mechanical properties are generally consistent with the function of these different cell types.

Non-NASA Center↗
DOE OSTI2025

Quadrupole strength in isobaric triplets

The dependence of the 𝐸⁢2 matrix elements on isospin projection 𝑇 𝑧 is linked to the conservation of the isospin symmetry. To study this conjecture, we calculated the 𝐵⁡(𝐸⁢2 : 2 + →0 + ) rates for the even-even 𝑇=1 mirror nuclei with 42 ≤ 𝐴 ≤ 98 within nuclear density functional theory, employing the generalized Bohr Hamiltonian, and carrying out angular momentum projection. We demonstrated that collective effects are crucial for describing experimental data near the 𝑁=𝑍 line without invoking explicit beyond-Coulomb isospin symmetry-breaking corrections. We also determined the 𝐵⁡(𝐸⁢2↓) values for odd-odd 𝑇 𝑧 =0 nuclei 70 Br and 78 Y in doubly blocked configurations. We discussed the requirements for accurately describing isobaric analog states and emphasized how current theoretical results should be interpreted within the study of isospin symmetry across isospin triplets.

collective levels↗
DOE OSTI2023

Uniform-density Bose-Einstein condensates of the Gross-Pitaevskii equation found by solving the inverse problem for the confining potential

Here, in this work, we consider a “reverse-engineering” approach to construct confining potentials that support exact, constant density kovaton solutions to the classical Gross-Pitaevskii equation (GPE) also known as the nonlinear Schr¨odinger equation (NLSE). In the one-dimensional case, the exact solution is the sum of stationary kink and anti-kink solutions, i.e. a kovaton, and in the overlapping region, the density is constant. In higher dimensions, the exact solutions are generalizations of this wave function. In the absence of self-interactions, the confining potential is similar to a smoothed out finite square well with minima also at the edges. When self-interactions are added, a term proportional to ±gψ*ψ gets added to the confining potential and ±gM, where M is the norm, gets added to the total energy. In the realm of stability analysis, we find (linearly) stable solutions in the case with repulsive self-interactions which also are stable to self-similar deformations. For attractive interactions, however, the minima at the edges of the potential get deeper and a barrier in the center forms as we increase the norm. This leads to instabilities at a critical value of M (related to the number of particles in the BEC). Comparing the stability criteria from Derrick’s theorem and Bogoliubov-de Gennes analysis stability results, we find that both predict stability for repulsive self-interactions and instability at a critical mass M for attractive interactions. However, the numerical analysis gives a much lower critical mass. The numerical analysis shows further that the initial instabilities violate the symmetry x → -x assumed by Derrick’s theorem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1989

Viscous flow in simple curved gaps. I - An asymptotic theory. II - Viscous stress and shape function

The present asymptotic theory for generalized incompressible two-dimensional steady flow in curved channels has been constructed in the limit when gas thickness approaches zero with its lateral dimensions fixed; successive asymptotic solution terms are analytically generated by quadratures. In the second part of this work, the curvature of the gap treated is arbitrary. It is established that each term in the series solution of velocity and pressure is the product of a scale factor and a universal shape functions. Various interaction modes between the volume rate-of-flow, curvature, and its variations, are identified and quantitatively characterized.

Fan, D.-N.↗
DOE OSTI2022

Interband Transitions and Critical Points of Single-Crystal Thoria Compared with Urania

The interband transitions of UO 2 are validated independently through cathode luminescence. A picture emerges consistent with density functional theory. While theory is generally consistent with experiment, it is evident from the comparison of UO 2 and ThO 2 that the choice of functional can significantly alter the bandgap and some details of the band structure, in particular at the conduction band minimum. Strictly ab initio predictions of the optical properties of the actinide compounds, based on density functional theory alone, continue to be somewhat elusive.

42 ENGINEERING↗
DOE OSTI2024

Amplitudes, supersymmetric black hole scattering at $\mathcal{O}\left({G}^5\right)$, and loop integration

We compute the potential-graviton contribution to the scattering amplitude, the radial action, and the scattering angle of two extremal black holes in $\mathcal{N}$ = 8 supergravity at the fifth post-Minkowskian order and to next-to-leading order in a large mass expansion (first self-force order). Properties of classical unitarity cuts allow us to focus on the integration-by-parts reduction of planar integrals, while nonplanar integrals at this order are obtained from the planar ones by straightforward manipulations. We present the solution to the differential equations for all master integrals necessary to evaluate the classical scattering amplitudes of massive scalar particles at this order in all gravitational theories, in particular in $\mathcal{N}$ = 8 supergravity, and in general relativity. Despite the appearance of higher-weight generalized polylogarithms and elliptic functions in the solution to the differential equation for master integrals, the final supergravity answer is remarkably simple and contains only (harmonic) polylogarithmic functions up to weight 2. The systematic analysis of elliptic integrals discussed here, as well as the particular organization of boundary integrals in $\mathcal{N}$ = 8 observables are independent of supersymmetry and may have wider applications, including to aspects of collider physics.

Black Holes↗
DOE OSTI2020

Electronic structure and thermoelectric properties of CoAsSb with post-DFT approaches

Here, we study the electronic structure and thermoelectric properties of recently synthesized CoAsSb. The calculated bandgap becomes more accurate for increasingly-complex electronic structure methods: generalized gradient approximation, hybrid functionals, selfconsistent linearized quasiparticle GW method (LQSGW), and LQSGW with simplified vertex corrections. By equating the bandgaps of each method from a rigid shift of the bands, we evaluate the contributions made to thermoelectric properties beyond the bandgap. In doing so, we evaluate the efficacy of a commonpractice: a rigid shift applied to less-costly electronic structure methods in order to gain some accuracy of the more-costly methods. We find that while the shift made the Seebeck coefficients much closer to one another than from the original electronic structures, there remain differences in the Goldsmid-Sharp (thermoelectric) bandgap between the methods and from the intended electronic bandgap. Additionally, some lasting differences in temperature dependence remain between the methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Physics-guided neural networks with engineering domain knowledge for hybrid process modeling

As neural networks are more frequently used to solve problems in science and engineering, the methods used to incorporate scientific knowledge into these networks are becoming increasingly complex. Here, this work breaks down these complicated techniques into a set of basic strategies which can easily be applied to diverse situations. Several novel neural networks are built using the categories laid out in this work. These networks are tested on simulated data from a continuous stirred tank reactor (CSTR) model to evaluate the advantages provided by each network. The three points demonstrated in this work are: (1) architectural hybrid models can speed up convergence and reduce the amount of data necessary to train a model; (2) adding a physics-guided loss function can improve model generalization and make models more physically consistent; (3) using physics-guided initialization and transfer learning improves accuracy and speeds up convergence, but can harm generalizability if used incorrectly.

42 ENGINEERING↗
DOE OSTI2023

Reviewing flexibility in industrial electrification: Focusing on green ammonia and steel in the United States

The renewable energy transition in the power sector involves a paradigm shift for flexibility. Flexible demand is more attainable than before, while supply flexibility faces new constraints due to the increased share of variable renewable resources. Increased demand flexibility can allow less use of peaking power plants and avoid need for additional capacity. Industrial flexibility could be especially interesting in this regard, as industrial customers are larger on average than other customers and have typically provided the largest share of demand response in the United States. We consider industrial demand, studying characteristics of flexible industrial loads. We then examine the dynamics of change occurring around industrial load flexibility by focusing on two case studies: green ammonia and steel production via electric arc furnaces. Electric arc furnace steel production is an important component of current demand response programs, whereas green ammonia and green fuels offer new paradigms for flexibility. We analyze the structure and functions of the technological innovation systems of load flexibility in those two industries via interviews with twenty-two stakeholders. Here, we conclude that in the United States these technological innovation systems are not well-functioning for industry in general or for steel, but do seem to be present for green ammonia. Additionally, explicit connections between scope two greenhouse gas emissions (those from purchased energy) and flexibility are lacking, and likewise industry stakeholders do not appear to make a connection between decarbonization and load flexibility, thus flexible demand is not viewed as a tool in industrial decarbonization.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2023

Fast model-based scenario optimization in NSTX-U enabled by analytic gradient computation

Model-based optimization offers a systematic approach to advanced scenario planning. In this case, the feedforward-control inputs (actuator trajectories) that are needed to attain and sustain a desired scenario are obtained by solving a nonlinear constrained optimization problem. This class of problems generally minimize a cost function that measures the difference between desired and actual plasma states. Several numerical optimization algorithms, such as sequential quadratic programming, require repeated calculation of the cost function gradients with respect to the input trajectories. Calculating these gradients numerically can be computationally intensive, increasing the time needed to solve the feedforward-control optimization problem. Here, this work introduces a method to analytically calculate these cost function gradients from the current profile evolution model. This can significantly reduce the computational time and allow for fast feedforward-control optimization, which would eventually enable optimal scenario planning between discharges. The performance of the feedforward optimizer with analytical gradients is compared to a traditional optimization algorithm based on numerical gradients for different NSTX-U scenarios. The plasma dynamics in the optimization algorithm are simulated using the Control Oriented Transport SIMulator (COTSIM). Results of the work show that analytical gradients consistently reduce the computation time while achieving trajectories that are comparable to those obtained by traditional optimization algorithms based on numerical gradients.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Numerical schemes for a multi-species BGK model with velocity-dependent collision frequency

Here, we consider a kinetic description of multi-species gas mixture modeled with Bhatnagar-Gross-Krook (BGK) collision operators, in which the collision frequency varies not only in time and space but also with the microscopic velocity. In this model, the Maxwellians typically used in standard BGK operators are replaced by a generalization of such target functions, which are defined by a variational procedure. In this paper we present a numerical method for simulating this model, which uses an Implicit-Explicit (IMEX) scheme to minimize a certain potential function, mimicking the Lagrange functional that appears in the theoretical derivation. We show that theoretical properties such as conservation of mass, total momentum and total energy as well as positivity of the distribution functions are preserved by the numerical method, and illustrate its usefulness and effectiveness with numerical examples.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Phase-field modeling of dislocation–interstitial interactions

The mechanical behavior of body-centered cubic (BCC) materials can be dramatically affected by the presence of interstitial solute atoms. Here we present a new phase-field dislocation dynamics formulation to include the diffusion of interstitials. Short-range interactions are accounted for by a concentration-dependent lattice energy, and long-range interactions are accounted for by modifications to the elastic energy. The interstitial diffusion law introduces gradients that require methods for minimizing Gibbs oscillations, which is done via a modified Green’s function. The formulation is general to any solute-solvent system and is applied here to Nb-O as a model system, whose interstitial parameters are obtained from ab initio calculations. The effect of O on the core structures of Nb edge and screw dislocations is calculated. The diffusion of O to form interstitial atmospheres around dislocation cores is simulated, as well as the critical stresses required for dislocations to break away or cross slip from these atmospheres. Future applications of the method to simulate complex interstitial embrittlement mechanisms are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Leptodermic corrections to the TOV equations and nuclear astrophysics within the effective surface approximation

The macroscopic model for a neutron star (NS) as a liquid drop at the equilibrium is used to extend the Tolman-Oppenheimer-Volkoff (TOV) equations taking into account the gradient terms responsible for the system surface. The parameters of the Schwarzschild metric in the spherical case are found with these surface corrections to the known leading (zero) order of the leptodermic approximation a/R << 1, where a is the NS effective-surface (ES) thickness, and R is the effective NS radius. The energy density $\mathscr{E}$ is considered in a general form including the functions of the particle number density and of its gradient terms. The macroscopic gravitational component $Φ$(ρ) of the energy density is taken into account in the simplest form as expansion in powers of $ρ$ – $\overline{ρ}$, where $\overline{ρ}$ is the saturation density, up to second order, in terms of its contributions to the separation particle energy and incompressibility. Density distributions ρ across the NS ES in the normal direction to the ES, which are derived in the simple analytical form at the same leading approximation, was used for the derivation of the modified TOV (MTOV) equations by accounting for their NS surface corrections. As a result, the MTOV equations are analytically solved at first order and the results are compared with the standard TOV approach of the zero order.

Magner, A. G. [Institute for Nuclear Research, Kyi↗
DOE OSTI2022

Exploring Electrosynthesis: Bulk Electrolysis and Cyclic Voltammetry Analysis of the Shono Oxidation

As electrochemistry continues to gain broader acceptance and use within the organic chemistry community, it is important that advanced undergraduate students are exposed to fundamental and practical knowledge of electrochemical applications for chemical synthesis. Herein, we describe the development of an undergraduate laboratory experience that introduces synthetic and analytical electrochemistry concepts to an advanced organic chemistry class. Experiments focus on the electrooxidative α-functionalization of carbamates, more generally known as the Shono oxidation, and include cyclic voltammetry analysis of two cyclic carbamates and a constant current bulk electrolysis reaction. Here, the exercise offers students an authentic experience in organic electrochemistry, lays a practical and theoretical foundation for future engagement with concepts in electrochemistry and redox chemistry, and strengthens fundamental organic chemistry skills.

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

Norm-Conserving Pseudopotentials and Basis Sets to Explore Actinide Chemistry in Complex Environments

We constructed and tested a new set of norm-conserving pseudopotentials, as well as companion Gaussian basis sets, for the actinide series (Ac–Lr) using the Goedecker, Teter and Hutter (GTH) formalism within the generalized gradient approximation (GGA) functional for exchange-correlation of Perdew, Burke and Ernzerhof (PBE). This fills a critical void in the literature for studying the chemistry of 5f–block elements in condensed phase. The accuracy and reliability of the newly parameterized An–GTH pseudopotentials and basis sets, a variety of tests on actinide-containing molecules were carried out and compared to all-electron and available experimental results. The new pseudopotentials include both medium ( [Xe]4f 14 ) and large core ( [Xe]4f 14 5d 10 ) options that successfully reproduced structure and energetics, particularly redox processes. The medium core size set, in particular, reproduced all electron calculations over multiple oxidation states from 0 to VII, whereas the large core set was suitable only for the early series elements and low oxidation states.

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

Computational Study of Triphosphine-Ligated Cu(I) Catalysts for Hydrogenation of CO 2 to Formate

The catalyzed hydrogenation of CO 2 to formate via a triphosphine-ligated Cu(I) was studied computationally. Two bases, DBU and TBD, were studied in the context of two proposed mechanisms in MeCN solvent. Of the four functionals benchmarked, M06 was generally in the best agreement with the experimentally estimated values. Activation of H 2 was explored by using LCu(DBU) + to form LCuH. Dissociation of a ligand arm results in higher barriers to form the key hydride complex, LCuH. There is no significant difference between the choice of base, DBU or TBD with respect to the proposed mechanisms. We propose that the experimentally observed differences between DBU and TBD reactivity for this mechanism are due to off-pathway changes.

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

Reversing the Irreversible: Thermodynamic Stabilization of LiAlH 4 Nanoconfined Within a Nitrogen-Doped Carbon Host

A general problem when designing functional nanomaterials for energy storage is the lack of control over the stability and reactivity of metastable phases. Using the high-capacity hydrogen storage candidate LiAlH4 as an exemplar, we demonstrate an alternative approach to the thermodynamic stabilization of metastable metal hydrides by coordination to nitrogen binding sites within the nanopores of N-doped CMK-3 carbon (NCMK-3). The resulting LiAlH 4 @NCMK-3 material releases H 2 at temperatures as low as 126 °C with full decomposition below 240 °C, bypassing the usual Li 3 AlH 6 intermediate observed in bulk. Moreover, >80% of LiAlH 4 can be regenerated under 100 MPa H 2 , a feat previously thought to be impossible. Nitrogen sites are critical to these improvements, as no reversibility is observed with undoped CMK-3. Density functional theory predicts a drastically reduced Al–H bond dissociation energy and supports the observed change in the reaction pathway. Finally, the calculations also provide a rationale for the solid-state reversibility, which derives from the combined effects of nanoconfinement, Li adatom formation, and charge redistribution between the metal hydride and the host.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗