采用耦合器组件形式的耦合技术现已成为数值模拟中不同模式物理过程间实现复杂相互作用的主流方法。以美国阿尔贡国家实验室(Argonne National Laboratory)开发的模式耦合工具包(MCT)为基础构建了耦合器组件,并利用其实现了海冰模式CSIM...采用耦合器组件形式的耦合技术现已成为数值模拟中不同模式物理过程间实现复杂相互作用的主流方法。以美国阿尔贡国家实验室(Argonne National Laboratory)开发的模式耦合工具包(MCT)为基础构建了耦合器组件,并利用其实现了海冰模式CSIM5与区域海洋环流模式ROMS3在北极地区的耦合。耦合模式以多任务多数据(MPMD)并行方式运行,在采用Linux系统和pgi编译器的高性能微机及银河集群计算机系统下已实现利用其开展冰—洋相互作用数值模拟研究的功能。展开更多
Hydraulic butterfly valves have been widely applied in marine engineering because of their large switching torque, low pressure loss and suitability for large and medium diameter pipelines. Due to control problems res...Hydraulic butterfly valves have been widely applied in marine engineering because of their large switching torque, low pressure loss and suitability for large and medium diameter pipelines. Due to control problems resulting from switching angular speeds of the hydraulic butterfly valve, a throttle-governing control mode has been widely adopted, and detailed analysis has been carried out worldwide on the structural principle concerning speed-regulation and the load torque on the shaft while opening or closing a hydraulic butterfly valve. However relevant reports have yet been published on the change law, the error and the influencing factors of the rotational angular velocity of the hydraulic butterfly valve while opening and closing. In this article, research was based on some common specifications of a hydraulic butterfly valve with a symmetrical valve flap existing in a marine environment. The throttle governing system supplied by the accumulator to achieve the switching of the hydraulic control valve was adopted, and the mathematical models of the system were established in the actual conditions while the numerical simulations took place. The simulation results and analysis show that the rotational angular velocity and the error of the hydraulic butterfly valve while switching is influenced greatly by the drainage amount of the accumulator, resulting in pressure loss in the pipeline, the temperature of hydraulic medium and the load of the hydraulic butterfly valve. The simulation results and analysis provide a theoretical basis for the choice of the total capacity of the accumulator and pipeline diameters in a throttle governing system with a hydraulic butterfly valve.It also determines the type and specification of the hydraulic butterfly valve and the design of motion parameters of the transported fluid.展开更多
The South China Sea contains tremendous oil and gas resources in deepwater areas. However, one of the keys for deepwater exploration, the investigation of deepwater floating platforms, is very inadequate. In this pape...The South China Sea contains tremendous oil and gas resources in deepwater areas. However, one of the keys for deepwater exploration, the investigation of deepwater floating platforms, is very inadequate. In this paper, the authors studied and compared the hydrodynamics and global motion behaviors of typical deepwater platforms in the South China Sea environment. The hydrodynamic models of three main types of floating platforms, e.g. the Semi-submersible, tension leg platform （TLP）, and Truss Spar, which could potentially be utilized in the South China Sea, were established by using the 3-D potential theory. Additionally, some important considerations which significantly influence the hydrodynamics were given. The RAOs in frequency domains as well as global motions in time domains under time-varying wind, random waves, and current in 100-y, 10-y, and 1-y return period environment conditions were predicted, compared, and analyzed. The results indicate that the heave and especially the pitch motion of the TLP are favorable. The heave response of the Truss Spar is perfect and comparable with that of the TLP when the peak period of random waves is low. However, the pitch motion of Truss Spar is extraordinarily lar^er than that of Semi-submersible and TLP.展开更多
Using the latest version of Mesoscale Modeling System （MMSv3）, we assimilated wind data from the scatterometer and built a model to assimilate the wind field over eastern China seas and adjacent waters and applied t...Using the latest version of Mesoscale Modeling System （MMSv3）, we assimilated wind data from the scatterometer and built a model to assimilate the wind field over eastern China seas and adjacent waters and applied the wave model WAVEWATCH-Ⅲ to test the sea area with assimilative wind and blended wind of QSCAT and NCEP as driving forces. High precision and resolution numerical wave results were obtained. Analysis indicated that if we replace the model wind result with the blended wind, better sea surface wind results and wave results could be obtained.展开更多
A non-hydrostatic, Boussinesq, and three-dimensional large eddy simulation(LES) model was used to study the impact of the Earth's rotation on turbulence and the redistribution of energy in turbulence kinetic energ...A non-hydrostatic, Boussinesq, and three-dimensional large eddy simulation(LES) model was used to study the impact of the Earth's rotation on turbulence and the redistribution of energy in turbulence kinetic energy(TKE) budget. A set of numerical simulations was conducted,(1) with and without rotation,(2) at different latitudes(10°N, 30°N, 45°N, 60°N, and 80°N),(3) with wave breaking and with Langmuir circulation, and(4) under different wind speeds(5, 10, 20, and 30 m/s). The results show that eddy viscosity decreases when rotation is included, indicating that rotation weakens the turbulence strength. The TKE budget become tight with rotation and the effects of rotation grow with latitude. However, rotation become less important under Langmuir circulation since the transport term is strong in the vertical direction. Finally, simulations were conducted based on field data from the Boundary Layer and Air-Sea Transfer Low Wind(CBLAST-Low) experiment. The results, although more complex, are consistent with the results obtained from earlier simulations using ideal numerical conditions.展开更多
Based on the Fourier transform, the analytical solution of boundary integral equations formulated for the complex velocity of a 2-D steady linear surface flow is derived. It has been found that before the radiation co...Based on the Fourier transform, the analytical solution of boundary integral equations formulated for the complex velocity of a 2-D steady linear surface flow is derived. It has been found that before the radiation condition is imposed,free waves appear both far upstream and downstream. In order to cancel the free waves in far upstream regions, the eigensolution of a specific eigenvalue, which satisfies the homogeneous boundary integral equation, is found and superposed to the analytical solution. An example, a submerged vortex, is used to demonstrate the derived analytical solution. Furthermore,an analytical approach to imposing the radiation condition in the numerical solution of boundary integral equations for 2-D steady linear wave problems is proposed.展开更多
In this paper, we discussed the features of atmospheric circulations over Eurasia as a response to sea surface temperature anomalies （SSTAs） over the tropical Indian Ocean, the equatorial Pacific, Kuroshio and the N...In this paper, we discussed the features of atmospheric circulations over Eurasia as a response to sea surface temperature anomalies （SSTAs） over the tropical Indian Ocean, the equatorial Pacific, Kuroshio and the North Atlantic. Our results are shown as follows： （1） CAM3.0, driven by the combined SSTAs over the four oceanic regions, can simulate well the features of anomalous atmospheric circulations over Eurasia in January 2008, indicating that the effects of the SSTAs over these four regions were one of the key causes of the anomalous systems over Eurasia. （2） The SSTAs over each key region contributed to the intensification of blocking over the Urals Mountains and a main East Asian trough. However, the influence of the SSTAs over individual oceanic regions differed from one another in other aspects. The SSTAs over the North Atlantic had an impact on the 500-hPa anomalous height （Z500A） over the middle-high latitudes and had a somewhat smaller effect over the low latitudes. For the warm SSTAs over Kuroshio, the subtropical high was much stronger, spread farther north than usual, and had an anomalous easterly that dominated the northwest Pacific Ocean. The warm SSTAs over the tropical Indian Ocean could have caused a negative Z500A from West Asia to Middle Asia, a remarkably anomalous southwesterly from the Indian Ocean to the south of China and an anomalous anticyclone circulation over the South China Sea-Philippine Sea region. Because of the La Nifia event, the winter monsoon was stronger than normal, with an anomalously cooler northerly over the southeastern coastal areas of China. （3） The combined effects of the SSTAs over the four key regions were likely more important to the atmospheric circulation anomalies of January 2008 over Eurasia than the effects of individual or partly combined SSTAS. This unique SSTA distribution possibly led to the circulation anomalies over Eurasia in January 2008, especially the atmospheric circulation anomalies over the subtropics, which were more similar to those of the winter E1 Nifio events than to the circulation anomalies following La Nifia.展开更多
The objective of the present investigation is to study the vortex-induced vibrations （VIV） for flow past a circular cylinder. The turbulent flow is simulated by using a 2-D standard k-ε model incorporating the fini...The objective of the present investigation is to study the vortex-induced vibrations （VIV） for flow past a circular cylinder. The turbulent flow is simulated by using a 2-D standard k-ε model incorporating the finite volume method （FVM） and the Semi-Implicit Method for the Pressure Linked Equations （SIMPLE） algorithm on non-orthogonal boundary-fitted collocated grids. The wall boundaries are approximated with wall functions. In the numerical cases, the turbulent wake patterns are studied by plotting the streamlines and the turbulent kinetic energy contours. The pressure distributions are investigated. Analyses of the vortex-induced force coefficients and the structural vibrations are carried out. The variations of the Strouhal number with the Reynolds number and of the vortex-induced force coefficients with the reduced velocity are obtained. The results show that this numerical approach is feasible and efficient in investigating the VIV problem for a circular cylinder.展开更多
Long-term evolution of the Black Sea dynamics （1980-2020） is reconstructed by means of numerical simulation. The model of the Black Sea circulation is z-coordinate model with 4.8 km horizontal space resolution and 4...Long-term evolution of the Black Sea dynamics （1980-2020） is reconstructed by means of numerical simulation. The model of the Black Sea circulation is z-coordinate model with 4.8 km horizontal space resolution and 40 levels in vertical direction. Mixing processes in the upper layer are parameterized with the Mellor-Yamada turbulent model. As for the boundary conditions on the sea surface, we used atmospheric forcing functions for the Black Sea region provided by CMCC using regional climate model COSMO-CLM. These data have a spatial resolution of 14km and a daily temporal resolution. To evaluate the quality of the Black Sea circulation dynamics, derived from simulation, the modeling results are compared with results of the Black Sea physical reanalysis. This reanalysis was performed by assimilating the temperature and salinity profiles from hydrographic surveys conducted during 1971-1993.展开更多
Recent satellite altimeter observations have indicated that internal tides （ITs） from the Luzon Strait （LS） propagate more than 2 500 km into the Western Pacific （WP）. This study used a high-resolution three-dim...Recent satellite altimeter observations have indicated that internal tides （ITs） from the Luzon Strait （LS） propagate more than 2 500 km into the Western Pacific （WP）. This study used a high-resolution three-dimensional numerical model to reproduce and examine the ITs radiation process. The propagation of diurnal and semidiurnal ITs showed different patterns and variations. Diurnal ITs with lower frequency were affected more by the earth＇s rotation and they were bent more toward the equator than semidiurnal ITs. ITs phase speeds are functions of latitude and diurnal ITs showed greater distinctions of phase speeds during propagation. For M2 ITs, the wavelength remained nearly unchanged but the beam width increased significantly during propagation away from the LS. For diurnal ITs （K1 and O0, the wavelength decreased noticeably with latitude, while the beam width varied little during propagation because of blocking by land. Baroclinic energy was also examined as a complement to satellite results reported by Zhao （2014）. The magnitude of the generated baroclinic energy flux reduced remarkably within 300 km from the generation site but it then decayed slowly when propagating into abyssal sea. Baroclinic energy of diurnal ITs was found to dissipate at a slower rate than semidiurnal ITs along the main propagation path in the WP.展开更多
After the damage of Fukushima Daiichi Nuclear Power Plant, a great number of radioactive materials were released into the Pacific Ocean. Therefore, it is necessary to research on the temporal and spatial distribution ...After the damage of Fukushima Daiichi Nuclear Power Plant, a great number of radioactive materials were released into the Pacific Ocean. Therefore, it is necessary to research on the temporal and spatial distribution of these radionuclides. We use Princeton Ocean Model to simulate the circulation of the coast water of Fukushima NPP and obtain the concentration of caesium-137 by solving the diffusion equations. We employ the Monte Carlo N-particle(MCNP) code to assess the external doses caused by these contaminated sea water. To improve the efficiency and effectiveness of volume source in MCNP code, we establish a transformation method between spot source and volume source, and determine an appropriate range of volume source. Finally, we calculate the absorbed doses of every organ/tissue and the effective dose of a human body.展开更多
文摘采用耦合器组件形式的耦合技术现已成为数值模拟中不同模式物理过程间实现复杂相互作用的主流方法。以美国阿尔贡国家实验室(Argonne National Laboratory)开发的模式耦合工具包(MCT)为基础构建了耦合器组件,并利用其实现了海冰模式CSIM5与区域海洋环流模式ROMS3在北极地区的耦合。耦合模式以多任务多数据(MPMD)并行方式运行,在采用Linux系统和pgi编译器的高性能微机及银河集群计算机系统下已实现利用其开展冰—洋相互作用数值模拟研究的功能。
文摘Hydraulic butterfly valves have been widely applied in marine engineering because of their large switching torque, low pressure loss and suitability for large and medium diameter pipelines. Due to control problems resulting from switching angular speeds of the hydraulic butterfly valve, a throttle-governing control mode has been widely adopted, and detailed analysis has been carried out worldwide on the structural principle concerning speed-regulation and the load torque on the shaft while opening or closing a hydraulic butterfly valve. However relevant reports have yet been published on the change law, the error and the influencing factors of the rotational angular velocity of the hydraulic butterfly valve while opening and closing. In this article, research was based on some common specifications of a hydraulic butterfly valve with a symmetrical valve flap existing in a marine environment. The throttle governing system supplied by the accumulator to achieve the switching of the hydraulic control valve was adopted, and the mathematical models of the system were established in the actual conditions while the numerical simulations took place. The simulation results and analysis show that the rotational angular velocity and the error of the hydraulic butterfly valve while switching is influenced greatly by the drainage amount of the accumulator, resulting in pressure loss in the pipeline, the temperature of hydraulic medium and the load of the hydraulic butterfly valve. The simulation results and analysis provide a theoretical basis for the choice of the total capacity of the accumulator and pipeline diameters in a throttle governing system with a hydraulic butterfly valve.It also determines the type and specification of the hydraulic butterfly valve and the design of motion parameters of the transported fluid.
基金Supported by the National Sci-Tech Major Special Item(No.2008ZX05056-03)
文摘The South China Sea contains tremendous oil and gas resources in deepwater areas. However, one of the keys for deepwater exploration, the investigation of deepwater floating platforms, is very inadequate. In this paper, the authors studied and compared the hydrodynamics and global motion behaviors of typical deepwater platforms in the South China Sea environment. The hydrodynamic models of three main types of floating platforms, e.g. the Semi-submersible, tension leg platform （TLP）, and Truss Spar, which could potentially be utilized in the South China Sea, were established by using the 3-D potential theory. Additionally, some important considerations which significantly influence the hydrodynamics were given. The RAOs in frequency domains as well as global motions in time domains under time-varying wind, random waves, and current in 100-y, 10-y, and 1-y return period environment conditions were predicted, compared, and analyzed. The results indicate that the heave and especially the pitch motion of the TLP are favorable. The heave response of the Truss Spar is perfect and comparable with that of the TLP when the peak period of random waves is low. However, the pitch motion of Truss Spar is extraordinarily lar^er than that of Semi-submersible and TLP.
基金Supported by the High-Tech Research and Development Program of China (863 Program, No. 2001AA633070 2003AA604040)the National Natural Science Foundation of China (No. 40476015).
文摘Using the latest version of Mesoscale Modeling System （MMSv3）, we assimilated wind data from the scatterometer and built a model to assimilate the wind field over eastern China seas and adjacent waters and applied the wave model WAVEWATCH-Ⅲ to test the sea area with assimilative wind and blended wind of QSCAT and NCEP as driving forces. High precision and resolution numerical wave results were obtained. Analysis indicated that if we replace the model wind result with the blended wind, better sea surface wind results and wave results could be obtained.
基金Supported by the National Natural Science Foundation of China(Nos.41206015,41106019)the National Basic Research Program of China(973 Program)(Nos.2011CB403501,2012CB417402)the Fund for Creative Research Groups by NSFC(No.41121064)
文摘A non-hydrostatic, Boussinesq, and three-dimensional large eddy simulation(LES) model was used to study the impact of the Earth's rotation on turbulence and the redistribution of energy in turbulence kinetic energy(TKE) budget. A set of numerical simulations was conducted,(1) with and without rotation,(2) at different latitudes(10°N, 30°N, 45°N, 60°N, and 80°N),(3) with wave breaking and with Langmuir circulation, and(4) under different wind speeds(5, 10, 20, and 30 m/s). The results show that eddy viscosity decreases when rotation is included, indicating that rotation weakens the turbulence strength. The TKE budget become tight with rotation and the effects of rotation grow with latitude. However, rotation become less important under Langmuir circulation since the transport term is strong in the vertical direction. Finally, simulations were conducted based on field data from the Boundary Layer and Air-Sea Transfer Low Wind(CBLAST-Low) experiment. The results, although more complex, are consistent with the results obtained from earlier simulations using ideal numerical conditions.
文摘Based on the Fourier transform, the analytical solution of boundary integral equations formulated for the complex velocity of a 2-D steady linear surface flow is derived. It has been found that before the radiation condition is imposed,free waves appear both far upstream and downstream. In order to cancel the free waves in far upstream regions, the eigensolution of a specific eigenvalue, which satisfies the homogeneous boundary integral equation, is found and superposed to the analytical solution. An example, a submerged vortex, is used to demonstrate the derived analytical solution. Furthermore,an analytical approach to imposing the radiation condition in the numerical solution of boundary integral equations for 2-D steady linear wave problems is proposed.
基金Special Funds for Public Welfare of China (GYHY(QX) 2008-06-005)Science and Technology Innovation Program of Jiangsu Province (CX09B_221Z)
文摘In this paper, we discussed the features of atmospheric circulations over Eurasia as a response to sea surface temperature anomalies （SSTAs） over the tropical Indian Ocean, the equatorial Pacific, Kuroshio and the North Atlantic. Our results are shown as follows： （1） CAM3.0, driven by the combined SSTAs over the four oceanic regions, can simulate well the features of anomalous atmospheric circulations over Eurasia in January 2008, indicating that the effects of the SSTAs over these four regions were one of the key causes of the anomalous systems over Eurasia. （2） The SSTAs over each key region contributed to the intensification of blocking over the Urals Mountains and a main East Asian trough. However, the influence of the SSTAs over individual oceanic regions differed from one another in other aspects. The SSTAs over the North Atlantic had an impact on the 500-hPa anomalous height （Z500A） over the middle-high latitudes and had a somewhat smaller effect over the low latitudes. For the warm SSTAs over Kuroshio, the subtropical high was much stronger, spread farther north than usual, and had an anomalous easterly that dominated the northwest Pacific Ocean. The warm SSTAs over the tropical Indian Ocean could have caused a negative Z500A from West Asia to Middle Asia, a remarkably anomalous southwesterly from the Indian Ocean to the south of China and an anomalous anticyclone circulation over the South China Sea-Philippine Sea region. Because of the La Nifia event, the winter monsoon was stronger than normal, with an anomalously cooler northerly over the southeastern coastal areas of China. （3） The combined effects of the SSTAs over the four key regions were likely more important to the atmospheric circulation anomalies of January 2008 over Eurasia than the effects of individual or partly combined SSTAS. This unique SSTA distribution possibly led to the circulation anomalies over Eurasia in January 2008, especially the atmospheric circulation anomalies over the subtropics, which were more similar to those of the winter E1 Nifio events than to the circulation anomalies following La Nifia.
基金The present research is financially supported by the National Natural Science Foundation of China(No.50479027)the Natural Science Foundation of Qingdao,China(No.04-3-jj-01).The support is gratefully appreciated.
文摘The objective of the present investigation is to study the vortex-induced vibrations （VIV） for flow past a circular cylinder. The turbulent flow is simulated by using a 2-D standard k-ε model incorporating the finite volume method （FVM） and the Semi-Implicit Method for the Pressure Linked Equations （SIMPLE） algorithm on non-orthogonal boundary-fitted collocated grids. The wall boundaries are approximated with wall functions. In the numerical cases, the turbulent wake patterns are studied by plotting the streamlines and the turbulent kinetic energy contours. The pressure distributions are investigated. Analyses of the vortex-induced force coefficients and the structural vibrations are carried out. The variations of the Strouhal number with the Reynolds number and of the vortex-induced force coefficients with the reduced velocity are obtained. The results show that this numerical approach is feasible and efficient in investigating the VIV problem for a circular cylinder.
文摘Long-term evolution of the Black Sea dynamics （1980-2020） is reconstructed by means of numerical simulation. The model of the Black Sea circulation is z-coordinate model with 4.8 km horizontal space resolution and 40 levels in vertical direction. Mixing processes in the upper layer are parameterized with the Mellor-Yamada turbulent model. As for the boundary conditions on the sea surface, we used atmospheric forcing functions for the Black Sea region provided by CMCC using regional climate model COSMO-CLM. These data have a spatial resolution of 14km and a daily temporal resolution. To evaluate the quality of the Black Sea circulation dynamics, derived from simulation, the modeling results are compared with results of the Black Sea physical reanalysis. This reanalysis was performed by assimilating the temperature and salinity profiles from hydrographic surveys conducted during 1971-1993.
基金Supported by the National Natural Science Foundation of China(Nos.41528601,41376029,U1406401,41421005)the Strategic Pioneering Research Program of CAS(Nos.XDA10020104,XDA10020101)the CAS Interdisciplinary Innovation Team“Ocean Mesoscale Dynamical Processes and ecological effect”
文摘Recent satellite altimeter observations have indicated that internal tides （ITs） from the Luzon Strait （LS） propagate more than 2 500 km into the Western Pacific （WP）. This study used a high-resolution three-dimensional numerical model to reproduce and examine the ITs radiation process. The propagation of diurnal and semidiurnal ITs showed different patterns and variations. Diurnal ITs with lower frequency were affected more by the earth＇s rotation and they were bent more toward the equator than semidiurnal ITs. ITs phase speeds are functions of latitude and diurnal ITs showed greater distinctions of phase speeds during propagation. For M2 ITs, the wavelength remained nearly unchanged but the beam width increased significantly during propagation away from the LS. For diurnal ITs （K1 and O0, the wavelength decreased noticeably with latitude, while the beam width varied little during propagation because of blocking by land. Baroclinic energy was also examined as a complement to satellite results reported by Zhao （2014）. The magnitude of the generated baroclinic energy flux reduced remarkably within 300 km from the generation site but it then decayed slowly when propagating into abyssal sea. Baroclinic energy of diurnal ITs was found to dissipate at a slower rate than semidiurnal ITs along the main propagation path in the WP.
基金supported by the National Natural Science Foundation of China(Grant No.71373140)the State Key Laboratory of NBC Protection for Civilian(Grant No.SKLNBC0308)Tsinghua University Initiative Scientific Research Program(Grant Nos.2012Z10137,2012THZ0124)
文摘After the damage of Fukushima Daiichi Nuclear Power Plant, a great number of radioactive materials were released into the Pacific Ocean. Therefore, it is necessary to research on the temporal and spatial distribution of these radionuclides. We use Princeton Ocean Model to simulate the circulation of the coast water of Fukushima NPP and obtain the concentration of caesium-137 by solving the diffusion equations. We employ the Monte Carlo N-particle(MCNP) code to assess the external doses caused by these contaminated sea water. To improve the efficiency and effectiveness of volume source in MCNP code, we establish a transformation method between spot source and volume source, and determine an appropriate range of volume source. Finally, we calculate the absorbed doses of every organ/tissue and the effective dose of a human body.
