A cooperative conjugate gradient method for linear systems permitting multithread implementation of low complexity

This paper proposes a generalization of the conjugate gradient (CG) method used to solve the equation $Ax=b$ for a symmetric positive definite matrix $A$ of large size $n$. The generalization consists of permitting the scalar control parameters (= stepsizes in gradient and conjugate gradient directions) to be replaced by matrices, so that multiple descent and conjugate directions are updated simultaneously. Implementation involves the use of multiple agents or threads and is referred to as cooperative CG (cCG), in which the cooperation between agents resides in the fact that the calculation of each entry of the control parameter matrix now involves information that comes from the other agents. For a sufficiently large dimension $n$, the use of an optimal number of cores gives the result that the multithread implementation has worst case complexity $O(n^{2+1/3})$ in exact arithmetic. Numerical experiments, that illustrate the interest of theoretical results, are carried out on a multicore computer.Comment: Expanded version of manuscript submitted to the IEEE-CDC 2012 (Conference on Decision and Control

The benefit of high-conductivity materials in film cooled turbine nozzles

This study presents an experimental and numerical investigation of the beneficial effect of higher conductivity materials in HP turbine nozzles. Most of the literature studies focus on the maximum temperature that a nozzle can withstand, whereas the effect of thermal gradients is often neglected. However thermal gradients have higher influence on the life of the components and they have to be given careful consideration. In this work it is shown that thermal gradients are reduced by using high conductivity materials and, as a consequence, the nozzles life is appreciably increased. A representative film cooled leading edge with an internal impingement plate was studied experimentally at Texas AM University. Two materials were used, namely polycarbonate and stainless steel, in order to highlight the impact of conduction on coolant effectiveness. Numerically conjugate heat transfer simulations have been carried out with an in house solver to analyse in detail the impact of conduction and internal convection. Both experimental and numerical results show that by increasing the conductivity in the solid region, the thermal gradients are strongly reduced. Numerically it is shown that using inserts of nickel-aluminide alloys in nozzles may reduce the thermal gradients from 3 to 4 times if compared to nowadays design. © 2012 Elsevier Inc

Direct Minimization Generating Electronic States with Proper Occupation Numbers

We carry out the direct minimization of the energy functional proposed by Mauri, Galli and Car to derive the correct self-consistent ground state with fractional occupation numbers for a system degenerating at the Fermi level. As a consequence, this approach enables us to determine the electronic structure of metallic systems to a high degree of accuracy without the aid of level broadening of the Fermi-distribution function. The efficiency of the method is illustrated by calculating the ground-state energy of C$_2$ and Si$_2$ molecules and the W(110) surface to which a tungsten adatom is adsorbed.Comment: 4 pages, 4 figure

Postprocessing of Non-Conservative Flux for Compatibility with Transport in Heterogeneous Media

A conservative flux postprocessing algorithm is presented for both steady-state and dynamic flow models. The postprocessed flux is shown to have the same convergence order as the original flux. An arbitrary flux approximation is projected into a conservative subspace by adding a piecewise constant correction that is minimized in a weighted $L^2$ norm. The application of a weighted norm appears to yield better results for heterogeneous media than the standard $L^2$ norm which has been considered in earlier works. We also study the effect of different flux calculations on the domain boundary. In particular we consider the continuous Galerkin finite element method for solving Darcy flow and couple it with a discontinuous Galerkin finite element method for an advective transport problem.Comment: 34 pages, 17 figures, 11 table

Computational alternatives to obtain time optimal jet engine control

Two computational methods to determine an open loop time optimal control sequence for a simple single spool turbojet engine are described by a set of nonlinear differential equations. Both methods are modifications of widely accepted algorithms which can solve fixed time unconstrained optimal control problems with a free right end. Constrained problems to be considered have fixed right ends and free time. Dynamic programming is defined on a standard problem and it yields a successive approximation solution to the time optimal problem of interest. A feedback control law is obtained and it is then used to determine the corresponding open loop control sequence. The Fletcher-Reeves conjugate gradient method has been selected for adaptation to solve a nonlinear optimal control problem with state variable and control constraints

Optimal switching of a nanomagnet assisted by microwaves

We develop an efficient and general method for optimizing the microwave field that achieves magnetization switching with a smaller static field. This method is based on optimal control and renders an exact solution for the 3D microwave field that triggers the switching of a nanomagnet with a given anisotropy and in an oblique static field. Applying this technique to the particular case of uniaxial anisotropy, we show that the optimal microwave field, that achieves switching with minimal absorbed energy, is modulated both in frequency and in magnitude. Its role is to drive the magnetization from the metastable equilibrium position towards the saddle point and then damping induces the relaxation to the stable equilibrium position. For the pumping to be efficient, the microwave field frequency must match at the early stage of the switching process the proper precession frequency of the magnetization, which depends on the magnitude and direction of the static field. We investigate the effect of the static field (in amplitude and direction) and of damping on the characteristics of the microwave field. We have computed the switching curves in the presence of the optimal microwave field. The results are in qualitative agreement with micro-SQUID experiments on isolated nanoclusters. The strong dependence of the microwave field and that of the switching curve on the damping parameter may be useful in probing damping in various nanoclusters.Comment: 9 pages, 8 figure