The influence of the non-Newtonian properties of blood on the flow in large arteries: unsteady flow in a 90 degrees curved tube.

A numerical and experimental investigation of unsteady entry flow in a 90 degrees curved tube is presented to study the impact of the non-Newtonian properties of blood on the velocity distribution. The time-dependent flow rate for the Newtonian and the non-Newtonian blood analog fluid were identical. For the numerical computation, a Carreau-Yasuda model was employed to accommodate the shear thinning behavior of the Xanthan gum solution. The viscoelastic properties were not taken into account. The experimental results indicate that significant differences between the Newtonian and non-Newtonian fluid are present. The numerical results for both the Newtonian and the non-Newtonian fluid agree well with the experimental results. Since viscoelasticity was not included in the numerical code, shear thinning behavior of the blood analog fluid seems to be the dominant non-Newtonian property, even under unsteady flow conditions. Finally, a comparison between the non-Newtonian fluid model and a Newtonian fluid at a rescaled Reynolds number is presented. The rescaled Reynolds number, based on a characteristic rather than the high-shear rate viscosity of the Xanthan gum solution, was about three times as low as the original Reynolds number. Comparison reveals that the character of flow of the non-Newtonian fluid is simulated quite well by using the appropriate Reynolds number.     

Mental representation of number in different numerical forms

How are numerical operations implemented within the human brain? It has been suggested that there are at least three different codes for representing number: a verbal code that is used to manipulate number words and perform mental numerical operations (e.g., multiplication), a visual code that is used to decode frequently used visual number forms (e.g., Arabic digits), and an abstract analog code that may be used to represent numerical quantities. Furthermore, each of these codes is associated with a different neural substrate. We extend these studies using dense-sensor event-related EEG recording techniques to investigate the temporal pattern of notation-specific effects observed in a parity judgement (odd versus even) task in which single numbers were presented in one of four different numerical notations. Contrasts between different notations demonstrated clear modulations in the visual evoked potentials (VEP) recorded. We observed increased amplitudes for the P1 and N1 components of the VEP that were specific to Arabic numerals and to dot configurations but differed for random and recognizable (die-face) dot configurations. These results demonstrate clear, notation-specific differences in the time course of numerical information processing and provide electrophysiological support for the triple-code model of numerical representation.     

Program Code Generator for Cardiac Electrophysiology Simulation with Automatic PDE Boundary Condition Handling

Clinical and experimental studies involving human hearts can have certain limitations. Methods such as computer simulations can be an important alternative or supplemental tool. Physiological simulation at the tissue or organ level typically involves the handling of partial differential equations (PDEs). Boundary conditions and distributed parameters, such as those used in pharmacokinetics simulation, add to the complexity of the PDE solution. These factors can tailor PDE solutions and their corresponding program code to specific problems. Boundary condition and parameter changes in the customized code are usually prone to errors and time-consuming. We propose a general approach for handling PDEs and boundary conditions in computational models using a replacement scheme for discretization. This study is an extension of a program generator that we introduced in a previous publication. The program generator can generate code for multi-cell simulations of cardiac electrophysiology. Improvements to the system allow it to handle simultaneous equations in the biological function model as well as implicit PDE numerical schemes. The replacement scheme involves substituting all partial differential terms with numerical solution equations. Once the model and boundary equations are discretized with the numerical solution scheme, instances of the equations are generated to undergo dependency analysis. The result of the dependency analysis is then used to generate the program code. The resulting program code are in Java or C programming language. To validate the automatic handling of boundary conditions in the program code generator, we generated simulation code using the FHN, Luo-Rudy 1, and Hund-Rudy cell models and run cell-to-cell coupling and action potential propagation simulations. One of the simulations is based on a published experiment and simulation results are compared with the experimental data. We conclude that the proposed program code generator can be used to generate code for physiological simulations and provides a tool for studying cardiac electrophysiology.     

GPUs, a new tool of acceleration in CFD: efficiency and reliability on smoothed particle hydrodynamics methods

Smoothed Particle Hydrodynamics (SPH) is a numerical method commonly used in Computational Fluid Dynamics (CFD) to simulate complex free-surface flows. Simulations with this mesh-free particle method far exceed the capacity of a single processor. In this paper, as part of a dual-functioning code for either central processing units (CPUs) or Graphics Processor Units (GPUs), a parallelisation using GPUs is presented. The GPU parallelisation technique uses the Compute Unified Device Architecture (CUDA) of nVidia devices. Simulations with more than one million particles on a single GPU card exhibit speedups of up to two orders of magnitude over using a single-core CPU. It is demonstrated that the code achieves different speedups with different CUDA-enabled GPUs. The numerical behaviour of the SPH code is validated with a standard benchmark test case of dam break flow impacting on an obstacle where good agreement with the experimental results is observed. Both the achieved speed-ups and the quantitative agreement with experiments suggest that CUDA-based GPU programming can be used in SPH methods with efficiency and reliability.     

Calculation of one-dimensional plasma flows with allowance for the kinetics of ion collisions

The influence of kinetic effects on the generation of line X radiation during the spherical implosion of a laser corona plasma with two ion species is studied under the conditions prevailing in experiments with thin-wall spherical targets in the Iskra-5 laser facility of the All-Russia Research Institute of Experimental Physics (Sarov, Russia). Kinetic processes occurring in a multicharged plasma are investigated using a specially devised code for solving one-dimensional Landau equations for a nondegenerate multicomponent plasma by the Monte Carlo method (the KIN-MC code). The code was developed using the quasineutral plasma approximation under the assumption that the electron distribution function is locally equilibrium. The model equations are presented, the scheme of numerical solution is described, and the calculated results are discussed.     

Empirical validation of SAR values predicted by FDTD modeling.

Rapid increase in the use of numerical techniques to predict current density or specific absorption rate (SAR) in sophisticated three dimensional anatomical computer models of man and animals has resulted in the need to understand how numerical solutions of the complex electrodynamics equations match with empirical measurements. This aspect is particularly important because different numerical codes and computer models are used in research settings as a guide in designing clinical devices, telecommunication systems, and safety standards. To ensure compliance with safety guidelines during equipment design, manufacturing and maintenance, realistic and accurate models could be used as a bridge between empirical data and actual exposure conditions. Before these tools are transitioned into the hands of health safety officers and system designers, their accuracy and limitations must be verified under a variety of exposure conditions using available analytical and empirical dosimetry techniques. In this paper, empirical validation of SAR values predicted by finite difference time domain (FDTD) numerical code on sphere and rat is presented. The results of this study show a good agreement between empirical and theoretical methods and, thus, offer a relatively high confidence in SAR predictions obtained from digital anatomical models based on the FDTD numerical code.     

A digital model for the venous junctions

The venous network in the lower limbs is composed of a considerable number of confluent junctions. Each of these singularities introduces some blood flow disturbances. Each physiological junction is unique, in terms of its geometry as well as the blood flow rate. In order to account for this great variability, we developed a numerical model based on the use of the N3S code (a software package for solving Navier-Stokes equations). To test the validity of the model, one of the numerical simulations is compared with the data obtained in the corresponding experimental configuration. The velocity measurements were carried out with an ultrasonic pulsed Doppler velocimeter. We also measured pressure differences using differential sensors. The numerical computations were then used to obtain the values of the flow variables at any point, with various geometrical and flow configurations. As far as the velocity field is concerned, a very marked three-dimensional pattern with swirls was observed. The pressure evolution was also strongly disturbed, with a non-linear decrease. All these data indicate that confluence effects cannot be neglected when evaluating pressure decreases. With a tool of this kind, it is possible to accurately predict the disturbances associated with any geometrical configuration or any flow rate.     

Analysis of the effect of swimmer's head position on swimming performance using computational fluid dynamics

The aim of this numerical work is to analyze the effect of the position of the swimmer's head on the hydrodynamic performances in swimming. In this initial study, the problem was modeled as 2D and in steady hydrodynamic state. The geometry is generated by the CAD software CATIA and the numerical simulation is carried out by the use of the CFD Fluent code. The standard k-epsilon turbulence model is used with a specific wall law. Three positions of the head were studied, for a range of Reynolds numbers about 10(6). The obtained numerical results revealed that the position of the head had a noticeable effect on the hydrodynamic performances, strongly modifying the wake around the swimmer. The analysis of these results made it possible to propose an optimal position of the head of a swimmer in underwater swimming.     

A Digital Model for the Venous Junctions

The venous network in the lower limbs is composed of a considerable number of confluent junctions. Each of these singularities introduces some blood flow disturbances. Each physiological junction is unique, in terms of its geometry as well as the blood flow rate. In order to account for this great variability, we developed a numerical model based on the use of the N3S code (a software package for solving Navier-Stokes equations). To test the validity of the model, one of the numerical simulations is compared with the data obtained in the corresponding experimental configuration. The velocity measurements were carried out with an ultrasonic pulsed Doppler velocimeter. We also measured pressure differences using differential sensors. The numerical computations were then used to obtain the values of the flow variables at any point, with various geometrical and flow configurations. As far as the velocity field is concerned, a very marked three-dimensional pattern with swirls was observed. The pressure evolution was also strongly disturbed, with a non-linear decrease. All these data indicate that confluence effects cannot be neglected when evaluating pressure decreases. With a tool of this kind, it is possible to accurately predict the disturbances associated with any geometrical configuration or any flow rate.     

Identification of the plasma boundary shape and position in the Damavand tokamak

A series of experiments and numerical calculations have been done on the Damavand tokamak for accurate determination of equilibrium parameters, such as the plasma boundary position and shape. For this work, the pickup coils of the Damavand tokamak were recalibrated and after that a plasma boundary shape identification code was developed for analyzing the experimental data, such as magnetic probes and coils currents data. The plasma boundary position, shape and other parameters are determined by the plasma shape identification code. A free-boundary equilibrium code was also generated for comparison with the plasma boundary shape identification results and determination of required fields to obtain elongated plasma in the Damavand tokamak.     

核素迁移的二维非平衡吸附数值模式研究

为了进一步研究核素在孔隙介质中的迁移规律,本文将水动力弥散方程与核素在孔隙介质中非平衡吸附关系式相耦合,建立了非平衡吸附数值模式,并给出了非平衡吸附数值模式中吸附解吸速率参数的获取方法．为了检验所建模式与编制的程序的正确性,利用试验场（８５）Ｓｒ迁移的野外试验资料进行了计算,并与核素迁移平衡吸附数值模式的计算结果进行了对比．结果表明,二维非平衡吸附核素迁移数值模式能够更好地描述核素迁移实际．     

A numerical simulation of intimal thickening under shear in arteries.

A model of intima thickening proposed by Friedman and his coworkers (1,2) is incorporated in our computer code to simulate the growth of intima under shear. The computer code is based on a finite volume method in a boundary-fitted coordinate system. It is found that the model yields an evenly-distributed thickening over a straight, smooth vessel wall. However, in a stenosed or a curved artery, thicker intima can be formed in preferential regions due to unevenly-distributed wall shear stresses. The results clearly demonstrate the correlations among the geometry, wall shear rate and the plaque localization in arteries. The model is applied to a straight artery with a stenosis or sinus, a smooth curved artery and a stenosed curved artery. The effects of stenosis/sinus and lumen curvature on the flows and the intimal thickening are studied. The simulation provides a numerical visualization of the intimal thickening in a dynamic way.     

Vascular wall flow-induced forces in a progressively enlarged aneurysm model

The current study is focused on the numerical investigation of the flow field induced by the unsteady flow in the vicinity of an abdominal aortic aneurysm model. The computational fluid dynamics code used is based on the finite volume method, and it has already been used in various bioflow studies. For modelling the rheological behaviour of blood, the Quemada non-Newtonian model is employed, which is suitable for simulating the two-phase character of blood namely a suspension of blood cells in plasma. For examining its non-Newtonian effects a comparison with a corresponding Newtonian flow is carried out. Furthermore, the investigation is focused on the distribution of the flow-induced forces on the interior wall of the aneurysm and in order to study the development of the distribution with the gradual enlargement of the aneurysm, three different degrees of aneurysm-growth have been assumed. Finally and for examining the effect of the distribution on the aneurysm growth, a comparison is made between the pressure and wall shear-stress distributions at the wall for each growth-degree.     

Vascular wall flow-induced forces in a progressively enlarged aneurysm model

The current study is focused on the numerical investigation of the flow field induced by the unsteady flow in the vicinity of an abdominal aortic aneurysm model. The computational fluid dynamics code used is based on the finite volume method, and it has already been used in various bioflow studies. For modelling the rheological behaviour of blood, the Quemada non-Newtonian model is employed, which is suitable for simulating the two-phase character of blood namely a suspension of blood cells in plasma. For examining its non-Newtonian effects a comparison with a corresponding Newtonian flow is carried out. Furthermore, the investigation is focused on the distribution of the flow-induced forces on the interior wall of the aneurysm and in order to study the development of the distribution with the gradual enlargement of the aneurysm, three different degrees of aneurysm-growth have been assumed. Finally and for examining the effect of the distribution on the aneurysm growth, a comparison is made between the pressure and wall shear-stress distributions at the wall for each growth-degree.     

Dependence of leukocyte capture on instantaneous pulsatile flow

Atherosclerosis, an artery disease, is currently the leading cause of death in the United States in both men and women. The first step in the development of atherosclerosis involves leukocyte adhesion to the arterial endothelium. It is broadly accepted that blood flow, more specifically wall shear stress (WSS), plays an important role in leukocyte capture and subsequent development of an atherosclerotic plaque. What is less known is how instantaneous WSS, which can vary by up to 5 Pa over one cardiac cycle, influences leukocyte capture. In this paper we use direct numerical simulations (DNS), performed using an in-house code, to illustrate that leukocyte capture is different whether as a function of instantaneous or time-averaged blood flow. Specifically, a stenotic plaque is modeled using a computational fluid dynamics (CFD) solver through fully three-dimensional Navier-Stokes equations and the immersed boundary method. Pulsatile triphasic inflow is used to simulate the cardiac cycle. The CFD is coupled with an agent-based leukocyte capture model to assess the impact of instantaneous hemodynamics on stenosis growth. The computed wall shear stress agrees well with the results obtained with a commercial software, as well as with theoretical results in the healthy region of the artery. The analysis emphasizes the importance of the instantaneous flow conditions in evaluating the leukocyte rate of capture. That is, the capture rate computed from mean flow field is generally underpredicted compared to the actual rate of capture. Thus, in order to obtain a reliable estimate, the flow unsteadiness during a cardiac cycle should be taken into account.     

Magnetic induction at 60 Hz in the human heart: a comparison between the in situ and isolated scenarios

Numerical modelling is used to estimate the electric fields and currents induced in the human heart and associated major blood vessels by 60 Hz external magnetic fields. The modelling is accomplished using a scalar-potential finite-difference code applied to a 3.6-mm resolution voxel-based model of the whole human body. The main goal of the present work is a comparison between the induced field levels in the heart located in situ and in isolation. This information is of value in assessing any health risks due to such fields, given that some existing protection standards consider the heart as an isolated conducting body. It is shown that the field levels differ significantly between these two scenarios. Consequently, data from more realistic and detailed numerical studies are required for the development of reliable standards.     

A numerical study of the shape of the surface separating flow into branches in microvascular bifurcations.

The shape of the separating surface formed by the streamlines entering the branches of microvascular bifurcations plays a major role in determining the distribution of red blood cells and other blood constituents downstream from the bifurcation. Using the finite element method, we determined the shape of the surface through numerical solution of three dimensional Navier-Stokes equations for fluid flow at low Reynolds numbers in a T-type bifurcation of circular tubes. Calculations were done for a wide range of daughter branch to parent vessel diameter ratios and flow ratios. The effect of Reynolds number was also studied. Our numerical results are in good agreement with previously reported experimental data of Rong and Carr (Microvascular Research, Vol. 39, pp. 186-202, 1990). The numerical results of this study will be used to predict the concentration of blood constituents downstream from microvascular bifurcations providing that the inlet concentration profile is known.     

Genetic code,hamming distance and stochastic matrices

In this paper we use the Gray code representation of the genetic code C = 00, U = 10, G = 11 and A = 01 (C pairs with G, A pairs with U) to generate a sequence of genetic code-based matrices. In connection with these code-based matrices, we use the Hamming distance to generate a sequence of numerical matrices. We then further investigate the properties of the numerical matrices and show that they are doubly stochastic and symmetric. We determine the frequency distributions of the Hamming distances, building blocks of the matrices, decomposition and iterations of matrices. We present an explicit decomposition formula for the genetic code-based matrix in terms of permutation matrices, which provides a hypercube representation of the genetic code. It is also observed that there is a Hamiltonian cycle in a genetic code-based hypercube.     

Study of plasma heating in discharges with neutral beam injection in the Globus-M spherical tokamak

Results from experimental studies on the injection of high-energy neutral hydrogen beams into the plasma of the Globus-M spherical tokamak are reviewed. In the Introduction, the importance of these studies for implementing the controlled fusion research program and constructing the ITER tokamak is proved. Some problems related to the use of neutral beam injection in small and low-aspect-ratio tokamaks is analyzed. Results are presented from numerical simulations of the experiment by using the ASTRA transport code. It is shown that the use of neutral beam injection in the Globus-M tokamak ensures efficient ion heating and increases the plasma stored energy. The greater part of the review is devoted to the survey of experiments on the injection of 22-to 30-keV hydrogen and deuterium beams with a power of 0.4–0.8 MW into the plasma of the Globus-M spherical tokamak in a wide range of plasma currents and densities. The experimental results are analyzed and compared with the results of numerical simulations. The achievement of top plasma parameters is highlighted.