Space–Time Chaos,Critical Phenomena,and Bifurcations of Solutions of Infinite-Dimensional Systems of Ordinary Differential Equations

We study infinite-dimensional systems of ordinary differential equations having applications in some popular and important physical problems. The appearance of infinite-dimensional space–time chaos is considered, namely, the bifurcations and critical phenomena that occur in the phase space of the systems and explain some physical problems are described.     

Optimal control of unsteady compressible viscous flows

The control of complex, unsteady flows is a pacing technology for advances in fluid mechanics. Recently, optimal control theory has become popular as a means of predicting best case controls that can guide the design of practical flow control systems. However, most of the prior work in this area has focused on incompressible flow which precludes many of the important physical flow phenomena that must be controlled in practice including the coupling of fluid dynamics, acoustics, and heat transfer. This paper presents the formulation and numerical solution of a class of optimal boundary control problems governed by the unsteady two‐dimensional compressible Navier–Stokes equations. Fundamental issues including the choice of the control space and the associated regularization term in the objective function, as well as issues in the gradient computation via the adjoint equation method are discussed. Numerical results are presented for a model problem consisting of two counter‐rotating viscous vortices above an infinite wall which, due to the self‐induced velocity field, propagate downward and interact with the wall. The wall boundary control is the temporal and spatial distribution of wall‐normal velocity. Optimal controls for objective functions that target kinetic energy, heat transfer, and wall shear stress are presented along with the influence of control regularization for each case. Copyright © 2002 John Wiley & Sons, Ltd.     

An algebraic multigrid solver for Navier-Stokes problems

An efficient numerical method is presented for solving the equations of motion for viscous fluids. The equations are discretized on the basis of unstructured finite element meshes and then solved by direct iteration. Advective fluxes are temporarily fixed at each iteration to provide a linearized set of coupled equations which are then also solved by iteration using a fully implicit algebraic multigrid (AMG) scheme. A rapid convergence to machine accuracy is achieved that is almost mesh-independent. The scaling of computing time with mesh size is therefore close to the optimum.     

On superlinear convergence of quasi-Newton methods for nonsmooth equations

We show that strong differentiability at solutions is not necessary for superlinear convergence of quasi-Newton methods for solving nonsmooth equations. We improve the superlinear convergence result of Ip and Kyparisis for general quasi-Newton methods as well as the Broyden method. For a special example, the Newton method is divergent but the Broyden method is superlinearly convergent.     

An algorithm for the computation of consistent initial values for differential–algebraic equations

We use boundary value methods to compute consistent initial values for fully implicit nonlinear differential-algebraic equations. The obtained algorithm uses variable order formulae and a deferred correction technique to evaluate the error. A rigorous theory is stated for nonlinear index 1, 2 and 3 DAEs of Hessenberg form. Numerical tests on classical index 1, 2 and 3 DAE problems are reported. This revised version was published online in June 2006 with corrections to the Cover Date.     

Energy inconsistency of Galerkin approximationsfor plane Couette flow

For classical solutions of the incompressible Navier-Stokes equations (NSE) the energybalance between kinetic energy, work done by external forces, and viscous dissipation holds rigorously true. It is shown in this paper that standard Galerkin approximations violate energy balance in the case of plane Couette flow, whereas Poiseuille flow turns out to be energy consistent at any cutoff. The main reason for this discrepancy is seen in the different boundary conditions between the stationary linear shear flow and its disturbances. In our analysis, essentially, we introduce an auxiliary external force field which enforces the finite dimensional Galerkin approximation to fulfil the NSE. It is exemplarily demonstrated how the energy discrepancy decreases when the number of disturbed modes is increased which couple to the stationary shear flow.     

Decay of confined,two-dimensional,spatially periodic arrays of vortices: A numerical investigation

The disarrangement of a perturbed lattice of vortices was studied numerically. The basic state is an exponentially decaying, exact solution of the Navier-Stokes equations. Square arrays of vortices with even numbers of vortex cells along each side were perturbed and their evolution was investigated. Whether the energy in the perturbation grows somewhat before it decays or decays monotonically depends on the initial strength of the vortices of the basic state, the extent of lateral confinement and the structure of the perturbation. The critical condition for temporally local instability, i.e. the critical amplitude of the basic state that must be exceeded to allow energy transfer from the basic state to the perturbation, is discussed. In the strongly confined case of a square lattice of four vortices the appearance of enchancement of global rotation is the result of energy transfer from the basic state to a temporally local unstable mode. Energy is transferred from the basic state to larger-scaled structures (inverse cascade) only if the scales of the larger structures are inherently contained in the initial structure of the perturbation. The initial structure of the double array of vortices is not maintained except for a very special form of perturbation. The facts that large scales decay more slowly than small scales and that, when non-linearities are sufficiently strong, energy is transferred from one scale to another explain the differences in the disarrangement process for different initial strengths of the vortices of the basic state. The stronger vortices, i.e. the vortices perturbed in a manner that increases their strength, tend to dominate the weaker vortices. The pairing and subsequent merging (or capture) of vortices of like sense into larger-scale vortices are described in terms of peaks in the evolution of the square root of the palinstrophy divided by the enstrophy.     

A NUMERICAL EMBEDDING METHOD FOR SOLVING THE NONLINEAR COMPLEMENTARITY PROBLEM(Ⅱ)-ALGORITHM AND ITS CONVERGENCE

In this paper,based on the results presented in part Ⅰ of this paper[18],we present a numerical embedding algorithm for solving the nonlinear complementarityproblem, and prove its convergence carefully. Numerical experiments show that thealgorithm is successful.     

Optimal inspection and perfect repair

Models are developed for decision making about monitoring andmaintenance of systems whose performance through time is describedby a general stochastic process. The system is monitored andpreventive and corrective maintenance actions are carried outin response to the observed system state. The decision processis simplified by using the maximum process as a decision variable.The models developed generalize age replacement models and othersimple maintenance strategies. The approach can deal with failuresthat prevent the system functioning further, and also failuresthat are defined by regulation or economic considerations. Attentionis restricted to perfect repair and inspection, but the structureprovides the framework for further developments.     

OSCILLATION FOR SOLUTIONS OF SYSTEMS OF N-TH ORDER PARTIAL FUNCTIONAL DIFFERENTIAL EQUATIONS

In this paper, some sufficient conditions for the oscillation for solutions to systems of n-th order partial functional differential equations are obtained.