discontinuous galerkin finite element methods for linear elasticity and quasistatic linear viscoelasticity we consider a finite element in space and quadrature in time discretization of a compressible linear quasistatic viscoelasticity problem the spatial discretization uses a discontinous galerkin finite element method based on polynomials of degree r termed dg r and the time discretization uses a trapezoidal rectangle rule approximation to the volterra history integral both semi and fully discrete a priori error estimates are derived without recourse to gronwall s inequality and therefore the error bounds do not show exponential growth in time moreover the convergence rates are optimal in both h and r providing that the finite element space contains a globally continuous interpolant to the exact solution e g when using the standard k polynomial basis on simplicies or tensor product polynomials k on quadrilaterals when this is not the case e g using k on quadri laterals the convergence rate is suboptimal in r but remains optimal in h we also consider a reduction of the problem to standard linear elasticity where similarly optimal a priori error estimates are derived for the dg r approximation

isogeometric analysis using t splines we explore t splines a generalization of nurbs enabling local refinement as a basis for isogeometric analysis we review t splines as a surface design methodology and then develop it for engineering analysis applications we test t splines on some elementary two dimensional and three dimensional fluid and structural analysis problems and attain good results in all cases we summarize the current status of t splines their limitations and future possibilities

isogeometric numerical dispersion analysis for two dimensional elastic wave propagation abstract in this paper we carry out a numerical dispersion analysis for the linear two dimensional elastodynamics equations approximated by means of nurbs based isogeometric analysis in the framework of the galerkin method specifically we consider the analysis of harmonic plane waves in an isotropic and homogeneous elastic medium we compare and discuss the errors associated with the compressional and shear wave velocities and we provide the anisotropic curves for numerical approximations obtained by considering b spline and nurbs basis functions of different regularity namely globally c 0 and c p 1 continuous p being the polynomial degree we conclude our analysis by numerically simulating the seismic wave propagation in a sinusoidal shaped valley with discontinuous elastic parameters across an internal interface

a mortar spectral finite element method for complex 2d and 3d elastodynamic problems abstract we present a hybrid spectral element finite element domain decomposition method for solving elastic wave propagation problems aim of the method is to exploit both the enhanced accuracy of spectral elements allowing significant reductions of the computational load and the flexibility of finite elements in treating irregular geometries and non linear media interfaces between finite and spectral elements are managed by the mortar method which enjoys optimal accuracy and allows to discretize the different regions independently several test cases illustrate the robustness of this method and its tolerance to a significant range of geometrical models characterized by different grid refinements

new fractional step runge kutta nystrom methods up to order three abstract fractional step runge kutta nystro fsrkn methods have been revealed to be an excellent option to integrate numerically many multidimensional evolution models governed by second order in time partial differential equations these methods combined with suitable spatial discretizations lead to strong computational cost reductions respect to many classical implicit time integrators in this paper we present the construction process of several implicit fsrkn methods of two and three levels which attain orders up to three and satisfy adequate stability properties we have also performed some numerical experiments in order to show the unconditionally convergent behavior of these schemes as well as their computational advantages

isogeometric analysis of the advective cahn hilliard equation spinodal decomposition under shear flow we present a numerical study of the spinodal decomposition of a binary fluid undergoing shear flow using the advective cahn hilliard equation a stiff nonlinear parabolic equation characterized by the presence of fourth order spatial derivatives our numerical solution procedure is based on isogeometric analysis an approximation technique for which basis functions of high order continuity are employed these basis functions allow us to directly discretize the advective cahn hilliard equation without resorting to a mixed formulation we present steady state solutions for rectangular domains in two dimensions and for the first time in three dimensions we also present steady state solutions for the two dimensional taylor couette cell to enforce periodic boundary conditions in this curved domain we derive and utilize a new periodic bezier extraction operator we present an extensive numerical study showing the effects of shear rate surface tension and the geometry of the domain on the phase evolution of the binary fluid theoretical and experimental results are compared with our simulations

discontinuous galerkin finite element method for the wave equation the symmetric interior penalty discontinuous galerkin finite element method is presented for the numerical discretization of the second order wave equation the resulting stiffness matrix is symmetric positive definite and the mass matrix is essentially diagonal hence the method is inherently parallel and leads to fully explicit time integration when coupled with an explicit time stepping scheme optimal a priori error bounds are derived in the energy norm and the l 2 norm for the semidiscrete formulation in particular the error in the energy norm is shown to converge with the optimal order cal o h min s ell with respect to the mesh size h the polynomial degree ell and the regularity exponent s of the continuous solution under additional regularity assumptions the l 2 error is shown to converge with the optimal order cal o h min s ell numerical results confirm the expected convergence rates and illustrate the versatility of the method

duality and unified analysis of discrete approximations in structural dynamics and wave propagation comparison of p method finite elements with k method nurbs we study the discretization behavior of classical finite element and nurbs approximations on problems of structural vibrations and wave propagation we find that on the basis of equal numbers of degrees of freedom and bandwidth nurbs have superior approximation properties in fact we observe that the high mode behavior of classical finite elements is divergent with the order of approximation a surprisingly negative result on the other hand nurbs offer almost spectral approximation properties and all modes converge with increasing order of approximation

on the stability of the rayleigh ritz method for eigenvalues this paper studies global stability properties of the rayleigh ritz approximation of eigenvalues of the laplace operator the focus lies on the ratios hat lambda k lambda k of the kth numerical eigenvalue hat lambda k and the kth exact eigenvalue lambda k in the context of classical finite elements the maximal ratio blows up with the polynomial degree for b splines of maximum smoothness the ratios are uniformly bounded with respect to the degree except for a few instable numerical eigenvalues which are related to the presence of essential boundary conditions these phenomena are linked to the inverse inequalities in the respective approximation spaces

optimal discontinuous galerkin methods for wave propagation we have developed and analyzed a new class of discontinuous galerkin methods dg which can be seen as a compromise between standard dg and the finite element fe method in the way that it is explicit like standard dg and energy conserving like fe in the literature there are many methods that achieve some of the goals of explicit time marching unstructured grid energy conservation and optimal higher order accuracy but as far as we know only our new algorithms satisfy all the conditions we propose a new stability requirement for our dg the stability analysis is based on the careful selection of the two fe spaces which verify the new stability condition the convergence rate is optimal with respect to the order of the polynomials in the fe spaces moreover the convergence is described by a series of numerical experiments

robustness of isogeometric structural discretizations under severe mesh distortion abstract this paper investigates higher order and higher continuity functions in isogeometric structural analysis under distortion of the control and physical meshes first the concepts behind isogeometric analysis are briefly reviewed the effect of perturbing control points upon the geometrical mapping is observed in the one dimensional setting the concept is extended to multiple dimensions with a series of cases where the interior of a cube is distorted in different ways and patch tests are performed to assess accuracy a circular fillet patch is introduced through the appropriate degeneration of a control mesh in order to address the problem of modeling small features in finite element models the fillet patch is tested with two problems of linear elasticity and is shown to be accurate the investigation is further extended to the setting of large deformation nearly incompressible nonlinear elasticity with two examples which make use of the recently proposed f projection method

a direct method for the computation of hopf bifurcation points we present a method for the numerical computation of a hopf bifurcation point of an evolution equation an augmented time independent system of equations has to be solved for which the hopf bifurcation point is an isolated solution in the case of a system of ordinary differential equations we propose an efficient implementation of newton s method for the solution of the resulting system of algebraic equations if the method is applied directly to a partial differential system with one space variable a two point boundary value problem is obtained we also compare this technique with a related method and results are given for two systems of partial differential equations describing chemical reacting systems

isogeometric analysis of geometric partial differential equations abstract we consider the numerical approximation of geometric partial differential equations pdes defined on surfaces in the 3d space in particular we focus on the geometric pdes deriving from the minimization of an energy functional by l 2 gradient flow we analyze two energy functionals the area which leads to the mean curvature flow a nonlinear second order pde and the willmore energy leading to the willmore flow a nonlinear fourth order pde we consider surfaces represented by single patch tensor product nurbs and discretize the pdes by means of nurbs based isogeometric analysis in the framework of the galerkin method to approximate the high order geometric pdes we use high order continuous nurbs basis functions for the time discretization of the nonlinear geometric pdes we use backward differentiation formulas bdf with extrapolation of the geometric quantities involved in the weak formulation of the problem in this manner we solve a linear problem at each time step we report numerical results concerning the mean curvature and willmore flows on different geometries of interest and we show the accuracy and efficiency of the proposed approximation scheme

isogeometric analysis of the electrophysiology in the human heart numerical simulation of the bidomain equations on the atria abstract we consider isogeometric analysis iga for the numerical solution of the electrophysiology of the atria which in this work is modeled by means of the bidomain equations on thin surfaces first we consider the bidomain equations coupled with the roger mcculloch ionic model on simple slabs here our goal is to evaluate the effects of the spatial discretization by iga and the use of different b spline basis functions on the accuracy of the approximation in particular regarding the accuracy of the front velocity and the dispersion error specifically we consider basis functions with high polynomial degree p and global high order continuity c p 1 in the computational domain our results show that the use of such basis functions is beneficial to the accurate approximation of the solution then we consider a realistic application of the bidomain equations coupled with the courtemanche ramirez nattel ionic model on the two human atria which are represented by means of two nurbs surfaces

non conforming high order approximations of the elastodynamics equation in this paper we formulate and analyze two non conforming high order strategies for the approximation of elastic wave problems in heterogeneous media namely the mortar spectral element method and the discontinuous galerkin spectral element method starting from a common variational formulation we make a full comparison of the two techniques from the points of view of accuracy convergence grid dispersion and stability

isogeometric analysis of the cahn hilliard phase field model abstract the cahn hilliard equation involves fourth order spatial derivatives finite element solutions are not common because primal variational formulations of fourth order operators are only well defined and integrable if the finite element basis functions are piecewise smooth and globally c 1 continuous there are a very limited number of two dimensional finite elements possessing c 1 continuity applicable to complex geometries but none in three dimensions we propose isogeometric analysis as a technology that possesses a unique combination of attributes for complex problems involving higher order differential operators namely higher order accuracy robustness two and three dimensional geometric flexibility compact support and most importantly the possibility of c 1 and higher order continuity a nurbs based variational formulation for the cahn hilliard equation was tested on two and three dimensional problems we present steady state solutions in two dimensions and for the first time in three dimensions to achieve these results an adaptive time stepping method is introduced we also present a technique for desensitizing calculations to dependence on mesh refinement this enables the calculation of topologically correct solutions on coarse meshes opening the way to practical engineering applications of phase field methodology

discrete dispersion relation for hp version finite element approximation at high wave number the dispersive properties of high order finite element schemes are analyzed in the setting of the helmholtz equation and an explicit form of the discrete dispersion relation is obtained for elements of arbitrary order it is shown that the numerical dispersion displays three different types of behavior depending on the size of the order of the method relative to the mesh size and the wave number quantitative estimates are obtained for the behavior and rates of decay of the dispersion error in the differing regimes all estimates are fully explicit and are shown to be sharp limits are obtained where transitions between the different regimes occur and used to provide guidelines for the selection of the mesh size and the polynomial order in terms of the wave number so that the dispersion error is controlled

studies of refinement and continuity in isogeometric structural analysis we investigate the effects of smoothness of basis functions on solution accuracy within the isogeometric analysis framework we consider two simple one dimensional structural eigenvalue problems and two static shell boundary value problems modeled with trivariate nurbs solids we also develop a local refinement strategy that we utilize in one of the shell analyses we find that increased smoothness that is the k method leads to a significant increase in accuracy for the problems of structural vibrations over the classical c0c0 continuous p method whereas a judicious insertion of c0c0 continuous surfaces about singularities in a mesh otherwise generated by the k method usually outperforms a mesh in which all basis functions attain their maximum level of smoothness we conclude that the potential for the k method is high but smoothness is an issue that is not well understood due to the historical dominance of c0c0 continuous finite elements and therefore further studies are warranted

finite dimensional approximation of nonlinear problems part ii limit points note univ paris 6 f 75230 paris 05 france ecole polytech ctr math appl cnrs equipe rech 747 f 91128 palaiseau france univ pavia cnr ist anal numer i 27100 pavia italy brezzi f univ pavia cnr ist math appl i 27100 pavia italy isi document delivery no lr761times cited 54cited reference count 13 reference asn article 1981 002doi 10 1007 bf01396184 record created on 2006 08 24 modified on 2017 05 12

3d stochastic gravity inversion using nonstationary covariances abstractthe flexibility of geostatistical inversions in geophysics is limited by the use of stationary covariances which implicitly and mostly for mathematical convenience assumes statistical homogeneity of the studied field for fields showing sharp contrasts due for example to faults or folds an approach based on the use of nonstationary covariances for cokriging inversion was developed the approach was tested on two synthetic cases and one real data set inversion results based on the nonstationary covariance were compared to the results from the stationary covariance for two synthetic models the nonstationary covariance better recovered the known synthetic models with the real data set the nonstationary assumption resulted in a better match with the known surface geology

an algorithm for the computation of hopf bifurcation points in comparison with other methods we present an algorithm for the computation of a hopf bifurcation point based on a direct method i e an augmented time independent system is solved the algorithm is designed to deal efficiently with problems representing a finite difference discretization of a partial differential system in this case the bandstructure of the jacobian matrix is exploited we also compare this technique with some other methods in regard to its efficiency and reliability

finite dimensional approximation of nonlinear problems in the first two papers of this series 4 5 we have studied a general method of approximation of nonsingular solutions and simple limit points of nonlinear equations in a banach space we derive here general approximation results of the branches of solutions in the neighborhood of a simple bifurcation point the abstract theory is applied to the galerkin approximation of nonlinear variational problems and to a mixed finite element approximation of the von karman equations

numerical approximation of a periodic linear parabolic problem we approximate a periodic linear parabolic problem using a galerkin method for the space variable and a spectral method for the time variable some optimal error estimates are derived

finite element and nurbs approximations of eigenvalue boundary value and initial value problems abstract we study the spectral approximation properties of finite element and nurbs spaces from a global perspective we focus on eigenfunction approximations and discover that the l 2 norm errors for finite element eigenfunctions exhibit pronounced spikes about the transition points between branches of the eigenvalue spectrum this pathology is absent in nurbs approximations by way of the pythagorean eigenvalue error theorem we determine that the squares of the energy norm errors of the eigenfunctions are the sums of the eigenvalue errors and the squares of the l 2 norm eigenfunction errors the spurious behavior of the higher eigenvalues for standard finite elements is well known and therefore inherited by the energy norm errors along with the spikes in the l 2 norm of the eigenfunction errors the eigenvalue pathology is absent for nurbs the implications of these results to the corresponding elliptic boundary value problem and parabolic and hyperbolic initial value problems are discussed

approximation of hopf bifurcation we make several assumptions on a nonlinear evolution problem ensuring the existence of a hopf bifurcation under a fairly general approximation condition we define a discrete problem which retains the bifurcation property and we prove an error estimate between the branches of exact and approximate periodic solutions

efficient quadrature for nurbs based isogeometric analysis abstract we initiate the study of efficient quadrature rules for nurbs based isogeometric analysis a rule of thumb emerges the half point rule indicating that optimal rules involve a number of points roughly equal to half the number of degrees of freedom or equivalently half the number of basis functions of the space under consideration the half point rule is independent of the polynomial order of the basis efficient rules require taking into account the precise smoothness of basis functions across element boundaries several rules of practical interest are obtained and a numerical procedure for determining efficient rules is presented we compare the cost of quadrature for typical situations arising in structural mechanics and fluid dynamics the new rules represent improvements over those used previously in isogeometric analysis

isogeometric collocation cost comparison with galerkin methods and extension to adaptive hierarchical nurbs discretizations abstract we compare isogeometric collocation with isogeometric galerkin and standard c 0 finite element methods with respect to the cost of forming the matrix and residual vector the cost of direct and iterative solvers the accuracy versus degrees of freedom and the accuracy versus computing time on this basis we show that isogeometric collocation has the potential to increase the computational efficiency of isogeometric analysis and to outperform both isogeometric galerkin and standard c 0 finite element methods when a specified level of accuracy is to be achieved with minimum computational cost we then explore an adaptive isogeometric collocation method that is based on local hierarchical refinement of nurbs basis functions and collocation points derived from the corresponding multi level greville abscissae we introduce the concept of weighted collocation that can be consistently developed from the weighted residual form and the two scale relation of b splines using weighted collocation in the transition regions between hierarchical levels we are able to reliably handle coincident collocation points that naturally occur for multi level greville abscissae the resulting method combines the favorable properties of isogeometric collocation and hierarchical refinement in terms of computational efficiency local adaptivity robustness and straightforward implementation which we illustrate by numerical examples in one two and three dimensions

a simple algorithm for obtaining nearly optimal quadrature rules for nurbs based isogeometric analysis we develop new quadrature rules for isogeometric analysis based on the solution of a local nonlinear problem a simple and robust algorithm is developed to determine the rules which are exact for important b spline spaces of uniform and geometrically stretched knot spacings we consider both periodic and open knot vector configurations and illustrate the efficiency of the rules on selected boundary value problems we find that the rules are almost optimally efficient but much easier to obtain than optimal rules which require the solution of global nonlinear problems that are often ill posed