KELVIN-HELMHOLTZ INSTABILITY BY SPH
|
|
- Elwin Page
- 6 years ago
- Views:
Transcription
1 II International Conference on Particle-based Methods - Fundamentals and Applications PARTICLES 2011 E. Oñate and D.R.J. Owen (Eds) KELVIN-HELMHOLTZ INSTABILITY BY SPH M. S. Shadloo, M. Yildiz Faculty of Engineering and Natural Sciences, Advanced Composites & Polymer Processing Laboratory, Sabanci University, Tuzla, Istanbul, Turkey. mostafa@sabanciuniv.edu, meyildiz@sabanciuniv.edu Key words: Smoothed Particle Hydrodynamics (SPH), Shear flow, Two-Phase flow, Interfacial flows, Surface tension, Kelvin-Helmholtz Instability (KHI). Abstract. In this paper, we have modeled the Kelvin-Helmholtz Instability (KHI) problem of an incompressible two-phase immiscible fluid in a stratified inviscid shear flow with interfacial tension using Smoothed Particle Hydrodynamics (SPH) method. The time dependent evolution of the two-fluid interface over a wide range of Richardson number (Ri) and for three different density ratios is numerically investigated. The simulation results are compared with analytical solutions in the linear regime. It was observed that the SPH method requires a Richardson number lower than unity (i.e., Ri = 0.8) for the onset of KHI, and that the artificial viscosity plays a significant role in obtaining physically correct simulation results that are in agreement with analytical solutions. The numerical algorithm presented in this work can easily handle a two-phase fluid flow with various density ratios. 1 INTRODUCTION Flow instability at the interface between two horizontal parallel streams of different velocities and densities, with the heavier fluid at the bottom, is called the Kelvin-Helmholtz Instability (KHI). The KHI is induced by either velocity shear within a continuous fluid or a sufficiently large velocity difference across the interface of a multiphase fluid. The instability kicks in when the destabilizing effect of shear across the interface overcomes the stabilizing effect of stratification due to gravity and/or surface tension if it exists. The KHI manifests itself as a row of horizontal eddies (in the form of waves) aligned across the interface. These eddies or waves are referred to as main billows. There are several well-known natural situations where the KHI can be observed such as wind blowing over the ocean or water surface, a meteor entering the Earth s atmosphere, the interface between the tails of comets and solar wind, or the interface between a liquid layer and a compressible gas, among others. 1
2 The Kelvin-Helmholtz instability problem was solved first for the ideal case of inviscid and incompressible fluids in 1871 by Lord Kelvin. It has been studied both theoretically [1] and experimentally [2], as well as numerically [3, 4]. Smoothed Particles Hydrodynamics (SPH) is one of the most successful Lagrangian methods which was introduced separately by Gingold and Monaghan and Lucy in 1977 initially to simulate astrophysical problems. Due to its flexibility and capability of modeling complex engineering problems, the SPH method has been extensively used for solving a wide variety of problems [5, 6]. As well, there are some recently reported works which have been conducted to investigate the feasibility and the ability of the SPH method to capture the physics behind the KHI [7, 8]. The aim of this work is to simulate the KHI by using a Weakly Compressible SPH (WCSPH) method, thereby showing the ability of the SPH technique to capture this hydrodynamic instability for various density ratios. The current presentation differs from earlier works in the following aspects: the previously reported relevant works treated the fluid to be compressible whereas here it is treated as an incompressible liquid. Unlike other pertinent SPH papers, this paper systematically examines the effects of the stabilizing parameters (e.g. surface tension (σ) and stratifying body force ( ρg)) individually or concurrently on the occurrence of the KHI for a wide range of Ri numbers and three density ratios. 2 SMOOTHED PARTICLE HYDRODYNAMICS SPH is a Lagrangian particle-based method which has been extensively used over a decade for solving partial differential equations (thermo mechanical balance laws integrated with relevant constitutive equations) for both fluids and solids dynamic problems.in comparison to mesh-dependent approaches, the SPH method significantly facilitates the handling of the problems with free or interfacial surface, whose locations are priori unknown and have to be computed during the solution of the problem since the position of particles determine the final configuration of the flow domain. Therefore, it is a particularly powerful tool to study multiphase fluid flow problems. The term particles merely refers to movable points that represent fluid elements in three dimensional space, and each particle is associated with mass, density, velocity and other relevant hydrodynamics properties. In another way of saying, the continuum under the consideration is discretized into finite number of field points which are known as particles. Field variables carried by each particle are calculated through an interpolation (weighting, smoothing, kernel) function, W(r ij, h) from those of surrounding neighbor particles. Here, r ij is the magnitude of the distance vector ( r ij = r i - r j ) between the particle of interest i (an interpolated particle) and its neighboring particles j (neighbors of the interpolated particle), and h is called smoothing length, which is an important input parameter that controls the number of active neighbor particles within the support domain of smoothing function. Here, r i and r j are the position vectors for particles i and j, respectively. The weighting function should have the following attributes; the normalization, the Dirac-delta 2
3 function, the compactness, and spherically symmetric even function properties. Due to the compactness property, W(r ij, h) = 0 when r ij > h, in computations, a given particle interacts with only its nearest neighbors contained within the radius of support domain h. For an arbitrary function (scalar (A), vectorial ( A), or tensorial (A), the kernel approximation in the form of the integral representation can be written as A( r i ) = A( r i ) A( r j )W(r ij, h)d 3 r j. (1) Space where the angle bracket denotes the kernel approximation. In the SPH method, it is common practice to replace the integral operation with a summation over all the particles within the support domain of the kernel function (i.e. r ij < h ). In what follows, one can write A i = A( r i ) = 1 A j W ij (2) ψ i where we have used a concise notation (i.e. W(r ij, h) = W ij ), and the differential volume element d 3 r j is approximated by the reciprocal of the number density ψ i, which is defined for the particle i as j ψ i = j W ij, (3) By using a differentiable kernel, the spatial derivative of the function A evaluated at the location of particle i can be simply reduced to an ordinary differentiation of the interpolant function W ij with respect to the particle i times the volume and the field variables of neighbor particles. It is due to this analytical approximation that the SPH method does not require a grid. For example, one of the widely used SPH gradient discretization scheme of the field A can be written in the form of i A i = ri A( r i ) = 1 (A j A i ) i W ij. (4) ψ i where ri A( r i ) or in a concise notation i A i denotes the spatial derivative of the field A with respect to i th particle coordinates. Since the kernel function determines the extent of the interaction between the particle of interest and its neighbors, its choice is related to the accuracy, efficiency, and stability of the resulting algorithm. In this work, the quintic kernel based on the fifth-order spline function is adopted [9]. j 3
4 2.1 Governing Equations The governing equations used to solve the fluid problems in this article are the mass and linear momentum balance equations which are expressed in the Lagrangian form and given in the direct notation as dρ dt = ρ v, (5) ρ d v dt = p + ρ g + F σ, (6) where p is the pressure, v is the fluid velocity, and F σ denotes a force per unit volume exerted on the fluid-fluid interface due to surface tension. Here, it should be pointed out that in reality, the surface tension force on the fluid-fluid interface is a surface force (i.e., force per unit area). For the sake of computational efficiency and simplifications, we here have followed a common practice whereby the surface tension force per unit area is included in the momentum balance equation as the force per unit volume. As written in eq.(12),this volume force is proportional to product of the interface gradient (gradient of the color field used to distinguish two fluids) and the surface curvature. This approach is referred to as the continuum surface force (CSF) formulation and is widely used in modeling surface tension force of multiphase fluid flows. Also, it should be noted that due to the inviscid fluid assumption, the linear momentum balance equation does not have a viscous force term. In the WCSPH method, for enforcing the incompressibility condition, the density can be evolved either using eq. (5) in combination with eq. (4) or directly utilizing the number density approach as given in eq. (7). ρ i = m i ψ i (7) As for the balance of linear momentum equation, the acceleration term can be discretized in the general form of a i = 1 m i Fi, (8) where F i is the total force acting on the particle i, which can be written in the discrete form as [ F i = m i g + F σ i j ( p i ψ 2 i ] + p j ) + Π ψj 2 ij i W ij. (9) Here, the first term in the summation sign (p i /ψ 2 i + p j/ψ 2 j ) is the discretized form of the pressure gradient force and Π ij is artificial viscosity. 4
5 The artificial viscosity term is included in the current model in order for circumventing numerical instabilities due to the meshless nature of the SPH method. This term introduces some numerical diffusion into the model, thus preventing non-physical oscillations. There are various forms of artificial viscosity terms in SPH literature. In this work, we have implemented the one suggested by Monaghan in [10] in the form of Π ij = 8µ ij( v ij r ij ) ψ i ψ j (r 2 ij + ǫh2 m) (10) where v ij = v i v j, r ij = r i r j, h m = h i+h j 2, ǫ is introduced to prevent singularity when r ij = 0 and µ ij is harmonic average of µ i and µ j defined by 2 µ ij = 1 µ i + 1, (11) µ j where µ i = 1 8 αh ic i ρ i in which c is sound speed and α is the artificial viscosity constant. In accordance with the continuum surface force (CSF) formulation, the surface tension force term in eq.(9) is formulated as F σ i = σ i κ i n i δ s, (12) where κ i is the interface curvature, and n i is the unit vector normal to interface and δ s is the interface dirac delta function(δ s = C ). The unit vector normal to the fluid-fluid interface n i is where i C i is computed as. n i = ic i i C i. (13) i C i = j C ji ψ ij i W ij. (14) where ψ ij = (ψ i + ψ j )/2 and C ji = C j C i. Finally, the interface curvature κ i is calculated from κ i = j 1 ψ ij n ij i W ij, (15) where n ij = n i n j. 5
6 2.2 Time Integration In order to increment the time-steps in the SPH algorithm, we have used a Velocity Verlet predictor corrector method [10]. This technique is an explicit time integration scheme, and is relatively simple to implement. The time integration scheme starts with the predictor step to compute the intermediate particle positions based on current particles velocities as follow r n+1 2 i = r n i + t 2 vn i. (16) where n is the temporal index. Afterward, new densities for each particles are calculated from eq.(3) and eq.(7) at current intermediate particle positions. Since fluid is assumed weakly compressible, the artificial equation of state of the form p p o = c 2 (ρ ρ o ), (17) is used to compute new pressure values at current particle positions. This state equation enforces the incompressibility condition on the flow such that a small variation in density produces a relatively large change in pressure hence preventing the dilatation of the fluid. Here, c is artificial sound velocity and the subscript o denotes the reference condition. The speed of sound for each particle must be chosen carefully to ensure that the fluid is very closely incompressible. In this work, upon selecting the sound speed twenty times larger than the maximum fluid velocity, the density variation is limited to one percent. Having computed the acceleration of each particle from eq.(8) using the current pressure and velocity values, one can obtain the velocity of new time-step through the time integration as v n+1 i = v n i + a n+1 2 i t. (18) Finally, the position of each particle is corrected at the new time-step based on this velocity. It is to be noted that in the SPH method, the order of particles affects the accuracy of interpolations for gradient and Laplacian computations. Therefore, for the computational stability and accuracy, it is preferable to move the SPH particles in a more orderly fashion, which can be achieved through using the XSPH technique proposed by Monaghan et.al [10]. The XSPH method includes the contribution from neighboring particles, thereby causing a fluid particle to move with an average velocity r n+1 i = r n+1 2 i + t n+1 2 v i, (19) where v n+1 i is the XSPH-averaged new velocity by the correction factor of
7 Figure 1: Configuration of the Kelvin-Helmholtz instability at the initial time step. 3 Definition of The Problem The KHI can occur at the interface between two horizontal parallel streams of different velocities and densities, with the heavier fluid at the bottom. For the simulation of this natural flow phenomenon, we here consider two immiscible fluids that are intervened between two infinite parallel horizontal plate with the height of H ( 0 < y < H) (Fig.1). For simplicity, the x-dimension of the computational domain L (0 < x < L) is chosen to be equal to the domain height H (L = H). The computational domain for the KHI problem is represented by a set of particles created on a Cartesian grid with an equidistant particle spacing. At the beginning of the simulations, the computational domain is halved by a horizontal midline (H 1 = H 2 = H 2 and H = 1m), where each half represents a different fluid region. The number of particles for each fluid region is the same. An initial sinusoidal perturbation is applied to the fluidfluid interface through switching the type of particles near the midline. The wavelength of the initial disturbance is set to be equal to the domain length (λ = L) so that the instability is to be confined into mid section of the model domain. The magnitude of the perturbation is ζ o /H 0.03 where ζ o is the initial amplitude of the applied disturbance. Let U 1 and be the velocity and density of the basic state of the upper layer and U 2 and ρ 2 be those of the bottom layer. Particles of two fluids initially at rest are set into motion in opposite directions with the same magnitude (i.e. U 1 = U 2 = U = 0.5 m ). s To be able to show the effect of density on the KH instability, in simulations we have used several density ratios ρ 2 = 2, 5 and 10 where is set to be = 1000 kg. When all m 3 modeling parameters are active, the surface tension force acts (σ) only on the interface particles in the direction of the unit normal whereas gravity (g) acts in the downward direction on all particles. For determining the number of SPH particles for the particle number independent results, we run several test cases where three different particle arrays, namely, 80 80, and were used. It was observed that the array of particles is sufficient for capturing primary wave as well as obtaining particle- 7
8 independent solutions. 4 RESULTS For unstable situations, the analytical non-dimensional growth rate G e in the linear regime can be written as [11] G e = Im(ω) = 2π ρ 2 + ρ 2 1 Ri. (20) where Ri is the dimensionless Richardson number Ri = + ρ 2 k ρ 2 (U 1 U 2 ) 2 ( g(ρ2 ) + k 2 σ ). (21) For numerical investigation, the numerical growth rate G n is calculated in the form of G n = ˆζ/ζ o 1 t, (22) where ˆζ is the amplitude of the disturbance at time t and t is the dimensionless time. t = t U 2 U 1, (23) H where t is real time and H is domain height. To be able to compare the analytical growth rate in eq.(20) with the numerical one in eq.(22), (because eq.(20) is only valid for the linear regime), ˆζ and t are calculated when the wave amplitude reaches to 10 percentage of the domain height (ˆζ = 0.1H). Having perturbed the fluid-fluid interface at the initial time (t = 0) by a small disturbance, under certain input parameters (i.e., surface tension, gravity, density, etc.), the interface disturbance grows and the flow system becomes unstable. Fig.2 illustrates the growth of the interface disturbance as a function of time in the two-dimensional KHI problem for the density ratios of ρ 2 = 2 and ρ 2 = 10 at Ri = For this simulation, the artificial viscosity coefficient in eq.(10) is set to be α = As a result of the interface disturbance, the heavier fluid starts moving in a positive vertical direction, while the lighter fluid in the opposite direction. Thus, both fluids begin to penetrate into each other. As the time progresses, the height of the instability gets larger, and due to the inertial effect, both fluids tend to gain horizontal velocity opposite to their initial bulk velocities. A small vortex appears and the flow regime is no longer linear. This process results in the formation of the main billow. It should be noted that the linear stability growth rate formula is valid only before this time step. At later times, the characteristic form of the KHI becomes much more obvious. Passing time, the non-linear flow regime results in the formation of a Cat s Eye vortex and the fingering phenomenon for the density ratios of ρ 2 = 2 and ρ 2 = 10, respectively. 8
9 M. S. Shadloo, M. Yildiz (I) (II) Figure 2: Time evolution of the interface in the two-dimensional KHI problem with the density ratios of ρ2 ρ2 ρ1 = 2 and ρ1 = 10, and the Ri = 0.01; (a)t = 1.0, (b)t = 2.0, (c)t = 3.0, (d)t = 4.0; (I) particle distributions and (II) velocity vectors. 9
10 In the simulation with the density ratio of 2, both fluids have relatively close inertial forces, and therefore, the vortex is not advected significantly by fluid streams. Consequently, as the simulation progresses, the flow circulation forms the Cat s Eye shape. On the other hand, due to the fact that there exists a relatively large difference in the inertial forces between the upper and the lower fluid layers for the density ratio of 10, and the heavier fluid at the bottom of the modeling domain has a greater inertial force than the lighter fluid at the top, the flow circulation is advected faster in the flow direction of the heavier fluid whereby it leaves the flow domain through the left side and re-enters it from the right side. Accordingly, the translational motion of the flow circulation along the interface brings about the elongation of the crest of the wave, or the fingering phenomenon. From eq.(20) it is also notable that for a given density ratio, the Ri number is the only parameter that controls the stability of the two fluid system in the KHI phenomena. Towards this end, it is important to determine the critical value for this number, which defines the border between stable and unstable flow regimes. The results of the simulations have shown that in the SPH method, the critical value for the Ri number is approximately 0.8 for all density ratios, which is slightly smaller than the one determined using the linear stability analysis. This difference might be attributed to the artificial viscosity utilized in the SPH method, numerical diffusion and the methodology used to perturb the initial fluid-fluid interface. Fig. 3a shows a plot of the growth rate versus the Ri number in which numerically and analytically computed growth rates are compared. One can see from the figure that the numerically computed growth rate decreases with increasing Ri number and/or density ratio, which is consistent with eq.( 20), and simulation results are in close agreement with those corresponding to analytical solutions. The comparison of results presented in figs.2 and 3 for the density ratios of 2 and 10 for a given Ri number reveals that the density ratio significantly affects the shape of the main billow as well as the growth rate. It is also important to note that with increasing density ratio, the transition from a linear to non-linear regime is delayed to later simulation times. The artificial viscosity is one of the reasons that may cause numerically obtained simulation results to deviate slightly from analytical ones. Fig. 3b illustrates the effect of the artificial viscosity on the time evolution of the interface in the two-dimensional KHI problem for one specific test case, which is chosen as a representative for the whole data. In this specific test case, ρ 2 = 10. As seen from the figure, upon choosing a low artificial viscosity coefficient, the numerical results are in better agreement with those of the linear stability analysis. One can also notice that the growth rate decreases as the utilized artificial viscosity coefficient increases. To have stable numerical simulations, the artificial viscosity coefficient can not be chosen to be too small (as an example, α and for ρ 2 = 2 and for ρ 2 = 10, respectively). Therefore, it should be selected carefully in order to have physically valid numerical results, which can predict the KHI phenomena accurately without loosing numerical stability. 10
11 (a) Figure 3: Growth rates G of KHI in the linear regime for various Richardson numbers and (a) the effect of density ratios (α = 0.001) and (b) the effect of artificial viscosities ( ρ2 = 10). 5 CONCLUSIONS The SPH method has been used to simulate Kelvin-Helmholtz instabilities in twodimensional incompressible two-phase flows under the inviscid assumption through taking the effect of surface tension and body forces into account. Simulations were performed for numerous Richardson numbers, density ratios and artificial viscosities. It was shown that under the influence of certain input parameters (body force, surface tension, and density ratios), a flow instability develops in a two-phase fluid system with an initially disturbed fluid-fluid interface. This instability grows in time and, subsequently, the flow system experiences a transition from a linear to a non-linear regime. It was also illustrated that all simulations results are in reasonably good agreement in terms of growth rate with those corresponding to analytical solutions of the linear stability analysis for KHI. The noted discrepancies between numerical analytical results might be attributed to numerical diffusions, to the inclusion of artificial viscosity and to the nature of the initial disturbance of the fluid-fluid interface. Additionally, it was observed that the growth rate is higher for lower density ratio and Richardson numbers, and reaches to free shear flow limit at Richardson numbers near zero. Based on the linear stability analysis, for Richardson numbers lower than unity (Ri < 1), the flow can be unstable, but for the numerical method used in this paper, the flow instability occurs at (Ri) number values less than roughly 0.8. Numerical results suggest that for a given density ratio, the growth rate of the instability is only controlled by the Richardson number, and independent of the nature of stabilizing forces. Also it is shown that the artificial viscosity plays a significant role in all simulation, therefore, it should be chosen such that it preserves the stability of the numerical method as well as captures all complex physics behind this phenomenon. (b) 11
12 ACKNOWLEDGMENT Funding provided by the European Commission Research Directorate General under Marie Curie International Reintegration Grant program with the Grant Agreement number of (PIRG03-GA ) is gratefully acknowledged. The first author also acknowledges the Yousef Jameel Scholarship. REFERENCES [1] Lamb H. Hydrodynamics. Dover Publications, [2] Lin PY, Hanrattty TJ. Prediction of the initiation of slugs with linear stability theory. International journal of multiphase flow 1986; 12(1): [3] Zhang R, He X, Doolen G, Chen S. Surface tension effects on two-dimensional two-phase Kelvin Helmholtz instabilities. Advances in Water Resources 2001; 24(3-4): [4] Herrmann M. A Eulerian level set/vortex sheet method for two-phase interface dynamics. Journal of Computational Physics 2005; 203(2): [5] Monaghan JJ. Simulating free surface flows with sph. Journal of Computational Physics 1994; 110: [6] Shadloo MS, Zainali A, Sadek SH, Yildiz M. Improved incompressible smoothed particle hydrodynamics method for simulating flow around bluff bodies. Computer Methods in Applied Mechanics and Engineering 2011; 200(9-12): [7] Agertz O, Moore B, Stadel J, Potter D, Miniati F, Read J, Mayer L, Gawryszczak A, Kravtsov A, Nordlund A, et al.. Fundamental differences between SPH and grid methods. Monthly Notices of the Royal Astronomical Society 2007; 380(3): [8] Price DJ. Modelling discontinuities and Kelvin-Helmholtz instabilities in SPH. Arxiv preprint arxiv: ;. [9] Morris JP. Analysis of Smoothed Particle Hydrodynamics with Applications. Monash University, [10] Monaghan JJ. Smoothed particle hydrodynamics. Reports on Progress in Physics 2005; 68: [11] Shadloo MS, Yildiz M. Numerical modeling of kelvinhelmholtz instability using smoothed particle hydrodynamics. International Journal for Numerical Methods in Engineering 2011; doi: /nme
Prototype Instabilities
Prototype Instabilities David Randall Introduction Broadly speaking, a growing atmospheric disturbance can draw its kinetic energy from two possible sources: the kinetic and available potential energies
More informationPHYS 432 Physics of Fluids: Instabilities
PHYS 432 Physics of Fluids: Instabilities 1. Internal gravity waves Background state being perturbed: A stratified fluid in hydrostatic balance. It can be constant density like the ocean or compressible
More information1. Comparison of stability analysis to previous work
. Comparison of stability analysis to previous work The stability problem (6.4) can be understood in the context of previous work. Benjamin (957) and Yih (963) have studied the stability of fluid flowing
More informationNonlinear shape evolution of immiscible two-phase interface
Nonlinear shape evolution of immiscible two-phase interface Francesco Capuano 1,2,*, Gennaro Coppola 1, Luigi de Luca 1 1 Dipartimento di Ingegneria Industriale (DII), Università di Napoli Federico II,
More informationFor example, for values of A x = 0 m /s, f 0 s, and L = 0 km, then E h = 0. and the motion may be influenced by horizontal friction if Corioli
Lecture. Equations of Motion Scaling, Non-dimensional Numbers, Stability and Mixing We have learned how to express the forces per unit mass that cause acceleration in the ocean, except for the tidal forces
More informationKelvin Helmholtz Instability
Kelvin Helmholtz Instability A. Salih Department of Aerospace Engineering Indian Institute of Space Science and Technology, Thiruvananthapuram November 00 One of the most well known instabilities in fluid
More informationPressure corrected SPH for fluid animation
Pressure corrected SPH for fluid animation Kai Bao, Hui Zhang, Lili Zheng and Enhua Wu Analyzed by Po-Ram Kim 2 March 2010 Abstract We present pressure scheme for the SPH for fluid animation In conventional
More informationSmooth Particle Hydrodynamic (SPH) Presented by: Omid Ghasemi Fare Nina Zabihi XU Zhao Miao Zhang Sheng Zhi EGEE 520
Smooth Particle Hydrodynamic (SPH) Presented by: Omid Ghasemi Fare Nina Zabihi XU Zhao Miao Zhang Sheng Zhi EGEE 520 OUTLINE Ø Introduction and Historical Perspective: Ø General Principles: Ø Governing
More informationCHAPTER 19. Fluid Instabilities. In this Chapter we discuss the following instabilities:
CHAPTER 19 Fluid Instabilities In this Chapter we discuss the following instabilities: convective instability (Schwarzschild criterion) interface instabilities (Rayleight Taylor & Kelvin-Helmholtz) gravitational
More informationGodunov methods in GANDALF
Godunov methods in GANDALF Stefan Heigl David Hubber Judith Ngoumou USM, LMU, München 28th October 2015 Why not just stick with SPH? SPH is perfectly adequate in many scenarios but can fail, or at least
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationInvestigation of an implicit solver for the simulation of bubble oscillations using Basilisk
Investigation of an implicit solver for the simulation of bubble oscillations using Basilisk D. Fuster, and S. Popinet Sorbonne Universités, UPMC Univ Paris 6, CNRS, UMR 79 Institut Jean Le Rond d Alembert,
More informationThe Evolution of SPH. J. J. Monaghan Monash University Australia
The Evolution of SPH J. J. Monaghan Monash University Australia planetary disks magnetic fields cosmology radiation star formation lava phase change dam break multiphase waves bodies in water fracture
More informationImprovement of Calculation Stability for Slow Fluid Flow Analysis Using Particle Method *
Materials Transactions, Vol. 58, No. 3 (2017) pp. 479 to 484 2017 Japan Foundry Engineering Society Improvement of Calculation Stability for Slow Fluid Flow Analysis Using Particle Method * Naoya Hirata
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More informationInstabilities due a vortex at a density interface: gravitational and centrifugal effects
Instabilities due a vortex at a density interface: gravitational and centrifugal effects Harish N Dixit and Rama Govindarajan Abstract A vortex placed at an initially straight density interface winds it
More informationFluid Dynamics. Part 2. Massimo Ricotti. University of Maryland. Fluid Dynamics p.1/17
Fluid Dynamics p.1/17 Fluid Dynamics Part 2 Massimo Ricotti ricotti@astro.umd.edu University of Maryland Fluid Dynamics p.2/17 Schemes Based on Flux-conservative Form By their very nature, the fluid equations
More information28. THE GRAVITATIONAL INTERCHANGE MODE, OR g-mode. We now consider the case where the magneto-fluid system is subject to a gravitational force F g
8. THE GRAVITATIONAL INTERCHANGE MODE, OR g-mode We now consider the case where the magneto-fluid system is subject to a gravitational force F g =!g, where g is a constant gravitational acceleration vector.
More informationEFFECT OF LIQUID PHASE COMPRESSIBILITY ON MODELING OF GAS-LIQUID TWO-PHASE FLOWS USING TWO-FLUID MODEL
EFFECT OF LIQUID PHASE COMPRESSIBILITY ON MODELING OF GAS-LIQUID TWO-PHASE FLOWS USING TWO-FLUID MODEL Vahid SHOKRI 1*,Kazem ESMAEILI 2 1,2 Department of Mechanical Engineering, Sari Branch, Islamic Azad
More information9 Fluid Instabilities
9. Stability of a shear flow In many situations, gaseous flows can be subject to fluid instabilities in which small perturbations can rapidly flow, thereby tapping a source of free energy. An important
More informationFluid Mechanics Prof. T.I. Eldho Department of Civil Engineering Indian Institute of Technology, Bombay. Lecture - 17 Laminar and Turbulent flows
Fluid Mechanics Prof. T.I. Eldho Department of Civil Engineering Indian Institute of Technology, Bombay Lecture - 17 Laminar and Turbulent flows Welcome back to the video course on fluid mechanics. In
More informationNumerical Studies of Droplet Deformation and Break-up
ILASS Americas 14th Annual Conference on Liquid Atomization and Spray Systems, Dearborn, MI, May 2001 Numerical Studies of Droplet Deformation and Break-up B. T. Helenbrook Department of Mechanical and
More informationTopics in Fluid Dynamics: Classical physics and recent mathematics
Topics in Fluid Dynamics: Classical physics and recent mathematics Toan T. Nguyen 1,2 Penn State University Graduate Student Seminar @ PSU Jan 18th, 2018 1 Homepage: http://math.psu.edu/nguyen 2 Math blog:
More informationModule 3: "Thin Film Hydrodynamics" Lecture 12: "" The Lecture Contains: Linear Stability Analysis. Some well known instabilities. Objectives_template
The Lecture Contains: Linear Stability Analysis Some well known instabilities file:///e /courses/colloid_interface_science/lecture12/12_1.htm[6/16/2012 1:39:16 PM] Linear Stability Analysis This analysis
More informationChapter 1. Continuum mechanics review. 1.1 Definitions and nomenclature
Chapter 1 Continuum mechanics review We will assume some familiarity with continuum mechanics as discussed in the context of an introductory geodynamics course; a good reference for such problems is Turcotte
More informationFluid Dynamics. Massimo Ricotti. University of Maryland. Fluid Dynamics p.1/14
Fluid Dynamics p.1/14 Fluid Dynamics Massimo Ricotti ricotti@astro.umd.edu University of Maryland Fluid Dynamics p.2/14 The equations of fluid dynamics are coupled PDEs that form an IVP (hyperbolic). Use
More informationFluid Animation. Christopher Batty November 17, 2011
Fluid Animation Christopher Batty November 17, 2011 What distinguishes fluids? What distinguishes fluids? No preferred shape Always flows when force is applied Deforms to fit its container Internal forces
More informationSPH for the Modeling of non-newtonian Fluids with Thermal-Dependent Rheology
SPH for the Modeling of non-newtonian Fluids with Thermal-Dependent Rheology G. Bilotta 1,2 1 Dipartimento di Matematica e Informatica, Università di Catania, Italy 2 Istituto Nazionale Geofisica e Vulcanologia,
More information2. FLUID-FLOW EQUATIONS SPRING 2019
2. FLUID-FLOW EQUATIONS SPRING 2019 2.1 Introduction 2.2 Conservative differential equations 2.3 Non-conservative differential equations 2.4 Non-dimensionalisation Summary Examples 2.1 Introduction Fluid
More informationThe effect of a background shear current on large amplitude internal solitary waves
The effect of a background shear current on large amplitude internal solitary waves Wooyoung Choi Dept. of Mathematical Sciences New Jersey Institute of Technology CAMS Report 0506-4, Fall 005/Spring 006
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationReduction of parasitic currents in the DNS VOF code FS3D
M. Boger a J. Schlottke b C.-D. Munz a B. Weigand b Reduction of parasitic currents in the DNS VOF code FS3D Stuttgart, March 2010 a Institut für Aerodynamik und Gasdynamik, Universität Stuttgart, Pfaffenwaldring
More informationEXTENDED FREE SURFACE FLOW MODEL BASED ON THE LATTICE BOLTZMANN APPROACH
METALLURGY AND FOUNDRY ENGINEERING Vol. 36, 2010, No. 2 Micha³ Szucki*, Józef S. Suchy***, Pawe³ ak*, Janusz Lelito**, Beata Gracz* EXTENDED FREE SURFACE FLOW MODEL BASED ON THE LATTICE BOLTZMANN APPROACH
More information12.1 Viscous potential flow (VPF)
1 Energy equation for irrotational theories of gas-liquid flow:: viscous potential flow (VPF), viscous potential flow with pressure correction (VCVPF), dissipation method (DM) 1.1 Viscous potential flow
More informationRoughness Sub Layers John Finnigan, Roger Shaw, Ned Patton, Ian Harman
Roughness Sub Layers John Finnigan, Roger Shaw, Ned Patton, Ian Harman 1. Characteristics of the Roughness Sub layer With well understood caveats, the time averaged statistics of flow in the atmospheric
More informationRecapitulation: Questions on Chaps. 1 and 2 #A
Recapitulation: Questions on Chaps. 1 and 2 #A Chapter 1. Introduction What is the importance of plasma physics? How are plasmas confined in the laboratory and in nature? Why are plasmas important in astrophysics?
More information2 The incompressible Kelvin-Helmholtz instability
Hydrodynamic Instabilities References Chandrasekhar: Hydrodynamic and Hydromagnetic Instabilities Landau & Lifshitz: Fluid Mechanics Shu: Gas Dynamics 1 Introduction Instabilities are an important aspect
More informationThree Dimensional Models of RR Lyrae Pulsation
Regional Variable Star Conference: Physics & Astronomy Department, Michigan State University: 40 Years of Variable Stars: A Celebration of Contributions by Horace A. Smith ed. K. Kinemuchi (Sunspot, NM:
More informationChapter 1. Governing Equations of GFD. 1.1 Mass continuity
Chapter 1 Governing Equations of GFD The fluid dynamical governing equations consist of an equation for mass continuity, one for the momentum budget, and one or more additional equations to account for
More informationFundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics
Fundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI (after: D.J. ACHESON s Elementary Fluid Dynamics ) bluebox.ippt.pan.pl/
More informationINTERNAL GRAVITY WAVES
INTERNAL GRAVITY WAVES B. R. Sutherland Departments of Physics and of Earth&Atmospheric Sciences University of Alberta Contents Preface List of Tables vii xi 1 Stratified Fluids and Waves 1 1.1 Introduction
More informationComputational Astrophysics
Computational Astrophysics Lecture 1: Introduction to numerical methods Lecture 2:The SPH formulation Lecture 3: Construction of SPH smoothing functions Lecture 4: SPH for general dynamic flow Lecture
More informationA semi-implicit finite volume implementation of the CSF method for treating surface tension in interfacial flows
A semi-implicit finite volume implementation of the CSF method for treating surface tension in interfacial flows M. Raessi, M. Bussmann*, and J. Mostaghimi Department of Mechanical and Industrial Engineering,
More informationA Momentum Exchange-based Immersed Boundary-Lattice. Boltzmann Method for Fluid Structure Interaction
APCOM & ISCM -4 th December, 03, Singapore A Momentum Exchange-based Immersed Boundary-Lattice Boltzmann Method for Fluid Structure Interaction Jianfei Yang,,3, Zhengdao Wang,,3, and *Yuehong Qian,,3,4
More informationChapter 3. Stability theory for zonal flows :formulation
Chapter 3. Stability theory for zonal flows :formulation 3.1 Introduction Although flows in the atmosphere and ocean are never strictly zonal major currents are nearly so and the simplifications springing
More informationmeters, we can re-arrange this expression to give
Turbulence When the Reynolds number becomes sufficiently large, the non-linear term (u ) u in the momentum equation inevitably becomes comparable to other important terms and the flow becomes more complicated.
More informationIntroduction to Marine Hydrodynamics
1896 1920 1987 2006 Introduction to Marine Hydrodynamics (NA235) Department of Naval Architecture and Ocean Engineering School of Naval Architecture, Ocean & Civil Engineering First Assignment The first
More informationStability of stratified two-phase flows in horizontal channels
Stability of stratified two-phase flows in horizontal channels I. Barmak a), A. Gelfgat, H. Vitoshkin, A. Ullmann, and N. Brauner School of Mechanical Engineering, Tel Aviv University, Tel Aviv 69978,
More informationcentrifugal acceleration, whose magnitude is r cos, is zero at the poles and maximum at the equator. This distribution of the centrifugal acceleration
Lecture 10. Equations of Motion Centripetal Acceleration, Gravitation and Gravity The centripetal acceleration of a body located on the Earth's surface at a distance from the center is the force (per unit
More informationInterfacial waves in steady and oscillatory, two-layer Couette flows
Interfacial waves in steady and oscillatory, two-layer Couette flows M. J. McCready Department of Chemical Engineering University of Notre Dame Notre Dame, IN 46556 Page 1 Acknowledgments Students: M.
More informationColloquium FLUID DYNAMICS 2012 Institute of Thermomechanics AS CR, v.v.i., Prague, October 24-26, 2012 p.
Colloquium FLUID DYNAMICS 212 Institute of Thermomechanics AS CR, v.v.i., Prague, October 24-26, 212 p. ON A COMPARISON OF NUMERICAL SIMULATIONS OF ATMOSPHERIC FLOW OVER COMPLEX TERRAIN T. Bodnár, L. Beneš
More information1/27/2010. With this method, all filed variables are separated into. from the basic state: Assumptions 1: : the basic state variables must
Lecture 5: Waves in Atmosphere Perturbation Method With this method, all filed variables are separated into two parts: (a) a basic state part and (b) a deviation from the basic state: Perturbation Method
More informationCHAPTER 5. RUDIMENTS OF HYDRODYNAMIC INSTABILITY
1 Lecture Notes on Fluid Dynamics (1.63J/.1J) by Chiang C. Mei, 00 CHAPTER 5. RUDIMENTS OF HYDRODYNAMIC INSTABILITY References: Drazin: Introduction to Hydrodynamic Stability Chandrasekar: Hydrodynamic
More informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More informationAn Overview of Fluid Animation. Christopher Batty March 11, 2014
An Overview of Fluid Animation Christopher Batty March 11, 2014 What distinguishes fluids? What distinguishes fluids? No preferred shape. Always flows when force is applied. Deforms to fit its container.
More informationContents. Parti Fundamentals. 1. Introduction. 2. The Coriolis Force. Preface Preface of the First Edition
Foreword Preface Preface of the First Edition xiii xv xvii Parti Fundamentals 1. Introduction 1.1 Objective 3 1.2 Importance of Geophysical Fluid Dynamics 4 1.3 Distinguishing Attributes of Geophysical
More informationShear instabilities. Chapter Energetics of shear instabilities
Chapter 7 Shear instabilities In this final Chapter, we continue our study of the stability of fluid flows by looking at another very common source of instability, shear. By definition, shear occurs whenever
More informationDynamics of Transient Liquid Injection:
Dynamics of Transient Liquid Injection: K-H instability, vorticity dynamics, R-T instability, capillary action, and cavitation William A. Sirignano University of California, Irvine -- Round liquid columns
More informationDiscussion panel 1: Enforcing incompressibility in SPH
Discussion panel 1: Enforcing incompressibility in SPH Weakly compressible or not weakly compressible, this is the question. Introduction: Andrea Colagrossi Chairman : Antonio Souto Iglesias, Physical/Mathematical
More informationGeneral introduction to Hydrodynamic Instabilities
KTH ROYAL INSTITUTE OF TECHNOLOGY General introduction to Hydrodynamic Instabilities L. Brandt & J.-Ch. Loiseau KTH Mechanics, November 2015 Luca Brandt Professor at KTH Mechanics Email: luca@mech.kth.se
More informationHeat Transfer Benchmark Problems Verification of Finite Volume Particle (FVP) Method-based Code
PROCEEDING OF 3 RD INTERNATIONAL CONFERENCE ON RESEARCH, IMPLEMENTATION AND EDUCATION OF MATHEMATICS AND SCIENCE YOGYAKARTA, 16 17 MAY 2016 Heat Transfer Benchmark Problems Verification of Finite Volume
More informationKelvin-Helmholtz and Rayleigh-Taylor instabilities in partially ionised prominences
Highlights of Spanish Astrophysics VII, Proceedings of the X Scientific Meeting of the Spanish Astronomical Society held on July 9-13, 2012, in Valencia, Spain. J. C. Guirado, L.M. Lara, V. Quilis, and
More informationDonald Slinn, Murray D. Levine
2 Donald Slinn, Murray D. Levine 2 Department of Civil and Coastal Engineering, University of Florida, Gainesville, Florida College of Oceanic and Atmospheric Sciences, Oregon State University, Corvallis,
More informationJet Stability: A computational survey
Jet Stability Galway 2008-1 Jet Stability: A computational survey Rony Keppens Centre for Plasma-Astrophysics, K.U.Leuven (Belgium) & FOM-Institute for Plasma Physics Rijnhuizen & Astronomical Institute,
More informationarxiv: v3 [astro-ph] 2 Jun 2008
Modelling discontinuities and Kelvin-Helmholtz instabilities in SPH Daniel J. Price arxiv:0709.2772v3 [astro-ph] 2 Jun 2008 Abstract School of Physics, University of Exeter, Exeter EX4 4QL, UK In this
More informationTurbulence models and excitation of solar oscillation modes
Center for Turbulence Research Annual Research Briefs Turbulence models and excitation of solar oscillation modes By L. Jacoutot, A. Wray, A. G. Kosovichev AND N. N. Mansour. Motivation and objectives
More informationTurbulence. 2. Reynolds number is an indicator for turbulence in a fluid stream
Turbulence injection of a water jet into a water tank Reynolds number EF$ 1. There is no clear definition and range of turbulence (multi-scale phenomena) 2. Reynolds number is an indicator for turbulence
More informationModeling, Simulating and Rendering Fluids. Thanks to Ron Fediw et al, Jos Stam, Henrik Jensen, Ryan
Modeling, Simulating and Rendering Fluids Thanks to Ron Fediw et al, Jos Stam, Henrik Jensen, Ryan Applications Mostly Hollywood Shrek Antz Terminator 3 Many others Games Engineering Animating Fluids is
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationChE 385M Surface Phenomena University of Texas at Austin. Marangoni-Driven Finger Formation at a Two Fluid Interface. James Stiehl
ChE 385M Surface Phenomena University of Texas at Austin Marangoni-Driven Finger Formation at a Two Fluid Interface James Stiehl Introduction Marangoni phenomena are driven by gradients in surface tension
More informationLevel Set and Phase Field Methods: Application to Moving Interfaces and Two-Phase Fluid Flows
Level Set and Phase Field Methods: Application to Moving Interfaces and Two-Phase Fluid Flows Abstract Maged Ismail Claremont Graduate University Level Set and Phase Field methods are well-known interface-capturing
More informationFluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition
Fluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition C. Pozrikidis m Springer Contents Preface v 1 Introduction to Kinematics 1 1.1 Fluids and solids 1 1.2 Fluid parcels and flow
More informationSimulation of Liquid Jet Breakup Process by Three-Dimensional Incompressible SPH Method
Seventh International Conference on Computational Fluid Dynamics (ICCFD7), Big Island, Hawaii, July 9-13, 212 ICCFD7-291 Simulation of Liquid Jet Breakup Process by Three-Dimensional Incompressible SPH
More informationSimulating Interfacial Tension of a Falling. Drop in a Moving Mesh Framework
Simulating Interfacial Tension of a Falling Drop in a Moving Mesh Framework Anja R. Paschedag a,, Blair Perot b a TU Berlin, Institute of Chemical Engineering, 10623 Berlin, Germany b University of Massachusetts,
More informationInternal boundary layers in the ocean circulation
Internal boundary layers in the ocean circulation Lecture 9 by Andrew Wells We have so far considered boundary layers adjacent to physical boundaries. However, it is also possible to find boundary layers
More informationFigure 1. adiabatically. The change in pressure experienced by the parcel is. dp = -ρ o gξ
6. Internal waves Consider a continuously stratified fluid with ρ o (z) the vertical density profile. z p' ξ p ρ ρ ο (z) Figure 1. Figure by MIT OpenCourseWare. At a point P raise a parcel of water by
More informationTowards Understanding Simulations of Galaxy Formation. Nigel Mitchell. On the Origin of Cores in Simulated Galaxy Clusters
Towards Understanding Simulations of Galaxy Formation Nigel Mitchell On the Origin of Cores in Simulated Galaxy Clusters Work published in the Monthly Notices of the Royal Astronomy Society Journal, 2009,
More informationEULERIAN DERIVATIONS OF NON-INERTIAL NAVIER-STOKES EQUATIONS
EULERIAN DERIVATIONS OF NON-INERTIAL NAVIER-STOKES EQUATIONS ML Combrinck, LN Dala Flamengro, a div of Armscor SOC Ltd & University of Pretoria, Council of Scientific and Industrial Research & University
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationSmoothed Particle Hydrodynamics (SPH) 4. May 2012
Smoothed Particle Hydrodynamics (SPH) 4. May 2012 Calculating density SPH density estimator Weighted summation over nearby particles: ρ(r) = N neigh b=1 m bw (r r b, h) W weight function with dimension
More informationStructure Formation and Particle Mixing in a Shear Flow Boundary Layer
Structure Formation and Particle Mixing in a Shear Flow Boundary Layer Matthew Palotti palotti@astro.wisc.edu University of Wisconsin Center for Magnetic Self Organization Ellen Zweibel University of Wisconsin
More informationDetailed 3D modelling of mass transfer processes in two phase flows with dynamic interfaces
Detailed 3D modelling of mass transfer processes in two phase flows with dynamic interfaces D. Darmana, N.G. Deen, J.A.M. Kuipers Fundamentals of Chemical Reaction Engineering, Faculty of Science and Technology,
More informationDirect Numerical Simulation of Single Bubble Rising in Viscous Stagnant Liquid
Direct Numerical Simulation of Single Bubble Rising in Viscous Stagnant Liquid Nima. Samkhaniani, Azar. Ajami, Mohammad Hassan. Kayhani, Ali. Sarreshteh Dari Abstract In this paper, direct numerical simulation
More informationParticle-based Fluids
Particle-based Fluids Particle Fluids Spatial Discretization Fluid is discretized using particles 3 Particles = Molecules? Particle approaches: Molecular Dynamics: relates each particle to one molecule
More informationPrinciples of Convection
Principles of Convection Point Conduction & convection are similar both require the presence of a material medium. But convection requires the presence of fluid motion. Heat transfer through the: Solid
More informationConvergence of AMR and SPH simulations I. Hydrodynamical resolution and convergence tests
MNRAS 432, 711 727 (2013) Advance Access publication 2013 April 22 doi:10.1093/mnras/stt509 Convergence of AMR and SPH simulations I. Hydrodynamical resolution and convergence tests D. A. Hubber, 1,2,3
More informationToy stars in one dimension
Mon. Not. R. Astron. Soc. 350, 1449 1456 (004) doi:10.1111/j.1365-966.004.07748.x Toy stars in one dimension J. J. Monaghan 1 and D. J. Price 1 School of Mathematical Sciences, Monash University, Clayton
More informationUniversity of Bristol - Explore Bristol Research. Link to publication record in Explore Bristol Research PDF-document.
Dobra, T., Lawrie, A., & Dalziel, S. B. (2016). Nonlinear Interactions of Two Incident Internal Waves. 1-8. Paper presented at VIIIth International Symposium on Stratified Flows, San Diego, United States.
More informationdynamics of f luids in porous media
dynamics of f luids in porous media Jacob Bear Department of Civil Engineering Technion Israel Institute of Technology, Haifa DOVER PUBLICATIONS, INC. New York Contents Preface xvii CHAPTER 1 Introduction
More informationQ.1 The most abundant gas in the atmosphere among inert gases is (A) Helium (B) Argon (C) Neon (D) Krypton
Q. 1 Q. 9 carry one mark each & Q. 10 Q. 22 carry two marks each. Q.1 The most abundant gas in the atmosphere among inert gases is (A) Helium (B) Argon (C) Neon (D) Krypton Q.2 The pair of variables that
More informationPAPER 331 HYDRODYNAMIC STABILITY
MATHEMATICAL TRIPOS Part III Thursday, 6 May, 016 1:30 pm to 4:30 pm PAPER 331 HYDRODYNAMIC STABILITY Attempt no more than THREE questions. There are FOUR questions in total. The questions carry equal
More informationHydrodynamic Instability. Agenda. Problem definition. Heuristic approach: Intuitive solution. Mathematical approach: Linear stability theory
Environmental Fluid Mechanics II Stratified Flow and Buoyant Mixing Hydrodynamic Instability Dong-Guan Seol INSTITUTE FOR HYDROMECHANICS National Research Center of the Helmholtz Association www.kit.edu
More informationDetailed Outline, M E 521: Foundations of Fluid Mechanics I
Detailed Outline, M E 521: Foundations of Fluid Mechanics I I. Introduction and Review A. Notation 1. Vectors 2. Second-order tensors 3. Volume vs. velocity 4. Del operator B. Chapter 1: Review of Basic
More informationarxiv:cond-mat/ v1 [cond-mat.stat-mech] 11 Dec 2002
arxiv:cond-mat/0212257v1 [cond-mat.stat-mech] 11 Dec 2002 International Journal of Modern Physics B c World Scientific Publishing Company VISCOUS FINGERING IN MISCIBLE, IMMISCIBLE AND REACTIVE FLUIDS PATRICK
More informationDirect numerical simulation of salt sheets and turbulence in a double-diffusive shear layer
Direct numerical simulation of salt sheets and turbulence in a double-diffusive shear layer Satoshi Kimura, William Smyth July 3, 00 Abstract We describe three-dimensional direct numerical simulations
More informationConvective Mass Transfer
Convective Mass Transfer Definition of convective mass transfer: The transport of material between a boundary surface and a moving fluid or between two immiscible moving fluids separated by a mobile interface
More informationGoals of this Chapter
Waves in the Atmosphere and Oceans Restoring Force Conservation of potential temperature in the presence of positive static stability internal gravity waves Conservation of potential vorticity in the presence
More informationViscous non-linear theory of Richtmyer-Meshkov Instability. Abstract
Viscous non-linear theory of Richtmyer-Meshkov Instability Pierre Carles and Stéphane Popinet Laboratoire de Modélisation en Mécanique, Université Pierre et Marie Curie, Case 162, 4 place Jussieu, 75252
More informationThe Evolution of Large-Amplitude Internal Gravity Wavepackets
The Evolution of Large-Amplitude Internal Gravity Wavepackets Sutherland, Bruce R. and Brown, Geoffrey L. University of Alberta Environmental and Industrial Fluid Dynamics Laboratory Edmonton, Alberta,
More informationSchool and Conference on Analytical and Computational Astrophysics November, Angular momentum transport in accretion disks
2292-13 School and Conference on Analytical and Computational Astrophysics 14-25 November, 2011 Angular momentum transport in accretion disks Gianluigi Bodo Osservatorio Astronomico, Torino Italy Angular
More information