SEARCH · Search NASA

Results for “problems”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,171 records · Page 65

All sources NASA NTRS DOE OSTI

NASA NTRS1994

The algebraic multigrid projection for eigenvalue problems; backrotations and multigrid fixed points

The periods of the theorem for the algebraic multigrid projection (MGP) for eigenvalue problems, and of the multigrid fixed point theorem for multigrid cycles combining MGP with backrotations, are presented. The MGP and the backrotations are central eigenvector separation techniques for multigrid eigenvalue algorithms. They allow computation on coarse levels of eigenvalues of a given eigenvalue problem, and are efficient tools in the computation of eigenvectors.

Costiner, Sorin↗

NASA NTRS1994

Knowledge based translation and problem solving in an intelligent individualized instruction system

An Intelligent Individualized Instruction I(sup 3) system is being built to provide computerized instruction. We present the roles of a translator and a problem solver in an intelligent computer system. The modular design of the system provides for easier development and allows for future expansion and maintenance. CLIPS modules and classes are utilized for the purpose of the modular design and inter module communications. CLIPS facts and rules are used to represent the system components and the knowledge base. CLIPS provides an inferencing mechanism to allow the I(sup 3) system to solve problems presented to it in English.

Jung, Namho↗

NASA NTRS1995

High-performance parallel analysis of coupled problems for aircraft propulsion

This research program deals with the application of high-performance computing methods to the numerical simulation of complete jet engines. The program was initiated in 1993 by applying two-dimensional parallel aeroelastic codes to the interior gas flow problem of a by-pass jet engine. The fluid mesh generation, domain decomposition and solution capabilities were successfully tested. Attention was then focused on methodology for the partitioned analysis of the interaction of the gas flow with a flexible structure and with the fluid mesh motion driven by these structural displacements. The latter is treated by an ALE technique that models the fluid mesh motion as that of a fictitious mechanical network laid along the edges of near-field fluid elements. New partitioned analysis procedures to treat this coupled 3-component problem were developed in 1994. These procedures involved delayed corrections and subcycling, and have been successfully tested on several massively parallel computers. For the global steady-state axisymmetric analysis of a complete engine we have decided to use the NASA-sponsored ENG10 program, which uses a regular FV-multiblock-grid discretization in conjunction with circumferential averaging to include effects of blade forces, loss, combustor heat addition, blockage, bleeds and convective mixing. A load-balancing preprocessor for parallel versions of ENG10 has been developed. It is planned to use the steady-state global solution provided by ENG10 as input to a localized three-dimensional FSI analysis for engine regions where aeroelastic effects may be important.

Felippa, C. A.↗

NASA NTRS1995

The initial-value problem for viscous channel flows

Plane viscous channel flows are perturbed and the ensuing initial-value problems are investigated in detail. Unlike traditional methods where traveling wave normal modes are assumed for solution, this works offers a means whereby completely arbitrary initial input can be specified without having to resort to eigenfunction expansions. The full temporal behavior, including both early time transients and the long time asymptotics, can be determined for any initial disturbance. Effects of three-dimensionality can be assessed. The bases for the analysis are: (a) linearization of the governing equations; (b) Fourier decomposition in the spanwise and streamwise directions of the flow; and (c) direct numerical integration of the resulting partial differential equations. All of the stability data that are known for such flows can be reproduced. Also, the optimal initial condition can be determined in a straight forward manner and such optimal conditions clearly reflect transient growth data that is easily determined by a rational choice of a basis for the initial conditions. Although there can be significant transient growth for subcritical values of the Reynolds number using this approach it does not appear possible that arbitrary initial conditions will lead to the exceptionally large transient amplitudes that have been determined by optimization of normal modes. The approach is general and can be applied to other classes of problems where only a finite discrete spectrum exists, such as the boundary layer for example.

Criminale, W. O.↗

NASA NTRS1994

Separation of crack extension modes in composite delamination problems

This work concerns fracture mechanics modeling of composite delamination problems. In order to predict delamination resistance, an applied stress intensity factor, K, or energy release rate, G, must be compared to a mode-dependent critical value of K or G from experiment. In the interfacial fracture analysis of most applications and some tests, the mode of crack extension is not uniquely defined. It is instead a function of distance from the crack tip due to the oscillating singularity existing at the tip. In this work, a consistent method is presented of extracting crack extension modes in such cases. In particular, use of the virtual crack closure technique (VCCT) to extract modes of crack extension is studied for cases of a crack along the interface between two in-plane orthotropic materials. Modes of crack extension extracted from oscillatory analyses using VCCT are a function of the virtual crack extension length, delta. Most existing efforts to obtain delta-independent modes of crack extension involve changing the analysis in order to eliminate its oscillatory nature. One such method involves changing one or more properties of the layers to make the oscillatory exponent parameter, epsilon, equal zero. Standardized application of this method would require consistent criteria for identifying which properties can be altered without changing the physical aspects of the problem. Another method involves inserting a thin homogeneous layer (typically referred to as a resin interlayer) along the interface and placing the crack within it. The drawbacks of this method are that it requires increased modeling effort and introduces the thickness of the interlayer as an additional length parameter. The approach presented here does not attempt to alter the interfacial fracture analysis to eliminate its oscillatory behavior. Instead, the argument is made that the oscillatory behavior is non-physical and that if its effects were separated from VCCT quantities, then consistent, delta-independent modes of crack extension could be defined. Knowledge of the near-tip fields in a planar orthotropic material interfacial fracture analysis is used to determine the explicit delta dependence of VCCT parameters. Once this delta dependence is determined, energy release rates are defined with this delta dependence factored out. This modified VCCT method is applied to results from two finite element test cases. It is shown that, as predicted, delta-independent modes of crack extension result. The modified VCCT approach shows potential as a consistent method of extracting crack extension modes. It uses the same information from a finite element analysis (i.e., nodal forces and displacements) as the traditional VCCT method does. The A-independent modes extracted using the modified VCCT approach can also be used as guides to test the convergence of finite element solutions.

Beuth, Jack L.↗

NASA NTRS1994

Theoretical studies of spectroscopic problems of importance for atmospheric radiation measurements

Many of the instruments used to deduce the physical parameters of the Earth's atmosphere necessary for climate studies or for pollution monitoring (for instance, temperature versus pressure or number densities of trace molecules) rely on the existence of accurate spectroscopic data and an understanding of the physical processes responsible for the absorption or emission of radiation. During the summer, research was either continued or begun on three distinct problems: (1) an improved theoretical framework for the calculation of the far-wing absorption of allowed spectral lines; (2) a refinement of the calculation of the collision-induced fundamental spectrum of N2; and (3) an investigation of possible line-mixing effects in the fundamental spectrum of CH4. Progress in these three areas is summarized below. During the past few years, we have developed a theoretical framework for the calculation of the absorption of radiation by the far wings of spectral lines. Such absorption due to water vapor plays a crucial role in the greenhouse effect as well as limiting the retrieval of temperature profiles from satellite data. Several improvements in the theory have been made and the results are being prepared for publication. Last year we published results for the theoretical calculation of the absorption of radiation due to the dipoles induced during binary collisions of N2 molecules using independently measured molecular parameters; the results were in reasonable agreement with experimental data. However, recent measurements have revealed new fine structure that has been attributed to line-mixing effects. We do not think that this is correct, rather that the structure results from short-range anisotropic dipoles. We are in the process of including this refinement in our theoretical calculation in order to compare with the new experimental data. Subtle changes in the spectra of CH4 measured by researchers at Langley have also been attributed to line-mixing effects. By analyzing the same spectral lines we have attempted to verify or rule out possible line-mixing mechanisms. Due to the complexity and richness of the spectrum of this highly symmetric molecule, as well as the small magnitude of the effects, a detailed first-principle calculation of the mixing is a difficult problem. Before such a program is undertaken it is important to glean as much information as possible concerning the possible mechanisms by a systematic analysis of the existing data.

Tipping, Richard H.↗

NASA NTRS1995

Multigrid solution strategies for adaptive meshing problems

This paper discusses the issues which arise when combining multigrid strategies with adaptive meshing techniques for solving steady-state problems on unstructured meshes. A basic strategy is described, and demonstrated by solving several inviscid and viscous flow cases. Potential inefficiencies in this basic strategy are exposed, and various alternate approaches are discussed, some of which are demonstrated with an example. Although each particular approach exhibits certain advantages, all methods have particular drawbacks, and the formulation of a completely optimal strategy is considered to be an open problem.

Mavriplis, Dimitri J.↗

NASA NTRS1995

A procedure for automating CFD simulations of an inlet-bleed problem

A procedure was developed to improve the turn-around time for computational fluid dynamics (CFD) simulations of an inlet-bleed problem involving oblique shock-wave/boundary-layer interactions on a flat plate with bleed into a plenum through one or more circular holes. This procedure is embodied in a preprocessor called AUTOMAT. With AUTOMAT, once data for the geometry and flow conditions have been specified (either interactively or via a namelist), it will automatically generate all input files needed to perform a three-dimensional Navier-Stokes simulation of the prescribed inlet-bleed problem by using the PEGASUS and OVERFLOW codes. The input files automatically generated by AUTOMAT include those for the grid system and those for the initial and boundary conditions. The grid systems automatically generated by AUTOMAT are multi-block structured grids of the overlapping type. Results obtained by using AUTOMAT are presented to illustrate its capability.

Chyu, Wei J.↗

NASA NTRS1995

Application of the space-time conservation element and solution element method to two-dimensional advection-diffusion problems

The existing 2-D alpha-mu scheme and alpha-epsilon scheme based on the method of space-time conservation element and solution element, which were constructed for solving the linear 2-D unsteady advection-diffusion equation and unsteady advection equation, respectively, are tested. Also, the alpha-epsilon scheme is modified to become the V-E scheme for solving the nonlinear 2-D inviscid Burgers equation. Numerical solutions of six test problems are presented in comparison with their exact solutions or numerical solutions obtained by traditional finite-difference or finite-element methods. It is demonstrated that the 2-D alpha-mu, alpha-epsilon, and nu-epsilon schemes can be used to obtain numerical results which are more accurate than those based on some of the traditional methods but without using any artificial tuning in the computation. Similar to the previous 1-D test problems, the high accuracy and simplicity features of the space-time conservation element and solution element method have been revealed again in the present 2-D test results.

Wang, Xiao-Yen↗

NASA NTRS1995

ICASE/LaRC Workshop on Benchmark Problems in Computational Aeroacoustics (CAA)

The proceedings of the Benchmark Problems in Computational Aeroacoustics Workshop held at NASA Langley Research Center are the subject of this report. The purpose of the Workshop was to assess the utility of a number of numerical schemes in the context of the unusual requirements of aeroacoustical calculations. The schemes were assessed from the viewpoint of dispersion and dissipation -- issues important to long time integration and long distance propagation in aeroacoustics. Also investigated were the effect of implementation of different boundary conditions. The Workshop included a forum in which practical engineering problems related to computational aeroacoustics were discussed. This discussion took the form of a dialogue between an industrial panel and the workshop participants and was an effort to suggest the direction of evolution of this field in the context of current engineering needs.

Hardin, Jay C.↗

NASA NTRS1995

Application of low dissipation and dispersion Runge-Kutta schemes to benchmark problems in computational aeroacoustics

We investigate accurate and efficient time advancing methods for computational aeroacoustics, where non-dissipative and non-dispersive properties are of critical importance. Our analysis pertains to the application of Runge-Kutta methods to high-order finite difference discretization. In many CFD applications, multi-stage Runge-Kutta schemes have often been favored for their low storage requirements and relatively large stability limits. For computing acoustic waves, however, the stability consideration alone is not sufficient, since the Runge-Kutta schemes entail both dissipation and dispersion errors. The time step is now limited by the tolerable dissipation and dispersion errors in the computation. In the present paper, it is shown that if the traditional Runge-Kutta schemes are used for time advancing in acoustic problems, time steps greatly smaller than that allowed by the stability limit are necessary. Low Dissipation and Dispersion Runge-Kutta (LDDRK) schemes are proposed, based on an optimization that minimizes the dissipation and dispersion errors for wave propagation. Optimizations of both single-step and two-step alternating schemes are considered. The proposed LDDRK schemes are remarkably more efficient than the classical Runge-Kutta schemes for acoustic computations. Numerical results of each Category of the Benchmark Problems are presented. Moreover, low storage implementations of the optimized schemes are discussed. Special issues of implementing numerical boundary conditions in the LDDRK schemes are also addressed.

Hu, F. Q.↗

NASA NTRS1995

Results of two new methods for aeroacoustics benchmark problems

Two new methods for the numerical solution of conservation laws (the Euler equations in particular) are presented: a uniformly second-order accurate upwind scheme and a third-order accurate centered scheme. Results of these schemes are shown for problems 1, 2, and 5 of this workshop's benchmark problems.

Huynh, H. T.↗

NASA NTRS1995

Application of traditional CFD methods to nonlinear computational aeroacoustics problems

This paper describes an implementation of a high order finite difference technique and its application to the category 2 problems of the ICASE/LaRC Workshop on Computational Aeroacoustics (CAA). Essentially, a popular Computational Fluid Dynamics (CFD) approach (central differencing, Runge-Kutta time integration and artificial dissipation) is modified to handle aeroacoustic problems. The changes include increasing the order of the spatial differencing to sixth order and modifying the artificial dissipation so that it does not significantly contaminate the wave solution. All of the results were obtained from the CM5 located at the Numerical Aerodynamic Simulation Laboratory. lt was coded in CMFortran (very similar to HPF), using programming techniques developed for communication intensive large stencils, and ran very efficiently.

Chyczewski, Thomas S.↗

NASA NTRS1994

Magnetocentrifugally driven flows from young stars and disks. 2: Formulation of the dynamical problem

We formulate the dynamical problem of a cool wind centrifugally driven from the magnetic interface of a young star and an adjoining Keplerian disk. We examine the situation for mildly accreting T Tauri stars that rotate slowly as well as rapidly accreting protostars that rotate near break-up. In both cases a wind can be driven from a small X-region just outside the stellar magnetopause, where the field lines assume an open geometry and are rooted to material that rotates at an angular speed equal both to the local Keplerian value and to the stellar angular speed. Assuming axial symmetry for the ideal magnetohydrodynamic flow, which requires us to postpone asking how the (lightly ionized) gas is loaded onto field lines, we can formally integrate all the governing equations analytically except for a partial equation that describes how streamlines spread in the meridional plane. Apart from the difficulty of dealing with PDEs of mixed type, finding the functional forms of the conserved quantities along streamlines - the ratio beta of magnetic field to mass flux, the specific energy H of the fluid in the rotating frame, and the total specific angular momentum J carried in the matter and the field - constitutes a standard difficulty in this kind of (Grad-Shafranov) formalism. Fortunately, because the ratio of the thermal speed of the mass-loss regions to the Keplerian speed of rotation of the interface constitutes a small parameter epsilon, we can attack the overall problem by the method of matched asymptotic expansions. This procedure leads to a natural and systematic technique for obtaining the relevant functional dependences of beta, H, and J. Moreover, we are able to solve analytically for the properties of the flow emergent from the small transsonic region driven by gas pressure without having to specify the detailed form of any of the conserved functions, beta, H, and J. This analytical solution provides inner boundary conditions for the numerical computation in a companion paper by Najita & Shu of the larger region where the main acceleration to terminal speeds occurs.

Shu, Frank H.↗

NASA NTRS1994

The Eady problem for a basic state with zero PV gradient but beta does not equal 0

The classic Eady problem is modified to include beta does not equal 0, but with the basic distributions of temperature and zonal flow adjusted to preserve zero meridional gradients of basic-state potential vorticity in the fluid interior. Much of the mathematical simplicity of the classic problem is retained; however, the results differ in important ways. Specifically, the instability now has a long-wave cutoff in addition to the traditional short-wave cutoff. The former is associated with the fact that the phase speeds of the edge waves begin to differ so much as wavenumber is reduced that the two edge waves can no longer interact in order to form unstable modes. For the unstable modes, this manifests itself in that the steering level for unstable modes is always below the middle of the fluid and approaches the lower boundary near the long-wave cutoff. Relatedly, the amplitude of the unstable geopotential perturbations is larger at the upper boundary than at the lower boundary. Finally, below the long-wave cutoff, one of the neutral waves has a phase speed that becomes increasingly easterly as wavenumber decreases. This allows a resonant response to planetary-scale stationary forcing.

Lindzen, Richard S.↗

NASA NTRS1994

A mapping method for the gravitational few-body problem with dissipation

Recently a new class of numerical integration methods - 'mixed variable symplectic integrators' - has been introduced for studying long-term evolution in the conservative gravitational few-body problem. These integrators are an order of magnitude faster than conventional ordinary differential equations (ODE) integration methods. Here we present a simple modification of this method to include small non-gravitational forces. The new scheme provides a similar advantage of computational speed for a larger class of problems in Solar System dynamics.

Malhotra, Renu↗

NASA NTRS1992

Analytic solution of the thickness problem of a rectangular wing in steady subsonic flow

An analytic solution of the thickness problem of a rectangular wing with parabolic airfoil section in three-dimensional flow is presented. The free-air solution is obtained by integrating the equation of the axial perturbation velocity. The Prandtl-Glauert rule can be used to derive the subsonic solution. Parts of the free-air solution are verified by taking the limit of the axial peturbation velocity on the model surface as the wing span goes to infinity. Pressure coefficients are studied for a selected wing geometry in the flow field. The solution of the thickness problem of a rectangular wing in a rectangular wind tunnel is derived using the free-air solution and the method of images.

Ulbrich, N.↗

NASA NTRS1995

What is the problem with He-3

We consider the galactic evolutionary history of He-3 in models which deplete deuterium by as much as a factor of 2 to approximately 15 from its primordial value to its present-day observed value in the interstellar medium (ISM). We show that when He-3 production in low-mass stars (1-3 solar mass) is included over the history of the galaxy, He-3 is greatly overproduced and exceeds the inferred solar values and the abundances determined in galactic H II regions. Furthermore, the ISM abundances show a disturbing dispersion which is difficult to understand from the point of view of standard chemical evolution models. In principle, resolution of the problem may lie in either (1) the calculated He-3 production in low-mass stars; (2) the observations of the He-3 abundance; or (3) an observational bias twoard regions of depleted He-3. Since He-3 observations in planetary nebula support the calculated He-3 production in low-mass stars, option (1) is unlikely. We will argue for option (3) since the He-3 interstellar observations are indeed made in regions dominated by massive stars in which He-3 is destroyed. In conclusion, we note that the problem with He-3 seems to be galactic and not cosmological.

Olive, Keith A.↗