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.

DIFFERENTIAL-ALGEBRAIC APPROACH TO COUPLED PROBLEMS OF DYNAMIC THERMOELASTICITY

An efficient numerical approach for the general thermomechanical problems was developed and it was tested for a two-dimensional thermoelasticity problem. The main idea of our numerical method is based on the reduction procedure of the original system of PDEs describing coupled thermomechanical behavior to a system of Differential Algebraic Equations （DAEs） where the stress-strain relationships are treated as algebraic equations. The resulting system of DAEs was then solved with a Backward Differentiation Formula （BDF） using a fully implicit algorithm. The described procedure was explained in detail, and its effectiveness was demonstrated on the solution of a transient uncoupled thermoelastic problem, for which an analytical solution is known, as well as on a fully coupled problem in the two-dimensional case.

SOlvaBILITY OF HIGHER INDEX TIME-VARYING LINEAR DIFFERENTIAL-ALGEBRAIC EQUATIONS

Linear differential-algebraic equations (DAEs) with time-varying coefficients A(t)x(1)(t) + B(t)x(t) = q(t), which are tractable with a higher index. are discussed. Their essential properties are investigated. Some equivalent system,,; are given. Using them the paper shows how to state properly initial and boundary conditions for these DAEs. The existence and uniqueness theory of the solution of the initial and boundary value problems for higher index DAEs are proposed.

A Class of Parallel Runge-Kutta Methods for Differential-Algebraic Systems of Index 2

A class of parallel Runge-Kutta Methods for differential-algebraic equations of index 2are constructed for multiprocessor system. This paper gives the order conditions and investigatesthe convergence theory for such methods.

Class of numerical methods for differential-algebraic systems with discontinuous right-hand sides

Numerical methods for Differential-Algebraic systems with discontinuous right-hand sides is discussed. A class of continuous Rosenbrock methods are constructed, and numerical experiments show that the continuous Rosenbrock methods are effective. Applying the methods, a fast and high-precision numerical algorithm is given to deal with typical discontinuous parts, which occur frequently in differential-algebraic systems(DAS).

Stability of Nonlinear Differential-Algebraic Systems Via Additive Identity

The stability analysis for nonlinear differentialalgebraic systems is addressed using tools from classical control theory. Sufficient stability conditions relying on matrix inequalities are established via Lyapunov Direct Method. In addition, a novel interpretation of differential-algebraic systems as feedback interconnection of a purely differential system and an algebraic system allows reducing the stability analysis to a smallgain-like condition. The study of stability properties for constrained mechanical systems, for a class of Lipschitz differential-algebraic systems and for an academic example is used to illustrate the theory.

NONNEGATIVITY OF SOLUTIONS OF NONLINEAR FRACTIONAL DIFFERENTIAL-ALGEBRAIC EQUATIONS

Nonlinear fractional differential-algebraic equations often arise in simulating integrated circuits with superconductors. How to obtain the nonnegative solutions of the equations is an important scientific problem. As far as we known, the nonnegativity of solutions of the nonlinear fractional differential-algebraic equations is still not studied. In this article, we investigate the nonnegativity of solutions of the equations. Firstly, we discuss the existence of nonnegative solutions of the equations, and then we show that the nonnegative solution can be approached by a monotone waveform relaxation sequence provided the initial iteration is chosen properly. The choice of initial iteration is critical and we give a method of finding it. Finally, we present an example to illustrate the efficiency of our method.

ON A REGULARIZATION OF INDEX 2 DIFFERENTIAL-ALGEBRAIC EQUATIONS WITH PROPERLY STATED LEADING TERM

In this article, linear regular index 2 DAEs A（t）[D（t）x（t）]＇ ＋ B（t）x（t） = q（t） are considered. Using a decoupling technique, initial condition and boundary condition are properly formulated. Regular index 1 DAEs are obtained by a regularization method. We study the behavior of the solution of the regularization system via asymptotic expansions. The error analysis between the solutions of the DAEs and its regularization system is given.

Delay and Its Time-Derivative Dependent Stable Criterion for Differential-Algebraic Systems

In this paper, the stable problem for differential-algebraic systems is investigated by a convex op-timization approach. Based on the Lyapunov functional method and the delay partitioning approach, some delay and its time-derivative dependent stable criteria are obtained and formulated in the form of simple linear matrix inequalities (LMIs). The obtained criteria are dependent on the sizes of delay and its time-derivative and are less conservative than those produced by previous approaches.

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.

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.

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.

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

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

Sparse Approximations of the Schur Complement for Parallel Algebraic Hybrid Solvers in 3D

In this paper we study the computational performance of variants of an algebraic additive Schwarz preconditioner for the Schur complement for the solution of large sparse linear systems.In earlier works,the local Schur complements were computed exactly using a sparse direct solver.The robustness of the preconditioner comes at the price of this memory and time intensive computation that is the main bottleneck of the approach for tackling huge problems.In this work we investigate the use of sparse approximation of the dense local Schur complements.These approximations are computed using a partial incomplete LU factorization.Such a numerical calculation is the core of the multi-level incomplete factorization such as the one implemented in pARMS. The numerical and computing performance of the new numerical scheme is illustrated on a set of large 3D convection-diffusion problems;preliminary experiments on linear systems arising from structural mechanics are also reported.

电力系统混合整数线性规划问题的运筹决策关键技术综述与展望

机组组合、检修计划、拓扑运行优化、电力系统规划等电力系统混合整数线性规划(MILP)问题旨在实现电力资源的最佳配置,应用广泛,其精准性与高效性直接影响了电力系统的安全性与经济性。随着“双碳”目标的提出,新型电力系统MILP问题模型复杂度更高、计算效率要求更严格,对当前运筹决策技术提出了更严峻的挑战。然而,现有依赖于国外进口求解器的电力系统运筹决策技术面临“组合爆炸”,且求解器依赖进口面临“卡脖子”困境,亟须实现技术突破。为此,该文系统地梳理了电力系统MILP问题的运筹决策技术,以及近年来通用MILP问题的最新进展,并展望了电力系统MILP问题运筹决策关键技术未来的研究方向,旨在为我国相关研究工作提供参考和思路。

大东湖有压深隧淤积数值模拟研究

基于Godunov型有限体积法建立了能模拟有压管道冲淤的一维水动力-固体物输移-管道冲淤数学模型,并基于该模型模拟预测了大东湖深隧有压管道内的固体物淤积过程,探讨了深隧淤积后的水力冲刷效果。研究结果表明,深隧按现状条件下运行时,管道长期使用基本不会发生淤积,但当遇到极端情况致使固体物浓度达到一定程度时,管道发生明显淤积。当管道最大淤积厚度达到30 cm时,采用水力冲刷将前期淤积物全部冲刷完所需的时间约为16~24 h。

SMTLOC:基于多源频谱的SMT求解器缺陷定位

SMT求解器作为重要的基础软件,其存在的缺陷可能会导致依赖于它的软件功能失效,甚至带来安全事故.然而,修复SMT求解器缺陷是一个十分耗时的任务,因为开发者需要花费大量的时间和精力来理解并找到缺陷的根本原因.虽然已有许多软件缺陷定位方面的研究,但尚未有系统的工作研究如何自动定位SMT求解器缺陷.因此,提出一种基于多源频谱的SMT求解器缺陷定位方法SMTLOC.首先,对于给定的SMT求解器缺陷,SMTLOC提出一种枚举算法,用以对触发该缺陷的公式进行变异,从而生成一组不触发缺陷,但与触发缺陷的公式具有相似执行路径的证人公式.然后,SMTLOC根据证人公式的执行路径以及SMT求解器的源码信息,提出一种融合覆盖频谱和历史频谱的文件可疑度计算方法,从而定位可能存在缺陷的文件.为了验证SMTLOC的有效性,收集60个SMT求解器缺陷.实验结果表明,SMTLOC的缺陷定位效果明显优于传统的频谱缺陷定位方法,SMTLOC可以将46.67%的缺陷定位在TOP-5的文件内,定位效果提升了133.33%.

TRIP2.0_SOLVER的开发与应用

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

基于多样性SAT求解器和新颖性搜索的软件产品线测试

软件产品线测试是一项非常具有挑战性的工作.基于相似性的测试方法通过提升测试集的多样性以达到提高测试覆盖率和缺陷检测率的目的.因其具有良好的可拓展性和较好的测试效果,目前已成为软件产品线测试的重要手段之一.在该测试方法中,如何产生多样化的测试用例和如何维护测试集的多样性是两个关键问题.针对以上问题,提出一种基于多样性可满足性(SAT)求解器和新颖性搜索(novelty search,NS)的软件产品线测试算法.具体地,所提算法同时采用两类多样性SAT求解器产生多样化的测试用例.特别地,为了改善随机局部搜索SAT求解器的多样性,提出一种基于概率向量的通用策略产生候选解.此外,为同时维护测试集的全局和局部多样性,设计并运用两种基于NS算法思想的归档策略.在50个真实软件产品线上的消融和对比实验验证多样性SAT求解器和两种归档策略的有效性,以及所提算法较其他主流算法的优越性.

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

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