Dimensional scaling and numerical similarity in hyperbolic method for diffusion
Explore the source record for details and available documents.
Hiroaki Nishikawa↗
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.
A Numerical and Experimental Study of Coflow Laminar Diffusion Flames: Effects of Gravity and Inlet Velocity
In this work, the influence of gravity, fuel dilution, and inlet velocity on the structure, stabilization, and sooting behavior of laminar coflow methane-air diffusion flames was investigated both computationally and experimentally. A series of flames measured in the Structure and Liftoff in Combustion Experiment (SLICE) was assessed numerically under microgravity and normal gravity conditions with the fuel stream CH4 mole fraction ranging from 0.4 to 1.0. Computationally, the MC-Smooth vorticity-velocity formulation of the governing equations was employed to describe the reactive gaseous mixture; the soot evolution process was considered as a classical aerosol dynamics problem and was represented by the sectional aerosol equations. Since each flame is axisymmetric, a two-dimensional computational domain was employed, where the grid on the axisymmetric domain was a nonuniform tensor product mesh. The governing equations and boundary conditions were discretized on the mesh by a nine-point finite difference stencil, with the convective terms approximated by a monotonic upwind scheme and all other derivatives approximated by centered differences. The resulting set of fully coupled, strongly nonlinear equations was solved simultaneously using a damped, modified Newton's method and a nested Bi-CGSTAB linear algebra solver. Experimentally, the flame shape, size, lift-off height, and soot temperature were determined by flame emission images recorded by a digital camera, and the soot volume fraction was quantified through an absolute light calibration using a thermocouple. For a broad spectrum of flames in microgravity and normal gravity, the computed and measured flame quantities (e.g., temperature profile, flame shape, lift-off height, and soot volume fraction) were first compared to assess the accuracy of the numerical model. After its validity was established, the influence of gravity, fuel dilution, and inlet velocity on the structure, stabilization, and sooting tendency of laminar coflow methane-air diffusion flames was explored further by examining quantities derived from the computational results.
microgravity↗
Effects of Structure and Hydrodynamics on the Sooting Behavior of Spherical Microgravity Diffusion Flames
Recent experimental, numerical and analytical work has shown that the stoichiometric mixture fraction (Z(sub st)) can have a profound effect on soot formation in diffusion flames. These findings were obtained at constant flame temperature (T(sub ad)), employing the approach described in Du and Axelbaum (1995, 1996). For example, a fuel mixture containing 1 mole of ethylene and 11.28 moles of nitrogen burning in pure oxygen ((Z(sub st)) = 0.78) has the same adiabatic flame temperature (2370 K) as that of pure ethylene burning in air ((Z(sub st)) = 0.064). An important finding of these works was that at sufficiently high (Z(sub st)), flames remain blue as strain rate approaches zero in counterflow flames, or as flame height and residence time approach infinity in coflowing flames. Lin and Faeth (1996a) coined the term permanently blue to describe such flames. Two theories have been proposed to explain the appearance of permanently-blue flames at high (Z(sub st)). They are based on (1) hydrodynamics and (2) flame structure. Previous experimental studies in normal gravity are not definitive as to which, if either, mechanism is dominant because both hydrodynamics and structure suppress soot formation at high (Z(sub st)) in coflowing and counterflowing diffusion flames. In counterflow flames with (Z(sub st)) < 0.5 streamlines at the flame sheet are directed toward the fuel. Newly formed soot is convected into richer regions, favoring soot growth over oxidation. For (Z(sub st)) > 0.5, convection at the flame is toward the oxidizer, thus enhancing soot oxidization. Thus, in counterflow flames, hydrodynamics causes soot to be convected towards the oxidizer at high (Z(sub st)) which suppresses soot formation. Axelbaum and co-workers maintain that while the direction of convection can impact soot growth and oxidation, these processes alone cannot cause permanently-blue flames. Soot growth and oxidation are dependent on the existence of soot particles and the presence of soot is invariably accompanied by yellow luminosity. Soot-particle inception, on the other hand, arises from gas-phase reactions and its dependence on flow direction is weak, similar to that of other gas-phase reactions in flames. For example, when the flame moves across the stagnation plane no significant changes in flame chemistry are observed. Furthermore, since the soot-inception zone has a finite thickness, soot has been produced in counterflow flames with (Z(sub st)) > 0.5. For large (Z(sub st)) the fuel concentration decreases and oxygen concentration increases in the soot forming regions of the flame. This yields a shift in the OH profile toward the fuel side of the flame, and this shift can dramatically influence soot inception because it essentially narrows the soot inception zone. Soot-free (permanently-blue) conditions can be realized when the structure of the flame is adjusted to the extent that significant oxidizing species exist on the fuel side of the flame at temperatures above the critical temperature for soot inception, ca. 1250 K. In previously considered flames it was impossible to independently vary flame structure and convection direction. In contrast, spherical diffusion flames (which generally require microgravity) allow both properties to be varied independently. We altered structure (Z(sub st)) by exchanging inert between the oxidizer and the fuel and we independently varied convection direction at the flame sheet by interchanging the injected and ambient gases. In this work we established four flames: (a) ethylene issuing into air, (b) diluted ethylene issuing into oxygen, (c) air issuing into ethylene, and (d) oxygen issuing into diluted ethylene. (Z(sub st)) is 0.064 in flames (a) and (c) and 0.78 in flames (b) and (d). The convection direction is from fuel to oxidizer in flames (a) and (b) and from oxidizer to fuel in flames (c) and (d). Under the assumption of equal diffusivities of all species and heat, the stoichiometric contours of these flames have identical temperatures and nitrogen concentrations.
Sunderland, P. B.↗
Numerical simulation of the interaction of transport, diffusion and chemical reactions in an urban plume
A model system is presented that takes into account the main physical and chemical processes on the regional scale here in an area of 100x100 sq km. The horizontal gridsize used is 2x2 sq km. For a case study, it is demonstrated how the model system can be used to separate the contributions of the processes advection, turbulent diffusion, and chemical reactions to the diurnal cycle of ozone. In this way, typical features which are visible in observations and are reproduced by the numerical simulations can be interpreted.
Vogel, Bernhard↗
Solutions for Reacting and Nonreacting Viscous Shock Layers with Multicomponent Diffusion and Mass Injection
Numerical solutions are presented for the viscous shocklayer equations where the chemistry is treated as being either frozen, equilibrium, or nonequilibrium. Also the effects of the diffusion model, surface catalyticity, and mass injection on surface transport and flow parameters are considered. The equilibrium calculations for air species using multicomponent: diffusion provide solutions previously unavailable. The viscous shock-layer equations are solved by using an implicit finite-difference scheme. The flow is treated as a mixture of inert and thermally perfect species. Also the flow is assumed to be in vibrational equilibrium. All calculations are for a 45 deg hyperboloid. The flight conditions are those for various altitudes and velocities in the earth's atmosphere. Data are presented showing the effects of the chemical models; diffusion models; surface catalyticity; and mass injection of air, water, and ablation products on heat transfer; skin friction; shock stand-off distance; wall pressure distribution; and tangential velocity, temperature, and species profiles.
Moss, J. N.↗
Reacting viscous-shock-layer solutions with multicomponent diffusion and mass injection
Numerical solutions of the viscous-shock-layer equation where the chemistry is treated as being either frozen, equilibrium, or nonequilibrium are presented. Also the effects of the diffusion model, surface catalysis, and mass injection on surface transport and flow parameters are considered. The flow is treated as a mixture of five inert and thermally perfect species. The viscous-shock-layer equations are solved by using an implicit-difference scheme. All calculations are for hyperboloids with included angles of 20 and 45. The flight conditions are those for various altitudes and velocities in the earth's atmosphere. Data are presented to show the effects of the chemical models; diffusion models; surface catalysis; and mass injection of air on heat transfer; skin friction; shock standoff distance; wall pressure distribution; and tangential velocity, temperature, and species profiles. The results show that an equilibrium analysis can substantially overpredict the heat-transfer rates for flow conditions experienced by earth-orbital entry vehicles. Moreover, at such conditions surface catalysis significantly influences heat-transfer and flow-field properties. If a binary rather than a multicomponent diffusion model is assumed, negligible errors in most flow properties result. Quantitative results are presented that show the effect of mass injection on flow properties within and downstream of the injection region.
Moss, J. N.↗
A finite volume Euler calculation of the aerodynamics of transonic airfoil-vortex interaction
Unsteady inviscid transonic airfoil-vortex interaction is numerically analyzed by solving the two-dimensional unsteady Euler equations in integral form using a finite volume scheme. The solution procedure is based on an explicit Runge-Kutta time-stepping scheme wherein the spatial terms are central-differenced and a combination of second- and fourth-differences in the flow variables is used to form the numerical dissipation terms to stabilize the scheme. A velocity decomposition technique is applied to alleviate the problem of vortex diffusion by the numerical dissipation terms and to treat the interaction of a Rankine vortex with an airfoil accurately. Results obtained are compared with available numerical data.
Damodaran, Murali↗
Modeling Microwave-Enhanced Chemical Vapor Infiltration Process for Preventing Premature Pore Closure
The chemical vapor infiltration (CVI) process involves infiltrating a porous preform with reacting gases that undergo chemical transformation at high temperatures to deposit the ceramic phase within the pores, ultimately leading to a dense composite. The conventional CVI process in composite manufacturing needs to follow an isothermal approach to minimize temperature differences between the external and internal surfaces of the preform, ensuring that reactive gases infiltrate internal pores before external surfaces seal. Here, this study addresses the challenge of premature pore closure in CVI processes through microwave heating. A frequency-domain microwave solver is developed in OpenFOAM to investigate volumetric heating mechanisms within the preform. Through numerical studies, we demonstrate the capability of microwave heating of creating an inside-out temperature inversion. This inversion accelerates reactions proximal to the preform center, effectively mitigating the risk of premature external pore closure and ensuring uniform densification. The results reveal a significant enhancement in temperature inversion when high-permittivity reflectors are incorporated to generate resonant waves. This microwave heating strategy is then coupled with high-fidelity direct numerical simulation (DNS) of reacting flow, enabling the analysis of resulting densification processes. The DNS includes detailed chemistry and realistic diffusion coefficients. The numerical results can be used to estimate the impact of microwave-induced temperature inversion on densification in productions.
42 ENGINEERING↗
Modeling Microwave-Enhanced Chemical Vapor Infiltration Process for Preventing Premature Pore Closure
The chemical vapor infiltration (CVI) process involves infiltrating a porous preform with reacting gases that undergo chemical transformation at high temperatures to deposit the ceramic phase within the pores, ultimately leading to a dense composite. The conventional CVI process in composite manufacturing needs to follow an isothermal approach to minimize temperature differences between the external and internal surfaces of the preform, ensuring that reactive gases infiltrate internal pores before external surfaces seal. This study addresses the challenge of premature pore closure in CVI processes through microwave heating. A frequency-domain microwave solver is developed in Open-FOAM to investigate volumetric heating mechanisms within the preform. Through numerical studies, we demonstrate the capability of microwave heating of creating an inside-out temperature inversion. This inversion accelerates reactions proximal to the preform center, effectively mitigating the risk of premature external pore closure and ensuring uniform densification. The results reveal a significant enhancement in temperature inversion when high-permittivity reflectors are incorporated to generate resonant waves. This microwave heating strategy is then coupled with high-fidelity direct numerical simulation (DNS) of reacting flow, enabling the analysis of resulting densification processes. The DNS simulation includes detailed chemistry and realistic diffusion coefficients. The numerical results can be used to estimate the impact of microwave-induced temperature inversion on densification in productions.
Ge, Wenjun↗
Analysis and design of numerical schemes for gas dynamics. 2: Artificial diffusion and discrete shock structure
The effect of artificial diffusion on discrete shock structures is examined for a family of schemes which includes scalar diffusion, convective upwind and split pressure (CUSP) schemes, and upwind schemes with characteristics splitting. The analysis leads to conditions on the diffusive flux such that stationary discrete shocks can contain a single interior point. The simplest formulation which meets these conditions is a CUSP scheme in which the coefficients of the pressure differences is fully determined by the coefficient of convective diffusion. It is also shown how both the characteristic and CUSP schemes can be modified to preserve constant stagnation enthalpy in steady flow, leading to four variants, the E and H-characteristic schemes, and the E and H-CUSP schemes. Numerical results are presented which confirm the properties of these schemes.
Jameson, Antony↗
Structure and Stability of Burke-Schumann Diffusion Flames
The general goal of this NASA Grant is twofold: to improve our understanding of (1) the influence of buoyancy on the stability and structure of Burke-Schumann type diffusion flames, and (2) the effects of buoyancy on vortex-flame interactions in co-flow diffusion flames. A numerical code with a higher order accuracy for spatial discretization is developed in this project for simulation of time-dependent diffusion flames by Sheu and Sheu and Chen, and an extended reduced mechanism is incorporated for prediction of methane oxidation and NO(x)(NO, NO2, and N2O) formation and emission from methane Burke-Schumann diffusion flame (BSDF) as reported in Sheu, and Sheu and Chen. Initial investigation of vortex and flame interaction within the context of fast chemistry is reported. Experiments are conducted in reduced pressure to study the lift-off and stabilization of methane-fueled BSDF in reduced buoyancy environments due to reduced pressure. Measurements of temperature and species concentrations are made in normal and reduced pressure environments to study the effects of buoyancy on the structure of BSDF, and will be reported in this paper. To study the buoyancy effects on the lift-off and stabilization of methane-fueled jet diffusion flames in coflowing air, a glovebox investigation, Enclosed Laminar Flames (ELF), has been proposed and approved for space-based testing on the fourth United States Microgravity Payload (USMP-4) mission, scheduled for October 1997. A brief description of the ELF investigation is also presented.
Lee, Yong G.↗
Simple Analytical Expressions for Electron Pitch Angle Diffusion Coefficients
In the present paper, we have calculated electron pitch angle diffusion coefficients due to resonant interactions with whistler mode lower band chorus (LBC), upper band chorus (UBC), and electrostatic electron cyclotron harmonic (ECH) waves. Calculations have been performed at two values of the ratio of electron plasma frequency to gyro-frequency and thirteen representative values of electron energies for the plasma sheet electrons. The numerical data of diffusion coefficients have been fitted to simple analytical expressions for each wave mode. These analytical expressions allow for simple evaluation of pitch angle diffusion coefficients for the arbitrary pitch angle, energy, the ambient magnetic field, the wave amplitude, and the ratio of plasma frequency to gyro-frequency. In the case of LBC waves, the analytical coefficients are generally within a factor of two to the numerical coefficients, except at higher pitch angles where the numerical coefficients drop to show negligible values. Likewise, also for UBC waves, the analytical coefficients are generally within a factor of two to the numerical coefficients. The analytical coefficients for ECH waves are generally in agreement with numerical coefficients. However, the analytical expressions developed for ECH waves do not reproduce sharp fluctuations, dips, and gaps, the so-called banded structure observed in the data of numerical coefficients. The effect of the shape variation on the profile of the wave spectral intensity and the effect of cold temperature on the numerical coefficients are also investigated. Applications of the analytical expressions of diffusion coefficients are discussed.
coefficients↗
Investigation of core transport changes in DIII-D discharges with off-axis T e profile peaks
Abstract DIII-D discharges that transition to H-mode solely with off-axis electron cyclotron heating (ECH) often exhibit strong off-axis peaking of electron temperature profiles at the heating location. Electron heat transport properties near these off-axis temperature peaks have been studied using modulated ECH. The Fourier analyzed electron temperature data have been used to infer electron thermal diffusivity. Comparisons with numerical solutions of the time-dependent electron thermal equation find that the data are consistent with a narrow region with electron diffusivity χ e an order of magnitude lower than the average value across the plasma, suggesting an electron internal transport barrier (ITB) near the ECH heating location. Detailed profile analysis and equilibrium reconstructions suggest that the formation of these ITBs are correlated with off-axis values of the safety factor q being near 1. Furthermore, the ECH driven H-mode discharges demonstrate more rapid electron heating rate near the ECH deposition location than L-mode discharges with higher auxiliary ECH heating power. Additional modeling attributes this difference to the modification of electron heat transport in the core at the L-H transition, which also sustains the off-axis electron temperature peaks.
Physics↗
Significance of Low‐Velocity Zones on Solute Retention in Rough Fractures
Natural fractures are characterized by high internal heterogeneity. This internal variability is the cause of flow channeling, which in turn leads to contaminant transport taking place primarily along the high-velocity channels. Mass exchange between the high-velocity channels and the low-velocity zones has the potential to enhance contaminant retention, due to solute diffusion into the low-velocity zones and subsequent exposure to additional surface area for diffusion into the bordering rock matrix. Here, we derive a random walk particle tracking method for heterogeneous fractures, which includes an additional term to account for the aperture gradient. The method takes into account advection, diffusion in the fracture and matrix diffusion. The developed numerical framework is applied to assess the effect of low-velocity zones in rough self-affine fractures. The results show that diffusion into low-velocity zones has a visible but modest impact on contaminant retention. The magnitude of this impact does not change considerably, regardless of whether diffusion into the rock matrix is considered in the model, and increases for a decreasing average Péclet number of the fracture.
58 GEOSCIENCES↗
Numerical thermalization in 2D PIC simulations: Practical estimates for low-temperature plasma simulations
The process of numerical thermalization in particle-in-cell (PIC) simulations has been studied extensively. It is analogous to Coulomb collisions in real plasmas, causing particle velocity distributions (VDFs) to evolve toward a Maxwellian as macroparticles experience polarization drag and resonantly interact with the fluctuation spectrum. This paper presents a practical tutorial on the effects of numerical thermalization in 2D PIC applications. Scenarios of interest include simulations, which must be run for many thousands of plasma periods and contain a population of cold electrons that leave the simulation space very slowly. This is particularly relevant to many low-temperature plasma discharges and materials processing applications. We present numerical drag and diffusion coefficients and their associated timescales for a variety of grid resolutions, discussing the circumstances under which the electron VDF is modified by numerical thermalization. Though the effects described here have been known for many decades, direct comparison of analytically derived, velocity-dependent numerical relaxation timescales to those of other relevant processes has not often been applied in practice due to complications that arise in calculating thermalization rates in 1D simulations. Using these comparisons, we estimate the impact of numerical thermalization in several examples of low-temperature plasma applications including capacitively coupled plasma discharges, inductively coupled plasma discharges, beam plasmas, and hollow cathode discharges. Finally, we discuss possible strategies for mitigating numerical relaxation effects in 2D PIC simulations.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
A high-order, localized-artificial-diffusivity method for Eulerian simulation of multi-material elastic-plastic deformation with strain hardening
A high-order method for Eulerian simulation of material undergoing large elastic–plastic deformation is developed. Thermodynamically consistent hyperelastic constitutive relations are assumed, facilitating the treatment of solids, liquids, and gases in a unified manner. Here, the method enables the simulation of multi-material interactions using a diffuse interface approach. Numerical capturing of material interfaces, shock waves, contact surfaces, and elastic-plastic strain discontinuities using high-order compact-difference schemes is assisted by Localized Artificial Diffusivity (LAD). In the new setting involving elastic–plastic deformation, the previously established terms for the artificial properties are verified to effectively regularize normal shocks. Additional LAD terms are introduced to the elastic and plastic kinematic equations to regularize shear shocks and other strain discontinuities, improving solution stability. Other important features of the method that improve robustness include the numerical treatment of compatibility terms in the kinematic equations, and the treatment of rotation. Particular emphasis is focused toward new advancements of the methods for plastic-deformation integration and the associated strain hardening of the material, including rate-dependent plasticity. The method is demonstrated on a variety of test problems, including 1-D impacts, a variant of the Shu-Osher problem, a Taylor impact, and a Richtmyer-Meshkov instability between two elastic–plastic solids with strain hardening.
42 ENGINEERING↗
Evidence for the 3D Radiative Effects of Boundary-Layer Clouds From Observations of Direct and Diffuse Surface Solar Fluxes
Numerical experiments have revealed the importance of horizontal transport of light in the presence of clouds (“3D effects”), with consequences for climate, weather, and solar resource availability predictions. Yet, analysis of 3D effects from observations remain sparse because of the difficulty to isolate the effect of horizontal transport in radiation measurements. In this study, we provide observational evidence for 3D effects based on the direct-diffuse partition of surface solar fluxes. It is compared to outputs from the ecRad radiative transfer scheme run on retrieved cloud profiles. Furthermore, the direct-beam calculation takes careful account of the field-of-view of the pyrheliometer to ensure consistency between observed and modeled direct fluxes. Only the solver that accounts for 3D effects is able to reproduce the observed mean direct-diffuse partition as a function of solar zenith angle and cloud cover, in particular at large solar zenith angles where cloud sides intercept most of the direct beam.
54 ENVIRONMENTAL SCIENCES↗
Numerical Simulation of an Enclosed Laminar Jet Diffusion Flame in Microgravity Environment: Comparison with ELF Data
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. The diffusion flame is found at the surface where the fuel jet and oxygen meet, react, and consume each other. 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. Flame stability is strongly dependent on the fuel jet velocity. When the fuel jet velocity is sufficiently low, the diffusion flame anchors at the burner rim. When the fuel jet velocity is increased, the flame base gradually moves downstream. However, when the fuel jet velocity increases beyond a critical value, the flame base abruptly jumps downstream. When this "jump" occurs, the flame is said to have reached its lift-off condition and the critical fuel jet velocity is called the lift-off velocity. While lifted, the flame is not attached to the burner and it appears to float in mid-air. Flow conditions are such that the flame cannot be maintained at the burner rim despite the presence of both fuel and oxygen. When the fuel jet velocity is further increased, the flame will eventually extinguish at its blowout condition. In contrast, if the fuel jet velocity of a lifted flame is reduced, the flame base moves upstream and abruptly returns to anchor at the burner rim. The fuel jet velocity at reattachment can be much lower than that at lift off, illustrating the hysteresis effect present in flame stability. Although there have been numerous studies of flame stability, the controlling mechanisms are not well understood. This uncertainty is described by Pitts in his review of various competing theories of lift off and blow out in turbulent jet diffusion flames. There has been some research on the stability of laminar flames, but most studies have focused on turbulent flames. It is also well known that the airflow around the fuel jet can significantly alter the lift off, reattachment and blow out of the jet diffusion flame. Buoyant convection is sufficiently strong in 1-g flames that it can dominate the flow-field, even at the burner rim. In normal-gravity testing, it is very difficult to delineate the effects of the forced airflow from those of the buoyancy-induced flow. Comparison of normal-gravity and microgravity flames provides clear indication of the influence of forced and buoyant flows on the flame stability. The overall goal of the Enclosed Laminar Flames (ELF) investigation (STS-87/USMP-4 Space Shuttle mission, November to December 1997) is to improve our understanding of the effects of buoyant convection on the structure and stability of co-flow diffusion flame, e.g., see http://zeta.lerc.nasa.gov/expr/elf.htm. The ELF hardware meets the experiment hardware limit of the 35-liter interior volume of the glovebox working area, and the 180x220-mm dimensions of the main door. The ELF experiment module is a miniature, fan-driven wind tunnel, equipped with a gas supply system. A 1.5-mm diameter nozzle is located on the duct's flow axis. The cross section of the duct is nominally a 76-mm square with rounded corners. The forced air velocity can be varied from about 0.2 to 0.9 m/s. The fuel flow can be set as high as 3 std. cubic centimeter (cc) per second, which corresponds to a nozzle exit velocity of up to 1.70 m/s. The ELF hardware and experimental procedure are discussed in detail in Brooker et al. The 1-g test results are repeated in several experiments following the STS-87 Mission. The ELF study is also relevant to practical systems because the momentum-dominated behavior of turbulent flames can be achieved in laminar flames in microgravity. The specific objectives of this paper are to evaluate the use reduced model for simulation of flame lift-off and blowout.
Jia, Kezhong↗