Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion equation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,603 records · Page 89

Theoretical model of vertical distributions of CO and CH4 in the mesosphere and upper stratosphere.

Carbon compounds (CO and CO2) and methane are added to an existing computer program (Shimazaki and Laird, 1970, 1972), which calculates diurnal variations in density distributions of pure oxygen compounds, hydrogen compounds, and nitrogen compounds. For the present model calculations, the original program was modified to calculate the steady-state profiles for all constituents, and then these solutions were used as the initial condition for solving the time-dependent equations. The effects of chemical reactions and vertical eddy diffusion transport on the CH4 and CO concentrations above 40 km, where one-dimensional treatment may not be too inappropriate, are discussed. Some of the 15 constituents included in the computer program interact with CO and CH4 chemically, and CH4 is more sensitive to the eddy diffusion coefficient than any other constituent. Therefore it is easier to determine the most appropriate profile of the eddy diffusion coefficient by comparing the model profile with the observations for CH4.

Shimazaki, T.↗

Electron precipitation in the postmidnight sector of the auroral zones

Measurements of the angular distributions and energy spectra of electron intensities within the energy range 50 eV to 15 keV with electrostatic analyzer arrays on board the low-altitude satellite Injun 5 are reported for the postmidnight sector of the auroral zones during the high-intensity events accompanying magnetic substorms. Precipitation features on closed terrestrial field lines well equatorward of the trapping boundary for electrons with energies greater than 45 keV are examined. Precipitation of low-energy electron intensities was characterized by isotropy for all pitch angles outside the atmospheric backscatter cone. The region of electron precipitation observed is associated with the diffuse aurora and with pulsating aurora in the postmidnight sector. Similar variations of the energetic electron intensities with energies above 45 keV were observed in the regions of fluctuating energy fluxes of low-energy electrons associated with auroral luminosity. The increases of energetic electron intensities were not coincident with those of the principal energy fluxes into the atmosphere, except when the average electron energy for the energy fluxes was unusually high, i.e., in the 10-keV range. Precipitation of electron intensities within these energy ranges is consistent with strong pitch angle diffusion of electron intensities near or at the magnetic equator by high-frequency wave turbulence, the effectiveness of which is modulated by perturbations attributable to micropulsations.

Frank, L. A.↗

Characteristics of the motions, turbulence intensity, diffusivity, flux of momentum and sensible heat in the upper atmosphere

Analyses of the meteorological rocket data obtained from an experiment conducted at 3-hour intervals at 8 western meridional rocket stations are presented. Large variations in the meridional wind contribute substantially to overall turbulence in the tropical stratosphere. The solar semidiurnal component of wind oscillations in the tropics was observed to be much higher than predicted by theory, often exceeding the magnitude of the diurnal amplitude throughout the stratosphere. The observed value of the solar diurnal amplitude in the stratosphere was in line with theoretical prediction. The solar terdiurnal amplitudes for temperature, meridional and zonal winds were non-negligible and must be considered in any harmonic analysis. Phase angle variation with height was rapid for all harmonics; however, there was general agreement between predicted and observed phase angles. Because of large changes in the mean winds in the mesosphere with season, harmonic determinations are difficult. There appear to be large zonal wind changes even within the same season as mentioned previously. Turbulence diffusivity in the upper stratosphere is greater near the equator than in the mid-latitudes.

Kao, S. K.↗

Determination of optimum sunlight concentration level in space for 3-5 cascade solar cells

Current-voltage curves were calculated for each cell in a cascade structure using a solar cell diode equation and superposition. Terms for the light generated current, diffusion current, space charge recombination current and series and shunt resistance are included. Individual current voltage curves are added in series with ohmic resistance losses for the cell interconnects to obtain the cascade cell performance. Temperature was varied with concentration, using several models, and ranged from 55 C at one Sun to between 80 and 200 C at 100 Suns. A variety of series resistance and internal resistances were used. Coefficients of the diffusion and recombination terms are strongly temperature dependent. The study indicates that maximum efficiency (30%) occurs in the 50 to 100X Sun concentration range, provided series resistance is below 0.015 ohm-sq cm and cell temperature is 80 C at 100 Suns.

Curtis, H. B.↗

Modelling of dimensional stability of fiber reinforced composite materials

Various methods of predicting the expansion and diffusion properties of composite laminates are reviewed. The prediction equations for continuous fiber composites can be applied to SMC composites as the effective fiber aspect ratio in the latter is large enough. The effect of hygrothermal expansion on the dimensional stability of composite laminates was demonstrated through the warping of unsymmetric graphite/epoxy laminates. The warping is very sensitive to the size of the panel, and to the moisture content which is in turn sensitive to the relative humidity in the environment. Thus, any long term creep test must be carried out in a humidity-controlled environment. Environmental effects in SMC composites and bulk polyester were studied under seven different environments. The SMC composites chosen are SMC-R25, SMC-R40, and SMC-R65.

Hahn, H. T.↗

The effect of convection upon charged particle transport in random magnetic fields

In a coordinate system moving with the plasma and random magnetic fields of a wind that blows with constant velocity in the direction of the guiding field, transport of energetic particles is described by a Boltzmann equation which is similar to the one that describes unconvected transport. Although this formulation is mathematically identical to that developed by Luhmann, which refers to the system where the guiding field is static, there are both practical and fundamental reasons to adopt the new approach. It leads to first-order approximate transport equations which are similar to those that apply in the absence of convection. However, these equations are more general than Parker's description of diffusion and convection, for they describe the coherent modes of transport that appear when the mean free path is large compared to the scale length for spatial variations of the guiding field, and they are valid for arbitrary wind velocity. The latter characteristic opens up new possibilities for analyzing particle transport in relativistic flows seen in some astronomical objects.

Earl, J. A.↗

The influence of gravity wave breaking on the general circulation of the middle atmosphere

The zonal mean solstice circulation of the global middle atmosphere is simulated using a semi-spectral numerical model. Radiative heating and cooling is computed by the algorithm of Wehrbein and Leovy. Mechanical dissipation is represented by the gravity wave breaking parameterization of Lindzen. An inertial adjustment parameterization is used to prevent the development of inertially unstable meridional shears near the equator. It is shown that gravity wave drag and diffusion in the mesosphere can account for the observed large departure from radiative equilibrium in both summer and winter. Experiments incorporating a forced stationary wavenumber 1 disturbance indicate that planetary wave EP flux convergences, although they may modify the mean flow profile significantly, cannot provide the major source of mechanical dissipation in the winter mesosphere. A simulated sudden warming is accompanied by an equally strong mesospheric cooling. This cooling is caused primarily by the relaxation of the polar mesosphere toward radiative equilibrium when the easterly mean winds in the polar stratosphere induced by the sudden warming reduce the transmission of gravity waves into the mesosphere.

Holton, J. R.↗

Numerical simulation and comparison with experiment for self-excited oscillations in a diffuser flow

This paper describes numerical simulations of self-excited oscillations in a two-dimensional transonic diffuser flow obtained by solving the Navier-Stokes equations with a two-equation turbulence model. Comparisons were made between the computational results and experimental data. For the mean flowfields, the agreement between computation and experiment is good for the wall pressures, shock location, and the separation and reattachment points. However, the thickness of the computed recirculation zone is about 50 percent of the measured thickness. For the fluctuating flowfields, a great deal of qualitative similarity exists between the computation and experiment; however, the predicted oscillation frequency is about 50 percent higher than the measured value. The formation of a succession of downstream-traveling counter-rotating vortices, as seen experimentally, is also vividly displayed in the numerical results.

Hsieh, T.↗

Effect of free stream turbulence on flow separation

An analysis is presented to study the influence of free stream turbulence on separation in the case of crossflow over a circular cylinder. An eddy diffusivity model has been formulated and the governing momentum equation has been solved numerically. A correlation parameter has been suggested to correlate the angle of separation. It has been found that the angle of separation increases with increasing free stream turbulence intensity.

Gorla, R. S. R.↗

Electron temperature and concentration in a thermal atomic oxygen source

A thermal atomic oxygen source for materials screening was built for NASA by Boeing Aerospace. The objective here was to use a microwave interferometer and Langmuir probe to characterize the electron concentration in this thermal atomic oxygen source. Typical operating conditions in the thermal atomic oxygen source were found to produce electron concentrations that were well below the detection threshold of the interferometer (10(exp 8) cm (sup -3)). The researchers calibrated (with the interferometer) the Langmuir probe at an artificially high plasma density and then used the circular and the square Langmuir probes to measure the low electron concentrations that exist during materials exposure tests. Electron concentration was measured as a function of power and position. The electrons were lost to the walls through ambipolar diffusion, and their concentration was accurately described by an equation. The electron concentration was proportional to power squared and decayed exponentially with distance.

Pedrow, Patrick Dennis↗

Validation of parameterization scheme for eddy diffusion from satellite data

The eddy diffusion coefficient K(yy) has been calculated usign LIMS for the months of December 1978 and January and February 1979. Two methods have been used. The first implements the suggestion made by Tung (1987) to parameterize the eddy transport as a diffusive process along isentropes. The second method integrates the equation relating the parcel displacements to the eddy velocity fields. The latter method uses a filtering on both space and time domains to isolate transients and is referred to as the 'spectral method'. Results from the first method are shown to be reliable only for quiescent periods, breaking down when the meridional gradient of potential vorticity is negligible. Results from the two methods are in agreement only for very disturbed conditions, when transience is readily isolated. It is concluded that the parameterizations suggested for eddy transport and calculated in this paper may be meaningful for quiet periods, but are not reliable for unsteady and very large amplitude disturbances.

Sassi, F.↗

Empirical Equation for Turbulent Forced-Convection Heat Transfer for Prandtl Numbers from 0.001 to 1000

A review is made of some of the experimental data and analyses applicable to convective heat transfer in fully turbulent flow in smooth tubes with liquid metals and viscous Newtonian fluids. An empirical equation is evolved that closely approximates heat-transfer values obtained from selected analyses and experimental data for Prandtl numbers from 0.001 to 1000. The terms included in the equation are Reynolds number, Prandtl number, and an empirical diffusivity ratio between heat and momentum.

vonGlahn, Uwe H.↗

A Self-Contained Mapping Closure Approximation for Scalar Mixing

Scalar turbulence exhibits interplays of coherent structures and random fluctuations over a broad range of spatial and temporal scales. This feature necessitates a probabilistic description of the scalar dynamics, which can be achieved comprehensively by using probability density functions (PDFs). Therefore, the challenge is to obtain the scalar PDFs (Lundgren 1967; Dopazo 1979). Generally, the evolution of a scalar is governed by three dynamical processes: advection, diffusion and reaction. In a PDF approach (Pope 1985), the advection and reaction can be treated exactly but the effect of molecular diffusion has to be modeled. It has been shown (Pope 1985) that the effect of molecular diffusion can be expressed as conditional dissipation rates or conditional diffusions. The currently used models for the conditional dissipation rates and conditional diffusions (Pope 1991) have resisted deduction from the fundamental equations and are unable to yield satisfactory results for the basic test cases of decaying scalars in isotropic turbulence, although they have achieved some success in a variety of individual cases. The recently developed mapping closure approach (Pope 1991; Chen, Chen & Kraichnan 1989; Kraichnan 1990; Klimenko & Pope 2003) provides a deductive method for conditional dissipation rates and conditional di usions, and the models obtained can successfully describe the shape relaxation of the scalar PDF from an initial double delta distribution to a Gaussian one. However, the mapping closure approach is not able to provide the rate at which the scalar evolves. The evolution rate has to be modeled. Therefore, the mapping closure approach is not closed. In this letter, we will address this problem.

He, Guo-Wei↗

Evaluation of a Multigrid Scheme for the Incompressible Navier-Stokes Equations

A fast multigrid solver for the steady, incompressible Navier-Stokes equations is presented. The multigrid solver is based upon a factorizable discrete scheme for the velocity-pressure form of the Navier-Stokes equations. This scheme correctly distinguishes between the advection-diffusion and elliptic parts of the operator, allowing efficient smoothers to be constructed. To evaluate the multigrid algorithm, solutions are computed for flow over a flat plate, parabola, and a Karman-Trefftz airfoil. Both nonlifting and lifting airfoil flows are considered, with a Reynolds number range of 200 to 800. Convergence and accuracy of the algorithm are discussed. Using Gauss-Seidel line relaxation in alternating directions, multigrid convergence behavior approaching that of O(N) methods is achieved. The computational efficiency of the numerical scheme is compared with that of Runge-Kutta and implicit upwind based multigrid methods.

Swanson, R. C.↗

Rotationally Symmetric Potential Flows

This paper includes the following topics: 1) Characteristic differential equations; 2) Treatment of practical examples; 3) First example: Diffuser; and 4) Second Example: Nozzle.

Schaefer, Manfred↗

Color Histogram Diffusion for Image Enhancement

Various color histogram equalization (CHE) methods have been proposed to extend grayscale histogram equalization (GHE) for color images. In this paper a new method called histogram diffusion that extends the GHE method to arbitrary dimensions is proposed. Ranges in a histogram are specified as overlapping bars of uniform heights and variable widths which are proportional to their frequencies. This diagram is called the vistogram. As an alternative approach to GHE, the squared error of the vistogram from the uniform distribution is minimized. Each bar in the vistogram is approximated by a Gaussian function. Gaussian particles in the vistoram diffuse as a nonlinear autonomous system of ordinary differential equations. CHE results of color images showed that the approach is effective.

Kim, Taemin↗

Quasi 1D Modeling of Mixed Compression Supersonic Inlets

The AeroServoElasticity task under the NASA Supersonics Project is developing dynamic models of the propulsion system and the vehicle in order to conduct research for integrated vehicle dynamic performance. As part of this effort, a nonlinear quasi 1-dimensional model of the 2-dimensional bifurcated mixed compression supersonic inlet is being developed. The model utilizes computational fluid dynamics for both the supersonic and subsonic diffusers. The oblique shocks are modeled utilizing compressible flow equations. This model also implements variable geometry required to control the normal shock position. The model is flexible and can also be utilized to simulate other mixed compression supersonic inlet designs. The model was validated both in time and in the frequency domain against the legacy LArge Perturbation INlet code, which has been previously verified using test data. This legacy code written in FORTRAN is quite extensive and complex in terms of the amount of software and number of subroutines. Further, the legacy code is not suitable for closed loop feedback controls design, and the simulation environment is not amenable to systems integration. Therefore, a solution is to develop an innovative, more simplified, mixed compression inlet model with the same steady state and dynamic performance as the legacy code that also can be used for controls design. The new nonlinear dynamic model is implemented in MATLAB Simulink. This environment allows easier development of linear models for controls design for shock positioning. The new model is also well suited for integration with a propulsion system model to study inlet/propulsion system performance, and integration with an aero-servo-elastic system model to study integrated vehicle ride quality, vehicle stability, and efficiency.

Kopasakis, George↗