Engineering topics

Kwak, Dochan

Publications and source records attributed to Kwak, Dochan.

At least 55 records · Page 3

NASA NTRS2002

Time-Dependent Simulations of Turbopump Flows

Unsteady flow simulations for RLV (Reusable Launch Vehicles) 2nd Generation baseline turbopump for one and half impeller rotations have been completed by using a 34.3 Million grid points model. MLP (Multi-Level Parallelism) shared memory parallelism has been implemented in INS3D, and benchmarked. Code optimization for cash based platforms will be completed by the end of September 2001. Moving boundary capability is obtained by using DCF module. Scripting capability from CAD (computer aided design) geometry to solution has been developed. Data compression is applied to reduce data size in post processing. Fluid/Structure coupling has been initiated.

Kiris, Cetin↗

NASA NTRS2002

Numerical Simulation of Rolling-Airframes Using a Multi-Level Cartesian Method

A supersonic rolling missile with two synchronous canard control surfaces is analyzed using an automated, inviscid, Cartesian method. Sequential-static and time-dependent dynamic simulations of the complete motion are computed for canard dither schedules for level flight, pitch, and yaw maneuver. The dynamic simulations are compared directly against both high-resolution viscous simulations and relevant experimental data, and are also utilized to compute dynamic stability derivatives. The results show that both the body roll rate and canard dither motion influence the roll-averaged forces and moments on the body. At the relatively, low roll rates analyzed in the current work these dynamic effects are modest, however the dynamic computations are effective in predicting the dynamic stability derivatives which can be significant for highly-maneuverable missiles.

Murman, Scott M.↗

NASA NTRS2002

Spectral (Finite) Volume Method for Conservation Laws on Unstructured Grids II: Extension to Two Dimensional Scalar Equation

The framework for constructing a high-order, conservative Spectral (Finite) Volume (SV) method is presented for two-dimensional scalar hyperbolic conservation laws on unstructured triangular grids. Each triangular grid cell forms a spectral volume (SV), and the SV is further subdivided into polygonal control volumes (CVs) to supported high-order data reconstructions. Cell-averaged solutions from these CVs are used to reconstruct a high order polynomial approximation in the SV. Each CV is then updated independently with a Godunov-type finite volume method and a high-order Runge-Kutta time integration scheme. A universal reconstruction is obtained by partitioning all SVs in a geometrically similar manner. The convergence of the SV method is shown to depend on how a SV is partitioned. A criterion based on the Lebesgue constant has been developed and used successfully to determine the quality of various partitions. Symmetric, stable, and convergent linear, quadratic, and cubic SVs have been obtained, and many different types of partitions have been evaluated. The SV method is tested for both linear and non-linear model problems with and without discontinuities.

Wang, Z. J.↗

NASA NTRS2002

Multilevel Error Estimation and Adaptive h-Refinement for Cartesian Meshes with Embedded Boundaries

This paper presents the development of a mesh adaptation module for a multilevel Cartesian solver. While the module allows mesh refinement to be driven by a variety of different refinement parameters, a central feature in its design is the incorporation of a multilevel error estimator based upon direct estimates of the local truncation error using tau-extrapolation. This error indicator exploits the fact that in regions of uniform Cartesian mesh, the spatial operator is exactly the same on the fine and coarse grids, and local truncation error estimates can be constructed by evaluating the residual on the coarse grid of the restricted solution from the fine grid. A new strategy for adaptive h-refinement is also developed to prevent errors in smooth regions of the flow from being masked by shocks and other discontinuous features. For certain classes of error histograms, this strategy is optimal for achieving equidistribution of the refinement parameters on hierarchical meshes, and therefore ensures grid converged solutions will be achieved for appropriately chosen refinement parameters. The robustness and accuracy of the adaptation module is demonstrated using both simple model problems and complex three dimensional examples using meshes with from 10(exp 6), to 10(exp 7) cells.

Aftosmis, M. J.↗

NASA NTRS2002

Progress in the Simulation of Steady and Time-Dependent Flows with 3D Parallel Unstructured Cartesian Methods

The proposed paper will present recent extensions in the development of an efficient Euler solver for adaptively-refined Cartesian meshes with embedded boundaries. The paper will focus on extensions of the basic method to include solution adaptation, time-dependent flow simulation, and arbitrary rigid domain motion. The parallel multilevel method makes use of on-the-fly parallel domain decomposition to achieve extremely good scalability on large numbers of processors, and is coupled with an automatic coarse mesh generation algorithm for efficient processing by a multigrid smoother. Numerical results are presented demonstrating parallel speed-ups of up to 435 on 512 processors. Solution-based adaptation may be keyed off truncation error estimates using tau-extrapolation or a variety of feature detection based refinement parameters. The multigrid method is extended to for time-dependent flows through the use of a dual-time approach. The extension to rigid domain motion uses an Arbitrary Lagrangian-Eulerlarian (ALE) formulation, and results will be presented for a variety of two- and three-dimensional example problems with both simple and complex geometry.

Aftosmis, M. J.↗

NASA NTRS2002

Theoretical Study of Tautomerization Reactions for the Ground and First Excited Electronic States of Adenine

Geometrical structures and energetic properties for different tautomers of adenine are calculated in this study, using multi-configurational wave functions. Both the ground and the lowest singlet excited state potential energy surfaces are studied. Four tautomeric forms are considered, and their energetic order is found to be different on the ground and the excited state potential energy surfaces. Minimum energy reaction paths are obtained for hydrogen atom transfer (tautomerization) reactions in the ground and the lowest excited electronic states. It is found that the barrier heights and the shapes of the reaction paths are different for the ground and the excited electronic states, suggesting that the probability of such tautomerization reaction is higher on the excited state potential energy surface. This tautomerization process should become possible in the presence of water or other polar solvent molecules and should play an important role in the photochemistry of adenine.

Salter, Latasha M.↗

NASA NTRS2002

Spectral (Finite) Volume Method for One Dimensional Euler Equations

Consider a mesh of unstructured triangular cells. Each cell is called a Spectral Volume (SV), denoted by Si, which is further partitioned into subcells named Control Volumes (CVs), indicated by C(sub i,j). To represent the solution as a polynomial of degree m in two dimensions (2D) we need N = (m+1)(m+2)/2 pieces of independent information, or degrees of freedom (DOFs). The DOFs in a SV method are the volume-averaged mean variables at the N CVs. For example, to build a quadratic reconstruction in 2D, we need at least (2+1)(3+1)/2 = 6 DOFs. There are numerous ways of partitioning a SV, and not every partition is admissible in the sense that the partition may not be capable of producing a degree m polynomial. Once N mean solutions in the CVs of a SV are given, a unique polynomial reconstruction can be obtained.

Wang, Z. J.↗

NASA NTRS2002

An Accurate Quartic Force Field, Fundamental Frequencies, and Binding Energy for the High Energy Density Material T(d)N4

The CCSD(T) method has been used to compute a highly accurate quartic force field and fundamental frequencies for all N-14 and N-15 isotopomers of the high energy density material T(sub d)N(sub 4). The computed fundamental frequencies show beyond doubt that the bands observed in a matrix isolation experiment by Radziszewski and coworkers are not due to different isotopomers of T(sub d)N(sub 4). The most sophisticated thermochemical calculations to date yield a N(sub 4) -> 2N(sub 2) heat of reaction of 182.22 +/- 0.5 kcal/mol at 0 K (180.64 +/- 0.5 at 298 K). It is hoped that the data reported herein will aid in the ultimate detection of T(sub d)N(sub 4).

Lee, Timothy J.↗

NASA NTRS2002

Blood Pump Development Using Rocket Engine Flow Simulation Technology

This viewgraph presentation provides information on the transfer of rocket engine flow simulation technology to work involving the development of blood pumps. Details are offered regarding the design and requirements of mechanical heart assist devices, or VADs (ventricular assist device). There are various computational fluid dynamics issues involved in the visualization of flow in such devices, and these are highlighted and compared to those of rocket turbopumps.

Kiris, Cetin C.↗

NASA NTRS2002

Theoretical Study of Infrared and Raman Spectra of Hydrated Magnesium Sulfate Salts

Harmonic and anharmonic vibrational frequencies, as well as infrared and Raman intensities, are calculated for MgSO4.nH20 (n=1-3). Electronic structure theory at the second order Moller-Plesset perturbation theory (MP2) level with a triple-zeta + polarization (TZP) basis set is used to determine the geometry, properties, and vibrational spectra of pure and hydrated MgSO4 salts. The direct vibrational self-consistent field (VSCF) method and its correlation corrected (CC-VSCF) extension are used to determine anharmonic corrections to vibrational frequencies and intensities for the pure MgSO4 and its complex with one water molecule. Very significant differences are found between vibrational of water molecules in complexes with MgSO4 and pure water. Some of the O-H stretching frequencies are shifted to the red very significantly (by up to 1500-2000/cm) upon complexation with magnesium sulfate. They should be observed between 1700 and 3000/cm in a region very different from the corresponding O-H stretch frequency region of pure water (3700-3800/cm). In addition, the SO2 stretching vibrations are found at lower frequency regions than the water vibrations. They can serve as unique identifiers for the presence of sulfate salts. The predicted infrared and Raman spectra should be of valuable help in the design of future missions and analysis of observed data from the ice surface of Jupiter's moon Europa that possibly contains hydrated MgSO4 salts.

Chaban, Galina M.↗

NASA NTRS2002

Parametric Time-Dependent Navier-Stokes Computations for a YAV-8B Harrier in Ground Effect

The Harrier Jump Jet has the distinction of being the only powered-lift aircraft in the free world to achieve operational status and to have flown in combat. This V/STOL aircraft can take-off and land vertically or utilize very short runways by directing its four exhaust nozzles towards the ground. Transition to forward flight is achieved by rotating these nozzles into a horizontal position. Powered-lift vehicles have certain advantages over conventional strike fighters. Their V/STOL capabilities allow for safer carrier operations, smaller carrier size, and quick reaction time for troop support. Moreover, they are not dependent on vulnerable land-based runways. The AV-8A Harrier first entered service in the British Royal Air Force (RAF) during 1969, and the U.S. Marine Corps (USMC) in 1971. The AV-8B was a redesign to achieve improved payload capacity, range, and accuracy. This modified design first entered service with the USMC and RAF in 1985. The success and unique capabilities of the Harrier has prompted the design of a powered-lift version of the Joint Strike Fighter (JSF). The flowfield for the Harrier near the ground during low-speed or hover flight operations is very complex and time-dependent. A sketch of this flowfield is shown. Warm air from the fan is exhausted from the front nozzles, while a hot air/fuel mixture from the engine is exhausted from the rear nozzles. These jets strike the ground and move out radially forming a ground jet-flow. The ambient freestream, due to low-speed forward flight or - headwind during hover, opposes the jet-flow. This interaction causes the flow to separate and form a ground vortex. The multiple jets also interact with each other near the ground and form an upwash or jet fountain, which strikes the underside of the fuselage. If the aircraft is sufficiently close to the ground, the inlet can ingest ground debris and hot gases from the fountain and ground vortex. This Hot Gas Ingestion (HGI) can cause a sudden loss of thrust (powered lift), and the vehicle may crash. The high-speed jet flow along the ground can also entrain the ambient flow, resulting in a low pressure region underneath the vehicle. The accompanied loss of lift is referred to as the suckdown effect. The ground vortex may also be highly unsteady, dramatically changing its size and position with time at low frequencies, e.g., 1 Hz.

Chaderjian, Neal M.↗

NASA NTRS2002

Definition of Contravariant Velocity Components

In this paper we have reviewed the basics of tensor analysis in an attempt to clarify some misconceptions regarding contravariant and covariant vector components as used in fluid dynamics. We have indicated that contravariant components are components of a given vector expressed as a unique combination of the covariant base vector system and, vice versa, that the covariant components are components of a vector expressed with the contravariant base vector system. Mathematically, expressing a vector with a combination of base vector is a decomposition process for a specific base vector system. Hence, the contravariant velocity components are decomposed components of velocity vector along the directions of coordinate lines, with respect to the covariant base vector system. However, the contravariant (and covariant) components are not physical quantities. Their magnitudes and dimensions are controlled by their corresponding covariant (and contravariant) base vectors.

Hung, Ching-moa↗

NASA NTRS2002

Use of High Fidelity Methods in Multidisciplinary Optimization-A Preliminary Survey

Multidisciplinary optimization is a key element of design process. To date multidiscipline optimization methods that use low fidelity methods are well advanced. Optimization methods based on simple linear aerodynamic equations and plate structural equations have been applied to complex aerospace configurations. However, use of high fidelity methods such as the Euler/ Navier-Stokes for fluids and 3-D (three dimensional) finite elements for structures has begun recently. As an activity of Multidiscipline Design Optimization Technical Committee (MDO TC) of AIAA (American Institute of Aeronautics and Astronautics), an effort was initiated to assess the status of the use of high fidelity methods in multidisciplinary optimization. Contributions were solicited through the members MDO TC committee. This paper provides a summary of that survey.

Guruswamy, Guru P.↗

NASA NTRS2002

PEGASUS 5: An Automated Pre-Processor for Overset-Grid CFD

An all new, automated version of the PEGASUS software has been developed and tested. PEGASUS provides the hole-cutting and connectivity information between overlapping grids, and is used as the final part of the grid generation process for overset-grid computational fluid dynamics approaches. The new PEGASUS code (Version 5) has many new features: automated hole cutting; a projection scheme for fixing gaps in overset surfaces; more efficient interpolation search methods using an alternating digital tree; hole-size optimization based on adding additional layers of fringe points; and an automatic restart capability. The new code has also been parallelized using the Message Passing Interface standard. The parallelization performance provides efficient speed-up of the execution time by an order of magnitude, and up to a factor of 30 for very large problems. The results of three example cases are presented: a three-element high-lift airfoil, a generic business jet configuration, and a complete Boeing 777-200 aircraft in a high-lift landing configuration. Comparisons of the computed flow fields for the airfoil and 777 test cases between the old and new versions of the PEGASUS codes show excellent agreement with each other and with experimental results.

Suhs, Norman E.↗

NASA NTRS2002

Parametric Study of a YAV-8B Harrier in Ground Effect using Time-Dependent Navier-Stokes Computations

A process is described which enables the generation of 35 time-dependent viscous solutions for a YAV-8B Harrier in ground effect in one week. Overset grids are used to model the complex geometry of the Harrier aircraft and the interaction of its jets with the ground plane and low-speed ambient flow. The time required to complete this parametric study is drastically reduced through the use of process automation, modern computational platforms, and parallel computing. Moreover, a dual-time-stepping algorithm is described which improves solution robustness. Unsteady flow visualization and a frequency domain analysis are also used to identify and correlated key flow structures with the time variation of lift.

Pandya, Shishir↗

NASA NTRS2002

Progress in Unsteady Turbopump Flow Simulations

This viewgraph presentation discusses unsteady flow simulations for a turbopump intended for a reusable launch vehicle (RLV). The simulation process makes use of computational grids and parallel processing. The architecture of the parallel computers used is discussed, as is the scripting of turbopump simulations.

Kiris, Cetin C.↗

NASA NTRS2002

Development and Applications of a Modular Parallel Process for Large Scale Fluid/Structures Problems

A modular process that can efficiently solve large scale multidisciplinary problems using massively parallel supercomputers is presented. The process integrates disciplines with diverse physical characteristics by retaining the efficiency of individual disciplines. Computational domain independence of individual disciplines is maintained using a meta programming approach. The process integrates disciplines without affecting the combined performance. Results are demonstrated for large scale aerospace problems on several supercomputers. The super scalability and portability of the approach is demonstrated on several parallel computers.

Guruswamy, Guru P.↗

NASA NTRS2002

Prebiotic Processing induced by Comet and Meteor Impact

In their study of organic synthesis from impact shocks using the laser-induced-plasma (LIP) technique, McKay and Borucki(l) found that organic synthesis preferentially occurred in a reducing gas mixture rich in methane, and not in a mixture rich in carbon dioxide. This result means chemical models based on the thermodynamical equilibrium approach do not apply to shock chemistry. In this study, we employ the technique of reacting flow, i.e., chemical kinetics in a fluid flow, to simulate the chemistry occurring in LIP and in the wake region from comet or meteor impact. Three different air compositions have been used: (1) 1/3 CO2 and 2/3 H2, (2) pure CH4, and (3) 1/4 CH4, 1/4 CO2, and 1/2 H2O. The stoichiometric ratio of gas mixtures (1) and (3) are kept the same. For (1) we obtain equal mole fractions of CO and H2O as the major products and for (2) C2H2 is the major product. In both cases our results are in agreement with Ref. (1). For (3) we find an interesting case where the nature of chemicals produced to be critically dependent on the flowfield temperature. At the higher temperature part of the wake region, CO and H2O are the dominant products, whereas in the cooler region C2H2 is the dominant product. Further studies of these reactions, as well as for the gas mixture including N2, are being pursued.

Dateo, Christopher E.↗