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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

NASA NTRS1992

Approximate Solutions Of Equations Of Steady Diffusion

Rigorous analysis yields reliable criteria for "best-fit" functions. Improved "curve-fitting" method yields approximate solutions to differential equations of steady-state diffusion. Method applies to problems in which rates of diffusion depend linearly or nonlinearly on concentrations of diffusants, approximate solutions analytic or numerical, and boundary conditions of Dirichlet type, of Neumann type, or mixture of both types. Applied to equations for diffusion of charge carriers in semiconductors in which mobilities and lifetimes of charge carriers depend on concentrations.

Edmonds, Larry D.↗

NASA NTRS1994

The formation of sharp structures by ambipolar diffusion

The effect of ambipolar diffusion is investigated using simple numerical models. Examples are shown where sharp structures develop around magnetic nulls. In contrast to the case of ordinary diffusion, the magnetic field topology is conserved by ambipolar diffusion. This is demonstrated in an example where differential rotation winds up an initially uniform magnetic field and brings oppositely oriented field lines close together. It is argued that ambipolar diffusion produces structures of scales small enough for reconnection to occur.

Brandenburg, Axel↗

DOE OSTI2020

Diffuser Braze Development

At the request of SRTE, through the PDRD project, SRNL is developing a vacuum brazing process to facilitate the fabrication of hydrogen purification systems. In particular, this effort is directed at tube-style hydrogen diffusers. The diffusers are comprised of numerous palladium-silver tubes that are brazed to a nickel fitting, which is subsequently welded to a stainless steel tube, fitting, or tube sheet. Historically, torch brazing has been used as a heat source; for this work vacuum furnace brazing was suggested to improve the uniformity of the process. This report describes the braze process development for braze selection and simple nickel fitting to tube brazing. Five braze alloys were selected and sessile drop testing was conducted over a range of temperatures from 25 to 75° C over the liquidus temperature. Samples were characterized for wetting angle, erosion, and brazability. The cup-to-tube braze joints were inspected using digital radiography, ultrasonic testing, and metallographic examination. Two of the five alloys exhibited excellent fabrication properties and are considered for further development.

36 MATERIALS SCIENCE↗

DOE OSTI2021

NRC Multiphysics Analysis Capability Deployment FY21: Part 3

This report details the progress and activities of Idaho National Laboratory (INL) on the Nuclear Regulatory Commission (NRC) project “Development and Modeling Support for Advanced Non-Light Water Reactors.” The deliverables completed for this report are: Deliverable 1c: the capability to model gas mixtures was added to Pronghorn. A test problem mimicking the conditions achieved in a depressurized loss of forced cooling (DLOFC) event was solved with both RELAP-5 and Pronghorn. Pronghorn employed a finite vol ume method with the Kurganov-Tadmor discretization. The comparison between the mass fraction spatial profiles computed with RELAP-5 and Pronghorn clearly shows the presence of numerical artifacts (i.e., overly diffusive behavior at low Mach numbers). We confirmed that the problem disappears at higher Mach numbers. We recommend future work on the implementation of a low Mach finite volume formulation to better treat low Mach number problems. Deliverable 2a: we demonstrated two approaches to model the radiation/conduction/natural convection heat transfer across a stagnant gas for the PBMR-400 design using Pronghorn. The first approach is based on the net radiation method, which relies on the computation of view factors with the Multiphysics Object-Oriented Simulation Environment (MOOSE) ray tracing capability. The second method is a traditional thermal resistance approach. The test problems include both 2D and 3D geometries. In all cases, the results show very good agreement during a DLOFC transient. This confirms that the faster thermal resistance method produces solutions that are equivalent to the net radiation method for this geometry. Deliverable 3d: we demonstrated the use of the advection kernel for the delayed neutron precursor equation in Griffin with a 2D MSFR model. The results appear physical but further verification is recommended. We also recommend the addition of conjugate heat transfer to compute the temperatures and model the thermomechanic behavior of the reflectors and other structures. Significant memory and performance issues were encountered in the 3D axisymmetric model. Future work is recommended in this area. Task 8g: this task allows multidimensional MOOSE applications to be coupled to system codes (RELAP-7 and SAM). We implemented a faster multiphysics iteration coupling algorithm, which provides an overall 6× acceleration of the 3D-1D coupling of the core multidi- mensional fluid flow solver and the 1D primary and secondary loop model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1971

Zero-gravity transient thermal mixing simulation

The experimental program described is an outgrowth of independent investigations into alternate redesign concepts for the Apollo SM cryogenic oxygen storage system. The experiments were continued, after the redesign was established, to provide physical insight into transient thermal mixing in zero-gravity and to aid in the characterization of the system performance in flight. Zero-gravity heat transfer and fluid mixing were simulated experimentally through an analogy between unsteady heat conduction and species diffusion. To further support numerical analyses of the cryogenic oxygen storage system, the experimental investigation was extended to include a cubical tank geometry, representative of existing numerical models. In general, the transient flow patterns in the cubical tank are far more complex than those of the spherical tank and the extent of fluid mixing is significantly greater but less repeatable.

Lands, J. F., Jr.↗

NASA NTRS1974

Scattering and absorption from poly-dispersed aerosols

Results from numerical computations of the diffuse sky radiances and degree of polarization are given. These results were obtained by solving the equation of radiative transfer in the formulation of Eschelbach. In these computations, multiple scattering and absorption by aerosol particles were considered. In addition, results from experimental measurements of the complex refractive index of aerosol samples at different relative humidities are presented. From the radiation flux divergences which were computed based on the determined properties of atmospheric aerosol particles, atmospheric heating rates were derived which were found to be comparable to the heating rates by water vapor.

Bullrich, K.↗

NASA NTRS1975

Turbulence in a gaseous hydrogen-liquid oxygen rocket combustion chamber

The intensity of turbulence and the Lagrangian correlation coefficient for a LOX-GH2 rocket combustion chamber was determined from experimental measurements of tracer gas diffusion. A combination of Taylor's turbulent diffusion theory and a numerical method for solving the conservation equations of fluid mechanics was used to calculate these quantities. Taylor's theory was extended to consider the inhomogeneity of the turbulence field in the axial direction of the combustion chamber, and an exponential function was used to represent the Lagrangian correlation coefficient. The results indicate that the value of the intensity of turbulence reaches a maximum of 14% at a location about 7" downstream from the injector. The Lagrangian correlation coefficient associated with this value is given by the above exponential expression where alpha = 10,000/sec.

Lebas, J.↗

NASA NTRS1980

Computer program documentation D1FLTD to drive SINDA boundary nodes: User's guide

The thermal model correlation process begins when measured thermocouple data is available from the orbital flight tests of the shuttle. For this effort, it is necessary to convert some of the system improved numerical differencing analyzer (SINDA) diffusion or arithmetic nodes to boundary nodes and then drive these boundary nodes to the temperature profile of a flight measurement. An efficient way to provide this capability within the SINDA and OFT software systems is to provide a new SINDA routine, D1FLTD, for use in VARIABLES 1 of SINDA, to access the processed (word-addressable) orbital data reduction center flight data and store the appropriate measurement temperature in the desired SINDA temperature location. The ODRC flight data that is to be used for driving the boundary nodes must be assigned a logical unit number and must reside on a word-addressable file. The user must also provide two SINDA constants for the word positions of the first and last words of the temperature record for each measurement identifier (MID), i.e. each call to D1FLTD, used in the model. D1FLTD is then called from the VARIABLES 1 block to obtain the SINDA boundary node temperature for any MID on the file at any time point.

Damico, S. J.↗

NASA NTRS1981

Mixing of a passive scalar in isotropic and sheared homogeneous turbulence

In order to calculate the velocity and scalar fields, the three dimensional, time-dependent equations of motion and the diffusion equation were solved numerically. The following cases were treated: isotropic, homogeneous turbulence with decay of a passive scalar; and homogeneous turbulent shear flow with a passive scalar whose mean varies linearly in the spanwise direction. The solutions were obtained at relatively low Reynolds numbers so that all of the turbulent scales could be resolved without modeling. Turbulent statistics such as integral length scales, Taylor microscales, Kolmogorov length scale, one- and two-point correlations of velocity-velocity and velocity-scalar, turbulent Prandtl/Schmidt number, r.m.s. values of velocities, the scalar quantity and pressure, skewness, decay rates, and decay exponents were calculated. The results are compared with the available expermental results, and good agreement is obtained.

Shirani, E.↗

NASA NTRS1982

Heat generation in aircraft tires under free rolling conditions

A method was developed for calculating the internal temperature distribution in an aircraft tire while free rolling under load. The method uses an approximate stress analysis of each point in the tire as it rolls through the contact patch, and from this stress change the mechanical work done on each volume element may be obtained and converted into a heat release rate through a knowledge of material characteristics. The tire cross-section is then considered as a body with internal heat generation, and the diffusion equation is solved numerically with appropriate boundary conditions of the wheel and runway surface. Comparison with data obtained with buried thermocouples in tires shows good agreement.

Clark, S. K.↗

NASA NTRS1983

Heat generation in aircraft tires

Clark, S. K.↗

NASA NTRS1985

The theoretical and experimental investigation of cosmic ray Forbush-effects

The theoretical results of analysis of the expected spatial distributions of density, gradients and anisotropy of cosmic rays, obtained on the basis of the numerical solution of anisotropic diffusion equation in the presence of the disturbances of shock wave type in the interplanetary space are presented. The theoretical calculations on the definition of the energy spectrum and anisotropy of galactic cosmic rays during Forbush effect are compared with the experimental data.

Alania, M. V.↗

NASA NTRS1985

Heat generation in aircraft tires

Clark, S. K.↗

NASA NTRS1986

Particle acceleration at a termination shock. I - Application to the solar wind and the anomalous component

The results of a numerical study of the diffusive acceleration of charged particles at the termination shock of the solar wind are reported. In the model a realistic magnetic field structure is employed which is similar to that observed in the solar wind. In addition to causing spatial variation of the diffusion tensor, a major effect of the magnetic field is to cause guiding-center drifts of the accelerated particles, both in the solar wind and at the shock. It is demonstrated that the inclusion of the drifts has a large effect on the acceleration. It is concluded, furthermore, that acceleration at the termination shock, in conjunction with drifts, can explain several observed features of the anomalous component.

Jokipii, J. R.↗

NASA NTRS1989

Solution techniques for incompressible flow problems

A three-step Petrov-Galerkin (PG)/operator spliting scheme for the time-dependent incompressible Navier-Stokes equations is proposed. Each time step is split into two Stokes problems and one nonlinear convection-diffusion problem. Using a PG technique on the two outer Stokes problems ensures a stable scheme despite equal-order interpolation, while using a streamline upwind PG scheme on the inner convection-diffusion problem ensures a numerically stable solution at high Reynolds numbers. Numerical tests of this method have been carried out.

Tezduyar, T. E.↗

NASA NTRS1989

Assessment of numerical methods for the solution of fluid dynamics equations for nonlinear resonance systems

The capability of accurate nonlinear flow analysis of resonance systems is essential in many problems, including combustion instability. Classical numerical schemes are either too diffusive or too dispersive especially for transient problems. In the last few years, significant progress has been made in the numerical methods for flows with shocks. The objective was to assess advanced shock capturing schemes on transient flows. Several numerical schemes were tested including TVD, MUSCL, ENO, FCT, and Riemann Solver Godunov type schemes. A systematic assessment was performed on scalar transport, Burgers' and gas dynamic problems. Several shock capturing schemes are compared on fast transient resonant pipe flow problems. A system of 1-D nonlinear hyperbolic gas dynamics equations is solved to predict propagation of finite amplitude waves, the wave steepening, formation, propagation, and reflection of shocks for several hundred wave cycles. It is shown that high accuracy schemes can be used for direct, exact nonlinear analysis of combustion instability problems, preserving high harmonic energy content for long periods of time.

Przekwas, A. J.↗

NASA NTRS1990

Application of implicit numerical techniques to the solution of the three-dimensional diffusion equation

Implicit techniques for calculating three-dimensional, time-dependent heat diffusion in a cube are tested with emphasis on storage efficiency, accuracy, and speed of calculation. For this purpose, a tensor product technique with both Chebyshev collocation and finite differences and a generalized conjugate gradient technique with finite differences are used in conjunction with Crank-Nicolson discretization. An Euler explicit finite difference calculation is performed for use as a benchmark. The implicit techniques are found to be competitive with the Euler explicit method in terms of storage efficiency and speed of calculation and offer advantages both in accuracy and stability. Mesh stretching in the finite difference calculations is shown to markedly improve the accuracy of the solution.

Peltier, Leonard Joel↗

NASA NTRS1990

A two-dimensional kinematic dynamo model of the ionospheric magnetic field at Venus

The results of a high-resolution, two-dimensional, time dependent, kinematic dynamo model of the ionospheric magnetic field of Venus are presented. Various one-dimensional models are considered and the two-dimensional model is then detailed. In this model, the two-dimensional magnetic induction equation, the magnetic diffusion-convection equation, is numerically solved using specified plasma velocities. Origins of the vertical velocity profile and of the horizontal velocities are discussed. It is argued that the basic features of the vertical magnetic field profile remain unaltered by horizontal flow effects and also that horizontal plasma flow can strongly affect the magnetic field for altitudes above 300 km.

Cravens, T. E.↗