Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “wave spectrum”

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 469 records · Page 26

HF-induced parametric decay in presence of field-aligned irregularities

A nonlinear theory for the saturation of the HF-induced parametric decay instability is formulated in terms of the normal modes of a field-aligned plasma wave duct. Coupling coefficients between modes are derived, and a saturation spectrum is computed numerically for a typical irregularity. The computed spectrum shows that unstable waves may propagate at larger angles to the magnetic field than would be the case in a regular medium. The relevance of this work to the experiment of Muldrew and Showen (1977), concerning the height of the HF-induced plasma line at Arecibo, is discussed. A possible explanation is also found for the 'decay line' feature in comparison with the 'broad bump' in the observed backscatter spectrum.

Rypdal, K.↗

Theoretical study of turbulent channel flow - Bulk properties, pressure fluctuations, and propagation of electromagnetic waves

In this paper, we apply two theoretical turbulence models, DIA and the recent GISS model, to study properties of a turbulent channel flow. Both models provide a turbulent kinetic energy spectral function E(k) as the solution of a non-linear equation; the two models employ the same source function but different closures. The source function is characterized by a rate n sub s (k) which is derived from the complex eigenvalues of the Orr-Sommerfeld (OS) equation in which the basic flow is taken to be of a Poiseuille type. The O-S equation is solved for a variety of Reynolds numbers corresponding to available experimental data. A physical argument is presented whereby the central line velocity characterizing the basic flow, U0 sup L, is not to be identified with the U0 appearing in the experimental Reynolds number. The theoretical results are compared with two types of experimental data: (1) turbulence bulk properties, and (2) properties that depend strongly on the structure of the turbulence spectrum at low wave numbers. The only existing analytical expression for Pi (k) cannot be used in the present case because it applies to the case of a flat plate, not a finite channel.

Canuto, V. M.↗

Tropical Waves and the Quasi-Biennial Oscillation in a 7-km Global Climate Simulation

This study investigates tropical waves and their role in driving a quasi-biennial oscillation (QBO)-like signal in stratospheric winds in a global 7-km-horizontal-resolution atmospheric general circulation model. The Nature Run (NR) is a 2-year global mesoscale simulation of the Goddard Earth Observing System Model, version 5 (GEOS-5). In the tropics, there is evidence that the NR supports a broad range of convectively generated waves. The NR precipitation spectrum resembles the observed spectrum in many aspects, including the preference for westward-propagating waves. However, even with very high horizontal resolution and a healthy population of resolved waves, the zonal force provided by the resolved waves is still too low in the QBO region and parameterized gravity wave drag is the main driver of the NR QBO-like oscillation (NRQBO). The authors suggest that causes include coarse vertical resolution and excessive dissipation. Nevertheless, the very-high-resolution NR provides an opportunity to analyze the resolved wave forcing of the NR-QBO. In agreement with previous studies, large-scale Kelvin and small-scale waves contribute to the NRQBO driving in eastward shear zones and small-scale waves dominate the NR-QBO driving in westward shear zones. Waves with zonal wavelength,1000 km account for up to half of the small-scale (,3300 km) resolved wave forcing in eastward shear zones and up to 70% of the small-scale resolved wave forcing in westward shear zones of the NR-QBO.

NR-QBO↗

Metaplectic geometrical optics for ray-based modeling of caustics: Theory and algorithms

The optimization of radio frequency-wave (RF) systems for fusion experiments is often performed using ray-tracing codes, which rely on the geometrical-optics (GO) approximation. However, GO fails at caustics such as cutoffs and focal points, erroneously predicting the wave intensity to be infinite. This is a critical shortcoming of GO, since the caustic wave intensity is often the quantity of interest, e.g., RF heating. Full-wave modeling can be used instead, but the computational cost limits the speed at which such optimizations can be performed. Here, we have developed a less expensive alternative called metaplectic geometrical optics (MGO). Instead of evolving waves in the usual x (coordinate) or k (spectral) representation, MGO uses a mixed X$\equiv$Ax+Bk representation. By continuously adjusting the matrix coefficients A and B along the rays, one can ensure that GO remains valid in the X coordinates without caustic singularities. The caustic-free result is then mapped back onto the original x space using metaplectic transforms. Here, we overview the MGO theory and review algorithms that will aid the development of an MGO-based ray-tracing code. We show how using orthosymplectic transformations leads to considerable simplifications compared to previously published MGO formulas. We also prove explicitly that MGO exactly reproduces standard GO when evaluated far from caustics (an important property that until now has only been inferred from numerical simulations), and we relate MGO to other semiclassical caustic-removal schemes published in the literature. Finally this discussion is then augmented by an explicit comparison of the computed spectrum for a wave bounded between two cutoffs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Radio wave scattering observations of the solar corona First-order measurements of expansion velocity and turbulence spectrum using Viking and Mariner 10 spacecraft

Radio wave scattering data were collected at 3.6 and 13 cm wavelengths by means of the radio link between the Viking orbiters and the earth during the Nov. 25, 1976 solar conjunction of Mars, which occurred near the beginning of solar cycle 21; Mariner 10 solar activity observations during 1974 are also used. It is found that the temporal frequency variance spectrum of amplitude fluctuations is useful for characterizing the bulk motion of the plasma, and the spectral index of electron density turbulence is obtained. The measurements of solar wind velocity and spectral index cover 78 days for Viking and 49 days for Mariner 10 and show the combined effects of changing heliocentric distance, solar latitude, and solar longitude as well as solar activity. It is concluded that the observational velocity profile differs significantly from the theoretical profiles in two ways: (1) the theoretical profile does not show the abrupt change in velocity at about 15 solar radii, and (2) the observational profile shows acceleration at larger radial distances than the model profiles. The observational profiles indicate velocities of less than about 150 km/sec out to 15 solar radii.

Tyler, G. L.↗

A New Paradigm for Supergranulation Derived from Large-Distance Time-Distance Helioseismology: Pancakes

With large separations (10-24 deg heliocentric), it has proven possible to cleanly separate the horizontal and vertical components of supergranular flow with time-distance helioseismology. These measurements require very broad filters in the k-$\omega$ power spectrum as apparently supergranulation scatters waves over a large area of the power spectrum. By picking locations of supergranulation as peaks in the horizontal divergence signal derived from f-mode waves, it is possible to simultaneously obtain average properties of supergranules and a high signal/noise ratio by averaging over many cells. By comparing ray-theory forward modeling with HMI measurements, an average supergranule model with a peak upflow of 240 m/s at cell center at a depth of 2.3 Mm and a peak horizontal outflow of 700 m/s at a depth of 1.6 Mm. This upflow is a factor of 20 larger than the measured photospheric upflow. These results may not be consistent with earlier measurements using much shorter separations (<5 deg heliocentric). With a 30 Mm horizontal extent and a few Mm in depth, the cells might be characterized as thick pancakes.

Duvall, Thomas L.↗

Three-wave coupling observed between a shear Alfvén wave and a kink-unstable magnetic flux rope

Results from a laboratory experiment are presented in which, for the first time, a shear Alfvén wave is launched using an antenna in a current-carrying plasma column that is tailored to be either stable or unstable to the kink oscillation. As the plasma is driven kink unstable, the frequency power spectrum of the Alfvén wave evolves from a single peak to a peak with multiple sidebands separated by integer multiples of the kink frequency. The main sidebands (one on either side of the launched wave peak in the power spectrum) are analyzed using azimuthal wavenumber matching, perpendicular and parallel wavenumber decomposition, and bispectral time series analysis. The dispersion relation and three-wave matching conditions are satisfied, given each sideband is a propagating Alfvén wave that results from the interaction of the pump Alfvén wave and the co-propagating component of a half-wavelength, standing kink mode. The interaction is shown to generate smaller perpendicular wavelength Alfvén waves that drive energy transport to scales that will approach the dissipation scale of k⊥ρs=1, with k⊥ being the perpendicular wavenumber and ρs being the ion gyroradius at the electron temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On Dirac equations for linear magnetoacoustic waves propagating in an isothermal atmosphere

A new analytical approach to study linear magnetoacoustic waves propagating in an isothermal, stratified, and uniformly magnetized atmosphere is presented. The approach is based on Dirac equations, and the theory of Sturm-Liouville operators is used to investigate spectral properties of the obtained Dirac Hamiltonians. Two cases are considered: (1) the background magnetic field is vertical, and the waves are separated into purely magnetic (transverse) and purely acoustic (longitudinal) modes; and (2) the field is tilted with respect to the vertical direction and the magnetic and acoustic modes become coupled giving magnetoacoustic waves. For the first case, the Dirac Hamiltonian possesses either a discrete spectrum, which corresponds to standing magnetic waves, or a continuous spectrum, which can be clearly identified with freely propagating acoustic waves. For the second case, the quantum mechanical perturbation calculus is used to study coupling and energy exchange between the magnetic and acoustic components of magnetoacoustic waves. It is shown that this coupling may efficiently prevent trapping of magnetoacoustic waves instellar atmospheres.

Alicki, R.↗

Continued RXTE Observations of Three Blazars

This project continued (during Cycle 3 of the RXTE mission) weekly observations of the quasars PKS 1510-089 and 3C 454.3 and the BL Lac object OJ 287. The data were analyzed by the PI and his graduate student. In the middle of Cycle 3, the PI changed the mode of observation of 3C 454.3, based on a report from another astronomer that observations with the Rosat satellite had found a strong confusing X-ray source about 45 arcmin away. The RXTE observations, in which the flux was measured while pointing directly at 3C 454.3 and then 15' from 3C 454.3, in the direction away from the confusing source. The flux decreased considerably, and so it appears that nearly all the X-ray flux measured by RXTE comes from the confusing source rather than the quasar. OJ 287 was detected only at one epoch of the Cycle 3 observations. This corresponded to a particularly quiet period for this source, when its flux at other wavebands was also low. From the single epoch of detection, no conclusions can be drawn concerning the general nature of the X-ray emission from this object. In contrast to the other two sources, the results of the observations of PKS 1510-089, when combined with data from cycles 3 and 4, are very interesting. The RXTE (2.4-20 keV) abd 14.5 GHz University of Michigan Radio Astronomy Observatory (UMRAO) light curves through the end of RXTE Cycle 4 are shown for PKS 1510-089. Visual inspection reveals an obvious correlation between the X-ray and 14.5 GHz light curves. While this is as expected in the SSC model in a source with simultaneous mm-wave (the photons thought to be scattered to X-ray energies) and cm-wave variations, the X-ray spectrum is flat ("energy" spectral index alpha = 0.2-0.5, where F (sub v) proportional to v (sup - alpha)) compared with the mm-wave synchrotron spectrum (alpha approximately equal to 0.8). The PI is currently engaged in theoretical work to determine whether this is possible for an SSC model.

Marscher, Alan P.↗

Preliminary investigation of acousto-ultrasonic evaluation of metal-matrix composite specimens

Acoustoultrasonic (AU) measurements were performed on a series of (0)8, (90)8, (+/-30)2s, (+/-60)2s, (0/90)2s, and (90/0)2s tensile specimens composed of eight laminated layers of continuous, SiC-fiber-reinforced Ti-15-3 matrix. The frequency spectrum was dominated by frequencies of longitudinal wave resonance through the thickness of the specimen at the sending transducer; signal propagation from the sender to the receiver was by shear waves. The magnitude of the frequency spectrum of the AU signal was used for calculating a stress-wave factor (SWF) based upon integrating the spectral distribution function. The SWF was sensitive to fiber-matrix debonding caused by mechanical strain and also to unconstrained thermal cycling. Damage in this metal-matrix composite structure due to tensile loading or thermal cycling was detected by this change in the magnitude of the frequency spectrum.

Kautz, Harold E.↗

Boson star normal modes

Boson stars are gravitationally bound objects that arise in ultralight dark matter models and form in the centers of galactic halos or axion miniclusters. We systematically study the excitations of a boson star, taking into account the mixing between positive and negative frequencies introduced by gravity. We show that the spectrum contains zero-energy modes in the monopole and dipole sectors resulting from spontaneous symmetry breaking by the boson star background. We analyze the general properties of the eigenmodes and derive their orthogonality and completeness conditions which have non-standard form due to the positive-negative frequency mixing. The eigenvalue problem is solved numerically for the first few energy levels in different multipole sectors and the results are compared to the solutions of the Schrödinger equation in fixed boson star gravitational potential. The two solutions differ significantly for the lowest modes, but get close for higher levels. We further confirm the normal mode spectrum in 3D wave simulations where we inject perturbations with different multipoles. As an application of the normal mode solutions, we compute the matrix element entering the evaporation rate of a boson star immersed in a hot axion gas. The computation combines the use of exact wavefunctions for the low-lying bound states and of the Schrödinger approximation for the high-energy excitations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Magnetic-field generation and its effect on ablative Rayleigh–Taylor instability in diffusive ablation fronts

The effects of self-generated magnetic fields on the ablative Rayleigh–Taylor (RT) instability are investigated in the linear regime. The main governing parameters are the Froude number (Fr), which stands for the ratio between ablative convection and acceleration of the target, and the Mach number at the ablation front (Ma), assumed to be small (isobaricity). During the development of the RT instability, magnetic fields are generated due to misalignment between pressure and density gradients (Biermann-battery effect). They accumulate at the section of the ablation front where the Nernst and the plasma velocities cancel each other. The magnetic field modifies the dynamics of the instability through the Righi–Leduc term, which acts as a heat source in the energy equation. It is found that the B fields affect perturbations with short wavelengths up to the most unstable wave in the spectrum. Here, the B field plays a destabilizing role for moderate Froude numbers and becomes stabilizing for large Froude numbers. For plastic ablators, the Fr threshold is found to be Fr = 5.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Excitation of toroidally localized harmonics of global Alfvén eigenmodes

Abstract The spherical tokamak NSTX and the upgraded, higher toroidal field, NSTX-U, routinely create plasmas with large un-thermalized populations of super-Alfvénic fast ions. This population of non-thermal energetic ions excites a broad spectrum of Alfvénic waves, from the lower frequency (≈30 kHz–≈200 kHz) toroidal Alfvén eigenmodes (Cheng and Chance 1986 Phys. Fluids 29 3695), to the higher frequency (400 kHz–3 MHz) global Alfvén eigenmodes (GAEs) and compressional Alfvén eigenmodes (Goedbloed 1975 Phys. Fluids 18 1258). In this paper we present evidence that the GAE non-linearly excite modes, presumably GAE, at frequencies consistent with non-linear or 3-wave coupling. The observation of the excitation of 2 nd harmonic GAE through the intrinsic non-linearity of Global Alfvén modes demonstrates that the non-linear terms can act as an exciter-antenna inside the plasma, broadcasting at harmonics of the mode frequency and with concomitantly shorter wavelengths. As with experiments using external antenna to excite otherwise weakly stable Alfvénic modes, this data can provide information about the stability of modes at harmonics of the GAE. It may also provide information on the nature of the non-linearities in the wave dispersion equation or a direct measure of the mode amplitude (Smith et al 2006 Phys. Plasmas 13 042504). We also report that the short wavelength 2 nd harmonic GAE can be strongly toroidally localized. The observation that the shorter wavelength ( n ≈ 20) 2 nd harmonic modes can be toroidally localized potentially has implications for the impact of fast-ion driven instabilities on fast-ion confinement in ITER and future fusion reactors (Gorelenkov et al 2014 Nucl. Fusion 54 125001).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Spectra of magnetic fields from electroweak symmetry breaking

We characterize magnetic fields produced during electroweak symmetry breaking by nondynamical numerical simulations based on the Kibble mechanism. The generated magnetic fields were thought to have an energy spectrum ∝ k3 for small wave numbers k, but here we show that it is actually a spectrum ∝ k4 along with characteristic fluctuations in the magnetic helicity. Using scaling results from magnetohydro- dynamics simulations for the evolution and assuming that the initial magnetic field is coherent on the electroweak Hubble scale, we estimate the magnetic field strength to be ∼10−13 G on kpc scales at the present epoch for nonhelical fields. For maximally helical fields we obtain ∼10−10 G on Mpc scales. We also give scalings of these estimates for partially helical fields.

Vachaspati, Tanmay [Arizona State University]↗

Nuclear Physics Exascale Requirements Review: An Office of Science Review sponsored jointly by Advanced Scientific Computing Research and Nuclear Physics, June 15 - 17, 2016, Gaithersburg, Maryland

Imagine being able to predict — with unprecedented accuracy and precision — the structure of the proton and neutron, and the forces between them, directly from the dynamics of quarks and gluons, and then using this information in calculations of the structure and reactions of atomic nuclei and of the properties of dense neutron stars (NSs). Also imagine discovering new and exotic states of matter, and new laws of nature, by being able to collect more experimental data than we dream possible today, analyzing it in real time to feed back into an experiment, and curating the data with full tracking capabilities and with fully distributed data mining capabilities. Making this vision a reality would improve basic scientific understanding, enabling us to precisely calculate, for example, the spectrum of gravity waves emitted during NS coalescence, and would have important societal applications in nuclear energy research, stockpile stewardship, and other areas. This review presents the components and characteristics of the exascale computing ecosystems necessary to realize this vision.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Exploring New Physics on Cosmological Scales

In recent work with her group, the PI developed a general formalism to compute from first principles the projected mass density (convergence) power spectrum of the substructure in galactic halos under different populations of dark matter sub halos. She constructed a halo model-based formalism, computing the 1-subhalo and the 2-subhalo terms from first principles for the first time. She found that the asymptotic slope of the substructure power spectrum at large wave number reflects the internal density profile of the sub halos, and proposed this as a key observable to discern between different dark matter scenarios.

79 ASTRONOMY AND ASTROPHYSICS↗

Quantum Monte Carlo Methods for Astrophysical Applications

In recent years, new astrophysical observations have provided a wealth of exciting input for nuclear physics. For example, the observations of two-solar-mass neutron stars put strong constraints on possible phase transitions to exotic phases in strongly interacting matter at high densities. Furthermore, the recent observation of a neutron-star merger in both the electromagnetic spectrum and gravitational waves has provided compelling evidence that neutron-star mergers are an important site for the production of extremely neutron-rich nuclei within the r-process. In the coming years, an abundance of new observations is expected, which will continue to provide crucial constraints on the nuclear physics of these events. To reliably analyze such astrophysical observations and extract information on nuclear physics, it is very important that a consistent approach to nuclear systems is used. Such an approach consists of a precise and accurate method to solve the nuclear many-body problem in nuclei and nuclear matter, combined with modern nuclear Hamiltonians that allow to estimate the theoretical uncertainties. Quantum Monte Carlo methods are ideally suited for such an approach and have been successfully used to describe atomic nuclei and nuclear matter. In this contribution, I will present a detailed description of Quantum Monte Carlo methods focusing on the application of these methods to astrophysical problems. In particular, I will discuss how to use Quantum Monte Carlo methods to describe nuclear matter of relevance to the physics of neutron stars.

79 ASTRONOMY AND ASTROPHYSICS↗

The ICON-A model for direct QBO simulations on GPUs (version icon-cscs:baf28a514)

Abstract. Classical numerical models for the global atmosphere, as used for numerical weather forecasting or climate research, have been developed for conventional central processing unit (CPU) architectures. This hinders the employment of such models on current top-performing supercomputers, which achieve their computing power with hybrid architectures, mostly using graphics processing units (GPUs). Thus also scientific applications of such models are restricted to the lesser computer power of CPUs. Here we present the development of a GPU-enabled version of the ICON atmosphere model (ICON-A), motivated by a research project on the quasi-biennial oscillation (QBO), a global-scale wind oscillation in the equatorial stratosphere that depends on a broad spectrum of atmospheric waves, which originates from tropical deep convection. Resolving the relevant scales, from a few kilometers to the size of the globe, is a formidable computational problem, which can only be realized now on top-performing supercomputers. This motivated porting ICON-A, in the specific configuration needed for the research project, in a first step to the GPU architecture of the Piz Daint computer at the Swiss National Supercomputing Centre and in a second step to the JUWELS Booster computer at the Forschungszentrum Jülich. On Piz Daint, the ported code achieves a single-node GPU vs. CPU speedup factor of 6.4 and allows for global experiments at a horizontal resolution of 5 km on 1024 computing nodes with 1 GPU per node with a turnover of 48 simulated days per day. On JUWELS Booster, the more modern hardware in combination with an upgraded code base allows for simulations at the same resolution on 128 computing nodes with 4 GPUs per node and a turnover of 133 simulated days per day. Additionally, the code still remains functional on CPUs, as is demonstrated by additional experiments on the Levante compute system at the German Climate Computing Center. While the application shows good weak scaling over the tested 16-fold increase in grid size and node count, making also higher resolved global simulations possible, the strong scaling on GPUs is relatively poor, which limits the options to increase turnover with more nodes. Initial experiments demonstrate that the ICON-A model can simulate downward-propagating QBO jets, which are driven by wave–mean flow interaction.

54 ENVIRONMENTAL SCIENCES↗