Simulation of three-dimensional viscous flow in turbomachinery geometries using a solution-adaptive unstructured mesh methodology

This paper presents a numerical method for the simulation of flow in turbomachinery blade rows using a solution-adaptive mesh methodology. The fully three-dimensional, compressible, Reynolds-averaged Navier-Stokes equations with k-ε turbulence modeling (and low Reynolds number damping terms) are solved on an unstructured mesh formed from tetrahedral finite volumes. At stages in the solution, mesh refinement is carried out based on flagging cell faces with either a fractional variation of a chosen variable (like Mach number) greater than a given threshold or with a mean value of the chosen variable within a given range. Several solutions are presented, including that for the highly three-dimensional flow associated with the corner stall and secondary flow in a transonic compressor cascade, to demonstrate the potential of the new method.

Influence of a bow of finite width on bowed string motion: numerical modelling and experimental evidence

An enhanced physical model of the bowed string presented previously [1] is explored. It takes into account: the width of the bow, the angular motion of the string, bow-hair elasticity and string bending stiffness. The results of an analytical investigation of a model system - an infinite string sticking to a bow of finite width and driven on one side of the bow - are compared with experimental results published by Cremer [2] and reinterpreted here. Comparison shows that both the width of the bow and the bow-hair elasticity have a large impact on the reflection and transmission behaviour. In general, bending stiffness plays a minor role. Furthermore, a method of numerical simulation of the stiff string bowed with a bow of finite width is presented along with some preliminary results.

Three-dimensional direct numerical simulations of autoignition in turbulent non-premixed flows with simple and complex chemistry

3D Direct Numerical Simulations (DNS) of autoignition in turbulent non-premixed flows between fuel and hotter air have been carried out using both 1-step and complex chemistry consisting of a 22 species n-heptane mechanism to investigate spontaneous ignition timing and location. The simple chemistry results showed that the previous findings from 2D DNS that ignition occurred at the most reactive mixture fraction (ξMR) and at small values of the conditional scalar dissipation rate (N|ξMR) are valid also for 3D turbulent mixing fields. Performing the same simulation many times with different realizations of the initial velocity field resulted in a very narrow statistical distribution of ignition delay time, consistent with a previous conjecture that the first appearance of ignition is correlated with the low-N content of the conditional probability density function of N. The simulations with complex chemistry for conditions outside the Negative Temperature Coefficient (NTC) regime show behaviour similar to the single-step chemistry simulations. However, in the NTC regime, the most reactive mixture fraction is very rich and ignition seems to occur at high values of scalar dissipation. Copyright © 2006 by ASME.

The use of object oriented modelling and simulation tools in constraint satisfaction for embodiment design

The software package Dymola, which implements the new, vendor-independent standard modelling language Modelica, exemplifies the emerging generation of object-oriented modelling and simulation tools. This paper shows how, in addition to its simulation capabilities, it may be used as an embodiment design tool, to size automatically a design assembled from a library of generic parametric components. The example used is a miniature model aircraft diesel engine. To this end, the component classes contain extra algebraic equations calculating the overload factor (or its reciprocal, the safety factor) for all the different modes of failure, such as buckling or tensile yield. Thus the simulation results contain the maximum overload or minimum safety factor for each failure mode along with the critical instant and the device state at which it occurs. The Dymola "Initial Conditions Calculation" function, controlled by a simple software script, may then be used to perform automatic component sizing. Each component is minimised in mass, subject to a chosen safety factor against failure, over a given operating cycle. Whilst the example is in the realm of mechanical design, it must be emphasised that the approach is equally applicable to the electrical or mechatronic domains, indeed to any design problem requiring numerical constraint satisfaction.

CFD-DEM simulation of propagation of sound waves in fluid particles fluidised medium

In this work, speed of sound in 2 phase mixture has been explored using CFD-DEM (Computational Fluid Dynamcis - Discrete Element Modelling). In this method volume averaged Navier Stokes, continuity and energy equations are solved for fluid. Particles are simulated as individual entities; their behaviour is captured by Newton's laws of motion and classical contact mechanics. Particle-fluid interaction is captured using drag laws given in literature.The speed of sound in a medium depends on physical properties. It has been found experimentally that speed of sound drops significantly in 2 phase mixture of fluidised particles because of its increased density relative to gas while maintaining its compressibility. Due to the high rate of heat transfer within 2 phase medium as given in Roy et al. (1990), it has been assumed that the fluidised gas-particle medium is isothermal.The similar phenomenon has been tried to be captured using CFD-DEM numerical simulation. The disturbance is introduced and fundamental frequency in the medium is noted to measure the speed of sound for e.g. organ pipe. It has been found that speed of sound is in agreement with the relationship given in Roy et al. (1990). Their assumption that the system is isothermal also appears to be valid.

Stress-strain response of hydrate-bearing sands: Numerical study using discrete element method simulations

Gas hydrate is a crystalline solid found within marine and subpermafrost sediments. While the presence of hydrates can have a profound effect on sediment properties, the stress-strain behavior of hydrate-bearing sediments is poorly understood due to inherent limitations in laboratory testing. In this study, we use numerical simulations to improve our understanding of the mechanical behavior of hydrate-bearing sands. The hydrate mass is simulated as either small randomly distributed bonded grains or as "ripened hydrate" forming patchy saturation, whereby sediment clusters with 100% pore-filled hydrate saturation are distributed within a hydrate-free sediment. Simulation results reveal that reduced sand porosity and higher hydrate saturation cause an increase in stiffness, strength, and dilative tendency, and the critical state line shifts toward higher void ratio and higher shear strength. In particular, the critical state friction angle increases in sands with patchy saturation, while the apparent cohesion is affected the most when the hydrate mass is distributed in pores. Sediments with patchy hydrate distribution exhibit a slightly lower strength than sediments with randomly distributed hydrate. Finally, hydrate dissociation under drained conditions leads to volume contraction and/or stress relaxation, and pronounced shear strains can develop if the hydrate-bearing sand is subjected to deviatoric loading during dissociation.

A priori assessment of algebraic flame surface density models in the context of large eddy simulation for nonunity lewis number flames in the thin reaction zones regime

The performance of algebraic flame surface density (FSD) models has been assessed for flames with nonunity Lewis number (Le) in the thin reaction zones regime, using a direct numerical simulation (DNS) database of freely propagating turbulent premixed flames with Le ranging from 0.34 to 1.2. The focus is on algebraic FSD models based on a power-law approach, and the effects of Lewis number on the fractal dimension D and inner cut-off scale η i have been studied in detail. It has been found that D is strongly affected by Lewis number and increases significantly with decreasing Le. By contrast, η i remains close to the laminar flame thermal thickness for all values of Le considered here. A parameterisation of D is proposed such that the effects of Lewis number are explicitly accounted for. The new parameterisation is used to propose a new algebraic model for FSD. The performance of the new model is assessed with respect to results for the generalised FSD obtained from explicitly LES-filtered DNS data. It has been found that the performance of the most existing models deteriorates with decreasing Lewis number, while the newly proposed model is found to perform as well or better than the most existing algebraic models for FSD. © 2012 Mohit Katragadda et al.

Direct simulation of turbulent entrainment due to a plume impinging on a density interface

The law for turbulent entrainment due to plumes and jets impinging on a density interface is subject to significant uncertainty, with reported differences in entrainment rates up to a factor of 10. We report preliminary results obtained by Direct Numerical Simulation which are part of a PRACE project on turbulent entrainment carried out on JUGENE at Jülich, Germany. Various interface tracking methods are discussed and the entrainment coefficient is determined.

Boundary integral technique for the numerical modelling of the air flow for the Aaberg exhaust system

A boundary integral technique has been developed for the numerical simulation of the air flow for the Aaberg exhaust system. For the steady, ideal, irrotational air flow induced by a jet, the air velocity is an analytical function. The solution of the problem is formulated in the form of a boundary integral equation by seeking the solution of a mixed boundary-value problem of an analytical function based on the Riemann-Hilbert technique. The boundary integral equation is numerically solved by converting it into a system of linear algebraic equations, which are solved by the process of the Gaussian elimination. The air velocity vector at any point in the solution domain is then computed from the air velocity on the boundary of the solution domains.

Large-eddy simulation of complex geometry jets

Computations are made for chevron and coflowing jet nozzles. The latter has a bypass ratio of 6:1. Also, unlike the chevron nozzle, the core flow is heated, making the inlet conditions reminiscent of those for a real engine. A large-eddy resolving approach is used with circa 12 × 10 6 cell meshes. Because the codes being used tend toward being dissipative the subgrid scale model is abandoned, giving what can be termed numerical large-eddy simulation. To overcome near-wall modeling problems a hybrid numerical large-eddy simulation-Reynolds-averaged Navier-Stokes related method is used. For y + ≤ 60 a Reynolds-averaged Navier-Stokes model is used. Blending between the two regions makes use of the differential Hamilton-Jabobi equation, an extension of the eikonal equation. For both nozzles, results show encouraging agreement with measurements of other workers. The eikonal equation is also used for ray tracing to explore the effect of the mean flow on acoustic ray trajectories, thus yielding a coherent solution strategy. © 2011 by Cambridge University.

Large eddy simulation of a hot, high speed coflowing jet

Computations are made of a short cowl coflowing jet nozzle with a bypass ratio 8 : 1. The core flow is heated, making the inlet conditions reminiscent of those for a real engine. A large eddy resolving approach is used with a 12 × 106 cell mesh. Since the code being used tends towards being dissipative the sub-grid scale (SGS) model is abandoned giving what can be termed Numerical Large Eddy Simulation (NLES). To overcome near wall modelling problems a hybrid NLES-RANS (Reynolds Averaged Navier-Stokes) related method is used. For y+ ≤ 60 a κ-l model is used. Blending between the two regions makes use of the differential Hamilton-Jabobi (HJ) equation, an extension of the eikonal equation. Results show encouraging agreement with existing measurements of other workers. The eikonal equation is also used for acoustic ray tracing to explore the effect of the mean flow on acoustic ray trajectories, thus yielding a coherent solution strategy. Copyright © 2011 by ASME.

Numerical study of soil heterogeneity effects on contaminant transport in unsaturated soil (model development and validation)

The movement of chemicals through soil to groundwater is a major cause of degradation of water resources. In many cases, serious human and stock health implications are associated with this form of pollution. The study of the effects of different factors involved in transport phenomena can provide valuable information to find the best remediation approaches. Numerical models are increasingly being used for predicting or analyzing solute transport processes in soils and groundwater. This article presents the development of a stochastic finite element model for the simulation of contaminant transport through soils with the main focus being on the incorporation of the effects of soil heterogeneity in the model. The governing equations of contaminant transport are presented. The mathematical framework and the numerical implementation of the model are described. The comparison of the results obtained from the developed stochastic model with those obtained from a deterministic method and some experimental results shows that the stochastic model is capable of predicting the transport of solutes in unsaturated soil with higher accuracy than deterministic one. The importance of the consideration of the effects of soil heterogeneity on contaminant fate is highlighted through a sensitivity analysis regarding the variance of saturated hydraulic conductivity as an index of soil heterogeneity. © 2011 John Wiley & Sons, Ltd.