Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “waveform relaxation”

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 19 records

A note on windowing for the waveform relaxation

The technique of windowing has been often used in the implementation of the waveform relaxations for solving ODE's or time dependent PDE's. Its efficiency depends upon problem stiffness and operator splitting. Using model problems, the estimates for window length and convergence rate are derived. The electiveness of windowing is then investigated for non-stiff and stiff cases respectively. lt concludes that for the former, windowing is highly recommended when a large discrepancy exists between the convergence rate on a time interval and the ones on its subintervals. For the latter, windowing does not provide any computational advantage if machine features are disregarded. The discussion is supported by experimental results.

Zhang, Hong↗

Optimal convolution SOR acceleration of waveform relaxation with application to semiconductor device simulation

In this paper we describe a novel generalized SOR (successive overrelaxation) algorithm for accelerating the convergence of the dynamic iteration method known as waveform relaxation. A new convolution SOR algorithm is presented, along with a theorem for determining the optimal convolution SOR parameter. Both analytic and experimental results are given to demonstrate that the convergence of the convolution SOR algorithm is substantially faster than that of the more obvious frequency-independent waveform SOR algorithm. Finally, to demonstrate the general applicability of this new method, it is used to solve the differential-algebraic system generated by spatial discretization of the time-dependent semiconductor device equations.

Reichelt, Mark↗

Overlapping localized exponential time differencing methods for diffusion problems

The localized exponential time differencing (ETD) based on overlapping domain decomposition has been recently introduced for extreme-scale phase field simulations of coarsening dynamics, which displays excellent parallel scalability in supercomputers. This paper serves as the first step toward building a solid mathematical foundation for this approach. We study the overlapping localized ETD schemes for a model time-dependent diffusion equation discretized in space by the standard central difference. Two methods are proposed and analyzed for solving the fully discrete localized ETD systems: the first one is based on Schwarz iteration applied at each time step and involves solving stationary problems in the subdomains at each iteration, while the second one is based on the Schwarz waveform relaxation algorithm in which time-dependent subdomain problems are solved at each iteration. The convergences of the associated iterative solutions to the corresponding fully discrete localized ETD solution and to the exact semidiscrete solution are rigorously proved. Finally, numerical experiments are also carried out to confirm theoretical results and to compare the performance of the two methods.

58 GEOSCIENCES↗

Relationship between early diastolic intraventricular pressure gradients, an index of elastic recoil, and improvements in systolic and diastolic function

BACKGROUND: Early diastolic intraventricular pressure gradients (IVPGs) have been proposed to relate to left ventricular (LV) elastic recoil and early ventricular "suction." Animal studies have demonstrated relationships between IVPGs and systolic and diastolic indices during acute ischemia. However, data on the effects of improvements in LV function in humans and the relationship to IVPGs are lacking. METHODS AND RESULTS: Eight patients undergoing CABG and/or infarct exclusion surgery had a triple-sensor high-fidelity catheter placed across the mitral valve intraoperatively for simultaneous recording of left atrial (LA), basal LV, and apical LV pressures. Hemodynamic data obtained before bypass were compared with those with similar LA pressures and heart rates obtained after bypass. From each LV waveform, the time constant of LV relaxation (tau), +dP/dt(max), and -dP/dt(max) were determined. Transesophageal echocardiography was used to determined end-diastolic (EDV) and end-systolic (ESV) volumes and ejection fractions (EF). At similar LA pressures and heart rates, IVPG increased after bypass (before bypass 1.64+/-0.79 mm Hg; after bypass 2.67+/-1.25 mm Hg; P<0.01). Significant improvements were observed in ESV, as well as in apical and basal +dP/dt(max), -dP/dt(max), and tau (each P<0.05). Overall, IVPGs correlated inversely with both ESV (IVPG=-0.027[ESV]+3.46, r=-0.64) and EDV (IVPG=-0.027[EDV]+4.30, r=-0.70). Improvements in IVPGs correlated with improvements in apical tau (Deltatau =5.93[DeltaIVPG]+4.76, r=0.91) and basal tau (Deltatau =2.41[DeltaIVPG]+5.13, r=-0.67). Relative changes in IVPGs correlated with changes in ESV (DeltaESV=-0.97[%DeltaIVPG]+23.34, r=-0.79), EDV (DeltaEDV=-1.16[%DeltaIVPG]+34.92, r=-0.84), and EF (DeltaEF=0.38[%DeltaIVPG]-8.39, r=0.85). CONCLUSIONS: Improvements in LV function also increase IVPGs. These changes in IVPGs, suggestive of increases in LV suction and elastic recoil, correlate directly with improvements in LV relaxation and ESV.

Non-NASA Center↗

Testing Gravity with Realistic Gravitational Waveforms in Pulsar Timing Arrays

We consider the effects of relaxing the assumption that gravitational waves composing the stochastic gravitational wave background (SGWB) are uncorrelated between frequencies in analyses of the data from Pulsar Timing Arrays (PTAs). While individual monochromatic plane waves are often a good approximation, a background composed of astrophysical sources cannot be monochromatic since an infinite plane wave carries no signal. We consider how relaxing this assumption allows us to extract potential information about modified dispersion relations and other fundamental physics questions, as both the group and phase velocity of waves become relevant. After developing the formalism we carry out simple Gaussian wavepacket examples and then consider more realistic waveforms, such as that from binary inspirals. When the frequency evolves only slowly across the PTA temporal baseline, the monochromatic assumption at an effective mean frequency remains a good approximation and we provide scaling relations that characterize its accuracy.

Hu, Wayne↗

Model experiment to study the effect of turbulence on risetime and waveform of N waves

Typical measured sonic boom rise times are two to five times longer than rise times calculated using molecular relaxation theory. The difference may be due to atmospheric turbulence. A model experiment was set up to study the influence of turbulence on waveform and rise time of spark-produced N waves. The N waves propagate through turbulence generated by a plane jet. The model turbulence is scaled down from atmospheric turbulence by approximately the same factor as the model N wave is scaled down from the sonic boom. Our experiments show that passage through the turbulence produces a wide variety of changes in the N waveform. Spiked and rounded N waves are observed, and average rise time is increased by a factor of about 2. A tentative observation based on data obtained so far is that rise time is always increased, never decreased, by turbulence.

Lipkens, Bart↗

High-Speed Research: 1994 Sonic Boom Workshop: Atmospheric Propagation and Acceptability Studies

The workshop proceedings include papers on atmospheric propagation and acceptability studies. Papers discussing atmospheric effects on the sonic boom waveform addressed several issues. It has long been assumed that the effects of molecular relaxation are adequately accounted for by assuming that a steady state balance between absorption and nonlinear wave steepening exists. It was shown that the unsteadiness induced by the nonuniform atmosphere precludes attaining this steady state. Further, it was shown that the random atmosphere acts as a filter, effectively filtering out high frequency components of the distorted waveform. Several different propagation models were compared, and an analysis of the sonic boom at the edge of the primary carpet established that the levels there are bounded. Finally, a discussion of the levels of the sonic boom below the sea surface was presented.

Mccurdy, David A.↗

Effect of quasiperiodic and random noise on many-body dynamical decoupling protocols

Symmetries (and their spontaneous rupturing) can be used to protect and engender novel quantum phases and lead to interesting collective phenomena. In the work of K. Agarwal and I. Martin, the authors described a general dynamical decoupling (polyfractal) protocol that can be used to engineer multiple discrete symmetries in many-body systems. This work expands on the former by studying the effect of quasiperiodic and random noise on such a dynamical scheme. We find generally that relaxation of engineered symmetry generators proceeds by (i) an initial relaxation on microscopic timescales to a prethermal plateau whose height is independent of noise, (ii) a linear relaxation regime with a noise-dependent rate, followed by (iii) a slow logarithmic relaxation regime that is only present for quasiperiodic noise. We glean the essential features of these regimes via scaling collapses and show that they can be generally explained by the spectral properties of the various noise waveforms considered. In particular, the quasiperiodic noise is characterized by highly time-dependent spectrum with a noise floor that mimics white noise, and peaks that grow sharper with time. We argue that both the noise floor and peaks contribute to the initial linear-in-time relaxation while the logarithmic regime is initiated when the peaks become sufficiently well resolved and cease to contribute to further relaxation. Here, we provide numerical evidence to justify these findings.

74 ATOMIC AND MOLECULAR PHYSICS↗

A numerical model for sonic boom propagation through an inhomogeneous, windy atmosphere

The ZEPHYRUS computer model calculates sonic boom distortion during propagation through the atmosphere. The model includes the effects of nonlinear distortion, attenuation, dispersion, and wind. Trial runs with the model indicate that, in general, stable shocks have not formed when the sonic boom reaches the ground. Also, the rise time of the lead shock may strongly depend on the overall waveform shape, and may be significantly increased by purely dispersive effects when matching occurs between the rise time and the characteristic oxygen molecular relaxation time.

Robinson, Leick D.↗

Analysis of sonic boom data to quantify distortions of shock profiles

Researchers at Penn State have been examining some sonic boom waveforms recorded during overflights by the Air Force which have become available to NASA and its contractors. The quality of the digitized data and the supporting meteorological data was such that one could test the applicability of molecular relaxation theories. In the late sixties, it had been supposed that the finite rise times in the absence of turbulence had neglected the vibrational relaxation of nitrogen molecules. Bass et al. have demonstrated that molecular relaxation definitely gives the correct order of magnitude of the observed rise times. However, the Air Force data in conjunction with the recent steady-state shock profile model theory of Kang and Pierce give the first opportunity to make a detailed quantitative assessment of the molecular relaxation hypothesis. Currently an investigation is ongoing to establish a method of quantifying the distortion of a sonic boom wave from a classic N-wave shape using the Air Force data taken at Edwards AFB in 1987. Using the premise that energy will be conserved approximately for a sonic boom wave both before and after the boom passes through the Earth's turbulent boundary layer, a classic undistorted waveform is constructed from the distorted signature received at the ground. A correlation between the mean-squared deviation of the distorted and undistorted waveforms and the distance the boom travels through the turbulence is sought.

Gionfriddo, Thomas A.↗

Influence of modulation frequency in rubidium cell frequency standards

The error signal which is used to control the frequency of the quartz crystal oscillator of a passive rubidium cell frequency standard is considered. The value of the slope of this signal, for an interrogation frequency close to the atomic transition frequency is calculated and measured for various phase (or frequency) modulation waveforms, and for several values of the modulation frequency. A theoretical analysis is made using a model which applies to a system in which the optical pumping rate, the relaxation rates and the RF field are homogeneous. Results are given for sine-wave phase modulation, square-wave frequency modulation and square-wave phase modulation. The influence of the modulation frequency on the slope of the error signal is specified. It is shown that the modulation frequency can be chosen as large as twice the non-saturated full-width at half-maximum without a drastic loss of the sensitivity to an offset of the interrogation frequency from center line, provided that the power saturation factor and the amplitude of modulation are properly adjusted.

Audoin, C.↗

Progress in modeling atmospheric propagation of sonic booms

The improved simulation of sonic boom propagation through the real atmosphere requires greater understanding of how the transient acoustic pulses popularly termed sonic booms are affected by humidity and turbulence. A realistic atmosphere is invariably somewhat turbulent, and may be characterized by an ambient fluid velocity v and sound speed c that vary from point to point. The absolute humidity will also vary from point to point, although possibly not as irregularly. What is ideally desired is a relatively simple scheme for predicting the probable spreads in key sonic boom signature parameters. Such parameters could be peak amplitudes, rise times, or gross quantities obtainable by signal processing that correlate well with annoyance or damage potential. The practical desire for the prediction scheme is that it require a relatively small amount of knowledge, possibly of a statistical nature, concerning the atmosphere along, the propagation path from the aircraft to the ground. The impact of such a scheme, if developed, implemented, and verified, would be that it would give the persons who make planning decisions a tool for assessing the magnitude of environmental problems that might result from any given overflight or sequence of overflights. The technical approach that has been followed by the author and some of his colleagues is to formulate a hierarchy of simple approximate models based on fundamental physical principles and then to test these models against existing data. For propagation of sonic booms and of other types of acoustic pulses in nonturbulent model atmospheres, there exists a basic overall theoretical model that has evolved as an outgrowth of geometrical acoustics. This theoretical model depicts the sound as propagating within ray tubes in a manner analogous to sound in a waveguide of slowly varying cross-section. Propagation along the ray tube is quasi-one-dimensional, and a wave equation for unidirectional wave propagation is used. A nonlinear term is added to this equation to account for nonlinear steepening, and the formulation has been carried through to allow for spatially varying sound speed, ambient density, and ambient wind velocities. The model intrinsically neglects diffraction, so it cannot take into account what has previously been mentioned in the literature as possibly important mechanisms for turbulence-related distortion. The model as originally developed could predict an idealized N-waveform which often agrees with data in terms of peak amplitude and overall positive phase duration. It is possible, moreover, to develop simple methods based on the physics of relaxation processes for incorporating molecular relaxation into the quasi-one-dimensional model of nonlinear propagation along ray tubes.

Pierce, Allan D.↗

CEBAF C100 Fault Classification based on Time Domain RF Signals

The CEBAF 12 GeV upgrade project, which was completed and commissioned in 2014, included the construction and installation of 80 new 7-cell superconducting cavities that were configured in 10 cryomodules. In 2018, the software and hardware in the digital low level RF systems were configured such that a fault would trigger an acquisition process which records waveform records of 17 of the RF signals for each of the 8 cavities within the cryomodule for subsequent analysis. These waveforms are especially useful in C100 cryomodules as there is a 10% mechanical coupling between adjacent cavities. When one cavity has a fault and the gradient is reduced quickly, it will mechanically deform due relaxation of the Lorentz force effects. This deformation change causes perturbations in the adjacent cavities which, in turn, causes a cascade of cavity faults that are difficult to understand without the time domain data. This contribution will describe the types of faults encountered during operation and their signatures in the time domain data, as well as how is being used to modify the setup of the machine and implement improvements to the cryomodules.

Powers, T.↗

Relaxation and turbulence effects on sonic boom signatures

The rudimentary theory of sonic booms predicts that the pressure signatures received at the ground begin with an abrupt shock, such that the overpressure is nearly abrupt. This discontinuity actually has some structure, and a finite time is required for the waveform to reach its peak value. This portion of the waveform is here termed the rise phase, and it is with this portion that this presentation is primarily concerned. Any time characterizing the duration of the rise phase is loosely called the 'rise time.' Various definitions are used in the literature for this rise time. In the present discussion the rise time can be taken as the time for the waveform to rise from 10 percent of its peak value to 90 percent of its peak value. The available data on sonic booms that appears in the open literature suggests that typical values of shock over-pressure lie in the range of 30 Pa to 200 Pa, typical values of shock duration lie in the range of 150 ms to 250 ms, and typical values of the rise time lie in the range of 1 ms to 5 ms. The understanding of the rise phase of sonic booms is important because the perceived loudness of a shock depends primarily on the structure of the rise phase. A longer rise time typically implies a less loud shock. A primary question is just what physical mechanisms are most important for the determination of the detailed structure of the rise phase.

Pierce, Allan D.↗

The making of an Alfvenic fluctuation: The resolution of a second-order analysis

Ulysses observations of the high speed polar streams show that they are largely occupied by very large amplitude Alfvenic fluctuations accompanied by many rotational discontinuities. These fluctuations have a nearly constant magnetic intensity or amplitude, and the magnetic field direction per wave cycle sweeps only through a limited arc, much as a car wiperblade would do. Barnes and Hollweg (JGR, 79, 2302, 1974) suggested that this unusual waveform could arise from an obliquely propagating and linearly polarized Alfven wave of finite amplitude. From a second-order analysis, they showed that the existence of a particular solution with a constant amplitude but could not resolve the outcome of the homogeneous solution which consisted of fast waves. They suggested that Landau damping of these fast waves may be needed to get the observed waveform. We present a 1 1/2 D hybrid simulation which is fully nonlinear and correctly describes the ion kinetics for an initially monochromatic and linearly polarized Alfven wave propagating obliquely to the background magnetic field. The wave has a large amplitude and a wavelength so long that it can be considered dispersionless for simulation times. At early times, the second harmonic in density and in magnetic field transverse to the initial wave magnetic field are generated and have more power than other harmonics. Steepening is observed with a weak fast shock emerging, but no rotational discontinuity is left behind, and instead a constant amplitude and an arc-shaped waveform is made. The compressional component which develops after the shocks have dissipated is to zeroth order better described as a pure acoustic wave than as a fast wave. This might be explained by the relaxing of the Alfven wave to a state where its ponderomotive force vanishes so that the compressional component can travel almost independently of it.

Vasquez, Bernard J.↗

Energy Budget of Brittle Fracturing in Granite Under Stress Relaxation and Creep

Creep and relaxation are the two major time-dependent fracturing processes in rocks. While a considerable amount of research has been done in understanding these two mechanisms, critical gaps remain regarding how different energy components evolve during time-dependent fracturing processes in rocks. In this study, a series of relaxation and creep experiments were conducted on prismatic Barre granite specimens in the laboratory to estimate the energy budget of brittle fracturing in granite. For the input energy, the work done by the machine (W) is calculated and for the output energy the radiated seismic energy (E R ), released in the form of acoustic emissions (AEs), is calculated as the only measurable output energy component in the conducted experiments. The low-frequency plateau (Ω 0 ) and corner frequency (f 0 ) for each AE waveform was estimated by fitting the observed AE spectra with the theoretical spectra using the Omega model. These parameters were used to estimate the seismic moments (M 0 ) based on the radiation pattern for the double couple (shear) and non-double-couple (non-shear) events. The range of f 0 and M 0 varied from 150 to 750 kHz and 10 –4 to 10 –1 N m, respectively. Moment magnitude (M w ) varied in a wider range from –9 to –6 in creep and –8.5 to –7 in relaxation. Stress drops (Δ σ ) and source radius (r) were estimated for the AEs using Brune’s model. Here, the results report on three primary observations: (1) the effects of different source mechanisms on the estimated source parameters showed that M 0 and Δ σ were higher for DC events as compared to NDC in both relaxation and creep. (2) The radiation efficiency in the case of creep is 70% higher as compared to relaxation and, (3) the stress drop estimated in relaxation and creep demonstrated a breakdown in scaling with the seismic moment.

58 GEOSCIENCES↗

Wave equations and computational models for sonic boom propagation through a turbulent atmosphere

The improved simulation of sonic boom propagation through the real atmosphere requires greater understanding of how the transient acoustic pulses popularly termed sonic booms are affected by atmospheric turbulence. A nonlinear partial differential equation that can be used to simulate the effects of smaller-scale atmospheric turbulence on sonic boom waveforms is described. The equation is first order in the time derivative and involves an extension of geometrical acoustics to include diffraction phenomena. Various terms in the equation are explained in physical terms. Such terms include those representing convection at the wave speed, diffraction, molecular relaxation, classical dissipation, and nonlinear steepening. The atmospheric turbulence enters through an effective sound speed, which varies with all three spatial coordinates, and which is the sum of the local sound speed and the component of the turbulent flow velocity projected along a central ray that connects the aircraft trajectory with the listener.

Pierce, Allan D.↗