Search NASA⌕ Search

SEARCH · Search NASA

Results for “Discontinuous fields”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Meshfree Methods for Modeling Chemo-Mechanical Cathode Cracking in Li-Ion Batteries

Through repeated charging and discharging, the electrodes of a Li-ion battery experience cyclic loading causing swelling and contraction due to the movement of lithium, also known as intercalation. This research focuses specifically on the chemo-mechanical cracking in the cathode. Because cathode particles are comprised of many randomly-oriented grains, which have highly anisotropic material properties, the expansion and contraction is very non-uniform. As a result, stress concentrations tend to form between grains, which necessitates the modeling of crack propagation largely along grain boundaries. Chemo-mechanical models are simulated mainly by two methods: the cohesive zone model (CZM) and the continuous damage model (CDM). The CZM is more accurate at capturing the sharp discontinuities of a crack but is very computationally expensive and intractable for large-scale models as a result. Conversely, the CDM is easily computed but not well-suited to easily allow for discontinuous field variables, which are inherent across a crack. This study aims to improve the CDM's ability to capture discontinuous cracks. Current versions of the CDM use the finite element method (FEM), which is one of the most widely used approaches for spatial discretization. This research investigates the use of the Reproducing Kernel Particle Method (RKPM), a meshfree method, for spatial discretization and aims to achieve a chemo-mechanical crack propagation model that achieves both intermediate accuracy and computational time between a CDM and CZM. To gain a better understanding of how the meshfree model compares to the current FEM model in capturing cathode crack propagation, further investigation is needed.

cathode cracking↗
DOE OSTI2023

Meshfree Multiphysics Damage Modeling of Li-ion Battery Materials

Through repeated charging and discharging cycles, the electrodes of a Li-ion battery experience significant swelling and contraction due to the movement of lithium, also known as intercalation. This research focuses specifically on the chemo-mechanical cracking in the cathode. Because cathode particles are comprised of many randomly-oriented grains, which have highly anisotropic material properties, the expansion and contraction is very non-uniform. As a result, stress concentrations tend to form between grains, which necessitates the modeling of crack propagation largely along grain boundaries and material interfaces. Chemo-mechanical damage models are generally simulated by one of two methods: the cohesive zone model (CZM) and the continuous damage model (CDM). The CZM is more accurate at capturing the sharp discontinuities of a crack but is very computationally expensive and intractable for large-scale models as a result. Conversely, the CDM is easily computed but not well-suited to easily allow for discontinuous field variables, which are inherent across a crack. This study aims to improve the CDM's ability to capture discontinuous cracks. Current versions of the CDM use the finite element method (FEM), which is one of the most widely used approaches for spatial discretization. This research investigates the use of the Reproducing Kernel Particle Method (RKPM), a meshfree method, for spatial discretization and aims to achieve a coupled chemo-mechanical crack propagation model that enhances accuracy while maintaining high computational efficiency. This model features a fully coupled, iterative electrochemistry solution, which informs the meshfree damage model of impending crack formation. To gain a better understanding of how the meshfree model compares to the current FEM model in capturing cathode crack propagation, further investigation is needed.

damage modeling↗
DOE OSTI2016

TURBULENCE-GENERATED PROTON-SCALE STRUCTURES IN THE TERRESTRIAL MAGNETOSHEATH

Recent results of numerical magnetohydrodynamic simulations suggest that in collisionless space plasmas, turbulence can spontaneously generate thin current sheets. These coherent structures can partially explain the intermittency and the non-homogenous distribution of localized plasma heating in turbulence. In this Letter, Cluster multi-point observations are used to investigate the distribution of magnetic field discontinuities and the associated small-scale current sheets in the terrestrial magnetosheath downstream of a quasi-parallel bow shock. It is shown experimentally, for the first time, that the strongest turbulence-generated current sheets occupy the long tails of probability distribution functions associated with extremal values of magnetic field partial derivatives. During the analyzed one-hour time interval, about a hundred strong discontinuities, possibly proton-scale current sheets, were observed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Exotic $\mathbb{Z}_N$ symmetries, duality, and fractons in 3+1-dimensional quantum field theory

Following our earlier analyses of nonstandard continuum quantum field theories, we study here gapped systems in 3+1 dimensions, which exhibit fractonic behavior. In particular, we present three dual field theory descriptions of the low-energy physics of the X-cube model. A key aspect of our constructions is the use of discontinuous fields in the continuum field theory. Spacetime is continuous, but the fields are not.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Neural Network-Enhanced Reproducing Kernel Particle Method for Image-Based Multiphysics Damage Modeling of Energy Storage Materials

Energy storage materials undergo significant stresses during charge/discharge cycling, which makes understanding their reliability and durability fundamental in predicting performance and service life. Strong electrochemical-mechanical coupling and highly anisotropic material properties contribute to the formation and propagation of micro-cracking, largely along material interfaces and grain boundaries. With microstructural images supplied by the National Renewable Energy Laboratory (NREL), image-based modeling techniques are used to represent the complex material microstructures that dictate the coupled physics of these systems. Traditional electrochemical-mechanical models rely on mesh-based finite element methods, which can lead to difficulties in capturing crack propagation due to mesh dependency. Additionally, commonly used damage models, such as the continuous damage model and the cohesive zone model, often have steep tradeoffs between discontinuous field accuracy and computational expense. In this work, a neural network-enhanced reproducing kernel particle method (NN-RKPM) [1] is leveraged to accurately capture damage and crack propagation throughout the material by learning the location, orientation, and sharpness of discontinuity while allowing for a coarser nodal distribution than that necessary for capturing sharp solution transitions using traditional mesh-based methods. NN-RKPM is used to inform how crack opening and closure in turn affect the coupled chemical equations and material microstructure. Reference: [1] Baek, J., Chen, J. S., Susuki, K., "Neural Network enhanced Reproducing Kernel Particle Method for Modeling Localizations," International Journal for Numerical Methods in Engineering, Vol. 123, pp 4422-4454, https://doi.org/10.1002/nme.7040, 2022.

damage modeling↗
DOE OSTI2014

Stark effect in ensembles of polar (0001) Al{sub 0.5}Ga{sub 0.5}N/GaN quantum dots and comparison with semipolar (11−22) ones

This work presents a continuous-wave photoluminescence study of Al{sub 0.5}Ga{sub 0.5}N/GaN quantum dots grown by ammonia-assisted molecular beam epitaxy on sapphire, either on the wurtzite polar (0001) or the semipolar (11−22) plane. Due to interface polarization discontinuities, the polar dots are strongly red-shifted by the Stark effect and emit in the visible range. Carrier injection screening of the polarization charges has been studied. A model relying on average dot heights and dot height variances, as measured by transmission electron microscopy, is proposed. It can account for the injection dependent luminescence energies and efficiencies. The electric field discontinuity deduced from the fittings is in good agreement with theoretical expectations for our barrier composition. On the contrary, semipolar quantum dot ensembles always emit above the gap of GaN strained to Al{sub 0.5}Ga{sub 0.5}N. Their luminescence linewidth is significantly lower than that of polar ones, and their energy does not shift with injection. Our study then confirms the expected strong decrease of the Stark effect for (11−22) grown (Al,Ga)N/GaN heterostructures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Phase field modeling of discontinuous dynamic recrystallization in hot deformation of magnesium alloys

We have developed a quantitative phase field model of discontinuous dynamic recrystallization (PF-DDRX) for investigating grain structure evolution and its impact on mechanical response in magnesium alloy during hot deformation. A set of isothermal compression tests were conducted by Gleeble thermo-mechanical simulator and used to determine PF-DDRX model parameters such as critical stresses of recrystallization, activation energies of dislocation, grain boundaries and recrystallization grain nucleation. In the simulations, the experimentally measured average grain size and distribution were used to generate initial microstructures. The results show that the typical ‘necklace’ microstructures observed in DDRX forms at the earlier stage of recrystallization. The inhomogeneity of microstructures first increases rapidly and then decreases gradually and reaches a constant with the increase of dynamic recrystallization volume fraction. Both predicted grain structure and mechanical response are in good agreement with experimental results. This demonstrates that the model parameters determined by experimental data are reliable for the developed PF-DDRX model to predict the grain structure evolution and flow behavior of AZ80 magnesium alloy under hot deformation. This work presents a path to combine experiments and modeling to give quantitative prediction and minimize experimental cost.

magnesium alloy, Phase field method, dynamic recry↗
DOE OSTI2025

On the numerical sensitivity of cellular automata grain structure predictions to large thermal gradients and cooling rates

Cellular automata (CA) models of as-solidified grain structure, originally developed and applied to casting, have become a common means of predicting grain structure resulting from Additive Manufacturing (AM) processes. The majority of these models are based on the decentered octahedron approach, which attempts to correct for the effect of grid anisotropy on the prediction of competitive solidification of dendritic grains. However, AM solidification occurs under cooling rates ($\dot{T}$) and thermal gradients (G) that are orders of magnitude larger than those encountered in casting, and no systematic investigation on the effect of the CA model cell size (Δx) and time step (Δt) on AM microstructure predictions has been performed. Here, in this study, such an investigation is first performed via simulation of individual grains of various crystallographic orientations with a fixed, unidirectional G, showing that CA prediction of the steady-state undercooling matched the expected values based on the interfacial response function at small G and deviated from the expected values at large G. Simulation of competitive growth of multiple grains showed a weakening of the predicted texture as G and Δx became large. Simulation of solidification under AM conditions, where G and $\dot{T}$ vary spatially across the melt pools, showed that not only does grain selection weaken and deviate from expectations at large Δx, but grains with crystallographic $\langle$100$\rangle$ aligned with the grid directions are more adversely affected by the temperature field discontinuities than grains with other crystallographic orientations. Despite the fact that the exact grain competition results depended on Δt, the overall texture development was notably less sensitive to Δt than Δx, provided that a reasonable value of Δt is selected based on the ratio of Δx to the maximum local solidification velocity in the simulation domain. Finally, from the directional solidification and AM simulation results, an analysis of computational cost compared to simulation resolution is performed based on an equation derived to quantify the relatively inaccuracy in grain selection based on the model and temperature field inputs. From this analysis, it is concluded that there is a need for algorithmic improvements to improve CA grain competition accuracy for large G processing conditions as sufficiently small Δx to resolve the necessary competition is intractable for many AM processing conditions.

36 MATERIALS SCIENCE↗
NASA NTRS1994

Introduction to Generalized Functions with Applications in Aerodynamics and Aeroacoustics

Generalized functions have many applications in science and engineering. One useful aspect is that discontinuous functions can be handled as easily as continuous or differentiable functions and provide a powerful tool in formulating and solving many problems of aerodynamics and acoustics. Furthermore, generalized function theory elucidates and unifies many ad hoc mathematical approaches used by engineers and scientists. We define generalized functions as continuous linear functionals on the space of infinitely differentiable functions with compact support, then introduce the concept of generalized differentiation. Generalized differentiation is the most important concept in generalized function theory and the applications we present utilize mainly this concept. First, some results of classical analysis, are derived with the generalized function theory. Other applications of the generalized function theory in aerodynamics discussed here are the derivations of general transport theorems for deriving governing equations of fluid mechanics, the interpretation of the finite part of divergent integrals, the derivation of the Oswatitsch integral equation of transonic flow, and the analysis of velocity field discontinuities as sources of vorticity. Applications in aeroacoustics include the derivation of the Kirchhoff formula for moving surfaces, the noise from moving surfaces, and shock noise source strength based on the Ffowcs Williams-Hawkings equation.

Farassat, F.↗
NASA NTRS1985

On the existence and stability conditions for mixed-hybrid finite element solutions based on Reissner's variational principle

The extensions of Reissner's two-field (stress and displacement) principle to the cases wherein the displacement field is discontinuous and/or the stress field results in unreciprocated tractions, at a finite number of surfaces ('interelement boundaries') in a domain (as, for instance, when the domain is discretized into finite elements), is considered. The conditions for the existence, uniqueness, and stability of mixed-hybrid finite element solutions based on such discontinuous fields, are summarized. The reduction of these global conditions to local ('element') level, and the attendant conditions on the ranks of element matrices, are discussed. Two examples of stable, invariant, least-order elements - a four-node square planar element and an eight-node cubic element - are discussed in detail.

Karlovitz, L. A.↗
DOE OSTI2026

Performance of Heterostructural TaC/AlGaN Schottky Diodes Based on First Principles Electronic Structure Properties

Advances in ultra-wide bandgap materials, such as high Al-content AlxGa1-xN (AlGaN), are essential for next generation power electronics, but the requirement for lattice matched substrates is currently a significant obstacle. Recently, conductive TaC has emerged as a promising virtual substrate for AlGaN heteroepitaxy, with wurtzite (0001) AlxGa1-xN lattice-matched to rocksalt (111) TaC at x ~ 0.5. Thus, understanding and controlling the electronic properties of the TaC/AlGaN interface is key for developing technological applications based on TaC/AlGaN devices. Using density functional theory and electronic structure calculations, we here investigate TaC/Al0.5Ga0.5N interfaces, where we include explicit alloy models in the slab calculations. We predict the Schottky barrier height and the electric field discontinuity resulting from interface charges. Considering all possible combinations of (Ta, C) substrate termination, (Al/Ga, N) nucleation, and (Al/Ga, N) polarity, we construct a chemical potential phase diagram to identify the stable interfaces that can be accessed through variation of the synthesis conditions. The predicted interface electronic properties are implemented in device performance simulations to demonstrate a practical design for a strain-free, high-efficiency TaC/AlGaN Schottky diode with a low barrier height and without interface charges, underscoring the potential of TaC as a substrate for ultra-wide bandgap devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1972

Transport of cosmic rays in the solar corona

A method by which energetic cosmic ray particles in the interplanetary medium can diffuse in solar longitude. The method consists of particle motion occurring along current sheets separating discontinuous field structures in the corona. These sheets can serve as pathways along which energetic particles drift at nearly their propagation speed. A model of the coronal magnetic field and the large scale field behavior are analyzed. A mechanism by which particles can travel across the coronal magnetic field and arrive at the earth is suggested.

Schatten, K. H.↗
NASA NTRS1974

Post-shock spikes: A new feature of proton and alpha enhancements associated with an interplanetary shock wave

Abrupt and prolonged enhancements in the intensities of 100 to approximately 2000 keV nucleon protons and alpha particles observed in interplanetary space are interpreted as particle populations confined between an interplanetary shock front and a magnetic field discontinuity. Prominent intensity spikes observed only below approximately 400 keV per charge for both protons and alpha particles several hours behind the shock front suggest that some fraction of the confined particles is accelerated by an energy per charge dependent process.

Gloeckler, G.↗
NASA NTRS1975

Modulation of low-energy cosmic rays

The relation between the diffusion coefficient of cosmic rays in the solar wind and the power spectrum of interplanetary magnetic field fluctuations, established in recent theories, is tested directly for low energy protons (below 80 MeV). In addition, an attempt is made to determine whether the particles are scattered by magnetic field discontinuities or by fluctuations between discontinuities. Predictions of a perturbation solution of the Fokker-Planck equation are compared with observations of the cosmic ray radial gradient. It is found that at energies between 40 and 80 MeV, galactic cosmic ray protons respond to changes in the predicted diffusion coefficients (i.e., the relationship under consideration holds at these low energies). The relation between changes in the proton flux and modulation parameters is best when the contribution of discontinuities is subtracted, which means that scattering is caused by fluctuations between discontinuities. There appears to be no distinct relation between changes in the modulation parameters and changes in the intensity of 20 to 40 MeV protons.

Sari, J. W.↗
NASA NTRS1977

Electromagnetic sounding of the moon using Apollo 16 and Lunokhod 2 surface magnetometer observations /preliminary results/

A technique of deep electromagnetic sounding of the moon using simultaneous magnetic-field measurements at two lunar surface sites is described. The method, used with the assumption that deep electrical conductivity is a function only of lunar radius, has the advantage of allowing calculation of the external driving field from two surface-site measurements only and therefore does not require data from a lunar orbiting satellite. A transient-response calculation is presented for the example of a magnetic-field discontinuity, measured simultaneously by Apollo 16 and Lunokhod 2 surface magnetometers.

Vanian, L. L.↗
NASA NTRS1978

Plasma electron observations in the wake of Venus

Two plasma regimes were observed in the Venus wake, one characterized by negative-going anisotropies accompanied by high bulk speeds and low magnetic variance, and the other characterized by positive-going anisotropies accompanied by low bulk speeds and high magnetic variance. The two modes alternate in time in a stochastic bimodal fashion as interplanetary magnetic field discontinuities propagate past Venus and switch on the positive anisotropies. The positive mode probably represents enhanced planetary ion pickup which results in instabilities, the damping of which energizes the electrons in the distant wake.

Yeates, C. M.↗
NASA NTRS1984

High-altitude observations of an intense inverted V event

Inverted-V events which generally occur in the pre-midnight sector over the auroral zone are frequently associated with reversals in the convection electric field. Such a reversal is observed by the University of Iowa quasispherical LEPEDEA on board ISEE 1 at an altitude of 13 RE on May 1, 1978. The bulk of the plasma shows a large shear over a five-minute interval. The associated change in the convection electric field is 5.1 mV/m. Large values of the field-aligned current are simultaneously detected. The potential structure appears to extend to the satellite location. Using a theoretical model, the field-aligned current due to the electric field discontinuity has been calculated. The magnitude of the parallel potential drop and width of the inverted-V region agree well with observation.

Huang, C. Y.↗