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 433 records · Page 24

Recent advances in numerical simulation and control of asymmetric flows around slender bodies

The problems of asymmetric flow around slender bodies and its control are formulated using the unsteady, compressible, thin-layer or full Navier-Stokes equations which are solved using an implicit, flux-difference splitting, finite-volume scheme. The problem is numerically simulated for both locally-conical and three-dimensional flows. The numerical applications include studies of the effects of relative incidence, Mach number and Reynolds number on the flow asymmetry. For the control of flow asymmetry, the numerical simulation cover passive and active control methods. For the passive control, the effectiveness of vertical fins placed in the leeward plane of geometric symmetry and side strakes with different orientations is studied. For the active control, the effectiveness of normal and tangential flow injection and surface heating and a combination of these methods is studied.

Kandil, Osama A.↗

High Order Approximations for Compressible Fluid Dynamics on Unstructured and Cartesian Meshes

The development of high-order accurate numerical discretization techniques for irregular domains and meshes is often cited as one of the remaining challenges facing the field of computational fluid dynamics. In structural mechanics, the advantages of high-order finite element approximation are widely recognized. This is especially true when high-order element approximation is combined with element refinement (h-p refinement). In computational fluid dynamics, high-order discretization methods are infrequently used in the computation of compressible fluid flow. The hyperbolic nature of the governing equations and the presence of solution discontinuities makes high-order accuracy difficult to achieve. Consequently, second-order accurate methods are still predominately used in industrial applications even though evidence suggests that high-order methods may offer a way to significantly improve the resolution and accuracy for these calculations. To address this important topic, a special course was jointly organized by the Applied Vehicle Technology Panel of NATO's Research and Technology Organization (RTO), the von Karman Institute for Fluid Dynamics, and the Numerical Aerospace Simulation Division at the NASA Ames Research Center. The NATO RTO sponsored course entitled "Higher Order Discretization Methods in Computational Fluid Dynamics" was held September 14-18, 1998 at the von Karman Institute for Fluid Dynamics in Belgium and September 21-25, 1998 at the NASA Ames Research Center in the United States. During this special course, lecturers from Europe and the United States gave a series of comprehensive lectures on advanced topics related to the high-order numerical discretization of partial differential equations with primary emphasis given to computational fluid dynamics (CFD). Additional consideration was given to topics in computational physics such as the high-order discretization of the Hamilton-Jacobi, Helmholtz, and elasticity equations. This volume consists of five articles prepared by the special course lecturers. These articles should be of particular relevance to those readers with an interest in numerical discretization techniques which generalize to very high-order accuracy. The articles of Professors Abgrall and Shu consider the mathematical formulation of high-order accurate finite volume schemes utilizing essentially non-oscillatory (ENO) and weighted essentially non-oscillatory (WENO) reconstruction together with upwind flux evaluation. These formulations are particularly effective in computing numerical solutions of conservation laws containing solution discontinuities. Careful attention is given by the authors to implementational issues and techniques for improving the overall efficiency of these methods. The article of Professor Cockburn discusses the discontinuous Galerkin finite element method. This method naturally extends to high-order accuracy and has an interpretation as a finite volume method. Cockburn addresses two important issues associated with the discontinuous Galerkin method: controlling spurious extrema near solution discontinuities via "limiting" and the extension to second order advective-diffusive equations (joint work with Shu). The articles of Dr. Henderson and Professor Schwab consider the mathematical formulation and implementation of the h-p finite element methods using hierarchical basis functions and adaptive mesh refinement. These methods are particularly useful in computing high-order accurate solutions containing perturbative layers and corner singularities. Additional flexibility is obtained using a mortar FEM technique whereby nonconforming elements are interfaced together. Numerous examples are given by Henderson applying the h-p FEM method to the simulation of turbulence and turbulence transition.

Barth, Timothy↗

Relative motion characteristics of 2 near-Earth Satellites

The stability of the nonlinear dynamical system of two GRAVSAT - type satellites was investigated by performing several numerical experiments which provide the simulations of the relative motion characteristics between the two satellites for various specified time intervals. The simulations included the relative range, range-rate, and relative acceleration magnitude. These simulations were generated with respect to appropriate initial orbital elements which were obtained such that the instantaneous separation distance between the two satellites has small fluctuations from a specified constant separation distance. The simulation results indicate that the behavior of the relative motions is very sensitive to the initial orbital elements of the satellites and that for a specified time interval of interest. A stable behavior is possible only with the use of an appropriate set of initial orbital elements compatible with the gravity field used to derive them.

Schutz, B. E.↗

Radiation Heat Transfer Procedures for Space-Related Applications

Over the last contract year, a numerical procedure for combined conduction-radiation heat transfer using unstructured grids has been developed. As a result of this research, one paper has been published in the Numerical Heat Transfer Journal. One paper has been accepted for presentation at the International Center for Heat and Mass Transfer's International Symposium on Computational Heat Transfer to be held in Australia next year. A journal paper is under review by my NASA's contact. A conference paper for the ASME National Heat Transfer conference is under preparation. In summary, a total of four (4) papers (two journal and two conference) have been published, accepted or are under preparation. There are two (2) to three (3) more papers to be written for the project. In addition to the above publications, one book chapter, one journal paper and six conference papers have been published as a result of this project. Over the last contract year, the research project resulted in one Ph.D. thesis and partially supported another Ph.D. student. My NASA contact and myself have formulated radiation heat transfer procedures for materials with different indices of refraction and for combined conduction-radiation heat transfer. We are trying to find other applications for the procedures developed under this grant.

Chai, John C.↗

Relating a Jet-Surface Interaction Experiment to a Commercial Supersonic Transport Aircraft Using Numerical Simulations

NASA and industry partners desire to reintroduce commercial supersonic airliners to the air transportation system. There are a number of technical challenges that must be overcome by future commercial supersonic airliners to make them viable solutions in society. NASA is specifically concerned with the challenges of reducing boom during supersonic cruise, maximizing range, and reducing airport community noise to acceptable levels. Concepts for commercial supersonic transports, such as the concept aircraft by Lockheed Martin pictured in Figure 1, place the engine nozzles in close proximity to wing and tail surfaces. However, the effects of noise shielding and noise radiation are not fully understood for installed propulsion systems. A series of acoustic tests were conducted on the NASA Glenn Research Centers Nozzle Acoustic Test Rig (NATR) to address the challenge of reducing airport community noise, which is often dominated by jet noise. To best represent the conceptual aircraft in the acoustic tests, noise measurements were taken of the jet in close proximity of simulated aerodynamic surfaces, not simply of an isolated jet.

nozzle flow↗

Volcanism in the Noachis-Hellas region of Mars, 2

Mariner 9 and Viking Orbiter imagery has been examined to investigate volcanism in the Noachis-Hellas region of Mars. A wide variety of probable and possible volcanic features, including cinder cones, flood lavas, wrinkel ridges, cratered domes on the floors of impact craters, shield volcanoes, and possible stratovolcanoes has been observed. The flood lavas of ridged plains and broad, low shield volcanoes are indicative of very fluid lavas of ultrabasic or very basic basaltic composition. Evidence for relatively recent pyroclastic volcanism is present in numerous probable cinder cones and some isolated cratered mountains tentatively identified as stratovolcanoes. The volcanic complex and ridged plains at 68 deg S, 323 deg W appear to be relatively young. The occurrence of the five major volcanic centers of the Noachis-Hellas region at the intersections of ancient Hellas impact basin rings with a postulated northeast-trending fracture system suggests that large-scale volcanic activity may have begun soon after formation of the basin.

Peterson, J. E.↗

External relations

Mechanisms for a variety of external relations were developed by NASA for programs that clearly would involve more people outside NASA than in, and in which numerous other agencies would have strong interests. The relation between NASA and the Department of Defense, the Academy of Science, and the President's Science Advisory Committe are summarized. The publication of results and technology transfer between NASA and universities and other scientific societies are discussed.

Source record↗

Water-Leaving Contribution to Polarized Radiation Field Over Ocean

The top-of-atmosphere (TOA) radiation field from a coupled atmosphere-ocean system (CAOS) includes contributions from the atmosphere, surface, and water body. Atmo-spheric correction of ocean color imagery is to retrieve water-leaving radiance from the TOA measurement, from which ocean bio-optical properties can be obtained. Knowledge of the ab-solute and relative magnitudes of water-leaving signal in the TOA radiation field is important for designing new atmospheric correction algorithms and developing retrieval algorithms for new ocean biogeochemical parameters. In this paper we present a systematic sensitivity study of water-leaving contribution to the TOA radiation field, from 340 nm to 865 nm, with polarization included. Ocean water inherent optical properties are derived from bio-optical models for two kinds of waters, one dominated by phytoplankton (PDW) and the other by non-algae particles (NDW). In addition to elastic scattering, Raman scattering and fluorescence from dissolved organic matter in ocean waters are included. Our sensitivity study shows that the polarized reflectance is minimized for both CAOS and ocean signals in the backscattering half plane, which leads to numerical instability when calculating water leaving relative contribution, the ratio between polarized water leaving and CAOS signals. If the backscattering plane is excluded, the water-leaving polarized signal contributes less than 9% to the TOA polarized reflectance for PDW in the whole spectra. For NDW, the polarized water leaving contribution can be as much as 20% in the wavelength range from 470 to 670 nm. For wavelengths shorter than 452 nm or longer than 865 nm, the water leaving contribution to the TOA polarized reflectance is in general smaller than 5% for NDW. For the TOA total reflectance, the water-leaving contribution has maximum values ranging from 7% to 16% at variable wavelengths from 400 nm to 550 nm from PDW. The water leaving contribution to the TOA total reflectance can be as large as 35%for NDW, which is in general peaked at 550 nm. Both the total and polarized reflectances from water-leaving contributions approach zero in the ultraviolet and near infrared bands. These facts can be used as constraints or guidelines when estimating the water leaving contribution to the TOA reflectance for new atmospheric correction algorithms for ocean color imagery.

inelastic scatterin↗

Theoretical overview and modeling of the sodium and potassium atmospheres of mercury

A general theoretical overview for the sources, sinks, gas-surface interactions, and transport dynamics of sodium and potassium in the exospheric atmsophere of Mercury is given. Information for these four factors, which control the spatial distribution of these two alkali-group gases about the planet, is incorporated in numerical models. The spatial nature and relative importance of the initial source atom atmosphere and the ambient (ballistic hopping) atom atmosphere are then examined and are shown to be controlled and coupled to a great extent by the extremely large and variable solar radiation acceleration experienced by sodium and potassium as they resonantly scatter solar photons. The lateral (antisunward) transport rate of thermally accommodated sodium and potassium ambient atoms is shown to be driven by the solar radiation acceleration and, over a significant portion of Mercury's orbit about the Sun, is sufficiently rapid to be competitive with the short photoionization lifetimes for these atoms when they are located on the summit surface near or within about 30 deg of the terminator. The lateral transport rate is characterized by a migration time determined by model calculations for an ensemble of atoms initially starting at a point source on the surface (i.e., a numerical spacetime dependent Green's function). Four animations for the spacetime evolution of the sodium (or potassium) atmosphere produced by a point source on the surface are presented on a videotape format. For extended surface sources for sodium and potassium, the local column density is determined by competition between the photoionization lifetimes and the lateral transport times of atoms originating from different surface source locations. Sodium surface source fluxes (referenced to Mercury at perihelion) that are required on the sunlit hemisphere to reproduce the typically observed several megarayleighs of D2 emission-line brightness and the inferred column densities of 1-2 x 10(exp 11) atoms per sq cm range from approximately 2-5 x 10(exp 7) atoms/sq cm/sec. The sodium model is applied to study observational data that document an anticorrelation in the average sodium column density and solar radiation acceleration. Lateral transport driven by the solar radiation acceleration is shown to produce this behavior for combinations of different sources and surface accomodation coefficients. The best fit model fits to the observational data require a significant degree of thermal accommodation of the ambient sodium atoms to the surface and a source rate that decreases as an inverse power of 1.5 to 2 in heliocentric distance.

Smyth, William H.↗

Ice-Over-Water Cloud Properties in an Artificial Neural Network Approach

Clouds are a crucial component of the atmospheric energy system, particularly the radiative balance within, above, and below the troposphere. The vertical distribution of cloud mass and phase determines layer heating rates, the loss of radiation to space, and the amount of radiative heating at the surface. Thus, it is important to know how clouds are distributed both vertically and horizontally at all times of day. Satellite remote sensing is the only approach available to monitor clouds day and night around the globe. In this paper several artificial neural network (ANN) algorithms, employing several Aqua MODIS infrared channels, profiles of relative humidity and temperature from GMAO numerical weather analyses, and the retrieved total cloud visible optical depth, are trained to detect multilayer ice-over-water cloud systems and to retrieve some of their properties as identified by a year of 2008 Aqua MODIS data matched with CloudSat and CALIPSO (CC) cloud profiles. The CC lidar and radar profiles provide the vertical structure that serves as output truth for the multilayer algorithm. The neural networks were trained using one year (2008) of cloud top height data from the CC dataset, with correlation around 0.94 (0.95) and MAE as low as 0.82 (0.81) km for nonpolar regions during the day (night). Applying the trained ANN to independent year 2009 MODIS data resulted in a combined ML and single layer hit rate of 86.4% (85.1%) for nonpolar regions during the day (night). Since the ANN is trained using near-nadir MODIS pixels, infrared radiance corrections were developed as a function of view zenith angle from MODIS and applied to off-nadir pixels when processing MODIS swath data. The multilayer amount derived with the ANN is relatively invariant with increasing view zenith angle compared to the multilayer amount without the corrections.

MODIS↗

A Numerical Climate Observing Network Design Study

This project was concerned with three related questions of an optimal design of a climate observing system: 1. The spatial sampling characteristics required from an ARGO system. 2. The degree to which surface observations from ARGO can be used to calibrate and test satellite remote sensing observations of sea surface salinity (SSS) as it is anticipated now. 3. The more general design of an climate observing system as it is required in the near future for CLIVAR in the Atlantic. An important question in implementing an observing system is that of the sampling density required to observe climate-related variations in the ocean. For that purpose this project was concerned with the sampling requirements for the ARGO float system, but investigated also other elements of a climate observing system. As part of this project we studied the horizontal and vertical sampling characteristics of a global ARGO system which is required to make it fully complementary to altimeter data with the goal to capture climate related variations on large spatial scales (less thanAttachment: 1000 km). We addressed this question in the framework of a numerical model study in the North Atlantic with an 1/6 horizontal resolution. The advantage of a numerical design study is the knowledge of the full model state. Sampled by a synthetic float array, model results will therefore allow to test and improve existing deployment strategies with the goal to make the system as optimal and cost-efficient as possible. Attachment: "Optimal observations for variational data assimilation".

Stammer, Detlef↗

Stability analysis for laminar flow control, part 2

Topics covered include: (1) optimization of the numerics of the SALLY stability analysis code; (2) relation between temporal and spatial stability theory; (3) compressible flow stability calculations; (4) spectral methods for the boundary layer equations; and (5) numerical study of nonlinear, nonparallel stability of incompressible flows.

Orszag, S. A.↗

Structural tailoring/analysis for hypersonic components - Executive system development

No direct analytical or integrated numerical tool exists today for the optimal design of a generic class of built-up actively cooled composite structure for applications in hypersonic propulsion ducts. The need exists for a numerical tool to perform the comprehensive design/analysis of a panel on the inlet wall under hypersonic flight conditions. Such a tool requires relatively complex multi-disciplinary analysis. One such numerical tool controlled by an executive system has been developed and is named as STAHYC (Structural Tailoring/Analysis for HYpersonic Components). A detailed account of the executive system development of STAHYC along with the results of one example inlet panel design problem is given in this paper.

Narayanan, G. V.↗

Multi-dimensional high order essentially non-oscillatory finite difference methods in generalized coordinates

The nonlinear stability of compact schemes for shock calculations is investigated. In recent years compact schemes were used in various numerical simulations including direct numerical simulation of turbulence. However to apply them to problems containing shocks, one has to resolve the problem of spurious numerical oscillation and nonlinear instability. A framework to apply nonlinear limiting to a local mean is introduced. The resulting scheme can be proven total variation (1D) or maximum norm (multi D) stable and produces nice numerical results in the test cases. The result is summarized in the preprint entitled 'Nonlinearly Stable Compact Schemes for Shock Calculations', which was submitted to SIAM Journal on Numerical Analysis. Research was continued on issues related to two and three dimensional essentially non-oscillatory (ENO) schemes. The main research topics include: parallel implementation of ENO schemes on Connection Machines; boundary conditions; shock interaction with hydrogen bubbles, a preparation for the full combustion simulation; and direct numerical simulation of compressible sheared turbulence.

Shu, Chi-Wang↗

An analysis of the flow field near the fuel injection location in a gas core reactor.

An analytical study is presented which shows the effects of large energy release and the concurrent high acceleration of inner stream fluid on the coaxial flow field in a gas core reactor. The governing equations include the assumptions of only radial radiative transport of energy represented as an energy diffusion term in the Euler equations. The method of integral relations is used to obtain the numerical solution. Results show that the rapidly accelerating, heat generating inner stream actually shrinks in radius as it expands axially.

Weinstein, H.↗

Outline of a new approach to the analysis of complex systems and decision processes.

Development of a conceptual framework for dealing with systems which are too complex or too ill-defined to admit of precise quantitative analysis. The approach outlined is based on the premise that the key elements in human thinking are not numbers, but labels of fuzzy sets - i.e., classes of objects in which the transition from membership to nonmembership is gradual rather than abrupt. The approach in question has three main distinguishing features - namely, the use of so-called 'linguistic' variables in place of or in addition to numerical variables, the characterization of simple relations between variables by conditional fuzzy statements, and the characterization of complex relations by fuzzy algorithms.

Zadeh, L. A.↗

Evaluation of analytical procedures for prediction of turbulent boundary layers on a porous wall

An analytical study has been made to determine how well current boundary layer prediction techniques work when there is mass transfer normal to the wall. The data that were considered in this investigation were for two-dimensional, incompressible, turbulent boundary layers with suction and blowing. Some of the bleed data were taken in an adverse pressure gradient. An integral prediction method was used three different porous wall skin friction relations, in addition to a solid-surface relation for the suction cases. A numerical prediction method was also used. Comparisons were made between theoretical and experimental skin friction coefficients, displacement and momentum thicknesses, and velocity profiles. The integral method with one of the porous wall skin friction laws gave very good agreement with data for most of the cases considered. The use of the solid-surface skin friction law caused the integral to overpredict the effectiveness of the bleed. The numerical techniques also worked well for most of the cases.

Towne, C. E.↗