Engineering topics

Shang, H. M.

Publications and source records attributed to Shang, H. M..

At least 19 records

NASA NTRS2001

Computation of Coupled Thermal-Fluid Problems in Distributed Memory Environment

The thermal-fluid coupling problems are very important to aerospace and engineering applications. Instead of analyzing heat transfer and fluid flow separately, this study merged two well-accepted engineering solution methods, SINDA for thermal analysis and FDNS for fluid flow simulation, into a unified multi-disciplinary thermal fluid prediction method. A fully conservative patched grid interface algorithm for arbitrary two-dimensional and three-dimensional geometry has been developed. The state-of-the-art parallel computing concept was used to couple SINDA and FDNS for the communication of boundary conditions through PVM (Parallel Virtual Machine) libraries. Therefore, the thermal analysis performed by SINDA and the fluid flow calculated by FDNS are fully coupled to obtain steady state or transient solutions. The natural convection between two thick-walled eccentric tubes was calculated and the predicted results match the experiment data perfectly. A 3-D rocket engine model and a real 3-D SSME geometry were used to test the current model, and the reasonable temperature field was obtained.

Wei, H.↗

NASA NTRS1996

The study of flow pattern and phase-change problem in die casting process

The flow pattern and solidification phenomena in die casting process have been investigated in the first phase study. The flow pattern in filling process is predicted by using a VOF (volume of fluid) method. A good agreement with experimental observation is obtained for filling the water into a die cavity with different gate geometry and with an obstacle in the cavity. An enthalpy method has been applied to solve the solidification problem. By treating the latent heat implicitly into the enthalpy instead of explicitly into the source term, the CPU time can be reduced at least 20 times. The effect of material properties on solidification fronts is tested. It concludes that the dependence of properties on temperature is significant. The influence of the natural convection over the diffusion has also been studied. The result shows that the liquid metal solidification phenomena is diffusion dominant, and the natural convection can affect the shape of the interface. In the second phase study, the filling and solidification processes will be considered simultaneously.

Wang, T. S.↗

NASA NTRS1996

Numerical Modeling of Spray Combustion with an Unstructured-Grid Method

The present unstructured-grid method follows strictly the basic finite volume forms of the conservation laws of the governing equations for the entire flow domain. High-order spatially accurate formulation has been employed for the numerical solutions of the Navier-Stokes equations. A two-equation k-epsilon turbulence model is also incorporated in the unstructured-grid solver. The convergence of the resulted linear algebraic equation is accelerated with preconditioned Conjugate Gradient method. A statistical spray combustion model has been incorporated into the present unstructured-grid solver. In this model, spray is represented by discrete particles, rather than by continuous distributions. A finite number of computational particles are used to predict a sample of total population of particles. Particle trajectories are integrated using their momentum and motion equations and particles exchange mass, momentum and energy with the gas within the computational cell in which they are located. The interaction calculations are performed simultaneously and eliminate global iteration for the two-phase momentum exchange. A transient spray flame in a high pressure combustion chamber is predicted and then the solution of liquid-fuel combusting flow with a rotating cup atomizer is presented and compared with the experimental data. The major conclusion of this investigation is that the unstructured-grid method can be employed to study very complicated flow fields of turbulent spray combustion. Grid adaptation can be easily achieved in any flow domain such as droplet evaporation and combustion zone. Future applications of the present model can be found in the full three-dimensional study of flow fields of gas turbine and liquid propulsion engine combustion chambers with multi-injectors.

Shang, H. M.↗

NASA NTRS1993

A numerical procedure for analysis of finite rate reacting flows

Combustion processes in rocket propulsion systems are characterized by the existence of multiple, vastly differing time and length scales, as well as flow-speeds at wide variation of Mach numbers. The chemical kinetics processes in the highly active reaction zone are characterized by much smaller scales compared to fluid convective and diffusive time scales. An operator splitting procedure for transient finite rate chemistry problems has been developed using a pressure based method, which can be applied to all speed flows without difficulties. The splitting of chemical kinetics terms formed the fluid-mechanical terms of the species equation ameliorated the difficulties associated with the disparate time scales and stiffness in the set of equations which describes highly exothermic combustion. A combined efficient ordinary differential equations (ODE) solver was used to integrate the effective chemical source terms over the residence time at each grid cell. One and two dimensional reacting flow situations were carried out to demonstrate and verify the current procedure. Different chemical kinetics with different degrees of nonlinearity have also been incorporated to test the robustness and generality of the proposed method.

Shang, H. M.↗

NASA NTRS1993

Development of an algebraic stress/two-layer model for calculating thrust chamber flow fields

Following the consensus of a workshop in Turbulence Modeling for Liquid Rocket Thrust Chambers, the current effort was undertaken to study the effects of second-order closure on the predictions of thermochemical flow fields. To reduce the instability and computational intensity of the full second-order Reynolds Stress Model, an Algebraic Stress Model (ASM) coupled with a two-layer near wall treatment was developed. Various test problems, including the compressible boundary layer with adiabatic and cooled walls, recirculating flows, swirling flows and the entire SSME nozzle flow were studied to assess the performance of the current model. Detailed calculations for the SSME exit wall flow around the nozzle manifold were executed. As to the overall flow predictions, the ASM removes another assumption for appropriate comparison with experimental data, to account for the non-isotropic turbulence effects.

Chen, C. P.↗

NASA NTRS1993

Liquid propellant rocket engine combustion simulation with a time-accurate CFD method

Time-accurate computational fluid dynamics (CFD) algorithms are among the basic requirements as an engineering or research tool for realistic simulations of transient combustion phenomena, such as combustion instability, transient start-up, etc., inside the rocket engine combustion chamber. A time-accurate pressure based method is employed in the FDNS code for combustion model development. This is in connection with other program development activities such as spray combustion model development and efficient finite-rate chemistry solution method implementation. In the present study, a second-order time-accurate time-marching scheme is employed. For better spatial resolutions near discontinuities (e.g., shocks, contact discontinuities), a 3rd-order accurate TVD scheme for modeling the convection terms is implemented in the FDNS code. Necessary modification to the predictor/multi-corrector solution algorithm in order to maintain time-accurate wave propagation is also investigated. Benchmark 1-D and multidimensional test cases, which include the classical shock tube wave propagation problems, resonant pipe test case, unsteady flow development of a blast tube test case, and H2/O2 rocket engine chamber combustion start-up transient simulation, etc., are investigated to validate and demonstrate the accuracy and robustness of the present numerical scheme and solution algorithm.

Chen, Y. S.↗

NASA NTRS1993

Advanced Multi-Phase Flow CFD Model Development for Solid Rocket Motor Flowfield Analysis

It is known that the simulations of solid rocket motor internal flow field with AL-based propellants require complex multi-phase turbulent flow model. The objective of this study is to develop an advanced particulate multi-phase flow model which includes the effects of particle dynamics, chemical reaction and hot gas flow turbulence. The inclusion of particle agglomeration, particle/gas reaction and mass transfer, particle collision, coalescence and breakup mechanisms in modeling the particle dynamics will allow the proposed model to realistically simulate the flowfield inside a solid rocket motor. The Finite Difference Navier-Stokes numerical code FDNS is used to simulate the steady-state multi-phase particulate flow field for a 3-zone 2-D axisymmetric ASRM model and a 6-zone 3-D ASRM model at launch conditions. The 2-D model includes aft-end cavity and submerged nozzle. The 3-D model represents the whole ASRM geometry, including additional grain port area in the gas cavity and two inhibitors. FDNS is a pressure based finite difference Navier-Stokes flow solver with time-accurate adaptive second-order upwind schemes, standard and extended k-epsilon models with compressibility corrections, multi zone body-fitted formulations, and turbulence particle interaction model. Eulerian/Lagrangian multi-phase solution method is applied for multi-zone mesh. To simulate the chemical reaction, penalty function corrected efficient finite-rate chemistry integration method is used in FDNS. For the AL particle combustion rate, the Hermsen correlation is employed. To simulate the turbulent dispersion of particles, the Gaussian probability distribution with standard deviation equal to (2k/3)(exp 1/2) is used for the random turbulent velocity components. The computational results reveal that the flow field near the juncture of aft-end cavity and the submerged nozzle is very complex. The effects of the turbulent particles affect the flow field significantly and provide better prediction of the ASRM performance. The multi-phase flow analysis using the FDNS code in the present research can be used as a design tool for solid rocket motor applications.

Liaw, Paul↗

NASA NTRS1993

Spray combustion model improvement study, 1

This study involves the development of numerical and physical modeling in spray combustion. These modeling efforts are mainly motivated to improve the physical submodels of turbulence, combustion, atomization, dense spray effects, and group vaporization. The present mathematical formulation can be easily implemented in any time-marching multiple pressure correction methodologies such as MAST code. A sequence of validation cases includes the nonevaporating, evaporating and_burnin dense_sprays.

Chen, C. P.↗

NASA NTRS1993

Turbulence modelling of flow fields in thrust chambers

Following the consensus of a workshop in Turbulence Modelling for Liquid Rocket Thrust Chambers, the current effort was undertaken to study the effects of second-order closure on the predictions of thermochemical flow fields. To reduce the instability and computational intensity of the full second-order Reynolds Stress Model, an Algebraic Stress Model (ASM) coupled with a two-layer near wall treatment was developed. Various test problems, including the compressible boundary layer with adiabatic and cooled walls, recirculating flows, swirling flows, and the entire SSME nozzle flow were studied to assess the performance of the current model. Detailed calculations for the SSME exit wall flow around the nozzle manifold were executed. As to the overall flow predictions, the ASM removes another assumption for appropriate comparison with experimental data to account for the non-isotropic turbulence effects.

Chen, C. P.↗

NASA NTRS1992

An efficient pressure-velocity procedure for gas-droplet two-phase flow calculations

This paper presents a noniterative numerical technique for computing time-dependent gas-droplet flows. The method is a fully interacting combination of Eulerian fluid and Lagrangian particle calculations. The interaction calculations between the two phases are formulated on a pressure-velocity-coupling procedure based on the operator-splitting technique. This procedure eliminates the global iterations required in the conventional particle-source-in-cell procedure. Turbulent dispersion calculations are treated by a stochastic procedure. Numerical calculations and comparisons with available experimental data as well as efficiency assessments are given for some sprays typical of spray combustion applications.

Chen, C. P.↗

NASA NTRS1992

Numerical modeling for dilute and dense sprays

We have successfully implemented a numerical model for spray-combustion calculations. In this model, the governing gas-phase equations in Eulerian coordinate are solved by a time-marching multiple pressure correction procedure based on the operator-splitting technique. The droplet-phase equations in Lagrangian coordinate are solved by a stochastic discrete particle technique. In order to simplify the calculation procedure for the circulating droplets, the effective conductivity model is utilized. The k-epsilon models are utilized to characterize the time and length scales of the gas phase in conjunction with turbulent modulation by droplets and droplet dispersion by turbulence. This method entails random sampling of instantaneous gas flow properties and the stochastic process requires a large number of computational parcels to produce the satisfactory dispersion distributions even for rather dilute sprays. Two major improvements in spray combustion modelings were made. Firstly, we have developed a probability density function approach in multidimensional space to represent a specific computational particle. Secondly, we incorporate the Taylor Analogy Breakup (TAB) model for handling the dense spray effects. This breakup model is based on the reasonable assumption that atomization and drop breakup are indistinguishable processes within a dense spray near the nozzle exit. Accordingly, atomization is prescribed by injecting drops which have a characteristic size equal to the nozzle exit diameter. Example problems include the nearly homogeneous and inhomogeneous turbulent particle dispersion, and the non-evaporating, evaporating, and burning dense sprays. Comparison with experimental data will be discussed in detail.

Chen, C. P.↗

NASA NTRS1992

Studies on fuel spray characteristics in high-pressure environment

The present study deals with several issues involving the improvement of physical submodels and the computational efficiency in modeling dense fuel sprays. To improve the computational efficiency, a parcel PDF approach is implemented which can account for turbulent dispersion within each computational parcel. The advantage of a parcel PDF tracking method is to reduce the number of computational parcels representing the spray dynamics as well as to obtain grid-independent solutions for two-phase flows. To account for the dense spray effects, an existing drop collision and coalescence model, two breakup models, and a Reitz's wave instability model were used. These models were incorporated into a state-of-the-art multiphase all-speed transient flow solution procedure. Comparative performance for two breakup models as well as the turbulence modulation effects are also studied. Validation cases include the nonevaporating and evaporating solid-cone dense sprays. The predictions show a reasonably good agreement with available experimental results in terms of spray penetration, drop sizes, gas and drop mean velocities, and gas and drop rms velocities. The numerical results indicate that the present parcel PDF model has the capability of accurately representing drop dispersion in dense sprays with manageable number of computational parcels.

Shang, H. M.↗

NASA NTRS1992

Computational fluid dynamics combustion analysis evaluation

This study involves the development of numerical modelling in spray combustion. These modelling efforts are mainly motivated to improve the computational efficiency in the stochastic particle tracking method as well as to incorporate the physical submodels of turbulence, combustion, vaporization, and dense spray effects. The present mathematical formulation and numerical methodologies can be casted in any time-marching pressure correction methodologies (PCM) such as FDNS code and MAST code. A sequence of validation cases involving steady burning sprays and transient evaporating sprays will be included.

Kim, Y. M.↗

NASA NTRS1992

Numerical studies of dilute and dense spray characteristics

Several issues involving the improvement of physical submodels and the computational efficiency in modeling dilute and dense spray combustion are discussed. First, the implementations of a dispersion width approach accounting for turbulent dispersion within each computational parcel is discussed. This is essentially a statistical transport model and the testings of this model confirm the capability of accurately representing dispersion in nearly-homogeneous and inhomogeneous turbulent flows with improved efficiency over the delta function stochastic separated flow model. To account for the dense spray effects, an existing drop collision and coalescence model and a Taylor analogy breakup (TAB) model were employed. These models were incorporated into a state-of-the-art multiphase all-speed transient flow solution procedure. Several examples including nonevaporating, evaporating, and burning dense spray cases were studied. The numerical results show reasonably good comparisons with available experimental data in terms of spray penetration, drop sizes, and overall configuration of a spray flame.

Kim, Y. M.↗

NASA NTRS1992

Characterization and modeling for dilute and dense sprays

The present numerical modeling of dilute and dense spray combustion solves the governing gas-phase equations in Eulerian coordinates by means of a time-marching multiple pressure correction procedure that is based on the operator-splitting technique. A sequence of validating cases involving (1) nonevaporating, (2) evaporating, (3) burning, (4) dilute, and (5) dense spray cases is undertaken; it is found that in the combusting dilute spray case, the present procedure correctly predicts the general features of flows and yields qualitative agreement with experimental data. In the other cases, reasonably good agreement with experimental results is obtained.

Chen, C. P.↗

NASA NTRS1991

A computer code for multiphase all-speed transient flows in complex geometries. MAST version 1.0

The operation of the MAST code, which computes transient solutions to the multiphase flow equations applicable to all-speed flows, is described. Two-phase flows are formulated based on the Eulerian-Lagrange scheme in which the continuous phase is described by the Navier-Stokes equation (or Reynolds equations for turbulent flows). Dispersed phase is formulated by a Lagrangian tracking scheme. The numerical solution algorithms utilized for fluid flows is a newly developed pressure-implicit algorithm based on the operator-splitting technique in generalized nonorthogonal coordinates. This operator split allows separate operation on each of the variable fields to handle pressure-velocity coupling. The obtained pressure correction equation has the hyperbolic nature and is effective for Mach numbers ranging from the incompressible limit to supersonic flow regimes. The present code adopts a nonstaggered grid arrangement; thus, the velocity components and other dependent variables are collocated at the same grid. A sequence of benchmark-quality problems, including incompressible, subsonic, transonic, supersonic, gas-droplet two-phase flows, as well as spray-combustion problems, were performed to demonstrate the robustness and accuracy of the present code.

Chen, C. P.↗

NASA NTRS1991

Non-isotropic turbulence effects on spray combustion

A numerical model for the prediction of local properties of statistically stationary spray-combusting flows is evaluated by comparison with experimental data. To appraise the relative performance of turbulence models, computations were carried out by the k-epsilon model and the algebraic stress model. The present numerical results show the qualitative agreement with experimental data. In terms of overall local flow properties, the algebraic stress model improves a degree of conformity to the experimental data due to its ability to introduce the nonisotropic turbulence effects. Two swirl numbers are considered to investigate the influence of swirl on the droplet evaporation and trajectories, and the effects of droplet/turbulence interactions in flow properties. It is found that the large swirl produces a higher evaporationn rate, and more intensive turbulent mixing and burning. The discrepancies observed in the results are attributed mainly to uncertainties in the initial spray size and velocity distributions, the droplet/wall impingement interaction, the combustion model with the fast chemistry and the turbulence models dealing with the strong streamline curvature and complex interactions between the dispersed droplets and the continuous gas-phase flows.

Kim, Y. M.↗

NASA NTRS1991

A novel gas-droplet numerical method for spray combustion

This paper presents a non-iterative numerical technique for computing time-dependent gas-droplet flows. The method is a fully-interacting combination of Eulerian fluid and Lagrangian particle calculation. The interaction calculations between the two phases are formulated on a pressure-velocity coupling procedure based on the operator-splitting technique. This procedure eliminates the global iterations required in the conventional particle-source-in-cell (PSIC) procedure. Turbulent dispersion calculations are treated by a stochastic procedure. Numerical calculations and comparisons with available experimental data, as well as efficiency assessments are given for some sprays typical of spray combustion applications.

Chen, C. P.↗