Optimal control of unsteady compressible viscous flows

The control of complex, unsteady flows is a pacing technology for advances in fluid mechanics. Recently, optimal control theory has become popular as a means of predicting best case controls that can guide the design of practical flow control systems. However, most of the prior work in this area has focused on incompressible flow which precludes many of the important physical flow phenomena that must be controlled in practice including the coupling of fluid dynamics, acoustics, and heat transfer. This paper presents the formulation and numerical solution of a class of optimal boundary control problems governed by the unsteady two‐dimensional compressible Navier–Stokes equations. Fundamental issues including the choice of the control space and the associated regularization term in the objective function, as well as issues in the gradient computation via the adjoint equation method are discussed. Numerical results are presented for a model problem consisting of two counter‐rotating viscous vortices above an infinite wall which, due to the self‐induced velocity field, propagate downward and interact with the wall. The wall boundary control is the temporal and spatial distribution of wall‐normal velocity. Optimal controls for objective functions that target kinetic energy, heat transfer, and wall shear stress are presented along with the influence of control regularization for each case. Copyright © 2002 John Wiley & Sons, Ltd.     

On estimating pressure build-up in severely accelerated,shocked flow

A new approach to the calculation of the high pressures characterizing the flow field in front of a piston undergoing severe acceleration over the short term is presented. In contrast with previous approaches where the computational domain is altered but which stop short of transforming velocities, here the problem is solved in an accelerating non-Euclidean co-ordinate system where the piston is stationary. The method is applied to a study of the problem of premature sabot separation. Through use of Harten's second-order-accurate TVD scheme, flow simulations are performed for both 1D and 3D axisymmetric geometries. The simple 1D model gives pressure profiles surprisingly close to those of the more physical 3D model.     

Influence of Relative Permeability on the Stability Characteristics of Immiscible Flow in Porous Media

Relative permeability functions for immiscible displacements in porous media show a wide range of profiles. Although, this behavior is well known, its impact on the stability of the displacement process is unexplored. Our analysis clearly demonstrates for the first time that the viscous instability characteristics of two-phase flows are governed not only by their end point values, but are strongly dependent on the actual profile of relative permeability functions. Linear stability analysis predicts the capacity of the flow to develop large scale fingers which can result in substantial bypassing of the resident fluid. It is observed that relative permeability functions attributed to drainage processes yield a more unstable displacement as compared to functions related to imbibition processes. Moreover, instability is observed to increase for those relative permeability functions which result from increased wettability of the wetting fluid. High accuracy numerical simulations show agreement with these predictions and demonstrate how large amplitude viscous fingers result in significant bypassing for certain relative permeability functions. In the nonlinear regime, the finger amplitude grows at a rate ∝ t^1/2 initially, drops to t^1/4 at a later time and finally grows ∝ t. The basic mechanisms of finger interaction, however, are not substantially influenced by relative permeability functions.     

Measurements and simulations of particle dispersion in a turbulent flow

A particle imaging technique has been used to collect droplet displacement statistics in a round turbulent jet of air. Droplets are injected on the jet axis, and a laser sheet and position-sensitive photomultiplier tube are used to track their radial displacement and time-of-flight. Dispersion statistics can be computed which are Lagrangian or Eulerian in nature. The experiments have been simulated numerically using a second-order closure scheme for the jet and a stochastic simulation for the particle trajectories. Results are presented for non-vaporizing droplets of sizes from 35 to 160 μm. The simulations have underscored the importance of initial conditions and early droplet displacement history on the droplet trajectory for droplets with large inertia relative to the turbulence. Estimates of initial conditions have been made and their effect on dispersion is quantified.     

Averaged topological equations for dispersed two-phase flows

A general relationship between the volume fraction and the specific interfacial area for averaged dispersed two-phase flows is proposed. This relationship, expressed as a basic set of two scalar evolution equations and two vectorial non-uniformity state equations, is an analytical result obtained by a systematic approach using the derivatives of some generalized functions and a local volume-averaging technique. The proposed set of equations was expressed for measurable macroscopic parameters of the system and has the same generality as the averaged transport equations of two-phase flows. By combination of the basic set of equations, called the averaged topological equations (ATEs), second-order ATEs for the volume fraction were found. The second-order ATEs were expressed both by a Lagrangian formulation and by a Eulerian formulation. The importance and physical meaning of the ATEs developed in this study were clarified within the framework of the theory of kinematic waves.     

Modélisation du champ des vitesses de l'urine dans un bolus urétéralModelling of urine flow in an ureteral bolus

Urine transport is made from the kidney to the bladder through the ureter by isolated pockets called bolus. To determine the urine flow in a bolus, we use an adherence condition on the interface urine/wall. It gives us an infinite linear system verified by a set of parameters. An iterative and convergent algorithm allows us to solve this system and to determine analytically the components of the velocity vector in the bolus. To cite this article: A. Vogel et al., C. R. Mecanique 332 (2004).     

Numerical method for calculation of the incompressible flow in general curvilinear co‐ordinates with double staggered grid

A solution methodology has been developed for incompressible flow in general curvilinear co‐ordinates. Two staggered grids are used to discretize the physical domain. The first grid is a MAC quadrilateral mesh with pressure arranged at the centre and the Cartesian velocity components located at the middle of the sides of the mesh. The second grid is so displaced that its corners correspond to the centre of the first grid. In the second grid the pressure is placed at the corner of the first grid. The discretized mass and momentum conservation equations are derived on a control volume. The two pressure grid functions are coupled explicitly through the boundary conditions and implicitly through the velocity of the field. The introduction of these two grid functions avoids an averaging of pressure and velocity components when calculating terms that are generated in general curvilinear co‐ordinates. The SIMPLE calculation procedure is extended to the present curvilinear co‐ordinates with double grids. Application of the methodology is illustrated by calculation of well‐known external and internal problems: viscous flow over a circular cylinder, with Reynolds numbers ranging from 10 to 40, and lid‐driven flow in a cavity with inclined walls are examined. The numerical results are in close agreement with experimental results and other numerical data. Copyright © 2003 John Wiley & Sons, Ltd.     

承受拉-扭组合载荷的薄壁圆筒屈服时产生塑性流动的条件

提出了在无外力作功的情况下，具有Bauschinger效应的弹塑性材料处于屈服状态产生自发的塑性流动时应满足的条件．这个条件不仅与材料的力学性能有关，而且还处决于材料的具体的载荷边界条件和变形．举例说明了承受拉一扭组合的薄壁圆筒中，采用组合强化模型时，产生塑性流动的具体条件．     

Effect of biotransformation on multispecies plume evolution and natural attenuation

Biological transformation of volatile organic compounds is one of the key factors that influence contaminant-plume evolution and thus natural attenuation. In this study we investigate the effect of biological transformation on the transport of contaminants in the aqueous and gaseous phases. The analysis includes the study of the effect of density-driven advection of contaminants in the gaseous phase on multiphase and multispecies flow, fate and transport modeling in the subsurface. Trichloroethylene (TCE) and its two byproducts, dichloroethylene and vinyl chloride, are analyzed as the target contaminants. Our results indicate that density-driven advection of the gaseous phase, which is initiated by evaporation of TCE as a nonaqueous phase liquid, increases the downward and also the lateral migration of TCE within the unsaturated zone. This process also influences the location of high-concentration zones of the byproducts of TCE in the unsaturated and the saturated zones. Biotransformation of TCE contributes to the reduction of dissolved TCE plume development as expected. The daughter byproducts, which are introduced into the subsurface system, show distinct transport patterns as they are affected by their independent degradation kinetics and density-driven advection. These observations, which are based on our simulation results for biotransformation and transport of TCE and its byproducts, are useful in evaluating the natural attenuation processes, its potential health hazards and also the evaluation of potential plume development at contaminated sites.