Prediction of steady flows passing fixed cylinders using deep learning

Considerable attention has been given to deep-learning and machine-learning techniques in an effort to reduce the computational cost of computational fluid dynamics simulation. The present paper addresses the prediction of steady flows passing many fixed cylinders using a deep-learning model and investigates the accuracy of the predicted velocity field. The deep-learning model outputs the x- and y-components of the flow velocity field when the cylinder arrangement is input. The accuracy of the predicted velocity field is investigated, focusing on the velocity profile of the fluid flow and the fluid force acting on the cylinders. The present model accurately predicts the flow when the number of cylinders is equal to or close to that set in the training dataset. The extrapolation of the prediction to a smaller number of cylinders results in error, which can be interpreted as internal friction of the fluid. The results of the fluid force acting on the cylinders suggest that the present deep-learning model has good generalization performance for systems with a larger number of cylinders.

Steady plasma flows in a periodic non-symmetric domain

Steady plasma flows have been studied almost exclusively in systems with continuous symmetry or in open domains. In the absence of continuous symmetry, the lack of a conserved quantity makes the study of flows intrinsically challenging. In a toroidal domain, the requirement of double periodicity for physical quantities adds to the complications. In particular, the magnetohydrodynamics (MHD) model of plasma steady state with the flow in a non-symmetric toroidal domain allows the development of singularities when the rotational transform of the magnetic field is rational, much like the equilibrium MHD model. In this work, we show that steady flows can still be maintained provided the rotational transform is close to rational and the magnetic shear is weak. We extend the techniques developed in carrying out perturbation methods to all orders for static MHD equilibrium by Weitzner (Phys. Plasmas, vol. 21, 2014, p. 022515) to MHD equilibrium with flows. We construct perturbative MHD equilibrium in a doubly periodic domain with nearly parallel flows by systematically eliminating magnetic resonances order by order. We then utilize an additional symmetry of the flow problem, first discussed by Hameiri (J. Math. Phys., vol. 22, 1981, pp. 2080–2088, § III), to obtain a generalized Grad–Shafranov equation for a class of non-symmetric three-dimensional MHD equilibrium with flows both parallel and perpendicular to the magnetic field. For this class of flows, we can obtain non-symmetric generalizations of integrals of motion, such as Bernoulli's function and angular momentum. Finally, we obtain the generalized Hamada conditions, which are necessary to suppress singular currents in such a system when the magnetic field lines are closed. We do not attempt to address the question of neoclassical damping of flows.

The effects of non-Newtonian blood modeling and pulsatility on hemodynamics in the food and drug administration’s benchmark nozzle model

BACKGROUND: Computational fluid dynamics (CFD) is an important tool for predicting cardiovascular device performance. The FDA developed a benchmark nozzle model in which experimental and CFD data were compared, however, the studies were limited by steady flows and Newtonian models. OBJECTIVE: Newtonian and non-Newtonian blood models will be compared under steady and pulsatile flows to evaluate their influence on hemodynamics in the FDA nozzle. METHODS: CFD simulations were validated against the FDA data for steady flow with a Newtonian model. Further simulations were performed using Newtonian and non-Newtonian models under both steady and pulsatile flows. RESULTS: CFD results were within the experimental standard deviations at nearly all locations and Reynolds numbers. The model differences were most evident at Re = 500, in the recirculation regions, and during diastole. The non-Newtonian model predicted blunter upstream velocity profiles, higher velocities in the throat, and differences in the recirculation flow patterns. The non-Newtonian model also predicted a greater pressure drop at Re = 500 with minimal differences observed at higher Reynolds numbers. CONCLUSIONS: An improved modeling framework and validation procedure were used to further investigate hemodynamics in geometries relevant to cardiovascular devices and found that accounting for blood’s non-Newtonian and pulsatile behavior can lead to large differences in predictions in hemodynamic parameters.

Onset of absolutely unstable behaviour in the Stokes layer: a Floquet approach to the Briggs method

For steady flows, the Briggs (Electron-Stream Interaction with Plasmas. MIT Press, 1964) method is a well-established approach for classifying disturbances as either convectively or absolutely unstable. Here, the framework of the Briggs method is adapted to temporally periodic flows, with Floquet theory utilised to account for the time periodicity of the Stokes layer. As a consequence of the antiperiodicity of the flow, symmetry constraints are established that are used to describe the pointwise evolution of the disturbance, with the behaviour governed by harmonic and subharmonics modes. On coupling the symmetry constraints with a cusp-map analysis, multiple harmonic and subharmonic cusps are found for each Reynolds number of the flow. Therefore, linear disturbances experience subharmonic growth about fixed spatial locations. Moreover, the growth rate associated with the pointwise development of the disturbance matches the growth rate of the disturbance maximum. Thus, the onset of the Floquet instability (Blennerhassett & Bassom, J. Fluid Mech., vol. 464, 2002, pp. 393–410) coincides with the onset of absolutely unstable behaviour. Stability characteristics are consistent with the spatio-temporal disturbance development of the family-tree structure that has hitherto only been observed numerically via simulations of the linearised Navier–Stokes equations (Thomas et al., J. Fluid Mech., vol. 752, 2014, pp. 543–571; Ramage et al., Phys. Rev. Fluids, vol. 5, 2020, 103901).

On the lifespan of recirculating suspensions with pulsatile flow

A Lagrangian analysis is performed to measure the rate at which recirculating fluid is replaced (depleted) in pulsatile flows. Based on this approach, we then investigate how depletion is affected in dense suspensions. Experiments are conducted for pure liquid as well as suspensions with volume fractions of $\varPhi =5\,\%$ , 10 % and 20 %. Using Lagrangian tracking and pathline extension techniques, the depletion of the recirculation region is quantified via the trajectories of individual fluid parcels exiting the domain. Pulsatile flows with varying concentrations of hydrogel beads, up to a volume fraction of 20 %, are compared at mean Reynolds numbers of $Re=4800$ , 9600 and 14 400, while the Strouhal number ( $St=0.04$ , 0.08 and 0.15) and amplitude ratio ( $\lambda =0.25$ , 0.50 and 0.95) are systematically varied. A so-called ‘depletion efficiency’ is calculated for each test case, which is shown to increase with increasing Strouhal number and amplitude ratio. For most pulsatile cases, periodic vortex formation significantly increases depletion efficiency through enhanced entrainment of recirculating fluid. Conversely, low-amplitude pulsatile flows are dominated by Kelvin–Helmholtz instabilities, which do not penetrate into the recirculation region, and thus their depletion efficiency is markedly lower as a result. The efficiency trends and depletion mechanisms remain virtually unchanged between the pure liquid and each of the suspension concentrations under almost all flow conditions, which forms an unexpected conclusion. The only exception is for low-amplitude and steady flows, where increasing the suspension volume fraction is shown to suppress fluid transport across the shear layer, which in turn slows depletion and decreases the overall depletion efficiency.

The Momentum Conserving Scheme for Two-Layer Shallow Flows

This paper confronts the numerical simulation of steady flows of fluid layers through channels of varying bed and width. The fluid consists of two immiscible fluid layers with constant density, and it is assumed to be of a one-dimensional shallow flow. The governing equation is a coupled system of two-layer shallow water models. In this paper, we apply a direct extension of the momentum conserving scheme previously used for solving the one layer shallow water equations. Computations of various steady-state solutions are used to demonstrate the performance of the proposed numerical scheme. Under the influence of a given flow rate, the numerical steady interface is generated in a channel topography with a hump. The results obtained confirm the analytic steady interface of the two-layer rigid-lid model. Furthermore, the same scheme was used with an additional artificial damping to simulate the maximal exchange flow in channels of varying width. The numerical steady interface agreed well with the analytical steady solutions.

Kink Oscillations of Coronal Loops

Kink oscillations of coronal loops, i.e., standing kink waves, is one of the most studied dynamic phenomena in the solar corona. The oscillations are excited by impulsive energy releases, such as low coronal eruptions. Typical periods of the oscillations are from a few to several minutes, and are found to increase linearly with the increase in the major radius of the oscillating loops. It clearly demonstrates that kink oscillations are natural modes of the loops, and can be described as standing fast magnetoacoustic waves with the wavelength determined by the length of the loop. Kink oscillations are observed in two different regimes. In the rapidly decaying regime, the apparent displacement amplitude reaches several minor radii of the loop. The damping time which is about several oscillation periods decreases with the increase in the oscillation amplitude, suggesting a nonlinear nature of the damping. In the decayless regime, the amplitudes are smaller than a minor radius, and the driver is still debated. The review summarises major findings obtained during the last decade, and covers both observational and theoretical results. Observational results include creation and analysis of comprehensive catalogues of the oscillation events, and detection of kink oscillations with imaging and spectral instruments in the EUV and microwave bands. Theoretical results include various approaches to modelling in terms of the magnetohydrodynamic wave theory. Properties of kink oscillations are found to depend on parameters of the oscillating loop, such as the magnetic twist, stratification, steady flows, temperature variations and so on, which make kink oscillations a natural probe of these parameters by the method of magnetohydrodynamic seismology.