Ecological assessment data of the marine flora and fauna of Point Lookout in 2014 (PLEA)

In 2014, UniDive (The University of Queensland Underwater Club) conducted an ecological assessment of the Point Lookout Dive sites for comparison with similar surveys conducted in 2001. Involvement in the project was voluntary. Members of UniDive who were marine experts conducted training for other club members who had no, or limited, experience in identifying marine organisms and mapping habitats. Since the 2001 detailed baseline study, no similar seasonal survey has been conducted. The 2014 data is particularly important given that numerous changes have taken place in relation to the management of, and potential impacts on, these reef sites. In 2009, Moreton Bay Marine Park was re-zoned, and Flat Rock was converted to a marine national park zone (Green zone) with no fishing or anchoring. In 2012, four permanent moorings were installed at Flat Rock. Additionally, the entire area was exposed to the potential effects of the 2011 and 2013 Queensland floods, including flood plumes which carried large quantities of sediment into Moreton Bay and surrounding waters. The population of South East Queensland has increased from 2.49 million in 2001 to 3.18 million in 2011 (BITRE, 2013). This rapidly expanding coastal population has increased the frequency and intensity of both commercial and recreational activities around Point Lookout dive sites (EPA 2008). Methodology used for the PLEA project was based on the 2001 survey protocols, Reef Check Australia protocols and Coral Watch methods. This hybrid methodology was used to monitor substrate and benthos, invertebrates, fish, and reef health impacts. Additional analyses were conducted with georeferenced photo transects. The PLEA marine surveys were conducted over six weekends in 2014 totaling 535 dives and 376 hours underwater. Two training weekends (February and March) were attended by 44 divers, whilst biological surveys were conducted on seasonal weekends (February, May, July and October). Three reefs were surveyed, with two semi-permanent transects at Flat Rock, two at Shag Rock, and one at Manta Ray Bommie. Each transect was sampled once every survey weekend, with the transect tapes deployed at a depth of 10 m below chart datum. Fish populations were assessed using a visual census along 3 x 20 m transects. Each transect was 5 m wide (2.5 m either side of the transect tape), 5 m high and 20 m in length. Fish families and species were chosen that are commonly targeted by recreational or commercial fishers, or targeted by aquarium collectors, and that were easily identified by their body shape. Rare or otherwise unusual species were also recorded. Target invertebrate populations were assessed using visual census along 3 x 20 m transects. Each transect was 5 m wide (2.5 m either side of the transect tape) and 20 m in length. The diver surveying invertebrates conducted a 'U-shaped' search pattern, covering 2.5 m on either side of the transect tape. Target impacts were assessed using a visual census along the 3 x 20 m transects. Each transect was 5 m wide (2.5 m either side of the transect tape) and 20 m in length. The transect was surveyed via a 'U-shaped' search pattern, covering 2.5 m on either side of the transect tape. Substrate surveys were conducted using the point sampling method, enabling percentage cover of substrate types and benthic organisms to be calculated. The substrate or benthos under the transect line was identified at 0.5m intervals, with a 5m gap between each of the three 20m segments. Categories recorded included various growth forms of hard and soft coral, key species/growth forms of algae, other living organisms (i.e. sponges), recently killed coral, and, non-living substrate types (i.e. bare rock, sand, rubble, silt/clay).

Spatially explicit estimates of length and biomass of Micromesistius poutassou (Blue Whiting) from acoustic and trawl surveys in the North-East Atlantic in May 2010

Acoustic estimates of herring and blue whiting abundance were obtained during the surveys using the Simrad ER60 scientific echosounder. The allocation of NASC-values to herring, blue whiting and other acoustic targets were based on the composition of the trawl catches and the appearance of echo recordings. To estimate the abundance, the allocated NASC -values were averaged for ICES-squares (0.5° latitude by 1° longitude). For each statistical square, the unit area density of fish (rA) in number per square nautical mile (N*nm-2) was calculated using standard equations (Foote et al., 1987; Toresen et al., 1998). To estimate the total abundance of fish, the unit area abundance for each statistical square was multiplied by the number of square nautical miles in each statistical square and then summed for all the statistical squares within defined subareas and over the total area. Biomass estimation was calculated by multiplying abundance in numbers by the average weight of the fish in each statistical square then summing all squares within defined subareas and over the total area. The Norwegian BEAM soft-ware (Totland and Godø 2001) was used to make estimates of total biomass.

Spatially explicit estimates of length and biomass of Clupea harengus (Norwegian spring spawning herring) from acoustic and trawl surveys in the North-East Atlantic in May 2011

Acoustic estimates of herring and blue whiting abundance were obtained during the surveys using the Simrad ER60 scientific echosounder. The allocation of NASC-values to herring, blue whiting and other acoustic targets were based on the composition of the trawl catches and the appearance of echo recordings. To estimate the abundance, the allocated NASC -values were averaged for ICES-squares (0.5° latitude by 1° longitude). For each statistical square, the unit area density of fish (rA) in number per square nautical mile (N*nm-2) was calculated using standard equations (Foote et al., 1987; Toresen et al., 1998). To estimate the total abundance of fish, the unit area abundance for each statistical square was multiplied by the number of square nautical miles in each statistical square and then summed for all the statistical squares within defined subareas and over the total area. Biomass estimation was calculated by multiplying abundance in numbers by the average weight of the fish in each statistical square then summing all squares within defined subareas and over the total area. The Norwegian BEAM soft-ware (Totland and Godø 2001) was used to make estimates of total biomass.

Spatially explicit estimates of length and biomass of Clupea harengus (Norwegian spring spawning herring) from acoustic and trawl surveys in the North-East Atlantic in May 2007

Acoustic estimates of herring and blue whiting abundance were obtained during the surveys using the Simrad ER60 scientific echosounder. The allocation of NASC-values to herring, blue whiting and other acoustic targets were based on the composition of the trawl catches and the appearance of echo recordings. To estimate the abundance, the allocated NASC -values were averaged for ICES-squares (0.5° latitude by 1° longitude). For each statistical square, the unit area density of fish (rA) in number per square nautical mile (N*nm-2) was calculated using standard equations (Foote et al., 1987; Toresen et al., 1998). To estimate the total abundance of fish, the unit area abundance for each statistical square was multiplied by the number of square nautical miles in each statistical square and then summed for all the statistical squares within defined subareas and over the total area. Biomass estimation was calculated by multiplying abundance in numbers by the average weight of the fish in each statistical square then summing all squares within defined subareas and over the total area. The Norwegian BEAM soft-ware (Totland and Godø 2001) was used to make estimates of total biomass.

A regularised continuum damage model based on the mesoscopic scale for soft tissue

Material properties of soft fibrous tissues are highly conditioned by the hierarchical structure of this kind of composites. Collagen based tissues present, at decreasing length scales, a complex framework of fibres, fibrils, tropocollagen molecules and amino-acids. Understanding the mechanical behaviour at nano-scale level is critical to accurately incorporate this structural information in phenomenological damage models. In this work we derive a relationship between the mechanical and geometrical properties of the fibril constituents and the soft tissue material parameters at macroscopic scale. A Hodge–Petruska two-dimensional model has been used to describe the fibrils as staggered arrays of tropocollagen molecules. After a mechanical characterisation of each of the fibril components, two fibril failures modes have been defined related with two planes of weakness. A phenomenological continuous damage model with regularised softening was presented along with meso-structurally based definitions for its material parameters. Finally, numerical analysis at fibril, fibre and tissue levels are presented to show the capabilities of the model

A short segment within the cytoplasmic domain of the neural cell adhesion molecule (N-CAM) is essential for N-CAM-induced NF-κB activity in astrocytes

The neural cell adhesion molecule (N-CAM) is expressed on the surface of astrocytes, where its homophilic binding leads to the activation of the transcription factor NF-κB. Transfection of astrocytes with a construct encompassing the transmembrane region and the cytoplasmic domain of N-CAM (designated Tm-Cyto, amino acids 685–839 in the full-length molecule) inhibited this activation up to 40%, and inhibited N-CAM-induced translocation of NF-κB to the nucleus. N-CAM also activated NF-κB in astrocytes from N-CAM knockout mice, presumably through binding to a heterophile. This activation, however, was not blocked by Tm-Cyto expression, indicating that the inhibitory effect of the Tm-Cyto construct is specific for cell surface N-CAM. Deletions and point mutations of the cytoplasmic portion of the Tm-Cyto construct indicated that the region between amino acids 780 and 800 were essential for inhibitory activity. This region contains four threonines (788, 793, 794, and 797). Mutation to alanine of T788, T794, or T797, but not T793, abolished inhibitory activity, as did mutation of T788 or T797 to aspartic acid. A Tm-Cyto construct with T794 mutated to aspartic acid retained inhibitory activity but did not itself induce a constitutive NF-κB response. This result suggests that phosphorylation of T794 may be necessary but is not the triggering event. Overall, these findings define a short segment of the N-CAM cytoplasmic domain that is critical for N-CAM-induced activation of NF-κB and may be important in other N-CAM-mediated signaling.

Invariant chain made in Escherichia coli has an exposed N-terminal segment that blocks antigen binding to HLA-DR1 and a trimeric C-terminal segment that binds empty HLA-DR1.

Invariant chain (Ii), a membrane glycoprotein, binds class II major histocompatibility complex (MHC) glycoproteins, probably via its class II-associated Ii peptide (CLIP) segment, and escorts them toward antigen-containing endosomal compartments. We find that a soluble, trimeric ectodomain of Ii expressed and purified from Escherichia coli blocks peptide binding to soluble HLA-DR1. Proteolysis indicates that Ii contains two structural domains. The C-terminal two-thirds forms an alpha-helical domain that trimerizes and interacts with empty HLA-DR1 molecules, augmenting rather than blocking peptide binding. The N-terminal one-third, which inhibits peptide binding, is proteolytically susceptible over its entire length. In the trimer, the N-terminal domains act independently with each CLIP segment exposed and free to bind an MHC class II molecule, while the C-terminal domains act as a trimeric unit.

Relationship among Corneal Biomechanics, Refractive Error, and Axial Length

Purpose: To evaluate the relationship between different ocular and corneal biomechanical parameters in emmetropic and ametropic healthy white children. Methods: This study included 293 eyes of 293 healthy Spanish children (135 boys and 158 girls), ranging in age from 6 to 17 years. Subjects were divided according to the refractive error: control (emmetropia, 99 children), myopia (100 children), and hyperopia (94 children) groups. In all cases, corneal hysteresis (CH) and corneal resistance factor (CRF) were evaluated with the Ocular Response Analyzer system. Axial length (AL) and mean corneal power were also measured by partial coherence interferometry (IOLMaster), and central corneal thickness (CCT) and anterior chamber depth were measured by anterior segment optical coherence tomography (Visante). Results: Mean (±SD) CH and CRF were 12.12 (±1.71) and 12.30 (±1.89) mm Hg, respectively. Mean (±SD) CCT was 542.68 (±37.20) μm and mean (±SD) spherical equivalent was +0.14 (±3.41) diopters. A positive correlation was found between CH and CRF (p < 0.001), and both correlated as well with CCT (p < 0.0001). Corneal resistance factor was found to decrease with increasing age (p = 0.01). Lower levels of CH were associated with longer AL and more myopia (p < 0.001 and p = 0.001, respectively). Higher values of CH were associated with increasing hyperopia. Significant differences in CH were found between emmetropic and myopic groups (p < 0.001) and between myopic and hyperopic groups (p = 0.011). There were also significant differences in CRF between emmetropic and myopic groups (p = 0.02). Multiple linear regression analysis showed that lower CH and CRF significantly associated with thinner CCT, longer AL, and flatter corneal curvature. Conclusions: The Ocular Response Analyzer corneal biomechanical properties seem to be compromised in myopia from an early age, especially in high myopia.

Spatially explicit estimates of length and biomass of Micromesistius poutassou (Blue Whiting) from acoustic and trawl surveys in the North-East Atlantic in May 2012

Acoustic estimates of herring and blue whiting abundance were obtained during the surveys using the Simrad ER60 scientific echosounder. The allocation of NASC-values to herring, blue whiting and other acoustic targets were based on the composition of the trawl catches and the appearance of echo recordings. To estimate the abundance, the allocated NASC -values were averaged for ICES-squares (0.5° latitude by 1° longitude). For each statistical square, the unit area density of fish (rA) in number per square nautical mile (N*nm-2) was calculated using standard equations (Foote et al., 1987; Toresen et al., 1998). To estimate the total abundance of fish, the unit area abundance for each statistical square was multiplied by the number of square nautical miles in each statistical square and then summed for all the statistical squares within defined subareas and over the total area. Biomass estimation was calculated by multiplying abundance in numbers by the average weight of the fish in each statistical square then summing all squares within defined subareas and over the total area. The Norwegian BEAM soft-ware (Totland and Godø 2001) was used to make estimates of total biomass.

Biomechanical aspects of the anterior segment in human myopia

The thesis investigates the relationship between the biomechanical properties of the anterior human sclera and cornea in vivo using Schiotz tonometry (ST), rebound tonometry (RBT, iCare) and the Ocular Response Analyser (ORA, Reichert). Significant differences in properties were found to occur between scleral quadrants. Structural correlates for the differences were examined using Partial Coherent Interferometry (IOLMaster, Zeiss), Optical Coherent tomography (Visante OCT), rotating Scheimpflug photography (Pentacam, Oculus) and 3-D Magnetic Resonance Imaging (MRI). Subject groups were employed that allowed investigation of variation pertaining to ethnicity and refractive error. One hundred thirty-five young adult subjects were drawn from three ethnic groups: British-White (BW), British-South-Asian (BSA) and Hong-Kong-Chinese (HKC) comprising non-myopes and myopes. Principal observations: ST demonstrated significant regional variation in scleral resistance a) with lowest levels at quadrant superior-temporal and highest at inferior-nasal; b) with distance from the limbus, anterior locations showing greater resistance. Variations in resistance using RBT were similar to those found with ST; however the predominantly myopic HKC group had a greater overall mean resistance when compared to the BW-BSA group. OCT-derived scleral thickness measurements indicated the sclera to be thinner superiorly than inferiorly. Thickness varied with distance from the corneolimbal junction, with a decline from 1 to 2 mm followed by a successive increase from 3 to 7 mm. ORA data varied with ethnicity and refractive status; whilst axial length (AL) was associated with corneal biometrics for BW-BSA individuals it was associated with IOP in the HKC individuals. Complex interrelationships were found between ORA Additional-Waveform-Parameters and biometric data provided by the Pentacam. OCT indicated ciliary muscle thickness to be greater in myopia and more directly linked to posterior ocular volume (from MRI) than AL. Temporal surface areas (SAs, from MRI) were significantly smaller than nasal SAs in myopic eyes; globe bulbosity (from MRI) was constant across quadrants.

Use of imaging technology to better understand soft contact lens fit dynamics

The principal theme of this thesis is the identification of additional factors affecting, and consequently to better allow, the prediction of soft contact lens fit. Various models have been put forward in an attempt to predict the parameters that influence soft contact lens fit dynamics; however, the factors that influence variation in soft lens fit are still not fully understood. The investigations in this body of work involved the use of a variety of different imaging techniques to both quantify the anterior ocular topography and assess lens fit. The use of Anterior-Segment Optical Coherence Tomography (AS-OCT) allowed for a more complete characterisation of the cornea and corneoscleral profile (CSP) than either conventional keratometry or videokeratoscopy alone, and for the collection of normative data relating to the CSP for a substantial sample size. The scleral face was identified as being rotationally asymmetric, the mean corneoscleral junction (CSJ) angle being sharpest nasally and becoming progressively flatter at the temporal, inferior and superior limbal junctions. Additionally, 77% of all CSJ angles were within ±50 of 1800, demonstrating an almost tangential extension of the cornea to form the paralimbal sclera. Use of AS-OCT allowed for a more robust determination of corneal diameter than that of white-to-white (WTW) measurement, which is highly variable and dependent on changes in peripheral corneal transparency. Significant differences in ocular topography were found between different ethnicities and sexes, most notably for corneal diameter and corneal sagittal height variables. Lens tightness was found to be significantly correlated with the difference between horizontal CSJ angles (r =+0.40, P =0.0086). Modelling of the CSP data gained allowed for prediction of up to 24% of the variance in contact lens fit; however, it was likely that stronger associations and an increase in the modelled prediction of variance in fit may have occurred had an objective method of lens fit assessment have been made. A subsequent investigation to determine the validity and repeatability of objective contact lens fit assessment using digital video capture showed no significant benefit over subjective evaluation. The technique, however, was employed in the ensuing investigation to show significant changes in lens fit between 8 hours (the longest duration of wear previously examined) and 16 hours, demonstrating that wearing time is an additional factor driving lens fit dynamics. The modelling of data from enhanced videokeratoscopy composite maps alone allowed for up to 77% of the variance in soft contact lens fit, and up to almost 90% to be predicted when used in conjunction with OCT. The investigations provided further insight into the ocular topography and factors affecting soft contact lens fit.

Measurement of scleral thickness in humans using anterior segment optical coherent tomography

Anterior segment optical coherent tomography (AS-OCT, Visante; Zeiss) is used to examine meridional variation in anterior scleral thickness (AST) and its association with refractive error, ethnicity and gender. Scleral cross-sections of 74 individuals (28 males; 46 females; aged between 18-40 years (27.7±5.3)) were sampled twice in random order in 8 meridians: [superior (S), inferior (I), nasal (N), temporal (T), superior-temporal (ST), superior-nasal (SN), inferior-temporal (IT) and inferior-nasal (IN)]. AST was measured in 1mm anterior-toposterior increments (designated the A-P distance) from the scleral spur (SS) over a 6mm distance. Axial length and refractive error were measured with a Zeiss IOLMaster biometer and an open-view binocular Shin-Nippon autorefractor. Intra- And inter-observer variability of AST was assessed for each of the 8 meridians. Mixed repeated measures ANOVAs tested meridional and A-P distance differences in AST with refractive error, gender and ethnicity. Only right eye data were analysed. AST (mean±SD) across all meridians and A-P distances was 725±46μm. Meridian SN was the thinnest (662±57μm) and I the thickest (806 ±60μm). Significant differences were found between all meridians (p<0.001), except S:ST, IT:IN, IT:N and IN:N. Significant differences between A-P distances were found except between SS and 6 mm and between 2 and 4mm. AST measurements at 1mm (682±48 μm) were the thinnest and at 6mm (818±49 μm) the thickest (p<0.001); a significant interaction occurred between meridians and A-P distances (p<0.001). AST was significantly greater (p<0.001) in male subjects but no significant differences were found between refractive error or ethnicity. Significant variations in AST occur with regard to meridian and distance from the SS and may have utility in selecting optimum sites for pharmaceutical or surgical intervention.

Evolução dos sistemas de falhas de borda do Rifte Potiguar com base em curvas de crescimento de falhas

The discussion about rift evolution in the Brazilian Equatorial margin during the South America-Africa breakup in the Jurassic/Cretaceous has been focused in many researches. But rift evolution based on development and growth of faults has not been well explored. In this sense, we investigated the Cretaceous Potiguar Basin in the Equatorial margin of Brazil to understand the geometry of major faults and the influence of crustal heterogeneity and preexisting structural fabric in the evolution of the basin internal architecture. Previous studies pointed out that the rift is an asymmetrical half-graben elongated along the NE-SW direction. We used 2D seismic, well logs and 3D gravity modeling to analyze four major border fault segments and determine their maximum displacement (Dmax) and length (L) ratio in the Potiguar Rift. We constrained the 3D gravity modeling with well data and the interpretation of seismic sections. The difference of the fault displacement measured in the gravity model is in the order of 10% compared to seismic and well data. The fault-growth curves allowed us to divide the faulted rift border into four main fault segments, which provide roughly similar Dmax/L ratios. Fault-growth curves suggest that a regional uniform tectonic mechanism influenced growth of the rift fault segments. The variation of the displacements along the fault segments indicates that the fault segments were formed independently during rift initiation and were linked by hard and soft linkages. The latter formed relay ramps. In the interconnection zones the Dmax/L ratios are highest due to interference of fault segment motions. We divided the evolution of the Potiguar Rift into five stages based on these ratios and correlated them with the major tectonic stages of the breakup between South America and Africa in Early Cretaceous.

Evolução dos sistemas de falhas de borda do Rifte Potiguar com base em curvas de crescimento de falhas

The discussion about rift evolution in the Brazilian Equatorial margin during the South America-Africa breakup in the Jurassic/Cretaceous has been focused in many researches. But rift evolution based on development and growth of faults has not been well explored. In this sense, we investigated the Cretaceous Potiguar Basin in the Equatorial margin of Brazil to understand the geometry of major faults and the influence of crustal heterogeneity and preexisting structural fabric in the evolution of the basin internal architecture. Previous studies pointed out that the rift is an asymmetrical half-graben elongated along the NE-SW direction. We used 2D seismic, well logs and 3D gravity modeling to analyze four major border fault segments and determine their maximum displacement (Dmax) and length (L) ratio in the Potiguar Rift. We constrained the 3D gravity modeling with well data and the interpretation of seismic sections. The difference of the fault displacement measured in the gravity model is in the order of 10% compared to seismic and well data. The fault-growth curves allowed us to divide the faulted rift border into four main fault segments, which provide roughly similar Dmax/L ratios. Fault-growth curves suggest that a regional uniform tectonic mechanism influenced growth of the rift fault segments. The variation of the displacements along the fault segments indicates that the fault segments were formed independently during rift initiation and were linked by hard and soft linkages. The latter formed relay ramps. In the interconnection zones the Dmax/L ratios are highest due to interference of fault segment motions. We divided the evolution of the Potiguar Rift into five stages based on these ratios and correlated them with the major tectonic stages of the breakup between South America and Africa in Early Cretaceous.

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.