Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical calculations: variational”

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 181 records · Page 10

Stability of compressible three-dimensional boundary-layer flows

For compressible three-dimensional flow, the method of multiple scales to formulate the three-dimensional stability problem and determine the partial-differential equations governing variations of the amplitude and complex wavenumbers is used. A method for following one specific wave along its trajectory to ascertain the characteristics of the most unstable disturbance is proposed. Numerical results using the flow over the X-21 wing as calculated from the Kaups-Cebeci code will be presented.

Reed, H. L.↗

Numerical calculations of mass transfer flow in semi-detached binary systems

The details of the mass transfer flow near the inner Lagrangian point in a semidetached binary system are numerically calculated. A polytropic equation of state with n = 3/2 is used. The dependence of the mass transfer rate on the degree to which the star overfills its Roche lobe is calculated, and good agreement with previous analytic estimates is found. The variation of mass transfer rate which occurs if the binary system has a small eccentricity is calculated and is used to cast doubt on the model for superhumps in dwarf novae proposed by Papaloizou and Pringle (1979).

Edwards, D. A.↗

A numerical study of the unsteady processes associated with the type IV shock interaction

A time-dependent ADI formulation of the two dimensional Navier-Stokes equations coupled with the total variation diminishing (TVD) scheme is used to examine the type IV shock interaction time accurately. The experimentally measured and computationally calculated unsteadiness of the type IV interaction is shown to be related to the formation and shedding of shear layers in the shock layer, and the frequency of oscillation of the jet is a strong function of the shedding rate of these shear layers. For one particular geometry, an oscillation frequency of 1.4 kHz was calculated. The effects of shock impingement angle and shock strength on the location and magnitude of the peak pressure are also identified. In one case it is shown that a 1-deg change in the impinging shock angle can alter the location of the associated jet by almost 11 deg with a change in the peak surface pressure of 32 percent.

Lind, Charles A.↗

Computational Study of Electron-Molecule Collisions Related to Low-Temperature Plasmas

Computational study of electron-molecule collisions not only complements experimental measurements, but can also be used to investigate processes not readily accessible experimentally. A number of ab initio computational methods are available for these types of calculations. Here we describe a recently developed technique, the finite element Z-matrix method, Analogous to the R-matrix, method, it partitions the space into regions and employs real matrix elements. However, unlike the implementation of the R-matrix method commonly used in atomic and molecular physics, the Z-matrix method is fully variational. In the present implementation, a mixed basis of finite elements and Gaussians is used to represent the continuum electron, thus offering full flexibility without imposing fixed boundary conditions. Numerical examples include the electron-impact dissociation of N2 via the metastable A3Su+ state, a process which may be important in the lower thermosphere, and the dissociation of the CF radical, a process of interest to plasma etching. To understand the dissociation pathways, large scale quantum chemical calculations have been carried out for all target states which dissociate to the lowest five limits in the case of N2, and to the lowest two limits in the case of CF. For N2, the structural calculations clearly show the preference for predissociation if the initial state is the ground X1 Sg+ state, but direct dissociation appears to be preferable if the initial state is the A3Su+ state. Multi-configuration SCF target functions are used in the collisional calculation.

Huo, Winifred M.↗

Numerical experiments with a symmetric high-resolution shock-capturing scheme

Characteristic-based explicit and implicit total variation diminishing (TVD) schemes for the two-dimensional compressible Euler equations have recently been developed. This is a generalization of recent work of Roe and Davis to a wider class of symmetric (non-upwind) TVD schemes other than Lax-Wendroff. The Roe and Davis schemes can be viewed as a subset of the class of explicit methods. The main properties of the present class of schemes are that they can be implicit, and, when steady-state calculations are sought, the numerical solution is independent of the time step. In a recent paper, a comparison of a linearized form of the present implicit symmetric TVD scheme with an implicit upwind TVD scheme originally developed by Harten and modified by Yee was given. Results favored the symmetric method. It was found that the latter is just as accurate as the upwind method while requiring less computational effort. Currently, more numerical experiments are being conducted on time-accurate calculations and on the effect of grid topology, numerical boundary condition procedures, and different flow conditions on the behavior of the method for steady-state applications. The purpose here is to report experiences with this type of scheme and give guidelines for its use.

Yee, H. C.↗

Hingeless rotor theory and experiment on vibration reduction by periodic variation of conventional controls

The reduction of the n per rev. pitch-, roll- and vertical vibrations of an n-bladed rotor by n per rev. sinusoidal variations of the collective and cyclic controls is investigated. The numerical results presented refer to a four-bladed, 7.5-foot model and are based on frequency response tests conducted under an Army-sponsored research program. The following subjects are treated: extraction of the rotor transfer functions (.073R hub flapping and model thrust versus servo valve command, amplitude and phase), calculation of servo commands (volts) required to compensate .073R hub flapping (3P and 5P) and model thrust (4P), evaluation of the effect of the vibratory control inputs on blade loads, and theoretical prediction of the root flapbending moments generated by 0 to 5P perturbations of the feathering angle and rotor angle of attack. Five operating conditions are investigated covering advance ratios from approximately 0.2 to 0.85. The feasibility of vibration reduction by periodic variation on conventional controls is evaluated.

Sissingh, G. J.↗

Simulation studies of a model of high-density metallic hydrogen

Upper bounds for the ground-state energies of liquid and solid phases of metallic hydrogen and metallic deuterium have been calculated with variational methods and Monte Carlo techniques. At four densities (0.8, 1.2, 1.36, and 1.488) crystalline phases are clearly preferred in the sense that the energy difference, when compared to the liquid, is in excess of the errors inherent in the numerical procedures. At a fifth density (1.6), the energy differences between solid and liquid phases are smaller than these errors.

Mon, K. K.↗

Computational analysis of the SSME fuel preburner flow

A computational fluid dynamics model which simulates the steady state operation of the SSME fuel preburner is developed. Specifically, the model will be used to quantify the flow factors which cause local hot spots in the fuel preburner in order to recommend experiments whereby the control of undesirable flow features can be demonstrated. The results of a two year effort to model the preburner are presented. In this effort, investigating the fuel preburner flowfield, the appropriate transport equations were numerically solved for both an axisymmetric and a three-dimensional configuration. Continuum's VAST (Variational Solution of the Transport equations) code, in conjunction with the CM-1000 Engineering Analysis Workstation and the NASA/Ames CYBER 205, was used to perform the required calculations. It is concluded that the preburner operational anomalies are not due to steady state phenomena and must, therefore, be related to transient operational procedures.

Wang, T. S.↗

Computational Study of Flow Establishment in a Ram Accelerator

The temporal evolution of the combustion process established during projectile transition from the launch tube into the ram accelerator section containing an explosive hydrogen-oxygen-argon gas mixture is studied. The Navier-Stokes equations for chemically reacting flow are solved in a fully coupled manner, using an implicit, time accurate algorithm. The solution procedure is based on a spatially second order total variation diminishing scheme and a temporally second order, variable-step, backward differentiation formula method. The hydrogen-oxygen chemistry is modeled with a 9-species, 19-step mechanism. The accuracy of the solution method is first demonstrated by several benchmark calculations. Numerical simulations of two ram accelerator configurations are then presented. In particular, the temporal developments of shock-induced combustion and thrust forces are followed. Positive thrust is established in both cases; however, in one of the ram accelerator configurations studied, combustion in the boundary layer enhances its separation, ultimately resulting in unstart.

Yungster, S.↗

Numerical solutions of three-dimensional Navier-Stokes equations for closed-bluff bodies

Methods of grid adaption are reviewed and a method is developed with the capability of adaption to several flow variables. This method is based on a variational approach and is an algebraic method which does not require the solution of partial differential equations. The method has also been formulated in such a way that there is no need for any matrix inversion. The method is used in conjunction with the calculation of hypersonic flow over a blunt nose. The equations of motion are the compressible Navier-Stokes equations where all viscous terms are retained. They are solved by the MacCormack time-splitting method and a movie has been produced which shows simultaneously the transient behavior of the solution and the grid adaption. The results are compared with the experimental and other numerical results.

Abolhassani, Jamshid S.↗

Application of finite difference techniques to noise propagation in jet engine ducts

A finite difference formulation is presented for wave propagation in a rectangular two-dimensional duct without steady flow. The difference technique, which should be useful in the study of acoustically treated inlet and exhausts ducts used in turbofan engines, can readily handle acoustical flow field complications such as axial variations in wall impedance and cross section area. In the numerical analysis, the continuous acoustic field is lumped into a series of grid points in which the pressure and velocity at each grid point are separated into real and imaginary terms. An example calculation is also presented for the sound attenuation in a two-dimensional straight soft-walled suppressor.

Baumeister, K. J.↗

Thermal microwave emission from half-space random media

Brightness temperatures resulting from microwave thermal emission from a half-space random medium are calculated. The random medium has a nonuniform temperature profile and is characterized by correlation functions that possess both vertical and lateral variations. Radiative transfer equations are derived. They are solved with an iterative integral equation approach for small scattering albedo and with a numerical approach for general cases. New results are illustrated, discussed, and compared with various special cases. Phenomena caused by resonant scattering, which is absent under Rayleigh approximations, are revealed. Scattering in the lateral direction contributes to the decrease of the equivalent reflectivity and the increase of the brightness temperature.

Tsang, L.↗

Collapse and fragmentation of rotating, adiabatic clouds

A numerical hydrodynamics code has been used to calculate the collapse of rotating, adiabatic clouds. The three-dimensional nature of the calculation allows the clouds to fragment in the dynamic collapse phase. Clouds with adiabatic exponent of 7/5 and initial cos(2 phi) density variations fragment into binary systems if the initial ratio of thermal to gravitational energy is small (about 0.05). Clouds with higher thermal energy, however, damp the density variation and form near-equilibrium ellipsoids, with ratios of rotational to gravitational energy less than the critical value for dynamic growth of nonaxisymmetry in Maclaurin spheroids. Even with an adiabatic pressure law, dynamic fragmentation of a collapsing cloud is possible, implying for star formation theory that the low thermal energy fragments produced in isothermal collapse calculations may undergo a subsequent dynamic fragmentation in the nonisothermal regime.

Boss, A. P.↗

A numerical study of the stability of radiative shocks

Attention is given to the oscillatory instability of optically thin radiative shocks in time-dependent numerical calculations of accretion flows onto degenerate dwarfs. The present nonlinear calculations yield good quantitative agreement with the linear results obtained for oscillation frequencies, damping rates, and critical alpha-values. The fundamental mode and the first overtone in the shock radius and luminosity variations can be clearly identified, and evidence is sometimes seen for the second overtone. Time-dependent calculations are also performed which include additional physics relevant to degenerate dwarf accretion, such as electron thermal conduction, unequal electron and ion temperatures, Compton cooling, and relativistic corrections to the bremsstrahlung cooling law. All oscillatory modes are found to be damped, and hence stable, in the case of a 1-solar mass white dwarf accreting in spherical symmetry.

Imamura, J. N.↗

Surface deformation, gravity and the geoid from a three-dimensional convection model at low Rayleigh numbers

The two-dimensional surface deformation, gravity field, and geoid are calculated from the temperature fields of a number of numerical models of constant-viscosity three-dimensional convective flows, heated from within and from below. The surface deformation and the geoid are insensitive to the short-wavelength features of the temperature variation; the gravity field is less smooth, although it still does not clearly indicate the narrowness of the upwelling and downwelling regions at large Rayleigh numbers. The results indicate that the convective upwelling beneath regions showing positive geoid and residual depth anomalies is more localized than the horizontal extent of these features would suggest. The results are relevant to attempts to understand the behavior of the earth mantle.

Craig, Claire Harvey↗

Computational Study of Flow Establishment in Hypersonic Pulse Facilities

This paper presents a study of the temporal evolution of the combustion flowfield established by the interaction of ram-accelerator-type projectiles with an explosive gas mixture accelerated to hypersonic speeds in an expansion tube. The Navier-Stokes equations for a chemically reacting gas are solved in a fully coupled manner using an implicit, time accurate algorithm. The solution procedure is based on a spatially second order, total variation diminishing (TVD) scheme and a temporally second order, variable-step, backward differentiation formula method. The hydrogen-oxygen chemistry is modeled with a 9-species, 19-step mechanism. The accuracy of the solution method is first demonstrated by several benchmark calculations. Numerical simulations of expansion tube flowfields are then presented for two different configurations. In particular, the development of the shock-induced combustion process is followed. In one case, designed to ensure ignition only in the boundary layer, the lateral extent of the combustion front during the initial transient phase was surprisingly large. The time histories of the calculated thrust and drag forces on the ram accelerator projectile are also presented.

Yungster, S.↗

Design of an Object-Oriented Turbomachinery Analysis Code: Initial Results

Performance prediction of turbomachines is a significant part of aircraft propulsion design. In the conceptual design stage, there is an important need to quantify compressor and turbine aerodynamic performance and develop initial geometry parameters at the 2-D level prior to more extensive Computational Fluid Dynamics (CFD) analyses. The Object-oriented Turbomachinery Analysis Code (OTAC) is being developed to perform 2-D meridional flowthrough analysis of turbomachines using an implicit formulation of the governing equations to solve for the conditions at the exit of each blade row. OTAC is designed to perform meanline or streamline calculations; for streamline analyses simple radial equilibrium is used as a governing equation to solve for spanwise property variations. While the goal for OTAC is to allow simulation of physical effects and architectural features unavailable in other existing codes, it must first prove capable of performing calculations for conventional turbomachines.OTAC is being developed using the interpreted language features available in the Numerical Propulsion System Simulation (NPSS) code described by Claus et al (1991). Using the NPSS framework came with several distinct advantages, including access to the pre-existing NPSS thermodynamic property packages and the NPSS Newton-Raphson solver. The remaining objects necessary for OTAC were written in the NPSS framework interpreted language. These new objects form the core of OTAC and are the BladeRow, BladeSegment, TransitionSection, Expander, Reducer, and OTACstart Elements. The BladeRow and BladeSegment consumed the initial bulk of the development effort and required determining the equations applicable to flow through turbomachinery blade rows given specific assumptions about the nature of that flow. Once these objects were completed, OTAC was tested and found to agree with existing solutions from other codes; these tests included various meanline and streamline comparisons of axial compressors and turbines at design and off-design conditions.

compressors↗

Polyatomic molecular Dirac-Hartree-Fock calculations with Gaussian basis sets

Numerical methods have been used successfully in atomic Dirac-Hartree-Fock (DHF) calculations for many years. Some DHF calculations using numerical methods have been done on diatomic molecules, but while these serve a useful purpose for calibration, the computational effort in extending this approach to polyatomic molecules is prohibitive. An alternative more in line with traditional quantum chemistry is to use an analytical basis set expansion of the wave function. This approach fell into disrepute in the early 1980's due to problems with variational collapse and intruder states, but has recently been put on firm theoretical foundations. In particular, the problems of variational collapse are well understood, and prescriptions for avoiding the most serious failures have been developed. Consequently, it is now possible to develop reliable molecular programs using basis set methods. This paper describes such a program and reports results of test calculations to demonstrate the convergence and stability of the method.

Dyall, Kenneth G.↗