Modified alternating direction methods for solving a two-dimensional non-continuous seepage flow with fractional derivatives

In this paper, a two-dimensional non-continuous seepage flow with fractional derivatives (2D-NCSF-FD) in uniform media is considered, which has modified the well known Darcy law. Using the relationship between Riemann-Liouville and Grunwald-Letnikov fractional derivatives, two modified alternating direction methods: a modified alternating direction implicit Euler method and a modified Peaceman-Rachford method, are proposed for solving the 2D-NCSF-FD in uniform media. The stability and consistency, thus convergence of the two methods in a bounded domain are discussed. Finally, numerical results are given.

Numerical simulation for the 3D seepage flow with fractional derivatives in porous media

In this paper, the numerical simulation of the 3D seepage flow with fractional derivatives in porous media is considered under two special cases: non-continued seepage flow in uniform media (NCSFUM) and continued seepage flow in non-uniform media (CSF-NUM). A fractional alternating direction implicit scheme (FADIS) for the NCSF-UM and a modified Douglas scheme (MDS) for the CSF-NUM are proposed. The stability, consistency and convergence of both FADIS and MDS in a bounded domain are discussed. A method for improving the speed of convergence by Richardson extrapolation for the MDS is also presented. Finally, numerical results are presented to support our theoretical analysis.

Modeling of seepage flow through layered solis

A two-dimensional finite difference model, which solves mixed type of Richards' equation, whose non-linearity is dealt with modified Picard's iteration and strongly implicit procedure to solve the resulting equations, is presented. Modeling of seepage flow through heterogeneous soils, which is common in the field is addressed in the present study. The present model can be applied to both unsaturated and saturated soils and can handle very dry initial condition and steep wetting fronts. The model is validated by comparing experimental results reported in the literature. Newness of this two dimensional model is its application on layered soils with transient seepage face development, which has not been reported in the literature. Application of the two dimensional model for studying unconfined drainage due to sudden drop of water table at seepage face in layered soils is demonstrated. In the present work different sizes of rectangular flow domain with different types of layering are chosen. Sensitivity of seepage height due to problem dimension of layered system is studied. The effect of aspect ratio on seepage face development in case of the flow through layered soil media is demonstrated. The model is also applied to random heterogeneous soils in which the randomness of the model parameters is generated using the turning band technique. The results are discussed in terms of phreatic surface and seepage height development and also flux across the seepage face. Such accurate modeling of seepage face development and quantification of flux moving across the seepage face becomes important while modeling transport problems in variably saturated media.

Fundamentals of steady state, non-Newtonian rimming flow

Rimming flow on the inner surface of a horizontal rotating cylinder is investigated. Using a scale analysis, a theoretical description is obtained for steady-state non-Newtonian flow. Simple lubrication theory is applied since the Reynolds number is small and the liquid film is thin. Since the Deborah number is very small the flow is viscometric. The shear-thinning number, which characterizes the shear-thinning effect, may be small or large. A general constitutive law for this kind of flow requires only a single function relating shear stress and shear rate that corresponds to a generalized Newtonian liquid. For this case the run-off condition for rimming flow is derived. Provided the run-off condition is satisfied, the existence of a continuous steady-state solution is proved. The rheological models, which show Newtonian behavior at low shear rates with transition to power-law shear thinning at moderate shear rates, are considered. Numerical results are carried out for the Carreau and Ellis models, which exhibit Newtonian behavior near the free surface and power-law behavior near the wall of the rotating cylinder.

Stability of the tunnel face reinforced by bolts under seepage flow conditions

Il presente lavoro tratta la stabilità del fronte di scavo, rinforzato con barre di consolidamento ed interessato da drenaggi in avanzamento, di gallerie sotto falda in rocce tenere o terreni. Tale studio è stato sviluppato dal progetto di Tesi attraverso l’analisi all’equilibrio limite che approssima il fronte di scavo con un rettangolo e considera un meccanismo di rottura composto da un cuneo, a tergo del fronte, caricato da un prisma. Il metodo descritto consente di tenere conto dell’effetto stabilizzante delle barre, mediante una distribuzione della pressione di supporto non uniforme. Nel caso di gallerie sotto falda, lo stesso metodo permette inoltre di considerare l’effetto destabilizzante dei gradienti idraulici. Sono state ricavate soluzioni analitiche per la valutazione della stabilità, ed implementate successivamente nel software di analisi numerica MATLAB. Dalle analisi condotte è emerso che il numero minimo di barre per garantire la stabilità del fronte di scavo è in molti casi elevato e risulta impossibile da porre in opera in terreni scarsamente coesivi o in gallerie sotto elevati battenti d’acqua. Per risolvere questa situazione si può prevedere l’inserimento di drenaggi in avanzamento, con lo scopo di diminuire i gradienti idraulici nei pressi del fronte della galleria. Il modello che descrive il nuovo andamento dei carichi idraulici, considerando la presenza di dreni, è stato realizzato con il software commerciale agli elementi finiti COMSOL. Una volta determinati gli andamenti dei carichi idraulici, sono stati condotti studi parametrici sull’effetto dei dreni combinato con gli elementi di rinforzo. Dopo tali analisi sono stati ricavati nomogrammi adimensionali che tengano conto della presenza contemporanea delle barre e dei dreni. Tali diagrammi costituiscono uno strumento utile e valido per la progettazione del rinforzo del fronte di scavo. Infine sono stati realizzati confronti fra casi di studio reali e risultati ottenuti dal modello.

Pre-crash and non-crash traffic flow trends analysis on motorways

Crashes that occur on motorways contribute to a significant proportion (40-50%) of non-recurrent motorway congestions. Hence, reducing the frequency of crashes assists in addressing congestion issues (Meyer, 2008). Crash likelihood estimation studies commonly focus on traffic conditions in a short time window around the time of a crash while longer-term pre-crash traffic flow trends are neglected. In this paper we will show, through data mining techniques that a relationship between pre-crash traffic flow patterns and crash occurrence on motorways exists. We will compare them with normal traffic trends and show this knowledge has the potential to improve the accuracy of existing models and opens the path for new development approaches. The data for the analysis was extracted from records collected between 2007 and 2009 on the Shibuya and Shinjuku lines of the Tokyo Metropolitan Expressway in Japan. The dataset includes a total of 824 rear-end and sideswipe crashes that have been matched with crashes corresponding to traffic flow data using an incident detection algorithm. Traffic trends (traffic speed time series) revealed that crashes can be clustered with regards to the dominant traffic patterns prior to the crash. Using the K-Means clustering method with Euclidean distance function allowed the crashes to be clustered. Then, normal situation data was extracted based on the time distribution of crashes and were clustered to compare with the “high risk” clusters. Five major trends have been found in the clustering results for both high risk and normal conditions. The study discovered traffic regimes had differences in the speed trends. Based on these findings, crash likelihood estimation models can be fine-tuned based on the monitored traffic conditions with a sliding window of 30 minutes to increase accuracy of the results and minimize false alarms.

Influence of spatial variability of permeability property on steady state seepage flow and slope stability analysis

In recent years, spatial variability modeling of soil parameters using random field theory has gained distinct importance in geotechnical analysis. In the present Study, commercially available finite difference numerical code FLAC 5.0 is used for modeling the permeability parameter as spatially correlated log-normally distributed random variable and its influence on the steady state seepage flow and on the slope stability analysis are studied. Considering the case of a 5.0 m high cohesive-frictional soil slope of 30 degrees, a range of coefficients of variation (CoV%) from 60 to 90% in the permeability Values, and taking different values of correlation distance in the range of 0.5-15 m, parametric studies, using Monte Carlo simulations, are performed to study the following three aspects, i.e., (i) effect ostochastic soil permeability on the statistics of seepage flow in comparison to the analytic (Dupuit's) solution available for the uniformly constant permeability property; (ii) strain and deformation pattern, and (iii) stability of the given slope assessed in terms of factor of safety (FS). The results obtained in this study are useful to understand the role of permeability variations in slope stability analysis under different slope conditions and material properties. (C) 2009 Elsevier B.V. All rights reserved.

Non-linear couette flow with heat transfer using the B-G-K model

In recent years a large number of investigators have devoted their efforts to the study of flow and heat transfer in rarefied gases, using the BGK [1] model or the Boltzmann kinetic equation. The velocity moment method which is based on an expansion of the distribution function as a series of orthogonal polynomials in velocity space, has been applied to the linearized problem of shear flow and heat transfer by Mott-Smith [2] and Wang Chang and Uhlenbeck [3]. Gross, Jackson and Ziering [4] have improved greatly upon this technique by expressing the distribution function in terms of half-range functions and it is this feature which leads to the rapid convergence of the method. The full-range moments method [4] has been modified by Bhatnagar [5] and then applied to plane Couette flow using the B-G-K model. Bhatnagar and Srivastava [6] have also studied the heat transfer in plane Couette flow using the linearized B-G-K equation. On the other hand, the half-range moments method has been applied by Gross and Ziering [7] to heat transfer between parallel plates using Boltzmann equation for hard sphere molecules and by Ziering [83 to shear and heat flow using Maxwell molecular model. Along different lines, a moment method has been applied by Lees and Liu [9] to heat transfer in Couette flow using Maxwell's transfer equation rather than the Boltzmann equation for distribution function. An iteration method has been developed by Willis [10] to apply it to non-linear heat transfer problems using the B-G-K model, with the zeroth iteration being taken as the solution of the collisionless kinetic equation. Krook [11] has also used the moment method to formulate the equivalent continuum equations and has pointed out that if the effects of molecular collisions are described by the B-G-K model, exact numerical solutions of many rarefied gas-dynamic problems can be obtained. Recently, these numerical solutions have been obtained by Anderson [12] for the non-linear heat transfer in Couette flow,

NON-COERCIVE RICCI FLOW INVARIANT CURVATURE CONES

This note is a study of nonnegativity conditions on curvature preserved by the Ricci flow. We focus on a specific class of curvature conditions which we call non-coercive: These are the conditions for which nonnegative curvature and vanishing scalar curvature does not imply flatness. We show, in dimensions greater than 4, that if a Ricci flow invariant nonnegativity condition is satisfied by all Einstein curvature operators with nonnegative scalar curvature, then this condition is just the nonnegativity of scalar curvature. As a corollary, we obtain that a Ricci flow invariant curvature condition, which is stronger than a nonnegative scalar curvature, cannot be strictly satisfied by curvature operators (other than multiples of the identity) of compact Einstein symmetric spaces. We also investigate conditions which are satisfied by all conformally flat manifolds with nonnegative scalar curvature.

Kinetic model of continuous-wave flow chemical lasers

In this paper an analysis of the kinetic theory of the continuous-wave flow chemical lasers(CWFCL) is presented with emphasis being laid on the effects of inhomogeneous broadeningon CWFCL's performance. The results obtained are applicable to the case where laser fre-quency is either coincident or incoincident with that of the eenter of the line shape. This rela-tion has been,compared with that of the rate model in common use. These two models are almostidentical as the broadening parameter η is larger than 1. The smaller the value of η, thegreater the difference between the results of these two models will be. For fixed η, the dif-ferences between fhe results of the two models increase with the increase of the frequencyshift parameter ξ. When η is about less than 0.2. the kinetic model can predict exactly the in-homogeneous broadening effects,while the rate model cannot.