Related Experiment Video
Updated: Jun 26, 2026

Mechanical Control of Relaxation Using Intact Cardiac Trabeculae
Published on: February 17, 2023
Integrating biosystem models using waveform relaxation.
Linzhong Li1, Robert M Seymour, Stephen Baigent
1Institute for Energy Technology, Kjeller, Norway.
The waveform relaxation (WR) method efficiently computes complex biological models by integrating submodels. This parallel computation approach handles diverse model couplings and time scales for systems biology.
Area of Science:
- Systems Biology
- Computational Biology
- Mathematical Modelling
Background:
- Integrating component models into larger composite models is a key challenge in systems biology.
- Coupling of submodels can be unidirectional or bidirectional with variable strengths, complicating integration.
- Efficient and systematic methods are needed for computing complex, linked biological systems.
Purpose of the Study:
- To adapt the waveform relaxation (WR) method for parallel computation of ordinary differential equations (ODEs).
- To present WR as a general methodology for computing systems of linked submodels in systems biology.
- To demonstrate the flexibility of WR in handling multitime-scale computation and model heterogeneity.
Main Methods:
- Adaptation of the waveform relaxation (WR) method for parallel computation of ODEs.
- Application of WR to four distinct test cases involving coupled harmonic oscillators and calcium oscillations.
- Testing WR on single-cell and multicellular models, including complex behaviors like bursting and chaotic dynamics.
Main Results:
- The WR method successfully computed systems of linked submodels across diverse scenarios.
- Demonstrated WR's capability in handling both deterministic and stochastic simulations of calcium oscillations.
- Validated WR's effectiveness for multitime-scale computations and heterogeneous model components.
Conclusions:
- The waveform relaxation method offers a flexible approach for computing complex biological systems.
- WR enables the capture of global solutions independent of individual component solution techniques.
- This methodology facilitates parallel computation and integration of diverse submodels in systems biology.
Related Concept Videos
BIBO stability of continuous and discrete -time systems
To determine the BIBO stability, the convolution integral is utilized when a bounded continuous-time input is applied to a Linear Time-Invariant (LTI) system.
Parameters Affecting Nonlinear Elimination: Zero-Order Input, First-Order Absorption and Two-Compartment Model
When a drug is administered through a constant intravenous infusion and eliminated via nonlinear pharmacokinetics, it follows zero-order input. For example, oral drugs undergo first-order absorption upon administration and are eliminated through nonlinear pharmacokinetics.
In the case of subcutaneously administered drugs,...
Model Approaches for Pharmacokinetic Data: Physiological Models
Linear time-invariant Systems
The input-output behavior of an LTI system can be fully defined by its response to an impulsive excitation at its input. Once this impulse response is known, the system's reaction to any other input can be calculated...
Multi-input and Multi-variable systems
In the absence of...
Classification of Systems-II