The automatic control of boundary-layer transition — Experiments and computation

In recent years there has been an increasing interest in the control of boundary-layer transition through the use of wall suction. In the current work suction is provided through one or more suction panels situated close to the leading edge of a plate. Experiments show that boundary-layer pressure fluctuation measurements can be used to identify the position of transition. Transition can be maintained at a desired location with minimum power consumption by employing an automatic adaptive feedback control loop which regulates the suction flow rates of two independent suction panels. This can be expressed as a constrained optimization problem. To allow the suction flow rates to be updated, a modified least mean squares algorithm is used within the control loop. Experimental measurements show that the control algorithm allows fast and stable convergence towards the optimum suction distribution for a double suction panel configuration. Numerical simulations have also been performed. The two-dimensional boundary layer was calculated allowing the viscous boundary layer to interact with the inviscid outer flow. Following linear stability theory the spatial growth rates are calculated by solving an Orr-Sommerfeld type eigenvalue problem, with the streamwise location of transition predicted via thee ^N -method. Applying the same optimization strategy as in the experiments, good qualitative agreement between computations and experiments was found. The optimization algorithm has been applied to computer models where the relation between suction flow rates and transition location is described by an empirical analytical function. This shows that the controller can in principle be applied to systems with more than two suction panels.Nomenclature b transition location with zero suction - d desired transition location - e(k) error signal - k iteration index - p rms pressure - p _ref reference rms pressure - r sum of the reference pressure - u streamwise velocity - u _e external velocity - inviscid external velocity - A wave amplitude - F( ) cost function - I identity matrix - N maximum amplification factor - P projection matrix - R Reynolds number - Re Reynolds number based on the boundary-layer thickness - R matrix of weights - Tu turbulence level - vector of suction flow rates - v normal velocity - v _wall suction velocity at the surface - x streamwise coordinates - x _m microphone location - x _T(k) measured transition location - y normal coordinate - y(k) sum of the measured pressures - w(k) noise - plate length - r +i _i - free stream velocity - * displacement thickness - gradient vector - Lagrange multiplier - controller gain - disturbance stream function - disturbance amplitude - wave frequency = complex wave number     

Active vibration isolation of machinery and sensitive equipment using H ∞ control criterion and particle swarm optimization method

Isolating the sensitive equipment from vibrating base or the foundation from machinery vibration is of practical importance in a number of engineering fields. With the development of the vibration control techniques and increasing requirements for the higher-performance vibration isolation in industry and everyday life, active vibration isolation exhibits the best performances. In this paper, active vibration isolation reducing vibration transmitted from vibrating base to sensitive equipment and from machinery to foundation was investigated. Controller as static output feedback was considered to design components of active isolation system. An active control is provided by using H _∞ control criteria to design this controller. This criterion is presented as a cost function and then optimized by Particle Swarm Optimization (PSO) algorithm. The approach is validated using numerical simulation. Results show that this static output feedback H _∞ controller using PSO algorithm can get good performance to reduce the effect of unwanted vibration and disturbances.     

Tracking control of uncertain switched nonlinear cascade systems: a nonlinear H ∞ sliding mode control method

This paper considers the tracking control problem for a class of uncertain switched nonlinear cascade systems via the multiple Lyapunov functions (MLFs) method. Each subsystem under consideration is composed of two cascade-connected parts: the null space dynamics part and the range space dynamics part. The two main robust control strategies, nonlinear H _∞ control (NHC) and the sliding mode control (SMC), are integrated to function in a complementary manner for tracking control tasks. Furthermore, sufficient conditions for the solvability of the tracking control problem of the switched system and design of both switching laws and controllers are presented. Finally, a simulation example is provided to demonstrate the feasibility of the developed method.     

On a Crystalline Variational Problem, Part I:¶First Variation and Global L∞ Regularity

Let φ:ℝ ⁿ → [0,+∞[ be a given positively one-homogeneous convex function, and let ?_φ≔{φ≤ 1 }. Pursuing our interest in motion by crystalline mean curvature in three dimensions, we introduce and study the class ?_φ (ℝ ⁿ ) of “smooth” boundaries in the relative geometry induced by the ambient Banach space (ℝ ⁿ , φ). It can be seen that, even when ?_φ is a polytope, ?_φ(ℝ ⁿ ) cannot be reduced to the class of polyhedral boundaries (locally resembling ∂?_φ). Curved portions must be necessarily included and this fact (as well as the nonsmoothness of ∂?_φ) is the source of several technical difficulties related to the geometry of Lipschitz manifolds. Given a boundary δE in the class ?_φ(ℝ ⁿ ), we rigorously compute the first variation of the corresponding anisotropic perimeter, which leads to a variational problem on vector fields defined on δE. It turns out that the minimizers have a uniquely determined (intrinsic) tangential divergence on δE. We define such a divergence to be the φ-mean curvature κ_φ of δE; the function κ_φ is expected to be the initial velocity of δE, whenever δE is considered as the initial datum for the corresponding anisotropic mean curvature flow. We prove that κ_φ is bounded on δE and that its sublevel sets are characterized through a variational inequality.     

Robust H ∞ output tracking control for fuzzy networked systems with stochastic sampling and multiplicative noise

In this paper, the problem of robust sampled-data H _?? output tracking control is investigated for a class of fuzzy networked systems with stochastic sampling, multiplicative noise and time-varying norm-bounded uncertainties. For the sake of technical simplicity, only two different sampling periods are considered, their occurrence probabilities are given constants and satisfy Bernoulli distribution, and they can be extended to the case with multiple stochastic sampling periods. By using an input-delay method, the probabilistic system is transformed into a stochastic continuous time-delay system. A?new linear matrix inequality(LMI)-based procedure is proposed for designing state-feedback controllers, which would guarantee that the closed-loop networked system with stochastic sampling tracks the output of a given reference model well in the sense of H _??. Conservatism is reduced by taking the probability into account. Both network-induced delays and packet dropouts have been considered. Finally, an illustrative example is given to show the usefulness and effectiveness of the proposed H _?? output tracking design.     

Control problems of Chen–Lee system by adaptive control method

Nonlinear Dynamics - This paper investigates the control problems of Chen–Lee system. Based on control theory of nonlinear systems, the adaptive controllers were proposed to achieve...     

A Variational Model for Plastic Slip and Its Regularization via Γ-Convergence

A variational model is presented able to interpret the onset of plastic deformations, here modeled as displacement jumps occurring along slip surfaces at constant yielding stress. The corresponding strain energy functional, leading to a free-discontinuity problem set in the space of SBV functions, is then approximated by a sequence of regularized elliptic functionals following the seminal work by Ambrosio and Tortorelli (Commun. Pure Appl. Math. 43, 999–1036, 1990) within the framework of Γ-convergence. Comparisons between the results obtainable with the free-discontinuity model and its regularized approximation, in terms of stability of the pure elastic phase, irreversibility of plastic slip and response under unloading, are presented, in general, for the 2-D case of antiplane shear and exemplified, in particular, for the 1-D case.     

Variational principles and conservation laws to the Burridge–Knopoff equation

We generate conservation laws for the Burridge–Knopoff equation which model nonlinear dynamics of earthquake faults by a new conservation theorem proposed recently by Ibragimov. One can employ this new general theorem for every differential equation (or systems) and derive new local and nonlocal conservation laws. Nonlocal conservation laws comprise nonlocal variables defined by the adjoint equations to the Burridge–Knopoff equation.     

Variational Principles for EIastoplastic Continua∗

ABSTRACT

A mixed variational principle is constrained by a homogeneous yield function using a Lagrange multiplier. The Lagrange factor corresponds to the scalar factor in Prager's normality rule for the plastic strain increments. Several reduced functionals and their associated constitutive equations are derived by eliminating some variables.     

Intrinsic Geometry and Analysis of Diffusion Processes and L ∞-Variational Problems

The aim of this paper is twofold. First, we obtain a better understanding of the intrinsic distance of diffusion processes. Precisely, (a) for all n ≧ 1, the diffusion matrix A is weak upper semicontinuous on Ω if and only if the intrinsic differential and the local intrinsic distance structures coincide; (b) if n = 1, or if n ≧ 2 and A is weak upper semicontinuous on Ω, the intrinsic distance and differential structures always coincide; (c) if n ≧ 2 and A fails to be weak upper semicontinuous on Ω, the (non-)coincidence of the intrinsic distance and differential structures depend on the geometry of the non-weak-upper-semicontinuity set of A. Second, for an arbitrary diffusion matrix A, we show that the intrinsic distance completely determines the absolute minimizer of the corresponding L ^∞-variational problem, and then obtain the existence and uniqueness for given boundary data. We also give an example of a diffusion matrix A for which there is an absolute minimizer that is not of class C ¹. When A is continuous, we also obtain the linear approximation property of the absolute minimizer.     

The moiré method and the evaluation of principal-moment and stress directions

This paper discusses, in relation to the moiré method as used for the solution of plate bending and two-dimensional stress problems, two graphical techniques for the determination of the directions of principal moments and stresses.The so-called isoclinic method and the point method are described.The application of these new techniques on three different models—a circular disk under diametrically opposite loads and two different circular plates subjected to a lateral load—are fully discussed.The graphically determined principal-stress and moment directions show excellent agreement with analytically determined comparable values.Paper was presented at 1963 SESA Annual Meeting held in Boston, Mass., on November 6–8.     

Hole-drilling method using grating rosette and Moiré interferometry

The hole-drilling method is one of the most wellknown methods for measuring residual stresses. To identify unknown plane stresses in a specimen, a circular hole is first drilled in the infinite plate under plane stress, then the strains resulting from the hole drilling is measured. The strains may be acquired from interpreting the Moire signature around the hole. In crossed grating Moire interferometry, the horizontal and vertical displacement fields （u and v） can be obtained to determinate two strain fields and one shearing strain field. In this paper, by means of Moire interferometry and three directions grating （grating rosette） developed by the authors, three displacement fields （u, v and s） are obtained to acquire three strain fields. As a practical application, the hole-drilling method is adopted to measure the relief strains for aluminum and fiber reinforced composite. It is a step by step method; in each step a single laminate or equivalent depth is drilled to find some relationships between the drilling depth and the residual strains relieved in the fiber reinforced composite materials.     

Various Reciprocal Theorems and Variational Principles in the Theories of Nonlocal Micropolar Linear Elastic Mediums

In the first part of our paper,we have extended the concepts of the classical convolution and the“convolution scalar product”given by I.Hlavácěk and presented the concepts of the“convolution vector”and the“convolution vector scalar product”,which enable us to extend the initial value as well as the initial-boundary value problems for the equation with the operator coefficients to those for the system of equations with the operator coefficients.In the second part of this paper,based on the concepts of the convolution vector and the con-volution vector scalar product,two fundamental types of reciprocal theorems of the non-local micro-polar linear elastodynamics for inhomogeneous and anisotropic solids are derived.In the third part of this paper,based on the concepts and results in the first and second parts as well as the Lagrange multiplies method which is presented by W.Z.Chien,four main types of variational principles are given for the nonlocal micropolar linear elastodynamics for inhomogeneous     

On a Volume‐Constrained Variational Problem

. Existence of minimizers for a volume-constrained energy $ E(u) := \int_{\Omega} W(\nabla u)\, dx Existence of minimizers for a volume-constrained energy E(u) : = ò_W W(?u) dx E(u) := \int_{\Omega} W(\nabla u)\, dx where L^N({u = z_i}) = a_i, i = 1, ?, P, {\cal L}^N(\{u = z_i\}) = \alpha_i, i = 1, \ldots, P, is proved for the case in which z_iz_i are extremal points of a compact, convex set in \Bbb R^d\Bbb R^d and under suitable assumptions on a class of quasiconvex energy densities W. Optimality properties are studied in the scalar-valued problem where d=1d=1, P=2P=2, W(x)=|x|²W(\xi)=|\xi|^2, and the &-limit as the sum of the measures of the 2 phases tends to \L(W)\L(\Omega) is identified. Minimizers are fully characterized when N=1N=1, and candidates for solutions are studied for the circle and the square in the plane.     

An extension of Neumann’s method for shape control of force-induced elastic vibrations by eigenstrains

This paper is concerned with the force-induced vibrations of linear elastic solids and structures. We seek a transient distribution of actuating stresses produced by additional eigenstrain, such that the vibrations produced by a given set of imposed forces are exactly compensated. This problem, known as dynamic shape control problem in structural engineering, or as dynamic displacement compensation problem in automatic control, is inverse to the usual direct problem of determining displacements due to imposed forces and actuation stresses. In the present paper, we extend a method, which was introduced by F.E. Neumann for demonstrating the uniqueness of direct elastodynamic problems. We use this extended Neumann method in order to show that the distribution of the actuating stresses for shape control must be equal to any statically admissible stress distribution that is in temporal equilibrium with the imposed forces. We furthermore discuss the role of stresses corresponding to this class of solutions in some detail, emphasizing the non-unique nature of a statically admissible stress. As an analytical justification of our formulations, we show that our method reveals some static results by J.M.C. Duhamel and by W. Voigt and D.E. Carlson. Particularly, our method can be interpreted as a dynamic extension of the Duhamel body-force analogy. We moreover present numerical results for a dynamically loaded, irregularly shaped domain in a state of plane strain. These finite element computations give excellent evidence for the validity of the presented method of shape control for both, the case of a step-input and the case of a harmonic excitation.     

Numerical simulation with a TVD–FVM method for circular cylinder wake control by a fairing

A water drop-shaped fairing is applied to control the wake behind a circular cylinder and to suppress the formation of Karman vortex street in this paper. The results are evaluated using high resolution CFD technique. A finite-volume total variation diminishing (TVD) approach based upon the recently proposed elemental velocity vector transformation (EVVT) method, which aims at solving the incompressible turbulent flow for irregular boundary conditions with renormalization group (RNG) turbulence model, is used to simulate the flow field around circular cylinder systems. The calculations are carried out with cylinder systems with and without fairings, while the fairings have different top shape angles within the range of $30°~90°$. The Reynolds number ranges from 1000 to 50 000. It is shown that the simulation results of present numerical method reaches good agreement with the available experimental and numerical simulation data of typical circular cylinder flow and a fixed fairing cylinder system flow. Compared with bare cylinder, the faired bluff structures can obviously reduce the lift and drag forces and alter the vortex shedding frequency. Overall, the mean drag coefficient can be reduced up to about (10–31)% and the RMS lift coefficient can be reduced up to (30–99)% for all faired systems at given Reynolds numbers. The influence of Reynolds number and attack angles on the flow field characters of bare cylinder and faired cylinders is also discussed. The faired structures with shape angles within $30°~45°$under zero-attack-angle-inflow case are considered as the optimal structures, with which the mean drag coefficient and the RMS lift coefficient can be reduced up to (26–31)% and (98–99)%, respectively. Considering the influence of attack angles on lift and drag coefficients reduction, 75° shaped faired structure may be taken as a proper option.     

The Euler Equation and¶Absolute Minimizers of L∞ Functionals

The Aronsson-Euler equation for the functional

Weighted L∞ bounds and uniqueness for the Boltzmann BGK model

We prove rather general L bounds for hydrodynamical fields in terms of weighted L norms of the kinetic density. We use these estimates to prove L bounds and uniqueness for the solution of the BGK Equation, which is a simple relaxation model introduced by Bhatnagar, Gross & Krook to mimic Boltzmann flows.