Stable multivariate rational interpolation for parameter-dependent aerospace models

A multivariate, robust, rational interpolation method for propagating uncertainties in several dimensions is presented. The algorithm for selecting numerator and denominator polynomial orders is based on recent work that uses a singular value decomposition approach. In this paper we extend this algorithm to higher dimensions and demonstrate its efficacy in terms of convergence and accuracy, both as a method for response suface generation and interpolation. To obtain stable approximants for continuous functions, we use an L2 error norm indicator to rank optimal numerator and denominator solutions. For discontinous functions, a second criterion setting an upper limit on the approximant value is employed. Analytical examples demonstrate that, for the same stencil, rational methods can yield more rapid convergence compared to pseudospectral or collocation approaches for certain problems. © 2012 AIAA.

Models for acoustic propagation through turbofan exhaust flows

Turbomachinery noise radiating into the rearward arc is an important problem. This noise is scattered by the trailing edges of the nacelle and the jet exhaust, and interacts with the shear layers between the external flow, bypass stream and jet, en route to the far field. In the past a range of relevant model problems involving semi-infinite cylinders have been solved. However, one limitation of previous solutions is that they do not allow for the jet nozzle to protrude a finite distance beyond the end of the nacelle (or in certain configurations being buried a finite distance upstream). In this paper we use the matrix Wiener-Hopf technique, which will allow precisely the finite nacelle-jet nozzle separation to be included. The crucial step in our work is to factorise a certain matrix as a product of terms analytic and invertible in the upper/lower halves of the complex plane. The way we do this matrix factorisation is quite different in the buried and protruding nozzle cases. In the buried case our solution method is the so-called pole-removal technique. In the technically more demanding protruding case, however, we must first use Pade approximants to generate a uniformly-valid, meromorphic representation of a certain function, before the same pole-removal method can be applied. Sample results are presented, investigating in particular the effects of exit plane stagger. © 2007 by B Veitch and N Peake.

Electrochemical hydrogen storage in (Ti1-xVx)(2)Ni (x=0.05-0.3) alloys comprising icosahedral quasicrystalline phase

For (Ti1-xVx)(2)Ni (x = 0.05,0.1,0.15,0.2 and 0.3) ribbons, synthesized by arc-melting and subsequent melt-spinning techniques, an icosahedral quasicrystalline phase was present, either in the amorphous matrix or together with the stable Ti2Ni-type phase. With increasing x values, the maximum discharge capacity of the alloy electrodes increased until reached 271.3 mAh/g when x = 0.3. The cycling capacity retention rates for these electrodes were approximately 80% after a preliminary test of 30 consecutive cycles of charging and discharging.

Electrochemical Hydrogen Storage in Ti1.6V0.4Ni1-xCox Icosahedral Quasicrystalline Alloys

The discovery of the icosahedral phase (i-phase) in rapidly quenched Ti1.6V0.4Ni1-xCox (x=0.02-01) alloys is described herein. The i-phase occurs in a similar amount relative to the coexisting beta Ti phase. The electron diffraction patterns show the distinct spot anisotropy, indicating that the i-phase is metastable. The electrochemical hydrogen storage performance of these five alloy electrodes are also reported herein. The hydrogen desorption of nonelectrochemical recombination in the cyclic voltammetric (CV) response exhibits the demand for electrocatalytic activity improvement.

Preparation and transport properties of a bulk icosahedral quasicrystalline Ti45Zr35Ni17Cu3 alloy

A bulk alloy which consists of the single icosahedral quasicrystalline phase (I-phase) in Ti45Zr35Ni17CU3 alloy has been fabricated by mechanical alloying and subsequent pulse discharge sintering technique. Crystallographic structure analyses show that the bulk alloy is an I-phase. The transport properties of the bulk alloy are examined, and the results show that the room-temperature thermal conductivity is 5.347 W K-(1) m(-1), and the electrical conductivity decreases with increasing the temperature from 300 to 450K. The Seebeck coefficient is negative at the temperature range from 300 to 360K, and changes to positive from 370 to 450K. Hall effect measurements indicate the bulk I-phase alloy has a high carrier concentration. The specific heat capacity increases when the temperature increases from 280 to 324 K.

Semiclassical complex angular momentum theory and Pade reconstruction for resonances, rainbows, and reaction thresholds

A semiclassical complex angular momentum theory, used to analyze atom-diatom reactive angular distributions, is applied to several well-known potential (one-particle) problems. Examples include resonance scattering, rainbow scattering, and the Eckart threshold model. Pade reconstruction of the corresponding matrix elements from the values at physical (integral) angular momenta and properties of the Pade approximants are discussed in detail.

Bose systems in spatially random or time-varying potentials

Bose systems, subject to the action of external random potentials, are considered. For describing the system properties, under the action of spatially random potentials of arbitrary strength, the stochastic mean-field approximation is employed. When the strength of disorder increases, the extended Bose-Einstein condensate fragments into spatially disconnected regions, forming a granular condensate. Increasing the strength of disorder even more transforms the granular condensate into the normal glass. The influence of time-dependent external potentials is also discussed. Fastly varying temporal potentials, to some extent, imitate the action of spatially random potentials. In particular, strong time-alternating potential can induce the appearance of a nonequilibrium granular condensate.

The Q-D algorithm for transforming series expansions into a corresponding continued fraction: an extension to cope with zero coefficients

An algorithm for deriving a continued fraction that corresponds to two series expansions simultaneously, when there are zero coefficients in one or both series, is given. It is based on using the Q-D algorithm to derive the corresponding fraction for two related series, and then transforming it into the required continued fraction. Two examples are given. (C) 2003 Elsevier B.V. All rights reserved.

Application of a meshless method in electromagnetics

An improved meshless method is presented with an emphasis on the detailed description of this new computational technique and its numerical implementations by investigating the usefulness of a commonly neglected parameter in this paper. Two approaches to enforce essential boundary conditions are also thoroughly investigated. Numerical tests on a mathematical function is carried out as a means of validating the proposed method. It will be seen that the proposed method is more robust than the conventional ones. Applications in solving electromagnetic problems are also presented.

The application of interpolating MLS approximations to the analysis of MHD flows

The element-free Galerkin method (EFGM) is a very attractive technique for solutions of partial differential equations, since it makes use of nodal point configurations which do not require a mesh. Therefore, it differs from FEM-like approaches by avoiding the need of meshing, a very demanding task for complicated geometry problems. However, the imposition of boundary conditions is not straightforward, since the EFGM is based on moving-least-squares (MLS) approximations which are not necessarily interpolants. This feature requires, for instance, the introduction of modified functionals with additional unknown parameters such as Lagrange multipliers, a serious drawback which leads to poor conditionings of the matrix equations. In this paper, an interpolatory formulation for MLS approximants is presented: it allows the direct introduction of boundary conditions, reducing the processing time and improving the condition numbers. The formulation is applied to the study of two-dimensional magnetohydrodynamic flow problems, and the computed results confirm the accuracy and correctness of the proposed formulation. (C) 2002 Elsevier B.V. All rights reserved.

Some Extensions of M-Fractions Related to Strong Stieltjes Distributions

The purpose of this paper is to show the symmetric relations that appear between the coefficients of some even and odd extensions of the M-fractions related to a certain kind of symmetric strong Stieltjes distribution.

On the two point Padé table for a distribution

Some additional recurrence relations for the denominator polynomials of two point Padé approximants are derived. An example in which the coefficients of one of the two series, from which the Padé approximants are derived, are moments of a distribution is considered. For this example, properties of the denominator polynomials, and their zeros, are described.

Frações contínuas que correspondem a séries de potências em dois pontos

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Meromorphic continuation of functions and arbitrary distribution of interpolation points

We characterize the region of meromorphic continuation of an analytic function ff in terms of the geometric rate of convergence on a compact set of sequences of multi-point rational interpolants of ff. The rational approximants have a bounded number of poles and the distribution of interpolation points is arbitrary.

Método da propagação de feixe de ângulo largo para análise de guias de ondas ópticos não-lineares

Este trabalho propõe uma extensão do método de propagação de feixe (BPM - Beam Propagation Method) para a análise de guias de ondas ópticos e acopladores baseados em materiais não-lineares do tipo Kerr. Este método se destina à investigação de estruturas onde a utilização da equação escalar de Helmholtz (EEH) em seu limite paraxial não mais se aplica. Os métodos desenvolvidos para este fim são denominados na literatura como métodos de propagação de feixe de ângulo largo. O formalismo aqui desenvolvido é baseado na técnica das diferenças finitas e nos esquemas de Crank-Nicholson (CN) e Douglas generalizado (GD). Estes esquemas apresentam como característica o fato de apresentarem um erro de truncamento em relação ao passo de discretização transversal, Δx, proporcional a O(Δx2) para o primeiro e O(Δx4). A convergência do método em ambos esquemas é otimizada pela utilização de um algoritmo interativo para a correção do campo no meio não-linear. O formalismo de ângulo largo é obtido pela expansão da EEH para os esquemas CN e GD em termos de polinômios aproximantes de Padé de ordem (1,0) e (1,1) para CN e GD, e (2,2) e (3,3) para CN. Os aproximantes de ordem superior a (1,1) apresentam sérios problemas de estabilidade. Este problema é eliminado pela rotação dos aproximantes no plano complexo. Duas condições de contorno nos extremos da janela computacional são também investigadas: 1) (TBC - Transparent Boundary Condition) e 2) condição de contorno absorvente (TAB - Transparent Absorbing Boundary). Estas condições de contorno possuem a facilidade de evitar que reflexões indesejáveis sejam transmitidas para dentro da janela computacional. Um estudo comparativo da influência destas condições de contorno na solução de guias de ondas ópticos não-lineares é também abordada neste trabalho.