An upwind-differencing scheme for the incompressible Navier-Stokes equations

The steady state incompressible Navier-Stokes equations in 2-D are solved numerically using the artificial compressibility formulation. The convective terms are upwind-differenced using a flux difference split approach that has uniformly high accuracy throughout the interior grid points. The viscous fluxes are differenced using second order accurate central differences. The numerical system of equations is solved using an implicit line relaxation scheme. Although the current study is limited to steady state problems, it is shown that this entire formulation can be used for solving unsteady problems. Characteristic boundary conditions are formulated and used in the solution procedure. The overall scheme is capable of being run at extremely large pseudotime steps, leading to fast convergence. Three test cases are presented to demonstrate the accuracy and robustness of the code. These are the flow in a square-driven cavity, flow over a backward facing step, and flow around a 2-D circular cylinder

Computation of incompressible viscous flows through turbopump components

A finite-difference, three-dimensional, incompressible Navier-Stokes formulation for calculating the flow through turbopump components is presented. The solution method is based on the pseudocompressibility approach and uses an implicit-upwind differencing scheme together with the Gauss-Seidel line-relaxation method. Both steady and unsteady flow calculations can be performed using the presented algorithm. In this paper, the equations are solved in steadily rotating reference frames by using the steady-state formulation in order to simulate the flow through a turbopump inducer. Eddy viscosity is computed by using the Baldwin-Lomax model. Numerical results are compared with experimental measurements and good agreement is found between the two. Time-accurate calculations will be reported in future publications

Numerical Simulations of a Nonlinear Transmission Line

This thesis investigates wave propagation along a nonlinear transmission line having a voltage-dependent capacitance. The telegrapher\u27s equations that model wave propagation along such a transmission line are derived and shown to represent a nonlinear hyperbolic system of balance laws. The Lax-Friedrichs, Lax-Wendroff, and hybrid numerical schemes for obtaining approximate solutions to nonlinear hyperbolic systems of balance laws are presented, analyzed, and applied to the nonlinear telegrapher\u27s equations. The Lax-Wendroff and hybrid schemes are invoked to numerically simulate wave propagation along a nonlinear transmission line. Simulations obtained via the hybrid scheme are used to briefly study the potential application of a nonlinear transmission line as a pulse generator. Since they may generate ultra short (nanosecond) rise time voltage waveforms, nonlinear transmission lines could have applicability to pulsed power systems and high power microwaves. Moreover, short duration, high intensity voltage pulses have been demonstrated to penetrate biological cells and affect the internal organelles. Thus, such voltage shaping could provide the technological tools for intracellular manipulation and several potential biomedical engineering applications. Though only one particular voltage-dependent capacitance was selected for the nonlinear transmission line analysis, the method described is general. Other voltage dependencies could be realized through suitable mate1ial tailoring, by using composites, and by sectioning the transmission line layout. These alterations would afford relative assessments of voltage wave propagation in nonlinear transmission lines, and also provide for analyses of shock wave generation

Efficient simulation of incompressible viscous flow over multi-element airfoils

The incompressible, viscous, turbulent flow over single and multi-element airfoils is numerically simulated in an efficient manner by solving the incompressible Navier-Stokes equations. The computer code uses the method of pseudo-compressibility with an upwind-differencing scheme for the convective fluxes and an implicit line-relaxation solution algorithm. The motivation for this work includes interest in studying the high-lift take-off and landing configurations of various aircraft. In particular, accurate computation of lift and drag at various angles of attack, up to stall, is desired. Two different turbulence models are tested in computing the flow over an NACA 4412 airfoil; an accurate prediction of stall is obtained. The approach used for multi-element airfoils involves the use of multiple zones of structured grids fitted to each element. Two different approaches are compared: a patched system of grids, and an overlaid Chimera system of grids. Computational results are presented for two-element, three-element, and four-element airfoil configurations. Excellent agreement with experimental surface pressure coefficients is seen. The code converges in less than 200 iterations, requiring on the order of one minute of CPU time (on a CRAY YMP) per element in the airfoil configuration

Numerical Simulation of Flow Through an Artificial Heart

A solution procedure was developed that solves the unsteady, incompressible Navier-Stokes equations, and was used to numerically simulate viscous incompressible flow through a model of the Pennsylvania State artificial heart. The solution algorithm is based on the artificial compressibility method, and uses flux-difference splitting to upwind the convective terms; a line-relaxation scheme is used to solve the equations. The time-accuracy of the method is obtained by iteratively solving the equations at each physical time step. The artificial heart geometry involves a piston-type action with a moving solid wall. A single H-grid is fit inside the heart chamber. The grid is continuously compressed and expanded with a constant number of grid points to accommodate the moving piston. The computational domain ends at the valve openings where nonreflective boundary conditions based on the method of characteristics are applied. Although a number of simplifing assumptions were made regarding the geometry, the computational results agreed reasonably well with an experimental picture. The computer time requirements for this flow simulation, however, are quite extensive. Computational study of this type of geometry would benefit greatly from improvements in computer hardware speed and algorithm efficiency enhancements

Computation of incompressible viscous flows through artificial heart devices with moving boundaries

The extension of computational fluid dynamics techniques to artificial heart flow simulations is illustrated. Unsteady incompressible Navier-Stokes equations written in 3-D generalized curvilinear coordinates are solved iteratively at each physical time step until the incompressibility condition is satisfied. The solution method is based on the pseudo compressibility approach and uses an implicit upwind differencing scheme together with the Gauss-Seidel line relaxation method. The efficiency and robustness of the time accurate formulation of the algorithm are tested by computing the flow through model geometries. A channel flow with a moving indentation is computed and validated with experimental measurements and other numerical solutions. In order to handle the geometric complexity and the moving boundary problems, a zonal method and an overlapping grid embedding scheme are used, respectively. Steady state solutions for the flow through a tilting disk heart valve was compared against experimental measurements. Good agreement was obtained. The flow computation during the valve opening and closing is carried out to illustrate the moving boundary capability

Tracking Debris Shed by a Space-Shuttle Launch Vehicle

The DEBRIS software predicts the trajectories of debris particles shed by a space-shuttle launch vehicle during ascent, to aid in assessing potential harm to the space-shuttle orbiter and crew. The user specifies the location of release and other initial conditions for a debris particle. DEBRIS tracks the particle within an overset grid system by means of a computational fluid dynamics (CFD) simulation of the local flow field and a ballistic simulation that takes account of the mass of the particle and its aerodynamic properties in the flow field. The computed particle trajectory is stored in a file to be post-processed by other software for viewing and analyzing the trajectory. DEBRIS supplants a prior debris tracking code that took .15 minutes to calculate a single particle trajectory: DEBRIS can calculate 1,000 trajectories in .20 seconds on a desktop computer. Other improvements over the prior code include adaptive time-stepping to ensure accuracy, forcing at least one step per grid cell to ensure resolution of all CFD-resolved flow features, ability to simulate rebound of debris from surfaces, extensive error checking, a builtin suite of test cases, and dynamic allocation of memory

Numerical solution of the incompressible Navier-Stokes equations

The current work is initiated in an effort to obtain an efficient, accurate, and robust algorithm for the numerical solution of the incompressible Navier-Stokes equations in two- and three-dimensional generalized curvilinear coordinates for both steady-state and time-dependent flow problems. This is accomplished with the use of the method of artificial compressibility and a high-order flux-difference splitting technique for the differencing of the convective terms. Time accuracy is obtained in the numerical solutions by subiterating the equations in psuedo-time for each physical time step. The system of equations is solved with a line-relaxation scheme which allows the use of very large pseudo-time steps leading to fast convergence for steady-state problems as well as for the subiterations of time-dependent problems. Numerous laminar test flow problems are computed and presented with a comparison against analytically known solutions or experimental results. These include the flow in a driven cavity, the flow over a backward-facing step, the steady and unsteady flow over a circular cylinder, flow over an oscillating plate, flow through a one-dimensional inviscid channel with oscillating back pressure, the steady-state flow through a square duct with a 90 degree bend, and the flow through an artificial heart configuration with moving boundaries. An adequate comparison with the analytical or experimental results is obtained in all cases. Numerical comparisons of the upwind differencing with central differencing plus artificial dissipation indicates that the upwind differencing provides a much more robust algorithm, which requires significantly less computing time. The time-dependent problems require on the order of 10 to 20 subiterations, indicating that the elliptical nature of the problem does require a substantial amount of computing effort