Computational Physics (physics.comp-ph)

  • PDF
    The present paper deals with the numerical solution of the incompressible Navier-Stokes equations using high-order discontinuous Galerkin (DG) methods for discretization in space. For DG methods applied to the dual splitting projection method, instabilities have recently been reported that occur for coarse spatial resolutions and small time step sizes. By means of numerical investigation we give evidence that these instabilities are related to the discontinuous Galerkin formulation of the velocity divergence term and the pressure gradient term that couple velocity and pressure. Integration by parts of these terms with a suitable definition of boundary conditions is required in order to obtain a stable and robust method. Since the intermediate velocity field does not fulfill the boundary conditions prescribed for the velocity, a consistent boundary condition is derived from the convective step of the dual splitting scheme to ensure high-order accuracy with respect to the temporal discretization. This new formulation is stable in the limit of small time steps for both equal-order and mixed-order polynomial approximations. Although the dual splitting scheme itself includes inf-sup stabilizing contributions, we demonstrate that spurious pressure oscillations appear for equal-order polynomials and small time steps highlighting the necessity to consider inf-sup stability explicitly.
  • PDF
    Using the method of transversal Laplace transform the field amplitude in the parabolic approximation is calculated in the two-dimensional free space using initial values of the amplitude specified on an arbitrary shaped monotonic curve. The obtained amplitude depends on one \it a priori unknown function, which can be found from a Volterra first kind integral equation. In a special case of field amplitude specified on a concave parabolic curve the exact solution is derived. Both solutions can be used to study the light propagation from surfaces of tilted objects of an arbitrary shape including grazing incidence X-ray mirrors. They can also find applications in the analysis of coherent imaging problems of the X-ray optics, in phase retrieval algorithms as well as to solve inverse problems in cases when the initial field amplitude is sought on a curved surface.
  • PDF
    Partial differential equations (p.d.e) equipped of spatial derivatives of fractional order capture anomalous transport behaviors observed in diverse fields of Science. A number of numerical methods approximate their solutions in dimension one. Focusing our effort on such p.d.e. in higher dimension with Dirichlet boundary conditions, we present an approximation based on Lattice Boltzmann Method with Bhatnagar-Gross-Krook (BGK) or Multiple-Relaxation-Time (MRT) collision operators. First, an equilibrium distribution function is defined for simulating space-fractional diffusion equations in dimensions 2 and 3. Then, we check the accuracy of the solutions by comparing with i) random walks derived from stable Lévy motion, and ii) exact solutions. Because of its additional freedom degrees, the MRT collision operator provides accurate approximations to space-fractional advection-diffusion equations, even in the cases which the BGK fails to represent because of anisotropic diffusion tensor or of flow rate destabilizing the BGK LBM scheme.
  • PDF
    In this work, we present a highly accurate spectral neighbor analysis potential (SNAP) model for molybdenum (Mo) developed through the rigorous application of machine learning techniques on large materials data sets. Despite Mo's importance as a structural metal, existing force fields for Mo based on the embedded atom and modified embedded atom methods still do not provide satisfactory accuracy on many properties. We will show that by fitting to the energies, forces and stress tensors of a large density functional theory (DFT)-computed dataset on a diverse set of Mo structures, a Mo SNAP model can be developed that achieves close to DFT accuracy in the prediction of a broad range of properties, including energies, forces, stresses, elastic constants, melting point, phonon spectra, surface energies, grain boundary energies, etc. We will outline a systematic model development process, which includes a rigorous approach to structural selection based on principal component analysis, as well as a differential evolution algorithm for optimizing the hyperparameters in the model fitting so that both the model error and the property prediction error can be simultaneously lowered. We expect that this newly developed Mo SNAP model will find broad applications in large-scale, long-time scale simulations.

Recent comments

Chris Ferrie May 01 2017 12:52 UTC

I've posted a public referee report on this paper here:

Salvatore Mandrà Nov 21 2014 13:50 UTC

The manuscript has been widely revised to focus the reader's attention on the proposed method and its application in presence of local disorder.


Salvatore, Gian Giacomo and Alán

Salvatore Mandrà Aug 01 2014 19:11 UTC

Thanks Dr. Hastings for your comment. It is true that the transverse field Ising model does not satisfy the requirements to apply our method with an exponential reduction. Indeed, the opposite would be quite impressive since the random Ising model is a NP-Hard problem and we ourselves would be prett

Matt Hastings Aug 01 2014 16:43 UTC

The "quite general" conditions seem not to include the transverse field Ising model, the subject of most of the intensive numerical work previously. Incidentally, the terms "Lanczos" and "Krylov subspace" might be helpful.

Steve Flammia Jan 10 2014 08:01 UTC

This paper is an important first step towards making numerical and computational studies more easily checkable and reproducible. Very interesting stuff.

Noon van der Silk Jan 10 2014 16:26 UTC

Interestingly, another thematically-similar paper popped up today - (posting the link here as I assume not many people are subbed to cs.DL)