An alternative to reduction of surface pressure to sea level

The pitfalls of the present method of reducing surface pressure to sea level are reviewed, and an alternative, adjusted pressure, P, is proposed. P is obtained from solution of a Poisson equation over a continental region, using the simplest boundary condition along the perimeter or coastline where P equals the sea level pressure. The use of P would avoid the empiricisms and disadvantages of pressure reduction to sea level, and would produce surface pressure charts which depict the true geostrophic wind at the surface

Vertical distribution of wind speed, temperature and humidity above a water surface

An observational program is described which has been used to obtain an accurate determination of vertical profiles of wind speed, temperature, and vapor pressure over a salt water inlet with an over-water air fetch of about five miles. The wind profiles show systematic anomalies of 1 or 2 %, which are not explainable as instrumental or observational error. The curvature of wind profiles over water shows the same dependence on Richardson number as that found by others over land. Temperature profiles are similar in this respect, but curvature of the vapor pressure profiles shows little dependence on stability. Values of the resistance coefficient computed from wind profiles are at the lower limit of those reported by most other investigators...

Matrix exponential-based closures for the turbulent subgrid-scale stress tensor

Two approaches for closing the turbulence subgrid-scale stress tensor in terms of matrix exponentials are introduced and compared. The first approach is based on a formal solution of the stress transport equation in which the production terms can be integrated exactly in terms of matrix exponentials. This formal solution of the subgrid-scale stress transport equation is shown to be useful to explore special cases, such as the response to constant velocity gradient, but neglecting pressure-strain correlations and diffusion effects. The second approach is based on an Eulerian-Lagrangian change of variables, combined with the assumption of isotropy for the conditionally averaged Lagrangian velocity gradient tensor and with the recent fluid deformation approximation. It is shown that both approaches lead to the same basic closure in which the stress tensor is expressed as the matrix exponential of the resolved velocity gradient tensor multiplied by its transpose. Short-time expansions of the matrix exponentials are shown to provide an eddy-viscosity term and particular quadratic terms, and thus allow a reinterpretation of traditional eddy-viscosity and nonlinear stress closures. The basic feasibility of the matrix-exponential closure is illustrated by implementing it successfully in large eddy simulation of forced isotropic turbulence. The matrix-exponential closure employs the drastic approximation of entirely omitting the pressure-strain correlation and other nonlinear scrambling terms. But unlike eddy-viscosity closures, the matrix exponential approach provides a simple and local closure that can be derived directly from the stress transport equation with the production term, and using physically motivated assumptions about Lagrangian decorrelation and upstream isotropy

The energy budget in Rayleigh-Benard convection

It is shown using three series of Rayleigh number simulations of varying aspect ratio AR and Prandtl number Pr that the normalized dissipation at the wall, while significantly greater than 1, approaches a constant dependent upon AR and Pr. It is also found that the peak velocity, not the mean square velocity, obeys the experimental scaling of Ra^{0.5}. The scaling of the mean square velocity is closer to Ra^{0.46}, which is shown to be consistent with experimental measurements and the numerical results for the scaling of Nu and the temperature if there are strong correlations between the velocity and temperature.Comment: 5 pages, 3 figures, new version 13 Mar, 200

Adaptively refined large eddy simulations of clusters

We present a numerical scheme for modelling unresolved turbulence in cosmological adaptive mesh refinement codes. As a first application, we study the evolution of turbulence in the intra-cluster medium and in the core of a galaxy cluster. Simulations with and without subgrid scale model are compared in detail. Since the flow in the ICM is subsonic, the global turbulent energy contribution at the unresolved length scales is smaller than 1% of the internal energy. We find that the production of turbulence is closely correlated with merger events occurring in the cluster environment, and its dissipation locally affects the cluster energy budget. Because of this additional source of dissipation, the core temperature is larger and the density is smaller in the presence of subgrid scale turbulence than in the standard adiabatic run, resulting in a higher entropy core value.Comment: Submitted to ApJ, 14 pages, 14 figures, 3 table

Analysis of the shearing instability in nonlinear convection and magnetoconvection

Numerical experiments on two-dimensional convection with or without a vertical magnetic field reveal a bewildering variety of periodic and aperiodic oscillations. Steady rolls can develop a shearing instability, in which rolls turning over in one direction grow at the expense of rolls turning over in the other, resulting in a net shear across the layer. As the temperature difference across the fluid is increased, two-dimensional pulsating waves occur, in which the direction of shear alternates. We analyse the nonlinear dynamics of this behaviour by first constructing appropriate low-order sets of ordinary differential equations, which show the same behaviour, and then analysing the global bifurcations that lead to these oscillations by constructing one-dimensional return maps. We compare the behaviour of the partial differential equations, the models and the maps in systematic two-parameter studies of both the magnetic and the non-magnetic cases, emphasising how the symmetries of periodic solutions change as a result of global bifurcations. Much of the interesting behaviour is associated with a discontinuous change in the leading direction of a fixed point at a global bifurcation; this change occurs when the magnetic field is introduced

GPU-Accelerated Large-Eddy Simulation of Turbulent Channel Flows

High performance computing clusters that are augmented with cost and power efficient graphics processing unit (GPU) provide new opportunities to broaden the use of large-eddy simulation technique to study high Reynolds number turbulent flows in fluids engineering applications. In this paper, we extend our earlier work on multi-GPU acceleration of an incompressible Navier-Stokes solver to include a large-eddy simulation (LES) capability. In particular, we implement the Lagrangian dynamic subgrid scale model and compare our results against existing direct numerical simulation (DNS) data of a turbulent channel flow at Reτ = 180. Overall, our LES results match fairly well with the DNS data. Our results show that the Reτ = 180 case can be entirely simulated on a single GPU, whereas higher Reynolds cases can benefit from a GPU cluster

Spiral valve parasites of blue and common thresher sharks as indicators of shark feeding behaviour and ecology

Open Access via the Jisc Wiley agreement Acknowledgements This work would not have been possible without the assistance and samples provided by the National Marine Fisheries Service (NMFS) Southwest Region Fishery Observer Program and the participating drift gillnet fishermen. A. Arevalo, E. Reed, H. Colley, J. Williams, J. Tamez and K. Tran assisted with spiral valve dissections and parasite sorting in the lab. D. Losey helped with library research. D. Sweetnam, A. Yau, A. Thompson, M. Craig, S. Stohs, G. DiNardo provided constructive critiques that helped improve the manuscript. This research was supported by the National Oceanographic Atmospheric Administration (NOAA).Peer reviewedPublisher PD

Stereodivergent Synthesis of Enantioenriched 4-Hydroxy-2- cyclopentenones

Protected 4-hydroxycyclopentenones (4-HCPs) constitute an important class of intermediates in chemical synthesis. A route to this class of compound has been developed. Key steps include Noyori reduction (which establishes the stereochemistry of the product), ring-closing metathesis, and simple functional group conversions to provide a set of substituted 4-HCPs in either enantiomeric form