Implicit Tau-Leaping Algorithm for Stiff Chemical System Simulation

Resolve Bottlenecks,
Find Innovative Solutions
Generate Solutions

Solution Overview

Problem

Current methods for simulating chemical and biological systems are time-consuming and costly, particularly in drug development, due to the need for extensive laboratory experiments and clinical trials, and struggle with accurately modeling stiff systems, which can lead to unforeseen side effects and inefficiencies.

Innovation Solution

A simulation tool that uses an implicit tau-leaping algorithm with a bracketing algorithm to maximize the time interval for simulating stiff chemical or biological systems, allowing for user-specified accuracy and reducing the need for iterative solving, thereby accelerating simulations while maintaining stability.

Engineering Contradictions & Design Principles

VSEngineering Contradiction Analysis

1Productivity

If the time interval size is increased to accelerate simulation, then productivity improves, but measurement precision deteriorates

Engineering Contradiction:
Improvesimulation speedVSAvoidsimulation accuracy
Core Design Contradiction:
ProductivityVSMeasurement precision

Solution Approach 1:

The patent implements dynamic time interval adjustment where the simulation time step automatically adapts based on system state. The algorithm monitors reaction rates and species concentrations, increasing the time interval when changes are slow and decreasing it when rapid changes occur, thereby optimizing both simulation speed and accuracy throughout the simulation process

Inventive Principle:
Principle #15Dynamics

Solution Approach 2:

The patent employs feedback mechanisms where simulation results from each time interval are used to determine the appropriate time interval for the next step. Error estimates and convergence criteria guide the adjustment of time step size, ensuring that accuracy requirements are met while maximizing simulation efficiency

Inventive Principle:
Principle #23Feedback

2Measurement precision

If iterative solving methods are used to solve nonlinear algebraic equations, then measurement precision improves, but loss of time increases

Engineering Contradiction:
Improvesolution accuracyVSAvoidcomputational time
Core Design Contradiction:
Measurement precisionVSLoss of time

Solution Approach 1:

The patent applies preliminary action by implementing effective truncation strategies that eliminate unnecessary iterative steps. By analyzing the system characteristics beforehand and using intelligent initial guesses based on previous time steps, the algorithm reduces the number of iterations required to achieve convergence, saving computational time while maintaining solution accuracy

Inventive Principle:
Principle #10Preliminary action

Solution Approach 2:

The patent utilizes parameter changes by dynamically adjusting solver tolerance levels and maximum iteration counts based on the current simulation state. For well-behaved systems or later time points where high precision is less critical, looser tolerances are applied to reduce computational burden, while stricter tolerances are used when high accuracy is required

Inventive Principle:
Principle #35Parameter changes

Data Source

PatentUS9811640B1Automatic selection of time interval size in implicit tau-leaping algorithm
Publication Date: 2017.11.07 MATHWORKS INC
  • US9811640B1 patent drawing
  • US9811640B1 patent drawing
  • US9811640B1 patent drawing

AI summary

A modeling and simulation environment enables a user to create a model of a stiff chemical or biological system. The stiff chemical or biological system refers to a system that includes at least two or more different time scales of the chemical or biochemical reactions. The modeling and simulation environment also enables the user to solve the stiff chemical or biological system using a leaping algorithm. The leaping algorithm may leap over simulation of the reactions occurring in the time interval to accelerate the simulation of the chemical or biological system. The algorithm maximizes the size of the time interval by adjusting the size of the time interval using a bracketing algorithm.