SEARCH · Search NASA

Results for “divided differences”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 433 records · Page 24

All sources NASA NTRS DOE OSTI

NASA NTRS2009

Hardware Implementation of a Bilateral Subtraction Filter

A bilateral subtraction filter has been implemented as a hardware module in the form of a field-programmable gate array (FPGA). In general, a bilateral subtraction filter is a key subsystem of a high-quality stereoscopic machine vision system that utilizes images that are large and/or dense. Bilateral subtraction filters have been implemented in software on general-purpose computers, but the processing speeds attainable in this way even on computers containing the fastest processors are insufficient for real-time applications. The present FPGA bilateral subtraction filter is intended to accelerate processing to real-time speed and to be a prototype of a link in a stereoscopic-machine- vision processing chain, now under development, that would process large and/or dense images in real time and would be implemented in an FPGA. In terms that are necessarily oversimplified for the sake of brevity, a bilateral subtraction filter is a smoothing, edge-preserving filter for suppressing low-frequency noise. The filter operation amounts to replacing the value for each pixel with a weighted average of the values of that pixel and the neighboring pixels in a predefined neighborhood or window (e.g., a 9 9 window). The filter weights depend partly on pixel values and partly on the window size. The present FPGA implementation of a bilateral subtraction filter utilizes a 9 9 window. This implementation was designed to take advantage of the ability to do many of the component computations in parallel pipelines to enable processing of image data at the rate at which they are generated. The filter can be considered to be divided into the following parts (see figure): a) An image pixel pipeline with a 9 9- pixel window generator, b) An array of processing elements; c) An adder tree; d) A smoothing-and-delaying unit; and e) A subtraction unit. After each 9 9 window is created, the affected pixel data are fed to the processing elements. Each processing element is fed the pixel value for its position in the window as well as the pixel value for the central pixel of the window. The absolute difference between these two pixel values is calculated and used as an address in a lookup table. Each processing element has a lookup table, unique for its position in the window, containing the weight coefficients for the Gaussian function for that position. The pixel value is multiplied by the weight, and the outputs of the processing element are the weight and pixel-value weight product. The products and weights are fed to the adder tree. The sum of the products and the sum of the weights are fed to the divider, which computes the sum of products the sum of weights. The output of the divider is denoted the bilateral smoothed image. The smoothing function is a simple weighted average computed over a 3 3 subwindow centered in the 9 9 window. After smoothing, the image is delayed by an additional amount of time needed to match the processing time for computing the bilateral smoothed image. The bilateral smoothed image is then subtracted from the 3 3 smoothed image to produce the final output. The prototype filter as implemented in a commercially available FPGA processes one pixel per clock cycle. Operation at a clock speed of 66 MHz has been demonstrated, and results of a static timing analysis have been interpreted as suggesting that the clock speed could be increased to as much as 100 MHz.

Huertas, Andres↗

NASA NTRS1999

Gravitational Acceleration Effects on Macrosegregation: Experiment and Computational Modeling

Experiments were performed under terrestrial gravity (1g) and during parabolic flights (10-2 g) to study the solidification and macrosegregation patterns of Al-Cu alloys. Alloys having 2% and 5% Cu were solidified against a chill at two different cooling rates. Microscopic and Electron Microprobe characterization was used to produce microstructural and macrosegregation maps. In all cases positive segregation occurred next to the chill because shrinkage flow, as expected. This positive segregation was higher in the low-g samples, apparently because of the higher heat transfer coefficient. A 2-D computational model was used to explain the experimental results. The continuum formulation was employed to describe the macroscopic transports of mass, energy, and momentum, associated with the solidification phenomena, for a two-phase system. The model considers that liquid flow is driven by thermal and solutal buoyancy, and by solidification shrinkage. The solidification event was divided into two stages. In the first one, the liquid containing freely moving equiaxed grains was described through the relative viscosity concept. In the second stage, when a fixed dendritic network was formed after dendritic coherency, the mushy zone was treated as a porous medium. The macrosegregation maps and the cooling curves obtained during experiments were used for validation of the solidification and segregation model. The model can explain the solidification and macrosegregation patterns and the differences between low- and high-gravity results.

Leon-Torres, J.↗

NASA NTRS1980

Predictions of mineral assemblages in planetary interiors

It is shown that mineral compatibilities in the model system CaO-MgO-Al2O3-SiO2 can be applied to deduce the mineral assemblages expected in planetary interiors and their variation with depth. In general, the available estimates of bulk composition of the terrestrial planets suggest that the terrestrial planets can be divided into two groups based on their predicted mineral assemblages. The terrestrial, Venusian, and lunar bulk compositions are expected to display the following sequence of mineral assemblages with increasing pressure: plagioclase lherzolite, spinel lherzolite, and garnet lherzolite. The sequences expected in Martian and Mercurian are different: spinel-plagioclase wehrlite, spinel lherzolite, and spinel-garnet wehrlite. These assemblages have a major influence on the compositions of liquids produced by melting of these planetary interiors, on the solidus temperatures, and thus on the nature of planetary differentiation and the types of magmas extruded at planetary surfaces.

Stolper, E.↗

NASA NTRS1985

Performance degradation of helicopter rotor in forward flight due to ice

This study addresses the analytical assessment of the degradation in the forward flight performance of the front rotor Boeing Vertol CH47D helicopter in a rime ice natural icing encounter. The front rotor disk was divided into 24 15-deg sections and the local Mach number and angle of attack were evaluated as a function of azimuthal and radial location for a specified flight condition. Profile drag increments were then calculated as a function of azimuthal and radial position for different times of exposure to icing, and the rotor performance was re-evaluated including these drag increments. The results of the analytical prediction method, such as horsepower required to maintain a specific flight condition, as a function of icing time have been generated. The method to illustrate the value of such an approach in assessing performance changes experienced by a helicopter rotor as a result of rime ice accretion is described.

Korkan, K. D.↗

NASA NTRS1986

Computation of transonic separated wing flows using an Euler/Navier-Stokes zonal approach

A computer program called Transonic Navier Stokes (TNS) has been developed which solves the Euler/Navier-Stokes equations around wings using a zonal grid approach. In the present zonal scheme, the physical domain of interest is divided into several subdomains called zones and the governing equations are solved interactively. The advantages of the Zonal Grid approach are as follows: (1) the grid for any subdomain can be generated easily; (2) grids can be, in a sense, adapted to the solution; (3) different equation sets can be used in different zones; and, (4) this approach allows for a convenient data base organization scheme. Using this code, separated flows on a NACA 0012 section wing and on the NASA Ames WING C have been computed. First, the effects of turbulence and artificial dissipation models incorporated into the code are assessed by comparing the TNS results with other CFD codes and experiments. Then a series of flow cases is described where data are available. The computed results, including cases with shock-induced separation, are in good agreement with experimental data. Finally, some futuristic cases are presented to demonstrate the abilities of the code for massively separated cases which do not have experimental data.

Kaynak, Uenver↗

NASA NTRS1999

Global Acceleration of Coronal Mass Ejections

Using the observed relation between speeds of coronal mass ejections (CMEs) near the Sun and in the solar wind, we estimate a global acceleration acting on the CMEs. Our study quantifies the qualitative results of Gosling [1997] and numerical simulations that CMEs at 1 AU with speeds closer to the solar wind. We found a linear relation between the global acceleration and the initial speed of the CMEs and the absolute value of the acceleration is similar to the slow solar wind acceleration. Our study naturally divides CMEs into fast and slow ones, the dividing line being the solar wind speed. Our results have important implications to space weather prediction models which need to incorporate this effect in estimating the CME arrival time at 1 AU. We show that the arrival times of CMEs at 1 AU are drastically different from the zero acceleration case.

Gopalswamy, Nat↗

NASA NTRS1971

Elevation-relief ratio, hypsometric integral, and geomorphic area-altitude analysis.

Mathematical proof establishes identity of hypsometric integral and elevation-relief ratio, two quantitative topographic descriptors developed independently of one another for entirely different purposes. Operationally, values of both measures are in excellent agreement for arbitrarily bounded topographic samples, as well as for low-order fluvial watersheds. By using a point-sampling technique rather than planimetry, elevation-relief ratio (defined as mean elevation minus minimum elevation divided by relief) is calculated manually in about a third of the time required for the hypsometric integral.

Pike, R. J.↗

NASA NTRS1972

Empirical evaluation of pump inlet compliance

Cavitation compliance was determined experimentally from pulsing tests on a number of rocket turbopumps. The primary test data used for this study are those for the Rocketdyne H-1, F-1, and J-2 oxidizer and fuel pumps employed on Saturn vehicles. The study shows that these data can be correlated by a particular form of nondimensionalization, the key feature of which is to divide the operating cavitation number or suction specific speed by its value at head breakdown. An expression is obtained for a best-fit curve for these data. Another set of test data for the Aerojet LR87 and 91 pumps can be correlated by a somewhat different nondimensional pump performance parameter, specifically by relating the cavitation number to its position between the head breakdown point and the point of zero slope of the head coefficient versus cavitation number. Recommendations are given for the estimation of the cavitation compliance for new designs in the Rocketdyne family of pumps.

Ghahremani, F. G.↗

NASA NTRS1993

Initial results from the NASA-Lewis wave rotor experiment

Wave rotors may play a role as topping cycles for jet engines, since by their use, the combustion temperature can be raised without increasing the turbine inlet temperature. In order to design a wave rotor for this, or any other application, knowledge of the loss mechanisms is required, and also how the design parameters affect those losses. At NASA LeRC, a 3-port wave rotor experiment operating on the flow-divider cycle, has been started with the objective of determining the losses. The experimental scheme is a three factor Box-Behnken design, with passage opening time, friction factor, and leakage gap as the factors. Variation of these factors is provided by using two rotors, of different length, two different passage widths for each rotor, and adjustable leakage gap. In the experiment, pressure transducers are mounted on the rotor, and give pressure traces as a function of rotational angle at the entrance and exit of a rotor passage. In addition, pitot rakes monitor the stagnation pressures for each port, and orifice meters measure the mass flows. The results show that leakage losses are very significant in the present experiment, but can be reduced considerably by decreasing the rotor to wall clearance spacing.

Wilson, Jack↗

NASA NTRS1993

Initial results from the NASA Lewis wave rotor experiment

Wilson, Jack↗

NASA NTRS1977

Spacelab output processing system architectural study

Two different system architectures are presented. The two architectures are derived from two different data flows within the Spacelab Output Processing System. The major differences between these system architectures are in the position of the decommutation function (the first architecture performs decommutation in the latter half of the system and the second architecture performs that function in the front end of the system). In order to be examined, the system was divided into five stand-alone subsystems; Work Assembler, Mass Storage System, Output Processor, Peripheral Pool, and Resource Monitor. The work load of each subsystem was estimated independent of the specific devices to be used. The candidate devices were surveyed from a wide sampling of off-the-shelf devices. Analytical expressions were developed to quantify the projected workload in conjunction with typical devices which would adequately handle the subsystem tasks. All of the study efforts were then directed toward preparing performance and cost curves for each architecture subsystem.

Source record↗

NASA NTRSDate not supplied

Implementation of A New Microwave Scattering Database and A Forward Model for Active Microwave Sensors in CRTM

Radiative transfer models are extensively used for the assimilation of satellite observations into NWP models as well as retrieving geophysical products from satellite measurements. CRTM is a community model developed by NOAA JCSDA and widely used for different purposes requiring RT calculations. CRTM requires bulk scattering lookup tables in order to perform all-sky RT calculations. However, the current CRTM lookup tables for microwave frequencies were generated based on the Mie theory by assuming spherical frozen particles. The scattering lookup tables generated using the DDA technique has shown to largely improve the RT scattering calculations in the MW region. This presentation targets (i) the implementation and validation of a DDA database that was originally developed for the ARTS RT model into CRTM, and (ii) developing the CRTM active sensor module that takes advantage of the backscattering coefficients computed using the DDA method. The DDA database only provides single scattering properties of different habits, while CRTM requires bulk scattering properties. The CRTM cloud coefficients were previously generated based on the effective radius for representing the size of the particles. However, effective radius is neither measurable nor provided by the NWP models, thus need to be estimated from other geophysical variables such as water content. Therefore, in addition to calculating the CRTM bulk scattering properties from the DDA single scattering database, the CRTM was also largely modified to use cloud water content (kg.m-3), instead of effective radius, for performing the interpolation over size/mass of the particles. CRTM already requires water content as input, thus no extra variables are required for performing scattering calculations using the new ARTS DDA database. The CRTM scattering modules search for effective radius in cloud coefficient files and will use the cloud water content if the effective radius dimension is not found in the cloud coefficient files. Figure 1 shows the CRTM simulated brightness temperatures computed using different cloud coefficients versus ATMS observed values over Hurricane Irma on September 7, 2017 at 18:00 UTC. We used all the cloud water content values included in ERA5 with default CRTM/DDA habits for water, rain, snow, ice, hail, and graupel. ERA5 does not provide separate water content values for ice, hail, and graupel, thus the ice water content values were divided between ice, hail, and graupel clouds similar to what was explained in the previous section. In channels with a frequency lower than 90 GHz, emission from water and rain clouds can compensate for cloud scattering so that cloud contaminated Tbs are larger than corresponding clear sky Tbs. The DDA simulations for channels 1-7 largely perform better than the Mie simulations. The DDA simulations show a mix of small negative and positive simulated minus observed values, while the Mie results show large negative biases. The weighting functions for some of the ATMS temperature sounding channels (channels 9-15) peak mostly above the clouds, therefore the measured Tbs become less sensitive to clouds so that the results of both Mie and DDA become very similar. The Mie lookup tables generate excessive scattering for channel 16, but not enough scattering for the water vapor channels. In the specific case of Hurricane Maria, the DDA lookup tables do not generate enough scattering for channel 16, but the DDA results are much more consistent with observations for water vapor channels than for channel 16. It should be noted that the results may vary if we use other habits to represent snow, hail, and graupel in the DDA simulations. Although these results clearly show the advantage of the DDA database over the Mie dataset, different error sources such as error in the observations, displacement of clouds in the ERA5 reanalysis, and also lack of convective clouds or in general errors in the input atmospheric and cloud profiles contribute to the differences between the simulated and observed values. Aside from the improvements in the simulations, a major advantage of the new dataset is a large number of habits that can be used to tune the data assimilation systems to perform well in different weather conditions.

Isaac Moradi↗

NASA NTRS2011

An Assessment of Differences Between Cloud Effective Particle Radius Retrievals for Marine Water Clouds from Three MODIS Spectral Bands

The Moderate Resolution Imaging Spectroradiometer (MODIS) cloud product provides three separate 1 km resolution retrievals of cloud particle effective radii (r (sub e)), derived from 1.6, 2.1 and 3.7 micron band observations. In this study, differences among the three size retrievals for maritime water clouds (designated as r (sub e), 1.6 r (sub e), 2.1 and r (sub e),3.7) were systematically investigated through a series of case studies and global analyses. Substantial differences are found between r (sub e),3.7 and r (sub e),2.1 retrievals (delta r (sub e),3.7-2.l), with a strong dependence on cloud regime. The differences are typically small, within +/- 2 micron, over relatively spatially homogeneous coastal stratocumulus cloud regions. However, for trade wind cumulus regimes, r (sub e),3.7 was found to be substantially smaller than r (sub e),2.1, sometimes by more than 10 micron. The correlation of delta r(sub e),3.7-2.1 with key cloud parameters, including the cloud optical thickness (tau), r (sub e) and a cloud horizontal heterogeneity index (H-sigma) derived from 250 m resolution MODIS 0.86 micron band observations, were investigated using one month of MODIS Terra data. It was found that differences among the three r (sub e) retrievals for optically thin clouds (tau <5) are highly variable, ranging from - 15 micron to 10 micron, likely due to the large MODIS retrieval uncertainties when the cloud is thin. The delta r (sub e),3.7-2.1 exhibited a threshold-like dependence on both r (sub e),2.l and H-sigma. The re,3.7 is found to agree reasonably well with re,2.! when re,2.l is smaller than about 15J-Lm, but becomes increasingly smaller than re,2.1 once re,2.! exceeds this size. All three re retrievals showed little dependence when H-sigma < 0.3 (defined as standard deviation divided by the mean for the 250 m pixels within a 1 km pixel retrieval). However, for H-=sigma >0.3, both r (sub e),1.6 and r (sub e),2.1 were seen to increase quickly with H-sigma. On the other hand, r (sub e),3.7 statistics showed little dependence on H-sigma and remained relatively stable over the whole range of H-sigma values. Potential contributing causes to the substantial r (sub e),3.7 and r (sub e),2.1 differences are discussed. In particular, based on both 1-D and 3-D radiative transfer simulations, we have elucidated mechanisms by which cloud heterogeneity and 3-D radiative effects can cause large differences between r (sub e),3.7 and r (sub e),2.l retrievals for highly inhomogeneous clouds. Our results suggest that the contrast in observed delta r (sub e)3.7-2.1 between cloud regimes is correlated with increases in both cloud r (sub e) and H-sigma. We also speculate that in some highly inhomogeneous drizzling clouds, vertical structure induced by drizzle and 3-D radiative effects might operate together to cause dramatic differences between r (sub e),3.7 and r (sub e),2.1 retrievals.

Platnick, Steven↗

NASA NTRS1979

Co-apsidal autonomous terminal rendezvous in Mars orbit

The autonomous (or unmanned) Mars Orbit Rendezvous may be divided into four phases: approach, terminal, station-keeping, and docking. The paper outlines a feasible concept for the terminal rendezvous phase. A co-apsidal terminal rendezvous technique is described which utilizes the advantages of aligning the lines of apsides of two coplanar elliptic orbits of slightly different periods but of the same apoapsis distance. An initial maneuver nulling the inclination between the orbits is followed by another maneuver to align the two lines of apsides. A proper phasing maneuver is then executed at the coapsidal point, and the final rendezvous maneuver is performed at the same coapsidal point to complete the terminal rendezvous transfer. Cases are formulated in order of complexity with a view to illustrate the basic ideas of this technique. Pertinent calculations are also presented.

Tang, C. C. H.↗

NASA NTRS1989

The nature of cometary materials

Because cometary surfaces are likely to be far colder and of a different composition than planetary surfaces, there are some new considerations that must be examined in regards to placing instrumented packages or sample return devices on their surfaces. The qualitative analysis of the problem of attaching hardware to a comet and not being ejected back into space can be divided into two parts. The first problem is to pierce the mantle and obtain access to the icy core. Drilling through the mantle requires that the drilling forces be reacted. Reacting such forces probably requires attachment to the icy core below. Therefore, some kinetic impact piercing device is likely to be required as the first act of attachment. The second problem for a piercing device to overcome is the force produced by the impact kinetic energy that tries to eject the piercing device back into space. The mantle and icy core can absorb some of the impact kinetic energy in the form of fracture formation and friction energy. The energy that is not absorbed in these two ways is stored by the core as elastic deformation of the mantle and icy core. It is concluded that because the cometary materials are almost certainly brittle and the icy core is likely to be self lubricating, the elastic rebound and gas pressure expulsion forces must be counteracted by forces greater than those that may be provided by a piercing device or its capture devices (barbs).

Stephens, James↗

NASA NTRS1991

The compressible mixing layer

The capability of turbulence modeling correctly to handle natural unsteadiness appearing in compressible turbulent flows is investigated. Physical aspects linked to the unsteadiness problem and the role of various flow parameters are analyzed. It is found that unsteady turbulent flows can be simulated by dividing these motions into an 'organized' part for which equations of motion are solved and a remaining 'incoherent' part represented by a turbulence model. Two-equation turbulence models and second-order turbulence models can yield reasonable results. For specific compressible unsteady turbulent flow, graphic presentations of different quantities may reveal complementary physical features. Strong compression zones are observed in rapid flow parts but shocklets do not yet occur.

Vandromme, Dany↗

NASA NTRS2000

Dynamic Load-Balancing for Distributed Heterogeneous Computing of Parallel CFD Problems

The developed methodology is aimed at improving the efficiency of executing block-structured algorithms on parallel, distributed, heterogeneous computers. The basic approach of these algorithms is to divide the flow domain into many sub- domains called blocks, and solve the governing equations over these blocks. Dynamic load balancing problem is defined as the efficient distribution of the blocks among the available processors over a period of several hours of computations. In environments with computers of different architecture, operating systems, CPU speed, memory size, load, and network speed, balancing the loads and managing the communication between processors becomes crucial. Load balancing software tools for mutually dependent parallel processes have been created to efficiently utilize an advanced computation environment and algorithms. These tools are dynamic in nature because of the chances in the computer environment during execution time. More recently, these tools were extended to a second operating system: NT. In this paper, the problems associated with this application will be discussed. Also, the developed algorithms were combined with the load sharing capability of LSF to efficiently utilize workstation clusters for parallel computing. Finally, results will be presented on running a NASA based code ADPAC to demonstrate the developed tools for dynamic load balancing.

Ecer, A.↗

NASA NTRS1973

Major and trace elements in igneous rocks from Apollo 15.

The concentrations of major and trace elements have been determined in igneous rocks from Apollo 15. All materials analyzed have typical depletions of Eu except for minerals separated from sample 15085. Four samples have concentrations of trace elements that are similar to those of KREEP. The samples of mare basalt from Apollo 15 have higher concentrations of FeO, MgO, Mn, and Cr and lower concentrations of CaO, Na2O, K2O, and rare-earth elements (REE) as compared to the samples of mare basalt from Apollos 11, 12, and 14. The samples can be divided into two groups on the basis of their normative compositions. One group is quartz normative and has low concentrations of FeO while the other is olivine normative and has high concentrations of FeO. The trace element data indicate that the samples of olivine normative basalt could be from different portions of a single lava flow.

Helmke, P. A.↗