Material constants for a finite element model of the intervertebral disk with a fiber composite annulus

A simple axisymmetric finite element model of a human spine segment containing two adjacent vertebrae and the intervening intervertebral disk was constructed. The model incorporated four substructures: one to represent each of the vertebral bodies, the annulus fibrosus, and the nucleus pulposus. A semi-analytic technique was used to maintain the computational economies of a two-dimensional analysis when nonaxisymmetric loads were imposed on the model. The annulus material was represented as a layered fiber-reinforced composite. This paper describes the selection of material constants to represent the anisotropic layers of the annulus. It shows that a single set of material constants can be chosen so that model predictions of gross disk behavior under compression, torsion, shear, and moment loading are in reasonable agreement with the mean and range of experimentally measured disk behaviors. It also examines the effects of varying annular material properties.     

Effect of nonaxisymmetric hematocrit distribution on non-Newtonian blood flow in small tubes

Hematocrit distribution and red blood cell aggregation are the major determinants of blood flow in narrow tubes at low flow rates. It has been observed experimentally that in microcirculation the hematocrit distribution is not uniform. This nonuniformity may result from plasma skimming and cell screening effects and also from red cell sedimentation. The goal of the present study is to understand the effect of nonaxisymmetric hematocrit distribution on the flow of human and cat blood in small blood vessels of the microcirculation. Blood vessels are modeled as circular cylindrical tubes. Human blood is described by Quemada's rheological model, in which local viscosity is a function of both the local hematocrit and a structural parameter that is related to the size of red blood cell aggregates. Cat blood is described by Casson's model. Eccentric hematocrit distribution is considered such that the axis of the cylindrical core region of red cell suspension is parallel to the axis of the blood vessel but not coincident. The problem is solved numerically by using finite element method. The calculations predict nonaxisymmetric distribution of velocity and shear stress in the blood vessel and the increase of apparent viscosity with increasing eccentricity of the core.     

Progress in theory of instabilities in a rotating plasma

A review is given of the basic results of modern theory of instabilities in a rotating plasma. Both axisymmetric and nonaxisymmetric perturbations are considered. Main attention is given to the magnetorotational instability (MRI), discovered earlier by Velikhov, and the rotational-convective instability (RCI) discussed in a number of papers of astrophysical trend. For qualitative explanation of the results, a local approach is used which, with equilibrium plasma pressure gradient and/or nonsymmetry of perturbations, requires operation with nonlocal azimuthal perturbed magnetic field. The gravity and effects of pressure anisotropy are taken into account. In addition to hydrodynamic, the electrodynamic approach is formulated. The drift effects are considered. Analyzed are the ideal instabilities and those depending on the dissipative effects: viscosity and heat conductivity. The MRI is considered at presence of the charged dust particles. Besides the local approach, the nonlocal approach is formulated for the plasma model with a steplike profile of angular rotation frequency. Alongside with perturbations which frequencies are small compared to the ion cyclotron frequency, the perturbations are analyzed with frequencies larger than the ion cyclotron frequency. The latter corresponds to the Hall regime and subregime of nonmagnetized plasma.     

Steady propagation of a liquid plug in a two-dimensional channel

In this study, we investigate the steady propagation of a liquid plug within a two-dimensional channel lined by a uniform, thin liquid film. The Navier-Stokes equations with free-surface boundary conditions are solved using the finite volume numerical scheme. We examine the effect of varying plug propagation speed and plug length in both the Stokes flow limit and for finite Reynolds number (Re). For a fixed plug length, the trailing film thickness increases with plug propagation speed. If the plug length is greater than the channel width, the trailing film thickness agrees with previous theories for semi-infinite bubble propagation. As the plug length decreases below the channel width, the trailing film thickness decreases, and for finite Re there is significant interaction between the leading and trailing menisci and their local flow effects. A recirculation flow forms inside the plug core and is skewed towards the rear meniscus as Re increases. The recirculation velocity between both tips decreases with the plug length. The macroscopic pressure gradient, which is the pressure drop between the leading and trailing gas phases divided by the plug length, is a function of U and U2, where U is the plug propagation speed, when the fluid property and the channel geometry are fixed. The U2 term becomes dominant at small values of the plug length. A capillary wave develops at the front meniscus, with an amplitude that increases with Re, and this causes large local changes in wall shear stresses and pressures.     

Wave Propagation through a Viscous Fluid Contained in a Tethered, Initially Stressed, Orthotropic Elastic Tube

To give a realistic representation of the pulse propagation in arteries a theoretical analysis of the wave propagation through a viscous incompressible fluid contained in an initially stressed elastic tube is considered. The tube is assumed to be orthotropic and its longitudinal motion is constrained by a uniformly distributed additional mass, a dashpot and a spring. The fluid is assumed to be Newtonian. The analysis is restricted to propagation of small amplitude harmonic waves whose wavelength is large compared to the radius of the vessel. Elimination of arbitrary constants from the general solutions of the equations of motion of the fluid and the wall gives a frequency equation to determine the velocity of propagation. Two roots of this equation give the velocity of propagation of two distinct outgoing waves. One of the waves propagates slower than the other. The propagation properties of s lower waves are very slightly affected by the degree of anisotropy of the wall. The velocity of propagation of faster waves decreases as the ratio of the longitudinal modulus of elasticity to the circumferential modulus decreases; transmission of these waves is very little affected. The influence of the tethering on the propagation velocity of slower waves is negligibly small; transmission of these waves is seriously affected. In tethered tubes faster waves are completely attenuated.     

Preemptive spatial competition under a reproduction-mortality constraint

Spatially structured ecological interactions can shape selection pressures experienced by a population's different phenotypes. We study spatial competition between phenotypes subject to antagonistic pleiotropy between reproductive effort and mortality rate. The constraint we invoke reflects a previous life-history analysis; the implied dependence indicates that although propagation and mortality rates both vary, their ratio is fixed. We develop a stochastic invasion approximation predicting that phenotypes with higher propagation rates will invade an empty environment (no biotic resistance) faster, despite their higher mortality rate. However, once population density approaches demographic equilibrium, phenotypes with lower mortality are favored, despite their lower propagation rate. We conducted a set of pairwise invasion analyses by simulating an individual-based model of preemptive competition. In each case, the phenotype with the lowest mortality rate and (via antagonistic pleiotropy) the lowest propagation rate qualified as evolutionarily stable among strategies simulated. This result, for a fixed propagation to mortality ratio, suggests that a selective response to spatial competition can extend the time scale of the population's dynamics, which in turn decelerates phenotypic evolution.     

Experimental validation of a time-domain-based wave propagation model of blood flow in viscoelastic vessels

Time-domain-based one-dimensional wave propagation models of the arterial system are preferable over one-dimensional wave propagation models in the frequency domain since the latter neglect the non-linear convection forces present in the physiological situation, especially when the vessel is tapered. Moreover, one-dimensional wave propagation models of the arterial system can be used to provide boundary conditions for fully three-dimensional fluid-structure interaction computations that are usually defined in the time domain. In this study, a time-domain-based one-dimensional wave propagation model in a cross-sectional area, flow and pressure (A,q,p)-formulation is developed. Using this formulation, a constitutive law that includes viscoelasticity based on the mechanical behaviour of a Kelvin body, is introduced. The resulting pressure and flow waves travelling through a straight and tapered vessel are compared to experimental data obtained from measurements in an in vitro setup. The model presented shows to be well suited to predict wave propagation through these straight and tapered vessels with viscoelastic wall properties and hereto can serve as a time-domain-based method to model wave propagation in the human arterial system.     

Study using the Noble model of excitation propagation through a branching node]

Excitation propagation through a branching node with 3 thin branches fused into a thick one has been studied on Noble model. When the ratio of the diameters of thin and thick branches equals 1 : 10 the propagation is blocked. When the diameter ratio is 1 : 8 and excitation comes asynchroneously (along one of input branches 10 msec later than along two others) a 30 msec delay of propagation is observed. The importance of the results obtained for explaining atrio-ventricular delay is discussed.     

Wave Propagation through a Newtonian Fluid Contained within a Thick-Walled,Viscoelastic Tube

The propagation of harmonic pressure waves through a Newtonian fluid contained within a thick-walled, viscoelastic tube is considered as a model of arterial blood flow. The fluid is assumed to be homogeneous and Newtonian, and its motion to be laminar and axisymmetric. The wall is assumed to be isotropic, incompressible, linear, and viscoelastic. It is also assumed that the motion is such that the convective acceleration is negligible. The motion of the fluid is described by the linearized form of the Navier-Stokes equations and the motion of the wall by classical elasticity theory. The frequency dependence of the wall mechanical properties are represented by a three parameter, relaxation-type model. Using boundary conditions describing the continuity of stress and velocity components in the fluid and the wall, explicit solutions for the system of equations of the model have been obtained. The longitudinal fluid impedance has been expressed in terms of frequency and the system parameters. The frequency equation has been solved and the propagation constant also expressed in terms of frequency and system parameters. The results indicate that the fluid impedance is smaller than predicted by the rigid tube model or by Womersley''s constrained elastic tube model. Also, the velocity of propagation is generally slower and the transmission per wavelength less than predicted by Womersley''s elastic tube model. The propagation constant is very sensitive to changes in the degree of wall viscoelasticity.     

Conditions for Loss Compensation of Surface Plasmon Polaritons Propagation on a Metal/Gain Medium Boundary

One way to compensate for the surface plasmon polariton (SPP) propagation losses is to use a gain medium. However, simply ensuring high enough gain is not sufficient because it may violate the bounded character of the wave. Therefore, a detailed theoretical analysis is needed for the determination of the conditions for lossless or amplified SPP propagation. Here presented is an exact theoretical analysis of the SPP propagation in the case of an infinite metal/gain medium boundary. It is shown that the conditions for lossless/amplified SPP propagation can be conveniently examined and presented as a simply connected region in the complex plane of the gain medium dielectric function. Effective and minimum gain parameters are introduced, which facilitates the simultaneous analyses of different gain media/metals combinations. The practical application of these results is illustrated for several gain media/metal (silver, gold and aluminium) systems.     

Epileptiform activity in a neocortical network: a mathematical model

 A simple mathematical model describing the generation and propagation of epileptiform activity in a cerebral cortical network is presented. The model consists of a system of nonlinear delay differential equations. Physiological properties are taken into account as nonlinear transmission of signals at the synapse, temporal and spatial summation of incoming signals at the soma, active membrane characteristics, and dendritic and axonal propagation times. The influence of the connectivity and the temporal parameters on the oscillatory properties of the model is studied. The computer simulations are in agreement with experimental observations in cortical networks: whereas a weak excitatory or strong inhibitory synaptic connection strength produces a stationary status with short-lasting responses to external stimuli, increases in excitation or decreases in inhibition induce spontaneous and stimulus-evoked rhythmic discharges. Synaptic burst-like activity is observed only for an intermediate range of excitatory and inhibitory connection strengths and external inputs. The form and duration of the bursts can also be controlled by the temporal parameters. The results demonstrate that relatively simple mathematical equations are sufficient to model some of the network properties underlying the generation and propagation of epileptiform activity. Received: 2 October 2000 / Accepted in revised form: 4 March 2001     

圆红冬孢酵母两阶段培养法生产微生物油脂

为缩短发酵时间,减少原料消耗,采用细胞增殖和油脂积累分离的两阶段模式,培养圆红冬孢酵母Rhodosporidium toruloides AS 2.1389生产微生物油脂。结果表明,细胞增殖阶段获得的R.toruloides AS 2.1389细胞,重悬接种在葡萄糖溶液中,可快速积累油脂,菌体油脂含量超过自身干重的55%。增殖阶段细胞的菌龄越高,产油能力越强。油脂积累阶段使用高浓度葡萄糖溶液或未灭菌的葡萄糖溶液,油脂合成可以有效进行。油脂中脂肪酸成分以含16和18个碳原子的长链脂肪酸为主,可作为制备生物柴油的新型原料。     

Electrical impedance tomography imaging using a priori ultrasound data

Transverse propagation was previously found to occur in a two-dimensional model of cardiac muscle using the PSpice software program for electronic circuit design and analysis. Longitudinal propagation within each chain, and transverse propagation between parallel chains, occurred even when there were no gap-junction (g-j) channels inserted between the simulated myocardial cells either longitudinally or transversely. In those studies, there were pronounced edge (boundary) effects and end-effects even within single chains. Transverse velocity increased with increase in model size. The present study was performed to examine boundary effects on transverse propagation velocity when the length of the chains was held constant at 10 cells and the number of parallel chains was varied from 3 to 5, to 7, to 10, and to 20. The number of g-j channels was either zero, both longitudinally and transversely (0/0), or 100/100. Some experiments were also made at 100/0, 1/1, and 10/10. Transverse velocity and overall velocity (both longitudinal and transverse components) was calculated from the measured total propagation time (TPT), i.e., the elapsed time between when the first action potential (AP) and the last AP crossed the zero potential level. The transverse g-j channels were placed only at the ends of each chain, such that propagation would occur in a zigzag pattern. Electrical stimulation was applied intracellularly between cells A1 and A2. It was found that, with no g-j channels (0/0), overall velocity increased almost linearly when more and more chains were placed in parallel. In contrast, with many g-j channels (100/100), there was a much flatter relationship between overall velocity and number of parallel chains. The difference in velocities with 0/0 channels and 100/100 channels was reduced as the number of chains was increased. In conclusion, edges have important effects on propagation velocity (overall and transverse) in cardiac muscle simulations.     

Motion of nonaxisymmetric red blood cells in cylindrical capillaries

We analyze theoretically the single-file flow of asymmetric red blood cells along cylindrical capillaries. Red cells in narrow capillaries are typically nonaxisymmetric, with the cell membrane moving continuously around the cell. In our analysis, cell shape and streamlines of membrane motion are prescribed. Lubrication theory is used to compute velocities and pressures in the fluid surrounding the cell. Conditions of zero lift, zero torque, zero drag, and energy conservation in the cell are imposed. Predicted tank-treading frequency, cell inclination and transverse displacement are small. Cell asymmetry and tank-treading are found to have little effect on the apparent viscosity of blood in capillaries with diameters up to 7 microns.     

Arterial pulse attenuation prediction using the decaying rate of a pressure wave in a viscoelastic material model

The present study examines the possibility of attenuating blood pulses by means of introducing prosthetic viscoelastic materials able to absorb energy and damp such pulses. Vascular prostheses made of polymeric materials modify the mechanical properties of blood vessels. The effect of these materials on the blood pulse propagation remains to be fully understood. Several materials for medical applications, such as medical polydimethylsiloxane or polytetrafluoroethylene, show viscoelastic behavior, modifying the original vessel stiffness and affecting the propagation of blood pulses. This study focuses on the propagation of pressure waves along a pipe with viscoelastic materials using the Maxwell and the Zener models. An expression of exponential decay has been obtained for the Maxwell material model and also for low viscous coefficient values in the Zener model. For relatively high values of the viscous term in the Zener model, the steepest part of the pulse can be damped quickly, leaving a smooth, slowly decaying wave. These mathematical models are critical to tailor those materials used in cardiovascular implants to the mechanical environment they are confronted with to repair or improve blood vessel function.     

Scalar equation for wave beams in a magnetized plasma

A heuristic procedure is proposed for deriving quasi-optical equations for wave beams in anisotropic gyrotropic media with allowance for aberrations, spatial dispersion, and absorption. To solve such equations numerically, a method is developed that generalizes the operator exponent method. The applicability limits of the aberration-free approximation for simulating the propagation of beams in absorbing media are determined. Numerical examples of the propagation of beams in the vicinity of the electron cyclotron resonance in plasmas in actual devices are presented.     

Theoretical analysis of the magnetocardiographic pattern for reentry wave propagation in a three-dimensional human heart model

We present a computational study of reentry wave propagation using electrophysiological models of human cardiac cells and the associated magnetic field map of a human heart. We examined the details of magnetic field variation and related physiological parameters for reentry waves in two-dimensional (2-D) human atrial tissue and a three-dimensional (3-D) human ventricle model. A 3-D mesh system representing the human ventricle was reconstructed from the surface geometry of a human heart. We used existing human cardiac cell models to simulate action potential (AP) propagation in atrial tissue and 3-D ventricular geometry, and a finite element method and the Galerkin approximation to discretize the 3-D domain spatially. The reentry wave was generated using an S1-S2 protocol. The calculations of the magnetic field pattern assumed a horizontally layered conductor for reentry wave propagation in the 3-D ventricle. We also compared the AP and magnetocardiograph (MCG) magnitudes during reentry wave propagation to those during normal wave propagation. The temporal changes in the reentry wave motion and magnetic field map patterns were also analyzed using two well-known MCG parameters: the current dipole direction and strength. The current vector in a reentry wave forms a rotating spiral. We delineated the magnetic field using the changes in the vector angle during a reentry wave, demonstrating that the MCG pattern can be helpful for theoretical analysis of reentry waves.     

Female sterility in Ulmus minor (Ulmaceae): a hypothesis invoking the cost of sex in a clonal plant

A high incidence of individuals with low seed set was found in two populations of the field elm Ulmus minor, a European tree that reproduces sexually and via vegetative propagation through root sprouting. One population was a seminatural stand, while the other was established by artificial propagation of genotypes sampled widely across Spain. The low seed set in both populations was due to both pre- and post-zygotic factors, the importance of which vary between genotypes. These factors included gynoecial malformations that produced a non-ovulated pistil, early gynoecial necrosis (i.e., necrosis before any opportunities for pollination), and seed abortion. Female sterility gave rise to two classes of individuals: trees that were largely female-sterile but dispersed normal quantities of viable pollen, and trees that dispersed both normal pollen and substantial numbers of seeds. Reduced production of protein-rich seeds may increase the resource availability for clonal propagation, helping to maintain female-sterile individuals with hermaphrodites.     

Replication of viruses in a growing plaque: a reaction-diffusion model.

An understanding of the viral replication process commonly referred to as "plaque growth" is developed in the context of a reaction-diffusion model. The interactions among three components: the virus, the healthy host, and the infected host are represented using rates of viral adsorption and desorption to the cell surface, replication and release by host lysis, and diffusion. The solution to the full model reveals a maximum in the dependence of the velocity of viral propagation on its equilibrium adsorption constant, suggesting that conditions can be chosen where viruses which adsorb poorly to their hosts will replicate faster in plaques than those which adsorb well. Analytic expressions for the propagation velocity as a function of the kinetic and diffusion parameters are presented for the limiting cases of equilibrated adsorption, slow adsorption, fast adsorption, and large virus yields. Hindered diffusion at high host concentrations must be included for quantitative agreement with experimental data.     

Experimental validation of a microcracking-based toughening mechanism for cortical bone

It has been proposed that cortical bone derives its toughness by forming microcracks during the process of crack propagation (J. Biomech. 30 (1997) 763; J. Biomech. 33 (2000) 1169). The purpose of this study was to experimentally validate the previously proposed microcrack-based toughening mechanism in cortical bone. Crack initiation and propagation tests were conducted on cortical bone compact tension specimens obtained from the antlers of red deer. For these tests, the main fracture crack was either propagated to a predetermined crack length or was stopped immediately after initiating from the notch. The microcracks produced in both groups of specimens were counted in the same surface area of interest around and below the notch, and crack growth resistance and crack propagation velocity were analyzed. There were more microcracks in the surface area of interest in the propagation than in initiation specimens showing that the formation of microcracks continued after the initiation of a fracture crack. Crack growth resistance increased with crack extension, and crack propagation velocity vs. crack extension curves demonstrated the characteristic jump increase and decrease pattern associated with the formation of microcracks. The scanning electron micrographs of crack initiation and propagation displayed the formation of a frontal process zone and a wake, respectively. These results support the microcrack-based toughening mechanism in cortical bone. Bone toughness is, therefore, determined by its ability to form microcracks during fracture.