Based on Least Square Method, this paper presents variational principle for handling various water gravity wave theories and the unified water gravity wave theory was given. By using this variational principle of unif...Based on Least Square Method, this paper presents variational principle for handling various water gravity wave theories and the unified water gravity wave theory was given. By using this variational principle of unified water wave theory, two kinds of improved linear waves were derived. The first one uses the same boundary conditions which were applied to derive 5-order Stokes wave. The second one uses the accurate boundary conditions (Eqs. 11 and 12). The two improved linear waves were compared with the existing linear wave.展开更多
A new theory is developed here for evaluating solitary waves on water, with results of high accuracy uniformly valid for waves of all heights, from the highest wave with a corner crest of 120<SUP></SUP> do...A new theory is developed here for evaluating solitary waves on water, with results of high accuracy uniformly valid for waves of all heights, from the highest wave with a corner crest of 120<SUP></SUP> down to very low ones of diminishing height. Solutions are sought for the Euler model by employing a unified expansion of the logarithmic hodograph in terms of a set of intrinsic component functions analytically determined to represent all the intrinsic properties of the wave entity from the wave crest to its outskirts. The unknown coefficients in the expansion are determined by minimization of the mean-square error of the solution, with the minimization optimized so as to take as few terms as needed to attain results as high in accuracy as attainable. In this regard, Stokess formula, F<SUP>2</SUP>= tan , relating the wave speed (the Froude number F) and the logarithmic decrement of its wave field in the outskirt, is generalized to establish a new criterion requiring (for minimizing solution error) the functional expansion to contain a finite power series in M terms of Stokess basic term (singular in ), such that 2M is just somewhat beyond unity, i.e. 2M1. This fundamental criterion is fully validated by solutions for waves of various amplitude-to-water depth ratio =a/h, especially about 0.01, at which M=10 by the criterion. In this pursuit, the class of dwarf solitary waves, defined for waves with 0.01, is discovered as a group of problems more challenging than even the highest wave. For the highest wave, a new solution is determined here to give the maximum height <SUB>hst</SUB>=0.8331990, and speed F<SUB>hst</SUB>=1.290890, accurate to the last significant figure, which seems to be a new record.展开更多
Based on theoretical system of current Maxwell’s equations, the Maxwell’s equations for LEM waves concealed in full current law and Faraday’s law of electromagnetic induction (Faraday’s law) are proposed. Then, ta...Based on theoretical system of current Maxwell’s equations, the Maxwell’s equations for LEM waves concealed in full current law and Faraday’s law of electromagnetic induction (Faraday’s law) are proposed. Then, taking them as the fundamental equations, the wave equation and energy equation of LEM waves are established, and a new electromagnetic wave propagation mode based on the mutual induction of scalar electromagnetic fields/vortex magneto-electric fields, which was overlooked in current Maxwell’s equations, are put forward. Moreover, through theoretical derivation based on vacuum LEM waves, the Maxwell’s equations of the gravitational field generated by vacuum LEM waves, the wave equations of the electromagnetic scalar potential/magnetic vector potential and the constraint equation governing the wave phase-velocities between LEM/TEM waves are discovered. Finally, on the basis of these theoretical research results, the electromagnetic properties of vacuum LEM waves are analyzed in detail, encompassing the speed of light, harmless penetrability to the human body, absorption and stable storage by water, the possibility of generating artificial gravitational fields, and the capability of extracting free energy. This reveals the medical functional mechanism of LEM waves and establishes a solid theoretical basis for the application of LEM waves in the fields of medicine and energy.展开更多
Very Large Floating Structures (VLFS) have drawn considerable attention recently due to their potential significance in the exploitation of ocean resources and in the utilization of ocean space. Efficient and accurate...Very Large Floating Structures (VLFS) have drawn considerable attention recently due to their potential significance in the exploitation of ocean resources and in the utilization of ocean space. Efficient and accurate estimation of their hydroelastic responses to waves is very important for the design. Recently, an efficient numerical algorithm was developed by Ertekin and Kim (1999). However, in their analysis, the linear Level I Green-Naghdi (GN) theory is employed to describe fluid dynamics instead of the conventional linear wave (LW) theory of finite water depth. They claimed that this linear level I GN theory provided better predictions of the hydroelastic responses of VLFS than the linear wave theory. In this paper, a detailed derivation is given in the conventional linear wave theory framework with the same quantity as used in the linear level I GN theory framework. This allows a critical comparison between the linear wave theory and the linear level I GN theory. It is found that the linear level I GN theory can be regarded as an approximation to the linear wave theory of finite water depth. The consequences of the differences between these two theories in the predicted hydroelastic responses are studied quantitatively. And it is found that the linear level I GN theory is not superior to the linear wave theory. Finally, various factors affecting the hydroelastic response of VLFS are studied with the implemented algorithm.展开更多
A unified perturbation theory is developed here for calculating solitary waves of all heights by series expansion of base flow variables in powers of a small base parameter to eighteenth order for the one-parameter fa...A unified perturbation theory is developed here for calculating solitary waves of all heights by series expansion of base flow variables in powers of a small base parameter to eighteenth order for the one-parameter family of solutions in exact form, with all the coefficients determined in rational numbers. Comparative studies are pursued to investigate the effects due to changes of base parameters on （i） the accuracy of the theoretically predicted wave properties and （ii） the rate of convergence of perturbation expansion. Two important results are found by comparisons between the theoretical predictions based on a set of parameters separately adopted for expansion in turn. First, the accuracy and the convergence of the perturbation expansions, appraised versus the exact solution provided by an earlier paper [1] as the standard reference, are found to depend, quite sensitively, on changes in base parameter. The resulting variations in the solution are physically displayed in various wave properties with differences found dependent on which property （e.g. the wave amplitude, speed, its profile, excess mass, momentum, and energy）, on what range in value of the base, and on the rank of the order n in the expansion being addressed. Secondly, regarding convergence, the present perturbation series is found definitely asymptotic in nature, with the relative error δ （n） （the relative mean-square difference between successive orders n of wave elevations） reaching a minimum, δm at a specific order, n = n both depending on the base adopted, e.g. nm,α= 11-12 based on parameter α （wave amplitude）, nm,δ = 15 on δ （amplitude-speed square ratio）, and nm.ε= 17 on ε （ wave number squared）. The asymptotic range is brought to completion by the highest order of n = 18 reached in this work.展开更多
Based on the linear wave, solitary wave and fifth order stokes wave derived by use of the Unified Variational Principle of Water Gravity Wave (UVPWGW), this paper derives stream function wave theory by using UVPWGW. T...Based on the linear wave, solitary wave and fifth order stokes wave derived by use of the Unified Variational Principle of Water Gravity Wave (UVPWGW), this paper derives stream function wave theory by using UVPWGW. This paper will handle the Kinematic Free Surface Boundary Condition (KFSBC) and Dynamic Free Surface Boundary Condition (DFSBC) directly and give the optimum solution, instead of the conditions Sigma(Q(av) - Q(i))(2) = min, and the related equations of stational condition. When the wave height H, period T and water depth D are given, the original stream function wave will be determined, and can not be adjusted if it does not agree with the real case; in the present method, the adjustment can be done by adding several constraint conditions, for example, the wave profile can be adjusted according to the condition of accurate peak position. The examples given in this paper show that for the original stream function wave, the DFSBC can be fairly well satisfied, but the KFSBC can not; however, the stream function wave derived by UVPWGW is better than the original one in the sense of minimum error squares in the aspect of the level at which KFSBC and DFSBC are satisfied.展开更多
This work considers initiation of nonlinear waves, their propagation, reflection, and their interactions in thermoelastic solids and thermoviscoelastic solids with and without memory. The conservation and balance laws...This work considers initiation of nonlinear waves, their propagation, reflection, and their interactions in thermoelastic solids and thermoviscoelastic solids with and without memory. The conservation and balance laws constituting the mathematical models as well as the constitutive theories are derived for finite deformation and finite strain using second Piola-Kirchoff stress tensor and Green’s strain tensor and their material derivatives [1]. Fourier heat conduction law with constant conductivity is used as the constitutive theory for heat vector. Numerical studies are performed using space-time variationally consistent finite element formulations derived using space-time residual functionals and the non-linear equations resulting from the first variation of the residual functional are solved using Newton’s Linear Method with line search. Space-time local approximations are considered in higher order scalar product spaces that permit desired order of global differentiability in space and time. Computed results for non-linear wave propagation, reflection, and interaction are compared with linear wave propagation to demonstrate significant differences between the two, the importance of the nonlinear wave propagation over linear wave propagation as well as to illustrate the meritorious features of the mathematical models and the space-time variationally consistent space-time finite element process with time marching in obtaining the numerical solutions of the evolutions.展开更多
The problem of oblique wave （internal wave） propagation over a small deformation in a channel flow consisting of two layers was considered. The upper fluid was assumed to be bounded above by a rigid lid, which is an...The problem of oblique wave （internal wave） propagation over a small deformation in a channel flow consisting of two layers was considered. The upper fluid was assumed to be bounded above by a rigid lid, which is an approximation for the free surface, and the lower one was bounded below by an impermeable bottom surface having a small deformation; the channel was unbounded in the horizontal directions. Assuming irrotational motion, the perturbation technique was employed to calculate the first-order corrections of the velocity potential in the two fluids by using Green＇s integral theorem suitably with the introduction of appropriate Green＇s functions. Those functions help in calculating the reflection and transmission coefficients in terms of integrals involving the shape ftmction c（x） representing the bottom deformation. Three-dimensional linear water wave theory was utilized for formulating the relevant boundary value problem. Two special examples of bottom deformation were considered to validate the results. Consideration of a patch of sinusoidal ripples （having the same wave number） shows that the reflection coefficient is an oscillatory function of the ratio of twice the x-component of the wave number to the ripple wave number. When this ratio approaches one, the theory predicts a resonant interaction between the bed and the interface, and the reflection coefficient becomes a multiple of the number of ripples. High reflection of incident wave energy occurs if this number is large. Similar results were observed for a patch of sinusoidal ripples having different wave numbers. It was also observed that for small angles of incidence, the reflected energy is greater compared to other angles of incidence up to π/ 4. These theoretical observations are supported by graphical results.展开更多
The distributions of the wave-induced radiation stress tensor over depth are studied by us- ing the linear wave theory, which are divided into three regions, i. e., above the mean water level, be- low the wave trough ...The distributions of the wave-induced radiation stress tensor over depth are studied by us- ing the linear wave theory, which are divided into three regions, i. e., above the mean water level, be- low the wave trough level, and between these two levels. The computational expressions of the wave-in- duced radiation stress tensor at the arbitrary wave angle are established by means of the Eulerian coordi- nate transformation, and the asymptotic forms for deep and shallow water are also presented. The verti- cal variations of a 30°incident wave-induced radiation stress tensor in deep water, intermediate water and shallow water are calculated respectively. The following conclusions are obtained from computations. The wave-induced radiation stress tensor below the wave trough level is induced by the water wave parti- cle velocities only, whereas both the water wave particle velocities and the wave pressure contribute to the tensor above the wave trough level. The vertical variations of the wave-induced radiation stress ten- sor are influenced substantially by the velocity component in the direction of wave propagation. The dis- tributions of the wave-induced radiation stress tensor over depth are nonuiniform and the proportion of the tensor below the wave trough level becomes considerable in the shallow water. From the water surface to the seabed, the reversed variations occur for the predominant tensor components.展开更多
Based on the previous study results, two higher accuracy explicit solutions to the dispersion equation for wave length are presented in this paper. These two solutions have an accuracy of 0. 1% over all wave lengths, ...Based on the previous study results, two higher accuracy explicit solutions to the dispersion equation for wave length are presented in this paper. These two solutions have an accuracy of 0. 1% over all wave lengths, which is sufficiently complete for practical application. At the same time, several previous explicit solutions also have been reviewed and compared herein. In comparison with accuracy, the results show that the present two solutions are as good as Wu and Thornton＇s solution （which has a good accuracy over all wave lengths, but its calculation formula is so complex that it is hard to be used with a hand calculator）, and are better than the other solutions, they may be rather useful in practical calculation with a hand calculator or computer.展开更多
Obliquely incident water wave scattering by an uneven channel-bed in the form of a small bottom undulation in a two-layer fluid is investigated within the frame work of three-dimensional linear water wave theory. The ...Obliquely incident water wave scattering by an uneven channel-bed in the form of a small bottom undulation in a two-layer fluid is investigated within the frame work of three-dimensional linear water wave theory. The upper fluid is assumed to be bounded above by a rigid lid, while the lower one is bounded below by a bottom surface having a small deformation and the channel is unbounded in the horizontal directions. Assuming irrotational motion, perturbation technique is employed to calculate the first-order corrections to the velocity potentials in the two fluids by using Fourier transform approximately, and also to calculate the reflection and transmission coefficients in terms of integrals involving the shape function representing the bottom deformation. Consideration of a patch of sinusoidal ripples shows that the reflection coefficient is an oscillatory function of the ratio of twice the component of the wave number along x-axis and the ripple wave number. When this ratio approaches one, the theory predicts a resonant interaction between the bed and interface, and the reflection coefficient becomes a multiple of the number of ripples. High reflection of incident wave energy occurs if this number is large.展开更多
Using linear water wave theory,three-dimensional problems concerning the interaction of waves with spherical structures in a fluid which contains a three-layer fluid consisting of a layer of finite depth bounded above...Using linear water wave theory,three-dimensional problems concerning the interaction of waves with spherical structures in a fluid which contains a three-layer fluid consisting of a layer of finite depth bounded above by freshwater of finite depth with free surface and below by an infinite layer of water of greater density are considered.In such a situation timeharmonic waves with a given frequency can propagate with three wavenumbers.The sphere is submerged in either of the three layers.Each problem is reduced to an infinite system of linear equations by employing the method of multipoles and the system of equations is solved numerically by standard technique.The hydrodynamic forces(vertical and horizontal forces)are obtained and depicted graphically against the wavenumber.When the density ratio of the upper and middle layer is made to approximately one,curves for vertical and horizontal forces almost coincide with the corresponding curves for the case of a two-layer fluid with a free surface.This means that in the limit,the density ratio of the upper and middle layer goes to approximately one,the solution agrees with the solution for the case of a two-layer fluid with a free surface.展开更多
A joint probability density function (PDF) for surface slopes in two arbitrary directions is derived on the basis of Longuet Higgins's linear model for three-dimensionol (3-D) random wave field. and the correlatio...A joint probability density function (PDF) for surface slopes in two arbitrary directions is derived on the basis of Longuet Higgins's linear model for three-dimensionol (3-D) random wave field. and the correlation moments of surface slopes. as parameters in the PDF, are expressed in terms of directional spectrum of ocean waves. So long as the directional spectrum model is given, these parameters are determined. Since the directional spectrum models proposed so far are mostly parameterized by the wind speed and fetch, this allows for substituting these parameters with thc wind speed and fetch. As an example, the wind speed and fetch are taken to be 14 m ' s and 200 km, and the Hasselmann and Donclan directional spectra are, respectively, use to compute these parameters. Some novel results a reobtained. One of the increasing interesting results is that the variances of surface slope in downwind and cross-wind directions determined by the Donclan directional spectra are close to those measured by Cox and Munk (1954). Some discussions are made on these results.展开更多
The scattering problem involving water waves by small undulation on the porous ocean-bed in a two-layer fluid,is investigated within the framework of the two-dimensional linear water wave theory where the upper layer ...The scattering problem involving water waves by small undulation on the porous ocean-bed in a two-layer fluid,is investigated within the framework of the two-dimensional linear water wave theory where the upper layer is covered by a thin uniform sheet of ice modeled as a thin elastic plate.In such a two-layer fluid there exist waves with two different modes,one with a lower wave number propagate along the ice-cover whilst those with a higher wave number propagate along the interface.An incident wave of a particular wave number gets reflected and transmitted over the bottom undulation into waves of both modes.Perturbation analysis in conjunction with the Fourier transform technique is used to derive the first-order corrections of reflection and transmission coefficients for both the modes due to incident waves of two different modes.One special type of bottom topography is considered as an example to evaluate the related coefficients in detail.These coefficients are depicted in graphical forms to demonstrate the transformation of wave energy between the two modes and also to illustrate the effects of the ice sheet and the porosity of the undulating bed.展开更多
When designing a wave power plant, reliable and fast simulation tools are required. Computational fluid dynamics （CFD） software provides high accuracy but with a very high computational cost, and in operational, mod...When designing a wave power plant, reliable and fast simulation tools are required. Computational fluid dynamics （CFD） software provides high accuracy but with a very high computational cost, and in operational, moderate sea states, linear potential flow theories may be sufficient to model the hydrodynamics. In this paper, a model is built in COMSOL Multiphysics to solve for the hydrodynamic parameters of a point-absorbing wave energy device. The results are compared with a linear model where the hydrodynamical parameters are computed using WAMIT, and to experimental results from the Lysekil research site. The agreement with experimental data is good for both numerical models.展开更多
Using the linear wave theory, the distributions of the wave induced excess momentum fluxes over depth at the arbitrary wave angle and their asymptotic forms for deep and shallow water are developed. Results indicate ...Using the linear wave theory, the distributions of the wave induced excess momentum fluxes over depth at the arbitrary wave angle and their asymptotic forms for deep and shallow water are developed. Results indicate that the distribution of the wave induced excess momentum fluxes over depth is non uniform and the contributions of the component below the wave trough to the total momentum fluxes become considerable in shallow water. On the basis of the Navier Stokes equations, the simplified three dimensional mathematical model is established by taking a phase average over a wavelength. It is found that there are the terms of the wave induced excess momentum fluxes varying over depth in the model, which illustrates the situation of wave current interactions and the vertical structure of current velocity are changed because of different wave induced excess momentum fluxes at various vertical location. The finite difference method is employed to solve the simplified model. Performances of the two dimensional vertically integrated equations are evaluated against available numerical and experimental results including the cases of wave set up on a plane beach, longshore current due to an oblique wave, wave induced nearshore circulation in a semi enclosed seas, and wave current interactions. All cases yield satisfactory agreements. The three dimensional mathematical model is applied to the numerical simulation of wave current interactions, and it performs well in predicting the vertical velocity structure and the plane flow field.展开更多
The origin of elementary particle mass is considered as a function of n-valued graviton quanta. To develop this concept we begin in a cold region of “empty space” comprised of only microscopic gravitons oscillating ...The origin of elementary particle mass is considered as a function of n-valued graviton quanta. To develop this concept we begin in a cold region of “empty space” comprised of only microscopic gravitons oscillating at angular frequency ω. From opposite directions enters a pair of stray protons. Upon colliding, heat and energy are released. Customarily, this phase and what follows afterward would be described by Quantum Chromodynamics (QCD). Instead, we argue for an intermediary step. One in which neighboring gravitons absorb discrete amounts of plane-wave energy. Captured by the graviton, the planewave becomes a standing wave, whereupon its electromagnetic energy densities are converted into gravitational quanta. Immediately thereafter an elementary particle is formed and emitted, having both mass and spin. From absorption to conversion to emission occurs in less than 3.7 × 10−16 s. During this basic unit of hybrid time, general relativity and quantum physics unite into a common set of physical laws. As additional stray protons collide the process continues. Over eons, vast regions of spacetime become populated with low-mass particles. These we recognize to be dark matter by its effects on large scale structures in the universe. Its counterpart, dark energy, arises when the conversion of gravitational quanta to particle emission is interrupted. This causes the gravitational quanta to be ejected. It is recognized by its large scale effects on the universe.展开更多
文摘Based on Least Square Method, this paper presents variational principle for handling various water gravity wave theories and the unified water gravity wave theory was given. By using this variational principle of unified water wave theory, two kinds of improved linear waves were derived. The first one uses the same boundary conditions which were applied to derive 5-order Stokes wave. The second one uses the accurate boundary conditions (Eqs. 11 and 12). The two improved linear waves were compared with the existing linear wave.
文摘A new theory is developed here for evaluating solitary waves on water, with results of high accuracy uniformly valid for waves of all heights, from the highest wave with a corner crest of 120<SUP></SUP> down to very low ones of diminishing height. Solutions are sought for the Euler model by employing a unified expansion of the logarithmic hodograph in terms of a set of intrinsic component functions analytically determined to represent all the intrinsic properties of the wave entity from the wave crest to its outskirts. The unknown coefficients in the expansion are determined by minimization of the mean-square error of the solution, with the minimization optimized so as to take as few terms as needed to attain results as high in accuracy as attainable. In this regard, Stokess formula, F<SUP>2</SUP>= tan , relating the wave speed (the Froude number F) and the logarithmic decrement of its wave field in the outskirt, is generalized to establish a new criterion requiring (for minimizing solution error) the functional expansion to contain a finite power series in M terms of Stokess basic term (singular in ), such that 2M is just somewhat beyond unity, i.e. 2M1. This fundamental criterion is fully validated by solutions for waves of various amplitude-to-water depth ratio =a/h, especially about 0.01, at which M=10 by the criterion. In this pursuit, the class of dwarf solitary waves, defined for waves with 0.01, is discovered as a group of problems more challenging than even the highest wave. For the highest wave, a new solution is determined here to give the maximum height <SUB>hst</SUB>=0.8331990, and speed F<SUB>hst</SUB>=1.290890, accurate to the last significant figure, which seems to be a new record.
文摘Based on theoretical system of current Maxwell’s equations, the Maxwell’s equations for LEM waves concealed in full current law and Faraday’s law of electromagnetic induction (Faraday’s law) are proposed. Then, taking them as the fundamental equations, the wave equation and energy equation of LEM waves are established, and a new electromagnetic wave propagation mode based on the mutual induction of scalar electromagnetic fields/vortex magneto-electric fields, which was overlooked in current Maxwell’s equations, are put forward. Moreover, through theoretical derivation based on vacuum LEM waves, the Maxwell’s equations of the gravitational field generated by vacuum LEM waves, the wave equations of the electromagnetic scalar potential/magnetic vector potential and the constraint equation governing the wave phase-velocities between LEM/TEM waves are discovered. Finally, on the basis of these theoretical research results, the electromagnetic properties of vacuum LEM waves are analyzed in detail, encompassing the speed of light, harmless penetrability to the human body, absorption and stable storage by water, the possibility of generating artificial gravitational fields, and the capability of extracting free energy. This reveals the medical functional mechanism of LEM waves and establishes a solid theoretical basis for the application of LEM waves in the fields of medicine and energy.
基金by the National Natural Science Foundation of China(50039010)the Science and Technology Development Foundation of Shanghai Municipal Government(00XD14015)
文摘Very Large Floating Structures (VLFS) have drawn considerable attention recently due to their potential significance in the exploitation of ocean resources and in the utilization of ocean space. Efficient and accurate estimation of their hydroelastic responses to waves is very important for the design. Recently, an efficient numerical algorithm was developed by Ertekin and Kim (1999). However, in their analysis, the linear Level I Green-Naghdi (GN) theory is employed to describe fluid dynamics instead of the conventional linear wave (LW) theory of finite water depth. They claimed that this linear level I GN theory provided better predictions of the hydroelastic responses of VLFS than the linear wave theory. In this paper, a detailed derivation is given in the conventional linear wave theory framework with the same quantity as used in the linear level I GN theory framework. This allows a critical comparison between the linear wave theory and the linear level I GN theory. It is found that the linear level I GN theory can be regarded as an approximation to the linear wave theory of finite water depth. The consequences of the differences between these two theories in the predicted hydroelastic responses are studied quantitatively. And it is found that the linear level I GN theory is not superior to the linear wave theory. Finally, various factors affecting the hydroelastic response of VLFS are studied with the implemented algorithm.
基金The project partly supported by the National Natural Science Foundation of China(19925414,10474045)
文摘A unified perturbation theory is developed here for calculating solitary waves of all heights by series expansion of base flow variables in powers of a small base parameter to eighteenth order for the one-parameter family of solutions in exact form, with all the coefficients determined in rational numbers. Comparative studies are pursued to investigate the effects due to changes of base parameters on （i） the accuracy of the theoretically predicted wave properties and （ii） the rate of convergence of perturbation expansion. Two important results are found by comparisons between the theoretical predictions based on a set of parameters separately adopted for expansion in turn. First, the accuracy and the convergence of the perturbation expansions, appraised versus the exact solution provided by an earlier paper [1] as the standard reference, are found to depend, quite sensitively, on changes in base parameter. The resulting variations in the solution are physically displayed in various wave properties with differences found dependent on which property （e.g. the wave amplitude, speed, its profile, excess mass, momentum, and energy）, on what range in value of the base, and on the rank of the order n in the expansion being addressed. Secondly, regarding convergence, the present perturbation series is found definitely asymptotic in nature, with the relative error δ （n） （the relative mean-square difference between successive orders n of wave elevations） reaching a minimum, δm at a specific order, n = n both depending on the base adopted, e.g. nm,α= 11-12 based on parameter α （wave amplitude）, nm,δ = 15 on δ （amplitude-speed square ratio）, and nm.ε= 17 on ε （ wave number squared）. The asymptotic range is brought to completion by the highest order of n = 18 reached in this work.
文摘Based on the linear wave, solitary wave and fifth order stokes wave derived by use of the Unified Variational Principle of Water Gravity Wave (UVPWGW), this paper derives stream function wave theory by using UVPWGW. This paper will handle the Kinematic Free Surface Boundary Condition (KFSBC) and Dynamic Free Surface Boundary Condition (DFSBC) directly and give the optimum solution, instead of the conditions Sigma(Q(av) - Q(i))(2) = min, and the related equations of stational condition. When the wave height H, period T and water depth D are given, the original stream function wave will be determined, and can not be adjusted if it does not agree with the real case; in the present method, the adjustment can be done by adding several constraint conditions, for example, the wave profile can be adjusted according to the condition of accurate peak position. The examples given in this paper show that for the original stream function wave, the DFSBC can be fairly well satisfied, but the KFSBC can not; however, the stream function wave derived by UVPWGW is better than the original one in the sense of minimum error squares in the aspect of the level at which KFSBC and DFSBC are satisfied.
文摘This work considers initiation of nonlinear waves, their propagation, reflection, and their interactions in thermoelastic solids and thermoviscoelastic solids with and without memory. The conservation and balance laws constituting the mathematical models as well as the constitutive theories are derived for finite deformation and finite strain using second Piola-Kirchoff stress tensor and Green’s strain tensor and their material derivatives [1]. Fourier heat conduction law with constant conductivity is used as the constitutive theory for heat vector. Numerical studies are performed using space-time variationally consistent finite element formulations derived using space-time residual functionals and the non-linear equations resulting from the first variation of the residual functional are solved using Newton’s Linear Method with line search. Space-time local approximations are considered in higher order scalar product spaces that permit desired order of global differentiability in space and time. Computed results for non-linear wave propagation, reflection, and interaction are compared with linear wave propagation to demonstrate significant differences between the two, the importance of the nonlinear wave propagation over linear wave propagation as well as to illustrate the meritorious features of the mathematical models and the space-time variationally consistent space-time finite element process with time marching in obtaining the numerical solutions of the evolutions.
文摘The problem of oblique wave （internal wave） propagation over a small deformation in a channel flow consisting of two layers was considered. The upper fluid was assumed to be bounded above by a rigid lid, which is an approximation for the free surface, and the lower one was bounded below by an impermeable bottom surface having a small deformation; the channel was unbounded in the horizontal directions. Assuming irrotational motion, the perturbation technique was employed to calculate the first-order corrections of the velocity potential in the two fluids by using Green＇s integral theorem suitably with the introduction of appropriate Green＇s functions. Those functions help in calculating the reflection and transmission coefficients in terms of integrals involving the shape ftmction c（x） representing the bottom deformation. Three-dimensional linear water wave theory was utilized for formulating the relevant boundary value problem. Two special examples of bottom deformation were considered to validate the results. Consideration of a patch of sinusoidal ripples （having the same wave number） shows that the reflection coefficient is an oscillatory function of the ratio of twice the x-component of the wave number to the ripple wave number. When this ratio approaches one, the theory predicts a resonant interaction between the bed and the interface, and the reflection coefficient becomes a multiple of the number of ripples. High reflection of incident wave energy occurs if this number is large. Similar results were observed for a patch of sinusoidal ripples having different wave numbers. It was also observed that for small angles of incidence, the reflected energy is greater compared to other angles of incidence up to π/ 4. These theoretical observations are supported by graphical results.
基金The project was supported by the Research Fund for the Doctoral Program of Higher Education of China under contractNo. 9802940
文摘The distributions of the wave-induced radiation stress tensor over depth are studied by us- ing the linear wave theory, which are divided into three regions, i. e., above the mean water level, be- low the wave trough level, and between these two levels. The computational expressions of the wave-in- duced radiation stress tensor at the arbitrary wave angle are established by means of the Eulerian coordi- nate transformation, and the asymptotic forms for deep and shallow water are also presented. The verti- cal variations of a 30°incident wave-induced radiation stress tensor in deep water, intermediate water and shallow water are calculated respectively. The following conclusions are obtained from computations. The wave-induced radiation stress tensor below the wave trough level is induced by the water wave parti- cle velocities only, whereas both the water wave particle velocities and the wave pressure contribute to the tensor above the wave trough level. The vertical variations of the wave-induced radiation stress ten- sor are influenced substantially by the velocity component in the direction of wave propagation. The dis- tributions of the wave-induced radiation stress tensor over depth are nonuiniform and the proportion of the tensor below the wave trough level becomes considerable in the shallow water. From the water surface to the seabed, the reversed variations occur for the predominant tensor components.
基金This study was financially supported by the Doctor Degree ProgramFoundation of the Ministry of Education of China(Grant No.20050294009)
文摘Based on the previous study results, two higher accuracy explicit solutions to the dispersion equation for wave length are presented in this paper. These two solutions have an accuracy of 0. 1% over all wave lengths, which is sufficiently complete for practical application. At the same time, several previous explicit solutions also have been reviewed and compared herein. In comparison with accuracy, the results show that the present two solutions are as good as Wu and Thornton＇s solution （which has a good accuracy over all wave lengths, but its calculation formula is so complex that it is hard to be used with a hand calculator）, and are better than the other solutions, they may be rather useful in practical calculation with a hand calculator or computer.
基金partially supported by a research grant from Department of Science and Technology(DST),India(No.SB/FTP/MS-003/2013)
文摘Obliquely incident water wave scattering by an uneven channel-bed in the form of a small bottom undulation in a two-layer fluid is investigated within the frame work of three-dimensional linear water wave theory. The upper fluid is assumed to be bounded above by a rigid lid, while the lower one is bounded below by a bottom surface having a small deformation and the channel is unbounded in the horizontal directions. Assuming irrotational motion, perturbation technique is employed to calculate the first-order corrections to the velocity potentials in the two fluids by using Fourier transform approximately, and also to calculate the reflection and transmission coefficients in terms of integrals involving the shape function representing the bottom deformation. Consideration of a patch of sinusoidal ripples shows that the reflection coefficient is an oscillatory function of the ratio of twice the component of the wave number along x-axis and the ripple wave number. When this ratio approaches one, the theory predicts a resonant interaction between the bed and interface, and the reflection coefficient becomes a multiple of the number of ripples. High reflection of incident wave energy occurs if this number is large.
文摘Using linear water wave theory,three-dimensional problems concerning the interaction of waves with spherical structures in a fluid which contains a three-layer fluid consisting of a layer of finite depth bounded above by freshwater of finite depth with free surface and below by an infinite layer of water of greater density are considered.In such a situation timeharmonic waves with a given frequency can propagate with three wavenumbers.The sphere is submerged in either of the three layers.Each problem is reduced to an infinite system of linear equations by employing the method of multipoles and the system of equations is solved numerically by standard technique.The hydrodynamic forces(vertical and horizontal forces)are obtained and depicted graphically against the wavenumber.When the density ratio of the upper and middle layer is made to approximately one,curves for vertical and horizontal forces almost coincide with the corresponding curves for the case of a two-layer fluid with a free surface.This means that in the limit,the density ratio of the upper and middle layer goes to approximately one,the solution agrees with the solution for the case of a two-layer fluid with a free surface.
基金This work is financially supported by the National Natural Science Foundation of China(No.49676277)863-818 Project(05-02)
文摘A joint probability density function (PDF) for surface slopes in two arbitrary directions is derived on the basis of Longuet Higgins's linear model for three-dimensionol (3-D) random wave field. and the correlation moments of surface slopes. as parameters in the PDF, are expressed in terms of directional spectrum of ocean waves. So long as the directional spectrum model is given, these parameters are determined. Since the directional spectrum models proposed so far are mostly parameterized by the wind speed and fetch, this allows for substituting these parameters with thc wind speed and fetch. As an example, the wind speed and fetch are taken to be 14 m ' s and 200 km, and the Hasselmann and Donclan directional spectra are, respectively, use to compute these parameters. Some novel results a reobtained. One of the increasing interesting results is that the variances of surface slope in downwind and cross-wind directions determined by the Donclan directional spectra are close to those measured by Cox and Munk (1954). Some discussions are made on these results.
基金Supprted by the ISIRD grant(Ref.No.16-3/10/IITRPR/Acad/116)
文摘The scattering problem involving water waves by small undulation on the porous ocean-bed in a two-layer fluid,is investigated within the framework of the two-dimensional linear water wave theory where the upper layer is covered by a thin uniform sheet of ice modeled as a thin elastic plate.In such a two-layer fluid there exist waves with two different modes,one with a lower wave number propagate along the ice-cover whilst those with a higher wave number propagate along the interface.An incident wave of a particular wave number gets reflected and transmitted over the bottom undulation into waves of both modes.Perturbation analysis in conjunction with the Fourier transform technique is used to derive the first-order corrections of reflection and transmission coefficients for both the modes due to incident waves of two different modes.One special type of bottom topography is considered as an example to evaluate the related coefficients in detail.These coefficients are depicted in graphical forms to demonstrate the transformation of wave energy between the two modes and also to illustrate the effects of the ice sheet and the porosity of the undulating bed.
基金supported by the Center for Natural Disaster Science(CNDS)the Swedish Research Council(VR,Grant Number 2015-04657)+1 种基金Lars Hiertas FoundationBengt Ingestrms scholarship fund
文摘When designing a wave power plant, reliable and fast simulation tools are required. Computational fluid dynamics （CFD） software provides high accuracy but with a very high computational cost, and in operational, moderate sea states, linear potential flow theories may be sufficient to model the hydrodynamics. In this paper, a model is built in COMSOL Multiphysics to solve for the hydrodynamic parameters of a point-absorbing wave energy device. The results are compared with a linear model where the hydrodynamical parameters are computed using WAMIT, and to experimental results from the Lysekil research site. The agreement with experimental data is good for both numerical models.
文摘Using the linear wave theory, the distributions of the wave induced excess momentum fluxes over depth at the arbitrary wave angle and their asymptotic forms for deep and shallow water are developed. Results indicate that the distribution of the wave induced excess momentum fluxes over depth is non uniform and the contributions of the component below the wave trough to the total momentum fluxes become considerable in shallow water. On the basis of the Navier Stokes equations, the simplified three dimensional mathematical model is established by taking a phase average over a wavelength. It is found that there are the terms of the wave induced excess momentum fluxes varying over depth in the model, which illustrates the situation of wave current interactions and the vertical structure of current velocity are changed because of different wave induced excess momentum fluxes at various vertical location. The finite difference method is employed to solve the simplified model. Performances of the two dimensional vertically integrated equations are evaluated against available numerical and experimental results including the cases of wave set up on a plane beach, longshore current due to an oblique wave, wave induced nearshore circulation in a semi enclosed seas, and wave current interactions. All cases yield satisfactory agreements. The three dimensional mathematical model is applied to the numerical simulation of wave current interactions, and it performs well in predicting the vertical velocity structure and the plane flow field.
文摘The origin of elementary particle mass is considered as a function of n-valued graviton quanta. To develop this concept we begin in a cold region of “empty space” comprised of only microscopic gravitons oscillating at angular frequency ω. From opposite directions enters a pair of stray protons. Upon colliding, heat and energy are released. Customarily, this phase and what follows afterward would be described by Quantum Chromodynamics (QCD). Instead, we argue for an intermediary step. One in which neighboring gravitons absorb discrete amounts of plane-wave energy. Captured by the graviton, the planewave becomes a standing wave, whereupon its electromagnetic energy densities are converted into gravitational quanta. Immediately thereafter an elementary particle is formed and emitted, having both mass and spin. From absorption to conversion to emission occurs in less than 3.7 × 10−16 s. During this basic unit of hybrid time, general relativity and quantum physics unite into a common set of physical laws. As additional stray protons collide the process continues. Over eons, vast regions of spacetime become populated with low-mass particles. These we recognize to be dark matter by its effects on large scale structures in the universe. Its counterpart, dark energy, arises when the conversion of gravitational quanta to particle emission is interrupted. This causes the gravitational quanta to be ejected. It is recognized by its large scale effects on the universe.
