Page 1

Displaying 1 – 15 of 15

Showing per page

A comparison of deterministic and Bayesian inverse with application in micromechanics

Radim Blaheta, Michal Béreš, Simona Domesová, Pengzhi Pan (2018)

Applications of Mathematics

The paper deals with formulation and numerical solution of problems of identification of material parameters for continuum mechanics problems in domains with heterogeneous microstructure. Due to a restricted number of measurements of quantities related to physical processes, we assume additional information about the microstructure geometry provided by CT scan or similar analysis. The inverse problems use output least squares cost functionals with values obtained from averages of state problem quantities...

A linear mixed finite element scheme for a nematic Ericksen–Leslie liquid crystal model

F. M. Guillén-González, J. V. Gutiérrez-Santacreu (2013)

ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique

In this work we study a fully discrete mixed scheme, based on continuous finite elements in space and a linear semi-implicit first-order integration in time, approximating an Ericksen–Leslie nematic liquid crystal model by means of a Ginzburg–Landau penalized problem. Conditional stability of this scheme is proved via a discrete version of the energy law satisfied by the continuous problem, and conditional convergence towards generalized Young measure-valued solutions to the Ericksen–Leslie problem...

Coarse-graining schemes and a posteriori error estimates for stochastic lattice systems

Markos A. Katsoulakis, Petr Plecháč, Luc Rey-Bellet, Dimitrios K. Tsagkarogiannis (2007)

ESAIM: Mathematical Modelling and Numerical Analysis

The primary objective of this work is to develop coarse-graining schemes for stochastic many-body microscopic models and quantify their effectiveness in terms of a priori and a posteriori error analysis. In this paper we focus on stochastic lattice systems of interacting particles at equilibrium. The proposed algorithms are derived from an initial coarse-grained approximation that is directly computable by Monte Carlo simulations, and the corresponding numerical error is calculated using the...

Combining stochastic and deterministic approaches within high efficiency molecular simulations

Bruno Escribano, Elena Akhmatskaya, Jon Mujika (2013)

Open Mathematics

Generalized Shadow Hybrid Monte Carlo (GSHMC) is a method for molecular simulations that rigorously alternates Monte Carlo sampling from a canonical ensemble with integration of trajectories using Molecular Dynamics (MD). While conventional hybrid Monte Carlo methods completely re-sample particle’s velocities between MD trajectories, our method suggests a partial velocity update procedure which keeps a part of the dynamic information throughout the simulation. We use shadow (modified) Hamiltonians,...

Convergence of gradient-based algorithms for the Hartree-Fock equations

Antoine Levitt (2012)

ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique

The numerical solution of the Hartree-Fock equations is a central problem in quantum chemistry for which numerous algorithms exist. Attempts to justify these algorithms mathematically have been made, notably in [E. Cancès and C. Le Bris, Math. Mod. Numer. Anal. 34 (2000) 749–774], but, to our knowledge, no complete convergence proof has been published, except for the large-Z result of [M. Griesemer and F. Hantsch, Arch. Rational Mech. Anal. (2011) 170]. In this paper, we prove the convergence of...

Convergence of gradient-based algorithms for the Hartree-Fock equations

Antoine Levitt (2012)

ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique

The numerical solution of the Hartree-Fock equations is a central problem in quantum chemistry for which numerous algorithms exist. Attempts to justify these algorithms mathematically have been made, notably in [E. Cancès and C. Le Bris, Math. Mod. Numer. Anal. 34 (2000) 749–774], but, to our knowledge, no complete convergence proof has been published, except for the large-Z result of [M. Griesemer and F. Hantsch, Arch. Rational Mech. Anal. (2011) 170]. In this paper, we prove the convergence of...

Convergence of gradient-based algorithms for the Hartree-Fock equations∗

Antoine Levitt (2012)

ESAIM: Mathematical Modelling and Numerical Analysis

The numerical solution of the Hartree-Fock equations is a central problem in quantum chemistry for which numerous algorithms exist. Attempts to justify these algorithms mathematically have been made, notably in [E. Cancès and C. Le Bris, Math. Mod. Numer. Anal. 34 (2000) 749–774], but, to our knowledge, no complete convergence proof has been published, except for the large-Z result of [M. Griesemer and F. Hantsch, Arch. Rational Mech. Anal. (2011) ...

Finite element approximation of finitely extensible nonlinear elastic dumbbell models for dilute polymers

John W. Barrett, Endre Süli (2012)

ESAIM: Mathematical Modelling and Numerical Analysis

We construct a Galerkin finite element method for the numerical approximation of weak solutions to a general class of coupled FENE-type finitely extensible nonlinear elastic dumbbell models that arise from the kinetic theory of dilute solutions of polymeric liquids with noninteracting polymer chains. The class of models involves the unsteady incompressible Navier–Stokes equations in a bounded domain Ω ⊂ ℝd, d = 2 or 3, for the velocity...

Finite element approximation of finitely extensible nonlinear elastic dumbbell models for dilute polymers

John W. Barrett, Endre Süli (2012)

ESAIM: Mathematical Modelling and Numerical Analysis

We construct a Galerkin finite element method for the numerical approximation of weak solutions to a general class of coupled FENE-type finitely extensible nonlinear elastic dumbbell models that arise from the kinetic theory of dilute solutions of polymeric liquids with noninteracting polymer chains. The class of models involves the unsteady incompressible Navier–Stokes equations in a bounded domain Ω ⊂ ℝd, d = 2 or 3, for the velocity...

Numerical analysis of the planewave discretization of some orbital-free and Kohn-Sham models

Eric Cancès, Rachida Chakir, Yvon Maday (2012)

ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique

In this article, we provide a priorierror estimates for the spectral and pseudospectral Fourier (also called planewave) discretizations of the periodic Thomas-Fermi-von Weizsäcker (TFW) model and for the spectral discretization of the periodic Kohn-Sham model, within the local density approximation (LDA). These models allow to compute approximations of the electronic ground state energy and density of molecular systems in the condensed phase. The TFW model is strictly convex with respect to the...

Numerical analysis of the planewave discretization of some orbital-free and Kohn-Sham models

Eric Cancès, Rachida Chakir, Yvon Maday (2011)

ESAIM: Mathematical Modelling and Numerical Analysis

In this article, we provide a priori error estimates for the spectral and pseudospectral Fourier (also called planewave) discretizations of the periodic Thomas-Fermi-von Weizsäcker (TFW) model and for the spectral discretization of the periodic Kohn-Sham model, within the local density approximation (LDA). These models allow to compute approximations of the electronic ground state energy and density of molecular systems in the condensed phase. The TFW model is strictly convex with respect to the...

Role of Molecular Chaos in Granular Fluctuating Hydrodynamics

G. Costantini, A. Puglisi (2011)

Mathematical Modelling of Natural Phenomena

We perform a numerical study of the fluctuations of the rescaled hydrodynamic transverse velocity field during the cooling state of a homogeneous granular gas. We are interested in the role of Molecular Chaos for the amplitude of the hydrodynamic noise and its relaxation in time. For this purpose we compare the results of Molecular Dynamics (MD, deterministic dynamics) with those from Direct Simulation Monte Carlo (DSMC, random process), where Molecular...

The continuous Coupled Cluster formulation for the electronic Schrödinger equation

Thorsten Rohwedder (2013)

ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique

Nowadays, the Coupled Cluster (CC) method is the probably most widely used high precision method for the solution of the main equation of electronic structure calculation, the stationary electronic Schrödinger equation. Traditionally, the equations of CC are formulated as a nonlinear approximation of a Galerkin solution of the electronic Schrödinger equation, i.e. within a given discrete subspace. Unfortunately, this concept prohibits the direct application of concepts of nonlinear numerical analysis...

Currently displaying 1 – 15 of 15

Page 1