Flow Field Characteristics of Multi-Trophic Artificial Reef Based on Computation Fluid Dynamics

On the basis of computational fluid dynamics,the flow field characteristics of multi-trophic artificial reefs,including the flow field distribution features of a single reef under three different velocities and the effect of spacing between reefs on flow scale and the flow state,were analyzed.Results indicate upwelling,slow flow,and eddy around a single reef.Maximum velocity,height,and volume of upwelling in front of a single reef were positively correlated with inflow velocity.The length and volume of slow flow increased with the increase in inflow velocity.Eddies were present both inside and backward,and vorticity was positively correlated with inflow velocity.Space between reefs had a minor influence on the maximum velocity and height of upwelling.With the increase in space from 0.5 L to 1.5 L(L is the reef lehgth),the length of slow flow in the front and back of the combined reefs increased slightly.When the space was 2.0 L,the length of the slow flow decreased.In four different spaces,eddies were present inside and at the back of each reef.The maximum vorticity was negatively correlated with space from 0.5 L to 1.5 L,but under 2.0 L space,the maximum vorticity was close to the vorticity of a single reef under the same inflow velocity.

Structural parameter optimization for novel internal-loop iron–carbon micro-electrolysis reactors using computational fluid dynamics

It is generally recognized that internal-loop reactors are well-developed mass and heat-transfer multiphase flow reactors. However, the internal flow field in the internal-loop reactor is influenced by the structure parameter of the reactor, which has a great effect on the reaction efficiency. In this study, the computational fluid dynamics simulation method was used to determine the influence of reactor structure on flow field, and a volume-offluid model was employed to simulate the gas–liquid, two-phase flow of the internal-loop micro-electrolysis reactor. Hydrodynamic factors were optimized when the height-to-diameter ratio was 4:1, diameter ratio was9:1, draft-tube axial height was 90 mm. Three-dimensional simulations for the water distributor were carried out, and the results suggested that the optimal conditions are as follows: the number of water distribution pipes was four, and an inhomogeneous water distribution was used. According to the results of the simulation,the suitable structure can be used to achieve good fluid mechanical properties, such as the good liquid circulation velocity and gas holdup, which provides a good theoretical foundation for the application of the reactor.

Application of Computational Fluid Dynamics and Fluid Structure Interaction Techniques for Calculating the 3D Transient Flow of Journal Bearings Coupled with Rotor Systems

Journal bearings are important parts to keep the high dynamic performance of rotor machinery. Some methods have already been proposed to analysis the flow field of journal bearings, and in most of these methods simplified physical model and classic Reynolds equation are always applied. While the application of the general computational fluid dynamics （CFD）-fluid structure interaction （FSI） techniques is more beneficial for analysis of the fluid field in a journal bearing when more detailed solutions are needed. This paper deals with the quasi-coupling calculation of transient fluid dynamics of oil film in journal bearings and rotor dynamics with CFD-FSI techniques. The fluid dynamics of oil film is calculated by applying the so-called ＂dynamic mesh＂ technique. A new mesh movement approacb is presented while the dynamic mesh models provided by FLUENT are not suitable for the transient oil flow in journal bearings. The proposed mesh movement approach is based on the structured mesh. When the joumal moves, the movement distance of every grid in the flow field of bearing can be calculated, and then the update of the volume mesh can be handled automatically by user defined function （UDF）. The journal displacement at each time step is obtained by solving the moving equations of the rotor-bearing system under the known oil film force condition. A case study is carried out to calculate the locus of the journal center and pressure distribution of the journal in order to prove the feasibility of this method. The calculating results indicate that the proposed method can predict the transient flow field of a journal bearing in a rotor-bearing system where more realistic models are involved. The presented calculation method provides a basis for studying the nonlinear dynamic behavior of a general rotor-bearing system.

Three-dimensional Computational Fluid Dynamics Modeling of Two-phase Flow in a Structured Packing Column

Characterizing the complex two-phase hydrodynamics in structured packed columns requires a powerful modeling tool.The traditional two-dimensional model exhibits limitations when one attempts to model the detailed two-phase flow inside the columns.The present paper presents a three-dimensional computational fluid dynamics(CFD) model to simulate the two-phase flow in a representative unit of the column.The unit consists of an entire corrugation channel and describes well the real liquid flow conditions.The detailed unsteady two-phase 3D CFD calculations on column packed with Flexipak 1Y were implemented within the volume of fluid(VOF) mathematical framework.The CFD model was validated by comparing the calculated thickness of liquid film with the available experimental data.Special attention was given to quantitative analysis of the effects of gravity on the hydrodynamics.Fluctuations in the liquid mass flow rate and the calculated pressure drop loss were found to be qualitatively in agreement with the experimental observations.

Computational dynamics of soft machines

Soft machine refers to a kind of mechanical system made of soft materials to complete sophisticated missions, such as handling a fragile object and crawling along a narrow tunnel corner, under low cost control and actuation. Hence, soft machines have raised great challenges to computational dynamics. In this review article, recent studies of the authors on the dynamic modeling, numerical simulation, and experimental validation of soft machines are summarized in the framework of multibody system dynamics. The dynamic modeling approaches are presented first for the geometric nonlinearities of coupled overall motions and large deformations of a soft component, the physical nonlinearities of a soft component made of hyperelastic or elastoplastic materials, and the frictional contacts/impacts of soft components, respectively. Then the computation approach is outlined for the dynamic simulation of soft machines governed by a set of differential-algebraic equations of very high dimensions, with an emphasis on the efficient computations of the nonlinear elastic force vector of finite elements. The validations of the proposed approaches are given via three case studies, including the locomotion of a soft quadrupedal robot, the spinning deployment of a solar sail of a spacecraft, and the deployment of a mesh reflector of a satellite antenna, as well as the corresponding experimental studies. Finally, some remarks are made for future studies.

调谐液柱阻尼器-结构系统风致振动响应的CFD/CSD耦合分析方法

针对调谐液柱阻尼器(tuned liquid column damper, TLCD)难以建立其精确的非线性理论分析模型,且其力学性能试验成本高和耗时长等问题,首先采用计算流体动力学(CFD)数值模拟方法,对TLCD系统的力学性能和动力特征进行仿真模拟,在此基础上进一步提出了基于计算流体动力学/计算结构动力学(CFD/CSD)耦合分析方法,求解带TLCD系统的高层建筑结构的风致动力响应。通过开展某一TLCD系统在特定底部激励下的力学性能和动力特性试验,得到其内液体晃荡的自由液面波高和晃动力时程,验证了CFD数值模拟方法可以准确地分析TLCD水箱内液体的非线性晃动特征。随后对风工程领域广泛采用的76层建筑结构振动控制Benchmark模型,假设其顶部设置TLCD系统时主体结构在三种风速重现期(10、50和100年)风速对应的横风向动力风荷载激励下的风致控制效率,采用提出的CFD/CSD耦合分析方法,进行了数值仿真模拟分析。耦合分析结果表明,TLCD系统对Benchmark模型的风致加速度、速度和位移响应均有一定的控制效果,对加速度响应的控制效果要优于对位移响应的控制效果。该研究方法可为复杂TLCD系统对高层建筑的风振控制分析提供有效的参考。

Structural Parameter Analyses on Rotor Airloads with New Type Blade-Tip Based on CFD/CSD Coupling Method

For accurate aeroelastic analysis,the unsteady rotor flowfield is solved by computational fluid dynamics(CFD)module based on RANS/Euler equations and moving-embedded grid system,while computational structural dynamics(CSD)module is introduced to handle blade flexibility.In CFD module,dual time-stepping algorithm is employed in temporal discretization,Jameson two-order central difference(JST)scheme is adopted in spatial discretization and B-L turbulent model is used to illustrate the viscous effect.The CSD module is developed based on Hamilton′s variational principles and moderate deflection beam theory.Grid deformation is implemented using algebraic method through coordinate transformations to achieve deflections with high quality and efficiency.A CFD/CSD loose coupling strategy is developed to transfer information between rotor flowfield and blade structure.The CFD and the CSD modules are verified seperately.Then the CFD/CSD loose coupling is adopted in airloads prediction of UH-60A rotor under high speed forward flight condition.The calculated results agree well with test data.Finally,effects of torsional stiffness properties on airloads of rotors with different tip swept angles(from 10° forward to 30° backward)are investigated.The results are evaluated through pressure distribution and airloads variation,and some meaningful conclusions are drawn the moderated shock wave strength and pressure gradient caused by varied tip swept angle and structural properties.

Numerical Simulation of ATPS Parachute Transient Dynamics Using Fluid-Structure Interaction Method

In order to simulate and analyze the dynamic characteristics of the parachute from advanced tactical parachute system(ATPS),a nonlinear finite element algorithm and a preconditioning finite volume method are employed and developed to construct three dimensional parachute fluid-structure interaction(FSI)model.Parachute fabric material is represented by membrane-cable elements,and geometrical nonlinear algorithm is employed with wrinkling technique embedded to simulate the large deformations of parachute structure by applying the NewtonRaphson iteration method.On the other hand,the time-dependent flow surrounding parachute canopy is simulated using preconditioned lower-upper symmetric Gauss-Seidel(LU-SGS)method.The pseudo solid dynamic mesh algorithm is employed to update the flow-field mesh based on the complex and arbitrary motion of parachute canopy.Due to the large amount of computation during the FSI simulation,massage passing interface(MPI)parallel computation technique is used for all those three modules to improve the performance of the FSI code.The FSI method is tested to simulate one kind of ATPS parachutes to predict the parachute configuration and anticipate the parachute descent speeds.The comparison of results between the proposed method and those in literatures demonstrates the method to be a useful tool for parachute designers.

An efficient technique for recovering responses of parameterized structural dynamic problems

In this article,an effective technique is developed to efficiently obtain the output responses of parameterized structural dynamic problems.This technique is based on the conception of reduced basis method and the usage of linear interpolation principle.The original problem is projected onto the reduced basis space by linear interpolation projection,and subsequently an associated interpolation matrix is generated.To ensure the largest nonsingularity,the interpolation matrix needs to go through a timenode choosing process,which is developed by applying the angle of vector spaces.As a part of this technique,error estimation is recommended for achieving the computational error bound.To ensure the successful performance of this technique,the offline-online computational procedures are conducted in practical engineering.Two numerical examples demonstrate the accuracy and efficiency of the presented method.

Computational Fluid Dynamics Approach for Predicting Pipeline Response to Various Blast Scenarios: A Numerical Modeling Study

Recent industrial explosions globally have intensified the focus in mechanical engineering on designing infras-tructure systems and networks capable of withstanding blast loading.Initially centered on high-profile facilities such as embassies and petrochemical plants,this concern now extends to a wider array of infrastructures and facilities.Engineers and scholars increasingly prioritize structural safety against explosions,particularly to prevent disproportionate collapse and damage to nearby structures.Urbanization has further amplified the reliance on oil and gas pipelines,making them vital for urban life and prime targets for terrorist activities.Consequently,there is a growing imperative for computational engineering solutions to tackle blast loading on pipelines and mitigate associated risks to avert disasters.In this study,an empty pipe model was successfully validated under contact blast conditions using Abaqus software,a powerful tool in mechanical engineering for simulating blast effects on buried pipelines.Employing a Eulerian-Lagrangian computational fluid dynamics approach,the investigation extended to above-surface and below-surface blasts at standoff distances of 25 and 50 mm.Material descriptions in the numerical model relied on Abaqus’default mechanical models.Comparative analysis revealed varying pipe performance,with deformation decreasing as explosion-to-pipe distance increased.The explosion’s location relative to the pipe surface notably influenced deformation levels,a key finding highlighted in the study.Moreover,quantitative findings indicated varying ratios of plastic dissipation energy(PDE)for different blast scenarios compared to the contact blast(P0).Specifically,P1(25 mm subsurface blast)and P2(50 mm subsurface blast)showed approximately 24.07%and 14.77%of P0’s PDE,respectively,while P3(25 mm above-surface blast)and P4(50 mm above-surface blast)exhibited lower PDE values,accounting for about 18.08%and 9.67%of P0’s PDE,respectively.Utilising energy-absorbing materials such as thin coatings of ultra-high-strength concrete,metallic foams,carbon fiber-reinforced polymer wraps,and others on the pipeline to effectively mitigate blast damage is recommended.This research contributes to the advancement of mechanical engineering by providing insights and solutions crucial for enhancing the resilience and safety of underground pipelines in the face of blast events.

Influence of Surrounding Structures upon the Aerodynamic and Acoustic Performance of the Outdoor Unit of a Split Air-Conditioner

DC-inverter split air-conditioner is widely used in Chinese homes as a result of its high-efficiency and energy-saving. Recently, the researches on its outdoor unit have focused on the influence of surrounding structures upon the aerodynamic and acoustic performance, however they are only limited to the influence of a few parameters on the performance, and practical design of the unit requires more detailed parametric analysis. Three-dimensional computational fluid dynamics（CFD） and computational aerodynamic acoustics（CAA） simulation based on FLUENT solver is used to study the influence of surrounding structures upon the aforementioned properties of the unit. The flow rate and sound pressure level are predicted for different rotating speed, and agree well with the experimental results. The parametric influence of three main surrounding structures（i.e. the heat sink, the bell-mouth type shroud and the outlet grille） upon the aerodynamic performance of the unit is analyzed thoroughly. The results demonstrate that the tip vortex plays a major role in the flow fields near the blade tip and has a great effect on the flow field of the unit. The inlet ring＇s size and throat＇s depth of the bell-mouth type shroud, and the through-flow area and configuration of upwind and downwind sections of the outlet grille are the most important factors that affect the aerodynamic performance of the unit. Furthermore, two improved schemes against the existing prototype of the unit are developed, which both can significantly increase the flow rate more than 6 %（i.e. 100 m3·h～（-1）） at given rotating speeds. The inevitable increase of flow noise level when flow rate is increased and the advantage of keeping a lower rotating speed are also discussed. The presented work could be a useful guideline in designing the aerodynamic and acoustic performance of the split air-conditioner in engineering practice.

Structure optimization of the organ-pipe cavitating nozzle and its erosion ability test on hydrate-bearing sediments

Cavitating jet is a promising drilling rate improvement technology in both the marine natural gas hydrate (NGH) fluidization exploitation method and the integrated radial jet drilling and completion method. In present study, we aim to improve the efficiency of jet erosion and extracting NGH. With a computational fluid dynamics (CFD) method, the pressure, velocity and cavitation field characteristics of organ-pipe cavitating jet (OPCJ) are analysed. The divergent angle, throat length, and divergent length of OPCJ nozzle are preferred to obtain stronger jet cavitation erosion effect. Laboratory experiments of gas hydrate-bearing sediments (GHBS) erosion by OPCJ and conical jet (CJ) are conducted to compare and validate the jet erosion performance. The impinging models of OPCJ and CJ are constructed to study the impact characteristics. Results show that the preferred values of divergent angle, throat length, and divergent length are 15°, 1d, and 3d, respectively, in present simulation conditions. For GHBS, the OPCJ possesses the advantages of high efficiency and low energy consumption. Moreover, the OPCJ has higher penetration efficiency, while showing equivalent penetration ability compared to CJ. During the impinging process, the OPCJ can induce stronger impact pressure and turbulence effect, and also shows stronger chambering effect and bottom cleaning ability compared to CJ. This study presents the erosion performance of OPCJ and CJ on GHBS, and provides preliminary insights on the potential field applications in NGH exploitation.

Crystal MD:Molecular Dynamic Simulation Software for Metal with BCC Structure

Material irradiation effect plays an important role in material science.However,it is lack of high-throughput irradiation facility and process of evolution and development,which lead to lack of basic scientific theory about atomic scale materials design and development guidance.High-performance computing for simulation makes deeply understanding of micro-level-material possible.In this paper,a new data structure is proposed for the parallel simulation of metal materials evolution with crystal structure under irradiation defects.Compared with LAMMPS and IMD,which are two popular molecular dynamic simulation versions,our method takes much less memory on multi-core clusters.

Adaptive Lagrange finite element methods for high precision vibrations and piezoelectric acoustic wave computations in SMT structures and plates with nano interfaces

This paper discusses the validity of (adaptive) Lagrange generalized plain finite element method (FEM) and plate element method for accurate analysis of acoustic waves in multi-layered piezoelectric structures with tiny interfaces between metal electrodes and surface mounted piezoelectric substrates. We have come to conclusion that the quantitative relationships between the acoustic and electric fields in a piezoelectric structure can be accurately determined through the proposed finite element methods. The higher-order Lagrange FEM proposed for dynamic piezoelectric computation is proved to be very accurate (prescribed relative error 0.02% - 0.04% ) and a great improvement in convergence accuracy over the higher order Mindlin plate element method for piezoelectric structural analysis due to the assumptions and corrections in the plate theories.The converged lagrange finite element methods are compared with the plate element methods and the computedresults are in good agreement with available exact and experimental data. The adaptive Lagrange finite elementmethods and a new FEA computer program developed for macro- and micro-scale analyses are reviewed, and recently extended with great potential to high-precision nano-scale analysis in this paper and the similarities between piezoelectric and seismic wave propagations in layered structures and plates are stressed.

穿孔对穿孔型Kenics静态混合器混合性能的影响

为了提升传统Kenics静态混合器的混合效果,在传统螺旋插入带结构中引入穿孔结构,设计出穿孔型Kenics静态混合器。采用浓度变异系数评估穿孔内径(d)、穿孔间距(S)、混合单元长径比(Ar)对穿孔型静态混合器混合性能的影响。模拟结果表明,设置的穿孔结构会显著改变流体的流动状态,在穿孔处区域产生加速流动,进而提升混合器的混合性能;随着d的增加,流体之间的对流和交换作用逐渐增强,混合效果逐渐增强,但随着d的持续增加,导致混合器内部螺旋通道被破坏,混合效果反而减弱;仅改变S对混合效果没有明显影响;对于相同的混合器长度,Ar的减小导致混合元件数量增加,液体之间的扭转、剪切效果增强,混合器混合效果更好。

喷水推进器进水流道格栅流固耦合特性研究

为了研究格栅部件对喷水推进船舶效率的影响,本文拟采用流固耦合的数值计算方法系统地对格栅表面受力及其强度特性进行计算,结合有限元方法和流固耦合理论,以STARCCM+和ABAQUS软件为工具,根据目前已知的格栅剖面设计方法,对目前格栅的剖面形式进行改进,形成新型对称剖面形式的格栅。基于上述计算流体力学与流固耦合计算方法对其进行周围流场、结构静力强度、结构动力强度进行计算与分析。结果表明:对称形式的格栅剖面对于降低格栅的变形与应力水平有一定的帮助,采用新型对称形式的剖面结构使格栅总体变形模式发生了变化,但是最大变形位置与最大应力出现在格栅的位置未发生明显变化,新型对称剖面形式格栅具有比较好的静强度与动强度特性,各栅条均匀承受流体载荷。

基于CEL方法的带结构物溃坝流数值模拟

随着全球气温的不断升高,极端天气现象增多,溃坝现象成为水工结构方向需要重点防范的自然灾害。溃坝流的发展会对下游结构物产生巨大的破坏作用,因此深入探究溃坝流现象的水动力学特性以及溃坝流下游结构物所受到的冲击力影响势在必行。本文针对带有结构物的溃坝流动问题现象,采用CEL方法(Coupled Eulerian-Lagrangian analysis,耦合欧拉-拉格朗日方法),对于溃坝流的发展、演化进行了数值模拟。研究结果表明,模拟结果中溃坝流自由表面及流体尖端时程变化与文献试验结果拟合较好,CEL方法可应用在溃坝流现象的模拟,并且具有较好的准确性和精度。结构物材料性能、高度及与结构物与溃坝流间的位置,均会影响溃坝流的运动状态以及结构物的受力变形状态,在不考虑结构物破坏的情况下,提高结构物弹性模量,增加结构物高度,以及减小结构物与溃坝流间的距离均可提高结构物对于溃坝流的阻挡作用。

基于GA-BPNN代理模型的混合式黏滞阻尼器构造优化

采取CFD(computational fluid dynamics)数值模拟和基于代理模型的优化方法,研究混合式黏滞阻尼器在多参数设计空间内的构造参数优化问题。首先,分析黏滞阻尼器参数对力学性能的影响规律,基于正交试验理论选取间隙宽度、缸体内径、活塞厚度、阻尼孔直径为关键参数,利用优化拉丁超立方设计和CFD数值模拟获得初始样本点耗能性能指标;然后,使用遗传算法(genetic algorithm,GA)和反向传播神经网络(backpropagation neural network,BPNN)的方法构建代理模型,其能够准确反映构造参数耗能性能指标和构造参数之间的繁复关系。进一步地,循环更新GA-BPNN代理模型并结合非支配排序遗传算法(non-dominated sorting genetic algorithm-II,NSGA-II),以最佳耗能能力为目标,实现在全设计空间内的寻优,确定在典型阻尼介质黏度下双出杆混合式黏滞阻尼器的最优构造参数。最后,探讨了构造参数对黏滞阻尼器流动性能的影响机理。研究结果表明,优化双出杆混合式黏滞阻尼器构造参数在选用基于GA-BPNN代理模型的优化方法时,可有效提高计算效率且具有较高预测精度。从性能指标提升角度,全黏度下最优构造参数阻尼孔直径、活塞厚度、缸体内径和间隙宽度分别为1.31 mm、120.78 mm、199.87 mm、0.5 mm,相较于基准构造,最大阻尼力和能量耗材比分别提升33.8%和64.8%。该方法可为黏滞阻尼器构造参数优化及相关研究提供参考。

高旋弹二维修正引信双旋结构气动特性及气动力表征

针对配二维弹道修正引信高旋弹具有弹体气动参数非对称、纵向和横向修正紧密耦合等问题,为了准确表征高旋弹气动参数、明确修正弹丸气动特性和产生机理,提出了基于CFD(computational fluid dynamics)仿真的双旋结构气动力计算分析方法.在构建双滚转域流场仿真模型的基础上,对比了二维修正引信不同控制状态下的弹丸受力情况;明确了高旋弹固有气动力和二维修正引信所引起的气动力;建立并推导了攻角与滚转角耦合情况下的舵片受力模型.研究表明:针对二维修正组件,需要考虑合攻角与舵滚转角的相对位置关系以计算诱导阻力;受迎背风和舵片绕流影响,舵片受力模型和弹体的横、纵向气动力均随攻角、滚转角及马赫数变化.