Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion coefficients”

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

MODIS Solar Diffuser Degradation Determination and Its Spectral Dependency

This study presents a modeling approach to improve solar diffuser (SD) degradation determination from SD stability monitor (SDSM) measurements. The MODIS instrument uses a SDto calibrate its reflective solar bands (RSBs) on-orbit. Due to the imperfectly designed SDSM sun view screen, the SD reflectance tracked by SDSM has large noise. The SDSM measurements noise is spectrally coherent and can be minimized by normalizing measurements to the least degradeddetector 9 (936 nm). In this study, a SD degradation model is used to determine the SDdegradation's wavelength dependency and the detector 9 degradation is estimated by the model solution.The results show the SD degradations measured at 6 SDSM detectors (554 _ 936 nm) have stable relationships, where the degradation is inversely proportion to 1/wavelength^4. The model estimated SD degradation at SDSM detector 9 wavelength (936 nm) is ~0.9% from 2002 to 2018.Based on the SD degradation model solution, the SD degradation at short/mid wave bands are estimated to improve short/mid wave bands calibration. The model can also be used to improve interpolating SD degradation at SDSM detectors to RSB wavelengths. Compared to linear interpolation, bands 9 and 10 show the largest differences of up to 0.3 and 0.4% respectively. These differences directly impact the calibration coefficients of these bands.

radiometric calibration↗

Parameterization of linear wave chemical transport in planetary atmospheres by eddy diffusion

Two-dimensional transport by linear waves is parameterized by wave velocity and temperature correlations and nonconservative properties of the atmosphere in a residual Eulerian framework. Explicit expressions derived for the eddy tracer flux vector are written in eddy diffusion tensor form, but it is demonstrated that strong coupling of chemistry and temperature through the rate constants, as is the case for ozone, requires additional terms proportional to the mixing ratio and not to its gradient. A variety of limiting cases are discussed. In the residual Eulerian framework the K(yz), K(zy) components are nonzero for a nonconservative atmosphere and/or a reactive tracer. Tracers have common transport coefficients in the limit where their chemistry is negligible in comparison to transience. It is shown that the usual one-dimensional vertical eddy diffusion continuity equation is rigorously valid only when the temperature dependence of the rate constants is negligible and the tracer has strong vertical stratification relative to its horizontal variability.

Strobel, D. F.↗

Effects of positive ion implantation into antireflection coating of silicon solar cells

The state of technological development of Si solar cells for highest obtained efficiency and radiation resistance is summarized. The various theoretical analyses of Si solar cells are reviewed. It is shown that factors controlling blue response are carrier diffusion length, surface recombination, impurity concentration profile in surface region, high level of surface impurity concentration (degeneracy), reflection coefficient of oxide, and absorption coefficient of Si. The theory of ion implantation of charge into the oxide antireflection coating is developed and side effects are discussed. The experimental investigations were directed at determining whether the blue response of Si solar cells could be improved by phosphorus ion charges introduced into the oxide antireflection coating.

Middleton, A. E.↗

Damage coefficients in low resistivity silicon

Electron and proton damage coefficients are determined for low resistivity silicon based on minority-carrier lifetime measurements on bulk material and diffusion length measurements on solar cells. Irradiations were performed on bulk samples and cells fabricated from four types of boron-doped 0.1 ohm-cm silicon ingots, including the four possible combinations of high and low oxygen content and high and low dislocation density. Measurements were also made on higher resistivity boron-doped bulk samples and solar cells. Major observations and conclusions from the investigation are discussed.

Srour, J. R.↗

Radiation damage in silicon NIP solar cells

The performance parameters of n(+) p p(+) silicon solar cells of varying thicknesses with boron doped p base resistivities of 1250 and 84 ohm centimeters were determined. High injection theory was used to analyze the experimental data. Results from an analysis of open circuit voltages show a much greater contribution to V(oc) from the back junction than is the case for the lower resistivity in common use. The base minority carrier distribution is seen to be significant in determining the contribution of V(B), the base contribution to V(oc). Although V(B) is small, its value increases with increasing radiation fluence. In this connection it is noted that, with illumination from the p(+) side, the sign of V(B) becomes positive, and V(B) itself becomes an additive term to V(oc). Diffusion lengths determined under high injection conditions are significantly greater than those obtained under low injection, while damage coefficients under low injection are higher than those obtained under high injection conditions.

Weinberg, I.↗

Finite element modeling of frictionally restrained composite interfaces

The use of special interface finite elements to model frictional restraint in composite interfaces is described. These elements simulate Coulomb friction at the interface, and are incorporated into a standard finite element analysis of a two-dimensional isolated fiber pullout test. Various interfacial characteristics, such as the distribution of stresses at the interface, the extent of slip and delamination, load diffusion from fiber to matrix, and the amount of fiber extraction or depression are studied for different friction coefficients. The results are compared to those obtained analytically using a singular integral equation approach, and those obtained by assuming a constant interface shear strength. The usefulness of these elements in micromechanical modeling of fiber-reinforced composite materials is highlighted.

Ballarini, Roberto↗

On Supersonic-Inlet Boundary-Layer Bleed Flow

Boundary-layer bleed in supersonic inlets is typically used to avoid separation from adverse shock-wave/boundary-layer interactions and subsequent total pressure losses in the subsonic diffuser and to improve normal shock stability. Methodologies used to determine bleed requirements are reviewed. Empirical sonic flow coefficients are currently used to determine the bleed hole pattern. These coefficients depend on local Mach number, pressure ratio, hole geometry, etc. A new analytical bleed method is presented to compute sonic flow coefficients for holes and narrow slots and predictions are compared with published data to illustrate the accuracy of the model. The model can be used by inlet designers and as a bleed boundary condition for computational fluid dynamic studies.

Harloff, Gary J.↗

Technical Report Series on Global Modeling and Data Assimilation: A Coupled Ocean-Atmosphere Radiative Model for Global Ocean Biogeochemical Models - Volume 22

An ocean-atmosphere radiative model (OARM) evaluates irradiance availability and quality in the water column to support phytoplankton growth and drive ocean thermodynamics. An atmospheric component incorporates spectral and directional effects of clear and cloudy skies as a function of atmospheric optical constituents, and spectral reflectance across the air-sea interface. An oceanic component evaluates the propagation of spectral and directional irradiance through the water column as a function of water, five phytoplankton groups, and chromophoric dissolved organic matter. It tracks the direct and diffuse streams from the atmospheric component, and a third stream, upwelling diffuse irradiance. The atmospheric component of OARM was compared to data sources at the ocean surface with a coefficient of determination (r2) of 0.97 and a root mean square of 12.1%.

Gregg, Watson W.↗

Parameters Impacting Columnated Granular Soil Pneumatic Seal Performance

The ability of a column of loose granular soil to form a pneumatic seal was investigated by varying the diameter of the soil column, the effective column height, and the level of compaction in the soil. The soil column diameter was tested at three levels using pipes with inner diameters measuring 5.08, 10.16, and 15.24 centimeters (2, 4, and 6 inches). Soil was filled in each pipe to form 15, 30, and 45 centimeter (6, 12, and 18 inch) tall soil columns. Each diameter/height configuration was also tested at three levels of soil compaction, compared by calculating the bulk density of the soil with mass and volume measurements. GRC-1a simulant was used, with approximate low/medium/high bulk densities of 1.6, 1.75, and 1.9 g/cc achieved with a combination of vibration and tamping. The order of tests for a given column diameter was randomized and repeated three times. With the top of the soil open to atmosphere at room temperature, compressed air was injected through a small diffuser at the column base with several small downward-facing holes. The number and size of these holes was scaled such that a constant total inlet orifice area to column cross sectional area ratio was maintained for each column diameter. Inlet air pressure was slowly increased via a precision regulator to preserve quasi-static equilibrium in the soil column to minimize the impact of dynamics. Air pressure was increased all the way through the static and bubbling regimes until slugging or turbulent behavior was observed in the soil to ensure that the entire static regime had been captured during data collection. A three-factor, three-level analysis of variance (ANOVA) statistical analysis was performed on the resulting data to determine the extent to which each physical parameter impacted the soil column seal performance. It was concluded that there is statistically significant evidence that column height, and the interaction between column height and diameter impact soil seal performance. In all other cases there was insufficient data to identify a statistically significant causal relationship. Additionally, plots were generated comparing experimental data to the predictive formula developed by Ogino et al. for fluidized beds in 1993. Because this model was developed for industrial spouted fluidized beds, the accuracy of its output prior to fluidization in the static seal ‘edge case’ is unknown, especially considering in this application the working gas was diffused across the column base rather than being injected through a spout. Further, the fidelity of the model had not yet been tested with lunar soil simulants. Plotting the Ogino et al. model alongside test data allows for a more intuitive sense of the impact of test parameters on soil seal performance, as well as providing a quick means to further tune this predictive model for more accurate use with static seals across granular lunar soil simulants. The model provided by Ogino et al. was further tuned using test data to determine the degree to which a spouted bed model could be applied to a slightly modified set of testing conditions: lunar soil simulants and a more diffuse, homogeneous application of pneumatic pressure. A Matlab script was created to test different values for the leading coefficient and exponents in Ogino’s formula. The script swept preset ranges, then iterated with higher resolutions across narrower ranges to converge on the optimal value for each parameter. The resulting adjusted model was compared to the original, as well as test data with noticeable improvements across the entire test domain.

Jack Stewart↗

Analytical investigation of ram-jet-engine performance in flight Mach number range from 3 to 7

An analytical investigation was made of the performance of isolated ram-jet engines in the flight Mach number range from 3 to 7 for two types of diffuser, a high-efficiency diffuser, and a normal-shock diffuser. The fuel was assumed to be a hydrocarbon similar to gasoline. The conclusions reached are: (1) a design altitude of about 100,000 feet is desirable for a high-efficiency high Mach number ram jet on the basis of engine construction and performance; and (2) although greater thrust could be obtained with other fuels, gasoline provides sufficient energy release for maximum engine efficiency in the flight Mach number range investigated. The maximum engine efficiency calculated was 0.47, which occurred at a Mach number of 5. At a Mach number of 7, the maximum propulsive-thrust coefficient was 0.57.

ENGINES, RAM JET↗

An unconditionally-stable central differencing scheme for high Reynolds number flows

The central difference scheme (CDS) is a second order accurate scheme which is free of numerical diffusion (in the second order sense) and is simple to implement: However, for grid Peclet numbers larger than 2, the CDS leads to over- and undershoots and is unstable. The present paper describes a method, called CONDIF, which retains the essential nature of the CDS but eliminates the over- and under-shoots. It leads to unconditionally positive coefficients and, in the limit, approaches the CDS for all values of grid Peclet numbers. The CONDIF modifies the CDS by introducing a controlled amount of numerical diffusion based on the local gradients. In the worst case the scheme yields results similar to those of the hybrid scheme. This paper reports the results obtained from CONDIF for a number of test problems which have been widely used for comparative study of numerical schemes in the published literature. For most of these problems, the CONDIF results are significantly more accurate than the hybrid scheme at high Peclet numbers. In particular, the CONDIF scheme depicts much lower level of numerical diffusion than the hybrid scheme even when the Peclet number is very high and the flow is at large angles to the grid.

Runchal, Akshai K.↗

Actinic Flux Calculations: A Model Sensitivity Study

calculate direct and diffuse surface irradiance and actinic flux (downwelling (2p) and total (4p)) for the reference model. Sensitivity analysis has shown that the accuracy of the radiative transfer flux calculations for a unit ETS (i.e. atmospheric transmittance) together with a numerical interpolation technique for the constituents' vertical profiles is better than 1% for SZA less than 70(sub o) and wavelengths longer than 310 nm. The differences increase for shorter wavelengths and larger SZA, due to the differences in pseudo-spherical correction techniques and vertical discretetization among the codes. Our sensitivity study includes variation of ozone cross-sections, ETS spectra and the effects of wavelength shifts between vacuum and air scales. We also investigate the effects of aerosols on the spectral flux components in the UV and visible spectral regions. The "aerosol correction factors" (ACFs) were calculated at discrete wavelengths and different SZAs for each flux component (direct, diffuse, reflected) and prescribed IPMMI aerosol parameters. Finally, the sensitivity study was extended to calculation of selected photolysis rates coefficients.

Krotkov, Nickolay A.↗

Three Stage Cool Flame Droplet Burning Behavior of n-Alkane Droplets at Elevated Pressure Conditions: Hot, Warm and Cool Flame

Transient, isolated n-alkane droplet combustion is simulated at elevated pressure for helium-diluent substituted-air mixtures. We report the presence of unique quasi-steady, three-stage burning behavior of large sphero-symmetric n-alkane droplets at these elevated pressures and helium substituted ambient fractions. Upon initiation of reaction, hot-flame diffusive burning of large droplets is initiated that radiatively extinguishes to establish cool flame burning conditions in nitrogen/oxygen “air” at atmospheric and elevated pressures. However, at elevated pressure and moderate helium substitution for nitrogen ( X He > 20%), the initiated cool flame burning proceeds through two distinct, quasi-steady-state, cool flame burning conditions. The classical “Hot flame ” ( ∼1500 K) radiatively extinguishes into a “Warm flame ” burning mode at a moderate maximum reaction zone temperature ( ∼ 970 K), followed by a transition to a lower temperature ( ∼765 K), quasi-steady “Cool flame ” burning condition. The reaction zone (“flame”) temperatures are associated with distinctly different yields in intermediate reaction products within the reaction zones and surrounding near-field, and the flame-standoff ratios characterizing each burning mode progressively decrease. The presence of all three stages first appears with helium substitution near 20%, and the duration of each stage is observed to be strongly dependent on helium substitutions level between 20–60%. For helium substitution greater than 60%, the hot flame extinction is followed by only the lower temperature cool flame burning mode. In addition to the strong coupling between the diffusive loss of both energy and species and the slowly evolving degenerate branching in the low and negative temperature coefficient (NTC) kinetic regimes, the competition between the low-temperature chain branching and intermediate-temperature chain termination reactions control the “Warm” and “Cool ” flame quasi-steady conditions and transitioning dynamics.

Droplet Combustion↗

Diffusion in the lower stratosphere as determined from lidar measurements of volcanic aerosol dispersion

Lidar measurements of the stratospheric aerosol layer from the Fuego volcanic eruption in 1974 are analyzed to yield estimates of effective vertical mixing coefficients K(z). The data at 19 deg N latitude give K(z) = 6.6 x 10 to the 2nd sq cm/sec for the altitude range of 19 + or - 1.5 km, while the data at 37 deg N yield K(z) = 13 x 10 to the 2nd sq cm/sec for 19 + or - 3.5 km. The variability in these values is about a factor of 2. These derived K(z) values are interpreted as being due mainly to diffusion.

Remsberg, E. E.↗

Parametric studies with an atmospheric diffusion model that assesses toxic fuel hazards due to the ground clouds generated by rocket launches

Parametric studies were made with a multilayer atmospheric diffusion model to place quantitative limits on the uncertainty of predicting ground-level toxic rocket-fuel concentrations. Exhaust distributions in the ground cloud, cloud stabilized geometry, atmospheric coefficients, the effects of exhaust plume afterburning of carbon monoxide CO, assumed surface mixing-layer division in the model, and model sensitivity to different meteorological regimes were studied. Large-scale differences in ground-level predictions are quantitatively described. Cloud alongwind growth for several meteorological conditions is shown to be in error because of incorrect application of previous diffusion theory. In addition, rocket-plume calculations indicate that almost all of the rocket-motor carbon monoxide is afterburned to carbon dioxide CO2, thus reducing toxic hazards due to CO. The afterburning is also shown to have a significant effect on cloud stabilization height and on ground-level concentrations of exhaust products.

Stewart, R. B.↗

The prediction of nozzle performance and heat transfer in hydrogen/oxygen rocket engines with transpiration cooling, film cooling, and high area ratios

An advanced engineering computational model has been developed to aid in the analysis and design of hydrogen/oxygen chemical rocket engines. The complete multi-species, chemically reacting and diffusing Navier-Stokes equations are modelled, finite difference approach that is tailored to be conservative in an axisymmetric coordinate system for both the inviscid and viscous terms. Demonstration cases are presented for a 1030:1 area ratio nozzle, a 25 lbf film cooled nozzle, and transpiration cooled plug-and-spool rocket engine. The results indicate that the thrust coefficient predictions of the 1030:1 nozzle and the film cooled nozzle are within 0.2 to 0.5 percent, respectively, of experimental measurements when all of the chemical reaction and diffusion terms are considered. Further, the model's predictions agree very well with the heat transfer measurements made in all of the nozzle test cases. The Soret thermal diffusion term is demonstrated to have a significant effect on the predicted mass fraction of hydrogen along the wall of the nozzle in both the laminar flow 1030:1 nozzle and the turbulent plug-and-spool rocket engine analysis cases performed. Further, the Soret term was shown to represent a significant fraction of the diffusion fluxes occurring in the transpiration cooled rocket engine.

Kacynski, Kenneth J.↗

Uncertainty Determination for Aeroheating in Uranus and Saturn Probe Entries by the Monte Carlo Method

The 2013-2022 Decaedal survey for planetary exploration has identified probe missions to Uranus and Saturn as high priorities. This work endeavors to examine the uncertainty for determining aeroheating in such entry environments. Representative entry trajectories are constructed using the TRAJ software. Flowfields at selected points on the trajectories are then computed using the Data Parallel Line Relaxation (DPLR) Computational Fluid Dynamics Code. A Monte Carlo study is performed on the DPLR input parameters to determine the uncertainty in the predicted aeroheating, and correlation coefficients are examined to identify which input parameters show the most influence on the uncertainty. A review of the present best practices for input parameters (e.g. transport coefficient and vibrational relaxation time) is also conducted. It is found that the 2(sigma) - uncertainty for heating on Uranus entry is no more than 2.1%, assuming an equilibrium catalytic wall, with the uncertainty being determined primarily by diffusion and H(sub 2) recombination rate within the boundary layer. However, if the wall is assumed to be partially or non-catalytic, this uncertainty may increase to as large as 18%. The catalytic wall model can contribute over 3x change in heat flux and a 20% variation in film coefficient. Therefore, coupled material response/fluid dynamic models are recommended for this problem. It was also found that much of this variability is artificially suppressed when a constant Schmidt number approach is implemented. Because the boundary layer is reacting, it is necessary to employ self-consistent effective binary diffusion to obtain a correct thermal transport solution. For Saturn entries, the 2(sigma) - uncertainty for convective heating was less than 3.7%. The major uncertainty driver was dependent on shock temperature/velocity, changing from boundary layer thermal conductivity to diffusivity and then to shock layer ionization rate as velocity increases. While radiative heating for Uranus entry was negligible, the nominal solution for Saturn computed up to 20% radiative heating at the highest velocity examined. The radiative heating followed a non-normal distribution, with up to a 3x variation in magnitude. This uncertainty is driven by the H(sub 2) dissociation rate, as H(sub 2) that persists in the hot non-equilibrium zone contributes significantly to radiation.

Palmer, Grant↗

Laser-Induced Fluorescence Measurements and Modeling of Nitric Oxide in Counterflow Diffusion Flames

The feasibility of making quantitative nonintrusive NO concentration ([NO]) measurements in nonpremixed flames has been assessed by obtaining laser-induced fluorescence (LIF) measurements of [NO] in counterflow diffusion flames at atmospheric and higher pressures. Comparisons at atmospheric pressure between laser-saturated fluorescence (LSF) and linear LIF measurements in four diluted ethane-air counterflow diffusion flames with strain rates from 5 to 48/s yielded excellent agreement from fuel-lean to moderately fuel-rich conditions, thus indicating the utility of a model-based quenching correction technique, which was then extended to higher pressures. Quantitative LIF measurements of [NO] in three diluted methane-air counterflow diffusion flames with strain rates from 5 to 35/s were compared with OPPDIF model predictions using the GRI (version 2.11) chemical kinetic mechanism. The comparisons revealed that the GRI mechanism underpredicts prompt-NO by 30-50% at atmospheric pressure. Based on these measurements, a modified reaction rate coefficient for the prompt-NO initiation reaction was proposed which causes the predictions to match experimental data. Temperature measurements using thin filament pyrometry (TFP) in conjunction with a new calibration method utilizing a near-adiabatic H2-air Hencken burner gave very good comparisons with model predictions in these counterflow diffusion flames. Quantitative LIF measurements of [NO] were also obtained in four methane-air counterflow partially-premixed flames with fuel-side equivalence ratios (phi(sub B)) of 1.45, 1.6, 1.8 and 2.0. The measurements were in excellent agreement with model predictions when accounting for radiative heat loss. Spatial separation between regions dominated by the prompt and thermal NO mechanisms was observed in the phi(sub B) = 1.45 flame. The modified rate coefficient proposed earlier for the prompt-NO initiation reaction improved agreement between code predictions and measurements in the region where prompt-NO dominates. Finally, LIF measurements of NO were obtained in counterflow diffusion flames at 2 to 5 atm. Comparisons between [NO] measurements and predictions show that the GRI mechanism underpredicts prompt-NO by a factor of two to three at all pressures. In general, the results indicate a need for refinement of the CH chemistry, especially the pressure-dependent CH formation and destruction reactions.

Ravikrishna, Rayavarapu V.↗