a semi implicit numerical scheme for reacting flow an additive semi implicit projection scheme for the simulation of unsteady combustion in two dimensions is constructed the scheme relies on a zero mach number formulation of the compressible conservation equations with detailed chemistry the governing equations are discretized in space using second order differences and integrated in time using a semi implicit approach time integration of the evolution equations for species mass fraction thermodynamic pressure and density is performed using a semi implicit nonsplit scheme that combines a second order predictor corrector treatment of convection and diffusion terms and a stiff integrator for the reaction source terms meanwhile the momentum equations are integrated using a second order projection scheme the projection scheme is based on a predictor corrector approach that couples the evolution of the velocity and density fields in order to stabilize computations of reacting flows with large density variations a pressure poisson equation is inverted following both the predictor and corrector steps using a fast solver the advantages of the stiff integration of reaction source terms are analyzed by comparing the performance of the scheme to that of a predictor corrector scheme in which reaction and diffusion are integrated in a similar nonstiff fashion the comparison in based on both one dimensional 1d unsteady tests of a premixed methane air flame and unsteady two dimensional tests of the same flame interacting with a counterrotating vortex pair in both cases the grimech1 2 reaction mechanism with 32 species and 177 elementary reactions is used computed results show that the stiff reaction scheme enables selection of larger time steps and thus leads to substantial improvement in the performance of the computations for the present reaction mechanism and flame conditions speedup factors of about 10 are achieved in the 1d tests and about five in two dimensions possible extensions of the present scheme to further improve efficiency are also discussed

a splitting method for incompressible flows with variable density based on a pressure poisson equation a new fractional time stepping technique for solving incompressible flows with variable density is proposed the main feature of this method is that as opposed to other known algorithms the pressure is determined by just solving one poisson equation per time step which greatly reduces the computational cost the stability of the method is proved and the performance of the method is numerically illustrated

enlarging the domain of attraction and maximising convergence rate for delta operator systems with actuator saturation in this paper we present some new approaches to improve the feedback property of delta operator systems with actuator saturation both enlarging the domain of attraction and maximising the convergence rate are the desired feedback properties the lifting technique is used to enlarge the domain of attraction for the delta operator systems subject to actuator saturation the maximisation of convergence rate is realised by designing control gain inside a given ellipsoid a necessary and sufficient condition is proposed for the contractive invariance of the given ellipsoid simulation results are provided to demonstrate the effectiveness of the developed techniques

dynamical system design from a control perspective finite frequency positive realness approach nonminimum phase zeros are well known to limit the best achievable control performance when the control gain is allowed to be arbitrarily high on the other hand the phase crossover appears to be a limiting factor for performance when high gain controllers are not allowed in particular the positive realness in a finite frequency range seems crucial for achieving good performance in the presence of control constraints this paper will first give multiple reasons to support this conjecture and then develop a systematic method for designing mechanical systems to achieve the finite frequency positive real ffpr property specifically we present a state space characterization of the ffpr property by generalizing the well known kalman yakubovich popov lemma to deal with a class of frequency domain inequalities that are required to hold within a finite frequency interval the result is further extended for uncertain systems to give a sufficient condition for satisfaction of a robust ffpr property the nominal ffpr result is interpreted in the time domain in terms of input output signals finally we show that certain sensor actuator placement problems to achieve the ffpr property can be reduced to finite dimensional convex problems involving linear matrix inequalities the method is applied to the shape design of a swing arm for magnetic storage devices with the objective of maximizing the control bandwidth achievable with a limited actuator power

a high order numerical algorithm for dns of low mach number reactive flows with detailed chemistry and quasi spectral accuracy a novel and efficient algorithm is presented in this paper to deal with dns of turbulent reacting flows under the low mach number assumption with detailed chemistry and a quasi spectral accuracy the temporal integration of the equations relies on an operating split strategy where chemical reactions are solved implicitly with a stiff solver and the convection diffusion operators are solved with a runge kutta chebyshev method the spatial discretisation is performed with high order compact schemes and a fft based constant coefficient spectral solver is employed to solve a variable coefficient poisson equation the numerical implementation takes advantage of the 2decomp fft libraries developed by 1 which are based on a pencil decomposition method of the domain and are proven to be computationally very efficient an enhanced pressure correction method is proposed to speed up the achievement of machine precision accuracy it is demonstrated that a second order accuracy is reached in time while the spatial accuracy ranges from fourth order to sixth order depending on the set of imposed boundary conditions the software developed to implement the present algorithm is called holomac and its numerical efficiency opens the way to deal with dns of reacting flows to understand complex turbulent and chemical phenomena in flames

an hybrid finite volume finite element method for variable density incompressible flows this paper is devoted to the numerical simulation of variable density incompressible flows modeled by the navier stokes system we introduce an hybrid scheme which combines a finite volume approach for treating the mass conservation equation and a finite element method to deal with the momentum equation and the divergence free constraint the breakthrough relies on the definition of a suitable footbridge between the two methods through the design of compatibility condition in turn the method is very flexible and allows to deal with unstructured meshes several numerical tests are performed to show the scheme capabilities in particular the viscous rayleigh taylor instability evolution is carefully investigated

delta operator formulated discrete time approximations of continuous time systems given a continuous time system a technique to directly obtain an approximate delta operator formulated discrete time system spl delta system is presented for this purpose the analog of the well known boxer thaler integrators q forms applicable to shift operator formulated discrete time systems q systems are derived for spl delta systems next using these spl delta forms a method to obtain an approximate spl delta system of a given continuous time system is derived this algorithm is easily implementable in a computer with little computational burden it is shown that as sampling time decreases the spl delta system thus obtained yields the given continuous time system further verifying the close equivalence between this formulation and continuous time systems two examples illustrating advantages that may be gained by utilizing these spl delta forms in digitizing analog systems are also included

regular article a semi implicit numerical scheme for reacting flow a stiff 1 operator split projection scheme is constructed for the simulation of unsteady two dimensional reacting flow with detailed kinetics the scheme is based on the compressible conservation equations for an ideal gas mixture in the zero mach number limit the equations of motion are spatially discretized using second order centered differences and are advanced in time using a new stiff predictor corrector approach the new scheme is a modified version of the additive stiff scheme introduced in a previous effort by h n najm p s wyckoff and o m knio 1998 j comput phys 143 381 the predictor updates the scalar fields using a strang type operator split integration step which combines several explicit diffusion sub steps with a single stiff step for the reaction terms such that the global time step may significantly exceed the critical diffusion stability limit convection terms are explicitly handled using a second order multi step scheme the velocity field is advanced using a projection scheme which consists of a partial convection diffusion update followed by a pressure correction step a split approach is also adopted for the convection diffusion step in the momentum update this splitting combines an explicit treatment of the convective terms at the global time step with several explicit fractional steps for diffusion finally a corrector step is implemented in order to couple the evolution of the density and velocity fields and to stabilize the computations the corrector acts only on the convective terms and the pressure field while the predicted updates due to diffusion and reaction are left unchanged the correction of the scalar fields is implemented using a single step non split non stiff second order time integration a similar procedure is used for the velocity field which is followed by a pressure projection step the performance and behavior of the operator split scheme are first analyzed based on tests for a nonlinear reaction diffusion equation in one space dimension followed by computations with a detailed c 1c 2 methane air mechanism in one and two dimensions the tests are used to verify that the scheme is effectively second order in time and to suggest guidelines for selecting integration parameters including the number of fractional diffusion steps and tolerance levels in the stiff integration for two dimensional simulations with the present reaction mechanism flame conditions and resolution parameters speedup factors of about 5 are achieved over the previous additive scheme and about 25 over the original explicit scheme

output feedback control for t s fuzzy delta operator systems with time varying delays via an input output approach in this paper the problem of induced l 2 disturbance attenuation control is investigated for discrete time takagi sugeno t s fuzzy delta operator systems with time varying delays and disturbance input via an input output approach a model transformation method is utilized to approximate the time varying delay in t s fuzzy delta operator systems by applying the scaled small gain ssg theorem and lyapunov krasovskii functional approach a sufficient condition is established to guarantee that the closed loop system is asymptotically stable and has an induced l 2 disturbance attenuation performance the existence condition of the dynamic output feedback controller can be solved via convex optimization problems finally simulation results are given to demonstrate the feasibility and effectiveness of the proposed method

suitable initial conditions the aim of this note is to draw attention on a computational problem related to the initial simulation of very large time dependant systems the underlying theoretical problem has been the object of many relevant works in theoretical mathematical and in the applied literature as well but it seems that its computational implications are not familiar to many we hope that this note will help

time accurate calculation of variable density flows with strong temperature gradients and combustion a time accurate algorithm is proposed for low mach number variable density flows with or without chemical reactions the algorithm is based on a predictor corrector time integration scheme that employs a projection method for the momentum equation a constant coefficient poisson equation is solved for the pressure following both the predictor and corrector steps to fully satisfy the continuity equation at each time step spatial discretization is performed on a collocated grid system that offers computational simplicity and straightforward extension to curvilinear coordinate systems to avoid the pressure odd even decoupling that is typically encountered in such grids a flux interpolation technique is introduced for the equations governing variable density flows an important characteristic of the proposed algorithm is that it can be applied to flows in both open and closed domains its robustness and accuracy are illustrated with a series of numerical experiments in particular we present simulations of non isothermal turbulent channel flow as well as simulations of a premixed flame vortex interaction

robust stability of discrete time systems using delta operators the robust stability of discrete time systems formulated in terms of the delta delta operator is discussed that is given the nominal characteristic equation p delta of a discrete time system it is of interest to know how much the coefficients can be perturbed while preserving stability a procedure to obtain the maximum intervals for a perturbed polynomial p delta to still be stable is presented

a projection fem for variable density incompressible flows abstract this work describes a new finite element projection method for the computation of incompressible viscous flows of nonuniform density one original idea of the proposed method consists in factorizing the density variable partly outside and partly inside the time evolution operator in the momentum equation to prevent spatial discretization errors in the mass conservation to affect the kinetic energy balance of the fluid it is shown that unconditional stability in the incremental version of the projection method is possible provided two projections are performed per time step in particular a second order accurate bdf projection method is presented and its numerical performance is illustrated by test computations and comparisons

multi constrained fault estimation observer design with finite frequency specifications for continuous time systems the design of a multi constrained full order fault estimation observer ffeo with finite frequency specifications is studied for continuous time systems by constructing an augmented system a multi constrained ffeo in finite frequency domain is proposed to achieve fault estimation meanwhile the presented ffeo can avoid the overdesign problem generated by the entire frequency domain by the generalised kalman yakubovich popov lemma furthermore by introducing slack variables improved results on ffeo design in different frequency domains are obtained such that different lyapunov matrices can be separately designed for each constraint simulation results are presented to demonstrate the effectiveness and potentials of the proposed techniques

boundary conditions for a divergence free velocity pressure formulation of the navier stokes equations new sets of boundary conditions for the velocity pressure formulation of the incompressible navier stokes equations are derived the boundary conditions have the same form on both inflow and outflow boundaries and lead to a divergence free solution moreover the specific form of the boundary condition makes it possible derive a symmetric positive definite equation system for the internal pressure numerical experiments support the theoretical conclusions

non fragile positive real control for uncertain linear neutral delay systems this paper is concerned with the problem of non fragile positive real control for uncertain neutral delay systems with time invariant norm bounded parameter uncertainty time delays are assumed to appear in both the state and the controlled output equations the state feedback gains are with norm bounded controller uncertainties for both the cases with additive and multiplicative controller uncertainties we address the problem of designing memoryless state feedback controllers such that for all admissible uncertainties the resulting closed loop system is stable and the closed loop transfer function is extended strictly positive real sufficient conditions for the existence of desired controllers are given in terms of linear matrix inequalities lmis when these lmis are feasible the expected memoryless state feedback controller can be easily constructed via convex optimization an illustrative example is given to demonstrate the validity and applicability of the proposed approach

analysis and design for delta operator systems with actuator saturation this paper studies the problem of estimating the domain of attraction for delta operator systems under a pre designed saturated linear feedback ellipsoids are used as the estimates of domain of attraction for delta operator systems a condition is derived in terms of an auxiliary feedback matrix for determining if a given ellipsoid is contractive invariant moreover some necessary and sufficient conditions are given for single input delta operator systems the results are also further extended to determine the invariant sets for delta operator systems with persistent disturbances linear matrix inequalities based on the methods are developed for constructing feedback laws that achieve disturbance rejection with guaranteed stability requirements some numerical examples are given to illustrate the effectiveness of the developed techniques

conservative high order finite difference schemes for low mach number flows three finite difference algorithms are proposed to solve a low mach number approximation for the navier stokes equations these algorithms exhibit fourth order spatial and second order temporal accuracy they are dissipation free and thus well suited for dns and les of turbulent flows the key ingredient common to each of the methods presented is a poisson equation with variable coefficient that is solved for the hydrodynamic pressure this feature ensures that the velocity field is constrained correctly it is shown that this approach is needed to avoid violation of the conservation of kinetic energy in the inviscid limit which would otherwise arise through the pressure term in the momentum equation an existing set of finite difference formulae for incompressible flow is generalized to handle arbitrary large density fluctuations with no violation of conservation through the non linear convective terms an algorithm which conserves mass momentum and kinetic energy fully is obtained when an approximate equation of state is used instead of the exact one results from a model problem are used to show both spatial and temporal convergence rates and several test cases are presented to illustrate the performance of the algorithms

gauge uzawa methods for incompressible flows with variable density two new gauge uzawa schemes are constructed for incompressible flows with variable density one is in the conserved form while the other is in the convective form it is shown that the first order versions of both schemes in their semi discretized form are unconditionally stable numerical experiments indicate that the first order resp second order versions of the two schemes lead to first order resp second order convergence rate for all variables and that these schemes are suitable for handling problems with large density ratios such as in the situation of air bubble rising in water

a robust colocated implicit algorithm for direct numerical simulation of compressible turbulent flows a non dissipative robust implicit algorithm is proposed for direct numerical and large eddy simulation of compressible turbulent flows the algorithm addresses the problems caused by low mach numbers and under resolved high reynolds numbers it colocates variables in space to allow easy extension to unstructured grids and discretely conserves mass momentum and total energy the navier stokes equations are non dimensionalized using an incompressible scaling for pressure and the energy equation is used to obtain an expression for the velocity divergence a pressure correction approach is used to solve the resulting equations such that the discrete divergence is constrained by the energy equation as a result the discrete equations analytically reduce to the incompressible equations at very low mach number i e the algorithm overcomes the acoustic time scale limit without preconditioning or solution of an implicit system of equations the algorithm discretely conserves kinetic energy in the incompressible inviscid limit and is robust for inviscid compressible turbulence on the convective time scale these properties make it well suited for dns les of compressible turbulent flows results are shown for acoustic propagation the incompressible taylor problem periodic shock tube problem and isotropic turbulence

improved finite word length characteristics in digital control using delta operators this paper examines some of the consequences of finite word lengths in digital control it is shown that in many cases of practical importance the usual shift operator formulation is inferior to an alternative formulation which we designate the delta operator approach this latter approach is shown to give better coefficient representation and less roundoff noise in many cases we thus argue that the shift operator and its associated z transform can be replaced by delta operators and their associated transform which we designate a transform an added advantage of this approach is that discrete designs and transforms converge to their continuous time counterparts as the sampling rate is increased

high order splitting methods for the incompressible navier stokes equations abstract a new pressure formulation for splitting methods is developed that results in high order time accurate schemes for the solution of the incompressible navier stokes equations in particular improved pressure boundary conditions of high order in time are introduced that minimize the effect of erroneous numerical boundary layers induced by splitting methods a new family of stiffly stable schemes is employed in mixed explicit implicit time intgration rules these schemes exhibit much broader stability regions as compared to adams family schemes typically used in splitting methods their stability properties remain almost constant as the accuracy of the integration increases so that robust third or higher order time accurate schemes can readily be constructed that remain stable at relatively large cfl number the new schemes are implemented within the framework of spectral element discretizations in space so that flexibility and accuracy is guaranteed in the numerical experimentation a model stokes problem is studied in detail and several examples of navier stokes solutions of flows in complex geometries are reported comparison is made with the previously used first order in time spectral element splitting and non splitting e g uzawa schemes high order splitting spectral element methods combine accuracy in space and time and flexibility in geometry and thus can be very efficient in direct simulations of turbulent flows in complex geometries

consensus protocol for multiple delta operator systems abstract this paper considers the consensus problem for multiple delta operator systems the delta operator system is used for describing continuous time processes at rapid rate sampling moreover it has been known to have a better model description than that of shift operator based approach in this work the consensus problem is first solved for multiple delta operator systems by using a distributed protocol to achieve consensus several important sufficient conditions are proposed for the fast sampling systems based on the parameters of sampling period and consensus strength in addition two distributed consensus protocols are constructed for both leaderless multiple delta operator systems and leader follower multiple delta operator systems finally the effectiveness of the theoretical analysis is illustrated by applying distributed consensus protocol to multi agent systems

generalized conservative approximations of split convective derivative operators we propose a strategy to design locally conservative finite difference approximations of convective derivatives for shock free compressible flows with arbitrary order of accuracy that generalizes the approach of ducros et al 2000 1 and that can be applied as a building block of low dissipative hybrid shock capturing methods the approximations stem from application of standard central difference formulas to split forms of the convective terms in the compressible euler equations which guarantee strong numerical stability and near energy preservation in the inviscid limit a convenient implementation of the high order fluxes is suggested which guarantees improved computational efficiency over existing methods numerical tests performed for isotropic turbulence at zero viscosity show stability of schemes with order of accuracy up to ten and effectiveness of convective splitting of kennedy and gruber 2008 2 in providing extra stability in the presence of strong density variations numerical simulations of compressible turbulent boundary layer flow indicate suitability of the method for non uniform grids and overall support superior computational efficiency of high order schemes

a second order projection method for variable density flows this paper describes a second order projection method for variabledensity incompressible flows the method is suitable for both finite amplitude density variations and for fluids that are modeled using a boussinesq approximation it is based on a second order fractional step scheme in which diffusion convection terms are advanced without enforcing the incompressibility condition and the resulting intermediate velocity field is then projected onto the space of discretely divergencefree vector fields the nonlinear convection terms are treated using a godunov type procedure that is second order for smooth flow and remain stable and non oscillatory for nonsmooth flows with low fluid viscosities the method is described for finite amplitude density variation and the simplifications for a boussinesq approximation are sketched numerical results are presented that validate the convergence properties of the method and demonstrate its performance on more realistic problems

time domain interpretations of frequency domain inequalities on semi finite ranges many of the significant results in systems and control literature rely on characterizations of system properties in terms of frequency domain inequalities fdis and or time domain inequalities tdis classical fdis required to hold on the entire frequency range have been interpreted by equivalent tdis so that satisfaction of one implies that of the other recent developments have addressed fdis within semi finite frequency ranges to increase flexibility in the system analysis and synthesis this paper provides necessary and sufficient conditions for a general fdi to hold within a restricted frequency range in terms of tdis to be satisfied by every inputs within a certain class specified by a matrix valued integral quadratic constraint

a numerical method for dns les of turbulent reacting flows a spatially non dissipative implicit numerical method to simulate turbulent reacting flows over a range of mach numbers is described the compressible navier stokes equations are rescaled so that the zero mach number equations are discretely recovered in the limit of zero mach number the dependent variables are co located in space and thermodynamic variables are staggered from velocity in time the algorithm discretely conserves kinetic energy in the incompressible inviscid non reacting limit the chemical source terms are implicit in time to allow for stiff chemical mechanisms the algorithm is readily extended to complex chemical mechanisms numerical examples using both simple and complex chemical mechanisms are presented

on mass conservation and solvability of the discretized variable density zero mach navier stokes equations abstract this paper presents a study of a variable density zero mach navier stokes discretization and subtle difficulties associated with the pressure a consistent compatibility condition between the temperature and velocity is first derived it is shown that one can accurately and efficiently solve all primitive variables with the uniform order of accuracy of the specified discretization and time marching schemes it is shown that higher density ratios than those commonly reported are achievable and it is in part traced to the solvability of the pressure poisson equation numerical properties are assessed with the method of manufactured solutions and several classical testing cases are carried out to validate the code ranging from variable density mixing problems to premixed flames with forcing

positive real control problem for uncertain linear time invariant systems abstract this paper focuses on positive real control of linear time invariant systems which are subjected to norm bounded uncertainties in the state equation we address the problem of designing a linear dynamic output feedback controller that robustly stabilizes the uncertain system and achieves the extended strict positive realness property for a given closed loop transfer function it is shown that a solution to the above problem can be obtained by solving a scaled strict positive real control problem for which no parameter uncertainty occurs

high order compact alternating direction implicit method for the generalized sine gordon equation high order compact alternating direction implicit scheme is given for solving the generalized sine gordon equation in a two dimensional rectangular domain we apply the compact finite difference operators to obtain a fourth order discretization for the second order space derivatives and we give a linearized three time level algorithm for solving the original nonlinear equation error estimate is given by the energy method numerical results are provided to verify the accuracy and efficiency of this algorithm

null controllable region of delta operator systems subject to actuator saturation abstractin this paper we give exact description of null controllable regions for delta operator systems subject to actuator saturation the null controllable region is in terms of a set of extremal trajectories of anti stable subsystems for the delta operator system with real eigenvalues or complex eigenvalues the description is simplified to an explicit formula which is used to characterise the boundary of a null controllable region the relations of null controllable regions are shown separately for continuous time systems discrete time systems and delta operator systems two numerical examples are given to illustrate the effectiveness of the proposed techniques on null controllable regions

numerical solution of the nonlinear klein gordon equation a numerical method is developed to solve the nonlinear one dimensional klein gordon equation by using the cubic b spline collocation method on the uniform mesh points we solve the problem for both dirichlet and neumann boundary conditions the convergence and stability of the method are proved the method is applied on some test examples and the numerical results have been compared with the exact solutions the l 2 l and root mean square errors rms in the solutions show the efficiency of the method computationally

direct numerical simulation of a turbulent reactive plume on a parallel computer a computational algorithm is described for direct numerical simulation dns of a reactive plume in spatially evolving grid turbulence the algorithm uses sixth order compact differencing in conjunction with a fifth order compact boundary scheme which has been developed and found to be stable a compact filtering method is discussed as a means of stabilizing calculations where the viscous diffusive terms are differenced in their conservative form this approach serves as an alternative to nonconservative differencing previously advocated as a means of damping the 2 waves in numerically solving the low mach number equations the time derivative of the density field in the pressure poisson equation was found to be the most destabilizing part of the calculation even ordered finite difference approximations to this derivative were found to be more stable allow for larger density gradients than odd ordered approximations turbulence at the inlet boundary is generated by scanning through an existing three dimensional field of fully developed turbulence in scanning through the inlet field it was found that a high order interpolation e g cubic spline interpolation is necessary in order to provide continuous velocity derivatives regarding pressure a neumann inlet condition combined with a dirichlet outlet condition was found to work well the chemistry follows the single step irreversible global reaction fuel r oxidizer 1 r product heat with parameters chosen to match experimental data as far as allowed by resolution constraints simulation results are presented for four different cases in order to examine the effects of heat release damkohler number and arrhenius kinetics on the flow physics statistical data from the dns are compared to theory and wind tunnel data and found in reasonable agreement with regard to growth of turbulent length scales decay of turbulent kinetic energy decay of centerline scalar concentration decrease in scalar rms and spread of plume profile reactive scalar statistics are consistent with expected behavior

delay dependent robust stability and stabilization of uncertain linear delay systems a linear matrix inequality approach this paper considers the problems of robust stability analysis and robust control design for a class of uncertain linear systems with a constant time delay the uncertainty is assumed to be norm bounded and appears in all the matrices of the state space model we develop methods for robust stability analysis and robust stabilization the proposed methods are dependent on the size of the delay and are given in terms of linear matrix inequalities

an overview of projection methods for incompressible flows a series of numerical issues related to the analysis and implementation of fractional step methods for incompressible flows are addressed in this paper these methods are often referred to in the literature as projection methods and can be classified into three classes namely the pressure correction methods the velocity correction methods and the consistent splitting methods for each class of schemes theoretical and numerical convergence results available in the literature are reviewed and open questions are discussed the essential results are summarized in a table which could serve as a useful reference to numerical analysts and practitioners published by elsevier b v msc 65n35 65n22 65f05 35j05

a solution for the odd even decoupling problem in pressure correction algorithms for variable density flows when the navier stokes equations are solved on a colocated mesh a spurious mode for the pressure can appear if no special attention is paid to the discretization of the pressure this pressure mode can be suppressed by a pressure weighted interpolation formula for the mass flux over a cell face in this paper a similar cure is presented in the framework of pressure correction methods in variable density flow special attention is given to the solvability condition for the resulting poisson like equation for the pressure it consists of two remedies a correction term for the cell face velocity is introduced and the stencil for the discrete laplacian is compacted we finally show the applicability of the method on general curvilinear coordinate systems in three dimensions

skew symmetric form of convective terms and fully conservative finite difference schemes for variable density low mach number flows the form of convective terms for compressible flow equations is discussed in the same way as for an incompressible one by morinishi et al y morinishi t s lund o v vasilyev p moin fully conservative higher order finite difference schemes for incompressible flow j comput phys 124 1998 90 and fully conservative finite difference schemes suitable for shock free unsteady compressible flow simulations are proposed commutable divergence advective and skew symmetric forms of convective terms are defined by including the temporal derivative term for compressible flow these forms are analytically equivalent if the continuity is satisfied and the skew symmetric form is secondary conservative without the aid of the continuity while the divergence form is primary conservative the relations between the present and existing quasi skew symmetric forms are also revealed commutable fully discrete finite difference schemes of convection are then derived in a staggered grid system and they are fully conservative provided that the corresponding discrete continuity is satisfied in addition a semi discrete convection scheme suitable for compact finite difference is presented based on the skew symmetric form the conservation properties of the present schemes are demonstrated numerically in a three dimensional periodic inviscid flow the proposed fully discrete fully conservative second order accurate scheme is also used to perform the dns of compressible isotropic turbulence and the simulation of open cavity flow

adaptive h control in finite frequency domain for uncertain linear systems this paper is concerned with the problem of adaptive h controller design in finite frequency domains for uncertain linear systems the uncertainties are assumed to be time invariant unknown but bounded which appear affinely in the matrices of system models an adaptive mechanism is introduced to construct a novel finite frequency h controller with time varying gains by using lyapunov theory and parseval s theorem the controller design conditions are given in terms of a set of linear matrix inequalities lmis it is shown that the proposed finite frequency controller with time varying gains can achieve better h performance than the traditional ones with fixed gains finally a numerical example of the f 18 aircraft model is given to illustrate the presented theoretical results

h sub spl infin and positive real control for linear neutral delay systems this note is concerned with the h sub spl infin and positive real control problems for linear neutral delay systems the purpose of h sub spl infin control is the design of a memoryless state feedback controller which stabilizes the neutral delay system and reduces the h sub spl infin norm of the closed loop transfer function from the disturbance to the controlled output to a prescribed level while the purpose of positive real control is to design a memoryless state feedback controller such that the resulting closed loop system is stable and the closed loop transfer function is extended strictly positive real sufficient conditions for the existence of the desired controllers are given in terms of a linear matrix inequality lmi when this lmi is feasible the expected memoryless state feedback controllers can be easily constructed via convex optimization

interactions among periodic waves and solitary waves of the n 1 dimensional sine gordon field exact solutions of the n 1 dimensional sine gordon field equation are studied with help of those of the cubic nonlinear klein gordon fields the mapping relations among the sine gordon field equation and the cubic nonlinear klein gordon fields are pure algebraic by solving the cubic nonlinear klein gordon equations many new types of exact explicit solutions such as the periodic periodic interaction waves periodic kink interaction waves periodic perturbed snake shape solitary waves etc are displayed both analytically and graphically

numerically stable formulations of convective terms for turbulent compressible flows abstract a systematic analysis of the discrete conservation properties of non dissipative central difference approximations of the compressible navier stokes equations is reported a generalized splitting of the nonlinear convective terms is considered and energy preserving formulations are fully characterized by deriving a two parameter family of split forms previously developed formulations reported in literature are shown to be particular members of this family novel splittings are introduced and discussed as well furthermore the conservation properties yielded by different choices for the energy equation i e total and internal energy entropy are analyzed thoroughly it is shown that additional preserved quantities can be obtained through a suitable adaptive selection of the split form within the derived family local conservation of primary invariants which is a fundamental property to build high fidelity shock capturing methods is also discussed in the paper numerical tests performed for the taylor green vortex at zero viscosity fully confirm the theoretical findings and show that a careful choice of both the splitting and the energy formulation can provide remarkably robust and accurate results

low frequency positive real control for delta operator systems a strictly positive real control problem for delta operator systems in a low frequency range is presented by using the generalized kalman yakubovic popov lemma the objective of the strictly positive real control problem is to design a controller such that the transfer function is strictly positive real and the resulting closed loop system is stable sufficient conditions for the low frequency strictly positive real controller of the closed loop delta operator systems are presented in terms of solutions to a set of linear matrix inequalities a numerical example is given to illustrate the effectiveness and potential for the developed techniques

a conservative adaptive projection method for the variable density incompressible navier stokes equations in this paper we present a method for solving the equations governing time dependent variable density incompressible flow in two or three dimensions on an adaptive hierarchy of grids the method is based on a projection formulation in which we first solve advection diffusion equations to predict intermediate velocities and then project these velocities onto a space of approximately divergence free vector fields our treatment of the first step uses a specialized second order upwind method for differencing the nonlinear convection terms that provides a robust treatment of these terms suitable for inviscid and high reynolds number flow density and other scalars are advected in such a way as to maintain conservation if appropriate and free stream preservation our approach to adaptive refinement uses a nested hierarchy of logically rectangular girds with simultaneous refinement of the girds in both space and time the integration algorithm on the grid hierarchy is a recursive procedure in which coarse grids are advanced in time fine grids are advanced multiple steps to reach the same time as the coarse grids and the data at different levels are then synchronized the single grid algorithm is described briefly but the emphasis here is on the time stepping procedure for the adaptive hierarchy numerical examples are presented to demonstrate the algorithms s accuracy and convergence properties and illustrate the behavior of the method an additional example demonstrates the performance of the method on a more realistic problem namely a three dimensional variable density shear layer

generalized kyp lemma unified frequency domain inequalities with design applications the celebrated kalman yakubovic spl caron popov kyp lemma establishes the equivalence between a frequency domain inequality fdi and a linear matrix inequality and has played one of the most fundamental roles in systems and control theory this paper first develops a necessary and sufficient condition for an s procedure to be lossless and uses the result to generalize the kyp lemma in two aspects the frequency range and the class of systems and to unify various existing versions by a single theorem in particular our result covers fdis in finite frequency intervals for both continuous discrete time settings as opposed to the standard infinite frequency range the class of systems for which fdis are considered is no longer constrained to be proper and nonproper transfer functions including polynomials can also be treated we study implications of this generalization and develop a proper interface between the basic result and various engineering applications specifically it is shown that our result allows us to solve a certain class of system design problems with multiple specifications on the gain phase properties in several frequency ranges the method is illustrated by numerical design examples of digital filters and proportional integral derivative controllers