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 937 records · Page 52

Integrated identification and robust control tuning for large space structures

System identification is studied for the explicit purpose of supporting modern H-infinity robust control design objectives. In the analysis, the true plant is not assumed to be in the identification model set. An integrated identification/robust control problem is posed in which the optimal solution guarantees the best robust performance relative to the system information contained in a given experimental data set. A numerical example demonstrating an approximate solution to the problem indicates the usefulness of the approach.

Yam, Y.↗

Multiblock grid generation for jet engine configurations

The goal was to create methods for generating grids with minimal human intervention that are applicable to a wide range of problems and compatible with existing numerical methods and with existing and proposed computers. The following topics that are related to multiblock grid generation are briefly covered in viewgraph form: finding a domain decomposition, dimensioning grids, grid smoothing, manipulating grids and decompositions, and some specializations for jet engine configurations.

Stewart, Mark E. M.↗

Flexure-torsion behavior of prismatic beams. I - Section properties via power series

The behavior of a tip-loaded cantilever beam with an arbitrary cross section is studied using Saint-Venant's semi-inverse method along with a power series solution for the out-of-plane flexure and torsion warping functions. The power series coefficients are determined by solving a set of variationally derived linear algebraic equations. For complex cross sections, the calculated coefficients represented a 'best-fit approximation' to the exact warping function. The resulting warping functions are used to determine the cross-sectional properties (torsion constant, shear correction factors, shear deformation coefficients, and shear center location). A new linear relation is developed for locating the shear center, where the twist rate is zero about the line of shear centers. Moreover, the kinematic relations for a new fully compatible one-dimensional beam theory are developed. Numerical results are presented first to verify the approach and second to provide section data on NACA four-series airfoils not currently found in the literature.

Kosmatka, J. B.↗

A global analytical representation of the magnetic field produced by the region 2 Birkeland currents and the partial ring current

A quantitative model is developed of the magnetic field produced by the electric current system of region 2 Birkeland currents, closed via the partial ring current. The distribution of j-perpendicular is computed from a given axially asymmetric spatial distribution of hot isotropic magnetospheric plasma over an infinitely thin L shell in an axisymmetric purely dipolar geomagnetic field, while the field-aligned current density is found from the continuity of the net electric current. The magnetic field distribution is derived by a Biot-Savart integral over the electric current system. An assumed cosine dependence of the plasma pressure on local time makes it possible to reduce the problem of analytical representation of the B field to two dimensions. The obtained numerical fits for the partial ring current/region 2 Birkeland current magnetic field are relatively simple, continuous, and valid throughout the whole extraterrestrial space from ionospheric heights up to tens of earth radii.

Tsyganenko, N. A.↗

A fast Poisson solver for unsteady incompressible Navier-Stokes equations on the half-staggered grid

In this paper, a fast Poisson solver for unsteady, incompressible Navier-Stokes equations with finite difference methods on the non-uniform, half-staggered grid is presented. To achieve this, new algorithms for diagonalizing a semi-definite pair are developed. Our fast solver can also be extended to the three dimensional case. The motivation and related issues in using this second kind of staggered grid are also discussed. Numerical testing has indicated the effectiveness of this algorithm.

Golub, G. H.↗

Conflict Probability Estimation for Free Flight

The safety and efficiency of free flight will benefit from automated conflict prediction and resolution advisories. Conflict prediction is based on trajectory prediction and is less certain the farther in advance the prediction, however. An estimate is therefore needed of the probability that a conflict will occur, given a pair of predicted trajectories and their levels of uncertainty. A method is developed in this paper to estimate that conflict probability. The trajectory prediction errors are modeled as normally distributed, and the two error covariances for an aircraft pair are combined into a single equivalent covariance of the relative position. A coordinate transformation is then used to derive an analytical solution. Numerical examples and Monte Carlo validation are presented.

Paielli, Russell A.↗

The Complex Stratigraphy of the Highland Crust in the Serenitatis Region of the Moon Inferred from Mineral Fragment Chemistry

Large impact basins are natural drill holes into the Moon, and their ejecta carries unique information about the rock types and stratigraphy of the lunar crust. We have conducted an electron microprobe study of mineral fragments in the poikilitic melt breccias collected from the Taurus Mountains at the Apollo 17 landing site. These breccias are virtually unanimously agreed to be impact melt produced in the Serenitatis impact event. They contain lithic fragments and much more abundant mineral fragments of crustal origin. We have made precise microprobe analyses of minor element abundances in fragments of olivine, pyroxene, and plagioclase to provide new information on the possible source rocks and the crustal stratigraphy in the Serenitatis region. These data were also intended to elucidate the nature of the cryptic geochemical component in breccias such as these with low-K Fra Mauro basalt compositions. We chose the finest-grained (i.e., most rapidly quenched) breccias for study, to avoid reacted and partly assimilated fragments as much as possible. Most of the mineral fragments appear to have been derived from rocks that would fall into the pristine igneous Mg-suite as represented by lithic fragments in the Apollo collection, or reasonable extensions of it. Gabbroic rocks were more abundant in the target stratigraphy than is apparent from the Apollo sample collection. Some pyroxene and plagiociase, but probably not much olivine, could be derived from feldspathic granulites, which are metamorphosed polymict breccias. Some mineral fragments are from previously unknown rocks. These include highly magnesian olivines (up to Fo(sub 94)), possibly volcanic in origin, that exacerbate the difficulty in explaining highly magnesian rocks in the lunar crust. It appears that some part of the lunar interior has an mg*(= 100 x Mg/(Mg/Fe) atomic) greater than the conventional bulk Moon value of 80-84. Other volcanic rocks, including mare basalts, and rapidly- cooled impact melt rocks do not contribute significantly to the fragment population. Nor do ferroan anorthosites contribute more than a tiny part of even the plagiociase fragment population. A few mineral fragments that are consistent with the cryptic low-K Fra Mauro chemical component were found, and these appear to be from gabbroic sources. The mineral fragment populations cannot be mixed in their observed proportions to produce the whole rock composition, because the fragments are more refractory and deficient in Ti, P, and alkalis. A preferential contribution to the melt from a rock similar to sodic ferrogabbro can partly resolve the discrepancy. The population of mineral fragments requires a very diverse population of igenous rocks that are not all related to each other, demonstrating the existence of a complex crust built of numerous separate igneous plutons. Many of these plutons may have crystallized at shallow depths. The chemical composition of the melt breccias, in combination with the mineral fragment data and an understanding of the cratering process, suggests that the deepest crust sampled by the Serenitatis impace (not necessarily the deepest crust) was basaltic in composition, including KREEP and gabbroic rocks like sodic ferrogabbro, and lacking abundant olivine-rich material. These were overlain by Mg-suite rocks of varied types, including norites and troctolites that supplied most of the olivine mineral fragments. Granulities, which are metamorphosed and more feldspathic breccias, were abundant near the surface. Remote sensing indicates that the entire Serenitatis region lacks ferroan anorthosite, consistent with the results of our study.

Ryder, Graham↗

Parametric Study of Three-Stringer-Panel Compression-After-Impact Strength

Damage tolerance requirements for integrally-stiffened composite wing skins are typically met using design allowables generated by testing impact-damaged subcomponents, such as three-stringer stiffened panels. To improve these structures, it is necessary to evaluate the critical design parameters associated with three-stringer stiffened-panel compressive behavior. During recent research and development programs, four structural parameters were identified as sources for strength variation: (a) material system, (b) stringer configuration, (c) skin layup, and (d) form of axial reinforcement (tape versus pultruded carbon rods). Relative effects of these parameters on damage resistance and damage tolerance were evaluated numerically and experimentally. Material system and geometric configuration had the largest influence on damage resistance; location and extent of the damage zone influenced the sublaminate buckling behavior, failure initiation site, and compressive ultimate strength. A practical global-local modeling technique captured observed experimental behavior and has the potential to identify critical damage sites and estimate failure loads prior to testing. More careful consideration should be given to accurate simulation of boundary conditions in numerical and experimental studies.

Rousseau, Carl Q.↗

A Real-Time Method for Estimating Viscous Forebody Drag Coefficients

This paper develops a real-time method based on the law of the wake for estimating forebody skin-friction coefficients. The incompressible law-of-the-wake equations are numerically integrated across the boundary layer depth to develop an engineering model that relates longitudinally averaged skin-friction coefficients to local boundary layer thickness. Solutions applicable to smooth surfaces with pressure gradients and rough surfaces with negligible pressure gradients are presented. Model accuracy is evaluated by comparing model predictions with previously measured flight data. This integral law procedure is beneficial in that skin-friction coefficients can be indirectly evaluated in real-time using a single boundary layer height measurement. In this concept a reference pitot probe is inserted into the flow, well above the anticipated maximum thickness of the local boundary layer. Another probe is servomechanism-driven and floats within the boundary layer. A controller regulates the position of the floating probe. The measured servomechanism position of this second probe provides an indirect measurement of both local and longitudinally averaged skin friction. Simulation results showing the performance of the control law for a noisy boundary layer are then presented.

Whitmore, Stephen A.↗

A Smoothed Boundary Condition for Reducing Nonphysical Field Effects

In this paper, we examine the problem associated with abruptly mixing boundary conditions in the context of a two-dimensional semiconductor device simulator. Explicitly, this paper addresses the transition between an ohmic-type Dirichlet condition and a passivated Neumann boundary. In the traditional setting, the details or the transition between the two boundary types are not addressed and an abrupt transition is assumed. Subsequently, the calculated observables (most notably the potential) exhibit discontinuous derivatives near the surface at the point where the boundary type switches. This paper proposes an alternative condition which models the progression between the two boundary types through the use of a finite length, smoothed boundary whereby the numerical discontinuities are eliminated. The physical and mathematical basis for this smoothed boundary condition is discussed and examples of the technique's implementation given. It is found that the proposed boundary condition is numerically efficient and can be implemented in pre-existing device simulators with relative ease.

Smith, Arlynn W.↗

Software for Refining or Coarsening Computational Grids

A computer program performs calculations for refinement or coarsening of computational grids of the type called 'structured' (signifying that they are geometrically regular and/or are specified by relatively simple algebraic expressions). This program is designed to facilitate analysis of the numerical effects of changing structured grids utilized in computational fluid dynamics (CFD) software. Unlike prior grid-refinement and -coarsening programs, this program is not limited to doubling or halving: the user can specify any refinement or coarsening ratio, which can have a noninteger value. In addition to this ratio, the program accepts, as input, a grid file and the associated restart file, which is basically a file containing the most recent iteration of flow-field variables computed on the grid. The program then refines or coarsens the grid as specified, while maintaining the geometry and the stretching characteristics of the original grid. The program can interpolate from the input restart file to create a restart file for the refined or coarsened grid. The program provides a graphical user interface that facilitates the entry of input data for the grid-generation and restart-interpolation routines.

Daines, Russell↗

A Cartesian Adaptive Level Set Method for Two-Phase Flows

In the present contribution we develop a level set method based on local anisotropic Cartesian adaptation as described in Ham et al. (2002). Such an approach should allow for the smallest possible Cartesian grid capable of resolving a given flow. The remainder of the paper is organized as follows. In section 2 the level set formulation for free surface calculations is presented and its strengths and weaknesses relative to the other free surface methods reviewed. In section 3 the collocated numerical method is described. In section 4 the method is validated by solving the 2D and 3D drop oscilation problem. In section 5 we present some results from more complex cases including the 3D drop breakup in an impulsively accelerated free stream, and the 3D immiscible Rayleigh-Taylor instability. Conclusions are given in section 6.

Ham, F.↗

World Ocean Circulation Experiment (WOCE) Young Investigator Workshops

The World Ocean Circulation Experiment (WOCE) Young Investigator Workshops goals and objectives are: a) to familiarize Young Investigators with WOCE models, datasets and estimation procedures; b) to offer intensive hands-on exposure to these models ard methods; c) to build collaborations among junior scientists and more senior WOCE investigators; and finally, d) to generate ideas and projects leading to fundable WOCE synthesis projects. To achieve these goals and objectives, the Workshop will offer a mixture of tutorial lectures on numerical models and estimation procedures, advanced seminars on current WOCE synthesis activities and related projects, and the opportunity to conduct small projects which put into practice the techniques advanced in the lectures.

Austin, Meg↗

Software for Refining or Coarsening Computational Grids

A computer program performs calculations for refinement or coarsening of computational grids of the type called "structured" (signifying that they are geometrically regular and/or are specified by relatively simple algebraic expressions). This program is designed to facilitate analysis of the numerical effects of changing structured grids utilized in computational fluid dynamics (CFD) software. Unlike prior grid-refinement and -coarsening programs, this program is not limited to doubling or halving: the user can specify any refinement or coarsening ratio, which can have a noninteger value. In addition to this ratio, the program accepts, as input, a grid file and the associated restart file, which is basically a file containing the most recent iteration of flow-field variables computed on the grid. The program then refines or coarsens the grid as specified, while maintaining the geometry and the stretching characteristics of the original grid. The program can interpolate from the input restart file to create a restart file for the refined or coarsened grid. The program provides a graphical user interface that facilitates the entry of input data for the grid-generation and restart-interpolation routines.

Daines, Russell↗

Hubble Space Telescope Imaging of the Circumnuclear Environments of the CfA Seyfert Galaxies: Nuclear Spirals and Fueling

We present archival Hubble Space Telescope (HST) images of the nuclear regions of 43 of the 46 Seyfert galaxies found in the volume limited,spectroscopically complete CfA Redshift Survey sample. Using an improved method of image contrast enhancement, we created detailed high-quality " structure maps " that allow us to study the distributions of dust, star clusters, and emission-line gas in the circumnuclear regions (100-1000 pc scales) and in the associated host galaxy. Essentially all of these Seyfert galaxies have circumnuclear dust structures with morphologies ranging from grand-design two-armed spirals to chaotic dusty disks. In most Seyfert galaxies there is a clear physical connection between the nuclear dust spirals on hundreds of parsec scales and large-scale bars and spiral arms in the host galaxies proper. These connections are particularly striking in the interacting and barred galaxies. Such structures are predicted by numerical simulations of gas flows in barred and interacting galaxies and may be related to the fueling of active galactic nuclei by matter inflow from the host galaxy disks. We see no significant differences in the circumnuclear dust morphologies of Seyfert 1s and 2s, and very few Seyfert 2 nuclei are obscured by large-scale dust structures in the host galaxies. If Sevfert 2s are obscured Sevfert Is, then the obscuration must occur on smaller scales than those probed by HST.

Pogge, Richard W.↗

Percentage Contributions from Atmospheric and Surface Features to Computed Brightness Temperatures

Over the past few years, there has become an increasing interest in the use of millimeter-wave (mm-wave) and sub-millimeter-wave (submm-wave) radiometer observations to investigate the properties of ice particles in clouds. Passive radiometric channels respond to both the integrated particle mass throughout the volume and field of view, and to the amount, location, and size distribution of the frozen (and liquid) particles with the sensitivity varying for different frequencies and hydrometeor types. One methodology used since the 1960's to discern the relationship between the physical state observed and the brightness temperature (TB) is through the temperature weighting function profile. In this research, the temperature weighting function concept is exploited to analyze the sensitivity of various characteristics of the cloud profile, such as relative humidity, ice water path, liquid water path, and surface emissivity. In our numerical analysis, we compute the contribution (in Kelvin) from each of these cloud and surface characteristics, so that the sum of these various parts equals the computed TB. Furthermore, the percentage contribution from each of these characteristics is assessed. There is some intermingling/contamination of the contributions from various components due to the integrated nature of passive observations and the absorption and scattering between the vertical layers, but all in all the knowledge gained is useful. This investigation probes the sensitivity over several cloud classifications, such as cirrus, blizzards, light snow, anvil clouds, and heavy rain. The focus is on mm-wave and submm-wave frequencies, however discussions of the effects of cloud variations to frequencies as low as 10 GHz and up to 874 GHz will also be presented. The results show that nearly 60% of the TB value at 89 GHz comes from the earth's surface for even the heaviest blizzard snow rates. On the other hand, a significant percentage of the TB value comes from the snow in the cloud for 166, and 183 plus or minus 7 GHz for the heavy and medium snow rates. For submm-wave channels, there is no contribution from the surface because these channels cannot probe through clouds, nor normal water vapor amounts in clear air regions. This work is extremely valuable in physically-based retrieval algorithm development research.

Gail Skofronick-Jackson↗

Thin Shell, Segmented X-Ray Mirrors

Thin foil mirrors were introduced as a means of achieving high throughput in an X-ray astronomical imaging system in applications for which high angular resolution were not necessary. Since their introduction, their high filling factor, modest mass, relative ease of construction, and modest cost have led to their use in numerous X-ray observatories, including the Broad Band X-ray Telescope, ASCA, and Suzaku. The introduction of key innovations, including epoxy replicated surfaces, multilayer coatings, and glass mirror substrates, has led to performance improvements, and in their becoming widely used for X-ray astronomical imaging at energies above 10 keV. The use of glass substrates has also led to substantial improvement in angular resolution, and thus their incorporation into the NASA concept for the International X-ray Observatory with a planned 3 in diameter aperture. This paper traces the development of foil mirrors from their inception in the 1970's through their current and anticipated future applications.

Petre, Robert↗

Software for Refining or Coarsening Computational Grids

A computer program performs calculations for refinement or coarsening of computational grids of the type called structured (signifying that they are geometrically regular and/or are specified by relatively simple algebraic expressions). This program is designed to facilitate analysis of the numerical effects of changing structured grids utilized in computational fluid dynamics (CFD) software. Unlike prior grid-refinement and -coarsening programs, this program is not limited to doubling or halving: the user can specify any refinement or coarsening ratio, which can have a noninteger value. In addition to this ratio, the program accepts, as input, a grid file and the associated restart file, which is basically a file containing the most recent iteration of flow-field variables computed on the grid. The program then refines or coarsens the grid as specified, while maintaining the geometry and the stretching characteristics of the original grid. The program can interpolate from the input restart file to create a restart file for the refined or coarsened grid. The program provides a graphical user interface that facilitates the entry of input data for the grid-generation and restart-interpolation routines.

Daines, Russell↗