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.

Gas kinetic flux solver based finite volume weighted essentially non-oscillatory scheme for inviscid compressible flows

A high-order gas kinetic flux solver(GKFS)is presented for simulating inviscid compressible flows.The weighted essentially non-oscillatory(WENO)scheme on a uniform mesh in the finite volume formulation is combined with the circular function-based GKFS(C-GKFS)to capture more details of the flow fields with fewer grids.Different from most of the current GKFSs,which are constructed based on the Maxwellian distribution function or its equivalent form,the C-GKFS simplifies the Maxwellian distribution function into the circular function,which ensures that the Euler or Navier-Stokes equations can be recovered correctly.This improves the efficiency of the GKFS and reduces its complexity to facilitate the practical application of engineering.Several benchmark cases are simulated,and good agreement can be obtained in comparison with the references,which demonstrates that the high-order C-GKFS can achieve the desired accuracy.

Parallel Iterative FEM Solver with Initial Guess for Frequency Domain Electromagnetic Analysis

The finite element method is a key player in computational electromag-netics for designing RF(Radio Frequency)components such as waveguides.The frequency-domain analysis is fundamental to identify the characteristics of the components.For the conventional frequency-domain electromagnetic analysis using FEM(Finite Element Method),the system matrix is complex-numbered as well as indefinite.The iterative solvers can be faster than the direct solver when the solver convergence is guaranteed and done in a few steps.However,such complex-numbered and indefinite systems are hard to exploit the merit of the iterative solver.It is also hard to benefit from matrix factorization techniques due to varying system matrix parts according to frequency.Overall,it is hard to adopt conventional iterative solvers even though the system matrix is sparse.A new parallel iterative FEM solver for frequency domain analysis is implemented for inhomogeneous waveguide structures in this paper.In this implementation,the previous solution of the iterative solver of Matlab(Matrix Laboratory)employ-ing the preconditioner is used for the initial guess for the next step’s solution process.The overlapped parallel stage using Matlab’s Parallel Computing Toolbox is also proposed to alleviate the cold starting,which ruins the convergence of early steps in each parallel stage.Numerical experiments based on waveguide structures have demonstrated the accuracy and efficiency of the proposed scheme.

Shared Memory Semi-Implicit Solver for Hydrodynamical Instability Processes

The Advection-Diffusion Reaction (ADR) equation appears in many problems in nature. This constitutes a general model that is useful in various scenarios, from porous media to atmospheric processes. Particularly, it is used at the interface between two fluids where different types of instabilities due to surface mobility may appear. Together with the ADR equation, the Darcy-Brinkman model describes the phenomena known as fingering that appear in different contexts. The study of this type of system gains in complexity when the number of chemical species dissolved in both fluids increases. With more solutes, the increasing complexity of this phenomenon generally requires much computational power. To face the need for more computational resources, we build a solver tool based on an Alternating Direction Implicit (ADI) scheme that can be run in Central Processing Unit (CPU) and Graphic Processing Unit (GPU) architectures on any notebook. The implementation is done using the MATLAB platform to compare both versions. It is shown that using the GPU version strongly saves both resources and calculation times.

基于SMSolver的铝合金附着式升降防护平台竖向主框架结构设计校验

附着式升降防护平台是现代建筑施工,特别是高层建筑施工中的重要装备。竖向主框架是附着式升降防护平台中承受和传递平台的竖向和水平荷载的主要结构,结构构件多并且复杂。本文针对一种铝合金附着式升降防护平台的竖向主框架结构进行研究,采用SM Solver软件静力学方法,计算出竖向主框架各构件的内力,再进行构件材料强度校验。

TRIP2.0_SOLVER的开发与应用

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

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

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

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

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

求解优化问题的混合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.

Aspen Open Solver接口技术研究

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

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）

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.

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

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

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.