Comparison of two-and three-dimensional modeling of invert trap for sewer solid management

In the present study, five different invert trap configurations (rectangular with and without lids on both sides; trapezoidal, trapezoidal with rectangular base and rectangular with trapezoidal base with lids on both sides) were simulated for both two-dimensional (2D) and three-dimensional (3D) flow conditions for three sediment types (sand, styrocell and plastic beads) at six flow rates (0.35, 0.70, 1.05, 1.35, 4.55 and 9.95L/s) for each trap. Computational fluid dynamics (CFD)-based modeling using FLUENT software with Renormalization Group (RNG) k-ε model along with discrete phase model (DPM) were used in the simulations. A hexagonal/tetrahedral and map-type non-uniform grid was chosen to discretize the entire computational domain and a control volume finite difference method was used to solve the governing equations. The flow rates selected in the present study cover the entire range of flow rate expected for dry weather and monsoon. The simulation is capable of differentiating between 2D and 3D modeling of particle trajectories, the effects of flow rate and trap geometry on flow patterns developed in the trap. The sediment retention ratio for 2D is higher than that for 3D modeling for all flow conditions, particle types and model geometry due to inclusion of lateral effects in 3D modeling. The invert trap having rectangular shape with trapezoidal base is found to be the most efficient configuration in both 2D and 3D modeling.     

Effective flow of a viscous liquid through a helical pipe

We study the flow of a viscous fluid through a pipe with helical shape parameterized with , where the small parameter stands for the distance between two coils of the helix. The pipe has small cross-section of size . Using the asymptotic analysis of the microscopic flow described by the Navier–Stokes system, with respect to the small parameter that tends to zero, we find the effective fluid flow described by an explicit formula of the Poisseuile type including a small distorsion due to the particular geometry of the pipe. To cite this article: E. Marušić-Paloka, I. Pažanin, C. R. Mecanique 332 (2004).

Résumé

On considère un écoulement dans un tube de section circulaire et de forme hélicoïdale paramétré par , où est la distance entre deux tours de la spirale. Le rayon de la section du tube est lui aussi supposé égal à . A partir de l'écoulement microscopique décrit par le système de Navier–Stokes et en utilisant l'analyse asymptotique par rapport à ce petit paramètre on obtient l'écoulemment effectif décrit par une formule explicite de type Poiseuille associée à une petite déviation due à la géometrie du tube. Pour citer cet article : E. Marušić-Paloka, I. Pažanin, C. R. Mecanique 332 (2004).     

CFD-simulation of a circulating fluidized bed riser

In the current work, a model of the fluid mechanics in the riser of a circulating fluidized bed （CFB） has been implemented using computational fluid dynamics （CFD）. The model developed shall be used in future as the basis of 3D-reactor model for the simulation of large scale CFB combustors. The two-fluid model （TFM） approach is used to represent the fluid mechanics involved in the flow. The computational implementation is accomplished by the commercial software FLUENT. Different closure formulations are tested on a simplified geometry. Two different turbulence formulations, namely the swirl modified RNG k-e model and the Realizable k-e model, are tested in combination with two different approaches to solid phase turbulence, namely the dispersion and per phase approach. One focus of the current work is put on the study of different drag correlations. Besides the drag correlations by Syamlal et al. [Syamlal, M., Rogers, W., ＆ O＇Brien, T. J. （1993）. MFIX documentation theory guide. Technical Report DOE/METC-9411004, U.S. Department of Energy （DOE）. Morgantown Energy Technology Center： Morgantown, WV] and Gidaspow [Gidaspow, D. （1994）. Multiphaseflow andfluidization. New York： Academic Press] the EMMS model has been used to determine the momentum exchange between the two phases. The resulting formulation is then used to simulate a 1-m × 0.3-m cold CFB setup and is validated by experimental results [Schlichtharle, P. （2000）. Fluid dynamics and mixing of solids and gas in the bottom zone of circulating fluidized beds. Unoublished doctoral dissertation, Technische Universitaet Hamburg-Harburg, Shaker Verlag： Aachen].     

A two-dimensional effective model describing fluid–structure interaction in blood flow: analysis, simulation and experimental validation

We derive a closed system of effective equations describing a time-dependent flow of a viscous incompressible Newtonian fluid through a long and narrow elastic tube. The 3D axially symmetric incompressible Navier–Stokes equations are used to model the flow. Two models are used to describe the tube wall: the linear membrane shell model and the linearly elastic membrane and the curved, linearly elastic Koiter shell model. We study the behavior of the coupled fluid–structure interaction problem in the limit when the ratio between the radius and the length of the tube, , tends to zero. We obtain the reduced equations that are of Biot type with memory. An interesting feature of the reduced equations is that the memory term explicitly captures the viscoelastic nature of the coupled problem. Our model provides significant improvement over the standard 1D approximations of the fluid–structure interaction problem, all of which assume an ad hoc closure assumption for the velocity profile. We performed experimental validation of the reduced model using a mock circulatory flow loop assembled at the Cardiovascular Research Laboratory at the Texas Heart Institute. Experimental results show excellent agreement with the numerically calculated solution. Major applications include blood flow through large human arteries. To cite this article: S. Čanić et al., C. R. Mecanique 333 (2005).     

Effect of dislocation structure on fracture toughness of strained BCC-metals

Mechanical properties of materials depend on their structure. Examined are the effects of dislocation structure on fracture toughness and mechanisms of fracture of BCC-metals. Fracture toughness was determined by depending specimens with cracks introduced into the plane perpendicular to the plane of rolling. Fracture toughness increases with decreasing yield stress (for =15–25%). This is due to instability of slightly misoriented cell structure under repeated loading. The peak of fracture toughness at the temperature 77 K was not observed. The increase of fracture toughness for high strained metals (>60% for Mo, and >85% for Cr) corresponds to cell size reduction and the change of fracture mechanisms.     

Domain perturbation method and shape of a bubble in a uniform flow of an inviscid liquid

In many problems with a free boundary there is defined a small parameter, , for which the solution is sometimes known for a particular value, =0, and the general solution is obtained as a series in the parameter. To find this solution, the equations can be written on a reference configuration and solved in a fixed domain. The purpose of this study is to show that this method of domain perturbation is a good one. The range of validity of this method will be studied on the model example of the irrotational flow of a perfect fluid around a bubble. The radius of convergence of the series solution will be determined, as will the nature of the solution in the neighbourhood of the first real singularity.     

Large amplitude oscillatory elongation flow

A filament stretching rheometer (FSR) was used for measuring the elongation flow with a large amplitude oscillative elongation imposed upon the flow. The large amplitude oscillation imposed upon the elongational flow as a function of the time t was defined as where ε is the Hencky strain, is a constant elongational rate for the base elongational flow, Λ the strain amplitude (Λ ≥ 0), and Ω the strain frequency. A narrow molecular mass distribution linear polystyrene with a molecular weight of 145 kg/mol was subjected to the oscillative flow. The onset of the steady periodic regime is reached at the same Hencky strain as the onset of the steady elongational viscosity ( Λ = 0). The integral molecular stress function formulation within the ‘interchain pressure’ concept agrees qualitatively with the experiments.     

Dynamical stability of the response of oscillators with discontinuous or steep first derivative of restoring characteristic

The influence of factors which can lead to incorrect prediction of dynamical stability of the periodic response of oscillators which contain a non-linear restoring characteristic with discontinuous or steep first derivative is considered in this paper. For that purpose, a simple one degree-of-freedom system with a piecewise-linear force-displacement relationship subjected to a harmonic excitation is analysed. Stability of the periodic response obtained in the frequency domain by the incremental harmonic balance method is determined by using the Floquet–Liapounov theorem. Responses in the time domain are obtained by digital simulation. The accuracy of determining the eigenvalues of the monodromy matrix (in the considered example) significantly depend on the corrective vector norm r , the accuracy of numerical determination of the times when the system undergoes a stiffness change, and on the number of step functions M (used in the Hsu's procedure), only for r >1×10⁻⁵,>1×10⁻⁵ and M<2000. Otherwise, except if the maximum modulus of the eigenvalues of the monodromy matrix is very close to unity, their influence on estimation of dynamical stability is minor. On the contrary, neglecting very small harmonic terms of the actual time domain response can cause a very large error in the evaluation of the eigenvalues of the monodromy matrix, and so they can lead to incorrect prediction of the dynamical stability of the solution, regardless of whether the maximum modulus of the eigenvalues of the monodromy matrix is close to unity or not.     

Spanwise effects on instabilities of compressible flow over a long rectangular cavity

The stability properties of two-dimensional (2D) and three-dimensional (3D) compressible flows over a rectangular cavity with length-to-depth ratio of \(L/D=6\) are analyzed at a free-stream Mach number of \(M_\infty =0.6\) and depth-based Reynolds number of \(Re_D=502\). In this study, we closely examine the influence of three-dimensionality on the wake mode that has been reported to exhibit high-amplitude fluctuations from the formation and ejection of large-scale spanwise vortices. Direct numerical simulation (DNS) and bi-global stability analysis are utilized to study the stability characteristics of the wake mode. Using the bi-global stability analysis with the time-averaged flow as the base state, we capture the global stability properties of the wake mode at a spanwise wavenumber of \(\beta =0\). To uncover spanwise effects on the 2D wake mode, 3D DNS are performed with cavity width-to-depth ratio of \(W/D=1\) and 2. We find that the 2D wake mode is not present in the 3D cavity flow with \(W/D=2\), in which spanwise structures are observed near the rear region of the cavity. These 3D instabilities are further investigated via bi-global stability analysis for spanwise wavelengths of \(\lambda /D=0.5{-}2.0\) to reveal the eigenspectra of the 3D eigenmodes. Based on the findings of 2D and 3D global stability analysis, we conclude that the absence of the wake mode in 3D rectangular cavity flows is due to the release of kinetic energy from the spanwise vortices to the streamwise vortical structures that develops from the spanwise instabilities.     

The Asymptotic Finite-dimensional Character of a Spectrally-hyperviscous Model of 3D Turbulent Flow

We obtain attractor and inertial-manifold results for a class of 3D turbulent flow models on a periodic spatial domain in which hyperviscous terms are added spectrally to the standard incompressible Navier–Stokes equations (NSE). Let P _m be the projection onto the first m eigenspaces of A =−Δ, let μ and α be positive constants with α ≥3/2, and let Q _m =I − P _m , then we add to the NSE operators μ A _φ in a general family such that A _φ≥Q _m A ^α in the sense of quadratic forms. The models are motivated by characteristics of spectral eddy-viscosity (SEV) and spectral vanishing viscosity (SVV) models. A distinguished class of our models adds extra hyperviscosity terms only to high wavenumbers past a cutoff λ _m0 where m ₀ ≤ m, so that for large enough m ₀ the inertial-range wavenumbers see only standard NSE viscosity. We first obtain estimates on the Hausdorff and fractal dimensions of the attractor (respectively and ). For a constant K _α on the order of unity we show if μ ≥ ν that and if μ ≤ ν that where ν is the standard viscosity coefficient, l ₀ = λ₁^−1/2 represents characteristic macroscopic length, and is the Kolmogorov length scale, i.e. where is Kolmogorov’s mean rate of dissipation of energy in turbulent flow. All bracketed constants and K _α are dimensionless and scale-invariant. The estimate grows in m due to the term λ _m /λ₁ but at a rate lower than m ^3/5, and the estimate grows in μ as the relative size of ν to μ. The exponent on is significantly less than the Landau–Lifschitz predicted value of 3. If we impose the condition , the estimates become for μ ≥ ν and for μ ≤ ν. This result holds independently of α, with K _α and c _α independent of m. In an SVV example μ ≥ ν, and for μ ≤ ν aspects of SEV theory and observation suggest setting for 1/c within α orders of magnitude of unity, giving the estimate where c _α is within an order of magnitude of unity. These choices give straight-up or nearly straight-up agreement with the Landau–Lifschitz predictions for the number of degrees of freedom in 3D turbulent flow with m so large that (e.g. in the distinguished-class case for m ₀ large enough) we would expect our solutions to be very good if not virtually indistinguishable approximants to standard NSE solutions. We would expect lower choices of λ _m (e.g. with a > 1) to still give good NSE approximation with lower powers on l ₀/l _ε, showing the potential of the model to reduce the number of degrees of freedom needed in practical simulations. For the choice , motivated by the Chapman–Enskog expansion in the case m = 0, the condition becomes , giving agreement with Landau–Lifschitz for smaller values of λ _m then as above but still large enough to suggest good NSE approximation. Our final results establish the existence of a inertial manifold for reasonably wide classes of the above models using the Foias/Sell/Temam theory. The first of these results obtains such an of dimension N > m for the general class of operators A _φ if α > 5/2. The special class of A _φ such that P _m A _φ = 0 and Q _m A _φ ≥ Q _m A ^α has a unique spectral-gap property which we can use whenever α ≥ 3/2 to show that we have an inertial manifold of dimension m if m is large enough. As a corollary, for most of the cases of the operators A _φ in the distinguished-class case that we expect will be typically used in practice we also obtain an , now of dimension m ₀ for m ₀ large enough, though under conditions requiring generally larger m ₀ than the m in the special class. In both cases, for large enough m (respectively m ₀), we have an inertial manifold for a system in which the inertial range essentially behaves according to standard NSE physics, and in particular trajectories on are controlled by essentially NSE dynamics.      

Strong Solutions for Generalized Newtonian Fluids

We consider the motion of a generalized Newtonian fluid, where the extra stress tensor is induced by a potential with p-structure (p = 2 corresponds to the Newtonian case). We focus on the three dimensional case with periodic boundary conditions and extend the existence result for strong solutions for small times from \tfrac{5}{3}$$ " align="middle" border="0"> (see [16]) to \tfrac{7}{5}.$$ " align="middle" border="0"> Moreover, for we improve the regularity of the velocity field and show that for all 0.$$ " align="middle" border="0"> Within this class of regularity, we prove uniqueness for all \tfrac{7}{5}.$$ " align="middle" border="0"> We generalize these results to the case when p is space and time dependent and to the system governing the flow of electrorheological fluids as long as      

Streamline Topologies in Stokes Flow Within Lid-Driven Cavities

Stokes flow in a rectangular cavity with two moving lids (with equal speed but in opposite directions) and aspect ratio A (height to width) is considered. An analytic solution for the streamfunction, , expressed as an infinite series of Papkovich–Fadle eigenfunctions is used to reveal changes in flow structures as A is varied. Reducing A from A=0.9 produces a sequence of flow transformations at which a saddle stagnation point changes to a centre (or vice versa) with the generation of two additional stagnation points. To obtain the local flow topology as A0, we expand the velocity field about the centre of the cavity and then use topological methods. Expansion coefficients depend on the cavity aspect ratio which is considered as a separation parameter. The normal-form transformations result in a much simplified system of differential equations for the streamlines encapsulating all features of the original system. Using the simplified system, streamline patterns and their bifurcations are obtained, as A0.     

Three-dimensional computational fluid dynamics (volume of fluid) modelling coupled with a stochastic discrete phase model for the performance analysis of an invert trap experimentally validated using field sewer solids

Invert traps are used to trap sewer solids flowing into a sewer drainage system. The performance of the invert trap in an open rectangular channel was experimentally and numerically analysed using field sewer solids collected from a sewer drain. Experiments showed that the free water surface rises over the central opening (slot) of the invert trap, which reduces the velocity near the slot and allows more sediment to be trapped in comparison with the case for the fixed-lid model (assuming closed conduit flow with a shear-free top wall) used by earlier investigators. This phenomenon cannot be modelled using a closed conduit model as no extra space is provided for the fluctuation of the water surface, whereas this space is provided in the volume of fluid (VOF) model in the form of air space in ANSYS Fluent 14.0 software. Additionally, the zero atmospheric pressure at the free water surface cannot be modelled in a fixed-lid model. In the present study, experimental trap efficiencies of the invert trap using field sewer solids were fairly validated using a three-dimensional computational fluid dynamics model (VOF model) coupled with a stochastic discrete phase model. The flow field (i.e., velocities) predicted by the VOF model were compared with experimental velocities obtained employing particle image velocimetry. The water surface profile above the invert trap predicted by the VOF model was found to be in good agreement with the experimentally measured profile. The present study thus showed that the VOF model can be used with the stochastic discrete phase model to well predict the performance of invert traps.     

On the Darcy–Brinkman Flow Through a Channel with Slightly Perturbed Boundary

The goal of this paper is to study the effects of a slightly perturbed boundary on the Darcy–Brinkman flow through a porous channel. We start from a rectangular domain and then perturb the upper part of its boundary by the product of the small parameter \(\epsilon \) and arbitrary smooth function h. Using asymptotic analysis with respect to \(\epsilon \), the effective model has been formally derived. Being in the form of the explicit formulae for the velocity and pressure, the asymptotic approximation clearly shows the nonlocal effects of the small boundary perturbation. The error analysis is also conducted providing the order of accuracy of the asymptotic solution.     

On the Errors Associated with Two-Dimensional Stochastic Solute Transport Models

Many subsurface solute transport studies employ numerical modeling techniques to estimate solute arrival times. Simplifying assumptions must be made to define the modeling domain within a mathematical framework. One common assumption is that the vertical flow is negligible such that the flow field can be simulated with a two-dimensional model. Reducing the vertical dimension reduces the number of flow paths that a solute can take. In a heterogenous medium, artificially removing the 3rd dimension may lead to erroneous results. We investigate the error in the simulated solute breakthrough associated with a two-dimensional model. We also use a stochastic solution of solute arrival time to derive a transform of a two-dimensional ln (k) field so that solute transport more closely resembles three-dimensional transport behavior. The moment equations for two- and three-dimensional domains were solved simultaneously to calculate this transform. The results indicate that the removal of the vertical variability (3D 2D) introduces a 5–10% error in the predicted solute breakthrough. The error tends to increase with increased hydraulic conductivity variance. Numerical experiments confirm that the transform developed herein decreases the relative error of particle breakthrough curves.     

Random Kick-Forced 3D Navier–Stokes Equations in a Thin Domain

We consider the Navier–Stokes equations in the thin 3D domain , where is a two-dimensional torus. The equation is perturbed by a non-degenerate random kick force. We establish that, firstly, when ε ≪ 1, the equation has a unique stationary measure and, secondly, after averaging in the thin direction this measure converges (as ε → 0) to a unique stationary measure for the Navier–Stokes equation on . Thus, the 2D Navier–Stokes equations on surfaces describe asymptotic in time, and limiting in ε, statistical properties of 3D solutions in thin 3D domains.     

Progress in subgrid-scale transport modelling for continuous hybrid non-zonal RANS/LES simulations

The partially integrated transport modelling (PITM) method can be viewed as a continuous approach for hybrid RANS/LES modelling allowing seamless coupling between the RANS and the LES regions. The subgrid turbulence quantities are thus calculated from spectral equations depending on the varying spectral cutoff location [Schiestel, R., Dejoan, A., 2005. Towards a new partially integrated transport model for coarse grid and unsteady turbulent flow simulations. Theoretical and Computational Fluid Dynamics 18, 443–468; Chaouat, B., Schiestel, R., 2005. A new partially integrated transport model for subgrid-scale stresses and dissipation rate for turbulent developing flows. Physics of Fluids, 17 (6)] The PITM method can be applied to almost all statistical models to derive its hybrid LES counterpart. In the present work, the PITM version based on the transport equations for the turbulent Reynolds stresses together with the dissipation transport rate equation is now developed in a general formulation based on a new accurate energy spectrum function E(κ) valid in both large and small eddy ranges that allows to calibrate more precisely the c_sgs2 function involved in the subgrid dissipation rate _sgs transport equation. The model is also proposed here in an extended form which remains valid in low Reynolds number turbulent flows. This is achieved by considering a characteristic turbulence length-scale based on the total turbulent energy and the total dissipation rate taking into account the subgrid and resolved parts of the dissipation rate. These improvements allow to consider a large range of flows including various free flows as well as bounded flows. The present model is first tested on the decay of homogeneous isotropic turbulence by referring to the well known experiment of Comte-Bellot and Corrsin. Then, initial perturbed spectra E(κ) with a peak or a defect of energy are considered for analysing the model capabilities in strong non-equilibrium flow situations. The second test case is the classical fully turbulent channel flow that allows to assess the performance of the model in non-homogeneous flows characterised by important anisotropy effects. Different simulations are performed on coarse and refined meshes for checking the grid independence of solutions as well as the consistency of the subgrid-scale model when the filter width is changed. A special attention is devoted to the sharing out of the energy between the subgrid-scales and the resolved scales. Both the mean velocity and the turbulent stress computations are compared with data from direct numerical simulations.     

Density-Dependent Incompressible Fluids in Bounded Domains

This paper is devoted to the study of the initial value problem for density dependent incompressible viscous fluids in a bounded domain of with boundary. Homogeneous Dirichlet boundary conditions are prescribed on the velocity. Initial data are almost critical in term of regularity: the initial density is in W^1,q for some q > N, and the initial velocity has fractional derivatives in L^r for some r > N and arbitrarily small. Assuming in addition that the initial density is bounded away from 0, we prove existence and uniqueness on a short time interval. This result is shown to be global in dimension N = 2 regardless of the size of the data, or in dimension N ≥ 3 if the initial velocity is small. Similar qualitative results were obtained earlier in dimension N = 2, 3 by O. Ladyzhenskaya and V. Solonnikov in [18] for initial densities in W^1,∞ and initial velocities in with q > N.     

Effect of lubricant inertia in bearings with a non-Newtonian lubricant

Summary The inertia effects in externally pressurized and squeeze film bearings with lubricants obeying a power law are considered. It is found that the inertia forces decrease the load capacity of the externally pressurized bearing with a given flow rate and the inertia effect increases with the flow behaviour index. At a given feeding pressure, on the other hand, the inertia increases or decreases the load capacity when the flow behaviour index is smaller than or greater than 3, respectively. For squeeze films between circular plates and rectangular plates, the rate of squeeze is slowed down by the inertia and the inertia effect is larger in dilatant lubricants than in pseudoplastic lubricants.Nomenclature 2a diameter of the bearing, width of rectangular plates - 2b diameter of the recess - 2h film thickness - 2h ₀ initial thickness of squeeze films - l length of the rectangular plates - m consistency index - n flow behaviour index - p pressure - p _e external pressure - p _i feeding pressure - q flow rate - r radial distance - t time - u velocity of the lubricant - v squeeze velocity - w load capacity - W dimensionless load capacity - axial distance - viscosity of the lubricant - density of the lubricant     

Cardiac valve prosthesis flow performances measured by 2D and 3D-stereo particle image velocimetry

The present paper reports the design and performance tests of a new artificial heart valve test bench that was specially devised to employ Particle Image Velocimetry (PIV) to perform flow analysis. Among the useful characteristics of this new test bench are the repeatable high quality of the developed flow, the generation of time-variable flow with a feedback-controlled actuator, high versatility in changing and controlling the flow parameters (average rate and the beat frequency), and good optical access for PIV measurements. Different chambers and flow conditions have been used to perform tests using 2D-PIV and 3D-StereoPIV. The aim of such tests were to study the fluid dynamical characteristics of mechanical and biological cardiac prostheses, and to evaluate the procedures used to reduce measurement uncertainty due to the 3D components of the pulsatilee flow through the cardiac valve prostheses. It has been possible to observe the evolution of the complete 2D and 3D flow disturbance induced by the valve prostheses for each phase of their cycle on the upstream and downstream volumes.Abbreviations Re Reynold Number - St Strouhal Number - Q volumetric flow rate - density - dynamic viscosity - pulsatile flow parameter - f pulse rate - d nominal pipe diameter - d _i internal pipe diameter - d _e external pipe diameter - s pipe thickness - t -time - v -velocity - Scheimpflug condition angle