Samuel Heroy, D. Taylor, F. Shi, M. Forest, P. Mucha
{"title":"Rigidity Percolation in Disordered 3D Rod Systems","authors":"Samuel Heroy, D. Taylor, F. Shi, M. Forest, P. Mucha","doi":"10.1137/21m1401206","DOIUrl":"https://doi.org/10.1137/21m1401206","url":null,"abstract":"In composite materials composed of soft polymer matrix and stiff, high-aspect-ratio particles, the composite undergoes a transition in mechanical strength when the inclusion phase surpasses a critical density. This phenomenon (rheological or mechanical percolation) is well-known to occur in many composites at a critical density that exceeds the conductivity percolation threshold. Conductivity percolation occurs as a consequence of contact percolation, which refers to the conducting particles' formation of a connected component that spans the composite. Rheological percolation, however, has evaded a complete theoretical explanation and predictive description. A natural hypothesis is that rheological percolation arises due to rigidity percolation, whereby a rigid component of inclusions spans the composite. We model composites as random isotropic dispersions of soft-core rods, and study rigidity percolation in such systems. Building on previous results for two-dimensional systems, we develop an approximate algorithm that identifies spanning rigid components through iteratively identifying and compressing provably rigid motifs -- equivalently, decomposing giant rigid components into rigid assemblies of successively smaller rigid components. We apply this algorithm to random rod systems to estimate a rigidity percolation threshold and explore its dependence on rod aspect ratio. We show that this transition point, like the contact percolation transition point, scales inversely with the average (aspect ratio-dependent) rod excluded volume. However, the scaling of the rigidity percolation threshold, unlike the contact percolation scaling, is valid for relatively low aspect ratio. Moreover, the critical rod contact number is constant for aspect ratio above some relatively low value; and lies below the prediction from Maxwell's isostatic condition.","PeriodicalId":313703,"journal":{"name":"Multiscale Model. Simul.","volume":"21 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2021-03-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"134012096","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Online adaptive algorithm for Constraint Energy Minimizing Generalized Multiscale Discontinuous Galerkin Method","authors":"Sai-Mang Pun, Siu Wun Cheung","doi":"10.1137/21m1402625","DOIUrl":"https://doi.org/10.1137/21m1402625","url":null,"abstract":"In this research, we propose an online basis enrichment strategy within the framework of a recently developed constraint energy minimizing generalized multiscale discontinuous Galerkin method (CEM-GMsDGM). Combining the technique of oversampling, one makes use of the information of the current residuals to adaptively construct basis functions in the online stage to reduce the error of multiscale approximation. A complete analysis of the method is presented, which shows the proposed online enrichment leads to a fast convergence from multiscale approximation to the fine-scale solution. The error reduction can be made sufficiently large by suitably selecting oversampling regions and the number of oversampling layers. Further, the convergence rate of the enrichment algorithm depends on a factor of exponential decay regarding the number of oversampling layers and a user-defined parameter. Numerical results are provided to demonstrate the effectiveness and efficiency of the proposed online adaptive algorithm.","PeriodicalId":313703,"journal":{"name":"Multiscale Model. Simul.","volume":"35 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2021-03-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"121307542","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Statistical Learning of Nonlinear Stochastic Differential Equations from Nonstationary Time Series using Variational Clustering","authors":"V. Boyko, S. Krumscheid, N. Vercauteren","doi":"10.1137/21m1403989","DOIUrl":"https://doi.org/10.1137/21m1403989","url":null,"abstract":"Parameter estimation for non-stationary stochastic differential equations (SDE) with an arbitrary nonlinear drift, and nonlinear diffusion is accomplished in combination with a non-parametric clustering methodology. Such a model-based clustering approach includes a quadratic programming (QP) problem with equality and inequality constraints. We couple the QP problem to a closed-form likelihood function approach based on suitable Hermite expansion to approximate the parameter values of the SDE model. The classification problem provides a smooth indicator function, which enables us to recover the underlying temporal parameter modulation of the one-dimensional SDE. The numerical examples show that the clustering approach recovers a hidden functional relationship between the SDE model parameters and an additional auxiliary process. The study builds upon this functional relationship to develop closed-form, non-stationary, data-driven stochastic models for multiscale dynamical systems in real-world applications.","PeriodicalId":313703,"journal":{"name":"Multiscale Model. Simul.","volume":"138 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2021-02-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"122063278","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Higher order phase averaging for highly oscillatory systems","authors":"W. Bauer, C. Cotter, B. Wingate","doi":"10.1137/21m1430546","DOIUrl":"https://doi.org/10.1137/21m1430546","url":null,"abstract":"We introduce a higher order phase averaging method for nonlinear oscillatory systems. Phase averaging is a technique to filter fast motions from the dynamics whilst still accounting for their effect on the slow dynamics. Phase averaging is useful for deriving reduced models that can be solved numerically with more efficiency, since larger timesteps can be taken. Recently, Haut and Wingate (2014) introduced the idea of computing finite window numerical phase averages in parallel as the basis for a coarse propagator for a parallel-in-time algorithm. In this contribution, we provide a framework for higher order phase averages that aims to better approximate the unaveraged system whilst still filtering fast motions. Whilst the basic phase average assumes that the solution independent of changes of phase, the higher order method expands the phase dependency in a basis which the equations are projected onto. In this new framework, the original numerical phase averaging formulation arises as the lowest order version of this expansion. Our new projection onto functions that are $k$th degree polynomials in the phase gives rise to higher order corrections to the phase averaging formulation. We illustrate the properties of this method on an ODE that describes the dynamics of a swinging spring due to Lynch (2002). Although idealized, this model shows an interesting analogy to geophysical flows as it exhibits a slow dynamics that arises through the resonance between fast oscillations. On this example, we show convergence to the non-averaged (exact) solution with increasing approximation order also for finite averaging windows. At zeroth order, our method coincides with a standard phase average, but at higher order it is more accurate in the sense that solutions of the phase averaged model track the solutions of the unaveraged equations more accurately.","PeriodicalId":313703,"journal":{"name":"Multiscale Model. Simul.","volume":"39 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2021-02-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"124865586","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Inverse Random Potential Scattering for Elastic Waves","authors":"Jianliang Li, Peijun Li, Xu Wang","doi":"10.1137/22m1497183","DOIUrl":"https://doi.org/10.1137/22m1497183","url":null,"abstract":"This paper is concerned with the inverse elastic scattering problem for a random potential in three dimensions. Interpreted as a distribution, the potential is assumed to be a microlocally isotropic Gaussian random field whose covariance operator is a classical pseudo-differential operator. Given the potential, the direct scattering problem is shown to be well-posed in the distribution sense by studying the equivalent Lippmann--Schwinger integral equation. For the inverse scattering problem, we demonstrate that the microlocal strength of the random potential can be uniquely determined with probability one by a single realization of the high frequency limit of the averaged compressional or shear backscattered far-field pattern of the scattered wave. The analysis employs the integral operator theory, the Born approximation in the high frequency regime, the microlocal analysis for the Fourier integral operators, and the ergodicity of the wave field.","PeriodicalId":313703,"journal":{"name":"Multiscale Model. Simul.","volume":"206 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2021-02-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"122666321","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Potential Singularity Formation of Incompressible Axisymmetric Euler Equations with Degenerate Viscosity Coefficients","authors":"T. Hou, D. Huang","doi":"10.1137/22m1470906","DOIUrl":"https://doi.org/10.1137/22m1470906","url":null,"abstract":"In this paper, we present strong numerical evidences that the incompressible axisymmetric Euler equations with degenerate viscosity coefficients and smooth initial data of finite energy develop a potential finite-time locally self-similar singularity at the origin. An important feature of this potential singularity is that the solution develops a two-scale traveling wave that travels towards the origin. The two-scale feature is characterized by the scaling property that the center of the traveling wave is located at a ring of radius $O((T-t)^{1/2})$ surrounding the symmetry axis while the thickness of the ring collapses at a rate $O(T-t)$. The driving mechanism for this potential singularity is due to an antisymmetric vortex dipole that generates a strong shearing layer in both the radial and axial velocity fields. Without the viscous regularization, the $3$D Euler equations develop a sharp front and some shearing instability in the far field. On the other hand, the Navier-Stokes equations with a constant viscosity coefficient regularize the two-scale solution structure and do not develop a finite-time singularity for the same initial data.","PeriodicalId":313703,"journal":{"name":"Multiscale Model. Simul.","volume":"6 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2021-02-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"125314446","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"High Order Numerical Homogenization for Dissipative Ordinary Differential Equations","authors":"Zeyu Jin, Ruo Li","doi":"10.1137/21m1397003","DOIUrl":"https://doi.org/10.1137/21m1397003","url":null,"abstract":"We propose a high order numerical homogenization method for dissipative ordinary differential equations (ODEs) containing two time scales. Essentially, only first order homogenized model globally in time can be derived. To achieve a high order method, we have to adopt a numerical approach in the framework of the heterogeneous multiscale method (HMM). By a successively refined microscopic solver, the accuracy improvement up to arbitrary order is attained providing input data smooth enough. Based on the formulation of the high order microscopic solver we derived, an iterative formula to calculate the microscopic solver is then proposed. Using the iterative formula, we develop an implementation to the method in an efficient way for practical applications. Several numerical examples are presented to validate the new models and numerical methods.","PeriodicalId":313703,"journal":{"name":"Multiscale Model. Simul.","volume":"7 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2021-02-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"124948762","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Multiscale Invariants of Floquet Topological Insulators","authors":"G. Bal, Daniel Massatt","doi":"10.1137/21m1392826","DOIUrl":"https://doi.org/10.1137/21m1392826","url":null,"abstract":"This paper analyzes Floquet topological insulators resulting from the time-harmonic irradiation of electromagnetic waves on two dimensional materials such as graphene. We analyze the bulk and edge topologies of approximations to the evolution of the light-matter interaction. Topologically protected interface states are created by spatial modulations of the drive polarization across an interface. In the high-frequency modulation regime, we obtain a sequence of topologies that apply to different time scales. Bulk-difference invariants are computed in detail and a bulk-interface correspondence is shown to apply. We also analyze a high-frequency high-amplitude modulation resulting in a large-gap effective topology topologically that remains valid only for moderately long times.","PeriodicalId":313703,"journal":{"name":"Multiscale Model. Simul.","volume":"56 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2021-01-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"128254732","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Homogenization of the Stokes System in a Non-Periodically Perforated Domain","authors":"S. Wolf","doi":"10.1137/21m1390815","DOIUrl":"https://doi.org/10.1137/21m1390815","url":null,"abstract":"In our recent work [8], we have studied the homogenization of the Poisson equation in a class of non periodically perforated domains. In this paper, we examine the case of the Stokes system. We consider a porous medium in which the characteristic distance between two holes, denoted by ε, is proportional to the characteristic size of the holes. It is well known (see [1],[17] and [19]) that, when the holes are periodically distributed in space, the velocity converges to a limit given by the Darcy’s law when the size of the holes tends to zero. We generalize these results to the setting of [8]. The non-periodic domains are defined as a local perturbation of a periodic distribution of holes. We obtain classical results of the homogenization theory in perforated domains (existence of correctors and regularity estimates uniform in ε) and we prove H−convergence estimates for particular force fields.","PeriodicalId":313703,"journal":{"name":"Multiscale Model. Simul.","volume":"6 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2021-01-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"130265961","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Continuum Model and Numerical Method for Dislocation Structure and Energy of Grain Boundaries","authors":"Xiaoxue Qin, Yejun Gu, Luchan Zhang, Y. Xiang","doi":"10.1137/20m1366782","DOIUrl":"https://doi.org/10.1137/20m1366782","url":null,"abstract":"We present a continuum model to determine the dislocation structure and energy of low angle grain boundaries in three dimensions. The equilibrium dislocation structure is obtained by minimizing the grain boundary energy that is associated with the constituent dislocations subject to the constraint of Frank's formula. The orientation-dependent continuous distributions of dislocation lines on grain boundaries are described conveniently using the dislocation density potential functions, whose contour lines on the grain boundaries represent the dislocations. The energy of a grain boundary is the total energy of the constituent dislocations derived from discrete dislocation dynamics model, incorporating both the dislocation line energy and reactions of dislocations. The constrained energy minimization problem is solved by the augmented Lagrangian method and projection method. Comparisons with atomistic simulation results show that our continuum model is able to give excellent predictions of the energy and dislocation densities of both planar and curved low angle grain boundaries.","PeriodicalId":313703,"journal":{"name":"Multiscale Model. Simul.","volume":"22 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2021-01-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"130989349","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}