A hybrid subcell-remapping algorithm for staggered multi-material arbitrary Lagrangian-Eulerian methods

A new flux-based hybrid subcell-remapping algorithm for staggered multimaterial arbitrary Lagrangian-Eulerian (MMALE) methods is presented. This new method is an effective generalization of the original subcell-remapping method to the multi-material regime (LOUBERE, R. and SHASHKOV,M. A subcell remapping method on staggered polygonal grids for arbitrary-Lagrangian-Eulerian methods. Journal of Computational Physics, 209, 105–138 (2005)). A complete remapping procedure of all fluid quantities is described detailedly in this paper. In the pure material regions, remapping of mass and internal energy is performed by using the original subcell-remapping method. In the regions near the material interfaces, remapping of mass and internal energy is performed with the intersection-based fluxes where intersections are performed between the swept regions and pure material polygons in the Lagrangian mesh, and an approximate approach is then introduced for constructing the subcell mass fluxes. In remapping of the subcell momentum, the mass fluxes are used to construct the momentum fluxes by multiplying a reconstructed velocity in the swept region. The nodal velocity is then conservatively recovered. Some numerical examples simulated in the full MMALE regime and several purely cyclic remapping examples are presented to prove the properties of the remapping method.

An efficient algorithm to compute transient pressure responses of slanted wells with arbitrary inclination in reservoirs

Compared with vertical and horizontal wells, the solution and computation of transient pressure responses of slanted wells are more complex. Vertical and horizontal wells are both simplified cases of slanted wells at particular inclination, so the model for slanted wells is more general and more complex than other models for vertical and horizontal wells. Many authors have studied unsteady-state flow of fluids in slanted wells and various solutions have been proposed. However, until now, few of the published results pertain to the computational efficiency. Whether in the time domain or in the Laplace domain, the computation of integration of complex functions is necessary in obtaining pressure responses of slanted wells, while the computation of the integration is complex and time-consuming. To obtain a perfect type curve the computation time is unacceptable even with an aid of high-speed computers. The purpose of this paper is to present an efficient algorithm to compute transient pressure distributions caused by slanted wells in reservoirs. Based on rigorous derivation, the transient pressure solution for slanted wells of any inclination angle is presented. Assuming an infinite-conductivity wellbore, the location of the equivalent-pressure point is determined. More importantly, according to the characteristics of the integrand in a transient pressure solution for slanted wells, the whole integral interval is partitioned into several small integral intervals, and then the method of variable substitution and the variable step-size piecewise numerical integration are employed. The amount of computation is significantly reduced and the computational efficiency is greatly improved. The algorithm proposed in this paper thoroughly solved the difficulty in the efficient and high-speed computation of transient pressure distribution of slanted wells with any inclination angle.

Single-Step Arbitrary Lagrangian-Eulerian Discontinuous Galerkin Method for 1-D Euler Equations

We propose an explicit,single-step discontinuous Galerkin method on moving grids using the arbitrary Lagrangian-Eulerian approach for one-dimensional Euler equations.The grid is moved with the local fluid velocity modified by some smoothing,which is found to con-siderably reduce the numerical dissipation introduced by Riemann solvers.The scheme preserves constant states for any mesh motion and we also study its positivity preservation property.Local grid refinement and coarsening are performed to maintain the mesh qual-ity and avoid the appearance of very small or large cells.Second,higher order methods are developed and several test cases are provided to demonstrate the accuracy of the proposed scheme.

Negative Norm Estimates for Arbitrary Lagrangian-Eulerian Discontinuous Galerkin Method for Nonlinear Hyperbolic Equations

In this paper,we present the negative norm estimates for the arbitrary Lagrangian-Eulerian discontinuous Galerkin(ALE-DG)method solving nonlinear hyperbolic equations with smooth solutions.The smoothness-increasing accuracy-conserving(SIAC)filter is a post-processing technique to enhance the accuracy of the discontinuous Galerkin(DG)solutions.This work is the essential step to extend the SIAC filter to the moving mesh for nonlinear problems.By the post-processing theory,the negative norm estimates are vital to get the superconvergence error estimates of the solutions after post-processing in the L2 norm.Although the SIAC filter has been extended to nonuniform mesh,the analysis of fil-tered solutions on the nonuniform mesh is complicated.We prove superconvergence error estimates in the negative norm for the ALE-DG method on moving meshes.The main dif-ficulties of the analysis are the terms in the ALE-DG scheme brought by the grid velocity field,and the time-dependent function space.The mapping from time-dependent cells to reference cells is very crucial in the proof.The numerical results also confirm the theoreti-cal proof.

A Novel Memory Compress Algorithm for Arbitrary Waveform Generator

A memory compress algorithm for 12\|bit Arbitrary Waveform Generator (AWG) is presented and optimized. It can compress waveform memory for a sinusoid to 16×13bits with a Spurious Free Dynamic Range (SFDR) 90.7dBc (1/1890 of uncompressed memory at the same SFDR) and to 8×12bits with a SFDR 79dBc. Its hardware cost is six adders and two multipliers. Exploiting this memory compress technique makes it possible to build a high performance AWG on a chip.

An Algorithm for Partitioning Polygons into Convex Parts

An algorithm for partitioning arbitrary simple polygons into a number of convex parts was presented. The concave vertices were determined first, and then they were moved by using the method connecting the concave vertices with the vertices of falling into its region B,so that the primary polygon could be partitioned into two subpolygons. Finally, this method was applied recursively to the subpolygons until all the concave vertices were removed. This algorithm partitions the polygon into O(l) convex parts, its time complexity is max(O(n),O(l 2)) multiplications, where n is the number of vertices of the polygon and l is the number of the concave vertices.

Arbitrary Lagrangian‑Eulerian Discontinuous Galerkin Methods for KdV Type Equations

In this paper,several arbitrary Lagrangian-Eulerian discontinuous Galerkin(ALE-DG)methods are presented for Korteweg-de Vries(KdV)type equations on moving meshes.Based on the L^(2) conservation law of KdV equations,we adopt the conservative and dissipative numerical fuxes for the nonlinear convection and linear dispersive terms,respectively.Thus,one conservative and three dissipative ALE-DG schemes are proposed for the equations.The invariant preserving property for the conservative scheme and the corresponding dissipative properties for the other three dissipative schemes are all presented and proved in this paper.In addition,the L^(2)-norm error estimates are also proved for two schemes,whose numerical fuxes for the linear dispersive term are both dissipative type.More precisely,when choosing the approximation space with the piecewise kth degree polynomials,the error estimate provides the kth order of convergence rate in L^(2)-norm for the scheme with the conservative numerical fuxes applied for the nonlinear convection term.Furthermore,the(k+1∕2)th order of accuracy can be proved for the ALE-DG scheme with dissipative numerical fuxes applied for the convection term.Moreover,a Hamiltonian conservative ALE-DG scheme is also presented based on the conservation of the Hamiltonian for KdV equations.Numerical examples are shown to demonstrate the accuracy and capability of the moving mesh ALE-DG methods and compare with stationary DG methods.

WinoNet:Reconfigurable look-up table-based Winograd accelerator for arbitrary precision convolutional neural network inference

To solve the hardware deployment problem caused by the vast demanding computational complexity of convolutional layers and limited hardware resources for the hardware network inference,a look-up table(LUT)-based convolution architecture built on a field-programmable gate array using integer multipliers and addition trees is used.With the help of the Winograd algorithm,the optimization of convolution and multiplication is realized to reduce the computational complexity.The LUT-based operator is further optimized to construct a processing unit(PE).Simultaneously optimized storage streams improve memory access efficiency and solve bandwidth constraints.The data toggle rate is reduced to optimize power consumption.The experimental results show that the use of the Winograd algorithm to build basic processing units can significantly reduce the number of multipliers and achieve hardware deployment acceleration,while the time-division multiplexing of processing units improves resource utilization.Under this experimental condition,compared with the traditional convolution method,the architecture optimizes computing resources by 2.25 times and improves the peak throughput by 19.3 times.The LUT-based Winograd accelerator can effectively solve the deployment problem caused by limited hardware resources.

基于流固耦合的救生筏海锚阻力性能分析

救生筏是提高航行安全的海上救生设备.作为救生筏附件,海锚能有效降低恶劣海况下救生筏的漂移速度,提高其稳定性.为对海锚进行阻力性能分析,基于Star-CCM+与LS-DYNA软件对其水中工作过程及阻力性能进行研究,建立了可透气柔性体海锚的流固耦合动力学模型与不考虑救生筏受风影响变形的刚性体模型.对救生筏在海上典型风速条件下的漂移速度及海锚的减速效果进行了数值模拟,对不同救生筏漂移速度下海锚的阻力性能和结构变形进行了综合分析.修正了救生筏筏体总阻力的前体尾流影响与刚-柔性体海锚转换影响.柔性体海锚的阻力拟合公式与文献中的相似性证明了数值模拟的可靠性.结果表明海锚可有效降低救生筏的漂移速度,提高救生筏被捕救的可能性.

A Vertex-Centered Arbitrary Lagrangian-Eulerian Finite Volume Method with Sub-Cells for Two-Dimensional Compressible Flow

In this paper,we present a new vertex-centered arbitrary LagrangianEulerian(ALE)finite volume scheme for two-dimensional compressible flow.In our scheme,the momentum equation is discretized on the vertex control volume,while the mass equation and the energy equation are discretized on the sub-cells which are included in the vertex control volume.We attain the average of the fluid velocity on the vertex control volume directly by solving the conservation equations.Then we can obtain the fluid velocity at vertex with the reconstructed polynomial of the velocity.This fluid velocity is chosen as the mesh velocity,which makes the mesh move in a Lagrangian manner.Two WENO(Weighted Essentially Non-Oscillatory)reconstructions for the density(the total energy)and the velocity are used to make our scheme achieve the anticipated accuracy.Compared with the general vertexcentered schemes,our scheme with the new approach for the space discretization can simulate some multi-material flows which do not involve large deformations.In addition,our scheme has good robustness,and some numerical examples are presented to demonstrate the anticipated accuracy and the good properties of our scheme.

不同截面肋柱-软尾结构单相流动传热比较

基于任意拉格朗日-欧拉(ALE)方法,利用动网格技术和重叠网格技术数值模拟研究了不同Reynolds数、截面形状以及长径比下的肋柱-软尾结构在通道中的双向流固耦合换热问题。模拟工况为:Reynolds数Re=200,275,351;肋柱截面形状:圆形和方形;长径比k=2,3,4。研究表明:Re=275时圆形截面肋柱-软尾结构的流动换热综合能力优于方形截面结构,k=3时圆形截面肋柱-软尾结构的流动换热综合能力最佳,肋柱周围的局部换热能力方形截面结构好于圆形截面结构;k=3时,随着Reynolds数的增大,不同截面形状的肋柱-软尾结构的流动换热综合能力逐渐升高,并且高Reynolds数、高长径比的圆形截面综合流动换热能力最佳;与不加软尾肋柱结构的流动换热能力进行对比发现,对于圆形截面形状的肋柱结构,增加软尾后综合换热能力增大了19.46%。

阵列起爆爆轰波阵面形状变化及其特征研究

为提升药包的爆炸冲击效能,针对爆轰冲击波传播过程,建立了水下爆炸药包内爆轰波与水中形成的冲击波流固耦合数值模型。引入了阵列角研究考虑攻角效应的单源药包内边缘阵列起爆,并基于爆轰波叠加原理模拟了阵列起爆后爆轰波耦合增强及跨介质传递初始冲击波过程。此外,通过改变阵列角探究了药包内边缘阵列起爆对爆轰波阵面形状、爆轰波压和爆轰波速及水中形成的初始冲击波特征参数的影响。研究表明:两侧对称爆轰波相互作用会导致其传播方向转变为垂直向上,耦合后波压升高,波速小幅度增加;阵列角为20°~45°时,阵列角越小,则爆轰波交汇耦合角度越大,交汇后的耦合爆轰波越早形成类平面波的波阵面。此外,考虑攻角效应的单源药包内边缘阵列起爆对冲击波载荷峰值的增益效果比对波阵面传播速度的增益效果更为显著。本文研究可为水下阵列爆炸的冲击效能优化提供参考。

Simulation of extrusion process of complicated aluminium profile and die trial

The arbitrary Lagrangian-Eulerian（ALE） adaptive remeshing technology and the HyperXtrude software of transient finite element simulations were used on analogue simulation of aluminium extrusion processing.The field distributions of strain rate,stress,temperature and velocity of metal flow were obtained.The results are basically consistent with the experiment,which indicates that this method may successfully predict the defects in the actual extrusion process.

ALE框架下几种不同Godunov型格式的数值比较

给出一种有限体积Godunov型的ALE方法,用于求解多介质大变形的可压缩流动问题.由于方法具有任意的网格移动速度,可在拉氏、欧氏和ALE之间切换,具有较强的适应性.通过数值算例对这3种框架下的数值特性进行了比较研究.同时还研究了几种不同Godunov型格式的数值行为特性,分析比较了它们对激波和接触间断的分辨效果.

基于DEM的任意方向坡度计算方法

坡度是数字地形分析中重要的地形因子,在水文分析、土壤侵蚀模拟、地貌类型划分等地学分析及工程上均有广泛的应用。由于地形的各向异性,坡度在各个方向上并不相同。目前,基于DEM所提取的坡度,均是位于最陡方向上的坡度,而任意给定方向坡度计算是地学分析和工程应用中较为常见的问题。在对目前DEM坡度提取算法研究的基础上,通过数学分析,给出了格网DEM上任意方向格网点的坡度计算模型和计算流程,并给出了相应的计算实例,从而将坡度计算模型从特定方向推向任意方向,使DEM坡度计算更具普适性。

任意形状工作面沉陷预测的概率积分法及其应用

概率积分法是煤矿开采沉陷预测的重要方法,可适用于缓倾斜、倾斜煤层开采地表移动分布曲线符合正态分布规律的地表移动和变形计算。针对原方法对任意形状工作面存在的不足,在原有走向和倾向概率积分函数的基础上,以走向和倾向为基准划分积分函数区域,在该区域内剖分非结构化的三角形单元作为二重积分的基本单位;通过坐标旋转变换的二重积分换元法,转化为以直线为上、下限的二重积分,计算旋转后的影响半径,对新的二重积分采用复合辛普森数值积分公式,积分得出地表任意点在基本三角形开采单元下的积分下沉值;最后,通过叠加计算完成任意点、任意形状工作面概率积分法的沉陷预测计算。算法在GIS系统下实现,实例应用效果良好,可为"三下"采煤提供地表移动预测计算的支持。

一种新型遗传算法及其在土坡任意滑动面确定中的应用

提出了一种便于工程应用的改进模拟任意滑动面策略。针对简单遗传算法容易陷入局部最优值的缺陷,借鉴和声算法产生新解的方式来产生遗传算法中的子代个体。该策略能综合利用整个父体的信息来构建子代。将该方法用于求解两个复杂土坡的最小安全系数,所得结果与已有文献结果比较证明,新型遗传算法能搜索到更小的安全系数,可应用于土坡稳定分析。

Chameleon算法的改进

结合Chameleon算法可以发现高质量的任意形状、大小和密度的自然簇及一趟聚类算法快速高效的特点,研究可以处理混合属性的高效聚类算法.首先简单改进Chameleon算法,使之可以处理含分类属性的数据;进而提出一种两阶段聚类算法.第一阶段使用一趟聚类算法对数据集进行初始划分,第二阶段利用改进的Chameleon算法归并初始划分而得到最终聚类.在真实数据集和人造数据集上的实验结果表明,提出的两阶段聚类算法是有效可行的.

几种任意步距步进相移算法的误差分析与对比

介绍了到目前为止的几种任意步距步进相移算法，并针对相移干涉仪的两种主要误差———移相误差和探测器非线性误差进行了计算仿真，进而比较分析了它们对这这些误差的抑制能力，其结果可为实际应用合理地选择算法提供理论依据。