Coupled Systems of Differential-Algebraic and Kinetic Equations with Application to the Mathematical Modelling of Muscle Tissue

We consider a coupled system composed of a linear differential-algebraic equation (DAE) and a linear large-scale system of ordinary differential equations where the latter stands for the dynamics of numerous identical particles. Replacing the discrete particles by a kinetic equation for a particle density, we obtain in the mean-field limit the new class of partially kinetic systems. We investigate the influence of constraints on the kinetic theory of those systems and present necessary adjustments. We adapt the mean-field limit to the DAE model and show that index reduction and the mean-field limit commute. As a main result, we prove Dobrushin's stability estimate for linear systems. The estimate implies convergence of the mean-field limit and provides a rigorous link between the particle dynamics and their kinetic description. Our research is inspired by mathematical models for muscle tissue where the macroscopic behaviour is governed by the equations of continuum mechanics, often discretised by the finite element method, and the microscopic muscle contraction process is described by Huxley's sliding filament theory. The latter represents a kinetic equation that characterises the state of the actin-myosin bindings in the muscle filaments. Linear partially kinetic systems are a simplified version of such models, with focus on the constraints.Comment: 32 pages, 18 figure

In silico assessment of biomedical products: the conundrum of rare but not so rare events in two case studies

In silico clinical trials, defined as “The use of individualized computer simulation in the development or regulatory evaluation of a medicinal product, medical device, or medical intervention,” have been proposed as a possible strategy to reduce the regulatory costs of innovation and the time to market for biomedical products. We review some of the the literature on this topic, focusing in particular on those applications where the current practice is recognized as inadequate, as for example, the detection of unexpected severe adverse events too rare to be detected in a clinical trial, but still likely enough to be of concern. We then describe with more details two case studies, two successful applications of in silico clinical trial approaches, one relative to the University of Virginia/Padova simulator that the Food and Drug Administration has accepted as possible replacement for animal testing in the preclinical assessment of artificial pancreas technologies, and the second, an investigation of the probability of cardiac lead fracture, where a Bayesian network was used to combine in vivo and in silico observations, suggesting a whole new strategy of in silico-augmented clinical trials, to be used to increase the numerosity where recruitment is impossible, or to explore patients’ phenotypes that are unlikely to appear in the trial cohort, but are still frequent enough to be of concern

Nonlinear Feed Forward Control of a Perturbed Satellite using ExtendedLeast Squares Adaptation and a Luenberger Observer

International audienc

Multiscale modeling of the neuromuscular system: Coupling neurophysiology and skeletal muscle mechanics

Mathematical models and computer simulations have the great potential to substantially increase our understanding of the biophysical behavior of the neuromuscular system. This, however, requires detailed multiscale, and multiphysics models. Once validated, such models allow systematic in silico investigations that are not necessarily feasible within experiments and, therefore, have the ability to provide valuable insights into the complex interrelations within the healthy system and for pathological conditions. Most of the existing models focus on individual parts of the neuromuscular system and do not consider the neuromuscular system as an integrated physiological system. Hence, the aim of this advanced review is to facilitate the prospective development of detailed biophysical models of the entire neuromuscular system. For this purpose, this review is subdivided into three parts. The first part introduces the key anatomical and physiological aspects of the healthy neuromuscular system necessary for modeling the neuromuscular system. The second part provides an overview on state-of-the-art modeling approaches representing all major components of the neuromuscular system on different time and length scales. Within the last part, a specific multiscale neuromuscular system model is introduced. The integrated system model combines existing models of the motor neuron pool, of the sensory system and of a multiscale model describing the mechanical behavior of skeletal muscles. Since many sub-models are based on strictly biophysical modeling approaches, it closely represents the underlying physiological system and thus could be employed as starting point for further improvements and future developments. This article is categorized under: Physiology > Mammalian Physiology in Health and Disease Analytical and Computational Methods > Computational Methods Models of Systems Properties and Processes > Organ, Tissue, and Physiological Models

A Comparison Between Active Strain and Active Stress in Transversely Isotropic Hyperelastic Materials

Active materials are media for which deformations can occur in absence of loads, given an external stimulus. Two approaches to the modeling of such materials are mainly used in literature, both based on the introduction of a new tensor: an additive stress Pact in the active stress case and a multiplicative strain Fa in the active strain one. Aim of this paper is the comparison between the two approaches on simple shears. Considering an incompressible and transversely isotropic material, we design constitutive relations for Pact and Fa so that they produce the same results for a uniaxial deformation along the symmetry axis. We then study the two approaches in the case of a simple shear deformation. In a hyperelastic setting, we show that the two approaches produce different stress components along a simple shear, unless some necessary conditions on the strain energy density are fulfilled. However, such conditions are very restrictive and rule out the usual elastic strain energy functionals. Active stress and active strain therefore produce different results in shear, even if they both fit uniaxial data. Our results show that experimental data on the stress-stretch response on uniaxial deformations are not enough to establish which activation approach can capture better the mechanics of active materials. We conclude that other types of deformations, beyond the uniaxial one, should be taken into consideration in the modeling of such materials

Distributed multi-scale muscle simulation in a hybrid MPI–CUDA computational environment

We present Mexie, an extensible and scalable software solution for distributed multi-scale muscle simulations in a hybrid MPI–CUDA environment. Since muscle contraction relies on the integration of physical and biochemical properties across multiple length and time scales, these models are highly processor and memory intensive. Existing parallelization efforts for accelerating multi-scale muscle simulations imply the usage of expensive large-scale computational resources, which produces overwhelming costs for the everyday practical application of such models. In order to improve the computational speed within a reasonable budget, we introduce the concept of distributed calculations of multi-scale muscle models in a mixed CPU–GPU environment. The concept is applied to a two-scale muscle model, in which a finite element macro model is coupled with the microscopic Huxley kinetics model. Finite element calculations of a continuum macroscopic model take place strictly on the CPU, while numerical solutions of the partial differential equations of Huxley’s cross-bridge kinetics are calculated on both CPUs and GPUs. We present a modular architecture of the solution, along with an internal organization and a specific load balancer that is aware of memory boundaries in such a heterogeneous environment. Solution was verified on both benchmark and real-world examples, showing high utilization of involved processing units, ensuring high scalability. Speed-up results show a boost of two orders of magnitude over any previously reported distributed multi-scale muscle models. This major improvement in computational feasibility of multi-scale muscle models paves the way for new discoveries in the field of muscle modeling and future clinical applications.Author's versio