A curve shortening flow rule for closed embedded plane curves with a prescribed rate of change in enclosed area

Motivated by a problem from fluid mechanics, we consider a generalization of the standard curve shortening flow problem for a closed embedded plane curve such that the area enclosed by the curve is forced to decrease at a prescribed rate. Using formal asymptotic and numerical techniques, we derive possible extinction shapes as the curve contracts to a point, dependent on the rate of decreasing area; we find there is a wider class of extinction shapes than for standard curve shortening, for which initially simple closed curves are always asymptotically circular. We also provide numerical evidence that self-intersection is possible for non-convex initial conditions, distinguishing between pinch-off and coalescence of the curve interior.

Evaluating Methods to Predict Streamflow at Ungauged Sites using Regional Flow Duration Curves: A Case Study

Precise information on streamflows is of major importance for planning and monitoring of water resources schemes related to hydro power, water supply, irrigation, flood control, and for maintaining ecosystem. Engineers encounter challenges when streamflow data are either unavailable or inadequate at target locations. To address these challenges, there have been efforts to develop methodologies that facilitate prediction of streamflow at ungauged sites. Conventionally, time intensive and data exhaustive rainfall-runoff models are used to arrive at streamflow at ungauged sites. Most recent studies show improved methods based on regionalization using Flow Duration Curves (FDCs). A FDC is a graphical representation of streamflow variability, which is a plot between streamflow values and their corresponding exceedance probabilities that are determined using a plotting position formula. It provides information on the percentage of time any specified magnitude of streamflow is equaled or exceeded. The present study assesses the effectiveness of two methods to predict streamflow at ungauged sites by application to catchments in Mahanadi river basin, India. The methods considered are (i) Regional flow duration curve method, and (ii) Area Ratio method. The first method involves (a) the development of regression relationships between percentile flows and attributes of catchments in the study area, (b) use of the relationships to construct regional FDC for the ungauged site, and (c) use of a spatial interpolation technique to decode information in FDC to construct streamflow time series for the ungauged site. Area ratio method is conventionally used to transfer streamflow related information from gauged sites to ungauged sites. Attributes that have been considered for the analysis include variables representing hydrology, climatology, topography, land-use/land- cover and soil properties corresponding to catchments in the study area. Effectiveness of the presented methods is assessed using jack knife cross-validation. Conclusions based on the study are presented and discussed. (C) 2015 The Authors. Published by Elsevier B.V.

Similarity solutions for flow and heat transfer of non-Newtonian fluid over a stretching surface

Similarity solutions are carried out for flow of power law non-Newtonian fluid film on unsteady stretching surface subjected to constant heat flux. Free convection heat transfer induces thermal boundary layer within a semi-infinite layer of Boussinesq fluid. The nonlinear coupled partial differential equations (PDE) governing the flow and the boundary conditions are converted to a system of ordinary differential equations (ODE) using two-parameter groups. This technique reduces the number of independent variables by two, and finally the obtained ordinary differential equations are solved numerically for the temperature and velocity using the shooting method. The thermal and velocity boundary layers are studied by the means of Prandtl number and non-Newtonian power index plotted in curves.

On Random Planar Curves and Their Scaling Limits

Planar curves arise naturally as interfaces between two regions of the plane. An important part of statistical physics is the study of lattice models. This thesis is about the interfaces of 2D lattice models. The scaling limit is an infinite system limit which is taken by letting the lattice mesh decrease to zero. At criticality, the scaling limit of an interface is one of the SLE curves (Schramm-Loewner evolution), introduced by Oded Schramm. This family of random curves is parametrized by a real variable, which determines the universality class of the model. The first and the second paper of this thesis study properties of SLEs. They contain two different methods to study the whole SLE curve, which is, in fact, the most interesting object from the statistical physics point of view. These methods are applied to study two symmetries of SLE: reversibility and duality. The first paper uses an algebraic method and a representation of the Virasoro algebra to find common martingales to different processes, and that way, to confirm the symmetries for polynomial expected values of natural SLE data. In the second paper, a recursion is obtained for the same kind of expected values. The recursion is based on stationarity of the law of the whole SLE curve under a SLE induced flow. The third paper deals with one of the most central questions of the field and provides a framework of estimates for describing 2D scaling limits by SLE curves. In particular, it is shown that a weak estimate on the probability of an annulus crossing implies that a random curve arising from a statistical physics model will have scaling limits and those will be well-described by Loewner evolutions with random driving forces.

Continuous-Flow Surface Aeration Systems

An aeration process in ail activated sludge plant is a continuous-flow system. In this system, there is a steady input flow (flow from the primary clarifier or settling tank with some part from the secondary clarifier or secondary settling tank) and output flow connection to the secondary clarifier or settling tank. The experimental and numerical results obtained through batch systems can not be relied on and applied for the designing of a continuous aeration tank. In order to scale up laboratory results for field application, it is imperative to know the geometric parameters of a continuous system. Geometric parameters have a greater influence on the mass transfer process of surface aeration systems. The present work establishes the optimal geometric configuration of a continuous-flow surface aeration system. It is found that the maintenance of these optimal geometric parameters systems result in maximum aeration efficiency. By maintaining the obtained optimal geometric parameters, further experiments are conducted in continuous-flow surface aerators with three different sizes in order to develop design curves correlating the oxygen transfer coefficient and power number with the rotor speed. The design methodology to implement the presently developed optimal geometric parameters and correlation equations for field application is discussed.

Aspiration-based utility functions in a planning model for timber flow management.

The study presents a theory of utility models based on aspiration levels, as well as the application of this theory to the planning of timber flow economics. The first part of the study comprises a derivation of the utility-theoretic basis for the application of aspiration levels. Two basic models are dealt with: the additive and the multiplicative. Applied here solely for partial utility functions, aspiration and reservation levels are interpreted as defining piecewisely linear functions. The standpoint of the choices of the decision-maker is emphasized by the use of indifference curves. The second part of the study introduces a model for the management of timber flows. The model is based on the assumption that the decision-maker is willing to specify a shape of income flow which is different from that of the capital-theoretic optimum. The utility model comprises four aspiration-based compound utility functions. The theory and the flow model are tested numerically by computations covering three forest holdings. The results show that the additive model is sensitive even to slight changes in relative importances and aspiration levels. This applies particularly to nearly linear production possibility boundaries of monetary variables. The multiplicative model, on the other hand, is stable because it generates strictly convex indifference curves. Due to a higher marginal rate of substitution, the multiplicative model implies a stronger dependence on forest management than the additive function. For income trajectory optimization, a method utilizing an income trajectory index is more efficient than one based on the use of aspiration levels per management period. Smooth trajectories can be attained by squaring the deviations of the feasible trajectories from the desired one.

Influence of Osmotic Suction on the Soil-Water Characteristic Curves of Compacted Expansive Clay

Unsaturated clays are subject to osmotic suction gradients in geoenvironmental engineering applications and it therefore becomes important to understand the effect of these chemical concentration gradients on soil-water characteristic curves (SWCCs). This paper brings out the influence of induced osmotic suction gradient on the wetting SWCCs of compacted clay specimens inundated with sodium chloride solutions/distilled water at vertical stress of 6.25 kPa in oedometer cells. The experimental results illustrate that variations in initial osmotic suction difference induce different magnitudes of osmotic induced consolidation and osmotic consolidation strains thereby impacting the wetting SWCCs and equilibrium water contents of identically compacted clay specimens. Osmotic suction induced by chemical concentration gradients between reservoir salt solution and soil-water can be treated as an equivalent net stress component, (p(pi)) that decreases the swelling strains of unsaturated specimens from reduction in microstructural and macrostructural swelling components. The direction of osmotic flow affects the matric SWCCs. Unsaturated specimens experiencing osmotic induced consolidation and osmotic consolidation develop lower equilibrium water content than specimens experiencing osmotic swelling during the wetting path. The findings of the study illustrate the need to incorporate the influence of osmotic suction in determination of the matric SWCCs.

Flow Instabilities and fracture in Ti-6Al-4V deformed in compression at 298 K to 673 K

Uniaxial compression tests were conducted on Ti-6Al-4V specimens in the strain-rate range df 0.001 to 1 s(-1) and temperature range of 298 to 673 K. The stress-strain curves exhibited a peak flow stress followed by flow softening. Up to 523 K, the specimens cracked catastrophically after the flow softening started. Adiabatic shear banding was observed in this regime. The fracture surface exhibited both mode I and II fracture features. The state of stress existing in a compression test specimen when bulging occurs is responsible for this fracture. The instabilities observed in the present tests are classified as ''geometric'' in nature and are state-of-stress dependant, unlike the ''intrinsic'' instabilities, which are dependant on the dynamic constitutive behavior of the material.

Stability of the viscous flow of a fluid through a flexible tube

The stability of Hagen-Poiseuille flow of a Newtonian fluid of viscosity eta in a tube of radius R surrounded by a viscoelastic medium of elasticity G and viscosity eta(s) occupying the annulus R < r < HR is determined using a linear stability analysis. The inertia of the fluid and the medium are neglected, and the mass and momentum conservation equations for the fluid and wall are linear. The only coupling between the mean flow and fluctuations enters via an additional term in the boundary condition for the tangential velocity at the interface, due to the discontinuity in the strain rate in the mean flow at the surface. This additional term is responsible for destabilizing the surface when the mean velocity increases beyond a transition value, and the physical mechanism driving the instability is the transfer of energy from the mean flow to the fluctuations due to the work done by the mean flow at the interface. The transition velocity Gamma(t) for the presence of surface instabilities depends on the wavenumber k and three dimensionless parameters: the ratio of the solid and fluid viscosities eta(r) = (eta(s)/eta), the capillary number Lambda = (T/GR) and the ratio of radii H, where T is the surface tension of the interface. For eta(r) = 0 and Lambda = 0, the transition velocity Gamma(t) diverges in the limits k much less than 1 and k much greater than 1, and has a minimum for finite k. The qualitative behaviour of the transition velocity is the same for Lambda > 0 and eta(r) = 0, though there is an increase in Gamma(t) in the limit k much greater than 1. When the viscosity of the surface is non-zero (eta(r) > 0), however, there is a qualitative change in the Gamma(t) vs. k curves. For eta(r) < 1, the transition velocity Gamma(t) is finite only when k is greater than a minimum value k(min), while perturbations with wavenumber k < k(min) are stable even for Gamma--> infinity. For eta(r) > 1, Gamma(t) is finite only for k(min) < k < k(max), while perturbations with wavenumber k < k(min) or k > k(max) are stable in the limit Gamma--> infinity. As H decreases or eta(r) increases, the difference k(max)- k(min) decreases. At minimum value H = H-min, which is a function of eta(r), the difference k(max)-k(min) = 0, and for H < H-min, perturbations of all wavenumbers are stable even in the limit Gamma--> infinity. The calculations indicate that H-min shows a strong divergence proportional to exp (0.0832 eta(r)(2)) for eta(r) much greater than 1.

Hidden order in serrated flow of metallic glasses

We report results of statistical and dynamic analysis of the serrated stress-time curves obtained from compressive constant strain-rate tests on two metallic glass samples with different ductility levels in an effort to extract hidden information in the seemingly irregular serrations. Two distinct types of dynamics are detected in these two alloy samples. The stress-strain curve corresponding to the less ductile Zr65Cu15Ni10Al10 alloy is shown to exhibit a finite correlation dimension and a positive Lyapunov exponent, suggesting that the underlying dynamics is chaotic. In contrast, for the more ductile Cu47.5Zr47.5Al5 alloy, the distributions of stress drop magnitudes and their time durations obey a power-law scaling reminiscent of a self-organized critical state. The exponents also satisfy the scaling relation compatible with self-organized criticality. Possible physical mechanisms contributing to the two distinct dynamic regimes are discussed by drawing on the analogy with the serrated yielding of crystalline samples. The analysis, together with some physical reasoning, suggests that plasticity in the less ductile sample can be attributed to stick-slip of a single shear band, while that of the more ductile sample could be attributed to the simultaneous nucleation of a large number of shear bands and their mutual interactions. (C) 2011 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Stability of non-parabolic flow in a flexible tube

Flows with velocity profiles very different from the parabolic velocity profile can occur in the entrance region of a tube as well as in tubes with converging/diverging cross-sections. In this paper, asymptotic and numerical studies are undertaken to analyse the temporal stability of such 'non-parabolic' flows in a flexible tube in the limit of high Reynolds numbers. Two specific cases are considered: (i) developing flow in a flexible tube; (ii) flow in a slightly converging flexible tube. Though the mean velocity profile contains both axial and radial components, the flow is assumed to be locally parallel in the stability analysis. The fluid is Newtonian and incompressible, while the flexible wall is modelled as a viscoelastic solid. A high Reynolds number asymptotic analysis shows that the non-parabolic velocity profiles can become unstable in the inviscid limit. This inviscid instability is qualitatively different from that observed in previous studies on the stability of parabolic flow in a flexible tube, and from the instability of developing flow in a rigid tube. The results of the asymptotic analysis are extended numerically to the moderate Reynolds number regime. The numerical results reveal that the developing flow could be unstable at much lower Reynolds numbers than the parabolic flow, and hence this instability can be important in destabilizing the fluid flow through flexible tubes at moderate and high Reynolds number. For flow in a slightly converging tube, even small deviations from the parabolic profile are found to be sufficient for the present instability mechanism to be operative. The dominant non-parallel effects are incorporated using an asymptotic analysis, and this indicates that non-parallel effects do not significantly affect the neutral stability curves. The viscosity of the wall medium is found to have a stabilizing effect on this instability.

A dynamical instability due to fluid–wall coupling lowers the transition Reynolds number in the flow through a flexible tube

A flow-induced instability in a tube with flexible walls is studied experimentally. Tubes of diameter 0.8 and 1.2 mm are cast in polydimethylsiloxane (PDMS) polymer gels, and the catalyst concentration in these gels is varied to obtain shear modulus in the range 17–550 kPa. A pressure drop between the inlet and outlet of the tube is used to drive fluid flow, and the friction factor $f$ is measured as a function of the Reynolds number $Re$. From these measurements, it is found that the laminar flow becomes unstable, and there is a transition to a more complicated flow profile, for Reynolds numbers as low as 500 for the softest gels used here. The nature of the $f$–$Re$ curves is also qualitatively different from that in the flow past rigid tubes; in contrast to the discontinuous increase in the friction factor at transition in a rigid tube, it is found that there is a continuous increase in the friction factor from the laminar value of $16\ensuremath{/} Re$ in a flexible tube. The onset of transition is also detected by a dye-stream method, where a stream of dye is injected into the centre of the tube. It is found that there is a continuous increase of the amplitude of perturbations at the onset of transition in a flexible tube, in contrast to the abrupt disruption of the dye stream at transition in a rigid tube. There are oscillations in the wall of the tube at the onset of transition, which is detected from the laser scattering off the walls of the tube. This indicates that the coupling between the fluid stresses and the elastic stresses in the wall results in an instability of the laminar flow.