Steady Flow Simulations Research Articles

Stratification of a population of intracranial aneurysms using blood flow metrics

Indices of the intra-aneurysm hemodynamic environment have been proposed as potentially indicative of their longitudinal outcome. To be useful, the indices need to be used to stratify large study populations and tested against known outcomes. The first objective was to compile the diverse hemodynamic indices reported in the literature. Furthermore, as morphology is often the only patient-specific information available in large population studies, the second objective was to assess how the ranking of aneurysms in a population is affected by the use of steady flow simulation as an approximation to pulsatile flow simulation, even though the former is clearly non-physiological. Sixteen indices of aneurysmal hemodynamics reported in the literature were compiled and refined where needed. It was noted that, in the literature, these global indices of flow were always time-averaged over the cardiac cycle. Steady and pulsatile flow simulations were performed on a population of 198 patient-specific and 30 idealised aneurysm models. All proposed hemodynamic indices were estimated and compared between the two simulations. It was found that steady and pulsatile flow simulations had a strong linear dependence (r ≥ 0.99 for 14 indices; r ≥ 0.97 for 2 others) and rank the aneurysms in an almost identical fashion (ρ ≥ 0.99 for 14 indices; ρ ≥ 0.96 for other 2). When geometry is the only measured piece of information available, stratification of aneurysms based on hemodynamic indices reduces to being a physically grounded substitute for stratification of aneurysms based on morphology. Under such circumstances, steady flow simulations may be just as effective as pulsatile flow simulation for estimating most key indices currently reported in the literature.

Pulsatile flow analysis in a bend square channel by modal decomposition

The article analyzes the results of numerical simulations of pulsatile flow in a square channel with a 90° bend using Proper Orthogonal Decomposition (POD). The simulations of laminar steady and unsteady flow obtained by ANSYS Fluent commercial system were used as input for modal decomposition. Several regimes for different values of individual parameters (stationary flow component, ) were simulated. Not only energy-related modes but also modes significant from the viewpoint of energy dissipation. Other evaluated modes are the modes corresponding to the pulsatile flow frequency and the modes reflecting the flow separation at the bend.

Coupled Fluid, Thermal and Structural Analysis of Nozzle in Solid Rocket Motor

A coupled fluid, thermal and structure analysis model was established for solid rocket motor nozzle. Based on a steady flow simulation and a transient thermal simulation, a transient structure analysis simulation was executed. The model took a frictional coefficient 0.3 into account to simulate a rub action between the parts of the nozzle. Result showed that with a frictional coefficient, Von Mises stress could be released which gives a constructive recommendation for a structure simulation of a nozzle. Also, the model was more close to the actual. A more dependable approach was introduced and it could contribute to designing department.

Approximating hemodynamics of cerebral aneurysms with steady flow simulations

Computational fluid dynamics (CFD) simulations can be employed to gain a better understanding of hemodynamics in cerebral aneurysms and improve diagnosis and treatment. However, introduction of CFD techniques into clinical practice would require faster simulation times. The aim of this study was to evaluate the use of computationally inexpensive steady flow simulations to approximate the aneurysm's wall shear stress (WSS) field. Two experiments were conducted. Experiment 1 compared for two cases the time-averaged (TA), peak systole (PS) and end diastole (ED) WSS field between steady and pulsatile flow simulations. The flow rate waveform imposed at the inlet was varied to account for variations in heart rate, pulsatility index, and TA flow rate. Consistently across all flow rate waveforms, steady flow simulations accurately approximated the TA, but not the PS and ED, WSS field. Following up on experiment 1, experiment 2 tested the result for the TA WSS field in a larger population of 20 cases covering a wide range of aneurysm volumes and shapes. Steady flow simulations approximated the space-averaged WSS with a mean error of 4.3%. WSS fields were locally compared by calculating the absolute error per node of the surface mesh. The coefficient of variation of the root-mean-square error over these nodes was on average 7.1%. In conclusion, steady flow simulations can accurately approximate the TA WSS field of an aneurysm. The fast computation time of 6min per simulation (on 64 processors) could help facilitate the introduction of CFD into clinical practice.

Computational Modeling of Non-Newtonian Blood Flow Through Stenosed Arteries in the Presence of Magnetic Field

Steady flow simulations of blood flow in an axisymmetric stenosed artery, subjected to a static magnetic field, are performed to investigate the influence of artery size, magnetic field strength, and non-Newtonian behavior on artery wall shear stress and pressure drop in the stenosed section. It is found that wall shear stress and pressure drop increase by decreasing artery size, assuming non-Newtonian fluid, and increasing magnetic field strength. In the computations, the shear thinning behavior of blood is accounted for by the Carreau-Yasuda model. Computational results are compared and found to be inline with available experimental data.

Implicit high-order discontinuous Galerkin method with HWENO type limiters for steady viscous flow simulations

Two types of implicit algorithms have been improved for high order discontinuous Galerkin (DG) method to solve compressible Navier-Stokes (NS) equations on triangular grids. A block lower-upper symmetric Gauss-Seidel (BLU-SGS) approach is implemented as a nonlinear iterative scheme. And a modified LU-SGS (LLU-SGS) approach is suggested to reduce the memory requirements while retain the good convergence performance of the original LU-SGS approach. Both implicit schemes have the significant advantage that only the diagonal block matrix is stored. The resulting implicit high-order DG methods are applied, in combination with Hermite weighted essentially non-oscillatory (HWENO) limiters, to solve viscous flow problems. Numerical results demonstrate that the present implicit methods are able to achieve significant efficiency improvements over explicit counterparts and for viscous flows with shocks, and the HWENO limiters can be used to achieve the desired essentially non-oscillatory shock transition and the designed high-order accuracy simultaneously.

OpenFOAM의 비압축성 유동 해석정밀도 평가

The goal of the research is to evaluate the open source code of OpenFOAM for the use of nuclear plant flow simulation objectively. Of the various incompressible flow solvers, simpleFoam, pimpelFoam are then tested under three validated cases (backward facing step, flow over circular cylinder and turbulent round jet flow). For the evaluation of steady state incompressible laminar flow simulation, low reynolds number of backward facing step flow was solved by simpleFoam. The resultant of the reattached lengths turned out to be similar with the other experimental and simulation results. For transient flow simulation, flow over circular cylinder and turbulent round jet flow were solved by pimpleFoam. The simulation accuracy was evaluated by comparing the resultant flow patterns with the description of the characteristics of the flow over the circular cylinder. The quantitative accuracy was evaluated for no more than 85% by comparing it to the decaying constants of the turbulent round jet velocity.

Flow prediction in cerebral aneurysms based on geometry reconstruction from 3D rotational angiography

We present an immersed boundary (IB) method for the simulation of steady blood flow inside a realistic cerebral aneurysm. We reconstruct a segment of the cerebrovascular system that contains an aneurysm, by using medical images obtained with three dimensional rotational angiography (3DRA). The main focus is on evaluating the sensitivity of flow predictions to the various steps of the vascular reconstruction process. Starting from the raw medical data, we analyze the fluid-mechanical consequences of the steps needed to generate the IB masking function for our simulations. We illustrate the IB method by applying it to a realistic aneurysm and investigate the role of (i) numerical resolution of the geometry; (ii) the selection of the specific vascular segment used in the simulations; and (iii) the influence of the smoothness of the periodic vessel extension to complete the computational model. Because of an unavoidable degree of uncertainty in the medical images, the geometry of the vessels and the aneurysm can be reconstructed only approximately. We also incorporate these slight uncertainties in the masking function by introducing inner and outer 'bounding' geometries and analyze the sensitivity of the flow predictions to these variations in the masking function. The numerical solutions computed in the inner and outer bounding geometries provide practical upper and lower bounds for basic flow properties, thus quantifying the reliability of the numerical solution, subject to uncertainties in the geometry of the flow domain.

Design of high performance badminton shuttlecocks: virtual and rapid prototyping approach

In this paper, a framework utilising virtual and rapid prototyping to aid design of a shuttlecock skirt was proposed and implemented. A standard template for flow simulation on ANSYS® CFX was developed to predict the aerodynamic performance of the virtual prototypes, with focus on reduced simulation time and complexity. This forms the basis of the computer aided design model – computational fluid dynamics – performance evaluation loop in the proposed framework. Steady state flow simulation was conducted for three simple computer aided design models of shuttlecocks, with different gap dimensions on the skirt. Wind tunnel drag study of the models reproduced through rapid prototyping validated the simulation result. The proposed methodology can also be applied to the existing design process for integration of virtual and rapid prototyping to reduce the burden of design iteration.

Application of entropy production theory to hydro-turbine hydraulic analysis

The understanding of hydraulic behavior in the hydro turbine requires the detailed study of fluid flow in the turbine. Previous methods of analyzing the numerical simulation results on the fluid machinery are short of intuitiveness on energy dissipation. In this paper, the energy dissipation was analyzed based on the entropy production theory. 3-D steady flow simulations and entropy production calculations of the reduced hydro turbine were carried out. The results indicated that the entropy production theory was suitable for evaluating the performance of the hydro turbine. The energy dissipation in the guide vanes area weighted nearly 25% of the whole flow passage, and mainly happened at the head and tail areas of the vanes. However, more than half the energy dissipation occurred in the runner, mostly at the leading edge of runner blade and the trailing edge of runner blade. Meanwhile, close to 20% of the energy dissipation occurred in the elbow. And it can be concluded that the method of entropy production analysis has the advantages of determining the quantity of energy dissipation and where the dissipation happens.

Hamiltonian thermostats fail to promote heat flow

Hamiltonian mechanics can be used to constrain temperature simultaneously with energy. We illustrate the interesting situations that develop when two different temperatures are imposed within a composite Hamiltonian system. The model systems we treat are ϕ4 chains, with quartic tethers and quadratic nearest-neighbor Hooke’s-law interactions. This model is known to satisfy Fourier’s law. Our prototypical problem sandwiches a Newtonian subsystem between hot and cold Hamiltonian reservoir regions. We have characterized four different Hamiltonian reservoir types. There is no tendency for any of these two-temperature Hamiltonian simulations to transfer heat from the hot to the cold degrees of freedom. Evidently steady heat flow simulations require energy sources and sinks, and are therefore incompatible with Hamiltonian mechanics.

Fluid Mechanics Analysis and Effect of Splitter Blades on the Characteristics of High-Speed Centrifugal Pump

Due to its single-stage high head, compact structure, convenient maintenance, low cost and high reliability, the high-speed centrifugal pump has been widely used in refining, petrochemical, chemical industry, and aerospace engineering fields. However, as the speed increases, the increasingly high demand for stability is showing and cavitation characteristics of the high-speed centrifugal pump is emerging. Reducing the mechanical vibration, hydraulic losses and improving the pressure distribution of the high-speed pump have become one of the main research fields for high-speed pumps. Three types of splitter blades were designed in this paper. Numerical simulations of the 3D steady flow in the high-speed pump have been carried out by using the software ANSYS CFX 12.0 and the comparative analysis of the velocity, pressure distributions and fluid mechanics in the high-speed pump under three different flow conditions were also carried out，Respectively. The effects of splitter blades on the high-speed pump performance were studied by analyzing both its inner flow field and hydraulic performance. The calculated H-Q curve of the high-speed pump was compared with the test results, and they agreed well with each other. The relative error of head at design conditions was lower than that at the non-design conditions.

Simulation of Back Pressure Effect on Behaviour of Convergent Divergent Nozzle

In this research a simulation of steady flow of a gas through a convergent divergent nozzle which has a varying cross sectional area will be considered. The nature of the flow can be explained by considering how the flow and its characteristics in the nozzle changes as nthe back pressure Pb is decreased.The characteristics of gas flow i.e.(Mach number, static pressure, density, velocity magnitude and static temperature) distributions for the convergent divergent nozzle are implemented by using the ANSYS Fluent 12.1 software to solve the quasi-one dimensional nozzle flow.The reductions in the back pressure cannot affect conditions upstream of the throat. The nozzle is, therefore, choked. The shock wave increases the pressure, density and temperature and reduces the velocity and Mach number to a subsonic value, and as back pressure is further reduced to a certain value, the extent of the supersonic flow region increases, the shock wave moving further down the divergent portion of the nozzle towards the exit plane.

A low-dimensional tool for predicting force decomposition coefficients for varying inflow conditions

We develop a low-dimensional tool to predict the effects of unsteadiness in the inflow on force coefficients acting on a circular cylinder using proper orthogonal decomposition (POD) modes from steady flow simulations. The approach is based on combining POD and linear stochastic estimator (LSE) techniques. We use POD to derive a reduced-order model (ROM) to reconstruct the velocity field. To overcome the difficulty of developing a ROM using Poisson’s equation, we relate the pressure field to the velocity field through a mapping function based on LSE. The use of this approach to derive force decomposition coefficients (FDCs) under unsteady mean flow from basis functions of the steady flow is illustrated. For both steady and unsteady cases, the final outcome is a representation of the lift and drag coefficients in terms of velocity and pressure temporal coefficients. Such a representation could serve as the basis for implementing control strategies or conducting uncertainty quantification.

J013021 テンソル分解による解析結果データのデータ圧縮*([J013-02]流体情報学とデータ解析(2))

With the increase in the size of simulation results database, it has become a big challenge to transfer large data between different locations. In this study, we have developed a data compression method based on tensor compression. Voxel simulation results data are represented as third-order tensors and a tensor decomposition algorithm based on HOOI is applied to the tensors. After tensor decompression, the original tensor is decomposed into a core tensor and a series of basis matrices, whose summed size is considerably smaller than that of the original tensor. As a result, a compression ratio of over 60:1 is achieved for steady flow simulation and the error is below 5%. A compression ratio of over 70:1 is achieved for unsteady flow simulation and the error is below 5%. We have confirmed that no significant information is lost during the data compression process.

English

In this work, the numerical simulation of steady heat transfer and fluid flow over a bank of flat tubes in staggered configurations for determining the constant surface temperature is presented. The results are attained using the finite volume method (FVM) and body fitted coordinates (BFC) technique. Transverse ratios (ST/Ds) of the pitch to small diameter of 3.0, 4.0 and 6.0 are also considered. The Reynolds numbers used are 10, 20, 60, 80 and 100, and the Prandtl number is taken as 0.7. The isothermal line, streamline and average Nusselt number were analyzed in this paper. It was found that the strength of the heat transfer between the surface of the tubes and the air flow increases with increasing Reynolds number and increasing pitch-to-diameter ratios. Also, the effect of the Reynolds number clear for the isothermal line, streamline and the average Nusselt number.

Frictional dependence of shallow-granular flows from discrete particle simulations

A shallow-layer model for granular flows down inclines is completed with a closure relation for the macroscopic bed friction obtained from micro-scale, discrete particle simulations of steady flows over geometrically rough bases with contact friction. Microscopic friction can be different between bulk particles and with particles at the base, where the latter is systematically varied. When extending the known friction closure relation to be a function of both bulk flow and bed properties, surprisingly, we find that the macroscopic bed friction is only weakly dependent on the contact friction of the bed particles and is predominantly determined by the properties of the flowing particles. Implications for constitutive modelling and possible experiments to better understand the bulk rheology are discussed.

Significance of plaque morphology in modifying flow characteristics in a diseased coronary artery: Numerical simulation using plaque measurements from intravascular ultrasound imaging

The paper deals with numerical investigation of the effect of plaque morphology on the flow characteristics in a diseased coronary artery using realistic plaque morphology. The morphological information of the lumen and the plaque is obtained from intravascular ultrasound imaging measurements of 42 patients performed at Cleveland Clinic Foundation, Ohio. For this data, study of Bhaganagar et al. (2010) [1] has revealed the stenosis for 42 patients can be categorized into four types – type I (peak-valley), type II (ascending), type III (descending), and type IV (diffuse). The aim of the present study is to isolate the effect of shape of the stenosis on the flow characteristics for a given degree of the stenosis. In this study, we conduct fluid dynamic simulations for the four stenosis types (type I–IV) and analyze the differences in the flow characteristics between these types. Finely refined tetrahedral mesh for the 3-D solid model of the artery with plaques has been generated. The 3-D steady flow simulations were performed using the turbulence (k–ε) model in a finite volume based computational fluid dynamics solver. The axial velocity, the radial velocity, turbulence kinetic energy and wall shear stress profiles of the plaque have been analyzed. From the axial and radial velocity profiles results the differences in the velocity patterns are significantly visible at proximal as well as distal to the throat, region of maximum stenosis. Turbulent kinetic energy and wall shear stress profiles have revealed significant differences in the vicinity of the plaque. Additional unsteady flow simulations have been performed to validate the hypothesis of the significance of plaque morphology in flow alterations in diseased coronary artery. The results revealed the importance of accounting for plaque morphology in addition to plaque height to accurately characterize the turbulent flow in a diseased coronary artery.

Investigation of skin porosity damping effects on free stream disturbance induced unsteady wing loads

AbstractNumerical simulations were performed to analyse the possibility of damping abrupt incoming free stream disturbances upon a porous aerofoil using an unsteady Reynolds-averaged Navier-Stokes (URANS) model. To mimic the turbulence disturbance levels that are typically encountered in the atmosphere, two flow configurations were considered. In the first configuration, the unsteadiness of the flow was created with vortices shed from a circular cylinder installed ahead of a WTEA-TE1 aerofoil. The continuous von Kármán shedding vortices contained within the cylinder wake were convected downstream and projected upon the aerofoil. In the second configuration, an instantaneous pair of discrete vortices was created by a rotational snapping of a flat plate, installed upstream of the aerofoil. Solid and porous aerofoil configurations, with porosity settings of 11 and 22%, were applied on 50% of the chord of the aerofoil starting from the leading edge. Both steady and unsteady flow simulations were performed to assess the performance of the porosity under steady and unsteady effects. The steady state flow simulations revealed a noticeable reduction in the aerofoil lift coefficient for the porous aerofoil. For unsteady solutions with a continuous or distinct series of vortices interacting with the aerofoil, the porosity showed insignificant damping of the lift coefficient amplitude. The porosity values investigated in the current exercise had indiscernible effect upon the unsteady lift-load alleviations caused by free stream disturbances.

고압터빈 익렬 주위 유동해석에서 난류모델의 영향 평가

Steady flow simulations through a high pressure turbine guide vanes were carried out. The main objective of the present work is to study the performance of turbulence models on the steady flow prediction from aerodynamic and aerothermal points of view. Three turbulence models were compared, namely SST, k-<TEX>${\omega}$</TEX> and <TEX>${\omega}$</TEX>-Reynolds stress models. The laminar results were also compared. The comparison was done with emphasis on the isentropic Mach number and heat transfer coefficient along the blade, and total pressure loss in the wake region. The calculated isentropic Mach number showed reasonable agreement with experimental data along the blade surface for all three turbulent models. For the total pressure loss in the wake region, <TEX>${\omega}$</TEX>-Reynolds stress model showed the best agreement with the experimental data. However, unless using an appropriate transition model, the heat transfer coefficients of all three turbulent models showed poor agreement with experimental data.