Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical diffusion”

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 955 records · Page 53

N-particle irreducible actions for stochastic fluids

We construct one- and two-particle irreducible (1PI and 2PI) effective actions for the stochastic fluid dynamics of a conserved density undergoing diffusive motion. We compute the 1PI action in one-loop order and the 2PI action in two-loop approximation. We derive a set of Schwinger-Dyson equations and regularize the resulting equations using Pauli-Villars fields. We numerically solve the Schwinger-Dyson equations for a non-critical fluid. We find that higher-loop effects summed by the Schwinger-Dyson renormalize the non-linear coupling. We also find indications of a diffuson-cascade, the appearance of n-loop correction with smaller and smaller exponential suppression.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Radiosensitizing the Vasculature of Primary Brainstem Gliomas Fails to Improve Tumor Response to Radiation Therapy

Diffuse intrinsic pontine gliomas (DIPGs) arise in the pons and are the leading cause of death from brain tumors in children. DIPGs are routinely treated with radiation therapy, which temporarily improves neurological symptoms but generally fails to achieve local control. Because numerous clinical trials have not improved survival from DIPG over standard radiation therapy alone, there is a pressing need to evaluate new therapeutic strategies for this devastating disease. Vascular damage caused by radiation therapy can increase the permeability of tumor blood vessels and promote tumor cell death.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Effect of thermal oxidation on helium implanted pure iron

The impact of thermal oxidation on Helium (He) implanted pure iron (Fe) was investigated to experimentally evaluate how thermal oxidation influences the diffusion and distribution of He within the material. In case of the sample with the lowest dose (2 × 10 17 ions/cm 2 ), the thinnest oxide layer was observed compared to the non-implanted pure Fe. It is due to He nano bubbles implanted in the material, which hinder the diffusion of Fe ions necessary for forming the oxide layer. For the medium dose sample (5 × 10 17 ions/cm 2 ), the oxide layer was slightly thicker than that of the lowest dose, however it had more porous structure due to the numerous nano He bubbles. In addition, a large bubble region around 150–200 nm depth was generated. Here, the highest dose sample (1 × 10 18 ions/cm 2 ) was observed that the escape of He was accelerated by the large amount of nano He bubbles from within the metal, forming a porous Fe matrix on the top surface and an oxide layer that is very porous but thicker than those of lower doses.

He implantation↗

Influence of non-solvent chemistry on polybenzimidazole hollow fiber membrane preparation

Polybenzimidazole membrane materials have attractive H 2 /CO 2 separation characteristics and high thermo-chemical stability for elevated temperature synthesis (syn) gas separations. The development of PBI membranes with a thin defect-free selective layer and porous support morphology is vital to achieving industrially attractive separation performance. This work is focused on developing a fundamental understanding of the liquid-liquid demixing-based phase inversion process for asymmetric PBI hollow fiber membrane (HFM) formation. The development of industrially attractive HFMs is a challenging process due to the complex interplay between phase equilibria, phase inversion kinetics, and interfacial mass transfer that exist during the liquid-liquid demixing process. Numerous parameters including the dope, bore, and outer coagulant chemistries and compositions significantly influence the HFM morphologies produced. Here, a systematic study is conducted to investigate the phase inversion process parameters including the roles of the non-solvent solubility and diffusivity parameters with respect to the solvent and PBI on the phase inversion process. Furthermore, the influence of dope, bore and coagulant chemistries and compositions on PBI HFM morphology are investigated. As a result, the fabricated PBI HFMs are evaluated for their ideal H 2 and CO 2 permeance and H 2 /CO 2 selectivity at 250 °C to benchmark their separation performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

State-variable model of porous electrode batteries 1-D electrochemical/thermal behavior

Compared to full-order numerical models, the state-variable model (SVM) achieves multiple orders of reduction in computation time for simulation of electrochemical/thermal transport and reaction processes occurring inside lithium-ion and similar batteries with porous electrode construction. The model efficiently describes Li diffusion dynamics in active material and electrolyte phases, charge transfer kinetics, charge balance, and electron transport processes. The model requires a separate pre-processing step to identify SVM parameters.

Kim, Gi-Heon↗

Direct numerical simulation of single-species and binary-species boundary layers at high pressure

Direct numerical simulations of single-species and binary-species temporal boundary layers at high pressure are performed. The main objective is to investigate flow physics in a binary-species boundary layer at high pressure, with special attention to mass diffusion. The working fluids are nitrogen in the single-species cases, and a mixture of nitrogen and methane in the binary-species cases. An investigation of mean profiles shows that velocity and temperature profiles have steep gradients near the wall, whereas the mass fraction profiles do not have a steep gradient. This result indicates that the similarity law between velocity, temperature and mass fraction does not hold true in the binary-species boundary layer. The comparison of turbulent fluctuations shows that the qualitative characteristics of velocity and temperature are similar to each other. In contrast, profiles of mass fraction fluctuations are largely different from those of velocity and temperature. This result indicates that mass diffusion is not similar to momentum and thermal diffusions in binary-species boundary layers.

Toki, Takahiko↗

Slant, Fan, and Narrow: The Response of Stellar Streams to a Tilting Galactic Disk

Stellar streams are sensitive tracers of the gravitational potential, which is typically assumed to be static in the inner Galaxy. However, massive mergers like Gaia-Sausage-Enceladus can impart torques on the stellar disk of the Milky Way that result in the disk tilting at rates of up to 10°–20° Gyr –1 . Here, we demonstrate the effects of disk tilting on the morphology and kinematics of stellar streams. Through a series of numerical experiments, we find that streams with nearby apocenters (r apo ≲ 20 kpc) are sensitive to disk tilting, with the primary effect being changes to the stream's on-sky track and width. Interestingly, disk tilting can produce both more diffuse streams and more narrow streams, depending on the orbital inclination of the progenitor and the direction in which the disk is tilting. Our model of Pal 5's tidal tails for a tilting rate of 15° Gyr –1 is in excellent agreement with the observed stream's track and width, and reproduces the extreme narrowing of the trailing tail. We also find that failure to account for a tilting disk can bias constraints on shape parameters of the Milky Way's local dark matter distribution at the level of 5%–10%, with the direction of the bias changing for different streams. Disk tilting could therefore explain discrepancies in the Milky Way's dark matter halo shape inferred using different streams.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Investigation of Gel Polymer Electrolyte Failure Mechanism in Lithium-Sulfur Batteries

Lithium sulfur (Li-S) batteries are a promising candidate for next-generation energy storage due to their high theoretical energy density. However, well-known challenges such as polysulfide crossover, large volume changes during cycling, lithium dendrite growth, and passivation of reaction sites limit the capacity and capacity retention of Li-S batteries. One approach to address these issues is to use gel polymer electrolyte (GPE) in place of the typical liquid electrolyte. It is believed that GPE can reduce the dissolution and diffusion of the intermediate polysulfides and enhance the capacity retention in Li-S batteries, but additional research is needed. In particular, the difference in failure mechanisms between GPE and liquid electrolyte cells is not well understood, which is important knowledge for the rational design of high performance GPE cells. We have investigated the effects of gel polymer electrolyte on lithium-sulfur coin cell failure mechanisms. In particular, we used a discharge test with pauses to probe how the choice of electrolyte affected the amount of capacity recovered due to the pause. After discharging normally until failure, the cell was paused for some time, after which discharge was restarted until the second and final failure. Gel polymer electrolyte was found to have a higher percentage capacity recovery after pausing than liquid electrolyte, and that recovery improved at higher C-rates. Large percentage capacities recoveries suggest that mass transport limitations are leading to cell failure rather than cathode passivation as is frequently proposed. To help explain this behavior, we use a 1-dimension + time numerical continuum model for the Li-S full cell that includes the pausing step. The model includes the electrochemical reactions and precipitation of sulfur and intermediate polysulfide species, transport of species through the battery, and the effect of mass transport limitations on cell failure. The simulations show that variations in the bulk polysulfide diffusion coefficients can explain many of the differences between GPE and liquid electrolyte experimental results, and further confirm that mass transport limitations are a critical obstacle in improving the performance of GPE Li-S batteries.

25 ENERGY STORAGE↗

Radial oscillations of viscous neutron stars: Zero diffusion case

In this article, the spectrum of radial oscillations of neutron stars is systematically studied within two frameworks of viscous relativistic hydrodynamics: the relativistic Navier-Stokes and Israel-Stewart theories. A correspondence is established between the discrete stellar eigenmodes and the continuous dispersion relation of perturbations around a homogeneous fluid, providing a basis for interpreting our numerical results. We analyze the Newtonian limit and assess the impact of relativistic corrections, such as the gravitational redshifting of microscopic relaxation timescales. We show that bulk viscosity can significantly affect the behavior of both hydrodynamic and nonhydrodynamic fundamental modes, and that, depending on the magnitude of the viscous effects, it is the nonhydrodynamic mode that becomes unstable beyond the turning point in a sequence of equilibrium configurations. These results provide a useful step toward systematic studies of neutron star quasinormal modes in the presence of viscosity.

79 ASTRONOMY AND ASTROPHYSICS↗

1-D coupled surface flow and transport equations revisited via the physics-informed neural network approach

The de Saint-Venant equation (SVE) and advection–diffusion equation (ADE) are commonly employed to solve solute transport problems in surface water. In this work, we propose a mesh-free method based on the physics-informed neural network (PINN) to solve the one dimensional (1-D) SVE, ADE, and the coupled SVE and ADE (SVE-ADE) under various initial and boundary conditions. The PINN model extends the architecture of deep neural networks (DNNs) with implementation of loss function, which are additionally subject to constraints imposed by the physical laws of SVE and ADE, along with their initial and boundary conditions. In such a manner, PINNs can be quickly steered to the true solution while obeying the physical laws. The results of PINN model are compared with the analytical and/or numerical solutions under various conditions to investigate its accuracy and efficiency in solving the SVE, ADE, and SVE-ADE. Our results indicate PINN can accurately simulate the shock wave morphology and avoid numerical dissipation in unsteady flow condition. The PINN method outweighed traditional numerical methods in several aspects, including its ability to function with small amounts of data, no grid discretization, and random selection of sampling points, etc. Additionally, the PINN method is also suitable for solving inverse problems with sparse and noisy data. With 1% noise and 2000 initial and boundary condition points (N u ), the errors of the estimated flow rate (v) and diffusion coefficient (D) are 0.003% and 0.105%, respectively, which indicate the accuracy and robustness of the proposed method. Finally, our results indicate the capability and robustness of the proposed PINN methodology for solving multi-physics problems, irrespective of the presence of sparse and noisy data in the training dataset.

54 ENVIRONMENTAL SCIENCES↗

A fractional model for anomalous diffusion with increased variability: Analysis, algorithms and applications to interface problems

Fractional equations have become the model of choice in several applications where heterogeneities at the microstructure result in anomalous diffusive behavior at the macroscale. Here, we introduce a new fractional operator characterized by a doubly-variable fractional order and possibly truncated interactions. Under certain conditions on the model parameters and on the regularity of the fractional order we show that the corresponding Poisson problem is well-posed. Additionally, we introduce a finite element discretization and describe an efficient implementation of the finite-element matrix assembly in the case of piecewise constant fractional order. Through several numerical tests, we illustrate the improved descriptive power of this new operator across media interfaces. Furthermore, we present one-dimensional and two-dimensional h-convergence results that show that the variable-order model has the same convergence behavior as the constant-order model.

97 MATHEMATICS AND COMPUTING↗

Models of primary runaway electron distribution in the runaway vortex regime

Generation of runaway electrons (RE) beams can possibly induce the most deleterious effect of tokamak disruptions. A number of recent numerical calculations have confirmed the formation of a RE bump in their energy distribution by taking into account Synchrontron radiational damping force due to RE’s gyromotions. Here, we present a detailed examination on how the bump location changes at different pitch-angle and the characteristics of the RE pitch-angle distribution. Although REs moving along the magnetic field are preferably accelerated and then populate the phase-space of larger pitch-angle mainly through diffusions, an off-axis peak can still form due to the presence of the vortex structure which causes accumulation of REs at low pitch-angle. A simplified Fokker- Planck model and its semi-analytical solutions based on local expansions around the O point is used to illustrate the characteristics of RE distribution around the O point of the runaway vortex in phase-space. The calculated energy location of the O point together with the local energy and pitch-angle distributions agree with the full numerical solution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Performance prediction of electrohydrodynamic thrusters by the perturbation method

In this paper, we present a novel method for analyzing electrohydrodynamic (EHD) thrusters. The method is based on a perturbation technique applied to a set of drift-diffusion equations, similar to the one introduced in our previous study on estimating breakdown voltage. The thrust-to-current ratio is generalized to represent the performance of EHD thrusters. We have compared the thrust-to-current ratio obtained theoretically with that obtained from the proposed method under atmospheric air conditions, and we have obtained good quantitative agreement. Also, we have conducted a numerical simulation in more complex thruster geometries, such as the dual-stage thruster developed by Masuyama and Barrett [Proc. R. Soc. A 469, 20120623 (2013)]. We quantitatively clarify the fact that if the magnitude of a third electrode voltage is low, the effective gap distance shortens, whereas if the magnitude of the third electrode voltage is sufficiently high, the effective gap distance lengthens.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Current-Assisted Diffusion Bonding of Extruded Ti-22Al-25Nb Alloy by Spark Plasma Sintering: Interfacial Microstructure and Mechanical Properties

Spark plasma sintering (SPS) technology was used to current-assisted bond extruded Ti-22Al-25Nb alloy. The effects of bonding temperature (920-980 °C) and bonding time (10-30 min) on the microstructure evolution and shear strength of this alloy were investigated systematically. The temperature distribution in the specimen during the current-assisted bonding process was also analyzed by numerical simulation. It is noted that the highest temperature was obtained at the bonding interface. As the bonding temperature and bonding time increased, the voids in the interface shrank increasingly until they vanished. A complete metallurgical bonding interface could be produced at 960 °C/20 min/10 MPa, exhibiting the highest shear strength of 269.3 MPa. In addition, the shear strength of the bonded specimen depended on its interfacial microstructure. With increased bonding temperature, the fracture mode transformed from the intergranular fracture at the bonding interface to the cleavage fracture in the substrate.

36 MATERIALS SCIENCE↗

The Role of Interfacial Thermal Resistance in Li-Ion Battery Thermal Management: Preprint

Temperature critically affects the performance, life and safety of lithium-ion batteries. Therefore, it is essential to understand heat generation and dissipation within individual battery cells and battery packs to plan a proper thermal management strategy. One of the key challenges is that interfacial heat transfer of a battery unit is difficult to quantify. The steady-state absolute method and the transient laser-flash-diffusivity method were employed to measure heat conductivities of battery layer stacks and individual battery layer separately. Results show flash diffusivity method gives higher thermal conductivity at both cross-plane and in-plane directions. The difference is primarily caused by interfacial thermal resistance so that it can be estimated by steady-state and transient measurements. To investigate the effects of interfacial thermal transport beyond individual cell level, a multiphysics battery model is used. The model is built upon a multi-scale multi-domain modeling framework for battery packs that accounts for the interplay across multiple physical phenomena. Benefits of a battery module using thermal management materials are quantified through numerical experiments. During a thermal runaway event, it is found interfacial thermal resistance can mitigate thermal runaway in a battery module by significantly reducing heat transfer between cells.

33 ADVANCED PROPULSION SYSTEMS↗

Water vapor oxidation of SiC layer in surrogate TRISO fuel particles

Under accidental conditions for high temperature gas-cooled reactors (HTGR), the SiC layer in tri-structural-isotropic (TRISO) fuel particles can be exposed to water vapor. In this study, oxidation behaviors of surrogate TRISO fuel particles were investigated in a He-20 vol% water vapor mixed atmosphere at temperatures up to 1600 °C. The growth of the crystalline oxide passivation layer with cracks and pores followed a parabolic law with time, where the maximum was 2.3 μm at 1600 °C. The oxide layer and the SiC surface under the oxide became flattened with increasing temperature, as a function of the silica viscosity and the diffusion path of water vapor. Volatilization of the oxide layer was analyzed using a mechanistic model that an inert gas in oxidizing atmospheres could influence the magnitude of volatilization. The fracture load and strength of the oxidized and thinned SiC layer were numerically estimated to decrease from 2.27 to 1.69 N and from 317 to 299 MPa, respectively, with the SiC thickness decrease from 35 to 32 μm. This prediction indicates that the oxidized SiC layer should retain fission products. Additionally, the mechanical integrity of each layer in the TRISO fuel particle after oxidation was evaluated. The results in this work provide important data for the safety analysis of accidental scenarios in HTGRs.

36 MATERIALS SCIENCE↗

Geologic stress modulates fluid mixing at fracture intersections

Fracture intersections are critical links that enable flow and transport in subsurface fracture networks, and their behavior strongly influences fluid mixing in a network. Although all subsurface fractures are subjected to geological stress, we lack a fundamental understanding of how fracture intersection geometry evolves under stress and how these changes influence fluid mixing. Here, we combine 3D printing, 3D X-ray tomographic imaging, and 3D pore-scale numerical simulations to reveal stress-induced changes in intersection geometry and their impact on mixing. Mixing is found to be strongly affected by partial closure of an intersection under stress. As an intersection closes, the void area for fluid flow and diffusion decreases leading to substantial deviations between conventional mixing models and full pore-scale modeling. To address this, we propose a modified mixing model that accounts for intersection deformation, which is essential for accurate modeling of solute transport and mixing through fracture networks.

15 GEOTHERMAL ENERGY↗