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.

SourceResearch connectionFrom yearThrough year

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

At least 703 records · Page 39

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

Towards Robust and Accurate Implicit Gradient Methods for Second- and Third-Order Nodal-Gradient Cell-Centered Finite-Volume Discretizations on Tetrahedral Grids

In this paper, we introduce implicit gradient methods as alternatives to conventional least-squares gradient methods for second- and third-order nodal-gradient cell-centered finite-volume discretizations, where solutions are stored at cells but gradients are stored at nodes. Because of the unique configuration of solutions and gradients, implicit gradient systems developed for the node-centered edge-based discretization method can be directly applied once the numerical solutions are interpolated from cells to nodes with sufficient accuracy. The resulting defect-correction solver can be loosely coupled with a flow-equation solver, and at convergence, solutions and gradients that satisfy the corresponding residual equations are obtained. Each iteration is relatively cheap compared with least-squares methods involving hundreds of neighbors. Numerical results are presented for accuracy verification studies and some simple but realistic flow problems.

Computational Fluid Dynamics↗

NASA NTRS1991

NASA DC-8 airborne research laboratory

Since the summer of 1987, NASA Ames Research Center has been operating a DC-8 equipped with CFM 56 engines as a flying research laboratory. In this relatively short time, the DC-8, with its tremendous capabilities, has made significant contributions to numerous scientific fields. Capable of staying aloft for over 12 hours, the DC-8 has flown directly over both the North and South Poles, gathering data relating to the ozone hole. Operating from a few thousand feet to over 40,000 feet above sea level the interchangeable payload capability of the DC-8 has made it a versatile scientific tool. The DC-8 also plays a vital role in the development of new satellite-borne sensors as very often those sensors are test-flown on the DC-8 before they are launched into space. The tremendous range and instrument carrying capability make the DC-8 an ideal flying laboratory. A few of the programs the DC-8 has participated in as well as a sampling of the instruments carried are outlined.

Degreef, Leo H.↗

NASA NTRS1958

On flow of electrically conducting fluids over a flat plate in the presence of a transverse magnetic field

The use of a magnetic field to control the motion of electrically conducting fluids is studied. The incompressible boundary-layer solutions are found for flow over a flat plate when the magnetic field is fixed relative to the plate or to the fluid. The equations are integrated numerically for the effect of the transverse magnetic field on the velocity and temperature profiles, and hence, the skin friction and rate of heat transfer. It is concluded that the skin friction and the heat-transfer rate are reduced when the transverse magnetic field is fixed relative to the plate and increased when fixed relative to the fluid. The total drag is increased in all of the areas.

Rossow, Vernon J↗

NASA NTRS2014

On the Computation of High Order Rys Quadrature Weights and Nodes

Since its introduction in 1976, the Rys Quadrature method has proven a very attractive method for evaluating electron repulsion integrals for calculations using Gaussian type orbitals. Since then, there have been considerable refinements of the method, but at it's core, Gaussian weights and nodes are used to exactly evaluate using a numerical approach to the transform integral. One of the powers of the Rys Quadrature method is the relative ease in evaluating integrals involving functions of high angular momentum. In this work we report on the complete resolution of these numerical difficulties, and we have easily computed accurate quadrature weights and nodes up to order 101. All calculations were carried out using 128-bit precision.

Computation↗

NASA NTRS1972

The relation between momentum transfer and capture and total scattering cross sections for ion-dipole collisions

Numerical values of momentum transfer cross sections sigma sub m for ion-dipole collisions are compared with the corresponding capture cross sections sigma sub c as a function of ion velocity and rotational temperature. For values of dipole moment mu from 1 to 4 Debyes the sigma sub m/sigma sub c ratio is in the range 1.2 to 2.0 (roughly). This is in contrast to the simple relation for Langevin collisions where sigma sub m/sigma sub c is equal to or approximately 1.10 independent of polarizability of the target atom. At low temperatures, the momentum transfer cross sections can be as large as 2000 A squared but they are only about 15 to 30 percent of the total scattering cross sections sigma sub S.

Dugan, J. V., Jr.↗

NASA NTRS2001

Numerical Simulation of Delamination Growth in Composite Materials

The use of decohesion elements for the simulation of delamination in composite materials is reviewed. The test methods available to measure the interfacial fracture toughness used in the formulation of decohesion elements are described initially. After a brief presentation of the virtual crack closure technique, the technique most widely used to simulate delamination growth, the formulation of interfacial decohesion elements is described. Problems related with decohesion element constitutive equations, mixed-mode crack growth, element numerical integration and solution procedures are discussed. Based on these investigations, it is concluded that the use of interfacial decohesion elements is a promising technique that avoids the need for a pre-existing crack and pre-defined crack paths, and that these elements can be used to simulate both delamination onset and growth.

Camanho, P. P.↗

NASA NTRS1989

A numerical method for predicting hypersonic flowfields

The flow about a body traveling at hypersonic speed is energetic enough to cause the atmospheric gases to chemically react and reach states in thermal nonequilibrium. The prediction of hypersonic flowfields requires a numerical method capable of solving the conservation equations of fluid flow, the chemical rate equations for specie formation and dissociation, and the transfer of energy relations between translational and vibrational temperature states. Because the number of equations to be solved is large, the numerical method should also be as efficient as possible. The proposed paper presents a fully implicit method that fully couples the solution of the fluid flow equations with the gas physics and chemistry relations. The method flux splits the inviscid flow terms, central differences of the viscous terms, preserves element conservation in the strong chemistry source terms, and solves the resulting block matrix equation by Gauss Seidel line relaxation.

Maccormack, Robert W.↗

NASA NTRS1983

The evolution of Tollmien-Schlichting waves near a leading edge. II - Numerical determination of amplitudes

In the first part of this investigation, Goldstein (1983) has shown that the amplitude of the spatially growing Tollmien-Schlichting wave generated by a time-harmonic free-stream disturbance is related to the coefficient multiplying the lowest-order asymptotic eigensolution of the unsteady boundary-layer equation. In the present study, a numerical solution of the unsteady boundary-layer equation is used to relate the amplitude of the asymptotic eigensolution, and consequently of the Tollmien-Schlichting wave, to that of the imposed free-stream disturbance for the special case of a uniformly pulsating stream. It is pointed out that the ideas of this study can be extended to other, more complex bodies and free-stream oscillations.

Goldstein, M. E.↗

NASA NTRS1972

Analysis of supersonic combustion flow fields with embedded subsonic regions

The viscous characteristic analysis for supersonic chemically reacting flows was extended to include provisions for analyzing embedded subsonic regions. The numerical method developed to analyze this mixed subsonic-supersonic flow fields is described. The boundary conditions are discussed related to the supersonic-subsonic and subsonic-supersonic transition, as well as a heuristic description of several other numerical schemes for analyzing this problem. An analysis of shock waves generated either by pressure mismatch between the injected fluid and surrounding flow or by chemical heat release is also described.

Dash, S.↗

NASA NTRS1979

General Potential Theory of Arbitrary Wing Sections

The problem of determining the two dimensional potential flow around wing sections of any shape is examined. The problem is condensed into the compact form of an integral equation capable of yielding numerical solutions by a direct process. An attempt is made to analyze and coordinate the results of earlier studies relating to properties of wing sections. The existing approximate theory of thin wing sections and the Joukowski theory with its numerous generalizations are reduced to special cases of the general theory of arbitrary sections, permitting a clearer perspective of the entire field. The method which permits the determination of the velocity at any point of an arbitrary section and the associated lift and moments is described. The method is also discussed in terms for developing new shapes of preassigned aerodynamical properties.

Theodorsen, T.↗

NASA NTRS1976

The axisymmetric jet counterflow problem

A numerical method for computing the flow field of an underexpanded axisymmetric supersonic jet flowing counter to a uniform supersonic stream is described. The strong shock waves and shear layers which occur in the flow are treated as normal and tangential discontinuities bounding regions of inviscid flow. Within the regions, the unsteady gasdynamic equations are solved in integral form on a moving computational mesh. As the numerical flow field evolves, the mesh is deformed according to the jump relations applied at the discontinuities, thus determining their positions as part of the solution. Also presented are computed results obtained for a wide range of jet total pressure. The standoff distances of the bow shock and jet shock are in good agreement with those measured experimentally, thus demonstrating the accuracy of the numerical method.

Schiff, L. B.↗

NASA NTRS1993

Numerical modeling of the energy storage and release in solar flares

This paper reports on investigation of the photospheric magnetic field-line footpoint motion (usually referred to as shear motion) and magnetic flux emerging from below the surface in relation to energy storage in a solar flare. These causality relationships are demonstrated by using numerical magnetohydrodynamic simulations. From these results, one may conclude that the energy stored in solar flares is in the form of currents. The dynamic process through which these currents reach a critical value is discussed as well as how these currents lead to energy release, such as the explosive events of solar flares.

Wu, S. T.↗

NASA NTRS2017

Validation of Quasi-Invariant Ice Cloud Radiative Quantities with MODIS Satellite-Based Cloud Property Retrievals

Similarity relations applied to ice cloud radiance calculations are theoretically analyzed and numerically validated. If t(1v) and t(1vg) are conserved where t is optical thickness, v the single-scattering albedo, and g the asymmetry factor, it is possible that substantially different phase functions may give rise to similar radiances in both conservative and non-conservative scattering cases, particularly in the case of large optical thicknesses. In addition to theoretical analysis, this study uses operational ice cloud optical thickness retrievals from the Moderate Resolution Imaging Spectroradiometer (MODIS) Level 2 Collection5 (C5) and Collection 6 (C6) cloud property products to verify radiative similarity relations. It is found that, if the MODIS C5 and C6 ice cloud optical thickness values are multiplied by their respective (1wg)factors, the resultant products referred to as the effective optical thicknesses become similar with their ratio values around unity. Furthermore, the ratios of the C5 and C6 ice cloud effective optical thicknesses display an angular variation pattern similar to that of the corresponding ice cloud phase function ratios. The MODIS C5 and C6 values of ice cloud similarity parameter, defined as [(1w)(1(exp. 1/2)wg)]12, also tend to be similar.

Radiative transfer↗

NASA NTRS2010

Roles of Engineering Correlations in Hypersonic Entry Boundary Layer Transition Prediction

Efforts to design and operate hypersonic entry vehicles are constrained by many considerations that involve all aspects of an entry vehicle system. One of the more significant physical phenomenon that affect entry trajectory and thermal protection system design is the occurrence of boundary layer transition from a laminar to turbulent state. During the Space Shuttle Return To Flight activity following the loss of Columbia and her crew of seven, NASA's entry aerothermodynamics community implemented an engineering correlation based framework for the prediction of boundary layer transition on the Orbiter. The methodology for this implementation relies upon the framework of correlation techniques that have been in use for several decades. What makes the Orbiter boundary layer transition correlation implementation unique is that a statistically significant data set was acquired in multiple ground test facilities, flight data exists to assist in establishing a better correlation and the framework was founded upon state of the art chemical nonequilibrium Navier Stokes flow field simulations. The basic tenets that guided the formulation and implementation of the Orbiter Return To Flight boundary layer transition prediction capability will be reviewed as a recommended format for future empirical correlation efforts. The validity of this approach has since been demonstrated by very favorable comparison of recent entry flight testing performed with the Orbiter Discovery, which will be graphically summarized. These flight data can provide a means to validate discrete protuberance engineering correlation approaches as well as high fidelity prediction methods to higher confidence. The results of these Orbiter engineering and flight test activities only serve to reinforce the essential role that engineering correlations currently exercise in the design and operation of entry vehicles. The framework of information-related to the Orbiter empirical boundary layer transition prediction capability will be utilized to establish a fresh perspective on this role, to illustrate how quantitative statistical evaluations of empirical correlations can and should be used to assess accuracy and to discuss what the authors' perceive as a recent heightened interest in the application of high fidelity numerical modeling of boundary layer transition. Concrete results will also be developed related to empirical boundary layer transition onset correlations. This will include assessment of the discrete protuberance boundary layer transition onset data assembled for the Orbiter configuration during post-Columbia Return To Flight. Assessment of these data will conclude that momentum thickness Reynolds number based correlations have superior coefficients and uncertainty in comparison to roughness height based Reynolds numbers, aka Re(sub k) or Re(sub kk). In addition, linear regression results from roughness height Reynolds number based correlations will be evaluated, leading to a hypothesis that non-continuum effects play a role in the processes associated with incipient boundary layer transition on discrete protuberances.

Campbell, Charles H.↗

NASA NTRS1999

Validity of the Relation Between Spontaneous and Stimulated Emissions in Semiconductors

The Einstein relation between spontaneous emission and absorption was originally derived for a system consists of a two-state subsystem representing matter and harmonic fields representing radiation. The derivation is based on the detailed balance between these two subsystems under thermal equilibrium. The relationship was later investigated in connection with the interactions between radiation field and solids or semiconductors. The simple derivation dose not hold for semiconductors in general. In certain limiting cases, simple relation was obtained. The validity of this relation is important not only because of its fundamental role connecting two of the most fundamental optical processes in semiconductors, but mostly also because of its wide use as a practical method to measure the optical gain of a semiconductor. The validity of this relation for semiconductors has been an issue of controversial for some time. In this paper we numerically examine the validity of this relationship for several different lineshapes including Lorentzian, Gaussian, Sech, and a convoluted double Lorentzians (CDL). We find out that at relatively low density above transparency level, all first three lineshapes violate the Einstein relation. The relation is approximately valid at high density. At very high density, the validity of the Einstein relation holds well for all three lineshapes. The reason behind this observation is explained. The CDL lineshape has been shown analytically to obey the Einstein relationship previously. We show that for a 2D semiconductor with parabolic bands, the CDL lineshape can be integrated analytically. This analytic lineshape is compared with a simple Lorentzian lineshape.

Ning, Cun-Zheng↗

NASA NTRS2005

Impact into Coarse Grained Spheres

Several experimental studies [1,2,3] indicate that differences in the grain size of the target relative to the projectile could influence the cratering process. Impacts into coarse sand grains of size comparable to the projectile show some discrepancies with existing relationships for crater growth [e.g. 4]. Similarly, targets of ne grained, uniform in diameter glass spheres show differences in crater depth, transient crater diameter, and volume of ejecta excavated as a function of grain size [2,3]. The purpose of this work is to continue investigating how the relative grain size may influence early time coupling between a projectile and target, with implications for subsequent ejecta excavation and crater growth. In previous efforts we used numerical techniques to focus on the propagation of shock waves in coarse, granular media emphasizing the influence of relative grain size on crater growth, ejecta production, cratering efficiency, target strength, and crater shape [5,6,7]. In this study, we use experimental techniques - in part as a reality check for the numerical studies - to report on how coarse grained targets might influence ejecta excavation and crater shape. This body of work possesses important implications for ejecta excavation and cratering efficiency on asteroids that may possess rubble pile-like structures, and on planets that may possess either pre-fractured surfaces or large-scale heterogeneities in shock impedance.

Barnouin-Jha, O. S.↗

NASA NTRS1973

Analysis of a two-dimensional type 6 shock-interference pattern using a perfect-gas code and a real-gas code

Numerical codes were developed to calculate the two dimensional flow field which results when supersonic flow encounters double wedge configurations whose angles are such that a type 4 pattern occurs. The flow field model included the shock interaction phenomena for a delta wing orbiter. Two numerical codes were developed, one which used the perfect gas relations and a second which incorporated a Mollier table to define equilibrium air properties. The two codes were used to generate theoretical surface pressure and heat transfer distributions for velocities from 3,821 feet per second to an entry condition of 25,000 feet per second.

Bertin, J. J.↗

NASA NTRS1973

Computational considerations in application of the finite element method for analysis of unsteady flow around airfoils.

Difficulties related to the application of the finite element method for numerical solution of the Navier-Stokes equation are discussed. The unsteady flow around a stationary and oscillating NACA 0012 airfoil was analyzed. In an effort to establish the range of applicability of the developed numerical method, the accuracy and stability problems in representing unsteady flow patterns and determining pressure distribution around the airfoil were investigated. Error involved in determining the velocity field and in the numerical integration of the discretized equations were analyzed in terms of geometry of the finite element gridwork, boundary conditions for stationary and oscillating airfoils, angles of attack and flow conditions.

Bratanow, T.↗