Multielemental Analysis of Various Kinds of Whisky.

Whisky (whiskey) consists of many trace elements coming from the raw materials used in its fermentation, distillation and maturation processes. These ingredients assure the exceptional organoleptic characteristics of the beverage. Their analysis is important to better control the stages of fermentation, distillation, taste repeatability and for product quality assurance as well as from the brand protection point of view. This article presents the usefulness of modern analytical techniques based on elemental analysis. ICP mass spectrometry and CV atomic absorption spectroscopy were applied to distinguish whisky produced in Scotland from whisky coming from Ireland and the United States. The collected semi-quantitative data were used for multivariate analysis performed using the Statistica 10.0 software. The results showed that Irish whiskey is characterized by quite a high amount of Ba and Ti compared with other samples, which made it possible to distinguish this sample from the others. No strict correlation was found between the type of whisky and the amount of trace elements, however, the projection of objects on the first two components revealed that single malt samples created one cluster.

Numerical flow experiment for assessing predictors for cerebrovascular accidents in patients with PHACES syndrome.

There is an increased risk of cerebrovascular accidents (CVA) in individuals with PHACES, yet the precise causes are not well understood. In this analysis, we aimed to examine the role of arteriopathy in PHACES syndrome as a potential contributor to CVA. We analyzed clinical and radiological data from 282 patients with suspected PHACES syndrome. We analyzed clinical features, including the presence of infantile hemangioma and radiological features based on magnetic resonance angiography or computed tomography angiography, in individuals with PHACES syndrome according to the Garzon criteria. To analyze intravascular blood flow, we conducted a simulation based on the Fluid-Structure Interaction (FSI) method, utilizing radiological data. The collected data underwent statistical analysis. Twenty patients with PHACES syndrome were included. CVAs were noted in 6 cases. Hypoplasia (p = 0.03), severe tortuosity (p < 0.01), absence of at least one main cerebral artery (p < 0.01), and presence of persistent arteries (p = 0.01) were associated with CVAs, with severe tortuosity being the strongest predictor. The in-silico analysis showed that the combination of hypoplasia and severe tortuosity resulted in a strongly thrombogenic environment. Severe tortuosity, combined with hypoplasia, is sufficient to create a hemodynamic environment conducive to thrombus formation and should be considered high-risk for cerebrovascular accidents (CVAs) in PHACES patients.

Intracranial bypass for giant aneurysms treatment assessed by computational fluid dynamics (CFD) analysis.

Unruptured giant intracranial aneurysms (GIA) are those with diameters of 25 mm or greater. As aneurysm size is correlated with rupture risk, GIA natural history is poor. Parent artery occlusion or trapping plus bypass revascularization should be considered to encourage intra-aneurysmal thrombosis when other treatment options are contraindicated. The mechanistic background of these methods is poorly studied. Thus, we assessed the potential of computational fluid dynamics (CFD) and fluid-structure interaction (FSI) analyses for clinical use in the preoperative stage. A CFD investigation in three patient-specific flexible models of whole arterial brain circulation was performed. A C6 ICA segment GIA model was created based on CT angiography. Two models were then constructed that simulated a virtual bypass in combination with proximal GIA occlusion, but with differing middle cerebral artery (MCA) recipient vessels for the anastomosis. FSI and CFD investigations were performed in three models to assess changes in flow pattern and haemodynamic parameters alternations (wall shear stress (WSS), oscillatory shear index (OSI), maximal time averaged WSS (TAWSS), and pressure). General flow splitting across the entire domain was affected by virtual bypass procedures, and any deficiency was partially compensated by a specific configuration of the circle of Willis. Following the implementation of bypass procedures, a reduction in haemodynamic parameters was observed within the aneurysm in both cases under analysis. In the case of the temporal MCA branch bypass, the decreases in the studied parameters were slightly greater than in the frontal MCA branch bypass. The reduction in the magnitude of the chosen area-averaged parameters (averaged over the aneurysm wall surface) was as follows: WSS 35.7%, OSI 19.0%, TAWSS 94.7%, and pressure 24.2%. FSI CFD investigation based on patient-specific anatomy models with subsequent stimulation of virtual proximal aneurysm occlusion in conjunction with bypass showed that this method creates a pro-thrombotic favourable environment whilst reducing intra-aneurysmal pressure leading to shrinking. MCA branch recipient selection for optimum haemodynamic conditions should be evaluated individually in the preoperative stage.

Computational fluid dynamics; a new diagnostic tool in giant intracerebral aneurysm treatment.

Giant intracerebral aneurysms (GIA) comprise up to 5 % of all intracranial aneurysms. The indirect surgical strategy, which leaves the GIA untouched but reverses the blood flow by performing a bypass in combination with proximal parent artery occlusion is a useful method to achieve spontaneous aneurysm occlusion. The goal of this study was to assess the utility of computational fluid dynamics (CFD) in preoperative GIA treatment planning. We hypothesise that CFD simulations will predict treatment results. A fluid-structure interaction (FSI) CFD investigation was performed for the entire arterial brain circulation. The analyses were performed in three patient-specific CT angiogram models. The first served as the reference geometry with a C6 internal carotid artery (ICA) GIA, the second a proximal parent artery occlusion (PAO) and virtual bypass to the frontal M2 branch of the middle cerebral artery (MCA), and the third a proximal PAO in combination with a temporal M2 branch bypass. The volume of "old blood", flow residence time (FRT), dynamic viscosity and haemodynamic changes were also analysed. The "old blood" within the aneurysm in the bypass models reached 41 % after 20 cardiac cycles while in the reference model it was fully washed out. In Bypass 2 "old blood" was also observed in the main trunk of the MCA after 20 cardiac cycles. Extrapolation of the results yielded a duration of 4 years required to replace the "old blood" inside the aneurysm after bypass revascularization. In both bypass models a 7-fold increase in mean blood viscosity in the aneurysm region was noted. Bypass revascularization combined with proximal PAO favours thrombosis. Areas prone to thrombus formation, and subsequently the treatment outcomes, were accurately identified in the preoperative model. Virtual surgical operations can give a remarkable insight into haemodynamics that could support operative decision-making.

The Elemental Fingerprints of Different Types of Whisky as Determined by ICP-OES and ICP-MS Techniques in Relation to Their Type, Age, and Origin.

A total of 170 samples of whisky from 11 countries were analysed in terms of their elemental profiles. The levels of 31 elements were determined by Inductively Coupled Plasma Mass Spectrometry (ICP-MS): Ag, Al, B, Ba, Be, Bi, Cd, Co, Cr, Cu, Li, Mn, Mo, Ni, Pb, Sb, Sn, Sr, Te, Tl, U, and V, Inductively Coupled Plasma Optical Emission Spectrometry (ICP-OES) Ca, Fe, K, Mg, P, S, Ti, and Zn and Cold Vapor-Atomic Absorption (CV-AAS): Hg techniques in those alcoholic samples. A comparative analysis of elemental profiles was made on the basis of the content of chosen elements with regard to selected parameters: country of origin, type of whisky (single malt and blended) and age of products. One of the elements which clearly distinguishes single malt and blended types of whisky is copper. Single malt Scotch whisky had a uniform concentration of copper, which is significantly higher for all malt whisky samples when compared with the blended type. Analysis of samples from the USA (n = 26) and Ireland (n = 15) clearly revealed that the objects represented by the same product but originating from independent bottles (e.g., JB, JDG, BUS brands) show common elemental profiles. On the other hand, comparative analysis of Scotch whisky with respect to aging time revealed that the longer the alcohol was aged, (i.e., the longer it stayed in the barrel), the higher the content of Cu and Mn that was recorded.

Assessment of the Authenticity of Whisky Samples Based on the Multi-Elemental and Multivariate Analysis.

Two hundred and five samples of whisky, including 170 authentic and 35 fake products, were analyzed in terms of their elemental profiles in order to distinguish them according to the parameter of their authenticity. The study of 31 elements (Ag, Al, B, Ba, Be, Bi, Cd, Co, Cr, Cu, Li, Mn, Mo, Ni, Pb, Sb, Sn, Sr, Te, Tl, U, V, Ca, Fe, K, Mg, P, S, Ti and Zn) was performed using the Inductively Coupled Plasma Mass Spectrometry (ICP-MS), Inductively Coupled Plasma Optical Emission Spectrometry (ICP-OES) and Cold Vapor-Atomic Absorption (CVAAS) techniques. Additionally, the pH values of all samples were determined by pH-meter, and their isotopic ratios of 88Sr/86Sr, 84Sr/86Sr, 87Sr/86Sr and 63Cu/65Cu were assessed, based on the number of counts by ICP-MS. As a result of conducted research, elements, such as Mn, K, P and S, were identified as markers of whisky adulteration related to the age of alcohol. The concentrations of manganese, potassium and phosphorus were significantly lower in the fake samples (which were not aged, or the aging period was much shorter than legally required), compared to the original samples (in all cases subjected to the aging process). The observed differences were related to the migration of these elements from wooden barrels to the alcohol contained in them. On the other hand, the sulfur concentration in the processed samples was much higher in the counterfeit samples than in the authentic ones. The total sulfur content, such as that of alkyl sulfides, decreases in alcohol with aging in the barrels. Furthermore, counterfeit samples can be of variable origin and composition, so they cannot be characterized as one group with identical or comparable features. Repeatedly, the element of randomness dominates in the production of these kinds of alcohols. However, as indicated in this work, the extensive elemental analysis supported by statistical tools can be helpful, especially in the context of detecting age-related adulteration of whisky. The results presented in this paper are the final part of a comprehensive study on the influence of selected factors on the elemental composition of whisky.

Porous Media Computational Fluid Dynamics and the Role of the First Coil in the Embolization of Ruptured Intracranial Aneurysms.

BACKGROUND: The objective of our project was to identify a late recanalization predictor in ruptured intracranial aneurysms treated with coil embolization. This goal was achieved by means of a statistical analysis followed by a computational fluid dynamics (CFD) with porous media modelling approach. Porous media CFD simulated the hemodynamics within the aneurysmal dome after coiling. METHODS: Firstly, a retrospective single center analysis of 66 aneurysmal subarachnoid hemorrhage patients was conducted. The authors assessed morphometric parameters, packing density, first coil volume packing density (1st VPD) and recanalization rate on digital subtraction angiograms (DSA). The effectiveness of initial endovascular treatment was visually determined using the modified Raymond-Roy classification directly after the embolization and in a 6- and 12-month follow-up DSA. In the next step, a comparison between porous media CFD analyses and our statistical results was performed. A geometry used during numerical simulations based on a patient-specific anatomy, where the aneurysm dome was modelled as a separate, porous domain. To evaluate hemodynamic changes, CFD was utilized for a control case (without any porosity) and for a wide range of porosities that resembled 1-30% of VPD. Numerical analyses were performed in Ansys CFX solver. RESULTS: A multivariate analysis showed that 1st VPD affected the late recanalization rate (p < 0.001). Its value was significantly greater in all patients without recanalization (p < 0.001). Receiver operating characteristic curves governed by the univariate analysis showed that the model for late recanalization prediction based on 1st VPD (AUC 0.94 (95%CI: 0.86-1.00) is the most important predictor of late recanalization (p < 0.001). A cut-off point of 10.56% (sensitivity-0.722; specificity-0.979) was confirmed as optimal in a computational fluid dynamics analysis. The CFD results indicate that pressure at the aneurysm wall and residual flow volume (blood volume with mean fluid velocity > 0.01 m/s) within the aneurysmal dome tended to asymptotically decrease when VPD exceeded 10%. CONCLUSIONS: High 1st VPD decreases the late recanalization rate in ruptured intracranial aneurysms treated with coil embolization (according to our statistical results > 10.56%). We present an easy intraoperatively calculable predictor which has the potential to be used in clinical practice as a tip to improve clinical outcomes.

Risk Factors for Recanalization after Coil Embolization.

The aim of our study was to identify risk factors for recanalization 6 months after coil embolization using clinical data followed by computational fluid dynamics (CFD) analysis. METHODS: Firstly, clinical data of 184 patients treated with coil embolization were analyzed retrospectively. Secondly, aneurysm models for high/low recanalization risk were generated based on ROC curves and their cut-off points. Afterward, CFD was utilized to validate the results. RESULTS: In multivariable analysis, aneurysm filling during the first embolization was an independent risk factor whilst packing density was a protective factor of recanalization after 6 months in patients with aSAH. For patients with unruptured aneurysms, packing density was found to be a protective factor whilst the aneurysm neck size was an independent risk factor. Complex flow pattern and multiple vortices were associated with aneurysm shape and were characteristic of the high recanalization risk group. CONCLUSIONS: Statistical analysis suggested that there are various factors influencing recanalization risk. Once certain values of morphometric parameters are exceeded, a complex flow with numerous vortices occurs. This phenomenon was revealed due to CFD investigations that validated our statistical research. Thus, the complex flow pattern itself can be treated as a relevant recanalization predictor.

Ultrafast optical circuit switching for data centers using integrated soliton microcombs.

Due to the slowdown of Moore's law, it will become increasingly challenging to efficiently scale the network in current data centers utilizing electrical packet switches as data rates grow. Optical circuit switches (OCS) represent an appealing option to overcome this issue by eliminating the need for expensive and power-hungry transceivers and electrical switches in the core of the network. In particular, optical switches based on tunable lasers and arrayed waveguide grating routers are quite promising due to the use of a passive core, which increases fault tolerance and reduces management overhead. Such an OCS-network can offer high bandwidth, low network latency and an energy-efficient and scalable data center network. To support dynamic data center workloads efficiently, however, it is critical to switch between wavelengths at nanosecond (ns) timescales. Here we demonstrate ultrafast OCS based on a microcomb and semiconductor optical amplifiers (SOAs). Using a photonic integrated Si3N4 microcomb, sub-ns (<520 ps) switching along with the 25-Gbps non-return-to-zero (NRZ) and 50-Gbps four-level pulse amplitude modulation (PAM-4) burst mode data transmission is achieved. Further, we use a photonic integrated circuit comprising an Indium phosphide based SOA array and an arrayed waveguide grating to show sub-ns switching (<900 ps) along with 25-Gbps NRZ burst mode transmission providing a path towards a more scalable and energy-efficient wavelength-switched network for data centers in the post Moore's Law era.

A two-stage model of an arteriovenous fistula maturation process.

PURPOSE: An arteriovenous fistula has been a widely accepted vascular access for hemodialysis, however, a fistula maturation process is still not fully understood. In the short period of time, right after vein and artery shunting, the physical and biological changes take place mainly in the venous wall. A two-stage modeling method of arteriovenous fistula maturation process was proposed and presented. METHODS: The first stage of the maturation was modeled with two-way coupled fluid structure interaction computer simulations. Whereas for the second, biological stage, a model was based on the change in the elasticity of the venous wall due to wall shear stress (WSS) modifications. RESULTS: The relation between stress and radial and circumferential strain, based on Lame's theory, makes possible to introduce a mathematical model defining modulus of elasticity, averaged WSS, and venous diameter as time functions. The presented model enables one to predict changes in the monitored parameters in the arteriovenous fistula taking place in the time longer than 90 days. CONCLUSIONS: We found that probably the majority of fistulas can be assessed to be mature too early, when the adequate blood flow rate is achieved but mean WSS still remains at the non-physiological level (>10 Pa).

Inhaled drug airflow patterns and particles deposition in the paediatric respiratory tract.

PURPOSE: The effectiveness of inhaled drugs is strictly related to areas reachable by drug particles. Unless particles reach the desired part of the bronchial tree, their influence might not meet the expectations. Consequently, the disease progress might not be stopped or even slowed down. Therefore, the primary objective of this research was to analyze the airflow patterns and particle deposition of a standard inhaled drug using computational fluid dynamics. METHODS: The study was devoted to the analysis of the particle diameter influence on their deposition areas within the entire respiratory tract. Two patient-specific respiratory tract models, for 6 and 12-year-old patients, were reconstructed based on the computed tomography examinations. Numerical analyses were carried out as stationary ones with the constant inflow of the particles of various diameters (within the range of 1-50 µm). It was proven that depending on the particle size, their deposition within the respiratory tract varies significantly. RESULTS: The vast majority of the particles with diameters over 20 µm is gathered on the walls of the throat, whereas particles of diameters 5-15 µm are accumulated mainly on the trachea walls, leaving the alveoli insufficiently supplied with the drug particles. CONCLUSIONS: The inhaled drug size cannot be treated as negligible factor during the drug spraying. An improper distribution of the particles might not inhibit the symptoms of the asthma. Numerical simulations may improve drugs selection and visualize their distribution along the airways, which might accelerate asthma treatment personalization.

Methods for determination of stagnation in pneumatic ventricular assist devices.

BACKGROUND: A pneumatic paediatric ventricular assist device developed at the Foundation of Cardiac Surgery Development, Zabrze, equipped with valves based on J. Moll's design, with later modifications introduced at the Institute of Turbomachinery, Lodz University of Technology, was tested numerically and experimentally. The main aim of those investigations was to detect stagnation zones within the ventricular assist device and indicate advantages and limitations of both approaches. METHODS: In the numerical transient test, a motion of the diaphragm and discs was simulated. Two different methods were used to illustrate stagnation zones in the ventricular assist device. The flow pattern inside the chamber was represented by velocity contours and vectors to validate the results using images obtained in the laser particle image velocimetry experiment. RESULTS: The experimental light-based method implied problems with proper illumination of regions in the wall vicinity. High-resolution flow data and other important parameters as stagnation regions or flow patterns in regions not accessible for light in the particle image velocimetry method are derived in the numerical solution. However, computations of a single case are much more time-consuming if compared to a single experiment conducted on a well-calibrated stand. CONCLUSION: The resulting main vortexes in the central part of the pump chamber and the velocity magnitudes are correlated in both methods, which are complementary and when used together offer better insight into the flow structure inside the ventricular assist device and enable a deeper analysis of the results.

Experimental investigations of the aerated polymethylmethacrylate-based vertebral cement flow in capillaries.

INTRODUCTION: A constant growth in the population suffering from osteoporotic vertebral weakening is observed. As a result, vertebroplasty procedures become more and more common. Unfortunately, they may be associated with several complications occurring during bone cement injection, including its leakage or overheating of tissues. Despite several experimental studies, there is a lack of data related to random aeration of the bone cement. Therefore, the main objective of the following investigations was to emphasize that random aeration of the bone cement, and, consequently, a compressibility factor, could not be treated as a negligible factor during the vertebroplasty procedure and had to be taken into account in the development of the mathematical model. MATERIALS AND METHODS: A special test rig to reproduce the vertebroplasty procedure was designed and built. The authors conducted numerous experiments on polymethylmethacrylate-based bone cement flows, analyzing different flow conditions, such as volume flow rate and flow channel diameter. Time periods of the flow front between characteristic sections and pressures (differential and gauge) were measured. RESULTS: All investigations revealed that bone cements mixed in special mixing kits were characterized by a random level of aeration which led to varied flow parameters. Comparing the experimental results with the theoretical values of the continuity equation, the highest difference in the flow duration reached 140%. DISCUSSION: It has been proven that the aeration of the bone cement alters the flow dynamics. Therefore, much more data are required for statistical analysis to validate a mathematical model of the bone cement flow.

Flow visualization through two types of aortic prosthetic heart valves using stereoscopic high-speed particle image velocimetry.

The number of candidates waiting for a heart valve replacement rises yearly. Even though there is a trend toward implantation of biological valves or reconstruction, the prosthetic heart valves (PHVs) are still commonly used for implantation or as a part of cardiac assist devices in many countries worldwide. However, the hemodynamic consequences of these valves are still not completely understood. Unfortunately, these devices currently do not perform sufficiently on a long-term basis and may lead to several complications, many of them are related to fluid mechanical aspects. A novel method, stereoscopic high-speed particle image velocimetry, was applied to quantify all three velocity components behind a PHV in detailed time domain. In this study, we compared clinically used bileaflet aortic prosthetic (ATS) valve and monoleaflet prototype of tilting disk PHV. The absolute velocities calculated out of two and three velocity components were compared to each other to estimate the overall difference in the desired region of interest. The most significant discrepancies between the two- and three-component absolute velocities were found at the regions of Valsalva sinuses and in a major jet stream of monoleaflet PHV.

Particle image velocimetry tests on pediatric 45-cc and 30-cc ventricle assist devices: effects of heart rate on VAD operation.

BACKGROUND: This study investigated flow analysis inside pediatric ventricle assist devices (VADs) designed and manufactured at the Foundation for Cardiac Surgery Development (FRK), Zabrze, Poland. The main goal of the experiment was to define the minimal heart rate admissible in clinical practice. METHODS: The flow was directed by mechanical, single-disc valves developed at the Lodz University of Technology, Institute of Turbomachinery in Lodz, Poland. VAD operation conditions under different heart rates were analyzed. Measurements were performed on Religa PED pediatric VADs (45 cm3 and 30 cm3) with a particle image velocimetry (PIV) system. RESULTS: Due to the PIV method used, the measurements were made without interference of the measuring system onto the flow structure in the investigated channel, as the measurement procedure is noninvasive. During the investigations conducted in different measurement planes, the majority of the flow volume in the chamber was observable. CONCLUSIONS: The measurements at different heart rates demonstrated a significant influence of this parameter on the flow nature in the heart ventricle. Additionally, it was found that the heart rate affected the operation of heart valves in the VAD.

Numerical investigations of the unsteady blood flow in the end-to-side arteriovenous fistula for hemodialysis.

PURPOSE: The aim of this study was to investigate the blood flow in the end-to-side arteriovenous (a-v) fistula, taking into account its pulsating nature and the patient-specific geometry of blood vessels. Computational Fluid Dynamics (CFD) methods were used for this analysis. METHODS: DICOM images of the fistula, obtained from the angio-computed tomography, were a source of the data applied to develop a 3D geometrical model of the fistula. The model was meshed, then the ANSYS CFX v. 15.0 code was used to perform simulations of the flow in the vessels under analysis. Mesh independence tests were conducted. The non-Newtonian rheological model of blood and the Shear Stress Transport model of turbulence were employed. Blood vessel walls were assumed to be rigid. RESULTS: Flow patterns, velocity fields, the volume flow rate, the wall shear stress (WSS) propagation on particular blood vessel walls were shown versus time. The maximal value of the blood velocity was identified in the anastomosis - the place where the artery is connected to the vein. The flow rate was calculated for all veins receiving blood. CONCLUSIONS: The blood flow in the geometrically complicated a-v fistula was simulated. The values and oscillations of the WSS are the largest in the anastomosis, much lower in the artery and the lowest in the cephalic vein. A strong influence of the mesh on the results concerning the maximal and area-averaged WSS was shown. The relation between simulations of the pulsating and stationary flow under time-averaged flow conditions was presented.

Simulations of the blood flow in the arterio-venous fistula for haemodialysis.

The Ciminio-Brescia arterio-venous fistula is a preferred vascular access for haemodialysis, but it is often associated with the development of vascular complications, due to changes in hemodynamic conditions. Computational fluid dynamics methods were involved to carry out seven simulations of the blood flow through the fistula for the patient specific (geometrical) case and various boundary conditions. The geometrical data, obtained from the angio-computed tomography, were used to create a 3-dimensional CAD model of the fistula. The blood flow patterns, blood velocity and the wall shear stress, thought to play a key role in the development of typical complications (stenoses, thromboses, aneurysms, etc.), have been analyzed in this study. The blood flow is reversed locally downstream the anastomosis (where the artery is connected to the vein) and downstream the stenosis in the cannulated vein. Blood velocity reaches abnormal value in the anastomosis during the systolic phase of the cardiac cycle (2.66 m/s). The wall shear stress changes in this place during a single cycle of the heart operation from 27.9 to 71.3 Pa (average 41.5 Pa). The results are compared with data found in the literature.

Numerical simulations of the blood flow in the patient-specific arterial cerebral circle region.

The Cerebral Circle Region, also known as the Circle of Willis (CoW), is a loop of arteries that form arterial connections between supply arteries to distribute blood throughout the cerebral mass. Among the population, only 25% to 50% have a complete system of arteries forming the CoW. 3D time-varying simulations for three different patient-specific artery anatomies of CoW were performed in order to gain a better insight into the phenomena existing in the cerebral blood flow. The models reconstructed on the basis of computer tomography images start from the aorta and include the largest arteries that supply the CoW and the arteries of CoW. Velocity values measured during the ultrasound examination have been compared with the results of simulations. It is shown that the flow in the right anterior artery in some cases may be supplied from the left internal carotid artery via the anterior communicating artery. The investigations conducted show that the computational fluid dynamic tool, which provides high resolution in both time and space domains, can be used to support physicians in diagnosing patients of different ages and various anatomical arterial structures.

Numerical simulations of the blood flow through vertebral arteries.

Vertebral arteries are two arteries whose structure and location in human body result in development of special flow conditions. For some of the arteries, one can observe a significant difference between flow rates in the left and the right arteries during ultrasonography diagnosis. Usually the reason of such a difference was connected with pathology of the artery in which a smaller flow rate was detected. Simulations of the flow through the selected type of the vertebral artery geometry for twenty five cases of artery diameters have been carried out. The main aim of the presented experiment was to visualize the flow in the region of vertebral arteries junction in the origin of the basilar artery. It is extremely difficult to examine this part of human circulation system, thus numerical experiments may be helpful in understanding the phenomena occurring when two relatively large arteries join together to form one vessel. The obtained results have shown that an individual configuration and diameters of particular arteries can exert an influence on the flow in them and affect a significant difference between flow rates for vertebral arteries. It has been assumed in the investigations that modelled arteries were absolutely normal, without any pathology. In the numerical experiment, the non-Newtonian model of blood was employed.