Engineering topics

Chung, T. J.

Publications and source records attributed to Chung, T. J..

At least 19 records

NASA NTRS2003

Finite Element Method for Capturing Ultra-relativistic Shocks

While finite element methods are used extensively by researchers solving computational fluid dynamics in fields other than astrophysics, their use in astrophysical fluid simulations has been predominantly overlooked. Current simulations using other methods such as finite difference and finite volume (based on finite difference) have shown remarkable results, but these methods are limited by their fundamental properties in aspects that are important for simulations with complex geometries and widely varying spatial and temporal scale differences. We have explored the use of finite element methods for astrophysical fluids in order to establish the validity of using such methods in astrophysical environments. We present our numerical technique applied to solving ultra-relativistic (Lorentz Factor Gamma >> 1) shocks which are prevalent in astrophysical studies including relativistic jets and gamma-ray burst studies. We show our finite element formulation applied to simulations where the Lorentz factor ranges up to 2236 and demonstrate its stability in solving ultra-relativistic flows. Our numerical method is based on the Flowfield Dependent Variation (FDV) Method, unique in that numerical diffusion is derived from physical parameters rather than traditional artificial viscosity methods. Numerical instabilities account for most of the difficulties when capturing shocks in this regime. Our method results in stable solutions and accurate results as compared with other methods.

Richardson, G. A.↗

NASA NTRS2002

Computational Relativistic Astrophysics Using the Flow Field-Dependent Variation Theory

We present our method for solving general relativistic nonideal hydrodynamics. Relativistic effects become pronounced in such cases as jet formation from black hole magnetized accretion disks which may lead to the study of gamma-ray bursts. Nonideal flows are present where radiation, magnetic forces, viscosities, and turbulence play an important role. Our concern in this paper is to reexamine existing numerical simulation tools as to the accuracy and efficiency of computations and introduce a new approach known as the flow field-dependent variation (FDV) method. The main feature of the FDV method consists of accommodating discontinuities of shock waves and high gradients of flow variables such as occur in turbulence and unstable motions. In this paper, the physics involved in the solution of relativistic hydrodynamics and solution strategies of the FDV theory are elaborated. The general relativistic astrophysical flow and shock solver (GRAFSS) is introduced, and some simple example problems for computational relativistic astrophysics (CRA) are demonstrated.

Richardson, G. A.↗

NASA NTRS2001

Parallelization of the Flow Field Dependent Variation Scheme for Solving the Triple Shock/Boundary Layer Interaction Problem

A parallelized version of the Flowfield Dependent Variation (FDV) Method is developed to analyze a problem of current research interest, the flowfield resulting from a triple shock/boundary layer interaction. Such flowfields are often encountered in the inlets of high speed air-breathing vehicles including the NASA Hyper-X research vehicle. In order to resolve the complex shock structure and to provide adequate resolution for boundary layer computations of the convective heat transfer from surfaces inside the inlet, models containing over 500,000 nodes are needed. Efficient parallelization of the computation is essential to achieving results in a timely manner. Results from a parallelization scheme, based upon multi-threading, as implemented on multiple processor supercomputers and workstations is presented.

Schunk, Richard Gregory↗

NASA NTRS2001

Computational Relativistic Astrophysics Using the Flowfield-Dependent Variation Theory

Theoretical models, observations and measurements have preoccupied astrophysicists for many centuries. Only in recent years, has the theory of relativity as applied to astrophysical flows met the challenges of how the governing equations can be solved numerically with accuracy and efficiency. Even without the effects of relativity, the physics of magnetohydrodynamic flow instability, turbulence, radiation, and enhanced transport in accretion disks has not been completely resolved. Relativistic effects become pronounced in such cases as jet formation from black hole magnetized accretion disks and also in the study of Gamma-Ray bursts (GRB). Thus, our concern in this paper is to reexamine existing numerical simulation tools as to the accuracy and efficiency of computations and introduce a new approach known as the flowfield-dependent variation (FDV) method. The main feature of the FDV method consists of accommodating discontinuities of shock waves and high gradients of flow variables such as occur in turbulence and unstable motions. In this paper, the physics involved in the solution of relativistic hydrodynamics and solution strategies of the FDV theory are elaborated. The general relativistic astrophysical flow and shock solver (GRAFSS) is introduced, and some simple example problems for Computational Relativistic Astrophysics (CRA) are demonstrated.

Richardson, G. A.↗

NASA NTRS2000

Parallel Processing in Combustion Analysis

The objective of this research is to demonstrate the application of the Flow-field Dependent Variation (FDV) method to a problem of current interest in supersonic chemical combustion. Due in part to the stiffness of the chemical reactions, the solution of such problems on unstructured three dimensional grids often dictates the use of parallel computers. Preliminary results for the injection of a supersonic hydrogen stream into vitiated air are presented.

Schunk, Richard Gregory↗

NASA NTRS2000

Airbreathing Propulsion System Analysis Using Multithreaded Parallel Processing

In this paper, parallel processing is used to analyze the mixing, and combustion behavior of hypersonic flow. Preliminary work for a sonic transverse hydrogen jet injected from a slot into a Mach 4 airstream in a two-dimensional duct combustor has been completed [Moon and Chung, 1996]. Our aim is to extend this work to three-dimensional domain using multithreaded domain decomposition parallel processing based on the flowfield-dependent variation theory. Numerical simulations of chemically reacting flows are difficult because of the strong interactions between the turbulent hydrodynamic and chemical processes. The algorithm must provide an accurate representation of the flowfield, since unphysical flowfield calculations will lead to the faulty loss or creation of species mass fraction, or even premature ignition, which in turn alters the flowfield information. Another difficulty arises from the disparity in time scales between the flowfield and chemical reactions, which may require the use of finite rate chemistry. The situations are more complex when there is a disparity in length scales involved in turbulence. In order to cope with these complicated physical phenomena, it is our plan to utilize the flowfield-dependent variation theory mentioned above, facilitated by large eddy simulation. Undoubtedly, the proposed computation requires the most sophisticated computational strategies. The multithreaded domain decomposition parallel processing will be necessary in order to reduce both computational time and storage. Without special treatments involved in computer engineering, our attempt to analyze the airbreathing combustion appears to be difficult, if not impossible.

Schunk, Richard Gregory↗

NASA NTRS1999

Unified CFD Methods Via Flowfield-Dependent Variation Theory

This paper addresses the flowfield-dependent variation (FDV) methods in which complex physical phenomena are taken into account in the final form of partial differential equations to be solved so that finite difference methods (FDM) or finite element methods (FEM) themselves will not dictate the physics, but rather are no more than simply the options how to discretize between adjacent nodal points or within an element. The variation parameters introduced in the formulation are calculated from the current flowfield based on changes of Mach numbers, Reynolds numbers, Peclet numbers, and Damkohler numbers between adjacent nodal points, which play many significant roles such as adjusting the governing equations (hyperbolic, parabolic, and/or e!liptic), resolving various physical phenomena, and controlling the accuracy and stability of the numerical solution. The theory is verified by a number of example problems addressing the physical implications of the variation parameters which resemble the flowfield itself, shock capturing mechanism, transitions and interactions between inviscid/viscous, compressibility/incompressibility, and laminar/turbulent flows.

Chung, T. J.↗

NASA NTRS1997

Flowfield-Dependent Mixed Explicit-Implicit (FDMEL) Algorithm for Computational Fluid Dynamics

Despite significant achievements in computational fluid dynamics, there still remain many fluid flow phenomena not well understood. For example, the prediction of temperature distributions is inaccurate when temperature gradients are high, particularly in shock wave turbulent boundary layer interactions close to the wall. Complexities of fluid flow phenomena include transition to turbulence, relaminarization separated flows, transition between viscous and inviscid incompressible and compressible flows, among others, in all speed regimes. The purpose of this paper is to introduce a new approach, called the Flowfield-Dependent Mixed Explicit-Implicit (FDMEI) method, in an attempt to resolve these difficult issues in Computational Fluid Dynamics (CFD). In this process, a total of six implicitness parameters characteristic of the current flowfield are introduced. They are calculated from the current flowfield or changes of Mach numbers, Reynolds numbers, Peclet numbers, and Damkoehler numbers (if reacting) at each nodal point and time step. This implies that every nodal point or element is provided with different or unique numerical scheme according to their current flowfield situations, whether compressible, incompressible, viscous, inviscid, laminar, turbulent, reacting, or nonreacting. In this procedure, discontinuities or fluctuations of an variables between adjacent nodal points are determined accurately. If these implicitness parameters are fixed to certain numbers instead of being calculated from the flowfield information, then practically all currently available schemes of finite differences or finite elements arise as special cases. Some benchmark problems to be presented in this paper will show the validity, accuracy, and efficiency of the proposed methodology.

Garcia, S. M.↗

NASA NTRS1995

Influence of Two-Phase Thermocapillary Flow on Liquid Retention in Microscopic Pores

An important feature of screened propellant acquisition devices is the retention capability or maximum maintainable pressure difference across the porous barrier separating the liquid and gas. Previous experiments with liquid hydrogen showed a marked reduction in retention when the tank containing the device was pressurized with hydrogen vapor. These tests, however did not indicate any appreciable degradation in retention with helium pressurization or direct heating through the screen. The objective of this article is to determine if the thermocapillary convection arising from phase change in the microscopic pores of such screens could cause these disparities in performance. A numerical model of flow in a single pore suggests that the thermocapillary-induced gradient in liquid pressure along the surface can strongly affect surface morphology. In an evaporative environment, this gradient exerts a stabilizing influence on surface curvature, and preserves the momentum balance between the liquid and gas. With condensation, it causes a force imbalance and a destabilizing suction in the middle of the pore that reduces retention. Results also indicate that introducing an inert gas, such as helium, suppresses this retention loss mechanism by lowering thermocapillary circulation and its associated interfacial pressure gradient.

Schmidt, George R.↗

NASA NTRS1995

Thermocapillary flow with evaporation and condensation at low gravity. Part 2: Deformable surface

The free surface behavior of a volatile wetting liquid at low gravity is studied using scaling and numerical techniques. An open cavity model, which was applied in part 1 to investigate fluid flow and heat transfer in non-deforming pores, is used to evaluate the influence of convection on surface morphology with length scales and subcooling/superheating limits of 1 less than or equal to D less than or equal to 10(exp 2) microns and approximately 1 K, respectively. Results show that the menisci shapes of highly wetting fluids are sensitive to thermocapillary flow and to a lesser extent the recoil force associated with evaporation and condensation. With subcooling, thermocapillarity produces a suction about the pore centerline that promotes loss of mechanical equilibrium, while condensation exerts an opposing force that under some conditions offsets this destabilizing influence. With superheating, thermocapillarity and evaporation act in the same direction and mutually foster surface stability. All of these trends are magnified by high capillary and Biot numbers, and the stronger circulation intensities associated with small contact angles. These phenomena strongly depend on the thermal and interfacial equilibrium between the liquid and vapor, and have important ramifications for systems designed to maintain a pressure differential across a porous surface.

Schmidt, G. R.↗

NASA NTRS1994

Influence of Two-Phase Thermocapillary Flow on Cryogenic Liquid Retention in Microscopic Pores

Previous experiments indicate that the bubble point pressure of spacecraft liquid hydrogen acquisition devices is reduced substantially when the ullage is pressurized with heated hydrogen vapor. The objective is to determine whether the two-phase thermocapillary convection arising from thermodynamic non-equilibrium along the porous surfaces of such devices could lead to this observed degradation in retention performance. We also examine why retention capability appears to be unaffected by pressurization with heated helium or direct heating through the porous structure. Computational assessments based on coupled solution of the flowfield and liquid free surface indicate that for highly wetting fluids in small pores, dynamic pressure and vapor recoil dictate surface morphology and drive meniscus deformation. With superheating, the two terms exert the same influence on curvature and promote mechanical equilibrium, but with subcooling, the pressure distribution produces a suction about the pore center-line that degrades retention. This result points to thermocapillary-induced deformation arising from condensation as the cause for retention loss. It also indicates that increasing the level of non-equilibrium by reducing accommodation coefficient restricts deformation and explains why retention failure does not occur with direct screen heating or helium pressurization.

Schmidt, G. R.↗

NASA NTRS1992

Thermocapillary flow about an evaporating meniscus

The steady motion and thermal behavior of an evaporating superheated liquid in a small cavity bounded by isothermal sidewalls is examined. Scaling analyses and a two-dimensional finite element model are used to investigate the influence of thermocapillarity, buoyancy, and temperature-dependent mass flux on flowfield, interfacial heat transfer, and meniscus morphology. Numerical investigations indicate the existence of two counter-rotating cells symmetric about the cavity center. Results also show that evaporation tends to counteract this circulation by directing flow toward the hotter sidewalls. Although thermocapillarity and evaporation yield different flowfield distributions, both effects tend to increase interfacial temperature and heat transfer.

Schmidt, G. R.↗

NASA NTRS1991

Combustion wave instabilities in turbulent flows

Navier-Stokes solutions are used to calculate oscillatory components of pressure, velocity, and density, which in turn provide necessary data to compute energy growth factors to determine combustion instability. It is shown that wave instabilities are associated with changes in entropy and the space and time averages of oscillatory components of pressure, velocity and density, together with the mean flow field in the energy equation. Compressible laminar and turbulent flows and reacting flows with hydrogen/oxygen combustion are considered in this paper. The space shuttle main engine combustion/thrust chamber is used for illustration of the theory. The analysis shows that the increase of mean pressure and disturbances consistently results in the increase in instability. It is shown that adequate combustion instability analysis requires at least third order nonlinearity in energy growth or decay.

Yoon, W. S.↗

NASA NTRS1990

CFD code evaluation for internal flow modeling

Research on the computational fluid dynamics (CFD) code evaluation with emphasis on supercomputing in reacting flows is discussed. Advantages of unstructured grids, multigrids, adaptive methods, improved flow solvers, vector processing, parallel processing, and reduction of memory requirements are discussed. As examples, researchers include applications of supercomputing to reacting flow Navier-Stokes equations including shock waves and turbulence and combustion instability problems associated with solid and liquid propellants. Evaluation of codes developed by other organizations are not included. Instead, the basic criteria for accuracy and efficiency have been established, and some applications on rocket combustion have been made. Research toward an ultimate goal, the most accurate and efficient CFD code, is in progress and will continue for years to come.

Chung, T. J.↗

NASA NTRS1990

Navier-Stokes Entropy Controlled Combustion Instability Analysis for Liquid Propellants

Navier-Stokes solutions are used to calculate oscillatory components of pressure, velocity, and density, which in turn provide necessary data to compute energy growth factors to determine combustion instability. It is shown that wave instabilities are associated with changes in entropy and the space and time averages of oscillatory components of pressure, velocity and density, together with the mean flow field in the energy equation. Compressible laminar and turbulent flows and reacting flows with hydrogen/oxygen combustion are considered. The SSME combustion/thrust chamber is used for illustration of the theory. The analysis shows that the increase of mean pressure and disturbances consistently results in the increase of instability. It is shown that adequate combustion instability analysis requires at least third order nonlinearity in energy growth or decay.

Chung, T. J.↗

NASA NTRS1990

Combustion instability analysis

A theory and computer program for combustion instability analysis are presented. The basic theoretical foundation resides in the concept of entropy-controlled energy growth or decay. Third order perturbation expansion is performed on the entropy-controlled acoustic energy equation to obtain the first order integrodifferential equation for the energy growth factor in terms of the linear, second, and third order energy growth parameters. These parameters are calculated from Navier-Stokes solutions with time averages performed on as many Navier-Stokes time steps as required to cover at least one peak wave period. Applications are made for a 1-D Navier-Stokes solution for the Space Shuttle Main Engine (SSME) thrust chamber with cross section area variations taken into account. It is shown that instability occurs when the mean pressure is set at 2000 psi with 30 percent disturbances. Instability also arises when the mean pressure is set at 2935 psi with 20 percent disturbances. The system with mean pressures and disturbances more adverse that these cases were shown to be unstable.

Chung, T. J.↗

NASA NTRS1990

Effects of turbulence mixing, variable properties, and vaporization on spray droplet combustion

Combustion of liquid fuels in the form of spray droplets is simulated numerically. Various vaporization models are examined as to their performance in finite element calculations involving a turbulent flow field. The Eulerian coordinate for the gas and Lagrangian coordinate for the liquid spray droplets are coupled through source terms being updated in the equations of continuity, momentum, and energy. The k-epsilon and modified eddy breakup models are used for simulating turbulent spray combustion flow field. Numerical results for the droplet trajectories, droplet heating, recirculation characteristics, and effects of evaporation models are evaluated. It is also shown that the finite element method is advantageous in dealing with complex geometries, complex boundary conditions, adaptive unstructured grids.

Kim, Y. M.↗

NASA NTRS1989

Finite element analysis in fluids; Proceedings of the Seventh International Conference on Finite Element Methods in Flow Problems, University of Alabama, Huntsville, Apr. 3-7, 1989

Recent advances in computational fluid dynamics are examined in reviews and reports, with an emphasis on finite-element methods. Sections are devoted to adaptive meshes, atmospheric dynamics, combustion, compressible flows, control-volume finite elements, crystal growth, domain decomposition, EM-field problems, FDM/FEM, and fluid-structure interactions. Consideration is given to free-boundary problems with heat transfer, free surface flow, geophysical flow problems, heat and mass transfer, high-speed flow, incompressible flow, inverse design methods, MHD problems, the mathematics of finite elements, and mesh generation. Also discussed are mixed finite elements, multigrid methods, non-Newtonian fluids, numerical dissipation, parallel vector processing, reservoir simulation, seepage, shallow-water problems, spectral methods, supercomputer architectures, three-dimensional problems, and turbulent flows.

Chung, T. J.↗