Search NASA⌕ Search

SEARCH · Search NASA

Results for “free energy perturbation”

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 19 records

Statistical mechanics of light elements at high pressure. VII - A perturbative free energy for arbitrary mixtures of H and He

A model free energy is presented which accurately represents results from 45 high-precision Monte Carlo calculations of the thermodynamics of hydrogen-helium mixtures at pressures of astrophysical and planetophysical interest. The free energy is calculated using free-electron perturbation theory (dielectric function theory), and is an extension of the expression given in an earlier paper in this series. However, it fits the Monte Carlo results more accurately, and is valid for the full range of compositions from pure hydrogen to pure helium. Using the new free energy, the phase diagram of mixtures of liquid metallic hydrogen and helium is calculated and compared with earlier results. Sample results for mixing volumes are also presented, and the new free energy expression is used to compute a theoretical Jovian adiabat and compare the adiabat with results from three-dimensional Thomas-Fermi-Dirac theory. The present theory gives slightly higher densities at pressures of about 10 megabars.

Hubbard, W. B.↗

Nongyrotropy as a source of instability and mode coupling

Nongyrotropic particle populations can bring about linear mode coupling in homogeneous media among the three eigen-modes of parallel propagation in gyrotropic magnetoplasmas. These interactions stimulate, in general, wave activity that does not occur in corresponding (random gyrophase) gyrotropic ambients. Solutions of the dispersion equation illustrate that simple introduction of gyrophase organization can (1) excite electrostatic (and electromagnetic) perturbations in media whose free energy sources are solely electromagnetic, and (3) drive hybrid (both electrostatic and electromagnetic) wave growth in thoroughly stable Maxwellian plasmas.

Brinca, A. L.↗

Thermophysical properties of simple liquid metals: A brief review of theory

In this paper, we review the current theory of the thermophysical properties of simple liquid metals. The emphasis is on thermodynamic properties, but we also briefly discuss the nonequilibrium properties of liquid metals. We begin by defining a 'simple liquid metal' as one in which the valence electrons interact only weakly with the ionic cores, so that the interaction can be treated by perturbation theory. We then write down the equilibrium Hamiltonian of a liquid metal as a sum of five terms: the bare ion-ion interaction, the electron-electron interaction, the bare electron-ion interaction, and the kinetic energies of electrons and ions. Since the electron-ion interaction can be treated by perturbation, the electronic part contributes in two ways to the Helmholtz free energy: it gives a density-dependent term which is independent of the arrangement of ions, and it acts to screen the ion-ion interaction, giving rise to effective ion-ion pair potentials which are density-dependent, in general. After sketching the form of a typical pair potential, we briefly enumerate some methods for calculating the ionic distribution function and hence the Helmholtz free energy of the liquid: monte Carlo simulations, molecular dynamics simulations, and thermodynamic perturbation theory. The final result is a general expression for the Helmholtz free energy of the liquid metal. It can be used to calculate a wide range of thermodynamic properties of simple metal liquids, which we enumerate. They include not only a range of thermodynamic coefficients of both metals and alloys, but also many aspects of the phase diagram, including freezing curves of pure elements and phase diagrams of liquid alloys (including liquidus and solidus curves). We briefly mention some key discoveries resulting from previous applications of this method, and point out that the same methods work for other materials not normally considered to be liquid metals (such as colloidal suspensions, in which the suspended microspheres behave like ions screened by the salt solution in which they are suspended). We conclude with a brief discussion of some non-equilibrium (i.e., transport) properties which can be treated by an extension of these methods. These include electrical resistivity, thermal conductivity, viscosity, atomic self-diffusion coefficients, concentration diffusion coefficients in alloys, surface tension and thermal emissivity. Finally, we briefly mention two methods by which the theory might be extended to non-simple liquid metals: these are empirical techniques (i.e., empirical two- and three-body potentials), and numerical many-body approaches. Both may be potentially applicable to extremely complex systems, such as nonstoichiometric liquid semiconductor alloys.

Stroud, David↗

Good Practices in Free-energy Calculations

As access to computational resources continues to increase, free-energy calculations have emerged as a powerful tool that can play a predictive role in drug design. Yet, in a number of instances, the reliability of these calculations can be improved significantly if a number of precepts, or good practices are followed. For the most part, the theory upon which these good practices rely has been known for many years, but often overlooked, or simply ignored. In other cases, the theoretical developments are too recent for their potential to be fully grasped and merged into popular platforms for the computation of free-energy differences. The current best practices for carrying out free-energy calculations will be reviewed demonstrating that, at little to no additional cost, free-energy estimates could be markedly improved and bounded by meaningful error estimates. In energy perturbation and nonequilibrium work methods, monitoring the probability distributions that underlie the transformation between the states of interest, performing the calculation bidirectionally, stratifying the reaction pathway and choosing the most appropriate paradigms and algorithms for transforming between states offer significant gains in both accuracy and precision. In thermodynamic integration and probability distribution (histogramming) methods, properly designed adaptive techniques yield nearly uniform sampling of the relevant degrees of freedom and, by doing so, could markedly improve efficiency and accuracy of free energy calculations without incurring any additional computational expense.

Pohorille, Andrew↗

The equation of state for stellar envelopes. II - Algorithm and selected results

A free-energy-minimization method for computing the dissociation and ionization equilibrium of a multicomponent gas is discussed. The adopted free energy includes terms representing the translational free energy of atoms, ions, and molecules; the internal free energy of particles with excited states; the free energy of a partially degenerate electron gas; and the configurational free energy from shielded Coulomb interactions among charged particles. Internal partition functions are truncated using an occupation probability formalism that accounts for perturbations of bound states by both neutral and charged perturbers. The entire theory is analytical and differentiable to all orders, so it is possible to write explicit analytical formulas for all derivatives required in a Newton-Raphson iteration; these are presented to facilitate future work. Some representative results for both Saha and free-energy-minimization equilibria are presented for a hydrogen-helium plasma with N(He)/N(H) = 0.10. These illustrate nicely the phenomena of pressure dissociation and ionization, and also demonstrate vividly the importance of choosing a reliable cutoff procedure for internal partition functions.

Mihalas, Dimitri↗

Kinetic instability associated with spherical shell distribution of cometary pickup ions

The stability of ion shell distribution is an interesting and important issue in view of the fact that it is stable under perturbations propagating parallel to the ambient magnetic field, despite large amount of free energy it contains. It was found, however, that if the shell distribution of ions has a sufficient amount of drift speed with respect to the solar wind hydromagnetic waves, it can become unstable. The obliquely propagating hydromagnetic instability associated with the ion shell distribution was recently discovered. This letter presents the unified stability analysis of the ion shell distribution, under what condition the drift-induced instability of the type first discussed by Gary and Sinha (1989), and the oblique hydromagnetic instability discussed by Wu and Yoon (1990) are excited, and how the two types of instabilities are related.

Yoon, Peter H.↗

Ion Acceleration and Transport in Solar Flares

The purpose of the work proposed for this grant was to develop a promising model for ion acceleration in impulsive solar flares. Solar flares are among the most energetic and interesting phenomena in the solar system, releasing up to 10(exp 32) ergs of energy over timescales ranging from a few tens of seconds to a few tens of minutes. Much of this energy appears as energetic electrons and ions, which produce a wide range of observable radiations. These radiations, in turn, are valuable diagnostics of the acceleration mechanism, the identification of which is the fundamental goal of solar flare research. The specific mechanism we proposed to investigate was based on cascading Alfven waves, the essence of which was as follows: During the primary flare energy release, it is widely believed that magnetic free energy is made available through the large-scale restructuring of the flare magnetic field. Any perturbation of a magnetic field will lead to the formation of MagnetoHydroDynamic (MHD) waves of wavelength comparable to the initial scale of the perturbation. Since the scalesize of a flare energy release region will likely be 10(exp 8)-10(exp 9) cm, the MHD waves will be of very long wavelength. However, it is well known that wave steepening will lead to a cascade of wave energy to smaller wavelengths. Now, MHD waves consist of two specific modes-the Alfven wave and the fast mode wave, and it is the Alfven wave which can interact with the ambient ions and accelerate them via cyclotron resonance. As the Alfven waves cascade to smaller wavenumbers, they can resonate with ions of progressively lower energy, until they eventually (actually, this is less than approx. 1 s) can resonate with ions in the thermal distribution. These ions are then energized out of the thermal background and, since lower-frequency waves are already present as a result of the cascading, to relativistic energies. Hence, cascading Alfven waves naturally accelerate ions from thermal to relativistic energies in one step with one basic mechanism.

Miller, James A.↗

Developing disturbances in shear

The origin and growth to moderate amplitude of disturbances in shear flow has been traditionally ascribed to the linear modal instability of the flow. Recent work on initial value problems has suggested that nonmodal growth of perturbations may be of equal and perhaps greater importance, at least in cases of rapid development. Examples of robust growth in model problems which support no instabilities and in baroclinic flows with realistic Ekman damping for which the exponential modes have zero or negative growth rates have been shown. Such examples have focused attention on the perturbations which are configured to tap the energy of mean flows. Here a critical examination of these favorably configured perturbations is given, making use of the simple constant free shear barotropic model which allows construction of exact two-dimensional isolated wave packet solutions.

Farrell, Brian↗

Impacts of free-floating objects: Unique space station experiments

The transfer of momentum and kinetic energy between planetary bodies forms the basis for wide ranging problems in planetary science ranging from the collective long term effects of minor perturbations to the catastrophic singular effect of a major collision. Although the collisional transfer of momentum and energy was discussed over the last two decades, major issues remain that largely reflect current limitations in Earth based experimental conditions and 3-D numerical codes. Two examples with potential applications in a Space Station laboratory, are presented: asteroid spin rates and orientations, and planetary disruption/spin rates. Asteroid spin rate and orientation experiments are needed wherein free floating nonspining and spining objects of varying strength, porosity, and volatility are impacted at varying velocities and angles. A space station platform also could provide an opportunity to test important facets of planetary disruption/spin rate models by allowing freely suspended spherical targets of varying viscosities, internal density gradients, and spin rates.

Schultz, P. H.↗

Dynamic magnetic reconnection at an X-type neutral point

The relaxation of a two-dimensional 'X-type' neutral point magnetic field disturbed from equilibrium is considered. Perturbations are shown to possess well-defined azimuthal modes which allow an exact determination of the magnetic annihilation rate. Free magnetic energy is dissipated by oscillatory reconnection which couples resistive diffusion at the neutral point to global advection of the outer field. The decay of azimuthally symmetric (m = 0) modes - the only modes associated with topological reconnection - is limited by the dissipation time scale of the 'fundamental' (n = 0) mode with no radial nodes. This mode decays over typically 100 Alfven times. An analytic treatment shows that the oscillation and decay time scales couple according to a given law.

Craig, I. J. D.↗

First Satellite-detected Perturbations of Outgoing Longwave Radiation Associated with Blowing Snow Events over Antarctica

We present the first satellite-detected perturbations of the outgoing longwave radiation (OLR) associated with blowing snow events over the Antarctic ice sheet using data from Cloud-Aerosol Lidar with Orthogonal Polarization and Clouds and the Earth's Radiant Energy System. Significant cloud-free OLR differences are observed between the clear and blowing snow sky, with the sign andmagnitude depending on season and time of the day. During nighttime, OLRs are usually larger when blowing snow is present; the average difference in OLRs between without and with blowing snow over the East Antarctic Ice Sheet is about 5.2 W/m2 for the winter months of 2009. During daytime, in contrast, the OLR perturbation is usually smaller or even has the opposite sign. The observed seasonal variations and day-night differences in the OLR perturbation are consistent with theoretical calculations of the influence of blowing snow on OLR. Detailed atmospheric profiles are needed to quantify the radiative effect of blowing snow from the satellite observations.

outgoing longwave radiation (OLR)↗

Selective Two-Photon Absorptive Resonance Femtosecond-Laser Electronic-Excitation Tagging (STARFLEET) Velocimetry in Flow and Combustion Diagnostics

Selective Two-Photon Absorptive Resonance Femtosecond-Laser Electronic-Excitation Tagging (STARFLEET), a non-seeded ultrafast-laser-based velocimetry technique, is demonstrated in reactive and non-reactive flows. STARFLEET is pumped via a two-photon resonance in N2 using 202.25-nm 100-fs light. STARFLEET greatly reduces the per-pulse energy required (30 μJ/pulse) to generate the signature FLEET emission compared to the conventional FLEET technique (1.1 mJ/pulse). This reduction in laser energy results in less energy deposited in the flow, which allows for reduced flow perturbations (reactive and non-reactive), increased thermometric accuracy, and less severe damage to materials. Velocity measurements conducted in a free jet of N2 and in a premixed flame show good agreement with theoretical velocities and further demonstrate the significantly less-intrusive nature of STARFLEET.

Jiang, Naibo↗

Impacts of free-floating objects: Unique Space Station experiments

The transfer of momentum and kinetic energy between planetary bodies forms the basis for wide-ranging problems in planetary science ranging from the collective long-term effects of minor perturbations to the catastrophic singular effect of a major collision. In the former case, the evolution of asteroid spin rates and orientations and planetary rotation rates are cited. In the latter case, the catastrophic angular momenta and the near-global disruption of partially molten planets are included. Although the collisional transfer of momentum and energy were discussed over the last two decades, major issues remain that largely reflect current limitations in earth-based experimental conditions and 3-D numerical codes. Two examples with potential applications in a Space Station laboratory are presented.

Schultz, Peter H.↗

Interaction of a Boundary Layer with a Turbulent Wake

The objective of this grant was to study the transition mechanisms on a flat-plate boundary layer interacting with the wake of a bluff body. This is a simplified configuration presented and designed to exemplify the phenomena that occur in multi-element airfoils, in which the wake of an upstream element impinges on a downstream one. Some experimental data is available for this configuration at various Reynolds numbers. The first task carried out was the implementation and validation of the immersed-boundary method. This was achieved by performing calculations of the flow over a cylinder at low and moderate Reynolds numbers. The low-Reynolds number results are discussed, which is enclosed as Appendix A. The high-Reynolds number results are presented in a paper in preparation for the Journal of Fluid Mechanics. We performed calculations of the wake-boundary-layer interaction at two Reynolds numbers, Re approximately equal to 385 and 1155. The first case is discussed and a comparison of the two calculations is reported. The simulations indicate that at the lower Reynolds number the boundary layer is buffeted by the unsteady Karman vortex street shed by the cylinder. This is shown: long streaky structures appear in the boundary layer in correspondence of the three-dimensionalities in the rollers. The fluctuations, however, cannot be self-sustained due to the low Reynolds-number, and the flow does not reach a turbulent state within the computational domain. In contrast, in the higher Reynolds-number case, boundary-layer fluctuations persist after the wake has decayed (due, in part, to the higher values of the local Reynolds number Re achieved in this case); some evidence could be observed that a self-sustaining turbulence generation cycle was beginning to be established. A third simulation was subsequently carried out at a higher Reynolds number, Re=3900. This calculation gave results similar to those of the Re=l155 case. Turbulence was established at fairly low Reynolds number, as a consequence of the high level of the free-stream perturbation. An instantaneous flow visualization for that case is shown. A detailed examination of flow statistics in the transitional and turbulent regions, including the evolution of the turbulent kinetic energy (TKE) budget and frequency spectra showed the formation and evolution of turbulent spots characteristic of the bypass transition mechanism. It was also observed that the turbulent eddies achieved an equilibrium, fully developed turbulent states first, as evidenced by the early agreement achieved by the terms in the TKE budget with those observed in turbulent flows. Once a turbulent Reynolds stress profile had been established, the velocity profile began to resemble a turbulent one, first in the inner region and later in the outer region of the wall layer. An extensive comparison of the three cases, including budgets, mean velocity and Reynolds stress profiles and flow visualization, is included. The results obtained are also presented.

Piomelli, Ugo↗

Laser Assisted Free-Free Transition in Electron - Atom Collision

Free-free transition is studied for electron-Hydrogen atom system in ground state at very low incident energies in presence of an external homogeneous, monochromatic and linearly polarized laser field. The incident electron is considered to be dressed by the laser in a non perturbative manner by choosing the Volkov solutions in both the channels. The space part of the scattering wave function for the electron is solved numerically by taking into account the effect of electron exchange, short range as well as of the long range interactions. Laser assisted differential as well as elastic total cross sections are calculated for single photon absorption/emission in the soft photon limit, the laser intensity being much less than the atomic field intensity. A strong suppression is noted in the laser assisted cross sections as compared to the field free situations. Significant difference is noted in the singlet and the triplet cross sections.

Sinha, C.↗

Arctic Cloud Response to a Perturbation in Sea Ice Concentration: The North Water Polynya

Surface and atmosphere energy exchanges play an important role in the Arctic climate system by influencing the lower atmospheric stability and humidity, sea ice melt and growth, and surface temperature. Sea ice significantly alters the character of these energy exchanges relative to ice-free ocean. The observed decline in Arctic sea ice since 1979 motivates questions related to the evolving role of surface-atmosphere coupling and potential feedbacks on the Arctic system. Due to the strong wintertime cloud warming effect, a critical question concerns the potential response of low clouds to Arctic sea ice decline. Previous approaches relied on interannual variability to investigate the cloud response to sea ice decline. However, the covariation between atmospheric conditions and sea ice makes it difficult to define an observational control when using interannual variability. To circumvent this difficulty, we exploit the recurring North Water polynya, an episodic opening in the northern Baffin Bay sea ice, as a natural laboratory to isolate the cloud response to a rapid, near-step perturbation in sea ice. Our results show that during the event, (a) low-cloud cover is 10%–33% larger over the polynya than nearby sea ice, (b) cloud liquid water content is up to 400% larger over the polynya than nearby sea ice, and (c) the surface cloud radiative effect is 18 W/sq. m larger over the polynya than nearby sea ice. Our results provide evidence that the low-cloud response during a polynya is a positive feedback lengthening the event.

Emily E. Monroe↗

Observations of Pc 1-2 waves in the outer magnetosphere

High time resolution measurements made inside the magnetopause from L = 7 to L = 14 by the flux gate magnetometer aboard the Ogo 5 satellite, which demonstrate that the Pc 1-2 magnetic field oscillations are present in the outer magnetosphere, are examined. The solar wind dynamic pressure and the magnetospheric cold plasma were enhanced, and the IMF had a southward component 2 hours before the events, but there was no significant correlation between the time of the events and either storm recovery phase or substorm onset. The magnetic field perturbations in each event were transverse to the ambient field with amplitudes from 2 to 8 gamma, and most of the transverse perturbations showed left-handed polarization. Excursions into the magnetosheath during two events revealed a disturbed magnetic field with significant power in frequencies of minimum 0.1 Hz, and although the magnetosheath was a source of free energy for waves at Pc 1-2 frequencies, the data within the magnetosphere suggest that the observed pulsations were ion cyclotron waves generated at the geomagnetic equator at large radial distances.

Kaye, S. M.↗