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.

SourceResearch connectionFrom yearThrough year

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

At least 1,261 records · Page 70

All sources NASA NTRS DOE OSTI

NASA NTRS1995

A turbulence module for the NPARC code

A turbulence module is developed for the 2D version of the NPARC code which is currently restricted to planar or axisymmetric flows without swirling. Four turbulence models have been built into the module: Baldwin-Lomax, Chien, Shih-Lumley and CMOTT models. The first is a mixing-length eddy-viscosity model which is mainly used for initialization of computational fields and the last three are the low Reynolds number two-equation models. Unlike chien's model, both the Shih-Lumley and CMOTT models do not involve the dimensionless wall distance y(sup +), an advantage for separated flow calculations. Contrary to the NPARC and most other compressible codes, the non-delta form of transport equations is used which leads to a simpler linearization and is more effective than using the delta form in ensuring the positiveness of the turbulent kinetic energy and its dissipation rate. To reduce numerical diffusion while maintaining necessary stability, a second-order accurate and bounded scheme is used for the convective terms of the turbulent transport equations. This scheme is implemented in a deferred correction manner so that the main coefficients of the resulting difference equations are always positive, thus making the numerical solutions process unconditionally stable. The system of equations is solved via a decoupled method and by the alternating direction TDMA of Thomas. The module can be easily linked to the NPARC code for turbulent flow calculations.

Zhu, J.↗

NASA NTRS1999

Efficient Multi-Stage Time Marching for Viscous Flows via Local Preconditioning

A new method has been developed to accelerate the convergence of explicit time-marching, laminar, Navier-Stokes codes through the combination of local preconditioning and multi-stage time marching optimization. Local preconditioning is a technique to modify the time-dependent equations so that all information moves or decays at nearly the same rate, thus relieving the stiffness for a system of equations. Multi-stage time marching can be optimized by modifying its coefficients to account for the presence of viscous terms, allowing larger time steps. We show it is possible to optimize the time marching scheme for a wide range of cell Reynolds numbers for the scalar advection-diffusion equation, and local preconditioning allows this optimization to be applied to the Navier-Stokes equations. Convergence acceleration of the new method is demonstrated through numerical experiments with circular advection and laminar boundary-layer flow over a flat plate.

Kleb, William L.↗

NASA NTRS1999

Efficient Multi-Stage Time Marching for Viscous Flows via Local Preconditioning

Kleb, William L.↗

NASA NTRS2019

Two-Stage Autoignition Dynamics of N-Dodecane Droplets Under Normal Gravity at High Pressures

Experimental observations of two-stage, autoignition characteristics of fiber-supported, n-dodecane fuel droplets under normal gravity are reported for a wide range of ambient temperatures, and pressures varying from 1 atm to 25 atm, which encompasses the fuel critical pressure (~ 18 atm). A high-speed camera is used to obtain backlit images of the autoignition process. During two-stage ignition, a cool-flame front forms below the droplet in the natural convective plume of the fuel vapor originating from the vaporizing droplet, and propagates upward, followed by the formation of a hot-flame kernel that expands rapidly and surrounds the droplet, eventually leading to the classical diffusion flame. Experimentally measured cool, and hot flame formation locations and propagation speeds are presented for varying initial temperatures and pressures. Measured first, and second stage ignition delay times are also reported. The influences of negative temperature coefficient (NTC) regime chemistry on flame starting location, propagation speeds, and ignition delay times are described. The first-stage cool flame characteristics are found to be controlled by both the chemistry, and physical aspects such as droplet evaporation and fuel-oxidizer mixing, while the second-stage hot flame ignition is primarily controlled by chemical kinetics. Further work is currently underway to develop models to quantitatively predict the observed behavior.

Rose, Evan↗

NASA NTRS2013

Direct Simulation of Extinction in a Slab of Spherical Particles

The exact multiple sphere superposition method is used to calculate the coherent and incoherent contributions to the ensemble-averaged electric field amplitude and Poynting vector in systems of randomly positioned nonabsorbing spherical particles. The target systems consist of cylindrical volumes, with radius several times larger than length, containing spheres with positional configurations generated by a Monte Carlo sampling method. Spatially dependent values for coherent electric field amplitude, coherent energy flux, and diffuse energy flux, are calculated by averaging of exact local field and flux values over multiple configurations and over spatially independent directions for fixed target geometry, sphere properties, and sphere volume fraction. Our results reveal exponential attenuation of the coherent field and the coherent energy flux inside the particulate layer and thereby further corroborate the general methodology of the microphysical radiative transfer theory. An effective medium model based on plane wave transmission and reflection by a plane layer is used to model the dependence of the coherent electric field on particle packing density. The effective attenuation coefficient of the random medium, computed from the direct simulations, is found to agree closely with effective medium theories and with measurements. In addition, the simulation results reveal the presence of a counter-propagating component to the coherent field, which arises due to the internal reflection of the main coherent field component by the target boundary. The characteristics of the diffuse flux are compared to, and found to be consistent with, a model based on the diffusion approximation of the radiative transfer theory.

spheres↗

NASA NTRS1996

Cross-Effects in Microgravity Flows

Film growth by chemical/physical vapor deposition is a process of considerable interest in microgravity experiments. The absence of natural convection should allow better control of film growth processes but, in highly non-isothermal ampoules, thermal slip (creep) can become a matter of significant concern. The reported research is a theoretical and experimental investigation of the flow of gas/vapor mixtures under non-continuum conditions. The Boltzmann equation has been solved for a monatomic gas under non-condensing conditions and the various phenomenological coefficients have been computed. Computations for realistic potentials as well as for velocity and creep slip have been completed and the creep slip has been found to be dependent on the type of gas confirming the accuracy of previous variational results. The variational technique has been extended and planar flows calculated via the Burnett solutions. Velocity, diffusion and creep slips have been computed for gas mixtures and previously unknown dependencies of the creep slip on the mixture properties have been observed. Also for gas mixtures, an integral representation of the linearized Boltzmann operator has been developed for use in numerical and variational calculations for all intermolecular force laws. Two, two-bulb capillary systems have been designed, built and tested for the measurements of cross-flows; one of glass for isothermal measurements and one of stainless steel for non-isothermal measurements. Extensive data have been collected for Ar-He and N2-He mixtures at a variety of pressures and mole ratios. Viscosity, velocity slip coefficients and tangential momentum accommodation coefficients have been obtained from measurements with a spinning rotor gauge via a new theory that has been formulated for the spinning rotor gauge in the slip regime. The FIDAP fluid dynamics code has been applied to condensing flows in ampoules in the continuum regime and agreement obtained with the earlier work of Duval.

Loyalka, Sudarshan K.↗

NASA NTRS2023

Numerical Investigation of Film Coefficient Approximation for Chemically Reacting Boundary-Layer Flows

Aerothermal analysis of spacecraft planetary entry is heavily dependent on heritage engineering models. The film coefficient heat transfer model examined in this paper estimates the convective heating to the vehicle for a laminar, dissociated, chemically reacting boundary layer for an Earth atmosphere. This model requires information about the vehicle and flowfield for a given trajectory point and estimates a proportional relationship between enthalpy potential and convective heat flux. In practice it is the aerothermal engineer who must decide which assumptions are appropriate for his/her application. This work looks at numerous CFD simulations for an arbitrary, axisymmetric flight vehicle to analyze the relative importance of both the mass and energy constraints imposed at the wall boundary, as well as the effect of various diffusion models. Within the subset of tested energy boundary conditions, it is found that the most desirable energy boundary condition is the radiative equilibrium boundary condition, which permits conservative estimates of convective heat flux, but also generates flowfield-dependent spatial thermal distributions along the surface. Other key findings are presented in an effort to make the film coefficient engineering model readily available to design engineers across industry.

Film coefficient↗

NASA NTRS1992

Impeller tandem blade study with grid embedding for local grid refinement

Flow non-uniformity at the discharge of high power density impellers can result in significant unsteady interactions between impeller blades and downstream diffuser vanes. These interactions result in degradation of both performance and pump reliability. The MSFC Pump Technology Team has recognized the importance of resolving this problem and has thus initiated the development and testing of a high head coefficient impeller. One of the primary goals of this program is to improve impeller performance and discharge flow uniformity. The objective of the present work is complimentary. Flow uniformity and performance gains were sought through the application of a tandem blade arrangement. The approach adopted was to numerically establish flow characteristics at the impeller discharge for the baseline MSFC impeller and then parametrically evaluate tandem blade configurations. A tandem design was sought that improves both impeller performance and discharge uniformity. The Navier-Stokes solver AEROVISC was used to conduct the study. Grid embedding is used to resolve local gradients while attempting to minimize model size. Initial results indicate that significant gains in flow uniformity can be achieved through the tandem blade concept and that blade clocking rather than slot location is the primary driver for flow uniformity.

Bache, George↗

NASA NTRS2024

Toward 2D Dynamo Models Calibrated by Global 3D Relativistic Accretion Disk Simulations

Two-dimensional models assuming axisymmetry are an economical way to explore the long-term evolution of black hole accretion disks, but they are only realistic if the feedback of the nonaxisymmetric turbulence on the mean momentum and magnetic fields is incorporated. Dynamo terms added to the 2D induction equation should be calibrated to 3D MHD simulations. For generality, the dynamo tensors should be calibrated as functions of local variables rather than explicit functions of spatial coordinates in a particular basis. In this paper, we study the feedback of non-axisymmetric features on the 2D mean fields using a global 3D, relativistic, Cartesian simulation from the Illinois GRMHD code. We introduce new methods for estimating overall dynamo alpha and turbulent diffusivity effects as well as measures of the dominance of non-axisymmetric components of energies and fluxes within the disk interior. We attempt closure models of the dynamo EMF using least squares fitting, considering both models where coefficient tensors are functions of space and more global, covariant models. None of these models are judged satisfactory, but we are able to draw conclusions on what sorts of generalizations are and are not promising.

Matthew D Duez↗

NASA NTRS1972

Analytical correlation of centrifugal compressor design geometry for maximum efficiency with specific speed

Centrifugal compressor performance was examined analytically to determine optimum geometry for various applications as characterized by specific speed. Seven specific losses were calculated for various combinations of inlet tip-exit diameter ratio, inlet hub-tip diameter ratio, blade exit backsweep, and inlet-tip absolute tangential velocity for solid body prewhirl. The losses considered were inlet guide vane loss, blade loading loss, skin friction loss, recirculation loss, disk friction loss, vaneless diffuser loss, and vaned diffuser loss. Maximum total efficiencies ranged from 0.497 to 0.868 for a specific speed range of 0.257 to 1.346. Curves of rotor exit absolute flow angle, inlet tip-exit diameter ratio, inlet hub-tip diameter ratio, head coefficient and blade exit backsweep are presented over a range of specific speeds for various inducer tip speeds to permit rapid selection of optimum compressor size and shape for a variety of applications.

Galvas, M. R.↗

NASA NTRS1983

Tribological and microstructural characteristics of ion-nitrided steels

Three steels AISI 4140, AISI 4340 and AISI 304 stainless steel were ion nitrided in a plasma consisting of a 75:25 mixture of H2:N2, sometimes with a trace of CH4. Their surface topography was characterized by SEM and two distinct compound phases were identified: the gamma and the epsilon. The core-case hardness profiles were also established. The low Cr alloy steels have an extended diffusion zone in contrast to the 304 stainless steels which have a sharp interface. The depth of ion-nitriding is increased as the Cr content is decreased. Friction tests reveal that the gamma surface phase has a lower coefficient of friction than the epsilon phase. The lowest coefficient of friction is achieved when both the rider and the specimen surface are ion nitrided.

Spalvins, T.↗

NASA NTRS1983

Tribological and microstructural characteristics of ion-nitrided steels

Three steels AISI 4140, AISI 4340 and AISI 304 stainless steel were ion nitrided in a plasma consisting of a 75:25 mixture of H2:N2, sometimes with a trace of CH4. Their surface topography was characterized by SEM and two distinct compound phases were identified: the gamma and the epsilon. The core-case hardness profiles were also established. The low Cr alloy steels have an extended diffusion zone in contrast to the 3034 stainless steels which have a sharp interface. The depth of ion-nitriding is increased as the Cr content is decreased. Friction tests reveal that the gamma surface phase has a lower coefficient of friction than the epsilon phase. The lowest coefficient of friction is achieved when both the rider and the specimen surface are ion nitrided. Previously announced in STAR as N83-24635

Spalvins, T.↗

NASA NTRS1985

Frictional and structural characterization of ion-nitrided low and high chromium steels

Low Cr steels AISI 41410, AISI 4340, and high Cr austenitic stainless steels AISI 304, AISI 316 were ion nitrided in a dc glow discharge plasma consisting of a 75 percent H2 - 25 percent N2 mixture. Surface compound layer phases were identified, and compound layer microhardness and diffusion zone microhardness profiles were established. Distinct differences in surface compound layer hardness and diffusion zone profiles were determined between the low and high Cr alloy steels. The high Cr stainless steels after ion nitriding displayed a hard compound layer and an abrupt diffusion zone. The compound layers of the high Cr stainless steels had a columnar structure which accounts for brittleness when layers are exposed to contact stresses. The ion nitrided surfaces of high and low Cr steels displayed a low coefficient of friction with respect to the untreated surfaces when examined in a pin and disk tribotester.

Spalvins, T.↗

NASA NTRS1985

Frictional and structural characterization of ion-nitrided low and high chromium steels

Spalvins, T.↗

NASA NTRS2014

Wind-US Results for the AIAA 2nd Propulsion Aerodynamics Workshop

This presentation contains Wind-US results presented at the 2nd Propulsion Aerodynamics Workshop. The workshop was organized by the American Institute of Aeronautics and Astronautics, Air Breathing Propulsion Systems Integration Technical Committee with the purpose of assessing the accuracy of computational fluid dynamics for air breathing propulsion applications. Attendees included representatives from government, industry, academia, and commercial software companies. Participants were encouraged to explore and discuss all aspects of the simulation process including the effects of mesh type and refinement, solver numerical schemes, and turbulence modeling. The first set of challenge cases involved computing the thrust and discharge coefficients for a 25deg conical nozzle for a range of nozzle pressure ratios between 1.4 and 7.0. Participants were also asked to simulate two cases in which the 25deg conical nozzle was bifurcated by a solid plate, resulting in vortex shedding (NPR=1.6) and shifted plume shock (NPR=4.0). A second set of nozzle cases involved computing the discharge and thrust coefficients for a convergent dual stream nozzle for a range of subsonic nozzle pressure ratios. The workshop committee also compared the plume mixing of these cases across various codes and models. The final test case was a serpentine inlet diffuser with an outlet to inlet area ratio of 1.52 and an offset of 1.34 times the inlet diameter. Boundary layer profiles, wall static pressure, and total pressure at downstream rake locations were examined.

Wind-US↗

NASA NTRS1977

Miniature drag force anemometer

A miniature drag force anemometer is described which is capable of measuring dynamic velocity head and flow direction. The anemometer consists of a silicon cantilevered beam 2.5 mm long, 1.5 mm wide, and 0.25 mm thick with an integrated diffused strain gage bridge, located at the base of the beam, as the force measuring element. The dynamics of the beam are like that of a second order system with a natural frequency of about 42 kHz and a damping coefficient of 0.007. The anemometer can be used in both forward and reversed flow. Measured flow characteristics up to Mach 0.6 are presented along with application examples including turbulence measurements.

Krause, L. N.↗

NASA NTRS1983

Cyclic oxidation of superalloys

The development of mathematical/computer models of cyclic oxidation in high-temperature alloys is described. The acceleration of oxidation attack caused by the spalling of protective oxide scales is discussed. The COREST paralinear model based on empirical data, designed by Barrett and Pressler, is found to be of limited value. Basic research on the spalling process in different alloy systems is reviewed, and the phenomenon is characterized: occurring mainly within the oxide, it is caused by the different thermal expansion coefficients of oxide and metal and is a linear function of the total retained oxide. This information has been used in the development of the discrete oxide-spall model COSP and in its modification by addition of a diffusion-zone-profile model to predict the failure of components (Heckel and Nesbitt, unpublished). The COSP model is shown to produce good predictions of cyclic oxidation under a variety of conditions, using data from simple isothermic tests.

Lowell, C. E.↗

NASA NTRS1973

Viscous flow over spinning cones at angle of attack.

A numerical finite-difference method is developed for evaluating the Magnus coefficients on spinning cones in laminar flow. The merged layer, the strong interaction region, and the downstream boundary layer are all considered. The numerical method is a predictor-corrector scheme developed for three-dimensional flows with or without crossflow diffusion. This method is particularly useful in problems in which a symmetry plane does not exist. Several contributions to the Magnus force and moments are considered. These include asymmetries in displacement thickness, centrifugal force and crossflow shear, and the effects of crossflow separation and vortex formation. Comparisons are made with experimental data and other analyses.

Lin, T. C.↗