An improved low-order boundary element method for breaking surface waves
|
|
- Byron Harris
- 6 years ago
- Views:
Transcription
1 Wave Motion 43 (2006) An improved low-order boundary element method for breaking surface waves N. Drimer, Y. Agnon Department of Civil and Environmental Engineering, Technion, Haifa 32000, Israel Received 21 February 2005; received in revised form 16 August 2005; accepted 27 September 2005 Available online 28 November 2005 Abstract A new compatibility condition, applied to the velocity of the free surface particles, enables the simulation of plunging breakers, up to the creation of a thin et, without numerical smoothing. An improved formulation of a boundary element method for the simulation of surface waves problems is derived and tested through some numerical examples. A numerical wave channel with a wavemaker and an energy absorber is developed and used for a comparison of the non-linear solution with first and second order steady solutions in the frequency domain. Finally, the shoaling of with published experimental results Elsevier B.V. All rights reserved. Keywords: Low-order boundary elements; Plunging breakers; Non-linear waves 1. Introduction Mathematical modeling of surface waves requires the solution of a non-linear boundary value problem. One way to treat the non-linearity is by perturbation series with respect to a small parameter, usually the wave slope. However, the perturbation method is limited to weak non-linearity and cannot be used to simulate highly non-linear flows, such as wave breaking. Another approach is a fully non-linear solution by a time stepping procedure. Several numerical models have been presented in the literature, which simulate highly non-linear free surface flows, most of them are based on the boundary element method (BEM) which has the advantage of reducing the spatial dimension involved in the numerical calculations, by one. Longuet-Higgins and Cokelet [1] presented the first numerical simulation of an overturning wave. Small instabilities, which arose in their calculations, were suppressed by a smoothing procedure. Since then, many different formulations were developed and a variety of applications were presented. Many of the models make use of complex potential or conformal mapping [1 5], hence, cannot be extended to three-dimensional problems. One formulation, which may be extended to three-dimensional problems, and is capable of accommodating a body of arbitrary shape, is the direct BEM formulation used by Jansen [6]. New et al. [7] have studied overturning waves, and their results still serve as a benchmark. A convenient method to simulate an open-ended channel is to add an artificial damping term, presented by Le Mehaute [8], to the dynamic boundary condition, in order to simulate an energy absorber. This artificial term is added only in an extension of the wave channel, hence, does not influence the solution in the physical region. Wave Corresponding author /$ see front matter 2005 Elsevier B.V. All rights reserved. doi: /.wavemoti
2 242 N. Drimer, Y. Agnon / Wave Motion 43 (2006) shoaling and wave breaking has been studied analytically and numerically in many studies, which are used to verify the results of the present model [9 15]. A review of most of the methods is given by Dold [2] and by Banner and Peregrine [16]. Perlin et al. [17] used particle image and tracking velocimetry in their experiments. Overturning et velocity reached 1.7 times the linear wave theory phase velocity and the acceleration was up to 1.1 g. Chen and Kharif [18] studied breaking via the two-dimensional Navier Stokes equation. The fluid motion was found to be irrotational up to the et re-entry, supporting the use of a potential flow model prior to the final stage of breaking. Schaffer et al. [19] used a roller concept in a high order Boussinesq model to incorporate breaking. Grilli et al. [20] developed a fully non-linear three-dimensional model using high order BEM. Two recent reviews are presented by Grilli and Subramanya [21] and by Kim et al. [22]. Donescu and Virgin [23] have developed an Eulerian Lagrangian implicit time integration formulation and applied it to 3D wave body interaction. Jansen [6] discusses the order of element, which may be used in the BEM formulation, and indicates several advantages of the linear element over lower or high order elements. However, a maor disadvantage of the linear element is the non-uniqueness of the normal velocity at the nodal points, which causes numerical instability when high curvature develops at the free surface. Such instability may be suppressed by smoothing but this may suppress physical instability as well. Smoothing procedures were used by Jansen as well as by many others [1,2,4,7]. While the use of high order elements is advanced and provides attractive results, it is best suited for a smooth boundary, and encounters difficulties when considering wave breaking and wave body interaction. The present work pursues an alternative course. The idea is to extend the capabilities of the simpler, linear element formulation. The linear element formulation is quite versatile, and suffers mostly from the difficulty in addressing corner flows, which arise, for example, at a breaking wave (see also Grilli et al. [20] and references therein for an evaluation of the difficulties arising at corners). In this paper we show that the instability associated with high curvature may be removed in a physical way by defining a unique normal velocity and relating it to the element variables by means of compatibility conditions. Following Jansen s model, we developed an improved formulation, which enables the simulation of plunging breakers up to the creation of a very thin et, without the use of any numerical smoothing. The angle between two adacent normals may sometimes reach almost 180 before the solution breaks down. The mathematical and numerical formulations are presented in Sections 2 and 3, respectively. In Sections 4.1 and 4.2 two verification problems of quite different nature are solved, and compared with previous results. Then the model is used to produce some new results. In Section 4.3, we simulate a wave channel with a piston wavemaker and an energy absorber. A comparison of the non-linear solution with first and second order steady solutions in the frequency domain is used for the tuning of the energy absorber, and then for investigating the validity of the frequency domain models. Several researchers utilized BEM models to investigate wave shoaling. Perhaps, the main difficulty in simulating waves shoaling is the absence of an exact mathematical expression for the radiation condition. In some works the breaker was generated by developing an initial solution of a steep sine wave, imposing periodicity conditions on both ends of the solution domain [3,6]. Such situations are useful for studying the local solution near the tip of the breaker. Recently, several works have been published, which modeled the hole shore profile, and simulate the whole shoaling process [5,9,10]. However, an initial solution of a solitary wave was used, which, as shown in the beginning of Section 4.4, caused no reflection problems. The design wave used in the design of offshore and coastal structures is usually of given deep-water height and length (or period). Hence, the simulation of waves of finite length is more practical. In Section 4.4, the numerical wave channel of Section 4.3 is utilized to simulate shoaling of waves of finite length, up to breaking. The results are compared with published experimental results. The model presented here has been recently used by Song and Banner [24 26] to study the onset and strength of wave breaking in deep- and shallow-water. 2. Mathematical formulation A two-dimensional flow in the yz plane is assumed, where y is a horizontal coordinate along the undisturbed free surface and z is a vertical coordinate pointed upwards. For irrotational and incompressible flow, there is a velocity potential Φ(y, z, t) defined by V = Φ (2.1)
3 N. Drimer, Y. Agnon / Wave Motion 43 (2006) Fig. 1. Solution domain definitions sketch. where F is the velocity potential, satisfies Laplace s equation: 2 Φ y Φ z 2 = 0 inω (2.2) The solution domain, Ω, is bounded by boundaries which may be free surface, solid or artificial boundaries, as shown in Fig. 1. A solid boundary may be fixed or moving. The free surface panicles satisfy the kinematic and dynamic boundary conditions: Dx Dt = Φ, on Γ f, (2.3) DΦ = 1 Dt 2 Φ Φ gz, on Γ f, (2.4) where x is the position vector (y, z) of a free surface particle and g is the gravity acceleration. The pressure above the free surface is assumed to be zero. On a solid boundary, an impenetration boundary condition will be satisfied Φ n = V s n, on Γ s (2.5) where n is a unit normal to Γ, positive out of Ω, and V s is the local velocity of the solid boundary. If the solution domain Ω is not completely bounded by the free surface Γ f, and the solid boundary, Γ s, additional boundaries must be defined, on which boundary conditions must be formulated. Those are the artificial boundaries and some of the possible are listed below. A simple and very popular possibility is to impose spatial periodicity, by the implementation of the boundary conditions: Φ(0,z) = Φ(λ, z), Φ n (0,z) = Φ n (λ,z), where λ is a pre-defined wavelength. Another possibility is to apply a linear radiation condition of the form ( ) Cn Φ = 0, (2.8) t on a vertical boundary. C is a given phase velocity and the sign defines a ± y propagation, respectively. The linear radiation condition is correct for a periodic wave with a known phase velocity C; but since for non-linear waves this is usually not the case; the artificial boundary will reflect some energy to the solution domain. A non-reflecting boundary condition, which does not require prior knowledge of a wave property, may be formulated by the use of artificial damping. One way to do it is to add an artificial damping term, vφ, to the right hand side of the dynamic boundary condition, where v is a positive coefficient (Le Mehaute [8]). In order to prevent non-physical effects on the solution, the artificial damping is used only along an extension to the region of interest (see Fig. 1). As mentioned by Le Mehaute [8], in order to minimize the energy reflection by the artificial boundary (which in this case is the vertical line, separating the region of interest and the region where the damping is applied), v must be slowly (2.6) (2.7)
4 244 N. Drimer, Y. Agnon / Wave Motion 43 (2006) increasing, from the artificial boundary along the damping region. One disadvantage of the artificial damping method is the increase of computation time by adding elements. In order to absorb satisfactorily steep waves, the extension length has to be about one wavelength. 3. Numerical formulation The initial boundary value problem is solved by a time stepping procedure, where in each time step a boundary value problem is solved in the instantaneous spatial domain. In an improved method for solving the boundary value problem is developed. The time stepping procedure is described in Section Solution of the boundary value problem The basis for the numerical evaluation of Φ is the familiar boundary element method. The direct formulation for a potential problem is used. According to the direct formulation a boundary integral equation is derived by the application of Green s theorem to the velocity potential Φ and a simple source, over a closed contour Γ enclosing Ω. The resultant integral equation is: αφ(ξ, t) = Γ [Φ(x, t) (ln r) ln r Φ(x, t)]dγ (3.1) where { 2π, ξ in Ω α = (3.2) Γ, ξ on Γ (inner angle of the contour Γ ) and r is the distance between x and ξ. Then, an algebraic formulation is derived, by means of discretization of the boundary Γ to finite elements, along which the unknown functions are approximated by their nodal values multiplying by interpolation functions (see Fig. 2): Φ = Φ 1 ψ 1 + φ 2 ψ 2, (3.3) Fig. 2. Boundary mesh definition sketch.
5 Φ = ( ) Φ ψ ( ) Φ ψ 2, 2 N. Drimer, Y. Agnon / Wave Motion 43 (2006) x = x 1 ψ 1 + x 2 ψ 2, (3.5) where the local nodes 1 and 2 of element are the global nodes, + 1, respectively, and ψ 1, ψ 2 are linear interpolation functions in the element region: ψ 1 = 1 (1 τ) 2 ψ 1 = 1 (3.6) (1 τ) 2 1 τ 1 The arc length dγ is, by definition: ( 2 2 ( 2 ) 2 dγ = (dxdx) 1/2 = y k ψ k) + z k ψ k k=1 k=1 1/2 (3.4) dτ = 1 2 S dτ (3.7) where S is the length of the element, and ( ) denotes derivative by τ. The boundary integral in (3.1) is then approximated as a sum of integrals of the interpolation functions ψ 1 and ψ 2 (3.6) over the boundary elements, where the unknown values of Φ and its gradient are left out of the integrals. Satisfying (3.1) for each of the N boundary nodes, results with a set of N equations with 3N unknowns [ ] N =1 A i φ B ( ) i ( Φ ) ( ) B (+) i ( Φ ) (+) = 0, i = 1, 2,...,N (3.8) where A i = 1 2 S 1 ψ 2 e 1 (ln r)dτ S ψ 1 e (ln r)dτ δ iα, (3.9) B ( ) i = 1 2 S 1 ψ 2 ln r dτ e 1 B (+) i = 1 (3.10) 2 S ψ 1 ln r dτ e r is the distance between node i and node and e indicates integration along the th element. In the geometry representation by linear elements, the unit normals, n are constants along each element and non-continuous between the elements. Hence, care must be taken to distinguish between the unknowns ( ) Φ ( ) = Φn 1, x = x ( Φ ) (+) = Φn, x = x (3.11) The index should be understood as an element index for S and n, and as a global node index for ( Φ ), Φ and x. The components of the unit normal n are related to the element nodes location by n = 1 (z +1 z,y y +1 ). (3.12) S In previous formulations, the coefficients B ( ) i and B (+) i were united to a single coefficient B i, which multiplied the single unknown, ( ) Φ, i.e. equality of ( ) Φ ( ) and ( ) Φ (+) was assumed. This assumption caused numerical instability, when large angles developed between two adacent normals, during the development of steep waves [6].
6 246 N. Drimer, Y. Agnon / Wave Motion 43 (2006) In order to formulate a solution method which distinguish between the normal velocities, a compatibility condition relating ( ) Φ ( ) and ( ) Φ (+) is required. A tangential derivative, of any function f following a free surface particle, may be computed by (finite difference central derivative) f s = f +1 f 1, (3.13a) DS DS = [ (y +1 y 1 ) 2 + (z +1 z 1 ) 2] 1/2, (3.13b) where f is the value of f at the location of the th particle. To the same degree of approximation the slope of the tangential velocity of the th particle will be ) y z = (3.14) ( y z z s = y +1 y 1 z +1 z 1. This is the slope of a line from node 1 to node + l shown in Fig. 2. The true normal velocity of the th particle is the velocity component, which is normal to the tangent defined by (3.14). Hence, the components of the true th unit normal are n = 1 (z +1 z 1,y 1 y +1 ). (3.15) DS Since the velocity vector of any particle is unique, the following compatibility conditions have to be satisfied ( Φ ( Φ where β ( ) ( Φ s ) ( ) ) (+) = cos β ( ) = cos β (+) ( ) Φ ( ) Φ is the angle between n ( ) ) = Φ +1 Φ 1 DS sin β ( ) + sin β (+) ( ) Φ s ( ) Φ, (3.16) s and n as shown in Fig. 2. The tangential velocity of the th particle is by (3.13a) If the boundary Γ is smooth at node ( = π), the compatibility condition is reduced to ( ) Φ ( ) ( ) Φ (+) ( ) Φ = =. (3.17) (3.18) Hence, the normal velocity ( Φ should not be split into ( Φ, ( Φ along the straight segments of the solid or artificial boundaries. At the corners the free surface and the solid or artificial boundaries, ( Φ will be split, but since at those corners ) two boundary conditions should be satisfied, the compatibility condition (3.16) is not needed. Substituting the compatibility conditions for the internal free surface particles (3.16) in (3.8) we obtain N [A i Φ ] [ ( ) ] (B ( ) i + B (+) Φ i ) [( ( ) ) B ( ) Φ ( ) ( Φ ) ] (+) i + B (+) i =1 Γ corners Γa s [ ( ) ( ) ] B (n) Φ i + B (s) Φ i = 0, i = 1, 2,...N, (3.19) Γ s ) ) ( ) ) (+)
7 where B (n) i B (s) i N. Drimer, Y. Agnon / Wave Motion 43 (2006) = B ( ) i cos B ( ) + B (+) i cos β (+), (3.20) = B ( ) i sin B (+) + B ( ) i sin β ( ). (3.21) Since ( Φ s ) is related to Φ 1 and Φ +1 by (3.17), (3.19) is a linear system of N algebraic equations with 2N unknowns. For a given Φ or ( Φ at each nodal point (by the boundary conditions), at each time step, the system ) may be solved to find ( Φ ) or Φ, respectively The time stepping procedure Similar procedures were used in previous works where the main differences are in the method of solution of the boundary value problem, the method of time integration, and the use of Lagrangian or Eulerian representation of the free surface. We use the formulation, a fourth order Runge Kutta time integrating procedure, and the boundary value solver described in Section 3.1. The time stepping procedure may be described by the following flow chart. 4. Results and discussion First, the method is tested through some numerical examples, by comparison with previous results. In Section 4.1, the evaluation of a steep propagating wave is considered while in Section 4.2 a steep but steady standing wave is followed. These are two problems with quite different nature. In the first, the program is required to follow an extremely unstable wave, which breaks after about half a period. A plunging breaker with a thin et is created which causes sharp angles between adacent free surface elements. The standing wave is smooth, also steep, and here the program is tested to follow an exact (25th order) steady solution, during 100 periods. Next, some new results are presented. In Section 4.3, a numerical simulation of a wave channel with a piston wavemaker and an energy absorber is demonstrated. Shoaling of waves of finite length, up to breaking, and a comparison with published experimental results, are presented in Section 4.4.
8 248 N. Drimer, Y. Agnon / Wave Motion 43 (2006) Fig. 3. Breaking of a steep sinusoidal wave. Comparing with Jansen [6]. Parameters of initial wave: wave steepness ka = , water depth ka = (linear wave theory) Breaking of steep sinusoidal wave One very popular problem, treated by developers of methods for the evaluation of non-linear free surface flow, is the development of a plunging breaker from very steep sine wave. This situation is not very physical. However, it is proper for examining the performances of the different methods. Our results have been compared with previous simulations of overturning waves made using two different methods. The first model, presented by Jansen [6] used a method which is similar to that used for the present model, but without the new compatibility condition for the particles velocity. The second model, presented by New et al. [7] use an alternative numerical method. Both of the models make use of smoothing procedures to prevent numerical instabilities. The problem is defined by Eqs. (2.2) (2.7), and an initial solution which is a sine Stokes wave defined by a given steepness ka, and non-dimensional water depth kh, where k, a and h are the wavenumber, wave amplitude and water depth, respectively. Fig. 3a presents the free surface profile at t = 0.38T, where T is the linear wave period, for an initial wave parameter ka = , kh = Jansen s results, using 60 free surface particles are shown by a solid line, while results of the present model, using 90 free surface particles, are shown by dots representing particles. Shortly after this time Jansen s calculations break down due to the high curvature at the tip of the et. The calculations continued after a special treatment of two points at the tip of the et, which, in addition to a smoothing procedure applied every five time steps, enabled to follow the free surface until t = 0.46T, shown at Fig. 3b. The improved method presented here enabled us to follow the free surface until t = 0.46T without any definition of special points during the solution and with no use of any smoothing procedure. The improvement, formulated here, was guided by Jansen s inspection of the instability problem, and his attempt to treat it by the special points. In Table 1, the error in energy conservation during the solution is presented for Jansen s and for the present solution. Fig. 4 presents a solution of a similar problem which was evaluated for comparison with New et al. [7]; the parameters of the initial wave are: ka = and kh = As it may be seen, better agreement between the two models was achieved at this simulation. New et al. [7] presented the free profile at t = 0.438T and 0.477T. We present one more profile at t = 0.567T, shortly before the numerical model broke down. Fifty-nine and 60 free surface particles were used by New et al. [7] and in the present simulation, respectively. The relative error in conservation of energy was 1.2%. 2.4% and 7% at t = 0.438, and 0.567, respectively. In both simulations the initial time step was T/160 and was shortened during the solution, such that no panicle moves more than the distance to its nearest neighbor, during one time step (following Longuet-Higgins and Cokelet [1]).
9 Table 1 Error in energy conservation t/t E T / ( 1 4 ρga2 0 ) N. Drimer, Y. Agnon / Wave Motion 43 (2006) Relative error [E T (t) E T (0)]/E T (0) Jansen Present model Jansen Present model Steep standing wave Another verification problem is the simulation of standing wave in a rectangular basin. The model is tested to follow an exact solution developed by Schwartz and Whitney [11] up to 25th order in the wave steepness, for infinite water depth. The wave steepness was taken as the limit of convergence of the analytical solution, ka = 0.3. The wavelength and the basin length were set to 2π. The basin depth in the numerical model was set to 2π in order to approximate infinite depth. The free surface, basin bottom and each of the two vertical walls was divided into 32, 8 and 16 elements, respectively, i.e. 72 elements total. The time step was T 0 /32, where T 0 = is the period calculated by the analytical solution. The free profiles at t = 0 (initial analytical solution), t = 0.5T, 100T and 100.5T, are presented in Fig. 5. The errors in energy conservation, run-up and period are presented in Table 2. As may be seen, the numerical model keeps reasonable accuracy along 100 periods, with quite a small number of elements A numerical wave channel One maor difficulty that arises in numerical simulations, as well as in physical experiments, in a wave channel, is the existence of reflections from the end of the channel. Le Mehaute [8] presented a method to simulate wave absorber by the addition of an artificial damping term, vφ, to the right-hand side of the dynamic boundary condition (2.4). Le Mehaute [8] demonstrated his method by means of comparison between an analytical linear solution and wave channel experiments. At this chapter, we use the artificial damping term to simulate a wave absorber at the end of a Fig. 4. Breaking of a steep sinusoidal wave. Comparison with New et al. [7]. Free surface profiles at t = 0.438T and 0.499T Parameters of initial wave: wave steepness ka = , water depth kh = (linear wave theory).
10 250 N. Drimer, Y. Agnon / Wave Motion 43 (2006) Fig. 5. Standing wave in a rectangular basin. Verification of the numerical solution s ability to follow on analytical solution, which was developed by Schwartz and Whitney [11] up-to 25th order in the wave steepness. Parameters: unit wavenumber k = 1, basin length Lk =2π, basin depth hk =2π (infinite), wave steepness ak = 0.3. numerical wave channel. The wave channel simulator is tested to develop a solution comparable to first and second order analytical solutions to the wavemaker problem. As was indicated by Le Mehaute [8] in order to decrease the reflections caused by the wave absorber, the artificial damping coefficient v must increase gently along the channel. After some numerical experiments we present the following example: The channel is of unit depth (at rest) h = 1, and a total length L = 32. The piston wavemaker starts with harmonic motion at t = 0. The piston displacement is given by { Ap, t 0 Y p = A p cos(ωt), t > 0 (4.1) where A p is the piston amplitude and ω is the angular frequency. ω corresponds to a wavelength λ = 8, by the linear dispersion relation. Table 2 Errors in the simulation of standing waves Time step t/t 0 No. of oscillations of the numerical solution Relative error Period Run-up Energy
11 N. Drimer, Y. Agnon / Wave Motion 43 (2006) Fig. 6. A simulation of a wave channel with a piston wavemaker and an energy absorber. The development of the free surface in time, (a e), up to a steady state, (f), which is compared to an analytical solution in the frequency domain, of second order, (g). Parameters: L/A p /λ(linear)/h = 32/0.05/8/1. Half of the channel, Y p < y < 16, is a physical region, while the other half, 16 < y < 32, is an energy absorbent region. The artificial damping coefficient v increases from v =0,aty =16tov = 0.5 ρω at y = 32, where ρ is the fluid density. Four figures are presented (6, 7, 8 and 9) for four wavemaker amplitudes, A p = 0.05, 0.1, 0.15 and 0.2, respectively. Each of the figures contains seven sub figures showing free surface profiles. The upper five show the wave development, after the initiation of the wavemaker, at T is the wavemaker period. At the sixth sub figure, the free surface profiles at 0.5T, 1.5T, 2.5T, 3.5T and 4.5T are shown by a dotted line (each dot represents a particle in the numerical model), and a steady solution results of a second order analytical solution in the frequency domain, is shown by a solid line. The lowest sub figure presents first and second order analytical solutions by a dotted and a solid line, respectively. The second order solution was calculated by means of an eigenfunctions expansion, as described by Drimer and Agnon [12]. As may be seen, the solution became steady after about four periods. The waves are almost completely damped while propagating one way through the damping region. Some amount of reflection by the end may be attenuated on the way back, from the end wall to the between the physical region and the damping region. For a small wavemaker amplitude, A p = 0.05, the linear solution is valid, and is almost not different from the second order solution. From the similarity between the non-linear solution at steady state (t = 19.5T) and the analytical solution,
12 252 N. Drimer, Y. Agnon / Wave Motion 43 (2006) Fig. 7. A simulation of a wave channel with a piston wavemaker and an energy absorber. The development of the free surface in time, (a e), up to a steady state, (f), which is compared to an analytical solution in the frequency domain, of second order, (g). Parameters: L/A p /λ(linear)/h = 32/0.10/8/1. it may be concluded that reflections from the damping region are quite small. While increasing the amplitude to 0.1, the non-linearity causes the crest to be and more pointed than the trough. Up to one wavelength away from the wavemaker the fully non-linear and the second order solution are very similar. Closer to the damping region small differences may be observed which may be caused by three reasons: some amount of reflections from the damping region, the analytical solution error by neglecting higher order terms and numerical error of the non-linear solution. Further increasing of the wavemaker amplitude results in increasing the difference between the solutions, also good agreement is still achieved near the wavemaker. For the largest amplitude, A p = 0.2 (Fig. 9), a small crest develops at the center of the trough, due to the second order solution. This small crest is not confirmed by the non-linear solution and it may be cancelled by higher order terms. The Shore Protection Manual presents a diagram by which a suitable wave theory may be chosen depending on the two parameters: H and h, where H and h are the wave height and the water depth, respectively. In our example, gt 2 gt 2 H gt 2 = , , and , for A p = 0.05, 0.1, 0.15 and 0.2, respectively. (H was calculated from the linear transfer function of the piston wavemaker.) According to the diagram, the first two cases (Figs. 6 and 7) fall in the region of second order Stokes theory, while the other two (Figs. 8 and 9) are recommended to be treated by a third order theory. Our examples
13 N. Drimer, Y. Agnon / Wave Motion 43 (2006) Fig. 8. A simulation of a wave channel with a piston wavemaker and an energy absorber. The development of the free surface in time, (a e), up to a steady state, (f), which is compared to an analytical solution in the frequency domain, of second order, (g). Parameters: L/A p /λ(linear)/h = 32/0.15/8/1. confirm the above recommendation except the possibility of using linear wave theory for the smallest amplitude (Fig. 6). The numerical wave channel was divided into 142 elements, 80 on the free surface, 8 on the wavemaker board, 6 on the reflecting wall and 48 on the bottom. The time step was T/ Wave shoaling A numerical solution of wave shoaling is very problematic. In a real shore, or in a physical model, the wave breaking causes the energy to be dissipated in a turbulent process, preventing wave reflection and generation of standing waves. In a numerical model a reflecting boundary condition is to be formulated in order to prevent waves reflection. One possible way to avoid the necessity of non-reflecting condition is the simulation of shoaling of a solitary wave. This single wave of infinite length may approach breaking before a significant of energy is reflected from the end of the channel. Besides this benefit, the solitary wave is a practical good approximation for shallow water wave, hence it was adopted by many researchers. Otta et al. [9], present a comprehensive parametric investigation of solitary wave shoaling, evaluated by a boundary element model. In Fig. 10, an example, evaluated by the present model, is presented.
14 254 N. Drimer, Y. Agnon / Wave Motion 43 (2006) Fig. 9. A simulation of a wave channel with a piston wavemaker and an energy absorber. The development of the free surface in time, (a e), up to a steady state, (f), which is compared to an analytical solution in the frequency domain, of second order, (g). Parameters: L/A p /λ(linear)/h = 32/0.20/8/1. The solitary wave is produced by a piston wavemaker, which processes the motion (Sugino et al. [13]) Y p A p tanh[ω(t t 0 )], where A p = 0.6h 0. g Ω = 0.5, h 0 t 0 = 5 h 0 g, and h 0 is the still water depth near the piston. The motion is stopped at t = 12 h 0 /g when the velocity approaches zero.
15 N. Drimer, Y. Agnon / Wave Motion 43 (2006) Fig. 10. Simulation of a solitary wave shoaling using a numerical wave channel. Free surface profiles at sequential times (t = 12, 14, 16, 18, 20, 21, 22, 23, 23.6, 24), until wave breaking. Units: length h 0 = 1, time (h 0 /g) 1/2 =1. The channel length, L =22h 0 is composed of a horizontal interval of length 2h 0 and an inclined interval of horizontal length 20h 0 and of slope S = The cannel is terminated by a reflecting wall. Fig. 10 presents successive free surface profiles up to a developed plunging breaker. The wave propagates until breaking while no significant energy is reflected by the quite close end wall. This example also demonstrates the ability of the new formulation to follow the free surface up to the development of an angle of almost 180 (a cusp) between elements, before a numerical instability occurs. Since the piston is stopped, the total energy in the channel is to be conserved. The relative error in energy conservation at t = 24 h 0 /g (the last profile presented) is E(t = 24) E(T = 12) E(t = 12) = For this example, the channel was divided into 121 elements, 80 on the free surface, 8 on the piston, 12 on the end wall, and 21 on the bottom. The time step was t = 30 1 h0 /g. The solitary wave is determined by two parameters wave height and water depth. Usually a design wave (in a deterministic approach) is determined by a given period and height at deep-water and its height variation versus the local depth relation is of interest. Since the design wave is one of a wave group the energy of the leading waves is to be absorbed. Numerical wave channel presented at Section 4.3 was used to demonstrate shoaling of a periodic wave group. The following parameters were chosen in order to compare with an experiment made by Hansen and Svendsen [14]: deep water wave height H = 37.5mm deep water wave period T = 1.67 s sea bottom slope S = 1:35 In order to reduce the time of calculations, the wave shoaling was evaluated up to a finite depth h 0 = 300 mm by linear wave theory according to which the wave height at depth h 0 is H 0 = 36.5 mm and the wavemaker period T 0 is equal T. The wavemaker motion was Y p = A p tanh(t/t 0 )cos(ωt)
16 256 N. Drimer, Y. Agnon / Wave Motion 43 (2006) Fig. 11. Simulation of wave shoaling, using o numerical waves channel. Longitudinal cross-sections of the channel with free surface profiles, at sequential times. where t 0 is a time constant which controls the amplitude increase, which was set to t 0 = T 0 in this example. Assuming the second wave will break in the physical region, and to a linear solution of the piston wavemaker the correspondence wavemaker amplitude is A p = 0.092h 0. The wavelength at h 0 is λ 0 = 8.868, by linear theory. The numerical channel is composed of a section of uniform depth h 0, of length λ 0, a section of sloped bottom 1:35 of length 35 (h 0 h min ), where h min = 0.21h 0, and a damping region of uniform depth h min of length λ 0. Along the damping region, the artificial damping coefficient v increases linearly from 0 to 0.3. Fig. 12. Simulation of wave shoaling, using o numerical waves channel. Longitudinal cross-sections of the channel with free surface profiles, at sequential times, focusing on the breaking zone.
17 N. Drimer, Y. Agnon / Wave Motion 43 (2006) Fig. 13. Wave shoaling. Wave s height vs. water depth. Comparison among Hansen and Svendsen s [14] experimental results, Stiassnie and Peregrine [15] theoretical results, and present model. Parameters: wave s height at deep-water 37.5 mm, wave period 1.67 s, bench slope 1:35. In scale longitudinal cross-sections of the channel, with free surface profiles at successive times are shown in Fig. 11. The first generated wave is absorbed by the damping region, while the second wave is broken in the physical region. Fig. 12 focuses on the breaker. In Fig. 13, the wave height versus water depth curve is presented and compared to experimental and theoretical results. Results of Hansen and Svendsen s [14] experiments were taken from Stiassnie and Peregrine [15] who used it to compare with theoretical results they developed by solving conservation equations, assuming gentle slope. Results of linear shoaling theory are also shown. Up to depth of 200 mm an agreement exists between the numerical results, the theoretical results and the experimental results. This confirms the use we made of linear theory for transformation from deep-water to depth h 0 = 300 mm. At about h = 200 mm the linear theory starts to differ from the reality while the finite amplitude theory of Stiassnie and Peregrine [15] and the present model continues to follow the experimental results. At h = 100 mm the symmetric wave assumption of the gentle slope theory does not hold and the wave height is overestimated. The results of the present model continue to follow the experimental results almost until wave breaking. Close to breaking the calculated wave is higher than the measured wave, possibly due to neglecting of viscosity. The breaking is characterized by a prostration of the curve. 5. Conclusions The application of a compatibility condition to the velocity of the free surface particles prevents numerical instabilities found in previous models, which used boundary element methods with linear elements. The source for the instabilities, as was detected by was non-uniqueness of the normal velocity at the nodal points, which breached the stability when the angle between the normals of two adacent elements became greater than about 40. A definition of a unique normal at each nodal point enables the development of almost 180 between two adacent normals. A new formulation, using the above compatibility condition, was derived and tested through some numerical examples. An artificial damping term, which is added to the dynamic free surface boundary condition, may be used for the simulation of an energy absorber at the end of a numerical wave channel. This is demonstrated by comparison of the non-linear solution with first and second order solutions of the piston wavemaker problem. Potential theory gives reasonable results for wave shoaling up to breaking. The use of an artificial energy absorber enables the simulation of shoaling of a periodic wave in a quite practical situation. This is demonstrated by a comparison of numerical results with published experimental results. References [1] M.S. Longuet-Higgins, E.D. Cokelet, The deformation of steep surface waves on water. I. A numerical method of computation, Proc. R. Soc. Lond. A 350 (1976) 1 26.
18 258 N. Drimer, Y. Agnon / Wave Motion 43 (2006) [2] J.W. Dold, An efficient surface-integral algorithm applied to unsteady gravity waves, J. Comp. Phys. 103 (1) (1992) [3] T. Vine, P. Brevig, Numerical simulation of breaking waves, Adv. Water Resour. 4 (1981) [4] J.W. Dold, D.H. Peregrine, An efficient boundary integral method for steep unsteady water waves, Numer. Methods Fluid Dyn. II (1986) [5] M.J. Cooker, A boundary-integral method for water wave motion over irregular beds, Eng. Anal. Boundary Elements 7 (4) (1990) [6] P.C.M. Jansen, A Boundary Element Model for Non-Linear Free Surface Phenomena, Delft University of Technology, Report No. 86-2, [7] A.L. New, P. McIver, D.H. Peregrine, Computations of overturning waves, J. Fluid Mech. 150 (1985) [8] B. Le Mehaute, Progressive wave absorber, J. Hydraulic Res. 10 (2) (1972) [9] A.K. Otta, I.A. Svendsen, S.T. Grilli, The breaking and run-up of solitary waves on beaches, Coastal Eng. 112 (1992) [10] J. Broeze, Computation of breaking waves with a panel method, J. Fluid Mech. 190 (1992) [11] L.W. Schwartz, A.K. Whitney, A semi-analytic solution for non-linear standing waves in deep water, J. Fluid Mech. 107 (1981) [12] N. Drimer, Y. Agnon, A hybrid boundary element method for second-order wave body interaction, Appl. Ocean Res. 16 (1994) [13] R. Sugino, H. Kawabata, N. Tosaka, Nonlinear free surface flow problems by boundary element Lagrangian solution procedure, in: Boundary Element Methods, Pergamon Press, [14] J.B. Hansen, I.A. Svendsen, Regular waves in shoaling water experimental data, Institute of Hydrodynamics and Hydraulic Engineering Technology, University of Denmark, Series Paper No. 21, [15] M. Stiassnie, D.H. Peregrine, Shoaling of finite-amplitude surface waves on water of slowly varying depth, J. Fluid Mech. 97 (1980) [16] M.L. Banner, D.H. Peregrine, Wave breaking in deep water, Ann. Rev. Fluid Mech. 25 (1993) [17] M. Perlin, J. He, L.P. Bernal, An experimental study of deep water plunging breakers, Phys. Fluids 8 (9) (1996) [18] G. Chen, C. Kharif, Two-dimensional Navier Stokes simulation of breaking waves, Phys. Fluids 11 (1) (1999) [19] H.A. Schaffer, P.A. Madsen, R. Deigaard, A Boussinesq model for waves breaking in shallow water, Coastal Eng. 20 (1993) [20] S.T. Grilli, P. Guyenne, F. Dias, A fully non-linear model for three-dimensional overturning waves over an arbitrary bottom, Int. J. Num. Methods Fluids 35 (7) (2001) [21] S.T. Grilli, R. Subramanya, Recent advances in the BEM modelling of nonlinear water waves, in: H. Power (Ed.), Boundary Element Applications in Fluid Mechanics, Advances in Fluid Mechanics Series, Computational Mechanics Publication, Southampton, UK, 1995, Chapter 4, pp [22] C.H. Kim, A.H. Clement, K. Tanizawa, Recent research and development of numerical wave tank a review, Int. J. Offshore Polar Eng. 9 (4) (1999) [23] P. Donescu, L.N. Virgin, An implicit boundary element solution with consistent linearization for free surface flows and non-linear fluid structure interaction of floating bodies, Int. J. Num. Methods Eng. 5 (4) (2001) [24] J.-B. Song, M.L. Banner, On determining the onset and strength of breaking for deep water waves. Part I: unforced irrotational wave groups, J. Phys. Oceanogr. 32 (9) (2002) [25] M.L. Banner, J.-B. Song, On determining the onset and strength of breaking for deep water waves. Part II: influence of wind forcing and surface shear, J. Phys. Oceanogr. 32 (9) (2002) [26] J.-B. Song, M.L. Banner, Influence of mean water depth and a subsurface sandbar on the onset and strength of wave breaking, J. Phys. Oceanogr. 34 (4) (2004)
A Two-Dimensional Numerical Wave Flume Part 1: Nonlinear Wave Generation, Propagation, and Absorption
S. F. Baudic Deepwater Systems Division, ABB Lummus Global Inc., Houston, TX 77242-2833 A. N. Williams Department of Civil and Environmental Engineering, University of Houston, Houston, TX 77204-4791 e-mail:
More informationGeneration and propagation of solitary wave over a steep sloping beach
Generation and propagation of solitary wave over a steep sloping beach Zouhaier Hafsia (1), Mehdi Ben Haj (2), Hedi Lamloumi (3), Khlifa Maalel (4) and Ridha Zgolli (5) (1-5) ENIT. Laboratoire d Hydraulique
More information1.2 Numerical Solutions of Flow Problems
1.2 Numerical Solutions of Flow Problems DIFFERENTIAL EQUATIONS OF MOTION FOR A SIMPLIFIED FLOW PROBLEM Continuity equation for incompressible flow: 0 Momentum (Navier-Stokes) equations for a Newtonian
More informationA Hybrid Technique Using Particle and Boundary-Element Methods for Wave-Body Interaction Problems
9th International Conference on Numerical Ship Hydrodynamics Ann Arbor, Michigan, August 5-8, 27 A Hybrid Technique Using Particle and Boundary-Element Methods for Wave-Body Interaction Problems Makoto
More informationExperimental Validation of the Computation Method for Strongly Nonlinear Wave-Body Interactions
Experimental Validation of the Computation Method for Strongly Nonlinear Wave-Body Interactions by Changhong HU and Masashi KASHIWAGI Research Institute for Applied Mechanics, Kyushu University Kasuga
More informationDevelopment of the Compliant Mooring Line Model for FLOW-3D
Flow Science Report 08-15 Development of the Compliant Mooring Line Model for FLOW-3D Gengsheng Wei Flow Science, Inc. October 2015 1. Introduction Mooring systems are common in offshore structures, ship
More informationThree-dimensional simulation of floating wave power device Xixi Pan 1, a, Shiming Wang 1, b, Yongcheng Liang 1, c
International Power, Electronics and Materials Engineering Conference (IPEMEC 2015) Three-dimensional simulation of floating wave power device Xixi Pan 1, a, Shiming Wang 1, b, Yongcheng Liang 1, c 1 Department
More information2D numerical simulation of ocean waves
2D numerical simulation of ocean waves Qingjie. Du,*, Y.C. Dennis. Leung Department of Mechanical Engineering, The University of Hong Kong, Hong Kong, China * Corresponding author. Tel: +852 51743593,
More informationNumerical Analysis of the Internal Kinematics and Dynamics of Three-dimensional Breaking Waves on Slopes ABSTRACT. Keyword: MATHEMATICAL FORMULATION
Proceedings of The Thirteenth (23) International Offshore and Polar Engineering Conference Honolulu, Hawaii, USA, May 25 3, 23 Copyright 23 by The International Society of Offshore and Polar Engineers
More informationNumerical Simulation of Coastal Wave Processes with the Use of Smoothed Particle Hydrodynamics (SPH) Method
Aristotle University of Thessaloniki Faculty of Engineering Department of Civil Engineering Division of Hydraulics and Environmental Engineering Laboratory of Maritime Engineering Christos V. Makris Dipl.
More informationMODELLING THE FLOW AROUND AN ISLAND AND A HEADLAND: APPLICATION OF A TWO MIXING LENGTH MODEL WITH TELEMAC3D. Nicolas Chini 1 and Peter K.
MODELLING THE FLOW AROUND AN ISLAND AND A HEADLAND: APPLICATION OF A TWO MIXING LENGTH MODEL WITH TELEMAC3D Nicolas Chini 1 and Peter K. Stansby 2 Numerical modelling of the circulation around islands
More informationWave Uplift on Platforms or Docks in Variable Depth Jiin-Jen Lee and C.P. Lai"
CHAPTER 149 ABSTRACT Wave Uplift on Platforms or Docks in Variable Depth Jiin-Jen Lee and C.P. Lai" Wave uplift forces on two dimensional platforms or docks in a variable depth medium has been studied
More informationLecture 1.1 Introduction to Fluid Dynamics
Lecture 1.1 Introduction to Fluid Dynamics 1 Introduction A thorough study of the laws of fluid mechanics is necessary to understand the fluid motion within the turbomachinery components. In this introductory
More informationCharacteristic Aspects of SPH Solutions
Characteristic Aspects of SPH Solutions for Free Surface Problems: Source and Possible Treatment of High Frequency Numerical Oscillations of Local Loads. A. Colagrossi*, D. Le Touzé & G.Colicchio* *INSEAN
More informationCS205b/CME306. Lecture 9
CS205b/CME306 Lecture 9 1 Convection Supplementary Reading: Osher and Fedkiw, Sections 3.3 and 3.5; Leveque, Sections 6.7, 8.3, 10.2, 10.4. For a reference on Newton polynomial interpolation via divided
More informationarxiv: v1 [physics.flu-dyn] 23 Oct 2018
Numerical models for evolution of extreme wave groups Eugeny Buldakov, Pablo Higuera 2 and Dimitris Stagonas 3 arxiv:8.9769v [physics.flu-dyn] 23 Oct 28 3 UCL, Department of Civil Engineering, Gower Street,
More informationQuantifying the Dynamic Ocean Surface Using Underwater Radiometric Measurement
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Quantifying the Dynamic Ocean Surface Using Underwater Radiometric Measurement Lian Shen Department of Mechanical Engineering
More informationCFD MODELING FOR PNEUMATIC CONVEYING
CFD MODELING FOR PNEUMATIC CONVEYING Arvind Kumar 1, D.R. Kaushal 2, Navneet Kumar 3 1 Associate Professor YMCAUST, Faridabad 2 Associate Professor, IIT, Delhi 3 Research Scholar IIT, Delhi e-mail: arvindeem@yahoo.co.in
More informationTable of contents for: Waves and Mean Flows by Oliver Bühler Cambridge University Press 2009 Monographs on Mechanics. Contents.
Table of contents for: Waves and Mean Flows by Oliver Bühler Cambridge University Press 2009 Monographs on Mechanics. Preface page 2 Part I Fluid Dynamics and Waves 7 1 Elements of fluid dynamics 9 1.1
More informationMESHLESS SOLUTION OF INCOMPRESSIBLE FLOW OVER BACKWARD-FACING STEP
Vol. 12, Issue 1/2016, 63-68 DOI: 10.1515/cee-2016-0009 MESHLESS SOLUTION OF INCOMPRESSIBLE FLOW OVER BACKWARD-FACING STEP Juraj MUŽÍK 1,* 1 Department of Geotechnics, Faculty of Civil Engineering, University
More informationCONFORMAL MAPPING POTENTIAL FLOW AROUND A WINGSECTION USED AS A TEST CASE FOR THE INVISCID PART OF RANS SOLVERS
European Conference on Computational Fluid Dynamics ECCOMAS CFD 2006 P. Wesseling, E. Oñate and J. Périaux (Eds) c TU Delft, The Netherlands, 2006 CONFORMAL MAPPING POTENTIAL FLOW AROUND A WINGSECTION
More informationCoupling of STAR-CCM+ to Other Theoretical or Numerical Solutions. Milovan Perić
Coupling of STAR-CCM+ to Other Theoretical or Numerical Solutions Milovan Perić Contents The need to couple STAR-CCM+ with other theoretical or numerical solutions Coupling approaches: surface and volume
More informationNumerical Wave Tank Modeling of Hydrodynamics of Permeable Barriers
ICHE 2014, Hamburg - Lehfeldt & Kopmann (eds) - 2014 Bundesanstalt für Wasserbau ISBN 978-3-939230-32-8 Numerical Wave Tank Modeling of Hydrodynamics of Permeable Barriers K. Rajendra & R. Balaji Indian
More informationAvailable online at ScienceDirect. Procedia Engineering 126 (2015 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 6 (05 ) 660 664 7th International Conference on Fluid Mechanics, ICFM7 Simulate wave body interaction based on the incompressible
More informationDYNAMIC ANALYSIS OF A GENERATOR ON AN ELASTIC FOUNDATION
DYNAMIC ANALYSIS OF A GENERATOR ON AN ELASTIC FOUNDATION 7 DYNAMIC ANALYSIS OF A GENERATOR ON AN ELASTIC FOUNDATION In this tutorial the influence of a vibrating source on its surrounding soil is studied.
More informationExample 13 - Shock Tube
Example 13 - Shock Tube Summary This famous experiment is interesting for observing the shock-wave propagation. Moreover, this case uses the representation of perfect gas and compares the different formulations:
More informationWave Hydro Dynamics Prof. V. Sundar Department of Ocean Engineering Indian Institute of Technology, Madras
Wave Hydro Dynamics Prof. V. Sundar Department of Ocean Engineering Indian Institute of Technology, Madras Module No. # 06 Wave Theories and Testing Facilities Lecture No. # 01 Finite Amplitude Wave Theories
More informationRealtime Water Simulation on GPU. Nuttapong Chentanez NVIDIA Research
1 Realtime Water Simulation on GPU Nuttapong Chentanez NVIDIA Research 2 3 Overview Approaches to realtime water simulation Hybrid shallow water solver + particles Hybrid 3D tall cell water solver + particles
More informationAcknowledgements. Prof. Dan Negrut Prof. Darryl Thelen Prof. Michael Zinn. SBEL Colleagues: Hammad Mazar, Toby Heyn, Manoj Kumar
Philipp Hahn Acknowledgements Prof. Dan Negrut Prof. Darryl Thelen Prof. Michael Zinn SBEL Colleagues: Hammad Mazar, Toby Heyn, Manoj Kumar 2 Outline Motivation Lumped Mass Model Model properties Simulation
More informationFree-surface flow under a sluice gate from deep water
Free-surface flow under a sluice gate from deep water 1,2 L.H. Wiryanto, 1 J. Widjaja, 3 B. Supriyanto 1 Faculty of Mathematics and Natural Sciences, 3 Faculty of Art and Design, Bandung Institute of Technology
More informationCHAPTER 3. Elementary Fluid Dynamics
CHAPTER 3. Elementary Fluid Dynamics - Understanding the physics of fluid in motion - Derivation of the Bernoulli equation from Newton s second law Basic Assumptions of fluid stream, unless a specific
More informationAalborg Universitet. Numerical 3-D Modelling of Overflows Larsen, Torben; Nielsen, L.; Jensen, B.; Christensen, E. D.
Aalborg Universitet Numerical 3-D Modelling of Overflows Larsen, Torben; Nielsen, L.; Jensen, B.; Christensen, E. D. Published in: Confernce Proceedings : 11th International Conference on Urban Drainage
More informationAurélien Thinat Stéphane Cordier 1, François Cany
SimHydro 2012:New trends in simulation - Hydroinformatics and 3D modeling, 12-14 September 2012, Nice Aurélien Thinat, Stéphane Cordier, François Cany Application of OpenFOAM to the study of wave loads
More informationNonlinear Potential Flow Solver Development in OpenFOAM
Nonlinear Potential Flow Solver Development in OpenFOAM A. Mehmood Plymouth University, UK April 19,2016 A. Mehmood Table of Contents 1 Motivation 2 Solution Methodology Mathematical Formulation Sequence
More informationC. Fochesato*, F. Dias**, S. Grilli*** * Mathématiques Appliquées de Bordeaux (Université de Bordeaux),
Wave EnergyFocusing in athree-dimensional Numerical WaveTank C. Fochesato* F. Dias** S. Grilli*** * Mathématiques Appliquées de Bordeau Université de Bordeau) Talence France ** Centre de Mathématiques
More informationDriven Cavity Example
BMAppendixI.qxd 11/14/12 6:55 PM Page I-1 I CFD Driven Cavity Example I.1 Problem One of the classic benchmarks in CFD is the driven cavity problem. Consider steady, incompressible, viscous flow in a square
More informationRotating Moving Boundary Analysis Using ANSYS 5.7
Abstract Rotating Moving Boundary Analysis Using ANSYS 5.7 Qin Yin Fan CYBERNET SYSTEMS CO., LTD. Rich Lange ANSYS Inc. As subroutines in commercial software, APDL (ANSYS Parametric Design Language) provides
More informationNumerical investigation of free surface elevation and celerity of solitary waves passing over submerged trapezoidal breakwaters
Int. Jnl. of Multiphysics Volume 9 Number 1 2015 61 Numerical investigation of free surface elevation and celerity of solitary waves passing over submerged trapezoidal breakwaters Parviz Ghadimi 1 *, Amin
More informationCoastal impact of a tsunami Review of numerical models
Coastal impact of a tsunami Review of numerical models Richard Marcer 2 Content Physics to simulate Different approaches of modelling 2D depth average Full 3D Navier-Stokes 3D model Key point : free surface
More informationA simple method for potential flow simulation of cascades
Sādhanā Vol. 35, Part 6, December 2010, pp. 649 657. Indian Academy of Sciences A simple method for potential flow simulation of cascades ARAVIND BHIMARASETTY and RAGHURAMAN N GOVARDHAN Department of Mechanical
More informationSPH: Towards the simulation of wave-body interactions in extreme seas
SPH: Towards the simulation of wave-body interactions in extreme seas Guillaume Oger, Mathieu Doring, Bertrand Alessandrini, and Pierre Ferrant Fluid Mechanics Laboratory (CNRS UMR6598) Ecole Centrale
More informationThe Level Set Method. Lecture Notes, MIT J / 2.097J / 6.339J Numerical Methods for Partial Differential Equations
The Level Set Method Lecture Notes, MIT 16.920J / 2.097J / 6.339J Numerical Methods for Partial Differential Equations Per-Olof Persson persson@mit.edu March 7, 2005 1 Evolving Curves and Surfaces Evolving
More informationKeywords: flows past a cylinder; detached-eddy-simulations; Spalart-Allmaras model; flow visualizations
A TURBOLENT FLOW PAST A CYLINDER *Vít HONZEJK, **Karel FRAŇA *Technical University of Liberec Studentská 2, 461 17, Liberec, Czech Republic Phone:+ 420 485 353434 Email: vit.honzejk@seznam.cz **Technical
More informationSloshing reduction effect of splitting wall in cylindrical tank
Sloshing reduction effect of splitting wall in cylindrical tank D.Takabatake 1, S.Sawada 2, N.Yoneyama 2 and M.Miura 3 1 Central Research Institute of Electric Power Industry, 1646,Abiko, Abiko, Chiba,
More informationLagrangian and Eulerian Representations of Fluid Flow: Kinematics and the Equations of Motion
Lagrangian and Eulerian Representations of Fluid Flow: Kinematics and the Equations of Motion James F. Price Woods Hole Oceanographic Institution Woods Hole, MA, 02543 July 31, 2006 Summary: This essay
More informationApplication of Finite Volume Method for Structural Analysis
Application of Finite Volume Method for Structural Analysis Saeed-Reza Sabbagh-Yazdi and Milad Bayatlou Associate Professor, Civil Engineering Department of KNToosi University of Technology, PostGraduate
More informationQUASI-3D SOLVER OF MEANDERING RIVER FLOWS BY CIP-SOROBAN SCHEME IN CYLINDRICAL COORDINATES WITH SUPPORT OF BOUNDARY FITTED COORDINATE METHOD
QUASI-3D SOLVER OF MEANDERING RIVER FLOWS BY CIP-SOROBAN SCHEME IN CYLINDRICAL COORDINATES WITH SUPPORT OF BOUNDARY FITTED COORDINATE METHOD Keisuke Yoshida, Tadaharu Ishikawa Dr. Eng., Tokyo Institute
More informationLevel-set and ALE Based Topology Optimization Using Nonlinear Programming
10 th World Congress on Structural and Multidisciplinary Optimization May 19-24, 2013, Orlando, Florida, USA Level-set and ALE Based Topology Optimization Using Nonlinear Programming Shintaro Yamasaki
More informationQuantifying the Dynamic Ocean Surface Using Underwater Radiometric Measurement
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Quantifying the Dynamic Ocean Surface Using Underwater Radiometric Measurement Lian Shen Department of Mechanical Engineering
More information1 Drevard, Marcer, Grilli, Fraunié and Rey
EXPERIMENTAL VALIDATION OF A COUPLED BEM-NAVIER-STOKES MODEL FOR SOLITARY WAVE SHOALING AND BREAKING D. Drevard 1,R. Marcer 2, S.T. Grilli 3, M. ASCE, P. Fraunié 1 and V. Rey 1 Abstract: This paper studies
More informationWAVE PATTERNS, WAVE INDUCED FORCES AND MOMENTS FOR A GRAVITY BASED STRUCTURE PREDICTED USING CFD
Proceedings of the ASME 2011 30th International Conference on Ocean, Offshore and Arctic Engineering OMAE2011 June 19-24, 2011, Rotterdam, The Netherlands OMAE2011-49593 WAVE PATTERNS, WAVE INDUCED FORCES
More informationA Direct Simulation-Based Study of Radiance in a Dynamic Ocean
1 DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. A Direct Simulation-Based Study of Radiance in a Dynamic Ocean LONG-TERM GOALS Dick K.P. Yue Center for Ocean Engineering
More informationComparison of a two-dimensional viscid and inviscid model for rotating stall analysis
Comparison of a two-dimensional viscid and inviscid model for rotating stall analysis S. LJEVAR, H.C. DE LANGE, A.A. VAN STEENHOVEN Department of Mechanical Engineering Eindhoven University of Technology
More informationJournal of Theoretical and Applied Mechanics, Sofia, 2015, vol. 45, No. 2, pp
Journal of Theoretical and Applied Mechanics, Sofia, 2015, vol. 45, No. 2, pp. 59 74 SOLID MECHANICS WAVE PROPAGATION DUE TO AN EMBEDDED SEISMIC SOURCE IN A GRADED HALF-PLANE WITH RELIEF PECULIARITIES.
More informationTechnical Report TR
Technical Report TR-2015-09 Boundary condition enforcing methods for smoothed particle hydrodynamics Arman Pazouki 1, Baofang Song 2, Dan Negrut 1 1 University of Wisconsin-Madison, Madison, WI, 53706-1572,
More informationNumerical and experimental investigations into liquid sloshing in a rectangular tank
The 2012 World Congress on Advances in Civil, Environmental, and Materials Research (ACEM 12) Seoul, Korea, August 26-30, 2012 Numerical and experimental investigations into liquid sloshing in a rectangular
More informationFLOWING FLUIDS AND PRESSURE VARIATION
Chapter 4 Pressure differences are (often) the forces that move fluids FLOWING FLUIDS AND PRESSURE VARIATION Fluid Mechanics, Spring Term 2011 e.g., pressure is low at the center of a hurricane. For your
More informationStudies of the Continuous and Discrete Adjoint Approaches to Viscous Automatic Aerodynamic Shape Optimization
Studies of the Continuous and Discrete Adjoint Approaches to Viscous Automatic Aerodynamic Shape Optimization Siva Nadarajah Antony Jameson Stanford University 15th AIAA Computational Fluid Dynamics Conference
More information2-D Tank Sloshing Using the Coupled Eulerian- LaGrangian (CEL) Capability of Abaqus/Explicit
2-D Tank Sloshing Using the Coupled Eulerian- LaGrangian (CEL) Capability of Abaqus/Explicit Jeff D. Tippmann, Sharat C. Prasad 2, and Parthiv N. Shah ATA Engineering, Inc. San Diego, CA 923 2 Dassault
More informationCommutative filters for LES on unstructured meshes
Center for Turbulence Research Annual Research Briefs 1999 389 Commutative filters for LES on unstructured meshes By Alison L. Marsden AND Oleg V. Vasilyev 1 Motivation and objectives Application of large
More informationNon-Newtonian Transitional Flow in an Eccentric Annulus
Tutorial 8. Non-Newtonian Transitional Flow in an Eccentric Annulus Introduction The purpose of this tutorial is to illustrate the setup and solution of a 3D, turbulent flow of a non-newtonian fluid. Turbulent
More informationITTC Recommended Procedures and Guidelines
Page 1 of 13 Table of Contents Codes in the Frequency Domain... 2 1. PURPOSE OF PROCEDURE... 2 2. SCOPE... 2 2.1 Present Theoretical Limitations... 4 2.2 Numerical Aspects... 4 2.3 Software Engineering
More informationGeometric Acoustics in High-Speed Boundary Layers
Accepted for presentation at the 9th International Symposium on Shock Waves. Madison, WI. July -9,. Paper #8 Geometric Acoustics in High-Speed Boundary Layers N. J. Parziale, J. E. Shepherd, and H. G.
More informationSupercomputing of Tsunami Damage Mitigation Using Offshore Mega-Floating Structures
International Innovation Workshop on Tsunami, Snow Avalanche and Flash Flood Energy Dissipation January 21-22, 2016, Maison Villemanzy in Lyon, France Supercomputing of Tsunami Damage Mitigation Using
More informationA Direct Simulation-Based Study of Radiance in a Dynamic Ocean
A Direct Simulation-Based Study of Radiance in a Dynamic Ocean Lian Shen Department of Civil Engineering Johns Hopkins University Baltimore, MD 21218 phone: (410) 516-5033 fax: (410) 516-7473 email: LianShen@jhu.edu
More informationTowards Real Time Simulation of Ship-Ship Interaction
Downloaded from orbit.dtu.dk on: Nov 12, 2018 Towards Real Time Simulation of Ship-Ship Interaction Lindberg, Ole; Bingham, Harry B.; Engsig-Karup, Allan Peter; Madsen, Per A. Published in: International
More informationFluid-structure Interaction by the mixed SPH-FE Method with Application to Aircraft Ditching
Fluid-structure Interaction by the mixed SPH-FE Method with Application to Aircraft Ditching Paul Groenenboom ESI Group Delft, Netherlands Martin Siemann German Aerospace Center (DLR) Stuttgart, Germany
More informationChapter 6. Semi-Lagrangian Methods
Chapter 6. Semi-Lagrangian Methods References: Durran Chapter 6. Review article by Staniford and Cote (1991) MWR, 119, 2206-2223. 6.1. Introduction Semi-Lagrangian (S-L for short) methods, also called
More informationMid-Year Report. Discontinuous Galerkin Euler Equation Solver. Friday, December 14, Andrey Andreyev. Advisor: Dr.
Mid-Year Report Discontinuous Galerkin Euler Equation Solver Friday, December 14, 2012 Andrey Andreyev Advisor: Dr. James Baeder Abstract: The focus of this effort is to produce a two dimensional inviscid,
More informationInviscid Flows. Introduction. T. J. Craft George Begg Building, C41. The Euler Equations. 3rd Year Fluid Mechanics
Contents: Navier-Stokes equations Inviscid flows Boundary layers Transition, Reynolds averaging Mixing-length models of turbulence Turbulent kinetic energy equation One- and Two-equation models Flow management
More informationA stable moving-particle semi-implicit method for free surface flows
Fluid Dynamics Research 38 (2006) 241 256 A stable moving-particle semi-implicit method for free surface flows B. Ataie-Ashtiani, Leila Farhadi Department of Civil Engineering, Sharif University of Technology,
More information2.7 Cloth Animation. Jacobs University Visualization and Computer Graphics Lab : Advanced Graphics - Chapter 2 123
2.7 Cloth Animation 320491: Advanced Graphics - Chapter 2 123 Example: Cloth draping Image Michael Kass 320491: Advanced Graphics - Chapter 2 124 Cloth using mass-spring model Network of masses and springs
More informationExample 24 Spring-back
Example 24 Spring-back Summary The spring-back simulation of sheet metal bent into a hat-shape is studied. The problem is one of the famous tests from the Numisheet 93. As spring-back is generally a quasi-static
More informationA Singular Example for the Averaged Mean Curvature Flow
To appear in Experimental Mathematics Preprint Vol. No. () pp. 3 7 February 9, A Singular Example for the Averaged Mean Curvature Flow Uwe F. Mayer Abstract An embedded curve is presented which under numerical
More information3D Finite Element Software for Cracks. Version 3.2. Benchmarks and Validation
3D Finite Element Software for Cracks Version 3.2 Benchmarks and Validation October 217 1965 57 th Court North, Suite 1 Boulder, CO 831 Main: (33) 415-1475 www.questintegrity.com http://www.questintegrity.com/software-products/feacrack
More informationNUMERICAL DIFFUSION AND DISSIPATION IN HYDROSTATIC MODELS OF INTERNAL WAVES
NUMERICAL DIFFUSION AND DISSIPATION IN HYDROSTATIC MODELS OF INTERNAL WAVES Ben R. Hodges 1 (A.M., ASCE) and Sarah Kelly Delavan 2 ABSTRACT Analysis of numerical diffusion and dissipation rates in a hydrostatic
More informationInvestigation of cross flow over a circular cylinder at low Re using the Immersed Boundary Method (IBM)
Computational Methods and Experimental Measurements XVII 235 Investigation of cross flow over a circular cylinder at low Re using the Immersed Boundary Method (IBM) K. Rehman Department of Mechanical Engineering,
More informationRevision of the SolidWorks Variable Pressure Simulation Tutorial J.E. Akin, Rice University, Mechanical Engineering. Introduction
Revision of the SolidWorks Variable Pressure Simulation Tutorial J.E. Akin, Rice University, Mechanical Engineering Introduction A SolidWorks simulation tutorial is just intended to illustrate where to
More informationComputation of Velocity, Pressure and Temperature Distributions near a Stagnation Point in Planar Laminar Viscous Incompressible Flow
Excerpt from the Proceedings of the COMSOL Conference 8 Boston Computation of Velocity, Pressure and Temperature Distributions near a Stagnation Point in Planar Laminar Viscous Incompressible Flow E. Kaufman
More informationcuibm A GPU Accelerated Immersed Boundary Method
cuibm A GPU Accelerated Immersed Boundary Method S. K. Layton, A. Krishnan and L. A. Barba Corresponding author: labarba@bu.edu Department of Mechanical Engineering, Boston University, Boston, MA, 225,
More informationMoving Interface Problems: Methods & Applications Tutorial Lecture II
Moving Interface Problems: Methods & Applications Tutorial Lecture II Grétar Tryggvason Worcester Polytechnic Institute Moving Interface Problems and Applications in Fluid Dynamics Singapore National University,
More informationON THE VELOCITY OF A WEIGHTED CYLINDER DOWN AN INCLINED PLANE
ON THE VELOCITY OF A WEIGHTED CYLINDER DOWN AN INCLINED PLANE Raghav Grover and Aneesh Agarwal RG (Grade 12 High School), AA (Grade 11 High School) Department of Physics, The Doon School, Dehradun. raghav.503.2019@doonschool.com,
More informationUnstructured Mesh Generation for Implicit Moving Geometries and Level Set Applications
Unstructured Mesh Generation for Implicit Moving Geometries and Level Set Applications Per-Olof Persson (persson@mit.edu) Department of Mathematics Massachusetts Institute of Technology http://www.mit.edu/
More informationIN A NUMERICAL WAVE TANK. Abstract: Tsunami generation by underwater landslides is simulated in a
LANDSLIDE TSUNAMI AMPLITUDE PREDICTION IN A NUMERICAL WAVE TANK Stephan T. Grilli 1, M.ASCE, Sylvia Vogelmann 1, and Philip Watts 2 Abstract: Tsunami generation by underwater landslides is simulated in
More informationA spectral boundary integral method for inviscid water waves in a finite domain
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 25; : 3 Published online in Wiley InterScience (www.interscience.wiley.com). DOI:.2/fld A spectral boundary integral method
More informationLATTICE-BOLTZMANN METHOD FOR THE SIMULATION OF LAMINAR MIXERS
14 th European Conference on Mixing Warszawa, 10-13 September 2012 LATTICE-BOLTZMANN METHOD FOR THE SIMULATION OF LAMINAR MIXERS Felix Muggli a, Laurent Chatagny a, Jonas Lätt b a Sulzer Markets & Technology
More informationNUMERICAL MODELING STUDY FOR FLOW PATTERN CHANGES INDUCED BY SINGLE GROYNE
NUMERICAL MODELING STUDY FOR FLOW PATTERN CHANGES INDUCED BY SINGLE GROYNE Jungseok Ho 1, Hong Koo Yeo 2, Julie Coonrod 3, and Won-Sik Ahn 4 1 Research Assistant Professor, Dept. of Civil Engineering,
More informationNumerical simulations of dam-break floods with MPM
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 00 (206) 000 000 www.elsevier.com/locate/procedia st International Conference on the Material Point Method, MPM 207 Numerical
More informationIncompressible Viscous Flow Simulations Using the Petrov-Galerkin Finite Element Method
Copyright c 2007 ICCES ICCES, vol.4, no.1, pp.11-18, 2007 Incompressible Viscous Flow Simulations Using the Petrov-Galerkin Finite Element Method Kazuhiko Kakuda 1, Tomohiro Aiso 1 and Shinichiro Miura
More informationDetermination of Free Surface in Steady-State Seepage through a Dam with Toe Drain by the Boundary Element Method
Determination of Free Surface in Steady-State Seepage through a Dam with Toe Drain by the Boundary Element Method Somchart Chantasiriwan Faculty of Engineering, Thammasat University, Rangsit Campus, Pathum
More informationVerification and Validation for Seismic Wave Propagation Problems
Chapter 26 Verification and Validation for Seismic Wave Propagation Problems (1989-2-24-25-28-29-21-211-217-) (In collaboration with Dr. Nima Tafazzoli, Dr. Federico Pisanò, Mr. Kohei Watanabe and Mr.
More informationOffshore Platform Fluid Structure Interaction (FSI) Simulation
Offshore Platform Fluid Structure Interaction (FSI) Simulation Ali Marzaban, CD-adapco Murthy Lakshmiraju, CD-adapco Nigel Richardson, CD-adapco Mike Henneke, CD-adapco Guangyu Wu, Chevron Pedro M. Vargas,
More information21. Efficient and fast numerical methods to compute fluid flows in the geophysical β plane
12th International Conference on Domain Decomposition Methods Editors: Tony Chan, Takashi Kako, Hideo Kawarada, Olivier Pironneau, c 2001 DDM.org 21. Efficient and fast numerical methods to compute fluid
More informationPhase-field simulation of two-phase micro-flows in a Hele-Shaw cell
Computational Methods in Multiphase Flow III 7 Phase-field simulation of two-phase micro-flows in a Hele-Shaw cell Y. Sun & C. Beckermann Department of Mechanical and Industrial Engineering, University
More informationIMPROVED WALL BOUNDARY CONDITIONS WITH IMPLICITLY DEFINED WALLS FOR PARTICLE BASED FLUID SIMULATION
6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 1115 June 2018, Glasgow, UK IMPROVED WALL BOUNDARY CONDITIONS WITH IMPLICITLY
More informationTransactions on Modelling and Simulation vol 10, 1995 WIT Press, ISSN X
Hydrodynamic coefficients and motions due to a floating cylinder in waves D.D. Bhatta, M. Rahman Department of Applied Mathematics, Technical University of Nova Scotia, Halifax, Nova Scotia, Canada B3J
More informationFEMLAB Exercise 1 for ChE366
FEMLAB Exercise 1 for ChE366 Problem statement Consider a spherical particle of radius r s moving with constant velocity U in an infinitely long cylinder of radius R that contains a Newtonian fluid. Let
More informationA Direct Simulation-Based Study of Radiance in a Dynamic Ocean
A Direct Simulation-Based Study of Radiance in a Dynamic Ocean Dick K.P. Yue Center for Ocean Engineering Massachusetts Institute of Technology Room 5-321, 77 Massachusetts Ave, Cambridge, MA 02139 phone:
More informationSimulation of Freak Wave Impact Using the Higher Order Spectrum
Simulation of Freak Wave Impact Using the Higher Order Spectrum The Naval Hydro Pack Hrvoje Jasak and Vuko Vukčević Faculty of Mechanical Engineering and Naval Architecture, Uni Zagreb, Croatia Wikki Ltd.
More informationPossibility of Implicit LES for Two-Dimensional Incompressible Lid-Driven Cavity Flow Based on COMSOL Multiphysics
Possibility of Implicit LES for Two-Dimensional Incompressible Lid-Driven Cavity Flow Based on COMSOL Multiphysics Masanori Hashiguchi 1 1 Keisoku Engineering System Co., Ltd. 1-9-5 Uchikanda, Chiyoda-ku,
More information