Viscoelastic behavior of erythrocyte membrane.

A nonlinear viscoelastic relation is developed to describe the viscoelastic properties of erythrocyte membrane. This constitutive equation is used in the analysis of the time-dependent aspiration of an erythrocyte membrane into a micropipette. Equations governing this motion are reduced to a nonlinear integral equation of the Volterra type. A numerical procedure based on a finite difference scheme is used to solve the integral equation and to match the experimental data. The data, aspiration length vs. time, is used to determine the relaxation function at each time step. The inverse problem of obtaining the time dependence of the aspiration length from a given relaxation function is also solved. Analytical results obtained are applied to the experimental data of Chien et al. 1978. Biophys. J. 24:463-487. A relaxation function similar to that of a four-parameter solid with a shear-thinning viscous term is proposed.     

Identification of Hammerstein model for bioreactors with input multiplicities

A closed loop identification method of Hammerstein model for continuous bioreactor with input multiplicity is proposed. Hammerstein model consists of nonlinear steady-state gain followed by a unity gain linear system. The method consists of first getting local first order plus time delay (FOPTD) models around the two input multiplicity values of the substrate feed concentration. The model parameters of the FOPTD is identified by a least square optimization method. The initial guess for the model parameters are obtained from the settling time, the initial delay in the closed loop servo response and using a simple proportional controller formula. From the local process gain values obtained for the several step changes around the two operating conditions, the nonlinear gain portion of the Hammerstein is then obtained. The actual nonlinear gain and the identified nonlinear gain is compared.     

A Nonlinearity in the Inhibitory Interactions in the Lateral Eye of Limulus

Receptor units in the eye of the horseshoe crab are more sensitive to lateral inhibition at some levels of excitation than they are at others. As a result, the steady-state inhibition of the response of a given unit is not directly proportional to the response levels of neighboring units. This effect may be represented by the introduction of a nonlinearity in the Hartline-Ratliff system of equations. The nonlinear inhibitory effect appears to increase the operating range of the receptor units.     

Random Search Algorithm for Solving the Nonlinear Fredholm Integral Equations of the Second Kind

In this paper, a randomized numerical approach is used to obtain approximate solutions for a class of nonlinear Fredholm integral equations of the second kind. The proposed approach contains two steps: at first, we define a discretized form of the integral equation by quadrature formula methods and solution of this discretized form converges to the exact solution of the integral equation by considering some conditions on the kernel of the integral equation. And then we convert the problem to an optimal control problem by introducing an artificial control function. Following that, in the next step, solution of the discretized form is approximated by a kind of Monte Carlo (MC) random search algorithm. Finally, some examples are given to show the efficiency of the proposed approach.     

A unified model for the combined temporal and spatial Broca-Sulzer effect

A briefly pulsed light is brighter than a pulse of longer duration. For a brief flash of fixed duration and small area, enlarging a target results initially in an increase in brightness. Beyond some critical area, the target dims with further increases in size. These two phenomena are the temporal and the spatial Broca-Sulzer effect respectively. Each effect can be modeled using a generalization of the Hartline-Ratliff equation. The resulting analysis produces a mathematically unified treatment of both effects.     

The identification of nonlinear biological systems: Wiener and Hammerstein cascade models

Systems that can be represented by a cascade of a dynamic linear subsystem preceded (Hammerstein cascade model) or followed (Wiener cascade model) by a static nonlinearity are considered. Various identification schemes that have been proposed for the Hammerstein and Wiener systems are critically reviewed with reference to the special problems that arise in the identification of nonlinear biological systems. Examples of Wiener and Hammerstein systems are identified from limited duration input-output data using an iterative identification scheme.     

Time-dependent ligand current into a saturating cell performing chemoreception

We determine the ligand current into a single spherical cell whose receptors become permanently blocked after binding ligands. Initially the cell is placed in a medium which contains ligands at uniform concentration. The analytical solution for the ligand distribution is obtained in terms of an integral over the solution at the cell surface. For the solution at the cell surface a nonlinear integral equation is derived which is solved numerically. We determine the time-dependent ligand current into the cell and the average number of free receptors in the cell surface as a function of time.     

Modeling of Coulomb collisions in a kinetic description of the electron cyclotron resonance plasma heating

The electron distribution function is modeled numerically with allowance for Coulomb collisions and quasilinear effects under cyclotron resonance conditions by solving a two-dimensional kinetic equation containing the quasilinear diffusion operator and the Coulomb collision operator in the Landau form. Two simplified model collision integrals that make it possible to describe electron heating by microwave radiation are considered. The first model collision operator is obtained by introducing the parametric time dependence of the temperature of the background Maxwellian electrons into the linear collision integral. It is shown that the heating of the bulk electrons can be described in a noncontradictory way if the temperature dynamics of the background electrons is calculated from the equation of energy balance, which is governed by the amount of the microwave power absorbed by the resonant electrons with the distribution function modified due to quasilinear effects. This conclusion is confirmed in a more rigorous fashion by comparing the solutions obtained using the first model Coulomb collision integral with those obtained using the second model integral, namely, the nonlinear operator derived by averaging the distribution function of the scattering electrons over pitch angles. The time-dependent linear collision integral is used to obtain analytic solutions describing quasi-steady electron heating with allowance for the quasilinear degradation of microwave power absorption.     

A model for the spatial spread of an epidemic

Summary We set up a deterministic model for the spatial spread of an epidemic. Essentially, the model consists of a nonlinear integral equation which has an unique solution. We show that this solution has a temporally asymptotic limit which describes the final state of the epidemic and is the minimal solution of another nonlinear integral equation. We outline the asymptotic behaviour of this minimal solution at a great distance from the epidemic's origin and generalize D. G. Kendall's pandemic threshold theorem (1957).     

Transient behavior of population density with competition for resources

In this paper we derive a perturbation expansion of some nonlinear diffusion equations proposed by J. G. Skellam and M. Kimura. The solution is derived by converting the differential equations into integral equations by means of the Green's function for the diffusion equation. This research was supported in part by the Office of Naval Research under Contract Nonr-595(17).     

Nonlinear feedforward control of bioreactors with input multiplicities

A steady-state nonlinear feedforward controller (FFC) for measurable disturbances is designed for a continuous bioreactor, which is represented by Hammerstein type nonlinear model wherein the nonlinearity is a polynomial with input multiplicities. The manipulated variable is the feed substrate concentration (S_f) and the disturbance variable is the dilution rate (D). The productivity (Q=DP) is considered as the controlled variable. The desired value of Q=3.73 gives two values of feed substrate concentration. The nonlinearity in the gain is considered for relating output to the manipulated variable and separately for the relation between output to disturbance variable. The FFC is also designed for the overall linearized system. The performance of the FFC is evaluated on the nonlinear differential equation model. The FFC is also designed for the model based on a single nonlinear steady-state equation containing both D and S_f. This nonlinear FFC gives the best performance. The nonlinear FFC is also designed by using only linear gain for the disturbance and nonlinear gain for the manipulated variable. Similarly, nonlinear FFC is also designed by using linear gain for the manipulated variable and the nonlinear gain for the disturbance variable. The performances of these FFC schemes are compared.     

EXISTENCE OF MULTIPLE PERIODIC SOLUTIONS TO A NONLINEAR INTEGRAL EQUATION FROM THE THEORY OF EPIDEMICS

Thenonlinearintegralequation(1.1)isageneralizedmodelforthespreadofdiseaseswithseasonaldependence.Inthispaper,wehaveprovedtheexistenceofatleastthreenontrivialnonnegativeperiodicsolutionstothisequation.     

A new asymptotic analysis of the nth order reaction-diffusion problem: analytical and numerical studies

We present a new formulation of the steady state, isothermal, nonlinear reaction-diffusion problem involving nth order reaction kinetics for slab geometry. This results in tractable expressions for the effectiveness factor as a function of the Thiele modulus, the Thiele modulus as a function of the centerline concentration, and the concentration profiles in the slab. The expressions are valid asymptotically in the limit of large orders n. We compare these results with the exact numerical solutions obtained by transforming the nonlinear differential equation into an integral form, using Green's function methods, and solving by successive approximations. The formulation for a membrane is also given, and the nature of the asymmetrical solution discussed. The analysis is facilitated through the introduction of pseudo-reaction orders. A comparison of the asymptotic Thiele modulus obtained herein with a previously given expression shows the present theory to be an improvement.     

Iterative algorithms for processing experimental data

The need to solve linear and nonlinear integral equations arise, e.g., in recovering plasma parameters from the data of multichannel diagnostics. The paper presents an iterative method for solving integral equations with a singularity at the upper limit of integration. The method consists in constructing successive approximations and calculating the integral by quadrature formulas in each integration interval. An example of application of the iterative algorithm to numerically solve an integral equation similar to those arising in recovering the plasma density profile from reflectometry data is presented.     

Thresholds and travelling waves for the geographical spread of infection

Summary A nonlinear integral equation of mixed Volterra-Fredholm type describing the spatio-temporal development of an epidemic is derived and analysed. Particular attention is paid to the hair-trigger effect and to the travelling wave problem.     

Stationary states of the Hartline-Ratliff model

For the Hartline-Ratliff model the exact conditions for the uniqueness of the stationary state are determined. Also sufficient conditions for dynamic stability are derived.     

Output regulation of a class of unstructured models of continuous bioreactors: steady-state approaches

This article deals with the output regulation of continuous bioreactors in the face of constant disturbances and inverse dynamics. Nonlinear controllers developed on the basis of approximate equilibrium manifolds can almost attenuate measurable or unmeasurable disturbances on the output. This nonlinear feed-forward/feedback control framework without any tuning parameters can be directly implemented to strictly nonlinear systems. Under dynamic actuator constraints and the availability of only output signals for use in the control law, closed-loop simulations demonstrate that the proposed control techniques are superior to a nonlinear PI control scheme based on the identified Hammerstein model.     

Storing covariance with nonlinearly interacting neurons

Summary A time-dependent, nonlinear model of neuronal interaction which was probabilistically analyzed in a previous article is shown here to be a natural generalization of the Hartline-Ratliff model of the Limulus retina. Although the primary physical variables in the model are the membrane potentials of neurons, the equations which govern the means and covariances of the membrane potentials are coupled through the average firing rates; as a consequence, the average firing rates control the selective storage and retrieval of covariance information. Motor learning in the cerebellar cortex is treated as a problem of covariance storage, and a prediction is made for the underlying synaptic plasticity: the change in synaptic strength between a parallel fiber and a Purkinje cell should be proportional to the covariance between discharges in the parallel fiber and the climbing fiber. Unlike previous proposals for synaptic plasticity, this prediction requires both facilitation and depression to occur (under different conditions) at the same synapse.