Simulated variability of the high‐latitude thermosphere induced by small‐scale gravity waves during a sudden stratospheric warming
Explore the source record for details and available documents.
Erdal Yiğit↗
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.
Laboratory simulation of pelagic asteroidal impact Atmospheric injection, benthic topography, and the surface wave radiation field
The hypervelocity impact (1.25-6 km/s) of projectiles into water overlying unconsolidated strata is reported for a variety of water depths. Variation of the background atmospheric pressure is found to be an important additional parameter. The relation of these exploratory experiments to pelagic impact of asteroidal-sized objects is discussed from the standpoint of atmospheric injection of sea water, the modification of bottom (benthic) topography, and the generation of a field of mega-amplitude sea waves and their propagation away from the source.
Gault, D. E.↗
Numerical simulation and analysis of a prefrontal squall line. I - Observations and basic simulation results. II - Propagation of the squall line as an internal gravity wave
Results of an observational and numerical study of the squall line that occurred on June 17-18, 1978, are presented. This squall line was initially triggered by the strong surface convergence along a cold front and stretched from Illinois to the Texas Panhandle. The squall line was aligned with the surface front during its initial development, but then propagated faster than the front, resulting in a separation of approximately 200 km by 0300 UTC and 300-400 km by 0600 UTC. The Colorado State University Regional Atmospheric Modeling System is used to model the squall-line development and propagation. Results are described from several experiments that tested the sensitivity to the use of the Kuo-type cumulus parameterization scheme and grid-scale microphysical processes. The movement of the squall line in the model is shown to be due to the propagation of a deep tropospheric internal gravity wave in a CISK-like process. The thermal and dynamic perturbations associated with the hypothesized wave are shown to be consistent with internal gravity wave theory.
Cram, Jennifer M.↗
Gravity Waves and Mesospheric Clouds in the Summer Middle Atmosphere: A Comparison of Lidar Measurements and Ray Modeling of Gravity Waves Over Sondrestrom, Greenland
We conducted gravity wave ray-tracing experiments within an atmospheric region centered near the ARCLITE lidar system at Sondrestrom, Greenland (67N, 310 deg E), in efforts to understand lidar observations of both upper stratospheric gravity wave activity and mesospheric clouds during August 1996 and the summer of 2001. The ray model was used to trace gravity waves through realistic three-dimensional daily-varying background atmospheres in the region, based on forecasts and analyses in the troposphere and stratosphere and climatologies higher up. Reverse ray tracing based on upper stratospheric lidar observations at Sondrestrom was also used to try to objectively identify wave source regions in the troposphere. A source spectrum specified by reverse ray tracing experiments in early August 1996 (when atmospheric flow patterns produced enhanced transmission of waves into the upper stratosphere) yielded model results throughout the remainder of August 1996 that agreed best with the lidar observations. The model also simulated increased vertical group propagation of waves between 40 km and 80 km due to intensifying mean easterlies, which allowed many of the gravity waves observed at 40 km over Sondrestrom to propagate quasi-vertically from 40-80 km and then interact with any mesospheric clouds at 80 km near Sondrestrom, supporting earlier experimentally-inferred correlations between upper stratospheric gravity wave activity and mesospheric cloud backscatter from Sondrestrom lidar observations. A pilot experiment of real-time runs with the model in 2001 using weather forecast data as a low-level background produced less agreement with lidar observations. We believe this is due to limitations in our specified tropospheric source spectrum, the use of climatological winds and temperatures in the upper stratosphere and mesosphere, and missing lidar data from important time periods.
LIDAR (LIGHT DETECTION AND RANGING)↗
Separating Gravitational Wave Signals from Instrument Artifacts
Central to the gravitational wave detection problem is the challenge of separating features in the data produced by astrophysical sources from features produced by the detector. Matched filtering provides an optimal solution for Gaussian noise, but in practice, transient noise excursions or "glitches" complicate the analysis. Detector diagnostics and coincidence tests can be used to veto many glitches which may otherwise be misinterpreted as gravitational wave signals. The glitches that remain can lead to long tails in the matched filter search statistics and drive up the detection threshold. Here we describe a Bayesian approach that incorporates a more realistic model for the instrument noise allowing for fluctuating noise levels that vary independently across frequency bands, and deterministic "glitch fitting" using wavelets as "glitch templates", the number of which is determined by a trans-dimensional Markov chain Monte Carlo algorithm. We demonstrate the method's effectiveness on simulated data containing low amplitude gravitational wave signals from inspiraling binary black hole systems, and simulated non-stationary and non-Gaussian noise comprised of a Gaussian component with the standard LIGO/Virgo spectrum, and injected glitches of various amplitude, prevalence, and variety. Glitch fitting allows us to detect significantly weaker signals than standard techniques.
Littenberg, Tyson B.↗
Revisiting Recent U.S. Heat Waves in a Warmer and More Humid Climate
The frequency and intensity of heat waves in the United States is projected to increase in the 21st century. We investigate dry and humid heat waves in a pair of high-resolution model simulations that constrain large-scale atmospheric circulation, to isolate the thermodynamic impacts on characteristics of present and future heat waves over the United States. The two kinds of heat waves show differences in mean intensity, amplitude, duration, and frequency over the Southeast, Northeast, and Midwest, while their characteristics are largely similar in the drier central and western United States. In a warmer climate, relative humidity is projected to decrease during dry heat waves, whereas it remains unchanged during humid heat waves. However, the overall increase in daily maximum temperature intensifies the heat stress during future humid and dry heat waves across all regions. With large-scale circulation constrained, these simulations emphasize the importance of thermodynamic drivers in determining future heat wave characteristics.
54 ENVIRONMENTAL SCIENCES↗
Results from Binary Black Hole Simulations in Astrophysics Applications
Present and planned gravitational wave observatories are opening a new astronomical window to the sky. A key source of gravitational waves is the merger of two black holes. The Laser Interferometer Space Antenna (LISA), in particular, is expected to observe these events with signal-to-noise ratio's in the thousands. To fully reap the scientific benefits of these observations requires a detailed understanding, based on numerical simulations, of the predictions of General Relativity for the waveform signals. New techniques for simulating binary black hole mergers, introduced two years ago, have led to dramatic advances in applied numerical simulation work. Over the last two years, numerical relativity researchers have made tremendous strides in understanding the late stages of binary black hole mergers. Simulations have been applied to test much of the basic physics of binary black hole interactions, showing robust results for merger waveform predictions, and illuminating such phenomena as spin-precession. Calculations have shown that merging systems can be kicked at up to 2500 km/s by the thrust from asymmetric emission. Recently, long lasting simulations of ten or more orbits allow tests of post-Newtonian (PN) approximation results for radiation from the last orbits of the binary's inspiral. Already, analytic waveform models based PN techniques with incorporated information from numerical simulations may be adequate for observations with current ground based observatories. As new advances in simulations continue to rapidly improve our theoretical understanding of the systems, it seems certain that high-precision predictions will be available in time for LISA and other advanced ground-based instruments. Future gravitational wave observatories are expected to make precision.
Baker, John G.↗
Laboratory measurements of the 7.5-9.38-mm absorption of gaseous ammonia (NH3) under simulated Jovian conditions
An attempt is made to infer the abundance and distribution of ammonia from RF emission measurements more accurately than heretofore, on the basis of the results of laboratory measurements for the mm-wave opacity of gaseous ammonia under simulated Jovian atmosphere conditions. The measurements were conducted at various frequencies in the 32-40 GHz range at 2 atm and 203 K; the atmospheric mixture was 88.34 percent H2, 9.81 percent He, and 1.85 percent NH3. Experimental results are found to be readily modeled by the Gross (1955) line-shape factor, rather than that of Van Vleck and Weisskopf (1945).
Joiner, Joanna↗
Gravitational waves from binary neutron star mergers with a spectral equation of state
In numerical simulations of binary neutron star systems, the equation of state of the dense neutron star matter is an important factor in determining both the physical realism and the numerical accuracy of the simulations. Some equations of state used in simulations are C 2 or smoother in the pressure/density relationship function, such as a polytropic equation of state, but may not have the flexibility to model stars or remnants of different masses while keeping their radii within known astrophysical constraints. Other equations of state, such as tabular or piece-wise polytropic, may be flexible enough to model additional physics and multiple stars' masses and radii within known constraints, but are not as smooth, resulting in additional numerical error. We will study in this paper a recently developed family of equation of state, using a spectral expansion with sufficient free parameters to allow for a larger flexibility than current polytropic equations of state, and with sufficient smoothness to reduce numerical errors compared to tabulated or piece-wise polytropic equations of state. We perform simulations at three mass ratios with a common chirp mass, using two distinct spectral equations of state, and at multiple numerical resolutions. We evaluate the gravitational waves produced from these simulations, comparing the phase error between resolutions and equations of state, as well as with respect to analytical models. From our simulations we estimate that the phase difference at merger for binaries with a dimensionless weighted tidal deformability difference greater than Δ$\tilde{Λ}$=55 can be captured by the SpEC code for these equations of state.
79 ASTRONOMY AND ASTROPHYSICS↗
Idealized Simulations of the Tropical Climate and Variability in the Single Column Atmosphere Model (SCAM): Radiative-Convective Equilibrium
To explore the interactions among column processes in the Community Atmosphere Model (CAM), the single-column version of CAM (SCAM) is integrated for 1000 days in radiative-convective equilibrium (RCE) with tropical values of boundary conditions, spanning a parameter or configuration space of model physics versions (v5 vs. v6), vertical resolution (standard and 60 levels), sea surface temperature (SST), and some interpretation-driven experiments. The simulated time-mean climate is reasonable, near observations and RCE of a cyclic cloud-resolving model. Updraft detrainment in the deep convection scheme produces distinctive grid-scale structures in humidity and cloud, which also interact with radiative transfer processes. These grid artifacts average out in multi-column RCE results reported elsewhere, illustrating the nuts-and-bolts interpretability that SCAM adds to the hierarchy of model configurations. Multi-day oscillations of precipitation arise from descent of warm convection-capping layers starting near the tropopause, eventually reset by a burst of convective deepening. Experiments reveal how these oscillations depend critically on an internal parameter that controls the number of neutral buoyancy levels allowed for determining cloud top and computing dilute convective available potential energy in the deep convection scheme, and merely modified a little by disabling cloud-base radiation (heating of cloud base). This strong dependence of transient behavior in 1D on this parameter will be tested in the second part of this work, in which SCAM is coupled to a parameterized dynamics of two-dimensional, linearized gravity wave, and in the 3D simulations in future study.
Meteorology & Atmospheric Sciences↗
Theory of Trotter Error with Commutator Scaling
The Lie-Trotter formula, together with its higher-order generalizations, provides a simple approach to decomposing the exponential of a sum of operators. Despite significant effort, the error scaling of such product formulas remains poorly understood. We develop a theory of Trotter error that overcomes the limitations of truncating the Baker-Campbell-Hausdorff expansion. Our analysis directly exploits the commutativity of operator summands, producing tighter error bounds for both real- and imaginary-time evolutions. Whereas previous work achieves similar goals for systems with geometric locality or Lie-algebraic structure, our approach holds in general. We give a host of improved algorithms for digital quantum simulation and quantum Monte Carlo methods, nearly matching or even outperforming the best previous results. Our applications include: (i) a simulation of second-quantized plane-wave electronic structure, nearly matching the interaction-picture algorithm of Low and Wiebe; (ii) a simulation of $k$-local Hamiltonians almost with induced one-norm scaling, faster than the qubitization algorithm of Low and Chuang; (iii) a simulation of rapidly decaying power-law interactions, outperforming the Lieb-Robinson-based approach of Tran et al.; (iv) a hybrid simulation of clustered Hamiltonians, dramatically improving the result of Peng, Harrow, Ozols, and Wu; and (v) quantum Monte Carlo simulations of the transverse field Ising model and quantum ferromagnets, tightening previous analyses of Bravyi and Gosset. We obtain further speedups using the fact that product formulas can preserve the locality of the simulated system. Specifically, we show that local observables can be simulated with complexity independent of the system size for power-law interacting systems, which implies a Lieb-Robinson bound nearly matching a recent result of Tran et al. Our analysis reproduces known tight bounds for first- and second-order formulas. We further investigate the tightness of our bounds for higher-order formulas. For quantum simulation of a one-dimensional Heisenberg model with an even-odd ordering of terms, our result overestimates the complexity by only a factor of $5$. Our bound is also close to tight for power-law interactions and other orderings of terms. This suggests that our theory can accurately characterize Trotter error in terms of both the asymptotic scaling and the constant prefactor.
quantum computing, numerical analysis↗
Shear Alfvén Waves Driven by Magnetic Reconnection as an Energy Source for the Aurora Borealis
Previous studies have concluded that the Hall magnetic field structures generated during magnetic reconnection are carried away by kinetic Alfvén waves. Here, we apply a kinetic simulation with an ion/electron mass ratio closer to its natural value and find that open boundary conditions and much reduced damping rates permit the shear Alfvén waves (SAWs) to become the main carrier of wave energy. Further, for magnetotail reconnection, these SAWs provide efficient transport of wave energy, enhancing the energy input for the Aurora Borealis by orders of magnitude above previous estimates.
58 GEOSCIENCES↗
Longitudinal and transverse modes of slosh wave excitation in rotating dewar associated with gravity jitters
Characteristics of slosh waves based on the dynamical behavior of oscillations at the liquid-vapor interface have been investigated. Twelve case studies of slosh wave excitation due to various frequencies of gravity jitters under different rotating speeds of the propellant tank and different levels of background gravity environment have been simulated. The study shows that slosh waves excited inside the spacecraft propellant tank are characterized by the lowest frequency of the waves initiated, frequencies of the gravity jitters imposed on the propellant system, the levels of background gravity environment, and dewar rotating speeds. Conditions for suppression and amplification of the slosh waves are discussed.
Hung, R. J.↗
Wave-current interaction study in the Gulf of Alaska for detection of eddies by synthetic aperture radar
High resolution Esa Remote Sensing Satellite-1 (ERS-1) Synthetic Aperture Radar (SAR) images are used to detect a mesoscale eddy. Such features limit dispersal of pollock larvae and therefore likely influence recruitment of fish in the Gulf of Alaska. During high sea states and high winds, the direct surface signature of the eddy was not clearly visible, but the wave refraction in the eddy area was observed. The rays of the wave field are traced out directly from the SAR image. The ray pattern gives information on the refraction pattern and on the relative variation of the wave energy along a ray through wave current interaction. These observations are simulated by a ray-tracing model which incorporates a surface current field associated with the eddy. The numerical results of the model show that the waves are refracted and diverge in the eddy field with energy density decreasing. The model-data comparison for each ray shows the model predictions are in good agreement with the SAR data.
Liu, Antony K.↗
Lumped Parameter Models of the Central Nervous System for VIIP Research
INTRODUCTION: Current long-duration missions to the International Space Station and future exploration-class missions beyond low-Earth orbit, such as to Mars and asteroids, expose astronauts to increased risk of Visual Impairment and Intracranial Pressure (VIIP) syndrome [1]. It has been hypothesized that the headward shift of cerebral spinal fluid (CSF) and blood in microgravity may cause significant elevation of intracranial pressure (ICP), which in turn induces VIIP syndrome through biomechanical pathways [1, 2]. However, there is insufficient evidence to confirm this hypothesis. In this light, we are developing lumped-parameter models of fluid transport in the central nervous system (CNS) as a means to simulate the influence of microgravity on ICP. The CNS models will also be used in concert with the lumped parameter and finite element models of the eye described in the realted IWS abstracts submitted by Nelson et al., Feola et al. and Ethier et al. METHODS: We have developed a nine compartment CNS model (Figure 1) capable of both time-dependent and steady state fluid transport simulations, based on the works of Stevens et al. [3]. The breakdown of compartments within the model includes: vascular (3), CSF (2), brain (1) and extracranial (3). The boundary pressure in the Central Arteries [A] node is prescribed using an oscillating pressure function PA(t) simulating the carotid pulsatile pressure wave as developed by Linninger et al. [4]. For each time step, pressures are integrated through time using an adaptive-timestep 4th and 5th order Runga-Kutta solver. Once pressures are found, constitutive equations are used to solve for flowrates (Q) between each compartment. In addition to fluid flow between the different compartments, compliance (C) interactions between neighboring compartments are represented. We are also developing a second CNS model based on the works of Linninger et al. [4] which takes a more granular approach to represent the interactions of the intracranial and spinal compartments with the inclusion of arteries, arterioles, capillaries, venules, veins, venous sinus, and ventricles. The flow through the arteries, veins and CSF compartments are governed by continuity, momentum and distensibility balance equations. Furthermore, unlike the Stevens et al. approach, the Monro-Kellie doctrine of constant cranial volume and the bi-phasic nature of the brain parenchyma are implemented. These features appear to be more consistent with the physiologic and anatomical behavior of the CNS, and follow a modeling philosophy similar to the lumped parameter eye model that is intended to be integrated with the CNS model. However, Linninger’s approach has never been implemented to include hydrostatic gradient and microgravity simulation capabilities. Therefore, we aim at implement this modeling approach for spaceflight simulations and assess its overall applicability to VIIP research. OBJECTIVES: We will present verification and validation test results for both models, as well as head-to-head comparison to explore their strengths and limitations with respect to mathematical implementation and physiological significance for VIIP research. In doing so, we hope to provide some guidance to the HRP research community on how to appropriately leverage lumped parameter models for space biomedical research.
Vera, J.↗
Process-oriented evaluation of quasi-stationary Rossby waves and their impact on surface air temperature extremes in dynamical downscaling over North America
Quasi-stationary Rossby waves are a crucial component of the general circulation and play a significant role in regional water and energy cycles, as well as in extreme events. However, process-oriented evaluation for Rossby waves is rarely performed for dynamical downscaling simulations. To close this gap, we evaluate three classes of dynamical downscaling approaches, with a focus on quasi-stationary Rossby waves and their impact on surface air temperature over North America during Northern Hemisphere summer. The three classes of models differ in the way large-scale forcing is provided: a limited-area model (LAM) constrained only by lateral boundary conditions, represented by RegCM4 from the North American branch of the Coordinated Regional Downscaling Experiment (NA-CORDEX), a LAM with spectral nudging to maintain consistency in large-scale dynamics with the forcing data, represented by the Weather Research and Forecasting (WRF) model simulation in NA-CORDEX, and a global variable-resolution model with smoothly varying grid spacings, represented by the Community Atmosphere Model version 5.4, with the Model for Prediction Across Scales (MPAS) as its dynamical core (CAM-MPAS). With no constraints on the atmospheric dynamics, CAM-MPAS exhibits several mean biases in the upper-level circulations over the Pacific Coast region: a weaker subtropical jet, a northward-shifted mid-latitude jet, and an overestimated southerly flow. With the lateral boundary constraint alone, RegCM4 also exhibits weaker jets and overestimated southerly winds off the West Coast. Rossby ray theory reveals that those wind biases direct incoming Rossby waves northward. The erroneously routed Rossby waves distort the relationship between the accumulation of wave activity over the US West Coast and surface temperature anomalies over the Southern Great Plains, which emerges approximately 4 d after the convergence of wave-activity flux in the ERA-Interim reanalysis. Furthermore, the response of heatwaves to the extreme wave activity flux is not reproduced by the two models, a serious drawback as a dynamical downscaling framework is expected to connect large-scale forcing to local-scale phenomena. The WRF model employing spectral nudging is largely free from the aforementioned problems. A pair of sensitivity simulations suggests that spectral nudging is the key to improving the dynamics of quasi-stationary Rossby waves and their impact on surface air temperature. Our results also demonstrate the effectiveness of Rossby wave diagnostics that allow for realistic background flows for assessing the credibility of dynamical downscaling over North America, where incoming Rossby waves propagate through complex circulation patterns before traveling across the continent.
54 ENVIRONMENTAL SCIENCES↗
Nonequilibrium, large-amplitude MHD fluctuations in the solar wind
Compressible MHD simulations in one dimension with three-dimensional vectors are used to investigate a number of processes relevant to problems in interplanetary physics. The simulations indicate that a large-amplitude nonequilibrium (e.g., linearly polarized) Alfvenic wave, which always starts with small relative fluctuations in the magnitude B of the magnetic field, typically evolves to flatten the magnetic profile in most regions. Under a wide variety of conditions B and the density rho become anticorrelated on average. If the mean magnetic field is allowed to decrease in time, the point where the transverse magnetic fluctuation amplitude delta B(sub T) is greater than the mean field B(sub 0) is not special, and large values of delta B(sub T)/B(sub 0) do not cause the compressive thermal energy to increase remarkably or the wave energy to dissipate at an unusually high rate. Nor does the 'backscatter' of the waves that occurs when the sound speed is less than the Alfven speed result, in itself, in substantial energy dissipation, but rather primarily in a phase change between the magnetic and velocity fields. For isolated wave packets the backscatter does not occur for any of the parameters examined; an initial radiation of acoustic waves away from the packet establishes a stable traveling structure. Thus these simulations, although greatly idealized compared to reality, suggest a picture in which the interplanetary fluctuations should have small deltaB and increasingly quasi-pressure balanced compressive fluctuations, as observed, and in which the dissipation and 'saturation' at delta B(sub T)/B(sub 0) approximately = 1 required by some theories of wave acceleration of the solar wind do not occur. The simulations also provide simple ways to understand the processes of nonlinear steepening and backscattering of Alfven waves and demonstrate the existence of previously unreported types of quasi-steady MHD states.
Roberts, D. Aaron↗
Global Distribution of EMIC Waves and Its Association to Subauroral Proton Precipitation During the 27 May 2017 Storm: Modeling and Multipoint Observations
Recent simulation studies using the RAM-SCB model showed that proton precipitation contributes significantly to the total energy flux deposited into the subauroral ionosphere thereby affecting the magnetosphere-ionosphere coupling. Here, in this study, we use the BATS-R-US + RAM-SCB model to understand the evolution of ElectroMagnetic Ion Cyclotron (EMIC) waves in the inner magnetosphere, their correspondence to the proton precipitation into the subauroral ionosphere, and to assess the performance of the model in reproducing the EMIC wave-particle interactions. During the 27 May 2017 storm, Arase and RBSP-A satellites observed typical signatures of EMIC waves in the inner magnetosphere. Within this interval, Defense Meteorological Satellite Program (DMSP) and National Oceanic and Atmospheric Administration (NOAA)/MetOp satellites observed significant proton precipitation in the dusk-midnight sector. Simulation results show that H- and He-band EMIC waves are excited within regions of strong temperature anisotropy near the plasmapause. The simulated growth rates of EMIC waves show a similar trend to that of the EMIC wave power observed by the Arase and RBSP-A satellites, suggesting that the model can reproduce the EMIC wave activity qualitatively. The simulated H-band waves in the dusk sector are stronger than He-band waves possibly due to the presence of excess protons in the boundary conditions obtained from the BATS-R-US code. The precipitating proton fluxes reproduced by the simulation with EMIC waves are found to agree reasonably well with the DMSP and NOAA/MetOp satellite observations. It is suggested that EMIC wave scattering of ring current ions can account for proton precipitation observed by the DMSP and MetOp satellites during the 27 May 2017 storm.
79 ASTRONOMY AND ASTROPHYSICS↗