Analysis, by simulation, of the impact of a technical default of a payment system participant.

Payment systems play a very important role in ensuring the safe and efficient transfer of deposits and financial instruments. Consequently, the failure of these systems may have a destabilising impact. Business continuity plans have thus been developed to ensure their robustness. However, their smooth functioning is also contingent on the capacity of participants to submit their payment orders. The Banque de France, in its role of overseer of the French payment systems, conducted a study with a view to enhancing its understanding of the consequences and the impact of the technical default of a participant in such systems. This study, carried out using a simulator of the functioning of the Paris Net Settlement (PNS) large-value payment system, operated by the CRI (Centrale des Règlements Interbancaires), shows that the technical default of a participant in this system has negative consequences on the smooth running of the system. Indeed, a situation in which a major participant, in the wake of a technical incident, is unable to submit its payment orders in a normal fashion to its counterparties in PNS, could further exacerbate congestion in the system and result in almost 10% of payments being rejected among non-defaulting participants. The consequences of a technical default could nevertheless be greatly reduced if the participants set their bilateral sender limits at a lower level than that currently observed and if they reacted rapidly to information indicating a technical default by reducing their bilateral limits with the defaulting participant (defaulter).

Holonomic constraints : an analytical result

Systems subjected to holonomic constraints follow quite complicated dynamics that could not be described easily with Hamiltonian or Lagrangian dynamics. The influence of holonomic constraints in equations of motions is taken into account by using Lagrange multipliers. Finding the value of the Lagrange multipliers allows to compute the forces induced by the constraints and therefore, to integrate the equations of motions of the system. Computing analytically the Lagrange multipliers for a constrained system may be a difficult task that is depending on the complexity of systems. For complex systems, it is most of the time impossible to achieve. In computer simulations, some algorithms using iterative procedures estimate numerically Lagrange multipliers or constraint forces by correcting the unconstrained trajectory. In this work, we provide an analytical computation of the Lagrange multipliers for a set of linear holonomic constraints with an arbitrary number of bonds of constant length. In the appendix of the paper, one would find explicit formulas for Lagrange multipliers for systems having 1, 2, 3, 4 and 5 bonds of constant length, linearly connected.Comment: 13 pages, no figures. To appear in J. Phys. A : Math. The

A 3D beam element analysis for R/C structural walls

To analyse the real 3D functioning of a structure under seismic loading the dialogue between tests and numerical simulations is needed. Within the framework of the TMR-ICONS research program, dynamic and cyclic tests on U-shaped shear walls have been performed at CEA Saclay and JRC Ispra respectively. More recently, for the French program ìCAMUS 2000î, shaking table tests have been performed on reinforced concrete structural walls. In order to simulate these tests, 3D multi-fiber beam elements are used. Comparison with the experimental results shows the well matching and the limitations of the approach

Yukawa potentials in systems with partial periodic boundary conditions I : Ewald sums for quasi-two dimensional systems

Yukawa potentials are often used as effective potentials for systems as colloids, plasmas, etc. When the Debye screening length is large, the Yukawa potential tends to the non-screened Coulomb potential ; in this small screening limit, or Coulomb limit, the potential is long ranged. As it is well known in computer simulation, a simple truncation of the long ranged potential and the minimum image convention are insufficient to obtain accurate numerical data on systems. The Ewald method for bulk systems, i.e. with periodic boundary conditions in all three directions of the space, has already been derived for Yukawa potential [cf. Y., Rosenfeld, {\it Mol. Phys.}, \bm{88}, 1357, (1996) and G., Salin and J.-M., Caillol, {\it J. Chem. Phys.}, \bm{113}, 10459, (2000)], but for systems with partial periodic boundary conditions, the Ewald sums have only recently been obtained [M., Mazars, {\it J. Chem. Phys.}, {\bf 126}, 056101 (2007)]. In this paper, we provide a closed derivation of the Ewald sums for Yukawa potentials in systems with periodic boundary conditions in only two directions and for any value of the Debye length. A special attention is paid to the Coulomb limit and its relation with the electroneutrality of systems.Comment: 40 pages, 5 figures and 4 table

3D Impacts Modeling of the Magnetic Pulse Welding Process and Comparison to Experimental Data

Magnetic Pulse Welding (MPW) is a solid state (cold) welding process known to present several advantages. When properly designed, such an assembly is stronger than the weakest base material even for multi-material joining. These high quality welds are due to an almost inexistent Heat Affected Zone which is not the case with fusion welding solutions. Another advantage is a welding time that is under a millisecond. In order to define the MPW parameters (mainly geometry, current and frequency), recent developments have made it possible to adapt welding windows from the Explosive Welding (EXW) for use in MPW. Until now, these welding windows have been simulated only in 2D geometries showing how the impact angle and the radial velocities progress in a welding window. The aim of this paper is to present our most recent development, which builds on this analysis to develop a 3D model in order to deal for example with local planar MPW. Simulation results will be presented and then compared to experimental data for a multimaterial join case

Numerical modelling for earthquake engineering: the case of lightly RC structural walls

Different types of numerical models exist to describe the non‐linear behaviour of reinforced concrete structures. Based on the level of discretization they are often classified as refined or simplified ones. The efficiency of two simplified models using beam elements and damage mechanics in describing the global and local behaviour of lightly reinforced concrete structural walls subjected to seismic loadings is investigated in this paper. The first model uses an implicit and the second an explicit numerical scheme. For each case, the results of the CAMUS 2000 experimental programme are used to validate the approaches

Structure and thermodynamics of a ferrofluid bilayer

We present extensive Monte Carlo simulations for the thermodynamic and structural properties of a planar bilayer of dipolar hard spheres for a wide range of densities, dipole moments and layer separations. Expressions for the stress and pressure tensors of the bilayer system are derived. For all thermodynamic states considered the interlayer energy is shown to be attractive and much smaller than the intralayer contribution to the energy. It vanishes at layer separations of the order of two hard sphere diameters. The normal pressure is negative and decays as a function of layer separation $h$ as $-1/h^5$. Intralayer and interlayer pair distribution functions and angular correlation functions are presented. Despite the weak interlayer energy strong positional and orientational correlations exist between particles in the two layers.Comment: 45 pages, 4 Tables, 9 Figure

Linearly scaling direct method for accurately inverting sparse banded matrices

In many problems in Computational Physics and Chemistry, one finds a special kind of sparse matrices, termed "banded matrices". These matrices, which are defined as having non-zero entries only within a given distance from the main diagonal, need often to be inverted in order to solve the associated linear system of equations. In this work, we introduce a new O(n) algorithm for solving such a system, being n X n the size of the matrix. We produce the analytical recursive expressions that allow to directly obtain the solution, as well as the pseudocode for its computer implementation. Moreover, we review the different options for possibly parallelizing the method, we describe the extension to deal with matrices that are banded plus a small number of non-zero entries outside the band, and we use the same ideas to produce a method for obtaining the full inverse matrix. Finally, we show that the New Algorithm is competitive, both in accuracy and in numerical efficiency, when compared to a standard method based in Gaussian elimination. We do this using sets of large random banded matrices, as well as the ones that appear when one tries to solve the 1D Poisson equation by finite differences.Comment: 24 pages, 5 figures, submitted to J. Comp. Phy

Efficient and reliable nonlocal damage models

We present an efficient and reliable approach for the numerical modelling of failure with nonlocal damage models. The two major numerical challenges––the strongly nonlinear, highly localized and parameter-dependent structural response of quasi-brittle materials, and the interaction between nonadjacent finite elements associated to nonlocality––are addressed in detail. Reliability of the numerical results is ensured by an h-adaptive strategy based on error estimation. We use a residual-type error estimator for nonlinear FE analysis based on local computations, which, at the same time, accounts for the nonlocality of the damage model. Efficiency is achieved by a proper combination of load-stepping control technique and iterative solver for the nonlinear equilibrium equations. A major issue is the computation of the consistent tangent matrix, which is nontrivial due to nonlocal interaction between Gauss points. With computational efficiency in mind, we also present a new nonlocal damage model based on the nonlocal average of displacements. For this new model, the consistent tangent matrix is considerably simpler to compute than for current models. The various ideas discussed in the paper are illustrated by means of three application examples: the uniaxial tension test, the three-point bending test and the single-edge notched beam test.Peer ReviewedPostprint (author’s final draft