NEW METHOD FOR IMPROVED CALCULATIONS OF UNSTEADY COMPLEX FLOWS IN LARGE ARTERIES

Using an improved computational fluid dynamics (CFD) method developed for highly unsteady three-dimensional flows, numerical simulations for oscillating flow cycles and detailed unsteady simulations of the flow and forces on the aortic vessels at the iliac bifurcation, for both healthy and diseased patients, are analyzed. Improvements in computational efficiency and acceleration in convergence are achieved by calculating both an unsteady pressure gradient which is due to fluid acceleration and a good global pressure field correction based on mass flow for the pressure Poisson equation. Applications of the enhanced method to oscillatory flow in curved pipes yield an order of magnitude increase in speed and efficiency, thus allowing the study of more complex flow problems such as flow through the mammalian abdominal aorta at the iliac arteries bifurcation. To analyze the large forces which can exist on stent graft of patients with abdominal aortic aneurysm (AAA) disease, a complete derivation of the force equations is presented. The accelerated numerical algorithm and the force equations derived are used to calculate flow and forces for two individuals whose geometry is obtained from CT data and whose respective blood pressure measurements are obtained experimentally. Although the use of endovascular stent grafts in diseased patients can alter vessel geometries, the physical characteristics of stents are still very different when compared to native blood vessels of healthy subjects. The geometry for the AAA stent graph patient studied in this investigation induced flows that resulted in large forces that are primarily caused by the blood pressure. These forces are also directly related to the flow cross-sectional area and the angle of the iliac arteries relative to the main descending aorta. Furthermore, the fluid flow is significantly disturbed in the diseased patient with large flow recirculation and stagnant regions which are not present for healthy subjects.     

On Distributed $H^1$ Shape Gradient Flows in Optimal Shape Design of Stokes Flows: Convergence Analysis and Numerical Applications

We consider optimal shape design in Stokes flow using $H^1$ shape gradient flows based on the distributed Eulerian derivatives. MINI element is used for discretizations of Stokes equation and Galerkin finite element is used for discretizations of distributed and boundary $H^1$ shape gradient flows. Convergence analysis with a priori error estimates is provided under general and different regularity assumptions. We investigate the performances of shape gradient descent algorithms for energy dissipation minimization and obstacle flow. Numerical comparisons in 2D and 3D show that the distributed $H^1$ shape gradient flow is more accurate than the popular boundary type. The corresponding distributed shape gradient algorithm is more effective.     

非牛顿流体在偏心环空中轴向层流的摄动解

本文应用摄动方法研究了非牛顿流体在偏心环形空间中轴向层流的流动规律.文中以相对偏心度ε作为摄动参数,得到了其流场、极值速度、压力梯度等流动参数的一阶解.     

On some unsteady flows of a non-Newtonian fluid

Some properties of unsteady unidirectional flows of a fluid of second grade are considered for flows impulsively started from rest by the motion of a boundary or two boundaries or by sudden application of a pressure gradient. Flows considered are: unsteady flow over a plane wall, unsteady Couette flow, flow between two parallel plates suddenly set in motion with the same speed, flow due to one rigid boundary moved suddenly and one being free, unsteady Poiseuille flow and unsteady generalized Couette flow. The results obtained are compared with those of the exact solutions of the Navier–Stokes equations. It is found that the stress at time zero on the stationary boundary for the flows generated by impulsive motion of a boundary or two boundaries is finite for a fluid of second grade and infinite for a Newtonian fluid. Furthermore, it is shown that for unsteady Poiseuille flow the stress at time zero on the boundary is zero for a Newtonian fluid, but it is not zero for a fluid of second grade.     

Open-surface MHD flow over a curved wall in the 3-D thin-shear-layer approximation

3-D thin-shear-layer equations for flows of conducting fluids in a magnetic field have been derived in orthogonal body-oriented coordinates and then applied to the analysis of MHD open-surface flows over a curved wall. Unlike the classic boundary-layer-type equations, present ones permit information to be propagated upstream through the induced magnetic field. Another departure from the classic theory is that the normal momentum equation keeps the balance between the pressure gradient term, and those related to gravity, centrifugal forces, and Lorentz force. Thus, the normal pressure variations are allowed. The model describes basic 3-D effects due to the wall curvature and spatial variations of the applied magnetic field. As a particular case, equations for flows with rotational symmetry have been derived. Numerical calculations were performed for open-surface flows over a body of revolution under conditions relevant to a fusion reactor (Hartmann number is 8500). Some specific flow patterns, such as flow thickening and spiral-type flows, have been observed and discussed. A special attention has been paid to the analysis of the magnetic propulsion as a tool for the active flow control by applying an electric current. It has been shown that depending on the applied current, the axial pressure gradient can act as an adverse pressure gradient or propulsion force.     

An adjoint-based pressure correction scheme

Incompressible flows are solenoidal, i. e. the divergence of the flow field is zero. This is an algebraic constraint on the solution in time. The pressure has to be determined, so that the constraint is fulfilled. To calculate the pressure often a Poisson equation is derived, which is then solved by an iterative method. Instead of this the constraint is here formulated as an optimization problem. The objective functional is taken as the square of the norm of the divergence. A gradient based optimization is performed to calculate the pressure in every time step of the simulation. By this an alternative iterative scheme is derived. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Peristaltic transport of a Newtonian fluid in a vertical asymmetric channel with heat transfer and porous medium

The problem of peristaltic flow of a Newtonian fluid with heat transfer in a vertical asymmetric channel through porous medium is studied under long-wavelength and low-Reynolds number assumptions. The flow is examined in a wave frame of reference moving with the velocity of the wave. The channel asymmetry is produced by choosing the peristaltic wave train on the walls to have different amplitudes and phase. The analytical solution has been obtained in the form of temperature from which an axial velocity, stream function and pressure gradient have been derived. The effects of permeability parameter, Grashof number, heat source/sink parameter, phase difference, varying channel width and wave amplitudes on the pressure gradient, velocity, pressure drop, the phenomenon of trapping and shear stress are discussed numerically and explained graphically.     

Localized asymptotic solutions of the navier-stokes equations and laminar wakes in an incompressible fluid

Equations are derived to describe the far-field laminar wake behind a body in incompressible fluid flow with an arbitrary distribution of the free-stream (unperturbed flow) velocity. For certain classes of free-stream flows, analysis of these equations enables various processes in narrow wakes or jets to be described (the interaction of the longitudinal transverse velocity components in a jet, cause it to accelerate or decelerate and conservation of the energy of the wake by distortion of its trajectory regardless of viscous dissipation). In particular, conditions are obtained for the wake growth in spiral flows, analogous to the Rayleigh conditions for the instability of two-dimensionally radially symmetric flows relative to three-dimensional short-wave perturbations.     

Non-linear peristaltic flow of a non-Newtonian fluid under effect of a magnetic field in a planar channel

This paper is devoted to the study of peristaltic flow of a fourth grade fluid in a channel under the considerations of long wavelength and low-Reynolds number. The flow is examined in a wave frame of reference moving with velocity of the wave. The analytic solution has been obtained in the form of a stream function from which the axial velocity and axial pressure gradient have been derived. The results for the pressure rise and frictional force per wavelength have also been computed numerically. The computational results indicate that the pressure rise and frictional force per wavelength are increased in case of non-Newtonian fluid when compared with Newtonian fluid. Several graphs of physical interest are displayed and discussed.     

Numerical analysis of 3D vortical cavity flow

The vortical flows of an incompressible fluid in a rectangular three-dimensional container with a large spanwise aspect ratio driven by a moving solid lid are studied using a combined compact finite difference (CCD) scheme with high accuracy and high resolution. The study focuses on the change of the steady flow structures in the cavity with Reynolds numbers ranging from 100 to 850. The results of the flow in the cavity with a spanwise aspect ratio 6.55 show that several stable closed streamlines localized near the symmetric plane are found at Re ≥500, while a closed stable streamline is found near the side wall at Re ≤300. The change of the flow pattern present in this system affects the diffusion properties in the flow but seems to have no qualitative effect on the global flow properties which include energy dissipation in the cavity. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Heat transfer in a Walter’s liquid B fluid over an impermeable stretching sheet with non-uniform heat source/sink and elastic deformation

In this present article an analysis is carried out to study the boundary layer flow behavior and heat transfer characteristics in Walter’s liquid B fluid flow. The stretching sheet is assumed to be impermeable, the effects of viscous dissipation, non-uniform heat source/sink in the presence and in the absence of elastic deformation (which was escaped from attention of researchers while formulating the viscoelastic boundary layer flow problems)on heat transfer are addressed. The basic boundary layer equations for momentum and heat transfer, which are non-linear partial differential equations, are converted into non-linear ordinary differential equations by means of similarity transformation. Analytical solutions are obtained for the resulting boundary value problems. The effects of viscous dissipation, Prandtl number, Eckert number and non-uniform heat source/sink on heat transfer (in the presence and in the absence of elastic deformation) are shown in several plots and discussed. Analytical expressions for the wall frictional drag coefficient, non-dimensional wall temperature gradient and non-dimensional wall temperature are obtained and are tabulated for various values of the governing parameters. The present study reveals that, the presence of work done by deformation in the energy equation yields an augment in the fluid’s temperature.     

Two-layer flows of generalized immiscible second grade fluids in a rectangular channel

Unsteady one-dimensional flows of two incompressible and immiscible generalized second grade fluids in a rectangular channel are studied. A constant pressure gradient acts in the flow direction, while the channel walls have oscillating translational motions in their planes. The generalization considered in this paper consists into a mathematical model based on constitutive equations of second grade fluid with Caputo time-fractional derivative in which the history of the shear stress influences the velocity gradient. The velocity and shear stress fields in the Laplace transform domain are obtained. Numerical solutions for the real velocity and shear stress have been found by employing the Stehfest numerical algorithm for the inverse Laplace transform. The influence of the fractional parameters on the velocity and shear stress has been studied by numerical simulations and graphical illustrations. It is found that the memory effects are significant only for small values of the time t.     

The linear stability of inviscid shear flow of a vibrationally excited diatomic gas

The problem of the linear stability of plane-parallel shear flows of a vibrationally excited compressible diatomic gas is investigated using a two-temperature gas dynamics model. The necessary and sufficient conditions for stability of the flows considered are obtained using the energy integrals of the corresponding linearized system for the perturbations. It is proved that thermal relaxation produces an additional dissipation factor, which enhances the flow stability. A region of eigenvalues of unstable perturbations is distinguished in the upper complex half-plane. Numerical calculations of the eigenvalues and eigenfunctions of the unstable inviscid modes are carried out. The dependence on the Mach number of the carrier stream, the vibrational relaxation time τ and the degree of non-equilibrium of the vibrational mode is analysed. The most unstable modes with maximum growth rate are obtained. It is shown that in the limit there is a continuous transition to well-known results for an ideal fluid as the Mach number and τ approach zero and for an ideal gas when τ → 0.     

Some analytical solutions for second grade fluid flows for cylindrical geometries

This paper deals with some unsteady flow problems of a second grade fluid. The flows are generated by the sudden application of a constant pressure gradient or by the impulsive motion of a boundary. The velocities of the flows are described by the partial differential equations. Exact analytic solutions of these differential equations are obtained. The well known solutions for a Navier–Stokes fluid in the hydrodynamic case appear as the limiting cases of our solutions.     

On the Relation between Enhanced Dissipation Timescales and Mixing Rates

We study diffusion and mixing in different linear fluid dynamics models, mainly related to incompressible flows. In this setting, mixing is a purely advective effect that causes a transfer of energy to high frequencies. When diffusion is present, mixing enhances the dissipative forces. This phenomenon is referred to as enhanced dissipation, namely the identification of a timescale faster than the purely diffusive one. We establish a precise connection between quantitative mixing rates in terms of decay of negative Sobolev norms and enhanced dissipation timescales. The proofs are based on a contradiction argument that takes advantage of the cascading mechanism due to mixing, an estimate of the distance between the inviscid and viscous dynamics, and an optimization step in the frequency cutoff. Thanks to the generality and robustness of our approach, we are able to apply our abstract results to a number of problems. For instance, we prove that contact Anosov flows obey logarithmically fast dissipation timescales. To the best of our knowledge, this is the first example of a flow that induces an enhanced dissipation timescale faster than polynomial. Other applications include passive scalar evolution in both planar and radial settings and fractional diffusion. © 2019 Wiley Periodicals, Inc.     

Unilateral gradient flow of the Ambrosio–Tortorelli functional by minimizing movements

Motivated by models of fracture mechanics, this paper is devoted to the analysis of a unilateral gradient flow of the Ambrosio–Tortorelli functional, where unilaterality comes from an irreversibility constraint on the fracture density. Solutions of such evolution are constructed by means of an implicit Euler scheme. An asymptotic analysis in the Mumford–Shah regime is then carried out. It shows the convergence towards a generalized heat equation outside a time increasing crack set. In the spirit of gradient flows in metric spaces, a notion of curve of maximal unilateral slope is also investigated, and analogies with the unilateral slope of the Mumford–Shah functional are also discussed.     

Gauging magnetorotational instability

Previously (Z. Angew. Math. Phys. 57:615–622, 2006), we examined the axisymmetric stability of viscous resistive magnetized Couette flow with emphasis on flows that would be hydrodynamically stable according to Rayleigh’s criterion: opposing gradients of angular velocity and specific angular momentum. A uniform axial magnetic field permeates the fluid. In this regime, magnetorotational instability (MRI) may occur. It was proved that MRI is suppressed, in fact no instability at all occurs, with insulating boundary conditions, when a term multipling the magnetic Prandtl number is neglected. Likewise, in the current work, including this term, when the magnetic resistivity is sufficiently large, MRI is suppressed. This shows conclusively that small magnetic dissipation is a feature of this instability for all magnetic Prandtl numbers. A criterion is provided for the onset of MRI.     

Gauging magnetorotational instability

Previously (Z. Angew. Math. Phys. 57:615–622, 2006), we examined the axisymmetric stability of viscous resistive magnetized Couette flow with emphasis on flows that would be hydrodynamically stable according to Rayleigh’s criterion: opposing gradients of angular velocity and specific angular momentum. A uniform axial magnetic field permeates the fluid. In this regime, magnetorotational instability (MRI) may occur. It was proved that MRI is suppressed, in fact no instability at all occurs, with insulating boundary conditions, when a term multipling the magnetic Prandtl number is neglected. Likewise, in the current work, including this term, when the magnetic resistivity is sufficiently large, MRI is suppressed. This shows conclusively that small magnetic dissipation is a feature of this instability for all magnetic Prandtl numbers. A criterion is provided for the onset of MRI.