Computing Residual Diffusivity by Adaptive Basis Learning via Spectral Method

Computing Residual Diffusivity by Adaptive Basis Learning via Spectral Method

Year:    2017

Numerical Mathematics: Theory, Methods and Applications, Vol. 10 (2017), Iss. 2 : pp. 351–372

Abstract

We study the residual diffusion phenomenon in chaotic advection computationally via adaptive orthogonal basis. The chaotic advection is generated by a class of time periodic cellular flows arising in modeling transition to turbulence in Rayleigh-Bénard experiments. The residual diffusion refers to the non-zero effective (homogenized) diffusion in the limit of zero molecular diffusion as a result of chaotic mixing of the streamlines. In this limit, the solutions of the advection-diffusion equation develop sharp gradients, and demand a large number of Fourier modes to resolve, rendering computation expensive. We construct adaptive orthogonal basis (training) with built-in sharp gradient structures from fully resolved spectral solutions at few sampled molecular diffusivities. This is done by taking snapshots of solutions in time, and performing singular value decomposition of the matrix consisting of these snapshots as column vectors. The singular values decay rapidly and allow us to extract a small percentage of left singular vectors corresponding to the top singular values as adaptive basis vectors. The trained orthogonal adaptive basis makes possible low cost computation of the effective diffusivities at smaller molecular diffusivities (testing). The testing errors decrease as the training occurs at smaller molecular diffusivities. We make use of the Poincaré map of the advection-diffusion equation to bypass long time simulation and gain accuracy in computing effective diffusivity and learning adaptive basis. We observe a non-monotone relationship between residual diffusivity and the amount of chaos in the advection, though the overall trend is that sufficient chaos leads to higher residual diffusivity.

You do not have full access to this article.

Already a Subscriber? Sign in as an individual or via your institution

Journal Article Details

Publisher Name:    Global Science Press

Language:    English

DOI:    https://doi.org/10.4208/nmtma.2017.s08

Numerical Mathematics: Theory, Methods and Applications, Vol. 10 (2017), Iss. 2 : pp. 351–372

Published online:    2017-01

AMS Subject Headings:   

Copyright:    COPYRIGHT: © Global Science Press

Pages:    22

Keywords:   

  1. Computing Effective Diffusivity of Chaotic and Stochastic Flows Using Structure-Preserving Schemes

    Wang, Zhongjian | Xin, Jack | Zhang, Zhiwen

    SIAM Journal on Numerical Analysis, Vol. 56 (2018), Iss. 4 P.2322

    https://doi.org/10.1137/18M1165219 [Citations: 10]
  2. An augmented subspace based adaptive proper orthogonal decomposition method for time dependent partial differential equations

    Dai, Xiaoying | Hu, Miao | Xin, Jack | Zhou, Aihui

    Journal of Computational Physics, Vol. 514 (2024), Iss. P.113231

    https://doi.org/10.1016/j.jcp.2024.113231 [Citations: 0]
  3. DeepParticle: Learning invariant measure by a deep neural network minimizing Wasserstein distance on data generated from an interacting particle method

    Wang, Zhongjian | Xin, Jack | Zhang, Zhiwen

    Journal of Computational Physics, Vol. 464 (2022), Iss. P.111309

    https://doi.org/10.1016/j.jcp.2022.111309 [Citations: 4]
  4. Two-Grid Based Adaptive Proper Orthogonal Decomposition Method for Time Dependent Partial Differential Equations

    Dai, Xiaoying | Kuang, Xiong | Xin, Jack | Zhou, Aihui

    Journal of Scientific Computing, Vol. 84 (2020), Iss. 3

    https://doi.org/10.1007/s10915-020-01288-9 [Citations: 2]