Three-dimensional Euler time accurate simulations of fan rotor-stator interactions

A numerical method useful to describe unsteady 3-D flow fields within turbomachinery stages is presented. The method solves the compressible, time dependent, Euler conservation equations with a finite volume, flux splitting, total variation diminishing, approximately factored, implicit scheme. Multiblock composite gridding is used to partition the flow field into a specified arrangement of blocks with static and dynamic interfaces. The code is optimized to take full advantage of the processing power and speed of the Cray Y/MP supercomputer. The method is applied to the computation of the flow field within a single stage, axial flow fan, thus reproducing the unsteady 3-D rotor-stator interaction

Eulerian method for multiphase interactions of soft solid bodies in fluids

We introduce an Eulerian approach for problems involving one or more soft solids immersed in a fluid, which permits mechanical interactions between all phases. The reference map variable is exploited to simulate finite-deformation constitutive relations in the solid(s) on the same fixed grid as the fluid phase, which greatly simplifies the coupling between phases. Our coupling procedure, a key contribution in the current work, is shown to be computationally faster and more stable than an earlier approach, and admits the ability to simulate both fluid--solid and solid--solid interaction between submerged bodies. The interface treatment is demonstrated with multiple examples involving a weakly compressible Navier--Stokes fluid interacting with a neo-Hookean solid, and we verify the method's convergence. The solid contact method, which exploits distance-measures already existing on the grid, is demonstrated with two examples. A new, general routine for cross-interface extrapolation is introduced and used as part of the new interfacial treatment

Heat Transfer Mechanism In Particle-Laden Turbulent Shearless Flows

Particle-laden turbulent flows are one of the complex flow regimes involved in a wide range of environmental, industrial, biomedical and aeronautical applications. Recently the interest has included also the interaction between scalars and particles, and the complex scenario which arises from the interaction of particle finite inertia, temperature transport, and momentum and heat feedback of particles on the flow leads to a multi-scale and multi-physics phenomenon which is not yet fully understood. The present work aims to investigate the fluid-particle thermal interaction in turbulent mixing under one-way and two-way coupling regimes. A recent novel numerical framework has been used to investigate the impact of suspended sub-Kolmogorov inertial particles on heat transfer within the mixing layer which develops at the interface of two regions with different temperature in an isotropic turbulent flow. Temperature has been considered a passive scalar, advected by the solenoidal velocity field, and subject to the particle thermal feedback in the two-way regime. A self-similar stage always develops where all single-point statistics of the carrier fluid and the suspended particles collapse when properly re-scaled. We quantify the effect of particle inertial, parametrized through the Stokes and thermal Stokes numbers, on the heat transfer through the Nusselt number, defined as the ratio of the heat transfer to the thermal diffusion. A scale analysis will be presented. We show how the modulation of fluid temperature gradients due to the statistical alignments of the particle velocity and the local carrier flow temperature gradient field, impacts the overall heat transfer in the two-way coupling regime

Improvement of Labyrinth Seal Leakage Rates Using Additive Manufacturing

The growing popularity of additive manufacturing in commercial applications has al- lowed for new ideas and new ways of thinking when designing components. Further optimization at the component level is possible, though powder metallurgy is in its infancy. This study explores the possibility of using additive manufacturing to develop better labyrinth seals in turbomachinery. Labyrinth seals have a torturous fluid path with high losses, thus limiting the amount of fluid leakage. These types of seals can be non-rotating, allowing them to better take advantage of the additive manufacturing process due to the absence of rotating stresses. Labyrinth seal performance is defined by its ability to limit leakage through a seal. Investigations on the ability to reduce this leakage using the inherent roughness from the additive manufacturing process and the addition of complex geometry only capable of being produced by additive manufacturing are explored. Incompressible and compressible fluid models are utilized in the study. Perfectly smooth seals with tooth counts of four, six, and eight are first simulated using ANSYS FLUENT and compared to theoretical models to determine accuracy. Roughness is then applied to the seals and leakage decreases of 0.5% to 1.5% are experienced for the incompressible model. Decreases of 1.0% to 3.5% are experienced for the compressible model. Flow visualization and line analysis are conducted for all seals tested to understand how fluid flow is behaving within the clearance region of the seal and seal chambers. Several additive manufacturing geometries are simulated against a control seal to determine geometries with the largest effect on leakage rates. These geometries are then adapted to a six tooth seal and simulated with roughness. This additively manufactured seal is then compared to the smooth and rough six tooth seal for both incompressible and compressible fluids. Leakage was decreased by 5% to 8% for the incompressible model and 5% to 7% reductions for the compressible model when compared to the smooth seal. Flow visualization and line analysis were also conducted for the additively manufactured six tooth seal. A basic outline for an experiment and test stand were developed for future work

Large Eddy Simulation of the Pulsating, Non-Reacting Flow in Combustion Chambers

In the present work a numerical investigation of the pulsating, non-reacting flow in a combustion chamber, and a coupled system of burner plenum and combustion chamber is presented. The results are compared with an analytical model and experimental data. The numerical simulations provide the damping ratio, an important input, for the analytical model, which was developed to predict the stability limits of technical combustion systems

Smoothed particle hydrodynamics and its applications for multiphase flow and reactive transport in porous media

Smoothed particle hydrodynamics (SPH) is a Lagrangian method based on a meshless discretization of partial differential equations. In this review, we present SPH discretization of the Navier-Stokes and advection-diffusion-reaction equations, implementation of various boundary conditions, and time integration of the SPH equations, and we discuss applications of the SPH method for modeling pore-scale multiphase flows and reactive transport in porous and fractured media.United States. Dept. of Energy. Office of Advanced Scientific Computing Research (Early Career Award, “New Dimension Reduction Methods and Scalable Algorithms for Multiscale Nonlinear Phenomena,” and Collaboratory on Mathematics for Mesoscopic Modeling of Materials (CM4)

Computational fluid dynamics

An overview of computational fluid dynamics (CFD) activities at the Langley Research Center is given. The role of supercomputers in CFD research, algorithm development, multigrid approaches to computational fluid flows, aerodynamics computer programs, computational grid generation, turbulence research, and studies of rarefied gas flows are among the topics that are briefly surveyed

Boundary layer simulator improvement

Boundary Layer Integral Matrix Procedure (BLIMPJ) has been identified by the propulsion community as the rigorous boundary layer program in connection with the existing JANNAF reference programs. The improvements made to BLIMPJ and described herein have potential applications in the design of the future Orbit Transfer Vehicle engines. The turbulence model is validated to include the effects of wall roughness and a way is devised to treat multiple smooth-rough surfaces. A prediction of relaminarization regions is examined as is the combined effects of wall cooling and surface roughness on relaminarization. A turbulence model to represent the effects of constant condensed phase loading is given. A procedure is described for thrust decrement calculation in thick boundary layers by coupling the T-D Kinetics Program and BLIMPJ and a way is provided for thrust loss optimization. Potential experimental studies in rocket nozzles are identified along with the required instrumentation to provide accurate measurements in support of the presented new analytical models

The compressible granular collapse in a fluid as a continuum: validity of a Navier-Stokes model with μ(J)\mu(J)-ϕ(J)\phi(J)-rheology

The incompressible $\mu(I)$-rheology has been used to study subaerial granular flows with remarkable success. For subaquatic granular flows, drag between grains and the pore fluid is substantially higher and the physical behaviour is more complex. High drag forces constrain the rearrangement of grains and dilatancy, leading to a considerable build-up of pore pressure. Its transient and dynamic description is the key to modelling subaquatic granular flows but out of the scope of incompressible models. In this work, we advance from the incompressible $\mu(I)$-rheology to the compressible $\mu(J)$-$\phi(J)$-rheology to account for pore pressure, dilatancy, and the scaling laws under subaquatic conditions. The model is supplemented with critical state theory to yield the correct properties in the quasi-static limit. The pore fluid is described by an additional set of conservation equations and the interaction with grains is described by a drag model. This new implementation enables us to include most of the physical processes relevant for submerged granular flows in a highly transparent manner. Both, the incompressible and compressible rheologies are implemented into OpenFOAM and various simulations at low and high Stokes numbers are conducted with both frameworks. We found a good agreement of the $\mu(J)$-$\phi(J)$-rheology with low Stokes number experiments, that incompressible models fail to describe. The combination of granular rheology, pore pressure, and drag model leads to complex phenomena such as apparent cohesion, remoulding, hydroplaning, and turbidity currents. The simulations give remarkable insights into these phenomena and increase our understanding of subaquatic mass transports