Calcolo最新文献

This work considers the numerical analysis of nonlinear Volterra integrodifferential equations that arise from, e.g., the theory of isotropic viscoelastic rods and plates. Novel properties of the kernel and its derivatives are derived via the Laplace transform, and the regularity of the solutions is proved. Then we apply the Crank–Nicolson method and the second-order convolution quadrature rule to develop the discrete-in-time scheme, and the energy argument is employed to analyze the stability and convergence of the proposed scheme. Numerical experiments are performed to substantiate the theoretical findings.

The (C^1)-(P_{k+1}) ((kge 4)) Argyris finite elements combined with the (C^0)-(P_k) Lagrange finite elements are locking-free with respect to the plate thickness, and quasi-optimal when solving the Reissner–Mindlin plate equation on triangular meshes. The method is truly conforming or consistent in the sense that no reduction operator is introduced to the formulation. Theoretical proof and numerical verification are presented.

In this manuscript, we aim to study the semi-analytical and the numerical solution of a semilinear time-fractional diffusion equation where the time-fractional term includes the combination of tempered fractional derivative and k-Caputo fractional derivative with a parameter (k ge 1). The application of the new integral transform, namely Elzaki transform of the tempered k-Caputo fractional derivative is shown here and thereafter the semi-analytical solution is obtained by using the Elzaki decomposition method. The model problem is linearized using Newton’s quasilinearization method, and then the quasilinearized problem is discretized by the difference scheme namely tempered (_kL2)-(1_sigma ) method. Stability and convergence analysis of the proposed scheme have been discussed in the (L_2)-norm using the energy method. In support of the theoretical results, numerical example has been incorporated.

In this paper, we study the lower and upper bounds for Stokes eigenvalues by finite element schemes. For the schemes studied here, roughly speaking, the loss of the local approximation property of the discrete velocity and pressure spaces may lead to different computed bounds of the eigenvalues. Formally theoretical analysis is constructed based on certain mathematical hypotheses, and numerical experiments are given to illustrate the validity of the theory.

We present an algorithm for Tchakaloff-like compression of quasi-Monte Carlo volume and surface integration on an arbitrary union of balls, via non-negative least squares. We also provide the corresponding Matlab codes together with several numerical tests.

In this study, we investigate a hybrid-type anisotropic weakly over-penalised symmetric interior penalty method for the Poisson equation on convex domains. Compared with the well-known hybrid discontinuous Galerkin methods, our approach is simple and easy to implement. Our primary contributions are the proposal of a new scheme and the demonstration of a proof for the consistency term, which allows us to estimate the anisotropic consistency error. The key idea of the proof is to apply the relation between the Raviart–Thomas finite element space and a discontinuous space. In numerical experiments, we compare the calculation results for standard and anisotropic mesh partitions.

The Douglas-Rachford method (DR) is one of the most computationally efficient iterative methods for the large scale linear systems of equations. Based on the randomized alternating reflection and relaxation strategy, we propose a randomized block Douglas–Rachford method for solving the matrix equation (AXB=C). The Polyak’s and Nesterov-type momentums are integrated into the randomized block Douglas–Rachford method to improve the convergence behaviour. The linear convergence of the resulting algorithms are proven. Numerical simulations and experiments of randomly generated data, real-world sparse data, image restoration problem and tensor product surface fitting in computer-aided geometry design are performed to illustrate the feasibility and efficiency of the proposed methods.

In this paper we study the convergence of a finite volume approximation of a convective diffusive elliptic problem with Neumann boundary conditions and (L^1) data. To deal with the non-coercive character of the equation and the low regularity of the right hand-side we mix the finite volume tools and the renormalized techniques. To handle the Neumann boundary conditions we choose solutions having a null median and we prove a convergence result. We present also some numerical experiments in dimension 2 to illustrate the rate of convergence.

This paper presents an extension of the transpose-free quasi-minimal residual (TFQMR) method for solving the generalized coupled Sylvester tensor equations. The new algorithm is based on the tensor format of the TFQMR process. We analyze the convergence behavior of this method and present a bound for the residual norm of the method depending on the specific parameter computed by the algorithm. The numerical experiments demonstrate the efficiency of the new method and confirm the theoretical results.