Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solve time”

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 577 records · Page 32

Dynamic grid deformation using Navier-displacement equation for deforming wings

For dynamic and aeroelastic applications of maneuvering wings, the solid boundaries undergo rigid-body motion and aeroelastic deformation. For rigid-body motion, the conservative fluid dynamics equation, in terms of the Eulerian description, is written in terms of a moving frame of reference, and the problem is solved on a time-independent body-conformed grid. For both rigid-body motion and aeroelastic deformation, the Navier-displacement equation, in terms of the Lagrangian coordinates, is modified for fluid-flow problems. It is used along with the Eulerian description of the conservative fluid dynamics equations to account for the grid deformation.

Kandil, Osama A.↗

Penetrative convective flows induced by internal heating and mantle compressibility

Penetrative convective flows induced in a spherical shell by combined effects of internal heating and mantle compressibility are investigated using mathematical and numerical formulations for compressible spherical shell convection. Isothermal stress-free boundary conditions applied at the top and the bottom of the shell are solved using a time-dependent finite difference code in a temperature, vorticity, stream function formulation for Rayleigh numbers ranging from the critical Rc up to 2000 Rc. Results indicate that compressibility, together with internal heating, could be a mechanism capable of generating spontaneously layered convection and local melting in the mantle and that non-Boussinesq effects must be considered in interpretations of geophysical phenomena.

Machetel, Philippe↗

Analysis of opposed jet hydrogen-air counter flow diffusion flame

A computational simulation of the opposed-jet diffusion flame is performed to study its structure and extinction limits. The present analysis concentrates on the nitrogen-diluted hydrogen-air diffusion flame, which provides the basic information for many vehicle designs such as the aerospace plane for which hydrogen is a candidate as the fuel. The computer program uses the time-marching technique to solve the energy and species equations coupled with the momentum equation solved by the collocation method. The procedure is implemented in two stages. In the first stage, a one-step forward overal chemical reaction is chosen with the gas phase chemical reaction rate determined by comparison with experimental data. In the second stage, a complete chemical reaction mechanism is introduced with detailed thermodynamic and transport property calculations. Comparison between experimental extinction data and theoretical predictions is discussed. The effects of thermal diffusion as well as Lewis number and Prandtl number variations on the diffusion flame are also presented.

Ho, Y. H.↗

Evolution of midplate hotspot swells: Numerical solutions

The evolution of midplate hotspot swells on an oceanic plate moving over a hot, upwelling mantle plume is numerically simulated. The plume supplies a Gaussian-shaped thermal perturbation and thermally-induced dynamic support. The lithosphere is treated as a thermal boundary layer with a strongly temperature-dependent viscosity. The two fundamental mechanisms of transferring heat, conduction and convection, during the interaction of the lithosphere with the mantle plume are considered. The transient heat transfer equations, with boundary conditions varying in both time and space, are solved in cylindrical coordinates using the finite difference ADI (alternating direction implicit) method on a 100 x 100 grid. The topography, geoid anomaly, and heat flow anomaly of the Hawaiian swell and the Bermuda rise are used to constrain the models. Results confirm the conclusion of previous works that the Hawaiian swell can not be explained by conductive heating alone, even if extremely high thermal perturbation is allowed. On the other hand, the model of convective thinning predicts successfully the topography, geoid anomaly, and the heat flow anomaly around the Hawaiian islands, as well as the changes in the topography and anomalous heat flow along the Hawaiian volcanic chain.

Liu, Mian↗

PROTEUS two-dimensional Navier-Stokes computer code, version 1.0. Volume 1: Analysis description

A new computer code was developed to solve the two-dimensional or axisymmetric, Reynolds averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The thin-layer or Euler equations may also be solved. Turbulence is modeled using an algebraic eddy viscosity model. The objective was to develop a code for aerospace applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The equations are written in nonorthogonal body-fitted coordinates, and solved by marching in time using a fully-coupled alternating direction-implicit procedure with generalized first- or second-order time differencing. All terms are linearized using second-order Taylor series. The boundary conditions are treated implicitly, and may be steady, unsteady, or spatially periodic. Simple Cartesian or polar grids may be generated internally by the program. More complex geometries require an externally generated computational coordinate system. The documentation is divided into three volumes. Volume 1 is the Analysis Description, and describes in detail the governing equations, the turbulence model, the linearization of the equations and boundary conditions, the time and space differencing formulas, the ADI solution procedure, and the artificial viscosity models.

Towne, Charles E.↗

PROTEUS two-dimensional Navier-Stokes computer code, version 1.0. Volume 2: User's guide

A new computer code was developed to solve the two-dimensional or axisymmetric, Reynolds averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The thin-layer or Euler equations may also be solved. Turbulence is modeled using an algebraic eddy viscosity model. The objective was to develop a code for aerospace applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The equations are written in nonorthogonal body-fitted coordinates, and solved by marching in time using a fully-coupled alternating direction-implicit procedure with generalized first- or second-order time differencing. All terms are linearized using second-order Taylor series. The boundary conditions are treated implicitly, and may be steady, unsteady, or spatially periodic. Simple Cartesian or polar grids may be generated internally by the program. More complex geometries require an externally generated computational coordinate system. The documentation is divided into three volumes. Volume 2 is the User's Guide, and describes the program's general features, the input and output, the procedure for setting up initial conditions, the computer resource requirements, the diagnostic messages that may be generated, the job control language used to run the program, and several test cases.

Towne, Charles E.↗

PROTEUS two-dimensional Navier-Stokes computer code, version 1.0. Volume 3: Programmer's reference

A new computer code was developed to solve the 2-D or axisymmetric, Reynolds-averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The thin-layer or Euler equations may also be solved. Turbulence is modeled using an algebraic eddy viscosity model. The objective was to develop a code for aerospace applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The equations are written in nonorthogonal body-fitted coordinates, and solved by marching in time using a fully-coupled alternating-direction-implicit procedure with generalized first- or second-order time differencing. All terms are linearized using second-order Taylor series. The boundary conditions are treated implicitly, and may be steady, unsteady, or spatially periodic. Simple Cartesian or polar grids may be generated internally by the program. More complex geometries require an externally generated computational coordinate system. The documentation is divided into three volumes. Volume 3 is the Programmer's Reference, and describes the program structure, the FORTRAN variables stored in common blocks, and the details of each subprogram.

Towne, Charles E.↗

Three-dimensional calculation of supersonic reacting flows using an LU scheme

A new three-dimensional numerical program incorporated with comprehensive real gas property models was developed to simulate supersonic reacting flows. The code employs an implicit finite volume, Lower-Upper (LU) time-marching method to solve the complete Navier-Stokes and species equations in a fully-coupled and very efficient manner. A chemistry model with nine species and eighteen reaction steps are adopted in the program to represent the chemical reaction of H2 and air. To demonstrate the capability of the program, flow fields of underexpanded hydrogen jets transversely injected into supersonic air stream inside the combustors of scramjets are calculated. Results clearly depict the flow characteristics, including the shock structure, separated flow regions around the injector, and the distribution of the combustion products.

Yu, Sheng-Tao↗

A general multiblock Euler code for propulsion integration. Volume 1: Theory document

A general multiblock Euler solver was developed for the analysis of flow fields over geometrically complex configurations either in free air or in a wind tunnel. In this approach, the external space around a complex configuration was divided into a number of topologically simple blocks, so that surface-fitted grids and an efficient flow solution algorithm could be easily applied in each block. The computational grid in each block is generated using a combination of algebraic and elliptic methods. A grid generation/flow solver interface program was developed to facilitate the establishment of block-to-block relations and the boundary conditions for each block. The flow solver utilizes a finite volume formulation and an explicit time stepping scheme to solve the Euler equations. A multiblock version of the multigrid method was developed to accelerate the convergence of the calculations. The generality of the method was demonstrated through the analysis of two complex configurations at various flow conditions. Results were compared to available test data. Two accompanying volumes, user manuals for the preparation of multi-block grids (vol. 2) and for the Euler flow solver (vol. 3), provide information on input data format and program execution.

Chen, H. C.↗

A general multiblock Euler code for propulsion integration. Volume 3: User guide for the Euler code

This manual explains the procedures for using the general multiblock Euler (GMBE) code developed under NASA contract NAS1-18703. The code was developed for the aerodynamic analysis of geometrically complex configurations in either free air or wind tunnel environments (vol. 1). The complete flow field is divided into a number of topologically simple blocks within each of which surface fitted grids and efficient flow solution algorithms can easily be constructed. The multiblock field grid is generated with the BCON procedure described in volume 2. The GMBE utilizes a finite volume formulation with an explicit time stepping scheme to solve the Euler equations. A multiblock version of the multigrid method was developed to accelerate the convergence of the calculations. This user guide provides information on the GMBE code, including input data preparations with sample input files and a sample Unix script for program execution in the UNICOS environment.

Chen, H. C.↗

Second mode interactions in supersonic boundary layers

The nonlinear evolution of a two-dimensional second mode unstable wave in a Mach 4.5 wall-bounded flow is computed by solving the full time-dependent compressible Navier-Stokes equations. A highly accurate solution is obtained using spectral collocation methods. It is shown that departure from linearity first occurs in the critical layer due to the cubic nonlinearities in the momentum equation. This is a direct result of the large density perturbations in this regime. Time evolution studies of the growth rate as a function of normal distance from the plate suggests that the mode is evolving toward a nonlinear saturated state, and that this problem is possibly amenable to standard weakly nonlinear perturbation methods.

Erlebacher, Gordon↗

Io's volcanic and sublimation atmospheres

Fully 3D axisymmetric gasdynamic equations simulating SO2 and H2S frost sublimation and SO2 dayside and nightside volcanic atmospheres on Io are numerically solved, using a time-explicit finite-volume formulation. Both the sublimation and volcanic atmospheres generate horizontal supersonic winds away from the subsolar point or the volcanic vent. While the sublimation atmosphere is primarily driven by horizontal pressure gradients determined by surface temperatures, the volcanic atmosphere is driven by pressure gradients that are determined by the source rate. Sublimation and condensation produce patterns of surface deposits which are characteristic of the two types of atmospheres. The volcanic model is quantitatively consistent with Voyager observations of ring deposits.

Moreno, Miguel A.↗

Autonomous space processor for orbital debris

The development of an Autonomous Space Processor for Orbital Debris (ASPOD) was the goal. The nature of this craft, which will process, in situ, orbital debris using resources available in low Earth orbit (LEO) is explained. The serious problem of orbital debris is briefly described and the nature of the large debris population is outlined. The focus was on the development of a versatile robotic manipulator to augment an existing robotic arm, the incorporation of remote operation of the robotic arms, and the formulation of optimal (time and energy) trajectory planning algorithms for coordinated robotic arms. The mechanical design of the new arm is described in detail. The work envelope is explained showing the flexibility of the new design. Several telemetry communication systems are described which will enable the remote operation of the robotic arms. The trajectory planning algorithms are fully developed for both the time optimal and energy optimal problems. The time optimal problem is solved using phase plane techniques while the energy optimal problem is solved using dynamic programming.

Ramohalli, Kumar↗

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.↗

Three-dimensional calculations of supersonic reacting flows using an LU scheme

A 3-D numerical program that incorporates comprehensive real gas property models was developed to simulate supersonic reacting flows. The code employs an implicit, finite volume, Lower-Upper (LU), time-marching method to solve the complete Navier-Stokes and species equations in a fully-coupled and efficient manner. A chemistry model with 9 species and 18 reaction steps is adopted in the program to represent the chemical reactions of H2 and air. To demonstrate the capability of the program, flow fields of underexpanded hydrogen jets transversely injected into the supersonic airstream inside the combustors of scramjets are calculated. Results clearly depict the flow characteristics, including the shock structure, the separated flow regions around the injector, and the distribution of the combustion products.

Yu, Sheng-Tao↗

Prediction of unsteady rotor-surface pressure and heat transfer from wake passings

The research described in this paper is a numerical investigation of the effects of unsteady flow on gas turbine heat transfer, particularly on a rotor blade surface. The unsteady flow in a rotor blade passage and the unsteady heat transfer on the blade surface as a result of wake/blade interaction are modeled by the inviscid flow/boundary layer approach. The Euler equations which govern the inviscid flow are solved using a time-accurate marching scheme. The unsteady flow in the blade passage is induced by periodically moving a wake model across the passage inlet. Unsteady flow solutions in the passage provide pressure gradients and boundary conditions for the boundary-layer equations which govern the viscous flow adjacent to the blade surface. Numerical solutions of the unsteady turbulent boundary layer yield surface heat flow values which can then be compared to experimental data. Comparisons with experimental data show that unsteady heat flux on the blade suction surface is well predicted, but the predictions of unsteady heat flux on the blade pressure surface do not agree.

Tran, Le T.↗

On waveform multigrid method

Waveform multigrid method is an efficient method for solving certain classes of time dependent PDEs. This paper studies the relationship between this method and the analogous multigrid method for steady-state problems. Using a Fourier-Laplace analysis, practical convergence rate estimates of the waveform multigrid iterations are obtained. Experimental results show that the analysis yields accurate performance prediction.

Taasan, Shlomo↗