SEARCH · Engineering Papers

Results for “secondary waves”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 487 records · Page 27

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Time-domain analysis of planar microstrip devices using a generalized Yee-algorithm based on unstructured grids

The generalized Yee-algorithm is presented for the temporal full-wave analysis of planar microstrip devices. This algorithm has the significant advantage over the traditional Yee-algorithm in that it is based on unstructured and irregular grids. The robustness of the generalized Yee-algorithm is that structures that contain curved conductors or complex three-dimensional geometries can be more accurately, and much more conveniently modeled using standard automatic grid generation techniques. This generalized Yee-algorithm is based on the the time-marching solution of the discrete form of Maxwell's equations in their integral form. To this end, the electric and magnetic fields are discretized over a dual, irregular, and unstructured grid. The primary grid is assumed to be composed of general fitted polyhedra distributed throughout the volume. The secondary grid (or dual grid) is built up of the closed polyhedra whose edges connect the centroid's of adjacent primary cells, penetrating shared faces. Faraday's law and Ampere's law are used to update the fields normal to the primary and secondary grid faces, respectively. Subsequently, a correction scheme is introduced to project the normal fields onto the grid edges. It is shown that this scheme is stable, maintains second-order accuracy, and preserves the divergenceless nature of the flux densities. Finally, for computational efficiency the algorithm is structured as a series of sparse matrix-vector multiplications. Based on this scheme, the generalized Yee-algorithm has been implemented on vector and parallel high performance computers in a highly efficient manner.

Gedney, Stephen D.↗

NASA NTRS1986

Residual stress measurements in carbon steel

External dc magnetic field-induced changes in natural velocity of Rayleigh surface waves were measured in steel specimens under various stress conditions. The low field slopes of curves representing the fractional changes of natural velocity were proved to provide correct stress information in steels with different metallurgical properties. The slopes of curves under uniaxial compression, exceeding about one third of the yield stress, fell below zero in all the specimens when magnetized along the stress axis. The slopes under tension varied among different steels but remained positive in any circumstances. The stress effect was observed for both applied and residual stress. A physical interpretation of these results is given based on the stress-induced domain structure changes and the delta epsilon effect. Most importantly, it is found that the influence of detailed metallurgical properties cause only secondary effects on the obtained stress information.

Heyman, J. S.↗

NASA NTRS2013

Silicon Carbide Telescope Investigations for the LISA Mission

Space-based gravitational wave (GW) detectors are conceived to detect GWs in the low frequency range (mili-Hertz) by measuring the distance between free-falling proof masses in spacecraft (SC) separated by 5 Gm. The reference in the last decade has been the joint ESA-NASA mission LISA. One of the key elements of LISA is the telescope since it simultaneously gathers the light coming from the far SC (approximately or equal to 100 pW) and expands, collimates and sends the outgoing beam (2 W) to the far SC. Demanding requirements have been imposed on the telescope structure: the dimensional stability of the telescope must be approximately or equal to 1pm Hz(exp−1/2) at 3 mHz and the distance between the primary and the secondary mirrors must change by less than 2.5 micrometer over the mission lifetime to prevent defocussing. In addition the telescope structure must be light, strong and stiff. For this reason a potential on-axis telescope structure for LISA consisting of a silicon carbide (SiC) quadpod structure has been designed, constructed and tested. The coefficient of thermal expansion (CTE) in the LISA expected temperature range has been measured with a 1% accuracy which allows us to predict the shrinkage/expansion of the telescope due to temperature changes, and pico-meter dimensional stability has been measured at room temperature and at the expected operating temperature for the LISA telescope (around -6[deg]C). This work is supported by NASA Grants NNX10AJ38G and NX11AO26G,

carbide↗

NASA NTRSDate not supplied

In-situ Inspection of Reflowable-Interface Composite Joints During Cure in an Autoclave

In structural bonds, the interface between adherend and adhesive is nearly two-dimensional making it susceptible to minute quantities of contamination, which can cause weak bonds. Regulatory organizations such as the Federal Aviation Administration (FAA) often require redundant load paths in secondary-bonded, primary-structures to alleviate concerns with bonded performance. Under NASA’s Convergent Aeronautics Solutions (CAS) project, the Adhesive Free Bonding of Composites (AERoBOND) project is investigating reformulated aerospace epoxy-matrix resins to enable reflow and diffusion of the resin at the joint interface during a secondary bonding and cure process. The reflow and intermixing of the matrix resin during assembly can eliminate the material discontinuity at the interface, thereby removing the dependence of bond performance on adhesion at a nearly two-dimensional boundary. AERoBOND process development has evaluated a multitude of parameters including the materials used, the stoichiometric offset of the epoxy resins, time and temperature of cure, and thickness of each ply. Without in-situ process monitoring, the condition of a test article is unknown until mechanical testing is completed. This presentation describes the use of an in-situ ultrasonic inspection system to monitor the joining of two composite parts assembled using the AERoBOND technique. This work quantifies the elimination of the interface by measuring the wave reflection, or lack thereof, at the joint throughout the cure cycle. In addition, the results indicate when reflow and cure of the epoxy resin at the joint occurred. By using the recently developed in-situ inspection with a mobile ultrasonic transducer, localized results can be obtained across a large portion of the joint with high resolution.

Tyler B Hudson↗

NASA NTRSDate not supplied

In-situ Inspection of Reflowable-Interface Composite Joints During Cure in an Autoclave

Tyler B Hudson↗

NASA NTRS1988

Turbulence characteristics inside a turbulent spot in plane Poiseuille flow

In wall-bounded shear flows the transition to turbulence through localized disturbances goes through a pattern starting with a development of shear layers. The localized normal velocity fluctuations induce normal vorticity through the lift-up effect. These shear layers become unstable to secondary disturbances, and if the amplitudes of the disturbances are large enough, a turbulent spot develops. Investigations of the spot in boundary layers has shown that the turbulent part of the spot is very similar to a fully developed boundary layer. Wygnanski et al. (1976) showed that the mean profile at the center-symmetry plane has a logarithmic region and Johansson et al. (1987) showed that both the higher-order statistics and flow structures in the spot were the same as in the corresponding fully developed flow. In what respects the turbulence inside the Poiseuille spot is similar to fully developed turbulent channel flow is studied. The numerically simulated spot is used, where the characteristics inside the spot are compared to those of the wave packet in the wingtip area. A recent experimental investigation of the velocity field associated with the Poiseuille spot by Klingmann et al. is used for comparison.

Henningson, D. S.↗

NASA NTRS1994

Parabolized stability equations

The parabolized stability equations (PSE) are a new approach to analyze the streamwise evolution of single or interacting Fourier modes in weakly nonparallel flows such as boundary layers. The concept rests on the decomposition of every mode into a slowly varying amplitude function and a wave function with slowly varying wave number. The neglect of the small second derivatives of the slowly varying functions with respect to the streamwise variable leads to an initial boundary-value problem that can be solved by numerical marching procedures. The PSE approach is valid in convectively unstable flows. The equations for a single mode are closely related to those of the traditional eigenvalue problems for linear stability analysis. However, the PSE approach does not exploit the homogeneity of the problem and, therefore, can be utilized to analyze forced modes and the nonlinear growth and interaction of an initial disturbance field. In contrast to the traditional patching of local solutions, the PSE provide the spatial evolution of modes with proper account for their history. The PSE approach allows studies of secondary instabilities without the constraints of the Floquet analysis and reproduces the established experimental, theoretical, and computational benchmark results on transition up to the breakdown stage. The method matches or exceeds the demonstrated capabilities of current spatial Navier-Stokes solvers at a small fraction of their computational cost. Recent applications include studies on localized or distributed receptivity and prediction of transition in model environments for realistic engineering problems. This report describes the basis, intricacies, and some applications of the PSE methodology.

Herbert, Thorwald↗

NASA NTRS2002

Spin-Up Instability of a Levitated Molten Drop in MHD-Flow Transition to Turbulence

When an alternating magnetic field interacts with induced eddy currents in a conducting body, there will be a repulsive force between the body and the driving coil system generating the field. This repulsive force is the basis of electromagnetic levitation, which allows containerless processing of different materials. The eddy currents in the conducting body also generate Joule heating. Axial rotation of electromagnetically levitated objects is a common observation in levitation systems and often an undesirable side effect of such experiments on 1-g and -g. There have been recent efforts to use magnetic damping and suppress this tendency of body rotation. The first report of rotation in EML drops was attributed to a slight asymmetry of the shape and location of the levitation coils could change the axis and speed of rotation. Other theories of sample rotation include a frequency difference in the traveling electromagnetic waves and a phase difference in two different applied fields of the same frequency. All of these different mechanisms share the following characteristics: the torque is small, constant for constant field strength, and very weakly dependent on the sample's temperature and phase (solid or liquid). During experiments on the MSL-1 (First Microgravity Science Laboratory) mission of the Space Shuttle (STS-83 and STS-94, April and July 1997), a droplet of palladium-silicon alloy was electromagnetically levitated for viscosity measurements. For the non-deforming droplet, the resultant MHD flow inside the drop is inferred from motion of impurities on the surface. These observations indicate formation of a pair of co-rotating toroidal flow structures inside the spheroidal levitated drop that undergo secondary flow instabilities. As rise in the fluid temperature rises, the viscosity falls and the internal flow accelerates and becomes oscillatory; and beyond a point in the experiments, the surface impurities exhibit non-coherent chaotic motion signifying emergence of turbulence inside the drop. In this work, a background of these set of observations will be given followed by a presentation of our results on the digital particle tracking analysis that has been performed on a number of available videos. The analysis indicates that the levitated drop attains a constant rotational speed during the melting phase and formation of the co-rotating axi-symmetric laminar toroidal structures. However, the rate of axial rotation increases dramatically during the deformation of the toroidal structures anti their breakup into chaotic entities. This new data suggests an interaction between the flow inside the levitated molten drop and the driving coils in the experiments. Possible mechanisms for this interaction will be reviewed. The data will also be used to make an assessment of existing theories on droplet rotation.

Abedian, B.↗

DOE OSTI2020

Electron ionization via dark matter-electron scattering and the Migdal effect

There are currently several existing and proposed experiments designed to probe sub-GeV dark matter (DM) using electron ionization in various materials. The projected signal rates for these experiments assume that this ionization yield arises only from DM scattering directly off electron targets, ignoring secondary ionization contributions from DM scattering off nuclear targets. We investigate the validity of this assumption and show that if sub-GeV DM couples with comparable strength to both protons and electrons, as would be the case for a dark photon mediator, the ionization signal from atomic scattering via the Migdal effect scales with the atomic number Z and 3-momentum transfer q as Z 2 q 2 . The result is that the Migdal effect is always subdominant to electron scattering when the mediator is light, but that Migdal-induced ionization can dominate over electron scattering for heavy mediators and DM masses in the hundreds of MeV range. We put these two ionization processes on identical theoretical footing, address some theoretical uncertainties in the choice of atomic wave functions used to compute rates, and discuss the implications for DM scenarios where the Migdal process dominates, including for XENON10, XENON100, and the recent XENON1T results on light DM scattering.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRSDate not supplied

Analog Environments for Venus Aerosol Instrument and Mission Concept Testing

The most common hypothesis discussed for life on Venus is a habitat in the dense, persistent clouds, which are primarily sulfuric acid with a secondary water component. Although Earth lacks a complete Venus analog environment, several partial analogs exist which can inform our understanding of the requirements for cloud habitability as well as serving as potential testbeds for future Venus instruments or missions. Earth’s tropospheric clouds have larger particles, more water, and shorter lifetimes than Venus’s; however, they have a partial overlap in pressure and temperature, and are our only example of an airborne habitat.They can inform how particle residence time and microbial generation time constrain habitable atmospheric regions. On the biological side, this includes available nutrients and energy, environmental stressors such as radiation, osmotic potential and acidity, and time spent in protective but inactive forms (endospores, cysts, akinetes); on the physical side, gravity, air density and viscosity, thermal lofting, gravity waves, scavenging due to precipitation, and other airflow dynamics, as well as effective particle radius and electrical charge. They can also enable sampling strategy tests for airborne life detection, ranging from10 2 to10 8 viable cells/mL (10 -3 to 10 -9 cells/particle). Conversely, Earth's stratospheric sulfate layer lacks an active microbial presence, and is at a colder and less dense altitude but has partially analogous values in terms of H2SO4 and water concentration, particle size, and number density. These aerosols are accessible to aircraft and high-altitude balloons. Even more accessible potential analogs include sea spray, which overlaps with the smaller (haze) particle sizes and densities at Venus and contains trace organic constituents; marine fog, which at higher densities can match Venus’s larger (Mode 2) particle properties and contains a range of trace compounds including ammonia, phosphates, sulfates, and nitrates; and spray from acidic hydrothermal systems, which can cover a wide variety of sizes, densities, and chemical compositions. These analogs could improve the development of Venus cloud sampling instruments or sondes, including field tests of material compatibility, particle capture efficiency, and analysis sensitivity.

Analog↗

NASA NTRS1987

Design, fabrication and performance of small, graphite electrode, multistage depressed collectors with 200-W, CW, 8- to 18-GHz traveling-wave tubes

Small multistage depressed collectors (MDC's) which used pyrolytic graphite, ion-beam-textured pyrolytic graphite, and isotropic graphite electrodes were designed, fabricated, and evaluated in conjuntion with 200-W, continuous wave (CW), 8- to 18-GHz traveling-wave tubes (TWT's). The design, construction, and performance of the MDC's are described. The bakeout performance of the collectors, in terms of gas evolution, was indistinguishable from that of typical production tubes with copper collectors. However, preliminary results indicate that some additional radiofrequency (RF) and dc beam processing time (and/or longer or higher temperature bakeouts) may be needed beyond that of typical copper electrode collectors. This is particularly true for pyrolytic graphite electrodes and for TWT's without appendage ion pumps. Extended testing indicated good long-term stability of the textured pyrolytic graphite and isotropic graphite electrode surfaces. The isotropic graphite in particular showed considerable promise as an MDC electrode material because of its high purity, low cost, simple construction, potential for very compact overall size, and relatively low secondary electron emission yield characteristics in the as-machined state. However, considerably more testing experience is required before definitive conclusions on its suitability for electronic countermeasure systems and space TWT's can be made.

Ebihara, Ben T.↗

DOE OSTI2020

Numerical study on spray collapse process of ECN spray G injector under flash boiling conditions

This work features a numerical study of the spray collapse process that occurs under flash-boiling, gasoline-direct-injection conditions. The Engine Combustion Network's Spray G injector was modeled, and iso-octane, n-hexane, and propane were used as fuels with the goal to investigate fuel property effects. The development of individual jets and jet-to-jet interaction were analyzed to reveal the physical mechanisms of spray collapse. The results agreed qualitatively well with the experimental data of propane available from the literature and the predictions revealed the model's ability to capture the near-nozzle shock structures witnessed in the experiments. Simulations showed that the flow development of flashing sprays was strongly tied with the fuel properties. Specifically, under-expanded jets were not observed with iso-octane due to the fuel's low vapor pressure. In the case of n-hexane, the jets were under-expanded, and a low-pressure core was observed inside the jets due to the flow expansion. With propane, the very strong expansion caused the flow to become transonic and generate a unique shock structure near the nozzle exit. The shock structure revealed the formation of primary shock cells within the jets and secondary shock cells induced by jet-to-jet interaction. Based on the results, the spray collapse under flash boiling conditions was mainly attributed to two effects, i.e., the decrease in pressure at the spray center induced by the low-pressure cores inside and between the jets, and the isolation effect due to the jet-to-jet interaction which contributed to maintain the pressure differential between the spray center and the ambient.

42 ENGINEERING↗

NASA NTRS1994

The high-state/low-state transition in V794 Aquilae

We present a V-magnitude light curve of the cataclysmic binary V794 Quilae covering an 800 day time span. The system shows variations of ups to 3 mag. In particular, there are two dips of approximately 1 mag which last for about 50 days, and are followed by a rapid recovery to the original brightness, and a third dip of approximately 3 mag lasting about 100 days. These fluctuations are thought to be caused by the response of the accretion disk in the system to the cessation of mass transfer from the mass-losing secondary star. We present computations using a time-dependent accretion disk code to delineate a range of allowed behavior for th accretion disk. To model the observed light curve, we require that the model parameters alpha(sub cold) and alpha(sub hot), which characterize the degree of coupling of the viscous dissipation to the orbital shear in the low and high states of the accretion disk, be smaller than they have been inferred to be in the dwarf novae, and that the ratio alpha(sub hot)/alpha(sub cold) also be smaller. The fact that alpha seems to vary with the rate of mass transfer may provide support for the model of Vishniac and Diamond, in which the impact of the mass transfer stream onto the outer edge of the accretion disk excites internal waves which transport angular momentum outward and provide the viscous dissipation.

Honeycutt, R. Kent↗

NASA NTRS2015

Attempt of Serendipitous Science During the Mojave Volatile Prospector Field Expedition

On 23 October a partial solar eclipse occurred across parts of the southwest United States between approximately 21:09 and 23:40 (UT), with maximum obscuration, 36%, occurring at 22:29 (UT). During 21-26 October 2014 the Mojave Volatile Prospector (MVP) field expedition deployed and operated the NASA Ames Krex2 rover in the Mojave desert west of Baker, California (Fig. 1, bottom). The MVP field expedition primary goal was to characterize the surface and sub-surface soil moisture properties within desert alluvial fans, and as a secondary goal to provide mission operations simulations of the Resource Prospector (RP) mission to a Lunar pole. The partial solar eclipse provided an opportunity during MVP operations to address serendipitous science. Science instruments on Krex2 included a neutron spectrometer, a near-infrared spectrometer with associated imaging camera, and an independent camera coupled with software to characterize the surface textures of the areas encountered. All of these devices are focused upon the surface and as a result are downward looking. In addition to these science instruments, two hazard cameras are mounted on Krex2. The chief device used to monitor the partial solar eclipse was the engineering development unit of the Near-Infrared Volatile Spectrometer System (NIRVSS) near-infrared spectrometer. This device uses two separate fiber optic fed Hadamard transform spectrometers. The short-wave and long-wave spectrometers measure the 1600-2400 and 2300-3400 nm wavelength regions with resolutions of 10 and 13 nm, respectively. Data are obtained approximately every 8 seconds. The NIRVSS stares in the opposite direction as the front Krex2.

soil moisture properties↗

NASA NTRS2006

The Dryline on 22 May 2002 during IHOP_2002: Convective-Scale Measurements at the Profiling Site

A detailed analysis of the structure of a double dryline observed over the Oklahoma panhandle during the first International H2O Project (IHOP_2002) convective initiation (CI) mission on 22 May 2002 is presented. A unique and unprecedented set of high temporal and spatial resolution measurements of water vapor mixing ratio, wind, and boundary layer structure parameters were acquired using the National Aeronautics and Space Administration (NASA) scanning Raman lidar (SRL), the Goddard Lidar Observatory for Winds (GLOW), and the Holographic Airborne Rotating Lidar Instrument Experiment (HARLIE), respectively. These measurements are combined with the vertical velocity measurements derived from the National Center for Atmospheric Research (NCAR) Multiple Antenna Profiler Radar (MAPR) and radar structure function from the high-resolution University of Massachusetts frequency-modulated continuous-wave (FMCW) radar to reveal the evolution and structure of the late afternoon double-dryline boundary layer. The eastern dryline advanced and then retreated over the Homestead profiling site in the Oklahoma panhandle, providing conditions ripe for a detailed observation of the small-scale variability within the boundary layer and the dryline. In situ aircraft data, dropsonde and radiosonde data, along with NCAR S-band dual-polarization Doppler radar (S-Pol) measurements, are also used to provide the larger-scale picture of the double-dryline environment. Moisture and temperature jumps of about 3 g kg(sup -1) and 1 -2 K, respectively, were observed across the eastern radar fine line (dryline), more than the moisture jumps (1-2 g kg(sup -1)) observed across the western radar fine line (secondary dryline). Most updraft plumes observed were located on the moist side of the eastern dryline with vertical velocities exceeding 3 m s(sup -1) and variable horizontal widths of 2-5 km, although some were as wide as 7-8 km. These updrafts were up to 1.5 g kg(sup -1) moister than the surrounding environment. Although models suggested deep convection over the Oklahoma panhandle and several cloud lines were observed near the dryline, the dryline itself did not initiate any storms over the intensive observation region (IOR). Possible reasons for this lack of convection are discussed. Strong capping inversion and moisture detrainment between the lifting condensation level and the level of free convection related to an overriding drier air, together with the relatively small near-surface moisture values (less than 10 g kg(sup -1)), were detrimental to CI in this case.

Demoz, Belay↗

DOE OSTI2022

Advances in physics and applications of 3D magnetic perturbations on the J-TEXT tokamak

In the last two years, three major technical improvements have been made on J-TEXT in supporting of the expanded operation regions and diagnostic capabilities. (1) The successful commission of the 105 GHz/500 kW/1 s electron cyclotron resonance heating (ECRH) system increasing the core electron temperature from 0.9 keV up to around 1.5 keV. (2) The poloidal divertor configuration with an X-point in the high-field side has been achieved. In particular, the 400 kW electron cyclotron wave has also been successfully injected into the diverted plasma. (3) A 256-channel electron cyclotron emission imaging diagnostic system and two sets of four-channel Doppler backscattering diagnostics have been successfully developed on J-TEXT, allowing detailed measurement of the electron temperature and density fluctuations for turbulence and MHD research. The locked mode (LM), especially the 2/1 LM, is one of the biggest threats to the plasma operation. Both the thresholds of 2/1 and 3/1 LM are observed to vary non-monotonically on electron density. The electrode biasing was applied successfully to unlock the LM from either a rotating or static resonant magnetic perturbation (RMP) field. In the presence of 2/1 LM, three kinds of standing wave (SW) structures have been observed to share a similar connection to the island structure, i.e. the nodes of the SWs locate around the O- or X-points of the 2/1 island. The control and mitigation of disruption is essential to the safe operation of ITER, and it has been systematically studied by applying a RMP field, massive gas injection (MGI) and shattered pellet injection on J-TEXT. When the RMP-induced 2/1 LM is larger than a critical width, the MGI shutdown process can be significantly influenced. If the phase difference between the O-point of LM and the MGI valve is +90° (or –90°), the penetration depth and the assimilation of impurities can be enhanced (or suppressed) during the pre-thermal quench (TQ) phase and result in a faster (or slower) TQ. A secondary MGI can also suppress the runaway electron (RE) generation, if the additional high-Z impurity gas arrives at the plasma edge before TQ. When the secondary MGI has been applied after the formation of the RE current plateau, the RE current can be dissipated, and the dissipation rate increases with the injected impurity quantity but saturates with a maximum of 28 MA s –1 . The non-local transport is experimentally observed in the ion transport channel. The electron thermal diffusivity significantly increases with the ECRH power. Theoretical work shows that significant intrinsic current can be driven by electromagnetic turbulence, and the robust formation mechanism of the E × B staircase is identified from the Hasegawa–Wakatani system.

divertor↗

NASA NTRS1986

Numerical Solution of the Incompressible Navier-Stokes Equations in Three-dimensional Generalized Curvilinear Coordinates

Numerically solving the incompressible Navier-Stokes equations is known to be time consuming and expensive. Testing of the INS3D computers code, which solves these equations with the use of the pseudocompressibility method, shows this method to be an efficient way to obtain the steady state solution. The effects of the waves introduced by the pseudocompressibility method are analyzed and criteria are set and tested for the choice of the pseudocompressibility parameter which governs the artificial sound speed. The code is tested using laminar flow over a two dimensional backward-facing step, and laminar flow over a two dimensional circular cylinder. The results of the computations over the backward-facing step are in excellent agreement with experimental results. The transient solution of the flow over the cylinder impulsively started from rest is in good agreement with experimental results. However, the computed frequency of periodic shedding of vortices behind the cylinder is not in agreement with the experimental value. For a three dimensional test case, computations were conducted for a cylinder end wall junction. The saddle point separation and horseshoe vortex system appear in the computed field. The solution also shows secondary vortex filaments which wrap around the cylinder and spiral up in the wake.

Rogers, S. E.↗

NASA NTRS2021

The NANOGrav 11yr Data Set: Limits on Supermassive Black Hole Binaries in Galaxies within 500 Mpc

Supermassive black hole binaries (SMBHBs) should form frequently in galactic nuclei as a result of galaxy mergers. At subparsec separations, binaries become strong sources of low-frequency gravitational waves (GWs), targeted by Pulsar Timing Arrays. We used recent upper limits on continuous GWs from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11 yr data set to place constraints on putative SMBHBs in nearby massive galaxies. We compiled a comprehensive catalog of ∼44,000 galaxies in the local universe (up to redshift ∼0.05) and populated them with hypothetical binaries, assuming that the total mass of the binary is equal to the SMBH mass derived from global scaling relations. Assuming circular equal-mass binaries emitting at NANOGrav's most sensitive frequency of 8 nHz, we found that 216 galaxies are within NANOGrav's sensitivity volume. We ranked the potential SMBHBs based on GW detectability by calculating the total signal-to-noise ratio such binaries would induce within the NANOGrav array. We placed constraints on the chirp mass and mass ratio of the 216 hypothetical binaries. For 19 galaxies, only very unequal-mass binaries are allowed, with the mass of the secondary less than 10% that of the primary, roughly comparable to constraints on an SMBHB in the Milky Way. However, we demonstrated that the (typically large) uncertainties in the mass measurements can weaken the upper limits on the chirp mass. Additionally, we were able to exclude binaries delivered by major mergers (mass ratio of at least 1/4) for several of these galaxies. We also derived the first limit on the density of binaries delivered by major mergers purely based on GW data.

Zaven Arzoumanian↗