two level method based on finite element and crank nicolson extrapolation for the time dependent navier stokes equations a fully discrete two level finite element method the two level method is presented for solving the two dimensional time dependent navier stokes problem the method requires a crank nicolson extrapolation solution u h tau 0 p h tau 0 on a spatial time coarse grid j h tau 0 and a backward euler solution u h tau p h tau on a space time fine grid j h tau the error estimates of optimal order of the discrete solution for the two level method are derived compared with the standard crank nicolson extrapolation method the one level method based on a space time fine grid j h tau the two level method is of the error estimates of the same order as the one level method in the h1 norm for velocity and the l2 norm for pressure however the two level method involves much less work than the one level method

discontinuous galerkin approximations of the stokes and navier stokes equations numerical schemes to compute approximate solutions of the evolutionary stokes and navier stokes equations are studied the schemes are discontinuous in time and conforming in space and of arbitrarily high order fully discrete error estimates are derived and dependence of the viscosity constant is carefully tracked it is shown that the errors are bounded by projection errors of the exact solution which exhibit optimal rates when the solutions are smooth

convergence analysis and error estimates for a parallel algorithm for solving the navier stokes equations this paper is concerned with the analysis of the convergence and the derivation of error estimates for a parallel algorithm which is used to solve the incompressible navier stokes equations as usual the main idea is to split the main differential operator this allows to consider independently the two main difficulties namely nonlinearity and incompressibility the results justify the observed accuracy of related numerical results

maximum l 2 convergence rates of the crank nicolson scheme to the stokes initial value problem let a denote the stokes operator and d a alpha the domain of its fractional powers a alpha we consider the homogeneous stokes initial value problem with initial data u 0 u 0 in d a 1 varepsilon varepsilon in 0 1 for stokes like equations the range varepsilon in 0 frac 1 4 is of special interest as any solution derived from varepsilon ge frac 1 4 would necessarily have to satisfy an additional in practice unverifiable compatibility condition at time t 0 approximating any strong solution u in c 0 0 infty d a 1 varepsilon in time direction on a finite time interval 0 t with a crank nicolson scheme we show convergence of order o frac tau 2 t 1 varepsilon which is maximal for the assumed regularity of the data and reflects the loss of regularity as t rightarrow 0 the error estimates are derived by energy and semigroup methods combined with a parabolic duality argument

on error estimates of the projection methods for the navier stokes equations second order schemes we present in this paper a rigorous error analysis of several projection schemes for the approximation of the unsteady incompressible navier stokes equations the error analysis is accomplished by interpreting the respective projection schemes as second order time discretizations of a perturbed system which approximates the navier stokes equations numerical results in agreement with the error analysis are also presented

on the convergence in distribution of measurable multifunctions random sets normal integrands stochastic processes and stochastic infima the concept of the distribution function of a closed valued measurable multifunction is introduced and used to study the convergence in distribution of sequences of multifunctions and the epi convergence in distribution of normal integrands and stochastic processes in particular various compactness criteria are exhibited the connections with the classical convergence theory for stochastic processes are analyzed and for purposes of illustration we apply the theory to sketch out a modified derivation of donsker s theorem brownian motion as a limit of normalized random walks we also suggest the potential application of the theory to the study of the convergence of stochastic infima

a staggered space time discontinuous galerkin method for the three dimensional incompressible navier stokes equations on unstructured tetrahedral meshes in this paper we propose a novel arbitrary high order accurate semi implicit space time discontinuous galerkin method for the solution of the three dimensional incompressible navier stokes equations on staggered unstructured curved tetrahedral meshes as is typical for space time dg schemes the discrete solution is represented in terms of space time basis functions this allows to achieve very high order of accuracy also in time which is not easy to obtain for the incompressible navier stokes equations similarly to staggered finite difference schemes in our approach the discrete pressure is defined on the primary tetrahedral grid while the discrete velocity is defined on a face based staggered dual grid while staggered meshes are state of the art in classical finite difference schemes for the incompressible navier stokes equations their use in high order dg schemes is still quite rare a very simple and efficient picard iteration is used in order to derive a space time pressure correction algorithm that achieves also high order of accuracy in time and that avoids the direct solution of global nonlinear systems formal substitution of the discrete momentum equation on the dual grid into the discrete continuity equation on the primary grid yields a very sparse five point block system for the scalar pressure which is conveniently solved with a matrix free gmres algorithm from numerical experiments we find that the linear system seems to be reasonably well conditioned since all simulations shown in this paper could be run without the use of any preconditioner even up to very high polynomial degrees for a piecewise constant polynomial approximation in time and if pressure boundary conditions are specified at least in one point the resulting system is in addition symmetric and positive definite this allows us to use even faster iterative solvers like the conjugate gradient method the flexibility and accuracy of high order space time dg methods on curved unstructured meshes allows to discretize even complex physical domains with very coarse grids in both space and time the proposed method is verified for approximation polynomials of degree up to four in space and time by solving a series of typical 3d test problems and by comparing the obtained numerical results with available exact analytical solutions or with other numerical or experimental reference data to the knowledge of the authors this is the first time that a space time discontinuous galerkin finite element method is presented for the three dimensional incompressible navier stokes equations on staggered unstructured tetrahedral grids high order staggered space time dg scheme for the 3d incompressible navier stokes equations primal grid based on curved tetrahedra dual grid based on non standard five point hexahedra new discretization of the viscous terms on the main grid using lifting operators on the dual grid very simple space time pressure correction algorithm providing also high order in time all linear systems have maximum sparsity and use the smallest possible optimal five point stencil

long term dissipativity of time stepping algorithms for an abstract evolution equation with applications to the incompressible mhd and navier stokes equations abstract this paper examines the long term dissipativity and unconditional non linear stability of time integration algorithms for the incompressible mhd equations the results apply as a particular case to the uncoupled incompressible navier stokes equations the analysis is performed in the setting defined by an abstract evolution equation possessing the same long term dissipative structure as the above physical systems a general class of time stepping schemes is analyzed in this context showing the non linear unconditional stability and long term dissipativity of particular members of this family of algorithms in particular algorithms showing these strong properties which additionally involve the solution of a linear problem at each time step are identified a particular linear algorithm is shown to exhibit optimal long term dissipative properties and to possess an algorithmic global attractor these abstract results apply rigorously to galerkin mixed finite element discretizations of the problem to this end it is shown that such galerkin type projection inherits the dissipative properties of the continuum problem computational aspects involved in the numerical implementation of these methods are addressed and illustrated in representative numerical simulations

on the transport diffusion algorithm and its applications to the navier stokes equations this paper deals with an algorithm for the solution of diffusion and or convection equations where we mixed the method of characteristics and the finite element method globally it looks like one does one step of transport plus one step of diffusion or projection but the mathematics show that it is also an implicit time discretization of thepde in lagrangian form we give an error bound h t h h t in the interesting case that holds also for the navier stokes equations even when the reynolds number is infinite euler equation

generation of vorticity near the boundary in planar navier stokes flows we construct the solutions of the planar navier stokes flow for a viscous incompressible fluid in the half plane by means of a boundary layer equation describing the production of vorticity on the boundary regularity properties are also discussed

convergence of euler stokes splitting of the navier stokes equations we consider approximation by partial time steps of a smooth solution of the navier stokes equations in a smooth domain in two or three space dimensions with no slip boundary condition for small k 0 we alternate the solution for time k of the inviscid euler equations with tangential boundary condition and the solution of the linear stokes equations for time k with the no slip condition imposed we show that this approximation remains bounded in h2 p and is accurate to order k in lp for p the principal difficulty is that the initial state for each stokes step has tangential velocity at the boundary generated during the euler step and thus does not satisfy the boundary condition for the stokes step the validity of such a fractional step method or splitting is an underlying principle for some computational methods 1994 john wiley sons inc

asymptotic behavior and time discretization analysis for the non stationary navier stokes problem the asymptotic behavior and the euler time discretization analysis are presented for the two dimensional non stationary navier stokes problem if the data and f t satisfy a uniqueness condition corresponding to the stationary navier stokes problem we then obtain the convergence of the non stationary navier stokes problem to the stationary navier stokes problem and the uniform boundedness stability and error estimates of the euler time discretization for the non stationary navier stokes problem

vibration and stability analyses of cross ply laminated circular cylindrical shells abstract an improved analytical procedure is introduced in the context of free vibration and stability problems of cross ply laminated circular cylindrical shells the donnell shear deformation type theory and donnell s classical theory are used to illustrate the procedure numerical results are presented for shells with a variety of boundary conditions special attention is given to the axisymmetric problems and new results not yet found in the literature for stability and vibration are reported it is recommended that the new technique be used for generating the levy type solutions in vibration and stability problems of laminated shell panels and plates when the number and order of equations are higher than those of the classical theories

stability and convergence of semi implicit time stepping algorithm for stationary incompressible magnetohydrodynamics abstract for the stationary incompressible magnetohydrodynamics mhd equations we provide a new uniqueness assumption a0 and show the exponential stability of the solution then the semi implicit time stepping algorithm is used to solve the stationary mhd equations the algorithm is proved to be unconditionally stable the discrete velocity and magnetic field are bounded in l 0 l for any space and time step sizes the error estimates for the algorithm are deduced under the uniqueness conditions finally numerical experiments are performed to testify our theoretical analysis

the postprocessed mixed finite element method for the navier stokes equations refined error bounds a postprocessing technique for mixed finite element methods for the incompressible navier stokes equations is analyzed the postprocess which amounts to solving a linear stokes problem is shown to increase the order of convergence of the method to which it is applied by one unit times the logarithm of the mesh diameter in proving the error bounds some superconvergence results are also obtained contrary to previous analysis of the postprocessing technique in the present paper we take into account the loss of regularity suffered by the solutions of the navier stokes equations at the initial time in the absence of nonlocal compatibility conditions of the data as in h g heywood and r rannacher siam j numer anal 25 1988 pp 489 512 where the same hypothesis is assumed no better than fifth order convergence is achieved

multilayered shell theories accounting for layerwise mixed description part 1 governing equations a reissner mixed variational equation is employed in this paper to derive the differential governing equations of multilayered double curved shells made of orthotropic laminae in linear static cases a layerwise description is referred to by assuming two independent fields in the thickness direction for the transverse stress both shear and normal components and displacement variables in each layer interlaminar values are used as the unknown variables of the introduced expansions the continuity conditions of displacements and transverse shear and normals stresses at the interfaces between two consecutive layers referred to as c 0 z requirements have been a priori fulfilled these have been used to drive the governing equations from a layer to a multilayered level classical displacement formulations and related equivalent single layer equations have been derived for comparison purposes no assumptions have been made concerning the terms of type thickness to radii shell ratio h r donnell s shallow shell type equations are given as particular cases for all of the considered theories indicial notations and arrays have been used extensively to handle the presented developments in a concise manner numerical evaluations and comparisons to exact and other available two dimensional solutions are given in a companion paper e carrera multilayered shell theories accounting for layerwise mixed description part 2 numerical evaluations aiaa journal vol 37 no 9 1999 pp 1117 1124

viscous splitting approximation of mixed hyperbolic parabolic convection diffusion equations we first analyse a semi discrete operator splitting method for nonlinear possibly strongly degenerate convection diffusion equations due to strong degeneracy solutions can be discontinuous and are in general not uniquely determined by their data hence weak solutions satisfying an entropy condition are sought we then propose and analyse a fully discrete splitting method which employs a front tracking scheme for the convection step and a finite difference scheme for the diffusion step numerical examples are presented which demonstrate that our method can be used to compute physically correct solutions to mixed hyperbolic parabolic convection diffusion equations

layerwise theory for dynamic buckling and postbuckling of laminated composite cylindrical shells a general layerwise theory for dynamic buckling and postbuckling of imperfect multilayered composite circular cylindrical shells is introduced in contrast to the conventional displacement based layerwise theories that are more suitable for semianalytical solutions and kinematic boundary conditions the present method includes the displacement and stress components at the same time thus the method provides an efficient formulation where both kinematic and force boundary conditions can be exactly satisfied to use the exact expressions of strain and stress components the full green s strain tensor for large deflections of circular cylinders is employed then almost no simplification is made in the development of the governing equations the applied loads can be composed of various mechanical loads axial compression external pressure external fluid pressure and torsion or a combination of them accuracy and convergence of the present theory in comparison with the three dimensional elasticity approaches proposed so far are higher and more effective finally a few examples of various references that have used different theories are reexamined for comparison purpose

on an efficient second order backward difference newton scheme for mhd system abstract for the 2d 3d time dependent magnetohydrodynamics mhd system we propose a new second order backward difference formula newton scheme sbdfn which is a combination of second order backward difference approximation of the time derivative terms and newton treatment of the nonlinear terms meanwhile the finite element method is applied as spatial discretization firstly the optimal convergence of spatially semi discrete form is deduced secondly the stability and well posedness of sbdfn scheme are provided under c where c is independent of h based on these results optimal error estimates of the scheme in time are proved by negative norm technique finally numerical experiments are carried out to validate the theoretical analysis

high order direct stokes solvers with or without temporal splitting numerical investigations of their comparative properties a recently proposed direct stokes solver which decouples the velocity and pressure operators without calling for a temporal scheme is numerically analyzed in comparison first with the splitting scheme proposed by g karniadakis m israeli and s orszag in j comput phys 97 1991 pp 414 443 and with the unique grid bbb p n bbb p n 2 uzawa approach for the space accuracy and computational costs the chebyshev collocation approximation is used to analyze the spectra of the continuous temporal evolution operators and their discrete time versions from the first to the fourth order an explicit boundary condition is also involved in the proposed stokes solver and it is numerically shown that the trace of bf delta bf u on the boundary must be evaluated through its bf nabla times bf nabla times contribution only otherwise the ellipticity is lost before proceeding to the time discretization the explicit evaluation of the rotational boundary term does not prevent the first and second order in time schemes from being unconditionally stable while the schemes built at the next two higher orders are limited by the usual explicit cal o n 4 stability criterion on the time step the proposed bf u p decoupling gives this limitation a less restrictive coefficient and supplies the expected temporal orders on the whole explored range of time step sizes the effective accuracy obtained for the navier stokes 2d solutions has been measured for the re 1000 lid driven cavity problem and has been found equivalent to what is supplied by the much more expensive uzawa decoupling method

projection method i convergence and numerical boundary layers this is the first of a series of papers on the subject of projection methods for viscous incompressible flow calculations the purpose of these papers is to provide a thorough understanding of the

analysis of a sequential regularization method for the unsteady navier stokes equations the incompressibility constraint makes navier stokes equations difficult a reformulation to a better posed problem is needed before solving it numerically the sequential regularization method srm is a reformulation which combines the penalty method with a stabilization method in the context of constrained dynamical systems and has the benefit of both methods in the paper we study the existence and uniqueness for the solution of the srm and provide a simple proof of the convergence of the solution of the srm to the solution of the navier stokes equations we also give error estimates for the time discretized srm formulation

general two dimensional theory of laminated cylindrical shells the theory accounts for a desired degree of approximation of the displacements through the thickness thus accounting for any discontinuities in their derivatives at the interface of laminae geometric nonlinearity in the sense of the von karman strains is also included navier type solutions of the linear theory are presented for simply supported boundary conditions

on the asymptotics of constrained m estimation limit theorems for an m estimate constrained to lie in a closed subset of r d are given under two different sets of regularity conditions a consistent sequence of global optimizers converges under chernoff regularity of the parameter set a n consistent sequence of local optimizers converges under clarke regularity of the parameter set in either case the asymptotic distribution is a projection of a normal random vector on the tangent cone of the parameter set at the true parameter value limit theorems for the optimal value are also obtained agreeing with chernoff s result in the case of maximum likelihood with global optimizers

flow rate boundary problems for an incompressible fluid in deformable domains formulations and solution methods in this paper we consider the numerical solution of the interaction of an incompressible fluid and an elastic structure in a truncated computational domain as well known in this case there is the problem of prescribing realistic boundary data on the artificial sections when only partial data are available this problem has been investigated extensively for the rigid case in this work we start considering the compliant case by focusing on the flow rate conditions for the fluid we propose three formulations of this problem different algorithms for its numerical solution and carry out several 2d numerical simulations with the aim of comparing the performances of the different algorithms

maximum a posteriori estimators as a limit of bayes estimators maximum a posteriori and bayes estimators are two common methods of point estimation in bayesian statistics it is commonly accepted that maximum a posteriori estimators are a limiting case of bayes estimators with 0 1 loss in this paper we provide a counterexample which shows that in general this claim is false we then correct the claim that by providing a level set condition for posterior densities such that the result holds since both estimators are defined in terms of optimization problems the tools of variational analysis find a natural application to bayesian point estimation

asymptotic stochastic programs consider a stochastic program with unique solution by the notion of epi convergence in distribution in local coordinates we define the asymptotic stochastic program associated to it stochastic programs may be classified according to the type of the pertaining asymptotic program in particular three such types are studied in detail the normal shift program the wiener type program and the poisson hyperplane program conditions for the convergence to each of the three types are given

a posteriori error estimations for mixed finite element approximations to the navier stokes equations a posteriori estimates for mixed finite element discretizations of the navier stokes equations are derived we show that the task of estimating the error in the evolutionary navier stokes equations can be reduced to the estimation of the error in a steady stokes problem as a consequence any available procedure to estimate the error in a stokes problem can be used to estimate the error in the nonlinear evolutionary problem a practical procedure to estimate the error based on the so called postprocessed approximation is also considered both the semidiscrete in space and the fully discrete cases are analyzed some numerical experiments are provided

finite element approximation of the nonstationary navier stokes problem i regularity of solutions and second order error estimates for spatial discretization this is the first part of a work dealing with the rigorous error analysis of finite element solutions of the nonstationary navier stokes equations second order error estimates are proven for spatial discretization using conforming or nonconforming elements the results indicate a fluid like behavior of the approximations even in the case of large data so long as the solution remains regular the analysis is based on sharp a priori estimates for the solution particularly reflecting its behavior as t to 0 and as t to infty it is shown that the regularity customarily assumed in the error analysis for corresponding parabolic problems cannot be realistically assumed in the case of the navier stokes equations as it depends on nonlocal compatibility conditions for the data the results which are presented here are independent of such compatibility conditions which cannot be verified in practice

euler implicit explicit iterative scheme for the stationary navier stokes equations in this paper a new uniqueness assumption a2 of the solution for the stationary navier stokes equations is presented under assumption a2 the exponential stability of the solution bar u bar p for the stationary navier stokes equations is proven moreover the euler implicit explicit scheme based on the mixed finite element is applied to solve the stationary navier stokes equations finally the almost unconditionally stability is proven and the optimal error estimates uniform in time are provided for the scheme

convergence of chorin marsden product formula in the half plane consider a viscous incompressible fluid in the half plane and letu t be a solution of the navier stokes equation in this paper we prove that the product formula e t n g t n u n u0 wheree t is the euler flow g t is the heat flow and is a suitable operator describing the vorticity production due to the boundary converges uniformly tou t in the limitn

error estimates for the spectral approximation of the nonstationary stokes problem this paper analyzes a discretization of the nonstationary stokes problem with periodic nonperiodic boundary conditions fourier legendre spectral approximations are considered with respect to the space variables while the time discretization is based on the second order crank nicolson scheme optimal convergence estimates are proved for the velocity field

a defect correction method for unsteady conduction convection problems ii time discretization in this paper a fully discrete defect correction mixed finite element method mfem for solving the non stationary conduction convection problems in two dimension which is leaded by combining the back euler time discretization with the two step defect correction in space is presented in this method we solve the nonlinear equations with an added artificial viscosity term on a finite element grid and correct these solutions on the same grid using a linearized defect correction technique the stability and the error analysis are derived the theory analysis shows that our method is stable and has a good convergence property some numerical results are also given which show that this method is highly efficient for the unsteady conduction convection problems

asymptotic theory for solutions in statistical estimation and stochastic programming new techniques of local sensitivity analysis for nonsmooth generalized equations are applied to the study of sequences of statistical estimates and empirical approximations to solutions of stochastic programs consistency is shown to follow from a certain local invertibility property and asymptotic distributions are derived from a generalized implicit function theorem that characterizes asymptotic behavior in situations where estimates are subjected to constraints and estimation functionals are nonsmooth

finite element approximations of the nonstationary navier stokes problem part iv error estimates for second order time discretization this paper provides an error analysis for the crank nicolson method of time discretization applied to spatially discrete galerkin approximations of the nonstationary navier stokes equations second order error estimates are proven locally in time under realistic assumptions about the regularity of the solution for approximations of an exponentially stable solution these local error estimates are extended uniformly in time as text t to infty

asymptotic behavior of statistical estimators and of optimal solutions of stochastic optimization problems we study the asymptotic behavior of the statistical estimators that maximize a not necessarily dieren tiable criterion function possibly subject to side constraints equalities and inequalities the consistency results generalize those of wald and huber conditions are also given under which one is still able to obtain asymptotic normality the analysis brings to the fore the relationship between the problem of nding statistical estimators and that of nding the optimal solutions of stochastic optimization problems with partial information the last section is devoted to the properties of the saddle points of the associated lagrangians

the stabilized extrapolated trapezoidal finite element method for the navier stokes equations abstract we consider a fully discrete stabilized finite element method for the navier stokes equations which is unconditionally stable and has second order temporal accuracy of o k 2 hk spatial error the method involves a simple artificial viscosity stabilization of the linear system for the approximation of the new time level connected to anti diffusion of its effects at the old time level lowering the cell reynolds number of the oseen problem to o 1 the method requires only the solution of one linear system arising from an oseen problem per time step the effectiveness of the method is illustrated in several numerical experiments

error analysis of fully discrete velocity correction methods for incompressible flows a fully discrete version of the velocity correction method proposed by guermond and shen 2003 for the time dependent navier stokes equations is introduced and analyzed it is shown that when accounting for space discretization additional consistency terms which vanish when space is not discretized have to be added to establish stability and optimal convergence error estimates are derived for both the standard version and the rotational version of the method these error estimates are consistent with those by guermond and shen 2003 as far as time discretiztion is concerned and are optimal in space for finite elements satisfying the inf sup condition

optimal error estimates for a semi implicit euler scheme for incompressible fluids with shear dependent viscosities certain rheological behaviors of fluids in engineering sciences are modeled by power law ansatz with p in 1 2 in the present paper a semi implicit time discretization scheme for such fluids is proposed the main result is the optimal mathcal o k error estimate where k is the time step size our results hold in the range p in 3 2 2 in the three dimensional setting for strong solutions of the continuous problem whose existence is guaranteed under appropriate assumptions on the data the estimates are uniform with respect to the degeneracy parameter delta in 0 delta 0 of the extra stress tensor additional regularity properties of the solution of the discrete problem are proved

analysis of numerical errors in large eddy simulation we consider the question of numerical errors in large eddy simulation it is often claimed that straightforward discretization and solution using centered methods of models for large eddy motion can simulate the motion of turbulent flows with complexity independent of the reynolds number and dependence only on the resolution delta of the eddies sought this report considers this question analytically is it possible to prove error estimates for discretizations of actually used large eddy models whose error constants depend only on delta but not re we consider the most common simplest and most mathematically tractable model and the most mathematically clear discretization in two cases we prove such an error estimate and try to explain why our technique of proof fails in the most general case our analysis aims to assume as little time regularity on the true solution as possible

calculation of incompressible viscous flows by an unconditionally stable projection fem this paper investigates the numerical performance of a finite element implementation of a new incremental fractional step method to compute steady and unsteady incompressible viscous flows under general boundary conditions and using unstructured meshes a variational framework is adopted which accommodates two different spaces for representing and approximating the velocity fields calculated respectively in the viscous and inviscid phases of the method but which leads to a very simple numerical scheme in terms of only one discrete velocity field an unconditionally stable semi implicit approximation of the nonlinear term is used to eliminate any time step restriction as far as the numerical stability is concerned numerical results for five test problems in two dimensions are reported to illustrate the flexibility of the proposed method

decoupled schemes for unsteady mhd equations ii finite element spatial discretization and numerical implementation in this article a decoupled fully discrete scheme for solving 2d magnetohydrodynamics mhd equations is proposed the decoupled scheme is used for time discretization and the finite element method is used for spatial discretization firstly the almost unconditional stability t c of this scheme is established then optimal l2 and h1 error estimates of numerical solution are provided finally a numerical example is presented to confirm our theoretical results and show the high efficiency of the decoupled scheme

first order decoupled finite element method of the three dimensional primitive equations of the ocean this paper is concerned with a first order fully discrete decoupled method for solving the three dimensional 3d primitive equations of the ocean with the dirichlet boundary conditions on the side where a decoupled semi implicit scheme is used for the time discretization and the p 1 p 1 p 1 p 1 p 1 finite element for velocity pressure and density is used for the spatial discretization of these equations the h 1 l 2 h 1 optimal error estimates for the numerical solution u h n p h n theta h n and the l 2 optimal error estimate for u n h theta h n are established under the restriction of 0 h le beta 1 and 0 tau le beta 2 for some positive constants beta 1 and beta 2 moreover numerical investigations are provided to show that the first order decoupled method is of almost unconditional convergence with accuracy mathcal o h tau in the h 1 norm and mathcal o h 2 tau in the l 2 norm for solving the 3d primitive equations of the ocean numerical results are gi

stability of difference schemes for hyperbolic parabolic equations the stable difference schemes approximately solving the nonlocal boundary value problem for hyperbolic parabolic equationd 2u t dt 2 au t f t 0 t 1 u 1 au m bu l fdu t dt au t g t 1 t 0 a b 1 0 m l 1in a hilbert space h with self adjoint positive definite operator a are presented the stability estimates for the solutions of the difference schemes of the mixed type boundary value problems for hyperbolic parabolic equations are obtained the theoretical statements for the solution of these difference schemes for hyperbolic parabolic equation are supported by the results of numerical experiments

a sequential regularization method for time dependent incompressible navier stokes equations the objective of the paper is to present a method called the sequential regularization method srm for the nonstationary incompressible navier stokes equations from the viewpoint of regularization of differential algebraic equations daes and to provide a way to apply a dae method to partial differential algebraic equations pdaes the srm is a functional iterative procedure it is proved that its convergence rate is o epsilon m where m is the number of the srm iterations and epsilon is the regularization parameter the discretization and implementation issues of the method are considered in particular a simple explicit difference scheme is analyzed and its stability is proved under the usual step size condition of explicit schemes it appears that the srm formulation is new in the navier stokes context unlike other regularizations or pseudocompressibility methods in the navier stokes context the regularization parameter epsilon in the srm need not be very small and the regularized problem in the sequence may be essentially nonstiff in time direction for any epsilon hence the stability condition is independent of epsilon even for explicit time discretization numerical experiments are given to verify our theoretical results

semi discrete stabilized finite element methods for navier stokes equations with nonlinear slip boundary conditions based on regularization procedure based on the pressure projection stabilized methods the semi discrete finite element approximation to the time dependent navier stokes equations with nonlinear slip boundary conditions is considered in this paper because this class of boundary condition includes the subdifferential property then the variational formulation is the navier stokes type variational inequality problem using the regularization procedure we obtain a regularized problem and give the error estimate between the solutions of the variational inequality problem and the regularized problem with respect to the regularized parameter varepsilon which means that the solution of the regularized problem converges to the solution of the navier stokes type variational inequality problem as the parameter varepsilon longrightarrow 0 moreover some regularized estimates about the solution of the regularized problem are also derived under some assumptions about the physical data the pressure projection stabilized finite element methods are used to the regularized problem and some optimal error estimates of the finite element approximation solutions are derived

a study of transverse normal stress effect on vibration of multilayered plates and shells this paper evaluates transverse normal stress zzeffect on vibration of multilayered structures to this purpose a mixed plate model initially introduced by toledano and murakami has been extended to dynamics analysis of double curved shells these models allow both continuous interlaminar transverse shear and normal stresses as well as the zigzag form of the displacement distribution in the shell thickness directions to be modelled governing equations have been derived by employing a reissner s mixed theorem classical models on the basis of standard displacement formulations have been considered for comparison purposes the evaluations of transverse stress effects have been conducted by comparing constant linear and higher order distributions of transverse displacement components in the plate thickness directions free vibrational response of layered simply supported plates cylindrical and spherical shells made of isotropic as well as orthotropic layers has been analyzed the numerical investigation carried out and comparison with earlier results has concluded that 1 the possibility of describing interlaminar continuous transverse normal stress makes the mixed theories more attractive with respect to other available modelling 2 any refinements of classical models are meaningless unless the effects of interlaminar continuous transverse shear and normal stresses are both taken into account in a multilayered shell theory

domains of fractional powers of the stokes operator in lr spaces soit d un domaine borne de r n n 2 a frontiere lisse s soit 1 r et soit x r la fermeture dans l r d n de tous les champs vectoriels solenoidaux lisses a support compact dans d soit w r m d l espace de sobolev d ordre m et soit le laplacien on considere l operateur de stokes a r p r de domaine d a r w w r 2 d n w 5 0 mx 2 soit a un nombre reel positif pour chaque positif il y a une constante ce a telle que a r z ce a exp im z pour tout z satisfaisant a rez 0 ou a r z est la norme operateur de a r z

a penalty model for the analysis of laminated composite shells abstract a theory for laminated composite shells and a finite element model based on this theory are presented due to a flexible layer wise description of the displacement field the theory accounts for an accurate description of the shear deformation and the stretching of transverse normals to the middle surface with respect to other layer wise theories which are available in the literature the present one distinguishes itself by the use of the penalty method to enforce the perfect bonding between adjacent layers this method offers the possibility to easily obtain accurate interlaminar stresses numerical results for simply supported laminated spherical shells are presented to show the features of the proposed theory and its accuracy

reconstruction of first order shear deformation shell theory for laminated and sandwich shells a first order shear deformation shell theory fsdt based on general curvilinear coordinates has been reconstructed through a strain energy transformation based on the definition of the mean displacements the transverse shear strains that are constant through the thickness are improved to vary through the thickness it is assumed that the displacements and in surface strains of three dimensional theory are approximated by those of fsdt relationship between fsdt and three dimensional theory has been systematically established in the averaged least square sense this relationship provides the closed form recovering relations for three dimensional variables expressed in terms of fsdt relations the present approach has been applied to a higher order zigzag shell theory comparisons of displacements and stresses of laminated and sandwich shells using the present theory are made with the original fsdt and three dimensional exact solutions

viscous splitting for the unbounded problem of the navier stokes equations the viscous splitting for the exterior initial boundary value problems of the navier stokes equations is considered it is proved that the approximate solutions are uniformly bounded in the space l 0 t hs l s 2 and converge with a rate of 0 k in the space l 0 t hi q where k is the length of the time steps

on convergent schemes for diffuse interface models for two phase flow of incompressible fluids with general mass densities we are concerned with convergence results for fully discrete finite element schemes suggested in grun and klingbeil j comput phys in press doi 10 1016 j jcp 2013 10 028 they were developed for the diffuse interface model in h abels h garcke and g grun math models methods appl sci 22 2012 1150013 which describes two phase flow of immiscible incompressible viscous fluids we formulate general conditions on discretization spaces and projection operators which allow us to prove compactness of discrete solutions with respect to both time and space and which hence permit us to establish convergence of the scheme to a generalized solution we identify a simple quantitative and physical criterion to decide whether this generalized solution is in fact a weak solution in this case our analysis provides another pathway to establish existence of weak solutions to the aforementioned model in two and in three space dimensions our argument is particularly based on higher regularity results

application of finite element method in aeroelasticity the subject of this paper is the numerical simulation of the interaction of two dimensional incompressible viscous flow and a vibrating airfoil which can rotate around the elastic axis and oscillate in the vertical direction the numerical simulation consists of the finite element approximation of the navier stokes equations coupled with the system of ordinary differential equations describing the airfoil motion the arbitrary lagrangian eulerian ale formulation of the navier stokes equations stabilization the finite element discretization and coupling of both models is discussed moreover the reynolds averaged navier stokes rans system of equations together with the spallart almaras turbulence model is also discussed the computational results of aeroelastic calculations are presented and compared with the nastran code solutions

a note on the second order of accuracy difference schemes for hyperbolic parabolic equations abstract the nonlocal boundary value problem for hyperbolic parabolic equation d 2 u t d t 2 au t f t 0 t 1 d u t d t au t g t 1 t 0 u 1 u 0 1 0 1 in a hilbert space h is considered the second order of accuracy difference schemes approximately solving this boundary value problem are presented the stability estimates for the solution of these difference schemes are established in applications the stability estimates for the solutions of the difference schemes of the mixed type boundary value problems for hyperbolic parabolic equations are obtained the theoretical statements for the solution of these difference schemes for hyperbolic parabolic equation are supported by the results of numerical experiments

efficient higher order shell theory for laminated composites with multiple delaminations a higher order zigzag theory has been developed for laminated composite shells with multiple delaminations general tensor based formulation is developed for arbitrary curved shell with exact geometric description a laminated shell theory with multiple delaminations for general lamination configurations is obtained by superposing a cubic varying displacement on a zigzag linearly varying displacement the von karman nonlinearity is included in the formulation for the potentialities of addressing problems requiring geometric nonlinearity such as large deflection and postbuckling problems when top and bottom surface transverse shear stress free conditions and interface transverse shear continuity conditions including delamination interfaces are imposed the displacement of the minimal degrees of freedom is obtained the proposed displacement field can systematically handle the number shape size and locations of the delaminations through the variational principle equilibrium equations and variationally consistent boundary conditions are obtained to assess the accuracy and efficiency of the present theory the linear buckling problem of cylindrical shell with multiple delaminations has been analyzed the higher order zigzag theory should work as an efficient tool for analyzing the behavior of composite laminated shells with multiple delaminations