ALBERT

Description: Accurate retrievals of aerosol size distribution are necessary to estimate aerosols' impact on climate and human health. The inversions of the Aerosol Robotic Network (AERONET) usually retrieve bimodal distributions. However, when the inversion is applied to Saharan and Sahelian dust, an additional mode of intermediate size between the coarse and fine modes is sometimes seen. This mode explains peculiarities in the behavior of the Angstrom exponent, along with the fine mode fraction retrieved using the spectral deconvolution algorithm, observed in a March 2006 dust storm. For this study, 15 AERONET sites in northern Africa and on the Atlantic are examined to determine the frequency and properties of the intermediate mode. The mode is observed most frequently at Ilorin in Nigeria. It is also observed at Capo Verde and multiple sites located within the Sahel but much less frequently at sites in the northern Sahara and the Canary Islands. The presence of the intermediate mode coincides with increases in Angstrom exponent, fine mode fraction, single-scattering albedo, and to a lesser extent percent sphericity. The Angstrom exponent decreases with increasing optical depth at most sites when the intermediate mode is present, but the fine mode fraction does not. Single-scattering albedo does not steadily decrease with fine mode fraction when the intermediate mode is present, as it does in typical mixtures of dust and biomass-burning aerosols. Continued investigation is needed to further define the intermediate mode's properties, determine why it differs from most Saharan dust, and identify its climate and health effects.

Description: This work uses a discontinuous-Galerkin spectral-element method (DGSEM) to solve the compressible Navier-Stokes equations [1{3]. The inviscid ux is computed using the approximate Riemann solver of Roe [4]. The viscous fluxes are computed using the second form of Bassi and Rebay (BR2) [5] in a manner consistent with the spectral-element approximation. The method of lines with the classical 4th-order explicit Runge-Kutta scheme is used for time integration. Results for polynomial orders up to p = 15 (16th order) are presented. The code is parallelized using the Message Passing Interface (MPI). The computations presented in this work are performed using the Sandy Bridge nodes of the NASA Pleiades supercomputer at NASA Ames Research Center. Each Sandy Bridge node consists of 2 eight-core Intel Xeon E5-2670 processors with a clock speed of 2.6Ghz and 2GB per core memory. On a Sandy Bridge node the Tau Benchmark [6] runs in a time of 7.6s.

Description: Forced periodic flows arise in a broad range of aerodynamic applications such as rotorcraft, turbomachinery, and flapping wing configurations. Standard practice involves solving the unsteady flow equations forward in time until the initial transient exits the domain and a statistically stationary flow is achieved. It is often required to simulate through several periods to remove the initial transient making unsteady design optimization prohibitively expensive for most realistic problems. An effort to reduce the computational cost of these calculations led to the development of the Harmonic Balance method [1, 2] which capitalizes on the periodic nature of the solution. The approach exploits the fact that forced temporally periodic flow, while varying in the time domain, is invariant in the frequency domain. Expanding the temporal variation at each spatial node into a Fourier series transforms the unsteady governing equations into a steady set of equations in integer harmonics that can be tackled with the acceleration techniques afforded to steady-state flow solvers. Other similar approaches, such as the Nonlinear Frequency Domain [3,4,5], Reduced Frequency [6] and Time-Spectral [7, 8, 9] methods, were developed shortly thereafter. Additionally, adjoint-based optimization techniques can be applied [10, 11] as well as frequency-adaptive methods [12, 13, 14] to provide even more flexibility to the method. The Fourier temporal basis functions imply spectral convergence as the number of harmonic modes, and correspondingly number of time samples, N, is increased. Some elect to solve the equations in the frequency domain directly, while others choose to transform the equations back into the time domain to simplify the process of adding this capability to existing solvers, but each harnesses the underlying steady solution in the frequency domain. These temporal projection methods will herein be collectively referred to as Time-Spectral methods. Time-Spectral methods have demonstrated marked success in reducing the computational costs associated with simulating periodic forced flows, but have yet to be fully applied to overset or Cartesian solvers for arbitrary motion with dynamic hole-cutting. Overset and Cartesian grid methodologies are versatile techniques capable of handling complex geometry configurations in practical engineering applications, and the combination of the Time-Spectral approach with this general capability potentially provides an enabling new design and analysis tool. In an arbitrary moving-body scenario for these approaches, a Lagrangian body moves through a fixed Eulerian mesh and mesh points in the Eulerian mesh interior to the solid body are removed (cut or blanked), leaving a hole in the Eulerian mesh. During the dynamic motion some gridpoints in the domain are blanked and do not have a complete set of time-samples preventing a direct implementation of the Time-Spectral method. Murman[6] demonstrated the Time-Spectral approach for a Cartesian solver with a rigid domain motion, wherein the hole cutting remains constant. Similarly, Custer et al. [15, 16] used the NASA overset OVERFLOW solver and limited the amount of relative motion to ensure static hole-cutting and interpolation. Recently, Mavriplis and Mundis[17] demonstrated a qualitative method for applying the Time-Spectral approach to an unstructured overset solver for arbitrary motion. The goal of the current work is to develop a robust and general method for handling arbitrary motion with the Time-Spectral approach within an overset or Cartesian mesh method, while still approaching the spectral convergence rate of the original Time-Spectral approach. The viscous OVERFLOW solver will be augmented with the new Time-Spectral algorithm and the capability of the method for benchmark problems in rotorcraft and turbomachinery will be demonstrated. This abstract begins with a brief synopsis of the Time-Spectral approach for overset grids and provides details of e current approach to allow for arbitrary motion. Model problem results in one and two dimensions are included to demonstrate the viability of the method and the convergence properties. Section IV briefly outlines the implementation into the OVERFLOW solver, and the abstract closes with a description of the benchmark test cases which will be included in the final paper.

Description: Management of water resources is a balancing act between temporally and spatially limited sources and competitive needs which can often exceed the supply. In order to manage water resources over a region such as the San Joaquin Valley or the Tennessee River Valley, it is pertinent to know the amount of water that has fallen in the watershed and where the water is going within it. Since rain gauge networks are typically sparsely spaced, it is typical that the majority of rainfall on the region may not be measured. To mitigate this under-sampling of rainfall, weather radar has long been employed to provide areal rainfall estimates. The Next-Generation Weather Radars (NEXRAD) make it possible to estimate rainfall over the majority of the conterminous United States. The NEXRAD Rainfall Estimation Processing System (NREPS) was developed specifically for the purpose of using weather radar to estimate rainfall for water resources management. The NREPS is tailored to meet customer needs on spatial and temporal scales relevant to the hydrologic or land-surface models of the end-user. It utilizes several techniques to mitigate artifacts in the NEXRAD data from contaminating the rainfall field. These techniques include clutter filtering, correction for occultation by topography as well as accounting for the vertical profile of reflectivity. This presentation will focus on improvements made to the NREPS system to map rainfall in the San Joaquin Valley for NASA s Water Supply and Management Project in California, but also ongoing rainfall mapping work in the Tennessee River watershed for the Tennessee Valley Authority and possible future applications in other areas of the continent.

Description: Direct numerical simulation (DNS) of turbulent compressible flows is performed using a higher-order space-time discontinuous-Galerkin finite-element method. The numerical scheme is validated by performing DNS of the evolution of the Taylor-Green vortex and turbulent flow in a channel. The higher-order method is shown to provide increased accuracy relative to low-order methods at a given number of degrees of freedom. The turbulent flow over a periodic array of hills in a channel is simulated at Reynolds number 10,595 using an 8th-order scheme in space and a 4th-order scheme in time. These results are validated against previous large eddy simulation (LES) results. A preliminary analysis provides insight into how these detailed simulations can be used to improve Reynoldsaveraged Navier-Stokes (RANS) modeling