Search NASA⌕ Search

SEARCH · Search NASA

Results for “FLUID DYNAMICS”

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 505 records · Page 28

Applying a Compact Porous Media Model to Numerically Derive Resistance Coefficients for Lattice Structures

Additive Manufacturing allows for exploring various geometries to achieve specific engineering criteria. Lattices are one geometry with unique properties, including being periodically repeating structures which allow flow through them to be represented as a porous media according to Darcy-Forchheimer equations. These equation’s coefficients are generally experimentally derived, but this work demonstrates the ability to numerically derive them with CFD. Simulations were performed using three-dimensional stead state Reynolds-averaged Navier-Stokes with a k-ω Shear Stress Transport turbulence model using Ansys Fluent. Three lattice geometries were investigated and drag coefficients were derived. The method was validated against externally published data for similar geometries demonstrating strong agreement, and grid convergence for all simulations was calculated with a Grid Convergence Index method. Wall roughness is demonstrated to have a non-negligible impact on results and roughness values are considered for the primary focus Octahedral geometry where both smooth wall and rough wall coefficients were derived. The porosity coefficients for the Octahedral geometry at 1.0 [m/s] were found to be 2.89×10 6 and 2.90×10 6 [1/(Pa*m*s)] for the permeability coefficients, 6.37×10 1 and 5.44×10 1 [m 2 /kg] for the inertial resistance coefficients, and with a max pressure drop of 5116.7 [Pa] and 4429.5 [Pa] for the smooth walls and rough walls, respectively. The derived numerical method enables rapid exploration and optimization of new lattice designs for diverse engineering applications.

42 ENGINEERING↗

Effect of soot concentration on radiative transport in a 2 m JP-8 pool fire

Hubbard et al. (2024) investigated conjugate heat transfer in a 7.9 m JP-8 pool fire, utilizing both heavily sooting and non-sooting fuel models. Their findings indicated that all models produced pool radiative heat fluxes within the measurement uncertainty for large-scale outdoor experiments. This study aims to further investigate the influence of soot concentration on radiative transport. To this end, Large Eddy Simulations (LES) of a 2 m diameter JP-8 pool fire were conducted, employing four turbulent combustion chemistry models: the Eddy Dissipation Concept (EDC) and three variants of flamelet model. Additionally, three variations of the one-equation subgrid scale turbulence kinetic energy model, and multiple mesh resolutions were utilized. The primary validation metrics were pool radiative heat flux and plume radiative fraction, with simulation data compared to experimental results published by Blanchat et al. (2010). The EDC model produced high soot concentrations near the pool surface, resulting in a significant shielding effect. In contrast, the flamelet combustion model has different kinetic limiters for soot production, leading to soot formation higher in the flame, where it can be more readily oxidized. Flamelet models incorporating both acetylene and benzene contributions to soot growth demonstrated better agreement with experimental data. To assess the sensitivity of pool radiative heat flux to soot mass concentration, we used a one-dimensional approximation over a wide range of soot concentrations and compared it to three-dimensional solutions obtained using Discrete Ordinates. Furthermore, this analysis illustrates the impacts of soot shielding and Turbulence Radiation Interactions (TRI), with TRI emission enhancement factors predicted to be in the range of 3–5 near the edges of the plume.

Computational fluid dynamics↗

Hybrid combustion modeling approach for turbulent jet ignition in natural-gas pre-chamber spark-ignition engines at high EGR

Here, this study presented a hybrid modeling approach for simulating turbulent jet ignition and combustion processes in a natural-gas pre-chamber spark-ignition engine operating under exhaust gas recirculation (EGR) diluted conditions. In-depth analyses of experimental data and simulation results from previous work [Chinnathambi et al., ICEF2021-67836; Kim et al., Fuel 409: 137815, 2026] revealed two key findings: (i) the magnitude of pressure difference between the pre-chamber and main chamber ($∆P_{PC-MC}$) was positively correlated with the combustion duration from the moment of $∆P_{PC-MC}=0$ to the point of 5% mass fraction burned, with larger $∆P_{PC-MC}$ associated with longer duration; and (ii) the turbulent combustion regime in the main chamber transitioned from the broken reaction zone to the corrugated flamelet regime, with the Karlovitz number exceeding 100 immediately after turbulent hot jets were ejected from nozzles, coinciding with observed local extinction events. To accurately simulate the entire combustion process, a hybrid approach was developed under Reynolds-Averaged Navier Stokes framework, combining the G-equation model for pre-chamber combustion with the multi-zone well-stirred reactor approach and a turbulence-chemistry interaction (TCI) submodel for main chamber combustion. The TCI submodel accounted for the attenuation of reaction rates due to turbulent strain and modeled local extinction by suppressing reaction rates under certain flow and flame conditions. When applied to three EGR rate conditions toward the dilution limit, the hybrid modeling approach accurately reproduced experimental data in terms of cylinder pressure, apparent heat release rate, and the observed positive correlation, including the delayed onset of main chamber combustion—a feature not captured by existing combustion models.

computational fluid dynamics simulation↗

Modeling and Experimental Validation of a Direct-Contact Counter-Flow Fluidized Bed Heat Exchanger for Thermal Energy Storage (TES) Applications

Particle-based thermal energy storage (TES) systems are an emerging energy storage technology. The technological advances have reduced costs, making TES more competitive and reliable in the marketplace but an efficient and reliable operation is heavily dependent on coherent heat transfer between air to particles or vice versa. The particle-based TES technologies provide an intermediate system that can store energy for short (0-10 h), long (10-200 h) and seasonal (> 200 h) timescales. The TES systems store energy by converting electricity to thermal energy; electricity can be directly sourced intermittent generation technologies and/or the grid, helping manage peak loads and other mismatches in supply and demand. The overall efficiency of the TES system depends on the performance of system components (particle storage silos and particle transfer mechanism etc.). The particle heat exchanger is one of the key system components that affects the system efficiency. The pressurized fluidized bed heat exchanger (PFB HX) performance is challenging to predict due to the chaotic behavior of particle and fluid interaction. This research presents a computational study of a novel direct-contact, counter-flow and air-to-particles PFB HX, that contributes in advancing the particle-based long-duration TES technologies. For the current analysis an unsteady Eulerian-Eulerian CFD model was developed and validated against experiments performed at the National Laboratory of the Rockies for two particle sizes (600 ..mu..m and 825 ..mu..m ). Following validation, parametric simulations were conducted to evaluate the effects of interphase drag models (Syamlal-O'Brien and Gidaspow), particle size, bed height and the influence of a frictional-viscosity term on hydrodynamics and heat transfer between the air & particles. The key findings from the analysis are: (1) for the studied operating window Syamlal-O'Brien provides superior agreement with measured gas temperatures (errors generally < 10%) while Gidaspow shows large deviations for the coarse particle case; (2) model predictions are most sensitive in the lower 0.2 m above the air distributor where bubble initiation and local mixing dominate interphase heat transfer; (3) representation of the distributor (number of inlet ports) materially affects predicted local mixing and temperature stratification; and (4) the Eulerian-Eulerian framework reproduces bulk thermal trends but shows regime dependent limitations for coarse particles, motivating mesoscale informed closures for scale-up analysis for future studies. These results provide validated guidance for drag selection and distributor design in particle-based thermal energy storage applications. Collectively, the validated model and parametric results quantify key drivers of PHB-HX performance and provide practical guidance for design and optimization. The results provide confidence in the model predictability and provide a step forward to improve on heat exchange performance. The demonstrated performance and modeling approach support the deployment and further development of this novel PHB-HX concept for robust, particle-based long-duration thermal energy storage systems.

25 ENERGY STORAGE↗

Derivation of physical equations for high-speed laser welding using large language models

It is challenging to formulate complex physical phenomena that occur in a manufacturing process, particularly when the available data are limited, rendering conventional data-driven approaches ineffective. This study aims to predict humping onset in high-speed laser welding by introducing a novel framework, namely text-to-equations generative pre-trained transformer (T2EGPT). This method leverages the capabilities of large language models (LLMs), in combination with sparse experimental data and enriched literature data, to derive an interpretable and generalizable equation for predicting humping initiation. By capturing key correlations among physical parameters, T2EGPT generates a compact and dimensionless expression that accurately predicts hump formation. The equation reveals that humping arises from the interplay between inertia-driven backward melt flow and capillary-driven surface stabilization, where inertial forces drive molten metal backward and capillary forces resist surface deformation. Furthermore, compared to traditional data-driven models, T2EGPT demonstrates enhanced predictive accuracy and cross-material transferability. More broadly, this study highlights the potential of LLMs to integrate textual information with data-driven discovery, enabling the extraction of physical laws in data-scarce scientific domains.

36 MATERIALS SCIENCE↗

A weighted shifted boundary method for immersed moving boundary simulations of Stokes' flow

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods, and relies on reformulating the original boundary value problem over a surrogate (approximate) computational domain. The surrogate domain is constructed so as to avoid cut cells and the associated problematic implementation and numerical integration issues. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions: hence the name of the method, that shifts the location and values of the boundary conditions. Here, in this article, we extend the SBM to the simulation of incompressible Stokes flow, by appropriately weighting its variational form with the elemental volume fraction of active fluid. This approach allows to drastically reduce spurious pressure oscillations in time, which are produced if the total volume of active fluid were to change abruptly over a time step. The proposed Weighted SBM (W-SBM) exactly preserves states of hydrostatic equilibrium, and induces small mass and momentum conservation errors, which converge as the grid is refined. This is in analogy to cutFEMs and related unfitted approaches, which rely on an affine representation of cut boundaries. We demonstrate the robustness and accuracy of the proposed method with an extensive suite of two-dimensional tests.

97 MATHEMATICS AND COMPUTING↗

Cost-efficient finite-volume high-order schemes for compressible magnetohydrodynamics

We present an efficient dimension-by-dimension finite-volume method which solves the adiabatic magnetohydrodynamics equations at high discretization order, using the constrained-transport approach on Cartesian grids. Results are presented up to tenth order of accuracy. The algorithmic architecture of this method is very close to that of commonly employed second-order schemes: it requires only one reconstructed value per face for each computational cell, independently of the scheme's order. This property is highly beneficial for the numerical efficiency. It results from reusing the required values already available in neighboring grid cells, in contrast to standard algorithms that require a number of reconstructions and evaluations which increases with the scheme's order of accuracy. At a given resolution, these high-order schemes present significantly less numerical dissipation than commonly employed lower-order approaches. Thus, results of comparable accuracy are achievable at a substantially coarser resolution, yielding overall performance gains. We also present a way to include physical dissipative terms: viscosity, magnetic diffusivity and cooling functions, respecting the finite-volume and constrained-transport frameworks. Benefits of this method are shown through applications in turbulent flows.

97 MATHEMATICS AND COMPUTING↗

Numerical simulation of vortex-induced vibration response of a single IEA 10-MW wind turbine blade

Three-dimensional simulation of vortex-induced vibration (VIV) of a single International Energy Agency (IEA) 10-MW reference wind turbine blade with a length of 97.325 m is performed using the ExaWind stack, an open-source suite of codes. This study aims to illustrate the spanwise VIV response characteristics and cross-validate the results with an existing commercial framework. Five near-body meshes and three time steps are selected for the convergence study. To improve computational efficiency, several VIV triggering methods are also compared to shorten the VIV development period. The ExaWind-based VIV simulation strategy for a single IEA 10-MW blade is determined. First, the modal shape is validated against published results. Then, spanwise VIV responses of four blade configurations under a fixed and varied incoming flow velocity are analyzed. Results show that the VIV response is dominated by the first edgewise (second overall) mode. Little first-mode contributions appear near the second-mode node, producing a pi phase jump, and a higher harmonics response occurs near the blade root. Rotational degrees of freedom are minor compared with translational motion. The response versus reduced velocity is analyzed, showing a two-branch behavior similar to that of VIV for a bluff cylinder. Across all tested cases, the dominant frequency remains locked to the natural frequency of the second mode with no observed desynchronization. A mild deviation is observed for the case of 90-degree pitch and 310-degree azimuth rotation near a reduced velocity of 6, which will be examined with additional cases in future work. These findings indicate that severe VIV responses can arise under specific configurations and flow conditions, thereby increasing the potential for VIV fatigue damage and requiring greater attention during operation.

17 WIND ENERGY↗

Defect-Mediated Diffusion Pathways in Spodumene Accelerate Lithium Transport

Lithium extraction from naturally occurring α-spodumene is hindered by poor lithium diffusivity, necessitating high-temperature phase transformation to a low-density β polymorph. Although β spodumene exhibits up to 5 orders of magnitude higher lithium-ion diffusivity, both phases have diffusion activation energies between 0.8 and 1 eV, indicating that polymorph density is not the controlling factor over diffusivity. We show that aluminum vacancies facilitate lithium-ion diffusion in α-spodumene by reducing the migration barrier from 2.4 to 0.9 eV. Bond valence site energy and nudged elastic band calculations show a new lithium local minimum site which promotes a one-dimensional percolation network by reducing the lithium intersite distance from 4.5 Å to 2.9 Å. However, aluminum vacancies are energetically unfavorable to percolate through the whole structure, resulting in very low net lithium diffusivity and highlighting the critical role of nonstoichiometric defects in facilitating lithium transport in rigid aluminosilicate structures.

Chemical structure↗

Enhanced Stability and Sensitivity for CA-125 Detection Under Microfluidic Shear Flow Using Polyethylene Glycol-Coated Biosensor

The microfluidic-based point-of-care (POC) diagnostic tool has garnered significant interest in recent years, offering rapid and cost-effective disease detection. There is a growing trend toward integrating microfluidic platforms with biosensors, aligning lab-on-a-chip technologies with POC diagnostic devices. Despite numerous efforts to incorporate biosensors into microfluidic systems, researchers have performed very limited investigations on the stability of biomarker detection when biosensors operate under microfluidic shear flow conditions. Gold nanoparticles (AuNPs) are a widely employed material in capacitive biosensors for antibody immobilization and sensitivity enhancement. However, AuNPs have limitations in providing stable detection of biomarkers within microfluidic shear flow due to their agglomeration nature. This study addresses these limitations by employing 2 kDa polyethylene glycol (PEG) as an intermediate biofunctional layer to immobilize CA-125 antibodies on gold-interdigitated electrodes for the stable and accurate detection of CA-125 antigens. The stabilities and sensitivities of AuNPs and PEG-coated biosensors are evaluated under both static drop and microfluidic shear flow conditions for CA-125 antigen detection. The experimental results demonstrate a capacitive signal response (5660 pF at 10 kHz) 2.2 times higher using the PEG-coated biosensor than the signal (2551 pF at 10 kHz) measured by the AuNP-coated biosensor in the detection of CA-125 antigen–antibody conjugation under static drop conditions, indicating the higher sensitivity of the PEG-coated biosensor. Additionally, the PEG-coated biosensor exhibits better consistency for the CA-125 antigen detection between static drop and microfluidic shear flow conditions (Cp decrease in percentage (ΔCp%↓) = 2.9% at 10 kHz) compared to the electrical signals measured using the AuNP-coated biosensor (ΔCp%↓ = 32.4% at 10 kHz), which suggests that the PEG-coated biosensor demonstrates higher stability for CA-125 antigen detection under microfluidic shear flow conditions. With these significant improvements brought by the PEG-coated biosensor, especially under microfluidic conditions, a substantial hurdle in developing electrical biosensors for POC diagnostic applications has been overcome, expediting further advancements in the field.

36 MATERIALS SCIENCE↗

Large Eddy Simulations of the Atmospheric Boundary Layer Over Satellite‐Sensed Sea Ice Maps

Surface heterogeneity in the marginal ice zone (MIZ) causes multiscale secondary atmospheric circulations that are challenging to model or observe. The absence or inadequate representation of these circulations in ocean‐atmosphere exchange schemes in climate models is partially responsible for the underestimation of Arctic sea ice loss. Observationally, such circulations obfuscate the interpretation of polar atmospheric chemistry measurements, among others. To address this open challenge, large‐eddy simulations are conducted over real‐world satellite‐sensed sea ice maps with an ice‐sea temperature contrast, as well as over idealizations of these maps that alter the ice pattern but conserve its fraction, showing that the ice fraction of a surface is not sufficient to predict the surface heat flux. In a second suite of simulations, three other heterogeneity metrics (representing the surface fragmentation, patch edge tortuosity, and patch size variability) are introduced to complement the ice fraction in describing the surface. Simulations varying these parameters suggest that they also significantly impact surface‐air interactions. A roughness contrast is then added to the surface temperature contrast, indicating that the contribution of roughness changes to the resulting atmospheric circulations is less pronounced than that of thermal heterogeneity. Based on these results, we illustrate, using a multi‐linear regression on these map features, that MIZ surface parameterizations in Earth Systems Models can be improved if they account for these various characteristics of the sea‐ice patterns.

54 ENVIRONMENTAL SCIENCES↗

Dust Direct Radiative Effect Including Large Particles and Component Minerals

The direct radiative effect (DRE) of dust aerosols in Earth system models (ESMs) remains highly uncertain, largely due to inadequate representations of particle size distribution (PSD) and mineral composition. Using NASA's Earth Surface Mineral Dust Source Investigation (EMIT) soil mineralogy data and observed PSD in an ESM that resolves dust mineral composition and emitted diameters from 0.1 to 70 μm, we find a near-neutral global dust net DRE (−0.057 W m −2 ), weaker than most previous estimates. Large dust (diameter >10 μm) contributes 30% of the global longwave dust optical depth, providing observational constraints on large-particle abundance, and offsets 20% of the dust shortwave cooling over major source regions. Incorporation of EMIT mineralogy reduces shortwave uncertainty by more than 50%. The remaining uncertainty mainly exists in processes controlling dust abundance, particularly the poorly understood transport of large dust, and longwave optical properties, which require additional observational constraints to more accurately quantify the dust DRE.

Li, Longlei [Cornell Univ., Ithaca, NY (United Sta↗

Role of Diurnal Cycle of Insolation on the MJO Propagation in the Maritime Continent

The diurnal cycle of convection in the Maritime Continent (MC) has been hypothesized to act as a barrier to the eastward propagation of the Madden‐Julian oscillation (MJO). To test this hypothesis, we use a regional model with realistic MJO to simulate an event from the boreal spring of 2013 that weakened and stalled over the MC. Two simulations are conducted: one that includes the diurnal cycle of insolation (CTL), and another without it (NO_DC). The MJO in the simulations was identified and tracked using a large‐scale precipitation tracking method that distinguishes propagation and non‐propagation unlike the usual Real‐time Multivariate MJO method. In the NO_DC simulation, the absence of diurnal heating reduces land precipitation, allowing more continuous eastward MJO propagation. An analysis of moist static energy budget reveals that MJO maintenance in NO_DC is due to increased longwave heating and reduced advection, whereas the persistent MJO propagation in NO_DC is due to increased advection and reduced longwave heating and surface latent heat flux. These processes, however, may vary across different parts of the MC, emphasizing the complexity of MJO propagation across the MC.

Zhou, Xin [National Center for Atmospheric Researc↗

Modeling Large Dust Aerosols in the Community Earth System Model Version 2 (CESM2)

Dust aerosols have a wide size distribution from less than 0.1 to over 100 μm and dominate Earth's atmospheric aerosol mass. However, most Earth system models (ESMs) inadequately represent dust aerosols larger than 10 μm in diameter, limiting the accuracy of the simulated dust cycle and climate impacts. Here, we introduce a new modeling framework that captures the full observed size distribution of dust aerosols, incorporating recent advances into a mineral-resolved version of the Community ESM, while addressing known issues in previous versions. Comprehensive evaluation against diverse observations of bulk dust and component minerals demonstrates that the model reproduces the observed dust cycle across particle sizes. Incorporating the previously unrepresented large-dust fractions substantially alters dust budget estimates, highlighting potential changes in simulated climate impacts and underscoring the importance of comprehensive size-resolved dust modeling. Despite these advancements, uncertainties persist. Our results indicate that a size-dependent reduction in settling velocity is required to reproduce the observed dust size distribution downwind of source regions. Specifically, in the new model, the gravitational settling velocity of dust particles larger than 10 μm in diameter must be reduced by as much as 85% to achieve agreement with observations. This empirical reduction serves as a constraint on physics-based models of dust settling. Future developments should address misrepresented physical processes that hinder accurate modeling of the large dust aerosol transport. Expanding observational data sets covering the full-size distribution is also essential to better constrain the dust cycle and improve the representation of dust optical properties and climate effects.

Li, Longlei [Cornell Univ., Ithaca, NY (United Sta↗

Atmospheric wind energization of ocean weather

Ocean weather comprises vortical and straining mesoscale motions, which play fundamentally different roles in the ocean circulation and climate system. Vorticity determines the movement of major ocean currents and gyres. Strain contributes to frontogenesis and the deformation of water masses, driving much of the mixing and vertical transport in the upper ocean. While recent studies have shown that interactions with the atmosphere damp the ocean’s mesoscale vortices O(100) km in size, the effect of winds on straining motions remains unexplored. Here, we derive a theory for wind work on the ocean’s vorticity and strain. Using satellite and model data, we discover that wind damps strain and vorticity at an equal rate globally, and unveil striking asymmetries based on their polarity. Subtropical winds damp oceanic cyclones and energize anticyclones outside strong current regions, while subpolar winds have the opposite effect. A similar pattern emerges for oceanic strain, where subtropical convergent flow is damped along the west-equatorward east-poleward direction and energized along the east-equatorward west-poleward direction. These findings reveal energy pathways through which the atmosphere shapes ocean weather.

54 ENVIRONMENTAL SCIENCES↗

Gaps and ways forward in atmospheric blocking and extreme weather research

Atmospheric blocking often results in significant weather extremes, such as heatwaves, droughts, cold spells, and floods in mid-latitude regions. However, the physical processes behind blocking and its response to climate change are not well understood, which undermines predictions and decision-making for climate mitigation and adaptation. As a phenomenon with a timescale at the interface of weather and climate, blocking interacts with various elements of the climate system and connects short-term weather events to long-term climate extremes. Understanding atmospheric blocking is crucial, given that climate change may impact its frequency, duration, and geographic distribution. This perspective discusses new experimental approaches to improve our understanding of this important subject.

Wang, Lei [Purdue University, West Lafayette, IN (↗

Aerosol effective radiative forcing increases Earth’s energy imbalance in recent decades

Earth’s energy imbalance (EEI), a key driver of climate change, has risen markedly over the last two decades. Greenhouse gas forcing, aerosol forcing, and cloud feedback all contribute to this increase. However, the role of aerosol forcing, particularly effective radiative forcing through aerosol-cloud interactions (ERF ACI ), remains highly uncertain and closely intertwined with cloud feedback. Here, we estimate shortwave ERF ACI using satellite observations and show its importance as a contributor to the EEI increase between 2003 and 2023. ERF ACI exhibits a significant increasing trend of 0.33 Wm −2 /decade averaged over oceans between 60°S and 60°N. The increasing trend of ERF ACI stems from a global decline in cloud droplet number concentration driven by decreasing anthropogenic aerosol emissions. It is of similar magnitude as the combined instantaneous forcing trend from greenhouse gases and aerosol-radiation interactions. Our results may reduce the gap between the simulated and observed shortwave contribution to EEI trends while implying a weak shortwave cloud feedback. They also suggest a stronger aerosol forcing than the multi-model mean.

Yuan, Tianle [NASA Goddard Space Flight Center (GS↗

Conformally coated scaffold design using water-tolerant Pr 1.8 Ba 0.2 NiO 4.1 for protonic ceramic electrochemical cells with 5,000-h electrolysis stability

Protonic ceramic electrochemical cells (PCECs) have potential as long-duration energy storage systems. However, their operational stability is limited under industrially relevant conditions due to the intrinsic chemical instability of doped barium cerate-based electrolytes and oxygen electrodes against H 2 O, as well as the poor electrode–electrolyte interfacial contact. Here, in this study, we present a conformally coated scaffold (CCS) design to comprehensively address these issues. A porous proton-conducting scaffold is constructed and conformally coated with Pr 1.8 Ba 0.2 NiO 4.1 electrocatalyst, which has high chemical stability against H 2 O, triple conductivity and hydration capability, and protects vulnerable electrolytes from H 2 O. The CCS structure consolidates the electrode–electrolyte interfacial bonding to enable fast proton transfer in the percolated network. This design enables PCECs to reach electrolysis stability for 5,000 h at −1.5 A cm−2 and 600 °C in 40% H 2 O. This work provides a general strategy to stabilize PCECs and offers guidance for designing resilient and stable solid-state energy storage systems.

Protonic ceramic electrochemical cell↗