Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constrained wave”

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 415 records · Page 23

Geosat crossover analysis in the tropical Pacific. Part 1: Constrained sinusoidal crossover adjustment

A new method (constrained sinusoidal crossover adjustment) for removing the orbit error in satellite altimetry is tested (using crossovers accumulated in the first 91 days of the Geosat non-repeat era in the tropical Pacific) and found to have excellent qualities. Two features distinguish the new method from the conventional bias-and-tilt crossover adjustment. First, a sine wave (with wavelength equaling the circumference of the Earth) is used to represent the orbit error for each satellite revolution, instead of the bias-and-tilt (and curvature, if necessary) approach for each segment of the satellite ground track. Secondly, the indeterminacy of the adjustment process is removed by a simple constraint minimizing the amplitudes of the sine waves, rather than by fixing selected tracks. Overall the new method is more accurate, more efficient, and much less cumbersome than the old. The idea of restricting the crossover adjustment to crossovers between tracks that are less than certain days apart in order to preserve the large-scale long-term oceanic variability is also tested with inconclusive results because the orbit error was unusually nonstationary in the initial 91 days of the GEOSAT mission.

Tai, Chang-Kou↗

Geosat crossover analysis in the tropical Pacific. I - Constrained sinusoidal crossover adjustment

A new method (constrained sinusoidal crossover adjustment) for removing the orbit error in satellite altimetry is tested (using crossovers accumulated in the first 91 days of the Geosat non-repeat era in the tropical Pacific) and found to have excellent qualities. Two features distinguish the new method from the conventional bias-and-tilt crossover adjustment. First, a sine wave (with wavelength equaling the circumference of the Earth) is used to represent the orbit error for each satellite revolution, instead of the bias-and-tilt (and curvature, if necessary) approach for each segment of the satellite ground track. Secondly, the indeterminacy of the adjustment process is removed by a simple constraint minimizing the amplitudes of the sine waves, rather than by fixing selected tracks. Overall the new method is more accurate, more efficient, and much less cumbersome than the old. The idea of restricting the crossover adjustment to crossovers between tracks that are less than certain days apart in order to preserve the large-scale long-term oceanic variability is also tested with inconclusive results because the orbit error was unusually nonstationary in the initial 91 days of the GEOSAT mission.

Tai, Chang-Kou↗

Gravitational Wave Sources as Timing References for LISA Data

In the megahertz gravitational-wave band, galactic ultra-compact binaries (UCBs) are continuous sources emitting at near-constant frequency. The signals from many of these galactic binaries will be sufficiently strong to be detectable by the Laser Interferometer Space Antenna (LISA) after approximately Omicron (1 week) of observing. In addition to their astrophysical value, these UCBs can be used to monitor the data quality of the observatory. This paper demonstrates the capabilities of galactic UCBs to be used as calibration sources for LISA by demanding signal coherence between adjacent week-long data segments separated by a gap in time of a priori unknown duration. A parameter for the gap duration is added to the UCB waveform model and used in a Markov-chain Monte Carlo algorithm simultaneously fitting for the astrophysical source parameters. Results from measurements of several UCBs are combined to produce a joint posterior on the gap duration. The measurement accuracy's dependence on how much is known about the UCBs through prior observing, and seasonal variations due to the LISA orbital motion, is quantified. The duration of data gaps in a two-week segment of data can be constrained to within approximately 0.2 s using Omicron (10) UCBs after one month of observing. The timing accuracy from UCBs improves to approximately or less than o.1 s after 1 year of mission operations. These results are robust to within a factor of approximately 2 when taking into account seasonal variations.

Gravitation Wave Sources↗

Experimental shock Hugoniot and melt boundary temperatures of aluminum

Aluminum is a ubiquitous component in dynamic compression, pulsed power, and other high energy density physics studies. Its high-pressure behavior and phase diagram are extensively studied standards in shockwave physics. While theoretical calculations and multiphase equations of state have been benchmarked to velocity measurements of loading and unloading waves, pressure and density under shock, and other mechanical data, experimental temperature data under these conditions have not been reported. We conducted a series of experiments shocking and releasing aluminum 6061 and 1100 samples into lithium fluoride windows. We measured temperature at the sample–window interface under steady compression and subsequent isentropic release. These results allow us to constrain the temperature of the solid Hugoniot and the boundary between the liquid and face-centered cubic solid phases.

Hartsfield, Thomas Murray [Sandia National Laborat↗

Using a 1D Radially Symmetric Coda Envelope Model for Robust Moment Magnitude ( M w ) Estimation in Iraq’s Tectonically Diverse Zones

Robust estimation of moment magnitude (M w ) can be challenging for Iraq due to the strong lateral heterogeneity across diverse tectonic zones. We aim to improve moment magnitude estimation by investigating the reliability of using a 1D coda envelope model in diverse tectonic zones of different lateral effects and offer a way forward for reliable estimates of M w for small events that are difficult to waveform model. Iraq comprises two main tectonic zones: (1) the Outer platform, consisting of the northwestern Zagros fold-thrust belt and the Mesopotamian foredeep, and (2) the Inner Arabian platform which is overlain by the Iraqi desert. A simple 1D coda envelope model was used because coda waves have a low sensitivity to the source and path heterogeneity. Three separate coda calibrations were conducted to investigate the robustness of a single 1D calibration to fit the country: Whole-region calibration, Zagros calibration, and Mesopotamia calibration. In the whole-region calibration, we used stations from both the Zagros and Mesopotamia zones. In the two other calibration models, we used only stations that were in those particular zones. Ground-truth reference spectra derived from the coda spectral ratio method were used to constrain high-frequency site terms. There was no drastic difference when comparing the moment magnitudes calculated from the waveform modeling and the three calibration models. The results show that the 1D coda envelope model is a reliable method even for a region with diverse tectonic zones. Therefore, we recommend using the whole-region calibration model for moment magnitude estimation that provides more complete path coverage and avoids biases introduced by path correction failures. The proposed calibration is a fundamental step in updating the comprehensive earthquake catalog and probabilistic seismic hazard assessments for Iraq.

58 GEOSCIENCES↗

Evaluation of thick wall wave guide element

The SPS transmitting array requires an architecture which will provide a low weight, high efficiency and high structural rigidity. As noted above, waveguide slot arrays constitute the most desirable option. Consequently, such an array has been chosen for the SPS. Waveguide slot arrays offer high efficiency, uniform illumination, are are fairly lighweight. Bandwidths of such arrays are narrow, typically 1/2-2%. Although this does not directly impact the SPS, which transmits power at a single frequency of 2.45 GHz, the narrow bandwidth does constrain the thermal and mechanical tolerances of the antenna. The purpose of this program is to better define the electronic aspects of an SPS specific waveguide slot array. The specific aims of the program are as follows: (1) To build a full-scale half-module, 10 stick, array, the design parameters for which are to be determined by analytical considerations tempered by experimental data on a single slotted radiating stick, (2) To experimentally evaluate the completed array with respect to antenna pattern, impedance and return loss, and (3) To measure swept transmission amplitude and phase to provide a data base for design of a receiving antenna.

Nalos, E. J.↗

Quantum Annealing with Inequality Constraints: The Set Cover Problem

Abstract Quantum annealing is a promising method for solving hard optimization problems by transforming them into quadratic unconstrained binary optimization (QUBO) problems. However, when constraints are involved, particularly multiple inequality constraints, incorporating them into the objective function poses challenges. In this paper, the authors present two novel approaches for solving problems with multiple inequality constraints on a quantum annealer and apply them to the set cover problem (SCP). The first approach uses the augmented Lagrangian method to represent the constraints, while the second approach employs a higher‐order binary optimization (HUBO) formulation. The experiments show that both approaches outperform the standard approach for solving the SCP on the D‐Wave Advantage quantum annealer. The HUBO formulation performs slightly better than the augmented Lagrangian method in solving the SCP, but its scalability in terms of embeddability in the quantum chip is worse. The results demonstrate that the proposed augmented Lagrangian and HUBO methods can successfully implement a large number of inequality constraints, making them applicable to a broad range of constrained problems beyond the SCP.

Djidjev, Hristo N.↗

On the Scaling and Growth Limit of Fluvial Dunes

Abstract Bedforms emerge in a variety of shapes and sizes when a granular erodible surface is subject to a strong enough shear flow, as observed in topographic data from submarine canyons, rivers, deserts, and planetary bodies. The two salient features of bedforms are the ability to collectively transport particles, and to generate form drag thus increasing flow resistance. These two mechanisms are in competition and contribute to force bedforms of increasing size migrating more slowly. In a dedicated large‐scale open channel flow facility, hierarchies of fluvial bedforms were generated and measured in equilibrium conditions. The corresponding scale‐dependent migration velocity and mass flux contributions were quantified in the wave number and frequency domains. Experimental results are paired with a validated set of theoretical models to demonstrate that ripples or dunes reach an equilibrium state when drag partitioning ensures high enough frictional drag to sustain the bedload transport of sediment, and low enough form drag to enable the migration of the largest bedform size. This mechanism is inferred to constrain the growth of bedforms when sediment supply is not the limiting factor.

58 GEOSCIENCES↗

A finite-volume scheme for modeling compressible magnetohydrodynamic flows at low Mach numbers in stellar interiors

Fully compressible magnetohydrodynamic (MHD) simulations are a fundamental tool for investigating the role of dynamo amplification in the generation of magnetic fields in deep convective layers of stars. The flows that arise in such environments are characterized by low (sonic) Mach numbers ($\mathscr{M}$ son ≲ 10 –2 ). In these regimes, conventional MHD codes typically show excessive dissipation and tend to be inefficient as the Courant–Friedrichs–Lewy (CFL) constraint on the time step becomes too strict. In this work we present a new method for efficiently simulating MHD flows at low Mach numbers in a space-dependent gravitational potential while still retaining all effects of compressibility. The proposed scheme is implemented in the finite-volume SEVEN-LEAGUE HYDRO (SLH) code, and it makes use of a low-Mach version of the five-wave Harten–Lax–van Leer discontinuities (HLLD) solver to reduce numerical dissipation, an implicit–explicit time discretization technique based on Strang splitting to overcome the overly strict CFL constraint, and a well-balancing method that dramatically reduces the magnitude of spatial discretization errors in strongly stratified setups. The solenoidal constraint on the magnetic field is enforced by using a constrained transport method on a staggered grid. We carry out five verification tests, including the simulation of a small-scale dynamo in a star-like environment at $\mathscr{M}$ son ~ 10 –3 . We demonstrate that the proposed scheme can be used to accurately simulate compressible MHD flows in regimes of low Mach numbers and strongly stratified setups even with moderately coarse grids.

79 ASTRONOMY AND ASTROPHYSICS↗

Machine learning framework for quantum sampling of highly constrained, continuous optimization problems

In recent years, there is growing interest in using quantum computers for solving combinatorial optimization problems. In this work, we developed a generic, machine learning-based framework for mapping continuous-space inverse design problems into surrogate quadratic unconstrained binary optimization (QUBO) problems by employing a binary variational autoencoder and a factorization machine. The factorization machine is trained as a low-dimensional, binary surrogate model for the continuous design space and sampled using various QUBO samplers. Using the D-Wave Advantage hybrid sampler and simulated annealing, we demonstrate that by repeated resampling and retraining of the factorization machine, our framework finds designs that exhibit figures of merit exceeding those of its training set. We showcase the framework’s performance on two inverse design problems by optimizing (i) thermal emitter topologies for thermophotovoltaic applications and (ii) diffractive meta-gratings for highly efficient beam steering. This technique can be further scaled to leverage future developments in quantum optimization to solve advanced inverse design problems for science and engineering applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Excitation of solar p-modes

We investigate the rates at which energy is supplied to individual p-modes as a function of their frequencies nu and angular degrees l. The observationally determined rates are compared with those calculated on the hypothesis that the modes are stochastically excited by turbulent convection. The observationally determined excitation rate is assumed to be equal to the product of the mode's energy E and its (radian) line width Gamma. We obtain E from the mode's mean square surface velocity with the aid of its velocity eigenfuction. We assume that Gamma measures the mode's energy decay rate, even though quasi-elastic scattering may dominate true absorption. At fixed l, E(Gamma) arises as nu(exp 7) at low nu, reaches a peak at nu approximately equal 3.5 mHz, and then declines as nu(exp 4.4) at higher nu . At fixed nu, E(Gamma) exhibits a slow decline with increasing l. To calculate energy input rates, P(sub alpha), we rely on the mixing-length model of turbulent convection. We find entropy fluctuations to be about an order of magnitude more effective than the Reynolds stress in exciting p-modes . The calculated P(sub alpha) mimic the nu(exp 7) dependence of E(Gamma) at low nu and the nu(exp -4.4) dependence at high nu. The break of 11.4 powers in the nu-dependence of E(Gamma) across its peak is attributed to a combination of (1) the reflection of high-frequency acoustic waves just below the photosphere where the scale height drops precipitously and (2) the absence of energy-bearing eddies with short enough correlation times to excite high-frequency modes. Two parameters associated with the eddy correlation time are required to match the location and shape of the break. The appropriate values of these parameters, while not unnatural, are poorly constrained by theory. The calculated P(sub alpha) can also be made to fit the magnitude of E(Gamma) with a reasonable value for the eddy aspect ratio. Our resutls suggest a possible explanation for the decline of mode energy with increasing l at fixed nu. Entropy fluctuations couple to changes in volume associated with the oscillation mode. These decrease with decreasing n at fixed nu, becoming almost zero for the f-mode.

Goldreich, Peter↗

Constraints on the Lunar Interior Through Comparison of Observed Apollo Seismic Travel Time With Simulated Lunar Interiors

Understanding the current state of the lunar interior is key to evaluating its formation history and evolution. The analyses of seismic waves recorded by 5 years of the Apollo passive seismic network, from 1972-1977, have provided some fundamental constraints on the lunar interior structure. Here, we show that comparisons between (a) calculated seismic velocities of simulated lunar interiors and (b) measured Apollo seismic velocities can resolve new solutions. Specifically, we show that comparison of seismic travel times can be used to constrain the lower bound of possible selenotherms.

Moon↗

Multifrequency survey of the intergalactic cloud in the M96 group

The intergalactic cloud of neutral hydrogen in the M96 group are examined for signs of emission over a wide range of frequencies, from radio waves to X rays. Past or present stellar activity in the gas might have been expected to produce detectable visual infrared, CO, OH, or radio recombination-line emission. None was detected. The limits are used to study physical conditions in the intergalactic gas. In particular, B and V band limits on starlight and IRAS limits on the presence of dust strongly constrain the presence of stars or stellar by-products. However, given the uncertainties about physical conditions in the intergalactic environment, it is difficult to rule out entirely the presence of stellar-processed materials. Results of neutral hydrogen mapping from a large-scale survey of the intergalactic cloud and surrounding region are also presented. These observations confirm that the gas is confined to a large ringlike structure. The simplest interpretation remains that the intergalactic gas in Leo is primordial.

Schneider, Stephen E.↗

Dynamical Signatures of Living Systems

One of the main challenges in modeling living systems is to distinguish a random walk of physical origin (for instance, Brownian motions) from those of biological origin and that will constitute the starting point of the proposed approach. As conjectured, the biological random walk must be nonlinear. Indeed, any stochastic Markov process can be described by linear Fokker-Planck equation (or its discretized version), only that type of process has been observed in the inanimate world. However, all such processes always converge to a stable (ergodic or periodic) state, i.e., to the states of a lower complexity and high entropy. At the same time, the evolution of living systems directed toward a higher level of complexity if complexity is associated with a number of structural variations. The simplest way to mimic such a tendency is to incorporate a nonlinearity into the random walk; then the probability evolution will attain the features of diffusion equation: the formation and dissipation of shock waves initiated by small shallow wave disturbances. As a result, the evolution never "dies:" it produces new different configurations which are accompanied by an increase or decrease of entropy (the decrease takes place during formation of shock waves, the increase-during their dissipation). In other words, the evolution can be directed "against the second law of thermodynamics" by forming patterns outside of equilibrium in the probability space. Due to that, a specie is not locked up in a certain pattern of behavior: it still can perform a variety of motions, and only the statistics of these motions is constrained by this pattern. It should be emphasized that such a "twist" is based upon the concept of reflection, i.e., the existence of the self-image (adopted from psychology). The model consists of a generator of stochastic processes which represents the motor dynamics in the form of nonlinear random walks, and a simulator of the nonlinear version of the diffusion equation which represents the mental dynamics. It has been demonstrated that coupled mental-motor dynamics can simulate emerging self-organization, prey-predator games, collaboration and competition, "collective brain," etc.

Zak, M.↗

Search for Gravitational Waves from Intermediate Mass Binary Black Holes

We present the results of a weakly modeled burst search for gravitational waves from mergers of non-spinning intermediate mass black holes (IMBH) in the total mass range 100-450 solar Mass and with the component mass ratios between 1:1 and 4:1. The search was conducted on data collected by the LIGO and Virgo detectors between November of 2005 and October of 2007. No plausible signals were observed by the search which constrains the astrophysical rates of the IMBH mergers as a function of the component masses. In the most efficiently detected bin centered on 88 + 88 solar Mass , for non-spinning sources, the rate density upper limit is 0.13 per Mpc(exp 3) per Myr at the 90% confidence level.

Blackburn, L.↗

The luminous, slow-rising orphan afterglow AT2019pim as a candidate moderately relativistic outflow

ABSTRACT Classical gamma-ray bursts (GRBs) have two distinct emission episodes: prompt emission from ultrarelativistic ejecta and afterglow from shocked circumstellar material. While both components are extremely luminous in known GRBs, a variety of scenarios predict the existence of luminous afterglow emission with little or no associated high-energy prompt emission. We present AT 2019pim, the first spectroscopically confirmed afterglow with no observed high-energy emission to be identified. Serendipitously discovered during follow-up observations of a gravitational-wave trigger and located in a contemporaneous TESS sector, it is hallmarked by a fast-rising ($t \approx 2$ h), luminous ($M_{\rm UV,peak} \approx -24.4$ mag) optical transient with accompanying luminous X-ray and radio emission. No gamma-ray emission consistent with the time and location of the transient was detected by Fermi-GBM or by Konus, placing constraining limits on an accompanying GRB. We investigate several independent observational aspects of the afterglow in the context of constraints on relativistic motion and find all of them are consistent with an initial Lorentz factor of $\Gamma _0 \approx$ 10–30 for the on-axis material, significantly lower than in any well-observed GRB and consistent with the theoretically predicted ‘dirty fireball’ scenario in which the high-energy prompt emission is stifled by pair production. However, we cannot rule out a structured jet model in which only the line-of-sight material was ejected at low-$\Gamma$, off-axis from a classical high-$\Gamma$ jet core, and an on-axis GRB with below-average gamma-ray efficiency also remains a possibility. This event represents a milestone in orphan afterglow searches, demonstrating that luminous optical afterglows lacking detected GRB counterparts can be identified and spectroscopically confirmed in real time.

Perley, Daniel A. (ORCID:0000000184721996)↗

A NICER View of the Massive Pulsar PSR J0740+6620 Informed by Radio Timing and XMM-Newton Spectroscopy

We report on Bayesian estimation of the radius, mass, and hot surface regions of the massive millisecond pulsar PSR J0740+6620, conditional on pulse-profile modeling of Neutron Star Interior Composition Explorer X-ray Timing Instrument event data. We condition on informative pulsar mass, distance, and orbital inclination priors derived from the joint North American Nanohertz Observatory for Gravitational Waves and Canadian Hydrogen Intensity Mapping Experiment/Pulsar wideband radio timing measurements of Fonseca et al. We use XMM-Newton European Photon Imaging Camera spectroscopic event data to inform our X-ray likelihood function. The prior support of the pulsar radius is truncated at 16 km to ensure coverage of current dense matter models. We assume conservative priors on instrument calibration uncertainty. We constrain the equatorial radius and mass of PSR J0740+6620 to be ${12.39}_{-0.98}^{+1.30}$ km and ${2.072}_{-0.066}^{+0.067}$ M ⊙ respectively, each reported as the posterior credible interval bounded by the 16% and 84% quantiles, conditional on surface hot regions that are non-overlapping spherical caps of fully ionized hydrogen atmosphere with uniform effective temperature; a posteriori, the temperature is ${\mathrm{log}}_{10}(T\,[{\rm{K}}])={5.99}_{-0.06}^{+0.05}$ for each hot region. All software for the X-ray modeling framework is open-source and all data, model, and sample information is publicly available, including analysis notebooks and model modules in the Python language. Furthermore, our marginal likelihood function of mass and equatorial radius is proportional to the marginal joint posterior density of those parameters (within the prior support) and can thus be computed from the posterior samples.

79 ASTRONOMY AND ASTROPHYSICS↗

Flow properties of the solar wind obtained from white light data and a two-fluid model

The flow properties of the solar wind from 1 R(sub s) to 1 AU were obtained using a two fluid model constrained by density and scale height temperatures derived from white light observations, as well as knowledge of the electron temperature in coronal holes. The observations were obtained with the white light coronographs on SPARTAN 201-1 and at Mauna Loa (Hawaii), in a north polar coronal hole from 1.16 to 5.5 R(sub s) on 11 Apr. 1993. By specifying the density, temperature, Alfven wave velocity amplitude and heating function at the coronal base, it was found that the model parameters fit well the constraints of the empirical density profiles and temperatures. The optimal range of the input parameters was found to yield a higher proton temperature than electron temperature in the inner corona. The results indicate that no preferential heating of the protons at larger distances is needed to produce higher proton than electron temperatures at 1 AU, as observed in the high speed solar wind.

Habbal, Shadia Rifai↗