Mathematical modeling of gate control theory

Agi, Egemen
The purpose of this thesis work is to model the gate control theory, which explains the modulation of pain signals, with a motivation of finding new possible targets for pain treatment and to find novel control algorithms that can be used in engineering practice. The difference of the current study from the previous modeling trials is that morphologies of neurons that constitute gate control system are also included in the model by which structure-function relationship can be observed. Model of an excitable neuron is constructed and the response of the model for different perturbations are investigated. The simulation results of the excitable cell model is obtained and when compared with the experimental findings obtained by using crayfish, it is found that they are in good agreement. Model encodes stimulation intensity information as firing frequency and also it can add sub-threshold inputs and fire action potentials as real neurons. Moreover, model is able to predict depolarization block. Absolute refractory period of the single cell model is found as 3.7 ms. The developed model, produces no action potentials when the sodium channels are blocked by tetrodotoxin. Also, frequency and amplitudes of generated action potentials increase when the reversal potential of Na is increased. In addition, propagation of signals along myelinated and unmyelinated fibers is simulated and input current intensity-frequency relationships for both type of fibers are constructed. Myelinated fiber starts to conduct when current input is about 400 pA whereas this minimum threshold value for unmyelinated fiber is around 1100 pA. Propagation velocity in the 1 cm long unmyelinated fiber is found as 0.43 m/s whereas velocity along myelinated fiber with the same length is found to be 64.35 m/s. Developed synapse model exhibits the summation and tetanization properties of real synapses while simulating the time dependency of neurotransmitter concentration in the synaptic cleft. Morphometric analysis of neurons that constitute gate control system are done in order to find electrophysiological properties according to dimensions of the neurons. All of the individual parts of the gate control system are connected and the whole system is simulated. For different connection configurations, results of the simulations predict the observed phenomena for the suppression of pain. If the myelinated fiber is dissected, the projection neuron generates action potentials that would convey to brain and elicit pain. However, if the unmyelinated fiber is dissected, projection neuron remains silent. In this study all of the simulations are preformed using Simulink.


Numerical simulation of laminar reaction flows
Tarhan, Tanıl; Selçuk, Nevin; Department of Chemical Engineering (2004)
Novel sequential and parallel computational fluid dynamic (CFD) codes based on method of lines (MOL) approach were developed for the numerical simulation of multi-component reacting flows using detailed transport and thermodynamic models. Both codes were applied to the prediction of a confined axisymmetric laminar co-flowing methane-air diffusion flame for which experimental data were available in the literature. Flame-sheet model for infinite-rate chemistry and one-, two-, and five- and ten-step reduced fi...
YAPICI, KERİM; Uludağ, Yusuf (FapUNIFESP (SciELO), 2013-10-01)
In this work, computer simulation results of steady incompressible flow in a 2-D square lid-driven cavity up to Reynolds number (Re) 65000 are presented and compared with those of earlier studies. The governing flow equations are solved by using the finite volume approach. Quadratic upstream interpolation for convective kinematics (QUICK) is used for the approximation of the convective terms in the flow equations. In the implementation of QUICK, the deferred correction technique is adopted. A non-uniform st...
Model Predictive control(MPC)performance for controlling reaction systems
Aşar, Işık; Özgen, Canan; Department of Chemical Engineering (2004)
In this study, the performance of the Model Predictive Controller (MPC) algorithm is investigated in two different reaction systems. The first case is a saponification reaction system where ethyl acetate reacts with sodium hydroxide to produce sodium acetate and ethanol in a CSTR. In the reactor, temperature and sodium acetate concentration are controlled by manipulating the flow rates of ethyl acetate and cooling water. The model of the reactor is developed considering first principal models. The experimen...
Numerical simulation of radiating flows
Karaismail, Ertan; Selçuk, Nevin; Department of Chemical Engineering (2005)
Predictive accuracy of the previously developed coupled code for the solution of the time-dependent Navier-Stokes equations in conjunction with the radiative transfer equation was first assessed by applying it to the prediction of thermally radiating, hydrodynamically developed laminar pipe flow for which the numerical solution had been reported in the literature. The effect of radiation on flow and temperature fields was demonstrated for different values of conduction to radiation ratio. It was found that ...
An empirical method for the second viral coefficients of non-standard fluids
Kis, Konrad; Orbey, Hasan (Elsevier BV, 1989-9)
A new empirical method is proposed for the extension of Pitzer-Curl type correlations of the second virial coefficient to non-standard fluids as define
Citation Formats
E. Agi, “Mathematical modeling of gate control theory,” M.S. - Master of Science, Middle East Technical University, 2009.