Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Conditional diffusion models”

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 757 records · Page 42

Energy balance in the solar transition region. I - Hydrostatic thermal models with ambipolar diffusion

The energy balance in the lower transition region is analyzed by constructing theoretical models which satisfy the energy balance constraint. The energy balance is achieved by balancing the radiative losses and the energy flowing downward from the corona. This energy flow is mainly in two forms: conductive heat flow and hydrogen ionization energy flow due to ambipolar diffusion. Hydrostatic equilibrium is assumed, and, in a first calculation, local mechanical heating and Joule heating are ignored. In a second model, some mechanical heating compatible with chromospheric energy-balance calculations is introduced. The models are computed for a partial non-LTE approach in which radiation departs strongly from LTE but particles depart from Maxwellian distributions only to first order. The results, which apply to cases where the magnetic field is either absent, or uniform and vertical, are compared with the observed Lyman lines and continuum from the average quiet sun. The approximate agreement suggests that this type of model can roughly explain the observed intensities in a physically meaningful way, assuming only a few free parameters specified as chromospheric boundary conditions.

Fontenla, J. M.↗

Efficient Low Dissipative High Order Schemes for Multiscale MHD Flows, I: Basic Theory

The objective of this paper is to extend our recently developed highly parallelizable nonlinear stable high order schemes for complex multiscale hydrodynamic applications to the viscous MHD equations. These schemes employed multiresolution wavelets as adaptive numerical dissipation controls t o limit the amount of and to aid the selection and/or blending of the appropriate types of dissipation to be used. The new scheme is formulated for both the conservative and non-conservative form of the MHD equations in curvilinear grids. The four advantages of the present approach over existing MHD schemes reported in the open literature are as follows. First, the scheme is constructed for long-time integrations of shock/turbulence/combustion MHD flows. Available schemes are too diffusive for long-time integrations and/or turbulence/combustion problems. Second, unlike exist- ing schemes for the conservative MHD equations which suffer from ill-conditioned eigen- decompositions, the present scheme makes use of a well-conditioned eigen-decomposition obtained from a minor modification of the eigenvectors of the non-conservative MHD equations t o solve the conservative form of the MHD equations. Third, this approach of using the non-conservative eigensystem when solving the conservative equations also works well in the context of standard shock-capturing schemes for the MHD equations. Fourth, a new approach to minimize the numerical error of the divergence-free magnetic condition for high order schemes is introduced. Numerical experiments with typical MHD model problems revealed the applicability of the newly developed schemes for the MHD equations.

Sjoegreen, Bjoern↗

Combustion of liquid fuel droplets in supercritical conditions

A comprehensive analysis of liquid-fuel droplet combustion in both sub- and super-critical environments has been conducted. The formulation is based on the complete conservation equations for both gas and liquid phases, and accommodates finite-rate chemical kinetics and a full treatment of liquid-vapor phase equilibrium at the droplet surface. The governing equations and the associated interface boundary conditions are solved numerically using a fully coupled, implicit scheme with the dual time-stepping integration technique. The model is capable of treating the entire droplet history, including the transition from the subcritical to the supercritical state. As a specific example, the combustion of n-pentane fuel droplets in air is studied for pressures of 5-140 atm. Results indicate that the ambient gas pressure exerts significant control of droplet gasification and burning processes through its influences on the fluid transport, gas/liquid interface thermodynamics, and chemical reactions. The droplet gasification rate increases progressively with pressure. However, the data for the overall burnout time exhibits a significant variation near the critical burning pressure, mainly as a result of reduced mass-diffusion rate and latent heat of vaporization with increased pressure. The influence of droplet size on the burning characteristics is also noted.

Shuen, J. S.↗

Navier-Stokes analysis and experimental data comparison of compressible flow within ducts

Many aircraft employ ducts with centerline curvature or changing cross-sectional shape to join the engine with inlet and exhaust components. S-ducts convey air to the engine compressor from the intake and often decelerate the flow to achieve an acceptable Mach number at the engine compressor by increasing the cross-sectional area downstream. Circular-to-rectangular transition ducts are used on aircraft with rectangular exhaust nozzles to connect the engine and nozzle. To achieve maximum engine performance, the ducts should minimize flow total pressure loss and total pressure distortion at the duct exit. Changes in the curvature of the duct centerline or the duct cross-sectional shape give rise to streamline curvature which causes cross stream pressure gradients. Secondary flows can be caused by deflection of the transverse vorticity component of the boundary layer. This vortex tilting results in counter-rotating vortices. Additionally, the adverse streamwise pressure gradient caused by increasing cross-sectional area can lead to flow separation. Vortex pairs have been observed in the exit planes of both duct types. These vortices are due to secondary flows induced by pressure gradients resulting from streamline curvature. Regions of low total pressure are produced when the vortices convect boundary layer fluid into the main flow. The purpose of the present study is to predict the measured flow field in a diffusing S-duct and a circular-to-rectangular transition duct with a full Navier-Stokes computer program, PARC3D, and to compare the numerical predictions with new detailed experimental measurements. The work was undertaken to extend previous studies and to provide additional CFD validation data needed to help model flows with strong secondary flow and boundary layer separation. The S-duct computation extends the study of Smith et al, and Harloff et al, which concluded that the computation might be improved by using a finer grid and more advanced turbulence models. The present study compares results for both the Baldwin-Lomas and k-epsilon turbulence models and is conducted with a refined grid. For the transition duct, two inlet conditions were considered, the first with straight flow and the second with swirling flow. The first case permits examination of the effects of the geometric transition on the flow field, while the second case includes the rotational flow effect characteristic of a gas turbine engine.

Harloff, G. J.↗

Physical electronics and surface physics

Current analytical and experimental investigations of the influence of surface processes on the rate of permeation of gases through solids are summarized. A brief description is presented of an analytical model that provides an approximate prediction of the dependence of the permeation rate on the rates of surface processes relative to the rate of the diffusion process. Results of experiments indicate that the desorption process of hydrogen from various materials is influenced by surface impurities to the degree that substantial changes in the spatial distribution of the desorbed molecules are observed when the impurities are varied. Except possibly for H-Nb, the permeation measurements indicate that the desorption process is not influenced to the degree that detectable changes in the permeation rates are observed when the impurities and temperatures are varied. On the basis of analysis, it is unlikely that the desorption process will be the rate-controlling step under the conditions of the present experiments.

Stickney, R. E.↗

Mass loading and velocity diffusion models for heavy pickup ions at comet Grigg-Skjellerup

We compare model predictions of cometary water group ion densities and the solar wind slow down with measurements made by the Giotto Johnstone plasma analyzer implanted ion sensor at the encounter with comet Grigg-Skjellerup (G-S) on July 10, 1992. The observed slope of the ion density profile on approach to the comet is unexpectedly steep. Possible explanations for this are discussed. We present also a preliminary investigation of the quasilinear velocity-space diffusion of the implanted heavy ion population at G-S using a transport equation including souce, convection, adiabatic compression, and velocity diffusion terms. Resulting distributions are anisotropic, in agreement with observations. We consider theoretically the waves that may be generated by the diffusion process for the observed solar wind conditions. At initial ion injections, waves are generated at omega approximately Omega(sub i) the ion gyrofrequency, and lower frequencies are predicted for diffusion toward a bispherical shell.

Huddleston, D. E.↗

Numerical Modeling of Flow Control in a Boundary-Layer-Ingesting Offset Inlet Diffuser at Transonic Mach Numbers

This paper will investigate the validation of a NASA developed, Reynolds-averaged Navier-Stokes (RANS) flow solver, OVERFLOW, for a boundary-layer-ingesting (BLI) offset (S-shaped) inlet in transonic flow with passive and active flow control devices as well as the baseline case. Numerical simulations are compared to wind tunnel results of a BLI inlet conducted at the NASA Langley 0.3-Meter Transonic Cryogenic Tunnel. Comparisons of inlet flow distortion, pressure recovery, and inlet wall pressures are performed. The numerical simulations are compared to the BLI inlet data at a freestream Mach number of 0.85 and a Reynolds number of approximately 2 million based on the length of the fan-face diameter. The numerical simulations with and without wind tunnel walls are performed, quantifying effects of the tunnel walls on the BLI inlet flow measurements. The wind tunnel test evaluated several different combinations of jet locations and mass flow rates as well as a vortex generator (VG) vane case. The numerical simulations will be performed on a single jet configuration for varying actuator mass flow rates at a fix inlet mass flow condition. Validation of the numerical simulations for the VG vane case will also be performed for varying inlet mass flow rates. Overall, the numerical simulations were able to predict the baseline circumferential flow distortion, DPCPavg, very well for comparisons made within the designed operating range of the BLI inlet. However the CFD simulations did predict a total pressure recovery that was 0.01 lower than the experiment. Numerical simulations of the baseline inlet flow also showed good agreement with the experimental inlet centerline surface pressures. The vane case showed that the CFD predicted the correct trends in the circumferential distortion for varying inlet mass flow but had a distortion level that was nearly twice as large as the experiment. Comparison to circumferential distortion measurements for a 15 deg clocked 40 probe rake indicated that the circumferential distortion levels are very sensitive to the symmetry of the flow and that a miss alignment of the vanes in the experiment could have resulted in this difference. The numerical simulations of the BLI inlet with jets showed good agreement with the circumferential inlet distortion levels for a range of jet actuator mass flow ratios at a fixed inlet mass flow rate. The CFD simulations for the jet case also predicted an average total pressure recovery that was 0.01 lower than the experiment as was seen in the baseline. Comparison of the flow features the jet case revealed that the CFD predicted a much larger vortex at the engine fan-face when compare to the experiment.

Allan Brian G.↗

Integration of the Radiation Belt Environment Model Into the Space Weather Modeling Framework

We have integrated the Fok radiation belt environment (RBE) model into the space weather modeling framework (SWMF). RBE is coupled to the global magnetohydrodynamics component (represented by the Block-Adaptive-Tree Solar-wind Roe-type Upwind Scheme, BATS-R-US, code) and the Ionosphere Electrodynamics component of the SWMF, following initial results using the Weimer empirical model for the ionospheric potential. The radiation belt (RB) model solves the convection-diffusion equation of the plasma in the energy range of 10 keV to a few MeV. In stand-alone mode RBE uses Tsyganenko's empirical models for the magnetic field, and Weimer's empirical model for the ionospheric potential. In the SWMF the BATS-R-US model provides the time dependent magnetic field by efficiently tracing the closed magnetic field-lines and passing the geometrical and field strength information to RBE at a regular cadence. The ionosphere electrodynamics component uses a two-dimensional vertical potential solver to provide new potential maps to the RBE model at regular intervals. We discuss the coupling algorithm and show some preliminary results with the coupled code. We run our newly coupled model for periods of steady solar wind conditions and compare our results to the RB model using an empirical magnetic field and potential model. We also simulate the RB for an active time period and find that there are substantial differences in the RB model results when changing either the magnetic field or the electric field, including the creation of an outer belt enhancement via rapid inward transport on the time scale of tens of minutes.

Glocer, A.↗

Use of GRACE Terrestrial Water Storage Retrievals to Evaluate Model Estimates by the Australian Water Resources Assessment System

Terrestrial water storage (TWS) estimates retrievals from the Gravity Recovery and Climate Experiment (GRACE) satellite mission were compared to TWS modeled by the Australian Water Resources Assessment (AWRA) system. The aim was to test whether differences could be attributed and used to identify model deficiencies. Data for 2003 2010 were decomposed into the seasonal cycle, linear trends and the remaining de-trended anomalies before comparing. AWRA tended to have smaller seasonal amplitude than GRACE. GRACE showed a strong (greater than 15 millimeter per year) drying trend in northwest Australia that was associated with a preceding period of unusually wet conditions, whereas weaker drying trends in the southern Murray Basin and southwest Western Australia were associated with relatively dry conditions. AWRA estimated trends were less negative for these regions, while a more positive trend was estimated for areas affected by cyclone Charlotte in 2009. For 2003-2009, a decrease of 7-8 millimeter per year (50-60 cubic kilometers per year) was estimated from GRACE, enough to explain 6-7% of the contemporary rate of global sea level rise. This trend was not reproduced by the model. Agreement between model and data suggested that the GRACE retrieval error estimates are biased high. A scaling coefficient applied to GRACE TWS to reduce the effect of signal leakage appeared to degrade quantitative agreement for some regions. Model aspects identified for improvement included a need for better estimation of rainfall in northwest Australia, and more sophisticated treatment of diffuse groundwater discharge processes and surface-groundwater connectivity for some regions.

van Dijk, A. I. J. M.↗

Decoupled Method for Reconstruction of Surface Conditions From Internal Temperatures On Ablative Materials With Uncertain Recession Model

Obtaining measurements of flight environments on ablative heat shields is both critical for spacecraft development and extremely challenging due to the harsh heating environment and surface recession. Thermocouples installed several millimeters below the surface are commonly used to measure the heat shield temperature response, but an ill-posed inverse heat conduction problem must be solved to reconstruct the surface heating environment from these measurements. Ablation can contribute substantially to the measurement response making solutions to the inverse problem strongly dependent on the recession model, which is often poorly characterized. To enable efficient surface reconstruction for recession model sensitivity analysis, a method for decoupling the surface recession evaluation from the inverse heat conduction problem is presented. The decoupled method is shown to provide reconstructions of equivalent accuracy to the traditional coupled method but with substantially reduced computational effort. These methods are applied to reconstruct the environments on the Mars Science Laboratory heat shield using diffusion limit and kinetically limited recession models.

Oliver, A. Brandon↗

Basal melt beneath whillans ice stream and ice streams A and C

We have used a recently derived map of the velocity of Whillans Ice Stream and Ice Streams A and C to help estimate basal melt. Temperature was modeled with a simple vertical advection-diffusion equation, 'tuned' to match temperature profiles. We find that most of the melt occurs beneath the tributaries where larger basal shear stresses and thicker ice favors greater melt (e.g., 10-20 mm/yr). The occurrence of basal freezing is predicted beneath much of the ice plains of Ice Stream C and Whillans Ice Stream. Modelled melt rates for when Ice Stream C was active suggest there was just enough melt water generated in its tributaries to balance basal freezing on its ice plain. Net basal melt for Whillans Ice Stream is positive due to smaller basal temperature gradients. Modelled temperatures on Whillans Ice Stream, however, were constrained by a single temperature profile at UpB. Basal temperature gradients for Whillans B1 and Ice Stream A may have conditions more similar to those beneath Ice Streams C and D, in which case, there may not be sufficient melt to sustain motion. This would be consistent with the steady deceleration of Whillans stream over the last few decades.

glaciers ice sheets↗

Computational fluid dynamics analysis in support of the simplex turbopump design

Simplex is a turbopump that is being developed at NASA/Marshall Space Flight Center (MSFC) by an in-house team. The turbopump consists of a single-stage centrifugal impeller, vaned-diffuser pump powered by a single-stage, axial, supersonic, partial admission turbine. The turbine is driven by warm gaseous oxygen tapped off of the hybrid motor to which it will be coupled. Rolling element bearings are cooled by the pumping fluid. Details of the configuration and operating conditions are given by Marsh. CFD has been used extensively to verify one-dimensional (1D) predictions, assess aerodynamic and hydrodynamic designs, and to provide flow environments. The complete primary flow path of the pump-end and the hot gas path of the turbine, excluding the inlet torus, have been analyzed. All CFD analyses conducted for the Simplex turbopump employed the pressure based Finite Difference Navier-Stokes (FDNS) code using a standard kappa-epsilon turbulence model with wall functions. More detailed results are presented by Garcia et. al. To support the team, loading and temperature results for the turbine rotor were provided as inputs to structural and thermal analyses, and blade loadings from the inducer were provided for structural analyses.

Garcia, Roberto↗

Investigation of Dynamic Oxygen Adsorption in Molten Solder Jetting Technology

Surface tension forces play a critical role in fluid dynamic phenomena that are important in materials processing. The surface tension of liquid metals has been shown to be very susceptible to small amounts of adsorbed oxygen. Consequently, the kinetics of oxygen adsorption can influence the capillary breakup of liquid-metal jets targeted for use in electronics assembly applications, where low-melting-point metals (such as tin-containing solders) are utilized as an attachment material for mounting of electronic components to substrates. By interpreting values of surface tension measured at various surface ages, adsorption and diffusion rates of oxygen on the surface of the melt can be estimated. This research program investigates the adsorption kinetics of oxygen on the surface of an atomizing molten-metal jet. A novel oscillating capillary jet method has been developed for the measurement of dynamic surface tension of liquids, and in particular, metal melts which are susceptible to rapid surface degradation caused by oxygen adsorption. The experimental technique captures the evolution of jet swells and necks continuously along the jet propagation axis and is used in conjunction with an existing linear, axisymmetric, constant-property model to determine the variation of the instability growth rate, and, in turn, surface tension of the liquid as a function of surface age measured from the exit orifice. The conditions investigated so far focus on a time window of 2-4ms from the jet orifice. The surface properties of the eutectic 63%Sn-37%Pb solder alloy have been investigated in terms of their variation due to O2 adsorption from a N2 atmosphere containing controlled amounts of oxygen (from 8 ppm to 1000 ppm). The method performed well for situations where the oxygen adsorption was low in that time window. The value of surface tension for the 63Sn-37Pb solder in pure nitrogen was found to be 0.49 N/m, in good agreement with previously published work. A characteristic time of O(1ms) or less was determined for the molten-metal surface to be saturated by oxygen at 1000 ppm concentration in N2.

Megaridis, Constantine M.↗

Survey of the Relationship Between Theory and Experiment for Convective Heat Transfer From Rocket Combustion Gases

The problem of predicting the distribution of heat flux in a rocket engine is complicated by a flow field and thermodynamics so perturbed by as yet indescribable combustion effects that it continues to be necessary to utilize a simplified theoretical model against which to compare experimental results in order to gain insight. The most successful model applied to this problem is that of a compressible, turbulent, accelerating boundary layer flow with an arbitrary point of origin, and chemical reaction restricted to equilibrium recombination or diffusion-limited recombination on a catalytic wall. Analysis of this model was not possible without making several assumptions about the nature of the skin friction coefficient relationship and Reynold's analogy in an accelerating boundary layer flow. The paper summarizes the analysis of this model (together with detailed treatment of the special problems of variable properties, chemical reaction, and sources of transport properties), and describes several closed-form equations that have also been developed and applied to the problem. Described are results of several experiments with hot air flow made for the purpose of checking the validity of the assumptions in the boundary layer analysis. Most significant, however, are presentations of sample distributions of local heat flux measured in rocket engines operating with a wide range of conditions, propellants, injectors and sizes. When compared with predictions, these data show the significant range of deviation in heat flux that results from the combustion effects and which precludes a generally accurate prediction theory or correlation equation until the combustion process itself is correlated.

D R Bartz↗

Combustion of liquid-fuel droplets in supercritical conditions

A comprehensive analysis of liquid-fuel droplet combustion in both subcritical and supercritical environments has been conducted. The formulation is based on the complete conservation equations for both gas and liquid phases, and accommodates variable thermophysical properties, finite-rate chemical kinetics, and a full treatment of liquid-vapor phase equilibrium at the drop surface. The governing equations and associated interfacial boundary conditions are solved numerically using a fully coupled, implicit scheme with the dual time-stepping integration technique. The model is capable of treating the entire droplet history, including the transition from the subcritical to supercritical state. As a specific example, the combustion of n-pentane fuel droplets in air is studied for pressures in the range of 5-140 atm. Results indicate that the ambient gas pressure exerts significant control of droplet gasification and burning processes through its influence on fluid transport, gas-liquid interfacial thermodynamics, and chemical reactions. The droplet gasification rate increases progressively with pressure. However, the data for the overall burnout time exhibit a considerable change in the combustion mechanism at the critical pressure, mainly as a result of reduced mass diffusivity and latent heat of vaporization with increased pressure.

Shuen, J. S.↗

Investigation of turbulent diffusion in the extreme lower atmosphere

Turbulent diffusion in the extreme lower layer of the atmosphere (up to 5 m) has been investigated. Turbulent flow was simulated under dry, stable and calm conditions by means of a 3.04 m diameter fan installed at a field site situated on flat grassland. The ambient wind was continuously monitored by means of a cup anemometer placed outside the wake, and the temperature distribution was measured by four thermometers placed on an 18 m tower, also outside the wake. Balloons and red smoke were used to visualize the wake flow and investigate the predominant sizes of turbulent eddies and their streamwise behavior. The mean and turbulent velocities along the turbulence line were measured using an array of hot-wire anomometers. Results provide substantial verification of a recently proposed model (Koper and Sadeh, 1975; Koper et al., 1978) relating the Lagrangian to the Eulerian turbulent velocity autocorrelation. In this model the Lagrangian autocorrelation is given by a domain integral over a set of ordinary Eulerian autocorrelations acquired simultaneously at all points within the flow field in question, which is viewed as a turbulence 'box'.

Koper, C. A., Jr.↗

Comprehensive modeling of turbulent particulate flows using Eulerian and Lagrangian schemes

This paper addresses turbulent particle dispersion and modulation effects in dilute gas-particle turbulent flows using the Eulerian and Lagrangian modeling approaches. Gradient diffusion approximations are employed in the Eulerian formulation, while a stochastic procedure is utilized to simulate turbulent dispersion in the Lagrangian formulation. The k-epsilon turbulence model is used to characterize the time and length scales of the continuous phase turbulence. For the particle size and loading considered, the turbulence transport equations must be modified to account for the modulation effects. Models are proposed for both Eulerian and Lagrangian schemes. Comparisons and predictions are made in fully developed gas-solid pipe flow and confined coaxial jets laden with particles. For the monodispersed system investigated, Eulerian approach is less expensive and gives more consistent results than the Lagrangian approach. The Lagrangian technique should be further developed to eliminate current inherent inconsistencies especially with regard to the symmetry boundary condition.

Adeniji-Fashola, A.↗

Diffusion mechanisms in Ir-coated Re for high-temperature, radiation-cooled rocket thrusters

Materials used for radiation-cooled rocket thrusters must be capable of surviving under extreme conditions of high temperatures and oxidizing environments. Thruster chambers were developed using chemical-vapor-deposited (CVD) Re coated with CVD Ir on the inside surface which is exposed to hot combustion gases. Ir serves as an oxidation barrier protecting the Re which maintains structural integrity at high temperatures. In order to predict and extend the performance limits of these Ir-coated Re thrusters, the diffusion kinetics of CVD materials at temperature are studied. Thruster end ring sections were examined using electron microprobe analysis both before and after exposure to high temperature vacuum environments. The resulting elemental maps for Re, Ir, and Mo in the near-surface region allow identification of diffusion mechanisms operating at these temperatures. Line scans for Ir and Re were fit using a diffusion model to extract relevant diffusion constants. The fastest diffusion process is seen to be grain boundary diffusion with Re diffusing down grain boundaries in the Ir overlayer. The measured dependence of the diffusion rate on temperature will allow prediction of operating lifetimes for these thrusters.

Hamilton, J. C.↗