Efficient implementation of the proper outlet flow conditions in blood flow simulations through asymmetric arterial bifurcations

We present an efficient implementation of the proper (in vivo) outlet boundary conditions in detailed, three‐dimensional (3D) and time‐periodic simulations of blood flow through arteries. This is achieved through the intermediate use of an approximate ‘simulant’ model of the outlet pressure/flow relationship corresponding to the full 3D and time‐dependent numerical simulation. This model allows us to efficiently couple the 3D outlet pressure/flow conditions to the equivalent relations due to the downstream arterial network, as obtained from a one‐dimensional approximate model in the form of Fourier frequency impedance coefficients. An adjustable time‐periodic function correction term in the simulant model requires input from the full 3D model that has to run iteratively until convergence. The advantage of the proposed numerical scheme is that it decouples the upstream detailed simulation from the downstream approximate network model offering exceptional versatility. This approach is demonstrated here in a series of detailed 3D simulations of blood flow, performed using the commercial software FLUENT?, through an asymmetric arterial bifurcation. Two cases are considered: first a healthy system patterned after the left main coronary arterial bifurcation, and second a diseased case where an occlusion has developed in one of the daughter vessels, resulting in strengthening the asymmetry of the bifurcation. Rapid convergence of the iterative process was achieved in both cases. Subtle changes occur in the shear patterns of the daughter vessels, whereas the flow distribution is quite different. In the presence of a stenosis additional regions of low shear develop due to inertial effects. Copyright © 2010 John Wiley & Sons, Ltd.     

Simulation of oscillatory flow in an aortic bifurcation using FVM and FEM: A comparative study of implementation strategies

Cardiovascular diseases are one of the major causes of long‐term morbidity and mortality in human beings. The nearly epidemic increase in prevalence of such diseases poses a serious threat to public health and calls for efficient methods of diagnosis and treatment. Non‐invasive diagnostic procedures such as MRI are often used in this context; however, these are limited in terms of spatial and temporal resolution and do not provide information on time‐dependent pressures and wall shear stresses—key quantities considered to be partially responsible for the formation and development of related pathologies. The present study is concerned with the numerical simulation of oscillatory flow through the abdominal aortic bifurcation. Computational fluid dynamics simulation of oscillatory flow in a branched geometry at high Reynolds numbers poses considerable challenges. The present study reports a detailed comparison of simulations performed with a finite volume and a finite element method, two approaches with significant differences in their discretization strategy, treatment of boundary conditions and other numerical aspects. Both solvers were parallelized, using loop parallelization of the BiCGStab linear solver for the finite volume and domain decomposition based on the Schur complement method for the finite element technique. The experience gained with these two approaches for the solution of flow in a bifurcation forms the focus of this study. Although similar results were obtained for both methods, the computation time required for convergence was found to be significantly smaller for the finite element approach. Copyright © 2010 John Wiley & Sons, Ltd.     

Slow Passage through a Saddle-Center Bifurcation

Received September 24, 1998; accepted August 9, 1999     

Anti-control of Hopf Bifurcation in a Delayed Predator-prey Gompertz Model

A delayed predator-prey Gompertz model is investigated. The stability is ana- lyzed. Anti-control of Hopf bifurcation for the model is presented. Numerical simulations are performed to confirm that the new feedback controller using time delay is efficient in creating Hopf bifurcation. Finally, main conclusions are included.     

Bifurcation Instability of sheet metal during spring-back

In automotive and home appliance industries, there are many complex-shaped sheet metal components which need to be fabricated in multiple stamping operations. For example, the manufacturing of an outer case of washing machine consists of stamping followed by a bending operation. After the first stage of the stamping process, a large amount of spring-back takes place, and therefore, it is difficult to proceed to the next stage of the bending process. In the stamping process of that kind of sheet component with low geometric constraint, the forming area is large compared to the forming depth. Therefore, the formed part is in an unstable state and is less geometrically constrained, which causes a large amount of spring-back. To investigate this phenomenon, finite element analyses are carried out. During a spring-back analysis after forming, bifurcation takes place and the finite element solution procedure using the Newton–Raphson scheme becomes unstable. To get a stable post-bifurcation solution, a bifurcation algorithm is introduced at the bifurcation point. The deformed shapes obtained from finite element analyses are in good agreement with the experimental data. From this study, it is shown that the bifurcation behaviour enlarges the spring-back and the degree of dimensional error. To obtain additional possible post-bifurcation solutions, non-bifurcation analyses using initial guesses obtained in a modal analysis are carried. For the initial guesses, lowed four eigenmodes are utilized. Finally, the post-bifurcation behaviour and spring-back amount are investigated for various process parameters including the forming depth, punch width and corner radius.     

基于L-S方法求解一类方程

基于Lyapunov-Schmidt方法求出给定方程的分岐方程,Newton迭代得到其在分岐点附近的近似非平凡解枝,得到了满意的结果.     

Period doubling toward chaos in a driven magnetic macrospin

The Landau-Lifshitz-Gilbert equation is analyzed in the case of a configuration involving easy plane isotropy under the influence of a sinusoidally oscillating magnetic field and a demagnetizing field. Through the use of numerical techniques, chaotic behavior is found and analyzed. By reducing the system to a discrete map (numerically), bifurcation diagrams for the system are computed. The system is found to exhibit a period doubling cascade route to chaos, and it obeys certain convergence rules for chaotic transitions outlined by Feigenbaum. A connection is drawn between the route to chaos and the geometry of the system, and comparisons are made with similar systems. Within the chaotic regime, windows of arbitrarily large period are suspected to exist, and explicitly illustrated and discussed for a period three window.     

A technique for flow-enhanced magnetic resonance angiography of the lower extremities

A two-dimensional, flow-enhanced gradient echo pulse sequence for nuclear magnetic resonance angiography is described. It employs interleaved, presaturated slices to acquire data efficiently on imagers which favor interleaved acquisition over sequential acquisition for multislice imaging. It is useful on any imager when the effective TR is extended to enhance the sensitivity to slow flow. The technique was applied to the region from aortic bifurcation to the iliac bifurcations of three normal volunteers. The right and left common iliac arteries and veins, the separation of the external and internal iliac arteries, and secondary branches were clearly depicted.     

Existence and bifurcation of integral manifolds with applications

In this paper a bifurcation theorem on the existence of integral manifolds is obtained by using contracting principle. As an application, sufficient conditions for a higher dimensional system to have an integral manifold are given. Especially the existence and uniqueness of a 3-dimensional invariant torus appearing in a 4-dimensional autonomous system with singularity of codimension two are proved.