The flow field around a freely swimming copepod in steady motion. Part II: Numerical simulation

Three-dimensional, numerical simulations of the flow field arounda freely swimming model-copepod were performed using a finite-volumecode. The model copepod had a realistic body shape representedby a curvilinear body-fitted coordinate system. The beatingmovement of the cephalic appendages was replaced by a distributedforce field acting on the water ventrally adjacent to the copepod'sbody. In the simulations, we took into account that freely swimmingcopepods are self-propelled bodies through properly couplingthe Navier–Stokes equations with the dynamic equationfor the copepod's body. Flow fields were calculated for fivesteady motions: (1) hovering, (2) sinking, (3) upwards swimming,(4) backwards swimming and (5) forwards swimming. The numericalresults confirm the conclusions drawn from the theoretical analysisusing Stokes flow models by Jiang et al. [in a companion paper(Jiang et al., 2002a)] for a spherical copepod shape and showthat the geometry of the flow field around a freely swimmingcopepod varies significantly with the different swimming behaviours.When a copepod hovers in the water, or swims very slowly, itgenerates a cone-shaped and wide flow field. In contrast, whena copepod sinks, or swims fast, the flow geometry is not cone-shaped,but cylindrical, narrow and long. The relationships betweencopepods' swimming behaviour and body orientation, hydrodynamicconspicuousness, energetics as well as feeding efficiency werediscussed, based on the simulation data. It is shown that thebehaviour of hovering or swimming slowly is more energeticallyefficient in terms of relative capture volume per energy expendedthan the behaviour of swimming fast, i.e. for a same amountof energy expended a hovering or slow-swimming copepod is ableto scan more water than a fast-swimming one. The numerical resultsalso suggest that the flow field generated by a fast-swimmingcopepod enables the copepod to use mechanoreception to perceivethe food/prey and therefore increases the food concentrationin the swept volume and that the flow field around a free-sinkingcopepod favours the copepod's mechanoreception while minimizingthe energy expense, so that the energy budget can still be maintainedfor both cases.     

Numerical simulation and comparative analysis of flow field in axial blood pumps

The objective study was to estimate the rheological properties and physiological compatibility of the blood pump by simulating the internal flow field of the blood pump. In this study we use computational fluid dynamics method to simulate and analyse two models of axial blood pumps with a three-blade diffuser and a six-blade diffuser, named pump I and pump II, respectively, and to compare the flow patterns of these two kinds of blood pumps while both of them satisfy the conditions of the normal human blood differential pressure and blood flow. Results indicate that (i) the high shear force occurs between the diffuser and the rotor in which the crucial place leads to haemolysis and (ii) under the condition of 100 mmHg pressure head and 5 l/min flow rate, the difference between the two kinds of blood pumps, as far as the haemolytic performance is concerned, is notable. The haemolysis index of the two pumps is 0.32% and 0.2%. In conclusion, the performance of the blood pump is influenced by the diffusers' blade number. Pump II performed better than pump I, which can be the basic model for blood pump option.     

Numerical flow simulation of pool-type fishways: new ways with well-known tools

Fish passage structures are built to restore the connectivity of rivers and allow the migration of aquatic fauna. In order to assess the functioning of a pool-type fishway, it is necessary, inter alia, to possess detailed knowledge of its flow structure, since observations of fishways, and in particular of the visible water surface, can only provide a rough idea of the actual conditions inside the pools. Numerical simulation has been used for many years to support engineering sciences. Especially, the modeling of flow processes in hydraulic machines can, on the one hand, help avoid major problems during the design stage of fish passage structures and, on the other, improve the structure’s hydraulic performance. To this end, two diploma theses within the framework of a research project of the local energy supplier Energie Baden-Württemberg AG (EnBW) employed modeling tools for 3D flow simulation, primarily for pool-slot fishways (PSF), and for traditional vertical slot fishways (VSF). Guest editors: R. L. Welcomme & G. Marmulla Hydropower, Flood Control and Water Abstraction: Implications for Fish and Fisheries     

Numerical simulation of the remodelling process of trabecular architecture around dental implants

Dental implants may alter the mechanical environment in the jawbone, thereby causing remodelling and adaptation of the surrounding trabecular bone tissues. To improve the efficacy of dental implant systems, it is necessary to consider the effect of bone remodelling on the performance of the prosthetic systems. In this study, finite element simulations were implemented to predict the evolution of microarchitecture around four implant systems using a previously developed model that combines both adaptive and microdamage-based mechano-sensory mechanisms in bone remodelling process. Changes in the trabecular architecture around dental implants were mainly focused. The simulation results indicate that the orientational and ladder-like architecture around the implants predicted herein is in good agreement with those observed in animal experiments and clinical observations. The proposed algorithms were shown to be effective in simulating the remodelling process of trabecular architecture around dental implant systems. In addition, the architectural features around four typical dental implant systems in alveolar bone were evaluated comparatively.     

Bacterial transformation and biodegradation processes simulation in horizontal subsurface flow constructed wetlands using CWM1-RETRASO

The performance and reliability of the CWM1-RETRASO model for simulating processes in horizontal subsurface flow constructed wetlands (HSSF CWs) and the relative contribution of different microbial reactions to organic matter (COD) removal in a HSSF CW treating urban wastewater were evaluated. Various different approaches with diverse influent configurations were simulated. According to the simulations, anaerobic processes were more widespread in the simulated wetland and contributed to a higher COD removal rate [72-79%] than anoxic [0-1%] and aerobic reactions [20-27%] did. In all the cases tested, the reaction that most contributed to COD removal was methanogenesis [58-73%]. All results provided by the model were in consonance with literature and experimental field observations, suggesting a good performance and reliability of CWM1-RETRASO. According to the good simulation predictions, CWM1-RETRASO is the first mechanistic model able to successfully simulate the processes described by the CWM1 model in HSSF CWs.     

Proposal of a thorax segment coordinate system for the 3D kinematical analysis of the cervical spine

The International Society of Biomechanics detailed the recommendations for 3D kinematics of intervertebral movements (Wu, et al. 2002. J Biomech. 35:543–548), but does not specify how to adapt this proposal to describe the kinematics of the cervical spine, between the head and the thorax. The analysis of the literature shows that no consensus exists at the present time on this subject. The objective of our study was to identify the reference points that formed the most rigid triplet allowing building an optimal thorax segment coordinate system (SCS). We thus measured the variations of distances between markers placed on various anatomical landmarks, and then the deformations of the combinations of three markers on different cervical movements of a sample of 10 asymptomatic subjects. The results show that the triplet formed by the sternum and both acromions undergoes less deformation on the flexion–extension movement. For all the other movements (lateral bending, axial rotation and complex movements), the triplet formed by sternum, T3 and TH (positioned on the thoracic spinal column, in a horizontal plane containing the sternal marker), undergoes less deformation. As a conclusion, the optimal triplet to define the thorax SCS for 3D kinematical analysis of the cervical spine is that formed by the markers: sternum, T3 and TH. This triplet makes it possible to define an orthonormal SCS, the axes of which coincide with anatomical directions, i.e. with the functional axes of the movement.     

Numerical simulation of birch pollen dispersion with an operational weather forecast system

We included a parameterisation of the emissions of pollen grains into the comprehensive model system COSMO-ART. In addition, a detailed density distribution of birch trees within Switzerland was derived. Based on these new developments, we carried out numerical simulations of the dispersion of pollen grains for an episode that occurred in April 2006 over Switzerland and the adjacent regions. Since COSMO-ART is based on the operational forecast model of the German Weather Service, we are presenting a feasibility study of daily pollen forecast based on methods which have been developed during the last two decades for the treatment of anthropogenic aerosol. A comparison of the model results and very detailed pollen counts documents the current possibilities and the shortcomings of the method and gives hints for necessary improvements.     

Numerical simulation of a stroke: computational problems and methodology

We discuss the difficulties of the numerical simulation of a stroke, and we describe the numerical methods which we have developed and used to obtain some realistic results. Nowadays, the computations are performed in two-dimensional slices of a brain, but the strategies to obtain full three-dimensional simulations are explored. This paper is written so as to be understandable by non-mathematicians.     

Analysis of the flow field of the krill,Euphausia pacifica

Velocity measurements were performed for the flow field generated by tethered krill Euphausia pacifica. The particle image velocimetry (PIV) technique was used to measure the velocity field in vertical planes aligned with the krill body axis. The krill generates a narrow jet-like flow behind and below the pleopods (roughly 25° below horizontal). The volume of fluid moving at greater than 10% of the maximum velocity near the pleopods is roughly 18 times larger than the volume of the krill. Thus, the hydrodynamic disturbance occupies a significantly larger region than the animal body. Other krill, sensing the flow disturbance, may take advantage of the flow induced by a neighbor to locate a mate or to draft for efficient propulsion.     

Analysis of the flow field of the krill, Euphausia pacifica

Velocity measurements were performed for the flow field generated by tethered krill Euphausia pacifica. The particle image velocimetry (PIV) technique was used to measure the velocity field in vertical planes aligned with the krill body axis. The krill generates a narrow jet-like flow behind and below the pleopods (roughly 25° below horizontal). The volume of fluid moving at greater than 10% of the maximum velocity near the pleopods is roughly 18 times larger than the volume of the krill. Thus, the hydrodynamic disturbance occupies a significantly larger region than the animal body. Other krill, sensing the flow disturbance, may take advantage of the flow induced by a neighbor to locate a mate or to draft for efficient propulsion.     

Impact of coincidence on granulocyte-platelet complex determination by flow cytometry is evaluated by a novel computer simulation model of coincidence

Cell complexes composed of two different cells labeled with different fluorophores can be detected as double positive events in the flow cytometer. Double positivity can originate not only from real complexes but from non-interacting coinciding cells as well. Coincidence has a high impact on the determination of the amount of platelet–granulocyte complexes since platelet concentration is in the orders of magnitude higher than that of the granulocytes. A computer model has been developed to simulate coincidence in the flow cytometer to reveal the contribution of coincidence to the overestimation of the total amount of platelet–granulocyte complexes. Mixtures of non-interacting fluorescent beads as well as EDTA-anticoagulated blood samples were analyzed in the flow cytometer.

An excellent fit was found between computer simulated and measured data pairs. Bead mixture in the flow cytometer and simulation of that resulted in 37.3 ± 1.3 and 35.7 ± 0.6% double positivity, respectively. 30.2 ± 4.3% double positivity was measured for EDTA-anticoagulated blood samples while simulation of that resulted in 28.3 ± 0.6%. Double positivity attributed to platelet–granulocyte complexes in slightly diluted blood samples might originate in coincidence and not from true complexes.     

Numerical simulation of blood pulsatile flow in a stenosed carotid artery using different rheological models

Symmetrical 30-60% stenosis in a common carotid artery under unsteady flow condition for Newtonian and six non-Newtonian viscosity models are investigated numerically. Results show power-law model produces higher deviations, in terms of velocity and wall shear stress in comparison with other models while generalized power-law and modified-Casson models are more prone to Newtonian state. Comparing separation length of recirculation region at different critical points of cardiac cycle confirms the necessity of considering blood flow in unsteady mode. Increasing stenosis intensity causes flow patterns more disturbed downstream of the stenosis and WSS appear to develop remarkably at the stenosis throat.     

Numerical simulation of trace element transport on subsurface environment pollution in coal mine spoil

An understanding of the dynamic behavior of trace elements leaching from coal mine spoil is important in predicting the groundwater quality. The relationship between trace element concentrations and leaching times, pH values of the media is studied. Column leaching tests conducted in the laboratory showed that there was a close correlation between pH value and trace element concentrations. The longer the leaching time, the higher the trace element concentrations. Different trace elements are differently affected by pH values of leaching media. A numerical model for water flow and trace element transport has been developed based on analyzing the characteristics of migration and transformation of trace elements leached from coal mine spoil. Solutions to the coupled model are accomplished by Eulerian-Lagrangian localized adjoint method. Numerical simulation shows that rainfall intensity determined maximum leaching depth. As rainfall intensity is 3.6ml/s, the outflow concentrations indicate a breakthrough of trace elements beyond the column base, with peak concentration at 90cm depth. And the subsurface pollution range has a trend of increase with time. The model simulations are compared to experimental results of trace element concentrations, with reasonable agreement between them. The analysis and modeling of trace elements suggested that the infiltration of rainwater through the mine spoil might lead to potential groundwater pollution. It provides theoretical evidence for quantitative assessment soil-water quality of trace element transport on environment pollution.     

Numerical simulations of surface plasmon resonance system for monitoring DNA hybridization and detecting protein-lipid film interactions

This paper presents a simple method to extract information about thin organic films from surface plasmon resonance (SPR) spectra. From numerical simulations it was found that a shift (Δθ _SPR) of an absorption peak in the SPR spectrum was directly proportional to the product of the thin organic film thickness and the refractive index difference between the thin organic film and a buffer soaking the sample. It was also found that Δθ _SPR was not sensitive to the thin organic film support of a gold film and a glass cover slip. Relationships between Δθ _SPR and distributions of macromolecule structures, in the thin organic films were theoretically established. Formulae were derived for a homemade SPR system to calculate length, transverse area, density and surface concentration of macromolecules in the thin organic film. The validity of these treatments was checked by precisely measuring the size of a single distearoylphosphatidylcholine molecule on a gold-supported phospholipid film; by quantitatively monitoring hybridization of synthesized oligonucleotides strands based on a biotin/avidin system; and by quantitatively detecting the steric hindrance of rabbit C-reactive protein specifically bound to phospholipid monolayers composed of synthesized lipids. Received: 4 May 1998 / Revised version: 27 July 1998 / Accepted: 27 August 1998     

Numerical simulation of blood flow in the left ventricle and aortic sinus using magnetic resonance imaging and computational fluid dynamics

Understanding cardiac blood flow patterns has many applications in analysing haemodynamics and for the clinical assessment of heart function. In this study, numerical simulations of blood flow in a patient-specific anatomical model of the left ventricle (LV) and the aortic sinus are presented. The realistic 3D geometry of both LV and aortic sinus is extracted from the processing of magnetic resonance imaging (MRI). Furthermore, motion of inner walls of LV and aortic sinus is obtained from cine-MR image analysis and is used as a constraint to a numerical computational fluid dynamics (CFD) model based on the moving boundary approach. Arbitrary Lagrangian–Eulerian finite element method formulation is used for the numerical solution of the transient dynamic equations of the fluid domain. Simulation results include detailed flow characteristics such as velocity, pressure and wall shear stress for the whole domain. The aortic outflow is compared with data obtained by phase-contrast MRI. Good agreement was found between simulation results and these measurements.     

Evaluation of the total antioxidant capacity by using a multipumping flow system with chemiluminescent detection

An automated flow-based procedure for assessment of total antioxidant capacity was developed. It involved a multipumping flow system, a recent approach to flow analysis, and exploited the ability of selected compounds to inhibit the chemiluminescence reactions of luminol or lucigenin with hydrogen peroxide. The system included several discretely actuated solenoid micropumps as the only active components of the flow manifold. This enabled the reproducible insertion and efficient mixing of very low volumes of sample and reagents as well as the transportation of the sample zone toward a flow-through luminometer, where the chemiluminometric response was monitored. With luminol as the chemiluminogenic reagent, linearity of the analytical curves was noted up to 3.2x10(-4), 1.1x10(-3), and 8.8x10(-8) molL-1 for Trolox, ascorbic acid, and resveratrol, respectively. With lucigenin, linear calibration plots up to 2x10(-5) molL-1 of Trolox and 5.7x10(-5)molL-1 of ascorbic acid were obtained. As favorable analytical figures of merit, the measurement precision (RSD typically between 0.2 and 2.0%, n=10), low operational costs, low reagent consumption, sampling rate (160 and 70 h-1), and versatility should be highlighted. The proposed system can be used in distinct analytical circumstances without requiring physical reconfiguration.     

手术前后三维动脉瘤模型的血液动力学模拟对比

本研究目的在于从临床患者的医学影像出发,在手术前后分别对三维动脉瘤模型进行了定常模拟,观察动脉瘤模型内的流场形态以及各血液动力学参数的变化并进行讨论。由于三维的数值模拟比二维的跟接近实际,而且更加直观形象,所以研究结果对于分析动脉瘤的破裂机理具有重要的临床应用价值。     

Numerical analysis of pulsatile blood flow and vessel wall mechanics in different degrees of stenoses

Hemodynamics factors and biomechanical forces play key roles in atherogenesis, plaque development and final rupture. In this paper, we investigated the flow field and stress field for different degrees of stenoses under physiological conditions. Disease is modelled as axisymmetric cosine shape stenoses with varying diameter reductions of 30%, 50% and 70%, respectively. A simulation model which incorporates fluid-structure interaction, a turbulence model and realistic boundary conditions has been developed. The results show that wall motion is constrained at the throat by 60% for the 30% stenosis and 85% for the 50% stenosis; while for the 70% stenosis, wall motion at the throat is negligible through the whole cycle. Peak velocity at the throat varies from 1.47 m/s in the 30% stenosis to 3.2m/s in the 70% stenosis against a value of 0.78 m/s in healthy arteries. Peak wall shear stress values greater than 100 Pa were found for > or =50% stenoses, which in vivo could lead to endothelial stripping. Maximum circumferential stress was found at the shoulders of plaques. The results from this investigation suggest that severe stenoses inhibit wall motion, resulting in higher blood velocities and higher peak wall shear stress, and localization of hoop stress. These factors may contribute to further development and rupture of plaques.     

二维细胞运动模型的数值模拟

研究了一个二维细胞运动模型,该模型是个动边界问题．我们运用运动网格的有限元算法对此模型进行了数值模拟,数值结果表明,细胞最终收敛到固定形状,并且以固定速度向前运动．     

Direct numerical simulation of transitional flow in a patient-specific intracranial aneurysm

In experiments turbulence has previously been shown to occur in intracranial aneurysms. The effects of turbulence induced oscillatory wall stresses could be of great importance in understanding aneurysm rupture. To investigate the effects of turbulence on blood flow in an intracranial aneurysm, we performed a high resolution computational fluid dynamics (CFD) simulation in a patient specific middle cerebral artery (MCA) aneurysm using a realistic, pulsatile inflow velocity. The flow showed transition to turbulence just after peak systole, before relaminarization occurred during diastole. The turbulent structures greatly affected both the frequency of change of wall shear stress (WSS) direction and WSS magnitude, which reached a maximum value of 41.5Pa. The recorded frequencies were predominantly in the range of 1-500Hz. The current study confirms, through properly resolved CFD simulations that turbulence can occur in intracranial aneurysms.