A marker-and-cell approach to free surface 2-D multiphase flows

This work describes a methodology to simulate free surface incompressible multiphase flows. This novel methodology allows the simulation of multiphase flows with an arbitrary number of phases, each of them having different densities and viscosities. Surface and interfacial tension effects are also included. The numerical technique is based on the GENSMAC front-tracking method. The velocity field is computed using a finite-difference discretization of a modification of the NavierStokes equations. These equations together with the continuity equation are solved for the two-dimensional multiphase flows, with different densities and viscosities in the different phases. The governing equations are solved on a regular Eulerian grid, and a Lagrangian mesh is employed to track free surfaces and interfaces. The method is validated by comparing numerical with analytic results for a number of simple problems; it was also employed to simulate complex problems for which no analytic solutions are available. The method presented in this paper has been shown to be robust and computationally efficient. Copyright (c) 2012 John Wiley & Sons, Ltd.

Experimental Comparison of Two Degrees-of-Freedom Vortex-Induced Vibration on High and Low Aspect Ratio Cylinders with Small Mass Ratio

Vortex-induced motion (VIM) is a specific way for naming the vortex-induced vibration (VIV) acting on floating units. The VIM phenomenon can occur in monocolumn production, storage and offloading system (MPSO) and spar platforms, structures presenting aspect ratio lower than 4 and unity mass ratio, i.e., structural mass equal to the displaced fluid mass. These platforms can experience motion amplitudes of approximately their characteristic diameters, and therefore, the fatigue life of mooring lines and risers can be greatly affected. Two degrees-of-freedom VIV model tests based on cylinders with low aspect ratio and small mass ratio have been carried out at the recirculating water channel facility available at NDF-EPUSP in order to better understand this hydro-elastic phenomenon. The tests have considered three circular cylinders of mass ratio equal to one and different aspect ratios, respectively L/D = 1.0, 1.7, and 2.0, as well as a fourth cylinder of mass ratio equal to 2.62 and aspect ratio of 2.0. The Reynolds number covered the range from 10 000 to 50 000, corresponding to reduced velocities from 1 to approximately 12. The results of amplitude and frequency in the transverse and in-line directions were analyzed by means of the Hilbert-Huang transform method (HHT) and then compared to those obtained from works found in the literature. The comparisons have shown similar maxima amplitudes for all aspect ratios and small mass ratio, featuring a decrease as the aspect ratio decreases. Moreover, some changes in the Strouhal number have been indirectly observed as a consequence of the decrease in the aspect ratio. In conclusion, it is shown that comparing results of small-scale platforms with those from bare cylinders, all of them presenting low aspect ratio and small mass ratio, the laboratory experiments may well be used in practical investigation, including those concerning the VIM phenomenon acting on platforms. [DOI: 10.1115/1.4006755]

Two-degree-of-freedom vortex-induced vibration of circular cylinders with very low aspect ratio and small mass ratio

The investigation of vortex-induced vibration on very short cylinders with two degrees of freedom has drawn the attention of a large number of researchers. Some investigations on such a problem are carried out in order to have a better understanding of the physics involved in vortex-induced motions of floating bodies such as offshore platforms. In this paper, experiments were carried out in a recirculating water channel over the range of Reynolds number 6000number decreases when the aspect ratio decreases, as also observed in previous works available in the literature, most of them for stationary cylinders.

Computational simulation of the flow-induced vibration of a circular cylinder subjected to wake interference

In this work, we considered the flow around two circular cylinders of equal diameter placed in tandem with respect to the incident uniform flow. The upstream cylinder was fixed and the downstream cylinder was completely free to move in the cross-stream direction, with no spring or damper attached to it. The centre-to-centre distance between the cylinders was four diameters, and the Reynolds number was varied from 100 to 645. We performed two- and three-dimensional simulations of this flow using a Spectral/hp element method to discretise the flow equations, coupled to a simple Newmark integration routine that solves the equation of the dynamics of the cylinder. The differences of the behaviours observed in the two- and three-dimensional simulations are highlighted and the data is analysed under the light of previously published experimental results obtained for higher Reynolds numbers.

Two-degree-of-freedom vortex-induced vibration of circular cylinders with very low aspect ratio and small mass ratio

The investigation of vortex-induced vibration on very short cylinders with two degrees of freedom has drawn the attention of a large number of researchers. Some investigations on such a problem are carried out in order to have a better understanding of the physics involved in vortex-induced motions of floating bodies such as offshore platforms. In this paper, experiments were carried out in a recirculating water channel over the range of Reynolds number 6000number decreases when the aspect ratio decreases, as also observed in previous works available in the literature, most of them for stationary cylinders.

Experimental analysis of coaxial jets: instability, flow and mixing characterization

The velocity and mixing field of two turbulent jets configurations have been experimentally characterized by means of cold- and hot-wire anemometry in order to investigate the effects of the initial conditions on the flow development. In particular, experiments have been focused on the effect of the separation wall between the two streams on the flow field. The results of the experiments have pointed out that the wake behind a thick wall separating wall has a strong influence on the flow field evolution. For instance, for nearly unitary velocity ratios, a clear vortex shedding from the wall is observable. This phenomenon enhances the mixing between the inner and outer shear layer. This enhancement in the fluctuating activity is a consequence of a local absolute instability of the flow which, for a small range of velocity ratios, behaves as an hydrodynamic oscillator with no sensibility to external perturbations. It has been suggested indeed that this absolute instability can be used as a passive method to control the flow evolution. Finally, acoustic excitation has been applied to the near field in order to verify whether or not the observed vortex shedding behind the separating wall is due to a global oscillating mode as predicted by the theory. A new scaling relationship has been also proposed to determine the preferred frequency for nearly unitary velocity ratios. The proposed law takes into account both the Reynolds number and the velocity ratio dependence of this frequency and, therefore, improves all the previously proposed relationships.

A LES study of a modified Ahmed body geometry

A numerical study using Large Eddy Simulation Coherent Structure Model (LES-CSM), of the flow around a simplified Ahmed body, has been done in this work of thesis. The models used are two salient geometries from the experimental investigation performed in [1], and consist, in particular, in two notch-back body geometries. Six simulation are carried out in total, changing Reynolds number and back-light angle of the model’s rear part. The Reynolds numbers used, based on the height of the models and the free stream velocity, are Re = 10000, Re = 30000 and Re = 50000. The back-light angles of the slanted surface with respect to the horizontal roof surface, that characterizes the vehicle, are taken as B = 31.8◦ and B = 42◦ respectively. The experimental results in [1] have shown that, depending on the parameter B, asymmetric and symmetric averaged flow over the back-light and in the wake for a symmetric geometry can be observed. The aims of the present work of master thesis are principally two. The first aim is to investigate and confirm the influence of the parameter B on the presence of the asymmetry of the averaged flow, and confirm the features described in the experimental results. The second important aspect is to investigate and observe the influence of the second variable, the Reynolds number, in the developing of the asymmetric flow itself. The results have shown the presence of the mentioned asymmetry as well as an influence of the Reynolds number on it.

Experimental analysis of a separation control device based on plasma technology

In this thesis effects of plasma actuators based on Dielectric Barrier Discharge (DBD) technology over a NACA 0015 bidimensional airfoil have been analyzed in an experimental way, at low Reynolds number. Work developed on thesis has been carried on in partnership with the Department of Electrical Engineering of Università di Bologna, inside Wind Tunnel of the Applied Aerodynamic Laboratory of Aerospace Engineering faculty. In order to verify the effectiveness of these active control devices, the analysis has shown how actuators succeed in prevent boundary layer separation only in certain conditions af angle of attack and Reynolds numbers. Moreover, in this thesis actuators’ chordwise position effect has been also analyzed, together with the influence of steady and unsteady operations.

Experimental and Numerical Analysis of Gas Forced Convection through Microtubes and Micro Heat Exchangers

The last decade has witnessed very fast development in microfabrication technologies. The increasing industrial applications of microfluidic systems call for more intensive and systematic knowledge on this newly emerging field. Especially for gaseous flow and heat transfer at microscale, the applicability of conventional theories developed at macro scale is not yet completely validated; this is mainly due to scarce experimental data available in literature for gas flows. The objective of this thesis is to investigate these unclear elements by analyzing forced convection for gaseous flows through microtubes and micro heat exchangers. Experimental tests have been performed with microtubes having various inner diameters, namely 750 m, 510 m and 170 m, over a wide range of Reynolds number covering the laminar region, the transitional zone and also the onset region of the turbulent regime. The results show that conventional theory is able to predict the flow friction factor when flow compressibility does not appear and the effect of fluid temperature-dependent properties is insignificant. A double-layered microchannel heat exchanger has been designed in order to study experimentally the efficiency of a gas-to-gas micro heat exchanger. This microdevice contains 133 parallel microchannels machined into polished PEEK plates for both the hot side and the cold side. The microchannels are 200 µm high, 200 µm wide and 39.8 mm long. The design of the micro device has been made in order to be able to test different materials as partition foil with flexible thickness. Experimental tests have been carried out for five different partition foils, with various mass flow rates and flow configurations. The experimental results indicate that the thermal performance of the countercurrent and cross flow micro heat exchanger can be strongly influenced by axial conduction in the partition foil separating the hot gas flow and cold gas flow.

Testing of subgrid scale (SGS) models for large-eddy simulation (LES) of turbulent channel flow

Sub-grid scale (SGS) models are required in order to model the influence of the unresolved small scales on the resolved scales in large-eddy simulations (LES), the flow at the smallest scales of turbulence. In the following work two SGS models are presented and deeply analyzed in terms of accuracy through several LESs with different spatial resolutions, i.e. grid spacings. The first part of this thesis focuses on the basic theory of turbulence, the governing equations of fluid dynamics and their adaptation to LES. Furthermore, two important SGS models are presented: one is the Dynamic eddy-viscosity model (DEVM), developed by \cite{germano1991dynamic}, while the other is the Explicit Algebraic SGS model (EASSM), by \cite{marstorp2009explicit}. In addition, some details about the implementation of the EASSM in a Pseudo-Spectral Navier-Stokes code \cite{chevalier2007simson} are presented. The performance of the two aforementioned models will be investigated in the following chapters, by means of LES of a channel flow, with friction Reynolds numbers $Re_\tau=590$ up to $Re_\tau=5200$, with relatively coarse resolutions. Data from each simulation will be compared to baseline DNS data. Results have shown that, in contrast to the DEVM, the EASSM has promising potentials for flow predictions at high friction Reynolds numbers: the higher the friction Reynolds number is the better the EASSM will behave and the worse the performances of the DEVM will be. The better performance of the EASSM is contributed to the ability to capture flow anisotropy at the small scales through a correct formulation for the SGS stresses. Moreover, a considerable reduction in the required computational resources can be achieved using the EASSM compared to DEVM. Therefore, the EASSM combines accuracy and computational efficiency, implying that it has a clear potential for industrial CFD usage.

Biaxial extensional motion of an inertially driven radially expanding liquid sheet

We consider the inertially driven, time-dependent biaxial extensional motion of inviscid and viscous thinning liquid sheets. We present an analytic solution describing the base flow and examine its linear stability to varicose (symmetric) perturbations within the framework of a long-wave model where transient growth and long-time asymptotic stability are considered. The stability of the system is characterized in terms of the perturbation wavenumber, Weber number, and Reynolds number. We find that the isotropic nature of the base flow yields stability results that are identical for axisymmetric and general two-dimensional perturbations. Transient growth of short-wave perturbations at early to moderate times can have significant and lasting influence on the long-time sheet thickness. For finite Reynolds numbers, a radially expanding sheet is weakly unstable with bounded growth of all perturbations, whereas in the inviscid and Stokes flow limits sheets are unstable to perturbations in the short-wave limit.

Effect Of Confined Impinging Jet Mixing Processing Parameters On Poly (Lactic Acid) Nanoparticle Morphology

Biodegradable nanoparticles are at the forefront of drug delivery research as they provide numerous advantages over traditional drug delivery methods. An important factor affecting the ability of nanoparticles to circulate within the blood stream and interact with cells is their morphology. In this study a novel processing method, confined impinging jet mixing, was used to form poly (lactic acid) nanoparticles through a solvent-diffusion process with Pluronic F-127 being used as a stabilizing agent. This study focused on the effects of Reynolds number (flow rate), surfactant presence in mixing, and polymer concentration on the morphology of poly (lactic acid) nanoparticles. In addition to looking at the parameters affecting poly (lactic acid) morphology, this study attempted to improve nanoparticle isolation and purification methods to increase nanoparticle yield and ensure specific morphologies were not being excluded during isolation and purification. The isolation and purification methods used in this study were centrifugation and a stir cell. This study successfully produced particles having pyramidal and cubic morphologies. Despite successful production of these morphologies the yield of non-spherical particles was very low, additionally great variability existed between redundant trails. Surfactant was determined to be very important for the stabilization of nanoparticles in solution but appears to be unnecessary for the formation of nanoparticles. Isolation and purification methods that produce a high yield of surfactant free particles have still not been perfected and additional testing will be necessary for improvement.¿

Hemodynamic Modeling of the Intrarenal Circulation

Three dimensional, time dependent numerical simulations of healthy and pathological conditions in a model kidney were performed. Blood flow in a kidney is not commonly investigated by computational approach, in contrast for example, to the flow in a heart. The flow in a kidney is characterized by relatively small Reynolds number (100 < Re < 0.01-laminar regime). The presented results give insight into the structure of such flow, which is hard to measure in vivo. The simulations have suggested that venous thrombosis is more likely than arterial thrombosis-higher shear rate observed. The obtained maximum velocity, as a result of the simulations, agrees with the observed in vivo measurements. The time dependent simulations show separation regimes present in the vicinity of the maximum pressure value. The pathological constriction introduced to the arterial geometry leads to the changes in separation structures. The constriction of a single vessel affects flow in the whole kidney. Pathology results in different flow rate values in healthy and affected branches, as well as, different pulsate cycle characteristic for the whole system.

Studies of spatial clustering of inertial particles in turbulence

We present studies of the spatial clustering of inertial particles embedded in turbulent flow. A major part of the thesis is experimental, involving the technique of Phase Doppler Interferometry (PDI). The thesis also includes significant amount of simulation studies and some theoretical considerations. We describe the details of PDI and explain why it is suitable for study of particle clustering in turbulent flow with a strong mean velocity. We introduce the concept of the radial distribution function (RDF) as our chosen way of quantifying inertial particle clustering and present some original works on foundational and practical considerations related to it. These include methods of treating finite sampling size, interpretation of the magnitude of RDF and the possibility of isolating RDF signature of inertial clustering from that of large scale mixing. In experimental work, we used the PDI to observe clustering of water droplets in a turbulent wind tunnel. From that we present, in the form of a published paper, evidence of dynamical similarity (Stokes number similarity) of inertial particle clustering together with other results in qualitative agreement with available theoretical prediction and simulation results. We next show detailed quantitative comparisons of results from our experiments, direct-numerical-simulation (DNS) and theory. Very promising agreement was found for like-sized particles (mono-disperse). Theory is found to be incorrect regarding clustering of different-sized particles and we propose a empirical correction based on the DNS and experimental results. Besides this, we also discovered a few interesting characteristics of inertial clustering. Firstly, through observations, we found an intriguing possibility for modeling the RDF arising from inertial clustering that has only one (sensitive) parameter. We also found that clustering becomes saturated at high Reynolds number.

A class of exact Navier-Stokes solutions for homogeneous flat-plate boundary layers and their linear stability

We introduce a new boundary layer formalism on the basis of which a class of exact solutions to the Navier–Stokes equations is derived. These solutions describe laminar boundary layer flows past a flat plate under the assumption of one homogeneous direction, such as the classical swept Hiemenz boundary layer (SHBL), the asymptotic suction boundary layer (ASBL) and the oblique impingement boundary layer. The linear stability of these new solutions is investigated, uncovering new results for the SHBL and the ASBL. Previously, each of these flows had been described with its own formalism and coordinate system, such that the solutions could not be transformed into each other. Using a new compound formalism, we are able to show that the ASBL is the physical limit of the SHBL with wall suction when the chordwise velocity component vanishes while the homogeneous sweep velocity is maintained. A corresponding non-dimensionalization is proposed, which allows conversion of the new Reynolds number definition to the classical ones. Linear stability analysis for the new class of solutions reveals a compound neutral surface which contains the classical neutral curves of the SHBL and the ASBL. It is shown that the linearly most unstable Görtler–Hämmerlin modes of the SHBL smoothly transform into Tollmien–Schlichting modes as the chordwise velocity vanishes. These results are useful for transition prediction of the attachment-line instability, especially concerning the use of suction to stabilize boundary layers of swept-wing aircraft.