Theoretical models and numerical studies of waves in a three-fluid medium
Mathematical models of wave propagation in three-fluid plasma and numerical solutions to dispersion and propagation properties
Murphree, D. L.↗
SEARCH · Engineering Papers
Results for “wave model”
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.
DETERMINATION OF THE DETONATION WAVE STRUCTURE
Laminar detonation wave structure in ozone for three models - coupled wave, von neuman-doring- zeldovich model, and decoupled wave
DETONATION WAVE↗
Equatorial waves simulated by the NCAR community climate model
The equatorial planetary waves simulated by the NCAR CCM1 general circulation model were investigated in terms of space-time spectral analysis (Kao, 1968; Hayashi, 1971, 1973) and energetic analysis (Hayashi, 1980). These analyses are particularly applied to grid-point data on latitude circles. In order to test some physical factors which may affect the generation of tropical transient planetary waves, three different model simulations with the CCM1 (the control, the no-mountain, and the no-cloud experiments) were analyzed.
Cheng, Xinhua↗
Long-term Global Morphology of Gravity Wave Activity Using UARS Data
An extensive body of research this quarter is documented. Further methodical analysis of temperature residuals in Cryogenic Limb Array Etalon Spectrometer (CLAES) Version 8 level 3AT data show signatures during December 1992 at middle and high northern latitudes that, when compared to Naval Research Laboratory/Mountain Wave Forecast Model (NRL)/(MWFM) mountain wave hindcasts, reveal evidence of long mountain waves in these data over Eurasia, Greenland, Scandinavia and North America. The explicit detection of gravity waves in limb-scanned Cryogenic Infrared Spectrometers and Telescopes for the Atmosphere (CRISTA) temperatures is modeled at length, to derive visibility functions. These insights are used to convert CRISTA gravity wave temperature residuals into data that more closely resemble gravity wave fluctuations detected in data from other satellite instruments, such as Microwave Limb Sounder (MLS), Limb Infrared Monitor of the Stratosphere (LIMS) and Global Positioning System/Meteorology (GPS)/(MET). Finally, newly issued mesospheric temperatures from inversion of CRISTA 15gin emissions are analyzed using a new method that uses separate Kalman fits to the ascending and descending node data. This allows us to study global gravity wave amplitudes at two local times, 12 hours apart. In the equatorial mesosphere, where a large diurnal tidal temperature signal exists, we see modulations of gravity wave activity that are consistent with gravity wave-tidal interactions produced by tidal temperature variability.
Eckermann, Stephen D.↗
Prediction of Complex Observed Shear Wave Splitting Patterns at Ryukyu Subduction Zone Using a Strong Intra-Slab Anisotropy Model
Complex shear wave splitting (SWS) patterns in subduction zones are often interpreted geodynamically as resulting from complex mantle flow; however, this may not always be necessary. We analyzed 7,093 high-quality SWS measurements from teleseismic S waves recorded by Hi-net stations across the Ryukyu arc in Japan. Our findings show a systematic rotation of the fast S polarization from trench-parallel to trench-perpendicular depending on the earthquake backazimuth. For the same earthquake, the measured splitting patterns also vary spatially across the southwest Japan. Using full-wave seismic modeling, we showed that a dipping slab with ~30% shear anisotropy of the tilted transverse isotropy (TTI) type, with a symmetry axis perpendicular to the slab interface, can predict the observed delay times and polarization rotation. Our results highlight the importance of considering dipping anisotropic slabs in interpreting SWS at subduction zones.
58 GEOSCIENCES↗
A comparison between numerically modelled and experimentally measured loss mechanisms in wave rotors
A numerical model has been developed which is capable of predicting the performance of a wave rotor (pressure exchanger) of specified geometry over a wide range of operating conditions. The model can account for the major loss mechanisms of leakage from the tube ends, fluid viscosity, heat transfer to the tube wails, finite tube opening time, shock waves, and non-uniform port flows. It is a one dimensional flow model which follows a single tube as it rotates past the various stationary ports. Since the model is relatively simple (i.e., one dimensional) it uses little computer time. This makes it suitable for design as well as analytical purposes. This paper will present a brief description of the model then discuss a comparison between the model predictions and several wave rotor experiments.
Paxson, Daniel E.↗
A Comparison Between Numerically Modelled and Experimentally Measured Loss Mechanisms in Wave Rotors
A numerical model has been developed which is capable of predicting the performance of a wave rotor (pressure exchanger) of specified geometry over a wide range of operating conditions. The model can account for the major loss mechanisms of leakage from the tube ends, fluid viscosity, heat transfer to the tube walls, finite tube opening time, shock waves, and non-uniform port flows. It is a one dimensional flow model which follows a single tube as it rotates past the various stationary ports. Since the model is relatively simple (i.e. one dimensional) it uses little computer time. This makes it suitable for design as well as analytical purposes. This paper will present a brief description of the model then discuss a comparison between the model predictions and several wave rotor experiments.
Paxson, Daniel E.↗
Implementation of a Three-Dimensional Planetary Boundary Layer Parameterization in a Coupled Modeling System and Evaluation of "Gray Zone" Simulations of a Wind-Wave Event Off the US California Coast Using Observations
A three-dimensional (3D) planetary boundary layer (PBL) parameterization was added to the Coupled Ocean Atmosphere Wave Sediment Transport (COAWST) model and the first coupled atmosphere/wave "gray zone" simulations were carried out using the 3D PBL parameterization. A case study of a significant wind-wave event from 4-8 November 2020 off the US west coast was examined, with a focus on the impacts of the event within an approximate 280x280 km central California coastal region covering the Bureau of Ocean Energy Management Morro Bay wind energy lease area. Simulations with both one-dimensional (1D; Yonsei University) and 3D PBL parameterizations were examined in this nest with ..delta..x=400 m. Two-way coupling was active, with near-surface winds feeding back to the wave model and bulk wave statistics feeding back to the atmospheric model. Both simulations compared favorably with buoy observations in capturing the timing and magnitude of wind speed, temperature, dewpoint, and significant wave height, as the front associated with the maritime weather system moved southward across coastal California. However, wind speed errors over land were larger. Time series of the vertical profile of winds below z=250 m from the simulations compared favorably with observations from the US Dept. of Energy Morro Bay Doppler lidar buoy. Differences between the coupled 1D and 3D PBL simulations were minor. The most notable difference was that the simulation using the 3D PBL parameterization had approximately 10 % weaker winds at the peak of the event causing 10 % lower significant wave heights. This difference was shown to be mostly due to differences in the vertical mixing treatment between the parameterizations, in particular the nonlocal downward mixing of higher momentum air in the 1D Yonsei University parameterization in convective conditions. Overall, this work demonstrates that the 3D PBL parameterization can be used in a coupled atmosphere/wave modeling framework with similar behaviors as traditional PBL parameterizations that do not provide horizontal turbulent variances and fluxes.
16 TIDAL AND WAVE POWER↗
Steady-state turbulence with a narrow inertial range
Coupled two-dimensional wave equations are solved on a computer to model Langmuir wave turbulence excited by a weak electron beam. The model includes wave growth due to beam-plasma interaction, and dissipation by Landau damping. The inertial range is limited to a relatively small number of modes such as could occur when the ratio of masses between the negative and positive ions is larger than in a hydrogen plasma, or when there is damping in long wavelength Langmuir waves. A steady state is found consisting of quasistable, collapsed wave packets. The effects of different beam parameters and the assumed narrow inertial range are considered. The results may be relevant to plasma turbulence observed in connection with type III solar bursts.
Weatherall, J. C.↗
Irreversible wave-mean flow interactions in a mechanistic model of the stratosphere
A wave-1 minor warming is simulated in a mechanistic, global, primitive equation model. The modification of the zonal flow by the rearrangement of potential vorticity on a middle stratopheric isentrope is compared in a fully nonlinear model and in a model with only one wave and the zonal flow (a quasi-linear model). The permanent rearrangement of potential vorticity during the wave episode is more intense and more localized meridionally in the fully nonlinear model, which is able to capture the process of planetary wave breaking in some detail. Additional experiments reveal that the differences between the quasi-linear and nonlinear models persist for a broad range of wave amplitudes, and that the quasi-linear model can qualitatively reproduce the modification of the zonal flow by the wave when the diffusive dissipation of the wave is enhanced. These results are discussed in the context of the theory of barotropic Rossby waves in shear flows, and in comparison with recent numerical simulations of the middle atmosphere.
Robinson, Walter A.↗
Origin of Intense Electron Heating in Relativistic Blast Waves
The modeling of gamma-ray burst afterglow emission bears witness to strong electron heating in the precursor of Weibel-mediated, relativistic collisionless shock waves propagating in unmagnetized electron–ion plasmas. In this Letter, we propose a theoretical model, which describes electron heating via a Joule-like process caused by pitch-angle scattering in the decelerating, self-induced microturbulence and the coherent charge-separation field induced by the difference in inertia between electrons and ions. The emergence of this electric field across the precursor of electron–ion shocks is confirmed by large-scale particle-in-cell (PIC) simulations. Integrating the model using a Monte Carlo-Poisson method, we compare the main observables to the PIC simulations to conclude that the above mechanism can indeed account for the bulk of electron heating.
79 ASTRONOMY AND ASTROPHYSICS↗
Hydrostatic calculations of axisymmetric flow and its stability for the AGCE model
Baroclinic waves in the atmospherics general circulation experiment (AGCE) apparatus by the use of numerical hydrostatic primitive equation models were determined. The calculation is accomplished by using an axisymmetric primitive equation model to compute, for a given set of experimental parameters, a steady state axisymmetric flow and then testing this axisymmetric flow for stability using a linear primitive equation model. Some axisymmetric flows are presented together with preliminary stability calculations.
Miller, T. L.↗
A numerical model study of long-term planetary wave predictability
A numerical model was developed to study planetary wave behavior over periods up to 4 months under the influences of stochastic and diffusive forcing by synoptic scale motions and of sea surface temperature anomalies. The model is suitable for Monte Carlo experiments which provide ensemble mean and variance predictions and permit the response to heating anomalies to be partly separated from random variability. The model incorporates the linear-balanced dynamical equations for wind at two tropospheric levels and temperature at one level. The model atmosphere is forced by semiempirical January heating linearly dependent on local air temperature. Topography and the lower boundary layer are represented by vertical motion and drag functions, respectively, of the lower level wind. The linear phase speeds of Rossby waves are modified to reflect the effect of the stratosphere.
Wobus, R. L.↗
Shock wave structure using nonlinear model Boltzmann equations
The structure of a strong plane shock wave in a monatomic rarefied perfect gas is one of the simplest problems able to be posed in kinetic theory, and one of the hardest to solve. Its simplicity lies in the absence of solid boundaries, geometrical complications, or internal molecular energy. Its difficulty arises from the great departure of the gas from equilibrium within the shock, which invalidates many of the techniques used successfully elsewhere in kinetic theory. In addition to this theoretical challenge, the modern development of ballistics and hypersonic flight has helped to stimulate extensive theoretical and experimental interest in the shock problem. The experimenters in turn have encountered great difficulties on account of the very small physical dimensions of shocks. In fact, until very recently indeed, any close comparisons of theoretical and experimental shock structure results have been rather unprofitable due to the inadequacies of both theory and experiment. During the last few years this situation has been appreciably improved by development of the Monte Carlo method. This allows idealized 'experiments' to be performed on large computers instead of in wind tunnels, using a known intermolecular force law. The most developed of these methods has been shown to be equivalent theoretically to the Boltzmann equation and to give results which agree extremely closely with measurements of high accuracy. Thus Monte Carlo results not only form the soundest basis for our present theoretical knowledge of shock wave structure, but, for purposes of developing other theories, can also be considered a very valuable experimental resource. However, such results remain very expensive to obtain. In this thesis we develop more economical kinetic theory methods for the approximate prediction of shock structure, and compare our results with those of the Monte Carlo method.
Segal, Ben Maurice↗
FullWave — A Full Wave Parallel Code for Modeling RF Fields in Hot Tokamak Plasma
FullWave is a computer code that simulates how radio-frequency (RF) waves travel and deposit energy in the hot plasma inside a fusion reactor. RF waves are used to heat the plasma and drive electrical current, which is essential for sustaining fusion reactions. The code uses a new algorithm that can handle much finer spatial detail than previous codes — more than 100 times finer — while running efficiently on national supercomputers. It incorporates a detailed physics model that captures subtle kinetic effects important for accurate prediction of wave behavior. Under this project, FullWave was extended to cover multiple RF frequency ranges relevant to present and future tokamaks, and validated against experimental parameters from the DIII-D tokamak at General Atomics. Results were published in peer-reviewed journal articles.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Role of the likelihood for elastic scattering uncertainty quantification
In the last decade, uncertainty quantification (UQ) for optical model potentials (OMPs) has become a focal point for nuclear reaction theory, and several competing approaches for OMP UQ have recently been developed. Here, we clarify recent efforts to compare frequentist and Bayesian approaches in the context of OMP UQ [G. B. King et al., Phys. Rev. Lett. 122, 232502 (2019)]. We replicate a portion of that OMP UQ study but use independent statistical tools. Specifically, we compare two methods for OMP parameter inference from elastic scattering data: the Levenberg-Marquardt algorithm for χ 2 minimization on one hand and Markov chain Monte Carlo (MCMC) sampling on the other. Separately, we assess the common practice of using a renormalized likelihood (χ 2 /N), N being the number of data points, instead of the canonical weighted-least-squares likelihood (χ 2 ), as a way of accounting for unknown data correlations. Here, we show that for a generic linear model and for a five-parameter OMP analysis, frequentist and uniform-prior Bayesian approaches recover the same optimum and uncertainty estimates—not systematically larger uncertainties for the Bayesian approach, as was concluded in G. B. King et al., Phys. Rev. Lett. 122, 232502 (2019). Further, we show that if an additional, near-degenerate parameter is introduced into the same OMP analysis such that the parameter posterior becomes non-Gaussian, then covariance-based estimates of uncertainty become unreliable. Finally, we show that regardless of optimization approach, if χ 2 /N is used for the likelihood, the resulting parametric uncertainties increase by $\sqrt{N}$, and that this is responsible for the conclusions drawn in the revisited study. Based on our replication results, we find that a fortuitous cancellation of unreported errors and the renormalization factor can lead to improvement in empirical coverages, as was the case in the original comparative study. We emphasize that developing and applying a realistic likelihood function is an essential task in a UQ analysis, and that several recent UQ studies that employed a renormalized likelihood (i.e., including a 1/N factor) may have yielded unrealistically large uncertainties for elastic-scattering observables. If the parameter posterior deviates from multivariate-normal, a sampling-based approach like MCMC has a clear advantage over methods that assume the Laplace approximation holds. We note that empirical coverage can serve as an important internal check for the analyst whose model or data may have additional, unaccounted-for uncertainties.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
The motion of ions specularly reflected off a quasi-parallel shock in the presence of large-amplitude, monochromatic MHD waves
A model is used to examine the motion of specularly reflected ions in the presence of large-amplitude, monochromatic, transverse MHD waves. The calculations of ion trajectories are described. The heating downstream from the quasi-parallel bow shock is analyzed. The relationship between the specularly reflected ions and their gyrospeeds and guiding center speeds is studied. The data reveal that the characteristics of the motion depend on the frequency, wavelength, phase, and amplitude of the wave that is converted into the shock.
Fuselier, S. A.↗
Reduced dimer production in solar-simulator-pumped continuous wave iodine lasers based on model simulations and scaling and pumping studies
A numerical rate equation model for a continuous wave iodine laser with longitudinally flowing gaseous lasant is validated by approximating two experiments that compare the perfluoroalkyl iodine lasants n-C3F7I and t-C4F9I. The salient feature of the simulations is that the production rate of the dimer (C4F9)2 is reduced by one order of magnitude relative to the dimer (C3F7)2. The model is then used to investigate the kinetic effects of this reduced dimer production, especially how it improves output power. Related parametric and scaling studies are also presented. When dimer production is reduced, more monomer radicals (t-C4F9) are available to combine with iodine ions, thus enhancing depletion of the laser lower level and reducing buildup of the principal quencher, molecular iodine. Fewer iodine molecules result in fewer downward transitions from quenching and more transitions from stimulated emission of lasing photons. Enhanced depletion of the lower level reduces the absorption of lasing photons. The combined result is more lasing photons and proportionally increased output power.
Costen, Robert C.↗