Search NASA⌕ Search

SEARCH · Search NASA

Results for “Directional Gaussian 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.

30 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Gravitational clustering in the expanding universe - Controlled high-resolution studies in two dimensions

Results are presented from a series of gravitational clustering simulations in two dimensions. These simulations are a significant departure from previous work, since in two dimensions one can have large dynamic range in both length scale and mass using present computer technology. Controlled experiments were conducted by varying the slope of power-law initial density fluctuation spectra and varying cutoffs at large k, while holding constant the phases of individual Fourier components and the scale of nonlinearity. Filaments are found in many different simulations, even with pure power-law initial conditions. By direct comparison, filaments, called 'second-generation pancakes' are shown to arise as a consequence of mild nonlinearity on scales much larger than the correlation length and are not relics of an initial lattice or due to sparse sampling of the Fourier components. Bumps of low amplitude in the two-point correlation are found to be generic but usually only statistical fluctuations. Power spectra are much easier to relate to initial conditions, and seem to follow a simple triangular shape (on log-log plot) in the nonlinear regime. The rms density fluctuation with Gaussian smoothing is the most stable indicator of nonlinearity.

Beacom, John Francis↗
NASA NTRS1998

Pursuit Latency for Chromatic Targets

The temporal dynamics of eye movement response to a change in direction of stimulus motion has been used to compare the processing speeds of different types of stimuli (Mulligan, ARVO '97). In this study, the pursuit response to colored targets was measured to test the hypothesis that the slow response of the chromatic system (as measured using traditional temporal sensitivity measures such as contrast sensitivity) results in increased eye movement latencies. Subjects viewed a small (0.4 deg) Gaussian spot which moved downward at a speed of 6.6 deg/sec. At a variable time during the trajectory, the dot's direction of motion changed by 30 degrees, either to the right or left. Subjects were instructed to pursue the spot. Eye movements were measured using a video ophthalmoscope with an angular resolution of approximately 1 arc min and a temporal sampling rate of 60 Hz. Stimuli were modulated in chrominance for a variety of hue directions, combined with a range of small luminance increments and decrements, to insure that some of the stimuli fell in the subjects' equiluminance planes. The smooth portions of the resulting eye movement traces were fit by convolving the stimulus velocity with an exponential having variable onset latency, time constant and amplitude. Smooth eye movements with few saccades were observed for all stimuli. Pursuit responses to stimuli having a significant luminance component are well-fit by exponentials having latencies and time constants on the order of 100 msec. Increases in pursuit response latency on the order of 100-200 msec are observed in response to certain stimuli, which occur in pairs of complementary hues, corresponding to the intersection of the stimulus section with the subjects' equiluminant plane. Smooth eye movements can be made in response to purely chromatic stimuli, but are slower than responses to stimuli with a luminance component.

Mulligan, Jeffrey B.↗
NASA NTRS2013

Downscaling Satellite Precipitation with Emphasis on Extremes: A Variational 1-Norm Regularization in the Derivative Domain

The increasing availability of precipitation observations from space, e.g., from the Tropical Rainfall Measuring Mission (TRMM) and the forthcoming Global Precipitation Measuring (GPM) Mission, has fueled renewed interest in developing frameworks for downscaling and multi-sensor data fusion that can handle large data sets in computationally efficient ways while optimally reproducing desired properties of the underlying rainfall fields. Of special interest is the reproduction of extreme precipitation intensities and gradients, as these are directly relevant to hazard prediction. In this paper, we present a new formalism for downscaling satellite precipitation observations, which explicitly allows for the preservation of some key geometrical and statistical properties of spatial precipitation. These include sharp intensity gradients (due to high-intensity regions embedded within lower-intensity areas), coherent spatial structures (due to regions of slowly varying rainfall),and thicker-than-Gaussian tails of precipitation gradients and intensities. Specifically, we pose the downscaling problem as a discrete inverse problem and solve it via a regularized variational approach (variational downscaling) where the regularization term is selected to impose the desired smoothness in the solution while allowing for some steep gradients(called 1-norm or total variation regularization). We demonstrate the duality between this geometrically inspired solution and its Bayesian statistical interpretation, which is equivalent to assuming a Laplace prior distribution for the precipitation intensities in the derivative (wavelet) space. When the observation operator is not known, we discuss the effect of its misspecification and explore a previously proposed dictionary-based sparse inverse downscaling methodology to indirectly learn the observation operator from a database of coincidental high- and low-resolution observations. The proposed method and ideas are illustrated in case studies featuring the downscaling of a hurricane precipitation field.

Hurricanes↗
NASA NTRS1990

NGC 1058: Gas motions in an extended, quiescent spiral disk

Researchers investigate in detail the motion of gas in the galaxy NGC 1058 using the very large array (VLA) to map the emission in the 21-cm line. This galaxy is so nearly face-on that the contribution to the line width due to the variation of the rotational velocity across the D-array beam is small compared with the random z-motion of the gas. Researchers confirm results of earlier studies (Lewis 1987, A. and A. Suppl., 63, 515; van der Kruit and Shostak 1984, A. and A., 134, 258) of the galaxy's total neutral hydrogen (HI) and kinematics, including the fact that the rotation curve drops faster than Keplerian at the outer edge of the disk, which is interpreted as a fortuitous twist of the plane of rotation in the outer disk. However, their very high velocity resolution (2.58 km s(exp -1) after Hanning smoothing) coupled with good spatial resolution, allows researchers to measure more accurately the line width, and even to some extent its shape, throughout the disk. One of the most interesting results of this study is the remarkable constancy of the line width in the outer disk. From radius 90 to 210 seconds the Gaussian velocity dispersion (sigma sub nu) of the 21-cm line has a mean value of 5.7 km s(exp -1) (after correcting for the spectral resolution) with a dispersion of less than 0.9 km s(exp -1) (after correcting for the spectral resolution) with a dispersion of less than 0.9 km s(exp -1). Translating this directly into a kinetic temperature (Doppler temperature): T sub Dopp equals 121K (sigma sub mu exp 2/(km s(exp -1) (exp 2) gives 4000 K, with a dispersion of less than 1500 K over the outer disk. This constancy is observed even when comparing the spiral arms versus inter-arm regions, which in the radius range from 100 to 150 seconds the surface density modulates (defined as the ratio N sub peak -N sub trough/N sub peak + N sub trough) from 0.5 to 0.25 in the range 150 to 200 seconds.

Hanson, Margaret Murray↗
DOE OSTI2022

Tracer-field cross-correlations with k -nearest neighbour distributions

ABSTRACT In astronomy and cosmology significant effort is devoted to characterizing and understanding spatial cross-correlations between points – e.g galaxy positions, high energy neutrino arrival directions, X-ray and AGN sources, and continuous fields – e.g. weak lensing meiand Cosmic Microwave Background maps. Recently, we introduced the k-nearest neighbour (kNN) formalism to better characterize the clustering of discrete (point) data sets. Here, we extend it to the point – field cross-correlations analysis. It combines kNN measurements of the point data set with measurements of the field smoothed at many scales. The resulting statistics are sensitive to all orders in the joint clustering of the points and the field. We demonstrate that this approach, unlike the 2-pt cross-correlation, can measure the statistical dependence of two data sets even when there are no linear (Gaussian) correlations between them. We further demonstrate that this framework is far more effective than the two point function in detecting cross-correlations when the continuous field is contaminated by high levels of noise. For a particularly high level of noise, the cross-correlation between haloes and the underlying matter field in a cosmological simulation, between 10 h−1 Mpc and 30 h−1 Mpc, is detected at >5σ significance using the technique presented here, when the two-point cross-correlation significance is ∼1σ. Finally, we show that kNN cross-correlations of haloes and the matter field can be well modelled on quasi-linear scales using the Hybrid Effective Field Theory (HEFT) framework, with the same set of bias parameters as are used for 2-pt cross-correlations. The substantial improvement in the statistical power of detecting cross-correlations using this method makes it a promising tool for various cosmological applications.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2003

Transition in a Supersonic Boundary-Layer Due to Roughness and Acoustic Disturbances

The transition process induced by the interaction of an isolated roughness with acoustic disturbances in the free stream is numerically investigated for a boundary layer over a flat plate with a blunted leading edge at a free stream Mach number of 3.5. The roughness is assumed to be of Gaussian shape and the acoustic disturbances are introduced as boundary condition at the outer field. The governing equations are solved using the 5'h-~rder accurate weighted essentially non-oscillatory (WENO) scheme for space discretization and using third- order total-variation-diminishing (TVD) Runge- Kutta scheme for time integration. The steady field induced by the two and three-dimensional roughness is also computed. The flow field induced by two-dimensional roughness exhibits different characteristics depending on the roughness heights. At small roughness heights the flow passes smoothly over the roughness, at moderate heights the flow separates downstream of the roughness and at larger roughness heights the flow separates upstream and downstream of the roughness. Computations also show that disturbances inside the boundary layer is due to the direct interaction of the acoustic waves and isolated roughness plays a minor role in generating instability waves.

Balakumar, P.↗
DOE OSTI2025

Statistical evaluation of microscale stress conditions leading to void nucleation in the weak shock regime

Here, we investigate the heterogeneity of the stress state driven by anisotropic deformation response at the single crystal level through five statistical volume element (SVE) calculations of polycrystalline BCC tantalum. This work focuses on grain boundaries as a prominent material defect type prone to void nucleation based upon experimental observations of predominantly intergranular void nucleation in this material. The SVEs are constructed to be statistically representative of larger volumes of material and are meshed such that mean and standard deviation of grain size and orientation information is reconstructed. The computational meshes feature hexahedral (brick) elements and smooth conformal grain boundaries where significant stress concentration is known to occur, a tail effect of interest in the extreme events process of dynamic ductile damage. An existing micromechanical crystallographic plasticity model shown to capture the single crystal behavior of BCC tantalum well is used to perform the polycrystal calculations. The model includes representation of the non-Schmid effect of non-planar screw dislocation kinetics in tantalum. A three-dimensional stress state time profile predicted by damage modeling of a flyer plate impact experiment is applied as boundary conditions to each SVE. Resulting grain boundary stress state statistics are strongly non-Gaussian. Significant structural evolution is observed within the compressive hold before unloading into tension in the stress profile. Strong angular dependence of grain boundary traction magnitude with shock direction is observed. Non-Schmid effects continue to suggest their influence on propensity of microstructural defect types to nucleate voids. A general void nucleation criterion is proposed using probability theory. The general framework is specified to polycrystalline BCC tantalum in the weak shock regime to include the SVE calculations and literature molecular dynamics calculations of grain boundary void nucleation strength. Probability density functions (PDFs) are used to describe the interaction between the local stress state heterogeneity and the distributed grain boundary void nucleation strength state. A causation entropy maximization procedure removes the requirement for ad hoc selection of a PDF functional form and provides a rigorous procedure for data-based PDF determination. The resulting physically informed PDF describes the spatial appearance frequency of nucleated voids as a function of applied macroscale pressure. Lower length scale physics are thus packaged in a precise and computationally efficient way to provide computational plasticity insight to macroscale dynamic ductile damage models.

36 MATERIALS SCIENCE↗
NASA NTRS2014

Transition Marshall Space Flight Center Wind Profiler Splicing Algorithm to Launch Services Program Upper Winds Tool

NASAs LSP customers and the future SLS program rely on observations of upper-level winds for steering, loads, and trajectory calculations for the launch vehicles flight. On the day of launch, the 45th Weather Squadron (45 WS) Launch Weather Officers (LWOs) monitor the upper-level winds and provide forecasts to the launch team via the AMU-developed LSP Upper Winds tool for launches at Kennedy Space Center (KSC) and Cape Canaveral Air Force Station. This tool displays wind speed and direction profiles from rawinsondes released during launch operations, the 45th Space Wing 915-MHz Doppler Radar Wind Profilers (DRWPs) and KSC 50-MHz DRWP, and output from numerical weather prediction models.The goal of this task was to splice the wind speed and direction profiles from the 45th Space Wing (45 SW) 915-MHz Doppler radar Wind Profilers (DRWPs) and KSC 50-MHz DRWP at altitudes where the wind profiles overlap to create a smooth profile. In the first version of the LSP Upper Winds tool, the top of the 915-MHz DRWP wind profile and the bottom of the 50-MHz DRWP were not spliced, sometimes creating a discontinuity in the profile. The Marshall Space Flight Center (MSFC) Natural Environments Branch (NE) created algorithms to splice the wind profiles from the two sensors to generate an archive of vertically complete wind profiles for the SLS program. The AMU worked with MSFC NE personnel to implement these algorithms in the LSP Upper Winds tool to provide a continuous spliced wind profile.The AMU transitioned the MSFC NE algorithms to interpolate and fill data gaps in the data, implement a Gaussian weighting function to produce 50-m altitude intervals in each sensor, and splice the data together from both DRWPs. They did so by porting the MSFC NE code written with MATLAB software into Microsoft Excel Visual Basic for Applications (VBA). After testing the new algorithms in stand-alone VBA modules, the AMU replaced the existing VBA code in the LSP Upper Winds tool with the new algorithms. They then tested the code in the LSP Upper Winds tool with archived data. The tool will be delivered to the 45 WS after the 50-MHz DRWP upgrade is complete and the tool is tested with real-time data. The 50-MHz DRWP upgrade is expected to be finished in October 2014.

Space launch↗
NASA NTRS1990

On three dimensional object recognition and pose-determination: An abstraction based approach

A method of computing reliable Gaussian and mean curvature sign-map descriptors from the polynomial approximation of surfaces was demonstrated. Such descriptors which are invariant under perspective variation are suitable for hypothesis generation. A means for determining the pose of constructed geometric forms whose algebraic surface descriptors are nonlinear in terms of their orienting parameters was developed. This was done by means of linear functions which are capable of approximating nonlinear forms and determining their parameters. It was shown that biquadratic surfaces are suitable companion linear forms for cylindrical approximation and parameter estimation. The estimates provided the initial parametric approximations necessary for a nonlinear regression stage to fine tune the estimates by fitting the actual nonlinear form to the data. A hypothesis-based split-merge algorithm for extraction and pose determination of cylinders and planes which merge smoothly into other surfaces was developed. It was shown that all split-merge algorithms are hypothesis-based. A finite-state algorithm for the extraction of the boundaries of run-length regions was developed. The computation takes advantage of the run list topology and boundary direction constraints implicit in the run-length encoding.

Quek, Kok How Francis↗
DOE OSTI2022

Cross-beam energy transfer in direct-drive ICF. II. Theory and simulation of mitigation through increased laser bandwidth

The response of the cross-beam energy transfer instability (CBET) to laser bandwidth is investigated through a combination of theory and simulation. Existing linear theory is generalized to treat broadband lasers, demonstrating that CBET is most effectively suppressed when the bandwidth exceeds the ion-acoustic wave (IAW) frequency. It is shown that for such bandwidths, reverse (seed to pump) transfer becomes possible, which reduces the net energy transfer rapidly as bandwidth is increased. The CBET gain exponent in this regime scales with bandwidth (Δω) as Δω−3 for Gaussian or Lorentzian laser spectra with different scalings possible for other spectra. Comparison of our theory with linearized fluid and particle-in-cell simulations, performed with the laser-plasma simulation environment (LPSE) and vector particle in cell (VPIC) codes, respectively, finds that the model is accurate in the absence of nonlinear processes. However, linear analysis also finds that the IAW energy density scales as Δω−1, implying that nonlinear effects may be more difficult to control than the CBET scaling would suggest. Indeed, nonlinear effects are found to be present in VPIC simulations with high-intensity lasers, despite minimal apparent CBET. Nonlinear processes in the VPIC cases include particle trapping, the two-ion wave decay, and ion wave self-focusing. In some high intensity VPIC cases, these effects lead to net energy transfer from seed to pump and increases to backscatter stimulated Brillouin scattering reflectivities. Finally, for a given bandwidth, we show that improved control of nonlinear processes can be achieved via smoothing by spectral dispersion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Non-Boussinesq subgrid-scale model with dynamic tensorial coefficients

A major drawback of Boussinesq-type subgrid-scale stress models used in large-eddy simulations is the inherent assumption of alignment between large-scale strain rates and filtered subgrid-stresses. A priori analyses using direct numerical simulation (DNS) data have shown that this assumption is invalid locally as subgrid-scale stresses are poorly correlated with the large-scale strain rates [J. Bardina, J. Ferziger, and W. Reynolds, Improved subgrid-scale models for large-eddy simulation, in Proceedings of the 13th Fluid and Plasmadynamics Conference, AIAA (1980); C. Meneveau and K. Katz, Scale-invariance and turbulence models for large-eddy simulation, Ann. Rev. Fluid Mech. 32, 1 (2000)]. In the present work, a new, non-Boussinesq subgrid-scale model is presented where the model coefficients are computed dynamically. Some previous non-Boussinesq models have observed issues in providing adequate dissipation of turbulent kinetic energy [e.g., Bardina et al., Proceedings of the 13th Fluid and Plasmadynamics Conference (1980); R. A. Clark, J. Ferziger, and W.C. Reynolds. Evaluation of subgrid-scale models using an accurately simulated turbulent flow, J. Fluid Mech. 91, 1 (1979); S. Stolz and N. A. Adams, An approximate deconvolution procedure for large-eddy simulation, Phys. Fluids 11, 1699 (1999)]; however, the present model is shown to provide sufficient dissipation using dynamic coefficients. Modeled subgrid-scale Reynolds stresses satisfy the consistency requirements of the governing equations for large-eddy simulation (LES), vanish in laminar flow and at solid boundaries, and have the correct asymptotic behavior in the near-wall region of a turbulent boundary layer. The new model, referred to as the dynamic tensor-coefficient Smagorinsky model (DTCSM), has been tested in simulations of canonical flows: decaying and forced homogeneous isotropic turbulence, and wall-modeled turbulent channel flow at high Reynolds numbers. The results show favorable agreement with DNS data. It has been shown that DTCSM offers similar predictive capabilities as the dynamic Smagorinsky model for canonical flows. In order to assess the performance of DTCSM in more complex flows, wall-modeled simulations of high Reynolds number flow over a Gaussian bump (Boeing speed bump) exhibiting smooth-body flow separation are performed. Predictions of surface pressure and skin friction, compared against DNS and experimental data, show improved accuracy from DTCSM in comparison to existing static coefficient (Vreman) and dynamic Smagorinsky model. The computational cost of performing LES with this model is up to 15% higher than the dynamic Smagorinsky model.

42 ENGINEERING↗
NASA NTRS2003

HEASARC Software Archive

(1) Chandra Archive: SAO has maintained the interfaces through which HEASARC gains access to the Chandra Data Archive. At HEASARC's request, we have implemented an anonymous ftp copy of a major part of the public archive and we keep that archive up-to- date. SAO has participated in the ADEC interoperability working group, establishing guidelines or interoperability standards and prototyping such interfaces. We have provided an NVO-based prototype interface, intending to serve the HEASARC-led NVO demo project. HEASARC's Astrobrowse interface was maintained and updated. In addition, we have participated in design discussions surrounding HEASARC's Caldb project. We have attended the HEASARC Users Group meeting and presented CDA status and developments. (2) Chandra CALDB: SA0 has maintained and expanded the Chandra CALDB by including four new data file types, defining the corresponding CALDB keyword/identification structures. We have provided CALDB upgrades for the public (CIAO) and for Standard Data Processing. Approximately 40 new files have been added to the CALDB in these version releases. There have been in the past year ten of these CALDB upgrades, each with unique index configurations. In addition, with the inputs from software, archive, and calibration scientists, as well as CIAO/SDP software developers, we have defined a generalized expansion of the existing CALDB interface and indexing structure. The purpose of this is to make the CALDB more generally applicable and useful in new and future missions that will be supported archivally by HEASARC. The generalized interface will identify additional configurational keywords and permit more extensive calibration parameter and boundary condition specifications for unique file selection. HEASARC scientists and developers from SAO and GSFC have become involved in this work, which is expected to produce a new interface for general use within the current year. (3) DS9: One of the decisions that came from last year's HEADCC meeting was to make the ds9 image display program the primary vehicle for displaying line graphics (as well as images). The first step required to make this possible was to enhance the line graphics capabilities of ds9. SAO therefore spent considerable effort upgrading ds9 to use Tcl 8.4 so that the BLT line graphics package could be built and imported into ds9 from source code, rather than from a pre-built (and generally outdated) shared library. This task, which is nearly complete, allows us to extend BLT as needed for the HEAD community. Following HEADCC discussion concerning archiving and the display of archived data, we extended ds9 to support full access to many astronomical Web-based archives sites, including HEASARC, MAST, CHANDRA, SKYVIEW, ADS, NED, SIMBAD, IRAS, NVRO, SAO TDC, and FIRST. Using ds9's new internal Web access capabilities, these archives can be accessed via their Web page. FITS images, plots, spectra, and journal abstracts can be referenced, down-loaded, and displayed directly and easily in ds9. For more information, see: http://hea-www.harvard.edu/saord/ds9. Also after the HEADCC discussion concerning region filtering, we extended the Funtools sample implementation of region filtering as described in: http://hea-www.harvard.edu/saord/funtools/regions.html. In particular, we added several new composite regions for event and image filtering, including elliptical and box annuli. We also extended the panda (Pie AND Annulus) region support to include box pandas and elliptical pandas. These new composite regions are especially useful in programs that need to count photons in each separate region using only a single pass through the data. Support for these new regions was added to ds9. In the same vein, we developed new region support for filtering images using simple FITS image masks, i.e. 8-bit or 16-bit FITS images where the value of a pixel is the region id number for that pixel. Other important enhancements to DS9 this year, include supporor multiple world coordinate systems, three dimensional event file binning, image smoothing, region groups and tags, the ability to save images in a number of image formats (such as JPEG, TIFF, PNG, FITS), improvements in support for integrating external analysis tools, and support for the virtual observatory. In particular, a full-featured web browser has been implemented within D S 9 . This provides support for full access to HEASARC archive sites such as SKYVIEW and W3BROWSE, in addition to other astronomical archives sites such as MAST, CHANDRA, ADS, NED, SIMBAD, IRAS, NVRO, SA0 TDC, and FIRST. From within DS9, the archives can be searched, and FITS images, plots, spectra, and journal abstracts can be referenced, downloaded and displayed The web browser provides the basis for the built-in help facility. All DS9 documentation, including the reference manual, FAQ, Know Features, and contact information is now available to the user without the need for external display applications. New versions of DS9 maybe downloaded and installed using this facility. Two important features used in the analysis of high energy astronomical data have been implemented in the past year. The first is support for binning photon event data in three dimensions. By binning the third dimension in time or energy, users are easily able to detect variable x-ray sources and identify other physical properties of their data. Second, a number of fast smoothing algorithms have been implemented in DS9, which allow users to smooth their data in real time. Algorithms for boxcar, tophat, and gaussian smoothing are supported.

White, Nicholas↗