Modeling coupled electro-chemo-mechanical phenomena within all-solid-state battery composite cathodes
Not Available
25 ENERGY STORAGE↗
SEARCH · Search NASA
Results for “phase field modeling”
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.
Extending Interfaces in 3D to Achieve Superior Nanoscale Strength in Ti/Nb Nanolaminates
Tuning the atomic-level structure of nanolaminates enables high strength, increased deformability, and the ability to absorb and mitigate damage due to varied crystalline defects, including dislocations. Here, we present the enhanced strength of Ti/Nb nanolaminates containing thick 3D interfaces (3DIs), relative to their chemically abrupt 2D counterparts. We examine the effects of crystallographic alignment and compositional gradients on mechanical behavior via experiments and phase-field-dislocation dynamics (PFDD) modeling. Mechanical testing reveals that nanolaminates containing thicker 3DIs demonstrate a 28% hardness enhancement compared with sharp-interface nanolaminates. PFDD modeling shows that the critical resolved shear stress (CRSS) increases with the 3DI thickness. Gradual compositional transitions in 3DIs were confirmed via scanning transmission electron microscopy and high-resolution transmission electron microscopy, showing sharp crystallographic transitions and a heightened interface topography. The findings establish a positive function between the 3DI thickness and mechanical robustness for hexagonal-closest-packed-containing composites, emphasizing the role of defect–interface interactions in tailoring the mechanical performance and providing a foundation for future interfacial engineering.
36 MATERIALS SCIENCE↗
Accelerating Phase Field Simulations Through Convolution-Only Deep Learning and Diffusion Models
Explore the source record for details and available documents.
Safta, Cosmin [Sandia National Laboratories (SNL-C↗
Quantum Filtering and Analysis of Multiplicities in Eigenvalue Spectra
Fine-grained spectral properties of quantum Hamiltonians, including both eigenvalues and their multiplicities, provide useful information for characterizing many-body quantum systems as well as for understanding phenomena such as topological order. Extracting such information with small additive error is #BQP-complete in the worst case. In this work, we introduce QFAMES (quantum filtering and analysis of multiplicities in eigenvalue spectra), a quantum algorithm that efficiently identifies clusters of closely spaced dominant eigenvalues and determines their multiplicities under physically motivated assumptions, which allows us to bypass worst-case complexity barriers. QFAMES also enables the estimation of observable expectation values within targeted energy clusters, providing a powerful tool for studying quantum phase transitions and other physical properties. We validate the effectiveness of QFAMES through numerical demonstrations, including its applications to characterizing quantum phases in the transverse-field Ising model and estimating the ground-state degeneracy of a topologically ordered phase in the two-dimensional toric code model. We also generalize QFAMES to the setting of mixed initial states. Our approach offers rigorous theoretical guarantees and significant advantages over existing subspace-based quantum spectral analysis methods, particularly in terms of the sample complexity and the ability to resolve degeneracies.
97 MATHEMATICS AND COMPUTING↗
Applicability of empirical data currently used in predicting solid propellant exhaust plumes
Theoretical and experimental approaches to exhaust plume analysis are compared. A two-phase model is extended to include treatment of reacting gas chemistry, and thermodynamical modeling of the gaseous phase of the flow field is considered. The applicability of empirical data currently available to define particle drag coefficients, heat transfer coefficients, mean particle size, and particle size distributions is investigated. Experimental and analytical comparisons are presented for subscale solid rocket motors operating at three altitudes with attention to pitot total pressure and stagnation point heating rate measurements. The mathematical treatment input requirements are explained. The two-phase flow field solution adequately predicts gasdynamic properties in the inviscid portion of two-phase exhaust plumes. It is found that prediction of exhaust plume gas pressures requires an adequate model of flow field dynamics.
Tevepaugh, J. A.↗
Reports From EXPORTS Modeling and Data-Mining Activities
EXPORTS -EXport Processes in the Ocean from RemoTe Sensing is NASA’s large field campaign focusing on development of a predictive understanding of the export, fate and carbon cycle impacts of the global net primary production. This co-funded program (NSF, private funding1, and international participation) was conceived in 2013, with EXPORTS Science plan published in 2016 (Siegel et al 2016, EXPORTS Writing Team 2015) and its implementation plan finalized in 2016 (EXPORTS Science Definition Team 2016). More details are presented in Siegel et al (2021). EXPORTS campaign is structured as a multiyear effort (Figure 1). It started with “Pre-EXPORTS” modeling and data-mining activity followed by a first phase with two major field programs and a second synthesis and modeling phase. The “Pre-EXPORTS” projects, total of 6 of them (Table 1), funded under A.3 Ocean Biology and Biogeochemistry 2015 call, helped to plan the field campaign (Resplandy et al 2019, Rousseaux & Gregg 2017), and supporting further global synthesis with datasets mined from the literature (e.g. Bisson et al (2020), Bisson et al (2018), Kramer and Siegel (2019)), directly responding to objectives outlined in Science Plan (see section 6 in Team (2015)) and Implementation team (see Figure 1 in Team (2016)). This document presents a compilation of the final reports of the Pre-EXPORTS funded projects, in hope of synthesizing the outcomes, and insuring the legacy of this program. Each of these reports contains a list of published papers, and reader should refer to them to see results in details.
Brandi J. McCarty↗
Parameterization of BISON model for fission gas release in U3Si2 nuclear fuel using phase-field simulations
Summary of mesoscale U3Si2 simulations to parametrize BISON, for submission to the CASL closeout technical symposium.
11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
High-Fidelity Generalization Method of Cells for Inelastic Periodic Multiphase Materials
An extension of a recently-developed linear thermoelastic theory for multiphase periodic materials is presented which admits inelastic behavior of the constituent phases. The extended theory is capable of accurately estimating both the effective inelastic response of a periodic multiphase composite and the local stress and strain fields in the individual phases. The model is presently limited to materials characterized by constituent phases that are continuous in one direction, but arbitrarily distributed within the repeating unit cell which characterizes the material's periodic microstructure. The model's analytical framework is based on the homogenization technique for periodic media, but the method of solution for the local displacement and stress fields borrows concepts previously employed by the authors in constructing the higher-order theory for functionally graded materials, in contrast with the standard finite-element solution method typically used in conjunction with the homogenization technique. The present approach produces a closed-form macroscopic constitutive equation for a periodic multiphase material valid for both uniaxial and multiaxial loading. The model's predictive accuracy in generating both the effective inelastic stress-strain response and the local stress said inelastic strain fields is demonstrated by comparison with the results of an analytical inelastic solution for the axisymmetric and axial shear response of a unidirectional composite based on the concentric cylinder model, and with finite-element results for transverse loading.
Aboudi, Jacob↗
An AeroCom Assessment of Black Carbon in Arctic Snow and Sea Ice
Though many global aerosols models prognose surface deposition, only a few models have been used to directly simulate the radiative effect from black carbon (BC) deposition to snow and sea ice. Here, we apply aerosol deposition fields from 25 models contributing to two phases of the Aerosol Comparisons between Observations and Models (AeroCom) project to simulate and evaluate within-snow BC concentrations and radiative effect in the Arctic. We accomplish this by driving the offline land and sea ice components of the Community Earth System Model with different deposition fields and meteorological conditions from 2004 to 2009, during which an extensive field campaign of BC measurements in Arctic snow occurred. We find that models generally underestimate BC concentrations in snow in northern Russia and Norway, while overestimating BC amounts elsewhere in the Arctic. Although simulated BC distributions in snow are poorly correlated with measurements, mean values are reasonable. The multi-model mean (range) bias in BC concentrations, sampled over the same grid cells, snow depths, and months of measurements, are −4.4 (−13.2 to +10.7) ng/g for an earlier phase of AeroCom models (phase I), and +4.1 (−13.0 to +21.4) ng/g for a more recent phase of AeroCom models (phase II), compared to the observational mean of 19.2 ng/g. Factors determining model BC concentrations in Arctic snow include Arctic BC emissions, transport of extra-Arctic aerosols, precipitation, deposition efficiency of aerosols within the Arctic, and meltwater removal of particles in snow. Sensitivity studies show that the model-measurement evaluation is only weakly affected by meltwater scavenging efficiency because most measurements were conducted in non-melting snow. The Arctic (60-90degN) atmospheric residence time for BC in phase II models ranges from 3.7 to 23.2 days, implying large inter-model variation in local BC deposition efficiency. Combined with the fact that most Arctic BC deposition originates from extra-Arctic emissions, these results suggest that aerosol removal processes are a leading source of variation in model performance. The multi-model mean (full range) of Arctic radiative effect from BC in snow is 0.15 (0.07-0.25) W/sq m and 0.18 (0.06-0.28) W/sq m in phase I and phase II models, respectively. After correcting for model biases relative to observed BC concentrations in different regions of the Arctic, we obtain a multi-model mean Arctic radiative effect of 0.17 W/sq m for the combined AeroCom ensembles. Finally, there is a high correlation between modeled BC concentrations sampled over the observational sites and the Arctic as a whole, indicating that the field campaign provided a reasonable sample of the Arctic.
carbon↗
Predicting fracture behavior in single crystal nickel using a coupled crystal plasticity phase field damage approach
Understanding the fracture behavior of single crystal metals is critical for predicting material performance under mechanical loading. Here, in this study, we investigate the fracture characteristics of single crystal nickel tensile bars using a crystal plasticity coupled phase field damage (CP-PFD) model. Experimental tensile tests were conducted on 15 specimens spanning five crystallographic orientations and three thickness variants per orientation. The results revealed two distinct fracture modes: brittle fractures with 45-degree angled surfaces and ductile fractures characterized by significant necking. The CP-PFD model successfully replicated these fracture behaviors, demonstrating strong agreement with experimental observations. The model effectively predicted the strain at which necking and fracture occurred, as well as the orientation-dependent fracture mechanisms. By comparing experimental and simulated fracture surfaces, we establish the CP-PFD model as a robust tool for predicting single crystal behavior and damage evolution. This work provides insight into the microstructural dependence of fracture behavior and establishes a predictive framework for modeling orientation-dependent damage evolution in single-crystal nickel.
Crystal plasticity↗
Coupled Monte Carlo Probability Density Function/ SPRAY/CFD Code Developed for Modeling Gas-Turbine Combustor Flows
The success of any solution methodology for studying gas-turbine combustor flows depends a great deal on how well it can model various complex, rate-controlling processes associated with turbulent transport, mixing, chemical kinetics, evaporation and spreading rates of the spray, convective and radiative heat transfer, and other phenomena. These phenomena often strongly interact with each other at disparate time and length scales. In particular, turbulence plays an important role in determining the rates of mass and heat transfer, chemical reactions, and evaporation in many practical combustion devices. Turbulence manifests its influence in a diffusion flame in several forms depending on how turbulence interacts with various flame scales. These forms range from the so-called wrinkled, or stretched, flamelets regime, to the distributed combustion regime. Conventional turbulence closure models have difficulty in treating highly nonlinear reaction rates. A solution procedure based on the joint composition probability density function (PDF) approach holds the promise of modeling various important combustion phenomena relevant to practical combustion devices such as extinction, blowoff limits, and emissions predictions because it can handle the nonlinear chemical reaction rates without any approximation. In this approach, mean and turbulence gas-phase velocity fields are determined from a standard turbulence model; the joint composition field of species and enthalpy are determined from the solution of a modeled PDF transport equation; and a Lagrangian-based dilute spray model is used for the liquid-phase representation with appropriate consideration of the exchanges of mass, momentum, and energy between the two phases. The PDF transport equation is solved by a Monte Carlo method, and existing state-of-the-art numerical representations are used to solve the mean gasphase velocity and turbulence fields together with the liquid-phase equations. The joint composition PDF approach was extended in our previous work to the study of compressible reacting flows. The application of this method to several supersonic diffusion flames associated with scramjet combustor flow fields provided favorable comparisons with the available experimental data. A further extension of this approach to spray flames, three-dimensional computations, and parallel computing was reported in a recent paper. The recently developed PDF/SPRAY/computational fluid dynamics (CFD) module combines the novelty of the joint composition PDF approach with the ability to run on parallel architectures. This algorithm was implemented on the NASA Lewis Research Center's Cray T3D, a massively parallel computer with an aggregate of 64 processor elements. The calculation procedure was applied to predict the flow properties of both open and confined swirl-stabilized spray flames.
Source record↗
Analysis of heat transfer and AuNPs-mediated photo-thermal inactivation of E. coli at varying laser powers using single-phase CFD modeling
In the wake of the COVID-19 pandemics, the demand for innovative and effective methods of bacterial inactivation has become a critical area of research, providing the impetus for this study. The purpose of this research is to analyze the AuNPs-mediated photothermal inactivation of E. coli. Gold nanoparticles irradiated by laser represent a promising technique for combating bacterial infection that combines high-tech and scientific progress. The intermediate aim of the work was to present the calibration of the model with respect to the gold nanorods experiment. The purpose of this work is to study the effect of initial concentration of E. coli bacteria, the design of the chamber and the laser power on heat transfer and inactivation of E. coli bacteria. Using the CFD simulation, the work combines three main concepts. 1. The conversion of laser light to heat has been described by a combination of three distinctive approximations: a- Discrete particle integration to take into account every nanoparticle within the system, b- Rayleigh-Drude approximation to determine the scattering and extinction coefficients and c- Lambert–Beer–Bourger law to describe the decrease in laser intensity across the AuNPs. 2. The contribution of the presence of E. coli bacteria to the thermal and fluid-dynamic fields in the microdevice was modeled by single-phase approach by determining the effective thermophysical properties of the water-bacteria mixture. 3. An approach based on a temperature threshold attained at which bacteria will be inactivated, has been used to predict bacterial response to temperature increases. The comparison of the thermal fields and temporal temperature changes obtained by the CFD simulation with those obtained experimentally confirms the accuracy of the light-heat conversion model derived from the aforementioned approximations. The results show a linear relationship between maximum temperature and variation in laser power over the range studied, which is in line with previous experimental results. It was also found that the temperature inside the microchamber can exceed 55 °C only when a laser power higher than 0.8 W is used, so bacterial inactivation begins. The experimental data allows to determinate the concentration of nanoparticles. This parameter is introduced into the mathematical model obtaining the same number of AuNPs. However, this assumption introduces a certain simplification, as in the mathematical model the distribution of nanoparticles is uniform. This work is directly connected to the use of gold nanoparticles for energy conversion, as well as the field of bacterial inactivation in microfluidic systems such as lab-on-a-chip. Presented mathematical and numerical models can be extended to the entire spectrum of wavelengths with particular use of white light in the inactivation of bacteria. This work represents a significant advancement in the field, as to the best of the authors’ knowledge, it is the first to employ a single-phase computational fluid dynamics (CFD) approach specifically combined with the thermal inactivation of bacteria. Moreover, this research pioneers the use of a numerical simulation to analyze the temperature threshold of photothermal inactivation of E. coli mediated by gold nanorods (AuNRs). The integration of these methodologies offers a new perspective on optimizing bacterial inactivation techniques, making this study a valuable contribution to both computational modeling and biomedical applications.
36 MATERIALS SCIENCE↗
Flux-Closure Domain Structures in Ferroelectric K 0.5 Na 0.5 NbO 3 Thin Films
Topological domain structures in ferroelectric materials have garnered increasing attention due to their intriguing physical properties and promising applications. While most existing topological structures in ferroelectric perovskite oxides originate from tetragonal or rhombohedral bulk phases, much less is understood about their counterparts in orthorhombic ferroelectrics. Here, in this work, we employ ferroelectric K 0.5 Na 0.5 NbO 3 (KNN) thin films as a model system and leverage phase-field simulations to theoretically predict the static structures and dynamic behaviors of three types of flux-closure domain configurations: in-plane (Type-I), out-of-plane (Type-II), and superdomain (Type-III) flux-closure structures. We systematically investigate the effects of finite size, misfit strains, and electrical boundary conditions on the formation and switching of these topological structures. For the Type-I structure, size reduction or small misfit strain facilitates a transition of the flux-closure pattern to polar vortices. Type-II structures emerge under open-circuit electrical boundary conditions of the film, forming at the junctions of specific domain walls with the film surface or the film–substrate interface. The formation mechanisms of these two flux-closure structures are rationalized from an energy perspective. We demonstrated switching capabilities of Type-II and Type-III structures by obtaining polarization–electric field hysteresis loops. Our simulations also reveal a reversible electric-field-induced transition between the orthorhombic and rhombohedral ferroelectric phases with a checkerboard domain pattern, during which the integrity of the flux-closure structure is preserved. These findings provide theoretical insights and practical guidance for identifying and manipulating topological structures in low-symmetry ferroelectrics, paving the way for developing energy-efficient microelectronic devices based on topological structures.
P-E loop↗
The Dynamics of Vapor Bubbles in Acoustic Pressure Fields
In spite of a superficial similarity with gas bubbles, the intimate coupling between dynamical and thermal processes confers to oscillating vapor bubbles some unique characteristics. This paper examines numerically the validity of some asymptotic-theory predictions such as the existence of two resonant radii and a limit size for a given sound amplitude and frequency. It is found that a small vapor bubble in a sound field of sufficient amplitude grows quickly through resonance and continues to grow thereafter at a very slow rate, seemingly indefinitely. Resonance phenomena therefore play a role for a few cycles at most, and reaching a limit size-if one exists at all-is found to require far more than several tens of thousands of cycles. It is also found that some small bubbles may grow or collapse depending on the phase of the sound field. The model accounts in detail for the thermo-fluid-mechanic processes in the vapor. In the second part of the paper, an approximate formulation valid for bubbles small with respect to the thermal penetration length in the vapor is derived and its accuracy examined, The present findings have implications for acoustically enhanced boiling heat transfer and other special applications such as boiling in microgravity.
Hao, Y.↗
A chain stretch-based gradient-enhanced model for damage and fracture in elastomers
Similar to quasi-brittle materials, it has been recently shown that elastomers can exhibit a macroscopically diffuse damage zone that accompanies the fracture process. In this study, we introduce a stretch-based gradient-enhanced damage (GED) model that allows the fracture to localize and also captures the development of a physically diffuse damage zone. This capability contrasts with the paradigm of the phase field method for fracture, where a sharp crack is numerically approximated in a diffuse manner. Capturing fracture localization and diffuse damage in our approach is achieved by considering nonlocal effects that encompass network topology, heterogeneity, and imperfections. These considerations motivate the use of a statistical damage function dependent upon the nonlocal deformation state. From this model, fracture toughness is realized as an output. While GED models have been classically utilized for damage modeling of structural engineering materials (e.g., concrete), they face challenges when trying to capture the cascade from damage to fracture, often leading to damage zone broadening (de Borst and Verhoosel, 2016). This deficiency contributed to the popularity of the phase-field method over the GED model for elastomers and other quasi-brittle materials. Other groups have proceeded with damage-based GED formulations that prove identical to the phase-field method (Lorentz et al., 2012), but these inherit the aforementioned limitations. To address this issue in a thermodynamically consistent framework, we implement two modeling features (a nonlocal driving force bound and a simple relaxation function) specifically designed to capture the evolution of a physically meaningful damage field and the simultaneous localization of fracture, thereby overcoming a longstanding obstacle in the development of these nonlocal strain- or stretch-based approaches. Here, we discuss several numerical examples to understand the features of the approach at the limit of incompressibility, and compare them to the phase-field method as a benchmark for the macroscopic response and fracture energy predictions.
Elastomers↗
Modeling Radiolysis and Chemical Reactions during Dry Storage of Aluminum-clad Spent Nuclear Fuel
After aluminum-clad spent nuclear fuel (ASNF) is removed from the reactor, it is initially stored in spent fuel pools, which are specially designed water-filled basins that provide temporary cooling to reduce the temperature of the fuel assemblies and provide radiation shielding. ASNF continues to generate heat due to the radioactive decay of elements within the fuel, which persists for many years post-shutdown as the residual radioactive products decay into more stable elements. During the wet storage period, an oxyhydroxide layer composed of boehmite/bayerite forms on the surfaces of the aluminum cladding from exposure to water in the pools. Road-ready packaging for long-term disposition of the ASNF involves dry storage in helium backfilled DOE standard canisters (DSCs). When the ASNF is removed from water storage and dried, most of the water is removed, but some physisorbed and chemisorbed water remains in the oxyhydroxide layers. This residual water can produce hydrogen when exposed to radiation from the ASNF during dry storage. Predicting hydrogen accumulation over time in the DSCs is critical for long-term storage considerations. Previous modeling efforts have developed coupled computational fluid dynamics (CFD)-chemical models to simulate temperature, pressure, and gas phase concentrations within the DSCs. These models use the thermal field predicted by CFD as input to a radiolysis model for the gas phase and the surface oxyhydroxide layer chemistry. Given the long storage period of the DSCs and the impracticality of long-term experiments, a simulation-based approach is necessary to assess chemical evolution within the canisters. This study advances the development of a modeling framework designed to simulate the chemical evolution of spent fuel canisters. Both thermal and radiation-driven reactions are considered, with radiation kinetics quantified using G-values. Sensitivity analysis identifies key parameters influencing species composition. Reaction pathway diagrams offer insight into dominant species formation routes, enabling more effective comparisons between model predictions and experimental observations, particularly regarding the production of hydrogen. Results show that the model predicts significant hydrogen gas production with minimal oxygen generation, primarily due to hydrogen formation via boehmite pathways. These findings underscore the importance of accurately characterizing surface-bound species and radiolysis kinetics. A deeper understanding of these mechanisms is critical for evaluating the long-term safety of nuclear waste storage.
12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗
Primordial black holes and magnetic fields in conformal neutrino mass models
Sufficiently strong and long-lasting first-order phase transitions can produce primordial black holes (PBHs) that contribute substantially to the dark matter abundance of the Universe, and can produce large-scale primordial magnetic fields. Here, we study these mechanisms in a generic class of conformal U(1)' models that also explain active neutrino oscillation data via the type-I seesaw mechanism. We find that phase transitions that occur at seesaw scales between 10 4 GeV and 10 11 GeV produce gravitational wave signals (from the dynamics of the phase transition and from the decay of cosmic string loops) at LISA/ET that can be correlated with microlensing signals of PBHs at the Roman Space Telescope, while scales near 10 11 GeV can be correlated with Hawking evaporation signals at future gamma-ray telescopes. LISA can probe the entire range of PBH masses between 1 × 10 -16 M ⊙ and 8 × 10 -11 M ⊙ if PBHs fully account for the dark matter abundance. For Z' masses between 40 TeV and 10 4 TeV, and 10 TeV right-handed neutrinos, helical magnetic fields can be produced with magnitudes ≳ 0.5 pG and coherence lengths ≳ 0.008 Mpc, above current blazar lower bounds.
cosmological phase transitions↗
Preliminary calibration plan for the Advanced Particles and Field Observatory (APAFO) magnetometer experiment
Prelaunch and postlaunch calibration plans for the APAFO magnetometer experiment are presented. A study of tradeoffs between boom length and spacecraft field is described; the results are summarized. The prelaunch plan includes: calibration of the Vector Fluxgate Magnetometer (VFM), Star Sensors, and Scalar Helium Magnetometer (SHM); optical bench integration; and acquisition of basic spacecraft field data. Postlaunch calibration has two phases. In phase one, SHM data are used to calibrate the VFM, total vector magnetic field data are used to calibrate a physical model of the spacecraft field, and both calibrations are refined by iteration. In phase two, corrected vector data are transformed into geocentric coordinates, previously undetected spacecraft fields are isolated, and initial geomagnetic field models are computed. Provided the SHM is accurate to the required 1.0 nT and can be used to calibrate the VFM to the required 3.0- nT accuracy, the tradeoff study indicates that a 12 m boom and a spacecraft field model uncertainty of 5 percent together allow the 1.0 nT spacecraft field error requirement to be met.
Voorhies, C. V.↗