A new method is presented for the computation of two-dimensional periodic progressive surface waves propagating under the combined influence of gravity and surface tension. The nonlinear surface is expressed by Fourie...A new method is presented for the computation of two-dimensional periodic progressive surface waves propagating under the combined influence of gravity and surface tension. The nonlinear surface is expressed by Fourier series with finite number of terms, after the computational domain is transformed into a unit circle. The dynamic boundary equation is used in its exact nonlinear form and the coefficients of Fourier series are found by the Newton-Raphson method successively. This is a neat method, Yielding high prescision with little computational effort.展开更多
Equivalent staggered-grid(ESG) as a new family of schemes has been utilized in seismic modeling,imaging,and inversion.Traditionally,the Taylor series expansion is often applied to calculate finite-difference(FD) coeff...Equivalent staggered-grid(ESG) as a new family of schemes has been utilized in seismic modeling,imaging,and inversion.Traditionally,the Taylor series expansion is often applied to calculate finite-difference(FD) coefficients on spatial derivatives,but the simulation results suffer serious numerical dispersion on a large frequency zone.We develop an optimized equivalent staggered-grid(OESG) FD method that can simultaneously suppress temporal and spatial dispersion for solving the second-order system of the 3 D elastic wave equation.On the one hand,we consider the coupling relations between wave speeds and spatial derivatives in the elastic wave equation and give three sets of FD coefficients with respect to the P-wave,S-wave,and converted-wave(C-wave) terms.On the other hand,a novel plane wave solution for the 3 D elastic wave equation is derived from the matrix decomposition method to construct the time-space dispersion relations.FD coefficients of the OESG method can be acquired by solving the new dispersion equations based on the Newton iteration method.Finally,we construct a new objective function to analyze P-wave,S-wave,and C-wave dispersion concerning frequencies.The dispersion analyses show that the presented method produces less modeling errors than the traditional ESG method.The synthetic examples demonstrate the effectiveness and superiority of the presented method.展开更多
The mixed covolume method for the regularized long wave equation is devel- oped and studied. By introducing a transfer operator γh, which maps the trial function space into the test function space, and combining the ...The mixed covolume method for the regularized long wave equation is devel- oped and studied. By introducing a transfer operator γh, which maps the trial function space into the test function space, and combining the mixed finite element with the finite volume method, the nonlinear and linear Euler fully discrete mixed covolume schemes are constructed, and the existence and uniqueness of the solutions are proved. The optimal error estimates for these schemes are obtained. Finally, a numerical example is provided to examine the efficiency of the proposed schemes.展开更多
Soil shear wave velocity (SWV) is an important parameter in geotechnical engineering. To measure the soil SWV, three methods are generally used in China, including the single-hole method, cross-hole method and the s...Soil shear wave velocity (SWV) is an important parameter in geotechnical engineering. To measure the soil SWV, three methods are generally used in China, including the single-hole method, cross-hole method and the surface-wave technique. An optimized approach based on a correlation function for single-hole SWV measurement is presented in this paper. In this approach, inherent inconsistencies of the artificial methods such as negative velocities, and too-large and too-small velocities, are eliminated from the single-hole method, and the efficiency of data processing is improved. In addition, verification using the cross-hole method of upper measuring points shows that the proposed optimized approach yields high precision in signal processing.展开更多
This Paper improves the existing fifth order Stokes wave theory by using least Square method, and givesthe optimum result in the meaning of minimum error Squares to satisfy the free surface boundary conditions, and th...This Paper improves the existing fifth order Stokes wave theory by using least Square method, and givesthe optimum result in the meaning of minimum error Squares to satisfy the free surface boundary conditions, and thewave profile can be adjusted according to the measured data. This paper also gives a simplified method for derivingthe parameters of the existing fifth order Stokes wave.展开更多
The nonlinear capillary-gravity wave produced by a vertically oscillating plate, in which the contact-angle model is considered, is studied by use of the Boundary Integral Equation Method (BIEM). The present numerica...The nonlinear capillary-gravity wave produced by a vertically oscillating plate, in which the contact-angle model is considered, is studied by use of the Boundary Integral Equation Method (BIEM). The present numerical experiment shows that the code is robust and efficient for modeling the generation and propagation of capillary-gravity waves. It is found that the wave heights of stationary periodic nonlinear waves radiated away from the plate are dependent on the parameters involved in the contact-angle model. The effect of the contact-angle hysteresis and the nonlinearity of capillary-gravity waves on the wave profile is discussed in the paper.展开更多
H1-Galerkin mixed methods are proposed for viscoelasticity wave equation.Depending on the physical quantities of interest,two methods are discussed.The optimal error estimates and the proof of the existence and unique...H1-Galerkin mixed methods are proposed for viscoelasticity wave equation.Depending on the physical quantities of interest,two methods are discussed.The optimal error estimates and the proof of the existence and uniqueness of semidiscrete solutions are derived for problems in one space dimension.And the methods don't require the LBB condition.展开更多
In this paper, the dispersive long wave equation is studied by Lie symmetry group theory. Firstly, the Lie symmetries of this system are calculated. Secondly, one dimensional optimal systems of Lie algebra and all the...In this paper, the dispersive long wave equation is studied by Lie symmetry group theory. Firstly, the Lie symmetries of this system are calculated. Secondly, one dimensional optimal systems of Lie algebra and all the symmetry reductions are obtained. Finally, based on the power series method and the extended Tanh function method, some new explicit solutions of this system are constructed.展开更多
Consider the inverse diffraction problem to determine a two-dimensional periodic structure from scattered elastic waves measured above the structure.We formulate the inverse problem as a least squares optimization pro...Consider the inverse diffraction problem to determine a two-dimensional periodic structure from scattered elastic waves measured above the structure.We formulate the inverse problem as a least squares optimization problem,following the two-step algorithm by G.Bruckner and J.Elschner[Inverse Probl.,19(2003),315–329]for electromagnetic diffraction gratings.Such a method is based on the Kirsch-Kress optimization scheme and consists of two parts:a linear severely ill-posed problem and a nonlinear well-posed one.We apply this method to both smooth(C2)and piecewise linear gratings for the Dirichlet boundary value problem of the Navier equation.Numerical reconstructions from exact and noisy data illustrate the feasibility of the method.展开更多
Structure health monitoring based on diagnostic Lamb waves has been found to be one of the most promising techniques recently. This paper has a brief review of the new developments on this method including the basic n...Structure health monitoring based on diagnostic Lamb waves has been found to be one of the most promising techniques recently. This paper has a brief review of the new developments on this method including the basic novel of the method, fundamentals and mathematics of Lamb wave propagation, narrowband and wideband Lamb wave excitation methods, optimization of excitation factors and diagnostic Lamb wave interpretation methods.展开更多
The hull form optimization concerns one of the most important applications of wave making resistance theories. In order to obtain a hull form with the minimum wave making resistance, an optimization design method base...The hull form optimization concerns one of the most important applications of wave making resistance theories. In order to obtain a hull form with the minimum wave making resistance, an optimization design method based on the CFD is proposed, which combines the Rankine source method with the nonlinear programming (NLP). The bow-body shape is optimized with the minimum wave making resistance as the objective function. A hull form modification function is introduced to represent an improved hull surface, which can be used to generate a new smooth hull surface by multiplying it by the offset data of the original hull surface. The parameters of the hull form modification function are taken as the design variables. Other constraint conditions can also be considered, for example, in optimizing the lines of the bow, appropriate displacements can be taken as the basic constraints. S60 hull form is selected as the original hull. Three improved hulls are obtained by optimal design. Rankine source method proves to be an effective method in ship form optimization based on analysis of the resistance performance and lines of the improved hull.展开更多
The ship hull surface optimization based on the wave resistance is an important issue in the ship engineering industry. The wavelet method may provide a convenient tool for the surface hull optimization. As a prelimin...The ship hull surface optimization based on the wave resistance is an important issue in the ship engineering industry. The wavelet method may provide a convenient tool for the surface hull optimization. As a preliminary study, we use the wavelet method to optimize the hull surface based on the Michel wave resistance for a Wigley model in this paper. Firstly, we express the model's surface by the wavelet decomposition expressions and obtain a reconstructed surface and then validate its accuracy. Secondly, we rewrite the Michel wave resistance formula in the wavelet bases, resulting in a simple formula containing only the ship hull surface's wavelet coefficients. Thirdly, we take these wavelet coefficients as optimization variables, and analyze the main wave resistance distribution in terms of scales and locations, to reduce the number of optimization variables. Finally, we obtain the optimal hull surface of the Wigley model through genetic algorithms, reducing the wave resistance almost by a half. It is shown that the wavelet method may provide a new approach for the hull optimization.展开更多
The hull form optimization concerns one of the most important applications of wave-making resistance theories.In recent years,scholars can determine the hull form by using the optimization method based on the computat...The hull form optimization concerns one of the most important applications of wave-making resistance theories.In recent years,scholars can determine the hull form by using the optimization method based on the computational fluid dynamics(CFD)and other mathematical techniques.In this paper,the hull form optimization method based on the Rankine source method and nonlinear programming(NLP)is discussed;in the optimization process,a hull form modification function is introduced to represent an improved hull surface and to generate a new smooth hull surface by changing its frame lines and bow stem profiles under the prescribed design constraints. Numerical example is given for a practical container hull form.Finally,shape optimization of bow bulls is shown for non-protruding and protruding bow bulls.This study presents a simplified and practical design method to the select frame lines of bow bulls.展开更多
The hydrodynamic shape of the heaving buoy is an important factor of the motion response in waves and thus concerns the energy conversion efficiency for the point absorbers(PAs).The current experience-based designs ar...The hydrodynamic shape of the heaving buoy is an important factor of the motion response in waves and thus concerns the energy conversion efficiency for the point absorbers(PAs).The current experience-based designs are time consuming and not very efficient,hence,faster and smarter methods are desirable.An automated optimization method based on a fully parametric modeling method and computational fluid dynamics(CFD),is proposed in this paper.Using this method,a benchmark buoy is screen designed and then optimized by maximizing the heave motion response.The geometry is described parametrically and deformed by means of the free-form deformation(FFD)method.During the optimization process,the expansion factor of control points is the basis for the variations.A combination of the Sobol and the non-dominated sorting genetic algorithm II(NSGA-II)is used to search for the solutions.After several iterations,the heaving buoy shape with optimal heave motion response is obtained.The analyses show that the heave motion response has increased 55.3%after optimization.The developed methodology is valid and seems to be a promising way to design a novel buoy that can significantly improve the wave energy conversion efficiency of the PAs in future.展开更多
Weak Galerkin finite element method is introduced for solving wave equation with interface on weak Galerkin finite element space(Pk(K),P_(k−1)(∂K),[P_(k−1)(K)]^(2)).Optimal order a priori error estimates for both spac...Weak Galerkin finite element method is introduced for solving wave equation with interface on weak Galerkin finite element space(Pk(K),P_(k−1)(∂K),[P_(k−1)(K)]^(2)).Optimal order a priori error estimates for both space-discrete scheme and implicit fully discrete scheme are derived in L1(L2)norm.This method uses totally discontinuous functions in approximation space and allows the usage of finite element partitions consisting of general polygonal meshes.Finite element algorithm presented here can contribute to a variety of hyperbolic problems where physical domain consists of heterogeneous media.展开更多
In order to reduce transmission loss of the optical waveguide in Mach-Zehnder (M-Z) electro-optical (EO) polymer modulator,the basic iterative formula of semi-vector finite-difference beam propagation method (FD-BPM) ...In order to reduce transmission loss of the optical waveguide in Mach-Zehnder (M-Z) electro-optical (EO) polymer modulator,the basic iterative formula of semi-vector finite-difference beam propagation method (FD-BPM) is obtained from the scalar wave equation. The transition waveguide is combined with S-type bend branch waveguide for the M-Z EO modulator in the branch waveguide. The effects of structure parameters such as ridge width,length of the branch waveguide and interferometer spacing on the transmission loss are systematically studied by using the semi-vector FD-BPM method. The structure is optimized as an S-sine bend branch waveguide,with rib width w=7μm,length of branch waveguide L=1200μm and interferometer spacing G=22 μm. The results show that the optimized structure can reduce transmission loss to 0.083 dB,which have a certain reference value to the design of optical waveguide in M-Z polymer modulator.展开更多
文摘A new method is presented for the computation of two-dimensional periodic progressive surface waves propagating under the combined influence of gravity and surface tension. The nonlinear surface is expressed by Fourier series with finite number of terms, after the computational domain is transformed into a unit circle. The dynamic boundary equation is used in its exact nonlinear form and the coefficients of Fourier series are found by the Newton-Raphson method successively. This is a neat method, Yielding high prescision with little computational effort.
文摘Equivalent staggered-grid(ESG) as a new family of schemes has been utilized in seismic modeling,imaging,and inversion.Traditionally,the Taylor series expansion is often applied to calculate finite-difference(FD) coefficients on spatial derivatives,but the simulation results suffer serious numerical dispersion on a large frequency zone.We develop an optimized equivalent staggered-grid(OESG) FD method that can simultaneously suppress temporal and spatial dispersion for solving the second-order system of the 3 D elastic wave equation.On the one hand,we consider the coupling relations between wave speeds and spatial derivatives in the elastic wave equation and give three sets of FD coefficients with respect to the P-wave,S-wave,and converted-wave(C-wave) terms.On the other hand,a novel plane wave solution for the 3 D elastic wave equation is derived from the matrix decomposition method to construct the time-space dispersion relations.FD coefficients of the OESG method can be acquired by solving the new dispersion equations based on the Newton iteration method.Finally,we construct a new objective function to analyze P-wave,S-wave,and C-wave dispersion concerning frequencies.The dispersion analyses show that the presented method produces less modeling errors than the traditional ESG method.The synthetic examples demonstrate the effectiveness and superiority of the presented method.
基金supported by the National Natural Science Fundation of China (No. 11061021)the Science Research of Inner Mongolia Advanced Education (Nos. NJ10006, NJ10016, and NJZZ12011)the National Science Foundation of Inner Mongolia (Nos. 2011BS0102 and 2012MS0106)
文摘The mixed covolume method for the regularized long wave equation is devel- oped and studied. By introducing a transfer operator γh, which maps the trial function space into the test function space, and combining the mixed finite element with the finite volume method, the nonlinear and linear Euler fully discrete mixed covolume schemes are constructed, and the existence and uniqueness of the solutions are proved. The optimal error estimates for these schemes are obtained. Finally, a numerical example is provided to examine the efficiency of the proposed schemes.
基金Pre-research Project of Yantai Unverity Under Project No. TM05B35Shandong Natural Science Foundation Under Project No. bs08003 Key Foundation of Ministry of Education Under Project No. 207062
文摘Soil shear wave velocity (SWV) is an important parameter in geotechnical engineering. To measure the soil SWV, three methods are generally used in China, including the single-hole method, cross-hole method and the surface-wave technique. An optimized approach based on a correlation function for single-hole SWV measurement is presented in this paper. In this approach, inherent inconsistencies of the artificial methods such as negative velocities, and too-large and too-small velocities, are eliminated from the single-hole method, and the efficiency of data processing is improved. In addition, verification using the cross-hole method of upper measuring points shows that the proposed optimized approach yields high precision in signal processing.
文摘This Paper improves the existing fifth order Stokes wave theory by using least Square method, and givesthe optimum result in the meaning of minimum error Squares to satisfy the free surface boundary conditions, and thewave profile can be adjusted according to the measured data. This paper also gives a simplified method for derivingthe parameters of the existing fifth order Stokes wave.
文摘The nonlinear capillary-gravity wave produced by a vertically oscillating plate, in which the contact-angle model is considered, is studied by use of the Boundary Integral Equation Method (BIEM). The present numerical experiment shows that the code is robust and efficient for modeling the generation and propagation of capillary-gravity waves. It is found that the wave heights of stationary periodic nonlinear waves radiated away from the plate are dependent on the parameters involved in the contact-angle model. The effect of the contact-angle hysteresis and the nonlinearity of capillary-gravity waves on the wave profile is discussed in the paper.
基金Supported by NNSF(10601022,11061021)Supported by NSF of Inner Mongolia Au-tonomous Region(200607010106)Supported by SRP of Higher Schools of Inner Mongolia(NJ10006)
文摘H1-Galerkin mixed methods are proposed for viscoelasticity wave equation.Depending on the physical quantities of interest,two methods are discussed.The optimal error estimates and the proof of the existence and uniqueness of semidiscrete solutions are derived for problems in one space dimension.And the methods don't require the LBB condition.
文摘In this paper, the dispersive long wave equation is studied by Lie symmetry group theory. Firstly, the Lie symmetries of this system are calculated. Secondly, one dimensional optimal systems of Lie algebra and all the symmetry reductions are obtained. Finally, based on the power series method and the extended Tanh function method, some new explicit solutions of this system are constructed.
基金the support by the German Research Foundation(DFG)under Grant No.EL 584/1-2.
文摘Consider the inverse diffraction problem to determine a two-dimensional periodic structure from scattered elastic waves measured above the structure.We formulate the inverse problem as a least squares optimization problem,following the two-step algorithm by G.Bruckner and J.Elschner[Inverse Probl.,19(2003),315–329]for electromagnetic diffraction gratings.Such a method is based on the Kirsch-Kress optimization scheme and consists of two parts:a linear severely ill-posed problem and a nonlinear well-posed one.We apply this method to both smooth(C2)and piecewise linear gratings for the Dirichlet boundary value problem of the Navier equation.Numerical reconstructions from exact and noisy data illustrate the feasibility of the method.
基金The authors acknowledge the financial supports from the National Natural Science Foundation of China under grant No.90305005,50135030
文摘Structure health monitoring based on diagnostic Lamb waves has been found to be one of the most promising techniques recently. This paper has a brief review of the new developments on this method including the basic novel of the method, fundamentals and mathematics of Lamb wave propagation, narrowband and wideband Lamb wave excitation methods, optimization of excitation factors and diagnostic Lamb wave interpretation methods.
文摘The hull form optimization concerns one of the most important applications of wave making resistance theories. In order to obtain a hull form with the minimum wave making resistance, an optimization design method based on the CFD is proposed, which combines the Rankine source method with the nonlinear programming (NLP). The bow-body shape is optimized with the minimum wave making resistance as the objective function. A hull form modification function is introduced to represent an improved hull surface, which can be used to generate a new smooth hull surface by multiplying it by the offset data of the original hull surface. The parameters of the hull form modification function are taken as the design variables. Other constraint conditions can also be considered, for example, in optimizing the lines of the bow, appropriate displacements can be taken as the basic constraints. S60 hull form is selected as the original hull. Three improved hulls are obtained by optimal design. Rankine source method proves to be an effective method in ship form optimization based on analysis of the resistance performance and lines of the improved hull.
基金Project supported by the Natural National Science Foundation of China(Grant Nos.51309040,51379033)the National Key Basic Research Development Program of China(973 Program,Grant No.2013CB036101)the Fundamental research fund for the Central Universities(Grant No.DMU3132015089)
文摘The ship hull surface optimization based on the wave resistance is an important issue in the ship engineering industry. The wavelet method may provide a convenient tool for the surface hull optimization. As a preliminary study, we use the wavelet method to optimize the hull surface based on the Michel wave resistance for a Wigley model in this paper. Firstly, we express the model's surface by the wavelet decomposition expressions and obtain a reconstructed surface and then validate its accuracy. Secondly, we rewrite the Michel wave resistance formula in the wavelet bases, resulting in a simple formula containing only the ship hull surface's wavelet coefficients. Thirdly, we take these wavelet coefficients as optimization variables, and analyze the main wave resistance distribution in terms of scales and locations, to reduce the number of optimization variables. Finally, we obtain the optimal hull surface of the Wigley model through genetic algorithms, reducing the wave resistance almost by a half. It is shown that the wavelet method may provide a new approach for the hull optimization.
基金the National Natural Science Foundation of China(No.51009087)
文摘The hull form optimization concerns one of the most important applications of wave-making resistance theories.In recent years,scholars can determine the hull form by using the optimization method based on the computational fluid dynamics(CFD)and other mathematical techniques.In this paper,the hull form optimization method based on the Rankine source method and nonlinear programming(NLP)is discussed;in the optimization process,a hull form modification function is introduced to represent an improved hull surface and to generate a new smooth hull surface by changing its frame lines and bow stem profiles under the prescribed design constraints. Numerical example is given for a practical container hull form.Finally,shape optimization of bow bulls is shown for non-protruding and protruding bow bulls.This study presents a simplified and practical design method to the select frame lines of bow bulls.
基金supported by the Key Area Research and Development Program of Guangdong Province(Grant Nos.2021B0101200002,2021B0202070002)the Natural Science Foundation of Guangdong Province(Grant Nos.2022A1515011285,2021A1515011771)Innovation Group Project of Southern Marine Science and Engineering Guangdong Laboratory(Zhuhai)(Grant No.SML2022008).
文摘The hydrodynamic shape of the heaving buoy is an important factor of the motion response in waves and thus concerns the energy conversion efficiency for the point absorbers(PAs).The current experience-based designs are time consuming and not very efficient,hence,faster and smarter methods are desirable.An automated optimization method based on a fully parametric modeling method and computational fluid dynamics(CFD),is proposed in this paper.Using this method,a benchmark buoy is screen designed and then optimized by maximizing the heave motion response.The geometry is described parametrically and deformed by means of the free-form deformation(FFD)method.During the optimization process,the expansion factor of control points is the basis for the variations.A combination of the Sobol and the non-dominated sorting genetic algorithm II(NSGA-II)is used to search for the solutions.After several iterations,the heaving buoy shape with optimal heave motion response is obtained.The analyses show that the heave motion response has increased 55.3%after optimization.The developed methodology is valid and seems to be a promising way to design a novel buoy that can significantly improve the wave energy conversion efficiency of the PAs in future.
文摘Weak Galerkin finite element method is introduced for solving wave equation with interface on weak Galerkin finite element space(Pk(K),P_(k−1)(∂K),[P_(k−1)(K)]^(2)).Optimal order a priori error estimates for both space-discrete scheme and implicit fully discrete scheme are derived in L1(L2)norm.This method uses totally discontinuous functions in approximation space and allows the usage of finite element partitions consisting of general polygonal meshes.Finite element algorithm presented here can contribute to a variety of hyperbolic problems where physical domain consists of heterogeneous media.
基金supported by the National High Technology Research and Development Program of China (No.2009AA03Z413)
文摘In order to reduce transmission loss of the optical waveguide in Mach-Zehnder (M-Z) electro-optical (EO) polymer modulator,the basic iterative formula of semi-vector finite-difference beam propagation method (FD-BPM) is obtained from the scalar wave equation. The transition waveguide is combined with S-type bend branch waveguide for the M-Z EO modulator in the branch waveguide. The effects of structure parameters such as ridge width,length of the branch waveguide and interferometer spacing on the transmission loss are systematically studied by using the semi-vector FD-BPM method. The structure is optimized as an S-sine bend branch waveguide,with rib width w=7μm,length of branch waveguide L=1200μm and interferometer spacing G=22 μm. The results show that the optimized structure can reduce transmission loss to 0.083 dB,which have a certain reference value to the design of optical waveguide in M-Z polymer modulator.