Search NASA⌕ Search

SEARCH · Search NASA

Results for “Local coordinate transform”

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

Time simulation of flutter with large stiffness changes

Time simulation of flutter, involving large local structural changes, is formulated with a state-space model that is based on a relatively small number of generalized coordinates. Free-free vibration modes are first calculated for a nominal finite-element model with relatively large fictitious masses located at the area of structural changes. A low-frequency subset of these modes is then transformed into a set of structural modal coordinates with which the entire simulation is performed. These generalized coordinates and the associated oscillatory aerodynamic force coefficient matrices are used to construct an efficient time-domain, state-space model for basic aeroelastic case. The time simulation can then be performed by simply changing the mass, stiffness and damping coupling terms when structural changes occur. It is shown that the size of the aeroelastic model required for time simulation with large structural changes at a few a priori known locations is similar to that required for direct analysis of a single structural case. The method is applied to the simulation of an aeroelastic wind-tunnel model. The diverging oscillations are followed by the activation of a tip-ballast decoupling mechanism that stabilizes the system but may cause significant transient overshoots.

Karpel, M.↗

Stability and transition on swept wings

This paper describes the extension and application of the Parabolized Stability Equations (PSE) to the stability and transition of the supersonic three-dimensional laminar boundary layer on a swept wing. The problem formulation uses a general coordinate transformation for arbitrary curvilinear body-fitted computational grids. Some testing using these coordinates is briefly described to help validate the software used for the investigation. The disturbance amplitude ratios as a function of chord position for supersonic (Mach 1.5) boundary layers on untapered, untwisted wings of different sweep angles are then presented and compared with those obtained from local parallel analyses.

Stuckert, Greg↗

Three-Dimensional Modeling of Aircraft High-Lift Components with Vehicle Sketch Pad

Vehicle Sketch Pad (OpenVSP) is a parametric geometry modeler that has been used extensively for conceptual design studies of aircraft, including studies using higher-order analysis. OpenVSP can model flap and slat surfaces using simple shearing of the airfoil coordinates, which is an appropriate level of complexity for lower-order aerodynamic analysis methods. For three-dimensional analysis, however, there is not a built-in method for defining the high-lift components in OpenVSP in a realistic manner, or for controlling their complex motions in a parametric manner that is intuitive to the designer. This paper seeks instead to utilize OpenVSP's existing capabilities, and establish a set of best practices for modeling high-lift components at a level of complexity suitable for higher-order analysis methods. Techniques are described for modeling the flap and slat components as separate three-dimensional surfaces, and for controlling their motion using simple parameters defined in the local hinge-axis frame of reference. To demonstrate the methodology, an OpenVSP model for the Energy-Efficient Transport (EET) AR12 wind-tunnel model has been created, taking advantage of OpenVSP's Advanced Parameter Linking capability to translate the motions of the high-lift components from the hinge-axis coordinate system to a set of transformations in OpenVSP's frame of reference.

Olson, Erik D.↗

Refinement of Methods for Evaluation of Near-Hypersingular Integrals in BEM Formulations

In this paper, we present advances in singularity cancellation techniques applied to integrals in BEM formulations that are nearly hypersingular. Significant advances have been made recently in singularity cancellation techniques applied to 1 R type kernels [M. Khayat, D. Wilton, IEEE Trans. Antennas and Prop., 53, pp. 3180-3190, 2005], as well as to the gradients of these kernels [P. Fink, D. Wilton, and M. Khayat, Proc. ICEAA, pp. 861-864, Torino, Italy, 2005] on curved subdomains. In these approaches, the source triangle is divided into three tangent subtriangles with a common vertex at the normal projection of the observation point onto the source element or the extended surface containing it. The geometry of a typical tangent subtriangle and its local rectangular coordinate system with origin at the projected observation point is shown in Fig. 1. Whereas singularity cancellation techniques for 1 R type kernels are now nearing maturity, the efficient handling of near-hypersingular kernels still needs attention. For example, in the gradient reference above, techniques are presented for computing the normal component of the gradient relative to the plane containing the tangent subtriangle. These techniques, summarized in the transformations in Table 1, are applied at the sub-triangle level and correspond particularly to the case in which the normal projection of the observation point lies within the boundary of the source element. They are found to be highly efficient as z approaches zero. Here, we extend the approach to cover two instances not previously addressed. First, we consider the case in which the normal projection of the observation point lies external to the source element. For such cases, we find that simple modifications to the transformations of Table 1 permit significant savings in computational cost. Second, we present techniques that permit accurate computation of the tangential components of the gradient; i.e., tangent to the plane containing the source element.

Fink, Patricia W.↗

Automatic blocking of nested loops

Blocked algorithms have much better properties of data locality and therefore can be much more efficient than ordinary algorithms when a memory hierarchy is involved. On the other hand, they are very difficult to write and to tune for particular machines. The reorganization is considered of nested loops through the use of known program transformations in order to create blocked algorithms automatically. The program transformations used are strip mining, loop interchange, and a variant of loop skewing in which invertible linear transformations (with integer coordinates) of the loop indices are allowed. Some problems are solved concerning the optimal application of these transformations. It is shown, in a very general setting, how to choose a nearly optimal set of transformed indices. It is then shown, in one particular but rather frequently occurring situation, how to choose an optimal set of block sizes.

Schreiber, Robert↗

Community Coordinated Modeling Center Support of Science Needs for Integrated Data Environment

Space science models are essential component of integrated data environment. Space science models are indispensable tools to facilitate effective use of wide variety of distributed scientific sources and to place multi-point local measurements into global context. The Community Coordinated Modeling Center (CCMC) hosts a set of state-of-the- art space science models ranging from the solar atmosphere to the Earth's upper atmosphere. The majority of models residing at CCMC are comprehensive computationally intensive physics-based models. To allow the models to be driven by data relevant to particular events, the CCMC developed an online data file generation tool that automatically downloads data from data providers and transforms them to required format. CCMC provides a tailored web-based visualization interface for the model output, as well as the capability to download simulations output in portable standard format with comprehensive metadata and user-friendly model output analysis library of routines that can be called from any C supporting language. CCMC is developing data interpolation tools that enable to present model output in the same format as observations. CCMC invite community comments and suggestions to better address science needs for the integrated data environment.

Kuznetsova, M. M.↗

H2, fixed architecture, control design for large scale systems

The H2, fixed architecture, control problem is a classic linear quadratic Gaussian (LQG) problem whose solution is constrained to be a linear time invariant compensator with a decentralized processing structure. The compensator can be made of p independent subcontrollers, each of which has a fixed order and connects selected sensors to selected actuators. The H2, fixed architecture, control problem allows the design of simplified feedback systems needed to control large scale systems. Its solution becomes more complicated, however, as more constraints are introduced. This work derives the necessary conditions for optimality for the problem and studies their properties. It is found that the filter and control problems couple when the architecture constraints are introduced, and that the different subcontrollers must be coordinated in order to achieve global system performance. The problem requires the simultaneous solution of highly coupled matrix equations. The use of homotopy is investigated as a numerical tool, and its convergence properties studied. It is found that the general constrained problem may have multiple stabilizing solutions, and that these solutions may be local minima or saddle points for the quadratic cost. The nature of the solution is not invariant when the parameters of the system are changed. Bifurcations occur, and a solution may continuously transform into a nonstabilizing compensator. Using a modified homotopy procedure, fixed architecture compensators are derived for models of large flexible structures to help understand the properties of the constrained solutions and compare them to the corresponding unconstrained ones.

Mercadal, Mathieu↗

Orbit determination by solving for gravity parameters with multiple arc data

The orbit of a satellite that repeats in the earth fixed coordinates is determined by combining GPS tracking data from multiple arcs. The satellite dynamics are modeled with the epoch state and a set of parameters, called the bin parameters, that account for the effect of the local gravitational field on the satellite current state. The epoch state is specific to each arc, and the bin parameters are common to all repeat arcs. The estimation algorithm is based on the Square Root Information Filter. It involves partitioning of the measurement matrix and use of the Householder transformation to combine multiple arc data and solve for the epoch states and the bin parameters. The bin parameters can then be converted into the earth's gravitational field with a modest amount of computation.

Wu, Jiun-Tsong↗

Solar Occultation Satellite Data and Derived Meteorological Products: Sampling Issues and Comparisons with Aura MLS

Derived Meteorological Products (DMPs, including potential temperature (theta), potential vorticity, equivalent latitude (EqL), horizontal winds and tropopause locations) have been produced for the locations and times of measurements by several solar occultation (SO) instruments and the Aura Microwave Limb Sounder (MLS). DMPs are calculated from several meteorological analyses for the Atmospheric Chemistry Experiment-Fourier Transform Spectrometer, Stratospheric Aerosol and Gas Experiment II and III, Halogen Occultation Experiment, and Polar Ozone and Aerosol Measurement II and III SO instruments and MLS. Time-series comparisons of MLS version 1.5 and SO data using DMPs show good qualitative agreement in time evolution of O3, N2O, H20, CO, HNO3, HCl and temperature; quantitative agreement is good in most cases. EqL-coordinate comparisons of MLS version 2.2 and SO data show good quantitative agreement throughout the stratosphere for most of these species, with significant biases for a few species in localized regions. Comparisons in EqL coordinates of MLS and SO data, and of SO data with geographically coincident MLS data provide insight into where and how sampling effects are important in interpretation of the sparse SO data, thus assisting in fully utilizing the SO data in scientific studies and comparisons with other sparse datasets. The DMPs are valuable for scientific studies and to facilitate validation of non-coincident measurements.

Manney, Gloria↗

A Potentially Useful for Airborne Separation in 4D-Trajectory ATM Operations

An aircraft equipped with Airborne Separation Assistance System functions and 4- dimensional trajectory management capabilities can have significant, potentially transforming, value to Air Traffic Management at the local and system levels. This paper discusses how certain vital characteristics envisioned in the Next Generation Air Transportation System enable some Air Traffic Management functions to be distributed to properly equipped aircraft, and it defines and illustrates this equipage level in a potential application. The new equipage level, perhaps the most capable of many levels permitted, enables an effective implementation of both near- and long-term 4-dimensional trajectory operations in complex airspace, with the aircraft providing the near-term tactical functions and conforming to the long-term trajectory attributes coordinated with ground-based Traffic Flow Management authorities. NASA s recent research and development of this proposed aircraft equipage for en-route and terminal-arrival operations is summarized. The role the equipage level may play in addressing key implementation challenges of reducing ground infrastructure cost, building in security and safety, and scaling to traffic demand is discussed.

Wing, David J.↗

High-Resolution Mapping of Lunar Crustal Magnetic Fields: Correlations with Albedo Markings of the Reiner Gamma Class

During the last eight months of the Lunar Prospector mission (December 1999-July 1999), the spacecraft was placed in a relatively low-altitude (15-30-km perapsis), near-polar orbit that allowed high-resolution mapping of crustal magnetic fields. We report here initial studies of the correlation of locally strong magnetic anomalies with unusual, swirl-like albedo markings of the Reiner Gamma class. Based on this correlation, which is known from earlier studies of Apollo subsatellite magnetometer data, it has been proposed that the swirls represent regions whose higher albedos have been preserved via deflection of the solar-wind ion bombardment by strong crustal fields. This model in turn depends on the hypothesis that solar-wind implanted H is at least one component of the process that optically matures exposed silicate surfaces in the inner solar system . Specifically, it is hypothesized that implanted H acts as an effective reducing agent to enhance the rate of production of nanophase metallic Fe particles from preexisting silicates during micrometeoroid impacts. According to the model, the curvilinear shapes of these albedo markings are caused, at least in part, by the geometry of ion deflections in a magnetic field. The improved resolution and coverage of the Prospector data allow more detailed mapping of the fields, especially on the lunar farside. This permits a more quantitative test of whether all albedo markings of this class are associated with strong local magnetic fields.Only if the latter condition is met can the solar-wind deflection hypothesis he valid. The basic procedure for mapping crustal magnetic fields using Lunar Prospector magnetometer data follows that developed for analysis of Apollo subsatellite magnetometer data. The specific mapping steps are (1) selection of mission time intervals suitable for mapping crustal fields; these are limited essentially either to times when the Moon is in a lobe of the geomagnetic tail or to times when the Moon is in the solar wind but the spacecraft is in the lunar wake; the data are transformed to a radial, east, and north coordinate system with measurements given as a function of spacecraft latitude, longitude, and altitude; (2) visual editing of individual orbit segments selected for minimal external field disturbances; (3) minimization of remaining low-frequency external fields for individual orbit data segments by quadratic detrending; and (4) two-dimensional filtering of individual orbit segments to produce a vector field map along the slightly curved surface defined by the spacecraft altitude; maps of the three field components (radial, east, and north), the field magnitude, and the spacecraft altitude are constructed. For data obtained at low to middle latitudes, the horizontal resolution of the field maps is limited by the orbit-track separation (about 30 km at the equator). Maps of the field magnitude have been constructed within limited selenographic regions based mainly on data acquired in March, April, and May of 1999. This was a time period when the orbit plane was nearly aligned with the Sun-Moon line so that field mapping was possible at times when the Moon was in the solar wind as well as when the Moon was in the geomagnetic tail. Most of the coverage is across the lunar farside. However, a shows an example of a field map produced from solar-wind wake data for a region including Reiner Gamm on western Oceanus Procellarum (location: 58.5W, 7.5N). The contour interval is 3 nT and the mean spacecraft altitude is 18 km to within the accuracy allowed by the resolution of the map (30 km or about 1 deg.); strong magnetic anomalies correlate closely with swirl locations. Individual orbit profiles (whose resolution along the orbit track is comparable to the spacecraft altitude of 18 km) also demonstrate a good correlation of field magnitude with surface albedo. In order to investigate the correlation of magnetic fields with the location of swirl features, we have reexamined available lunar imagery (Lunar Orbiter, Apollo, and Clementine) to identify and map swirl locations within regions where swirls have previously been mapped. In these images, swirls were distinguished from other high-albedo features such as crater rays by their curvilinear shapes and increased visibility in forward-scattered light. Digital maps of swirls identified by all available imagery were then superposed on maps of the field magnitude at the spacecraft altitude. Based upon analysis of these composite magnetic/geologic maps, we draw the preliminary conclusion that swirl features are associated with magnetic anomalies revealed by Lunar Prospector. Detailed maps of these swirl features are currently being constructed for the magnetically strong regions antipodal to the Imbrium, Serenitatis, and Crisium Basins. Additional information contained in the original,

Hood, L. L.↗

The Zircon-Reidite Shock Transformation: A Tem Study of Reidite Produced at A Newly Established 21.2 GPa Shock Pressure Threshold

The high-pressure polymorph of zircon (Zr2SiO4) is reidite, a phase with the scheelite structure that was first discovered to form at ca. 12 GPa in static high-pressure experiments [1]. Reidite was subsequently found to occur in naturally shocked zircons from several terrestrial impact localities [2,3], as well as in shock recovery-reverberation experiments performed above 20 GPa [4]. The nature and distribution of reidite occurrences in terrestrial impact localities has pointed to its potential use as a shock pressure indicator, as well as to the importance of understanding its age-reset effects in zircon geochronology [3]. In this study we have used coordinated SEM/Electron Backscatter Diffraction (EBSD) and Field-Emission Scanning Transmission Electron Microscopy (FE-STEM) to study the micro- and nanostructure of reidite and its shocked zircon host produced at 21.2 GPa under experimental shock recovery/reverberation conditions [5]. Our results document the TEM microstructure of reidite formed at a 21.2 GPa experimental shock pressure that is significantly lower than that observed in previous studies [4]. At these conditions the sample developed a diverse set of microstructures that have not been previously documented on the TEM scale.

shock metamorphism↗

Coordination in Large Collectives

Finding the subset of a set of imperfect devices (e.g., nano or micro devices) that results in the best aggregate device is a challenging problem. It is an abstraction of what will likely be a major difficulty in designing and controlling systems of nano or micro-scale components, particularly when a large fraction of those components may be unreliable. Rather than approaching this as a passive search problem, we transform the problem into one of coordination in a complex system by imbuing each device with simple decision making ability. In doing so we face the challenge of determining what each component should attempt to do so that the collective behavior solves the overall problem. Furthermore, in this padicular instance, we face problems of scaling (number of components in the thousands to tens of thousands), observability (components have limited sensing capabilities), and reliability (the components are faulty). We present an approach based on deriving component goals that are aligned with the overall system goal (e.g., forming best aggregate device), and can be computed using information readily (e.g., locally) available to the components. Then, each component in such a collective uses a simple reinforcement learning algorithm to selfishly pursue its own goals. Because those goals are derived in a principled manner, there is no need to use external mechanisms to force collaboration or coordination among the components to ensure that the system reaches a globally desirable solution. The results show that not only this approach provides improvements of over an order of magnitude over both traditional search methods and traditional multi-agent methods, but that the gains increase with the size of the system. This latter result makes this method ideal for domains where the number of components is currently in the thousands and will reach millions in the near future.

Tumer, Kagan↗

Aeronautical satellite antenna steering using magnetic field sensors

Designers of aeronautical satellite terminals are often faced with the problem of steering a directive antenna from an airplane or helicopter. This problem is usually solved by using aircraft orientation information derived from inertial sensors on-board the aircraft in combination with satellite ephemeris information calculated from geographic coordinates. This procedure works well but relies heavily on avionics that are external to the terminal. For the majority of small aircraft and helicopters which will form the bulk of future aeronautical satcom users, such avionics either do not exist or are difficult for the satellite terminal to interface with. At the Communications Research Center (CRC), work has been undertaken to develop techniques that use the geomagnetic field and satellite antenna pointing vectors (both of which are stationary in a local geographical area) to track the position of a satellite relative to a moving platform such as an aircraft. The performance of this technique is examined and a mathematical steering transformation is developed within this paper. Details are given regarding the experimental program that will be undertaken to test the concepts proposed herein.

Sydor, John↗

Statistical study of ULF wave occurrence in the dayside magnetosphere

Ultralow-frequency (ULF) waves are observed almost everywhere in the dayside magnetosphere. The mechanism by which these waves are generated and transformed in the dayside magnetosphere is still not understood. Here we report a statistical study of these waves based on magnetic field data from the International Sun-Earth Explorer 1 (ISEE 1) spacecraft. Data from the first traversal of the spacecraft through the entire dayside magnetosphere have been examined to determine the spatial distribution of wave occurrence. Successive 20-min segments of data were transformed to a field-aligned coordinate system. The parallel component was detrended and all three components of the field spectrally analyzed. Wave occurrence was defined by the presence of significant peaks in the power spectra. Wave events were categorized by three wave frequency bands: Pc 3 with T approximately 10-45 s; Pc 4 with T approximately 45-150 s; the short-period part of the Pc 5 wave band with T approximately 150-324 s. Properties of the spectral peaks were then entered into a data base. The data base was next sorted to determine the spatial occurrence pattern for the waves. Our results show that Pc 3 waves most frequently occur just outside synchronous orbit and are approximately centered on local noon. Pc 4 waves have a similar distribution with its peak further out. Pc 5 waves have high occurrence rate at the two flanks of the magnetosphere. Peaks in spectra obtained near the magnetopause are less clearly defined than those deeper in the magnetosphere.

Cao, M.↗

Changes in measured vector magnetic fields when transformed into heliographic coordinates

The changes that occur in measured magnetic fields when they are transformed into a heliographic coordinate system are investigated. To carry out this investigation, measurements of the vector magnetic field of an active region that was observed at 1/3 the solar radius from disk center are taken, and the observed field is transformed into heliographic coordinates. Differences in the calculated potential field that occur when the heliographic normal component of the field is used as the boundary condition rather than the observed line-of-sight component are also examined. The results of this analysis show: (1) that the observed fields of sunspots more closely resemble the generally accepted picture of the distribution of umbral fields if they are displayed in heliographic coordinates; (2) that the differences in the potential calculations are less than 200 G in field strength and 20 deg in field azimuth outside sunspots; and (3) that differences in the two potential calculations in the sunspot areas are no more than 400 G in field strength but range from 60 to 80 deg in field azimuth in localized umbral areas.

Hagyard, M. J.↗

Ocean Turbulence V: Mesoscale Modeling in Level Coordinates. The Effect of Random Nature of Density

The main result of this paper is the derivation of a new expression for the tracer subgrid term in level coordinates S(l) to be employed in O-GCM. The novel feature is the proper account of the random nature of the density field which strongly affects the transformation from isopycnal to level coordinates of the variables of interest, velocity and tracer fields, their correlation functions and ultimately the subgrid terms. In deriving our result we made use of measured properties of vertical ocean turbulence. The major new results are: 1) the new subgrid expression is different from that of the heuristic GM model, 2) u++(tracer)=1/2u+(thickness), where u++ and u+ are the tracer and thickness bolus velocities. In previous models, u++ = u+, 2) the subgrid for a tracer tau is not the same as that for the density rho even when one accounts for the obvious absence of a diffusion term in the latter. The difference stems from a new treatment of the stochastic nature of the density, 3) the mesoscale diffusivity enters both locally and non-locally, as the integral over all z's from the bottom of the ocean to the level z.

Canuto, V. M.↗

A new technique for calculating reentry base heating

The laminar base flow field of a two-dimensional reentry body has been studied by Telenin's method. The flow domain was divided into strips along the x-axis, and the flow variations were represented by Lagrange interpolation polynomials in the transformed vertical coordinate. The complete Navier-Stokes equations were used in the near wake region, and the boundary layer equations were applied elsewhere. The boundary conditions consisted of the flat plate thermal boundary layer in the forebody region and the near wake profile in the downstream region. The resulting two-point boundary value problem of 33 ordinary differential equations was then solved by the multiple shooting method. The detailed flow field and thermal environment in the base region are presented in the form of temperature contours, Mach number contours, velocity vectors, pressure distributions, and heat transfer coefficients on the base surface. The maximum heating rate was found on the centerline, and the two-dimensional stagnation point flow solution was adquate to estimate the maximum heating rate so long as the local Reynolds number could be obtained.

Meng, J. C. S.↗