SEARCH · Engineering Papers

Results for “unsteady mixing”

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 235 records · Page 13

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Flux-difference split parabolized Navier-Stokes algorithm for non-equilibrium chemically reacting flows

A flux-difference split explicit finite-difference algorithm is presented for solving the parabolized form of the equations governing three-dimensional nonequilibrium chemically reacting flows. The algorithm is based on an explicit noniterative, upwind space-marching scheme developed by Korte, but differs in that the unsteady Riemann problem, rather than the steady Riemann problem, is solved. The algorithm allows either a second or an approximately third-order accurate upwind treatment of the convection terms by employing the unsteady approximate Riemann solver of Roe. The source terms of the species transport equations are treated in either an explicit or implicit manner, and the species diffusion terms are modeled with either a Fickian or a multicomponent model. A validation of the algorithm is performed by comparing computational results with the 2-D Mach 14, 15 degree compression-corner data of Holden. The three-dimensional capability of the algorithm is demonstrated by computing Mach 2.7 flow over a swept wedge scramjet fuel injector, and three-dimensional reacting flow capability is demonstrated by a computing a shock-jet interaction concept for mixing and combustion enhancement.

White, J. A.↗

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 NTRS1984

Recovery of burner acoustic source structure from far-field sound spectra

A method is presented that permits the thermal-acoustic efficiency spectrum in a long turbulent burner to be recovered from the corresponding far-field sound spectrum. An acoustic source/propagation model is used based on the perturbation solution of the equations describing the unsteady one-dimensional flow of an inviscid ideal gas with a distributed heat source. The technique is applied to a long cylindrical hydrogen-flame burner operating over power levels of 4.5-22.3 kW. The results show that the thermal-acoustic efficiency at a given frequency, defined as the fraction of the total burner power converted to acoustic energy at that frequency, is rather insensitive to burner power, having a maximum value on the order of 10 to the -4th at 150 Hz and rolling off steeply with increasing frequency. Evidence is presented that acoustic agitation of the flame at low frequencies enhances the mixing of the unburned fuel and air with the hot products of combustion. The paper establishes the potential of the technique as a useful tool for characterizing the acoustic source structure in any burner, such as a gas turbine combustor, for which a reasonable acoustic propagation model can be postulated.

Mahan, J. R.↗

NASA NTRS2001

Structure and Variability of Water Vapor in the Upper Troposphere and Lower Stratosphere

Upper-tropospheric humidity (UTH) has been synoptically mapped via an algorithm that rejects small-scale undersampled variance, which is intrinsic to asymptotic measurements of water vapor, cloud, and other convective properties. Mapped distributions of UTH have been used, jointly with high-resolution Global Cloud Imagery (GCI), to study how the upper troposphere is humidified. The time-mean distribution of UTH is spatially correlated to the time-mean distribution of cold cloud fraction (eta)(sub c) (T < than 230 K). Regions of large UTH coincide with regions of large eta(sub c), which mark deep convection. They also coincide with regions of reduced vertical stability, in which the vertical gradient of theta is weakened by convective mixing. Coldest cloud cover is attended convective overshoots above the local tropopause, which is simultaneously coldest and highest. Together, these features reflect the upper-troposphere being ventilated by convection, which mixes in moist air from lower levels. Histograms of UTH and eta(sub c) have been applied to construct the joint probability density function, which quantifies the relationship between these properties. The expected value of UTH in convective regions is strongly correlated to the expected value of eta(sub c). In ensembles of asymptotic samples, the correlation between epsilon[UTH] and epsilon[eta(sub c)] exceeds 0.80. As these expectations reflect the most likely values, the strong correlation between epsilon[UTH] and epsilon[eta(sub c)] indicates that the large-scale organization of UTH is strongly shaped by convective pumping of moisture from lower levels. The same relationship holds for unsteady fields - even though, instantaneously, those fields are comprised almost entirely of small-scale convective structure. The spatial autocorrelation of UTH, constructed at high resolution from overpass data along ascending and descending tracks of the orbit, is limited to only a couple of degrees in the horizontal. This mirrors the spatial autocorrelation of eta(sub c), which likewise operates coherently on short scales. The short correlation scale of UTH, which reflects the scale of individual convective systems, is comparable to the spacing of retrievals from MLS. These scales are undersampled in the asynoptic measurements. Despite their prevalence, the mapping algorithm described above successfully recovers synoptic behavior operating coherently on large scales. It reveals eastward migration of anomalous UTH from the Indian ocean to the central Pacific, in association with the modulation of convection by the Madden-Julian oscillation. Additional information is contained in the original extended abstract.

Salby, Murry L.↗

NASA NTRS2020

Implicit Large-Eddy Simulations of Hot and Cold Supersonic Jets in Loci-CHEM

Turbulence is a critical factor in most aero- and propulsion-related flows. For many years, the industry-standard approach to incorporating the effect of turbulence in computational fluid dynamics (CFD)calculations has been the Reynolds Averaged Navier Stokes (RANS) modeling approach. RANS implicitly time-averages turbulent motion, and models its effect on the mixing of species, momentum, and energy in a flow field. The approach is economical and often allows a CFD calculation to proceed to steady-state. However, it makes several simplifying assumptions, and the results will only be as good as the RANS turbulence model. Meanwhile, the last two decades have witnessed rapid and explosive growth in computational power, the development of CFD codes which scale well over thousands of processors, and the refinement of unstructured grid-generation tools which facilitate rapid surface and volume gridding of complex geometries. CFD engineering calculations of 108– 109 finite-volume cells have become routine. Thus, large eddy simulation (LES), which attempts to directly resolve the unsteady motion of the largest scales of turbulence, is increasingly of interest for many fluids engineering problems. Additionally, hybrid RANS-LES approaches such as Spalart’s Detached Eddy Simulation (DES), which attempt to marry the best strengths of both methods, are in active use and development. Because it seeks to accurately resolve the unsteady motion of the largest scales of turbulence, LES places certain requirements on the numerical method used to discretize and advance the Navier Stokes equations. In particular, LES requires high time accuracy and minimal intrinsic numerical dispersion and dissipation over a wide range of length scales.

C I Morris↗

DOE OSTI2023

Experiments and Simulations of Droplet Burning of Isobutanol Mixed with a Tier II Gasoline Certification Fuel and its Surrogate, Paper Number: SHTC2023-115267

This study presents an experimental and numerical results of combustion dynamics of a Tier II certification fuel/iso-butanol mixtures and its surrogate(S5)/iso-butanol mixtures. The experiments were carried out where gas transport was determined entirely by fuel evaporation to promote spherical symmetry during the burning process. The droplet initial diameters were in the range 0.55 mm to 0.62 mm. The combustion characteristics such as droplet diameters, flame standoff ratios and soot shell standoff ratios of indolene/iso-butanol mixtures are compared with S5/iso-butanol mixtures and the implications of iso-butanol addition to indolene and S5 on combustion are discussed. The presence of iso-butanol is experimentally found to reduce the formation of soot as the iso-butanol concentration was increased and the iso-butanol addition insignificantly influence the burning rate. Simulations were carried out with the detailed kinetic chemistry (324 species and 17608 reactions) from POLIMI. The detailed numerical model included unsteady gas and liquid transport, soot chemistry and radiative transport. The impact of iso-butanol on the burning rate, flame structure and soot formation are presented. Droplet burn rates were well predicted while flame diameters showed more variability depending on how the flame was defined in the numerical simulation (peak gas temperature and OH concentration). Simulations also show the influence of iso-butanol concentration on flame temperature, greenhouse gas emissions and soot formations.

Guo, Songtao↗

DOE OSTI2023

Combustion of petroleum-based transportation fuels and their blends with biofuels: a new approach for developing surrogates and understanding the effects of blending

Accurate and versatile reaction mechanisms are necessary to simulate combustion of transportation fuels blended with biofuels for improving fuel efficiency, reducing fuel consumption and mitigating the formation and emissions of combustion products harmful to the atmosphere. This project involved a combined effort of simulation and experimentation of liquid fuel burning using a combustion configuration amenable to numerical modeling with a level of detail not previously achieved for blends of biofuels with gasoline certification and surrogate fuels. The configuration of an isolated droplet burning under conditions where gas transport arises solely from fuel evaporation was selected as the platform for experiments and numerical modeling. The spherical symmetry and one-dimensional gas transport that results enabled simulating the droplet burning process with an precedented level of detail. Processes associated with unsteady gas and liquid transport, formation of particulate and gaseous products, radiation, multicomponent phase equilibrium at the droplet surface and moving boundary effects from droplet evaporation, were incorporated in a single numerical framework. The numerical model was based on the open source code OPENSmoke++ (OS) adapted to incorporate the effects noted above. In addition to simulation, an experimental design was developed for the isolated droplet configuration to acquire the data needed for validating the numerical model. Several broad accomplishments of the project were the following: demonstrating generally excellent agreement between measured and simulated combustion parameters for the fuel blends examined that included mixtures of heptane and isobutanol as a model system, and surrogates comprised of up to seven miscible components mixed with ethanol or isobutanol as representative biofuel additives; using a new approach to validate reaction mechanisms of biofuel blends which incorporated fuel evaporation into the process along with developing an experimental design for acquiring data to compare with simulations; and showing that a wealth of information could be obtained on the combustion physics of biofuels from experiments requiring volumes on the order of only nanoliters at a time thus opening the way to evaluating biofuels synthesized by new processes early in development.

02 PETROLEUM↗

NASA NTRS2016

Femtosecond Laser Tagging Characterization of a Sweeping Jet Actuator Operating in the Compressible Regime

A sweeping jet (SWJ) actuator operating over a range of nozzle pressure ratios (NPRs) was characterized with femtosecond laser electronic excitation tagging (FLEET), single hot-wire anemometry (HWA) and high-speed/phase-averaged schlieren. FLEET velocimetry was successfully demonstrated in a highly unsteady, oscillatory flow containing subsonic through supersonic velocities. Qualitative comparisons between FLEET and HWA (which measured mass flux since the flow was compressible) showed relatively good agreement in the external flow profiles. The spreading rate was found to vary from 0.5 to 1.2 depending on the pressure ratio. The precision of FLEET velocity measurements in the external flow field was poorer (is approximately equal to 25 m/s) than reported in a previous study due to the use of relatively low laser fluences, impacting the velocity fluctuation measurements. FLEET enabled velocity measurements inside the device and showed that choking likely occurred for NPR ≥ 2.0, and no internal shockwaves were present. Qualitative oxygen concentration measurements using FLEET were explored in an effort to gauge the jet's mixing with the ambient. The jet was shown to mix well within roughly four throat diameters and mix fully within roughly eight throat diameters. Schlieren provided visualization of the internal and external flow fields and showed that the qualitative structure of the internal flow does not vary with pressure ratio and the sweeping mechanism observed for incompressible NPRs also probably holds for compressible NPRs.

Peters, Christopher J.↗

NASA NTRS1995

An Extended Lagrangian Method

A unique formulation of describing fluid motion is presented. The method, referred to as 'extended Lagrangian method,' is interesting from both theoretical and numerical points of view. The formulation offers accuracy in numerical solution by avoiding numerical diffusion resulting from mixing of fluxes in the Eulerian description. The present method and the Arbitrary Lagrangian-Eulerian (ALE) method have a similarity in spirit-eliminating the cross-streamline numerical diffusion. For this purpose, we suggest a simple grid constraint condition and utilize an accurate discretization procedure. This grid constraint is only applied to the transverse cell face parallel to the local stream velocity, and hence our method for the steady state problems naturally reduces to the streamline-curvature method, without explicitly solving the steady stream-coordinate equations formulated a priori. Unlike the Lagrangian method proposed by Loh and Hui which is valid only for steady supersonic flows, the present method is general and capable of treating subsonic flows and supersonic flows as well as unsteady flows, simply by invoking in the same code an appropriate grid constraint suggested in this paper. The approach is found to be robust and stable. It automatically adapts to flow features without resorting to clustering, thereby maintaining rather uniform grid spacing throughout and large time step. Moreover, the method is shown to resolve multi-dimensional discontinuities with a high level of accuracy, similar to that found in one-dimensional problems.

Liou, Meng-Sing↗

NASA NTRS2003

Predicting Modes of the Unsteady Vorticity Field near the Trailing Edge of a Blade

Progress on predicting modes of the unsteady velocity/vorticity field of a turbulent boundary layer from Reynolds stress statistics is described. Prediction of these modes, that provide the source terms for trailing edge noise predictions in aircraft engine fans and other configurations, will allow for the first time detailed viscous flow effects to be included in such noise calculations. The key accomplishments of this work in FY02 are: (1) The development of a Matlab code for the prediction of modes in two- and three-dimensional boundary layers, previously applied to plane wakes; (2) Predictions with the code using a constant lengthscale formulation in a fully developed turbulence channel flow. Comparison of these boundary layer predictions with available DNS simulation results; and (3) Formulation of an improved model using a variable lengthscale proportional to mixing length. Turbulent channel flow predictions and comparison with DNS results. This work is being carried out in continuous communication and collaboration with the Glegg research group at Florida Atlantic University, which will be incorporating mode predictions into engine noise calculations.

Devenport, William J.↗

DOE OSTI2026

Modeling and Experimental Validation of a Direct-Contact Counter-Flow Fluidized Bed Heat Exchanger for Thermal Energy Storage (TES) Applications

Particle-based thermal energy storage (TES) systems are an emerging energy storage technology. The technological advances have reduced costs, making TES more competitive and reliable in the marketplace but an efficient and reliable operation is heavily dependent on coherent heat transfer between air to particles or vice versa. The particle-based TES technologies provide an intermediate system that can store energy for short (0-10 h), long (10-200 h) and seasonal (> 200 h) timescales. The TES systems store energy by converting electricity to thermal energy; electricity can be directly sourced intermittent generation technologies and/or the grid, helping manage peak loads and other mismatches in supply and demand. The overall efficiency of the TES system depends on the performance of system components (particle storage silos and particle transfer mechanism etc.). The particle heat exchanger is one of the key system components that affects the system efficiency. The pressurized fluidized bed heat exchanger (PFB HX) performance is challenging to predict due to the chaotic behavior of particle and fluid interaction. This research presents a computational study of a novel direct-contact, counter-flow and air-to-particles PFB HX, that contributes in advancing the particle-based long-duration TES technologies. For the current analysis an unsteady Eulerian-Eulerian CFD model was developed and validated against experiments performed at the National Laboratory of the Rockies for two particle sizes (600 ..mu..m and 825 ..mu..m ). Following validation, parametric simulations were conducted to evaluate the effects of interphase drag models (Syamlal-O'Brien and Gidaspow), particle size, bed height and the influence of a frictional-viscosity term on hydrodynamics and heat transfer between the air & particles. The key findings from the analysis are: (1) for the studied operating window Syamlal-O'Brien provides superior agreement with measured gas temperatures (errors generally < 10%) while Gidaspow shows large deviations for the coarse particle case; (2) model predictions are most sensitive in the lower 0.2 m above the air distributor where bubble initiation and local mixing dominate interphase heat transfer; (3) representation of the distributor (number of inlet ports) materially affects predicted local mixing and temperature stratification; and (4) the Eulerian-Eulerian framework reproduces bulk thermal trends but shows regime dependent limitations for coarse particles, motivating mesoscale informed closures for scale-up analysis for future studies. These results provide validated guidance for drag selection and distributor design in particle-based thermal energy storage applications. Collectively, the validated model and parametric results quantify key drivers of PHB-HX performance and provide practical guidance for design and optimization. The results provide confidence in the model predictability and provide a step forward to improve on heat exchange performance. The demonstrated performance and modeling approach support the deployment and further development of this novel PHB-HX concept for robust, particle-based long-duration thermal energy storage systems.

25 ENERGY STORAGE↗

NASA NTRS1986

Transition and mixing in axisymmetric jets and vortex rings

A class of impulsively started, axisymmetric, laminar jets produced by a time dependent joint source of momentum are considered. These jets are different flows, each initially at rest in an unbounded fluid. The study is conducted at three levels of detail. First, a generalized set of analytic creeping flow solutions are derived with a method of flow classification. Second, from this set, three specific creeping flow solutions are studied in detail: the vortex ring, the round jet, and the ramp jet. This study involves derivation of vorticity, stream function, entrainment diagrams, and evolution of time lines through computer animation. From entrainment diagrams, critical points are derived and analyzed. The flow geometry is dictated by the properties and location of critical points which undergo bifurcation and topological transformation (a form of transition) with changing Reynolds number. Transition Reynolds numbers were calculated. A state space trajectory was derived describing the topological behavior of these critical points. This state space derivation yielded three states of motion which are universal for all axisymmetric jets. Third, the axisymmetric round jet is solved numerically using the unsteady laminar Navier Stokes equations. These equations were shown to be self similar for the round jet. Numerical calculations were performed up to a Reynolds number of 30 for a 60x60 point mesh. Animations generated from numerical solution showed each of the three states of motion for the round jet, including the Re = 30 case.

Allen, G. A., Jr.↗

DOE OSTI2024

A numerical method for simulating variable density flows in membrane desalination systems

Here, we present a novel method for simulating unsteady, variable density, fluid flows in membrane desalination systems. By assuming the density varies only with concentration and temperature, the scheme decouples the solution of the governing equations into two sequential blocks. The first solves the governing equations for the temperature and concentration fields, which are used to compute all thermophysical properties. The second block solves the conservation of mass and momentum equations for the velocity and pressure. We show that this is computationally more efficient than schemes that iterate over the full coupled equations in one block. We verify that the method achieves second-order spatial-temporal accuracy, and we use the method to investigate buoyancy-driven convection in a desalination process called vacuum membrane distillation. Specifically, we show that with gravity properly oriented, variations in temperature and concentration can trigger a double-diffusive instability that enhances mixing and improves water recovery. We also show that the instability can be strengthened by providing external heating.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2001

Multiple Mode Actuation of a Turbulent Jet

The effects of multiple mode periodic excitation on the evolution of a circular turbulent jet were studied experimentally. A short, wide-angle diffuser was attached to the jet exit. Streamwise and cross-stream excitations were introduced at the junction between the jet exit and the diffuser inlet on opposing sides of the jet. The introduction of high amplitude, periodic excitation in the streamwise direction enhances the mixing and promotes attachment of the jet shear-layer to the diffuser wall. Cross-stream excitation applied over a fraction of the jet circumference can deflect the jet away from the excitation slot. The two modes of excitation were combined using identical frequencies and varying the relative phase between the two actuators in search of an optimal response. It is shown that, for low and moderate periodic momentum input levels, the jet deflection angles depend strongly on the relative phase between the two actuators. Optimum performance is achieved when the phase difference is pi +/- pi/6. The lower effectiveness of the equal phase excitation is attributed to the generation of an azimuthally symmetric mode that does not produce the required non-axisymmetric vectoring. For high excitation levels, identical phase becomes more effective, while phase sensitivity decreases. An important finding was that with proper phase tuning, two unsteady actuators can be combined to obtain a non-linear response greater than the superposition of the individual effects.

Pack, LaTunia G.↗

NASA NTRS2001

Multiple Mode Actuation of a Turbulent Jet

Pack, LaTunia G.↗

NASA NTRS1993

A time-accurate finite volume method valid at all flow velocities

A finite volume method to solve the Navier-Stokes equations at all flow velocities (e.g., incompressible, subsonic, transonic, supersonic and hypersonic flows) is presented. The numerical method is based on a finite volume method that incorporates a pressure-staggered mesh and an incremental pressure equation for the conservation of mass. Comparison of three generally accepted time-advancing schemes, i.e., Simplified Marker-and-Cell (SMAC), Pressure-Implicit-Splitting of Operators (PISO), and Iterative-Time-Advancing (ITA) scheme, are made by solving a lid-driven polar cavity flow and self-sustained oscillatory flows over circular and square cylinders. Calculated results show that the ITA is the most stable numerically and yields the most accurate results. The SMAC is the most efficient computationally and is as stable as the ITA. It is shown that the PISO is the most weakly convergent and it exhibits an undesirable strong dependence on the time-step size. The degenerated numerical results obtained using the PISO are attributed to its second corrector step that cause the numerical results to deviate further from a divergence free velocity field. The accurate numerical results obtained using the ITA is attributed to its capability to resolve the nonlinearity of the Navier-Stokes equations. The present numerical method that incorporates the ITA is used to solve an unsteady transitional flow over an oscillating airfoil and a chemically reacting flow of hydrogen in a vitiated supersonic airstream. The turbulence fields in these flow cases are described using multiple-time-scale turbulence equations. For the unsteady transitional over an oscillating airfoil, the fluid flow is described using ensemble-averaged Navier-Stokes equations defined on the Lagrangian-Eulerian coordinates. It is shown that the numerical method successfully predicts the large dynamic stall vortex (DSV) and the trailing edge vortex (TEV) that are periodically generated by the oscillating airfoil. The calculated streaklines are in very good comparison with the experimentally obtained smoke picture. The calculated turbulent viscosity contours show that the transition from laminar to turbulent state and the relaminarization occur widely in space as well as in time. The ensemble-averaged velocity profiles are also in good agreement with the measured data and the good comparison indicates that the numerical method as well as the multipletime-scale turbulence equations successfully predict the unsteady transitional turbulence field. The chemical reactions for the hydrogen in the vitiated supersonic airstream are described using 9 chemical species and 48 reaction-steps. Consider that a fast chemistry can not be used to describe the fine details (such as the instability) of chemically reacting flows while a reduced chemical kinetics can not be used confidently due to the uncertainty contained in the reaction mechanisms. However, the use of a detailed finite rate chemistry may make it difficult to obtain a fully converged solution due to the coupling between the large number of flow, turbulence, and chemical equations. The numerical results obtained in the present study are in good agreement with the measured data. The good comparison is attributed to the numerical method that can yield strongly converged results for the reacting flow and to the use of the multiple-time-scale turbulence equations that can accurately describe the mixing of the fuel and the oxidant.

Kim, S.-W.↗

NASA NTRS1992

Surface tension and buoyancy-driven flow in a non-isothermal liquid bridge

The Navier-Stokes-Boussinesq equations governing the transport of momentum, mass and heat in a nonisothermal liquid bridge with a temperature-dependent surface tension are solved using a vorticity-stream-function formulation together with a nonorthogonal coordinate transformation. The equations are discretized using a pseudo-unsteady semi-implicit finite difference scheme and are solved by the ADI method. A Picard-type iteration is adopted which consists of inner and outer iterative processes. The outer iteration is used to update the shape of the free surface. Two schemes have been used for the outer iteration; both use the force balance normal to the free surface as the distinguished boundary condition. The first scheme involves successive approximation by the direct solution of the distinguished boundary condition. The second scheme uses the artificial force imbalance between the fluid pressure, viscous and capillary forces at the free surface which arises when the boundary condition for force balance normal to the surface is not satisfied. This artificial imbalance is then used to change the surface shape until the distinguished boundary condition is satisfied. These schemes have been used to examine a variety of model liquid bridge situations including purely thermocapillary-driven flow situations and mixed thermocapillary- and bouyancy-driven flow.

Zhang, Yiqiang↗

NASA NTRS1991

Three-Dimensional Analysis and Modeling of a Wankel Engine

A new computer code, AGNI-3D, has been developed for the modeling of combustion, spray, and flow properties in a stratified-charge rotary engine (SCRE). The mathematical and numerical details of the new code are described by the first author in a separate NASA publication. The solution procedure is based on an Eulerian-Lagrangian approach where the unsteady, three-dimensional Navier-Stokes equations for a perfect gas-mixture with variable properties are solved in generalized, Eulerian coordinates on a moving grid by making use of an implicit finite-volume, Steger-Warming flux vector splitting scheme. The liquid-phase equations are solved in Lagrangian coordinates. The engine configuration studied was similar to existing rotary engine flow-visualization and hot-firing test rigs. The results of limited test cases indicate a good degree of qualitative agreement between the predicted and measured pressures. It is conjectured that the impulsive nature of the torque generated by the observed pressure nonuniformity may be one of the mechanisms responsible for the excessive wear of the timing gears observed during the early stages of the rotary combustion engine (RCE) development. It was identified that the turbulence intensities near top-dead-center were dominated by the compression process and only slightly influenced by the intake and exhaust processes. Slow mixing resulting from small turbulence intensities within the rotor pocket and also from a lack of formation of any significant recirculation regions within the rotor pocket were identified as the major factors leading to incomplete combustion. Detailed flowfield results during exhaust and intake, fuel injection, fuel vaporization, combustion, mixing and expansion processes are also presented. The numerical procedure is very efficient as it takes 7 to 10 CPU hours on a CRAY Y-MP for one entire engine cycle when the computations are performed over a 31 x16 x 20 grid.

Raju, M. S.↗