海量散乱点云out-of-core快速均匀精简算法

海量散乱点云的精简对其显示及交互性操作具有重要的意义。为快速有效地精简海量散乱点云,改进CR树结点分裂算法,将上溢结点子结点包围盒集转化为包围盒的中心点集,利用CR树与数据库SQLite构建主存-辅存分级存储机制,从而实现海量散乱点云的out-of-core管理,计算CR树目标结点层中每个结点所包含点集的均值点,将距离均值点最近的点作为该点集的精简结果,根据目标结点层的不同实现海量散乱点云不同程度的精简。实验证明,基于分级存储机制的快速均匀精简算法能够快速有效地精简海量散乱点云。

Out-of-core clustering of volumetric datasets

In this paper we present a novel method for dividing and clustering large volumetric scalar out-of-core datasets. This work is based on the Ordered Cluster Binary Tree (OCBT) structure created using a top-down or divisive clustering method. The OCBT structure allows fast and efficient sub volume queries to be made in combination with level of detail (LOD) queries of the tree. The initial partitioning of the large out-of-core dataset is done by using non-axis aligned planes calculated using Principal Component Analysis (PCA). A hybrid OCBT structure is also proposed where an in-core cluster binary tree is combined with a large out-of-core file.

On the efficient non-linear solver for hydraulic fracturing and well cementing simulations based on Anderson acceleration

The aim of this study is to create a fast and stable iterative technique for numerical solution of a quasi-linear elliptic pressure equation. We developed a modified version of the Anderson acceleration(AA)algorithm to fixed-point(FP) iteration method. It computes the approximation to the solutions at each iteration based on the history of vectors in extended space, which includes the vector of unknowns, the discrete form of the operator, and the equation's right-hand side. Several constraints are applied to AA algorithm, including a limitation of the time step variation during the iteration process, which allows switching to the base FP iterations to maintain convergence. Compared to the base FP algorithm, the improved version of the AA algorithm enables a reliable and rapid convergence of the iterative solution for the quasi-linear elliptic pressure equation describing the flow of particle-laden yield-stress fluids in a narrow channel during hydraulic fracturing, a key technology for stimulating hydrocarbon-bearing reservoirs. In particular, the proposed AA algorithm allows for faster computations and resolution of unyielding zones in hydraulic fractures that cannot be calculated using the FP algorithm. The quasi-linear elliptic pressure equation under consideration describes various physical processes, such as the displacement of fluids with viscoplastic rheology in a narrow cylindrical annulus during well cementing,the displacement of cross-linked gel in a proppant pack filling hydraulic fractures during the early stage of well production(fracture flowback), and multiphase filtration in a rock formation. We estimate computational complexity of the developed algorithm as compared to Jacobian-based algorithms and show that the performance of the former one is higher in modelling of flows of viscoplastic fluids. We believe that the developed algorithm is a useful numerical tool that can be implemented in commercial simulators to obtain fast and converged solutions to the non-linear problems described above.

面向倾斜摄影测量数据的out-of-core调度机制

对于实时加载无人机倾斜摄影采集的大规模场景数据,场景LOD模型out-of-core调度机制能否快速加载显示用户期望的LOD节点至关重要。本文提出跨越式LOD模型out-of-core调度机制。基本思想是,在限制时间内无法完全加载用户期望的LOD节点时,利用多选择背包问题方法,求解生成替代效果最优的LOD节点集合,并且跨越所有中间层次的节点进行加载显示。实验结果表明,在相同的限制加载时间内,本文算法加载的LOD节点数量为OSG的30%左右,节点数据总加载量为OSG的50%左右,并且本文算法的最终渲染效果比OSG更加精细。

Implementation of a particle-in-cell method for the energy solver in 3D spherical geodynamic modeling

The thermal evolution of the Earth’s interior and its dynamic effects are the focus of Earth sciences.However,the commonly adopted grid-based temperature solver is usually prone to numerical oscillations,especially in the presence of sharp thermal gradients,such as when modeling subducting slabs and rising plumes.This phenomenon prohibits the correct representation of thermal evolution and may cause incorrect implications of geodynamic processes.After examining several approaches for removing these numerical oscillations,we show that the Lagrangian method provides an ideal way to solve this problem.In this study,we propose a particle-in-cell method as a strategy for improving the solution to the energy equation and demonstrate its effectiveness in both one-dimensional and three-dimensional thermal problems,as well as in a global spherical simulation with data assimilation.We have implemented this method in the open-source finite-element code CitcomS,which features a spherical coordinate system,distributed memory parallel computing,and data assimilation algorithms.

TRIP2.0_SOLVER的开发与应用

本文介绍了“亚跨超CFD软件平台”2.0版本流场解算器部分(TRIP2.0_SOLVER)的开发工作和应用情况,重点是结构网格部分的开发工作,包括对接网格、拼接网格、重叠网格等多种网格拓扑结构。TRIP2.0_SOLVER采用有限体积法离散雷诺平均的N-S方程,数值模拟绕流飞行器的空间流场和飞行器的气动特性,本文从软件已经具备的基本功能,用户界面的开发、软件测试和和多种网格拓扑结构的应用等四个方面总结了亚跨超CFD软件平台流场解算器(TRIP2.0_SOLVER)的工作进展情况。

基于ILOG SOLVER的Job-Shop调度算法实现

通过使用约束规划方法对Job-Shop调度问题进行描述和建模,设计用于求解Job-Shop调度问题的禁忌搜索算法,在此基础上基于先进的约束规划系统ILOG对算法进行实现。实践证明基于约束规划将ILOG优化组件应用于对Job-Shop调度问题的求解中,不仅可以大大提高编程效率而且最后结果也有显著提高。

SM—Solver在土木工程专业课教学中的应用

SM-Slover是一个基于结构力学的软件包，既直观又快捷。就其在结构力学、桥梁工程以及钢筋混凝土结构等专业课程教学方面的应用进行了举例说明，表明正确合理地利用SM-Slover，不仅有利于培养学生的专业素养，而且教师也能从中受益，起到了教学相长的效果。

Aspen Open Solver接口技术研究

以Aspen Open Solver接口集中的非线性代数方程组(NLA)部分作为研究对象,在对接口集进行系统地分析之后,利用AspenTech提供的接口代码将分别基于梯度和非基于梯度的四种求解算法嵌入生成solver组件,并实现用Aspen Plus调用该solver组件观察各种算法嵌入的结果。

求解优化问题的混合PSO-Solver算法

融合了粒子群算法(PSO)和Solver加载宏,形成混合PSO-Solver算法进行优化问题的求解。PSO作为全局搜索算法首先给出问题的全局可行解,Solver则是基于梯度信息的局部搜索工具,对粒子群算法得出的解再进行改进,二者互相结合,既加快了全局搜索的速度,又有效地避免了陷入局部最优。算法用VBA语言进行编程,简单且易于实现。通过对无约束优化问题和约束优化问题的求解,以及和标准PSO、其他一些混合算法的比较表明,PSO-Solver算法能够有效地提高求解过程的收敛速度和解的精确性。

约束满足问题求解及ILOG SOLVER系统简介

首先综述求解约束满足问题的基本算法和搜索策略 ,然后介绍ILOGSOLVER求解系统提供的类和函数的基本组成 。

An Irregular Wave Generating Approach Based on naoe-FOAM-SJTU Solver

In this paper,a wave generating approach for long-crest irregular waves in a numerical tank by our in-house solver naoe-FOAM-SJTU is presented.The naoe-FOAM-SJTU solver is developed using an open source tool kit,Open FOAM.Reynolds-averaged Navier？Stokes（RANS） equations are chosen as governing equations and the volume of fluid（VOF） is employed to capture the two phases interface.Incoming wave group is generated by imposing the boundary conditions of the tank inlet.A spectrum based correction procedure is developed to make the measured spectrum approaching to the target spectrum.This procedure can automatically adjust the wave generation signal based on the measured wave elevation by wave height probe in numerical wave tank.After 3 to 4 iterations,the measured spectrum agrees well with the target one.In order to validate this method,several wave spectra are chosen and validated in the numerical wave tank,with comparison between the final measured and target spectra.In order to investigate a practical situation,a modified Wigley hull is placed in the wave tank with incoming irregular waves.The wave-induced heave and pitch motions are treated by Fourier analysis to obtain motion responses,showing good agreements with the measurements.

Lattice Boltzmann Flux Solver:An Efficient Approach for Numerical Simulation of Fluid Flows

A lattice Boltzmann flux solver(LBFS)is presented for simulation of fluid flows.Like the conventional computational fluid dynamics(CFD)solvers,the new solver also applies the finite volume method to discretize the governing differential equations,but the numerical flux at the cell interface is not evaluated by the smooth function approximation or Riemann solvers.Instead,it is evaluated from local solution of lattice Boltzmann equation(LBE)at cell interface.Two versions of LBFS are presented in this paper.One is to locally apply one-dimensional compressible lattice Boltzmann(LB)model along the normal direction to the cell interface for simulation of compressible inviscid flows with shock waves.The other is to locally apply multi-dimensional LB model at cell interface for simulation of incompressible viscous and inviscid flows.The present solver removes the drawbacks of conventional lattice Boltzmann method(LBM)such as limitation to uniform mesh,tie-up of mesh spacing and time interval,limitation to viscous flows.Numerical examples show that the present solver can be well applied to simulate fluid flows with non-uniform mesh and curved boundary.

Efficient Elliptic Solver for the Mild Slope Equation Using BI-CGSTAB Method

The elliptic mild slope equation is used to simulate linear wave propagation over variable sea bed topography with mild slopes. The governing equation is discretized by the finite difference method. Based on the BI-CGSTAB technique, an attractive variant bf BI-Conjugate Gradients (BI-CG) method, the obtained linear algebraic system of equations is solved. Numerical experiments show that the BI-CGSTAB method is efficient for solving the elliptic mild slope equation. The results obtained by the BI-CGSTAB-Based method are much the same as those obtained by other authors with different solution methods, but the convergence rate is much faster than that of other methods.

RANS-VOF Solver for Solitary Wave Run-up on A Circular Cylinder

Simulation of solitary wave run-up on a vertical circular cylinder is carried out in a viscous numerical wave tank developed based on the open source codes Open FOAM. An incompressible two-phase flow solver naoe-FOAM-SJTU is used to solve the Reynolds-Averaged Navier–Stokes（RANS） equations with the SST k ？？ turbulence model. The PISO algorithm is utilized for the pressure-velocity coupling. The air-water interface is captured via Volume of Fluid（VOF） technique. The present numerical model is validated by simulating the solitary wave run-up and reflected against a vertical wall, and solitary wave run-up on a vertical circular cylinder. Comparisons between numerical results and available experimental data show satisfactory agreement. Furthermore, simulations are carried out to study the solitary wave run-up on the cylinder with different incident wave height H and different cylinder radius a. The relationships of the wave run-up height with the incident wave height H, cylinder radius a are analyzed. The evolutions of the scattering free surface and vortex shedding are also presented to give a better understanding of the process of nonlinear wave-cylinder interaction.

Harten-Lax-van Leer-contact(HLLC)approximation Riemann solver with elastic waves for one-dimensional elastic-plastic problems

A Harten-Lax-van Leer-contact （HLLC） approximate Riemann solver is built with elastic waves （HLLCE） for one-dimensional elastic-plastic flows with a hypo- elastic constitutive model and the von Mises＇ yielding criterion. Based on the HLLCE, a third-order cell-centered Lagrangian scheme is built for one-dimensional elastic-plastic problems. A number of numerical experiments are carried out. The numerical results show that the proposed third-order scheme achieves the desired order of accuracy. The third-order scheme is used to the numerical solution of the problems with elastic shock waves and elastic rarefaction waves. The numerical results are compared with a reference solution and the results obtained by other authors. The comparison shows that the pre- sented high-order scheme is convergent, stable, and essentially non-oscillatory. Moreover, the HLLCE is more efficient than the two-rarefaction Riemann solver with elastic waves （TRRSE）

结合SM Solver和MATLAB的结构内力计算分析方法

本文以主梁为结构内力计算为例,应用结构力学求解器(SM Solver)对工程结构进行内力计算分析,将其结果与系数表人工求解方法所得结果的主要数据进行对比分析,表明采用SM Solver计算的结果满足工程计算精度要求。将SM Solver计算结果导入数学分析软件MATLAB绘制出弯矩内力包络图,很好地解决了系数表人工求解方法中无法得到主梁跨中各控制截面最不利内力的问题。将SM Solver和MATLAB结合起来,提供了一种简便易行的结构内力计算和内力包络图绘制方法。

A numerical study for WENO scheme-based on different lattice Boltzmann flux solver for compressible flows

In this paper,the finite difference weighted essentially non-oscillatory (WENO) scheme is incorporated into the recently developed four kinds of lattice Boltzmann flux solver (LBFS) to simulate compressible flows,including inviscid LBFS Ⅰ,viscous LBFS Ⅱ,hybrid LBFS Ⅲ and hybrid LBFS Ⅳ.Hybrid LBFS can automatically realize the switch between inviscid LBFS Ⅰ and viscous LBFS Ⅱ through introducing a switch function.The resultant hybrid WENO-LBFS scheme absorbs the advantages of WENO scheme and hybrid LBFS.We investigate the performance of WENO scheme based on four kinds of LBFS systematically.Numerical results indicate that the devopled hybrid WENO-LBFS scheme has high accuracy,high resolution and no oscillations.It can not only accurately calculate smooth solutions,but also can effectively capture contact discontinuities and strong shock waves.

Kernel polynomial representation for imaginary-time Green's functions in continuous-time quantum Monte Carlo impurity solver

Inspired by the recently proposed Legendre orthogonal polynomial representation for imaginary-time Green s functions G（τ）,we develop an alternate and superior representation for G（τ） and implement it in the hybridization expansion continuous-time quantum Monte Carlo impurity solver.This representation is based on the kernel polynomial method,which introduces some integral kernel functions to filter the numerical fluctuations caused by the explicit truncations of polynomial expansion series and can improve the computational precision significantly.As an illustration of the new representation,we re-examine the imaginary-time Green＇s functions of the single-band Hubbard model in the framework of dynamical mean-field theory.The calculated results suggest that with carefully chosen integral kernel functions,whether the system is metallic or insulating,the Gibbs oscillations found in the previous Legendre orthogonal polynomial representation have been vastly suppressed and remarkable corrections to the measured Green＇s functions have been obtained.

Fast impurity solver for dynamical mean field theory based on second order perturbation around the atomic limit

This paper proposes an impurity solver for the dynamical mean field theory （DMFT） study of the Mott insulators, which is based on the second order perturbation of the hybridization function. After careful benchmarking with quantum Monte Carlo results on the anti-ferromagnetic phase of the Hubbard model, it concludes that this impurity solver can capture the main physical features in the strong coupling regime and can be a very useful tool for the LDA （local density approximation） ＋ DMFT studies of the Mort insulators with long range order.