Finite-difference approximation of energies in fracture mechanics
The superconsistent collocation method, which is based on a collocation grid different from the one used to represent the solution, has proven to be very accurate in the resolution of various functional equations. Excellent results can be also obtained for what concerns preconditioning. Some analysis and numerous experiments, regarding the use of finite-differences preconditioners, for matrices arising from pseudospectral approximations of advection-diffusion boundary value problems, are presented...
We propose and analyze several finite-element schemes for solving a grade-two fluid model, with a tangential boundary condition, in a two-dimensional polygon. The exact problem is split into a generalized Stokes problem and a transport equation, in such a way that it always has a solution without restriction on the shape of the domain and on the size of the data. The first scheme uses divergence-free discrete velocities and a centered discretization of the transport term, whereas the other schemes...
We propose and analyze several finite-element schemes for solving a grade-two fluid model, with a tangential boundary condition, in a two-dimensional polygon. The exact problem is split into a generalized Stokes problem and a transport equation, in such a way that it always has a solution without restriction on the shape of the domain and on the size of the data. The first scheme uses divergence-free discrete velocities and a centered discretization of the transport term, whereas the other schemes...
We present, in a uniform manner, several integral equations of the first kind for the solution of the two-dimensional interior Dirichlet boundary value problem. We apply a general numerical collocation method to the various equations, and thereby we compare the various integral equations, and recommend two of them. We give a survey of the various numerical methods, and present a simple method for the numerical solution of the recommended integral equations.
In this work we consider the dual-primal Discontinuous Petrov–Galerkin (DPG) method for the advection-diffusion model problem. Since in the DPG method both mixed internal variables are discontinuous, a static condensation procedure can be carried out, leading to a single-field nonconforming discretization scheme. For this latter formulation, we propose a flux-upwind stabilization technique to deal with the advection-dominated case. The resulting scheme is conservative and satisfies a discrete maximum...
In this work we consider the dual-primal Discontinuous Petrov–Galerkin (DPG) method for the advection-diffusion model problem. Since in the DPG method both mixed internal variables are discontinuous, a static condensation procedure can be carried out, leading to a single-field nonconforming discretization scheme. For this latter formulation, we propose a flux-upwind stabilization technique to deal with the advection-dominated case. The resulting scheme is conservative and satisfies a discrete...
The paper is concerned with the graph formulation of forced anisotropic mean curvature flow in the context of the heteroepitaxial growth of quantum dots. The problem is generalized by including anisotropy by means of Finsler metrics. A semi-discrete numerical scheme based on the method of lines is presented. Computational results with various anisotropy settings are shown and discussed.
We propose and analyse a abstract framework for augmented mixed formulations. We give a priori error estimate in the general case: conforming and nonconforming approximations with or without numerical integration. Finally, a posteriori error estimator is given. An example of stabilized formulation for Stokes problem is analysed.
There is a growing interest in high-order finite and spectral/hp element methods using continuous and discontinuous Galerkin formulations. In this paper we investigate the effect of h- and p-type refinement on the relationship between runtime performance and solution accuracy. The broad spectrum of possible domain discretisations makes establishing a performance-optimal selection non-trivial. Through comparing the runtime of different implementations...
The Poisson-Boltzmann (PB) model is an effective approach for the electrostatics analysis of solvated biomolecules. The nonlinearity associated with the PB equation is critical when the underlying electrostatic potential is strong, but is extremely difficult to solve numerically. In this paper, we construct two operator splitting alternating direction implicit (ADI) schemes to efficiently and stably solve the nonlinear PB equation in a pseudo-transient continuation approach. The operator splitting...
We consider incremental problem arising in elasto-plastic models with isotropic hardening. Our goal is to derive computable and guaranteed bounds of the difference between the exact solution and any function in the admissible (energy) class of the problem considered. Such estimates are obtained by an advanced version of the variational approach earlier used for linear boundary-value problems and nonlinear variational problems with convex functionals [24, 30]. They do no contain mesh-dependent constants...