超声速流动中非线性EASM湍流模式应用研究

针对超声速复杂流动区域精确模拟的需要，发展了基于k-ω可压缩修正形式的非线性显式代数雷诺应力模式（EASM），提高了该模式对超声速复杂流动的数值模拟精度。通过对二维超声速凹槽和三维双椭球的数值计算表明，与SA和SST常规线性涡黏性湍流模式比较，非线性的EASM模式对大分离以及剪切层流动结构的刻画能力更精细，对剪切层再附区的压力及摩擦系数分布模拟更加精确；EASM模式能够准确地模拟二次激波引起的压强和热流分布情况。     

水下对转螺旋桨流致辐射噪声机理与预报方法*

水下对转螺旋桨流致辐射噪声的预报对于水下目标的特征提取和分类识别具有重要意义。由桨叶的旋转引起的湍流场是水下对转螺旋桨流致辐射噪声的源场。分述了水下对转螺旋桨湍流边界层脉动、旋转干涉效应和空化效应引发的水动力噪声机制和研究进展,比较了目前工程应用中的3种对转螺旋桨流致辐射噪声预报方法的特点。在分析对转螺旋桨流致辐射噪声数值预报难点的基础上,综述了对转螺旋桨流致辐射噪声计算方法的研究进展,指出间接数值模拟方法是工程中进行对转螺旋桨流致辐射噪声预报的有效方法。     

A C0-WEAK GALERKIN FINITE ELEMENT METHOD FOR THE TWO-DIMENSIONAL NAVIER-STOKES EQUATIONS IN STREAM-FUNCTION FORMULATION

We propose and analyze a $C^0$-weak Galerkin (WG) finite element method for the numerical solution of the Navier-Stokes equations governing 2D stationary incompressible flows. Using a stream-function formulation, the system of Navier-Stokes equations is reduced to a single fourth-order nonlinear partial differential equation and the incompressibility constraint is automatically satisfied. The proposed method uses continuous piecewise-polynomial approximations of degree $k+2$ for the stream-function $\psi$ and discontinuous piecewise-polynomial approximations of degree $k+1$ for the trace of $\nabla\psi$ on the interelement boundaries. The existence of a discrete solution is proved by means of a topological degree argument, while the uniqueness is obtained under a data smallness condition. An optimal error estimate is obtained in $L^2$-norm, $H^1$-norm and broken $H^2$-norm. Numerical tests are presented to demonstrate the theoretical results.     

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.     

Application of affinity capillary electrophoresis for charge heterogeneity profiling of biopharmaceuticals

Charge heterogeneity profiling is important for the quality control (QC) of biopharmaceuticals. Because of the increasing complexity of these therapeutic entities [1], the development of alternative analytical techniques is needed. In this work, flow‐through partial‐filling affinity capillary electrophoresis (FTPFACE) has been established as a method for the analysis of a mixture of two similar monoclonal antibodies (mAbs). The addition of a specific ligand results in the complexation of one mAb in the co‐formulation, thus changing its migration time in the electric field. This allows the characterization of the charged variants of the non‐shifted mAb without interferences. Adsorption of proteins to the inner capillary wall has been circumvented by rinsing with guanidine hydrochloride before each injection. The presented FTPFACE approach requires only very small amounts of ligands and provides complete comparability with a standard CZE of a single mAb.     

正交异性复合材料储液容器入水过程数值分析

采用多物质ALE流固耦合算法,对柔性储液容器的入水问题进行了仿真研究。该方法可以方便地考虑多种流体同时与结构耦合。首先,通过模拟平板恒速入水试验,验证了该方法的有效性,并确定了空气的影响。进而,仿真分析了复合材料柔性容器入水时的动态响应行为,结果表明,底部"空气垫"及容器内空气对冲击载荷具有一定的缓冲作用,容器底部形状、储液量等因素也具有较大影响。仿真得到的响应规律为容器的优化设计提供了重要的理论参考。     

关于不定方程x3-1=749y2

利用同余式、平方剩余、Pell方程的解的性质、递归序列证明了:不定方程x³-1=749y²仅有整数解(x,y)=(1,0).     

Shape morphing and topology optimization of fluid channels by explicit boundary tracking

An integrated shape morphing and topology optimization approach based on the deformable simplicial complex methodology is developed to address Stokes and Navier‐Stokes flow problems. The optimized geometry is interpreted by a set of piecewise linear curves embedded in a well‐formed triangular mesh, resulting in a physically well‐defined interface between fluid and impermeable regions. The shape evolution is realized by deforming the curves while maintaining a high‐quality mesh through adaption of the mesh near the structural boundary, rather than performing global remeshing. Topological changes are allowed through hole merging or splitting of islands. The finite element discretization used provides smooth and stable optimized boundaries for simple energy dissipation objectives. However, for more advanced problems, boundary oscillations are observed due to conflicts between the objective function and the minimum length scale imposed by the meshing algorithm. A surface regularization scheme is introduced to circumvent this issue, which is specifically tailored for the deformable simplicial complex approach. In contrast to other filter‐based regularization techniques, the scheme does not introduce additional control variables, and at the same time, it is based on a rigorous sensitivity analysis. Several numerical examples are presented to demonstrate the applicability of the approach.     

Implicit–explicit schemes for nonlinear nonlocal equations with a gradient flow structure in one space dimension

Nonlinear convection–diffusion equations with nonlocal flux and possibly degenerate diffusion arise in various contexts including interacting gases, porous media flows, and collective behavior in biology. Their numerical solution by an explicit finite difference method is costly due to the necessity of discretizing a local spatial convolution for each evaluation of the convective numerical flux, and due to the disadvantageous Courant–Friedrichs–Lewy (CFL) condition incurred by the diffusion term. Based on explicit schemes for such models devised in the study of Carrillo et al. a second‐order implicit–explicit Runge–Kutta (IMEX‐RK) method can be formulated. This method avoids the restrictive time step limitation of explicit schemes since the diffusion term is handled implicitly, but entails the necessity to solve nonlinear algebraic systems in every time step. It is proven that this method is well defined. Numerical experiments illustrate that for fine discretizations it is more efficient in terms of reduction of error versus central processing unit time than the original explicit method. One of the test cases is given by a strongly degenerate parabolic, nonlocal equation modeling aggregation in study of Betancourt et al. This model can be transformed to a local partial differential equation that can be solved numerically easily to generate a reference solution for the IMEX‐RK method, but is limited to one space dimension.