Search NASA⌕ Search

SEARCH · Search NASA

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 487 records · Page 27

Progress of High Efficiency Centrifugal Compressor Simulations Using TURBO

Three-dimensional, time-accurate, and phase-lagged computational fluid dynamics (CFD) simulations of the High Efficiency Centrifugal Compressor (HECC) stage were generated using the TURBO solver. Changes to the TURBO Parallel Version 4 source code were made in order to properly model the no-slip boundary condition along the spinning hub region for centrifugal impellers. A startup procedure was developed to generate a converged flow field in TURBO. This procedure initialized computations on a coarsened mesh generated by the Turbomachinery Gridding System (TGS) and relied on a method of systematically increasing wheel speed and backpressure. Baseline design-speed TURBO results generally overpredicted total pressure ratio, adiabatic efficiency, and the choking flow rate of the HECC stage as compared with the design-intent CFD results of Code Leo. Including diffuser fillet geometry in the TURBO computation resulted in a 0.6 percent reduction in the choking flow rate and led to a better match with design-intent CFD. Diffuser fillets reduced annulus cross-sectional area but also reduced corner separation, and thus blockage, in the diffuser passage. It was found that the TURBO computations are somewhat insensitive to inlet total pressure changing from the TURBO default inlet pressure of 14.7 pounds per square inch (101.35 kilopascals) down to 11.0 pounds per square inch (75.83 kilopascals), the inlet pressure of the component test. Off-design tip clearance was modeled in TURBO in two computations: one in which the blade tip geometry was trimmed by 12 mils (0.3048 millimeters), and another in which the hub flow path was moved to reflect a 12-mil axial shift in the impeller hub, creating a step at the hub. The one-dimensional results of these two computations indicate non-negligible differences between the two modeling approaches.

turbomachinery↗

A quasi-linear kinetic equation for cosmic rays in the interplanetary medium

A kinetic equation for interplanetary cosmic rays is set up with the aid of weak-plasma-turbulence theory for an idealized radially symmetric model of the interplanetary magnetic field. As a starting point, this treatment invokes the Vlasov equation instead of the traditional Fokker-Planck equation. Quasi-linear theory is applied to obtain a momentum diffusion equation for the heliocentric frame of reference which describes the interaction of cosmic rays with convecting magnetic irregularities in the solar-wind plasma. Under restricted conditions, the well-known equation of solar modulation can be obtained from this kinetic equation.

Luhmann, J. G.↗

Thermofluid analysis of the SSME preburner using a gas-gas diffusion model for oxygen and hydrogen combustion at supercritical pressures

The paper discusses the thermofluid analysis of the Space Shuttle Main Engine (SSME) fuelside preburner. The governing equations have been solved numerically to predict flow, heat transfer, mixing, and combustion. A two-fluid approach is adopted in which oxygen is regarded as one fluid and hydrogen is regarded as the other fluid. The chemical kinetics is assumed to be very fast so that combustion is primarily controlled by the rate of mixing between oxygen and hydrogen. The preburner pressure is much greater than the critical pressures of oxygen and hydrogen; hence, a gas-gas diffusion model (rather than an evaporation model) has been developed to compute the rate of interphase mixing. Empirical correlations have been incorporated to account for the effect of slip on the interphase exchange. A sensitivity study has been performed with various model parameters. It is observed that the model can predict possibility of incomplete combustion and local regions of high temperatures under steady operating conditions. Some of these anomalies have been observed in actual tests, and the numerical model is useful for understanding possible causes and remedies. At least some measurements are needed for quantitative verification of the model.

Prakash, C.↗

Computation of diffuse sky irradiance from multidirectional radiance measurements

Accurate determination of the diffuse solar spectral irradiance directly above the land surface is important in characterizing the reflectance properties of these surfaces, especially vegetation canopies. This determination is also needed to infer the net radiation budget of the earth-atmosphere system above these surfaces. An algorithm is developed here for the computation of hemispheric diffuse irradiance using the measurements from an instrument called PARABOLA, which rapidly measures upwelling and downwelling radiances in three selected wavelength bands. The validity of the algorithm is established from simulations. The standard reference data set of diffuse radiances of Dave (1978), obtained by solving the radiative transfer equation numerically for realistic atmospheric models, is used to simulate PARABOLA radiances. Hemispheric diffuse irradiance is estimated from a subset of simulated radiances by using the algorithm described. The algorithm is validated by comparing the estimated diffuse irradiance with the true diffuse irradiance of the standard data set. The validations include sensitivity studies for two wavelength bands (visible, 0.65-0.67 micron; near infrared, 0.81-0.84 micron), different atmospheric conditions, solar elevations, and surface reflectances. In most cases the hemispheric diffuse irradiance computed from simulated PARABOLA radiances and the true irradiance obtained from radiative transfer calculations agree within 1-2 percent. This technique can be applied to other sampling instruments designed to estimate hemispheric diffuse sky irradiance.

Ahmad, Suraiya P.↗

Effects of buoyancy-driven flow and thermal boundary conditions on physical vapor transport

A 2D numerical model was developed in order to ascertain if reduced gravity conditions are beneficial to physical vapor transport (PVT) and to determine its tolerance limits to residual accelerations. This was solved using the PHOENICS finite-volume code. Reduction of gravitational accelerations to less than 0.1 g0 was found to be sufficient to suppress buoyancy-driven convection to an extent that diffusion was the dominant transport mode, whence a greater uniformity in the growth rate could be obtained. It is shown that a uniform temperature gradient on the ampoule walls causes the vapor to be supersaturated throughout the ampoule, potentially resulting in undesirable nucleation at the walls. A 'hump' in the wall temperature profile can be used to avoid this. The prevailing transport conditions determine the size of the hump needed.

Nadarajah, Arunan↗

The Effects of Buoyancy and Dilution on the Structure and Lift-off of Coflow Laminar Diffusion Flames

The ability to predict the coupled effects of complex transport phenomena with detailed chemical kinetics in diffusion flames is critical in the modeling of turbulent reacting flows and in understanding the processes by which soot formation and radiative transfer take place. In addition, an understanding of the factors that affect flame extinction in diffusion flames is critical in the suppression of fires and in improving engine efficiency. The goal of our characterizations of coflow laminar diffusion flames is to bring to microgravity the multidimensional diagnostic tools available in normal gravity, and in so doing provide a broader understanding of the successes and limitations of current combustion models. This will lead to a more detailed understanding of the interaction of convection, diffusion and chemistry in both buoyant and nonbuoyant environments. As a sensitive marker of changes in the flame shape, the number densities of excited-state CH (A(exp 2)delta, denoted CH*), and excited-state OH (A(exp 2)Sigma, denoted OH*) are measured in mu-g and normal gravity. Two-dimensional CH* and OH* number densities are deconvoluted from line-of-sight chemiluminescence measurements made on the NASA KC-135 reduced-gravity aircraft. Measured signal levels are calibrated, post-flight, with Rayleigh scattering. Although CH* and OH* kinetics are not well understood, the CH*, OH*, and ground-state CH distributions are spatially coincident in the flame anchoring region. Therefore, the ground-state CH distribution, which is easily computed, and the readily measured CH*/OH* distributions can be used to provide a consistent and convenient way of measuring lift-off height and flame shape in the diffusion flame under investigation. Given that the fuel composition affects flame chemistry and that buoyancy influences the velocity profile of the flow, we have the opportunity to computationally and experimentally study the roles of fluids and chemistry. In performing this microgravity study, improvements to the computational model have been made and new calculations performed for a range of gravity and flow conditions. Furthermore, modifications to the experimental approach were required as a consequence of the constraints imposed by existing microgravity facilities. Results from the computations and experiments are presented in the following sections.

Walsh, Kevin T.↗

The Role of Kinetic Processes in Magnetospheric Dynamics

The question of whether the micro scale controls the macro scale or vice-versa remains one of the most challenging problems in plasmas. A particular topic of interest within this context is collisionless magnetic reconnection, where both points of views are espoused by different groups of researchers. This presentation will focus on this topic. We will begin by analyzing the global properties of magnetic reconnection. Following that, we discuss properties of electron diffusion region dynamics both for guide field and antiparallel reconnection, and how they can be scaled to different inflow conditions. As a next step, we will study typical temporal variations of the microscopic dynamics with the objective of understanding the potential for secular changes to the macroscopic system. The research will be based on a combination of analytical theory and numerical modeling.

Hesse, Michael↗

Modeling of a Non-Dilute Alloy Solidification Under Terrestrial and Microgravity Conditions

A general finite element model of non-dilute alloy solidification was used to simulate growth of HgCdTe from the melt. 2-D and 3-D models with and without Oberbeck-Boussinesq (O-B) approximation were used. Verification of the 3-D model was based on experimental data. The location and the shape of the interface were determined using melting temperatures obtained from the HgCdTe liquidus curve. Double-diffusive convection in the melt is caused by radial temperature gradients and by rejection at the solid-melt interface of the heavier HgTe-rich solute. Because of these complicating factors, dimensional rather than non-dimensional modeling was performed. Models based on complete set of Navier-Stokes equations for compressible media are superior to the O-B approximation in microgravity, because they provide more accurate estimation of convection.

Bune, Andris V.↗

Airfoil deposition model

The methodology to predict deposit evolution (deposition rate and subsequent flow of liquid deposits) as a function of fuel and air impurity content and relevant aerodynamic parameters for turbine airfoils is developed in this research. The spectrum of deposition conditions encountered in gas turbine operations includes the mechanisms of vapor deposition, small particle deposition with thermophoresis, and larger particle deposition with inertial effects. The focus is on using a simplified version of the comprehensive multicomponent vapor diffusion formalism to make deposition predictions for: (1) simple geometry collectors; and (2) gas turbine blade shapes, including both developing laminar and turbulent boundary layers. For the gas turbine blade the insights developed in previous programs are being combined with heat and mass transfer coefficient calculations using the STAN 5 boundary layer code to predict vapor deposition rates and corresponding liquid layer thicknesses on turbine blades. A computer program is being written which utilizes the local values of the calculated deposition rate and skin friction to calculate the increment in liquid condensate layer growth along a collector surface.

Kohl, F. J.↗

Simulation of SEU Cross-sections using MRED under Conditions of Limited Device Information

This viewgraph presentation reviews the simulation of Single Event Upset (SEU) cross sections using the membrane electrode assembly (MEA) resistance and electrode diffusion (MRED) tool using "Best guess" assumptions about the process and geometry, and direct ionization, low-energy beam test results. This work will also simulate SEU cross-sections including angular and high energy responses and compare the simulated results with beam test data for the validation of the model. Using MRED, we produced a reasonably accurate upset response model of a low-critical charge SRAM without detailed information about the circuit, device geometry, or fabrication process

Lauenstein, J. M.↗

Two phase choke flow in tubes with very large L/D

Data were obtained for two phase and gaseous choked flow nitrogen in a long constant area duct of 16200 L/D with a diverging diffuser attached to the exit. Flow rate data were taken along five isotherms (reduced temperature of 0.81, 0.96, 1.06, 1.12, and 2.34) for reduced pressures to 3. The flow rate data were mapped in the usual manner using stagnation conditions at the inlet mixing chamber upstream of the entrance length. The results are predictable by a two phase homogeneous equilibrium choking flow model which includes wall friction. A simplified theory which in essence decouples the long tube region from the high acceleration choking region also appears to predict the data resonably well, but about 15 percent low.

Hendricks, R. C.↗

Simulation of the thermospheric tides by use of the NCAR thermospheric general circulation model

Numerical calculations of the thermospheric tidal winds and temperatures at equinox are presented. The calculations were made using the National Center for Atmospheric Research Thermospheric General Circulation Model (TGCM), which includes the effects of viscosity, conductivity, diffusion, ion drag, winds, and temperature gradients. The thermospheric diurnal and semidiurnal tides are excited in situ by solar heating and by ion-neutral momentum coupling. The semidiurnal tidal calculations also include the effects of upward propagating waves generated by heating in the lower atmosphere. This semidiurnal propagating component is modeled by use of the classical tidal perturbations as lower boundary conditions. The model is tuned by adjusting the propagating tidal forcing term until calculated semidiurnal wind and temperature fields best approximate incoherent scatter observations. The tidal TGCM results are consistent with previous theoretical work and successfully reproduce high altitude temperature and meridional velocity data, but they give significantly lower magnitudes for velocities and temperatures near 160 km than are seen by observations.

Fensen, C. G.↗

Modeling of outgassing and matrix decomposition in carbon-phenolic composites

A new release rate equation to model the phase change of water to steam in composite materials was derived from the theory of molecular diffusion and equilibrium moisture concentration. The new model is dependent on internal pressure, the microstructure of the voids and channels in the composite materials, and the diffusion properties of the matrix material. Hence, it is more fundamental and accurate than the empirical Arrhenius rate equation currently in use. The model was mathematically formalized and integrated into the thermostructural analysis code CHAR. Parametric studies on variation of several parameters have been done. Comparisons to Arrhenius and straight-line models show that the new model produces physically realistic results under all conditions.

Mcmanus, Hugh L.↗

Two phase choke flow in tubes with very large L/D

Two phase and gaseous choked flow data for fluid nitrogen were obtained for a test section which was a long constant area duct of 16 200 L/D with a diverging diffuser attached to the exit. Flow rate data were taken along five isotherms (reduced temperature of 0.81, 0.96, 1.06, 1.12, and 2.34) for reduced pressures to 3. The flow rate data were mapped in the usual manner using stagnation conditions at the inlet mixing chamber upstream of the entrance length. The results are predictable by a two-phase homogeneous equilibrium choking flow model which includes wall fraction. A simplified theory which in essence decouples the long tube region from the high acceleration choking region also appears to predict the data reasonably well, but about 15 percent low.

Hendricks, R. C.↗

Stratospheric eddy diffusion coefficients from tracer data

Global distributions of nitrous oxide, methane, ozone, and carbon 14 are used to estimate four sets of stratospheric eddy diffusion coefficients. A photochemical equilibrium model calculates O(3P), O(1D), H, HO2, OH, H2O2, NO, and NO2 densities, as a function of altitude, latitude, and time. The calculated O(1D), OH, and observed Cl densities are used to obtain the eddy profiles associated with the methane and nitrous oxide distributions, for altitudes between 10 and 40 km. Application of a constant flux condition to the seasonally averaged ozone data yields eddy values below 20 km. Time-dependent carbon 14 calculations produce eddy coefficients between 13 and 27 km. A composite profile is obtained by comparing the four sets of coefficients. Further, carbon 14 computations are used to test these profiles as well as those recommended in reports issued by the National Academy of Sciences in 1976 and 1979. The composite eddy profile produces the best agreement.

Massie, S. T.↗

Chemoviscosity modeling for thermosetting resins - I

A new analytical model for chemoviscosity variation during cure of thermosetting resins was developed. This model is derived by modifying the widely used WLF (Williams-Landel-Ferry) Theory in polymer rheology. Major assumptions involved are that the rate of reaction is diffusion controlled and is linearly inversely proportional to the viscosity of the medium over the entire cure cycle. The resultant first order nonlinear differential equation is solved numerically, and the model predictions compare favorably with experimental data of EPON 828/Agent U obtained on a Rheometrics System 4 Rheometer. The model describes chemoviscosity up to a range of six orders of magnitude under isothermal curing conditions. The extremely non-linear chemoviscosity profile for a dynamic heating cure cycle is predicted as well. The model is also shown to predict changes of glass transition temperature for the thermosetting resin during cure. The physical significance of this prediction is unclear at the present time, however, and further research is required. From the chemoviscosity simulation point of view, the technique of establishing an analytical model as described here is easily applied to any thermosetting resin. The model thus obtained is used in real-time process controls for fabricating composite materials.

Hou, T. H.↗

3-D Modeling of Directional Solidification of a Non-Dilute Alloy with Temperature and Concentration Fields Coupling via Materials Properties Dependence and via Double Diffusive Convection

Numerical simulation of the HgCdTe growth by the vertical Bridgman method was performed using FIDAP finite element code. Double-diffusive melt convection is analyzed, as the primary factor at controls inhomogeneity of the solidified material. Temperature and concentration fields in the model are also coupled via material properties, such as thermal and solutal expansion coefficients with the dependence on both temperature and concentration, and melting temperature evaluation from pseudobinary CdTe-HgTe phase diagram. Experimental measurements were used to obtain temperature boundary conditions. Parametric study of the melt convection dependence on the gravity conditions was undertaken. It was found, that the maximum convection velocity in the melt can be reduced under certain conditions. Optimal conditions to obtain a near flat solidified interface are discussed. The predicted interface shape is in agreement with one obtained experimentally by quenching. The results of 3-D calculations are compared with previous 2- D findings. A video film featuring 3-D melt convection will be presented.

Bune, Andris V.↗

Multi-Scale Physics in the Magnetosphere: The Role of Magnetic Reconnection

The question of whether the microscale controls the macroscale or vice-versa remains one of the most challenging problems in plasmas. A particular topic of interest within this context is collision less magnetic reconnection, where both pOints of views are espoused by different groups of researchers. This presentation will focus on this topic. We will begin by analyzing the properties of electron diffusion region dynamics both for guide field and anti-parallel reconnection, and how they can be scaled to different inflow conditions. As a next step, we will study typical temporal variations of the microscopic dynamics with the objective of understanding the potential for secular changes to the macroscopic system. The research will be based on a combination of analytical theory and numerical modeling.

Hesse, Michael↗