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,531 records · Page 85

All sources NASA NTRS DOE OSTI

NASA NTRS1974

Nitric oxide formation by meteoroids in the upper atmosphere

The process of nitric oxide formation during atmospheric entry of meteoroids is analyzed theoretically. An ablating meteoroid is assumed to be a point source in a uniform flow with a continuum regime evolving in its wake. The amount of nitric oxide produced by high-temperature reactions of air in the continuum regime is calculated by numerical integration of chemical-rate equations. This is accomplished by assuming that flow properties are constant across the reacting region, its radius being determined from considerations of shock-wave formation and molecular diffusion. The results, when summed over the observed mass, velocity, and entry-angle distributions of meteoroids, provide annual global production rates of nitric oxide as a function of altitude. The peak production is found to occur between 90- and 100-km altitude, the total annual rate being around 40 million kg.

Menees, G. P.↗

NASA NTRS1975

Further numerical calculations of the circulation of the atmosphere of Venus

Results of two-dimensional simulations of the deep circulation of Venus are presented. They prove that the high surface temperature can only be explained by the greenhouse effect, and that Goody and Robinson's (1966) dynamical model is not valid. Analytical radiative equilibrium solutions for a semi-grey atmosphere, both with and without an internal heat source, are presented. It is shown that the greenhouse effect is sufficient to produce the high surface temperature if the thermal optical depth is much greater than 100 and the ratio of solar to thermal optical depths does not exceed 0.005. A two-dimensional version of a three-dimensional model is used to test the validity of the new mechanism proposed by Gierasch (1975) to explain the 4-day circulation. Numerical experiments with large horizontal viscosities failed to show strong zonal velocities even for the case of large Prandtl numbers. It is observed that the dissipation of angular momentum introduced by the strong horizontal diffusion more than compensates for the upward transport of angular momentum due to the Hadley cell.

De Rivas, E. K.↗

NASA NTRS1978

A theoretical study of a laminar diffusion flame

Theoretical models of an axisymmetric laminar diffusion flame are discussed, with an emphasis on the behavior of such flames at increasing pressures. The flame-sheet or Burke-Schumann model (in terms of Bessel functions) and various boundary layer numerical solutions are presented and their results compared with experimental data. The most promising theoretical model combines the numerical flow field solution of the Patankar-Spalding computer code with the Pratt-Wormeck chemical reaction subroutine. The flame shapes for pressures of 1, 5, 10, 20, and 50 atmospheres were computed and agree remarkably well with experimental data. There is a noticeable shape change with pressure, believed to be a result of buoyancy effects. The chemical concentration profiles do not exhibit much dependence on pressure, a reflection of the fact that only one chemical mechanism was utilized at all pressures.

Frair, K. L.↗

NASA NTRS1984

A generalized behavioral model for rotating short period comets with spectral orbital elements and axial orientation

A generalized model for short period comets is developed which integrates in a fairly rigorous manner the isolation history of regions on rotating comets with specified axial orientation and the complex feedback processes involving heat, gas and dust transport, dust mantle development and coma opacity. Attention is focused on development, reconfiguration and partial or complete launching of dust mantles and the reciprocal effects of these three processes on ice surface temperature and gas and dust production. The dust mantle controls the H2O flux not only by its effect on the temperature at the ice interface but (dominantly) by its dynamic stability which strongly influences vapor diffusivity. The model includes the effects of latitude, rotation and spin axis orientation are included and applied to an initially homogeneous sphere of H2O ice and silicate using the orbital parameters of comet Encke. Numerous variations of the model, using combinations of grain size distribution, dust-to-ice ratio, latitude and spin axis orientation, are presented and discussed. Resulted for a similar nonrotating, constant Sun orientation models are also included.

Fanale, F. P.↗

NASA NTRS1987

Approximation techniques for domain identification in two-dimensional parabolic systems under boundary observations

The identification of the geometrical structure of the system boundary for a two-dimensional diffusion system is considered. The domain identification treated is converted into an optimization problem based on a fit-to-data criterion and theoretical convergence results for approximate identification techniques are discussed. Results of numerical experiments to demonstrate the efficiency of the theoretical ideas are reported.

Banks, H. T.↗

NASA NTRS1993

Laser diagnostics for microgravity droplet studies

Rapid advances have recently been made in numerical simulation of droplet combustion under microgravity conditions, while experimental capabilities remain relatively primitive. Calculations can now provide detailed information on mass and energy transport, complex gas-phase chemistry, multi-component molecular diffusion, surface evaporation and heterogeneous reaction, which provides a clearer picture of both quasi-steady as well as dynamic behavior of droplet combustion. Experiments concerning these phenomena typically result in pictures of the burning droplets, and the data therefrom describe droplet surface regression along with flame and soot shell position. With much more precise, detailed, experimental diagnostics, significant gains could be made on the dynamics and flame structural changes which occur during droplet combustion. Since microgravity experiments become increasingly more expensive as they progress from drop towers and flights to spaceborne experiments, there is a great need to maximize the information content from these experiments. Sophisticated measurements using laser diagnostics on individual droplets and combustion phenomena are now possible. These include measuring flow patterns and temperature fields within droplets, vaporization rates and vaporization enhancement, radical species profiling in flames and gas-phase flow-tagging velocimetry. Although these measurements are sophisticated, they have undergone maturation to the degree where with some development, they are applicable to studies of microgravity droplet combustion. This program beginning in September of 1992, will include a series of measurements in the NASA Learjet, KC-135 and Drop Tower facilities for investigating the range of applicability of these diagnostics while generating and providing fundamental data to ongoing NASA research programs in this area. This program is being conducted in collaboration with other microgravity investigators and is aimed toward supplementing their experimental efforts.

Winter, Michael↗

NASA NTRS1996

Ignition, Burning, and Extinction of a Strained Fuel Strip

Flame structure and ignition and extinction processes associated with a strained fuel strip are explored numerically using detailed transport and complex kinetics for a propane-air reaction. Ignition modes are identified that are similar to those predicted by one-step activation energy asymptotics, i.e., modes in which diffusion flames can ignite as independent or dependent interfaces, and modes in which single premixed or partially premixed flames ignite. These ignition modes are found to be dependent on critical combinations of strain rate, fuel strip thickness, and initial reactant temperatures. Extinction in this configuration is seen to occur due to fuel consumption by adjacent flames, although viscosity is seen to have the effect of delaying extinction by reducing the effective strain rate and velocity field experienced by the flames.

Selerland, T.↗

NASA NTRS1998

Ignition, Burning and Extinction of a Strained Fuel Strip with Complex Kinetics

Flame structure and ignition and extinction processes associated with a strained fuel strip are explored numerically using detailed transport and complex kinetics for a propane-air reaction. Ignition modes are identified that are similar to those predicted by one-step activation energy asymptotics, i.e., modes in which diffusion flames can ignite as independent or dependent interfaces and modes in which single premixed or partially premixed flames ignite and burn. These ignition modes have been found to be dependent on critical combinations of strain rate, fuel strip thickness and initial reactant temperatures. The formation of NO/NO2 is found to be strongly dependent on strain rate and the local molecular mixing of reactants which occurs as a consequence of strain. Extinction in this configuration is seen to occur due to fuel consumption by adjacent flames, although viscosity is seen to have the effect of delaying extinction by reducing the effective strain rate experienced by the flames. Response of the flames to oscillatory strain rates is seen to be strongly dependent on the amplitude and frequency of the oscillation.

Selerland, T.↗

NASA NTRS2001

Electro-Thermal Simulation Studies of SiC Junction Diodes Containing Screw Dislocations Under High Reverse-Bias Operation

The objective of this work was to conduct a modeling study of SiC P-N junction diodes operating under high reverse biased conditions. Analytical models and numerical simulation capabilities were to be developed for self-consistent electro-thermal analysis of the diode current-voltage (I-V) characteristics. Data from GRC indicate that screw dislocations are unavoidable in large area SiC devices, and lead to changes in the SiC diode electrical response characteristics under high field conditions. For example, device instability and failures linked to internal current filamentation have been observed. The physical origin of these processes is not well understood, and quantitative projections of the electrical behavior under high field and temperature conditions are lacking. Thermal calculations for SiC devices have not been reported in the literature either. So estimates or projections of peak device temperatures and power limitations do not exist. This numerical study and simulation analysis was aimed at resolving some of the above issues. The following tasks were successfully accomplished: (1) Development of physically based models using one- and two-dimensional drift-diffusion theory for the transport behavior and I-V characteristics; (2) One- and two-dimensional heat flow to account for internal device heating. This led to calculations of the internal temperature profiles, which in turn, were used to update the electrical transport parameters for a self-consistent analysis. The temperature profiles and the peak values were thus obtainable for a given device operating condition; (3) Inclusion of traps assumed to model the presence of internal screw dislocations running along the longitudinal direction; (4) Predictions of the operating characteristics with and without heating as a function of applied bias with and without traps. Both one and two-dimensional cases were implemented; (5) Assessment of device stability based on the operating characteristics. The presence of internal non-uniformities, particularly filamentary structures, was probed and demonstrated; (6) Cause and physical origins of filamentary behavior and unstable I-V characteristics were made transparent; (7) It was demonstrated that diodes containing defects would be more prone to thermal breakdown associated with the temperature dependent decrease in the thermal conductivity; and (8) Finally, negative differential resistance (S-shaped NDR) which can potential lead to device instability and filamentary behavior was shown to occur for diodes containing a line of defects such as could be associated with a screw dislocation line.

Joshi, R. P.↗

NASA NTRS2003

Cosmic: Carbon Monoxide And Soot In Microgravity Inverse Combustion

Almost seventy percent of fire related deaths are caused by the inhalation of toxins such as CO and soot that are produced when fires become underventilated.(1) Although studies have established the importance of CO formation during underventilated burning,(2) the formation processes of CO (and soot) in underventilated fires are not well understood. The goal of the COSMIC project is to study the formation processes of CO and soot in underventilated flames. A potential way to study CO and soot production in underventilated flames is the use of inverse diffusion flames (IDFs). An IDF forms between a central air jet and a surrounding fuel jet. IDFs are related to underventilated flames because they may allow CO and soot to escape unoxidized. Experiments and numerical simulations of laminar IDFs of CH4 and C2H4 were conducted in 1-g and micro-g to study CO and soot formation. Laminar flames were studied because turbulent models of underventilated fires are uncertain. Microgravity was used to alter CO and soot pathways. A IDF literature survey, providing background and establishing motivation for this research, was presented at the 5th IWMC.(3) Experimental results from 1-g C2H4 IDFs and comparisons with simulations, demonstrating similarities between IDFs and underventilated fires, were presented at the 6th IWMC.(4) This paper will present experimental results from micro-g and 1-g IDFs of CH4 and C2H4 as well as comparisons with simulations, further supporting the relation between IDFs and underventilated flames.

Mikofski, M. A.↗

NASA NTRS2014

Sound Transmission Loss Through a Corrugated-Core Sandwich Panel with Integrated Acoustic Resonators

The goal of this study is to better understand the effect of structurally integrated resonators on the transmission loss of a sandwich panel. The sandwich panel has facesheets over a corrugated core, which creates long aligned chambers that run parallel to the facesheets. When ports are introduced through the facesheet, the long chambers within the core can be used as low-frequency acoustic resonators. By integrating the resonators within the structure they contribute to the static load bearing capability of the panel while also attenuating noise. An analytical model of a panel with embedded resonators is derived and compared with numerical simulations. Predictions show that acoustic resonators can significantly improve the transmission loss of the sandwich panel around the natural frequency of the resonators. In one configuration with 0.813 m long internal chambers, the diffuse field transmission loss is improved by more than 22 dB around 104 Hz. The benefit is achieved with no added mass or volume relative to the baseline structure. The embedded resonators are effective because they radiate sound out-of-phase with the structure. This results in destructive interference, which leads to less transmitted sound power.

Schiller, Noah H.↗

NASA NTRS2016

The Impacts of Dry Dynamic Cores on Asymmetric Hurricane Intensification

The fundamental pathways for tropical cyclone (TC) intensification are explored by considering axisymmetric and asymmetric impulsive thermal perturbations to balanced, TC-like vortices using the dynamic cores of three different nonlinear numerical models. Attempts at reproducing the results of previous work, which used the community WRF Model, revealed a discrepancy with the impacts of purely asymmetric thermal forcing. The current study finds that thermal asymmetries can have an important, largely positive role on the vortex intensification, whereas other studies find that asymmetric impacts are negligible. Analysis of the spectral energetics of each numerical model indicates that the vortex response to asymmetric thermal perturbations is significantly damped in WRF relative to the other models. Spectral kinetic energy budgets show that this anomalous damping is primarily due to the increased removal of kinetic energy from the vertical divergence of the vertical pressure flux, which is related to the flux of inertia-gravity wave energy. The increased kinetic energy in the other two models is shown to originate around the scales of the heating and propagate upscale with time from nonlinear effects. For very large thermal amplitudes (50 K), the anomalous removal of kinetic energy due to inertia-gravity wave activity is much smaller, resulting in good agreement between models. The results of this paper indicate that the numerical treatment of small-scale processes that project strongly onto inertia-gravity wave energy can lead to significant differences in asymmetric TC intensification. Sensitivity tests with different time integration schemes suggest that diffusion entering into the implicit solution procedure is partly responsible for the anomalous damping of energy.

Guimond, Stephen R.↗

NASA NTRSDate not supplied

Large-Eddy Simulations of a Single-Injector Cooling Flow Using the High-Order Flux Reconstruction Method

A single-injector cooling flow into a heated crossflow was used as a validation case for large eddy simulations (LES) from two high-order CFD codes with comparisons to a state-of-the-art Reynolds averaged Navier-Stokes (RANS) turbulence model. The focus of this paper is on the validation of LES for one of the high-order CFD codes, namely the GFR (Glenn Flux Reconstruction) code that is being developed at NASA Glenn Research Center. Fourth-order LES were performed for the single-injector cooling flow configuration at blowing ratios of both 1.0 and 2.0, and a fifth-order LES was also performed with a blowing ratio of 1.0. The GFR simulations agree very well with the experiment mean quantities and reasonably well with the experiment turbulence quantities. The LES solutions from GFR and the other high-order CFD code, FDL3DI, agreed very closely for nearly every flow quantity examined, despite the fact that these two codes use completely different numerical methods, different grid geometries and domains, and different inflow boundary conditions. Finally, both LES show a significant accuracy improvement over RANS turbulence models for flows with large temperature gradients, where the standard gradient diffusion approximation is insufficient for temperature predictions.

Large Eddy Simulations↗

NASA NTRSDate not supplied

Large Eddy Simulations of a Single-Injector Cooling Flow Using the High-Order Flux Reconstruction Method

Large Eddy Simulation↗

DOE OSTI2021

On the Equilibrium Electrostatic Potential and Light-Induced Charge Redistribution in Halide Perovskite Structures

Lead halide perovskites are semiconductor materials which are employed as nonintentionally doped absorbers inserted between two selective carrier transport layers (SCTL), realizing a p-i-n or n-i-p heterojunction. In our study, we have developed and investigated a lateral device, based on methylammonium lead iodide (MAPbI3) in which the p-i-n heterojunction develops in the horizontal direction. Our research suggests that the effective doping level in the MAPbI3 film should be very low, below 1012 cm-3. Along the vertical direction, this doping level is not enough to screen the electric field of the buried heterojunction with the SCTL. The perovskite work function is therefore affected by the work function of the SCTL underneath. From drift-diffusion simulations, we show that intrinsic perovskite-SCTL structures develop mV range surface photovoltages (SPVs) under continuous illumination. However, perovskite-SCTL structures can develop SPVs of hundreds of mV, as confirmed by our measurements. We therefore analyzed the compatibility between low doping and low defect densities in the perovskite layer and such high SPV values using numerical modeling. It is shown that these high SPV values could originate from electronic processes due to large band offsets in the buried perovskite-SCTL heterojunctions, or at the SCTL-transparent conductive oxide (TCO) buried heterojunction. However, such electronic processes can hardly explain the long SPV persistence after switching off the illumination.

defects↗

NASA NTRS1984

The dayside Venus ionosphere. II - Combined numerical model of ion and neutral composition above 120 km

The median vertical profile of ion densities for the dayside Venusian ionosphere obtained by the orbiter retarding potential analyzer (ORPA) is simulated by one-dimensional model calculations. The model includes both neutral and ionic chemistry, eddy and molecular diffusion for neutral constituents, and ion-plasma diffusion for ionic constituents. The electron and ion temperatures measured by the ORPA are used to calculate the plasma diffusion coefficients and scale heights for ions. The predicted O2(+) densities below about 200 km agree particularly well with observations by the ORPA, but the model values are significantly less than those measured by the orbiter ion mass spectrometer. The observed ion composition is interpreted in terms of densities of the neutral atmosphere and its composition.

Shimazaki, T.↗

NASA NTRS1999

Experimental Observations on a Low Strain Counter-Flow Diffusion Flame: Flow and Bouyancy Effects

Diffusion flames are of great interest in fire safety and many industrial processes. The counter-flow configuration provides a constant strain flow, and therefore is ideal to study the structure of diffusion flames. Most studies have concentrated on the high velocity, high strain limit, since buoyantly induced instabilities will disintegrate the planar flame as the velocity decreases. Only recently, experimental studies in microgravity conditions have begun to explore the low strain regimes. Numerical work has shown the coupling between gas phase reaction rates, soot reaction rates, and radiation. For these programs, size, geometry and experimental conditions have been chosen to keep the flame unaffected by the physical boundaries. When the physical boundaries can not be considered infinitely far from the reaction zone discrepancies arise. A computational study that includes boundary effects and accounts for the deviations occurring when the major potential flow assumptions are relaxed was presented by Borlik et al. This development properly incorporates all heat loss terms and shows the possibility of extinction in the low strain regime. A major constraint of studying the low strain regime is buoyancy. Buoyant instabilities have been shown to have a significant effect on the nature of reactants and heat transport, and can introduce instabilities on the flow that result in phenomena such as flickering or fingering. The counter-flow configuration has been shown to provide a flame with no symmetry disrupting instabilities for inlet velocities greater than 50 mm/s. As the velocity approaches this limit, the characteristic length of the experiment has to be reduced to a few millimetres so as to keep the Rayleigh number (Ra L = (β g 0 L 3 ∆T)/(𝑎v)) below 2000. In this work, a rectangular counter-flow burner was used to study a two-dimensional counter-flow diffusion flame. Flow visualisation and Particle Image Velocimetry served to describe the nature of the stagnation plane for strain rates smaller than 100 (1/s). These experiments were conducted with a non-reacting flow. Video images of a propane air diffusion flame were used to describe the behaviour of a diffusion flame in this regime. Flame geometry and pulsation frequency are described.

J. A. Sutula↗

NASA NTRS2016

Tracking Blade Tip Vortices for Numerical Flow Simulations of Hovering Rotorcraft

Blade tip vortices generated by a helicopter rotor blade are a major source of rotor noise and airframe vibration. This occurs when a vortex passes closely by, and interacts with, a rotor blade. The accurate prediction of Blade Vortex Interaction (BVI) continues to be a challenge for Computational Fluid Dynamics (CFD). Though considerable research has been devoted to BVI noise reduction and experimental techniques for measuring the blade tip vortices in a wind tunnel, there are only a handful of post-processing tools available for extracting vortex core lines from CFD simulation data. In order to calculate the vortex core radius, most of these tools require the user to manually select a vortex core to perform the calculation. Furthermore, none of them provide the capability to track the growth of a vortex core, which is a measure of how quickly the vortex diffuses over time. This paper introduces an automated approach for tracking the core growth of a blade tip vortex from CFD simulations of rotorcraft in hover. The proposed approach offers an effective method for the quantification and visualization of blade tip vortices in helicopter rotor wakes. Keywords: vortex core, feature extraction, CFD, numerical flow visualization

Numerical↗