A Phase Transition for the Uniform Distribution in the Pattern Maximum Likelihood Problem

In this paper, we consider the setting of the pattern maximum likelihood (PML) problem studied by Orlitsky et al. We present a well-motivated heuristic algorithm for deciding the question of when the PML distribution of a given pattern is uniform. The algorithm is based on the concept of a ``uniform threshold''. This is a threshold at which the uniform distribution exhibits an interesting phase transition in the PML problem, going from being a local maximum to being a local minimum.

Evaluation of airline exercises prescribed to avoid deep vein thrombosis using fiber Bragg grating sensors

Various leg exercises have been recommended to prevent deep vein thrombosis (DVT), a condition where a blood clot forms in the deep veins, especially during long-haul flights. Accessing the benefit of each of these exercises in avoiding the DVT, which can be fatal, is important in the context of suggesting the correct and the most beneficial exercises. Present work aims at demonstrating the fiber Bragg grating (FBG)-based sensing methodology for measuring surface strains generated on the skin of the calf muscle to evaluate the suggested airline exercises to avoid DVT. As the dataset in the experiment involves multiple subjects performing these exercises, an inertial measurement unit has been used to validate the repetitiveness of each of the exercises. The surface strain on the calf muscle obtained using the FBG sensor, which is a measure of the calf muscle deformation, has been compared against the variation of blood velocity in the femoral vein of the thigh measured using a commercial electronic-phased array color Doppler ultrasound system. Apart from analyzing the effectiveness of suggested exercises, a new exercise which is more effective in terms of strain generated to avoid DVT is proposed and evaluated. (C) 2013 Society of Photo-Optical Instrumentation Engineers (SPIE)

A Kushner-Stratonovich Monte Carlo filter applied to nonlinear dynamical system identification

A Monte Carlo filter, based on the idea of averaging over characteristics and fashioned after a particle-based time-discretized approximation to the Kushner-Stratonovich (KS) nonlinear filtering equation, is proposed. A key aspect of the new filter is the gain-like additive update, designed to approximate the innovation integral in the KS equation and implemented through an annealing-type iterative procedure, which is aimed at rendering the innovation (observation prediction mismatch) for a given time-step to a zero-mean Brownian increment corresponding to the measurement noise. This may be contrasted with the weight-based multiplicative updates in most particle filters that are known to precipitate the numerical problem of weight collapse within a finite-ensemble setting. A study to estimate the a-priori error bounds in the proposed scheme is undertaken. The numerical evidence, presently gathered from the assessed performance of the proposed and a few other competing filters on a class of nonlinear dynamic system identification and target tracking problems, is suggestive of the remarkably improved convergence and accuracy of the new filter. (C) 2013 Elsevier B.V. All rights reserved.

Channel estimation relying on the minimum bit-error-ratio criterion for BPSK and QPSK signals

The authors consider the channel estimation problem in the context of a linear equaliser designed for a frequency selective channel, which relies on the minimum bit-error-ratio (MBER) optimisation framework. Previous literature has shown that the MBER-based signal detection may outperform its minimum-mean-square-error (MMSE) counterpart in the bit-error-ratio performance sense. In this study, they develop a framework for channel estimation by first discretising the parameter space and then posing it as a detection problem. Explicitly, the MBER cost function (CF) is derived and its performance studied, when transmitting binary phase shift keying (BPSK) and quadrature phase shift keying (QPSK) signals. It is demonstrated that the MBER based CF aided scheme is capable of outperforming existing MMSE, least square-based solutions.

Curvature based mobility analysis and form closure of smooth planar curves with multiple contacts

This paper presents a simple second-order, curvature based mobility analysis of planar curves in contact. The underlying theory deals with penetration and separation of curves with multiple contacts, based on relative configuration of osculating circles at points of contact for a second-order rotation about each point of the plane. Geometric and analytical treatment of mobility analysis is presented for generic as well as special contact geometries. For objects with a single contact, partitioning of the plane into four types of mobility regions has been shown. Using point based composition operations based on dual-number matrices, analysis has been extended to computationally handle multiple contacts scenario. A novel color coded directed line has been proposed to capture the contact scenario. Multiple contacts mobility is obtained through intersection of the mobility half-spaces. It is derived that mobility region comprises a pair of unbounded or a single bounded convex polygon. The theory has been used for analysis and synthesis of form closure configurations, revolute and prismatic kinematic pairs. (C) 2013 Elsevier Ltd. All rights reserved.

Spatially Adaptive Kernel Regression Using Risk Estimation

An important question in kernel regression is one of estimating the order and bandwidth parameters from available noisy data. We propose to solve the problem within a risk estimation framework. Considering an independent and identically distributed (i.i.d.) Gaussian observations model, we use Stein's unbiased risk estimator (SURE) to estimate a weighted mean-square error (MSE) risk, and optimize it with respect to the order and bandwidth parameters. The two parameters are thus spatially adapted in such a manner that noise smoothing and fine structure preservation are simultaneously achieved. On the application side, we consider the problem of image restoration from uniform/non-uniform data, and show that the SURE approach to spatially adaptive kernel regression results in better quality estimation compared with its spatially non-adaptive counterparts. The denoising results obtained are comparable to those obtained using other state-of-the-art techniques, and in some scenarios, superior.

Free turbulent shear layer in a point vortex gas as a problem in nonequilibrium statistical mechanics

This paper attempts to unravel any relations that may exist between turbulent shear flows and statistical mechanics through a detailed numerical investigation in the simplest case where both can be well defined. The flow considered for the purpose is the two-dimensional (2D) temporal free shear layer with a velocity difference Delta U across it, statistically homogeneous in the streamwise direction (x) and evolving from a plane vortex sheet in the direction normal to it (y) in a periodic-in-x domain L x +/-infinity. Extensive computer simulations of the flow are carried out through appropriate initial-value problems for a ``vortex gas'' comprising N point vortices of the same strength (gamma = L Delta U/N) and sign. Such a vortex gas is known to provide weak solutions of the Euler equation. More than ten different initial-condition classes are investigated using simulations involving up to 32 000 vortices, with ensemble averages evaluated over up to 10(3) realizations and integration over 10(4)L/Delta U. The temporal evolution of such a system is found to exhibit three distinct regimes. In Regime I the evolution is strongly influenced by the initial condition, sometimes lasting a significant fraction of L/Delta U. Regime III is a long-time domain-dependent evolution towards a statistically stationary state, via ``violent'' and ``slow'' relaxations P.-H. Chavanis, Physica A 391, 3657 (2012)], over flow time scales of order 10(2) and 10(4)L/Delta U, respectively (for N = 400). The final state involves a single structure that stochastically samples the domain, possibly constituting a ``relative equilibrium.'' The vortex distribution within the structure follows a nonisotropic truncated form of the Lundgren-Pointin (L-P) equilibrium distribution (with negatively high temperatures; L-P parameter lambda close to -1). The central finding is that, in the intermediate Regime II, the spreading rate of the layer is universal over the wide range of cases considered here. The value (in terms of momentum thickness) is 0.0166 +/- 0.0002 times Delta U. Regime II, extensively studied in the turbulent shear flow literature as a self-similar ``equilibrium'' state, is, however, a part of the rapid nonequilibrium evolution of the vortex-gas system, which we term ``explosive'' as it lasts less than one L/Delta U. Regime II also exhibits significant values of N-independent two-vortex correlations, indicating that current kinetic theories that neglect correlations or consider them as O(1/N) cannot describe this regime. The evolution of the layer thickness in present simulations in Regimes I and II agree with the experimental observations of spatially evolving (3D Navier-Stokes) shear layers. Further, the vorticity-stream-function relations in Regime III are close to those computed in 2D Navier-Stokes temporal shear layers J. Sommeria, C. Staquet, and R. Robert, J. Fluid Mech. 233, 661 (1991)]. These findings suggest the dominance of what may be called the Kelvin-Biot-Savart mechanism in determining the growth of the free shear layer through large-scale momentum and vorticity dispersal.

Cubic Sieve Congruence of the Discrete Logarithm Problem, and fractional part sequences

The Cubic Sieve Method for solving the Discrete Logarithm Problem in prime fields requires a nontrivial solution to the Cubic Sieve Congruence (CSC) x(3) equivalent to y(2)z (mod p), where p is a given prime number. A nontrivial solution must also satisfy x(3) not equal y(2)z and 1 <= x, y, z < p(alpha), where alpha is a given real number such that 1/3 < alpha <= 1/2. The CSC problem is to find an efficient algorithm to obtain a nontrivial solution to CSC. CSC can be parametrized as x equivalent to v(2)z (mod p) and y equivalent to v(3)z (mod p). In this paper, we give a deterministic polynomial-time (O(ln(3) p) bit-operations) algorithm to determine, for a given v, a nontrivial solution to CSC, if one exists. Previously it took (O) over tilde (p(alpha)) time in the worst case to determine this. We relate the CSC problem to the gap problem of fractional part sequences, where we need to determine the non-negative integers N satisfying the fractional part inequality {theta N} < phi (theta and phi are given real numbers). The correspondence between the CSC problem and the gap problem is that determining the parameter z in the former problem corresponds to determining N in the latter problem. We also show in the alpha = 1/2 case of CSC that for a certain class of primes the CSC problem can be solved deterministically in <(O)over tilde>(p(1/3)) time compared to the previous best of (O) over tilde (p(1/2)). It is empirically observed that about one out of three primes is covered by the above class. (C) 2013 Elsevier B.V. All rights reserved.

Enhanced Gas Adsorption on Graphitic Substrates via Defects and Local Curvature: A Density Functional Theory Study

Using van-der-Waals-corrected density functional theory calculations, we explore the possibility of engineering the local structure and morphology of high-surface-area graphene-derived materials to improve the uptake of methane and carbon dioxide for gas storage and sensing. We test the sensitivity of the gas adsorption energy to the introduction of native point defects, curvature, and the application of strain. The binding energy at topological point defect sites is inversely correlated with the number of missing carbon atoms, causing Stone-Wales defects to show the largest enhancement with respect to pristine graphene (similar to 20%). Improvements of similar magnitude are observed at concavely curved surfaces in buckled graphene sheets under compressive strain, whereas tensile strain tends to weaken gas binding. Trends for CO2 and CH4 are, similar, although CO2 binding is generally stronger by similar to 4 to 5 kJ mol(-1). However, the differential between the adsorption of CO2 and CH4 is much higher on folded graphene sheets and at concave curvatures; this could possibly be leveraged for CH4/CO2 flow separation and gasselective sensors.

Error analysis of discontinuous Galerkin methods for the Stokes problem under minimal regularity

In this article, we analyse several discontinuous Galerkin (DG) methods for the Stokes problem under minimal regularity on the solution. We assume that the velocity u belongs to H-0(1)(Omega)](d) and the pressure p is an element of L-0(2)(Omega). First, we analyse standard DG methods assuming that the right-hand side f belongs to H-1(Omega) boolean AND L-1(Omega)](d). A DG method that is well defined for f belonging to H-1(Omega)](d) is then investigated. The methods under study include stabilized DG methods using equal-order spaces and inf-sup stable ones where the pressure space is one polynomial degree less than the velocity space.

Fluctuating micro-heterogeneity in water-tert-butyl alcohol mixtures and lambda-type divergence of the mean cluster size with phase transition-like multiple anomalies

Water-tert-butyl alcohol (TBA) binary mixture exhibits a large number of thermodynamic and dynamic anomalies. These anomalies are observed at surprisingly low TBA mole fraction, with x(TBA) approximate to 0.03-0.07. We demonstrate here that the origin of the anomalies lies in the local structural changes that occur due to self-aggregation of TBA molecules. We observe a percolation transition of the TBA molecules at x(TBA) approximate to 0.05. We note that ``islands'' of TBA clusters form even below this mole fraction, while a large spanning cluster emerges above that mole fraction. At this percolation threshold, we observe a lambda-type divergence in the fluctuation of the size of the largest TBA cluster, reminiscent of a critical point. Alongside, the structure of water is also perturbed, albeit weakly, by the aggregation of TBA molecules. There is a monotonic decrease in the tetrahedral order parameter of water, while the dipole moment correlation shows a weak nonlinearity. Interestingly, water molecules themselves exhibit a reverse percolation transition at higher TBA concentration, x(TBA) approximate to 0.45, where large spanning water clusters now break-up into small clusters. This is accompanied by significant divergence of the fluctuations in the size of largest water cluster. This second transition gives rise to another set of anomalies around. Both the percolation transitions can be regarded as manifestations of Janus effect at small molecular level. (C) 2014 AIP Publishing LLC.

Sensitivity of simulated climate to latitudinal distribution of solar insolation reduction in solar radiation management

Solar radiation management (SRM) geoengineering has been proposed as a potential option to counteract climate change. We perform a set of idealized geoengineering simulations using Community Atmosphere Model version 3.1 developed at the National Center for Atmospheric Research to investigate the global hydrological implications of varying the latitudinal distribution of solar insolation reduction in SRM methods. To reduce the solar insolation we have prescribed sulfate aerosols in the stratosphere. The radiative forcing in the geoengineering simulations is the net forcing from a doubling of CO2 and the prescribed stratospheric aerosols. We find that for a fixed total mass of sulfate aerosols (12.6 Mt of SO4), relative to a uniform distribution which nearly offsets changes in global mean temperature from a doubling of CO2, global mean radiative forcing is larger when aerosol concentration is maximum at the poles leading to a warmer global mean climate and consequently an intensified hydrological cycle. Opposite changes are simulated when aerosol concentration is maximized in the tropics. We obtain a range of 1 K in global mean temperature and 3% in precipitation changes by varying the distribution pattern in our simulations: this range is about 50% of the climate change from a doubling of CO2. Hence, our study demonstrates that a range of global mean climate states, determined by the global mean radiative forcing, are possible for a fixed total amount of aerosols but with differing latitudinal distribution. However, it is important to note that this is an idealized study and thus not all important realistic climate processes are modeled.

Low-frequency modes in an equatorial shallow-water model with moisture gradients

Motivated by observations of the mean state of tropical precipitable water (PW), a moist, first baroclinic mode, shallow-water system on an equatorial beta-plane with a background saturation profile that depends on latitude and longitude is studied. In the presence of a latitudinal moisture gradient, linear analysis of the non-rotating problem reveals large-scale, symmetric, eastward and westward propagating unstable modes. The introduction of a zonal moisture gradient breaks the east-west symmetry of the unstable modes. The effects of rotation are then included by numerically solving the resulting eigenvalue problem on an equatorial beta-plane. With a purely meridional moisture gradient, the system supports large-scale, low-frequency, eastward and westward moving neutral modes. Some of the similarities, and some of the discrepancies of these modes with intraseasonal tropical waves are pointed out. Finally, a zonal moisture gradient in the presence of rotation renders some of the aforementioned neutral modes unstable. In particular, according to observations of large-scale, low-frequency tropical variability, it is seen that regions where the background saturation profile increases (decreases) to the east favour eastward (westward) moving moist modes.

On the stability of the L-p-norm of the Riemannian curvature tensor

We consider the Riemannian functional defined on the space of Riemannian metrics with unit volume on a closed smooth manifold M where R(g) and dv (g) denote the corresponding Riemannian curvature tensor and volume form and p a (0, a). First we prove that the Riemannian metrics with non-zero constant sectional curvature are strictly stable for for certain values of p. Then we conclude that they are strict local minimizers for for those values of p. Finally generalizing this result we prove that product of space forms of same type and dimension are strict local minimizer for for certain values of p.