Fluid Dynamics - Math 6750 Vorticity and Potential Flows - Week 8

Total Page:16

File Type:pdf, Size:1020Kb

Fluid Dynamics - Math 6750 Vorticity and Potential Flows - Week 8 Fluid Dynamics - Math 6750 Vorticity and Potential Flows - Week 8 1 Vorticity and potential flow The nonlinear term of the convective derivative can be expressed using a vector identity as 1 u · ru = r juj2 − u × (r × u): 2 We recall that the vorticity field is defined by ! = r × u and that if the vorticity field is identically zero, then the fluid is irrotational. From Calculus, we know that a field is irrotational if and only if it is a potential flow, i.e. ! = 0 , u = rΦ: Definition 1 (Potential flow.). A potential flow is a flow that can be represented as u = rΦ, where Φ is a scalar function called the velocity potential. Given an irrotational vector field u defined on some fluid region D, we may define a potential φ: D −! R as the line integral Z P φ(P ) = u · dx; O where O 2 D is some arbitrary reference point. For simply-connected fluid region D, the potential φ is independent of the path between O and P . Consequently, φ is a well-defined function and we have u = rφ. If u is irrotational and incompressible, then the corresponding velocity potential φ must be harmonic since r · u = r · rφ = ∆φ = 0: Definition 2 (Circulation.). The circulation of u on C (C simple, smooth, oriented, closed contour) is I ΓC = udγ: C Here the integral is understood as a line integral along C. Remark 1. Stokes' theorem applied to the velocity field yields Z ΓC = n · !dS; S where S is the surface enclosed by C. The right-hand side is sometimes referred to as the flux of vorticity through S. Example 1. If u is a potential flow, then ΓC = 0 on any well behaved closed curve. Remark 2. Since r × rΦ = 0 for any gradient field, the vorticity of a potential flow vanishes and irrotational flow are most always refer to as potential flow. 1 Example 2. Consider a rigid body rotation with angular velocity Ω. The flow field is given by u = (0; 0; Ω) × (x; y; z) = (−Ωy; Ωx; 0) and ! = r × u = 2(0; 0; Ω) = 2Ω 6= 0: Example 3. Consider the two-dimensional irrotational vortex with Ω = (0; 0; α=r2). Then αy αx u = − ; ; 0 r2 r2 and one can verify that ! = 0. Example 4. Consider the stagnation point flow u = (αx; −αy; 0): One can verify that ! = 0 and there exists a velocity potential φ such that u = rφ, with α φ(x; y) = (x2 + y2); 2 up to some additive constant. k Example 5. Consider the line vortex flow u = r eθ, with eθ the unit vector in polar coordinates defined by eθ = (− sin θ; cos θ; 0): Using the gradient operator in polar coordinates, the velocity potential φ, if it exists, satisfies 1 k @ φ = 0; @ φ = ;@ φ = 0: r r θ r z Up to some additive constant, we find that φ = kθ. This doesn't contradict the fact we 3 mentioned above since u is defined on R n f0g which is not simply-connected. For closed curves that don't go around the origin, the line integral is zero. But for the unit circle centered at the origin, the line integral is I u · dx = 2πk: C 2 Inviscid vortical flow A barotropic fluid is an example of a compressible fluid. Definition 3. A barotropic fluid is defined by specifying the pressure as a given function of the density, i.e. p = p(ρ). Theorem 1 (Kelvin). Let C(t) be a simple closed material curve in an inviscid fluid with body force −ρrΦ (conservative). If either ρ =constant or the fluid is barotropic, then the circulation is invariant under the flow: D Γ = 0 (1) Dt C(t) 2 Proof. We parameterise C(t) as C(t) = fx(s; t): 0 ≤ s ≤ 1g: Then ! d d Z d Z 1 @x(s; t) Γ(t) = u(x; t) dx = u(x(s; t); t) · ds dt dt C(t) dt 0 @s Z 1 Du @x(s; t) Z 1 @u(x(s; t); t) = · ds + u(x(s; t); t) · ds 0 Dt @x 0 @s Z 1 Du @x(s; t) Z 1 1 @ = · ds + fu(x(s; t); t) · u(x(s; t); t)g ds: 0 Dt @x 0 2 @s The second integral vanishes since C(t) is a closed curve. Rewriting the material derivative of u using Euler equation, we see that d Z 1 p @x(s; t) Γ(t) = −∇ + χ · ds dt 0 ρ @s Z p = −∇ + χ · dx = 0; C(t) ρ since we are integrating a conservative vector field over a closed curve. Remark 3. Kelvin's theorem allows for a direct Lagrangian description of the dynamics of vorticity, even though vorticity is defined as an Eulerian variable. Remark 4. If a body is moving through an inviscid, incompressible fluid, or if a uniform flow of such a fluid passes around a body, then the vorticity far from the body will be zero. Then, from Kelvin's theorem, the vorticity in the fluid will be zero everywhere. Remark 5. The theorem does not require the fluid region to be simply connected. The Euler equation enters the proof in evaluating the line integral around C, so the theorem continues to hold even if the viscous effect happens to be important elsewhere in the flow. Definition 4. A vortex line is a curve which is everywhere tangent to the vorticity !(x; t). The vortex line which pass through some simple closed curve in space are said to form the boundary of a vortex tube. A consequence of the Kelvin's theorem is the Helmholtz's theorem. Theorem 2 (Helmholtz). 1. The strength of a vortex tube does not vary with time. A vortex tube is a collection of vortex line (a line whose tangent is parallel to the voracity) through a closed contour moving through the fluid. 2. Vortex lines move with the fluid. A vortex cannot end in the fluid. It must form a ring, end at a boundary or go to infinity. 3. Fluid elements initially free of vorticity remain fluid of vorticity. 3 Proof. Let us define a vortex surface as a surface such that ! is everywhere tangent to the surface. Suppose that at t = 0, the vortext line is the intersection of two vortex surfaces S1 and S2. On each of these surfaces, !(x; 0) · ni = 0 where ni is normal to Si, i=1,2: We claim that both S1;S2 remain as vortex surfaces as time advances, i.e. as the points on Si move with the fluid, the surface they comprise at any time t is a vortex surface. Choose ∗ any arbitrary closed material curve C1(t) that encloses a surface S1 (t) ⊂ S1(t). It follows from Stokes theorem that Z Z u · dx = ! · n dS = 0: ∗ C1(0) S1 (0) ∗ This quantity remains zero by virtue of Kelvin's circulation theorem. Since C1(t) and so S1 (t) was arbitrary, this shows that ! · n1 = 0 on S1(t) and an identical argument shows that ! · n2 = 0 on S2(t). We proceed to prove the second statement. Consider any two arbitrary cross sections A1;A2 of a vortex tube and consider the vortex tube bounded by these two cross sections. Since u is incompressible, r · ! = r · (r · u) = 0 and it follows from the divergence theorem that Z Z Z 0 = r · ! dx = ! · n dS = ! · n dS: V @V A1[A2 Let C(t) be a closed material curve which lie on the wall of the vortex tube, enclosing a closed region A(t). Then Z Z !(x; t) · n(x; t) dS = u(x; t) · dx: A(t) C(t) It follows from Kelvin's circulation theorem that Γ remains constant as time proceeds. Example 6 (Vorticity generation near a wall). Consider a one dimensional shear flow near a wall. The viscous shear stress is given by τx = µ (@xv − @yu) = −µωz: Recall that for a Newtonian fluid and incompressible flow, the vorticity equation has the form (ν = µ/ρ) D ! = ! · ru + νr2!: (2) Dt Remark 6. If the flow is steady, then in dimensionless form the vorticity equation is 1 u · r! = ! · ru + r2!: Re The left-hand side represents the advection of vorticity by u, and the second term on the right- hand side represents the transport of vorticity by diffusion. The first time on the right-hand side is the production term associated with the stretching of vortex lines. 4 We now state some other versions of the vorticity equation for inviscid flows. Theorem 3. For an inviscid fluid of constant density, the vorticity equation is D ! = ! · ru: (3) Dt For a barotropic inviscid fluid, the vorticity equation is D ! Dρ ! = ! · ru + : (4) Dt ρ Dt Proof. Assuming that the body force is a potential force, f = −∇Φ, the conservation of linear momentum equation can be written as Du 1 + rp + rΦ = 0: Dt ρ Taking the curl of the previous equation and using the fact that the curl of a gradient field is zero yield Du r × = 0: Dt Expressing the nonlinear term with the vorticity and using the conservation of mass equation Dρ Dt + ρr · u = 0, we find 1 0 = r × @ u + r × (u · ru) = @ ! + r × rjuj2 − u × ! t t 2 = @t! − ! · ru + u · r! − ur · ! + !r · u D! ! Dρ = − ! · ru − : Dt ρ Dt Here we used the vector identity r × (a × b) = b · ra − a · rb + ar · b − br · a and the facts that the divergence of vorticity is zero as well as the curl of a gradient field.
Recommended publications
  • Future Directions of Computational Fluid Dynamics
    Future Directions of Computational Fluid Dynamics F. D. Witherden∗ and A. Jameson† Stanford University, Stanford, CA, 94305 For the past fifteen years computational fluid dynamics (CFD) has been on a plateau. Due to the inability of current-generation approaches to accurately predict the dynamics of turbulent separated flows, reliable use of CFD has been restricted to a small regionof the operating design space. In this paper we make the case for large eddy simulations as a means of expanding the envelope of CFD. As part of this we outline several key challenges which must be overcome in order to enable its adoption within industry. Specific issues that we will address include the impact of heterogeneous massively parallel computing hardware, the need for new and novel numerical algorithms, and the increasingly complex nature of methods and their respective implementations. I. Introduction Computational fluid dynamics (CFD) is a relatively young discipline, having emerged during the last50 years. During this period advances in CFD have been paced by advances in the available computational hardware, which have enabled its application to progressively more complex engineering and scientific prob- lems. At this point CFD has revolutionized the design process in the aerospace industry, and its use is pervasive in many other fields of engineering ranging from automobiles to ships to wind energy. Itisalso a key tool for scientific investigation of the physics of fluid motion, and in other branches of sciencesuch as astrophysics. Additionally, throughout its history CFD has been an important incubator for the formu- lation and development of numerical algorithms which have been seminal to advances in other branches of computational physics.
    [Show full text]
  • Potential Flow Theory
    2.016 Hydrodynamics Reading #4 2.016 Hydrodynamics Prof. A.H. Techet Potential Flow Theory “When a flow is both frictionless and irrotational, pleasant things happen.” –F.M. White, Fluid Mechanics 4th ed. We can treat external flows around bodies as invicid (i.e. frictionless) and irrotational (i.e. the fluid particles are not rotating). This is because the viscous effects are limited to a thin layer next to the body called the boundary layer. In graduate classes like 2.25, you’ll learn how to solve for the invicid flow and then correct this within the boundary layer by considering viscosity. For now, let’s just learn how to solve for the invicid flow. We can define a potential function,!(x, z,t) , as a continuous function that satisfies the basic laws of fluid mechanics: conservation of mass and momentum, assuming incompressible, inviscid and irrotational flow. There is a vector identity (prove it for yourself!) that states for any scalar, ", " # "$ = 0 By definition, for irrotational flow, r ! " #V = 0 Therefore ! r V = "# ! where ! = !(x, y, z,t) is the velocity potential function. Such that the components of velocity in Cartesian coordinates, as functions of space and time, are ! "! "! "! u = , v = and w = (4.1) dx dy dz version 1.0 updated 9/22/2005 -1- ©2005 A. Techet 2.016 Hydrodynamics Reading #4 Laplace Equation The velocity must still satisfy the conservation of mass equation. We can substitute in the relationship between potential and velocity and arrive at the Laplace Equation, which we will revisit in our discussion on linear waves.
    [Show full text]
  • Chapter 4 Differential Relations for a Fluid Particle
    Inviscid potential flow past an array of cylinders. The mathematics of potential theory, pre- sented in this chapter, is both beautiful and manageable, but results may be unrealistic when there are solid boundaries. See Figure 8.13b for the real (viscous) flow pattern. (Courtesy of Tecquipment Ltd., Nottingham, England) 214 Chapter 4 Differential Relations for a Fluid Particle Motivation. In analyzing fluid motion, we might take one of two paths: (1) seeking an estimate of gross effects (mass flow, induced force, energy change) over a finite re- gion or control volume or (2) seeking the point-by-point details of a flow pattern by analyzing an infinitesimal region of the flow. The former or gross-average viewpoint was the subject of Chap. 3. This chapter treats the second in our trio of techniques for analyzing fluid motion, small-scale, or differential, analysis. That is, we apply our four basic conservation laws to an infinitesimally small control volume or, alternately, to an infinitesimal fluid sys- tem. In either case the results yield the basic differential equations of fluid motion. Ap- propriate boundary conditions are also developed. In their most basic form, these differential equations of motion are quite difficult to solve, and very little is known about their general mathematical properties. However, certain things can be done which have great educational value. First, e.g., as shown in Chap. 5, the equations (even if unsolved) reveal the basic dimensionless parameters which govern fluid motion. Second, as shown in Chap. 6, a great number of useful so- lutions can be found if one makes two simplifying assumptions: (1) steady flow and (2) incompressible flow.
    [Show full text]
  • Potential Flow
    1.0 POTENTIAL FLOW One of the most important applications of potential flow theory is to aerodynamics and marine hydrodynamics. Key assumption. 1. Incompressibility – The density and specific weight are to be taken as constant. 2. Irrotationality – This implies a nonviscous fluid where particles are initially moving without rotation. 3. Steady flow – All properties and flow parameters are independent of time. (a) (b) Fig. 1.1 Examples of complicated immersed flows: (a) flow near a solid boundary; (b) flow around an automobile. In this section we will be concerned with the mathematical description of the motion of fluid elements moving in a flow field. A small fluid element in the shape of a cube which is initially in one position will move to another position during a short time interval as illustrated in Fig.1.1. Fig. 1.2 1.1 Continuity Equation v u = velocity component x direction y v y y y v = velocity component y direction u u x u x x x v x Continuity Equation Flow inwards = Flow outwards u v uy vx u xy v yx x y u v 0 - 2D x y u v w 0 - 3D x y z 1.2 Stream Function, (psi) y B Stream Line B A A u x -v The stream is continuity d vdx udy if (x, y) d dx dy x y u and v x y Integrated the equations dx dy C x y vdx udy C Continuity equation in 2 2 y x 0 or x y xy yx if 0 not continuity Vorticity equation, (rotational flow) 2 1 v u ; angular velocity (rad/s) 2 x y v u x y or substitute with 2 2 x 2 y 2 Irrotational flow, 0 Rotational flow.
    [Show full text]
  • Introduction to Computational Fluid Dynamics by the Finite Volume Method
    Introduction to Computational Fluid Dynamics by the Finite Volume Method Ali Ramezani, Goran Stipcich and Imanol Garcia BCAM - Basque Center for Applied Mathematics April 12–15, 2016 Overview on Computational Fluid Dynamics (CFD) 1. Overview on Computational Fluid Dynamics (CFD) 2 / 110 Overview on Computational Fluid Dynamics (CFD) What is CFD? I Fluids: mainly liquids and gases I The governing equations are known, but not their analytical solution: thus, we approximate it I By CFD we typically denote the set of numerical techniques used for the approximate solution (prevision) of the motion of fluids and the associated phenomena (heat exchange, combustion, fluid-structure interaction . ) I The solution of the governing differential (or integro-differential) equations is approximated by a discretization of space and time I From the continuum we move to the discrete level I The CFD is deeply connected to the improvement of computers in the last decades 3 / 110 Overview on Computational Fluid Dynamics (CFD) Applications of CFD I Any field where the fluid motion plays a relevant role: Industry Physics Medicine Meteorology Architecture Environment 4 / 110 Overview on Computational Fluid Dynamics (CFD) Applications of CFD II Moreover, the research front is particularly active: Basic research on fluid New numerical methods mechanics (e.g. transition to turbulence) Application oriented (e.g. renewable energy, competition . ) 5 / 110 Overview on Computational Fluid Dynamics (CFD) CFD: limits and potential I Method Advantages Disadvantages Experimental 1. More realistic 1. Need for instrumentation 2. Allows “complex” problems 2. Scale effects 3. Difficulty in measurements & perturbations 4. Operational costs Theoretical 1. Simple information 1.
    [Show full text]
  • Eindhoven University of Technology MASTER Ghost Counts of Gas Turbine Meters Araujo, S
    Eindhoven University of Technology MASTER Ghost counts of gas turbine meters Araujo, S. Award date: 2004 Link to publication Disclaimer This document contains a student thesis (bachelor's or master's), as authored by a student at Eindhoven University of Technology. Student theses are made available in the TU/e repository upon obtaining the required degree. The grade received is not published on the document as presented in the repository. The required complexity or quality of research of student theses may vary by program, and the required minimum study period may vary in duration. General rights Copyright and moral rights for the publications made accessible in the public portal are retained by the authors and/or other copyright owners and it is a condition of accessing publications that users recognise and abide by the legal requirements associated with these rights. • Users may download and print one copy of any publication from the public portal for the purpose of private study or research. • You may not further distribute the material or use it for any profit-making activity or commercial gain technische Department of Applied Physics u~iversiteit Fluid Dynamics Laberatory TU e emdhoven I Building: Cascade P.O. Box 513 Eindhoven University of Technology 5600MB Eindhoven Title Ghost counts of gas turbine meters Author S.B.Araujo Report number R-1581-A Date April2002 Masters thesis of the period April 2001 - April2002 Group Vortex Dynamics Advisors prof. dr. A. Hirschberg (TUle) dr. H. Riezebos (Gasunie) Abstract Thrbine meters are often used to measure volume flow through pipes. The company responsible for natural gas transportation in the Netherlands ( Gasunie) has recently observed what is called "ghost counts".
    [Show full text]
  • POTENTIAL FLOW: in Which the Vorticity Is Zero
    Chapter 1 Potential flow: 1.1 General formulation Inviscid and irrotational flows in the limit of high Reynolds number are referred to as ‘potential’ or ‘ideal’ flows. The term ‘inviscid’ refers to flows where viscous forces are small compared to inertial forces, so that the fluid viscosity can be neglected in comparison to fluid inertia. ‘Potential’ or ‘ideal’ flows are a class of inviscid flows in which the vorticity ω, which is the curl of the velocity vector, is zero, i. e. ∂uk ωi = ǫijk = 0 (1.1) ∂xj Since the curl of the velocity is zero, the velocity can be expressed as the gradient of a potential φ, ∂φ ui = (1.2) ∂xi and hence the name ‘potential flow’. Using equation 1.2 for the velocity field, the Navier-Stokes mass and mo- mentum equations can be written in terms of the velocity potential. The mass conservation equation, expressed in terms of the velocity potential, is, 2 ∂ui ∂ φ = 2 = 0 (1.3) ∂xi ∂xi Thus, the mass conservation equation simply states that the velocity potential satisfies the Laplace equation, 2φ = 0. Therefore, for potential flow we can use all the techiques developed∇ for solving the Laplace equation earlier. The momentum conservation equation an inviscid flow is, ∂ui ∂ui 1 ∂p fi + uj = + (1.4) ∂t ∂xj −ρ ∂xi ρ where fi is the body force per unit volume acting on the fluid. The second term on the right side of the above equation can be simplified for an irrotational flow 1 2 CHAPTER 1. POTENTIAL FLOW: in which the vorticity is zero.
    [Show full text]
  • A Hamiltonian Boussinesq Model with Horizontally Sheared Currents
    A Hamiltonian Boussinesq model with horizontally sheared currents Elena Gagarina1, Jaap van der Vegt, Vijaya Ambati, and Onno Bokhove2 Department of Applied Mathematics, University of Twente, Enschede, Netherlands 1Email: [email protected] 2Email: [email protected] Abstract We are interested in the numerical modeling of wave-current interactions around beaches’ surf zones. Any model to predict the onset of wave breaking at the breaker line needs to capture both the nonlinearity of the wave and its dispersion. We have formulated the Hamiltonian dynamics of a new water wave model. This model incorporates both the shallow water model and the potential flow model as limiting systems. The variational model derived by Cotter and Bokhove (2010) is such a model, but the variables used have been difficult to work with. Our new model has a three–dimensional velocity field consisting of the full three– dimensional potential field plus horizontal velocity components, such that the vertical component of vorticity is nonzero. Our aims are to augment the new model locally with bores and to derive a numerical finite element discretization of the new model including the capturing of bores. As a preliminary step, the variational finite element discretization of Miles’ variational principle coupled to an elliptic mesh generator is shown. 1. Introduction The beach surf zone is defined as the region of wave breaking and white capping between the moving shore line and the (generally time-dependent) breaker line. Consider nonlinear waves in the deeper water outside the surf zone approaching the beach, before any significant wave breaking occurs. The start of the surf zone on the offshore side is at the breaker line at which sustained wave breaking begins.
    [Show full text]
  • Volume 22 No
    International Journal of Pure and Applied Mathematics ————————————————————————– Volume 22 No. 3 2005, 327-336 DISTRIBUTION OF SPHEROIDAL FOCAL SINGULARITIES IN STOKES FLOW Panayiotis Vafeas Division of Applied Mathematics and Mechanics Department of Engineering Sciences University of Patras Patras, GR-265 04, GREECE e-mail: [email protected] Abstract: Stokes flow for the steady, non-axisymmetric motion of viscous, incompressible fluids in small Reynolds numbers (creeping flow), around small particles embedded within simply connected and bounded flow domains, is de- scribed by a pair of partial differential equations, which evolve the vector bihar- monic velocity and the scalar harmonic total pressure fields. There exist many representations of the solutions of those kinds of flow, in three-dimensional do- mains, appearing in the form of differential operators acting on harmonic and biharmonic potentials. On the other hand, the development of Stokes theory for two-dimensional flows has the advantage that uses only one potential func- tion (stream function) for the representation of the flow fields, but refers to axisymmetric flows. The effect of a distribution of sources – singularities, on the surface of a spheroidal particle or marginally on the focal segment, to the basic flow fields, is the goal of the present work. In particular, the proper confrontation of the problem is ensured by the introduction of the well-known Havelock’s Theorem for the presence of singularities, which provides us with the necessary integral representations of the velocity and the pressure. More- over, the interrelation of the eigenforms of the Papkovich-Neuber differential representation with those that arise from Stokes theory, in two-dimensional Received: May 13, 2005 c 2005, Academic Publications Ltd.
    [Show full text]
  • Horizontal Circulation and Jumps in Hamiltonian Wave Models
    EGU Journal Logos (RGB) Open Access Open Access Open Access Nonlin. Processes Geophys., 20, 483–500, 2013 Advanceswww.nonlin-processes-geophys.net/20/483/2013/ in Annales Nonlinear Processes doi:10.5194/npg-20-483-2013 Geosciences© Author(s) 2013. CC AttributionGeophysicae 3.0 License. in Geophysics Open Access Open Access Natural Hazards Natural Hazards and Earth System and Earth System Sciences Sciences Horizontal circulation and jumpsDiscussions in Hamiltonian wave models Open Access Open Access Atmospheric Atmospheric E. Gagarina1, J. van der Vegt1, and O. Bokhove1,2 Chemistry Chemistry 1Department of Applied Mathematics, University of Twente, Enschede, the Netherlands and Physics2School of Mathematics, University of Leeds,and Leeds, Physics UK Discussions Open Access CorrespondenceOpen Access to: O. Bokhove ([email protected]) Atmospheric Atmospheric Received: 10 January 2013 – Revised: 16 May 2013 – Accepted: 17 May 2013 – Published: 12 July 2013 Measurement Measurement Techniques Techniques Abstract. We are interested in the modelling of wave-current cal model that can predict the onset of wave breaking at the Discussions interactions around surf zones at beaches. Any model that breaker line will need to capture both the nonlinearity of the Open Access Open Access aims to predict the onset of wave breaking at the breaker line waves and their dispersion. Moreover the model has to in- needs to capture both the nonlinearity of the wave and its dis- clude vorticity effects to simulate wave–current interactions. Biogeosciences Biogeosciences persion. We have therefore formulated the HamiltonianDiscussions dy- Various mathematical models are used to describe water namics of a new water wave model, incorporating both the waves.
    [Show full text]
  • An Efficient Spectral-Projection Method for the Navier–Stokes
    JOURNAL OF COMPUTATIONAL PHYSICS 139, 308–326 (1998) ARTICLE NO. CP975872 An Efficient Spectral-Projection Method for the Navier–Stokes Equations in Cylindrical Geometries I. Axisymmetric Cases J. M. Lopez , and Jie Shen † Department of Mathematics and Earth System Science Center, Pennsylvania State University, † University Park, Pennsylvania 16802 E-mail: shen [email protected] Received June 6, 1996; revised August 7, 1997 An efficient and accurate numerical scheme is presented for the axisymmetric Navier–Stokes equations in primitive variables in a cylinder. The scheme is based on a new spectral-Galerkin approximation for the space variables and a second- order projection scheme for the time variable. The new spectral-projection scheme is implemented to simulate the unsteady incompressible axisymmetric flow with a singular boundary condition which is approximated to within a desired accuracy by using a smooth boundary condition. A sensible comparison is made with a standard second-order (in time and space) finite difference scheme based on a stream function- vorticity formulation and with available experimental data. The numerical results indicate that both schemes produce very reliable results and that despite the singular boundary condition, the spectral-projection scheme is still more accurate (in terms of a fixed number of unknowns) and more efficient (in terms of CPU time required for resolving the flow at a fixed Reynolds number to within a prescribed accuracy) than the finite difference scheme. More importantly, the spectral-projection scheme can be readily extended to three-dimensional nonaxisymmetric cases. c 1998 Academic Press 1. INTRODUCTION The main purpose of this paper and its sequel is to develop and validate an efficient and ac- curate numerical scheme for the Navier–Stokes equations (NSE) in cylindrical geometries.
    [Show full text]
  • Numerical Analysis of Combined Wave Radiation and Diffraction On
    water Article Numerical Analysis of Combined Wave Radiation and Diffraction on a Floating Barge Yajie Li 1 , Bin Xu 1,* , Desheng Zhang 1, Xi Shen 1 and Weibin Zhang 2 1 Research Center of Fluid Machinery Engineering and Technology, Jiangsu University, Zhen Jiang 212013, China; [email protected] (Y.L.); [email protected] (D.Z.); [email protected] (X.S.) 2 Key Laboratory (Fluid Machinery and Engineering Research Base) of Sichuan Province, Xihua University, Sichuan 610039, China; [email protected] * Correspondence: [email protected] Received: 21 October 2019; Accepted: 7 January 2020; Published: 11 January 2020 Abstract: A two-dimensional boundary element method (BEM) based on the potential flow theory is adopted to study the combined wave radiation and diffraction by a single barge. The wave-body interaction problems are simulated using a mixed Euler-Lagrangian scheme, with fully nonlinear boundary conditions. The numerical schemes are verified through comparing with existing results, which show that both the wave runups on the barge and hydrodynamic forces can be calculated with sufficient accuracy. Cases of a single barge subjected to sway motion and regular waves are studied. The real contribution of this study is the outcomes of the spectral analysis conducted for test cases when wave radiation effects are considered in addition to pure wave diffraction. The cases of sway motion with the same frequency as incident wave are simulated first. It is found that sway motion will reduce the overall horizontal force when the frequency is lower than a critical frequency. After that, the higher the frequency, the bigger the horizontal force increasing effect.
    [Show full text]