high order methods for low reynolds number flows around moving obstacles based on universal meshes summary we propose a family of methods for simulating two dimensional incompressible low reynolds number flow around a moving obstacle whose motion is prescribed the methods make use of a universal mesh a fixed background mesh that adapts to the geometry of the immersed obstacle at all times by adjusting a few elements in the neighborhood of the obstacle s boundary the resulting mesh provides a conforming triangulation of the fluid domain over which discretizations of any desired order of accuracy in space and time can be constructed using standard finite element spaces together with off the shelf time integrators we demonstrate the approach by using taylor hood elements to approximate the fluid velocity and pressure to integrate in time we consider implicit runge kutta schemes as well as a fractional step scheme we illustrate the methods and study their convergence numerically via examples that involve flow around obstacles that undergo prescribed deformations copyright 2015 john wiley sons ltd

a three field stabilized finite element method for fluid structure interaction elastic solid and rigid body limit summary we propose a full eulerian framework for solving fluid structure interaction fsi problems based on a unified formulation in which the fsis are modelled by introducing an extra stress in the momentum equation the obtained three field velocity pressure and stress system is solved using a stabilized finite element method the key feature of this unified formulation is the ability to describe different kind of interactions between the fluid and the structure which can be either elastic or a perfect rigid body without the need of treating this last case via penalization the level set method combined with a dynamic anisotropic mesh adaptation is used to track the fluid solid interface copyright 2015 john wiley sons ltd

multiscale phenomena green s functions the dirichlet to neumann formulation subgrid scale models bubbles and the origins of stabilized methods abstract an approach is developed for deriving variational methods capable of representing multiscale phenomena the ideas are first illustrated on the exterior problem for the helmholtz equation this leads to the well known dirichlet to neumann formulation next a class of subgrid scale models is developed and the relationships to bubble function methods and stabilized methods are established it is shown that both the latter methods are approximate subgrid scale models the identification for stabilized methods leads to an analytical formula for the intrinsic time scale whose origins have been a mystery heretofore

the fixed mesh ale approach for the numerical approximation of flows in moving domains in this paper we propose a method to approximate flow problems in moving domains using always a given grid for the spatial discretization and therefore the formulation to be presented falls within the category of fixed grid methods even though the imposition of boundary conditions is a key ingredient that is very often used to classify the fixed grid method our approach can be applied together with any technique to impose approximately boundary conditions although we also describe the one we actually favor our main concern is to properly account for the advection of information as the domain boundary evolves to achieve this we use an arbitrary lagrangian eulerian framework the distinctive feature being that at each time step results are projected onto a fixed background mesh that is where the problem is actually solved

improved fractional step method for simulating fluid structure interaction using the pfem summary numerical difficulties are present in the particle finite element method even though it has been shown to be a powerful and effective approach to simulating fluid structure interaction to overcome problems of mass loss on the free surface and the added mass effect an improved fractional step method fsm that handles added mass terms in a mathematically exact way is developed a further benefit is that no assumptions regarding the structural response are made in handling added mass terms thus it is straightforward to incorporate material nonlinearity in fluid structure interaction fsi under this approach patch tests and comparisons with experimental data are presented in order to verify and validate the improved fsm for fsi applications the computational cost of this approach is shown to be negligible compared with the other aspects of the fsm particularly when the size of the structure and the fluid structure interface is small relative to the volume of fluid copyright 2014 john wiley sons ltd

a complete method for the computations of mathieu characteristic numbers of integer orders a complete method and expressions for computing mathieu characteristic numbers mcns of integer orders are given the paper reviews available algorithms and their shortcomings and then suitable methods are selected on the basis of simplicity and accuracy for small orders n 4 new approximate formulas for mcns are developed these have much larger range than available expressions a new normalization technique for the mcns is proposed with the aid of this normalization a new approximate recurrence relation for mcns is obtained using the same normalization general formulas are obtained for regional limits where the approximate and asymptotic expressions can be used to compute mcns with these expressions the algorithms for computing mcns can handle very high orders up to n 1000 and beyond

adaptively weighted numerical integration in the finite cell method abstract with adaptively weighted aw numerical integration for a given set of quadrature nodes order and domain of integration the quadrature weights are obtained by solving a system of suitable moment fitting equations in least square sense the moments in the moment equations are approximated over a simplified domain that is homeomorphic to the original domain and then are corrected for the deviation from the original domain using shape sensitivity analysis in this paper we demonstrate the application of aw integration scheme in the context of the finite cell method which must perform numerical integration over arbitrary domains without meshing the standard integration technique employed in fcm is the characteristic function method that converts the continuous integrand over a complex domain into a discontinuous integrand over a simple box domain then well known integrand adaptivity techniques are employed to integrate the resulting discontinuous integrand over the box domain although this method is simple to implement it becomes computationally very expensive for realistic complex 3d domains such as sculptures bones and engines in contrast in aw scheme the quadrature weights directly adapt to the complex geometric domain without the need to making the integrand discontinuous leading to superior computational properties in this paper we demonstrate the computational efficiency of aw over the characteristic function method as it requires fewer subdivisions and less time to achieve a given accuracy in both two and three dimensions in addition aw offers a number of advantages including flexibility in the choice of quadrature points and basis functions

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

an immersed boundary method with formal second order accuracy and reduced numerical viscosity a formally second order accurate immersed boundary method is presented and tested in this paper we apply this new scheme to simulate the flow past a circular cylinder and study the effect of numerical viscosity on the accuracy of the computation by comparing the numerical results with those of a first order method the numerical evidence shows that the new scheme has less numerical viscosity and is therefore a better choice for the simulation of high reynolds number flows with immersed boundaries

error analysis for the computation of zeros of regular coulomb wave function and its first derivative in 1975 one of the coauthors ikebe showed that the problem of computing the zeros of the regular coulomb wave functions and their derivatives may be reformulated as the eigenvalue problem for infinite matrices approximation by truncation is justified but no error estimates are given there the class of eigenvalue problems studied there turns out to be subsumed in a more general problem studied by ikebe et al in 1993 where an extremely accurate asymptotic error estimate is shown in this paper we apply this error formula to the former case to obtain error formulas in a closed explicit form

universal meshes a method for triangulating planar curved domains immersed in nonconforming meshes summary we introduce a new method to triangulate planar curved domains that transforms a specific collection of triangles in a background mesh to conform to the boundary in the process no new vertices are introduced and connectivities of triangles are left unaltered the method relies on a novel way of parameterizing an immersed boundary over a collection of nearby edges with its closest point projection to guarantee its robustness we require that the domain be c2 regular the background mesh be sufficiently refined near the boundary and that specific angles in triangles near the boundary be strictly acute the method can render both straight edged and curvilinear triangulations for the immersed domain the latter includes curved triangles that conform exactly to the immersed boundary and ones constructed with isoparametric mappings to interpolate the boundary at select points high order finite elements constructed over these curved triangles achieve optimal accuracy which has customarily proven difficult in numerical schemes that adopt nonconforming meshes aside from serving as a quick and simple tool for meshing planar curved domains with complex shapes the method provides significant advantages for simulating problems with moving boundaries and in numerical schemes that require iterating over the geometry of domains with no conformity requirements the same background mesh can be adopted to triangulate a large family of domains immersed in it including ones realized over several updates during the coarse of simulating problems with moving boundaries we term such a background mesh as a universal mesh for the family of domains it can be used to triangulate universal meshes hence facilitate a framework for finite element calculations over evolving domains while using only fixed background meshes furthermore because the evolving geometry can be approximated with any desired order numerical solutions can be computed with high order accuracy we present demonstrative examples using universal meshes to simulate the interaction of rigid bodies with stokesian fluids copyright 2014 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

approximate imposition of boundary conditions in immersed boundary methods we analyze several possibilities to prescribe boundary conditions in the context of immersed boundary methods as basic approximation technique we consider the finite element method with a mesh that does not match the boundary of the computational domain and therefore dirichlet boundary conditions need to be prescribed in an approximate way as starting variational approach we consider nitsche s methods and we then move to two options that yield non symmetric problems but that turned out to be robust and efficient the essential idea is to use the degrees of freedom of certain nodes of the finite element mesh to minimize the difference between the exact and the approximated boundary condition

weak imposition of essential boundary conditions in the finite element approximation of elliptic problems with non matching meshes this is the accepted version of the following article codina r and baiges j 2015 weak imposition of essential boundary conditions in the finite element approximation of elliptic problems with non matching meshes int j numer meth engng 104 624 654 doi 10 1002 nme 4815 which has been published in final form at http onlinelibrary wiley com doi 10 1002 nme 4815 abstract

immersed interface methods for stokes flow with elastic boundaries or surface tension a second order accurate interface tracking method for the solution of incompressible stokes flow problems with moving interfaces on a uniform cartesian grid is presented the interface may consist of an elastic boundary immersed in the fluid or an interface between two different fluids the interface is represented by a cubic spline along which the singularly supported elastic or surface tension force can be computed the stokes equations are then discretized using the second order accurate finite difference methods for elliptic equations with singular sources developed in our previous paper siam j numer anal 31 1994 pp 1019 1044 the resulting velocities are interpolated to the interface to determine the motion of the interface an implicit quasi newton method is developed that allows reasonable time steps to be used

the eigenvalue problem for infinite compact complex symmetric matrices with application to the numerical computation of complex zeros of j0 z ij1 z and of bessel functions jm z of any real order m abstract consider computing simple eigenvalues of a given compact infinite matrix re garded as operating in the complex hilbert space l 2 by computing the eigenvalues of the truncated finite matrices and taking an obvious limit in this paper we deal with a special case where the given matrix is compact complex and symmetric but not necessarily hermitian two examples of application are studied the first is con cerned with the equation j 0 z ij 1 z 0 appearing in the analysis of the solitary wave runup on a sloping beach and the second with the zeros of the bessel function j m z of any real order m in each case the problem is reformulated as an eigenvalue problem for a compact complex symmetric tridiagonal matrix operator in l 2 whose eigenvalues are all simple a complete error analysis for the numerical solution by truncation is given based on the general theorems proved in this paper where the usefulness of the seldom used generalized rayleigh quotient is demonstrated

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

an immersed interface method for simulating the interaction of a fluid with moving boundaries in the immersed interface method boundaries are represented as singular force in the navier stokes equations which enters a numerical scheme as jump conditions recently we systematically derived all the necessary spatial and temporal jump conditions for simulating incompressible viscous flows subject to moving boundaries in 3d with second order spatial and temporal accuracy near the boundaries sheng xu z jane wang systematic derivation of jump conditions for the immersed interface method in three dimensional flow simulation siam j sci comput 2006 in press in this paper we implement the immersed interface method to incorporate these jump conditions in a 2d numerical scheme we study the accuracy efficiency and robustness of our method by simulating taylor couette flow flow induced by a relaxing balloon flow past single and multiple cylinders and flow around a flapping wing our results show that 1 our code has second order accuracy in the infinity norm for both the velocity and the pressure 2 the addition of an object introduces relatively insignificant computational cost 3 the method is equally effective in computing flow subject to boundaries with prescribed force or boundaries with prescribed motion

computation of multiple eigenvalues of infinite tridiagonal matrices in this paper it is first given as a necessary and sufficient condition that infinite matrices of a certain type have double eigenvalues the computation of such double eigenvalues is enabled by the newton method of two variables the three term recurrence relations obtained from its eigenvalue problem evp subsume the well known relations of a the zeros of j v z b the zeros of zj f v z hj v z c the evp of the mathieu differential equation and d the evp of the spheroidal wave equation the results of experiments are shown for the three cases a c for the computation of their double pairs

a fictitious domain method for dirichlet problem and applications abstract in this article we discuss the solution of the dirichlet problem for a class of elliptic operators by a lagrange multiplier fictitious domain method this approach allows the use of regular grids and therefore of fast specialized solvers for problems on complicated geometries the resulting saddle point system can be solved by an uzawa conjugate gradient algorithm in the case of two dimensional problems a quasi optimal preconditioner has been found by fourier analysis and numerical experiments confirm its nice scaling properties the resulting methodology is applied to a nonlinear time dependent problem namely the flow of a viscous plastic medium in a cylindrical pipe showing the potential of this methodology for some classes of nonlinear problems

the immersed interface method for elliptic equations with discontinuous coefficients and singular sources the authors develop finite difference methods for elliptic equations of the form nabla cdot beta x nabla u x kappa x u x f x in a region omega in one or two space dimension

the eigenvalue problem for infinite complex symmetric tridiagonal matrices with application abstract we consider an infinite complex symmetric not necessarily hermitian tridiagonal matrix t whose diagonal elements diverge to in modulus and whose off diagonal elements are bounded we regard t as a linear operator mapping a maximal domain in the hilbert space l2 into l2 assuming the existence of t 1 we consider the problem of approximating a given simple eigenvalue of t by an eigenvalue n of tn and nth order principal submatrix of t let x x 1 x 2 t be an eigenvector corresponding to assuming xtx 0 and f n 1 x n 1 x n 0 as n we show that there exists a sequence n of tn such that n f n 1 x n 1 1 o 1 x t x 0 where fn 1 represents the n n 1 element of t application to the following problems is included a solve jv z 0 for v given z 0 and b compute the eigenvalues of the mathieu equation fortunately the existence of t 1 need not be verified for these examples since we may show that t i with taken appropriately has an inverse

a topos foundation for theories of physics iii the representation of physical quantities with arrows this paper is the third in a series whose goal is to develop a fundamentally new way of viewing theories of physics our basic contention is that constructing a theory of physics is equivalent to finding a representation in a topos of a certain formal language that is attached to the system in paper ii we studied the topos representations of the propositional language pl s for the case of quantum theory and in the present paper we do the same thing for the more extensive local language l s one of the main achievements is to find a topos representation for self adjoint operators this involves showing that for any physical quantity a there is an arrow o a r where r is the quantity value object for this theory the construction of o a is an extension of the daseinisation of projection operators that was discussed in paper ii the object r is a monoid object only in the topos setsv h op of the theory and to enhance the applicability of the formalism we apply to r a topos a

comparison of some finite element methods for solving the diffusion convection reaction equation abstract in this paper we describe several finite element methods for solving the diffusion convection reaction equation none of them is new although the presentation is non standard in an effort to emphasize the similarities and differences between them in particular it is shown that the classical supg method is very similar to an explicit version of the characteristic galerkin method whereas the taylor galerkin method has a stabilization effect similar to a sub grid scale model which is in turn related to the introduction of bubble functions

efficient implementation of finite element methods on nonmatching and overlapping meshes in three dimensions in recent years a number of finite element methods have been formulated for the solution of partial differential equations on complex geometries based on nonmatching or overlapping meshes examples of such methods are the fictitious domain method the extended finite element method and nitsche s method in all these methods integrals must be computed over cut cells or subsimplices which is challenging to implement especially in three space dimensions in this note we address the main challenges of such an implementation and demonstrate good performance of a fully general code for automatic detection of mesh intersections and integration over cut cells and subsimplices as a canonical example of an overlapping mesh method we consider nitsche s method which we apply to poisson s equation and a linear elastic problem

stabilized finite element approximation of transient incompressible flows using orthogonal subscales abstract in this paper we present a stabilized finite element method to solve the transient navier stokes equations based on the decomposition of the unknowns into resolvable and subgrid scales the latter are approximately accounted for so as to end up with a stable finite element problem which in particular allows to deal with convection dominated flows and the use of equal velocity pressure interpolations three main issues are addressed the first is a method to estimate the behavior of the stabilization parameters based on a fourier analysis of the problem for the subscales secondly the way to deal with transient problems discretized using a finite difference scheme is discussed finally the treatment of the nonlinear term is also analyzed a very important feature of this work is that the subgrid scales are taken as orthogonal to the finite element space in the transient case this simplifies considerably the numerical scheme

the fixed mesh ale approach applied to solid mechanics and fluid structure interaction problems in this paper we propose a method to solve solid mechanics and fluid structure interaction problems using always a fixed background mesh for the spatial discretization the main feature of the method is that it properly accounts for the advection of information as the domain boundary evolves to achieve this we use an arbitrary lagrangian eulerian ale framework the distinctive characteristic being that at each time step results are projected onto a fixed background mesh for solid mechanics problems subject to large strains the fixed mesh fm ale method avoids the element stretching found in fully lagrangian approaches for fsi problems fm ale allows for the use of a single background mesh to solve both the fluid and the structure copyright 2009 john wiley sons ltd

an eulerian lagrangian method for fluid structure interaction based on level sets abstract we describe a method for fluid structure interaction with the interface and free surfaces defined by level sets the fluid is treated by an eulerian mesh whereas the solid or structure by a lagrangian mesh the level set method provides a compact framework for writing the weak form and for developing and computing the discretized equations lagrange multiplier methods and penalty methods are used to couple the fluid and structure results are presented for a few compressible fluid problems representative of high pressure waves

an adaptive fixed mesh ale method for free surface flows in this work we present a fixed mesh ale method for the numerical simulation of free surface flows capable of using an adaptive finite element mesh covering a background domain this mesh is successively refined and unrefined at each time step in order to focus the computational effort on the spatial regions where it is required some of the main ingredients of the formulation are the use of an arbitrary lagrangian eulerian formulation for computing temporal derivatives the use of stabilization terms for stabilizing convection stabilizing the lack of compatibility between velocity and pressure interpolation spaces and stabilizing the ill conditioning introduced by the cuts on the background finite element mesh and the coupling of the algorithm with an adaptive mesh refinement procedure suitable for running on distributed memory environments algorithmic steps for the projection between meshes are presented together with the algebraic fractional step approach used for improving the condition number of the linear systems to be solved the method is tested in several numerical examples the expected convergence rates both in space and time are observed smooth solution fields for both velocity and pressure are obtained as a result of the contribution of the stabilization terms finally a good agreement between the numerical results and the reference experimental data is obtained

an arbitrary lagrangian eulerian finite element method for finite strain plasticity this paper presents a new arbitrary lagrangian eulerian ale finite element formulation for finite strain plasticity in non linear solid mechanics we consider the models of finite strain plasticity defined by the multiplicative decomposition of the deformation gradient in an elastic and a plastic part f fefp with the stresses given by a hyperelastic relation in contrast with more classical ale approaches based on plastic models of the hypoelastic type the ale formulation presented herein considers the direct interpolation of the motion of the material with respect to the reference mesh together with the motion of the spatial mesh with respect to this same reference mesh this aspect is shown to be crucial for a simple treatment of the advection of the plastic internal variables and dynamic variables in fact this advection is carried out exactly through a particle tracking in the reference mesh a calculation that can be accomplished very efficiently with the use of the connectivity graph of the fixed reference mesh a staggered scheme defined by three steps the smoothing the advection and the lagrangian steps leads to an efficient method for the solution of the resulting equations we present several representative numerical simulations that illustrate the performance of the newly proposed methods both quasi static and dynamic conditions are considered in these model examples copyright 2003 john wiley sons ltd