Model studies of blood flow in basilar artery with 3D laser Doppler anemometer
NASA Astrophysics Data System (ADS)
Frolov, S. V.; Sindeev, S. V.; Liepsch, D.; Balasso, A.; Proskurin, S. G.; Potlov, A. Y.
2015-03-01
It is proposed an integrated approach to the study of basilar artery blood flow using 3D laser Doppler anemometer for identifying the causes of the formation and development of cerebral aneurysms. Feature of the work is the combined usage of both mathematical modeling and experimental methods. Described the experimental setup and the method of measurement of basilar artery blood flow, carried out in an interdisciplinary laboratory of Hospital Rechts der Isar of Technical University of Munich. The experimental setup used to simulate the blood flow in the basilar artery and to measure blood flow characteristics using 3D laser Doppler anemometer (3D LDA). Described a method of numerical studies carried out in Tambov State Technical University and the Bakoulev Center for Cardiovascular Surgery. Proposed an approach for sharing experimental and numerical methods of research to identify the causes of the basilar artery aneurysms.
Combination 3D TOP with 2D PC MRA Techique for cerebral blood flow volume measurement.
Guo, G; Wu, R H; Zhang, Y P; Guan, J T; Guo, Y L; Cheng, Y; terBrugge, K; Mikulis, D J
2006-01-01
To demonstrate the discrepancy of cerebral blood flow volume (BFV) estimation with 2D phase-contrast (2D PC) MRA guided with 3D time-of-flight (3D TOF) MR localization by using an "internal" standard. 20 groups of the common (CCA), internal (ICA), and external (ECA) carotid arteries in 10 healthy subjects were examined with 2D PC MRA guided by 3D TOF MR angiograms. The sum BFV of the internal and external carotid arteries was then compared with the ipsilateral common carotid artery flow. An accurate technique would demonstrate no difference. The difference was therefore a measure of accuracy of the method. 3D TOF MRA localization is presented to allow the determination of a slice orientation to improve the accuracy of 2D PC MRA in estimate the BFV. By using the combined protocols, there was better correlation in BFV estimate between the sum of ICA+ECA with the ipsilateral CCA (R2=0.729, P=0.000). The inconsistency (mean +/- SD) was found to be 6.95 +/- 5.95% for estimate the BFV in ICA+ECA and ipsilateral CCA. The main inconsistency was contributed to the ECA and its branches. Guided with 3D TOF MRA localization, 2D PC MRA is more accurate in the determination of blood flow volume in the carotid arteries. PMID:17946401
Grinberg, L; Cheever, E; Anor, T; Madsen, J R; Karniadakis, G E
2011-01-01
We compare results from numerical simulations of pulsatile blood flow in two patient-specific intracranial arterial networks using one-dimensional (1D) and three-dimensional (3D) models. Specifically, we focus on the pressure and flowrate distribution at different segments of the network computed by the two models. Results obtained with 1D and 3D models with rigid walls show good agreement in massflow distribution at tens of arterial junctions and also in pressure drop along the arteries. The 3D simulations with the rigid walls predict higher amplitude of the flowrate and pressure temporal oscillations than the 1D simulations with compliant walls at various segments even for small time-variations in the arterial cross-sectional areas. Sensitivity of the flow and pressure with respect to variation in the elasticity parameters is investigated with the 1D model. PMID:20661645
Parallel Adaptive Computation of Blood Flow in a 3D ``Whole'' Body Model
NASA Astrophysics Data System (ADS)
Zhou, M.; Figueroa, C. A.; Taylor, C. A.; Sahni, O.; Jansen, K. E.
2008-11-01
Accurate numerical simulations of vascular trauma require the consideration of a larger portion of the vasculature than previously considered, due to the systemic nature of the human body's response. A patient-specific 3D model composed of 78 connected arterial branches extending from the neck to the lower legs is constructed to effectively represent the entire body. Recently developed outflow boundary conditions that appropriately represent the downstream vasculature bed which is not included in the 3D computational domain are applied at 78 outlets. In this work, the pulsatile blood flow simulations are started on a fairly uniform, unstructured mesh that is subsequently adapted using a solution-based approach to efficiently resolve the flow features. The adapted mesh contains non-uniform, anisotropic elements resulting in resolution that conforms with the physical length scales present in the problem. The effects of the mesh resolution on the flow field are studied, specifically on relevant quantities of pressure, velocity and wall shear stress.
NASA Astrophysics Data System (ADS)
Carlsohn, Matthias F.; Kemmling, André; Petersen, Arne; Wietzke, Lennart
2016-04-01
Cerebral aneurysms require endovascular treatment to eliminate potentially lethal hemorrhagic rupture by hemostasis of blood flow within the aneurysm. Devices (e.g. coils and flow diverters) promote homeostasis, however, measurement of blood flow within an aneurysm or cerebral vessel before and after device placement on a microscopic level has not been possible so far. This would allow better individualized treatment planning and improve manufacture design of devices. For experimental analysis, direct measurement of real-time microscopic cerebrovascular flow in micro-structures may be an alternative to computed flow simulations. An application of microscopic aneurysm flow measurement on a regular basis to empirically assess a high number of different anatomic shapes and the corresponding effect of different devices would require a fast and reliable method at low cost with high throughout assessment. Transparent three dimensional 3D models of brain vessels and aneurysms may be used for microscopic flow measurements by particle image velocimetry (PIV), however, up to now the size of structures has set the limits for conventional 3D-imaging camera set-ups. On line flow assessment requires additional computational power to cope with the processing large amounts of data generated by sequences of multi-view stereo images, e.g. generated by a light field camera capturing the 3D information by plenoptic imaging of complex flow processes. Recently, a fast and low cost workflow for producing patient specific three dimensional models of cerebral arteries has been established by stereo-lithographic (SLA) 3D printing. These 3D arterial models are transparent an exhibit a replication precision within a submillimeter range required for accurate flow measurements under physiological conditions. We therefore test the feasibility of microscopic flow measurements by PIV analysis using a plenoptic camera system capturing light field image sequences. Averaging across a sequence of
3-D phantom to simulate cerebral blood flow and metabolic images for PET
Hoffman, E.J.; Cutler, P.D.; Digby, W.M.; Mazziotta, J.C. . Nuclear Medicine Lab.)
1990-04-01
A 3-dimensional brain phantom has been developed to simulate the activity distributions found in the human brain in the cerebral blood flow and metabolism studies currently employed in PET. The phantom has a single contiguous chamber and utilizes thin layers of lucite to provide apparent relative concentrations of 5, 1 and 0 for gray matter, white matter and ventricles, respectively, in the brain. The phantom and an ideal image set were created from the same set of data. Thus, the user has a basis for comparing measured images with an ideal image set which enables the user to make quantitative evaluation of the errors in PET studies with a data set similar to that obtained in patient studies.
3D Flow Visualization Using Texture Advection
NASA Technical Reports Server (NTRS)
Kao, David; Zhang, Bing; Kim, Kwansik; Pang, Alex; Moran, Pat (Technical Monitor)
2001-01-01
Texture advection is an effective tool for animating and investigating 2D flows. In this paper, we discuss how this technique can be extended to 3D flows. In particular, we examine the use of 3D and 4D textures on 3D synthetic and computational fluid dynamics flow fields.
Rigano, S; Ferrazzi, E; Boito, S; Pennati, G; Padoan, A; Galan, H
2010-01-01
The primary aim of this pilot study was to study uterine artery (UtA) blood flow volume in uneventful human pregnancies delivered at term, at mid and late gestation by means of 3D and bi-dimensional ultrasound imaging with angio-Doppler combined with fluid-dynamic modeling. Secondary aims were to correlate flow volume to placental site and to UtA Pulsatility Index (PI). Women with singleton, low-risk pregnancies were examined at mid and late gestation. The structure and course of the uterine artery (UtA) was studied in each patient by means of 3D-angio-Doppler and included vessel diameter D, blood flow velocity and PI (measured along the UtA). Fetal weight estimation and placental insertion site were assessed by ultrasound. A robust fluid-dynamic modeling was applied to calculate absolute flow and flow per unit fetal weight. Mean UtA diameter and blood flow velocity increased significantly (p < 0.0001) from mid-gestation to late gestation from 2.6 mm and 67.5 cm/s, to 3.0 mm and 85.3 cm/s, respectively, yielding an increasing absolute flow troughout gestation. h coefficient, derived by fluid-dynamic modeling to calculate mean velocity, increased significantly from 0.52 at mid-gestation to 0.57 at late gestation. UtA blood flow volume ml/min/kg-fetal weight was significantly higher at mid-gestation than at late gestation (535 ml/min/kg vs 193 ml/min/kg; p < 0.0001). In cases with strictly lateral placentas the ipsilateral UtA accommodates at mid and late gestation 63% and 67% of the total UtA flow. In central placentas UtA flow was evenly distributed between the two vessels. An inverse correlation was observed between PI and blood flow volume ml/min/kg (Pearson's coefficient r = -0.54). Our work confirms the technological and methodological limitations in the measurement of uterine artery blood flow. However, Doppler measurements supported by three-dimensional angio imaging of the uterine vessel, high resolution imaging and diameter measurement, and a robust
3D flow focusing for microfluidic flow cytometry with ultrasonics
NASA Astrophysics Data System (ADS)
Gnyawali, Vaskar; Strohm, Eric M.; Daghighi, Yasaman; van de Vondervoort, Mia; Kolios, Michael C.; Tsai, Scott S. H.
2015-11-01
We are developing a flow cytometer that detects unique acoustic signature waves generated from single cells due to interactions between the cells and ultrasound waves. The generated acoustic waves depend on the size and biomechanical properties of the cells and are sufficient for identifying cells in the medium. A microfluidic system capable of focusing cells through a 10 x 10 μm ultrasound beam cross section was developed to facilitate acoustic measurements of single cells. The cells are streamlined in a hydro-dynamically 3D focused flow in a 300 x 300 μm channel made using PDMS. 3D focusing is realized by lateral sheath flows and an inlet needle (inner diameter 100 μm). The accuracy of the 3D flow focusing is measured using a dye and detecting its localization using confocal microscopy. Each flowing cell would be probed by an ultrasound pulse, which has a center frequency of 375 MHz and bandwidth of 250 MHz. The same probe would also be used for recording the scattered waves from the cells, which would be processed to distinguish the physical and biomechanical characteristics of the cells, eventually identifying them. This technique has potential applications in detecting circulating tumor cells, blood cells and blood-related diseases.
Sun, Yawen; Cao, Wenwei; Ding, Weina; Wang, Yao; Han, Xu; Zhou, Yan; Xu, Qun; Zhang, Yong; Xu, Jianrong
2016-01-01
Abnormal reductions in cortical cerebral blood flow (CBF) have been identified in subcortical vascular cognitive impairment (SVCI). However, little is known about the pattern of CBF reduction in relation with the degree of cognitive impairment. CBF measured with three-dimensional (3D) Arterial Spin Labeling (ASL) perfusion magnetic resonance imaging (MRI) helps detect functional changes in subjects with SVCI. We aimed to compare CBF maps in subcortical ischemic vascular disease (SIVD) subjects with and without cognitive impairment and to detect the relationship of the regions of CBF reduction in the brain with the degree of cognitive impairment according to the z-score. A total of 53 subjects with SVCI and 23 matched SIVD subjects without cognitive impairment (controls), underwent a whole-brain 3D ASL MRI in the resting state. Regional CBF (rCBF) was compared voxel wise by using an analysis of variance design in a statistical parametric mapping program, with patient age and sex as covariates. Correlations were calculated between the rCBF value in the whole brain and the z-score in the 53 subjects with SVCI. Compared with the control subjects, SVCI group demonstrated diffuse decreased CBF in the brain. Significant positive correlations were determined in the rCBF values in the left hippocampus, left superior temporal pole gyrus, right superior frontal orbital lobe, right medial frontal orbital lobe, right middle temporal lobe, left thalamus and right insula with the z-scores in SVCI group. The noninvasively quantified resting CBF demonstrated altered CBF distributions in the SVCI brain. The deficit brain perfusions in the temporal and frontal lobe, hippocampus, thalamus and insula was related to the degree of cognitive impairment. Its relationship to cognition indicates the clinical relevance of this functional marker. Thus, our results provide further evidence for the mechanisms underlying the cognitive deficit in patients with SVCI.
Complex flow dynamics around 3D microbot prototypes.
Martínez-Aranda, Sergio; Galindo-Rosales, Francisco J; Campo-Deaño, Laura
2016-02-28
A new experimental setup for the study of the complex flow dynamics around 3D microbot prototypes in a straight microchannel has been developed and assessed. The ultimate aim of this work is focused on the analysis of the morphology of different microbot prototypes to get a better insight into their efficiency when they swim through the main conduits of the human circulatory system. The setup consists of a fused silica straight microchannel with a 3D microbot prototype fastened in the center of the channel cross-section by an extremely thin support. Four different prototypes were considered: a cube, a sphere and two ellipsoids with aspect ratios of 1 : 2 and 1 : 4, respectively. Flow visualization and micro-particle image velocimetry (μPIV) measurements were performed using Newtonian and viscoelastic blood analogue fluids. An efficiency parameter, ℑ, to discriminate the prototypes in terms of flow disturbance has been proposed. PMID:26790959
XML3D and Xflow: combining declarative 3D for the Web with generic data flows.
Klein, Felix; Sons, Kristian; Rubinstein, Dmitri; Slusallek, Philipp
2013-01-01
Researchers have combined XML3D, which provides declarative, interactive 3D scene descriptions based on HTML5, with Xflow, a language for declarative, high-performance data processing. The result lets Web developers combine a 3D scene graph with data flows for dynamic meshes, animations, image processing, and postprocessing. PMID:24808080
Magnetosheath Flow Anomalies in 3-D
NASA Technical Reports Server (NTRS)
Vaisberg, O. L.; Burch, J. L.; Smirnov, V. N.; Avanov, L. A.; Moore, T. E.; Waite, J. H., Jr.; Skalsky, A. A.; Borodkova, N. L.; Coffey, V. N.; Gallagher, D. L.; Rose, M. Franklin (Technical Monitor)
2000-01-01
Measurements of the plasma and magnetic field with high temporal resolution on the Interball Tail probe reveal many flow anomalies in the magnetosheath. They are usually seen as flow direction and number density variations, accompanied by magnetic field discontinuities. Large flow anomalies with number density variations of factor of 2 or more and velocity variations of 100 km/s or more are seen with periodicity of about I per hour. The cases of flow anomalies following in succession are also observed, and suggest their decay while propagating through the magnetosheath. Some magnetospheric disturbances observed in the outer magnetosphere after the satellite has crossed the magnetopause on the inbound orbit suggest their association with magnetosheath flow anomalies observed in the magnetosheath prior to magnetopause crossing.
Recent Advances in Visualizing 3D Flow with LIC
NASA Technical Reports Server (NTRS)
Interrante, Victoria; Grosch, Chester
1998-01-01
Line Integral Convolution (LIC), introduced by Cabral and Leedom in 1993, is an elegant and versatile technique for representing directional information via patterns of correlation in a texture. Although most commonly used to depict 2D flow, or flow over a surface in 3D, LIC methods can equivalently be used to portray 3D flow through a volume. However, the popularity of LIC as a device for illustrating 3D flow has historically been limited both by the computational expense of generating and rendering such a 3D texture and by the difficulties inherent in clearly and effectively conveying the directional information embodied in the volumetric output textures that are produced. In an earlier paper, we briefly discussed some of the factors that may underlie the perceptual difficulties that we can encounter with dense 3D displays and outlined several strategies for more effectively visualizing 3D flow with volume LIC. In this article, we review in more detail techniques for selectively emphasizing critical regions of interest in a flow and for facilitating the accurate perception of the 3D depth and orientation of overlapping streamlines, and we demonstrate new methods for efficiently incorporating an indication of orientation into a flow representation and for conveying additional information about related scalar quantities such as temperature or vorticity over a flow via subtle, continuous line width and color variations.
3D Flow reconstruction using ultrasound PIV
NASA Astrophysics Data System (ADS)
Poelma, C.; Mari, J. M.; Foin, N.; Tang, M.-X.; Krams, R.; Caro, C. G.; Weinberg, P. D.; Westerweel, J.
2011-04-01
Ultrasound particle image velocimetry (PIV) can be used to obtain velocity fields in non-transparent geometries and/or fluids. In the current study, we use this technique to document the flow in a curved tube, using ultrasound contrast bubbles as flow tracer particles. The performance of the technique is first tested in a straight tube, with both steady laminar and pulsatile flows. Both experiments confirm that the technique is capable of reliable measurements. A number of adaptations are introduced that improve the accuracy and applicability of ultrasound PIV. Firstly, due to the method of ultrasound image acquisition, a correction is required for the estimation of velocities from tracer displacements. This correction accounts for the fact that columns in the image are recorded at slightly different instances. The second improvement uses a slice-by-slice scanning approach to obtain three-dimensional velocity data. This approach is here demonstrated in a strongly curved tube. The resulting flow profiles and wall shear stress distribution shows a distinct asymmetry. To meaningfully interpret these three-dimensional results, knowledge of the measurement thickness is required. Our third contribution is a method to determine this quantity, using the correlation peak heights. The latter method can also provide the third (out-of-plane) component if the measurement thickness is known, so that all three velocity components are available using a single probe.
Slope instability in complex 3D topography promoted by convergent 3D groundwater flow
NASA Astrophysics Data System (ADS)
Reid, M. E.; Brien, D. L.
2012-12-01
Slope instability in complex topography is generally controlled by the interaction between gravitationally induced stresses, 3D strengths, and 3D pore-fluid pressure fields produced by flowing groundwater. As an example of this complexity, coastal bluffs sculpted by landsliding commonly exhibit a progression of undulating headlands and re-entrants. In this landscape, stresses differ between headlands and re-entrants and 3D groundwater flow varies from vertical rainfall infiltration to lateral groundwater flow on lower permeability layers with subsequent discharge at the curved bluff faces. In plan view, groundwater flow converges in the re-entrant regions. To investigate relative slope instability induced by undulating topography, we couple the USGS 3D limit-equilibrium slope-stability model, SCOOPS, with the USGS 3D groundwater flow model, MODFLOW. By rapidly analyzing the stability of millions of potential failures, the SCOOPS model can determine relative slope stability throughout the 3D domain underlying a digital elevation model (DEM), and it can utilize both fully 3D distributions of pore-water pressure and material strength. The two models are linked by first computing a groundwater-flow field in MODFLOW, and then computing stability in SCOOPS using the pore-pressure field derived from groundwater flow. Using these two models, our analyses of 60m high coastal bluffs in Seattle, Washington showed augmented instability in topographic re-entrants given recharge from a rainy season. Here, increased recharge led to elevated perched water tables with enhanced effects in the re-entrants owing to convergence of groundwater flow. Stability in these areas was reduced about 80% compared to equivalent dry conditions. To further isolate these effects, we examined groundwater flow and stability in hypothetical landscapes composed of uniform and equally spaced, oscillating headlands and re-entrants with differing amplitudes. The landscapes had a constant slope for both
Lattice Boltzmann Method for 3-D Flows with Curved Boundary
NASA Technical Reports Server (NTRS)
Mei, Renwei; Shyy, Wei; Yu, Dazhi; Luo, Li-Shi
2002-01-01
In this work, we investigate two issues that are important to computational efficiency and reliability in fluid dynamics applications of the lattice, Boltzmann equation (LBE): (1) Computational stability and accuracy of different lattice Boltzmann models and (2) the treatment of the boundary conditions on curved solid boundaries and their 3-D implementations. Three athermal 3-D LBE models (D3QI5, D3Ql9, and D3Q27) are studied and compared in terms of efficiency, accuracy, and robustness. The boundary treatment recently developed by Filippova and Hanel and Met et al. in 2-D is extended to and implemented for 3-D. The convergence, stability, and computational efficiency of the 3-D LBE models with the boundary treatment for curved boundaries were tested in simulations of four 3-D flows: (1) Fully developed flows in a square duct, (2) flow in a 3-D lid-driven cavity, (3) fully developed flows in a circular pipe, and (4) a uniform flow over a sphere. We found that while the fifteen-velocity 3-D (D3Ql5) model is more prone to numerical instability and the D3Q27 is more computationally intensive, the 63Q19 model provides a balance between computational reliability and efficiency. Through numerical simulations, we demonstrated that the boundary treatment for 3-D arbitrary curved geometry has second-order accuracy and possesses satisfactory stability characteristics.
Li, H-O; Huo, R; Xu, G-Q; Duan, Y-H; Nie, P; Ji, X-P; Cheng, Z-P; Xu, Z-D
2015-01-01
Objective: To evaluate the feasibility of high-pitch spiral CT in imaging vascular anomalies (VAs) with affluent blood flow in the head and neck of infants and children. Methods: For patients with suspected VAs and affluent blood flow pre-detected by ultrasound, CT was performed with high-pitch mode, individualized low-dose scan protocol and three-dimensional (3D) reformation. A five-point scale was used for image quality evaluation. Diagnostic accuracy was calculated with clinical diagnosis with/without pathological results as the reference standard. Radiation exposure and single-phase scan time were recorded. Treatment strategies were formulated based on CT images and results and were monitored through follow-up results. Results: 20 lesions were identified in 15 patients (median age of 11 months). The mean score of image quality was 4.13 ± 0.74. 7 patients (7/15, 46.67%) were diagnosed with haemangiomas, 6 patients (6/15, 40%) were diagnosed with venous malformations and 2 patients (2/15, 13.33%) were diagnosed with arteriovenous malformations. The average effective radiation doses of a single phase and of the total procedure were 0.27 ± 0.08 and 0.86 ± 0.21 mSv. The average scanning time of a single phase was 0.46 ± 0.09 s. After treatment, 13 patients (13/15, 86.67%) achieved excellent results, and 2 patients (2/15, 13.33%) showed good results in follow-up visits. Conclusion: High-pitch spiral CT with an individualized low-dose scan protocol and 3D reformation is an effective modality for imaging VAs with affluent blood flow in the head and neck of infants and children when vascular details are needed and ultrasound and MRI could not provide the complete information. Advances in knowledge: This study proposes an alternative modality for imaging VAs with affluent blood flow. PMID:26055504
Customizable engineered blood vessels using 3D printed inserts.
Pinnock, Cameron B; Meier, Elizabeth M; Joshi, Neeraj N; Wu, Bin; Lam, Mai T
2016-04-15
Current techniques for tissue engineering blood vessels are not customizable for vascular size variation and vessel wall thickness. These critical parameters vary widely between the different arteries in the human body, and the ability to engineer vessels of varying sizes could increase capabilities for disease modeling and treatment options. We present an innovative method for producing customizable, tissue engineered, self-organizing vascular constructs by replicating a major structural component of blood vessels - the smooth muscle layer, or tunica media. We utilize a unique system combining 3D printed plate inserts to control construct size and shape, and cell sheets supported by a temporary fibrin hydrogel to encourage cellular self-organization into a tubular form resembling a natural artery. To form the vascular construct, 3D printed inserts are adhered to tissue culture plates, fibrin hydrogel is deposited around the inserts, and human aortic smooth muscle cells are then seeded atop the fibrin hydrogel. The gel, aided by the innate contractile properties of the smooth muscle cells, aggregates towards the center post insert, creating a tissue ring of smooth muscle cells. These rings are then stacked into the final tubular construct. Our methodology is robust, easily repeatable and allows for customization of cellular composition, vessel wall thickness, and length of the vessel construct merely by varying the size of the 3D printed inserts. This platform has potential for facilitating more accurate modeling of vascular pathology, serving as a drug discovery tool, or for vessel repair in disease treatment. PMID:26732049
Flow over bio-inspired 3D herringbone wall riblets
NASA Astrophysics Data System (ADS)
Chen, Huawei; Rao, Fugang; Shang, Xiaopeng; Zhang, Deyuan; Hagiwara, Ichiro
2014-03-01
Under the inspiration of small riblets of shark skin, the microgroove drag reduction riblets whose direction set along fluid flow have been widely investigated. Herringbone-type riblets of bird flight feather are seldom exploited although bird also has excellent flight performance. Inspired from the flight feather, novel bio-inspired plane-3D (p-3D) and spatial-3D (s-3D) herringbone wall riblets are proposed. Through experiment measurement of drag reduction in water tunnel, maximum drag reduction of p-3D and s-3D herringbone riblets was about 17 and 20 %, higher than traditional microgroove riblets. Moreover, significant change of drag reduction was also found by change of the angle between herringbone riblets. In particular, maximum drag reduction occurred as angle between herringbone riblets was about 60° close to real flight feather, which indicates that microstructure of bird flight feather has great impact on flight performance.
3D-reconstruction of blood vessels by ultramicroscopy
Jährling, Nina; Becker, Klaus
2009-01-01
As recently shown, ultramicroscopy (UM) allows 3D-visualization of even large microscopic structures with µm resolution. Thus, it can be applied to anatomical studies of numerous biological and medical specimens. We reconstructed the three-dimensional architecture of tomato-lectin (Lycopersicon esculentum) stained vascular networks by UM in whole mouse organs. The topology of filigree branches of the microvasculature was visualized. Since tumors require an extensive growth of blood vessels to survive, this novel approach may open up new vistas in neurobiology and histology, particularly in cancer research. PMID:20539742
Multigrid calculations of 3-D turbulent viscous flows
NASA Technical Reports Server (NTRS)
Yokota, Jeffrey W.
1989-01-01
Convergence properties of a multigrid algorithm, developed to calculate compressible viscous flows, are analyzed by a vector sequence eigenvalue estimate. The full 3-D Reynolds-averaged Navier-Stokes equations are integrated by an implicit multigrid scheme while a k-epsilon turbulence model is solved, uncoupled from the flow equations. Estimates of the eigenvalue structure for both single and multigrid calculations are compared in an attempt to analyze the process as well as the results of the multigrid technique. The flow through an annular turbine is used to illustrate the scheme's ability to calculate complex 3-D flows.
Numerical analysis of 3-D potential flow in centrifugal turbomachines
NASA Astrophysics Data System (ADS)
Daiguji, H.
1983-09-01
A numerical method is developed for analysing a three-dimensional steady incompressible potential flow through an impeller in centrifugal turbomachines. The method is the same as the previous method which was developed for the axial flow turbomachines, except for some treatments in the downstream region. In order to clarify the validity and limitation of the method, a comparison with the existing experimental data and numerical results is made for radial flow compressor impellers. The calculated blade surface pressure distributions almost coincide with the quasi-3-D calculation by Krimerman and Adler (1978), but are different partly from the quasi-3-D calculation using one meridional flow analysis. It is suggested from this comparison that the flow through an impeller with high efficiency near the design point can be predicted by this fully 3-D numerical method.
USM3D Predictions of Supersonic Nozzle Flow
NASA Technical Reports Server (NTRS)
Carter, Melissa B.; Elmiligui, Alaa A.; Campbell, Richard L.; Nayani, Sudheer N.
2014-01-01
This study focused on the NASA Tetrahedral Unstructured Software System CFD code (USM3D) capability to predict supersonic plume flow. Previous studies, published in 2004 and 2009, investigated USM3D's results versus historical experimental data. This current study continued that comparison however focusing on the use of the volume souring to capture the shear layers and internal shock structure of the plume. This study was conducted using two benchmark axisymmetric supersonic jet experimental data sets. The study showed that with the use of volume sourcing, USM3D was able to capture and model a jet plume's shear layer and internal shock structure.
An elliptic calculation procedure for 3-D viscous flow
NASA Astrophysics Data System (ADS)
Moore, J. G.
1985-05-01
The computation of 3-D internal transonic flows by means of a 3-D Euler Code is discussed. A multidomain approach for time hyperbolic system is presented. This technique, based on the decomposition of the computational domain into several subdomains which may overlap one another, makes it possible to simplify some mesh generation problems and to fit discontinuities such as shocks and slip surfaces. A description of the 3-D Euler Code is given. The space discretization method and the treatment of boundary conditions are emphasized. Various applications of this code in turbomachinery are discussed.
3D Printed Micro Free-Flow Electrophoresis Device.
Anciaux, Sarah K; Geiger, Matthew; Bowser, Michael T
2016-08-01
The cost, time, and restrictions on creative flexibility associated with current fabrication methods present significant challenges in the development and application of microfluidic devices. Additive manufacturing, also referred to as three-dimensional (3D) printing, provides many advantages over existing methods. With 3D printing, devices can be made in a cost-effective manner with the ability to rapidly prototype new designs. We have fabricated a micro free-flow electrophoresis (μFFE) device using a low-cost, consumer-grade 3D printer. Test prints were performed to determine the minimum feature sizes that could be reproducibly produced using 3D printing fabrication. Microfluidic ridges could be fabricated with dimensions as small as 20 μm high × 640 μm wide. Minimum valley dimensions were 30 μm wide × 130 μm wide. An acetone vapor bath was used to smooth acrylonitrile-butadiene-styrene (ABS) surfaces and facilitate bonding of fully enclosed channels. The surfaces of the 3D-printed features were profiled and compared to a similar device fabricated in a glass substrate. Stable stream profiles were obtained in a 3D-printed μFFE device. Separations of fluorescent dyes in the 3D-printed device and its glass counterpart were comparable. A μFFE separation of myoglobin and cytochrome c was also demonstrated on a 3D-printed device. Limits of detection for rhodamine 110 were determined to be 2 and 0.3 nM for the 3D-printed and glass devices, respectively. PMID:27377354
3-D Flow Visualization with a Light-field Camera
NASA Astrophysics Data System (ADS)
Thurow, B.
2012-12-01
Light-field cameras have received attention recently due to their ability to acquire photographs that can be computationally refocused after they have been acquired. In this work, we describe the development of a light-field camera system for 3D visualization of turbulent flows. The camera developed in our lab, also known as a plenoptic camera, uses an array of microlenses mounted next to an image sensor to resolve both the position and angle of light rays incident upon the camera. For flow visualization, the flow field is seeded with small particles that follow the fluid's motion and are imaged using the camera and a pulsed light source. The tomographic MART algorithm is then applied to the light-field data in order to reconstruct a 3D volume of the instantaneous particle field. 3D, 3C velocity vectors are then determined from a pair of 3D particle fields using conventional cross-correlation algorithms. As an illustration of the concept, 3D/3C velocity measurements of a turbulent boundary layer produced on the wall of a conventional wind tunnel are presented. Future experiments are planned to use the camera to study the influence of wall permeability on the 3-D structure of the turbulent boundary layer.Schematic illustrating the concept of a plenoptic camera where each pixel represents both the position and angle of light rays entering the camera. This information can be used to computationally refocus an image after it has been acquired. Instantaneous 3D velocity field of a turbulent boundary layer determined using light-field data captured by a plenoptic camera.
MedlinePlus Videos and Cool Tools
As the heart pumps, the arteries carry oxygen-rich blood (shown in red) away from the heart and toward the body’s tissues and ... returns to the heart from the lungs, which pumps it throughout the body.
Vertical Flow Lithography for Fabrication of 3D Anisotropic Particles.
Habasaki, Shohei; Lee, Won Chul; Yoshida, Shotaro; Takeuchi, Shoji
2015-12-22
A microfluidics-based method for the 3D fabrication of anisotropic particles is reported. The method uses a vertical microchannel where tunable light patterns solidify photocurable resins for stacking multiple layers of the resins, thus enabling an application of stereolithography concepts to conventional flow lithography. Multilayered, tapered, and angular compartmental microparticles are demonstrated. PMID:26551590
Joint 3d Estimation of Vehicles and Scene Flow
NASA Astrophysics Data System (ADS)
Menze, M.; Heipke, C.; Geiger, A.
2015-08-01
driving. While much progress has been made in recent years, imaging conditions in natural outdoor environments are still very challenging for current reconstruction and recognition methods. In this paper, we propose a novel unified approach which reasons jointly about 3D scene flow as well as the pose, shape and motion of vehicles in the scene. Towards this goal, we incorporate a deformable CAD model into a slanted-plane conditional random field for scene flow estimation and enforce shape consistency between the rendered 3D models and the parameters of all superpixels in the image. The association of superpixels to objects is established by an index variable which implicitly enables model selection. We evaluate our approach on the challenging KITTI scene flow dataset in terms of object and scene flow estimation. Our results provide a prove of concept and demonstrate the usefulness of our method.
Two-equation turbulence modeling for 3-D hypersonic flows
NASA Technical Reports Server (NTRS)
Bardina, J. E.; Coakley, T. J.; Marvin, J. G.
1992-01-01
An investigation to verify, incorporate and develop two-equation turbulence models for three-dimensional high speed flows is presented. The current design effort of hypersonic vehicles has led to an intensive study of turbulence models for compressible hypersonic flows. This research complements an extensive review of experimental data and the current development of 2D turbulence models. The review of experimental data on 2D and 3D flows includes complex hypersonic flows with pressure profiles, skin friction, wall heat transfer, and turbulence statistics data. In a parallel effort, turbulence models for high speed flows have been tested against flat plate boundary layers, and are being tested against the 2D database. In the present paper, we present the results of 3D Navier-Stokes numerical simulations with an improved k-omega two-equation turbulence model against experimental data and empirical correlations of an adiabatic flat plate boundary layer, a cold wall flat plate boundary layer, and a 3D database flow, the interaction of an oblique shock wave and a thick turbulent boundary layer with a free stream Mach number = 8.18 and Reynolds number = 5 x 10 to the 6th.
An annotation system for 3D fluid flow visualization
NASA Technical Reports Server (NTRS)
Loughlin, Maria M.; Hughes, John F.
1995-01-01
Annotation is a key activity of data analysis. However, current systems for data analysis focus almost exclusively on visualization. We propose a system which integrates annotations into a visualization system. Annotations are embedded in 3D data space, using the Post-it metaphor. This embedding allows contextual-based information storage and retrieval, and facilitates information sharing in collaborative environments. We provide a traditional database filter and a Magic Lens filter to create specialized views of the data. The system has been customized for fluid flow applications, with features which allow users to store parameters of visualization tools and sketch 3D volumes.
MedlinePlus Videos and Cool Tools
As the heart pumps, the arteries carry oxygen-rich blood (shown in red) away from the heart and toward the body's tissues and vital organs. ... brain, liver, kidneys, stomach, and muscles, including the heart muscle itself. At the same time, the veins ...
The 3D Flow Field Around an Embedded Planet
NASA Astrophysics Data System (ADS)
Fung, Jeffrey; Artymowicz, Pawel; Wu, Yanqin
2015-10-01
3D modifications to the well-studied 2D flow topology around an embedded planet have the potential to resolve long-standing problems in planet formation theory. We present a detailed analysis of the 3D isothermal flow field around a 5 Earth-mass planet on a fixed circular orbit, simulated using our graphics processing unit hydrodynamics code PEnGUIn. We find that, overall, the horseshoe region has a columnar structure extending vertically much beyond the Hill sphere of the planet. This columnar structure is only broken for some of the widest horseshoe streamlines, along which high altitude fluid descends rapidly into the planet’s Bondi sphere, performs one horseshoe turn, and exits the Bondi sphere radially in the midplane. A portion of this flow exits the horseshoe region altogether, which we refer to as the “transient” horseshoe flow. The flow continues as it rolls up into a pair of up-down symmetric horizontal vortex lines shed into the wake of the planet. This flow, unique to 3D, affects both planet accretion and migration. It prevents the planet from sustaining a hydrostatic atmosphere due to its intrusion into the Bondi sphere, and leads to a significant corotation torque on the planet, unanticipated by 2D analysis. In the reported simulation, starting with a {{Σ }}˜ {r}-3/2 radial surface density profile, this torque is positive and partially cancels with the negative differential Lindblad torque, resulting in a factor of three slower planet migration rate. Finally, we report 3D effects can be suppressed by a sufficiently large disk viscosity, leading to results similar to 2D.
Optimization of a Hybrid Magnetic Bearing for a Magnetically Levitated Blood Pump via 3-D FEA.
Cheng, Shanbao; Olles, Mark W; Burger, Aaron F; Day, Steven W
2011-10-01
In order to improve the performance of a magnetically levitated (maglev) axial flow blood pump, three-dimensional (3-D) finite element analysis (FEA) was used to optimize the design of a hybrid magnetic bearing (HMB). Radial, axial, and current stiffness of multiple design variations of the HMB were calculated using a 3-D FEA package and verified by experimental results. As compared with the original design, the optimized HMB had twice the axial stiffness with the resulting increase of negative radial stiffness partially compensated for by increased current stiffness. Accordingly, the performance of the maglev axial flow blood pump with the optimized HMBs was improved: the maximum pump speed was increased from 6000 rpm to 9000 rpm (50%). The radial, axial and current stiffness of the HMB was found to be linear at nominal operational position from both 3-D FEA and empirical measurements. Stiffness values determined by FEA and empirical measurements agreed well with one another. The magnetic flux density distribution and flux loop of the HMB were also visualized via 3-D FEA which confirms the designers' initial assumption about the function of this HMB. PMID:22065892
Optimization of a Hybrid Magnetic Bearing for a Magnetically Levitated Blood Pump via 3-D FEA
Cheng, Shanbao; Olles, Mark W.; Burger, Aaron F.; Day, Steven W.
2011-01-01
In order to improve the performance of a magnetically levitated (maglev) axial flow blood pump, three-dimensional (3-D) finite element analysis (FEA) was used to optimize the design of a hybrid magnetic bearing (HMB). Radial, axial, and current stiffness of multiple design variations of the HMB were calculated using a 3-D FEA package and verified by experimental results. As compared with the original design, the optimized HMB had twice the axial stiffness with the resulting increase of negative radial stiffness partially compensated for by increased current stiffness. Accordingly, the performance of the maglev axial flow blood pump with the optimized HMBs was improved: the maximum pump speed was increased from 6000 rpm to 9000 rpm (50%). The radial, axial and current stiffness of the HMB was found to be linear at nominal operational position from both 3-D FEA and empirical measurements. Stiffness values determined by FEA and empirical measurements agreed well with one another. The magnetic flux density distribution and flux loop of the HMB were also visualized via 3-D FEA which confirms the designers’ initial assumption about the function of this HMB. PMID:22065892
Measuring Actin Flow in 3D Cell Protrusions
Chiu, Chi-Li; Digman, Michelle A.; Gratton, Enrico
2013-01-01
Actin dynamics is important in determining cell shape, tension, and migration. Methods such as fluorescent speckle microscopy and spatial temporal image correlation spectroscopy have been used to capture high-resolution actin turnover dynamics within cells in two dimensions. However, these methods are not directly applicable in 3D due to lower resolution and poor contrast. Here, we propose to capture actin flow in 3D with high spatial-temporal resolution by combining nanoscale precise imaging by rapid beam oscillation and fluctuation spectroscopy techniques. To measure the actin flow along cell protrusions in cell expressing actin-eGFP cultured in a type I collagen matrix, the laser was orbited around the protrusion and its trajectory was modulated in a clover-shaped pattern perpendicularly to the protrusion. Orbits were also alternated at two positions closely spaced along the protrusion axis. The pair cross-correlation function was applied to the fluorescence fluctuation from these two positions to capture the flow of actin. Measurements done on nonmoving cellular protrusion tips showed no pair-correlation at two orbital positions indicating a lack of flow of F-actin bundles. However, in some protrusions, the pair-correlation approach revealed directional flow of F-actin bundles near the protrusion surface with flow rates in the range of ∼1 μm/min, comparable to results in two dimensions using fluorescent speckle microscopy. Furthermore, we found that the actin flow rate is related to the distance to the protrusion tip. We also observed collagen deformation by concomitantly detecting collagen fibers with reflectance detection during these actin motions. The implementation of the nanoscale precise imaging by rapid beam oscillation method with a cloverleaf-shaped trajectory in conjunction with the pair cross-correlation function method provides a quantitative way of capturing dynamic flows and organization of proteins during cell migration in 3D in conditions of
Finite element solver for 3-D compressible viscous flows
NASA Technical Reports Server (NTRS)
Reddy, K. C.; Reddy, J. N.
1986-01-01
The space shuttle main engine (SSME) has extremely complex internal flow structure. The geometry of the flow domain is three-dimensional with complicated topology. The flow is compressible, viscous, and turbulent with large gradients in flow quantities and regions of recirculations. The analysis of the flow field in SSME involves several tedious steps. One is the geometrical modeling of the particular zone of the SSME being studied. Accessing the geometry definition, digitalizing it, and developing surface interpolations suitable for an interior grid generator require considerable amount of manual labor. There are several types of grid generators available with some general-purpose finite element programs. An efficient and robust computational scheme for solving 3D Navier-Stokes equations has to be implemented. Post processing software has to be adapted to visualize and analyze the computed 3D flow field. The progress made in a project to develop software for the analysis of the flow is discussed. The technical approach to the development of the finite element scheme and the relaxation procedure are discussed. The three dimensional finite element code for the compressible Navier-Stokes equations is listed.
A finite element solver for 3-D compressible viscous flows
NASA Technical Reports Server (NTRS)
Reddy, K. C.; Reddy, J. N.; Nayani, S.
1990-01-01
Computation of the flow field inside a space shuttle main engine (SSME) requires the application of state of the art computational fluid dynamic (CFD) technology. Several computer codes are under development to solve 3-D flow through the hot gas manifold. Some algorithms were designed to solve the unsteady compressible Navier-Stokes equations, either by implicit or explicit factorization methods, using several hundred or thousands of time steps to reach a steady state solution. A new iterative algorithm is being developed for the solution of the implicit finite element equations without assembling global matrices. It is an efficient iteration scheme based on a modified nonlinear Gauss-Seidel iteration with symmetric sweeps. The algorithm is analyzed for a model equation and is shown to be unconditionally stable. Results from a series of test problems are presented. The finite element code was tested for couette flow, which is flow under a pressure gradient between two parallel plates in relative motion. Another problem that was solved is viscous laminar flow over a flat plate. The general 3-D finite element code was used to compute the flow in an axisymmetric turnaround duct at low Mach numbers.
Unsteady 3D Turbulent Flow Separation around a ROV Body
NASA Astrophysics Data System (ADS)
Ungureanu, Costel; Lungu, Adrian
2009-09-01
Turbulent separated flows around ellipsoids of various aspect ratios are investigated using a numerical method. The Reynolds averaged equations for continuity and momentum are solved by cell-centered finite-volume method for the primitive variables to describe the 3D turbulent incompressible flow. The objectives of the study are: (a) to investigate the propulsive performances of a moving Remotely Operated Vehicle (ROV hereafter) ellipsoidal body; (b) to validate the computational solutions through comparisons with the experimental data; (c) to investigate the effects of the angle of attack on the separation pattern as well as on the hydrodynamic forces and moments.
Challenges in Lagrangian transport and predictability in 3D flows
NASA Astrophysics Data System (ADS)
Branicki, M.; Wiggins, S.; Kirwan, A. D.; Malek-Madani, R.
2011-12-01
The interplay between the geometrical theory of dynamical systems and the trajectory-based description of aperiodically time-dependent fluid flows has led to significant advances in understanding the role of chaotic transport in geophysical flows at scales dominated by advection. Lagrangian transport analysis utilizing either the time-dependent geometry of intersecting stable and unstable manifolds of the so-called Distinguished Hyperbolic Trajectories (DHT), or ridges of finite-time Lyapunov exponent fields (LCS), provide a much needed and complementary insight into ephemeral mechanisms responsible for the existence of `leaky' transport barriers and 'leaky' mesoscale eddies. However, to date most oceanic applications have been confined to 2D analysis of near surface regions in 'perfect' flows not accounting for model or measurement error, and with the tacit assumption of negligible vertical velocities. I will systematically address issues concerning the regimes of applicability of two-dimensional analysis in 3D aperiodically time-dependent flows, as well as outstanding challenges in fully 3D Lagrangian transport analysis. Even for perfect horizontal velocities, little is known about the vertical extent of stable/unstable manifolds associated with DHTs and/or other special structures relevant to stratified 3D flows. In particular, their sensitivity to errors in the vertical velocities and data assimilation methods has been little studied. Rigorous results regarding the above issues will be illustrated by revealing and mathematically tractable toy models, as well as examples from a detailed study in an eddy-rich region from the Gulf of Mexico and the Mediterranean. New ways of quantifying the uncertainty in Lagrangian predictions will also be presented.
3D vector flow using a row-column addressed CMUT array
NASA Astrophysics Data System (ADS)
Holbek, Simon; Christiansen, Thomas Lehrmann; Engholm, Mathias; Lei, Anders; Stuart, Mathias Bo; Beers, Christopher; Moesner, Lars Nordahl; Bagge, Jan Peter; Thomsen, Erik Vilain; Jensen, Jørgen Arendt
2016-04-01
This paper presents an in-house developed 2-D capacitive micromachined ultrasonic transducer (CMUT) applied for 3-D blood flow estimation. The probe breaks with conventional transducers in two ways; first, the ultrasonic pressure field is generated from thousands of small vibrating micromachined cells, and second, elements are accessed by row and/or column indices. The 62+62 2-D row-column addressed prototype CMUT probe was used for vector flow estimation by transmitting focused ultrasound into a flow-rig with a fully developed parabolic flow. The beam-to-flow angle was 90°. The received data was beamformed and processed offline. A transverse oscillation (TO) velocity estimator was used to estimate the 3-D vector flow along a line originating from the center of the transducer. The estimated velocities in the lateral and axial direction were close to zero as expected. In the transverse direction a characteristic parabolic velocity profile was estimated with a peak velocity of 0.48 m/s +/- 0.02 m/s in reference to the expected 0.54 m/s. The results presented are the first 3-D vector flow estimates obtained with a row-column CMUT probe, which demonstrates that the CMUT technology is feasible for 3-D flow estimation.
Using flow information to support 3D vessel reconstruction from rotational angiography
Waechter, Irina; Bredno, Joerg; Weese, Juergen; Barratt, Dean C.; Hawkes, David J.
2008-07-15
For the assessment of cerebrovascular diseases, it is beneficial to obtain three-dimensional (3D) morphologic and hemodynamic information about the vessel system. Rotational angiography is routinely used to image the 3D vascular geometry and we have shown previously that rotational subtraction angiography has the potential to also give quantitative information about blood flow. Flow information can be determined when the angiographic sequence shows inflow and possibly outflow of contrast agent. However, a standard volume reconstruction assumes that the vessel tree is uniformly filled with contrast agent during the whole acquisition. If this is not the case, the reconstruction exhibits artifacts. Here, we show how flow information can be used to support the reconstruction of the 3D vessel centerline and radii in this case. Our method uses the fast marching algorithm to determine the order in which voxels are analyzed. For every voxel, the rotational time intensity curve (R-TIC) is determined from the image intensities at the projection points of the current voxel. Next, the bolus arrival time of the contrast agent at the voxel is estimated from the R-TIC. Then, a measure of the intensity and duration of the enhancement is determined, from which a speed value is calculated that steers the propagation of the fast marching algorithm. The results of the fast marching algorithm are used to determine the 3D centerline by backtracking. The 3D radius is reconstructed from 2D radius estimates on the projection images. The proposed method was tested on computer simulated rotational angiography sequences with systematically varied x-ray acquisition, blood flow, and contrast agent injection parameters and on datasets from an experimental setup using an anthropomorphic cerebrovascular phantom. For the computer simulation, the mean absolute error of the 3D centerline and 3D radius estimation was 0.42 and 0.25 mm, respectively. For the experimental datasets, the mean absolute
Optic flow aided navigation and 3D scene reconstruction
NASA Astrophysics Data System (ADS)
Rollason, Malcolm
2013-10-01
An important enabler for low cost airborne systems is the ability to exploit low cost inertial instruments. An Inertial Navigation System (INS) can provide a navigation solution, when GPS is denied, by integrating measurements from inertial sensors. However, the gyrometer and accelerometer biases of low cost inertial sensors cause compound errors in the integrated navigation solution. This paper describes experiments to establish whether (and to what extent) the navigation solution can be aided by fusing measurements from an on-board video camera with measurements from the inertial sensors. The primary aim of the work was to establish whether optic flow aided navigation is beneficial even when the 3D structure within the observed scene is unknown. A further aim was to investigate whether an INS can help to infer 3D scene content from video. Experiments with both real and synthetic data have been conducted. Real data was collected using an AR Parrot quadrotor. Empirical results illustrate that optic flow provides a useful aid to navigation even when the 3D structure of the observed scene is not known. With optic flow aiding of the INS, the computed trajectory is consistent with the true camera motion, whereas the unaided INS yields a rapidly increasing position error (the data represents ~40 seconds, after which the unaided INS is ~50 metres in error and has passed through the ground). The results of the Monte Carlo simulation concur with the empirical result. Position errors, which grow as a quadratic function of time when unaided, are substantially checked by the availability of optic flow measurements.
Patterns of 3D flow in a rotating cylinder array
NASA Astrophysics Data System (ADS)
Craig, Anna; Dabiri, John; Koseff, Jeffrey
2015-11-01
Experimental data are presented for large arrays of rotating, finite-height cylinders, which show that the three-dimensional flows are strongly dependent on the geometric and rotational configurations of the array. Two geometric configurations of the cylinders, each with two rotational configurations, were examined for a total of four arrays. 2D PIV was conducted in multiple intersecting horizontal and vertical sheets at a location far downstream of the leading edge of the array in order to build up a picture of the 3D developed flow patterns. It was found that the rotation of the cylinders drives the formation of streamwise and transverse flow patterns between cylinders. These horizontal flow patterns, by conservation of mass, drive vertical flows through the top of the array. As the array of rotating cylinders may provide insight into the flow kinematics of an array of vertical axis wind turbines, this planform flux is of particular interest as it would bring down into the array high kinetic energy fluid from above the array, thus increasing the energy resource available to turbines far downstream of the leading edge of the array.
Supersonic turbulence in 3D isothermal flow collision
NASA Astrophysics Data System (ADS)
Folini, Doris; Walder, Rolf; Favre, Jean M.
2014-02-01
Large scale supersonic bulk flows are present in a wide range of astrophysical objects, from O-star winds to molecular clouds, galactic sheets, accretion, or γ-ray bursts. Associated flow collisions shape observable properties and internal physics alike. Our goal is to shed light on the interplay between large scale aspects of such collision zones and the characteristics of the compressible turbulence they harbor. Our model setup is as simple as can be: 3D hydrodynamical simulations of two head-on colliding, isothermal, and homogeneous flows with identical upstream (subscript u) flow parameters and Mach numbers 2 < Mu < 43. The turbulence in the collision zone is driven by the upstream flows, whose kinetic energy is partly dissipated and spatially modulated by the shocks confining the zone. Numerical results are in line with expectations from self-similarity arguments. The spatial scale of modulation grows with the collision zone. The fraction of energy dissipated at the confining shocks decreases with increasing Mu. The mean density is ρm ≈ 20ρu, independent of Mu. The root mean square Mach number is Mrms ≈ 0.25Mu. Deviations toward weaker turbulence are found as the collision zone thickens and for small Mu. The density probability function is not log-normal. The turbulence is inhomogeneous, weaker in the center of the zone than close to the confining shocks. It is also anisotropic: transverse to the upstream flows Mrms is always subsonic. We argue that uniform, head-on colliding flows generally disfavor turbulence that is at the same time isothermal, supersonic, and isotropic. The anisotropy carries over to other quantities like the density variance - Mach number relation. Line-of-sight effects thus exist. Structure functions differ depending on whether they are computed along a line-of-sight perpendicular or parallel to the upstream flow. Turbulence characteristics generally deviate markedly from those found for uniformly driven, supersonic, isothermal
MPSalsa 3D Simulations of Chemically Reacting Flows
Many important scientific and engineering applications require a detailed analysis of complex systems with coupled fluid flow, thermal energy transfer, mass transfer and nonequilibrium chemical reactions. Currently, computer simulations of these complex reacting flow problems are limited to idealized systems in one or two spatial dimensions when coupled with a detailed, fundamental chemistry model. The goal of our research is to develop, analyze and implement advanced MP numerical algorithms that will allow high resolution 3D simulations with an equal emphasis on fluid flow and chemical kinetics modeling. In our research, we focus on the development of new, fully coupled, implicit solution strategies that are based on robust MP iterative solution methods (copied from http://www.cs.sandia.gov/CRF/MPSalsa/). These simulations are needed for scientific and technical areas such as: combustion research for transportation, atmospheric chemistry modeling for pollution studies, chemically reacting flow models for analysis and control of manufacturing processes, surface catalytic reactors for methane to methanol conversion and chemical vapor deposition (CVD) process modeling for production of advanced semiconductor materials (http://www.cs.sandia.gov/CRF/MPSalsa/).
This project website provides six QuickTime videos of these simulations, along with a small image gallery and slideshow animations. A list of related publications and conference presentations is also made available.
Modeling Electric Current Flow in 3D Fractured Media
NASA Astrophysics Data System (ADS)
Demirel, S.; Roubinet, D.; Irving, J.
2014-12-01
The study of fractured rocks is extremely important in a variety of research fields and applications such as hydrogeology, hydrocarbon extraction and long-term storage of toxic waste. As fractures are highly conductive structures in comparison to the surrounding rock, their presence can be either an advantage or a drawback. For hydrocarbon extraction, fractures allow for quick and easy access to the resource whereas for toxic waste storage their presence increases the risk of leakage and migration of pollutants. In both cases, the identification of fracture network characteristics is an essential step. Recently, we have developed an approach for modeling electric current flow in 2D fractured media. This approach is based on a discrete-dual-porosity model where fractures are represented explicitly, the matrix is coarsely discretized into blocks, and current flow exchange between the fractures and matrix is analytically evaluated at the fracture-scale and integrated at the block-scale [1]. Although this approach has shown much promise and has proven its efficiency for 2D simulations, its extension to 3D remains to be addressed. To this end, we assume that fractures can be represented as two-dimensional finite planes embedded in the surrounding matrix, and we express analytically the distribution of electric potential at the fracture scale. This fracture-scale expression takes into account the electric-current-flow exchange with the surrounding matrix and flow conservation is enforced at the fracture intersections. The fracture-matrix exchange is then integrated at the matrix-block scale where the electric current flow conservation at the block boundaries is formulated with a modified finite volume method. With the objective of providing a low-computational-cost modeling approach adapted to 3D simulations in fractured media, our model is (i) validated and compared to existing modeling approaches and, (ii) used to evaluate the impact of the presence of fractures on
Energy flow in passive and active 3D cochlear model
Wang, Yanli; Steele, Charles; Puria, Sunil
2015-12-31
Energy flow in the cochlea is an important characteristic of the cochlear traveling wave, and many investigators, such as von Békésy and Lighthill, have discussed this phenomenon. Particularly after the discovery of the motility of the outer hair cells (OHCs), the nature of the power gain of the cochlea has been a fundamental research question. In the present work, direct three-dimensional (3D) calculations of the power on cross sections of the cochlea and on the basilar membrane are performed based on a box model of the mouse cochlea. The distributions of the fluid pressure and fluid velocity in the scala vestibuli are presented. The power output from the OHCs and the power loss due to fluid viscous damping are calculated along the length of the cochlea. This work provides a basis for theoretical calculations of the power gain of the OHCs from mechanical considerations.
Energy flow in passive and active 3D cochlear model
NASA Astrophysics Data System (ADS)
Wang, Yanli; Puria, Sunil; Steele, Charles
2015-12-01
Energy flow in the cochlea is an important characteristic of the cochlear traveling wave, and many investigators, such as von Békésy and Lighthill, have discussed this phenomenon. Particularly after the discovery of the motility of the outer hair cells (OHCs), the nature of the power gain of the cochlea has been a fundamental research question. In the present work, direct three-dimensional (3D) calculations of the power on cross sections of the cochlea and on the basilar membrane are performed based on a box model of the mouse cochlea. The distributions of the fluid pressure and fluid velocity in the scala vestibuli are presented. The power output from the OHCs and the power loss due to fluid viscous damping are calculated along the length of the cochlea. This work provides a basis for theoretical calculations of the power gain of the OHCs from mechanical considerations.
Slat Cove Unsteadiness Effect of 3D Flow Structures
NASA Technical Reports Server (NTRS)
Choudhari, Meelan M.; Khorrami, Mehdi R.
2006-01-01
Previous studies have indicated that 2D, time accurate computations based on a pseudo-laminar zonal model of the slat cove region (within the framework of the Reynolds-Averaged Navier-Stokes equations) are inadequate for predicting the full unsteady dynamics of the slat cove flow field. Even though such computations could capture the large-scale, unsteady vorticity structures in the slat cove region without requiring any external forcing, the simulated vortices were excessively strong and the recirculation zone was unduly energetic in comparison with the PIV measurements for a generic high-lift configuration. To resolve this discrepancy and to help enable physics based predictions of slat aeroacoustics, the present paper is focused on 3D simulations of the slat cove flow over a computational domain of limited spanwise extent. Maintaining the pseudo-laminar approach, current results indicate that accounting for the three-dimensionality of flow fluctuations leads to considerable improvement in the accuracy of the unsteady, nearfield solution. Analysis of simulation data points to the likely significance of turbulent fluctuations near the reattachment region toward the generation of broadband slat noise. The computed acoustic characteristics (in terms of the frequency spectrum and spatial distribution) within short distances from the slat resemble the previously reported, subscale measurements of slat noise.
Simulations of soluble surfactants in 3D multiphase flow
NASA Astrophysics Data System (ADS)
Muradoglu, Metin; Tryggvason, Gretar
2014-10-01
A finite-difference/front-tracking method is developed for simulations of soluble surfactants in 3D multiphase flows. The interfacial and bulk surfactant concentration evolution equations are solved fully coupled with the incompressible Navier-Stokes equations. A non-linear equation of state is used to relate interfacial surface tension to surfactant concentration at the interface. Simple test cases are designed to validate different parts of the numerical algorithm and the computational results are found to be in a good agreement with the analytical solutions. The numerical algorithm is parallelized using a domain-decomposition method. It is then applied to study the effects of soluble surfactants on the motion of buoyancy-driven bubbles in a straight square channel in nearly undeformable (spherical) and deformable (ellipsoidal) regimes. Finally the method is used to examine the effects of soluble surfactants on the lateral migration of bubbles in a pressure-driven channel flow. It is found that surfactant-induced Marangoni stresses counteract the shear-induced lift force and can reverse the lateral bubble migration completely, i.e., the contaminated bubble drifts away from the channel wall and stabilizes at the center of the channel when the surfactant-induced Marangoni stresses are sufficiently large.
ICEd-ALE Treatment of 3-D Fluid Flow.
Energy Science and Technology Software Center (ESTSC)
1999-09-13
Version: 00 SALE3D calculates three-dimensional fluid flow at all speeds, from the incompressible limit to highly supersonic. An implicit treatment of the pressure calculation similar to that in the Implicit Continuous-fluid Eulerian (ICE) technique provides this flow speed flexibility. In addition, the computing mesh may move with the fluid in a typical Lagrangian fashion, be held in an Eulerian manner, or move in some arbitrarily specified way to provide a continuous rezoning capability. This latitudemore » results from use of an Arbitrary Lagrangian-Eulerian (ALE) treatment of the mesh. The partial differential equations solved are the Navier-Stokes equations and the mass and internal energy equations. The fluid pressure is determined from an equation of state and supplemented with an artificial viscous pressure for the computation of shock waves. The computing mesh consists of a three-dimensional network of arbitrarily shaped, six-sided deformable cells, and a variety of user-selectable boundary conditions are provided in the program.« less
Antibody-free isolation of rare cancer cells from blood based on 3D lateral dielectrophoresis.
Cheng, I-Fang; Huang, Wei-Lun; Chen, Tzu-Ying; Liu, Chien-Wei; Lin, Yu-De; Su, Wu-Chou
2015-07-21
We present an antibody-free approach for the high-purity and high-throughput dielectrophoretic (DEP) isolation of circulating tumour cells (CTCs) from blood in a microfluidic chip. A hydrodynamic sheath flow is designed upstream in the chip to direct the suspension samples to the channel side walls, thus providing a queue to allow DEP-induced lateral displacements. High-throughput continuous cancer cell sorting (maximum flow rate: ~2.4 mL h(-1), linear velocity: ~4 mm s(-1)) is achieved with a sustained 3D lateral DEP (LDEP) particle force normal to the continuous through-flow. This design allows the continuous fractionation of micro/nanosized particles into different downstream subchannels based on the differences in their different critical negative DEP strengths/mobilities. The main advantage of this separation strategy is that increasing the channel length can effectively increase the throughput proportionally. The effective separation of rare cancer cells (<0.001%) from diluted human blood in a handheld chip is demonstrated. An enrichment factor of 10(5) and a recovery rate of ~85% from a 0.001% cancer cell sample are achieved at an optimal flow rate of 20 μL min(-1) passing through a 6 cm long LDEP channel with an appropriate voltage at a frequency of 10 kHz. A higher throughput of 2.4 mL h(-1) is also achieved with a 13 cm long metal-based microchannel. PMID:26085231
Free-Breathing 3D Whole Heart Black Blood Imaging with Motion Sensitized Driven Equilibrium
Srinivasan, Subashini; Hu, Peng; Kissinger, Kraig V.; Goddu, Beth; Goepfert, Lois; Schmidt, Ehud J.; Kozerke, Sebastian; Nezafat, Reza
2012-01-01
Purpose To assess the efficacy and robustness of motion sensitized driven equilibrium (MSDE) for blood suppression in volumetric 3D whole heart cardiac MR. Materials and Methods To investigate the efficacy of MSDE on blood suppression and myocardial SNR loss on different imaging sequences. 7 healthy adult subjects were imaged using 3D ECG-triggered MSDE-prep T1-weighted turbo spin echo (TSE), and spoiled gradient echo (GRE), after optimization of MSDE parameters in a pilot study of 5 subjects. Imaging artifacts, myocardial and blood SNR were assessed. Subsequently, the feasibility of isotropic spatial resolution MSDE-prep black-blood was assessed in 6 subjects. Finally, 15 patients with known or suspected cardiovascular disease were recruited to be imaged using conventional multi-slice 2D DIR TSE imaging sequence and 3D MSDE-prep spoiled GRE. Results The MSDE-prep yields significant blood suppression (75-92%), enabling a volumetric 3D black-blood assessment of the whole heart with significantly improved visualization of the chamber walls. The MSDE-prep also allowed successful acquisition of black-blood images with isotropic spatial resolution. In the patient study, 3D black-blood MSDE-prep and DIR resulted in similar blood suppression in LV and RV walls but the MSDE prep had superior myocardial signal and wall sharpness. Conclusion MSDE-prep allows volumetric black-blood imaging of the heart. PMID:22517477
Real-time multispectral 3-D photoacoustic imaging of blood phantoms
NASA Astrophysics Data System (ADS)
Kosik, Ivan; Carson, Jeffrey J. L.
2013-03-01
Photoacoustic imaging is exquisitely sensitive to blood and can infer blood oxygenation based on multispectral images. In this work we present multispectral real-time 3D photoacoustic imaging of blood phantoms. We used a custom-built 128-channel hemispherical transducer array coupled to two Nd:YAG pumped OPO laser systems synchronized to provide double pulse excitation at 680 nm and 1064 nm wavelengths, all during a triggered series of ultrasound pressure measurements lasting less than 300 μs. The results demonstrated that 3D PAI is capable of differentiating between oxygenated and deoxygenated blood at high speed at mm-level resolution.
Provost, J.; Papadacci, C.; Demene, C.; Gennisson, J-L.; Tanter, M.; Pernot, M.
2016-01-01
Ultrafast Doppler Imaging was introduced as a technique to quantify blood flow in an entire 2-D field of view, expanding the field of application of ultrasound imaging to the highly sensitive anatomical and functional mapping of blood vessels. We have recently developed 3-D Ultrafast Ultrasound Imaging, a technique that can produce thousands of ultrasound volumes per second, based on three-dimensional plane and diverging wave emissions, and demonstrated its clinical feasibility in human subjects in vivo. In this study, we show that non-invasive 3-D Ultrafast Power Doppler, Pulsed Doppler, and Color Doppler Imaging can be used to perform quantitative imaging of blood vessels in humans when using coherent compounding of three-dimensional tilted plane waves. A customized, programmable, 1024-channel ultrasound system was designed to perform 3-D Ultrafast Imaging. Using a 32X32, 3-MHz matrix phased array (Vermon, France), volumes were beamformed by coherently compounding successive tilted plane wave emissions. Doppler processing was then applied in a voxel-wise fashion. 3-D Ultrafast Power Doppler Imaging was first validated by imaging Tygon tubes of varying diameter and its in vivo feasibility was demonstrated by imaging small vessels in the human thyroid. Simultaneous 3-D Color and Pulsed Doppler Imaging using compounded emissions were also applied in the carotid artery and the jugular vein in one healthy volunteer. PMID:26276956
3d Modelling of Convective Flow In The Rhine Graben
NASA Astrophysics Data System (ADS)
Bächler, D.; Kohl, T.; Rybach, L.
Detailed investigations of the temperature distribution in the Rhine Graben indi- cate regular pattern of thermal anomalies following major north-south striking faults. These anomalies remain unexplained by conventional Rhine Graben studies based on 2D east-west striking sections. First analytical solutions for convective flow in vertical faults are applied for a clearly observable anomalous temperature pattern along ma- jor Rhine Graben faults. By these calculations the fault height, fault aperture, minimal fault permeability and time to convective onset is derived from the observed distances. Since analytical solutions are limited to simple model geometries further improvement was achieved by numerical model simulations, which allow to assume more com- plex initial and boundary conditions. Using the finite volume code TOUGH2 series of anomalies following the same fault were simulated by a 3D numerical model. Fo- cussing on the predominant north-south permeability structure the model consists of a vertical north-south striking fault and surrounding matrix. The fault geometries are based on the analytically predicted fault geometries (aperture=200m, height=3500m) and on the observed temperatures. Comparison of simulation results with observed temperatures shows that the fault is situated between 500 to 600m and 4200m. The fault permeability is taken as 5*10-13m2 and the fluid velocity in the fault is calcu- lated as 10-9 to 10-10 m/s. These results indicate the importance of our considerations since mass flux is much higher in the faults than across them. The minimal age of the anomaly is considered to be 77'000 years, since steady state is reached after this time span. The study proves that the observed temperature anomaly pattern along the gamma fault at Landau can be explained by north-south striking convection systems within fault zones. Similar situations have been found at Soultz. This may be a hint on a general feature of the major north-south striking
RV functional imaging: 3-D echo-derived dynamic geometry and flow field simulations.
Pasipoularides, Ares D; Shu, Ming; Womack, Michael S; Shah, Ashish; Von Ramm, Olaf; Glower, Donald D
2003-01-01
We describe a novel functional imaging approach for quantitative analysis of right ventricular (RV) blood flow patterns in specific experimental animals (or humans) using real-time, three-dimensional (3-D) echocardiography (RT3D). The method is independent of the digital imaging modality used. It comprises three parts. First, a semiautomated segmentation aided by intraluminal contrast medium locates the RV endocardial surface. Second, a geometric scheme for dynamic RV chamber reconstruction applies a time interpolation procedure to the RT3D data to quantify wall geometry and motion at 400 Hz. A volumetric prism method validated the dynamic geometric reconstruction against simultaneous sonomicrometric canine measurements. Finally, the RV endocardial border motion information is used for mesh generation on a computational fluid dynamics solver to simulate development of the early RV diastolic inflow field. Boundary conditions (tessellated endocardial surface nodal velocities) for the solver are directly derived from the endocardial geometry and motion information. The new functional imaging approach may yield important kinematic information on the distribution of instantaneous velocities in the RV diastolic flow field of specific normal or diseased hearts. PMID:12388220
Label-free characterization of white blood cells by measuring 3D refractive index maps
Yoon, Jonghee; Kim, Kyoohyun; Park, HyunJoo; Choi, Chulhee; Jang, Seongsoo; Park, YongKeun
2015-01-01
The characterization of white blood cells (WBCs) is crucial for blood analyses and disease diagnoses. However, current standard techniques rely on cell labeling, a process which imposes significant limitations. Here we present three-dimensional (3D) optical measurements and the label-free characterization of mouse WBCs using optical diffraction tomography. 3D refractive index (RI) tomograms of individual WBCs are constructed from multiple two-dimensional quantitative phase images of samples illuminated at various angles of incidence. Measurements of the 3D RI tomogram of WBCs enable the separation of heterogeneous populations of WBCs using quantitative morphological and biochemical information. Time-lapse tomographic measurements also provide the 3D trajectory of micrometer-sized beads ingested by WBCs. These results demonstrate that optical diffraction tomography can be a useful and versatile tool for the study of WBCs. PMID:26504637
Skalak, R; Keller, S R; Secomb, T W
1981-05-01
The historical development of the mechanics of blood flow can be traced from ancient times, to Leonardo da Vinci and Leonhard Euler and up to the present times with increasing biological knowledge and mathematical analysis. In the last two decades, quantitative and numerical methods have steadily given more complete and precise understanding. In the arterial system wave propagation computations based on nonlinear one-dimensional modeling have given the best representation of pulse wave propagation. In the veins, the theory of unsteady flow in collapsible tubes has recently been extensively developed. In the last decade, progress has been made in describing the blood flow at junctions, through stenoses, in bends and in capillary blood vessels. The rheological behavior of individual red blood cells has been explored. A working model consists of an elastic membrane filled with viscous fluid. This model forms a basis for understanding the viscous and viscoelastic behavior of blood. PMID:7024641
Effects of Presence, Copresence, and Flow on Learning Outcomes in 3D Learning Spaces
ERIC Educational Resources Information Center
Hassell, Martin D.; Goyal, Sandeep; Limayem, Moez; Boughzala, Imed
2012-01-01
The level of satisfaction and effectiveness of 3D virtual learning environments were examined. Additionally, 3D virtual learning environments were compared with face-to-face learning environments. Students that experienced higher levels of flow and presence also experienced more satisfaction but not necessarily more effectiveness with 3D virtual…
Quasi 3D modeling of water flow in vadose zone and groundwater
NASA Astrophysics Data System (ADS)
Kuznetsov, M.; Yakirevich, A.; Pachepsky, Y. A.; Sorek, S.; Weisbrod, N.
2012-07-01
SummaryThe complexity of subsurface flow systems calls for a variety of concepts leading to the multiplicity of simplified flow models. One habitual simplification is based on the assumption that lateral flow and transport in unsaturated zone are not significant unless the capillary fringe is involved. In such cases the flow and transport in the unsaturated zone above groundwater level can be simulated as a 1D phenomenon, whereas the flow and transport through groundwater are viewed as 2D or 3D phenomena. A new approach for a numerical scheme for 3D variably saturated flow using quasi 3D Richards' equation and finite difference scheme is presented. The corresponding numerical algorithm and the QUASI-3D computer code were developed. Results of the groundwater level simulations were compared with transient laboratory experimental data for 2D data constant-flux infiltration, quasi-3D HYDRUS-MODFLOW numerical model and a FULL-3D numerical model using Richards' equation. Hypothetical 3D examples of infiltration, pumping and groundwater mound dissipation for different spatial-time scales are presented. Water flow simulation for the Alto Piura aquifer (Peru) demonstrates the QUASI-3D model application at the regional scale. Computationally the QUASI-3D code was found to be more efficient by an order of 10-300%, while being accurate with respect to the benchmark fully 3D variable saturation code, when the capillary fringe was considered.
Recent Enhancements to USM3D Unstructured Flow Solver for Unsteady Flows
NASA Technical Reports Server (NTRS)
Pandya, Mohagna J.; Frink, Neal T.; Abdol-Hamid, Khaled S.; Chung, James J.
2004-01-01
The NASA USM3D unstructured flow solver is undergoing extensions to address dynamic flow problems in support of NASA and NAVAIR efforts to study the applicability of Computational Fluid Dynamics tools for the prediction of aircraft stability and control characteristics. The initial extensions reported herein include two second-order time stepping schemes, Detached-Eddy Simulation, and grid motion. This paper reports the initial code verification and validation assessment of the dynamic flow capabilities of USM3D. The cases considered are the classic inviscid shock-tube problem, low Reynolds number wake shedding from a NACA 0012 airfoil, high Reynolds number DES-based wake shedding from a 4-to-1 length-to-diameter cylinder, and forced pitch oscillation of a NACA 0012 airfoil with inviscid and turbulent flow.
Development of an accurate 3D blood vessel searching system using NIR light
NASA Astrophysics Data System (ADS)
Mizuno, Yoshifumi; Katayama, Tsutao; Nakamachi, Eiji
2010-02-01
Health monitoring system (HMS) and drug delivery system (DDS) require accurate puncture by needle for automatic blood sampling. In this study, we develop a miniature and high accurate automatic 3D blood vessel searching system. The size of detecting system is 40x25x10 mm. Our searching system use Near-Infrared (NIR) LEDs, CMOS camera modules and image processing units. We employ the stereo method for searching system to determine 3D blood vessel location. Blood vessel visualization system adopts hemoglobin's absorption characterization of NIR light. NIR LED is set behind the finger and it irradiates Near Infrared light for the finger. CMOS camera modules are set in front of the finger and it captures clear blood vessel images. Two dimensional location of the blood vessel is detected by luminance distribution of the image and its depth is calculated by the stereo method. 3D blood vessel location is automatically detected by our image processing system. To examine the accuracy of our detecting system, we carried out experiments using finger phantoms with blood vessel diameters, 0.5, 0.75, 1.0mm, at the depths, 0.5 ~ 2.0 mm, under the artificial tissue surface. Experimental results of depth obtained by our detecting system showed good agreements with given depths, and the availability of this system is confirmed.
Brien, Dianne L.; Reid, Mark E.
2007-01-01
Landslides are a common problem on coastal bluffs throughout the world. Along the coastal bluffs of the Puget Sound in Seattle, Washington, landslides range from small, shallow failures to large, deep-seated landslides. Landslides of all types can pose hazards to human lives and property, but deep-seated landslides are of significant concern because their large areal extent can cause extensive property damage. Although many geomorphic processes shape the coastal bluffs of Seattle, we focus on large (greater than 3,000 m3), deepseated, rotational landslides that occur on the steep bluffs along Puget Sound. Many of these larger failures occur in advance outwash deposits of the Vashon Drift (Qva); some failures extend into the underlying Lawton Clay Member of the Vashon Drift (Qvlc). The slope stability of coastal bluffs is controlled by the interplay of three-dimensional (3-D) variations in gravitational stress, strength, and pore-water pressure. We assess 3-D slope-stability using SCOOPS (Reid and others, 2000), a computer program that allows us to search a high-resolution digital-elevation model (DEM) to quantify the relative stability of all parts of the landscape by computing the stability and volume of thousands of potential spherical failures. SCOOPS incorporates topography, 3-D strength variations, and 3-D pore pressures. Initially, we use our 3-D analysis methods to examine the effects of topography and geology by using heterogeneous material properties, as defined by stratigraphy, without pore pressures. In this scenario, the least-stable areas are located on the steepest slopes, commonly in Qva or Qvlc. However, these locations do not agree well with observations of deep-seated landslides. Historically, both shallow colluvial landslides and deep-seated landslides have been observed near the contact between Qva and Qvlc, and commonly occur in Qva. The low hydraulic conductivity of Qvlc impedes ground-water flow, resulting in elevated pore pressures at the
Navier-Stokes solutions for rotating 3-D duct flows
NASA Astrophysics Data System (ADS)
Srivastava, B. N.
1988-07-01
This paper deals with the computation of three-dimensional viscous turbulent flow in a rotating rectangular duct of low aspect ratio using thin-layer Navier-Stokes equations. Scalar form of an approximate factorization implicit scheme along with a modified q-omega turbulence model has been utilized for mean flow predictions. The predicted mean flow behavior has been favorably compared with the experimental data for mean axial velocity, channel pressure and cross-flow velocities at a flow Mach number of 0.05 and a rotational speed of 300 rpm.
Verification of internal flow analyses in complex 3-D geometries
NASA Astrophysics Data System (ADS)
Choi, S. K.; Buggeln, R. C.
1992-11-01
Analysis of internal flow in advanced rocket propulsion systems is complicated by hardware geometry, high Reynolds numbers, rotation, high frequency phenomena, and near incompressibility. Typical of such a problem is the Space Shuttle Main Engine (SSME) hot gas manifold (HGM). Previous analyses of flow in the SSME HGM have been compared to air flow data and found to be inaccurate with respect to system losses, outer wall static pressures, and transfer duct environments. Such discrepancies could arise from flow measurement methodology, low order algorithms, turbulence modeling, and/or inadequate grid resolution. The objective of this work is to compare internal flow computational analyses to LDV flow measurements for the MSFC HGM pilot model configuration using two grids of different node density in the near wall region. Grids were generated with the EAGLE grid generator and calculations were made with the SRA MINT code. The calculated results were compared with HGM experimental data obtained in the MSFC water flow facility.
3D topographic correction of the BSR heat flow and detection of focused fluid flow
NASA Astrophysics Data System (ADS)
He, Tao; Li, Hong-Lin; Zou, Chang-Chun
2014-06-01
The bottom-simulating reflector (BSR) is a seismic indicator of the bottom of a gas hydrate stability zone. Its depth can be used to calculate the seafloor surface heat flow. The calculated BSR heat flow variations include disturbances from two important factors: (1) seafloor topography, which focuses the heat flow over regions of concave topography and defocuses it over regions of convex topography, and (2) the focused warm fluid flow within the accretionary prism coming from depths deeper than BSR. The focused fluid flow can be detected if the contribution of the topography to the BSR heat flow is removed. However, the analytical equation cannot solve the topographic effect at complex seafloor regions. We prove that 3D finite element method can model the topographic effect on the regional background heat flow with high accuracy, which can then be used to correct the topographic effect and obtain the BSR heat flow under the condition of perfectly flat topography. By comparing the corrected BSR heat flow with the regional background heat flow, focused fluid flow regions can be detected that are originally too small and cannot be detected using present-day equipment. This method was successfully applied to the midslope region of northern Cascadia subducting margin. The results suggest that the Cucumber Ridge and its neighboring area are positive heat flow anomalies, about 10%-20% higher than the background heat flow after 3D topographic correction. Moreover, the seismic imaging associated the positive heat flow anomaly areas with seabed fracture-cavity systems. This suggests flow of warm gas-carrying fluids along these high-permeability pathways, which could result in higher gas hydrate concentrations.
Implementation of Flow Tripping Capability in the USM3D Unstructured Flow Solver
NASA Technical Reports Server (NTRS)
Pandya, Mohagna J.; Abdol-Hamid, Khaled S.; Campbell, Richard L.; Frink, Neal T.
2006-01-01
A flow tripping capability is added to an established NASA tetrahedral unstructured parallel Navier-Stokes flow solver, USM3D. The capability is based on prescribing an appropriate profile of turbulence model variables to energize the boundary layer in a plane normal to a specified trip region on the body surface. We demonstrate this approach using the k-epsilon two-equation turbulence model of USM3D. Modification to the solution procedure primarily consists of developing a data structure to identify all unstructured tetrahedral grid cells located in the plane normal to a specified surface trip region and computing a function based on the mean flow solution to specify the modified profile of the turbulence model variables. We leverage this data structure and also show an adjunct approach that is based on enforcing a laminar flow condition on the otherwise fully turbulent flow solution in user-specified region. The latter approach is applied for the solutions obtained using other one-and two-equation turbulence models of USM3D. A key ingredient of the present capability is the use of a graphical user-interface tool PREDISC to define a trip region on the body surface in an existing grid. Verification of the present modifications is demonstrated on three cases, namely, a flat plate, the RAE2822 airfoil, and the DLR F6 wing-fuselage configuration.
Implementation of Flow Tripping Capability in the USM3D Unstructured Flow Solver
NASA Technical Reports Server (NTRS)
Pandya, Mohagna J.; Abdol-Harrid, Khaled S.; Campbell, Richard L.; Frink, Neal T.
2006-01-01
A flow tripping capability is added to an established NASA tetrahedral unstructured parallel Navier-Stokes flow solver, USM3D. The capability is based on prescribing an appropriate profile of turbulence model variables to energize the boundary layer in a plane normal to a specified trip region on the body surface. We demonstrate this approach using the k-e two-equation turbulence model of USM3D. Modification to the solution procedure primarily consists of developing a data structure to identify all unstructured tetrahedral grid cells located in the plane normal to a specified surface trip region and computing a function based on the mean flow solution to specify the modified profile of the turbulence model variables. We leverage this data structure and also show an adjunct approach that is based on enforcing a laminar flow condition on the otherwise fully turbulent flow solution in user specified region. The latter approach is applied for the solutions obtained using other one- and two-equation turbulence models of USM3D. A key ingredient of the present capability is the use of a graphical user-interface tool PREDISC to define a trip region on the body surface in an existing grid. Verification of the present modifications is demonstrated on three cases, namely, a flat plate, the RAE2822 airfoil, and the DLR F6 wing-fuselage configuration.
Gas flow environment and heat transfer nonrotating 3D program
NASA Technical Reports Server (NTRS)
Schulz, R. J.
1982-01-01
A complete set of benchmark quality data for the flow and heat transfer within a large rectangular turning duct is provided. These data are to be used to evaluate, and verify, three-dimensional internal viscous flow models and computational codes. The analytical contract objective is to select a computational code and define the capabilities of this code to predict the experimental results obtained. Details of the proper code operation will be defined and improvements to the code modeling capabilities will be formulated. Internal flow in a large rectangular cross-sectioned 90 deg. bend turning duct was studied. The duct construction was designed to allow detailed measurements to be made for the following three duct wall conditions: (1) an isothermal wall with isothermal flow; (2) an adiabatic wall with convective heat transfer by mixing between an unheated surrounding flow; and (3) an isothermal wall with heat transfer from a uniformly hot inlet flow.
A 3D velocimetry study of the flow through prosthetic heart valves
NASA Astrophysics Data System (ADS)
Ledesma, R.; Zenit, R.; Pulos, G.; Sanchez, E.; Juarez, A.
2006-11-01
Blood damage commonly appears in medical valve prothesis. It is a mayor concern for the designers and surgeons. It is well known that this damage and other complications result from the modified fluid dynamics through the replacement valve. To evaluate the performance of prosthetic heart valves, it is necessary to study the flow through them. To conduct this study , we have built a flow channel that emulates cardiac conditions and allows optical access such that a 3D-PIV velocimetry system could be used. The experiments are aimed to reconstruct the downstream structure of the flow through a mechanical and a bio-material tricuspid heart valve prothesis. Preliminary results show that the observed coherent structures can be related with haemolysis and trombosis, illnesses commonly found in valve prothesis recipients. The mean flow, the levels of strain rate and the turbulence intensity generated by the valves can also be directly related to blood damage. In general, bio-material made valves tend to reduce these complications.
3-D Particle Tracking Velocimetry: Development and Applications in Small Scale Flows
NASA Astrophysics Data System (ADS)
Tien, Wei-Hsin
valve design a shorter leaflets is advantageous because it prevents flow stasis and reduce the energy loss. Valve pairing tests showed that an orthogonal configuration of the valve pair result in a complicated 3-D flow around the valve, which can increase the mixing of the blood flow and prevent reversed flow in between the valves. The tests of different valve separation distance showed that the coupling effect of two valves was weakened as the separation distance increased, suggesting the existence of a separation distance between the two valves to maximize the coupling effect and keep the flow structure stable.
Engineering an in vitro air-blood barrier by 3D bioprinting
Horváth, Lenke; Umehara, Yuki; Jud, Corinne; Blank, Fabian; Petri-Fink, Alke; Rothen-Rutishauser, Barbara
2015-01-01
Intensive efforts in recent years to develop and commercialize in vitro alternatives in the field of risk assessment have yielded new promising two- and three dimensional (3D) cell culture models. Nevertheless, a realistic 3D in vitro alveolar model is not available yet. Here we report on the biofabrication of the human air-blood tissue barrier analogue composed of an endothelial cell, basement membrane and epithelial cell layer by using a bioprinting technology. In contrary to the manual method, we demonstrate that this technique enables automatized and reproducible creation of thinner and more homogeneous cell layers, which is required for an optimal air-blood tissue barrier. This bioprinting platform will offer an excellent tool to engineer an advanced 3D lung model for high-throughput screening for safety assessment and drug efficacy testing. PMID:25609567
NASA Astrophysics Data System (ADS)
Nilsson, G. E.
1997-01-01
The operating principles of Laser Doppler Perfusion Imaging (LDPI) for visualization of the tissue blood perfusion are explained. Using this emerging technology skin perfusion has been investigated in healthy volunteers and in patients with various conditions that affect skin blood flow. LDPI is anticipated to be particularly useful in evaluation of peripheral circulation in diabetics, as an objective tool in irritancy patch testing, assessment of burnt skin and visualization of spot-wise hyperperfusion in breast skin in association with carcinoma.
Gas flow environmental and heat transfer nonrotating 3D program
NASA Technical Reports Server (NTRS)
Geil, T.; Steinhoff, J.
1983-01-01
A complete set of benchmark quality data for the flow and heat transfer within a large rectangular turning duct is being compiled. These data will be used to evaluate and verify three dimensional internal viscous flow models and computational codes. The analytical objective is to select such a computational code and define the capabilities of this code to predict the experimental results. Details of the proper code operation will be defined and improvements to the code modeling capabilities will be formulated.
New 3-D flow interpolation method on moving ADCP data
NASA Astrophysics Data System (ADS)
Tsubaki, R.; Kawahara, Y.; Muto, Y.; Fujita, I.
2012-05-01
A simple but accurate interpolation procedure for obtaining the three-dimensional distribution of three-component velocity data, from moving acoustic doppler current profiler (ADCP) observation data, is proposed. For understanding actual flow structure within a river with complex bathymetry, the three-dimensional mean velocity field provides a basic picture of the flow. For obtaining the three-dimensional distribution of three-component velocity data, in this work, anisotropic gridding was introduced in order to remove the random component of measured velocity data caused by the turbulence of the flow and measurement error. A continuity correction based on the pressure equation was used to reduce both random and systematic errors. The accuracy of the developed method was evaluated using three-dimensional flow simulation data from a detached-eddy simulation (DES). By using the procedure developed, the complex flow structure surrounding the spur dikes section in the Uji River was successfully visualized and explored. The proposed method shows superiorities in both accuracy and consistency for the interpolated velocity field, as compared to the kriging and inverse-distance weighted (IDW) methods.
Segmented Domain Decomposition Multigrid For 3-D Turbomachinery Flows
NASA Technical Reports Server (NTRS)
Celestina, M. L.; Adamczyk, J. J.; Rubin, S. G.
2001-01-01
A Segmented Domain Decomposition Multigrid (SDDMG) procedure was developed for three-dimensional viscous flow problems as they apply to turbomachinery flows. The procedure divides the computational domain into a coarse mesh comprised of uniformly spaced cells. To resolve smaller length scales such as the viscous layer near a surface, segments of the coarse mesh are subdivided into a finer mesh. This is repeated until adequate resolution of the smallest relevant length scale is obtained. Multigrid is used to communicate information between the different grid levels. To test the procedure, simulation results will be presented for a compressor and turbine cascade. These simulations are intended to show the ability of the present method to generate grid independent solutions. Comparisons with data will also be presented. These comparisons will further demonstrate the usefulness of the present work for they allow an estimate of the accuracy of the flow modeling equations independent of error attributed to numerical discretization.
Lattice Boltzmann Model of 3D Multiphase Flow in Artery Bifurcation Aneurysm Problem
Abas, Aizat; Mokhtar, N. Hafizah; Ishak, M. H. H.; Abdullah, M. Z.; Ho Tian, Ang
2016-01-01
This paper simulates and predicts the laminar flow inside the 3D aneurysm geometry, since the hemodynamic situation in the blood vessels is difficult to determine and visualize using standard imaging techniques, for example, magnetic resonance imaging (MRI). Three different types of Lattice Boltzmann (LB) models are computed, namely, single relaxation time (SRT), multiple relaxation time (MRT), and regularized BGK models. The results obtained using these different versions of the LB-based code will then be validated with ANSYS FLUENT, a commercially available finite volume- (FV-) based CFD solver. The simulated flow profiles that include velocity, pressure, and wall shear stress (WSS) are then compared between the two solvers. The predicted outcomes show that all the LB models are comparable and in good agreement with the FVM solver for complex blood flow simulation. The findings also show minor differences in their WSS profiles. The performance of the parallel implementation for each solver is also included and discussed in this paper. In terms of parallelization, it was shown that LBM-based code performed better in terms of the computation time required. PMID:27239221
Lattice Boltzmann Model of 3D Multiphase Flow in Artery Bifurcation Aneurysm Problem.
Abas, Aizat; Mokhtar, N Hafizah; Ishak, M H H; Abdullah, M Z; Ho Tian, Ang
2016-01-01
This paper simulates and predicts the laminar flow inside the 3D aneurysm geometry, since the hemodynamic situation in the blood vessels is difficult to determine and visualize using standard imaging techniques, for example, magnetic resonance imaging (MRI). Three different types of Lattice Boltzmann (LB) models are computed, namely, single relaxation time (SRT), multiple relaxation time (MRT), and regularized BGK models. The results obtained using these different versions of the LB-based code will then be validated with ANSYS FLUENT, a commercially available finite volume- (FV-) based CFD solver. The simulated flow profiles that include velocity, pressure, and wall shear stress (WSS) are then compared between the two solvers. The predicted outcomes show that all the LB models are comparable and in good agreement with the FVM solver for complex blood flow simulation. The findings also show minor differences in their WSS profiles. The performance of the parallel implementation for each solver is also included and discussed in this paper. In terms of parallelization, it was shown that LBM-based code performed better in terms of the computation time required. PMID:27239221
Atherosclerosis imaging using 3D black blood TSE SPACE vs 2D TSE
Wong, Stephanie K; Mobolaji-Iawal, Motunrayo; Arama, Leron; Cambe, Joy; Biso, Sylvia; Alie, Nadia; Fayad, Zahi A; Mani, Venkatesh
2014-01-01
AIM: To compare 3D Black Blood turbo spin echo (TSE) sampling perfection with application-optimized contrast using different flip angle evolution (SPACE) vs 2D TSE in evaluating atherosclerotic plaques in multiple vascular territories. METHODS: The carotid, aortic, and femoral arterial walls of 16 patients at risk for cardiovascular or atherosclerotic disease were studied using both 3D black blood magnetic resonance imaging SPACE and conventional 2D multi-contrast TSE sequences using a consolidated imaging approach in the same imaging session. Qualitative and quantitative analyses were performed on the images. Agreement of morphometric measurements between the two imaging sequences was assessed using a two-sample t-test, calculation of the intra-class correlation coefficient and by the method of linear regression and Bland-Altman analyses. RESULTS: No statistically significant qualitative differences were found between the 3D SPACE and 2D TSE techniques for images of the carotids and aorta. For images of the femoral arteries, however, there were statistically significant differences in all four qualitative scores between the two techniques. Using the current approach, 3D SPACE is suboptimal for femoral imaging. However, this may be due to coils not being optimized for femoral imaging. Quantitatively, in our study, higher mean total vessel area measurements for the 3D SPACE technique across all three vascular beds were observed. No significant differences in lumen area for both the right and left carotids were observed between the two techniques. Overall, a significant-correlation existed between measures obtained between the two approaches. CONCLUSION: Qualitative and quantitative measurements between 3D SPACE and 2D TSE techniques are comparable. 3D-SPACE may be a feasible approach in the evaluation of cardiovascular patients. PMID:24876923
Laser direct writing 3D structures for microfluidic channels: flow meter and mixer
NASA Astrophysics Data System (ADS)
Lin, Chih-Lang; Liu, Yi-Jui; Lin, Zheng-Da; Wu, Bo-Long; Lee, Yi-Hsiung; Shin, Chow-Shing; Baldeck, Patrice L.
2015-03-01
The 3D laser direct-writing technology is aimed at the modeling of arbitrary three-dimensional (3D) complex microstructures by scanning a laser-focusing point along predetermined trajectories. Through the perspective technique, the details of designed 3D structures can be properly fabricated in a microchannel. This study introduces a direct reading flow meter and a 3D passive mixer fabricated by laser direct writing for microfluidic applications. The flow meter consists of two rod-shaped springs, a pillar, an anchor, and a wedge-shaped indicator, installed inside a microfluidic channel. The indicator is deflected by the flowing fluid while restrained by the spring to establish an equilibrium indication according to the flow rate. The measurement is readily carried out by optical microscopy observation. The 3D passive Archimedes-screw-shaped mixer is designed to disturb the laminar flow 3D direction for enhancing the mixing efficiency. The simulation results indicate that the screw provides 3D disturbance of streamlines in the microchannel. The mixing demonstration for fluids flowing in the micrchannel approximately agrees with the simulation result. Thanks to the advantage of the laser direct writing technology, this study performs the ingenious applications of 3D structures for microchannels.
Modeling blood flow heterogeneity.
King, R B; Raymond, G M; Bassingthwaighte, J B
1996-01-01
It has been known for some time that regional blood flows within an organ are not uniform. Useful measures of heterogeneity of regional blood flows are the standard deviation and coefficient of variation or relative dispersion of the probability density function (PDF) of regional flows obtained from the regional concentrations of tracers that are deposited in proportion to blood flow. When a mathematical model is used to analyze dilution curves after tracer solute administration, for many solutes it is important to account for flow heterogeneity and the wide range of transit times through multiple pathways in parallel. Failure to do so leads to bias in the estimates of volumes of distribution and membrane conductances. Since in practice the number of paths used should be relatively small, the analysis is sensitive to the choice of the individual elements used to approximate the distribution of flows or transit times. Presented here is a method for modeling heterogeneous flow through an organ using a scheme that covers both the high flow and long transit time extremes of the flow distribution. With this method, numerical experiments are performed to determine the errors made in estimating parameters when flow heterogeneity is ignored, in both the absence and presence of noise. The magnitude of the errors in the estimates depends upon the system parameters, the amount of flow heterogeneity present, and whether the shape of the input function is known. In some cases, some parameters may be estimated to within 10% when heterogeneity is ignored (homogeneous model), but errors of 15-20% may result, even when the level of heterogeneity is modest. In repeated trials in the presence of 5% noise, the mean of the estimates was always closer to the true value with the heterogeneous model than when heterogeneity was ignored, but the distributions of the estimates from the homogeneous and heterogeneous models overlapped for some parameters when outflow dilution curves were
Segmentation of Blood Vessels and 3D Representation of CMR Image
NASA Astrophysics Data System (ADS)
Jiji, G. W.
2013-06-01
Current cardiac magnetic resonance imaging (CMR) technology allows the determination of patient-individual coronary tree structure, detection of infarctions, and assessment of myocardial perfusion. The purpose of this work is to segment heart blood vessels and visualize it in 3D. In this work, 3D visualisation of vessel was performed into four phases. The first step is to detect the tubular structures using multiscale medialness function, which distinguishes tube-like structures from and other structures. Second step is to extract the centrelines of the tubes. From the centreline radius the cylindrical tube model is constructed. The third step is segmentation of the tubular structures. The cylindrical tube model is used in segmentation process. Fourth step is to 3D representation of the tubular structure using Volume . The proposed approach is applied to 10 datasets of patients from the clinical routine and tested the results with radiologists.
Approximate analytic solutions to 3D unconfined groundwater flow within regional 2D models
NASA Astrophysics Data System (ADS)
Luther, K.; Haitjema, H. M.
2000-04-01
We present methods for finding approximate analytic solutions to three-dimensional (3D) unconfined steady state groundwater flow near partially penetrating and horizontal wells, and for combining those solutions with regional two-dimensional (2D) models. The 3D solutions use distributed singularities (analytic elements) to enforce boundary conditions on the phreatic surface and seepage faces at vertical wells, and to maintain fixed-head boundary conditions, obtained from the 2D model, at the perimeter of the 3D model. The approximate 3D solutions are analytic (continuous and differentiable) everywhere, including on the phreatic surface itself. While continuity of flow is satisfied exactly in the infinite 3D flow domain, water balance errors can occur across the phreatic surface.
Digital holography for recovering 3D shape of red blood cells
NASA Astrophysics Data System (ADS)
Memmolo, P.; Miccio, L.; Merola, F.; Gennari, O.; Netti, P.; Ferraro, Pietro
2015-07-01
Full morphometric data analysis and 3D rendering of Red Blood Cells (RBCs) is provided by means of Digital Holography (DH) in combination with Optical Tweezers (OT). The proposed method is compared with a geometrical model of RBC in order to evaluate its accuracy and tested for many kinds of RBCs, from healthy ones with double-concavity to that with abnormal shapes. Applications in diagnostics are foreseen.
3D-Flow processor for a programmable Level-1 trigger (feasibility study)
Crosetto, D.
1992-10-01
A feasibility study has been made to use the 3D-Flow processor in a pipelined programmable parallel processing architecture to identify particles such as electrons, jets, muons, etc., in high-energy physics experiments.
NASA Technical Reports Server (NTRS)
Parikh, Paresh; Pirzadeh, Shahyar; Loehner, Rainald
1990-01-01
A set of computer programs for 3-D unstructured grid generation, fluid flow calculations, and flow field visualization was developed. The grid generation program, called VGRID3D, generates grids over complex configurations using the advancing front method. In this method, the point and element generation is accomplished simultaneously, VPLOT3D is an interactive, menudriven pre- and post-processor graphics program for interpolation and display of unstructured grid data. The flow solver, VFLOW3D, is an Euler equation solver based on an explicit, two-step, Taylor-Galerkin algorithm which uses the Flux Corrected Transport (FCT) concept for a wriggle-free solution. Using these programs, increasingly complex 3-D configurations of interest to aerospace community were gridded including a complete Space Transportation System comprised of the space-shuttle orbitor, the solid-rocket boosters, and the external tank. Flow solutions were obtained on various configurations in subsonic, transonic, and supersonic flow regimes.
Quaini, A.; Canic, S.; Glowinski, R.; Igo, S.; Hartley, C.J.; Zoghbi, W.; Little, S.
2011-01-01
This work presents a validation of a fluid-structure interaction computational model simulating the flow conditions in an in vitro mock heart chamber modeling mitral valve regurgitation during the ejection phase during which the trans-valvular pressure drop and valve displacement are not as large. The mock heart chamber was developed to study the use of 2D and 3D color Doppler techniques in imaging the clinically relevant complex intra-cardiac flow events associated with mitral regurgitation. Computational models are expected to play an important role in supporting, refining, and reinforcing the emerging 3D echocardiographic applications. We have developed a 3D computational fluid-structure interaction algorithm based on a semi-implicit, monolithic method, combined with an arbitrary Lagrangian-Eulerian approach to capture the fluid domain motion. The mock regurgitant mitral valve corresponding to an elastic plate with a geometric orifice, was modeled using 3D elasticity, while the blood flow was modeled using the 3D Navier-Stokes equations for an incompressible, viscous fluid. The two are coupled via the kinematic and dynamic conditions describing the two-way coupling. The pressure, the flow rate, and orifice plate displacement were measured and compared with numerical simulation results. In-line flow meter was used to measure the flow, pressure transducers were used to measure the pressure, and a Doppler method developed by one of the authors was used to measure the axial displacement of the orifice plate. The maximum recorded difference between experiment and numerical simulation for the flow rate was 4%, the pressure 3.6%, and for the orifice displacement 15%, showing excellent agreement between the two. PMID:22138194
Quaini, A; Canic, S; Glowinski, R; Igo, S; Hartley, C J; Zoghbi, W; Little, S
2012-01-10
This work presents a validation of a fluid-structure interaction computational model simulating the flow conditions in an in vitro mock heart chamber modeling mitral valve regurgitation during the ejection phase during which the trans-valvular pressure drop and valve displacement are not as large. The mock heart chamber was developed to study the use of 2D and 3D color Doppler techniques in imaging the clinically relevant complex intra-cardiac flow events associated with mitral regurgitation. Computational models are expected to play an important role in supporting, refining, and reinforcing the emerging 3D echocardiographic applications. We have developed a 3D computational fluid-structure interaction algorithm based on a semi-implicit, monolithic method, combined with an arbitrary Lagrangian-Eulerian approach to capture the fluid domain motion. The mock regurgitant mitral valve corresponding to an elastic plate with a geometric orifice, was modeled using 3D elasticity, while the blood flow was modeled using the 3D Navier-Stokes equations for an incompressible, viscous fluid. The two are coupled via the kinematic and dynamic conditions describing the two-way coupling. The pressure, the flow rate, and orifice plate displacement were measured and compared with numerical simulation results. In-line flow meter was used to measure the flow, pressure transducers were used to measure the pressure, and a Doppler method developed by one of the authors was used to measure the axial displacement of the orifice plate. The maximum recorded difference between experiment and numerical simulation for the flow rate was 4%, the pressure 3.6%, and for the orifice displacement 15%, showing excellent agreement between the two. PMID:22138194
NASA Astrophysics Data System (ADS)
Saris, Anne E. C. M.; Nillesen, Maartje M.; Lopata, Richard G. P.; de Korte, Chris L.
2013-03-01
Automated segmentation of 3D echocardiographic images in patients with congenital heart disease is challenging, because the boundary between blood and cardiac tissue is poorly defined in some regions. Cardiologists mentally incorporate movement of the heart, using temporal coherence of structures to resolve ambiguities. Therefore, we investigated the merit of temporal cross-correlation for automated segmentation over the entire cardiac cycle. Optimal settings for maximum cross-correlation (MCC) calculation, based on a 3D cross-correlation based displacement estimation algorithm, were determined to obtain the best contrast between blood and myocardial tissue over the entire cardiac cycle. Resulting envelope-based as well as RF-based MCC values were used as additional external force in a deformable model approach, to segment the left-ventricular cavity in entire systolic phase. MCC values were tested against, and combined with, adaptive filtered, demodulated RF-data. Segmentation results were compared with manually segmented volumes using a 3D Dice Similarity Index (3DSI). Results in 3D pediatric echocardiographic images sequences (n = 4) demonstrate that incorporation of temporal information improves segmentation. The use of MCC values, either alone or in combination with adaptive filtered, demodulated RF-data, resulted in an increase of the 3DSI in 75% of the cases (average 3DSI increase: 0.71 to 0.82). Results might be further improved by optimizing MCC-contrast locally, in regions with low blood-tissue contrast. Reducing underestimation of the endocardial volume due to MCC processing scheme (choice of window size) and consequential border-misalignment, could also lead to more accurate segmentations. Furthermore, increasing the frame rate will also increase MCC-contrast and thus improve segmentation.
Determining 3D Flow Fields via Multi-camera Light Field Imaging
Truscott, Tadd T.; Belden, Jesse; Nielson, Joseph R.; Daily, David J.; Thomson, Scott L.
2013-01-01
In the field of fluid mechanics, the resolution of computational schemes has outpaced experimental methods and widened the gap between predicted and observed phenomena in fluid flows. Thus, a need exists for an accessible method capable of resolving three-dimensional (3D) data sets for a range of problems. We present a novel technique for performing quantitative 3D imaging of many types of flow fields. The 3D technique enables investigation of complicated velocity fields and bubbly flows. Measurements of these types present a variety of challenges to the instrument. For instance, optically dense bubbly multiphase flows cannot be readily imaged by traditional, non-invasive flow measurement techniques due to the bubbles occluding optical access to the interior regions of the volume of interest. By using Light Field Imaging we are able to reparameterize images captured by an array of cameras to reconstruct a 3D volumetric map for every time instance, despite partial occlusions in the volume. The technique makes use of an algorithm known as synthetic aperture (SA) refocusing, whereby a 3D focal stack is generated by combining images from several cameras post-capture 1. Light Field Imaging allows for the capture of angular as well as spatial information about the light rays, and hence enables 3D scene reconstruction. Quantitative information can then be extracted from the 3D reconstructions using a variety of processing algorithms. In particular, we have developed measurement methods based on Light Field Imaging for performing 3D particle image velocimetry (PIV), extracting bubbles in a 3D field and tracking the boundary of a flickering flame. We present the fundamentals of the Light Field Imaging methodology in the context of our setup for performing 3DPIV of the airflow passing over a set of synthetic vocal folds, and show representative results from application of the technique to a bubble-entraining plunging jet. PMID:23486112
Reacting Multi-Species Gas Capability for USM3D Flow Solver
NASA Technical Reports Server (NTRS)
Frink, Neal T.; Schuster, David M.
2012-01-01
The USM3D Navier-Stokes flow solver contributed heavily to the NASA Constellation Project (CxP) as a highly productive computational tool for generating the aerodynamic databases for the Ares I and V launch vehicles and Orion launch abort vehicle (LAV). USM3D is currently limited to ideal-gas flows, which are not adequate for modeling the chemistry or temperature effects of hot-gas jet flows. This task was initiated to create an efficient implementation of multi-species gas and equilibrium chemistry into the USM3D code to improve its predictive capabilities for hot jet impingement effects. The goal of this NASA Engineering and Safety Center (NESC) assessment was to implement and validate a simulation capability to handle real-gas effects in the USM3D code. This document contains the outcome of the NESC assessment.
ODTLES : a model for 3D turbulent flow based on one-dimensional turbulence modeling concepts.
McDermott, Randy; Kerstein, Alan R.; Schmidt, Rodney Cannon
2005-01-01
This report describes an approach for extending the one-dimensional turbulence (ODT) model of Kerstein [6] to treat turbulent flow in three-dimensional (3D) domains. This model, here called ODTLES, can also be viewed as a new LES model. In ODTLES, 3D aspects of the flow are captured by embedding three, mutually orthogonal, one-dimensional ODT domain arrays within a coarser 3D mesh. The ODTLES model is obtained by developing a consistent approach for dynamically coupling the different ODT line sets to each other and to the large scale processes that are resolved on the 3D mesh. The model is implemented computationally and its performance is tested and evaluated by performing simulations of decaying isotropic turbulence, a standard turbulent flow benchmarking problem.
Strategies for Effectively Visualizing a 3D Flow Using Volume Line Integral Convolution
NASA Technical Reports Server (NTRS)
Interrante, Victoria; Grosch, Chester
1997-01-01
This paper discusses strategies for effectively portraying 3D flow using volume line integral convolution. Issues include defining an appropriate input texture, clarifying the distinct identities and relative depths of the advected texture elements, and selectively highlighting regions of interest in both the input and output volumes. Apart from offering insights into the greater potential of 3D LIC as a method for effectively representing flow in a volume, a principal contribution of this work is the suggestion of a technique for generating and rendering 3D visibility-impeding 'halos' that can help to intuitively indicate the presence of depth discontinuities between contiguous elements in a projection and thereby clarify the 3D spatial organization of elements in the flow. The proposed techniques are applied to the visualization of a hot, supersonic, laminar jet exiting into a colder, subsonic coflow.
Retinal blood flow measurement by using optical coherence tomography
NASA Astrophysics Data System (ADS)
Makita, Shuichi; Fabritius, Tapio; Miura, Masahiro; Yatagai, Toyohiko; Yasuno, Yoshiaki
2008-02-01
Quantification of the three-dimensional (3D) retinal vessel structure and blood flow is demonstrated. 3D blood flow distribution is obtained by Doppler optical coherence angiography (D-OCA). Vessel parameters, i.e. diameter, orientation, and position, are determined in an en face vessel image. The Doppler angle is estimated as the angle between the retinal vessel and the incident probing beam in representative cross-sectional flow image which extracted from the 3D flow distribution according to the vessel parameters. Blood flow velocity and volume rate can be quantified with these vessel parameters. The retinal blood flow velocity and volume rate are measured in the retinal vessels around the optic nerve head.
Three-dimensional potential flows from functions of a 3D complex variable
NASA Technical Reports Server (NTRS)
Kelly, Patrick; Panton, Ronald L.; Martin, E. D.
1990-01-01
Potential, or ideal, flow velocities can be found from the gradient of an harmonic function. An ordinary complex valued analytic function can be written as the sum of two real valued functions, both of which are harmonic. Thus, 2D complex valued functions serve as a source of functions that describe two-dimensional potential flows. However, this use of complex variables has been limited to two-dimensions. Recently, a new system of three-dimensional complex variables has been developed at the NASA Ames Research Center. As a step toward application of this theory to the analysis of 3D potential flow, several functions of a three-dimensional complex variable have been investigated. The results for two such functions, the 3D exponential and 3D logarithm, are presented in this paper. Potential flows found from these functions are investigated. Important characteristics of these flows fields are noted.
Correlations of Surface Deformation and 3D Flow Field in a Compliant Wall Turbulent Channel Flow.
NASA Astrophysics Data System (ADS)
Wang, Jin; Zhang, Cao; Katz, Joseph
2015-11-01
This study focuses on the correlations between surface deformation and flow features, including velocity, vorticity and pressure, in a turbulent channel flow over a flat, compliant Polydimethylsiloxane (PDMS) wall. The channel centerline velocity is 2.5 m/s, and the friction Reynolds number is 2.3x103. Analysis is based on simultaneous measurements of the time resolved 3D velocity and surface deformation using tomographic PIV and Mach-Zehnder Interferometry. The volumetric pressure distribution is calculated plane by plane by spatially integrating the material acceleration using virtual boundary, omni-directional method. Conditional sampling based on local high/low pressure and deformation events reveals the primary flow structures causing the deformation. High pressure peaks appear at the interface between sweep and ejection, whereas the negative deformations peaks (dent) appear upstream, under the sweeps. The persistent phase lag between flow and deformations are presumably caused by internal damping within the PDMS. Some of the low pressure peaks and strong ejections are located under the head of hairpin vortices, and accordingly, are associated with positive deformation (bump). Others bumps and dents are correlated with some spanwise offset large inclined quasi-streamwise vortices that are not necessarily associated with hairpins. Sponsored by ONR.
3D-printed devices for continuous-flow organic chemistry
Dragone, Vincenza; Sans, Victor; Rosnes, Mali H; Kitson, Philip J
2013-01-01
Summary We present a study in which the versatility of 3D-printing is combined with the processing advantages of flow chemistry for the synthesis of organic compounds. Robust and inexpensive 3D-printed reactionware devices are easily connected using standard fittings resulting in complex, custom-made flow systems, including multiple reactors in a series with in-line, real-time analysis using an ATR-IR flow cell. As a proof of concept, we utilized two types of organic reactions, imine syntheses and imine reductions, to show how different reactor configurations and substrates give different products. PMID:23766811
The computation of steady 3-D separated flows over aerodynamic bodies at incidence and yaw
NASA Technical Reports Server (NTRS)
Pulliam, T. H.; Pan, D.
1986-01-01
This paper describes the implementation of a general purpose 3-D NS code and its application to simulated 3-D separated vortical flows over aerodynamic bodies. The thin-layer Reynolds-averaged NS equations are solved by an implicit approximate factorization scheme. The pencil data structure enables the code to run on very fine grids using only limited incore memories. Solutions of a low subsonic flow over an inclined ellipsoid are compared with experimental data to validate the code. Transonic flows over a yawed elliptical wing at incidence are computed and separations occurred at different yaw angles are discussed.
Inverse cascades sustained by the transfer rate of angular momentum in a 3D turbulent flow.
López-Caballero, Miguel; Burguete, Javier
2013-03-22
The existence of energy cascades as signatures of conserved magnitudes is one of the universal characteristics of turbulent flows. In homogeneous 3D turbulence, the energy conservation produces a direct cascade from large to small scales, although in 2D, it produces an inverse cascade pointing towards small wave numbers. In this Letter, we present the first evidence of an inverse cascade in a fully developed 3D experimental turbulent flow where the conserved magnitude is the angular momentum. Two counterrotating flows collide in a central region where very large fluctuations are produced, generating a turbulent drag that transfers the external torque between different fluid layers. PMID:25166809
The performance & flow visualization studies of three-dimensional (3-D) wind turbine blade models
NASA Astrophysics Data System (ADS)
Sutrisno, Prajitno, Purnomo, W., Setyawan B.
2016-06-01
Recently, studies on the design of 3-D wind turbine blades have a less attention even though 3-D blade products are widely sold. In contrary, advanced studies in 3-D helicopter blade tip have been studied rigorously. Studies in wind turbine blade modeling are mostly assumed that blade spanwise sections behave as independent two-dimensional airfoils, implying that there is no exchange of momentum in the spanwise direction. Moreover, flow visualization experiments are infrequently conducted. Therefore, a modeling study of wind turbine blade with visualization experiment is needed to be improved to obtain a better understanding. The purpose of this study is to investigate the performance of 3-D wind turbine blade models with backward-forward swept and verify the flow patterns using flow visualization. In this research, the blade models are constructed based on the twist and chord distributions following Schmitz's formula. Forward and backward swept are added to the rotating blades. Based on this, the additional swept would enhance or diminish outward flow disturbance or stall development propagation on the spanwise blade surfaces to give better blade design. Some combinations, i. e., b lades with backward swept, provide a better 3-D favorable rotational force of the rotor system. The performance of the 3-D wind turbine system model is measured by a torque meter, employing Prony's braking system. Furthermore, the 3-D flow patterns around the rotating blade models are investigated by applying "tuft-visualization technique", to study the appearance of laminar, separated, and boundary layer flow patterns surrounding the 3-dimentional blade system.
Park, HyunJoo; Hong, Sung-Hee; Kim, Kyoohyun; Cho, Shin-Hyeong; Lee, Won-Ja; Kim, Youngchan; Lee, Sang-Eun; Park, YongKeun
2015-01-01
Babesia microti causes “emergency” human babesiosis. However, little is known about the alterations in B. microti invaded red blood cells (Bm-RBCs) at the individual cell level. Through quantitative phase imaging techniques based on laser interferometry, we present the simultaneous measurements of structural, chemical, and mechanical modifications in individual mouse Bm-RBCs. 3-D refractive index maps of individual RBCs and in situ parasite vacuoles are imaged, from which total contents and concentration of dry mass are also precisely quantified. In addition, we examine the dynamic membrane fluctuation of Bm-RBCs, which provide information on cell membrane deformability. PMID:26039793
3D cultured immortalized human hepatocytes useful to develop drugs for blood-borne HCV
Aly, Hussein Hassan; Shimotohno, Kunitada; Hijikata, Makoto
2009-02-06
Due to the high polymorphism of natural hepatitis C virus (HCV) variants, existing recombinant HCV replication models have failed to be effective in developing effective anti-HCV agents. In the current study, we describe an in vitro system that supports the infection and replication of natural HCV from patient blood using an immortalized primary human hepatocyte cell line cultured in a three-dimensional (3D) culture system. Comparison of the gene expression profile of cells cultured in the 3D system to those cultured in the existing 2D system demonstrated an up-regulation of several genes activated by peroxisome proliferator-activated receptor alpha (PPAR{alpha}) signaling. Furthermore, using PPAR{alpha} agonists and antagonists, we also analyzed the effect of PPAR{alpha} signaling on the modulation of HCV replication using this system. The 3D in vitro system described in this study provides significant insight into the search for novel anti-HCV strategies that are specific to various strains of HCV.
NASA Astrophysics Data System (ADS)
Yang, L. M.; Shu, C.; Wang, Y.; Sun, Y.
2016-08-01
The sphere function-based gas kinetic scheme (GKS), which was presented by Shu and his coworkers [23] for simulation of inviscid compressible flows, is extended to simulate 3D viscous incompressible and compressible flows in this work. Firstly, we use certain discrete points to represent the spherical surface in the phase velocity space. Then, integrals along the spherical surface for conservation forms of moments, which are needed to recover 3D Navier-Stokes equations, are approximated by integral quadrature. The basic requirement is that these conservation forms of moments can be exactly satisfied by weighted summation of distribution functions at discrete points. It was found that the integral quadrature by eight discrete points on the spherical surface, which forms the D3Q8 discrete velocity model, can exactly match the integral. In this way, the conservative variables and numerical fluxes can be computed by weighted summation of distribution functions at eight discrete points. That is, the application of complicated formulations resultant from integrals can be replaced by a simple solution process. Several numerical examples including laminar flat plate boundary layer, 3D lid-driven cavity flow, steady flow through a 90° bending square duct, transonic flow around DPW-W1 wing and supersonic flow around NACA0012 airfoil are chosen to validate the proposed scheme. Numerical results demonstrate that the present scheme can provide reasonable numerical results for 3D viscous flows.
Implementation of Advanced Two Equation Turbulence Models in the USM3D Unstructured Flow Solver
NASA Technical Reports Server (NTRS)
Wang, Qun-Zhen; Massey, Steven J.; Abdol-Hamid, Khaled S.
2000-01-01
USM3D is a widely-used unstructured flow solver for simulating inviscid and viscous flows over complex geometries. The current version (version 5.0) of USM3D, however, does not have advanced turbulence models to accurately simulate complicated flow. We have implemented two modified versions of the original Jones and Launder k-epsilon "two-equation" turbulence model and the Girimaji algebraic Reynolds stress model in USM3D. Tests have been conducted for three flat plate boundary layer cases, a RAE2822 airfoil and an ONERA M6 wing. The results are compared with those from direct numerical simulation, empirical formulae, theoretical results, and the existing Spalart-Allmaras one-equation model.
Farooqi, Kanwal M; Lengua, Carlos Gonzalez; Weinberg, Alan D; Nielsen, James C; Sanz, Javier
2016-08-01
The method of cardiac magnetic resonance (CMR) three-dimensional (3D) image acquisition and post-processing which should be used to create optimal virtual models for 3D printing has not been studied systematically. Patients (n = 19) who had undergone CMR including both 3D balanced steady-state free precession (bSSFP) imaging and contrast-enhanced magnetic resonance angiography (MRA) were retrospectively identified. Post-processing for the creation of virtual 3D models involved using both myocardial (MS) and blood pool (BP) segmentation, resulting in four groups: Group 1-bSSFP/MS, Group 2-bSSFP/BP, Group 3-MRA/MS and Group 4-MRA/BP. The models created were assessed by two raters for overall quality (1-poor; 2-good; 3-excellent) and ability to identify predefined vessels (1-5: superior vena cava, inferior vena cava, main pulmonary artery, ascending aorta and at least one pulmonary vein). A total of 76 virtual models were created from 19 patient CMR datasets. The mean overall quality scores for Raters 1/2 were 1.63 ± 0.50/1.26 ± 0.45 for Group 1, 2.12 ± 0.50/2.26 ± 0.73 for Group 2, 1.74 ± 0.56/1.53 ± 0.61 for Group 3 and 2.26 ± 0.65/2.68 ± 0.48 for Group 4. The numbers of identified vessels for Raters 1/2 were 4.11 ± 1.32/4.05 ± 1.31 for Group 1, 4.90 ± 0.46/4.95 ± 0.23 for Group 2, 4.32 ± 1.00/4.47 ± 0.84 for Group 3 and 4.74 ± 0.56/4.63 ± 0.49 for Group 4. Models created using BP segmentation (Groups 2 and 4) received significantly higher ratings than those created using MS for both overall quality and number of vessels visualized (p < 0.05), regardless of the acquisition technique. There were no significant differences between Groups 1 and 3. The ratings for Raters 1 and 2 had good correlation for overall quality (ICC = 0.63) and excellent correlation for the total number of vessels visualized (ICC = 0.77). The intra-rater reliability was good for Rater A (ICC = 0.65). Three models were successfully printed
High fidelity digital inline holographic method for 3D flow measurements.
Toloui, Mostafa; Hong, Jiarong
2015-10-19
Among all the 3D optical flow diagnostic techniques, digital inline holographic particle tracking velocimetry (DIH-PTV) provides the highest spatial resolution with low cost, simple and compact optical setups. Despite these advantages, DIH-PTV suffers from major limitations including poor longitudinal resolution, human intervention (i.e. requirement for manually determined tuning parameters during tracer field reconstruction and extraction), limited tracer concentration, and expensive computations. These limitations prevent this technique from being widely used for high resolution 3D flow measurements. In this study, we present a novel holographic particle extraction method with the goal of overcoming all the major limitations of DIH-PTV. The proposed method consists of multiple steps involving 3D deconvolution, automatic signal-to-noise ratio enhancement and thresholding, and inverse iterative particle extraction. The entire method is implemented using GPU-based algorithm to increase the computational speed significantly. Validated with synthetic particle holograms, the proposed method can achieve particle extraction rate above 95% with fake particles less than 3% and maximum position error below 1.6 particle diameter for holograms with particle concentration above 3000 particles/mm^{3}. The applicability of the proposed method for DIH-PTV has been further validated using the experiment of laminar flow in a microchannel and the synthetic tracer flow fields generated using a DNS turbulent channel flow database. Such improvements will substantially enhance the implementation of DIH-PTV for 3D flow measurements and enable the potential commercialization of this technique. PMID:26480377
Time-lapse 3-D seismic imaging of shallow subsurface contaminant flow.
McKenna, J; Sherlock, D; Evans, B
2001-12-01
This paper presents a physical modelling study outlining a technique whereby buoyant contaminant flow within water-saturated unconsolidated sand was remotely monitored utilizing the time-lapse 3-D (TL3-D) seismic response. The controlled temperature and pressure conditions, along with the high level of acquisition repeatability attainable using sandbox physical models, allow the TL3-D seismic response to pore fluid movement to be distinguished from all other effects. TL3-D seismic techniques are currently being developed to monitor hydrocarbon reserves within producing reservoirs in an endeavour to improve overall recovery. However, in many ways, sandbox models under atmospheric conditions more accurately simulate the shallow subsurface than petroleum reservoirs. For this reason, perhaps the greatest application for analogue sandbox modelling is to improve our understanding of shallow groundwater and environmental flow mechanisms. Two fluid flow simulations were conducted whereby air and kerosene were injected into separate water-saturated unconsolidated sand models. In both experiments, a base 3-D seismic volume was recorded and compared with six later monitor surveys recorded while the injection program was conducted. Normal incidence amplitude and P-wave velocity information were extracted from the TL3-D seismic data to provide visualization of contaminant migration. Reflection amplitudes displayed qualitative areal distribution of fluids when a suitable impedance contrast existed between pore fluids. TL3-D seismic reflection tomography can potentially monitor the change in areal distribution of fluid contaminants over time, indicating flow patterns. However, other research and this current work have not established a quantifiable relationship between either normal reflection amplitudes and attenuation and fluid saturation. Generally, different pore fluids will have unique seismic velocities due to differences in compressibility and density. The predictable
A moving mesh algorithm for 3-D regional groundwater flow with water table and seepage face
NASA Astrophysics Data System (ADS)
Knupp, P.
A numerical algorithm is described for solving the free-surface groundwater flow equations in 3-D large-scale unconfined aquifers with strongly heterogeneous conductivity and surface recharge. The algorithm uses a moving mesh to track the water-table as it evolves according to kinematic and seepage face boundary conditions. Both steady-state and transient algorithms are implemented in the SECO-Flow 3-D code and demonstrated on stratigraphy based on the Delaware Basin of south-eastern New Mexico.
3-D model of a radial flow sub-watt methanol fuel processor
Holladay, J. D.; Wang, Y.
2015-10-01
A 3-D model is presented for a novel sub-watt packed bed reactor. The reactor uses an annular inlet flow combined with a radial flow packed bed reactor. The baseline reactor is compared to a reactor with multiple outlets and a reactor with 3 internal fins. Increasing the outlets from 1 to 4 did improve the flow distribution, but did not increase the performance in the simulation. However, inserting fins allowed a decrease in temperature with same inlet flow of approximately 35K. Or the inlet flow rate could be increased by a factor of 2.8x while maintaining >99% conversion.
Numerical Optimization Strategy for Determining 3D Flow Fields in Microfluidics
NASA Astrophysics Data System (ADS)
Eden, Alex; Sigurdson, Marin; Mezic, Igor; Meinhart, Carl
2015-11-01
We present a hybrid experimental-numerical method for generating 3D flow fields from 2D PIV experimental data. An optimization algorithm is applied to a theory-based simulation of an alternating current electrothermal (ACET) micromixer in conjunction with 2D PIV data to generate an improved representation of 3D steady state flow conditions. These results can be used to investigate mixing phenomena. Experimental conditions were simulated using COMSOL Multiphysics to solve the temperature and velocity fields, as well as the quasi-static electric fields. The governing equations were based on a theoretical model for ac electrothermal flows. A Nelder-Mead optimization algorithm was used to achieve a better fit by minimizing the error between 2D PIV experimental velocity data and numerical simulation results at the measurement plane. By applying this hybrid method, the normalized RMS velocity error between the simulation and experimental results was reduced by more than an order of magnitude. The optimization algorithm altered 3D fluid circulation patterns considerably, providing a more accurate representation of the 3D experimental flow field. This method can be generalized to a wide variety of flow problems. This research was supported by the Institute for Collaborative Biotechnologies through grant W911NF-09-0001 from the U.S. Army Research Office.
3D positional tracking of ellipsoidal particles in a microtube flow using holographic microscopy
NASA Astrophysics Data System (ADS)
Byeon, Hyeok Jun; Seo, Kyung Won; Lee, Sang Joon
2014-11-01
Understanding of micro-scale flow phenomena is getting large attention under advances in micro-scale measurement technologies. Especially, the dynamics of particles suspended in a fluid is essential in both scientific and industrial fields. Moreover, most particles handled in research and industrial fields have non-spherical shapes rather than a simple spherical shape. Under various flow conditions, these non-spherical particles exhibit unique dynamic behaviors. To analyze these dynamic behaviors in a fluid flow, 3D positional information of the particles should be measured accurately. In this study, digital holographic microscopy (DHM) is employed to measure the 3D positional information of non-spherical particles, which are fabricated by stretching spherical polystyrene particles. 3D motions of those particles are obtained by interpreting the holograms captured from particles. Ellipsoidal particles with known size and shape are observed to verify the performance of the DHM technique. In addition, 3D positions of particles in a microtube flow are traced. This DHM technique exhibits promising potential in the analysis of dynamic behaviors of non-spherical particles suspended in micro-scale fluid flows.
One-layer microfluidic device for hydrodynamic 3D self-flow-focusing operating in low flow speed
NASA Astrophysics Data System (ADS)
Daghighi, Yasaman; Gnyawali, Vaskar; Strohm, Eric M.; Tsai, Scott S. H.; Kolios, Michael C.
2016-03-01
Hydrodynamic 3D flow-focusing techniques in microfluidics are categorized as (a) sheathless techniques which require high flow rates and long channels, resulting in high operating cost and high flow rates which are inappropriate for applications with flow rate limitations, and (b) sheath-flow based techniques which usually require excessive sheath flow rate to achieve hydrodynamic 3D flow-focusing. Many devices based on these principles use complicated fabrication methods to create multi-layer microchannels. We have developed a sheath-flow based microfluidic device that is capable of hydrodynamic 3D self-flow-focusing. In this device the main flow (black ink) in a low speed, and a sheath flow, enter through two inlets and enter a 180 degree curved channel (300 × 300 μm cross-section). Main flow migrates outwards into the sheath-flow due to centrifugal effects and consequently, vertical focusing is achieved at the end of the curved channel. Then, two other sheath flows horizontally confine the main flow to achieve horizontal focusing. Thus, the core flow is three-dimensionally focused at the center of the channel at the downstream. Using centrifugal force for 3D flow-focusing in a single-layer fabricated microchannel has been previously investigated by few groups. However, their demonstrated designs required high flow speed (>1 m/s) which is not suitable for many applications that live biomedical specie are involved. Here, we introduce a new design which is operational in low flow speed (<0.05 m/s) and is suitable for applications involving live cells. This microfluidic device can be used in detecting, counting and isolating cells in many biomedical applications.
3D separation over a wall-mounted hemisphere in steady and pulsatile flow
NASA Astrophysics Data System (ADS)
Carr, Ian A.; Plesniak, Michael W.
2014-11-01
Flow separation over a surface-mounted hemispheriod is prevalent in countless applications, both under steady (constant freestream velocity) and unsteady flow over the protuberance. Previous studies of 3D separation have been limited to steady inflow conditions. In biological and geophysical flows, pulsatile flow conditions are much more commonly observed, yet such conditions have not been well studied. Primarily motivated by previous studies of the flow observed in various human vocal fold pathologies, such as polyps, our research aims to fill the knowledge gap in unsteady 3D flow separation. This is achieved by characterizing surface pressure fields and velocity fields, focused primarily on the vortical flow structures and dynamics that occur around a hemispheroid protuberance under pulsatile flow conditions. Surface static pressure and two-dimensional, instantaneous and phase-averaged, particle image velocimetry data in steady and pulsatile flow are presented and compared. Coherent vortical flow structures have been identified using the λci vortex identification criterion. This analysis has revealed a novel set of flow structures dependent on the pulsatile flow forcing function. This material is based in part upon work supported by the National Science Foundation under Grant Number CBET-1236351.
3D Simulation of Velocity Profile of Turbulent Flow in Open Channel with Complex Geometry
NASA Astrophysics Data System (ADS)
Kamel, Benoumessad; Ilhem, Kriba; Ali, Fourar; Abdelbaki, Djebaili
Simulation of open channel flow or river flow presents unique challenge to numerical simulators, which is widely used in the applications of computational fluid dynamics. The prediction is extremely difficult because the flow in open channel is usually transient and turbulent, the geometry is irregular and curved, and the free-surface elevation is varying with time. The results from a 3D non-linear k- ɛ turbulence model are presented to investigate the flow structure, the velocity distribution and mass transport process in a meandering compound open channel and a straight open channel. The 3D numerical model for calculating flow is set up in cylinder coordinates in order to calculate the complex boundary channel. The finite volume method is used to disperse the governing equations and the SIMPLE algorithm is applied to acquire the coupling of velocity and pressure. The non-linear k- ɛ turbulent model has good useful value because of taking into account the anisotropy and not increasing the computational time. The main contributions of this study are developing a numerical method that can be applied to predict the flow in river bends with various bend curvatures and different width-depth ratios. This study demonstrates that the 3D non-linear k- ɛ turbulence model can be used for analyzing flow structures, the velocity distribution and pollutant transport in the complex boundary open channel, this model is applicable for real river and wetland problem.
Commercial turbofan engine exhaust nozzle flow analyses using PAB3D
NASA Technical Reports Server (NTRS)
Abdol-Hamid, Khaled S.; Uenishi, K.; Carlson, John R.; Keith, B. D.
1992-01-01
Recent developments of a three-dimensional (PAB3D) code have paved the way for a computational investigation of complex aircraft aerodynamic components. The PAB3D code was developed for solving the simplified Reynolds Averaged Navier-Stokes equations in a three-dimensional multiblock/multizone structured mesh domain. The present analysis was applied to commercial turbofan exhaust flow systems. Solution sensitivity to grid density is presented. Laminar flow solutions were developed for all grids and two-equation k-epsilon solutions were developed for selected grids. Static pressure distributions, mass flow and thrust quantities were calculated for on-design engine operating conditions. Good agreement between predicted surface static pressures and experimental data was observed at different locations. Mass flow was predicted within 0.2 percent of experimental data. Thrust forces were typically within 0.4 percent of experimental data.
Finite volume and finite element methods applied to 3D laminar and turbulent channel flows
Louda, Petr; Příhoda, Jaromír; Sváček, Petr; Kozel, Karel
2014-12-10
The work deals with numerical simulations of incompressible flow in channels with rectangular cross section. The rectangular cross section itself leads to development of various secondary flow patterns, where accuracy of simulation is influenced by numerical viscosity of the scheme and by turbulence modeling. In this work some developments of stabilized finite element method are presented. Its results are compared with those of an implicit finite volume method also described, in laminar and turbulent flows. It is shown that numerical viscosity can cause errors of same magnitude as different turbulence models. The finite volume method is also applied to 3D turbulent flow around backward facing step and good agreement with 3D experimental results is obtained.
Numerical simulation of a combined oxidation ditch flow using 3D k-epsilon turbulence model.
Luo, Lin; Li, Wei-min; Deng, Yong-sen; Wang, Tao
2005-01-01
The standard three dimensional(3D) k-epsilon turbulence model was applied to simulate the flow field of a small scale combined oxidation ditch. The moving mesh approach was used to model the rotor of the ditch. Comparison of the computed and the measured data is acceptable. A vertical reverse flow zone in the ditch was found, and it played a very important role in the ditch flow behavior. The flow pattern in the ditch is discussed in detail, and approaches are suggested to improve the hydrodynamic performance in the ditch. PMID:16313008
NASA Astrophysics Data System (ADS)
Timmerman, B. H.; Watt, D. W.; Bryanston-Cross, P. J.
1999-02-01
Using holographic interferometry the three-dimensional structure of unsteady and large-scale motions within subsonic and transonic turbulent jet flows has been studied. The instantaneous 3D flow structure is obtained by tomographic reconstruction techniques from quantitative phase maps recorded using a rapid-switching, double reference beam, double pulse laser system. The reconstruction of the jets studied here reveal a three-dimensional nature of the flow. In particular an increasing complexity can be seen in the turbulence as the flow progresses from the jet nozzle. Furthermore, a coherent three-dimensional, possibly rotating, structure can be seen to exist within these jets. The type of flow features illustrated here are not just of fundamental importance for understanding the behavior of free jet flows, but are also common to a number of industrial applications, ranging from the combustion flow within an IC engine to the transonic flow through the stages of a gas turbine.
Modeling of blood flow in arterial trees.
Anor, Tomer; Grinberg, Leopold; Baek, Hyoungsu; Madsen, Joseph R; Jayaraman, Mahesh V; Karniadakis, George E
2010-01-01
Advances in computational methods and medical imaging techniques have enabled accurate simulations of subject-specific blood flows at the level of individual blood cell and in complex arterial networks. While in the past, we were limited to simulations with one arterial bifurcation, the current state-of-the-art is simulations of arterial networks consisting of hundreds of arteries. In this paper, we review the advances in methods for vascular flow simulations in large arterial trees. We discuss alternative approaches and validity of various assumptions often made to simplify the modeling. To highlight the similarities and discrepancies of data computed with different models, computationally intensive three-dimensional (3D) and inexpensive one-dimensional (1D) flow simulations in very large arterial networks are employed. Finally, we discuss the possibilities, challenges, and limitations of the computational methods for predicting outcomes of therapeutic interventions for individual patients. PMID:20836052
ERIC Educational Resources Information Center
Clifford, Philip S.
2011-01-01
Organ blood flow is determined by perfusion pressure and vasomotor tone in the resistance vessels of the organ. Local factors that regulate vasomotor tone include myogenic and metabolic autoregulation, flow-mediated and conducted responses, and vasoactive substances released from red blood cells. The relative importance of each of these factors…
An improved parallel SPH approach to solve 3D transient generalized Newtonian free surface flows
NASA Astrophysics Data System (ADS)
Ren, Jinlian; Jiang, Tao; Lu, Weigang; Li, Gang
2016-08-01
In this paper, a corrected parallel smoothed particle hydrodynamics (C-SPH) method is proposed to simulate the 3D generalized Newtonian free surface flows with low Reynolds number, especially the 3D viscous jets buckling problems are investigated. The proposed C-SPH method is achieved by coupling an improved SPH method based on the incompressible condition with the traditional SPH (TSPH), that is, the improved SPH with diffusive term and first-order Kernel gradient correction scheme is used in the interior of the fluid domain, and the TSPH is used near the free surface. Thus the C-SPH method possesses the advantages of two methods. Meanwhile, an effective and convenient boundary treatment is presented to deal with 3D multiple-boundary problem, and the MPI parallelization technique with a dynamic cells neighbor particle searching method is considered to improve the computational efficiency. The validity and the merits of the C-SPH are first verified by solving several benchmarks and compared with other results. Then the viscous jet folding/coiling based on the Cross model is simulated by the C-SPH method and compared with other experimental or numerical results. Specially, the influences of macroscopic parameters on the flow are discussed. All the numerical results agree well with available data, and show that the C-SPH method has higher accuracy and better stability for solving 3D moving free surface flows over other particle methods.
Intuitive Visualization of Transient Flow: Towards a Full 3D Tool
NASA Astrophysics Data System (ADS)
Michel, Isabel; Schröder, Simon; Seidel, Torsten; König, Christoph
2015-04-01
Visualization of geoscientific data is a challenging task especially when targeting a non-professional audience. In particular, the graphical presentation of transient vector data can be a significant problem. With STRING Fraunhofer ITWM (Kaiserslautern, Germany) in collaboration with delta h Ingenieurgesellschaft mbH (Witten, Germany) developed a commercial software for intuitive 2D visualization of 3D flow problems. Through the intuitive character of the visualization experts can more easily transport their findings to non-professional audiences. In STRING pathlets moving with the flow provide an intuition of velocity and direction of both steady-state and transient flow fields. The visualization concept is based on the Lagrangian view of the flow which means that the pathlets' movement is along the direction given by pathlines. In order to capture every detail of the flow an advanced method for intelligent, time-dependent seeding of the pathlets is implemented based on ideas of the Finite Pointset Method (FPM) originally conceived at and continuously developed by Fraunhofer ITWM. Furthermore, by the same method pathlets are removed during the visualization to avoid visual cluttering. Additional scalar flow attributes, for example concentration or potential, can either be mapped directly to the pathlets or displayed in the background of the pathlets on the 2D visualization plane. The extensive capabilities of STRING are demonstrated with the help of different applications in groundwater modeling. We will discuss the strengths and current restrictions of STRING which have surfaced during daily use of the software, for example by delta h. Although the software focusses on the graphical presentation of flow data for non-professional audiences its intuitive visualization has also proven useful to experts when investigating details of flow fields. Due to the popular reception of STRING and its limitation to 2D, the need arises for the extension to a full 3D tool
A fast and accurate method to predict 2D and 3D aerodynamic boundary layer flows
NASA Astrophysics Data System (ADS)
Bijleveld, H. A.; Veldman, A. E. P.
2014-12-01
A quasi-simultaneous interaction method is applied to predict 2D and 3D aerodynamic flows. This method is suitable for offshore wind turbine design software as it is a very accurate and computationally reasonably cheap method. This study shows the results for a NACA 0012 airfoil. The two applied solvers converge to the experimental values when the grid is refined. We also show that in separation the eigenvalues remain positive thus avoiding the Goldstein singularity at separation. In 3D we show a flow over a dent in which separation occurs. A rotating flat plat is used to show the applicability of the method for rotating flows. The shown capabilities of the method indicate that the quasi-simultaneous interaction method is suitable for design methods for offshore wind turbine blades.
3D flow past transonic turbine cascade SE 1050 — Experiment and numerical simulations
NASA Astrophysics Data System (ADS)
Šimurda, D.; Fürst, J.; Luxa, M.
2013-08-01
This paper is concerned with experimental and numerical research on 3D flow past prismatic turbine cascade SE1050 (known in QNET network as open test case SE1050). The primary goal was to assess the influence of the inlet velocity profile on the flow structures in the interblade channel and on the flow field parameters at the cascade exit and to compare these findings to results of numerical simulations. Investigations of 3D flow past the cascade with non-uniform inlet velocity profile were carried out both experimentally and numerically at subsonic ( M 2is = 0.8) and at transonic ( M 2is = 1.2) regime at design angle of incidence. Experimental data was obtained using a traversing device with a five-hole conical probe. Numerically, the 3D flow was simulated by open source code OpenFOAM and in-house code. Analyses of experimental data and CFD simulations have revealed the development of distinctive vortex structures resulting from non-uniform inlet velocity profile. Origin of these structures results in increased loss of kinetic energy and spanwise shift of kinetic energy loss coefficient distribution. Differences found between the subsonic and the transonic case confirm earlier findings available in the literature. Results of CFD and experiments agree reasonably well.
Subramanian, K R; Thubrikar, M J; Fowler, B; Mostafavi, M T; Funk, M W
2000-01-01
We present a technique that accurately reconstructs complex three dimensional blood vessel geometry from 2D intravascular ultrasound (IVUS) images. Biplane x-ray fluoroscopy is used to image the ultrasound catheter tip at a few key points along its path as the catheter is pulled through the blood vessel. An interpolating spline describes the continuous catheter path. The IVUS images are located orthogonal to the path, resulting in a non-uniform structured scalar volume of echo densities. Isocontour surfaces are used to view the vessel geometry, while transparency and clipping enable interactive exploration of interior structures. The two geometries studied are a bovine artery vascular graft having U-shape and a constriction, and a canine carotid artery having multiple branches and a constriction. Accuracy of the reconstructions is established by comparing the reconstructions to (1) silicone moulds of the vessel interior, (2) biplane x-ray images, and (3) the original echo images. Excellent shape and geometry correspondence was observed in both geometries. Quantitative measurements made at key locations of the 3D reconstructions also were in good agreement with those made in silicone moulds. The proposed technique is easily adoptable in clinical practice, since it uses x-rays with minimal exposure and existing IVUS technology. PMID:11105284
Numerical modelling of gravel unconstrained flow experiments with the DAN3D and RASH3D codes
NASA Astrophysics Data System (ADS)
Sauthier, Claire; Pirulli, Marina; Pisani, Gabriele; Scavia, Claudio; Labiouse, Vincent
2015-12-01
Landslide continuum dynamic models have improved considerably in the last years, but a consensus on the best method of calibrating the input resistance parameter values for predictive analyses has not yet emerged. In the present paper, numerical simulations of a series of laboratory experiments performed at the Laboratory for Rock Mechanics of the EPF Lausanne were undertaken with the RASH3D and DAN3D numerical codes. They aimed at analysing the possibility to use calibrated ranges of parameters (1) in a code different from that they were obtained from and (2) to simulate potential-events made of a material with the same characteristics as back-analysed past-events, but involving a different volume and propagation path. For this purpose, one of the four benchmark laboratory tests was used as past-event to calibrate the dynamic basal friction angle assuming a Coulomb-type behaviour of the sliding mass, and this back-analysed value was then used to simulate the three other experiments, assumed as potential-events. The computational findings show good correspondence with experimental results in terms of characteristics of the final deposits (i.e., runout, length and width). Furthermore, the obtained best fit values of the dynamic basal friction angle for the two codes turn out to be close to each other and within the range of values measured with pseudo-dynamic tilting tests.
Delft3D-FLOW on PRACE infrastructures for real life hydrodynamic applications.
NASA Astrophysics Data System (ADS)
Donners, John; Genseberger, Menno; Jagers, Bert; de Goede, Erik; Mourits, Adri
2013-04-01
PRACE, the Partnership for Advanced Computing in Europe, offers access to the largest high-performance computing systems in Europe. PRACE invites and helps industry to increase their innovative potential through the use of the PRACE infrastructure. This poster describes different efforts to assist Deltares with porting the open-source simulation software Delft3D-FLOW to PRACE infrastructures. Analysis of the performance on these infrastructures has been done for real life flow applications. Delft3D-FLOW is a 2D and 3D shallow water solver which calculates non-steady flow and transport phenomena resulting from tidal and meteorological forcing on a curvilinear, boundary fitted grid in Cartesian or spherical coordinates. It also includes a module which sediment transport (both suspended and bed total load) and morphological changes for an arbitrary number of cohesive and non-cohesive fractions. As Delft3D-FLOW has been developed over several decades, with a variety of functionality and over 350k lines of source code, porting to PRACE infrastructures needs some effort. At the moment Delft3D-FLOW uses MPI with domain decomposition in one direction as its parallellisation approach. Because it is hard to identify scaling issues if one immediately starts with a complex case with many features enabled, different cases with increasing complexity have been used to investigate scaling of this parallellisation approach on several PRACE platforms. As a base reference case we started with a schematic high-resolution 2D hydrodynamic model of the river Waal that turned out to be surprisingly well-suited to the highly-parallel PRACE machines. Although Delft3D-FLOW employs a sophisticated build system, several modifications were required to port it to most PRACE systems due to the use of specific, highly-tuned compilers and MPI-libraries. After this we moved to a 3D hydrodynamic model of Rotterdam harbour that includes sections of the rivers Rhine and Meuse and a part of the North
NASA Technical Reports Server (NTRS)
Fleming, J. L.; Simpson, R. L.
1997-01-01
Laser Doppler velocimetry (LDV) measurements and hydrogen bubble flow visualization techniques were used to examine the near-wall flow structure of 2D and 3D turbulent boundary layers (TBLs) over a range of low Reynolds numbers. The goals of this research were (1) an increased understanding of the flow physics in the near wall region of turbulent boundary layers,(2) to observe and quantify differences between 2D and 3D TBL flow structures, and (3) to document Reynolds number effects for 3D TBLs. The LDV data have provided results detailing the turbulence structure of the 2D and 3D TBLs. These results include mean Reynolds stress distributions, flow skewing results, and U and V spectra. Effects of Reynolds number for the 3D flow were also examined. Comparison to results with the same 3D flow geometry but at a significantly higher Reynolds number provided unique insight into the structure of 3D TBLs. While the 3D mean and fluctuating velocities were found to be highly dependent on Reynolds number, a previously defined shear stress parameter was discovered to be invariant with Reynolds number. The hydrogen bubble technique was used as a flow visualization tool to examine the near-wall flow structure of 2D and 3D TBLs. Both the quantitative and qualitative results displayed larger turbulent fluctuations with more highly concentrated vorticity regions for the 2D flow.
3-D High-Lift Flow-Physics Experiment - Transition Measurements
NASA Technical Reports Server (NTRS)
McGinley, Catherine B.; Jenkins, Luther N.; Watson, Ralph D.; Bertelrud, Arild
2005-01-01
An analysis of the flow state on a trapezoidal wing model from the NASA 3-D High Lift Flow Physics Experiment is presented. The objective of the experiment was to characterize the flow over a non-proprietary semi-span three-element high-lift configuration to aid in assessing the state of the art in the computation of three-dimensional high-lift flows. Surface pressures and hot-film sensors are used to determine the flow conditions on the slat, main, and flap. The locations of the attachments lines and the values of the attachment line Reynolds number are estimated based on the model surface pressures. Data from the hot-films are used to determine if the flow is laminar, transitional, or turbulent by examining the hot-film time histories, statistics, and frequency spectra.
Adaptive volume rendering of cardiac 3D ultrasound images: utilizing blood pool statistics
NASA Astrophysics Data System (ADS)
Åsen, Jon Petter; Steen, Erik; Kiss, Gabriel; Thorstensen, Anders; Rabben, Stein Inge
2012-03-01
In this paper we introduce and investigate an adaptive direct volume rendering (DVR) method for real-time visualization of cardiac 3D ultrasound. DVR is commonly used in cardiac ultrasound to visualize interfaces between tissue and blood. However, this is particularly challenging with ultrasound images due to variability of the signal within tissue as well as variability of noise signal within the blood pool. Standard DVR involves a global mapping of sample values to opacity by an opacity transfer function (OTF). While a global OTF may represent the interface correctly in one part of the image, it may result in tissue dropouts, or even artificial interfaces within the blood pool in other parts of the image. In order to increase correctness of the rendered image, the presented method utilizes blood pool statistics to do regional adjustments of the OTF. The regional adaptive OTF was compared with a global OTF in a dataset of apical recordings from 18 subjects. For each recording, three renderings from standard views (apical 4-chamber (A4C), inverted A4C (IA4C) and mitral valve (MV)) were generated for both methods, and each rendering was tuned to the best visual appearance by a physician echocardiographer. For each rendering we measured the mean absolute error (MAE) between the rendering depth buffer and a validated left ventricular segmentation. The difference d in MAE between the global and regional method was calculated and t-test results are reported with significant improvements for the regional adaptive method (dA4C = 1.5 +/- 0.3 mm, dIA4C = 2.5 +/- 0.4 mm, dMV = 1.7 +/- 0.2 mm, d.f. = 17, all p < 0.001). This improvement by the regional adaptive method was confirmed through qualitative visual assessment by an experienced physician echocardiographer who concluded that the regional adaptive method produced rendered images with fewer tissue dropouts and less spurious structures inside the blood pool in the vast majority of the renderings. The algorithm has been
Simulation of 3-D Nonequilibrium Seeded Air Flow in the NASA-Ames MHD Channel
NASA Technical Reports Server (NTRS)
Gupta, Sumeet; Tannehill, John C.; Mehta, Unmeel B.
2004-01-01
The 3-D nonequilibrium seeded air flow in the NASA-Ames experimental MHD channel has been numerically simulated. The channel contains a nozzle section, a center section, and an accelerator section where magnetic and electric fields can be imposed on the flow. In recent tests, velocity increases of up to 40% have been achieved in the accelerator section. The flow in the channel is numerically computed us ing a 3-D parabolized Navier-Stokes (PNS) algorithm that has been developed to efficiently compute MHD flows in the low magnetic Reynolds number regime: The MHD effects are modeled by introducing source terms into the PNS equations which can then be solved in a very efficient manner. The algorithm has been extended in the present study to account for nonequilibrium seeded air flows. The electrical conductivity of the flow is determined using the program of Park. The new algorithm has been used to compute two test cases that match the experimental conditions. In both cases, magnetic and electric fields are applied to the seeded flow. The computed results are in good agreement with the experimental data.
Evaluation of the 3-D channeling flow in a fractured type of oil/gas reservoir
NASA Astrophysics Data System (ADS)
Ishibashi, T.; Watanabe, N.; Tsuchiya, N.; Tamagawa, T.
2013-12-01
An understanding of the flow and transport characteristics through rock fracture networks is of critical importance in many engineering and scientific applications. These include effective recovery of targeted fluid such as oil/gas, geothermal, or potable waters, and isolation of hazardous materials. Here, the formation of preferential flow path (i.e. channeling flow) is one of the most significant characteristics in considering fluid flow through rock fracture networks; however, the impact of channeling flow remains poorly understood. In order to deepen our understanding of channeling flow, the authors have developed a novel discrete fracture network (DFN) model simulator, GeoFlow. Different from the conventional DFN model simulators, we can characterize each fracture not by a single aperture value but by a heterogeneous aperture distribution in GeoFlow [Ishibashi et al., 2012]. As a result, the formation of 3-D preferential flow paths within fracture network can be considered by using this simulator. Therefore, we would challenge to construct the precise fracture networks whose fractures have heterogeneous aperture distributions in field scale, and to analyze fluid flows through the fracture networks by GeoFlow. In the present study, the Yufutsu oil/gas field in Hokkaido, Japan is selected as the subject area for study. This field is known as the fractured type of reservoir, and reliable DFN models can be constructed for this field based on the 3-D seismic data, well logging, in-situ stress measurement, and acoustic emission data [Tamagawa et al., 2012]. Based on these DFN models, new DFN models for 1,080 (East-West) × 1,080 (North-South) × 1,080 (Depth) m^3, where fractures are represented by squares of 44-346 m on a side, are re-constructed. In these new models, scale-dependent aperture distributions are considered for all fractures constructing the fracture networks. Note that the multi-scale modeling of fracture flow has been developed by the authors
Quasi 3D modeling of water flow and solute transport in vadose zone and groundwater
NASA Astrophysics Data System (ADS)
Yakirevich, A.; Kuznetsov, M.; Weisbrod, N.; Pachepsky, Y. A.
2013-12-01
The complexity of subsurface flow systems calls for a variety of concepts leading to the multiplicity of simplified flow models. One commonly used simplification is based on the assumption that lateral flow and transport in unsaturated zone is insignificant unless the capillary fringe is involved. In such cases the flow and transport in the unsaturated zone above groundwater level can be simulated as a 1D phenomenon, whereas through groundwater they are viewed as 2D or 3D phenomena. A new approach for a numerical scheme for 3D variably saturated flow and transport is presented. A Quasi-3D approach allows representing flow in the 'vadose zone - aquifer' system by a series of 1D Richards' equations solved in variably-saturated zone and by 3D-saturated flow equation in groundwater (modified MODFLOW code). The 1D and 3D equations are coupled at the phreatic surface in a way that aquifer replenishment is calculated using the Richards' equation, and solving for the moving water table does not require definition of the specific yield parameter. The 3D advection-dispersion equation is solved in the entire domain by the MT3D code. Using implicit finite differences approximation to couple processes in the vadose zone and groundwater provides mass conservation and increase of computational efficiency. The above model was applied to simulate the impact of irrigation on groundwater salinity in the Alto Piura aquifer (Northern Peru). Studies on changing groundwater quality in arid and semi-arid lands show that irrigation return flow is one of the major factors contributing to aquifer salinization. Existing mathematical models do not account explicitly for the solute recycling during irrigation on a daily scale. Recycling occurs throughout the unsaturated and saturated zones, as function of the solute mass extracted from pumping wells. Salt concentration in irrigation water is calculated at each time step as a function of concentration of both surface water and groundwater
3-D segmentation of retinal blood vessels in spectral-domain OCT volumes of the optic nerve head
NASA Astrophysics Data System (ADS)
Lee, Kyungmoo; Abràmoff, Michael D.; Niemeijer, Meindert; Garvin, Mona K.; Sonka, Milan
2010-03-01
Segmentation of retinal blood vessels can provide important information for detecting and tracking retinal vascular diseases including diabetic retinopathy, arterial hypertension, arteriosclerosis and retinopathy of prematurity (ROP). Many studies on 2-D segmentation of retinal blood vessels from a variety of medical images have been performed. However, 3-D segmentation of retinal blood vessels from spectral-domain optical coherence tomography (OCT) volumes, which is capable of providing geometrically accurate vessel models, to the best of our knowledge, has not been previously studied. The purpose of this study is to develop and evaluate a method that can automatically detect 3-D retinal blood vessels from spectral-domain OCT scans centered on the optic nerve head (ONH). The proposed method utilized a fast multiscale 3-D graph search to segment retinal surfaces as well as a triangular mesh-based 3-D graph search to detect retinal blood vessels. An experiment on 30 ONH-centered OCT scans (15 right eye scans and 15 left eye scans) from 15 subjects was performed, and the mean unsigned error in 3-D of the computer segmentations compared with the independent standard obtained from a retinal specialist was 3.4 +/- 2.5 voxels (0.10 +/- 0.07 mm).
3-D trajectory model for MDT using micro-spheres implanted within large blood vessels
NASA Astrophysics Data System (ADS)
Choomphon-anomakhun, Natthaphon; Natenapit, Mayuree
2016-09-01
Implant assisted magnetic drug targeting (IA-MDT) using ferromagnetic spherical targets implanted within large blood vessels and subjected to a uniform externally applied magnetic field (H0) has been investigated and reported for the first time. The capture areas (As) of magnetic drug carrier particles (MDCPs) were determined from the analysis of particle trajectories simulated from equations of motion. Then, the effects of various parameters, such as types of ferromagnetic materials in the targets and MDCPs, blood flow rates, mass fraction of the ferromagnetic material in the MDCPs, average radii of MDCPs (Rp) and the strength of H0 on the As were obtained. Furthermore, the effects of saturation magnetization of the ferromagnetic materials in the MDCPs and within the targets on the As were analyzed. After this, the suitable strengths of H0 and Rp for IA-MDT designs were reported. Dimensionless As, ranging from 2 to 7, was obtained with Rp ranging from 500 to 2500 nm, μ0H0 less than 0.8 T and a blood flow rate of 0.1 m s-1. The target-MDCP materials considered are iron-iron, iron-magnetite and SS409-magnetite, respectively.
NASA Astrophysics Data System (ADS)
Kincaid, C. R.; MacDougall, J. G.; Druken, K. A.; Fischer, K. M.
2010-12-01
Understanding patterns in plate scale mantle flow in subduction zones is key to models of thermal structure, dehydration reactions, volatile distributions and magma generation and transport in convergent margins. Different patterns of flow in the mantle wedge can generate distinct signatures in seismological observables. Observed shear wave fast polarization directions in several subduction zones are inconsistent with predictions of simple 2-D wedge corner flow. Geochemical signatures in a number of subduction zones also indicate 3-D flow and entrainment patterns in the wedge. We report on a series of laboratory experiments on subduction driven flow to characterize spatial and temporal variability in 3-D patterns in flow and shear-induced finite strain. Cases focus on how rollback subduction, along-strike dip changes in subducting plates and evolving gaps or tears in subduction zones control temporal-spatial patterns in 3-D wedge flow. Models utilize a glucose working fluid with a temperature dependent viscosity to represent the upper 2000 km of the mantle. Subducting lithosphere is modeled with two rubber-reinforced continuous belts. Belts pass around trench and upper/lower mantle rollers. The deeper rollers can move laterally to allow for time varying dip angle. Each belt has independent speed control and dip adjustment, allowing for along-strike changes in convergence rate and the evolution of slab gaps. Rollback is modeled using a translation system to produce either uniform and asymmetric lateral trench motion. Neutral density finite strain markers are distributed throughout the fluid and used as proxies for tracking the evolution of anisotropy through space and time in the evolving flow fields. Particle image velocimetry methods are also used to track time varying 3-D velocity fields for directly calculating anisotropy patterns. Results show that complex plate motions (rollback, steepening) and morphologies (gaps) in convergent margins produce flows with
3D Particle image velocimetry test of inner flow in a double blade pump impeller
NASA Astrophysics Data System (ADS)
Liu, Houlin; Wang, Kai; Yuan, Shouqi; Tan, Minggao; Wang, Yong; Ru, Weimin
2012-05-01
The double blade pump is widely used in sewage treatment industry, however, the research on the internal flow characteristics of the double blade pump with particle image velocimetry (PIV) technology is very little at present. To reveal inner flow characteristics in double blade pump impeller under off-design and design conditions, inner flows in a double blade pump impeller, whose specific speed is 111, are measured under the five off-design conditions and design condition by using 3D PIV test technology. In order to ensure the accuracy of the 3D PIV test, the external trigger synchronization system which makes use of fiber optic and equivalent calibration method are applied. The 3D PIV relative velocity synthesis procedure is compiled by using Visual C++ 2005. Then absolute velocity distribution and relative velocity distribution in the double blade pump impeller are obtained. Test results show that vortex exists in each condition, but the location, size and velocity of vortex core are different. Average absolute velocity value of impeller outlet increases at first, then decreases, and then increases again with increase of flow rate. Again average relative velocity values under 0.4, 0.8, and 1.2 design condition are higher than that under 1.0 design condition, while under 0.6 and 1.4 design condition it is lower. Under low flow rate conditions, radial vectors of absolute velocities at impeller outlet and blade inlet near the pump shaft decrease with increase of flow rate, while that of relative velocities at the suction side near the pump shaft decreases. Radial vectors of absolute velocities and relative velocities change slightly under the two large flow rate conditions. The research results can be applied to instruct the hydraulic optimization design of double blade pumps.
Numerical simulation of internal and external inviscid and viscous 3-D flow fields
NASA Astrophysics Data System (ADS)
Leicher, Stefan
1986-11-01
A numerical method for solving the 3-D Euler equations in geometrical complex domains was developed. The approach divides the computational space into multiple blocks whose structure follows the natural lines of the conficuration. A systematic, multi-block grid generation scheme is used to produce the grid. The flow solutions are obtained by solving the Euler equations by a finite volume discretization and a Runge-Kutta time stepping scheme. The main advantage of this method is the applicability to complex geometries, for example complete aircraft configurations including wing, fuselage, canard and tail. The coupling with a 3-D boundary layer method allows to account for viscous effects. Another application for the method was the simulation of flows in the presence of a propeller.
Simulation of a 3D unsteady flow in an axial turbine stage
NASA Astrophysics Data System (ADS)
Straka, Petr
2012-04-01
The contribution deals with a numerical simulation of an unsteady flow in an axial turbine stage. The solution is performed using an in-house numerical code developed in the Aeronautical and Test Institute, Plc. in Prague. The numerical code is based on a finite volume discretization of governing equations (Favre averaged Navier-Stokes equations) and a two-equations turbulence model. The temporal integration is based on the implicit second-order backward Euler formula, which is realized through the iteration process in dual time. The proposed numerical method is used for solution of the 3D, unsteady, viscous turbulent flow of a perfect gas in the axial turbine stage. The flow path consists of an input nozzle, stator blade-wheel, rotor blade-wheel, a shroud-seal gap and a diffuser. Attention is paid to the influence of a secondary flow structures, such as generated vortices and flow in shroud-seal gap.
Numerical analysis of the aeroelastic behaviour for the last turbine stage in 3D transonic flow
NASA Astrophysics Data System (ADS)
Gnesin, Vitaly; Kolodyazhnaya, Lyubov
2004-11-01
An understanding of the physics of the mutual interaction between gas flow and oscillating blades, and the development of predictive capabilities is essential for improving overall efficiency, durability and reliability. In this study presented the algorithm proposed involving the coupled solution of 3D unsteady flow through a turbine stage and dynamic problem for rotor blades motion by action of aerodynamic forces without separating outer and inner flow fluctuations. There has been performed the calculations for the last stage of the steam turbine under design and off-design regimes. It has investigated the mutual influence of both outer flow non-uniformity and blades oscillations. It has shown that amplitude-frequency spectrum of blade oscillations contains the high frequency harmonics, corresponding to rotor moving one stator blade pitch, and low frequency harmonics caused by blade oscillations and flow non-uniformity downstream from the blade row.
Ultrarapid detection of pathogenic bacteria using a 3D immunomagnetic flow assay.
Lee, Wonjae; Kwon, Donghoon; Chung, Boram; Jung, Gyoo Yeol; Au, Anthony; Folch, Albert; Jeon, Sangmin
2014-07-01
We developed a novel 3D immunomagnetic flow assay for the rapid detection of pathogenic bacteria in a large-volume food sample. Antibody-functionalized magnetic nanoparticle clusters (AbMNCs) were magnetically immobilized on the surfaces of a 3D-printed cylindrical microchannel. The injection of a Salmonella-spiked sample solution into the microchannel produced instant binding between the AbMNCs and the Salmonella bacteria due to their efficient collisions. Nearly perfect capture of the AbMNCs and AbMNCs-Salmonella complexes was achieved under a high flow rate by stacking permanent magnets with spacers inside the cylindrical separator to maximize the magnetic force. The concentration of the bacteria in solution was determined using ATP luminescence measurements. The detection limit was better than 10 cfu/mL, and the overall assay time, including the binding, rinsing, and detection steps for a 10 mL sample took less than 3 min. To our knowledge, the 3D immunomagnetic flow assay described here provides the fastest high-sensitivity, high-capacity method for the detection of pathogenic bacteria. PMID:24856003
Ultrarapid Detection of Pathogenic Bacteria Using a 3D Immunomagnetic Flow Assay
Lee, Wonjae; Kwon, Donghoon; Chung, Boram; Jung, Gyoo Yeol; Au, Anthony; Folch, Albert; Jeon, Sangmin
2015-01-01
We developed a novel 3D immunomagnetic flow assay for the rapid detection of pathogenic bacteria in a large-volume food sample. Antibody-functionalized magnetic nanoparticle clusters (AbMNCs) were magnetically immobilized on the surfaces of a 3D-printed cylindrical microchannel. The injection of a Salmonella-spiked sample solution into the microchannel produced instant binding between the AbMNCs and the Salmonella bacteria due to their efficient collisions. Nearly perfect capture of the AbMNCs and AbMNCs-Salmonella complexes was achieved under a high flow rate by stacking permanent magnets with spacers inside the cylindrical separator to maximize the magnetic force. The concentration of the bacteria in solution was determined using ATP luminescence measurements. The detection limit was better than 10 cfu/mL, and the overall assay time, including the binding, rinsing, and detection steps for a 10 mL sample took less than 3 min. To our knowledge, the 3D immunomagnetic flow assay described here provides the fastest high-sensitivity, high-capacity method for the detection of pathogenic bacteria. PMID:24856003
Ababou, R.
1996-12-31
Subsurface flow processes are inherently three-dimensional and heterogeneous over many scales. Taking this into account, for instance assuming random heterogeneity in 3-D space, puts heavy constraints on numerical models. An efficient numerical code has been developed for solving the porous media flow equations, appropriately generalized to account for 3-D, random-like heterogeneity. The code is based on implicit finite differences (or finite volumes), and uses specialized versions of pre-conditioned iterative solvers that take advantage of sparseness. With Diagonally Scaled Conjugate Gradients, in particular, large systems on the order of several million equations, with randomly variable coefficients, have been solved efficiently on Cray-2 and Cray-Y/MP8 machines, in serial mode as well as parallel mode (autotasking). The present work addresses, first, the numerical aspects and computational issues associated with detailed 3-D flow simulations, and secondly, presents a specific application related to the conductivity homogenization problem (identifying a macroscale conduction law, and an equivalent or effective conductivity). Analytical expressions of effective conductivities are compared with empirical values obtained from several large scale simulations conducted for single realizations of random porous media.
Endovascular blood flow measurement system
NASA Astrophysics Data System (ADS)
Khe, A. K.; Cherevko, A. A.; Chupakhin, A. P.; Krivoshapkin, A. L.; Orlov, K. Yu
2016-06-01
In this paper an endovascular measurement system used for intraoperative cerebral blood flow monitoring is described. The system is based on a Volcano ComboMap Pressure and Flow System extended with analogue-to-digital converter and PC laptop. A series of measurements performed in patients with cerebrovascular pathologies allows us to introduce “velocity-pressure” and “flow rate-energy flow rate” diagrams as important characteristics of the blood flow. The measurement system presented here can be used as an additional instrument in neurosurgery for assessment and monitoring of the operation procedure. Clinical data obtained with the system are used for construction of mathematical models and patient-specific simulations. The monitoring of the blood flow parameters during endovascular interventions was approved by the Ethics Committee at the Meshalkin Novosibirsk Research Institute of Circulation Pathology and included in certain surgical protocols for pre-, intra- and postoperative examinations.
3-D Numerical Modeling of MHD Flows in Variable Magnetic Field
NASA Astrophysics Data System (ADS)
Abdullina, K. I.; Bogovalov, S. V.
3-D numerical simulation of the liquid metal flow affected by the electromagnetic field in the magnetohydrodynamic (MHD) devices is performed. Software package ANSYS has been used for the numerical calculations. The non-stationary problem has been solved taking into account the influence of the metal flow on the electromagnetic field and nonlinear magnetic permeability of the ferromagnetic cores. Simplified calculations with constant magnetic permeability of the ferromagnetic cores have been performed as well. Comparison of these calculations shows that the simulation of the MHD pump can be performed in the linear approximation. The pump performance curve has been derived in this approximation.
Nonhydrostatic granular flow over 3-D terrain: New Boussinesq-type gravity waves?
NASA Astrophysics Data System (ADS)
Castro-Orgaz, Oscar; Hutter, Kolumban; Giraldez, Juan V.; Hager, Willi H.
2015-01-01
granular mass flow is a basic step in the prediction and control of natural or man-made disasters related to avalanches on the Earth. Savage and Hutter (1989) pioneered the mathematical modeling of these geophysical flows introducing Saint-Venant-type mass and momentum depth-averaged hydrostatic equations using the continuum mechanics approach. However, Denlinger and Iverson (2004) found that vertical accelerations in granular mass flows are of the same order as the gravity acceleration, requiring the consideration of nonhydrostatic modeling of granular mass flows. Although free surface water flow simulations based on nonhydrostatic depth-averaged models are commonly used since the works of Boussinesq (1872, 1877), they have not yet been applied to the modeling of debris flow. Can granular mass flow be described by Boussinesq-type gravity waves? This is a fundamental question to which an answer is required, given the potential to expand the successful Boussinesq-type water theory to granular flow over 3-D terrain. This issue is explored in this work by generalizing the basic Boussinesq-type theory used in civil and coastal engineering for more than a century to an arbitrary granular mass flow using the continuum mechanics approach. Using simple test cases, it is demonstrated that the above question can be answered in the affirmative way, thereby opening a new framework for the physical and mathematical modeling of granular mass flow in geophysics, whereby the effect of vertical motion is mathematically included without the need of ad hoc assumptions.
ERIC Educational Resources Information Center
Lassen, Niels A.; And Others
1978-01-01
Discusses the use of radioactive isotopes to graphically represent changes in the amount of blood flowing in areas of the human cerebral cortex, reflecting changes in the activity of those areas. Numerous illustrations are included. (Author/MA)
NASA Astrophysics Data System (ADS)
Harlander, U.; Wright, G. B.; Egbers, C.
2012-04-01
In the earth's atmosphere baroclinic instability is responsible for the heat and momentum transport from low to high latitudes. In the fifties, Raymond Hide used a rather simple laboratory experiment to study such vortices in the lab. The experiment is comprised by a cooled inner and heated outer cylinder mounted on a rotating platform, which mimics the heated tropical and cooled polar regions of the earth's atmosphere. The experiment shows rich dynamics that have been studied by varying the radial temperature difference and the rate of annulus revolution. At the Brandenburg University of Technology (BTU) Cottbus the differentially heated rotating annulus is a reference experiment of the DFG priority program 'MetStröm'. The 3D structure of the annulus flow field has been numerically simulated but, to our knowledge, has not been measured in the laboratory. In the present paper we use novel interpolation techniques to reconstruct the 3D annulus flow field from synchronous Particle Image Velocimetry (PIV) and Infrared Thermography (IRT) measurements. The PIV system is used to measure the horizontal velocity components at 40, 60, 80, 100, and 120 mm above the bottom. The uppermost level is thus 15 mm below the fluid's surface. The surface temperature is simultaneously measured by an infrared (IR) camera. The PIV and infrared cameras have been mounted above the annulus and they co-rotate with the annulus. From the PIV observations alone a coherent 3D picture of the flow cannot be constructed since the PIV measurements have been taken at different instants of time. Therefore a corresponding IR image has been recorded for each PIV measurement. These IR images can be used to reconstruct the correct phase of the measured velocity fields. Each IR and PIV image for which t>0 is rotated back to the position at t=0. Then all surface waves have the same phase. In contrast, the PIV velocity fields generally have different phases since they have been taken at different vertical
Numerical Simulation of 3-D Supersonic Viscous Flow in an Experimental MHD Channel
NASA Technical Reports Server (NTRS)
Kato, Hiromasa; Tannehill, John C.; Gupta, Sumeet; Mehta, Unmeel B.
2004-01-01
The 3-D supersonic viscous flow in an experimental MHD channel has been numerically simulated. The experimental MHD channel is currently in operation at NASA Ames Research Center. The channel contains a nozzle section, a center section, and an accelerator section where magnetic and electric fields can be imposed on the flow. In recent tests, velocity increases of up to 40% have been achieved in the accelerator section. The flow in the channel is numerically computed using a new 3-D parabolized Navier-Stokes (PNS) algorithm that has been developed to efficiently compute MHD flows in the low magnetic Reynolds number regime. The MHD effects are modeled by introducing source terms into the PNS equations which can then be solved in a very e5uent manner. To account for upstream (elliptic) effects, the flowfield can be computed using multiple streamwise sweeps with an iterated PNS algorithm. The new algorithm has been used to compute two test cases that match the experimental conditions. In both cases, magnetic and electric fields are applied to the flow. The computed results are in good agreement with the available experimental data.
Flow measurements in a model centrifugal pump by 3-D PIV
NASA Astrophysics Data System (ADS)
Yang, H.; Xu, H. R.; Liu, C.
2012-11-01
PIV (Particle Image Velocimetry), as an non-intrusive flow measurements technology, is widely used to investigate the flow fields in many areas. 3-D (three Dimensional) PIV has seldom been used to measure flow field in rotational impeller of centrifugal pump due to the difficulty of calibration in samll space. In this article, a specially manufactured water tank was used to perform the calibration for 3-D PIV measurement. The instantaneous absolute velocity in one impeller passage was obtained by merging of three sub zones and the relative velocity was acquired by velocity decomposition. The result shows that, when the pump runs at the condition of design flow rate, the radial component velocity Wr appears a concave distribution except the condition of R=45 mm. With the increase of radius, the circumference location of the minimum radial component velocity Wr moves from the pressure side to the suction side. At the same time, the tangential component velocity Wθ on the suction side decreases gradually with the increase of radius, while the component on the pressure side increases gradually. The secondary flow in different radius section has also been shown. At last, the error of PIV measurements was analyzed, which shows that the test results are accurate and the measured data is reliable.
Variational formulation of hybrid problems for fully 3-D transonic flow with shocks in rotor
NASA Technical Reports Server (NTRS)
Liu, Gao-Lian
1991-01-01
Based on previous research, the unified variable domain variational theory of hybrid problems for rotor flow is extended to fully 3-D transonic rotor flow with shocks, unifying and generalizing the direct and inverse problems. Three variational principles (VP) families were established. All unknown boundaries and flow discontinuities (such as shocks, free trailing vortex sheets) are successfully handled via functional variations with variable domain, converting almost all boundary and interface conditions, including the Rankine Hugoniot shock relations, into natural ones. This theory provides a series of novel ways for blade design or modification and a rigorous theoretical basis for finite element applications and also constitutes an important part of the optimal design theory of rotor bladings. Numerical solutions to subsonic flow by finite elements with self-adapting nodes given in Refs., show good agreement with experimental results.
3D-PTV measurement of the phototactic movement of algae in shear flow
NASA Astrophysics Data System (ADS)
Maeda, Tatsuyuki; Ishikawa, Takuji; Ueno, Hironori; Numayama-Tsuruta, Keiko; Imai, Yosuke; Yamaguchi, Takami
2012-11-01
Recently, swimming motion of algae cells is researched actively, because algae fuel is one of the hottest topic in engineering. It is known that algae swim toward the light for photosynthesis however, the effect of a background flow on the unidirectional swimming is unclear. In this study, we used Volvox as a model alga and placed them in a simple shear flow with or without light stimulus. The shear flow was generated by moving two flat sheets in the opposite direction tangentially. A red LED light (wave length 660 nm) was used as an observation light source, and a white LED light was used to stimulate cells for the phototaxis. The trajectories of individual cells were measured by a 3D-PTV system, consists of a pair of high-speed camera with macro lenses. The results were analyzed to understand the effect of the background shear flow on the phototaxis of cells.
Incorporating preferential flow into a 3D model of a forested headwater catchment
NASA Astrophysics Data System (ADS)
Glaser, Barbara; Jackisch, Conrad; Hopp, Luisa; Pfister, Laurent; Klaus, Julian
2016-04-01
Preferential flow plays an important role for water flow and solute transport. The inclusion of preferential flow, for example with dual porosity or dual permeability approaches, is a common feature in transport simulations at the plot scale. But at hillslope and catchment scales, incorporation of macropore and fracture flow into distributed hydrologic 3D models is rare, often due to limited data availability for model parameterisation. In this study, we incorporated preferential flow into an existing 3D integrated surface subsurface hydrologic model (HydroGeoSphere) of a headwater region (6 ha) of the forested Weierbach catchment in western Luxembourg. Our model philosophy was a strong link between measured data and the model setup. The model setup we used previously had been parameterised and validated based on various field data. But existing macropores and fractures had not been considered in this initial model setup. The multi-criteria validation revealed a good model performance but also suggested potential for further improvement by incorporating preferential flow as additional process. In order to pursue the data driven model philosophy for the implementation of preferential flow, we analysed the results of plot scale bromide sprinkling and infiltration experiments carried out in the vicinity of the Weierbach catchment. Three 1 sqm plots were sprinkled for one hour and excavated one day later for bromide depth profile sampling. We simulated these sprinkling experiments at the soil column scale, using the parameterisation of the base headwater model extended by a second permeability domain. Representing the bromide depth profiles was successful without changing this initial parameterisation. Moreover, to explain the variability between the three bromide depth profiles it was sufficient to adapt the dual permeability properties, indicating the spatial heterogeneity of preferential flow. Subsequently, we incorporated the dual permeability simulation in the
3D imaging of particle-scale rotational motion in cyclically driven granular flows
NASA Astrophysics Data System (ADS)
Harrington, Matt; Powers, Dylan; Cooper, Eric; Losert, Wolfgang
Recent experimental advances have enabled three-dimensional (3D) imaging of motion, structure, and failure within granular systems. 3D imaging allows researchers to directly characterize bulk behaviors that arise from particle- and meso-scale features. For instance, segregation of a bidisperse system of spheres under cyclic shear can originate from microscopic irreversibilities and the development of convective secondary flows. Rotational motion and frictional rotational coupling, meanwhile, have been less explored in such experimental 3D systems, especially under cyclic forcing. In particular, relative amounts of sliding and/or rolling between pairs of contacting grains could influence the reversibility of both trajectories, in terms of both position and orientation. In this work, we apply the Refractive Index Matched Scanning technique to a granular system that is cyclically driven and measure both translational and rotational motion of individual grains. We relate measured rotational motion to resulting shear bands and convective flows, further indicating the degree to which pairs and neighborhoods of grains collectively rotate.
On the growth of enstrophy in axisymmetric 3D Euler flows with swirl
NASA Astrophysics Data System (ADS)
Ayala, Diego; Doering, Charles
2015-11-01
By numerically solving suitable PDE-constrained optimization problems, we study a family of axisymmetric vector fields, having the structure of a pair of vortex rings with swirl, that maximize the instantaneous production of enstrophy in the context of 3-dimensional (3D) incompressible Euler flows. The axisymmetric fields are parametrized by their energy K , enstrophy E and helicity H . The imposed symmetry is justified by the results from the seminal work of Doering & Lu (2008), recently confirmed independently by Ayala & Protas (2015), where highly localized pairs of colliding vortex rings were found to be instantaneously optimal for enstrophy production in 3D Navier-Stokes flows. The axial symmetry allows for an exhaustive exploration of the parameter space (K , E , H , as the 3D problem is effectively reduced to a 2-dimensional system of partial differential equations for the modified azimuthal vorticity and the azimuthal circulation density, with the corresponding reduction in computational complexity. Possible connections between these optimal axisymmetric fields with swirl and the ``blow-up'' problem are discussed.
Intrathoracic tumour motion estimation from CT imaging using the 3D optical flow method
NASA Astrophysics Data System (ADS)
Guerrero, Thomas; Zhang, Geoffrey; Huang, Tzung-Chi; Lin, Kang-Ping
2004-09-01
The purpose of this work was to develop and validate an automated method for intrathoracic tumour motion estimation from breath-hold computed tomography (BH CT) imaging using the three-dimensional optical flow method (3D OFM). A modified 3D OFM algorithm provided 3D displacement vectors for each voxel which were used to map tumour voxels on expiration BH CT onto inspiration BH CT images. A thoracic phantom and simulated expiration/inspiration BH CT pairs were used for validation. The 3D OFM was applied to the measured inspiration and expiration BH CT images from one lung cancer and one oesophageal cancer patient. The resulting displacements were plotted in histogram format and analysed to provide insight regarding the tumour motion. The phantom tumour displacement was measured as 1.20 and 2.40 cm with full-width at tenth maximum (FWTM) for the distribution of displacement estimates of 0.008 and 0.006 cm, respectively. The maximum error of any single voxel's motion estimate was 1.1 mm along the z-dimension or approximately one-third of the z-dimension voxel size. The simulated BH CT pairs revealed an rms error of less than 0.25 mm. The displacement of the oesophageal tumours was nonuniform and up to 1.4 cm, this was a new finding. A lung tumour maximum displacement of 2.4 cm was found in the case evaluated. In conclusion, 3D OFM provided an accurate estimation of intrathoracic tumour motion, with estimated errors less than the voxel dimension in a simulated motion phantom study. Surprisingly, oesophageal tumour motion was large and nonuniform, with greatest motion occurring at the gastro-oesophageal junction. Presented at The IASTED Second International Conference on Biomedical Engineering (BioMED 2004), Innsbruck, Austria, 16-18 February 2004.
3D FEM Geometry and Material Flow Optimization of Porthole-Die Extrusion
Ceretti, Elisabetta; Mazzoni, Luca; Giardini, Claudio
2007-05-17
The aim of this work is to design and to improve the geometry of a porthole-die for the production of aluminum components by means of 3D FEM simulations. In fact, the use of finite element models will allow to investigate the effects of the die geometry (webs, extrusion cavity) on the material flow and on the stresses acting on the die so to reduce the die wear and to improve the tool life. The software used to perform the simulations was a commercial FEM code, Deform 3D. The technological data introduced in the FE model have been furnished by METRA S.p.A. Company, partner in this research. The results obtained have been considered valid and helpful by the Company for building a new optimized extrusion porthole-die.
Chaotic advection in blood flow.
Schelin, A B; Károlyi, Gy; de Moura, A P S; Booth, N A; Grebogi, C
2009-07-01
In this paper we argue that the effects of irregular chaotic motion of particles transported by blood can play a major role in the development of serious circulatory diseases. Vessel wall irregularities modify the flow field, changing in a nontrivial way the transport and activation of biochemically active particles. We argue that blood particle transport is often chaotic in realistic physiological conditions. We also argue that this chaotic behavior of the flow has crucial consequences for the dynamics of important processes in the blood, such as the activation of platelets which are involved in the thrombus formation. PMID:19658798
Calculation by the finite element method of 3-D turbulent flow in a centrifugal pump
NASA Astrophysics Data System (ADS)
Combes, J. F.
1992-02-01
In order to solve industrial flow problems in complex geometries, a finite element code, N3S, was developed. It allows the computation of a wide variety of 2-D or 3-D unsteady incompressible flows, by solving the Reynolds averaged Navier-Stokes equations together with a k-epsilon turbulence model. Some recent developments of this code concern turbomachinery flows, where one has to take into account periodic boundary conditions, as well as Coriolis and centrifugal forces. The numerical treatment is based on a fractional step method: at each time step, an advection step is solved successively by means of a characteristic method; a diffusion step for the scalar terms; and finally, a Generalized Stokes Problem by using a preconditioned Uzawa algorithm. The space discretization uses a standard Galerkin finite element method with a mixed formulation for the velocity and pressure. An application is presented of this code to the flow inside a centrifugal pump which was extensively tested on several air and water test rigs, and for which many quasi-3-D or Euler calculations were reported. The present N3S calculation is made on a finite element mesh comprising about 28000 tetrahedrons and 43000 nodes.
Cauchy's almost forgotten Lagrangian formulation of the Euler equation for 3D incompressible flow
NASA Astrophysics Data System (ADS)
Frisch, Uriel; Villone, Barbara
2014-09-01
Two prized papers, one by Augustin Cauchy in 1815, presented to the French Academy and the other by Hermann Hankel in 1861, presented to Göttingen University, contain major discoveries on vorticity dynamics whose impact is now quickly increasing. Cauchy found a Lagrangian formulation of 3D ideal incompressible flow in terms of three invariants that generalize to three dimensions the now well-known law of conservation of vorticity along fluid particle trajectories for two-dimensional flow. This has very recently been used to prove analyticity in time of fluid particle trajectories for 3D incompressible Euler flow and can be extended to compressible flow, in particular to cosmological dark matter. Hankel showed that Cauchy's formulation gives a very simple Lagrangian derivation of the Helmholtz vorticity-flux invariants and, in the middle of the proof, derived an intermediate result which is the conservation of the circulation of the velocity around a closed contour moving with the fluid. This circulation theorem was to be rediscovered independently by William Thomson (Kelvin) in 1869. Cauchy's invariants were only occasionally cited in the 19th century - besides Hankel, foremost by George Stokes and Maurice Lévy - and even less so in the 20th until they were rediscovered via Emmy Noether's theorem in the late 1960, but reattributed to Cauchy only at the end of the 20th century by Russian scientists.
NASA Astrophysics Data System (ADS)
Bocanegra, Humberto; Gorumlu, Seder; Aksak, Burak; Castillo, Luciano; Sheng, Jian
2015-11-01
Understanding how fluid flow interacts with micro-textured surfaces is crucial for a broad range of key biological processes and engineering applications including particle dispersion, pathogenic infections, and drag manipulation by surface topology. Existing methods, such as μPIV, suffers from low spatial resolution and fail to track tracer particle motion very close to a rough surface and within roughness elements. In this paper, we present a technique that combines high speed digital holographic microscopy (DHM) with a correlation based de-noising algorithm to overcome the optical interference generated by surface roughness and to capture a large number of 3D particle trajectories. It allows us to obtain a 3D velocity field with an uncertainty of 0.01% and 2D wall shear stress distribution at the resolution of ~ 65 μPa. Applying the technique to a microfluidics with a surface textured by microfibers, we find that the flow is three-dimensional and complex. While the microfibers affect the velocity flow field locally, their presence is felt globally in terms of wall shear stresses. The study of effect of microfiber patterns and flow characteristics on skin frictions are ongoing and will be reported.
A 3-D nonisothermal flow simulation and pulling force model for injection pultrusion processes
NASA Astrophysics Data System (ADS)
Mustafa, Ibrahim
1998-12-01
Injected Pultrusion (IP) is an efficient way of producing high quality, low cost, high volume and constant cross-section polymeric composites. This process has been developed recently, and the efforts to optimize it are still underway. This work is related to the development of a 3-D non-isothermal flow model for the IP processes. The governing equations for transport of mass, momentum and, energy are formulated by using a local volume averaging approach, and the Finite Element/Control Volume method is used to solve the system of equations numerically. The chemical species balance equation is solved in the Lagrangian frame of reference whereas the energy equation is solved using Galerkin, SU (Streamline Upwind), and SUPG (Streamline Upwind Petrov Galerkin) approaches. By varying degrees of freedom and the flow rates of the resin, it is shown that at high Peclet numbers the SUPG formulation performs better than the SU and the Galerkin methods in all cases. The 3-D model predictions for degree of cure and temperature are compared with a one dimensional analytical solution and the results are found satisfactory. Moreover, by varying the Brinkman Number, it is shown that the effect of viscous dissipation is insignificant. The 3-D flow simulations have been carried out for both thin and thick parts and the results are compared with the 2-D model. It is shown that for thick parts 2-D simulations render erroneous results. The effect of changing permeability on the flow fronts is also addressed. The effect of increasing taper angle on the model prediction is also investigated. A parametric study is conducted to isolate optimum conditions for both isothermal and non-isothermal cases using a straight rectangular die and a die with a tapered inlet. Finally, a simple pulling force model is developed and the pulling force required to pull the carbon-epoxy fiber resin system is estimated for dies of varying tapered inlet.
North Cascadia heat flux and fluid flow from gas hydrates: Modeling 3-D topographic effects
NASA Astrophysics Data System (ADS)
Li, Hong-lin; He, Tao; Spence, George D.
2014-01-01
The bottom-simulating reflector (BSR) of gas hydrate is well imaged from two perpendicular seismic grids in the region of a large carbonate mound, informally called Cucumber Ridge off Vancouver Island. We use a new method to calculate 3-D heat flow map from the BSR depths, in which we incorporate 3-D topographic corrections after calibrated by the drilling results from nearby (Integrated) Ocean Drilling Program Site 889 and Site U1327. We then estimate the associated fluid flow by relating it to the topographically corrected heat flux anomalies. In the midslope region, a heat flux anomaly of 1 mW/m2 can be associated with an approximate focused fluid flow rate of 0.09 mm/yr. Around Cucumber Ridge, high rates of focused fluid flow were observed at steep slopes with values more than double the average regional diffusive fluid discharge rate of 0.56 mm/yr. As well, in some areas of relatively flat seafloor, the focused fluid flow rates still exceeded 0.5 mm/yr. On the seismic lines the regions of focused fluid flow were commonly associated with seismic blanking zones above the BSR and sometimes with strong reflectors below the BSR, indicating that the faults/fractures provide high-permeability pathways for fluids to carry methane from BSR depths to the seafloor. These high fluid flow regions cover mostly the western portion of our area with gas hydrate concentration estimations of ~6% based on empirical correlations from Hydrate Ridge in south off Oregon, significantly higher than previously recognized values of ~2.5% in the eastern portion determined from Site U1327.
Predictions of bubbly flows in vertical pipes using two-fluid models in CFDS-FLOW3D code
Banas, A.O.; Carver, M.B.; Unrau, D.
1995-09-01
This paper reports the results of a preliminary study exploring the performance of two sets of two-fluid closure relationships applied to the simulation of turbulent air-water bubbly upflows through vertical pipes. Predictions obtained with the default CFDS-FLOW3D model for dispersed flows were compared with the predictions of a new model (based on the work of Lee), and with the experimental data of Liu. The new model, implemented in the CFDS-FLOW3D code, included additional source terms in the {open_quotes}standard{close_quotes} {kappa}-{epsilon} transport equations for the liquid phase, as well as modified model coefficients and wall functions. All simulations were carried out in a 2-D axisymmetric format, collapsing the general multifluid framework of CFDS-FLOW3D to the two-fluid (air-water) case. The newly implemented model consistently improved predictions of radial-velocity profiles of both phases, but failed to accurately reproduce the experimental phase-distribution data. This shortcoming was traced to the neglect of anisotropic effects in the modelling of liquid-phase turbulence. In this sense, the present investigation should be considered as the first step toward the ultimate goal of developing a theoretically sound and universal CFD-type two-fluid model for bubbly flows in channels.
A Computational Model for Suspended Large Rigid Bodies in 3D Unsteady Viscous Flows
NASA Astrophysics Data System (ADS)
Xiao, Feng
1999-11-01
A 3D numerical model for computing large rigid objects suspended in fluid flow has been developed. Rather than calculating the surface pressure upon the solid body, we evaluate the net force and torque based on a volume force formulation. The total effective force is obtained by summing up the forces at the Eulerian grids occupied by the rigid body. The effects of the moving bodies are coupled to the fluid flow by imposing the velocity field of the bodies to the fluid. A Poisson equation is used to compute the pressure over the whole domain. The objects are identified by color functions and calculated by the PPM scheme and a tangent function transformation which scales the transition region of the computed interface to a compact thickness. The model is then implemented on a parallel computer of distributed memory and validated with Stokes and low Reynolds number flows.
The 3-D CFD modeling of gas turbine combustor-integral bleed flow interaction
NASA Technical Reports Server (NTRS)
Chen, D. Y.; Reynolds, R. S.
1993-01-01
An advanced 3-D Computational Fluid Dynamics (CFD) model was developed to analyze the flow interaction between a gas turbine combustor and an integral bleed plenum. In this model, the elliptic governing equations of continuity, momentum and the k-e turbulence model were solved on a boundary-fitted, curvilinear, orthogonal grid system. The model was first validated against test data from public literature and then applied to a gas turbine combustor with integral bleed. The model predictions agreed well with data from combustor rig testing. The model predictions also indicated strong flow interaction between the combustor and the integral bleed. Integral bleed flow distribution was found to have a great effect on the pressure distribution around the gas turbine combustor.
Flow control on a 3D backward facing ramp by pulsed jets
NASA Astrophysics Data System (ADS)
Joseph, Pierric; Bortolus, Dorian; Grasso, Francesco
2014-06-01
This paper presents an experimental study of flow separation control over a 3D backward facing ramp by means of pulsed jets. Such geometry has been selected to reproduce flow phenomena of interest for the automotive industry. The base flow has been characterised using PIV and pressure measurements. The results show that the classical notchback topology is correctly reproduced. A control system based on magnetic valves has been used to produce the pulsed jets whose properties have been characterised by hot wire anemometry. In order to shed some light on the role of the different parameters affecting the suppression of the slant recirculation area, a parametric study has been carried out by varying the frequency and the momentum coefficient of the jets for several Reynolds numbers. xml:lang="fr"
Mimicking Natural Laminar to Turbulent Flow Transition: A Systematic CFD Study Using PAB3D
NASA Technical Reports Server (NTRS)
Pao, S. Paul; Abdol-Hamid, Khaled S.
2005-01-01
For applied aerodynamic computations using a general purpose Navier-Stokes code, the common practice of treating laminar to turbulent flow transition over a non-slip surface is somewhat arbitrary by either treating the entire flow as turbulent or forcing the flow to undergo transition at given trip locations in the computational domain. In this study, the possibility of using the PAB3D code, standard k-epsilon turbulence model, and the Girimaji explicit algebraic stresses model to mimic natural laminar to turbulent flow transition was explored. The sensitivity of flow transition with respect to two limiters in the standard k-epsilon turbulence model was examined using a flat plate and a 6:1 aspect ratio prolate spheroid for our computations. For the flat plate, a systematic dependence of transition Reynolds number on background turbulence intensity was found. For the prolate spheroid, the transition patterns in the three-dimensional boundary layer at different flow conditions were sensitive to the free stream turbulence viscosity limit, the reference Reynolds number and the angle of attack, but not to background turbulence intensity below a certain threshold value. The computed results showed encouraging agreements with the experimental measurements at the corresponding geometry and flow conditions.
Numerical simulation of unsteady flow characteristics for cavitation around a 3-D hydrofoil
NASA Astrophysics Data System (ADS)
Ahn, S. H.; Xiao, Y. X.; Wang, Z. W.
2015-01-01
At present it is possible to predict more accurately by various numerical methods established for cavitation simulation around a hydrofoil. However, for the solution of the complex unsteady cavity flow, it is still marginal. In this paper, numerical method is adopted to simulate cavitation around 3-D NACA0015 hydrofoil with homogeneous two-phase flow calculation using commercial code CFX-solver with two turbulence models, the standard RNG k-epsilon turbulence model and the modified RNG k-epsilon turbulence model respectively. First, pressure coefficient for non-cavitating flow, time averaged values of unsteady cavity flow around a hydrofoil are verified to simulate more closely to an actual cavity flow. And then frequency analysis is performed with Fast Fourier Transform. The results show that the calculation results with modified RNG k-epsilon turbulence model agree with experimental results in terms of mean cavity length and pressure drop, but the unsteady flow characteristics of oscillating cavitation still deviate slightly in terms of unsteady cavity flow.
Reconstruction and Visualization of Coordinated 3D Cell Migration Based on Optical Flow.
Kappe, Christopher P; Schütz, Lucas; Gunther, Stefan; Hufnagel, Lars; Lemke, Steffen; Leitte, Heike
2016-01-01
Animal development is marked by the repeated reorganization of cells and cell populations, which ultimately determine form and shape of the growing organism. One of the central questions in developmental biology is to understand precisely how cells reorganize, as well as how and to what extent this reorganization is coordinated. While modern microscopes can record video data for every cell during animal development in 3D+t, analyzing these videos remains a major challenge: reconstruction of comprehensive cell tracks turned out to be very demanding especially with decreasing data quality and increasing cell densities. In this paper, we present an analysis pipeline for coordinated cellular motions in developing embryos based on the optical flow of a series of 3D images. We use numerical integration to reconstruct cellular long-term motions in the optical flow of the video, we take care of data validation, and we derive a LIC-based, dense flow visualization for the resulting pathlines. This approach allows us to handle low video quality such as noisy data or poorly separated cells, and it allows the biologists to get a comprehensive understanding of their data by capturing dynamic growth processes in stills. We validate our methods using three videos of growing fruit fly embryos. PMID:26529743
Simulation of the 3D viscoelastic free surface flow by a parallel corrected particle scheme
NASA Astrophysics Data System (ADS)
Jin-Lian, Ren; Tao, Jiang
2016-02-01
In this work, the behavior of the three-dimensional (3D) jet coiling based on the viscoelastic Oldroyd-B model is investigated by a corrected particle scheme, which is named the smoothed particle hydrodynamics with corrected symmetric kernel gradient and shifting particle technique (SPH_CS_SP) method. The accuracy and stability of SPH_CS_SP method is first tested by solving Poiseuille flow and Taylor-Green flow. Then the capacity for the SPH_CS_SP method to solve the viscoelastic fluid is verified by the polymer flow through a periodic array of cylinders. Moreover, the convergence of the SPH_CS_SP method is also investigated. Finally, the proposed method is further applied to the 3D viscoelastic jet coiling problem, and the influences of macroscopic parameters on the jet coiling are discussed. The numerical results show that the SPH_CS_SP method has higher accuracy and better stability than the traditional SPH method and other corrected SPH method, and can improve the tensile instability. Project supported by the Natural Science Foundation of Jiangsu Province, China (Grant Nos. BK20130436 and BK20150436) and the Natural Science Foundation of the Higher Education Institutions of Jiangsu Province, China (Grant No. 15KJB110025).
Segmentation of bone structures in 3D CT images based on continuous max-flow optimization
NASA Astrophysics Data System (ADS)
Pérez-Carrasco, J. A.; Acha-Piñero, B.; Serrano, C.
2015-03-01
In this paper an algorithm to carry out the automatic segmentation of bone structures in 3D CT images has been implemented. Automatic segmentation of bone structures is of special interest for radiologists and surgeons to analyze bone diseases or to plan some surgical interventions. This task is very complicated as bones usually present intensities overlapping with those of surrounding tissues. This overlapping is mainly due to the composition of bones and to the presence of some diseases such as Osteoarthritis, Osteoporosis, etc. Moreover, segmentation of bone structures is a very time-consuming task due to the 3D essence of the bones. Usually, this segmentation is implemented manually or with algorithms using simple techniques such as thresholding and thus providing bad results. In this paper gray information and 3D statistical information have been combined to be used as input to a continuous max-flow algorithm. Twenty CT images have been tested and different coefficients have been computed to assess the performance of our implementation. Dice and Sensitivity values above 0.91 and 0.97 respectively were obtained. A comparison with Level Sets and thresholding techniques has been carried out and our results outperformed them in terms of accuracy.
Computing 3-D steady supersonic flow via a new Lagrangian approach
NASA Technical Reports Server (NTRS)
Loh, C. Y.; Liou, M.-S.
1993-01-01
The new Lagrangian method introduced by Loh and Hui (1990) is extended for 3-D steady supersonic flow computation. Details of the conservation form, the implementation of the local Riemann solver, and the Godunov and the high resolution TVD schemes are presented. The new approach is robust yet accurate, capable of handling complicated geometry and reactions between discontinuous waves. It keeps all the advantages claimed in the 2-D method of Loh and Hui, e.g., crisp resolution for a slip surface (contact discontinuity) and automatic grid generation along the stream.
Exact Relativistic Ideal Hydrodynamical Solutions in (1+3)D with Longitudinal and Transverse Flows
Liao, Jinfeng; Koch, Volker
2009-05-20
A new method for solving relativistic ideal hydrodynamics in (1+3)D is developed. Longitudinal and transverse radial flows are explicitly embedded into the ansatz for velocity field and the hydrodynamic equations are reduced to a single equation for the transverse velocity field only, which is analytically more tractable as compared with the full hydrodynamic equations. As an application we use the method to find analytically all possible solutions whose transverse velocity fields have power dependence on proper time and transverse radius. Possible application to the Relativistic Heavy Ion Collisions and possible generalizations of the method are discussed.
Reactive Flow Modeling of Liquid Explosives via ALE3D/Cheetah Simulations
Kuo, I W; Bastea, S; Fried, L E
2010-03-10
We carried out reactive flow simulations of liquid explosives such as nitromethane using the hydrodynamic code ALE3D coupled with equations of state and reaction kinetics modeled by the thermochemical code Cheetah. The simulation set-up was chosen to mimic cylinder experiments. For pure unconfined nitromethane we find that the failure diameter and detonation velocity dependence on charge diameter are in agreement with available experimental results. Such simulations are likely to be useful for determining detonability and failure behavior for a wide range of experimental conditions and explosive compounds.
Rotary slanted single wire CTA - a useful tool for 3D flows investigations
NASA Astrophysics Data System (ADS)
Jonáš, P.
2013-04-01
The procedure is described of experimental investigation of a statistically stationary generally nonisothermal 3D flow by means of a constant temperature anemometer (CTA) using single slanted heated wire, rotary round the fixed axis. The principle of this procedure is quite clear. The change of the heated wire temperature modifies ratio of CTA sensitivities to temperature and velocity fluctuations. Turning the heated wire through a proper angle changes the sensitivity to components of the instantaneous velocity vector. Some recommendations are presented based on long time experiences, e.g. on the choice of probe, on the probe calibration, to the measurement organization and to the evaluation of results.
Large-eddy simulation of 3D turbulent flow past a complete marine hydrokinetic turbine
NASA Astrophysics Data System (ADS)
Kang, S.; Sotiropoulos, F.
2011-12-01
A high-resolution computational framework was recently developed by Kang et al (Adv. Water Resour., submitted) for simulating three-dimensional (3D), turbulent flow past real-life, complete marine hydrokinetic (MHK) turbine configurations. In this model the complex turbine geometry is resolved by employing the curvilinear immersed boundary (CURVIB) method, which solves the 3D unsteady incompressible Navier-Stokes equations in generalized curvilinear domains with embedded arbitrarily complex, moving and/or stationary immersed boundaries (Ge and Sotiropoulos, 2007). Turbulence is simulated using the large-eddy simulation (LES) approach adapted in the context of the CURVIB method, with a wall model based on solving the simplified boundary layer equations used to reconstruct boundary conditions near all solid surfaces (Kang et al., 2011). The model can resolve the flow patterns generated by the rotor and all stationary components of the turbine as well as the interactions of the flow structures with the channel bed. We apply this model to carry out LES of the flow past the model-size hydrokinetic turbine deployed in the St. Anthony Falls Laboratory main channel. The mean velocities and second-order turbulence statistics measured in the downstream wake using acoustic Doppler velocimetry (ADV) are compared with the LES results. The comparisons show that the computed mean velocities and turbulent stresses are in good agreement with the measurements. The high-resolution LES data are used to explore physically important downstream flow characteristics such as the time-averaged wake structure, recovery of cross-sectionally averaged power potential, near-bed scour potential, etc. This work is supported by Verdant Power.
Mitiche, Amar; Sekkati, Hicham
2006-11-01
This study investigates a variational, active curve evolution method for dense three-dimentional (3D) segmentation and interpretation of optical flow in an image sequence of a scene containing moving rigid objects viewed by a possibly moving camera. This method jointly performs 3D motion segmentation, 3D interpretation (recovery of 3D structure and motion), and optical flow estimation. The objective functional contains two data terms for each segmentation region, one based on the motion-only equation which relates the essential parameters of 3D rigid body motion to optical flow, and the other on the Horn and Schunck optical flow constraint. It also contains two regularization terms for each region, one for optical flow, the other for the region boundary. The necessary conditions for a minimum of the functional result in concurrent 3D-motion segmentation, by active curve evolution via level sets, and linear estimation of each region essential parameters and optical flow. Subsequently, the screw of 3D motion and regularized relative depth are recovered analytically for each region from the estimated essential parameters and optical flow. Examples are provided which verify the method and its implementation. PMID:17063686
NASA Astrophysics Data System (ADS)
Morgan, J. P.; Hasenclever, J.; Shi, C.
2009-12-01
Computational studies of mantle convection face large challenges to obtain fast and accurate solutions for variable viscosity 3d flow. Recently we have been using parallel (MPI-based) MATLAB to more thoroughly explore possible pitfalls and algorithmic improvements to current ‘best-practice’ variable viscosity Stokes and D’Arcy flow solvers. Here we focus on study of finite-element solvers based on a decomposition of the equations for incompressible Stokes flow: Ku + Gp = f and G’u = 0 (K-velocity stiffness matrix, G-discretized gradient operator, G’=transpose(G)-discretized divergence operator) into a single equation for pressure Sp==G’K^-1Gp =G’K^-1f, in which the velocity is also updated as part of each pressure iteration. The outer pressure iteration is solved with preconditioned conjugate gradients (CG) (Maday and Patera, 1989), with a multigrid-preconditioned CG solver for the z=K^-1 (Gq) step of each pressure iteration. One fairly well-known pitfall (Fortin, 1985) is that constant-pressure elements can generate a spurious non-zero flow under a constant body force within non-rectangular geometries. We found a new pitfall when using an iterative method to solve the Kz=y operation in evaluating each G’K^-1Gq product -- even if the residual of the outer pressure equation converges to zero, the discrete divergence of this equation does not correspondingly converge; the error in the incompressibility depends on roughly the square of the tolerance used to solve each Kz=y velocity-like subproblem. Our current best recipe is: (1) Use flexible CG (cf. Notay, 2001) to solve the outer pressure problem. This is analogous to GMRES for a symmetric positive definite problem. It allows use of numerically unsymmetric and/or inexact preconditioners with CG. (2) In this outer-iteration, use an ‘alpha-bar’ technique to find the appropriate magnitude alpha to change the solution in each search direction. This improvement allows a similar iterative tolerance of
Flow-Through Stream Modeling with MODFLOW and MT3D: Certainties and Limitations.
Ben Simon, Rose; Bernard, Stéphane; Meurville, Charles; Rebour, Vincent
2015-01-01
This paper aims to assess MODFLOW and MT3D capabilities for simulating the spread of contaminants from a river exhibiting an unusual relationship with an alluvial aquifer, with the groundwater head higher than the river head on one side and lower on the other (flow-through stream). A series of simulation tests is conducted using a simple hypothetical model so as to characterize and quantify these limitations. Simulation results show that the expected contaminant spread could be achieved with a specific configuration composed of two sets of parameters: (1) modeled object parameters (hydraulic groundwater gradient, hydraulic conductivity values of aquifer and streambed), and (2) modeling parameters (vertical discretization of aquifer, horizontal refinement of stream modeled with River [RIV] package). The influence of these various parameters on simulation results is investigated, and potential complications and errors are identified. Contaminant spread from stream to aquifer is not always reproduced by MT3D due to the RIV package's inability to simulate lateral exchange fluxes between stream and aquifer. This paper identifies the need for a MODFLOW streamflow package allowing lateral stream-aquifer interactions and streamflow routine calculations. Such developments could be of particular interest for modeling contaminated flow-through streams. PMID:25557038
Coupled aeroelastic oscillations of a turbine blade row in 3D transonic flow
NASA Astrophysics Data System (ADS)
Gnesin, Vitaly; Kolodyazhnaya, Lyubov; Rzadkowski, Romuald
2001-10-01
This paper presents the mutual time - marching method to predict the aeroelastic stability of an oscillating blade row in 3D transonic flow. The ideal gas flow through a blade row is governed by the time dependent Euler equations in conservative form which are integrated by using the explicit monotonous second order accurate Godunov-Kolgan finite volume scheme and moving hybrid H-O grid. The structure analysis uses the modal approach and 3D finite element dynamic model of blade. The blade movement is assumed as a linear combination of the first modes of blade natural oscillations with the modal coefficients depending on time. To demonstrate the capability and correctness of the method, two experimentally investigated test cases have been selected, in which the blades had performed tuned harmonic bending or torsional vibrations (The 1st and 4th standard configurations of the “Workshop on Aeroelasticity in Turbomachines” by Bolcs and Fransson, 1986). The calculated results of aeroelastic behaviour of the blade row (4th standard configuration), are presented over a wide frequency range under different start regimes of interblade phase angle.
Fluid flow pathways study from the 3D seismic data offshore southwestern Taiwan
NASA Astrophysics Data System (ADS)
Chen, L.; Chi, W. C.; Chiang, H. T.; Lin, S.
2014-12-01
3D seismic reflection data provide detailed information on the physical properties of the crust, which can be used for hydrocarbon exploration. Recently, scientists from Taiwan and Germany are collaborating on a project to use a portable 3D seismic system, called P-Cable, to study gas hydrates offshore southwest Taiwan. We have collected 3 cubes, covering the active and passive margins. At these three sites, there is a wide-spread bottom-simulating reflector (BSR). We use the BSR to study the shallow thermal structures of these prospect sites, and use the temperature field information to study fluid migration patterns. We have also done in-situ heat flow measurements, and found similar results, showing focused fluid flow migrations in some pathways. Some of the high temperature fields also correlate with gas chimneys found through seismic attribute analyses. Preliminary results show that there might be active fluid migration above the BSR in the gas hydrate stability zone. In September and October of 2014, we will collect additional P-Cable datasets to be incorporated into this study. Such results will be used to evaluate some proposed sites for future drilling programs.
Real-time tracking with a 3D-Flow processor array
Crosetto, D.
1993-06-01
The problem of real-time track-finding has been performed to date with CAM (Content Addressable Memories) or with fast coincidence logic, because the processing scheme was thought to have much slower performance. Advances in technology together with a new architectural approach make it feasible to also explore the computing technique for real-time track finding thus giving the advantages of implementing algorithms that can find more parameters such as calculate the sagitta, curvature, pt, etc., with respect to the CAM approach. The report describes real-time track finding using new computing approach technique based on the 3D-Flow array processor system. This system consists of a fixed interconnection architecture scheme, allowing flexible algorithm implementation on a scalable platform. The 3D-Flow parallel processing system for track finding is scalable in size and performance by either increasing the number of processors, or increasing the speed or else the number of pipelined stages. The present article describes the conceptual idea and the design stage of the project.
Brightness-compensated 3-D optical flow algorithm for monitoring cochlear motion patterns
NASA Astrophysics Data System (ADS)
von Tiedemann, Miriam; Fridberger, Anders; Ulfendahl, Mats; de Monvel, Jacques Boutet
2010-09-01
A method for three-dimensional motion analysis designed for live cell imaging by fluorescence confocal microscopy is described. The approach is based on optical flow computation and takes into account brightness variations in the image scene that are not due to motion, such as photobleaching or fluorescence variations that may reflect changes in cellular physiology. The 3-D optical flow algorithm allowed almost perfect motion estimation on noise-free artificial sequences, and performed with a relative error of <10% on noisy images typical of real experiments. The method was applied to a series of 3-D confocal image stacks from an in vitro preparation of the guinea pig cochlea. The complex motions caused by slow pressure changes in the cochlear compartments were quantified. At the surface of the hearing organ, the largest motion component was the transverse one (normal to the surface), but significant radial and longitudinal displacements were also present. The outer hair cell displayed larger radial motion at their basolateral membrane than at their apical surface. These movements reflect mechanical interactions between different cellular structures, which may be important for communicating sound-evoked vibrations to the sensory cells. A better understanding of these interactions is important for testing realistic models of cochlear mechanics.
Measuring the orientation and rotation rate of 3D printed particles in turbulent flow
NASA Astrophysics Data System (ADS)
Voth, Greg; Marcus, Guy G.; Parsa, Shima; Kramel, Stefan; Ni, Rui; Cole, Brendan
2014-11-01
The orientation distribution and rotations of anisotropic particles plays a key role in many applications ranging from icy clouds to papermaking and drag reduction in pipe flow. Experimental access to time resolved orientations of anisotropic particles has not been easy to achieve. We have found that 3D printing technology can be used to fabricate a wide range of particle shapes with smallest dimension down to 300 μm. So far we have studied rods, crosses, jacks, tetrads, and helical shapes. We extract the particle orientations from stereoscopic video images using a method of least squares optimization in Euler angle space. We find that in turbulence the orientation and rotation rate of many particles can be understood using a simple picture of alignment of both the vorticity and a long axis of the particle with the Lagrangian stretching direction of the flow. This research is supported by NSF Grant DMR-1208990.
Measuring the orientation and rotation rate of 3D printed particles in turbulent flow
NASA Astrophysics Data System (ADS)
Voth, Greg; Kramel, Stefan; Cole, Brendan
2015-03-01
The orientation distribution and rotations of anisotropic particles plays a key role in many applications ranging from icy clouds to papermaking and drag reduction in pipe flow. Experimental access to time resolved orientations of anisotropic particles has not been easy to achieve. We have found that 3D printing technology can be used to fabricate a wide range of particle shapes with smallest dimension down to 300 ?m. So far we have studied rods, crosses, jacks, tetrads, and helical shapes. We extract the particle orientations from stereoscopic video images using a method of least squares optimization in Euler angle space. We find that in turbulence the orientation and rotation rate of many particles can be understood using a simple picture of alignment of both the vorticity and a long axis of the particle with the Lagrangian stretching direction of the flow.
Decay of the 3D inviscid liquid-gas two-phase flow model
NASA Astrophysics Data System (ADS)
Zhang, Yinghui
2016-06-01
We establish the optimal {Lp-L2(1 ≤ p < 6/5)} time decay rates of the solution to the Cauchy problem for the 3D inviscid liquid-gas two-phase flow model and analyze the influences of the damping on the qualitative behaviors of solution. Compared with the viscous liquid-gas two-phase flow model (Zhang and Zhu in J Differ Equ 258:2315-2338, 2015), our results imply that the friction effect of the damping is stronger than the dissipation effect of the viscosities and enhances the decay rate of the velocity. Our proof is based on Hodge decomposition technique, the {Lp-L2} estimates for the linearized equations and an elaborate energy method.
3D conformation of a flexible fiber in a turbulent flow
NASA Astrophysics Data System (ADS)
Verhille, Gautier; Bartoli, Adrien
2016-07-01
A growing number of studies is devoted to anisotropic particles in turbulent flows. In most cases, the particles are assumed to be rigid and their deformations are neglected. We present an adaptation of classical computer vision tools to reconstruct from two different images the 3D conformation of a fiber distorted by the turbulent fluctuations in a von Kármán flow. This technique allows us notably to characterize the fiber deformation by computing the correlation function of the orientation of the tangent vector. This function allows us to tackle the analogy between polymers and flexible fibers proposed by Brouzet et al. (Phys Rev Lett 112(7):074501, 2014). We show that this function depends on an elastic length ℓ _e which characterizes the particle flexibility, as is the case for polymers, but also on the fiber length L, contrary to polymers.
A numerical solution of 3D inviscid rotational flow in turbines and ducts
NASA Astrophysics Data System (ADS)
Oktay, Erdal; Akmandor, Sinan; Üçer, Ahmet
1998-04-01
The numerical solutions of inviscid rotational (Euler) flows were obtained using an explicit hexahedral unstructured cell vertex finite volume method. A second-order-accurate, one-step Lax-Wendroff scheme was used to solve the unsteady governing equations discretized in conservative form. The transonic circular bump, in which the location and the strength of the captured shock are well predicted, was used as the first test case. The nozzle guide vanes of the VKI low-speed turbine facility were used to validate the Euler code in highly 3D environment. Despite the high turning and the secondary flows which develop, close agreements have been obtained with experimental and numerical results associated with these test cases.
A 3-D implicit finite-volume model of shallow water flows
NASA Astrophysics Data System (ADS)
Wu, Weiming; Lin, Qianru
2015-09-01
A three-dimensional (3-D) model has been developed to simulate shallow water flows in large water bodies, such as coastal and estuarine waters. The eddy viscosity is determined using a newly modified mixing length model that uses different mixing length functions for the horizontal and vertical shear strain rates. The 3-D shallow water flow equations with the hydrostatic pressure assumption are solved using an implicit finite-volume method based on a quadtree (telescoping) rectangular mesh on the horizontal plane and the sigma coordinate in the vertical direction. The quadtree technique can locally refine the mesh around structures or in high-gradient regions by splitting a coarse cell into four child cells. The grid nodes are numbered with a one-dimensional index system that has unstructured grid feature for better grid flexibility. All the primary variables are arranged in a non-staggered grid system. Fluxes at cell faces are determined using a Rhie and Chow-type momentum interpolation, to avoid the possible spurious checkerboard oscillations caused by linear interpolation. Each of the discretized governing equations is solved iteratively using the flexible GMRES method with ILUT preconditioning, and coupling of water level and velocity among these equations is achieved by using the SIMPLEC algorithm with under-relaxation. The model has been tested in four cases, including steady flow near a spur-dyke, tidal flows in San Francisco Bay and Gironde Estuary, and wind-induced current in a flume. The calculated water levels and velocities are in good agreement with the measured values.
Eulerian and Lagrangian methods for vortex tracking in 2D and 3D flows
NASA Astrophysics Data System (ADS)
Huang, Yangzi; Green, Melissa
2014-11-01
Coherent structures are a key component of unsteady flows in shear layers. Improvement of experimental techniques has led to larger amounts of data and requires of automated procedures for vortex tracking. Many vortex criteria are Eulerian, and identify the structures by an instantaneous local swirling motion in the field, which are indicated by closed or spiral streamlines or pathlines in a reference frame. Alternatively, a Lagrangian Coherent Structures (LCS) analysis is a Lagrangian method based on the quantities calculated along fluid particle trajectories. In the current work, vortex detection is demonstrated on data from the simulation of two cases: a 2D flow with a flat plate undergoing a 45 ° pitch-up maneuver and a 3D wall-bounded turbulence channel flow. Vortices are visualized and tracked by their centers and boundaries using Γ1, the Q criterion, and LCS saddle points. In the cases of 2D flow, saddle points trace showed a rapid acceleration of the structure which indicates the shedding from the plate. For channel flow, saddle points trace shows that average structure convection speed exhibits a similar trend as a function of wall-normal distance as the mean velocity profile, and leads to statistical quantities of vortex dynamics. Dr. Jeff Eldredge and his research group at UCLA are gratefully acknowledged for sharing the database of simulation for the current research. This work was supported by the Air Force Office of Scientific Research under AFOSR Award No. FA9550-14-1-0210.
Computed Tomography 3-D Imaging of the Metal Deformation Flow Path in Friction Stir Welding
NASA Technical Reports Server (NTRS)
Schneider, Judy; Beshears, Ronald; Nunes, Arthur C., Jr.
2005-01-01
In friction stir welding (FSW), a rotating threaded pin tool is inserted into a weld seam and literally stirs the edges of the seam together. To determine optimal processing parameters for producing a defect free weld, a better understanding of the resulting metal deformation flow path is required. Marker studies are the principal method of studying the metal deformation flow path around the FSW pin tool. In our study, we have used computed tomography (CT) scans to reveal the flow pattern of a lead wire embedded in a FSW weld seam. At the welding temperature of aluminum, the lead becomes molten and is carried with the macro-flow of the weld metal. By using CT images, a 3-dimensional (3D) image of the lead flow pattern can be reconstructed. CT imaging was found to be a convenient and comprehensive way of collecting and displaying tracer data. It marks an advance over previous more tedious and ambiguous radiographic/metallographic data collection methods.
Test Problems for Reactive Flow HE Model in the ALE3D Code and Limited Sensitivity Study
Gerassimenko, M.
2000-03-01
We document quick running test problems for a reactive flow model of HE initiation incorporated into ALE3D. A quarter percent change in projectile velocity changes the outcome from detonation to HE burn that dies down. We study the sensitivity of calculated HE behavior to several parameters of practical interest where modeling HE initiation with ALE3D.
Viscous Incompressible Flow Computations for 3-D Steady and Unsteady Flows
NASA Technical Reports Server (NTRS)
Kwak, Dochan
2001-01-01
This viewgraph presentation gives an overview of viscous incompressible flow computations for three-dimensional steady and unsteady flows. Details are given on the use of computational fluid dynamics (CFD) as an engineering tool, solution methods for incompressible Navier-Stokes equations, numerical and physical characteristics of the primitive variable approach, and the role of CFD in the past and in current engineering and research applications.
Regulation of pulpal blood flow
Kim, S.
1985-04-01
The regulation of blood flow of the dental pulp was investigated in dogs and rats anesthetized with sodium pentobarbital. Pulpal blood flow was altered by variations of local and systemic hemodynamics. Macrocirculatory blood flow (ml/min/100 g) in the dental pulp was measured with both the /sup 133/Xe washout and the 15-microns radioisotope-labeled microsphere injection methods on the canine teeth of dogs, to provide a comparison of the two methods in the same tooth. Microcirculatory studies were conducted in the rat incisor tooth with microscopic determination of the vascular pattern, RBC velocity, and intravascular volumetric flow distribution. Pulpal resistance vessels have alpha- and beta-adrenergic receptors. Activation of alpha-receptors by intra-arterial injection of norepinephrine (NE) caused both a reduction in macrocirculatory Qp in dogs and decreases in arteriolar and venular diameters and intravascular volumetric flow (Qi) in rats. These responses were blocked by the alpha-antagonist PBZ. Activation of beta-receptors by intra-arterial injection of isoproterenal (ISO) caused a paradoxical reduction of Qp in dogs. In rats, ISO caused a transient increase in arteriolar Qi followed by a flow reduction; arteriolar dilation was accompanied by venular constriction. These macrocirculatory and microcirculatory responses to ISO were blocked by the alpha-antagonist propranolol.
The mantle wedge's transient 3-D flow regime and thermal structure
NASA Astrophysics Data System (ADS)
Davies, D. R.; Le Voci, G.; Goes, S.; Kramer, S. C.; Wilson, C. R.
2016-01-01
Arc volcanism, volatile cycling, mineralization, and continental crust formation are likely regulated by the mantle wedge's flow regime and thermal structure. Wedge flow is often assumed to follow a regular corner-flow pattern. However, studies that incorporate a hydrated rheology and thermal buoyancy predict internal small-scale-convection (SSC). Here, we systematically explore mantle-wedge dynamics in 3-D simulations. We find that longitudinal "Richter-rolls" of SSC (with trench-perpendicular axes) commonly occur if wedge hydration reduces viscosities to Pa s, although transient transverse rolls (with trench-parallel axes) can dominate at viscosities of Pa s. Rolls below the arc and back arc differ. Subarc rolls have similar trench-parallel and trench-perpendicular dimensions of 100-150 km and evolve on a 1-5 Myr time-scale. Subback-arc instabilities, on the other hand, coalesce into elongated sheets, usually with a preferential trench-perpendicular alignment, display a wavelength of 150-400 km and vary on a 5-10 Myr time scale. The modulating influence of subback-arc ridges on the subarc system increases with stronger wedge hydration, higher subduction velocity, and thicker upper plates. We find that trench-parallel averages of wedge velocities and temperature are consistent with those predicted in 2-D models. However, lithospheric thinning through SSC is somewhat enhanced in 3-D, thus expanding hydrous melting regions and shifting dehydration boundaries. Subarc Richter-rolls generate time-dependent trench-parallel temperature variations of up to K, which exceed the transient 50-100 K variations predicted in 2-D and may contribute to arc-volcano spacing and the variable seismic velocity structures imaged beneath some arcs.
Simulation of abrasive flow machining process for 2D and 3D mixture models
NASA Astrophysics Data System (ADS)
Dash, Rupalika; Maity, Kalipada
2015-12-01
Improvement of surface finish and material removal has been quite a challenge in a finishing operation such as abrasive flow machining (AFM). Factors that affect the surface finish and material removal are media viscosity, extrusion pressure, piston velocity, and particle size in abrasive flow machining process. Performing experiments for all the parameters and accurately obtaining an optimized parameter in a short time are difficult to accomplish because the operation requires a precise finish. Computational fluid dynamics (CFD) simulation was employed to accurately determine optimum parameters. In the current work, a 2D model was designed, and the flow analysis, force calculation, and material removal prediction were performed and compared with the available experimental data. Another 3D model for a swaging die finishing using AFM was simulated at different viscosities of the media to study the effects on the controlling parameters. A CFD simulation was performed by using commercially available ANSYS FLUENT. Two phases were considered for the flow analysis, and multiphase mixture model was taken into account. The fluid was considered to be a
Simulation of a Synthetic Jet in Quiescent Air Using TLNS3D Flow Code
NASA Technical Reports Server (NTRS)
Vatsa, Veer N.; Turkel, Eli
2007-01-01
Although the actuator geometry is highly three-dimensional, the outer flowfield is nominally two-dimensional because of the high aspect ratio of the rectangular slot. For the present study, this configuration is modeled as a two-dimensional problem. A multi-block structured grid available at the CFDVAL2004 website is used as a baseline grid. The periodic motion of the diaphragm is simulated by specifying a sinusoidal velocity at the diaphragm surface with a frequency of 450 Hz, corresponding to the experimental setup. The amplitude is chosen so that the maximum Mach number at the jet exit is approximately 0.1, to replicate the experimental conditions. At the solid walls zero slip, zero injection, adiabatic temperature and zero pressure gradient conditions are imposed. In the external region, symmetry conditions are imposed on the side (vertical) boundaries and far-field conditions are imposed on the top boundary. A nominal free-stream Mach number of 0.001 is imposed in the free stream to simulate incompressible flow conditions in the TLNS3D code, which solves compressible flow equations. The code was run in unsteady (URANS) mode until the periodicity was established. The time-mean quantities were obtained by running the code for at least another 15 periods and averaging the flow quantities over these periods. The phase-locked average of flow quantities were assumed to be coincident with their values during the last full time period.
Flow properties along field lines in a 3-D tilted-dipole geometry
NASA Technical Reports Server (NTRS)
Pizzo, V. J.
1995-01-01
A 3-D MHD simulation of a global, tilted-dipole solar wind flow pattern is analyzed to determine flow properties along individual magnetic field lines. In the model, flow conditions near the Sun are chosen to provide a reasonable match to the interplanetary configuration prevailing during the recent south polar passage by Ulysses, i.e., a streamer belt inclined approximately 30 deg to the solar equator and speeds ranging from 325-800 km/s. Field lines all across the stream pattern are traced from 1 to 10 AU by following the motion of marker particles embedded in the flow. It is found that those field lines threading the core of the interaction region are subject to significant latitudinal and relative longitudinal displacement over this range of heliocentric distance. Thus, observations taken at a fixed latitude in the inner solar system sample, over the course of a solar rotation, field lines which connect to a range of latitudes in the outer heliosphere. Maps of the field line displacements are presented to help visualize these connections. In addition, it is found that depending upon the location relative to the CIR structure, the radial evolution of fluid elements frozen to different field lines can deviate considerably from that of the canonical solar wind. That is, for selected subsets of field lines, large speed changes (not just at shocks) can be experienced; the density variation can be far from 1/r(exp 2), and the magnetic field intensity need not decay monotonically with distance.
Simulation of 3D flows past hypersonic vehicles in FlowVision software
NASA Astrophysics Data System (ADS)
Aksenov, A. A.; Zhluktov, S. V.; Savitskiy, D. V.; Bartenev, G. Y.; Pokhilko, V. I.
2015-11-01
A new implicit velocity-pressure split method is discussed in the given presentation. The method implies using conservative velocities, obtained at the given time step, for integration of the momentum equation and other convection-diffusion equations. This enables simulation of super- and hypersonic flows with account of motion of solid boundaries. Calculations of known test cases performed in the FlowVision software are demonstrated. It is shown that the method allows one to carry out calculations at high Mach numbers with integration step essentially exceeding the explicit time step.
NASA Astrophysics Data System (ADS)
Peng, Zhangli; Pak, On Shun; Young, Yuan-Nan; Liu, Allen; Stone, Howard
2015-11-01
We investigate the gating of mechanosensing channels (Mscls) on vesicles and cell membranes under different flow conditions using a multiscale approach. At the cell level (microns), the membrane tension is calculated using a 3D two-component whole-cell membrane model based on dissipative particle dynamics (DPD), including the cortex cytoskeleton and its interactions with the lipid bilayer. At the Mscl level (nanometers), we predict the relation between channel gating and the membrane tension obtained from a cell-level model using a semi-analytical model based on the bilayer hydrophobic mismatch energy. We systematically study the gating of Mscls of vesicles and cell membranes in constricted channel flows and shear flows, and explore the dependence of the gating on flow rate, cell shape and size. The results provide guidance for future experiments in inducing Mscl opening for various purposes such as drug delivery.
3D Structure and Internal Circulation of Pancake Vortices in Rotating Stratified Flows
NASA Astrophysics Data System (ADS)
Hassanzadeh, Pedram; Marcus, Philip; Aubert, Oriane; Le Bars, Michael; Le Gal, Patrice
2011-11-01
Jovian vortices, Atlantic meddies, and vortices of the protoplanetrary disks are examples of weakly-forced or unforced long-lived vortices in rotating stratified flows. Knowing the 3D structure and internal circulation of these vortices is essential in understanding their physics, which is not well-understood. For example, the aspect ratio of these vortices has been long thought to be f / N where f is the Coriolis parameter and N is the Brunt-Vaisala frequency. However, our recent theoretical and experimental study has shown that the aspect ratio in fact depends not only on f and N but also on the Rossby number and density mixing inside the vortex. The new scaling law also agrees with the available measurements of the meddies and Jupiter's Great Red Spot. High resolution 3D numerical simulations of the Navier-Stokes equation are carried out to confirm this new scaling law for a slowly (viscously) decaying anticyclonic vortex in which the Rossby number and stratification inside the vortex evolve in time. For a wide range of parameters and different distributions of density anomaly, the secondary circulations within the vortices are studied. The effect of a non-uniform background stratification is investigated, and the small cyclonic vortices that form above and below the anticyclone are studied.
Moon, Inkyu; Yi, Faliu; Lee, Yeon H; Javidi, Bahram; Boss, Daniel; Marquet, Pierre
2013-12-16
Quantitative phase (QP) images of red blood cells (RBCs), which are obtained by off-axis digital holographic microscopy, can provide quantitative information about three-dimensional (3D) morphology of human RBCs and the characteristic properties such as mean corpuscular hemoglobin (MCH) and MCH surface density (MCHSD). In this paper, we investigate modifications of the 3D morphology and MCH in RBCs induced by the period of storage time for the purpose of classification of RBCs with different periods of storage by using off-axis digital holographic microscopy. The classification of RBCs based on the duration of storage is highly relevant because a long storage of blood before transfusion may alter the functionality of RBCs and, therefore, cause complications in patients. To analyze any changes in the 3D morphology and MCH of RBCs due to storage, we use data sets from RBC samples stored for 8, 13, 16, 23, 27, 30, 34, 37, 40, 47, and 57 days, respectively. The data sets consist of more than 3,300 blood cells in eleven classes, with more than 300 blood cells per class. The classes indicate the storage period of RBCs and are listed in chronological order. Using the RBCs donated by healthy persons, the off-axis digital holographic microscopy reconstructs several quantitative phase images of RBC samples stored for eleven different periods. We employ marker-controlled watershed transform to remove the background in the RBC quantitative phase images obtained by the off-axis digital holographic microscopy. More than 300 single RBCs are extracted from the segmented quantitative phase images for each class. Such a large number of RBC samples enable us to obtain statistical distributions of the characteristic properties of RBCs after a specific period of storage. Experimental results show that the 3D morphology of the RBCs, in contrast to MCH, is essentially related to the aging of the RBCs. PMID:24514667
Potential for 3-D hyporheic exchange flow along a succession of pool-riffle sequences
NASA Astrophysics Data System (ADS)
Käser, Daniel; Binley, Andrew; Krause, Stefan; Heathwaite, Louise
2010-05-01
Pool-riffle sequences are key geomorphological features that can influence the ecology of streams by inducing a flow exchange between surface water and groundwater - a process called hyporheic exchange flow (HEF). The objective of this research was to test the suitability of a simple 3-D groundwater model for characterizing HEF induced by pool-riffle sequences that had been the focus of experimental study. Three reaches of 20 m were modelled separately. While the bed topography was surveyed and represented at a high resolution, the permeability distribution referred to a simple conceptual model consisting of two superposed layers. One hypothesis was that, despite its simplicity, the calibrated model would produce an acceptable fit between observed and simulated heads because its permeability structure resembled the natural system. The potential complexity of hyporheic flow patterns is well-known, yet this study highlights the usefulness of a simple conceptual model coupled to mechanistic flow equations for describing HEF in 3-D. The error structure of the calibrated model provides insight into various site-specific features. The root mean square error between computed and observed hydraulic heads (relative to the head drop over the structure) is comparable to other studies with more elaborate permeability structures. After calibration, a sensitivity analysis was conducted in order to determine the influence of permeability contrast between the layers, depth of the permeability interface, and basal flux on three HEF characteristics: residence time, lateral and vertical extent, and total flux. Results indicate that permeability characteristics can affect HEF in different ways. For example, the vertical extent is deepest in homogeneous conditions, whereas the lateral extent is not significantly affected by permeability contrast, or by the depth of the interface between the two layers. Thus bank piezometers may be insufficient to calibrate groundwater models of HEF
An explicit Runge-Kutta method for 3D turbulent incompressible flows
NASA Technical Reports Server (NTRS)
Sung, Chao-Ho; Lin, Cheng-Wen; Hung, C. M.
1988-01-01
A computer code has been developed to solve for the steady-state solution of the 3D incompressible Reynolds-averaged Navier-Stokes equations. The approach is based on the cell-center, central-difference, finite-volume formulation and an explicit one-step, multistage Runge-Kutta time-stepping scheme. The Baldwin-Lomax turbulence model is used. Techniques to accelerate the rate of convergence to a steady-state solution include the preconditioned method, the local time stepping, and the implicit residual smoothing. Improvements in computational efficiency have been demonstrated in several areas. This numerical procedure has been used to simulate the turbulent horseshoe vortex flow around an airfoil/flat-plate juncture.
NASA Astrophysics Data System (ADS)
Matsakos, T.; Chièze, J.-P.; Stehlé, C.; González, M.; Ibgui, L.; de Sá, L.; Lanz, T.; Orlando, S.; Bonito, R.; Argiroffi, C.; Reale, F.; Peres, G.
2014-08-01
The structure and dynamics of young stellar object (YSO) accretion shocks depend strongly on the local magnetic field strength and configuration, as well as on the radiative transfer effects responsible for the energy losses. We present the first 3D YSO shock simulations of the interior of the stream, assuming a uniform background magnetic field, a clumpy infalling gas, and an acoustic energy flux flowing at the base of the chromosphere. We study the dynamical evolution and the post-shock structure as a function of the plasma-beta (thermal pressure over magnetic pressure). We find that a strong magnetic field (~hundreds of Gauss) leads to the formation of fibrils in the shocked gas due to the plasma confinement within flux tubes. The corresponding emission is smooth and fully distinguishable from the case of a weak magnetic field (~tenths of Gauss) where the hot slab demonstrates chaotic motion and oscillates periodically.
Electroosmotic flow through a microparallel channel with 3D wall roughness.
Chang, Long; Jian, Yongjun; Buren, Mandula; Sun, Yanjun
2016-02-01
In this paper, a perturbation method is introduced to study the EOF in a microparallel channel with 3D wall roughness. The corrugations of the two walls are periodic sinusoidal waves of small amplitude in two directions either in phase or half-period out of phase. Based on linearized Poisson-Boltzmann equation, Laplace equation, and the Navier-Stokes equations, the perturbation solutions of velocity, electrical potential, and volume flow rate are obtained. By using numerical computation, the influences of the wall corrugations on the mean velocity are analyzed. The variations of electrical potential, velocity profile, mean velocity, and their dependences on the wave number α and β of wall corrugations in two directions, the nondimensional electrokinetic width K, and the zeta potential ratio between the lower wall and the upper wall ς are analyzed graphically. PMID:26333852
Progress Toward Overset-Grid Moving Body Capability for USM3D Unstructured Flow Solver
NASA Technical Reports Server (NTRS)
Pandyna, Mohagna J.; Frink, Neal T.; Noack, Ralph W.
2005-01-01
A static and dynamic Chimera overset-grid capability is added to an established NASA tetrahedral unstructured parallel Navier-Stokes flow solver, USM3D. Modifications to the solver primarily consist of a few strategic calls to the Donor interpolation Receptor Transaction library (DiRTlib) to facilitate communication of solution information between various grids. The assembly of multiple overlapping grids into a single-zone composite grid is performed by the Structured, Unstructured and Generalized Grid AssembleR (SUGGAR) code. Several test cases are presented to verify the implementation, assess overset-grid solution accuracy and convergence relative to single-grid solutions, and demonstrate the prescribed relative grid motion capability.
Generation of 3D Spatially Variable Anisotropy for Groundwater Flow Simulations.
Borghi, Andrea; Renard, Philippe; Courrioux, Gabriel
2015-01-01
Sedimentary units generally present anisotropy in their hydraulic properties, with higher hydraulic conductivity along bedding planes, rather than perpendicular to them. This common property leads to a modeling challenge if the sedimentary structure is folded. In this paper, we show that the gradient of the geological potential used by implicit geological modeling techniques can be used to compute full hydraulic conductivity tensors varying in space according to the geological orientation. For that purpose, the gradient of the potential, a vector normal to the bedding, is used to construct a rotation matrix that allows the estimation of the 3D hydraulic conductivity tensor in a single matrix operation. A synthetic 2D cross section example is used to illustrate the method and show that flow simulations performed in such a folded environment are highly influenced by this rotating anisotropy. When using the proposed method, the streamlines follow very closely the folded formation. This is not the case with an isotropic model. PMID:25648610
Multilevel local refinement and multigrid methods for 3-D turbulent flow
Liao, C.; Liu, C.; Sung, C.H.; Huang, T.T.
1996-12-31
A numerical approach based on multigrid, multilevel local refinement, and preconditioning methods for solving incompressible Reynolds-averaged Navier-Stokes equations is presented. 3-D turbulent flow around an underwater vehicle is computed. 3 multigrid levels and 2 local refinement grid levels are used. The global grid is 24 x 8 x 12. The first patch is 40 x 16 x 20 and the second patch is 72 x 32 x 36. 4th order artificial dissipation are used for numerical stability. The conservative artificial compressibility method are used for further improvement of convergence. To improve the accuracy of coarse/fine grid interface of local refinement, flux interpolation method for refined grid boundary is used. The numerical results are in good agreement with experimental data. The local refinement can improve the prediction accuracy significantly. The flux interpolation method for local refinement can keep conservation for a composite grid, therefore further modify the prediction accuracy.
Multigrid direct numerical simulation of the whole process of flow transition in 3-D boundary layers
NASA Technical Reports Server (NTRS)
Liu, Chaoqun; Liu, Zhining
1993-01-01
A new technology was developed in this study which provides a successful numerical simulation of the whole process of flow transition in 3-D boundary layers, including linear growth, secondary instability, breakdown, and transition at relatively low CPU cost. Most other spatial numerical simulations require high CPU cost and blow up at the stage of flow breakdown. A fourth-order finite difference scheme on stretched and staggered grids, a fully implicit time marching technique, a semi-coarsening multigrid based on the so-called approximate line-box relaxation, and a buffer domain for the outflow boundary conditions were all used for high-order accuracy, good stability, and fast convergence. A new fine-coarse-fine grid mapping technique was developed to keep the code running after the laminar flow breaks down. The computational results are in good agreement with linear stability theory, secondary instability theory, and some experiments. The cost for a typical case with 162 x 34 x 34 grid is around 2 CRAY-YMP CPU hours for 10 T-S periods.
Ergodic theory and experimental visualization of chaos in 3D flows
NASA Astrophysics Data System (ADS)
Sotiropoulos, Fotis; Mezic, Igor
2000-11-01
In his motivation for the ergodic hypothesis Gibbs invoked an analogy with fluid mixing: “…Yet no fact is more familiar to us than that stirring tends to bring a liquid to a state of uniform mixture, or uniform densities of its components…”. Although proof of the ergodic hypothesis is possible only for the simplest of systems using methods from ergodic theory, the use of the hypothesis has led to many accurate predictions in statistical mechanics. The problem of fluid mixing, however, turned out to be considerably more complicated than Gibbs envisioned. Chaotic advection can indeed lead to efficient mixing even in non-turbulent flows, but many non-mixed islands are known to persist within well-mixed regions. In numerical studies, Poincaré maps can be used to reveal the structure of such islands but their visualization in the laboratory requires laborious experimental procedures and is possible only for certain types of flows. Here we propose the first non-intrusive, simple to implement, and generally applicable technique for constructing experimental Poincaré maps and apply it to a steady, 3D, vortex breakdown bubble. We employ standard laser-induced fluorescence (LIF) and construct Poincaré maps by time averaging a sufficiently long sequence of instantaneous LIF images. We also show that ergodic theory methods provide a rigorous theoretical justification for this approach whose main objective is to reveal the non-ergodic regions of the flow.
Magnetic Damping of g-Jitter Driven Flows: 3-D Calculations
NASA Technical Reports Server (NTRS)
Shang, D. Y.; Li, B. Q.; deGroh, H. C.
1997-01-01
A 3-D numerical model is developed to represent the oscillating natural convection induced in a cylindrical cavity filled with Ga-doped germanium with and without the presence of an external magnetic field. The model is developed based on the penalty-finite element solution of the equations describing the transport of momentum, heat and solutal element as well as the electromagnetic field distribution in the melt pool. Automatic time step control is applied to help speed up the calculations. Numerical simulations are conducted to study the convection and magnetic damping effects as a function of frequency, directions and amplitudes of g-jitter and also the direction and magnitudes of the applied magnetic fields. The results show that the g-jitter driven flow is time dependent and exhibits a complex recirculating convection pattern in three dimensions and that an applied magnetic field can be employed to suppress this deleterious convective flow and both magnitude and orientation of the applied field are important in magnetic damping of the g-jitter induced convective flows.
Dynamic coupling between fluid flow and vein growth in fractures: a 3D numerical model
NASA Astrophysics Data System (ADS)
Schwarz, J.-O.; Enzmann, F.
2012-04-01
Fluid flow is one of the main mass transport mechanisms in the Earth's crust and abundant mineral vein networks are important indicators for fluid flow and fluid rock interaction. These systems are dynamic and part of the so called RTM processes (reaction-transport-mechanics). Understanding of mineral vein systems requires coupling of these processes. Here we present a conceptional model for dynamic vein growth of syntaxial, posttectonic veins generated by advective fluid flow and show first results of a numerical model for this scenario. Vein generation requires three processes to occur: (i) fracture generation by mechanical stress e.g. hydro-fracturing, (ii) flow of a supersaturated fluid on that fracture and (iii) crystallization of phase(s) on or in the fracture. 3D synthetic fractures are generated with the SynFrac code (Ogilvie, et al. 2006). Subsequently solutions of the Navier-Stokes equation for this fracture are computed by a computational fluid dynamics code called GeoDict (Wiegmann 2007). Transport (advective and diffusive) of chemical species to growth sites in the fracture and vein growth are computed by a self-written MATLAB script. The numerical model discretizes the wall rock and fracture geometry by volumetric pixels (voxels). Based on this representation, the model computes the three basic functions for vein generation: (a) nucleation, (b) fluid flow with transport of chemical species and (c) growth. The following conditions were chosen for these three modules. Nucleation is heterogeneous and occurs instantaneously at the wall rock/fracture interface. Advective and diffusive flow of a supersaturated fluid and related transport of chemical species occurs according to the computed fluid flow field by GeoDict. Concentration of chemical species at the inflow is constant, representing external fluid buffering. Changes/decrease in the concentration of chemical species occurs only due to vein growth. Growth of nuclei is limited either by transport of
Maji, Debnath; Suster, Michael A; Stavrou, Evi; Gurkan, Umut A; Mohseni, Pedram
2015-08-01
This paper reports on the design, fabrication, and testing of a microfluidic sensor for dielectric spectroscopy (DS) of human whole blood during coagulation. The sensor employs a three-dimensional (3D), parallel-plate, capacitive sensing structure with a floating electrode integrated into a microfluidic channel. Using an impedance analyzer and after a 5-point calibration, the sensor is shown to measure the real part of complex relative dielectric permittivity of human whole blood in a frequency range of 10kHz to 100MHz. The temporal variation of dielectric permittivity at 1MHz for human whole blood from three different healthy donors shows a peak in permittivity at ~ 4 to 5 minutes, which also corresponds to the onset of CaCl2-initiated coagulation of the blood sample verified visually. PMID:26737635
Enhancement of USM3D Unstructured Flow Solver for High-Speed High-Temperature Shear Flows
NASA Technical Reports Server (NTRS)
Pandya, Mohagna J.; Abdol-Hamid, Khaled S.; Frink, Neal T.
2009-01-01
Large temperature and pressure fluctuations have a profound effect on turbulence development in transonic and supersonic jets. For high-speed, high-temperature jet flows, standard turbulence models lack the ability to predict the observed mixing rate of a shear layer. Several proposals to address this deficiency have been advanced in the literature to modify the turbulence transport equations in a variety of ways. In the present study, some of the most proven and simple modifications to two-equation turbulence models have been selected and implemented in NASA's USM3D tetrahedral Navier-Stokes flow solver. The modifications include the addition of compressibility correction and pressure dilatation terms in the turbulence transport equations for high-speed flows, and the addition of a simple modification to the Boussinesq's closure model coefficient for high-temperature jets. The efficacy of the extended models is demonstrated by comparison with experimental data for two supersonic axisymmetric jet test cases at design pressure ratio.
Investigation of the effect of wall friction on the flow rate in 2D and 3D Granular Flow
NASA Astrophysics Data System (ADS)
Carballo-Ramirez, Brenda; Pleau, Mollie; Easwar, Nalini; Birwa, Sumit; Shah, Neil; Tewari, Shubha
We have measured the mass flow rate of spherical steel spheres under gravity in vertical, straight-walled 2 and 3-dimensional hoppers, where the flow velocity is controlled by the opening size. Our measurements focus on the role of friction and its placement along the walls of the hopper. In the 2D case, an increase in the coefficient of static friction from μ = 0.2 to 0.6 is seen to decrease the flow rate significantly. We have changed the placement of frictional boundaries/regions from the front and back walls of the 2D hopper to the side walls and floor to investigate the relative importance of the different regions in determining the flow rate. Fits to the Beverloo equation show significant departure from the expected exponent of 1.5 in the case of 2D flow. In contrast, 3D flow rates do not show much dependence on wall friction and its placement. We compare the experimental data to numerical simulations of gravity driven hopper granular flow with varying frictional walls constructed using LAMMPS*. *http://lammps.sandia.gov Supported by NSF MRSEC DMR 0820506.
Secondary flow structure in a model curved artery: 3D morphology and circulation budget analysis
NASA Astrophysics Data System (ADS)
Bulusu, Kartik V.; Plesniak, Michael W.
2015-11-01
In this study, we examined the rate of change of circulation within control regions encompassing the large-scale vortical structures associated with secondary flows, i.e. deformed Dean-, Lyne- and Wall-type (D-L-W) vortices at planar cross-sections in a 180° curved artery model (curvature ratio, 1/7). Magnetic resonance velocimetry (MRV) and particle image velocimetry (PIV) experiments were performed independently, under the same physiological inflow conditions (Womersley number, 4.2) and using Newtonian blood-analog fluids. The MRV-technique performed at Stanford University produced phase-averaged, three-dimensional velocity fields. Secondary flow field comparisons of MRV-data to PIV-data at various cross-sectional planes and inflow phases were made. A wavelet-decomposition-based approach was implemented to characterize various secondary flow morphologies. We hypothesize that the persistence and decay of arterial secondary flow vortices is intrinsically related to the influence of the out-of-plane flow, tilting, in-plane convection and diffusion-related factors within the control regions. Evaluation of these factors will elucidate secondary flow structures in arterial hemodynamics. Supported by the National Science Foundation under Grant Number CBET-0828903, and GW Center for Biomimetics and Bioinspired Engineering (COBRE). The MRV data were acquired at Stanford University in collaboration with Christopher Elkins and John Eaton.
NASA Astrophysics Data System (ADS)
Meheust, Y.; De Dreuzy, J.; Pichot, G.
2011-12-01
Flow channeling and permeability scaling in fractured media have been classically addressed either at the fracture- or at the network- scales. In the latter case they are linked to the topological structure of the network, while at the fracture scale they are controlled by the variability of the local aperture distribution inside individual fractures. In this study we analyze these two combined effects, investigating how flow localization below the scale of individual fractures influences that at the network scale and the resulting medium permeability. This is done by use of a new highly-resolved 3D discrete fracture network model (DFN). The local apertures of individual fractures are distributed according to a truncated Gaussian law, and exhibit self-affine spatial correlations that are bounded by an upper cutoff scale Lc; Lc and the fracture closure, defined as the ratio of the aperture fluctuations at scale Lc to the mean aperture, are considered homogeneous over the DFN. The network topology is controlled by a homogeneous scalar fracture density and a power law fracture length distribution. We have varied these features to investigate a large variety of DFN topologies, from sparse networks with varying degrees of fracture interconnections, flow bottlenecks and dead-ends (Fig. 1a), to dense well-connected networks (Fig. 1b). We have also investigated a large range of fracture closures, performing extensive simulations of about 105 different DFN realizations. At the fracture scale, accounting for local aperture fluctuations leads to a monotical deviation (which can exceed 50%) of the equivalent fracture transmissivity from the parallel plate behavior. At the network scale we observe a complex interaction between flow channeling within fracture planes and flow localization in the network. This interaction is controlled by the location of fracture interactions with respect to that of low local transmissivity zones (particularly the closed zones), in the fracture
3D Markov Process for Traffic Flow Prediction in Real-Time.
Ko, Eunjeong; Ahn, Jinyoung; Kim, Eun Yi
2016-01-01
Recently, the correct estimation of traffic flow has begun to be considered an essential component in intelligent transportation systems. In this paper, a new statistical method to predict traffic flows using time series analyses and geometric correlations is proposed. The novelty of the proposed method is two-fold: (1) a 3D heat map is designed to describe the traffic conditions between roads, which can effectively represent the correlations between spatially- and temporally-adjacent traffic states; and (2) the relationship between the adjacent roads on the spatiotemporal domain is represented by cliques in MRF and the clique parameters are obtained by example-based learning. In order to assess the validity of the proposed method, it is tested using data from expressway traffic that are provided by the Korean Expressway Corporation, and the performance of the proposed method is compared with existing approaches. The results demonstrate that the proposed method can predict traffic conditions with an accuracy of 85%, and this accuracy can be improved further. PMID:26821025
CFD Simulation of 3D Flow field in a Gas Centrifuge
Dongjun Jiang; Shi Zeng
2006-07-01
A CFD method was used to study the whole flow field in a gas centrifuge. In this paper, the VSM (Vector Splitting Method) of the FVM (Finite Volume Method) was used to solve the 3D Navier-Stokes equations. An implicit second-order upwind scheme was adopted. The numerical simulation was successfully performed on a parallel cluster computer and a convergence result was obtained. The simulation shows that: in the withdrawal chamber, a strong detached shock wave is formed in front of the scoop; as the radial position increases, the shock becomes stronger and the distance to scoop front surface is smaller. An oblique shock forms in the clearance between the scoop and the centrifuge wall; behind the shock-wave, the radially-inward motion of gas is induced because of the imbalance of the pressure gradient and the centrifugal force. In the separation chamber, a countercurrent is introduced. This indicates that CFD method can be used to study the complex three-dimensional flow field of gas centrifuges. (authors)
3D Markov Process for Traffic Flow Prediction in Real-Time
Ko, Eunjeong; Ahn, Jinyoung; Kim, Eun Yi
2016-01-01
Recently, the correct estimation of traffic flow has begun to be considered an essential component in intelligent transportation systems. In this paper, a new statistical method to predict traffic flows using time series analyses and geometric correlations is proposed. The novelty of the proposed method is two-fold: (1) a 3D heat map is designed to describe the traffic conditions between roads, which can effectively represent the correlations between spatially- and temporally-adjacent traffic states; and (2) the relationship between the adjacent roads on the spatiotemporal domain is represented by cliques in MRF and the clique parameters are obtained by example-based learning. In order to assess the validity of the proposed method, it is tested using data from expressway traffic that are provided by the Korean Expressway Corporation, and the performance of the proposed method is compared with existing approaches. The results demonstrate that the proposed method can predict traffic conditions with an accuracy of 85%, and this accuracy can be improved further. PMID:26821025
Adhikari, Deepak; Gemmell, Brad J; Hallberg, Michael P; Longmire, Ellen K; Buskey, Edward J
2015-11-01
We describe an automated, volumetric particle image velocimetry (PIV) and tracking method that measures time-resolved, 3D zooplankton trajectories and surrounding volumetric fluid velocity fields simultaneously and non-intrusively. The method is demonstrated for groups of copepods flowing past a wall-mounted cylinder. We show that copepods execute escape responses when subjected to a strain rate threshold upstream of a cylinder, but the same threshold range elicits no escape responses in the turbulent wake downstream. The method was also used to document the instantaneous slip velocity of zooplankton and the resulting differences in trajectory between zooplankton and non-inertial fluid particles in the unsteady wake flow, showing the method's capability to quantify drift for both passive and motile organisms in turbulent environments. Applications of the method extend to any group of organisms interacting with the surrounding fluid environment, where organism location, larger-scale eddies and smaller-scale fluid deformation rates can all be tracked and analyzed. PMID:26486364
USM3D Simulations of Saturn V Plume Induced Flow Separation
NASA Technical Reports Server (NTRS)
Deere, Karen; Elmlilgui, Alaa; Abdol-Hamid, K. S.
2011-01-01
The NASA Constellation Program included the Ares V heavy lift cargo vehicle. During the design stage, engineers questioned if the Plume Induced Flow Separation (PIFS) that occurred along Saturn V rocket during moon missions at some flight conditions, would also plague the newly proposed rocket. Computational fluid dynamics (CFD) was offered as a tool for initiating the investigation of PIFS along the Ares V rocket. However, CFD best practice guidelines were not available for such an investigation. In an effort to establish a CFD process and define guidelines for Ares V powered simulations, the Saturn V vehicle was used because PIFS flight data existed. The ideal gas, computational flow solver USM3D was evaluated for its viability in computing PIFS along the Saturn V vehicle with F-1 engines firing. Solutions were computed at supersonic freestream conditions, zero degree angle of attack, zero degree sideslip, and at flight Reynolds numbers. The effects of solution sensitivity to grid refinement, turbulence models, and the engine boundary conditions on the predicted PIFS distance along the Saturn V were discussed and compared to flight data from the Apollo 11 mission AS-506.
Ghost Particle Velocimetry: Accurate 3D Flow Visualization Using Standard Lab Equipment
NASA Astrophysics Data System (ADS)
Buzzaccaro, Stefano; Secchi, Eleonora; Piazza, Roberto
2013-07-01
We describe and test a new approach to particle velocimetry, based on imaging and cross correlating the scattering speckle pattern generated on a near-field plane by flowing tracers with a size far below the diffraction limit, which allows reconstructing the velocity pattern in microfluidic channels without perturbing the flow. As a matter of fact, adding tracers is not even strictly required, provided that the sample displays sufficiently refractive-index fluctuations. For instance, phase separation in liquid mixtures in the presence of shear is suitable to be directly investigated by this “ghost particle velocimetry” technique, which just requires a microscope with standard lamp illumination equipped with a low-cost digital camera. As a further bonus, the peculiar spatial coherence properties of the illuminating source, which displays a finite longitudinal coherence length, allows for a 3D reconstruction of the profile with a resolution of few tenths of microns and makes the technique suitable to investigate turbid samples with negligible multiple scattering effects.
A digital holography set-up for 3D vortex flow dynamics
NASA Astrophysics Data System (ADS)
Lebon, Benoît; Perret, Gaële; Coëtmellec, Sébastien; Godard, Gilles; Gréhan, Gérard; Lebrun, Denis; Brossard, Jérôme
2016-06-01
In the present paper, a digital in-line holography (DIH) set-up, with a converging beam, is used to take three-dimensional (3D) velocity measurements of vortices. The vortices are formed periodically at the edges of a submerged horizontal plate submitted to regular waves. They take the form of vortex filaments that extend from side to side of the channel. They undergo strongly three-dimensional instability mechanisms that remain very complicated to characterize experimentally. The experiments are performed in a 10 × 0.3 × 0.3 m3 wave flume. The DIH set-up is performed using a modulated laser diode emitting at the wavelength of 640 nm and a lensless CCD camera. The beam crosses the channel side to side. To reveal the flow dynamics, 30-μm hydrogen bubbles are generated at the edge of the plate to serve as tracers. Their locations are recorded on the holograms multiple times to access the dynamics of the flow. This method leads to an accuracy in the order of 100 μm on the axial location. Those measurements have been validated with stereo-PIV measurements. A very good agreement is found on time-averaged velocity fields between the two techniques.
Numerical Calculations of 3-D High-Lift Flows and Comparison with Experiment
NASA Technical Reports Server (NTRS)
Compton, William B, III
2015-01-01
Solutions were obtained with the Navier-Stokes CFD code TLNS3D to predict the flow about the NASA Trapezoidal Wing, a high-lift wing composed of three elements: the main-wing element, a deployed leading-edge slat, and a deployed trailing-edge flap. Turbulence was modeled by the Spalart-Allmaras one-equation turbulence model. One case with massive separation was repeated using Menter's two-equation SST (Menter's Shear Stress Transport) k-omega turbulence model in an attempt to improve the agreement with experiment. The investigation was conducted at a free stream Mach number of 0.2, and at angles of attack ranging from 10.004 degrees to 34.858 degrees. The Reynolds number based on the mean aerodynamic chord of the wing was 4.3 x 10 (sup 6). Compared to experiment, the numerical procedure predicted the surface pressures very well at angles of attack in the linear range of the lift. However, computed maximum lift was 5% low. Drag was mainly under predicted. The procedure correctly predicted several well-known trends and features of high-lift flows, such as off-body separation. The two turbulence models yielded significantly different solutions for the repeated case.
High-Resolution 3D Seismic Imaging of Fluid Flow Anomalies in the Southwest Barents Sea
NASA Astrophysics Data System (ADS)
Planke, S.; Eriksen, F. N.; Eriksen, O. K.; Assad, M.; Stokke, H. H.
2014-12-01
Fluid flow features imaged as gas flares in the water column, pockmarks and mud volcanoes on the seabed, and high-amplitude cross-cutting reflections and bright spots in the sub-surface are abundant in the SW Barents Sea offshore northern Norway. This region is covered by extensive conventional 2D and 3D deep penetration seismic reflection data and multibeam bathymetry. High-resolution 3D P-Cable seismic data have been acquired in the SW Barents Sea over the past few years to image the uppermost ca. 500 m of the sub-surface. The P-Cable system consist of 12 to 16 short streamers (25 m) that are towed on a cross-cable perpendicular to the vessel's steaming direction. This configuration allows for acquisition of seismic data with high trace density, typically with 6 m in-line separation. The vertical resolution is a good as 1-2 m using conventional site survey air gun configurations. The sedimentary succession in the SW Barents Sea consists of upper Paleozoic evaporites overlaid by Mesozoic and Cenozoic clastic sediments. There are several organic-rich intervals in the sequence, including Paleozoic coals and Triassic and Jurassic marine source rocks that are locally in the oil or gas maturation windows. Glacial erosion has locally removed kilometer thick Cenozoic and Mesozoic successions, leaving the Mesozoic and Paleozoic strata in shallow sub-surface. The new high-resolution 3D surveys have targeted shallow fluid anomalies in the subsurface. These are imaged as high-amplitude reflections in fault blocks and structural highs, locally cross-cutting well-defined Mesozoic reflections. Commonly, disturbed reflections are present in overlying sequences, or high-amplitude reflections are imaged in the glacial overburden sediments. Locally, hundreds of pockmarks are imaged by the seafloor reflection. The deep cross-cutting reflections are interpreted as hydrocarbon accumulations that locally migrate towards the surface. The fluids are stored in shallow gas pockets or
Simulation of bacteria transport processes in a river with Flow3D
NASA Astrophysics Data System (ADS)
Schwarzwälder, Kordula; Bui, Minh Duc; Rutschmann, Peter
2014-05-01
Water quality aspects are getting more and more important due to the European water Framework directive (WFD). One problem related to this topic is the inflow of untreated wastewater due to combined sewer overflows into a river. The wastewater mixture contains even bacteria like E. coli and Enterococci which are markers for water quality. In our work we investigated the transport of these bacteria in river Isar by using a large-scale flume in the outside area of our lab (Oskar von Miller Institute). Therefor we could collect basic data and knowledge about the processes which occur during bacteria sedimentation and remobilisation. In our flume we could use the real grain with the exact size distribution curve as in the river Isar which we want to simulate and we had the chance to nurture a biofilm which is realistic for the analysed situation. This biofilm plays an important role in the remobilisation processes, because the bacteria are hindered to be washed out back into the bulk phase as fast and in such an amount as this would happen without biofilm. The results of our experiments are now used for a module in the 3D software Flow3D to simulate the effects of a point source inlet of raw wastewater on the water quality. Therefor we have to implement the bacteria not as a problem of concentration with advection and diffusion but as single particles which can be inactivated during the process of settling and need to be hindered from remobilisation by the biofilm. This biofilm has special characteristic, it is slippery and has a special thickness which influences the chance of bacteria being removed. To achieve realistic results we have to include the biofilm with more than a probabilistic-tool to make sure that our module is transferable. The module should be as flexible as possible to be improved step by step with increasing quality of dataset.
Ocular Blood Flow Autoregulation Mechanisms and Methods
Luo, Xue; Shen, Yu-meng; Jiang, Meng-nan; Lou, Xiang-feng; Shen, Yin
2015-01-01
The main function of ocular blood flow is to supply sufficient oxygen and nutrients to the eye. Local blood vessels resistance regulates overall blood distribution to the eye and can vary rapidly over time depending on ocular need. Under normal conditions, the relation between blood flow and perfusion pressure in the eye is autoregulated. Basically, autoregulation is a capacity to maintain a relatively constant level of blood flow in the presence of changes in ocular perfusion pressure and varied metabolic demand. In addition, ocular blood flow dysregulation has been demonstrated as an independent risk factor to many ocular diseases. For instance, ocular perfusion pressure plays key role in the progression of retinopathy such as glaucoma and diabetic retinopathy. In this review, different direct and indirect techniques to measure ocular blood flow and the effect of myogenic and neurogenic mechanisms on ocular blood flow are discussed. Moreover, ocular blood flow regulation in ocular disease will be described. PMID:26576295
High-resolution 3D seismic data characterize fluid flow systems in the SW Barents Sea
NASA Astrophysics Data System (ADS)
Bünz, Stefan; Mienert, Jürgen; Rajan, Anupama
2010-05-01
The flow of fluids through marine sediments is one of the most dominant and pervasive processes in continental margins. These processes control the evolution of a sedimentary basin and its seafloor environment, and have implications for hydrocarbon exploration and seabed ecosystems. Many seep sites at the seafloor are associated with large but complex faunal communities that have received significant attention in recent years. However, there is a need for a better understanding of the driving mechanism of fluid flow in various geological settings, the accumulation of fluids in the subsurface and their focused flow through conduits and/or faults to the seabed. The Barents Sea is a large hydrocarbon-prone basin of the Norwegian Arctic region. A significant portion of the hydrocarbons has leaked or migrated into the shallow subsurface and is now trapped in gas-hydrate and shallow-gas reservoirs. Furthermore, there are few places in the Barents Sea, where methane gas is leaking from the seafloor into the oceanosphere. Accumulations of free gas in the shallow subsurface are considered a geohazard. They constitute a risk for safe drilling operations and they may pose a threat to global climate if the seal that is trapping them is breached. P-Cable 3D high-resolution seismic data from the Ringvassøya Fault Complex and the Polheim Sub-Platform provide new and detailed insight into fluid flow controls and accumulation mechanisms. The data shows a wide variety of fluid flow features, mostly in the form of pockmarks, bright spots, wipe-out zones or vertical zones of disturbed reflectivity. Fluids migrate by both diapiric mechanism and channelized along sedimentary layers. Glacigenic sediments generally form a strong boundary for fluid flow in the very shallow section. However, we can recognize pockmarks not only at the seafloor but also at one subsurface layer approximately 50 m below sea floor indicating a former venting period in the SW Barents Sea. At few locations high
NASA Technical Reports Server (NTRS)
Hardy, Terry L.; Tomsik, Thomas M.
1990-01-01
As part of the National Aero-Space Plane (NASP) project, the multi-dimensional effects of gravitational force, initial tank pressure, initial ullage temperature, and heat transfer rate on the 2-D temperature profiles were studied. FLOW-3D, a commercial finite difference fluid flow model, was used for the evaluation. These effects were examined on the basis of previous liquid hydrogen experimental data with gaseous hydrogen pressurant. FLOW-3D results were compared against an existing 1-D model. In addition, the effects of mesh size and convergence criteria on the analytical results were investigated. Suggestions for future modifications and uses of FLOW-3D for modeling of a NASP tank are also presented.
Turbulence modeling for subsonic separated flows over 2-D airfoils and 3-D wings
NASA Astrophysics Data System (ADS)
Rosen, Aaron M.
Accurate predictions of turbulent boundary layers and flow separation through computational fluid dynamics (CFD) are becoming more and more essential for the prediction of loads in the design of aerodynamic flight components. Standard eddy viscosity models used in many commercial codes today do not capture the nonequilibrium effects seen in a separated flow and thus do not generally make accurate separation predictions. Part of the reason for this is that under nonequilibrium conditions such as a strong adverse pressure gradient, the history effects of the flow play an important role in the growth and decay of turbulence. More recent turbulence models such as Olsen and Coakley's Lag model and Lillard's lagRST model seek to simulate these effects by lagging the turbulent variables when nonequilibrium effects become important. The purpose of the current research is to assess how these nonequilibrium turbulence models capture the separated regions on various 2-D airfoils and 3-D wings. Nonequilibrium models including the Lag model and the lagRST model are evaluated in comparison with three baseline models (Spalart-Allmaras, Wilcox's k-omega, and Menter's SST) using a modified version of the OVERFLOW code. Tuning the model coefficients of the Lag and lagRST models is also explored. Results show that the various lagRST formulations display an improvement in velocity profile predictions over the standard RANS models, but have trouble capturing the edge of the boundary layer. Experimental separation location measurements were not available, but several trends are noted which may be useful to tuning the model coefficients in the future.
Experimental Investigation of Material Flows Within FSWs Using 3D Tomography
Charles R. Tolle; Timothy A. White; Karen S. Miller; Denis E. Clark; Herschel B. Smartt
2008-06-01
There exists significant prior work using tracers or pre-placed hardened markers within friction stir welding (FSWing) to experimentally explore material flow within the FSW process. Our experiments replaced markers with a thin sheet of copper foil placed between the 6061 aluminum lap and butt joints that were then welded. The absorption characteristics of x-rays for copper and aluminum are significantly different allowing for non-destructive evaluation (NDE) methods such as x-ray computed tomography (CT) to be used to demonstrate the material movement within the weldment on a much larger scale than previously shown. 3D CT reconstruction of the copper components of the weldment allows for a unique view into the final turbulent state of the welding process as process parameters are varied. The x-ray CT data of a section of the weld region was collected using a cone-beam x-ray imaging system developed at the INL. Six-hundred projections were collected over 360-degrees using a 160-kVp Bremsstrahlung x-ray generator (25-micrometer focal spot) and amorphoussilicon x-ray detector. The region of the object that was imaged was about 3cm tall and 1.5cm x 1cm in cross section, and was imaged at a magnification of about 3.6x. The data were reconstructed on a 0.5x0.5x0.5 mm3 voxel grid. After reconstruction, the aluminum and copper could be easily discriminated using a gray level threshold allowing visualization of the copper components. Fractal analysis of the tomographic reconstructed material topology is investigated as a means to quantify macro level material flow based on process parameters. The results of multi-pass FSWs show increased refinement of the copper trace material. Implications of these techniques for quantifying process flow are discussed.
Heat Flow Partitioning Between Continents and Oceans - from 2D to 3D
NASA Astrophysics Data System (ADS)
Moresi, L. N.; Cooper, C. M.; Lenardic, A.
2010-12-01
Scalings derived from thermal network theory explain how the presence of continents can influence the Earth’s overall heat loss. Intuitively, it may seem that increasing the proportion of a planet’s surface area covered by continents would decrease the efficiency of heat transfer given that continents do not participate in convective overturn. However, this ignores the potential feedback between the insulating effect of continents and the temperature-dependent viscosity of the mantle (Lenardic et al, 2005, Cooper et al, 2007). When this feedback is considered, a clear regime exists in which the partial stagnation and insulation of the surface by buoyant continental crust can lead to an increase in heat flow compared to the uninsulated case. The numerical results used to verify the scalings have mostly been conducted in two dimensions in order to cover a very wide range of Rayleigh number, fraction of continental coverage, and continental thickness. However as more recent results show that the configuration of the crust also plays a role in determining the heat flow partitioning and global heat flow (See Lenardic et al, “Continents, Super-Continents, Mantle Thermal Mixing, and Mantle Thermal Isolation” in this session), we have begun to repeat this exhaustive and exhausting 2D study in 3D. Cooper, C.M., A. Lenardic, and L.-N. Moresi "Effects of continental insulation and the partioning of heat producing elements on the Earth's heat loss." Geophys. Res. Lett., 33 ,10.1029, 2006. Lenardic, A., L.-N. Moresi, A.M. Jellinek, and M. Manga "Continental insulation, mantle cooling, and the surface area of oceans and continents." Earth Planet. Sci. Lett., 234 ,317-333, 2005.
NASA Astrophysics Data System (ADS)
Leigh, S. J.; Purssell, C. P.; Billson, D. R.; Hutchins, D. A.
2014-09-01
Flow sensing is an essential technique required for a wide range of application environments ranging from liquid dispensing to utility monitoring. A number of different methodologies and deployment strategies have been devised to cover the diverse range of potential application areas. The ability to easily create new bespoke sensors for new applications is therefore of natural interest. Fused deposition modelling is a 3D printing technology based upon the fabrication of 3D structures in a layer-by-layer fashion using extruded strands of molten thermoplastic. The technology was developed in the late 1980s but has only recently come to more wide-scale attention outside of specialist applications and rapid prototyping due to the advent of low-cost 3D printing platforms such as the RepRap. Due to the relatively low-cost of the printers and feedstock materials, these printers are ideal candidates for wide-scale installation as localized manufacturing platforms to quickly produce replacement parts when components fail. One of the current limitations with the technology is the availability of functional printing materials to facilitate production of complex functional 3D objects and devices beyond mere concept prototypes. This paper presents the formulation of a simple magnetite nanoparticle-loaded thermoplastic composite and its incorporation into a 3D printed flow-sensor in order to mimic the function of a commercially available flow-sensing device. Using the multi-material printing capability of the 3D printer allows a much smaller amount of functional material to be used in comparison to the commercial flow sensor by only placing the material where it is specifically required. Analysis of the printed sensor also revealed a much more linear response to increasing flow rate of water showing that 3D printed devices have the potential to at least perform as well as a conventionally produced sensor.
Effects of flow control over a 3D turret -- Part II
NASA Astrophysics Data System (ADS)
Wallace, Ryan; Andino, Marlyn; Schmit, Ryan; Camphouse, Chris; Myatt, James; Glauser, Mark
2007-11-01
Building upon the 3D turret work done at Syracuse University an extended study was conducted in the Air Force Research Laboratory's Subsonic Aerodynamic Research Laboratory (SARL) wind tunnel at Wright-Patterson Air Force Base. The SARL experiments were performed at higher Reynolds and Mach numbers and therefore present a more complex, more challenging flow. Synthetic jets mounted upstream of the aperture were used to generate multiple actuation cases in order to provide a rich ensemble for plant model development based on the split POD method of Camphouse (2007). PIV velocity data was acquired along with simultaneous surface pressure data at various planes across the turret with and without open-loop control. In addition, a simple proportional closed-loop control was performed using the bandpass filtered first POD mode coefficient of the surface pressure as the feedback signal. The amplitude of the feedback signal was calibrated using the open-loop results which were the most effective in reducing the separation zone of the turret.
A Quasi-3-D Theory for Impedance Eduction in Uniform Grazing Flow
NASA Technical Reports Server (NTRS)
Watson, W. R.; Jones, M. G.; Parrott, T. L.
2005-01-01
A 2-D impedance eduction methodology is extended to quasi-3-D sound fields in uniform or shearing mean flow. We introduce a nonlocal, nonreflecting boundary condition to terminate the duct and then educe the impedance by minimizing an objective function. The introduction of a parallel, sparse, equation solver significantly reduces the wall clock time for educing the impedance when compared to that of the sequential band solver used in the 2-D methodology. The accuracy, efficiency, and robustness of the methodology is demonstrated using two examples. In the first example, we show that the method reproduces the known impedance of a ceramic tubular test liner. In the second example, we illustrate that the approach educes the impedance of a four-segment liner where the first, second, and fourth segments consist of a perforated face sheet bonded to honeycomb, and the third segment is a cut from the ceramic tubular test liner. The ability of the method to educe the impedances of multisegmented liners has the potential to significantly reduce the amount of time and cost required to determine the impedance of several uniform liners by allowing them to be placed in series in the test section and to educe the impedance of each segment using a single numerical experiment. Finally, we probe the objective function in great detail and show that it contains a single minimum. Thus, our objective function is ideal for use with local, inexpensive, gradient-based optimizers.
A 3-D Vortex Code for Parachute Flow Predictions: VIPAR Version 1.0
STRICKLAND, JAMES H.; HOMICZ, GREGORY F.; PORTER, VICKI L.; GOSSLER, ALBERT A.
2002-07-01
This report describes a 3-D fluid mechanics code for predicting flow past bluff bodies whose surfaces can be assumed to be made up of shell elements that are simply connected. Version 1.0 of the VIPAR code (Vortex Inflation PARachute code) is described herein. This version contains several first order algorithms that we are in the process of replacing with higher order ones. These enhancements will appear in the next version of VIPAR. The present code contains a motion generator that can be used to produce a large class of rigid body motions. The present code has also been fully coupled to a structural dynamics code in which the geometry undergoes large time dependent deformations. Initial surface geometry is generated from triangular shell elements using a code such as Patran and is written into an ExodusII database file for subsequent input into VIPAR. Surface and wake variable information is output into two ExodusII files that can be post processed and viewed using software such as EnSight{trademark}.
Thermocapillary bubble flow and coalescence in a rotating cylinder: A 3D study
NASA Astrophysics Data System (ADS)
Alhendal, Yousuf; Turan, A.; Al-mazidi, M.
2015-12-01
The process of thermocapillary bubbles rising in a rotating 3D cylinder in zero gravity was analysed and presented numerically with the aid of computational fluid dynamics (CFD) by means of the volume of fluid (VOF) method. Calculations were carried out to investigate in detail the effect of the rotational speed of the hosted liquid on the trajectory of both single and group bubbles driven by the Marangoni force in zero-gravity conditions. For rotational speeds from 0.25 to 2 rad/s, bubble displacement with angular motion was found to be directed between the hotter surface and the rotational axis. This is contrary to the conventional bubble flow from areas of high pressure to low pressure, radial direction, or from cold to hot regions, axial direction. The results demonstrate that for the ratio of rotational speeds to the thermocapillary bubble velocity larger than unity, the surface tension gradient is the dominant force and the bubble motion towards the hotter. On the other hand, for ratio less than 1, the bubble motion is dominated and is significantly affected by centrifugal force. As rotation speed increases, the amount of deflection increases and the Marangoni effect vanishes. The current study is novel in the sense that single- and multi-bubble motion incorporating thermocapillary forces in a rotating liquid in a zero-gravity environment has never been numerically investigated.
Permanganate oxidation of DNAPL in a large 3-D flow tank
NASA Astrophysics Data System (ADS)
Lee, E.; Seol, Y.; Fang, Y. C.; Schwartz, F. W.
2002-05-01
Potassium permanganate (KMnO4), as a metal-oxo reagent, can attack a double carbon-carbon bond and therefore oxidize common chlorinated ethylenes, such as trichloroethylene (TCE) and tetrachloroethylene (PCE). This feature of metal-oxo reagents facilitates the use of permanganate to remediation of chlorinated solvents in soil and groundwater. In this study, we evaluated the efficiency of TCE removal by permanganate oxidation in large three-dimensional flooding schemes. We constructed a large 3-D flow tank (L x W x D = 180 cm x 60 cm x 90 cm) where TCE source zone was installed in a saturated porous sandy medium. The tank was flushed at a flow rate of 51 L/day with permanganate solution (1,250 mg/L) for 63 days. Using chemical, electrical, and optical monitoring techniques we estimated temporal and areal variations in TCE, permanganate, MnO2 precipitates, conductivity, and chloride concentrations. TCE emplaced as DNAPL in a upstream source zone gradually moved downstream forming a TCE plume of about 120 cm long, 30 cm wide, and 55 cm deep. This TCE plume diminished considerably over time due to the in situ oxidation of the DNAPL. However, TCE was not completely destroyed and TCE concentration remained high (63 to 228 mg/L) in the shrunken TCE plume downstream after 63 days of permanganate flushing. Mass balance calculation indicated about 28% of TCE still remained in the system. This was attributed to the precipitation of low-permeability reaction by-product, i.e., MnO2, which caused flushing to become less efficient with time. Findings of this study are useful for developing a practical technique for enhancing the efficacy of the oxidative treatment of TCE using permanganate in the field conditions.
Quantitative OCT angiography of optic nerve head blood flow
Jia, Yali; Morrison, John C.; Tokayer, Jason; Tan, Ou; Lombardi, Lorinna; Baumann, Bernhard; Lu, Chen D.; Choi, WooJhon; Fujimoto, James G.; Huang, David
2012-01-01
Optic nerve head (ONH) blood flow may be associated with glaucoma development. A reliable method to quantify ONH blood flow could provide insight into the vascular component of glaucoma pathophysiology. Using ultrahigh-speed optical coherence tomography (OCT), we developed a new 3D angiography algorithm called split-spectrum amplitude-decorrelation angiography (SSADA) for imaging ONH microcirculation. In this study, a method to quantify SSADA results was developed and used to detect ONH perfusion changes in early glaucoma. En face maximum projection was used to obtain 2D disc angiograms, from which the average decorrelation values (flow index) and the percentage area occupied by vessels (vessel density) were computed from the optic disc and a selected region within it. Preperimetric glaucoma patients had significant reductions of ONH perfusion compared to normals. This pilot study indicates OCT angiography can detect the abnormalities of ONH perfusion and has the potential to reveal the ONH blood flow mechanism related to glaucoma. PMID:23243564
NASA Astrophysics Data System (ADS)
Izham, Mohamad Yusoff; Muhamad Uznir, Ujang; Alias, Abdul Rahman; Ayob, Katimon; Wan Ruslan, Ismail
2011-04-01
Existing 2D data structures are often insufficient for analysing the dynamism of saturation excess overland flow (SEOF) within a basin. Moreover, all stream networks and soil surface structures in GIS must be preserved within appropriate projection plane fitting techniques known as georeferencing. Inclusion of 3D volumetric structure of the current soft geo-objects simulation model would offer a substantial effort towards representing 3D soft geo-objects of SEOF dynamically within a basin by visualising saturated flow and overland flow volume. This research attempts to visualise the influence of a georeference system towards the dynamism of overland flow coverage and total overland flow volume generated from the SEOF process using VSG data structure. The data structure is driven by Green-Ampt methods and the Topographic Wetness Index (TWI). VSGs are analysed by focusing on spatial object preservation techniques of the conformal-based Malaysian Rectified Skew Orthomorphic (MRSO) and the equidistant-based Cassini-Soldner projection plane under the existing geodetic Malaysian Revised Triangulation 1948 (MRT48) and the newly implemented Geocentric Datum for Malaysia (GDM2000) datum. The simulated result visualises deformation of SEOF coverage under different georeference systems via its projection planes, which delineate dissimilar computation of SEOF areas and overland flow volumes. The integration of Georeference, 3D GIS and the saturation excess mechanism provides unifying evidence towards successful landslide and flood disaster management through envisioning the streamflow generating process (mainly SEOF) in a 3D environment.
Analysis of the TREAT loss-of-flow tests L6 and L7 using SAS3D
Morris, E.E.; Simms, R.; Gruber, E.E.
1985-01-01
The TREAT loss-of-flow tests L6 and L7 have been analyzed using the SAS3D accident analysis code. The impetus for the analysis was the need for experimentally supported fuel motion modeling in whole core accident studies performed in support of licensing of the Clinch River Breeder Reactor Project. The input prescription chosen for the SAS3D/SLUMPY fuel motion model gave reasonable agreement with the test results. Tests L6 and L7, each conducted with a cluster of three fuel pins, were planned to simulate key events in the loss-of-flow accident scenario for the Clinch River homogeneous reactor.
NASA Astrophysics Data System (ADS)
Wada, Ikuko; He, Jiangheng; Hasegawa, Akira; Nakajima, Junichi
2015-09-01
We develop a 3-D thermal model for the Northeast Japan subduction margin, using a realistic slab geometry for the subducting Pacific plate, and investigate the effects of oblique subduction and 3-D slab geometry on the mantle wedge flow pattern and the thermal structure. In the Tohoku region, the mantle wedge flow pattern is nearly two-dimensional resulting in a thermal structure similar to those obtained by a 2-D model, owing to the simple slab geometry and subduction nearly perpendicular to the margin. However, in Hokkaido, oblique subduction leads to 3-D mantle wedge flow with northerly inflow and west-northwestward outflow and also results in lower temperatures in the shallow part of the mantle wedge than in Tohoku due to lower sinking rate of the slab. Between Hokkaido and Tohoku, the slab has a hinge-like shape due to a relatively sharp change in the dip direction. In this hinge zone, northerly mantle inflow from Hokkaido and westerly mantle inflow from Tohoku converge, discouraging inflow from northwest and resulting in a cooler mantle wedge. The model-predicted mantle wedge flow patterns are consistent with observed seismic anisotropy and may explain the orientations of volcanic cross-chains. The predicted 3-D thermal structure correlates well with the along-arc variations in the location of the frontal arc volcanoes and help to provide new insights into the surface heat flow pattern and the down-dip extent of interplate earthquakes.
NASA Astrophysics Data System (ADS)
Jang, Jaeseong; Ahn, Chi Young; Jeon, Kiwan; Choi, Jung-il; Lee, Changhoon; Seo, Jin Keun
2015-03-01
A reconstruction method is proposed here to quantify the distribution of blood flow velocity fields inside the left ventricle from color Doppler echocardiography measurement. From 3D incompressible Navier- Stokes equation, a 2D incompressible Navier-Stokes equation with a mass source term is derived to utilize the measurable color flow ultrasound data in a plane along with the moving boundary condition. The proposed model reflects out-of-plane blood flows on the imaging plane through the mass source term. For demonstrating a feasibility of the proposed method, we have performed numerical simulations of the forward problem and numerical analysis of the reconstruction method. First, we construct a 3D moving LV region having a specific stroke volume. To obtain synthetic intra-ventricular flows, we performed a numerical simulation of the forward problem of Navier-Stokes equation inside the 3D moving LV, computed 3D intra-ventricular velocity fields as a solution of the forward problem, projected the 3D velocity fields on the imaging plane and took the inner product of the 2D velocity fields on the imaging plane and scanline directional velocity fields for synthetic scanline directional projected velocity at each position. The proposed method utilized the 2D synthetic projected velocity data for reconstructing LV blood flow. By computing the difference between synthetic flow and reconstructed flow fields, we obtained the averaged point-wise errors of 0.06 m/s and 0.02 m/s for u- and v-components, respectively.
NASA Astrophysics Data System (ADS)
Shaked, Natan Tzvi
2016-03-01
I present multimodal wide-field interferometric microscopy platform for label-free 3-D imaging of live cells during fast flow. Using holographic optical tweezers, multiple cells can be optically trapped and rapidity rotated on all axes, while acquired using an external off-axis wide-field interferometric module developed in our lab. The interferometric projections are rapidly processed into the 3-D refractive-index profile of the cells using a tomographic phase microscopy algorithms that take into consideration optical diffraction effects. The algorithms for the 3-D refractive-index reconstruction, and for calculating various morphological parameters that should serve for online sorting of cells, are efficiently implemented in a nearly real-time manner. The potential of this new high-throughput imaging technique is for label-free image analysis and sorting of cells during flow, to substitute current cell sorting devices, which are based on external labeling that eventually damages the cell sample.