The stability of a two-dimensional vorticity filament under uniform strain

The quantitative effects of uniform strain and background rotation on the stability of a strip of constant vorticity (a simple shear layer) are examined. The thickness of the strip decreases in time under the strain, so it is necessary to formulate the linear stability analysis for a time-dependent basic flow. The results show that even a strain rate γ (scaled with the vorticity of the strip) as small as 0.25 suppresses the conventional Rayleigh shear instability mechanism, in the sense that the r.m.s. wave steepness cannot amplify by more than a certain factor, and must eventually decay. For γ < 0.25 the amplification factor increases as γ decreases; however, it is only 3 when γ e 0.065. Numerical simulations confirm the predictions of linear theory at small steepness and predict a threshold value necessary for the formation of coherent vortices. The results help to explain the impression from numerous simulations of two-dimensional turbulence reported in the literature that filaments of vorticity infrequently roll up into vortices. The stabilization effect may be expected to extend to two- and three-dimensional quasi-geostrophic flows.

Interface dynamics in planar neural field models

Neural field models describe the coarse-grained activity of populations of interacting neurons. Because of the laminar structure of real cortical tissue they are often studied in two spatial dimensions, where they are well known to generate rich patterns of spatiotemporal activity. Such patterns have been interpreted in a variety of contexts ranging from the understanding of visual hallucinations to the generation of electroencephalographic signals. Typical patterns include localized solutions in the form of traveling spots, as well as intricate labyrinthine structures. These patterns are naturally defined by the interface between low and high states of neural activity. Here we derive the equations of motion for such interfaces and show, for a Heaviside firing rate, that the normal velocity of an interface is given in terms of a non-local Biot-Savart type interaction over the boundaries of the high activity regions. This exact, but dimensionally reduced, system of equations is solved numerically and shown to be in excellent agreement with the full nonlinear integral equation defining the neural field. We develop a linear stability analysis for the interface dynamics that allows us to understand the mechanisms of pattern formation that arise from instabilities of spots, rings, stripes and fronts. We further show how to analyze neural field models with linear adaptation currents, and determine the conditions for the dynamic instability of spots that can give rise to breathers and traveling waves.

Multiple stationary solutions of an irradiated slab

A mathematical model describing the heat budget of an irradiated medium is introduced. The one-dimensional form of the equations and boundary conditions are presented and analysed. Heat transport at one face of the slab occurs by absorption (and reflection) of an incoming beam of short-wave radiation with a fraction of this radiation penetrating into the body of the slab, a diffusive heat flux in the slab and a prescribed incoming heat flux term. The other face of the slab is immersed in its own melt and is considered to be a free surface. Here, temperature continuity is prescribed and evolution of the surface is determined by a Stefan condition. These boundary conditions are flexible enough to describe a range of situations such as a laser shining on an opaque medium, or the natural environment of polar sea ice or lake ice. A two-stream radiation model is used which replaces the simple Beer’s law of radiation attenuation frequently used for semi-infinite domains. The stationary solutions of the governing equations are sought and it is found that there exists two possible stationary solutions for a given set of boundary conditions and a range of parameter choices. It is found that the existence of two stationary solutions is a direct result of the model of radiation absorption, due to its effect on the albedo of the medium. A linear stability analysis and numerical calculations indicate that where two stationary solutions exist, the solution corresponding to a larger thickness is always stable and the solution corresponding to a smaller thickness is unstable. Numerical simulations reveal that when there are two solutions, if the slab is thinner than the smaller stationary thickness it will melt completely, whereas if the slab is thicker than the smaller stationary thickness it will evolve toward the larger stationary thickness. These results indicate that other mechanisms (e.g. wave-induced agglomeration of crystals) are necessary to grow a slab from zero initial thickness in the parameter regime that yields two stationary solutions.

A mathematical model of crystallization in an emulsion

A mathematical model incorporating many of the important processes at work in the crystallization of emulsions is presented. The model describes nucleation within the discontinuous domain of an emulsion, precipitation in the continuous domain, transport of monomers between the two domains, and formation and subsequent growth of crystals in both domains. The model is formulated as an autonomous system of nonlinear, coupled ordinary differential equations. The description of nucleation and precipitation is based upon the Becker–Döring equations of classical nucleation theory. A particular feature of the model is that the number of particles of all species present is explicitly conserved; this differs from work that employs Arrhenius descriptions of nucleation rate. Since the model includes many physical effects, it is analyzed in stages so that the role of each process may be understood. When precipitation occurs in the continuous domain, the concentration of monomers falls below the equilibrium concentration at the surface of the drops of the discontinuous domain. This leads to a transport of monomers from the drops into the continuous domain that are then incorporated into crystals and nuclei. Since the formation of crystals is irreversible and their subsequent growth inevitable, crystals forming in the continuous domain effectively act as a sink for monomers “sucking” monomers from the drops. In this case, numerical calculations are presented which are consistent with experimental observations. In the case in which critical crystal formation does not occur, the stationary solution is found and a linear stability analysis is performed. Bifurcation diagrams describing the loci of stationary solutions, which may be multiple, are numerically calculated.

A mathematical model of the sterol regulatory element binding protein 2 cholesterol biosynthesis pathway

Cholesterol is one of the key constituents for maintaining the cellular membrane and thus the integrity of the cell itself. In contrast high levels of cholesterol in the blood are known to be a major risk factor in the development of cardiovascular disease. We formulate a deterministic nonlinear ordinary differential equation model of the sterol regulatory element binding protein 2 (SREBP-2) cholesterol genetic regulatory pathway in an hepatocyte. The mathematical model includes a description of genetic transcription by SREBP-2 which is subsequently translated to mRNA leading to the formation of 3-hydroxy-3-methylglutaryl coenzyme A reductase (HMGCR), a main precursor of cholesterol synthesis. Cholesterol synthesis subsequently leads to the regulation of SREBP-2 via a negative feedback formulation. Parameterised with data from the literature, the model is used to understand how SREBP-2 transcription and regulation affects cellular cholesterol concentration. Model stability analysis shows that the only positive steady-state of the system exhibits purely oscillatory, damped oscillatory or monotic behaviour under certain parameter conditions. In light of our findings we postulate how cholesterol homestasis is maintained within the cell and the advantages of our model formulation are discussed with respect to other models of genetic regulation within the literature.

Phytoplankton–zooplankton dynamics in the ‘presence’ or ‘absence’ of toxic phytoplankton

The inhibitory effects of toxin-producing phytoplankton (TPP) on zooplankton modulate the dynamics of marine plankton. In this article, we employ simple mathematical models to compare theoretically the dynamics of phytoplankton–zooplankton interaction in situations where the TPP are present with those where TPP are absent. We consider two sets of three-component interaction models: one that does not include the effect of TPP and the other that does. The negative effects of TPP on zooplankton is described by a non-linear interaction term. Extensive theoretical analyses of the models have been performed to understand the qualitative behaviour of the model systems around every possible equilibria. The results of local-stability analysis and numerical simulations demonstrate that the two model-systems differ qualitatively with regard to oscillations and stability. The model system that does not include TPP is asymptotically stable around the coexisting equilibria, whereas, the system that includes TPP oscillates for a range of parametric values associated with toxin-inhibition rate and competition coefficients. Our analysis suggests that the qualitative dynamics of the plankton–zooplankton interactions are very likely to alter due to the presence of TPP species, and therefore the effects of TPP should be considered carefully while modelling plankton dynamics.

Modeling migratory grazing of zooplankton on toxic and non-toxic phytoplankton

Migratory grazing of zooplankton between non-toxic phytoplankton (NTP) and toxic phytoplankton (TPP) is a realistic phenomena unexplored so far. The present article is a first step in this direction. A mathematical model of NTP–TPP-zooplankton with constant and variable zooplankton migration is proposed and analyzed. The asymptotic dynamics of the model system around the biologically feasible equilibria is explored through local stability analysis. The dynamics of the proposed system is explored and displayed for different combination of migratory parameters and toxin inhibition parameters. Our analysis suggests that the migratory grazing of zooplankton has a significant role in determining the dynamic stability and oscillation of phytoplankton zooplankton systems.

Numerical analyses of Runge–Kutta implicit–explicit schemes for horizontally explicit, vertically implicit solutions of atmospheric models

With the prospect of exascale computing, computational methods requiring only local data become especially attractive. Consequently, the typical domain decomposition of atmospheric models means horizontally-explicit vertically-implicit (HEVI) time-stepping schemes warrant further attention. In this analysis, Runge-Kutta implicit-explicit schemes from the literature are analysed for their stability and accuracy using a von Neumann stability analysis of two linear systems. Attention is paid to the numerical phase to indicate the behaviour of phase and group velocities. Where the analysis is tractable, analytically derived expressions are considered. For more complicated cases, amplification factors have been numerically generated and the associated amplitudes and phase diagnosed. Analysis of a system describing acoustic waves has necessitated attributing the three resultant eigenvalues to the three physical modes of the system. To do so, a series of algorithms has been devised to track the eigenvalues across the frequency space. The result enables analysis of whether the schemes exactly preserve the non-divergent mode; and whether there is evidence of spurious reversal in the direction of group velocities or asymmetry in the damping for the pair of acoustic modes. Frequency ranges that span next-generation high-resolution weather models to coarse-resolution climate models are considered; and a comparison is made of errors accumulated from multiple stability-constrained shorter time-steps from the HEVI scheme with a single integration from a fully implicit scheme over the same time interval. Two schemes, “Trap2(2,3,2)” and “UJ3(1,3,2)”, both already used in atmospheric models, are identified as offering consistently good stability and representation of phase across all the analyses. Furthermore, according to a simple measure of computational cost, “Trap2(2,3,2)” is the least expensive.

Analysis of void volumes in proteins and application to stability of the p53 tumour suppressor protein

We have developed a new method for the analysis of voids in proteins (defined as empty cavities not accessible to solvent). This method combines analysis of individual discrete voids with analysis of packing quality. While these are different aspects of the same effect, they have traditionally been analysed using different approaches. The method has been applied to the calculation of total void volume and maximum void size in a non-redundant set of protein domains and has been used to examine correlations between thermal stability and void size. The tumour-suppressor protein p53 has then been compared with the non-redundant data set to determine whether its low thermal stability results from poor packing. We found that p53 has average packing, but the detrimental effects of some previously unexplained mutations to p53 observed in cancer can be explained by the creation of unusually large voids. (C) 2004 Elsevier Ltd. All rights reserved.

Long-term cultivation-independent microbial diversity analysis demonstrates that bacterial communities infecting the adult cystic fibrosis lung show stability and resilience.

Whilst not true in all cases, the microbial communities that chronically infect the airways of patients with CF can vary little over a year despite antibiotic perturbation. The species present tended to vary more between than within subjects, suggesting that each CF airway infection is unique, with relatively stable and resilient bacterial communities. The inverse relationship between community richness and disease severity is similar to findings reported in other mucosal infections.

Comparative analysis of performance and stability among composite cross populations, variety mixtures and pure lines of winter wheat in organic and conventional cropping systems

This study investigated the effects of increased genetic diversity in winter wheat (Triticum aestivum L.), either from hybridization across genotypes or from physical mixing of lines, on grain yield, grain quality, and yield stability in different cropping environments. Sets of pure lines (no diversity), chosen for high yielding ability or high quality, were compared with line mixtures (intermediate level of diversity), and lines crossed with each other in composite cross populations (CCPn, high diversity). Additional populations containing male sterility genes (CCPms) to increase outcrossing rates were also tested. Grain yield, grain protein content, and protein yield were measured at four sites (two organically-managed and two conventionally-managed) over three years, using seed harvested locally in each preceding year. CCPn and mixtures out-yielded the mean of the parents by 2.4% and 3.6%, respectively. These yield differences were consistent across genetic backgrounds but partly inconsistent across cropping environments and years. Yield stability measured by environmental variance was higher in CCPn and CCPms than the mean of the parents. An index of yield reliability tended to be higher in CCPn, CCPms and mixtures than the mean of the parents. Lin and Binns’ superiority values of yield and protein yield were consistently and significantly lower (i.e. better) in the CCPs than in the mean of the parents, but not different between CCPs and mixtures. However, CCPs showed greater early ground cover and plant height than mixtures. When compared with the (locally non-predictable) best-yielding pure line, CCPs and mixtures exhibited lower mean yield and somewhat lower yield reliability but comparable superiority values. Thus, establishing CCPs from smaller sets of high-performing parent lines might optimize their yielding ability. On the whole, the results demonstrate that using increased within-crop genetic diversity can produce wheat crops with improved yield stability and good yield reliability across variable and unpredictable cropping environments.

Omega-3 fatty acids, cardiovascular disease and stability of atherosclerotic plaques

Long-chain n-3 polyunsaturated fatty acids are found in oily fish and in fish oils and similar preparations. Substantial evidence from epidemiological and case-control studies indicates that consumption of fish, oily fish and long-chain n-3 fatty acids reduces risk of cardiovascular mortality. Secondary prevention studies using long-chain n-3 fatty acids in patients post-myocardial infarction have shown a reduction in total and cardiovascular mortality with an especially potent effect on sudden death. Long-chain n-3 fatty acids have been shown to beneficially modify a range of cardiovascular risk factors, which may result in primary cardiovascular prevention. However, reduced non-fatal and fatal events and a reduction in sudden death probably involve other mechanisms. Reduced thrombosis following long-chain n-3 fatty acids may play a role. A decrease in arrhythmias is a favoured mechanism of action of long-chain n-3 fatty acids and is supported by cell culture and animal studies. However human trials using implantable cardiac defibrillators have produced inconsistent findings and a recent meta-analysis does not support this mechanism of action. An alternative mechanism of action may be stabilisation of atherosclerotic plaques by long-chain n-3 fatty acids. This is suggested by one published human study which showed that incorporation of long-chain n-3 fatty acids into plaques collected at carotid endarterectomy resulted in fewer macrophages in the plaque and a morphology indicative of increased stability. These findings are supported from observations in an animal model and suggest that the primary effect of long-chain n-3 fatty acids might be on macrophages within the plaque.

Association of n-3 polyunsaturated fatty acids with stability of atherosclerotic plaques: a randomised controlled trial

Background N-3 polyunsaturated fatty acids (PUFAs) from oily fish protect against death from cardiovascular disease. We aimed to assess the hypothesis that incorporation of n-3 and n-6 PUFAs into advanced atherosclerotic plaques increases and decreases plaque stability, respectively. Methods We did a randomised controlled trial of patients awaiting carotid endarterectomy. We randomly allocated patients control, sunflower oil (n-6), or fish-oil (n-3) capsules until surgery. Primary outcome was plaque morphology indicative of stability or instability, and outcome measures were concentrations of EPA, DHA, and linoleic acid in carotid plaques; plaque morphology; and presence of macrophages in plaques. Analysis was per protocol. Findings 188 patients were enrolled and randomised; 18 withdrew and eight were excluded. Duration of oil treatment was 7-189 days (median 42) and did not differ between groups. The proportions of EPA and DHA were higher in carotid plaque fractions in patients receiving fish oil compared with those receiving control (absolute difference 0.5 [95% CI 0.3-0.7], 0.4 [0.1-0.6], and 0.2 [0.1-0.4] g/100 g total fatty acids for EPA; and 0.3 [0.0-0.8], 0.4 [0.1-0.7], and 0.3 [0.1-0.6] g/100 g total fatty acids for DHA; in plaque phospholipids, cholesteryl esters, and triacylglycerols, respectively). Sunflower oil had little effect on the fatty acid composition of lipid fractions. Fewer plaques from patients being treated with fish oil had thin fibrous caps and signs of inflammation and more plaques had thick fibrous caps and no signs of inflammation, compared with plaques in patients in the control and sunflower oil groups (odds ratio 0.52 [95% CI 0.24-0.89] and 1.19 [1.02-1.57] vs control; 0.49 [0.23-0.90] and 1.16 [1.01-1.53] vs sunflower oil). The number of macrophages in plaques from patients receiving fish oil was lower than in the other two groups. Carotid plaque morphology and infiltration by macrophages did not differ between control and sunflower oil groups. Interpretation Atherosclerotic plaques readily incorporate n-3 PUFAs from fish-oil supplementation, inducing changes that can enhance stability of atherosclerotic plaques. By contrast, increased consumption of n-6 PUFAs does not affect carotid plaque fatty-acid composition or stability over the time course studied here. Stability of plaques could explain reductions in non-fatal and fatal cardiovascular events associated with increased n-3 PUFA intake.

Enhancing the oxidative stability of rice crackers by addition of the ethanolic extract of phytochemicals from Cratoxylum formosum Dyer

Craloxylum formosum Dyer is consumed throughout the year as food and medicine in Thailand. It contains large amounts of chlorogenic acid and quinic acid derivatives. The antioxidative activity of the extract was studied in refined soybean oil coating on rice crackers without any seasoning. They were stored in accelerated oxidation conditions at 40 degrees C, 80% relative humidity (RH) in the dark for 18 days. The oxidative state of each sample was monitored by analyzing of the peroxide value (PV) and thiobarbituric acid reactive substances (TBARS) as well as by odor analysis by quantitative descriptive analysis (QDA). The C formosum extract was more effective than alpha-tocopherol due to metal ions present in the crackers, which resulted in alpha-tocopherol being less effective as an antioxidant. Sensory odor attributes of rice crackers were related more closely to TBARS than to PV values by linear regression analysis. The present study indicated that C. formosum extract was a promising source of a natural food antioxidant and was effective in inhibiting lipid oxidation in rice crackers.