Search NASA⌕ Search

SEARCH · Search NASA

Results for “matched field processing”

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 235 records · Page 13

FAST CALCULATION OF THE LOMB-SCARGLE PERIODOGRAM USING GRAPHICS PROCESSING UNITS

I introduce a new code for fast calculation of the Lomb-Scargle periodogram that leverages the computing power of graphics processing units (GPUs). After establishing a background to the newly emergent field of GPU computing, I discuss the code design and narrate key parts of its source. Benchmarking calculations indicate no significant differences in accuracy compared to an equivalent CPU-based code. However, the differences in performance are pronounced; running on a low-end GPU, the code can match eight CPU cores, and on a high-end GPU it is faster by a factor approaching 30. Applications of the code include analysis of long photometric time series obtained by ongoing satellite missions and upcoming ground-based monitoring facilities, and Monte Carlo simulation of periodogram statistical properties.

79 ASTRONOMY AND ASTROPHYSICS↗

The evolution of strongly modulated, low-frequency, moderate-amplitude wave packets in a dispersive plasma

The evolution of strongly modulated wave packets in a dispersive plasma that propagate parallel to the magnetic field is studied. Modulation effects are shown to reduce significantly (about 30 percent) the rate of spreading from that due to dispersion alone. For fluidlike behavior, nonlinearity has its greatest impact on evolution when the linear sound speed and initial wave packet speeds are well matched, resulting in a strong coupling between the wave magnetic and sonic components. Ion kinetic processes reduce the impact of nonlinearity and cause the rate of spreading to approach that expected from dispersion alone as the ratio of ion and electron temperatures, Ti/Te approaches 4. For Beta equal to or greater than 1 and Ti/Te of about 1, the coupled waveforms correspond qualitatively to kinetic treatments of the derivative nonlinear Schroedinger equation.

Vasquez, Bernard J.↗

Excitation of a Parallel Plate Waveguide by an Array of Rectangular Waveguides

This work addresses the problem of excitation of a parallel plate waveguide by an array of rectangular waveguides that arises in applications such as the continuous transverse stub (CTS) antenna and dual-polarized parabolic cylindrical reflector antennas excited by a scanning line source. In order to design the junction region between the parallel plate waveguide and the linear array of rectangular waveguides, waveguide sizes have to be chosen so that the input match is adequate for the range of scan angles for both polarizations. Electromagnetic wave scattered by the junction of a parallel plate waveguide by an array of rectangular waveguides is analyzed by formulating coupled integral equations for the aperture electric field at the junction. The integral equations are solved by the method of moments. In order to make the computational process efficient and accurate, the method of weighted averaging was used to evaluate rapidly oscillating integrals encountered in the moment matrix. In addition, the real axis spectral integral is evaluated in a deformed contour for speed and accuracy. The MoM results for a large finite array have been validated by comparing its reflection coefficients with corresponding results for an infinite array generated by the commercial finite element code, HFSS. Once the aperture electric field is determined by MoM, the input reflection coefficients at each waveguide port, and coupling for each polarization over the range of useful scan angles, are easily obtained. Results for the input impedance and coupling characteristics for both the vertical and horizontal polarizations are presented over a range of scan angles. It is shown that the scan range is limited to about 35 for both polarizations and therefore the optimum waveguide is a square of size equal to about 0.62 free space wavelength.

Rengarajan, Sembiam↗

Coiled transmission line pulse generators

Methods and apparatus are provided for fabricating and constructing solid dielectric "Coiled Transmission Line" pulse generators in radial or axial coiled geometries. The pour and cure fabrication process enables a wide variety of geometries and form factors. The volume between the conductors is filled with liquid blends of monomers, polymers, oligomers, and/or cross-linkers and dielectric powders; and then cured to form high field strength and high dielectric constant solid dielectric transmission lines that intrinsically produce ideal rectangular high voltage pulses when charged and switched into matched impedance loads. Voltage levels may be increased by Marx and/or Blumlein principles incorporating spark gap or, preferentially, solid state switches (such as optically triggered thyristors) which produce reliable, high repetition rate operation. Moreover, these Marxed pulse generators can be DC charged and do not require additional pulse forming circuitry, pulse forming lines, transformers, or an a high voltage spark gap output switch. The apparatus accommodates a wide range of voltages, impedances, pulse durations, pulse repetition rates, and duty cycles. The resulting mobile or flight platform friendly cylindrical geometric configuration is much more compact, light-weight, and robust than conventional linear geometries, or pulse generators constructed from conventional components. Installing additional circuitry may accommodate optional pulse shape improvements. The Coiled Transmission Lines can also be connected in parallel to decrease the impedance, or in series to increase the pulse length.

McDonald, Kenneth Fox↗

Fourier transform methods in local gravity modeling

New algorithms were derived for computing terrain corrections, all components of the attraction of the topography at the topographic surface and the gradients of these attractions. These algoriithms utilize fast Fourier transforms, but, in contrast to methods currently in use, all divergences of the integrals are removed during the analysis. Sequential methods employing a smooth intermediate reference surface were developed to avoid the very large transforms necessary when making computations at high resolution over a wide area. A new method for the numerical solution of Molodensky's problem was developed to mitigate the convergence difficulties that occur at short wavelengths with methods based on a Taylor series expansion. A trial field on a level surface is continued analytically to the topographic surface, and compared with that predicted from gravity observations. The difference is used to compute a correction to the trial field and the process iterated. Special techniques are employed to speed convergence and prevent oscillations. Three different spectral methods for fitting a point-mass set to a gravity field given on a regular grid at constant elevation are described. Two of the methods differ in the way that the spectrum of the point-mass set, which extends to infinite wave number, is matched to that of the gravity field which is band-limited. The third method is essentially a space-domain technique in which Fourier methods are used to solve a set of simultaneous equations.

Harrison, J. C.↗

A Combined Very Large Telescope and Gemini Study of the Atmosphere of the Directly Imaged Planet, Beta Pictoris b

We analyze new/archival VLT/NaCo and Gemini/NICI high-contrast imaging of the young, self-luminous planet Beta Pictoris b in seven near-to-mid IR photometric filters, using advanced image processing methods to achieve high signal-to-noise, high precision measurements. While Beta Pic b's near-IR colors mimic those of a standard, cloudy early-to-mid L dwarf, it is overluminous in the mid-infrared compared to the field L/T dwarf sequence. Few substellar/planet-mass objects-i.e., κ And b and 1RXJ 1609B-match Beta Pic b's JHKsL photometry and its 3.1 micron and 5 micron photometry are particularly difficult to reproduce. Atmosphere models adopting cloud prescriptions and large (approx. 60 micron)dust grains fail to reproduce the Beta Pic b spectrum. However, models incorporating thick clouds similar to those found forHR8799 bcde, but also with small (a fewmicrons) modal particle sizes, yield fits consistent with the data within the uncertainties. Assuming solar abundance models, thick clouds, and small dust particles (a = 4 micron), we derive atmosphere parameters of log(g) = 3.8 +/- 0.2 and Teff = 1575-1650 K, an inferred mass of 7+4 −3 MJ, and a luminosity of log(L/L) approx. −3.80 +/- 0.02. The best-estimated planet radius, is approx. equal to 1.65 +/- 0.06 RJ, is near the upper end of allowable planet radii for hot-start models given the host star's age and likely reflects challenges constructing accurate atmospheric models. Alternatively, these radii are comfortably consistent with hot-start model predictions if Beta Pic b is younger than is approx. equal to 7 Myr, consistent with a late formation well after its host star's birth approx. 12+8 −4 Myr ago.

overluminous↗

Paramagnetic Liquid Bridge in a Gravity-Compensating Magnetic Field

Magnetic levitation was used to stabilize cylindrical columns of a paramagnetic liquid in air between two solid supports. The maximum achievable length to diameter ratio R(sub max) was approx. (3.10 +/- 0.07), very close to the Rayleigh-Plateau limit of pi. For smaller R, the stability of the column was measured as a function of the Bond number, which could be continuously varied by adjusting the strength of the magnetic field. Liquid bridges supported by two solid surfaces have been attracting scientific attention since the time of Rayleigh and Plateau. For a cylindrical bridge of length L and diameter d, it was shown theoretically that in zero gravity the maximum slenderness ratio R (identically = L/d) is pi. The stability and ultimate collapse of such bridges is of interest because of their importance in a number of industrial processes and their potential for low gravity applications. In the presence of gravity, however, the cylindrical shape of an axisymmetric bridge tends to deform, limiting its stability and decreasing the maximum achievable value of R. Theoretical studies have discussed the stability and possible shapes of axisymmetric bridges. Experiments typically are performed in either a Plateau tank, in which the bridge is surrounded by a density-matched immiscible fluid, or in a space-borne microgravity environment. It has been shown, for example, that the stability limit R can be pushed beyond pi by using flow stabilization, by acoustic radiation pressure, or by forming columns in the presence of an axial electric field. In this work, magnetic levitation was used to simulate a low gravity environment and create quasi-cylindrical liquid columns in air. Use of a magnetic field permits us to continuously vary the Bond number B identically equal to (g)(rho)d(exp 2)/4(sigma), where g is the gravitational acceleration, rho is the density of the liquid, and sigma is the surface tension of the liquid in air. The dimensionless Bond number represents the relative importance of external forces acting on the liquid column to those due to surface tension. Our central result is that in a large magnetic field gradient we could create and stabilize columns of mixtures of water and paramagnetic manganese chloride tetrahydrate (MnCl2.4H2O), achieving a length to diameter ratio very close to pi.

Mahajan, Milind P.↗

Mapping TpPa-1 covalent organic framework (COF) molecular interactions in mixed solvents via atomistic modeling and experimental study

Complex solvent environments continue to limit the widespread adoption of organic solvent nanofiltration (OSN) in many chemical industry applications. In this paper we employ a commercially available covalent organic framework (COF), TpPa-1, and force field models to molecularly map separation performance of TpPa-1 membrane in mixed solvents. To minimize time and length scale mismatch between atomistic modeling and experiments, solvent permeance was normalized with water in modeling and experimental results to enable direct comparison. Model outputs, such as organic solvent permeance and solute rejection rate, matched well with filtration results. Since the atomistic models assume that all mass transfer is via through-pore transport, the discrepancies between modeling and experimental results provide insights on the effect of linear polymer defects, adsorption and interstitial mass transfer on polycrystalline COF membrane performance. In sum, force field models can serve as digital twins of COF membranes to simulate separation processes while capturing the effects of COF structure, chemistry, and crystallinity on membrane performance in complex organic solvent environments. Finally, this approach will provide insight into future COF design and synthesis for persisting separation challenges.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Parametric amplification of angularly multiplexed waves for application to beam smoothing

The parametric amplification of angularly multiplexed waves is studied via modeling and experiments. The phase-matching properties of the amplification process with a single pump, including parasitic interactions of pairs of signal/idler, are investigated. An experimental demonstration with up to six signal waves is presented, with a quantification of the output properties, including stability relative to the input pump and seeds. A proof-of-concept study of beam smoothing induced by combining up to 12 speckle fields generated from six signal waves and six idler waves is shown.

Dorrer, C. [University of Rochester, NY (United St↗

Conceptual Design and Dynamics Testing and Modeling of a Mars Tumbleweed Rover

The NASA Langley Research Center has been developing a novel concept for a Mars planetary rover called the Mars Tumbleweed. This concept utilizes the wind to propel the rover along the Mars surface, bringing it the potential to cover vast distances not possible with current Mars rover technology. This vehicle, in its deployed configuration, must be large and lightweight to provide the ratio of drag force to rolling resistance necessary to initiate motion from rest on the Mars surface. One Tumbleweed design concept that satisfies these considerations is called the Eggbeater-Dandelion. This paper describes the basic design considerations and a proposed dynamics model of the concept for use in simulation studies. It includes a summary of rolling/bouncing dynamics tests that used videogrammetry to better understand, characterize, and validate the dynamics model assumptions, especially the effective rolling resistance in bouncing/rolling dynamic conditions. The dynamics test used cameras to capture the motion of 32 targets affixed to a test article s outer structure. Proper placement of the cameras and alignment of their respective fields of view provided adequate image resolution of multiple targets along the trajectory as the test article proceeded down the ramp. Image processing of the frames from multiple cameras was used to determine the target positions. Position data from a set of these test runs was compared with results of a three dimensional, flexible dynamics model. Model input parameters were adjusted to match the test data for runs conducted. This process presented herein provided the means to characterize the dynamics and validate the simulation of the Eggbeater-Dandelion concept. The simulation model was used to demonstrate full scale Tumbleweed motion from a stationary condition on a flat-sloped terrain using representative Mars environment parameters.

Calhoun Philip C.↗

Simulations of Bluff Body Flow Interaction for Noise Source Modeling

The current study is a continuation of our effort to characterize the details of flow interaction between two cylinders in a tandem configuration. This configuration is viewed to possess many of the pertinent flow features of the highly interactive unsteady flow field associated with the main landing gear of large civil transports. The present effort extends our previous two-dimensional, unsteady, Reynolds Averaged Navier-Stokes computations to three dimensions using a quasilaminar, zonal approach, in conjunction with a two-equation turbulence model. Two distinct separation length-to-diameter ratios of L/D = 3.7 and 1.435, representing intermediate and short separation distances between the two cylinders, are simulated. The Mach 0.166 simulations are performed at a Reynolds number of Re = 1.66 105 to match the companion experiments at NASA Langley Research Center. Extensive comparisons with the measured steady and unsteady surface pressure and off-surface particle image velocimetry data show encouraging agreement. Both prominent and some of the more subtle trends in the mean and fluctuating flow fields are correctly predicted. Both computations and the measured data reveal a more robust and energetic shedding process at L/D = 3.7 in comparison with the weaker shedding in the shorter separation case of L/D = 1.435. The vortex shedding frequency based on the computed surface pressure spectra is in reasonable agreement with the measured Strouhal frequency.

Khorrami, Medi R.↗

Spin polarisation of tt¯γγ production at NLO+PS with GoSam interfaced to MadGraph5_aMC@NLO

Here, we present an interface between the multipurpose Monte Carlo tool MadGraph5_aMC@NLO and the automated amplitude generator GoSam. As a first application of this novel framework, we compute the NLO corrections to pp→ tt¯H and pp→ tt¯γγ matched to a parton shower. In the phenomenological analyses of these processes, we focus our attention on observables which are sensitive to the polarisation of the top quarks.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The velocity fields of elliptical galaxies: Steps toward a solution of the intrinsic shape problem

One of the few surviving signatures, at low redshift, of the process of galaxy formation should be the distribution of shapes of elliptical galaxies. Yet the problem of inferring this distribution from the observed ellipticals is still unsolved, because insufficient use has been made of kinematic information. The levels of 'sophistication' of the theoretical models and of the observations have up to now been poorly matched. The kinematic data available for most ellipticals consists of only major and minor axis spectra; and Franx et al. (1991) find, using simple geometric models, that the addition of only one kinematic parameter (the ratio of minor axis to major axis rotation velocity) to the photometry is just not enough to finely constrain the intrinsic shape distribution. On the other hand, the more elaborate self-consistent models (e.g., Levison and Richstone 1987, Statler 1987) have made only infrequent and model-dependent predictions of complicated velocity patterns, mostly at small radii, and have not discussed how they change with shape.

Statler, Thomas S.↗

Computationally efficient and error aware surrogate construction for numerical solutions of subsurface flow through porous media

Limiting the injection rate to restrict the pressure below a threshold at a critical location can be an important goal of simulations that model the subsurface pressure between injection and extraction wells. The pressure is approximated by the solution of Darcy’s partial differential equation for a given permeability field. The subsurface permeability is modeled as a random field since it is known only up to statistical properties. This induces uncertainty in the computed pressure. Solving the partial differential equation for an ensemble of random permeability simulations enables estimating a probability distribution for the pressure at the critical location. These simulations are computationally expensive, and practitioners often need rapid online guidance for real-time pressure management. An ensemble of numerical partial differential equation solutions is used to construct a Gaussian process regression model that can quickly predict the pressure at the critical location as a function of the extraction rate and permeability realization. The Gaussian process surrogate analyzes the ensemble of numerical pressure solutions at the critical location as noisy observations of the true pressure solution, enabling robust inference using the conditional Gaussian process distribution. Our first novel contribution is to identify a sampling methodology for the random environment and matching kernel technology for which fitting the Gaussian process regression model scales as O ( n log n ) instead of the typical O ( n 3 ) rate in the number of samples n used to fit the surrogate. The surrogate model allows almost instantaneous predictions for the pressure at the critical location as a function of the extraction rate and permeability realization. Our second contribution is a novel algorithm to calibrate the uncertainty in the surrogate model to the discrepancy between the true pressure solution of Darcy’s equation and the numerical solution. Finally, although our method is derived for building a surrogate for the solution of Darcy’s equation with a random permeability field, the framework broadly applies to solutions of other partial differential equations with random coefficients.

54 ENVIRONMENTAL SCIENCES↗

Petascale turbulence simulation using a highly parallel fast multipole method on GPUs

This paper reports large-scale direct numerical simulations of homogeneous-isotropic fluid turbulence, achieving sustained performance of 1.08 petaflop/s on gpu hardware using single precision. The simulations use a vortex particle method to solve the Navier–Stokes equations, with a highly parallel fast multipole method (fmm) as numerical engine, and match the current record in mesh size for this application, a cube of computational points solved with a spectral method. The standard numerical approach used in this field is the pseudo-spectral method, relying on the fft algorithm as the numerical engine. The particle-based simulations presented in this paper quantitatively match the kinetic energy spectrum obtained with a pseudo-spectral method, using a trusted code. In terms of parallel performance, weak scaling results show the fmm-based vortex method achieving 74% parallel efficiency on 4096 processes (one gpu per mpi process, 3 gpus per node of the tsubame-2.0 system). The fft-based spectral method is able to achieve just 14% parallel efficiency on the same number of mpi processes (using only cpu cores), due to the all-to-all communication pattern of the fft algorithm. The calculation time for one time step was 108 s for the vortex method and 154 s for the spectral method, under these conditions. Computing with 69 billion particles, this work exceeds by an order of magnitude the largest vortex-method calculations to date.

Computer Science↗

Saturn's Weather‐Driven Aurorae Modulate Oscillations in the Magnetic Field and Radio Emissions

Abstract The Cassini spacecraft revealed that Saturn's magnetic field displayed oscillations at a period originally thought to match the planetary rotation rate but later found not to. One of many proposed theories predicts that a polar twin‐cell neutral weather system drives this variation, producing observable differences in flows within Saturn's ionosphere. Here, using spectral observations of auroral emission lines taken by the Keck Observatory's Near Infrared Echelle Spectrograph (Keck‐NIRSPEC) in 2017, we derive ion line‐of‐sight velocity maps after grouping spectra into rotational quadrants matching phases of the planetary magnetic field. We measure 0.5 km s −1 wind systems in the ionosphere consistent with predicted neutral twin‐vortex flow patterns. These findings demonstrate that neutral winds in Saturn's polar regions cause the rotational period, as determined via the magnetic field, to exhibit differences and time variabilities relative to the planet's true period of rotation in a process never before seen within planetary atmospheres.

Chowdhury, M. N.↗

Technical Note: Development and performance of a software tool for quality assurance of online replanning with a conventional Linac or MR-Linac

Purpose: To develop an integrated quality assurance (QA) software tool for online replanning capable of efficiently and automatically checking radiation treatment (RT) planning parameters and gross plan quality, verifying treatment plan data transfer from treatment planning system (TPS) to record and verify (R&V) system, performing a secondary monitor unit (MU) calculation with or without a presence of a magnetic field from MR-Linac, and validating the delivery record consistency with the plan. Methods: The software tool, named ArtQA, was developed to obtain and compare plan and treatment parameters from both the TPS and the R&V system database. The TPS data are accessed via direct file reading and the R&V data are retrieved via open database connectivity and structured query language. Plan quality is evaluated with both the logical consistency of planning parameters and the achieved dose–volume histograms. Beams in between the TPS and R&V system are matched based on geometry configurations. To consider the effect of a 1.5 T transverse magnetic field from MR-Linac in the secondary MU calculation, a method based on modified Clarkson integration algorithm was developed and tested for a series of clinical situations. Results: ArtQA has been used in their clinic and can quickly detect inconsistencies and deviations in the entire RT planning process. With the use of the ArtQA tool, the efficiency for plan check including plan quality, data transfer, and delivery check can be improved by at least 60%. The newly developed independent MU calculation tool for MR-Linac reduces the difference between the plan and calculated MUs by 10%. Conclusions: The software tool ArtQA can be used to perform a comprehensive QA check from planning to delivery with conventional Linac or MR-Linac and is an essential tool for online replanning where the QA check needs to be performed rapidly.

60 APPLIED LIFE SCIENCES↗

Numerical experiments in boundary-layer stability

Numerical solution of the three-dimensional incompressible Navier-Stokes equations is used to study the instability of a flat-plate boundary layer in a manner analogous to the vibrating-ribbon experiments. Flow-field structures are observed which are very similar to those found in the vibrating-ribbon experiment to which computational initial conditions have been matched. Streamwise periodicity is assumed in the simulation so that the evolution occurs in time, but the events which constitute the instability are so similar to the spatially occurring ones of the laboratory that it seems clear the physical processes involved are the same. A spectral and finite difference numerical algorithm is employed in the simulation.

Wray, A.↗