An efficient parallel algorithm of variational nodal method for heterogeneous neutron transport problems

The heterogeneous variational nodal method(HVNM)has emerged as a potential approach for solving high-fidelity neutron transport problems.However,achieving accurate results with HVNM in large-scale problems using high-fidelity models has been challenging due to the prohibitive computational costs.This paper presents an efficient parallel algorithm tailored for HVNM based on the Message Passing Interface standard.The algorithm evenly distributes the response matrix sets among processors during the matrix formation process,thus enabling independent construction without communication.Once the formation tasks are completed,a collective operation merges and shares the matrix sets among the processors.For the solution process,the problem domain is decomposed into subdomains assigned to specific processors,and the red-black Gauss-Seidel iteration is employed within each subdomain to solve the response matrix equation.Point-to-point communication is conducted between adjacent subdomains to exchange data along the boundaries.The accuracy and efficiency of the parallel algorithm are verified using the KAIST and JRR-3 test cases.Numerical results obtained with multiple processors agree well with those obtained from Monte Carlo calculations.The parallelization of HVNM results in eigenvalue errors of 31 pcm/-90 pcm and fission rate RMS errors of 1.22%/0.66%,respectively,for the 3D KAIST problem and the 3D JRR-3 problem.In addition,the parallel algorithm significantly reduces computation time,with an efficiency of 68.51% using 36 processors in the KAIST problem and 77.14% using 144 processors in the JRR-3 problem.

Weak and Strong Convergence of Self Adaptive Inertial Subgradient Extragradient Algorithms for Solving Variational Inequality Problems

Many solutions of variational inequalities have been proposed,among which the subgradient extragradient method has obvious advantages.Two different algorithms are given for solving variational inequality problem in this paper.The problem we study is defined in a real Hilbert space and has L-Lipschitz and pseudomonotone condition.Two new algorithms adopt inertial technology and non-monotonic step size rule,and their convergence can still be proved when the value of L is not given in advance.Finally,some numerical results are designed to demonstrate the computational efficiency of our two new algorithms.

A RELAXED INERTIAL FACTOR OF THE MODIFIED SUBGRADIENT EXTRAGRADIENT METHOD FOR SOLVING PSEUDO MONOTONE VARIATIONAL INEQUALITIES IN HILBERT SPACES

In this paper,we investigate pseudomonotone and Lipschitz continuous variational inequalities in real Hilbert spaces.For solving this problem,we propose a new method that combines the advantages of the subgradient extragradient method and the projection contraction method.Some very recent papers have considered different inertial algorithms which allowed the inertial factor is chosen in[0;1].The purpose of this work is to continue working in this direction,we propose another inertial subgradient extragradient method that the inertial factor can be chosen in a special case to be 1.Under suitable mild conditions,we establish the weak convergence of the proposed algorithm.Moreover,linear convergence is obtained under strong pseudomonotonicity and Lipschitz continuity assumptions.Finally,some numerical illustrations are given to confirm the theoretical analysis.

Variational Iteration Method for Solving Time Fractional Burgers Equation Using Maple

The Time Fractional Burger equation was solved in this study using the Mabel software and the Variational Iteration approach. where a number of instances of the Time Fractional Burger Equation were handled using this technique. Tables and images were used to present the collected numerical results. The difference between the exact and numerical solutions demonstrates the effectiveness of the Mabel program’s solution, as well as the accuracy and closeness of the results this method produced. It also demonstrates the Mabel program’s ability to quickly and effectively produce the numerical solution.

A Dimension-Splitting Variational Multiscale Element-Free Galerkin Method for Three-Dimensional Singularly Perturbed Convection-Diffusion Problems

By introducing the dimensional splitting(DS)method into the multiscale interpolating element-free Galerkin(VMIEFG)method,a dimension-splitting multiscale interpolating element-free Galerkin(DS-VMIEFG)method is proposed for three-dimensional(3D)singular perturbed convection-diffusion(SPCD)problems.In the DSVMIEFG method,the 3D problem is decomposed into a series of 2D problems by the DS method,and the discrete equations on the 2D splitting surface are obtained by the VMIEFG method.The improved interpolation-type moving least squares(IIMLS)method is used to construct shape functions in the weak form and to combine 2D discrete equations into a global system of discrete equations for the three-dimensional SPCD problems.The solved numerical example verifies the effectiveness of the method in this paper for the 3D SPCD problems.The numerical solution will gradually converge to the analytical solution with the increase in the number of nodes.For extremely small singular diffusion coefficients,the numerical solution will avoid numerical oscillation and has high computational stability.

Regularization Methods to Approximate Solutions of Variational Inequalities

In this paper, we study the regularization methods to approximate the solutions of the variational inequalities with monotone hemi-continuous operator having perturbed operators arbitrary. Detail, we shall study regularization methods to approximate solutions of following variational inequalities: and with operator A being monotone hemi-continuous form real Banach reflexive X into its dual space X*, but instead of knowing the exact data (y0, A), we only know its approximate data satisfying certain specified conditions and D is a nonempty convex closed subset of X;the real function f defined on X is assumed to be lower semi-continuous, convex and is not identical to infinity. At the same time, we will evaluate the convergence rate of the approximate solution. The regularization methods here are different from the previous ones.

A Variational Multiscale Method for Particle Dispersion Modeling in the Atmosphere

A LES model is proposed to predict the dispersion of particles in the atmosphere in the context of Chemical,Biological,Radiological and Nuclear(CBRN)applications.The code relies on the Finite Element Method(FEM)for both the fluid and the dispersed solid phases.Starting from the Navier-Stokes equations and a general description of the FEM strategy,the Streamline Upwind Petrov-Galerkin(SUPG)method is formulated putting some emphasis on the related assembly matrix and stabilization coefficients.Then,the Variational Multiscale Method(VMS)is presented together with a detailed illustration of its algorithm and hierarchy of computational steps.It is demonstrated that the VMS can be considered as a more general version of the SUPG method.The final part of the work is used to assess the reliability of the implemented predictor/multicorrector solution strategy.

Aerodynamic Design of Airfoils Based on Variable-Domain Variational Finite Element Method

Designing airfoils according to given pressure (or velocity) distribution is one kind of free boundary problems. Free boundary condition can be coupled with the flow governing equations by variable-domain variational calculus, which makes it possible to calculate simultaneously the flow field and the free boundary. An accurate deduction of the variable-domain variational principles is taken herein to design airfoils in compressible and incompressible flows. Furthermore, two grid types (H and O) are used in the calculation with better results for the O-type grid. It is shown that convergence is accelerated and good results can be obtained even if the initial guessed airfoil shape is a triangle, demonstrating the strong adaptability of this method.

Variational iteration method for solving the mechanism of the Equatorial Eastern Pacific El Nino-Southern Oscillation

A class of coupled system for the E1 Nifio-Southern Oscillation （ENSO） mechanism is studied. Using the method of variational iteration for perturbation theory, the asymptotic expansions of the solution for ENSO model are obtained and the asymptotic behaviour of solution for corresponding problem is considered.

Variational iteration solving method for El Nio phenomenon atmospheric physics of nonlinear model

A class of E1 Niйo atmospheric physics oscillation model is considered. The E1 Niйo atmospheric physics oscillation is an abnormal phenomenon involved in the tropical Pacific ocean-atmosphere interactions. The conceptual oscillator model should consider the variations of both the eastern and western Pacific anomaly patterns. An E1 Niйo atmospheric physics model is proposed using a method for the variational iteration theory. Using the variational iteration method, the approximate expansions of the solution of corresponding problem are constructed. That is, firstly, introducing a set of functional and accounting their variationals, the Lagrange multiplicators are counted, and then the variational iteration is defined, finally, the approximate solution is obtained. From approximate expansions of the solution, the zonal sea surface temperature anomaly in the equatorial eastern Pacific and the thermocline depth anomaly of the sea-air oscillation for E1 Niйo atmospheric physics model can be analyzed. E1 Niйo is a very complicated natural phenomenon. Hence basic models need to be reduced for the sea-air oscillator and are solved. The variational iteration is a simple and valid approximate method.

Variational regularization method of solving the Cauchy problem for Laplace's equation: Innovation of the Grad–Shafranov(GS) reconstruction

The simplified linear model of Grad-Shafranov （GS） reconstruction can be reformulated into an inverse boundary value problem of Laplace＇s equation. Therefore, in this paper we focus on the method of solving the inverse boundary value problem of Laplace＇s equation. In the first place, the variational regularization method is used to deal with the ill- posedness of the Cauchy problem for Laplace＇s equation. Then, the ＇L-Curve＇ principle is suggested to be adopted in choosing the optimal regularization parameter. Finally, a numerical experiment is implemented with a section of Neumann and Dirichlet boundary conditions with observation errors. The results well converge to the exact solution of the problem, which proves the efficiency and robustness of the proposed method. When the order of observation error δ is 10-1, the order of the approximate result error can reach 10-3.

A NEW ITERATIVE METHOD FOR FINDING COMMON SOLUTIONS OF GENERALIZED EQUILIBRIUM PROBLEM,FIXED POINT PROBLEM OF INFINITE k-STRICT PSEUDO-CONTRACTIVE MAPPINGS,AND QUASI-VARIATIONAL INCLUSION PROBLEM

In this article, we introduce a hybrid iterative scheme for finding a common element of the set of solutions for a generalized equilibrium problems, the set of common fixed point for a family of infinite k-strict pseudo-contractive mappings, and the set of solutions of the variational inclusion problem with multi-valued maximal monotone mappings and inverse-strongly monotone mappings in Hilbert space. Under suitable conditions, some strong convergence theorems are proved. Our results extends the recent results in G.L.Acedo and H.K.Xu [2], Zhang, Lee and Chan [8], Wakahashi and Toyoda [9], Takahashi and Takahashi [I0] and S. S. Chang, H. W. Joseph Lee and C. K. Chan [II], S.Takahashi and W.Takahashi [12]. Moreover, the method of proof adopted in this article is different from those of [4] and [12].

SEMI-INVERSE METHOD AND GENERALIZED VARIATIONAL PRINCIPLES WITH MULTI-VARIABLES IN ELASTICITY

Semi-inverse method, which is an integration and an extension of Hu's try-and-error method, Chien's veighted residual method and Liu's systematic method, is proposed to establish generalized variational principles with multi-variables without arty variational crisis phenomenon. The method is to construct an energy trial-functional with an unknown function F, which can be readily identified by making the trial-functional stationary and using known constraint equations. As a result generalized variational principles with two kinds of independent variables (such as well-known Hellinger-Reissner variational principle and Hu-Washizu principle) and generalized variational principles with three kinds of independent variables (such as Chien's generalized variational principles) in elasticity have been deduced without using Lagrange multiplier method. By semi-inverse method, the author has also proved that Hu-Washizu principle is actually a variational principle with only two kinds of independent variables, stress-strain relations are still its constraints.

Numerical manifold method for thermo-mechanical coupling simulation of fractured rock mass

As a calculation method based on the Galerkin variation,the numerical manifold method(NMM)adopts a double covering system,which can easily deal with discontinuous deformation problems and has a high calculation accuracy.Aiming at the thermo-mechanical(TM)coupling problem of fractured rock masses,this study uses the NMM to simulate the processes of crack initiation and propagation in a rock mass under the influence of temperature field,deduces related system equations,and proposes a penalty function method to deal with boundary conditions.Numerical examples are employed to confirm the effectiveness and high accuracy of this method.By the thermal stress analysis of a thick-walled cylinder(TWC),the simulation of cracking in the TWC under heating and cooling conditions,and the simulation of thermal cracking of the SwedishÄspöPillar Stability Experiment(APSE)rock column,the thermal stress,and TM coupling are obtained.The numerical simulation results are in good agreement with the test data and other numerical results,thus verifying the effectiveness of the NMM in dealing with thermal stress and crack propagation problems of fractured rock masses.

Modified Lagrange Multiplier Method and Generalized Variational Principle in Fluid Mechanics

The Lagrange multiplier method plays an important role in establishing generalized variational principles notonly in tluid mechallics. but also in elasticity. Sometimes, however, one may come across variational crisis(somemultipliers vanish identically). failing to achieve his aim. The crisis is caused by the fact that the Inultipliers are treatedas independent variables in the process of variatioll. but after identification they become functions of the originalindependent variables. To overcome it, a Inodified Lagrange multiplier method or semi-inverse method has beenproposed to deduce generalized varistional principles. Some e-camples are given to illustrate its convenience andeffectiveness of the novel method.

AN ITERATIVE METHOD FOR THE DISCRETE PROBLEMS OF A CLASS OF ELLIPTICAL VARIATIONAL INEQUALITIES

Based on the nonlinear characiers of the discrete problems of some ellipticalvariational inequalities, this paper presents a numerical iterative method, the schemesof which are pithy and converge rapidly The new method possesses a high efficiency. insolving such applied engineering problems as obstacle problems and .free boundary.problems arising in fluid lubrications.

GENERALIZED VARIATIONAL DATA ASSIMILATION METHOD AND NUMERICAL EXPERIMENT FOR NON-DIFFERENTIAL SYSTEM

The generalized variational data assimilation for non-differential dynamical systems is studied.There is no tangent linear model for non-differential systems and thus the general adjoint model can not be derived in the traditional way.The weak form of the original system was introduced, and then the generalized adjoint model was derived. The generalized variational data assimilation methods were developed for non-differential low dimensional system and non-differential high dimensional system with global and local observations. Furthermore, ideas in inverse problems are introduced to 4DVAR (Four-dimensional variational) of non-differential partial differential system with local observations.

Inertial Subgradient Extragradient Algorithm for Solving Variational Inequality Problems with Pseudomonotonicity

In order to solve variational inequality problems of pseudomonotonicity and Lipschitz continuity in Hilbert spaces, an inertial subgradient extragradient algorithm is proposed by virtue of non-monotone stepsizes. Moreover, weak convergence and R-linear convergence analyses of the algorithm are constructed under appropriate assumptions. Finally, the efficiency of the proposed algorithm is demonstrated through numerical implementations.

Methods of Variations and Their Applications

In this review article, we begin with reviewing Calculus of variations giving few examples on its use to solve a large number of problems in geometry, physics, and other branches of knowledge. Afterwards, we direct our attention to different methods of variations which evolved during the last century and which include their use in eigenvalue problems and in finite difference methods and those adopted in classical and quantum mechanics. The methods used in evaluating products and quotients of functionals are also discussed along with variational iteration methods. Later on, a good number of applications in different areas are presented and discussed;then a concluding discussion is given.

A Unification of the Concepts of the Variational Iteration, Adomian Decomposition and Picard Iteration Methods;and a Local Variational Iteration Method

This paper compares the variational iteration method(VIM),the Adomian decomposition method(ADM)and the Picard iteration method(PIM)for solving a system of first o rder n onlinear o rdinary d ifferential e quations(ODEs).A unification of the concepts underlying these three methods is attempted by considering a very general iterative algorithm for VIM.It is found that all the three methods can be regarded as special cases of using a very general matrix of Lagrange multipliers in the iterative algorithm of VIM.The global variational iteration method is briefly reviewed,and further recast into a Local VIM,which is much more convenient and capable of predicting long term complex dynamic responses of nonlinear systems even if they are chaotic.