The development and stability analysis of a nonlinear growth model for microorganisms

A nonlinear dynamic model of microbial growth is established based on the theories of the diffusion response of thermodynamics and the chemotactic response of biology. Except for the two traditional variables, i.e. the density of bacteria and the concentration of attractant, the pH value, a crucial influencing factor to the microbial growth, is also considered in this model. The pH effect on the microbial growth is taken as a Gaussian function G0e-(f- fc)2/G1, where G0, G1 and fc are constants, f represents the pH value and fc represents the critical pH value that best fits for microbial growth. To study the effects of the reproduction rate of the bacteria and the pH value on the stability of the system, three parameters a, G0 and G1 are studied in detail, where a denotes the reproduction rate of the bacteria, G0 denotes the impacting intensity of the pH value to microbial growth and G1 denotes the bacterial adaptability to the pH value. When the effect of the pH value of the solution which microorganisms live in is ignored in the governing equations of the model, the microbial system is more stable with larger a. When the effect of the bacterial chemotaxis is ignored, the microbial system is more stable with the larger G1 and more unstable with the larger G0 for f0 > fc. However, the stability of the microbial system is almost unaffected by the variation G0 and G1 and it is always stable for f0 < fc under the assumed conditions in this paper. In the whole system model, it is more unstable with larger G1 and more stable with larger G0 for f0 < fc. The system is more stable with larger G1 and more unstable with larger G0 for f0 > fc. However, the system is more unstable with larger a for f0 < fc and the stability of the system is almost unaffected by a for f0 > fc. The results obtained in this study provide a biophysical insight into the understanding of the growth and stability behavior of microorganisms.

Linear stability of flow past two circular cylinders in a side-by-side arrangement

The linear stability of flow past two circular cylinders in a side-by-side arrangement is investigated theoretically, numerically and experimentally under the assumption of a two-dimensional flow field, in order to explore the origin of in-phase and antiphase oscillatory flows. Steady symmetric flow is realized at a small Reynolds number, but becomes unstable above a critical Reynolds number though the solution corresponding to the flow still satisfies the basic equations irrespective of the magnitude of the Reynolds number. We obtained the solution numerically and investigated its linear stability. We found that there are two kinds of unstable modes, i.e., antisymmetric and symmetric modes, which lead to in-phase and antiphase oscillatory flows, respectively. We determined the critical Reynolds numbers for the two modes and evaluated the critical distance at which the most unstable disturbance changes from the antisymmetric to the symmetric mode, or vice versa. ©2005 The Physical Society of Japan.

Stability analysis of tree structured decision functions

In multicriteria decision problems many values must be assigned, such as the importance of the different criteria and the values of the alternatives with respect to subjective criteria. Since these assignments are approximate, it is very important to analyze the sensitivity of results when small modifications of the assignments are made. When solving a multicriteria decision problem, it is desirable to choose a decision function that leads to a solution as stable as possible. We propose here a method based on genetic programming that produces better decision functions than the commonly used ones. The theoretical expectations are validated by case studies. © 2003 Elsevier B.V. All rights reserved.

Recent Results on Stability Analysis of an Optimal Assembly Line Balance

Two assembly line balancing problems are addressed. The first problem (called SALBP-1) is to minimize number of linearly ordered stations for processing n partially ordered operations V = {1, 2, ..., n} within the fixed cycle time c. The second problem (called SALBP-2) is to minimize cycle time for processing partially ordered operations V on the fixed set of m linearly ordered stations. The processing time ti of each operation i ∈V is known before solving problems SALBP-1 and SALBP-2. However, during the life cycle of the assembly line the values ti are definitely fixed only for the subset of automated operations V\V . Another subset V ⊆ V includes manual operations, for which it is impossible to fix exact processing times during the whole life cycle of the assembly line. If j ∈V , then operation times tj can differ for different cycles of the production process. For the optimal line balance b of the assembly line with operation times t1, t2, ..., tn, we investigate stability of its optimality with respect to possible variations of the processing times tj of the manual operations j ∈ V .

Stability of membrane bound reactions

We present a novel approach to the dynamics of reactions of diffusing chemical species with species fixed in space e.g. by binding to a membrane. The non-diffusing reaction partners are clustered in areas with a diameter smaller than the diffusion length of the diffusing partner. The activated fraction of the fixed species determines the size of an active sub-region of the cluster. Linear stability analysis reveals that diffusion is one of the ma jor determinants of the stability of the dynamics. We illustrate the model concept with Ca²⁺ dynamics in living cells, which has release channels as fixed reaction partners. Our results suggest that spatial and temporal structures in intracellular Ca²⁺ dynamics are caused by fluctuations due to the small number of channels per cluster.

SLOPE STABILITY ANALYSIS OF MOUNT MEAGER, SOUTH-WESTERN BRITISH COLUMBIA, CANADA

The Mount Meager Volcanic Complex (MMVC) in south-western British Columbia is a potentially active, hydrothermally altered massif comprising a series of steep, glaciated peaks. Climatic conditions and glacial retreat has led to the further weathering, exposure and de-buttressing of steep slopes composed of weak, unconsolidated material. This has resulted in an increased frequency of landslide events over the past few decades, many of which have dammed the rivers bordering the Complex. The breach of these debris dams presents a risk of flooding to the downstream communities. Preliminary mapping showed there are numerous sites around the Complex where future failure could occur. Some of these areas are currently undergoing progressive slope movement and display features to support this such as anti-scarps and tension cracks. The effect of water infiltration on stability was modelled using the Rocscience program Slide 6.0. The main site of focus was Mount Meager in the south- east of the Complex where the most recent landslide took place. Two profiles through Mount Meager were analysed along with one other location in the northern section of the MMVC, where instability had been detected. The lowest Factor of Safety (FOS) for each profile was displayed and an estimate of the volume which could be generated was deduced. A hazard map showing the inundation zones for various volumes of debris flows was created from simulations using LAHARZ. Results showed the massif is unstable, even before infiltration. Varying the amount of infiltration appears to have no significant impact on the FOS annually implying that small changes of any kind could also trigger failure. Further modelling could be done to assess the impact of infiltration over shorter time scales. The Slide models show the volume of material that could be delivered to the Lillooet River Valley to be of the order of 109 m3 which, based on the LAHARZ simulations, would completely inundate the valley and communities downstream. A major hazard of this is that the removal of such a large amount of material has the potential to trigger an explosive eruption of the geothermal system and renew volcanic activity. Although events of this size are infrequent, there is a significant risk to the communities downstream of the complex.

STABILITY ANALYSIS OF THE SLOPE ALONG US-2 BETWEEN EPOUFETTE BAY AND THE CUT RIVER BRIDGE

Recently, water was observed flowing from a section of steep slope along US-2 near St. Ignace, Michigan in addition to soil sloughing in the area where the water is flowing from the slope. An inspection of the area also showed the presence of sinkholes. The original construction drawing for US-2 also indicated that sinkholes were present in this area prior to road construction in 1948. An investigation was conducted to determine the overall stability of the slope. The slope consists primarily of aeolian sand deposits. Laboratory testing determined the shear strength of the slope material to have a friction angle around 30°, which is also the slope angle. Thus, the slope is at its maximum angle for stability—however, the slope is also heavily wooded which provides additional support to the slope. Although the area surrounding the water flow has been sloughing, the remaining slope remains intact.

STABILITY ANALYSIS AND HAZARD ASSESSMENT OF THE NORTHERN SLOPES OF SAN VICENTE VOLCANO IN CENTRAL EL SALVADOR

Geologic hazards affect the lives of millions of people worldwide every year. El Salvador is a country that is regularly affected by natural disasters, including earthquakes, volcanic eruptions and tropical storms. Additionally, rainfall-induced landslides and debris flows are a major threat to the livelihood of thousands. The San Vicente Volcano in central El Salvador has a recurring and destructive pattern of landslides and debris flows occurring on the northern slopes of the volcano. In recent memory there have been at least seven major destructive debris flows on San Vicente volcano. Despite this problem, there has been no known attempt to study the inherent stability of these volcanic slopes and to determine the thresholds of rainfall that might lead to slope instability. This thesis explores this issue and outlines a suggested method for predicting the likelihood of slope instability during intense rainfall events. The material properties obtained from a field campaign and laboratory testing were used for a 2-D slope stability analysis on a recent landslide on San Vicente volcano. This analysis confirmed that the surface materials of the volcano are highly permeable and have very low shear strength and provided insight into the groundwater table behavior during a rainstorm. The biggest factors on the stability of the slopes were found to be slope geometry, rainfall totals and initial groundwater table location. Using the results from this analysis a stability chart was created that took into account these main factors and provided an estimate of the stability of a slope in various rainfall scenarios. This chart could be used by local authorities in the event of a known extreme rainfall event to help make decisions regarding possible evacuation. Recommendations are given to improve the methodology for future application in other areas as well as in central El Salvador.

Hidden vorticity in binary Bose-Einstein condensates

We consider a binary Bose-Einstein condensate (BEC) described by a system of two-dimensional (2D) Gross-Pitaevskii equations with the harmonic-oscillator trapping potential. The intraspecies interactions are attractive, while the interaction between the species may have either sign. The same model applies to the copropagation of bimodal beams in photonic-crystal fibers. We consider a family of trapped hidden-vorticity (HV) modes in the form of bound states of two components with opposite vorticities S(1,2) = +/- 1, the total angular momentum being zero. A challenging problem is the stability of the HV modes. By means of a linear-stability analysis and direct simulations, stability domains are identified in a relevant parameter plane. In direct simulations, stable HV modes feature robustness against large perturbations, while unstable ones split into fragments whose number is identical to the azimuthal index of the fastest growing perturbation eigenmode. Conditions allowing for the creation of the HV modes in the experiment are discussed too. For comparison, a similar but simpler problem is studied in an analytical form, viz., the modulational instability of an HV state in a one-dimensional (1D) system with periodic boundary conditions (this system models a counterflow in a binary BEC mixture loaded into a toroidal trap or a bimodal optical beam coupled into a cylindrical shell). We demonstrate that the stabilization of the 1D HV modes is impossible, which stresses the significance of the stabilization of the HV modes in the 2D setting.

A review of drainage and spontaneous rupture in free standing thin films with tangentially immobile interfaces

A review of spontaneous rupture in thin films with tangentially immobile interfaces is presented that emphasizes the theoretical developments of film drainage and corrugation growth through the linearization of lubrication theory in a cylindrical geometry. Spontaneous rupture occurs when corrugations from adjacent interfaces become unstable and grow to a critical thickness. A corrugated interface is composed of a number of waveforms and each waveform becomes unstable at a unique transition thickness. The onset of instability occurs at the maximum transition thickness, and it is shown that only upper and lower bounds of this thickness can be predicted from linear stability analysis. The upper bound is equivalent to the Freakel criterion and is obtained from the zeroth order approximation of the H-3 term in the evolution equation. This criterion is determined solely by the film radius, interfacial tension and Hamaker constant. The lower bound is obtained from the first order approximation of the H-3 term in the evolution equation and is dependent on the film thinning velocity A semi-empirical equation, referred to as the MTR equation, is obtained by combining the drainage theory of Manev et al. [J. Dispersion Sci. Technol., 18 (1997) 769] and the experimental measurements of Radoev et al. [J. Colloid Interface Sci. 95 (1983) 254] and is shown to provide accurate predictions of film thinning velocity near the critical thickness of rupture. The MTR equation permits the prediction of the lower bound of the maximum transition thickness based entirely on film radius, Plateau border radius, interfacial tension, temperature and Hamaker constant. The MTR equation extrapolates to Reynolds equation under conditions when the Plateau border pressure is small, which provides a lower bound for the maximum transition thickness that is equivalent to the criterion of Gumerman and Homsy [Chem. Eng. Commun. 2 (1975) 27]. The relative accuracy of either bound is thought to be dependent on the amplitude of the hydrodynamic corrugations, and a semiempirical correlation is also obtained that permits the amplitude to be calculated as a function of the upper and lower bound of the maximum transition thickness. The relationship between the evolving theoretical developments is demonstrated by three film thickness master curves, which reduce to simple analytical expressions under limiting conditions when the drainage pressure drop is controlled by either the Plateau border capillary pressure or the van der Waals disjoining pressure. The master curves simplify solution of the various theoretical predictions enormously over the entire range of the linear approximation. Finally, it is shown that when the Frenkel criterion is used to assess film stability, recent studies reach conclusions that are contrary to the relevance of spontaneous rupture as a cell-opening mechanism in foams. (C) 2003 Elsevier Science B.V. All rights reserved.

Fluctuations in domain growth: Ginzburg-Landau equations with multiplicative noise

Ginzburg-Landau equations with multiplicative noise are considered, to study the effects of fluctuations in domain growth. The equations are derived from a coarse-grained methodology and expressions for the resulting concentration-dependent diffusion coefficients are proposed. The multiplicative noise gives contributions to the Cahn-Hilliard linear-stability analysis. In particular, it introduces a delay in the domain-growth dynamics.

Interfacial growth in driven diffusive systems

We have studied the growth of interfaces in driven diffusive systems well below the critical temperature by means of Monte Carlo simulations. We consider the region beyond the linear regime and of large values of the external field which has not been explored before. The simulations support the existence of interfacial traveling waves when asymmetry is introduced in the model, a result previously predicted by a linear-stability analysis. Furthermore, the generalization of the Gibbs-Thomson relation is discussed. The results provide evidence that the external field is a stabilizing effect which can be considered as effectively increasing the surface tension.

Effects of external noise on the Swift-Hohenberg Equation

The Swift-Hohenberg equation is studied in the presence of a multiplicative noise. This stochastic equation could describe a situation in which a noise has been superimposed on the temperature gradient between the two plates of a Rayleigh-Bnard cell. A linear stability analysis and numerical simulations show that, in constrast to the additive-noise case, convective structures appear in a regime in which a deterministic analysis predicts a homogeneous solution.

Viscous fingering as a paradigm of interfacial pattern formation: Recent results and new challenges

We review recent results on dynamical aspects of viscous fingering. The Saffman¿Taylor instability is studied beyond linear stability analysis by means of a weakly nonlinear analysis and the exact determination of the subcritical branch. A series of contributions pursuing the idea of a dynamical solvability scenario associated to surface tension in analogy with the traditional selection theory is put in perspective and discussed in the light of the asymptotic theory of Tanveer and co-workers. The inherently dynamical singular effects of surface tension are clarified. The dynamical role of viscosity contrast is explored numerically. We find that the basin of attraction of the Saffman¿Taylor finger depends on viscosity contrast, and that the sensitivity to this parameter is maximal in the usual limit of high viscosity contrast. The competing attractors are identified as closed bubble solutions. We briefly report on recent results and work in progress concerning rotating Hele-Shaw flows, topological singularities and wetting effects, and also discuss future directions in the context of viscous fingering

Interfacial instabilities of a fluid annulus in a rotating Hele-Shaw cell

We have studied the interfacial instabilities experienced by a liquid annulus as it moves radially in a circular Hele-Shaw cell rotating with angular velocity Omega. The instability of the leading interface (oil displacing air) is driven by the density difference in the presence of centrifugal forcing, while the instability of the trailing interface (air displacing oil) is driven by the large viscosity contrast. A linear stability analysis shows that the stability of the two interfaces is coupled through the pressure field already at a linear level. We have performed experiments in a dry cell and in a cell coated with a thin fluid layer on each plate, and found that the stability depends substantially on the wetting conditions at the leading interface. Our experimental results of the number of fingers resulting from the instability compare well with the predictions obtained through a numerical integration of the coupled equations derived from a linear stability analysis. Deep in the nonlinear regime we observe the emission of liquid droplets through the formation of thin filaments at the tip of outgrowing fingers.