An official website of the United States government
Abstract Unconditionally stable time stepping schemes are useful and often practically necessary for advancing parabolic operators in multi-scale systems. However, serious accuracy problems may emerge when taking time steps that far exceed the explicit stability limits. In our previous work, we compared the accuracy and performance of advancing parabolic operators in a thermodynamic MHD model using an implicit method and an explicit super time-stepping (STS) method. We found that while the STS method outperformed the implicit one with overall good results, it was not able to damp oscillatory behavior in the solution efficiently, hindering its practical use. In this follow-up work, we evaluate an easy-to-implement method for selecting a practical time step limit (PTL) for unconditionally stable schemes. This time step is used to ‘cycle’ the operator-split thermal conduction and viscosity parabolic operators. We test the new time step with both an implicit and STS scheme for accuracy, performance, and scaling. We find that, for our test cases here, the PTL dramatically improves the STS solution, matching or improving the solution of the original implicit scheme, while retaining most of its performance and scaling advantages. The PTL shows promise to allow more accurate use of unconditionally stable schemes for parabolic operators and reliable use of STS methods.
more »
« less
Stability analysis of explicit MPM
Abstract In this paper we analyze the stability of the explicit material point method (MPM). We focus on PIC, APIC, and CPIC transfers using quadratic and cubic splines in two and three dimensions. We perform a fully three‐dimensional Von Neumann stability analysis to study the behavior within the bulk of a material. This reveals the relationship between the sound speed, CFL number, and actual time step restriction and its dependence on discretization options. We note that boundaries are generally less stable than the interior, with stable time steps generally decreasing until the limit when particles become isolated. We then analyze the stability of a single particle to derive a novel time step restriction that stabilizes simulations at their boundaries. Finally, we show that for explicit MPM with APIC or CPIC transfers, there are pathological cases where growth is observed at arbitrarily small time steps sizes. While these cases do not necessarily pose a problem for practical usage, they do suggest that a guarantee of stability may be theoretically impossible and that necessary but not sufficient time step restrictions may be a necessary and practical compromise.
more »
« less
- Award ID(s):
- 2006570
- PAR ID:
- 10402386
- Publisher / Repository:
- Wiley-Blackwell
- Date Published:
- Journal Name:
- Computer Graphics Forum
- Volume:
- 41
- Issue:
- 8
- ISSN:
- 0167-7055
- Page Range / eLocation ID:
- p. 19-30
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract BackgroundAnimal movement data are increasingly used to make ecological inferences, as well as to inform conservation and management actions. While advanced statistical methods to estimate behavioral states from these datasets have become widely available, the large number to choose from may make it difficult for practitioners to decide which method best addresses their needs. To guide decisions, we compared the behavioral state estimates and inferences from three methods (movement persistence models [MPM], hidden Markov models [HMM], and mixed-membership method for movement [M4]) when analyzing animal telemetry data. Tracks of post-breeding adult male green sea turtles (Chelonia mydas) were treated as an empirical example for this method comparison. The effect of temporal scale on behavioral state estimates was also investigated (at 1, 4, and 8 h time steps). ResultsThe HMM and M4 models produced relatively similar behavioral state estimates (compared to the MPM) and estimated anywhere from three to five states depending on the time interval of the tracks and the method used. Likewise, for all three methods, sampling movement at coarser time scales smoothed estimates of behavioral transitions. Additionally, the selection of movement metrics for analysis by the HMM and M4 also appeared to be a critical decision regarding state estimation and interpretation. At the longest time step (8 h), all three models were able to distinguish area-restricted search (ARS) behavior from migratory behavior, with greater nuance estimated by the HMM and M4 methods. By comparison, the MPM was the only model that was able to identify fine-scale behavioral patterns when analyzing the shortest time step (1 h). Moreover, the analysis of tracks with short time steps via MPM identified likely periods of resting during long-distance migration, which had only previously been hypothesized in green turtles. ConclusionsWhile there is no single best method to estimate behavioral states, our findings demonstrate that results can vary widely among different statistical methods and that model assumptions should be thoroughly checked during the model fitting process to reduce any potential biases. Thus, practitioners should carefully consider which methods best address their needs while also accounting for the inherent properties of their telemetry dataset.more » « less
-
A bstract We find a new on-shell replica wormhole in a computation of the generating functional of JT gravity coupled to matter. We show that this saddle has lower action than the disconnected one, and that it is stable under restriction to real Lorentzian sections, but can be unstable otherwise. The behavior of the classical generating functional thus may be strongly dependent on the signature of allowed perturbations. As part of our analysis, we give an LM-style construction for computing the on-shell action of replicated manifolds even as the number of boundaries approaches zero, including a type of one-step replica symmetry breaking that is necessary to capture the contribution of the new saddle. Our results are robust against quantum corrections; in fact, we find evidence that such corrections may sometimes stabilize this new saddle.more » « less
-
Achieving stable stress solutions at large strains using the Material Point Method (MPM) is challenging due to the accumulation of errors associated with geometry discretization, cell-crossing noise, and volumetric locking. Several simplified attempts exist in the literature to mitigate these errors, including higher-order frameworks. However, the stability of the MPM solution in such frameworks has been limited to simple geometries and the single-phase formulation (i.e., neglecting pore fluid). Although never explored, multipatch isogeometric analysis offers desirable qualities to simulate complex geometries while mitigating errors in the MPM. The degree of required high-order spatial integration has also never been investigated to infer a minimum limit for the stability of the stress solution in MPM. This paper presents a general-purpose numerical framework for simulating stable stresses in porous media, capturing both near incompressibility and multiphase interactions. First, the numerical framework is presented considering Non-Uniform Rational B-splines (NURBS) to perform isogeometric analysis (IGA) in MPM. Additionally, a volumetric strain smoothing algorithm is used to alleviate errors associated with volumetric locking. Second, the manifestation of cell-crossing errors is assessed via a series of problems with orders ranging from linear to cubic interpolation functions. Third, the use of NURBS is investigated and verified for problems with circular geometries. Finally, multipatch analysis is deployed to simulate plane strain and 3D penetration in soils, considering nearly incompressible elastoplastic (total stress) analysis and fully-coupled hydro-mechanical (effective stress) analysis. The stability of the solution is also analyzed for different constitutive models. From the results, it can be concluded that the framework using cubic interpolation functions with strain smoothing is the most convenient, presenting stable stress solutions for a broad range of multiphase geotechnical applications.more » « less
-
Abstract We study diagonally implicit Runge-Kutta (DIRK) schemes when applied to abstract evolution problems that fit into the Gelfand-triple framework. We introduce novel stability notions that are well-suited to this setting and provide simple, necessary and sufficient, conditions to verify that a DIRK scheme is stable in our sense and in Bochner-type norms. We use several popular DIRK schemes in order to illustrate cases that satisfy the required structural stability properties and cases that do not. In addition, under some mild structural conditions on the problem we can guarantee compactness of families of discrete solutions with respect to time discretization.more » « less
