Numerical Simulation of Beam Network Models

Network models are used as efficient representation of materials with complex, interconnected locally one-dimensional structures. They typically accurately capture the mechanical properties of a material, while substantially reducing computational cost by avoiding full three-dimensional resolution. Applications include the simulation of fiber-based materials, porous media, and biological systems such as vascular networks. This article focuses on two representative problems: a stationary formulation describing the elastic deformation of beam networks, and a time-dependent formulation modeling elastic wave propagation in such materials. We propose a two-level additive domain decomposition method to efficiently solve the linear system associated with the stationary problem, as well as the linear systems that arise at each time step of the time-dependent problem through implicit time discretization. We present a rigorous convergence analysis of the domain decomposition method when used as a preconditioner, quantifying the convergence rate with respect to network connectivity and heterogeneity. The efficiency and robustness of the proposed approach are demonstrated through numerical simulations of the mechanical properties of commercial-grade paperboard.

💡 Research Summary

This research presents an advanced numerical framework for the efficient simulation of complex, interconnected one-dimensional structures, such as fiber-based materials, porous media, and biological vascular networks. A significant challenge in simulating these systems is the prohibitive computational cost associated with full three-dimensional (3D) resolutions. To address this, the authors utilize beam network models, which provide a highly efficient representation of materials by capturing essential mechanical properties while avoiding the massive overhead of 3D discretization.

The study focuses on two fundamental physical problems: a stationary formulation for describing the elastic deformation of beam networks and a time-dependent formulation for modeling elastic wave propagation. For the time-dependent problem, the authors employ implicit time discretization, which necessitates solving large-scale linear systems at each time step. To tackle the computational complexity of these systems, the paper proposes a two-level additive domain decomposition method (DDM) used as an effective preconditioner.

A major scientific contribution of this work is the rigorous convergence analysis of the proposed domain decomposition method. The authors do not merely propose an algorithm but mathematically quantify the convergence rate in relation to the network’s connectivity and the degree of heterogeneity within the material. This provides critical insights into how the topological complexity and material variations affect the numerical stability and speed of the solver, ensuring the method’s reliability for large-scale applications.

The practical utility and robustness of the proposed approach are validated through numerical simulations of the mechanical properties of commercial-grade paperboard. As a real-world example of a complex fiber network, the paperboard simulation demonstrates that the proposed method can accurately and efficiently capture the mechanical behavior of industrial-grade materials. This research offers a powerful computational tool for material scientists and engineers, promising significant reductions in simulation time for studying complex structural networks in fields ranging from tissue engineering to advanced composite material design.