砾石单元法与SPH耦合模型在数值波浪水槽中的应用研究

为研究以流体粒子描述波浪运动,以固体单元描述砾石运动的两相介质大变形运动,在港口、海岸工程科学研究中具有重要意义。本文提出砾石单元法（GEM）,介绍了光滑粒子动力学方法（SPH）和GEM的基本原理,阐述了GEM与离散单元法（DEM）的异同之处,说明了采用SPH方法与GEM构建波浪砾石耦合运动数学模型的方法和过程。应用SPH方法建立数值波浪水槽,用GEM模拟波浪作用下堆积砾石的滚落、坍塌变形,构建了SPH方法与GEM耦合数学模型。模拟了水槽造波和波浪生成过程和波浪作用下砾石的滚落、坍塌变形,并与物理模型试验成果进行了比较,结果基本吻合。本文提出的GEM法具有模拟单相堆积砾石运动和堆积砾石与流体粒子耦合多相介质运动的功能,是对DEM法的补充和改善。本文提出的堆积力学球概念和拟序排列求解方法是砾石单元法的重要组成部分。     

波浪对透空式结构物冲击作用的光滑粒子流体动力学数值模拟

基于光滑粒子流体动力学(SPH)方法的二维数值波浪水槽模拟了规则波对透空式结构物的冲击作用。通过黎曼解和CSPM相结合的方法对连续方程和动量方程进行了修正。对造波边界采用虚粒子法模拟。提出了一种耦合计算方法来消除造波边界附近的压力波动现象,在结构物边界处设置了适合的冲击边界条件。应用修正的SPH模型模拟了规则波对浪溅区结构物的冲击作用,给出了结构物附近流场和压力场的变化特征,并应用物理模型试验结果对数值模型进行了验证。     

波浪与大孔隙多孔介质相互作用的耦合数学模型

建立了波浪与大孔隙多孔介质相互作用的耦合数学模型,波浪域的控制方程为雷诺时均方程和k-ε紊流模型。对于计算域的入射波采用推板式造波,它可以是线性波、椭圆余弦波和孤立波。采用PLIC-VOF法追踪波浪自由表面。对于多孔介质内的孔隙流场采用非线性Forchheimer方程,两区域共享连续方程,最后导出的波浪域与孔隙流域的压力修正方程具有完全相同的形式,利用这个方程能够同时而不是分别求解波浪场和孔隙流场,避免了在内部边界上给定匹配条件,实现了波浪场与孔隙流场的同步耦合。波浪与粗颗粒海床、平底床面上抛石潜堤及斜坡上抛石潜堤相互作用的验证计算结果表明该模型可用于研究波浪与大孔隙多孔介质相互作用的问题。     

波浪作用下锚链系泊浮式防波堤动力响应的数值模型研究

浮式防波堤充分利用波能在水深方向的分布特性,在满足工程消浪要求的同时对海域水沙交换影响较小,且能够快速布置,在某些实际工程有一定应用前景。为了深入了解波浪作用下浮式防波堤的动力响应,基于OpenFOAM标准求解器olaFlow,在刚体运动求解计算中植入锚链求解模块MOODY(mooring cable dynamics),实现了基于重叠网格方法的浮体运动与锚链受力耦合求解,建立了锚链系泊浮式防波堤动力响应的二维数值模型。利用该数值模型对锚链系泊单方箱浮式防波堤在波浪作用下的透射系数、运动响应、锚链张力进行了模拟,并和相关试验结果进行了比较。结果表明,模型能够准确模拟二维波浪和浮式防波堤的相互作用,并用于三维模型的改进。     

波浪与多孔介质结构相互作用SPH模拟

应用经CSPM法和黎曼解修正后的光滑粒子流体动力学(SPH)方法,建立了主动吸收无反射数值波浪水槽,研究波浪作用下多孔介质结构的水动力特性。流体运动控制方程采用N-S方程,多孔渗水结构内流体的运动控制方程考虑渗流力的影响。数值计算结果给出了水槽内不同位置测点的波面历时曲线和越浪量随时间变化曲线,并同试验结果和Philip Liu的数值计算结果进行了比较。并对一个波浪周期内斜坡堤多孔介质结构内外的速度场和压力分布进行了讨论分析。计算分析表明,数值计算波面较Philip Liu的计算结果与试验结果吻合更好。说明应用SPH方法建立的二维数学模型能够较好地模拟破碎波在多孔渗水斜坡上的爬坡和越浪。     

陡坡上孤立波爬高的数值模拟

波浪爬高是海岸工程中重要的水动力学问题之一,其数值模拟方法通常是通过离散Navier-Stokes方程或Boussinesq方程实现的,其中基于光滑粒子流体动力学方法是近年发展起来的。本文应用该方法模拟相同水深下,不同波高的孤立波在45(°)陡坡上的爬高,模拟结果与理论计算结果及已有物理模型试验结果进行了对比,并模拟出孤立波激散破碎过程及粒子分布和速度场的变化过程。结果表明,对密度近似方程进行重新初始化保持了流场内的质量守恒,同时整个计算域内的压力分布更加规则,说明光滑粒子流体动力学法在波浪爬高计算中的有效性。     

波浪与起伏水平板防波堤相互作用数值模拟

利用自主研发的基于紧致插值曲线CIP(constrained interpolation profile)方法的数学模型,开展规则波与起伏水平板防波堤相互作用的数值模拟研究。模型在笛卡尔直角坐标下建立,以CIP方法为流场基本求解器,分步求解Navier-Stokes方程,利用高精度的流体体积类型的THINC/SW (tangent of hyperbola for interface capturing with slope weighting)方法重构自由液面,采用浸入边界IBM(immersed boundary method)方法处理波浪与起伏板防波堤的耦合作用问题,通过动量源项造波方法模拟波浪的产生。重点关注波浪的浅水变形和板两端涡旋脱落的非线性现象,分析不同潜深、波要素下的板周围流场分布、板的运动响应和波浪的反透射系数。结果表明:起伏水平板主要通过能量反射、板上浅水变形和板两端的涡脱落消能,能有效减小板后波高,具有作为防波堤的可行性。     

基于SPH方法的开孔沉箱比尺效应研究

由于在前壁上设置了尺寸较小的孔,开孔沉箱受流体黏性力作用显著,依照弗劳德数相似准则设计模型存在比尺效应。为揭示比尺效应,建立了模拟波浪与开孔沉箱相互作用的光滑粒子流体动力学(SPH)模型。其中流体运动由连续性方程和Navier-Stokes方程控制,固壁边界由改进的动力边界粒子施加。模型收敛性通过分析不同粒子分辨率下的波浪反射系数得到,模型精度通过比较计算与理论波浪反射系数证明。使用经过验证的SPH模型,计算并比较了不同几何比尺和开孔率下开孔沉箱附近的涡量场、箱体外侧的波面时程曲线和波浪反射系数。结果表明,随着模型几何比尺的减小,开孔沉箱受到偏大的流体黏性力,致使更多波能在湍流运动中耗散,进而减小了波浪反射系数并降低了箱体外侧的波面高度。     

低雷诺数下圆柱强迫振动特性光滑粒子流体动力学模拟

基于光滑粒子流体动力学(SPH)方法,开发了能够准确描述水流作用下圆柱强迫振动特性的数学模型。通过引入适合于无网格粒子法的开边界算法,来模拟出入流边界条件,建立了具有造流功能的SPH数值水槽。圆柱及计算域的上下边界均采用修正的动力边界条件进行模拟。借助于粒子位移矫正和压力修正算法,避免了圆柱周围流体粒子压力大幅震荡以及结构下游区域出现空腔等非物理性现象。使用典型的圆柱绕流数据,验证了所建SPH模型的计算性能,研究了固定圆柱在低雷诺数情况下的尾涡脱落模式和升阻力变化规律。明确了低雷诺数下强迫振动圆柱在频率锁定以及非锁定区间内的升力变化规律,量化了升力与外界激励频率之间的关系。     

深水斜坡式防波堤人工垫层块体稳定性及反射系数研究

船舶大型化的趋势导致港口设施向深水发展。而对于深水防波堤,在较大设计波浪作用下,护面块体和垫层块石的重量往往超出现行技术规范的要求,针对大型天然块石垫层选料和施工困难,大量的开山选石损害自然环境等问题,本文总结了国内外人工块体的研究成果,阐明了人工块体垫层代替大型块石垫层的可行性,遴选了3种类型人工块体作为垫层,进行护面块体稳定性试验,经对比分析得出扭王字护面下带槽方块垫层的新结构,并确定其重量比值关系对护面稳定性的影响;对比传统垫层块石,给出遴选人工块体垫层的K_D值、堤前反射系数经验公式及波浪爬高糙渗系数推荐值。通过物理模型试验证明了深水斜坡防波堤通过采用人工垫层块体代替大型天然垫层块石的方法是可行的,为港口工程技术规范的修订提供了参考。     

海水中痕量DMS和DMSP分析方法的研究

二甲基硫（DMS）是海洋排放的占优势地位的生源硫气体,其在大气中的氧化产物能够影响到环境酸化和世界的气候变化.因此, 测定海水中的DMS对于准确地评价其在全球硫循环所起的重要作用具有重要意义.本文中作者研究了海水中DMS的痕量分析技术.海水中的DMS首先采用气提-冷阱捕集技术进行预浓缩, 然后用带有火焰光度检测器的气相色谱（GC-FPD）进行分析.该方法的精确度在5%以内, 平均回收率为85.6% （82.8%-90.5%）, 最小检出限为0.15 ng S.β-二甲基巯基丙酸内盐（DMSP）的分析是通过将其在碱性溶液中分解成DMS来进行.作者采用此方法实测了黄海中DMS和DMSP的含量, 获得了理想的结果.     

Study on the axial exciting force characteristics of marine propellers considered the effect of the shaft and blade elasticity

Marine propellers usually operate in a spatially nonuniform wake and then the propeller exciting forces are produced. These exciting forces will cause serious hull vibration and noise radiation. So, there are many researches on the exciting forces of propellers. However, the effects of the shaft and blade elasticity are ignored in most studies. Therefore, firstly, considered the effects of the shaft and blade elasticity, a fluid–structure interaction dynamic model of the fluid-propeller-shaft system is established by coupled BEM and FEM. Then, based on this model, the characteristics of axial exciting force and theirs transmission mechanism to the hull via the shaft are studied. The research results show that first, there are two kinds of vibration modes for blade bending vibration: global mode and local mode. The elastic coupling effect between the blade and the shafting only affects the global mode, but hardly affects the local mode of the blade. Secondly, during the transmission of axial exciting force to the hull through the blades and shafting, only the global mode of the blades can amplify it, while the local mode cannot. These studies could provide a guideline for the optimal design of the propeller-shaft system to make the exciting force transmitted to the ship hull via the shafting be the smallest.     

光子晶体光纤带隙及模场的一种综合模拟分析方法

提出直接由MATLAB读入光纤端面图实现复杂结构光子晶体光纤(PCF)模型的快速建立,并综合平面波法(PWM)和频域有限差分(FDFD)法,模拟分析带隙型光子晶体光纤(PBG-PCF)的带隙和模场分布.利用PWM计算得到了在PBG-PCF中传输的光波频率及模式有效折射率范围;基于FDFD法在给定波长及模式有效折射率范围情况下,模拟得到了PBG-PCF中可能存在的模场分布及其他特性.以市售的PBG-PCF为例,验证了数值模拟的正确性,随后系统地分析了结构参数(晶格结构、"原子"占有率、背景材料折射率及"原子     

Material selection for the gas containment system of a compressed natural gas carrier fleet

In this work, the possible exploitation of fiber-reinforced composites in the context of maritime transportation of compressed natural gas (CNG) is investigated. In addition to a more conventional steel configuration, two different fiber materials, carbon and glass, are considered as construction materials for pressure vessels (PVs) to be stored on board ships, with thickness optimized by FEM analysis.The considered scenario is represented by the transportation of CNG from an offshore well to a terminal on shore. Fleets of ships carrying CNG in pressure vessels manufactured with the investigated materials are generated by means of a ship synthesis model (SSM) software and compared on the basis of technical and economical indicators.The choice of the construction material influences considerably the weight of the PVs, which represent a major item of total ship weight and reflects directly on the general transport performances in terms of resistance, seakeeping and reliability in the service. On the other hand, capital as well as operating expenditures are considerably affected by the choice. When exploring the design space, the ship synthesis model is able, at a preliminary stage of the design, to account for the various technical and economical aspects, their implications and relationships. Results are presented of computations carried out in a specific case, identified by the annual gas production and other characteristics of the well terminal and a cruising route for the ships. The comparison is carried out on the basis of the cost per transported unit of gas and of the percentage of success in the transportation process. The computations show that the choice of the PV material has a key influence on the results in terms of optimal number, dimensions and speed of the ships.     

Vibrations of a rectangular orthotropic plate with a free edge: a comparison of analytical and numerical results

Transverse vibrations of a rectangular orthotropic plate with a free edge while the others are supported according to four different combinations of boundary conditions are considered in the present paper. The presence of a concentrated mass is also taken into account. In view of the difficulty of satisfying exactly the natural boundary conditions at the free edge, the fundamental frequency of vibration is determined by means of the optimized Rayleigh-Ritz method and a pseudo-Fourier expansion which contains, in its argument, the optimization parameters γ_i. The results are in excellent agreement with those obtained by means of a finite element algorithmic procedure.     

Assessment of CPTU and static load test methods for predicting ultimate bearing capacity of pile

For the static pressure pile, the most important is to determine the standard value of ultimate bearing capacity of single pile. At present, the bearing capacity of pile foundation is usually determined by the cone penetration test (CPT) test. The empirical formula was used in practice, but the effect of excess pore water pressure generated in the penetration on the measured cone-tip resistance and side friction is neglected. In this study, based on the field test results at Fengyu playground at Yancheng Institute of Technology by CPT and CPTU methods, the bearing capacity of pile was predicted by the standardized methods, the LCPC (France method) and CPTU direct prediction methods. The prediction was also compared with the results by the static load test method. The prediction accuracy of the CPTU method was then discussed as well. The results reveal that the accuracy of the CPTU method was the highest, which was consistent with the results obtained by the static load test method. It is the best method and worthy of being applied for predicting bearing capacity of piles in engineering applications.     

温度和盐度对华贵栉孔扇贝免疫相关酶的联合效应

为研究温度和盐度对华贵栉孔扇贝免疫相关酶活力的联合效应,采用中心复合设计(CCD)和响应曲面分析法(RSM)进行试验,在实验室条件下研究了温度(19~31℃),盐度(22~38)对华贵栉孔扇贝酸性磷酸酶(ACP)、碱性磷酸酶(AKP)、谷胱甘肽还原酶(GR)三种免疫相关酶活性的联合效应,旨在考察温度与盐度对ACP,AKP和GR活性的影响及华贵栉孔扇贝免疫水平较高时的最佳温度和盐度组合,为华贵栉孔扇贝人工养殖所需条件提供理论基础。结果表明:温度的一次效应和二次效应对ACP,AKP和GR活力影响显著(P0.05);盐度的一次效应和二次效应对ACP,AKP和GR的活力影响显著(P0.05);温度和盐度的交互作用对ACP,AKP和GR活力的影响显著(P0.05)。采用响应曲面法,建立了温度和盐度对华贵栉孔扇贝ACP,AKP和GR活性影响的模型方程,该模型方程决定系数分别为0.986 6,0.981 3和0.957 6(P0.01),预测系数分别为0.934 5,0.888 0和0.738 4,表明该模型可以用于预测华贵栉孔扇贝ACP,AKP和GR活力变化。通过模型优化和验证实验,得到在温度24.30℃、盐度30.55时,ACP和AKP活力达到最小值,分别为170.43,181.74U/mg,而GR活力最大为0.34U/mg,满意度达0.95。     

A Maximum-Entropy Method for Estimating the Spectrum

Based on the maximum-entropy(ME)principle,a new power spectral estimator for random waves is derived in the form of S～(ω)=a/8H~2(2π)~(d 1)ω~-~((d 2))exp[-b(2π/ω)~n],by solving a variational problem subject to some quite general constraints.This robust method is comprehensive enough to describe the wave spectra even in extreme wave conditions and is superior to periodogram method that is not suitable to process comparatively short or intensively unsteady signals for its tremendous boundary effect and some inherent defects of FFT.Fortunately,the newly derived method for spectral estimation works fairly well,even though the sample data sets are very short and unsteady,and the reliability and efficiency of this spectral estimator have been preliminarily proved.     

Numerical experiments on the Rayleigh-Ritz method when applied to doubly connected plates in the case of free edge holes

Numerical experiments are performed on the determination of the fundamental frequency of transverse vibration of simply supported rectangular plates having rectangular holes with free edges. This constitutes a rather common technological situation since holes are present in plates or slabs due to operational conditions, namely passage of conduits or ducts, electric conductors, etc., Satisfying exactly the governing natural boundary conditions at the hole edges is practically an impossible task. This study reviews numerical experiments where the displacement function is expanded into a double Fourier series which constitutes the exact solution when the plate is simply connected. Satisfactory convergence is achieved when the plate is doubly connected.     

大气气溶胶中多环芳烃的分析

本文对国内外近年来关于大气气溶胶中多环芳烃的样品的采集、提取、分离及分析方法作了较为系统的介绍和比较 ,对各种方法的适用性作了评述