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 865 records · Page 48

Error propagation in the numerical solutions of the differential equations of orbital mechanics

The relationship between the eigenvalues of the linearized differential equations of orbital mechanics and the stability characteristics of numerical methods is presented. It is shown that the Cowell, Encke, and Encke formulation with an independent variable related to the eccentric anomaly all have a real positive eigenvalue when linearized about the initial conditions. The real positive eigenvalue causes an amplification of the error of the solution when used in conjunction with a numerical integration method. In contrast an element formulation has zero eigenvalues and is numerically stable.

Bond, V. R.↗

Global tropospheric chemistry: A plan for action

Prompted by an increasing awareness of the influence of human activity on the chemistry of the global troposphere, a panel was formed to (1) assess the requirement for a global study of the chemistry of the troposphere; (2) develop a scientific strategy for a comprehensive plan taking into account the existing and projected programs of the government; (3) assess the requirements of a global study in terms of theoretical knowledge, numerical modeling, instrumentation, observing platforms, ground-level observational techniques, and other related needs; and (4) outline the appropriate sequence and coordination required to achieve the most effective utilization of available resources. Part 1 presents a coordinated national blueprint for scientific investigations of biogeochemical cycles in the global troposphere. part 2 presents much of the background information of the present knowledge and gaps in the understanding of tropospheric chemical cycles and processes from which the proposed program was developed.

Source record↗

Acoustic emission monitoring crack propagation in single crystal silicon

The feasibility of acoustic emission (AE) monitoring of cracking and crack propagation in Si semiconductor materials was evaluated experimentally. A double torsion load relaxation method was employed wherein the propagation velocity and the AE levels in precracked (but not notched) boron-doped wafers were recorded simultaneously. A numerical model for the critical stress intensity factor (KIC) was used to relate the crack growth velocity, the instantaneous load and the load relaxation rate. All specimens were monitored with acoustic transducers at six points and examined with SEM after failure. The AE levels reached a peak amplitude of 70 dB at a KIC of 0.997 MNm to the -3/2 for cracking in the 111 plane. No AE was detected before the load reached the KIC, indicating that no subcritical crack growth occurs in Si. The results support the use of AE for monitoring crack propagation in crystal Si.

Chen, C. P.↗

Radiation and temperature effects in gallium arsenide, indium phosphide and silicon solar cells

The effects of radiation on performance are determined for both n(+)p and p(+)n GaAs and InP cells and for silicon n(+)p cells. It is found that the radiation resistance of InP is greater than that of both GaAs and Si under 1 MeV electron irradiation. For silicon, the observed decreased radiation resistance with decreased resistivity is attributed to the presence of a radiation induced boron-oxygen defect. Comparison of radiation damage in both p(+)n and n(+)p GaAs cells yields a decreased radiation resistance for the n(+)p cell attributable to increased series resistance, decreased shunt resistance, and relatively greater losses in the cell's p-region. For InP, the n(+)p configuration is found to have greater radiation resistance than the p(+)n cell. The increased loss in this latter cell is attributed to losses in the cell's emitter region. Temperature dependency results are interpreted using a theoretical relation for dVoc/cT which predicts that increased Voc should results in decreased numerical values for dPm/dT. The predicted correlation is observed for GaAs but not for InP a result which is attributed to variations in cell processing.

Weinberg, I.↗

Radiation and temperature effects in gallium arsenide, indium phosphide, and silicon solar cells

The effects of radiation on performance are determined for both n+p and p+n GaAs and InP cells and for silicon n+p cells. It is found that the radiation resistance of InP is greater than that of both GaAs and Si under 1-MeV electron irradiation. For silicon, the observed decreased radiation resistance with decreased resistivity is attributed to the presence of a radiation-induced boron-oxygen defect. Comparison of radiation damage in both p+n and n+p GaAs cells yields a decreased radiation resistance for the n+p cell attributable to increased series resistance, decreased shunt resistance, and relatively greater losses in the cell's p-region. For InP, the n+p configuration is found to have greater radiation resistance than the p+n cell. The increased loss in this latter cell is attributed to losses in the cell's emitter region. Temperature dependency results are interpreted using a theoretical relation for dVoc/dT, which predicts that increased Voc should result in decreased numerical values for dPm/dT. The predicted correlation is observed for GaAs but not for InP, a result which is attributed to variations in cell processing.

Weinberg, I.↗

Quotient-difference type generalizations of the power method and their analysis

The recursion relations that were proposed by W. F. Ford and A. Sidi (Appl. Numer. Math, 4 (1988), pp. 477-489) for implementing vector extrapolation methods are used for devising generalizations of the power method for linear operators. These generalizations are shown to produce approximations to largest eigenvalues of a linear operator under certain conditions. They are similar in form to the quotient-difference algorithm and share similar convergence properties with the latter. These convergence properties also resemble those obtained for the basic LR and QR algorithms. Finally, it is shown that the convergence rate produced by one fo these generalizations is twice as fast for normal operators as it is for nonnormal operators.

Sidi, Avram↗

Acceleration, gamma, and theta guidance for abort landing in a windshear

This paper is concerned with the guidance of abort landing trajectories in a windshear. First, optimal trajectories are determined by minimizing the peak value of the altitude drop. Then, two guidance schemes, approximating the optimal trajectors, are developed: acceleration guidance (based on the relative acceleration) and gamma guidance (based on the absolute path inclination). From numerical experiments, it appears that both the acceleration guidance and the gamma guidance yield trajectories that are close to the optimal trajectory. In addition, a theta guidance scheme (modified constant pitch guidance) is developed that is superior to the constant pitch guidance in terms of the altitude loss and the survival capability in severe windshears.

Miele, A.↗

Reliability analysis of structural ceramic components using a three-parameter Weibull distribution

Described here are nonlinear regression estimators for the three-Weibull distribution. Issues relating to the bias and invariance associated with these estimators are examined numerically using Monte Carlo simulation methods. The estimators were used to extract parameters from sintered silicon nitride failure data. A reliability analysis was performed on a turbopump blade utilizing the three-parameter Weibull distribution and the estimates from the sintered silicon nitride data.

Duffy, Stephen F.↗

An Eulerian/Lagrangian method for computing blade/vortex impingement

A combined Eulerian/Lagrangian approach to calculating helicopter rotor flows with concentrated vortices is described. The method computes a general evolving vorticity distribution without any significant numerical diffusion. Concentrated vortices can be accurately propagated over long distances on relatively coarse grids with cores only several grid cells wide. The method is demonstrated for a blade/vortex impingement case in 2D and 3D where a vortex is cut by a rotor blade, and the results are compared to previous 2D calculations involving a fifth-order Navier-Stokes solver on a finer grid.

Steinhoff, John↗

Three-dimensional boundary element thermal shape sensitivity analysis

A computationally efficient and accurate shape design sensitivity analysis approach for the thermal response of 3D solid objects is presented which utilizes a direct, singular, boundary element analysis formulation. The theoretical formulation for the primary response (the surface temperature and normal heat flux) sensitivities and the secondary response (tangential components of the heat flux, vector, interior point temperature and heat flux vector components) sensitivities are presented. A number of computational issues related to the overall efficiency of implementation of these formulations are discussed. Numerical results are presented to demonstrate the accuracy and efficiency of this approach.

Prasad, K. G.↗

An algebraic turbulence model for three-dimensional viscous flows

An algebraic turbulence model is proposed for use with three-dimensional Navier-Stokes analyses. It incorporates features of both the Baldwin-Lomax and Cebeci-Smith models. The Baldwin-Lomax model uses the maximum of a function f(y) to determine length and velocity scales. An analysis of the Baldwin-Lomax model shows that f(y) can have a spurious maximum close to the wall, causing numerical problems and non-physical results. The proposed model uses integral relations to determine delta(*) u(sub e) and delta used in the Cebeci-Smith mode. It eliminates a constant in the Baldwin-Lomax model and determines the two remaining constants by comparison to the Cebeci-Smith formulation. Pressure gradient effects, a new wake model, and the implementation of these features in a three-dimensional Navier-Stokes code are also described. Results are shown for a flat plate boundary layer, an annular turbine cascade, and endwall heat transfer in a linear turbine cascade. The heat transfer results agree well with experimental data which shows large variations in endwall Stanton number contours with Reynolds number.

Chima, R. V.↗

Laser diagnostics for microgravity droplet studies

Rapid advances have recently been made in numerical simulation of droplet combustion under microgravity conditions, while experimental capabilities remain relatively primitive. Calculations can now provide detailed information on mass and energy transport, complex gas-phase chemistry, multi-component molecular diffusion, surface evaporation and heterogeneous reaction, which provides a clearer picture of both quasi-steady as well as dynamic behavior of droplet combustion. Experiments concerning these phenomena typically result in pictures of the burning droplets, and the data therefrom describe droplet surface regression along with flame and soot shell position. With much more precise, detailed, experimental diagnostics, significant gains could be made on the dynamics and flame structural changes which occur during droplet combustion. Since microgravity experiments become increasingly more expensive as they progress from drop towers and flights to spaceborne experiments, there is a great need to maximize the information content from these experiments. Sophisticated measurements using laser diagnostics on individual droplets and combustion phenomena are now possible. These include measuring flow patterns and temperature fields within droplets, vaporization rates and vaporization enhancement, radical species profiling in flames and gas-phase flow-tagging velocimetry. Although these measurements are sophisticated, they have undergone maturation to the degree where with some development, they are applicable to studies of microgravity droplet combustion. This program beginning in September of 1992, will include a series of measurements in the NASA Learjet, KC-135 and Drop Tower facilities for investigating the range of applicability of these diagnostics while generating and providing fundamental data to ongoing NASA research programs in this area. This program is being conducted in collaboration with other microgravity investigators and is aimed toward supplementing their experimental efforts.

Winter, Michael↗

Reliability analysis of structural ceramic components using a three-parameter Weibull distribution

Described here are nonlinear regression estimators for the three-parameter Weibull distribution. Issues relating to the bias and invariance associated with these estimators are examined numerically using Monte Carlo simulation methods. The estimators were used to extract parameters from sintered silicon nitride failure data. A reliability analysis was performed on a turbopump blade utilizing the three-parameter Weibull distribution and the estimates from the sintered silicon nitride data.

Duffy, Stephen F.↗

Boundary formulations for three-dimensional continuum structural shape sensitivity analysis

The direct, singular, boundary element analysis formulation is shown to provide a basis for a computationally efficient and accurate shape design sensitivity analysis approach for the structural response of three-dimensional solid objects. The theoretical formulation for surface displacement and traction component sensitivities, and all components of the stress tensor is presented along with a formulation for the recovery of displacement and stress components in the interior of the object under consideration. Discussion of computational issues related to the overall efficiency of these formulations is given, along with numerical results to demonstrate the accuracy and efficiency of this approach.

Kane, J. H.↗

An algebraic turbulence model for three-dimensional viscous flows

An algebraic turbulence model is proposed for use with three-dimensional Navier-Stokes analyses. It incorporates features of both the Baldwin-Lomax and Cebeci-Smith models. The Baldwin-Lomax model uses the maximum of a function f(v) to determine length and velocity scales. An analysis of the Baldwin-Lomax model shows that f(v) can have a spurious maximum close to the wall, causing numerical problems and non-physical results. The proposed model uses integral relations to determine delta(*) u(sub e) and delta used in the Cebeci-Smith mode. It eliminates a constant in the Baldwin-Lomax model and determines the two remaining constants by comparison to the Cebeci-Smith formulation. Pressure gradient effects, a new wake model, and the implementation of these features in a three-dimensional Navier-Stokes code are also described. Results are shown for a flat plate boundary layer, an annular turbine cascade, and endwall heat transfer in a linear turbine cascade. The heat transfer results agree well with experimental data which shows large variations in endwall Stanton number contours with Reynolds number.

Chima, R. V.↗

Induced Velocities Near a Lifting Rotor with Nonuniform Disk Loading

A method is given for converting known uniformly loaded rotor induced velocities to correspond with arbitrary axisymmetric nonuniform disk load distributions. Numerical results for two specific distributions are given in chart form. Symmetry relations and relations between radial disk loading and wake velocities are developed. Experimental flow measurements are presented and compared with theory. Reasonable agreement is shown in the forward part of the flow when nonuniform loading is assumed, but far behind the rotor the flow is more like that of a wing.

Heyson, Harry H↗

Observational and numerical studies of extreme frontal scale contraction

The general objective of this effort is to increase understanding of how frontal scale contraction processes may create and sustain intense mesoscale precipitation along intensifying cold fronts. The five-part project (an expansion of the originally proposed two-part project) employed conventional meteorological data, special mesoscale data, remote sensing measurements, and various numerical models. First an idealized hydrostatic modeling study of the scale contraction effects of differential cloud cover on low-level frontal structure and dynamics was completed and published in a peer-reviewed journal. The second objective was to complete and publish the results from a three dimensional numerical model simulation of a cold front in which differential sensible heating related to cloud coverage patterns was apparently crucial in the formation of a severe frontal squall line. The third objective was to use a nonhydrostatic model to examine the nonlinear interactions between the transverse circulation arising from inhomogeneous cloud cover, the adiabatic frontal circulation related to semi-geostrophic forcing, and diabatic effects related to precipitation processes, in the development of a density current-like microstructure at the leading edge of cold fronts. Although the development of a frontal model that could be used to initialize such a primitive equation model was begun, we decided to focus our efforts instead on a project that could be successfully completed in this short time, due to the lack of prospects for continued NASA funding beyond this first year (our proposal was not accepted for future funding). Thus, a fourth task was added, which was to use the nonhydrostatic model to test tentative hypotheses developed from the most detailed observations ever obtained on a density current (primarily sodar and wind profiler data). These simulations were successfully completed, the findings were reported at a scientific conference, and the results have recently been submitted to a peer-reviewed journal. The fifth objective was to complete the analysis of data collected during the Cooperative Oklahoma Profiler Studies (COPS-91) field project, which was supported by NASA. The analysis of the mesoscale surface and sounding data, Doppler radar imagery, and other remote sensing data from multi frequency wind profiler, microwave radiometer, and the Radio Acoustic Sounding System has been completed. This study is a unique investigation of processes that caused the contraction of a cold front to a microscale zone exhibiting an undular bore-like structure. Results were reported at a scientific conference and are being prepared for publication. In summary, considerable progress has been achieved under NASA funding in furthering our understanding of frontal scale contraction and density current - gravity wave interaction processes, and in utilizing models and remotely sensed data in such studies.

Koch, Steven E.↗

Non-Coalescence Effects in Microgravity

Non-coalescence of two bodies of the same liquid and the suppression of contact between liquid drops and solid surfaces is being studied through a pair of parallel investigations being conducted at the Georgia Institute of Technology and the Microgravity Research and Support (MARS) Center in Naples, Italy. Both non-coalescence and contact suppression are achieved by exploiting the mechanism of thermocapillary convection to drive a lubricating film of surrounding gas (air) into the space between the two liquid free surfaces (non-coalescence) or between the drop free surface and the solid (contact suppression). Earlier experiments performed included flow-visualization experiments in both axisymmetric and (nearly) two-dimensional geometries and quantitative measurements of film thickness in the contact-suppression case in both geometries. Work done in the second year has focused on obtaining quantitative results relating to the effects of variable air pressure, development of analytical and numerical models of non-coalescing droplets and to pursuing potential applications of these self-lubricated systems.

Neitzel, G. Paul↗