Research and Application of a Multi-Field Co-Simulation Data Extraction Method Based on Adaptive Infinitesimal Element

Regarding the spatial profile extraction method of a multi-field co-simulation dataset,different extraction directions,locations,and numbers of profileswill greatly affect the representativeness and integrity of data.In this study,a multi-field co-simulation data extractionmethod based on adaptive infinitesimal elements is proposed.Themultifield co-simulation dataset based on related infinitesimal elements is constructed,and the candidate directions of data profile extraction undergo dimension reduction by principal component analysis to determine the direction of data extraction.Based on the fireworks algorithm,the data profile with optimal representativeness is searched adaptively in different data extraction intervals to realize the adaptive calculation of data extraction micro-step length.The multi-field co-simulation data extraction process based on adaptive microelement is established and applied to the data extraction process of the multi-field co-simulation dataset of the sintering furnace.Compared with traditional data extraction methods for multi-field co-simulation,the approximate model constructed by the data extracted from the proposed method has higher construction efficiency.Meanwhile,the relative maximum absolute error,root mean square error,and coefficient of determination of the approximationmodel are better than those of the approximation model constructed by the data extracted from traditional methods,indicating higher accuracy,it is verified that the proposed method demonstrates sound adaptability and extraction efficiency.

Multi-field dynamic modeling and numerical simulation of aluminum alloy resistance spot welding

In order to explore the influence of welding parameters and to investigate the Al alloy (AA) nugget formation process, a comprehensive model involving electrical-thermal-mechanical and metallurgical analysis was established to numerically display the resistance spot welding (RSW) process within multiple fields and understand the AA-RSW physics. A multi-disciplinary finite element method (FEM) framework and a empirical sub-model were built to analyze the affecting factors on weld nugget and the underlying nature of welding physics with dynamic simulation procedure. Specifically, a counter-intuitive phenomenon of the resistance time-variation caused by the transient inverse virtual variation (TIVV) effect was highlighted and analyzed on the basis of welding current and temperature distribution simulation. The empirical model describing the TIVV phenomenon was used for modifying the dynamic resistance simulation during the AA spot welding process. The numerical and experimental results show that the proposed multi-field FEM model agrees with the measured AA welding feature, and the modified dynamic resistance model captures the physics of nugget growth and the electrical-thermal behavior under varying welding current and fluctuating heat input.

Experimental research and numerical simulation of the multi-field performance of cemented paste backfill:Review and future perspectives

Cemented paste backfill(CPB)technology is a green mining method used to control underground goaves and tailings ponds.The curing process of CPB in the stope is the product of a thermo-hydro-mechanical-chemical multi-field performance interaction.At present,research on the multi-field performance of CPB mainly includes indoor similar simulation experiments,in-situ multi-field performance monitoring experiments,multi-field performance coupling model construction of CPB,and numerical simulation of the multi-field performance of CPB.Because it is hard to study the in-situ multi-field performance of CPB in the real stope,most current research on in-situ multi-field performance adopts the numerical simulation method.By simulating the conditions of CPB in the real stope(e.g.,maintenance environment,stope geometry,drainage conditions,and barricade and backfilling rates),the multi-field performance of CPB is further studied.This paper summarizes the mathematical models employed in the numerical simulation and lists the engineering application cases of numerical simulation in the in-situ multi-field performance of CPB.Finally,it proposes that the multi-field performance of CPB needs to strengthen the theoretical study of multi-field performance,form the strength design criterion based on the multi-field performance of CPB,perform a full-range numerical simulation of the multi-field performance of CPB,develop a pre-warning technology for the CPB safety of CPB,develop automatic and wireless sensors for the multi-field performance monitoring of CPB,and realize the application and popularization of CPB monitoring technology.

On the Transient Numerical Simulation of Solid Rocket Motor by Coupling Quasi One⁃Dimension Internal Flow with Three⁃Dimension Propellant Grain Burnback

The quasi one⁃dimension compressible flowfield coupled to the three⁃dimension propellant grain regression solved by the level⁃set method was used to simulate the transient internal ballistics of solid rocket motor.One⁃dimension flowfield instead of three⁃dimension can save computational cost on the premise of calculation accuracy because the radial and azimuthal variations parameters have little contribution to the internal flowfield.The grain regression in real⁃time could provide accurate geometrical information for simulation.A combination of flowfluid solver and grain regression can reappear in a relatively real internal ballistic flowfield,so it is good for further studying the instability of solid rocket motor.For level⁃set equations,the total variation⁃diminishing second⁃order Runge⁃Kutta method for temporal derivatives and a fifth⁃order weighted⁃essentially⁃non⁃oscillatory scheme for spatial derivatives were used.The total variation⁃diminishing MacCormack method was used to discrete the Euler equations in flowfield solver.Two modules of this code were tested in this study:one is the burning rate module and the other is the nozzle erosion module.Results show that the burning rate influenced the solid rocket motor efficiency,and the velocity profile in the chamber was affected by the nozzle shape,and the nozzle erosion could influence the head⁃end pressure spike.

Multi-field Coupling Simulation of Impact of Temperature and Density of Heat Injection Well on Carbon Budget during Hydrate Exploitation in Qilian Mountain Permafrost Region

Permafrost regions of Qilian Mountains in China are rich in gas hydrate resources.Once greenhouse gases in deep frozen layer are released into the atmosphere during hydrate mining,a series of negative consequences occur.This study aims to evaluate the impact of hydrate thermal exploitation on regional permafrost and carbon budgets based on a multi-physical field coupling simulation.The results indicate that the permeability of the frozen soil is anisotropic,and the low permeability frozen layer can seal the methane gas in the natural state.Heat injection mining of hydrates causes the continuous melting of permafrost and the escape of methane gas,which transforms the regional permafrost from a carbon sink to a carbon source.A higher injection temperature concentrates the heat and causes uneven melting of the upper frozen layer,which provides a dominant channel for methane gas and results in increased methane emissions.However,dense heat injection wells cause more uniform melting of the lower permafrost layer,and the melting zone does not extend to the upper low permeability formation,which cannot provide advantageous channels for methane gas.Therefore,a reasonable and dense number of heat injection wells can reduce the risk of greenhouse gas emissions during hydrate exploitation.

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.

基于GLC多项式谱元法的大地电磁场二维正演模拟

为了提高大地电磁场数值模拟的精度和效率,提出基于Gauss-Lobatto-Chebyshev(GLC)基函数谱元法的大地电磁场正演模拟方法。该方法首先从麦克斯韦方程组出发,推导了二维大地电磁场边值问题;然后基于Galerkin加权余量法将微分形式的边值问题转换成积分弱形式;最后采用GLC正交多项式插值基函数对全局问题进行离散化,利用Pardiso求解器求解大型稀疏线性方程组得到大地电磁场,实现了二维大地电磁场数值模拟。为了提高数值模拟效率,文中算法采用变密度规则网格剖分技术,即在电性复杂区域使用细网格,在电性均匀区域使用粗网格,采用OpenMP编程模式实现了多个频点的并行计算,以此达到缩短计算时间的目的。通过一维层状介质模型数值模拟结果验证了文中算法的正确性和精度,相比于Gauss-Lobatto-Legendre(GLL)多项式谱元法的数值模拟结果,GLC多项式谱元法的模拟结果精度更高。针对国际标准模型COMMEMI 2D-1和带地形模型,进行了基于GLC多项式谱元法、有限差分法和基于三角网格有限元法的正演模拟,三者数值结果对比表明,GLC多项式谱元法计算精度更高、网格依赖性更低。

基于黎曼求解器的SPH算法及其应用

基于传统光滑粒子流体动力学方法(SPH)进行数值模拟时,在边界处存在压力和速度震荡,甚至导致求解失败。为准确求解SPH方程式中的粒子的物理粘性,减小数值耗散,将黎曼求解器纳入SPH格式中,通过重整化过程对连续性方程中的梯度算子进行了校正。通过一系列经典算例的数值模拟结果表明,基于黎曼求解器的SPH算法能够准确的模拟液体的流动,且具有良好的准确性和鲁棒性。

DSMC: Fast direct simulation Monte Carlo solver for the Boltzmann equation by Multi-Chain Markov Chain and multicore programming

Direct Simulation Monte Carlo(DSMC)solves the Boltzmann equation with large Knudsen number.The Boltzmann equation generally consists of three terms:the force term,the diffusion term and the collision term.While the first two terms of the Boltzmann equation can be discretized by numerical methods such as the finite volume method,the third term can be approximated by DSMC,and DSMC simulates the physical behaviors of gas molecules.However,because of the low sampling efficiency of Monte Carlo Simulation in DSMC,this part usually occupies large portion of computational costs to solve the Boltzmann equation.In this paper,by Markov Chain Monte Carlo(MCMC)and multicore programming,we develop Direct Simulation Multi-Chain Markov Chain Monte Carlo(DSMC3):a fast solver to calculate the numerical solution for the Boltzmann equation.Computational results show that DSMC3 is significantly faster than the conventional method DSMC.

基于EES的太阳能光热发电系统性能模拟与优化

随着全球能源结构的持续转型和对可持续能源技术的不断追求,太阳能光热发电作为一种高效利用太阳能的方式受到了广泛关注。基于工程方程求解软件(Engineering Equation Solver,EES),构建高契合度的太阳能光热发电系统数学模型,对系统中的热力过程、能量转换效率以及关键影响因素进行模拟分析,进而提出一系列具有针对性的优化措施。实践证明,基于EES的太阳能光热发电系统的性能模拟与优化有助于提高系统光热发电效率,提升系统的经济性和环境友好性,降低能源损耗,具有广阔的应用前景。

优化和并行一个油藏数值模拟软件中的解法器

在当前共享存储对称多处理 ( SMP)并行机上 ,基于指导语句的并行程序设计模式 ,讨论了多功能油藏数值模拟软件中求解超过百万节点规模的解法器 ( MFS)的并行和优化技术。首先 ,结合当前微处理器的高性能特征 ,为了提高 Cache命中率 ,改进了数据和循环结构 ,并组织了 MFS的性能优化 ,在 R50 0 0上获得了 2 0 %的性能提高 ,并消除了并行化将可能引入的 Cache一致性冲突 ;然后 ,基于循环合并、区域分解和大粒度流水线并行技术 ,实现了 MFS的并行化 ;最后 ,在 POW-ER CHALL ENGE R80 0 0的 6台处理机和 R1 0 0 0 0的 8台处理机上 ,对三维三相 50万和 1 0 0万节点规模问题 ,分别组织了数值实验 ,并取得了超过 60 %的并行效率。

航空公司机组重调度问题建模和算法研究

不正常航班机组重调度问题是航空公司航班重调度决策研究的重要环节,针对机组重调度问题建模和算法进行研究。采用约束规划建模,基于ILOG Solver优化组件提供的框架,将模拟退火算法嵌入到ILOG Solver框架中对模型进行求解,并给出邻域生成搜索策略和移动的选择策略,使邻域搜索效率获得改善。案例计算结果显示,新的求解方法在CPU时间消耗和目标函数优化两个方面比直接采用ILOG Solver优化组件获得较大改善。

空间碎片超高速碰撞的数值方法研究

应用非线性动力学数值仿真软件AUTODYN-2D对空间碎片的超高速碰撞问题进行模拟,研究不同数值计算方法和材料模型的适用范围,空间碎片形状对碎片云的影响,冲击相变等问题。结果表明,在研究超高速碰撞中普遍采用的SPH方法可以很好地模拟碎片云的情况,但是难以对冲击相变问题进行模拟,而采用侵蚀模型的拉格朗日方法可较好地反映冲击相变现象。Tilloston状态方程适用于描述冲击相变问题,包括汽化问题。Johnson-Cook强度模型用来对应变硬化、应变率硬化以及熔化等问题进行模拟。因此,在当速度足以使材料发生相变时,Tilloston状态方程和Johnson-Cook强度模型是比较好的选择。

单级感应线圈炮工作过程的动态仿真

为了更好地对多级同步感应线圈炮进行设计和试验研究,为储能脉冲电容器驱动的单级感应线圈炮建立了描述其动态工作过程的数学模型(包括电路方程和电枢运动方程),并利用Ansoft有限元分析软件和2-D瞬态求解器进行了仿真研究。这种方法可以避免驱动线圈自感、电枢自感以及驱动线圈与电枢之间互感等参数计算的难题。仿真结果给出了驱动电流波形、电枢受到的电磁力、位移以及速度等参数随时间变化的波形。这对于预测单级感应线圈炮的加速性能具有十分重要的意义。仿真结果表明:感应线圈炮的主要加速作用发生在驱动电流的第1个半周期内(主要加速时间小于半周期),驱动电流有时对电枢起制动作用,制动作用范围与感应线圈炮的初始工作状况有关。研究结果表明:该仿真方法可用于对单级感应线圈炮的最佳初始位置进行优化,也可用于多级同步感应线圈炮工作过程的动态仿真。

卫星姿控发动机混合物羽流场分区耦合计算研究

研究求解喷管内流场N-S方程数值计算方法,发展基于N-S方程物面边界滑移流理论计算技术。提出求解羽流核心区轴对称DSMC模拟方法与远场三维DSMC仿真方案,发展多组元混合物羽流DSMC仿真方法。研究求解卫星姿控发动机内外近场、远场、倒流区和物面相互作用影响区多流域流场分区耦合计算技术,建立了一套用于求解混合物燃气羽流及对太阳电池帆板与卫星体表面撞击污染影响数值模拟方法。通过对分别安装于某在轨卫星不同位置两个典型姿控发动机燃气五组元混合物羽流计算研究及相关结果对比分析,证实本文数值方法可靠性。

全过程动态仿真中大型线性方程组的分块求解算法

电力系统全过程动态仿真能够将机电暂态、中期和长期动态过程有机地统一起来进行数字仿真,仿真过程中需要多次求解大型稀疏线性方程组。该方程组由电力系统设备模型的微分—代数方程式差分后的代数方程和输电网络模型的代数方程形成,其快速求解算法是电力系统全过程动态仿真的难点之一。文中提出一种利用仿真中矩阵结构特点的分块快速直接求解算法,并开发实现了大型电力系统线性方程组稀疏求解器(ESS)。该算法首先将稀疏矩阵分为4个分块矩阵,然后将其中规模最大的对角块进一步细分为多个更小的对角分块矩阵,并利用部分小分块具有相同结构的特点进行矩阵LU符号分解和数值分解,最后根据分块矩阵进行前代和回代求解计算。与现有其他求解器进行的算例对比表明,ESS具有较为明显的整体求解速度优势,特别是在矩阵LU分解方面。

航行器水中运动数学模型及数值仿真

根据水下垂直发射航行器水中运动的特点,利用质心运动定理和绕质心的动量矩定理,推导建立航行器水中运动的运动学模型和动力学模型。采用Matlab软件中的ODE解算器,对建立的数学模型进行数值仿真,得到航行器水中运动的速度和姿态角。

飞机安全性设计的外挂物投放数值方法

研究了飞机外挂物投放问题的数值模拟方法,并应用于安全性设计.通过采用运动嵌套网格方法,结合三维Euler方程对飞机外挂物投放过程进行了数值模拟,并针对该投放模拟方法提出了一套合理可行的并行计算方案.由不同弹射力作用方式对于外挂物投放安全性设计的影响,证实所提出方法对投放过程的数值模拟是有效的.