scholarly journals Model Reduction for Kinetic Models of Biological Systems

Symmetry ◽  
2020 ◽  
Vol 12 (5) ◽  
pp. 863
Author(s):  
Neveen Ali Eshtewy ◽  
Lena Scholz

High dimensionality continues to be a challenge in computational systems biology. The kinetic models of many phenomena of interest are high-dimensional and complex, resulting in large computational effort in the simulation. Model order reduction (MOR) is a mathematical technique that is used to reduce the computational complexity of high-dimensional systems by approximation with lower dimensional systems, while retaining the important information and properties of the full order system. Proper orthogonal decomposition (POD) is a method based on Galerkin projection that can be used for reducing the model order. POD is considered an optimal linear approach since it obtains the minimum squared distance between the original model and its reduced representation. However, POD may represent a restriction for nonlinear systems. By applying the POD method for nonlinear systems, the complexity to solve the nonlinear term still remains that of the full order model. To overcome the complexity for nonlinear terms in the dynamical system, an approach called the discrete empirical interpolation method (DEIM) can be used. In this paper, we discuss model reduction by POD and DEIM to reduce the order of kinetic models of biological systems and illustrate the approaches on some examples. Additional computational costs for setting up the reduced order system pay off for large-scale systems. In general, a reduced model should not be expected to yield good approximations if different initial conditions are used from that used to produce the reduced order model. We used the POD method of a kinetic model with different initial conditions to compute the reduced model. This reduced order model is able to predict the full order model for a variety of different initial conditions.

2004 ◽  
Vol 126 (4) ◽  
pp. 496-513 ◽  
Author(s):  
Deman Tang ◽  
Earl H. Dowell

Dynamic numerical simulation of a protein-ligand molecular chain connected to a moving atomic force microscope (AFM) has been studied. A sinusoidal base excitation of the cantilevered beam of the AFM is considered in some detail. A comparison between results for a single molecule and those for multiple molecules has been made. For a small number of molecules, multiple stable static equilibrium positions are observed and chaotic behavior may be generated via a period-doubling cascade for harmonic base excitation of the AFM. For many molecules in the chain, only a single static equilibrium position exists. To enable these calculations, reduced-order (dynamic) models are constructed for fully linear, combined linear/nonlinear and fully nonlinear systems. Several distinct reduced-order models have been developed that offer the option of increased computational efficiency at the price of greater effort to construct the particular reduced-order model. The agreement between the original and reduced-order models (ROM) is very good even when only one mode is included in the ROM for either the fully linear or combined linear/nonlinear systems provided the excitation frequency is lower than the fundamental natural frequency of the linear system. The computational advantage of the reduced-order model is clear from the results presented.


2012 ◽  
Vol 6-7 ◽  
pp. 135-142
Author(s):  
Xue Song Han ◽  
Yu Bo Duan

This paper extends the results obtained for one-dimensional Markovian jump systems to investigate the problem of H∞model reduction for a class of linear discrete time 2D Markovian jump systems with state delays in Roesser model which is time-varying and mode-independent. The reduced-order model with the same randomly jumping parameters is proposed which can make the error systems stochastically stable with a prescribed H∞ performance. A sufficient condition in terms of linear matrix inequalitiesSubscript text(LMIs) plus matrix inverse constraints are derived for the existence of a solution to the reduced-order model problems. The cone complimentarity linearization (CCL) method is exploited to cast them into nonlinear minimization problems subject to LMI constraints. A numerical example is given to illustrate the design procedures.


Author(s):  
Saurabh Biswas ◽  
Anindya Chatterjee

Hysteresis in material behaviour includes both signum nonlinearities as well as high dimensionality. Available models for component-level hysteretic behaviour are empirical. Here, we derive a low-order model for rate-independent hysteresis from a high-dimensional massless frictional system. The original system, being given in terms of signs of velocities, is first solved incrementally using a linear complementarity problem formulation. From this numerical solution, to develop a reduced-order model, basis vectors are chosen using the singular value decomposition. The slip direction in generalized coordinates is identified as the minimizer of a dissipation-related function. That function includes terms for frictional dissipation through signum nonlinearities at many friction sites. Luckily, it allows a convenient analytical approximation. Upon solution of the approximated minimization problem, the slip direction is found. A final evolution equation for a few states is then obtained that gives a good match with the full solution. The model obtained here may lead to new insights into hysteresis as well as better empirical modelling thereof.


Author(s):  
Sangram Redkar ◽  
S. C. Sinha

In this work, some techniques for order reduction of nonlinear systems with periodic coefficients subjected to external periodic excitations are presented. The periodicity of the linear terms is assumed to be non-commensurate with the periodicity of forcing vector. The dynamical equations of motion are transformed using the Lyapunov-Floquet (L-F) transformation such that the linear parts of the resulting equations become time-invariant while the forcing and/or nonlinearity takes the form of quasiperiodic functions. The techniques proposed here; construct a reduced order equivalent system by expressing the non-dominant states as time-varying functions of the dominant (master) states. This reduced order model preserves stability properties and is easier to analyze, simulate and control since it consists of relatively small number of states in comparison with the large scale system. Specifically, two methods are outlined to obtain the reduced order model. First approach is a straightforward application of linear method similar to the ‘Guyan reduction’, the second novel technique proposed here, utilizes the concept of ‘invariant manifolds’ for the forced problem to construct the fundamental solution. Order reduction approach based on invariant manifold technique yields unique ‘reducibility conditions’. If these ‘reducibility conditions’ are satisfied only then an accurate order reduction via ‘invariant manifold’ is possible. This approach not only yields accurate reduced order models using the fundamental solution but also explains the consequences of various ‘primary’ and ‘secondary resonances’ present in the system. One can also recover ‘resonance conditions’ associated with the fundamental solution which could be obtained via perturbation techniques by assuming weak parametric excitation. This technique is capable of handing systems with strong parametric excitations subjected to periodic and quasi-periodic forcing. These methodologies are applied to a typical problem and results for large-scale and reduced order models are compared. It is anticipated that these techniques will provide a useful tool in the analysis and control system design of large-scale parametrically excited nonlinear systems subjected to external periodic excitations.


Author(s):  
M. Ramasubramanian ◽  
M. Thirumarimurugan ◽  
P. Ananthi

Design of controller and analyzing the response of higher order system in real time environment would be very complex and expensive. Therefore, an attempt has been made in this paper to obtain the reduced order model of single-shaft Heavy duty gas turbine plants ranging from 18.2 to 106.7 MW by using various model order reduction techniques. The step response of Heavy duty gas turbine model using the reduced order models are compared with that of the original MATLAB/ Simulink model. Various time domain specifications and performance index criteria have been considered for analyzing the responses. The simulation results show that the response obtained by Routh approximation-Pade approximation technique based reduced order model mimics the original, higher order Heavy Duty gas turbine response. It is also proposed in this paper to improve the response by optimizing the co-efficients of reduced order model using Particle Swarm Optimization technique. On comparing the simulation results, Particle Swarm Optimization technique based reduced order model yield better transient and steady state response as close to original higher order system and hence it is identified as an optimal reduced order model for all Heavy Duty gas turbine plants in grid connected operation


2014 ◽  
Vol 2014 ◽  
pp. 1-5 ◽  
Author(s):  
Avadh Pati ◽  
Awadhesh Kumar ◽  
Dinesh Chandra

A Padé approximation based technique for designing a suboptimal controller is presented. The technique uses matching of both time moments and Markov parameters for model order reduction. In this method, the suboptimal controller is first derived for reduced order model and then implemented for higher order plant by partial feedback of measurable states.


Author(s):  
Yuqing Xie ◽  
Lin Li ◽  
Shuaibing Wang

Purpose To reduce the computational scale for quasi-magnetostatic problems, model order reduction is a good option. Reduced-order modelling techniques based on proper orthogonal decomposition (POD) and centroidal Voronoi tessellation (CVT) have been used to solve many engineering problems. The purpose of this paper is to investigate the computational principle, accuracy and efficiency of the POD-based and the CVT-based reduced-order method when dealing with quasi-magnetostatic problems. Design/methodology/approach The paper investigates computational features of the reduced-order method based on POD and CVT methods for quasi-magnetostatic problems. Firstly the construction method for the POD and the CVT reduced-order basis is introduced. Then, a reduced model is constructed using high-fidelity finite element solutions and a Galerkin projection. Finally, the transient quasi-magnetostatic problem of the TEAM 21a model is studied with the proposed reduced-order method. Findings For the TEAM 21a model, the numerical results show that both POD-based and CVT-based reduced-order approaches can greatly reduce the computational time compared with the full-order finite element method. And the results obtained from both reduced-order models are in good agreement with the results obtained from the full-order model, while the computational accuracy of the POD-based reduced-order model is a little higher than the CVT-based reduced-order model. Originality/value The CVT method is introduced to construct the reduced-order model for a quasi-magnetostatic problem. The computational accuracy and efficiency of the presented approaches are compared.


Author(s):  
Coşku Kasnakoğlu ◽  
R. Chris Camphouse ◽  
Andrea Serrani

In this paper, we consider a boundary control problem governed by the two-dimensional Burgers’ equation for a configuration describing convective flow over an obstacle. Flows over obstacles are important as they arise in many practical applications. Burgers’ equations are also significant as they represent a simpler form of the more general Navier–Stokes momentum equation describing fluid flow. The aim of the work is to develop a reduced-order boundary control-oriented model for the system with subsequent nonlinear control law design. The control objective is to drive the full order system to a desired 2D profile. Reduced-order modeling involves the application of an L2 optimization based actuation mode expansion technique for input separation, demonstrating how one can obtain a reduced-order Galerkin model in which the control inputs appear as explicit terms. Controller design is based on averaging and center manifold techniques and is validated with full order numerical simulation. Closed-loop results are compared to a standard linear quadratic regulator design based on a linearization of the reduced-order model. The averaging∕center manifold based controller design provides smoother response with less control effort and smaller tracking error.


Author(s):  
H Mallesam Dora ◽  

In this paper the Modified Routh Approximation (MRA) and Affine Arithmetic (AA) methods are investigates for obtaining the reduced order model (ROM) of SISO, discrete & MIMO uncertain systems into lower order system. Rigorous study and analysis of physical system direct to the outcome of systems with uncertainty instead of certain coefficients. Thus, systems having uncertain but bounded parameters known as uncertain systems are under consideration in this paper. An effective algorithm to determine the reduced order model is proposed here. This proposed methodology is verified using numerical examples available from the literature.


Sign in / Sign up

Export Citation Format

Share Document