Search results for "Numerical"

Energy-stable linear schemes for polymer-solvent phase field models

2017

We present new linear energy-stable numerical schemes for numerical simulation of complex polymer-solvent mixtures. The mathematical model proposed by Zhou, Zhang and E (Physical Review E 73, 2006) consists of the Cahn-Hilliard equation which describes dynamics of the interface that separates polymer and solvent and the Oldroyd-B equations for the hydrodynamics of polymeric mixtures. The model is thermodynamically consistent and dissipates free energy. Our main goal in this paper is to derive numerical schemes for the polymer-solvent mixture model that are energy dissipative and efficient in time. To this end we will propose several problem-suited time discretizations yielding linear scheme…

1982

The molecular weight distribution (MWD) of a high polymer is calculated from a weakly perturbed Zimm-plot of the classical light scattering on dilute solutions of Gaussian polymer coils (theta state). A typical Zimm-plot is simulated corresponding to the measurements of high accuracy as would be obtained by using the laser photometer described by Hack and Meyerhoff. The accuracy as published by these authors for small dissymmetries is used. Two inversion procedures are described in detail and tested, both of them avoiding the use of an empirical formula for the Laplace image of the calculated MWD. A complete numerical analysis of the results is given. The results are compared with those of …

GPU accelerated Monte Carlo simulations of lattice spin models

2011

We consider Monte Carlo simulations of classical spin models of statistical mechanics using the massively parallel architecture provided by graphics processing units (GPUs). We discuss simulations of models with discrete and continuous variables, and using an array of algorithms ranging from single-spin flip Metropolis updates over cluster algorithms to multicanonical and Wang-Landau techniques to judge the scope and limitations of GPU accelerated computation in this field. For most simulations discussed, we find significant speed-ups by two to three orders of magnitude as compared to single-threaded CPU implementations.

On the Complexity of the Bernstein Combinatorial Problem

2012

International audience

Mechanical properties of carbon nanotube reinforced nanocomposites

2013

One way to take advantage of the marvelous properties of the carbon nanotubes consists in incorporating them into a matrix to build composite materials. The best candidates for this task are undoubtedly polymers thanks to their strength, toughness, low weight, and easy processing. A mixed model, numerical-analytical, is presented that allows one to predict the elastic properties of carbon nanotube (CNT)/polymer composites containing a random distribution of CNTs, while taking account of the curvature that they show when immersed in the polymer. A three-dimensional model, using a single CNT immersed in an infinite matrix, allows us to numerically evaluate the concentration tensor by evaluati…

Computable majorants of the limit load in Hencky’s plasticity problems

2018

Abstract We propose a new method for analyzing the limit (safe) load of elastoplastic media governed by the Hencky plasticity law and deduce fully computable bounds of this load. The main idea of the method is based on a combination of kinematic approach and new estimates of the distance to the set of divergence free fields. We show that two sided bounds of the limit load are sharp and the computational efficiency of the method is confirmed by numerical experiments.

Sur la prise en compte de quelques conditions aux limites avec la méthode des éléments finis

2017

Ce cours a pour principal sujet la discrétisation de certaines conditions limites sur le bord d'un domaine à l'aide de la méthode des éléments finis. Nous étudions en particulier des problèmes elliptiques avec conditions de Dirichlet non-homogènes, puis avec conditions de Signorini. Les conditions de Signorini sont des conditions de contact unilatéral entre un solide et un support rigide. Différentes méthodes sont passées en revue, et nous détaillons l'analyse pour les méthodes classiques et celle de Nitsche.

Attractor as a convex combination of a set of attractors

2021

This paper presents an effective approach to constructing numerical attractors of a general class of continuous homogenous dynamical systems: decomposing an attractor as a convex combination of a set of other existing attractors. For this purpose, the convergent Parameter Switching (PS) numerical method is used to integrate the underlying dynamical system. The method is built on a convergent fixed step-size numerical method for ODEs. The paper shows that the PS algorithm, incorporating two binary operations, can be used to approximate any numerical attractor via a convex combination of some existing attractors. Several examples are presented to show the effectiveness of the proposed method.…

New Indicators of Approximation Errors for Problems in Continuum Mechanics

2010

In this paper we present a new error indicator for approximate solutions of elliptic problems. We discuss error indication with the paradigm of the diffusion problem, however the techniques are easily adaptable to more complicated elliptic problems, for example to linear elasticity, viscous flow models and electromagnetic models. The proposed indicator does not contain mesh dependent constants and it admits parallelization. nonPeerReviewed

Predictive control of convex polyhedron LPV systems with Markov jumping parameters

2012

The problem of receding horizon predictive control of stochastic linear parameter varying systems is discussed. First, constant coefficient matrices are obtained at each vertex in the interior of linear parameter varying system, and then, by considering semi-definite programming constraints, weight coefficients between each vertex are calculated, and the equal coefficients matrices for the time variable system are obtained. Second, in the given receding horizon, for each mode sequence of the stochastic convex polyhedron linear parameter varying systems, the optimal control input sequences are designed in order to make the states into a terminal invariant set. Outside of the receding horizon…