1 Introduction
The numerical tools of Partial Differential Equations (PDE) play a key role in the accuracy and acceleration of the solution process in engineering problems. In the last few years, significant progress has been achieved in this field due to the development of computer technology. Despite this progress, numerical methods are often generalized to solve differential equations in a heuristic manner [1]. Numerical methods consist of two types based on the discretization method. The first employs a mesh (mesh-based method) and the second uses a local approximation by polynomials (meshless method).
Mesh-based numerical methods are divided into three classical groups for solving PDEs: Finite Difference Methods (FDM), Finite Volume Methods (FVM), and Finite Element Methods (FEM) [2]. These methods have been used to solve engineering problems for more than 60 years and cover a wide scope of problems such as material science, aerospace engineering, computational physics, safety, reliability analysis, solid mechanic, fluid mechanical interactions, molecular interaction, etc. Mesh-based methods are powerful numerical tools for solving problems and evaluating, predicting, and simulating the behavior of systems. However, mesh-based methods have some limitations in that they cannot provide fully acceptable results. Moreover, despite the great progress in mesh generation, creating an accurate mesh remains difficult. In addition, mesh-based methods are typically time consuming [1].
Therefore, there is an interest in developing alternative methods that remove or reduce the need for meshing. A review of these methods shows that meshless methods have achieved considerable success in recent decades and have attracted substantial interest from researchers because of their applications for solving computational problems [3–9]. The most outstanding of these methods are introduced below.
The advent of the meshless concept dates back to the model of astrophysical phenomena without boundaries (Smooth Particle Hydrodynamics (SPH)); for example, exploding stars and dust clouds. This idea was proposed as a theory for non-spherical stars by Gingold and Monaghan [10]. The SPH scheme has been successfully applied to many problems such as heat conduction, explosion phenomena, and free surface [11, 12]. Then, Nayroles et al. (1992) proposed the diffuse element method that has been employed by many researchers. Belytschko et al. (2004) introduced an extended version of Nayroles's method, which included a series of improvements over the diffuse element method formulation (the element-less Galerkin methods) [14]. Liu et al. (1995, 1993) introduced a new method for correcting the lack of consistency in the SPH scheme, known as the reproducing kernel particle method.
Oñate et al. (1996) studied the finite point manner to model elasticity and plate bending and fluid flow problems using the collocation point technique, least squares approximation, and weighted least squares approximation. Since then, research has tended to use the radial basis function approximation technique [18–21]. Several studies on radial basis functions have been applied to solve physical problems such as transport phenomena [22], heat conduction [23], neutron diffusion [24, 25], and analysis of Kirchhoff plates [26]. De and Bathe (2000) proposed the finite spheres approach for constructing the approximation function. Moradipour and Yousefi (2016) used a meshless kernel-based method to solve the Black–Scholes variational inequality of American options. In addition to the analysis of two-dimensional elasticity problems, Ebrahimnejad et al. (2015) presented three schemes of the 2D meshless finite volume method. One of the main methods of meshless schemes is the Meshless Local Petrov-Galerkin (MLPG) method proposed by [30]. The MLPG method has frequently been used to successfully solve many engineering problems [31–35].
This method is based on the idea of the local weak form, which eliminates the requirement of the background cell. In addition, this method provides numerical integration in a meshless sense to simplify the integrand of the weak form [4, 6, 30, 36–38]. In this global weak form method, a rational basis is provided to construct meshless methods with a larger degree of flexibility [36, 39–44]. The MLPG method is an effective method for solving many problems, which is based on a local weak form and the Moving Least Squares Approximation (MLSA) to approximate the shape functions [45, 46]. As this method does not need a "finite element mesh", it is a truly meshless method [2]. All integrals in the formulation can be easily evaluated over regularly shaped domains and their boundaries [30]. The shape functions must be effectively constructed to achieve a desirable order of continuity.
Accordingly, the capability and flexibility of meshless methods provide the motivation for using this approach to solve the neutron diffusion equation in this study.
2 Methodology
In this research, the MLPG method was used to solve the neutron diffusion equation in two-dimensional geometry. First, the solution was approximated by using a local approximation in the problem domain. This approximation scheme was used to formulate the residual of the governing equations on the problem and find the shape function. The obtained shape function was employed in the MLPG algorithm to solve the problem. Then, the problem was modeled by the analytical solutions and the Galerkin finite element method (GFEM) method to evaluate the proposed method. In the following subsections, all steps required to derive the shape function and the weak form of the problem equation are explained. A schematic representation of the local support domains used to create the shape function is demonstrated in Fig. 1. The domain representation of the assumed problem of the mesh-based and meshless approaches is shown in Fig. 2.
-201811/1001-8042-29-11-018/alternativeImage/1001-8042-29-11-018-F001.jpg)
-201811/1001-8042-29-11-018/alternativeImage/1001-8042-29-11-018-F002.jpg)
2.1 Determining the shape function by MLSA
Selection of an approximation function has an important role for reaching a sufficient number of equations to determine the unknowns [47, 48]. Some of the most important functions consist of Shepard functions [40, 49], least squares [14, 50], local radial point interpolation [51], and local point interpolation [52]. One of the best schemes for interpolating data with an appropriate accuracy is the moving least squares method, which has more compatibility with MLPG method [53]. The moving least squares approximation was created by mathematic researchers for data fitting and surface construction [54]. The steps for determination of the shape function are as follows:
1) First, to demonstrate the trial function, the meshless approach should use local interpolation. Therefore, a function of the field unknown variables, U(X), is defined in a subdomain, Ωs, which indicates the support domain of the MLS approximation for U(X). It should be mentioned that the support domain is located in the problem domain Ωp. X=(x, y), which indicates a two-dimensional problem. To find an approximation distribution of function U in the support domain, the local moving squares approximation of function U(X) for all X in the support domain should be defined as follows [42]:
PT(X) is a complete monomial basis of order m. m is the number of monomial terms, here the linear basis m=3.
a(X) is a coefficient vector as follow:
The approximated values of the field function are:
2) In the second step, the weighted discrete L2 norm is defined as follows:
Ui is the nodal parameter of the field variable at node i. n indicates the number of neighborhood nodes of node i. The weight function
3) In the third step, for determining the coefficients matrix, a(X), J should be minimized with respect to a(X) [42]:
This equation leads to a linear relationship (7) between a(X) and
Solving Eq.(7) for a(X) gives [54]:
where Us is a vector for choosing the field function nodal parameters for interior nodes of the support domain as follows:
The matrices A(X) and B(X) are given as [54, 55]:
where
Hence, for this two-dimensional problem, A(X) is a symmetric matrix according to the linear basis (m=3) that can be explicitly written as:
-201811/1001-8042-29-11-018/media/1001-8042-29-11-018-M001.jpg)
少公式(12)
4) The fourth step is the derivation of the shape function. If the above equations are substituted into Eq. (4), the form of an interpolation function that is MLSA of u(X) at X is written as:
where
A flowchart of the algorithm for the above steps is shown in Fig. 3.
-201811/1001-8042-29-11-018/alternativeImage/1001-8042-29-11-018-F003.jpg)
2.2 Solving the neutron diffusion equation by MLPG
To solve the neutron diffusion equation, the nodal parameters uj and its derivatives should be defined. uj is a function of time, while the shape function
2.2.1 The temporal derivation of
2.2.2 The spatial derivation of
2.2.3 In addition, the second-order spatial domain is:
The derivatives are denoted by subscripts. In addition, a local quadrature domain,
As mentioned earlier,
2.3 Determining the residual of the neutron diffusion equation
In this step, two-dimensional neutron diffusion equations are first introduced. For arbitrary volume, v, the balance neutron equation is considered. This well-known equation is shown below.
[Time rate of change of the number of neutron in a volume]
=[Production rate in the volume]
- [Absorptions in the volume]
- [Net leakage from the surface of volume]
The first term is mathematically expressed in Eq. (21):
The production rate is considered as a neutron source in Eq. (22):
The neutron absorption term is presented in Eq. (23):
The leakage term is shown in Eq. (24):
The mathematical formulation of the neutron diffusion theory is obtained as follows:
Equation (25) can be rewritten as equation (26):
Where u is the solution, D is the diffusion coefficient, s is a neutron source, and
Equation (27) is substituted into equation (20) to give a system of n linear equations for n nodes:
Then, each of the neutron diffusion equation terms should be substituted by Eq.(16) into Eq.(19):
a) The first term is obtained by substituting Eq.(17) into Eq. (28):
b) The second term is obtained by substituting Eq. (19) into Eq. (28):
c) The third term is obtained by substituting Eq. (16) into Eq. (28):
d) The fourth term is obtained by substituting the fission effective factors into Eq. (22):
Therefore, Eq. (28) can be rewritten as:
For clarity, the sums are moved outside the integrals as shown in equation (34):
This equation can be rearranged into its final form as shown in Eq.(35):
It should be noted that the integrals obtained in discretized equations are regular integrals. To calculate the integrals, the numerical integration method based on the 20 point Gauss–Legendre quadrature was applied [33, 45, 58].
2.4 Solution of the time-dependent neutron diffusion equation
The first stage in the time-dependent solution is the discretization of space and time variables. The space and time variables are essentially discretized to gain a series of (nonlinear) algebraic systems. Therefore, it is necessary to decouple the space and time variables to solve this problem.
In this study, three different approaches to time discretization were applied, the forward Euler, Crank-Nicolson, and backward Euler scheme. It is notable that the results of the second scheme were more acceptable than other methods because the method is unconditionally stable. In addition, the time required for this process was shorter than other methods. In the method, time is replaced by the time derivative at the half-step
In this discretization method,
By collecting all the nodal parameters into a vector
where M and N are coefficients of
All approaches were calculated by a code written entirely in FORTRAN 95. To enforce the boundary conditions in the problem, the following provisions should be considered. Each element of matrix L, Lij, is defined based on the three conditions. If xj is in the support domain of boundary node xi, Lij is
In addition, each element of matrix K, Kij, is defined based on the three conditions. If xj is in the support domain of boundary node xi, Lij is zero.
As mentioned previously, the shape function is an important requirement for applying meshless methods. A qualified approach should satisfy several conditions to provide the shape functions of meshless methods [48]. These conditions consist of having an arbitrary nodal distribution, being numerically stable, satisfying a certain order of consistency, and being compactly supported [48]. The accuracy of interpolation for the arbitrary points depends on the number of nodes in the support domain [53]. The number of nodes in each interpolation domain should be used to determine the appropriate shape function and local interpolation schemes to match the sensitivity of any meshless interpolation methods to a variable. Therefore, an appropriate support domain should be selected to apply an accurate approximation [4]. It should be noted that, in order to evaluate the performance of the proposed method, several weight functions such as the cubic spline, the quadratic spline, the compact support radial basis, and Gaussian weight functions were studied. These weight functions are presented below:
1. Cubic spline weight function
2. Quadratic spline weight function
3. Compact support radial basis function (CSRBF)
4. Gaussian weight function
Where
The main steps of the proposed algorithm include:
1. Defining the initial condition of the problem;
2. Defining the boundary condition of the problem;
3. Setting up the cloud points;
4. Finding the affect domain of each node;
5. Looping over the time steps;
6. Appling the MLSA algorithm:
6.1 Selecting the neighboring nodes;
6.2 Deriving the shape functions and its derivative for the nodes;
6.3 Evaluating the nodal matrices and vectors;
6.4 Setting the nodal portion to the global matrices and vectors;
6.5 End loop (see Fig. 3);
7. Solving the government equations at each node;
8. Appling the moving least squares shape function for interior nodes of a local domain for recalculation;
9. Determining the coefficients at each node to obtain higher accuracy;
10. Determining and evaluating the obtained results;
11. Evaluating the effects of increasing the amount of nodes in the problem domain;
12. Recording the history of the state variables and their derivatives;
13. Back to 5;
14. End.
The flowchart of the MLPG algorithm is demonstrated in Fig. 4.
-201811/1001-8042-29-11-018/alternativeImage/1001-8042-29-11-018-F004.jpg)
3 Results and Discussion
Several example problems were studied to validate the proposed method. The results of the MLPG method were then compared to the results of the analytical solution and the Galerkin finite element method.
3.1 Convergence degree test
In order to calculate the convergence degree of the MLPG method, a slab is considered. The boundary conditions are a perfect reflector on the left and a bare surface on the right of the slab. The model geometry is presented in Fig. 5. The logarithm of the error versus the number of nodes is presented in Fig. 6. It is clear that increasing the number of nodes decreases the error. Fig. 7 illustrates the logarithm of the error versus the logarithm of the inverse of the number of nodes. It is well known that the slope of this diagram depicts the convergence rate. For further clarification, the convergence rate versus the number of nodes is presented in Fig. 7b, which demonstrates that the behavior of the proposed method is stable for less than 500 nodes.
-201811/1001-8042-29-11-018/alternativeImage/1001-8042-29-11-018-F005.jpg)
-201811/1001-8042-29-11-018/alternativeImage/1001-8042-29-11-018-F006.jpg)
-201811/1001-8042-29-11-018/alternativeImage/1001-8042-29-11-018-F007.jpg)
3.2 Constant source problem
In this test problem, a distributed source problem was considered. This problem consists of two models, a single region (Fig. 8) and two regions (Fig. 11) including a source in region 1 and a reflector in region 2 with dimensions of 150 cm × 150 cm. The boundary conditions include a perfect reflector on the left side and a bare surface on the right side. In this section, three types of problem were considered:
-201811/1001-8042-29-11-018/alternativeImage/1001-8042-29-11-018-F008.jpg)
-201811/1001-8042-29-11-018/alternativeImage/1001-8042-29-11-018-F011.jpg)
CASE-I: The geometry of the single region is presented in Fig. 8. The parameters required for this calculation are shown in Table 1. The fast and thermal flux of the MLPG method compared to that of GFEM and the analytical solution is presented in Figs. 9 and 10, respectively.
Applied models | Energygroup | Absorptioncross section | Nu-Fissioncross section | Diffusion coefficient |
---|---|---|---|---|
Slab I | 1 | 0.03051 | 0.000256 | 1.440130 |
2 | 0.10582 | 0.154726 | 0.015723 |
-201811/1001-8042-29-11-018/alternativeImage/1001-8042-29-11-018-F009.jpg)
-201811/1001-8042-29-11-018/alternativeImage/1001-8042-29-11-018-F010.jpg)
CASE-II: The geometry of the two-region problem is presented in Fig. 11 and the required parameters are shown in Table 2. A comparison of the fast and the thermal flux of the MLPG method with that of the GFEM and the analytical solution is presented in Figs. 12 and 13, respectively. As shown, the results of the proposed method exhibit good agreement with those of the analytical solution.
Applied models | Energy group | Absorption cross section | Nu-Fission cross section | Diffusion coefficient | Region |
---|---|---|---|---|---|
Slab-II | 1 | 0.03051 | 0.000256 | 1.440130 | Fuel |
2 | 0.10582 | 0.154726 | 0.015723 | ||
1 | 0.035128 | 0.000000 | 1.852524 | Reflector | |
2 | 0.032540 | 0.000000 | 0.294517 |
-201811/1001-8042-29-11-018/alternativeImage/1001-8042-29-11-018-F012.jpg)
-201811/1001-8042-29-11-018/alternativeImage/1001-8042-29-11-018-F013.jpg)
The average error of these results obtained by applying various weight functions is presented in Table 3, which indicates that the accuracy of the solution is improved by using Gaussian weight functions. In Table 4, the average error and CPU time of this problem are demonstrated for the Gaussian and CSRBF weight functions. Fig. 14 presents a comparison of the fast neutron flux in CASE-I and CASE-II obtained by inserting 50×50 nodes.
Error | MLPG Applied Weight Functions | |||
---|---|---|---|---|
Gaussian | CSRBF | Cubic Spline | Quadratic Spline | |
CASE-I | 0.0035 | 0.073 | 0.094 | 0.085 |
CASE-II | 0.0038 | 0.070 | 0.091 | 0.083 |
Problem No. | MLPG Applied Weight Functions | |||
---|---|---|---|---|
Gaussian | CSRBF | |||
Error | CPU Time (s) | Error | CPU Time(s) | |
CASE-I | 0.0035 | 2 | 0.073 | 2.5 |
CASE-II | 0.0038 | 3.5 | 0.070 | 6 |
-201811/1001-8042-29-11-018/alternativeImage/1001-8042-29-11-018-F014.jpg)
In addition, to evaluate the performance of the number of inserted nodes, MLPG methods with 20×20 nodes and 40×40 nodes were compared. As shown, an accurate value of the thermal flux was achieved by the MLPG with 40×40 nodes for CASE-I (Fig. 15). Fig. 15 indicates the performance of the MLPG method was improved by increasing the number of nodes. Fig. 16 demonstrates the contour of the fast and the thermal neutron flux. The best results for CASE-II were obtained with 40×40 nodes.
-201811/1001-8042-29-11-018/alternativeImage/1001-8042-29-11-018-F015.jpg)
-201811/1001-8042-29-11-018/alternativeImage/1001-8042-29-11-018-F016.jpg)
CASE-III: In addition, in order to evaluate the performance of the proposed method for dealing with curve boundaries, a problem with a curved domain was considered. The geometry of this problem is presented in Fig. 17. A comparison of the contours of the thermal neutron flux of this problem is shown in Fig. 18. Two MLPG methods with various numbers of nodes were employed in this problem. The results show good agreement between the proposed method and analytical solutions. In addition, the obtained results illustrate that, to achieve the required accuracy, a large number of nodes can be inserted in the problem domain. Furthermore, increasing the number of nodes in a curved domain by the MLPG method is much simpler than mesh generation for mesh-based methods. However, the integration process for a curved domain is challenging. As mentioned, a parallel algorithm was used in this study to improve the calculation.
-201811/1001-8042-29-11-018/alternativeImage/1001-8042-29-11-018-F017.jpg)
-201811/1001-8042-29-11-018/alternativeImage/1001-8042-29-11-018-F018.jpg)
CASE-IV: Moreover, in order to evaluate the performance of the proposed method, two types of node distribution, normal and scatter models, were applied to CASE-I. The average error of these results obtained by applying various weight functions is presented in Table 4. Table 5 indicates that the accuracy of the solution was improved by using Gaussian weight functions and the scatter distribution. The patterns of the node distributions are presented in Fig. 19.
Distribution Pattern | MLPG Applied Weight Functions | |||
---|---|---|---|---|
Gaussian | CSRBF | Cubic Spline | Quadratic Spline | |
Normal Distribution | 0.0035 | 0.073 | 0.094 | 0.085 |
Scatter Distribution | 0.0030 | 0.059 | 0.082 | 0.065 |
-201811/1001-8042-29-11-018/alternativeImage/1001-8042-29-11-018-F019.jpg)
3.3 Void test problem
The void test problem was first studied by Abuzaid and Gashut [61]. The problem involves a slab including multiple regions with a perfect reflector on the left side and a bare surface on the right side. Fig. 20 presents the geometry and properties of the problem. A comparison between the results of the proposed method and the analytical solution is provided in Fig. 21, which indicates a good agreement.
-201811/1001-8042-29-11-018/alternativeImage/1001-8042-29-11-018-F020.jpg)
-201811/1001-8042-29-11-018/alternativeImage/1001-8042-29-11-018-F021.jpg)
3.4 Neutronic calculations for a SMR reactor
The proposed method was applied to calculate the neutronic parameters in a Small Modular Reactor (SMR). The core used for the analysis was Westinghouse’s SMR, with an electrical output of 200 MW and an active core height of 2.4 m. The core consists of 89 assemblies contained within a core barrel and reactor vessel. These assemblies include 52 fuel assemblies and 37 control rod assemblies. The SMR core map is presented in Fig. 22. The reactor vessel components were based on an AP1000 design but modified to a reduced diameter and height of 3.5 m and 24.7 m, respectively [62].
-201811/1001-8042-29-11-018/alternativeImage/1001-8042-29-11-018-F022.jpg)
The fuel assembly was a square lattice in a standard 17×17 layout with 264 fuel rod locations, 24 guide tube locations, and 1 central location for instrumentation, incorporating the standard Westinghouse design specifications. The fuel of the core was uranium-oxide (UOX) with less than 5% U-235 enrichment (2.5w/o, 3.5w/o, 4.2%w/o U-235 fuel) and a density of 10.36 g/cm3. The fuel rod consisted of a uniform cylindrical pellet stacked together within a Zircaloy clad tube. Between the fuel stack and the clad, clearance was provided to accommodate fuel swelling due to accumulation of fission products, thereby preventing clad rupture. The gap was filled with helium gas to improve heat conduction from fuel to cladding. The guide tubes in the fuel assembly served as a location for the insertion of the rod cluster control assembly (RCCA). The RCCA is a spider assembly consisting of evenly spaced control rods of either silver indium cadmium or boron carbide rods based on the type of fuel used. Detailed specifications for the fuel rod, cladding, structure, control rod, and burnable poison (i.e. discrete and integral) were taken from CASL (Consortium for Advanced Simulation of LWRs) VERA core physics benchmark specifications [62].
This study reveals the ability of the MLPG method to solve the neutronic diffusion equation and evaluates the performance of this method when different selectable parameters were incorporated into the solution. In order to automate the procedure of the proposed method for this problem, a computer program was developed. The first stage of computational procedures consists of creating the cross section database by WIMSD4 code in three-group energy. Then, some of the main neutronic parameters of the core were calculated by the proposed method. The results of the MLPG method were compared with the results of the citation code and the reference values. In addition, as the selection of the weight function and the number of nodes in the problem domain plays an important role in solution accuracy, the performance of the MLPG method was examined with respect to these parameters. Table 6 demonstrates the effective multiplication factor as a main parameter of the core calculation. Tables 7 and 8 present the maximum to average flux ratio and the control rod worth, respectively. The choice of weight function is clearly a significant factor for obtaining an accurate solution. The most important frequently used weight functions are the cubic spline, quadratic spline, CSRBF, and Gaussian. Table 9 compares the average errors of results obtained from the various weight functions and reveals that the Gaussian function produces better accuracy than other functions.
Parameter | Citation | Referenced value | MLPG Applied weight functions | Nodes Insertion | |||
---|---|---|---|---|---|---|---|
Gaussian | CSRBF | Cubic Spline | Quadratic Spline | ||||
Effective Multiplication Factor | 1.12414 | 1.12511 | 1.12332 | 1.11032 | 1.11054 | 1.11902 | 37×37 |
1.12415 | 1.11940 | 1.11891 | 1.11974 | 91×91 | |||
1.12496 | 1.12439 | 1.12475 | 1.11996 | 289×289 |
Parameter | Citation | Referenced value | MLPG Applied weight functions | Nodes Insertion | |||
---|---|---|---|---|---|---|---|
Gaussian | CSRBF | Cubic Spline | Quadratic Spline | ||||
Maximum to Average Flux Ratio | 2.48520 | 2.4980 | 2.48213 | 2.43191 | 2.44521 | 2.39120 | 37×37 |
2.48500 | 2.45225 | 2.46378 | 2.44032 | 91×91 | |||
2.49739 | 2.48821 | 2.49707 | 2.45318 | 289×289 |
Parameter | Citation | Referenced value | MLPG Applied weight functions | Nodes Insertion | |||
---|---|---|---|---|---|---|---|
Gaussian | CSRBF | Cubic spline | Quartic spline | ||||
Control Rod Worth | 0.15613 | 0.1571 | 0.15243 | 0.14254 | 0.1550 | 0.14382 | 37×37 |
0.15615 | 0.14367 | 0.1561 | 0.14630 | 91×91 | |||
0.15694 | 0.14975 | 0.15679 | 0.14587 | 289×289 |
Error | MLPG Applied weight functions | |||
---|---|---|---|---|
Gaussian | CSRBF | Cubic spline | Quartic spline | |
Effective Multiplication Factor | 0.001 | 0.0064 | 0.032 | 0.045 |
Maximum to Average Flux Ratio | 0.002 | 0.041 | 0.037 | 0.179 |
Control Rod Worth | 0.001 | 0.046 | 0.020 | 0.071 |
In addition, the performance of the MLPG method with respect to the number of nodes was evaluated for the above parameters. The values obtained for these parameters are promisingly close to the citation code results and reference values for 289×289 nodes. Moreover, a comparison between the average errors and the CPU time by applying Gaussian weight functions and various node distributions to derive an effective multiplication factor are listed in Table 10. The errors decreases with increasing node number. In addition, the scatter distribution results are more accurate than those of the normal distribution.
Parameter | MLPG Gaussian weight function | |||||
---|---|---|---|---|---|---|
Nodes Normal Distribution | Nodes Scatter Distribution | |||||
Error | CPU Time(s) | Nodes | Error | CPU Time(s) | Nodes | |
Effective Multiplication Factor | 0.012 | 4.0 | 37×37 | 0.010 | 3.5 | 37×37 |
0.087 | 5.5 | 91×91 | 0.063 | 4.0 | 91×91 | |
0.001 | 10 | 289×289 | 0.001 | 7.5 | 289×289 |
The fast flux and thermal flux of the reactor core is demonstrated in Fig. 23 for the MLPG method with various node numbers. In addition, to evaluate the performance of the MLPG method with respect to the number of nodes, 37×37, 91× 91, and 289×289 nodes in the radial direction and 21 layers in the axial direction were employed in the problem domain for any of the fuel assembly types.
-201811/1001-8042-29-11-018/alternativeImage/1001-8042-29-11-018-F023.jpg)
Fig. 24 shows the power factor distributions for a 1/8 symmetry of this core. In this figure, the reference values are compared with the values of the proposed method (MLPG, 289×289 nodes). The reference value exhibits good agreement with the calculated value. The differences between these parameters are also demonstrated in Fig. 24. The comparison of these methods with the citation code and reference values indicates that increasing the number of nodes leads to a more accurate performance of the MLPG method.
-201811/1001-8042-29-11-018/alternativeImage/1001-8042-29-11-018-F024.jpg)
4 Conclusion
In this work, the meshless local Petrov-Galerkin method and moving least squares approximation were coupled to solve the neutron diffusion equation. Several problems were examined to demonstrate the degree of accuracy that can be obtained by employing this method as well as the applicability of the proposed method with respect to effective factors. The results of the proposed method were compared with the analytical solution and the GFEM method, indicating good agreement. Additionally, the proposed method was used to solve neutronic calculations in the SMR reactor. The results were compared to those obtained from the citation code and reference values.
The accuracy and precision of the proposed method were both acceptable and further improved by increasing the number of nodes and through selection of an appropriate weight function. It should be noted that the run-time varied from 0.5 to 10 s, depending mainly on the weight function type, the number of nodes per cell domain, the number of nodes in the boundary, the pattern of node distribution, the type of boundary conditions applied, and the particular computer system on which the code is run. As mentioned above, the numerical proposed method removes the mesh generation process by inserting nodes into the problem domain. Therefore, this study eliminated the mesh generation time plus the interactive and CPU time and decreased the time costs.
In summary, using the MLPG method to solve the neutron diffusion equations is promising due to improved computational accuracy and execution time. Accordingly, the outstanding performance of the MLPG method has been confirmed for nuclear engineering calculations, indicating that it can be successfully applied to the desired nuclear computational codes.
Mesh free methods
. in Encyclopedia of Computational Mechanics, edited byMesh free and particle methods and their applications
. Appl. Mech. Rev. 55, 1-34 (2002). doi: 10.1115/1.1431547Mesh free methods and their comparisons
. Int. J. Comput. Methods 2, 477-515 (2005). doi: 10.1142/S0219876205000673A meshless method for solving two-dimensional variable-order time fractional advection-diffusion equation
. J. Comput. Phys. 340, 655-669 (2017). doi: 10.1016/j.jcp.2017.03.061Transient heat conduction analysis using the MLPG method and modified precise time step integration method
. J. Comput. Phys. 230, 2736-2750 (2011). doi: 10.1016/j.jcp.2011.01.019A stabilized MLPG method for steady state incompressible fluid flow simulation
. J. Comput. Phys. 229, 8564-8577 (2010). doi: 10.1016/j.jcp.2010.08.001A moving least square reproducing polynomial meshless method
. Appl. Numer. Math. 69, 34-58 (2013). doi: 10.1016/j.apnum.2013.03.001Smoothed particle hydrodynamics: theory and application to non-spherical stars
. Mon. Not. R. Astron. Soc. 181, 375-389 (1977). doi: 10.1093/mnras/181.3.375Modeling confined multi-material heat and mass flows using SPH
. Appl. Math. Model. 22, 981-993 (1998). doi: 10.1016/S0307-904X(98)10031-8A point interpolation mesh free method for static and frequency analysis of two-dimensional piezoelectric structures
. Comput. Mech. 29, 510-519 (2002). doi: 10.1007/s00466-002-0360-9Generalizing the finite element method: Diffuse approximation and diffuse elements
. Comput. Mech. 10, 307-318 (1992). doi: 10.1007/BF00364252Moving least-square reproducing kernel methods: part I: methodology and convergence
. Comput. Method, Appl. Mech. 143, 113-154 (1977). doi: 10.1016/S0045-7825(96)01132-2Reproducing kernel particle methods for elastic and plastic problems
. in Advanced Computational Methods for Material Modeling, edited by D. A. Siginer, W. E. VanArsdale, C. M. Altan, and A. N. Alexandrou (ASME, New Orleans, 1993), pp. 175-189Reproducing kernel particle methods
. Numer. Method Fluids 20, 1081-1106 (1995). doi: 10.1002/fld.1650200824A finite point method in computational mechanics - Applications to convective transport and fluid flow
. Int. J. Numer. Methods Eng. 39, 3839-3866 (1996). doi: 10.1002/(SICI)1097-0207(19961130)39:22<3839::AID-NME27>3.0.CO;2-RUsing radial function basis in a boundrary-type meshless method, boundrary point method
. inRadial basis functions: theory and implementions
. (On estimating the curvature attributes and strain invariants of deformed surface through radial basis functions
. Comput. Appl. Math. 1-18 (2016). doi: 10.1007/s40314-016-0380-2A locking-free meshless local Petrov-Galerkin formulation for thick and thin plates
. J. Comput. Phys. 208, 116-133 (2005). doi: 10.1016/j.jcp.2005.02.008Formulation and analysis of a diffusion-velocity particle model for transport-dispersion equations
. 35, 447-473 (2014). doi: 10.1007/s40314-014-0200-5MLPG method for two-dimensional diffusion equation with Neumann’s and non-classical boundary conditions
. Appl. Numer. Math. 61, 170-180 (2011). doi: 10.1016/j.apnum.2010.09.002Application of radial point interpolation method to neutron diffusion field
. Trend Appl. Sci. Res. 7, 18-31 (2012). doi: 10.3923/tasr.2012.18.31A comparison of the meshless RBF collocation method with finite element and boundary element methods in neutron diffusion calculations
. Eng. Anal. Bound. Elem. 46, 30-40 (2014). doi: 10.1016/j.enganabound.2014.05.005Analysis of thick plates by using a higher-order shear and normal deformable plate theory and MLPG method with radial basis functions
. Comput. Methods Appl. Mech. Eng. 196, 979-987 (2007). doi: 10.1016/j.cma.2006.08.002The method of finite spheres
. Comput. Mech. 25, 329-345 (2000). doi: 10.1007/s004660050481Using a meshless kernel-based method to solve the Black - Scholes variational inequality of American options
. Comput. Appl. Math. 1-13 (2016). doi: 10.1007/s40314-016-0351-7Three types of meshless finite volume method for the analysis of two-dimensional elasticity problems
. Comput. Appl. Math. 36, 971 (2015). doi: 10.1007/s40314-015-0273-9New concepts in meshless methods
. Computational Mechanics. Int. J. Numer. Methods Eng. 47, 117-127 (2000). doi: 10.1002/(SICI)1097-0207(20000110/30)47:1/3<537::AID-NME783>3.0.CO;2-ENumerical investigation based on direct meshless local Petrov Galerkin (direct MLPG) method for solving generalized Zakharov system in one and two dimensions and generalized Gross-Pitaevskii equation
. Eng. Comput. 33, 983-996 (2017). doi: 10.1007/s00366-017-0510-5Direct meshless local Petrov-Galerkin method for elliptic interface problems with applications in electrostatic and elastostatic
. Comput. Methods Appl. Mech. Eng. 278, 479-498 (2014). doi: 10.1016/j.cma.2014.05.016A meshless local Petrov-Galerkin method for the time-dependent Maxwell equations
. J. Comput. Appl. Math. 268, 93-110 (2014). doi: 10.1016/j.cam.2014.02.013An adaptive meshless local Petrov-Galerkin method based on a posteriori error estimation for the boundary layer problems
. Appl. Numer. Math. 111, 181-196 (2017). doi: 10.1016/j.apnum.2016.09.007Meshless local Petrov-Galerkin and RBFs collocation methods for solving 2D fractional Klein-Kramers dynamics equation on irregular domains
. C. - Comput. Model. Eng. Sci. 107, 481-516 (2015).Element Free Galerkin Method for threedimensional structural analysis
. Comput. Model. Eng. Sci. 4, 497-508 (2001). doi: 10.3970/cmes.2001.002.497Meshless local Petrov-Galerkin method for two-dimensional nonlinear water wave problems
. J. Comput. Phys. 205, 611-625 (2005). doi: 10.1016/j.jcp.2004.11.010Numerical pricing of American options under two stochastic factor models with jumps using a meshless local Petrov - Galerkin method
. Appl. Numer. Math. 115, 252-274 (2017). doi: 10.1016/j.apnum.2017.01.015The Meshless Local Petrov-Galerkin (MLPG) method for solving incompressible Navier-Stokes equations
. Comput. Model. Eng. Sci. 2, 117-142 (2001). doi: 10.3970/cmes.2001.002.117A meshless local Petrov-Galerkin (MLPG) method for solving the bending problems of a thin plate
. CMES 3, 53-64 (2002). doi: 10.3970/cmes.2002.003.053Element-free galerkin method
. Int. J. Numer. Methods Eng. 37, 229-256 (1994). doi: 10.1002/nme.1620370205Meshless local Petrov-Galerkin (MLPG) method in combination with finite element and boundary element approaches
. Comput. Mech. 26, 536-546 (2000). doi: 10.1007/s004660000203Boundary node Petrov-Galerkin method in solid structures
. Comput. Appl. Math. 1-25 (2016). doi: 10.1007/s40314-016-0335-7Analysis of thermo-elastic problems using the improved element-free Galerkin method
. Comput. Appl. Math. 1-16 (2016). doi: 10.1007/s40314-016-0401-1Direct Meshless Local Petrov - Galerkin (DMLPG) method : A generalized MLS approximation
. Appl. Numer. Math. 68, 73-82 (2013). doi: 10.1016/j.apnum.2013.01.002Meshless Local Petrov - Galerkin (MLPG) method for the unsteady magnetohydrodynamic (MHD) flow through pipe with arbitrary wall conductivity
. Appl. Numer. Math. 59, 1043-1058 (2009). doi: 10.1016/j.apnum.2008.05.001The basis of meshless domain discretization: the meshless local Petrov-Galerkin (MLPG) method
. Adv. Comput. Math. 23, 73-93 (2005). doi: 10.1007/s10444-004-1813-9A meshless local Petrov-Galerkin (MLPG) method for free and forced vibration analyses for solids
. Comput. Mech. 27, 188-198 (2001). doi: 10.1007/s004660100237Surfaces generation by moving least squares methods
. Math. Comput. 37, 141-158 (1981). doi: 10.1090/S0025-5718-1981-0616367-1Local radial point interpolation (MLRPI) method for solving time fractional diffusion-wave equation with damping
. J. Comput. Phys. 312, 307-332 (2016). doi: 10.1016/j.jcp.2016.02.030A local point interpolation method for static and dynamic analysis of thin beams
. Comput. Methods Appl. Mech. Eng. 190, 5515-5528 (2001). doi: 10.1016/S0045-7825(01)00180-3Numerical solution to the unsteady two-dimensional Schrödinger equation using meshless local boundary integral equation method
. Int. J. Numer. Methods Eng. 76, 501-520 (2008). doi: 10.1002/nme.2338Topology optimization of free vibrating continuum structures based on the element free Galerkin method
. Struct. Multidiscip. Optim. 45, 119-127 (2012). doi: 10.1007/s00158-011-0667-2Mesh free method applied to the diffusion equation
. Parallel Numer. 5, 57-66 (2005). doi: 10.1007/s00158-011-0667-2Numerical solution of the multigroup neutron diffusion equation by the meshless RBF collocation method
. Math. Comput. Appl. 18, 399-407 (2013). doi: 10.3390/mca18030399Meshless local Petrov-Galerkin (MLPG) approximation to the two dimensional sine-Gordon equation
. J. Comput. Appl. Math. 233, 2737-2754 (2010). doi: 10.1016/j.cam.2009.11.022A Crank-Nicolson scheme for the Landau-Lifshitz equation without damping
. J. Comput. Appl. Math. 234, 613-623 (2010). doi: 10.1016/j.cam.2010.01.002Identification of exonic regions in DNA sequences using cross-correlation and noise suppression by discrete wavelet transform
. BMC Bioinformatics 12, 430 (2011). doi: 10.1186/1471-2105-12-430Discontinuous Finite Element Methods for Reactor Calculations
. in