A study of the aerodynamic characteristics of captor hoods in local exhaust ventilation systems

The research objectives were:- 1.To review the literature to establish the factors which have traditionally been regarded as most crucial to the design of effectlve exhaust ventilation systems. 2. To design, construct, install and calibrate a wind tunnel. 3. To develop procedures for air velocity measurement followed by a comprehensive programme of aerodvnamic data collection and data analysis for a variety of conditions. The major research findings were:- a) The literature in the subject is inadequate. There is a particular need for a much greater understanding of the aerodynamics of the suction flow field. b) The discrepancies between the experimentally observed centre-line velocities and those predicted by conventional formulae are unacceptably large. c) There was little agreement between theoretically calculated and observed velocities in the suction zone of captor hoods. d) Improved empirical formulae for the prediction of centre-line velocity applicable to the classical geometrically shaped suction openings and the flanged condition could be (and were) derived. Further analysis of data revealed that: - i) Point velocity is directly proportional to the suction. flow rate and the ratio of the point velocity to the average face velocity is constant. ii) Both shape, and size of the suction opening are significant factors as the coordinates of their points govern the extent of the effect of the suction flow field. iii) The hypothetical ellipsoidal potential function and hyperbolic streamlines were found experimentally to be correct. iv) The effect of guide plates depends on the size, shape and the angle of fitting. The effect was to very approximately double the suction velocity but the exact effect is difficult to predict. v) The axially symmetric openings produce practically symmetric flow fields. Similarity of connection pieces between the suction opening and the main duct in each case is essential in order to induce a similar suction flow field. Additionally a pilot study was made in which an artificial extraneous air flow was created, measured and its interaction with the suction flow field measured and represented graphically.

Flow processes in a radiant tube burner: isothermal flow

This paper presents the first part of a study of the combustion processes in an industrial radiant tube burner (RTB). The RTB is used typically in heat-treating furnaces. The work was initiated because of the need for improvements in burner lifetime and performance. The present paper is concerned with the flow of combustion air; a future paper will address the combusting flow. A detailed three-dimensional computational fluid dynamics model of the burner was developed, validated with experimental air flow velocity measurements using a split-film probe. Satisfactory agreement was achieved using the k-e turbulence model. Various features along the air inlet passage were subsequently analysed. The effectiveness of the air recuperator swirler was found to be significantly compromised by the need for a generous assembly tolerance. Also, a substantial circumferential flow maldistribution introduced by the swirler is effectively removed by the positioning of a constriction in the downstream passage.

Flow across microvessel walls through the endothelial surface glycocalyx and the interendothelial cleft

A mathematical model is presented for steady fluid flow across microvessel walls through a serial pathway consisting of the endothelial surface glycocalyx and the intercellular cleft between adjacent endothelial cells, with junction strands and their discontinuous gaps. The three-dimensional flow through the pathway from the vessel lumen to the tissue space has been computed numerically based on a Brinkman equation with appropriate values of the Darcy permeability. The predicted values of the hydraulic conductivity Lp, defined as the ratio of the flow rate per unit surface area of the vessel wall to the pressure drop across it, are close to experimental measurements for rat mesentery microvessels. If the values of the Darcy permeability for the surface glycocalyx are determined based on the regular arrangements of fibres with 6nm radius and 8nm spacing proposed recently from the detailed structural measurements, then the present study suggests that the surface glycocalyx could be much less resistant to flow compared to previous estimates by the one-dimensional flow analyses, and the intercellular cleft could be a major determinant of the hydraulic conductivity of the microvessel wall.

Lattice Boltzmann method for flow and heat transfer in microgeometries

Recent technological developments have made it possible to design various microdevices where fluid flow and heat transfer are involved. For the proper design of such systems, the governing physics needs to be investigated. Due to the difficulty to study complex geometries in micro scales using experimental techniques, computational tools are developed to analyze and simulate flow and heat transfer in microgeometries. However, conventional numerical methods using the Navier-Stokes equations fail to predict some aspects of microflows such as nonlinear pressure distribution, increase mass flow rate, slip flow and temperature jump at the solid boundaries. This necessitates the development of new computational methods which depend on the kinetic theory that are both accurate and computationally efficient. In this study, lattice Boltzmann method (LBM) was used to investigate the flow and heat transfer in micro sized geometries. The LBM depends on the Boltzmann equation which is valid in the whole rarefaction regime that can be observed in micro flows. Results were obtained for isothermal channel flows at Knudsen numbers higher than 0.01 at different pressure ratios. LBM solutions for micro-Couette and micro-Poiseuille flow were found to be in good agreement with the analytical solutions valid in the slip flow regime (0.01 < Kn < 0.1) and direct simulation Monte Carlo solutions that are valid in the transition regime (0.1 < Kn < 10) for pressure distribution and velocity field. The isothermal LBM was further extended to simulate flows including heat transfer. The method was first validated for continuum channel flows with and without constrictions by comparing the thermal LBM results against accurate solutions obtained from analytical equations and finite element method. Finally, the capability of thermal LBM was improved by adding the effect of rarefaction and the method was used to analyze the behavior of gas flow in microchannels. The major finding of this research is that, the newly developed particle-based method described here can be used as an alternative numerical tool in order to study non-continuum effects observed in micro-electro-mechanical-systems (MEMS).

Development of a dynamic headspace concentration technique for the non-contact sampling of human odor samples and the creation of canine training aids

Human scent and human remains detection canines are used to locate living or deceased humans under many circumstances. Human scent canines locate individual humans on the basis of their unique scent profile, while human remains detection canines locate the general scent of decomposing human remains. Scent evidence is often collected by law enforcement agencies using a Scent Transfer Unit, a dynamic headspace concentration device. The goals of this research were to evaluate the STU-100 for the collection of human scent samples, and to apply this method to the collection of living and deceased human samples, and to the creation of canine training aids. The airflow rate and collection material used with the STU-100 were evaluated using a novel scent delivery method. Controlled Odor Mimic Permeation Systems were created containing representative standard compounds delivered at known rates, improving the reproducibility of optimization experiments. Flow rates and collection materials were compared. Higher air flow rates usually yielded significantly less total volatile compounds due to compound breakthrough through the collection material. Collection from polymer and cellulose-based materials demonstrated that the molecular backbone of the material is a factor in the trapping and releasing of compounds. The weave of the material also affects compound collection, as those materials with a tighter weave demonstrated enhanced collection efficiencies. Using the optimized method, volatiles were efficiently collected from living and deceased humans. Replicates of the living human samples showed good reproducibility; however, the odor profiles from individuals were not always distinguishable from one another. Analysis of the human remains samples revealed similarity in the type and ratio of compounds. Two types of prototype training aids were developed utilizing combinations of pure compounds as well as volatiles from actual human samples concentrated onto sorbents, which were subsequently used in field tests. The pseudo scent aids had moderate success in field tests, and the Odor pad aids had significant success. This research demonstrates that the STU-100 is a valuable tool for dog handlers and as a field instrument; however, modifications are warranted in order to improve its performance as a method for instrumental detection.

Lattice Boltzmann Method for Flow and Heat Transfer in Microgeometries

Recent technological developments have made it possible to design various microdevices where fluid flow and heat transfer are involved. For the proper design of such systems, the governing physics needs to be investigated. Due to the difficulty to study complex geometries in micro scales using experimental techniques, computational tools are developed to analyze and simulate flow and heat transfer in microgeometries. However, conventional numerical methods using the Navier-Stokes equations fail to predict some aspects of microflows such as nonlinear pressure distribution, increase mass flow rate, slip flow and temperature jump at the solid boundaries. This necessitates the development of new computational methods which depend on the kinetic theory that are both accurate and computationally efficient. In this study, lattice Boltzmann method (LBM) was used to investigate the flow and heat transfer in micro sized geometries. The LBM depends on the Boltzmann equation which is valid in the whole rarefaction regime that can be observed in micro flows. Results were obtained for isothermal channel flows at Knudsen numbers higher than 0.01 at different pressure ratios. LBM solutions for micro-Couette and micro-Poiseuille flow were found to be in good agreement with the analytical solutions valid in the slip flow regime (0.01 < Kn < 0.1) and direct simulation Monte Carlo solutions that are valid in the transition regime (0.1 < Kn < 10) for pressure distribution and velocity field. The isothermal LBM was further extended to simulate flows including heat transfer. The method was first validated for continuum channel flows with and without constrictions by comparing the thermal LBM results against accurate solutions obtained from analytical equations and finite element method. Finally, the capability of thermal LBM was improved by adding the effect of rarefaction and the method was used to analyze the behavior of gas flow in microchannels. The major finding of this research is that, the newly developed particle-based method described here can be used as an alternative numerical tool in order to study non-continuum effects observed in micro-electro-mechanical-systems (MEMS).

Monitoring of down-hole parameters for early kick detection

The sudden hydrocarbon influx from the formation into the wellbore poses a serious risk to the safety of the well. This sudden influx is termed a kick, which, if not controlled, may lead to a blowout. Therefore, early detection of the kick is crucial to minimize the possibility of a blowout occurrence. There is a high probability of delay in kick detection, apart from other issues when using a kick detection system that is exclusively based on surface monitoring. Down-hole monitoring techniques have a potential to detect a kick at its early stage. Down-hole monitoring could be particularly beneficial when the influx occurs as a result of a lost circulation scenario. In a lost circulation scenario, when the down-hole pressure becomes lower than the formation pore pressure, the formation fluid may starts to enter the wellbore. The lost volume of the drilling fluid is compensated by the formation fluid flowing into the well bore, making it difficult to identify the kick based on pit (mud tank) volume observations at the surface. This experimental study investigates the occurrence of a kick based on relative changes in the mass flow rate, pressure, density, and the conductivity of the fluid in the down-hole. Moreover, the parameters that are most sensitive to formation fluid are identified and a methodology to detect a kick without false alarms is reported. Pressure transmitter, the Coriolis flow and density meter, and the conductivity sensor are employed to observe the deteriorating well conditions in the down-hole. These observations are used to assess the occurrence of a kick and associated blowout risk. Monitoring of multiple down-hole parameters has a potential to improve the accuracy of interpretation related to kick occurrence, reduces the number of false alarms, and provides a broad picture of down-hole conditions. The down-hole monitoring techniques have a potential to reduce the kick detection period. A down-hole assembly of the laboratory scale drilling rig model and kick injection setup were designed, measuring instruments were acquired, a frame was fabricated, and the experimental set-up was assembled and tested. This set-up has the necessary features to evaluate kick events while implementing down-hole monitoring techniques. Various kick events are simulated on the drilling rig model. During the first set of experiments compressed air (which represents the formation fluid) is injected with constant pressure margin. In the second set of experiments the compressed air is injected with another pressure margin. The experiments are repeated with another pump (flow) rate as well. This thesis consists of three main parts. The first part gives the general introduction, motivation, outline of the thesis, and a brief description of influx: its causes, various leading and lagging indicators, and description of the several kick detection systems that are in practice in the industry. The second part describes the design and construction of the laboratory scale down-hole assembly of the drilling rig and kick injection setup, which is used to implement the proposed methodology for early kick detection. The third part discusses the experimental work, describes the methodology for early kick detection, and presents experimental results that show how different influx events affect the mass flow rate, pressure, conductivity, and density of the fluid in the down-hole, and the discussion of the results. The last chapter contains summary of the study and future research.

Time series of air chemistry measurements at Georg-von-Neumayer station, Dronning Maud Land, Antarctica in the years 1982 to 1991

The first Air Chemistry Observatory at the German Antarctic station Georg von Neumayer (GvN) was operated for 10 years from 1982 to 1991. The focus of the established observational programme was on characterizing the physical properties and chemical composition of the aerosol, as well as on monitoring the changing trace gas composition of the background atmosphere, especially concerning greenhouse gases. The observatory was designed by the Institut für Umweltphysik, University of Heidelberg (UHEIIUP). The experiments were installed inside the bivouac lodge, mounted on a sledge and put upon a snow hill to prevent snow accumulation during blizzards. All experiments were under daily control and daily performance protocols were documented. A ventilated stainless steel inlet stack (total height about 3-4 m above the snow surface) with a 50% aerodynamic cut-off diameter around 7-10 µm at wind velocities between 4-10 m/s supplied all experiments with ambient air. Contamination free sampling was realized by several means: (i) The Air Chemistry Observatory was situated in a clean air area about 1500 m south of GvN. Due to the fact that northern wind directions are very rare, contamination from the base can be excluded for most of the time. (ii) The power supply (20 kW) is provided by a cable from the main station, thus no fuel-driven generator is operated in the very vicinity. (iii) Contamination-free sampling is controlled by the permanently recorded wind velocity, wind direction and by condensation particle concentration. Contamination was indicated if one of the following criteria were given: Wind direction within a 330°-30° sector, wind velocity <2.2 m/s or >17.5 m/s, or condensation particle concentrations >2500/cm**3 during summer, >800/cm**3 during spring/autumn and >400/cm**3 during winter. If one or a definable combination of these criteria were given, high volume aerosol sampling and part of the trace gas sampling were interrupted. Starting at 1982 through 1991-01-14 surface ozone was measured with an electrochemical concentration cell (ECC). Surface ozone mixing ratio are given in ppbv = parts per 10**9 by volume. The averaging time corresponds to the given time intervals in the data sheet. The accuracy of the values are better than ±1 ppbv and the detection limit is around 1.0 ppbv. Aerosols were sampled on two Whatman 541 cellulose filters in series and analyzed by ion chromatography at the UHEI-IUP. Generally, the sampling period was seven days but could be up to two weeks on occasion. The air flow was around 100 m**3/h and typically 10000-20000 m**3 of ambient air was forced through the filters for one sample. Concentration values are given in nanogram (ng) per 1 m**3 air at standard pressure and temperature (1013 mbar, 273.16 K). Uncertainties of the values were approximately ±10% to ±15% for the main components MSA, chloride, nitrate, sulfate and sodium, and between ±20% and ±30% for the minor species bromide, ammonium, potassium, magnesium and calcium.

Seawater carbonate chemistry and biological processes during experiments with postlarvae of barnacle of Semibalanus balanoides and Elminius modestus, 2010

Ocean acidification and global warming are occurring concomitantly, yet few studies have investigated how organisms will respond to increases in both temperature and CO2. Intertidal microcosms were used to examine growth, shell mineralogy and survival of two intertidal barnacle post-larvae, Semibalanus balanoides and Elminius modestus, at two temperatures (14 and 19°C) and two CO2 concentrations (380 and 1,000 ppm), fed with a mixed diatom-flagellate diet at 15,000 cells ml-1 with flow rate of 10 ml-1 min-1. Control growth rates, using operculum diameter, were 14 ± 8 µm day-1 and 6 ± 2 µm day-1 for S. balanoides and E. modestus, respectively. Subtle, but significant decreases in E. modestus growth rate were observed in high CO2 but there were no impacts on shell calcium content and survival by either elevated temperature or CO2. S. balanoides exhibited no clear alterations in growth rate but did show a large reduction in shell calcium content and survival under elevated temperature and CO2. These results suggest that a decrease by 0.4 pH(NBS) units alone would not be sufficient to directly impact the survival of barnacles during the first month post-settlement. However, in conjunction with a 4-5°C increase in temperature, it appears that significant changes to the biology of these organisms will ensue.

Seawater carbonate chemistry and growth rate during experiments with dinoflagellates, 2007

The effects of dissolved inorganic carbon (DIC) on the growth of 3 red-tide dinoflagellates (Ceratium lineatum, Heterocapsa triquetra and Prorocentrum minimum) were studied at pH 8.0 and at higher pH levels, depending upon the pH tolerance of the individual species. The higher pH levels chosen for experiments were 8.55 for C. lineatum and 9.2 for the other 2 species. At pH 8.0, which approximates the pH found in the open sea, the maximum growth in all species was maintained until the total DIC concentration was reduced below ~0.4 and 0.2 mM for C. lineatum and the other 2 species, respectively. Growth compensation points (concentration of inorganic carbon needed for maintenance of cells) were reached at ~0.18 and 0.05 mM DIC for C. lineatum and the other 2 species, respectively. At higher pH levels, maximum growth rates were lower compared to growth at pH 8, even at very high DIC concentrations, indicating a direct pH effect on growth. Moreover, the concentration of bio-available inorganic carbon (CO2 + HCO3-) required for maintenance as well as the half-saturation constants were increased considerably at high pH compared to pH 8.0. Experiments with pH-drift were carried out at initial concentrations of 2.4 and 1.2 mM DIC to test whether pH or DIC was the main limiting factor at a natural range of DIC. Independent of the initial DIC concentrations, growth rates were similar in both incubations until pH had increased considerably. The results of this study demonstrated that growth of the 3 species was mainly limited by pH, while inorganic carbon limitation played a minor role only at very high pH levels and low initial DIC concentrations.

Gas, water and volume characteristics during controlled degassing experiments on ODP Leg 164 PCS cores

A pressurized core with CH4 hydrate or dissolved CH4 should evolve gas volumes in a predictable manner as pressure is released over time at isothermal conditions. Incremental gas volumes were collected as pressure was released over time from 29 pressure core sampler (PCS) cores from Sites 994, 995, 996, and 997 on the Blake Ridge. Most of these cores were kept at or near 0ºC with an ice bath, and many of these cores yielded substantial quantities of CH4. Volume-pressure plots were constructed for 20 of these cores. Only five plots conform to expected volume and pressure changes for sediment cores with CH4 hydrate under initial pressure and temperature conditions. However, other evidence suggests that sediment in these five and at least five other PCS cores contained CH4 hydrate before core recovery and gas release. Detection of CH4 hydrate in a pressurized sediment core through volume-pressure relationships is complicated by two factors. First, significant quantities of CH4-poor borehole water fill the PCS and come into contact with the core. This leads to dilution of CH4 concentration in interstitial water and, in many cases, decomposition of CH4 hydrate before a degassing experiment begins. Second, degassing experiments were conducted after the PCS had equilibrated in an ice-water bath (0ºC). This temperature is significantly lower than in situ values in the sediment formation before core recovery. Our results and interpretations for PCS cores collected on Leg 164 imply that pressurized containers formerly used by the Deep Sea Drilling Project (DSDP) and currently used by ODP are not appropriately designed for direct detection of gas hydrate in sediment at in situ conditions through volume-pressure relationships.

Interactive effects of ocean acidification and neighboring corals on the growth of Pocillopora verrucosa

The physical and chemical environment around corals, as well as their physiology, can be affected by interactions with neighboring corals. This study employed small colonies (4 cm diameter) of Pocillopora verrucosa and Acropora hyacinthus configured in spatial arrays at 7 cm/s flow speed to test the hypothesis that ocean acidification (OA) alters interactions among them. Interaction effects were quantified for P. verrucosa using three measures of growth: calcification (i.e., weight), horizontal growth, and vertical growth. The study was carried out in May-June 2014 using corals from 10 m depth on the outer reef of Moorea, French Polynesia. Colonies of P. verrucosa were placed next to conspecifics or heterospecifics (A. hyacinthus) in arrangements of two or four colonies (pairs and aggregates) that were incubated at ambient and high pCO2 (1000 µatm) for 28 days. There was an effect of pCO2, and arrangement type on multivariate growth (utilizing the three measures of growth), but no interaction between the main effects. Conversely, arrangement and pCO2 had an interactive effect on calcification, with an overall 23 % depression at high pCO2 versus ambient pCO2 (i.e., pooled among arrangements). Within arrangements, there was a 34-45 % decrease in calcification for solitary and paired conspecifics, but no effect in conspecific aggregates, heterospecific pairs, or heterospecific aggregates. Horizontal growth was negatively affected by pCO2 and arrangement type, while vertical growth was positively affected by arrangement type. Together, our results show that conspecific aggregations can mitigate the negative effects of OA on calcification of colonies within an aggregation.

A Novel Methodology to Characterize the Constitutive Behaviour of Polyethylene terephthalate for the Stretch Blow Moulding process

The stretch blow moulding (SBM) process is the main method for the mass production of PET containers. And understanding the constitutive behaviour of PET during this process is critical for designing the optimum product and process. However due to its nonlinear viscoelastic behaviour, the behaviour of PET is highly sensitive to its thermomechanical history making the task of modelling its constitutive behaviour complex. This means that the constitutive model will be useful only if it is known to be valid under the actual conditions of interest to the SBM process. The aim of this work was to develop a new material characterization method providing new data for the deformation behaviour of PET relevant to the SBM process. In order to achieve this goal, a reliable and robust characterization method was developed based on an instrumented stretch rod and a digital image correlation system to determine the stress-strain relationship of material in deforming preforms during free stretch-blow tests. The effect of preform temperature and air mass flow rate on the deformation behaviour of PET was also investigated.

Modelling of water removal during a paper vacuum dewatering process using a Level-Set method

Water removal in paper manufacturing is an energy-intensive process. The dewatering process generally consists of four stages of which the first three stages include mechanical water removal through gravity filtration, vacuum dewatering and wet pressing. In the fourth stage, water is removed thermally, which is the most expensive stage in terms of energy use. In order to analyse water removal during a vacuum dewatering process, a numerical model was created by using a Level-Set method. Various different 2D structures of the paper model were created in MATLAB code with randomly positioned circular fibres with identical orientation. The model considers the influence of the forming fabric which supports the paper sheet during the dewatering process, by using volume forces to represent flow resistance in the momentum equation. The models were used to estimate the dry content of the porous structure for various dwell times. The relation between dry content and dwell time was compared to laboratory data for paper sheets with basis weights of 20 and 50 g/m2 exposed to vacuum levels between 20 kPa and 60 kPa. The comparison showed reasonable results for dewatering and air flow rates. The random positioning of the fibres influences the dewatering rate slightly. In order to achieve more accurate comparisons, the random orientation of the fibres needs to be considered, as well as the deformation and displacement of the fibres during the dewatering process.

Computer simulations of soft particles in flow

Understanding the dynamics of blood cells is a crucial element to discover biological mechanisms, to develop new efficient drugs, design sophisticated microfluidic devices, for diagnostics. In this work, we focus on the dynamics of red blood cells in microvascular flow. Microvascular blood flow resistance has a strong impact on cardiovascular function and tissue perfusion. The flow resistance in microcirculation is governed by flow behavior of blood through a complex network of vessels, where the distribution of red blood cells across vessel cross-sections may be significantly distorted at vessel bifurcations and junctions. We investigate the development of blood flow and its resistance starting from a dispersed configuration of red blood cells in simulations for different hematocrits, flow rates, vessel diameters, and aggregation interactions between red blood cells. Initially dispersed red blood cells migrate toward the vessel center leading to the formation of a cell-free layer near the wall and to a decrease of the flow resistance. The development of cell-free layer appears to be nearly universal when scaled with a characteristic shear rate of the flow, which allows an estimation of the length of a vessel required for full flow development, $l_c \approx 25D$, with vessel diameter $D$. Thus, the potential effect of red blood cell dispersion at vessel bifurcations and junctions on the flow resistance may be significant in vessels which are shorter or comparable to the length $l_c$. The presence of aggregation interactions between red blood cells lead in general to a reduction of blood flow resistance. The development of the cell-free layer thickness looks similar for both cases with and without aggregation interactions. Although, attractive interactions result in a larger cell-free layer plateau values. However, because the aggregation forces are short-ranged at high enough shear rates ($\bar{\dot{\gamma}} \gtrsim 50~\text{s}^{-1}$) aggregation of red blood cells does not bring a significant change to the blood flow properties. Also, we develop a simple theoretical model which is able to describe the converged cell-free-layer thickness with respect to flow rate assuming steady-state flow. The model is based on the balance between a lift force on red blood cells due to cell-wall hydrodynamic interactions and shear-induced effective pressure due to cell-cell interactions in flow. We expect that these results can also be used to better understand the flow behavior of other suspensions of deformable particles such as vesicles, capsules, and cells. Finally, we investigate segregation phenomena in blood as a two-component suspension under Poiseuille flow, consisting of red blood cells and target cells. The spatial distribution of particles in blood flow is very important. For example, in case of nanoparticle drug delivery, the particles need to come closer to microvessel walls, in order to adhere and bring the drug to a target position within the microvasculature. Here we consider that segregation can be described as a competition between shear-induced diffusion and the lift force that pushes every soft particle in a flow away from the wall. In order to investigate the segregation, on one hand, we have 2D DPD simulations of red blood cells and target cell of different sizes, on the other hand the Fokker-Planck equation for steady state. For the equation we measure force profile, particle distribution and diffusion constant across the channel. We compare simulation results with those from the Fokker-Planck equation and find a very good correspondence between the two approaches. Moreover, we investigate the diffusion behavior of target particles for different hematocrit values and shear rates. Our simulation results indicate that diffusion constant increases with increasing hematocrit and depends linearly on shear rate. The third part of the study describes development of a simulation model of complex vascular geometries. The development of the model is important to reproduce vascular systems of small pieces of tissues which might be gotten from MRI or microscope images. The simulation model of the complex vascular systems might be divided into three parts: modeling the geometry, developing in- and outflow boundary conditions, and simulation domain decomposition for an efficient computation. We have found that for the in- and outflow boundary conditions it is better to use the SDPD fluid than DPD one because of the density fluctuations along the channel of the latter. During the flow in a straight channel, it is difficult to control the density of the DPD fluid. However, the SDPD fluid has not that shortcoming even in more complex channels with many branches and in- and outflows because the force acting on particles is calculated also depending on the local density of the fluid.