A Wave Superposition-Finite Element Method for Calculating the Radiated Noise Generated by Volumetric Targets in Shallow Water

A combined method of wave superposition and finite element is proposed to solve the radiation noise of targets in shallow sea.Taking the sound propagation of spherical sound source in shallow sea as an example,the radiation sound field of the spherical sound source is equivalent to the linear superposition of the radiation sound field of several internal point sound sources,and then the radiated noise induced by spherical sound source can be predicted quickly.The accuracy and efficiency of the method are verified by comparing with the numerical results of finite element method,and the rapid prediction of underwater radiated noise of cylindrical shell is carried out based on the method.The results show that compared with the finite element method,the relative error of the calculation results under different simulation conditions does not exceed 0.1%,and the calculation time is about 1/10 of the finite element method,so this method can be used to solve the radiated noise of shallow underwater targets.

NUMERICAL INVESTIGATION OF TOROIDAL SHOCK WAVES FOCUSING USING DISCONTINUOUS GALERKIN FINITE ELEMENT METHOD

A numerical simulation of the toroidal shock wave focusing in a co-axial cylindrical shock tube is inves- tigated by using discontinuous Galerkin （DG） finite element method to solve the axisymmetric Euler equations. For validating the numerical method, the shock-tube problem with exact solution is computed, and the computed results agree well with the exact cases. Then, several cases with higher incident Mach numbers varying from 2.0 to 5.0 are simulated. Simulation results show that complicated flow-field structures of toroidal shock wave diffraction, reflection, and focusing in a co-axial cylindrical shock tube can be obtained at different incident Mach numbers and the numerical solutions appear steep gradients near the focusing point, which illustrates the DG method has higher accuracy and better resolution near the discontinuous point. Moreover, the focusing peak pres- sure with different grid scales is compared.

PARTITION OF UNITY FINITE ELEMENT METHOD FOR SHORT WAVE PROPAGATION IN SOLIDS

A partition of unity finite element method for numerical simulation of short wave propagation in solids is presented. The finite element spaces were constructed by multiplying the standard isoparametric finite element shape functions, which form a partition of unity, with the local subspaces defined on the corresponding shape functions, which include a priori knowledge about the wave motion equation in trial spaces and approximately reproduce the highly oscillatory properties within a single element. Numerical examples demonstrate the performance of the proposed partition of unity finite element in both computational accuracy and efficiency.

IN-PLANE WAVE MOTION IN FINITE ELEMENT MODEL

The analysis method of lattice dynamics in classical physics is extended to study the properties of in-plane wave motion in the hybrid-mass finite element model in this paper. The dispersion equations of P and SV waves in the discrete model are first obtained by means of separating the characteristic equation of the motion equation, and then used to analyse the properties of P-and SV-homogeneous, inhomogeneous waves and other types of motion in the model. The dispersion characters, cut-off frequencies of P and SV waves, the polarization drift and appendent anisotropic property of wave motion caused by the discretization are finally discussed.

Nonlinear finite element analysis of effect of seismic waves on dynamic response of Shiziping dam

Many high earth-rockfill dams are constructed in the west of China. The seismic intensity at the dam site is usually very high, thus it is of great importance to ensure the safety of the dam in meizoseismal area. A 3D FEM model is established to analyze the seismic responses of Shiziping earth-rockfill dam. The nonlinear elastic Duncan-Chang constitutive model and the equivalent viscoelastic constitutive model are used to simulate the static and dynamic stress strain relationships of the dam materials, respectively. Four groups of seismic waves are inputted from the top of the bedrock to analyze the dynamic responses of the dam. The numerical results show that the calculated dynamic magnification factors display a good consistency with the specification values. The site spectrum results in larger acceleration response than the specification spectrum. The analysis of relative dynamic displacement indicates that the displacement at the downstream side of the dam is larger than that at the upstream side. The displacement response reduces from the center of river valley to two banks. The displacement responses corresponding to the specification spectrum are a little smaller than those corresponding to the site spectrum. The analysis of shear stress indicates that a large shear stress area appears in the upstream overburden layer, where the shear stress caused by site waves is larger than that caused by specification waves. The analysis of dynamic principal stress indicates that the minimum dynamic stresses in corridor caused by specification and site waves have little difference. The maximum and minimum dynamic stresses are relatively large at two sides. The largest tensile stress occurs at two sides of the floor of grouting corridor, which may result in the crack near the corridor side. The numerical results present good consistency with the observation data of the grouting corridor in Wenchuan earthquake.

Propagations of Rayleigh and Love waves in ZnO films/glass substrates analyzed by three-dimensional finite element method

Propagation characteristics of surface acoustic waves（SAWs） in ZnO films/glass substrates are theoretically investigated by the three-dimensional（3D） finite element method. At first, for（11ˉ20） ZnO films/glass substrates, the simulation results confirm that the Rayleigh waves along the [0001] direction and Love waves along the [1ˉ100] direction are successfully excited in the multilayered structures. Next, the crystal orientations of the ZnO films are rotated, and the influences of ZnO films with different crystal orientations on SAW characterizations, including the phase velocity, electromechanical coupling coefficient, and temperature coefficient of frequency, are investigated. The results show that at appropriate h/λ, Rayleigh wave has a maximum k^2 of 2.4% in（90°, 56.5°, 0°） ZnO film/glass substrate structure; Love wave has a maximum k^2 of 3.81% in（56°, 90°, 0°） ZnO film/glass substrate structure. Meantime, for Rayleigh wave and Love wave devices, zero temperature coefficient of frequency（TCF） can be achieved at appropriate ratio of film thickness to SAW wavelength. These results show that SAW devices with higher k^2 or lower TCF can be fabricated by flexibly selecting the crystal orientations of ZnO films on glass substrates.

Scaled Boundary Finite Element Analysis of Wave Passing A Submerged Breakwater

The scaled boundary finite element method （SBFEM） is a novel semi-analytical technique combining the advantage of the finite element method （FEM） and the boundary element method （BEM） with its unique properties. In this paper, the SBFEM is used for computing wave passing submerged breakwaters, and the reflection coeffcient and transmission coefficient are given for the case of wave passing by a rectangular submerged breakwater, a rigid submerged barrier breakwater and a trapezium submerged breakwater in a constant water depth. The results are compared with the analytical solution and experimental results. Good agreement is obtained. Through comparison with the results using the dual boundary element method （DBEM）, it is found that the SBFEM can obtain higher accuracy with fewer elements. Many submerged breakwaters with different dimensions are computed by the SBFEM, and the changing character of the reflection coeffcient and the transmission coefficient are given in the current study.

A Finite Element Solution of Wave Forces on Submerged Horizontal Circular Cylinders

In this paper, a numerical model is established for estimating the wave forces on a submerged horizontal circular cylinder. For predicting the wave motion, a set of two dimensional Navier Stokes equations is solved numerically with a finite element method. In order to track the moving non linear wave surface boundary, the Navier Stokes equations are discretized in a moving mesh system. After each computational time step, the mesh is modified according to the changed wave surface boundary. In order to stabilize the numerical procedure, a three step finite element method is applied in the time integration. The water sloshing in a tank and wave propagation over a submerged bar are simulated for the first time to validate the present model. The computational results agree well with the analytical solution and the experimental data. Finally, the model is applied to the simulation of interaction between waves and a submerged horizontal circular cylinder. The effects of the KC number and the cylinder depth on the wave forces are studied.

Dynamic Adaptive Finite Element Analysis of Acoustic Wave Propagation Due to Underwater Explosion for Fluid-structure Interaction Problems

In this paper, an investigation into the propagation of far field explosion waves in water and their effects on nearby structures are carried out. For the far field structure, the motion of the fluid surrounding the structure may be assumed small, allowing linearization of the governing fluid equations. A complete analysis of the problem must involve simultaneous solution of the dynamic response of the structure and the propagation of explosion wave in the surrounding fluid. In this study, a dynamic adaptive finite element procedure is proposed. Its application to the solution of a 2D fluid-structure interaction is investigated in the time domain. The research includes：a） calculation of the far-field scatter wave due to underwater explosion including solution of the time-depended acoustic wave equation, b） fluid-structure interaction analysis using coupled Euler-Lagrangian approach, and c） adaptive finite element procedures employing error estimates, and re-meshing. The temporal mesh adaptation is achieved by local regeneration of the grid using a time-dependent error indicator based on curvature of pressure function. As a result, the overall response is better predicted by a moving mesh than an equivalent uniform mesh. In addition, the cost of computation for large problems is reduced while the accuracy is improved.

PARTITION OF UNITY FINITE ELEMENT METHOD FOR SHORT WAVE PROPAGATION IN SOLIDS

A partition of unity finite element method for numerical simulation of short wave propagation in solids is presented. The finite element spaces were constructed by multiplying the standard isoparametric finite element shape functions, which form a partition of unity, with the local subspaces defined on the corresponding shape functions, which include a priori knowledge about the wave motion equation in trial spaces and approximately reproduce the highly oscillatory properties within a single element. Numerical examples demonstrate the performance of the proposed partition of unity finite element in both computational accuracy and efficiency.

FINITE ELEMENT ANALYSIS OF WAVE PROPAGATION IN FLUID-SATURATED POROUS MEDIA

With the porous media model based on mixture theory, a finite element formulation for dynamic transient analysis of fluid_saturated two_phase porous media is presented. Time integration of the equation, deduced with penalty method, can be performed by using implicit or explicit method. One_dimensional wave propagation in column under step loading and impulsive loading are analyzed with the developed finite element program. The obtained curves of displacements, velocities, effective stresses and pore pressures against time demonstrate the existence of wave propagation phenomena, which coincide with the theoretical results.

Second-generation wavelet finite element based on the lifting scheme for GPR simulation

Ground-penetrating radar(GPR)is a highly efficient,fast and non-destructive exploration method for shallow surfaces.High-precision numerical simulation method is employed to improve the interpretation precision of detection.Second-generation wavelet finite element is introduced into the forward modeling of the GPR.As the finite element basis function,the second-generation wavelet scaling function constructed by the scheme is characterized as having multiple scales and resolutions.The function can change the analytical scale arbitrarily according to actual needs.We can adopt a small analysis scale at a large gradient to improve the precision of analysis while adopting a large analytical scale at a small gradient to improve the efficiency of analysis.This approach is beneficial to capture the local mutation characteristics of the solution and improve the resolution without changing mesh subdivision to realize the efficient solution of the forward GPR problem.The algorithm is applied to the numerical simulation of line current radiation source and tunnel non-dense lining model with analytical solutions.Result show that the solution results of the secondgeneration wavelet finite element are in agreement with the analytical solutions and the conventional finite element solutions,thereby verifying the accuracy of the second-generation wavelet finite element algorithm.Furthermore,the second-generation wavelet finite element algorithm can change the analysis scale arbitrarily according to the actual problem without subdividing grids again.The adaptive algorithm is superior to traditional scheme in grid refinement and basis function order increase,which makes this algorithm suitable for solving complex GPR forward-modeling problems with large gradient and singularity.

Finite Elements Based on Deslauriers-Dubuc Wavelets for Wave Propagation Problems

This paper presents the formulation of finite elements based on Deslauriers-Dubuc interpolating scaling functions, also known as Interpolets, for their use in wave propagation modeling. Unlike other wavelet families like Daubechies, Interpolets possess rational filter coefficients, are smooth, symmetric and therefore more suitable for use in numerical methods. Expressions for stiffness and mass matrices are developed based on connection coefficients, which are inner products of basis functions and their derivatives. An example in 1-D was formulated using Central Difference and Newmark schemes for time differentiation. Encouraging results were obtained even for large time steps. Results obtained in 2-D are compared with the standard Finite Difference Method for validation.

Finite element analysis of controlling the TC4 thin plate weldment wave-like deformation by welding with impacting rotation

The new technology of welding with impacting rotation is put forward to decrease the wave-like deformation of the TC4 thin plate weldment. The thermal stress and strain are vital to understand the mechanism of controlling the wave-like deformation. In order to know the development of internal thermal stress and strain, finite element method is utilized for- the stress and strain are difficult to be investigated by experimental methods during the welding process. Temperature field, thermal stress evolution and distortion of thin plate are compared with the test results such as weld thermal cycle, residual stress sectioning measurement, and the deflection of the thin plate respectively. By the finite element analysis and test results verification, the meehaaism of the technology to control the wave-like deformation is brought forward, non-uniform thermal elastic strain between compressive plastic region and elastic extensive region is diminished by a certain amount of extensive plastic deformation by welding with impacting rotation process.

GENERALIZED WAVE EQUATION FINITE ELEMENT METHOD FOR SOLVING TWO-DIMENSIONAL TIDAL WAVES

The study of tidal circulation has a long history . The numerical simulation of tidal flow has been developed greatly with the development of computer techniques in the past two decades. The generalized wave equation finite-element method is a relatively new numerical model for studying shallow water flow . This method was used to simulate tidal waves of the Gulf of St. Lawrence in Canada . The very good agreement of the numerical results with the field data indicated that the model is an effective and promising numerical method for solving two-dimensional tidal wave problems .

Finite element equations and numerical simulation of elastic wave propagation in two-phase anisotropic media

Based on Biot theory of two-phase anisotropic media and Hamilton theory about dynamic problem,finite element equations of elastic wave propagation in two-phase anisotropic media are derived in this paper.Numerical solution of finite element equations is given.Finally,Properties of elastic wave propagation are observed and analyzed through FEM modeling.

Finite element simulation of steady state SH wave motion

A discrete local transmitting boundary is combined with a lumped-mass finite element technique for simulating steady SH wave motion in a layered and unbounded medium.The combination decouples any node of the finite element model from the others except those in its neighborhood,and an effective algorithm of the Gauss elimination in implementation of the method is then devised to reduce the main memory and computing time.The method and its accuracy are first analysed via simulating steady SH wave motion in layered elastic media,the effective algorithm is then introduced and illustrated by simple examples.

A Brief Summary of Finite Element Method Applications to Nonlinear Wave-structure Interactions

We review recent advances in the finite element method (FEM) simulations of interactions between waves and structures. Our focus is on the potential theory with the fully nonlinear or second-order boundary condition. The present paper has six sections. A review of previous work on interactions between waves and ocean structures is presented in Section one. Section two gives the mathematical formulation. In Section three, the finite element discretization, mesh generation and the finite element linear system solution methods are described. Section four presents numerical methods including time marching schemes, computation of velocity, remeshing and smoothing techniques and numerical radiation conditions. The application of the FEM to the wave-structure interactions are presented in Section five followed by the concluding remarks in Section six.

Adaptive Lagrange finite element methods for high precision vibrations and piezoelectric acoustic wave computations in SMT structures and plates with nano interfaces

This paper discusses the validity of (adaptive) Lagrange generalized plain finite element method (FEM) and plate element method for accurate analysis of acoustic waves in multi-layered piezoelectric structures with tiny interfaces between metal electrodes and surface mounted piezoelectric substrates. We have come to conclusion that the quantitative relationships between the acoustic and electric fields in a piezoelectric structure can be accurately determined through the proposed finite element methods. The higher-order Lagrange FEM proposed for dynamic piezoelectric computation is proved to be very accurate (prescribed relative error 0.02% - 0.04% ) and a great improvement in convergence accuracy over the higher order Mindlin plate element method for piezoelectric structural analysis due to the assumptions and corrections in the plate theories.The converged lagrange finite element methods are compared with the plate element methods and the computedresults are in good agreement with available exact and experimental data. The adaptive Lagrange finite elementmethods and a new FEA computer program developed for macro- and micro-scale analyses are reviewed, and recently extended with great potential to high-precision nano-scale analysis in this paper and the similarities between piezoelectric and seismic wave propagations in layered structures and plates are stressed.

An ADI Finite Volume Element Method for a Viscous Wave Equation with Variable Coefficients

Based on rectangular partition and bilinear interpolation,we construct an alternating-direction implicit(ADI)finite volume element method,which combined the merits of finite volume element method and alternating direction implicit method to solve a viscous wave equation with variable coefficients.This paper presents a general procedure to construct the alternating-direction implicit finite volume element method and gives computational schemes.Optimal error estimate in L2 norm is obtained for the schemes.Compared with the finite volume element method of the same convergence order,our method is more effective in terms of running time with the increasing of the computing scale.Numerical experiments are presented to show the efficiency of our method and numerical results are provided to support our theoretical analysis.