Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “smoothing”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Potential Mars Surveyor 2001 Landing Sites: Low-Elevation Cratered "Highlands" in Central and Eastern Sinus Meridiani and Near Amenthes Fossae

The main scientific goal for the Mars Surveyor Program 2001 (MSP 01) landed mission is to collect and characterize 91 rock and 13 soil core samples using an integrated instrument suite onboard the Athena Rover. If possible, these samples will be retrieved and returned to Earth via a MSP 05 or MSP 07 mission. Preliminary engineering constraints for the MSP 01 landing site call for a location that lies between 15 S and 30 N, and below about 2 km elevation (based on Viking-era topography). Desirable landing sites for MSP 01 are to be located "in the ancient highlands where the environmental conditions may have been favorable to the preservation of evidence of possible prebiotic or biotic processes including the emergence (and, potentially, the persistence) of life". We interpret this to mean that the desirable sites include those that have evidence of aqueous sediments that might have been deposited during the Noachian and/or Hesperian Epochs of Mars' history. In addition to the search for subaqueous sedimentary deposits, we took into consideration the fact that the rover, Athena, will need to be able to access these materials. Thus, a site where aeolian deflation has occurred might be desirable because it might expose, in situ, layered sedimentary deposits. Deflated areas, of course, might include potential landing hazards in the form of meterscale buttes and mesas (e.g., Christmas Lake Valley, OR), thus careful study of such sites with high resolution images will be required before a decision is made to land. We have been examining three regions that have potential to be considered for MSP 01 landing sites. This work is based on Viking (VIS, IRTM) and Phobos 2 (Termoskan) observations and should be regarded as preliminary because we believe that the final site selection should also be based upon analysis of Mars Global Surveyor observations that help constrain mineralogy (TES) and local geomorphology (MOC, MOLA). (1) Eastern Sinus Meridiani Region (proposed by K. S. Edgett) Sinus Meridiani is a persistent low-albedo (< 0. 16) region on the martian equator that has been recognized for about 400 years. All of Sinus Meridiani is below the 2 km elevation constraint for MSP 01. The region includes cratered highlands, valley networks, aeolian dunes, and possible aqueous sedimentary deposits. The landing site study region is located in the eastern portion of Sinus Meridiani. It is bounded by latitudes 10 S to 2 N, longitudes 355 W to 345 W, and the elevations are mostly 1-2 km. The center of this area contains a medium-albedo (0.19-0.21), relatively smooth-surfaced deposit that was suggested by Rice to be a lacustrine deposit. Several potential landing areas can be suggested within this region. Based on Viking and Phobos 2 data, the favored sites so far are centered at 7.6 S, 346.9 W and 0.8 S, 349 W. Thermal inertias are 3.2-7.0 x 10(exp -3) cal /sq cm s(exp -0.5)/K; and rock abundances are around 2-6%. The site at 7.6 S, 346.9 W is at the southern end of the smooth, medium-albedo unit that might have a lacustrine origin. At this location, numerous channels appear to have drained toward the smooth unit. Viking images from orbit 747A show this area at about 15 m/pixel ground resolution. The images reveal that aeolian deflation has occurred along the deposit's margins. Bright (i.e., albedo >= 0.21) aeolian dunes are present on the channel floors and in some of the depressions on the smooth unit. The bright, apparently active dunes might consist of material (perhaps lakedeposited sands) that has been eroded from the smooth unit. The site at approximately 0.8 S, 349 W is selected because it offers an opportunity to solve a long-standing puzzle about Mars remote sensing. There are three main "color" units on Mars: "dark red, dark gray, and bright red". This landing site would allow the Athena rover an opportunity to investigate all three materials within close proximity (the best place on Mars to do so). There are no high resolution (better than 100 m/pixel) Viking or Mariner images of this site. (2) Central Sinus Meridiani Region (proposed by K. S. Edgett and T. J. Parker) Central Sinus Meridiani is characterized by two types of surfaces [4]. One is like typical martian cratered highlands elsewhere- there are old valley networks and old impact craters. The other is relatively smooth and flat. These two units are in contact around 3.1 S between 5 W and 4 E longitudes. Valley networks- including one at 6'S, 358 W that rivals the Grand Canyon of Arizona- once drained toward the smooth unit. Edgett and Parker [4] proposed that the smooth unit might consist of sediments laid down in a large Noachian-aged sea/ocean that would have covered much of the northern hemisphere. Schultz and Lutz [I I I suggested that it is a paleopolar layered deposit. Regardless, the smooth unit where it contacts the cratered terrain would make an excellent site for Athena rover to investigate. The site is best seen in Viking high resolution images from orbits 408B (about 30 m/pixel) and 746A (about 12 m/pixel). These images suggest that aeolian deflation has occurred along the margin of the smooth unit, and this deflation has exposed horizontal larrs of material. The elevation is about 0.5 km; thermal inertias are 6.5-8.0 x 10(exp -3) cal /sq cm s(exp -0.5) / K; rock abundances are 2-4%; and the surface is probably sandy with dark drifts and ripples but almost no actual dunes. We suggest a landing around 3.2 S, 3.0 W would test the aqueous sediment hypothesis and provide a potentially smooth surface on which to land. (3) Amenthes Fossae Region (proposed by S. N. Huntwork and K. S. Edgett) The Amenthes Fossae are a series of graben/fissures that are circumferential to the southeast side of Isidis Planitia. These fissures cross a variety of ancient, heavily cratered Noachian terrain and younger, Hesperian and Amazonian terrain. We focused our search on a region 0-15 N, 250 - 270 W. Elevations are -0.5 to 2 km. Depending upon whether Isidis Planitia was ever a water-rich environment, this region might have been influenced by aqueous sedimentation. Valley networks are common, and they drained toward the north and northwest. We focused our work on a set of Viking orbiter high-resolution images, 719A 1-48. These have resolutions 16-24 rri/pixel. We examined images 20-23, centered at 2 N, 258 W. This site, on the plains just southwest of a 42 km-diameter crater, includes a valley network channel, a relatively young crater ejecta deposit, a few buttes composed of presumably ancient, Noachian bedrock, and a "plains" unit. The plains might make an ideal landing surface, except for the presence of some fine-scale ridges (oriented approximately N-S). The ridges are probably yardangs, thus this site offers a place where aeolian deflation has probably exposed some of the layered rock units that comprise the "plains". Thermal inertias are 7.9-8.3 x 10(exp -3) cal/sq cm s(exp -0.5) /K and rock abundances are 10-15%. A rover traverse might include the opportunity to go down to the floor (and sample along the walls) of the valley network channel at 2 N 258.1 W.

Edgett, K. S.↗
NASA NTRS1994

Accurate interlaminar stress recovery from finite element analysis

The accuracy and robustness of a two-dimensional smoothing methodology is examined for the problem of recovering accurate interlaminar shear stress distributions in laminated composite and sandwich plates. The smoothing methodology is based on a variational formulation which combines discrete least-squares and penalty-constraint functionals in a single variational form. The smoothing analysis utilizes optimal strains computed at discrete locations in a finite element analysis. These discrete strain data are smoothed with a smoothing element discretization, producing superior accuracy strains and their first gradients. The approach enables the resulting smooth strain field to be practically C1-continuous throughout the domain of smoothing, exhibiting superconvergent properties of the smoothed quantity. The continuous strain gradients are also obtained directly from the solution. The recovered strain gradients are subsequently employed in the integration o equilibrium equations to obtain accurate interlaminar shear stresses. The problem is a simply-supported rectangular plate under a doubly sinusoidal load. The problem has an exact analytic solution which serves as a measure of goodness of the recovered interlaminar shear stresses. The method has the versatility of being applicable to the analysis of rather general and complex structures built of distinct components and materials, such as found in aircraft design. For these types of structures, the smoothing is achieved with 'patches', each patch covering the domain in which the smoothed quantity is physically continuous.

Tessler, Alexander↗
DOE OSTI2022

Measurement of shock roughness due to phase plate speckle imprinting relevant for x-ray diffraction experiments on 3rd and 4th generation light sources

Laser-shock compression experiments at 3rd and 4th generation light sources generally employ phase plates, which are inserted into the beamline to achieve a repeatable intensity distribution at the focal plane. Here, the laser intensity profile is characterized by a high-contrast, high-frequency laser speckle. Without sufficient smoothing, these laser non-uniformities can translate to a significant pressure distribution within the sample layer and can affect data interpretation in x-ray diffraction experiments. Here, we use a combination of one- and two-dimensional velocity interferometry to directly measure the extent to which spatial frequencies within the laser focal spot intensity pattern are smoothed out during propagation within the laser plasma and a polyimide ablator. We find that the use of thicker polyimide layers results in spatially smoother shock fronts, with the greatest degree of smoothing associated with the highest spatial frequencies. Focal spots with the smallest initial speckle separation produce the most rapid smoothing. Laser systems that employ smoothing by spectral dispersion techniques to rapidly modulate the focal plane intensity distribution are shown to be the most effective ones in producing a spatially smooth shock front. We show that a simple transport model combined with the known polyimide Hugoniot adequately describes the extent of shock smoothness as a function of polyimide thickness. Finally, our results provide a description of spatial structure smoothing across a shock front, which can be used to design targets on x-ray free electron laser facilities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1991

Investigation of the influence of a step change in surface roughness on turbulent heat transfer

The use is studied of smooth heat flux gages on the otherwise very rough SSME fuel pump turbine blades. To gain insights into behavior of such installations, fluid mechanics and heat transfer data were collected and are reported for a turbulent boundary layer over a surface with a step change from a rough surface to a smooth surface. The first 0.9 m length of the flat plate test surface was roughened with 1.27 mm hemispheres in a staggered, uniform array spaced 2 base diameters apart. The remaining 1.5 m length was smooth. The effect of the alignment of the smooth surface with respect to the rough surface was also studied by conducting experiments with the smooth surface aligned with the bases or alternatively with the crests of the roughness elements. Stanton number distributions, skin friction distributions, and boundary layer profiles of temperature and velocity are reported and are compared to previous data for both all rough and all smooth wall cases. The experiments show that the step change from rough to smooth has a dramatic effect on the convective heat transfer. It is concluded that use of smooth heat flux gages on otherwise rough surfaces could cause large errors.

Taylor, Robert P.↗
NASA NTRS1994

Topological analysis of the CfA redshift survey

We study the topology of large-scale structure in the Center for Astrophysics Redshift Survey, which now includes approximately 12,000 galaxies with limiting magnitude m(sub B) is less than or equal to 15.5. The dense sampling and large volume of this survey allow us to compute the topology on smoothing scales from 6 to 20/h Mpc; we thus examine the topology of structure in both 'nonlinear' and 'linear' regimes. On smoothing scales less than or equal to 10/h Mpc this sample has 3 times the number of resolution elements of samples examined in previous studies. Isodensity surface of the smoothed galaxy density field demonstrate that coherent high-density structures and large voids dominate the galaxy distribution. We compute the genus-threshold density relation for isodensity surfaces of the CfA survey. To quantify phase correlation in these data, we compare the CfA genus with the genus of realizations of Gaussian random fields with the power spectrum measured for the CfA survey. On scales less than or equal to 10/h Mpc the observed genus amplitude is smaller than random phase (96% confidence level). This decrement reflects the degree of phase coherence in the observed galaxy distribution. In other words the genus amplitude on these scales is not good measure of the power spectrum slope. On scales greater than 10/h Mpc, where the galaxy distribution is rougly in the 'linear' regime, the genus ampitude is consistent with the random phase amplitude. The shape of the genus curve reflects the strong coherence in the observed structure; the observed genus curve appears broader than random phase (94% confidence level for smoothing scales less than or equal to 10/h Mpc) because the topolgoy is spongelike over a very large range of density threshold. This departre from random phase consistent with a distribution like a filamentary net of 'walls with holes.' On smoothing scales approaching approximately 20/h Mpc the shape of the CfA genus curve is consistent with random phase. There is very weak evidence for a shift of the genus toward a 'bubble-like' topology. To test cosmological models, we compute the genus for mock CfA surveys drawn from large (L greater than or approximately 400/h Mpc) N-body simulations of three variants of the cold dark matter (CDM) cosmogony. The genus amplitude of the 'standard' CDM model (omega h = 0.5, b = 1.5) differs from the observations (96% confidence level) on smoothing scales is less than or approximately 10/h Mpc. An open CDM model (omega h = 0.2) and a CDM model with nonzero cosmological constant (omega h = 0.24, lambda (sub 0) = 0.6) are consistent with the observed genus amplitude over the full range of smoothing scales. All of these models fail (97% confidence level) to match the broadness of the observed genus curve on smoothing scales is less than or equal to 10/h Mpc.

Vogeley, Michael S.↗
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↗
DOE OSTI2026

Overlapping Schwarz Methods Are Not Anisotropy‐Robust Multigrid Smoothers

We analyze overlapping multiplicative Schwarz methods as smoothers in the geometric multigrid solution of two-dimensional anisotropic diffusion problems. For diffusion equations, it is well known that the smoothing properties of point-wise smoothers, such as Gauss Seidel, rapidly deteriorate as the strength of anisotropy increases. On the other hand, global smoothers based on line smoothing are known to generally provide good smoothing for diffusion problems, independent of the anisotropy strength. Here, a natural question is whether global methods are really necessary to achieve good smoothing in such problems, or whether it can be obtained with locally overlapping block smoothers using sufficiently large blocks and overlap. Through local Fourier analysis and careful numerical experimentation, we show that global methods are indeed necessary to achieve anisotropy-robust smoothing. Specifically, for any fixed block size bounded sufficiently far away from the global domain size, we find that the smoothing properties of overlapping multiplicative Schwarz rapidly deteriorate with increasing anisotropy, irrespective of the amount of overlap between blocks. Moreover, our results indicate that anisotropy-robust smoothing requires blocks of diameter 𝒪⁡(𝜖 −1/2 ) for anisotropy ratio 𝜖 ∈(0,1] .

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Illuminating dark matter halo density profiles without subhaloes

Cold dark matter haloes consist of a relatively smooth dark matter component as well as a system of bound subhaloes. It is the prevailing practice to include all mass, including mass in subhaloes, in studies of halo density profiles in simulations. However, often in observational studies satellites are treated as having their own distinct dark matter density profiles in addition to the profile of the host. This difference can make comparisons between theoretical and observed results difficult. In this work, we investigate density profiles of the smooth components of host haloes by excluding mass contained within subhaloes. We find that the density profiles of the smooth halo component (without subhaloes) differ substantially from the conventional halo density profile, declining more rapidly at large radii. We also find that concentrations derived from smooth density profiles exhibit less scatter at fixed mass and a weaker mass dependence than standard concentrations. Both smooth and standard halo profiles can be described by a generalized Einasto profile, an Einasto profile with a modified central slope, with smaller residuals than either a Navarro–Frenk–White or Einasto profile. Furthermore, these results hold for both Milky Way-mass and cluster-mass haloes. This new characterization of smooth halo profiles can be useful for many analyses, such as lensing and dark matter annihilation, in which the smooth and clumpy components of a halo should be accounted for separately.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2012

The Re-Analysis of Ozone Profile Data from a 41-Year Series of SBUV Instruments

In this study we present the validation of ozone profiles from a number of Solar Back Scattered Ultra Violet (SBUV) and SBUV/2 instruments that were recently reprocessed using an updated (Version 8.6) algorithm. The SBUV dataset provides the longest available record of global ozone profiles, spanning a 41-year period from 1970 to 2011 (except a 5-year gap in the 1970s) and includes ozone profile records obtained from the Nimbus-4 BUV and Nimbus-7 SBUV instruments, and a series of SBUV(/2) instruments launched on NOAA operational satellites (NOAA 09, 11, 14, 16, 17, 18, 19). Although modifications in instrument design were made in the evolution from the BUV instrument to the modern SBUV(/2) model, the basic principles of the measurement technique and retrieval algorithm remain the same. The long term SBUV data record allows us to create a consistent, calibrated dataset of ozone profiles that can be used for climate studies and trend analyses. In particular, we focus on estimating the various sources of error in the SBUV profile ozone retrievals using independent observations and analysis of the algorithm itself. For the first time we include in the metadata a quantitative estimate of the smoothing error, defined as the error due to profile variability that the SBUV observing system cannot inherently measure. The magnitude of the smoothing error varies with altitude, latitude, season and solar zenith angle. Between 10 and 1 hPa the smoothing errors for the SBUV monthly zonal mean retrievals are of the order of 1 %, but start to increase above and below this layer. The largest smoothing errors, as large as 15-20%, were detected in in the troposphere. The SBUV averaging kernels, provided with the ozone profiles in version 8.6, help to eliminate the smoothing effect when comparing the SBUV profiles with high vertical resolution measurements, and make it convenient to use the SBUV ozone profiles for data assimilation and model validation purposes. The smoothing error can also be minimized by combining layers of data, and we will discuss recommendations for this approach as well. The SBUV ozone profiles have been intensively validated against satellite profile measurements obtained from the Microwave Limb Sounders (MLS) (on board the UARS and AURA satellites), Stratospheric Aerosol and Gas Experiment (SAGE) and Michelson Interferometer for Passive Atmospheric Sounding (MIPAS). Also, we compare coincident and collocated SBUV ozone retrievals with observations made by ground-based instruments, such as microwave spectrometers, lidars, Umkehr instruments and balloon-borne ozonosondes. Finally, we compare the SBUV ozone profiles with output from the NASA GSFC GEOS-CCM model. In the stratosphere between 25 and 1 hPa the mean biases and standard deviations are within 5% for monthly mean ozone profiles. Above and below this layer the vertical resolution of the SBUV algorithm decreases and the effects of vertical smoothing should be taken into account. Though the SBUV algorithm has a coarser vertical resolution in the lower stratosphere and troposphere, it is capable of precisely estimating the integrated ozone column between the surface and 25 hPa. The time series of the tropospheric - lower stratospheric ozone column derived from SBUV agrees within 5% with the corresponding values observed by an ensemble of ozone sonde stations in North Hemisphere. Drift of the ozone time series obtained from each SBUV(/2) instrument relative to ground based and satellite measurements are evaluated and some features of individual SBUV(l2) instruments are discussed. In addition to evaluating individual instruments against independent observations, we also focus on the instrument to instrument consistency in the series. Overall, Version 8.6 ozone profiles obtained from two different SBUV(l2) instruments compare within a couple of percent during overlap periods and are consistently varying in time, with some exceptions. Some of the noted discrepancies might bssociated with ozone diurnal variations, since the difference in the local time of the observations for a pair of SBUV(l2) instruments could be several hours. Other issues include the potential short-term drift in measurements as the instrument orbit drifts, and measurements are obtained at high solar zenith angles (>85 ). Based on the results of the validation, a consistent, calibrated dataset of SBUV ozone profiles has been created based on internal calibration only.

Kramarova, Natalya↗
DOE OSTI2023

Protection Against Graph-Based False Data Injection Attacks on Power Systems

Graph signal processing (GSP) has emerged as a powerful tool for practical network applications, including power system monitoring. By representing power system voltages as smooth graph signals, recent research has focused on developing GSP-based methods for state estimation, attack detection, and topology identification. Included, efficient methods have been developed for detecting false data injection (FDI) attacks, which until now were perceived as non-smooth with respect to the graph Laplacian matrix. Consequently, these methods may not be effective against smooth FDI attacks. In this paper, we propose a graph FDI (GFDI) attack that minimizes the Laplacian-based graph total variation (TV) under practical constraints. In addition, we develop a low-complexity algorithm that solves the non-convex GDFI attack optimization problem using ell_1-norm relaxation, the projected gradient descent (PGD) algorithm, and the alternating direction method of multipliers (ADMM). We then propose a protection scheme that identifies the minimal set of measurements necessary to constrain the GFDI output to high graph TV, thereby enabling its detection by existing GSP-based detectors. Our numerical simulations on the IEEE-57 bus test case reveal the potential threat posed by well-designed GSP-based FDI attacks. Moreover, we demonstrate that integrating the proposed protection design with GSP-based detection can lead to significant hardware cost savings compared to previous designs of protection methods against FDI attacks.

Morgenstern, Gal↗
NASA NTRS1991

Coordination of human upper arm and forearm motion in rapid extensions

In many movements of the upper limb such as reaching, positioning, and displacing the objects, the hand moves smoothly along a uni-directional planar trajectory with a bell-shaped speed profile. Because of variability of the load, gravitational and velocity interaction forces between its segments, the dynamics of the arm during the motion is very complex. Motion of the upper arm and the forearm are coordinated to produce smooth movements despite such complex dynamics. This coordination constitutes an important organizational feature of arm movement. The present paper describes some experimental results related to the above and interprets their role in producing smooth motion of the hand. Trajectories of the right upper limb in vertical plane and the simultaneous muscles' EMG activities were recorded for extensions of various amplitudes under four loads. The forearm trajectories are smooth with bell-shaped speed profiles. The upper arm trajectories may have bimodal speed profiles and three segments, resembling an inverted 'Z.' Segmentation is sharper for the points near the shoulder joint, and is accentuated by load. As one moves from the central points near the shoulder joint to the peripheral points near the hand, the three segments merge and result in smooth curves with single-peaked speed profiles. The three segments have a coherent time course and can be identified rather accurately and non-ambiguously. The observed trajectories and the EMG patterns reveal an effective coordination strategy which utilizes the structure and the dynamics of the moving arm to produce smooth movement.

Nahvi, Mahmood↗
NASA NTRS2006

Automated, on-board terrain analysis for precision landings

Advances in space robotics technology hinge to a large extent upon the development and deployment of sophisticated new vision-based methods for automated in-space mission operations and scientific survey. To this end, we have developed a new concept for automated terrain analysis that is based upon a generic image enhancement platform|multi-scale retinex (MSR) and visual servo (VS) processing. This pre-conditioning with the MSR and the vs produces a "canonical" visual representation that is largely independent of lighting variations, and exposure errors. Enhanced imagery is then processed with a biologically inspired two-channel edge detection process, followed by a smoothness based criteria for image segmentation. Landing sites can be automatically determined by examining the results of the smoothness-based segmentation which shows those areas in the image that surpass a minimum degree of smoothness. Though the msr has proven to be a very strong enhancement engine, the other elements of the approach|the vs, terrain map generation, and smoothness-based segmentation|are in early stages of development. Experimental results on data from the Mars Global Surveyor show that the imagery can be processed to automatically obtain smooth landing sites. In this paper, we describe the method used to obtain these landing sites, and also examine the smoothness criteria in terms of the imager and scene characteristics. Several examples of applying this method to simulated and real imagery are shown.

Rahman, Zia-ur↗
NASA NTRS2006

Automated On-board Terrain Analysis for Precision Landings

Advances in space robotics technology hinge to a large extent upon the development and deployment of sophisticated new vision-based methods for automated in-space mission operations and scientific survey. To this end, we have developed a new concept for automated terrain analysis that is based upon a generic image enhancement platform-multi-scale Retinex (MSR) and visual servo (VS) processing. This pre-conditioning with the MSR and the VS produces a "canonical" visual representation that is largely independent of lighting variations, and exposure errors. Enhanced imagery is then processed with a biologically inspired two-channel edge detection process, followed by a smoothness based criteria for image segmentation. Landing sites can be automatically determined by examining the results of the smoothness-based segmentation which shows those areas in the image that surpass a minimum degree of smoothness. Though the MSR has proven to be a very strong enhancement engine, the other elements of the approach, the VS, terrain map generation, and smoothness-based segmentation, are in early stages of development. Experimental results on data from the Mars Global Surveyor show that the imagery can be processed to automatically obtain smooth landing sites. In this paper, we describe the method used to obtain these landing sites, and also examine the smoothness criteria in terms of the imager and scene characteristics. Several examples of applying this method to simulated and real imagery are shown.

Rahman, Zia-ur↗
NASA NTRSDate not supplied

Effect of Geometric Fidelity on the Aerodynamics of a Swept Wing with Scalloped Ice Accretion

The effect of geometric fidelity on the aerodynamics of a swept wing with a “scallop” ice shape was studied. Three geometric fidelity versions of the ice shape were studied. The High Fidelity ice shape maintained all of the highly 3D features of the scallop shape. The 3D Smooth ice shape was smoothly lofted over the High Fidelity shape in order to eliminate all of the locally 3D feature. The third ice shape was the 3D Smooth ice shape with roughness grit attached to the surface. The two 3D Smooth versions of the ice shape exhibited a flowfield characterized by a leading edge separation bubble that rolled into spanwise running vortex. The surface pressure data showed classic leading edge separation bubbles that are observed on 2D airfoils with leading edge horn ice shapes. The High Fidelity ice shape exhibited a flowfield characterized by streamwise vortices that formed downstream of the ice shape. The streamwise vortices observed downstream of the High Fidelity ice shape appeared to reduce the size of the separation suction peak observed downstream of the 3D Smooth ice shapes. This reduced the lift for the High Fidelity ice shape, when compared to the 3D Smooth ice shapes. However, these streamwise vortices may have allowed the flow to remain attached longer and slightly increased the stall angle of attack, even though the maximum lift was lower.

Aircraft icing↗
DOE OSTI2022

Improved Data Interpretation through Identification of Time Series Periodicity Changes

Analysis and interpretation of time series data is easiest when the data values occur at uniform intervals in time, but actual data may have differing data sampling frequencies, such as monthly and daily readings. Applying data analysis techniques, such as smoothing, to such a data set may not give a representative result between time segments. The ability to automatically distinguish time segments of differing data frequency would provide a means for applying data analysis independently to each segment, though a suitable blending at segment boundaries would be required. A method for detecting frequency changes was developed and applied to Gaussian and median smoothing of hydraulic head data from groundwater wells at the U.S. Department of Energy Hanford Site in southeastern Washington state. The process identifies time segments of high-frequency (daily) or low-frequency (greater than daily) data using adjusted-bandwidth Gaussian kernel density estimation and a threshold value, which are further refined to address small blocks of low-frequency data within larger blocks of high-frequency data. User-selectable levels of smoothing are then applied independently to the time segments prior to combining the segment results for a single smoothed data set. This time segment identification approach provides effective low- and high-frequency data separation, which provides a method to apply data analysis independently to each time segment.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

On the influence of over-parameterization in manifold based surrogates and deep neural operators

Constructing accurate and generalizable approximators (surrogate models) for complex physico-chemical processes exhibiting highly non-smooth dynamics is challenging. The main question is what type of surrogate models we should construct and should these models be under-parameterized or over-parameterized. In this work, we propose new developments and perform comparisons for two promising approaches: manifold-based polynomial chaos expansion (m-PCE) and the deep neural operator (DeepONet), and we examine the effect of over-parameterization on generalization. While m-PCE enables the construction of a mapping by first identifying low-dimensional embeddings of the input functions, parameters, and quantities of interest (QoIs), a neural operator learns the nonlinear mapping via the use of deep neural networks. Here, we demonstrate the performance of these methods in terms of generalization accuracy by solving the 2D time-dependent Brusselator reaction-diffusion system with uncertainty sources, modeling an autocatalytic chemical reaction between two species. We first propose an extension of the m-PCE by constructing a mapping between latent spaces formed by two separate embeddings of the input functions and the output QoIs. To further enhance the accuracy of the DeepONet, we introduce weight self-adaptivity in the loss function. We demonstrate that the performance of m-PCE and DeepONet is comparable for cases of relatively smooth input-output mappings. However, when highly non-smooth dynamics is considered, DeepONet shows higher approximation accuracy. We also find that for m-PCE, modest over-parameterization leads to better generalization, both within and outside of distribution, whereas aggressive over-parameterization leads to over-fitting. In contrast, an even highly over-parameterized DeepONet leads to better generalization for both smooth and non-smooth dynamics. Furthermore, we compare the performance of the above models with another recently proposed operator learning model, the Fourier Neural Operator, and show that its over-parameterization also leads to better generalization. Taken together, our studies show that m-PCE can provide very good accuracy at very low training cost, whereas a highly over-parameterized DeepONet can provide better accuracy and robustness to noise but at higher training cost. In both methods, the inference cost is negligible.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Integrability, normal forms, and magnetic axis coordinates

Integrable or near-integrable magnetic fields are prominent in the design of plasma confinement devices. Such a field is characterized by the existence of a singular foliation entirely consisting of invariant submanifolds. A compact regular leaf (a flux surface) of this foliation must be diffeomorphic to the two-torus. In a neighborhood of a flux surface, it is known that the magnetic field admits several exact smooth normal forms in which the field lines are straight. However, these normal forms break down near singular leaves, including elliptic and hyperbolic magnetic axes. In this work, the existence of exact smooth normal forms for integrable magnetic fields near elliptic and hyperbolic magnetic axes is established. In the elliptic case, smooth near-axis Hamada and Boozer coordinates are defined and constructed. Ultimately, these results establish previously conjectured smoothness properties for smooth solutions of the magnetohydrodynamic equilibrium equations. The key arguments are a consequence of a geometric reframing of integrability and magnetic fields: they are presymplectic systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Non-parametric Lagrangian biasing from the insights of neural nets

We present a Lagrangian model of galaxy clustering bias in which we train a neural net using the local properties of the smoothed initial density field to predict the late-time mass-weighted halo field.By fitting the mass-weighted halo field in the AbacusSummit simulations at z = 0.5, we find that including three coarsely spaced smoothing scales gives the best recovery of the halo power spectrum. Adding more smoothing scales may lead to 2–5% underestimation of the large-scale power and can cause the neural net to overfit.We find that the fitted halo-to-mass ratio can be well described by two directions in the original high-dimension feature space.Projecting the original features into these two principal components and re-training the neural net either reproduces the original training result, or outperforms it with a better match of the halo power spectrum. The elements of the principal components are unlikely to be assigned physical meanings, partly owing to the features being highly correlated between different smoothing scales.Our work illustrates a potential need to include multiple smoothing scales when studying galaxy bias, and this can be done easily with machine-learning methods that can take in high dimensional input feature space.

79 ASTRONOMY AND ASTROPHYSICS↗