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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

Space Weathering of Olivine: Samples, Experiments and Modeling

Olivine is a major constituent of chondritic bodies and its response to space weathering processes likely dominates the optical properties of asteroid regoliths (e.g. S- and many C-type asteroids). Analyses of olivine in returned samples and laboratory experiments provide details and insights regarding the mechanisms and rates of space weathering. Analyses of olivine grains from lunar soils and asteroid Itokawa reveal that they display solar wind damaged rims that are typically not amorphized despite long surface exposure ages, which are inferred from solar flare track densities (up to 10 (sup 7 y)). The olivine damaged rim width rapidly approaches approximately 120 nm in approximately 10 (sup 6 y) and then reaches steady-state with longer exposure times. The damaged rims are nanocrystalline with high dislocation densities, but crystalline order exists up to the outermost exposed surface. Sparse nanophase Fe metal inclusions occur in the damaged rims and are believed to be produced during irradiation through preferential sputtering of oxygen from the rims. The observed space weathering effects in lunar and Itokawa olivine grains are difficult to reconcile with laboratory irradiation studies and our numerical models that indicate that olivine surfaces should readily blister and amorphize on relatively short time scales (less than 10 (sup 3 y)). These results suggest that it is not just the ion fluence alone, but other variable, the ion flux that controls the type and extent of irradiation damage that develops in olivine. This flux dependence argues for caution in extrapolating between high flux laboratory experiments and the natural case. Additional measurements, experiments, and modeling are required to resolve the discrepancies among the observations and calculations involving solar wind processing of olivine.

Keller, L. P.↗
NASA NTRS1986

Application of shuttle imaging radar to geologic mapping

Images from the Shuttle Imaging Radar - B (SIR-B) experiment covering the area of the Panamint Mountains, Death Valley, California, were examined in the field and in the laboratory to determine their usefulness as aids for geologic mapping. The covered area includes the region around Wildrose Canyon where rocks ranging in age from Precambrian to Cenozoic form a moderately rugged portion of the Panamint Mountains, including sharp ridges, broad alluviated upland valleys, and fault-bounded grabens. The results of the study indicate that the available SIR-B images of this area primarily illustrate variations in topography, except in the broadly alluviated areas of Panamint Valley and Death Valley where deposits of differing reflectivity can be recognized. Within the mountainous portion of the region, three textures can be discerned, each representing a different mode of topographic expression related to the erosion characteristics of the underlying bedrock. Regions of Precambrian bedrock have smooth slopes and sharp ridges with a low density of gullies. Tertiary monolithologic breccias have smooth, steep slopes with an intermediate density of gullies with rounded ridges. Tertiary fanglomerates have steep rugged slopes with numerous steep-sided gullies and knife-sharp ridges. The three topographic types reflect the consistancy and relative susceptibility to erosion of the bedrock; the three types can readily be recognized on topographic maps. At present, it has not been possible to distinguish on the SIR-B image of the mountainous terrain the type of bedrock, independent of the topographic expression.

Labotka, T. C.↗
NASA NTRS1989

Nonlinear strain-displacement relations and flexible multibody dynamics

Dynamics of chains of flexible bodies undergoing large rigid body motions, but small elastic deflections are considered. The role of nonlinear strain-displacement relations in the development of the motion equations correct to first order in elastic deflections is investigated. The general form of these equations linearized only in the small elastic deflections is presented, and the relative significance of various nonlinear terms is studied both analytically and through the use of the numerical simulations. Numerical simulations are performed for a two link chain constrained to move in the plane, subject to hinge torques. Each link is modeled as a thin beam. Slew maneuver simulation results are compared for models with and without properly modeled kinematics of deformation. The goal of this case study is to quantify the importance of the terms in the equations of motion which arise from the inclusion of nonlinear strain-displacement relations. It is concluded that unless the consistently linearized equations in elastic deflections and speeds are available and necessary, the inconsistently (prematurely) linearized equations should be replaced in all cases by ruthlessly linearized equations: equations in which all nonlinear terms involving the elastic deflections and speeds are ignored.

Padilla, Carlos E.↗
NASA NTRS1985

Periodic orbits of the general three-body problem for the sun-Jupiter-Saturn system

Two families of symmetric periodic orbits of the planar, general, three-body problem are presented. The masses of the three bodies include ratios equal to the sun-Jupiter-Saturn system and the periods of the orbits of Jupiter and Saturn are in a 2:5 resonance. The (linear) stability of the orbits are studied in relation to eccentricity and mass variations. The generation of the two families of periodic orbits follows a systematic approach and employs (numerical) continuation from periodic orbits of the first and second kind in the circular restricted problem to the elliptic restricted problem and from the circular and elliptic problems to the general problem through bifurcation phenomena relating the three dynamical systems. The approach also provides insight into the evolutionary process of periodic orbits continued from the restricted problems to the general problem.

Kwok, J. H.↗
NASA NTRS1997

Turbulent Premixed Flame Propagation in Microgravity

A facility in which turbulent Couette flow could be generated in a microgravity environment was designed and built. To fit into the NASA Lewis drop tower the device had to be very compact. This means that edge effects and flow re-circulation were expected to affect the flow. The flow was thoroughly investigated using LDV and was found to be largely two dimensional away from the edges with constant turbulence intensities in the core. Slight flow asymmetries are introduced by the non symmetric re-circulation of the fluid outside the test region. Belt flutter problems were remedied by adding a pair of guide plates to the belt. In general, the flow field was found to be quite similar to previously investigated Couette flows. However, turbulence levels and associated shear stresses were higher. This is probably due to the confined re-circulation zone reintroducing turbulence into the test section. An estimate of the length scales in the flow showed that the measurements were able to resolve nearly all the length scales of interest. Using a new LES method for subgrid combustion it has been demonstrated that the new procedure is computational feasible even on workstation type environment. It is found that this model is capable of capturing the propagation of the premixed names by resolving the flame in the LES grid within 2-3 grid points. In contrast, conventional LES results in numerical smearing of the flame and completely inaccurate estimate of the turbulent propagation speed. Preliminary study suggests that there is observable effect of buoyancy in the 1g environment suggesting the need for microgravity experiments of the upcoming experimental combustion studies. With the cold flow properties characterized, an identical hot flow facility is under construction. It is assumed that the turbulence properties ahead of the flame in this new device will closely match the results obtained here. This is required since the hot facility will not enable LDV measurements. The reacting flow facility is also being constructed with planned drop tower experiments in mind. Therefore, issues related to safety and structural integrity are being take into account. Further development of the numerical model will also be carried out to include finite-rate kinetics for representative premixed cases. More detail analysis of the flame structure and propagation nature will be investigated. Simulations will also be compared to the flame properties observed in the experiments.

Menon, S.↗
NASA NTRS2016

The Impacts of Dry Dynamic Cores on Asymmetric Hurricane Intensification

The fundamental pathways for tropical cyclone (TC) intensification are explored by considering axisymmetric and asymmetric impulsive thermal perturbations to balanced, TC-like vortices using the dynamic cores of three different nonlinear numerical models. Attempts at reproducing the results of previous work, which used the community WRF Model, revealed a discrepancy with the impacts of purely asymmetric thermal forcing. The current study finds that thermal asymmetries can have an important, largely positive role on the vortex intensification, whereas other studies find that asymmetric impacts are negligible. Analysis of the spectral energetics of each numerical model indicates that the vortex response to asymmetric thermal perturbations is significantly damped in WRF relative to the other models. Spectral kinetic energy budgets show that this anomalous damping is primarily due to the increased removal of kinetic energy from the vertical divergence of the vertical pressure flux, which is related to the flux of inertia-gravity wave energy. The increased kinetic energy in the other two models is shown to originate around the scales of the heating and propagate upscale with time from nonlinear effects. For very large thermal amplitudes (50 K), the anomalous removal of kinetic energy due to inertia-gravity wave activity is much smaller, resulting in good agreement between models. The results of this paper indicate that the numerical treatment of small-scale processes that project strongly onto inertia-gravity wave energy can lead to significant differences in asymmetric TC intensification. Sensitivity tests with different time integration schemes suggest that diffusion entering into the implicit solution procedure is partly responsible for the anomalous damping of energy.

Guimond, Stephen R.↗
NASA NTRS1988

Design of an expert system for estimating the cost of new knowledge in high energy astrophysics

The High Energy Astrophysics Costing Tool (HEACT) is a combined Expert System/Numerical Simulation to evaluate sensors and experiments proposed for space platforms in order to determine their expected performance and to assess the amount and quality of scientific information likely to be gathered by such experiments. The end product of this tool will be both a cost estimate for the experiment package as deployed and also an estimate of the amount of new scientific knowledge (both new data and reductions in the uncertainties of previous measurements) which would result from the actual deployment and operation of that experiment package. The numerical simulation part of HEACT will contain components to: (1) calculate the high energy signatures of both actual and proposed classes of astrophysical objects, and (2) model the response and resolution of the candidate sensors for observing those objects. The rule based part of the HEACT will be based on knowledge in the areas of experiment design, instrument selection, and system integration developed from experience with existing high energy astrophysics observatories, both orbiting and lofted by rockets or balloons. This paper presents mainly results from the design of the rule-based part of HEACT. It focuses on the structure of the rules and the structures used to represent the knowledge related to the design and integration of such experiments. Methods for integrating the results from numerical simulations into the expert system will also be discussed.

Bosworth, Edward L., Jr.↗
NASA NTRS1980

Investigation of transonic flow in a cascade using an adaptive mesh

The solution of two-dimensional full potential equation for the analysis of steady transonic flow through cascades is investigated. Finite element method is employed in the analysis. Accuracy and efficiency of the obtained numerical solutions are discussed in terms of the employed computational grid. Accurate modeling of subsonic and supersonic flow regions together with the shock is discussed. The choice of artificial viscosity and relaxation factors are examined and related to the design of a computational grid. Shock capturing and shock fitting procedures are compared for improved accuracy and efficiency. Numerical results include cascades of Gostelow and NACA 0012 airfoils.

Ecer, A.↗
NASA NTRS1981

Interaction of upstream flow distortions with high Mach number cascades

Features of the interaction of flow distortions, such as gusts and wakes with blade rows of advance type fans and compressors having high tip Mach numbers are modeled. A typical disturbance was assumed to have harmonic time dependence and was described, at a far upstream location, in three orthogonal spatial coordinates by a double Fourier series. It was convected at supersonic relative to a linear cascade described as an unrolled annulus. Conditions were selected so that the component of this velocity parallel to the axis of the turbomachine was subsonic, permitting interaction between blades through the upstream as well as downstream flow media. A strong, nearly normal shock was considered in the blade passages which was allowed curvature and displacement. The flows before and after the shock were linearized relative to uniform mean velocities in their respective regions. Solution of the descriptive equations was by adaption of the Wiener-Hopf technique, enabling a determination of distortion patterns through and downstream of the cascade as well as pressure distributions on the blade and surfaces. Details of interaction of the disturbance with the in-passage shock were discussed. Infuences of amplitude, wave length, and phase of the disturbance on lifts and moments of cascade configurations are presented. Numerical results are clarified by reference to an especially orderly pattern of upstream vertical motion in relation to the cascade parameters.

Englert, G. W.↗
NASA NTRS1997

Numerical Flow Visualization in Basic- and Hyper-Cluster Spheres

Packed spherical particle beds have wide applications throughout the process industry and are usually analyzed using an appropriate combination of laminar and turbulent flows with empirically derived coefficients of which the Ergun (1952) relation is probably the best known. The 3-D complexity of the void distribution within the beds has precluded detailed studies of sphere clusters. Numerical modeling and flow vector visualization for the basic tetra- and hexa-sphere clusters and two hyper-sphere clusters are presented at two Reynolds numbers, 400 and 1200. Cutting planes are used to enable visualization of the complex flows generated within the sphere clusters and are discussed herein. The boundary conditions and flow fields for the simple clusters are also compared to the hyper-clusters with larger variations found for hexa-clusters.

Hendricks, R. C.↗
NASA NTRS1982

Cloud cover sensor

An apparatus is described which provides a numerical indication of the cloudiness at a particular time of a day. The apparatus includes a frame holding several light sensors such as photovoltaic cells, with a direct sensor mounted to directly face the Sun and indirect sensors mounted to face different portions of the sky not containing the Sun. A light shield shields the direct sensor from most of the sky except a small portion containing the Sun, and also shields each of the indirect sensors from direct sunlight. The relative values of the outputs from the direct and indirect sensors, enables the generation of a numerical indication of the degree of cloudiness at a particular time of day.

Laue, E. G.↗
NASA NTRS2000

Generation and Radiation of Acoustic Waves from a 2-D Shear Layer using the CE/SE Method

In the present work, the generation and radiation of acoustic waves from a 2-D shear layer problem is considered. An acoustic source inside of a 2-D jet excites an instability wave in the shear layer, resulting in sound Mach radiation. The numerical solution is obtained by solving the Euler equations using the space time conservation element and solution element (CE/SE) method. Linearization is achieved through choosing a small acoustic source amplitude. The Euler equations are nondimensionalized as instructed in the problem statement. All other conditions are the same except that the Crocco's relation has a slightly different form. In the following, after a brief sketch of the CE/SE method, the numerical results for this problem are presented.

Loh, Ching Y.↗
NASA NTRS1982

Numerical grid generation; Proceedings of the Symposium on Numerical Generation of Curvilinear Coordinate Systems and Their Use in the Numerical Solution of Partial Differential Equations, Nashville, TN, April 13-16, 1982

General curvilinear coordinate systems are considered along with the error induced by coordinate systems, basic differential models for coordinate generation, elliptic grid generation, conformal grid generation, algebraic grid generation, orthogonal grid generation, patched coordinate systems, and solid mechanics applications of boundary fitted coordinate systems. Attention is given to coordinate system control and adaptive meshes, the application of body conforming curvilinear grids for finite difference solution of external flow, the use of solution adaptive grids in solving partial differential equations, adaptive gridding for finite difference solutions to heat and mass transfer problems, and the application of curvilinear coordinate generation techniques to the computation of internal flows. Other topics explored are related to the solution of nonlinear water wave problems using boundary-fitted coordinate systems, the numerical modeling of estuarine hydrodynamics on a boundary-fitted coordinate system, and conformal grid generation for multielement airfoils.

Thompson, J. F.↗
NASA NTRS1985

Electromagnetic electron beam instabilities - Hot, isotropic beams

This paper considers the linear theory of electromagnetic instabilities driven by an electron beam in a homogeneous, nonrelativistic, Vlasov plasma. The beam is relatively hot, isotropic in its own frame, and streams parallel or antiparallel to a magnetic field B. Numerical solutions of the full dispersion equation for propagation parallel or antiparallel to B are presented, and the linear properties of the whistler heat flux and electron beam firehose instabilities are exhibited and compared. Under a broad range of parmameters the former mode has the lower beam speed threshold, and the larger maximum growth rate. In addition, it is demonstrated that, for a sufficiently large relative beam density, relative beam temperature, and plasma beta the whistler heat flux instability has a much lower beam speed threshold than the electrostatic electron beam instability. The application of these instabilities to first-order Fermi acceleration of electrons at space plasma shocks is discussed.

Gary, S. P.↗
NASA NTRS1992

Vorticity generation mechanisms in parallel injection schemes for supersonic mixing

A numerical study is reported of the three-dimensional nonreacting supersonic flow field produced by three parallel fuel injection schemes. Such injections are being considered as a means of enhancing fuel mixing in future Scramjet engine designs. The strength of vortex structures produced by ramp injectors and by rocket/jet interactions are quantified and their relative effectiveness in enhancing the fuel mixing process is addressed. An experimental setup for validation of the above numerical results has been constructed and preliminary results are presented. For the flow field parameters chosen, the ramp generated vorticity is found to be considerably larger than that generated by the shock. The unsteady recirculating flow field in the wake of the injector appears also to give rise to important additional vorticity and mixing.

Donohue, James M.↗
NASA NTRS1988

Studies on pressure response of gas bubbles contributions of condensed droplets in bubbles generated by a uniform nucleation

The response of a tiny gas bubble under reduced pressure is investigated in its relation to cavitation. Equations of motion are formulated for gas mixtures inside the bubble and numerical calculations performed for several examples. The conclusions are as follows: (1) at the onset of bubble growth, the gas mixture inside it adiabatically expands and the temperature decreases. Condensed droplets appear inside the gas mixture due to a uniform nucleation and the temperature recovers, thus the motion of the bubble is apparently isothermal; (2) the evaporation and condensation coefficient largely affects bubble motions (maximum radius, period and rate of attenuation of the bubble oscillation) including the uniform contraction; (3) the oscillation period of the bubble is longer as the equilibrium bubble radius is larger when the surrounding pressure decreases stepwise. In this circumstance the temperature inside the bubble is kept constant due to condensation evaporation phenomena and is nearly isothermal; and (4) when the surrounding pressure decreases in a stepwise fashion, the critical pressure bubble radius relation becomes closer to that for the isothermal process if the bubble radius is larger than 8 microns.

Matsumoto, Y.↗