Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “potential energy”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Nonaxisymmetric secular instabilities driven by star/disk coupling

We determine conditions for the onset of nonaxisymmetric secular instabilities in polytropes with a wide range of angular momentum distributions using Lagrangian techniques, and then calculate the growth rate of such instabilities when driven by the coupling of the perturbed star to a circumstellar disk. We use Langrangian displacement vectors with azimuthal dependence proportional to exp (im phi), where m is an integer and phi is the azimuthal coordinate. The onset of secular instability in terms of the quantity T/absolute value of W, the ratio of rotational kinetic energy to gravitational potential energy, is affected by both the compressibility and the angular momentum distribution of the polytrope. The largest effects occur when the angular momentum distribution is varied. For polytropic index n = 3/2, the onset of secular instability for the m = 2 mode (the bar mode), as determined by its neutral point, shifts from T/absolute value of W = 0.141 to 0.093, while the m = 5 mode neutral point shifts from T/absolute value of W = 0.088 to 0.031 over the range of angular momentum distributions we consider. The smallest critical T/absolute value of W-values occur for the angular momentum distributions which are the most peaked toward the equator. For the angular momentum distribution of a Maclaurin spheroid, as the polytropic index n is increased from 3/2 to 5/2, the neutral point for m = 2 shifts from T/absolute value of W = 0.141 to 0.144 and the netural point for m = 5 shifts from T/absolute value of W = 0.069 to 0.078. The netural points for m = 2 and 5 for the Maclaurin sequence (n = 0) are 0.137 and 0.0629, respectively. As the angular momentum distribution becomes more peaked toward the equatorial radius of the polytropes, the critical T/absolute value of W-values generally become less sensitive to the compressibility of the polytrope. Star/disk coupling can drive the secular instability in systems where the star is surrounded by a massive disk and, if the instability can grow to moderate amplitude, then the coupling can transport significant amounts of angular momentum from the star into the circumstellar disk. We find, for the particular case of rotating protostars during the accretion phase, that the instability growth time can be shorter than the accretion time. Further, if the instability can grow to amplitudes on the order of several percent, the star/disk coupling can remove angular momentum from the forming star faster than it is added by accretion.

Imamura, James, N.↗
NASA NTRS1996

A Search for Star Formation in the Translucent Cloud MBM 40

The star formation status of the translucent high-latitude molecular cloud, MBM 40, is explored through analysis of radio, infrared, optical, and X-ray data. With a peak visual extinction of 1 to 2 mag, MBM 40 is an example of a high-latitude cloud near the diffuse/translucent demarcation. However, unlike most translucent clouds, MBM 40 exhibits a compact morphology and a kinetic energy-to gravitational potential energy ratio near unity. Our radio data, encompassing the CO (J = 1-0), CS (J = 2-1), and H2CO 1(sub 11-1(sub 10), spectral line transitions, reveal that the cloud contains a ridge of molecular gas with n greater than or equal to 10(exp 3)/ cc. In addition, the molecular data, together with IRAS data, indicate that the mass of MBM 40 is approx. 40 solar mass. In light of the ever-increasing number of recently formed stars far from any dense molecular clouds or cores, we searched the environs of MBM 40 for any trace of recent star formation. We used the ROSAT All-Sky Survey X-ray data and a ROSAT PSPC pointed observation toward MBM 40 to identify 33 stellar candidates with properties consistent with pre-main-sequence (PMS) stars. Follow-up optical spectroscopy of the candidates with V less than 15.5 was conducted with the 1.5 m Fred Lawrence Whipple Observatory telescope in order to identify signatures of T Tauri or pre-main- sequence stars (such as the Li 6708 A resonance line). Since none of our optically observed candidates display standard PMS signatures, we conclude that MBM 40 displays no evidence of recent or ongoing star formation. The absence of high-density molecular cores in the cloud and the relatively low column density compared to star-forming interstellar clouds may be the principal reasons that MBM 40 is devoid of star formation. More detailed comparison between this cloud and other, higher extinction translucent and dark clouds may elucidate the necessary initial conditions for the onset of low-mass star formation.

Magnani, Loris↗
DOE OSTI2020

An extended methodology for automated calculations of non-Boltzmann kinetic sequences: H + C 2 H 2 + X and combustion impact

It is generally assumed in phenomenological kinetic models that bimolecular reactions only occur between species whose rovibrational energy follows a Boltzmann (thermal) distribution. That is, any complexes initially formed in non-Boltzmann distributions are assumed to be thermalized by energy-transferring collisions prior to bimolecular reactions. Given the high mole fractions of reactive species, X, in combustion environments, reactive collisions of the complexes with X often occur on the same timescale as energy-transferring collisions – yielding sequences proceeding through non-Boltzmann intermediates across multiple potential energy surfaces. Recent studies have shown that such non-Boltzmann kinetic sequences can have substantial impact on the global reactivity in combustion systems. Simulations of these non-Boltzmann reaction sequences, which can be described in phenomenological kinetic models via chemically termolecular reactions, require that rovibrational excitation from one potential energy surface be carried over to the next. Here this paper presents an extended theoretical and computational methodology that couples multiple master equations and derives rate constants for phenomenological reactions describing the conversion of thermal reactants to thermal products for use in phenomenological kinetic schemes. The methodology is then implemented using in-house scripts for non-Boltzmann sequences involving C2H$^*_3$ + X (with X = O 2 , H, and OH) where C 2 H$^*_3$ is formed via H + C 2 H 2 association – which were identified as having strong potential for influencing combustion predictions in a recent study. The results reveal that non-Boltzmann reaction sequences for X = O 2 (the primary focus of this paper) significantly alters the total conversion rate from H + C 2 H 2 to products and product branching fractions from those of thermal sequential pathways. Furthermore, the present results demonstrate that non-Boltzmann reaction sequences have significant impact – as high as an order of magnitude – on predicted ignition delay times. Similarly, they yield significantly different dependence of ignition delay times with temperature and O 2 mole fraction – yielding signatures that are likely observable experimentally.

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

The electron-transfer intermediates of the oxygen evolution reaction (OER) as polarons by in situ spectroscopy

The conversion of diffusive forms of energy (electrical and light) into short, compact chemical bonds by catalytic reactions regularly involves moving a carrier from an environment that favors delocalization to one that favors localization. While delocalization lowers the energy of the carrier through its kinetic energy, localization creates a polarization around the carrier that traps it in a potential energy minimum. The trapped carrier and its local distortion—termed a polaron in solids—can play a role as a highly reactive intermediate within energy-storing catalytic reactions but is rarely discussed as such. Here, we present this perspective of the polaron as a catalytic intermediate through recent in situ and time-resolved spectroscopic investigations of photo-triggered electrochemical reactions at material surfaces. The focus is on hole-trapping at metal–oxygen bonds, denoted M–OH*, in the context of the oxygen evolution reaction (OER) from water. The potential energy surface for the hole-polaron defines the structural distortions from the periodic lattice and the resulting “active” site of catalysis. This perspective will highlight how current and future time-resolved, multi-modal probes can use spectroscopic signatures of M–OH* polarons to obtain kinetic and structural information on the individual reaction steps of OER. A particular motivation is to provide the background needed for eventually relating this information to relevant catalytic descriptors by free energies. Lastly, the formation of the O–O chemical bond from the consumption of M–OH*, required to release O 2 and store energy in H 2 , will be discussed as the next target for experimental investigations.

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

Stochastic evaluation of four-component relativistic second-order many-body perturbation energies: A potentially quadratic-scaling correlation method

A second-order many-body perturbation correction to the relativistic Dirac-Hartree-Fock energy is evaluated stochastically by integrating 13-dimensional products of four-component spinors and Coulomb potentials. The integration in the real space of electron coordinates is carried out by the Monte Carlo (MC) method with the Metropolis sampling, whereas the MC integration in the imaginary-time domain is performed by the inverse-CDF (cumulative distribution function) method. The computational cost to reach a given relative statistical error for spatially compact but heavy molecules is observed to be no worse than cubic and possibly quadratic with the number of electrons or basis functions. This is a vast improvement over the quintic scaling of the conventional, deterministic second-order many-body perturbation method. The algorithm is also easily and efficiently parallelized with demonstrated 92% strong scalability going from 64 to 4096 processors for a fixed job size.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Real-World Evaluation of National Energy Efficiency Potential of Cold Storage Evaporator Technology in the Context of Engine Start-Stop Systems

National concerns over energy consumption and emissions from the transportation sector have prompted regulatory agencies to implement aggressive fuel economy targets for light-duty vehicles through the U.S. National Highway Traffic Safety Administration/Environmental Protection Agency (EPA) Corporate Average Fuel Economy (CAFE) program. Automotive manufacturers have responded by bringing competitive technologies to market that maximize efficiency while meeting or exceeding consumer performance and comfort expectations. In a collaborative effort among Toyota Motor Corporation, Argonne National Laboratory (ANL), and the National Renewable Energy Laboratory (NREL), the real-world savings of one such technology is evaluated. A commercially available Toyota Highlander equipped with two-phase cold storage technology was tested at ANL’s chassis dynamometer testing facility. The cold storage technology maintains the thermal state of air-conditioning evaporators to enable longer and more frequent engine-off operation in vehicles equipped with start-stop functionality. Test results were analyzed and provided to NREL where a novel simulation framework was developed and calibrated to the test data. The vehicle model was then exercised over a large set of real-world drive cycle and ambient condition data to estimate national-level fuel economy benefits. Results indicate that the cold storage evaporator provided national fuel consumption reductions of 0.113% relative to a conventional evaporator in the same vehicle. In addition, when the cold storage evaporator engine stop/start was enabled for any temperature and the baseline was limited to the EPA menu, Start and Stop credit assumption of 27°C, a national fuel savings of 0.497% was found. Fuel savings resulted from a combination of extended engine-off duration during idle events and increased frequency of deceleration fuel cutoff, both enabled by the ability of the cold storage evaporator to maintain thermal state in situations where air conditioning is active.

DIRECT ENERGY CONVERSION↗
DOE OSTI2024

Distinguishing homolytic vs heterolytic bond dissociation of phenylsulfonium cations with localized active space methods

Modeling chemical reactions with quantum chemical methods is challenging when the electronic structure varies significantly throughout the reaction and when electronic excited states are involved. Multireference methods, such as complete active space self-consistent field (CASSCF), can handle these multiconfigurational situations. However, even if the size of the needed active space is affordable, in many cases, the active space does not change consistently from reactant to product, causing discontinuities in the potential energy surface. The localized active space SCF (LASSCF) is a cheaper alternative to CASSCF for strongly correlated systems with weakly correlated fragments. The method is used for the first time to study a chemical reaction, namely the bond dissociation of a mono-, di-, and triphenylsulfonium cation. LASSCF calculations generate smooth potential energy scans more easily than the corresponding, more computationally expensive CASSCF calculations while predicting similar bond dissociation energies. Furthermore, our calculations suggest a homolytic bond cleavage for di- and triphenylsulfonium and a heterolytic pathway for monophenylsulfonium.

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

Atmospheric energetics in regions of intense convective activity

Synoptic-scale budgets of kinetic and total potential energy are computed using 3- and 6-h data at nine times from NASA's fourth Atmospheric Variability Experiment (AVE IV). Two intense squall lines occurred during the period. Energy budgets for areas that enclose regions of intense convection are shown to have systematic changes that relate to the life cycles of the convection. Some of the synoptic-scale energy processes associated with the convection are found to be larger than those observed in the vicinity of mature cyclones. Volumes enclosing intense convection are found to have large values of cross-contour conversion of potential to kinetic energy and large horizontal export of kinetic energy. Although small net vertical transport of kinetic energy is observed, values at individual layers indicate large upward transport. Transfer of kinetic energy from grid to subgrid scales of motion occurs in the volumes. Latent heat release is large in the middle and upper troposphere and is thought to be the cause of the observed cyclic changes in the budget terms. Total potential energy is found to be imported horizontally in the lower half of the atmosphere, transported aloft, and then exported horizontally. Although local changes of kinetic energy and total potential energy are small, interaction between volumes enclosing convection with surrounding larger volumes is quite large.

Fuelberg, H. E.↗
NASA NTRSDate not supplied

REACTER: A Versatile Tool for Large-Scale Reactive Molecular Dynamics

Accurately describing reactive events over long length and time scales remains a grand challenge of computational materials science. REACTER is a general protocol for modeling chemical reactions using classical force fields, and is implemented in the popular molecular dynamics software LAMMPS. REACTER has a growing user base and has been used as a model-building tool for a variety of materials, including thermoplastics, thermosets, glassy materials and composites. The method has also been applied to accelerated modeling of reversible chemical reactions, such as the formation of electrochemical components for batteries. Recently, the REACTER protocol has received some major upgrades to enhance its ability to predict when reactions occur and to make it easier to use. Force field parameters can now be automatically assigned to newly created bonds, angles and other interactions. Advanced reaction constraints have also been added, including an Arrhenius constraint to enforce an effective activation energy, a root-mean-square-deviation option for complex geometrical constraints, and a constraint based on the total potential energy of the atoms involved in a reactive site. This potential energy constraint allows for the accurate reproduction of DFT-based tight-binding (DFTB3) predicted bond dissociation curves for mechanically induced bond breaking.

polymer simulations, molecular dynamics↗
DOE OSTI2021

Taylor dispersion in osmotically driven laminar flows in phloem

Sucrose is among the main products of photosynthesis that are deemed necessary for plant growth and survival. It is produced in the mesophyll cells of leaves and translocated to different parts of the plant through the phloem. Progress in understanding this transport process remains fraught with experimental difficulties, thereby prompting interest in theoretical approaches and laboratory studies. The Münch pressure and mass flow model is one of the accepted hypotheses describing the physics of sucrose transport in the phloem. It is based on osmosis creating an energy potential difference between the source and the sink. The flow responding to this energy potential is assumed laminar and described by the Hagen–Poiseuille equation. This study revisits such osmotically driven flows in tubes with membrane walls by including the effects of Taylor dispersion on mass transport. This effect has been overlooked in phloem flow studies. Taylor dispersion can increase the effective transport of solutes by increasing the apparent diffusion coefficient. It is shown that, in addition to the conventional diffusive correction derived for impermeable tube walls, a new adjustment to the mean advective terms arises because of osmotic effects. Because the molecular Schmidt number is very large for sucrose in water, the sucrose front speed and travel times have a direct dependence on the Péclet number for different ranges of the Münch number. Fially, this study establishes upper limits on expected Taylor dispersion enhancement of sucrose transport.

42 ENGINEERING↗
NASA NTRS1986

Effects of eddy initial conditions on nonlinear forcing of planetary scale waves by amplifying baroclinic eddies

The previous study of Young and Villere concerning growth of planetary scale waves forced by wave-wave interactions of amplifying intermediate scale baroclinic eddies is extended to investigate effects of different eddy initial conditions. A global, spectral, primitive equation model is used for the calculations. For every set of eddy initial conditions considered, growth rates of planetary modes are considerably greater than growth rates computed from linear instability theory for a fixed zonally independent basic state. However, values of growth rates ranged over a factor of 3 depending on the particular set of eddy initial conditions used. Nonlinear forcing of planetary modes via wave-wave coupling becomes more important than baroclinic growth on the basic state at small values of the intermediate-scale modal amplitudes. The relative importance of direct transfer of kinetic energy from intermediate scales of motion to a planetary mode, compared to baroclinic conversion of available potential energy to kinetic energy within that planetary mode, depends on the individual case. In all cases, however, the transfer of either kinetic or available potential energy to the planetary modes was accomplished principally by wave-wave transfer from intermediate scale eddies, rather than from the zonally averaged state. The zonal wavenumber 2 planetary mode was prominent in all solutions, even in those for which eddy initial conditions were such that a different planetary mode was selectively forced at the start. General characteristics of the structural evolution of the planetary wave components of total heat and momentum flux, and modal structures themselves, were relatively insensitive to variations in eddy initial conditions, even though quantitative details varied from case to case.

Young, Richard E.↗
NASA NTRS1999

Interaction Between Tropical Convection and its Embedding Environment: An Energetics Analysis of a 2-D Cloud Resolving Simulation

The phase relation between the perturbation kinetic energy (K') associated with the tropical convection and the horizontal-mean moist available potential energy (bar-P) associated with environmental conditions is investigated by an energetics analysis of a numerical experiment. This experiment is performed using a 2-D cloud resolving model forced by the TOGA-COARE derived vertical velocity. The imposed upward motion leads to a decrease of bar-P directly through the associated vertical advective cooling, and to an increase of K' directly through cloud related processes, feeding the convection. The maximum K' and its maximum growth rate lags and leads, respectively, the maximum imposed large-scale upward motion by about 1-2 hours, indicating that convection is phase locked with large-scale forcing. The dominant life cycle of the simulated convection is about 9 hours, whereas the time scales of the imposed large-scale forcing are longer than the diurnal cycle. In the convective events, maximum growth of K' leads maximum decay of the perturbation moist available potential energy (P') by about 3 hours through vertical heat transport by perturbation circulation, and perturbation cloud heating. Maximum decay of P' leads maximum decay of bar-P by about one hour through the perturbation radiative, processes, the horizontal-mean cloud heating, and the large-scale vertical advective cooling. Therefore, maximum gain of K' occurs about 4-5 hours before maximum decay of bar-P.

Li, Xiaofan↗
NASA NTRS2022

Washington, D.C. & Maryland Energy: Estimating Solar Potential Using NASA POWER Data to Inform Renewable Energy Policy for Washington, D.C.

In line with the Sustainable D.C. 2.0 plan to combat climate change, Washington D.C. aims to decrease its greenhouse gas emissions by 100% by 2050. As solar energy is a clean, renewable energy form, its integration into the region’s power grids lowers energy costs and incentivizes sustainable development. We partnered with the Washington D.C. Department of Energy & Environment (DOEE) to determine how urban areas surrounding D.C. can better be incorporated into decisions regarding renewable energy policy. The team used NASA’s Prediction of Worldwide Energy Resources (POWER) solar data and a Light Detection and Ranging (LiDAR) derived digital surface model, to estimate and visualize rooftop solar potential for Maryland’s Prince George’s and Montgomery counties. POWER provided solar irradiance data adjusted for tilt angle while the digital surface model contributed aspect and slope data. This methodology factored out areas that were unsuitable for solar panel installation while displaying areas that possess a high potential for energy return. The team found the total rooftop solar potential for the study area to be almost 32 million kW, which is equivalent to roughly 660 kW per building. The methodology used to generate the solar potential maps can be applied to other regions of the country seeking to efficiently utilize solar energy. The end users at the DOEE can use our resulting solar potential map and data table to effectively target buildings that have the highest potential to generate solar energy.

Edward Cronin↗
DOE OSTI2021

Global Optimization of Chemical Cluster Structures: Methods, Applications, and Challenges

Chemical clusters are relevant to many applications in catalysis, separations, materials, and energy sciences. Experimentally, the structure of clusters is difficult to determine, but it is very important in understanding their chemistry and properties. Computational methods can be used to examine cluster structure, however finding the most stable structure is not simple, particularly as the cluster size increases. Global optimization techniques have long been used to tackle the problem of the most stable structure, but such approaches would have to look for a global minimum, while sampling local minima over the whole potential energy surface as well. In this review, the state-of-the-art theory of global optimization theory is summarized. First, the definition, significance, relation to experiments, and a brief history of global optimization is presented. We then discuss, in more detail, three versatile global optimization methods: the basin hopping, the artificial bee colony algorithm, and the genetic algorithm. We close with some representative application examples of global optimization of clusters since 2016 and the challenges, open questions and opportunities in this field.

Global optimization, Chemical clusters, Artificial↗
NASA NTRS1994

A detailed evaluation of heating processes in the middle atmosphere

A fundamental problem in the study of the terrestrial middle atmosphere is to calculate accurately the local heating due to the absorption of solar radiation. Knowledge of the heat budget is essential to understanding the atmospheric thermal structure, atmospheric motions, atmospheric chemistry, and their coupling. The evaluation of heating rates is complicated (especially above the stratopause) by the fact that the heating is not a simple one-step process. That is, the absorbed solar energy does not all immediately appear as heat. Rather, substantial portions of the incident energy may appear as internal energy of excited photolysis products (e.g., O(1D) or O2(1 delta)) or as chemical potential energy of product species such as atomic oxygen. The ultimate disposition of the internal and chemical energy possessed by the photolysis products determines the efficiency and thus the rate at which the middle atmosphere is heated. In studies of the heat budget, it is also vitally important to consider transport of long lived chemical species such as atomic oxygen above approximately 80 km. In such cases, the chemical potential energy may be transported great distances (horizontally or vertically) before undergoing a reaction to release the heat. Atomic oxygen influences the heating not only by reactions with itself and with O2 but also by reactions with odd-hydrogen species, especially those involving OH (Mlynczak and Solomon, 1991a). Consequently, absorbed solar energy may finally by converted to heat a long time after and at a location far from the original deposition. The purpose of this paper is to examine the solar and chemical heating processes and to present parameterizations for the heating efficiencies readily applicable for use in numerical models and heat budget studies. In the next two sections the processes relevant to the heating efficiencies for ozone and molecular oxygen will be reviewed. In section 4 the processes for the exothermic reactions will be reviewed and parameterizations for the heating efficiencies for both the solar and chemical processes will be presented in Section 5.

Mlynczak, Martin↗
DOE OSTI2024

Full‐dimensional coupled‐channel statistical approach to atom‐triatom systems and applications to H/D + O 3 reaction

Abstract The statistical quantum model (SQM), which assumes that the reactivity is controlled by entrance/exit channel quantum capture probabilities, is well suited for chemical reactions with a long‐lived intermediate complex. In this work, a time‐independent coupled‐channel implementation of the SQM approach is developed for atom‐triatom systems in full dimensionality. As SQM treats the capture dynamics quantum mechanically, it is capable of handling quantum effects such as tunneling. A detailed study of the H/D + O 3 capture dynamics was performed by applying the newly developed SQM method on an accurate global potential energy surface. Agreement with previous ring polymer molecular dynamics (RPMD) results on the same potential energy surface is excellent except for very low temperatures. The SQM results are also in reasonably good agreement with available experimental rate coefficients. The strong H/D kinetic isotope effect underscores the dominant role of quantum tunneling under an entrance channel barrier at low temperatures.

Yang, Dongzheng↗
DOE OSTI2022

A benchmark dataset for Hydrogen Combustion

The generation of reference data for deep learning models is challenging for reactive systems, and more so for combustion reactions due to the extreme conditions that create radical species and alternative spin states during the combustion process. Here, we extend intrinsic reaction coordinate (IRC) calculations with ab initio MD simulations and normal mode displacement calculations to more extensively cover the potential energy surface for 19 reaction channels for hydrogen combustion. A total of ~290,000 potential energies and ~1,270,000 nuclear force vectors are evaluated with a high quality range-separated hybrid density functional, ωB97X-V, to construct the reference data set, including transition state ensembles, for the deep learning models to study hydrogen combustion reaction.

08 HYDROGEN↗
NASA NTRS1982

Calculation of SF6-/SF6 and Cl-/CFCl3 electron attachment cross sections in the energy range 0-100 meV

Electron attachment cross sections for the processes SF6-/SF6 and Cl-/CFCl3 are calculated in a local theory using a model in which diatomic-like potential energy curves for the normal modes are constructed from available spectroscopic data. Thermally populated vibrational and rotational levels are included. Good agreement is found with experimental cross sections in the energy range 5-100 meV for a particular choice of potential energy curve parameters.

Chutjian, A.↗