Legendre-Weighted Residual Methods for System of Fractional Order Differential Equations

The numerical approach for finding the solution of fractional order systems of boundary value problems (BPVs) is derived in this paper. The implementation of the weighted residuals such as Galerkin, Least Square, and Collocation methods are included for solving fractional order differential equations, which is broadened to acquire the approximate solutions of fractional order systems with differentiable polynomials, namely Legendre polynomials, as basis functions. The algorithm of the residual formulations of matrix form can be coded efficiently. The interpretation of Caputo fractional derivatives is employed here. We have demonstrated these methods numerically through a few examples of linear and nonlinear BVPs. The results in absolute errors show that the present method efficiently finds the numerical solutions of fractional order systems of differential equations.

A WEIGHTED RESIDUAL METHOD FOR ELASTIC-PLASTIC ANALYSIS NEAR A CRACK TIP AND THE CALCULATION OF THE PLASTIC STRESS INTENSITY FACTORS

In this paper, a weighted residual method for the elastic-plastic analysis near a crack tip is systematically given by taking the model of power-law hardening under plane strain condition as a sample. The elastic-plastic solutions of the crack lip field and an approach based on the superposition of the nonlinear finite element method on the complete solution in the whole crack body field, to calculate the plastic stress intensity factors, are also developed. Therefore, a complete analvsis based on the calculation both for the crack tip field and for the whole crack body field is provided.

Improvement of the minimal residual method for solving nonsymmetric linear systems

In this paper, the minimal residual （MRES） method for solving nonsymmetric equation systems was improved, the recurrence relation was deduced between the approximate solutions of the linear equation system Ax = b, and a more effective method was presented, which can reduce the operational count and the storage.

Method of Semi Analytic Perturbation Weighted Residuals for Nonlinear Bending of Shallow Shells

The semi? analytic perturbation weighted residuals method was used to solve the nonlinear bending problem of shallow shells, and the fifth order B spline was taken as trial function to seek an efficient method for nonlinear bending problem of shallow shells. The results from the present method are in good agreement with those derived from other methods. The present method is of higher accuracy, lower computing time and wider adaptability. In addition, the design of computer program is simple and it is easy to be programmed.

A Combination of Residual Distribution and the Active Flux Formulations or a New Class of Schemes That Can Combine Several Writings of the Same Hyperbolic Problem:Application to the 1D Euler Equations

We show how to combine in a natural way(i.e.,without any test nor switch)the conservative and non-conservative formulations of an hyperbolic system that has a conservative form.This is inspired from two different classes of schemes:the residual distribution one(Abgrall in Commun Appl Math Comput 2(3):341–368,2020),and the active flux formulations(Eyman and Roe in 49th AIAA Aerospace Science Meeting,2011;Eyman in active flux.PhD thesis,University of Michigan,2013;Helzel et al.in J Sci Comput 80(3):35–61,2019;Barsukow in J Sci Comput 86(1):paper No.3,34,2021;Roe in J Sci Comput 73:1094–1114,2017).The solution is globally continuous,and as in the active flux method,described by a combination of point values and average values.Unlike the“classical”active flux methods,the meaning of the point-wise and cell average degrees of freedom is different,and hence follow different forms of PDEs;it is a conservative version of the cell average,and a possibly non-conservative one for the points.This new class of scheme is proved to satisfy a Lax-Wendroff-like theorem.We also develop a method to perform nonlinear stability.We illustrate the behaviour on several benchmarks,some quite challenging.

High-precision downward continuation of potential fi elds algorithm utilizing adaptive damping coeffi cient of generalized minimal residuals

The downward continuation of potential fields is a process of calculating their values in a lower plane based on those of a certain plane.This technology is not only a data processing method for resource exploration but also plays an extremely important role in military applications.However,the downward continuation of potential fields is a typical linear inverse problem that is ill-posed.Generalized minimal residuals(GMRES)is an eff ective solution to ill-posed inverse problems,but it is unstable under the condition wherein the GMRES is directly applied in the calculation process.Moreover,the long-term behavior of its iterative computation is a disordered,divergent result.Therefore,to obtain stable solutions,GMRES is applied to solve the normal equations of the downward continuation of potential fields;it is also used to prequalify for occasional interruptions in the operation process by adding the damping coefficient,thus strengthening the stability conditions of the equations of residual minimization.Finally,the stable downward continuation of the potential fields method is proposed.As indicated by the theoretical data and the measured testing data,the method proposed in this paper has the advantages of high-precision and excellent stability.Furthermore,compared with the Tikhonov iteration method,the proposed method avoids the need to choose regularization parameters.

An analysis of the incompressible viscous flows problem by meshless method

Generally the incompressible viscous flow problem is described by the Navier-Stokes equation. Based on the weighted residual method the discrete formulation of element-free Galerkin is inferred in this paper. By the step-bystep computation in the field of time, and adopting the least-square estimation of the-same-order shift, this paper has calculated both velocity and pressure from the decoupling independent equations. Each time fraction Newton-Raphson iterative method is applied for the velocity and pressure. Finally, this paper puts the method into practice of the shear-drive cavity flow, verifying the validity, high accuracy and stability.

AN ITERATIVE PARALLEL ALGORITHM OF FINITE ELEMENT METHOD

In this paper, a parallel algorithm with iterative form for solving finite element equation is presented. Based on the iterative solution of linear algebra equations, the parallel computational steps are introduced in this method. Also by using the weighted residual method and choosing the appropriate weighting functions, the finite element basic form of parallel algorithm is deduced. The program of this algorithm has been realized on the ELXSI-6400 parallel computer of Xi'an Jiaotong University. The computational results show the operational speed will be raised and the CPU time will be cut down effectively. So this method is one kind of effective parallel algorithm for solving the finite element equations of large-scale structures.

Fully coupled flow-induced vibration of structures under small deformation with GMRES method

Lagrangian-Eulerian formulations based on a generalized variational principle of fluid-solid coupling dynamics are established to describe flow-induced vibration of a structure under small deformation in an incompressible viscous fluid flow. The spatial discretization of the formulations is based on the multi-linear interpolating functions by using the finite element method for both the fluid and solid structures. The generalized trapezoidal rule is used to obtain apparently non-symmetric linear equations in an incremental form for the variables of the flow and vibration. The nonlinear convective term and time factors are contained in the non-symmetric coefficient matrix of the equations. The generalized minimum residual （GMRES） method is used to solve the incremental equations. A new stable algorithm of GMRES-Hughes-Newmark is developed to deal with the flow-induced vibration with dynamical fluid-structure interaction in complex geometries. Good agreement between the simulations and laboratory measurements of the pressure and blade vibration accelerations in a hydro turbine passage was obtained, indicating that the GiViRES-Hughes-Newmark algorithm presented in this paper is suitable for dealing with the flow-induced vibration of structures under small deformation.

TECHNOLOGY OF QUANTITATIVE ANALYSIS ON VIBRATORY STRESS RELIEF BASED ON ULTRASONIC TIME-OF-ARRIVAL METHOD

The effect of vibratory stress relief （VSR） is usually evaluated with the indirect method of observing the change of amplitude frequency response characteristics of structures. A new kind of evaluating method of VSR based on the ultrasonic time-of-arrival method （UTM）, which can obtain the residual stress directly through measuring the propagation time of ultrasonic wave in the material, is presented. At first, the principle of the measuring method of residual stress based on UTM is analyzed. Then the measuring system of the method is described, which is in virtue of ultrasonic flaw detector and high-sampling-rate digital oscillograph. And a set of calibration system that contains a piece of standard specimen is also introduced. Experimental results prove the relation between the residual stress and the propagation time of ultrasonic in workpieces. Finally, the measuring and calibration systems are applied in evaluating the effect of VSR. The final test results show that the method is effective.

APPLICATION OF SUBREGION FUNCTION METHOD FOR SOLVING BEAM-BOARD STRUCTURE

In this paper, the solution to the structure consisting of a bead and a board is given as a result of the application of the subregion function method which was suggested in ref. [1]. The same problem is also computed with finite element method. The comparison between the two results shows that the application of the subregion function in the method of weighted residuals is practical and effective, especially for solving compound structures.

AN APPROXIMATE METHOD FOR DETERMING THE VELOCITY PROFILE IN A LAMINAR BOUNDARY-LAYER ON FLAT PLATE

In this paper, using the integration method, it is sought to solve the problem for the laminar boundary_layer on a flat plate. At first, a trial function of the velocity profile which satisfies the basical boundary conditions is selected. The coefficients in the trial function awaiting decision are decided by using some numerical results of the boundary_layer differential equations. It is similar to the method proposed by Peng Yichuan, but the former is simpler. According to the method proposed by Peng, when the awaiting decision coefficients of the trial function are decided, it is sought to solve a third power algebraic equation. On the other hand, in this paper, there is only need for solving a linear algebraic equation. Moreover, the accuracy of the results of this paper is higher than that of Peng.

Vortex-Source Method for the 3D Incompressible Irrotational Flow

The weighted residuals method was used for obtaining the boundary integral representation of the velocity of the three-dimensional inviscid irrotational flow. It is shown that velocity in an arbitrary point of domain can be expressed through its values on the boundary. Boundary integral equations of the second kind for solving boundary-valued problems of the first and second kinds are developed. The result has been also generalised to the case of solenoidal vector fields with potential vorticity. It is shown that the resulting integral equations are Fredholm integral equations of the second kind and allow effective numerical solving of corresponding boundary-valued problems. Examples of numerical solutions for a sphere and an ellipsoid are given for demonstration of efficiency of the offered method.

Analytical Solution for Evaluation of Voltage Surge Distribution along Transformer Windings through the Method of Residues

Computer programs have definitely become indispensable for designing power transformer. Among several applications, computer programs are mostly used for electric field calculation and thus electrical insulation concerns. In consequence, studies based on analytical approach to basic studies of correlated problems have become even more important because they form the very basis of knowledge that is necessary to every transformer designer in view of taking all the advantages of computational analyses. On the other hand, one of the most important basic studies consists in the evaluation of voltage surge distribution along transformer windings for which the method of separation of variables has been extensively used thanks to some simplifying assumptions. With this aim, authors have developed and previously published works that show the applicability of an alternative and useful analytical method that is the method of the residues, which requires no simplification to be assumed. In this work, another important step is taken towards proofing the total applicability of this promising method that is through a practical problem. A comparison to the numerical method TLM （transmission line method） is also performed and concordance with TLM and experimental data confirms the proposal of the method of residues can be also applicable to several others problems of electromagnetism.

Paleotectonic residual strain energy in Southwest China

Based on the orthotropic elastic theory of rock masses, the X-ray method was used to measure the distribution of macro-residual strain energy density along a depth profile,using core samples taken from 47 large-aperture deep boreholes in four regions of Southwest China： the Longmenshan, Anninghe, Honghe, and Xianshuihe fault zones.Then, the vertical gradients of the macro-residual strain energy density and the macroresidual strain energy contained in high-energy cuboid block segments along each fault zone were determined. The results demonstrate that the macro-residual strain energy stored at shallow levels in the rock mass in these fault zones may be partly responsible for generating many large earthquakes and may explain why the large earthquakes in this region are typically shallow.

The Fractional Investigation of Fornberg-Whitham Equation Using an Efficient Technique

In the last few decades,it has become increasingly clear that fractional calculus always plays a very significant role in various branches of applied sciences.For this reason,fractional partial differential equations(FPDEs)are of more importance to model the different physical processes in nature more accurately.Therefore,the analytical or numerical solutions to these problems are taken into serious consideration and several techniques or algorithms have been developed for their solution.In the current work,the idea of fractional calculus has been used,and fractional FornbergWhithamequation(FFWE)is represented in its fractional view analysis.Awell-knownmethod which is residual power series method(RPSM),is then implemented to solve FFWE.TheRPSMresults are discussed through graphs and tables which conform to the higher accuracy of the proposed technique.The solutions at different fractional orders are obtained and shown to be convergent toward an integer-order solution.Because the RPSM procedure is simple and straightforward,it can be extended to solve other FPDEs and their systems.

The Fractional Investigation of Some Nonlinear Partial Differential Equations by Using an Efficient Procedure

The nonlinearity inmany problems occurs because of the complexity of the given physical phenomena.The present paper investigates the non-linear fractional partial differential equations’solutions using the Caputo operator with Laplace residual power seriesmethod.It is found that the present technique has a direct and simple implementation to solve the targeted problems.The comparison of the obtained solutions has been done with actual solutions to the problems.The fractional-order solutions are presented and considered to be the focal point of this research article.The results of the proposed technique are highly accurate and provide useful information about the actual dynamics of each problem.Because of the simple implementation,the present technique can be extended to solve other important fractional order problems.

Residue Based Open-Loop Modal Analysis Method for Detecting LFMR of PMSG-WFs Penetrated Power Systems

This paper proposes a residue based open-loop modal analysis method to detect low frequency modal resonance(LFMR),including asymmetric low frequency modal attraction(ALFMA)and asymmetric low frequency modal repulsion(ALFMR),of permanent magnetic synchronous generator based wind farms(PMSG-WFs)penetrated power systems.The formation of ALFMA and ALFMR caused by two open-loop low frequency oscillation(LFO)modes moving close and apart is analyzed in detail.Via predicting the trajectories of closed-loop LFO modes based on calculation of residue of open-loop LFO modes,both ALFMA and ALFMR can be detected.The proposed method can select LFO modes which move to the right half complex plane as control parameters vary.Simulation studies are carried out on a three-machine power system and a four-machine 11-bus power system to verify the properties of the proposed method.

Expansion modelling of discrete grey model based on multi-factor information aggregation

This paper aims to study a novel expansion discrete grey forecasting model, which could aggregate input information more effectively. In general, existing multi-factor grey forecasting models, such as one order and h variables grey forecasting model （GM （1, h））, always aggregate the main system variable and independent variables in a linear form rather than a nonlinear form, while a nonlinear form could be used in more cases than the linear form. And the nonlinear form could aggregate collinear independent factors, which widely lie in many multi-factor forecasting problems. To overcome this problem, a new approach, named as the Solow residual method, is proposed to aggregate independent factors. And a new expansion model, feedback multi-factor discrete grey forecasting model based on the Solow residual method （abbreviated as FDGM （1, h））, is proposed accordingly. Then the feedback control equation and the parameters＇ solution of the FDGM （1, h） model are given. Finally, a real application is used to test the modelling accuracy of the FDGM （1, h） model. Results show that the FDGM （1, h） model is much better than the nonhomogeneous discrete grey forecasting model （NDGM） and the GM （1, h） model.

Application of A Fast Multipole BIEM for Flow Diffraction from A 3D Body

A Fast Multipole Method (FMM) is developed as a numerical approach to the reduction of the computational cost and requirement memory capacity for a large in solving large-scale problems. In this paper it is applied to the boundary integral equation method (BIEM) for current diffraction from arbitrary 3D bodies. The boundary integral equation is discretized by higher order elements, the FMM is applied to avoid the matrix/vector product, and the resulting algebraic equation is solved by the Generalized Conjugate Residual method (GCR). Numerical examination shows that the FMM is more efficient than the direct evaluation method in computational cost and storage of computers.