Ellipse Fitting Using the Finite Rate of Innovation Sampling Principle

Standard approaches for ellipse fitting are based on the minimization of algebraic or geometric distance between the given data and a template ellipse. When the data are noisy and come from a partial ellipse, the state-of-the-art methods tend to produce biased ellipses. We rely on the sampling structure of the underlying signal and show that the x- and y-coordinate functions of an ellipse are finite-rate-of-innovation (FRI) signals, and that their parameters are estimable from partial data. We consider both uniform and nonuniform sampling scenarios in the presence of noise and show that the data can be modeled as a sum of random amplitude-modulated complex exponentials. A low-pass filter is used to suppress noise and approximate the data as a sum of weighted complex exponentials. The annihilating filter used in FRI approaches is applied to estimate the sampling interval in the closed form. We perform experiments on simulated and real data, and assess both objective and subjective performances in comparison with the state-of-the-art ellipse fitting methods. The proposed method produces ellipses with lesser bias. Furthermore, the mean-squared error is lesser by about 2 to 10 dB. We show the applications of ellipse fitting in iris images starting from partial edge contours, and to free-hand ellipses drawn on a touch-screen tablet.

Characteristics of Micronozzle Gas Flows

The fluid characteristics of gas flows in the micronozzle whose throat height is 20 μm were investigated by the direct simulation Monte Carlo (DSMC) method. In a series of cases, the dependence of mass flux on the pressure difference was gained, and the DSMC's results show good agreement with the experimental data. The comparison of mass flux and the Mach number contours between the DSMC and Navier-Stokes equations adding slip boundary also reveals quantitatively that the continuum model will be invalid gradually even when the average Knudsen number is smaller than 0.01. As one focus of the present paper, the phenomenon of the multiple expansion-compression waves that comes from the nozzle's divergent part was analyzed in detailed.

MILU-CG method and the numerical study on the flow around a rotating circular cylinder

A hybrid finite difference method and vortex method (HDV), which is based on domain decomposition and proposed by the authors (1992), is improved by using a modified incomplete LU decomposition conjugate gradient method (MILU-CG), and a high order implicit difference algorithm. The flow around a rotating circular cylinder at Reynolds number R-e = 1000, 200 and the angular to rectilinear speed ratio alpha is an element of (0.5, 3.25) is studied numerically. The long-time full developed features about the variations of the vortex patterns in the wake, and drag, lift forces on the cylinder are given. The calculated streamline contours agreed well with the experimental visualized flow pictures. The existence of critical states and the vortex patterns at the states are given for the first time. The maximum lift to drag force ratio can be obtained nearby the critical states.

Surface waves in a vertically excited circular cylindrical container

The nonlinear free surface amplitude equation, which has been derived from the inviscid fluid by solving the potential equation of water waves with a singular perturbation theory in a vertically oscillating rigid circular cylinder, is investigated successively in the fourth-order Runge-Kutta approach with an equivalent time-step. Computational results include the evolution of the amplitude with time, the characteristics of phase plane determined by the real and imaginary parts of the amplitude, the single-mode selection rules of the surface waves in different forced frequencies, contours of free surface displacement and corresponding three-dimensional evolution of surface waves, etc. In addition, the comparison of the surface wave modes is made between theoretical calculations and experimental measurements, and the results are reasonable although there are some differences in the forced frequency.

Experiencia social y tradición en el derecho penal

Resumen: El delito penal es aquello que la tradición, vivida en la experiencia de la familia y de la comunidad, permite a cada uno reconocer como un grave alejamiento de lo verdadero, lo bueno y lo correcto. El fundamento de la punibilidad penal es la imputación, el reconocimiento de la pertenencia del delito al sujeto como a su causa. Se intenta mostrar cómo este vínculo del derecho penal con la tradición jurídica sufrió dos rupturas: con el iluminismo jurídico y el kantismo se separó la imputación jurídica de su fundamento moral y con la codificación, se quebró la unidad del derecho penal universal fundado en lo bonum et aequum otorgando prevalencia a la idea de la legalidad esclava de los intereses de los Estados. Al mismo tiempo, se produjo una segunda ruptura: se predica una responsabilidad que se atribuye desde afuera, en forma objetiva, a centros de imputación –con frecuencia colectivos– que realizan la producción industrial. Así, según exigencias de la seguridad y la salud, el derecho penal se transforma en un instrumento de la política criminal y los contornos del tipo penal se operan a través de jueces y fiscales en prevención de las consecuencias futuras del “riesgo” progresivo de la producción industrial. La nota dominante es el “riesgo” colectivo y no el “hecho”.

Effects of temperature on the threshold of phosphorus for algal blooms

Algal blooms, worsening marine ecosystems and causing great economic loss, have been paid much attention to for a long time. Such environmental factors as light penetration, water temperature, and nutrient concentration are crucial in blooms processes. Among them, only nutrients can be controlled. Therefore, the threshold of nutrients for algal blooms is of great concern. To begin with, a dynamic eutrophication model has been constructed to simulate the algal growth and phosphorus cycling. The model encapsulates the essential biological processes of algal growth and decay, and phosphorus regeneration due to algal decay. The nutrient limitation is based upon commonly used Monod's kinetics. The effects of temperature and phosphorus limitation are particularly addressed. Then, we have endeavored to elucidate the threshold of phosphorus at different temperature for algal blooms. Based on the numerical simulation, the isoquant contours of change rate of alga as shown in the figure are obtained, which obviously demonstrate the threshold of nutrient at an arbitrary reasonable temperature. The larger the change rate is, the more rapidly the alga grows. If the phosphorus concentration at a given temperature remains larger than the threshold the algal biomass may increase monotonically, leading to the algal blooming. With the rising of temperature, the threshold is apparently reduced, which may explain why likely red tide disasters occur in a fine summer day. So, high temperature and sufficient phosphorus supply are the major factors which result in algal growth and blowout of red tide.

Space-time correlations of fluctuating velocities in turbulent shear flows

Space-time correlations or Eulerian two-point two-time correlations of fluctuating velocities are analytically and numerically investigated in turbulent shear flows. An elliptic model for the space-time correlations in the inertial range is developed from the similarity assumptions on the isocorrelation contours: they share a uniform preference direction and a constant aspect ratio. The similarity assumptions are justified using the Kolmogorov similarity hypotheses and verified using the direct numerical simulation DNS of turbulent channel flows. The model relates the space-time correlations to the space correlations via the convection and sweeping characteristic velocities. The analytical expressions for the convection and sweeping velocities are derived from the Navier-Stokes equations for homogeneous turbulent shear flows, where the convection velocity is represented by the mean velocity and the sweeping velocity is the sum of the random sweeping velocity and the shearinduced velocity. This suggests that unlike Taylor’s model where the convection velocity is dominating and Kraichnan and Tennekes’ model where the random sweeping velocity is dominating, the decorrelation time scales of the space-time correlations in turbulent shear flows are determined by the convection velocity, the random sweeping velocity, and the shear-induced velocity. This model predicts a universal form of the spacetime correlations with the two characteristic velocities. The DNS of turbulent channel flows supports the prediction: the correlation functions exhibit a fair good collapse, when plotted against the normalized space and time separations defined by the elliptic model.

Generalization of the JTZ model to open plane wakes

The JTZ model [C. Jung, T. T¶el and E. Ziemniak, Chaos 3, (1993) 555], as a theoretical model of a plane wake behind a circular cylinder in a narrow channel at a moderate Reynolds number, has previously been employed to analyze phenomena of chaotic scattering. It is ex- tended here to describe an open plane wake without the con¯ned nar- row channel by incorporating a double row of shedding vortices into the intermediate and far wake. The extended JTZ model is found in qualitative agreement with both direct numerical simulations and ex- perimental results in describing streamlines and vorticity contours. To further validate its applications to particle transport processes, the in- teraction between small spherical particles and vortices in an extended JTZ model °ow is studied. It is shown that the particle size has signif- icant in°uences on the features of particle trajectories, which have two characteristic patterns: one is rotating around the vortex centers and the other accumulating in the exterior of vortices. Numerical results based on the extended JTZ model are found in qualitative agreement with experimental ones in the normal range of particle sizes.

Numerical Study on Critical Wedge Angle of Cellular Detonation Reflections

The critical wedge angle (CWA) for the transition from regular reflection (RR) to Mach reflection (MR) of a cellular detonation wave is studied numerically by an improved space-time conservation element and solution element method together with a two-step chemical reaction model. The accuracy of that numerical way is verified by simulating cellular detonation reflections at a 19.3∘ wedge. The planar and cellular detonation reflections over 45∘–55∘ wedges are also simulated. When the cellular detonation wave is over a 50∘ wedge, numerical results show a new phenomenon that RR and MR occur alternately. The transition process between RR and MR is investigated with the local pressure contours. Numerical analysis shows that the cellular structure is the essential reason for the new phenomenon and the CWA of detonation reflection is not a certain angle but an angle range.

A dynamical systems analysis of vortex pinch-off

Vortex rings constitute the main structure in the wakes of a wide class of swimming and flying animals, as well as in cardiac flows and in the jets generated by some moss and fungi. However, there is a physical limit, determined by an energy maximization principle called the Kelvin-Benjamin principle, to the size that axisymmetric vortex rings can achieve. The existence of this limit is known to lead to the separation of a growing vortex ring from the shear layer feeding it, a process known as `vortex pinch-off', and characterized by the dimensionless vortex formation number. The goal of this thesis is to improve our understanding of vortex pinch-off as it relates to biological propulsion, and to provide future researchers with tools to assist in identifying and predicting pinch-off in biological flows.

To this end, we introduce a method for identifying pinch-off in starting jets using the Lagrangian coherent structures in the flow, and apply this criterion to an experimentally generated starting jet. Since most naturally occurring vortex rings are not circular, we extend the definition of the vortex formation number to include non-axisymmetric vortex rings, and find that the formation number for moderately non-axisymmetric vortices is similar to that of circular vortex rings. This suggests that naturally occurring vortex rings may be modeled as axisymmetric vortex rings. Therefore, we consider the perturbation response of the Norbury family of axisymmetric vortex rings. This family is chosen to model vortex rings of increasing thickness and circulation, and their response to prolate shape perturbations is simulated using contour dynamics. Finally, the response of more realistic models for vortex rings, constructed from experimental data using nested contours, to perturbations which resemble those encountered by forming vortices more closely, is simulated using contour dynamics. In both families of models, a change in response analogous to pinch-off is found as members of the family with progressively thicker cores are considered. We posit that this analogy may be exploited to understand and predict pinch-off in complex biological flows, where current methods are not applicable in practice, and criteria based on the properties of vortex rings alone are necessary.

Experimental Investigation of Thrust Fault Rupture Mechanics

Thrust fault earthquakes are investigated in the laboratory by generating dynamic shear ruptures along pre-existing frictional faults in rectangular plates. A considerable body of evidence suggests that dip-slip earthquakes exhibit enhanced ground motions in the acute hanging wall wedge as an outcome of broken symmetry between hanging and foot wall plates with respect to the earth surface. To understand the physical behavior of thrust fault earthquakes, particularly ground motions near the earth surface, ruptures are nucleated in analog laboratory experiments and guided up-dip towards the simulated earth surface. The transient slip event and emitted radiation mimic a natural thrust earthquake. High-speed photography and laser velocimeters capture the rupture evolution, outputting a full-field view of photo-elastic fringe contours proportional to maximum shearing stresses as well as continuous ground motion velocity records at discrete points on the specimen. Earth surface-normal measurements validate selective enhancement of hanging wall ground motions for both sub-Rayleigh and super-shear rupture speeds. The earth surface breaks upon rupture tip arrival to the fault trace, generating prominent Rayleigh surface waves. A rupture wave is sensed in the hanging wall but is, however, absent from the foot wall plate: a direct consequence of proximity from fault to seismometer. Signatures in earth surface-normal records attenuate with distance from the fault trace. Super-shear earthquakes feature greater amplitudes of ground shaking profiles, as expected from the increased tectonic pressures required to induce super-shear transition. Paired stations measure fault parallel and fault normal ground motions at various depths, which yield slip and opening rates through direct subtraction of like components. Peak fault slip and opening rates associated with the rupture tip increase with proximity to the fault trace, a result of selective ground motion amplification in the hanging wall. Fault opening rates indicate that the hanging and foot walls detach near the earth surface, a phenomenon promoted by a decrease in magnitude of far-field tectonic loads. Subsequent shutting of the fault sends an opening pulse back down-dip. In case of a sub-Rayleigh earthquake, feedback from the reflected S wave re-ruptures the locked fault at super-shear speeds, providing another mechanism of super-shear transition.

Computational modeling and psychophysics in low- and mid-level vision

This thesis addresses a series of topics related to the question of how people find the foreground objects from complex scenes. With both computer vision modeling, as well as psychophysical analyses, we explore the computational principles for low- and mid-level vision.

We first explore the computational methods of generating saliency maps from images and image sequences. We propose an extremely fast algorithm called Image Signature that detects the locations in the image that attract human eye gazes. With a series of experimental validations based on human behavioral data collected from various psychophysical experiments, we conclude that the Image Signature and its spatial-temporal extension, the Phase Discrepancy, are among the most accurate algorithms for saliency detection under various conditions.

In the second part, we bridge the gap between fixation prediction and salient object segmentation with two efforts. First, we propose a new dataset that contains both fixation and object segmentation information. By simultaneously presenting the two types of human data in the same dataset, we are able to analyze their intrinsic connection, as well as understanding the drawbacks of today’s “standard” but inappropriately labeled salient object segmentation dataset. Second, we also propose an algorithm of salient object segmentation. Based on our novel discoveries on the connections of fixation data and salient object segmentation data, our model significantly outperforms all existing models on all 3 datasets with large margins.

In the third part of the thesis, we discuss topics around the human factors of boundary analysis. Closely related to salient object segmentation, boundary analysis focuses on delimiting the local contours of an object. We identify the potential pitfalls of algorithm evaluation for the problem of boundary detection. Our analysis indicates that today’s popular boundary detection datasets contain significant level of noise, which may severely influence the benchmarking results. To give further insights on the labeling process, we propose a model to characterize the principles of the human factors during the labeling process.

The analyses reported in this thesis offer new perspectives to a series of interrelating issues in low- and mid-level vision. It gives warning signs to some of today’s “standard” procedures, while proposing new directions to encourage future research.

透镜自重变形引起波像差的有限元分析

在用大口径、长焦距平行光管模拟激光远场特性时，其光组透镜在重力作用下的变形不能忽略，为了分析对出射光束质量的影响，采用有限元分析软件“ANSYS”建立了平行光管光组中声400mm平凸透镜的有限元模型，给出一种分析透镜轴向变形引起的波像差的方法，在不同工况下，计算了平凸透镜在重力作用下轴向变形的峰谷值和均方根值，对轴向变形量均方根值最小工况画出了透镜表面变形的等值线图，计算了声350mm通光口径内的波像差峰谷值和均方根值，对平行光管光组的波像差做出估计，验证了设计的合理性。

Absolute left-handed behaviors in a triangular elliptical-rod photonic crystal

In this paper we theoretically study the left-handed behaviors in a two-dimensional triangular photonic crystal made of elliptical rods in air. An absolute left-handed region is found in the second photonic band by using the plane wave expansion method to analyze the photonic band structure and equifrequency contours. Typical left-handed behaviors such as negative refraction, flat superlensing and plano-concave lensing are demonstrated by the finite-difference time-domain simulations. These behaviors are also compared with the quasi-negative refraction and the resulted focusing effects in a square-lattice two-dimensional photonic crystal. (c) 2005 Optical Society of America

Scattering by two degenerate anisotropic modes in square-lattice dielectric photonic crystals

We systematically investigate the square-lattice dielectric photonic crystals that have been used to demonstrate flat slab imaging experimentally. A right-handed Bloch mode is found in the left-handed frequency region by using the plane wave expansion method to analyze the photonic band structure and equifrequency contours. Using the multiple scattering theory, numerical simulations demonstrate that the left-handed mode and the right-handed mode are excited simultaneously by a point source and result in two kinds of transmitted waves. Impacted by the evanescent waves, superposition of these transmitted waves brings on complicated near field distributions such as the so-called imaging and its disappearance.