Local Null Controllability of the N-Dimensional Navier–Stokes System with N − 1 Scalar Controls in an Arbitrary Control Domain

In this paper we deal with the local null controllability of the N-dimensional Navier–Stokes system with internal controls having one vanishing component. The novelty of this work is that no condition is imposed on the control domain.     

Numerical Simulation of Delft-Jet-in-Hot-Coflow (DJHC) Flames Using the Eddy Dissipation Concept Model for Turbulence–Chemistry Interaction

In this paper, we report results of a numerical investigation of turbulent natural gas combustion for a jet in a coflow of lean combustion products in the Delft-Jet-in-Hot-Coflow (DJHC) burner which emulates MILD (Moderate and Intense Low Oxygen Dilution) combustion behavior. The focus is on assessing the performance of the Eddy Dissipation Concept (EDC) model in combination with two-equation turbulence models and chemical kinetic schemes for about 20 species (Correa mechanism and DRM19 mechanism) by comparing predictions with experimental measurements. We study two different flame conditions corresponding to two different oxygen levels (7.6% and 10.9% by mass) in the hot coflow, and for two jet Reynolds number (Re = 4,100 and Re = 8,800). The mean velocity and turbulent kinetic energy predicted by different turbulence models are in good agreement with data without exhibiting large differences among the model predictions. The realizable k-ε model exhibits better performance in the prediction of entrainment. The EDC combustion model predicts too early ignition leading to a peak in the radial mean temperature profile at too low axial distance. However the model correctly predicts the experimentally observed decreasing trend of lift-off height with jet Reynolds number. A detailed analysis of the mean reaction rate of the EDC model is made and as possible cause for the deviations between model predictions and experiments a low turbulent Reynolds number effect is identified. Using modified EDC model constants prediction of too early ignition can be avoided. The results are weakly sensitive to the sub-model for laminar viscosity and laminar diffusion fluxes.     

On a Generalized Nonlinear K–ε Model and the Use of Extended Thermodynamics in Turbulence

A resent extension of the nonlinear K–ε model is critically discussed from a basic theoretical standpoint. While it was said in the paper that this model was formulated to incorporate relaxation effects, it will be shown that the model is incapable of describing one of the most basic such turbulent flows as is obvious but is described for clarity. It will be shown in detail that this generalized nonlinear K–ε model yields erroneous results for the Reynolds stress tensor when the mean strains are set to zero in a turbulent flow – the return-to-isotropy problem which is one of the most elementary relaxational turbulent flows. It is clear that K–ε type models cannot describe relaxation effects. While their general formalism can describe relaxation effects, the nonlinear K–ε model – which the paper is centered on – cannot. The deviatoric part of the Reynolds stress tensor is predicted to be zero when it actually only gradually relaxes to zero. Since this model was formulated by using the extended thermodynamics, it too will be critically assessed. It will be argued that there is an unsubstantial physical basis for the use of extended thermodynamics in turbulence. The role of Material Frame-Indifference and the implications for future research in turbulence modeling are also discussed. Received 19 February 1998 and accepted 23 October 1998     

Experimental and Numerical Study of a Hypersonic Separated Flow in the Vicinity of a Cone‐Flare Model

An axisymmetric laminar separated flow in the vicinity of a coneflare model is studied experimentally and numerically for a Mach number M = 6. The distributions of pressure and Stanton numbers along the model surface and velocity profiles in the region of shock wave–boundary layer interaction are measured and compared with the calculated data. The influence of the laminar–turbulent transition on flow parameters is studied numerically.     

Local Exact Controllability for the One-Dimensional Compressible Navier–Stokes Equation

In this paper we deal with the isentropic (compressible) Navier-Stokes equation in one space dimension and we adress the problem of the boundary controllability for this system. We prove that we can drive initial conditions which are sufficiently close to some constant states to those constant states. This is done under some natural hypotheses on the time of control and on the regularity on the initial conditions.     

Complete Global and Bifurcation Analysis of a Stoichiometric Predator–Prey Model

Loladze et al. (Bull Math Biol 62:1137–1162, 2000) proposed a highly cited stoichiometric predator–prey system, which is nonsmooth, and thus it is extremely difficult to analyze its global dynamics. The main challenge comes from the phase plane fragmentation and parameter space partitioning in order to perform a detailed and complete global stability and bifurcation analysis. Li et al. (J Math Biol 63:901–932, 2011) firstly discussed its global dynamical behavior with Holling type I functional response and found that the system has no limit cycles, and the internal equilibrium is globally asymptotically stable if it exists. Secondly, for the system with Holling type II functional response, Li et al. (2011) fixed all parameters (with realistic values) except K to perform the bifurcation analysis and obtained some interesting phenomena, for instance, the appearance of bistability and many bifurcation types. The aim of this paper is to provide a complete global analysis for the system with Holling type II functional response without fixing any parameter. Our analysis shows that the model has far richer dynamics than those found in the previous paper (Li et al. 2011), for example, four types of bistability appear: besides the bistability between an internal equilibrium and a limit cycle as shown in Li et al. (2011), the other three bistabilities occur between an internal equilibrium and a boundary equilibrium, between two internal equilibria, or between a boundary equilibrium and a limit cycle. In addition, this paper rigorously provides all possible bifurcation passways of this stoichiometric model with Holling type II functional response.     

Development and Validation of a Numerical Model of Flow Through Embankment Dams – Comparisons with Experimental Data and Analytical Solutions

The development and validation of a numerical simulation model of the flow through embankment dams is described. The paper focuses on basic verification studies, that is, comparisons with analytical solutions and data from laboratory experiments. Two experimental studies, one dealing with the flow in a Hele–Shaw cell and the other with the flow through a bed of packed glass beads, are also described. Comparisons are carried out with respect to the phreatic surfaces, pressure profiles, seepage levels and discharges. It is concluded that the agreement between experimental, analytical and numerical results is generally satisfactory.     

Numerical simulation of a compressible vortex–wall interaction

The wall interaction of isolated compressible vortices generated from a short driver section shock tube has been simulated numerically by solving the Navier–Stokes equations in axisymmetric form. The dynamics of shock-free (incident shock Mach number \(M = 1.36\)) and shock-embedded \((M = 1.57)\) compressible vortices near the wall has been studied in detail. The AUSM+ scheme with a fifth-order upwind interpolation formula is used for the convective fluxes. Time integration is performed using a low dissipative and dispersive fourth-order six-stage Runge–Kutta scheme. The evolution of primary and wall vortices has been shown using the velocity field, vorticity field, and numerical schlierens. The vortex impingement, shocklets, wall vortices, and their lift-off are clearly identified from the wall pressure time history. It has been observed that the maximum vorticity of the wall vortices reaches close to 30 % of the primary vortex for \(M = 1.36\) and it reaches up to 60 % for \(M = 1.57\). The net pressure force on the wall due to incident shock impingement is dominant compared to the compressible vortex impingement and their evolution.     

Numerical study of gas–solid two-phase flow and erosion in a cavity with a slope

Gate valve is mainly used to turn on or turn off the pipeline in pneumatic conveying. When the gate valve is fully open, the particles are easy to collide with the cavity rear wall and enter into the cavity, resulting in particles’ accumulation in the cavity. The particles in cavity will accumulate between the cavity bottom and the flashboard bottom wall and prevent the gate from turning off normally. Meanwhile, the particles’ collision with cavity rear wall will cause serious erosion. Both the particles’ accumulation and erosion will cause the poor sealing of the gate valve, further resulting in the leakage of the pipeline system. To reduce the particles’ accumulation in cavity and erosion on cavity when the gate valve is fully open, we simplify the gate valve into a cavity structure and study it. We find that adding a slope upstream the cavity can effectively reduce the particles’ accumulation in the cavity and the erosion on the cavity rear wall. In this work, Eulerian–Lagrangian method in commercial code (FLUENT) was used to study the gas–solid two-phase flow and erosion characteristics of a cavity with a slope. The particle distribution shows that the particles with Stokes number St = 1.3 and St = 13 cannot enter the cavity due to the slope, but the particles with St = 0.13 enter the cavity following the gas. For St = 13, the particles collide with the wall many times in the ideal cavity. Erosion results show that the slope can transfer the erosion on cavity rear wall to the slope and reduce the maximum erosion rate of the wall near the cavity to some degrees.     

Reflection of a Dam–Break Wave at a Vertical Wall. Numerical Modeling and Experiment

Results of the experimental study and numerical modeling of the reflection of a dam–break wave at the vertical end wall of a channel are given. A wave forms with distance from a partition creating the initial level difference of the liquid. It is shown that a numerical calculation based on the Zheleznyak—Pelinovskii nonlinear dispersion model satisfactorily describes the height of the splash–up, the amplitude of reflected waves, and the wave velocity in front of the wall for smooth and dam–break waves. It is also shown that, for smooth and weakly breaking (without significant entrainment of air) incoming waves, the experimental values of the height of the splash–up at the wall agree well with relevant experimental and calculated data for solitary waves.     

Numerical Study of Hydrodynamics and Heat and Mass Transfer of a Ducted Gas–Vapor‐Droplet Flow

A computational model has been developed to predict heat and mass transfer and hydrodynamic characteristics of a turbulent gas–vapor–droplet flow. Turbulent characteristics of the gas phase are computed using the k– model of turbulence. It is shown that, with increasing inlet droplet diameter, the rate of heat transfer between the duct surface and the vapor–gas mixture decreases appreciably, whereas the wall friction increases only insignificantly. The predicted values agree fairly well with available experimental and numerical data     

Numerical simulation and stability analysis of a novel reaction–diffusion COVID-19 model

Nonlinear Dynamics - In this study, a novel reaction–diffusion model for the spread of the new coronavirus (COVID-19) is investigated. The model is a spatial extension of the recent COVID-19...     

Numerical modeling and experimental investigation of gas–liquid slug formation in a microchannel T-junction

Gas–liquid two-phase flow in a microfluidic T-junction with nearly square microchannels of 113 μm hydraulic diameter was investigated experimentally and numerically. Air and water superficial velocities were 0.018–0.791 m/s and 0.042–0.757 m/s, respectively. Three-dimensional modeling was performed with computational fluid dynamics (CFD) software FLUENT and the volume of fluid (VOF) model. Slug flow (snapping/breaking/jetting) and stratified flow were observed experimentally. Numerically predicted void fraction followed a linear relationship with the homogeneous void fraction, while experimental values depended on the superficial velocity ratio U_G/U_L. Higher experimental velocity slip caused by gas inlet pressure build-up and oscillation caused deviation from numerical predictions. Velocity slip was found to depend on the cross-sectional area coverage of the gas slug, the formation of a liquid film and the presence of liquid at the channel corners. Numerical modeling was found to require improvement to treat the contact angle and contact line slip, and could benefit from the use of a dynamic boundary condition to simulate the compressible gas phase inlet reservoir.     

Behaviors of Solutions for a Singular Prey–Predator Model and its Shadow System

We study the asymptotic behaviors and quenching of the solutions for a two-component system of reaction–diffusion equations modeling prey–predator interactions in an insular environment. First, we give a global existence result for the solutions to the corresponding shadow system. Then, by constructing some suitable Lyapunov functionals, we characterize the asymptotic behaviors of global solutions to the shadow system. Also, we give a finite time quenching result for the shadow system. Finally, some global existence results for the original reaction–diffusion system are given.     

Numerical simulations of gas–liquid–solid flows in a hydrocyclone separator

The flow behavior in hydrocyclones is quite complex. In this study, the computational fluid dynamics (CFD) method was used to simulate the flow fields inside a hydrocyclone in order to investigate its separation efficiency. In the computational fluid dynamics study of hydrocyclones, the air-core dimension is a key to predicting the mass split between the underflow and overflow. In turn, the mass split influences the prediction of the size classification curve. Three models, the model, the Reynolds stress model (RSM) without considering the air-core, and the Reynolds stress turbulence model with the volume of fluid (VOF) multiphase model for simulating the air-core, were compared in terms of their predictions of velocity, axial and tangential velocity distributions, and separation proportion. The RSM with air-core simulation model, since it reproduces some detailed features of the turbulence and multiphase, clearly predicted the experimental data more closely than did the other two models.     

Spatiotemporal Patterns of a Reaction–Diffusion Substrate–Inhibition Seelig Model

In this paper, the spatiotemporal patterns of a reaction–diffusion substrate–inhibition chemical Seelig model are considered. We first prove that this parabolic Seelig model has an invariant rectangle in the phase plane which attracts all the solutions of the model regardless of the initial values. Then, we consider the long time behaviors of the solutions in the invariant rectangle. In particular, we prove that, under suitable “lumped parameter assumption” conditions, these solutions either converge exponentially to the unique positive constant steady states or to the spatially homogeneous periodic solutions. Finally, we study the existence and non-existence of Turing patterns. To find parameter ranges where system does not exhibit Turing patterns, we use the properties of non-constant steady states, including obtaining several useful estimates. To seek the parameter ranges where system possesses Turing patterns, we use the techniques of global bifurcation theory. These two different parameter ranges are distinguished in a delicate bifurcation diagram. Moreover, numerical experiments are also presented to support and strengthen our analytical analysis.     

A Spalart–allmaras Turbulence Model Implementation for High-order Discontinuous Galerkin Solution of the Reynolds-averaged Navier-stokes Equations

A novel and robust approach has been proposed for the high-order discontinuous Galerkin (DG) discretization of the Reynolds-averaged Navier-Stokes (RANS) equations with the turbulence model of Spalart-Allmaras (SA). The solution polynomials of the SA equation are reconstructed by the Hermite weighted essentially non-oscillatory (HWENO) scheme. Several practical techniques are suggested to simplify and extend a positivity-preserving limiter to further guarantee the positivity of SA working variable. The resulting positivity-preserving HWENO limiting method is compact and easy to implement on arbitrary meshes. Typical turbulent flows are conducted to assess the accuracy and robustness of the present method. Numerical experiments demonstrate that with the increasing grid or order resolution, the limited results of the working variable are getting closer to the unlimited ones. And the most obvious improvement with proposed method is on the computation of the working variable field in wake regions.     

Numerical and Experimental Study of Rayleigh–Bénard–Kelvin Convection

We performed experimental and numerical studies of combined effects of thermal buoyancy and magnetization force applied on a cubical enclosure of a paramagnetic fluid heated from below and cooled from top. The temperature difference between the hot and cold wall was kept constant. After considering neutral situation (i.e. a pure natural convection case), magnetic fields of different intensity were imposed. The magnetization force produced significant changes in flow (transition from laminar to turbulent regimes), wall-heat transfer (enhancement) and turbulence (turbulence structures reorganization). The strong magnetic field and its gradients were generated by a superconducting magnet which can generate magnetic field up to 10 T and where gradients of the magnetic induction can reach up to 900 T²/m. A good agreement between experiments and numerical simulations was obtained in predicting the integral wall heat transfer over entire range of considered working parameters. Numerical simulations provided a detailed insights into changes of the local wall-heat transfer and long-term time averaged first and second moments for different strengths of the imposed magnetic induction.     

Experimental Investigation of Stability of a Hypersonic Boundary Layer on a Cone–Flare Model

Stability of a hypersonic flow in the regions of laminar separation of the boundary layer on a cone–flare model is experimentally studied for a Mach number M = 5.92. Development of natural disturbances and artificial wave packets in the boundary layer and separation region is examined. It is shown that highfrequency disturbances are predominantly amplified in the separation region; the most unstable waves are those propagating with an angle close to 60° to the freestream direction. It is found that separation and reattachment lines are generators of twodimensional disturbances.