Gibbs sampling for time-delay-and amplitude estimation in underwater acoustics

Multipath arrivals at a receiving sensor are frequently encountered in many signal-processing areas, including sonar, radar, and communication problems. In underwater acoustics, numerous approaches to source localization, geoacoustic inversion, and tomography rely on accurate multipath arrival extraction. A novel method for estimation of time delays and amplitudes of arrivals with maximum a posteriori (MAP) estimation is presented here. MAP estimation is optimal if appropriate statistical models are selected for the data; implementation, requiring maximization of a multidimensional function, is computationally demanding. Gibbs sampling is proposed as an efficient means for estimating necessary posterior probability distributions, bypassing analytical calculations. The Gibbs sampler includes as unknowns time delays, amplitudes, noise variance, and number of arrivals. Through Monte Carlo simulations, the method is shown to have a performance very close to that of analytical MAP estimation. The method is also shown to be superior to expectation-maximization, which is often applied to time-delay estimation. The Gibbs sampling approach is demonstrated to be more informative than other time-delay estimation methods, providing complete posterior distributions compared to just point estimates; the distributions capture the uncertainty in the problem, presenting likely values of the unknowns that are different from simple point estimates.     

Properties of balanced and perfect matrices

In this paper we define wheel matrices and characterize some properties of matrices that are perfect but not balanced. A consequence of our results is that if a matrixA is perfect then we can construct a polynomial number of submatricesA _I,,A _n ofA such thatA is balanced if and only if all the submatricesA ₁,,A _n ofA are perfect. This implies that if the problem of testing perfection is polynomially solvable, then so is the problem of testing balancedness.Partial support under NSF Grants DMS-8606188, ECS-8800281 and DDM-8800281.     

Ab initio path-integral molecular dynamics

A path-integral molecular dynamics technique for strongly interacting atoms using ab initio potentials derived from density functional theory is implemented. This allows the efficient inclusion of nuclear quantum dispersion in ab initio simulations at finite temperatures. We present an application to the quantum cluster H ₅ ⁺ .     

Integer linear programming models for 2-staged two-dimensional Knapsack problems

 We are given a unique rectangular piece of stock material S, with height H and width W, and a list of m rectangular shapes to be cut from S. Each shape's type i (i = 1, ..., m) is characterized by a height , a width , a profit , and an upper bound ub _i indicating the maximum number of items of type i which can be cut. We refer to the Two-Dimensional Knapsack (TDK) as the problem of determining a cutting pattern of S maximizing the sum of the profits of the cut items. In particular, we consider the classical variant of TDK in which the maximum number of cuts allowed to obtain each item is fixed to 2, and we refer to this problem as 2-staged TDK (2TDK). For the 2TDK problem we present two new Integer Linear Programming models, we discuss their properties, and we compare them with other formulations in terms of the LP bound they provide. Finally, both models are computationally tested within a standard branch-and-bound framework on a large set of instances from the literature by reinforcing them with the addition of linear inequalities to eliminate symmetries. Received: October 17, 2000 / Accepted: December 19, 2001 Published online: September 27, 2002 Key words. packing – cutting – integer linear programming     

Strongly correlated superconductivity and pseudogap phase near a multiband Mott insulator

Near a Mott transition, strong electron correlations may enhance Cooper pairing. This is demonstrated in the dynamical mean field theory solution of a twofold-orbital degenerate Hubbard model with an inverted on-site Hund rule exchange, favoring local spin-singlet configurations. Close to the Mott insulator (which here is a local version of a valence bond insulator) a pseudogap non-Fermi-liquid metal, a superconductor, and a normal metal appear, in striking similarity with the physics of cuprates. The strongly correlated s-wave superconducting state has a larger Drude weight than the corresponding normal state. The role of the impurity Kondo problem is underscored.     

Superprism Behaviour of an Array of Photonic Crystal Waveguides

The behaviour of an array of photonic crystal waveguides is numerically investigated. It is shown that high dispersion may be achieved in the telecommunication window around 1550 nm, with a device whose dimensions are in the order of half a mm².     

Delayed feedback control of phase synchronisation in a neuronal network model

The human cerebral cortex can be separated into cortical areas forming a clustered network structure. We build two different clustered networks, where one network is based on a healthy brain and the other according to a brain affected by a neurodegenerative process. Each cortical area has a subnetwork with small-world properties. We verify that both networks exhibit rich-club organisation and phase synchronisation. Due to the fact that neuronal synchronisation can be related to brain diseases, we consider the delayed feedback control as a method to suppress synchronous behaviours. In this work, it is presented that depending on the feedback parameters, intensity and time delay, phase synchronisation in both networks can be suppressed. Therefore, one of our main results is to show that delayed feedback control can be used to suppress undesired synchronous behaviours not only in the healthy brain, but also in the brain marked by neurodegenerative processes.     

Chemical modification of chitosan under high-intensity ultrasound

Chitosan (CTS), a biocompatible, biodegradable, non-toxic polymer, dissolves in water only if pH is lowered under 6.5, when a substantial fraction of the amino groups is protonated. Its range of application has been much extended by partially depolymerising it or converting it to water-soluble derivatives. Working under high-intensity ultrasound at 17.8-18.5 kHz, using either a simple horn or a cup horn, we achieved a controlled depolymerization of CTS, also prepared in high yields several derivatives that can be useful intermediates for further chemical modification, as well as several water-soluble derivatives that lend themselves to a host of industrial applications. Compared to conventional methods, all these reactions went to completion in considerably shorter times at lower temperatures.     

Experimental and numerical investigation of a stagnation pulverised coal flame

A multi-stream Flamelet Progress Variable (FPV) model, specifically developed for coal combustion, is proposed. The model accounts for the different fuel streams associated with the volatile and char burnout products. The applicability of the new FPV model is investigated in a laminar stagnation pulverized coal flame. The flame considered is a premixed mixture of CH₄, O₂ and N₂, carrying pulverized coal particles, stabilized in an impinging wall. Spontaneous emissions of OH*, CH* and C${}_2^{*}$ are measured to identify the flame. The 1D numerical simulations of the experimental conditions are able to reproduce the main features of the flame. The applicability of the multi-stream FPV model to coal combustion is further evaluated with the a posteriori analysis of the FPV results, comparing the results with a reference model, where the species are fully transported and the chemistry directly evaluated. Then, with the budget analysis, the influence of the control variables used to build the look-up table is assessed by examining the conditional contributions to the overall transport terms of scalar quantities (e.g. species, temperature). The results of both analyses show that the proposed multi-stream FPV model can accurately predict the main features of coal combustion, with only minor issues related to the manifold used to build the look-up table.