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

Non-meanfield deterministic limits in chemical reaction kinetics far from equilibrium

A general mechanism is proposed by which small intrinsic fluctuations in a system far from equilibrium can result in nearly deterministic dynamical behaviors which are markedly distinct from those realized in the meanfield limit. The mechanism is demonstrated for the kinetic Monte-Carlo version of the Schnakenberg reaction where we identified a scaling limit in which the global deterministic bifurcation picture is fundamentally altered by fluctuations. Numerical simulations of the model are found to be in quantitative agreement with theoretical predictions.Comment: 4 pages, 4 figures (submitted to Phys. Rev. Lett.

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

Numerical Evaluation of the Accuracy and Stability Properties of High-order Direct Stokes Solvers with or without Temporal Splitting

The temporal stability and effective order of two different direct high-order Stokes solvers are examined. Both solvers start from the primitive variables formulation of the Stokes problem, but are distinct by the numerical uncoupling they apply on the Stokes operator. One of these solvers introduces an intermediate divergence free velocity for performing a temporal splitting (J. Comput. Phys. [1991] 97, 414-443) while the other treats the whole Stokes problem through the evaluation of a divergence free acceleration field (C.R. Acad. Sci. Paris [1994] 319 Serie I, 1455-1461; SIAM J. Scient. Comput. [2000] 22(4), 1386-1410). The second uncoupling is known to be consistent with the harmonicity of the pressure field (SIAM J. Scient. Comput. [2000] 22(4), 1386-1410). Both solvers proceed by two steps, a pressure evaluation based on an extrapolated in time (of theoretical order Je) Neumann condition, and a time implicit (of theoretical order Ji) diffusion step for the final velocity. These solvers are implemented with a Chebyshev mono-domain and a Legendre spectral element collocation schemes. The numerical stability of these four options is investigated for the sixteen combinations of (Je,Ji), 1 ≤ Je, Ji ≤

Large-eddy simulation of the lid-driven cubic cavity flow by the spectral element method

This paper presents the large-eddy simulation of the lid-driven cubic cavity flow by the spectral element method (SEM) using the dynamic model. Two spectral filtering techniques suitable for these simulations have been implemented. Numerical results for Reynolds number $\text{Re}=12'000$ are showing very good agreement with other experimental and DNS results found in the literature

Freezing-induced self-assembly of amphiphilic molecules

The self-assembly of amphiphilic molecules usually takes place in a liquid phase, near room temperature. Here, using small angle X-ray scattering (SAXS) experiments performed in real time, we show that freezing of aqueous solutions of copolymer amphiphilic molecules can induce self-assembly below 0{\deg}C.Comment: 10 pages, 3 figure

Three-dimensional flow instability in a lid-driven isosceles triangular cavity

Linear three-dimensional modal instability of steady laminar two-dimensional states developing in a lid-driven cavity of isosceles triangular cross-section is investigated theoretically and experimentally for the case in which the equal sides form a rectangular corner. An asymmetric steady two-dimensional motion is driven by the steady motion of one of the equal sides. If the side moves away from the rectangular corner, a stationary three-dimensional instability is found. If the motion is directed towards the corner, the instability is oscillatory. The respective critical Reynolds numbers are identified both theoretically and experimentally. The neutral curves pertinent to the two configurations and the properties of the respective leading eigenmodes are documented and analogies to instabilities in rectangular lid-driven cavities are discussed

Adaptive mesh refinement with spectral accuracy for magnetohydrodynamics in two space dimensions

We examine the effect of accuracy of high-order spectral element methods, with or without adaptive mesh refinement (AMR), in the context of a classical configuration of magnetic reconnection in two space dimensions, the so-called Orszag-Tang vortex made up of a magnetic X-point centered on a stagnation point of the velocity. A recently developed spectral-element adaptive refinement incompressible magnetohydrodynamic (MHD) code is applied to simulate this problem. The MHD solver is explicit, and uses the Elsasser formulation on high-order elements. It automatically takes advantage of the adaptive grid mechanics that have been described elsewhere in the fluid context [Rosenberg, Fournier, Fischer, Pouquet, J. Comp. Phys. 215, 59-80 (2006)]; the code allows both statically refined and dynamically refined grids. Tests of the algorithm using analytic solutions are described, and comparisons of the Orszag-Tang solutions with pseudo-spectral computations are performed. We demonstrate for moderate Reynolds numbers that the algorithms using both static and refined grids reproduce the pseudo--spectral solutions quite well. We show that low-order truncation--even with a comparable number of global degrees of freedom--fails to correctly model some strong (sup--norm) quantities in this problem, even though it satisfies adequately the weak (integrated) balance diagnostics.Comment: 19 pages, 10 figures, 1 table. Submitted to New Journal of Physic

Progress report and first operation of the GANIL injector

http://accelconf.web.cern.ch/AccelConf/c81/papers/abp-07.pdfInternational audienc