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.     

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.     

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.     

Preparative separation of immunoglobulins from bovine colostrum by continuous divergent-flow electrophoresis

Immunoglobulins in bovine colostrum were separated and fractionated from other proteins using the method and instrumentation developed in our laboratory. The proposed separation was based on bidirectional isotachophoresis/moving boundary electrophoresis with electrofocusing of the analytes in a pH gradient from 3.9 to 10.1. The preparative instrumentation included the trapezoidal non-woven fabric that served as separation space with divergent continuous flow. The defatted and casein precipitate-free colostrum supernatant was loaded directly into the instrument without any additional colostrum pre-preparation. Immunoglobulin G was fractionated from other immune proteins such as bovine serum albumin, β-lactoglobulin, and α-lactalbumin, and was continuously collected in separated fractions over 3 h. The fractions were further processed, and isolated immunoglobulin G in the liquid fractions was confirmed by sodium dodecyl sulfate-polyacrylamide gel electrophoresis and by re-focusing in gel isoelectric focusing. Separated immunoglobulin G was detected in seven fractions by sodium dodecyl sulfate-polyacrylamide gel electrophoresis with a gradually decreased concentration in the fractions. Re-focusing of the proteins in the fractions by gel isoelectric focusing revealed multiple separated zones of immunoglobulin G with the isoelectric point values covering the range from 5.4 to 7.2. Each fraction contained distinct zones with gradually increased isoelectric point values and decreased concentrations from fraction to fraction.     

基于高通量解析算法的复杂样品重叠气相色谱-质谱信号的快速分析

化学计量学算法为重叠气相色谱-质谱(GC-MS)信号的解析提供了有效手段,但其在计算过程中一般需要将数据进行分段处理,然后只对信号的某些区间进行解析,难以实现真正意义上的高通量分析。该文结合移动窗口目标转换因子分析(MWTTFA)和非负免疫算法(NNIA),建立了一种高通量解析方法。首先,根据所有可能存在的目标组分的标准质谱信息,利用MWTTFA检验复杂信号中存在的组分,并确定目标组分的质谱信息和洗脱时间区域。以得到的质谱信息作为后续计算的输入值,利用NNIA解析得到相应的色谱信息。采用快速升温程序对17种和42种农药混合标准样品的GC-MS信号进行分析,利用所建立的方法可在10 min内得到全部组分的色谱和质谱信息。     

具有N-各向异性拉普拉斯算子的超线性问题的先验界

通过改进Brezis和Merle的方法,结合Moser-Trudinger不等式,移动平面方法及比较原理,得到了方程-Q_Nu=f(u),u∈W_0~(1,N)(Ω)的正解的先验界,其中Ω是R~N中的一个有界光滑区域,非线性项f至多具有指数型增长.