A Multiresolution Census Algorithm for Calculating Vortex Statistics in Turbulent Flows

The fundamental equations that model turbulent flow do not provide much insight into the size and shape of observed turbulent structures. We investigate the efficient and accurate representation of structures in two-dimensional turbulence by applying statistical models directly to the simulated vorticity field. Rather than extract the coherent portion of the image from the background variation, as in the classical signal-plus-noise model, we present a model for individual vortices using the non-decimated discrete wavelet transform. A template image, supplied by the user, provides the features to be extracted from the vorticity field. By transforming the vortex template into the wavelet domain, specific characteristics present in the template, such as size and symmetry, are broken down into components associated with spatial frequencies. Multivariate multiple linear regression is used to fit the vortex template to the vorticity field in the wavelet domain. Since all levels of the template decomposition may be used to model each level in the field decomposition, the resulting model need not be identical to the template. Application to a vortex census algorithm that records quantities of interest (such as size, peak amplitude, circulation, etc.) as the vorticity field evolves is given. The multiresolution census algorithm extracts coherent structures of all shapes and sizes in simulated vorticity fields and is able to reproduce known physical scaling laws when processing a set of voriticity fields that evolve over time

Properties of continuous Fourier extension of the discrete cosine transform and its multidimensional generalization

A versatile method is described for the practical computation of the discrete Fourier transforms (DFT) of a continuous function $g(t)$ given by its values $g_{j}$ at the points of a uniform grid $F_{N}$ generated by conjugacy classes of elements of finite adjoint order $N$ in the fundamental region $F$ of compact semisimple Lie groups. The present implementation of the method is for the groups SU(2), when $F$ is reduced to a one-dimensional segment, and for $SU(2)\times ... \times SU(2)$ in multidimensional cases. This simplest case turns out to result in a transform known as discrete cosine transform (DCT), which is often considered to be simply a specific type of the standard DFT. Here we show that the DCT is very different from the standard DFT when the properties of the continuous extensions of these two discrete transforms from the discrete grid points $t_j; j=0,1, ... N$ to all points $t \in F$ are considered. (A) Unlike the continuous extension of the DFT, the continuous extension of (the inverse) DCT, called CEDCT, closely approximates $g(t)$ between the grid points $t_j$. (B) For increasing $N$, the derivative of CEDCT converges to the derivative of $g(t)$. And (C), for CEDCT the principle of locality is valid. Finally, we use the continuous extension of 2-dimensional DCT to illustrate its potential for interpolation, as well as for the data compression of 2D images.Comment: submitted to JMP on April 3, 2003; still waiting for the referee's Repor

Gamow-Teller strength distributions for nuclei in pre-supernova stellar cores

Electron-capture and $\beta$-decay of nuclei in the core of massive stars play an important role in the stages leading to a type II supernova explosion. Nuclei in the f-p shell are particularly important for these reactions in the post Silicon-burning stage of a presupernova star. In this paper, we characterise the energy distribution of the Gamow-Teller Giant Resonance (GTGR) for mid-fp-shell nuclei in terms of a few shape parameters, using data obtained from high energy, forward scattering (p,n) and (n,p) reactions. The energy of the GTGR centroid $E_{GT}$ is further generalised as function of nuclear properties like mass number, isospin and other shell model properties of the nucleus. Since a large fraction of the GT strength lies in the GTGR region, and the GTGR is accessible for weak transitions taking place at energies relevant to the cores of presupernova and collapsing stars, our results are relevant to the study of important $e^-$-capture and $\beta$-decay rates of arbitrary, neutron-rich, f-p shell nuclei in stellar cores. Using the observed GTGR and Isobaric Analog States (IAS) energy systematics we compare the coupling coefficients in the Bohr-Mottelson two particle interaction Hamiltonian for different regions of the Isotope Table.Comment: Revtex, 28 pages +7 figures (PostScript Figures, uuencoded, filename: Sutfigs.uu). If you have difficulty printing the figures, please contact [email protected]. Accepted for publication in Phys. Rev. C, Nov 01, 199

Progressive Meshes in an Operational Rate-Distortion Sense with Application to Terrain Data

This paper presents an efficient simplification method for regular meshes obtained with a binary subdivision scheme. Our mesh connectivity is constrained with a quadtree data structure. We propose a quadtree built especially for this class of meshes having a constant-time traversal property. We introduce a rate-distortion (RD) framework to decimate the mesh and build a progressive representation for the model. We propose to achieve the RD-optimal solutions for our quadtree-restricted setting: to obtain the optimal solutions, we show how to find the vertex in the quadtree yielding the best RD trade-off and then perform optimizations at variable rate, where the rate is given by a cost function (for example the number of triangles). All previous methods are restricted to constant rate optimization only. We compare the optimal approach to its greedy counterpart. We give computationally optimal formulations for all our algorithms on the quadtree. We apply our technique to a large dataset of terrains and give extensive experimental results

Estimation of urban sensible heat flux using a dense wireless network of observations

The determination of the sensible heat flux over urban terrain is challenging due to irregular surface geometry and surface types. To address this, in 2006-07, a major field campaign (LUCE) took place at the École Polytechnique Fédérale de Lausanne campus, a moderately occupied urban site. A distributed network of 92 wireless weather stations was combined with routine atmospheric profiling, offering high temporal and spatial resolution meteorological measurements. The objective of this study is to estimate the sensible heat flux over the built environment under convective conditions. Calculations were based on Monin-Obukhov similarity for temperature in the surface layer. The results illustrate a good agreement between the sensible heat flux inferred from the thermal roughness length approach and independent calibrated measurements from a scintillometer located inside the urban canopy. It also shows that using only one well-selected station can provide a good estimate of the sensible heat flux over the campus for convective conditions. Overall, this study illustrates how an extensive network of meteorological measurements can be a useful tool to estimate the sensible heat flux in complex urban environment

Practical solution to the Monte Carlo sign problem: Realistic calculations of 54Fe

We present a practical solution to the "sign problem" in the auxiliary field Monte Carlo approach to the nuclear shell model. The method is based on extrapolation from a continuous family of problem-free Hamiltonians. To demonstrate the resultant ability to treat large shell-model problems, we present results for 54Fe in the full fp-shell basis using the Brown-Richter interaction. We find the Gamow-Teller beta^+ strength to be quenched by 58% relative to the single-particle estimate, in better agreement with experiment than previous estimates based on truncated bases.Comment: 11 pages + 2 figures (not included

Electron capture on iron group nuclei

We present Gamow-Teller strength distributions from shell model Monte Carlo studies of fp-shell nuclei that may play an important role in the pre-collapse evolution of supernovae. We then use these strength distributions to calculate the electron-capture cross sections and rates in the zero-momentum transfer limit. We also discuss the thermal behavior of the cross sections. We find large differences in these cross sections and rates when compared to the naive single-particle estimates. These differences need to be taken into account for improved modeling of the early stages of type II supernova evolution

Complete 0 hbar omega calculations of Gamow-Teller strengths for nuclei in the iron region

Gamow-Teller strengths for selected nuclei in the iron region (A~56) have been investigated via shell-model Monte Carlo calculations with realistic interactions in the complete fp basis. Results for all cases show significant quenching relative to single-particle estimates, in quantitative agreement with (n,p) data. The J=1,T=0 residual interaction and the f_{7/2}-f_{5/2} spin-orbit splitting are shown to play major roles in the quenching mechanism. Calculated B(E2, 2^+_1 -> 0^+_1) values are in fair agreement with experiment using effective charges of e_p=1.1e and e_n=0.1e.Comment: 13 pages + 1 postscript file, Caltech preprint MAP-16