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 451 records · Page 25

On upstream blocking in a viscous diffusive stratified flow

The effect of diffusion of specie upon the flow about a transverse flat plate moving horizontally in a viscous stratified medium is considered. Asymptotic expansions are used to define a parameter regime where a viscous-diffusive-buoyancy balance is dominant. The solution, expressed in terms of an inverse Fourier transform, is numerically integrated. The results show that, as in the non-diffusive problem, a region of closed streamlines exists ahead of the body. However, unlike the case where diffusion is neglected, the density field within this recirculating region is uniquely determined and found to be statically stable. It is also found that varying the relative amount of diffusion affects not only the density distribution, but the velocity profile as well, indicating a strong coupling between the vorticity and specie equation.

Koop, C. G.↗

Prediction of the diffuse far infrared flux from the galactic plane

A basic model and simple numerical relations useful for future far infrared studies of the galaxy are presented. Making use of recent CO and other galactic surveys, the diffuse far infrared flux distribution from the galactic plane as a function of galactic longitude alternate theta for 4 deg or = alternate theta or = 90 and the far infrared emissivity as a function of galactocentric distance is predicted. Future measurements of the galactic far infrared flux would yield valuable information on the physical properties and distribution of dust and molecular clouds in the galaxy, particularly the inner region.

Fazio, G. G.↗

Nonlinearity Effects of Lateral Density Diffusion Coefficient on Gain-Guided VCSEL Performance

Electron and hole diffusions in the plane of semiconductor quantum wells play an important part in the static and dynamic operations of semiconductor lasers. In this paper, we apply a hydrodynamic model developed from the semiconductor Bloch equations to numerically study the effects of nonlinearity in the diffusion coefficient on single mode operation and direct modulation of a gain-guided InGaAs/GaAs multiple quantum well laser, operating not too far from threshold. We found that a small diffusion coefficient is advantageous for lowering the threshold current and increasing the modulation bandwidth. Most importantly, the effects of nonlinearity in the coefficient can be approximately reproduced by replacing the coefficient with an effective constant diffusion coefficient, which corresponds roughly to the half height density of the density distribution.

Li, Jian-Zhong↗

Coupled Lattice Boltzmann Modeling Framework for Pore-Scale Fluid Flow and Reactive Transport

In this paper, we propose a modeling framework for pore-scale fluid flow and reactive transport based on a coupled lattice Boltzmann model (LBM). We develop a modeling interface to integrate the LBM modeling code parallel lattice Boltzmann solver and the PHREEQC reaction solver using multiple flow and reaction cell mapping schemes. The major advantage of the proposed workflow is the high modeling flexibility obtained by coupling the geochemical model with the LBM fluid flow model. Consequently, the model is capable of executing one or more complex reactions within desired cells while preserving the high data communication efficiency between the two codes. Meanwhile, the developed mapping mechanism enables the flow, diffusion, and reactions in complex pore-scale geometries. We validate the coupled code in a series of benchmark numerical experiments, including 2D single-phase Poiseuille flow and diffusion, 2D reactive transport with calcite dissolution, as well as surface complexation reactions. The simulation results show good agreement with analytical solutions, experimental data, and multiple other simulation codes. In addition, we design an AI-based optimization workflow and implement it on the surface complexation model to enable increased capacity of the coupled modeling framework. Compared to the manual tuning results proposed in the literature, our workflow demonstrates fast and reliable model optimization results without incorporating pre-existing domain knowledge.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Galactic cosmic ray modulation from 1965-1970.

Numerical solutions of the cosmic-ray equation of transport within the solar cavity and including the effects of diffusion, convection, and energy losses due to adiabatic deceleration, have been used to reproduce the modulation of galactic electrons, protons, and helium nuclei observed during the period from 1965 to 1970. Kinetic energies between 10 and 10,000 MeV/nucleon are considered. Computed and observed spectra are given for the years 1965, 1968, 1969, and 1970 together with the diffusion coefficients. These diffusion coefficients are assumed to be of separable form in rigidity and radial dependence, and are consistent with the available magnetic-field power spectra. The force-field solutions are given for these diffusion coefficients and galactic spectra and are compared with the numerical solutions. It is shown that the energy losses and convection lead to near-earth nuclei spectra at kinetic energies less than or equal to 100 MeV/nucleon in which the differential intensity is proportional to the kinetic energy with little dependence on the form of the galactic spectrum. This dependence is in agreement with the observed spectra of all species of atomic nuclei and it is argued that this provides strong observational evidence for the presence of energy losses in the propagation process, and for the exclusion of low-energy galactic nuclei from near earth.

Urch, I. H.↗

Turbulence and stall in plane diffusers - Computational study

The effect of free-stream turbulence level on stall conditions in plane diffusers is investigated with a finite-volume numerical method which utilizes a k-epsilon turbulence model or a differential Reynolds stress model for closure. Computed results show qualitative agreement with measured data. The ability to eliminate stall and thereby increase pressure recovery in diffusers at large included angles by utilizing high free-stream turbulence is correctly predicted.

Demuren, A. O.↗

A detailed numerical investigation of Burke-Schumann gaseous and spray flames

The classical Burke-Schumann gaseous diffusion flame, and the effect of introducing fuel in the form of liquid sprays are investigated numerically. The time-dependent two-dimensional, axisymmetric conservation equations are solved for the gas phase. A Lagrangian approach is used for the dispersed phase. The chemical reactions are modeled through the one-step global reaction scheme. The numerical results show that the computed flame structure is significantly different from that given by the analytical solution. The computed flame is about 15 percent longer and 5 percent narrower than the classical Burke-Schumann flame. The predicted flame shape is, however, in better agreement with experimental observations than the Burke-Schumann flame. The effects of introducing liquid sprays is found to produce thinner and longer flames than its counterpart gaseous flame. A detailed comparison of the structures of gaseous and spray flames is made, and the differences are discussed.

Mawid, M. A.↗

Annual Research Briefs, 1998

The topics contained in this progress report are direct numerical simulation of turbulent non-premixed combustion with realistic chemistry; LES of non-premixed turbulent reacting flows with conditional source term estimation; measurements of the three-dimensional scalar dissipation rate in gas-phase planar turbulent jets; direct simulation of a jet diffusion flame; on the use of interpolating wavelets in the direct numerical simulation of combustion; on the use of a dynamically adaptive wavelet collocation algorithm in DNS (direct numerical simulation) of non-premixed turbulent combustion; 2D simulations of Hall thrusters; computation of trailing-edge noise at low mach number using LES and acoustic analogy; weakly nonlinear modeling of the early stages of bypass transition; interactions between freestream turbulence and boundary layers; interfaces at the outer boundaries of turbulent motions; largest scales of turbulent wall flows; the instability of streaks in near-wall turbulence; an implementation of the v(sup 2) - f model with application to transonic flows; heat transfer predictions in cavities; a structure-based model with stropholysis effects; modeling a confined swirling coaxial jet; subgrid-scale models based on incremental unknowns for large eddy simulations; subgrid scale modeling taking the numerical error into consideration; towards a near-wall model for LES of a separated diffuser flow; on the feasibility of merging LES with RANS (Reynolds Averaging Numerical simulation) for the near-wall region of attached turbulent flows; large-eddy simulation of a separated boundary layer; numerical study of a channel flow with variable properties; on the construction of high order finite difference schemes on non-uniform meshes with good conservation properties; development of immersed boundary methods for complex geometries; and particle methods for micro and macroscale flow simulations.

Spinks, Debra↗

Effects of confinement on partially premixed flames

Partially premixed combustion is an intermediate regime between the limiting cases of premixed and nonpremixed combustion. Although combustion problems are generally approached from one of these two limiting cases, there are many practical situations where flames cannot be considered as purely premixed or nonpremixed, and thus the partially premixed approach must be used. In partially premixed combustion, mechanisms from the premixed and nonpremixed regimes can coexist, and as a result some interesting new phenomena can arise. One such phenomenon is the flame stabilization in laminar mixing layers by triple flames. One of the first observations of triple flames was made by Phillips (1965), who investigated a triple flame propagating in a methane mixing layer. Kioni et al. (1993) also examined triple flames both experimentally and numerically. There have also been numerous analytical studies on the shape and propagation of triple flames under various assumptions by Dold (1989), Dold et al. (1991), and Hartley and Dold (1991). In terms of modeling, Muller et al. (1994) have combined the flamelet formulations for premixed and nonpremixed combustion in order to treat lifted diffusion flames. One common feature in the analytical and numerical studies mentioned above is the assumption of zero heat release, which is necessary to make the problem tractable. The effect of heat release on triple flames was investigated by Ruetsch et al. (1995), where for the unconfined case, flame speeds larger than their premixed counterparts were found. One of the most important practical situations in which these conditions arise is in lifted turbulent jet diffusion flames. At a critical velocity the burning zone of a fuel jet lifts off from the nozzle, moves to increasing distances as the jet velocity increases, and finally blows off. The mechanisms that control these phenomena, i.e. that determine the stability of these flames, are still not understood. In addition to regions where diffusion flame stabilization takes place, partially premixed conditions also exist during the ignition process in nonpremixed systems. Numerical simulations by Reveillon et al. (1994) of the ignition process in a weakly stirred mixture of fuel and oxidizer show that triple flames propagate along lines of stoichiometric mixture fraction throughout the fluid. In addition, Peters (1994) notes that NO(x) emissions are likely to be large in such transient cases, and therefore an understanding of triple flames can provide information concerning pollutant formation. This study extends the work previously done and examines the effects of lateral confinement on partially premixed flames. Once again, we study both the flame structure and propagation.

Ruetsch, G. R.↗

Application of a finite difference technique to thermal wave propagation

A finite difference formulation is presented for thermal wave propagation resulting from periodic heat sources. The numerical technique can handle complex problems that might result from variable thermal diffusivity, such as heat flow in the earth with ice and snow layers. In the numerical analysis, the continuous temperature field is represented by a series of grid points at which the temperature is separated into real and imaginary terms. Computer routines previously developed for acoustic wave propagation are utilized in the solution for the temperatures. The calculation procedure is illustrated for the case of thermal wave propagation in a uniform property semi-infinite medium.

Baumeister, K. J.↗

Application of a finite difference technique to thermal wave propagation

A finite difference formulation is presented for thermal wave propagation resulting from periodic heat sources. The numerical technique can handle complex problems that might result from variable thermal diffusivity, such as heat flow in the earth with ice and snow layers. In the numerical analysis, the continuous temperature field is represented by a series of grid points at which the temperature is separated into real and imaginary terms. Next, computer routines previously developed for acoustic wave propagation are utilized in the solution for the temperatures. The calculation procedure is illustrated for the case of thermal wave propagation in a uniform property semi-infinite medium.

Baumeister, K. J.↗

Plasma Diffusion in Self-Consistent Fluctuations

The problem of particle diffusion in position space, as a consequence ofeleclromagnetic fluctuations is addressed. Numerical results obtained with a self-consistent hybrid code are presented, and a method to calculate diffusion coefficient in the direction perpendicular to the mean magnetic field is proposed. The diffusion is estimated for two different types of fluctuations. The first type (resuiting from an agyrotropic in itiai setting)is stationary, wide band white noise, and associated to Gaussian probability distribution function for the magnetic fluctuations. The second type (result ing from a Kelvin-Helmholtz instability) is non-stationary, with a power-law spectrum, and a non-Gaussian probabi lity distribution function. The results of the study allow revisiting the question of loading particles of solar wind origin in the Earth magnetosphere.

Smets, R.↗

Diffusion quantum Monte Carlo approach to the polaritonic ground state

Making and using polaritonic states (i.e., hybrid electron-photon states) for chemical applications has recently become one of the most prominent and active fields that connects the communities of chemistry and quantum optics. Modeling of such polaritonic phenomena using ab initio approaches calls for new methodologies, leading to the reinvention of many commonly used electronic structure methods, such as Hartree-Fock, density functional, and coupled cluster theories. Here in this work, we explore the formally exact diffusion quantum Monte Carlo approach to obtain numerical solutions to the polaritonic ground state during the dissociation of the H 2 molecular system. Here, we examine various electron-nuclear-photon properties throughout the dissociation, such as changes to the minimum of the cavity Born-Oppenheimer surface, the localization of the electronic wave function, and the average mode occupation. Finally, we directly compare our results to that obtained with state-of-the-art, yet approximate, polaritonic coupled cluster approaches.

74 ATOMIC AND MOLECULAR PHYSICS↗

Self-Consistent Relativistic Electron Scattering using the Sherlock Scattering Model for X-ray Diagnostics

We present on a new, self-consistent, arbitrary-temperature Romberg integration scheme for modeling electron scattering in materials in a LANL Lagrangian Shock Hydro (LSH) code. Electron beam-target interactions are fundamental to a wide range of scientific and technological applications. When high-energy electron beams hit their target, they may scatter, deposit energy, or ionize the source. These processes govern the behavior and outcomes in nanotechnology manufacturing, electron microscopy, and modern X-ray diagnostics. Simulating these interactions is essential for interpreting experimental results, predicting material responses, and designing efficient tools and experiments. At Los Alamos, this is done using a LSH code, which is a multi-dimension, multi-material, massively parallel, multi-physics code used to simulate applications from asteroid impacts to electron beam interactions. By effectively and efficiently modeling the way that electrons scatter from the beam we can bolster these simulations and more accurately predict experimental outcomes. The model currently implemented in the LSH of interest is based on work by Papp and does not self-consistently preserve momentum in the slightly relativistic regime; here we adopt a model proposed by Braams and Karney and implement a Romberg integration scheme to compute the diffusion tensor. In this paper we will provide background on the Braams-Karney diffusion tensor as well as the Romberg integration scheme we employed to numerically solve for it. We will show that our integration scheme is accurate in solving for the set of scalar potentials used to re-express the diffusion tensor in differential form, and in solving for the diffusion coefficients in the larger LSH code. By using this diffusion tensor rather than the existing Papp one, and numerically integrating it with a Romberg method, we produce much more accurate, self-consistent results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effects of free stream turbulence intensity and integral length scale on heat transfer from a circular cylinder in crossflow

A model has been proposed for the momentum eddy diffusivity induced by free stream turbulence intensity and integral length scale. The eddy diffusivity model is applied to a cylinder situated in a uniform crossflow in the presence of free stream turbulence. Numerical solution of the governing momentum and energy equations with the proposed eddy diffusivity model yielded results for the skin friction coefficient and the Nusselt number. The numerical predictions of the present work are compared with experimental data and the agreement between the two is seen to be very good.

Gorla, R. S. R.↗

Solar-Tracker Diffuse-Response Algorithm for Balancing Energy Gain and Mechanical Wear in Photovoltaic Systems

Single-axis solar tracking maximizes photovoltaic energy production under clear-sky conditions; however, its effectiveness decreases under cloudy and overcast skies, where diffuse irradiance dominates and the optimal module orientation changes. Conventional tracking algorithms either neglect sky conditions or rely on simplified diffuse-response strategies that may trigger frequent tracker repositioning under variable cloud cover, leading to increased mechanical wear with marginal energy gains. This work proposes an enhanced diffuse-response tracking algorithm that explicitly accounts for both the intensity and temporal persistence of cloudiness. By requiring overcast conditions to persist for a minimum duration before reorienting the tracker to a diffuse-stow position, the proposed approach reduces unnecessary movements while preserving the benefits of diffuse-response operation. The algorithm is evaluated through numerical simulations based on historical meteorological data and validated using field measurements on monofacial and bifacial photovoltaic strings. The results show that the proposed strategy reduces excess tracker movement from 114% to 0.16% while maintaining nearly the same energy yield. Compared to a conventional diffuse-response algorithm, the associated energy reduction is minimal (≈0.17%) relative to the ≈0.37% yield gain observed at the studied location. These findings demonstrate that incorporating cloudiness duration enables a practical compromise between energy performance and tracker durability, particularly for monofacial photovoltaic systems.

14 SOLAR ENERGY↗

A numerical model of gravity wave breaking and stress in the mesosphere

The goal of the study is to calculate numerically the deceleration and heating caused by breaking gravity waves. The effect of the radiative dissipation of the wave is included as vertical-wavelength-dependent Newtonian cooling. The parameterization for zonal deceleration is extended by breaking gravity waves (Lindzen, 1981) to include the turbulent diffusion of heat and momentum. After describing the numerical model, the numerical results are presented and compared with the parameterizations in a noninteractive model of the mean zonal wind. Attention is then given to the transport of constituents by gravity waves and the attendant turbulent zone. It is noted that if gravity wave breaking were not an intermittent process, gravity wave stresses would produce an adiabatic mesosphere with a zonal mean velocity close to the phase speed of the breaking wave.

Schoeberl, M. R.↗

Analytic approximations in the study of the solar modulation of electrons

Numerical solutions to the transport equation of galactic cosmic rays in the interplanetary medium have been used to investigate the applicability of commonly used approximate analytic solutions for electrons (10 MeV to 10 GeV). We find that for a given cosmic-ray diffusion coefficient the force-field approximation is in reasonable agreement with the numerical solution at energies above 200 MeV, but deviates significantly at lower energies, depending on the shape of the interstellar electron spectrum. The diffusion-convection approximation agrees generally with the numerical solution within a factor of 2 over the entire energy range.

Cummings, A. C.↗