To predict the wave loads of a flexible trimaran in different wave fields,a one-way interaction numerical simulation method is proposed by integrating the fluid solver(Star-CCM+)and structural solver(Abaqus).Differing...To predict the wave loads of a flexible trimaran in different wave fields,a one-way interaction numerical simulation method is proposed by integrating the fluid solver(Star-CCM+)and structural solver(Abaqus).Differing from the existing coupled CFD-FEA method for monohull ships in head waves,the presented method equates the mass and stiffness of the whole ship to the hull shell so that any transverse and longitudinal section stress of the hull in oblique waves can be obtained.Firstly,verification study and sensitivity analysis are carried out by comparing the trimaran motions using different mesh sizes and time step schemes.Discussion on the wave elevation of uni-and bi-directional waves is also carried out.Then a comprehensive analysis on the structural responses of the trimaran in different uni-directional regular wave and bi-directional cross sea conditions is carried out,respectively.Finally,the differences in structural response characteristics of trimaran in different wave fields are studied.The results show that the present method can reduce the computational burden of the two-way fluid-structure interaction simulations.展开更多
Based on the analysis of the advantages and disadvantages of some vertical coordinates applied in the calculation of the Changjiang diluted water (CDW), a new hybrid vertical coordinate is designed, which uses σ co...Based on the analysis of the advantages and disadvantages of some vertical coordinates applied in the calculation of the Changjiang diluted water (CDW), a new hybrid vertical coordinate is designed, which uses σ coordinate for current and σ-z coordinate for salinity. To combine the current and salinity, the Eulerian-Lagrangian method is used for the salinity calculation, and the baroclinic pressure gradient (BPG) is calculated on the salinity sited layers. The new hybrid vertical coordinate is introduced to the widely used model of POM (Princeton Ocean Model) to make a new model of POM-σ-z. The BPG calculations of an ideal case show that POM-σ-z model brings smaller error than POM model does. The simulations of CDW also show that POM-σ-z model is better than POM model on simulating the salinity and its front.展开更多
Typhoon-generated waves are simulated with two numerical wave models, the SWAN model for the coastal and Yangtze Estuary domain, nested within the WAVEWATCHIII (WW3) for the basin-scale East China Sea domain. Typhoo...Typhoon-generated waves are simulated with two numerical wave models, the SWAN model for the coastal and Yangtze Estuary domain, nested within the WAVEWATCHIII (WW3) for the basin-scale East China Sea domain. Typhoon No. 8114 is chosen because it was very strong, and generated high waves in the Estuary. WW3 was implemented for the East China Sea coarse-resolution computational domain, to simulate the waves over a large spatial scale and provide boundary conditions for SWAN model simulations, implemented on a fine-resolution nested domain for the Yangtze Estuary area. The Takahashi wind model is applied to the simulation of the East China Sea scale (3-hourly) and Yangtze Estuary scale (1-hourly) winds. Simulations of significant wave heights in the East China Sea show that the highest waves are on the right side of the storm track, and maxima tend to occur at the eastern deep-water open boundary of the Yangtze Estuary. In the Yangtze Estuary, incoming swell is dominant over locally generated waves before the typhoon approaches the Estuary. As the typhoon approaches the Estuary, wind waves and swell coexist, and the wave direction is mainly influenced by the swell direction and the complex topography.展开更多
The 4.2 ka event that occurred during the period from 4 500–3 900 a BP was characterized by cold and dry climates and resulted in the collapse of civilizations around the world. The cause of this climatic event, howe...The 4.2 ka event that occurred during the period from 4 500–3 900 a BP was characterized by cold and dry climates and resulted in the collapse of civilizations around the world. The cause of this climatic event, however, has been under debate. We collected four corals(Porites lutea) from Yongxing Island, Xisha Islands, South China Sea, dated them with the U-series method, and measured the annual coral growth rates using X-ray technology. The dating results showed that the coral growth ages were from 4 500–3 900 a BP, which coincide well with the period of the4.2 ka event. We then reconstructed annual sea surface temperature anomaly(SSTA) variations based on the coral growth rates. The growth rate-based SSTA results showed that the interdecadal SSTA from 4 500–3 900 a BP was lower than that during modern times(1961–2008 AD). A spectral analysis showed that the SSTA variations from4 500–3 900 a BP were under the influence of El Nino-Southern Oscillation(ENSO) activities. From 4 500–4 100 a BP, the climate exhibited La Nina-like conditions with weak ENSO intensity and relatively stable and lower SSTA amplitudes. From 4 100–3 900 a BP, the climate underwent a complicated period of ENSO variability and showed alternating El Nino-or La Nina-like conditions at interdecadal time scales and large SSTA amplitudes. We speculate that during the early and middle stages of the 4.2 ka event, the cold climate caused by weak ENSO activities largely weakened social productivity. Then, during the end stages of the 4.2 ka event, the repeated fluctuations in the ENSO intensity caused frequent extreme weather events, resulting in the collapse of civilizations worldwide. Thus, the new evidence obtained from our coral records suggests that the 4.2 ka event as well as the related collapse of civilizations were very likely driven by ENSO variability.展开更多
Two kinds of regression equations are used to reproduce the sediment flux of the 26 small coastal watersheds in southeastern China.The first kind is the global equations suggested by Milliman and Syvitski(1992),Mulder...Two kinds of regression equations are used to reproduce the sediment flux of the 26 small coastal watersheds in southeastern China.The first kind is the global equations suggested by Milliman and Syvitski(1992),Mulder and Syvitski(1996),Syvitski et al.(2003),and Syvitski and Milliman(2007).The second kind is the modified equations revised by the characteristics of the coastal watersheds,including the drainage area,mean water discharge,and mean sediment discharge.Compared with the observations of the hydrometric stations,the global equations overestimate the sediment flux by 1–2 orders of magnitude.By using the modified equations,the accuracy of the estimated sediment flux is significantly improved,with the relative error in the range of 7%–24%.The reason for the overestimation mainly caused by different parameters’domain and regression coefficients between global rivers and study coastal watersheds.This study demonstrates that modification needs to be considered when using global regression equations to reproduce the sediment flux of the small coastal watersheds in southeastern China.展开更多
In the preliminary design stage of the full form ships, in order to obtain a hull form with low resistance and maximum propulsion efficiency, an optimization design program for a full form ship with the minimum thrust...In the preliminary design stage of the full form ships, in order to obtain a hull form with low resistance and maximum propulsion efficiency, an optimization design program for a full form ship with the minimum thrust deduction factor has been developed, which combined the potential flow theory and boundary layer theory with the optimization technique. In the optimization process, the Sequential Unconstrained Minimization Technique(SUMT) interior point method of Nonlinear Programming(NLP) was proposed with the minimum thrust deduction factor as the objective function. An appropriate displacement is a basic constraint condition, and the boundary layer separation is an additional one. The parameters of the hull form modification function are used as design variables. At last, the numerical optimization example for lines of after-body of 50000 DWT product oil tanker was provided, which indicated that the propulsion efficiency was improved distinctly by this optimal design method.展开更多
On average,five to six storms occur in the Qiongzhou Strait every year,causing significant damage to coastal geomorphology and several property losses.Tropical Storm Bebinca is the most unusual and complex storm event...On average,five to six storms occur in the Qiongzhou Strait every year,causing significant damage to coastal geomorphology and several property losses.Tropical Storm Bebinca is the most unusual and complex storm event that has occurred in this region over the last 10 years.To detect the high-frequency beachface responses to the storm,a pressure sensor was deployed in the surf zone to record the free sea surface height,and the heights of grid pile points on the beachface were measured manually to determine beach elevation changes during this storm.Empirical Mode Decomposition and related analysis techniques were used to analyze the high-frequency topography and wave data.The results showed that:(1)the beachface response process occurred in three stages.The first stage was the rapid response stage,wherein the spring tide berm began to erode significantly,and the front edge of the beach berm reacted closely.The two beach sections resisted the harmful energy of the main storm.In the second stage,the beach slope increased after a large sediment loss on the beach berm and its front edge.To adapt to the storm energy,the beach at the low tide line began to erode,and the beach slope decreased.In the third stage,after the storm turned,the wave energy was significantly attenuated,and the beach berm eroded to resist the residual wave energy.The beachface began to oscillate and recover.(2)The main wave surface was the superimposed product of a few internal mode functions.Similar results were observed in beachface changes.High-frequency driving factors determine the local characteristics of beach evolution,and low-frequency driving factors determine the beach evolution trend.(3)The response of sediment to the storm was not a single sea-transportation,but a single-or two-way conversion driven by factors such as wave energy,swash flow,and secondary wave breaking.(4)TheΩ-RTR model is not completely applicable to beach states that undergo rapid changes during storms.Therefore,it is necessary to carry out further research on beach state identification during storms.展开更多
Much has been written of the error in computing the baroclinic pressure gradient (BPG) with sigma coordinates in ocean or atmospheric numerical models. The usual way to reduce the error is to subtract area-averaged de...Much has been written of the error in computing the baroclinic pressure gradient (BPG) with sigma coordinates in ocean or atmospheric numerical models. The usual way to reduce the error is to subtract area-averaged density stratification of the whole computation region. But if there is great difference between the area-averaged and the local averaged density stratification, the error will be obvious. An example is given to show that the error from this method may be larger than that from no correction sometimes. The definition of local area is put forward. Then, four improved BPG difference schemes of subtracting the local averaged density stratification are designed to reduce the error. Two of them are for diagnostic calculation (density field is fixed), and the others are for prognostic calculation (density field is not fixed). The results show that the errors from these schemes all significantly decrease.展开更多
The resonant motion of the fluid inside a narrow gap between two fixed boxes induced by incident regular waves with various wave heights is investigated by adopting a two-dimensional numerical wave flume based on an o...The resonant motion of the fluid inside a narrow gap between two fixed boxes induced by incident regular waves with various wave heights is investigated by adopting a two-dimensional numerical wave flume based on an opensourced CFD package,OpenFOAM.The two boxes have identical draft and height,but the upstream box has a variable breadth.This article focuses on the influences of the breadth ratio,defined as the ratio of the breadth of the upstream box to that of the downstream box,on the following three aspects of hydrodynamic characteristics of gap resonance:(1)the wave height amplifications inside the gap,and in front and at the rear of the structure system,(2)the reflection,transmission,and energy loss coefficients of the structure system,and(3)the response and damping time of the fluid resonance.It is found that the fluid resonant frequency,the amplification factor of the resonant wave height inside the gap and the maximum energy loss coefficient of the structure system are shown to gradually decrease with the increase of the breadth ratio.The response time of gap resonance is shown to first increase and then decrease with the breadth ratio overall,regardless of the incident wave height,and the configuration that the two boxes have the same breadth would bring the largest response time of gap resonance.展开更多
This paper explores an analytical model for Elastic Ring Squeeze Film Damper(ERSFD) with thin-walled ring and turbulent-jet orifices, and uncovers its Oil Film Pressure Performance(OFPP). Firstly, the ring deformation...This paper explores an analytical model for Elastic Ring Squeeze Film Damper(ERSFD) with thin-walled ring and turbulent-jet orifices, and uncovers its Oil Film Pressure Performance(OFPP). Firstly, the ring deformation is addressed by using the Fourier series expansion approach and the orifice outflow rate is characterized with the Prandtl boundary layer theory. Secondly, applying finite difference scheme, the influence of elastic ring flexibility, orifice diameter, and attitude angle on the OFPP is analyzed. Finally, Outer chamber pressure was measured experimentally at different rotor speeds. The results indicate that the outer chamber pressure coats an individual load-carrying region and spreads symmetrically pertaining to the attitude angle. Its amplitude drops as the elastic ring flexibility decreases but boosts with the reduction of the orifice diameter.For inner chamber pressure, the orifice diameter effects a similar trend to the outer cavity, but exhibits more stable distribution regarding the attitude angle. Minimizing the elastic ring flexibility causes an increase in amplitude. The model is validated by the test results giving that the outer chamber pressure shifts synchronously and periodically with the variation of the attitude angle,while the pressure amplitude increases slightly at higher rotor speeds.展开更多
基金financially supported by the State Key Laboratory of Structural Analysis,Optimization and CAE Software for Industrial Equipment,Dalian University of Technology(Grant No.GZ23112)the Shandong Provincial Natural Science Foundation,China(Grant No.ZR2021ME146).
文摘To predict the wave loads of a flexible trimaran in different wave fields,a one-way interaction numerical simulation method is proposed by integrating the fluid solver(Star-CCM+)and structural solver(Abaqus).Differing from the existing coupled CFD-FEA method for monohull ships in head waves,the presented method equates the mass and stiffness of the whole ship to the hull shell so that any transverse and longitudinal section stress of the hull in oblique waves can be obtained.Firstly,verification study and sensitivity analysis are carried out by comparing the trimaran motions using different mesh sizes and time step schemes.Discussion on the wave elevation of uni-and bi-directional waves is also carried out.Then a comprehensive analysis on the structural responses of the trimaran in different uni-directional regular wave and bi-directional cross sea conditions is carried out,respectively.Finally,the differences in structural response characteristics of trimaran in different wave fields are studied.The results show that the present method can reduce the computational burden of the two-way fluid-structure interaction simulations.
文摘Based on the analysis of the advantages and disadvantages of some vertical coordinates applied in the calculation of the Changjiang diluted water (CDW), a new hybrid vertical coordinate is designed, which uses σ coordinate for current and σ-z coordinate for salinity. To combine the current and salinity, the Eulerian-Lagrangian method is used for the salinity calculation, and the baroclinic pressure gradient (BPG) is calculated on the salinity sited layers. The new hybrid vertical coordinate is introduced to the widely used model of POM (Princeton Ocean Model) to make a new model of POM-σ-z. The BPG calculations of an ideal case show that POM-σ-z model brings smaller error than POM model does. The simulations of CDW also show that POM-σ-z model is better than POM model on simulating the salinity and its front.
基金This project is supported bythe Canadian Panel on Energy Research and Development (Offshore Environmental Fac-tors Program) , ONR (US Office of Naval Research) via GoMOOS-the Gulf of Maine Ocean Observing System,Petroleum Research Atlantic Canada (PRAC) ,and the CFCAS (Canada Foundation for Climate and AtmosphericStudies) ,Canadian Panel on Energy Research and Development (Offshore Environmental Factors Program) .It is al-so supported bythe Advanced Doctoral Fund of the Ministry of Education of China (Grant No.20030294010)
文摘Typhoon-generated waves are simulated with two numerical wave models, the SWAN model for the coastal and Yangtze Estuary domain, nested within the WAVEWATCHIII (WW3) for the basin-scale East China Sea domain. Typhoon No. 8114 is chosen because it was very strong, and generated high waves in the Estuary. WW3 was implemented for the East China Sea coarse-resolution computational domain, to simulate the waves over a large spatial scale and provide boundary conditions for SWAN model simulations, implemented on a fine-resolution nested domain for the Yangtze Estuary area. The Takahashi wind model is applied to the simulation of the East China Sea scale (3-hourly) and Yangtze Estuary scale (1-hourly) winds. Simulations of significant wave heights in the East China Sea show that the highest waves are on the right side of the storm track, and maxima tend to occur at the eastern deep-water open boundary of the Yangtze Estuary. In the Yangtze Estuary, incoming swell is dominant over locally generated waves before the typhoon approaches the Estuary. As the typhoon approaches the Estuary, wind waves and swell coexist, and the wave direction is mainly influenced by the swell direction and the complex topography.
基金The National Natural Science Foundation of China under contract No.91428203the Guangxi Scientific Projects under contract Nos AD17129063 and AA17204074the Bagui Fellowship from Guangxi of China
文摘The 4.2 ka event that occurred during the period from 4 500–3 900 a BP was characterized by cold and dry climates and resulted in the collapse of civilizations around the world. The cause of this climatic event, however, has been under debate. We collected four corals(Porites lutea) from Yongxing Island, Xisha Islands, South China Sea, dated them with the U-series method, and measured the annual coral growth rates using X-ray technology. The dating results showed that the coral growth ages were from 4 500–3 900 a BP, which coincide well with the period of the4.2 ka event. We then reconstructed annual sea surface temperature anomaly(SSTA) variations based on the coral growth rates. The growth rate-based SSTA results showed that the interdecadal SSTA from 4 500–3 900 a BP was lower than that during modern times(1961–2008 AD). A spectral analysis showed that the SSTA variations from4 500–3 900 a BP were under the influence of El Nino-Southern Oscillation(ENSO) activities. From 4 500–4 100 a BP, the climate exhibited La Nina-like conditions with weak ENSO intensity and relatively stable and lower SSTA amplitudes. From 4 100–3 900 a BP, the climate underwent a complicated period of ENSO variability and showed alternating El Nino-or La Nina-like conditions at interdecadal time scales and large SSTA amplitudes. We speculate that during the early and middle stages of the 4.2 ka event, the cold climate caused by weak ENSO activities largely weakened social productivity. Then, during the end stages of the 4.2 ka event, the repeated fluctuations in the ENSO intensity caused frequent extreme weather events, resulting in the collapse of civilizations worldwide. Thus, the new evidence obtained from our coral records suggests that the 4.2 ka event as well as the related collapse of civilizations were very likely driven by ENSO variability.
基金The National Natural Science Foundation of China under contract Nos 41625021,41676079 and 41906021the Project of Enhancing School with Innovation of Guangdong Ocean University under contract No.Q18307the Program for Scientific Research Start-up Funds of Guangdong Ocean University under contract No.060302112010
文摘Two kinds of regression equations are used to reproduce the sediment flux of the 26 small coastal watersheds in southeastern China.The first kind is the global equations suggested by Milliman and Syvitski(1992),Mulder and Syvitski(1996),Syvitski et al.(2003),and Syvitski and Milliman(2007).The second kind is the modified equations revised by the characteristics of the coastal watersheds,including the drainage area,mean water discharge,and mean sediment discharge.Compared with the observations of the hydrometric stations,the global equations overestimate the sediment flux by 1–2 orders of magnitude.By using the modified equations,the accuracy of the estimated sediment flux is significantly improved,with the relative error in the range of 7%–24%.The reason for the overestimation mainly caused by different parameters’domain and regression coefficients between global rivers and study coastal watersheds.This study demonstrates that modification needs to be considered when using global regression equations to reproduce the sediment flux of the small coastal watersheds in southeastern China.
基金financially supported by the National Natural Science Foundation of China(Grant No.51009087)
文摘In the preliminary design stage of the full form ships, in order to obtain a hull form with low resistance and maximum propulsion efficiency, an optimization design program for a full form ship with the minimum thrust deduction factor has been developed, which combined the potential flow theory and boundary layer theory with the optimization technique. In the optimization process, the Sequential Unconstrained Minimization Technique(SUMT) interior point method of Nonlinear Programming(NLP) was proposed with the minimum thrust deduction factor as the objective function. An appropriate displacement is a basic constraint condition, and the boundary layer separation is an additional one. The parameters of the hull form modification function are used as design variables. At last, the numerical optimization example for lines of after-body of 50000 DWT product oil tanker was provided, which indicated that the propulsion efficiency was improved distinctly by this optimal design method.
基金The National Natural Science Foundation of China under contract Nos 42176167 and 41676079the Project of Enhancing School with Innovation,Guangdong Ocean University under contract No.Q18307.
文摘On average,five to six storms occur in the Qiongzhou Strait every year,causing significant damage to coastal geomorphology and several property losses.Tropical Storm Bebinca is the most unusual and complex storm event that has occurred in this region over the last 10 years.To detect the high-frequency beachface responses to the storm,a pressure sensor was deployed in the surf zone to record the free sea surface height,and the heights of grid pile points on the beachface were measured manually to determine beach elevation changes during this storm.Empirical Mode Decomposition and related analysis techniques were used to analyze the high-frequency topography and wave data.The results showed that:(1)the beachface response process occurred in three stages.The first stage was the rapid response stage,wherein the spring tide berm began to erode significantly,and the front edge of the beach berm reacted closely.The two beach sections resisted the harmful energy of the main storm.In the second stage,the beach slope increased after a large sediment loss on the beach berm and its front edge.To adapt to the storm energy,the beach at the low tide line began to erode,and the beach slope decreased.In the third stage,after the storm turned,the wave energy was significantly attenuated,and the beach berm eroded to resist the residual wave energy.The beachface began to oscillate and recover.(2)The main wave surface was the superimposed product of a few internal mode functions.Similar results were observed in beachface changes.High-frequency driving factors determine the local characteristics of beach evolution,and low-frequency driving factors determine the beach evolution trend.(3)The response of sediment to the storm was not a single sea-transportation,but a single-or two-way conversion driven by factors such as wave energy,swash flow,and secondary wave breaking.(4)TheΩ-RTR model is not completely applicable to beach states that undergo rapid changes during storms.Therefore,it is necessary to carry out further research on beach state identification during storms.
基金The Major State Basic Research Program of China under contract No. 2002412403the National Natural Science Foundation of China un-der contract No. 40306014.
文摘Much has been written of the error in computing the baroclinic pressure gradient (BPG) with sigma coordinates in ocean or atmospheric numerical models. The usual way to reduce the error is to subtract area-averaged density stratification of the whole computation region. But if there is great difference between the area-averaged and the local averaged density stratification, the error will be obvious. An example is given to show that the error from this method may be larger than that from no correction sometimes. The definition of local area is put forward. Then, four improved BPG difference schemes of subtracting the local averaged density stratification are designed to reduce the error. Two of them are for diagnostic calculation (density field is fixed), and the others are for prognostic calculation (density field is not fixed). The results show that the errors from these schemes all significantly decrease.
基金financially supported by the National Natural Science Foundation of China (Grant Nos.52371277 and 52001071)the State Key Laboratory of Hydraulic Engineering Intelligent Construction and Operation (Tianjin University) (Grant No.HESS-2323)+1 种基金the Postgraduate Research&Practice Innovation Program of Jiangsu Province (Grant No.KYCX23_3902)the Guangdong Basic and Applied Basic Research Foundation (Grant No.:2023A1515010890)。
文摘The resonant motion of the fluid inside a narrow gap between two fixed boxes induced by incident regular waves with various wave heights is investigated by adopting a two-dimensional numerical wave flume based on an opensourced CFD package,OpenFOAM.The two boxes have identical draft and height,but the upstream box has a variable breadth.This article focuses on the influences of the breadth ratio,defined as the ratio of the breadth of the upstream box to that of the downstream box,on the following three aspects of hydrodynamic characteristics of gap resonance:(1)the wave height amplifications inside the gap,and in front and at the rear of the structure system,(2)the reflection,transmission,and energy loss coefficients of the structure system,and(3)the response and damping time of the fluid resonance.It is found that the fluid resonant frequency,the amplification factor of the resonant wave height inside the gap and the maximum energy loss coefficient of the structure system are shown to gradually decrease with the increase of the breadth ratio.The response time of gap resonance is shown to first increase and then decrease with the breadth ratio overall,regardless of the incident wave height,and the configuration that the two boxes have the same breadth would bring the largest response time of gap resonance.
基金supported by the National Natural Science Foundation of China(No.52005158)。
文摘This paper explores an analytical model for Elastic Ring Squeeze Film Damper(ERSFD) with thin-walled ring and turbulent-jet orifices, and uncovers its Oil Film Pressure Performance(OFPP). Firstly, the ring deformation is addressed by using the Fourier series expansion approach and the orifice outflow rate is characterized with the Prandtl boundary layer theory. Secondly, applying finite difference scheme, the influence of elastic ring flexibility, orifice diameter, and attitude angle on the OFPP is analyzed. Finally, Outer chamber pressure was measured experimentally at different rotor speeds. The results indicate that the outer chamber pressure coats an individual load-carrying region and spreads symmetrically pertaining to the attitude angle. Its amplitude drops as the elastic ring flexibility decreases but boosts with the reduction of the orifice diameter.For inner chamber pressure, the orifice diameter effects a similar trend to the outer cavity, but exhibits more stable distribution regarding the attitude angle. Minimizing the elastic ring flexibility causes an increase in amplitude. The model is validated by the test results giving that the outer chamber pressure shifts synchronously and periodically with the variation of the attitude angle,while the pressure amplitude increases slightly at higher rotor speeds.