spectral hp least squares finite element formulation for the navier stokes equations we consider the application of least squares finite element models combined with spectral hp methods for the numerical solution of viscous flow problems the paper presents the formulation validation and application of a spectral hp algorithm to the numerical solution of the navier stokes equations governing two and three dimensional stationary incompressible and low speed compressible flows the navier stokes equations are expressed as an equivalent set of first order equations by introducing vorticity or velocity gradients as additional independent variables and the least squares method is used to develop the finite element model high order element expansions are used to construct the discrete model the discrete model thus obtained is linearized by newton s method resulting in a linear system of equations with a symmetric positive definite coefficient matrix that is solved in a fully coupled manner by a preconditioned conjugate gradient method spectral convergence of the l2 least squares functional and l2 error norms is verified using smooth solutions to the two dimensional stationary poisson and incompressible navier stokes equations numerical results for flow over a backward facing step steady flow past a circular cylinder three dimensional lid driven cavity flow and compressible buoyant flow inside a square enclosure are presented to demonstrate the predictive capability and robustness of the proposed formulation

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

an upwinding embedded boundary method for maxwell s equations in media with material interfaces 2d case we propose a new upwinding embedded boundary method to solve time dependent maxwell s equations in media with material interfaces a global second order finite difference method is obtained by combining central difference schemes away from the interfaces and upwinding technique with jump conditions near the interfaces the proposed finite difference method allows time step based on a uniform mesh independent of the locations and shapes of the interfaces moreover the scheme is simple to implement in multidimensional cases numerical tests of wave equations with various types of material interfaces and electromagnetic scattering of 2d cylinders confirm the stability uniform accuracy and ease of implementation of the method

hp dgfem for second order elliptic problems in polyhedra ii exponential convergence the goal of this paper is to establish exponential convergence of hp version interior penalty ip discontinuous galerkin dg finite element methods for the numerical approximation of linear second order elliptic boundary value problems with homogeneous dirichlet boundary conditions and piecewise analytic data in three dimensional polyhedral domains more precisely we shall analyze the convergence of the hp ip dg methods considered in d schotzau c schwab t p wihler siam j numer anal 51 2013 pp 1610 1633 based on axiparallel sigma geometric anisotropic meshes and bm s linear anisotropic polynomial degree distributions

dsc analysis of free edged beams by an iteratively matched boundary method this communication introduces a novel scheme for the treatment of free edge supports in the analysis of beams by using the discrete singular convolution dsc algorithm accommodating free edges has been a challenging issue in the dsc analysis of beams plates and shells an iteratively matched boundary imb method is proposed to overcome the difficulty numerical experiments are carried out to demonstrate that the proposed imb method works very well in dealing with arbitrary combinations of beam edge supports

a galerkin formulation of the mib method for three dimensional elliptic interface problems we develop a three dimensional 3d galerkin formulation of the matched interface and boundary mib method for solving elliptic partial differential equations pdes with discontinuous coefficients i e the elliptic interface problem the present approach builds up two sets of elements respectively on two extended subdomains which both include the interface as a result two sets of elements overlap each other near the interface fictitious solutions are defined on the overlapping part of the elements so that the differentiation operations of the original pdes can be discretized as if there was no interface the extra coefficients of polynomial basis functions which furnish the overlapping elements and solve the fictitious solutions are determined by interface jump conditions consequently the interface jump conditions are rigorously enforced on the interface the present method utilizes cartesian meshes to avoid the mesh generation in conventional finite element methods fems we implement the proposed mib galerkin method with three different elements namely rectangular prism element five tetrahedron element and six tetrahedron element which tile the cartesian mesh without introducing any new node the accuracy stability and robustness of the proposed 3d mib galerkin are extensively validated over three types of elliptic interface problems in the first type interfaces are analytically defined by level set functions these interfaces are relatively simple but admit geometric singularities in the second type interfaces are defined by protein surfaces which are truly arbitrarily complex the last type of interfaces originates from multiprotein complexes such as molecular motors near second order accuracy has been confirmed for all of these problems to our knowledge it is the first time for an fem to show a near second order convergence in solving the poisson equation with realistic protein surfaces additionally the present work offers the first known near second order accurate method for c 1 continuous or h 2 continuous solutions associated with a lipschitz continuous interface in a 3d setting

an observation of quadratic algebra dual family of nonlinear coherent states and their non classical properties in the generalized isotonic oscillator in this paper we construct nonlinear coherent states for the generalized isotonic oscillator and study their non classical properties in detail by transforming the deformed ladder operators suitably which generate the quadratic algebra we obtain heisenberg algebra from the algebra we define two non unitary and one unitary displacement type operators while the action of one of the non unitary type operators reproduces the original nonlinear coherent states the other one fails to produce a new set of nonlinear coherent states dual pair we show that these dual states are not normalizable for the nonlinear coherent states we evaluate the parameter a3 and examine the non classical nature of the states through quadratic and amplitude squared squeezing effect further we derive analytical formula for the p function q function and the wigner function for the nonlinear coherent states all of them confirm the non classicality of the nonlinear coherent states in addition to the above we obtain the

spectral element methods for elliptic problems in nonsmooth domains engineering applications frequently require the numerical solution of elliptic boundary value problems in irregularly shaped domains for smooth problems spectral element methods have proved very successful since they can accommodate fairly complicated geometries while retaining a rapid rate of convergence geometric singularities however often give rise to singular solutions the accuracy of the spectral element methods is then degraded and they offer no apparent advantage over low order finite element methods in many cases however the singular structure of the solution is known and its form may be exploited by the spectral element method among the various ways of doing so through supplementary basis functions eigenfunction expansions and graded meshes the method of auxiliary mapping proves to be particularly effective for certain simple cases the problem is transformed to a coordinate system in which the solution is analytic and exponential convergence is recovered even when this is not possible the singularity is usually much weaker after mapping so that other treatments are more effective in the new coordinate system in this paper we study different ways of treating singularities and in particular the method of auxiliary mapping coupled with the use of supplementary basis functions error estimates are presented explaining why the combined approach is more effective and these estimates are confirmed through a number of numerical experiments for the laplace poisson and helmholtz equations

an immersed interface method for solving anisotropic elliptic boundary value problems in three dimensions we have developed a first order stable cartesian grid discretization that uses only interior grid points for inhomogeneous anisotropic elliptic operators subject to neumann boundary conditions on a bounded nonrectangular geometry in three dimensions for this discretization method a necessary and sufficient condition depending on the mesh size h for the existence of this first order stable scheme at a regular i e interior grid point is found in terms of the anisotropy matrix for this discretization method a way to analyze the existence of a first order stable scheme at an irregular i e boundary grid point is also given the arguments are identical to those for the two dimensional case m a dumett and j p keener a numerical method for solving anisotropic elliptic boundary value problems on an irregular domain in 2d manuscript only the details change unlike in m a dumett and j p keener a numerical method for solving anisotropic elliptic boundary value problems on an irregular domain in 2d manuscript a discussion of dirichlet and robin boundary conditions is also included in particular it is shown that the gerschgorin condition does not impose sign restrictions on irregular grid points stencil coefficients as in the neumann case

analysis of first order system least squares fosls for elliptic problems with discontinuous coefficients part i first order system least squares fosls is a recently developed methodology for solving partial differential equations among its advantages are that the finite element spaces are not restricted by the inf sup condition imposed for example on mixed methods and that the least squares functional itself serves as an appropriate error measure this paper studies the fosls approach for scalar second order elliptic boundary value problems with discontinuous coefficients irregular boundaries and mixed boundary conditions a least squares functional is defined and ellipticity is established in a natural norm of an appropriately scaled least squares bilinear form for some geometries this ellipticity is independent of the size of the jumps in the coefficients the occurrence of singularities at interface corners cross points reentrant corners and irregular boundary points is discussed and a basis of singular functions with local support around singular points is established a companion paper shows that the singular basis functions can be added at little extra cost and lead to optimal performance of standard finite element discretization and multilevel solver techniques also independent of the size of coefficient jumps for some geometries

a polylogarithmic bound for an iterative substructuring method for spectral elements in three dimensions iterative substructuring methods form an important family of domain decomposition algorithms for elliptic finite element problems a p version finite element method based on continuous piecewise q p functions is considered for second order elliptic problems in three dimensions this special method can also be viewed as a conforming spectral element method an iterative method is designed for which the condition number of the relevant operator grows only in proportion to 1 log p 2 this bound is independent of jumps in the coefficient of the elliptic problem across the interfaces between the subregions numerical results are also reported which support the theory

preconditioning for feature selection and regression in high dimensional problems we consider regression problems where the number of predictors greatly exceeds the number of observations we propose a method for variable selection that first estimates the regression function yielding a preconditioned response variable the primary method used for this initial regression is supervised principal components then we apply a standard procedure such as forward stepwise selection or the lasso to the preconditioned response variable in a number of simulated and real data examples this two step procedure outperforms forward stepwise selection or the usual lasso applied directly to the raw outcome we also show that under a certain gaussian latent variable model application of the lasso to the preconditioned response variable is consistent as the number of predictors and observations increases moreover when the observational noise is rather large the suggested procedure can give a more accurate estimate than lasso we illustrate our method on some real problems including survival analysis with microarray data

immersed interface methods for neumann and related problems in two and three dimensions we develop and apply a finite difference method to discretize the laplacian operator with neumann boundary conditions on an irregular domain using a regular cartesian grid the method is an extension of the immersed interface method developed by leveque and li siam j numer anal 31 1994 1019 1044 with careful selection of stencils the method is second order accurate and produces a matrix that is stable diagonally semidominant the method is illustrated on several two dimensional problems and one three dimensional problem

fourth order accurate evaluation of integrals in potential theory on exterior 3d regions we present two methods for the rapid high order accurate evaluation of integrals in potential theory on general unbounded 3d regions our methods allow for direct calculation of derivatives of the integrals as well one of the methods uses a fourth order compact stencil and the other uses a nonstandard variant of richardson extrapolation both methods involve calculation of discontinuities in high order derivatives of the integrals across the boundary of the integration region the extrapolation method in addition involves correction for the discontinuities in truncation error the number of operations required for the methods is essentially equal to twice the number of operations needed to solve poisson s equation on a regular grid both methods avoid problems associated with using quadrature methods to evaluate integrals with singular kernels numerical results are presented for experiments on a variety of geometries in free space

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

a recovery based a posteriori error estimator for h curl interface problems abstract this paper introduces a new recovery based a posteriori error estimator for the lowest order nedelec finite element approximation to the h curl interface problem the error estimator is analyzed by establishing both the reliability and the efficiency bounds and is supported by numerical results under certain assumptions it is proved that the reliability and efficiency constants are independent of the jumps of the coefficients

a boundary condition capturing method for poisson s equation on irregular domains interfaces have a variety of boundary conditions or jump conditions that need to be enforced the ghost fluid method gfm was developed to capture the boundary conditions at a contact discontinuity in the inviscid euler equations and has been extended to treat more general discontinuities such as shocks detonations and deflagrations and compressible viscous flows in this paper a similar boundary condition capturing approach is used to develop a new numerical method for the variable coefficient poisson equation in the presence of interfaces where both the variable coefficients and the solution itself may be discontinuous this new method is robust and easy to implement even in three spatial dimensions furthermore the coefficient matrix of the associated linear system is the standard symmetric matrix for the variable coefficient poisson equation in the absence of interfaces allowing for straightforward application of standard black box solvers

on the design of discontinuous galerkin methods for elliptic problems based on hybrid formulations abstract the objective of this paper is to present a new framework for the design of discontinuous galerkin dg methods for elliptic problems the idea is to start from a hybrid formulation of the problem involving as unknowns the main field in the interior of the element domains and its fluxes and traces on the element boundaries rather than working with this three field formulation fluxes are modeled using finite difference expressions and then the traces are determined by imposing continuity of fluxes although other strategies could be devised this procedure is applied to four elliptic problems namely the convection diffusion equation in the diffusion dominated regime the stokes problem the darcy problem and the maxwell problem we justify some well known dg methods with some modifications that in fact allow to improve the performance of the original methods particularly when the physical properties are discontinuous

a fictitious domain approach with spread interface for elliptic problems with general boundary conditions the aim of this article is to solve second order elliptic problems in an original physical domain using a fictitious domain method with a spread interface approach the main idea of the fictitious domain approach consists in immersing the original domain of study into a geometrically bigger and simpler one called fictitious domain as the spatial discretization is being performed in the fictitious domain this method allows the use of structured meshes the discretization is not boundary fitted to the original physical domain this paper describes several ways to impose dirichlet robin or neumann boundary conditions on a spread immersed interface without locally modifying the numerical scheme and without using lagrange multipliers the numerical applications focus on difiusion and convection problems in the unit disk with dirichlet or robin boundary conditions for such problems analytical solutions can be determined for a correctly chosen source term the numerical resolution is performed using a q1 finite element scheme the spread interface approach is then combined with a local adaptive mesh refinement algorithm in order to increase the precision in the vinicity of the immersed boundary the l2 norm of the errors is computed in order to evaluate the capability of the method immersed boundaries are found in many industrial applications like two phase flow simulations fluid structure interaction etc this article represents a first step towards the simulation of these kinds of applications

preconditioners for spectral element methods for elliptic and parabolic problems in this paper we propose preconditioners for spectral element methods for elliptic and parabolic problems these preconditioners are constructed using separation of variables and are easy to invert moreover they are spectrally equivalent to the quadratic forms which they are used to approximate