Search NASA⌕ Search

SEARCH · Search NASA

Results for “pseudo two-dimensional model”

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.

16 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Simulation of Impedance Changes with Aging in Lithium Titanate-based Cells Using Physics-Based Dimensionless Modeling

Quantifying aging effects in lithium-ion cells with chemistries that have a flat open circuit potential is challenging. We implement a physics-based electrochemical model to track changes in the electrochemical impedance response of lithium titanate-based cells. Frequency domain equations of a pseudo two-dimensional model are made dimensionless, and the corresponding non-dimensional parameters are estimated using a Levenberg-Marquardt routine. The model weighs the relative contributions of changes in diffusion, ionic conduction within the electrolyte phase against solid phase electronic conduction towards cell aging. Solid-phase diffusion, charge transfer resistance and double layer capacitance at the solid-liquid interface are accounted for in the particle impedance. The estimation routine tracks dimensionless parameters using accelerated cycling data from full cells over 1000 cycles. The model can be deployed within a short time for state estimation using physics-based models without requiring prior knowledge of the battery chemistry, format, or capacity.

25 ENERGY STORAGE↗
DOE OSTI2024

Rapid Inverse Parameter Inference Using Physics-Informed Neural Network

As Li-ion batteries become more essential in today's economy, tools need to be developed to accurately and rapidly diagnose a battery's internal state-of-health. Using a Li-ion battery's (high-rate) voltage response, it is proposed to determine a battery's internal state through Bayesian calibration. However, Bayesian calibration is notoriously slow and requires thousands of model runs. To accelerate parameter inference using Bayesian calibration, a surrogate model is developed to replace the underlying physics-based Li-ion model. Developing a surrogate model for rapid Bayesian calibration analysis is discussed for both the single particle model (SPM) and the pseudo two-dimensional (P2D) model. Surrogate models are constructed using physics-informed neural networks (PINNs) that encode the influence of internal properties on observed voltage responses. In practice, a neural network can be trained by: 1) using simulation results of the physics-based model (i.e., a data-loss approach); 2) using the residuals of the governing equations themselves (i.e., a physics-loss approach); or 3) using a combination of simulation results and governing equation residuals. In the present work, PINNs are developed using a variety of training losses and neural network architectures. In this analysis, it is shown that a PINN surrogate model can be reliably trained with only physics-informed loss. However, using a coupled data-informed and physics-loss approach produced the most accurate PINNs.

Bayesian calibration↗
DOE OSTI2022

A Robust Numerical Treatment of Solid-Phase Diffusion in Pseudo Two-Dimensional Lithium-Ion Battery Models

Solid-phase diffusion in active materials of lithium-ion batteries significantly affects charging and safety-related behavior of lithium-ion batteries. Therefore, it is essential to develop an efficient and robust numerical algorithm for solving solid-phase diffusion equations in physics-based battery models. In this work, we discuss the origins of numerical instabilities that can occur when solving the solid-phase diffusion equations using iterative methods. Then, in order to resolve such issues, we propose a simple numerical treatment to the surface flux term of discretized solid-phase diffusion equations. To demonstrate its numerical robustness, the proposed method is implemented into a pseudo two-dimensional (P2D) physics-based battery model and simulations are conducted at wide ranges of operating conditions. Even with extremely poor initial guesses for the Li+ concentrations of the active materials, computations using the proposed method do not diverge and the their computational speeds are comparable to those with conventional initial guesses. Comprehensive tests of the proposed method are also performed with a dynamic current profile based on US06 driving profile and a multi-stage charging profile with very high initial C-rate (12C).

battery modeling↗
DOE OSTI2022

A Tanks-in-Series Approach to Estimate Parameters for Lithium-Ion Battery Models

Advanced Battery Management Systems (BMS) play a vital role in monitoring, predicting, and controlling the performance of lithium-ion batteries. BMS employing sophisticated electrochemical models can help increase battery cycle life and minimize charging time. However, in order to realize the full potential of electrochemical model-based BMS, it is critical to ensure accurate predictions and proper model parameterization. The accuracy of the predictions of an electrochemical model is dependent on the accuracy of its parameters, the values of which might change with battery cycling and aging. Parameter estimation for an electrochemical model is generally challenging due to the nonlinear nature and computational complexity of the model equations. To this end, this work utilizes the recently proposed Tanks-in-Series model for Li-ion batteries (J.Electrochem. Soc., 167, 013534 (2020)) to perform parameter estimation. The Tanks-in-Series approach allows for substantially faster parameter estimation compared to the original pseudo two-dimensional (p2D) model. The objective of this work is thus to demonstrate the gain in computational efficiency from the Tanks-in-Series approach. A sensitivity analysis of model parameters is also performed to benchmark the fidelity of the Tanks-in-Series model.

25 ENERGY STORAGE↗
DOE OSTI2022

Physics-Informed Neural Network Modeling of Li-Ion Batteries: Preprint

Li-ion batteries (LIB) are a promising solution to enable storage of intermittent energy sources due to their high energy density. However, LIBs are known to significantly degrade after about 1000 charge-discharge cycles. LIBs degrade following different degradation modes and at a rate that depends on the operating conditions (external temperature, load). To plan the installation of batteries, appropriate understanding and prediction capabilities of their lifecycle is needed. High-fidelity numerical models of LIBs such as the pseudo-two-dimensional (P2D) model have been shown to accurately represent the charge-discharge-cycle of an LIB given ac- curate choice of the physical parameters. Given the large number of P2D parameters, adjusting them using forward runs is intractable. This work describes the development of a physics-informed neural network (PINN) as a surrogate substitute of the P2D model that captures parameter dependence. The PINN is advantageous as it needs little to no data, and can naturally encode the dependence of every model parameter. Here, a specific training procedure is adopted to efficiently cover parameter space, handle model stiffness and enforce boundary conditions. The trained PINN is validated against numerical solutions of the P2D model, and its applicability to battery degradation modeling is discussed.

battery degradation↗
DOE OSTI2024

Theoretical design and performance of three-dimensional, pillared FeS 2 cathodes

Three-dimensional (3D) electrode design can provide improved capacities and rate capabilities over conventional two-dimensional electrodes by enhancing electrical and ionic transport. Here, expanding upon our previous modeling efforts for conversion chemistry lithium-ion batteries, we develop a pseudo-four-dimensional (P4D) approach that is subsequently used to investigate the design of a pillared FeS 2 electrode. The model considers transport in three dimensions with an additional “fourth” dimension corresponding to the solid-state lithium transport within the active material particles. Additionally, we allow for expansion of the active material during the conversion reaction to understand how internal stresses impact the electrochemical performance of the cell. By optimizing the model with respect to areal capacity, we are able to predict areal capacities up to 16.8 mAh/cm 2 for an areal current density of 1.78 mA/cm 2 and a 103% improvement for the three-dimensional electrodes over planar electrodes of equal volume. Despite the promising results, our simulations suggest that 3D design may be difficult for conversion cathode materials due to the large internal stresses that arise during conversion. Nevertheless, the model is robust and adaptable to other materials that may be more suitable for 3D electrodes due to a lesser change in volume during discharge.

Conversion cathode materials↗
DOE OSTI2024

Analytical Model for Atomic Relaxation in Twisted Moiré Materials

By virtue of being atomically thin, the electronic properties of heterostructures built from two-dimensional materials are strongly influenced by atomic relaxation. The atomic layers behave as flexible membranes rather than rigid crystals. Here we develop an analytical theory of lattice relaxation in twisted moiré materials. We obtain analytical results for the lattice displacements and corresponding pseudo gauge fields, as a function of twist angle. We benchmark our results for twisted bilayer graphene and twisted WSe 2 bilayers using large-scale molecular dynamics simulations. Our single-parameter theory is valid in graphene bilayers for twist angles 𝜃 ≳ 0.7°, and in twisted WSe 2 for 𝜃 ≳ 1.6°. Furthermore, we also investigate how relaxation alters the electronic structure in twisted bilayer graphene, providing a simple extension to the continuum model to account for lattice relaxation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Magnetism in the two-dimensional dipolar XY model

Here, motivated by a recent experiment on a square-lattice Rydberg atom array realizing a long-range dipolar XY model [C. Chen et al., Nature (London) 616, 691 (2023)], we numerically study the model's equilibrium properties. We obtain the phase diagram, critical properties, entropies, variance of the magnetization, and site-resolved correlation functions. We consider both ferromagnetic and antiferromagnetic interactions and apply quantum Monte Carlo and pseudo-Majorana functional renormalization group techniques, generalizing the latter to a U⁡(1) symmetric setting. Our simulations perform extensive thermometry in dipolar Rydberg atom arrays and establish conditions for adiabaticity and thermodynamic equilibrium. On the ferromagnetic side of the experiment, we determine the entropy per particle S/N≈0.5, close to the one at the critical temperature, S c /N=0.585⁢(15). The simulations suggest the presence of an out-of-equilibrium plateau at large distances in the correlation function, thus motivating future studies on the nonequilibrium dynamics of the system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Two-Dimensional Ferroelectric Altermagnets: From Model to Material Realization

Multiferroic altermagnets offer new opportunities for magnetoelectric coupling and electrically tunable spintronics. However, due to intrinsic symmetry conflicts between altermagnetism and ferroelectricity, achieving their coexistence, known as ferroelectric altermagnets (FEAM), remains an outstanding challenge, especially in two-dimensional (2D) systems. Here, we propose a universal, symmetry-based design principle for 2D FEAM, supported by tight-binding models and first-principles calculations. We show that lattice distortions can break spin equivalence and introduce the necessary rotation-related symmetry, enabling altermagnetism with electrically reversible spin splitting. Guided by this framework, we identify a family of 2D vanadium oxyhalides and sulfide halides as promising FEAM candidates. In these compounds, pseudo Jahn-Teller distortions and Peierls-like dimerization cooperatively establish the required symmetry conditions. Here, we further propose the magneto-optical Kerr effect as an experimental probe to confirm FEAM and its electric spin reversal. Furthermore, our findings provide a practical framework for 2D FEAM and advancing electrically controlled spintronic devices.

2D materials↗
DOE OSTI2022

A Simulation-based Method for Correcting Mode Coupling in CMB Angular Power Spectra

Modern cosmic microwave background (CMB) analysis pipelines regularly employ complex time-domain filters, beam models, masking, and other techniques during the production of sky maps and their corresponding angular power spectra. However, these processes can generate couplings between multipoles from the same spectrum and from different spectra, in addition to the typical power attenuation. Within the context of pseudo-C ℓ based, MASTER-style analyses, the net effect of the time-domain filtering is commonly approximated by a multiplicative transfer function, F ℓ , that can fail to capture mode mixing and is dependent on the spectrum of the signal. To address these shortcomings, we have developed a simulation-based spectral correction approach that constructs a two-dimensional transfer matrix, ${J}_{{\ell }{\ell }^{\prime} }$, which contains information about mode mixing in addition to mode attenuation. We demonstrate the application of this approach on data from the first flight of the Spider balloon-borne CMB experiment.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Chemical state mapping of simulant Chernobyl lava-like fuel containing material using micro-focused synchrotron X-ray spectroscopy

Uranium speciation and redox behaviour is of critical importance in the nuclear fuel cycle. X-ray absorption near-edge spectroscopy (XANES) is commonly used to probe the oxidation state and speciation of uranium, and other elements, at the macroscopic and microscopic scale, within nuclear materials. Two-dimensional (2D) speciation maps, derived from microfocus X-ray fluorescence and XANES data, provide essential information on the spatial variation and gradients of the oxidation state of redox active elements such as uranium. In the present work, we elaborate and evaluate approaches to the construction of 2D speciation maps, in an effort to maximize sensitivity to the U oxidation state at the U L 3 -edge, applied to a suite of synthetic Chernobyl lava specimens. Our analysis shows that calibration of speciation maps can be improved by determination of the normalized X-ray absorption at excitation energies selected to maximize oxidation state contrast. The maps are calibrated to the normalized absorption of U L 3 XANES spectra of relevant reference compounds, modelled using a combination of arctangent and pseudo-Voigt functions (to represent the photoelectric absorption and multiple-scattering contributions). We validate this approach by microfocus X-ray diffraction and XANES analysis of points of interest, which afford average U oxidation states in excellent agreement with those estimated from the chemical state maps. This simple and easy-to-implement approach is general and transferrable, and will assist in the future analysis of real lava-like fuel-containing materials to understand their environmental degradation, which is a source of radioactive dust production within the Chernobyl shelter.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Application of mesh refinement to relativistic magnetic reconnection

During relativistic magnetic reconnection, antiparallel magnetic fields undergo a rapid change in topology, releasing a large amount of energy in the form of non-thermal particle acceleration. This work explores the application of mesh refinement to 2D reconnection simulations to efficiently model the inherent disparity in length-scales. We have systematically investigated the effects of mesh refinement and determined necessary modifications to the algorithm required to mitigate non-physical artifacts at the coarse–fine interface. We have used the ultrahigh-order pseudo-spectral analytical time-domain Maxwell solver to analyze how its use can mitigate the numerical dispersion that occurs with the finite-difference time-domain (or “Yee”) method. Absorbing layers are introduced at the coarse–fine interface to eliminate spurious effects that occur with mesh refinement. We also study how damping the electromagnetic fields and current density in the absorbing layer can help prevent the non-physical accumulation of charge and current density at the coarse–fine interface. Using a mesh refinement ratio of 8 for two-dimensional magnetic reconnection simulations, we obtained good agreement with the high-resolution baseline simulation, using only 36% of the macroparticles and 71% of the node-hours needed for the baseline. The methods presented here are especially applicable to 3D systems where higher memory savings are expected than in 2D, enabling comprehensive, computationally efficient 3D reconnection studies in the future.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Spatially Quasi-Periodic Water Waves of Infinite Depth

Abstract We formulate the two-dimensional gravity-capillary water wave equations in a spatially quasi-periodic setting and present a numerical study of solutions of the initial value problem. We propose a Fourier pseudo-spectral discretization of the equations of motion in which one-dimensional quasi-periodic functions are represented by two-dimensional periodic functions on a torus. We adopt a conformal mapping formulation and employ a quasi-periodic version of the Hilbert transform to determine the normal velocity of the free surface. Two methods of time-stepping the initial value problem are proposed, an explicit Runge–Kutta (ERK) method and an exponential time-differencing (ETD) scheme. The ETD approach makes use of the small-scale decomposition to eliminate stiffness due to surface tension. We perform a convergence study to compare the accuracy and efficiency of the methods on a traveling wave test problem. We also present an example of a periodic wave profile containing vertical tangent lines that is set in motion with a quasi-periodic velocity potential. As time evolves, each wave peak evolves differently, and only some of them overturn. Beyond water waves, we argue that spatial quasi-periodicity is a natural setting to study the dynamics of linear and nonlinear waves, offering a third option to the usual modeling assumption that solutions either evolve on a periodic domain or decay at infinity.

Wilkening, Jon (ORCID:0000000327827596)↗
DOE OSTI2024

Probabilistic Context Neighborhood model for lattices

Here we present the Probabilistic Context Neighborhood model designed for two-dimensional lattices as a variation of a Markov random field assuming discrete values. In this model, the neighborhood structure has a fixed geometry but a variable order, depending on the neighbors’ values. Our model extends the Probabilistic Context Tree model, originally applicable to one-dimensional space. It retains advantageous properties, such as representing the dependence neighborhood structure as a graph in a tree format, facilitating an understanding of model complexity. Furthermore, we adapt the algorithm used to estimate the Probabilistic Context Tree to estimate the parameters of the proposed model. We illustrate the accuracy of our estimation methodology through simulation studies. Additionally, we apply the Probabilistic Context Neighborhood model to spatial real-world data, showcasing its practical utility.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Influence of Terminal Carboxyl Groups on the Structure and Reactivity of Functionalized m-Carboranethiolate Self-Assembled Monolayers

The structure and function of self-assembled monolayers (SAMs) at the nanoscale are determined by the steric and electronic effects of their building blocks. Carboranethiol molecules form pristine monolayers that provide tunable two-dimensional systems to probe lateral and interfacial interactions. Additional ω-functionality, such as carboxyl groups, can be introduced to change the properties of the exposed surfaces. Here, two geometrically similar isomeric m-carborane analogues of m-mercaptobenzoic acid, 1-COOH-7-SH-1,7-C 2 B 10 H 10 and racem-1-COOH-9-SH-1,7-C 2 B 10 H 10 , are characterized and their SAMs on Au{111} are examined. The latter isomer belongs to the rare group of chiral cage molecules and becomes, to our knowledge, the first example assembled on Au{111}. Although different in symmetry, molecules of both isomers assemble into similar hexagonal surface patterns. The nearest-neighbor spacing of 8.4 ± 0.4 Å is larger than that of non-carboxylated isomers, consistent with the increased steric demands of the carboxyl groups. Computational modeling reproduced this spacing and suggests a tilt relative to the surface normal. However, tilt domains are not observed experimentally, suggesting the presence of strong lateral interactions. Analyses of the influence of the functional groups through the pseudo-aromatic m-carborane skeleton showed that the thiol group attached to either carbon or boron atoms increases the carboxyl group acidity in solution. In contrast, the acidity of the exposed carboxyl group in the SAMs decreases upon surface attachment; computational analyses suggest that the driving force of this shift is the dielectric of the environment in the monolayer as a result of confined intermolecular interactions, proximity to the Au surface, and partial desolvation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Higher-order topological corner and bond-localized modes in magnonic insulators

We theoretically investigate a two-dimensional decorated honeycomb lattice framework to realize a second-order topological magnon insulator (SOTMI) phase featuring distinct corner-localized modes. Here, our study emphasizes the pivotal role of spin-magnon mapping in characterizing bosonic topological properties, which exhibit differences from their fermionic counterparts. We employ a symmetry indicator topological invariant to identify and characterize this SOTMI phase, particularly for systems respecting time-reversal and rotational symmetry. Using a spin model defined on a honeycomb lattice geometry, we demonstrate that introducing “kekulé” type distortions yields a topological phase. In contrast, “antikekulé” distortions result in a nontopological magnonic phase. The presence of kekulé distortions manifests in two distinct topologically protected bosonic corner modes—an intrinsic and a pseudo, based on the specific edge terminations. On the other hand, antikekulé distortions give rise to Tamm/Shockley type bond-localized boundary modes, which are nontopological and reliant on particular edge termination. We further investigate the effects of random out-of-plane exchange anisotropy disorder on the robustness of these bosonic corner modes. The distinction between SOTMIs and their fermionic counterparts arises due to the system-specific magnonic onsite energies, a crucial feature often overlooked in prior literature. Our study unveils exciting prospects for engineering higher-order topological phases in magnon systems and enhances our understanding of their unique behavior within decorated honeycomb lattices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗