A multiphase first order model for non-equilibrium sand erosion,transport and sedimentation

Three phenomena are involved in sand movement: erosion, wind transport, and sedimentation. This paper presents a comprehensive easy-to-use multiphase model that include all three aspects with a particular attention to situations in which erosion due to wind shear and sedimentation due to gravity are not in equilibrium. The interest is related to the fact that these are the situations leading to a change of profile of the sand bed.     

热扩散和电场对风沙流发展过程的影响

建立了描述在来流水平风场和由于地表热扩散产生的垂向风场联合作用下的风沙流发展过程的基本方程．通过定量分析表明:地表热扩散和风沙电场对风沙流发展过程的影响十分明显．在此基础上,给出地表热扩散和风沙电场对风沙流中的沙粒跃移轨迹、输沙率以及风沙流达到充分发展的时间等的影响规律．     

基流风速包含垂直变化时的三维无旋Non-Boussinesq流的地形重力波拖曳解

用WKB近似方法建立了表达三维地形重力波拖曳的解析Non-Boussinesq扰动模型,其中在大Richardson数条件下给出了(静力和非静力模型的)重力波拖曳及其地表扰动气压的二阶表达式.通过针对经典的理想化三维圆钟型山体的一个算例证明,当基流风速切变为线性时,重力波拖曳随着切变的增强而减弱;并且前向垂直切变(forward-shear,风速随高度增加)所对应的重力波拖曳比反向切变(backward-shear,风速随高度减小)所对应的重力波拖曳减弱得更快.这种现象与模型是否采用静力近似无关.     

Steady state numerical simulation of the particle collection efficiency of a new urban sustainable gravity settler using design of experiments by FVM

The aim of this work is to analyze the efficiency of a new sustainable urban gravity settler to avoid the solid particle transport, to improve the water waste quality and to prevent pollution problems due to rain water harvesting in areas with no drainage pavement. In order to get this objective, it is necessary to solve particle transport equations along with the turbulent fluid flow equations since there are two phases: solid phase (sand particles) and fluid phase (water). In the first place, the turbulent flow is modelled by solving the Reynolds-averaged Navier-Stokes (RANS) equations for incompressible viscous flows through the finite volume method (FVM) and then, once the flow velocity field has been determined, representative particles are tracked using the Lagrangian approach. Within the particle transport models, a particle transport model termed as Lagrangian particle tracking model is used, where particulates are tracked through the flow in a Lagrangian way. The full particulate phase is modelled by just a sample of about 2,000 individual particles. The tracking is carried out by forming a set of ordinary differential equations in time for each particle, consisting of equations for position and velocity. These equations are then integrated using a simple integration method to calculate the behaviour of the particles as they traverse the flow domain. The entire FVM model is built and the design of experiments (DOE) method was used to limit the number of simulations required, saving on the computational time significantly needed to arrive at the optimum configuration of the settler. Finally, conclusions of this work are exposed.     

SOIL EROSION AND DEGRADATION BASED ON SAND PARTICLES TRANSPORT CAUSED BY WIND BLOWING

Abstract In this paper, the effect of sand particles transport caused by wind blowing and its role in the land degradation and desertification process is considered. For the modeling of the 3D landscape, a grayscale height map has been used, the vegetation has been modeled using a Lindenmayer system, and the sand particles have been modeled as a 3D mesh‐free particles system. It was assumed that both the sand motion and the wind motion are incompressible continuum systems and their behavior follows the Navier–Stokes equations. To simulate the sand transport, the Navier–Stokes equations are discretized using the moving particle Semi‐implicit (MPS) method. Different types of revegetation patterns (windbreakers) have been used to show some effective measures preventing soils from erosion.     

Solvability of a supercritical free surface flow problem under gravity-capillarity

In this paper, we consider a problem of a supercritical free surface flow over an obstacle lying on the bottom of a channel in 2D. The flow is irrotational, stationary and the fluid is ideal and incompressible. We take into account both the gravity and the effects of the superficial tension. The problem is nonlinear, it is formulated by the Laplace operator and the dynamic condition defined on the free surface of the fluid domain (Bernoulli equation). Using the perturbation stream function, we linearize the problem and we give a priori properties of the solution. These a priori properties allow us to construct a space where we can use the Lax–Milgram’s theorem to prove the existence and the uniqueness of the solution of the problem.     

Fluid accumulation in thin-film flows driven by surface tension and gravity

In this note we derive a model describing the two-dimensional viscous flow driven by surface tension and gravity of a thin liquid film near a stagnation point. In the thin-film approximation of such a flow, accumulation takes place where the combined effects of gravity and surface tension stop the flow. These stagnation points are characterised to leading order by the geometry of the substrate. We first derive the thin-film approximation that describes the flow away from such accumulation regions. Then, assuming the existence of isolated stagnation points, we derive the boundary layer equation describing the inner structure of solutions describing accumulation. The existence of these solutions has been proved by the authors elsewhere. Finally, in order to justify the model we prove the existence of curves that give a substrate with an isolated stagnation point.     

Stabilization of periodic Stokesian Hele–Shaw flows of ferrofluids

We consider an incompressible ferrofluid in a vertical Hele–Shaw cell and develop a proper analytic framework for the free interface and the velocity potential of the fluid in a periodic geometry. The flow is assumed to obey a non-Newtonian Darcy law. The forces influencing the fluid are gravity, surface tension and the response to a magnetic field induced by a current. In addition, the flow is stabilized at the lower boundary component by an external source b. We prove a well-posedness result for the flow near flat solutions. Moreover, we find conditions on the parameters and on the slope of b for the exponential stability and instability of flat interfaces. Furthermore, we identify values for the current's intensity ι where critical bifurcation of nontrivial finger-shaped solutions from the branch of trivial (flat) solutions takes place.     

Modeling and control of surface gravity waves in a model of a copper converter

Undesirable splashing appears in copper converters when air is injected into the molten matte to trigger the conversion process. We consider here a cylindrical container horizontally placed and containing water, where gravity waves on the liquid surface are generated due to water injection through a lateral submerged nozzle. The fluid dynamics in a transversal section of the converter is modeled by a 2-D inviscid potential flow involving a gravity wave equation with local damping on the liquid surface. Once the model is established, using a finite element method, the corresponding natural frequencies and normal modes are numerically computed in the absence of injection, and the solution of the system with injection is obtained using the spectrum. If a finite number of modes is considered, this approximation leads to a system of ordinary differential equations where the input is represented by the fluid injection. The dynamics is simulated as perturbations around a constant fluid injection solution, which is the desired operating state of the system, considering that the conversion process does not have to be stopped or seriously affected by the control. The solution is naturally unstable without control and the resulting increase of amplitude of the surface waves are assimilable to the splashing inside the converter. We show numerically that a variable flow around the operating injection is able to sensibly reduce these waves. This control is obtained by a LQG feedback law by measuring the elevation of the free surface at the point corresponding to the opposite extreme to where the nozzle injection is placed.     

On a mathematical model for laser-induced thermotherapy

We study a mathematical model for laser-induced thermotherapy, a minimally invasive cancer treatment. The model consists of a diffusion approximation of the radiation transport equation coupled to a bio-heat equation and a model to describe the evolution of the coagulated zone. Special emphasis is laid on a refined model of the applicator device, accounting for the effect of coolant flow inside. Comparisons between experiment and simulations show that the model is able to predict the experimentally achieved temperatures reasonably well.     

明渠层流失稳与沙纹成因机理研究

动床水流中,泥沙起动之后,往往要出现沙纹,沙纹成因各家的解释不一。白玉川,罗纪生的观点是:沙纹的尺度较小,主要是由于明渠层流不稳定性波或床面近壁流层中小尺度拟序结构发展演化所致。当床面边界附近扰动波或拟序结构以及水流自身所产生的床面底部切应力大于Shields切应力后,床面即产生响应,形成沙纹;如果扰动所产生的扰动切应力频率接近床面泥沙固有频率,则产生与泥沙颗粒的共振,这种现象也称之为“泥沙的检波性质”,此时床面发生最大响应,沙纹发展速度也最快。     

Resonant wind-driven oceanic motions

We are interested here in describing the linear response of a rapidly rotating fluid to some surface stress, possibly due to the wind. The distinctive feature of the model considered here lies in the fact that the stress admits fast time oscillations and may be resonant with the Coriolis force. In addition to the usual Ekman layer, the model exhibits another – much larger – boundary layer, and some global vertical profile. We prove, in particular, that for large times, the wind effect is no longer localized in the vicinity of the surface. To cite this article: A.-L. Dalibard, L. Saint-Raymond, C. R. Acad. Sci. Paris, Ser. I 347 (2009).     

A fully Eulerian multiphase model of windblown sand coupled with morphodynamic evolution: Erosion,transport, deposition,and avalanching

Modeling unsteady windblown sand dynamics requires not only treatment of the sand present in the air as a suspended constituent of a mixture but also consideration of erosion and sedimentation phenomena and consequently of the morphodynamic evolution of the sand-bed surface, including avalanching, especially in the presence of natural or human-built obstacles, artifacts, and infrastructures. With this aim in mind, we present a comprehensive multiphase model capable of accurately simulating all the physical phenomena mentioned above, producing satisfactory results, with reasonable computational effort. As test cases, two- and three-dimensional simulations of dune evolution are reported, as is windblown sand transport over a straight vertical wall. Examples of sand transport around other obstacles are given to show the flexibility of the model and its usefulness for such engineering applications.     

Micropolar fluid flow through the membrane composed of impermeable cylindrical particles coated by porous layer under the effect of magnetic field

In the present work, the magnetohydrodynamic flow of a micropolar fluid through the membrane composed of impermeable cylindrical particles coated by porous layer is considered. The flow of a fluid is taken parallel to an axis of cylinder and a uniform magnetic field is applied in transverse direction of the flow. The problem is solved by using the cell model technique for the flow through assemblage of cylindrical particles. The solution of the problem has been obtained by using no-slip condition, continuity of velocity and stresses at interfaces along with Happle's no-couple stress condition as the boundary conditions. The expressions for the linear velocity, micro-rotational velocity, flow rate and hydrodynamic permeability of the membrane are achieved in this work. The obtained solution for velocities is used to plot the graph against various transport parameters such as, Hartmann number, coupling parameter, porosity, scaling parameter etc. The effect of these transport parameters on the flow velocity, micro-rotational velocity, and the hydrodynamic permeability of the membrane have been presented and discussed in this work.     

Numerical modeling of hydrodynamic and hydrothermal characteristics in subtropical alpine lake

A three-dimensional, time-dependent hydrodynamic and hydrothermal model was performed and applied to the subtropical alpine Yuan-Yang Lake (YYL) in northeastern region of Taiwan. The model was driven with discharge inflow, heat, and wind stress to simulate the hydrodynamic and hydrothermal in the lake. The model was validated with measured water surface elevation, current, and temperature in 2008. The overall model simulation results are in quantitative agreement with the available field data. The validated model was then used to investigate wind-driven current, mean circulation, and residence time in the YYL. The modeling results reveal that the velocity field along the wind axis present the variations over depth with return current where the velocity at the surface layer is along the wind direction while it is opposite near 1 m below water surface. The simulated mean current indicates that the surface currents flow towards the southwest direction and form a clock-wise rotation. The calculated residence time is strongly dependent on the inflows and wind effects. Regression analysis of model results reveals that an exponential regression equation can be employed to correlate the residence time to change of discharge input. The residence time without wind stress is higher than that with wind effect, indicating that wind plays an important role in lake mixing. The calculated residence time is approximately 2-2.5 days under low inflow with wind effect.     

Optimization methods for pipeline transportation of natural gas with variable specific gravity and compressibility

In this paper, the problem of flow maximization in pipeline systems for transmission of natural gas is addressed. We extend previously suggested models by incorporating the variation in pipeline flow capacities with gas specific gravity and compressibility. Flow capacities are modeled as functions of pressure, compressibility and specific gravity by the commonly-used Weymouth equation, and the California Natural Gas Association method is used to model compressibility as a function of specific gravity and pressure. The sources feeding the transmission network do not necessarily supply gas with equal specific gravity. In our model, it is assumed that when different flow streams enter a junction point, the specific gravity of the resulting flow is a weighted average of the specific gravities of entering flows. We also assume the temperature to be constant, and the system to be in steady state. Since the proposed model is non-convex, and global optimization hence can be time consuming, we also propose a heuristic method based on an iterative scheme in which a simpler NLP model is solved in each iteration. Computational experiments are conducted in order to assess the computability of the model by applying a global optimizer, and to evaluate the performance of the heuristic approach. When applied to a wide set of test instances, the heuristic method provides solutions with deviations less than 10% from optimality, and in many instances turns out to be exact. We also report several experiments demonstrating that letting the compressibility and the specific gravity be global constants can lead to significant errors in the estimates of the total network capacity.     

Instability mechanisms in buoyant-thermocapillary liquid pools

A cylindrical volume of fluid, with a free surface on top, is heated by a parabolic heat flux from above. Two physical effects drive a flow: thermocapillary effects due to free-surface temperature gradients introduced by the non-uniform heat flux and buoyancy forces due to gravity. The basic axisymmetric flow is computed by finite volumes and its stability is investigated by a linear-stability analysis. It is found that the critical stability boundaries and modes are similar to those known from the half-zone model of crystal growth. For low Prandtl numbers the critical mode is steady and three-dimensional. We find an asymptotic critical value in the limit of vanishing Prandtl number. For increasing Prandtl number the critical Reynolds number increases. Near unit Prandtl number no threshold could be found with the present computational limitations. For Prandtl numbers larger than unity, the critical mode is oscillatory and the critical Reynolds number decreases with the Prandtl number. We present evidence that the low- and high-Prandtl-number instabilities are essentially centrifugal respectively due to the hydrothermal-wave mechanism. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical modelling of salinity variations due to wind and thermohaline forcing in the Persian Gulf

Salinity is an important component of the marine system. Due to shallow nature of the Persian Gulf, the salinity has been influenced by both wind driven and surface thermohaline fluxes (heat and moisture fluxes). In this study, the seasonal distribution of salinity and its variations due to wind stress and thermohaline forcing are investigated by using a three-dimensional hydrodynamic model, Coupled Hydrodynamical–Ecological Model for Regional and Shelf Seas (COHERENS). The simulation results show that the salinity in the Persian Gulf experiences dramatic spatial and temporal variations. The influence of the thermohaline forcing is considerably more than the wind stress on the salinity. The effect of the surface thermohaline fluxes over the salinity field is generally to increase the salinity for almost all the water column during the year. This effect is high during September–November where the evaporative surface salinity flux dominates over inflow of low-salinity values of Indian Ocean Surface Water. The wind forcing at the most regions of the Persian Gulf, in particular at the United Arab Emirate (UAE) coast and Bahrain–Qatar shelf, freshens the water all the year round. The wind and thermohaline forcing in March–June have strong potential to generate stratification in salinity structure. The model predictions, which are successful in simulating many features of observed pattern, indicate that the surface water of the Gulf is saltier in winter than that in spring and early summer. Both heat fluxes and wind stress play an important role for this seasonal cycle of the surface salinity.     

Well‐posedness of a generalized time‐harmonic transport equation for acoustics in flow

We study a generalized time‐harmonic transport equation, which appears in the Goldstein equations and allows us to model the acoustic radiation in a flow. We investigate the well‐posedness of this transport problem. The result will be established under the assumption of a Ω‐filling flow, which, in 2D, is simply equivalent to a flow that does not vanish. The approach relies on the method of characteristics, which leads to the resolution of the transport equation along the streamlines, and on general results of functional analysis. The theoretical results are illustrated with numerical results obtained with a Streamline Upwind Petrov‐Galerkin finite element scheme.     

Resonant Interactions between Vortical Flows and Water Waves. Part I: Deep Water

Any weak, steady vortical flow is a solution to leading order of the inviscid fluid equations with a free surface, so long as this flow has horizontal streamlines coinciding with the undisturbed free surface. This work considers the propagation of irrotational surface gravity waves when such a vortical flow is present. In particular, when the vortical flow and the irrotational surface waves are both periodic, resonant interactions can occur between the various components of the flow. The periodic vortical component of the flow is proposed as a model for more complicated vortical flows that would affect surface waves in the ocean, such as the turbulence in the wake of a ship. These resonant interactions are studied in two dimensions, both in the limit of deep water (Part I) and shallow water (Part II). For deep water, the resonant set of surface waves is governed by “triad-like” ordinary differential equations for the wave amplitudes, whose coefficients depend on the underlying rotational flow. These coefficients are calculated explicitly and the stability of various configurations of waves is discussed. The effect of three dimensionality is also briefly mentioned.