Journal of Computational Dynamics
April 2022 , Volume 9 , Issue 2
Special issue on structural dynamical systems
Select all articles
We examine a variety of numerical methods that arise when considering dynamical systems in the context of physics-based simulations of deformable objects. Such problems arise in various applications, including animation, robotics, control and fabrication. The goals and merits of suitable numerical algorithms for these applications are different from those of typical numerical analysis research in dynamical systems. Here the mathematical model is not fixed a priori but must be adjusted as necessary to capture the desired behaviour, with an emphasis on effectively producing lively animations of objects with complex geometries. Results are often judged by how realistic they appear to observers (by the "eye-norm") as well as by the efficacy of the numerical procedures employed. And yet, we show that with an adjusted view numerical analysis and applied mathematics can contribute significantly to the development of appropriate methods and their analysis in a variety of areas including finite element methods, stiff and highly oscillatory ODEs, model reduction, and constrained optimization.
In this work we propose a new numerical scheme for solving Richards' equation within Gardner's framework and accomplishing mass conservation. In order to do so, we resort to Kirchhoff transformation of Richards' equation in mixed form, so to exploit specific Gardner model features, obtaining a linear second order partial differential equation. Then, leveraging the mass balance condition, we integrate both sides of the equation over a generic grid cell and discretize integrals using trapezoidal rule. This approach provides a linear non-homogeneous initial value problem with respect to the Kirchhoff transform variable, whose solution yields the sought numerical scheme. Such a scheme is proven to be
We explore the applicability of splitting methods involving complex coefficients to solve numerically the time-dependent Schrödinger equation. We prove that a particular class of integrators are conjugate to unitary methods for sufficiently small step sizes when applied to problems defined in the group
Periodic solutions of delay equations are usually approximated as continuous piecewise polynomials on meshes adapted to the solutions' profile. In practical computations this affects the regularity of the (coefficients of the) linearized system and, in turn, the effectiveness of assessing local stability by approximating the Floquet multipliers. To overcome this problem when computing multipliers by collocation, the discretization grid should include the piecewise adapted mesh of the computed periodic solution. By introducing a piecewise version of existing pseudospectral techniques, we explain why and show experimentally that this choice is essential in presence of either strong mesh adaptation or nontrivial multipliers whose eigenfunctions' profile is unrelated to that of the periodic solution.
This paper analyzes conservation issues in the discretization of certain stochastic dynamical systems by means of stochastic
Recent experimental observations reveal that local cellular contraction pulses emerge via a combination of fast positive and slow negative feedbacks based on a signal network composed of Rho, GEF and Myosin interactions [
We are interested in the numerical solution of coupled semilinear partial differential equations (PDEs) in two and three dimensions. Under certain assumptions on the domain, we take advantage of the Kronecker structure arising in standard space discretizations of the differential operators and illustrate how the resulting system of ordinary differential equations (ODEs) can be treated directly in matrix or tensor form. Moreover, in the framework of the proper orthogonal decomposition (POD) and the discrete empirical interpolation method (DEIM) we derive a two- and three-sided model order reduction strategy that is applied directly to the ODE system in matrix and tensor form respectively. We discuss how to integrate the reduced order model and, in particular, how to solve the tensor-valued linear system arising at each timestep of a semi-implicit time discretization scheme. We illustrate the efficiency of the proposed method through a comparison to existing techniques on classical benchmark problems such as the two- and three-dimensional Burgers equation.
In this study, the formation of the adult sea urchin shape is rationalized within the Turing's theory paradigm. The emergence of protrusions from the expanding underlying surface is described through a reaction-diffusion model with Gray-Scott kinetics on a growing oblate spheroid. The case of slow exponential isotropic growth is considered. The model is first studied in terms of the spatially homogenous equilibria and of the bifurcations involved. Turing diffusion-driven instability is shown to occur and the impact of the slow exponential growth on the resulting Turing regions adequately discussed. Numerical investigations validate the theoretical results showing that the combination between an inhibitor and an activator can result in a distribution of spot concentrations that underlies the development of ambulacral tentacles in the sea urchin's adult stage. Our findings pave the way for a model-driven experimentation that could improve the current biological understanding of the gene control networks involved in patterning.
The Virtual Element Method (VEM) is a Galerkin approximation method that extends the Finite Element Method (FEM) to polytopal meshes. In this paper, we present a conforming formulation that generalizes the Scott-Vogelius finite element method for the numerical approximation of the Stokes problem to polygonal meshes in the framework of the virtual element method. In particular, we consider a straightforward application of the virtual element approximation space for scalar elliptic problems to the vector case and approximate the pressure variable through discontinuous polynomials. We assess the effectiveness of the numerical approximation by investigating the convergence on a manufactured solution problem and a set of representative polygonal meshes. We numerically show that this formulation is convergent with optimal convergence rates except for the lowest-order case on triangular meshes, where the method coincides with the
We propose a numerical method for approximating integro-differential equations arising in age-of-infection epidemic models. The method is based on a non-standard finite differences approximation of the integral term appearing in the equation. The study of convergence properties and the analysis of the qualitative behavior of the numerical solution show that it preserves all the basic properties of the continuous model with no restrictive conditions on the step-length
In this paper, we introduce a rigorous computational approach to prove existence of rotation invariant patterns for a nonlinear Laplace-Beltrami equation posed on the 2-sphere. After changing to spherical coordinates, the problem becomes a singular second order boundary value problem (BVP) on the interval
The global structure of the atmosphere and the oceans is a continuous source of intriguing challenges in geophysical fluid dynamics (GFD). Among these, jets are determinant in the air and water circulation around the Earth. In the last fifty years, thanks to the development of more and more precise and extensive observations, it has been possible to study in detail the atmospheric formations of the giant-gas planets in the solar system. For those planets, jets are the dominant large scale structure. Starting from the 70s, various theories combining observations and mathematical models have been proposed in order to describe their formation and stability. In this paper, we propose a purely algebraic approach to describe the spontaneous formation of jets on a spherical domain. Analysing the algebraic properties of the 2D Euler equations, we give a characterization of the different jets' structures. The calculations are performed starting from the discrete Zeitlin model of the Euler equations. For this model, the classification of the jets' structures can be precisely described in terms of reductive Lie algebras decomposition. The discrete framework provides a simple tool for analysing both from a theoretical and and a numerical perspective the jets' formation. Furthermore, it allows to extend the results to the original Euler equations.
Classical symplectic partitioned Runge–Kutta methods can be obtained from a variational formulation where all the terms in the discrete Lagrangian are treated with the same quadrature formula. We construct a family of symplectic methods allowing the use of different quadrature formulas (primary and secondary) for different terms of the Lagrangian. In particular, we study a family of methods using Lobatto quadrature (with corresponding Lobatto IIIA-B symplectic pair) as a primary method and Gauss–Legendre quadrature as a secondary method. The methods have the same implicitness as the underlying Lobatto IIIA-B pair, and, in addition, they are P-stable, therefore suitable for application to highly oscillatory problems.
Add your name and e-mail address to receive news of forthcoming issues of this journal:
[Back to Top]