Types of Potential Flows

Total Page:16

File Type:pdf, Size:1020Kb

Types of Potential Flows Types of Potential Flows Simple Approximations of Fluid Flows – Lesson 4 Solving Potential Flow Problems • We now have all the tools for solving complete flow problems using the potential flow model. The essential steps are as follows: 1. Solve Laplace’s equation for either the velocity potential or stream function (in 2D), subject to appropriate boundary conditions. 2. Compute the velocity components from the derivatives of velocity potential or stream function. 3. Compute the pressure field using Bernoulli’s equation. Note that we must define the constant for Bernoulli’s equation from a known boundary condition (e.g., freestream conditions). • Question: How to solve Laplace’s Equation for flows of interest? ‐ Analytical Solutions: Employ exact, analytical solutions to create a flow field of interest. ‐ Numerical Solutions: Use well-known numerical algorithms to compute solutions to Laplace’s equations on a computer. 2 Exact Potential Flow Solutions • We will now examine some selected exact solutions to Laplace’s equation which we will later combine (using superposition) to create various basic flows. • The exact solutions we will work with are: ‐ Uniform Flow ‐ Source/Sink Flow ‐ Vortex Flow ‐ Doublet Flow • In all cases, the validity of the solutions can be proven by substituting them into Laplace’s equation. • The flow configuration we will assume in all cases is an unbounded domain, with one or more features of interest within the domain (e.g., a solid object, vortex, source, etc.). • We will also confine our attention to 2D solutions. 3 2D Laplace’s Equation for 휙 and 휓 휕2휙 휕2휙 휕휙 휕휓 + = 0 푢 = = − y 휕푥2 휕푦2 휕푥 휕푦 Cartesian 휕2휓 휕2휓 휕휙 휕휓 + = 0 푣 = = 휕푥2 휕푦2 휕푦 휕푥 풓 휽 2 휕 휕휙 1 휕 휙 휕휙 1 휕휓 x 푟 + 2 = 0 푣 = − = − origin 휕푟 휕푟 푟 휕휃 푟 휕푟 푟 휕휃 Cylindrical-Polar 휕 휕휓 1 휕2휓 1 휕휙 휕휓 푟 + = 0 푣휃 = = 휕푟 휕푟 푟 휕휃2 푟 휕휃 휕푥 4 Uniform Flow Equipotential Lines • A solution which describes a uniform flow field of magnitude 푉∞ is given by the velocity potential Uniform Flow 휙 푥, 푦 = 푢∞푥 + 푣∞푦 푽∞ 휓 푥, 푦 = −푢∞푦 + 푣∞푥 2 2 Streamlines 푉∞ = 푢∞ + 푣∞ • In cylindrical-polar coordinates 휙 푟, 휃 = 푉∞푟 cos 휃 휓 푟, 휃 = 푉∞푟 sin 휃 5 Potential Flow image courtesy – www.potentialflow.com Point Source/Sink Flow • Consider a 2D point source/sink flow of prescribed strength 푆ሶ (units m3/s per unit depth). We can think of this flow as an idealized fluid volumetric source/sink emerging from a point and issuing uniformly into the domain. Source Flow Sink Flow 퐒ሶ > ퟎ 퐒ሶ < ퟎ • The exact potential flow solution centered at point 푥0, 푦0 is given by: ሶ 푆ሶ 푦 − 푦 푆 2 2 −1 0 휙 푥, 푦 = ln 푥 − 푥0 + 푦 − 푦0 휓 푥, 푦 = tan 4휋 2휋 푥 − 푥0 6 Point Source/Sink Flow (cont.) • In cylindrical-polar coordinates, for a source/sink located at the origin: 푆ሶ 푆ሶ 휙 푟, 휃 = ln 푟 휓 푟, 휃 = 휃 2휋 2휋 Equipotential Lines Streamlines 7 Potential Flow image courtesy – www.potentialflow.com Superposition Example: The Rankine Oval • We can combine solutions to Laplace’s equation (by adding and subtracting them) to create a new solution. This concept is called superposition. • Using the exact solutions for the uniform and source/sink flows, we can create a flow-field around an oval shaped object by superposition. This object is called a Rankine oval. • The recipe for a Rankine oval solution consists of summing: ‐ Uniform flow = 푢∞ 푥 ‐ Source flow of strength +푆ሶ at 푥0 = −푎, 푦0 = 0 ‐ Sink flow of strength −푆ሶ at 푥0 = +푎, 푦0 = 0 • The resulting solutions for the velocity potential and stream function are: 푆ሶ 푥 + 푎 2 + 푦2 ሶ 푆 −1 푦 −1 푦 휙 푥, 푦 = 푢∞푥 + ln 휓 푥, 푦 = −푢 푦 + tan − tan 4휋 푥 − 푎 2 + 푦2 ∞ 2휋 푥 − 푎 푥 + 푎 8 Rankine Oval Solution (푆ሶ = 20, 푎 = ±2) Wall streamline Streamlines Equipotential Flows • Streamline forms a closed curve, which can be thought of as a “body” of the ellipse; streamlines outside of this curve represent potential flows around this body. 9 Potential Flow image courtesy – www.potentialflow.com Potential Doublet (Dipole) • The doublet is a special flow pattern that emerges if we take a source and sink of equal strength and constrain the product of source/sink strength and distance to be a constant (휅), as the distance approaches zero. The resulting solutions for velocity potential and stream function are: Equipotential Flows Streamlines 푑휃 휃2 휃1 푙 푆ሶ 푙 −푆ሶ 푆ሶ −푆ሶ 휅 = 푆ሶ푙 = 푐표푛푠푡푎푛푡 휅 휙 푟, 휃 = cos 휃 2휋푟 휅 휓 푟, 휃 = − sin 휃 2휋푟 10 Potential Flow image courtesy – www.potentialflow.com Potential Vortex Flow • A potential vortex is a rotational flow around a point. • The velocity is prescribed as a pure tangential velocity which is Potential Vortex inversely proportional to the radial distance from the center 퐶 푣 = 휃 푟 where 퐶 is a constant, known as the vortex strength. • It is customary to use the definition of circulation to define the constant as follows: 2휋 2휋 ퟐ Γ = ර 푉휃 푟푑휃 = ර 퐶푑휃 = 2휋퐶 휞 has units of 풎 /풔 0 0 hence Γ 푣 = 휃 2휋푟 11 Vortex Flow • The exact potential flow solution centered at the origin is given by: Streamlines Equipotential Flows Γ 휙 푟, 휃 = − 휃 2휋 Γ 휓 푟, 휃 = ln 푟 2휋 • In Cartesian form: Γ 푦 휙 푟, 휃 = − tan−1 2휋 푥 • Note the potential vortex formulation does not contradict the Γ 2 2 irrotational nature of the flow. The irrotational flow requires that fluid 휓 푟, 휃 = ln 푥 + 푦 elements do not rotate about their axis, which is satisfied in this case, 2휋 as fluid elements move along circular lines around the vortex center but do not rotate about their axes. 12 Potential Flow image courtesy – www.potentialflow.com Potential Flow Past a Circular Cylinder • Using superposition of a doublet and a uniform flow in the 푥 direction, we can construct a potential flow solution for flow around a circular cylinder of radius R: 휅 휙 푟, 휃 = 푉 푟 + cos 휃 ∞ 2휋푟 푦 휅 휓 푟, 휃 = 푉∞ 푟 − sin 휃 2휋푟 푅 • The velocity components are thus 푉∞ 푥 푅2 푣 푟, 휃 = 푉 1 − cos 휃 푟 ∞ 푟2 휅 Where 푅2 = 푅2 2휋푉 푣 푟, 휃 = −푉 1 + sin 휃 ∞ 휃 ∞ 푟2 13 Potential Flow Past a Circular Cylinder (cont.) • We can now apply Bernoulli’s equation to evaluate the pressure along the cylinder’s surface (푟 = 푅). Doing this yields: 1 푝 = 푝 + 휌푉2 1 − 4 sin2휃 푠 ∞ 2 ∞ • Note that the velocity reaches a maximum at the top/bottom of the cylinder and thus pressure is a minimum at those locations. • If we integrate the pressure force in the 푥-direction (drag) around the cylinder, we find that the total force is zero! ‐ This result is known as D’Alembert’s Paradox. A cylinder in potential flow experiences no drag. ‐ The paradox is resolved when we consider that we are not including viscous effects in our model. 14 Comparison with Real Flow • This figure shows a comparison of the potential flow surface pressure distribution for a circular cylinder with the actual flow observed in nature ‐ The agreement is reasonable only near the leading edge of the cylinder Potential flow Potential flow solutions do NOT predict separation. Flow separation is what causes drag on real objects. 15 Potential Flow Past a Rotating Circular Cylinder • Using superposition of a vortex of strength Γ , a doublet and a uniform flow in the 푥 direction, we can construct potential flow solutions for flow around a rotating circular cylinder of radius R: 푅2 Γ 휙 푟, 휃 = 푉∞ 1 + cos 휃 − 휃 푟2 2휋 푦 푅2 Γ 푟 휓 푟, 휃 = 푉 1 − sin 휃 + ln 푅 ∞ 푟2 2휋 푅 푉∞ 푥 • The velocity components are thus 푅2 푣 푟, 휃 = 푉 1 − cos 휃 푟 ∞ 푟2 푅2 Γ 푣 푟, 휃 = −푉 1 + sin 휃 − 휃 ∞ 푟2 2휋푟 16 Potential Flow Past a Rotating Circular Cylinder (cont.) • The velocity on the surface of the cylinder can be obtained by substituting 푟 = 푅 Γ 푣 푎, 휃 = −2푉 sin 휃 − 휃 ∞ 2휋푅 • Using this expression, we can now analyze the pressure coefficient for the rotating cylinder 2 2 푣 2 2Γ sin 휃 Γ 퐶푝 = 1 − = 1 − 4 sin 휃 + + 푉∞ 휋푅푉∞ 2휋푅푉∞ • Just like the previous case, if we integrate the pressure force in the 푥-direction (drag) around the cylinder, we find that the total force is zero! • However, this time due to the asymmetry induced along the 푦-axis, we get a net force along the positive 푦-direction, 퐿’ (Lift per unit span) ′ Kutta-Joukowski Theorem 퐿 = 휌∞푉∞Γ 17 Potential Flow and Lifting Bodies • Recall the definition of circulation given earlier Γ = − ර 푉 ∙ 푑푟Ԧ 퐶 • Kutta-Joukowski Theorem – The lift force exerted on an arbitrary 2D body in a potential flow is proportional to the circulation. C 퐹 퐿 = 휌푉 Γ 푏 ∞ ‐ Body shape is arbitrary ‐ Applications: airfoils, lifting surfaces, struts. ‐ Thin airfoil theory is based on using potential flow methods to determine the lift of an arbitrary thin airfoil shape. 18 Summary • Potential flow solutions are reasonable approximations to real flows if viscous effects can be ignored and flow upstream is undisturbed (irrotational). • Potential flow equations are linear and relatively easy to solve. • The solution techniques discussed in this lesson can be extended to compressible flows, unsteady flows and flows with body forces • Lift forces can be computed quite accurately, but drag forces will be totally missed, since drag is inherently a viscous phenomenon.
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]
  • 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]
  • 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]
  • 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]
  • Variational Water-Wave Models and Pyramidal Freak Waves
    Variational water-wave models and pyramidal freak waves Floriane Marie Pauline Gidel Submitted in accordance with the requirements for the degree of Doctor of Philosophy The University of Leeds Department of Applied Mathematics June 2018 ii iii Declaration The candidate confirms that the work submitted is her own, except where work which has formed part of jointly authored publications has been included. The contribution of the candidate and the other authors to this work has been explicitly indicated below. The candidate confirms that appropriate credit has been given within the thesis where reference has been made to the work of others. The work in Chapter 2 of the thesis has appeared in publication as follows: [57] F. Gidel, O. Bokhove and A. Kalogirou, Variational modelling of extreme waves through oblique interaction of solitary waves: application to Mach reflection, Nonlinear Processes in Geophysics, 24 (2017), 43-60, doi: 10.5194/npg-24-43-2017. The original idea was given by O. Bokhove as an extension of his previous collaboration with A. Kalogirou [20]. Therefore, sections 2.2.2, 2.4.1 and 2.4.2 follow the method initially introduced in [20]. The other sections are directly attributable to the candidate’s work. In addition, the candidate wrote the publication fully, with proofreading from A. Kalogirou, O.Bokhove, M. Kelmanson, G. Kapsenberg and T. Bunnik. The work in Chapters 3 and 4 of the thesis has appeared in publication as follows: [58] F. Gidel, O.Bokhove and M. Kelmanson, Driven nonlinear potential flow with wave breaking at shallow-water beaches, Proc.
    [Show full text]
  • A Coupled Variational Principle for 2D Interactions Between Water-Waves and a Rigid-Body Containing fluid
    Under consideration for publication in JFM RAPIDS 1 A coupled variational principle for 2D interactions between water-waves and a rigid-body containing fluid Hamid Alemi Ardakani1y 1Department of Mathematics, College of Engineering, Mathematics and Physical Sciences, University of Exeter, Penryn Campus, Penryn, Cornwall TR10 9EZ, UK (Received xx; revised xx; accepted xx) New variational principles are given for the two-dimensional interactions between gravity driven water-waves and a rotating and translating rectangular vessel dynamically coupled to its interior potential flow with uniform vorticity. The complete set of equations of motion for the exterior water waves, and the exact nonlinear hydrodynamic equations of motion for the vessel in the roll/pitch, sway/surge, and heave directions, and also the full set of equations of motion for the vessel's interior fluid, relative to the body coordinate system attached to the rotating-translating vessel, are derived from two Lagrangian functionals. 1. Introduction Luke (1967) presented a variational principle for the classical water-wave problem described by equations ∆Φ := ΦXX + ΦYY = 0 for − H (X) < Y < Γ (X; t) ; 1 Φt + rΦ · rΦ + gY = 0 on Y = Γ (X; t) ; 2 (1.1) ΦY = Γt + ΦX ΓX on Y = Γ (X; t) ; ΦY + ΦX HX = 0 on Y = −H (X) ; where (X; Y ) is the spatial coordinate system, Φ (X; Y; t) is the velocity potential of an irrotational fluid lying between Y = −H (X) and Y = Γ (X; t) with the gravity acceleration g acting in the negative Y direction. In the horizontal direction X, the fluid domain is cut off by a vertical surface Σ which extends from the bottom to the free surface.
    [Show full text]
  • Potential Flow Theory
    Potential Flow Theory “When a flow is both frictionless and irrotational, pleasant things happen.” 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.. 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. Potential Flow Theory The Stream Function Is a clever device which allows us to wipe out the continuity equation and solve the momentum equation directly for the single variable. Continuity equation Potential Flow Theory The most common application is incompressible flow in the xy plane This equation is satisfied identically if a function ψ(x, y) is defined such that and the above equation becomes The Eq. shows that this new function must be defined such that Potential Flow Theory The vorticity, or curl V, is an interesting function Where One important application is inviscid irrotational flow in the xy plane, where ωZ =0 and the above Eq. is reduced to This is the second-order Laplace equation for which many solutions and analytical techniques are known. Also, boundary conditions like Eq. reduce to Potential Flow Theory Geometric Interpretation of Stream Function The fancy mathematics above would serve by itself to make the stream function immortal and always useful to engineers. Even better, though, stream function has a beautiful geometric interpretation: Lines of constant stream function are streamlines of the flow.
    [Show full text]