Phase-winding Solutions for Axisymmetric Convection between Rotating Planes Uniformly Heated from Below

The concept of phase winding, where the argument of a complexamplitude function changes slowly with position, may providea theoretical explanation of observed changes in the wavelengthof concentric roll patterns in shallow circular cylinders uniformlyheated from below. Here it is shown that in an infinite fluidlayer, constraints on phase winding associated with the motionnear the centre of the roll pattern are weaker than, for instance,those associated with a two-dimensional pattern confined bya lateral wall. The theoretical treatment, which is based onthe methods of matched asymptotic expansions and multiple scales,includes the effect of rotation of the layer. The solution forthe weakly nonlinear axisymmetric state at slightly supercriticalRayleigh numbers is found by matching an outer solution, expressedin terms of an amplitude function, to an inner solution validnear the axis of rotation. Some observations are also made,concerning modifications to the solution caused by the presenceof an outer cylindrical wall     

Analytical model for deformable roll coating with nip feed

The presented work deals with an analytical solution for a roll coater with deformable rolls, as commenly used in industry. The focus is on the calculation of the nip feed system in a forward coating mode. The calculation is done with thin film theory and includes Hook's law for elasticity for modelling the elasticity behaviour of the deformable roll. The new model has been validated by experimental data from literature. The agreement is sufficiently good. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A theory for synchronization of dynamical systems

In this paper, a theory for synchronization of multiple dynamical systems under specific constraints is developed from a theory of discontinuous dynamical systems. The concepts on synchronization of two or more dynamical systems to specific constraints are presented. The synchronization, desynchronization and penetration of multiple dynamical systems to multiple specified constraints are discussed, and the necessary and sufficient conditions for such synchronicity are developed. The synchronicity of two dynamical systems to a single specific constraint and to multiple specific constraints is investigated. Finally, the synchronization and the corresponding complexity for multiple slave systems with multiple master systems are discussed briefly. The meaning of synchronization for dynamical systems with constraints is extended as a generalized, universal concept. The theory presented in this paper may be as a universal theory for dynamical systems. The paper provides a theoretic frame work in order to control slave systems which can be synchronized with master systems through specific constraints in a general sense.     

Multiple solutions of a boundary-value problem in enzyme kinetics

In some immobilized enzyme systems the steady state of substrate concentration may suddenly change from a low profile to a high profile or vice versa when the physical parameters of the systems pass through certain critical values. This phenomenon is due to the transition from a unique solution to multiple solutions (or vice versa) of the enzyme reaction equation. This problem is studied by considering two physical parameters which represent the internal reaction mechanism and the external influence on the boundary of the reaction-diffusion medium. Both analytical and numerical results for the problem are presented. The analytical results include some sufficient conditions for the existence of multiple steady-state solutions as well as a unique solution. Various numerical results of the problem including time-dependent solutions and their convergence to steady-state solutions are given.     

Modelling and experimental validation of the deflection of a leveller for hot heavy plates

In order to successfully automate levelling processes, in particular for heavy plates, the deflection of the leveller has to be compensated based on a deflection model. In this work, a detailed mathematical deflection model of a hot leveller with bending mechanism and its experimental validation are presented. The roll intermesh profiles are calculated based on the deflection of the work rolls that are elastically supported by support rolls, frames, posts and adjustment screws. The deflection model is suited to compensate the effect of deflection on the roll intermesh and the plate flatness as well as to assess the loads of critical parts, for example the support rolls. A new experimental design to measure the deflection of a leveller is presented and successfully applied for model validation. The work roll deflection is measured directly by means of displacement sensors that are inserted in cut-outs of test plates. These test plates are modelled as linear elastic stripes. For normal load levels, the relative accuracy (repeatability) of the roll intermesh prediction of the model is better than 0.08 mm.     

Computation of the normal forms for general M-DOF systems using multiple time scales. Part I: autonomous systems

This paper is concerned with the symbolic computation of the normal forms of general multiple-degree-of-freedom oscillating systems. A perturbation technique based on the method of multiple time scales, without the application of center manifold theory, is generalized to develop efficient algorithms for systematically computing normal forms up to any high order. The equivalence between the perturbation technique and Poincaré normal form theory is proved, and general solution forms are established for solving ordered perturbation equations. A number of cases are considered, including the non-resonance, general resonance, resonant case containing 1:1 primary resonance, and combination of resonance with non-resonance. “Automatic” Maple programs have been developed which can be executed by a user without knowing computer algebra and Maple. Examples are presented to show the efficiency of the perturbation technique and the convenience of symbolic computation. This paper is focused on autonomous systems, and non-autonomous systems are considered in a companion paper.     

Application of a block modified Chebyshev algorithm to the iterative solution of symmetric linear systems with multiple right hand side vectors

Summary. An adaptive Richardson iteration method is described for the solution of large sparse symmetric positive definite linear systems of equations with multiple right-hand side vectors. This scheme ``learns' about the linear system to be solved by computing inner products of residual matrices during the iterations. These inner products are interpreted as block modified moments. A block version of the modified Chebyshev algorithm is presented which yields a block tridiagonal matrix from the block modified moments and the recursion coefficients of the residual polynomials. The eigenvalues of this block tridiagonal matrix define an interval, which determines the choice of relaxation parameters for Richardson iteration. Only minor modifications are necessary in order to obtain a scheme for the solution of symmetric indefinite linear systems with multiple right-hand side vectors. We outline the changes required. Received April 22, 1993     

Some remarks on period doubling in systems with symmetry

Period doubling of periodic solutions in systems with symmetry leads to certain group theoretical difficulties, if a periodic solution possesses a mixed spatio-temporal symmetry. Based on a result of Vanderbauwhede [11] on period doubling with symmetry a method is presented to determine systematically the bifurcations that one may expect in such a system. The results are used to analyse multiple period doublings of periodic solutions with dihedral group symmetry.     

Localized Roll-Solutions in the Ekman-Couette-System

In the Ekman-Couette-System, where the usual Couette-System is additionally rotated about its normal axis, localized single roll solutions have been known for some time. By varying different system parameters, a new localized solution emerging by saddle-node bifurcations from these known solutions has now been found. This new solution is of the same localized nature like the old solution but with an additional roll. Further bifurcations lead then to an increasing number of rolls, still localized. This behavior is kind of analogous to the so called ‘homoclinic snaking’ which has recently been investigated in conjunction with the Swift-Hohenberg equation and binary fluid convection. It might link the unstable localized single roll solutions with stable multi roll solutions or even with stable periodic roll solutions, which has to be shown yet. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Investigation of the starting conditions of rolling thixotropic polymeric systems

The solution of the problem of flow in the roll gap of a polymeric system whose rheological behavior is described by Leonov's thixotropic model is presented. Analytical expressions are obtained for specific pressures, thrusts, and torques under transient rolling conditions, which permitted developing a method of calculating the maximum values of the energy and force characteristics.Lensovet Leningrad Technological Institute. Translated from Mekhanika Polimerov, No. 1, pp. 133–140, January–February, 1972.     

Nonlinear Spatial Evolution of Small Disturbances from Boundary Conditions in Flat Inclined-Channel Flow

The asymptotic behavior of small disturbances as they evolve spatially from boundary conditions in a flat inclined channel is determined. These small disturbances develop into traveling waves called roll waves, first discussed by Dressler in 1949. Roll waves exist if the Froude number F exceeds 2, which consist of a periodic pattern of bores, or discontinuities. After confirming the instability condition for   F > 2  for the linearized equations in the boundary value case, the nonlinear boundary value problem for the weakly unstable region of F slightly larger than 2 is studied. Multiple scales and the Fredholm alternative theorem are applied to determine the evolution of the solution in space. It is found that the solution is dominated by the evolution of the disturbance along one characteristic. The shock conditions governing the asymptotic solution are determined and these conditions are used to determine the approximate shape of the resulting traveling wave from the solution. Both asymptotic and numerical results for periodic disturbances are presented.     

A numerical evaluation of solvers for the periodic Riccati differential equation

Efficient and accurate structure exploiting numerical methods for solving the periodic Riccati differential equation (PRDE) are addressed. Such methods are essential, for example, to design periodic feedback controllers for periodic control systems. Three recently proposed methods for solving the PRDE are presented and evaluated on challenging periodic linear artificial systems with known solutions and applied to the stabilization of periodic motions of mechanical systems. The first two methods are of the type multiple shooting and rely on computing the stable invariant subspace of an associated Hamiltonian system. The stable subspace is determined using either algorithms for computing an ordered periodic real Schur form of a cyclic matrix sequence, or a recently proposed method which implicitly constructs a stable deflating subspace from an associated lifted pencil. The third method reformulates the PRDE as a convex optimization problem where the stabilizing solution is approximated by its truncated Fourier series. As known, this reformulation leads to a semidefinite programming problem with linear matrix inequality constraints admitting an effective numerical realization. The numerical evaluation of the PRDE methods, with focus on the number of states (n) and the length of the period (T) of the periodic systems considered, includes both quantitative and qualitative results.     

Multiple right-hand side techniques for the numerical simulation of quasistatic electric and magnetic fields

The simulation of slowly varying transient electric high-voltage fields and magnetic fields requires the repeated and successive solution of high-dimensional linear algebraic systems of equations with identical or near-identical system matrices and different right-hand side vectors. For these solution processes which are required within implicit time integration schemes and nonlinear (quasi-)Newton–Raphson methods an iterative multiple right-hand side (mrhs) scheme is used which recycles vector subspaces resulting from previous preconditioned conjugate gradient iteration runs. The combination of this scheme with a subspace projection extrapolation start value generation scheme is discussed. Numerical results for three-dimensional electric and magnetic field simulations are presented and the efficiency of the new schemes re-using eigenvector information from previous iteration processes with different tolerance criteria are compared to those of standard conjugate gradient iterations.     

Deterministic Global Optimization in Nonlinear Optimal Control Problems

The accurate solution of optimal control problems is crucial in many areas of engineering and applied science. For systems which are described by a nonlinear set of differential-algebraic equations, these problems have been shown to often contain multiple local minima. Methods exist which attempt to determine the global solution of these formulations. These algorithms are stochastic in nature and can still get trapped in local minima. There is currently no deterministic method which can solve, to global optimality, the nonlinear optimal control problem. In this paper a deterministic global optimization approach based on a branch and bound framework is introduced to address the nonlinear optimal control problem to global optimality. Only mild conditions on the differentiability of the dynamic system are required. The implementa-tion of the approach is discussed and computational studies are presented for four control problems which exhibit multiple local minima.     

Computing singular solutions to nonlinear analytic systems

Summary A method to generate an accurate approximation to a singular solution of a system of complex analytic equations is presented. Since manyreal systems extend naturally tocomplex analytic systems, this porvides a method for generating approximations to singular solutions to real systems. Examples include systems of polynomials and systems made up of trigonometric, exponential, and polynomial terms. The theorem on which the method is based is proven using results from several complex variables. No special conditions on the derivatives of the system, such as restrictions on the rank of the Jacobian matrix at the solution, are required. The numerical method itself is developed from techniques of homotopy continuation and 1-dimensional quadrature. A specific implementation is given, and the results of numerical experiments in solving five test problems are presented.     

On the approximate solution of some integral equations of the theory of elasticity and mathematical physics PMM vol. 31, no. 6, 1967, pp. 1117–1131

In this paper the development of the method presented in [1] is carried out with application to two types of integral equations encountered in mathematical physics in the investigation of many mixed problems with circular separation line of boundary conditions and in the investigation of plane mixed problems.

The algorithm is given for reducing these integral equations to solution of equivalent infinite linear algebraic systems. It is proved that the resulting infinite systems are quasi completely regular for sufficiently large values of dimensionless parameter A which enters into the systems. It is shown that reduction (truncation) of infinite systems results in finite systems of linear algebraic equations with almost triangular matrices. The last circumstance simplifies considerably the solution of these finite systems after which the solution of initial integral equations is found from simple equations. For given accuracy of the approximate solution and decrease of parameter λ the number of equations in reduced systems increases.

As an example the solution is presented for the axisymmetric problem of a die acting on an elastic layer lying without friction on a rigid foundation.     

Switching dynamics of multiple linear oscillators

In this paper, the switching dynamics of linear oscillators with arbitrary discontinuous forcing are investigated through the concept of switching systems, and such switching systems consist of countable prescribed linear oscillators with different external excitations. The traditional treatments are to smoothen the discontinuity at switching points of two subsystems in a switching system, which can provide an approximate solution only. Therefore, an alternative method is presented to obtain an exact solution of the resultant switching linear system. Under periodic piecewise forcing and random forcing, the corresponding exact solutions and stochastic responses of switching linear systems are developed. For any periodic forcing, the periodic responses and stability of the resultant system composed of multiple linear oscillators in different time intervals are presented. In addition, the resultant switching system consisting of two oscillators are discussed, and the corresponding stability analysis is carried out.     

Direct Shooting Method for the Numerical Solution of Higher-Index DAE Optimal Control Problems

A method for the numerical solution of state-constrained optimal control problems subject to higher-index differential-algebraic equation (DAE) systems is introduced. For a broad and important class of DAE systems (semiexplicit systems with algebraic variables of different index), a direct multiple shooting method is developed. The multiple shooting method is based on the discretization of the optimal control problem and its transformation into a finite-dimensional nonlinear programming problem (NLP). Special attention is turned to the mandatory calculation of consistent initial values at the multiple shooting nodes within the iterative solution process of (NLP). Two different methods are proposed. The projection method guarantees consistency within each iteration, whereas the relaxation method achieves consistency only at an optimal solution. An illustrative example completes this article.     

Planforms in two and three dimensions

Summary When solving systems of PDE with two space dimensions it is often assumed that the solution is spatially doubly periodic. This assumption is usually made in systems such as the Boussinesq equation or reaction-diffusion equations where the equations have Euclidean invariance. In this article we use group theoretic techniques to determine a large class of spatially doubly periodic solutions that are forced to existence near a steady-state bifurcation from a translation-invariant equilibrium.This type of bifurcation problem has been considered by many authors when studying a number of different systems of PDE. Typically, these studies focus at the beginning on equilibria that are spatially periodic with respect to a fixed planar lattice type-such as square or hexagonal. Our focus is different in that we attempt to find all spatially periodic equilibria that bifurcate on all lattices. This point of view leads to some technical simplifications such as being able to restrict to translation free irreducible representations.Of course, many of the types of solutions that we find are well-known-such as hexagon and roll solutions on a hexagonal lattice. This coordinated group theoretic approach does lead, however, to solutions which seem not to have been discussed previously (antisquare solutions on a square lattice) as well as to a more complete classification of the symmetry types of possible solutions. Moreover, our methods extend to triply periodic solutions of PDE with three spatial variables. Some of these results, namely those concerned with primitive cubic lattices, are presented here. The complete results on triply periodic solutions may be found in [6, 7].In honor of Klaus Kirchgässner on the occasion of his sixtieth birthdayResearch supported in part by NSF/DARPA (DMS-8700897) and by the Texas Advanced Research Program (ARP-1100).     

Almost Sure Stability Analysis of Parametric Roll in Random Seas Based on Top Lyapunov Exponent

Stability analysis of the upright position of a ship in random head or following seas is presented. Such seas lead to parametric excitation of roll motion due to periodic variations of the righting lever. The development of simple criteria for the occurrence of parametric induced roll motion in random seas is of major interest for improvement of the international code on intact stability provided by the International Maritime Organization. The stability analysis in random seas is based on the calculation of the top Lyapunov exponent using the fact, that a negative top Lyapunov exponent yields no roll motion. With this findings, roll motion can be excluded for specific sea states. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)