Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical relativity”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

The numerical viscosity of entropy stable schemes for systems of conservation laws

Discrete approximations to hyperbolic systems of conservation laws are studied. The amount of numerical viscosity present in such schemes, is quantified and related to their entropy stability by means of comparison. To this end, conservative schemes which are also entropy conservative are constructed. These entropy conservative schemes enjoy second-order accuracy; moreover, they admit a particular interpretation within the finite-element frameworks, and hence can be formulated on various mesh configurations. It is then shown that conservative schemes are entropy stable if and only if they contain more viscosity than the mentioned above entropy conservative ones.

Tadmor, E.↗
NASA NTRS1986

Middle Atmosphere Program. Handbook for MAP, volume 20

Various topics related to investigations of the middle atmosphere are discussed. Numerical weather prediction, performance characteristics of weather profiling radars, determination of gravity wave and turbulence parameters, case studies of gravity-wave propagation, turbulence and diffusion due to gravity waves, the climatology of gravity waves, mesosphere-stratosphere-troposphere radar, antenna arrays, and data management techniques are among the topics discussed.

Bowhill, S. A.↗
NASA NTRS1986

Satellite attitude motion models for capture and retrieval investigations

The primary purpose of this research is to provide mathematical models which may be used in the investigation of various aspects of the remote capture and retrieval of uncontrolled satellites. Emphasis has been placed on analytical models; however, to verify analytical solutions, numerical integration must be used. Also, for satellites of certain types, numerical integration may be the only practical or perhaps the only possible method of solution. First, to provide a basis for analytical and numerical work, uncontrolled satellites were categorized using criteria based on: (1) orbital motions, (2) external angular momenta, (3) internal angular momenta, (4) physical characteristics, and (5) the stability of their equilibrium states. Several analytical solutions for the attitude motions of satellite models were compiled, checked, corrected in some minor respects and their short-term prediction capabilities were investigated. Single-rigid-body, dual-spin and multi-rotor configurations are treated. To verify the analytical models and to see how the true motion of a satellite which is acted upon by environmental torques differs from its corresponding torque-free motion, a numerical simulation code was developed. This code contains a relatively general satellite model and models for gravity-gradient and aerodynamic torques. The spacecraft physical model for the code and the equations of motion are given. The two environmental torque models are described.

Cochran, John E., Jr.↗
NASA NTRS1987

The numerical viscosity of entropy stable schemes for systems of conservation laws. I

Discrete approximations to hyperbolic systems of conservation laws are studied. The amount of numerical viscosity present in such schemes is quantified and related to their entropy stability by means of comparison.To this end, conservative schemes which are also entropy-conservative are constructed. These entropy-conservative schemes enjoy second-order accuracy; moreover, they can be interpreted as piecewise-linear finite-element methods, and hence can be formulated on various mesh configurations. It is then shown that conservative schemes are entropy stable, if and (for three-point schemes) only they contain more viscosity than that present in the above-mentioned entropy-conservative ones.

Tadmor, Eitan↗
NASA NTRS1989

Transonic aerodynamic design experience

Advancements have occurred in transonic numerical simulation that place aerodynamic performance design into a relatively well developed status. Efficient broad band operating characteristics can be reliably developed at the conceptual design level. Recent aeroelastic and separated flow simulation results indicate that systematic consideration of an increased range of design problems appears promising. This emerging capability addresses static and dynamic structural/aerodynamic coupling and nonlinearities associated with viscous dominated flows.

Bonner, E.↗
NASA NTRS1989

The length scale for sub-grid-scale parameterization with anisotropic resolution

Use of the Smagorinsky eddy-viscosity formulation and related schemes for subgrid-scale parameterization of large eddy simulation models requires specification of a single length scale, earlier related by Lilly to the scale of filtering and/or numerical resolution. An anisotropic integration of the Kolmogoroff enstrophy spectrum allows generalization of that relationship to anisotropic resolution. It is found that the Deardorff assumption is reasonably accurate for small anisotropies and can be simply improved for larger values.

Lilly, Douglas K.↗
NASA NTRS1988

Zonal techniques for flowfield simulation about aircraft

A technique for performing conservative flowfield calculations on zonal meshes is described. The underlying flow solver is an implicit, upwind finite volume scheme which can incorporate either a perfect gas or an equilibrium air equation of state. Two different approaches which yield identical results, in terms of performing a conservative flux calculation on a zonal interface, are described and compared in terms of numerical efficiency. The capability of the method to handle relatively complex geometries is demonstrated by considering the flowfield about a model SR71 aircraft.

Walters, Robert W.↗
NASA NTRS1989

Protostellar formation in rotating interstellar clouds. VIII - Inner core formation

The results are presented of a variety of spherically symmetric one-dimensional (1D) calculations intended to determine the robustness of the dynamical hiccup phenomenon in protostellar cores. The 1D models show that the phenomenon is relatively insensitive to changes in the equations of state, numerical resolution, initial density and temperature, and the radiative transfer approximation. In 1D, the hiccup results in an explosive destruction of the entire inner protostellar core. Inner core formation is studied with a sequence of three-dimensional models which show that rapid inner core rotation stabilizes the hiccup instability. Instead, the inner core becomes quite flat and undergoes a cycle of binary fragmentation, binary decay into a single object surrounded by a bar, breakup of the bar into a binary, etc. When lesser amounts of rotation are involved, the inner core does hiccup somewhat, but mass is ejected in only a few directions, leading to several broad streams of ejecta.

Boss, Alan P.↗
NASA NTRS1990

Computation of three-dimensional turbulent boundary layers with heat transfer in a plane of symmetry using embedded wall-layer functions

In the calculation of turbulent boundary layers, a large number of mesh points are required to adequately resolve the intense variation in the velocity and enthalpy in the near-wall region. A substantial reduction in computational effort may be realized by representing the velocity and enthalpy profiles in the wall layer by analytical embedded functions. The effectively inviscid flow in the outer part of the boundary layer may then be resolved by employing a relatively coarser mesh. To obtain complete profiles, the outer numerical solution is matched asymptotically to the inner wall-layer analytical solution. To date, this approach has been restricted to two-dimensional flows; in the present study, a method which may be utilized for turbulent boundary layers with heat transfer in a plane of symmetry is developed as a first step in the application of the embedded-function method to full three-dimensional flows. The present method uses only about half as many mesh points as that required in a conventional procedure, which calculates the flow all the way to the wall, but there is no degradation in accuracy of the computed results.

Degani, A. T.↗
NASA NTRS1990

Energy stability of thermocapillary convection in a model of the float-zone crystal-growth process

Energy stability theory has been applied to a basic state of thermocapillary convection occurring in a cylindrical half-zone of finite length to determine conditions under which the flow will be stable. Because of the finite length of the zone, the basic state must be determined numerically. Instead of obtaining stability criteria by solving the related Euler-Lagrange equations, the variational problem is attacked directly by discretization of the integrals in the energy identity using finite differences. Results of the analysis are values of the Marangoni number below which axisymmetric disturbances to the basic state will decay, for various values of the other parameters governing the problem.

Shen, Y.↗
NASA NTRS1990

Two dimensional sound propagation over a ridge

The sound propagation over an arbitrarily-shaped ridge situated on an infinite plane is treated using the Boundary Element Method, for which the Helmholtz integral equation in two-dimensional half space is formulated. A previously developed computational scheme for the acoustic problem at wave length shorther than the characteristic dimension of the radiator or scatterer is applied to the propagation of relatively low frequency sound over large terrain features. The numerical results using the Boundary Element Method are compared to those using algorithms for the Matched Asymptotic Expansion technique.

Eversman, W.↗
NASA NTRS1992

CFD lends the government a hand

The present survey of important and novel CFD applications being developed and implemented by U.S. Government contractors gives attention to naval vessel flow-modeling, Army ballistic and rotary wing aerodynamics, and NASA hypersonic vehicle related applications of CFD. CFD-generated knowledge of numerical algorithms, fluid motion, and supercomputer use is being incorporated into such additional areas as computational electromagnetics and acoustics. Attention is presently given to CFD methods' development status in such fields as submarine boundary layers, hypersonic kinetic energy projectile shock structures, helicopter main rotor tip flows, and National Aerospace Plane aerothermodynamics.

Lekoudis, Spiro↗
NASA NTRS1993

A hybrid structured-unstructured grid method for unsteady turbomachinery flow computations

A hybrid grid technique for the solution of 2D, unsteady flows is developed. This technique is capable of handling complex, multiple component geometries in relative motion, such as those encountered in turbomachinery. The numerical approach utilizes a mixed structured-unstructured zonal grid topology along with modeling equations and solution methods that are most appropriate in the individual domains, therefore combining the advantages of both structured and unstructured grid techniques.

Mathur, Sanjay R.↗
NASA NTRS1993

Turbulent structure of a wingtip vortex in the near field

The turbulent rollup of a vortex generated by a rectangular wing has been investigated. Extensive mean and turbulence measurements of the flowfield on a wingtip and in the near field have been completed. Velocity fluctuation measurements show that the near-field core is not laminar. A large axial velocity excess was found to exist in the core of the vortex. A momentum balance in the near-field of the wingtip showed that the magnitude of the core Reynolds-stress gradient terms are the same order as the largest terms in the governing equations. Navier-Stokes computations of the identical configuration, including wind tunnel walls and using measured inflow and outflow boundary conditions, reproduced many of the features of the experiment. Inherent limitations of the Baldwin-Barth turbulence model combined with limited grid resolution caused the computed vortex core to be more diffuse than desired. The momentum balance also demonstrated that the level of numerically generated false diffusion in the vortex core is relatively high.

Zilliac, Gregory G.↗
NASA NTRS1995

Viscous-flow analysis of a subsonic transport aircraft high-lift system and correlation with flight data

High-lift system aerodynamics has been gaining attention in recent years. In an effort to improve aircraft performance, comprehensive studies of multi-element airfoil systems are being undertaken in wind-tunnel and flight experiments. Recent developments in Computational Fluid Dynamics (CFD) offer a relatively inexpensive alternative for studying complex viscous flows by numerically solving the Navier-Stokes (N-S) equations. Current limitations in computer resources restrict practical high-lift N-S computations to two dimensions, but CFD predictions can yield tremendous insight into flow structure, interactions between airfoil elements, and effects of changes in airfoil geometry or free-stream conditions. These codes are very accurate when compared to strictly 2D data provided by wind-tunnel testing, as will be shown here. Yet, additional challenges must be faced in the analysis of a production aircraft wing section, such as that of the NASA Langley Transport Systems Research Vehicle (TSRV). A primary issue is the sweep theory used to correlate 2D predictions with 3D flight results, accounting for sweep, taper, and finite wing effects. Other computational issues addressed here include the effects of surface roughness of the geometry, cove shape modeling, grid topology, and transition specification. The sensitivity of the flow to changing free-stream conditions is investigated. In addition, the effects of Gurney flaps on the aerodynamic characteristics of the airfoil system are predicted.

Potter, R. C.↗
NASA NTRS1998

Bi-Level Integrated System Synthesis (BLISS)

BLISS is a method for optimization of engineering systems by decomposition. It separates the system level optimization, having a relatively small number of design variables, from the potentially numerous subsystem optimizations that may each have a large number of local design variables. The subsystem optimizations are autonomous and may be conducted concurrently. Subsystem and system optimizations alternate, linked by sensitivity data, producing a design improvement in each iteration. Starting from a best guess initial design, the method improves that design in iterative cycles, each cycle comprised of two steps. In step one, the system level variables are frozen and the improvement is achieved by separate, concurrent, and autonomous optimizations in the local variable subdomains. In step two, further improvement is sought in the space of the system level variables. Optimum sensitivity data link the second step to the first. The method prototype was implemented using MATLAB and iSIGHT programming software and tested on a simplified, conceptual level supersonic business jet design, and a detailed design of an electronic device. Satisfactory convergence and favorable agreement with the benchmark results were observed. Modularity of the method is intended to fit the human organization and map well on the computing technology of concurrent processing.

Sobieszczanski-Sobieski, Jaroslaw↗
NASA NTRS1999

Potential Flow Interactions With Directional Solidification

The effect of convective melt motion on the growth of morphological instabilities in crystal growth has been the focus of many studies in the past decade. While most of the efforts have been directed towards investigating the linear stability aspects, relatively little attention has been devoted to experimental and numerical studies. In a pure morphological case, when there is no flow, morphological changes in the solid-liquid interface are governed by heat conduction and solute distribution. Under the influence of a convective motion, both heat and solute are redistributed, thereby affecting the intrinsic morphological phenomenon. The overall effect of the convective motion could be either stabilizing or destabilizing. Recent investigations have predicted stabilization by a flow parallel to the interface. In the case of non-parallel flows, e.g., stagnation point flow, Brattkus and Davis have found a new flow-induced morphological instability that occurs at long wavelengths and also consists of waves propagating against the flow. Other studies have addressed the nonlinear aspects (Konstantinos and Brown, Wollkind and Segel)). In contrast to the earlier studies, our present investigation focuses on the effects of the potential flow fields typically encountered in Hele-Shaw cells. Such a Hele-Shaw cell can simulate a gravity-free environment in the sense that buoyancy-driven convection is largely suppressed, and hence negligible. Our interest lies both in analyzing the linear stability of the solidification process in the presence of potential flow fields, as well as in performing high-accuracy nonlinear simulations. Linear stability analysis can be performed for the flow configuration mentioned above. It is observed that a parallel potential flow is stabilizing and gives rise to waves traveling downstream. We have built a highly accurate numerical scheme which is validated at small amplitudes by comparing with the analytically predicted results for the pure morphological case. We have been able to observe nonlinear effects at larger times. Preliminary results for the case when flow is imposed also provide good validation at small amplitudes.

Buddhavarapu, Sudhir S.↗
NASA NTRS2001

Effect of Under-Resolved Grids on High Order Methods

There has been much discussion on verification and validation processes for establishing the credibility of CFD simulations. Since the early 1990s, many of the aeronautical and mechanical engineering related reference journals mandated that any accepted articles in numerical simulations (without known solutions to compared with) need to perform a minimum of one level of grid refinement and time step reduction. Due to the difficulty in analysis, the effect of under-resolved grids and the nonlinear behavior of available spatial discretizations, are scarcely discussed in the literature. Here, an under-resolved numerical simulation is one where the grid spacing being used is too coarse to resolve the smallest physically relevant scales of the chosen continuum governing equations that are of interest to the numerical modeler. With the advent of new developments in fourth-order or higher spatial schemes, it has become common to regard high order schemes as more accurate, reliable and require less grid points. The danger comes when one tries to perform computations with the coarsest grid possible while still hoping to maintain numerical results sufficiently accurate for complex flows, and especially, data-limited problems. On one hand, high order methods when applies to highly coupled multidimensional complex nonlinear problems might have different stability, convergence and reliability behavior than their well studied low order counterparts, especially for nonlinear schemes such as TVD, MUSCL with limiters, ENO, WENO and discrete Galerkin. On the other hand, high order methods involve more operation counts and systematic grid convergence study can be time consuming and prohibitively expansive. At the same time it is difficult to fully understand or categorize the different nonlinear behavior of finite discretizations, especially at the limits of under-resolution when different types of bifurcation phenomena might occur, depending on the combination of grid spacings, time steps, initial conditions and numerical treatments of boundary conditions.

Yee, H. C.↗