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.     

抽吸和压力梯度在层流边界层转捩过程中的作用

用空间模式的二次稳定性理论研究了抽吸和压力梯度对边界层三维亚谐扰动流动稳定性的影响．数值结果表明，固体边界上的抽吸有明显的层流控制作用，逆压梯度则有较强的不稳定作用．     

Optimality conditions and duality theory for minimizing sums of the largest eigenvalues of symmetric matrices

The sum of the largest eigenvalues of a symmetric matrix is a nonsmooth convex function of the matrix elements. Max characterizations for this sum are established, giving a concise characterization of the subdifferential in terms of a dual matrix. This leads to a very useful characterization of the generalized gradient of the following convex composite function: the sum of the largest eigenvalues of a smooth symmetric matrix-valued function of a set of real parameters. The dual matrix provides the information required to either verify first-order optimality conditions at a point or to generate a descent direction for the eigenvalue sum from that point, splitting a multiple eigenvalue if necessary. Connections with the classical literature on sums of eigenvalues and eigenvalue perturbation theory are discussed. Sums of the largest eigenvalues in the absolute value sense are also addressed.This paper is dedicated to Phil Wolfe on the occasion of his 65th birthday.The work of this author was supported by the National Science Foundation under grants CCR-8802408 and CCR-9101640.The work of this author was supported in part during a visit to Argonne National Laboratory by the Applied Mathematical Sciences subprogram of the Office of Energy Research of the U.S. Department of Energy under contract W-31-109-Eng-38, and in part during a visit to the Courant Institute by the U.S. Department of Energy under Contract DEFG0288ER25053.     

A note on an SOR-like method for augmented systems

Golub et al. (2001, BIT, 41, 71–85) gave a generalizedsuccessive over-relaxation method for the augmented systems.In this paper, the connection between the SOR-like method andthe preconditioned conjugate gradient (PCG) method for the augmentedsystems is investigated. It is shown that the PCG method isat least as accurate (fast) as the SOR-like method. Numericalexamples demonstrate that the PCG method is much faster thanthe SOR-like method.     

精确搜索下的非线性共轭梯度法

该文提出一种无约束优化非线性共轭梯度法，证明了精确线性 搜索下的全局收敛性。当目标函数为一致凸函数时，证明了算法具有线性收敛速度。数值实验表明算法对于求解实际问题是有效的。     

Optimal shape design for packaging containing heating elements by inverse heat transfer method

The present study is concerned with the determination of the optimal shape for a package containing multiple heating elements. The optimization tool has been developed based on the inverse heat transfer (IHT) approach, incorporating a direct problem solver, a numerical grid generator, a direct-differentiation sensitivity analyzer, and the conjugate gradient method. Shape design that leads to a specified outer surface temperature distribution is predicted by the approach. In this study, the effects of internal heat generation on optimal shapes of the packagings have also been evaluated. Several practical cases with various imbedded heating elements and thermal conditions are studied. Results show that the approach provides an efficient computer-aided design scheme for the shape profile determination.     

一类共轭梯度法的全局收敛性结果

本文证明了在Grippo-Lucidi线搜索下当βk取βk=σ1βPRPk+σ2βnewk,其中σ1≥0,σ2≥0,σ1+σ2＞0,βnewk=gTk(gk-gk-1)/-dTk-1gk-1时一类共轭梯度法的全局收敛性,并给出了此类方法良好的数值效果.     

Flow Search Approach and New Bounds for the m-Step Linear Conjugate Gradient Algorithm

A flow search approach is presented in this paper. In the approach, each iterative process involves a subproblem, whose variables are the stepsize parameters. Every feasible solution of the subproblem corresponds to some serial search stages, the stepsize parameters in different search stages may interact mutually, and their optimal values are determined by evaluating the total effect of the interaction. The main idea of the flow search approach is illustrated via the minimization of a convex quadratic function. Based on the flow search approach, some properties of the m-step linear conjugate gradient algorithm are analyzed and new bounds on its convergence rate are also presented. Theoretical and numerical results indicate that the new bounds are better than the well-known ones.     

Elliptic systems with nonlinearity q greater or equal than two with controlled growth. Global Hölder continuity of solutions to the Dirichlet problem

Hölder regularity up to the boundary of the solutions to the Dirichlet problem for second order elliptic systems with nonlinearity q>2 and with controlled growth is proved when n?q+2.     

Shifting the Superparamagnetic Limit of Nanosized Copper Iron Spinel

In this paper we analyze trend of EFG values measured at Cd impurity in a group of semiconducting delafossites with chemical formula CuBO₂ (B = Al, Fe, Cr, Nd). We conclude that this trend reveals one of the most subtle details in electronic spectrum of the compounds: if impurity states are formed within or out of the band gap. In CuAlO₂ and CuFeO₂ the Cd EFG exhibits larger value than in CuCrO₂ and CuNdO₂, when Cd substitutes the Cu atom. This occurs because in the first two compounds the Cd forms shallow band within the gap, and in the second two compounds does not. When Cd occupies the B position it exhibits almost the same EFG in all delafossites. In this case, Cd does not form its states within the gap in none of the compounds. To arrive to these conclusions we analyzed and calculated various systems (Cd-doped CuAlO₂ and CuCrO₂ compounds, fictitious molecules), using the FP-LAPW method.