基于二阶锥规划的宽带波束形成器设计

宽带波束形成器有两种典型的实现方式,分别基于FIR滤波器和基于长方形阵列,将这两种波束形成器的波束响应表达为统一的形式,并根据参考波束的选择情况,为该类波束优化问题的求解建立了数学模型,该数学模型属于二阶锥约束问题,可以用二阶锥规划方法进行求解。仿真结果表明,提出的求解宽带恒定束宽响应波束形成器权系数的数学模型具有通用性,并且把二阶锥规划方法运用到宽带波束形成器优化设计中,其约束控制灵活,问题求解方便,设计结果精确。     

新型线性先行GA算法在油品调和中的应用

油品调和是生产成品油的关键环节,调和配方决定了利润的最终大小。配方优化模型中物料平衡约束、产品产量约束属于线性约束,但质量指标约束中存在的非线性约束,使整个求解问题成为非线性约束求解问题。针对炼油厂油品调和复杂非线性约束配方优化问题,提出了一种新的线性先行GA算法。通过将线性约束与非线性约束拆分开,先求解凸多面体的全部顶点和极方向,解出满足线性约束的所有可行解所在的可行区域,进而用GA算法在此可行域内,通过优化凸组合的系数,从而达到搜索整个线性解区间的目的,完成非线性优化模型的求解。仿真结果表明,新算法大大缩小了GA算法的搜索空间和需要处理的约束,能快速的获得理想的调和配方,使利润最大并保证对质量卡边的要求。     

三维非定常不可压涡量—速度Navier—Stokes方程组的有限差分法

提出了一种数值求解三维非定常涡量—速度形式的不可压Navier-Stokes方程组的有限差分方法,该方法在空间方向上具有二阶精度,并且系数矩阵具有对角占优性,因此适合高雷诺数问题的数值求解．同时,给出了适合的二阶涡量边界条件．通过对有精确解的狄利克雷边值问题和典型的驱动方腔流问题的数值实验,验证了本文格式的精确性、稳定性和有效性．     

基于神经网络的冗余度TT-VGT机器人的运动学求解

应用BP神经网络对冗余度TT-VGT机器人的位姿正解进行训练学习,进而求解机器人 的位姿反解问题．根据网络模型求得机器人的一、二阶影响系数,应用神经网络求解雅可比 矩阵的伪逆．并对七重四面体的变几何桁架机器人进行了仿真计算．     

伯格谱估计算法的一种改进

在实数据情况下，由伯格算法得到的谱有时会出现谱分裂和谱峰偏移。已有文献对这种现象产生的原因作了各种分析和猜测。其中之一认为由于一阶反射系数的计算存在误差，导致预测误差滤波器（PEF）其余系数计算的误差。根据这种分析，本文提出了一种改进算法。该算法不是直接利用伯格算法计算一阶反射系数，而是按二阶PEF输出总功率最小的原则，求解二阶PEF系数，再求得一附反射系数。在获得二阶PEF系数之后，仍用伯格递推算法求其他高阶系数。结果，在几乎不增加计算量的情况下，使谱估计的性能得到明显改善。对信号长度为1/4信号周期的奇数倍且初相为л/4时的最坏情况，在计算机上进行了计算。结果表明，谱峰偏移大为减小。     

基于多目标优化的汽油调和优化模型及其应用研究

汽油的调和方案是提高炼油厂经济效益的关键因素.通过现有汽油调和优化模型的分析,本文结合汽油辛烷值的预测模型1DPE模型,提出了1种以调和成本和辛烷值过剩最少为目标、将经验知识也作为约束条件的多目标优化汽油调和优化模型,该模型可应用于不同标号的汽油产品的调和,根据模型的非线性特征,采用遗传算法进行求解,最后,通过某炼油厂的1个应用实例验证了上述成品油调和优化模型及其求解方法的有效性,为炼油厂减少企业成本、保证产品合格率提供了新途径.     

采用渐近法研究负载钢丝绳的固有振动

结合Liouville-Green变换,改进了求解变系数二阶线性齐次方程的渐近法.并采用改进后的渐近法研究了负载钢丝绳的固有振动问题,推导出了其固有振动的近似频率特征方程.实例计算表明,改进后的渐近法不但比Bessel函数法计算简便,而且计算精度也非常高.     

数值求解优化问题在活动轮廓模型上的应用

针对活动轮廓模型图像分割过程中迭代次数多,分割速度慢的问题,提出一种高阶的数值求解方法。分析活动轮廓模型中基于全局信息的CV模型,以及基于局部信息的LBF模型,LIF模型。使用二阶、三阶Runge-Kutta方法,原始Euler方法对模型进行数值求解实验对比分析。并对LBF模型中平滑项系数,时间步长的选取进行讨论。通过对非同质图像、同质图像的实验结果分析表明,所采用的数值方法能够有效地提高数值收敛精度、减少迭代次数、计算效率高。对不同系数和时间步长,数值方法也能表现出较好的稳定性。     

二阶离散线性时变系统的一种稳定性判据

本文对二阶离散线性时变稳定性问题进行了研究。在离散线性时变系统Lyapunov稳定理论的基础上，给出二阶离散线性时变系统稳定的充分条件。即系统系数满足系数冻结条件和系数变化条件。仿真验证了结论的正确性。     

一阶、二阶导数耦合的紧致差分格式及其应用

在数值计算中可能遇到求解一阶导数和二阶导数耦合的微分方程,为了能用紧致差分格式进行计算,针对这样的方程,建立了考虑一阶、二阶导数耦合的紧致差分格式,利用这一方法可以直接对方程进行离散求解。通过具体算例,验证该类紧致差分格式的优越性,还将这类紧致差分格式运用到求解二维偏微分方程中。     

A simple implementation of the immersed interface methods for Stokes flows with singular forces

In this paper, we introduce a simple version of the immersed interface method (IIM) for Stokes flows with singular forces along an interface. The numerical method is based on applying the Taylor’s expansions along the normal direction and incorporating the solution and its normal derivative jumps into the finite difference approximations. The fluid variables are solved in a staggered grid, and a new accurate interpolating scheme for the non-smooth velocity has been developed. The numerical results show that the scheme is second-order accurate.     

Unconditionally stable ADI scheme of higher-order for linear hyperbolic equations

An unconditionally stable alternating direction implicit (ADI) method of higher-order in space is proposed for solving two- and three-dimensional linear hyperbolic equations. The method is fourth-order in space and second-order in time. The solution procedure consists of a multiple use of one-dimensional matrix solver which produces a computational cost effective solver. Numerical experiments are conducted to compare the new scheme with the existing scheme based on second-order spatial discretization. The effectiveness of the new scheme is exhibited from the numerical results.     

A direct-forcing immersed boundary method for the thermal lattice Boltzmann method

In this study, a direct-forcing immersed boundary method (IBM) for thermal lattice Boltzmann method (TLBM) is proposed to simulate the non-isothermal flows. The direct-forcing IBM formulas for thermal equations are derived based on two TLBM models: a double-population model with a simplified thermal lattice Boltzmann equation (Model 1) and a hybrid model with an advection–diffusion equation of temperature (Model 2). As an interface scheme, which is required due to a mismatch between boundary and computational grids in the IBM, the sharp interface scheme based on second-order bilinear and linear interpolations (instead of the diffuse interface scheme, which uses discrete delta functions) is adopted to obtain the more accurate results. The proposed methods are validated through convective heat transfer problems with not only stationary but also moving boundaries – the natural convection in a square cavity with an eccentrically located cylinder and a cold particle sedimentation in an infinite channel. In terms of accuracy, the results from the IBM based on both models are comparable and show a good agreement with those from other numerical methods. In contrast, the IBM based on Model 2 is more numerically efficient than the IBM based on Model 1.     

A numerical model for multiphase flow based on the GMPPS formulation. Part I: Kinematics

The CFD modeling of two-dimensional multiphase flows is a useful tool in industry, although accurate modeling itself remains a difficult task. One of the difficulties is to track the complicated topological deformations of the interfaces between different phases. This paper describes a marker-particle method designed to track fluid interfaces for fluid flows of at least three phases. The interface-tracking scheme presented in this paper is the first part of a series of papers presenting our complete model based on a one-field Godunov marker-particle projection scheme (GMPPS). In this part, we shall focus on the presentation of the interface-tracking scheme and the kinematic tests we conducted to examine the scheme’s ability to accurately track interfacial movements typified by vorticity-induced stretching and tearing of the interface. Our test results show that for a set of carefully designed and commonly used error measures, relative percentage errors never exceed 2% for all of the tests and grid sizes considered, provided a sufficient number of marker particles are used. We shall also demonstrate that the method is of second-order accuracy and the interface transition width remains constant never exceeding three cell widths.     

Conservative compact difference scheme for the Zakharov–Rubenchik equations

In this article, a compact difference scheme is investigated to solve the Zakharov–Rubenchik equations in one dimension. The new scheme is proved to conserve the total mass and energy in the discrete sense. Rigorous error estimates are established for the new method with the help of an induction argument in energy space which show that the new scheme has second-order accuracy in time and fourth-order accuracy in space. Extensive numerical results are provided to verify our theoretical analysis, and show the accuracy and efficiency of the new scheme.     

An Explicit Implicit Scheme for Cut Cells in Embedded Boundary Meshes

We present a new mixed explicit implicit time stepping scheme for solving the linear advection equation on a Cartesian cut cell mesh. We use a standard second-order explicit scheme on Cartesian cells away from the embedded boundary. On cut cells, we use an implicit scheme for stability. This approach overcomes the small cell problem—that standard schemes are not stable on the arbitrarily small cut cells—while keeping the cost fairly low. We examine several approaches for coupling the schemes in one dimension. For one of them, which we refer to as flux bounding, we can show a TVD result for using a first-order implicit scheme. We also describe a mixed scheme using a second-order implicit scheme and combine both mixed schemes by an FCT approach to retain monotonicity. In the second part of this paper, extensions of the second-order mixed scheme to two and three dimensions are discussed and the corresponding numerical results are presented. These indicate that this mixed scheme is second-order accurate in \(L^1\) and between first- and second-order accurate along the embedded boundary in two and three dimensions.     

GIFT密码算法的二阶门限实现及其安全性评估

目前已知GIFT算法的防护方案仅能抗一阶功耗攻击。为了使该算法能抵御高阶功耗攻击,利用GIFT算法的数学结构,结合门限实现方法,构造了GIFT算法的二阶门限实现方案;针对非线性部件S盒的数学特征,分别构建（3,9）,（6,7）以及（5,10）等二阶门限实现方案,并选取所需硬件资源最少的（3,9）方案,将其在FPGA平台下实现。在实现时,为了更好地分析对比所耗硬件资源,使用并行和串行的硬件实现方法。结果表明,在NanGate 45nm工艺库下,并行实现消耗的总面积为12 043 GE,串行实现消耗的总面积为6 373 GE。通过采集实际功耗曲线,利用T-test对该二阶门限实现方案进行侧信息泄露评估,证实了该二阶门限实现方案的安全性。     

Study on a grid interface algorithm for solutions of incompressible Navier–Stokes equations

Grid interface treatment is a crucial issue in solving unsteady, three-dimensional, incompressible Navier–Stokes equations by domain decomposition methods. Recently, a mass flux based interpolation (MFBI) interface algorithm was proposed for Chimera grids [Tang HS, Jones SC, Sotiropoulos F. An overset grid method for 3D unsteady incompressible flows. J Comput Phys, 2003;191:567–600] and it has been successfully applied to a variety of flows. MFBI determines velocity and pressure at grid interfaces by mass conservation and interpolation, and it is easy to implement. Compared with the commonly used standard interpolation, which directly interpolates velocity as well as pressure, the proposed interface algorithm gives fewer solution oscillations and faster convergence rates. This paper makes a study on MFBI. Starting with discussions about grid connectivity, it is shown that MFBI is second-order accurate for mass flux across grid interface. It is also derived that the scheme provides second-order accuracy for momentum flux. In addition, another version of MFBI is presented. At last, numerical examples are presented to demonstrate that MFBI honors mass flux balance at grid interfaces and it leads to second-order accurate solutions.     

A new ADI technique for two-dimensional parabolic equation with an integral condition

Nonclassical parabolic initial-boundary value problems arise in the study of several important physical phenomena. This paper presents a new approach to treat complicated boundary conditions appearing in the parabolic partial differential equations with nonclassical boundary conditions. A new fourth-order finite difference technique, based upon the Noye and Hayman (N-H) alternating direction implicit (ADI) scheme, is used as the basis to solve the two-dimensional time dependent diffusion equation with an integral condition replacing one boundary condition. This scheme uses less central processor time (CPU) than a second-order fully implicit scheme based on the classical backward time centered space (BTCS) method for two-dimensional diffusion. It also has a larger range of stability than a second-order fully explicit scheme based on the classical forward time centered space (FTCS) method. The basis of the analysis of the finite difference equations considered here is the modified equivalent partial differential equation approach, developed from the 1974 work of Warming and Hyeet. This allows direct and simple comparison of the errors associated with the equations as well as providing a means to develop more accurate finite difference methods. The results of numerical experiments for the new method are presented. The central processor times needed are also reported. Error estimates derived in the maximum norm are tabulated.     

Parameter-uniform hybrid numerical scheme for time-dependent convection-dominated initial-boundary-value problems

In this paper, we propose a numerical scheme which is almost second-order spatial accurate for a one-dimensional singularly perturbed parabolic convection-diffusion problem exhibiting a regular boundary layer. The proposed numerical scheme consists of classical backward-Euler method for the time discretization and a hybrid finite difference scheme for the spatial discretization. We analyze the scheme on a piecewise-uniform Shishkin mesh for the spatial discretization to establish uniform convergence with respect to the perturbation parameter. Numerical results are presented to validate the theoretical results.