Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffuse layer model”

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 73 records · Page 4

Effect of Gravity Level on the Particle Shape and Size During Zeolite Crystal Growth

A microscopic diffusion model is developed to represent solute transport in the boundary layer of a growing zeolite crystal. This model is used to describe the effect of gravity on particle shape and solute distribution. Particle dynamics and crystal growth kinetics serve as the boundary conditions of flow and convection-diffusion equations. A statistical rate theory is used to obtain the rate of solute transport across the growing interface, which is expressed in terms of concentration and velocity of solute species. Microgravity can significantly decrease the solute velocity across the growing interface compared to its earth-based counterpart. The extent of this reduction highly depends on solute diffusion constant in solution. Under gravity, the flow towards the crystal enhances solute transport rate across the growing interface while the flow away from crystals reduces this rate, suggesting a non-uniform growth rate and thus an elliptic final shape. However, microgravity can significantly reduce the influence of flow and obtain a final product with perfect spherical shape. The model predictions compare favorably with the data of space experiment of zeolites grown in space.

Song, Hong-Wei↗

Numerical Study of the Transient Vaporization of an Oxygen Droplet at Sub- and Super-Critical Conditions

Unsteady vaporization of a droplet in a high-pressure quiescent environment has been studied. With spherical symmetry and constant pressure, the process is diffusion controlled. At low pressures droplet heating is significant for most of the droplet lifetime, and an unsteady analysis of the gas phase is required. Then, the vaporization of a liquid oxygen droplet in gaseous hydrogen at moderate and high pressures is considered. At super-critical pressures, the surface temperature reaches the computed critical mixture value where a model for super-critical combustion is needed. The details of the diffusion layer of dissolved hydrogen do not significantly affect the results. Various methods are explored to compute the gas-phase density. Under the model's assumptions, homogeneous nucleation is not likely to occur.

Delplanque, J.-P.↗

Modeling the Uptake of Hydrogen Chloride onto Interior Spacecraft Materials

Hydrogen chloride (HCl) is a major combustion product from the pyrolysis of polyvinyl chloride (PVC) insulated electrical wires, a common spacecraft fire safety concern. Models at two different scales were developed to predict HCl uptake on anodized, chromate conversion coated (Iridite), and bare aluminum surfaces, as well as on Nomex fabric: a macroscopic one-step global surface reaction model where all the active sites are on the exterior surface, and a pore model where the interior active sites deeper into the oxide layer can also be accessed by HCl. Experiments were performed to calibrate kinetic and diffusion constants in the models. A cast acrylic test cell was used to measure the differences between the inlet and outlet concentration of HCl after inserting a sample rod of the test material. For the materials with a thin (< 200 µm) or no oxide layer, the macroscopic surface reaction model adequately predicts the experimental measurements. For the anodized aluminum with a thicker oxide layer, the pore model provided a better match to experimental results. The results will be discussed with respect to the spacecraft fire safety project (Saffire).

modeling↗

Properties of the +70 kilometers per second cloud toward HD 203664

I present high-resolution International Ultraviolet Explorer (IUE) spectra of the ultraviolet absorption in an intermediate-velocity interstellar cloud (nu(sub LSR) approximately equal to +70 km/s) toward HD 203664. The combined, multiple IUE images result in spectra with S/N = 15-40 and resolutions of approximately 20-25 km/s. The intermediate-velocity cloud absorption is present in ultraviolet lines of C II, C II(sup *), C IV, N I, O I, Mg I, Mg II, Al II, Al III, Si II, Si III, Si IV, S II, Cr II, Mn II, Fe II, and Zn II. The relative abundances of low-ionization species suggest an electron density of 0.15-0.34/cu cm and a temperature of 5300-6100 K in the neutral and weakly ionized gas. Given the presence of high-ionization gas tracers such as Si IV and C IV, ionized portions of the cloud probably contribute to the relatively large values of n(sub e) derived from measurements of the lower ionization species. The high-ionization species in the cloud have an abundance ratio, N(C IV)/N(Si IV) approximately equal to 4.5, similar to that inferred for collisionally ionized cloud interfaces at temperatures near 10(exp 5) K along other sight lines. When referenced to sulfur, the abundances of most elements in the cloud are within a factor of 5 of their solar values, which suggests that the +70 km/s gas has a previous origin in the Galactic disk despite a recent determination by Little et al. that the cloud lies at a distance of 200-1500 pc below the Galactic plane. I have checked this result against a model of the ionization for the diffuse ionized gas layer of the Galaxy and find that this conclusion is essentially unchanged as long as the ionization parameter is low as implied by the abundances of adjoining ionization states of aluminum and silicon. The processes responsible for the production of highly ionized gas in the +70 km/s cloud appear to be able to account for the inferred dust grain destruction as well.

Sembach, Kenneth R.↗

Multiple Scattering of Laser Pulses in Snow Over Ice: Modeling the Potential Bias in ICESat Altimetry

The primary goal of NASA's current ICESat and future ICESat2 missions is to map the altitude of the Earth's land ice with high accuracy using laser altimetry technology, and to measure sea ice freeboard. Ice however is a highly transparent optical medium with variable scattering and absorption properties. Moreover, it is often covered by a layer of snow with varying depth and optical properties largely dependent on its age. We describe a modeling framework for estimating the potential altimetry bias caused by multiple scattering in the layered medium. We use both a Monte Carlo technique and an analytical diffusion model valid for optically thick media. Our preliminary numerical results are consistent with estimates of the multiple scattering delay from laboratory measurements using snow harvested in Greenland, namely, a few cm. Planned refinements of the models are described.

Davis, A. B.↗

ACC Subduction by Mesoscales

The mesoscale contribution to Subduction in the Southern Ocean was recently studied by Salle'e and Rintoul (2011, SR) using the following mesoscale model. The adiabatic A-regime was modeled with the GM stream function, the diabatic D-regime was modeled with tapering functions,the D-A interface was taken to be at the mixed layer depth and the mesoscale diffusivity was either a constant or given by a 2D model. Since the resulting subductions were an order of magnitude smaller than the data of ±200 m/yr (Mazloff et al., 2010), SR showed that if instead of the above model-dependent mesoscale diffusivities, they employed the ones by Salle'e et al. (2008) from surface drifter observations, the subductions compared significantly better with the data. On those grounds, SR suggested a tenfold increase of the diffusivity. In this work, we suggest that since the mesoscale diffusivity is but one component of a muc large mesoscale parameterization, one should first assess the latter's overall performance followed by the assessment of the predicted ACC subduction. We employ the mesoscale model formulated31 in Canuto et al. (2018; 2019, that includes recent theoretical and observational advances and that was assessed against a variety of data including the output of 17 other OGCMs. The ACC diffusivities compare well with drifter data by Salle'e et al. (2008) and the ACC subduction rates are in agreement with the data.

V M Canuto↗

Three dimensional flow field inside compressor rotor, including blade boundary layers

The Reynolds stress equation, pressure strain correlation, and dissipative terms and diffusion are discussed in relation to turbulence modelling using the Reynolds stress model. Algebraic modeling of Reynolds stresses and calculation of the boundary layer over an axial cylinder are examined with regards to the kinetic energy model for turbulence modelling. The numerical analysis of blade and hub wall boundary layers, and an experimental study of rotor blade boundary layer in an axial flow compressor rotor are discussed. The Patankar-Spalding numerical method for two dimensional boundary layers is included.

Galmes, J. M.↗

Analytic modeling of a spray diffusion flame

A detailed model for a spray diffusion flame is described. The model is based on the boundary layer form of the equations of motion, with droplet transport accounted for using a discretized droplet size distribution function. Interphase transport of mass and energy are accounted for, with a flame-sheet model used to describe the combustion process on a droplet scale. Near dynamic equilibrium is assumed for the description of droplet transport; droplets can diffuse relative to the gas phase. Gas-phase mixing is accounted for using a two-equation turbulence model; buoyancy effects are included, with a temperature fluctuation equation used to account for buoyancy effects on turbulence structure. Thermal radiation from gas-phase CO2 and H2O is included. Gas-phase chemical kinetics are modeled using a 20-reaction, 10-species version of the advanced quasi-global chemical kinetics formulation. Results are compared with data for a vaporizing Freon spray and a pentane spray flame. It is shown that the computational approach provides a reasonably valid picture of the overall development of a spray diffusion flame, and, furthermore, provides a useful tool for the parametric examination of the spray combustion process.

Harsha, P. T.↗

Evaluation of Full Reynolds Stress Turbulence Models in FUN3D

Full seven-equation Reynolds stress turbulence models are a relatively new and promising tool for todays aerospace technology challenges. This paper uses two stress-omega full Reynolds stress models to evaluate challenging flows including shock-wave boundary layer interactions, separation and mixing layers. The Wilcox and the SSGLRR full second-moment Reynolds stress models are evaluated for four problems: a transonic two-dimensional diffuser, a supersonic axisymmetric compression corner, a compressible planar shear layer, and a subsonic axisymmetric jet. Simulation results are compared with experimental data and results using the more commonly used Spalart-Allmaras (SA) one-equation and the Menter Shear Stress Transport (SST) two-equation models.

computational fluid dynamics↗

Wind enhanced planetary escape: Collisional modifications

The problem of thermal escape is considered in which both the effects of thermospheric winds at the exobase and collisions below the exobase are included in a Monte Carlo calculation. The collisions are included by means of a collisional relaxation layer of a background gas which models the transition region between the exosphere and the thermosphere. The wind effects are considered in the limiting cases of vertical and horizontal flows. Two species are considered: terrestrial hydrogen and terrestrial helium. In the cases of terrestrial hydrogen the escape fluxes were found to be strongly filtered or throttled by collisions at high exospheric temperatures. The model is applied to molecular hydrogen diffusing through a methane relaxation layer under conditions possible on Titan. The results are similar to the case of terrestrial hydrogen with wind enhanced escape being strongly suppressed by collisions. It is concluded that wind enhanced escape is not an important process on Titan.

Curtis, S. A.↗

Diurnal forcing of planetary atmospheres

A free convection parameterization has been introduced into the Mars Planetary Boundary Layer Model (MPBL). Previously, the model would fail to generate turbulence under conditions of zero wind shear, even when statically unstable. This in turn resulted in erroneous results at the equator, for example, when the lack of Coriolis forcing allowed zero wind conditions. The underlying cause of these failures was the level 2 second-order turbulence closure scheme which derived diffusivities as algebraic functions of the Richardson number (the ratio of static stability to wind shear). In the previous formulation, the diffusivities were scaled by the wind shear--a convenient parameter since it is non-negative. This was the drawback that all diffusivities are zero under conditions of zero shear (viz., the free convection case). The new scheme tests for the condition of zero shear in conjunction with static instability and recalculates the diffusivities using a static stability scaling. The results for a simulation of the equatorial boundary layer at autumnal equinox are presented. (Note that after some wind shear is generated, the model reverts to the traditional diffusivity calculation.)

Houben, Howard C.↗

A three-dimensional axisymmetric calculation procedure for turbulent flows in a radial vaneless diffuser

An analytical model is proposed to calculate the three-dimensional axisymmetric turbulent flowfield in a radial vneless diffuser. The model assumes that the radial and tangential boundary layer profiles are approximated by power law profiles. Then, using the integrated radial and tangential momentum and continuity equations for the boundary layer and corresponding inviscid equations for the core flow, there results six ordinary differential equation in six unknowns which are easily solved using a Runge-Kutta technique. A model is also proposed for fully developed flow. The results using this technique were compared with the results from a three-dimensional viscous, axisymmetric duct code and with experimental data and good quantitative agreement was obtained.

Schumann, L. F.↗

Verification of Unstructured Grid Adaptation Components

Adaptive unstructured grid techniques have made limited impact on production analysis workflows where the control of discretization error is critical to obtaining reliable simulation results. Recent progress has matured a number of independent implementations of flow solvers, error estimation methods, and anisotropic grid adaptation mechanics. Known differences and previously unknown differences in grid adaptation components and their integrated processes are identified here for study. Unstructured grid adaptation tools are verified using analytic functions and the Code Comparison Principle. Three analytic functions with different smoothness properties are adapted to show the impact of smoothness on implementation differences. A scalar advection-diffusion problem with an analytic solution that models a boundary layer is adapted to test individual grid adaptation components. Laminar flow over a delta wing and turbulent flow over an ONERA M6 wing are verified with multiple, independent grid adaptation procedures to show consistent convergence to fine-grid forces and a moment. The scalar problems illustrate known differences in a grid adaptation component implementation and a previously unknown interaction between components. The wing adaptation cases in the current study document a clear improvement to existing grid adaptation procedures. The stage is set for the infusion of verified grid adaptation into production fluid flow simulations.

Park, Michael A.↗

Verification of Unstructured Grid Adaptation Components

Adaptive unstructured grid techniques have made limited impact on production analysis workflows where the control of discretization error is critical to obtaining reliable simulation results. Recent progress has matured a number of independent implementations of flow solvers, error estimation methods, and anisotropic grid adaptation mechanics. Known differences and previously unknown differences in grid adaptation components and their integrated processes are identified here for study. Unstructured grid adaptation tools are verified using analytic functions and the Code Comparison Principle. Three analytic functions with different smoothness properties are adapted to show the impact of smoothness on implementation differences. A scalar advection-diffusion problem with an analytic solution that models a boundary layer is adapted to test individual grid adaptation components. The scalar problems illustrate known differences in a grid adaptation component implementation and a previously unknown interaction between components. Laminar flow over a delta wing is verified with multiple, independent grid adaptation procedures to show consistent convergence to fine-grid forces and pitching moment.

Park, Michael A.↗

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.↗

Effects of turbulence model selection on the prediction of complex aerodynamic flows

Numerical simulations of viscous transonic flow over a circular-arc airfoil and in a diffuser are described. The simulations are made with a new computer program designed to serve as a tool in the development of improved turbulence models for complex flows. The program incorporates zero-, one-, and two-equation eddy viscosity models and includes a variety of subsonic and supersonic boundary conditions. The airfoil flow contains a shock-separated boundary-layer interaction that has resisted previous attempts at simulation. The diffuser flow also contains a shock-boundary-layer interaction, which has not been simulated previously. Calculations using standard turbulence models, developed originally for incompressible unseparated flows, are described. Results indicate that although there are interesting differences in predictions between the various models, none of them predict the flows accurately. Suggestions for improved turbulence models are discussed.

Coakley, T. J.↗

The stability of ground ice in the equatorial region of Mars

The lifetime of an unreplenished layer of ground ice lying within 30 deg of the Martian equator was examined within the context of the existing data base on Martian regolith and climate. Data on the partial pressure of H2O in the Martian atmosphere and the range of mean annual temperatures indicated the ground ice would be restricted to latitudes poleward of 40 deg. However, the ground ice near the poles may be a relic from early Martian geologic times held in place by a thin layer of regolith. Consideration of twelve model pore size distributions, similar to silt- and clay-type earth soils, was combined with a parallel pore model of gaseous diffusion to calculate the flux of H2O molecules escaping from the subsurface ground ice layer. Martian equatorial ground ice was found to be influenced by the soil structure, the magnitude of the geothermal gradient, the climatic desorption of CO2 from the regolith. It is concluded that equatorial ground ice is present on Mars only if a process of replenishment is active.

Clifford, S. M.↗

Downstream boundary effects on the frequency of self-excited oscillations in transonic diffuser flows

An investigation of downstream boundary effects on the frequency of self-excited oscillations in two-dimensional, separated transonic diffuser flows has been conducted numerically by solving the compressible, Reynolds-averaged, thin-layer Navier-Stokes equation with a two-equation turbulence model. It was found that the unsteady diffuser flowfields are very sensitive to the location of the downstream boundary. Extension of the diffuser downstream boundary significantly reduces the frequency and amplitude of oscillations for pressure, velocity and shock. Computational results suggest that the mechanism causing the self-excited oscillation changes from viscous convective wave dominated oscillations to inviscid acoustic wave dominated oscillations when the location of downstream boundary varies from 8.66 to 134.7 throat height. The existence of a suction slot in the experimental setup obscures the physical downstream boundary and, therefore, presents a difficulty for quantitative comparisons between computation and experiment.

Hsieh, T.↗