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 1,171 records · Page 65

All sources NASA NTRS DOE OSTI

NASA NTRS1988

Finite element procedures for time-dependent convection-diffusion-reaction systems

New finite element procedures based on the streamline-upwind/Petrov-Galerkin formulations are developed for time-dependent convection-diffusion-reaction equations. These procedures minimize spurious oscillations for convection-dominated and reaction-dominated problems. The results obtained for representative numerical examples are accurate with minimal oscillations. As a special application problem, the single-well chemical tracer test (a procedure for measuring oil remaining in a depleted field) is simulated numerically. The results show the importance of temperature effects on the interpreted value of residual oil saturation from such tests.

Tezduyar, T. E.↗

NASA NTRS1990

Spectral determinations for discrete sources with EGRET

The ability of the EGRET (Energetic Gamma-Ray Experimental Telescope) to determine the spectral parameters of point sources in 14-day exposures, as planned for the initial survey phase of the GRO (Gamma Ray Observatory) mission, is explored by numerical simulation. Results are given for both galactic and extragalactic objects as a function of source strength and for representative levels of diffuse background emission.

Hughes, E. B.↗

NASA NTRS1991

Navier-Stokes simulation of inlet flow fields-experimental verification

The objectives were to develop a state-of-the-art numerical capability for the prediction of internal and external flow fields, validate the Navier-Stokes solution for internal flows, treat realistic diffuser inlet systems, and provide a technical base for inlet system design. Progress is reported in viewgraph form.

Chyu, Wei J.↗

NASA NTRS1997

Extinguishment of a Diffusion Flame Over a PMMA Cylinder by Depressurization in Low-Gravity

The behavior of flames in low-speed flows in low-gravity is relevant to spacecraft fire safety. Previous work has shown that flames in the presence of low-speed forced flows in low-gravity may be more flammable than in a forced flow of the same magnitude in normal gravity. Additionally, fire suppression plans for the International Space Station include the use of venting (depressurization) as an emergency option for extinguishing fires. This procedure would induce flows in the affected compartment that could temporarily intensify the fire, as was observed in flammability tests of solids conducted on board Skylab. Despite a general Understanding, current knowledge of the combined effects of reduced pressure and forced flow on a burning solid in low-gravity is inadequate for the design of a venting extinguishment system. Previous studies in low-g have examined flammability limits for thermally thin solids. However, there are differences when burning thick materials because the interior solid-phase temperature continuously changes, which affects the percentage of gas-phase heat feedback to the solid-phase. Changes in the heat feedback to the solid-phase can affect the flammability characteristics of the material. In the current work, the extinction of a diffusion flame burning over PMMA (polymethyl methacrylate) cylinders during depressurization with a low-speed cross flow was examined experimentally and via numerical simulations.

Goldmeer, Jeffrey S.↗

NASA NTRS1999

Diffusion Flame Tip Instabilities of a Wide Sample in Microgravity

This work is a study of diffusion flame tip instabilities of a wide thermally thin sample in a microgravity environment. The purpose of this work is to determine the thermodiffusive and hydrodynamic instability mechanisms which cause a spreading diffusion flame to become corrugated and later fragment into smaller separate flames that we call flamelets. These thermodiffusive and hydrodynamic mechanisms are prominent forces in the near limit because buoyancy is quashed in the microgravity environment. A fundamental knowledge of the underlying mechanisms controlling the extinction limit of a diffusion flame can be applied to any real world device involving a spreading diffusion flame as well as fire safety applications on earth and in space. Our work, which combines microgravity experiments, theoretical models, and numerical computation, is a joint effort between NASA and Michigan State University. Experiments in the 2.2 second drop tower at NASA Lewis Research Center in the Microgravity Combustion Branch have been a preliminary step in locating a region where these instabilities are prevalent. These positive results will be discussed as well as the second phase of our project and future plans and recommendations.

Oravecz, Lisa↗

NASA NTRS1997

Studying Room Acoustics using a Monopole-Dipole Microphone Array

The use of a soundfield microphone for examining the directional nature of a room impulse response was reported recently. By cross-correlating monopole and co-located dipole microphone signals aligned with left-right, up-down, and front-back axes, a sense of signal direction of arrival is revealed. The current study is concerned with the array's ability to detect individual reflections and directions of arrival, as a function of the cross-correlation window duration. If is window is too long, weak reflections are overlooked; if too short, spurious detections result. Guidelines are presented for setting the window width according to perceptual criteria. Formulas are presented describing the accuracy with which direction of arrival can be estimated as a function of room specifics and measurement noise. The direction of arrival of early reflections is more accurately determined than that of later reflections which are quieter and more numerous. The transition from a fairly directional sound field at the beginning of the room impulse response to a uni-directional diffuse field is examined. Finally, it is shown that measurements from additional dipole orientations can significantly improve the ability to detect reflections and estimate their directions of arrival.

Begault, Durand R.↗

NASA NTRS2003

Particle Engulfment and Pushing By Solidifying Interfaces

The study of particle behavior at solid/liquid interfaces (SLI s) is at the center of the Particle Engulfment and Pushing (PEP) research program. Interactions of particles with SLI s have been of interest since the 1960 s, starting with geological observations, i.e., frost heaving. Ever since, this field of research has become significant to such diverse areas as metal matrix composite materials, fabrication of superconductors, and inclusion control in steels. The PEP research effort is geared towards understanding the fundamental physics of the interaction between particles and a planar SLI. Experimental work including 1-g and mu-g experiments accompany the development of analytical and numerical models. The experimental work comprised of substantial groundwork with aluminum (Al) and zinc (Zn) matrices containing spherical zirconia particles, mu-g experiments with metallic Al matrices and the use of transparent organic metal-analogue materials. The modeling efforts have grown from the initial steady-state analytical model to dynamic models, accounting for the initial acceleration of a particle at rest by an advancing SLI. To gain a more comprehensive understanding, numerical models were developed to account for the influence of the thermal and solutal field. Current efforts are geared towards coupling the diffusive 2-D front tracking model with a fluid flow model to account for differences in the physics of interaction between 1-g and -g environments. A significant amount of this theoretical investigation has been and is being performed by co-investigators at NASA MSFC.

Source record↗

NASA NTRS2015

High-Order Hyperbolic Residual-Distribution Schemes on Arbitrary Triangular Grids

In this paper, we construct high-order hyperbolic residual-distribution schemes for general advection-diffusion problems on arbitrary triangular grids. We demonstrate that the second-order accuracy of the hyperbolic schemes can be greatly improved by requiring the scheme to preserve exact quadratic solutions. We also show that the improved second-order scheme can be easily extended to third-order by further requiring the exactness for cubic solutions. We construct these schemes based on the LDA and the SUPG methodology formulated in the framework of the residual-distribution method. For both second- and third-order-schemes, we construct a fully implicit solver by the exact residual Jacobian of the second-order scheme, and demonstrate rapid convergence of 10-15 iterations to reduce the residuals by 10 orders of magnitude. We demonstrate also that these schemes can be constructed based on a separate treatment of the advective and diffusive terms, which paves the way for the construction of hyperbolic residual-distribution schemes for the compressible Navier-Stokes equations. Numerical results show that these schemes produce exceptionally accurate and smooth solution gradients on highly skewed and anisotropic triangular grids, including curved boundary problems, using linear elements. We also present Fourier analysis performed on the constructed linear system and show that an under-relaxation parameter is needed for stabilization of Gauss-Seidel relaxation.

Mazaheri, Alireza↗

NASA NTRS1999

Turbulent Diffusivities for Momentum, Heat, Salt and Passive Scalars

A program began ten years ago to build a turbulence model to describe high Reynolds numbers flows. Specifically, the aim was to devise a model that would satisfy two basic conditions: reproduce available turbulence data (laboratory, numerical simulations, etc.) concerning shear driven flows, buoyancy driven flows, 2D turbulence, freely decaying turbulence etc., and yet be manageable so as to be used, for example, in General Circulation Models (GCM'S). The model was presented in a series of papers that have appeared in Physics of Fluids since 1996. A total of about 80 turbulence statistics were reproduced. The model has no adjustable parameters. The next step was to apply the model to construct the vertical diffusities K for moment K(sub m), temperature K(sub h), salt K(sub s). and passive scalars K(sub c). First, we took K(sub s) = K(sub h) and tested the model using the GFDL ocean model. The results for the profiles of T and S vs. depth are indistinguishable from those derived using the latest model, the KPP model by the NCAR group. Presently, we are running the same Geophysical Fluid Dynamics Laboratory (GFDL) code relaxing the assumptive K(sub s) = K(sub h). Indeed, the turbulence model yields a salt diffusivity that depends on Ri and R rho (= Beta delta S/delta z/alpha delta T/delta z) in such a way that K(sub s) may be quite different from K(sub h). Salt fingers and double diffusivity laboratory data are reproduced. Results from the ocean model will be available shortly. Finally, we are trying to derive the horizontal diffusivities with the goal of providing a physically acceptable representation of mesoscale eddies. The recently suggested GMW parameterization has improved several O-GCM results and the goal here is to try to derive/justify it from a turbulence model and/or propose improvements/modifications. Theoretical work is in progress.

Canuto, V. M.↗

DOE OSTI2018

A front-fixing finite element method for pricing American options under regime-switching jump-diffusion models

A front-fixing finite element method is applied to solve the partial integro-differential equations (PIDEs) arising in pricing American options under Markov-modulated jump-diffusion models with the free boundaries feature. For this purpose, we used a front-fixing method to transform the pricing problem into a nonlinear parabolic integro-differential equation on a fixed domain. Then the variational form of the resulting problem is solved by a finite element method. Under some appropriate assumptions, we establish the stability of the method and illustrate some numerical results to examine the rate of convergence of the proposed method for the pricing problem and compare its accuracy to some recent works on pricing American options under regime-switching jump-diffusion models.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1975

Heating of solid targets with laser pulses

Analytical and numerical solutions to the heat-conduction equation are obtained for the heating of absorbing media with pulsed lasers. The spatial and temporal form of the temperature is determined using several different models of the laser irradiance. Both surface and volume generation of heat are discussed. It is found that if the depth of thermal diffusion for the laser-pulse duration is large compared to the optical-attenuation depth, the surface- and volume-generation models give nearly identical results. However, if the thermal-diffusion depth for the laser-pulse duration is comparable to or less than the optical-attenuation depth, the surface-generation model can give significantly different results compared to the volume-generation model. Specific numerical results are given for a tungsten target irradiated by pulses of different temporal durations and the implications of the results are discussed with respect to the heating of metals by picosecond laser pulses.

Bechtel, J. H.↗

DOE OSTI2017

The Weakly Nonlinear Magnetorotational Instability in a Local Geometry

The magnetorotational instability (MRI) is a fundamental process of accretion disk physics, but its saturation mechanism remains poorly understood despite considerable theoretical and computational effort. We present a multiple-scales analysis of the non-ideal MRI in the weakly nonlinear regime—that is, when the most unstable MRI mode has a growth rate asymptotically approaching zero from above. Here, we develop our theory in a local, Cartesian channel. Our results confirm the finding by Umurhan et al. that the perturbation amplitude follows a Ginzburg–Landau equation. We further find that the Ginzburg–Landau equation will arise for the local MRI system with shear-periodic boundary conditions, when the effects of ambipolar diffusion are considered. A detailed force balance for the saturated azimuthal velocity and vertical magnetic field demonstrates that, even when diffusive effects are important, the bulk flow saturates via the combined processes of reducing the background shear and rearranging and strengthening the background vertical magnetic field. We directly simulate the Ginzburg–Landau amplitude evolution for our system, and demonstrate the pattern formation our model predicts on long scales of length- and timescales. We compare the weakly nonlinear theory results to a direct numerical simulation of the MRI in a thin-gap Taylor Couette flow.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1995

Numerical modeling of Bridgman growth of PbSnTe in a magnetic field

In this work we study heat and mass transport, fluid motion, and solid/liquid phase change in the process of steady Bridgman growth of Pb(.8)Sn(.2)Te (LTT) in an axially-imposed uniform magnetic field under terrestrial and microgravity conditions. In particular, this research is concerned with the interrelationships among segregation, buoyancy-driven convection, and magnetic damping in the LTT melt. The main objectives are to provide a quantitative understanding of the complex transport phenomena during solidification of the nondilute binary of LTT, to provide estimates of the strength of magnetic field required to achieve the desired diffusion-dominated growth, and to assess the role of magnetic damping for space and earth based control of the buoyancy-induced convection. The problem was solved by using FIDAP and numerical results for both vertical and horizontal growth configurations with respect to the acceleration of gravity vector are presented.

Yao, Minwu↗

DOE OSTI2021

An asymptotically compatible approach for Neumann-type boundary condition on nonlocal problems

In this paper we consider 2D nonlocal diffusion models with a finite nonlocal horizon parameter δ characterizing the range of nonlocal interactions, and consider the treatment of Neumann-like boundary conditions that have proven challenging for discretizations of nonlocal models. We propose a new generalization of classical local Neumann conditions by converting the local flux to a correction term in the nonlocal model, which provides an estimate for the nonlocal interactions of each point with points outside the domain. While existing 2D nonlocal flux boundary conditions have been shown to exhibit at most first order convergence to the local counter part as δ → 0, the proposed Neumann-type boundary formulation recovers the local case as O(δ 2 ) in the L∞(Ω) norm, which is optimal considering the O(δ 2 ) convergence of the nonlocal equation to its local limit away from the boundary. We analyze the application of this new boundary treatment to the nonlocal diffusion problem, and present conditions under which the solution of the nonlocal boundary value problem converges to the solution of the corresponding local Neumann problem as the horizon is reduced. To demonstrate the applicability of this nonlocal flux boundary condition to more complicated scenarios, we extend the approach to less regular domains, numerically verifying that we preserve second-order convergence for non-convex domains with corners. Finally, based on the new formulation for nonlocal boundary condition, we develop an asymptotically compatible meshfree discretization, obtaining a solution to the nonlocal diffusion equation with mixed boundary conditions that converges with O(δ 2 ) convergence.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1999

Influence of Buoyant Convection on the Stability of Enclosed Laminar Flames

Enclosed diffusion flames are commonly found in practical combustion systems, such as the power-plant combustor, gas turbine combustor, and jet engine after-burner. In these systems, fuel is injected into a duct with a co-flowing or cross-flowing air stream. In combustors, this flame is anchored at the burner (i.e., fuel jet inlet) unless adverse conditions cause the flame to lift off or blow out. Investigations of burner stability study the lift off, reattachment, and blow out of the flame. There have been numerous studies of flame stability. Relatively few studies have investigated the stability of flames with an oxidizer co-flow, compared with the number of studies on (nearly) free jet diffusion flames. The air flow around the fuel jet can significantly alter the lift off, reattachment and blow out of the jet diffusion flame. In normal gravity, however, the effects of the air flow on flame stability are often complicated by the presence of buoyant convection. A comparison of normal-gravity and microgravity flames can provide clear indication of the influence of forced and buoyant flows on the flame stability. The overall goal of the Enclosed Laminar Flames (ELF) research, described at the following URL site: http://zeta.lerc.nasa.gov/expr/elf.htm, is to improve our understanding of the effects of buoyant convection on the structure and stability of co-flow diffusion flames.

Brooker, John E.↗

NASA NTRS1993

Evaporation and combustion of LOX under supercritical and subcritical conditions

The objective is to study the evaporation and combustion of LOX under supercritical and subcritical conditions both experimentally and theoretically. In the evaporation studies, evaporation rate and surface temperature were measured when LOX vaporizing in helium environments at pressures ranging from 5 to 68 atm. A Varian 3700 gas chromatograph was employed to measure the oxygen concentration above the LOX surface. For the combustion tests, high-magnification video photography was used to record direct images of the flame shape of a LOX/H2/He laminar diffusion flame. The gas composition in the post-flame region is also being measured with the gas sampling and chromatography analysis. These data are being used to validate the theoretical model. A comprehensive theoretical model with the consideration of the solubility of ambient gases as well as variable thermophysical properties was formulated and solved numerically to study the gasification and burning of LOX at elevated pressures. The calculated flame shape agreed reasonably well with the edge of the observed luminous flame surface. The effect of gravity on the flame structure of laminar diffusion flames was found to be significant. In addition, the predicted results using the flame-sheet model were compared with those based upon full equilibrium calculations (which considered the formation of intermediate species) at supercritical pressures. Except at the flame front where temperature exceeded 2,800 K, the flame-sheet and equilibrium solutions in terms of temperature distributions were in very close agreement. The temperature deviation in the neighborhood of the flame front is caused by the effect of high-temperature dissociation.

Yang, A. S.↗

DOE OSTI2020

Affine Approximation of Parametrized Kernels and Model Order Reduction for Nonlocal and Fractional Laplace Models

In this work, we consider parametrized problems driven by spatially nonlocal integral operators with parameter-dependent kernels. In particular, kernels with varying nonlocal interaction radius $\delta > 0$ and fractional Laplace kernels, parametrized by the fractional power $s\in(0,1)$, are studied. Furthermore, in order to provide an efficient and reliable approximation of the solution for different values of the parameters, we develop the reduced basis method as a parametric model order reduction approach. Major difficulties arise since the kernels are not affine in the parameters, singular, and discontinuous. Moreover, the spatial regularity of the solutions depends on the varying fractional power $s$. To address this, we derive regularity and differentiability results with respect to $\delta$ and $s$, which are of independent interest for other applications such as optimization and parameter identification. We then use these results to construct affine approximations of the kernels by local polynomials. Finally, we certify the method by providing reliable a posteriori error estimators, which account for all approximation errors, and support the theoretical findings by numerical experiments.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Error Estimates for the Optimal Control of a Parabolic Fractional PDE

In this work, we consider the integral definition of the fractional Laplacian and analyze a linear-quadratic optimal control problem for the so-called fractional heat equation; control constraints are also considered. We derive existence and uniqueness results, first order optimality conditions, and regularity estimates for the optimal variables. To discretize the state equation we propose a fully discrete scheme that relies on an implicit finite difference discretization in time combined with a piecewise linear finite element discretization in space. We derive stability results and a novel $L^2(0,T;L^2(\Omega))$ a priori error estimate. On the basis of the aforementioned solution technique, we propose a fully discrete scheme for our optimal control problem that discretizes the control variable with piecewise constant functions, and we derive a priori error estimates for it. We illustrate the theory with one- and two-dimensional numerical experiments.

97 MATHEMATICS AND COMPUTING↗