Search NASA⌕ Search

SEARCH · Search NASA

Results for “phase field 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 361 records · Page 20

Precipitation Processes Derived from TRMM Satellite Data, Cloud Resolving Model and Field Campaigns

Rainfall is a key link in the hydrologic cycle and is a primary heat source for the atmosphere. The vertical distribution of latent-heat release, which is accompanied by rainfall, modulates the large-scale circulations of the tropics and in turn can impact midlatitude weather. This latent heat release is a consequence of phase changes between vapor, liquid. and solid water. Present large-scale weather and climate models can simulate cloud latent heat release only crudely thus reducing their confidence in predictions on both global and regional scales. In this paper, NASA Tropical Rainfall Measuring (TRMM) precipitation radar (PR) derived rainfall information and the Goddard Convective and Stratiform Heating (CSH) algorithm used to estimate the four-dimensional structure of global monthly latent heating and rainfall profiles over the global tropics from December 1997 to October 2000. Rainfall latent heating and radar reflectively structure between ENSO (1997-1998 winter) and non-ENSO (1998-1999 winter) periods are examined and compared. The seasonal variation of heating over various geographic locations (i.e. Indian ocean vs west Pacific; Africa vs S. America) are also analyzed. In addition, the relationship between rainfall latent heating maximum heating level), radar reflectively and SST are examined.

Tao, W.-K.↗

Application of Advanced Very High Resolution Radiometer vegetation index to study atmosphere-biosphere exchange of CO2

Normalized difference vegetation indices derived from radiances measured by the Advanced Very High Resolution Radiometer were used to prescribe the phasing of terrestrial photosynthesis. The satellite data were combined with field data on soil respiration and a global map of net primary productivity to obtain the seasonal exchange of CO2 between the atmosphere and the terrestrial biosphere. The monthly fluxes of CO2 thus obtained were employed as source/sink functions in a global three-dimensional atmospheric tracer transport model to simulate the annual oscillations of CO2 in the atmosphere. The results demonstrate that satellite data of high spatial and temporal resolution can be used to provide quantitative information about seasonal and longer-term variations of photosynthetic activity on a global scale.

Fung, I. Y.↗

Particle tracing in the magnetosphere - New algorithms and results

A fast and efficient method is employed to trace charged particles through realistic magnetospheric electric and magnetic fields, greatly reducing computer simulation times. The method is applicable for particles having arbitrary charge, energy, or pitch angle but which conserve the first two adiabatic invariants. An efficient method is used to classify drift orbits, which greatly facilitates identification of the last closed drift path or other drift boundaries. The time-independent evolution of the bounce-averaged phase space density along convective drift orbits is calculated. These three tools can be used to quantitatively describe convective evolution of the particle distribution from the tail, an essential step in understanding the production of unstable distributions in the magnetosphere.

Sheldon, R. B.↗

The Effect of Magnetic Topology on the Escape of Flare Particles

Magnetic reconnection in the solar atmosphere is believed to be the driver of most solar explosive phenomena. Therefore, the topology of the coronal magnetic field is central to understanding the solar drivers of space weather. Of particular importance to space weather are the impulsive Solar Energetic particles that are associated with some CME/eruptive flare events. Observationally, the magnetic configuration of active regions where solar eruptions originate appears to agree with the standard eruptive flare model. According to this model, particles accelerated at the flare reconnection site should remain trapped in the corona and the ejected plasmoid. However, flare-accelerated particles frequently reach the Earth long before the CME does. We present a model that may account for the injection of energetic particles onto open magnetic flux tubes connecting to the Earth. Our model is based on the well-known 2.5D breakout topology, which has a coronal null point (null line) and a four-flux system. A key new addition, however, is that we include an isothermal solar wind with open-flux regions. Depending on the location of the open flux with respect to the null point, we find that the flare reconnection can consist of two distinct phases. At first, the flare reconnection involves only closed field, but if the eruption occurs close to the open field, we find a second phase involving interchange reconnection between open and closed. We argue that this second reconnection episode is responsible for the injection of flare-accelerated particles into the interplanetary medium. We will report on our recent work toward understanding how flare particles escape to the heliosphere. This work uses high-resolution 2.5D MHD numerical simulations performed with the Adaptively Refined MHD Solver (ARMS).

Antiochos, S. K.↗

A Model for SEP Escape

Magnetic reconnection in the solar atmosphere is believed to be the driver of most solar explosive phenomena. Therefore, the topology of the coronal magnetic field is central to understanding the solar drivers of space weather. Of particular importance to space weather are the impulsive Solar Energetic particles that are associated with some CME/eruptive flare events. Observationally, the magnetic configuration of active regions where solar eruptions originate appears to agree with the standard eruptive flare model. According to this model, particles accelerated at the flare reconnection site should remain trapped in the corona and the ejected plasmoid. However, flare-accelerated particles frequently reach the Earth long before the CME does. We present a model that may account for the injection of energetic particles onto open magnetic flux tubes connecting to the Earth. Our model is based on the well-known 2.5D breakout topology, which has a coronal null point (null line) and a four-flux system. A key new addition, however, is that we include an isothermal solar wind with open-flux regions. Depending on the location of the open flux with respect to the null point, we find that the flare reconnection can consist of two distinct phases. At first, the flare reconnection involves only closed field, but if the eruption occurs close to the open field, we find a second phase involving interchange reconnection between open and closed. We argue that this second reconnection episode is responsible for the injection of flare-accelerated particles into the interplanetary medium. We will report on our recent work toward understanding how flare particles escape to the heliosphere. This work uses high-resolution 2.5D MHD numerical simulations performed with the Adaptively Refined MHD Solver (ARMS).

Antiochos, S. K.↗

CFD Analysis of the 24-inch JIRAD Hybrid Rocket Motor

A series of multispecies, multiphase computational fluid dynamics (CFD) analyses of the 24-inch diameter joint government industry industrial research and development (JIRAD) hybrid rocket motor is described. The 24-inch JIRAD hybrid motor operates by injection of liquid oxygen (LOX) into a vaporization plenum chamber upstream of ports in the hydroxyl-terminated polybutadiene (HTPB) solid fuel. The injector spray pattern had a strong influence on combustion stability of the JIRAD motor so a CFD study was initiated to define the injector end flow field under different oxidizer spray patterns and operating conditions. By using CFD to gain a clear picture of the flow field and temperature distribution within the JIRAD motor, it is hoped that the fundamental mechanisms of hybrid combustion instability may be identified and then suppressed by simple alterations to the oxidizer injection parameters such as injection angle and velocity. The simulations in this study were carried out using the General Algorithm for Analysis of Combustion SYstems (GALACSY) multiphase combustion codes. GALACSY consists of a comprehensive set of droplet dynamic submodels (atomization, evaporation, etc.) and a computationally efficient hydrocarbon chemistry package built around a robust Navier-Stokes solver optimized for low Mach number flows. Lagrangian tracking of dispersed particles describes a closely coupled spray phase. The CFD cases described in this paper represent various levels of simplification of the problem. They include: (A) gaseous oxygen with combusting fuel vapor blowing off the walls at various oxidizer injection angles and velocities, (B) gaseous oxygen with combusting fuel vapor blowing off the walls, and (C) liquid oxygen with combusting fuel vapor blowing off the walls. The study used an axisymmetric model and the results indicate that the injector design significantly effects the flow field in the injector end of the motor. Markedly different recirculation patterns are observed in the vaporization chamber as the oxygen velocity and/or spray pattern is varied. The ability of these recirculation patterns to stabilize the diffusion flame above the surface of the solid fuel gives a plausible explanation for the experimentally determined combustion stability characteristics of the JIRAD motor, and suggests how combustion stability can be assured by modifications to the injector design.

Liang, Pak-Yan↗

Buoyancy Driven Shear Flows of Bubble Suspensions

In this work the gas volume fraction and the root-mean-squared fluid velocity are measured in buoyancy driven shear flows of bubble suspensions in a tall, inclined, rectangular channel. The experiments are performed under conditions where We << 1a nd Re >> 1, for which comparisons are made with kinetic theory and numerical simulations. Here Re = gamma(a(exp 2)/nu is the Reynolds number and We = rho(gamma(exp 2))a(exp 3)/sigma is the Weber number; gamma is the shear rate, a is the bubble radius, nu is the kinematic viscosity of the liquid, rho is the density of the liquid, and sigma is the surface tension of the gas/liquid interface. Kang et al. calculated the bubble phase pressure and velocity variance of sheared bubble suspensions under conditions where the bubbles are spherical and the liquid phase velocity field can be approximated using potential flow theory, i.e. We= 0 and Re >> 1. Such conditions can be achieved in an experiment using gas bubbles, with a radius of O(0.5mm), in water. The theory requires that there be no average relative motion of the gas and liquid phases, hence the motivation for an experimental program in microgravity. The necessity of performing preliminary, Earth based experiments, however, requires performing experiments where the gas phase rises in the liquid, which significantly complicates the comparison of experiments with theory. Rather than comparing experimental results with theory for a uniform, homogeneous shear flow, experiments can be compared directly with solutions of the averaged equations of motion for bubble suspensions. This requires accounting for the significant lift force acting on the gas phase when the bubbles rise parallel to the average velocity of the sheared suspension. Shear flows can be produced in which the bubble phase pressure gradient, arising from shear induced collisions amongst the bubbles, balances a body force (centrifugal or gravitational) on the gas phase. A steady, non-uniform gas volume fraction can be measured, from which the bubble phase pressure gradient can be obtained and compared to theory and numerical simulations. The presence of bounding walls further complicates the experiments, since the detailed interactions of the bubbles with bounding walls is not well understood, especially in the presence of gravity, where the momentum and energy exchange depends on the inclination of the wall.

Koch, D. L.↗

Simulation study of Type 2 counterstreaming electrons along auroral field lines

The production of counterstreaming electrons associated with parallel fields along auroral field lines is examined through the use of computer simulation. A 2 1/2-dimensional (two spatial and three velocity dimensions) electrostatic particle algorithm and auroral boundary conditions are used to set up a self-consistent V potential structure. The simulation produces signatures of counterstreaming electrons resembling those observed by the Dynamics Explorer 1 satellite. The main signatures are as follows: (1) the phase space contours of the electron distribution function are elongated along the V-parallel axis, and (2) the energy of electrons streaming in the upward direction is comparable to the energy of the accelerated electron beam. The simulation indicates that a portion of the accelerated electron beam is trapped by large amplitude electrostatic waves produced through the two-stream instability. Strong wave-particle interactions then thermalize the trapped electrons to produce suprathermal electrons streaming in the direction opposite to that of the accelerated electron beam. These results suggest a possible mechanism of producing counterstreaming electron fluxes through nonlinear processes of the two-stream instability.

Wagner, J. S.↗

Progress of High Efficiency Centrifugal Compressor Simulations Using TURBO

Three-dimensional, time-accurate, and phase-lagged computational fluid dynamics (CFD) simulations of the High Efficiency Centrifugal Compressor (HECC) stage were generated using the TURBO solver. Changes to the TURBO Parallel Version 4 source code were made in order to properly model the no-slip boundary condition along the spinning hub region for centrifugal impellers. A startup procedure was developed to generate a converged flow field in TURBO. This procedure initialized computations on a coarsened mesh generated by the Turbomachinery Gridding System (TGS) and relied on a method of systematically increasing wheel speed and backpressure. Baseline design-speed TURBO results generally overpredicted total pressure ratio, adiabatic efficiency, and the choking flow rate of the HECC stage as compared with the design-intent CFD results of Code Leo. Including diffuser fillet geometry in the TURBO computation resulted in a 0.6 percent reduction in the choking flow rate and led to a better match with design-intent CFD. Diffuser fillets reduced annulus cross-sectional area but also reduced corner separation, and thus blockage, in the diffuser passage. It was found that the TURBO computations are somewhat insensitive to inlet total pressure changing from the TURBO default inlet pressure of 14.7 pounds per square inch (101.35 kilopascals) down to 11.0 pounds per square inch (75.83 kilopascals), the inlet pressure of the component test. Off-design tip clearance was modeled in TURBO in two computations: one in which the blade tip geometry was trimmed by 12 mils (0.3048 millimeters), and another in which the hub flow path was moved to reflect a 12-mil axial shift in the impeller hub, creating a step at the hub. The one-dimensional results of these two computations indicate non-negligible differences between the two modeling approaches.

turbomachinery↗

Auroral, Ionospheric and Ground Magnetic Signatures of Magnetopause Surface Modes

Surface waves on Earth's magnetopause have a controlling effect upon global magnetospheric dynamics. Since spacecraft provide sparse in situ observation points, remote sensing these modes using ground-based instruments in the polar regions is desirable. However, many open conceptual questions on the expected signatures remain. Therefore, we provide predictions of key qualitative features expected in auroral, ionospheric, and ground magnetic observations through both magnetohydrodynamic theory and a global coupled magnetosphere-ionosphere simulation of a magnetopause surface eigenmode. These show monochromatic oscillatory field-aligned currents (FACs), due to both the surface mode and its non-resonant Alfvén coupling, are present throughout the magnetosphere. The currents peak in amplitude at the equatorward edge of the magnetopause boundary layer, not the open-closed boundary as previously thought. They also exhibit slow poleward phase motion rather than being purely evanescent. We suggest the upward FAC perturbations may result in periodic auroral brightenings. In the ionosphere, convection vortices circulate the poleward moving FAC structures. Finally, surface mode signals are predicted in the ground magnetic field, with ionospheric Hall currents rotating perturbations by approximately (but not exactly) 90° compared to the magnetosphere. Thus typical dayside magnetopause surface modes should be strongest in the East-West ground magnetic field component. Overall, all ground-based signatures of the magnetopause surface mode are predicted to have the same frequency across L-shells, amplitudes that maximize near the magnetopause's equatorward edge, and larger latitudinal scales than for field line resonance. Implications in terms of ionospheric Joule heating and geomagnetically induced currents are discussed.

M. O. Archer↗

Numerical Simulations and Diagnostic Studies Relating Meteorology To Atmospheric Chemistry during TRACE-P

The Florida State University (FSU) team participated extensively in the pre-mission planning for TRACE-P through meetings, telephone calls, and e-mails. During Spring 2001, Prof. Fuelberg served as DC-8 Mission Meteorologist during the field campaign. He prepared meteorological guidance for each flight of the DC-8 and flew on each mission. After the field phase, FSU prepared various meteorological products, included backward air trajectories, for each flight of the DC-8 and P-3B. These were posted on the FSU and NASA-GTE web sites for use by all the Science Team. During the two-year post mission period, FSU conducted research relating meteorology to atmospheric chemistry during TRACE-P. This led to three journal articles in the Journal of Geophysical Research. FSU personnel were the lead authors on each of these articles. Abstracts of these articles are attached. In addition, the FSU team collaborated with other members of the TRACE-P Science Team to incorporate meteorological factors into their research. A list of publications resulting from these interactions is included.

Fuelberg, Henry E.↗

Biogeochemical Protocols and Diagnostics for the CMIP6 Ocean Model Intercomparison Project (OMIP)

The Ocean Model Intercomparison Project (OMIP) focuses on the physics and biogeochemistry of the ocean component of Earth system models participating in the sixth phase of the Coupled Model Intercomparison Project (CMIP6). OMIP aims to provide standard protocols and diagnostics for ocean models, while offering a forum to promote their common assessment and improvement. It also offers to compare solutions of the same ocean models when forced with reanalysis data (OMIP simulations) vs. when integrated within fully coupled Earth system models (CMIP6). Here we detail simulation protocols and diagnostics for OMIP's biogeochemical and inert chemical tracers. These passive-tracer simulations will be coupled to ocean circulation models, initialized with observational data or output from a model spin-up, and forced by repeating the 1948-2009 surface fluxes of heat, fresh water, and momentum. These so-called OMIP-BGC simulations include three inert chemical tracers (CFC-11, CFC-12, SF [subscript] 6) and biogeochemical tracers (e.g., dissolved inorganic carbon, carbon isotopes, alkalinity, nutrients, and oxygen). Modelers will use their preferred prognostic BGC model but should follow common guidelines for gas exchange and carbonate chemistry. Simulations include both natural and total carbon tracers. The required forced simulation (omip1) will be initialized with gridded observational climatologies. An optional forced simulation (omip1-spunup) will be initialized instead with BGC fields from a long model spin-up, preferably for 2000 years or more, and forced by repeating the same 62-year meteorological forcing. That optional run will also include abiotic tracers of total dissolved inorganic carbon and radiocarbon, CTabio and 14CTabio, to assess deep-ocean ventilation and distinguish the role of physics vs. biology. These simulations will be forced by observed atmospheric histories of the three inert gases and CO2 as well as carbon isotope ratios of CO2. OMIP-BGC simulation protocols are founded on those from previous phases of the Ocean Carbon-Cycle Model Intercomparison Project. They have been merged and updated to reflect improvements concerning gas exchange, carbonate chemistry, and new data for initial conditions and atmospheric gas histories. Code is provided to facilitate their implementation.

Carbon↗

The flight and impact dispersion of the Comet Nucleus Penetrator - A six-degree-of-freedom simulation

The flight dynamics of a Penetrator probe which, as part of a Comet Rendezvous and Asteroid Flyby mission will impact a known periodic comet, possibly Hempel 2, and penetrate beneath the surface to obtain measurements, is discussed. The mission is described, including the penetrator-separation dynamics, the dynamics during initial coast and rocket thrust phase, and the final coast to impact. The effects on the results that may be caused by the comet nucleus gravity field variations, uncertainty in the spacecraft state vector, the spacecraft pointing accuracy, initial separation errors, and rocket motor thrust errors are examined.

Mascy, A. C.↗

NASA Tech Briefs, July 2009

Topics covered include: Dual Cryogenic Capacitive Density Sensor; Hail Monitor Sensor; Miniature Six-Axis Load Sensor for Robotic Fingertip; Improved Blackbody Temperature Sensors for a Vacuum Furnace; Wrap-Around Out-the-Window Sensor Fusion System; Wide-Range Temperature Sensors with High-Level Pulse Train Output; Terminal Descent Sensor Simulation; A Robust Mechanical Sensing System for Unmanned Sea Surface Vehicles; Additive for Low-Temperature Operation of Li-(CF)n Cells; Li/CFx Cells Optimized for Low-Temperature Operation; Number Codes Readable by Magnetic-Field-Response Recorders; Determining Locations by Use of Networks of Passive Beacons; Superconducting Hot-Electron Submillimeter-Wave Detector; Large-Aperture Membrane Active Phased-Array Antennas; Optical Injection Locking of a VCSEL in an OEO; Measuring Multiple Resistances Using Single-Point Excitation; Improved-Bandwidth Transimpedance Amplifier; Inter-Symbol Guard Time for Synchronizing Optical PPM; Novel Materials Containing Single-Wall Carbon Nanotubes Wrapped in Polymer Molecules; Light-Curing Adhesive Repair Tapes; Thin-Film Solid Oxide Fuel Cells; Zinc Alloys for the Fabrication of Semiconductor Devices; Small, Lightweight, Collapsible Glove Box; Radial Halbach Magnetic Bearings; Aerial Deployment and Inflation System for Mars Helium Balloons; Steel Primer Chamber Assemblies for Dual Initiated Pyrovalves; Voice Coil Percussive Mechanism Concept for Hammer Drill; Inherently Ducted Propfans and Bi-Props; Silicon Nanowire Growth at Chosen Positions and Orientations; Detecting Airborne Mercury by Use of Gold Nanowires; Detecting Airborne Mercury by Use of Palladium Chloride; Micro Electron MicroProbe and Sample Analyzer; Nanowire Electron Scattering Spectroscopy; Electron-Spin Filters Would Offer Spin Polarization Greater than 1; Subcritical-Water Extraction of Organics from Solid Matrices; A Model for Predicting Thermoelectric Properties of Bi2Te3; Integrated Miniature Arrays of Optical Biomolecule Detectors; A Software Rejuvenation Framework for Distributed Computing; Kurtosis Approach to Solution of a Nonlinear ICA Problem; Robust Software Architecture for Robots; R4SA for Controlling Robots; Bio-Inspired Neural Model for Learning Dynamic Models; Evolutionary Computing Methods for Spectral Retrieval; Monitoring Disasters by Use of Instrumented Robotic Aircraft; Complexity for Survival of Living Systems; Using Drained Spacecraft Propellant Tanks for Habitation; Connecting Node; and Electrolytes for Low-Temperature Operation of Li-CFx Cells.

Source record↗

Delay and Disruption Tolerant Networking MACHETE Model

To verify satisfaction of communication requirements imposed by unique missions, as early as 2000, the Communications Networking Group at the Jet Propulsion Laboratory (JPL) saw the need for an environment to support interplanetary communication protocol design, validation, and characterization. JPL's Multi-mission Advanced Communications Hybrid Environment for Test and Evaluation (MACHETE), described in Simulator of Space Communication Networks (NPO-41373) NASA Tech Briefs, Vol. 29, No. 8 (August 2005), p. 44, combines various commercial, non-commercial, and in-house custom tools for simulation and performance analysis of space networks. The MACHETE environment supports orbital analysis, link budget analysis, communications network simulations, and hardware-in-the-loop testing. As NASA is expanding its Space Communications and Navigation (SCaN) capabilities to support planned and future missions, building infrastructure to maintain services and developing enabling technologies, an important and broader role is seen for MACHETE in design-phase evaluation of future SCaN architectures. To support evaluation of the developing Delay Tolerant Networking (DTN) field and its applicability for space networks, JPL developed MACHETE models for DTN Bundle Protocol (BP) and Licklider/Long-haul Transmission Protocol (LTP). DTN is an Internet Research Task Force (IRTF) architecture providing communication in and/or through highly stressed networking environments such as space exploration and battlefield networks. Stressed networking environments include those with intermittent (predictable and unknown) connectivity, large and/or variable delays, and high bit error rates. To provide its services over existing domain specific protocols, the DTN protocols reside at the application layer of the TCP/IP stack, forming a store-and-forward overlay network. The key capabilities of the Bundle Protocol include custody-based reliability, the ability to cope with intermittent connectivity, the ability to take advantage of scheduled and opportunistic connectivity, and late binding of names to addresses.

Segui, John S.↗

Comparison of INP Parameterizations for Dust Minerals in Climatological Simulations With a Global Model

The effect of aerosol particles on ice nucleation and, in turn, the formation of ice and mixed phase clouds is recognized as one of the largest sources of uncertainty in weather and climate prediction. We utilize an improved sectional dust module in NASA GISS Earth System ModelE2.1, which distinguishes eight different mineral species and accretions between iron oxides and the other minerals. Simulations over a period of 20 years have been carried out with this model, and the mineral fields and other model variables (temperature, relative humidity) are used to calculate the ice nucleating particle (INP) number concentration, applying time-independent and time-dependent INP parameterizations, such as active site parameterization and water activity based immersion freezing model (ABIFM). We study how the dependence of the parameterizations on different model variables affects the mean INP number concentration. The sensitivity of the INP number concentration to fundamental dust properties such as emitted mineral size distributions and mixing state between minerals is also investigated. Results show that the sensitivity of the total INP number concentration to the emitted dust size distribution is rather small, but the sensitivity over the whole size range obscures offsetting differences in the magnitude and the sign of the sensitivity between smaller and larger particles.

INP Parameterizations↗

Broadband Performance of TPF's High-contrast Imaging Testbed: Modeling and Simulations

The broadband performance of the high-contrast imaging testbed (HCIT) at JPL is investigated through optical modeling and simulations. The analytical tool is an optical simulation algorithm developed by combining the HCIT's optical model with a speckle-nulling algorithm that operates directly on coronagraphic images, an algorithm identical to the one currently being used on the HCIT to actively suppress scattered light via a deformable mirror. It is capable of performing full three-dimensional end-to-end near-field diffraction analysis on the HCIT's optical system. By conducting speckle-nulling optimization, we clarify the HCIT's capability and limitations in terms of its broadband contrast performance under various realistic conditions. Considered cases include non-ideal occulting masks, such as a mask with optical density and wavelength dependent parasitic phase-delay errors (i.e., a not band-limited occulting mask) and the one with an optical-density profile corresponding to a measured, non-standard profile, as well as the independently measured phase errors of all optics. Most of the information gathered on the HCIT's optical components through measurement and characterization over the last several years at JPL has been used in this analysis to make the predictions as accurate as possible. The best contrast values predicted so far by our simulations obtainable on the HCIT illuminated with a broadband light having a bandwidth of 80nm and centered at 800nm wavelength are Cm=1.1x10-8 (mean) and C4=4.9x10-8 (at 4(lamda)/D), respectively. In this paper we report our preliminary findings about the broadband light performance of the HCIT.

integrated modeling↗

Li-Doped Ionic Liquid Electrolytes: From Bulk Phase to Interfacial Behavior

Ionic liquids have been proposed as candidate electrolytes for high-energy density, rechargeable batteries. We present an extensive computational analysis supported by experimental comparisons of the bulk and interfacial properties of a representative set of these electrolytes as a function of Li-salt doping. We begin by investigating the bulk electrolyte using quantum chemistry and ab initio molecular dynamics to elucidate the solvation structure of Li(+). MD simulations using the polarizable force field of Borodin and coworkers were then performed, from which we obtain an array of thermodynamic and transport properties. Excellent agreement is found with experiments for diffusion, ionic conductivity, and viscosity. Combining MD simulations with electronic structure computations, we computed the electrochemical window of the electrolytes across a range of Li(+)-doping levels and comment on the role of the liquid environment. Finally, we performed a suite of simulations of these Li-doped electrolytes at ideal electrified interfaces to evaluate the differential capacitance and the equilibrium Li(+) distribution in the double layer. The magnitude of differential capacitance is in good agreement with our experiments and exhibits the characteristic camel-shaped profile. In addition, the simulations reveal Li(+) to be highly localized to the second molecular layer of the double layer, which is supported by additional computations that find this layer to be a free energy minimum with respect to Li(+) translation.

molecular dynamics↗