stability and convergence of two level iterative methods for the stationary incompressible magnetohydrodynamics with different reynolds numbers abstract in the paper we develop some two level finite element iterative methods and use these methods to solve the stationary incompressible magnetohydrodynamics mhd with different reynolds numbers under some different uniqueness conditions firstly we use the stokes type iterative method newton type iterative method and oseen type iterative method of m times on a coarse mesh with mesh size h and then solve a linear problem with the stokes type newton type and oseen type correction of one time on a fine grid with mesh sizes h h furthermore we analyze the uniform stability and convergence of these methods with respect to reynolds numbers mesh sizes h and h and iterative times m finally the stationary incompressible mhd equations with large reynolds number are researched by the one level finite element method based on the oseen type iteration on a fine mesh under a weak uniqueness condition

on a two grid finite element scheme combined with crank nicolson method for the equations of motion arising in the kelvin voigt model in this paper we present a second order accurate crank nicolson scheme applied to a system arising from galerkin finite element two grid approximations to the equations of motion described by the kelvin voigt viscoelastic fluid flow model optimal error estimates in l l 2 norm and l h 1 norm for velocity and l l 2 norm for pressure are established these estimates are shown to be uniform in time numerical results are provided to verify the theoretical estimates

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

optimal error bounds for two grid schemes applied to the navier stokes equations abstract we consider two grid mixed finite element schemes for the spatial discretization of the incompressible navier stokes equations a standard mixed finite element method is applied over the coarse grid to approximate the nonlinear navier stokes equations while a linear evolutionary problem is solved over the fine grid the previously computed galerkin approximation to the velocity is used to linearize the convective term for the analysis we take into account the lack of regularity of the solutions of the navier stokes equations at the initial time in the absence of nonlocal compatibility conditions of the data optimal error bounds are obtained

time domain finite difference and finite element methods for maxwell equations in complex media extensions of finite difference time domain fdtd and finite element time domain fetd algorithms are reviewed for solving transient maxwell equations in complex media also provided are a few representative examples to illustrate the modeling capabilities of fdtd and fetd for complex media the term complex media refers here to media with dispersive bi anisotropic inhomogeneous and or nonlinear properties present in the constitutive tensors

a computational approach to approximate and plausible reasoning with applications to expert systems the intended purpose of this paper is twofold proposing a common basis for the modeling of uncertainty and imprecision and discussing various kinds of approximate and plausible reasoning schemes in this framework together with probability different kinds of uncertainty measures credibility and plausibility functions in the sense of shafer possibility measures in the sense of zadeh and the dual measures of necessity sugeno s g fuzzy measures are introduced in a unified way the modeling of imprecision in terms of possibility distribution is then presented and related questions such as the measure of the uncertainty of fuzzy events the probability and possibility qualification of statements the concept of a degree of truth and the truth qualification of propositions are discussed at length deductive inference from premises weighted by different kinds of measures by uncertainty or by truth values in the framework of various multivalued logics is fully investigated then deductive inferences from imprecise or fuzzy premises are dealt with patterns of reasoning where both uncertainty and imprecision are present are also addressed the last section is devoted to the combination of uncertain or imprecise pieces of information given by different sources on the whole this paper is a tentative survey of quantitative approaches in the modeling of uncertainty and imprecision including recent theoretical proposals as well as more empirical techniques such as the ones developed in expert systems such as mycin or prospector the management of uncertainty and imprecision in reasoning patterns being a key issue in artificial intelligence

using fuzzy sets in operational research abstract during the last four decades a large number of mainly mathematical tools have been developed in and for operational research they are primarily devices to find optimal solutions after a problem has been modelled formally hardly any new tools have been developed for proper modelling of problems which are not well structured and which can not easily be cast into classical mathematical models fuzzy set theory as proposed by l zadeh in 1965 promises to bridge part of this gap it is much more adaptable to different problem structures and better suited to model human evaluation and decision making processes than classical mathematics this paper will primarily focus attention on four questions why fuzzy sets in or what are fuzzy sets where can they be used and how are they being used

on a higher order accurate fully discrete galerkin approximation to the navier stokes equations we consider approximating the solution of the initial and boundary value problem for the navier stokes equations in bounded twoand three dimensional domains using a nonstandard galerkin finite element method for the space discretization and the third order accurate three step backward differentiation method coupled with extrapolation for the nonlinear terms for the time stepping the resulting scheme requires the solution of one linear system per time step plus the solution of five linear systems for the computation of the required initial conditions all these linear systems have the same matrix the resulting approximations of the velocity are shown to have optimal rate of convergence in l2 under suitable restrictions on the discretization parameters of the problem and the size of the solution in an appropriate function space

a spectral element method for the navier stokes equations with improved accuracy we present a procedure to enhance the accuracy of spectral element methods for the evolutionary navier stokes equations the approximations to the velocity and the pressure obtained using a spectral element method are postprocessed by solving a discrete stokes problem only once when the time evolution is completed the analysis shows that this procedure improves the order of convergence for the approximations to the velocity and the pressure the choice of a discrete pressure space endowed with an inf sup condition independent of the discretization parameter leads to optimal error bounds for the pressure

two level method and some a priori estimates in unsteady navier stokes calculations a two level method proposed for quasielliptic problems is adapted in this paper to the simulation of unsteady incompressible navier stokes flows the method requires a solution of a nonlinear problem on a coarse grid and a solution of linear symmetric problem on a fine grid the scaling between these two grids is superlinear approximation stability and convergence aspects of a fully discrete scheme are considered stability properties of the two level scheme are compared with those for a commonly used semi implicit scheme some new estimates are also proved for the latter

an approach for solving a fuzzy multiobjective programming problem abstract in this paper we present a new approach based on the nearest interval approximation operator for dealing with a multiobjective programming problem with fuzzy valued objective functions by the way we have established a karush kuhn tucker k k t kind of pareto optimality conditions for the resulting interval multiobjective program to this end we made use of gh differentiability of involved interval valued functions two algorithms play a pivotal role in the proposed method the first one returns a nearest interval approximation to a given fuzzy number the other one makes use of k k t conditions to deliver a pareto optimal solution of the above mentioned resulting interval program

distinction between several subsets of fuzzy measures abstract in this paper we place several fuzzy measure subsets in relation one with the other the subsets under study are those corresponding to the definitions of probability measure sugeno s g measure shafer s belief function and zadeh s possibility measure we study the intersection of these subsets and we show the particular role of dirac s measures in this comparison we limit ourself to the case of mappins whose domain is the collection of all subsets of a finite set finally the obtained partial results are summarized in only one figure which shoul clarify the specificity of each of the above definitions

multi level spectral galerkin method for the navier stokes problem i spatial discretization a multi level spectral galerkin method for the two dimensional non stationary navier stokes equations is presented the method proposed here is a multiscale method in which the fully nonlinear navier stokes equations are solved only on a low dimensional space subsequent approximations are generated on a succession of higher dimensional spaces j 2 j by solving a linearized navier stokes problem around the solution on the previous level error estimates depending on the kinematic viscosity 0 1 are also presented for the j level spectral galerkin method the optimal accuracy is achieved when we demonstrate theoretically that the j level spectral galerkin method is much more efficient than the standard one level spectral galerkin method on the highest dimensional space

postprocessing finite element methods for the navier stokes equations the fully discrete case an accuracy enhancing postprocessing technique for finite element discretizations of the navier stokes equations is analyzed the technique had been previously analyzed only for semidiscretizations and fully discrete methods are addressed in the present paper we show that the increased spatial accuracy of the postprocessing procedure is not affected by the errors arising from any convergent time stepping procedure further refined bounds are obtained when the time stepping procedure is either the backward euler method or the two step backward differentiation formula

analysis of two grid method for semi linear elliptic equations by new mixed finite element scheme this article considers two grid stable mixed finite element method based on the less regularity of flux velocity in practice for the semi linear elliptic equations approximated by the p 0 p 1 velocity pressure pair which satisfies the inf sup condition this method involves solving a small semi linear system on a coarse mesh with mesh size h and a linear system problem on a fine mesh with mesh size h o h which can still maintain the asymptotically optimal accuracy it provides an approximate solution with the convergence rate of the same order as the usual mixed finite element solution solving the semi linear elliptic equations on a fine mesh hence the two grid scheme can reduce the computational cost finally numerical tests confirm the theoretical results of the present method

a simplified two level method for the steady navier stokes equations abstract a simplified two level method for solving the steady two dimensional incompressible navier stokes equations is investigated in this article the simplified two level method involves solving one small nonlinear navier stokes problem on the coarse mesh and one linear stokes problem hence linear with positive definite symmetric part on the fine mesh the uniform optimal error estimates with respect to large for the standard galerkin method and the two level method are proven finally some numerically tests to confirm the theoretical results of the simplified two level method and the standard galerkin method are provided

nonlinear galerkin methods this article presents a new method of integrating evolution differential equations the non linear galerkin method that is well adapted to the long term integration of such equations while the usual galerkin method can be interpreted as a projection of the considered equation on a linear space the methods considered here are related to the projection of the equation on a nonlinear manifold from the practical point of view some terms have been identified as small and sometimes but not always disregarded

effective preconditioning of uzawa type schemes for a generalized stokes problem the schur complement of a model problem is considered as a preconditioner for the uzawa type schemes for the generalized stokes problem the stokes problem with the additional term alpha boldmath u in the motion equation the implementation of the preconditioned method requires for each iteration only one extra solution of the poisson equation with neumann boundary conditions for a wide class of 2d and 3d domains a theorem on its convergence is proved in particular it is established that the method converges with a rate that is bounded by some constant independent of alpha some finite difference and finite element methods are discussed numerical results for finite difference mac scheme are provided

bi fdtd a novel finite difference time domain formulation for modeling wave propagation in bi isotropic media this paper presents a newly developed finite difference time domain fdtd technique referred to as bi fdtd for modeling electromagnetic wave interactions with bi isotropic bi media the theoretical foundation for the bi fdtd method will be developed based on a wavefield decomposition the main advantage of this approach is that the two sets of wavefields are uncoupled and can be viewed as propagating in an equivalent isotropic medium which makes it possible to readily apply conventional fdtd analysis techniques the bi fdtd scheme will also be extended to include the dispersive nature of chiral media an important subclass of bi isotropic media this extension represents the first of its kind in the fdtd community validations of this new model are demonstrated for a chiral half space and a chiral slab

an extension to possibilistic linear programming abstract this paper proposes an extension to the solution of possibilistic linear programming problems with fuzzy number parameters introduced by buckley fuzzy sets and systems31 1989 329 341 the method relies on cuts of the fuzzy number parameters to generate a succession of pairs of classical linear programs the solutions of these linear programs produce a convex possibility distribution of the optimal objective values with the corresponding levels of the decision variables the approach is illustrated with an imprecise water quality management problem

the postprocessed mixed finite element method for the navier stokes equations a postprocessing technique for mixed finite element methods for the incompressible navier stokes equations is studied the technique was earlier developed for spectral and standard finite element methods for dissipative partial differential equations the postprocessing amounts to solving a stokes problem on a finer grid or higher order space once the time integration on the coarser mesh is completed the analysis presented here shows that this technique increases the convergence rate of both the velocity and the pressure approximations numerical experiments are presented that confirm both this increase in the convergence rate and the corresponding improvement in computational efficiency

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

modelling of wave propagation in wire media using spatially dispersive finite difference time domain method numerical aspects the finite difference time domain fdtd method is applied for modelling of wire media as artificial dielectrics both frequency dispersion and spatial dispersion effects in wire media are taken into account using the auxiliary differential equation method according to the authors knowledge this is the first time when the spatial dispersion effect is considered in the fdtd modelling the stability of developed spatially dispersive fdtd formulations is analyzed through the use of von neumann method combined with the routh hurwitz criterion the results show that the conventional stability courant limit is preserved using standard discretization scheme for wire media modelling flat sub wavelength lenses formed by wire media are chosen for validation of proposed spatially dispersive fdtd formulation results of the simulations demonstrate excellent sub wavelength imaging capability of the wire medium slabs the size of the simulation domain is significantly reduced using the modified perfectly matched layer mpml which can be placed in close vicinity of the wire medium it is demonstrated that the reflections from the mpml wire medium interface are less than 70 db that lead to dramatic improvement of convergence compared to conventional simulations

finite element approximation of the nonstationary navier stokes problem part iii smoothing property and higher order error estimates for spatial discretization this paper continues our error analysis of finite element galerkin approximation of the nonstationary navier stokes equations optimal order error estimates both local and global are derived for higher order finite elements under appropriate assumptions about the smoothness and stability of the solution these assumptions take into account the loss of regularity at t 0 that one generally has to expect in the absence of higher order nonlocal compatibility conditions for the data of the problem

a second order accuracy for a full discretized time dependent navier stokes equations by a two grid scheme we study a second order two grid scheme fully discrete in time and space for solving the navier stokes equations the two grid strategy consists in discretizing in the first step the fully non linear problem in space on a coarse grid with mesh size h and time step t and in the second step in discretizing the linearized problem around the velocity u h computed in the first step in space on a fine grid with mesh size h and the same time step the two grid method has been applied for an analysis of a first order fully discrete in time and space algorithm and we extend the method to the second order algorithm this strategy is motivated by the fact that under suitable assumptions the contribution of u h to the error in the non linear term is measured in the l 2 norm in space and time and thus has a higher order than if it were measured in the h 1 norm in space we present the following results if h 2 h 3 t 2 then the global error of the two grid algorithm is of the order of h 2 the same as would have been obtained if the non linear problem had been solved directly on the fine grid

postprocessing fourier galerkin method for the navier stokes equations a full discrete two level postprocessing fourier galerkin scheme for the unsteady navier stokes equations with a periodic boundary condition is proposed in this paper by defining a new projection the interaction between the large and small eddies is reflected by the associated space splitting to some extent therefore a weakly coupled system of the large and small eddies is obtained stability and error estimates for the weakly coupled two level scheme are established in the paper the proposed scheme is an effective algorithm because nonlinearity is treated only in the coarse level subspace by solving the coarse level standard galerkin equation and the matrix of the linear algebraic equations arising in each fine level time stepping is a sparse matrix especially for large scale computations as a result the new scheme saves a lot of cpu time and memory in comparison with the standard fourier galerkin method for deriving an approximation of prescribed accuracy

investigations on two kinds of two level stabilized finite element methods for the stationary navier stokes equations abstract this article considers two kinds of two level stabilized finite element methods based on local gauss integral technique for the two dimensional stationary navier stokes equations approximated by the lowest equal order elements which do not satisfy the inf sup condition the two level methods consist of solving a small non linear system on the coarse mesh and then solving a linear system on the fine mesh the error analysis shows that the two level stabilized finite element methods provide an approximate solution with the convergence rate of the same order as the usual stabilized finite element solution solving the navier stokes equations on a fine mesh for a related choice of mesh widths therefore the two level methods are of practical importance in scientific computation finally the performance of two kinds of two level stabilized methods are compared in efficiency and precision aspects by a series of numerical experiments the conclusion is that simple two level stabilized method is a viable choice for the lowest equal order approximations of the stationary navier stokes problem

finite difference and pseudospectral time domain methods applied to backward wave metamaterials backward wave bw materials that have simultaneously negative real parts of their electric permittivity and magnetic permeability can support waves where phase and power propagation occur in opposite directions these materials were predicted to have many unusual electromagnetic properties among them amplification of the near field of a point source which could lead to the perfect reconstruction of the source field in an image j pendry phys rev lett vol 85 pp 3966 2000 often systems containing bw materials are simulated using the finite difference time domain technique we show that this technique suffers from a numerical artifact due to its staggered grid that makes its use in simulations involving bw materials problematic the pseudospectral time domain technique on the other hand uses a collocated grid and is free of this artifact it is also shown that when modeling the dispersive bw material the linear frequency approximation method introduces error that affects the frequency of vanishing reflection while the auxiliary differential equation the z transform and the bilinear frequency approximation method produce vanishing reflection at the correct frequency the case of vanishing reflection is of particular interest for field reconstruction in imaging applications

nonlinear galerkin methods the finite elements case with the increase in the computing power and the advent of supercomputers the approximation of evolution equations on large intervals of time is emerging as a new type of numerical problem in this article we consider the approximation of evolution equations on large intervals of time when the space discretization is accomplished by finite elements the algorithm that we propose called the nonlinear galerkin method stems from the theory of dynamical systems and amounts to some approximation of the attractor in the discrete finite elements space essential here is the utilization of incremental unknown which is accomplished in finite elements by using hierarchical bases beside a detailed description of the algorithm the article includes some technical results on finite elements spaces and a full study of the stability and convergence of the method

stability of higher order hood taylor methods the stability of a higher order hood taylor method for the approximation of the stationary stokes equations using continuous piecewise polynomials of degree 3 to approximate velocities and continuo

the postprocessing galerkin and nonlinear galerkin methods a truncation analysis point of view we revisit the postprocessing algorithm and give a justification from a classical truncation analysis point of view we assume a perturbation expansion for the high frequency mode component of solutions to the underlying equation keeping terms to certain orders we then generate approximate systems which correspond to numerical schemes we show that the first two leading order methods are in fact the postprocessed galerkin and postprocessed nonlinear galerkin methods respectively hence postprocessed galerkin is a natural leading order method more natural than the standard galerkin method for approximating solutions of parabolic dissipative pdes the analysis is presented in the framework of the two dimensional navier stokes equation nse however similar analysis may be done for any parabolic dissipative nonlinear pde the truncation analysis is based on asymptotic estimates in time for the low and high mode components we also introduce and investigate an alternative postprocessing scheme whic

fuzzy optimization an appraisal abstract this paper takes a general look at core ideas that make up the burgeoning body of fuzzy mathematical programming emphasizing the methodological view although fuzzy mathematical programming has enjoyed a rapidly increasing acceptance within the scientific community some technical hurdles exist to hinder a unanimity reasons for this as well as possible ways for improvement are also discussed

postprocessing the linear finite element method we extend the idea of the postprocessing galerkin methods for dissipative evolution equations to the case of the linear finite element method the postprocessing technique has been developed earlier for spectral methods and for higher order finite element methods the analysis shows that this procedure improves the order of convergence of the piecewise linear galerkin finite element approximation in the h1 norm we show by means of numerical experiments that there is no improvement in the order of convergence in the l2 norm

a multilevel mesh independence principle for the navier stokes equations multilevel finite element discretization methods for the navier stokes equations are considered in contrast to usual multilevel methods a superlinear scaling of the consecutive meshwidths h j 1 mathcal o h j alpha j is used on the coarsest mesh the discretized system is solved which is small and nonlinear on each subsequent mesh only one or two newton correction steps are performed so that only one or two larger linear systems are solved the scalings of the meshwidths that lead to optimal accuracy of the approximate solution in both the h 1 and l 2 norm are investigated an error analysis of independent interest is also presented for the basic finite element method using the conservation form of the nonlinear term this formulation leads to a nonlinear system whose nonlinearity is easier to resolve than the systems arising when either the convective or explicitly skew symmetric corms are used for the nonlinear term

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

water quality management with imprecise information abstract the application of fuzzy set and possibility theories for the representation of imprecise information in water quality management problems is investigated imprecision is associated with fuzziness which is non statistical in nature and refers to the absence of sharp boundaries in the information imprecise parameters in water quality decision making can be represented by possibility distributions defining maximum achievable probabilities the corresponding possibilistic programming problem is viewed as an alternative to the stochastic one where the parameters are modeled as fuzzy variables instead of random variables the possibilistic problem can be solved through a succession of classical linear programs the resulting possibility distributions of the objective value provide a possibilistic assessment of the risk based on possibility levels

nonlinear galerkin methods and mixed finite elements two grid algorithms for the navier stokes equations a nonlinear galerkin method using mixed finite elements is presented for the two dimensional incompressible navier stokes equations the scheme is based on two finite element spaces x h and x h for the approximation of the velocity defined respectively on one coarse grid with grid size h and one fine grid with grid size h ll h and one finite element space m h for the approximation of the pressure nonlinearity and time dependence are both treated on the coarse space we prove that the difference between the new nonlinear galerkin method and the standard galerkin solution is of the order of h 2 both in velocity h 1 omega and pressure l 2 omega norm we also discuss a penalized version of our algorithm which enjoys similar properties

a finite element variational multiscale method for incompressible flow in this paper we present a numerical scheme prove stability existence of solutions uniqueness and convergence of the incompressible navier stokes equations it has the advantage of being defined from strictly algebraic considerations a significant feature of the present method is that the structure of the stabilization term based on the multiscale enrichment and derived from the navier stokes problem itself ample numerical experiments are carried out to confirm the theory and illustrate the effectiveness of the scheme on incompressible fluid

fuzzy sets as a basis for a theory of possibility the theory of possibility described in this paper is related to the theory of fuzzy sets by defining the concept of a possibility distribution as a fuzzy restriction which acts as an elastic constraint on the values that may be assigned to a variable more specifically if f is a fuzzy subset of a universe of discourse u u which is characterized by its membership function f then a proposition of the form x is f where x is a variable taking values in u induces a possibility distribution tx which equates the possibility of x taking the value u to f u the compatibility of u with f in this way x becomes a fuzzy variable which is associated with the possibility distribution tx in much the same way as a random variable is associated with a probability distribution in general a variable may be associated both with a possibility distribution and a probability distribution with the weak connection between the two expressed as the possibility probability consistency principle a thesis advanced in this paper is that the imprecision that is intrinsic in natural languages is in the main possibilistic rather than probabilistic in nature thus by employing the concept of a possibility distribution a proposition p in a natural language may be translated into a procedure which computes the probability distribution of a set of attributes which are implied by p several types of conditional translation rules are discussed and in particular a translation rule for propositions of the form x is f is possible where is a number in the interval 0 1 is formulated and illustrated by examples

error estimates on a new nonlinear galerkin method based on two grid finite elements a new nonlinear galerkin method based on finite element discretization is presented in this paper for semilinear parabolic equations the new scheme is based on two different finite element spaces defined respectively on one coarse grid with grid size h and one fine grid with grid size h ll h nonlinearity and time dependence are both treated on the coarse space and only a fixed stationary equation needs to be solved on the fine space at each time with linear finite element discretizations it is proved that the difference between the new nonlinear galerkin solution and the standard galerkin solution in h 1 omega norm is of the order of h 3

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

joint solution to fuzzy programming problems we propose a new solution concept for fuzzy programming problems it is based on our new method of solving fuzzy equations 101 for simplicity we discuss in detail only fuzzy linear programming in this paper we define and obtain the basic properties of the joint solution a fuzzy vector in rn and the optimal value of the objective function a fuzzy number three examples are presented illustrating these concepts

error estimates of two level finite element method for smagorinsky model two level finite element method for simulating smagorinsky model in large eddy simulation is investigated in this two level algorithm a linearized discrete problem is solved by using the variational multiscale method on the coarse mesh corresponding to newton linearization method the linearized smagorinsky model is solved on the fine mesh the error estimates derived in this paper implies that by choosing appropriate mesh sizes and the radius of the spatial filter used in smagorinsky model the two level method provides the optimal convergence rates for the velocity in h1 norm and the pressure in l2 norm meanwhile the two different numerical experiments are given to support the optimal convergence rates and the high efficiency of two level algorithm

galerkin method and optimal error algorithms in hilbert spaces we consider operator equations of the formlu f wherel belongs to the class of linear bounded by a constantm and coercive with a constantm operators from a hilbert spacev onto its dualv andf belongs to a hilbert spacew v we study optimality of the galerkin methodp n lu n p n f whereu n v n v n is subspace ofv p n is the orthogonal projector ontov n andp n is dual top n we show that the galerkin method is quasi optimal independently of the choice of the subspacev n and the spacew ifm m in the casem m optimality of the method depends strongly on the choice ofv n andw therefore we define a new algorithm which is always optimal independently of the choice ofv n andw and relations betweenm andm

a two level method with backtracking for the navier stokes equations we consider a two level method for resolving the nonlinearity in finite element approximations of the equilibrium navier stokes equations the method yields l2 and h1 optimal velocity approximations and an l2 optimal pressure approximation the two level method involves solving one small nonlinear coarse mesh system one oseen problem hence linear with positive definite symmetric part on the fine mesh and one linear correction problem on the coarse mesh the algorithm we study produces an approximate solution with the optimal asymptotic in h accuracy for any fixed reynolds number we do not consider the behavior of the error for fixed h as re rightarrow infty i e for flows in transition to turbulence

a postprocessed galerkin method with chebyshev or legendre polynomials we present an approximate inertial manifold based postprocess to enhance chebyshev or legendre spectral galerkin methods we prove that the postprocess improves the order of convergence of the galerkin solution yielding the same accuracy as the nonlinear galerkin method numerical experiments show that the new method is computationally more efficient than galerkin and nonlinear galerkin methods new approximation results for chebyshev polynomials are presented

on the rate of convergence of the nonlinear galerkin methods in this paper we provide estimates to the rate of convergence of the nonlinear galerkin approximation method in particular and by means of an illustrative example we show that the nonlinear galerkin method converges faster than the usual galerkin method

stability analysis of the nonlinear galerkin method our object in this article is to describe some numerical schemes for the approximation of nonlinear evolution equations and to study the stability of the schemes spatial discretization can be performed by either spectral or pseudospectral methods finite elements or finite differences time discretization is done by two level schemes partly or fully explicit the algorithms that we present stem from the study of the evolution equations from the dynamical systems point of view they are based on a differentiated treatment of the small and large wave lengths and they are particularly adapted to the integration of such equations on large intervals of time introduction the long time integration of dissipative evolution equations has emerged as a new type of problem in numerical analysis and its appearance is due to the considerable increase in computing power that we have seen during the last years much effort has been devoted in the past to the approximation of evolution equations on finite intervals of time or in the dynamically simple case where the solution converges to a steady stationary state as time goes to infinity however in many physically relevant situations the solutions to a dissipative evolution equation do not converge to a stationary solution they rather remain time dependent and they converge to a compact attractor that encompasses the nonlinear dynamics hence and whether this is addressed implicitly or explicitly the numerical integration of the evolution equation is then closely related to the approximation of the attractor that may be a complicated set even a fractal inertial algorithms are new integration algorithms that are based on the approximation of the attractor by smooth simple manifolds namely the approximate inertial manifolds these algorithms will be recalled below but at this point let us point out that one of their main features is a differentiated treatment of small and large wave lengths an inertial manifold is an exact law of dependence of small wave lengths in terms of the large ones the small wave received april 3 1990 1980 mathematics subject classification 1985 revision primary 65f10 65f35 65n20 35a40 1991 american mathematical society 0025 5718 91 1 00 25 perpage

possibility and necessity in optimization abstract we employ possibility and necessity in fuzzy programming and in other optimization problems that traditionally use random variables the possibility approach to multiobjective programming produces the same results as fuzzy programming a detailed example is presented to illustrate the use of possibility and necessity in optimization problems that usually employ random variables to describe uncertainty

fdtd for nth order dispersive media previously a method for applying the finite difference time domain fdtd method to dispersive media with complex permittivity described by a function with a single first order pole was presented this method involved the recursive evaluation of a discrete convolution and was therefore relatively efficient in this work the recursive convolution approach is extended to media with dispersions described by multiple second order poles the significant change from the first order implementation is that the single backstore variable for each second order pole is complex the approach is demonstrated for a pulsed plane wave incident on a medium with a complex permittivity described by two second order poles and excellent agreement is obtained with the exact solution

fuzzy programming and the pareto optimal set it is shown how fuzzy programming using either the min or product operator may be used to generate the whole pareto optimal set for nonlinear concave or convex multiobjective programming problems also a new solution procedure for a fuzzy program when the min operator is employed is presented

postprocessing the galerkin method a novel approach to approximate inertial manifolds a postprocess of the standard galerkin method for the discretization of dissipative equations is presented the postprocessed galerkin method uses the same approximate inertial manifold phi app to approximate the high wave number modes of the solution as in the nonlinear galerkin nlg method however in this postprocessed galerkin method the value of phi app is calculated only once after the time integration of the standard galerkin method is completed contrary to the nlg in which phi app evolves with time and affects the time evolution of the lower wave number modes the postprocessed galerkin method which is much cheaper to implement computationally than the nlg is shown in the case of fourier modes to possess the same rate of convergence accuracy as the simplest version of the nlg which is based on either the foias manley temam approximate inertial manifold or the euler galerkin approximate inertial manifold this is proved for some problems in one and two spatial dimensions including the navier stokes equations under periodic boundary conditions the advantages of postprocessing that we present here apply not only to the standard galerkin method but also to the computationally more efficient pseudospectral method

stability and phase error analysis of fd td in dispersive dielectrics four fd td extensions for the modeling of pulse propagation in debye or lorentz dispersive media are analyzed through studying the stability and phase error properties of the coupled difference equations corresponding to maxwell s equations and to the equations for the dispersion for good overall accuracy the author shows that all schemes should be run at their courant stability limit and that the timestep should finely resolve the medium timescales particularly for debye schemes it should be at least spl delta t 10 sup spl minus 3 spl tau while for lorentz schemes it should be spl delta t 10 sup spl minus 2 spl tau where spl tau is a typical medium relaxation time a numerical experiment with a debye medium confirms this the author has determined that two of the discretizations for debye media are totally equivalent in the lorentz medium case the author establishes that the method that uses the polarization differential equation to model dispersion is stable for all wavenumbers and that the method using the local in time constitutive relation is weakly unstable for modes with wavenumber k such that k spl delta x spl pi 2

an approximate inertial manifolds approach to postprocessing the galerkin method for the navier stokes equation in a recent paper we have introduced a postprocessing procedure for the galerkin method for dissipative evolution partial differential equations with periodic boundary conditions the postprocessing technique uses approximate inertial manifolds to approximate the high modes the small scale components in the exact solutions in terms of the galerkin approximations which in this case play the role of the lower modes large scale components this procedure can be seen as a defect correction technique but contrary to standard procedures the correction is computed only when the time evolution is completed here we extend these results to more realistic boundary conditions specifically we study in detail the two dimensional navier stokes equations subject to homogeneous nonslip dirichlet boundary conditions we also discuss other equations such as reaction diffusion systems and the cahn hilliard equations

frequency dependent fdtd methods using z transforms the frequency dependent finite difference time domain method fd sup 2 td method has been shown to be capable of correctly calculating electromagnetic propagation through media whose dielectric properties are frequency dependent however as researchers search for more elaborate applications the formulation of the fd sup 2 td methods becomes more complex in this work the mathematics of the fd sup 2 td method is developed using z transform theory this has the advantages of presenting a clearer formulation and allowing researchers to draw on the literature of systems analysis and signal processing disciplines

long time unconditional stability of a two level hybrid method for nonstationary incompressible navier stokes equations abstract this paper presents a two level maccormack rapid solver method for solving a 2d time dependent navier stokes problem that models an incompressible fluid flow a decoupling approach based on the coupling term via spatial extrapolation is proposed for devising decoupled marching algorithms for the nonstationary incompressible model in this hybrid method an explicit maccormack scheme provides a h 1 optimal velocity approximation and a l 2 optimal pressure approximation by solving a nonlinear navier stokes problem on a coarse grid with mesh size h while an implicit crank nicolson algorithm consists in dealing with the fully discrete linear generalized stokes problem on a fine grid with mesh size h h the theoretical result suggests that our method is unconditionally stable over long time intervals while the numerical evidences both confirm the theoretic analysis and show that the algorithm is convergent the tests also suggest that our method is both cheaper and faster than the two level finite element galerkin scheme

a novel two grid method for semilinear elliptic equations a new finite element discretization technique based on two coarse and fine subspaces is presented for a semilinear elliptic boundary value problem the solution of a nonlinear system on the fine space is reduced to the solution of two small one linear and one nonlinear systems on the coarse space and a linear system on the fine space it is shown both theoretically and numerically that the coarse space can be extremely coarse and still achieve asymptotically optimal approximation as a result the numerical solution of such a nonlinear equation is not significantly more expensive than the solution of one single linearized equation

modality constrained programming problems a unified approach to fuzzy mathematical programming problems in the setting of possibility theory abstract in this paper fuzzy mathematical programming problems are formulated based on the idea analogous with the chance constrained programming problem the difference in meaning between the ambiguity of the coefficients and that of the decision maker s preference is emphasized the constraints with fuzzy coefficients are treated as the restriction that should be satisfied properly rather than perfectly the objective functions with fuzzy coefficients are treated variously depending on the interpretations i e the optimization of the modalities the optimization of the fractile the minimization of the ambiguity and so forth the deterministic equivalent constraints and the deterministic equivalent problems are shown when the constraints and the objective functions are linear a numerical example is given to illustrate the proposed formulations

ranking fuzzy numbers in the setting of possibility theory the arithmetic manipulation of fuzzy numbers or fuzzy intervals is now well understood equally important for application purposes is the problem of ranking fuzzy numbers or fuzzy intervals which is addressed in this paper a complete set of comparison indices is proposed in the framework of zadeh s possibility theory it is shown that generally four indices enable one to completely describe the respective locations of two fuzzy numbers moreover this approach is related to previous ones and its possible extension to the ranking of n fuzzy numbers is discussed at length finally it is shown that all the information necessary and sufficient to characterize the respective locations of two fuzzy numbers can be recovered from the knowledge of their mutual compatibilities

two level picard and modified picard methods for the navier stokes equations abstract iterative methods of picard type for the navier stokes equations are known to converge only for quite small reynolds numbers however we study methods involving just one such iteration at general reynolds numbers for the initial approximation a coarse mesh of width h0 is used the corrected approximation is computed by just one picard or modified picard step on a fine mesh of width h1 for example h1 may be of order o h 0 2 when linear velocity elements are used the resulting method requires the solution of a small system of nonlinear equations on the coarse mesh and only one larger linear system on the fine mesh this two level picard method is proven to converge for fixed reynolds number as h 0 further the fine mesh solution satisfies a quasi optimal error bound the error constants grow as re as for the usual finite element method one very heuristic explanation why one step of the divergent picard method might work when beginning with a coarse mesh approximation is that the terms neglected involve lower order derivatives thus they are approximated with higher accuracy on the coarse mesh this is linked with a smoothing property of the fine mesh step

on a hierarchical basis multilevel method with nonconforming p1 elements we develop a hierarchical basis multilevel method for symmetric second order elliptic boundary value problems in two dimensional polygonal domains based on nonconforming p1 triangular elements the main result is that the condition number of the hierarchical discretization is bounded byo k wherek is the number of refinement levels for comparison yserentant s result yields the sharp boundo k 2 for the corresponding hierarchical discretization with conforming linear elements