HOPF BIFURCATION IN A CALCIUM OSCILLATION MODEL AND ITS CONTROL: FREQUENCY DOMAIN APPROACH

2013 ◽  
Vol 23 (01) ◽  
pp. 1350012 ◽  
Author(s):  
YU CHANG ◽  
LILI ZHOU ◽  
JINLIANG WANG

The Hopf bifurcation in a calcium oscillation model is theoretically analyzed by Hopf bifurcation theory in frequency domain. Approximation expressions for frequencies and amplitudes of periodic orbits arising from Hopf bifurcation are provided by using second-order harmonic balance method. In addition, a new method is proposed to control the amplitudes of the periodic orbits. Numerical simulations show the effectiveness of the method for suppressing periodic oscillations.

2007 ◽  
Vol 17 (04) ◽  
pp. 1355-1366 ◽  
Author(s):  
WENWU YU ◽  
JINDE CAO

In this paper, a general two-neuron model with time delay is considered, where the time delay is regarded as a parameter. It is found that Hopf bifurcation occurs when this delay passes through a sequence of critical value. By analyzing the characteristic equation and using the frequency domain approach, the existence of Hopf bifurcation is determined. The stability of bifurcating periodic solutions are determined by the harmonic balance approach, Nyquist criterion and the graphic Hopf bifurcation theorem. Numerical results are given to justify the theoretical analysis.


2014 ◽  
Vol 2014 ◽  
pp. 1-8
Author(s):  
Hongliang Yao ◽  
Qian Zhao ◽  
Qi Xu ◽  
Bangchun Wen

The efficiency and accuracy of common time and frequency domain methods that are used to simulate the response of a rotor system with malfunctions are compared and analyzed. The Newmark method and the incremental harmonic balance method are selected as typical representatives of time and frequency domain methods, respectively. To improve the simulation efficiency, the fixed interface component mode synthesis approach is combined with the Newmark method and the receptance approach is combined with the incremental harmonic balance method. Numerical simulations are performed for rotor systems with single and double frequency excitations. The inherent characteristic that determines the efficiency of the two methods is analyzed. The results of the analysis indicated that frequency domain methods are suitable single and double frequency excitation rotor systems, whereas time domain methods are more suitable for multifrequency excitation rotor systems.


Author(s):  
Anna Engels-Putzka ◽  
Jan Backhaus ◽  
Christian Frey

This paper describes the development and initial application of an adjoint harmonic balance solver. The harmonic balance method is a numerical method formulated in the frequency domain which is particularly suitable for the simulation of periodic unsteady flow phenomena in turbomachinery. Successful applications of this method include unsteady aerodynamics as well as aeroacoustics and aeroelasticity. Here we focus on forced response due to the interaction of neighboring blade rows. In the CFD-based design and optimization of turbomachinery components it is often helpful to be able to compute not only the objective values — e.g. performance data of a component — themselves, but also their sensitivities with respect to variations of the geometry. An efficient way to compute such sensitivities for a large number of geometric changes is the application of the adjoint method. While this is frequently used in the context of steady CFD, it becomes prohibitively expensive for unsteady simulations in the time domain. For unsteady methods in the frequency domain, the use of adjoint solvers is feasible, but still challenging. The present approach employs the reverse mode of algorithmic differentiation (AD) to construct a discrete adjoint of an existing harmonic balance solver in the framework of an industrially applied CFD code. The paper discusses implemen-tational issues as well as the performance of the adjoint solver, in particular regarding memory requirements. The presented method is applied to compute the sensitivities of aeroelastic objectives with respect to geometric changes in a turbine stage.


Author(s):  
Loi¨c Salles ◽  
Laurent Blanc ◽  
Fabrice Thouverez ◽  
Alexander M. Gouskov ◽  
Pierrick Jean

Contact interfaces with dry friction are frequently used in turbomachinery. Dry friction damping produced by the sliding surfaces of these interfaces reduces the amplitude of bladed-disk vibration. The relative displacements at these interfaces lead to fretting-wear which reduces the average life expectancy of the structure. Frequency response functions are calculated numerically by using the multi-Harmonic Balance Method (mHBM). The Dynamic Lagrangian Frequency-Time method is used to calculate contact forces in the frequency domain. A new strategy for solving non-linear systems based on dual time stepping is applied. This method is faster than using Newton solvers. It was used successfully for solving Nonlinear CFD equations in the frequency domain. This new approach allows identifying the steady state of worn systems by integrating wear rate equations a on dual time scale. The dual time equations are integrated by an implicit scheme. Of the different orders tested, the first order scheme provided the best results.


Author(s):  
Laura Junge ◽  
Graham Ashcroft ◽  
Peter Jeschke ◽  
Christian Frey

Due to the relative motion between adjacent blade rows the aerodynamic flow fields within turbomachinery are normally dominated by deterministic, periodic phenomena. In the numerical simulation of such unsteady flows (nonlinear) frequency-domain methods are therefore attractive as they are capable of fully exploiting the given spatial and temporal periodicity, as well as capturing or modelling flow nonlinearity. Central to the efficiency and accuracy of such frequency-domain methods is the selection of the frequencies and the circumferential modes to be resolved in simulations. Whilst trivial in the context of the simulation of a single compressor- or turbine-stage, the choice of solution modes becomes substantially more involved in multi-stage configurations. In this work the importance of mode scattering, in the context of the unsteady aerodynamic field, is investigated and quantified. It is shown that scattered modes can substantially impact the unsteady flow field and are essential for the accurate modelling of wake propagation within multistage configurations. Furthermore, an iterative approach is outlined, based on the spectral analysis of the circumferential modes at the interfaces between blade rows, to identify the dominant solution modes that should be resolved in the adjacent blade row. To demonstrate the importance of mode scattering and validate the approach for their identification the unsteady blade row interaction within a 4.5 stage axial compressor is computed using both the harmonic balance method and, based on a full annulus midspan simulation, a time-domain method. Through the inclusion of scattered modes it is shown that the solution quality of the harmonic balance results is comparable to that of the nonlinear time-domain simulation.


2018 ◽  
Vol 140 (9) ◽  
Author(s):  
Elham Shamsara ◽  
Zahra Afsharnezhad ◽  
Elham Javidmanesh

In this paper, we present a discontinuous cytotoxic T cells (CTLs) response for HTLV-1. Moreover, a delay parameter for the activation of CTLs is considered. In fact, a system of differential equation with discontinuous right-hand side with delay is defined for HTLV-1. For analyzing the dynamical behavior of the system, graphical Hopf bifurcation is used. In general, Hopf bifurcation theory will help to obtain the periodic solutions of a system as parameter varies. Therefore, by applying the frequency domain approach and analyzing the associated characteristic equation, the existence of Hopf bifurcation by using delay immune response as a bifurcation parameter is determined. The stability of Hopf bifurcation periodic solutions is obtained by the Nyquist criterion and the graphical Hopf bifurcation theorem. At the end, numerical simulations demonstrated our results for the system of HTLV-1.


Author(s):  
Christian Frey ◽  
Graham Ashcroft ◽  
Hans-Peter Kersken

This paper compares various approaches to simulate unsteady blade row interactions in turbomachinery. Unsteady simulations of turbomachinery flows have gained importance over the last years since increasing computing power allows the user to consider 3D unsteady flows for industrially relevant configurations. Furthermore, for turbomachinery flows, the last two decades have seen considerable efforts in developing adequate CFD methods which exploit the rotational symmetries of blade rows and are therefore up to several orders of magnitude more efficient than the standard unsteady approach for full wheel configurations. This paper focusses on the harmonic balance method which has been developed recently by the authors. The system of equations as well as the iterative solver are formulated in the frequency domain. The aim of this paper is to compare the harmonic balance method with the time-linearized as well as the non-linear unsteady approach. For the latter the unsteady flow fields in a fan stage are compared to reference results obtained with a highly resolved unsteady simulation. Moreover the amplitudes of the acoustic modes which are due to the rotor stator interaction are compared to measurement data available for this fan stage. The harmonic balance results for different sets of harmonics in the blade rows are used to explain the minor discrepancies between the time-linearized and unsteady results published by the authors in previous publications. The results show that the differences are primarily due to the neglection of the two-way coupling in the time-linearized simulations.


Sign in / Sign up

Export Citation Format

Share Document