a domain decomposition method based on weighted interior penalties for advection diffusion reaction problems we propose a domain decomposition method for advection diffusion reaction equations based on nitsche s transmission conditions the advection dominated case is stabilized using a continuous interior penalty approach based on the jumps in the gradient over element boundaries we prove the convergence of the finite element solutions of the discrete problem to the exact solution and propose a parallelizable iterative method the convergence of the resulting domain decomposition method is proved and this result holds true uniformly with respect to the diffusion parameter the numerical scheme that we propose here can thus be applied straightforwardly to diffusion dominated advection dominated and hyperbolic problems some numerical examples are presented in different flow regimes showing the influence of the stabilization parameter on the performance of the iterative method and we compare our method with some other domain decomposition techniques for advection diffusion equations

a robust nitsche s formulation for interface problems abstract in this work we propose a novel weighting for the interfacial consistency terms arising in a nitsche variational form we demonstrate through numerical analysis and extensive numerical evidence that the choice of the weighting parameter has a great bearing on the stability of the method consequently we propose a weighting that results in an estimate for the stabilization parameter such that the method remains well behaved in varied settings ranging from the configuration of embedded interfaces resulting in arbitrarily small elements to such cases where a large contrast in material properties exists an important consequence of this weighting is that the bulk as well as the interfacial fields remain well behaved in the presence of a elements with arbitrarily small volume fractions b large material heterogeneities and c both large heterogeneities as well as arbitrarily small elements we then highlight the accuracy and efficiency of the proposed formulation through numerical examples focusing particular attention on interfacial quantities of interest

robust computational algorithms for dynamic interface tracking in three dimensions front tracking provides sharp resolution of wave fronts through the active tracking of interfaces between distinct materials a major challenge to this method is to handle changes in the interface topology we describe two algorithms implemented in the front tracking code frontier to model dynamic changes in three dimensional interfaces the two methods can be combined to give a hybrid method that is superior to each individual method the success of these algorithms is shown by simulations of rayleigh taylor instability which is an interfacial instability driven by an acceleration directed across a material interface our numerical results are validated by comparing the numerical computation of the velocity of a single rising bubble with an analytic model for the bubble velocity

variational and projection methods for the volume constraint in finite deformation elasto plasticity this paper focuses on the treatment of volume constraints which in the context of elasto plasticity typically arise as a result of assuming volume preserving plastic flow projection methods based on the modification of the discrete gradient operator b often proposed on an ad hoc basis are systematically obtained in the variational context furnished by a three field hu washizu principle the fully nonlinear formulation proposed here is based on a local multiplicative split of the deformation gradient into volume preserving and dilatational parts without relying on rate forms of the weak form of momentum balance this approach fits naturally in a formulation of plasticity based on the multiplicative decomposition of the deformation gradient and enables one to exactly enforce the condition of volume preserving plastic flow within the framework proposed in this paper rate forms and incrementally objective algorithms are entirely bypassed

total positivity and its applications preface i spline functions total positivity of the spline kernel and its applications b bojanov total positivity properties of lb splines t lyche l l schumaker total positivity and splines k morken totally positive systems and interpolation by subspaces of generalized splines m sommer h strauss ii matrix theory preservation of total nonnegativity under the hadamard product and related topics j garloff d g wagner vertex implications for totally nonnegative under the hadamard product and related topics j garloff on factorizations of totally positive matrices j garloff on factorizations of totally positive matrices m gasca j m pena iii geometric modeling total positivity and optimal bases j m carnicer j m pena total positivity and the shape of curves t n t goodman tchebycheff curves m l mazure h pottmann iv probability and mathematical biology analytical and probabilistic methods for totally positive matrices g s goodman some inequalities of total positivity in pure and applied mathematics s karlin z ziegler new directions and problems of multivariate total positivity motivated by genetic recombination probabilities s karlin v approximation theory tchebycheff spaces and total positivity j m carnicer j m pena cebysev and weak cebysev systems r a zalik vi complex analysis zero sets of entire fourier transforms of multiply positive density functions o m katkova vii statistics total nonnegativity of moment matrices and its application to e optimal designs in polynomial regression b heiligers analysis of regression dependence in 2 xn bivariate distributions and some applications in contingency tables s kasala m b rao viii real analysis total positivity and nonlinear analysis r q jia solutions of certain functional equations s karlin alternance and snakes in the nonlinear finite moment problem k tatalian ix combinatorics the applications of total positivity to combinatorics and conversely f brenti x integral equations spectral properties of totally positive kernels and matrices a pinkus index

a finite element method for the simulation of strong and weak discontinuities in solid mechanics in this paper we introduce and analyze a finite element method for elasticity problems with interfaces the method allows for discontinuities internal to the elements in the approximation across the interface we propose a general approach that can handle both perfectly and imperfectly bonded interfaces without modifications of the code for the case of linear elasticity we show that optimal order of convergence holds without restrictions on the location of the interface relative to the mesh we present numerical examples for the linear case as well as for contact and crack propagation model problems

on some techniques for approximating boundary conditions in the finite element method we discuss the stabilization of finite element methods in which essential boundary conditions are approximated by babuska s method of lagrange multipliers and we show that there is a close connection with this technique and a classical method by nitsche

contact with friction in multi material arbitrary lagrangian eulerian formulations using x fem a contact method with friction for the multi dimensional lagrangian step in multi material arbitrary lagrangian eulerian ale formulations is presented in our previous research the extended finite element method x fem was used to create independent fields i e velocity strain rate force mass etc for each material in the problem to model contact without friction the research presented here includes the extension to friction and improvements to the accuracy and robustness of our previous study the accelerations of the multi material nodes are obtained by coupling the material force and mass fields as a function of the prescribed contact similarly the velocities of the multi material nodes are recalculated using the conservation of momentum when the prescribed contact requires it the coupling procedures impose the same nodal velocity on the coupled materials in the direction normal to their interface during the time step update as a result the overlap of materials is prevented and unwanted separation does not occur three different types of contacts are treated perfectly bonded frictionless slip and slip with friction example impact problems are solved and the numerical solutions are presented copyright 2008 john wiley sons ltd

an eulerian finite element formulation for modelling stationary finite strain elastic deformation processes a new finite element formulation aimed at the modelling of stationary finite strain elastic deformation processes is presented the new formulation is based on an eulerian description of motion and on the transport of the deformation gradient tensor copyright 2004 john wiley sons ltd

an extended finite element formulation for contact in multi material arbitrary lagrangian eulerian calculations a contact method for the multi dimensional lagrangian step in multi material arbitrary lagrangian eulerian ale formulations is presented the extended finite element method x fem is used to create independent fields i e velocity strain rate force mass etc for each material in the problem the accelerations of the multi material nodes are obtained by coupling the material force and mass fields as a function of the prescribed contact similarly the velocities of the multi material nodes are recalculated using the conservation of momentum when the prescribed contact requires it the coupling procedures impose the same nodal velocity on the coupled materials in the direction normal to their interface during the time step update as a result the overlap of materials is prevented and unwanted separation does not occur three different contacts are treated perfectly bonded frictionless slip and slip with friction example impact problems are solved and the numerical solutions are presented

generalized bezier curves and surfaces based on lupas q analogue of bernstein operator in this paper a new generalization of bezier curves with one shape parameter is constructed it is based on the lupas q analogue of bernstein operator which is the first generalized bernstein operator based on the q calculus the new curves have some properties similar to classical bezier curves moreover we establish degree evaluation and de casteljau algorithms for the generalization furthermore we construct the corresponding tensor product surfaces over the rectangular domain and study the properties of the surfaces as well as the degree evaluation and de casteljau algorithms compared with q bezier curves and surfaces based on phillips q bernstein polynomials our generalizations show more flexibility in choosing the value of q and superiority in shape control of curves and surfaces the shape parameters provide more convenience for the curve and surface modeling

current algebra aks theorem and new super evolution equations the authors have deduced a new class of integrable super evolution equations by using a supersymmetric version of the aks adler kostant sym theorem in conjunction with the homogeneous space reduction technique of marshall 1988 and current algebra of samenov tian shansky 1983

extending fundamental formulas from classical b splines to quantum b splines we derive a collection of fundamental formulas for quantum b splines analogous to known fundamental formulas for classical b splines starting from known recursive formulas for evaluation and quantum differentiation along with quantum analogues of the marsden identity we derive quantum analogues of the de boor fix formula for the dual functionals explicit formulas for the quantum b splines in terms of divided differences of truncated power functions formulas for computing divided differences of arbitrary functions by quantum integrating certain quantum derivatives of these functions with respect to the quantum b splines closed formulas for the quantum integral of the quantum b splines over their support and finally a 1 q convolution formula for uniform q b splines

optimal bases for a class of mixed spaces and their associated spline spaces we consider spaces of the form span where the functions u i i 1 4 are algebraic polynomials or trigonometric or hyperbolic functions we find intervals 0 a where we can guarantee that the spaces possess normalized totally positive bases and so shape preserving representations we construct their normalized b bases with optimal shape preserving properties we also present a unified approach to deal with the associated spline spaces a recursive procedure to obtain bases with optimal shape preserving and stability properties is presented

an unfitted finite element method based on nitsche s method for elliptic interface problems in this paper we propose a method for the finite element solution of elliptic interface problem using an approach due to nitsche the method allows for discontinuities internal to the elements in the approximation across the interface we show that optimal order of convergence holds without restrictions on the location of the interface relative to the mesh further we derive a posteriori error estimates for the purpose of controlling functionals of the error and present some numerical examples

a bubble stabilized finite element method for dirichlet constraints on embedded interfaces we examine a bubble stabilized finite element method for enforcing dirichlet constraints on embedded interfaces by embedded we refer to problems of general interest wherein the geometry of the interface is assumed independent of some underlying bulk mesh as such the robust imposition of dirichlet constraints using a lagrange multiplier field is not trivial to focus issues we consider a simple one sided problem that is representative of a wide class of evolving interface problems the bulk field is decomposed into coarse and fine scales giving rise to coarse scale and fine scale one sided sub problems the fine scale solution is approximated with bubble functions permitting static condensation and giving rise to a stabilized form bearing strong analogy with a classical method importantly the method is simple to implement readily extends to multiple dimensions obviates the need to specify any free stabilization parameters and can lead to a symmetric positive definite system of equations the performance of the method is then examined through several numerical examples the accuracy of the lagrange multiplier is compared to results obtained using a local version of the domain integral method the variational multiscale approach proposed herein is shown to both stabilize the lagrange multiplier and improve the accuracy of the post processed fluxes copyright 2006 john wiley sons ltd

robust imposition of dirichlet boundary conditions on embedded surfaces summary we develop both stable and stabilized methods for imposing dirichlet constraints on embedded three dimensional surfaces in finite elements the stable method makes use of the vital vertex algorithm of bechet et al bechet et al 2009 int j numer meth engng 78 8 931 954 albeit with a modified set of discontinuous basis functions for the lagrange multiplier field on the embedded surface the stabilized method extends the work of dolbow and harari dolbow and harari 2009 int j numer meth engng 78 2 229 252 to three dimensional surfaces algorithmic and implementational details of both methods are provided several three dimensional benchmark problems are studied to compare and contrast the accuracy of the two approaches the results indicate that both methods yield optimal rates of convergence in various quantities of interest with the primary differences being in the surface flux the utility of the domain integral for extracting accurate surface fluxes is demonstrated for both techniques published by john wiley sons ltd copyright c 2000 john wiley sons ltd

discontinuous galerkin finite element methods for interface problems a priori and a posteriori error estimations discontinuous galerkin dg finite element methods were studied by many researchers for second order elliptic partial differential equations and a priori error estimates were established when the solution of the underlying problem is piecewise h 3 2 epsilon smooth with epsilon 0 however elliptic interface problems with intersecting interfaces do not possess such a smoothness in this paper we establish a quasi optimal a priori error estimate for interface problems whose solutions are only h 1 alpha smooth with alpha in 0 1 and hence fill a theoretical gap of the dg method for elliptic problems with low regularity the second part of the paper deals with the design and analysis of robust residual and recovery based a posteriori error estimators theoretically we show that the residual and recovery estimators studied in this paper are robust with respect to the dg norm i e their reliability and efficiency bounds do not depend on the jump provided that the distribution of coefficients is locally quasi monotone

skew symmetric nitsche s formulation in isogeometric analysis dirichlet and symmetry conditions patch coupling and frictionless contact abstract a simple skew symmetric nitsche s formulation is introduced into the framework of isogeometric analysis iga to deal with various problems in small strain elasticity essential boundary conditions symmetry conditions for kirchhoff plates patch coupling in statics and in modal analysis as well as signorini contact conditions for linear boundary or interface conditions the skew symmetric formulation is parameter free for contact conditions it remains stable and accurate for a wide range of the stabilization parameter several numerical tests are performed to illustrate its accuracy stability and convergence performance we investigate particularly the effects introduced by nitsche s coupling including the convergence performance and condition numbers in statics as well as the extra outlier frequencies and corresponding eigenmodes in structural dynamics we present the hertz test the block test and a 3d self contact example showing that the skew symmetric nitsche s formulation is a suitable approach to simulate contact problems in iga

a semi lagrangean time integration approach for extended finite element methods many computational problems incorporate discontinuities that evolve in time the extendend finite element method xfem is able to represent discontinuities sharply on fixed arbitrary meshes but numerical difficulties arise if these discontinuities move in time we point out that this issue is crucial for interface problems with strongly discontinuous fields on fixed grids a method using semi lagrangean techniques is proposed to adequately handle time integration based on finite difference schemes in the context of the xfem the basic idea is to adapt previous numerical solutions to the current interface position by tracking back virtual lagrangean particles to their previous positions where an appropriate solution can be extrapo lated from a smooth field convergence properties of the proposed method in time and space are thoroughly studied for two one dimensional model problems finally the method is applied to the particularly challenging problem of premixed combustion where the discontinuity appears at the flame front separating the burnt from the unburnt gases a two dimensional and a three dimensional expanding flame demonstrates that the method is sufficiently accurate to retain the properties of the overall nitsche type formulation for interface problems with embedded strong discontinuities

the extended generalized finite element method an overview of the method and its applications an overview of the extended generalized finite element method gefm xfem with emphasis on methodological issues is presented this method enables the accurate approximation of solutions that involve jumps kinks singularities and other locally non smooth features within elements this is achieved by enriching the polynomial approximation space of the classical finite element method the gefm xfem has shown its potential in a variety of applications that involve non smooth solutions near interfaces among them are the simulation of cracks shear bands dislocations solidification and multi field problems copyright 2010 john wiley sons ltd

stable imposition of stiff constraints in explicit dynamics for embedded finite element methods summary we investigate various strategies to enforce the kinematics at an embedded interface for transient problems within the extended finite element method in particular we focus on explicit time integration of the semi discrete equations of motion and extend both dual and primal variational frameworks for constraint enforcement to a transient regime we reiterate the incompatibility of the dual formulation with purely explicit time integration and the severe restrictions placed by the courant friedrichs levy condition on primal formulations we propose an alternate consistent formulation for the primal method and derive an estimate for the stabilization parameter which is more amenable in an explicit dynamics framework importantly the use of the new estimate circumvents the need for any tolerances as an interface approaches an element boundary we also show that with interfacial constraints existing mass lumping schemes can lead to prohibitively small critical time steps accordingly we propose a mass lumping procedure which provides a more favorable estimate these techniques are then demonstrated on several benchmark numerical examples where we compare and contrast the accuracy of the primal methods against the dual methods in enforcing the constraints copyright 2012 john wiley sons ltd

eulerian formulation using stabilized finite element method for large deformation solid dynamics this paper describes an eulerian formulation for large deformation solid dynamics in the present eulerian approach an advective equation is solved using the stream upwind petrov galerkin finite element method the eulerian finite element method is applied to path dependent solid analyses such as impact bar and ductile necking problems these computational results using the eulerian finite element method are compared with the results obtained from using the lagrangian finite element method and an eulerian formulation based on a finite difference method copyright 2007 john wiley sons ltd

convexity of rational curves and total positivity we analyze convexity preserving properties of curves from a geometric point of view we also characterize totally positive systems of functions in terms of geometric convexity preserving properties of the rational curves rational bezier and nonuniform rational b spline curves are included in this setting

on strategies for enforcing interfacial constraints and evaluating jump conditions with the extended finite element method we consider a problem stemming from recent models of phase transitions in stimulus responsive hydrogels wherein a sharp interface separates swelled and collapsed phases extended finite element methods that approximate the local solution with an enriched basis such that the mesh need not explicitly fit the interface geometry are emphasized attention is focused on the weak enforcement of the normal configurational force balance and various options for evaluating the jump in the normal component of the solute flux at the interface we show that as the reciprocal interfacial mobility vanishes it plays the role of a penalty parameter enforcing a pure dirichlet constraint eventually triggering oscillations in the interfacial velocity we also examine alternative formulations employing a lagrange multiplier to enforce this constraint it is shown that the most convenient choice of basis for the lagrange multiplier results in oscillations in the multiplier field and a decrease in accuracy and rate of convergence in local error norms suggesting a lack of stability in the discrete formulation under such conditions neither the direct evaluation of the gradient of the approximation at the phase interface nor the interpretation of the lagrange multiplier field provide a robust means to obtain the jump in the normal component of solute flux fortunately the adaptation and use of local domain integral methodologies considerably improves the flux evaluations several example problems are presented to compare and contrast the various techniques and methods

residual free bubbles for embedded dirichlet problems abstract we examine a stabilized method employing residual free bubbles for enforcing dirichlet constraints on embedded finite element interfaces in particular we focus attention on problems where the underlying mesh is not explicitly fitted to the geometry of the interface the residual free bubble problem is derived for a simple case and extensions are discussed we show that under certain conditions stabilization only requires knowledge of the residual free bubble on the interface we then examine methods to approximate the residual free bubble a series of benchmark tests are used to demonstrate the accuracy and robustness of the method comparisons are made to nitsche s method employing an eigenvalue estimate for the global stability parameter particular emphasis is placed on the accuracy of flux evaluations on the interface finally we employ the method to simulate an evolving interface problem motivated by resin transfer molding

a fully eulerian formulation for fluid structure interaction problems in this work we present a fully eulerian framework for fluid structure interaction fsi problems coupling the incompressible navier stokes equations with a hyperelastic solid the fully eulerian framework is a monolithic implicit variational formulation for the coupled problem in contrast to the well established arbitrary lagrangian eulerian ale coordinates the fully eulerian framework formulates both subproblems fluid and solid in eulerian coordinates this concept circumvents various difficulties connected to ale coordinates since no artificial domain mapping is used the formulation is an interface capturing method and uses an extension of the solid s deformation the initial point set to detect the interface location by construction very large deformation as well as topology changes like contact of the solid to the domain boundary or other solid parts are possible

total positivity and neville elimination abstract neville elimination is described in terms of schur complements of matrices and used to improve some well known characterizations of totally positive and strictly totally positive matrices

a grid based particle method for moving interface problems we propose a novel algorithm for modeling interface motions the interface is represented and is tracked using quasi uniform meshless particles these particles are sampled according to an underlying grid such that each particle is associated to a grid point which is in the neighborhood of the interface the underlying grid provides an eulerian reference and local sampling rate for particles on the interface it also renders neighborhood information among the meshless particles for local reconstruction of the interface the resulting algorithm which is based on lagrangian tracking using meshless particles with eulerian reference grid can naturally handle control topological changes moreover adaptive sampling of the interface can be achieved easily through local grid refinement with simple quad oct tree data structure extensive numerical examples are presented to demonstrate the capability of our new algorithm

imposing essential boundary conditions in mesh free methods imposing essential boundary conditions is a key issue in mesh free methods the mesh free interpolation does not verify the kronecker delta property and therefore the imposition of prescribed values is not as straightforward as for the finite element method the aim of this paper is to present a general overview on the existing techniques to enforce essential boundary conditions in galerkin based mesh free methods special attention is paid to the mesh free coupling with finite elements for the imposition of prescribed values and to methods based on a modification of the galerkin weak form particular examples are used to analyze and compare their performance in different situations

on methods for stabilizing constraints over enriched interfaces in elasticity summary enriched finite element approaches such as the extended finite element method provide a framework for constructing approximations to solutions of non smooth problems internal features such as boundaries are represented in such methods by using discontinuous enrichment of the standard finite element basis within such frameworks however imposition of interface constraints and or constitutive relations can cause unexpected difficulties depending upon how relevant fields are interpolated on un gridded interfaces this work address the stabilized treatment of constraints in an enriched finite element context both the lagrange multiplier and penalty enforcement of tied constraints for an arbitrary boundary represented in an enriched finite element context can lead to instabilities and artificial oscillations in the traction fields we demonstrate two alternative variational methods that can be used to enforce the constraints in a stable manner in a bubble stabilized approach fine scale degrees of freedom are added over elements supporting the interface the variational form can be shown to have a similar form to a second approach we consider nitsche s method with the exception that the stabilization terms follow directly from the bubble functions in this work we examine alternative variational methods for enforcing a tied constraint on an enriched interface in the context of two dimensional elasticity we examine several benchmark problems in elasticity and show that only nitsche s method and the bubble stabilization approach produce stable traction fields over internal boundaries we also demonstrate a novel difference between the penalty method and nitsche s method in that the latter passes the patch test exactly regardless of the stabilization parameter s magnitude results for more complicated geometries and triple interface junctions are also presented copyright 2008 john wiley sons ltd

a weighted nitsche stabilized method for small sliding contact on frictional surfaces abstract we propose a weighted nitsche framework for small sliding frictional contact problems on three dimensional interfaces the proposed method inherits the advantages of both augmented lagrange multiplier and penalty methods while also addressing their shortcomings algorithmic details of the traction update and consistent linearization in the presence of nitsche terms are provided several benchmark numerical experiments are conducted and the results are compared with existing studies the results are encouraging and indicate accurate satisfaction of the non interpenetration constraint stable tractions and asymptotic quadratic convergence of the newton raphson method

new nonlinear evolution equations related to some superloop algebras constructs hierarchies of nonlinear evolution equations on the superloop algebras related to the super lie algebras b 0 1 and sl 2 1 the hierarchy associated with sl 2 1 contains a generalisation of the super nonlinear schrodinger equation with a cubic non derivative interaction term in complex anticommuting fields in all cases the author gives a recursive procedure for constructing involutive integrals of local polynomials

a robust preconditioning technique for the extended finite element method the extended finite element method enhances the approximation properties of the finite element space by using additional enrichment functions but the resulting stiffness matrices can become ill conditioned in that case iterative solvers need a large number of iterations to obtain an acceptable solution in this paper a procedure is described to obtain stiffness matrices whose condition number is close to the one of the finite element matrices without any enrichments a domain decomposition is employed and the algorithm is very well suited for parallel computations the method was tested in numerical experiments to show its effectiveness the experiments have been conducted for structures containing cracks and material interfaces we show that the corresponding enrichments can result in arbitrarily ill conditioned matrices the method proposed here however provides well conditioned matrices and can be applied to any sort of enrichment the complexity of this approach and its relation to the domain decomposition is discussed computation times have been measured for a structure containing multiple cracks for this structure the computation times could be decreased by a factor of 2

a finite strain eulerian formulation for compressible and nearly incompressible hyperelasticity using high order b spline finite elements summary we present a numerical formulation aimed at modeling the nonlinear response of elastic materials using large deformation continuum mechanics in three dimensions this finite element formulation is based on the eulerian description of motion and the transport of the deformation gradient when modeling a nearly incompressible solid the transport of the deformation gradient is decomposed into its isochoric part and the jacobian determinant as independent fields a homogeneous isotropic hyperelastic solid is assumed and b splines based finite elements are used for the spatial discretization a variational multiscale residual based approach is employed to stabilize the transport equations the performance of the scheme is explored for both compressible and nearly incompressible applications the numerical results are in good agreement with theory illustrating the viability of the computational scheme copyright 2011 john wiley sons ltd

some remarks on weighted lupas q bezier curves abstract recently han et al 2016 introduced weighted lupas q bezier curves as a novel generalization of the classical rational bezier curves in this note we make some remarks on their article in particular we revisit the rational quadratic case for an illustrative exposition

the mortar finite element method with lagrange multipliers the present paper deals with a variant of a non conforming domain decomposition technique the mortar finite element method in the opposition to the original method this variant is never conforming because of the relaxation of the matching constraints at the vertices and the edges in 3d of subdomains it is shown that written under primal hybrid formulation the approximation problem issued from a discretization of a second order elliptic equation in 2d is nonetheless well posed and provides a discrete solution that satisfies optimal error estimates with respect to natural norms finally the parallelization advantages consequence of this variant are also addressed

the shifted boundary method for embedded domain computations part ii linear advection diffusion and incompressible navier stokes equations abstract we propose a new embedded finite element method for the linear advection diffusion equation and the laminar and turbulent incompressible navier stokes equations the proposed method belongs to the class of surrogate approximate boundary algorithms and is based on the idea of shifting the location where boundary conditions are applied from the true to a surrogate boundary accordingly boundary conditions enforced weakly are appropriately modified to preserve optimal error convergence rates we include the full analysis of stability and convergence of the method in the linear advection diffusion equation and a battery of tests for the case of laminar and turbulent incompressible navier stokes equations we also discuss the conservation properties of the proposed method in all cases

improved implementation and robustness study of the x fem for stress analysis around cracks numerical crack propagation schemes were augmented in an elegant manner by the x fem method the use of special tip enrichment functions as well as a discontinuous function along the sides of the crack allows one to do a complete crack analysis virtually without modifying the underlying mesh which is of industrial interest especially when a numerical model for crack propagation is desired this paper improves the implementation of the x fem method for stress analysis around cracks in three ways first the enrichment strategy is revisited the conventional approach uses a topological enrichment only the elements touching the front are enriched we suggest a geometrical enrichment in which a given domain size is enriched the improvements obtained with this enrichment are discussed second the conditioning of the x fem both for topological and geometrical enrichments is studied a preconditioner is introduced so that off the shelf iterative solver packages can be used and perform as well on x fem matrices as on standard fem matrices the preconditioner uses a local nodal cholesky based decomposition third the numerical integration scheme to build the x fem stiffness matrix is dramatically improved for tip enrichment functions by the use of an ad hoc integration scheme a 2d benchmark problem is designed to show the improvements and the robustness

nitsche s method for general boundary conditions we introduce a method for treating general boundary conditions in the finite element method generalizing an approach due to nitsche 1971 for approximating dirichlet boundary conditions we use poisson s equations as a model problem and prove a priori and a posteriori error estimates the method is also compared with the traditional galerkin method the theoretical results are verified numerically

cutfem discretizing geometry and partial differential equations we discuss recent advances on robust unfitted finite element methods on cut meshes these methods are designed to facilitate computations on complex geometries obtained for example from computer

spectral element method for three dimensional elliptic problems with smooth interfaces abstract in this paper we propose a least squares spectral element method for three dimensional elliptic interface problems the differentiability estimates and the main stability theorem using non conforming spectral element functions are proven the proposed method is free from any kind of first order reformulation a suitable preconditioner is constructed with help of the regularity estimate and proposed stability estimates which is used to control the condition number we show that these preconditioners are spectrally equivalent to the quadratic forms by which we approximate them we obtain the error estimates which show the exponential accuracy of the method numerical results are obtained for both straight and curved interfaces to show the efficiency of the proposed method

on the backlund transformation and hamiltonian properties of superevaluation equations backlund and darboux backlund transformations are deduced for the superevaluation equations recently deduced by kupershmidt b kuperschmidt a super k dv equation an integrable system preprint utsi tullahoma 1984 phys lett a 102 213 1983 j phys a 17 l 863 1984 and gurses h gurses and o oguz phys lett a 108 437 1985 by a extension of the technique of bpt m boiti f pempinelli and g z tu nuovo cimento b 79 231 1984 to anticommuting variables the bi hamiltonian structure and hence the form of the recursion operator for the lie backlund symmetry for such equations are deduced incidentally some explicit forms of the lie backlund symmetry are also deduced

on the lupas q transform the lupasq transform emerges in the study of the limit q lupas operator the latter comes out naturally as a limit for a sequence of the lupas q analogues of the bernstein operator lately it has been studied by several authors from different perspectives in mathematical analysis and approximation theory this operator is closely related to the q deformed poisson probability distribution which is used widely in the q boson operator calculus given email protected 0 1 email protected c 0 1 the q lupas transform of f is defined by l qf z 1 z q email protected email protected k 0 f 1 q k q k k 1 2 q q kz k in this paper we study some analytic properties of l qf z in particular we examine the conditions under which l qf can either be an entire function or a rational one

the shifted boundary method for hyperbolic systems embedded domain computations of linear waves and shallow water flows abstract we propose a new computational approach for embedded boundary simulations of hyperbolic systems and in particular the linear wave equations and the nonlinear shallow water equations the proposed approach belongs to the class of surrogate approximate boundary algorithms and is based on the idea of shifting the location where boundary conditions are applied from the true to a surrogate boundary accordingly boundary conditions enforced weakly are appropriately modified to preserve optimal error convergence rates this framework is applied here in the setting of a stabilized finite element method even though other spatial discretization techniques could have been employed accuracy stability and robustness of the proposed method are tested by means of an extensive set of computational experiments for the acoustic wave propagation equations and shallow water equations comparisons with standard weak boundary conditions imposed on body fitted grids which conform to the geometry of the computational domain boundaries are also presented

normalized b basis of the space of trigonometric polynomials and curve design a normalized b basis of the space of trigonometric polynomials of degree n is presented some interesting properties of the basis functions are given based on the basis symmetric trigonometric polynomial curves like bezier curves are constructed the trigonometric polynomial curves present the shape of their control polygons well thus the theoretics and methods are proposed for curve representation of the trigonometric polynomial space by adding additional control points the given graph examples show that the trigonometric polynomial curves are nearer their control polygons than the bezier curves for the same parametric variable and the same degree

q bernstein polynomials and bezier curves we define q bernstein polynomials which generalize the classical bernstein polynomials and show that the difference of two consecutive q bernstein polynomials of a function f can be expressed in terms of second order divided differences of f it is also shown that the approximation to a convex function by its q bernstein polynomials is one sided a parametric curve is represented using a generalized bernstein basis and the concept of total positivity is applied to investigate the shape properties of the curve we study the nature of degree elevation and degree reduction for this basis and show that degree elevation is variation diminishing as for the classical bernstein basis

weighted lupas q bezier curves this paper is concerned with a new generalization of rational bernstein bezier curves involving q integers as shape parameters a one parameter family of rational bernstein bezier curves weighted lupas q bezier curves is constructed based on a set of lupas q analogue of bernstein functions which is proved to be a normalized totally positive basis the generalized rational bezier curve is investigated from a geometric point of view the investigation provides the geometric meaning of the weights and the representation for conic sections we also obtain degree evaluation and de casteljau algorithms by means of homogeneous coordinates numerical examples show that weighted lupas q bezier curves have more modeling flexibility than classical rational bernstein bezier curves and lupas q bezier curves and meanwhile they provide better approximations to the control polygon than rational phillips q bezier curves

an approach to generate superextensions of integrable systems an algorithm to generate integrable systems is extended to the super case some new examples of superextensions of integrable systems are illustrated we also generalize the trace identity due to tu to the super case and use it to establish hamiltonian structures of superextensions of integrable systems under consideration

a powerful approach to generate new integrable systems in this paper an effective algorithm to generate integrable systems is given as a result many new integrable equations are derived in a systematic way

analysis of backward euler extended finite element discretization of parabolic problems with moving interfaces abstract numerical experiments performed by several authors repeatedly confirmed that the extended finite element method combined with classical time stepping schemes such as backward euler or crank nicolson provides optimal rates of convergence when applied to problems with moving interfaces in spite of the rapid spread of this discretization approach a complete space and time error analysis is not available yet the objective of this work is to contribute to closing the gap between observations and rigorous analysis more precisely we preform a thorough space and time error analysis of different variants of the extended finite element method combined with backward euler time advancing scheme when applied to problems with moving interfaces our results prove that the discretization error measured in the natural norms decreases with the expected rate when space and time discretization steps are refined we also show that the space and time approximation properties are fully decoupled as a consequence high order approximation schemes could be developed and analyzed within the theoretical framework proposed here numerical experiments are finally addressed for the verification of the algorithms besides the interest of providing rigorous error bounds we believe that a general theoretical framework is extremely helpful as a guide for further developing and refining extended finite element methods for free interface problems

a mixed eulerian lagrangian displacement model for large deformation analysis in solid mechanics abstract a new mixed kinematic description for large deformation solid mechanics is presented the novel aspect of the new model is that both the initial and current configurations are treated as unknown large deformation statements of the principles of virtual work and stationary potential energy are formulated with the new description and presented in incremental form suitable for iterative finite element solutions the well known total lagrangian and updated lagrangian descriptions are shown to be special cases of the mixed formulation the new description is advantageous in situations where it is difficult to select an appropriate material discretization prior to analysis applications to large deformation frictional contact and fracture mechanics are presented to demonstrate the utility of the new kinematic description

bernstein schurer kantorovich operators based on q integers in this paper we introduce a new kantorovich type generalization of the q bernstein schurer operators defined in muraru 2011 first we give the basic convergence theorem and then obtain the local direct results for these operators estimating the rate of convergence by using the modulus of smoothness and the lipschitz type maximal function respectively we also obtain a voronovskaja type theorem and investigate the statistical approximation properties of these operators with the help of a korovkin type statistical approximation theorem given in duman 2008

a stable lagrange multiplier space for stiff interface conditions within the extended finite element method this paper introduces a new algorithm to define a stable lagrange multiplier space to impose stiff interface conditions within the context of the extended finite element method in contrast to earlier approaches we do not work with an interior penalty formulation as e g for nitsche techniques but impose the constraints weakly in terms of lagrange multipliers roughly speaking a stable and optimal discrete lagrange multiplier space has to satisfy two criteria a best approximation property and a uniform inf sup condition owing to the fact that the interface does not match the edges of the mesh the choice of a good discrete lagrange multiplier space is not trivial here we propose a new algorithm for the local construction of the lagrange multiplier space and show that a uniform inf sup condition is satisfied a counterexample is also presented i e the inf sup constant depends on the mesh size and degenerates as it tends to zero numerical results in two dimensional confirm the theoretical ones copyright

stable generalized finite element method sgfem abstract the generalized finite element method gfem is a partition of unity method pum where the trial space of standard finite element method fem is augmented with non polynomial shape functions with compact support these shape functions which are also known as the enrichments mimic the local behavior of the unknown solution of the underlying variational problem gfem has been successfully used to solve a variety of problems with complicated features and microstructure however the stiffness matrix of gfem is badly conditioned much worse compared to the standard fem and there could be a severe loss of accuracy in the computed solution of the associated linear system in this paper we address this issue and propose a modification of the gfem referred to as the stable gfem sgfem we show that sgfem retains the excellent convergence properties of gfem does not require a ramp function in the presence of blending elements and the conditioning of the associated stiffness matrix is not worse than that of the standard fem moreover sgfem is very robust with respect to the parameters of the enrichments we show these features of sgfem on several examples

imposing dirichlet boundary conditions with nitsche s method and spline based finite elements a key challenge while employing non interpolatory basis functions in finite element methods is the robust imposition of dirichlet boundary conditions the current work studies the weak enforcement of such conditions for b spline basis functions with application to both second and fourth order problems this is achieved using concepts borrowed from nitsche s method which is a stabilized method for imposing constraints on surfaces conditions for the stability of the system of equations are derived for each class of problem stability parameters in the nitsche weak form are then evaluated by solving a local generalized eigenvalue problem at the dirichlet boundary the approach is designed to work equally well when the grid used to build the splines conforms to the physical boundary of interest as well as to the more general case when it does not through several numerical examples the approach is shown to yield optimal rates of convergence copyright 2010 john wiley sons ltd

a coupled eulerian lagrangian extended finite element formulation for simulating large deformations in hyperelastic media with moving free boundaries we present a coupled eulerian lagrangian cel formulation aimed at modeling the moving interface of hyperelastic materials undergoing large to extreme deformations this formulation is based on an eulerian description of kinematics of deformable bodies together with an updated lagrangian formulation for the transport of the deformation gradient tensor the extended finite element method xfem is used to discretize the mechanical equilibrium and deformation gradient transport equations in a two phase domain a mixed interpolation scheme biquadratic for the velocity and bilinear for the deformation gradient is adopted to improve the accuracy of the numerical formulation the interface describing the deformed shape of the body is represented by the level set function and is evolved using the grid based particle method the performance of the scheme is explored in two dimensions in the compressible regime for an adequate spatial and temporal discretization our numerical results are in good agreement with theory and with numerical results from the traditional lagrangian formulation in abaqus the advantage of the proposed formulation is that material motion is not coupled with that of the mesh this eliminates the issues of mesh distortion and the need for remeshing associated with lagrangian formulations when bodies undergo very large distortions it is therefore well adapted to describe the motion of complex fluids and soft matter whose physical properties are intermediate between conventional liquids and solids