Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ADI”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

New Techniques for High-Contrast Imaging with ADI: The ACORNS-ADI SEEDS Data Reduction Pipeline

We describe Algorithms for Calibration, Optimized Registration, and Nulling the Star in Angular Differential Imaging (ACORNS-ADI), a new, parallelized software package to reduce high-contrast imaging data, and its application to data from the Strategic Exploration of Exoplanets and Disks (SEEDS) survey. We implement seyeral new algorithms, includbg a method to centroid saturated images, a trimmed mean for combining an image sequence that reduces noise by up to approx 20%, and a robust and computationally fast method to compute the sensitivitv of a high-contrast obsen-ation everywhere on the field-of-view without introducing artificial sources. We also include a description of image processing steps to remove electronic artifacts specific to Hawaii2-RG detectors like the one used for SEEDS, and a detailed analysis of the Locally Optimized Combination of Images (LOCI) algorithm commonly used to reduce high-contrast imaging data. ACORNS-ADI is efficient and open-source, and includes several optional features which may improve performance on data from other instruments. ACORNS-ADI is freely available for download at www.github.com/t-brandt/acorns_-adi under a BSD license

Brandt, Timothy D.↗

ADI on staggered mesh - A method for the calculation of compressible convection

An alternating direction implicit (ADI) method has been applied to a staggered grid for the computation of convection in a highly stratified fluid. Since artificial viscosity is not needed, subtle effects like the onset of convection can be studied. These computations compare well with the 2-D results by Graham (1975) and also agree with standard Boussinesq results when taken to that limit. Good efficiency has been achieved with a time step hundreds of times larger than the stability limit imposed by the explicit treatment of diffusion and the Courant number is not restricted to be below 1. The Navier-Stokes equation contains cross spatial derivatives which are treated explicitly in most ADI schemes. The destabilizing effect of such a practice on a 2-D model system with second-order spatial derivative terms only was analyzed and found to be not excessive.

Chan, K. L.↗

Improving the convergence rate of parabolic ADI methods

The rate of convergence to steady state of parabolic Alternating Direction Implicit (ADI) solvers is analyzed in terms of the L(2)-norms of the residuals. The analysis allows one to predict the number of iterations necessary for convergence as function of the Courant number, Lambda. A simple modification of existing ADI codes is devised. It improves the convergence rate substantially and is insensitive to the Courant number in a large range of Lambda.

Abarbanel, S. S.↗

Improving the convergence rate to steady state of parabolic ADI methods

The present, residuals' L(2)-norms analysis of the rate of convergence to steady state for parabolic ADI solvers allows the prediction of the number of iterations required for convergence, as a function of the Courant number alpha. A modification of current ADI codes is presented which significantly improves the convergence rate and is insensitive to the Courant number over a large range of alpha. This corrected algorithm is tested for the cases of Dirichlet problems for uniform grids of many mesh sizes, mixed Dirichlet-Neumann problems, and problems defined on stretched grids and/or problems with variable coefficients.

Abarbanel, Saul S.↗

Second-order non-iterative ADI solution of non-linear partial differential equations

A new method for the solution of non-linear partial differential equations by an ADI procedure is described. Although the method is second order accurate in time, it does not require either iterations or predictor corrector methods to overcome the nonlinearity of the equations. Thus the computational effort required for the solution of the non-linear problem becomes similar to that required for the linear case. The method is applied to a two-dimensional 'extended Burgers equation'. Linear stability is studied, and some numerical solutions obtained. The improved accuracy obtained by the 2nd order truncation error is clearly manifested.

Wolfshtein, M.↗

Intermediate boundary conditions for LOD, ADI and approximate factorization methods

A general approach to determining the correct intermediate boundary conditions for dimensional splitting methods is presented. The intermediate solution U is viewed as a second order accurate approximation to a modified equation. Deriving the modified equation and using the relationship between this equation and the original equation allows us to determine the correct boundary conditions for U*. This technique is illustrated by applying it to locally one dimensional (LOD) and alternating direction implicit (ADI) methods for the heat equation in two and three space dimensions. The approximate factorization method is considered in slightly more generality.

Leveque, R. J.↗

Implementation of an ADI method on parallel computers

The implementation of an ADI method for solving the diffusion equation on three parallel/vector computers is discussed. The computers were chosen so as to encompass a variety of architectures. They are: the MPP, an SIMD machine with 16K bit serial processors; FLEX/32, an MIMD machine with 20 processors; and CRAY/2, an MIMD machine with four vector processors. The Gaussian elimination algorithm is used to solve a set of tridiagonal systems on the FLEX/32 and CRAY/2 while the cyclic elimination algorithm is used to solve these systems on the MPP. The implementation of the method is discussed in relation to these architectures and measures of the performance on each machine are given. Simple performance models are used to describe the performance. These models highlight the bottlenecks and limiting factors for this algorithm on these architectures. Finally, conclusions are presented.

Fatoohi, Raad A.↗

Navier-Stokes simulation of wind-tunnel flow using LU-ADI factorization algorithm

The three dimensional Navier-Stokes solution code using the LU-ADI factorization algorithm was employed to simulate the workshop test cases of transonic flow past a wing model in a wind tunnel and in free air. The effect of the tunnel walls is well demonstrated by the present simulations. An Amdahl 1200 supercomputer having 128 Mbytes main memory was used for these computations.

Obayashi, Shigeru↗

Implementation of an ADI method on parallel computers

In this paper the implementation of an ADI method for solving the diffusion equation on three parallel/vector computers is discussed. The computers were chosen so as to encompass a variety of architectures. They are the MPP, an SIMD machine with 16-Kbit serial processors; Flex/32, an MIMD machine with 20 processors; and Cray/2, an MIMD machine with four vector processors. The Gaussian elimination algorithm is used to solve a set of tridiagonal systems on the Flex/32 and Cray/2 while the cyclic elimination algorithm is used to solve these systems on the MPP. The implementation of the method is discussed in relation to these architectures and measures of the performance on each machine are given. Simple performance models are used to describe the performance. These models highlight the bottlenecks and limiting factors for this algorithm on these architectures. Finally conclusions are presented.

Fatoohi, Raad A.↗

Navier-Stokes simulation of wind-tunnel flow using LU-ADI factorization algorithm

The three dimensional Navier-Stokes solution code using the LU-ADI factorization algorithm was employed to simulate the workshop test cases of transonic flow past a wing model in a wind tunnel and in free air. The effect of the tunnel walls is well demonstrated by the present simulations. An Amdahl 1200 supercomputer having 128 Mbytes main memory was used for these computations.

Obayashi, Shigeru↗

Implementation of ADI: Schemes on MIMD parallel computers

In order to simulate the effects of the impingement of hot exhaust jets of High Performance Aircraft on landing surfaces a multi-disciplinary computation coupling flow dynamics to heat conduction in the runway needs to be carried out. Such simulations, which are essentially unsteady, require very large computational power in order to be completed within a reasonable time frame of the order of an hour. Such power can be furnished by the latest generation of massively parallel computers. These remove the bottleneck of ever more congested data paths to one or a few highly specialized central processing units (CPU's) by having many off-the-shelf CPU's work independently on their own data, and exchange information only when needed. During the past year the first phase of this project was completed, in which the optimal strategy for mapping an ADI-algorithm for the three dimensional unsteady heat equation to a MIMD parallel computer was identified. This was done by implementing and comparing three different domain decomposition techniques that define the tasks for the CPU's in the parallel machine. These implementations were done for a Cartesian grid and Dirichlet boundary conditions. The most promising technique was then used to implement the heat equation solver on a general curvilinear grid with a suite of nontrivial boundary conditions. Finally, this technique was also used to implement the Scalar Penta-diagonal (SP) benchmark, which was taken from the NAS Parallel Benchmarks report. All implementations were done in the programming language C on the Intel iPSC/860 computer.

Vanderwijngaart, Rob F.↗

Flutter analysis of a supersonic cascade in time domain using an ADI Euler solver

The aeroelastic stability of a two-dimensional cascade oscillating in supersonic axial flow is analyzed in the time domain. The aeroelastic model consists of a single degree of freedom typical section structural model for each blade of the cascade and an unsteady two-dimensional cascade aerodynamic model based on the Euler equations. The Euler equations are solved using a time accurate Alternating Direction Implicit (ADI) solution scheme. The aeroelastic equations are integrated in time. The effect of interblade phase angle is included in the aeroelastic analysis by an appropriate choice of initial and boundary conditions. Flutter predictions are obtained from the time response of a flat plate cascade in single degree of freedom pitching motion. The results correlate well with those obtained from a separate frequency domain flutter analysis for all values of interblade phase angles considered. Flutter results are then presented for cascades having airfoil sections representative of a supersonic throughflow fan. The validity of the time integration method for a cascade of airfoils at various interblade phase angles is demonstrated.

Reddy, T. S. R.↗

Multi-partitioning for ADI-schemes on message passing architectures

A kind of discrete-operator splitting called Alternating Direction Implicit (ADI) has been found to be useful in simulating fluid flow problems. In particular, it is being used to study the effects of hot exhaust jets from high performance aircraft on landing surfaces. Decomposition techniques that minimize load imbalance and message-passing frequency are described. Three strategies that are investigated for implementing the NAS Scalar Penta-diagonal Parallel Benchmark (SP) are transposition, pipelined Gaussian elimination, and multipartitioning. The multipartitioning strategy, which was used on Ethernet, was found to be the most efficient, although it was considered only a moderate success because of Ethernet's limited communication properties. The efficiency derived largely from the coarse granularity of the strategy, which reduced latencies and allowed overlap of communication and computation.

Vanderwijngaart, Rob F.↗

Gradient-Based Aerodynamic Shape Optimization Using ADI Method for Large-Scale Problems

A gradient-based shape optimization methodology, that is intended for practical three-dimensional aerodynamic applications, has been developed. It is based on the quasi-analytical sensitivities. The flow analysis is rendered by a fully implicit, finite volume formulation of the Euler equations.The aerodynamic sensitivity equation is solved using the alternating-direction-implicit (ADI) algorithm for memory efficiency. A flexible wing geometry model, that is based on surface parameterization and platform schedules, is utilized. The present methodology and its components have been tested via several comparisons. Initially, the flow analysis for for a wing is compared with those obtained using an unfactored, preconditioned conjugate gradient approach (PCG), and an extensively validated CFD code. Then, the sensitivities computed with the present method have been compared with those obtained using the finite-difference and the PCG approaches. Effects of grid refinement and convergence tolerance on the analysis and shape optimization have been explored. Finally the new procedure has been demonstrated in the design of a cranked arrow wing at Mach 2.4. Despite the expected increase in the computational time, the results indicate that shape optimization, which require large numbers of grid points can be resolved with a gradient-based approach.

Pandya, Mohagna J.↗

Shared Memory Parallelization of an Implicit ADI-type CFD Code

A parallelization study designed for ADI-type algorithms is presented using the OpenMP specification for shared-memory multiprocessor programming. Details of optimizations specifically addressed to cache-based computer architectures are described and performance measurements for the single and multiprocessor implementation are summarized. The paper demonstrates that optimization of memory access on a cache-based computer architecture controls the performance of the computational algorithm. A hybrid MPI/OpenMP approach is proposed for clusters of shared memory machines to further enhance the parallel performance. The method is applied to develop a new LES/DNS code, named LESTool. A preliminary DNS calculation of a fully developed channel flow at a Reynolds number of 180, Re(sub tau) = 180, has shown good agreement with existing data.

Hauser, Th.↗

Six-degree-of-freedom missile simulation using the ADI AD 100 digital computer and ADSIM simulation language

The use of an AD 100 computer and the ADSIM language in the six-degree-of-freedom digital simulation of an air-to-ground missile is illustrated. The missile is launched from a moving platform, typically a helicopter, and is capable of striking a mobile target up to 10 kilometers away. The missile could be any tactical missile. The performance numbers of the AD 100 show that it is possible to implement a high performance missile model in a real-time simulation without the problems associated with an implementation on a general purpose computer using FORTRAN.

Zwaanenburg, Koos↗

Scalability and Portability of Two Parallel Implementations of ADI

Two domain decompositions for the implementation of the NAS Scalar Penta-diagonal Parallel Benchmark on MIMD systems are investigated, namely transposition and multi-partitioning. Hardware platforms considered are the Intel iPSC/860 and Paragon XP/S-15, and clusters of SGI workstations on ethernet, communicating through PVM. It is found that the multi-partitioning strategy offers the kind of coarse granularity that allows scaling up to hundreds of processors on a massively parallel machine. Moreover, efficiency is retained when the code is ported verbatim (save message passing syntax) to a PVM environment on a modest size cluster of workstations.

Phung, Thanh↗