Rayleigh–Lamb Waves in a Compound Wave Guide. Reflection from and Transmission Through the Vertical Interface between Media with Different Impedances

The method of superposition is used to study the first normal wave reflecting from and transmitting through the interface in a compound waveguide consisting of two rigidly joined elastic half-strips with equal width and different mechanical properties. We study how the impedances of the contacting media influence the transformation of the energy of the incident wave to those of the reflected and transmitted waves. Two cases are considered — propagating waves of higher orders appear in the reflected wavefield earlier than in the transmitted wavefield and propagating waves of higher orders appear in the transmitted wavefield earlier than in the reflected wave field. For both cases, the impedances vary so that the incident wave can propagate in both more rigid and softer media. It is shown that by increasing the impedances of the contacting media, the interface can be made more transparent     

Reflection and transmission of regular waves at a surface-pitching slotted barrier

The interactions between regular surface waves and a surface-pitching slotted barrier are investigated both analytically and experimentally.A quasi-linear theory is developed using the eigenfunction expansion method.The energy dissipation within the barriers is modeled by a quadratic friction factor,and an equivalent linear dissipation coefficient,which is depth-varying,wave-height dependent,is introduced to linearize the matching condition at the surface-pitching barrier.By comparing the theoretical results with laboratory experiments,it is shown that the present method can satisfactorily predict the variation of the reflection and transmission coefficients with wave height.     

Reflection and Transmission of Plane Waves at a Water–Porous Sediment Interface with a Double-Porosity Substrate

This paper investigates the wave propagation at the interface between the ocean and the ocean floor. The ocean floor is assumed to be composed of covered porous sediment with an underlying double-porosity substrate. For this purpose, plane wave reflection and transmission in the coupled water–porous sediment–double-porosity substrate system are analytically solved in terms of displacement potentials. Using numerical examples, the effects of the material properties of the underlying double-porosity substrate on the reflection coefficients are discussed in detail. Variations in pore and fracture fluid, fracture volume fraction, and permeability coefficients are considered. In addition, two cases of boundary conditions at the porous sediment–double-porosity substrate interface, i.e., sealed-pore boundary and open-pore boundary, are compared in the numerical calculations. Results show that material property variations in the double-porosity substrate may significantly affect the reflected wave in the overlying water if the sandwiched sediment depth is less than the critical value.     

Reflection of an Oblique Shock Wave in a Reacting Gas with a Finite Relaxation–Zone Length

Reflection of an oblique shock wave in a reacting gas with a finite length of the chemical–reaction zone is studied. Shock polars for an arbitrary heat release behind the oblique shock wave are constructed. Transition criteria from regular to Mach reflection and back are obtained. It is shown that transition criteria are significantly changed if the reaction–zone length is taken into account.     

Numerical Modeling of Unsteady Radiative–Convective Heat Transfer on a Flat–Plate Boundary Layer in a Selectively Emitting and Scattering Medium

A conjugation problem for radiative–convective heat transfer in a turbulent flow of a high–temperature gas—particle medium around a thermally thin plate is considered. The plate experiences intense heating from an outside source that emits radiation in a restricted spectral range. Unsteady temperature fields and heat–flux distributions along the plate are calculated. The results permit prediction of the effect of the type and concentration of particles on the dynamics of the thermal state of both the medium in the boundary layer and the plate itself under conditions of its outside heating by a high–temperature source of radiation.     

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.     

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 simulation of the fluid–structure interaction in a tube array under cross flow at moderate and high Reynolds number

The unsteady loads in a tube bundle are studied at moderate and high Reynolds number by means of URANS and hybrid (DDES) modelling. The onset of fluid-elastic instability is analysed for different structural parameters, Scruton number and reduced velocity. The simulations have been carried out with the code NSMB (Navier–Stokes Multi Block) by using turbulence modelling methods URANS and DDES (Delayed Detached Eddy Simulation). The CEA-DIVA configuration is considered for the cylinders array for an inter-tube Reynolds number 60 000. The study is carried out for a configuration of (4×5) cylinders in static conditions as well as for the vertical free motion of one of the central cylinders in one DOF (Degree Of Freedom).The inter-tube Reynolds number is 60 000. It is found that this cylinder spontaneously displays an oscillatory motion which first corresponds to Vortex Induced Vibration (VIV), associated to a lock-in mechanism for low values of the reduced velocity and secondly develops Movement Induced Vibration, MIV, for higher values of the reduced velocity. The variation of the cylinder׳s oscillations frequency, of the unsteady loads and the structure׳s displacement are studied as a function of the reduced velocity for low and high values of the Scruton number. The increase of the phase-lag between the forces and the displacement is predicted and discussed for different Scruton number values and reduced velocities.     

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.     

Wave features of a hyperbolic reaction–diffusion model for Chemotaxis

The Extended Thermodynamic theory is used to derive a hyperbolic reaction–diffusion model for Chemotaxis. Linear stability analysis is performed to study the nature of the equilibrium states against uniform and nonuniform perturbations. A particular emphasis is given to the occurrence of the Turing bifurcation. The existence of traveling wave solutions connecting the two steady states is investigated and the governing equations are numerically integrated to validate the analytical results. The propagation of plane harmonic waves is analyzed and the stability regions in terms of the model parameters are shown. The frequency dependence of the phase velocity and of the attenuation is also illustrated. Finally, in order to have a measure of the non linear stability, the propagation of acceleration waves is studied, the wave amplitude is derived and the critical time is evaluated.     

Existence of Strong Solutions for a System Coupling the Navier–Stokes Equations and a Damped Wave Equation

We consider a fluid–structure interaction problem coupling the Navier–Stokes equations with a damped wave equation which describes the displacement of a part of the boundary of the fluid domain. The system is considered first in the two-dimensional setting and in a second part it is adapted to the three-dimensional setting.     

Traveling Wave Solutions for a Class of Predator–Prey Systems

We use a shooting method to show the existence of traveling wave fronts and to obtain an explicit expression of minimum wave speed for a class of diffusive predator?Cprey systems. The existence of traveling wave fronts indicates the existence of a transition zone from a boundary equilibrium to a co-existence steady state and the minimum wave speed measures the asymptotic speed of population spread in some sense. Our approach is a significant improvement of techniques introduced by Dunbar. The advantage of our method is that it does not need the notion of Wazewski??s set and LaSalle??s invariance principle used in Dunbar??s approach. In our approach, we convert the equations for traveling wave solutions to a system of first order equations by a ??non-traditional transformation??. With this converted new system, we are able to construct a Liapunov function, which gives an immediate implication of the boundedness and convergence of the relevant class of heteroclinic orbits. Our method provides a more efficient way to study the existence of traveling wave solutions for general predator?Cprey systems.     

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.     

One-Dimensional Interaction of a Cylindrical Unloading Wave with a Moving Elastic–Plastic Boundary

The one-dimensional dynamic problem of the theory of large elastic–plastic deformations is considered for the interaction of an unloading wave with an elastic–plastic boundary. It is shown that before the occurrence of the unloading wave, the increasing pressure gradient leads to quasistatic deformation of the elasti©viscoplastic material filling the round tube, which is retained in the tube due to friction on its wall, resulting in the formation of near-wall viscoplastic flow and an elastic core. The unloading wave is initiated at the moment of the onset of slippage of the material along the inner wall of the tube. Calculations were conducted using the ray method of constructing approximate solutions behind strong discontinuity surfaces, and ray expansions of the solutions behind the cylindrical surfaces of discontinuities were obtained.     

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.