Search NASA⌕ Search

SEARCH · Search NASA

Results for “Gaussian bump”

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.

28 records · Page 2

Turbulence Model Validation Through Joint Experimental /Computational Studies of Separated Flow Over A Three-Dimensional Tapered Bump: Part I - Experimental Investigation

The results of a three-year experimental campaign aimed at comprehensively documenting the separated flow over a three-dimensional bump are presented with the purpose of generating a benchmark experimental database useful in validating computational fluid dynamics (CFD) flow simulations. The bump model geometry was designed to provide well-defined and repeatable smooth-body flow separation conditions that were suitable for both experiments and simulations. The bump had a Gaussian streamwise profile with a constant height equal to 8.5% of its width over the central 60% of the test section width. The remaining 40% were outboard spanwise portions that gradually taper to zero using an error function profile to minimize side-wall boundary layer effects. The model was immersed in a turbulent boundary layer that was developed on a suspended flat plate in the Notre Dame Mach 0.6 Wind Tunnel. In order to document the effect of the incoming boundary layer thickness on the flow separation, the bump model could be located at two streamwise positions. The mean velocity and turbulence intensity of the wind tunnel freestream flow field and approaching turbulent boundary layer were fully documented. The measurements of the flow separation region included surface visualization, wall shear stress using oil-film interferometry, mean and dynamic surface pressure, and planar and stereoscopic particle image velocimetry. The experiments were conducted over a range of Mach numbers from 0.05 to 0.2 corresponding to a range of Reynolds numbers based on the test section spanwise dimension (0.914m) of 1.0 × 10 6 ≤ Re L = U ∞ L/ν ≤ 4.0 × 10 6 . The bulk of the results are presented for the higher Mach number conditions of 0.1 and 0.2 with Re L = 2.0 × 10 6 and 4.0 × 10 6 , respectively. Extensive uncertainty analysis of the data was performed. The data is archived in the NASA Langley Turbulence Modeling Resource website at https://turbmodels.larc.nasa.gov/Other_exp_Data/speedbump_sep_exp.html. In addition to the experiments, a computational effort was made in parallel by the CFD group at Boeing Research & Technology highlighting the usefulness of the data set, which is outlined in the accompanying CFD report.

Flow separation↗

Simulations of a Turbulent Flow Subjected to Favorable and Adverse Pressure Gradients

This paper reports the results from a direct numerical simulation of an initially turbulent boundary layer passing over a wall-mounted “speed bump” geometry. The speed bump, represented in the form of a Gaussian distribution profile, generates a favorable pressure gradient region over the upstream half of the geometry, followed by an adverse pressure gradient over the downstream half. The boundary layer approaching the bump undergoes strong acceleration in the favorable pressure gradient region before experiencing incipient or very weak separation within the adverse pressure gradient region. These types of flows have proven to be particularly challenging to predict using lower-fidelity simulation tools based on various turbulence modeling approaches and warrant the use of the highest-fidelity simulation techniques. Simulation results are utilized to examine the key phenomena present in the flowfield, such as relaminarization/stabilization in the strong acceleration region succeeded by retransition to turbulence near the onset of adverse pressure gradient, incipient/weak separation, and development of internal layers where the sense of streamwise pressure gradient changes at the foot, apex and tail of the bump. The present direct numerical simulation is performed using a flow solver developed exclusively for graphics processing units, which is found to provide a significant speedup compared to an earlier solver optimized for central processing unit architectures.

Ali Uzun↗

Deriving Deformable Mirror Performance Requirements in Simulation with Experimental Verification

Coronagraph instruments rely on predictable and stable deformable mirror (DM) surface displacement to achieve the contrast required to detect Earth-sized exoplanets orbiting nearby Solar-type stars. Anomalous DM behavior, such as unstable or pinned actuators, can limit contrast in coronagraphs. Simulating how these undesired behaviors affect the performance of a high contrast imaging architecture is important for developing requirements on their associated hardware. Simulating a vortex coronagraph (VC) with two deformable mirrors, this study quantifies how the number of pinned actuators affect the performance of Focal Plane Wavefront Sensing and Control algorithms using both Grid Search Electric Field Conjugation (EFC) and Planned EFC, which uses Beta-Bumping. The simulation also quantifies how various types of voltage noise such as zero-mean Gaussian noise, zero-mean periodic noise, and drift can affect the contrast of a VC during an observation run. A tolerance of a change in the Mean Normalized Intensity of ${1\times10^{-11}}$ is allocated to both types of error. If Planned EFC is used, only 1 pinned actuator on both DMs can be tolerated. If only pure Grid Search EFC is used the DMs cannot have any pinned actuators. For the case of zero-mean Gaussian noise and zero-mean periodic noise, one can tolerate a noise standard deviation of no more than ${\sigma = 0.45 \text{ mV}}$. For drift, we can only tolerate ${\sigma = 0.30 \text{ mV}}$ or less. These results show that the DM electronics and the DM themselves need to be nearly defect free to avoid having more than 1 pinned actuator. The electronics need to be tested for different types of noise statistics and that both the average and standard deviation of the noise should be measured.

Dean L. Palmer↗

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↗

On the Nature of the Compact Object in SS 433. Observational Evidence of X-Ray Photon Index Saturation

We present an analysis of the X-ray spectral properties observed from black hole , candidate (BHC) binary SS 433. We have analyzed Rossi X-ray Time Explorer (RXTE) data from this source, coordinated with Green Bank Interferometer/RATAN-600. We show that SS 433 undergoes a X-ray spectral transition from the low hard state (LHS) to the intermediate state (IS). We show that the X-ray broad-band energy spectra during all spectral states are well fit by a sum of so called "Bulk Motion Comptonization (BMC) component" and by two (broad and narrow) Gaussians for the continuum and line emissions respectively. In addition to these spectral model components we also find a strong feature that we identify as a" blackbody-like (BB)" component which color temperature is in the range of 4-5 keV in 24 IS spectra during the radio outburst decay in SS 433. Our observational results on the "high temperature BB" bump leads us to suggest the presence of gravitationally redshifted annihilation line emission in this source. In fact this spectral feature has been recently reproduced in Monte Carlo simulations by Laurent and Titarchuk. We have also established the photon index saturation at about 2.3 in index vs mass accretion correlation. This index-mass accretion correlation allows us to evaluate the low limit of black hole (BH) mass of compact object in SS 433, M(sub bh) approximately > 2 solar masses, using the scaling method using BHC GX 339-4 as a reference source. Our estimate of the BH mass in SS 433 is consistent with recent BH mass measurement using the radial-velocity measurements of the binary system by Hillwig & Gies who find that M(sub x)( = (4.3 +/- 0.8) solar masses. This is the smallest BH mass found up to now among all BH sources. Moreover, the index saturation effect versus mass accretion rate revealed in SS 433, like in a number of other BH candidates, is the strong observational evidence for the presence of a BH in SS 433.

Seifina, Elena↗

Smooth Dispersion Is Physically Appropriate: Assessing and Amending the D4 Dispersion Model

The addition of dispersion corrections to density functionals is essential for accurate energy and geometry predictions. Among them, the D4 scheme is popular due to its low computational cost and high accuracy. However, due to its design, the D4 correction can occasionally lead to anomalies, such as unphysical curvature and bumps in the potential energy surface. We find these anomalies are common in the D4 model, although observable consequences are rarer than in the D3 model for reasons we explain. Nevertheless, we uncover instances of unphysical local minima and stationary points with the D4 scheme and propose two solutions that yield smoother dispersion energy as a function of nuclear position. One is trivial to implement, based on a smoother reparametrization of Gaussian weighting (D4S) to find the effective coordination number. The other replaces Gaussian weighting with soft linear interpolation (D4SL). These new approaches usually remove artificial extremum points, while maintaining accuracy.

Tkachenko, Nikolay V↗

Testing Surrogate-Based Optimization with the Fortified Branin-Hoo Extended to Four Dimensions

Some popular functions used to test global optimization algorithms have multiple local optima, all with the same value, making them all global optima. It is easy to make them more challenging by fortifying them via adding a localized bump at the location of one of the optima. In previous work the authors illustrated this for the Branin-Hoo function and the popular differential evolution algorithm, showing that the fortified Branin-Hoo required an order of magnitude more function evaluations. This paper examines the effect of fortifying the Branin-Hoo function on surrogate- based optimization, which usually proceeds by adaptive sampling. Two algorithms are considered. The EGO algorithm, which is based on a Gaussian process (GP) and an algorithm based on radial basis functions (RBF). EGO is found to be more frugal in terms of the number of required function evaluations required to identify the correct basin, but it is expensive to run on a desktop, limiting the number of times the runs could be repeated to establish sound statistics on the number of required function evaluations. The RBF algorithm was cheaper to run, providing more sound statistics on performance. A four-dimensional version of the Branin-Hoo function was introduced in order to assess the effect of dimensionality. Furthermore, it was found that the difference between the ordinary function and the fortified one was much more pronounced for the four-dimensional function compared to the two dimensional one.

97 MATHEMATICS AND COMPUTING↗

Blanco DECam Bulge Survey (BDBS)

The Blanco DECam Bulge Survey (BDBS) has imaged more than 200 square degrees of the southern Galactic bulge, providing photometry in the ugrizy filters for ~250 million unique stars. The presence of a strong foreground disk population, along with complex reddening and extreme image crowding, has made it difficult to constrain the presence of young and intermediate age stars in the bulge population. We employed an accurate cross-match of BDBS with the latest data release (EDR3) from the Gaia mission, matching more than 140 million sources with BDBS photometry and Gaia EDR3 photometry and astrometry. We relied on Gaia EDR3 astrometry, without any photometric selection, to produce clean BDBS bulge colour-magnitude diagrams (CMDs). Gaia parallaxes were used to filter out bright foreground sources, and a Gaussian mixture model fit to Galactic proper motions could identify stars kinematically consistent with bulge membership. We applied this method to 127 different bulge fields of 1 deg 2 each, with |$\ell$| ≤ 9.5 ° and –9.5 ° ≤ b ≤ –2.5 °. The astrometric cleaning procedure removes the majority of blue stars in each field, especially near the Galactic plane, where the ratio of blue to red stars is .10%, increasing to values ~20% at higher Galactic latitudes. We rule out the presence of a widespread population of stars younger than 2 Gyr. The vast majority of blue stars brighter than the turnoff belong to the foreground population, according to their measured astrometry. We introduce the distance between the observed red giant branch bump and the red clump as a simple age proxy for the dominant population in the field, and we confirm the picture of a predominantly old bulge. Further work is needed to apply the method to estimate ages to fields at higher latitudes, and to model the complex morphology of the Galactic bulge. We also produce transverse kinematic maps, recovering expected patterns related to the presence of the bar and of the X-shaped nature of the bulge.

79 ASTRONOMY AND ASTROPHYSICS↗

Computationally efficient subglacial drainage modelling using Gaussian process emulators: GlaDS-GP v1.0

Subglacial drainage models represent water flow at the ice–bed interface through coupled distributed and channelized systems to determine water pressure, discharge, and drainage system geometry. While they are used to understand processes such as the relationship between surface melt and ice flow, the number of uncertain model parameters and the computational cost of running models makes it difficult to adequately explore the high-dimensional parameter space and evaluate uncertainty in model predictions. Here, we develop Gaussian process (GP) emulators that make fast predictions with associated uncertainty of subglacial drainage model outputs. Using a truncated principal component (PC) basis representation, we construct a GP emulator for diurnally averaged subglacial water pressure. We also explore emulation of scalar variables describing drainage efficiency and configuration. We train the emulators using ensembles of up to 512 simulations varying eight parameters of the Glacier Drainage System (GlaDS) model on a synthetic domain intended to represent an ice-sheet margin. The emulators make predictions ∼ 1000 times faster than GlaDS simulations, with errors <3 % for the water pressure field and ∼ 5 %–9 % for drainage efficiency and configuration. We apply the emulators to explore the eight-dimensional parameter space by computing variance-based parameter sensitivity indices, finding that three parameters (ice flow coefficient, bed bump aspect ratio, and the subglacial cavity system conductivity) explain 90 % of the variance in modelled water pressure in response to parameter changes. The GP emulator approach described here is well suited to integrating observational data with models to make calibrated, credible predictions of subglacial drainage.

58 GEOSCIENCES↗