Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “energy simulation”

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 451 records · Page 25

Ultrarelativistic electromagnetic pulses in plasmas

The physical processes of a linearly polarized electromagnetic pulse of highly relativistic amplitude in an underdense plasma accelerating particles to very high energies are studied through computer simulation. An electron-positron plasma is considered first. The maximum momenta achieved scale as the square of the wave amplitude. This acceleration stops when the bulk of the wave energy is converted to particle energy. The pulse leaves behind as a wake a vacuum region whose length scales as the amplitude of the wave. The results can be explained in terms of a snow plow or piston-like action of the radiation on the plasma. When a mass ratio other than unity is chosen and electrostatic effects begin to play a role, first the ion energy increases faster than the electron energy and then the electron energy catches up later, eventually reaching the same value.

Ashour-Abdalla, M.↗

ECUT (Energy Conversion and Utilization Technologies Program). Biocatalysis Project

Presented are the FY 1985 accomplishments, activities, and planned research efforts of the Biocatalysis Project of the U.S. Department of Energy, Energy Conversion and Utilization Technologies (ECUT) Program. The Project's technical activities were organized as follows: In the Molecular Modeling and Applied Genetics work element, research focused on (1) modeling and simulation studies to establish the physiological basis of high temperature tolerance in a selected enzyme and the catalytic mechanisms of three species of another enzyme, and (2) determining the degree of plasmid amplification and stability of several DNA bacterial strains. In the Bioprocess Engineering work element, research focused on (1) studies of plasmid propagation and the generation of models, (2) developing methods for preparing immobilized biocatalyst beads, and (3) developing an enzyme encapsulation method. In the Process Design and Analysis work element, research focused on (1) further refinement of a test case simulation of the economics and energy efficiency of alternative biocatalyzed production processes, (2) developing a candidate bioprocess to determine the potential for reduced energy consumption and facility/operating costs, and (3) a techno-economic assessment of potential advancements in microbial ammonia production.

Source record↗

Development of a Simulation Capability for the Space Station Active Rack Isolation System

To realize quality microgravity science on the International Space Station, many microgravity facilities will utilize the Active Rack Isolation System (ARIS). Simulation capabilities for ARIS will be needed to predict the microgravity environment. This paper discusses the development of a simulation model for use in predicting the performance of the ARIS in attenuating disturbances with frequency content between 0.01 Hz and 10 Hz. The derivation of the model utilizes an energy-based approach. The complete simulation includes the dynamic model of the ISPR integrated with the model for the ARIS controller so that the entire closed-loop system is simulated. Preliminary performance predictions are made for the ARIS in attenuating both off-board disturbances as well as disturbances from hardware mounted onboard the microgravity facility. These predictions suggest that the ARIS does eliminate resonant behavior detrimental to microgravity experimentation. A limited comparison is made between the simulation predictions of ARIS attenuation of off-board disturbances and results from the ARIS flight test. These comparisons show promise, but further tuning of the simulation is needed.

Johnson, Terry L.↗

How Should Life Support Be Modeled and Simulated?

Why do most space life support research groups build and investigate large models for systems simulation? The need for them seems accepted, but are we asking the right questions and solving the real problems? The modeling results leave many questions unanswered. How then should space life support be modeled and simulated? Life support system research and development uses modeling and simulation to study dynamic behavior as part of systems engineering and analysis. It is used to size material flows and buffers and plan contingent operations. A DoD sponsored study used the systems engineering approach to define a set of best practices for modeling and simulation. These best practices describe a systems engineering process of developing and validating requirements, defining and analyzing the model concept, and designing and testing the model. Other general principles for modeling and simulation are presented. Some specific additional advice includes performing a static analysis before developing a dynamic simulation, applying the mass and energy conservation laws, modeling on the appropriate system level, using simplified subsystem representations, designing the model to solve a specific problem, and testing the model on several different problems. Modeling and simulation is necessary in life support design but many problems are outside its scope.

Jones, Harry W.↗

The NASA Space Radiation Laboratory

The Space Radiation Element (SRE) of the Human Research Program (HRP) aims to establish a robust biospecimen and tissue sharing collection (BTSC) to improve sample collection, tracking, access, distribution, and usability to maximize scientific return with intention. Future missions to the Moon, Mars, and beyond will require the appropriate characterization and mitigation of the many of the health risks associated with cosmic radiation exposure. An ongoing partnership between the National Aeronautics and Space Agency (NASA) and the Office of Science for the U.S. Department of Energy (DOE) has established the NASA Space Radiation Laboratory (NSRL) at the DOE Brookhaven National Laboratory (BNL) to create a space radiation ground-analog. NSRL provides the capability to use simulated space radiation with high energy ion beams for assessment of biological outcomes relevant space radiation risks as well as electronics testing for avionics safety. Current NSRL irradiation capabilities include a variety of exposure variables available for experimental needs: - Up to approximately 10 Gy/min dose-rates with a 20 × 20 cm 2 beam area. - Down to several hundred particles within a 60 × 60 cm 2 beam area. - Dose fractionation across multiple days is also possible. - For most radiobiology experiments the exposure structure, referred to as “spills”, has a 4 second repetition time where ions are extracted uniformly during a 0.3-0.4 second spill, followed by a ~3.6 second recovery time when the beam is off. - Both single or mixed-beam field irradiations are possible including a full 33-beam GCRSim and 6-beam SimGCRSim. Biological experiment capabilities: Two labs for Cell/tissue culture procedures can be conducted in two lab spaces with a separate lab available for sample processing that contains two chemical fume hoods. Sixteen fully equipped cell culture laboratories are available in the NSRL’s Long Term Support area in the Medical Building. The current vivarium configuration can support approximately 8,000 rodents/year for animal-based experiments. The NSRL area animal husbandry system includes five caging racks with 48 spaces in each coupled to a ventilation/filtration system in largest room, and two smaller rooms with conventional caging systems. Schedule: The NSRL is nominally open for experiment scheduling during three runs per year; Spring (March – May), Summer (May - June), and Fall (Sept – Nov). NASA funding is not required for biological and electronic testing experiments at the NSRL.

Jason M Weeks↗

Simulation of SEU Cross-sections using MRED under Conditions of Limited Device Information

This viewgraph presentation reviews the simulation of Single Event Upset (SEU) cross sections using the membrane electrode assembly (MEA) resistance and electrode diffusion (MRED) tool using "Best guess" assumptions about the process and geometry, and direct ionization, low-energy beam test results. This work will also simulate SEU cross-sections including angular and high energy responses and compare the simulated results with beam test data for the validation of the model. Using MRED, we produced a reasonably accurate upset response model of a low-critical charge SRAM without detailed information about the circuit, device geometry, or fabrication process

Lauenstein, J. M.↗

Prototype Demonstration of an Integrated Solar Concentrator System and Carbothermal Reactor Using Solar Energy to Extract Oxygen from Regolith

The Carbothermal Reduction Demonstration (CaRD) project was an effort to develop a prototype system to demonstrate the extraction of oxygen from simulated lunar regolith using concentrated solar energy and a carbothermal reaction. The prototype consisted of a deployable solar concentrator capable of tracking the sun, carbothermal reactor, fluid system, gas analysis, and a solar concentrator control system consisting of avionics and software. These subsystems were developed by multiple NASA centers and a private industry partner, Sierra Space. The various teams worked together to define requirements and interfaces to successfully assemble the complex system and demonstrate an integrated solar carbothermal process. The solar concentrator developed at Glenn Research Center (GRC) was designed to be stowed for a launch environment then deployed on the lunar surface. It utilized a crossed dragone configuration of composite mirrors to direct horizontal sunlight onto a target 90° from the incoming sunlight. The key performance parameters for the solar concentrator were efficiency and power density. The carbothermal reactor was developed by Sierra Space through a separate project called the Carbothermal Oxygen Production Reactor (COPR) where it successfully demonstrated a fully automated process in a thermal vacuum environment. The fluid system needed for the carbothermal reaction was also developed by Sierra Space and successfully demonstrated in the same thermal vacuum test. The gas analysis system was developed at Kennedy Space Center (KSC) and was required to determine the amount of oxygen extracted during each test. The gas analysis system was based on the Mass Spectrometer Observing Lunar Operations (MSOLO) instrument. Avionics and software for the CaRD prototype were also developed at KSC and based on experience with MSOLO avionics and software. The control system was designed to stow, deploy, track the sun, and perform beam alignment of the concentrated light. The prototype subsystems were integrated and tested at Johnson Space Center’s (JSC) Energy Systems Test Area. A heliostat was used to direct sunlight toward the prototype in a way that is representative of the sunlight conditions at the south pole of the Moon. When concentrated sunlight was focused on simulated lunar regolith within the reactor, the gas analysis team confirmed the presence of carbon monoxide gas, which confirmed that a solar carbothermal reaction took place. The key performance parameter for the integrated prototype was grams of oxygen extracted per kilowatt hour of energy arriving at the concentrator primary mirror. The prototype design successfully demonstrated end-to-end capability and further steps to achieve a flight capable system have been defined. With lunar data, engineers would be able to design a scaled-up system capable of extracting oxygen from regolith at useful quantities for crew life support and rocket propellant. On the long term, this method of In-Situ Resource Utilization could be used to drastically reduce the cost and risk of a sustained human presence on the Moon by reducing the amount of oxygen that would have to be delivered.

Oxygen from Regolith↗

The simulated features of heliospheric cosmic-ray modulation with a time-dependent drift model. II - On the energy dependence of the onset of new modulation in 1987

Time-dependent heliospheric cosmic-ray modulation for the period 1985-1989 is simulated by means of a time-dependent axially symmetric drift model with an emulated wavy heliospheric neutral sheet (HNS). The model is used to extend previous calculations to other energies in order to study a possible energy-dependence of the onset of new modulation at various radial distances in 1987. The model, with the outward propagating changes of the HNS as the only time-dependent parameters, is found to predict essentially no energy-dependence in the time when new modulation started in the simulated 1987. When a more practical approach in defining 'constant' modulation in 1987 is followed, the present calculations can be interpreted to indicate that the end of the recovery period in 1987 happened progressively earlier and the onset of new modulation progressively later the higher the rigidity of the cosmic rays. This period of relatively unchanged modulation is predicted to last longer with increasing radial distance.

Le Roux, J. A.↗

Voyager 2 plasma ion observations in the magnetosphere of Uranus

Positive ion measurements in the magnetosphere of Uranus have been made by the Voyager 2 plasma science experiment. The paper presents an overview of the entire data set and a detailed analysis of the observations from the inner magnetosphere which complements and extends results reported elsewhere. Densities and temperatures are obtained from an analysis which incorporates details of the instrumental response. These results are then used to calculate flux tube particle and energy content to support the hypothesis that the plasma transport is controlled by a solar wind-driven magnetospheric convection system. Variations in the flux tube content suggest both a local source of plasma, produced from the neutral hydrogen corona of Uranus, and a nonlocal source, convected inwared and heated by adiabatic compression. In each case a proton composition is inferred. Sharp boundaries in the high-energy (approximately 1 keV) plasma population are interpreted in terms of the spatial extent of the magnetospheric convection, with significant shielding of the convection electric field. The convection theory is also used in a simulation of the low-energy (approximately 10 eV) ion component using the neutral hydrogen source, resulting in distribution functions which qualitatively agree with the observations.

Selesnick, Richard S.↗

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↗

Radiation Protection Effectiveness of Polymeric Based Shielding Materials at Low Earth Orbit

Correlations of limited ionizing radiation measurements onboard the Space Transportation System (STS; shuttle) and the International Space Station (ISS) with numerical simulations of charged particle transport through spacecraft structure have indicated that usage of hydrogen rich polymeric materials improves the radiation shielding performance of space structures as compared to the traditionally used aluminum alloys. We discuss herein the radiation shielding correlations between measurements on board STS-81 (Atlantis, 1997) using four polyethylene (PE) spheres of varying radii, and STS-89 (Endeavour, 1998) using aluminum alloy spheres; with numerical simulations of charged particle transport using the Langley Research Center (LaRC)-developed High charge (Z) and Energy TRaNsport (HZETRN) algorithm. In the simulations, the Galactic Cosmic Ray (GCR) component of the ionizing radiation environment at Low Earth Orbit (LEO) covering ions in the 1< or equals Z< or equals 28 range is represented by O'Neill's (2004) model. To compute the transmission coefficient for GCR ions at LEO, O'Neill's model is coupled with the angular dependent LaRC cutoff model. The trapped protons/electrons component of LEO environment is represented by a LaRC-developed time dependent procedure which couples the AP8min/AP8max, Deep River Neutron Monitor (DRNM) and F10.7 solar radio frequency measurements. The albedo neutron environment resulting from interaction of GCR ions with upper atmosphere is modeled through extrapolation of the Atmospheric Ionizing Radiation (AIR) measurements. With the validity of numerical simulations through correlation with PE and aluminum spheres measurements established, we further present results from the expansion of the simulations through the selection of high hydrogen content commercially available polymeric constituents such as PE foam core and Spectra fiber(Registered TradeMark) composite face sheet to assess their radiation shield properties as compared to generic PE.

Badavi, Francis F.↗

Energy storage and dissipation in the magnetotail during substorms. I - Particle simulations. II - MHD simulations

2D electromagnetic particle simulations are used to investigate the dynamics of the tail during development of substorms under the influence of the pressure in the magnetospheric boundary layer and the dawn-to-dusk electric field. It is shown that pressure pulses result in thinning of the tail current sheet as the magnetic field becomes pinched near the region where the pressure pulse is applied. The pinching leads to the tailward flow of the current sheet plasma and the eventual formation and injection of a plasmoid. Surges in the dawn-to-dusk electric field cause plasma on the flanks to convect into the center of the current sheet, thereby thinning the current sheet. The pressure in the magnetospheric boundary laser is coupled to the dawn-to-dusk electric field through the conductivity of the tail. Changes in the predicted evolution of the magnetosphere during substorms due to changes in the resistivity are investigated under the assumption that MHD theory provides a suitable representation of the global or large-scale evolution of the magnetotail to changes in the solar wind and to reconnection at the dayside magnetopause. It is shown that the overall evolution of the magnetosphere is about the same for three different resistivity distributions with plasmoid formation and ejection in each case.

Winglee, R. M.↗

Micro-Energy Rates for Damage Tolerance and Durability of Composite Structures

In this paper, the adhesive bond strength of lap-jointed graphite/aluminum composites is examined by computational simulation. Computed micro-stress level energy release rates are used to identify the damage mechanisms associated with the corresponding acoustic emission (AE) signals. Computed damage regions are similarly correlated with ultrasonically scanned damage regions. Results show that computational simulation can be used with suitable NDE methods for credible in-service monitoring of composites.

Chamis, Christos C.↗

Forward-looking Assimilation of MODIS-derived Snow Covered Area into a Land Surface Model

Snow cover over land has a significant impact on the surface radiation budget, turbulent energy fluxes to the atmosphere, and local hydrological fluxes. For this reason, inaccuracies in the representation of snow covered area (SCA) within a land surface model (LSM) can lead to substantial errors in both offline and coupled simulations. Data assimilation algorithms have the potential to address this problem. However, the assimilation of SCA observations is complicated by an information deficit in the observation SCA indicates only the presence or absence of snow, and not snow volume and by the fact that assimilated SCA observations can introduce inconsistencies with atmospheric forcing data, leading to non-physical artifacts in the local water balance. In this paper we present a novel assimilation algorithm that introduces MODIS SCA observations to the Noah LSM in global, uncoupled simulations. The algorithm utilizes observations from up to 72 hours ahead of the model simulation in order to correct against emerging errors in the simulation of snow cover while preserving the local hydrologic balance. This is accomplished by using future snow observations to adjust air temperature and, when necessary, precipitation within the LSM. In global, offline integrations, this new assimilation algorithm provided improved simulation of SCA and snow water equivalent relative to open loop integrations and integrations that used an earlier SCA assimilation algorithm. These improvements, in turn, influenced the simulation of surface water and energy fluxes both during the snow season and, in some regions, on into the following spring.

Zaitchik, Benjamin F.↗

A Computational Fluid Dynamic and Heat Transfer Model for Gaseous Core and Gas Cooled Space Power and Propulsion Reactors

A computational model based on the axisymmetric, thin-layer Navier-Stokes equations is developed to predict the convective, radiation and conductive heat transfer in high temperature space nuclear reactors. An implicit-explicit, finite volume, MacCormack method in conjunction with the Gauss-Seidel line iteration procedure is utilized to solve the thermal and fluid governing equations. Simulation of coolant and propellant flows in these reactors involves the subsonic and supersonic flows of hydrogen, helium and uranium tetrafluoride under variable boundary conditions. An enthalpy-rebalancing scheme is developed and implemented to enhance and accelerate the rate of convergence when a wall heat flux boundary condition is used. The model also incorporated the Baldwin and Lomax two-layer algebraic turbulence scheme for the calculation of the turbulent kinetic energy and eddy diffusivity of energy. The Rosseland diffusion approximation is used to simulate the radiative energy transfer in the optically thick environment of gas core reactors. The computational model is benchmarked with experimental data on flow separation angle and drag force acting on a suspended sphere in a cylindrical tube. The heat transfer is validated by comparing the computed results with the standard heat transfer correlations predictions. The model is used to simulate flow and heat transfer under a variety of design conditions. The effect of internal heat generation on the heat transfer in the gas core reactors is examined for a variety of power densities, 100 W/cc, 500 W/cc and 1000 W/cc. The maximum temperature, corresponding with the heat generation rates, are 2150 K, 2750 K and 3550 K, respectively. This analysis shows that the maximum temperature is strongly dependent on the value of heat generation rate. It also indicates that a heat generation rate higher than 1000 W/cc is necessary to maintain the gas temperature at about 3500 K, which is typical design temperature required to achieve high efficiency in the gas core reactors. The model is also used to predict the convective and radiation heat fluxes for the gas core reactors. The maximum value of heat flux occurs at the exit of the reactor core. Radiation heat flux increases with higher wall temperature. This behavior is due to the fact that the radiative heat flux is strongly dependent on wall temperature. This study also found that at temperature close to 3500 K the radiative heat flux is comparable with the convective heat flux in a uranium fluoride failed gas core reactor.

Anghaie, S.↗

Magnetohydrodynamic MACH Code Used to Simulate Magnetoplasmadynamic Thrusters

The On-Board Propulsion program at the NASA Glenn Research Center is utilizing a state of-the-art numerical simulation to model the performance of high-power electromagnetic plasma thrusters. Such thrusters are envisioned for use in lunar and Mars cargo transport, piloted interplanetary expeditions, and deep-space robotic exploration of the solar system. The experimental portion of this program is described in reference 1. This article describes the numerical modeling program used to guide the experimental research. The synergistic use of numerical simulations and experimental research has spurred the rapid advancement of high-power thruster technologies for a variety of bold new NASA missions. From its inception as a U.S. Department of Defense code in the mid-1980's, the Multiblock Arbitrary Coordinate Hydromagnetic (MACH) simulation tool has been used by the plasma physics community to model a diverse range of plasma problems--including plasma opening switches, inertial confinement fusion concepts, compact toroid formation and acceleration, z-pinch implosion physics, laser-target interactions, and a variety of plasma thrusters. The MACH2 code used at Glenn is a time-dependent, two-dimensional, axisymmetric, multimaterial code with a multiblock structure. MACH3, a more recent three-dimensional version of the code, is currently undergoing beta tests. The MACH computational mesh moves in an arbitrary Lagrangian-Eulerian (ALE) fashion that allows the simulation of diffusive-dominated and dispersive-dominated problems, and the mesh can be refined via a variety of adaptive schemes to capture regions of varying characteristic scale. The mass continuity and momentum equations model a compressible viscous fluid, and three energy equations are used to simulate nonthermal equilibrium between electrons, ions, and the radiation field. Magnetic fields are modeled by an induction equation that includes resistive diffusion, the Hall effect, and a thermal source for magnetic fields. Various models of plasma resistivity are included, along with ablation models and multiport circuit solvers. The set of equations is closed using either an ideal gas or real equation of state.

Mikellides, Pavlos G.↗

A Technique for Assimilating GOES-Derived Land Surface Products into Regional Models to Improve the Representation of Land Surface Forcing

As the parameterizations of surface energy budgets in regional models have become more complete physically, models have the potential to be much more realistic in simulations of coupling between surface radiation, hydrology, and surface energy transfer. Realizing the importance of properly specifying the surface energy budget, many institutions are using land-surface models to represent the lower boundary forcing associated with biophysical processes and soil hydrology. However, the added degrees of freedom due to inclusion of such land-surface schemes require the specification of additional parameters within the model system such as vegetative resistances, green vegetation fraction, leaf area index, soil physical and hydraulic characteristics, stream flow, runoff, and the vertical distribution of soil moisture. A technique has been developed for assimilating GOES-IR skin temperature tendencies into the surface energy budget equation of a mesoscale model so that the simulated rate of temperature change closely agrees with the satellite observations. A critical assumption of the technique is that the availability of moisture (either from the soil or vegetation) is the least known term in the model's surface energy budget. Therefore, the simulated latent heat flux, which is a function of surface moisture availability, is adjusted based upon differences between the modeled and satellite-observed skin temperature tendencies. An advantage of this technique is that satellite temperature tendencies are assimilated in an energetically consistent manner that avoids energy imbalances and surface stability problems that arise from direct assimilation of surface shelter temperatures. The fact that the rate of change of the satellite skin temperature is used rather than the absolute temperature means that sensor calibration is not as critical. An advantage of this technique for short-range forecasts is that it does not require a complex land-surface formulation within the atmospheric model. As a result, the need to specify poorly known soil and vegetative characteristics is eliminated. The GOES assimilation technique has been incorporated into the PSU/NCAR MM5. Results will be presented to demonstrate the ability of the assimilation scheme to improve short-term simulations of near-surface air temperature and mixing ratio during the warm season for several selected cases during June of 1995 and June 1997 which exhibit a variety of atmospheric and land-surface conditions. In addition, the simulations produced by the assimilation technique will be compared with those produced by the MM5/BATS coupled model initialized with soil water content derived from a soil hydrology model.

Lapenta, William M.↗

Calculating Free Energies Using Scaled-Force Molecular Dynamics Algorithm

One common objective of molecular simulations in chemistry and biology is to calculate the free energy difference between different states of the system of interest. Examples of problems that have such an objective are calculations of receptor-ligand or protein-drug interactions, associations of molecules in response to hydrophobic, and electrostatic interactions or partition of molecules between immiscible liquids. Another common objective is to describe evolution of the system towards a low energy (possibly the global minimum energy), 'native' state. Perhaps the best example of such a problem is folding of proteins or short RNA molecules. Both types of problems share the same difficulty. Often, different states of the system are separated by high energy barriers, which implies that transitions between these states are rare events. This, in turn, can greatly impede exploration of phase space. In some instances this can lead to 'quasi non-ergodicity', whereby a part of phase space is inaccessible on timescales of the simulation. A host of strategies has been developed to improve efficiency of sampling the phase space. For example, some Monte Carlo techniques involve large steps which move the system between low-energy regions in phase space without the need for sampling the configurations corresponding to energy barriers (J-walking). Most strategies, however, rely on modifying probabilities of sampling low and high-energy regions in phase space such that transitions between states of interest are encouraged. Perhaps the simplest implementation of this strategy is to increase the temperature of the system. This approach was successfully used to identify denaturation pathways in several proteins, but it is clearly not applicable to protein folding. It is also not a successful method for determining free energy differences. Finally, the approach is likely to fail for systems with co-existing phases, such as water-membrane systems, because it may lead to spontaneous mixing. A similar difficulty may be encountered in any method relying on global modifications of phase space.

Darve, Eric↗