Search NASA⌕ Search

SEARCH · Search NASA

Results for “MID POINT ALGORITHM”

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 19 records

Derivation and experimental verification of clock synchronization theory

The objective of this work is to validate mathematically derived clock synchronization theories and their associated algorithms through experiment. Two theories are considered, the Interactive Convergence Clock Synchronization Algorithm and the Mid-Point Algorithm. Special clock circuitry was designed and built so that several operating conditions and failure modes (including malicious failures) could be tested. Both theories are shown to predict conservative upper bounds (i.e., measured values of clock skew were always less than the theory prediction). Insight gained during experimentation led to alternative derivations of the theories. These new theories accurately predict the clock system's behavior. It is found that a 100% penalty is paid to tolerate worst case failures. It is also shown that under optimal conditions (with minimum error and no failures) the clock skew can be as much as 3 clock ticks. Clock skew grows to 6 clock ticks when failures are present. Finally, it is concluded that one cannot rely solely on test procedures or theoretical analysis to predict worst case conditions. conditions.

INTERACTIVE CONVERGENCE CLOCK↗

Notes and comments on computational elastoplasticity - Some new models and their numerical implementation

The following topics are discussed in this paper: (1) the basic interactive nature of classical elasto-plasticity and a redefinition of elastic and plastic processes that facilitates numerical calculations, (2) generalized mid-point or end-point algorithms to determine the stress increment in an elastic-plastic solid from a given strain increment, (3) an endochronic (internal time) rate theory of time-independent plasticity which encompasses various multiple-yield-surface theories and nonlinear kinematic hardening theories as its specializations, and (4) comments on finite element and boundary element methods for solving boundary value problems in elasto-plasticity.

Atluri, S. N.↗

An optical systems analysis approach to image resampling

All types of image registration require some type of resampling, either during the registration or as a final step in the registration process. Thus the image(s) must be regridded into a spatially uniform, or angularly uniform, coordinate system with some pre-defined resolution. Frequently the ending resolution is not the resolution at which the data was observed with. The registration algorithm designer and end product user are presented with a multitude of possible resampling methods each of which modify the spatial frequency content of the data in some way. The purpose of this paper is threefold: (1) to show how an imaging system modifies the scene from an end to end optical systems analysis approach, (2) to develop a generalized resampling model, and (3) empirically apply the model to simulated radiometric scene data and tabulate the results. A Hanning windowed sinc interpolator method will be developed based upon the optical characterization of the system. It will be discussed in terms of the effects and limitations of sampling, aliasing, spectral leakage, and computational complexity. Simulated radiometric scene data will be used to demonstrate each of the algorithms. A high resolution scene will be "grown" using a fractal growth algorithm based on mid-point recursion techniques. The result scene data will be convolved with a point spread function representing the optical response. The resultant scene will be convolved with the detection systems response and subsampled to the desired resolution. The resultant data product will be subsequently resampled to the correct grid using the Hanning windowed sinc interpolator and the results and errors tabulated and discussed.

Lyon, Richard G.↗

A discrete momentum-conserving explicit algorithm for multibody dynamics analysis

A discrete momentum-conserving, explicit time integration is presented. The accurate feature and simplicity of the present algorithm are realized by a mid-point implicit formula for integrating the Euler parameters and a second-order discrete momentum-conserving form of the central difference algorithm, respectively. The accuracy and robustness of the algorithm is demonstrated by example problems which exhibit large overall rigid motions under holonomic constraints.

Park, K. C.↗

A discrete momentum-conserving explicit algorithm for rigid body dynamics analysis

A discrete momentum-conserving explicit time integration is presented. The accurate feature and simplicity of the present algorithm are realized by a mid-point implicit formula for integrating the Euler parameters and a second-order discrete momentum-conserving form of the central difference algorithm, respectively. The accuracy and robustness of the algorithm is demonstrated by example problems which exhibit large overall rigid motions under holonomic constraints.

Park, K. C.↗

Staggered solution procedures for multibody dynamics simulation

The numerical solution procedure for multibody dynamics (MBD) systems is termed a staggered MBD solution procedure that solves the generalized coordinates in a separate module from that for the constraint force. This requires a reformulation of the constraint conditions so that the constraint forces can also be integrated in time. A major advantage of such a partitioned solution procedure is that additional analysis capabilities such as active controller and design optimization modules can be easily interfaced without embedding them into a monolithic program. After introducing the basic equations of motion for MBD system in the second section, Section 3 briefly reviews some constraint handling techniques and introduces the staggered stabilized technique for the solution of the constraint forces as independent variables. The numerical direct time integration of the equations of motion is described in Section 4. As accurate damping treatment is important for the dynamics of space structures, we have employed the central difference method and the mid-point form of the trapezoidal rule since they engender no numerical damping. This is in contrast to the current practice in dynamic simulations of ground vehicles by employing a set of backward difference formulas. First, the equations of motion are partitioned according to the translational and the rotational coordinates. This sets the stage for an efficient treatment of the rotational motions via the singularity-free Euler parameters. The resulting partitioned equations of motion are then integrated via a two-stage explicit stabilized algorithm for updating both the translational coordinates and angular velocities. Once the angular velocities are obtained, the angular orientations are updated via the mid-point implicit formula employing the Euler parameters. When the two algorithms, namely, the two-stage explicit algorithm for the generalized coordinates and the implicit staggered procedure for the constraint Lagrange multipliers, are brought together in a staggered manner, they constitute a staggered explicit-implicit procedure which is summarized in Section 5. Section 6 presents some example problems and discussions concerning several salient features of the staggered MBD solution procedure are offered in Section 7.

Park, K. C.↗

A New, More Physically Based Algorithm, for Retrieving Aerosol Properties over Land from MODIS

The MOD Imaging Spectrometer (MODIS) has been successfully retrieving aerosol properties, beginning in early 2000 from Terra and from mid 2002 from Aqua. Over land, the retrieval algorithm makes use of three MODIS channels, in the blue, red and infrared wavelengths. As part of the validation exercises, retrieved spectral aerosol optical thickness (AOT) has been compared via scatterplots against spectral AOT measured by the global Aerosol Robotic NETwork (AERONET). On one hand, global and long term validation looks promising, with two-thirds (average plus and minus one standard deviation) of all points falling between published expected error bars. On the other hand, regression of these points shows a positive y-offset and a slope less than 1.0. For individual regions, such as along the U.S. East Coast, the offset and slope are even worse. Here, we introduce an overhaul of the algorithm for retrieving aerosol properties over land. Some well-known weaknesses in the current aerosol retrieval from MODIS include: a) rigid assumptions about the underlying surface reflectance, b) limited aerosol models to choose from, c) simplified (scalar) radiative transfer (RT) calculations used to simulate satellite observations, and d) assumption that aerosol is transparent in the infrared channel. The new algorithm attempts to address all four problems: a) The new algorithm will include surface type information, instead of fixed ratios of the reflectance in the visible channels to the mid-IR reflectance. b) It will include updated aerosol optical properties to reflect the growing aerosol retrieved from eight-plus years of AERONE". operation. c) The effects of polarization will be including using vector RT calculations. d) Most importantly, the new algorithm does not assume that aerosol is transparent in the infrared channel. It will be an inversion of reflectance observed in the three channels (blue, red, and infrared), rather than iterative single channel retrievals. Thus, this new formulation of the MODIS aerosol retrieval over land includes more physically based surface, aerosol and radiative transfer with fewer potentially erroneous assumptions.

Levy, Robert C.↗

High Resolution Observations of Drop Size Distribution for GPM Ground Validation

During the Mid-latitude Convective Cloud Experiment (MC3E), NASA's GPM GV Disdrometer and Radar Observations of Precipitation (DROP) Facility deployed an array of disdrometers and rain gauges in northern Oklahoma to sample, with high resolution, the drop size distribution (DSD) for use in development of precipitation retrieval algorithms for the GPM core satellite. The DROP Facility instruments deployed during MC3E consisted of 16 autonomous Parsivel units, 5 two-dimensional video disdrometers (2DVD), a vertical pointing K-band radar and 32 tipping bucket rain gauges. There were several rainfall events during MC3E in which rain drops exceeding 5 mm in diameter were recorded. During a convective rainfall event one of the 2DVDs measured a large rain drop of 8 mm in diameter falling from a convective cell under which several other disdrometers revealed a broad DSD spectrum across this cell and median drop diameter of 3.1 mm. The NPOL radar, which was scanning over the disdrometer array in high resolution RHI mode (every 40 sec), during a stratiform rainfall event revealed a 1 km thick bright band with locally thicker segments that produced rainfall rates of 6-10 mm/hr and indications of DSD modulations as enhanced concentrations of mid-sized rain drops (1.5-3.0 mm) descended below cloud base. These observations are key in characterizing the DSD variability within the 4 km footprint of the GPM Dual-wavelength Precipitation Radar that will provide a global map of precipitation as well as DSD.

Gatlin, Patrick N.↗

Accounting for the photochemical variation in stratospheric NO2 in the SAGE III/ISS solar occultation retrieval

The Stratospheric Aerosol and Gas Experiment (SAGE) III has been operating on the International Space Station (ISS) since mid-2017. Nitrogen dioxide (NO2) number density profiles are routinely retrieved from SAGE III/ISS solar occultation measurements in the middle atmosphere. Although NO2 density varies throughout the day due to photochemistry, the standard SAGE NO2 retrieval algorithm neglects these variations along the instrument’s line of sight by assuming that the number density has a constant gradient within a given vertical layer of the atmosphere. This assumption will result in a retrieval bias for a species like NO2 that changes rapidly across the terminator. In this work we account for diurnal variations in retrievals of NO2 from the SAGE III/ISS measurements, and we determine the impact of this algorithm improvement on the resulting NO2 number densities. The first step in applying the diurnal correction is to use publicly available SAGE III/ISS products to convert the retrieved number density profiles to optical depth profiles. The retrieval is then re-performed with a new matrix that applies photochemical scale factors for each point along the line of sight according to the changing solar zenith angle. In general NO2 that is retrieved by accounting for these diurnal variations is more than 10% lower than the standard algorithm below 30 km. This effect is greatest in winter at high latitudes and generally greater for sunrise occultations than sunset. Comparisons with coincident profiles from the Optical Spectrograph and InfraRed Imager System (OSIRIS) show that NO2 from SAGE III/ISS is generally biased high; however the agreement improves by up to 20% in the mid-stratosphere when diurnal variations are accounted for in the retrieval. We conclude that diurnal variations along the SAGE III/ISS line of sight are an important term to consider for NO2 analyses at altitudes below 30 km.

Stratospheric NO2↗

High Resolution Observations of Drop Size Distribution for GPM Ground Validation

During the Mid-latitude Continental Convective Cloud Experiment (MC3E), NASA's GPM GV Disdrometer and Radar Observations of Precipitation (DROP) Facility deployed an array of disdrometers and rain gauges in northern Oklahoma to sample, with high resolution, the drop size distribution for use in development of precipitation retrieval algorithms for the GPM core satellites. The DROP Facility instruments deployed during MC3E consisted of 16 autonomous Parsivel units, 5 two-dimensional video disdrometers (2dvds), a vertically pointing K band radar, and 32 tipping bucket rain gauges. There were several rainfall events during MC3E in which rain drops exceeding 6 mm in diameter were recorded. The disdrometer array revealed large rain drops with diameters exceeding 6 mm and 8 mm during two separate stratiform and convective rainfall events, respectively. The NPOL radar, which was scanning in high resolution RHI mode (every 40 sec) over the disdrometer array during the stratiform event, indicated a 1 km thick bright band with a differential reflectivity column of 2-3 dB extending below the melting layer to the surface where the large drops were recorded by the 2dvds. These large drops are important for GPM since they can have a great impact upon satellite precipitation retrieval, especially near the ground and below heavy convective rainfall cores where satellites have had problems depicting the rainfall.

Gatlin, Patrick N.↗

MODIS Status

The Moderate Resolution Imaging Spectroradiometer (MODIS) was sent into orbit on the Earth Observing System (EOS) Terra spacecraft in December 1999. The MODIS began taking observations on February 24, 2000. Instrument checkout and characterization progressed to the point in the Fall of 2000 along with the checkout of algorithms so that approximately 40 products for land, ocean and atmosphere studies on global and regional scales were and are now being produced systematically. These products are now designated "beta", products indicating that they are still being examined and validated, but in most cases are amenable to examination by the scientific and applications community to further assess their utility and provide feedback to the validation process and the MODIS Science Team. It is expected that many products will progress to "provisionally useful" or fully "validated" status by mid-year 2001. one goal is to produce a systematically processed data set by the end of 2001 from MODIS extends from November 2000 through October 2001. Overall the MODIS instrument and the associated data processing systems are performing well. The many examples of MODIS observations indicate that the prospect for highly useful and exciting studies of the Earth-atmosphere system using MODIS data looks very good.

Salomonson, Vincent V.↗

Acoustic Calibration of the Exterior Effects Room at the NASA Langley Research Center

The Exterior Effects Room (EER) at the NASA Langley Research Center is a 39-seat auditorium built for psychoacoustic studies of aircraft community noise. The original reproduction system employed monaural playback and hence lacked sound localization capability. In an effort to more closely recreate field test conditions, a significant upgrade was undertaken to allow simulation of a three-dimensional audio and visual environment. The 3D audio system consists of 27 mid and high frequency satellite speakers and 4 subwoofers, driven by a real-time audio server running an implementation of Vector Base Amplitude Panning. The audio server is part of a larger simulation system, which controls the audio and visual presentation of recorded and synthesized aircraft flyovers. The focus of this work is on the calibration of the 3D audio system, including gains used in the amplitude panning algorithm, speaker equalization, and absolute gain control. Because the speakers are installed in an irregularly shaped room, the speaker equalization includes time delay and gain compensation due to different mounting distances from the focal point, filtering for color compensation due to different installations (half space, corner, baffled/unbaffled), and cross-over filtering.

Faller, Kenneth J., II↗

Global electric field determination in the Earth's outer magnetosphere using energetic charged particles

Although many properties of the Earth's magnetosphere have been measured and quantified in the past 30 years since it was discovered, one fundamental measurement (for zeroth order MHD equilibrium) has been made infrequently and with poor spatial coverage - the global electric field. This oversight is due in part to the neglect of theorists. However, there is renewed interest in the convection electric field because it is now realized to be central to many magnetospheric processes, including the global MHD equilibrium, reconnection rates, Region 2 Birkeland currents, magnetosphere ionosphere coupling, ring current and radiation belt transport, substorm injections, and several acceleration mechanisms. Unfortunately the standard experimental methods have not been able to synthesize a global field (excepting the pioneering work of McIlwain's geostationary models) and we are left with an overly simplistic theoretical field, the Volland-Stern electric field model. Single point measurements of the plasmapause were used to infer the appropriate amplitudes of this model, parameterized by K(sub p). Although this result was never intended to be the definitive electric field model, it has gone nearly unchanged for 20 years. The analysis of current data sets requires a great deal more accuracy than can be provided by the Volland-Stern model. The variability of electric field shielding has not been properly addressed although effects of penetrating magnetospheric electric fields has been seen in mid-and low-latitude ionospheric data sets. The growing interest in substorm dynamics also requires a much better assessment of the electric fields responsible for particle injections. Thus we proposed and developed algorithms for extracting electric fields from particle data taken in the Earth's magnetosphere. As a test of the effectiveness of these new techniques, we analyzed data taken by the AMPTE/CCE spacecraft in equatorial orbit from 1984 to 1989.

Eastman, Timothy E.↗

Global electric field determination in the Earth's outer magnetosphere using charged particles

Although many properties of the Earth's magnetosphere have been measured and quantified in the past 30 years since it was discovered, one fundamental (for a zeroeth order magnetohydrodynamic (MHD) equilibrium) measurement was made infrequently and with poor spatial coverage: the global electric field. This oversight is in part due to the difficulty of measuring a plasma electric field, and in part due to the difficulty of measuring a plasma electric field, and in part due to the neglect of theorists. However, there is renewed interest in the convection electric field, since it has been realized that it is vital for understanding many aspects of the magnetosphere: the global MHD equilibrium, reconnection rates, Region 2 Birkeland currents, magnetosphere-ionosphere coupling, ring current and radiation belt transport, substorm injections, acceleration mechanisms, etc. Unfortunately the standard experimental methods have not been able to synthesize a global field (excepting the pioneering work of McIlwain's geostationary models), and we are left with an overly simplistic theoretical field, the Volland-Stern electric field mode. Again, single point measurements of the plasma pause were used to infer the appropriate amplitudes of the model, parameterized by Kp (Maynard & Chen, JGR 1975). Although this result was never intended to be the definitive electric field model, it has gone nearly unchanged for 15 years. However, the data sets being taken today require a great deal more accuracy than can be provided by the Volland-Stern model. Nor has the variability of the electric field shielding been properly addressed, although effects of penetrating magnetospheric electric fields has been seen in mid- and low-latitude ionospheric data sets. The growing interests in substorm dynamics also requires a much better assessment of the electric fields responsible for particle injections. Thus, we proposed and developed algorithms for extracting electric fields from particle data taken in the earth's magnetosphere. As a test of the effectiveness of these techniques, we analyzed data taken by the AMPTE/CCE spacecraft in equatorial orbit between 1984-1988. Some analytic tools had to be developed before construction of computer algorithms, and they are discussed.

Eastman, T.↗

Comparison of Aircraft Models and Integration Schemes for Interval Management in the TRACON

Reusable models of common elements for communication, computation, decision and control in air traffic management are necessary in order to enable simulation, analysis and assurance of emergent properties, such as safety and stability, for a given operational concept. Uncertainties due to faults, such as dropped messages, along with non-linearities and sensor noise are an integral part of these models, and impact emergent system behavior. Flight control algorithms designed using a linearized version of the flight mechanics will exhibit error due to model uncertainty, and may not be stable outside a neighborhood of the given point of linearization. Moreover, the communication mechanism by which the sensed state of an aircraft is fed back to a flight control system (such as an ADS-B message) impacts the overall system behavior; both due to sensor noise as well as dropped messages (vacant samples). Additionally simulation of the flight controller system can exhibit further numerical instability, due to selection of the integration scheme and approximations made in the flight dynamics. We examine the theoretical and numerical stability of a speed controller under the Euler and Runge-Kutta schemes of integration, for the Maintain phase for a Mid-Term (2035-2045) Interval Management (IM) Operational Concept for descent and landing operations. We model uncertainties in communication due to missed ADS-B messages by vacant samples in the integration schemes, and compare the emergent behavior of the system, in terms of stability, via the boundedness of the final system state. Any bound on the errors incurred by these uncertainties will play an essential part in a composable assurance argument required for real-time, flight-deck guidance and control systems,. Thus, we believe that the creation of reusable models, which possess property guarantees, such as safety and stability, is an innovative and essential requirement to assessing the emergent properties of novel airspace concepts of operation.

Neogi, Natasha↗

Effect of MODIS Terra Radiometric Calibration Improvements on Collection 6 Deep Blue Aerosol Products: Validation and Terra/Aqua Consistency

The Deep Blue (DB) algorithm's primary data product is midvisible aerosol optical depth (AOD). DB applied to Moderate Resolution Imaging Spectroradiometer (MODIS) measurements provides a data record since early 2000 for MODIS Terra and mid-2002 for MODIS Aqua. In the previous data version (Collection 5, C5), DB production from Terra was halted in 2007 due to sensor degradation; the new Collection 6 (C6) has both improved science algorithms and sensor radiometric calibration. This includes additional calibration corrections developed by the Ocean Biology Processing Group to address MODIS Terra's gain, polarization sensitivity, and detector response versus scan angle, meaning DB can now be applied to the whole Terra record. Through validation with Aerosol Robotic Network (AERONET) data, it is shown that the C6 DB Terra AOD quality is stable throughout the mission to date. Compared to the C5 calibration, in recent years the RMS error compared to AERONET is smaller by approximately 0.04 over bright (e.g., desert) and approximately 0.01-0.02 over darker (e.g., vegetated) land surfaces, and the fraction of points in agreement with AERONET within expected retrieval uncertainty higher by approximately 10% and approximately 5%, respectively. Comparisons to the Aqua C6 time series reveal a high level of correspondence between the two MODIS DB data records, with a small positive (Terra-Aqua) average AOD offset <0.01. The analysis demonstrates both the efficacy of the new radiometric calibration efforts and that the C6 MODIS Terra DB AOD data remain stable (to better than 0.01 AOD) throughout the mission to date, suitable for quantitative scientific analyses.

Collection 6 Deep Blue↗

Preparing the JWST Observatory for Science Observations

The James Webb Space Telescope (JWST) is a large (6.5 m) near-and mid-infrared telescope scheduled for launch in 2021. JWST will be used to explore fundamental questions in astrophysics and planetary science, including the evolution of galaxies, the first light of stars, the formation of stars and planets as well as the characterization of exoplanets. To achieve these scientific goals, NASA and its international partners, the Canadian Space Agency (CSA) and the European Space Agency (ESA), have developed four science instruments for JWST and have prepared detailed commissioning plans for each. The NASA Project Science team has also outlined activities to characterize the performance of the Observatory as a whole. These activities include: 1) monitoring of the instrument and telescope cooldown for contamination mitigation; 2) measurements of straylight and other backgrounds for validation of the JWST stray light models as well as instrument background subtraction algorithms; 3) characterization of optical thermal distortion between its hot and cold telescope pointings within the observable field of regard and 4) trending of observatory performance parameters. Each of these activities has been planned in close collaboration with the observatory wavefront, commissioning planning, and science instrument teams. Here we present the plans for these activities, the expected results and how they will impact future guest observer (GO) proposals and JWST science community

Erin C Smith↗

P69 Using the NASA-Unified WRF to Assess the Impacts of Real-Time Vegetation on Simulations of Severe Weather

Since June 2010, the NASA Short-term Prediction Research and Transition (SPoRT; Goodman et al. 2004; Darden et al. 2010; Stano et al. 2012; Fuell et al. 2012) Center has been generating a real-time Normalized Difference Vegetation Index (NDVI) and corresponding Green Vegetation Fraction (GVF) composite based on reflectances from NASA s Moderate Resolution Imaging Spectroradiometer (MODIS) instrument. This dataset is generated at 0.01 resolution across the Continental United States (CONUS), and updated daily. The goal of producing such a vegetation dataset is to improve over the default climatological GVF dataset in land surface and numerical weather prediction models, in order to have better simulations of heat and moisture exchange between the land surface and the planetary boundary layer. Details on the SPoRT/MODIS vegetation composite algorithm are presented in Case et al. (2011). Vegetation indices such as GVF and Leaf Area Index (LAI) are used by land surface models (LSMs) to represent the horizontal and vertical density of plant vegetation (Gutman and Ignatov 1998), in order to calculate transpiration, interception and radiative shading. Both of these indices are related to the NDVI; however, there is an inherent ambiguity in determining GVF and LAI simultaneously from NDVI, as described in Gutman and Ignatov (1998). One practice is to specify the LAI while allowing the GVF to vary both spatially and temporally, as is done in the Noah LSM (Chen and Dudhia 2001; Ek et al. 2003). Operational versions of Noah within several of the National Centers for Environmental Prediction (NCEP) global and regional modeling systems hold the LAI fixed, while the GVF varies according to a global monthly climatology. This GVF climatology was derived from NDVI data on the NOAA Advanced Very High Resolution Radiometer (AVHRR) polar orbiting satellite, using information from 1985 to 1991 (Gutman and Ignatov 1998; Jiang et al. 2010). Representing data at the mid-point of every month, the climatological dataset is on a grid with 0.144 (~16 km) spatial resolution and is distributed with the community WRF model (Ek et al. 2003; Jiang et al. 2010; Skamarock et al. 2008).

Case, Jonathan L.↗