Search results for "Galerkin method"
showing 10 items of 71 documents
Numerical analysis of the Oseen-type Peterlin viscoelastic model by the stabilized Lagrange-Galerkin method, Part II: A linear scheme
2017
This is the second part of our error analysis of the stabilized Lagrange-Galerkin scheme applied to the Oseen-type Peterlin viscoelastic model. Our scheme is a combination of the method of characteristics and Brezzi-Pitk\"aranta's stabilization method for the conforming linear elements, which leads to an efficient computation with a small number of degrees of freedom especially in three space dimensions. In this paper, Part II, we apply a semi-implicit time discretization which yields the linear scheme. We concentrate on the diffusive viscoelastic model, i.e. in the constitutive equation for time evolution of the conformation tensor a diffusive effect is included. Under mild stability condi…
Computational aspects in 2D SBEM analysis with domain inelastic actions
2009
The Symmetric Boundary Element Method, applied to structures subjected to temperature and inelastic actions, shows singular domain integrals. In the present paper the strong singularity involved in the domain integrals of the stresses and tractions is removed, and by means of a limiting operation, this traction is evaluated on the boundary. First the weakly singular domain integral in the Somigliana Identity (S.I.) of the displacements is regularized and the singular integral is transformed into a boundary one using the Radial Integration Method; subsequently, using the differential operator applied to the displacement field, the S.I. of the tractions inside the body is obtained and through…
Incremental elastoplastic analysis for active macro-zones
2012
SUMMARY In this paper a strategy to perform incremental elastoplastic analysis using the symmetric Galerkin boundary element method for multidomain type problems is shown. The discretization of the body is performed through substructures, distinguishing the bem-elements characterizing the so-called active macro-zones, where the plastic consistency condition may be violated, and the macro-elements having elastic behaviour only. Incremental analysis uses the well-known concept of self-equilibrium stress field here shown in a discrete form through the introduction of the influence matrix (self-stress matrix). The nonlinear analysis does not use updating of the elastic response inside each plas…
Well-balanced bicharacteristic-based scheme for multilayer shallow water flows including wet/dry fronts
2013
The aim of this paper is to present a new well-balanced finite volume scheme for two-dimensional multilayer shallow water flows including wet/dry fronts. The ideas, presented here for the two-layer model, can be generalized to a multilayer case in a straightforward way. The method developed here is constructed in the framework of the Finite Volume Evolution Galerkin (FVEG) schemes. The FVEG methods couple a finite volume formulation with evolution operators. The latter are constructed using the bicharacteristics of multidimensional hyperbolic systems. However, in the case of multilayer shallow water flows the required eigenstructure of the underlying equations is not readily available. Thus…
Comparison between adaptive and uniform discontinuous Galerkin simulations in dry 2D bubble experiments
2013
Accepted by the Journal of Computational Physics Adaptive mesh refinement generally aims to increase computational efficiency without compromising the accuracy of the numerical solution. However it is an open question in which regions the spatial resolution can actually be coarsened without affecting the accuracy of the result. This question is investigated for a specific example of dry atmospheric convection, namely the simulation of warm air bubbles. For this purpose a novel numerical model is developed that is tailored towards this specific application. The compressible Euler equations are solved with a Discontinuous Galerkin method. Time integration is done with an IMEXmethod and the dy…
Active macro-zone approach for incremental elastoplastic-contact analysis
2013
SUMMARY The symmetric boundary element method, based on the Galerkin hypotheses, has found an application in the nonlinear analysis of plasticity and in contact-detachment problems, but both dealt with separately. In this paper, we want to treat these complex phenomena together as a linear complementarity problem. A mixed variable multidomain approach is utilized in which the substructures are distinguished into macroelements, where elastic behavior is assumed, and bem-elements, where it is possible that plastic strains may occur. Elasticity equations are written for all the substructures, and regularity conditions in weighted (weak) form on the boundary sides and in the nodes (strong) betw…
Monotonic solution of heterogeneous anisotropic diffusion problems
2013
Anisotropic problems arise in various areas of science and engineering, for example groundwater transport and petroleum reservoir simulations. The pure diffusive anisotropic time-dependent transport problem is solved on a finite number of nodes, that are selected inside and on the boundary of the given domain, along with possible internal boundaries connecting some of the nodes. An unstructured triangular mesh, that attains the Generalized Anisotropic Delaunay condition for all the triangle sides, is automatically generated by properly connecting all the nodes, starting from an arbitrary initial one. The control volume of each node is the closed polygon given by the union of the midpoint of…
Galerkin Scheme-Based Determination of Survival Probability of Oscillators With Fractional Derivative Elements
2016
In this paper, an approximate semi-analytical approach is developed for determining the first-passage probability of randomly excited linear and lightly nonlinear oscillators endowed with fractional derivative elements. The amplitude of the system response is modeled as one-dimensional Markovian process by employing a combination of the stochastic averaging and the statistical linearization techniques. This leads to a backward Kolmogorov equation which governs the evolution of the survival probability of the oscillator. Next, an approximate solution of this equation is sought by resorting to a Galerkin scheme. Specifically, a convenient set of confluent hypergeometric functions, related to …
Analysis of a Parabolic Cross-Diffusion Semiconductor Model with Electron-Hole Scattering
2007
The global-in-time existence of non-negative solutions to a parabolic strongly coupled system with mixed Dirichlet–Neumann boundary conditions is shown. The system describes the time evolution of the electron and hole densities in a semiconductor when electron-hole scattering is taken into account. The parabolic equations are coupled to the Poisson equation for the electrostatic potential. Written in the quasi-Fermi potential variables, the diffusion matrix of the parabolic system contains strong cross-diffusion terms and is only positive semi-definite such that the problem is formally of degenerate type. The existence proof is based on the study of a fully discretized version of the system…
Mauro Picone, Sandro Faedo, and the numerical solution of partial differential equations in Italy (1928-1953)
2013
In this paper we revisit the pioneering work on the numerical analysis of partial differential equations (PDEs) by two Italian mathematicians, Mauro Picone (1885-1977) and Sandro Faedo (1913-2001). We argue that while the development of constructive methods for the solution of PDEs was central to Picone's vision of applied mathematics, his own work in this area had relatively little direct influence on the emerging field of modern numerical analysis. We contrast this with Picone's influence through his students and collaborators, in particular on the work of Faedo which, while not the result of immediate applied concerns, turned out to be of lasting importance for the numerical analysis of …