Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulated 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.

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

At least 721 records · Page 40

Modeling of wave-coherent pressures in the turbulent boundary layer above water waves

The behavior of air pressure fluctuations induced by progressive water waves generated mechanically in a laboratory tank was simulated by solving a modified Orr-Sommerfeld equation in a transformed Eulerian wave-following frame of reference. Solution is obtained by modeling the mean and wave-coherent turbulent Reynolds stresses, the behavior of which in the turbulent boundary layer above the waves was simulated using a turbulent kinetic energy-dissipation model, properly modified to account for free-surface proximity and favorable pressure gradient effects. The distribution of both the wave-coherent turbulent Reynolds stress and pressure amplitudes and their corresponding phase lags was found to agree reasonably well with available laboratory data.

Papadimitrakis, Yiannis ALEX.↗

Spectral broadening from turbulence in multiscale lower hybrid current drive simulations

Abstract The scattering of lower hybrid (LH) waves due to scrape-off layer filaments is investigated. It is revealed that scattering can account for the LH spectral gap without any ad hoc modification to the wave-spectrum. This is shown using a multiscale simulation approach which allows, for the first time, the inclusion of full-wave scattering physics in ray-tracing/Fokker–Planck calculations. In this approach, full-wave scattering probabilities are calculated for a wave interacting with a statistical ensemble of filaments. These probabilities are coupled to ray-tracing equations using radiative transfer theory. This allows the modeling of scattering along the entire ray-trajectory, which can be important in the multi-pass regime. Simulations are conducted for LH current drive (LHCD) in Alcator C-Mod, resulting in excellent agreement with experimental current and hard x-ray profiles. A region in filament parameter space is identified in which the impact of scattering on LHCD is saturated. Such a state coincides with experimental LHCD measurements, suggesting saturation indeed occurs in C-Mod, and therefore the exact statistical properties of the filaments are not important.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A comparative study of three techniques for using the spectral matrix in wave analysis

Three techniques - power spectral analysis, coherency analysis, and wave analysis - for studying ULF waves in the magnetosphere are described. These techniques, while mathematically dissimilar, are all based on the assumption that, at any single frequency, a single plane wave, propagating along the normal to the plane, is present. Wave polarization parameters provided by these techniques are important in the development of models of wave generation. A comparative study of the three techniques was carried out using simulated data involving known wave and noise properties and real ULF event data from ATS 1 and 6.

Arthur, C. W.↗

WEC-GRID

WEC-GRID is an open-source Python library crafted to simulate the integration of Wave Energy Converters (WECs) and Current Energy Converters (CECs) into renowned power grid simulators like PSS®E & PyPSA.

Barajas-Ritchie, Alexander↗

Estimating gravity wave parameters from oblique high-frequency backscatter: Modeling and analysis

A new technique for estimating electron density perturbation amplitudes of traveling ionospheric disturbances (TIDs), using HF radar data, is presented. TIDs are observed in HF radar data as enhancements of the ground-scattered power which propagate through the radar's field of view. These TIDs are the ionospheric manifestation of atmospheric acoustic-gravity waves. TID electron density perturbation amplitudes were estimated by simulating the radar returns, using HF ray tracing through a model ionosphere perturbed by a model gravity wave. The simulation determined the return power in the ground-scattered portion of the signal as a function of range, and this was compared to HF radar data from the Goose Bay HF radar at a time when evidence of gravity waves was present in the data. By varying the amplitude of the electron density perturbation in the model it was possible to estimate the perturbation of the actual wave. It was found that the perturbations that are observed by the Goose Bay HF radar are of the order of 20% to 35%. It was also found that the number of observable power enhancements, and the relative amplitudes of these enhancements, depended on the vertical thickness of the gravity wave's source region. From the simulations and observations it was estimated that the source region for the case presented here was approximately 20 km thick. In addition, the energy in the wave packet was calculated and compared to an estimate of the available energy in the source region. It was found that the wave energy was about 0.2% of the estimated available source region energy.

Bristow, W. A.↗

Collaborative Research: Parametric Instabilities of Alfven Waves in Low‐beta Plasmas

Alfven waves are of fundamental importance in magnetized plasmas. This project aims to advance our understanding of the basic physics of the parametric instabilities of Alfven waves using advanced computer simulations. The results will have implications for important processes in the heliosphere involving Alfven waves, such as coronal heating and solar wind heating.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Heart Fibrillation and Parallel Supercomputers

The Luo and Rudy 3 cardiac cell mathematical model is implemented on the parallel supercomputer CRAY - T3D. The splitting algorithm combined with variable time step and an explicit method of integration provide reasonable solution times and almost perfect scaling for rectilinear wave propagation. The computer simulation makes it possible to observe new phenomena: the break-up of spiral waves caused by intracellular calcium and dynamics and the non-uniformity of the calcium distribution in space during the onset of the spiral wave.

Kogan, B. Y.↗

The effect of SST and soil moisture anomalies on GLA model simulations of the 1988 U.S. summer drought

A series of simulations of the late spring and early summer of 1988 were conducted in order to study the relative importance of different boundary forcings to the Goddard Laboratory for Atmospheres (GLA) model's simulation of the heat wave and drought over the Great Plains of the United States during this time period. Separate 60-day simulations were generated from 10, 20, and 30 May 1988 with a variety of boundary condition datasets. For the control experiment, climatological boundary conditions were used. This was followed by experiments in which either the observed 1988 sea surface temperatures (SST) or derived 1988 soil moisture values, or both, were used in place of the climatological fields. Additional experiments were conducted in which only tropical or midlatitude SST anomalies were used. The impact of the different boundary forcings was evaluated relative to the control simulations of the precipitation and surface air temperature over the Great Plains. It was found that the tropical SST anomalies had a significant effect in reducing precipitation in this area, while the midlatitude anomalies did not. Due to the prescribed climatological soil moistures for the SST experiments, a significant increase in surface temperature did not occur in these simulations. In contrast, the simulations with the anomalous 1988 soil moistures produced both a larger reduction of precipitation and a significant increase in surface temperature over the Great Plains. The simulations with both anomalous SST and soil moisture showed only a slight augmentation of the heat wave and drought relative to the experiments with anomalous soil moisture alone.

Atlas, R.↗

Kinetic properties of mirror waves in magnetosheath plasmas

Linear and nonlinear properties of waves excited by the mirror instability in high beta, low anisotropy plasmas characteristic of the magnetosheath are investigated using linear theory and one-dimensional hybrid simulations. The mechanisms for wave growth and saturation at low amplitudes are discussed. A new method is considered for generating the large amplitude mirror waves observed in the magnetosheath based on external compression of magnetic flux tubes. Simulations in which the anisotropy is maintained by recycling the ions shows this process can inhibit the growth of ion cyclotron waves and enhances the growth of mirror waves.

Mckean, M. E.↗

Scattering wave packets of hadrons in gauge theories: Preparation on a quantum computer

Quantum simulation holds promise of enabling a complete description of high-energy scattering processes rooted in gauge theories of the Standard Model. A first step in such simulations is preparation of interacting hadronic wave packets. To create the wave packets, one typically resorts to adiabatic evolution to bridge between wave packets in the free theory and those in the interacting theory, rendering the simulation resource intensive. In this work, we construct a wave-packet creation operator directly in the interacting theory to circumvent adiabatic evolution, taking advantage of resource-efficient schemes for ground-state preparation, such as variational quantum eigensolvers. By means of an ansatz for bound mesonic excitations in confining gauge theories, which is subsequently optimized using classical or quantum methods, we show that interacting mesonic wave packets can be created efficiently and accurately using digital quantum algorithms that we develop. Specifically, we obtain high-fidelity mesonic wave packets in the Z 2 and U(1) lattice gauge theories coupled to fermionic matter in 1+1 dimensions. Our method is applicable to both perturbative and non-perturbative regimes of couplings. The wave-packet creation circuit for the case of the Z 2 lattice gauge theory is built and implemented on the Quantinuum H1-1 trapped-ion quantum computer using 13 qubits and up to 308 entangling gates. The fidelities agree well with classical benchmark calculations after employing a simple symmetry-based noise-mitigation technique. This work serves as a step toward quantum computing scattering processes in quantum chromodynamics.

97 MATHEMATICS AND COMPUTING↗

Atmospheric Energy Spectra in Global Kilometre-Scale Models

Eleven 40-day long integrations of five different global models with horizontal resolutions of less than 9 km are compared in terms of their global energy spectra. The method of normal-mode function decomposition is used to distinguish between balanced (Rossby wave; RW) and unbalanced (inertia-gravity wave; IGW) circulation. The simulations produce the expected canonical shape of the spectra, but their spectral slopes at mesoscales, and the zonal scale at which RW and IGW spectra intersect differ significantly. The partitioning of total wave energies into RWs an IGWs is most sensitive to the turbulence closure scheme and this partitioning is what determines the spectral crossing scale in the simulations, which differs by a factor of up to two. It implies that care must be taken when using simple spatial filtering to compare gravity wave phenomena in storm-resolving simulations, even when the model horizontal resolutions are similar. In contrast to the energy partitioning between the RWs and IGWs, changes in turbulence closure schemes do not seem to strongly affect spectral slopes, which only exhibit major differences at mesoscales. Despite their minor contribution to the global (horizontal kinetic plus potential available) energy, small scales are important for driving the global mean circulation. Our results support the conclusions of previous studies that the strength of convection is a relevant factor for explaining discrepancies in the energies at small scales. The models studied here produce the major large-scale features of tropical precipitation patterns. However, particularly at large horizontal wavenumbers, the spectra of upper tropospheric vertical velocity, which is a good indicator for the strength of deep convection, differ by factors of three or more in energy. High vertical kinetic energies at small scales are mostly found in those models that do not use any convective parameterisation.

DYAMOND↗

Field structure of collapsing wave packets in 3D strong Langmuir turbulence

A simple model is constructed for the electric fields in the collapsing wave packets found in 3D simulations of driven and damped isotropic strong Langmuir turbulence. This model, based on a spherical-harmonic decomposition of the electrostatic potential, accounts for the distribution of wave-packet shapes observed in the simulations, particularly the predominance of oblate wave packets. In contrast with predictions for undamped and undriven subsonic collapse of scalar fields, oblate vector-field wave packets do not flatten during collapse but, instead, remain approximately self-similar and rigid.

Newman, D. L.↗

Millimeter Wave Doppler Radar Testing Using the Lunar Regolith Ejecta Simulator

A linear conveyor system has been devised to maintain lunar regolith simulant at a range of velocities, distances, angles, particle densities, and particle sizes. The Lunar Regolith Ejecta Simulator (LuRES) system shall test prototype and potential flight instrumentation for measuring the velocities of regolith particles accelerated by a plume-surface interaction (PSI) and determining effects of particle number densities and particle size distributions on return signals. Here we present the results of initial LuRES testing of the Millimeter Wave Doppler Radar (MWDR) instrument. The LuRES system held a single, ~900 µm BP-1 lunar regolith simulant particle at a continuous ~5 m/s velocity. The return signal of the MWDR from this LuRES configuration was analyzed to calibrate the MWDR instrument. Further testing is planned for the remaining of 2023 to calibrate the MWDR for various particle sizes and densities at a continuous velocity. These results will be presented in the final paper.

PSI↗

The Role of Vapor Pressure Deficit in the CLM Simulated Interaction Between Urban Heat Islands and Heat Waves Over CONUS

Abstract Although both urban and rural temperatures are expected to increase under heat waves (HWs), whether the urban heat island (UHI) intensity becomes stronger under HWs remains unknown especially at the daily mean and large spatial scales. Using an urbanized land surface model, we quantify the interactions between UHIs and HWs over the Contiguous United States (CONUS). Synergistic interactions (i.e., increased UHI intensities under HWs) are observed over the eastern and western U.S. However, negative interactions are found in the Central U.S. due to the stronger inhibition of rural evapotranspiration by vapor pressure deficit (VPD) stresses. The interactions between UHIs and HWs in the Central U.S. will be further reduced along with the elevated VPD stresses in a hotter future. The results highlight the importance of properly parameterizing the sensitivity of urban and rural evapotranspiration to various environmental stresses in climate and earth system models.

Wang, Linying [National Key Laboratory of Earth Sy↗

Simulations for x-ray imaging of wave-packet dynamics

Previous work on imaging wave-packet dynamics with x-ray scattering revealed that the scattering patterns deviate substantially from the notion of instantaneous momentum density of the wave packet. Here we show that scattering patterns can provide clear insights into the electron wave-packet dynamics if the final state of the scattered electron and the scattered photon momentum are determined simultaneously. The scattering probability is shown to be proportional to the modulus square of the Fourier transform of the instantaneous electronic spatial wave function weighted by the final state of the electron. Several cases for the choice of final state of the electron are explored. First, the case where the final state can be measured up to a given principal quantum number n and orbital angular momentum l are presented. Next, the case where the final states can only be determined up to a given energy is discussed. Lastly, the case of an initial wave packet consisting of a large amount of a known stationary state and a small amount of an unknown stationary state is examined. The scattering profile is used to determine the properties of the unknown state in the wave packet.

74 ATOMIC AND MOLECULAR PHYSICS↗

Simulation of the inertial-conductive subrange

The inertial-conductive subrange spectrum of a passive temperature field of a fluid of small Prandtl number is determined by large-eddy simulation. Results for simulations of both freely decaying and forced turbulence are presented. In the simulations of freely decaying turbulence, a subgrid model is used to simulate a decaying turbulent velocity field with a k exp -5/3 inertial subrange energy spectrum convecting eight different decaying temperature fields with well-resolved conductive subranges. In the simulations of forced turbulence, a subgrid model is again used to simulate an inertial subrange velocity field; however, the velocity field is now forced at the lowest wave numbers of the simulation, and an external uniform mean temperature gradient is imposed. Statistically stationary velocity and temperature fluctuations are generated. The results of the decaying and forced simulations are in excellent agreement with the Batchelor et al. (1959) k exp -17/3 spectrum in the far inertial-conductive subrange, whereas significant departures from the Batchelor spectrum are observed in the near inertial-conductive subrange.

Chasnov, Jeffrey R.↗

Simulation of Blunt-Fin-Induced Shock-Wave and Turbulent Boundary-Layer Interaction

The Reynolds-averaged Navier-Stokes equations are solved numerically for supersonic flow over a blunt fin mounted on a flat plate. The fin shock causes the boundary layer to separate. which results in a complicated, three-dimensional shock-wave and boundary-layer interaction. The computed results are in good agreement with the mean static pressure measured on the fin and the flat plate. The main features, such as peak pressure on the fin leading edge and a double peak on the plate. are predicted well. The role of the horseshoe vortex is discussed. This vortex leads to the development of high-speed flow and, hence, low-pressure regions on the fin and the plate. Different thicknesses of the incoming boundary layer have been studied. Varying the thicknesses by an order of magnitude shows that the size of the horseshoe vortex and, therefore, the spatial extent of the interaction are dominated by inviscid flow and only weakly dependent on the Reynolds number. Colored graphics are used to show details of the interaction flow field.

Hung, Ching-Mao↗