Particle systems on flows

We start with a stochastic flow of diffeomorphisms of the space. Particles enter the space at random times and places. Each particle is carried by the flow for some random amount of time. We examine the point process formed by the particles at a fixed time, on the evolution of that point process as time varies, and on the equilibrium law of the point process.     

Application of Smoothed Particle Hydrodynamics for modelling gated spillway flows

Computational models of spillways are important for evaluating and improving dam safety, optimising spillway design and updating operating conditions. Traditionally, scaled down physical models have been used for validation and to collect hydraulic data. Computational fluid dynamics (CFD) models however provide advantages in time, cost and resource reduction. CFD models also provide greater efficiency when evaluating a range of spillway designs or operating conditions. Within the present literature, most studies of computational spillway models utilise a mesh-based method. In this work we use the particle based method of Smoothed Particle Hydrodynamics (SPH) to model weir flow through a four bay, gated, spillway system. Advantages of SPH for such modelling include automatic representation of the free surface flow behaviour due to the Lagrangian nature of the method, and the ability to incorporate complex and dynamic boundary objects such as gate structures or debris. To validate the SPH model, the reservoir water depth simulated is compared with a related physical study. The effect of SPH resolution on the predicted water depth is evaluated. The change in reservoir water level with discharge rates for weir flow conditions is also investigated, with the difference in simulated and experimental water depths found to range from 0.16% to 11.48%. These results are the first quantitative validation of the SPH method to capture spillway flow in three dimensions. The agreement achieved demonstrates the capability of the SPH method for modelling spillway flows.     

Particle trajectories beneath small amplitude shallow water waves in constant vorticity flows

We investigate the particle trajectories in a constant vorticity shallow water flow over a flat bed as periodic waves propagate on the water’s free surface. Within the framework of small amplitude waves, we find the solutions of the nonlinear differential equations system which describes the particle motion in the considered case, and we describe the possible particle trajectories. Depending on the relation between the initial data and the constant vorticity, some particle trajectories are undulating curves to the right, or to the left, others are loops with forward drift, or with backward drift, others can follow some peculiar shapes.     

Multi-objective analysis of dam release flows in rivers downstream from hydropower reservoirs

This paper assesses dam releases from hydropower reservoirs in order to optimize power production and fish habitat protection. A multi-objective programming model includes output from 2-D hydraulic simulation for habitat assessment to optimize power production and fish habitat suitability as a Pareto set. To identify the optimal Pareto set two different approaches are used and compared: ε-constraint methods and non-dominant-sorting genetic algorithm (NSGA II). To formulate the ecological objective the river habitat quality is quantified by the weighted usable area (WUA). The relation between the WUA and the river flow-rate is obtained by using a 2D hydraulic model in which the hydraulic characteristics of river current – depth and velocity – are calculated by a finite difference numerical integration of two-dimensional shallow water equations on a boundary fitted non orthogonal curvilinear mesh. This approach allows the integration of motion equations on geometrically complex domains as those representing the morphology of natural watercourses. The performance of the proposed methodology is analyzed in a case study of a stretch of the Piave river downstream of the dam of the Pieve di Cadore reservoir (Belluno, Italy).     

Nonlinear mechanics of accretion

We formulate a geometric nonlinear theory of the mechanics of accretion. In this theory, the material manifold of an accreting body is represented by a time-dependent Riemannian manifold with a time-independent metric that at each point depends on the state of deformation at that point at its time of attachment to the body, and on the way the new material is added to the body. We study the incompatibilities induced by accretion through the analysis of the material metric and its curvature in relation to the foliated structure of the accreted body. Balance laws are discussed and the initial boundary value problem of accretion is formulated. The particular cases where the growth surface is either fixed or traction-free are studied and some analytical results are provided. We numerically solve several accretion problems and calculate the residual stresses in nonlinear elastic bodies induced from accretion.

     

A study of truly incompressible and weakly compressible Smoothed Particle Hydrodynamics methods to model incompressible flows with free surfaces

The Smoothed Particle Hydrodynamics (SPH) method is a meshless discretization method for solving, e.g., the Navier-Stokes equations. By now, it is used for hydraulic problems as well as for solid bodies. In general, there are two distinguishable approaches for incompressible fluid flows. One is called weakly compressible SPH (WCSPH) and the other is called truly incompressible SPH (ISPH). The main difference between these two approaches is the way of pressure evaluation. In WCSPH, a state equation is used, while in ISPH the pressure Poisson equation is solved. Each approach has its advantages as well as its disadvantages, for example the complexity of the numerical algorithm for WCSPH is smaller than for ISPH, but the pressure field is more accurate for ISPH. In this work, two representative examples are studied. The simulations were performed with two different codes, one developed at the Institute of Engineering and Computational Mechanics and one at the Institute of Chemical Process Engineering. It is the aim of this work to show some properties of WCSPH and ISPH as well as to compare two different implementations that, in detail, are very complex. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Structures of magnetized thin accretion disks

We investigate the magnetohydrodynamic (MHD) process in thin accretion disks. The relevant momentum as well as magnetic reduction equations in the thin disk approximation are included. On the basis of these equations, we examine numerically the stationary structures, including distributions of the surface mass density, temperature and flow velocities of a disk around a young stellar object (YSO). The numerical results are as follows: (i) There should be an upper limit to the magnitude of magnetic field, such an upper limit corresponds to the equipartition field. For relevant magnitude of magnetic field of the disk’s interior the disk remains approximately Keplerian. (ii) The distribution of effective temperature T(r) is a smoothly decreasing function of radius with power index γ -1/2, corresponding to the observed radiation flux density, provided that the magnetic field is suitably chosen.     

Multicommodity flows in graphs

Suppose that G is a graph, and (s_i,t_i) (1≤i≤k) are pairs of vertices; and that each edge has a integer-valued capacity (≥0), and that q_i≥0 (1≤i≤k) are integer-valued demands. When is there a flow for each i, between s_i and t_i and of value q_i, such that the total flow through each edge does not exceed its capacity? Ford and Fulkerson solved this when k=1, and Hu when k=2. We solve it for general values of k, when G is planar and can be drawn so that s₁,…, s_l, t₁, …, t_l,…,t_l are all on the boundary of a face and s_l+1, …,S_k, t_l+1,…,t_k are all on the boundary of the infinite face or when t₁=?=t_l and G is planar and can be drawn so that s_l+1,…,s_k, t₁,…,t_k are all on the boundary of the infinite face. This extends a theorem of Okamura and Seymour.     

Nilpotent and recursive flows

In this paper we introduce a class of nonlinear vector fields on infinite dimensional manifolds such that the corresponding evolution equations can be solved with the same method one uses to solve ordinary differential equations with constant coeficients. Mostly, these equations are nonlinear partial differential equations. It is shown that these flows are characterized by a generalization of the ‘method of variation of constants’ which is widely used for second order problems to find general solutions out of particular ones. Invariant densities are constructed for these flows in a natural way. These invariant densities are providing an essential tool for solving initial value and boundary value problems for the equations under consideration. Many applications are presented     

One-point suspensions and wreath products of polytopes and spheres

It is known that the suspension of a simplicial complex can be realized with only one additional point. Suitable iterations of this construction generate highly symmetric simplicial complexes with various interesting combinatorial and topological properties. In particular, infinitely many non-PL spheres as well as contractible simplicial complexes with a vertex-transitive group of automorphisms can be obtained in this way.     

Analyticity and flows in von Neumann algebras

Let $\text{B}$ be a von Neumann algebra, let {α_t}_tεR be an ultraweakly continuous one-parameter group of ¹-automorphisms of $\text{B}$, and let $\text{U}$ be the set of all A such that for each ? in $\text{B}$₁, the function t → ?(α_t(A)) lies in H^∞($\text{R}$. Then $\text{U}$ is an ultraweakly closed subalgebra of $\text{B}$ containing the identity which is proper and non-self-adjoint if {α_t}_tεR is not trivial. In this paper, a systematic investigation into the structure theory of $\text{U}$ is begun. Two of the more note-worthy developments are these. First of all, conditions under which $\text{U}$ is a subdiagonal algebra in $\text{B}$, in the sense of Arveson, are determined. The analysis provides a common perspective from which to view a large number of hitherto unrelated algebras. Second, the invariant subspace structure of $\text{U}$ is determined and conditions under which $\text{U}$ is a reductive subalgebra of $\text{B}$ are found. These results are then used to produce examples where $\text{U}$ is a proper, non-self-adjoint, reductive subalgebra of $\text{B}$. The examples do not answer the reductive algebra question, however, because although ultraweakly closed, the subalgebras are weakly dense in $\text{B}$.