Matching Items (21)

132633-Thumbnail Image.png

A High Order Numerical Scheme for Wave Propagation Problems

Description

The physics of waves control most of the world, in multiple forms, such as electromagnetic waves. Mathematicians and physicists have developed equations which describe the patterns in which waves evolve over time, while moving through space. Due to their partial

The physics of waves control most of the world, in multiple forms, such as electromagnetic waves. Mathematicians and physicists have developed equations which describe the patterns in which waves evolve over time, while moving through space. Due to their partial differential form, solutions to these equations must be approximated. This study introduces a new numerical scheme to perform the approximation which is highly stable and computationally efficient. This numerical scheme is formulated with respect to Maxwell’s equations, employing spatial and temporal staggering to implement a fourth-order phase accuracy. It is then compared to the traditional Yee scheme and the Runge-Kutta 3 scheme in one-dimensional applications, revealing a similar accuracy to the Runge-Kutta 3 scheme while requiring less computations per time step. Simulations are then performed in the two-dimensional case. First, no boundary conditions are implemented, causing reflection at the edge of the spatial domain. Next, the simulation is conducted while employing absorbing boundary conditions, simulating wave propagation over an infinite spatial domain. These results are compared to the results of a large domain simulation, in which the wave propagation does not reach the boundaries. Comparing the simulations, it is concluded that the numerical scheme is stable and highly accurate when employing absorbing boundary conditions. Finally, the scheme is tested in two dimensions with wave propagation through nonlinear media, as opposed to the prior simulations which were performed as if in a vacuum. After performing spectral analysis on the resulting waves after a long-time domain simulation, the resulting angular frequencies match those expected from theory. Therefore, the scheme is concluded to be powerful in one-dimensional, two-dimensional, and nonlinear simulations, all while being computationally efficient.

Contributors

Agent

Created

Date Created
2019-05

134565-Thumbnail Image.png

A Numerical and Analytical Study of Wave Reflection and Transmission across the Tropopause

Description

A numerical study of wave-induced momentum transport across the tropopause in the presence of a stably stratified thin inversion layer is presented and discussed. This layer consists of a sharp increase in static stability within the tropopause. The wave propagation

A numerical study of wave-induced momentum transport across the tropopause in the presence of a stably stratified thin inversion layer is presented and discussed. This layer consists of a sharp increase in static stability within the tropopause. The wave propagation is modeled by numerically solving the Taylor-Goldstein equation, which governs the dynamics of internal waves in stably stratified shear flows. The waves are forced by a flow over a bell shaped mountain placed at the lower boundary of the domain. A perfectly radiating condition based on the group velocity of mountain waves is imposed at the top to avoid artificial wave reflection. A validation for the numerical method through comparisons with the corresponding analytical solutions will be provided. Then, the method is applied to more realistic profiles of the stability to study the impact of these profiles on wave propagation through the tropopause.

Contributors

Created

Date Created
2017-05

129251-Thumbnail Image.png

Simulating Meteorological Profiles to Study Noise Propagation From Freeways

Description

Forecasts of noise pollution from a highway line segment noise source are obtained from a sound propagation model utilizing effective sound speed profiles derived from a Numerical Weather Prediction (NWP) limited area forecast with 1 km horizontal resolution and near-ground

Forecasts of noise pollution from a highway line segment noise source are obtained from a sound propagation model utilizing effective sound speed profiles derived from a Numerical Weather Prediction (NWP) limited area forecast with 1 km horizontal resolution and near-ground vertical resolution finer than 20 m. Methods for temporal along with horizontal and vertical spatial nesting are demonstrated within the NWP model for maintaining forecast feasibility. It is shown that vertical nesting can improve the prediction of finer structures in near-ground temperature and velocity profiles, such as morning temperature inversions and low level jet-like features. Accurate representation of these features is shown to be important for modeling sound refraction phenomena and for enabling accurate noise assessment. Comparisons are made using the parabolic equation model for predictions with profiles derived from NWP simulations and from field experiment observations during mornings on November 7 and 8, 2006 in Phoenix, Arizona. The challenges faced in simulating accurate meteorological profiles at high resolution for sound propagation applications are highlighted and areas for possible improvement are discussed.

Contributors

Created

Date Created
2015-05-01

129252-Thumbnail Image.png

Multiscale Modeling and Evaluation of Urban Surface Energy Balance in the Phoenix Metropolitan Area

Description

Physical mechanisms of incongruency between observations and Weather Research and Forecasting (WRF) Model predictions are examined. Limitations of evaluation are constrained by (i) parameterizations of model physics, (ii) parameterizations of input data, (iii) model resolution, and (iv) flux observation resolution.

Physical mechanisms of incongruency between observations and Weather Research and Forecasting (WRF) Model predictions are examined. Limitations of evaluation are constrained by (i) parameterizations of model physics, (ii) parameterizations of input data, (iii) model resolution, and (iv) flux observation resolution. Observations from a new 22.1-m flux tower situated within a residential neighborhood in Phoenix, Arizona, are utilized to evaluate the ability of the urbanized WRF to resolve finescale surface energy balance (SEB) when using the urban classes derived from the 30-m-resolution National Land Cover Database. Modeled SEB response to a large seasonal variation of net radiation forcing was tested during synoptically quiescent periods of high pressure in winter 2011 and premonsoon summer 2012. Results are presented from simulations employing five nested domains down to 333-m horizontal resolution. A comparative analysis of model cases testing parameterization of physical processes was done using four configurations of urban parameterization for the bulk urban scheme versus three representations with the Urban Canopy Model (UCM) scheme, and also for two types of planetary boundary layer parameterization: the local Mellor–Yamada–Janjić scheme and the nonlocal Yonsei University scheme. Diurnal variation in SEB constituent fluxes is examined in relation to surface-layer stability and modeled diagnostic variables. Improvement is found when adapting UCM for Phoenix with reduced errors in the SEB components. Finer model resolution is seen to have insignificant (<1 standard deviation) influence on mean absolute percent difference of 30-min diurnal mean SEB terms.

Contributors

Created

Date Created
2015-06-11

Achieving Accurate Simulations of Urban Impacts on Ozone at High Resolution

Description

The effects of urbanization on ozone levels have been widely investigated over cities primarily located in temperate and/or humid regions. In this study, nested WRF-Chem simulations with a finest grid resolution of 1 km are conducted to investigate ozone concentrations

The effects of urbanization on ozone levels have been widely investigated over cities primarily located in temperate and/or humid regions. In this study, nested WRF-Chem simulations with a finest grid resolution of 1 km are conducted to investigate ozone concentrations O3 due to urbanization within cities in arid/semi-arid environments. First, a method based on a shape preserving Monotonic Cubic Interpolation (MCI) is developed and used to downscale anthropogenic emissions from the 4 km resolution 2005 National Emissions Inventory (NEI05) to the finest model resolution of 1 km. Using the rapidly expanding Phoenix metropolitan region as the area of focus, we demonstrate the proposed MCI method achieves ozone simulation results with appreciably improved correspondence to observations relative to the default interpolation method of the WRF-Chem system. Next, two additional sets of experiments are conducted, with the recommended MCI approach, to examine impacts of urbanization on ozone production: (1) the urban land cover is included (i.e., urbanization experiments) and, (2) the urban land cover is replaced with the region's native shrubland. Impacts due to the presence of the built environment on O3 are highly heterogeneous across the metropolitan area. Increased near surface O3 due to urbanization of 10–20 ppb is predominantly a nighttime phenomenon while simulated impacts during daytime are negligible. Urbanization narrows the daily O3 range (by virtue of increasing nighttime minima), an impact largely due to the region's urban heat island. Our results demonstrate the importance of the MCI method for accurate representation of the diurnal profile of ozone, and highlight its utility for high-resolution air quality simulations for urban areas.

Contributors

Created

Date Created
2014-11-01

128203-Thumbnail Image.png

Regional-Scale Transport of Air Pollutants: Impacts of Southern California Emissions on Phoenix Ground-Level Ozone Concentrations

Description

In this study, WRF-Chem is utilized at high resolution (1.333 km grid spacing for the innermost domain) to investigate impacts of southern California anthropogenic emissions (SoCal) on Phoenix ground-level ozone concentrations ([O3]) for a pair of recent exceedance episodes. First,

In this study, WRF-Chem is utilized at high resolution (1.333 km grid spacing for the innermost domain) to investigate impacts of southern California anthropogenic emissions (SoCal) on Phoenix ground-level ozone concentrations ([O3]) for a pair of recent exceedance episodes. First, WRF-Chem control simulations, based on the US Environmental Protection Agency (EPA) 2005 National Emissions Inventories (NEI05), are conducted to evaluate model performance. Compared with surface observations of hourly ozone, CO, NOX, and wind fields, the control simulations reproduce observed variability well. Simulated [O3] are comparable with the previous studies in this region. Next, the relative contribution of SoCal and Arizona local anthropogenic emissions (AZ) to ozone exceedances within the Phoenix metropolitan area is investigated via a trio of sensitivity simulations: (1) SoCal emissions are excluded, with all other emissions as in Control; (2) AZ emissions are excluded with all other emissions as in Control; and (3) SoCal and AZ emissions are excluded (i.e., all anthropogenic emissions are eliminated) to account only for Biogenic emissions and lateral boundary inflow (BILB). Based on the USEPA NEI05, results for the selected events indicate the impacts of AZ emissions are dominant on daily maximum 8 h average (DMA8) [O3] in Phoenix. SoCal contributions to DMA8 [O3] for the Phoenix metropolitan area range from a few ppbv to over 30 ppbv (10–30 % relative to Control experiments). [O3] from SoCal and AZ emissions exhibit the expected diurnal characteristics that are determined by physical and photochemical processes, while BILB contributions to DMA8 [O3] in Phoenix also play a key role.

Contributors

Created

Date Created
2015-08-21

127861-Thumbnail Image.png

Observation and Simulation of Wave Breaking in the Southern Hemispheric Stratosphere During VORCORE

Description

An interesting occurrence of a Rossby wave breaking event observed during the VORCORE experiment is presented and explained. Twenty-seven balloons were launched inside the Antarctic polar vortex. Almost all of these balloons evolved in the stratosphere around 500K within the

An interesting occurrence of a Rossby wave breaking event observed during the VORCORE experiment is presented and explained. Twenty-seven balloons were launched inside the Antarctic polar vortex. Almost all of these balloons evolved in the stratosphere around 500K within the vortex, except the one launched on 28 October 2005. In this case, the balloon was caught within a tongue of high potential vorticity (PV), and was ejected from the polar vortex. The evolution of this event is studied for the period between 19 and 25 November 2005. It is found that at the beginning of this period, the polar vortex experienced distortions due to the presence of Rossby waves. Then, these waves break and a tongue of high PV develops. On 25 November, the tongue became separated from the vortex and the balloon was ejected into the surf zone. Lagrangian simulations demonstrate that the air masses surrounding the balloon after its ejection were originating from the vortex edge. The wave breaking and the development of the tongue are confined within a region where a planetary Quasi-Stationary Wave 1 (QSW1) induces wind speeds with weaker values. The QSW1 causes asymmetry in the wind speed and the horizontal PV gradient along the edge of the polar vortex, resulting in a localized jet. Rossby waves with smaller scales propagating on top of this jet amplify as they enter the jet exit region and then break. The role of the QSW1 on the formation of the weak flow conditions that caused the non-linear wave breaking observed near the vortex edge is confirmed by three-dimensional numerical simulations using forcing with and without the contribution of the QSW1.

Contributors

Created

Date Created
2013-04-16

135327-Thumbnail Image.png

A New Numerical Method Based on Leapfrog for Atmospheric and Ocean Modeling

Description

A semi-implicit, fourth-order time-filtered leapfrog numerical scheme is investigated for accuracy and stability, and applied to several test cases, including one-dimensional advection and diffusion, the anelastic equations to simulate the Kelvin-Helmholtz instability, and the global shallow water spectral model to

A semi-implicit, fourth-order time-filtered leapfrog numerical scheme is investigated for accuracy and stability, and applied to several test cases, including one-dimensional advection and diffusion, the anelastic equations to simulate the Kelvin-Helmholtz instability, and the global shallow water spectral model to simulate the nonlinear evolution of twin tropical cyclones. The leapfrog scheme leads to computational modes in the solutions to highly nonlinear systems, and time-filters are often used to damp these modes. The proposed filter damps the computational modes without appreciably degrading the physical mode. Its performance in these metrics is superior to the second-order time-filtered leapfrog scheme developed by Robert and Asselin.

Contributors

Created

Date Created
2016-05

136845-Thumbnail Image.png

Visualizing Lagrangian Coherent Structures from Realistic Flow Data

Description

The goal of this project was to examine the separatricies that define regions of distinct flow behaviors in realistic time-dependent dynamical systems. In particular, we adapted previously available methods for computing the Finite-Time Lyapunov Exponent (FTLE) to a set of

The goal of this project was to examine the separatricies that define regions of distinct flow behaviors in realistic time-dependent dynamical systems. In particular, we adapted previously available methods for computing the Finite-Time Lyapunov Exponent (FTLE) to a set of measured wind velocity data in order to visualize the separatricies as ridges of the FTLE field in a section of the atmosphere. This visualization required a number of alterations to the original methods, including interpolation techniques and two different adaptive refinement schemes for producing more detailed results. Overall, there were two computations performed with the wind velocity data: once along a single spherical surface, on which the separatricies could be visualized as material lines, and then along a three-dimensional section of the atmosphere, for which the separatricies were material surfaces. The resulting figures provide an image of the Antarctic polar vortex from the wind velocity data, which is consistent with other data gathered on the same date.

Contributors

Created

Date Created
2014-05

137818-Thumbnail Image.png

Lagrangian Skeletons in Hurricane Katrina

Description

This thesis shows analyses of mixing and transport patterns associated with Hurricane Katrina as it hit the United States in August of 2005. Specifically, by applying atmospheric velocity information from the Weather Research and Forecasting System, finite-time Lyapunov exponents have

This thesis shows analyses of mixing and transport patterns associated with Hurricane Katrina as it hit the United States in August of 2005. Specifically, by applying atmospheric velocity information from the Weather Research and Forecasting System, finite-time Lyapunov exponents have been computed and the Lagrangian Coherent Structures have been identified. The chaotic dynamics of material transport induced by the hurricane are results from these structures within the flow. Boundaries of the coherent structures are highlighted by the FTLE field. Individual particle transport within the hurricane is affected by the location of these boundaries. In addition to idealized fluid particles, we also studied inertial particles which have finite size and inertia. Basing on established Maxey-Riley equations of the dynamics of particles of finite size, we obtain a reduced equation governing the position process. Using methods derived from computer graphics, we identify maximizers of the FTLE field. Following and applying these ideas, we analyze the dynamics of inertial particle transport within Hurricane Katrina, through comparison of trajectories of dierent sized particles and by pinpointing the location of the Lagrangian Coherent Structures.

Contributors

Agent

Created

Date Created
2012-12