Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational molecular design”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Concept for Inclusion of Analytical and Computational Capability in Optical Plume Anomaly Detection (OPAD) for Measurement of Neutron Flux

Researchers are working on many fronts to make possible high-speed, automated classification and quantification of constituent materials in numerous environments. NASA's Marshall Space Flight Center has implemented a system for rocket engine flowfields/plumes. The Optical Plume Anomaly Detector (OPAD) system was designed to utilize emission and absorption spectroscopy for monitoring molecular and atomic particulates in gas plasma. An accompanying suite of tools and analytical package designed to utilize information collected by OPAD is known as the Engine Diagnostic Filtering System (EDiFiS). The current combination of these systems identifies atomic and molecular species and quantifies mass loss rates in H2/O2 rocket plumes. Capabilities for real-time processing are being advanced on several fronts, including an effort to hardware encode components of the EDiFiS for health monitoring and management. This paper addresses the OPAD with its tool suites, and discusses what is considered a natural progression: a concept for taking OPAD to the next logical level of high energy physics, incorporating fermion and boson particle analyses in measurement of neutron flux.

Patrick, Marshall Clint↗
NASA NTRS2004

Concept for Inclusion of Analytical and Computational Capability in Optical Plume Anomaly Detection (OPAD) for Measurement of Neutron Flux

Researchers are working on many konts to make possible high speed, automated classification and quantification of constituent materials in numerous environments. NASA's Marshall Space Flight Center has implemented a system for rocket engine flow fields/plumes; the Optical Plume Anomaly Detection (OPAD) system was designed to utilize emission and absorption spectroscopy for monitoring molecular and atomic particulates in gas plasma. An accompanying suite of tools and analytical package designed to utilize information collected by OPAD is known as the Engine Diagnostic Filtering System (EDIFIS). The current combination of these systems identifies atomic and molecular species and quantifies mass loss rates in H2/O2 rocket plumes. Additionally, efforts are being advanced to hardware encode components of the EDIFIS in order to address real-time operational requirements for health monitoring and management. This paper addresses the OPAD with its tool suite, and discusses what is considered a natural progression: a concept for migrating OPAD towards detection of high energy particles, including neutrons and gamma rays. The integration of these tools and capabilities will provide NASA with a systematic approach to monitor space vehicle internal and external environment.

Patrick, M. Clinton↗
NASA NTRS2019

Aerodynamics for the ADEPT SR-1 Flight Experiment

Adaptable, Deployable, Entry, and Placement Technology (ADEPT) is a combination of a heatshield and an aerodynamic decelerator for atmospheric entry applications. The ADEPT Sounding Rocket (SR)-1 mission was a suborbital flight experiment of an 0.7 m-diameter ADEPT to verify system-level performance and to characterize dynamic stability behavior. The aerodynamic database for ADEPT SR-1 was constructed from non-continuum and continuum flowfield computations, along with data from recent ADEPT ground testing and the IRVE-3 flight test vehicle. High-altitude (free-molecular and transitional regimes) data were generated using DSMC methods. Pre-flight predictions of continuum static aerodynamics coefficients were derived from Reynolds-Averaged Navier-Stokes solutions at conditions along a design trajectory, with comparisons to available ground test data of the nano-ADEPT geometry. Dynamic pitch damping characteristics were taken from functional forms developed for the IRVE-3 flight test vehicle through ballistic range testing. Comparison of pre-flight predictions to post-flight reconstruction of aerodynamic force and moment coefficients is presented.

Korzun, Ashley M.↗
NASA NTRS2001

Experiments and Model Development for the Investigation of Sooting and Radiation Effects in Microgravity Droplet Combustion

Today, despite efforts to develop and utilize natural gas and renewable energy sources, nearly 97% of the energy used for transportation is derived from combustion of liquid fuels, principally derived from petroleum. While society continues to rely on liquid petroleum-based fuels as a major energy source in spite of their finite supply, it is of paramount importance to maximize the efficiency and minimize the environmental impact of the devices that burn these fuels. The development of improved energy conversion systems, having higher efficiencies and lower emissions, is central to meeting both local and regional air quality standards. This development requires improvements in computational design tools for applied energy conversion systems, which in turn requires more robust sub-model components for combustion chemistry, transport, energy transport (including radiation), and pollutant emissions (soot formation and burnout). The study of isolated droplet burning as a unidimensional, time dependent model diffusion flame system facilitates extensions of these mechanisms to include fuel molecular sizes and pollutants typical of conventional and alternative liquid fuels used in the transportation sector. Because of the simplified geometry, sub-model components from the most detailed to those reduced to sizes compatible for use in multi-dimensional, time dependent applied models can be developed, compared and validated against experimental diffusion flame processes, and tested against one another. Based on observations in microgravity experiments on droplet combustion, it appears that the formation and lingering presence of soot within the fuel-rich region of isolated droplets can modify the burning rate, flame structure and extinction, soot aerosol properties, and the effective thermophysical properties. These observations led to the belief that perhaps one of the most important outstanding contributions of microgravity droplet combustion is the observation that in the absence of asymmetrical forced and natural convection, a soot shell is formed between the droplet surface and the flame, exerting an influence on the droplet combustion response far greater than previously recognized. The effects of soot on droplet burning parameters, including burning rate, soot shell dynamics, flame structure, and extinction phenomena provide significant testing parameters for studying the structure and coupling of soot models with other sub-model components.

Choi, Mun Young↗
NASA NTRS2013

Raster-Based Approach to Solar Pressure Modeling

An algorithm has been developed to take advantage of the graphics processing hardware in modern computers to efficiently compute high-fidelity solar pressure forces and torques on spacecraft, taking into account the possibility of self-shading due to the articulation of spacecraft components such as solar arrays. The process is easily extended to compute other results that depend on three-dimensional attitude analysis, such as solar array power generation or free molecular flow drag. The impact of photons upon a spacecraft introduces small forces and moments. The magnitude and direction of the forces depend on the material properties of the spacecraft components being illuminated. The parts of the components being lit depends on the orientation of the craft with respect to the Sun, as well as the gimbal angles for any significant moving external parts (solar arrays, typically). Some components may shield others from the Sun. The purpose of this innovation is to enable high-fidelity computation of solar pressure and power generation effects of illuminated portions of spacecraft, taking self-shading from spacecraft attitude and movable components into account. The key idea in this innovation is to compute results dependent upon complicated geometry by using an image to break the problem into thousands or millions of sub-problems with simple geometry, and then the results from the simpler problems are combined to give high-fidelity results for the full geometry. This process is performed by constructing a 3D model of a spacecraft using an appropriate computer language (OpenGL), and running that model on a modern computer's 3D accelerated video processor. This quickly and accurately generates a view of the model (as shown on a computer screen) that takes rotation and articulation of spacecraft components into account. When this view is interpreted as the spacecraft as seen by the Sun, then only the portions of the craft visible in the view are illuminated. The view as shown on the computer screen is composed of up to millions of pixels. Each of those pixels is associated with a small illuminated area of the spacecraft. For each pixel, it is possible to compute its position, angle (surface normal) from the view direction, and the spacecraft material (and therefore, optical coefficients) associated with that area. With this information, the area associated with each pixel can be modeled as a simple flat plate for calculating solar pressure. The vector sum of these individual flat plate models is a high-fidelity approximation of the solar pressure forces and torques on the whole vehicle. In addition to using optical coefficients associated with each spacecraft material to calculate solar pressure, a power generation coefficient is added for computing solar array power generation from the sum of the illuminated areas. Similarly, other area-based calculations, such as free molecular flow drag, are also enabled. Because the model rendering is separated from other calculations, it is relatively easy to add a new model to explore a new vehicle or mission configuration. Adding a new model is performed by adding OpenGL code, but a future version might read a mesh file exported from a computer-aided design (CAD) system to enable very rapid turnaround for new designs

Wright, Theodore W. II↗
NASA NTRSDate not supplied

“Cold-Flow” Experiments Supporting CFD of Mixing Flowfields for High-Speed Fuel Injectors for Scramjet Applications

Recent flight demonstrations of supersonic combustion ramjet (scramjet) vehicles prove their increasing promise for military (rapid response and strike capability on global scale), aerospace (safer and more affordable access to space), and civil aviation (hypersonic point-to-point transport) applications. Currently, these technologies are still in their early development stages with commercial interest and investment at only a fraction of that of government organizations such as NASA and DoD. To advance hypersonic air-breathing propulsion technologies to the technology readiness levels necessary for access to space or widespread commercialization, further government investments and university engagement are needed over the next decade and likely beyond. This is because designing scramjet propulsion devices capable of robust high-speed air-breathing operation, characterized by rapid fuel and air mixing, short combustion times, and ensuring stable flame, has proven difficult. In this presentation, after a brief engineering-level introduction to scramjets technology, we will discuss one of the key challenges in high-speed propulsion design, namely the fuel injector design. Attempts at improving fuel injection to enhance fuel-air mixing while simultaneously reducing engine thrust losses have received a great deal of attention over the years. Although some amount of loss is thermodynamically unavoidable and occurs due to the desired effect of molecular mixing of fuel and air, any losses beyond this minimum required amount reduces the thrust potential of the engine. The Enhanced Injection and Mixing Project (EIMP) at NASA Langley Research Center aims at addressing this design challenge by analyzing the performance of a number of baseline and novel fuel injectors for high-speed applications. The project leverages computational and experimental capabilities with the goal to investigate scramjet fuel injection and mixing physics, improve our understanding of the underlying physical process, and develop enhancement strategies relevant to hypersonic flight Mach numbers. The talk will discuss the current computational and experimental research approaches using one of the baseline injectors considered by the EIMP as an example.

Scramjet↗
NASA NTRS1982

Electronic levels and charge distribution near the interface of nickel

The energy levels in clusters of nickel atoms were investigated by means of a series of cluster calculations using both the multiple scattering and computational techniques (designated SSO) which avoids the muffin-tin approximation. The point group symmetry of the cluster has significant effect on the energy of levels nominally not occupied. This influences the electron transfer process during chemisorption. The SSO technique permits the approaching atom or molecule plus a small number of nickel atoms to be treated as a cluster. Specifically, molecular levels become more negative in the O atom, as well as in a CO molecule, as the metal atoms are approached. Thus, electron transfer from the nickel and bond formation is facilitated. This result is of importance in understanding chemisorption and catalytic processes.

Waber, J. T.↗
NASA NTRS1991

The solution of radiative transfer problems in molecular bands without the LTE assumption by accelerated lambda iteration methods

An iterative method based on the use of approximate transfer operators, which was designed initially to solve multilevel NLTE line formation problems in stellar atmospheres, is adapted and applied to the solution of the NLTE molecular band radiative transfer in planetary atmospheres. The matrices to be constructed and inverted are much smaller than those used in the traditional Curtis matrix technique, which makes possible the treatment of more realistic problems using relatively small computers. This technique converges much more rapidly than straightforward iteration between the transfer equation and the equations of statistical equilibrium. A test application of this new technique to the solution of NLTE radiative transfer problems for optically thick and thin bands (the 4.3 micron CO2 band in the Venusian atmosphere and the 4.7 and 2.3 micron CO bands in the earth's atmosphere) is described.

Kutepov, A. A.↗
NASA NTRS2005

Atom Resonance Lines for Modeling Atmosphere: Studies of Pressure-Broadening of Alkali Atom Resonance Lines for Modeling Atmospheres of Extrasolar Giant Planets and Brown Dwarfs

We report on progress made in a joint program of theoretical and experimental research to study the line-broadening of alkali atom resonance lines due to collisions with species such as helium and molecular hydrogen. Accurate knowledge of the line profiles of Na and K as a function of temperature and pressure will allow such lines to serve as valuable diagnostics of the atmospheres of brown dwarfs and extra-solar giant planets. A new experimental apparatus has been designed, built and tested over the past year, and we are poised to begin collecting data on the first system of interest, the potassium resonance lines perturbed by collisions with helium. On the theoretical front, calculations of line-broadening due to sodium collisions with helium are nearly complete, using accurate molecular potential energy curves and transition moments just recently computed for this system. In addition we have completed calculations of the three relevant potential energy curves and associated transition moments for K - He, using the MOLPRO quantum chemistry codes. Currently, calculations of the potential surfaces describing K-H2 are in progress.

Hasan, Hashima↗
NASA NTRS1969

Models of Earth's Atmosphere (120 to 1000 Km)

Atmospheric conditions encountered by a spacecraft in orbit about the Earth are important factors in space vehicle design, mission planning, and mission operations. Density is the primary atmospheric property that affects the spacecraft's orbital altitude, lifetime, and motion in the altitude range of 120 to 1000 kilometers. Near the lower limit of this range where density is greatest, a spacecraft will generally remain in orbit for a very short time; near the upper limit, the density effect on orbital lifetime is almost negligible. Density directly affects the torques which result from aerodynamic interaction between the space vehicle and the atmosphere; such torques must be considered in design of spacecraft attitude control systems. Density scale height is required in heating calculations for space vehicles re-entering the Earth's upper atmosphere. Density as well as chemical composition and temperature are needed in calculating a spacecraft's drag coefficient. Chemical composition and temperature also are required in the design of experiment sensors to be flown in this altitude range. Because of variability of atmospheric conditions with spatial location and solar condition, invariant models of the Earth's atmosphere (120 to 1000 kilometers) would not be useful for most engineering applications. Therefore, this monograph presents a computerized version of Jacchia's prediction method to provide models of the Earth's atmosphere which vary with solar condition and location. The resulting atmospheric models, which are predicted for particular times and locations, provide atmospheric density, chemical composition, temperature, molecular mass, and density scale height between 120 and 1000 kilometers altitude. In addition to the computerized method, a quick-look prediction method is given that may be used to obtain an estimate of atmospheric density for any time and spatial location without the use of a computer. A sample problem illustrates this method. Both methods provide models of mean density and models having reasonable upper extremes for density. The analytical approaches in both methods are considered to be the best available, but they could be refined considerably by additional data and study. Information contained in this monograph applies to altitudes between 120 and 1000 kilometers; other design criteria monographs in this series will provide atmospheric information below this altitude region.

Weidner, D. K.↗
NASA NTRS2014

Modeling Materials: Design for Planetary Entry, Electric Aircraft, and Beyond

NASA missions push the limits of what is possible. The development of high-performance materials must keep pace with the agency's demanding, cutting-edge applications. Researchers at NASA's Ames Research Center are performing multiscale computational modeling to accelerate development times and further the design of next-generation aerospace materials. Multiscale modeling combines several computationally intensive techniques ranging from the atomic level to the macroscale, passing output from one level as input to the next level. These methods are applicable to a wide variety of materials systems. For example: (a) Ultra-high-temperature ceramics for hypersonic aircraft-we utilized the full range of multiscale modeling to characterize thermal protection materials for faster, safer air- and spacecraft, (b) Planetary entry heat shields for space vehicles-we computed thermal and mechanical properties of ablative composites by combining several methods, from atomistic simulations to macroscale computations, (c) Advanced batteries for electric aircraft-we performed large-scale molecular dynamics simulations of advanced electrolytes for ultra-high-energy capacity batteries to enable long-distance electric aircraft service; and (d) Shape-memory alloys for high-efficiency aircraft-we used high-fidelity electronic structure calculations to determine phase diagrams in shape-memory transformations. Advances in high-performance computing have been critical to the development of multiscale materials modeling. We used nearly one million processor hours on NASA's Pleiades supercomputer to characterize electrolytes with a fidelity that would be otherwise impossible. For this and other projects, Pleiades enables us to push the physics and accuracy of our calculations to new levels.

Supercomputing↗
NASA NTRS2002

Equivalent-Continuum Modeling With Application to Carbon Nanotubes

A method has been proposed for developing structure-property relationships of nano-structured materials. This method serves as a link between computational chemistry and solid mechanics by substituting discrete molecular structures with equivalent-continuum models. It has been shown that this substitution may be accomplished by equating the vibrational potential energy of a nano-structured material with the strain energy of representative truss and continuum models. As important examples with direct application to the development and characterization of single-walled carbon nanotubes and the design of nanotube-based devices, the modeling technique has been applied to determine the effective-continuum geometry and bending rigidity of a graphene sheet. A representative volume element of the chemical structure of graphene has been substituted with equivalent-truss and equivalent continuum models. As a result, an effective thickness of the continuum model has been determined. This effective thickness has been shown to be significantly larger than the interatomic spacing of graphite. The effective thickness has been shown to be significantly larger than the inter-planar spacing of graphite. The effective bending rigidity of the equivalent-continuum model of a graphene sheet was determined by equating the vibrational potential energy of the molecular model of a graphene sheet subjected to cylindrical bending with the strain energy of an equivalent continuum plate subjected to cylindrical bending.

Odegard, Gregory M.↗
NASA NTRS1971

Dynamic upper atmospheric force model on stabilized vehicles for a high-precision trajectory computer program

The upper atmosphere model draws heavily on the behavior of the earth's upper atmosphere which exhibits cyclic as well as irregular variations in density profile, temperature, pressure, and composition in unison with solar activities as deduced from the more recent land-based and satellite observations. The lift and drag model is designed specifically for inertially stabilized vehicles of the Mariner class, with possible extension to gravity gradient stabilized vehicles of the GEOS class. The model considers operation in the free molecular flow regimes with large Knudsen numbers. The vehicle is considered a composite structure with basic components having well-defined shapes, each with its own surface characteristics in terms of temperature, reflectivity, and accommodation of free stream molecules. The model takes into account both the calculation of precise aerodynamic force coefficients in terms of expansion of modified Bessel functions in speed ratios and angle of attack, and approximate force coefficients when the speed ratios approach infinity. Other considerations include specular and diffused reflectivity, shielding, and shadow effects.

Khatib, A. R.↗
NASA NTRSDate not supplied

Advanced Materials for the Lunar Surface: Multiscale Computational Design of Refractory Alloys and Carbides

Emerging operational environments, such as the lunar surface, present novel challenges for NASA and drive the need for advanced materials in applications like fission surface power systems. To address these demands, computational materials science is rapidly evolving to augment or replace costly and hazardous empirical testing. Although materials selection at NASA remains predominantly experimentally driven, advanced simulation methodologies are being steadily integrated into the engineering lifecycle. This work details the application of multiscale simulation techniques—including first-principles calculations, CALPHAD, dislocation dynamics, and molecular dynamics—at NASA's Ames Research Center to evaluate advanced materials for extreme environments. First, we present contributions to the Space Nuclear Propulsion Project. Be-cause propellant channel coatings in nuclear thermal rockets must withstand high-pressure, high-temperature hydro-gen, optimizing these materials is critical. First-principles calculations were employed to establish a rigorous quantitative and qualitative understanding of the behavior of the refractory carbides ZrC, NbC, and their mixtures in high-enthalpy hydrogen environments. This necessitated the generation of high-fidelity thermodynamic models for both stoichiometric and carbon-depleted carbides, both with and without the presence of hydrogen. Furthermore, we highlight efforts under the Refractory Alloy Additive Manufacturing Build Optimization (RAAMBO) project, where existing and novel alloy compositions were assessed for additive manufacturing printability and subsequent performance in applications such as heat pipes and rocket nozzle extensions. This was accomplished through a comprehensive multiscale simulation framework that bridged the gap from the nanometer to the millimeter scale. Across both initiatives, rigorous validation against empirical data was prioritized. By systematically employing a verified and validated computational frame-work, we demonstrate how simulation effectively supports multidisciplinary engineering efforts, builds project-wide confidence, and drives critical materials development.

computational materials↗
NASA NTRS2008

Redesigned Human Metabolic Simulator

A design has been formulated for a proposed improved version of an apparatus that simulates atmospheric effects of human respiration by introducing controlled amounts of carbon dioxide, water vapor, and heat into the air. Denoted a human metabolic simulator (HMS), the apparatus is used for testing life-support equipment when human test subjects are not available. The prior version of the HMS, to be replaced, was designed to simulate the respiratory effects of as many as four persons. It exploits the catalytic combustion of methyl acetate, for which the respiratory quotient (the molar ratio of carbon dioxide produced to oxygen consumed) is very close to the human respiratory quotient of about 0.86. The design of the improved HMS provides for simulation of the respiratory effects of as many as eight persons at various levels of activity. The design would also increase safety by eliminating the use of combustion. The improved HMS (see figure) would include a computer that would exert overall control. The computer would calculate the required amounts of oxygen removal, carbon dioxide addition, water addition, and heat addition by use of empirical equations for metabolic profiles of respiration and heat. A blower would circulate air between the HMS and a chamber containing a life-support system to be tested. With the help of feedback from a mass flowmeter, the blower speed would be adjusted to regulate the rate of flow according to the number of persons to be simulated and to a temperature-regulation requirement (the air temperature would indirectly depend on the rate of flow, among other parameters). Oxygen would be removed from the circulating air by means of a commercially available molecular sieve configured as an oxygen concentrator. Oxygen, argon, and trace amounts of nitrogen would pass through a bed in the molecular sieve while carbon dioxide, the majority of nitrogen, and other trace gases would be trapped by the bed and subsequently returned to the chamber. If, as recommended, the oxygen concentrator were of a rotating twelve-bed design, then variations in the product stream could be made very small. Carbon dioxide would be added directly to the circulating air by simple injection from a supply tank. The rate of injection would be maintained at the required rate by use of a mass flowmeter/controller. In the same way, nitrogen would be added to make up for the small amount of nitrogen lost through the oxygen concentrator. Water vapor would be added to the circulating air by heating the corresponding required flow of water to steam in a heat exchanger. More heat, required to complete the simulation of the thermal effect of respiration, would be added through another heat exchanger. Heat would be supplied to both heat exchangers via a hot-oil loop.

Duffield, Bruce↗
NASA NTRS1991

Thermal modeling of nickel-hydrogen battery cells operating under transient orbital conditions

An analytical study of the thermal operating characteristics of nickel-hydrogen battery cells is presented. Combined finite-element and finite-difference techniques are employed to arrive at a computationally efficient composite thermal model representing a series-cell arrangement operating in conjunction with a radiately coupled baseplate and coldplate thermal bus. An aggressive, low-mass design approach indicates that thermal considerations can and should direct the design of the thermal bus arrangement. Special consideration is given to the potential for mixed conductive and convective processes across the hydrogen gap. Results of a compressible flow model are presented and indicate the transfer process is suitably represented by molecular conduction. A high-fidelity thermal model of the cell stack (and related components) indicates the presence of axial and radial temperature gradients. A detailed model of the thermal bus reveals the thermal interaction of individual cells and is imperative for assessing the intercell temperature gradients.

Schrage, Dean S.↗
NASA NTRS2017

Standard Isotherm Fit Information for Dry CO2 on Sorbents for 4-Bed Molecular Sieve

Onboard the ISS, one of the systems tasked with removal of metabolic carbon dioxide (CO2) is a 4-bed molecular sieve (4BMS) system. In order to enable a 4-person mission to succeed, systems for removal of metabolic CO2 must reliably operate for several years while minimizing power, mass, and volume requirements. This minimization can be achieved through system redesign and/or changes to the separation material(s). A material screening process has identified the most reliable sorbent materials for the next 4BMS. Sorbent characterization will provide the information necessary to guide system design by providing inputs for computer simulations.

Cmarik, G. E.↗
NASA NTRS2015

Free-Energy Calculations. A Mathematical Perspective

Ion channels are pore-forming assemblies of transmembrane proteins that mediate and regulate ion transport through cell walls. They are ubiquitous to all life forms. In humans and other higher organisms they play the central role in conducting nerve impulses. They are also essential to cardiac processes, muscle contraction and epithelial transport. Ion channels from lower organisms can act as toxins or antimicrobial agents, and in a number of cases are involved in infectious diseases. Because of their important and diverse biological functions they are frequent targets of drug action. Also, simple natural or synthetic channels find numerous applications in biotechnology. For these reasons, studies of ion channels are at the forefront of biophysics, structural biology and cellular biology. In the last decade, the increased availability of X-ray structures has greatly advanced our understanding of ion channels. However, their mechanism of action remains elusive. This is because, in order to assist controlled ion transport, ion channels are dynamic by nature, but X-ray crystallography captures the channel in a single, sometimes non-native state. To explain how ion channels work, X-ray structures have to be supplemented with dynamic information. In principle, molecular dynamics (MD) simulations can aid in providing this information, as this is precisely what MD has been designed to do. However, MD simulations suffer from their own problems, such as inability to access sufficiently long time scales or limited accuracy of force fields. To assess the reliability of MD simulations it is only natural to turn to the main function of channels - conducting ions - and compare calculated ionic conductance with electrophysiological data, mainly single channel recordings, obtained under similar conditions. If this comparison is satisfactory it would greatly increase our confidence that both the structures and our computational methodologies are sufficiently accurate. Channel conductance, defined as the ratio of ionic current through the channel to applied voltage, can be calculated in MD simulations by way of applying an external electric field to the system and counting the number of ions that traverse the channel per unit time. If the current is small, a voltage significantly higher than the experimental one needs to be applied to collect sufficient statistics of ion crossing events. Then, the calculated conductance has to be extrapolated to the experimental voltage using procedures of unknown accuracy. Instead, we propose an alternative approach that applies if ion transport through channels can be described with sufficient accuracy by the one-dimensional diffusion equation in the potential given by the free energy profile and applied voltage. Then, it is possible to test the assumptions of the equation, recover the full voltage/current dependence, determine the reliability of the calculated conductance and reconstruct the underlying (equilibrium) free energy profile, all from MD simulations at a single voltage. We will present the underlying theory, model calculations that test this theory and simulations on ion conductance through a channel that has been extensively studied experimentally. To our knowledge this is the first case in which the complete, experimentally measured dependence of the current on applied voltage has been reconstructed from MD simulations.

free energy↗