Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Distributed Generation”

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 919 records · Page 51

Introduction to Microgrids

A microgrid is simply a typical power system that has distributed energy resources (DERs), like small diesel generators, with loads in the network. The resiliency offered by microgrids to retain service during system events makes them a suitable support for the electric grid during natural disasters. One key factor that differentiates microgrids from a typical distribution grid is its ability to disconnect itself in the event of a blackout and sustain its loads through DERs. Microgrid systems will reconnect to the main grid when the grid comes back up. The resiliency offered by microgrids to retain service during system events makes them a suitable support for the electric grid during natural disasters. This presentation discusses the typical key steps involved in microgrid planning, installation, operation in a interoperable and interconnectable fashion.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Analysis of Electricity Price Differentials for the Utah Office of Energy Development

Operation Gigawatt is a vision for Utah to meet its growing energy needs with a focus on doubling energy generation. One of the greatest challenges to meet this goal will be the transmission and distribution system. Utah seeks analysis to inform discussions on new generation, the likely network congestion based on new power flows, and modeling of the state's energy infrastructure. To support this need, LBNL developed a baseline understanding of recent historical congestion in the existing transmission system based on electricity market prices. The resulting products visualize and provide metrics on congestion within Utah and between Utah and its neighbors during 2015-2023.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A Novel Technique for Distributed Measuring of Strain in the Vadose Zone

Characterizing water content and pressure changes in the vadose zone is important to understanding a variety of geologic processes, ranging from permeability, to evapotranspiration, and aquifer recharge. Changes in water content, or pressure, cause strain in the solid porous medium of the vadose zone, and it is possible that measuring those strains could be used a characterization tool. The Coherence-length-gated Microwave Photonics Interferometry (CMPI) technique measures strain at high resolution along many intervals defined by pairs of reflectors distributed along an optical fiber. This technique has recently been developed at Clemson University, and it has the spatial, and temporal resolution to characterize the strains that are expected to occur with hydrologic changes in the vadose zone. However, the technique has never been used to measure strain in porous media, so its capabilities remain uncertain. The objective of this thesis is to evaluate the ability of using CMPI to measure strain changes in the vadose zone. The research approach consists of conducting laboratory tests using a column filled with sand that was subjected to changes in water content and pressure. An optical fiber with CMPI reflectors was deployed in a high surface area ribbon and used to measure strain along the axis of the column. The column was made from 8-inch, Schedule 40 PVC pipe (75 cm tall, 20 cm inner diameter) and filled with medium-grained sand (K = 2.5x10-6 m/s, porosity = 0.28, Coefficient of Uniformity = 2, van Genuchten (n = 4.34 and  = 0.00039 m-1), Young’s Modulus = 18 – 47 MPa). Five pressure and four temperature sensors attached to probes inserted into the wall. The optical fiber includes five reflectors spaced 10 cm apart on the inside and outside of the column along the vertical axis. Strain measurements are made between pairs of reflectors. The reflectors are created in 250-micron-diameter acrylate-coated single mode Corning SMF28e+ optical fiber using a femtosecond laser. The optical fiber was laminated between two pieces of polyester film creating a large surface area to transfer strain from the porous media to the optical fiber. The experiments were conducted by filling the column from the bottom or infiltrating water from the top. The water was allowed to equilibrate to room temperature prior to each test in order to limit thermoelastic strain. Five injection tests with a rate of 250 ml/min and three infiltration tests at varying rates were conducted, and the results show patterns of strain and pressure are generally similar. Hydrologic conditions define three zones based on the pressure magnitude and distribution. 1.) Ambient Zone where the pressure heads are quasi-static and the pressure gradient is roughly unity (head gradient of zero). This is the uppermost zone and is characterized by negative pressures. 2.) Transition Zone where the pressures increase from ambient to zero, are changing relatively rapidly and the pressure gradient is relatively steep (pressure head gradients of 2). The Transition zone is 10 to 15 cm thick. 3.) Positive Pressure Zone where the pressure is positive, the rate of change is slower than in the transition zone and the gradient is flatter (pressure head gradient 1.1 to 1.2). This is the lowest zone in the column. Injection of water causes the pressure to increase and the three zones to move upward. The Transition zone moves at a rate of approximately 0.0002 m/s +/- .00005, according to analyses of pressure profiles. This is the same as the average velocity of the water calculated as volumetric flux/effective porosity =5.4x10-5 m/s/0.28. Strain signals in the range of 10s  were observed with a noise level of generally less than 0.1  (signal to noise ratio of greater than 100) during injection and drainage. The spatial and temporal distributions of strain caused by injection depend on the location of the moving pressure zones. The locations of the different pressure zones were determined from pressure profiles at different times and these data were transferred to strain time series. This showed that strain in the Ambient pressure zone is either unchanged or slightly compressive, whereas strain in Transition zone is tensile and roughly proportional to the pressure change, and strain ranges from tensile to compressive in the Positive Pressure zone. In general, the strain is variable at the top of the Positive Pressure zone, but it appears to be consistently tensile lower in the zone. The spatial distribution of strain causes three distinct stages in the strain times series measured between pairs of CMPI reflectors. The strain is unchanged or slightly compressive during Stage 1, it increases (positive strain is tensile) approximately linearly at a rate of 0.02 to 0.03 μ/s during the Stage 2, and then flattens out in Stage 3. The strain increases again during Stage 4 of the time series. The pressure also changes in stages. It is unchanged during Stage 1 and then increases relatively rapidly at a rate of between 4 and 10 Pa/s during Stage 2 and continues to increase during Stages 3 and 4, but at a rate slightly slower than during Stage 2 (from 2.5 to 3.5 Pa/s). The rates of pressure and strain change are consistent with basic analyses. The rate of pressurization is similar to the calculated velocity of flow in the saturated zone, and this is similar to the measured velocity of the strain increase. The ratio of the pressurization rate and strain rate during Stage 2 (2.5 (Pa/s) /0.03 (μ/s)) is 80 MPa, which is approximately equal to the upper range expected considering the uniaxial compression and the Young’s Modulus of the sand. An interesting effect occurs when the upper surface of the sand becomes saturated. Significant compression (several 10s of ) occurs as the pressure and saturation increase at the upper surface. Compression occurs throughout the column, but the effect is greatest at the top of the column. The rate of compression is fastest slightly before ponding occurs, but it slows markedly and nearly stops when water starts to accumulate at the surface (ponding). This effect reverses when the surface of the sand is drained, resulting in tension throughout the column. The compression caused by this effect can be as large or larger than the tensile strain that accumulated during filling. This effect was unexpected because increasing pressure is normally associated with tensile strain. Nevertheless, this effect was observed consistently in all tests where the pressure changes at the upper surface, including tests where water was injected from below or infiltrated from above. This effect behaves as if the pore pressure at the upper surface of the sand exerts a normal force on the boundary (increasing pore pressure exerts a downward compression on the boundary). The results of these laboratory experiments indicate that CMPI can measure strain caused by fluid pressure changes in the vadose zone with a signal to noise ratio of 100 or more. Injection and drainage cause a strain signal that is complex, but repeatable. The magnitudes of the strain signal are consistent with magnitudes that are expected based on poroelastic calculations using independently measured properties of the sand. These results indicate that the CMPI technique with an optical fiber laminated in a polyester ribbon generates data that represent the strain distribution during hydrologic processes of imbibition and drainage in variably saturated sand. This suggests that distributed strain measurements using CMPI could be a viable approach for evaluating processes in the vadose zone, laying the groundwork for future field implementation.

47 OTHER INSTRUMENTATION↗

Development of Supervisory Control System for Thermal Energy Distribution System

The integrated energy system (IES) refers to the combination of nuclear energy generation with other energy sources to enable the efficiency and reliability of power generations. To create technologically viable and economically competitive systems, supervisory control strategies are critical for optimizing performance and ensuring stability across different energy generation, transportation, and utilization. This work focuses on the control strategies for the thermal energy distribution system (TEDS), which is a cornerstone of the Dynamic Energy Transport and Integration Laboratory at Idaho National Laboratory. TEDS currently relies on operators to coordinate across different components to manage energy storage and ensure efficiency. This work demonstrates the use of model predictive control (MPC) with surrogate models in determining optimal setpoints for major TEDS components. The capability of MPC-based supervisory control system is evaluated by autonomously matching the power outputs from a Dymola-based TEDS with target heat demands.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Record magnetic field generation by short-pulse laser-driven capacitor-coil targets

Magnetic fields generated by capacitor-coil targets driven by intense short-pulse lasers have been characterized using ultrafast proton radiography. A 1-kJ, 15-ps laser at a center wavelength of 1053 nm irradiated the back plate of the capacitor with an intensity of ∼8.3 × 10 18 W/cm 2 , creating ultra large currents in the connecting coils. High-quality proton data obtained in the axial probing geometry show definitive signatures of magnetic field generation, allowing precision measurement of the field distribution and strength. The data show a coil current of 120 ± 10 kA producing 200 ± 20 Tesla magnetic fields at the coil center at 1.127 ns after the laser drive. In conclusion, this sets a record for magnetic field generation by the short-pulse-powered capacitor-coil targets.

Gao, Lan [Princeton Plasma Physics Laboratory (PPP↗

Multi-Probe 23A Pre-Shot report

The Multi-Probe Radiography campaign will field a full day of experiments at the Omega EP laser facility. The experiments will consist of shots that alternate between Omega EP’s two short-pulse laser beams to generate proton and deuteron beams from a variety of film and foam targets. The proton beam is intended to be used as a radiographic source, while the deuteron beam will be used as the pitcher in a pitcher-catcher neutron-generation concept. For the proton beam, our goal is to characterize the angular distribution and energy spectra of unwanted x-rays that are generated during proton acceleration and were found to be problematic for simultaneous radiography in MP-22A. As part of this primary goal, we will test whether an x-ray shield, XBLK, can be used to prevent these x-rays from impacting other x-ray diagnostics. For the deuteron beam, our primary goal is to characterize the energy spectra and beam profile of deuterons accelerated from CD foils and foams.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A deep generative model enables automated structure elucidation of novel psychoactive substances

Over the past decade, the illicit drug market has been reshaped by the proliferation of clandestinely produced designer drugs. These agents, referred to as new psychoactive substances (NPSs), are designed to mimic the physiological actions of better-known drugs of abuse while skirting drug control laws. The public health burden of NPS abuse obliges toxicological, police and customs laboratories to screen for them in law enforcement seizures and biological samples. However, the identification of emerging NPSs is challenging due to the chemical diversity of these substances and the fleeting nature of their appearance on the illicit market. Here, in this study, we present DarkNPS, a deep learning-enabled approach to automatically elucidate the structures of unidentified designer drugs using only mass spectrometric data. Our method employs a deep generative model to learn a statistical probability distribution over unobserved structures, which we term the structural prior. We show that the structural prior allows DarkNPS to elucidate the exact chemical structure of an unidentified NPS with an accuracy of 51% and a top-10 accuracy of 86%. Our generative approach has the potential to enable de novo structure elucidation for other types of small molecules that are routinely analysed by mass spectrometry.

Cheminformatics↗

Combining Experimental and Theoretical Techniques to Gain an Atomic Level Understanding of the Defect Binding Mechanism in Hard Carbon Anodes for Sodium Ion Batteries

Sodium ion batteries (NIBs) are an attractive alternative to lithium–ion batteries in applications that require large–scale energy storage due to sodium's high natural abundance and low cost. Hard carbon (HC) is the most promising anode material for NIBs; however, there is a knowledge gap in the understanding of the sodium binding mechanism that prevents a rational design of HC. This study tunes sucrose–derived HC via synthesis temperature then evaluates the structural, physical, and electrochemical properties. Neutron total scattering is used to generate structural models by fitting pair distribution functions (PDF) with a combination of molecular dynamics and reverse Monte Carlo methods. From this model, the number and type of structural features are identified, quantified, and correlated to the galvanostatic charge/discharge. A method of PDF “fingerprinting” binding sites using Na probe atoms is developed and analyzing these PDFs reveals an atomistic view of ion binding sites responsible for “defect” storage mechanisms. Combining these techniques results in an atomic–level study that provides a big picture of the Na–binding mechanism in NIBs, which allows for more precise tuning of the structure–property relationships in the future. Finally, the methodologies developed will also enable new strategies for the analysis of amorphous functional materials.

25 ENERGY STORAGE↗

The Kimberlina synthetic multiphysics dataset for CO 2 monitoring investigations

Abstract We present a synthetic multi‐scale, multi‐physics dataset constructed from the Kimberlina 1.2 CO 2 reservoir model based on a potential CO 2 storage site in the Southern San Joaquin Basin of California. Among 300 models, one selected reservoir‐simulation scenario produces hydrologic‐state models at the onset and after 20 years of CO 2 injection. Subsequently, these models were transformed into geophysical properties, including P‐ and S‐wave seismic velocities, saturated density where the saturating fluid can be a combination of brine and supercritical CO 2 , and electrical resistivity using established empirical petrophysical relationships. From these 3D distributions of geophysical properties, we have generated synthetic time‐lapse seismic, gravity and electromagnetic responses with acquisition geometries that mimic realistic monitoring surveys and are achievable in actual field situations. We have also created a series of synthetic well logs of CO 2 saturation, acoustic velocity, density and induction resistivity in the injection well and three monitoring wells. These were constructed by combining the low‐frequency trend of the geophysical models with the high‐frequency variations of actual well logs collected at the potential storage site. In addition, to better calibrate our datasets, measurements of permeability and pore connectivity have been made on cores of Vedder Sandstone, which forms the primary reservoir unit. These measurements provide the range of scales in the otherwise synthetic dataset to be as close to a real‐world situation as possible. This dataset consisting of the reservoir models, geophysical models, simulated time‐lapse geophysical responses and well logs forms a multi‐scale, multi‐physics testbed for designing and testing geophysical CO 2 monitoring systems as well as for imaging and characterization algorithms. The suite of numerical models and data have been made publicly available for downloading on the National Energy Technology Laboratory's (NETL) Energy Data Exchange (EDX) website.

58 GEOSCIENCES↗

Visualization of Hemispherical Post‐Detonation Fireball Internal Structures

Hemispherical charges are initiated above a transparent plate in a half‐plane configuration, allowing optical access to the internal regions of the luminous fireball. High‐speed visualization enables characterization of the internal luminous structure, showing clearly the bright shell and dark core regions. Spectroscopic and pyrometric diagnostics are applied to these flows, generating quantitative information on the spatial distribution of temperature inside the fireball. Both cased and uncased charges can be examined using this methodology, and several different high explosives are tested. This approach can be useful in validation of detailed explosive fireball models. Initial comparisons with such models are presented.

detonation↗

Error-mitigated data-driven circuit learning on noisy quantum hardware

Application-level benchmarks measure how well a quantum device performs meaningful calculations. In the case of parameterized circuit training, the computational task is the preparation of a target quantum state via optimization over a loss landscape. This is complicated by various sources of noise, fixed hardware connectivity, and generative modeling, the choice of target distribution. Gradient-based training has become a useful benchmarking task for noisy intermediate-scale quantum computers because of the additional requirement that the optimization step uses the quantum device to estimate the loss function gradient. In this work, we use gradient-based data-driven circuit learning to qualitatively evaluate the performance of several superconducting platform devices and present results that show how error mitigation can improve the training of quantum circuit Born machines with 28 tunable parameters.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Reactor and fuel cycle performance of light water reactor fuel with 235 U enrichments above 5%

Recent advances in nuclear fuel materials research, particularly on the topic of accident-tolerant fuels, have brought up potential opportunities for expanding the operating envelope of existing light water reactors. As many of the performance improvements offered by these technologies may be most fully realized by increasing fuel enrichment beyond the standard 5% limit, this paper examines the potential reactor performance and fuel cycle performance of low-enriched uranium oxide fueled light water reactors by generically considering pressurized water reactors with 235 U enrichment from 5 to 7%. Advanced cladding, including accident-tolerant cladding, has the potential to increase fuel burnup limits related to hydrogen in the cladding that coincide with those limits associated with end-of-life reactivity. Therefore, higher enrichment will be necessary in order to realize the higher fuel burnups. This work includes evaluation of the fuel cycle length, discharge burnup, reactivity coefficients, and fuel cycle performance, including radioactive waste and environmental impact metrics per unit energy generated. The analysis was performed using the evaluation metrics from the US Department of Energy Office of Nuclear Energy Fuel Cycle Evaluation and Screening Study. The reactor performance and safety analysis show that enrichments between 5 and 7% would have similar fuel temperature and moderator temperature coefficients. However, the soluble boron coefficient would decrease in magnitude, requiring more corrosive boric acid in the coolant or other methods of reactivity control during the fuel cycle. At these higher enrichments the maximum burnup at the rim of the fuel pellet would increase by almost a factor of two, which is expected to impact the formation of high-burnup structure in the fuel and the corresponding thermo-mechanical fuel properties. The fuel cycle performance assessment shows that increasing enrichment reduces the quantity of high-level waste disposed per unit energy generated, but it increases the natural resource requirements normalized to a gigawatt-electricity-per-year basis. Another impact is the slightly higher discharge burnup, resulting in somewhat different activity levels of the spent nuclear fuel and high-level waste radioactivity at 100 and 100,000 years after fuel discharge. The environmental impacts—including land use, water use, carbon emission, and radiological exposure—are of the same magnitude per unit energy generated. Finally, however, the impacts are distributed differently. Less than 5% enrichment has marginally more impact on the back-end of the fuel cycle, and greater than 5% enrichment has marginally more impact on the front-end of the fuel cycle. Ultimately, no neutronic or reactor safety hindrances to employing light water reactor fuel with enrichments greater than 5% are identified; given the achievable reactor performance benefits with advanced fuels, further practical exploration of increased enrichment fuel is recommended.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

In situ multi-tier auto-ignition detection applied to dual-fuel combustion simulations

Here we use an anomaly detection methodology that is centered on analyzing fourth-order joint moments (co-kurtosis), particularly focusing on its application in auto-ignition of combustion problems with large numbers of species. Unsupervised anomaly detection is challenging to generalize across problem types and domains. A recent technique, centered on analyzing information in the fourth-order joint moment co-kurtosis, has shown promise, especially for high-dimensional scientific data. In this work we present developments to the co-kurtosis based anomaly detection method needed to make it effective and scalable for large-scale distributed scientific data, such as those generated by massively parallel simulations. An in situ co-kurtosis algorithm is employed as the anomaly detection method for identifying ignition kernels in simulations of turbulent combustion. Here, we extend an existing methodology which identifies regions of the domain where anomalies are present, and add another tier of anomaly detection where the individual samples contributing to the anomaly are identified. We apply this algorithm on-the-fly to a variety of turbulent reacting flow problems and compare it to the widely used (but significantly more expensive) chemical explosive mode analysis (CEMA). We demonstrate the ability of the method to detect and identify the onset of low and high temperature ignition which can be used for computational steering, as chemical and combustion anomalies occur intermittently at spatio-temporal locations unknown a priori. Finally, we apply our lightweight in situ algorithm to an exascale high-fidelity simulation with a total of 2.4 Trillion degrees of freedom, performed using an adaptive mesh refinement solver. Furthermore, through a scalability analysis, we show that the relative computational cost of this in-situ anomaly detection algorithm compared to an iteration of the reacting flow solver is negligible.

97 MATHEMATICS AND COMPUTING↗

Large eddy simulation of wind farm performance in horizontally and vertically staggered layouts

This numerical investigation employs Large Eddy Simulation (LES) coupled with Actuator Disk Model (ADM) to evaluate wind farm layout optimization strategies. The study presents a systematic analysis of aligned, horizontal staggering, vertical staggering, and mixed (combination of horizontal and vertical) staggering configurations, aiming to establish optimal design parameters for enhanced power production. The investigation examines key performance metrics including mean velocity distributions, turbulence intensity characteristics, and power generation efficiency. Results demonstrate better performance of both horizontal and vertical staggering patterns compared to conventional aligned configurations, with horizontal staggering exhibiting notably higher power output than vertical arrangements. Our findings also suggest that mixed configurations, incorporating both horizontal and vertical staggering, can offer optimal performance characteristics. As a result, this research advances the understanding of wake interactions in complex wind farm layouts and provides design guidelines for maximizing wind farm power generation efficiency through strategic turbine positioning.

17 WIND ENERGY↗

In situ characterization of phase transformations in petroleum pitch by high temperature X-ray diffraction

The phase transformation from isotropic to mesophase of a petroleum-derived pitch was monitored by in situ X-ray diffraction at 410 °C. The kinetics of the transformation were characterized by monitoring the growth of a broad peak at 2θ~24° and could be described by Avrami’s equation. The volatilization of low molecular weight components originally present in the pitch, or produced during the condensation reaction, was assessed by thermal gravimetry analysis, following the same thermal history. Ex situ characterizations of the pitch precursor before and after thermal treatment provided insight into the evolution of the molecular structure associated with mesophase generation (optical polarized microscopy, molecular weight distribution by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry, X-ray photoelectron spectroscopy, Auger spectroscopy, and thermal stability). Here, this investigation demonstrates the value of using complementary in situ and ex situ experimental techniques to characterize complex mechanisms during the processing of organic materials.

02 PETROLEUM↗

Healing of nanocracks by collision cascades in nickel

We use molecular dynamics simulations to investigate the interaction between collision cascades and pre-existing nano-cracks in nickel (Ni). We find that collision cascades cause the cracks to heal whenever a cascade initiates a thermal spike whose core overlaps with the crack. Moreover, crack healing affects the distribution of radiation-induced defects generated during a collision cascade: defect clusters, such as dislocations and stacking fault tetrahedra, dominate in the presence of crack healing, while isolated point defects dominate in its absence. Furthermore, this work indicates that irradiation may reduce the number of pre-existing crack-like flaws in metallic components, potentially improving their mechanical resilience.

36 MATERIALS SCIENCE↗

TransformerG2G: Adaptive time-stepping for learning temporal graph embeddings using transformers

Dynamic graph embedding has emerged as a very effective technique for addressing diverse temporal graph analytic tasks (i.e., link prediction, node classification, recommender systems, anomaly detection, and graph generation) in various applications. Such temporal graphs exhibit heterogeneous transient dynamics, varying time intervals, and highly evolving node features throughout their evolution. Hence, incorporating long-range dependencies from the historical graph context plays a crucial role in accurately learning their temporal dynamics. In this paper, we develop a graph embedding model with uncertainty quantification, TransformerG2G, by exploiting the advanced transformer encoder to first learn intermediate node representations from its current state (t) and previous context (over timestamps [t–1,t–l], l is the length of context). Moreover, we employ two projection layers to generate lower-dimensional multivariate Gaussian distributions as each node's latent embedding at timestamp t. We consider diverse benchmarks with varying levels of "novelty" as measured by the TEA (Temporal Edge Appearance) plots. Here, our experiments demonstrate that the proposed TransformerG2G model outperforms conventional multi-step methods and our prior work (DynG2G) in terms of both link prediction accuracy and computational efficiency, especially for high degree of novelty. Furthermore, the learned time-dependent attention weights across multiple graph snapshots reveal the development of an automatic adaptive time stepping enabled by the transformer. Importantly, by examining the attention weights, we can uncover temporal dependencies, identify influential elements, and gain insights into the complex interactions within the graph structure. For example, we identified a strong correlation between attention weights and node degree at the various stages of the graph topology evolution.

97 MATHEMATICS AND COMPUTING↗

The Donnan potential revealed

Selective transport of solutes across a membrane is critical for many biological, water treatment and energy conversion and storage systems. When a charged membrane is equilibrated with an electrolyte, an unequal distribution of ions arises between phases, generating the so-called Donnan electrical potential at the solution/membrane interface. The Donnan potential results in the partial exclusion of co-ion, providing the basis of permselectivity. Although there are well-established ways to indirectly estimate the Donnan potential, it has been widely reported that it cannot be measured directly. Here we report the first direct measurement of the Donnan potential of an ion exchange membrane equilibrated with salt solutions. Our results highlight the dependence of the Donnan potential on external salt concentration and counter-ion valence, and show a reasonable agreement with current theoretical models of IEMs, which incorporate ion activity coefficients. By directly measuring the Donnan potential, we eliminate ambiguities that arise from limitations inherent in current models.

25 ENERGY STORAGE↗