Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solids flow”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,639 records · Page 91

Incompressible viscous flow simulations of the NFAC wind tunnel

The capabilities of an existing 3-D incompressible Navier-Stokes flow solver, INS3D, are extended and improved to solve turbulent flows through the incorporation of zero- and two-equation turbulence models. The two-equation model equations are solved in their high Reynolds number form and utilize wall functions in the treatment of solid wall boundary conditions. The implicit approximate factorization scheme is modified to improve the stability of the two-equation solver. Applications to the 3-D viscous flow inside the 80 by 120 feet open return wind tunnel of the National Full Scale Aerodynamics Complex (NFAC) are discussed and described.

Champney, Joelle Milene↗

Direct simulation of compressible wall-bounded turbulence

When analyzing many turbulent flows, the effects of compressibility can be neglected. Even some relatively high-speed flows, such as boundary layers generated by a supersonic aircraft, produce turbulent statistics that are similar to those found for the incompressible case. There are other situations, however, in which the non-zero divergence of the turbulence leads to behavior that is fundamentally different from that found at constant density. Examples include flows created by internal combustion engines, hypersonic flight, and supersonic combustion. It is with instances such as these that this project is concerned. In particular, we are interested in the effects of compressibility on turbulence near a smooth solid constant temperature surface; our primary objective is an increased physical understanding that can be used to improve turbulence models of wall-bounded compressible flows. With this in mind, we have begun a direct numerical simulation (DNS) study of turbulence in a plane channel. Because all of the relevant spatial and temporal scales are to be resolved, the simulations require no subgrid scale parameterization. The DNS code developed by Buell to study compressible plane Couette flow was modified to solve the compressible Navier-Stokes equations in the plane channel. The channel was chosen over the Couette flow for two reasons: (1) to avoid using the very large streamwise domains needed to adequately capture the large Couette vortical structures; and (2) to make use of previous experience by considering the compressible version of a well established case and isolate finite Mach number effects by comparing to the incompressible channel.

Coleman, Gary N.↗

Simulation of transonic viscous flow over a fighter-like configuration including inlet

The simulation of transonic viscous flow over a modified F-16A including inlet is presented. A zonal approach is utilized which allows appropriate clustering suitable for viscous calculations on all solid surfaces. Computational efficiency is enhanced by solving the thin-layer Navier-Stokes equations in viscous zones adjacent to the aircraft and the Euler equations on those coarse zones away from the aircraft. The flow conditions for this transonic case are freestream Mach number of 0.9, angle of attack of 4.12 deg, and a Reynolds number (based on root chord) of 4.5 million. A total of 19 zones are utilized yielding a total of 350,000 grid points. This case required about 3000 iterations to reduce the residual by three orders, which takes about 10 hr of CPU time on the Cray X-MP/48 computer. Pressure distributions on the wing and on cross sections through the inlet region compare favorably with the experimental data for this transonic case.

Flores, Jolen↗

A fourth-order Lagrangian discontinuous Galerkin method using a hierarchical orthogonal basis on curvilinear grids

The existing high-order Lagrangian discontinuous Galerkin (DG) hydrodynamic methods are restricted to using quadratic meshes with quadratic polynomials (P2), which in turn, yield up to third-order accuracy. Additionally, these existing DG hydrodynamic schemes, when extended to work with cubic meshes and cubic polynomials (P3), can be unstable on strong-shock problems. Therefore, this paper presents a new Lagrangian discontinuous Galerkin (DG) hydrodynamic method to simulate compressible material dynamics ( e.g., gasses, fluids, and solids) with strong-shocks using cubic meshes and cubic polynomials, and delivers up to fourth-order accuracy on smooth flows. The stability on shock problems is achieved using new hierarchical orthogonal basis functions and a new subcell mesh stabilization (SMS) scheme for cubic meshes. The accuracy and robustness of the new high-order accurate Lagrangian DG hydrodynamic method is demonstrated by simulating a diverse suite of challenging test problems covering gas and solid dynamic problems on curvilinear meshes.

97 MATHEMATICS AND COMPUTING↗

Reactive Hydrodynamics in Rotating Spherical and Cylindrical Geometry

In turbulent spray combustion among many complex interactions between local flow structures called turbulent eddies and droplets are those associated with rotation of droplets. In general, for a complete statistical description of turbulent sprays, consideration of at least four degrees of freedom respectively associated with translational, rotational, vibrational (pulsational), and internal motions of the droplet are needed. Clearly the interactions between all degrees of freedom of the droplets and those for the gaseous background field will be exceedingly complex. For example, one type of interaction between the translational and the rotational velocity of droplets results in droplet helicity, H(d) = w(d).v(d), the significance of which in turbulent spray combustion is yet to be recognized. The role of droplet rotation in turbulent spray combustion modeling and its impact on the evaporation of liquid fuel droplets was recently investigated. Also, the impact of rotation on combustion of solid particles such as is encountered in pulverized coal combustion has been emphasized. The problem of viscous flow around a rotating sphere discussed above also occurs in other areas of physical sciences such as astrophysics and geophysics. Consequently, the subject has been addressed in many classical as well as more recent investigations. According to these investigations, the rotation of a rigid sphere in an otherwise quiescent, unconfined environment results in the motion of the fluid towards the poles. The polar flows from the northern and southern hemispheres move along helical trajectories towards the equatorial plane. Eventually, the polar flows collide at the equatorial plane, thus producing a sheet of rotating fluid that is radially ejected outward on this plane. Therefore, a droplet induces a strained flow field as a result of its rotation. Since the spatial extent of equatorial jets could easily exceed many droplet diameters, interactions between neighboring droplets are enhanced as a result of their rotation. Also, the equatorial jet substantially alters the spherical geometry of the diffusion flame surface that surrounds a rotating droplet. The objective of the research is to gain more knowledge about the hydrodynamics within and around rotating spherical and cylindrical body of fluid, and the behavior of diffusion or premixed flame surfaces that could surround such symmetric body of rotating fluid.

Sohrab, Siavash H.↗

A Numerical Method for Incompressible Flow with Heat Transfer

A numerical method for the convective heat transfer problem is developed for low speed flow at mild temperatures. A simplified energy equation is added to the incompressible Navier-Stokes formulation by using Boussinesq approximation to account for the buoyancy force. A pseudocompressibility method is used to solve the resulting set of equations for steady-state solutions in conjunction with an approximate factorization scheme. A Neumann-type pressure boundary condition is devised to account for the interaction between pressure and temperature terms, especially near a heated or cooled solid boundary. It is shown that the present method is capable of predicting the temperature field in an incompressible flow.

Sa, Jong-Youb↗

Thermal interaction between an impinging hot jet and a conducting solid surface

Powered-lift aircraft may produce severe high-temperature environments which are potentially damaging to a landing surface or the aircraft. The interaction betweean the high temperature flow field and a nonadiabatic landing surface is analyzed with a coupled computational fluid dynamics/solid thermal conduction computer code, HOTJET. The HOTJET code couples time-accurate, implicit, factored solution schemes for the governing fluid dynamics equations (Reynolds-averaged Navier-Stokes equations) to the unsteady thermal conduction equation, which governs heat flux within a solid. HOTJET is validated against exact solutions to the thermal conduction and Navier-Stokes equations. First-of-a-kind results are included which show the impact of surface material properties on the fluid physics and the coupled fluid/material thermal fields.

Abeloff, P. A.↗

Extinction of the Stagnation Point Diffusion Flame: Effect of Conductive Heat Loss into Solid Interior

The flammability boundary of a solid is U-shaped with the ordinate as ambient oxygen percentage (or total pressure) and the abscissa represents flow velocity, gravity level, stretch rate, etc. The left side of the flammability boundary is the quenching branch, and the right side is the blow-off branch. Quenching is due to a weak flame whose heat release from combustion cannot overcome the excessive heat losses from, for example, radiation. Blow-off is due to a short gas residence time that is smaller than the chemical reaction time (traditionally referred to as the small Damköhler number effect). For thick solids, the heat conduction into the solid interior has a long transient after ignition before the solid temperature reaches the steady state. During this slow transient, there is excessive (i.e. more than that at the steady state) gas phase heat conduction into the solid interior. From the viewpoint of the quasi-static gaseous flame, this is an additional heat loss that will affect the gas flame's flammability. In this work, a previously developed one-dimensional axisymmetric stagnation-point diffusion flame model is extended to account for the excessive heat loss into the solid beyond the steady burning state. A non-dimensional excessive conductive heat loss parameter Ψ is defined in the paper to account for the transient solid effect. For each Ψ, the gas phase is considered reaching a quasi-steady state. Extensive numerical computations have been performed to determine the extinction boundary and the characteristics of near-limit flames. The U-shaped flammability boundary is now a series of boundaries with Ψ as a parameter. Flame structures under the same environment condition but with different Ψ are compared. For a given oxygen ambient, the extinction boundary Ψ vs stretch rate is an inverted U-shape. The boundary consists of a blowoff branch and a radiative quenching branch. Extensive computed data including flame and pyrolysis temperatures, burning rate, flame standoff distance, species concentration, and reaction rate are presented along the extinction boundaries. Heat balance analysis is performed for both the gas phase and the gas-solid interface. Along the extinction boundary, the relative importance of conduction, convection, radiation, and chemical reaction are analyzed in detail.

Stagnation point diffusion flame↗

Mechanism of Macrosegregation Formation in Continuous Casting Slab: A Numerical Simulation Study

Solidified shell bulging is supposed to be the main reason for slab center segregation, while the influence of thermal shrinkage rarely has been considered. In this article, a thermal shrinkage model coupled with the multiphase solidification model is developed to investigate the effect of the thermal shrinkage, solidification shrinkage, grain sedimentation, and thermal flow on solute transport in the continuous casting slab. In this model, the initial equiaxed grains contract freely with the temperature decrease, while the coherent equiaxed grains and columnar phase move directionally toward the slab surface. The results demonstrate that the center positive segregation accompanied by negative segregation in the periphery zone is mainly caused by thermal shrinkage. During the solidification process, liquid phase first transports toward the slab surface to compensate for thermal shrinkage, which is similar to the case considering solidification shrinkage, and then it moves opposite to the slab center near the solidification end. It is attributed to the sharp decrease of center temperature and the intensive contract of solid phase, which cause the enriched liquid to be squeezed out. With the effect of grain sedimentation and thermal flow, the negative segregation at the external arc side (zone A1) and the positive segregation near the columnar-to-equiaxed transition at the inner arc side (position B1) come into being. Besides, it is found that the grain sedimentation and thermal flow only influence solute transport before equiaxed grains impinge with each other, while the solidification and thermal shrinkage still affect solute redistribution in the later stage.

36 MATERIALS SCIENCE↗

Numerical simulation of the Space Shuttle Launch Vehicle flowfield with real gas solid rocket plume effects

A numerical simulation of the external ascent flowfield of the Space Shuttle Launch Vehicle (SSLV) was carried out at the freestream Mach number 1.25, the angle of attack -5.1 deg, and the flight Reynolds number 3.25 x 10 exp 6/ft. The simulation is an extension of the solution by Kandula et al. (1991) and incorporates variable gamma effects with a high fidelity grid appropriate for a flight Reynolds number flow solution. Three-dimensional SSLV steady-state simulations with either perfect-gas or variable-gamma-gas Redesigned-Solid-Rocket-Motor (RSRM) plumes are computed on a 14-grid overlapping body-conforming grid system, and the influence of the RSRM exhaust plumes on the SSLV component pressure distributions and integrated loads is examined.

Slotnick, J. P.↗

Aerodynamic detuning for aeroelastic control of stability and forced response of supersonic rotors

An unsteady aerodynamic model is developed to analyze flutter and aerodynamically forced response of aerodynamically detuned supersonic axial flow rotors. Alternate blade aerodynamic detuning is considered, accomplished by alternating the circumferential spacing of adjacent blades as small solidity variations which do not have a dominant effect on the steady performance of a rotor. The unsteady aerodynamics are determined by developing an influence coefficient technique which is appropriate for both aerodynamically tuned and detuned rotor configurations. Torsion mode rotor stability and aerodynamically forced response are then analyzed with this unsteady aerodynamic model by combining it with a single-degree-of-freedom structural model. The effects of this detuning on the flutter and forced response characteristics of supersonic axial flow rotors is then demonstrated by considering baseline twelve bladed rotors.

Spara, Karen M.↗

Effect of including variable gas properties and entrained particles in the flow analysis of the ASRM nozzle

The objective of this study is to determine the effect of using constant gas properties for the analysis of the Advanced Solid Rocket Motor (ASRM) nozzle and to gain more understanding concerning those types of analysis which might require this additional complexity. Kinetics data for viscosity, thermal conductivity, specific heat, and specific heat ratio are extracted from the Solid Propellant Rocket Motor Performance Prediction Computer Program (SPP) and tabulated as a function of temperature. These tables are added to the Aerovisc computational fluid dynamics (CFD) code in place of the constant gas property values. The results of a CFD analysis of the ASRM 48 inch motor with constant gas properties will be compared with an analysis which uses variable gas properties. Mach number, surface pressure, and torque plots are presented. A full scale ASRM analysis using SPP with and without particle flow is presented.

Clayton, Curtis D.↗

Pore-scale modelling of multiphase reactive flow: application to mineral dissolution with production of CO 2

A micro-continuum approach is proposed to simulate the dissolution of solid minerals at the pore scale in the presence of multiple fluid phases. The approach employs an extended Darcy–Brinkman–Stokes formulation that accounts for the interfacial tension between the two immiscible fluid phases and the moving contact line at the mineral surface. The simulation framework is validated using an experimental microfluidic device that provides time-lapse images of the dissolution dynamics. The set-up involves a single-calcite crystal and the subsequent generation of CO 2 bubbles in the domain. The dissolution of the calcite crystal and the production of gas during the acidizing process are analysed. We then show that the production of CO 2 bubbles during the injection of acid in a carbonate formation may limit the overall dissolution rate and prevent the emergence of wormholes.

bio-inspired, mechanical behavior, carbon sequestr↗

Wall-Modeled Large Eddy Simulation Method for Unstructured-Grid Navier-Stokes Solvers

This paper reports on the implementation and assessment of a Wall-Modeled Large-Eddy Simulation (WMLES) methodology in an unstructured-grid, node-centered flow solver, FUN3D that is developed and supported at the NASA Langley Research Center. Finite-volume (FV) and finite-element (FE) discretization schemes considered in the study provide formal second-order spatial accuracy. Large-Eddy Simulations (LES) resolve large-scale turbulent-flow features and filter out small-scale effects using the Vreman subgrid-scale model. At solid-wall boundaries, a shear-stress model is employed to provide a proper boundary-flux closure. The nonlinear equations are integrated in time using either an optimized backward difference formula or an implicit multistage Runge-Kutta temporal scheme. The implicit equations at each time step are solved by strong nonlinear iteration schemes. WMLES demonstrations are shown for two high-lift configurations, namely, the McDonnell Douglas 30P30N multielement airfoil and a NASA High-Lift Common Research Model. Results show that the WMLES approaches implemented in the FV and FE discretization methods produce consistent solutions and are capable of capturing key aerodynamic characteristics and flow structures for high-lift configurations at a wide range of angles of attack including maximum-lift conditions. In the 30P30N example, correct trends in the variations of integrated aerodynamic forces and moments, surface pressure distributions, and boundary-layer profiles are captured as the Reynolds number is increased.

CFD; turbulence modeling; High-Lift flow simulatio↗

Heat transfer evaluation in a plasma core reactor

Numerical evaluations of heat transfer in a fissioning uranium plasma core reactor cavity, operating with seeded hydrogen propellant, was performed. A two-dimensional analysis is based on an assumed flow pattern and cavity wall heat exchange rate. Various iterative schemes were required by the nature of the radiative field and by the solid seed vaporization. Approximate formulations of the radiative heat flux are generally used, due to the complexity of the solution of a rigorously formulated problem. The present work analyzes the sensitivity of the results with respect to approximations of the radiative field, geometry, seed vaporization coefficients and flow pattern. The results present temperature, heat flux, density and optical depth distributions in the reactor cavity, acceptable simplifying assumptions, and iterative schemes. The present calculations, performed in cartesian and spherical coordinates, are applicable to any most general heat transfer problem.

Smith, D. E.↗

Skin Friction and Transition Location Measurement on Supersonic Transport Models

Flow visualization techniques were used to obtain both qualitative and quantitative skin friction and transition location data in wind tunnel tests performed on two supersonic transport models at Mach 2.40. Oil-film interferometry was useful for verifying boundary layer transition, but careful monitoring of model surface temperatures and systematic examination of the effects of tunnel start-up and shutdown transients will be required to achieve high levels of accuracy for skin friction measurements. A more common technique, use of a subliming solid to reveal transition location, was employed to correct drag measurements to a standard condition of all-turbulent flow on the wing. These corrected data were then analyzed to determine the additional correction required to account for the effect of the boundary layer trip devices.

Kennelly, Robert A., Jr.↗

Solid Surface Combustion Experiment: Thick Fuel Results

The results of experiments for spread over polymethylmethacrylate, PMMA, samples in the microgravity environment of the Space Shuttle are described. The results are coupled with modelling in an effort to describe the physics of the spread process for thick fuels in a quiescent, microgravity environment and uncover differences between thin and thick fuels. A quenching phenomenon not present for thin fuels is delineated, namely the fact that for thick fuels the possibility exists that, absent an opposing flow of sufficient strength to press the flame close enough to the fuel surface to allow the heated layer in the solid to develop, the heated layer fails to become 'fully developed.' The result is that the flame slows, which in turn causes an increase in the relative radiative loss from the flame, leading eventually to extinction. This potential inability of a thick fuel to develop a steady spread rate is not present for a thin fuel because the heated layer is the fuel thickness, which reaches a uniform temperature across the thickness relatively rapidly.

Altenkirch, Robert A.↗

Performance of a single-stage transonic compressor with a blade-tip solidity of 1.5 and comparison with 1.3 and 1.7 solidity stages

The overall and blade-element performance of a transonic compressor stage with a tip solidity of 1.5 is presented over the stable operating range at rotative speeds from 50 to 100 percent of design speed. State peak efficiency of 0.82 was obtained at a weight flow of 29.4 kg.sec (200.4 (kg/sec)/m2 of annulus area) and a pressure ratio of 1.71. Stall margin at design speed was 14 percent. A comparison of three stages in a solidity study showed that the performance of the 1.5 solidity stage and the 1.3 solidity stage were nearly identical but that the performance of the 1.7 solidity stage was significantly lower.

Osborn, W. M.↗