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 991 records · Page 55

Efficient data-driven regression for reduced-order modeling of spatial pattern formation

We present an efficient data-driven regression approach for constructing reduced-order models (ROMs) of reaction-diffusion systems exhibiting pattern formation. The ROMs are learned non-intrusively from available training data of physically accurate numerical simulations. The method can be applied to general nonlinear systems through the use of polynomial model form, while not requiring knowledge of the underlying physical model, governing equations, or numerical solvers. The process of learning ROMs is posed as a low-cost least-squares problem in a reduced-order subspace identified via Proper Orthogonal Decomposition (POD). Numerical experiments on classical pattern-forming systems–including the Schnakenberg and Mimura–Tsujikawa models–demonstrate that higher-order surrogate models significantly improve prediction accuracy while maintaining low computational cost. The proposed method provides a flexible, non-intrusive model reduction framework, well suited for the analysis of complex spatio-temporal pattern formation phenomena.

Data-driven modeling↗

Universality aspects of quantum corrections to transverse momentum broadening in QCD media

We study non-linear quantum corrections to transverse momentum broadening (TMB) of a fast parton propagating in dense QCD matter in the leading logarithmic approximation. These non-local corrections yield an anomalous super-diffusive behavior characterized by a heavy tailed distribution which is associated with Lévy random walks. Using a formal analogy with the physics of traveling waves, we show that at late times the transverse momentum distribution tends to a universal scaling regime. We derive analytic solutions in terms of an asymptotic expansion around the scaling limit for both fixed and running coupling. We note that our analytic approach yields a good agreement with the exact numerical solutions down to realistic values of medium length. Finally, we discuss the interplay between system size and energy dependence of the diffusion coefficient $\hat{q}$ and its connection with the gluon distribution function that is manifest at large transverse momentum transfer.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Turbulent solutions of the equations of fluid motion

Some turbulent solutions of the unaveraged Navier-Stokes equations (equations of fluid motion) are reviewed. Those equations are solved numerically in order to study the nonlinear physics of incompressible turbulent flow. Initial three-dimensional cosine velocity fluctuations and periodic boundary conditions are used in most of the work considered. The three components of the mean-square velocity fluctuations are initially equal for the conditions chosen. The resulting solutions show characteristics of turbulence such as the linear and nonlinear excitation of small-scale fluctuations. For the stronger fluctuations, the initially nonrandom flow develops into an apparently random turbulence. Thus randomness or turbulence can arise as a consequence of the structure of the Navier-Stokes equations. The cases considered include turbulence which is statistically homogeneous or inhomogeneous and isotropic or anisotropic. A mean shear is present in some cases. A statistically steady-state turbulence is obtained by using a spatially periodic body force. Various turbulence processes, including the transfer of energy between eddy sizes and between directional components, and the production, dissipation, and spatial diffusion of turbulence, are considered. It is concluded that the physical processes occurring in turbulence can be profitably studied numerically.

Deissler, R. G.↗

Direct numerical simulation and Reynolds-averaged Navier-Stokes modeling of the sudden viscous dissipation for multicomponent turbulence

Simulations of a turbulent multicomponent fluid mixture undergoing isotropic deformations are carried out to investigate the sudden viscous dissipation. This dissipative mechanism was originally demonstrated using simulations of an incompressible single-component fluid [S. Davidovits and N. J. Fisch, Phys. Rev. Lett. 116, 105004 (2016)]. By accounting for the convective and diffusive transfer of various species, the current work aims to increase the physical fidelity of previous simulations and their relevance to inertial confinement fusion applications. Direct numerical simulations of the compressed fluid show that the sudden viscous dissipation of turbulent kinetic energy is unchanged from the single-component scenario. More importantly, the simulations demonstrate that the mass fraction variance and covariance for the various species also exhibit a sudden viscous decay. Reynolds-averaged Navier-Stokes simulations were carried out using the $k - l$ model to assess its ability to reproduce the sudden viscous dissipation. Results show that the standard $k - l$ formulation does not capture the sudden decay of turbulent kinetic energy, mass-fraction variance, and mass-fraction covariance for simulations with various compression and expansion rates, or different exponents for the power-law model of viscosity. A new formulation of the $k - l$ model that is based on previous improvements to the $k - ε$ family of models is proposed, which leads to consistently good agreement with the direct numerical simulations for all the isotropic deformations under consideration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Coupled thermal–hydrological–mechanical modeling of CO 2 -enhanced coalbed methane recovery

CO 2 -enhanced coalbed methane recovery, also known as CO 2 -ECBM, is a potential win-win approach for enhanced methane production while simultaneously sequestering injected anthropogenic CO 2 to decrease CO 2 emissions into the atmosphere. Here, CO 2 -ECBM is simulated using a coupled thermal–hydrological–mechanical (THM) numerical model that considers multiphase (gas and water) flow and solubility, multicomponent (CO 2 and CH 4 ) diffusion and adsorption, heat transfer and coal deformation. The coupled model is based on the TOUGH-FLAC simulator, which is applied here for the first time to model CO 2 -ECBM. The capacity of the simulator for modeling methane production is verified by a code-to-code comparison with the general-purpose finite-element solver COMSOL. Then, the TOUGH-FLAC simulator is applied in an isothermal simulation to study the variations in permeability evolution during a CO 2 -ECBM operation while considering four different stress-dependent permeability models that have been implemented into the simulator. Finally, the TOUGH-FLAC simulator is applied in non-isothermal simulations to model THM responses during a CO 2 -ECBM operation.Our simulations show that the permeability evolution, mechanical stress, and deformation are all affected by changes in pressure, temperature and adsorption swelling, with adsorption swelling having the largest effect. The calculated stress changes do not induce any mechanical failure in the coal seam, except near the injection well in one case of a very unfavorable stress field.

01 COAL, LIGNITE, AND PEAT↗

Analytical population dynamics underlying harmful algal blooms triggered by prey avoidance

A modified version of the NPZD ecosystem model is used to analytically examine the effects of predation avoidance, a possible mechanism for triggering harmful algal blooms (HAB). To resolve HAB development caused by predation avoidance, an additional phytoplankton functional group is considered, one that has slower nutrient uptake and better predation avoidance characteristics than the non-harmful phytoplankton group used in traditional NPZD models. Because the two phytoplankton groups (one non-harmful and one HAB) compete for only one resource within the same system, steady state (equilibrium) conditions cannot occur without the presence of zooplankton; only the non-harmful phytoplankton group, which defeats the HAB group in the resource competition, can survive in the equilibrium. The presence of sufficient zooplankton effectively acts to replenish the nutrient pool by consuming the non-harmful phytoplankton. When this occurs, two equilibrium states are found: one with both phytoplankton groups coexisting, and one that only includes the HAB group. Further, the condition required for equilibrium is that the total nitrogen within the system should be larger than a threshold determined by model coefficients. The threshold and feasibility of the equilibrium are sensitive to the relative HAB predation avoidance coefficient. If the coefficient is larger than the ratio of net growth rates between the HAB and non-harmful phytoplankton group, the threshold becomes infinite, and an equilibrium is not feasible. The time scale for the system to reach an equilibrium state that includes a HAB group is determined asymptotically. The dependence of a threshold condition as a controlling factor may explain the regime shift of dominant species causing HABs. The ecosystem model is fully implemented into the Regional Ocean Modeling System and applied to an idealized coastal embayment (with depths and geometry taken from San Francisco Bay) to show numerically the dominance of prey avoidance dynamics in a natural shallow water environment that includes advection and diffusion. The analytical results improve strategies for HAB modeling and provide guidance for setting model coefficients necessary to resolve a HAB event.

54 ENVIRONMENTAL SCIENCES↗

Pulsating Poiseuille flow of a cement slurry

In this paper, we investigate the pulsatile flow of a cement slurry in a pipe. The constitutive relation for the viscous stress tensor is based on the power-law model, where the shear viscosity not only depends on the shear rate but also the volume fraction of the cement particles. To solve for the volume fraction field, a convection-diffusion equation is used. The dimensionless form of the governing equations and the boundary conditions are solved numerically at different pulsatile cycles. Further, a parametric study is conducted to study the effect of different dimensionless numbers at different cycles on the velocity and the volume fraction profiles.

36 MATERIALS SCIENCE↗

Unravelling Photocarrier Dynamics beyond the Space Charge Region for Photoelectrochemical Water Splitting

Semiconductor photoelectrodes for photoelectrochemical (PEC) water splitting require efficient carrier generation, separation, and transport at and beyond the space charge region (SCR) formed at the aqueous interface. The trade-off between photon collection and minority carrier delivery governs the photoelectrode design and implies maximum water splitting efficiency at an electrode thickness equivalent to the light absorption depth. Here, using planar ZnO thin films as a model system, we identify the photocarriers beyond the SCR as another significant source to substantially enhance the PEC performance. The high-quality ZnO films synthesized by pulsed laser deposition feature very few deep trap states and support a long photocarrier lifetime. Combined with photoelectrochemical characterization, ultrafast spectroscopy, and numerical calculations, it is revealed that engineering the exciton concentration gradient by film thickness facilitates the inward diffusion of photocarriers from the neighboring illuminated region to the SCR and, therefore, achieves a record high quantum efficiency over 80% at a thickness far beyond its light absorption depth and the SCR width. Furthermore, these results elucidate the important role of the photocarriers beyond SCR for the PEC process and provide new insight into exploring the full potential for efficient photoelectrode materials with large exciton diffusivity.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Structural magnetic glassiness in the spin ice Dy 2 Ti 2 O 7

The origin and nature of glassy dynamics presents one of the central enigmas of condensed-matter physics across a broad range of systems ranging from window glass to spin glasses. The spin-ice compound Dy 2 Ti 2 O 7 , which is perhaps best known as hosting a three-dimensional Coulomb spin liquid with magnetically charged monopole excitations, also falls out of equilibrium at low temperature. How and why it does so remains an open question. Based on an analysis of low-temperature diffuse neutron-scattering experiments employing different cooling protocols alongside recent magnetic noise studies, combined with extensive numerical modeling, we argue that upon cooling, the spins freeze into what may be termed a “structural magnetic glass,” without an a priori need for chemical or structural disorder. Specifically, our model indicates the presence of frustration on two levels, first producing a near-degenerate constrained manifold inside which phase ordering kinetics is in turn frustrated. A remarkable feature is that monopoles act as sole annealers of the spin network and their pathways and history encode the development of glass dynamics, allowing the glass formation to be visualized. Our results suggest that spin ice Dy 2 Ti 2 O 7 provides one prototype of magnetic glass formation specifically and a setting for the study of kinetically constrained systems more generally.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Characterizing crystalline defects in single nanoparticles from angular correlations of single-shot diffracted X-rays

Characterizing and controlling the uniformity of nanoparticles is crucial for their application in science and technology because crystalline defects in the nanoparticles strongly affect their unique properties. Recently, ultra-short and ultra-bright X-ray pulses provided by X-ray free-electron lasers (XFELs) opened up the possibility of structure determination of nanometre-scale matter with Å spatial resolution. However, it is often difficult to reconstruct the 3D structural information from single-shot X-ray diffraction patterns owing to the random orientation of the particles. This report proposes an analysis approach for characterizing defects in nanoparticles using wide-angle X-ray scattering (WAXS) data from free-flying single nanoparticles. The analysis method is based on the concept of correlated X-ray scattering, in which correlations of scattered X-ray are used to recover detailed structural information. WAXS experiments of xenon nanoparticles, or clusters, were conducted at an XFEL facility in Japan by using the SPring-8 Ångstrom compact free-electron laser (SACLA). Bragg spots in the recorded single-shot X-ray diffraction patterns showed clear angular correlations, which offered significant structural information on the nanoparticles. The experimental angular correlations were reproduced by numerical simulation in which kinematical theory of diffraction was combined with geometric calculations. We also explain the diffuse scattering intensity as being due to the stacking faults in the xenon clusters.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Modeling of non-equilibrium partitioning in direct metal deposited copper–iron concentrated alloys

The direct metal deposition (DMD) additive manufacturing process produces high cooling rates within a small melt pool and can lead to high amounts of solute trapping. These high cooling rates limit diffusion and lead to the formation of non-equilibrium phases. In this work, we utilize a numerical model to calculate the degree of solute trapping, defined as non-equilibrium partitioning. A theoretical case with overall composition fixed to 50Cu-50Fe at.%. was performed to observe the influence of increasing solidification rates. We then simulate DMD of equimolar Cu-Fe powder printed on mild steel substrate and the calculated non-equilibrium phase compositions were consistent with experimental observations reported earlier on this alloy composition. For a single deposited track, cooling rates are high enough to yield significant solute trapping. Here, the degree of solute trapping is highest near the free surface and has a gradient that correlates with the cooling rate gradient.

36 MATERIALS SCIENCE↗

Hot oxygen atoms in the upper atmosphere of Venus

Using data obtained by an ultraviolet spectrometer on the Pioneer-Venus Orbiter, energy and altitude distribution of nonthermal oxygen atoms in the Venus atmosphere for typical daytime and nighttime conditions are calculated using two different numerical methods. Agreement is found between results obtained by the two-stream transport method and the diffusion equation technique. Dissociative recombination of molecular oxygen ions and charge exchange of atomic oxygen ions with the neutral hydrogen and oxygen gas are the source terms, and calculations show that the dissociative recombination term is the dominant one. It is shown that measured hot atom concentrations exceed estimates of the thermal hydrogen atom concentrations derived from Brinton et al. (1980) over most of the daytime exosphere, and the calculated concentrations do so to an even greater extent. It is concluded that hot oxygen atoms are an important component of the dayside exosphere of Venus, their most important effect being the absorption of solar wind.

Nagy, A. F.↗

Numerical analysis of SSME preburner injector atomization and combustion processes

The coaxial spray injection and combustion flowfields of a Space Shuttle Main Engine preburner injector element have been analyzed using a three-phase numerical code. The processes of atomization, evaporation, secondary droplet breakup, and multispecies chemistry, as well as turbulent diffusion, are included. The model produced realistic pictures of the complex internal flowfield, including liquid jet length, spray shape, flame-zone size and characteristics, and predicted temperatures that seem to be in agreement with test data envelopes. It predicted an external group combustion type of flame. Salient combustion and mixing features are discussed and sources of uncertainty are pointed out for future studies.

Liang, P. Y.↗

Navier-Stokes analysis of cold scramjet-afterbody flows

The progress of two efforts in coding solutions of Navier-Stokes equations is summarized. The first effort concerns a 3-D space marching parabolized Navier-Stokes (PNS) code being modified to compute the supersonic mixing flow through an internal/external expansion nozzle with multicomponent gases. The 3-D PNS equations, coupled with a set of species continuity equations, are solved using an implicit finite difference scheme. The completed work is summarized and includes code modifications for four chemical species, computing the flow upstream of the upper cowl for a theoretical air mixture, developing an initial plane solution for the inner nozzle region, and computing the flow inside the nozzle for both a N2/O2 mixture and a Freon-12/Ar mixture, and plotting density-pressure contours for the inner nozzle region. The second effort concerns a full Navier-Stokes code. The species continuity equations account for the diffusion of multiple gases. This 3-D explicit afterbody code has the ability to use high order numerical integration schemes such as the 4th order MacCormack, and the Gottlieb-MacCormack schemes. Changes to the work are listed and include, but are not limited to: (1) internal/external flow capability; (2) new treatments of the cowl wall boundary conditions and relaxed computations around the cowl region and cowl tip; (3) the entering of the thermodynamic and transport properties of Freon-12, Ar, O, and N; (4) modification to the Baldwin-Lomax turbulence model to account for turbulent eddies generated by cowl walls inside and external to the nozzle; and (5) adopting a relaxation formula to account for the turbulence in the mixing shear layer.

Baysal, Oktay↗

Low-energy ion acceleration at quasi-perpendicular shocks: Transverse diffusion

The problem of ion injection and acceleration at quasi perpendicular shocks has been the subject of some debate over the past two decades. It is widely known that these shocks efficiently accelerate particles that are well in the high-energy tail of the distribution. However, the issue of injection, or the acceleration of low-energy ions, has yet to reach a consensus. The fundamental issue is whether there is enough diffusion normal to the magnetic field for the particles to remain near the shock. Since transverse diffusion is a physical process that is not well understood in space plasmas, this is an important, and difficult issue to address. In this report, we will investigate the ion injection problem by performing test particle orbit integrations using synthesized turbulent fields. These fields are fully three-dimensional so that transverse diffusion is possible (cross-field diffusion is not possible in geometries where the electromagnetic fields are less than three dimensional). The synthesized fields are produced by superimposing a three-dimensional wave field on a background field. For completeness, we will compare the results from this model with the more well-established theories, such as the diffusive approximation and scatter-free shock drift acceleration. We will also compare these results with other numerical simulation techniques such as the well known hybrid simulation, and other test-particle calculations in which the shock fields are specified to have less than three dimensions. We will also discuss some recent relevant observations and how these compare with our results.

Giacalone, J.↗

The Effect of Velocity on the Extinction Behavior of a Diffusion Flame during Transient Depressurization

Current fire suppression plans for the International Space Station include the use of venting (depressurization) as a method for extinguishing a fire. Until recently this process had only been examined as part of a material flammability experiment performed on Skylab in the early 1970's. Due to the low initial pressure (0.35 Atm) and high oxygen concentration (65%), the Skylab experimental results are not applicable for understanding the effects of venting on a fire in a space station environment (21%O2, 1 Atm). Recent research examined the extinction behavior of a diffusion flame over a polymethyl methacrylate (PMMA) cylinder during a transient depressurization in low-gravity. The numerical model was used to examine extinction limits as a function of depressurization rate, forced flow velocity, and initial solid phase temperature. The experimental and numerically predicted extinction data indicated that as the solid phase temperature increased the pressure required to extinguish the flame decreased. The numerical model was also used to examine conditions not obtainable in the low-gravity experiments. From these simulations, a series of extinction boundaries were generated that showed a region of increased flammability existed at a forced flow of 10 cm/s. Analysis of these extinction boundaries indicated that they were quasi-steady in nature, and that the final extinction conditions were independent of the transient process. The velocity range in the previous study was limited and thus the results did not examine the effects of velocities less than 1 cm/s or greater than 20 cm/s. This study utilized low-gravity experiments performed on NASA's Reduced-gravity Research Aircraft Laboratory and numerical simulations to examine conditions applicable to the Space Station environment. This paper extends the analysis of the previous study to a comprehensive examination of the effect of increased velocity on extinction behavior and extinction limits during a transient depressurization in low-gravity. This is achieved by examining extinction data from buoyant (normal-gravity) and low-buoyant (low-gravity) depressurization. experiments, as well as from numerical predictions of flame behavior during depressurization in a non-buoyant (zero-gravity) environment.

Goldmeer, Jeffrey S.↗

NASA Galactic Cosmic Radiation Environment Model: Badhwar-O'Neill (2014)

The Badhwar‐O'Neill (BON) Galactic Cosmic Ray (GCR) flux model is used by NASA to certify microelectronic systems and in the analysis of radiation health risks for human space flight missions. Of special interest to NASA is the kinetic energy region below 4.0 GeV/n due to the fact that exposure from GCR behind shielding (e.g., inside a space vehicle) is heavily influenced by the GCR particles from this energy domain. The BON model numerically solves the Fokker‐Planck differential equation to account for particle transport in the heliosphere due to diffusion, convection, and adiabatic deceleration under the assumption of a spherically symmetric heliosphere. The model utilizes a GCR measurements database from various particle detectors to determine the boundary conditions. By using an updated GCR database and improved model fit parameters, the new BON model (BON14) is significantly improved over the previous BON models for describing the GCR radiation environment of interest to human space flight.

O'Neill, P. M.↗

Diffused morphotropic phase boundary in relaxor-PbTiO 3 crystals: High piezoelectricity with improved thermal stability

Solid solution ferroelectrics are the most widely used piezoelectric material for numerous electromechanical applications, including sensors, actuators, and transducers. A milestone in ferroelectric research is the discovery of the morphotropic phase boundary that was first reported in Pb(Zr x Ti 1–x )O 3 , which has been extensively solicited to improve the performance of various solid solution ferroelectrics, including those having the highest piezoelectricity known today. However, due to the inherent correlation between phase transition and thermodynamic imbalance, the efforts of building the phase boundary encounter the challenge that high performance materials are generally accompanied by property instability. In this paper, we report a comprehensive study on the crystalline symmetry and polar nano-regions in relaxor-PbTiO 3 crystals by synchrotron measurements and property characterizations. In contrast to the common belief that the morphotropic phase boundary is a narrow composition region, the morphotropic phase boundary in ternary Pb(In 1/2 Nb 1/2 )O 3 -Pb(Mg 1/3 Nb 2/3 )O 3 -PbTiO 3 is demonstrated as an extensive region covering at least 9% of the PbTiO 3 content. Such a diffused morphotropic phase boundary is associated with the intensive interaction of polar nano-regions, leading to high piezoelectricity (>1500 pC/N) with greatly improved thermal stability, where the piezoelectric variation is ~90% over the temperature range of 273–373 K, which is about a factor of 3 lower compared to its binary counterpart Pb(Mg 1/3 Nb 2/3 )O 3 -PbTiO 3 . This work sheds light on the fundamental understanding of nanoscale inhomogeneity and macroscopic symmetry in relaxor-PbTiO 3 , and successfully interlinks the structure and properties in complex solid solutions. The existence of a diffused morphotropic phase boundary is also expected to benefit other ferroic systems beyond ferroelectrics, such as ferromagnetic and multiferroic materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗