Miscellaneous Weighted Residual Methods

Fragile points method for Euler–Bernoulli beams

This paper briefly describes the development of generalized collocation methods for solving coupled systems of nonlinear, one-dimensional, parabolic, partial differential equations. parabolic, partial differential equations. The methods are applied to several problems that are representative of fluid flow in a porous medium. The numerical results for the problems solved indicate that collocation methods using C piecewise polynomials with Gauss-Legendre collocation points polynomials with Gauss-Legendre collocation points are more efficient than conventional, second-order, finite-difference methods. In particular, for problems whose solutions are smooth, one obtains more accuracy per unit of computer time with increasingly higher-order collocation methods. For problems whose solutions are characterized by a steep, shock-like, transient wave front, one can conclude that increasing the order of the collocation method does not necessarily result in a more efficient method, but high-accuracy collocation methods still appear to be more efficient than conventional, second-order, finite-difference methods. Also, for a given accuracy, collocation methods generally require less computer storage than conventional, second-order, finite-difference methods. Finally, the collocation technique presented here is shown to be not directly applicable to the saturation equation for immiscible flow. Introduction The purpose of this paper is to investigate collocation methods for the numerical solution of nonlinear equations similar to those that describe the flow of a fully compressible fluid in a porous medium and those that describe the process by which one miscible liquid displaces another liquid in a porous medium. The objective is to determine if collocation methods are applicable and economically justified for solving these types of petroleum engineering problems. The conclusions petroleum engineering problems. The conclusions reached are based on the amount of computer time and computer storage that is required to obtain a specified accuracy. Since only two problems are solved, the results may not be extendable to more general problems; however, experience in solving a wide variety of problems indicates that this is the general behavior that might be expected when using collocation methods. Previously published results indicate how Galerkin, collocation, and finite-element methods can be used to solve various problems that arise in petroleum engineering. Most of these papers present a variety of arguments in favor of these present a variety of arguments in favor of these methods. The principal argument is that better answers can be obtained for the same computational effort than by finite-difference methods. Culham and Varga have presented the most convincing evidence in favor of Galerkin methods. However, they only considered piecewise linear and piecewise cubic-basis functions and they only considered one problem. Because of these limitations, it is not problem. Because of these limitations, it is not clear whether their conclusions can be extended to higher-order Galerkin methods or to other types of problems. Also, they did not consider collocation problems. Also, they did not consider collocation methods in any detail. This paper considers the use of increasingly higher-order collocation methods and compares the results with those obtained with the usual second-order, finite-difference method. The Galerkin method is not considered, even though that is the method with which collocation must ultimately compete. This was done for several reasons. A preliminary operation count indicates that Galerkin methods will not be competitive with collocation methods for most problems. Also, it is not clear how to implement the higher-order Galerkin methods to achieve maximum efficiency. These remarks indicate that the relative merits of collocation and Galerkin methods are still uncertain. PROBLEMS SOLVED PROBLEMS SOLVED Two problems are solved to serve as a basis for the comparisons. The first problem, called the "gas flow problem," is the same problem considered by Culham and Varga. This problem is nonlinear and includes a volumetric source term. SPEJ P. 345

The use of mesh-based methods for modeling elongated composite structures with inhomogeneities leads to an increase in computational costs when discretizing the waveguide part of greater linear dimensions, while the semi-analytical numerical methods do not allow one to directly describe the structures with local inhomogeneities of an arbitrary shape. To compensate for the shortcomings of these two classes of numerical methods, we propose a hybrid scheme based on the spectral element method (SEM) and the semi-analytical finite element method (SAFEM) for studying the anti-plane vibrations of a composite structure in the frequency domain. Thus, for a waveguide, this scheme makes it possible to represent the solution via the SAFEM as a sum of modes or guided waves, and the adjacent regions are discretized using the SEM. The displacement and stress continuity conditions are imposed on the common boundary of two domains. To couple the solutions, we introduce an auxiliary function of displacement, which is approximated by applying the same basis functions as those used in SEM and SAFEM (Lagrange interpolation polynomials on the Gauss–Legendre–Lobatto nodal points). The unknown expansion coefficients of this function are determined by the Galerkin and collocation methods. It has been established that both methods provide the same accuracy. The results obtained by the hybrid approach employing the Galerkin and collocation methods are compared with the results calculated in the standard finite element software. It is shown that they are in good agreement. The presented hybrid approach can be straightforwardly extended to the case of in-plane motion, but it requires significant refinement for a three-dimensional case.

Multiscalets in the multiwavelet family are used as the basis and testing functions in Galerkin's method. Since the multiscalets are orthogonal to their translations under the Sobolev inner product, the resulting Galerkin's method behaves like a collocation method but possesses the ability of derivative tracking for unknown functions in solving integral equations. The former makes the method simple in implementation and the latter allows to use coarse meshes in discretization. These robust features have been demonstrated in solving two‐dimensional (2D) electromagnetic (EM) problems, but have not been exploited in three‐dimensional (3D) scenarios. For 3D problems, the unknown functions in the integral equations are dependent on two coordinate variables. In order to preserve the use of coarse meshes for 3D cases, we realize the omnidirectional derivative tracking by tracking the directional derivatives along two orthogonal directions, or equivalently tracking the gradient. This process yields a nonsquare matrix equation and we use the least‐squares method (LSM) to solve it. Numerical examples show that the multiscalet‐based Galerkin's method is also robust in solving for 3D EM integral equations with a minor cost increase from LSM. Copyright © 2007 John Wiley & Sons, Ltd.

This paper discusses spectral and spectral element methods with Legendre-Gauss-Lobatto nodal basis for general 2nd-order elliptic eigenvalue problems. The special work of this paper is as follows. (1) We prove a priori and a posteriori error estimates for spectral and spectral element methods. (2) We compare between spectral methods, spectral element methods, finite element methods and their derivedp-version,h-version, andhp-version methods from accuracy, degree of freedom, and stability and verify that spectral methods and spectral element methods are highly efficient computational methods.

The finite volume method in the context of the finite element method

Exponentially accurate spectral and spectral element methods for fractional ODEs

This study presents a rigorous and in-depth comparative analysis of the Chebyshev Collocation and Spectral Galerkin methods for solving nonlinear singular initial value problems, which are fundamental in modeling astrophysical phenomena, with a focus on the Lane-Emden equation as a model test case, by applying both methodologies, performance was evaluated based on precise metrics including accuracy, rate of convergence, computational efficiency, the results unequivocally demonstrate the superiority of the Spectral Galerkin method across all critical performance aspects, the method achieved true exponential convergence for smooth solutions, reaching a precision approaching machine accuracy using significantly lower approximation degrees than the collocation method. Furthermore, the Galerkin method exhibited superior robustness and a higher algebraic rate of convergence in cases where the solution is less regular, while maintaining excellent numerical stability and better-conditioned algebraic systems, we conclude that the integral nature of the orthogonal projection in the Galerkin method endows it with a structural advantage

On variationally consistent versus heuristic mass formulations in cut and extended finite element methods

Time/wave domain analysis for axially moving pre-stressed nanobeam by wavelet-based spectral element method

Spectral element method in time for rapidly actuated systems

Spectral methods are evaluated for the solution of population balance problems. Both a simplified population balance equation (PBE) with an analytical solution available and a rigorous PBE are considered in this numerical analysis. In comparison with the orthogonal collocation, tau, and Galerkin methods, the least-squares method does not produce the same favorable results. On the other hand, the least-squares method with a direct minimization solver is capable of producing equally accurate results of the model equations as the orthogonal collocation, tau, and Galerkin methods. Compared to the conventional least-squares method, the direct minimization formulation is better conditioned but does not produce a symmetric positive-definite system matrix.

One of the well-known methods for the determination of wind energy potential is the two-parameter Weibull distribution. It is clear that the success of the Weibull distribution for wind energy applications depends on the estimation of the parameters which can be determined by using various numerical methods. In the present study, Monte Carlo simulation method is performed by using six parameters estimation method that is used in the estimation of Weibull distribution parameters such as Maximum Likelihood Estimation (MLE), Least Squares Method (LSM), Method of Moments (MOM), Method of Logarithmic Moments (MLM), Percentile Method (PM), and L-Moment Method (LM), and is compared to Mean Squared Error (MSE) and Mean Absolute Percentage Error (MAPE). In this study, the wind energy potential of the Meşelik region in Eskişehir was modeled with two-parameter Weibull distribution. The average wind speed (m/s) data, which are gathered in 10-minute intervals from the measuring device installed 10 meters about the ground in Meşelik Campus of Eskişehir Osmangazi University, is used. As a result of the simulation study, it has been determined that MLE is the best parameter estimation method for two-parameter Weibull distribution in large sample sizes, and LM has the closest performance to MLE. The wind speed (m/h) data of the region has been successfully modeled with two-parameter Weibull distribution and the highest average wind power density has been obtained in July as 49.38295 (W/m2) while the lowest average wind power density has been obtained in October as 19.30044 (W/m2).

On the comparison of semi-analytical methods for the stability analysis of delay differential equations

The Lagrange–Galerkin Spectral Element Method on Unstructured Quadrilateral Grids