Fluid-Structure Interaction by the Spectral Element Method

Viscous fluid-structure interaction is treated with an arbitrary Lagrangian- Eulerian formulation. The spatial discretization is performed by the spectral element method for the fluid part where the Navier-Stokes equations are integrated and in the solid part where transient linear elasticity is described by the Navier equations. Time marching algorithms are second-order accurate in time in both the fluid and the solid. The algorithm is applied to the flow in a plane channel partially obstructed by a solid component able to move under the action of the fluid flo

On the role of fine-sand dune dynamics in controlling water depth changes in Rio Parapeti, Serrania Borebigua (Southern sub-Andean zone of Bolivia)

The role of the fine-dune sand dynamics in controlling the natural regeneration of the upper layer of a riverbed used for filtration is studied at the Choreti test reach of Rio Parapeti, in the Southern sub-Andean zone of Bolivia. Local production of drinking water relies on Riverbed Filtration, the delivery of which depends on the river water depth and the riverbed permeability. There is a strong, natural, declamation process of the upper layer maintained by dune bed-forms migrating downstream. It is thus essential to understand and represent local water depth changes as a function of the incoming discharge. We show the vortex-drag model can be used to correctly calculate the stream velocity in natural environment. Then we study the sand dunes characteristic (wavelength and celerity) in the Rio Parapeti. Because of the shallow-flow configuration the dominant dune length can be easily extracted from satellite images taken at various dates. We also show that it is more than likely that dune movement can be followed by the simple deployment of a pressure probe into the water under stable discharge condition, even if further data and investigation are necessary to confirm this

Improving Prolog Programs: Refactoring for Prolog

Refactoring is an established technique from the OO-community to restructure code: it aims at improving software readability, maintainability and extensibility. Although refactoring is not tied to the OO-paradigm in particular, its ideas have not been applied to Logic Programming until now. This paper applies the ideas of refactoring to Prolog programs. A catalogue is presented listing refactorings classified according to scope. Some of the refactorings have been adapted from the OO-paradigm, while others have been specifically designed for Prolog. Also the discrepancy between intended and operational semantics in Prolog is addressed by some of the refactorings. In addition, ViPReSS, a semi-automatic refactoring browser, is discussed and the experience with applying \vipress to a large Prolog legacy system is reported. Our main conclusion is that refactoring is not only a viable technique in Prolog but also a rather desirable one.Comment: To appear in ICLP 200

Quantitative Stability of Linear Infinite Inequality Systems under Block Perturbations with Applications to Convex Systems

The original motivation for this paper was to provide an efficient quantitative analysis of convex infinite (or semi-infinite) inequality systems whose decision variables run over general infinite-dimensional (resp. finite-dimensional) Banach spaces and that are indexed by an arbitrary fixed set $J$. Parameter perturbations on the right-hand side of the inequalities are required to be merely bounded, and thus the natural parameter space is $l_{\infty}(J)$. Our basic strategy consists of linearizing the parameterized convex system via splitting convex inequalities into linear ones by using the Fenchel-Legendre conjugate. This approach yields that arbitrary bounded right-hand side perturbations of the convex system turn on constant-by-blocks perturbations in the linearized system. Based on advanced variational analysis, we derive a precise formula for computing the exact Lipschitzian bound of the feasible solution map of block-perturbed linear systems, which involves only the system's data, and then show that this exact bound agrees with the coderivative norm of the aforementioned mapping. In this way we extend to the convex setting the results of [3] developed for arbitrary perturbations with no block structure in the linear framework under the boundedness assumption on the system's coefficients. The latter boundedness assumption is removed in this paper when the decision space is reflexive. The last section provides the aimed application to the convex case

On Distributive Subalgebras of Qualitative Spatial and Temporal Calculi

Qualitative calculi play a central role in representing and reasoning about qualitative spatial and temporal knowledge. This paper studies distributive subalgebras of qualitative calculi, which are subalgebras in which (weak) composition distributives over nonempty intersections. It has been proven for RCC5 and RCC8 that path consistent constraint network over a distributive subalgebra is always minimal and globally consistent (in the sense of strong $n$-consistency) in a qualitative sense. The well-known subclass of convex interval relations provides one such an example of distributive subalgebras. This paper first gives a characterisation of distributive subalgebras, which states that the intersection of a set of $n\geq 3$ relations in the subalgebra is nonempty if and only if the intersection of every two of these relations is nonempty. We further compute and generate all maximal distributive subalgebras for Point Algebra, Interval Algebra, RCC5 and RCC8, Cardinal Relation Algebra, and Rectangle Algebra. Lastly, we establish two nice properties which will play an important role in efficient reasoning with constraint networks involving a large number of variables.Comment: Adding proof of Theorem 2 to appendi

Magneto-shear modes and a.c. dissipation in a two-dimensional Wigner crystal

The a.c. response of an unpinned and finite 2D Wigner crystal to electric fields at an angular frequency $\omega$ has been calculated in the dissipative limit, $\omega \tau \ll 1$, where $\tau ^{-1}$ is the scattering rate. For electrons screened by parallel electrodes, in zero magnetic field the long-wavelength excitations are a diffusive longitudinal transmission line mode and a diffusive shear mode. A magnetic field couples these modes together to form two new magneto-shear modes. The dimensionless coupling parameter $\beta =2(c_{t}/c_{l})|\sigma_{xy}/\sigma_{xx}|$ where $c_{t}$ and $c_{l}$ are the speeds of transverse and longitudinal sound in the collisionless limit and $\sigma_{xy}$ and $\sigma_{xx}$ are the tensor components of the magnetoconductivity. For $\beta \geqslant 1$, both the coupled modes contribute to the response of 2D electrons in a Corbino disk measurement of magnetoconductivity. For $\beta \gg 1$, the electron crystal rotates rigidly in a magnetic field. In general, both the amplitude and phase of the measured a.c. currents are changed by the shear modulus. In principle, both the magnetoconductivity and the shear modulus can be measured simultaneously.Comment: REVTeX, 7 pp., 4 eps figure

From nonwetting to prewetting: the asymptotic behavior of 4He drops on alkali substrates

We investigate the spreading of 4He droplets on alkali surfaces at zero temperature, within the frame of Finite Range Density Functional theory. The equilibrium configurations of several 4He_N clusters and their asymptotic trend with increasing particle number N, which can be traced to the wetting behavior of the quantum fluid, are examined for nanoscopic droplets. We discuss the size effects, inferring that the asymptotic properties of large droplets correspond to those of the prewetting film

Low-Temperature Mobility of Surface Electrons and Ripplon-Phonon Interaction in Liquid Helium

The low-temperature dc mobility of the two-dimensional electron system localized above the surface of superfluid helium is determined by the slowest stage of the longitudinal momentum transfer to the bulk liquid, namely, by the interaction of surface and volume excitations of liquid helium, which rapidly decreases with temperature. Thus, the temperature dependence of the low-frequency mobility is \mu_{dc} = 8.4x10^{-11}n_e T^{-20/3} cm^4 K^{20/3}/(V s), where n_e is the surface electron density. The relation T^{20/3}E_\perp^{-3} << 2x10^{-7} between the pressing electric field (in kV/cm) and temperature (in K) and the value \omega < 10^8 T^5 K^{-5}s^{-1} of the driving-field frequency have been obtained, at which the above effect can be observed. In particular, E_\perp = 1 kV/cm corresponds to T < 70 mK and \omega/2\pi < 30 Hz.Comment: 4 pages, 1 figur

Shaping bursting by electrical coupling and noise

Gap-junctional coupling is an important way of communication between neurons and other excitable cells. Strong electrical coupling synchronizes activity across cell ensembles. Surprisingly, in the presence of noise synchronous oscillations generated by an electrically coupled network may differ qualitatively from the oscillations produced by uncoupled individual cells forming the network. A prominent example of such behavior is the synchronized bursting in islets of Langerhans formed by pancreatic \beta-cells, which in isolation are known to exhibit irregular spiking. At the heart of this intriguing phenomenon lies denoising, a remarkable ability of electrical coupling to diminish the effects of noise acting on individual cells. In this paper, we derive quantitative estimates characterizing denoising in electrically coupled networks of conductance-based models of square wave bursting cells. Our analysis reveals the interplay of the intrinsic properties of the individual cells and network topology and their respective contributions to this important effect. In particular, we show that networks on graphs with large algebraic connectivity or small total effective resistance are better equipped for implementing denoising. As a by-product of the analysis of denoising, we analytically estimate the rate with which trajectories converge to the synchronization subspace and the stability of the latter to random perturbations. These estimates reveal the role of the network topology in synchronization. The analysis is complemented by numerical simulations of electrically coupled conductance-based networks. Taken together, these results explain the mechanisms underlying synchronization and denoising in an important class of biological models