Engineering topics

Kwak, Dochan

Publications and source records attributed to Kwak, Dochan.

At least 73 records · Page 4

NASA NTRS2002

Time Dependent Density Functional Theory Calculations of Large Compact PAH Cations: Implications for the Diffuse Interstellar Bands

We investigate the electronic absorption spectra of several maximally pericondensed polycyclic aromatic hydrocarbon radical cations with time dependent density functional theory calculations. We find interesting trends in the vertical excitation energies and oscillator strengths for this series containing pyrene through circumcoronene, the largest species containing more than 50 carbon atoms. We discuss the implications of these new results for the size and structure distribution of the diffuse interstellar band carriers.

Weisman, Jennifer L.↗

NASA NTRS2002

Electronic Absorption Spectra of Neutral Perylene (C20H12), Terrylene (C30H16), and Quaterrylene (C40H20) and their Positive and Negative Ions: Ne Matrix-Isolation Spectroscopy and Time Dependent Density Functional Theory Calculations

We present a full experimental and theoretical study of an interesting series of polycyclic aromatic hydrocarbons, the oligorylenes. The absorption spectra of perylene, terrylene and quaterrylene in neutral, cationic and anionic charge states are obtained by matrix-isolation spectroscopy in Ne. The experimental spectra are dominated by a bright state that red shifts with growing molecular size. Excitation energies and state symmetry assignments are supported by calculations using time dependent density functional theory methods. These calculations also provide new insight into the observed trends in oscillator strength and excitation energy for the bright states: the oscillator strength per unit mass of carbon increases along the series.

Halasinski, Thomas M.↗

NASA NTRS2002

Definition of Contravariant Velocity Components

This is an old issue in computational fluid dynamics (CFD). What is the so-called contravariant velocity or contravariant velocity component? In the article, we review the basics of tensor analysis and give the contravariant velocity component a rigorous explanation. For a given coordinate system, there exist two uniquely determined sets of base vector systems - one is the covariant and another is the contravariant base vector system. The two base vector systems are reciprocal. The so-called contravariant velocity component is really the contravariant component of a velocity vector for a time-independent coordinate system, or the contravariant component of a relative velocity between fluid and coordinates, for a time-dependent coordinate system. The contravariant velocity components are not physical quantities of the velocity vector. Their magnitudes, dimensions, and associated directions are controlled by their corresponding covariant base vectors. Several 2-D (two-dimensional) linear examples and 2-D mass-conservation equation are used to illustrate the details of expressing a vector with respect to the covariant and contravariant base vector systems, respectively.

Hung, Ching-Mao↗

NASA NTRS2002

Quantum Dynamics Study of the Isotopic Effect on Capture Reactions: HD, D2 + CH3

Time-dependent wave-packet-propagation calculations are reported for the isotopic reactions, HD + CH3 and D2 + CH3, in six degrees of freedom and for zero total angular momentum. Initial state selected reaction probabilities for different initial rotational-vibrational states are presented in this study. This study shows that excitations of the HD(D2) enhances the reactivities; whereas the excitations of the CH3 umbrella mode have the opposite effects. This is consistent with the reaction of H2 + CH3. The comparison of these three isotopic reactions also shows the isotopic effects in the initial-state-selected reaction probabilities. The cumulative reaction probabilities (CRP) are obtained by summing over initial-state-selected reaction probabilities. The energy-shift approximation to account for the contribution of degrees of freedom missing in the six dimensionality calculation is employed to obtain approximate full-dimensional CRPs. The rate constant comparison shows H2 + CH3 reaction has the biggest reactivity, then HD + CH3, and D2 + CH3 has the smallest.

Wang, Dunyou↗

NASA NTRS2002

The Overgrid Interface for Computational Simulations on Overset Grids

Computational simulations using overset grids typically involve multiple steps and a variety of software modules. A graphical interface called OVERGRID has been specially designed for such purposes. Data required and created by the different steps include geometry, grids, domain connectivity information and flow solver input parameters. The interface provides a unified environment for the visualization, processing, generation and diagnosis of such data. General modules are available for the manipulation of structured grids and unstructured surface triangulations. Modules more specific for the overset approach include surface curve generators, hyperbolic and algebraic surface grid generators, a hyperbolic volume grid generator, Cartesian box grid generators, and domain connectivity: pre-processing tools. An interface provides automatic selection and viewing of flow solver boundary conditions, and various other flow solver inputs. For problems involving multiple components in relative motion, a module is available to build the component/grid relationships and to prescribe and animate the dynamics of the different components.

Chan, William M.↗

NASA NTRS2002

Automated CFD Parameter Studies on Distributed Parallel Computers

The objective of the current work is to build a prototype software system which will automated the process of running CFD jobs on Information Power Grid (IPG) resources. This system should remove the need for user monitoring and intervention of every single CFD job. It should enable the use of many different computers to populate a massive run matrix in the shortest time possible. Such a software system has been developed, and is known as the AeroDB script system. The approach taken for the development of AeroDB was to build several discrete modules. These include a database, a job-launcher module, a run-manager module to monitor each individual job, and a web-based user portal for monitoring of the progress of the parameter study. The details of the design of AeroDB are presented in the following section. The following section provides the results of a parameter study which was performed using AeroDB for the analysis of a reusable launch vehicle (RLV). The paper concludes with a section on the lessons learned in this effort, and ideas for future work in this area.

Rogers, Stuart E.↗

NASA NTRS2002

Unsteady Computations of a Jet in a Crossflow with Ground Effect

A numerical study of a jet in crossflow with ground effect is conducted using OVERFLOW with dual time-stepping and low Mach number preconditioning. The results of the numerical study are compared to an experiment to show that the numerical methods are capable of capturing the dominant features of the flow field as well as the unsteadiness associated with the ground vortex.

Pandya, Shishir↗

NASA NTRS2002

Quantal Study of the Exchange Reaction for N + N2 using an ab initio Potential Energy Surface

The N + N2 exchange rate is calculated using a time-dependent quantum dynamics method on a newly determined ab initio potential energy surface (PES) for the ground A" state. This ab initio PES shows a double barrier feature in the interaction region with the barrier height at 47.2 kcal/mol, and a shallow well between these two barriers, with the minimum at 43.7 kcal/mol. A quantum dynamics wave packet calculation has been carried out using the fitted PES to compute the cumulative reaction probability for the exchange reaction of N + N2(J=O). The J - K shift method is then employed to obtain the rate constant for this reaction. The calculated rate constant is compared with experimental data and a recent quasi-classical calculation using a LEPS PES. Significant differences are found between the present and quasiclassical results. The present rate calculation is the first accurate 3D quantal dynamics study for N + N2 reaction system and the ab initio PES reported here is the first such surface for N3.

Wang, Dunyou↗

NASA NTRS2002

Detection of HCN and C2H2 in ISO Spectra of Oxygen-Rich AGB Stars

Cool oxygen-rich AGB stars were not expected to have organic molecules like HCN in either their photospheres or circumstellar envelopes (CSEs). The discovery of HCN and CS microwave emission from the shallowest CSE layers of these stars was a considerable surprise and much theoretical effort has been expended in explaining the presence of such organics. To further explore this problem, we have undertaken a systematic search of oxygen-rich AGB stellar spectra in the Infrared Space Observatory (ISO) data archive. Our purposes are to find evidence regarding critical molecular species that could be of value in choosing among the proposed theoretical models, to locate spectral features which might give clues to conditions deeper in the CSEs, and to lay the groundwork for future SIRTF (Space Infrared Telescope Facility) and SOFIA (Stratospheric Observatory for Infrared Astronomy) observations. Using carefully reduced observations, we have detected weak absorption features arising from HCN and possibly C2H2 in a small number of oxygen-rich AGB stars. The most compelling case is NML Cyg which shows both HCN (14 microns) and CO2 (15 microns). VY CMa, a similar star, shows evidence for HCN, but not CO2. Two S-type stars show evidence for the C-H bending transitions: W Aql at 14 microns (HCN) and both W Aql and S Cas at 13.7 microns (C2H2). Both W Aql and S Cas as well as S Lyr, a SC-type star, show 3 micron absorption which may arise from the C-H stretch of HCN and C2H2. In the case of NML Cyg, we show that the HCN and CO2 spectral features are formed in the CSE at temperatures well above those of the outermost CSE layers and derive approximate column densities. In the case of the S-stars, we discuss the evidence for the organic features and their photospheric origin.

Carbon, Duane F.↗

NASA NTRS2002

Dissociative Ionization of Pyridine by Electron Impact

In order to understand the damage of biomolecules by electrons, a process important in radiation damage, we undertake a study of the dissociative ionization (DI) of pyridine (C5H5N) from the low-lying ionization channels. The methodology used is the same as in the benzene study. While no experimental DI data are available, we compare the dissociation products from our calculations with the dissociative photoionization measurements of Tixier et al. using dipole (e, e(+) ion) coincidence spectroscopy. Comparisons with the DI of benzene is also made so as to understand the difference in DI between a heterocyclic and an aromatic molecule.

Dateo, Christopher↗

NASA NTRS2002

Progress Toward Generation of a Navier-Stokes Database for a Harrier in Ground Effect

The Harrier YAV-8B aircraft is capable of vertical and short-field take-off and landing (V/STOL) by directing its four exhaust nozzles toward the ground, or conventional flight by rotating its nozzles into a horizontal position. The British Royal Air Force and the United States Marine Corps have used this aircraft for more than 30 years to provide a quick reaction time for troop support, and reduce the need for long runways. The success of this powered-lift (PL) vehicle has also prompted the more recent design of the Joint Strike Fighter (JSF). However there are significant safety issues that must be addressed when operating a PL vehicle in close proximity to the ground. Hot Gas Ingestion (HGI) by the inlets can result in a rapid loss of powered lift; and high-speed jet flows along the ground plane can induce low pressures underneath the vehicle, causing a 'suck-down' effect. Under these conditions, departure from controlled flight may occur. Moreover, unsteady ground vortices and jet fountains can affect the aircraft,s controllability and its proximity to ground troops. The viscous, time-dependent flow fields of PL vehicles are difficult to accurately and efficiently predict using Computational Fluid Dynamics (CFD). A number of researchers have used the time-dependent Reynolds-averaged Navier-Stokes (RANS) equations to compute flows for single and multiple jets in a cross-flow. A few have added some geometric complexity to the problem by computing flows for jet-augmented delta wings near a ground plane. Smith et.al. computed for the first time a single RANS solution about a simplified Harrier. This geometry included a fuselage, wing, leading edge root extension (LERX), inlets, and exhaust nozzles. All of these investigations cite two practical problems with computing these flows: 1) the need for improved solution accuracy; and, 2) the need for faster solution methods. We view the need for faster solution methods as key to improving the solution accuracy and making this class of computation more routine. One can hardly refine grids, explore the use of advanced turbulence models, and generate databases when it takes weeks of dedicated computer time for a single solution. Chaderjian, Ahmad, Pandya, and Murman have focused on reducing the time-to-solution for this very difficult and complex problem through process automation and exploitation of parallel computing. They began with the Harrier geometry reported, and added a deflected wing flap and empennage for greater realism. To date more than 80 solutions have been carried out. This paper will describe this process and progress made in reducing the time required to generate a simple longitudinal force and moment database for a Harrier in ground effect. It shows a typical snap-shot from an unsteady streakline animation, where fluid particles are colored by temperature. The ground vortex and a jet-fountain vortex are highlighted. It also shows a similar streakline image, where HGI occurs due to the vehicle in close proximity to the ground. It is show the mean lift coefficient as a function of angle of attack and height. The angle of attack range was 4 deg less than or = alpha less than or = 10 deg with an increment of 1 degree, and the height range was 10 ft less than or = h less than or = 30ft with an increment of 5 feet. This 35 solution database was extended to over 2500 cases using a monotone cubic-spline interpolation procedure. The suck-down effect (reduction of lift near the ground) is highlighted in the figure. The "cushion effect," the conventional reduction of lift as the vehicle moves out of ground effect, is also indicated. All 35 RANS solutions were obtained using 952 Silicon Graphics Origin 2000 and 3000 processors in dedicated mode for one week. Typically, 112 processors were assigned to each case. Some other cases used fewer processors to utilize all available CPUS. The final paper will report on the automation of the solution process, including: grid generation, job monitoring, solution completion criteria, and post processing. Moreover, improvements in parallel efficiency for a dual time-step algorithm for the RANS equations will also be presented. Results will be discussed in detail using unsteady streakline flow visualization to correlate unsteady flow structures with dominant aerodynamic frequencies. The stability derivatives, CL, and CL, will also be presented.

Chaderjian, Neal M.↗

NASA NTRS2002

Best Practices In Overset Grid Generation

Grid generation for overset grids on complex geometry can be divided into four main steps: geometry processing, surface grid generation, volume grid generation and domain connectivity. For each of these steps, the procedures currently practiced by experienced users are described. Typical problems encountered are also highlighted and discussed. Most of the guidelines are derived from experience on a variety of problems including space launch and return vehicles, subsonic transports with propulsion and high lift devices, supersonic vehicles, rotorcraft vehicles, and turbomachinery.

Chan, William M.↗

NASA NTRS2002

Progress in Unsteady Turbopump Flow Simulations Using Overset Grid Systems

This viewgraph presentation provides information on unsteady flow simulations for the Second Generation RLV (Reusable Launch Vehicle) baseline turbopump. Three impeller rotations were simulated by using a 34.3 million grid points model. MPI/OpenMP hybrid parallelism and MLP shared memory parallelism has been implemented and benchmarked in INS3D, an incompressible Navier-Stokes solver. For RLV turbopump simulations a speed up of more than 30 times has been obtained. Moving boundary capability is obtained by using the DCF module. Scripting capability from CAD geometry to solution is developed. Unsteady flow simulations for advanced consortium impeller/diffuser by using a 39 million grid points model are currently underway. 1.2 impeller rotations are completed. The fluid/structure coupling is initiated.

Kiris, Cetin C.↗

NASA NTRS2002

Comparison of the T(sub 1) and D(sub 1) Diagnostics: A New Definition for the Open-Shell D(sub 1) Diagnostic

It is shown that the T(sub 1) operator used in a previous study to define the open-shell D(sub 1) diagnostic is invalid, and leads to an arbitrary definition of the open-shell D(sub 1) diagnostic. A new definition is proposed that eliminates this ambiguity and approximately restores the mathematical relationship previously noted between the closed-shell D(sub 1) and T(sub 1) diagnostics. Statistical comparison of the T(sub 1) and D(sub 1) diagnostics shows a very high degree of correlation between them, although it is argued that both diagnostics used together can provide more information than either can separately.

Lee, Timothy J.↗

NASA NTRS2001

Total Electron-Impact Ionization Cross-Sections of CFx and NFx (x = 1 - 3)

The discrepancy between experimental and theoretical total electron-impact ionization cross sections for a group of fluorides, CFx, and NFx, (x = 1 - 3), is attributed to the inadequacies in previous theoretical models. Cross-sections calculated using a recently developed siBED (simulation Binary-Encounter-Dipole) model that takes into account the shielding of the long-range dipole potential between the scattering electron and target are in agreement with experimentation. The present study also carefully reanalyzed the previously reported experimental data to account for the possibility of incomplete collection of fragment ions and the presence of ion-pair formation channels. For NF3, our experimental and theoretical cross-sections compare well with the total ionization cross-sections recently reported by Haaland et al. in the region below dication formation.

Huo, Winifred M.↗

NASA NTRS2001

Towards Accurate Ab Initio Predictions of the Spectrum of Methane

We have carried out extensive ab initio calculations of the electronic structure of methane, and these results are used to compute vibrational energy levels. We include basis set extrapolations, core-valence correlation, relativistic effects, and Born- Oppenheimer breakdown terms in our calculations. Our ab initio predictions of the lowest lying levels are superb.

Schwenke, David W.↗

NASA NTRS2001

Lifetime of a Chemically Bound Helium Compound

The rare-gas atoms are chemically inert, to an extent unique among all elements. This is due to the stable electronic structure of the atoms. Stable molecules with chemically bound rare-gas atoms are, however, known. A first such compound, XePtF6, W2S prepared in 1962 and since then a range of molecules containing radon, xenon and krypton have been obtained. Most recently, a first stable chemically bound compound of argon was prepared, leaving neon and helium as the only elements for which stable chemically bound molecules are not yet known. Electronic structure calculations predict that a metastable species HHeF exists, but significance of the result depends on the unknown lifetime. Here we report quantum dynamics calculations of the lifetime of HHeF, using accurate interactions computed from electronic structure theory. HHeF is shown to disintegrate by tunneling through energy barriers into He + HF and H + He + F the first channel greatly dominating. The lifetime of HHeF is more than 120 picoseconds, that of DHeF is 14 nanoseconds. The relatively long lifetimes are encouraging for the preparation prospects of this first chemically bound helium compound.

Chaban, Galina M.↗

NASA NTRS2001

Time-Accurate Solutions of Incompressible Navier-Stokes Equations for Potential Turbopump Applications

Two numerical procedures, one based on artificial compressibility method and the other pressure projection method, are outlined for obtaining time-accurate solutions of the incompressible Navier-Stokes equations. The performance of the two method are compared by obtaining unsteady solutions for the evolution of twin vortices behind a at plate. Calculated results are compared with experimental and other numerical results. For an un- steady ow which requires small physical time step, pressure projection method was found to be computationally efficient since it does not require any subiterations procedure. It was observed that the artificial compressibility method requires a fast convergence scheme at each physical time step in order to satisfy incompressibility condition. This was obtained by using a GMRES-ILU(0) solver in our computations. When a line-relaxation scheme was used, the time accuracy was degraded and time-accurate computations became very expensive.

Kiris, Cetin↗