Application of Computational Fluid Dynamics in Modeling Blood Flow in Human Thoracic Aorta

The aim of this study was to simulate blood flow in thoracic human aorta and understand the role of flow dynamics in the initialization and localization of atherosclerotic plaque in human thoracic aorta. The blood flow dynamics in idealized and realistic models of human thoracic aorta were numerically simulated in three idealized and two realistic thoracic aorta models. The idealized models of thoracic aorta were reconstructed with measurements available from literature, and the realistic models of thoracic aorta were constructed by image processing Computed Tomographic (CT) images. The CT images were made available by South Karelia Central Hospital in Lappeenranta. The reconstruction of thoracic aorta consisted of operations, such as contrast adjustment, image segmentations, and 3D surface rendering. Additional design operations were performed to make the aorta model compatible for the numerical method based computer code. The image processing and design operations were performed with specialized medical image processing software. Pulsatile pressure and velocity boundary conditions were deployed as inlet boundary conditions. The blood flow was assumed homogeneous and incompressible. The blood was assumed to be a Newtonian fluid. The simulations with idealized models of thoracic aorta were carried out with Finite Element Method based computer code, while the simulations with realistic models of thoracic aorta were carried out with Finite Volume Method based computer code. Simulations were carried out for four cardiac cycles. The distribution of flow, pressure and Wall Shear Stress (WSS) observed during the fourth cardiac cycle were extensively analyzed. The aim of carrying out the simulations with idealized model was to get an estimate of flow dynamics in a realistic aorta model. The motive behind the choice of three aorta models with distinct features was to understand the dependence of flow dynamics on aorta anatomy. Highly disturbed and nonuniform distribution of velocity and WSS was observed in aortic arch, near brachiocephalic, left common artery, and left subclavian artery. On the other hand, the WSS profiles at the roots of branches show significant differences with geometry variation of aorta and branches. The comparison of instantaneous WSS profiles revealed that the model with straight branching arteries had relatively lower WSS compared to that in the aorta model with curved branches. In addition to this, significant differences were observed in the spatial and temporal profiles of WSS, flow, and pressure. The study with idealized model was extended to study blood flow in thoracic aorta under the effects of hypertension and hypotension. One of the idealized aorta models was modified along with the boundary conditions to mimic the thoracic aorta under the effects of hypertension and hypotension. The results of simulations with realistic models extracted from CT scans demonstrated more realistic flow dynamics than that in the idealized models. During systole, the velocity in ascending aorta was skewed towards the outer wall of aortic arch. The flow develops secondary flow patterns as it moves downstream towards aortic arch. Unlike idealized models, the distribution of flow was nonplanar and heavily guided by the artery anatomy. Flow cavitation was observed in the aorta model which was imaged giving longer branches. This could not be properly observed in the model with imaging containing a shorter length for aortic branches. The flow circulation was also observed in the inner wall of the aortic arch. However, during the diastole, the flow profiles were almost flat and regular due the acceleration of flow at the inlet. The flow profiles were weakly turbulent during the flow reversal. The complex flow patterns caused a non-uniform distribution of WSS. High WSS was distributed at the junction of branches and aortic arch. Low WSS was distributed at the proximal part of the junction, while intermedium WSS was distributed in the distal part of the junction. The pulsatile nature of the inflow caused oscillating WSS at the branch entry region and inner curvature of aortic arch. Based on the WSS distribution in the realistic model, one of the aorta models was altered to induce artificial atherosclerotic plaque at the branch entry region and inner curvature of aortic arch. Atherosclerotic plaque causing 50% blockage of lumen was introduced in brachiocephalic artery, common carotid artery, left subclavian artery, and aortic arch. The aim of this part of the study was first to study the effect of stenosis on flow and WSS distribution, understand the effect of shape of atherosclerotic plaque on flow and WSS distribution, and finally to investigate the effect of lumen blockage severity on flow and WSS distributions. The results revealed that the distribution of WSS is significantly affected by plaque with mere 50% stenosis. The asymmetric shape of stenosis causes higher WSS in branching arteries than in the cases with symmetric plaque. The flow dynamics within thoracic aorta models has been extensively studied and reported here. The effects of pressure and arterial anatomy on the flow dynamic were investigated. The distribution of complex flow and WSS is correlated with the localization of atherosclerosis. With the available results we can conclude that the thoracic aorta, with complex anatomy is the most vulnerable artery for the localization and development of atherosclerosis. The flow dynamics and arterial anatomy play a role in the localization of atherosclerosis. The patient specific image based models can be used to diagnose the locations in the aorta vulnerable to the development of arterial diseases such as atherosclerosis.

Predictors of fluid intravasation during operative hysteroscopy: a preplanned prospective observational study with 200 cases

PURPOSE: To verify the predictors of intravasation rate during hysteroscopy.METHODS: Prospective observational study (Canadian Task Force classification II-1). All cases (n=200 women; 22 to 86 years old) were treated in an operating room setting. Considering respective bag overfill to calculate water balance, we tested two multiple linear regression models: one for total intravasation (mL) and the other for absorption rate (mL.min-1). The predictors tested (independent variables) were energy (mono/bipolar), tube patency (with/without tubal ligation), hysterometry (cm), age≤50 years, body surface area (m2), surgical complexity (with/without myomectomy) and duration (min).RESULTS: Mean intravasation was significantly higher when myomectomy was performed (442±616 versus 223±332 mL; p<0.01). In the proposed multiple linear regression models for total intravasation (adjusted R2=0.44; p<0.01), the only significant predictors were myomectomy and duration (p<0.01).In the proposed model for intravasation rate (R2=0.39; p<0.01), only myomectomy and hysterometry were significant predictors (p=0.02 and p<0.01, respectively).CONCLUSIONS: Not only myomectomy but also hysterometry were significant predictors of intravasation rate during operative hysteroscopy.

Peritoneal fluid changes in horses subjected to small colon distension

Intestinal devitalization in cases of small colon obstruction may be difficult to detect based only in clinical signs. The purpose was to serially evaluate blood and peritoneal fluid of horses subjected to small colon distension. Seventeen adult horses were allotted in three groups. In the small colon-distended group (DG, n=7) a surgically-implanted latex balloon was inflated to promote intraluminal small colon distension. In the shamoperated group (SG, n=5), the balloon was implanted but not inflated, and no surgery was done in the control group (CG, n=5). Blood and peritoneal fluid were sampled before and after (6 samples with a 30-minute interval) intestinal obstruction for cytological and biochemical analyses. No significant changes in clinical signs occurred within groups or across time during the experimental period. There were no statistical differences among SG and SG groups in hematologic and blood chemistry variables. Although total protein concentration and lactate dehydrogenase (LDH) activity in peritoneal fluid remained most of the time within reference values during the experimental period in all groups, increases from baseline values were detected in SG and DG groups. Such increases occurred earlier, progressively and with greater magnitude in the DG when compared with the SG (P<0.05). Increases from baselines values were also observed in total nucleated cells and neutrophils counts in the DG (P<0.05). In conclusion, distension of the equine small colon induced progressive subtle increases in total protein and LDH concentrations in the peritoneal fluid during the first hours. Serial evaluation of these variables in peritoneal fluid may be useful for early detection of intestinal devitalization in clinical cases of equine small colon obstruction.

Enteral fluid therapy through nasogastric tube in rumen cannulated goats

This paper reports the effects of fluid therapy in goats through nasogastric route with an electrolyte solution composed by concentrations of sodium, potassium and chloride similar to goat plasma (140mmol/L of Na+, 4.5mmol/L of K+, 110mmol/L of Cl-). Four Alpine Chamoisee goats, two of them with evident leakage of the rumen cannulas, were used in a crossover experimental design of two periods and two groups. In one group the two goats were submitted to a treatment protocol to induce dehydration before the fluid therapy, whereas the other group was not. Fluid therapy consisted supplying 10mL/kg/h of the electrolyte solution during 8 hours. No signs of discomfort or stress were observed. The dehydration model employed caused a mild dehydration indicated by decrease in feces humidity, body weight and abdominal circumference, and increase in plasma total solids concentration. During fluid therapy globular volume and plasma total solids decreased, whereas % body weight and abdominal circumference increased. No signs of hyperhydration were observed and serum electrolytes (Na+, Cl-, K+) presented no significant alterations in both groups. Fluid therapy proposed in this study was efficient to treat dehydration, even for rumen cannulated animals with evident leakage, and can be administrated safely with no electrolyte imbalance.

The influence of selected pathological states on the somatostatin-like immunoreactive (SOM-LI) endocrine cells in the mucosal layer of the porcine descending colon

This study reports on changes in the number of somatostatin-like immunoreactive (SOM-LI) endocrine cells in the porcine descending colon, caused by chemically driven inflammation, axotomy and proliferative enteropathy (PE). The distribution pattern of SOM-LI endocrine cells has been studied using the routine single-labelling immunofluorescence technique. Semi-quantitative evaluation of the number of the SOM-immunostained endocrine cells within the mucosal layer of the porcine descending colon has been based on counting of all endocrine cells immunoreactive to SOM per unit area (0,1 mm²). Under physiological conditions the number of SOM-LI endocrine cells has been shown to constitute 3,30±0,22. All applied pathological processes resulted in changes in the SOM-like immunoreactivity, which varied in particular processes studied. The number of SOM-LI endocrine cells increased to 6,28±0,31 and 4,43±0,35 during chemically driven inflammation and proliferative enteropathy, respectively, and decreased to 1,17%±0,16 after axotomy. The obtained results suggest that SOM-LI endocrine cells may participate in various pathological states within porcine descending colon and their functions probably depend on the type of pathological factor.

Outbreak of infectious laryngotracheitis in large multi-age egg layer chicken flocks in Minas Gerais, Brazil

A recent (November 2010) outbreak of infectious laryngotracheitis (ILT) in a multi-age laying hen facility in Minas Gerais state, Brazil, is described. Previous ILT outbreak in laying hens was only notified in São Paulo state, Brazil, in 2002. In the outbreak described here, the affected population was approximately eight million hens, with flock sizes ranging from 100,000 to 2,900,000 chickens. The average mortality ranged from 1 to 6%, and morbidity was around 90% (most of the twenty seven farms of the area were positive for ILT virus). Three multi-age laying farms from one company were selected for this report. Clinical signs included prostration, dyspnea, conjunctivitis, occasional swelling of the paranasal sinuses and bloody mucous nasal discharge. Severely affected chickens presented with dyspnea, gasping and became cyanotic before death. At necropsy, these chickens had fibrinous exudate blocking the larynx and the lumen of cranial part of the trachea. In addition, conjunctivitis with intense hyperemia, edema and sinuses with caseous exudate were present. On histopathology, there were marked necrosis and desquamation of respiratory ephitelium and conjunctiva with numerous syncytial cells formation and fibrinous exudate. Moderate to marked non suppurative (especially lymphocytes and plasma cells) infiltration in the lamina propria also was observed. Sixteen out of 20 examined chickens, eosinophilic intranuclear inclusion bodies were observed in the syncytial cells. The DNA extracted from larynx and trachea produced positive PCR results for ILT virus (ILTV) DNA using formalin-fixed, paraffin embedded (FFPE) samples. Amplicons from a small region of ICP4 gene were submitted to sequencing and showed 100% identity with ILTV EU104910.1 (USA strain), 99% with ILTV JN596963.1 (Australian strain) and 91% with ILTV JN580316.1 (Gallid herpesvirus 1 CEO vaccine strain) and JN580315.1 (Gallid herpesvirus 1 TCO vaccine strain).

Porcine follicular fluid concentration of free insulin-like growth factor-I collected from different diameter ovarian follicles

The study aimed to quantify the concentrations of free IGF-I in serum and fluid of ovarian follicles in pre-pubertal gilts and describe the ovarian morphology by measuring the size of the ovaries and counting the number of surface follicles. Ovaries (n=1,000) from pre-pubertal gilts were obtained immediately after slaughter. A total of 10 samplings were performed, with ovaries obtained from 50 females for each collection. The follicles situated on the surface of each ovary were classified as small (SFs, 2 to 5mm in diameter) or large (LFs 6 to 10mm in diameter) and the follicular fluid was obtained by follicle aspiration. The collection of serum samples was performed after the gilts exsanguination using sterile tubes. From the pool of serum and follicular fluid obtained from 50 females, the concentration of free IGF-I was determined in each sample using an enzyme immunoassay kit (ELISA). The description of ovarian morphometry was performed in 100 ovaries from randomly selected gilts. The larger and smaller lengths of ovaries were measured, and the total number of SFs and LFs present on the surface of each ovary were also counted. The IGF-I concentration was greater (P<0.05) in LFs (170.92±88.29 ng/mL) compared with SFs (67.39±49.90ng/mL) and serum (73.48±34.63ng/mL). The largest and smallest length of the ovaries was 26.0±3.0 and 19.0mm ±2.0mm, respectively. The number of SFs (70.86±25.76) was greater (P<0.01) than LFs (6.54±5.26). The study concluded that LFs present greater levels of IGF-I when compared with SFs and blood, which is related to increased activity of the LFs and its differentiation to ovulation. In addition, ovaries of pre-pubertal gilts have a higher number of SFs compared to LFs. Therefore, our study demonstrated unique data regarding the physiological concentration of free IGF-I in ovarian follicles, that can be used in future research to evaluate the addition of this hormone in the in vitro production media of porcine embryos with the goal to improve the technique efficiency.

Lactated Ringer's solution or 0.9% sodium chloride as fluid therapy in pigeons (Columba livia) submitted to humerus osteosynthesis

The study aimed to compare the effects of intraosseous infusion of lactated Ringer's and 0.9% sodium chloride solutions on the electrolytes and acid-base balance in pigeons submitted to humerus osteosynthesis. Eighteen pigeons were undergoing to isoflurane anesthesia by an avalvular circuit system. They were randomly assigned into two groups (n=9) receiving lactated Ringer's solution (LR) or 0.9% sodium chloride (SC), in a continuous infusion rate of 20mL/kg/h, by using an intraosseous catheter into the tibiotarsus during 60-minute anesthetic procedure. Heart rate (HR), and respiratory rate (RR) were measured every 10 min. Venous blood samples were collected at 0, 30 and 60 minutes to analyze blood pH, PvCO2, HCO3 -, Na+ and K+. Blood gases and electrolytes showed respiratory acidosis in both groups during induction, under physical restraint. This acidosis was evidenced by a decrease of pH since 0 min, associated with a compensatory response, observed by increasing of HCO3 - concentration, at 30 and 60 min. It was not observed any changes on Na+ and K+ serum concentrations. According to the results, there is no reason for choosing one of the two solutions, and it could be concluded that both fluid therapy solutions do not promote any impact on acid-base balance and electrolyte concentrations in pigeons submitted to humerus osteosynthesis.

Effect of different surface roughnesses on a turbulent boundary layer

The classical treatment of rough wall turbulent boundary layers consists in determining the effect the roughness has on the mean velocity profile. This effect is usually described in terms of the roughness function delta U+. The general implication is that different roughness geometries with the same delta U+ will have similar turbulence characteristics, at least at a sufficient distance from the roughness elements. Measurements over two different surface geometries (a mesh roughness and spanwise circular rods regularly spaced in the streamwise direction) with nominally the same delta U+ indicate significant differences in the Reynolds stresses, especially those involving the wall-normal velocity fluctuation, over the outer region. The differences are such that the Reynolds stress anisotropy is smaller over the mesh roughness than the rod roughness. The Reynolds stress anisotropy is largest for a smooth wall. The small-scale anisotropy and interniittency exhibit much smaller differences when the Taylor microscale Reynolds number and the Kolmogorov-normalized mean shear are nominally the same. There is nonetheless evidence that the small-scale structure over the three-dimensional mesh roughness conforms more closely with isotropy than that over the rod-roughened and smooth walls.

An experimental/numerical study of the turbulent boundary layer development along a surface with a sudden change in roughness

A theory for the description of turbulent boundary layer flows over surfaces with a sudden change in roughness is considered. The theory resorts to the concept of displacement in origin to specify a wall function boundary condition for a kappa-epsilon model. An approximate algebraic expression for the displacement in origin is obtained from the experimental data by using the chart method of Perry and Joubert(J.F.M., vol. 17, pp. 193-122, 1963). This expression is subsequently included in the near wall logarithmic velocity profile, which is then adopted as a boundary condition for a kappa-epsilon modelling of the external flow. The results are compared with the lower atmospheric observations made by Bradley(Q. J. Roy. Meteo. Soc., vol. 94, pp. 361-379, 1968) as well as with velocity profiles extracted from a set of wind tunnel experiments carried out by Avelino et al.( 7th ENCIT, 1998). The measurements are found to be in good agreement with the theoretical computations. The skin-friction coefficient was calculated according to the chart method of Perry and Joubert(J.F.M., vol. 17, pp. 193-122, 1963) and to a balance of the integral momentum equation. In particular, the growth of the internal boundary layer thickness obtained from the numerical simulation is compared with predictions of the experimental data calculated by two methods, the "knee" point method and the "merge" point method.

Fluid dynamic modelling and simulation of circulating fluidized bed reactors: importance of the interface turbulence transfer

The main objective of this work is to analyze the importance of the gas-solid interface transfer of the kinetic energy of the turbulent motion on the accuracy of prediction of the fluid dynamic of Circulating Fluidized Bed (CFB) reactors. CFB reactors are used in a variety of industrial applications related to combustion, incineration and catalytic cracking. In this work a two-dimensional fluid dynamic model for gas-particle flow has been used to compute the porosity, the pressure, and the velocity fields of both phases in 2-D axisymmetrical cylindrical co-ordinates. The fluid dynamic model is based on the two fluid model approach in which both phases are considered to be continuous and fully interpenetrating. CFB processes are essentially turbulent. The model of effective stress on each phase is that of a Newtonian fluid, where the effective gas viscosity was calculated from the standard k-epsilon turbulence model and the transport coefficients of the particulate phase were calculated from the kinetic theory of granular flow (KTGF). This work shows that the turbulence transfer between the phases is very important for a better representation of the fluid dynamics of CFB reactors, especially for systems with internal recirculation and high gradients of particle concentration. Two systems with different characteristics were analyzed. The results were compared with experimental data available in the literature. The results were obtained by using a computer code developed by the authors. The finite volume method with collocated grid, the hybrid interpolation scheme, the false time step strategy and SIMPLEC (Semi-Implicit Method for Pressure Linked Equations - Consistent) algorithm were used to obtain the numerical solution.

Non-Linear Unsteady Aerodynamic Response Approximation Using Multi-Layer Functionals

Non-linear functional representation of the aerodynamic response provides a convenient mathematical model for motion-induced unsteady transonic aerodynamic loads response, that accounts for both complex non-linearities and time-history effects. A recent development, based on functional approximation theory, has established a novel functional form; namely, the multi-layer functional. For a large class of non-linear dynamic systems, such multi-layer functional representations can be realised via finite impulse response (FIR) neural networks. Identification of an appropriate FIR neural network model is facilitated by means of a supervised training process in which a limited sample of system input-output data sets is presented to the temporal neural network. The present work describes a procedure for the systematic identification of parameterised neural network models of motion-induced unsteady transonic aerodynamic loads response. The training process is based on a conventional genetic algorithm to optimise the network architecture, combined with a simplified random search algorithm to update weight and bias values. Application of the scheme to representative transonic aerodynamic loads response data for a bidimensional airfoil executing finite-amplitude motion in transonic flow is used to demonstrate the feasibility of the approach. The approach is shown to furnish a satisfactory generalisation property to different motion histories over a range of Mach numbers in the transonic regime.

Modeling of the Fluid-Structure Interaction in Inelastic Piping Systems

In many engineering applications, compliant piping systems conveying liquids are subjected to inelastic deformations due to severe pressure surges such as plastic tubes in modern water supply transmission lines and metallic pipings in nuclear power plants. In these cases the design of such systems may require an adequate modeling of the interactions between the fluid dynamics and the inelastic structural pipe motions. The reliability of the prediction of fluid-pipe behavior depends mainly on the adequacy of the constitutive equations employed in the analysis. In this paper it is proposed a systematic and general approach to consistently incorporate different kinds of inelastic behaviors of the pipe material in a fluid-structure interaction analysis. The main feature of the constitutive equations considered in this work is that a very simple numerical technique can be used for solving the coupled equations describing the dynamics of the fluid and pipe wall. Numerical examples concerning the analysis of polyethylene and stainless steel pipe networks are presented to illustrate the versatility of the proposed approach.

Wind tunnel simulation of atmospheric boundary layer flows

The present work shows how thick boundary layers can be produced in a short wind tunnel with a view to simulate atmospheric flows. Several types of thickening devices are analysed. The experimental assessment of the devices was conducted by considering integral properties of the flow and the spectra: skin-friction, mean velocity profiles in inner and outer co-ordinates and longitudinal turbulence. Designs based on screens, elliptic wedge generators, and cylindrical rod generators are analysed. The paper describes in detail the experimental arrangement, including the features of the wind tunnel and of the instrumentation. The results are compared with experimental data published by other authors and with naturally developed flows.