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.

INFLUENTIAL FACTORS IN SURFACE WATER QUALITY IN CATCHMENTS WITHIN THE PAMPA BIOME WITH DIFFERENT LAND USE1

The aim of this study was to identify, by multivariate statistical technique, the physic, chemical and biological variables that best characterize the quality of surface waters in two small rural catchments with different land uses (eucalyptus silviculture (SC) vs. pasture and extensive livestock (LC)) located in Rosário do Sul, RS - Brazil. Monitoring was conducted during the months of August 2011 to August 2012 and the following parameters were analyzed: Ca2+, Mg2+, K+, SO42-, Cl-, pH, electrical conductivity, turbidity, alkalinity, suspended and dissolved solids, biochemical oxygen demand , total coliforms, Escherichia coli and temperature, flow and rainfall. Through the use of FA/PCA, it was found that the model best fit to express water quality of in LC that was composed of five factors which account for 83.5% of the total variance, while for SC, four factors accounted for 85.12% of the variance. In LC, the five main factors were, respectively, soluble salts, diffuse pollution, solid, and both anthropogenic and organic factors. In SC, the four factors were namely: soluble salts, mineral, nutritional and diffuse pollution factors. The results of this study showed that by replacing the traditional soil usage (pasture and livestock) with planted forest, diffuse pollution was attenuated but, however, it did not result in major changes in the physical-chemical and biological characteristics of the water. Another point to note is that factorial analysis did not result in a large reduction in the number of variables, once the best model fit occurred with the addition of 15 of 18 analyzed variables (LC) and 17 of 18 analyzed variables (SC).

Seasonality of net radiation in two sub-basins of Paracatu by the use of modis sensor products

The net radiation (Rn) represents the main source of energy for physical and chemical processes that occur in the surface-atmosphere interface, and it is used for air and soil heating, water transfer, in the form of vapor from the surface to the atmosphere, and for the metabolism of plants, especially photosynthesis. If there is no record of net radiation in certain areas, the use of information is important to help determine it. Among them we can highlight those provided by remote sensing. In this context, this work aims to estimate the net radiation, with the use of products of MODIS sensor, in the sub-basins of Entre Ribeiros creek and Preto River, located between the Brazilian states of Goiás and Minas Gerais. The SEBAL (Surface Energy Balance Algorithm for Land) was used to obtain the Rn in four different days in the period of July to October, 2007. The Rn results obtained were consistent with others cited in the literature and are important because the orbital information can help determine the Rn in areas where there are not automatic weather stations to record the net radiation.

Horizontal pressures in cylindrical metal silos and comparison with different international standards

The focus of this research was to evaluate the horizontal pressures on a cylindrical metal silo of corrugated walls and flat bottom with 1.82m diameter and 5.4m high, and to compare the values with those obtained theoretically by the ISO 11697, EP 433 and AS 3774 standards. The silo was symmetrically filled and constant speed with wheat cv. soft red for two different height/diameter ratios (H/D) and was unloaded through three orifices with a diameter of 71.6mm, one concentric and two eccentrics. Horizontal pressures were measured on the walls of the silo at three positions using hydraulic type pressure cells. The results showed that shortly after the start of the unloading, there was a mass flow above the quota of H/D = 1.2, whereas below this quota funnel flow occurred. It can be said that the EP 433 standard was more appropriate to predict horizontal pressures in silos in H/D ratio = 1.0, with eccentric unloading. For the H/D ratio = 1.5, AS 3774 standard was the one that produced values closer to the experimental.

Application of response surface methodology to study the biological removal of nitrogen from effluent of cattle slaughterhouse in a sequencing batch reactor

The aim of this study was to evaluate the efficiency of a sequencing batch reactor (SBR) on biological removal of nitrogen from cattle slaughterhouse wastewater by nitrification/denitrification processes. The effects of initial concentration of ammoniacal nitrogen were investigated at 100; 150 and 200 mg L-1 and air flow rate at 0.125; 0.375 and 0.625 L min¹ Lreactor-1 on the nitrogen compounds removal, by a Central Composite Rotational Design (CCRD) configuration. There were variations from 9.2 to 94.9%, 4.0 to 19.6% and 20.8 to 92.0% in the conversion of ammoniacal nitrogen to nitrate and nitrite concentration and removal of total nitrogen, respectively. The increase of air flow rate and decrease of the initial concentration of ammoniacal nitrogen resulted in higher efficiencies of total nitrogen removal, as well as the conversion of ammoniacal nitrogen to nitrate. During the pre-established intervals of this study, the removal and conversion efficiencies of nitrogen compounds above 85% were achieved in air flow rate variations from 0.375 to 0.725 L min-1 Lreactor-1 and initial concentration of ammoniacal nitrogen from 80 to 200 mg L-1. On denitrification process, we obtained efficiencies from 91.5 to 96.9% on the removal of nitrite/nitrate and from 78.3 to 87.9% on the removal of organic matter.

Impact of the substitution of reference annual streamflow by monthly streamflow on the potential use of water resources

The consideration of the streamflow seasonality has a high potential to improve the water use. In order to give subsidies to the optimization of water use, it was evaluated the impact of the change of reference annual streamflow by the monthly streamflows in the potential water use throughout the hydrography of Paracatu sub-Basin. It was evaluated the impact on Q7,10 (lowest average streamflow during a 7-day period with an average recurrence of 10 years) and on Q95 (permanent flow present 95% of the time). The use of monthly streamflow to substitute the annual streamflow had a high potential of improvement of water resources use in the sub-Basin studied. The use of monthly Q 7,10 in substitution of annual Q 7,10 increases the potential water use that vary from about 10% in the months of lower water availability to values exceeding 200% in the months with higher availability of surface water resources. The use of monthly Q95 in substitution of the annual Q95 implies in changes oscillating from reduction of 37% in months of higher water restriction to values exceeding 100% in the months of higher availability, so the use of monthly Q95 instead of the annual Q95 enables the more rational and safe use of water resources.

Internal ambience of beehives Apis mellifera with different colors and roofing materials in the sub middle of the São Francisco Valley

ABSTRACTThe current study aimed to evaluate the influence of three colors and two types of roofing materials under the internal temperature of bee colonies Apis mellifera. The experiment was conducted at the Agricultural Sciences Campus at the Federal University of Sao Francisco Valley located in Petrolina-PE, in November and December 2013, using 24 colonies housed in Langstroth hives. The experiment was a completely randomized factorial design (3x2) with three colors of box (blue, white, and traditional) and two types of cover (with and without the use of plaster) with six treatments and four replications. The internal temperature dates of the colonies were hourly recorded, during 24 hours, and surface temperatures were hourly recorded between 08h00 and 17h00. The highest values for surface and internal temperature were registered in the blue painted boxes without the use of plasterboard, and the blue painted boxes covered with plasterboard respectively. However, the lowest values ​​were found in the white painted hives and hives that have not received the plasterboard. It is recommended to paint boxes with bright colors, and the use of plasterboard had no effect in reducing the internal temperature.

Backpressure effects on the flow-pressure relation of driplines

ABSTRACT For drip irrigation design and management, it is necessary to know the relation between flow and pressure acting on emitters. In the case of subsurface drip irrigation, the backpressure phenomenon may change the hydraulic characteristics of emitters. Thus, this study aimed at determining such relationship between flow and pressure of different driplines in surface and subsurface conditions; aiming to find possible differences in hydraulic behavior. We tested four emitter types; two pressure compensating (D5000 and Hydro PCND) and two non-pressure compensating (TalDrip and Jardiline). Emitter flow rates were attained in atmospheric conditions and submerged in water, in which submergence levels represented backpressure. Assays were performed using inlet pressures of 80, 100, 120, and 150 kPa for the Hydro PCND dripline and 25, 50, 100, and 150 kPa for the other ones; the backpressures were of 0.49, 1.47, 2.45, 4.41 and 6.37 kPa with four replications. The emitters had their proportionality constants and discharge exponents changed in submerged applications, representing backpressure effect. Non-pressure compensating emitters had their discharge exponent decreased, while in pressure compensating ones, it was increased. Backpressure reduced emitter flow rates at all evaluated pressures.

Modeling of turbulent flow through radial diffuser

The work considers the modeling of turbulent flow in radial diffuser with axial feeding. Due to its claimed capability to predict flow including features such as separation, curvature and adverse pressure gradient, the RNG k-epsilon model of Orzag et al. (1993) is applied in the present analysis. The governing equations are numerically solved using the finite volume methodology. Experiments were conducted to assess the turbulence model. Numerical results of pressure distribution on the front disk surface for different flow conditions when compared to the experimental data indicated that the RNG k-epsilon model is adequate to predict this class of flow.

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.

Comparison between one-dimensional uncoupled and convection-conduction conjugated approaches in finned surface heat transfer

This work studies the forced convection problem in internal flow between concentric annular ducts, with radial fins at the internal tube surface. The finned surface heat transfer is analyzed by two different approaches. In the first one, it is assumed one-dimensional heat conduction along the internal tube wall and fins, with the convection heat transfer coefficient being a known parameter, determined by an uncoupled solution. In the other way, named conjugated approach, the mathematical model (continuity, momentum, energy and K-epsilon equations) applied to tube annuli problem was numerically solved using finite element technique in a coupled formulation. At first time, a comparison was made between results obtained for the conjugated problem and experimental data, showing good agreement. Then, the temperature profiles under these two approaches were compared to each other to analyze the validity of the one-dimensional classical formulation that has been utilized in the heat exchanger design.

Surface Tension Implementation for Gensmac2d

In the present work we describe a method which allows the incorporation of surface tension into the GENSMAC2D code. This is achieved on two scales. First on the scale of a cell, the surface tension effects are incorporated into the free surface boundary conditions through the computation of the capillary pressure. The required curvature is estimated by fitting a least square circle to the free surface using the tracking particles in the cell and in its close neighbors. On a sub-cell scale, short wavelength perturbations are filtered out using a local 4-point stencil which is mass conservative. An efficient implementation is obtained through a dual representation of the cell data, using both a matrix representation, for ease at identifying neighbouring cells, and also a tree data structure, which permits the representation of specific groups of cells with additional information pertaining to that group. The resulting code is shown to be robust, and to produce accurate results when compared with exact solutions of selected fluid dynamic problems involving surface tension.

Horizontal projections of area 17 in Cebus monkeys: metric features, and modular and laminar distribution

Metric features and modular and laminar distributions of intrinsic projections of area 17 were studied in Cebus apella. Anterogradely and retrogradely labeled cell appendages were obtained using both saturated pellets and iontophoretic injections of biocytin into the operculum. Laminar and modular distributions of the labeled processes were analyzed using Nissl counterstaining, and/or cytochrome oxidase and/or NADPH-diaphorase histochemistry. We distinguished three labeled cell types: pyramidal, star pyramidal and stellate cells located in supragranular cortical layers (principally in layers IIIa, IIIb a, IIIb ß and IIIc). Three distinct axon terminal morphologies were found, i.e., Ia, Ib and II located in granular and supragranular layers. Both complete and partial segregation of group I axon terminals relative to the limits of the blobs of V1 were found. The results are compatible with recent evidence of incomplete segregation of visual information flow in V1 of Old and New World primates

Respiratory effects of kynurenic acid microinjected into the ventromedullary surface of the rat

Several studies demonstrate that, within the ventral medullary surface (VMS), excitatory amino acids are necessary components of the neural circuits involved in the tonic and reflex control of respiration and circulation. In the present study we investigated the cardiorespiratory effects of unilateral microinjections of the broad spectrum glutamate antagonist kynurenic acid (2 nmol/200 nl) along the VMS of urethane-anesthetized rats. Within the VMS only one region was responsive to this drug. This area includes most of the intermediate respiratory area, partially overlapping the rostral ventrolateral medulla (IA/RVL). When microinjected into the IA/RVL, kynurenic acid produced a respiratory depression, without changes in mean arterial pressure or heart rate. The respiratory depression observed was characterized by a decrease in ventilation, tidal volume and mean inspiratory flow and an increase in respiratory frequency. Therefore, the observed respiratory depression was entirely due to a reduction in the inspiratory drive. Microinjections of vehicle (200 nl of saline) into this area produced no significant changes in breathing pattern, blood pressure or heart rate. Respiratory depression in response to the blockade of glutamatergic receptors inside the rostral VMS suggests that neurons at this site have an endogenous glutamatergic input controlling the respiratory cycle duration and the inspiratory drive transmission.

Two-phase flow measurements with spacer grid in a rod bundle geometry

Heat transfer effectiveness in nuclear rod bundles is of great importance to nuclear reactor safety and economics. An important design parameter is the Critical Heat Flux (CHF), which limits the transferred heat from the fuel to the coolant. The CHF is determined by flow behaviour, especially the turbulence created inside the fuel rod bundle. Adiabatic experiments can be used to characterize the flow behaviour separately from the heat transfer phenomena in diabatic flow. To enhance the turbulence, mixing vanes are attached to spacer grids, which hold the rods in place. The vanes either make the flow swirl around a single sub-channel or induce cross-mixing between adjacent sub-channels. In adiabatic two-phase conditions an important phenomenon that can be investigated is the effect of the spacer on canceling the lift force, which collects the small bubbles to the rod surfaces leading to decreased CHF in diabatic conditions and thus limits the reactor power. Computational Fluid Dynamics (CFD) can be used to simulate the flow numerically and to test how different spacer configurations affect the flow. Experimental data is needed to validate and verify the used CFD models. Especially the modeling of turbulence is challenging even for single-phase flow inside the complex sub-channel geometry. In two-phase flow other factors such as bubble dynamics further complicate the modeling. To investigate the spacer grid effect on two-phase flow, and to provide further experimental data for CFD validation, a series of experiments was run on an adiabatic sub-channel flow loop using a duct-type spacer grid with different configurations. Utilizing the wire-mesh sensor technology, the facility gives high resolution experimental data in both time and space. The experimental results indicate that the duct-type spacer grid is less effective in canceling the lift force effect than the egg-crate type spacer tested earlier.