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 469 records · Page 26

Finite-element analysis of turbulent flow in annular exhaust diffusers of gas turbine engines

A modified version of the Petrov-Galerkin weighted residual method coupled with a biquadratic finite element of the Lagrangian type was used to develop a finite-element model of the turbulent flow field in the annular exhaust diffuser of a gas turbine engine. The swirling flow field was analyzed with emphasis on the diffuser off-design operation. A comparison of the numerical results with experimental data shows that the model is applicable to moderately separating flows of the kind that are typically associated with the off-design performance of diffusing passages in gas turbines.

Baskharone, E. A.↗

Momentum deposition on Wolf-Rayet winds: Nonisotropic diffusion with effective gray opacity

We derive the velocity and mass-loss rate of a steady state Wolf-Rayet (WR) wind, using a nonisotropic diffusion approximation applied to the transfer between strongly overlapping spectral lines. Following the approach of Friend & Castor (1983), the line list is assumed to approximate a statistically parameterized Poisson distribution in frequency, so that photon transport is controlled by an angle-dependent, effectively gray opacity. We show the nonisotropic diffusion approximation yields good agreement with more accurate numerical treatments of the radiative transfer, while providing analytic insight into wind driving by multiple scattering. We illustrate, in particular, that multiple radiative momentum deposition does not require that potons be repeatedly reflected across substantial distances within the spherical envelope, but indeed is greatest when photons undergo a nearly local diffusion, e.g., through scattering by many lines closely spaced in frequency. Our results reiterate the view that the so-called 'momentum problem' of Wolf-Rayet winds is better characterized as an 'opacity problem' of simply identfying enough lines. One way of increasing the number of thick lines in Wolf-Rayet winds is to transfer opacity from saturated to unsaturated lines, yielding a steeper opacity distribution than that found in OB winds. We discuss the implications of this perspective for extending our approach to W-R wind models that incorporate a more fundamental treatment of the ionization and excitation processes that determine the line opacity. In particular, we argue that developing statistical descriptions of the lines to allow an improved effective opacity for the line ensemble would offer several advantages for deriving such more fundamental W-R wind models.

Gayley, Kenneth G.↗

An evaluation of computer codes for simulating the Galileo Probe aerothermal entry environment

The approaches of three computer flow field codes (HYVIS, COLTS, and RASLE), used to determine the Galileo Probe aerothermal environment and its effect on the design of the thermal protection system, are analyzed in order to resolve differences in their predicted results. All three codes account for the hypersonic, massively blown, radiation shock layers, characteristic of Jupiter entry. Significant differences, however, are evident in their solution procedures: the governing conservation equations, the numerical differencing methods, the governing physics (chemical, radiation, diffusion, and turbulence models), and the basic physical data (thermodynamic, transport, chemical, and spectral properties for atomic and molecular species). Solutions are compared for two near peak heating entry conditions for a Galileo Probe baseline configuration, having an initial mass of 242 kg and simulating entry into the Orton nominal atmosphere. The modern numerical methodology of COLTS and RASLE appear to provide an improved capability for coupled flow-field solutions.

Menees, G. P.↗

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 is predicted as a function of galactic longitude for longitudes between 4 and 90 deg; the far-infrared emissivity is predicted as a function of galactocentric distance. 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.↗

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 planes is predicted as a function of galactic longitude for the region between 4 and 90 deg longitude; the far-infrared emissivity is predicted as a function of galactocentric distance. 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.↗

Turbulent Flame Processes Via Diffusion Flame-Vortex Ring Interactions

Flame-vortex interactions are canonical configurations that can be used to study the underlying processes occurring in turbulent reacting flows. This configuration contains many of the fundamental aspects of the coupling between fluid dynamics and combustion that could be investigated with more controllable conditions than are possible under direct investigations of turbulent flames. Diffusion flame-vortex ring interaction contains many of the fundamental elements of flow, transport, combustion, and soot processes found in turbulent diffusion flames. Some of these elements include concentrated vorticity, entrainment and mixing, strain and nonequilibrium phenomena, diffusion and differential diffusion, partial premixing and diluent effects, soot formation and oxidation, and heat release effects. Such simplified flowfield allows the complex processes to be examined more closely and yet preserving the physical processes present in turbulent reacting flows. Furthermore, experimental results from the study of flame-vortex interactions are useful for the validation of numerical simulations and more importantly to deepen our understanding of the fundamental processes present in reacting flows. Experimental and numerical results obtained under microgravity conditions of the diffusion flame-vortex ring interaction are summarized in this paper. Results are obtained using techniques that include Flame Luminosity Imaging (FLI), Laser Soot-Mie Scattering (LSMS), Computational Fluid Dynamics and Combustion (CFDC), and Diode Laser Spectroscopy/Iterative Temperature with Assumed Chemistry (DLS/ITAC).

Dahm, Werner J. A.↗

Tackling the curse of dimensionality in fractional and tempered fractional PDEs with physics-informed neural networks

Fractional and tempered fractional partial differential equations (PDEs) are effective models of long-range interactions, anomalous diffusion, and non-local effects. Traditional numerical methods for these problems are mesh-based, thus struggling with the curse of dimensionality (CoD). Physics-informed neural networks (PINNs) offer a promising solution due to their universal approximation, generalization ability, and mesh-free training. In principle, Monte Carlo fractional PINN (MC-fPINN) estimates fractional derivatives using Monte Carlo methods and thus could lift CoD. However, this may cause significant variance and errors, hence affecting convergence; in addition, MC-fPINN is sensitive to hyperparameters. In general, numerical methods and specifically PINNs for tempered fractional PDEs are under-developed. Herein, we extend MC-fPINN to tempered fractional PDEs to address these issues, resulting in the Monte Carlo tempered fractional PINN (MC-tfPINN). To reduce possible high variance and errors from Monte Carlo sampling, we replace the one-dimensional (1D) Monte Carlo with 1D Gaussian quadrature, applicable to both MC-fPINN and MC-tfPINN. We validate our methods on various forward and inverse problems of fractional and tempered fractional PDEs, scaling up to 100,000 dimensions. Our improved MC-fPINN/MC-tfPINN using quadrature consistently outperforms the original versions in accuracy and convergence speed in very high dimensions.

42 ENGINEERING↗

Two-color pyrometry based flame to fuel surface radiative heat flux diagnostic using flamelets

In this study, a non-intrusive two-color pyrometry based radiative heat flux diagnostic is presented that accounts for the self-absorption effects. This diagnostic relies on looking up the flame temperature and soot volume fraction using the ratio of color intensities obtained from digital cameras and numerical solutions of 1D steady-state diffusion flames. Virtual two color pyrometry of computed flame show a unique one-to-one mapping of strained flames with intensity ratio; thereby allowing flames to be indexed and creation of a two-color flamelet manifold (TCFM). The TCFM is a function of two variables; the ratio of red to green intensities (flame index) and physical distance to the stoichiometric surface defined using a level-set function. The developed diagnostic is used in upward flame spread experiments where an approximate flame hull reconstruction method is applied to provide the 3D level-set function and intensity ratio for use of TCFM to provide radiative absorption properties for 3D ray tracing. Accounting for self- absorption effects provide improved soot volume fraction estimates. In conclusion, the radiative heat flux estimates are found to be in very good agreement with literature and improved results are obtained near the flame base.

42 ENGINEERING↗

Material-dependent photon ionizing radiation effects in Si and GaAs PIN diodes: A numerical investigation

We present a finite-element drift-diffusion-Poisson model in the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework to compare the radiation response of silicon (Si) and gallium arsenide (GaAs) PIN diodes under high-energy photon irradiation. The model solves coupled carrier continuity and Poisson’s equations with Shockley-Read-Hall recombination, and is verified against standard analytical J-V behavior. Using a simplified 1D geometry with ideal Ohmic contacts, we quantify device response under forward and reverse bias with a 100 MeV photon flux. Under forward bias, Si exhibits markedly greater radiation sensitivity than GaAs, including larger increases in current density, stronger local field and carrier-product perturbations, and higher recombination. Under reverse bias, GaAs shows larger radiation-induced photocurrent and broader current-density peaks near junctions, indicating an advantage for photodetection. Integrated steady-state recombination is consistently higher in Si across voltages. Under periodic photon pulses, GaAs produces higher-amplitude photoresponse and settles more rapidly than Si. These results highlight material-dependent trade-offs for radiation-tolerant, high-speed optoelectronics and provide guidance for selecting PIN architectures in aerospace, nuclear, and high-energy physics environments.

36 MATERIALS SCIENCE↗

Continuous polyamorphic transition in high-entropy metallic glass

Polyamorphic transition (PT) is a compelling and pivotal physical phenomenon in the field of glass and materials science. Understanding this transition is of scientific and technological significance, as it offers an important pathway for effectively tuning the structure and property of glasses. In contrast to the PT observed in conventional metallic glasses (MGs), which typically exhibit a pronounced first-order nature, herein we report a continuous PT (CPT) without first-order characteristics in high-entropy MGs (HEMGs) upon heating. This CPT behavior is featured by the continuous structural evolution at the atomic level and an increasing chemical concentration gradient with temperature, but no abrupt reduction in volume and energy. The continuous transformation is associated with the absence of local favorable structures and chemical heterogeneity caused by the high configurational entropy, which limits the distance and frequency of atomic diffusion. As a result of the CPT, numerous glass states can be generated, which provides an opportunity to understand the nature, atomic packing, formability, and properties of MGs. Moreover, this discovery highlights the implication of configurational entropy in exploring polyamorphic glasses with an identical composition but highly tunable structures and properties.

36 MATERIALS SCIENCE↗

Thermal interfaces in dynamic compression experiments

Temperature is a difficult thermodynamic variable to measure in dynamic compression experiments. Optical pyrometry is a general-purpose technique for measuring temperature from a radiant surface, but that surface is often the interface between distinct materials with temperatures that vary spatially along the loading direction. This leads to a fundamental problem: how does the measured interface temperature relate to this temperature profile along the compression axis? Numerical analysis of loading history and thermal diffusion at these interfaces shows that seemingly subtle changes in experiment geometry can lead to very different temperature profiles. We compare these results to laboratory temperature measurements of shock-compressed tin.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

High Temperature Copper Metallization: Demand, Hurdles and Reliability

As newer cells structures come online, the pressing need to replace silver in the metallization pastes has renewed interest in alternative technologies employing base metals. Copper typically leads the charge with its abundance and lower cost but has faced numerous obstacles from relatively higher oxidation and diffusion rates which can damage the lifetime of the devices. In this study, a low-cost alternative to silver metallization pastes has been shown on PERC cells. The screen printable copper paste can be fired in air at temperatures >500 degrees C, and the impact of processing conditions and equipment on the performance and reliability of 274 cm2 cells have been evaluated. Through damp heat testing of micro-modules using 16 cm2 cells, routes that can lead to both the failure and success of durable contacts have been demonstrated.

copper↗

Atlas of lunar radar maps at 70-cm wavelength

The intensity distribution of lunar radar echoes has been mapped for two-thirds of the earth-visible lunar surface at a wavelength of 70 cm. The depolarizing effects of the lunar surface were observed by simultaneously receiving the radar echoes in opposite polarizations. These echoes were mapped with areal resolutions of 25-100 sq km. Mapping with this resolution confirmed that the young craters have enhanced returns. A few craters were found to have enhanced echoes only from their rims. Backscattering differences were also observed between various areas within a mare, between different highland areas, and between maria and adjacent highlands. These scattering differences were interpreted with a simple model, which assumed that the surface backscattered with varying amounts of quasi-specular and diffuse power. Only an increase in the diffuse power was needed to give the numerical values of the enhancements.

Thompson, T. W.↗

A finite difference method of solving anisotropic scattering problems

A new method of solving radiative transfer problems is described including a comparison of its speed with that of the doubling method, and a discussion of its accuracy and suitability for computations involving variable optical properties. The method uses a discretization in angle to produce a coupled set of first-order differential equations which are integrated between discrete depth points to produce a set of recursion relations for symmetric and anti-symmetric angular sums of the radiation field at alternate depth points. The formulation given here includes depth-dependent anisotropic scattering, absorption, and internal sources, and allows arbitrary combinations of specular and non-Lambertian diffuse reflection at either or both boundaries. Numerical tests of the method show that it can return accurate emergent intensities even for large optical depths. The method is also shown to conserve flux to machine accuracy in conservative atmospheres

Barkstrom, B. R.↗

Solar particle propagation from 1 to 5 AU

The radial variation of the interplanetary propagation of solar cosmic rays is investigated on a grand scale by performing a statistical analysis of charged-particle data obtained by Pioneers 10 and 11 from launch to Jovian encounter (i.e., at distances of 1 to 5 AU). A numerical propagation model is adopted which includes diffusion, convection, and adiabatic deceleration, as well as a variable power-law diffusion coefficient. The study of solar particle events is carried out by analyzing individually each of five key parameters that are uniquely defined in each solar particle event: the time to maximum flux starting from the initial release at the sun, the anisotropy at the time to maximum flux, the flux at that time, the time width of the event at half the flux at the time to maximum, and the decay time. Combined results for the proton energy ranges from 3.4 to 5.2 MeV and 24 to 30 MeV are presented, the average radial diffusion coefficients between 1 and 5 AU are determined for each energy range, and it is found that the combination of a near-impusive injection with a temporally constant spectral index produces a reasonably good fit to the high-energy data but not as good a fit to the low-energy data. Implications of the results obtained are discussed for cosmic-ray propagation theory, solar modulation studies, and interplanetary acceleration of solar particles.

Zwickl, R. D.↗

Systems improved numerical differencing analyzer

Program solves physical problems governed by diffusion-type equations, provided that equations can be modeled by lumped-parameter representation. Program is used for thermal analysis, and could be adapted to solve Fourier, Poisson, and Laplace differential equations. Program is in FORTRAN IV and Assembler for execution on UNIVAC 1100-series or CYBER 175.

Source record↗

Application of a Galerkin finite element method to atmospheric transport problems

Numerical simulation of the movement of a contaminant within the atmosphere presents difficulties due to the multidimensionality of the problem and the fact that the horizontal transport is usually convection dominated, that the boundary conditions are mixed, and that both slow and fast atmospheric chemical reactions can be important. In this study, numerical experiments using a Crank-Nicolson Galerkin finite element method to solve the time-dependent partial differential equations demonstrate the applicability and accuracy of this method for the variety of conditions encountered in atmospheric pollutant modeling. The Crank-Nicolson Galerkin method using piecewise linear, piecewise cubic Hermite polynomials, and upwind finite elements is shown to accurately model the pure convection of initial wave forms. Numerical results studying the interactions of convection, diffusion, chemical reaction, pollutant removal, and the effects of contaminant emission source strength, source location and multiple sources are also presented.

Carmichael, G. R.↗