Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical Analysis and Computation”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Analysis of internal ablation for the thermal control of aerospace vehicles

A new method of thermal protection for transatmospheric vehicles is introduced. The method involves the combination of radiation, ablation and transpiration cooling. By placing an ablating material behind a fixed-shape, porous outer shield, the effectiveness of transpiration cooling is made possible while retaining the simplicity of a passive mechanism. A simplified one-dimensional approach is used to derive the governing equations. Reduction of these equations to non-dimensional form yields two parameters which characterize the thermal protection effectiveness of the shield and ablator combination for a given trajectory. The non-dimensional equations are solved numerically for a sample trajectory corresponding to glide re-entry. Four typical ablators are tested and compared with results obtained by using the thermal properties of water. For the present level of analysis, the numerical computations adequately support the analytical model.

Camberos, Jose A.↗
NASA NTRS1990

Spectral Analysis Of Linear, Shift-Invariant Interpolants

Method of analysis provides quantitative measure of reconstruction and interpolation performances of linear, shift-invariant interpolants. Criterion of performance based upon mean-square error of difference between sampled and reconstructed functions. Applicable to reconstruction algorithms used in processing of signals and images and to types of interpolants used in numerical analysis, computer-aided design, and computer graphics.

Lansing, Donald L.↗
NASA NTRS1993

An instability associated with a magnetosphere-disk interaction

The evolution of a thin accretion disk surrounding a rapidly rotating magnetosphere is considered. By taking account of the variations of the magnetospheric boundary in response to the conditions at the inner edge of the disk, we find from linear analysis and numerical computation that the accretion disk can become unstable. Mass can be accreted by the central object in a cyclic fashion, with the cycle involving the storage and release of mass in the inner parts of the disk. The physical origin of the instability is associated with the variations of the magnetospheric boundary about corotation. The recurrence time scale of the cycle can vary by several orders of magnitude depending on the details of the conditions at the magnotosphere. The possible applicability of this instability process to the 'rapid burster' MXB 1730-335 is briefly discussed.

Spruit, H. C.↗
NASA NTRS1996

Analysis of condensation on a horizontal cylinder with unknown wall temperature and comparison with the Nusselt model of film condensation

Theoretical analysis and numerical computations are performed to set forth a new model of film condensation on a horizontal cylinder. The model is more general than the well-known Nusselt model of film condensation and is designed to encompass all essential features of the Nusselt model. It is shown that a single parameter, constructed explicitly and without specification of the cylinder wall temperature, determines the degree of departure from the Nusselt model, which assumes a known and uniform wall temperature. It is also known that the Nusselt model is reached for very small, as well as very large, values of this parameter. In both limiting cases the cylinder wall temperature assumes a uniform distribution and the Nusselt model is approached. The maximum deviations between the two models is rather small for cases which are representative of cylinder dimensions, materials and conditions encountered in practice.

Bahrami, Parviz A.↗
NASA NTRS2006

Large Array Channel Capacity in the Presence of Interference

We develop a model for a large array ground receiver system for use in deep-space communications, and analyze the resulting array channel capacity. The model includes effects of array geometry, time-dependent spacecraft orbital trajectory, point and extended interference sources, and elevation-dependent noise and tropospheric channel variations. Channel capacity is expressed as the ratio of determinants of covariance matrices characterizing source, interference, and additive noise, and then reduced to a simpler quadratic form more amenable to analysis and numerical computation. This formulation facilitates inclusion of array and channel characteristics into the model, as well as comparison of optimal, suboptimal, and equivalent single antenna configurations on achievable throughput. Realistic examples of ground array channel capacity calculations are presented, demonstrating the impact of array geometry, planetary interference sources, and array combining algorithm design upon the achievable data throughput.

Vilnrotter, V.↗
NASA NTRS2015

2nd-Order CESE Results For C1.4: Vortex Transport by Uniform Flow

The Conservation Element and Solution Element (CESE) method was used as implemented in the NASA research code ez4d. The CESE method is a time accurate formulation with flux-conservation in both space and time. The method treats the discretized derivatives of space and time identically and while the 2nd-order accurate version was used, high-order versions exist, the 2nd-order accurate version was used. In regards to the ez4d code, it is an unstructured Navier-Stokes solver coded in C++ with serial and parallel versions available. As part of its architecture, ez4d has the capability to utilize multi-thread and Messaging Passage Interface (MPI) for parallel runs.

Space-time CE/SE method↗
NASA NTRS1982

Numerical analysis and FORTRAN program for the computation of the turbulent wakes of turbomachinery rotor blades, isolated airfoils and cascade of airfoils

Turbulent wakes of turbomachinery rotor blades, isolated airfoils, and a cascade of airfoils were investigated both numerically and experimentally. Low subsonic and incompressible wake flows were examined. A finite difference procedure was employed in the numerical analysis utilizing the continuity, momentum, and turbulence closure equations in the rotating, curvilinear, and nonorthogonal coordinate system. A nonorthogonal curvilinear coordinate system was developed to improve the accuracy and efficiency of the numerical calculation. Three turbulence models were employed to obtain closure of the governing equations. The first model was comprised to transport equations for the turbulent kinetic energy and the rate of energy dissipation, and the second and third models were comprised of equations for the rate of turbulent kinetic energy dissipation and Reynolds stresses, respectively. The second model handles the convection and diffusion terms in the Reynolds stress transport equation collectively, while the third model handles them individually. The numerical results demonstrate that the second and third models provide accurate predictions, but the computer time and memory storage can be considerably saved with the second model.

Hah, C.↗
NASA NTRS1993

Plasma heating by collisionless magnetic reconnection - Analysis and computation

Analytic and numerical results on particle acceleration in 2D collisionless magnetic reconnection are presented. The particles are followed until they reach an outgoing flux surface at the same distance from the origin as the starting surface. The magnetic moment is not conserved for particles passing through the unmagnetized region around the X line at the origin. Other particles cross the separatrix without passing near the X line. The magnetic moment of the first class of outgoing particles is randomized, whereas it can be considered for the second class. The analytic model is based upon the observation of the final kinetic energy as a function of the initial conditions. Analytic results are shown to predict a Maxwellian tail for the distribution function in the perpendicular kinetic energy, with this energy much greater than the parallel kinetic energy. Numerical results showing that the predicted tail temperature agrees with the numerically computed temperature to within 10 percent over 4 orders of magnitude in the electric field are presented.

Moses, R. W.↗
NASA NTRS2012

Error Estimates of the Ares I Computed Turbulent Ascent Longitudinal Aerodynamic Analysis

Numerical predictions of the longitudinal aerodynamic characteristics for the Ares I class of vehicles, along with the associated error estimate derived from an iterative convergence grid refinement, are presented. Computational results are based on an unstructured grid, Reynolds-averaged Navier-Stokes analysis. The validity of the approach to compute the associated error estimates, derived from a base grid to an extrapolated infinite-size grid, was first demonstrated on a sub-scaled wind tunnel model at representative ascent flow conditions for which the experimental data existed. Such analysis at the transonic flow conditions revealed a maximum deviation of about 23% between the computed longitudinal aerodynamic coefficients with the base grid and the measured data across the entire roll angles. This maximum deviation from the wind tunnel data was associated with the computed normal force coefficient at the transonic flow condition and was reduced to approximately 16% based on the infinite-size grid. However, all the computed aerodynamic coefficients with the base grid at the supersonic flow conditions showed a maximum deviation of only about 8% with that level being improved to approximately 5% for the infinite-size grid. The results and the error estimates based on the established procedure are also presented for the flight flow conditions.

Abdol-Hamid, Khaled S.↗
NASA NTRS1981

Manufacturing in space: Fluid dynamics numerical analysis

Natural convection in a spherical container with cooling at the center was numerically simulated using the Lockheed-developed General Interpolants Method (GIM) numerical fluid dynamic computer program. The numerical analysis was simplified by assuming axisymmetric flow in the spherical container, with the symmetry axis being a sphere diagonal parallel to the gravity vector. This axisymmetric spherical geometry was intended as an idealization of the proposed Lal/Kroes growing experiments to be performed on board Spacelab. Results were obtained for a range of Rayleigh numbers from 25 to 10,000. For a temperature difference of 10 C from the cooling sting at the center to the container surface, and a gravitional loading of 0.000001 g a computed maximum fluid velocity of about 2.4 x 0.00001 cm/sec was reached after about 250 sec. The computed velocities were found to be approximately proportional to the Rayleigh number over the range of Rayleigh numbers investigated.

Robertson, S. J.↗
NASA NTRS1950

Calculation of the aerodynamic loading of swept and unswept flexible wings of arbitrary stiffness

A method is presented for calculating the aerodynamic loading, the divergence speed, and certain stability derivatives of swept and unswept wings and tail surfaces of arbitrary stiffness. Provision is made for using either stiffness curves and root rotation constants or structural influence coefficients in the analysis. Computing forms, tables of numerical constants required in the analysis, and an illustrative example are included to facilitate calculations by means of the method.

Diederich, Franklin W↗
NASA NTRS1988

Hybrid transfinite element modeling/analysis of nonlinear heat conduction problems involving phase change

The present paper describes the applicability of hybrid transfinite element modeling/analysis formulations for nonlinear heat conduction problems involving phase change. The methodology is based on application of transform approaches and classical Galerkin schemes with finite element formulations to maintain the modeling versatility and numerical features for computational analysis. In addition, in conjunction with the above, the effects due to latent heat are modeled using enthalpy formulations to enable a physically realistic approximation to be dealt computationally for materials exhibiting phase change within a narrow band of temperatures. Pertinent details of the approach and computational scheme adapted are described in technical detail. Numerical test cases of comparative nature are presented to demonstrate the applicability of the proposed formulations for numerical modeling/analysis of nonlinear heat conduction problems involving phase change.

Tamma, Kumar K.↗
NASA NTRS1981

Numerical simulation of natural convection in a spherical container due to cooling at the center (idealization of the Lal/Kroes experiment)

Natural convection in a spherical container with cooling at the center was numerically simulated using a numerical fluid dynamics computer program. The numerical analysis was simplified by assuming axisymmetric flow in the spherical container, with the symmetry axis being a sphere diagonal parallel to the gravity vector. This axisymmetric spherical geometry was intended as an idealization of the proposed Lal/Kroes crystal growing experiment to be performed on Spacelab. Results were obtained for a range of Rayleigh numbers from 25 to 10,000. The computed velocities were found to be approximately proportional to the Rayleigh number over the range of Rayleigh numbers investigated.

Robertson, S. J.↗
NASA NTRS1972

Thermal performance of 625-kg/cu m elastomeric ablative materials on spherically blunted 0.44-radian cones

Spherically blunted 0.44-radian (25 deg) half-angle conical models coated with elastomeric ablative materials were tested in supersonic arc-heated wind tunnels to evaluate performance of the ablators over a range of conditions typical of lifting entry. Four test conditions were combinations of stagnation point-heat transfer rates of 2.3 and 4.5 MW/m2 and stagnation pressures of 20 and 2kN/m2. Afterbody values of heat transfer rate and pressure were 0.05 to 0.20 of stagnation point values. Stagnation enthalpy varied from 4.4 to 25 MJ/kg (1900 to 11000 Btu/lbm) and free-stream Mach number was in a range from 3.5 to 4. Ablative materials retained the spherical nose shape throughout tests at the lower heat transfer level, but receded, assuming a flattened nose shape, during tests at the high heat transfer level. The residue layer that formed on the conical after-body was weak, friable, and extensively cracked. The reference ablative material, which contained phenolic microspheres, generally retained the conical shape on the model afterbody. However, a modified ablator, in which phenolic microspheres were replaced with silica microspheres, deformed and separated from the undegraded material, and thereby produced a very uneven surface. Substrate temperatures and ablator recession were in good agreement with values computed by a numerical analysis.

Champman, A. J.↗