A linear theory for global electrocortical activity and its control by the lateral hypothalamus

A linear model for electrocortical waves and their control by the lateral hypothalamus is proposed. It is argued that such a linear model is not in contradition to non-linearity of neural elements on the microscopic scale. Telencephalic structures are treated as a mass of linked oscillators generating activity with a number of resonant modes. The lateral hypothalamus is regarded as controlling damping of activity in the telencephalic mass, and therefore exerting a specific parametric control over all signal processing in the cortical networks. An initial test is proposed to assess the constancy of telencephalic natural frequencies, with variation in lateral hypothalamic damping.     

Effects of peristaltic and longitudinal wave motion of the channel wall on movement of micro-organisms: application to spermatozoa transport

A mathematical model is presented to study the motion of the spermatozoa in the cervical canal by considering the transverse waves along its tail and the transverse and longitudinal motions of the cervical wall. In an attempt to control fertility by reducing the speed of sperm, the transverse waves have been considered in the direction opposite to the motion of the spermatozoa. It has been shown that by having appropriate transverse wave motion and longitudinal velocity, the sperm may not be able to move towards the oviduct even if it could continue to have its own propelling velocity. A particular case of the motion of a thin plane sheet in a channel under peristaltic motion of its walls has also been obtained and studied.     

Capturing the Dynamics of a Hybrid Multiscale Cancer Model with a Continuum Model

Cancer is a complex disease involving processes at spatial scales from subcellular, like cell signalling, to tissue scale, such as vascular network formation. A number of multiscale models have been developed to study the dynamics that emerge from the coupling between the intracellular, cellular and tissue scales. Here, we develop a continuum partial differential equation model to capture the dynamics of a particular multiscale model (a hybrid cellular automaton with discrete cells, diffusible factors and an explicit vascular network). The purpose is to test under which circumstances such a continuum model gives equivalent predictions to the original multiscale model, in the knowledge that the system details are known, and differences in model results can be explained in terms of model features (rather than unknown experimental confounding factors). The continuum model qualitatively replicates the dynamics from the multiscale model, with certain discrepancies observed owing to the differences in the modelling of certain processes. The continuum model admits travelling wave solutions for normal tissue growth and tumour invasion, with similar behaviour observed in the multiscale model. However, the continuum model enables us to analyse the spatially homogeneous steady states of the system, and hence to analyse these waves in more detail. We show that the tumour microenvironmental effects from the multiscale model mean that tumour invasion exhibits a so-called pushed wave when the carrying capacity for tumour cell proliferation is less than the total cell density at the tumour wave front. These pushed waves of tumour invasion propagate by triggering apoptosis of normal cells at the wave front. Otherwise, numerical evidence suggests that the wave speed can be predicted from linear analysis about the normal tissue steady state.     

Theory of longitudinal plasma waves with allowance for ion motion

A study is made of the propagation of steady-state large-amplitude longitudinal plasma waves in a cold collisionless plasma with allowance for both electron and ion motion. Conditions for the existence of periodic potential waves are determined. The electric field, potential, frequency, and wavelength are obtained as functions of the wave phase velocity and ion-to-electron mass ratio. Taking into account the ion motion results in the nonmonotonic dependence of the frequency of the waves with the maximum possible amplitudes on the wave phase velocity. Specifically, at low phase velocities, the frequency is equal to the electron plasma frequency for linear waves. As the phase velocity increases, the frequency first decreases insignificantly, reaches its minimum value, and then increases. As the phase velocity increases further, the frequency continues to increase and, at relativistic phase velocities, again becomes equal to the plasma frequency. Finally, as the phase velocity approaches the speed of light, the frequency increases without bound.     

The equation of motion for sperm flagella.

The equation of motion for sperm flagella, in which the elastic bending moment and the active contractile moment are balanced by the moment from the viscous resistance of the surrounding fluid, is solved for a wave solution that superimposes partial solutions. Substitution of the expression for the wave solution into the equation leads to an expression for the active contractile moment. This active moment can be decomposed into two parts. The first part describes an active moment that travels over the flagellum with the mechanical flagellar wave, the second part represents a moment in phase over the entire length of the flagellum, which decreases linearly towards the distal tip. The linear synchronous moment, to which an amount of traveling moment has been added as a perturbation, leads to wave solutions that closely resemble flagellar waves. Properties such as wavelength and wave amplitudes and also the shape of the waves in sea urchin sperm flagella at different frequencies are accurately described by the theory. The change in wave shape in sea urchin sperm flagella at raised viscosity is predicted well by the theory. The different wave properties caused in bull sperm flagella by different boundary conditions at the proximal junction are explained. When only a traveling active moment is present in a flagellum, the wave solutions describe waves of a small wave length in a long flagellum. Some properties of the wave motion of sperm flagella are derived from the theory and verified experimentally.     

Effects of Shear Flows on Nonlinear Waves in Excitable Media

If an excitable medium is moving with relative shear, the waves of excitation may be broken by the motion. We consider such breaks for the case of a constant linear shear flow. The mechanisms and conditions for the breaking of solitary waves and wavetrains are essentially different: the solitary waves require the velocity gradient to exceed a certain threshold, whilst the breaking of repetitive wavetrains happens for arbitrarily small velocity gradients. Since broken waves evolve into new spiral wave sources, this leads to spatio-temporal irregularity.     

Dynamic reconstruction of the fish locomotor wave

Theoretical modeling substantiates the bionic solution of an optimal wave propulsor for water transport, which must have the amplitude and phase characteristics of the fish bending locomotor waves. Use is made of a computer model of multilink chain bending waves where arbitrary distributions can be set for amplitudes and phases of separate segments. With the linear increase in the amplitude of segmental oscillations, a numerical experiment determines the optimal phasing whereby each segment provides a maximal longitudinal component of the motive force. Two versions of hydrodynamic interaction of the segments with water are compared, implying (i) linear and (ii) quadratic drag: the optimal phasing of the transverse and longitudinal oscillations proves to be (i) orthogonal and (ii) nearly orthogonal. The computed bending wave shapes corresponding to dynamic optimization are consistent with the experimental data.     

Nonlinear Planetary Electromagnetic Vortex Structures in the Ionospheric F-Layer

A study is made of the dynamics of planetary-scale electromagnetic waves in the F-layer of the ionosphere. It is shown that, in this layer, a new branch of large-scale magneto-ionospheric wave perturbations is generated under the action of the latitudinal variations of the geomagnetic field, which are a constant property of the ionosphere. The waves propagate along the parallels with phase velocities of tens to hundreds of km/s. The pulsations of the geomagnetic field in the waves can be as strong as several tens of nT. A possible self-localization effect is revealed: these waves may form nonlinear localized solitary vortices moving either westward or eastward along the parallels with velocities much higher than the phase velocities of the linear waves. The characteristic dimension of a vortex is about 10⁴ km or even larger. The magnetic fields generated by vortex structures are one order of magnitude stronger than those in linear waves. The vortices are long-lived formations and may be regarded as elements of strong structural turbulence in the ionosphere. The properties of the wave structures under investigation are very similar to those of ultralow-frequency perturbations observed experimentally in the ionosphere at middle latitudes.     

A test for constant natural frequencies in electrocortical activity under lateral hypothalamic control

An initial test for a theory of lateral hypothalamic regulation of electrocortical activity is undertaken. The theory supposes lateral hypothalamic input directly or indirectly damps telencephalic resonances involving linear wave phenomena, enabling this pathway to act as parametric control of information processing in cortical neural networks. Relative changes in left and right electrocortical power spectra are used to test for the presence of resonant modes with constant natural frequencies in conditions of asymmetrical damping, following unilateral lesion of the lateral hypothalamus. Natural frequency values for the modes clustered about center frequencies in the EEG band are obtained. This method has the advantage of minimising the effects of time-variation and the recorded signal's distortion from the electrocortical local spatial average, but limits consideration to five dominant modes of resonance. The uncertainty of true model order, and errors in curve-fitting impose limitations on the test.     

Flutter in collapsible tubes: a theoretical model of wheezes

A mathematical analysis of flow through a flexible channel is examined as a model of flow-induced flutter oscillations that pertain to the production of wheezing breath sounds. The model provides predictions for the critical fluid speed that will initiate flutter waves of the wall, as well as their frequency and wavelength. The mathematical results are separated into linear theory (small oscillations) and nonlinear theory (larger oscillations). Linear theory determines the onset of the flutter, whereas nonlinear theory determines the relationships between the fluid speed and both the wave amplitudes and frequencies. The linear theory predictions correlate well with data taken at the onset of flutter and flow limitation during experiments of airflow in thick-walled collapsible tubes. The nonlinear theory predictions correlate well with data taken as these flows are forced to higher velocities while keeping the flow rate constant. Particular ranges of the parameters are selected to investigate and discuss the applications to airway flows. According to this theory, the mechanism of generation of wheezes is based in the interactions of fluid forces and friction and wall elastic-restoring forces and damping. In particular, a phase delay between the fluid pressure and wall motion is necessary. The wave speed theory of flow limitation is discussed with respect to the specific data and the flutter model.     

Formation and evolution of flapping and ballooning waves in magnetospheric plasma sheet

By adopting Lembége & Pellat’s 2D plasma-sheet model, we investigate the flankward flapping motion and Sunward ballooning propagation driven by an external source (e.g., magnetic reconnection) produced initially at the sheet center. Within the ideal MHD framework, we adopt the WKB approximation to obtain the Taylor–Goldstein equation of magnetic perturbations. Fourier spectral method and Runge–Kutta method are employed in numerical simulations, respectively, under the flapping and ballooning conditions. Studies expose that the magnetic shears in the sheet are responsible for the flapping waves, while the magnetic curvature and the plasma gradient are responsible for the ballooning waves. In addition, the flapping motion has three phases in its temporal development: fast damping phase, slow recovery phase, and quasi-stabilized phase; it is also characterized by two patterns in space: propagating wave pattern and standing wave pattern. Moreover, the ballooning modes are gradually damped toward the Earth, with a wavelength in a scale size of magnetic curvature or plasma inhomogeneity, only 1–7% of the flapping one; the envelops of the ballooning waves are similar to that of the observed bursty bulk flows moving toward the Earth.     

Separation of arterial pressure into a nonlinear superposition of solitary waves and a windkessel flow

A simplified model of arterial blood pressure intended for use in model-based signal processing applications is presented. The main idea is to decompose the pressure into two components: a travelling wave which describes the fast propagation phenomena predominating during the systolic phase and a windkessel flow that represents the slow phenomena during the diastolic phase. Instead of decomposing the blood pressure pulse into a linear superposition of forward and backward harmonic waves, as in the linear wave theory, a nonlinear superposition of travelling waves matched to a reduced physical model of the pressure, is proposed. Very satisfactory experimental results are obtained by using forward waves, the N-soliton solutions of a Korteweg–de Vries equation in conjunction with a two-element windkessel model. The parameter identifiability in the practically important 3-soliton case is also studied. The proposed approach is briefly compared with the linear one and its possible clinical relevance is discussed.     

Stimulus motion propels traveling waves in binocular rivalry

State transitions in the nervous system often take shape as traveling waves, whereby one neural state is replaced by another across space in a wave-like manner. In visual perception, transitions between the two mutually exclusive percepts that alternate when the two eyes view conflicting stimuli (binocular rivalry) may also take shape as traveling waves. The properties of these waves point to a neural substrate of binocular rivalry alternations that have the hallmark signs of lower cortical areas. In a series of experiments, we show a potent interaction between traveling waves in binocular rivalry and stimulus motion. The course of the traveling wave is biased in the motion direction of the suppressed stimulus that gains dominance by means of the wave-like transition. Thus, stimulus motion may propel the traveling wave across the stimulus to the extent that the stimulus motion dictates the traveling wave's direction completely. Using a computational model, we show that a speed-dependent asymmetry in lateral inhibitory connections between retinotopically organized and motion-sensitive neurons can explain our results. We argue that such a change in suppressive connections may play a vital role in the resolution of dynamic occlusion situations.     

Storm Resistant Boat Designing Based on the Geometry and Movement of Water Strider

A novel conceptual model of boat is developed based on the bionic properties of water strider. This four-legged model is called storm resistance boat. Concentrated balancing forces, asymmetry of arms and smooth motion against waves and the effect of arms length on the reduction of drag and effects of wave are the major characteristics of water strider which are considered for designing this model. The results indicate the characteristic improvement of the small boat and its resistance against strong waves as well as marine ill-conditions. This boat can be considered as a high speed rescue boat in marine traffic.     

Early development and quorum sensing in bacterial biofilms

 We develop mathematical models to examine the formation, growth and quorum sensing activity of bacterial biofilms. The growth aspects of the model are based on the assumption of a continuum of bacterial cells whose growth generates movement, within the developing biofilm, described by a velocity field. A model proposed in Ward et al. (2001) to describe quorum sensing, a process by which bacteria monitor their own population density by the use of quorum sensing molecules (QSMs), is coupled with the growth model. The resulting system of nonlinear partial differential equations is solved numerically, revealing results which are qualitatively consistent with experimental ones. Analytical solutions derived by assuming uniform initial conditions demonstrate that, for large time, a biofilm grows algebraically with time; criteria for linear growth of the biofilm biomass, consistent with experimental data, are established. The analysis reveals, for a biologically realistic limit, the existence of a bifurcation between non-active and active quorum sensing in the biofilm. The model also predicts that travelling waves of quorum sensing behaviour can occur within a certain time frame; while the travelling wave analysis reveals a range of possible travelling wave speeds, numerical solutions suggest that the minimum wave speed, determined by linearisation, is realised for a wide class of initial conditions. Received: 10 February 2002 / Revised version: 29 October 2002 / Published online: 19 March 2003 Key words or phrases: Bacterial biofilm – Quorum sensing – Mathematical modelling – Numerical solution – Asymptotic analysis – Travelling wave analysis     

Instabilities of a circularly polarized wave with trapped particles in an isotropic plasma

The structure and stability of a transverse electromagnetic wave propagating with a velocity lower than the speed of light in an unmagnetized plasma are considered. The stationary finite-amplitude wave is described by exact solutions to the Vlasov-Maxwell equations. However, unlike the well-known electrostatic analog, the Bernstein-Greene-Kruskal wave, the wave structure is determined to a large extent by the presence of trapped particles with a shear of transverse velocities, without which the existence of waves with a refraction index larger than unity is impossible. It is shown that the main origin of the wave instability is the longitudinal motion of trapped particles relative to the background plasma. Expressions for the growth rates in the main instability regimes are found under definite restrictions on the wave parameters.     

High-frequency surface waves at a plasma-metal interface: I. Linear model

A study is made of the dispersion properties of surface waves at a plasma-metal interface under thermodynamically nonequilibrium conditions such that a space charge sheath forms at the plasma boundary. In the simplest model, the sheath is described as a dielectric with a given permittivity. The wave parameters in a highly collisional plasma are discussed. The effect of interaction between waves propagating near the opposite plasma boundaries is considered, in particular, for space charge sheaths of different thicknesses. Conditions are determined under which the parameters of surface waves are substantially altered by the plasma-sheath geometric resonance.     

Numerical simulation of circulation in a Caribbean-type backreef lagoon

A numerical model based on the vertically intergrated equations of motion and continuity was used to simulate circulation in a shallow, well-mixed, Caribbean-type backreef lagoon. As the focus of previous research, Great Pond Bay, St. Croix, provided a suitable study site. Tides, wind, and the effect of waves impinging on the reef were incorporated within the model. Direction of simulated flow under various wind and wave regimes agrees well with patterns found during current meter and drogue tracking experiments. When momentum transfer over the reef due to wave breaking and set-up is added, the magnitude of flow within the lagoon, ranging 5–30 cm/s, also compares favorably to in situ measurements. Model results indicate a relatively rapid and realistic flushing period of less than one day. Results of this study indicate that numerical modeling techniques potentially offer an accurate and cost-effective means to predict the pattern of water movement within coral-reef lagoons.     

Excitation of growing wake waves by an electron bunch in a plasma in an intense pump wave

A study is made of the excitation of wake waves by a one-dimensional electron bunch in an electron plasma in the presence of an intense monochromatic pump wave with circular polarization. In the main state (in the absence of a bunch), the interaction between a pump wave and a plasma is described by Maxwell's equations and the nonlinear relativistic hydrodynamic equations for a cold plasma. The excitation of linear waves by a one-dimensional bunch is investigated against a cold plasma background. It is shown that, in a certain range of parameter values of the bunch, pump wave, and plasma, the excitation is resonant in character and the amplitude of the excited wake waves increases with distance from the bunch.     

The movement of spermatozoa with helical head: theoretical analysis and experimental results.

The present work is concerned with the study of the swimming of flagellated microscopic organisms with a helical head and a helical pattern of flagellar beating, such as Xenopus sperms. The theoretical approach is similar to that taken by Chang and Wu (1971) in the study of helical flagellar movement. The model used in the present study allows us to determine the velocity of propulsion (U) and the frequency of rotation of the sperm head (fh) as a function of the frequency of the wave of motion (ft) traveling along the tail. The results relative to the case of helical and planar flagellar waves are compared. Our main finding is that the helical shape of the head seems to increase the efficiency of propulsion of the spermatozoon when compared with the more commonly shaped spherical head. Experimentally measured values of fh versus U may be fitted by a linear plot whose slope is much higher than that corresponding to the case of planar flagellar beating. This fact is consistent with an effectively three-dimensional (nonplanar) movement of the flagellar tail. However, the results do not fit those predicted from a circular helix, suggesting that a different shape of the flagellar beating should be considered.