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
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.
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.
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
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
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.
Ashby, S.F.; Falgout, R.D.; Smith, S.G.; Fogwell, T.W.
1994-09-01
This paper discusses the numerical simulation of groundwater flow through heterogeneous porous media. The focus is on the performance of a parallel multigrid preconditioner for accelerating convergence of conjugate gradients, which is used to compute the hydraulic pressure head. The numerical investigation considers the effects of enlarging the domain, increasing the grid resolution, and varying the geostatistical parameters used to define the subsurface realization. The results were obtained using the PARFLOW groundwater flow simulator on the Cray T3D massively parallel computer.
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
Barlebo, H.C.; Rosbjerg, D.; Hill, M.C.
1996-01-01
An extensive amount of data including hydraulic heads, hydraulic conductivities and concentrations of several solutes from controlled injections have been collected during the MADE 1 and MADE 2 experiments at a heterogeneous site near Columbus, Mississippi. In this paper the use of three-dimensional inverse groundwater models including simultaneous estimation of flow and transport parameters is proposed to help identify the dominant characteristics at the site. Simulations show that using a hydraulic conductivity distribution obtained from 2187 borehole flowmeter tests directly in the model produces poor matches to the measured hydraulic heads and tritium concentrations. Alternatively, time averaged hydraulic head maps are used to define zones of constant hydraulic conductivity to be estimated. Preliminary simulations suggest that in the case of conservative transport many, but not all, of the major plume characteristics can be explained by large-scale heterogeneity in recharge and hydraulic conductivity.
Quasi 3D modeling of water flow in vadose zone and groundwater
Technology Transfer Automated Retrieval System (TEKTRAN)
The 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 ...
SABAE-HW3D: a Meteor-Hydrological Model Coupling the Land Surface to Groundwater Flow
NASA Astrophysics Data System (ADS)
Loukili, Y.; Woodbury, A. D.
2007-12-01
A new coupled model linking the land surface scheme SABAE-HW with the saturated groundwater flow is introduced. SABAE-HW stands for: Soil Atmosphere Boundary, Accurate Evaluations of Heat and Water. It was presented (Loukili et al., 2007) as a parallel extension to the Canadian LAnd Surface Scheme (CLASS) which permits a free choice of soil column depth and layers. This development is important as it allows for the coupling of the atmosphere, land surface and subsurface zones. The generic domain under study is subdivided into soil columns whose surface areas represent GCM grid squares. The depths of columns extend to the saturated zone where a water table lower boundary condition is prescribed. In each column the unsaturated flow is forced by the corresponding meteorological data, and the moisture (liquid and frozen) and temperature profiles are computed by SABAE-HW using the half hour time step. The horizontal flow in the saturated zone is described by the vertically integrated model incorporating a storage term, and discretized using a finite volume scheme operating on the same GCM quadrilateral mesh. Since the saturated flow model affords larger time steps, bottom drainages from soil columns are summed up and input as cell recharges. Meanwhile, as the water table fluctuates, the individual column's mesh is allowed to deform. This physically based coupling strategy was selected for its superiority in providing stable and consistent results. In fact, the convergence to steady state situations in both the unsaturated and saturated zones is realized in many validation numerical tests. Moreover, our SABAE-HW3D code was successfully benchmarked against the finite element code Seep/W (Geo-Slope International, 2002) for steady and transient groundwater flows through different soil types. Even when handling coarse grids, SABAE-HW3D solutions are free of moisture oscillation and under or overshoot near the water table. The model is used to understand and assess properly
NASA Astrophysics Data System (ADS)
Gräbe, A. C.; Guttman, J.; Rödiger, T.; Siebert, C.; Merz, R.; Kolditz, O.
2012-12-01
Semi-arid to arid regions are usually characterized by a scarcity of precipitation and a lack of stream flow. Especially in desert environments, groundwater is one of the most important fresh water sources and its recharge is basically controlled by two main mechanisms: the direct regional infiltration of precipitation in the mountains and interdrainage areas in the first place and secondly the flood water infiltration through ephemeral channel beds (transmission loss). Due to extensive spatio-temporal data scarcity, direct quantitative estimations of groundwater recharge are often difficult to perform, and numerical models simulating the water fluxes, have to be applied to enable a quantitative approximation of the groundwater recharge. We made an assumption about the quantity of recharge for the subsurface catchment of the western Dead Sea escarpment, which is at the same time the input for the complex groundwater flow model of the Judea Group Aquifer. This can only be suggested if the hydrogeological situation in the tectonically complex region is fully understood. A number of simplified models of the Judea Group aquifer have been formulated and employed using a two-dimensional (one horizontal layered) numerical simulation of groundwater flow (Baida et al. 1978; Goldschtoff & Shachnai, 1980; Guttman, 2000; Laronne Ben-Itzhak & Gvirtzmann, 2005). However, all previous approaches focused only on a limited area of the Judea Group aquifer. We developed a high resolution regional groundwater flow model for the entire western basin of the Dead Sea. Whereas the structural model could be defined using a large geological dataset, the challenge was to generate the groundwater flow model with only limited well data. With the help of the scientific software OpenGeoSys (OGS) the challenge was reliably solved resulting in a simulation of the hydraulic characteristics (hydraulic conductivity and hydraulic head) of the cretaceous aquifer system, which was calibrated using PEST.
NASA Astrophysics Data System (ADS)
Zurek, A. J.; Witczak, S.; Dulinski, M.; Wachniew, P.; Rozanski, K.; Kania, J.; Postawa, A.; Karczewski, J.; Moscicki, W. J.
2014-08-01
A dedicated study was launched in 2010 with the main aim to better understand the functioning of groundwater dependent terrestrial ecosystem (GDTE) located in southern Poland. The GDTE consists of a valuable forest stand (Niepolomice Forest) and associated wetland (Wielkie Bloto fen). A wide range of tools (environmental tracers, geochemistry, geophysics, 3-D flow and transport modeling) was used. The research was conducted along three major directions: (i) quantification of the dynamics of groundwater flow in various parts of the aquifer associated with GDTE, (ii) quantification of the degree of interaction between the GDTE and the aquifer, and (iii) 3-D modeling of groundwater flow in the vicinity of the studied GDTE and quantification of possible impact of enhanced exploitation of the aquifer on the status of GDTE. Environmental tracer data (tritium, stable isotopes of water) strongly suggest that upward leakage of the aquifer contributes significantly to the present water balance of the studied wetland and associated forest. Physico-chemical parameters of water (pH, conductivity, Na / Cl ratio) confirm this notion. Model runs indicate that prolonged groundwater abstraction through the newly-established network of water supply wells, conducted at maximum permitted capacity (ca. 10 000 m3 d-1), may trigger drastic changes in the ecosystem functioning, eventually leading to its degradation.
NASA Astrophysics Data System (ADS)
Zurek, A. J.; Witczak, S.; Dulinski, M.; Wachniew, P.; Rozanski, K.; Kania, J.; Postawa, A.; Karczewski, J.; Moscicki, W. J.
2015-02-01
Groundwater-dependent ecosystems (GDEs) have important functions in all climatic zones as they contribute to biological and landscape diversity and provide important economic and social services. Steadily growing anthropogenic pressure on groundwater resources creates a conflict situation between nature and man which are competing for clean and safe sources of water. Such conflicts are particularly noticeable in GDEs located in densely populated regions. A dedicated study was launched in 2010 with the main aim to better understand the functioning of a groundwater-dependent terrestrial ecosystem (GDTE) located in southern Poland. The GDTE consists of a valuable forest stand (Niepolomice Forest) and associated wetland (Wielkie Błoto fen). It relies mostly on groundwater from the shallow Quaternary aquifer and possibly from the deeper Neogene (Bogucice Sands) aquifer. In July 2009 a cluster of new pumping wells abstracting water from the Neogene aquifer was set up 1 km to the northern border of the fen. A conceptual model of the Wielkie Błoto fen area for the natural, pre-exploitation state and for the envisaged future status resulting from intense abstraction of groundwater through the new well field was developed. The main aim of the reported study was to probe the validity of the conceptual model and to quantify the expected anthropogenic impact on the studied GDTE. A wide range of research tools was used. The results obtained through combined geologic, geophysical, geochemical, hydrometric and isotope investigations provide strong evidence for the existence of upward seepage of groundwater from the deeper Neogene aquifer to the shallow Quaternary aquifer supporting the studied GDTE. Simulations of the groundwater flow field in the study area with the aid of a 3-D flow and transport model developed for Bogucice Sands (Neogene) aquifer and calibrated using environmental tracer data and observations of hydraulic head in three different locations on the study area
NASA Astrophysics Data System (ADS)
Spottke, Ina; Zechner, Eric; Huggenberger, Peter
2005-09-01
A 3D geological model of the area east of Basel on the southeastern border of the Upper Rhine Graben, consisting of 47 faults and six stratigraphic horizons relevant for groundwater flow, was developed using borehole data, geological maps, geological cross sections, and outcrop data. This model provides new insight into the discussions about the kinematics of the area between the southeastern border of the Upper Rhine Graben and the Tabular Jura east of Basel. A 3D analysis showed that both thin-skinned and thick-skinned tectonic elements occur in the modeled area and that the Anticline and a series of narrow graben structures developed simultaneously during an extensional stress-field varying from E W to SSE NNW, which lasted from the Middle Eocene to Late Oligocene. In a new approach the faults and horizons of the 3D geological model were transferred into discrete elements with distributed hydrogeological properties in order to simulate the 3D groundwater flow regime within the modeled aquifers. A three-layer approach with a horizontal regularly spaced grid combined with an irregular property distribution of transmissivity in depth permitted the piezometric head of the steady-state model to be automatically calibrated to corresponding measurements using more than 200 piezometers. Groundwater modeling results demonstrated that large-scale industrial pumping affected the groundwater flow field in the Upper Muschelkalk aquifer at distances of up to 2 km to the south. The results of this research will act as the basis for further model developments, including salt dissolution and solute transport in the area, and may ultimately help to provide predictions for widespread land subsidence risks.
3D Groundwater Flow Model in the Arid Region of Tafilalet Oasis System (South East of Morocco)
NASA Astrophysics Data System (ADS)
Bouaamlat, I.; Larabi, A.; Faouzi, M.
2013-05-01
The plain of Tafilalet contains an important oasis located in the Southeast of Morocco in a pre-Saharan area, characterized by an arid climate with a large deficit water budget. It has a behavior of a large depression resulting from erosion of a set of geological coverage during the Quaternary period. It also forms a small Mesopotamia crossed by two main rivers from the mountains of the High Atlas: Ziz and Rheris. The oasis of Tafilalet is an area of old traditions irrigation where agriculture is the main activity of the region that represents approximately 37% of the total area (637 km2). In this study, a three-dimensional model of groundwater flow was developed for the aquifer system of Tafilalet, to assist the decision makers as a "management tool" in order to assess alternative schemes for development and exploitation of groundwater resources in the Tafilalet plain, using Modflow2000 code. It is the first mathematical model performed for this oasis plain, taking into account the most possible real hydrogeological conditions and using the geographical information system (GIS) for the organisation and treatment of data and applying a multidisciplinary approach combining geostatistical and hydrogeological modeling. The conceptual model, in terms of hydrogeological modeling was therefore considered as a monolayer model and the aquifer system is mainly heterogeneous with lateral different hydraulic conductivities, which are ranging from 3.10-7 to 5.10-2 m/s, but most of them are located between 2.10-4 and 8.10-3 m/s. The results of the model calibration under steady state (1960) and transient state conditions, starting from this time, show reasonable agreement between observed and simulated water levels for the observation wells. After calibration, the model contributed to better groundwater characterization, the hydrodynamic parameters obtained from the model are much representative of reality. As a management tool, this model can help the manager to take
NASA Astrophysics Data System (ADS)
Zurek, Anna; Witczak, Stanislaw; Kania, Jaroslaw; Wachniew, Przemyslaw; Rozanski, Kazimierz; Dulinski, Marek; Jench, Olga
2013-04-01
Niepolomice Forest. There is a growing concern that continued exploitation of those wells may lead to lowering water table in the Niepolomice Forest area and, as a consequence, may trigger drastic changes in this unique ecosystem. A dedicated study was launched with the main aim to quantify the interaction between Niepolomice Forest, with the focus the Wielkie Bloto fen, and the underlying Bogucice Sands aquifer. The work was pursued along three major lines: (i) vertical profiling of the Wielkie Bloto fen aimed at characterizing chemical and isotope contrast in the shallow groundwater occupying the Quaternary cover in order to identify upward leakage of deeper groundwater in the investigated area, (ii) regular monitoring of flow rate, chemistry and environmental isotopes of the Dluga Woda stream draining the Wielkie Bloto fen, and (iii) 3D modeling of groundwater flow in the vicinity of the Wielkie Bloto fen focusing on quantifying the impact of the Wola Batorska well field on the regional groundwater flow patterns. The results of isotope and chemical analyses confirmed existence of upward seepage of groundwater from the Bogucice Sands aquifer in the area of Wielkie Bloto fen. Preliminary assessment of the water balance of Dluga Woda catchment indicates that the baseflow originating from groundwater seepage is equal approximately 16% of the annual precipitation. Results of 3D flow model applied to the study area indicate that prolonged operation of the well-field Wola Batorska at maximum capacity may lead to substantial lowering of water table in the Niepolomice Forest area and, as a consequence, endanger further existence of this unique GDTE. Acknowledgements. Partial financial support of this work through GENESIS project (http:/www.thegenesisproject.eu) funded by the European Commission 7FP contract 226536, and through statutory funds of the AGH University of Science and Technology (projects No.11.11.140.026 and 11.11.220.01) is kindly acknowledged.
Voss, Clifford I.; Provost, A.M.
2002-01-01
SUTRA (Saturated-Unsaturated Transport) is a computer program that simulates fluid movement and the transport of either energy or dissolved substances in a subsurface environment. This upgraded version of SUTRA adds the capability for three-dimensional simulation to the former code (Voss, 1984), which allowed only two-dimensional simulation. The code employs a two- or three-dimensional finite-element and finite-difference method to approximate the governing equations that describe the two interdependent processes that are simulated: 1) fluid density-dependent saturated or unsaturated ground-water flow; and 2) either (a) transport of a solute in the ground water, in which the solute may be subject to: equilibrium adsorption on the porous matrix, and both first-order and zero-order production or decay; or (b) transport of thermal energy in the ground water and solid matrix of the aquifer. SUTRA may also be used to simulate simpler subsets of the above processes. A flow-direction-dependent dispersion process for anisotropic media is also provided by the code and is introduced in this report. As the primary calculated result, SUTRA provides fluid pressures and either solute concentrations or temperatures, as they vary with time, everywhere in the simulated subsurface system. SUTRA flow simulation may be employed for two-dimensional (2D) areal, cross sectional and three-dimensional (3D) modeling of saturated ground-water flow systems, and for cross sectional and 3D modeling of unsaturated zone flow. Solute-transport simulation using SUTRA may be employed to model natural or man-induced chemical-species transport including processes of solute sorption, production, and decay. For example, it may be applied to analyze ground-water contaminant transport problems and aquifer restoration designs. In addition, solute-transport simulation with SUTRA may be used for modeling of variable-density leachate movement, and for cross sectional modeling of saltwater intrusion in
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.
Kipp, K.L.
1987-01-01
The Heat- and Soil-Transport Program (HST3D) simulates groundwater flow and associated heat and solute transport in three dimensions. The three governing equations are coupled through the interstitial pore velocity, the dependence of the fluid density on pressure, temperature, the solute-mass fraction , and the dependence of the fluid viscosity on temperature and solute-mass fraction. The solute transport equation is for only a single, solute species with possible linear equilibrium sorption and linear decay. Finite difference techniques are used to discretize the governing equations using a point-distributed grid. The flow-, heat- and solute-transport equations are solved , in turn, after a particle Gauss-reduction scheme is used to modify them. The modified equations are more tightly coupled and have better stability for the numerical solutions. The basic source-sink term represents wells. A complex well flow model may be used to simulate specified flow rate and pressure conditions at the land surface or within the aquifer, with or without pressure and flow rate constraints. Boundary condition types offered include specified value, specified flux, leakage, heat conduction, and approximate free surface, and two types of aquifer influence functions. All boundary conditions can be functions of time. Two techniques are available for solution of the finite difference matrix equations. One technique is a direct-elimination solver, using equations reordered by alternating diagonal planes. The other technique is an iterative solver, using two-line successive over-relaxation. A restart option is available for storing intermediate results and restarting the simulation at an intermediate time with modified boundary conditions. This feature also can be used as protection against computer system failure. Data input and output may be in metric (SI) units or inch-pound units. Output may include tables of dependent variables and parameters, zoned-contour maps, and plots of the
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
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.
NASA Astrophysics Data System (ADS)
Kupila, Juho
2015-04-01
3D modeling of groundwater areas is an important research method in groundwater surveys. Model of geological soil structure improves the knowledge of linkage between land use planning and groundwater protection. Results can be used as base information when developing the water supply services and anticipating and performing the measures needed in case of environmental accidents. Also, collected information is utilized when creating the groundwater flow model. In Finland, structure studies have been conducted in cooperation (among others) with the municipalities and local water suppliers and with the authorities from the Centre for Economic Development, Transport and the Environment. Geological Survey of Finland carries out project "Structure studies in Kolpene groundwater area" in Rovaniemi, Finnish Lapland. Study site is located in northern Finland, in the vicinity of the city center of Rovaniemi. Extent of the area is about 13 square kilometers and there are lots of urban residential areas and other human activities. The objective of this project is to determine the geological structure of the Kolpene groundwater area so that the results can be used to estimate the validity of the present exclusion area and possible risks to the groundwater caused by the land use. Soil layers of the groundwater area are studied by means of collecting information by heavy drilling, geophysical surveying (ground penetrating radar and gravimeter measurements) and water sampling from the installed observation pipes. Also the general geological and hydrological mappings are carried out. Main results which will be produced are: 1) the model of the bedrock surface, 2) the model of the surface of the ground water and flow directions, 3) the thickness of ground water saturated soil layers and 4) location and main characteristics of the soil layers which are significant to the ground water conditions. The preparing studies have been started at the end of 2013 and the results will be
NASA Astrophysics Data System (ADS)
San Nicolo, Lorenz; Burger, Ulrich; Zurlo, Raffaele
2014-05-01
The Brenner Base Tunnel (BBT) will cross the Isarco Valley near the village of Fortezza (BZ) at a depth of approximately 20 m below the riverbed of the Isarco river. The design of this roughly 1 km long stretch through alluvial sediments and below groundwater level required detailed knowledge of the prevailing hydrogeological conditions. In particular, it was necessary to determine if dewatering procedures were feasible and what the impacts on natural water flows in the aquifer after completion of the infrastructure will be. The study area is a typical Alpine valley, filled with alluvial sediments to a maximum depth of approximately 120m. The valley is bounded by granitic rocks with regional, water saturated main fault zones. In addition to the Isarco River, the area is shaped by two lateral rivers. The deposits of these lateral rivers form main alluvial fans. A 3D flow model of the aquifer was elaborated as part of the framework of the study. The model was calibrated and validated on the basis of two subsequent analyses: - a preliminary hydrogeological characterization of the area by means of a pumping test carried out in a well. - an experimental pumping test in five wells with a significant drawdown of the aquifer level in a broad area of the study zone. In addition to providing basic data for the planning of eventual dewatering procedures and for the simulation of the impact of the completed infrastructure, the model highlighted the complex interaction between the Isarco River and the aquifer. Significant variations in the rate of infiltration of water from the river into the groundwater were ascertained that are caused by changes in the permeability of the riverbed over very short distances. The Isarco River, which in the simulation domain extends over 1.5 km, was divided into nine segments with variations in riverbed permeability of more than an order of magnitude. The causes of these variations were not analysed in detail during this study. However, this
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.
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.
Brien, Dianne L.; Reid, Mark E.
2008-01-01
In Seattle, Washington, deep-seated landslides on bluffs along Puget Sound have historically caused extensive damage to land and structures. These large failures are controlled by three-dimensional (3-D) variations in strength and pore-water pressures. We assess the slope stability of part of southwestern Seattle using a 3-D limit-equilibrium analysis coupled with a 3-D groundwater flow model. Our analyses use a high-resolution digital elevation model (DEM) combined with assignment of strength and hydraulic properties based on geologic units. The hydrogeology of the Seattle area consists of a layer of permeable glacial outwash sand that overlies less permeable glacial lacustrine silty clay. Using a 3-D groundwater model, MODFLOW-2000, we simulate a water table above the less permeable units and calibrate the model to observed conditions. The simulated pore-pressure distribution is then used in a 3-D slope-stability analysis, SCOOPS, to quantify the stability of the coastal bluffs. For wet winter conditions, our analyses predict that the least stable areas are steep hillslopes above Puget Sound, where pore pressures are elevated in the outwash sand. Groundwater flow converges in coastal reentrants, resulting in elevated pore pressures and destabilization of slopes. Regions predicted to be least stable include the areas in or adjacent to three mapped historically active deep-seated landslides. The results of our 3-D analyses differ significantly from a slope map or results from one-dimensional (1-D) analyses.
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.
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.
NASA Astrophysics Data System (ADS)
Daughney, C.; Toews, M. W.; Morgenstern, U.; Cornaton, F. J.; Jackson, B. M.
2013-12-01
Lake Rotorua is a focus of culture and tourism in New Zealand. The lake's water quality has declined since the 1970s, partly due to nutrient inputs that reach the lake via the groundwater system. Improved land use management within the catchment requires prediction of the spatial variations of groundwater transit time from land surface to the lake, and from this the prediction of current and future nutrient inflows to the lake. This study combines the two main methods currently available for determination of water age: numerical groundwater models and hydrological tracers. A steady-state 3D finite element model was constructed to simulate groundwater flow and transport of tritium and age at the catchment scale (555 km2). The model materials were defined using a 3D geologic model and included ignimbrites, rhyolites, alluvial and lake bottom sediments. The steady-state saturated groundwater flow model was calibrated using observed groundwater levels in boreholes (111 locations) and stream flow measurements from groundwater-fed streams and springs (61 locations). Hydraulic conductivities and Cauchy boundary conditions associated with the streams, springs and lake were parameterized. The transport parameters for the model were calibrated using 191 tritium samples from 105 locations (springs, streams and boreholes), with most locations having two sample dates. The transport model used steady-state flow, but simulated the transient transport and decay of tritium from rainfall recharge between 1945 and 2012. An additional 1D unsaturated sub-model was added to account for tritium decay from the ground surface to the water table. The sub-model is linked on top of the 3D model, and uses the water table depths and material properties from the 3D model. The adjustable calibration parameters for the transport model were porosity and van Genuchten parameters related to the unsaturated sub-models. Calibration of the flow model was achieved using a combination of automated least
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.
Nicole Lautze
2015-01-01
Groundwater flow model for the island of Oahu. Data is from the following sources: Rotzoll, K., A.I. El-Kadi. 2007. Numerical Ground-Water Flow Simulation for Red Hill Fuel Storage Facilities, NAVFAC Pacific, Oahu, Hawaii - Prepared TEC, Inc. Water Resources Research Center, University of Hawaii, Honolulu.; Whittier, R.B., K. Rotzoll, S. Dhal, A.I. El-Kadi, C. Ray, G. Chen, and D. Chang. 2004. Hawaii Source Water Assessment Program Report – Volume VII – Island of Oahu Source Water Assessment Program Report. Prepared for the Hawaii Department of Health, Safe Drinking Water Branch. University of Hawaii, Water Resources Research Center. Updated 2008.; and Whittier, R. and A.I. El-Kadi. 2009. Human and Environmental Risk Ranking of Onsite Sewage Disposal Systems – Final. Prepared by the University of Hawaii, Dept. of Geology and Geophysics for the State of Hawaii Dept. of Health, Safe Drinking Water Branch. December 2009.
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.
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
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.
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.
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
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.
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.
Groundwater hydrology--coastal flow
Sanford, Ward E.
2010-01-01
How groundwater flow varies when long-term external conditions change is little documented. Geochemical evidence shows that sea-level rise at the end of the last glacial period led to a shift in the flow patterns of coastal groundwater beneath Florida.
Nicole Lautze
2015-01-01
Groundwater flow model for Kauai. Data is from the following sources: Whittier, R. and A.I. El-Kadi. 2014. Human and Environmental Risk Ranking of Onsite Sewage Disposal Systems For the Hawaiian Islands of Kauai, Molokai, Maui, and Hawaii – Final. Prepared by the University of Hawaii, Dept. of Geology and Geophysics for the State of Hawaii Dept. of Health, Safe Drinking Water Branch. September 2014.; and Whittier, R.B., K. Rotzoll, S. Dhal, A.I. El-Kadi, C. Ray, G. Chen, and D. Chang. 2004. Hawaii Source Water Assessment Program Report – Volume IV – Island of Kauai Source Water Assessment Program Report. Prepared for the Hawaii Department of Health, Safe Drinking Water Branch. University of Hawaii, Water Resources Research Center. Updated 2015.
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
NASA Astrophysics Data System (ADS)
Janssen, G.; Del Val Alonso, L.; Groenendijk, P.; Griffioen, J.
2012-12-01
We developed an on-line coupling between the 1D/quasi-2D nutrient transport model ANIMO and the 3D groundwater transport model code MT3DMS. ANIMO is a detailed, process-oriented model code for the simulation of nitrate leaching to groundwater, N- and P-loads on surface waters and emissions of greenhouse gasses. It is the leading nutrient fate and transport code in the Netherlands where it is used primarily for the evaluation of fertilization related legislation. In addition, the code is applied frequently in international research projects. MT3DMS is probably the most commonly used groundwater solute transport package worldwide. The on-line model coupling ANIMO-MT3DMS combines the state-of-the-art descriptions of the biogeochemical cycles in ANIMO with the advantages of using a 3D approach for the transport through the saturated domain. These advantages include accounting for regional lateral transport, considering groundwater-surface water interactions more explicitly, and the possibility of using MODFLOW to obtain the flow fields. An additional merit of the on-line coupling concept is that it preserves feedbacks between the saturated and unsaturated zone. We tested ANIMO-MT3DMS by simulating nutrient transport for the period 1970-2007 in a Dutch agricultural polder catchment covering an area of 118 km2. The transient groundwater flow field had a temporal resolution of one day and was calculated with MODFLOW-MetaSWAP. The horizontal resolution of the model grid was 100x100m and consisted of 25 layers of varying thickness. To keep computation times manageable, we prepared MT3DMS for parallel computing, which in itself is a relevant development for a large community of groundwater transport modelers. For the parameterization of the soil, we applied a standard classification approach, representing the area by 60 units with unique combinations of soil type, land use and geohydrological setting. For the geochemical parameterization of the deeper subsurface, however, we
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
Monitoring probe for groundwater flow
Looney, B.B.; Ballard, S.
1994-08-23
A monitoring probe for detecting groundwater migration is disclosed. The monitor features a cylinder made of a permeable membrane carrying an array of electrical conductivity sensors on its outer surface. The cylinder is filled with a fluid that has a conductivity different than the groundwater. The probe is placed in the ground at an area of interest to be monitored. The fluid, typically saltwater, diffuses through the permeable membrane into the groundwater. The flow of groundwater passing around the permeable membrane walls of the cylinder carries the conductive fluid in the same general direction and distorts the conductivity field measured by the sensors. The degree of distortion from top to bottom and around the probe is precisely related to the vertical and horizontal flow rates, respectively. The electrical conductivities measured by the sensors about the outer surface of the probe are analyzed to determine the rate and direction of the groundwater flow. 4 figs.
Monitoring probe for groundwater flow
Looney, Brian B.; Ballard, Sanford
1994-01-01
A monitoring probe for detecting groundwater migration. The monitor features a cylinder made of a permeable membrane carrying an array of electrical conductivity sensors on its outer surface. The cylinder is filled with a fluid that has a conductivity different than the groundwater. The probe is placed in the ground at an area of interest to be monitored. The fluid, typically saltwater, diffuses through the permeable membrane into the groundwater. The flow of groundwater passing around the permeable membrane walls of the cylinder carries the conductive fluid in the same general direction and distorts the conductivity field measured by the sensors. The degree of distortion from top to bottom and around the probe is precisely related to the vertical and horizontal flow rates, respectively. The electrical conductivities measured by the sensors about the outer surface of the probe are analyzed to determine the rate and direction of the groundwater flow.
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.
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
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
Development of a geological 3D-model for improved calculations of groundwater vulnerability
NASA Astrophysics Data System (ADS)
Lerch, C.; Hoppe, A.
2007-06-01
The role and function of the unsaturated zone in groundwater protection is a basic question in land use planning. To quantify this protective ability, different approaches have been adopted which focus on various issues including the water retention period, behavior of organic components or amount of infiltration water available. Most methods use information from drill logs and interpolate the point data to delineate a groundwater protection map. In some cases the dimensions of the geological structures are considered, but often the spatial structure is disregarded. In this paper, a method is introduced which combines the development of a 3D geological model and the application of a standard method for computing groundwater vulnerability. Thus a spatially-continuous evaluation of groundwater protection is possible. As a test case, a spatial geological model was set up at a site southeast of Frankfurt am Main, Germany. The result was used to compute and map the local groundwater vulnerability. The result was compared with interpolated maps from drill data.
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
NASA Astrophysics Data System (ADS)
Stirewalt, G. L.; Shepherd, J. C.
2003-12-01
/liter) around the MPB and elevated nitrate (> 2000 milligrams/ liter) around storage ponds. Vertical connectivity was suggested between the TGWS and SGWS, while the DGWS appeared relatively isolated from the overlying aquifers. Lateral movement of uranium was also suggested over time. For example, lateral migration in the TGWS is suggested along a shallow depression in the bedrock surface trending south-southwest from the southwest corner of the MPB. Another pathway atop the buried bedrock surface, trending west-northwest from the MPB and partially reflected by current surface topography, suggested lateral migration of nitrate in the SGWS. Lateral movement of nitrate in the SGWS was also indicated north, south, and west of the largest storage pond. Definition of contaminant plume movement over time is particularly important for assessing direction and rate of migration and the potential need for preventive measures to control contamination of groundwater outside facility property lines. The 3D geospatial property models proved invaluable for visualizing and analyzing variations in subsurface uranium and nitrate contamination in space and time within and between the three aquifers at the site. The models were an exceptional visualization tool for illustrating extent, volume, and quantitative amounts of uranium and nitrate contamination in the subsurface to regulatory decision-makers in regard to site decommissioning issues, including remediation concerns, providing a perspective not possible to achieve with traditional 2D maps. The geohydrologic framework model provides a conceptual model for consideration in flow and transport analyses.
NASA Astrophysics Data System (ADS)
Hassen, Imen; Gibson, Helen; Hamzaoui-Azaza, Fadoua; Negro, François; Rachid, Khanfir; Bouhlila, Rachida
2016-08-01
The challenge of this study was to create a 3D geological and structural model of the Kasserine Aquifer System (KAS) in central Tunisia and its natural extension into north-east Algeria. This was achieved using an implicit 3D method, which honors prior geological data for both formation boundaries and faults. A current model is presented which provides defendable predictions for the spatial distribution of geology and water resources in aquifers throughout the model-domain. This work has allowed validation of regional scale geology and fault networks in the KAS, and has facilitated the first-ever estimations of groundwater resources in this region by a 3D method. The model enables a preliminary assessment of the hydraulic significance of the major faults by evaluating their influence and role on groundwater flow within and between four compartments of the multi-layered, KAS hydrogeological system. Thus a representative hydrogeological model of the study area is constructed. The possible dual nature of faults in the KAS is discussed in the context that some faults appear to be acting both as barriers to horizontal groundwater flow, and simultaneously as conduits for vertical flow. Also discussed is the possibility that two flow directions occur within the KAS, at a small syncline area of near Feriana. In summary, this work evaluates the influence of aquifer connectivity and the role of faults and geology in groundwater flow within the KAS aquifer system. The current KAS geological model can now be used to guide groundwater managers on the best placement for drilling to test and further refine the understanding of the groundwater system, including the faults connectivity. As more geological data become available, the current model can be easily edited and re-computed to provide an updated model ready for the next stage of investigation by numerical flow modeling.
NASA Astrophysics Data System (ADS)
Li, S.
2002-05-01
Taking advantage of the recent developments in groundwater modeling research and computer, image and graphics processing, and objected oriented programming technologies, Dr. Li and his research group have recently developed a comprehensive software system for unified deterministic and stochastic groundwater modeling. Characterized by a new real-time modeling paradigm and improved computational algorithms, the software simulates 3D unsteady flow and reactive transport in general groundwater formations subject to both systematic and "randomly" varying stresses and geological and chemical heterogeneity. The software system has following distinct features and capabilities: Interactive simulation and real time visualization and animation of flow in response to deterministic as well as stochastic stresses. Interactive, visual, and real time particle tracking, random walk, and reactive plume modeling in both systematically and randomly fluctuating flow. Interactive statistical inference, scattered data interpolation, regression, and ordinary and universal Kriging, conditional and unconditional simulation. Real-time, visual and parallel conditional flow and transport simulations. Interactive water and contaminant mass balance analysis and visual and real-time flux update. Interactive, visual, and real time monitoring of head and flux hydrographs and concentration breakthroughs. Real-time modeling and visualization of aquifer transition from confined to unconfined to partially de-saturated or completely dry and rewetting Simultaneous and embedded subscale models, automatic and real-time regional to local data extraction; Multiple subscale flow and transport models Real-time modeling of steady and transient vertical flow patterns on multiple arbitrarily-shaped cross-sections and simultaneous visualization of aquifer stratigraphy, properties, hydrological features (rivers, lakes, wetlands, wells, drains, surface seeps), and dynamically adjusted surface flooding area
Bailey, Ryan T.; Morway, Eric D.; Niswonger, Richard G.; Gates, Timothy K.
2013-01-01
A numerical model was developed that is capable of simulating multispecies reactive solute transport in variably saturated porous media. This model consists of a modified version of the reactive transport model RT3D (Reactive Transport in 3 Dimensions) that is linked to the Unsaturated-Zone Flow (UZF1) package and MODFLOW. Referred to as UZF-RT3D, the model is tested against published analytical benchmarks as well as other published contaminant transport models, including HYDRUS-1D, VS2DT, and SUTRA, and the coupled flow and transport modeling system of CATHY and TRAN3D. Comparisons in one-dimensional, two-dimensional, and three-dimensional variably saturated systems are explored. While several test cases are included to verify the correct implementation of variably saturated transport in UZF-RT3D, other cases are included to demonstrate the usefulness of the code in terms of model run-time and handling the reaction kinetics of multiple interacting species in variably saturated subsurface systems. As UZF1 relies on a kinematic-wave approximation for unsaturated flow that neglects the diffusive terms in Richards equation, UZF-RT3D can be used for large-scale aquifer systems for which the UZF1 formulation is reasonable, that is, capillary-pressure gradients can be neglected and soil parameters can be treated as homogeneous. Decreased model run-time and the ability to include site-specific chemical species and chemical reactions make UZF-RT3D an attractive model for efficient simulation of multispecies reactive transport in variably saturated large-scale subsurface systems.
Bailey, Ryan T; Morway, Eric D; Niswonger, Richard G; Gates, Timothy K
2013-01-01
A numerical model was developed that is capable of simulating multispecies reactive solute transport in variably saturated porous media. This model consists of a modified version of the reactive transport model RT3D (Reactive Transport in 3 Dimensions) that is linked to the Unsaturated-Zone Flow (UZF1) package and MODFLOW. Referred to as UZF-RT3D, the model is tested against published analytical benchmarks as well as other published contaminant transport models, including HYDRUS-1D, VS2DT, and SUTRA, and the coupled flow and transport modeling system of CATHY and TRAN3D. Comparisons in one-dimensional, two-dimensional, and three-dimensional variably saturated systems are explored. While several test cases are included to verify the correct implementation of variably saturated transport in UZF-RT3D, other cases are included to demonstrate the usefulness of the code in terms of model run-time and handling the reaction kinetics of multiple interacting species in variably saturated subsurface systems. As UZF1 relies on a kinematic-wave approximation for unsaturated flow that neglects the diffusive terms in Richards equation, UZF-RT3D can be used for large-scale aquifer systems for which the UZF1 formulation is reasonable, that is, capillary-pressure gradients can be neglected and soil parameters can be treated as homogeneous. Decreased model run-time and the ability to include site-specific chemical species and chemical reactions make UZF-RT3D an attractive model for efficient simulation of multispecies reactive transport in variably saturated large-scale subsurface systems. PMID:23131109
Global scale groundwater flow model
NASA Astrophysics Data System (ADS)
Sutanudjaja, Edwin; de Graaf, Inge; van Beek, Ludovicus; Bierkens, Marc
2013-04-01
As the world's largest accessible source of freshwater, groundwater plays vital role in satisfying the basic needs of human society. It serves as a primary source of drinking water and supplies water for agricultural and industrial activities. During times of drought, groundwater sustains water flows in streams, rivers, lakes and wetlands, and thus supports ecosystem habitat and biodiversity, while its large natural storage provides a buffer against water shortages. Yet, the current generation of global scale hydrological models does not include a groundwater flow component that is a crucial part of the hydrological cycle and allows the simulation of groundwater head dynamics. In this study we present a steady-state MODFLOW (McDonald and Harbaugh, 1988) groundwater model on the global scale at 5 arc-minutes resolution. Aquifer schematization and properties of this groundwater model were developed from available global lithological model (e.g. Dürr et al., 2005; Gleeson et al., 2010; Hartmann and Moorsdorff, in press). We force the groundwtaer model with the output from the large-scale hydrological model PCR-GLOBWB (van Beek et al., 2011), specifically the long term net groundwater recharge and average surface water levels derived from routed channel discharge. We validated calculated groundwater heads and depths with available head observations, from different regions, including the North and South America and Western Europe. Our results show that it is feasible to build a relatively simple global scale groundwater model using existing information, and estimate water table depths within acceptable accuracy in many parts of the world.
: Belcher, Wayne R., (Edited By); Sweetkind, Donald S.
2010-01-01
A numerical three-dimensional (3D) transient groundwater flow model of the Death Valley region was developed by the U.S. Geological Survey for the U.S. Department of Energy programs at the Nevada Test Site and at Yucca Mountain, Nevada. Decades of study of aspects of the groundwater flow system and previous less extensive groundwater flow models were incorporated and reevaluated together with new data to provide greater detail for the complex, digital model. A 3D digital hydrogeologic framework model (HFM) was developed from digital elevation models, geologic maps, borehole information, geologic and hydrogeologic cross sections, and other 3D models to represent the geometry of the hydrogeologic units (HGUs). Structural features, such as faults and fractures, that affect groundwater flow also were added. The HFM represents Precambrian and Paleozoic crystalline and sedimentary rocks, Mesozoic sedimentary rocks, Mesozoic to Cenozoic intrusive rocks, Cenozoic volcanic tuffs and lavas, and late Cenozoic sedimentary deposits of the Death Valley regional groundwater flow system (DVRFS) region in 27 HGUs. Information from a series of investigations was compiled to conceptualize and quantify hydrologic components of the groundwater flow system within the DVRFS model domain and to provide hydraulic-property and head-observation data used in the calibration of the transient-flow model. These studies reevaluated natural groundwater discharge occurring through evapotranspiration (ET) and spring flow; the history of groundwater pumping from 1913 through 1998; groundwater recharge simulated as net infiltration; model boundary inflows and outflows based on regional hydraulic gradients and water budgets of surrounding areas; hydraulic conductivity and its relation to depth; and water levels appropriate for regional simulation of prepumped and pumped conditions within the DVRFS model domain. Simulation results appropriate for the regional extent and scale of the model were provided
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
Calmbacher, C.W. )
1992-12-01
Visualizing and delineating subsurface geological features, groundwater contaminant plumes, soil contamination, geological faults, shears and other features can prove invaluable to environmental consultants, engineers, geologists and hydrogeologists. Three-dimensional modeling is useful for a variety of applications from planning remediation to site planning design. The problem often is figuring out how to convert drilling logs, map lists or contaminant levels from soil and groundwater into a 3-D model. Three-dimensional subsurface modeling is not a new requirement, but a flexible, easily applied method of developing such models has not always been readily available. LYNX Geosystems Inc. has developed the Geoscience Modeling System (GMS) in answer to the needs of those regularly having to do three-dimensional geostatistical modeling. The GMS program has been designed to allow analysis, interpretation and visualization of complex geological features and soil and groundwater contamination. This is a powerful program driven by a 30 volume modeling technology engine. Data can be entered, stored, manipulated and analyzed in ways that will present very few limitations to the user. The program has selections for Geoscience Data Management, Geoscience Data Analysis, Geological Modeling (interpretation and analysis), Geostatistical Modeling and an optional engineering component.
: Belcher, Wayne R., (Edited By)
2004-01-01
A numerical three-dimensional (3D) transient ground-water flow model of the Death Valley region was developed by the U.S. Geological Survey for the U.S. Department of Energy programs at the Nevada Test Site and at Yucca Mountain, Nevada. Decades of study of aspects of the ground-water flow system and previous less extensive ground-water flow models were incorporated and reevaluated together with new data to provide greater detail for the complex, digital model. A 3D digital hydrogeologic framework model (HFM) was developed from digital elevation models, geologic maps, borehole information, geologic and hydrogeologic cross sections, and other 3D models to represent the geometry of the hydrogeologic units (HGUs). Structural features, such as faults and fractures, that affect ground-water flow also were added. The HFM represents Precambrian and Paleozoic crystalline and sedimentary rocks, Mesozoic sedimentary rocks, Mesozoic to Cenozoic intrusive rocks, Cenozoic volcanic tuffs and lavas, and late Cenozoic sedimentary deposits of the Death Valley Regional Ground-Water Flow System (DVRFS) region in 27 HGUs. Information from a series of investigations was compiled to conceptualize and quantify hydrologic components of the ground-water flow system within the DVRFS model domain and to provide hydraulic-property and head-observation data used in the calibration of the transient-flow model. These studies reevaluated natural ground-water discharge occurring through evapotranspiration and spring flow; the history of ground-water pumping from 1913 through 1998; ground-water recharge simulated as net infiltration; model boundary inflows and outflows based on regional hydraulic gradients and water budgets of surrounding areas; hydraulic conductivity and its relation to depth; and water levels appropriate for regional simulation of prepumped and pumped conditions within the DVRFS model domain. Simulation results appropriate for the regional extent and scale of the model were
Sweetkind, Donald S.; Faunt, Claudia C.; Hanson, Randall T.
2013-01-01
Groundwater is the sole source of water supply in Cuyama Valley, a rural agricultural area in Santa Barbara County, California, in the southeasternmost part of the Coast Ranges of California. Continued groundwater withdrawals and associated water-resource management concerns have prompted an evaluation of the hydrogeology and water availability for the Cuyama Valley groundwater basin by the U.S. Geological Survey, in cooperation with the Water Agency Division of the Santa Barbara County Department of Public Works. As a part of the overall groundwater evaluation, this report documents the construction of a digital three-dimensional geologic framework model of the groundwater basin suitable for use within a numerical hydrologic-flow model. The report also includes an analysis of the spatial variability of lithology and grain size, which forms the geologic basis for estimating aquifer hydraulic properties. The geologic framework was constructed as a digital representation of the interpreted geometry and thickness of the principal stratigraphic units within the Cuyama Valley groundwater basin, which include younger alluvium, older alluvium, and the Morales Formation, and underlying consolidated bedrock. The framework model was constructed by creating gridded surfaces representing the altitude of the top of each stratigraphic unit from various input data, including lithologic and electric logs from oil and gas wells and water wells, cross sections, and geologic maps. Sediment grain-size data were analyzed in both two and three dimensions to help define textural variations in the Cuyama Valley groundwater basin and identify areas with similar geologic materials that potentially have fairly uniform hydraulic properties. Sediment grain size was used to construct three-dimensional textural models that employed simple interpolation between drill holes and two-dimensional textural models for each stratigraphic unit that incorporated spatial structure of the textural data.
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
Age, double porosity, and simple reaction modifications for the MOC3D ground-water transport model
Goode, Daniel J.
1999-01-01
This report documents modifications for the MOC3D ground-water transport model to simulate (a) ground-water age transport; (b) double-porosity exchange; and (c) simple but flexible retardation, decay, and zero-order growth reactions. These modifications are incorporated in MOC3D version 3.0. MOC3D simulates the transport of a single solute using the method-ofcharacteristics numerical procedure. The age of ground water, that is the time since recharge to the saturated zone, can be simulated using the transport model with an additional source term of unit strength, corresponding to the rate of aging. The output concentrations of the model are in this case the ages at all locations in the model. Double porosity generally refers to a separate immobilewater phase within the aquifer that does not contribute to ground-water flow but can affect solute transport through diffusive exchange. The solute mass exchange rate between the flowing water in the aquifer and the immobile-water phase is the product of the concentration difference between the two phases and a linear exchange coefficient. Conceptually, double porosity can approximate the effects of dead-end pores in a granular porous media, or matrix diffusion in a fractured-rock aquifer. Options are provided for decay and zero-order growth reactions within the immobilewater phase. The simple reaction terms here extend the original model, which included decay and retardation. With these extensions, (a) the retardation factor can vary spatially within each model layer, (b) the decay rate coefficient can vary spatially within each model layer and can be different for the dissolved and sorbed phases, and (c) a zero-order growth reaction is added that can vary spatially and can be different in the dissolved and sorbed phases. The decay and growth reaction terms also can change in time to account for changing geochemical conditions during transport. The report includes a description of the theoretical basis of the model, a
Modeling groundwater flow on MPPs
Ashby, S.F.; Falgout, R.D.; Smith, S.G.; Tompson, A.F.B.
1993-10-01
The numerical simulation of groundwater flow in three-dimensional heterogeneous porous media is examined. To enable detailed modeling of large contaminated sites, preconditioned iterative methods and massively parallel computing power are combined in a simulator called PARFLOW. After describing this portable and modular code, some numerical results are given, including one that demonstrates the code`s scalability.
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…
Groundwater flow and transport modeling
Konikow, L.F.; Mercer, J.W.
1988-01-01
Deterministic, distributed-parameter, numerical simulation models for analyzing groundwater flow and transport problems have come to be used almost routinely during the past decade. A review of the theoretical basis and practical use of groundwater flow and solute transport models is used to illustrate the state-of-the-art. Because of errors and uncertainty in defining model parameters, models must be calibrated to obtain a best estimate of the parameters. For flow modeling, data generally are sufficient to allow calibration. For solute-transport modeling, lack of data not only limits calibration, but also causes uncertainty in process description. Where data are available, model reliability should be assessed on the basis of sensitivity tests and measures of goodness-of-fit. Some of these concepts are demonstrated by using two case histories. ?? 1988.
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.
NASA Astrophysics Data System (ADS)
Wang, Weishi; Munz, Matthias; Oswald, Sascha; Strasser, Daniel; Lensing, Hermann
2016-04-01
Bank filtration, by its effective improvement of water quality is widely used in many countries for water supply, and its major characteristics, the interaction between groundwater and surface water has been a hot topic for decades. As a key parameter, the travel time of the infiltrating river water to the wells is considered to be highly correlated with its water quality and has always been used as a main reference for estimating the filtering performance. As a periodic environmental tracer, heat has been used for estimating travel times by comparing the attenuation and the phase shift for temperature patterns in both the river and groundwater observation points. In most cases, the methods applied are analytical time series analysis, or 2D and 3D groundwater models with homogeneous attributes, in which many details of geological discontinuity and heterogeneity might be missed and further decrease the reliability of model result. However in our study, the transient heat transport model was set up based on a calibrated transient groundwater model with complex and discontinuous geological structures referenced by available geological information. At the study area, a water work is placed hundreds of meters from a river. By the pumping induced hydraulic gradient, river water flows into pumping wells through the river bank and shallow aquifers. The unconsolidated impermeable glacial deposits of different glacial periods showed discontinuities in forms of geological windows and lenses. Referenced by 145 drillings and 7 geological cross-sections, a geological model was set up and further translated into a groundwater model in FEFLOW. The model was first calibrated by FEPEST in steady state referenced by 104 observation wells and then it was adapted into a transient model. Influenced by an excavation at the channel bottom, a substantial water head rise happened. And in the model this could be simulated well by introducing an increasing hydraulic conductivity 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.
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.
NASA Astrophysics Data System (ADS)
Zanini, A.; Tanda, M.
2007-12-01
The groundwater in Italy plays an important role as drinking water; in fact it covers about the 30% of the national demand (70% in Northern Italy). The mineral water distribution in Italy is an important business with an increasing demand from abroad countries. The mineral water Companies have a great interest in order to increase the water extraction, but for the delicate and complex geology of the subsoil, where such very high quality waters are contained, a particular attention must be paid in order to avoid an excessive lowering of the groundwater reservoirs or great changes in the groundwater flow directions. A big water Company asked our University to set up a numerical model of the groundwater basin, in order to obtain a useful tool which allows to evaluate the strength of the aquifer and to design new extraction wells. The study area is located along Appennini Mountains and it covers a surface of about 18 km2; the topography ranges from 200 to 600 m a.s.l.. In ancient times only a spring with naturally sparkling water was known in the area, but at present the mineral water is extracted from deep pumping wells. The area is characterized by a very complex geology: the subsoil structure is described by a sequence of layers of silt-clay, marl-clay, travertine and alluvial deposit. Different groundwater layers are present and the one with best quality flows in the travertine layer; the natural flow rate seems to be not subjected to seasonal variations. The water age analysis revealed a very old water which means that the mineral aquifers are not directly connected with the meteoric recharge. The Geologists of the Company suggest that the water supply of the mineral aquifers comes from a carbonated unit located in the deep layers of the mountains bordering the spring area. The valley is crossed by a river that does not present connections to the mineral aquifers. Inside the area there are about 30 pumping wells that extract water at different depths. We built a 3
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.
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.
NASA Astrophysics Data System (ADS)
Goderniaux, Pascal; Davy, Philippe; Bresciani, Etienne; Dreuzy, Jean-Raynald; Borgne, Tanguy
2013-04-01
The distribution of groundwater fluxes in aquifers is strongly influenced by topography, and organized between hillslope and regional scales. The objective of this study is to provide new insights regarding the compartmentalization of aquifers at the regional scale and the partitioning of recharge between shallow/local and deep/regional groundwater transfers. A finite-difference flow model was implemented, and the flow structure was analyzed as a function of recharge (from 20 to 500 mm/yr), at the regional-scale (1400 km2), in three dimensions, and accounting for variable groundwater discharge zones; aspects which are usually not considered simultaneously in previous studies. The model allows visualizing 3-D circulations, as those provided by Tothian models in 2-D, and shows local and regional transfers, with 3-D effects. The probability density function of transit times clearly shows two different parts, interpreted using a two-compartment model, and related to regional groundwater transfers and local groundwater transfers. The role of recharge on the size and nature of the flow regimes, including groundwater pathways, transit time distributions, and volumes associated to the two compartments, have been investigated. Results show that topography control on the water table and groundwater compartmentalization varies with the recharge rate applied. When recharge decreases, the absolute value of flow associated to the regional compartment decreases, whereas its relative value increases. The volume associated to the regional compartment is calculated from the exponential part of the two-compartment model, and is nearly insensitive to the total recharge fluctuations.
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.
NASA Astrophysics Data System (ADS)
Fuchsluger, Martin; Götzl, Gregor
2014-05-01
flow has been realized. In addition the effects of the basement of the building to the groundwater flow have been analyzed. The results of the 2D model show an underestimation of more than 10 % of the performance of the groundwater utilization facility and a considerable smaller groundwater table drawdown compared to the 3D simulations. This is due to the possibility of 3D modeling to consider (i) the heat distribution and storage in the adjacent layers, (ii) the climatic surface effect and (iii) vertical groundwater flow.
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
Flow partitioning in regional groundwater flow systems as a function of recharge and topography
NASA Astrophysics Data System (ADS)
Goderniaux, P.; Davy, P.; Bresciani, E.; De Dreuzy, J.; Le Borgne, T.
2013-12-01
The distribution of groundwater fluxes in aquifers is strongly influenced by topography, and organized between hillslope and regional scales. In this study, we use a finite-difference flow model to quantify the partitioning of recharge and compartmentalization of aquifers between shallow/local and deep/regional groundwater transfers. The flow structure is analyzed for a regional aquifers, as a function of recharge (from 20 to 500 mm/yr), in 3-dimensions, and accounting for variable groundwater discharge zones. The Probability Density Function of transit times shows two different parts, interpreted using a two-compartment model, related to regional and local groundwater flows. The role of recharge on the size and nature of the flow regimes, including groundwater pathways, transit time distributions, and volumes associated to the two compartments is investigated. Results show that topography control on the water table and groundwater compartmentalization varies with the recharge rate applied. The volume associated to the regional compartment is calculated from the exponential part of the two-compartment model, and is nearly insensitive to the total recharge fluctuations. The model also allows visualizing 3D circulations, as those provided by Tothian models in 2D, and shows local and regional transfers, with 3D effects. Results are presented for a specific basin (1400 km2) in Brittany (France). Preliminary results using different kinds of topography are presented and compared.
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.
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.
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.
West Maui Groundwater Flow Model
Nicole Lautze
2015-01-01
Groundwater flow model for West Maui. Data is from the following sources: Whittier, R. and A.I. El-Kadi. 2014. Human and Environmental Risk Ranking of Onsite Sewage Disposal Systems For the Hawaiian Islands of Kauai, Molokai, Maui, and Hawaii – Final. Prepared by the University of Hawaii, Dept. of Geology and Geophysics for the State of Hawaii Dept. of Health, Safe Drinking Water Branch. September 2014; and Whittier, R.B., K. Rotzoll, S. Dhal, A.I. El-Kadi, C. Ray, G. Chen, and D. Chang. 2004. Hawaii Source Water Assessment Program Report – Volume V – Island of Maui Source Water Assessment Program Report. Prepared for the Hawaii Department of Health, Safe Drinking Water Branch. University of Hawaii, Water Resources Research Center. Updated 2008.
Hawaii Island Groundwater Flow Model
Nicole Lautze
2015-01-01
Groundwater flow model for Hawaii Island. Data is from the following sources: Whittier, R.B., K. Rotzoll, S. Dhal, A.I. El-Kadi, C. Ray, G. Chen, and D. Chang. 2004. Hawaii Source Water Assessment Program Report – Volume II – Island of Hawaii Source Water Assessment Program Report. Prepared for the Hawaii Department of Health, Safe Drinking Water Branch. University of Hawaii, Water Resources Research Center. Updated 2008; and Whittier, R. and A.I. El-Kadi. 2014. Human and Environmental Risk Ranking of Onsite Sewage Disposal Systems For the Hawaiian Islands of Kauai, Molokai, Maui, and Hawaii – Final. Prepared by the University of Hawaii, Dept. of Geology and Geophysics for the State of Hawaii Dept. of Health, Safe Drinking Water Branch. September 2014.
East Maui Groundwater Flow Model
Nicole Lautze
2015-01-01
Groundwater flow model for East Maui. Data is from the following sources: Whittier, R. and A.I. El-Kadi. 2014. Human and Environmental Risk Ranking of Onsite Sewage Disposal Systems For the Hawaiian Islands of Kauai, Molokai, Maui, and Hawaii – Final. Prepared by the University of Hawaii, Dept. of Geology and Geophysics for the State of Hawaii Dept. of Health, Safe Drinking Water Branch. September 2014; and Whittier, R.B., K. Rotzoll, S. Dhal, A.I. El-Kadi, C. Ray, G. Chen, and D. Chang. 2004. Hawaii Source Water Assessment Program Report – Volume V – Island of Maui Source Water Assessment Program Report. Prepared for the Hawaii Department of Health, Safe Drinking Water Branch. University of Hawaii, Water Resources Research Center. Updated 2008.
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.
NASA Astrophysics Data System (ADS)
James, A. L.; McDonnell, J. J.; Tromp-van Meerveld, H.
2006-12-01
Recent experimental work at the Panola Mountain Experimental Watershed trenched hillslope (Atlanta, GA) has shown that formation of transient groundwater in bedrock depressions during storm events is a precondition for the generation of subsurface stormflow. While there are clear threshold relations between subsurface stormflow and storm total precipitation, resolving the linkages between rainfall input, soil water recharge, transient groundwater formation and resulting subsurface stormflow has been difficult. Part of the problem is that we lack observational capabilities to interrogate the distributed internal slope response to storm rainfall. Even more problematic is our inability to quantify key boundary conditions like the permeability contrast at the soil-bedrock interface and its spatial variability. This paper uses a 3-D physically-based distributed model of the Panola hillslope (the TOUGH2 simulator) to model and visualize the linkages between storm rainfall, soil water recharge, transient groundwater development and resulting subsurface stormflow. We explore the effect of the bedrock permeability on the development of transient saturation within the hillslope in response to observed storm events. Preliminary results indicate that an increase in the estimated characteristic bedrock permeability at the Panola hillslope can deteriorate and even eliminate the connected patterns of transient saturation developed during a storm event that result in subsurface storm runoff. Our 3-D visualizations and virtual experiments with different bedrock permeability values provide insight into how bedrock permeability, antecedent soil moisture and storm conditions conspire to create the patterns of transient groundwater and resulting subsurface stormflow.
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
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.
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.
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.
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.
Ducci, Daniela; Sellerino, Mariangela
2013-03-01
The aim of this paper is to apply a methodology in order to reconstruct a lithostratigraphic 3D model of an aquifer so as to define some parameters involved in the evaluation of the aquifer vulnerability to contamination of porous aquifers. The DRASTIC, SINTACS and AVI methods have been applied to an alluvial coastal aquifer of southern Italy. The stratigraphic reconstruction has been obtained by interpolating stratigraphic data from more than one borehole per 2 km. The lithostratigraphic reconstruction of a 3D model has been applied and used for three-dimensional or two-dimensional representations. In the first two methods, the layers of the vadose zone and the aquifer media have been evaluated not only by the interpolation of the single boreholes and piezometers, but also by the 3D model, assigning the scores of the parameters of each layer of the 3D model. The comparison between the maps constructed from the weighted values in each borehole and the maps deriving from the attribution of the values of each layer of the 3D model, highlights that the second representation avoids or minimizes the "bullseye" effect linked to the presence of boreholes with higher or lower values. The study has demonstrated that it is possible to integrate a 3D lithostratigraphic model of an aquifer in the assessment of the parameters involved in the evaluation of the aquifer vulnerability to contamination by Point Count System methods. PMID:23391897
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.
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
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
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.
Influence of perched groundwater on base flow
Niswonger, R.G.; Fogg, G.E.
2008-01-01
Analysis with a three-dimensional variably saturated groundwater flow model provides a basic understanding of the interplay between streams and perched groundwater. A simplified, layered model of heterogeneity was used to explore these relationships. Base flow contribution from perched groundwater was evaluated with regard to varying hydrogeologic conditions, including the size and location of the fine-sediment unit and the hydraulic conductivity of the fine-sediment unit and surrounding coarser sediment. Simulated base flow was sustained by perched groundwater with a maximum monthly discharge in excess of 15 L/s (0.6 feet3/s) over the length of the 2000-m stream reach. Generally, the rate of perched-groundwater discharge to the stream was proportional to the hydraulic conductivity of sediment surrounding the stream, whereas the duration of discharge was proportional to the hydraulic conductivity of the fine-sediment unit. Other aspects of the perched aquifer affected base flow, such as the depth of stream penetration and the size of the fine-sediment unit. Greater stream penetration decreased the maximum base flow contribution but increased the duration of contribution. Perched groundwater provided water for riparian vegetation at the demand rate but reduced the duration of perched-groundwater discharge nearly 75%. Copyright 2008 by the American Geophysical Union.
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
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.
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
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.
STRING 3: An Advanced Groundwater Flow Visualization Tool
NASA Astrophysics Data System (ADS)
Schröder, Simon; Michel, Isabel; Biedert, Tim; Gräfe, Marius; Seidel, Torsten; König, Christoph
2016-04-01
The visualization of 3D groundwater flow is a challenging task. Previous versions of our software STRING [1] solely focused on intuitive visualization of complex flow scenarios for non-professional audiences. STRING, developed by Fraunhofer ITWM (Kaiserslautern, Germany) and delta h Ingenieurgesellschaft mbH (Witten, Germany), provides the necessary means for visualization of both 2D and 3D data on planar and curved surfaces. In this contribution we discuss how to extend this approach to a full 3D tool and its challenges in continuation of Michel et al. [2]. This elevates STRING from a post-production to an exploration tool for experts. In STRING moving pathlets provide an intuition of velocity and direction of both steady-state and transient flows. The visualization concept is based on the Lagrangian view of the flow. To capture every detail of the flow an advanced method for intelligent, time-dependent seeding is used building on the Finite Pointset Method (FPM) developed by Fraunhofer ITWM. Lifting our visualization approach from 2D into 3D provides many new challenges. With the implementation of a seeding strategy for 3D one of the major problems has already been solved (see Schröder et al. [3]). As pathlets only provide an overview of the velocity field other means are required for the visualization of additional flow properties. We suggest the use of Direct Volume Rendering and isosurfaces for scalar features. In this regard we were able to develop an efficient approach for combining the rendering through raytracing of the volume and regular OpenGL geometries. This is achieved through the use of Depth Peeling or A-Buffers for the rendering of transparent geometries. Animation of pathlets requires a strict boundary of the simulation domain. Hence, STRING needs to extract the boundary, even from unstructured data, if it is not provided. In 3D we additionally need a good visualization of the boundary itself. For this the silhouette based on the angle of
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.
NASA Astrophysics Data System (ADS)
Cao, Guoliang; Han, Dongmei; Currell, Matthew J.; Zheng, Chunmiao
2016-09-01
Groundwater flow in deep sedimentary basins results from complex evolution processes on geological timescales. Groundwater flow systems conceptualized according to topography and/or groundwater table configuration generally assume a near-equilibrium state with the modern landscape. However, the time to reach such a steady state, and more generally the timescales of groundwater flow system evolution are key considerations for large sedimentary basins. This is true in the North China Basin (NCB), which has been studied for many years due to its importance as a groundwater supply. Despite many years of study, there remain contradictions between the generally accepted conceptual model of regional flow, and environmental tracer data. We seek to reconcile these contractions by conducting simulations of groundwater flow, age and heat transport in a three dimensional model, using an alternative conceptual model, based on geological, thermal, isotope and historical data. We infer flow patterns under modern hydraulic conditions using this new model and present the theoretical maximum groundwater ages under such a flow regime. The model results show that in contrast to previously accepted conceptualizations, most groundwater is discharged in the vicinity of the break-in-slope of topography at the boundary between the piedmont and central plain. Groundwater discharge to the ocean is in contrast small, and in general there are low rates of active flow in the eastern parts of the basin below the central and coastal plain. This conceptualization is more compatible with geochemical and geothermal data than the previous model. Simulated maximum groundwater ages of ∼1 Myrs below the central and coastal plain indicate that residual groundwater may be retained in the deep parts of the basin since being recharged during the last glacial period or earlier. The groundwater flow system has therefore probably not reached a new equilibrium state with modern-day hydraulic conditions. The
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.
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.
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.
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.
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.
Simulation of groundwater flow around a pilot waste tire site
Lin, Y.; Chyi, L.L. . Dept. of Geology)
1994-04-01
The groundwater system in the area contains several aquifers divided by underclays as confining layers. The three aquifers addressed under this research include the upper aquifer above the No. 6 underclay, the middle aquifer above the No. 5 underclay and the lower aquifer above the No. 4 underclay. A quasi-3D conceptual model is used in simulating groundwater flow and MODFLOW is used for the actual modeling. The boundary of the model area is defined by a local divide on the north, the Hugle Run on the west and the Temple Ditch on the east and the Sandy Creek on the south. The model consists of three aquifer layers and two confining layers. The upper one is unconfined, the middle one is a hybrid, and the lower one is confined. In the simulation of the upper aquifer, drain boundary condition is defined to include seepage. In the simulation of the middle aquifer, the Hugle Run and the Temple ditch associated with the alluvium represent a drain boundary and the Sand Creek to the south represents a constant head. The lower aquifer does not have a natural hydrogeologic boundary. The equipotential contours are assigned as constant head boundary. Simulation has indicated that minimal groundwater flow existed in the upper aquifer and it is flowing toward the southwest. Groundwater flow in the middle aquifer near the pilot waste site is somewhat divergent to the south. The velocity is slow near the site but increases quickly toward the southwest. Leakage from the upper aquifer appears to be predominant over the regional flow in the recharge of the middle aquifer. The lower aquifer is recharged mainly by the regional flow and is flowing almost due south.
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.
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
Considering Barometric Pressure in Groundwater Flow Investigations
Spane, Frank A. )
2002-06-18
Well water-level elevation measurements are commonly used as the basis to delineate groundwater-flow patterns (i.e., flow direction and hydraulic gradient). Barometric pressure fluctuations, however, can have a discernable impact on well water-levels. These barometric effects may lead to erroneous indications of hydraulic head within the aquifer. Total hydraulic head within the aquifer, not well water-level elevation, is the hydrologic parameter for determining groundwater-flow direction and hydraulic gradient conditions. For low-gradient, unconfined aquifer sites exhibiting variable vadose zone characteristics (e.g., thickness, pneumatic diffusivity), barometric pressure fluctuations can also cause actual, temporal changes in lateral flow direction and flow velocity. Discrete water-level measurements used to determine the average or long-term groundwater-flow conditions, therefore, may provide non-representative results. Calculation of the barometric response characteristics for individual wells provides the basis to account for the temporal effects of barometric pressure fluctuations from monitor well measurements, so that average, long-term groundwater-flow pattern behavior can be determined.
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.
Effects of wetlands creation on groundwater flow
Hensel, B.R.; Miller, M.V.
1991-01-01
Changes in groundwater flow were observed near four Experimental Wetland Areas (EWAs) constructed along a reach of the Des Plaines River in northeastern Illinois. These changes were observed during monthly monitoring of groundwater elevation in nested piezometers and shallow observation wells before and after the wetlands were filled with water. A numerical model was calibrated with observed data and used to estimate seepage from the wetlands into the Des Plaines River. After the wetlands became operational, groundwater levels in adjacent wells increased by about 0.5m, while water levels in wells distant from the wetlands decreased. The increase in groundwater levels near the wetlands is a result of seepage from the wetlands. Numerical predictions of seepage from the wetlands are 60-150 m3 day-1 for two wetlands situated over sand and gravel and less than 1 m3 day-1 for two wetlands situated over clayey till. The difference in seepage rates is attributed to two factors. First, the hydraulic conductivity of the sand and gravel unit is greater than that of the till, and thus there is less mounding and a greater capacity for transmitting water beneath the wetlands overlying this deposit. Secondly, the wetlands located over till are groundwater flow-through ponds, whereas the wetlands over the sand and gravel are primarily groundwater recharge areas. The model was used to estimate that seepage from the wetlands will double groundwater discharge into the Des Plaines River and a tributary relative to pre-operational discharge from the study area. Overall, the wetlands have acted as a constant head boundary, stabilizing groundwater flow patterns. ?? 1991.
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.
NASA Astrophysics Data System (ADS)
Tóth, Ádám; Havril, Tímea; Simon, Szilvia; Galsa, Attila; Monteiro Santos, Fernando A.; Müller, Imre; Mádl-Szőnyi, Judit
2016-08-01
Groundwater flow, driven, controlled and determined by topography, geology and climate, is responsible for several natural surface manifestations and affected by anthropogenic processes. Therefore, flowing groundwater can be regarded as an environmental agent. Numerical simulation of groundwater flow could reveal the flow pattern and explain the observed features. In complex geologic framework, where the geologic-hydrogeologic knowledge is limited, the groundwater flow model could not be constructed based solely on borehole data, but geophysical information could aid the model building. The integrated model construction was presented via the case study of the Tihany Peninsula, Hungary, with the aims of understanding the background and occurrence of groundwater-related environmental phenomena, such as wetlands, surface water-groundwater interaction, slope instability, and revealing the potential effect of anthropogenic activity and climate change. The hydrogeologic model was prepared on the basis of the compiled archive geophysical database and the results of recently performed geophysical measurements complemented with geologic-hydrogeologic data. Derivation of different electrostratigraphic units, revealing fracturing and detecting tectonic elements was achieved by systematically combined electromagnetic geophysical methods. The deduced information can be used as model input for groundwater flow simulation concerning hydrostratigraphy, geometry and boundary conditions. The results of numerical modelling were interpreted on the basis of gravity-driven regional groundwater flow concept and validated by field mapping of groundwater-related phenomena. The 3D model clarified the hydraulic behaviour of the formations, revealed the subsurface hydraulic connection between groundwater and wetlands and displayed the groundwater discharge pattern, as well. The position of wetlands, their vegetation type, discharge features and induced landslides were explained as
Regression modeling of ground-water flow
Cooley, R.L.; Naff, R.L.
1985-01-01
Nonlinear multiple regression methods are developed to model and analyze groundwater flow systems. Complete descriptions of regression methodology as applied to groundwater flow models allow scientists and engineers engaged in flow modeling to apply the methods to a wide range of problems. Organization of the text proceeds from an introduction that discusses the general topic of groundwater flow modeling, to a review of basic statistics necessary to properly apply regression techniques, and then to the main topic: exposition and use of linear and nonlinear regression to model groundwater flow. Statistical procedures are given to analyze and use the regression models. A number of exercises and answers are included to exercise the student on nearly all the methods that are presented for modeling and statistical analysis. Three computer programs implement the more complex methods. These three are a general two-dimensional, steady-state regression model for flow in an anisotropic, heterogeneous porous medium, a program to calculate a measure of model nonlinearity with respect to the regression parameters, and a program to analyze model errors in computed dependent variables such as hydraulic head. (USGS)
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.
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.
GROUNDWATER FLOW IN LOW-PERMEABILITY ENVIRONMENTS.
Neuzil, C.E.
1986-01-01
Certain geologic media are known to have small permeability; subsurface environments composed of these media and lacking well developed secondary permeability have groundwater flow systems with many distinctive characteristics. Moreover, groundwater flow in these environments appears to influence the evolution of certain hydrologic, geologic, and geochemical systems, may affect the accumulation of petroleum and ores, and probably has a role in the structural evolution of parts of the crust. Such environments are also important in the context of waste disposal. This review attempts to synthesize the diverse contributions of various disciplines to the problem of flow in low-permeability environments. Problems hindering analysis are enumerated together with suggested approaches to overcoming them. A common thread running through the discussion is the significance of size- and time-scale limitations of the ability to directly observe flow behavior and make significance of size- and time-scale limitations of the ability to directly observe flow behavior and make measurements of parameters.
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.
NASA Astrophysics Data System (ADS)
Iwasaki, Yumi; Nakamura, Kimihito; Horino, Haruhiko; Kawashima, Shigeto
2014-12-01
Assessing factors that influence groundwater levels such as land use and pumping strategy, is essential to adequately manage groundwater resources. A transient numerical model for groundwater flow with infiltration was developed for the Tedori River alluvial fan (140 km2), Japan. The main water input into the groundwater body in this area is irrigation water, which is significantly influenced by land use, namely paddy and upland fields. The proposed model consists of two models, a one-dimensional (1-D) unsaturated-zone water flow model (HYDRUS-1D) for estimating groundwater recharge and a 3-D groundwater flow model (MODFLOW). Numerical simulation of groundwater flow from October 1975 to November 2009 was performed to validate the model. Simulation revealed seasonal groundwater level fluctuations, affected by paddy irrigation management. However, computational accuracy was limited by the spatiotemporal data resolution of the groundwater use. Both annual groundwater levels and recharge during the irrigation periods from 1975 to 2009 showed long-term decreasing trends. With the decline in rice-planted paddy field area, groundwater recharge cumulatively decreased to 61 % of the peak in 1977. A paddy-upland crop-rotation system could decrease groundwater recharge to 73-98 % relative to no crop rotation.
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.
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
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.
Groundwater Flow and Arsenic Biogeochemistry in Bangladesh
NASA Astrophysics Data System (ADS)
Harvey, C. F.
2004-12-01
Although groundwater in Bangladesh is severely contaminated by arsenic, little is known about the complex transient patterns of groundwater flow that flush solutes from aquifers and carry solutes into the subsurface. Hydrologic modeling results for our field site in the Munshiganj district indicate that groundwater flow is vigorous, flushing the aquifer over time-scales of decades and also introducing solute loads into the aquifer with recharge from rice fields, ponds and rivers. The combined hydrologic and biogeochemical results from our field site imply that the biogeochemistry of the aquifer system may not be in steady-state, and that the net effect of competing processes could either increase or decrease arsenic concentrations over the next decades. Modeling results suggest that irrigation has greatly changed the location, timing and chemical content of recharge to the aquifer, drawing large fluxes of anoxic water into the aquifer during the dry season that may mobilize arsenic from oxides in near-surface sediments.
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.
Geomorphic aspects of groundwater flow
NASA Astrophysics Data System (ADS)
LaFleur, Robert G.
The many roles that groundwater plays in landscape evolution are becoming more widely appreciated. In this overview, three major categories of groundwater processes and resulting landforms are considered: (1) Dissolution creates various karst geometries, mainly in carbonate rocks, in response to conditions of recharge, geologic setting, lithology, and groundwater circulation. Denudation and cave formation rates can be estimated from kinetic and hydraulic parameters. (2) Groundwater weathering generates regoliths of residual alteration products at weathering fronts, and subsequent exhumation exposes corestones, flared slopes, balanced rocks, domed inselbergs, and etchplains of regional importance. Groundwater relocation of dissolved salts creates duricrusts of various compositions, which become landforms. (3) Soil and rock erosion by groundwater processes include piping, seepage erosion, and sapping, important agents in slope retreat and headward gully migration. Thresholds and limits are important in many chemical and mechanical groundwater actions. A quantitative, morphometric approach to groundwater landforms and processes is exemplified by selected studies in carbonate and clastic terrains of ancient and recent origins. Résumé Les rôles variés joués par les eaux souterraines dans l'évolution des paysages deviennent nettement mieux connus. La revue faite ici prend en considération trois grandes catégories de processus liés aux eaux souterraines et les formes associées: (1) La dissolution crée des formes karstiques variées, surtout dans les roches carbonatées, en fonction des conditions d'alimentation, du cadre géologique, de la lithologie et de la circulation des eaux souterraines. Les taux d'érosion et de formation des grottes peuvent être estimés à partir de paramètres cinétiques et hydrauliques. (2) L'érosion par les eaux souterraines donne naissance à des régolites, résidus d'altération sur des fronts d'altération, et l'exhumation r
Geomorphic aspects of groundwater flow
NASA Astrophysics Data System (ADS)
LaFleur, Robert G.
The many roles that groundwater plays in landscape evolution are becoming more widely appreciated. In this overview, three major categories of groundwater processes and resulting landforms are considered: (1) Dissolution creates various karst geometries, mainly in carbonate rocks, in response to conditions of recharge, geologic setting, lithology, and groundwater circulation. Denudation and cave formation rates can be estimated from kinetic and hydraulic parameters. (2) Groundwater weathering generates regoliths of residual alteration products at weathering fronts, and subsequent exhumation exposes corestones, flared slopes, balanced rocks, domed inselbergs, and etchplains of regional importance. Groundwater relocation of dissolved salts creates duricrusts of various compositions, which become landforms. (3) Soil and rock erosion by groundwater processes include piping, seepage erosion, and sapping, important agents in slope retreat and headward gully migration. Thresholds and limits are important in many chemical and mechanical groundwater actions. A quantitative, morphometric approach to groundwater landforms and processes is exemplified by selected studies in carbonate and clastic terrains of ancient and recent origins. Résumé Les rôles variés joués par les eaux souterraines dans l'évolution des paysages deviennent nettement mieux connus. La revue faite ici prend en considération trois grandes catégories de processus liés aux eaux souterraines et les formes associées: (1) La dissolution crée des formes karstiques variées, surtout dans les roches carbonatées, en fonction des conditions d'alimentation, du cadre géologique, de la lithologie et de la circulation des eaux souterraines. Les taux d'érosion et de formation des grottes peuvent être estimés à partir de paramètres cinétiques et hydrauliques. (2) L'érosion par les eaux souterraines donne naissance à des régolites, résidus d'altération sur des fronts d'altération, et l'exhumation r
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.
3-D Particle Tracking Velocimetry: Development and Applications in Small Scale Flows
NASA Astrophysics Data System (ADS)
Tien, Wei-Hsin
The thesis contains two parts of studies. In part I, a novel volumetric velocimetry technique is developed to measure the 3-D flow field of small-scale flows. The technique utilizes a color-coded pinhole plate with multiple light sources aligned to each pinhole to achieve high particle image density and large measurable depth on a single lens microscope system. A color separation algorithm and an improved particle identification algorithm are developed to identify individual particle images from each pinhole view. Furthermore, a calibration-based technique based on epi-polar line search method is developed to reconstruct the spatial coordinates of the particle, and a new two-frame tracking particle-tracking algorithm is developed to calculate the velocity field. The system was setup to achieve a magnification of 2.69, resulting in an imaging volume of 3.35 x 2.5 x 1.5 mm3 and showed satisfactory measurement accuracy. The technique was then further miniaturized to achieve a magnification of 10, resulting in a imaging volume of 600 x 600 x 600 microm3. The system was applied to a backward-facing step flow to test its ability to reconstruct the unsteady flow field with two-frame tracking. Finally, this technique was applied to a steady streaming flow field in a microfluidic device used to trap particles. The results revealed the three-dimensional flow structure that has not been observed in previous studies, and provided insights to the design of a more efficient trapping device. In part II, an in-vitro study was carried out to investigate the flow around a prosthetic venous valve. Using 2-D PIV, the dynamics of the valve motion was captured and the velocity fields were measured to investigate the effect of the sinus pocket and the coupling effect of a pair of valves. The PIV and hemodynamic results showed that the sinus pocket around the valve functioned as a flow regulator to smooth the entrained velocity profile and suppress the jet width. For current prosthetic
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.
IN SITU FLOW METERS AROUND A GROUNDWATER CIRCULATION WELL (ABSTRACT)
The primary benefit of groundwater circulation well (GCW) technology is the development of strong vertical flows surrounding the treatment well. The extent of significant vertical flow surrounding a circulation well is difficult to establish from traditional groundwater elevation...
IN SITU FLOW METERS AROUND A GROUNDWATER CIRCULATION WELL (PAPER)
The primary benefit of groundwater circulation well (GCW) technology is the development of strong vertical flows surrounding the treatment well. The extent of significant vertical flow surrounding a circulation well is difficult to establish from traditional groundwater elevation...
3PE: A Tool for Estimating Groundwater Flow Vectors
Evaluation of hydraulic gradients and the associated groundwater flow rates and directions is a fundamental aspect of hydrogeologic characterization. Many methods, ranging in complexity from simple three-point solution techniques to complex numerical models of groundwater flow, ...
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.
Modeling groundwater flow on massively parallel computers
Ashby, S.F.; Falgout, R.D.; Fogwell, T.W.; Tompson, A.F.B.
1994-12-31
The authors will explore the numerical simulation of groundwater flow in three-dimensional heterogeneous porous media. An interdisciplinary team of mathematicians, computer scientists, hydrologists, and environmental engineers is developing a sophisticated simulation code for use on workstation clusters and MPPs. To date, they have concentrated on modeling flow in the saturated zone (single phase), which requires the solution of a large linear system. they will discuss their implementation of preconditioned conjugate gradient solvers. The preconditioners under consideration include simple diagonal scaling, s-step Jacobi, adaptive Chebyshev polynomial preconditioning, and multigrid. They will present some preliminary numerical results, including simulations of groundwater flow at the LLNL site. They also will demonstrate the code`s scalability.
Multiphase groundwater flow near cooling plutons
Hayba, D.O.; Ingebritsen, S.E.
1997-01-01
We investigate groundwater flow near cooling plutons with a computer program that can model multiphase flow, temperatures up to 1200??C, thermal pressurization, and temperature-dependent rock properties. A series of experiments examines the effects of host-rock permeability, size and depth of pluton emplacement, single versus multiple intrusions, the influence of a caprock, and the impact of topographically driven groundwater flow. We also reproduce and evaluate some of the pioneering numerical experiments on flow around plutons. Host-rock permeability is the principal factor influencing fluid circulation and heat transfer in hydrothermal systems. The hottest and most steam-rich systems develop where permeability is of the order of 10-15 m2. Temperatures and life spans of systems decrease with increasing permeability. Conduction-dominated systems, in which permeabilities are ???10-16m2, persist longer but exhibit relatively modest increases in near-surface temperatures relative to ambient conditions. Pluton size, emplacement depth, and initial thermal conditions have less influence on hydrothermal circulation patterns but affect the extent of boiling and duration of hydrothermal systems. Topographically driven groundwater flow can significantly alter hydrothermal circulation; however, a low-permeability caprock effectively decouples the topographically and density-driven systems and stabilizes the mixing interface between them thereby defining a likely ore-forming environment.
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
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.
Quasi 3D modelling of water flow in the sandy soil
NASA Astrophysics Data System (ADS)
Rezaei, Meisam; Seuntjens, Piet; Joris, Ingeborg; Boënne, Wesley; De Pue, Jan; Cornelis, Wim
2016-04-01
Monitoring and modeling tools may improve irrigation strategies in precision agriculture. Spatial interpolation is required for analyzing the effects of soil hydraulic parameters, soil layer thickness and groundwater level on irrigation management using hydrological models at field scale. We used non-invasive soil sensor, a crop growth (LINGRA-N) and a soil hydrological model (Hydrus-1D) to predict soil-water content fluctuations and crop yield in a heterogeneous sandy grassland soil under supplementary irrigation. In the first step, the sensitivity of the soil hydrological model to hydraulic parameters, water stress, crop yield and lower boundary conditions was assessed after integrating models at one soil column. Free drainage and incremental constant head conditions were implemented in a lower boundary sensitivity analysis. In the second step, to predict Ks over the whole field, the spatial distributions of Ks and its relationship between co-located soil ECa measured by a DUALEM-21S sensor were investigated. Measured groundwater levels and soil layer thickness were interpolated using ordinary point kriging (OK) to a 0.5 by 0.5 m in aim of digital elevation maps. In the third step, a quasi 3D modelling approach was conducted using interpolated data as input hydraulic parameter, geometric information and boundary conditions in the integrated model. In addition, three different irrigation scenarios namely current, no irrigation and optimized irrigations were carried out to find out the most efficient irrigation regime. In this approach, detailed field scale maps of soil water stress, water storage and crop yield were produced at each specific time interval to evaluate the best and most efficient distribution of water using standard gun sprinkler irrigation. The results show that the effect of the position of the groundwater level was dominant in soil-water content prediction and associated water stress. A time-dependent sensitivity analysis of the hydraulic
Connections between groundwater flow and transpiration partitioning.
Maxwell, Reed M; Condon, Laura E
2016-07-22
Understanding freshwater fluxes at continental scales will help us better predict hydrologic response and manage our terrestrial water resources. The partitioning of evapotranspiration into bare soil evaporation and plant transpiration remains a key uncertainty in the terrestrial water balance. We used integrated hydrologic simulations that couple vegetation and land-energy processes with surface and subsurface hydrology to study transpiration partitioning at the continental scale. Both latent heat flux and partitioning are connected to water table depth, and including lateral groundwater flow in the model increases transpiration partitioning from 47 ± 13 to 62 ± 12%. This suggests that lateral groundwater flow, which is generally simplified or excluded in Earth system models, may provide a missing link for reconciling observations and global models of terrestrial water fluxes. PMID:27463671
Connections between groundwater flow and transpiration partitioning
NASA Astrophysics Data System (ADS)
Maxwell, Reed M.; Condon, Laura E.
2016-07-01
Understanding freshwater fluxes at continental scales will help us better predict hydrologic response and manage our terrestrial water resources. The partitioning of evapotranspiration into bare soil evaporation and plant transpiration remains a key uncertainty in the terrestrial water balance. We used integrated hydrologic simulations that couple vegetation and land-energy processes with surface and subsurface hydrology to study transpiration partitioning at the continental scale. Both latent heat flux and partitioning are connected to water table depth, and including lateral groundwater flow in the model increases transpiration partitioning from 47 ± 13 to 62 ± 12%. This suggests that lateral groundwater flow, which is generally simplified or excluded in Earth system models, may provide a missing link for reconciling observations and global models of terrestrial water fluxes.
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.
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.
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.
Sneddon, Kristen W.; Powers, Michael H.; Johnson, Raymond H.; Poeter, Eileen P.
2002-01-01
Dense nonaqueous phase liquids (DNAPLs) are a pervasive and persistent category of groundwater contamination. In an effort to better understand their unique subsurface behavior, a controlled and carefully monitored injection of PCE (perchloroethylene), a typical DNAPL, was performed in conjunction with the University of Waterloo at Canadian Forces Base Borden in 1991. Of the various geophysical methods used to monitor the migration of injected PCE, the U.S. Geological Survey collected 500-MHz ground penetrating radar (GPR) data. These data are used in determining calibration parameters for a multiphase flow simulation. GPR data were acquired over time on a fixed two-dimensional surficial grid as the DNAPL was injected into the subsurface. Emphasis is on the method of determining DNAPL saturation values from this time-lapse GPR data set. Interactive full-waveform GPR modeling of regularized field traces resolves relative dielectric permittivity versus depth profiles for pre-injection and later-time data. Modeled values are end members in recursive calculations of the Bruggeman-Hanai-Sen (BHS) mixing formula, yielding interpreted pre-injection porosity and post-injection DNAPL saturation values. The resulting interpreted physical properties of porosity and DNAPL saturation of the Borden test cell, defined on a grid spacing of 50 cm with 1-cm depth resolution, are used as observations for calibration of a 3-D multiphase flow simulation. Calculated values of DNAPL saturation in the subsurface at 14 and 22 hours after the start of injection, from both the GPR and the multiphase flow modeling, are interpolated volumetrically and presented for visual comparison.
Groundwater flow as a cooling agent of the continental lithosphere
NASA Astrophysics Data System (ADS)
Kooi, Henk
2016-03-01
Groundwater that flows through the outer shell of the Earth as part of the hydrologic cycle influences the distribution of heat and, thereby, the temperature field in the Earth’s crust. Downward groundwater flow in recharge areas lowers crustal temperatures, whereas upward flow in discharge areas tends to raise temperatures relative to a purely conductive geothermal regime. Here I present numerical simulations of generalized topography-driven groundwater flow. The simulations suggest that groundwater-driven convective cooling exceeds groundwater-driven warming of the Earth’s crust, and hence that groundwater flow systems cause net temperature reductions of groundwater basins. Moreover, the simulations demonstrate that this cooling extends into the underlying crust and lithosphere. I find that horizontal components of groundwater flow play a central role in this net subsurface cooling by conveying relatively cold water to zones of upward groundwater flow. The model calculations suggest that the crust and lithosphere beneath groundwater basins can cool by several tens of degrees Celsius where groundwater flows over large distances in basins that consist of crustal rock. In contrast, groundwater-induced cooling is small in unconsolidated sedimentary settings, such as deltas.
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.
Coupled Groundwater and Heat Flow in the Tahoe Basin Region
NASA Astrophysics Data System (ADS)
Trask, J. C.
2002-12-01
We propose that recent developments in available temperature probe technology and improvements in appropriate modeling software, together with the advent of desktop high-speed computing, have enabled the thermal analysis approach to be an inexpensive, robust, and practical way to constrain groundwater flow fields in a wide variety of environments. We present field measurements and numerical models of coupled heat and groundwater flow in the Tahoe Basin region. In montane regions such as the Tahoe Basin, steep topography provides a driving force for deep groundwater flow. Deep groundwater flow re-routes subsurface heat flow, impacting temperature gradients to depth, including the shallow subsurface (<100m depth). In the Tahoe Basin region, the magnitude of deep groundwater flow on the areal or regional scale has been largely unknown. We present examples of borehole temperature profiles that constrain possible areal groundwater flow patterns, including the magnitude of flow beneath the bottom of boreholes probed.
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.
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.
Patterns in groundwater chemistry resulting from groundwater flow
NASA Astrophysics Data System (ADS)
Stuyfzand, Pieter J.
Groundwater flow influences hydrochemical patterns because flow reduces mixing by diffusion, carries the chemical imprints of biological and anthropogenic changes in the recharge area, and leaches the aquifer system. Global patterns are mainly dictated by differences in the flux of meteoric water passing through the subsoil. Within individual hydrosomes (water bodies with a specific origin), the following prograde evolution lines (facies sequence) normally develop in the direction of groundwater flow: from strong to no fluctuations in water quality, from polluted to unpolluted, from acidic to basic, from oxic to anoxic-methanogenic, from no to significant base exchange, and from fresh to brackish. This is demonstrated for fresh coastal-dune groundwater in the Netherlands. In this hydrosome, the leaching of calcium carbonate as much as 15m and of adsorbed marine cations (Na+, K+, and Mg2+) as much as 2500m in the flow direction is shown to correspond with about 5000yr of flushing since the beach barrier with dunes developed. Recharge focus areas in the dunes are evidenced by groundwater displaying a lower prograde quality evolution than the surrounding dune groundwater. Artificially recharged Rhine River water in the dunes provides distinct hydrochemical patterns, which display groundwater flow, mixing, and groundwater ages. Résumé Les écoulements souterrains influencent les différents types hydrochimiques, parce que l'écoulement réduit le mélange par diffusion, porte les marques chimiques de changements biologiques et anthropiques dans la zone d'alimentation et lessive le système aquifère. Ces types dans leur ensemble sont surtout déterminés par des différences dans le flux d'eau météorique traversant le sous-sol. Dans les "hydrosomes" (masses d'eau d'origine déterminée), les lignes marquant une évolution prograde (séquence de faciès) se développent normalement dans la direction de l'écoulement souterrain : depuis des fluctuations fortes de la
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.
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
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.
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
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
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.
Groundwater Flow Model for Taos, New Mexico
NASA Astrophysics Data System (ADS)
Burck, P. W.; Barroll, P. W.; Core, A. B.; Rappuhn, D.
2003-12-01
The New Mexico Office of the State Engineer - Hydrology Bureau (OSE) has developed a regional groundwater flow model for Taos, New Mexico. The MODFLOW 2000 model will serve as a tool to evaluate alternatives in settlement negotiations in an on-going water rights adjudication. If current settlement negotiations fail, it is conceivable that the model might be used in support of litigation. OSE produced the model in cooperation with technical representatives of the various parties to the adjudication. Regional hydrogeologic data including well records, aquifer test results, stream flow measurements and seepage studies have been shared relatively freely among the parties. A recent deep drilling program conducted in conjunction with the negotiation effort has added substantially to the hydrogeologic data set. Among the hydrologic processes simulated by the model are mountain front recharge; areal recharge from precipitation; evapotranspiration; discharge from springs; river and stream flow; accretions to groundwater from irrigation return flow, seepage from acequias, canals, and ditches, and deep percolation; and pumping by municipal entities and mutual domestic water users associations. The resulting model files are available for all parties to review and evaluate. Comments are assessed and many have resulted in significant improvements to the model. At this stage, however, it is unclear whether adopting this cooperative approach will increase the likelihood of model acceptance by the parties.