7,909 research outputs found
Application of the Contouring Method to Extended Microlensed Sources
The method devised by Lewis et al. (1993) for calculating the light curve of
a microlensed point source is expanded to two dimensions to enable the
calculation of light curves of extended sources. This method is significantly
faster than the ray shooting method that has been used in the past. The
increased efficiency is used to obtain much higher resolution light curves over
increased timescales. We investigate the signatures arising from different
source geometries in a realistic microlensing model. We show that a large
fraction of high magnification events (HMEs) in image A of Q2237+0305 involve
only one caustic, and could therefore yield information on the structure of the
quasar continuum through the recognition of a characteristic event shape. In
addition, the cataloguing of HMEs into morphological type will, in theory,
enable the direction of the transverse motion, as well as the source size to be
obtained from long term monitoring.Comment: 10 pages including 4 figures. Accepted for publication in M.N.R.A.
The combined Lagrangian advection method
We present and test a new hybrid numerical method for simulating layerwise-two-dimensional geophysical flows. The method radically extends the original Contour-Advective Semi-Lagrangian (CASL) algorithm by combining three computational elements for the advection of general tracers (e.g. potential vorticity, water vapor, etc.): (1) a pseudospectral method for large scales, (2) Lagrangian contours for intermediate to small scales, and (3) Lagrangian particles for the representation of general forcing and dissipation. The pseudo-spectral method is both efficient and highly accurate at large scales, while contour advection is efficient and accurate at small scales, allowing one to simulate extremely finescale structure well below the basic grid scale used to represent the velocity field. The particles allow one to efficiently incorporate general forcing and dissipation
Towards Real-Time Detection and Tracking of Spatio-Temporal Features: Blob-Filaments in Fusion Plasma
A novel algorithm and implementation of real-time identification and tracking
of blob-filaments in fusion reactor data is presented. Similar spatio-temporal
features are important in many other applications, for example, ignition
kernels in combustion and tumor cells in a medical image. This work presents an
approach for extracting these features by dividing the overall task into three
steps: local identification of feature cells, grouping feature cells into
extended feature, and tracking movement of feature through overlapping in
space. Through our extensive work in parallelization, we demonstrate that this
approach can effectively make use of a large number of compute nodes to detect
and track blob-filaments in real time in fusion plasma. On a set of 30GB fusion
simulation data, we observed linear speedup on 1024 processes and completed
blob detection in less than three milliseconds using Edison, a Cray XC30 system
at NERSC.Comment: 14 pages, 40 figure
Contouring with uncertainty
As stated by Johnson [Joh04], the visualization of uncertainty remains one of the major challenges for the visualization community. To achieve this, we need to understand and develop methods that allow us not only to
consider uncertainty as an extra variable within the visualization process, but to treat it as an integral part. In this paper, we take contouring, one of the most widely used visualization techniques for two dimensional data, and
focus on extending the concept of contouring to uncertainty. We develop special techniques for the visualization of uncertain contours. We illustrate the work through application to a case study in oceanography
On the non-local geometry of turbulence
A multi-scale methodology for the study of the non-local geometry of eddy structures in turbulence is developed. Starting from a given three-dimensional field, this consists of three main steps: extraction, characterization and classification of structures. The extraction step is done in two stages. First, a multi-scale decomposition based on the curvelet transform is applied to the full three-dimensional field, resulting in a finite set of component three-dimensional fields, one per scale. Second, by iso-contouring each component field at one or more iso-contour levels, a set of closed iso-surfaces is obtained that represents the structures at that scale. The characterization stage is based on the joint probability density function (p.d.f.), in terms of area coverage on each individual iso-surface, of two differential-geometry properties, the shape index and curvedness, plus the stretching parameter, a dimensionless global invariant of the surface. Taken together, this defines the geometrical signature of the iso-surface. The classification step is based on the construction of a finite set of parameters, obtained from algebraic functions of moments of the joint p.d.f. of each structure, that specify its location as a point in a multi-dimensional âfeature spaceâ. At each scale the set of points in feature space represents all structures at that scale, for the specified iso-contour value. This then allows the application, to the set, of clustering techniques that search for groups of structures with a common geometry. Results are presented of a first application of this technique to a passive scalar field obtained from 5123 direct numerical simulation of scalar mixing by forced, isotropic turbulence (ReÎť = 265). These show transition, with decreasing scale, from blob-like structures in the larger scales to blob- and tube-like structures with small or moderate stretching in the inertial range of scales, and then toward tube and, predominantly, sheet-like structures with high level of stretching in the dissipation range of scales. Implications of these results for the dynamical behaviour of passive scalar stirring and mixing by turbulence are discussed
- âŚ