SEARCH · Search NASA
Results for “convention equation”
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.
Asymmetries in ionization of atomic superposition states by ultrashort laser pulses
Abstract Progress in ultrafast science allows for probing quantum superposition states with ultrashort laser pulses in the new regime where several linear and nonlinear ionization pathways compete. Interferences of pathways can be observed in the photoelectron angular distribution and in the past they have been analyzed for atoms and molecules in a single quantum state via anisotropy and asymmetry parameters. Those conventional parameters, however, do not provide comprehensive tools for probing superposition states in the emerging research area of bright and ultrashort light sources, such as free-electron lasers and high-order harmonic generation. We propose a new set of generalized asymmetry parameters which are sensitive to interference effects in the photoionization and the interplay of competing pathways as the laser pulse duration is shortened and the laser intensity is increased. The relevance of the parameters is demonstrated using results of state-of-the-art numerical solutions of the time-dependent Schrödinger equation for ionization of helium atom and neon atom.
Port of New York and New Jersey Drayage Electrification Analysis
The National Renewable Energy Laboratory (NREL) evaluated the potential for drayage electrification in the Port of New York and New Jersey (PoNYNJ), with a focus on operators: Harbor Freight Transport (HF), Safeway Trucking (SWT), and International Motor Freight Inc (IMF). This report summarizes the data collection and electrification evaluation of all three drayage operators, includes detailed operational data, and identifies the performance requirements for battery electric tractors (BETs) and corresponding infrastructure operated within the context of PoNYNJ drayage operation. This report also details a methodology to evaluate opportunities, strategies, and challenges associated with future expansions of BETs in meeting PANYNJ emissions goals. The Port Authority has established a goal of achieving Net Zero carbon emissions by 2050 across all facilities, including from tenant and stakeholder sources such as drayage trucks. NREL used real-world performance data collected on the three PoNYNJ drayage operations, along with modeling and analysis tools to compare BET to diesel trucks. From March to July 2021, NREL collected 1Hz vehicle and engine data from 46 drayage trucks at the three operators totaling nearly 121,000 miles of operation, providing enough information to assess vehicle operations for electrification potential. A Future Automotive Systems Technology Simulator (FASTSim) electric truck powertrain model was validated using PoNYNJ data and scenarios were run to evaluate drayage truck electrification requirements over the real-world cycles. The first scenario examined BET viability with minimal changes to existing operations. This assumes the trucks charge when stopped for two hours or longer, have a functional battery size of 375 kWh, and can charge at 270 kilowatts (kW) average which are the specification of the commercially available Freightliner eCascadia. The second scenario looked at what operational, charging infrastructure, and BET technology changes would be needed to fully electrify. Finally, detailed analysis was run on charging rate structure to understand operational costs to the fleets. The studied drayage trucks averaged 5.1 MPG, spent roughly 9% of their energy at idle, and drove an average of 140 miles per day with a maximum daily distance of 573 miles. The FASTSim model results indicate a comparable BET would use 417 kWh of energy per day on average accounting for cargo weight, which is close to the full usable capacity of the eCascadia currently available on the market. Based on the daily average operating data, partial fleet electrification is possible with current technology. However, some specific days of operation would require over 1,600 kWh of energy due to longer distances traveled by the trucks and more intense operation. Trucks used for long distance and intense operation cannot be readily electrified with current technology without operational changes. Full adoption of BETs could reduce CO 2 emissions from these fleets by roughly 75% today, eliminating 76 metric tons of CO 2 (MTCO 2 ) per vehicle each year, which equates to 24,100 MTCO 2 per year for all three operators. Commercially available direct current fast chargers (DCFC) have charge rates up to 350 kW. Based on the average daily modeled energy use for each operator, current industrial rate structures, and the assumption of 350 kW peak charging, full drayage electrification would increase electricity consumption. In addition, peak demand usage would increase with unmanaged charging along with cost of electricity having a direct impact on cost per mile for electric vehicles. The resulting cost per mile for BETs along with comparable cost per mile for conventional diesel trucks are also examined at $\$$4.00 per gallon of diesel. It will be important for PANYNJ and the drayage operators within the PoNYNJ to consider these load impacts to their existing electrical infrastructure and devise operational strategies that avoid coincident charging of vehicles to mitigate demand charges. Despite these electricity cost increases, savings from reductions in diesel consumption will help offset the costs of this increased electricity consumption. However, prices of both electricity and diesel are subject to change based on various factors meaning the realized savings will vary over time. This shows BETs could be cost-competitive on an energy cost per mile basis for all scenarios while diesel is above $\$$3.00/gal. Further, if diesel prices dropped to the 15-year low of $2.33/gal, it would still be cost competitive to operate the EVs with electricity costs of 16.3 ¢/kWh or less.
Machine learning of superconducting critical temperature from Eliashberg theory
Abstract The Eliashberg theory of superconductivity accounts for the fundamental physics of conventional superconductors, including the retardation of the interaction and the Coulomb pseudopotential, to predict the critical temperature T c . McMillan, Allen, and Dynes derived approximate closed-form expressions for the critical temperature within this theory, which depends on the electron–phonon spectral function α 2 F ( ω ). Here we show that modern machine-learning techniques can substantially improve these formulae, accounting for more general shapes of the α 2 F function. Using symbolic regression and the SISSO framework, together with a database of artificially generated α 2 F functions and numerical solutions of the Eliashberg equations, we derive a formula for T c that performs as well as Allen–Dynes for low- T c superconductors and substantially better for higher- T c ones. This corrects the systematic underestimation of T c while reproducing the physical constraints originally outlined by Allen and Dynes. This equation should replace the Allen–Dynes formula for the prediction of higher-temperature superconductors.
Dynamic Properties of FeCrMnNi, a High Entropy Alloy
High entropy alloys (HEAs) represent a new class of materials that potentially possess an ideal combination of high strength and ductility and may be candidate materials that can withstand high amounts of damage before failure, especially under extreme conditions. To this end, a series of shock/spall loading experiments were conducted on an Fe/Cr/Mn/Ni HEA to assess not only its equation of state (EOS) but also its resistance to damage and failure. A single-stage light gas gun was used to conduct a series of flyer plate symmetric impact experiments coupled with recovery of the samples for post-mortem microstructural analysis to obtain fundamental understanding regarding the damage morphology. Our results show that the measured EOS for this alloy was in reasonable agreement with an estimated Hugoniot value. Additionally, the spall strength of these alloys was found to vary between 1.4 and 1.9 GPa within the same plate but was comparable to other iron-based conventional alloys. DOE/NV/03624--1059
Spectral properties and enhanced superconductivity in renormalized Migdal-Eliashberg theory
Migdal-Eliashberg theory describes the properties of the normal and superconducting states of electron-phonon-mediated superconductors based on a perturbative treatment of the electron-phonon interactions. It is necessary to include both electron and phonon self-energies self-consistently in Migdal-Eliashberg theory in order to match numerically exact results from determinantal quantum Monte Carlo in the adiabatic limit. Here in this work we provide a method to obtain the real-axis solutions of the Migdal-Eliashberg equations with electron and phonon self-energies calculated self-consistently. Our method avoids the typical challenge of computing cumbersome singular integrals on the real axis and is numerically stable and exhibits fast convergence. Analyzing the resulting real-frequency spectra and self-energies of the two-dimensional Holstein model, we find that self-consistently including the lowest-order correction to the phonon self-energy significantly affects the solution of the Migdal-Eliashberg equations. The calculation captures the broadness of the spectral function, renormalization of the phonon dispersion, enhanced effective electron-phonon coupling strength, minimal increase in the electron effective mass, and the enhancement of superconductivity which manifests as a superconducting ground state despite strong competition with charge-density-wave order. We discuss surprising differences in two common definitions of the electron-phonon coupling strength derived from the electron mass and the density of states, quantities which are accessible through experiments such as angle-resolved photoemission spectroscopy and electron tunneling. An approximate upper bound on 2Δ/T c for conventional superconductors mediated by retarded electron-phonon interactions is proposed.
A CONDUCTION-BASED HEAT PIPE MODEL FOR ANALYZING THE ENTIRE PROCESS OF LIQUID-METAL HEAT PIPE STARTUP
The Heat Pipe-cooled Microreactor (HPM) is one of the micro nuclear reactor designs under active study at the U.S. Idaho National Laboratory. Among the major concerns of HPM research is to understand the startup behavior of the heat pipe-cooled system associated with the startup of the liquid-metal heat pipes initialing from frozen state. The startup of liquid-metal heat pipes typically involves a number of nonlinear mass and heat transport processes including the phase change from solid to liquid and vapor. Hence, it is still a huge challenge to simulate the liquid-metal heat pipe startup using conventional CFD methods and software. The major difficulties of numerical CFD modeling come from the phase-change process, multiphase interaction, microporous wick flow, and compressible gas dynamics that occur during startup of the liquid-metal heat pipes. This paper proposes a simplified conduction-based method to provide practical insights into the entire startup process of the liquid-metal heat pipes while mitigating the challenges of addressing all the complex physics. We discuss the theoretical basis and modeling assumptions to analyze the liquid-metal heat pipe startup from frozen state based solely on heat-conduction equations. Then, the proposed model is implemented into the commercial CFD software to verify the model performance. The model prediction results are discussed via the comparison with the experimental data obtained from sodium heat-pipe startup experiments.
Adaptive Data-Driven Deep-Learning Surrogate Model for Frontal Polymerization in Dicyclopentadiene
Frontal polymerization (FP) is a self-sustaining curing process that enables rapid and energy-efficient manufacturing of thermoset polymers and composites. Computational methods conventionally used to simulate the FP process are time-consuming, and repeating simulations are required for sensitivity analysis, uncertainty quantification, or optimization of the manufacturing process. Here, in this work, we develop an adaptive surrogate deep-learning model for FP of dicyclopentadiene (DCPD), which predicts the evolution of temperature and degree of cure orders of magnitude faster than the finite-element method (FEM). The adaptive algorithm provides a strategy to select training samples efficiently and save computational costs by reducing the redundancy of FEM-based training samples. The adaptive algorithm calculates the residual error of the FP governing equations using automatic differentiation of the deep neural network. A probability density function expressed in terms of the residual error is used to select training samples from the Sobol sequence space. The temperature and degree of cure evolution of each training sample are obtained by a 2D FEM simulation. The adaptive method is more efficient and has a better prediction accuracy than the random sampling method. With the well-trained surrogate neural network, the FP characteristics (front speed, shape, and temperature) can be extracted quickly from the predicted temperature and degree-of-cure fields.
Advancing simulations of coupled electron and phonon nonequilibrium dynamics using adaptive and multirate time integration
Electronic structure calculations in the time domain provide a deeper understanding of nonequilibrium dynamics in materials. The real-time Boltzmann equation (rt-BTE), used in conjunction with accurate interactions computed from first principles, has enabled reliable predictions of coupled electron and lattice dynamics. However, the timescales and system sizes accessible with this approach are still limited, with two main challenges being the different timescales of electron and phonon interactions and the cost of computing collision integrals. As a result, only a few examples of these calculations exist, mainly for two-dimensional (2D) materials. Here we leverage adaptive and multirate time integration methods to achieve a major step forward in solving the coupled rt-BTEs for electrons and phonons. Relative to conventional (non-adaptive) time-stepping, our approach achieves a 10x speedup for a target accuracy, or greater accuracy by 3–6 orders of magnitude for the same computational cost, enabling efficient calculations in both 2D and bulk materials. This efficiency is showcased by computing the coupled electron and lattice dynamics in graphene up to ~100 ps, as well as modeling ultrafast lattice dynamics and thermal diffuse scattering maps in bulk materials (silicon and gallium arsenide). In addition to improved efficiency, our adaptive method can resolve the characteristic rates of different physical processes, thus naturally bridging different timescales. This enables simulations of longer timescales and provides a framework for modeling multiscale dynamics of coupled degrees of freedom in matter. Our work opens new opportunities for quantitative studies of nonequilibrium physics in materials, including driven lattice dynamics with phonons coupled to electrons, spin, and other degrees of freedom.
Analytical model of plasma response to external magnetic perturbation in absence of no-slip condition
Recent simulation and experimental results suggest that the magnetic island and flow on a resonant surface often do not satisfy the “no-slip” condition in the steady state. A new theoretical model on a nonlinear plasma response to external magnetic perturbation in the absence of a no-slip condition is proposed. The model is composed of the equations for the evolution of both the width and phase of the magnetic island due to forced reconnection driven by the external magnetic perturbation, and the force-balance equation for the plasma flow. When the island width is much less than the resistive layer width, the island growth is governed by the linear Hahm–Kulsrud–Taylor solution in the presence of time-dependent plasma flow. In the other regime when the island width is much larger than the resistive layer width, the evolution of both the island width and phase can be described using the Rutherford theory. The island solution is used to construct the quasilinear electromagnetic force, which, together with the viscous one, contributes to the nonlinear variation in plasma flow. The no-slip condition assumed in the conventional error field theory is not imposed here, where the island oscillation frequency depends on but does not necessarily equal to the plasma flow frequency at the rational surface.
Data Science Enabled Enabled Discovery of Superconductors (Final Progress Report)
This Final Technical Report describes efforts by 4 PIs at the University of Florida (Peter Hirschfeld, Richard Hennig, Greg Stewart and James Hamlin), over the period September 2019-August 2023, to use data science and machine learning techniques to discover new conventional superconductors. The PIs constructed a discovery loop with two theorists and two experimentalists to: develop algorithms to machine learn descriptors correlating strongly with the critical temperature Tc (PI's Peter Hirschfeld, UF Physics and Richard Hennig, UF Materials Science and En), synthesize and measure properties of promising materials, and feed back the knowledge gained into the prediction algorithm. This work was motivated by the theoretical prediction and experimental discovery of high-pressure, high-pressure hydride superconductors, and to find ways to recreate the high critical temperatures in these systems at ambient pressure. Highlights from the grant include: 1) a new equation for Tc in terms of moments of the electron-phonon spectral function, improving on the so-called Allen-Dynes equation (1975); 2) study of the metastable A15 superconductor Nb3Si, formed under explosive compression at ~1000GPa to determine the kinetic barrier to the ground state structure; 3) the development of ultra-fast machine-learned atomic potentials for molecular dynamics, and 4) the discovery of superconductivity at 19K in WB2 arising from metastable defect structures in the crystal.
PCET‐Driven Reactivity of Neptunyl(VI) Yields Oxo‐Bridged Np(V) and Np(IV) Species
Two unconventional polynuclear complexes of neptunium (Np) featuring mono-mathematical equation -oxo motifs have been accessed by proton-coupled electron transfer (PCET) reactivity involving the dissolution of neptunyl(VI) diacetate dihydrate (NpO 2 (OAc) 2 (H 2 O) 2 ∙ HOAc) in methanol followed by addition of a pentadentate Schiff-base ligand. One complex is a mixed-valent [Np V ,Np IV , Np V ] trimer with two bridging mathematical equation μ 2 -oxos and the other is a [Np V , Np V ] dimer featuring a single mathematical equation μ 2 -oxo. In both complexes the outer Np centers are capped with terminal oxo ligands. Spectroscopic and spectrokinetic studies aimed at elucidating mechanistic details of complex formation in this system show that intermediate multinuclear [Np V O 2 (OAc)] n species form prior to metal chelation by the ligand; electrolysis experiments demonstrate that production of Np(V) gives rise to asynchronous proton transfer that does not occur otherwise (in the Np(VI) state) as well as condensation with loss of H 2 O and formation of the polynuclear complexes. We attribute the oxo-deficient nature of these products, with respect to conventional actinyl ([AnO 2 ] m+ ) species, to the reduction/condensation reaction sequence of PCET.
Stress intensity factor models using mechanics-guided decomposition and symbolic regression
The finite element method can be used to compute accurate stress intensity factors (SIFs) for cracks with complex geometries and boundary conditions. In contrast, handbook solutions act as surrogate SIF models that provide significantly faster evaluation times. However, the development of conventional surrogate SIF models relies on manual development based on low-order parameterizations. This limits surrogate model accuracy and generalizability. Here, in this paper, we develop a framework for the automated development of mechanics-guided handbook SIF solutions by using interpretable machine learning via genetic programming for symbolic regression (GPSR). Formalizing the mechanics-based approach of Raju and Newman, SIF training data is decomposed into multiple subsets. This decomposition enables parallel GPSR model development of subfunctions, each of which accounts for specific geometrical corrections with respect to a known analytical model. Using this mechanics-based approach with GPSR allows for equations to be learned with improved accuracy and reduced complexity relative to the Raju Newman equations while maintaining the inherent interpretability of mathematical expressions. In this paper, we present equations that match the complexity of the Raju Newman equations while having reduced error, as well as equations with similar errors and reduced complexity.
First-principles free energies by hybrid thermodynamic integration for phase equilibria and fission product solubility in molten salts
The overarching goal of this research project is quantitative microscopic modeling of phase diagrams and solubility of actinides and lanthanide-fission products in molten salts for molten salt reactors. Application of molten salts as coolants and liquid fuels for next-generation reactors requires quantitative understanding of phase diagrams, equations of state and transport properties of binary and ternary molten salt mixtures, as well as solubility of fuel, fission products and reactor-wall corrosion products, but these remain beyond the reach of current first-principles computational techniques. We will address this fundamental capability gap using a novel hybrid framework for thermodynamic integration to directly predict absolute free energies of solid and molten salts at any specified conditions of composition, temperature and pressure. This framework will combine a hierarchy of first-principles, molecular dynamics and continuum techniques to deliver first-principles-level accuracy at a fraction of the computational cost of conventional ab initio molecular dynamics approaches. In particular, this technique is expected to provide first-principles accuracy, while off-loading most of the computational work on much less expensive classical molecular dynamics and continuum techniques.
Continued Verification of MOOSE Structural Mechanics Tools for Modeling Core Bowing Phenomena in Fast Reactors
Under the U.S. Department of Energy Office of Nuclear Energy’s Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism whereby increased power (which leads to temperature and flux gradients) influences the core to bow into less reactive configurations when the restraint system is properly designed. The phenomenon includes a complex interplay of radiation transport, duct temperature calculations involving fluid flow and heat transfer, and thermo-mechanical responses to the induced temperature and flux gradients. Structural material properties are also important to determining inelastic response to longer term flux gradients which cause irradiation creep and swelling. While core bowing provides a strong negative reactivity feedback when the restraint system is designed properly, it also results in additional forces between assemblies which increase the loads required to extricate them during refueling or control rod movement. Therefore, the restraint system must be designed with these tradeoffs in mind. The first stage of the work, which commenced in FY21 and continues through FY22, assesses thermo-mechanical modeling tools for producing core bowing predictions consistent with conventional tools. The Multiphysics Object Oriented Simulation Environment (MOOSE) Tensor Mechanics and Contact Modules are employed. This status report describes work on additional thermo-mechanical benchmark verification problems with increased complexity from the examples demonstrated in FY21. Several benchmark verification examples were selected from the IAEA verification and validation report. These examples involve clusters of ducts representative of a sector of a hexagonal reactor core which bow into each other and cause contact and load pad elevations, as well as single ducts subjected to irradiation fields undergoing swelling and subsequent bowing. The MOOSE-based results were compared to both IAEA benchmark participants’ results, analytic equations as available, and NUBOW-3D, a beam model code developed by Argonne National Laboratory. In every case, the MOOSE results agreed with other simulations results, providing additional verification basis of the tools for this particular physics application.
Ensemble Learning Based Convex Approximation of Three-Phase Power Flow
Though the convex optimization has been widely used in power systems, it still cannot guarantee to yield a tight (accurate) solution to some problems. To mitigate this issue, this paper proposes an ensemble learning based convex approximation for alternating current (AC) power flow equations that differs from the existing convex relaxations. The proposed approach is based on three-phase quadratic power flow equations in rectangular coordinates. To develop this data-driven convex approximation of power flows, the polynomial regression (PR) is first deployed as a basic learner to fit convex relationships between the independent and dependent variables. Then, ensemble learning algorithms such as gradient boosting (GB) and bagging are introduced to combine learners to boost model performance. Based on the learned convex approximation of power flow, optimal power flow (OPF) is formulated as a convex quadratic programming problem. The simulation results on IEEE standard cases of both balanced and unbalanced systems show that, in the context of solving OPF, the proposed data-driven convex approximation outperforms the conventional semi-definite programming (SDP) relaxation in both accuracy and computational efficiency, especially in the cases that the conventional SDP relaxation fails
Driving rotamak currents with minimal power dissipation
The Rotamak is a proposed thermonuclear fusion device which employs rotating magnetic fields (RMF) to generate an azimuthal current to produce a field-reversed configuration. The efficiency of the currents that produce the field reversal by RMFs was debated some 40 years ago. The debate revolved around whether the currents would incur dissipation by the conventional Spitzer perpendicular resistivity, or whether some other relation between current and dissipation would be more appropriate. By employing an electron–ion pitch-angle scattering model, we find that the dissipation is non-Spitzer in nature. However, curiously, there appears to exist a regime where the power dissipated to maintain the current becomes vanishingly small.
Dielectric screening and vacancy formation for large neutral and charged Si n H m ( n > 1500 ) nanocrystals using real-space pseudopotentials
A commonly used procedure for computing the properties of defects in crystalline materials is to consider a large supercell that includes the defect of interest. This is a straightforward technique as standard energy band codes can be used for such computations. For neutral defects, the only impediment of such an approach is to avoid defect-defect interactions between adjoining cells. However, this procedure can be complex if the defect of interest is charged as the system at large contains Coulombic divergences. Moreover, some have recently argued that the conventional definition of formation energies for charged defects cannot be reconciled with statistical mechanics. Here, we focus on an alternative approach. We consider large nanocrystals wherein a charged defect can be placed. Since the system is confined, a charged defect within the nanocrystal does not result in a Coulombic divergence. The chief impediment is computational, i.e., while no defect-defect or Coulombic divergences are present, the nanocrystal must be sufficiently large to allow the system to properly replicate a bulklike configuration. With the development of new algorithms and hardware advances, computations for systems of sufficient size to address this issue are feasible. In particular, we solve the Kohn-Sham equation in real space using pseudopotential-density-functional theory for large silicon nanocrystals, which contain thousands of atoms. Further, we focus on (i) the screening of a point charge and (ii) the formation of a charged vacancy in hydrogen-terminated silicon nanocrystals. This approach allows us to examine the role of quantum confinement in addition to exploring the bulk limit. Comparisons to other methods confirm the viability of this approach.