Nonlinear evolution analysis of T-S disturbance wave at finite amplitude in nonparallel boundary layers

ＩｎｔｒｏｄｕｃｔｉｏｎＮｏｎｌｉｎｅａｒｓｔａｂｉｌｉｔｙｐｒｏｂｌｅｍｉｎｎｏｎｐａｒａｌｌｅｌｂｏｕｎｄａｒｙｌａｙｅｒｓｎｅｅｄｓｍａｉｎｌｙｃｏｎｓｉｄｅｒｔｈｅｅｖｏｌｕｔｉｏｎｏｆｕｎｓｔａｂｉｌｉｔｙＴ＿Ｓｗａｖｅａｎｄｉｎｔｅｒａｃｔｉｏｎｂｅｔｗｅｅｎｗａｖｅｓｏｆｄｉｆｆｅｒｅｎｔｆｒｅｑｕｅｎｃｉｅｓ.Ｉｔｉｓｖｅｒｙｉｍｐｏｒｔａｎｔｔｏｓｉｍｕｌａｔｅａｃｃｕｒａｔｅｌｙｓｐａｔｉａｌｅｖｏｌｕｔｉｏｎｒｏｌｅｏｆｄｉｓｔｕｒｂａｎｃｅｉｎｂｏｕｎｄａｒｙｌａｙｅｒｓ,…     

A horizontally curvilinear non‐hydrostatic model for simulating nonlinear wave motion in curved boundaries

A horizontally curvilinear non‐hydrostatic free surface model that embeds the second‐order projection method, the so‐called θ scheme, in fractional time stepping is developed to simulate nonlinear wave motion in curved boundaries. The model solves the unsteady, Navier–Stokes equations in a three‐dimensional curvilinear domain by incorporating the kinematic free surface boundary condition with a top‐layer boundary condition, which has been developed to improve the numerical accuracy and efficiency of the non‐hydrostatic model in the standard staggered grid layout. The second‐order Adams–Bashforth scheme with the third‐order spatial upwind method is implemented in discretizing advection terms. Numerical accuracy in terms of nonlinear phase speed and amplitude is verified against the nonlinear Stokes wave theory over varying wave steepness in a two‐dimensional numerical wave tank. The model is then applied to investigate the nonlinear wave characteristics in the presence of dispersion caused by reflection and diffraction in a semicircular channel. The model results agree quantitatively with superimposed analytical solutions. Finally, the model is applied to simulate nonlinear wave run‐ups caused by wave‐body interaction around a bottom‐mounted cylinder. The numerical results exhibit good agreement with experimental data and the second‐order diffraction theory. Overall, it is shown that the developed model, with only three vertical layers, is capable of accurately simulating nonlinear waves interacting within curved boundaries. Copyright © 2011 John Wiley & Sons, Ltd.     

Outflow boundary conditions for one-dimensional finite element modeling of blood flow and pressure waves in arteries

Flow and pressure waves, originating due to the contraction of the heart, propagate along the deformable vessels and reflect due to tapering, branching, and other discontinuities. The size and complexity of the cardiovascular system necessitate a “multiscale” approach, with “upstream” regions of interest (large arteries) coupled to reduced-order models of “downstream” vessels. Previous efforts to couple upstream and downstream domains have included specifying resistance and impedance outflow boundary conditions for the nonlinear one-dimensional wave propagation equations and iterative coupling between three-dimensional and one-dimensional numerical methods. We have developed a new approach to solve the one-dimensional nonlinear equations of blood flow in elastic vessels utilizing a space-time finite element method with GLS-stabilization for the upstream domain, and a boundary term to couple to the downstream domain. The outflow boundary conditions are derived following an approach analogous to the Dirichlet-to-Neumann (DtN) method. In the downstream domain, we solve simplified zero/one-dimensional equations to derive relationships between pressure and flow accommodating periodic and transient phenomena with a consistent formulation for different boundary condition types. In this paper, we also present a new boundary condition that accommodates transient phenomena based on a Green’s function solution of the linear, damped wave equation in the downstream domain.     

A meshless numerical wave tank for simulation of nonlinear irregular waves in shallow water

Time domain simulation of the interaction between offshore structures and irregular waves in shallow water becomes a focus due to significant increase of liquefied natural gas (LNG) terminals. To obtain the time series of irregular waves in shallow water, a numerical wave tank is developed by using the meshless method for simulation of 2D nonlinear irregular waves propagating from deep water to shallow water. Using the fundamental solution of Laplace equation as the radial basis function (RBF) and locating the source points outside the computational domain, the problem of water wave propagation is solved by collocation of boundary points. In order to improve the computation stability, both the incident wave elevation and velocity potential are applied to the wave generation. A sponge damping layer combined with the Sommerfeld radiation condition is used on the radiation boundary. The present model is applied to simulate the propagation of regular and irregular waves. The numerical results are validated by analytical solutions and experimental data and good agreements are observed. Copyright © 2008 John Wiley & Sons, Ltd.     

Application of a perfectly matched layer to the nonlinear wave equation

We derive a perfectly matched layer-like damping layer for the nonlinear wave equation. In the layer, only two auxiliary variables are needed. In the linear case the layer is perfectly matched, but in the nonlinear case it is not. Well posedness is established for the linear case. We also prove various energy estimates which can be used as a starting point for establishing stability of more general cases. In particular, we are able to show estimates for a special type of nonlinearity.Numerical experiments that show the effectiveness of the layer are presented both for nonlinear and linear problems. In the computations, we use an eighth order summation-by-parts discretization in space and implement the boundary conditions using a penalty procedure. We present new stability results for this discretization applied to the second order wave equation in the case with Dirichlet boundary conditions.     

Generation of linear and nonlinear waves in numerical wave tank using clustering technique-volume of fluid method

A two-dimensional(2D) numerical model is developed for the wave simulation and propagation in a wave flume.The fluid flow is assumed to be viscous and incompressible,and the Navier-Stokes and continuity equations are used as the governing equations.The standard k-ε model is used to model the turbulent flow.The NavierStokes equations are discretized using the staggered grid finite difference method and solved by the simplified marker and cell(SMAC) method.Waves are generated and propagated using a piston type wave maker.An open boundary condition is used at the end of the numerical flume.Some standard tests,such as the lid-driven cavity,the constant unidirectional velocity field,the shearing flow,and the dam-break on the dry bed,are performed to valid the model.To demonstrate the capability and accuracy of the present method,the results of generated waves are compared with available wave theories.Finally,the clustering technique(CT) is used for the mesh generation,and the best condition is suggested.     

Finite element method for shallow water equation including open boundary condition

A method to deal with an open boundary condition in the analysis of water surface waves, the tide, etc. by means of the finite element method is proposed in this paper. The present method has two important features relating to the treatment of the open boundary condition. The first feature is to consider the non-reflective virtual boundary condition which has been developed in the numerical wave analysis method. The incident wave conditions without spurious reflected waves can be imposed at the open boundary. The second feature is to identify the amplitude of the components of incident waves in terms of observed water elevations in the field of standing waves. This can be done as a parameter identification based on an optimization problem by applying the conjugate gradient method. The applicability of this method to wave propagation problems is verified by several numerical computations.     

Elastic wave scattering by a three-dimensional inhomogeneity in an elastic half space

In this paper we will consider scattering of elastic waves in a half space. The half space is an isotropic, linear and homogeneous medium except for a finite inhomogeneity. The T-matrix method (also called the “extended boundary condition method” or “null field approach”) is extended to derive expressions for the elastic field inside the half space and the surface field on the interface. The assumptions on the source that excites the half space are fairly weak. In the numerical applications found in this paper we assume a Rayleigh surface wave to be the incoming field, and we only compute the surface displacements. We make illustrations on some simple types of scatterers (spheres and spheroids; the latter ones can be arbitrarily oriented).     

Three dimensional shock wave/boundary layer interactions

An investigation into a three-dimensional, curved shock wave interacting with a three-dimensional, curved boundary layer on a slender body is presented. Three different nose profiles mounted on a cylindrical body were tested in a supersonic wind tunnel and numerically simulated by solving the Navier–Stokes equations. The conical and hemispherical nose profiles tested were found to generate shock waves of sufficient strength to separate the boundary layer on the cylinder, while the shock wave generated by the ogival profile did not separate the boundary layer. For the separated flow, separation was found to occur predominantly on the windward side of the cylinder with the lee-side remaining shielded from the direct impact of the incident shock wave. A thickening of the boundary layer on the lee-side of all the profiles was observed, and in the conical and hemispherical cases this leads to the re-formation of the incident shock wave some distance away from the surface of the cylinder. A complex reflection pattern off the shock wave/boundary layer interaction (SWBLI) was also identified for the separated flow cases. For comparative purposes, an inviscid simulation was performed using the hemispherical profile. Significant differences between the viscous and inviscid results were noted including the absence of a boundary layer leading to a simplified shock wave reflection pattern forming. The behaviour of the incident shock wave on the lee-side of the cylinder was also affected with the shock wave amalgamating on the surface of the cylinder instead of away from the surface as per the viscous case. Test data from the wind tunnel identified two separation lines present on the cylindrical surface of the hemispherical SWBLI generator. The pair of lines were not explicitly evident in the original CFD simulations run, but were later identified in a high-resolution simulation.     

Nonlinear development of a wave in a boundary layer

In recent years definite progress has been achieved in the construction of theoretical models of nonlinear wave processes which lead to a transition from laminar to turbulent flow [1, 2]. At the same time, there is a shortage of actual experimental material, especially for flows in a boundary layer. Fairly thorough experimental studies have been carried out only on the initial stage of the development of disturbances in a boundary layer, which is satisfactorily describable by the linear theory of hydrodynamic stability. In evaluating the theoretical models of subsequent stages of the transition, investigators have been forced to turn chiefly to much earlier experiments carried out by the United States National Bureau of Standards [3, 4], in which the main attention was concentrated on the three-dimensional structure of the transition region. The present investigation was undertaken for the purpose of obtaining detailed data on the structure of the flow in the transition region when there is disturbance in the laminar boundary layer of a two-dimensional wave. In order to make the two-dimensional nonlinear effects stand out more clearly, the amplitude of the wave was specified to be fairly large from the very outset. In contrast to earlier investigations, the main attention was centered on the study of the spectral composition of the disturbance field.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 49–58, May–June, 1977.     

Solving a fully nonlinear highly dispersive Boussinesq model with mesh‐less least square‐based finite difference method

Combining mesh‐less finite difference method and least square approximation, a new numerical model is developed for water wave propagation model in two horizontal dimensions. In the numerical formulation of the method, the approximation of the unknown functions and their derivatives are constructed on a set of nodes in a local circular‐shaped region. The Boussinesq equations studied in this paper is a fully nonlinear and highly dispersive model, which is composed of the exact boundary conditions and the truncated series expansion solution of the Laplace equation. The resultant system involves a sparse, unsymmetrical matrix to be solved at each time step of the simulation. Matrix solutions are studied to reduce the computing resource requirements and improve the efficiency and accuracy. The convergence properties of the present numerical method are investigated. Preliminary verifications are given for nonlinear wave shoaling problems; the numerical results agree well with experimental data available in the literature. Copyright © 2006 John Wiley & Sons, Ltd.     

Simulation of supersonic turbulent flow in the vicinity of an inclined backward-facing step

Large eddy simulation (LES) is a viable and powerful tool to analyse unsteady three-dimensional turbulent flows. In this article, the method of LES is used to compute a plane turbulent supersonic boundary layer subjected to different pressure gradients. The pressure gradients are generated by allowing the flow to pass in the vicinity of an expansion–compression ramp (inclined backward-facing step with leeward-face angle of 25°) for an upstream Mach number of 2.9. The inflow boundary condition is the main problem for all turbulent wall-bounded flows. An approach to solve this problem is to extract instantaneous velocities, temperature and density data from an auxiliary simulation (inflow generator). To generate an appropriate realistic inflow condition to the inflow generator itself the rescaling technique for compressible flows is used. In this method, Morkovin's hypothesis, in which the total temperature fluctuations are neglected compared with the static temperature fluctuations, is applied to rescale and generate the temperature profile at inlet. This technique was successfully developed and applied by the present author for an LES of subsonic three-dimensional boundary layer of a smooth curved ramp. The present LES results are compared with the available experimental data as well as numerical data. The positive impact of the rescaling formulation of the temperature is proven by the convincing agreement of the obtained results with the experimental data compared with published numerical work and sheds light on the quality of the developed compressible inflow generator.     

Boundary element method for wave trapping by a multi-layered trapezoidal breakwater near a sloping rigid wall

This study examines the multiple layers in a rubble mound breakwater and their effect on reflection and dissipation of incoming ocean waves. The numerical model is developed using multi-domain boundary element method for oblique water wave trapping near a sloping wall by a multi-layered trapezoidal porous structure, which is utilized to model armour, filter and core layers while examining the hydrodynamics in different configurations. Both, the constant element and linear element approaches to boundary element method are discussed. The cases of bottom-standing porous structures as being submerged and fully extended are considered. The wave hydrodynamics over the structure is described by the reflection and dissipation coefficients along with the forces acting on the sloping wall, and is influenced by wave and structural parametrics of the system. The influence of armour layer in different configurations is highlighted for various structural and wave parameters.

     

Global volume conservation in unsteady free surface flows with energy absorbing far‐end boundaries

A wave absorption filter for the far‐end boundary of semi‐infinite large reservoirs is developed for numerical simulation of unsteady free surface flows. Mathematical model is based on finite volume solution of the Navier–Stokes equations and depth‐integrated continuity equation to track the free surface. The Sommerfeld boundary condition is applied at the far‐end of the truncated computational domain. A dissipation zone is formed by applying artificial pressure on water surface to dissipate the kinetic energy of the outgoing waves. The computational scheme is tested to verify the conservation of total fluid volume in the domain for long simulation durations. Combination of the Sommerfeld boundary and dissipation zone can effectively minimize reflections and prevent cumulative changes in total fluid volume in the domain. Solitary wave, nonlinear periodic waves and irregular waves are simulated to illustrate the numerical developments. Earthquake excited surface waves and nonlinear hydrodynamic pressures in a dam–reservoir are computed. Copyright © 2009 John Wiley & Sons, Ltd.     

Numerical simulation of fully nonlinear irregular wave tank in three dimension

A fully nonlinear irregular wave tank has been developed using a three‐dimensional higher‐order boundary element method (HOBEM) in the time domain. The Laplace equation is solved at each time step by an integral equation method. Based on image theory, a new Green function is applied in the whole fluid domain so that only the incident surface and free surface are discretized for the integral equation. The fully nonlinear free surface boundary conditions are integrated with time to update the wave profile and boundary values on it by a semi‐mixed Eulerian–Lagrangian time marching scheme. The incident waves are generated by feeding analytic forms on the input boundary and a ramp function is introduced at the start of simulation to avoid the initial transient disturbance. The outgoing waves are sufficiently dissipated by using a spatially varying artificial damping on the free surface before they reach the downstream boundary. Numerous numerical simulations of linear and nonlinear waves are performed and the simulated results are compared with the theoretical input waves. Copyright © 2006 John Wiley & Sons, Ltd.     

Dispersion and attenuation of torsional wave in a viscoelastic layer bonded between a layer and a half-space of dry sandy media

Propagation of a torsional wave in a doubly-layered half-space structure of an initially stressed heterogeneous viscoelastic layer sandwiched between a layer and a half-space of heterogeneous dry sandy media is studied. A closed form complex expression for the velocity profile is obtained under effective boundary conditions. The real part of the complex expression provides a dispersion equation, and the imaginary part yields a damping equation. The derived dispersion and damped equations are in well agreement with the classical Love wave condition. In addition, to study the effect of the dissipation factor, the attenuation coefficient, the sandy parameters, the initial stress, the heterogeneity parameters, and the thickness ratio parameter, some noteworthy contemplations are made by numerical calculations and graphical visuals. The results of this paper may present a deeper insight into the behaviour of propagation phenomena in heterogeneous viscoelastic and heterogeneous dry sandy materials that can provide a theoretical guide for the design and optimization in the field of earthquake engineering. The study also reveals that the presence of a damping part due to viscoelasticity affects the torsional wave propagation significantly.     

Linear water wave propagation through multiple floating elastic plates of variable properties

We investigate the problem of linear water wave propagation under a set of elastic plates of variable properties. The problem is two-dimensional, but we allow the waves to be incident from an angle. Since the properties of the elastic plates can be set arbitrarily, the solution method can also be applied to model regions of open water as well as elastic plates. We assume that the boundary conditions at the plate edges are the free boundary conditions, although the method could be extended straightforwardly to cover other possible boundary conditions. The solution method is based on an eigenfunction expansion under each elastic plate and on matching these expansions at each plate boundary. We choose the number of matching conditions so that we have fewer equations than unknowns. The extra equations are found by applying the free-edge boundary conditions. We show that our results agree with previous work and that they satisfy the energy balance condition. We also compare our results with a series of experiments using floating elastic plates, which were performed in a two-dimensional wave tank.