Elastic Analysis of Anisotropic Functionally Graded Rotating Disks With Non-Uniform Thicknesses北大核心CSCD

研究了任意梯度变化的变厚度各向异性转动圆盘的弹性问题.假设圆盘绕刚性轴匀速转动,其材料性能和厚度沿径向任意梯度变化.考虑圆盘在中心转轴处受位移约束,外侧自由,根据各向异性转动圆盘的平衡微分方程,得到关于径向应力的Fredholm积分方程,继而通过对Fredholm积分方程进行数值求解,得到结构的位移场和应力场.对具体梯度变化情况仅需代入相应梯度变化进行求解即可.数值算例部分,通过假设厚度、弹性模量等参数为特殊的幂函数形式,将由Fredholm积分方程求出的数值解与对应的精确解进行对比,以及针对常见的Voigt模型,将由该方法算得的数值解和ANSYS有限元计算结果进行对比,验证了该方法的准确性和精度.其次,针对Voigt模型,重点分析了厚度变化、材料性能梯度参数、各向异性度等对应力场和位移场的影响.提出了针对材料性能和厚度沿径向呈任意梯度变化的圆盘结构弹性分析方法,将为优化功能梯度圆盘的结构和材料参数、有效调整构件应力分布、提高结构安全性,提供强有力的工具;算例分析结果对功能梯度圆盘在复杂条件下的结构安全设计有重要的理论指导意义.     

Numerical Investigation of Particle Focusing Patterns in Laminar Pipe Flow With High Reynolds Numbers北大核心CSCD

The inertial focusing characteristics of particles in laminar flow pipes with high Re numbers were studied based on the “relative motion model”. In order to solve the problem of long pipes with high Re number flow, periodic boundary conditions were imposed on the inlet and outlet of the pipe. The research results show that the use of periodic boundary conditions can effectively reduce the computational, and the mechanical properties of particles in high Re flow can be calculated by using L=4D pipe. The difference from the low Re number is that as the Re number continues to increase,the lift force of the particles in the radial direction is no longer distributed as a parabola. The lift curve has a concave area between r+ =0.5 ~ 0.7, and there is a tendency for a new inertial focus point to appear in this section. By means of particles of a+ =1/17 for Re > 1 000, this new focus point position is solvable. In addition, in the analysis of the flow field, a secondary flow occurs around the particle, and its intensity gradually increases with the Re number and the closeness of the particle to the wall. The generation of the secondary flow affects the spatial distribution of the particle lift. © 2023 Editorial Office of Applied Mathematics and Mechanics. All rights reserved.     

Mesoscopic Numerical Study on Flow Boiling Heat Transfer Performance in Channels With Multiple Rectangular Heaters北大核心CSCD

The flow boiling phenomenon in a channel with multiple rectangular heaters under a constant wall temperature was numerically studied with the lattice Boltzmann method. The effects of spacings between heaters, heater lengths and heater surface wettabilities on the bubble morphology, the bubble area and the heat flux on the heater surface, were studied. The results show that, the bubble growth rate increases with the spacing between heaters. The larger the bubble area is, the earlier the nucleated bubbles will leave the heater surface. The corresponding boiling heat transfer performance increases by 12% with the spacing between heaters growing from 250 lattices to 1 000 lattices. On the other hand, the longer the heater length is, the earlier the bubble will nucleate and leave the heater surface, and the better the boiling heat transfer performance will be. The boiling heat transfer performance increases by 13% with the heater length rising from 16 lattices to 22 lattices. In addition, the bubble nucleates later on the hydrophilic surface than on the hydrophobic surface. Compared with the hydrophilic surface, the hydrophobic surface retains residual bubbles after the leaving of bubbles from the heater. The average heat flux and the bubble area of the hydrophilic surface are less than those of the hydrophobic surface. With the contact angle changing from 77° to 120°, the heat transfer performance increases by 26%. Finally, the orthogonal test results indicate that, the wettability of the heat exchanger surface has the greatest influence on the flow boiling heat transfer performance, while the heater length has the least influence. © 2022 Editorial Office of Applied Mathematics and Mechanics. All rights reserved.     

Simulation Study on Dam Break Flow Based on the B-Spline Material Point Method北大核心CSCD

The dam break flow poses a common free surface flow problem in hydraulic engineering, and its ac⁃ curate simulation is of great engineering significance. The B⁃spline material point method (BSMPM), as an im⁃ proved algorithm of the material point method (MPM), has optimized accuracy and convergence in material point calculations and unique algorithmic advantages in free surface flow problems. Based on the BSMPM, a weakly compressible BSMPM (WC⁃BSMPM) was developed through introduction of an artificial equation of state. The simulation of the dam break flow problem was carried out, with the effects of the order of the B⁃ spline interpolation basis function on the simulation results analyzed. The results show that, the simulated fluid wavefront position, the wavefront velocity and the elevation variation at a given position are basically consistent with the existing experimental results. As the order of the basis function increases, the computation time will lengthen for about 1.5 times. However, the computation times of the BSMPM of different orders will uniformly increase approximately linearly with the background grid size. The validity of the WC⁃BSMPM simulation of the dam break flow problem was verified. The research provides a new idea and method for the simulation of dam break flow problems. © 2023 Editorial Office of Applied Mathematics and Mechanics. All rights reserved.     

Lattice Boltzmann Study on the Motion of Dual Droplets in Microchannels With Contact Angle Hysteresis北大核心CSCD

接触角滞后表现为流体在非理想固体表面上运动时前进接触角和后退接触角不同,是两相流体在润湿表面上流动的重要现象.该文采用改进的伪势格子Boltzmann(LB)多组分模型,并与几何润湿边界条件相结合,研究了两个液滴在具有接触角滞后性微通道表面上的运动行为,主要研究了通道内特征数、通道表面性质以及液滴初始参数的影响.研究结果表明：毛细数的增大有助于液滴的移动,然而并不利于液滴的排出,且毛细数的增加对上游液滴的影响大于其对下游液滴的影响;另一方面,接触角滞后性窗口越大,液滴运动和形变更迟缓,但形变程度更明显,两液滴更早地发生合并,但更晚地排出管道;液滴间距的增加使液滴的运动行为在不同阶段表现为不同的模式,但都导致通道中残留小液滴,使得液滴排出通道的时间增加.研究结果还表明：上游液滴和下游液滴的相对尺寸差距越大,越不利于液滴排出管道.     

Dynamics Analysis of Cannibalistic Model With Density Dependence in Mature Stage北大核心CSCD

在考虑成熟阶段具有密度制约的基础上,建立了一类具有卵-成熟阶段的同类相食模型.该文从两个方面讨论了模型的动力学性态：当种群不存在同类相食时,构造Lyapunov函数证明平衡点的全局渐近稳定性;当种群存在同类相食时,利用中心流形定理证明同类相食使模型产生鞍结点分支,通过构造Dulac函数说明在二维自治系统中不存在极限环,得到了平衡点的全局稳定性.最后,利用数值模拟验证了所得相应结果的正确性.     

Simulation of Electroosmotic and Pressure-Driven Mixed Flow of Viscoelastic Fluids in Converging-Diverging Tubes北大核心CSCD

The electroosmotic and pressure⁃driven mixed flow was widely used in various biochemical microflu⁃ idic fields, where the elastic instability of the viscoelastic fluid cannot be ignored. A viscoelastic fluid was used to numerically simulate electroosmotic and pressure⁃driven mixed flow in a 10 ∶ 1 ∶ 10 microchannel converging⁃ diverging tube. The effects of different pressures and different polymer concentrations on the fluid flow were studied, and the superposition principle for the velocity distributions of Newtonian fluids and viscoelastic fluids in converging⁃diverging tubes was analyzed. The results show that, the reverse pressure brings the viscoelastic fluid into higher instability, which makes the inlet vortex larger by 25 μm for every 1 Pa pressure increase. The positive pressure makes the eddy current smaller. For a relatively small reverse pressure, the inlet vortex in⁃ creases with the polymer concentration and tends to be stable gradually. For a relatively large reverse pressure, the vortex size first increases and then decreases with the polymer concentration. © 2023 Editorial Office of Applied Mathematics and Mechanics. All rights reserved.     

UGV Path Programming Based on the DQN With Noise in the Output Layer北大核心CSCD

在DQN算法的框架下,研究了无人车路径规划问题.为提高探索效率,将处理连续状态的DQN算法加以变化地应用到离散状态,同时为平衡探索与利用,选择仅在DQN网络输出层添加噪声,并设计了渐进式奖励函数,最后在Gazebo仿真环境中进行实验.仿真结果表明:①该策略能快速规划出从初始点到目标点的无碰撞路线,与Q-learning算法、DQN算法和noisynet_DQN算法相比,该文提出的算法收敛速度更快;②该策略关于初始点、目标点、障碍物具有泛化能力,验证了其有效性与鲁棒性.     

Numerical Simulation of Hydraulic Fractures Intersecting Natural Fractures in Shale With Plastic Deformation北大核心CSCD

The plastic deformation and numerous natural joints of shale pose a great challenge for the predic⁃ tion of the hydraulic fracture geometry extension. Based on the finite element method, a fully coupled numeri⁃ cal model for elastoplastic hydraulic fractures was established with natural fractures and bedding planes consid⁃ ered. The numerical model was validated with the KGD analytical solution and Blanton’s curve. The numerical results show that, compared with the numerical model solution of linear elasticity, the hydraulic fractures are prone to enter the natural weak interface due to the rock plastic deformation. The rock plastic deformation area mainly lies in the reservoir layer during the fracture propagation. In the case of rock ductile damage, the hy⁃ draulic fracture is more likely to penetrate the bedding plane. Hydraulic fractures can directly penetrate natural fractures and bedding planes at high injection rates due to large driving forces. The study provides new insights in terms of hydraulic fracture extension in elastoplastic formations. © 2023 Editorial Office of Applied Mathematics and Mechanics. All rights reserved.     

Fundamental Solution of Rotating Generalized Stokes Problem in R~3

In this paper,the fundamental solution of rotating generalized Stokes problem in R 3 is established.To obtain it,some fundamental solutions of other problems also are established,such as generalized Laplace problem,generalized Stokes problem and rotating Stokes problem.     

旋转功能梯度圆柱壳的固有频率计算

Love的一阶理论被用来计算旋转功能梯度圆柱壳的固有频率, 为了验证该方法的有效性,计算了简支不旋转各向同性和功能梯度圆柱壳的固有频率,计算结果与相关文献中的结果具有较好的一致性．通过几个数值算例,研究了幂指数、x和θ方向的波数、厚度-半径比对简支旋转功能梯度圆柱壳固有频率的影响．结果表明:后向波的固有频率随着转速的增加而增加,前向波的固有频率随着转速的增加而减小,前向波和后向波的固有频率随着厚度-半径比增加而增加．     

A Hybrid Immersed Interface Method for Driven Stokes Flow in an Elastic Tube

We present a hybrid numerical method for simulating fluid flow through a compliant, closed tube, driven by an internal source and sink. Fluid is assumed to be highly viscous with its motion described by Stokes flow. Model geometry is assumed to be axisymmetric, and the governing equations are implemented in axisymmetric cylindrical coordinates, which capture 3D flow dynamics with only 2D computations. We solve the model equations using a hybrid approach: we decompose the pressure and velocity fields into parts due to the surface forcings and due to the source and sink, with each part handled separately by means of an appropriate method. Because the singularly-supported surface forcings yield an unsmooth solution, that part of the solution is computed using the immersed interface method. Jump conditions are derived for the axisymmetric cylindrical coordinates. The velocity due to the source and sink is calculated along the tubular surface using boundary integrals. Numerical results are presented that indicate second-order accuracy of the method.     

粘滑混合边界条件下平面边界前的Stokes流动

对具有粘滑混合边界条件的平面边界,建立一个Stokes流动的一般性定理,利用双调和函数A与调和函数B,表示了3维Stokes流动的速度场和压力场.关于无滑动平面边界前Stokes流动的早期定理,成为该一般性定理的一个特例.进一步地,从一般性定理导出了一个推论,根据该Stokes流函数,给出了粘滑边界条件时刚性平面轴对称Stokes流动问题的解,得到了流体作用在边界上的牵引力和扭矩公式.给出了一个说明性的例子.     

计及Hall和离子滑移电流的旋转流体在与时间相关初值条件下的MHDCouette流动

考虑Hall和离子滑移电流的影响,在旋转系统中研究导电流体非稳定的MHD Couette流动．在小数值磁场Reynolds数假定下,推导出基本的控制方程,使用著名的Laplace变换技术,数值地求解该基本方程．分两种情况:磁场相对于流体或者移动平板固定时,得到速度和表面摩擦力统一的闭式表达式．用图形讨论了问题的不同参数,对速度和表面摩擦力的影响．所得结果显示,主流速度u和次生速度v随着Hall电流而增大．离子滑移电流的增大,也会导致主流速度u的增大,但会使次生速度v减小．还给出了旋转、Hall和离子滑移参数的综合影响,确定了次生运动对流体流动的贡献．     

Inversion transformation of Stokes flows bounded by parallel planes

The positive inversion transformation applied to a two-dimensional Stokes flow around bodies leads alike to a Stokes flow. This fact can be exploited to find new two-dimensional Stokes flow solutions around inverse bodies. Some features of this method, such as the relations between the reference and inverse fluid velocity fields, are presented followed by an application to examples of cellular flow between two parallel plates induced by rotating or translating cylinder. Thus hydrodynamic characteristics of flow around circular bodies obtained by inversion of the plates are straightforward deduced. Typical fluid flow patterns around two circular cylinders in contact placed in the centre of a rotating or a translating circular cylinder are thus illustrated.     

三维Stokes绕流问题的无限元逼近(英文)

设计了一种最少自由度的无限元方法来实现三维Stokes绕流问题的求解.通过验证强制性和inf-sup条件,我们证明了相应的离散混合变分问题解的适定性,并在加权Sobolev空间中得到了误差的先验估计.数值实验结果验证了解的收敛性.     

柱状固粒两相边界层流场稳定性研究

从运动方程和本构方程出发,推导得到了含柱状粒子两相流场的修正Orr-Sommerfeld方程,然后在边界层流场中,采用数值计算方法,得到了含柱状粒子流场的稳定性中性曲线,给出了流场失稳的临界雷诺数.结果表明在所述情况下,柱状粒子对流场起着抑制失稳的作用,而且抑制的程度随着柱状粒子体积分数和长径比的增加而提高.     

平面边界前Stokes流动基本的奇异性

利用双调和函数A和调和函数B，给出了三维Stokes流动速度场和压力场的描述．由此建立了计算区域边界为固定无滑移平面边界Stokes流动基本奇异性的一般定理．刚性平面前轴对称Stokes流动的Collins定理成为本定理的特例．给出的几个例证说明了方法的有效性．     

微极流体蠕动泵经由滑移边界管道输送的Stokes流动

计及管道边界条件滑移的影响,研究微极流体蠕动泵,经由圆柱形管道输运的Stokes流动.壁面运动的控制方程为正弦波方程.使用润滑理论,得到了轴向速度、微转动向量、流函数、压力梯度、摩擦力和机械效率的解析数值解.用图形表示出构成参数,如像耦合参数、微极参数和表征蠕流泵特性的滑移参数、摩擦力和俘获现象的影响.数值计算表明,当耦合参数较大时,需要蠕动泵的压力更大,而微极参数和滑移参数正相反.俘获团块的大小随耦合参数和微极参数的减小而缩小,而随滑移参数的增大而缩小.