Stability of Viscous St. Venant Roll Waves: From Onset to Infinite Froude Number Limit

We study the spectral stability of roll wave solutions of the viscous St. Venant equations modeling inclined shallow water flow, both at onset in the small Froude number or “weakly unstable” limit \(F\rightarrow 2^+\) and for general values of the Froude number F, including the limit \(F\rightarrow +\infty \). In the former, \(F\rightarrow 2^+\), limit, the shallow water equations are formally approximated by a Korteweg-de Vries/Kuramoto–Sivashinsky (KdV–KS) equation that is a singular perturbation of the standard Korteweg-de Vries (KdV) equation modeling horizontal shallow water flow. Our main analytical result is to rigorously validate this formal limit, showing that stability as \(F\rightarrow 2^+\) is equivalent to stability of the corresponding KdV–KS waves in the KdV limit. Together with recent results obtained for KdV–KS by Johnson–Noble–Rodrigues–Zumbrun and Barker, this gives not only the first rigorous verification of stability for any single viscous St. Venant roll wave, but a complete classification of stability in the weakly unstable limit. In the remainder of the paper, we investigate numerically and analytically the evolution of the stability diagram as Froude number increases to infinity. Notably, we find transition at around \(F=2.3\) from weakly unstable to different, large-F behavior, with stability determined by simple power-law relations. The latter stability criteria are potentially useful in hydraulic engineering applications, for which typically \(2.5\le F\le 6.0\).     

极限分析的广义界限定理

本文研究了非连续流动场中,刚塑性介质极限分析完全解的界限问题.提出了一个包括界面条件及间断面条件在内的混合边值问题的广义变分原理,建立了极限载荷乘子的变分解析公式.并证明了一个新的界限定理,其中的场变量将不再受到屈服条件、不可压缩条件等约束的限制.此定理的推论给出了变分解与完全解之间的关系.初步应用表明,对于简单选取的场变量,由本文公式可以得到准确解的较佳界限值,结果具有较好的稳定性.     

Steady flows of compressible fluids in a rigid container with upper free boundary

We consider fluid in a smooth rigid container whose lateral boundary is a piece of vertical cylinder, bounded above by a free upper surface. As basic flow we consider the non homogeneous rest state in the presence of gravity, and of a surface tension. Under these assumptions, we study the existence of a steady free boundary and a steady motion in of an isothermal viscous gas, resulting as perturbation to the rest state in correspondence of small non potential perturbations to the (large potential) gravitational force. We linearize the problem by prescribing the unknown domain , then we make use of the iterative scheme introduced by Heywood and Padula. Our method is based on an iteration between the Neumann problem for a non homogeneous Stokes system for the velocity, the Neumann problem for an elliptic problem on for height, and a steady transport equation for the perturbation to the density. The difference of boundary condition between lateral boundary and free upper surfaces causes a singularity at the intersection (contact line). To avoid singularities on the contact line, we adopt weighted Sobolev spaces.     

Variational FORMLATION OF STEADY FLOW IN A HOPPER

A variational problem about maximal stable arches in a hopper is formulated.This problem idealizes an industrial problem related to guaranteeing reliable flow of material out of a storage silo. To obtain existence, the generalized function spaces are introduced and studied.Specifically, functions in the spaces can''be discontinuous in the ihteriOr of the''domain as well as along the boundaries. For the von Mises type of material in two dimension, the limit load is, estimated and its asymptotic behavior is investigated.     

The stress intensity factor for non-smooth fractures in antiplane elasticity

Motivated by some questions arising in the study of quasistatic growth in brittle fracture, we investigate the asymptotic behavior of the energy of the solution u of a Neumann problem near a crack in dimension 2. We consider non smooth cracks K that are merely closed and connected. At any point of density 1/2 in K, we show that the blow-up limit of u is the usual “cracktip” function ${C\sqrt{r}\sin(\theta/2)}$ , with a well-defined coefficient (the “stress intensity factor” or SIF). The method relies on Bonnet’s monotonicity formula (Bonnet, Variational methods for discontinuous structures, pp. 93–103. Birkhäuser, Basel, 1996) together with Γ-convergence techniques.     

Formulation and well-posedness of the Cauchy problem for a diffusion equation with discontinuous degenerating coefficients

The choice of a differential diffusion operator with discontinuous coefficients that corresponds to a finite flow velocity and a finite concentration is substantiated. For the equation with a uniformly elliptic operator and a nonzero diffusion coefficient, conditions are established for the existence and uniqueness of a solution to the corresponding Cauchy problem. For the diffusion equation with degeneration on a half-line, it is proved that the Cauchy problem with an arbitrary initial condition has a unique solution if and only if there is no flux from the degeneration domain to the ellipticity domain of the operator. Under this condition, a sequence of solutions to regularized problems is proved to converge uniformly to the solution of the degenerate problem in L ₁(R) on each interval.     

Winding behaviour of finite-time singularities of the harmonic map heat flow<Superscript>*</Superscript>

We settle a number of questions about the possible behaviour of the harmonic map heat flow at finite-time singularities. In particular, we show that a type of nonuniqueness of bubbles can occur at finite time, we show that the weak limit of the flow at the singular time can be discontinuous, we determine exactly the (polynomial) rate of blow-up in one particular example, and we show that winding behaviour of the flow can lead to an unexpected failure of convergence when the flow is (locally) lifted to the universal cover of the target manifold.Partly supported by an EPSRC Advanced Research Fellowship.     

A new method to exploit the Entropy Principle and galilean invariance in the macroscopic approach of Extended Thermodynamics

Abstract In extended thermodynamic the entropy principle and the Galilean invariance dictate respectively constraints for the constitutive equations and the velocity dependence. The entropy principle in particular requires the existence of a privileged field, the main field u^′, such that the original system becomes symmetric hyperbolic and is generated by four potentials. It is not easy to solve the restrictions of both principles, if we use as field the non convective main field and the velocity v. This is due to the fact that are not independent. Rather its components satisfy three scalar constraints. The aim of this paper is to solve the full problem using as new strategy to consider as independent variables and requiring an appropriate differential constraint. This new procedure is very efficient and we are able to solve the problem of 13 moments in the full non linear case (far from equilibrium). It turns out that the knowledge of only the equilibrium state function is sufficient to close the system. Keywords: Extended Thermodynamics, Entropy Principle, Galilean invariance, Rarefied Gas, Hyperbolic systems Mathematics Subject Classification (2000): 74A20, 76P05, 35l60     

Numerical and asymptotic study of non‐axisymmetric magnetohydrodynamic boundary layer stagnation‐point flows

Both numerical and asymptotic analyses are performed to study the similarity solutions of three‐dimensional boundary‐layer viscous stagnation point flow in the presence of a uniform magnetic field. The three‐dimensional boundary‐layer is analyzed in a non‐axisymmetric stagnation point flow, in which the flow is developed because of influence of both applied magnetic field and external mainstream flow. Two approaches for the governing equations are employed: the Keller‐box numerical simulations solving full nonlinear coupled system and a corresponding linearized system that is obtained under a far‐field behavior and in the limit of large shear‐to‐strain‐rate parameter (λ). From these two approaches, the flow phenomena reveals a rich structure of new family of solutions for various values of the magnetic number and λ. The various results for the wall stresses and the displacement thicknesses are presented along with some velocity profiles in both directions. The analysis discovered that the flow separation occurs in the secondary flow direction in the absence of magnetic field, and the flow separation disappears when the applied magnetic field is increased. The flow field is divided into a near‐field (due to viscous forces) and far‐field (due to mainstream flows), and the velocity profiles form because of an interaction between two regions. The magnetic field plays an important role in reducing the thickness of the boundary‐layer. A physical explanation for all observed phenomena is discussed. Copyright © 2017 John Wiley & Sons, Ltd.     

Strong convergence of the discontinuous Galerkin scheme for the low regularity miscible displacement equations

Strong convergence of the numerical solution to a weak solution is proved for a nonlinear coupled flow and transport problem arising in porous media. The method combines a mixed finite element method for the pressure and velocity with an interior penalty discontinuous Galerkin method in space for the concentration. Using functional tools specific to broken Sobolev spaces, the convergence of the broken gradient of the numerical concentration to the weak solution is obtained in the L² norm. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 489–513, 2017     

New stress and velocity fields for highly frictional granular materials

The idealized theory for the quasi-static flow of granular materialswhich satisfy the Coulomb–Mohr hypothesis is considered.This theory arises in the limit as the angle of internal frictionapproaches /2, and accordingly these materials may be referredto as being ‘highly frictional’. In this limit,the stress field for both two-dimensional and axially symmetricflows may be formulated in terms of a single nonlinear second-orderpartial differential equation for the stress angle. To obtainan accompanying velocity field, a flow rule must be employed.Assuming the non-dilatant double-shearing flow rule, a furtherpartial differential equation may be derived in each case, thistime for the streamfunction. Using Lie symmetry methods, a completeset of group-invariant solutions is derived for both systems,and through this process new exact solutions are constructed.Only a limited number of exact solutions for gravity-drivengranular flows are known, so these results are potentially importantin many practical applications. The problem of mass flow througha two-dimensional wedge hopper is examined as an illustration.     

Finite element approximation of viscoelastic fluid flow: Existence of approximate solutions and error bounds

Summary We study a finite element approximation of viscoelastic fluid flow obeying an Oldroyd B type constitutive law. The approximate stress, velocity and pressure are respectivelyP ₁ discontinuous,P ₂ continuous,P ₁ continuous. We use the method of Lesaint-Raviart for the convection of the extra stress tensor. We suppose that the continuous problem admits a sufficiently smooth and sufficiently small solution. We show by a fixed point method that the approximate problem has a solution and we give an error bound.This work has been supported in part by the GDR CNRS 901 Rhéologie der polymères fondus.     

Self-Similarity and the Evaluation of Long-Time Nonhydrostatic Effects in Compositionally Driven Gravity Flows in Deep Surroundings

In the study of compositionally driven gravity currents involving one or more homogeneous fluid layers, it has been customary to adopt the hydrostatic assumption for the pressure field in each layer which, in turn, leads to a depth‐independent horizontal velocity field in each of these layers and significant simplifications to the governing equations. Under this hydrostatic paradigm, each layer will then have its motion governed by the well‐known reduced dimension shallow‐water equations. For the so‐called ‐layer or reduced gravity shallow‐water equations, similarity solutions for fixed volume gravity currents released in rectangular geometry have been found. Very few attempts have been made to evaluate contributions arising from the possible loss of hydrostatic balance in the context of the problems treated using the classic shallow‐water approach. Where such attempts have been pursued, they have usually been carried out in a time‐independent context or using layer‐averaged equations and very small amplitude disturbances. The vast majority of these studies into nonhydrostatic effects do not include any relevant numerical work to assess these effects. In this paper, we develop an approach for evaluating nonhydrostatic contributions to the flow field for bottom gravity currents in deep surroundings and rectangular geometry. Our approach makes no assumptions on the amplitudes of the disturbances and does not depend on layer‐averaging in the governing equations. We seek asymptotic expansions of the solutions to the Euler equations for a shallow fluid by using the small parameter δ², where δ is the aspect ratio of the flow regime. At leading order the equations enforce hydrostatic balance while those obtained at first order retain certain nonhydrostatic effects which we evaluate. Our method for evaluation of these first‐order contributions employs the self‐similar nature of the solution to the leading‐order equations in the new first‐order equations without any vertical averaging procedures being employed.     

A new shallow water model with polynomial dependence on depth

In this paper, we study two‐dimensional Euler equations in a domain with small depth. With this aim, we introduce a small non‐dimensional parameter ε related to the depth and we use asymptotic analysis to study what happens when ε becomes small. We obtain a model for ε small that, after coming back to the original domain, gives us a shallow water model that considers the possibility of a non‐constant bottom, and the horizontal velocity has a dependence on z introduced by the vorticity when it is not zero. This represents an interesting novelty with respect to shallow water models found in the literature. We stand out that we do not need to make a priori assumptions about velocity or pressure behaviour to obtain the model. The new model is able to approximate the solutions to Euler equations with dependence on z (reobtaining the same velocities profile), whereas the classic model just obtains the average velocity. Copyright © 2007 John Wiley & Sons, Ltd.     

Nonlinear dynamic response for functionally graded shallow spherical shell under low velocity impact in thermal environment

Based on Giannakopoulos’s 2-D functionally graded material (FGM) contact model, a modified contact model is put forward to deal with impact problem of the functionally graded shallow spherical shell in thermal environment. The FGM shallow spherical shell, having temperature dependent material property, is subjected to a temperature field uniform over the shell surface but varying along the thickness direction due to steady-state heat conduction. The displacement field and geometrical relations of the FGM shallow spherical shell are established on the basis of Timoshenko–Midlin theory. And the nonlinear motion equations of the FGM shallow spherical shell under low velocity impact in thermal environment are founded in terms of displacement variable functions. Using the orthogonal collocation point method and the Newmark method to discretize the unknown variable functions in space and in time domain, the whole problem is solved by the iterative method. In numerical examples, the contact force and nonlinear dynamic response of the FGM shallow spherical shell under low velocity impact are investigated and effects of temperature field, material and geometrical parameters on contact force and dynamic response of the FGM shallow spherical shell are discussed.     

A combined hybrid mixed element method for incompressible miscible displacement problem with local discontinuous Galerkin procedure

In this article, we propose a combined hybrid discontinuous mixed finite element method for miscible displacement problem with local discontinuous Galerkin method. Here, to obtain more accurate approximation and deal with the discontinuous case, we use the hybrid mixed element method to approximate the pressure and velocity, and use the local discontinuous Galerkin finite element method for the concentration. Compared with other combined methods, this method can improve the efficiency of computation, deal with the discontinuous problem well and keep local mass balance. We study the convergence of this method and give the corresponding optimal error estimates in L^∞(L²) for velocity and concentration and the super convergence in L^∞(H¹) for pressure. Finally, we also present some numerical examples to confirm our theoretical analysis.     

Biomagnetic fluid flow over a stretching sheet with non linear temperature dependent magnetization

The flow of a heated ferrofluid over a linearly stretching sheet is studied in the pres- ence of an applied magnetic field due to a magnetic dipole. It is assumed that the applied magnetic field is sufficiently strong to saturate the ferrofluid and the variation of magnetization with temperature can be approximated by a non linear function of temperature difference. By introducing appropriate non dimensional variables the problem is described by a coupled and non linear system of ordinary differential equations with its boundary conditions which is solved numerically by applying an efficient numerical technique based on the common finite difference method. The obtained results are presented graphically for different values of the parameters entering into the problem under consideration and the dependence of the flow field from these parameters is discussed. A comparative study, with a similar problem which has already been solved and documented in literature, is also made wherever necessary, emphasizing the impor- tance of the non-linear variation of magnetization with temperature. Emphasis is also given in the obtained results for Prandtl number equal to 21 and critical exponent = 0.368 which are important and interesting in Biomagnetic Fluid Dynamics.     

A weak Galerkin finite element method for the Oseen equations

In this paper, a weak Galerkin finite element method for the Oseen equations of incompressible fluid flow is proposed and investigated. This method is based on weak gradient and divergence operators which are designed for the finite element discontinuous functions. Moreover, by choosing the usual polynomials of degree i ≥ 1 for the velocity and polynomials of degree i ? 1 for the pressure and enhancing the polynomials of degree i ? 1 on the interface of a finite element partition for the velocity, this new method has a lot of attractive computational features: more general finite element partitions of arbitrary polygons or polyhedra with certain shape regularity, fewer degrees of freedom and parameter free. Stability and error estimates of optimal order are obtained by defining a weak convection term. Finally, a series of numerical experiments are given to show that this method has good stability and accuracy for the Oseen problem.     

Vertical slot fishways: Mathematical modeling and optimal management

Fishways are the main type of hydraulic devices currently used to facilitate migration of fish past obstructions (dams, waterfalls, rapids,…$\mathrm{rapids},\dots$) in rivers. In this paper we present a mathematical formulation of an optimal control problem related to the optimal management of a vertical slot fishway, where the state system is given by the shallow water equations, the control is the flux of inflow water, and the cost function reflects the need of rest areas for fish and of a water velocity suitable for fish leaping and swimming capabilities. We give a first-order optimality condition for characterizing the optimal solutions of this problem. From a numerical point of view, we use a characteristic-Galerkin method for solving the shallow water equations, and we use an optimization algorithm for the computation of the optimal control. Finally, we present numerical results obtained for the realistic case of a standard nine pools fishway.     

Biomagnetic fluid flow over a stretching sheet with non linear temperature dependent magnetization

The flow of a heated ferrofluid over a linearly stretching sheet is studied in the pres- ence of an applied magnetic field due to a magnetic dipole. It is assumed that the applied magnetic field is sufficiently strong to saturate the ferrofluid and the variation of magnetization with temperature can be approximated by a non linear function of temperature difference. By introducing appropriate non dimensional variables the problem is described by a coupled and non linear system of ordinary differential equations with its boundary conditions which is solved numerically by applying an efficient numerical technique based on the common finite difference method. The obtained results are presented graphically for different values of the parameters entering into the problem under consideration and the dependence of the flow field from these parameters is discussed. A comparative study, with a similar problem which has already been solved and documented in literature, is also made wherever necessary, emphasizing the impor- tance of the non-linear variation of magnetization with temperature. Emphasis is also given in the obtained results for Prandtl number equal to 21 and critical exponent = 0.368 which are important and interesting in Biomagnetic Fluid Dynamics.