localisation in a cosserat continuum under static and dynamic loading conditions abstract localisation studies have been carried out for an unconstrained elasto plastic strain softening cosserat continuum because of the presence of an internal length scale in this continuum model a perfect convergence is found upon mesh refinement a finite constant width of the localisation zone and a finite energy dissipation are computed under static as well as under transient loading conditions because of the existence of rotational degrees of freedom in a cosserat continuum additional wave types arise and wave propagation becomes dispersive this has been investigated analytically and numerically for an elastic cosserat continuum and an excellent agreement has been found between both solutions

complete analysis of phase transitions and ensemble equivalence for the curie weiss potts model using the theory of large deviations we analyze the phase transition structure of the curie weiss potts spin model which is a mean field approximation to the nearest neighbor potts model it is equivalent to the potts model on the complete graph on n vertices the analysis is carried out both for the canonical ensemble and the microcanonical ensemble besides giving explicit formulas for the microcanonical entropy and for the equilibrium macrostates with respect to the two ensembles we analyze ensemble equivalence and nonequivalence at the level of equilibrium macrostates relating these to concavity and support properties of the microcanonical entropy the curie weiss potts model is the first statistical mechanical model for which such a detailed and rigorous analysis has been carried out

a discontinuous galerkin formulation of a model of gradient plasticity at finite strains abstract the work presented here constitutes an extension to the finite strain regime of a discontinuous galerkin based strain gradient plasticity formulation presented in djoko et al j k djoko f ebobisse a t mcbride b d reddy a discontinuous galerkin formulation for classical and gradient plasticity part 1 formulation and analysis comput methods appl mech engrg 196 2007 3881 3897 j k djoko f ebobisse a t mcbride b d reddy a discontinuous galerkin formulation for classical and gradient plasticity part 2 algorithms and numerical analysis comput methods appl mech engrg 197 2007 1 21 the focus here is on algorithmic and computational aspects of the formulation at finite strains the adoption of a logarithmic hyperelastic plastic formulation preserves the essential features of the infinitesimal formulation this key ingredient allows the predictor corrector solution algorithms developed for the infinitesimal gradient formulation to be extended readily to the finite strain regime the use of low order elements is essential to contain the computational expense of the formulation but these elements are prone to locking the method of enhanced assumed strains for geometrically nonlinear problems is utilised to circumvent this limitation the form of the consistent tangent modulus is derived for the case of gradient plasticity two numerical examples are presented to illustrate aspects of the approximation scheme and the algorithm as well as features of the model of gradient plasticity

an adaptive finite element algorithm for gradient theory of plasticity with coupling to damage in this paper an adaptive finite element algorithm is proposed for resolving localization zones in plastic solids with coupling to damage whereby the solution is regularized via the introduction of gradient theory it is possible to avoid the mesh dependence that otherwise would become inevitable when the localization narrows progressively due to the development of damage the adaptive method is based on an a posteriori error computation in energy norm for the plastic multiplier based on the constitutive subproblem which is solved via a mixed finite element formulation the performance of the proposed method is investigated for a few simple numerical examples in 1d and 2d

a high order accurate parallel solver for maxwell s equations on overlapping grids a scheme for the solution of the time dependent maxwell s equations on composite overlapping grids is described the method uses high order accurate approximations in space and time for maxwell s equations written as a second order vector wave equation high order accurate symmetric difference approximations to the generalized laplace operator are constructed for curvilinear component grids the modified equation approach is used to develop high order accurate approximations that use only three time levels and have the same time stepping restriction as the second order scheme discrete boundary conditions for perfect electrical conductors and for material interfaces are developed and analyzed the implementation is optimized for component grids that are cartesian resulting in a fast and efficient method the solver runs on parallel machines with each component grid distributed across one or more processors numerical results in two and three dimensions are presented for the fourth order accurate version of the method these results demonstrate the accuracy and efficiency of the approach

a conservative treatment of zonal boundaries for euler equation calculations finite difference calculations require the generation of a grid for the region of interest a zonal approach wherein the given region is subdivided into zones and the grid for each zone is generated independently makes the grid generation process for complicated topologies and for regions requiring selective grid refinement a fairly simple task this approach results in new boundaries within the given region that is zonal boundaries at the interfaces of the various zones the zonal boundary scheme the integration scheme used to update the points on the zonal boundary for the euler equations must be conservative accurate stable and applicable to general curvilinear coordinate systems a zonal boundary scheme with these desirable properties is developed in this study the scheme is designed for explicit first order accurate integration schemes but can be modified to accommodate second order accurate explicit and implicit integration schemes results for inviscid flow including supersonic flow over a cylinder blast wave diffraction by a ramp and one dimensional shock tube flow are obtained on zonal grids the conservative nature of the zonal boundary scheme permits the smooth transition of the discontinuities associated with these flows from one zone to another the calculations also demonstrate the continuity of contour lines across zonal boundaries that can be achieved with the present zonal scheme

the patch test for mixed formulations this paper presents a simple extension of the patch test to mixed formulations to provide the necessary and sufficient conditions for convergence the general algebraic conditions of babuska and brezzi are given a simple form and a conceptual application of the patch tests serves to point out the instability of several well known formulations for incompressible problems the test is also applied to mixed displacement strain formulations of elasticity and limiting formulations are noted

large deviations for the empirical field of a gibbs measure on applique le principe de contraction pour obtenir un principe de grande deviation pour la distribution des distributions empiriques on donne une description detaillee de la fonction de taux resultante

computational framework for coupling compressible and low mach number codes a hybrid multicode computational method is developed that allows combining compressible and low mach number codes into a single computational solver the proposed methodology can be used for integrated simulations of multicomponent engineering problems the unsteady coupling between the two codes is performed via exchanging time dependent state information through the interfaces using overlapping meshes the proposed coupling procedure is validated on laminar steady and unsteady test problems both constant density and variable density regimes of the low mach number code are investigated results of the test cases performed in frameworks of one way coupling and two way coupling are documented several methods of implementation of interface conditions for the compressible code are compared by looking at the levels of dilatation as well as numerical errors the injection method when all the variables are simply transformed from one code to another through an interpolation procedure is chosen for its superior stability accuracy and simplicity of implementation numerical results for the coupled calculations obtained with the injection method show good agreement with both standalone calculations and analytical solutions

an overlapping schwarz method for spectral element solution of the incompressible navier stokes equations efficient solution of the navier stokes equations in complex domains is dependent upon the availability of fast solvers for sparse linear systems for unsteady incompressible flows the pressure operator is the leading contributor to stiffness as the characteristic propagation speed is infinite in the context of operator splitting formulations it is the pressure solve which is the most computationally challenging despite its elliptic origins we examine several preconditioners for the consistent l2poisson operator arising in the pn pn 2spectral element formulation of the incompressible navier stokes equations we develop a finite element based additive schwarz preconditioner using overlapping subdomains plus a coarse grid projection operator which is applied directly to the pressure on the interior gauss points for large two dimensional problems this approach can yield as much as a fivefold reduction in simulation time over previously employed methods based upon deflation

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

a fourth order accurate method for the incompressible navier stokes equations on overlapping grids a method is described to solve the time dependent incompressible navier stokes equations with finite differences on curvilinear overlapping grids in two or three space dimensions the scheme is fourth order accurate in space and uses the momentum equations for the velocity coupled to a poisson equation for the pressure the boundary condition for the pressure is taken as u 0 extra numerical boundary conditions are chosen to make the scheme accurate and stable the velocity is advanced explicitly in time any standard time stepping scheme such as runge kutta can be used the poisson equation is solved using direct or iterative sparse matrix solvers or by the multigrid algorithm computational results in two and three space dimensions are given

a staggered grid multidomain spectral method for the compressible navier stokes equations we present a flexible non conforming staggered grid chebyshev spectral multidomain method for the solution of the compressible navier stokes equations in this method subdomain corners are not included in the approximation thereby simplifying the subdomain connectivity to allow for local refinement of the polynomial order or subdomain subdivision non conforming subdomains are treated by a mortar method spectral accuracy is shown in one and two dimensional linear and non linear problems application is made to four compressible flow problems the couette flow a steady boundary layer over a flat plate steady transonic flow in a nozzle and unsteady flow over a cylinder at a reynolds number of 75

large deviations for gibbs random fields a large deviation principle for gibbs random fields on zd is proven and a corresponding large deviations proof of the gibbs variational formula is given a generalization of the lanford theory of large deviations is also obtained

moving overlapping grids with adaptive mesh refinement for high speed reactive and non reactive flow we consider the solution of the reactive and non reactive euler equations on two dimensional domains that evolve in time the domains are discretized using moving overlapping grids in a typical grid construction boundary fitted grids are used to represent moving boundaries and these grids overlap with stationary background cartesian grids block structured adaptive mesh refinement amr is used to resolve fine scale features in the flow such as shocks and detonations refinement grids are added to base level grids according to an estimate of the error and these refinement grids move with their corresponding base level grids the numerical approximation of the governing equations takes place in the parameter space of each component grid which is defined by a mapping from fixed parameter space to moving physical space the mapped equations are solved numerically using a second order extension of godunov s method the stiff source term in the reactive case is handled using a runge kutta error control scheme we consider cases when the boundaries move according to a prescribed function of time and when the boundaries of embedded bodies move according to the surface stress exerted by the fluid in the latter case the newton euler equations describe the motion of the center of mass of the each body and the rotation about it and these equations are integrated numerically using a second order predictor corrector scheme numerical boundary conditions at slip walls are described and numerical results are presented for both reactive and non reactive flows that demonstrate the use and accuracy of the numerical approach

multiplicative schwarz methods for parabolic problems based on domain decomposition we give a multiplicative schwarz domain decomposition method for semi linear parabolic problems we consider the relation between the convergence rate and discretization parameters including the diameter of the sub domain we give the error estimate which tells us that the convergence of the approximate solution is independent of the iteration number at each time level finally we give a numerical example

effects of softening in elastic plastic structural dynamics softening understood as unstable behaviour of material or structural components is considered herein for its possible consequences on the overall behaviour of discrete dynamic models of elastic plastic beam structures in the absence of geometric effects it is shown that multiplicity of incremental solutions response bifurcations and manifestations of overall instability may occur the bifurcated responses may exhibit different scenarios for the same excitation e g shakedown or local damage up to failure an insight into and criteria for such occurrences are achieved by formulating the rate and the finite increment problem as linear complementarity problems period doubling and deterministic chaos under harmonic excitation are observed when reversible holonomic softening behaviour is assumed the strong sensitivity of the phenomena investigated with respect to the choice of the structural model is pointed out the findings are elucidated by illustrative examples

preconditioned spectral multi domain discretization of the incompressible navier stokes equations in this paper we present a preconditioned multi domain algorithm applied to the elliptic kernels arising from the spectral collocation of the incompressible navier stokes equations in three space dimensions with one homogeneous direction the technique based on the iterative solution of the schur complement matrix allows for efficient numerical solution of the operators in complex geometries consisting of a collection of non overlapping rectangular subdomains the method is shown to be nearly optimal in terms of condition number behavior in a double path of refinement strategy i e whenever both the number of legendre modes and the number of subdomains are significantly increased it is thus well suited for engineering applications in the fields of direct numerical simulation and large eddy simulation of turbulence

large deviations for markov processes with resetting markov processes restarted or reset at random times to a fixed state or region in space have been actively studied recently in connection with random searches foraging and population dynamics here we study the large deviations of time additive functions or observables of markov processes with resetting by deriving a renewal formula linking generating functions with and without resetting we are able to obtain the rate function of such observables characterizing the likelihood of their fluctuations in the long time limit we consider as an illustration the large deviations of the area of the ornstein uhlenbeck process with resetting other applications involving diffusions random walks and jump processes with resetting or catastrophes are discussed

a spectrally accurate method for overlapping grid solution of incompressible navier stokes equations an overlapping mesh methodology that is spectrally accurate in space and up to third order accurate in time is developed for solution of unsteady incompressible flow equations in three dimensional domains the ability to decompose a global domain into separate but overlapping subdomains eases mesh generation procedures and increases flexibility of modeling flows with complex geometries the methodology employs implicit spectral element discretization of equations in each subdomain and explicit treatment of subdomain interfaces with spectrally accurate spatial interpolation and high order accurate temporal extrapolation and requires few if any iterations yet maintains the global accuracy and stability of the underlying flow solver the overlapping mesh methodology is thoroughly validated using two dimensional and three dimensional benchmark problems in laminar and turbulent flows the spatial and temporal convergence is documented and is in agreement with the nominal order of accuracy of the solver the influence of long integration times as well as inflow outflow global boundary conditions on the performance of the overlapping grid solver is assessed in a turbulent benchmark of fully developed turbulent pipe flow the turbulent statistics with the overlapping grids is validated against published available experimental and other computation data scaling tests are presented that show near linear strong scaling even for moderately large processor counts a novel overlapping grid method for incompressible flow equations is developed the overlapping grid method is based on spectral element discretization the method maintains global spectral accuracy with polynomial refinement an explicit interface extrapolation scheme for temporal coupling is proposed the scheme allows for high order temporal accuracy with only a few iterations

adaptive local overlapping grid methods for parabolic systems in two space dimensions abstract adaptive mesh refinement techniques are described for two dimensional systems of parabolic partial differential equations solutions are calculated using galerkin s method with a piecewise bilinear basis in space and backward euler integration in time a posteriori estimates of the local discretization error of piecewise bilinear finite element solutions are obtained by a p refinement technique these error estimates are used to control a local h refinement strategy where finer grids are recursively introduced in regions where a prescribed tolerance is exceeded fine grids at a given level of refinement may overlap each other and independent solutions are generated on each of them a version of the schwarz alternating principle is used to coordinate solutions between overlapping fine grids computational results demonstrating the performance of the adaptive procedure on linear and nonlinear problems and apparent convergence of the error estimate for linear heat conduction problem and uniform global refinement is presented

theory and numerics of a thermodynamically consistent framework for geometrically linear gradient plasticity the paper presents the theory and the numerics of a thermodynamically consistent formulation of gradient plasticity at small strains starting from the classical local continuum formulation which fails to produce physically meaningful and numerically converging results within localization computations a thermodynamically motivated gradient plasticity formulation is envisioned the model is based on an assumption for the helmholtz free energy incorporating the gradient of the internal history variable a yield condition and the postulate of maximum dissipation resulting in an associated structure as a result the driving force conjugated to the hardening evolution is identified as the quasi non local drag stress which incorporates besides the strictly local drag stress essentially the divergence of a vectorial hardening flux at the numerical side besides the balance of linear momentum the algorithmic consistency condition has to be solved in weak form thereby the crucial issue is the determination of the active constraints exhibiting plastic loading which is solved by an active set search algorithm borrowed from convex non linear programming moreover different discretization techniques are proposed in order to compare the fe performance in local plasticity with the advocated gradient formulation both for hardening and softening copyright 2001 john wiley sons ltd

finite elements with embedded displacement discontinuity a generalized variable formulation this paper is intended to bring a contribution towards a satisfactory simulation of those fracture phenomena which result in the appearance and development of discrete cracks to this purpose a general mixed finite element formulation is proposed based on the concept of generalized variables in prager s sense the displacement field inside an element is modelled by the sum of two contributions a regular continuous part which is governed by standard shape functions and a possibly discontinuous one which is introduced soon after a suitable criterion is satisfied the formulation is first specialized to a one dimensional case then a triangular element for two dimensional problems is described in detail analytical and numerical examples are presented in order to clarify the formulation and to point out the essential role of inter element conformity copyright 2000 john wiley sons ltd

exact partition function zeros of the wako saito munoz eaton protein model i compute exact partition function zeros of the wako saito munoz eaton model for various secondary structural elements and for two proteins 1bbl and 1i6c by using both analytic and numerical methods two state and barrierless downhill folding transitions can be distinguished by a gap in the distribution of zeros at the positive real axis

a finite difference domain decomposition algorithm for numerical solution of the heat equation a domain decomposition algorithm for numerically solving the heat equation in one and two space dimensions is presented in this procedure interface values between subdomains are found by an explicit finite difference formula once these values are calculated interior values are determined by backward differencing in time a natural extension of this method allows for the use of different time steps in different subdomains maximum norm error estimates for these procedures are derived which demonstrate that the error incurred at the interfaces is higher order in the discretization parameters

large deviations and maximum entropy principle for interacting random fields on mathbb z d we present a new approach to the principle of large deviations for the empirical field of a gibbsian random field on the integer lattice zd this approach has two main features first we can replace the traditional weak topology by the finer topology of convergence of cylinder probabilities and thus obtain estimates which are finer and more widely applicable second we obtain as an immediate consequence a limit theorem for conditional distributions under conditions on the empirical field the limits being those predicted by the maximum entropy principle this result implies a general version of the equivalence of gibbs ensembles stating that every microcanonical limiting state is a grand canonical equilibrium state we also prove a converse to the last statement and discuss some applications 0 introduction as is well known the study of the asymptotic probabilities of large fluctuations of time averages or space averages away from the mean is based on two fundamental principles the principle of large deviations and the maximum entropy principle the former provides the exact rate of exponential decay of the fluctuation probabilities whereas the latter predicts the limiting conditional distribution under the condition that the fluctuations are large it is obvious that these principles are intimately related in this paper we investigate these principles in the case of interacting random fields on the integer lattice 7zd the setup is the following first we let e 4 be any measurable space we shall assume throughout that e a is standard borel but we shall avoid making any explicit topological assumptions on e so we do not assume that e is polish next we let s zd be the d dimensional integer lattice and q f e gas the associated product space in lattice models of statistical mechanics e is called the state space or single spin space and q f the configuration space we let j s f denote the set of all probability

a discrete formulation for elastic solids with damaging interfaces abstract an elastic solid with embedded interfaces loci of possible displacement discontinuities is considered here decohesion and quasi brittle fracture processes are simulated by making use of softening interface laws which relate tractions to displacements jumps a discrete formulation in terms of interface variables only is obtained the space discretization is carried out by means of a mixed finite element approach in which all interface variables are modelled study of uniqueness of the rate problem formulated in terms of interface variables and of stability of the equilibrium states is presented some examples are shown in order to clarify the formulation

a uniform strain hexahedron and quadrilateral with orthogonal hourglass control the treatment of zero energy modes which arise due to one point integration of first order isoparametric finite elements is addressed a method for precisely isolating these modes for arbitrary geometry is developed two hourglass control schemes viscous and elastic are presented and examined in addition a convenient one point integration scheme which analytically integrates the element volume and uniform strain modes is presented

on the use of a high order overlapping grid method for coupling in cfd caa this paper presents a theoretical analysis and two applications of a high order overlapping grid method for coupling cartesian and curvilinear grids developed in order to simulate aerodynamic noise first the overlapping grid method based on lagrange interpolating polynomials is described and a theoretical analysis of the interpolation operator is then carried out it shows that the interpolation generates spurious modes that depend on the wavenumbers of the signal besides it also gives the optimal conditions in which interpolation can be applied then an application to the simulation of the aeroacoustic noise generated by the vortex shedding behind a cylinder is presented during this simulation it appears that interpolation can create some spurious acoustic modes in regions where hydrodynamic fluctuations are significant as predicted by the theoretical analysis it is shown that these spurious modes disappear when a refined cartesian grid is used 26 points per wavelength of the vortex shedding were found to be adequate in this study at last the simulation of the aerodynamic noise of a three element high lift wing profile has then been carried out for this application the main acoustic source at the slat trailing edge is represented analytically the propagation of the generated acoustic wave is simulated with a mean flow at rest and with a steady turbulent mean flow computed by rans the first application allows us to assess the method by comparing the results to a reference solution the second one shows that the influence of a non uniform mean flow on the directivity of an acoustic source can be observed in complex geometries this application therefore shows that the proposed coupling method is well adapted to complex geometries that are usually met in industrial applications

high order compact finite difference methods on general overset grids this work investigates the coupling of a very high order finite difference algorithm for the solution of conservation laws on general curvilinear meshes with overset grid techniques originally developed to address complex geometric configurations the solver portion of the algorithm is based on pade type compact finite differences of up to sixth order with up to 10th order filters employed to remove spurious waves generated by grid non uniformities boundary conditions and flow non linearities the overset grid approach is utilized as both a domain decomposition paradigm for implementation of the algorithm on massively parallel machines and as a means for handling geometric complexity in the computational domain two key features have been implemented in the current work the ability of the high order algorithm to accommodate holes cut in grids by the overset grid approach and the use of high order interpolation at non coincident grid overlaps several high order high accuracy interpolation methods were considered and a high order explicit non optimized lagrangian method was found to be the most accurate and robust for this application several two dimensional benchmark problems were examined to validate the interpolation methods and the overall algorithm these included grid to grid interpolation of analytic test functions the inviscid convection of a vortex laminar flow over single and double cylinder configurations and the scattering of acoustic waves from one and three cylinder configurations the employment of the overset grid techniques coupled with high order interpolation at overset boundaries was found to be an effective way of employing the high order algorithm for more complex geometries than was previously possible

a chimera method based on a dirichlet neumann robin coupling for the navier stokes equations we present a chimera method for the numerical solution of incompressible flows past objects in relative motion the chimera method is implemented as an iteration by subdomain method based on dirichlet neumann robin coupling the dd method we propose is not only geometric but also algorithmic for the solution on each subdomain is obtained on separate processes and the exchange of information between the subdomains is carried out by a master code this strategy is very flexible as it requires almost no modification to the original numerical code therefore only the master code has to be adapted to the numerical codes and the strategies used on each subdomain as a basic flow solver we a use stabilized finite element method

a large deviation perspective on ratio observables in reset processes robustness of rate functions we study large deviations of a ratio observable in discrete time reset processes the ratio takes the form of a current divided by the number of reset steps and as such it is not extensive in time a large deviation rate function can be derived for this observable via contraction from the joint probability density function of current and number of reset steps the ratio rate function is differentiable and we argue that its qualitative shape is robust i e it is generic for reset processes regardless of whether they have short or long range correlations we discuss similarities and differences with the rate function of the efficiency in stochastic thermodynamics

a domain decomposition method for incompressible viscous flow a method for using domain decomposition to solve the equations of incompressible viscous flow is presented the method is described in detail and test results are given for two test problems a notable feature of the method is that the incompressibility constraint is never explicitly imposed the domain decomposition uses finite difference and spectral methods on overlapping domains with second order accurate interpolation of the velocity relating the solutions on the different domains the method is shown to be globally second order accurate by the test results

on nonconservative algorithms for grid interfaces in computations of fluid flows by domain decomposition methods the necessity of conservation at grid interfaces is now widely claimed in this paper we consider nonconservative algorithms for the interfaces our investigation begins with discussions about typical nonconservative interface treatments for one dimensional calculations the analysis shows that the conservation error of a numerical solution caused by a nonconservative interface matching has an upper bound when the solution itself is bounded furthermore if the numerical solution converges as the mesh size goes to zero it converges to a weak solution of the problem under certain conditions that may be detected numerically also we acquire similar results for two dimensional calculations on grids intersecting with each other in an arbitrary way in order to illustrate the theoretical results we present numerical examples and demonstrate that under those conditions conservation error reduces and accuracy for jumps as well as locations of discontinuities improves as the mesh size decreases

on uniqueness of the dynamic finite step problem in gradient dependent softening plasticity the dynamic evolution of an elastoplastic softening solid is considered a material model including in the yield function the laplacian of the plastic multiplier is used to regularize the problem the dynamic finite step problem is formulated according to a generalized mid point integration scheme space discretization is carried out by a mixed finite element technique based on generalized variables a sufficient uniqueness condition of the finite step solution is proved for a one dimensional problem also a necessary and sufficient condition is presented a simple numerical test shows the regularizing properties mesh independence of the proposed model and the positive influence of the gradient term also on the time step amplitude ensuring uniqueness of solution

generalized midpoint finite element dynamic analysis of elastoplastic systems the application of a four parameter generalized midpoint integration method to small displacement elastoplastic dynamic problems is considered the integration of the governing relations over a time step is formulated as a non linear algebraic problem finite step problem by making use of generalized variables finite element modelling uniqueness and an extremum characterization of the solution of the finite step problem are discussed sufficient conditions for convergence of various iterative algorithms including newton raphson are given the b stability of the proposed procedure is assessed

a wavenumber based extrapolation and interpolation method for use in conjunction with high order finite difference schemes the errors incurred in using extrapolation and interpolation in large scale computations are analyzed and quantified in the wavenumber space if a large extrapolation stencil is used the errors in the low wavenumbers can be significantly reduced however the errors in the high wavenumbers are at the same time greatly increased the opposite is true if the stencil size is reduced based on the wavenumber analysis an optimized extrapolation and interpolation method is proposed the optimization is carried out over a selected band of wavenumbers it is known that extrapolation often leads to numerical instability the instability is the result of large error amplification in the high wavenumber range to reduce the tendency to trigger numerical instability it is proposed that an extra constraint be imposed on the optimized extrapolation method the added constraint aims to reduce error amplification over the high wavenumbers numerical examples are provided to illustrate that accurate and stable numerical results can be obtained in large scale simulation using a high order finite difference scheme and the proposed optimized extrapolation method when the same problems are recomputed using the familiar high order polynomials extrapolation method in the lagrange form in one case the numerical results are plagued by large errors and ultimately instability in another problem it is found that the use of the lagrange polynomials extrapolation method would lead immediately to numerical instability

composite overlapping meshes for the solution of partial differential equations abstract we discuss the generation of curvilinear composite overlapping grids and the numerical solution of partial differential equations on them a composite overlapping grid consists of a set of curvilinear component grids that cover a region and overlap where they meet continuity conditions interpolation are imposed at the overlapping boundaries the principal advantage of composite grids is in the generation of grids for regions with complicated geometries the grid construction program cmpgrd is used to create composite grids with any number of component grids we describe some techniques for the solution of elliptic and time dependent pdes on composite meshes applications to the solution of the compressible navier stokes equations are presented

theoretical and computational aspects of a thermodynamically consistent framework for geometrically linear gradient damage this paper presents the theory and the numerics of an isotropic gradient damage formulation within a thermodynamical background the main motivation is provided by localization computations whereby classical local continuum formulations fail to produce physically meaningful and numerically converging results we propose a formulation in terms of the helmholtz free energy incorporating the gradient of the damage field a dissipation potential and the postulate of maximum dissipation as a result the driving force conjugated to damage evolution is identified as the quasi nonlocal energy release rate which essentially incorporates the divergence of a vectorial damage flux besides the strictly local energy release rate on the numerical side besides balance of linear momentum the algorithmic consistency condition must be solved in weak form thereby the crucial issue is the selection of active constraints which is solved by an active set search algorithm borrowed from convex nonlinear programming in the examples we compare the behavior in local damage with the performance of the gradient formulation

solvable models of classical lattice gases we study classical lattice gases at fixed temperature but variable fugacity it is shown how the thermodynamic functions may be calculated exactly provided the boltzmann weights are representable as principal minors of a convolution operator we explicitly construct this operator for the cluster models of fisher and felderhof

chebyshev collocation method and multi domain decomposition for navier stokes equations in complex curved geometries a general multidomain decomposition is proposed for the numerical solution of the 2d incompressible stationary navier stokes equations the solution technique consists in a chebyshev orthogonal collocation method preconditioned by a standard galerkin finite element technique the preconditioned system is then solved through a richardson procedure the domain of interest is decomposed into quadrilaterals curved when needed a gordon transfinite interpolation performs the curvilinear grid generation of the obtained simply connected planar subdomains the interface conditions naturally incorporated into the finite element approach relate neighbour subdomains through the normal jump of appropriate fluxes across internal boundaries where an integral form ofc 1 continuity is consequently achieved at convergence of the iterative processes the study of model stokes problems demonstrates that the current method still behaves spectrally in distorted geometries for curvilinear distortion a loss of several orders of magnitude is observed in the solution accuracy even when the distortion is very limited finally some results of flow simulation in a constricted channel are proposed to illustrate the abilities of the present method to treat navier stokes problems

gradient dependent plasticity formulation and algorithmic aspects a plasticity theory is proposed in which the yield strength not only depends on an equivalent plastic strain measure hardening parameter but also on the laplacian thereof the consistency condition now results in a differential equation instead of an algebraic equation as in conventional plasticity to properly solve the set of non linear differential equations the plastic multiplier is discretized in addition to the usual discretization of the displacements for appropriate boundary conditions this formulation can also be derived from a variational principle accordingly the theory is complete

parameter identification in explicit structural dynamics performance of the extended kalman filter abstract an extended kalman filter ekf approach is adopted in this paper for structural systems subject to dynamic loadings to simultaneously estimate the state and calibrate constitutive parameters to pursue this aim after space and time discretizations a joint system state vector is introduced gathering nodal displacements and constitutive parameters to be identified because of the linearization of the discretized equations governing filter updates the ekf can lead to inaccurate model calibration it is shown that even though the state of the system is always followed with a high level of accuracy unsatisfactory parameter estimations can be obtained in the case of degrading strength of the structure that is in the case of softening both single degree of freedom dof and multi dof dynamic systems are analyzed to detect the possible sources of this inaccuracy

a spectral element method for fluid dynamics laminar flow in a channel expansion a spectral element method that combines the generality of the finite element method with the accuracy of spectral techniques is proposed for the numerical solution of the incompressible navier stokes equations in the spectral element discretization the computational domain is broken into a series of elements and the velocity in each element is represented as a high order lagrangian interpolant through chebyshev collocation points the hyperbolic piece of the governing equations is then treated with an explicit collocation scheme while the pressure and viscous contributions are treated implicitly with a projection operator derived from a variational principle the implementation of the technique is demonstrated on a one dimensional inflow outflow advection diffusion equation and the method is then applied to laminar two dimensional separated flow in a channel expansion comparisons are made with experiment and previous numerical work

a unified approach for variationally consistent finite elements in elastoplasticity abstract a general variationally consistent approach for elastoplastic finite element analysis is presented a weak formulation of the finite step i e in terms of finite increments elastoplastic boundary value problem is obtained by enforcing the stationarity of a mixed functional of hu washizu type a peculiarity of the present formulation is that the constitutive equations are expressed in weak form also all fields are then independently discretized in terms of so called generalized variables the result is a set of non linear algebraic equations a subset of which can be interpreted as the backward difference time integrated constitutive law of the finite element an iterative procedure where all calculations are at the element and not at the gauss point level resting on newton raphson scheme is proposed for the solution of the governing equations it is shown how the customary displacement approach and also some other mixed formulations in elastoplasticity can be regarded as special cases of the present formulation

a variational principle for gradient plasticity abstract we elaborate on a generalized plasticity model which belongs to the class of gradient models suggested earlier by aifantis and co workers the generalization of the conventional theory of plasticity has been accomplished by the inclusion of higher order spatial gradients of the equivalent plastic strain in the yield condition first it is shown how these gradients affect the critical condition for the onset of localization and allow for a wavelength selection analysis leading to estimates for the width and or spacing of shear bands due to the presence of higher order gradients additional boundary conditions for the equivalent plastic strain are required this question and also the associated problem of the formulation and solution of general boundary value problems were left open in the previous work we demonstrate here that upon assuming a certain type of additional boundary conditions the structural symmetries of the gradient dependent constitutive model are such that there exists a variational principle for the displacement rates and the rate of the equivalent plastic strain the variational principle can serve as a basis for the numerical solution of boundary value problems in the sense of the finite element method explicit expressions for the tangent stillness matrix and the generalized nodal point forces are given

stability analysis of interface temporal discretization in grid overlapping methods we investigate the stability of a temporal discretization of interface terms in grid overlapping methods a matrix stability analysis is performed on a model problem of the one dimensional diffusion equation on overlapping grids the scheme stability is first analyzed theoretically and a proof of the unconditional stability of the first order interface extrapolation scheme with the first and second order time integration for any overlap size is presented for the higher order schemes we obtain explicit estimates of the spectral radius of the corresponding discrete matrix operator and document the values of the stability threshold depending on the number of grid points and the size of overlap the influence of iterations on stability properties is also investigated numerical experiments are then presented relating the obtained stability bounds to the observed numerical values semidiscrete analysis confirms the derived scaling for the stability bounds

analysis of a fractional step method on overset grids a fractional step method for solving the incompressible navier stokes equations on overset grids is derived as a matrix factorization of the spatially and temporally discretized system of equations the algorithm is applied to several test problems using second order accurate finite volume flux differencing on staggered grid systems and a hybrid implicit explicit time advancement scheme spatial order of accuracy is shown to depend on the behavior of the overset grid overlap during grid refinement the temporal order of accuracy of the time advancement algorithm on a single grid is maintained on the overset grid

a discontinuous galerkin formulation for classical and gradient plasticity part 2 algorithms and numerical analysis this work is the second of a two part investigation into the use of discontinuous galerkin methods for obtaining approximate solutions to problems of classical and gradient plasticity part i j k djoko f ebobisse a t mcbride b d reddy a discontinuous galerkin formulation for classical and gradient plasticity part 1 formulation and analysis comput methods appl mech engrg 196 2007 3881 3897 presented the formulation and analysis of such problems this part focusses on algorithmic and computational aspects of the problem in particular it is shown that the predictor corrector algorithms of classical plasticity are readily extended to the case of gradient plasticity and to discontinuous galerkin formulations conditions for convergence of the algorithms are presented for the elastic secant and consistent tangent predictors the form of the consistent tangent modulus is established for the case of gradient plasticity a selection of numerical examples is presented and discussed with a view to illustrating aspects of the approximation scheme and algorithms as well as features of the model of gradient plasticity adopted here

on multigrid for overlapping grids the solution of elliptic partial differential equations on composite overlapping grids using multigrid is discussed an approach is described that provides a fast and memory efficient scheme for the solution of boundary value problems in complex geometries the key aspects of the new scheme are an automatic coarse grid generation algorithm an adaptive smoothing technique for adjusting residuals on different component grids and the use of local smoothing near interpolation boundaries other important features include optimizations for cartesian component grids the use of over relaxed red black smoothers and the generation of coarse grid operators through galerkin averaging numerical results in two and three dimensions show that very good multigrid convergence rates can be obtained for dirichlet neumann and mixed boundary conditions a comparison to krylov based solvers shows that the multigrid solver can be much faster and require significantly less memory