On linear instability mechanisms in incompressible open cavity flow

A theoretical study of linear global instability of incompressible flow over a rectangular spanwise-periodic open cavity in an unconfined domain is presented. Comparisons with the limited number of results available in the literature are shown. Subsequently, the parameter space is scanned in a systematic manner, varying Reynolds number, incoming boundary-layer thickness and length-to-depth aspect ratio. This permits documenting the neutral curves and leading eigenmode characteristics of this flow. Correlations constructed using the results obtained collapse all available theoretical data on the three-dimensional instabilities.

Global instability analysis and control of three-dimensional long open cavity flow

The three-dimensional wall-bounded open cavity may be considered as a simplified geometry found in industrial applications such as leading gear or slotted flats on the airplane. Understanding the three-dimensional complex flow structure that surrounds this particular geometry is therefore of major industrial interest. At the light of the remarkable former investigations in this kind of flows, enough evidences suggest that the lateral walls have a great influence on the flow features and hence on their instability modes. Nevertheless, even though there is a large body of literature on cavity flows, most of them are based on the assumption that the flow is two-dimensional and spanwise-periodic. The flow over realistic open cavity should be considered. This thesis presents an investigation of three-dimensional wall-bounded open cavity with geometric ratio 6:2:1. To this aim, three-dimensional Direct Numerical Simulation (DNS) and global linear instability have been performed. Linear instability analysis reveals that the onset of the first instability in this open cavity is around Recr 1080. The three-dimensional shear layer mode with a complex structure is shown to be the most unstable mode. I t is noteworthy that the flow pattern of this high-frequency shear layer mode is similar to the observed unstable oscillations in supercritical unstable case. DNS of the cavity flow carried out at different Reynolds number from steady state until a nonlinear saturated state is obtained. The comparison of time histories of kinetic energy presents a clearly dominant energetic mode which shifts between low-frequency and highfrequency oscillation. A complete flow patterns from subcritical cases to supercritical case has been put in evidence. The flow structure at the supercritical case Re=1100 resembles typical wake-shedding instability oscillations with a lateral motion existed in the subcritical cases. Also, This flow pattern is similar to the observations in experiments. In order to validate the linear instability analysis results, the topology of the composite flow fields reconstructed by linear superposition of a three-dimensional base flow and its leading three-dimensional global eigenmodes has been studied. The instantaneous wall streamlines of those composited flows display distinguish influence region of each eigenmode. Attention has been focused on the leading high-frequency shear layer mode; the composite flow fields have been fully recognized with respect to the downstream wave shedding. The three-dimensional shear layer mode is shown to give rise to a typical wake-shedding instability with a lateral motions occurring downstream which is in good agreement with the experiment results. Moreover, the spanwise-periodic, open cavity with the same length to depth ratio has been also studied. The most unstable linear mode is different from the real three-dimensional cavity flow, because of the existence of the side walls. Structure sensitivity of the unstable global mode is analyzed in the flow control context. The adjoint-based sensitivity analysis has been employed to localized the receptivity region, where the flow is more sensible to momentum forcing and mass injection. Because of the non-normality of the linearized Navier-Stokes equations, the direct and adjoint field has a large spatial separation. The strongest sensitivity region is locate in the upstream lip of the three-dimensional cavity. This numerical finding is in agreement with experimental observations. Finally, a prototype of passive flow control strategy is applied.

Compressible modes in a square lid-driven cavity

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Using the GPU to design complex profile extrusion dies

In the present work the benefits of using graphics processing units (GPU) to aid the design of complex geometry profile extrusion dies, are studied. For that purpose, a3Dfinite volume based code that employs unstructured meshes to solve and couple the continuity, momentum and energy conservation equations governing the fluid flow, together with aconstitutive equation, was used. To evaluate the possibility of reducing the calculation time spent on the numerical calculations, the numerical code was parallelized in the GPU, using asimple programing approach without complex memory manipulations. For verificationpurposes, simulations were performed for three benchmark problems: Poiseuille flow, lid-driven cavity flow and flow around acylinder. Subsequently, the code was used on the design of two real life extrusion dies for the production of a medical catheter and a wood plastic composite decking profile. To evaluate the benefits, the results obtained with the GPU parallelized code were compared, in terms of speedup, with a serial implementation of the same code, that traditionally runs on the central processing unit (CPU). The results obtained show that, even with the simple parallelization approach employed, it was possible to obtain a significant reduction of the computation times.

Simulations of incompressible fluid flows by a least squares finite element method

In this work simulations of incompressible fluid flows have been done by a Least Squares Finite Element Method (LSFEM) using velocity-pressure-vorticity and velocity-pressure-stress formulations, named u-p-ω) and u-p-τ formulations respectively. These formulations are preferred because the resulting equations are partial differential equations of first order, which is convenient for implementation by LSFEM. The main purposes of this work are the numerical computation of laminar, transitional and turbulent fluid flows through the application of large eddy simulation (LES) methodology using the LSFEM. The Navier-Stokes equations in u-p-ω and u-p-τ formulations are filtered and the eddy viscosity model of Smagorinsky is used for modeling the sub-grid-scale stresses. Some benchmark problems are solved for validate the numerical code and the preliminary results are presented and compared with available results from the literature. Copyright © 2005 by ABCM.

Simulação numérica de escoamentos de fluidos pelo método de elementos finitos baseado em volumes de controle em malhas não estruturadas

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Simulação numérica de escoamentos de fluídos pelo método de elementos finitos baseado em volumes de controle com a técnica de passo fracionado

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Numerical study on mixed convection in a horizontal flow past a square porous cavity using UNIFAES scheme

Mixed convection on the flow past a heated length and past a porous cavity located in a horizontal wall bounding a saturated porous medium is numerically simulated. The cavity is heated from below. The steady-state regime is studied for several intensities of the buoyancy effects due to temperature variations. The influences of Péclet and Rayleigh numbers on the flow pattern and the temperature distributions are examined. Local and global Nusselt numbers are reported for the heated surface. The convective-diffusive fluxes at the volume boundaries are represented using the UNIFAES, Unified Finite Approach Exponential-type Scheme, with the Power-Law approximation to reduce the computing time. The conditions established by Rivas for the quadratic order of accuracy of the central differencing to be maintained in irregular grids are shown to be extensible to other quadratic schemes, including UNIFAES, so that accuracy estimates could be obtained.

Influence of water flow rate on shear bond strength of resin composite to Er : YAG cavity preparation

Purpose: To evaluate in vitro the influence of water flow rate on shear bond strength of a resin composite to enamel and dentin after Er:YAG cavity preparation. Methods: Ten bovine incisors were selected and roots removed. Crowns were sectioned in four pieces, resulting in 40 samples that were individually embedded in polyester resin (n=10), and ground to plane the enamel and expose the dentin. The bonding site was delimited and samples were randomly assigned according to cavity preparation: (1) Er:YAG/1.0 mL/minute; (2) Er:YAG/1.5 mL/minute; (3) Er:YAG/2.0 mL/minute and (4) High speed handpiece/bur (control group). Samples were fixed to a metallic device, where composite resin cylinders were prepared. Subsequently, they were stored for 24 hours and subjected to a shear bond strength test (500N at 0.5 mm/minute). Results: Means (MPa) were: enamel: 1: 12.8; 2: 16.8; 3: 17.5; 4: 36.0 and Dentin: 1: 13.6; 2: 18.7; 3: 12.1; 4: 21.3. Data were submitted to ANOVA and Tukey`s test. Adhesion to enamel was more efficient than for dentin. The cavities prepared with conventional bur (control) presented higher statistically significant bond strength values (P<0.05) than for Er:YAG laser for both enamel and dentin. No significant differences were observed between water flow rates employed during enamel ablation. For dentin, the shear bond strength of 2.0 mL/minute water flow rate was lower than for 1.5 mL/minute and 1.0 mL/minute rates. The Er:YAG laser adversely affected shear bond strength of resin composite to both enamel and dentin, regardless of the water flow rate used.

Three-dimensional flow instability in a lid-driven isosceles triangular cavity

Linear three-dimensional modal instability of steady laminar two-dimensional states developing in a lid-driven cavity of isosceles triangular cross-section is investigated theoretically and experimentally for the case in which the equal sides form a rectangular corner. An asymmetric steady two-dimensional motion is driven by the steady motion of one of the equal sides. If the side moves away from the rectangular corner, a stationary three-dimensional instability is found. If the motion is directed towards the corner, the instability is oscillatory. The respective critical Reynolds numbers are identified both theoretically and experimentally. The neutral curves pertinent to the two configurations and the properties of the respective leading eigenmodes are documented and analogies to instabilities in rectangular lid-driven cavities are discussed.

Linear instability analysis of incompressible flow over a cuboid cavity

Direct numerical simulations are performed to analyze the three-dimensional instability of flows over three-dimensional cavities. The flow structures at different Reynolds numbers are investigated by using the spectral-element solver nek5000. As the Reynolds number increasing, the lateral wall effects become more important, the recirculation zone shrinks, the front vortex increases and the flow structure inside of the cavity becomes more complex. Results show that the flow bifurcates from a steady state to an oscillatory regime beyond a value of Reynolds number Re = 1100.

Zonal flow induced by longitudinal librations of a rotating cylindrical cavity

Longitudinal librations represent oscillations about the axis of a rotating axisymmetric fluid filled cavity. An analytical theory is developed for the case of a cylindrical cavity in the limit when the libration frequency is small in comparison with the rotation rate, but large in comparison with the inverse of the spin-up time. It is shown that through the nonlinear advection in the Ekman layers the librations cause the fluid to rotate more slowly. © 2010 Elsevier B.V. All rights reserved.

Cellular Renewal and Improvement of Local Cell Effector Activity in Peritoneal Cavity in Response to Infectious Stimuli

The peritoneal cavity (PerC) is a singular compartment where many cell populations reside and interact. Despite the widely adopted experimental approach of intraperitoneal (i.p.) inoculation, little is known about the behavior of the different cell populations within the PerC. To evaluate the dynamics of peritoneal macrophage (Mempty set) subsets, namely small peritoneal Mempty set (SPM) and large peritoneal Mempty set (LPM), in response to infectious stimuli, C57BL/6 mice were injected i.p. with zymosan or Trypanosoma cruzi. These conditions resulted in the marked modification of the PerC myelo-monocytic compartment characterized by the disappearance of LPM and the accumulation of SPM and monocytes. In parallel, adherent cells isolated from stimulated PerC displayed reduced staining for beta-galactosidase, a biomarker for senescence. Further, the adherent cells showed increased nitric oxide (NO) and higher frequency of IL-12-producing cells in response to subsequent LPS and IFN-gamma stimulation. Among myelo-monocytic cells, SPM rather than LPM or monocytes, appear to be the central effectors of the activated PerC; they display higher phagocytic activity and are the main source of IL-12. Thus, our data provide a first demonstration of the consequences of the dynamics between peritoneal Mempty set subpopulations by showing that substitution of LPM by a robust SPM and monocytes in response to infectious stimuli greatly improves PerC effector activity.

Rheological, mechanical and mucoadhesive properties of thermoresponsive, bioadhesive binary mixtures composed of poloxamer 407 and carbopol 974P designed as platforms for implantable drug delivery systems for use in the oral cavity

This study described the formulation and characterisation of the viscoelastic, mechanical and mucoadhesive properties of thermoresponsive, binary polymeric systems composed of poloxamer (P407) and poly(acrylic acid, C974P) that were designed for use as a drug delivery platform within the oral cavity. Monopolymeric and binary polymeric formulations were prepared containing 10, 15 and 20% (w/w) poloxamer (407) and 0.10-0.25% (w/w) poly(acrylic acid, 934P). The flow theological and viscoelastic properties of the formulations were determined using controlled stress and oscillatory rheometry, respectively, the latter as a function of temperature. The mechanical and mucoadhesive properties (namely the force required to break the bond between the formulation and a pre-hydrated mucin disc) were determined using compression and tensile analysis, respectively. Binary systems composed of 10% (w/w) P407 and C934P were elastoviscous, were easily deformed under stress and did not exhibit mucoadhesion. Formulations containing 15 or 20% (w/w) Pluronic P407 and C934P exhibited a sol-gel temperature T(sol/gel), were viscoelastic and offered high elasticity and resistance to deformation at 37 degrees C. Conversely these formulations were elastoviscous and easily deformed at temperatures below the sol-gel transition temperature. The sol-gel transition temperatures of systems containing 15% (w/w) P407 were unaffected by the presence of C934P; however, increasing the concentration of C934P decreased the T(sol/gel) in formulations containing 20%(w/w) P407. Rheological synergy between P407 and C934P at 37 degrees C was observed and was accredited to secondary interactions between these polymers, in addition to hydrophobic interactions between P407 micelles. Importantly, formulations composed of 20% (w/w) P407 and C934P exhibited pronounced mucoadhesive properties. The ease of administration (below the T(sol/gel)) in conjunction with the viscoelastic (notably high elasticity) and mucoadhesive properties (at body temperature) render the formulations composed of 20% (w/w) P407 and C934P as potentially useful platforms for mucoadhesive, controlled topical drug delivery within the oral cavity. (c) 2009 Published by Elsevier B.V.