This study investigated the interaction between irregular waves and vertical breakwater.The main goal was to determine the wave force on the breakwater for different depths,along with the evolution of the secondary wa...This study investigated the interaction between irregular waves and vertical breakwater.The main goal was to determine the wave force on the breakwater for different depths,along with the evolution of the secondary wave generated by overtopping.The open source code DualSPHysics was used in the simulation.Wavelet transform was employed to remove the acoustic components from the weakly compressible smoothed particle hydrodynamics pressure solution.The results of the experiments and simulations showed that although the wave height at the low water level(case 2 was smaller than that at the high water level(case 1,the horizontal and uplifting force of case 2 was higher than that of case 1.In case 2,a large impact pressure occurred because of the plunging wave and wave breaking effects.As the water level was close to the breakwater top in case 1,the secondary wave generated by overtopping was studied.Results indicated that wave and breakwater interactions transfer wave energy from the leading wave component to higher harmonics.Two main harmonics were identified in the secondary wave.The first harmonic was the incident peak frequency,which was still the dominant component.During the wave propagation behind the breakwater,the frequencies of harmonics remained unchanged,but the amplitude of the first harmonic showed an obvious change.Under the same wave condition,the first and second harmonics did not change with depth.展开更多
To address the problem of subdividing inflexible rectangular grid models and their poor definition of velocity interfaces,we propose a complex structure triangular net for a minimum traveltime ray tracing global algor...To address the problem of subdividing inflexible rectangular grid models and their poor definition of velocity interfaces,we propose a complex structure triangular net for a minimum traveltime ray tracing global algorithm.Our procedure is:(1) Subdivide a triangle grid based on the Delaunay triangular subdivision criterion and the relationships of the points,lines,and the surfaces in the subdividing area.(2) Define the topology relationships and related concepts of triangular unit ray tracing.(3) The source point and wave arrival points at any time compose the propagating plane wave and the minimum traveltime and secondary source positions are calculated during the plane wave propagation.We adopt the hyperbolic approximation global algorithm for secondary source retrieving.(4) By minimum traveltime ray tracing,collect the path from receiver to source points with the neighborhood point's traveltime and the direction of the secondary source.Numerical simulation examples are given to test the algorithm.The results show that the triangular net ray tracing method demonstrates model subdivision flexibility,precise velocity discontinuity interfaces,and accurate computations.展开更多
A DTM map of the study area in the Taiwan Shoal was drawn based on precise and high- density data acquired in a field survey by a multi-beam sounding system (R2Sonic2024). We identified sand waves in the study area ...A DTM map of the study area in the Taiwan Shoal was drawn based on precise and high- density data acquired in a field survey by a multi-beam sounding system (R2Sonic2024). We identified sand waves in the study area at water depths of 13.89-49.12 m; the main sand waves had heights of 5- 25 m, lengths of 0.1-2.0 km, and crest lines 0.1-6.5 km long. The spatial distribution of the sand waves on the seabed is dense in the north and sparse in the south and the directions range between 50°-80° and 90°-135°. Between the main sand waves, secondary sand waves develop with heights of 0.1-5.0 m and lengths of 10-100 m, which are difficult to detect by satellite remote sensing. By comparing the evolution structures of the secondary and main sand waves, we identified three evolution modes of the secondary sand waves: parallel, oblique, and divergent modes according to the relative crest directions. Suitable water depth, reciprocating current speeds between 40 and 100 cm/s, and abundant sediment supply create favorable conditions for the formation of linear sand waves. Comparing the DTM maps and profiles of the June 2012 and June 2013 surveys of the same area, we found that the shape and morphology of the sand waves remained mostly unchanged under normal hydrodynamic conditions.展开更多
An oscillating buoy wave power device (OD) is a device extracting wave power by an oscillating buoy. Being excited by waves, the buoy heaves up and down to convert wave energy into electricity by means of a mechanical...An oscillating buoy wave power device (OD) is a device extracting wave power by an oscillating buoy. Being excited by waves, the buoy heaves up and down to convert wave energy into electricity by means of a mechanical or hydraulic device. Compared with an Oscillating Water Column (OWC) wave power device, the OD has the same capture width ratio as the OWC does, but much higher secondary conversion efficiency. Moreover, the chamber of the OWC, which is the most expensive and difficult part to be built, is not necessary for the OD, so it is easier to construct an OD. In this paper, a numerical calculation is conducted for an optimal design of the OD firstly, then a model of the device is built and, a model test is carried out in a wave tank. The results show that the total efficiency of the OD is much higher than that of the OWC and that the OD is a promising wave power device.展开更多
The equivalent damage calculation formulae of fatigue crack formation andgrowth are established. In order to compile the fatigue crack formation and growth accelerated loadspectra, the main wave shapes and load sequen...The equivalent damage calculation formulae of fatigue crack formation andgrowth are established. In order to compile the fatigue crack formation and growth accelerated loadspectra, the main wave shapes and load sequence of the actual load spectrum are kept constant, andthe carrier waves are cut off. And secondary waves are put together into new secondary waves toshorten the test time according to the equivalent damage calculation formulae respectively. Then bythe fatigue cumulative damage calculation of the fatigue crack formation and growth accelerated loadspectra, the one corresponding to the bigger damage is determined as the fatigue accelerated testload spectrum. Therefore in the test process, the fatigue accelerated test spectrum may be appliedtill fatigue failure, the engineering fatigue crack length of full-scale structure need not beinspected, and the fatigue crack formation accelerated load spectrum need not be transferred intothe fatigue crack growth accelerated load spectrum. Finally, it is verified by tests of two kinds ofspecimens that the damages of the specimens caused by the accelerated load spectra are near tothose by the actual load spectra; namely, the tested life of actual load spectra is similar to thatof accelerated load spectra. But the test time of accelerated load spectra is shortened by aboutthree-quarters that of actual load spectra. From these tests, it is also found that the fatigueaccelerated test spectrum has an advantage over FALSTAFF spectra.展开更多
Parabolized stability equations (PSE) approach is used to investigate problems of secondary instability in supersonic boundary layers. The results show that the mechanism of secondary instability does work, whether ...Parabolized stability equations (PSE) approach is used to investigate problems of secondary instability in supersonic boundary layers. The results show that the mechanism of secondary instability does work, whether the 2-D fundamental disturbance is of the first mode or second mode T-S wave. The variation of the growth rates of the 3-D sub-harmonic wave against its span-wise wave number and the amplitude of the 2-D fundamental wave is found to be similar to those found in incompressible boundary layers. But even as the amplitude of the 2-D wave is as large as the order 2%, the maximum growth rate of the 3-D sub-harmonic is still much smaller than the growth rate of the most unstable second mode 2-D T-S wave. Consequently, secondary instability is unlikely the main cause leading to transition in supersonic boundary layers.展开更多
To investigate the mechanism of secondary circulations in rip current systems, and to explore the relationship between wave conditions and secondary circulation intensity, a series of numerical experiments is performe...To investigate the mechanism of secondary circulations in rip current systems, and to explore the relationship between wave conditions and secondary circulation intensity, a series of numerical experiments is performed using coupled nearshore wave model and circulation model. In these experiments, the rip currents and secondary circulations generated above barred beaches with rip channels are simulated. A comparison experiment is conducted to investigate the formation and hydrodynamics of the secondary circulations. Model results indicate that the secondary circulations consist of alongshore flows driven by wave set-up near the shoreline, part of the feeder currents driven by the wave set-up over the bars, and onshore flows at the end of the rip channel driven by wave breaking and convection. The existence of the secondary circulation barely affects the rip current, but narrows and intensifies the feeder currents. Three groups of experiments of varying incident wave conditions are performed to investigate the relationship between wave conditions and secondary circulation intensity. The velocity of the alongshore flow of the secondary circulation is sensitive to the variation of the incident wave height and water depth. It is also found that the alongshore flow intensity is in direct proportion to the alongshore variation of the wave height gradient between the bars and the shoreline.展开更多
A preliminary investigation is conducted to study the characteristics of sheath damping near a dielectric wall with secondary electron emission (SEE). Making use of the linear analysis of the sheath stability, we ha...A preliminary investigation is conducted to study the characteristics of sheath damping near a dielectric wall with secondary electron emission (SEE). Making use of the linear analysis of the sheath stability, we have found two major contributions to the sheath damping, one similar to the Landau damping in uniform plasmas and another determined by local electric field and electron density of the steady-state sheath. It indicates that in a classical sheath regime the damping in the sheath region monotonically increases towards the wall and decreases with the enhancement of SEE effect. In order to verify the theoretical analysis, sheath oscillation processes induced by an initial disturbance are simulated with a time-dependent one-dimensional (1D) sheath model. Numerical results obtained are consistent with the theoretical analysis qualitatively.展开更多
By direct numerical simulation, a supersonic boundary layer was studied to see whether the mechanism for the generation of sub-harmonic waves, similar to those for the incompressible flows, existed in the process of l...By direct numerical simulation, a supersonic boundary layer was studied to see whether the mechanism for the generation of sub-harmonic waves, similar to those for the incompressible flows, existed in the process of laminar-turbulent transition. The results showed that mechanisms of both resonant triad and secondary instability did exist. Discussions were made on whether these two mechanisms are really important in laminarturbulent transition.展开更多
A fixed oscillating water column(OWC)-type wave energy converter consists of an air chamber,an air turbine and a generator.The energy conversion processes are the primary conversion in an air chamber and the secondary...A fixed oscillating water column(OWC)-type wave energy converter consists of an air chamber,an air turbine and a generator.The energy conversion processes are the primary conversion in an air chamber and the secondary conversion of the turbine.For the practical use,it is necessary to develop a design method which can consider the incident wave motion,the motion of the internal free surface affected in the structure such as a partly submerged wall,the fluctuation of air pressure in an air chamber,and the rotation of the air turbine.At here,the authors carried out the wave tank tests using the model OWC equipped with the impulse turbine and a generator to obtain the experimental data needed to make this design method.As the result,the efficiencies of the three cases with different speed ratio between generator and turbine,and the effects of the curtain wall depth and the wave length on the energy conversion performance were clarified.展开更多
基金the National Natural Sci-ence Foundation of China-Shandong Joint Fund(No.U1706226)the National Natural Science Foundation of China(No.51779236)the Key Science and Tech-nology Projects of International Cooperation of Ministry of Transport,China(Nos.2019-GH-002,TKS190304).
文摘This study investigated the interaction between irregular waves and vertical breakwater.The main goal was to determine the wave force on the breakwater for different depths,along with the evolution of the secondary wave generated by overtopping.The open source code DualSPHysics was used in the simulation.Wavelet transform was employed to remove the acoustic components from the weakly compressible smoothed particle hydrodynamics pressure solution.The results of the experiments and simulations showed that although the wave height at the low water level(case 2 was smaller than that at the high water level(case 1,the horizontal and uplifting force of case 2 was higher than that of case 1.In case 2,a large impact pressure occurred because of the plunging wave and wave breaking effects.As the water level was close to the breakwater top in case 1,the secondary wave generated by overtopping was studied.Results indicated that wave and breakwater interactions transfer wave energy from the leading wave component to higher harmonics.Two main harmonics were identified in the secondary wave.The first harmonic was the incident peak frequency,which was still the dominant component.During the wave propagation behind the breakwater,the frequencies of harmonics remained unchanged,but the amplitude of the first harmonic showed an obvious change.Under the same wave condition,the first and second harmonics did not change with depth.
基金the National Natural Science Foundation of China(Grant Nos.50974081,50774051, 51034003)the Foundation for the Author of National Excellent Doctoral Dissertation of PR China(Grant No.200958)+1 种基金the Program for Changjiang Scholars and Innovative Research Team in University(Grant No.IRT0843)the Open Research Fund Program of the State Key Laboratory of Coal Resources and Safe Mining(China University of Mining and Technology, Beijing)(No.2007-07)
文摘To address the problem of subdividing inflexible rectangular grid models and their poor definition of velocity interfaces,we propose a complex structure triangular net for a minimum traveltime ray tracing global algorithm.Our procedure is:(1) Subdivide a triangle grid based on the Delaunay triangular subdivision criterion and the relationships of the points,lines,and the surfaces in the subdividing area.(2) Define the topology relationships and related concepts of triangular unit ray tracing.(3) The source point and wave arrival points at any time compose the propagating plane wave and the minimum traveltime and secondary source positions are calculated during the plane wave propagation.We adopt the hyperbolic approximation global algorithm for secondary source retrieving.(4) By minimum traveltime ray tracing,collect the path from receiver to source points with the neighborhood point's traveltime and the direction of the secondary source.Numerical simulation examples are given to test the algorithm.The results show that the triangular net ray tracing method demonstrates model subdivision flexibility,precise velocity discontinuity interfaces,and accurate computations.
基金Funding was provided by the Marine Public Special Project (201105001 and 201205004)the Natural Science Foundation (China project, 41276058)
文摘A DTM map of the study area in the Taiwan Shoal was drawn based on precise and high- density data acquired in a field survey by a multi-beam sounding system (R2Sonic2024). We identified sand waves in the study area at water depths of 13.89-49.12 m; the main sand waves had heights of 5- 25 m, lengths of 0.1-2.0 km, and crest lines 0.1-6.5 km long. The spatial distribution of the sand waves on the seabed is dense in the north and sparse in the south and the directions range between 50°-80° and 90°-135°. Between the main sand waves, secondary sand waves develop with heights of 0.1-5.0 m and lengths of 10-100 m, which are difficult to detect by satellite remote sensing. By comparing the evolution structures of the secondary and main sand waves, we identified three evolution modes of the secondary sand waves: parallel, oblique, and divergent modes according to the relative crest directions. Suitable water depth, reciprocating current speeds between 40 and 100 cm/s, and abundant sediment supply create favorable conditions for the formation of linear sand waves. Comparing the DTM maps and profiles of the June 2012 and June 2013 surveys of the same area, we found that the shape and morphology of the sand waves remained mostly unchanged under normal hydrodynamic conditions.
基金This work was financially supported by the National Natural Science Foundation of China(Grant No.59979025),the High Tech Research and Development(863)Program(Grant No.2001AA516010)and the Science&Technology Program of Guangdong(C32004)
文摘An oscillating buoy wave power device (OD) is a device extracting wave power by an oscillating buoy. Being excited by waves, the buoy heaves up and down to convert wave energy into electricity by means of a mechanical or hydraulic device. Compared with an Oscillating Water Column (OWC) wave power device, the OD has the same capture width ratio as the OWC does, but much higher secondary conversion efficiency. Moreover, the chamber of the OWC, which is the most expensive and difficult part to be built, is not necessary for the OD, so it is easier to construct an OD. In this paper, a numerical calculation is conducted for an optimal design of the OD firstly, then a model of the device is built and, a model test is carried out in a wave tank. The results show that the total efficiency of the OD is much higher than that of the OWC and that the OD is a promising wave power device.
基金This project is supported by National Natural Science and Aeronautic Science Foundation of China (No. 50005003, No. 01A51011)
文摘The equivalent damage calculation formulae of fatigue crack formation andgrowth are established. In order to compile the fatigue crack formation and growth accelerated loadspectra, the main wave shapes and load sequence of the actual load spectrum are kept constant, andthe carrier waves are cut off. And secondary waves are put together into new secondary waves toshorten the test time according to the equivalent damage calculation formulae respectively. Then bythe fatigue cumulative damage calculation of the fatigue crack formation and growth accelerated loadspectra, the one corresponding to the bigger damage is determined as the fatigue accelerated testload spectrum. Therefore in the test process, the fatigue accelerated test spectrum may be appliedtill fatigue failure, the engineering fatigue crack length of full-scale structure need not beinspected, and the fatigue crack formation accelerated load spectrum need not be transferred intothe fatigue crack growth accelerated load spectrum. Finally, it is verified by tests of two kinds ofspecimens that the damages of the specimens caused by the accelerated load spectra are near tothose by the actual load spectra; namely, the tested life of actual load spectra is similar to thatof accelerated load spectra. But the test time of accelerated load spectra is shortened by aboutthree-quarters that of actual load spectra. From these tests, it is also found that the fatigueaccelerated test spectrum has an advantage over FALSTAFF spectra.
基金Project supported by the National Natural Science Foundation of China(Nos.10632050,90716007)the Foundation of LIU Hui Center of Applied Mathematics of Nankai University and Tianjin University
文摘Parabolized stability equations (PSE) approach is used to investigate problems of secondary instability in supersonic boundary layers. The results show that the mechanism of secondary instability does work, whether the 2-D fundamental disturbance is of the first mode or second mode T-S wave. The variation of the growth rates of the 3-D sub-harmonic wave against its span-wise wave number and the amplitude of the 2-D fundamental wave is found to be similar to those found in incompressible boundary layers. But even as the amplitude of the 2-D wave is as large as the order 2%, the maximum growth rate of the 3-D sub-harmonic is still much smaller than the growth rate of the most unstable second mode 2-D T-S wave. Consequently, secondary instability is unlikely the main cause leading to transition in supersonic boundary layers.
基金supported by China’s Public Science and Technology Research Funds Projects of Ocean (No. 200905013-4)by Ministry of Science and Technology of China (No. 2011BAC03B01)
文摘To investigate the mechanism of secondary circulations in rip current systems, and to explore the relationship between wave conditions and secondary circulation intensity, a series of numerical experiments is performed using coupled nearshore wave model and circulation model. In these experiments, the rip currents and secondary circulations generated above barred beaches with rip channels are simulated. A comparison experiment is conducted to investigate the formation and hydrodynamics of the secondary circulations. Model results indicate that the secondary circulations consist of alongshore flows driven by wave set-up near the shoreline, part of the feeder currents driven by the wave set-up over the bars, and onshore flows at the end of the rip channel driven by wave breaking and convection. The existence of the secondary circulation barely affects the rip current, but narrows and intensifies the feeder currents. Three groups of experiments of varying incident wave conditions are performed to investigate the relationship between wave conditions and secondary circulation intensity. The velocity of the alongshore flow of the secondary circulation is sensitive to the variation of the incident wave height and water depth. It is also found that the alongshore flow intensity is in direct proportion to the alongshore variation of the wave height gradient between the bars and the shoreline.
基金Project supported by the National Natural Science Foundation for Distinguished Young Scholars (Grant No. 50925625)the National Natural Science Foundation of China (Grant Nos. 10975026 and 10875024)the Scientific Research Foundation of the Higher Educational Institutions of Liaoning Education Bureau of China (Grant No. 2009A047)
文摘A preliminary investigation is conducted to study the characteristics of sheath damping near a dielectric wall with secondary electron emission (SEE). Making use of the linear analysis of the sheath stability, we have found two major contributions to the sheath damping, one similar to the Landau damping in uniform plasmas and another determined by local electric field and electron density of the steady-state sheath. It indicates that in a classical sheath regime the damping in the sheath region monotonically increases towards the wall and decreases with the enhancement of SEE effect. In order to verify the theoretical analysis, sheath oscillation processes induced by an initial disturbance are simulated with a time-dependent one-dimensional (1D) sheath model. Numerical results obtained are consistent with the theoretical analysis qualitatively.
基金Project supported by the National Natural Science Foundation of China (No.90205021), the Special Foundation for Doctoral Dissertations Research (No.200328) and the Science Foundation of Liuhui Center of Applied Mathematics, Nankai University and Tianjin University
文摘By direct numerical simulation, a supersonic boundary layer was studied to see whether the mechanism for the generation of sub-harmonic waves, similar to those for the incompressible flows, existed in the process of laminar-turbulent transition. The results showed that mechanisms of both resonant triad and secondary instability did exist. Discussions were made on whether these two mechanisms are really important in laminarturbulent transition.
基金This investigation was carried out as a continuation of the“Program for the Promotion of New Energy Infrastructure Development”,supported by the Mitsubishi Research Institute(MRI)/the Ministry of Economy,Trade and Industry(METI),Japan.
文摘A fixed oscillating water column(OWC)-type wave energy converter consists of an air chamber,an air turbine and a generator.The energy conversion processes are the primary conversion in an air chamber and the secondary conversion of the turbine.For the practical use,it is necessary to develop a design method which can consider the incident wave motion,the motion of the internal free surface affected in the structure such as a partly submerged wall,the fluctuation of air pressure in an air chamber,and the rotation of the air turbine.At here,the authors carried out the wave tank tests using the model OWC equipped with the impulse turbine and a generator to obtain the experimental data needed to make this design method.As the result,the efficiencies of the three cases with different speed ratio between generator and turbine,and the effects of the curtain wall depth and the wave length on the energy conversion performance were clarified.