Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “grid stability”

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 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Improving the geological interpretation of magnetic and gravity satellite anomalies

Quantitative analysis of the geologic component of observed satellite magnetic and gravity fields requires accurate isolation of the geologic component of the observations, theoretically sound and viable inversion techniques, and integration of collateral, constraining geologic and geophysical data. A number of significant contributions were made which make quantitative analysis more accurate. These include procedures for: screening and processing orbital data for lithospheric signals based on signal repeatability and wavelength analysis; producing accurate gridded anomaly values at constant elevations from the orbital data by three-dimensional least squares collocation; increasing the stability of equivalent point source inversion and criteria for the selection of the optimum damping parameter; enhancing inversion techniques through an iterative procedure based on the superposition theorem of potential fields; and modeling efficiently regional-scale lithospheric sources of satellite magnetic anomalies. In addition, these techniques were utilized to investigate regional anomaly sources of North and South America and India and to provide constraints to continental reconstruction. Since the inception of this research study, eleven papers were presented with associated published abstracts, three theses were completed, four papers were published or accepted for publication, and an additional manuscript was submitted for publication.

Hinze, William J.↗
NASA NTRS1997

Evaluating and Understanding Parameterized Convective Processes and their Role in the Development of Mesoscale Precipitation Systems

Research efforts during the second year have centered on improving the manner in which convective stabilization is achieved in the Penn State/NCAR mesoscale model MM5. Ways of improving this stabilization have been investigated by (1) refining the partitioning between the Kain-Fritsch convective parameterization scheme and the grid scale by introducing a form of moist convective adjustment; (2) using radar data to define locations of subgrid-scale convection during a dynamic initialization period; and (3) parameterizing deep-convective feedbacks as subgrid-scale sources and sinks of mass. These investigations were conducted by simulating a long-lived convectively-generated mesoscale vortex that occurred during 14-18 Jul. 1982 and the 10-11 Jun. 1985 squall line that occurred over the Kansas-Oklahoma region during the PRE-STORM experiment. The long-lived vortex tracked across the central Plains states and was responsible for multiple convective outbreaks during its lifetime.

Fritsch, J. Michael↗
NASA NTRS1998

Comments on the Diffusive Behavior of Two Upwind Schemes

The diffusive characteristics of two upwind schemes, multi-dimensional fluctuation splitting and locally one-dimensional finite volume, are compared for scalar advection-diffusion problems. Algorithms for the two schemes are developed for node-based data representation on median-dual meshes associated with unstructured triangulations in two spatial dimensions. Four model equations are considered: linear advection, non-linear advection, diffusion, and advection-diffusion. Modular coding is employed to isolate the effects of the two approaches for upwind flux evaluation, allowing for head-to-head accuracy and efficiency comparisons. Both the stability of compressive limiters and the amount of artificial diffusion generated by the schemes is found to be grid-orientation dependent, with the fluctuation splitting scheme producing less artificial diffusion than the finite volume scheme. Convergence rates are compared for the combined advection-diffusion problem, with a speedup of 2.5 seen for fluctuation splitting versus finite volume when solved on the same mesh. However, accurate solutions to problems with small diffusion coefficients can be achieved on coarser meshes using fluctuation splitting rather than finite volume, so that when comparing convergence rates to reach a given accuracy, fluctuation splitting shows a speedup of 29 over finite volume.

Wood, William A.↗
NASA NTRS1998

Diffusion Characteristics of Upwind Schemes on Unstructured Triangulations

The diffusive characteristics of two upwind schemes, multi-dimensional fluctuation splitting and dimensionally-split finite volume, are compared for scalar advection-diffusion problems. Algorithms for the two schemes are developed for node-based data representation on median-dual meshes associated with unstructured triangulations in two spatial dimensions. Four model equations are considered: linear advection, non-linear advection, diffusion, and advection-diffusion. Modular coding is employed to isolate the effects of the two approaches for upwind flux evaluation, allowing for head-to-head accuracy and efficiency comparisons. Both the stability of compressive limiters and the amount of artificial diffusion generated by the schemes is found to be grid-orientation dependent, with the fluctuation splitting scheme producing less artificial diffusion than the dimensionally-split finite volume scheme. Convergence rates are compared for the combined advection-diffusion problem, with a speedup of 2-3 seen for fluctuation splitting versus finite volume when solved on the same mesh. However, accurate solutions to problems with small diffusion coefficients can be achieved on coarser meshes using fluctuation splitting rather than finite volume, so that when comparing convergence rates to reach a given accuracy, fluctuation splitting shows a 20-25 speedup over finite volume.

Wood, William A.↗
NASA NTRS2021

Drought Relief and Reversal Over North America from 1500 to 2016

Season-to-season persistence of soil moisture drought varies across North America. Such inter-seasonal autocorrelation can have modest skill in forecasting future conditions several months in advance. Because robust instrumental observations of precipitation span less than 100 years, the temporal stability of the relationship between seasonal moisture anomalies is uncertain. The North American Seasonal Precipitation Atlas (NASPA) is a gridded network of separately reconstructed cool (December-April) and warm (May-July) season precipitation series and offers new insights on the intra-annual changes in drought for up to 2,000 years. Here, the NASPA precipitation reconstructions are rescaled to represent the long-term soil moisture balance during the cool season and three-month long atmospheric moisture during the warm season. These rescaled seasonal reconstructions are then used to quantify the frequency, magnitude, and spatial extent of cool season drought that was relieved or reversed during the following summer months. The adjusted seasonal reconstructions reproduce the general patterns of large-scale drought amelioration and termination in the instrumental record during the 20th century and are used to estimate relief and reversals for the most skillfully reconstructed past 500-years. Sub-continental to continental scale reversals of cool season drought in the following warm season have been rare, but the reconstructions display periods prior to the instrumental data of increased reversal probabilities for the Mid-Atlantic region and the US Southwest. Drought relief at the continental scale may arise in part from macroscale ocean-atmosphere processes, while the smaller scale regional reversals may reflect land-surface feedbacks and stochastic variability.

Drought↗
DOE OSTI2013

Global analysis of night marine air temperature and its uncertainty since 1880

An updated version of the Met Office Hadley Centre's monthly night marine air temperature data set is presented. It is available on a 5° latitude-longitude grid from 1880 as anomalies relative to 1961–1990 calendar-monthly climatological average night marine air temperature (NMAT). Adjustments are made for changes in observation height; these depend on estimates of the stability of the near surface atmospheric boundary layer. In previous versions of the data set, ad hoc adjustments were also made for three periods and regions where poor observational practice was prevalent. These adjustments are re-examined. Estimates of uncertainty are calculated for every grid box and result from measurement errors, uncertainty in adjustments applied to the observations, uncertainty in the measurement height, and under-sampling. The new data set is a clear improvement over previous versions in terms of coverage because of the recent digitization of historical observations from ships' logbooks. However, the periods prior to about 1890 and around World War II remain particularly uncertain, and sampling is still sparse in some regions in other periods. A further improvement is the availability of uncertainty estimates for every grid box and every month. Previous versions required adjustments that were dependent on contemporary measurements of sea surface temperature (SST); to avoid these, the new data set starts in 1880 rather than 1856. Overall agreement with variations of SST is better for the updated data set than for previous versions, supporting existing estimates of global warming and increasing confidence in the global record of temperature variability and change.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1995

On the Importance of the Dynamics of Discretizations

It has been realized recently that the discrete maps resulting from numerical discretizations of differential equations can possess asymptotic dynamical behavior quite different from that of the original systems. This is the case not only for systems of Ordinary Differential Equations (ODEs) but in a more complicated manner for Partial Differential Equations (PDEs) used to model complex physics. The impact of the modified dynamics may be mild and even not observed for some numerical methods. For other classes of discretizations the impact may be pronounced, but not always obvious depending on the nonlinear model equations, the time steps, the grid spacings and the initial conditions. Non-convergence or convergence to periodic solutions might be easily recognizable but convergence to incorrect but plausible solutions may not be so obvious - even for discretized parameters within the linearized stability constraint. Based on our past four years of research, we will illustrate some of the pathology of the dynamics of discretizations, its possible impact and the usage of these schemes for model nonlinear ODEs, convection-diffusion equations and grid adaptations.

Sweby, Peter K.↗
DOE OSTI2023

Efficient learning of power grid voltage control strategies via model-based deep reinforcement learning

Here this article proposes a model-based deep reinforcement learning (DRL) method to design emergency control strategies for short-term voltage stability problems in power systems. Recent advances show promising results for model-free DRL-based methods in power systems control problems. But in power systems applications, these model-free methods have certain issues related to training time (clock time) and sample efficiency; both are critical for making state-of-the-art DRL algorithms practically applicable. DRL-agent learns an optimal policy via a trial-and-error method while interacting with the real-world environment. It is also desirable to minimize the direct interaction of the DRL agent with the real-world power grid due to its safety-critical nature. Additionally, the state-of-the-art DRL-based policies are mostly trained using a physics-based grid simulator where dynamic simulation is computationally intensive, lowering the training efficiency. We propose a novel model-based DRL framework where a deep neural network (DNN)-based dynamic surrogate model (SM), instead of a real-world power grid or physics-based simulation, is utilized within the policy learning framework, making the process faster and more sample efficient. However, having stable training in model-based DRL is challenging because of the complex system dynamics of large-scale power systems. We addressed these issues by incorporating imitation learning to have a warm start in policy learning, reward-shaping, and multi-step loss in surrogate model training. Finally, we achieved 97.5% reduction in samples and 87.7% reduction in training time for an application to the IEEE 300-bus test system.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2011

Phase-field investigation on the non-equilibrium interface dynamics of rapid alloy solidification

The research program reported here is focused on critical issues that represent conspicuous gaps in current understanding of rapid solidification, limiting our ability to predict and control microstructural evolution (i.e. morphological dynamics and microsegregation) at high undercooling, where conditions depart significantly from local equilibrium. More specifically, through careful application of phase-field modeling, using appropriate thin-interface and anti-trapping corrections and addressing important details such as transient effects and a velocity-dependent (i.e. adaptive) numerics, the current analysis provides a reasonable simulation-based picture of non-equilibrium solute partitioning and the corresponding oscillatory dynamics associated with single-phase rapid solidification and show that this method is a suitable means for a self-consistent simulation of transient behavior and operating point selection under rapid growth conditions. Moving beyond the limitations of conventional theoretical/analytical treatments of non-equilibrium solute partitioning, these results serve to substantiate recent experimental findings and analytical treatments for single-phase rapid solidification. The departure from the equilibrium solid concentration at the solid-liquid interface was often observed during rapid solidification, and the energetic associated non-equilibrium solute partitioning has been treated in detail, providing possible ranges of interface concentrations for a given growth condition. Use of these treatments for analytical description of specific single-phase dendritic and cellular operating point selection, however, requires a model for solute partitioning under a given set of growth conditions. Therefore, analytical solute trapping models which describe the chemical partitioning as a function of steady state interface velocities have been developed and widely utilized in most of the theoretical investigations of rapid solidification. However, these solute trapping models are not rigorously verified due to the difficulty in experimentally measuring under rapid growth conditions. Moreover, since these solute trapping models include kinetic parameters which are difficult to directly measure from experiments, application of the solute trapping models or the associated analytic rapid solidification model is limited. These theoretical models for steady state rapid solidification which incorporate the solute trapping models do not describe the interdependency of solute diffusion, interface kinetics, and alloy thermodynamics. The phase-field approach allows calculating, spontaneously, the non-equilibrium growth effects of alloys and the associated time-dependent growth dynamics, without making the assumptions that solute partitioning is an explicit function of velocity, as is the current convention. In the research described here, by utilizing the phase-field model in the thin-interface limit, incorporating the anti-trapping current term, more quantitatively valid interface kinetics and solute diffusion across the interface are calculated. In order to sufficiently resolve the physical length scales (i.e. interface thickness and diffusion boundary length), grid spacings are continually adjusted in calculations. The full trajectories of transient planar growth dynamics under rapid directional solidification conditions with different pulling velocities are described. As a validation of a model, the predicted steady state conditions are consistent with the analytic approach for rapid growth. It was confirmed that rapid interface dynamics exhibits the abrupt acceleration of the planar front when the effect of the non-equilibrium solute partitioning at the interface becomes signi ficant. This is consistent with the previous linear stability analysis for the non-equilibrium interface dynamics. With an appropriate growth condition, the continuous oscillation dynamics was able to be simulated using continually adjusting grid spacings. This oscillatory dynamics including instantaneous jump of interface velocities are consistent with a previous phenomenological model by and a numerical investigation, which may cause the formation of banded structures. Additionally, the selection of the steady state growth dynamics in the highly undercooled melt is demonstrated. The transition of the growth morphology, interface velocity selection, and solute trapping phenomenon with increasing melt supersaturations was described by the phase-field simulation. The tip selection for the dendritic growth was consistent with Ivantsov's function, and the non-equilibrium chemical partitioning behavior shows good qualitative agreement with the Aziz's solute trapping model even though the model parameter(V D ) remains as an arbitrary constant. This work is able to show the possibility of comprehensive description of rapid alloy growth over the entire time-dependent non-equilibrium phenomenon.

36 MATERIALS SCIENCE↗
NASA NTRS2001

Effect of Under-Resolved Grids on High Order Methods

There has been much discussion on verification and validation processes for establishing the credibility of CFD simulations. Since the early 1990s, many of the aeronautical and mechanical engineering related reference journals mandated that any accepted articles in numerical simulations (without known solutions to compared with) need to perform a minimum of one level of grid refinement and time step reduction. Due to the difficulty in analysis, the effect of under-resolved grids and the nonlinear behavior of available spatial discretizations, are scarcely discussed in the literature. Here, an under-resolved numerical simulation is one where the grid spacing being used is too coarse to resolve the smallest physically relevant scales of the chosen continuum governing equations that are of interest to the numerical modeler. With the advent of new developments in fourth-order or higher spatial schemes, it has become common to regard high order schemes as more accurate, reliable and require less grid points. The danger comes when one tries to perform computations with the coarsest grid possible while still hoping to maintain numerical results sufficiently accurate for complex flows, and especially, data-limited problems. On one hand, high order methods when applies to highly coupled multidimensional complex nonlinear problems might have different stability, convergence and reliability behavior than their well studied low order counterparts, especially for nonlinear schemes such as TVD, MUSCL with limiters, ENO, WENO and discrete Galerkin. On the other hand, high order methods involve more operation counts and systematic grid convergence study can be time consuming and prohibitively expansive. At the same time it is difficult to fully understand or categorize the different nonlinear behavior of finite discretizations, especially at the limits of under-resolution when different types of bifurcation phenomena might occur, depending on the combination of grid spacings, time steps, initial conditions and numerical treatments of boundary conditions.

Yee, H. C.↗
DOE OSTI2015

Novel Sensor for the In Situ Measurement of Uranium Fluxes

The goal of this project was to develop a sensor that incorporates the field-tested concepts of the passive flux meter to provide direct in situ measures of flux for uranium and groundwater in porous media. Measurable contaminant fluxes [J] are essentially the product of concentration [C] and groundwater flux or specific discharge [q ]. The sensor measures [J] and [q] by changes in contaminant and tracer amounts respectively on a sorbent. By using measurement rather than inference from static parameters, the sensor can directly advance conceptual and computational models for field scale simulations. The sensor was deployed in conjunction with DOE in obtaining field-scale quantification of subsurface processes affecting uranium transport (e.g., advection) and transformation (e.g., uranium attenuation) at the Rifle IFRC Site in Rifle, Colorado. Project results have expanded our current understanding of how field-scale spatial variations in fluxes of uranium, groundwater and salient electron donor/acceptors are coupled to spatial variations in measured microbial biomass/community composition, effective field-scale uranium mass balances, attenuation, and stability. The coupling between uranium, various nutrients and micro flora can be used to estimate field-scale rates of uranium attenuation and field-scale transitions in microbial communities. This research focuses on uranium (VI), but the sensor principles and design are applicable to field-scale fate and transport of other radionuclides. Laboratory studies focused on sorbent selection and calibration, along with sensor development and validation under controlled conditions. Field studies were conducted at the Rifle IFRC Site in Rifle, Colorado. These studies were closely coordinated with existing SBR (formerly ERSP) projects to complement data collection. Small field tests were conducted during the first two years that focused on evaluating field-scale deployment procedures and validating sensor performance under controlled field conditions. In the third and fourth year a suite of larger field studies were conducted. For these studies, the uranium flux sensor was used with uranium speciation measurements and molecular-biological tools to characterize microbial community and active biomass at synonymous wells distributed in a large grid. These field efforts quantified spatial changes in uranium flux and field-scale rates of uranium attenuation (ambient and stimulated), uranium stability, and quantitatively assessed how fluxes and effective reaction rates were coupled to spatial variations in microbial community and active biomass. Analyses of data from these field experiments were used to generate estimates of Monod kinetic parameters that are ‘effective’ in nature and optimal for modeling uranium fate and transport at the field-scale. This project provided the opportunity to develop the first sensor that provides direct measures of both uranium (VI) and groundwater flux. A multidisciplinary team was assembled to include two geochemists, a microbiologist, and two quantitative contaminant hydrologists. Now that the project is complete, the sensor can be deployed at DOE sites to evaluate field-scale uranium attenuation, source behavior, the efficacy of remediation, and off-site risk. Because the sensor requires no power, it can be deployed at remote sites for periods of days to months. The fundamental science derived from this project can be used to advance the development of predictive models for various transport and attenuation processes in aquifers. Proper development of these models is critical for long-term stewardship of contaminated sites in the context of predicting uranium source behavior, remediation performance, and off-site risk.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Predicting the potentials, solubilities and stabilities of metal-acetylacetonates for non-aqueous redox flow batteries using density functional theory calculations

New active materials are needed to improve the performance and reduce the cost of non-aqueous redox flow batteries (RFBs) for grid-scale energy storage applications. Efforts to develop better performing materials, which have largely been empirical, would benefit from a better understanding of relationships between structural, electronic and RFB-relevant functional properties. This paper focuses on metal-acetylacetonates, a class of metal coordination complexes that has shown promise for use in RFBs, and describes correlations between their experimentally measured standard potentials, solubilities, and stabilities (cycle lifes), and selected chemical, structural and electronic properties determined from Density Functional Theory (DFT) calculations. The training set consisted of 16 complexes including 5 different metals and 11 different substituents on the acetylacetonate ligand. Standard potentials for those compounds were calculated and are in good agreement with experimentally measured results. A predictive equation based on the solvation energies and dipole moments, two easily computed properties, reasonably modeled the experimentally determined solubilities. Importantly, we were able to identify a descriptor for the stability of acetylacetonates. The experimentally determined stability, quantified as the cycle life to a given degree of degradation, correlated with the percentage of the highest occupied (HOMO) or lowest unoccupied molecular orbital (LUMO) on the metal of the complex. This percentage is influenced by the degree of ligand innocence (irreducibility), and complexes with the most innocent ligands yielded the most stable redox reactions. To this end, VO(acetylacetonate)(2) and Fe(acetylacetonate)(3), with nearly 80% of the HOMO and LUMO on the metal, possessed the most stable oxidation and reduction half-reactions, respectively. The structure-function relationships and correlations presented in this paper could be used to predict new, highly soluble and stable complexes for RFB applications.

Kucharyson, J. F.↗
NASA NTRS1993

Fuel optimal reorientation of axisymmetric spin-stabilized satellites

The fuel optimal spin axis reorientation problem for large and small maneuver times is solved using a new numerical technique. The approach, which employs the adaptive grid bisection search, is shown to be useful in determining the fuel optimal impulsive control strategy, particularly when rapid maneuvers are desired. For large enough maneuver times, the method produces the two-impulse reorientation maneuver. For even larger maneuver times, the existence of multiple optimal pulsing schemes is demonstrated.

Redmond, Jim↗
DOE OSTI2025

Co-Simulation of PSS/E, OpenDSS, and PSCAD for Power Systems Stability Analysis With Inverter-Based Resources: Preprint

The increasing penetration of inverter-based resources (IBRs) is reshaping the dynamic behavior of power systems. IEEE standard 1547-2018 suggests that distributed energy resources (DERs) should provide grid services such as voltage and frequency supports. On the other hand, dynamic events caused by IBRs such as sub-synchronous oscillation have been reported. These developments necessitate improvement in the current modeling capabilities to better understand the interdependencies within power systems. These include interactions between transmission and distribution systems, among IBRs themselves, and between IBRs and conventional resources. In this paper, we present a co-simulation model integrating PSS/E, OpenDSS, and PSCAD to analyze IBR impacts on the stability of transmission and distribution systems. A key challenge in developing a co-simulation model is ensuring interoperability among different simulators (interfacing and data flow) while maintaining accurate results. Using the developed model, we simulate the impact of IBRs on power systems in two test cases: 1) fault ride through (FRT) capability during a generation trip contingency; 2) IBR-induced sub-synchronous oscillation. The results show the effectiveness of the co-simulation model.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2016

Stem Inc. SunShot Incubator Program Final Technical Report

In this Energy Storage Control Algorithms project, Stem sought to develop tools and control algorithms to increase the value and reduce balance-of-system and grid integration costs associated with adding distributed solar generation to the grid. These advances fell under the headings SolarScope and SolarController. Stem sought to create initial market traction with a fully commercialized product for the solar industry to size storage systems (SolarScope) as well as a solar intermittency-mitigation framework for utilities (SolarController) in the course of the project. The company sought to align strategic growth plans and enable the rollout of the products to broader audiences in multiple geographic regions by leveraging the major solar companies in the national market as partners. Both final products were both intended to be commercialized. They are: SolarScope: Analysis tool to identify viable PV + storage projects and thereby expedite the sales and interconnection processes. SolarScope combines customer load data, PV production estimates, utility rate tariff, and simulated storage into a simple user interface for PV developers. Developers can easily identify viable solar + storage sites without the need for complex and time consuming, site-by-site spreadsheet modeling. SolarContoller: Tool to autonomously dispatch distributed storage in order to mitigate voltage fluctuation and reduce curtailment. SolarController co-optimizes, in real time, storage dispatch for circuit stability and curtailment reduction, enabling higher penetrations of PV. SolarController is automated, not requiring utility dispatch or management, as Stem hardware senses grid voltage, frequency, customer load, PV production, and power factor. In the end the two products met with different outcomes. SolarScope was tested by potential users, and continues to be used as a foundational platform for partnership with key solar industry partners. SolarController, on the other hand, was successful in lab testing but was not commercialized due to a lack of marketability and lack of interested customer base. Together the development of these two products marked a material step forward for Stem; and a new milestone along the pathway of integration for the solar and storage industries. SolarScope is leading to real, out-of-the-lab project development in storage + solar for the commercial customer sector. Meanwhile SolarController has opened the eyes of regulators and utility executives alike to the potential of distributed solar and by doing so, has moved the conversation forward for the integration of distributed energy resources more broadly on the grid.

25 ENERGY STORAGE↗
DOE OSTI2020

Analysis of the SBP-SAT Stabilization for Finite Element Methods Part I: Linear Problems

In the hyperbolic community, discontinuous Galerkin (DG) approaches are mainly applied when finite element methods are considered. As the name suggested, the DG framework allows a discontinuity at the element interfaces, which seems for many researchers a favorable property in case of hyperbolic balance laws. On the contrary, continuous Galerkin methods appear to be unsuitable for hyperbolic problems and there exists still the perception that continuous Galerkin methods are notoriously unstable. To remedy this issue, stabilization terms are usually added and various formulations can be found in the literature. However, this perception is not true and the stabilization terms are unnecessary, in general. In this paper, we deal with this problem, but present a different approach. We use the boundary conditions to stabilize the scheme following a procedure that are frequently used in the finite difference community. Here, the main idea is to impose the boundary conditions weakly and specific boundary operators are constructed such that they guarantee stability. This approach has already been used in the discontinuous Galerkin framework, but here we apply it with a continuous Galerkin scheme. No internal dissipation is needed even if unstructured grids are used. Further, we point out that we do not need exact integration, it suffices if the quadrature rule and the norm in the differential operator are the same, such that the summation-by-parts property is fulfilled meaning that a discrete Gauss Theorem is valid. This contradicts the perception in the hyperbolic community that stability issues for pure Galerkin scheme exist. In numerical simulations, we verify our theoretical analysis.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

Surface engineered porous silicon for stable, high performance electrochemical supercapacitors

Silicon materials remain unused for supercapacitors due to extreme reactivity of silicon with electrolytes. However, doped silicon materials boast a low mass density, excellent conductivity, a controllably etched nanoporous structure, and combined earth abundance and technological presence appealing to diverse energy storage frameworks. Here, we demonstrate a universal route to transform porous silicon (P-Si) into stable electrodes for electrochemical devices through growth of an ultra-thin, conformal graphene coating on the P-Si surface. This graphene coating simultaneously passivates surface charge traps and provides an ideal electrode-electrolyte electrochemical interface. This leads to 10–40X improvement in energy density, and a 2X wider electrochemical window compared to identically-structured unpassivated P-Si. This work demonstrates a technique generalizable to mesoporous and nanoporous materials that decouples the engineering of electrode structure and electrochemical surface stability to engineer performance in electrochemical environments. Specifically, we demonstrate P-Si as a promising new platform for grid-scale and integrated electrochemical energy storage.

36 MATERIALS SCIENCE↗