@article{CotterHamPainetal.2009, author = {Cotter, Colin J. and Ham, David A. and Pain, Christopher C. and Reich, Sebastian}, title = {LBB stability of a mixed Galerkin finite element pair for fluid flow simulations}, issn = {0021-9991}, doi = {10.1016/j.jcp.2008.09.014}, year = {2009}, abstract = {We introduce a new mixed finite element for solving the 2- and 3-dimensional wave equations and equations of incompressible flow. The element, which we refer to as P1(D)-P2, uses discontinuous piecewise linear functions for velocity and continuous piecewise quadratic functions for pressure. The aim of introducing the mixed formulation is to produce a new flexible element choice for triangular and tetrahedral meshes which satisfies the LBB stability condition and hence has no spurious zero-energy modes. The advantage of this particular element choice is that the mass matrix for velocity is block diagonal so it can be trivially inverted; it also allows the order of the pressure to be increased to quadratic whilst maintaining LBB stability which has benefits in geophysical applications with Coriolis forces. We give a normal mode analysis of the semi-discrete wave equation in one dimension which shows that the element pair is stable, and demonstrate that the element is stable with numerical integrations of the wave equation in two dimensions, an analysis of the resultant discrete Laplace operator in two and three dimensions on various meshes which shows that the element pair does not have any spurious modes. We provide convergence tests for the element pair which confirm that the element is stable since the convergence rate of the numerical solution is quadratic.}, language = {en} } @article{AkhmatskayaBouRabeeReich2009, author = {Akhmatskaya, Elena and Bou-Rabee, Nawaf and Reich, Sebastian}, title = {Erratum to "A comparison of generalized hybrid Monte Carlo methods with and without momentum flip" [J. Comput. Phys. 228 (2009), S. 2256 - 2265]}, issn = {0021-9991}, doi = {10.1016/j.jcp.2009.06.039}, year = {2009}, abstract = {The generalized hybrid Monte Carlo (GHMC) method combines Metropolis corrected constant energy simulations with a partial random refreshment step in the particle momenta. The standard detailed balance condition requires that momenta are negated upon rejection of a molecular dynamics proposal step. The implication is a trajectory reversal upon rejection, which is undesirable when interpreting GHMC as thermostated molecular dynamics. We show that a modified detailed balance condition can be used to implement GHMC without momentum flips. The same modification can be applied to the generalized shadow hybrid Monte Carlo (GSHMC) method. Numerical results indicate that GHMC/GSHMC implementations with momentum flip display a favorable behavior in terms of sampling efficiency, i.e., the traditional GHMC/GSHMC implementations with momentum flip got the advantage of a higher acceptance rate and faster decorrelation of Monte Carlo samples. The difference is more pronounced for GHMC. We also numerically investigate the behavior of the GHMC method as a Langevin-type thermostat. We find that the GHMC method without momentum flip interferes less with the underlying stochastic molecular dynamics in terms of autocorrelation functions and it to be preferred over the GHMC method with momentum flip. The same finding applies to GSHMC.}, language = {en} } @article{BergemannGottwaldReich2009, author = {Bergemann, Kay and Gottwald, Georg and Reich, Sebastian}, title = {Ensemble propagation and continuous matrix factorization algorithms}, issn = {0035-9009}, doi = {10.1002/qj.457}, year = {2009}, abstract = {We consider the problem of propagating an ensemble of solutions and its characterization in terms of its mean and covariance matrix. We propose differential equations that lead to a continuous matrix factorization of the ensemble into a generalized singular value decomposition (SVD). The continuous factorization is applied to ensemble propagation under periodic rescaling (ensemble breeding) and under periodic Kalman analysis steps (ensemble Kalman filter). We also use the continuous matrix factorization to perform a re-orthogonalization of the ensemble after each time-step and apply the resulting modified ensemble propagation algorithm to the ensemble Kalman filter. Results from the Lorenz-96 model indicate that the re-orthogonalization of the ensembles leads to improved filter performance.}, language = {en} } @article{LeimkuhlerReich2009, author = {Leimkuhler, Benedict and Reich, Sebastian}, title = {A metropolis adjusted Nos{\´e}-Hoover thermostat}, issn = {0764-583X}, doi = {10.1051/M2an/2009023}, year = {2009}, abstract = {We present a Monte Carlo technique for sampling from the canonical distribution in molecular dynamics. The method is built upon the Nose-Hoover constant temperature formulation and the generalized hybrid Monte Carlo method. In contrast to standard hybrid Monte Carlo methods only the thermostat degree of freedom is stochastically resampled during a Monte Carlo step.}, language = {en} } @article{AkhmatskayaBouRabeeReich2009, author = {Akhmatskaya, Elena and Bou-Rabee, Nawaf and Reich, Sebastian}, title = {A comparison of generalized hybrid Monte Carlo methods with and without momentum flip}, issn = {0021-9991}, doi = {10.1016/j.jcp.2008.12.014}, year = {2009}, abstract = {The generalized hybrid Monte Carlo (GHMC) method combines Metropolis corrected constant energy simulations with a partial random refreshment step in the particle momenta. The standard detailed balance condition requires that momenta are negated upon rejection of a molecular dynamics proposal step. The implication is a trajectory reversal upon rejection, which is undesirable when interpreting GHMC as thermostated molecular dynamics. We show that a modified detailed balance condition can be used to implement GHMC without momentum flips. The same modification can be applied to the generalized shadow hybrid Monte Carlo (GSHMC) method. Numerical results indicate that GHMC/GSHMC implementations with momentum flip display a favorable behavior in terms of sampling efficiency, i.e., the traditional GHMC/GSHMC implementations with momentum flip got the advantage of a higher acceptance rate and faster decorrelation of Monte Carlo samples. The difference is more pronounced for GHMC. We also numerically investigate the behavior of the GHMC method as a Langevin-type thermostat. We find that the GHMC method without momentum flip interferes less with the underlying stochastic molecular dynamics in terms of autocorrelation functions and it to be preferred over the GHMC method with momentum flip. The same finding applies to GSHMC.}, language = {en} }