Interface‐fitted moving mesh method for axisymmetric two‐phase flow in microchannels

A boundary‐fitted moving mesh scheme is presented for the simulation of two‐phase flow in two‐dimensional and axisymmetric geometries. The incompressible Navier‐Stokes equations are solved using the finite element method, and the mini element is used to satisfy the inf‐sup condition. The interface between the phases is represented explicitly by an interface adapted mesh, thus allowing a sharp transition of the fluid properties. Surface tension is modelled as a volume force and is discretized in a consistent manner, thus allowing to obtain exact equilibrium (up to rounding errors) with the pressure gradient. This is demonstrated for a spherical droplet moving in a constant flow field. The curvature of the interface, required for the surface tension term, is efficiently computed with simple but very accurate geometric formulas. An adaptive moving mesh technique, where smoothing mesh velocities and remeshing are used to preserve the mesh quality, is developed and presented. Mesh refinement strategies, allowing tailoring of the refinement of the computational mesh, are also discussed. Accuracy and robustness of the present method are demonstrated on several validation test cases. The method is developed with the prospect of being applied to microfluidic flows and the simulation of microchannel evaporators used for electronics cooling. Therefore, the simulation results for the flow of a bubble in a microchannel are presented and compared to experimental data.     

无人机视角下的林木爆燃火特性研究

森林火灾“爆燃”现象的特征是突然发生的高强度、高蔓延速度的燃烧。目前为止，关于“爆燃火”的原因还没有达成共识。以无人机视角下对林木爆燃火特性研究，以四川木里特大森林火灾为研究对象，通过分析凉山州某森林扑火部队3月31日木里森林火灾当天KWT（科卫泰）无人机航拍火场画面，结合无人机实时影像及实地调研数据，分析了峡谷地形林火蔓延时空特征，探讨了峡谷中风向风速变化时空分布规律，研究了地形变化条件下，不同海拔高度风速特征，建立了无人机倾角测量风速模型（其中为风速m·s-1，为无人机倾角°）。结果表明，高山峻岭特殊地形环境下每天4:00—12:00时间段为静风期，为峡谷林火扑救最佳时期；午后15:00—17:00和晚上20:00—22:00为山谷地形风速活跃期；仿真软件数据显示山顶、谷底与山腰不同海拔位置的风速风向不统一，谷底会产生乱流现象，且风速与海拔不存在正相关关系，小气候在复杂地形中占主导影响地位；在谷口至山谷深处的中间位置会出现气流速度的波峰状态，并易形成乱流，为爆燃火发生提供了客观必要条件。该研究可为复杂地形环境下，森林草原火灾扑救安全提供数据和技术支撑。     

Invariant discrete flows

In this paper, we investigate the evolution of joint invariants under invariant geometric flows using the theory of equivariant moving frames and the induced invariant discrete variational complex. For certain arc length preserving planar curve flows invariant under the special Euclidean group , the special linear group , and the semidirect group , we find that the induced evolution of the discrete curvature satisfies the differential‐difference mKdV, KdV, and Burgers' equations, respectively. These three equations are completely integrable, and we show that a recursion operator can be constructed by precomposing the characteristic operator of the curvature by a certain invariant difference operator. Finally, we derive the constraint for the integrability of the discrete curvature evolution to lift to the evolution of the discrete curve itself.     

Dynamics modeling and robotic-assist,leader-follower control of tractor convoys

This paper proposes a generalized dynamics model and a leader-follower control architecture for skid-steered tracked vehicles towing polar sleds. The model couples existing formulations in the literature for the powertrain components with the vehicle-terrain interaction to capture the salient features of terrain trafficability and predict the vehicles response. This coupling is essential for making realistic predictions of the vehicles traversing capabilities due to the power-load relationship at the engine output. The objective of the model is to capture adequate fidelity of the powertrain and off-road vehicle dynamics while minimizing the computational cost for model based design of leader-follower control algorithms. The leader-follower control architecture presented proposes maintaining a flexible formation by using a look-ahead technique along with a way point following strategy. Results simulate one leader-follower tractor pair where the leader is forced to take an abrupt turn and experiences large oscillations of its drawbar arm indicating potential payload instability. However, the follower tractor maintains the flexible formation but keeps its payload stable. This highlights the robustness of the proposed approach where the follower vehicle can reject errors in human leader driving.     

A Runge-Kutta–based WENO reconstruction on moving mesh for compressible fluids

In this paper, we construct a high-order moving mesh method based on total variation diminishing Runge-Kutta and weighted essential nonoscillatory reconstruction for compressible fluid system. Beginning with the integral form of fluid system, we get the semidiscrete system with an arbitrary mesh velocity. We use weighted essential nonoscillatory reconstruction to get the space accuracy on moving meshes, and the time accuracy is obtained by modified Runge-Kutta method; the mesh velocity is determined by moving mesh method. One- and two-dimensional numerical examples are presented to demonstrate the efficient and accurate performance of the scheme.     

A proper orthogonal decomposition variational multiscale meshless interpolating element-free Galerkin method for incompressible magnetohydrodynamics flow

In the recent decade, the meshless methods have been handled for solving most of PDEs due to easiness of the meshless methods. One of the popular meshless methods is the element-free Galerkin (EFG) method that was first proposed for solving some problems in the solid mechanics. The test and trial functions of the EFG are based on the special basis. Recently, some modifications have been developed to improve the EFG method. One of these improvements is the variational multiscale EFG procedure. In the current article, the shape functions of interpolation moving least squares approximation have been applied to the variational multiscale EFG technique for solving the incompressible magnetohydrodynamics flow. In order to reduce the elapsed CPU time of simulation, we employ a reduced-order model based on the proper orthogonal decomposition technique. The current combination can be referred to as the reduced-order variational multiscale EFG technique. To illustrate the reduction in CPU time used as well as the efficiency of the proposed method, we applied it for the two-dimensional cases.     

基于粒子滤波的运动目标光电定位仿真研究

运动目标的光电定位不能像静止目标那样简单做均值滤波，鉴于此，引入粒子滤波算法，它不仅可以应用于线性系统，而且还适用于非线性系统。结合光电定位需求，详细推导了计算公式及初值和参数选取公式，对只含测量噪声以及含有测量和运动噪声等的海面运动目标光电无源定位算法进行了仿真计算，验证了算法的有效性，讨论了噪声强度对滤波效果的影响，滤波参数选择对滤波效果的影响，目标运动方式对滤波跟随性的影响，重采样算法对滤波效果的影响等。所得结论为:粒子滤波可用于运动目标光电定位过程，可有效降低定位误差；粒子滤波算法具有较强鲁棒性，适用于噪声较大、目标运动形态变化大等情况。     

Meshless acoustic analysis using a weakly singular Burton-Miller boundary integral formulation

The Burton-Miller boundary integral formulation is solved by a complex variable boundary element-free method (CVBEFM) for the boundary-only meshless analysis of acoustic problems with arbitrary wavenumbers. To regularize both strongly singular and hypersingular integrals and to avoid the computation of the solid angle and its normal derivative, a weakly singular Burton-Miller formulation is derived by considering the normal derivative of the solid angle and adopting the singularity subtraction procedures. To facilitate the implementation of the CVBEFM and the approximation of gradients of the boundary variables, a stabilized complex variable moving least-square approximation is selected in the meshless discretization procedure. The results show the accuracy and efficiency of the present CVBEFM and reveal that the method can produce satisfactory results for all wavenumbers, even for extremely large wavenumbers such as k = 10 000.     

Thermoelastic solutions for thermal distributions moving over thin slim rod under memory-dependent three-phase lag magneto-thermoelasticity

Abstract

Enlightened by the Caputo fractional derivative, the present study deals with a novel mathematical model of generalized thermoelasticity to investigate the transient phenomena due to the influence of magnetic field and moving heat source in a rod in the context of three-phase lag (TPL) theory of thermoelasticity. Both ends of the rod are fixed and heat insulated. Employing Laplace transform as a tool, the problem has been transformed into the space-domain and solved analytically. Finally, solutions in the real-time domain are obtained by applying the inverse Laplace transform. Numerical calculation for stress, displacement, and temperature within the rod is carried out and displayed graphically. The effect of moving heat source speed on temperature, stress, and temperature is studied. It is found from the distributions that the temperature, thermally induced displacement and stress of the rod are found to decrease at large source speed. For the better understanding of the effect of moving heat source on all the distributions, three animations are added.     

基于无人机多光谱图像的土壤水分检测方法研究

以表层土壤为对象,探究土壤的多光谱反射率与土壤水分含量相关性,进行基于无人机多光谱图像的土壤水分含量预测模型方法的探究。选取中国农业大学通州实验站为研究区域,实地采集试验田的土壤样本100组,按照一定梯度配制土壤含水量,配成的土壤含水率为10%~50%之间,土壤含量的真实值采用土壤烘干法进行测定。多光谱相机灵巧便捷,可搭载在无人机上对土壤进行监测。将RedEdged-M型多光谱相机搭载在Phantom 3型无人机上,选择阳光充足的采集环境,实时采集土壤样本的多光谱图像,建立土壤多光谱信息与水分含量之间的模型。利用处理光谱数据的ENVI5.3软件提取土壤样本多光谱信息,以多光谱相机自带的标准白板反射率为100%,计算出土壤样本在蓝、绿、红、红边、近红外五个波段的光谱反射率。采用BP神经网络算法、支持向量机算法、偏最小二乘算法分别建立基于无人机多光谱图像的土壤水分含量的预测模型。以80组土壤样本数据作为训练集,建立基于多光谱图像的土壤水分含量预测模型。采用莱文贝格-马夸特算法对BPNN进行改进,提高了其训练速度,当网络结构为5-10-1时,训练效果最好,本文选择该网络结构;SVM采取高斯核函数,当参数为0.56时,模型效果最好。本研究采用归一化均方根误差(NRMSE)和决策系数(R 2)对三种土壤水分含量的预测模型进行定量对比。以20组土壤样本数据作为测试集,结果可知,基于BP神经网络土壤水分含量预测模型的NRMSE为0.268,R 2为0.872;基于支持向量机的土壤水分含量预测模型的NRMSE为0.298,R 2为0.821;基于偏最小二乘土壤水分含量预测模型的NRMSE为0.316,R 2为0.789。对三种模型分析可知,基于BPNN的土壤水分含量预测模型效果均较好。结果可知,土壤的光谱反射率与含水率间存在较密切的相关性,将多光谱相机搭载在无人机上可以对土壤水分含量进行有效的实时监测,对监测土壤墒情提供技术支持和理论支撑。