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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Numerical analysis of natural convection in liquid droplets by phase change

A numerical analysis is performed on thermocapillary buoyancy convection induced by phase change in a liquid droplet. A finite-difference code is developed using an alternating-direction implicit (ADI) scheme. The intercoupling relation between thermocapillary force, buoyancy force, fluid property, heat transfer, and phase change, along with their effects on the induced flow patterns, are disclosed. The flow is classified into three types: thermocapillary, buoyancy, and combined convection. Among the three mechanisms, the combined convection simulates the experimental observations quite well, and the basic mechanism of the observed convection inside evaporating sessile drops is thus identified. It is disclosed that evaporation initiates unstable convection, while condensation always brings about a stable density distribution which eventually damps out all fluid disturbances. Another numerical model is presented to study the effect of boundary recession due to evaporation, and the 'peeling-off' effect (the removal of the surface layer of fluid by evaporation) is shown to be relevant.

Duh, J. C.↗
NASA NTRS1991

Navier-Stokes Simulation Of Wind-Tunnel Flow

Report describes three-dimensional Navier-Stokes simulation of transonic flow past model wing in wind tunnel and in free air. Numerical integration performed by LU-ADI factorization algorithm. Simulation shows effects of walls of tunnel.

Obayashi, S.↗
NASA NTRS1989

Numerical analysis of flow through oscillating cascade sections

The design of turbomachinery blades requires the prevention of flutter for all operating conditions. However, flow field predictions used for aeroelastic analysis are not well understood for all flow regimes. The present research focuses on numerical solutions of the Euler and Navier-Stokes equations using an ADI procedure to model two-dimensional, transonic flow through oscillating cascades. The model prescribes harmonic pitching motions for the blade sections for both zero and non-zero inter-blade phase angles. The code introduces the use of a deforming grid technique for convenient specification of the periodic boundary conditions. Approximate nonreflecting boundary conditions have been coded for the inlet and exit boundary conditions. Sample unsteady solutions have been performed for an oscillating cascade and compared to experimental data. Also, test cases were fun for a flat plate cascade to compare with an unsteady, small-perturbation, subsonic analysis. The predictions for oscillating cascades with non-zero inter-blade phase angles are in good agreement with experimental data and small-perturbation theory. The zero degree inter-blade phase angle cases, which were near a resonant condition, differ from the experiment and theory. Studies on reflecting versus non-reflecting inlet and exit boundary conditions show that the treatment of the boundary can have a significant effect on the first harmonic, unsteady pressure distributions for certain flow conditions. This code is expected to be used as a tool for reviewing simpler models that do not include the full nonlinear aerodynamics or as a final check for designs against flutter in turbomachinery.

Huff, Dennis L.↗
NASA NTRS1989

Aerodynamic stability and heating analyses for the aeroassist flight experiment vehicle

Since ground based flow simulations are presently unable to model flight conditions expected for AOTVs (Aeroassist Orbital Transfer Vehicle) and other hypersonic space vehicles, computer codes are being developed to provide design parameters necessary for structure, guidance, and control aspects. Over the past four years, VRFLO (Viscous Reactive Flow) has been written to model finite-rate chemistry and viscous effects for a variety of aerobrake bodies. VRFLO includes a number of unique features that are summarized as follows: (1) Grid generation is an integral part of the code for several aerobrake configurations which includes the wake flow region; (2) The formulation is valid for three air chemical models; (3) An ADI central difference technique is used to solve the Navier-Stokes and species continuity equations in split groups; and (4) Grid density and numerical damping are minimized by shock-fitting and conformal mapping of body points.

Mcgary, J.↗
NASA NTRS1990

Transonic blade-vortex interactions noise: A parametric study

Transonic Blade-Vortex Interactions (BVI) are simulated numerically and the noise mechanisms are investigated. The 2-D high frequency transonic small disturbance equation is solved numerically (VTRAN2 code). An Alternating Direction Implicit (ADI) scheme with monotone switches is used; viscous effects are included on the boundary and the vortex is simulated by the cloud-in-cell method. The Kirchoff method is used for the extension of the numerical 2-D near field aerodynamic results to the linear acoustic 3-D far field. The viscous effect (shock/boundary layer interaction) on BVI is investigated. The different types of shock motion are identified and compared. Two important disturbances with different directivity exist in the pressure signal and are believed to be related to the fluctuating lift and drag forces. Noise directivity for different cases is shown. The maximum radiation occurs at an angle between 60 and 90 deg below the horizontal for an airfoil fixed coordinate system and depends on the details of the airfoil shape. Different airfoil shapes are studied and classified according to the BVI noise produced.

Lyrintzis, A. S.↗
NASA NTRS1990

Computational results for flows over 2-D ramp and 3-D obstacle with an upwind Navier-Stokes solver

An implicit, finite-difference, upwind, full Navier-Stokes solver was applied to supersonic/hypersonic flows over two-dimensional ramps and three-dimensional obstacle. Some of the computed results are presented. The numerical scheme used in the study is an implicit, spacially second order accurate, upwind, LU-ADI scheme based on Roe's approximate Reimann solver with MUSCL differencing of Van Leer. An algebraic grid generation scheme based on generalized interpolation scheme was used in generating the grids for the various 2-D and 3-D problems.

Venkatapathy, Ethiraj↗
NASA NTRS1990

A study of the noise mechanisms of transonic blade-vortex interactions

Transonic blade-vortex interactions (BVI) are simulated numerically and the noise mechanisms are investigated. The two-dimensional high frequency transonic small disturbance equation is solved numerically (VTRAN2 code). An ADI scheme with monotone switches is used; viscous effects are included on the boundary, and the vortex is simulated by the cloud in cell method. The Kirchhoff method is used for the extension of the numerical two-dimensional near-field aerodynamic results to the linear acoustic three dimensional far field. The viscous effects (shock/boundary layer interactions) on BVI is investigated. The different types of shock motion are identified and compared. Two important disturbances with different directivity exist in the pressure signal and are believed to be related to the fluctuating lift and drag forces. Noise directivity for different cases is shown. The maximum radiation occurs at an angle between 60 and 90 degrees below the horizontal for an airfoil-fixed coordinate system and depends on the details of the airfoil shape. Different airfoil shapes are studied and classified according to the BVI noise produced.

Lyrintzis, Anastasios S.↗
NASA NTRS1989

Compact finite volume methods for the diffusion equation

The paper describes an approach to treating initial-boundary-value problems by finite volume methods in which the parallel between differential and difference arguments is closely maintained. By using intrinsic geometrical properties of the volume elements, it is possible to describe discrete versions of the div, curl, and grad operators which lead, using summation-by-parts techniques, to familiar energy equations as well as the div curl = 0 and curl grad = 0 identities. For the diffusion equation, these operators describe compact schemes whose convergence is assured by the energy equations and which yield both the potential and the flux vector with second-order accuracy. A simplified potential form is especially useful for obtaining numerical results by multigrid and ADI methods.

Rose, Milton E.↗
NASA NTRS1989

A pseudo-spectral collocation method applied to the problem of convective diffusive transport in fluids subject to unsteady residual accelerations

The problem of determining the sensitivity of Bridgman-Stockbarger directional solidification experiments to residual accelerations of the type associated with spacecraft in low earth orbit is analyzed numerically using a pseudo-spectral collocation method. The approach employs a novel iterative scheme combining the method of artificial compressibility and a generalized ADI method. The results emphasize the importance of the consideration of residual accelerations and careful selection of the operating conditions in order to take full advantages of the low gravity conditions.

Alexander, J. Iwan↗
NASA NTRS1991

Navier-Stokes computation of wing/rotor interaction for a tilt rotor in hover

A method has been developed to analyze the wing/rotor interaction of tilt rotor aircraft in hover. The unsteady, thin-layer compressible Navier-Stokes equations are solved using an implicit, finite difference scheme that employs LU-ADI factorization. The rotor is modeled as an actuator disk which imparts a radial and azimuthal distribution of pressure rise and swirl to the flowfield. The 'chimera' approach of grid point blanking is used to update the rotor boundary conditions. Results are presented for both a rotor alone and for wing/rotor interaction where the thrust coefficient is 0.0164 and wing flap deflection is 67 degrees. Many of the complex flow features are captured including the fountain effect, leading and trailing edge separation, and the unsteady wake beneath the wing. Wing surface pressures compare fairly well with experimental data although the time-averaged download is about twenty percent higher than the measured value. This discrepancy is due to a combination of factors that are discussed.

Fejtek, Ian↗
NASA NTRS1991

Computational aspects of chemically reacting flows

An implicit, time-dependent approach is taken to solve the inviscid conservative equations in a vertex-based finite-volume formulation. The convective fluxes are estimated by a modified Riemann solver, and higher-order accuracy is obtained by means of the MUSCL. Implicit techniques are adopted to deal with the equation stiffness as a result of the clustered grid and chemical source term. The real-gas Roe matrix and the ADI and LU schemes have been tested on the sample problems: hypersonic flow over blunt configurations, such as the forebody of a double ellipse, and the near wake of a conic aerobrake. Temperature distribution and standoff distance are validated against the shock-fitting, central-differencing results. Equilibrium temperature may be recovered by artificially increasing the chemical rates at the wall. The present method satisfactorily resolves the flow structure, including the strong expansion behind the flat base. The relative merits and performance of the real-gas TVD and implicit schemes are discussed briefly.

Li, C. P.↗
NASA NTRS1991

A numerical study of the flow establishment time in hypersonic shock tunnels

This paper presents a numerical study of unsteady hypersonic viscous flows in 2D hypersonic shock tunnels. The numerical experiments on the combined nozzle/model flows, utilizing the configurations mounted in the test section of a 14.7-long shock tunnel, which has a design Mach number of Me = 6 and a stagnation temperature of T(0) = 3063 R, have been performed by solving the 2D Navier-Stokes equations with an ADI-type TVD scheme. The flow establishment times around aerodynamic models, such as a flat plate, a biconvex airfoil, thin and thick double wedges, and a circular cylinder, are calculated and compared. The analysis shows that, for the most severely separated case of a circular cylinder, the flow becomes stabilized within about 8 msec including the starting process, and its nondimensional flow establishment time is 46. The flow establishment time increases with increasing length of the separated region and decreases with increasing freestream velocity.

Lee, Jang-Yeon↗
NASA NTRS1991

Effects of leading-edge flap oscillation on unsteady delta wing flow and rock control

The isolated and interdisciplinary problems of unsteady fluid dynamics and rigid-body dynamics and control of delta wings with and without leading-edge flap oscillation are considered. For the fluid dynamics problem, the unsteady, compressible, thin-layer Navier-Stokes (NS) equations, which are written relative to a moving frame of reference, are solved along with the unsteady, linearized, Navier-displacement (ND) equations. The NS equations are solved for the flowfield using an implicit finite-volume scheme. The ND equations are solved for the grid deformation, if the leading-edge flaps oscillate, using an ADI scheme. For the dynamics and control problem, the Euler equation of rigid-body rolling motion for a wing and its flaps are solved interactively with the fluid dynamics equations for the wing-rock motion and subsequently for its control. A four-stage Runge-Kutta scheme is used to explicitly integrate the dynamics equation.

Kandil, Osama A.↗
NASA NTRS1991

A Fourier-Chebyshev pseudospectral method for solving steady 3-D Navier-Stokes and heat equations in cylindrical cavities

A Fourier-Chebyshev pseudospectral method for solving steady 3D Navier-Stokes equations in cylindrical cavities is presented and discussed. The general method is pseudo-unsteady and uses a semi-implicit finite difference scheme for the time discretization. The generalized ADI (Alternating Direction Implicit) procedure is then applied to reduce the problem to a successive solution of one-dimensional problems. The spatial approximation uses a Fourier-Galerkin approximation in the periodic direction and a Chebyshev-collocation approximation in the other directions. Difficulties related to the pressure are surmounted by using the artifical compressibility method. A suitable variable change was chosen to avoid the problem of singularity at the axis generated by cylindrical coordinates. The method is first tested on an advection-diffusion equation and then on the Navier-Stokes equations. Finally, the method is illustrated by a convection problem of a differentially heated fluid.

Pulicani, J. P.↗
NASA NTRS1992

Implementation/validation of a low Reynolds number two-equation turbulence model in the Proteus Navier-Stokes code: Two-dimensional/axisymmetric

The implementation and validation of the Chien low Reynolds number k-epsilon turbulence model in the two dimensional axisymmetric version Proteus, a compressible Navier-Stokes computer code, are presented. The set of k-epsilon equations are solved by marching in time using a coupled alternating direction implicit (ADI) solution procedure with generalized first or second order time differencing. To validate Proteus and the k-epsilon turbulence model, laminar and turbulent computations were done for several benchmark test cases: incompressible fully developed 2-D channel flow; fully developed axisymmetric pipe flow; boundary layer flow over a flat plate; and turbulent Sajben subsonic transonic diffuser flows. Proteus results from these test cases showed good agreement with analytical results and experimental data. Detailed comparisons of both mean flow and turbulent quantities showed that the Chien k-epsilon turbulence model given good results over a wider range of turbulent flow than the Baldwin-Lomax turbulence model in the Proteus code with no significant CPU time penalty for more complicated flow cases.

Bui, Trong T.↗
NASA NTRS1992

Development of a multigrid transonic potential flow code for cascades

Finite-volume methods for discretizing transonic potential flow equations have proven to be very flexible and accurate for both two and three dimensional problems. Since they only use local properties of the mapping, they allow decoupling of the grid generation from the rest of the problem. A very effective method for solving the discretized equations and converging to a solution is the multigrid-ADI technique. It has been successfully applied to airfoil problems where O type, C type and slit mappings have been used. Convergence rates for these cases are more than an order of magnitude faster than with relaxation techniques. In this report, we describe a method to extend the above methods, with the C type mappings, to airfoil cascade problems.

Steinhoff, John↗
NASA NTRS1992

Fast methods to numerically integrate the Reynolds equation for gas fluid films

The alternating direction implicit (ADI) method is adopted, modified, and applied to the Reynolds equation for thin, gas fluid films. An efficient code is developed to predict both the steady-state and dynamic performance of an aerodynamic journal bearing. An alternative approach is shown for hybrid journal gas bearings by using Liebmann's iterative solution (LIS) for elliptic partial differential equations. The results are compared with known design criteria from experimental data. The developed methods show good accuracy and very short computer running time in comparison with methods based on an inverting of a matrix. The computer codes need a small amount of memory and can be run on either personal computers or on mainframe systems.

Dimofte, Florin↗
NASA NTRS1992

The impact of time step definition on code convergence and robustness

We have implemented preconditioning for multi-species reacting flows in two independent codes, an implicit (ADI) code developed in-house and the RPLUS code (developed at LeRC). The RPLUS code was modified to work on a four-stage Runge-Kutta scheme. The performance of both the codes was tested, and it was shown that preconditioning can improve convergence by a factor of two to a hundred depending on the problem. Our efforts are currently focused on evaluating the effect of chemical sources and on assessing how preconditioning may be applied to improve convergence and robustness in the calculation of reacting flows.

Venkateswaran, S.↗