Search

This paper presents image velocimetry measurements on turbulent flows adjacent to a permeable bed made of randomly packed glass particles. For measuring flow velocities inside the bed, the refractive index of the glass particles was matched with that of the fluid. By continuously scanning in the transverse direction, we measured the streamwise and vertical velocity components within a three-dimensional domain (3D2C-PIV), including first- and second-order turbulent statistics. We established how the scanning travel speed is associated with the laser sheet thickness and the space-time velocity fluctuations for collecting reliable measurements. The methodology was applied to free-surface flows over a sloping bed under low relative submergence and supercritical conditions. Space- and time-averaged profiles were obtained in a representative elementary volume as defined by thedouble-averagingprocedure (Nikora et al. in J Hydraulic Eng.127(2):123–133, 2001). A turbulent boundary layer over the rough bed was observed when experiments were run at intermediate Reynolds numbers Re =$$O(1000)$$O(1000). Apart from measuring subsurface velocities, this method shed light on the part played by the rough bed in the overall flow dynamics: the roughness layer was a buffer region within which porosity varied sharply and turbulent stress was rapidly dampened.Graphic abstract

This work predicts the evolution of self-excited thermo-acoustic instabilities in a gas turbine model combustor using large eddy simulation. The applied flow solver is fully compressible and comprises a transported sub-grid probability density function approach in conjunction with the Eulerian stochastic fields method. An unstable operating condition in the PRECCINSTA test case—known to exhibit strong flame oscillations driven by thermo-acoustic instabilities—is the chosen target configuration. Good results are obtained in a comparison of time-averaged flow statistics against available measurement data. The flame’s self-excited oscillatory behaviour is successfully captured without any external forcing. Power spectral density analysis of the oscillation reveals a dominant thermo-acoustic mode at a frequency of 300 Hz; providing remarkable agreement with previous experimental observations. Moreover, the predicted limit-cycle amplitude is found to closely match its respective measured value obtained from experiments with rigid metal combustion chamber side walls. Finally, a phase-resolved study of the oscillation cycle is carried out leading to a detailed description of the physical mechanisms that sustain the closed feedback loop.

In this study, free vibration analyses of embedded carbon and silica carbide nanotubes lying on an elastic matrix are performed based on Eringen’s nonlocal elasticity theory. These nanotubes are modeled as nanobeam and nanorod. Elastic matrix is considered as Winkler–Pasternak elastic foundation and axial elastic media for beam and rod models, respectively. The vibration formulations of the beam and rod are derived by utilizing Hamilton’s principle. The obtained equations of motions are solved by the method of separation of variables and finite element-based Hermite polynomials for various boundary conditions. The effects of boundary conditions, system modeling, structural sizes such as length, cross-sectional sizes, elastic matrix, mode number, and nonlocal parameters on the natural frequencies of these nanostructures are discussed in detail. Moreover, the availability of size-dependent finite element formulation is investigated in the vibration problem of nanobeams/rods resting on an elastic matrix.

A backward facing micro-step is a building block for many microfluidic devices. Due to micron sized characteristic dimensions, the gas flow in such a geometry is rarefied in nature. Such rarefied gas flows are widely solved using the Direct Simulation Monte Carlo (DSMC) technique. Flow separation, circulation and re-attachment are some of the basic characteristics of step flows. The objective of this study is to analyze the effect of rarefaction on the flow properties and the separation of the flow. The range of selected Knudsen number (Kn) covers the slip and transition regime from a value of 0.0311–13.25. The pressure ratios employed are 3 and 5. It is observed that the slip velocity continuously increases while the centre-line velocity first decreases, then remains constant and finally increases with increase in Kn. At the step, separation of the flow is seen for Kn 0.1325 while no such separation is observed in the range of Kn from 0.198 to 13.25. The corresponding Re for these ranges are 6.43 to 0.67 and 0.392 to 0.012 respectively. The re-attachment length decreases with increase in Kn whereas it increases with increase in Re. A stronger pressure force and a weaker diffusion effect leads to flow separation in the slip regime whereas stronger diffusion and weaker pressure force lead to an absence of flow separation in the transition regime. Finally, this work presents for the first time the existence of the Knudsen minimum for such a backward step geometry.

This paper presents an experimental study of heat transfer in a slightly rarefied gas flowing in a circular tube with constant wall temperature boundary condition. Local temperature measurements are carried out for the first time in rarefied gas flows to investigate into the anomalous values of Nusselt number (Nu) obtained by a previous experimental study (Demsis et al., 2009). The present measurements agree well with the low Nu reported by Demsis et al. (2009) when the Nu is obtained using their procedure; additionally, the measurements reveal the importance of end effects in determining the Nusselt number in rarefied gases. The Nusselt number obtained in the present experiments tends to zero with increasing axial conduction. Nu shows a weak dependence on Knudsen number for the range investigated here (0.001 < Kn < 0.012).

This study concerns with the free vibration analyses for centrifugally stiffened cracked beams modelled as functionally graded material in the axial direction. The variation of material properties along the cracked beam is expressed in terms of the power law distribution. Analyses are carried out by the Rayleigh–Ritz method that uses shape functions and energy expressions written for centrifugally stiffened Euler-Bernoulli beams. The influence of the hub radius, angular velocity, crack location to beam length ratio and crack depth to thickness ratio and cross-section size ratio and inhomogeneity on natural frequencies are examined for the beams with the clamped-free boundary condition. The results are verified by using the data available in the open literature and/or the outputs of finite element analyses performed for an axially functionally graded solid beam. The results show that the fundamental frequency parameters reduce with the increasing the ratio of crack depth to thickness. The highest natural frequency droop exists close to the root of the beam for flap-wise and chord-wise vibrations. Another important finding is that reduction of natural frequency parameters in chord-wise modes is larger than that in flap-wise modes when the crack location is varied. This new research should help to improve predictions of the impact of the crack on free vibrations for the rotating cracked axially functionally graded beams.

Numerical blood flow simulations are typically set up from anatomical medical images and calibrated using velocity measurements. However, the accuracy of the computational geometry itself is limited by the resolution of the anatomical image. We first show that applying standard no‐slip boundary conditions on inaccurately extracted boundaries can cause large errors in the results, in particular the pressure gradient. In this work, we therefore propose to augment the flow model calibration by slip/transpiration boundary conditions, whose parameters are then estimated using velocity measurements. Numerical experiments show that this methodology can considerably improve the accuracy of the estimated pressure gradients and 3D velocity fields when the vessel geometry is uncertain.

This paper presents a comparative study of numerical, experimental and empirical techniques on the effect of localised air blast loads on mild steel and armour steel plates. The blast load effects on monolithic plates have been accounted for by using different approaches provided in the Finite Element hydrocode ABAQUS 6.13, namely an Eulerian Lagrangian and a Coupled Eulerian Lagrangian model. In the first model, the air and the explosive were modelled using multi-material Eulerian grids while the plate was modelled using a rigid Lagrangian mesh, while in the second model the rigid target was replaced with deformable plate.\ud \ud \ud \ud The transient deformation of the plate, strain localisation, pressure distribution on the plate have been investigated in the FE models, which have been validated against small scale experimental data for a limited range of charge sizes for both the mild steel and armoured steel. Despite the lower deflection of armour steel compared to mild steel plates, both plates were shown to undergo rupture upon similar charge mass and stand-off. For this purpose, a non-dimensional analysis was carried out with consideration of stand-off distance and slenderness ratio to predict the rupture impulse.

3D Numerical study of temperature variation for subsonic rarefied gas flow in a microchannel is carried out using an in-house MPI-based parallelized DSMC code. The temperature drop in the microchannel decreases with an increase in the aspect ratio whereas it increases with an increase in the pressure ratio, the cross-aspect ratio (CAR), and the Knudsen number. 3D and 2D simulations results are compared and effect of the CAR and Knudsen number are brought out. Finally, a correlation that predicts the temperature drop is formulated along with a list of conditions that ensures a near isothermal flow.

High temperature gradients in heat sinks of micro devices cause substantial properties variations in gas-micro-convection; which constitutes an important strategic research area in non-rarefaction scaling effects. In present analysis, compressible Navier-Stokes equations incorporating temperature dependent density and thermal conductivity variations, for steady, viscous flow of gas at subsonic speeds (Mach number << 1) are numerically solved. Circular micro-pipe geometry, subjected to constant wall heat flux was chosen, The key observations are: In addition to known effects, two additional physical mechanisms increasingly surface at micro scale; and also determine the micro convection characteristics within continuum regime. They are, induced radial convection due to density variations in radial and axial directions and axial conduction induced due to thermal conductivity variations along the flow. High density variations in radial and axial directions, cause velocity gradients and increase fluid friction irreversibility. High thermal conductivity variations cause flattening of temperature profile, induce axial conduction and substantially increase in entropy generation. The results highlight the need for incorporating fluid properties variation in thermal designs of heat sinks for micro devices. Two-way link between velocity profile and temperature, influence entropy generation. Corner effects at entrance and exit influence entropy generation in all cases studied. (C) 2014 Elsevier Ltd. All rights reserved.

This article presents a three-dimensional numerical framework for the simulation of fluid-fluid immiscible compounds in complex geometries, based on the multiple-relaxation-time lattice Boltzmann method to model the fluid dynamics and the color-gradient approach to model multicomponent flow interaction. New lattice weights for the lattices D3Q15, D3Q19, and D3Q27 that improve the Galilean invariance of the color-gradient model as well as for modeling the interfacial tension are derived and provided in the Appendix. The presented method proposes in particular an approach to model the interaction between the fluid compound and the solid, and to maintain a precise contact angle between the two-component interface and the wall. Contrarily to previous approaches proposed in the literature, this method yields accurate solutions even in complex geometries and does not suffer from numerical artifacts like nonphysical mass transfer along the solid wall, which is crucial for modeling imbibition-type problems. The article also proposes an approach to model inflow and outflow boundaries with the color-gradient method by generalizing the regularized boundary conditions. The numerical framework is first validated for three-dimensional (3D) stationary state (Jurin's law) and time-dependent (Washburn's law and capillary waves) problems. Then, the usefulness of the method for practical problems of pore-scale flow imbibition and drainage in porous media is demonstrated. Through the simulation of nonwetting displacement in two-dimensional random porous media networks, we show that the model properly reproduces three main invasion regimes (stable displacement, capillary fingering, and viscous fingering) as well as the saturating zone transition between these regimes. Finally, the ability to simulate immiscible two-component flow imbibition and drainage is validated, with excellent results, by numerical simulations in a Berea sandstone, a frequently used benchmark case used in this field, using a complex geometry that originates from a 3D scan of a porous sandstone. The methods presented in this article were implemented in the open-source PALABOS library, a general C++ matrix-based library well adapted for massive fluid flow parallel computation.

International audience; In this work, a nonlocal quasi-3D trigonometric plate theory for micro/nanoscale plates is proposed. In order to introduce the size influences, the Eringen's nonlocal elasticity theory is utilized. In addition, the theory considers both shear deformation and thickness stretching effects by a trigonometric variation of all displacements within the thickness, and respects the stress-free boundary conditions on the top and bottom surfaces of the plate without considering the shear correction factor. The advantage of this theory is that, in addition to considering the small scale and thickness stretching effects (epsilon(z)not equal 0), the displacement field is modelled with only 5 unknowns as the first order shear deformation theory (FSDT). Analytical solutions for vibration of simply supported micro/nanoscale plates are illustrated, and the computed results are compared with the available solutions in the literature and finite element model using ABAQUS software package. The influences of the nonlocal parameter, shear deformation and thickness stretching on the vibration behaviors of the micro/ nanoscale plates are examined

International audience; The analysis of self-excited combustion instabilies encountered in a laboratory-scale, swirl-stabilized combustion system is presented. The instability is successfully captured by reactive large-eddy simulation (LES) and analyzed by using a global acoustic energy equation. This energy equation shows how the source term due to combustion (equivalent to the Rayleigh criterion) is balanced by the acoustic fluxes at the boundaries when reaching the limit cycle. Additionally, an Helmholtz-equation solver including flame-acoustics interaction modeling is used to predict the stability characteristics of the system. Feeding the flame-transfer function from the LES into this solver allows to predict an amplification rate for each mode. The unstable mode encountered in the LES compares well with the mode of the highest amplification factor in the Helmholtz-equation solver, in terms of mode shape as well as in frequency.

In this paper, some aspects of water impact and green water loading are considered by numerically investigating a dambreak problem and water entry problems. The numerical method is based on the Navier-Stokes equations that describe the flow of an incompressible viscous fluid. The equations are discretised on a fixed Cartesian grid using the finite volume method. Even though very small cut cells can appear when moving an object through the fixed grid, the method is stable. The free surface is displaced using the Volume-of-Fluid method together with a local height function, resulting in a strictly mass conserving method. The choice of boundary conditions at the free surface appears to be crucial for the accuracy and robustness of the method. For validation, results of a dambreak simulation are shown that can be compared with measurements. A box has been placed in the flow, as a model for a container on the deck of an offshore floater on which forces are calculated. The water entry problem has been investigated by dropping wedges with different dead-rise angles, a cylinder and a cone into calm water with a prescribed velocity. The resulting free surface dynamics, with the sideways jets, has been compared with photographs of experiments. Also a comparison of slamming coefficients with theory and experimental results has been made. Finally, a drop test with a free falling wedge has been simulated. (c) 2005 Elsevier Inc. All rights reserved.

Investigation of combustion instabilities in gas turbine combustors require the knowledge of flame transfer functions. Those can be obtained by experimental measurement or by Large Eddy Simulations (LES). Because calculations are usually limited to a portion of the whole combustor, boundary conditions are of crucial importance. It is common practice to inject acoustic perturbations for the flame transfer function measurement in form of velocity perturbations (u(t)). We present an alternative method based on a characteristic treatment of the Euler Equations. It consists of injecting sound waves traveling into the computational inlet while letting outgoing waves leave the domain without reflection. This method has several advantages concerning the study of flame transfer functions compared to injecting velocity perturbations. Both techniques are compared for cases where analytical solutions may be derived (a duct without flame and a planar laminar flame) and for one case where a CFD code is necessary (a laminar Bunsen-type flame). © 2002 by The Combustion Institute

This paper presents an analysis of the energy transported by disturbances in gaseous combustion. It extends the previous work of Myers (J. Fluid Mech., vol. 226, 1991, 383–400) and so includes non-zero mean-flow quantities, large-amplitude disturbances, varying specific heats and chemical non-equilibrium. This extended form of Myers’ ‘disturbance energy’ then enables complete identification of the conditions under which the famous Rayleigh source term can be derived from the equations governing combusting gas motion. These are: small disturbances in an irrotational, homentropic, non-diffusive (in terms of species, momentum and energy) and stationary mean flow at chemical equilibrium. Under these assumptions, the Rayleigh source term becomes the sole source term in a conservation equation for the classical acoustic energy. It is also argued that the exact disturbance energy flux should become an acoustic energy flux in the far-field surrounding a (reacting or non-reacting) jet. In this case, the volume integral of the disturbance energy source terms are then directly related to the area-averaged far-field sound produced by the jet. This is demonstrated by closing the disturbance energy budget over a set of aeroacoustic, direct numerical simulations of a forced, low-Mach-number, laminar, premixed flame. These budgets show that several source terms are significant, including those involving the mean-flow and entropy fields. This demonstrates that the energetics of sound generation cannot be examined by considering the Rayleigh source term alone.

We introduce a novel multiple-relaxation time ( modified MRT) Lattice Boltzmann scheme for simulation of confined suspension flow. Via careful tuning of the free eigenvalues of the collision operator we can substantially reduce the error in the so-called hydrodynamic radius. Its performance has been compared to that of the TRT scheme for several benchmark problems. We have found that the optimal value of the free eigenvalue depends on the curvature of the solid–fluid interfaces. Hence, we have investigated suspension flow problems, with confining boundaries of different curvatures. We have found that the modified MRT scheme is better suited for suspension flow in curved confining walls, while the TRT scheme is better for suspension flow confined between planar walls. With both schemes we have investigated problems for confined suspension flows, namely 1) drag forces experienced by spheres flowing in confining flow channels of different cross sections, and 2) the lubrication force between a sedimenting sphere and the end cap of a confining cylindrical capillary.

A thorough study of shear stress within the lattice Boltzmann method is provided. Via standard multiscale Chapman-Enskog expansion we investigate the dependence of the error in shear stress on grid resolution showing that the shear stress obtained by the lattice Boltzmann method is second order accurate. This convergence, however, is usually spoiled by the boundary conditions. It is also investigated which value of the relaxation parameter minimizes the error. Furthermore, for simulations using velocity boundary conditions, an artificial mass increase is often observed. This is a consequence of the compressibility of the lattice Boltzmann fluid. We investigate this issue and derive an analytic expression for the time-dependence of the fluid density in terms of the Reynolds number, Mach number and a geometric factor for the case of a Poiseuille flow through a rectangular channel in three dimensions. Comparison of the analytic expression with results of lattice Boltzmann simulations shows excellent agreement., Comment: 15 pages, 4 figures, 2 tables

International audience; An adiabatic homogeneous model to account for multiperforated liners in combustion chamber flow simulations is described. It is based on a suction and an injection model to reproduce the average effect of effusion cooling on both sides of the plate. The coupled suction/injection model has been specifically designed to be used in industrial full-scale computations of gas turbine combustion chambers, where effusion cooling is commonly used for controlling the temperature of the liners. Notably, it can be used with a coarse grid, the real perforated plate being replaced by a homogeneous boundary condition where the model is applied. The new modeled boundary condition conserves the inviscid part of the wall fluxes, which are shown to be the main contribution, as evidenced by the analysis of former wall-resolved simulations. Conserving the wall fluxes allows reproduction of the global structure of the flow and leads to reasonable comparisons with experimental data. Hence, the proposed new model provides a practical way to account for multiperforated plates with inclined perforations without resolving the flow in the perforations.

International audience; Large-eddy simulations of a generic turbulent flow with discrete effusion are reported. The computational domain is periodic in both streamwise and spanwise directions and contains both the injection and the suction sides. The blowing ratio is close to 1.2 while the Reynolds number in the aperture is of order 2600. The numerical results for this fully developed bi-periodic turbulent flow with effusion are compared to available experimental data from a large-scale spatially evolving isothermal configuration. It is shown that many features are shared by the two flow configurations. The main difference is related to the mean streamwise velocity profile, which is more flat for the bi-periodic situation where the cumulative effect of an infinite number of upstream jets is accounted for. The necessity of considering both sides of the plate is also established by analysing the vortical structure of the flow and some differences with the classical jet-in-crossflow case are highlighted. Finally, the numerical results are analysed in terms of wall modelling for full-coverage film cooling. For the operating point considered, it is demonstrated that the streamwise momentum flux is dominated by non-viscous effects, although the area where only the viscous shear stress contributes is very large given the small porosity value (4%).

We consider the coupling between three-dimensional (3D) and one-dimensional (1D) fluid- structure interaction (FSI) models describing blood flow inside compliant vessels. The 1D model is a hyperbolic system of partial differential equations. The 3D model consists of the Navier-Stokes equations for incompressible Newtonian fluids coupled with a model for the vessel wall dynamics. A non standard formulation for the Navier-Stokes equations is adopted to have suitable boundary conditions for the coupling of the models. With this we derive an energy estimate for the fully 3D- 1D FSI coupling. We consider several possible models for the mechanics of the vessel wall in the 3D problem and show how the 3D-1D coupling depends on them. Several comparative numerical tests illustrating the coupling are presented.

The mathematical formulation of large-eddy simulation is introduced and the closure problem (subgrid scale modeling) is briefly reviewed. Some fields which need further investigation, in order that large-eddy simulation becomes a completely reliable tool, in particular for the analysis of complex flows, are identified and discussed. Specific problems for application to jet flows are also illustrated. Finally, examples of large-eddy simulations of jets in cross flow are presented.