Abstract

Proper orthogonal decomposition (POD) and its extension based on time-windows have been shown to greatly improve the effectiveness of recovering smooth ensemble solutions from noisy particle data. However, to successfully de-noise any molecular system, a large number of measurements still need to be provided. In order to achieve a better efficiency in processing time-dependent fields, we have combined POD with a well-established signal processing technique, wavelet-based thresholding. In this novel hybrid procedure, the wavelet filtering is applied within the POD domain and referred to as WAVinPOD. The algorithm exhibits promising results when applied to both synthetically generated signals and particle data. In this work, the simulations compare the performance of our new approach with standard POD or wavelet analysis in extracting smooth profiles from noisy velocity and density fields. Numerical examples include molecular dynamics and dissipative particle dynamics simulations of unsteady force- and shear-driven liquid flows, as well as phase separation phenomenon. Simulation results confirm that WAVinPOD preserves the dimensionality reduction obtained using POD, while improving its filtering properties through the sparse representation of data in wavelet basis. This paper shows that WAVinPOD outperforms the other estimators for both synthetically generated signals and particle-based measurements, achieving a higher signal-to-noise ratio from a smaller number of samples. The new filtering methodology offers significant computational savings, particularly for multi-scale applications seeking to couple continuum informations with atomistic models. It is the first time that a rigorous analysis has compared de-noising techniques for particle-based fluid simulations.

Highlights

  • Particle-based simulations, e.g. molecular dynamics, are well suited to investigate the effects of fluid–solid interactions and are widely used to study a broad range of complex physical phenomena [1]

  • We report the results of applying WPOD, 2D wavelet thresholding and WAVinPOD to velocity measurements from time-dependent channel flows and density profiles from a simulation of phase separation phenomenon, performed with either molecular dynamics (MD) or Dissipative particle dynamics (DPD)

  • Wavelet thresholding within Proper orthogonal decomposition (POD) extracted similar quality velocity profiles but WPOD required 10× more snapshots; WAVinPOD produced an ensemble with 18% higher signal-to-noise ratio (SNR) than WPOD, if we assume that the true solution is a set of profiles obtained with WPOD for NPOD = 5000

Read more

Summary

Introduction

Particle-based simulations, e.g. molecular dynamics, are well suited to investigate the effects of fluid–solid interactions and are widely used to study a broad range of complex physical phenomena [1]. / Journal of Computational Physics 321 (2016) 169–190 can be viewed as a coarse-graining of molecular dynamics (a DPD particle is a collection of MD molecules), allowing for mesoscale modelling of complex fluids, e.g. surfactant solutions. Uncertainty in the measurements is increased through thermal fluctuations introduced by thermostats and sampling with a finite number of particles. These effects are generally referred to as noise, and a major challenge in particle simulations is to filter, or de-noise, the fluctuations to obtain an accurate ensemble prediction. The development of an efficient filtering technique that provides clean particle distribution functions and smooth gradients, when coupling across different length and time-scales (multi-scale modelling), is highly desirable

Objectives
Methods
Findings
Conclusion
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call