Chapter 3 FEM on Nonlinear Free-Surface Flow
Total Page:16
File Type:pdf, Size:1020Kb
Numerical Models in Fluid-Structure Interaction Chapter 3 FEM on nonlinear free-surface flow K.J. Bai1 & J.W. Kim2 1Seoul National University, Korea 2American Bureau of Shipping, USA 3.1 Introduction Recently there has been a growing need for the treatment of the nonlinear water- wave problems. In treating the nonlinearity, there are basically two different approaches: One is based on the long-wave approximations or a shallow-water theory that gives an approximate formulation including the effect of the nonlinearity. The other is to treat the Laplace equation defined in the fluid domain with the exact nonlinear free-surface boundary conditions and the exact body-boundary condition if a body is present. The former is to treat a simplified nonlinear formulation defined only in a horizontal free-surface plane, which is one dimension less than the original fluid domain. The free-surface flow problems have been of interest to many naval architects and ocean and coastal engineers for a long time. The generation and evolution of the water waves and their interaction with the man-made structures are the main concerns in the problem. The theoretical investigations on the topics have been usually made in the scope of the potential theory by assuming the fluid is inviscid and incompressible. The most distinctive feature of this problem is the 83 WIT Transactions on State of the Art in Science and Engineering, Vol 18, © 2005 WIT Press www.witpress.com, ISSN 1755-8336 (on-line) doi:10.2495/978-1-85312-837-0/03 Numerical Models in Fluid-Structure Interaction presence of the free surface. This is one of the typical free-boundary problems where the part of the boundary should be obtained as a part of the solution. Due to this complexity the theoretical investigation on the problem has been restricted until recently. In the past, the linearized problem has been mainly treated. Sometimes, the nonlinear effect is considered by treating the second-order problem obtained in a systematic perturbation expansion. Another approach, where the leading nonlinearity is included, is to use approximate theoretical models such as the well known KdV and Boussinesq equations based on the long-wave approximation or shallow-water theory. Recently, the application of the Green–Naghdi equation is also proved to be an efficient method for treating nonlinear free-surface flow problems (Ertekin et al., 1986). In this chapter we will present the finite-element method applied to the nonlinear water-wave problems, i.e., in ship motions, wave resistance, lifting problems, and initial-value problems. The present numerical scheme can be used for the validation of the existing approximate theories, i.e., the KdV and the Boussinesq equations or Green–Naghdi equation as well as for a better prediction for more realistic physical models. In Section 3.2, we will describe the finite-element method applied to a two- dimensional hydrofoil problem. In this numerical scheme we introduced a buffer zone where the forced damping is introduced to reduce the wave elevation in matching the nonlinear numerical solution with the far-field linear solution. However, in three dimensions the nonlinear wave amplitude will reduce as it propagates to infinity. Thus the introduction of the fictitious damping term is not essential. In Section 3.3, we treat an initial-value problem for a physical model of an axi-symmetric container filled with water freely falling into a flat solid surface. In this problem we included the surface tension. In the numerical procedures, the finite-element subdivisions are made such that the free-surface elevation can be described as the multi-valued functions, since the surface tension as well as the gravity is important in the flow. In Section 3.4, the generation of solitons in a shallow-water towing-tank near the critical speed as well as a numerical towing- tank simulation for arbitrary tank conditions is discussed. In Section 3.5 a preliminary result of a sloshing problem by the same finite-element method used in Section 3.4 is discussed. In Section 3.6 a nonlinear diffraction problem is discussed. In Section 3.7, a brief concluding remark is given. 84 WIT Transactions on State of the Art in Science and Engineering, Vol 18, © 2005 WIT Press www.witpress.com, ISSN 1755-8336 (on-line) Numerical Models in Fluid-Structure Interaction 3.2 Nonlinear steady waves due to a two-dimensional hydrofoil An application is described of the localized finite-element method to a steady nonlinear free-surface flow past a submerged two-dimensional hydrofoil at an arbitrary angle of attack. The earlier investigations with the linear free-surface boundary condition in Bai (1978) have shown some disagreement between the computed results and the experimental measurements for the cases of shallow submergence. The computational method of solution is the localized finite- element method based on the classical Hamilton’s principle. In the present study, a notable step is introduced in the matching procedure between the fully nonlinear and the linear sub-domains. Additional numerical results of wave resistance, lift force, and circulation strength are presented. Details can be found in Bai and Han (1994). It should be noted that in the earlier treatment of the linear hydrofoil problems, the velocity-potential formulation was used. However, in the following treatment of the nonlinear problem, a stream-function formulation is used since it has a distinct advantage in the numerical treatment. 3.2.1 Mathematical formulation We consider here a steady uniform flow past a fixed two-dimensional hydrofoil submerged in a fluid. The coordinate system is right-handed and rectangular as shown in Fig. 3.1. The y-axis is directed opposite to the force of gravity, and the x-axis coincides with the undisturbed free surface. The unit normal vector nr is always directed outward from the fluid domain. We neglect surface tension and assume that the fluid is inviscid, incompressible, and that the motions are irrotational. Here we assume the water depth, H, is constant. The steady two-dimensional flow is described by a total stream function Ψ(x, y) = Uy +ψ (x, y) , (3.1) where U is the incoming uniform flow velocity in the upstream and ψ (x, y) is the perturbation stream function. The total stream function Ψ of the incoming uniform flow is set to 0 on the undisturbed free surface and −UH on the bottom. The perturbation stream function ψ must satisfy 85 WIT Transactions on State of the Art in Science and Engineering, Vol 18, © 2005 WIT Press www.witpress.com, ISSN 1755-8336 (on-line) Numerical Models in Fluid-Structure Interaction y S F 0 x S U 0 H c S TE S R1 R2 D S B Figure 3.1: Sketch of boundary configurations. ∇2ψ (x, y) = 0 (3.2) in the fluid domain D. The kinematic and dynamic boundary conditions on the free surface SF are as follows: ψ = −Uη or Uη +ψ = 0 x x 1 on S , y =η(x) , (3.3a,b) Uψ + gη + ()∇ψ ⋅∇ψ = 0 F y 2 where y =η(x) and g are the free-surface elevation and the gravitational acceleration, respectively. The boundary condition on the hydrofoil surface S0 can be Ψ(x, y) = C on S0 , (3.4) where the constant C will be determined as a part of the solution. To determine the constant C, we require at the trailing edge denoted by TE an additional condition, i.e., the Kutta condition, which states that the pressure on the upper and lower surfaces is continuous. Therefore, the tangential velocities on the upper and lower surfaces have the same magnitude at the trailing edge and bounded, that is, 86 WIT Transactions on State of the Art in Science and Engineering, Vol 18, © 2005 WIT Press www.witpress.com, ISSN 1755-8336 (on-line) Numerical Models in Fluid-Structure Interaction ∂ ∂ ()Uy +ψ + = (Uy +ψ − ) at TE, (3.5) ∂n ∂n where ∂ / ∂n is the normal derivative and ψ + and ψ − are the perturbation stream function ψ on the upper and lower surfaces, respectively. The boundary condition on the bottom, SB , i.e., y = −H is ψ = 0 on SB . (3.6) y x x x 1 S 2 3 F 0 x S R1 S U 0 S H J J R2 c TE 2 3 J 1 D 1.1 D N DB D1.2 SB Figure 3.2: Subdomains in nonlinear calculation. As the radiation condition, we require that no disturbances exist far upstream, that is, ∇ψ → 0 as x → −∞ , (3.7) and that the perturbed flow velocity is bounded in the far downstream, i.e., ∇ψ is bounded as x → +∞ . (3.8) We hereby assume that the solution of eqns (3.2~3.8) can be determined uniquely. 3.2.2 Variational principle A mathematical formulation of the nonlinear free-surface flow phenomena of an incompressible ideal fluid can be described by variational principles based on the 87 WIT Transactions on State of the Art in Science and Engineering, Vol 18, © 2005 WIT Press www.witpress.com, ISSN 1755-8336 (on-line) Numerical Models in Fluid-Structure Interaction classical Hamilton’s principle as given in Kim and Bai (1990). Therefore, the formulation given in the previous section can be replaced by an equivalent variational functional formulation as follows. Now the nonlinear problem defined by eqns (3.2–3.8) is replaced by an equivalent Hamilton’s principle for the Lagrangian. The Lagrangian L is defined as L = T – V + W, (3.9) where T, V, and W are the kinetic energy, the potential energy, and the work done by external pressure, respectively.