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.

At least 19 records

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↗

General analytic results for nonlinear waves and solitons in molecular clouds

We study nonlinear wave phenomena in self-gravitating fluid systems, with a particular emphasis on applications to molecular clouds. This paper presents analytical results for one spatial dimension. We show that a large class of physical systems can be described by theories with a 'charge density' q(rho); this quantity replaces the density on the right-hand side of the Poisson equation for the gravitational potential. We use this formulation to prove general results about nonlinear wave motions in self-gravitating systems. We show that in order for stationary waves to exist, the total charge (the integral of the charge density over the wave profile) must vanish. This 'no-charge' property for solitary waves is related to the capability of a system to be stable to gravitational perturbations for arbitrarily long wavelengths. We find necessary and sufficient conditions on the charge density for the existence of solitary waves and stationary waves. We study nonlinear wave motions for Jeans-type theories (where q(rho) = rho-rho(sub 0)) and find that nonlinear waves of large amplitude are confined to a rather narrow range of wavelengths. We also study wave motions for molecular clouds threaded by magnetic fields and show how the allowed range of wavelengths is affected by the field strength. Since the gravitational force in one spatial dimension does not fall off with distance, we consider two classes of models with more realistic gravity: Yukawa potentials and a pseudo two-dimensional treatment. We study the allowed types of wave behavior for these models. Finally, we discuss the implications of this work for molecular cloud structure. We argue that molecular clouds can support a wide variety of wave motions and suggest that stationary waves (such as those considered in this paper) may have already been observed.

Adams, Fred C.↗

Radiative Transfer Simulations of the Two-Dimensional Ocean Glint Reflectance and Determination of the Sea Surface Roughness

An optimized discrete-ordinate radiative transfer model (DISORT3) with a pseudo-two-dimensional bidirectional reflectance distribution function (BRDF) is used to simulate and validate ocean glint reflectances at an infrared wavelength (1036 nm) by matching model results with a complete set of BRDF measurements obtained from the NASA cloud absorption radiometer (CAR) deployed on an aircraft. The surface roughness is then obtained through a retrieval algorithm and is used to extend the simulation into the visible spectral range where diffuse reflectance becomes important. In general, the simulated reflectances and surface roughness information are in good agreement with the measurements, and the diffuse reflectance in the visible, ignored in current glint algorithms, is shown to be important. The successful implementation of this new treatment of ocean glint reflectance and surface roughness in DISORT3 will help improve glint correction algorithms in current and future ocean color remote sensing applications.

glint↗

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↗

Investgation of radiative interactions in high-speed entrance region flows

The influence of radiative energy transfer on the entrance region flow is investigated under supersonic and subsonic flow conditions in a channel. Two-dimensional compressible Navier-Stokes equations are solved numerically in conjunction with the radiative flux equations. The channel walls are assumed to be black. Nongray as well as pseudo gray gas models are used to represent the absorption-emission characteristic of the medium. The participating species considered are different amounts of water vapor in water vapor-air mixtures. Results obtained for different flow conditions indicate that the radiative interaction can moderately influence the overall energy transfer, but the flowfield is not changed significantly.

Tiwari, S. N.↗

On the Two-Dimensionalization of Homogeneous Rotating Turbulence

Large-eddy simulation of the incompressible Navier-Stokes equations has been used to examine the long-time development of initially isotropic turbulence subjected to solid-body rotation. The simulations were carried out using a pseudo-spectral method with 128 x 128 x 512 collocation points in a computational domain that is four times larger along the rotation axis than in the other directions; subgrid-scale motions were parameterized using a spectral eddy viscosity model modified for system rotation. Simulation results show that the correlation length along the rotation am's of velocities orthogonal to the rotation vector exhibits rapid growth while the integral length-scale of velocities aligned with the rotation axis is relatively unaffected by rotation. Examination of the energy spectrum of two-dimensional, two-component motions indicates the presence of an inverse cascade of energy. System rotation also causes an alignment of vorticity along the rotation axis with relatively stronger cyclonic vorticity than anticyclonic. The onset of anisotropic effects are well characterized by Rossby numbers defined in terms of both macroscopic and microscopic quantities.

Squires, K. D.↗

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↗

Computational vortex dynamics in two and three dimensions

A discussion of some solutions and outstanding problems of dissipationless and nearly dissipationless flows in two and three dimensions is presented. In particular, two-dimensional steady-state solutions, linear stability and nonlinear evolutions including axisymmetrization of isolated vortex regions and merger or condensation of like-signed vorticity regions, are reviewed. These results are obtained from pseudo-spectral codes, contour dynamical (CD) codes and a new Hamiltonian moment-model of the Euler equations

Melander, Mogens V.↗

Numerical Analysis of the Sea State Bias for Satellite Altimetry

Theoretical understanding of the dependence of sea state bias (SSB) on wind wave conditions has been achieved only for the case of a unidirectional wind-driven sea. Recent analysis of Geosat and TOPEX altimeter data showed that additional factors, such as swell, ocean currents, and complex directional properties of realistic wave fields, may influence SSB behavior. Here we investigate effects of two-dimensional multimodal wave spectra using a numerical model of radar reflection from a random, non-Gaussian surface. A recently proposed ocean wave spectrum is employed to describe sea surface statistics. The following findings appear to be of particular interest: (1) Sea swell has an appreciable effect in reducing the SSB coefficient compared with the pure wind sea case but has less effect on the actual SSB owing to the corresponding increase in significant wave height. (2) Hidden multimodal structure (the two-dimensional wavenumber spectrum contains separate peaks, for swell and wind seas, while the frequency spectrum looks unimodal) results in an appreciable change of SSB. (3) For unimodal, purely wind-driven seas, the influence of the angular spectral width is relatively unimportant; that is, a unidirectional sea provides a good qualitative model for SSB if the swell is absent. (4) The pseudo wave age is generally much better fo parametrization the SSB coefficient than the actual wave age (which is ill-defined for a multimodal sea) or wind speed. (5) SSB can be as high as 5% of the significant wave height, which is significantly greater than predicted by present empirical model functions tuned on global data sets. (6) Parameterization of SSB in terms of wind speed is likely to lead to errors due to the dependence on the (in practice, unknown) fetch.

Glazman, R. E.↗

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↗

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↗

Thermochemical nonequilibrium and radiative interactions in supersonic hydrogen-air combustion

The two-dimensional, elliptic Navier-Stokes equations are used to investigate supersonic flows with nonequilibrium chemistry and thermodynamics, coupled with radiation, for hydrogen-air systems. The chemistry source term in the species equation is treated implicitly to alleviate the stiffness associated with fast reactions. The explicit, unsplit MacCormack finite-difference scheme is used to advance the governing equations in time, until convergence is achieved. The specific problem considered is the premixed, expanding flow in a supersonic nozzle. The reacting flow consists of seven species, one of which is the inert N2 molecule. The thermal state of the gas is modeled with one translational-rotational temperature and five vibrational temperatures. The harmonic oscillator model is used in the formulation for vibrational relaxation. The tangent slab approximation is used in the radiative flux formulation. A pseudo-gray model is used to represent the absorption-emission characteristics of the participating species. Results obtained for specific conditions indicate the presence of nonequilibrium in the expansion region. This reduces the radiative interactions and can have a significant influence on the flowfield.

Chandrasekhar, R.↗

Investigation of chemically reacting and radiating supersonic internal flows

The two-dimensional spatially elliptic Navier-Stokes equations are used to investigate the chemically reacting and radiating supersonic flow of the hydrogen-air system between two parallel plates and in a channel with a ten degree compression-expansion ramp at the lower boundary. The explicit unsplit finite-difference technique of MacCormack is used to advance the governing equations in time until convergence is achieved. The chemistry source term in the species equation is treated implicitly to alleviate the stiffness associated with fast reactions. The tangent slab approximation is employed in the radiative flux formation. Both pseudo-gray and nongray models are used to represent the absorption characteristics of the participating species. Results obtained for specific conditions indicate that the radiative interaction can have a significant influence on the flow field.

Mani, M.↗

Numerical solutions of chemically reacting and radiating flows

The two-dimensional spatially elliptic Navier-Stokes equations are used to investigate the chemically reacting and radiating supersonic flow of the hydrogen-air system between two parallel plates and in a channel with a ten-degree compression-expansion ramp at the lower boundary. The explicit unsplit finite-difference technique of MacCormack is used to advance the governing equations in time until convergence is achieved. The chemistry source term in the species equation is treated implicitly to alleviate the stiffness associated with fast reactions. The tangent slab approximation is employed in the radiative flux formulation. Both pseudo-gray and nongray models are used to represent the absorption-emission characteristics of the participating species. Results obtained for specific conditions indicate that the radiative interaction can have a significant influence on the flow field.

Mani, M.↗

An explicit Runge-Kutta method for turbulent reacting flow calculations

The paper presents a numerical method for the solution of the conservation equations governing steady, reacting, turbulent viscous flow in two-dimensional geometries, in both Cartesian and axisymmetric coordinates. These equations are written in Favre-averaged form and closed with a first order model. A two-equation K-epsilon model, where low Reynolds number and compressibility effects are included, and a modified eddy-break up model are used to simulate fluid mechanics turbulence, chemistry and turbulence-combustion interaction. The solution is obtained by using a pseudo-unsteady method with improved perturbation propagation properties. The equations are discretized in space by using a finite volume formulation. An explicit multi-stage dissipative Runge-Kutta algorithm is then used to advance the flow equations in the pseudo-time. The method is applied to the computation of both diffusion and premixed turbulent reacting flows. The computed temperature distributions compare favorably with experimental data.

Boretti, A. A.↗

The relation between ferroelasticity and superconductivity

The high-temperature superconductivity is explained widely by the layered crystal structure. The one- and two-dimensional subsystems and their interaction are investigated here. It is assumed that the high-T(sub c) superconductivity takes place in the two-dimensional subsystem and the increase of the phase transition temperature from 60 K up to 90 K is the consequence of turning on the influence of one-dimensional chains. The interaction between the two subsystems is transferred along the c axis by the phonons of breathing mode, which causes the hybridization of the electronic bonds between these subsystems. The experimental works indicate that the existence of both the chains Cu(1)-O and their interaction with the superconducting plane of Cu(2)-O modify the temperature of the transition to the superconducting state. It is seen from the neutron scattering data that the rates of the interatomic distance dependencies on temperature are changed around 140 K and 90 K. The 'zig-zag' order in Cu(1)-O chains has been postulated but, on the other hand, the vibrations with a large amplitude only were reported. The bi-stabilized situation of the oxygen ions can be caused by the change of distance between these ions and the Ba ions. It leads to the appearance of a two-well potential. Its parameters depend on temperature and the dynamics of the oxygen ions' movement. They can induce the antipolar order, which can be, however, more or less chaotic. The investigation of the ferroelastic properties of Y-Ba-Cu-O samples lead to the conclusion that they are related to jumps of ions inside the given chain and not to a diffusion between different sites in the ab plane. Researchers deduce thus that the fluctuating oxygen ions from these chains create dipoles in the ab plane. They can be described with the pseudo-spin formalism/ - Pauli matrices/. The system can be described with the Ising model. The pseudo-spins interact with phonons and influence the superconductivity in the second subsystem.

Molak, A.↗

Comparison between analytical models and finite-difference simulations in transmission-line tap resistors and L-type cross-Kelvin resistors

Approximate analytical models of the transmission-line tap resistor and the cross-Kelvin resistor are compared with computer-simulated pseudo-three-dimensional resistor network models. The analytical formulas are in good agreement with the simulations over a useful range of parameters and are readily applied to the extraction of the contact resistivity and the sheet resistances of the semiconducting layer under and outside the contacts. The extraction procedure, which is easily implemented on a personal computer, is carried out for the case of alloyed AuGeNi/GaAs contacts, illustrating the importance of (1) distinguishing between layer sheet resistance under and outside the contacts and (2) considering two-dimensional current flow in the semiconducting layer.

Scorzoni, Andrea↗

Radiative interactions in chemically reacting supersonic internal flows

The two-dimensional, elliptic Navier-Stokes equations are used to investigate supersonic flows with finite-rate chemistry and radiation for hydrogen-air systems. The chemistry source terms in the species equation is treated implicitly to alleviate the stiffness associated with fast reactions. The explicit, unsplit MacCormack finite-difference scheme is used to advance the governing equations in time, until convergence is achieved. The specific problem considered is the premixed flow in a channel with a ten-degree compression ramp. Three different chemistry models are used, accounting for increasing number of reactions and participating species. Two chemistry models assure nitrogen as inert, while the third model accounts for nitrogen reactions and NO(x) formation. The tangent slab approximation is used in the radiative flux formulation. A pseudo-gray model is used to represent the absorption-emission characteristics of the participating species. Results obtained for specific conditions indicate that the radiative interactions vary substantially, depending on reactions involving HO2 and NO species and that this can have a significant influence on the flowfield.

Tiwari, S. N.↗