The goal of this study is to solve the neutron diffusion equation by using a meshless method and evaluate its performance compared to traditional methods. This paper proposes a novel method based on coupling the meshl...The goal of this study is to solve the neutron diffusion equation by using a meshless method and evaluate its performance compared to traditional methods. This paper proposes a novel method based on coupling the meshless local Petrov–Galerkin approach and the moving least squares approximation. This computational procedure consists of two main steps. The first involved applying the moving least squares approximation to construct the shape function based on the problem domain. Then, the obtained shape function was used in the meshless local Petrov–Galerkin method to solve the neutron diffusion equation.Because the meshless method is based on eliminating the mesh-based topologies, the problem domain was represented by a set of arbitrarily distributed nodes. There is no need to use meshes or elements for field variable interpolation. The process of node generation is simply and fully automated, which can save time. As this method is a local weak form, it does not require any background integration cells and all integrations are performed locally over small quadrature domains. To evaluate the proposed method,several problems were considered. The results were compared with those obtained from the analytical solution and a Galerkin finite element method. In addition, the proposed method was used to solve neutronic calculations in thesmall modular reactor. The results were compared with those of the citation code and reference values. The accuracy and precision of the proposed method were acceptable. Additionally, adding the number of nodes and selecting an appropriate weight function improved the performance of the meshless local Petrov–Galerkin method. Therefore, the proposed method represents an accurate and alternative method for calculating core neutronic parameters.展开更多
In this paper, a meshfree boundary integral equation （BIE） method, called the moving Kriging interpolation- based boundary node method （MKIBNM）, is developed for solving two-dimensional potential problems. This st...In this paper, a meshfree boundary integral equation （BIE） method, called the moving Kriging interpolation- based boundary node method （MKIBNM）, is developed for solving two-dimensional potential problems. This study combines the DIE method with the moving Kriging interpolation to present a boundary-type meshfree method, and the corresponding formulae of the MKIBNM are derived. In the present method, the moving Kriging interpolation is applied instead of the traditional moving least-square approximation to overcome Kronecker＇s delta property, then the boundary conditions can be imposed directly and easily. To verify the accuracy and stability of the present formulation, three selected numerical examples are presented to demonstrate the efficiency of MKIBNM numerically.展开更多
In this paper,we apply Fokas unified method to study the initial boundary value(IBV)problems for nonlinear integrable equation with 3×3 Lax pair on the finite interval[0,L].The solution can be expressed by the so...In this paper,we apply Fokas unified method to study the initial boundary value(IBV)problems for nonlinear integrable equation with 3×3 Lax pair on the finite interval[0,L].The solution can be expressed by the solution of a 3×3 Riemann-Hilbert(RH)problem.The relevant jump matrices are written in terms of matrix-value spectral functions s(k),S(k),S_(l)(k),which are determined by initial data at t=0,boundary values at x=0 and boundary values at x=L,respectively.What's more,since the eigenvalues of 3×3 coefficient matrix of k spectral parameter in Lax pair are three different values,search for the path of analytic functions in RH problem becomes a very interesting thing.展开更多
间接边界积分方程法IBIEM(indirect boundary integral equation method)求解波动问题时控制方程基本解构造依赖经验判断和试算,导致宽频散射求解不够稳定。本文通过粒子群优化-人工神经网络建立IBIEM控制方程基本解构造模型,以数据驱...间接边界积分方程法IBIEM(indirect boundary integral equation method)求解波动问题时控制方程基本解构造依赖经验判断和试算,导致宽频散射求解不够稳定。本文通过粒子群优化-人工神经网络建立IBIEM控制方程基本解构造模型,以数据驱动代替经验判断,处理基本解构造过程中的不确定性。以二维峡谷对平面SH波散射IBIEM模拟为例验证所建模型的可靠性。结果表明,所建IBIEM控制方程基本解构造模型可对虚拟波源位置和数量的最优设置进行有效预测,兼顾计算效率和精度,大幅提高IBIEM求解波动问题时的稳定性和高效性;虚拟波源位置和数量最优设置方案受入射波频率和场地几何条件影响显著,且表现出非单调变化特征,依据经验设置基本解可靠性较差,以数据驱动的预测模型具有明显优势。本文所建方法可为IBIEM求解其他类型场地地震波动问题提供参考。展开更多
Traditional geoacoustic inversions are generally solved by matched-field processing in combination with metaheuristic global searching algorithms which usually need massive computations. This paper proposes a new phys...Traditional geoacoustic inversions are generally solved by matched-field processing in combination with metaheuristic global searching algorithms which usually need massive computations. This paper proposes a new physical framework for geoacoustic retrievals. A parabolic approximation of wave equation with non-local boundary condition is used as the forward propagation model. The expressions of the corresponding tangent linear model and the adjoint operator are derived, respectively, by variational method. The analytical expressions for the gradient of the cost function with respect to the control variables can be formulated by the adjoint operator, which in turn can be used for optimization by the gradient-based method.展开更多
针对声学边界元法中解的非唯一性和奇异积分问题,基于组合亥姆霍兹积分方程公式(combined helmholtz integral equation formulation,CHIEF)法思想,将常规边界元方程和等效源方程进行联立,并利用两者方程系数矩阵间的耦合等价关系,间接...针对声学边界元法中解的非唯一性和奇异积分问题,基于组合亥姆霍兹积分方程公式(combined helmholtz integral equation formulation,CHIEF)法思想,将常规边界元方程和等效源方程进行联立,并利用两者方程系数矩阵间的耦合等价关系,间接替换计算常规边界元法中的奇异系数矩阵,进而提出一种具有全频域唯一解、高计算精度和高稳定性的耦合CHIEF法。该方法将等效源方程作为补充方程,不仅解决了传统CHIEF法内点补充方程失效的问题,而且矩阵的间接替换计算避免了直接计算奇异积分,显著提高了计算效率和精度。通过声辐射和声散射的典型算例对比了所提方法、常规边界元法、常规Burton-Miller法和等效源法的计算效果。结果表明,所提方法不仅在全波数域内均能获得唯一解,且其计算精度和效率均优于常规边界元法和常规Burton-Miller方法,其系数矩阵条件数远低于等效源法。展开更多
无网格局部彼得洛夫-伽辽金(meshless local Petrov-Galerkin,MLPG)法是一种具有代表性的无网格方法,在计算力学领域得到广泛应用.然而,这种方法在边界上需执行积分运算,通常很难处理不规则求解域问题.为了克服MLPG法的这种局限性,提出...无网格局部彼得洛夫-伽辽金(meshless local Petrov-Galerkin,MLPG)法是一种具有代表性的无网格方法,在计算力学领域得到广泛应用.然而,这种方法在边界上需执行积分运算,通常很难处理不规则求解域问题.为了克服MLPG法的这种局限性,提出了无网格局部强弱(meshless local strong-weak,MLSW)法.MLSW法采用MLPG法离散内部求解域,采用无网格介点(meshless intervention-point,MIP)法施加自然边界条件,并采用配点法施加本质边界条件,避免执行边界积分运算,可适用于求解各类复杂的不规则域问题.从理论上讲,这种结合式方法,既保持了MLPG法稳定而精确计算的优势,同时兼备配点型方法在处理复杂结构问题时简洁而灵活的优势,实现了弱式法和强式法的优势互补.此外,MLSW法采用移动最小二乘核(moving least squares core,MLSc)近似法来构造形函数,是对传统移动最小二乘(moving least squares,MLS)近似法的一种改进.MLSc使用核基函数代替通常的基函数,有利于数值求解的精确性和稳定性,而且其导数近似计算变得更为简单.数值算例结果初步表明:这种新方法实施简单,求解稳定、精确,表现出适合工程运用的潜力.展开更多
文摘The goal of this study is to solve the neutron diffusion equation by using a meshless method and evaluate its performance compared to traditional methods. This paper proposes a novel method based on coupling the meshless local Petrov–Galerkin approach and the moving least squares approximation. This computational procedure consists of two main steps. The first involved applying the moving least squares approximation to construct the shape function based on the problem domain. Then, the obtained shape function was used in the meshless local Petrov–Galerkin method to solve the neutron diffusion equation.Because the meshless method is based on eliminating the mesh-based topologies, the problem domain was represented by a set of arbitrarily distributed nodes. There is no need to use meshes or elements for field variable interpolation. The process of node generation is simply and fully automated, which can save time. As this method is a local weak form, it does not require any background integration cells and all integrations are performed locally over small quadrature domains. To evaluate the proposed method,several problems were considered. The results were compared with those obtained from the analytical solution and a Galerkin finite element method. In addition, the proposed method was used to solve neutronic calculations in thesmall modular reactor. The results were compared with those of the citation code and reference values. The accuracy and precision of the proposed method were acceptable. Additionally, adding the number of nodes and selecting an appropriate weight function improved the performance of the meshless local Petrov–Galerkin method. Therefore, the proposed method represents an accurate and alternative method for calculating core neutronic parameters.
基金Project supported by the Young Scientists Fund of the National Natural Science Foundation of China(Grant No.10902076)the Natural Science Foundation of Shanxi Province of China(Grant No.2007011009)+1 种基金the Scientific Research and Development Program of the Shanxi Higher Education Institutions(Grant No.20091131)the Doctoral Startup Foundation of Taiyuan University of Science and Technology(Grant No.200708)
文摘In this paper, a meshfree boundary integral equation （BIE） method, called the moving Kriging interpolation- based boundary node method （MKIBNM）, is developed for solving two-dimensional potential problems. This study combines the DIE method with the moving Kriging interpolation to present a boundary-type meshfree method, and the corresponding formulae of the MKIBNM are derived. In the present method, the moving Kriging interpolation is applied instead of the traditional moving least-square approximation to overcome Kronecker＇s delta property, then the boundary conditions can be imposed directly and easily. To verify the accuracy and stability of the present formulation, three selected numerical examples are presented to demonstrate the efficiency of MKIBNM numerically.
基金supported by the National Natural Science Foundation of China(11901167,11971313 and 51879045)Key scientific research projects of higher education institutions in Henan,China(18B110008).
文摘In this paper,we apply Fokas unified method to study the initial boundary value(IBV)problems for nonlinear integrable equation with 3×3 Lax pair on the finite interval[0,L].The solution can be expressed by the solution of a 3×3 Riemann-Hilbert(RH)problem.The relevant jump matrices are written in terms of matrix-value spectral functions s(k),S(k),S_(l)(k),which are determined by initial data at t=0,boundary values at x=0 and boundary values at x=L,respectively.What's more,since the eigenvalues of 3×3 coefficient matrix of k spectral parameter in Lax pair are three different values,search for the path of analytic functions in RH problem becomes a very interesting thing.
文摘间接边界积分方程法IBIEM(indirect boundary integral equation method)求解波动问题时控制方程基本解构造依赖经验判断和试算,导致宽频散射求解不够稳定。本文通过粒子群优化-人工神经网络建立IBIEM控制方程基本解构造模型,以数据驱动代替经验判断,处理基本解构造过程中的不确定性。以二维峡谷对平面SH波散射IBIEM模拟为例验证所建模型的可靠性。结果表明,所建IBIEM控制方程基本解构造模型可对虚拟波源位置和数量的最优设置进行有效预测,兼顾计算效率和精度,大幅提高IBIEM求解波动问题时的稳定性和高效性;虚拟波源位置和数量最优设置方案受入射波频率和场地几何条件影响显著,且表现出非单调变化特征,依据经验设置基本解可靠性较差,以数据驱动的预测模型具有明显优势。本文所建方法可为IBIEM求解其他类型场地地震波动问题提供参考。
基金Project supported by the Foundation of Key Laboratory of Marine Intelligent Equipment and System of Ministry of Education,China(Grant No.SJTUMIES1908)the National Natural Science Foundation of China(Grant No.41775027)
文摘Traditional geoacoustic inversions are generally solved by matched-field processing in combination with metaheuristic global searching algorithms which usually need massive computations. This paper proposes a new physical framework for geoacoustic retrievals. A parabolic approximation of wave equation with non-local boundary condition is used as the forward propagation model. The expressions of the corresponding tangent linear model and the adjoint operator are derived, respectively, by variational method. The analytical expressions for the gradient of the cost function with respect to the control variables can be formulated by the adjoint operator, which in turn can be used for optimization by the gradient-based method.
文摘针对声学边界元法中解的非唯一性和奇异积分问题,基于组合亥姆霍兹积分方程公式(combined helmholtz integral equation formulation,CHIEF)法思想,将常规边界元方程和等效源方程进行联立,并利用两者方程系数矩阵间的耦合等价关系,间接替换计算常规边界元法中的奇异系数矩阵,进而提出一种具有全频域唯一解、高计算精度和高稳定性的耦合CHIEF法。该方法将等效源方程作为补充方程,不仅解决了传统CHIEF法内点补充方程失效的问题,而且矩阵的间接替换计算避免了直接计算奇异积分,显著提高了计算效率和精度。通过声辐射和声散射的典型算例对比了所提方法、常规边界元法、常规Burton-Miller法和等效源法的计算效果。结果表明,所提方法不仅在全波数域内均能获得唯一解,且其计算精度和效率均优于常规边界元法和常规Burton-Miller方法,其系数矩阵条件数远低于等效源法。
文摘无网格局部彼得洛夫-伽辽金(meshless local Petrov-Galerkin,MLPG)法是一种具有代表性的无网格方法,在计算力学领域得到广泛应用.然而,这种方法在边界上需执行积分运算,通常很难处理不规则求解域问题.为了克服MLPG法的这种局限性,提出了无网格局部强弱(meshless local strong-weak,MLSW)法.MLSW法采用MLPG法离散内部求解域,采用无网格介点(meshless intervention-point,MIP)法施加自然边界条件,并采用配点法施加本质边界条件,避免执行边界积分运算,可适用于求解各类复杂的不规则域问题.从理论上讲,这种结合式方法,既保持了MLPG法稳定而精确计算的优势,同时兼备配点型方法在处理复杂结构问题时简洁而灵活的优势,实现了弱式法和强式法的优势互补.此外,MLSW法采用移动最小二乘核(moving least squares core,MLSc)近似法来构造形函数,是对传统移动最小二乘(moving least squares,MLS)近似法的一种改进.MLSc使用核基函数代替通常的基函数,有利于数值求解的精确性和稳定性,而且其导数近似计算变得更为简单.数值算例结果初步表明:这种新方法实施简单,求解稳定、精确,表现出适合工程运用的潜力.
