Ito, Kazumasa; Yongkoo, Seol
2003-04-09
Water fluxes in unsaturated, fractured rock involve the physical processes occurring at fracture-matrix interfaces within fracture networks. Modeling these water fluxes using a discrete fracture network model is a complicated effort. Existing preprocessors for TOUGH2 are not suitable to generate grids for fracture networks with various orientations and inclinations. There are several 3-D discrete-fracture-network simulators for flow and transport, but most of them do not capture fracture-matrix interaction. We have developed a new 3-D discrete-fracture-network mesh generator, FRACMESH, to provide TOUGH2 with information about the fracture network configuration and fracture-matrix interactions. FRACMESH transforms a discrete fracture network into a 3 dimensional uniform mesh, in which fractures are considered as elements with unique rock material properties and connected to surrounding matrix elements. Using FRACMESH, individual fractures may have uniform or random aperture distributions to consider heterogeneity. Fracture element volumes and interfacial areas are calculated from fracture geometry within individual elements. By using FRACMESH and TOUGH2, fractures with various inclinations and orientations, and fracture-matrix interaction, can be incorporated. In this paper, results of flow and transport simulations in a fractured rock block utilizing FRACMESH are presented.
Efficacy of 3-Dimensional plates over Champys miniplates in mandibular anterior fractures
Barde, Dhananjay H; Mudhol, Anupama; Ali, Fareedi Mukram; Madan, R S; Kar, Sanjay; Ustaad, Farheen
2014-01-01
Background: Mandibular fractures are treated surgically by either rigid or semi-rigid fixation, two techniques that reflect almost opposite concept of craniomaxillofacial osteosynthesis. The shortcomings of these fixations led to the development of 3 dimensional (3D) miniplates. This study was designed with the aim of evaluating the efficiency of 3D miniplate over Champys miniplate in anterior mandibular fractures. Materials & Methods: This study was done in 40 patients with anterior mandibular fractures. Group I consisting of 20 patients in whom 3D plates were used for fixation while in Group II consisting of other 20 patients, 4 holes straight plates were used. The efficacy of 3D miniplate over Champy’s miniplate was evaluated in terms of operating time, average pain, post operative infection, occlusion, wound dehiscence, post operative mobility and neurological deficit. Results: The mean operation time for Group II was more compared to Group I (statistically significant).There was significantly greater pain on day of surgery and at 2nd week for Group II patients but there was no significant difference between the two groups at 4th week. The post operative infection, occlusal disturbance, wound dehiscence, post operative mobility at facture site, neurological deficit was statistically insignificant (chi square test). Conclusion: The results of this study suggest that fixation of anterior mandibular fractures with 3D plates provides three dimensional stability and carries low morbidity and infection rates. The only probable limitation of these 3D plates may be excessive implant material, but they seem to be easy alternative to champys miniplate. How to cite the article: Barde DH, Mudhol A, Ali FM, Madan RS, Kar S, Ustaad F. Efficacy of 3-Dimensional plates over Champys miniplates in mandibular anterior fractures. J Int Oral Health 2014;6(1):20-6. PMID:24653598
Compartmentalization analysis using discrete fracture network models
La Pointe, P.R.; Eiben, T.; Dershowitz, W.; Wadleigh, E.
1997-08-01
This paper illustrates how Discrete Fracture Network (DFN) technology can serve as a basis for the calculation of reservoir engineering parameters for the development of fractured reservoirs. It describes the development of quantitative techniques for defining the geometry and volume of structurally controlled compartments. These techniques are based on a combination of stochastic geometry, computational geometry, and graph the theory. The parameters addressed are compartment size, matrix block size and tributary drainage volume. The concept of DFN models is explained and methodologies to compute these parameters are demonstrated.
Compartmentalization analysis using discrete fracture network models
La Pointe, P.R.; Eiben, T.; Dershowitz, W.; Wadleigh, E.
1997-12-31
This paper illustrates how Discrete Fracture Network (DFN) technology can serve as a basis for the calculation of reservoir engineering parameters for the development of fractured reservoirs. It describes the development of quantitative techniques for defining the geometry and volume of structurally controlled compartments. These techniques are based on a combination of stochastic geometry, computational geometry, and graph theory. The parameters addressed are compartment size, matrix block size and tributary drainage volume. The concept of DFN models is explained and methodologies to compute these parameters are demonstrated.
Discrete modeling of hydraulic fracturing processes in a complex pre-existing fracture network
NASA Astrophysics Data System (ADS)
Kim, K.; Rutqvist, J.; Nakagawa, S.; Houseworth, J. E.; Birkholzer, J. T.
2015-12-01
Hydraulic fracturing and stimulation of fracture networks are widely used by the energy industry (e.g., shale gas extraction, enhanced geothermal systems) to increase permeability of geological formations. Numerous analytical and numerical models have been developed to help understand and predict the behavior of hydraulically induced fractures. However, many existing models assume simple fracturing scenarios with highly idealized fracture geometries (e.g., propagation of a single fracture with assumed shapes in a homogeneous medium). Modeling hydraulic fracture propagation in the presence of natural fractures and homogeneities can be very challenging because of the complex interactions between fluid, rock matrix, and rock interfaces, as well as the interactions between propagating fractures and pre-existing natural fractures. In this study, the TOUGH-RBSN code for coupled hydro-mechanical modeling is utilized to simulate hydraulic fracture propagation and its interaction with pre-existing fracture networks. The simulation tool combines TOUGH2, a simulator of subsurface multiphase flow and mass transport based on the finite volume approach, with the implementation of a lattice modeling approach for geomechanical and fracture-damage behavior, named Rigid-Body-Spring Network (RBSN). The discrete fracture network (DFN) approach is facilitated in the Voronoi discretization via a fully automated modeling procedure. The numerical program is verified through a simple simulation for single fracture propagation, in which the resulting fracture geometry is compared to an analytical solution for given fracture length and aperture. Subsequently, predictive simulations are conducted for planned laboratory experiments using rock-analogue (soda-lime glass) samples containing a designed, pre-existing fracture network. The results of a preliminary simulation demonstrate selective fracturing and fluid infiltration along the pre-existing fractures, with additional fracturing in part
Effects of using a continuum representation of discrete fracture networks
Hull, L.C.; Clemo, T.M.
1987-01-01
The substitution of matrix or continuum permeability for discrete fracture permeability in the simulation of complex fracture systems requires a radically different treatment of transport in the matrix. The spatial distribution of pressure is reasonably well described by inclusion of only the major fractures. Transport of tracer and heat, however, depends on a detailed knowledge of fluid velocities. Two factors are involved. First, the velocities are dependent on the active porosity of the system. Because fractures channel flow, the active porosity may be much smaller than the total porosity of the system. Secondly, the distribution of velocities is generally not normally distributed precluding the use of a Gaussian dispersion model. Characterization of the active porosity and velocity distribution are necessary to quantify tracer and heat movement.
Towards effective flow simulations in realistic discrete fracture networks
NASA Astrophysics Data System (ADS)
Berrone, Stefano; Pieraccini, Sandra; Scialò, Stefano
2016-04-01
We focus on the simulation of underground flow in fractured media, modeled by means of Discrete Fracture Networks. Focusing on a new recent numerical approach proposed by the authors for tackling the problem avoiding mesh generation problems, we further improve the new family of methods making a step further towards effective simulations of large, multi-scale, heterogeneous networks. Namely, we tackle the imposition of Dirichlet boundary conditions in weak form, in such a way that geometrical complexity of the DFN is not an issue; we effectively solve DFN problems with fracture transmissivities spanning many orders of magnitude and approaching zero; furthermore, we address several numerical issues for improving the numerical solution also in quite challenging networks.
Radionuclide migration analysis using a discrete fracture network model
Ijiri, Y.; Sawada, A.; Webb, E.K.; Watari, S.; Hatanaka, K.; Uchida, M.; Ishiguro, K.; Umeki, H.; Dershowitz, W.S.
1999-07-01
This paper describes an approach for assessing the geosphere performance of nuclear waste disposal in fractured rock. In this approach, a three-dimensional heterogeneous channel-network model is constructed using a stochastic discrete fracture network (DFN) code. Radionuclide migration in the channel-network model is solved using the Laplace transform Galerkin finite element method, taking into account advection-dispersion in a fracture network, matrix diffusion, sorption in the rock matrix as well as radioactive chain decay. Preliminary radionuclide migration analysis was performed for fifty realizations of a synthetic block-scale DFN model. The total radionuclide release from all packages in the repository was estimated from the statistics of the results of fifty realizations under the hypothesis of ergodicity. The interpretation of the result of the three-dimensional network model by a combination of simpler one-dimensional parallel plate models is also discussed.
Dershowitz, William S.; Einstein, Herbert H.; LaPoint, Paul R.; Eiben, Thorsten; Wadleigh, Eugene; Ivanova, Violeta
1998-12-01
This report summarizes research conducted for the Fractured Reservoir Discrete Feature Network Technologies Project. The five areas studied are development of hierarchical fracture models; fractured reservoir compartmentalization, block size, and tributary volume analysis; development and demonstration of fractured reservoir discrete feature data analysis tools; development of tools for data integration and reservoir simulation through application of discrete feature network technologies for tertiary oil production; quantitative evaluation of the economic value of this analysis approach.
NASA Astrophysics Data System (ADS)
Meheust, Y.; De Dreuzy, J.; Pichot, G.
2011-12-01
Flow channeling and permeability scaling in fractured media have been classically addressed either at the fracture- or at the network- scales. In the latter case they are linked to the topological structure of the network, while at the fracture scale they are controlled by the variability of the local aperture distribution inside individual fractures. In this study we analyze these two combined effects, investigating how flow localization below the scale of individual fractures influences that at the network scale and the resulting medium permeability. This is done by use of a new highly-resolved 3D discrete fracture network model (DFN). The local apertures of individual fractures are distributed according to a truncated Gaussian law, and exhibit self-affine spatial correlations that are bounded by an upper cutoff scale Lc; Lc and the fracture closure, defined as the ratio of the aperture fluctuations at scale Lc to the mean aperture, are considered homogeneous over the DFN. The network topology is controlled by a homogeneous scalar fracture density and a power law fracture length distribution. We have varied these features to investigate a large variety of DFN topologies, from sparse networks with varying degrees of fracture interconnections, flow bottlenecks and dead-ends (Fig. 1a), to dense well-connected networks (Fig. 1b). We have also investigated a large range of fracture closures, performing extensive simulations of about 105 different DFN realizations. At the fracture scale, accounting for local aperture fluctuations leads to a monotical deviation (which can exceed 50%) of the equivalent fracture transmissivity from the parallel plate behavior. At the network scale we observe a complex interaction between flow channeling within fracture planes and flow localization in the network. This interaction is controlled by the location of fracture interactions with respect to that of low local transmissivity zones (particularly the closed zones), in the fracture
Characterization of fracture processes by continuum and discrete modelling
NASA Astrophysics Data System (ADS)
Kaliske, M.; Dal, H.; Fleischhauer, R.; Jenkel, C.; Netzker, C.
2012-09-01
A large number of methods to describe fracture mechanical features of structures on basis of computational algorithms have been developed in the past due to the importance of the topic. In this paper, current and promising numerical approaches for the characterization of fracture processes are presented. A fracture phenomenon can either be depicted by a continuum formulation or a discrete notch. Thus, starting point of the description is a micromechanically motivated formulation for the development of a local failure situation. A current, generalized method without any restriction to material modelling and loading situation in order to describe an existing crack in a structure is available through the material force approach. One possible strategy to simulate arbitrary crack growth is based on an adaptive implementation of cohesive elements in combination with the standard discretization of the body. In this case, crack growth criteria and the determination of the crack propagation direction in combination with the modification of the finite element mesh are required. The nonlinear structural behaviour of a fibre reinforced composite material is based on the heterogeneous microstructure. A two-scale simulation is therefore an appropriate and effective way to take into account the scale differences of macroscopic structures with microscopic elements. In addition, fracture mechanical structural properties are far from being sharp and deterministic. Moreover, a wide range of uncertainties influence the ultimate load bearing behaviour. Therefore, it is evident that the deterministic modelling has to be expanded by a characterization of the uncertainty in order to achieve a reliable and realistic simulation result. The employed methods are illustrated by numerical examples.
Uncertainties in Parameterizing Faults for Discrete Fracture Network Models (Invited)
NASA Astrophysics Data System (ADS)
Doe, T.
2013-12-01
Discrete Fracture Network (DFN) models represent conducting faults and fractures as two dimensional features with transmissivity that defines the flow capacity, storativity that represents the open volume and compressibility, and aperture that determines the flow velocity. While these properties work well for capturing the pressure and flow behaviors they are not complete for understanding solute transport in groundwater or production in fault-dominated oil and gas reservoirs. For oil and gas production, water breakthrough will be controlled by the small portion of the fault that accounts for most of the fault's transmissivity, while the major portion in the fault's storage will be released more slowly. Similarly, in contaminant hydrology, the multiple porosities of the faults may enhance matrix diffusion effects. This presentation focuses on faults in granitic basement rock with low matrix permeability and porosity. Although these rocks contain joints and simple fractures, the major portion of the flow and storage may appear in faults and shear fractures that span a large range of scales. Unlike simple fractures or joints, faults and shear fractures have complex internal geometries including a range of materials including altered rock, breccia, mylonite, cataclasite, and gouge. These materials may have significant porosities ranging from a few percent to over thirty percent in heavily weathered or altered rocks. The core materials may act as either barriers or conduits. Altered granitic rock may also have porosities in the range of several percent. As faults may have thicknesses from a few centimeters to several tens of meters, the pore volumes in fault rock can be very significant. A key to understanding fault zone behaviors lies in having a range for the relative pore volumes of the conducting portion of the fault versus that of the larger damage zone. This presentation summarizes a literature review focusing on these fault properties. The conducting portion of
Kijima, Hiroaki; Konishi, Natsuo; Kubota, Hitoshi; Tazawa, Hiroshi; Tani, Takayuki; Suzuki, Norio; Kamo, Keiji; Okudera, Yoshihiko; Sasaki, Ken; Kawano, Tetsuya; Shimada, Yoichi
2014-01-01
The reliability of proximal femoral fracture classifications using 3DCT was evaluated, and a comprehensive “area classification” was developed. Eleven orthopedists (5–26 years from graduation) classified 27 proximal femoral fractures at one hospital from June 2013 to July 2014 based on preoperative images. Various classifications were compared to “area classification.” In “area classification,” the proximal femur is divided into 4 areas with 3 boundary lines: Line-1 is the center of the neck, Line-2 is the border between the neck and the trochanteric zone, and Line-3 links the inferior borders of the greater and lesser trochanters. A fracture only in the first area was classified as a pure first area fracture; one in the first and second area was classified as a 1-2 type fracture. In the same way, fractures were classified as pure 2, 3-4, 1-2-3, and so on. “Area classification” reliability was highest when orthopedists with varying experience classified proximal femoral fractures using 3DCT. Other classifications cannot classify proximal femoral fractures if they exceed each classification's particular zones. However, fractures that exceed the target zones are “dangerous” fractures. “Area classification” can classify such fractures, and it is therefore useful for selecting osteosynthesis methods. PMID:25610659
NASA Astrophysics Data System (ADS)
Donado-Garzon, L. D.; Pardo, Y.
2013-12-01
Fractured media are very heterogeneous systems where occur complex physical and chemical processes to model. One of the possible approaches to conceptualize this type of massifs is the Discrete Fracture Network (DFN). Donado et al., modeled flow and transport in a granitic batholith based on this approach and found good fitting with hydraulic and tracer tests, but the computational cost was excessive due to a gigantic amount of elements to model. We present in this work a methodology based on percolation theory for reducing the number of elements and in consequence, to reduce the bandwidth of the conductance matrix and the execution time of each network. DFN poses as an excellent representation of all the set of fractures of the media, but not all the fractures of the media are part of the conductive network. Percolation theory is used to identify which nodes or fractures are not conductive, based on the occupation probability or percolation threshold. In a fractured system, connectivity determines the flow pattern in the fractured rock mass. This volume of fluid is driven through connection paths formed by the fractures, when the permeability of the rock is negligible compared to the fractures. In a population of distributed fractures, each of this that has no intersection with any connected fracture do not contribute to generate a flow field. This algorithm also permits us to erase these elements however they are water conducting and hence, refine even more the backbone of the network. We used 100 different generations of DFN that were optimized in this study using percolation theory. In each of the networks calibrate hydrodynamic parameters as hydraulic conductivity and specific storage coefficient, for each of the five families of fractures, yielding a total of 10 parameters to estimate, at each generation. Since the effects of the distribution of fault orientation changes the value of the percolation threshold, but not the universal laws of classical
Management of mandibular angle fractures using a 1.7 mm 3-dimensional strut plate
Pandey, Varnika; Bhutia, Ongkila; Nagori, Shakil Ahmed; Seith, Ashu; Roychoudhury, Ajoy
2015-01-01
Aim We report our experience with the use of 1.7 mm 3-dimentional (3D) strut plate for the management of mandibular angle fractures. Methods This prospective study enrolled 15 patients in whom mandibular angle fractures were treated with 1.7 mm 3D plate using trans-buccal trochar. Patients were evaluated at 72 h, 2 weeks, 6 weeks and 12 weeks for fracture stability, occlusion, soft-tissue swelling, infection and post-operative inferior alveolar nerve damage. Other complications like wound dehiscence, non-union, mal-union and hardware failure were also assessed. Results In the immediate post-operative period, fracture instability was seen in 1 (6.7%) patient which resolved by 2 weeks. Mild occlusal discrepancy was also noted in 1 (6.7%) patient. Wound dehiscence was seen in 5 (33.3%) patients and all resolved by local measures. 1 (6.7%) patient developed post-operative nerve paraesthesia. Immediate post-operative radiographic evaluation demonstrated optimal reduction in all cases with no inferior border gaping. No case of infection, hardware failure, non-union and mal-union was noted. Conclusion Within the limitations of the study, 1.7 mm 3D strut plate was found to be effective for management of non-communited mandibular angle fractures. PMID:26937367
He, Jiliang; Tan, Guoqing; Zhou, Dongsheng; Sun, Liang; Li, Qinghu; Yang, Yongliang; Liu, Ping
2016-01-01
Abstract Percutaneous screw insertion for minimally displaced or reducible acetabular fracture using x-ray fluoroscopy and computer-assisted navigation system has been advocated by some authors. The purpose of this study was to compare intraoperative conditions and clinical results between isocentric C-arm 3-dimensional (Iso-C 3D) fluoroscopy and conventional fluoroscopy for percutaneous retrograde screwing of acetabular anterior column fracture. A prospective cohort study was conducted. A total of 22 patients were assigned to 2 different groups: 10 patients in the Iso-C 3D navigation group and 12 patients in the conventional group. The operative time, fluoroscopic time, time of screw insertion, blood loss, and accuracy were analyzed between the 2 groups. There were significant differences in operative time, screw insertion time, fluoroscopy time, and mean blood loss between the 2 groups. Totally 2 of 12 (16.7%) screws were misplaced in the conventional fluoroscopy group, and all 10 screws were in safe zones in the navigation group. Percutaneous screw fixation using the Iso-C 3D computer-assisted navigation system significantly reduced the intraoperative fluoroscopy time and blood loss in percutaneous screwing for acetabular anterior column fracture. The Iso-C 3D computer-assisted navigation system provided a reliable and effective method for percutaneous screw insertion in acetabular anterior column fractures compared to conventional fluoroscopy. PMID:26765448
He, Jiliang; Tan, Guoqing; Zhou, Dongsheng; Sun, Liang; Li, Qinghu; Yang, Yongliang; Liu, Ping
2016-01-01
Percutaneous screw insertion for minimally displaced or reducible acetabular fracture using x-ray fluoroscopy and computer-assisted navigation system has been advocated by some authors. The purpose of this study was to compare intraoperative conditions and clinical results between isocentric C-arm 3-dimensional (Iso-C 3D) fluoroscopy and conventional fluoroscopy for percutaneous retrograde screwing of acetabular anterior column fracture.A prospective cohort study was conducted. A total of 22 patients were assigned to 2 different groups: 10 patients in the Iso-C 3D navigation group and 12 patients in the conventional group. The operative time, fluoroscopic time, time of screw insertion, blood loss, and accuracy were analyzed between the 2 groups.There were significant differences in operative time, screw insertion time, fluoroscopy time, and mean blood loss between the 2 groups. Totally 2 of 12 (16.7%) screws were misplaced in the conventional fluoroscopy group, and all 10 screws were in safe zones in the navigation group. Percutaneous screw fixation using the Iso-C 3D computer-assisted navigation system significantly reduced the intraoperative fluoroscopy time and blood loss in percutaneous screwing for acetabular anterior column fracture.The Iso-C 3D computer-assisted navigation system provided a reliable and effective method for percutaneous screw insertion in acetabular anterior column fractures compared to conventional fluoroscopy. PMID:26765448
Huang, Hai; Plummer, Mitchell; Podgorney, Robert
2013-02-01
Advancement of EGS requires improved prediction of fracture development and growth during reservoir stimulation and long-term operation. This, in turn, requires better understanding of the dynamics of the strongly coupled thermo-hydro-mechanical (THM) processes within fractured rocks. We have developed a physically based rock deformation and fracture propagation simulator by using a quasi-static discrete element model (DEM) to model mechanical rock deformation and fracture propagation induced by thermal stress and fluid pressure changes. We also developed a network model to simulate fluid flow and heat transport in both fractures and porous rock. In this paper, we describe results of simulations in which the DEM model and network flow & heat transport model are coupled together to provide realistic simulation of the changes of apertures and permeability of fractures and fracture networks induced by thermal cooling and fluid pressure changes within fractures. Various processes, such as Stokes flow in low velocity pores, convection-dominated heat transport in fractures, heat exchange between fluid-filled fractures and solid rock, heat conduction through low-permeability matrices and associated mechanical deformations are all incorporated into the coupled model. The effects of confining stresses, developing thermal stress and injection pressure on the permeability evolution of fracture and fracture networks are systematically investigated. Results are summarized in terms of implications for the development and evolution of fracture distribution during hydrofracturing and thermal stimulation for EGS.
An implicit finite element method for discrete dynamic fracture
Jobie M. Gerken
1999-12-01
A method for modeling the discrete fracture of two-dimensional linear elastic structures with a distribution of small cracks subject to dynamic conditions has been developed. The foundation for this numerical model is a plane element formulated from the Hu-Washizu energy principle. The distribution of small cracks is incorporated into the numerical model by including a small crack at each element interface. The additional strain field in an element adjacent to this crack is treated as an externally applied strain field in the Hu-Washizu energy principle. The resulting stiffness matrix is that of a standard plane element. The resulting load vector is that of a standard plane element with an additional term that includes the externally applied strain field. Except for the crack strain field equations, all terms of the stiffness matrix and load vector are integrated symbolically in Maple V so that fully integrated plane stress and plane strain elements are constructed. The crack strain field equations are integrated numerically. The modeling of dynamic behavior of simple structures was demonstrated within acceptable engineering accuracy. In the model of axial and transverse vibration of a beam and the breathing mode of vibration of a thin ring, the dynamic characteristics were shown to be within expected limits. The models dominated by tensile forces (the axially loaded beam and the pressurized ring) were within 0.5% of the theoretical values while the shear dominated model (the transversely loaded beam) is within 5% of the calculated theoretical value. The constant strain field of the tensile problems can be modeled exactly by the numerical model. The numerical results should therefore, be exact. The discrepancies can be accounted for by errors in the calculation of frequency from the numerical results. The linear strain field of the transverse model must be modeled by a series of constant strain elements. This is an approximation to the true strain field, so some
Discrete fracture simulations of the hydrogeology at Koongarra, Northern Territory, Australia
Smoot, J.L.
1992-04-01
The US Department of Energy is studying the Alligator Rivers Natural Analogue Project site at Koongarra, Northern Territory, Australia to investigate and simulate radionuclide migration in fractured rocks. Discrete fracture simulations were conducted within a cubic volume (180-m edge length) of fractured Cahill Formation schist oriented with one major axis parallel to the trend of the Koongarra Fault. Five hundred fractures are simulated within this domain. The fractures have a mean orientation parallel to the idealized plane of the Koongarra Fault dipping 55{degrees} SE. Simple flow modeling of this fracture network was conducted by assigning constant head boundaries to upgradient and downgradient vertical faces of the cube, which trend parallel to the fault. No-flow boundaries were assigned to all other faces. The fracture network allows hydraulic communication across the block, in spite of relatively low fracture density across the block.
A new computer code for discrete fracture network modelling
NASA Astrophysics Data System (ADS)
Xu, Chaoshui; Dowd, Peter
2010-03-01
The authors describe a comprehensive software package for two- and three-dimensional stochastic rock fracture simulation using marked point processes. Fracture locations can be modelled by a Poisson, a non-homogeneous, a cluster or a Cox point process; fracture geometries and properties are modelled by their respective probability distributions. Virtual sampling tools such as plane, window and scanline sampling are included in the software together with a comprehensive set of statistical tools including histogram analysis, probability plots, rose diagrams and hemispherical projections. The paper describes in detail the theoretical basis of the implementation and provides a case study in rock fracture modelling to demonstrate the application of the software.
NASA Astrophysics Data System (ADS)
Makedonska, Nataliia; Hyman, Jeffrey D.; Karra, Satish; Painter, Scott L.; Gable, Carl W.; Viswanathan, Hari S.
2016-08-01
The apertures of natural fractures in fractured rock are highly heterogeneous. However, in-fracture aperture variability is often neglected in flow and transport modeling and individual fractures are assumed to have uniform aperture distribution. The relative importance of in-fracture variability in flow and transport modeling within kilometer-scale field-scale fracture networks has been under a matter of debate for a long time because the flow in each single fracture is controlled not only by in-fracture variability but also by boundary conditions. Computational limitations have previously prohibited researchers from investigating the relative importance of in-fracture variability in flow and transport modeling within large-scale fracture networks. We address this question by incorporating internal heterogeneity of individual fractures into flow simulations within kilometer scale three-dimensional fracture networks, where fracture intensity, P32 (ratio between total fracture area and domain volume) is between 0.027 and 0.031 [1/m]. A recently developed discrete fracture network (DFN) simulation capability, dfnWorks, is used to generate DFNs that include in-fracture aperture variability represented by a stationary log-normal stochastic field with various correlation lengths and variances. The Lagrangian transport parameters, non-reacting travel time and cumulative retention, are calculated along particles streamlines. It is observed that due to local flow channeling early particle travel times are more sensitive to in-fracture variability than the tails of travel time distributions, where no significant effect of the in-fracture transmissivity variations and spatial correlation length is observed.
Makedonska, Nataliia; Hyman, Jeffrey D.; Karra, Satish; Painter, Scott L.; Gable, Carl W.; Viswanathan, Hari S.
2016-06-17
The apertures of natural fractures in fractured rock are highly heterogeneous. However, in-fracture aperture variability is often neglected in flow and transport modeling and individual fractures are assumed to have uniform aperture distribution. The relative importance of in-fracture variability in flow and transport modeling within kilometer18 scale field–scale fracture networks has been under a matter of debate for a long time because the flow in each single fracture is controlled not only by in-fracture variability but also by boundary conditions. Computational limitations have previously prohibited researchers from investigating the relative importance of in-fracture variability in flow and transport modeling withinmore » large-scale fracture networks. We address this question by incorporating internal heterogeneity of individual fractures into 23 flow simulations within kilometer scale three-dimensional fracture networks, where fracture intensity, P32 (ratio between total fracture area and domain volume) is between 0.027 and 0.031 [1/m]. A recently developed discrete fracture network (DFN) simulation capability, dfnWorks, is used to generate DFNs that include in-fracture aperture variability represented by a stationary log-normal stochastic field with various correlation lengths and variances. The Lagrangian transport parameters, non-reacting travel time and cumulative retention, are calculated along particles streamlines. It is observed that due to local flow channeling early particle travel times are more sensitive to in-fracture variability than the tails of travel time distributions, where no significant effect of the in-fracture transmissivity variations and spatial correlation length is observed.« less
Zhou, Q.; Salve, R.; Liu, H.-H.; Wang, J.S.Y.; Hudson, D.
2006-01-01
A mesoscale (21??m in flow distance) infiltration and seepage test was recently conducted in a deep, unsaturated fractured rock system at the crossover point of two underground tunnels. Water was released from a 3??m ?? 4??m infiltration plot on the floor of an alcove in the upper tunnel, and seepage was collected from the ceiling of a niche in the lower tunnel. Significant temporal and (particularly) spatial variabilities were observed in both measured infiltration and seepage rates. To analyze the test results, a three-dimensional unsaturated flow model was used. A column-based scheme was developed to capture heterogeneous hydraulic properties reflected by these spatial variabilities observed. Fracture permeability and van Genuchten ?? parameter [van Genuchten, M.T., 1980. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci. Soc. Am. J. 44, 892-898] were calibrated for each rock column in the upper and lower hydrogeologic units in the test bed. The calibrated fracture properties for the infiltration and seepage zone enabled a good match between simulated and measured (spatially varying) seepage rates. The numerical model was also able to capture the general trend of the highly transient seepage processes through a discrete fracture network. The calibrated properties and measured infiltration/seepage rates were further compared with mapped discrete fracture patterns at the top and bottom boundaries. The measured infiltration rates and calibrated fracture permeability of the upper unit were found to be partially controlled by the fracture patterns on the infiltration plot (as indicated by their positive correlations with fracture density). However, no correlation could be established between measured seepage rates and density of fractures mapped on the niche ceiling. This lack of correlation indicates the complexity of (preferential) unsaturated flow within the discrete fracture network. This also indicates that continuum
Towards improved 3D cross-borehole electrical resistivity imaging of discrete fracture networks
NASA Astrophysics Data System (ADS)
Robinson, J.; Slater, L. D.; Johnson, T. J.; Ntarlagiannis, D.; Lacombe, P.; Johnson, C. D.; Tiedeman, C. R.; Goode, D.; Day-Lewis, F. D.; Shapiro, A. M.; Lane, J. W.
2012-12-01
There is a need to better characterize discrete fractures in contaminated bedrock aquifers to determine the migration of injected remediation amendments away from boreholes. A synthetic cross-borehole electrical resistivity study was conducted assuming a discrete fracture model of an existing contaminated site with known fracture locations. Four boreholes and two discrete fracture zones, assumed to be the dominant electrical and hydraulically conductive pathways, were explicitly modeled within an unstructured tetrahedral finite-element mesh. To simulate field conditions, 5% random Gaussian noise was added to all synthetic datasets. We first evaluated different regularization constraints starting with an uninformed smoothness-constrained inversion, to which a priori information was incrementally added. We found major improvements when (1) smoothness regularization constraints were relaxed (or disconnected) along boreholes and fractures, (2) a homogeneous conductivity was assumed along boreholes, and (3) borehole conductivity constraints, which could be determined from a fluid specific-conductance log, were applied. We also evaluated the effect of including borehole packers on the fracture-zone model recovery. We found the estimated fracture-zone conductivities with the inclusion of packers were comparable to similar trials excluding the use of packers regardless of electrical potential changes. The misplacement of fracture regularization disconnects easily can be misinterpreted as actual fracture locations. Conductivities within misplaced disconnects were near the starting model value and removing smoothing between boreholes and assumed fracture locations helped in identifying incorrectly located fracture regularization disconnects. Model sensitivity structure improved when regularization disconnects were (1) applied along the boreholes and fracture zones, and (2) fracture-zone regularization disconnects were placed where actual fractures existed. A field study
Combined Finite-Discrete Element Method for Simulation of Hydraulic Fracturing
NASA Astrophysics Data System (ADS)
Yan, Chengzeng; Zheng, Hong; Sun, Guanhua; Ge, Xiurun
2016-04-01
Hydraulic fracturing is widely used in the exploitation of unconventional gas (such as shale gas).Thus, the study of hydraulic fracturing is of particular importance for petroleum industry. The combined finite-discrete element method (FDEM) proposed by Munjiza is an innovative numerical technique to capture progressive damage and failure processes in rock. However, it cannot model the fracturing process of rock driven by hydraulic pressure. In this study, we present a coupled hydro-mechanical model based on FDEM for the simulation of hydraulic fracturing in complex fracture geometries, where an algorithm for updating hydraulic fracture network is proposed. The algorithm can carry out connectivity searches for arbitrarily complex fracture networks. Then, we develop a new combined finite-discrete element method numerical code (Y-flow) for the simulation of hydraulic fracturing. Finally, several verification examples are given, and the simulation results agree well with the analytical or experimental results, indicating that the newly developed numerical code can capture hydraulic fracturing process correctly and effectively.
Incorporating Discrete Irregular Fracture Zone Networks into 3D Paleohydrogeologic Simulations
NASA Astrophysics Data System (ADS)
Normani, S. D.
2015-12-01
Dual continuum computational models which include both porous media and discrete fracture zones are valuable tools in assessing groundwater migration and pathways in fractured rock systems. Fracture generation models can produce stochastic realizations of fracture networks which honor geological structures and fracture propagation behaviors. Surface lineament traces can be propagated to depth based on fracture zone statistics to produce representations of geological structures in rock. The generated discrete, complex and irregular fracture zone networks, represented as a triangulated mesh, are embedded using orthogonal quadrilateral elements within a three-dimensional hexahedral finite element mesh. A detailed coupled density-dependent paleohydrogeologic groundwater analysis of a hypothetical 104 km2 portion of the Canadian Shield has been conducted using the discrete-fracture dual continuum finite element model FRAC3DVS to investigate the characterization of large-scale fracture zone networks on groundwater and tracer movement during a 120,000 year paleoclimate cycle. Permeability reduction due to permafrost was also applied. Time series data for the depth of permafrost, along with ice thickness and lake depth, were provided by the University of Toronto (UofT) Glacial Systems Model. The crystalline rock between fracture zones was assigned properties characteristic of those reported for the Canadian Shield. Total dissolved solids concentrations of 300 g/L are encountered at depth. Surface water features and a Digital Elevation Model (DEM) were used in a GIS framework to define the watershed boundaries at surface water divides and to populate the finite element mesh. This work will illustrate the long-term evolution and stability of the geosphere and groundwater systems to external perturbations caused by glaciation through the use of performance measures such as Mean Life Expectancy and the migration of a unit tracer to depth over a paleoclimate cycle.
NASA Astrophysics Data System (ADS)
Zhang, Qi-Hua
2015-10-01
Finite element generation of complicated fracture networks is the core issue and source of technical difficulty in three-dimensional (3-D) discrete fracture network (DFN) flow models. Due to the randomness and uncertainty in the configuration of a DFN, the intersection lines (traces) are arbitrarily distributed in each face (fracture and other surfaces). Hence, subdivision of the fractures is an issue relating to subdivision of two-dimensional (2-D) domains with arbitrarily-distributed constraints. When the DFN configuration is very complicated, the well-known approaches (e.g. Voronoi Delaunay-based methods and advancing-front techniques) cannot operate properly. This paper proposes an algorithm to implement end-to-end connection between traces to subdivide 2-D domains into closed loops. The compositions of the vertices in the common edges between adjacent loops (which may belong to a single fracture or two connected fractures) are thus ensured to be topologically identical. The paper then proposes an approach for triangulating arbitrary loops which does not add any nodes to ensure consistency of the meshes at the common edges. In addition, several techniques relating to tolerance control and improving code robustness are discussed. Finally, the equivalent permeability of the rock mass is calculated for some very complicated DFNs (the DFN may contain 1272 fractures, 633 connected fractures, and 16,270 closed loops). The results are compared with other approaches to demonstrate the veracity and efficiency of the approach proposed in this paper.
Discrete Dual Porosity Modeling of Electrical Current Flow in Fractured Media
NASA Astrophysics Data System (ADS)
Roubinet, D.; Irving, J.
2013-12-01
The study of fractured rocks is highly important in a variety of research fields and applications such as hydrogeology, geothermal energy, hydrocarbon extraction, and the long-term storage of toxic waste. Fractured media are characterized by a large contrast in permeability between the fractures and the rock matrix. For hydrocarbon extraction, the presence of highly conductive fractures is an advantage as they allow for quick and easy access to the resource. For toxic waste storage, however, the fractures represent a significant drawback as there is an increased risk of leakage and migration of pollutants deep into the subsurface. In both cases, the identification of fracture network characteristics is a critical, challenging, and required step. A number of previous studies have indicated that the presence of fractures in geological materials can have a significant impact on geophysical electrical resistivity measurements. It thus appears that, in some cases, geoelectrical surveys might be used to obtain useful information regarding fracture network characteristics. However, existing geoelectrical modeling tools and inversion methods are not properly adapted to deal with the specific challenges of fractured media. This prevents us from fully exploring the potential of the method to characterize fracture network properties. We thus require, as a first step, the development of accurate and efficient numerical modeling tools specifically designed for fractured domains. Building on the discrete fracture network (DFN) approach that has been widely used for modeling groundwater flow in fractured rocks, we have developed a discrete dual-porosity model for electrical current flow in fractured media. Our novel approach combines an explicit representation of the fractures with fracture-matrix electrical flow exchange at the block-scale. Tests in two dimensions show the ability of our method to deal with highly heterogeneous fracture networks in a highly computationally
NASA Astrophysics Data System (ADS)
Ryerson, F. J.; Ezzedine, S. M.; Glascoe, L. G.; Antoun, T. H.
2011-12-01
Fractures and fracture networks are the principle pathways for migration of water, heat and mass in enhanced geothermal systems, oil and gas reservoirs, CO2 leakage from saline aquifers, and radioactive and toxic industrial wastes from underground storage repositories. A major issue to overcome when characterizing a fractured reservoir is that of data limitation due to accessibility and affordability. Moreover, the ability to map discontinuities in the rock with available geological and geophysical tools tends to decrease particularly as the scale of the discontinuity goes down. Data collected are often reduced to probability distribution functions for predictive modeling and simulation in a stochastic framework such as stochastic discrete fracture network. Stochastic discrete fracture network models enable probabilistic assessment of flow, transport and geomechanical phenomena that are not adequately captured using continuum models. Despite the fundamental uncertainties inherited within the probabilistic reduction of the sparse data collected, very little work has been conducted on quantifying uncertainty on the reduced probabilistic distribution functions. In the current study, we investigate the impact of parameter uncertainties of the distribution functions that characterize discrete fracture networks on the flow, heat and mass transport and geomechanics. Numerical results of first, second and third moments, normalized to a base case scenario, are presented and compared to theoretical results extended from percolation theory. (Prepared by LLNL under Contract DE-AC52-07NA27344)
Dershowitz, W.S.; La Pointe, P.R.; Einstein, H.H.; Ivanova, V.
1998-01-01
This report describes progress on the project, {open_quotes}Fractured Reservoir Discrete Feature Network Technologies{close_quotes} during the period March 7, 1996 to February 28, 1997. The report presents summaries of technology development for the following research areas: (1) development of hierarchical fracture models, (2) fractured reservoir compartmentalization and tributary volume, (3) fractured reservoir data analysis, and (4) integration of fractured reservoir data and production technologies. In addition, the report provides information on project status, publications submitted, data collection activities, and technology transfer through the world wide web (WWW). Research on hierarchical fracture models included geological, mathematical, and computer code development. The project built a foundation of quantitative, geological and geometrical information about the regional geology of the Permian Basin, including detailed information on the lithology, stratigraphy, and fracturing of Permian rocks in the project study area (Tracts 17 and 49 in the Yates field). Based on the accumulated knowledge of regional and local geology, project team members started the interpretation of fracture genesis mechanisms and the conceptual modeling of the fracture system in the study area. Research on fractured reservoir compartmentalization included basic research, technology development, and application of compartmentalized reservoir analyses for the project study site. Procedures were developed to analyze compartmentalization, tributary drainage volume, and reservoir matrix block size. These algorithms were implemented as a Windows 95 compartmentalization code, FraCluster.
Prediction of Fracture Behavior in Rock and Rock-like Materials Using Discrete Element Models
NASA Astrophysics Data System (ADS)
Katsaga, T.; Young, P.
2009-05-01
The study of fracture initiation and propagation in heterogeneous materials such as rock and rock-like materials are of principal interest in the field of rock mechanics and rock engineering. It is crucial to study and investigate failure prediction and safety measures in civil and mining structures. Our work offers a practical approach to predict fracture behaviour using discrete element models. In this approach, the microstructures of materials are presented through the combination of clusters of bonded particles with different inter-cluster particle and bond properties, and intra-cluster bond properties. The geometry of clusters is transferred from information available from thin sections, computed tomography (CT) images and other visual presentation of the modeled material using customized AutoCAD built-in dialog- based Visual Basic Application. Exact microstructures of the tested sample, including fractures, faults, inclusions and void spaces can be duplicated in the discrete element models. Although the microstructural fabrics of rocks and rock-like structures may have different scale, fracture formation and propagation through these materials are alike and will follow similar mechanics. Synthetic material provides an excellent condition for validating the modelling approaches, as fracture behaviours are known with the well-defined composite's properties. Calibration of the macro-properties of matrix material and inclusions (aggregates), were followed with the overall mechanical material responses calibration by adjusting the interfacial properties. The discrete element model predicted similar fracture propagation features and path as that of the real sample material. The path of the fractures and matrix-inclusion interaction was compared using computed tomography images. Initiation and fracture formation in the model and real material were compared using Acoustic Emission data. Analysing the temporal and spatial evolution of AE events, collected during the
Shouchun Deng; Robert Podgorney; Hai Huang
2011-02-01
Key challenges associated with the EGS reservoir development include the ability to reliably predict hydraulic fracturing and the deformation of natural fractures as well as estimating permeability evolution of the fracture network with time. We have developed a physics-based rock deformation and fracture propagation simulator by coupling a discrete element model (DEM) for fracturing with a network flow model. In DEM model, solid rock is represented by a network of discrete elements (often referred as particles) connected by various types of mechanical bonds such as springs, elastic beams or bonds that have more complex properties (such as stress-dependent elastic constants). Fracturing is represented explicitly as broken bonds (microcracks), which form and coalesce into macroscopic fractures when external and internal load is applied. The natural fractures are represented by a series of connected line segments. Mechanical bonds that intersect with such line segments are removed from the DEM model. A network flow model using conjugate lattice to the DEM network is developed and coupled with the DEM. The fluid pressure gradient exerts forces on individual elements of the DEM network, which therefore deforms the mechanical bonds and breaks them if the deformation reaches a prescribed threshold value. Such deformation/fracturing in turn changes the permeability of the flow network, which again changes the evolution of fluid pressure, intimately coupling the two processes. The intimate coupling between fracturing/deformation of fracture networks and fluid flow makes the meso-scale DEM- network flow simulations necessary in order to accurately evaluate the permeability evolution, as these methods have substantial advantages over conventional continuum mechanical models of elastic rock deformation. The challenges that must be overcome to simulate EGS reservoir stimulation, preliminary results, progress to date and near future research directions and opportunities will be
Anna, L.O.
1998-09-01
A three-dimensional discrete fracture model was completed to investigate the potential effects of fractures on the flow of water at Yucca Mountain, Nye County, Nevada. A fracture network of the Exploratory Studies Facility starter tunnel area was simulated and calibrated with field data. Two modeled volumes were used to simulate three-dimensional fracture networks of the Tiva Canyon tuff. One volume had a width and length of 150 meters, and the other had a width and length of 200 meters; both volumes were 60 meters thick. The analysis shows that the fracture system in the Exploratory Studies Facility starter tunnel area has numerous connected fractures that have relatively large permeabilities. However, pathway analysis between three radial boreholes indicated there were few pathways and little connection, which is consistent with results of cross-boreholes pressure testing. Pathway analysis also showed that at the scales used there was only one pathway connecting one end of the flow box to the opposite end. The usual vertical pathway was along one large fracture, whereas in four horizontal directions the pathway was from multiple fracture connections. As a result, the fracture network can be considered sparse. The fracture network was refined by eliminating nonconductive fractures determined from field-derived permeabilities. Small fractures were truncated from the simulated network without any effect on the overall connectivity. Fractures as long as 1.25 meters were eliminated (a large percentage of the total number of fractures) from the network without altering the number of pathways. Five directional permeabilities were computed for the 150- and 200-meter-scale flow box areas. Permeabilities for the 150-meter scale vary by almost two orders of magnitude, with the principal permeability direction being easterly. At the 200-meter scale, however, the flow box permeabilities only vary by a factor of four, with the principal permeability direction being vertical.
Discrete fracture patterns of virus shells reveal mechanical building blocks.
Ivanovska, Irena L; Miranda, Roberto; Carrascosa, Jose L; Wuite, Gijs J L; Schmidt, Christoph F
2011-08-01
Viral shells are self-assembled protein nanocontainers with remarkable material properties. They combine simplicity of construction with toughness and complex functionality. These properties make them interesting for bionanotechnology. To date we know little about how virus structure determines assembly pathways and shell mechanics. We have here used atomic force microscopy to study structural failure of the shells of the bacteriophage Φ29. We observed rigidity patterns following the symmetry of the capsid proteins. Under prolonged force exertion, we observed fracture along well-defined lines of the 2D crystal lattice. The mechanically most stable building block of the shells was a trimer. Our approach of "reverse engineering" the virus shells thus made it possible to identify stable structural intermediates. Such stable intermediates point to a hierarchy of interactions among equal building blocks correlated with distinct next-neighbor interactions. The results also demonstrate that concepts from macroscopic materials science, such as fracture, can be usefully employed in molecular engineering. PMID:21768340
Hyman, Jeffrey De'Haven; Aldrich, Garrett Allen; Viswanathan, Hari S.; Makedonska, Nataliia; Karra, Satish
2016-08-25
We characterize how different fracture size-transmissivity relationships influence flow and transport simulations through sparse three-dimensional discrete fracture networks. Although it is generally accepted that there is a positive correlation between a fracture's size and its transmissivity/aperture, the functional form of that relationship remains a matter of debate. Relationships that assume perfect correlation, semicorrelation, and noncorrelation between the two have been proposed. To study the impact that adopting one of these relationships has on transport properties, we generate multiple sparse fracture networks composed of circular fractures whose radii follow a truncated power law distribution. The distribution of transmissivities are selected somore » that the mean transmissivity of the fracture networks are the same and the distributions of aperture and transmissivity in models that include a stochastic term are also the same. We observe that adopting a correlation between a fracture size and its transmissivity leads to earlier breakthrough times and higher effective permeability when compared to networks where no correlation is used. While fracture network geometry plays the principal role in determining where transport occurs within the network, the relationship between size and transmissivity controls the flow speed. Lastly, these observations indicate DFN modelers should be aware that breakthrough times and effective permeabilities can be strongly influenced by such a relationship in addition to fracture and network statistics.« less
Discrete fracture modeling of hydro-mechanical damage processes in geological systems
NASA Astrophysics Data System (ADS)
Kim, K.; Rutqvist, J.; Houseworth, J. E.; Birkholzer, J. T.
2014-12-01
This study presents a modeling approach for investigating coupled thermal-hydrological-mechanical (THM) behavior, including fracture development, within geomaterials and structures. In the model, the coupling procedure consists of an effective linkage between two codes: TOUGH2, a simulator of subsurface multiphase flow and mass transport based on the finite volume approach; and an implementation of the rigid-body-spring network (RBSN) method, a discrete (lattice) modeling approach to represent geomechanical behavior. One main advantage of linking these two codes is that they share the same geometrical mesh structure based on the Voronoi discretization, so that a straightforward representation of discrete fracture networks (DFN) is available for fluid flow processes. The capabilities of the TOUGH-RBSN model are demonstrated through simulations of hydraulic fracturing, where fluid pressure-induced fracturing and damage-assisted flow are well represented. The TOUGH-RBSN modeling methodology has been extended to enable treatment of geomaterials exhibiting anisotropic characteristics. In the RBSN approach, elastic spring coefficients and strength parameters are systematically formulated based on the principal bedding direction, which facilitate a straightforward representation of anisotropy. Uniaxial compression tests are simulated for a transversely isotropic material to validate the new modeling scheme. The model is also used to simulate excavation fracture damage for the HG-A microtunnel in the Opalinus Clay rock, located at the Mont Terri underground research laboratory (URL) near Saint-Ursanne, Switzerland. The Opalinus Clay has transversely isotropic material properties caused by natural features such as bedding, foliation, and flow structures. Preferential fracturing and tunnel breakouts were observed following excavation, which are believed to be strongly influenced by the mechanical anisotropy of the rock material. The simulation results are qualitatively
NASA Astrophysics Data System (ADS)
Ezzedine, S. M.
2010-12-01
Fractures and fracture networks are the principle pathways for migration of water, heat and mass in enhanced geothermal systems, oil and gas reservoirs, CO2 leakage from saline aquifers, and radioactive and toxic industrial wastes from underground storage repositories. A major issue to overcome when characterizing a fractured reservoir is that of data limitation due to accessibility and affordability. Moreover, the ability to map discontinuities in the rock with available geological and geophysical tools tends to decrease particularly as the scale of the discontinuity goes down. Geological characterization data include measurements of fracture density, orientation, extent, and aperture, and are based on analysis of outcrops, borehole optical and acoustic televiewer logs, aerial photographs, and core samples among others. All of these measurements are taken at the field scale through a very sparse limited number of deep boreholes. These types of data are often reduced to probability distributions function for predictive modeling and simulation in a stochastic framework such as stochastic discrete fracture network. Stochastic discrete fracture network models enable, through Monte Carlo realizations and simulations, for probabilistic assessment of flow and transport phenomena that are not adequately captured using continuum models. Despite the fundamental uncertainties inherited within the probabilistic reduction of the sparse data collected, very little work has been conducted on quantifying uncertainty on the reduced probabilistic distribution functions. In the current study, using nested Monte Carlo simulations, we present the impact of parameter uncertainties of the distribution functions that characterize discrete fracture networks on the flow, heat and mass transport. Numerical results of first, second and third moments, normalized to a base case scenario, are presented and compared to theoretical results extended from percolation theory.
Colloid facilitated transport of lanthanides through discrete fractures in chalk
NASA Astrophysics Data System (ADS)
Tran, Emily; Klein Ben-David, Ofra; Teutsch, Nadya; Weisbrod, Noam
2015-04-01
Geological disposal of high-level radioactive waste is the internationally agreed-upon, long term solution for the disposal of long lived radionuclides and spent fuel. Eventually, corrosion of the waste canisters may lead to leakage of their hazardous contents, and the radionuclides can ultimately make their way into groundwater and pose a threat to the biosphere. Engineered bentonite barriers placed around nuclear waste repositories are generally considered sufficient to impede the transport of radionuclides from their storage location to the groundwater. However, colloidal-sized mobile bentonite particles eroding from these barriers have come under investigation as a potential transport vector for radionuclides sorbed to them. In addition, the presence of organic matter in groundwater has been shown to additionally facilitate the uptake of radionuclides by the clay colloids. This study aims to evaluate the transport behaviors of radionuclides in colloid-facilitated transport through a fractured chalk matrix and under geochemical conditions representative of the Negev desert, Israel. Lanthanides are considered an acceptable substitute to actinides for research on radionuclide transportation due to their similar chemical behavior. In this study, the migration of Ce both with and without colloidal particles was explored and compared to the migration of a conservative tracer (bromide). Tracer solutions containing known concentrations of Ce, bentonite colloids, humic acid and bromide were prepared in a matrix solution containing salt concentrations representative of that of the average rain water found in the Negev. These solutions were then injected into a flow system constructed around a naturally fractured chalk core. Samples were analyzed for Ce and Br using ICP-MS, and colloid concentrations were determined using spectrophotographic analysis. Breakthrough curves comparing the rates of transportation of each tracer were obtained, allowing for comparison of
Development of discrete flow paths in unsaturated fractures at Yucca Mountain.
Bodvarsson, G S; Wu, Yu-Shu; Zhang, Keni
2003-01-01
We have carried out numerical modeling studies to investigate the development of discrete-fracture flow paths and flow-focusing phenomena in the unsaturated rock of the potential repository horizon at Yucca Mountain, Nevada. These studies are based on two- and three-dimensional (2-D and 3-D) numerical models using site-specific parameters. The 2-D and 3-D models use high-resolution spatial discretization to explicitly include effects of discrete fractures with stochastically developed fracture permeabilities and a continuum approach. The permeability field is generated based on air permeability measurements at various scales. For most of the cases considered, uniform infiltration with different average rates (1-500 mm/year) is prescribed at the top of the model, while variability in outflow at the bottom of the model is used to evaluate the degree of flow focusing. In addition, scenarios involving nonuniform infiltration at the top boundary, different permeability correlation lengths and different flow-allocation schemes were analyzed. The modeling results obtained from all of the cases showed a remarkably similar flow-focusing pattern at the repository horizon. Furthermore, tracer transport simulation results also revealed additional features of focused flow and transport through the fracture network. PMID:12714283
NASA Astrophysics Data System (ADS)
Ezzedine, S. M.
2009-12-01
Fractures and fracture networks are the principal pathways for transport of water and contaminants in groundwater systems, enhanced geothermal system fluids, migration of oil and gas, carbon dioxide leakage from carbon sequestration sites, and of radioactive and toxic industrial wastes from underground storage repositories. A major issue to overcome when characterizing a fractured reservoir is that of data limitation due to accessibility and affordability. Moreover, the ability to map discontinuities in the rock with available geological and geophysical tools tends to decrease particularly as the scale of the discontinuity goes down. Geological characterization data include measurements of fracture density, orientation, extent, and aperture, and are based on analysis of outcrops, borehole optical and acoustic televiewer logs, aerial photographs, and core samples, among other techniques. All of these measurements are taken at the field scale through a very sparse limited number of deep boreholes. These types of data are often reduced to probability distribution functions for predictive modeling and simulation in a stochastic framework such as a stochastic discrete fracture network. Stochastic discrete fracture network models enable, through Monte Carlo realizations and simulations, probabilistic assessment of flow and transport phenomena that are not adequately captured using continuum models. Despite the fundamental uncertainties inherited within the probabilistic reduction of the sparse data collected, very little work has been conducted on quantifying uncertainty on the reduced probabilistic distribution functions. In the current study, using nested Monte Carlo simulations, we present the impact of parameter uncertainties of the distribution functions of fracture density, orientation, aperture and size on the flow and transport using topological measures such as fracture connectivity, physical characteristics such as effective hydraulic conductivity tensors, and
The Effect of Loading Rate on Hydraulic Fracturing in Synthetic Granite - a Discrete Element Study
NASA Astrophysics Data System (ADS)
Tomac, I.; Gutierrez, M.
2015-12-01
Hydraulic fracture initiation and propagation from a borehole in hard synthetic rock is modeled using the two dimensional Discrete Element Method (DEM). DEM uses previously established procedure for modeling the strength and deformation parameters of quasi-brittle rocks with the Bonded Particle Model (Itasca, 2004). A series of simulations of laboratory tests on granite in DEM serve as a reference for synthetic rock behavior. Fracturing is enabled by breaking parallel bonds between DEM particles as a result of the local stress state. Subsequent bond breakage induces fracture propagation during a time-stepping procedure. Hydraulic fracturing occurs when pressurized fluid induces hoop stresses around the wellbore which cause rock fracturing and serves for geo-reservoir permeability enhancement in oil, gas and geothermal industries. In DEM, a network of fluid pipes and reservoirs is used for mathematical calculation of fluid flow through narrow channels between DEM particles, where the hydro-mechanical coupling is fully enabled. The fluid flow calculation is superimposed with DEM stress-strain calculation at each time step. As a result, the fluid pressures during borehole pressurization in hydraulic fracturing, as well as, during the fracture propagation from the borehole, can be simulated. The objective of this study is to investigate numerically a hypothesis that fluid pressurization rate, or the fluid flow rate, influences upon character, shape and velocity of fracture propagation in rock. The second objective is to better understand and define constraints which are important for successful fracture propagation in quasi-brittle rock from the perspective of flow rate, fluid density, viscosity and compressibility relative to the rock physical properties. Results from this study indicate that not only too high fluid flow rates cause fracture arrest and multiple fracture branching from the borehole, but also that the relative compressibility of fracturing fluid and
Zhou, Jing; Huang, Hai; Deo, Milind
2015-10-01
The interaction between hydraulic fractures (HF) and natural fractures (NF) will lead to complex fracture networks due to the branching and merging of natural and hydraulic fractures in unconventional reservoirs. In this paper, a newly developed hydraulic fracturing simulator based on discrete element method is used to predict the generation of complex fracture network in the presence of pre-existing natural fractures. By coupling geomechanics and reservoir flow within a dual lattice system, this simulator can effectively capture the poro-elastic effects and fluid leakoff into the formation. When HFs are intercepting single or multiple NFs, complex mechanisms such as direct crossing, arresting, dilating and branching can be simulated. Based on the model, the effects of injected fluid rate and viscosity, the orientation and permeability of NFs and stress anisotropy on the HF-NF interaction process are investigated. Combined impacts from multiple parameters are also examined in the paper. The numerical results show that large values of stress anisotropy, intercepting angle, injection rate and viscosity will impede the opening of NFs.
Identification of Fracture Toughness for Discrete Damage Mechanics Analysis of Glass-Epoxy Laminates
NASA Astrophysics Data System (ADS)
Barbero, E. J.; Cosso, F. A.; Martinez, X.
2014-08-01
A methodology for determination of the intralaminar fracture toughness is presented, based on fitting discrete damage mechanics (DDM) model predictions to available experimental data. DDM is constitutive model that, when incorporated into commercial finite element software via user material subroutines, is able to predict intralaminar transverse and shear damage initiation and evolution in terms of the fracture toughness of the composite. The applicability of the DDM model is studied by comparison to available experimental data for Glass-Epoxy laminates. Sensitivity of the DDM model to h- and p-refinement is studied. Also, the effect of in-situ correction of strength is highlighted.
NASA Astrophysics Data System (ADS)
Makedonska, N.; Karra, S.; Painter, S. L.; Viswanathan, H.; Gable, C. W.
2014-12-01
Natural gas from unconventional fossil energy sources such as shale and tight gasformations has a profound impact on US energy independence. The current state ofproduction of methane and other hydrocarbons from low permeability shale involvesprocesses such as hydraulic fracturing of rock, multiphase flow, and recovery of the gasvia these fractures. Although hydraulic fracturing has been used for the past couple ofdecades, little is known about the underlying mechanisms behind the production curvesthat are seen in the field, such as, reasons for 50-60% decline after the first year in typicalproduction curves.Numerical experiments on a realistic fractured shale system are presented toidentify the effect of complex flow of gas in fractures and matrix diffusion on theproduction curve. For characterizing flow, including the characteristics and geometriesfor the fracture networks, we use a methodology that incorporates a recently developeddiscrete fracture network meshing approach [1], which is combined with the highlyparallel PFLOTRAN subsurface flow and reactive transport code [2] and a new particletracking capability [3]. The results of this reservoir-scale methodology for analyzing thedecline in gas production rates indicate dominant flow in fractures in the initial highproduction rate. Increase in matrix diffusivity improves production recovery after theinitial production of gas from fractures. Moreover, it is observed that increasing aperturevariability within a single fracture has little effect on the production compared tovariations of the mean fracture aperture from fracture to fracture in a fracture network. [1] Hyman, J.D., Gable C.W., Painter S.L., and Makedonska N., ConformingDelaunay Triangulation of Stochastically Generated Three Dimensional DiscreteFracture Networks: a Feature Rejection Algorithm f or Meshing Strategy, SIAM J.Sci. Comput, 2014 (in press). [2] Lichtner, P.C., Hammond G.E., Lu C., Karra S., Bisht G., Mills R.T., and KumarJ., PFLOTRAN User
Shale Fracture Analysis using the Combined Finite-Discrete Element Method
NASA Astrophysics Data System (ADS)
Carey, J. W.; Lei, Z.; Rougier, E.; Knight, E. E.; Viswanathan, H.
2014-12-01
Hydraulic fracturing (hydrofrac) is a successful method used to extract oil and gas from highly carbonate rocks like shale. However, challenges exist for industry experts estimate that for a single $10 million dollar lateral wellbore fracking operation, only 10% of the hydrocarbons contained in the rock are extracted. To better understand how to improve hydrofrac recovery efficiencies and to lower its costs, LANL recently funded the Laboratory Directed Research and Development (LDRD) project: "Discovery Science of Hydraulic Fracturing: Innovative Working Fluids and Their Interactions with Rocks, Fractures, and Hydrocarbons". Under the support of this project, the LDRD modeling team is working with the experimental team to understand fracture initiation and propagation in shale rocks. LANL's hybrid hydro-mechanical (HM) tool, the Hybrid Optimization Software Suite (HOSS), is being used to simulate the complex fracture and fragment processes under a variety of different boundary conditions. HOSS is based on the combined finite-discrete element method (FDEM) and has been proven to be a superior computational tool for multi-fracturing problems. In this work, the comparison of HOSS simulation results to triaxial core flooding experiments will be presented.
Discrete Modeling of Early-Life Thermal Fracture in Ceramic Nuclear Fuel
Spencer, Benjamin W.; Huang, Hai; Dolbow, John E.; Hales, Jason D.
2015-03-01
Fracturing of ceramic fuel pellets heavily influences performance of light water reactor (LWR) fuel. Early in the life of fuel, starting with the initial power ramp, large thermal gradients cause high tensile hoop and axial stresses in the outer region of the fuel pellets, resulting in the formation of radial and axial cracks. Circumferential cracks form due to thermal gradients that occur when the power is ramped down. These thermal cracks cause the fuel to expand radially, closing the pellet/cladding gap and enhancing the thermal conductance across that gap, while decreasing the effective conductivity of the fuel in directions normal to the cracking. At lower length scales, formation of microcracks is an important contributor to the decrease in bulk thermal conductivity that occurs over the life of the fuel as the burnup increases. Because of the important effects that fracture has on fuel performance, a realistic, physically based fracture modeling capability is essential to predict fuel behavior in a wide variety of normal and abnormal conditions. Modeling fracture within the context of the finite element method, which is based on continuous interpolations of solution variables, has always been challenging because fracture is an inherently discontinuous phenomenon. Work is underway at Idaho National Laboratory to apply two modeling techniques model fracture as a discrete displacement discontinuity to nuclear fuel: The extended finite element method (XFEM), and discrete element method (DEM). XFEM is based on the standard finite element method, but with enhancements to represent discontinuous behavior. DEM represents a solid as a network of particles connected by bonds, which can arbitrarily fail if a fracture criterion is reached. This paper presents initial results applying the aforementioned techniques to model fuel fracturing. This work has initially focused on early life behavior of ceramic LWR fuel. A coupled thermal-mechanical XFEM method that includes
3D Discrete Element Model with 1 Million Particles: an Example of Hydro-fracturing
NASA Astrophysics Data System (ADS)
Liu, C.; Pollard, D. D.
2013-12-01
The Discrete Element Method (DEM) permits large relative motion and breakage of elements, and does not require re-meshing, for example as would the Finite Element Method. DEM has a wide range of applications in the fields of solid-earth geophysics, geomechanics, mining engineering, and structural geology. However, due to the computational cost, particle numbers of discrete element models are generally less than a few tens of thousands, which limits the applications. A new 3D DEM system 'MatDEM' can complete dynamic simulations of one million particles. The conversion formulas between particle parameters and model mechanical properties were derived, and the conversion of energy in DEM can be simulated. In a recent paper (Liu et al., 2013, JGR), the analytical solutions of elastic properties and failure modes of a 2D close-packed discrete element model were proposed. Based on these theoretical results, it is easy to create materials using DEM, which have similar mechanical properties to rock. Given the mechanical properties and state of stress, geologists and engineers can investigate the characteristics of rock deformation and failure under different conditions. MatDEM provides an alternative way to study the micro-macro relationships of rock and soil, and the evolution of geologic structures. As an example, MatDEM was used to investigate the generation and development of fluid driven fractures around a micro pore. The simulation result of fractures of an anisotropic 3D model, which includes 1 million particles, is demonstrated. Via parallel computing technology, MatDEM may handle tens of millions of particles in near future. Left: Fluid pressure is applied in the pore to generate fractures. Right: Simulation results (black segments represent fractures).
Discrete fracture hydromechanical model for the disturbed rock zone in a clay rock
NASA Astrophysics Data System (ADS)
Asahina, D.; Houseworth, J. E.; Birkholzer, J. T.
2013-12-01
We have developed a coupled thermal-hydrological-mechanical (THM) fracture damage model, TOUGH-RBSN, to investigate the behavior of fracture generation and evolution in rock in the presence of perturbations to THM conditions. This model combines the capabilities of the TOUGH2 simulator to represent thermal-hydrological processes with a rigid-body-spring-network (RBSN) model, a type of discrete modeling, to treat geomechanical and fracture-damage processes. In particular, the development and evolution of fractures in the excavation damaged zone (EDZ) of a clay rock, with application to high-level nuclear waste disposal, is a focus for this model development. Previously, the TOUGH-RBSN approach has been used to model fracture damage under tensile conditions as a result of desiccation shrinkage. The next phase of model testing will be application to the HG-A test being conducted at the Mont Terri underground research laboratory (URL) near Saint-Ursanne, Switzerland. This test is being conducted in a 13-m long, 1-m diameter microtunnel in the Opalinus clay rock in which a test section at the far end of the microtunnel is isolated using a packer. The test is specifically targeted to observe how fluids injected into the test section penetrate into the rock, with particular emphasis on the EDZ. The HG-A microtunnel was excavated in 2005 and subsequent mapping of the tunnel surface shows preferential fracturing and tunnel breakouts along zones where bedding planes are tangential to the tunnel wall and where faults intercept the tunnel. It appears that the EDZ fracture damage can be attributed to both tensile and shear fracturing mechanisms. A series of injection tests with water and gas have been performed which also show preferential invasion of the fluid pressure along the observed damage zones, as well as fracture self-sealing over time. The TOUGH-RBSN approach has been successfully applied to modeling fracture driven by predominately tensile loading, whereas only
Nitao, J.J.; Buscheck, T.A.
1995-11-09
An in situ heater test was performed at G-Tunnel, Nevada Nuclear Test Site, to investigate the thermal-hydrological response of unsaturated, fractured volcanic tuff under conditions similar to those at Yucca Mountain. The NUFT flow and transport code was used to model the test using discrete-fracture and equivalent-continuum approaches. Nonequilibrium fracture flow and thermal buoyant gas-phase convection were found to be the likely causes for observed lack of condensate imbibition into the matrix. The potential repository at Yucca Mountain was also modeled. Disequilibrium fracture flow is predicted to occur for less than a hundred years after emplacement followed by a period of fracture-matrix equilibrium, during which the equivalent-continuum and discrete-fracture models give almost identical results.
Particle tracking approach for transport in three-dimensional discrete fracture networks
Makedonska, Nataliia; Painter, Scott L; Bui, Quan M; Gable, Carl; Karra, Satish
2015-01-01
The discrete fracture network (DFN) model is a method to mimic discrete pathways for fluid flow through a fractured low-permeable rock mass, and may be combined with particle tracking simulations to address solute transport. However, experience has shown that it is challenging to obtain accurate transport results in three-dimensional DFNs because of the high computational burden and difficulty in constructing a high-quality unstructured computational mesh on simulated fractures. We present a new particle tracking capability, which is adapted to control volume (Voronoi polygons) flow solutions on unstructured grids (Delaunay triangulations) on three-dimensional DFNs. The locally mass-conserving finite-volume approach eliminates massmore » balance related problems during particle tracking. The scalar fluxes calculated for each control volume face by the flow solver are used to reconstruct a Darcy velocity at each control volume centroid. The groundwater velocities can then be continuously interpolated to any point in the domain of interest. The control volumes at fracture intersections are split into four pieces and the velocity is reconstructed independently on each piece, which results in multiple groundwater velocities at the intersection, one for each fracture on each side of the intersection line. This technique enables detailed particle transport representation through a complex DFN structure. Verified for small DFNs, the new simulation capability enables numerical experiments on advective transport in large DFNs to be performed. We demonstrate this particle transport approach on a DFN model using parameters similar to those of crystalline rock at a proposed geologic repository for spent nuclear fuel in Forsmark, Sweden.« less
Particle tracking approach for transport in three-dimensional discrete fracture networks
Makedonska, Nataliia; Painter, Scott L; Bui, Quan M; Gable, Carl; Karra, Satish
2015-01-01
The discrete fracture network (DFN) model is a method to mimic discrete pathways for fluid flow through a fractured low-permeable rock mass, and may be combined with particle tracking simulations to address solute transport. However, experience has shown that it is challenging to obtain accurate transport results in three-dimensional DFNs because of the high computational burden and difficulty in constructing a high-quality unstructured computational mesh on simulated fractures. We present a new particle tracking capability, which is adapted to control volume (Voronoi polygons) flow solutions on unstructured grids (Delaunay triangulations) on three-dimensional DFNs. The locally mass-conserving finite-volume approach eliminates mass balance related problems during particle tracking. The scalar fluxes calculated for each control volume face by the flow solver are used to reconstruct a Darcy velocity at each control volume centroid. The groundwater velocities can then be continuously interpolated to any point in the domain of interest. The control volumes at fracture intersections are split into four pieces and the velocity is reconstructed independently on each piece, which results in multiple groundwater velocities at the intersection, one for each fracture on each side of the intersection line. This technique enables detailed particle transport representation through a complex DFN structure. Verified for small DFNs, the new simulation capability enables numerical experiments on advective transport in large DFNs to be performed. We demonstrate this particle transport approach on a DFN model using parameters similar to those of crystalline rock at a proposed geologic repository for spent nuclear fuel in Forsmark, Sweden.
NASA Astrophysics Data System (ADS)
Ţene, Matei; Al Kobaisi, Mohammed Saad; Hajibeygi, Hadi
2016-09-01
This paper introduces an Algebraic MultiScale method for simulation of flow in heterogeneous porous media with embedded discrete Fractures (F-AMS). First, multiscale coarse grids are independently constructed for both porous matrix and fracture networks. Then, a map between coarse- and fine-scale is obtained by algebraically computing basis functions with local support. In order to extend the localization assumption to the fractured media, four types of basis functions are investigated: (1) Decoupled-AMS, in which the two media are completely decoupled, (2) Frac-AMS and (3) Rock-AMS, which take into account only one-way transmissibilities, and (4) Coupled-AMS, in which the matrix and fracture interpolators are fully coupled. In order to ensure scalability, the F-AMS framework permits full flexibility in terms of the resolution of the fracture coarse grids. Numerical results are presented for two- and three-dimensional heterogeneous test cases. During these experiments, the performance of F-AMS, paired with ILU(0) as second-stage smoother in a convergent iterative procedure, is studied by monitoring CPU times and convergence rates. Finally, in order to investigate the scalability of the method, an extensive benchmark study is conducted, where a commercial algebraic multigrid solver is used as reference. The results show that, given an appropriate coarsening strategy, F-AMS is insensitive to severe fracture and matrix conductivity contrasts, as well as the length of the fracture networks. Its unique feature is that a fine-scale mass conservative flux field can be reconstructed after any iteration, providing efficient approximate solutions in time-dependent simulations.
Discrete fracture network modelling of Groß Schönebeck stimulation treatment
NASA Astrophysics Data System (ADS)
Urpi, L.; Zimmermann, G.; Blöcher, G.; van Wees, J. D. A. M.; Wassing, B.
2012-04-01
Microseismic events associated to geothermal reservoir are recorded before, during and after the establishment of an Enhanced Geothermal System (EGS). Differentiating recorded seismicity between natural and induced can be ambiguous, but reservoir response to stimulations treatment can be modelled and give useful insights in designing treatment. Our model reproduces the stimulation treatment done at Groß Schönebeck reservoir. The stimulation target was the volcanic layer of the Rotliegend formation at 4km depth. The treatment increased by a factor 22 the productivity index (volume of fluid produced per unit of time per drawdown) of the reservoir, while recorded seismicity was lower than expected (<100 events recorded, max magnitude mL=-1). The microevents signal properties are compatible with shearing events and their locations along a plane give us little information on the processes occurring during the water injection. Bottom hole injection pressure was over the minimum stress, therefore we expect fracture opening due to both tensile and shearing opening. The interaction between a primary fracture initiating from the wellbore due to the injection and the secondary discrete fracture network (DFN) has been modelled with a hydraulic fracturing simulator. The DFN has been populated on the basis of available borehole data and lithological properties determined from rock sample analysis. Heat transfer has been computed since the temperature difference between injected fluid and reservoir rocks is over 100°C, inducing thermoelastic stresses around the the fractures. The model shows the importance of the natural fracture network, perturbing the reservoir state at distance in direction parallel to the minimum stress. Results are compared with commercial tensile and shear fracture models, which are compatible in term of predicting increased productivity of the well. We justify the absence of larger magnitude events after shut-in with the low natural seismicity and the
Fracture and impulse based finite-discrete element modeling of fragmentation
NASA Astrophysics Data System (ADS)
Paluszny, A.; Tang, X. H.; Zimmerman, R. W.
2013-11-01
A numerical method for fragmentation is presented that combines the finite element method with the impulse-based discrete element method (impulse-based FDEM). In contrast to existing methods, fragments are not represented as a conglomeration of spheres; instead, their shapes are represented using solid modeling techniques, and are the result of multiple fracture growth. Fracture growth within each three-dimensional fragment is controlled by stress intensity factors computed using the finite element method and the reduced virtual integration technique. Non-convex fragment interaction and movement is modeled using impulse dynamics, rather than a penalty-based method. Collisions leading to fracture are handled individually by propagating pre-existing internal flaws and cracks. The method utilizes decoupled geometry and mesh representation, and local failure and propagation criteria. Fractures that reach volume boundaries lead to further fragmentation. The approach is demonstrated by the fragmentation of a sphere, which exhibits a velocity-dependent fragment size distribution. The distribution is characterized by a two-parameter Weibull distribution, an emergent property of the simulation. Results are in good agreement with experimental data.
Discrete Fracture Network Models for Risk Assessment of Carbon Sequestration in Coal
Jack Pashin; Guohai Jin; Chunmiao Zheng; Song Chen; Marcella McIntyre
2008-07-01
A software package called DFNModeler has been developed to assess the potential risks associated with carbon sequestration in coal. Natural fractures provide the principal conduits for fluid flow in coal-bearing strata, and these fractures present the most tangible risks for the leakage of injected carbon dioxide. The objectives of this study were to develop discrete fracture network (DFN) modeling tools for risk assessment and to use these tools to assess risks in the Black Warrior Basin of Alabama, where coal-bearing strata have high potential for carbon sequestration and enhanced coalbed methane recovery. DFNModeler provides a user-friendly interface for the construction, visualization, and analysis of DFN models. DFNModeler employs an OpenGL graphics engine that enables real-time manipulation of DFN models. Analytical capabilities in DFNModeler include display of structural and hydrologic parameters, compartmentalization analysis, and fluid pathways analysis. DFN models can be exported to third-party software packages for flow modeling. DFN models were constructed to simulate fracturing in coal-bearing strata of the upper Pottsville Formation in the Black Warrior Basin. Outcrops and wireline cores were used to characterize fracture systems, which include joint systems, cleat systems, and fault-related shear fractures. DFN models were constructed to simulate jointing, cleating, faulting, and hydraulic fracturing. Analysis of DFN models indicates that strata-bound jointing compartmentalizes the Pottsville hydrologic system and helps protect shallow aquifers from injection operations at reservoir depth. Analysis of fault zones, however, suggests that faulting can facilitate cross-formational flow. For this reason, faults should be avoided when siting injection wells. DFN-based flow models constructed in TOUGH2 indicate that fracture aperture and connectivity are critical variables affecting the leakage of injected CO{sub 2} from coal. Highly transmissive joints
Senapati, Rajeev; Zhang Jianmei
2010-02-22
Advanced ceramic materials have been extensively applied in aerospace, automobile and other industries. However, the reliability of the advanced ceramics is a major concern because of the brittle nature of the materials. In this paper, combination of nondestructive testing and numerical modeling Discrete Element Method is proposed to identify the fracture origin in ceramics. The nondestructive testing--laser scattering technology is first performed on the ceramic components to reveal the machining-induced damage such as cracks and the material-inherent flaws such as voids, then followed by the four point bending test. Discrete Element software package PFC{sup 2D} is used to simulate the four point bending test and try to identify where the fractures start. The numerical representation of the ceramic materials is done by generating a densely packed particle system using the specimen genesis procedure and then applying the suitable microparameters to the particle system. Simulation of four point bending test is performed on materials having no defects, materials having manufacturing-induced defects like cracks, and materials having material-inherent flaws like voids. The initiation and propagation of defects is modeled and the mean contact force on the loading ball is also plotted. The simulation prediction results are well in accordance with the nondestructive testing results.
NASA Astrophysics Data System (ADS)
Zhang, Na; Yao, Jun; Huang, Zhaoqin; Wang, Yueying
2013-06-01
Numerical simulation in naturally fractured media is challenging because of the coexistence of porous media and fractures on multiple scales that need to be coupled. We present a new approach to reservoir simulation that gives accurate resolution of both large-scale and fine-scale flow patterns. Multiscale methods are suitable for this type of modeling, because it enables capturing the large scale behavior of the solution without solving all the small features. Dual-porosity models in view of their strength and simplicity can be mainly used for sugar-cube representation of fractured media. In such a representation, the transfer function between the fracture and the matrix block can be readily calculated for water-wet media. For a mixed-wet system, the evaluation of the transfer function becomes complicated due to the effect of gravity. In this work, we use a multiscale finite element method (MsFEM) for two-phase flow in fractured media using the discrete-fracture model. By combining MsFEM with the discrete-fracture model, we aim towards a numerical scheme that facilitates fractured reservoir simulation without upscaling. MsFEM uses a standard Darcy model to approximate the pressure and saturation on a coarse grid, whereas fine scale effects are captured through basis functions constructed by solving local flow problems using the discrete-fracture model. The accuracy and the robustness of MsFEM are shown through several examples. In the first example, we consider several small fractures in a matrix and then compare the results solved by the finite element method. Then, we use the MsFEM in more complex models. The results indicate that the MsFEM is a promising path toward direct simulation of highly resolution geomodels.
DFNWorks. A discrete fracture network framework for modeling subsurface flow and transport
Hyman, Jeffrey D.; Karra, Satish; Makedonska, Nataliia; Gable, Carl W.; Painter, Scott L.; Viswanathan, Hari S.
2015-08-10
DFNWorks is a parallalized computational suite to generate three-dimensional discrete fracture networks (DFN) and simulate flow and transport. Developed at Los Alamos National Laboratory over the past five years, it has been used to study flow and transport in fractured media at scales ranging from millimeters to kilometers. The networks are created and meshed using dfnGen, which combines fram (the feature rejection algorithm for meshing) methodology to stochastically generate three-dimensional DFNs on the basis of site specific data with the LaGriT meshing toolbox to create a high-quality computational mesh representation, specifically a conforming Delaunay triangulation suitable for high performance computingmore » finite volume solvers, of the DFN in an intrinsically parallel fashion. Flow through the network is simulated in dfnFlow, which utilizes the massively parallel subsurface flow and reactive transport finite volume code pflotran. A Lagrangian approach to simulating transport through the DFN is adopted within dfnTrans, which is an extension of the walkabout particle tracking method to determine pathlines through the DFN. Example applications of this suite in the areas of nuclear waste repository science, hydraulic fracturing and CO2 sequestration are also included.« less
dfnWorks: A discrete fracture network framework for modeling subsurface flow and transport
Hyman, Jeffrey D.; Karra, Satish; Makedonska, Nataliia; Gable, Carl W.; Painter, Scott L.; Viswanathan, Hari S.
2015-11-01
DFNWORKS is a parallelized computational suite to generate three-dimensional discrete fracture networks (DFN) and simulate flow and transport. Developed at Los Alamos National Laboratory over the past five years, it has been used to study flow and transport in fractured media at scales ranging from millimeters to kilometers. The networks are created and meshed using DFNGEN, which combines FRAM (the feature rejection algorithm for meshing) methodology to stochastically generate three-dimensional DFNs with the LaGriT meshing toolbox to create a high-quality computational mesh representation. The representation produces a conforming Delaunay triangulation suitable for high performance computing finite volume solvers in an intrinsically parallel fashion. Flow through the network is simulated in dfnFlow, which utilizes the massively parallel subsurface flow and reactive transport finite volume code PFLOTRAN. A Lagrangian approach to simulating transport through the DFN is adopted within DFNTRANS to determine pathlines and solute transport through the DFN. Example applications of this suite in the areas of nuclear waste repository science, hydraulic fracturing and CO_{2} sequestration are also included.
DFNWORKS: A discrete fracture network framework for modeling subsurface flow and transport
NASA Astrophysics Data System (ADS)
Hyman, Jeffrey D.; Karra, Satish; Makedonska, Nataliia; Gable, Carl W.; Painter, Scott L.; Viswanathan, Hari S.
2015-11-01
DFNWORKS is a parallelized computational suite to generate three-dimensional discrete fracture networks (DFN) and simulate flow and transport. Developed at Los Alamos National Laboratory over the past five years, it has been used to study flow and transport in fractured media at scales ranging from millimeters to kilometers. The networks are created and meshed using DFNGEN, which combines FRAM (the feature rejection algorithm for meshing) methodology to stochastically generate three-dimensional DFNs with the LAGRIT meshing toolbox to create a high-quality computational mesh representation. The representation produces a conforming Delaunay triangulation suitable for high performance computing finite volume solvers in an intrinsically parallel fashion. Flow through the network is simulated in DFNFLOW, which utilizes the massively parallel subsurface flow and reactive transport finite volume code PFLOTRAN. A Lagrangian approach to simulating transport through the DFN is adopted within DFNTRANS to determine pathlines and solute transport through the DFN. Example applications of this suite in the areas of nuclear waste repository science, hydraulic fracturing and CO2 sequestration are also included.
DFNWorks. A discrete fracture network framework for modeling subsurface flow and transport
Hyman, Jeffrey D.; Karra, Satish; Makedonska, Nataliia; Gable, Carl W.; Painter, Scott L.; Viswanathan, Hari S.
2015-08-10
DFNWorks is a parallalized computational suite to generate three-dimensional discrete fracture networks (DFN) and simulate flow and transport. Developed at Los Alamos National Laboratory over the past five years, it has been used to study flow and transport in fractured media at scales ranging from millimeters to kilometers. The networks are created and meshed using dfnGen, which combines fram (the feature rejection algorithm for meshing) methodology to stochastically generate three-dimensional DFNs on the basis of site specific data with the LaGriT meshing toolbox to create a high-quality computational mesh representation, specifically a conforming Delaunay triangulation suitable for high performance computing finite volume solvers, of the DFN in an intrinsically parallel fashion. Flow through the network is simulated in dfnFlow, which utilizes the massively parallel subsurface flow and reactive transport finite volume code pflotran. A Lagrangian approach to simulating transport through the DFN is adopted within dfnTrans, which is an extension of the walkabout particle tracking method to determine pathlines through the DFN. Example applications of this suite in the areas of nuclear waste repository science, hydraulic fracturing and CO2 sequestration are also included.
dfnWorks: A discrete fracture network framework for modeling subsurface flow and transport
Hyman, Jeffrey D.; Karra, Satish; Makedonska, Nataliia; Gable, Carl W.; Painter, Scott L.; Viswanathan, Hari S.
2015-11-01
DFNWORKS is a parallelized computational suite to generate three-dimensional discrete fracture networks (DFN) and simulate flow and transport. Developed at Los Alamos National Laboratory over the past five years, it has been used to study flow and transport in fractured media at scales ranging from millimeters to kilometers. The networks are created and meshed using DFNGEN, which combines FRAM (the feature rejection algorithm for meshing) methodology to stochastically generate three-dimensional DFNs with the LaGriT meshing toolbox to create a high-quality computational mesh representation. The representation produces a conforming Delaunay triangulation suitable for high performance computing finite volume solvers in anmore » intrinsically parallel fashion. Flow through the network is simulated in dfnFlow, which utilizes the massively parallel subsurface flow and reactive transport finite volume code PFLOTRAN. A Lagrangian approach to simulating transport through the DFN is adopted within DFNTRANS to determine pathlines and solute transport through the DFN. Example applications of this suite in the areas of nuclear waste repository science, hydraulic fracturing and CO2 sequestration are also included.« less
Birkholzer, J.; Karasaki, K.
1996-07-01
Fracture network simulators have extensively been used in the past for obtaining a better understanding of flow and transport processes in fractured rock. However, most of these models do not account for fluid or solute exchange between the fractures and the porous matrix, although diffusion into the matrix pores can have a major impact on the spreading of contaminants. In the present paper a new finite element code TRIPOLY is introduced which combines a powerful fracture network simulator with an efficient method to account for the diffusive interaction between the fractures and the adjacent matrix blocks. The fracture network simulator used in TRIPOLY features a mixed Lagrangian-Eulerian solution scheme for the transport in fractures, combined with an adaptive gridding technique to account for sharp concentration fronts. The fracture-matrix interaction is calculated with an efficient method which has been successfully used in the past for dual-porosity models. Discrete fractures and matrix blocks are treated as two different systems, and the interaction is modeled by introducing sink/source terms in both systems. It is assumed that diffusive transport in the matrix can be approximated as a one-dimensional process, perpendicular to the adjacent fracture surfaces. A direct solution scheme is employed to solve the coupled fracture and matrix equations. The newly developed combination of the fracture network simulator and the fracture-matrix interaction module allows for detailed studies of spreading processes in fractured porous rock. The authors present a sample application which demonstrate the codes ability of handling large-scale fracture-matrix systems comprising individual fractures and matrix blocks of arbitrary size and shape.
NASA Astrophysics Data System (ADS)
Zeeb, Conny; Frühwirt, Thomas; Konietzky, Heinz
2015-04-01
Key to a successful exploitation of deep geothermal reservoirs in a petrothermal environment is the hydraulic stimulation of the host rock to increase permeability. The presented research investigates the fracture propagation and interaction during hydraulic stimulation of multiple fractures in a highly anisotropic stress field. The presented work was conducted within the framework of the OPTIRISS project, which is a cooperation of industry partners and universities in Thuringia and Saxony (Federal States of Germany) and was funded by the European Fond for Regional Development. One objective was the design optimization of the subsurface geothermal heat exchanger (SGHE) by means of numerical simulations. The presented simulations were conducted applying 3DEC (Itasca™), a software tool based on the discrete element method. The simulation results indicate that the main direction of fracture propagation is towards lower stresses and thus towards the biosphere. Therefore, barriers might be necessary to limit fracture propagation to the designated geological formation. Moreover, the hydraulic stimulation significantly alters the stresses in the vicinity of newly created fractures. Especially the change of the minimum stress component affects the hydraulic stimulation of subsequent fractures, which are deflected away from the previously stimulated fractures. This fracture deflection can render it impossible to connect all fractures with a second borehole for the later production. The results of continuative simulations indicate that a fracture deflection cannot be avoided completely. Therefore, the stage alignment was modified to minimize fracture deflection by varying (1) the pauses between stages, (2) the spacing's between adjacent stages, and (3) the angle between stimulation borehole and minimum stress component. An optimum SGHE design, which implies that all stimulated fractures are connected to the production borehole, can be achieved by aligning the stimulation
Anna, L.O.
1998-09-01
Discrete-fracture modeling is part of site characterization for evaluating Yucca Mountain, Nye County, Nevada, as a potential high-level radioactive-waste repository site. Because most of the water and gas flow may be in fractures in low-porosity units, conventional equivalent-continuum models do not adequately represent the flow system. Discrete-fracture modeling offers an alternative to the equivalent-continuum method. This report describes how discrete-fracture networks can be constructed and used to answer concerns about the flow system at Yucca Mountain, including quantifying fracture connectivity, deriving directional-permeability distributions for one-and two-phase flow, determining parameters of anisotropy at different scales, and determining at what scale the rock functions as an equivalent continuum. A three-dimensional discrete-fracture model was developed to investigate the effects of fractures on flow of water and gas in the Topopah Spring tuff of Miocene age in the Exploratory Studies Facility at Yucca Mountain. Fracture data, used as model input, were taken exclusively from detailed line surveys in the Exploratory Studies Facility and converted into input parameters for simulation. A simulated fracture network was calibrated to field data. The simulated discrete fracture network was modified by eliminating nonconductive fractures determined from field-derived permeabilities. Small fractures also were removed from the simulated network without affecting the overall connectivity. Fractures, as much as 1.50 meters in length, were eliminated (a large percentage of the total number of fractures) from the network without altering the number of connected pathways. The analysis indicates that the fracture system in the Exploratory Studies Facility has numerous connected fractures that have relatively large permeabilities, but there are relatively few connected pathways across the simulated region. The fracture network was, therefore, sparse.
NASA Astrophysics Data System (ADS)
Keralavarma, Shyam Mohan
behavior of a large number of nano and micro scale defects such as dislocations, vacancies and grain boundaries. Continuum models relate macroscopically observable quantities such as stress and strain by coarse graining the discrete defect microstructure. While continuum models provide a good approximation for the effective behavior of bulk materials, several deviations have been observed in experiments at small scales such as an intrinsic size dependence of the material strength. Discrete dislocation dynamics (DD) is a mesoscale method for obtaining the mechanical response of a material by direct simulation of the motion and interactions of dislocations. The model incorporates an intrinsic length scale in the dislocation Burgers vector and potentially allows for size dependent mechanical behavior to emerge naturally from the dynamics of the dislocation ensemble. In the second part of this dissertation, a simplified twodimensional DD model is employed to study several phenomena of practical interest such as strain hardening under homogeneous deformation, growth of microvoids in a crystalline matrix and creep of single crystals at elevated temperatures. These studies have been enabled by several recent enhancements to the existing two-dimensional DD framework described in Chapter V. The main contributions from this research are: (i) development of a fully anisotropic continuum model of void growth for use in ductile fracture simulations and (ii) enhancing the capabilities of an existing two-dimensional DD framework for large scale simulations in complex domains and at elevated temperatures.
NASA Astrophysics Data System (ADS)
Griffith, W. A.; Ghaffari, H.; Barber, T. J.; Borjas, C.
2015-12-01
The motions of Earth's tectonic plates are typically measured in millimeters to tens of centimeters per year, seemingly confirming the generally-held view that tectonic processes are slow, and have been throughout Earth's history. In line with this perspective, the vast majority of laboratory rock mechanics research focused on failure in the brittle regime has been limited to experiments utilizing slow loading rates. On the other hand, many natural processes that pose significant risk for humans (e.g., earthquakes and extraterrestrial impacts), as well as risks associated with human activities (blow-outs, explosions, mining and mine failures, projectile penetration), occur at rates that are hundreds to thousands of times faster than those typically simulated in the laboratory. Little experimental data exists to confirm or calibrate theoretical models explaining the connection between these dramatic events and the pulverized rocks found in fault zones, impacts, or explosions; however the experimental data that does exist is thought-provoking: At the earth's surface, the process of brittle fracture passes through a critical transition in rocks at high strain rates (101-103s-1) between regimes of discrete fracture and distributed fragmentation, accompanied by a dramatic increase in strength. Previous experimental works on this topic have focused on key thresholds (e.g., peak stress, peak strain, average strain rate) that define this transition, but more recent work suggests that this transition is more fundamentally dependent on characteristics (e.g., shape) of the loading pulse and related microcrack dynamics, perhaps explaining why for different lithologies different thresholds more effectively define the pulverization transition. In this presentation we summarize some of our work focused on this transition, including the evolution of individual defects at the microscopic, microsecond scale and the energy budget associated with the brittle fragmentation process as a
NASA Astrophysics Data System (ADS)
Steefel, C. I.; Lichtner, P. C.
1998-08-01
A numerical multicomponent reactive transport model described fully in Steefel and Lichtner (1998)[Steefel, C.I., Lichtner, P.C., 1998. Multicomponent reactive transport in discrete fractures, I. Controls on reaction front geometry. J. Hydrol. (in press)] is used to simulate the infiltration of hyperalkaline groundwater along discrete fractures at Maqarin, Jordan, a site considered as a natural analogue to cement-bearing nuclear waste repositories. In the Eastern Springs area at Maqarin, two prominent sets of sub-parallel fractures trending NW-SE are approximately perpendicular to the local water table contours, with the slope of the water table indicating north-westward flow. Extensive mineralogic investigations [Alexander W.R. (Ed.), 1992. A natural analogue study of cement-buffered, hyperalkaline groundwaters and their interaction with a sedimentary host rock. NAgrA Technical Report (NTB 91-10), Wettingen, Switzerland; Milodowski, A.E., Hyslop, E.K., Pearce, J.M., Wetton, P.D., Kemp, S.J., Longworth, G., Hodginson, E., and Hughes, C.R., 1998. Mineralogy and geochemistry of the western springs area. In: Smellie, J.A.T. (ed.), 1998: Maqarin Natural Analogue Study: Phase III. SKB Technical Report TR98-04, Stockholm, Sweden] indicate that the width of intense rock alteration zone bordering the fractures changes from about 4 mm at one locality (the M1 sampling site) to approximately 1 mm 100 m to the north-west in the flow direction (the M2 site), suggesting a lessening of alteration intensity in that direction. Using this information, the dimensionless parameter δ v/φ D' (φ=porosity, D'=effective diffusion coefficient in rock matrix, δ=fracture aperture, and v=fluid velocity in the fracture) and measurements of the local hydraulic head gradient and effective diffusion coefficient in the rock matrix, a mean fracture aperture of 0.194 mm is calculated assuming the cubic law applies. This information, in combination with measured groundwater compositions at the
NASA Astrophysics Data System (ADS)
Makedonska, N.; Painter, S. L.; Hyman, J.; Karra, S.; Gable, C. W.; Viswanathan, H. S.
2015-12-01
Aperture variability within individual fractures is usually neglected in modeling flow and transport through fractured media. Typically, individual fractures are assumed to be homogeneous. However, in reality, individual fractures are heterogeneous, which may affect flow and transport in fractured media. The relative importance of including in-fracture variability in flow and transport modeling has been under debate for a long time. Previous studies have shown flow channeling on an individual fracture with internal variability, where the fracture is considered isolated from the rest of the fracture network. Although these studies yield some clear insights into the process, the boundary conditions are impractical for field-scale networks, where the realistic boundary conditions are determined by fracture connections in the network. Therefore, flow in a single fracture is controlled not only by in-fracture variability but also by boundary conditions. In order to address the question of the importance of in-fracture variability, the internal heterogeneity of every individual fracture is incorporated into a three-dimensional fracture network, represented by a composition of intersecting fractures. The new DFN simulation capability, dfnWorks, is used to generate a kilometer scale DFNs similar to the Forsmark, Sweden site. In our DFN model, the in-fracture aperture variability is scattered over each cell of the computational mesh along the fracture, representing by a stationary Gaussian random field with various correlation lengths. The Lagrangian particle tracking is conducted in multiple DFN realizations and the flow-dependent Lagrangian parameters, non-reacting travel time, τ, and cumulative reactivity parameter, β, are obtained along particles streamlines. It is shown that early particle travel times are more sensitive to in-fracture aperture variability than tails of travel time distributions, where no significant effect of the aperture variations and spatial
NASA Astrophysics Data System (ADS)
Leonardi, Christopher R.; McCullough, Jon W. S.; Jones, Bruce D.; Williams, John R.
2016-04-01
This paper describes the development of a computational framework that can be used to describe the electromagnetic excitation of rigid, spherical particles in suspension. In this model the mechanical interaction and kinematic behaviour of the particles is modelled using the discrete element method, while the surrounding fluid mechanics is modelled using the lattice Boltzmann method. Electromagnetic effects are applied to the particles as an additional set of discrete element forces, and the implementation of these effects was validated by comparison to the theoretical equations of point charges for Coulomb's law and the Lorentz force equation. Oscillating single and multiple particle tests are used to investigate the sensitivity of particle excitation to variations in particle charge, field strength, and frequency. The further capabilities of the model are then demonstrated by a numerical illustration, in which a hydraulic fracture fluid is excited and monitored within a hydraulic fracture. This modelling explores the feasibility of using particle vibrations within the fracture fluid to aid in the monitoring of fracture propagation in unconventional gas reservoirs.
NASA Astrophysics Data System (ADS)
Profit, Matthew; Dutko, Martin; Yu, Jianguo; Cole, Sarah; Angus, Doug; Baird, Alan
2016-04-01
This paper presents a novel approach to predict the propagation of hydraulic fractures in tight shale reservoirs. Many hydraulic fracture modelling schemes assume that the fracture direction is pre-seeded in the problem domain discretisation. This is a severe limitation as the reservoir often contains large numbers of pre-existing fractures that strongly influence the direction of the propagating fracture. To circumvent these shortcomings, a new fracture modelling treatment is proposed where the introduction of discrete fracture surfaces is based on new and dynamically updated geometrical entities rather than the topology of the underlying spatial discretisation. Hydraulic fracturing is an inherently coupled engineering problem with interactions between fluid flow and fracturing when the stress state of the reservoir rock attains a failure criterion. This work follows a staggered hydro-mechanical coupled finite/discrete element approach to capture the key interplay between fluid pressure and fracture growth. In field practice, the fracture growth is hidden from the design engineer and microseismicity is often used to infer hydraulic fracture lengths and directions. Microseismic output can also be computed from changes of the effective stress in the geomechanical model and compared against field microseismicity. A number of hydraulic fracture numerical examples are presented to illustrate the new technology.
Makedonska, Nataliia; Painter, Scott L.; Bui, Quan M.; Gable, Carl W.; Karra, Satish
2015-09-16
The discrete fracture network (DFN) model is a method to mimic discrete pathways for fluid flow through a fractured low-permeable rock mass, and may be combined with particle tracking simulations to address solute transport. However, experience has shown that it is challenging to obtain accurate transport results in three-dimensional DFNs because of the high computational burden and difficulty in constructing a high-quality unstructured computational mesh on simulated fractures. We present a new particle tracking capability, which is adapted to control volume (Voronoi polygons) flow solutions on unstructured grids (Delaunay triangulations) on three-dimensional DFNs. The locally mass-conserving finite-volume approach eliminates massmore » balance-related problems during particle tracking. The scalar fluxes calculated for each control volume face by the flow solver are used to reconstruct a Darcy velocity at each control volume centroid. The groundwater velocities can then be continuously interpolated to any point in the domain of interest. The control volumes at fracture intersections are split into four pieces, and the velocity is reconstructed independently on each piece, which results in multiple groundwater velocities at the intersection, one for each fracture on each side of the intersection line. This technique enables detailed particle transport representation through a complex DFN structure. Verified for small DFNs, the new simulation capability enables numerical experiments on advective transport in large DFNs to be performed. As a result, we demonstrate this particle transport approach on a DFN model using parameters similar to those of crystalline rock at a proposed geologic repository for spent nuclear fuel in Forsmark, Sweden.« less
Makedonska, Nataliia; Painter, Scott L.; Bui, Quan M.; Gable, Carl W.; Karra, Satish
2015-09-16
The discrete fracture network (DFN) model is a method to mimic discrete pathways for fluid flow through a fractured low-permeable rock mass, and may be combined with particle tracking simulations to address solute transport. However, experience has shown that it is challenging to obtain accurate transport results in three-dimensional DFNs because of the high computational burden and difficulty in constructing a high-quality unstructured computational mesh on simulated fractures. We present a new particle tracking capability, which is adapted to control volume (Voronoi polygons) flow solutions on unstructured grids (Delaunay triangulations) on three-dimensional DFNs. The locally mass-conserving finite-volume approach eliminates mass balance-related problems during particle tracking. The scalar fluxes calculated for each control volume face by the flow solver are used to reconstruct a Darcy velocity at each control volume centroid. The groundwater velocities can then be continuously interpolated to any point in the domain of interest. The control volumes at fracture intersections are split into four pieces, and the velocity is reconstructed independently on each piece, which results in multiple groundwater velocities at the intersection, one for each fracture on each side of the intersection line. This technique enables detailed particle transport representation through a complex DFN structure. Verified for small DFNs, the new simulation capability enables numerical experiments on advective transport in large DFNs to be performed. As a result, we demonstrate this particle transport approach on a DFN model using parameters similar to those of crystalline rock at a proposed geologic repository for spent nuclear fuel in Forsmark, Sweden.
NASA Astrophysics Data System (ADS)
Seok, E.; Gale, J.
2012-12-01
Fractured crystalline rocks, that have low matrix permeability, are characterized by strong spatial variability in permeability and porosity. However, most 3-D flow and transport models are constructed by developing a model mesh at the surface and then projecting this mesh as a series of layers to the assigned bottom of the model. The input parameters, especially permeability, are often assigned to individual layers as average properties within specific zones. These average properties vary from layer to layer based on the mapped hydrogeological variations but rarely embed the underlying spatial variability of the key parameters of permeability and porosity. To compare the impact on simulated flow in a discrete fracture of assigning average versus spatially varying permeability values, a one cubic metre of granitic rock with a single natural fracture was instrumented, loaded in a structural loading frame and a series of flow tests were completed at different normal loads. Water injection tests were completed in seventeen boreholes, both individually and in groups, that were drilled into this fracture plane providing a series of local permeability measurements and hydraulic heads distributed over the fracture plane as a function of normal stress. The fracture apertures were measured where the fracture intersected the four side of the granitic rock sample. The 3-D finite element model that gave a good match between the measured and computed flowrates used input parameters that were generated using the individual borehole measurements as constraints, such that the small scale variability, provided by the apertures measured around the outside of the block, was superimposed on the large scale variability. This laboratory and numerical model study demonstrated that flow models will provide a more representative measure of flow, and hence transport, if the input parameters properly reflect the small and large scale spatial variability in the key flow properties.
NASA Astrophysics Data System (ADS)
Cil, Mehmet B.; Alshibli, Khalid A.
2015-02-01
The constitutive behavior and deformation characteristics of uncemented granular materials are to a large extent derived from the fabric or geometry of the particle structure and the interparticle friction resulting from normal forces acting on particles or groups of particles. Granular materials consist of discrete particles with a fabric (microstructure) that changes under loading. Synchrotron micro-computed tomography (SMT) has emerged as a powerful non-destructive 3D scanning technique to study geomaterials. In this paper, SMT was used to acquire in situ scans of the oedometry test of a column of three silica sand particles. The sand is known as ASTM 20-30 Ottawa sand, and has a grain size between US sieves #20 (0.841 mm) and #30 (0.595 mm). The characteristics and evolution of particle fracture in sand were examined using SMT images, and a 3D discrete element method (DEM) was used to model the fracture behavior of sand particles. It adopts the bonded particle model to generate a crushable agglomerate that consists of a large number of small spherical sub-particles. The agglomerate shape matches the 3D physical shape of the tested sand particles by mapping the particle morphology from the SMT images. The paper investigates and discusses the influence of agglomerate packing (i.e., the number and size distribution of spherical sub-particles that constitute the agglomerate) and agglomerate shape on the fracture behavior of crushable particles.
Fu, P; Johnson, S M; Hao, Y; Carrigan, C R
2011-01-18
The primary objective of our current research is to develop a computational test bed for evaluating borehole techniques to enhance fluid flow and heat transfer in enhanced geothermal systems (EGS). Simulating processes resulting in hydraulic fracturing and/or the remobilization of existing fractures, especially the interaction between propagating fractures and existing fractures, represents a critical goal of our project. To this end, we are continuing to develop a hydraulic fracturing simulation capability within the Livermore Distinct Element Code (LDEC), a combined FEM/DEM analysis code with explicit solid-fluid mechanics coupling. LDEC simulations start from an initial fracture distribution which can be stochastically generated or upscaled from the statistics of an actual fracture distribution. During the hydraulic stimulation process, LDEC tracks the propagation of fractures and other modifications to the fracture system. The output is transferred to the Non-isothermal Unsaturated Flow and Transport (NUFT) code to capture heat transfer and flow at the reservoir scale. This approach is intended to offer flexibility in the types of analyses we can perform, including evaluating the effects of different system heterogeneities on the heat extraction rate as well as seismicity associated with geothermal operations. This paper details the basic methodology of our approach. Two numerical examples showing the capability and effectiveness of our simulator are also presented.
Discrete fracture modeling-ESF North Portal Area, Yucca Mountain Nevada
Anna, Lawrence O.
1995-01-01
Statistical parameters from three-dimensional fracture network and hydraulic parameters were developed to be used in site scale models. This approach utilizes geometric fracture models and assess their impact on flow characteristics and parameters. Laboratory and field-testing data will be integrated to calibrate the flow models and to determine sensitivities of the system.
Rain water filtering down through the soil will provide recharge of the saturated fractured rock aquifer. he computer model FRACNET has been designed to distribute areal recharge into linear fracture zones in order to complete the regional water balance. n this presentation, a te...
Hyman, Jeffrey De'Haven; Painter, S. L.; Viswanathan, H.; Makedonska, N.; Karra, S.
2015-09-12
We investigate how the choice of injection mode impacts transport properties in kilometer-scale three-dimensional discrete fracture networks (DFN). The choice of injection mode, resident and flux-weighted, is designed to mimic different physical phenomena. It has been hypothesized that solute plumes injected under resident conditions evolve to behave similarly to solutes injected under flux-weighted conditions. Previously, computational limitations have prohibited the large-scale simulations required to investigate this hypothesis. We investigate this hypothesis by using a high-performance DFN suite, dfnWorks, to simulate flow in kilometer-scale three-dimensional DFNs based on fractured granite at the Forsmark site in Sweden, and adopt a Lagrangian approach to simulate transport therein. Results show that after traveling through a pre-equilibrium region, both injection methods exhibit linear scaling of the first moment of travel time and power law scaling of the breakthrough curve with similar exponents, slightly larger than 2. Lastly, the physical mechanisms behind this evolution appear to be the combination of in-network channeling of mass into larger fractures, which offer reduced resistance to flow, and in-fracture channeling, which results from the topology of the DFN.
Hyman, Jeffrey De'Haven; Painter, S. L.; Viswanathan, H.; Makedonska, N.; Karra, S.
2015-09-12
We investigate how the choice of injection mode impacts transport properties in kilometer-scale three-dimensional discrete fracture networks (DFN). The choice of injection mode, resident and flux-weighted, is designed to mimic different physical phenomena. It has been hypothesized that solute plumes injected under resident conditions evolve to behave similarly to solutes injected under flux-weighted conditions. Previously, computational limitations have prohibited the large-scale simulations required to investigate this hypothesis. We investigate this hypothesis by using a high-performance DFN suite, dfnWorks, to simulate flow in kilometer-scale three-dimensional DFNs based on fractured granite at the Forsmark site in Sweden, and adopt a Lagrangian approachmore » to simulate transport therein. Results show that after traveling through a pre-equilibrium region, both injection methods exhibit linear scaling of the first moment of travel time and power law scaling of the breakthrough curve with similar exponents, slightly larger than 2. Lastly, the physical mechanisms behind this evolution appear to be the combination of in-network channeling of mass into larger fractures, which offer reduced resistance to flow, and in-fracture channeling, which results from the topology of the DFN.« less
Hyman, Jeffrey; Painter, Scott L; Viswanathan, Hari; Makedonska, Nataliia; Karra, Satish
2015-01-01
We investigate how the choice of injection mode impacts transport properties in kilometer-scale three-dimensional discrete fracture networks (DFN). The choice of injection mode, resident or flux-weighted, is designed to mimic different physical phenomena. It has been hypothesized that solute plumes injected under resident conditions evolve to behave similarly to solutes injected under flux-weighted conditions. Previously, computational limitations have prohibited the large scale simulations required to investigate this hypothesis. We investigate this hypothesis by using a high performance DFN suite, dfnWorks, to simulate flow in kilometer-scale three-dimensional DFNs based on fractured granite at the Forsmark site in Sweden, and adopt amore » Lagrangian approach to simulate transport therein. Results show that a after traveling through a pre-equilibrium region both injection methods exhibit linear scaling of the first moment of travel time and power law scaling of the breakthrough curve with similar exponents, slightly larger than two. The physical mechanisms behind this evolution appear to be the combination of in-network channeling of mass into larger fractures which o er reduced resistance to flow and in-fracture channeling that results from the topology of the DFN.« less
Hyman, Jeffrey; Painter, Scott L; Viswanathan, Hari; Makedonska, Nataliia; Karra, Satish
2015-01-01
We investigate how the choice of injection mode impacts transport properties in kilometer-scale three-dimensional discrete fracture networks (DFN). The choice of injection mode, resident or flux-weighted, is designed to mimic different physical phenomena. It has been hypothesized that solute plumes injected under resident conditions evolve to behave similarly to solutes injected under flux-weighted conditions. Previously, computational limitations have prohibited the large scale simulations required to investigate this hypothesis. We investigate this hypothesis by using a high performance DFN suite, dfnWorks, to simulate flow in kilometer-scale three-dimensional DFNs based on fractured granite at the Forsmark site in Sweden, and adopt a Lagrangian approach to simulate transport therein. Results show that a after traveling through a pre-equilibrium region both injection methods exhibit linear scaling of the first moment of travel time and power law scaling of the breakthrough curve with similar exponents, slightly larger than two. The physical mechanisms behind this evolution appear to be the combination of in-network channeling of mass into larger fractures which o er reduced resistance to flow and in-fracture channeling that results from the topology of the DFN.
Stafford, P.L.
1996-05-01
Simulations of a tritium tracer experiment in fractured shale saprolite, conducted at the Oak Ridge National Laboratory, were performed using 1D and 2D equivalent porous medium (EPM) and discrete-fracture/matrix-diffusion (DFMD) models. The models successfully reproduced the general shape of the breakthrough curves in down-gradient monitoring wells which are characterized by rapid first arrival, a slow-moving center of mass, and a persistent ``tail`` of low concentration. In plan view, the plume shows a large degree of transverse spreading with the width almost as great as the length. EPM models were sensitive to dispersivity coefficient values which had to be large (relative to the 3.7m distance between the injection and monitoring wells) to fit the tail and transverse spreading. For example, to fit the tail a longitudinal dispersivity coefficient, {alpha}{sub L}, of 0.8 meters for the 2D simulations was used. To fit the transverse spreading, a transverse dispersivity coefficient, {alpha}{sub T}, of 0.8 to 0.08 meters was used indicating an {alpha}{sub L}/{alpha}{sub T} ratio between 10 and 1. Transverse spreading trends were also simulated using a 2D DFMD model using a few larger aperture fractures superimposed onto an EPM. Of the fracture networks studied, only those with truncated fractures caused transverse spreading. Simulated tritium levels in all of the cases were larger than observed values by a factor of approximately 100. Although this is partly due to input of too much tritium mass by the models it appears that dilution in the wells, which were not purged prior to sampling, is also a significant factor. The 1D and 2D EPM models were fitted to monitoring data from the first five years of the experiment and then used to predict future tritium concentrations.
The Evolution of Fracture Systems in Rocks with Veins: Insights from 3D Discrete Element Models
NASA Astrophysics Data System (ADS)
Virgo, S.; Urai, J. L.; Abe, S.
2014-12-01
Observations from natural vein systems suggest that preexisting veins can strongly influence orientation, continuity and connectivity of fractures in a rock even in cases where the orientation of the veins is incompatible with the orientation of the stress field. We present a numerical method to model cycles of fracturing and sealing in a rotating stress field to simulate such systems, for different strength ratios of host rock and vein. We study a layered model under vertical stress and uniaxial horizontal extension. This represents common conditions in sedimentary basins with layers of varying composition. The model with fractures that form during the first deformation phase is sealed and deformed again in a different direction to model the effect of a changing horizontal stress field. We find different types of fracture interaction with veins, depending on the strength contrast between veins and host rock and amount of rotation. The crack-seal and crack-jump mechanisms ensue naturally from the models as a result of the strength of the vein material relative to the host rock. Weak veins localize fracturing and reactivate, even in high misorientation to the extension direction. Connecting fractures between reactivated veins form at a higher angle to the veins than expected. In these systems, the connectivity of the fracture network is dramatically increased. Veins stronger than the host rock have less influence on the new fractures. Most fractures crosscut the veins by the step-over mechanism. Deflection occurs for favorable vein orientations but the deflection length is very short. The results are in good agreement with natural crack seal vein networks found in carbonate rocks of the Oman Mountains. We find that preexisting veins can change the fracture behavior of a rock in a way that new fractures do not necessarily align with the principle extension direction and form a highly connected network with reactivated veins that dramatically enhances lateral
NASA Astrophysics Data System (ADS)
Gable, C. W.; Hyman, J.; Karra, S.; Makedonska, N.; Painter, S. L.; Viswanathan, H. S.
2015-12-01
dfnWorks generates discrete fracture networks (DFN) of planar polygons, creates a high quality conforming Delaunay triangulation of the intersecting DFN polygons, assigns properties (aperture, permeability) using geostatistics, sets boundary and initial conditions, solves pressure/flow in single or multi-phase fluids (water, air, CO2) using the parallel PFLOTRAN or serial FEHM, and solves for transport using Lagrangian particle tracking. We outline the dfnWorks workflow and present applications from a range of fractured rock systems. dfnWorks (http://www.lanl.gov/expertise/teams/view/dfnworks) is composed of three main components, all of which are freely available. dfnGen generates a distribution of fracture polygons from site characterization data (statistics or deterministic fractures) and utilizes the FRAM (Feature Rejection Algorithm for Meshing) to guarantee the mesh generation package LaGriT (lagrit.lanl.gov) will generate a high quality conforming Delaunay triangular mesh. dfnWorks links the mesh to either PFLOTRAN (pflotran.org) or FEHM (fehm.lanl.gov) for solving flow and transport. The various physics options available in FEHM and PFLOTRAN such as single and multi-phase flow and reactive transport are all available with appropriate initial and boundary conditions and material property models. dfnTrans utilizes explicit Lagrangian particle tracking on the DFN using a velocity field reconstructed from the steady state pressure/flow field solution obtained in PFLOTRAN or FEHM. Applications are demonstrated for nuclear waste repository in fractured granite, CO2 sequestration and extraction of unconventional hydrocarbon resources.
The Combined Finite-Discrete Element Method applied to the Study of Rock Fracturing Behavior in 3D
Rougier, Esteban; Bradley, Christopher R.; Broom, Scott T.; Knight, Earl E.; Munjiza, Ante; Sussman, Aviva J.; Swift, Robert P.
2011-01-01
Since its introduction the combined finite-discrete element method (FEM/DEM), has become an excellent tool to address a wide range of problems involving fracturing and fragmentation of solids. Within the context of rock mechanics, the FEM/DEM method has been applied to many complex industrial problems such as block caving, deep mining techniques, rock blasting, seismic waves, packing problems, rock crushing problems, etc. In the real world most of the problems involving fracture and fragmentation of solids are three dimensional problems. With the aim of addressing these problems an improved 2D/3D FEM/DEM capability has been developed at Los Alamos National Laboratory (LANL). These capabilities include state of the art 3D contact detection, contact interaction, constitutive material models, and fracture models. In this paper, Split Hopkinson Pressure Bar (SHPB) Brazilian experiments are simulated using this improved 2D/3D FEM/DEM approach which is implemented in LANL's MUNROU (Munjiza-Rougier) code. The results presented in this work show excellent agreement with both the SHPB experiments and previous 2D numerical simulations performed by other FEM/DEM research groups.
NASA Astrophysics Data System (ADS)
Cvetkovic, V.; Frampton, A.; Painter, S.; Selroos, J.
2008-12-01
An important challenge in subsurface hydrology is predictive modeling of tracer transport in sparsely fractured rock. A particular issue relevant for applications is how to accurately account for retention processes that are due to exchange (diffusion-sorption) of tracers between mobile fluid in fractures and immobile fluid in the rock matrix. Typically, tracers are subject to decay processes which may involve chains and in-growth (e.g., for radionuclides and some classes of hydrocarbons). Recently, a comprehensive particle-based methodology for upscaling transport with emphasis on tracer retention has been presented and applied to stochastic 2D discrete fracture networks (Frampton and Cvetkovic 2007, WRR, 43, W10429). Furthermore, a time domain random walk method has also recently been presented that effectively accounts for different exchange mechanisms and in-growth (Painter et al. 2008, WRR, 44, W01406). Now we present further advances in coupling these novel methodologies for solving radionuclide transport, and apply them to realistic 3D fracture networks, based on comprehensive data sets obtained from site characterization of the Laxemar area in south-east Sweden. Site measurements have revealed at least five fracture sets based on statistically significant orientation data, exhibiting power-law behaviour for fracture size and inferred transmissivity distributions. A few equally probable DFN realizations are generated based on these interpretations of the field data, in which advective fluid flow is solved using boundary conditions that mimic natural conditions. Thereafter, many particles are injected and tracked through the system, providing first- passage distributions of particle residence time and of the transport resistance parameter (quantifying the hydrodynamic control of retention). These distributions are then used as a basis for implementing the particle time-domain random walk model for radionuclide transport with retention and in-growth. Also, an
Hai Huang; Ben Spencer; Jason Hales
2014-10-01
A discrete element Model (DEM) representation of coupled solid mechanics/fracturing and heat conduction processes has been developed and applied to explicitly simulate the random initiations and subsequent propagations of interacting thermal cracks in a ceramic nuclear fuel pellet during initial rise to power and during power cycles. The DEM model clearly predicts realistic early-life crack patterns including both radial cracks and circumferential cracks. Simulation results clearly demonstrate the formation of radial cracks during the initial power rise, and formation of circumferential cracks as the power is ramped down. In these simulations, additional early-life power cycles do not lead to the formation of new thermal cracks. They do, however clearly indicate changes in the apertures of thermal cracks during later power cycles due to thermal expansion and shrinkage. The number of radial cracks increases with increasing power, which is consistent with the experimental observations.
2014-01-01
Background Osteoporotic fractures cause a large health burden and substantial costs. This study estimated the expected fracture numbers and costs for the remaining lifetime of postmenopausal women in Germany. Methods A discrete event simulation (DES) model which tracks changes in fracture risk due to osteoporosis, a previous fracture or institutionalization in a nursing home was developed. Expected lifetime fracture numbers and costs per capita were estimated for postmenopausal women (aged 50 and older) at average osteoporosis risk (AOR) and for those never suffering from osteoporosis. Direct and indirect costs were modeled. Deterministic univariate and probabilistic sensitivity analyses were conducted. Results The expected fracture numbers over the remaining lifetime of a 50 year old woman with AOR for each fracture type (% attributable to osteoporosis) were: hip 0.282 (57.9%), wrist 0.229 (18.2%), clinical vertebral 0.206 (39.2%), humerus 0.147 (43.5%), pelvis 0.105 (47.5%), and other femur 0.033 (52.1%). Expected discounted fracture lifetime costs (excess cost attributable to osteoporosis) per 50 year old woman with AOR amounted to €4,479 (€1,995). Most costs were accrued in the hospital €1,743 (€751) and long-term care sectors €1,210 (€620). Univariate sensitivity analysis resulted in percentage changes between -48.4% (if fracture rates decreased by 2% per year) and +83.5% (if fracture rates increased by 2% per year) compared to base case excess costs. Costs for women with osteoporosis were about 3.3 times of those never getting osteoporosis (€7,463 vs. €2,247), and were markedly increased for women with a previous fracture. Conclusion The results of this study indicate that osteoporosis causes a substantial share of fracture costs in postmenopausal women, which strongly increase with age and previous fractures. PMID:24981316
NASA Astrophysics Data System (ADS)
Follin, Sven; Hartley, Lee; Rhén, Ingvar; Jackson, Peter; Joyce, Steven; Roberts, David; Swift, Ben
2014-03-01
The large-scale geological structure of the crystalline rock at the proposed high-level nuclear waste repository site at Forsmark, Sweden, has been classified in terms of deformation zones of elevated fracture frequency. The rock between deformation zones was divided into fracture domains according to fracture frequency. A methodology to constrain the geometric and hydraulic parameters that define a discrete fracture network (DFN) model for each fracture domain is presented. The methodology is based on flow logging and down-hole imaging in cored boreholes in combination with DFN realizations, fracture connectivity analysis and pumping test simulations. The simulations suggest that a good match could be obtained for a power law size distribution where the value of the location parameter equals the borehole radius but with different values for the shape parameter, depending on fracture domain and fracture set. Fractures around 10-100 m in size are the ones that typically form the connected network, giving inflows in the simulations. The report also addresses the issue of up-scaling of DFN properties to equivalent continuous porous medium (ECPM) bulk flow properties. Comparisons with double-packer injection tests provide confidence that the derived DFN formulation of detailed flows within individual fractures is also suited to simulating mean bulk flow properties and their spatial variability.
NASA Astrophysics Data System (ADS)
Weisbrod, N.; Tran, E. L.; Klein-BenDavid, O.; Teutsch, N.
2015-12-01
Geological disposal of high-level radioactive waste is the long term solution for the disposal of long lived radionuclides and spent fuel. However, some radionuclides might be released from these repositories into the subsurface as a result of leakage, which ultimately make their way into groundwater. Engineered bentonite barriers around nuclear waste repositories are generally considered sufficient to impede the transport of radionuclides from their source to the groundwater. However, colloidal-sized mobile bentonite particles ("carrier" colloids) originating from these barriers have come under investigation as a potential transport vector for radionuclides sorbed to them. As lanthanides are generally accepted to have the same chemical behaviors as their more toxic actinide counterparts, lanthanides are considered an acceptable substitute for research on radionuclide transportation. This study aims to evaluate the transport behaviors of lanthanides in colloid-facilitated transport through a fractured chalk matrix and under geochemical conditions representative the Negev desert, Israel. The migration of Ce both with and without colloidal particles was explored and compared to the migration of a conservative tracer (bromide) using a flow system constructed around a naturally fractured chalk core. Results suggest that mobility of Ce as a solute is negligible. In experiments conducted without bentonite colloids, the 1% of the Ce that was recovered migrated as "intrinsic" colloids in the form of carbonate precipitates. However, the total recovery of the Ce increased to 9% when it was injected into the core in the presence of bentonite colloids and 13% when both bentonite and precipitate colloids were injected. This indicates that lanthanides are essentially immobile in chalk as a solute but may be mobile as carbonate precipitates. Bentonite colloids, however, markedly increase the mobility of lanthanides through fractured chalk matrices.
NASA Astrophysics Data System (ADS)
Virgo, Simon; Abe, Steffen; Urai, Janos L.
2016-03-01
We present the results of a comparative study of loading conditions on the interactions between extension fractures and veins. We model the fracture behavior of brittle discrete element materials each containing a tabular vein body of variable orientation and strength in two different loading conditions. The first is uniaxial tension, applied with servo-controlled sidewalls. The second is a boudinage boundary condition in which a tensile triaxial stress state is induced in the brittle model volume by quasi-viscous extensional deformation in the adjacent layers. Most of the fracture- vein interactions observed in uniaxial tension also exists in boudinage boundary conditions. However, the importance of each interaction mechanism for a given configuration of relative strength and misorientation of the vein may differ according to the loading mechanism. Nucleation and internal deflection is under both boundary conditions the dominating fracture-vein interaction style in weak veins. In uniaxial tension models, strong veins tend to alter the fracture path by external deflection, while under boudinage loading these veins are more likely overcome by the fracture step over mechanism. Dynamic bifurcation of fractures was observed in uniaxial tension models but never for boudinage boundary conditions. This is because the acceleration of fracture tips in these conditions is suppressed by interaction with distributed fractures as well as viscous damping by the neighboring layers.
NASA Astrophysics Data System (ADS)
Hodgetts, David; Seers, Thomas
2015-04-01
-deterministic, outcrop constrained discrete fracture network modeling code to derive volumetric fault intensity measures (fault area per unit volume / fault volume per unit volume). Producing per-vertex measures of volumetric intensity; our method captures the spatial variability in 3D fault density across a surveyed outcrop, enabling first order controls to be probed. We demonstrate our approach on pervasively faulted exposures of a Permian aged reservoir analogue from the Vale of Eden Basin, UK.
NASA Astrophysics Data System (ADS)
Novakowski, K. S.
2015-12-01
The development of conceptual models for solute migration in discrete fracture networks has typically been based on a combination of core logs, borehole geophysics, and some form of single-well hydraulic test using discrete zones. More rarely, interwell hydraulic tests and interwell tracer experiments are utilised to directly explore potential transport pathways. The latter methods are less widely employed simply due to potentially significant increases in the cost and effort in site characterization. To date however there is a paucity of literature comparing the efficacy of the standard procedure with what should be more definitive identification of transport pathways using interwell methods. In the present study, a detailed comparison is conducted by developing conceptual models from three separate data sets, the first based on core logs, geology and single-well hydraulic tests, the second based on a large suite of pulse interference tests, and the third based on a series of radially-divergent and injection-withdrawal tracer experiments. The study was conducted in an array of five HQ-sized wells, 28-32 m in depth and arranged in a five star pattern, 10 m on a side. The wells penetrate the contact between a Cambrian-aged limestone, and underlying Precambrian gneiss. The core was logged for potentially open fractures using a ranking system, and 87 contiguous hydraulic tests were conducted using a 0.85-m packer spacing. A total of 57 pulse interference tests were conducted using two wells as injection points, and 11 tracer experiments were conducted using either sample collection or in-situ detection via a submersible fluorometer. The results showed very distinct conceptual models depending on the data set, with the model based on the single-well testing significantly over-predicting the number and connection of solute transport pathways. The results of the pulse interference tests also over predict the transport pathways, but to a lesser degree. Quantification of
NASA Astrophysics Data System (ADS)
Watanabe, N.; Wong, L.; Bloecher, G.; Cacace, M.; Kolditz, O.
2012-12-01
We present our recent development of the finite element method (FEM) for simulating coupled thermo-hydro-mechanical (THM) processes in discretely fractured porous media and an application to geothermal reservoir modeling for the research test site Gross Schoenebeck in Germany operated by the GFZ German Research Centre for Geosciences. Numerical analysis of multi-physics problems in fractured rocks is important for various geotechnical applications. In particular for enhanced geothermal reservoirs where induced fractures and possibly natural fault systems dominate the system behavior, explicit modeling of those characteristic fractures (i.e. discrete fracture models) is essential to get more detailed understanding of in-situ processes and reliable estimations of heat extraction from those deep reservoirs. However, as fractures are mechanical discontinuities, it is difficult to solve the problems using continuity based numerical methods such as the FEM. Currently, equivalent porous medium or multiple continuum model approaches are often only the way to model fractured rocks with the FEM. The authors have recently developed lower-dimensional interface elements (LIEs) for modeling mechanics-involved coupled processes with pre-existing fractures (Watanabe et al. 2012 IJNME). The method does not require any double nodes unlike conventional interface elements. Moreover, for coupled problems, the approach allows for the use of a single mesh for both mechanical and other related processes such as flow and transport. All the code developments have been carried out within the scientific open source project OpenGeoSys (www.opengeosys.net) (Kolditz et al. 2012 EES). Using both traditional and new simulation techniques, a geothermal reservoir model for the research test site Gross Schoenebeck has been developed. Unstructured meshing of the complex faulted reservoir including both rock matrix and fracture elements has been conducted using recently developed automatic
Hall, Michael C.
1963-01-01
Recent studies on the epidemiology and repair of fractures are reviewed. The type and severity of the fracture bears a relation to the age, sex and occupation of the patient. Bone tissue after fracture shows a process of inflammation and repair common to all members of the connective tissue family, but it repairs with specific tissue. Cartilage forms when the oxygen supply is outgrown. After a fracture, the vascular bed enlarges. The major blood supply to healing tissue is from medullary vessels and destruction of them will cause necrosis of the inner two-thirds of the cortex. Callus rapidly mineralizes, but full mineralization is achieved slowly; increased mineral metabolism lasts several years after fracture. PMID:13952119
... commonly happen because of car accidents, falls, or sports injuries. Other causes are low bone density and osteoporosis, which cause weakening of the bones. Overuse can cause stress fractures, which are very small cracks in the ...
A fracture is a break, usually in a bone. If the broken bone punctures the skin, it is called an open ... falls, or sports injuries. Other causes are low bone density and osteoporosis, which cause weakening of the ...
Teleportation of a 3-dimensional GHZ State
NASA Astrophysics Data System (ADS)
Cao, Hai-Jing; Wang, Huai-Sheng; Li, Peng-Fei; Song, He-Shan
2012-05-01
The process of teleportation of a completely unknown 3-dimensional GHZ state is considered. Three maximally entangled 3-dimensional Bell states function as quantum channel in the scheme. This teleportation scheme can be directly generalized to teleport an unknown d-dimensional GHZ state.
NASA Astrophysics Data System (ADS)
Bruel, D.; Baujard, C.
2005-05-01
Heat extraction from deep engineered fractured formations is currently under investigation at the Soultz sous Forêts site with the support of the European Commission. The challenge is to develop a reservoir at great depth and to circulate a fluid in order to recover heat and produce electricity. The pilot project evolved toward a three well system at 5 km in depth with temperatures close to 200 C. Massive hydraulic tests performed to develop the reservoir have shown from the recorded micro-seismic signature, that fractures can easily be re-activated. The discussion now focusses on the hydraulic significance of the shear failure mechanism, considered as the source of the accoustic emissions. To improve our understanding of these coupled hydrau-mechanical processes, a numerical model was presented [1], based on a 3D random description of fracture networks. Local flow rules along equivalent 1D channels connecting the fractures can account for (i) a normal closure versus effective stress law together with (ii) a dilatant behaviour during shearing motion when a Mohr-Coulomb failure criterion is met. The purpose of the present work is to simulate injection tests in some synthetic fracture network using power law distributions for the fracture size, and to analyse the spatio-temporal growth of the sheared zones. Assuming that this process is analogue to the triggering of the microseismicity, we then perform an evaluation of the so called SBRC reservoir characterisation method [2] stating that the spatial position of the triggering front in an homogeneous isotropic poroelastic medium with a hydraulic diffusivity Dh is at time t given by √4 π Dh t. We conclude to its validity, although it is found sensitive to the hypothesis of critically stressed pre-existing fractures. The connectivity of the sub-set of subcritically oriented fractures plays a major role in the succes of a stimulation treatment and controls an equivalent macro-cohesion behaviour at the reservoir scale
Gureghian, A.B.
1990-08-01
Analytical solutions based on the Laplace transforms are presented for the one-dimensional, transient, advective-dispersive transport of a reacting radionuclide through a discrete planar fracture with constant aperture subject to diffusion in the surrounding rock matrix where both regions of solute migration display residual concentrations. The dispersion-free solutions, which are of closed form, are also reported. The solution assumes that the ground-water flow regime is under steady-state and isothermal conditions and that the rock matrix is homogeneous, isotropic, and saturated with stagnant water. The verification of the solution was performed by means of related analytical solutions dealing with particular aspects of the transport problem under investigation on the one hand, and a numerical solution capable of handling the complete problem on the other. The integrals encountered in the general solution are evaluated by means of a composite Gauss-Legendre quadrature scheme. 9 refs., 8 figs., 32 tabs.
NASA Astrophysics Data System (ADS)
Selroos, J. O.; Appleyard, P.; Bym, T.; Follin, S.; Hartley, L.; Joyce, S.; Munier, R.
2015-12-01
In 2011 the Swedish Nuclear Fuel and Waste Management Company (SKB) applied for a license to start construction of a final repository for spent nuclear fuel at Forsmark in Northern Uppland, Sweden. The repository is to be built at approximately 500 m depth in crystalline rock. A stochastic, discrete fracture network (DFN) concept was chosen for interpreting the surface-based (incl. boreholes) data, and for assessing the safety of the repository in terms of groundwater flow and flow pathways to and from the repository. Once repository construction starts, also underground data such as tunnel pilot borehole and tunnel trace data will become available. It is deemed crucial that DFN models developed at this stage honors the mapped structures both in terms of location and geometry, and in terms of flow characteristics. The originally fully stochastic models will thus increase determinism towards the repository. Applying the adopted probabilistic framework, predictive modeling to support acceptance criteria for layout and disposal can be performed with the goal of minimizing risks associated with the repository. This presentation describes and illustrates various methodologies that have been developed to condition stochastic realizations of fracture networks around underground openings using borehole and tunnel trace data, as well as using hydraulic measurements of inflows or hydraulic interference tests. The methodologies, implemented in the numerical simulators ConnectFlow and FracMan/MAFIC, are described in some detail, and verification tests and realistic example cases are shown. Specifically, geometric and hydraulic data are obtained from numerical synthetic realities approximating Forsmark conditions, and are used to test the constraining power of the developed methodologies by conditioning unconditional DFN simulations following the same underlying fracture network statistics. Various metrics are developed to assess how well the conditional simulations compare to
3-Dimensional Topographic Models for the Classroom
NASA Technical Reports Server (NTRS)
Keller, J. W.; Roark, J. H.; Sakimoto, S. E. H.; Stockman, S.; Frey, H. V.
2003-01-01
We have recently undertaken a program to develop educational tools using 3-dimensional solid models of digital elevation data acquired by the Mars Orbital Laser Altimeter (MOLA) for Mars as well as a variety of sources for elevation data of the Earth. This work is made possible by the use of rapid prototyping technology to construct solid 3-Dimensional models of science data. We recently acquired rapid prototyping machine that builds 3-dimensional models in extruded plastic. While the machine was acquired to assist in the design and development of scientific instruments and hardware, it is also fully capable of producing models of spacecraft remote sensing data. We have demonstrated this by using Mars Orbiter Laser Altimeter (MOLA) topographic data and Earth based topographic data to produce extruded plastic topographic models which are visually appealing and instantly engage those who handle them.
Fuller, Sam M; Butz, Daniel R; Vevang, Curt B; Makhlouf, Mansour V
2014-09-01
Three-dimensional printing is being rapidly incorporated in the medical field to produce external prosthetics for improved cosmesis and fabricated molds to aid in presurgical planning. Biomedically engineered products from 3-dimensional printers are also utilized as implantable devices for knee arthroplasty, airway orthoses, and other surgical procedures. Although at first expensive and conceptually difficult to construct, 3-dimensional printing is now becoming more affordable and widely accessible. In hand surgery, like many other specialties, new or customized instruments would be desirable; however, the overall production cost restricts their development. We are presenting our step-by-step experience in creating a bone reduction clamp for finger fractures using 3-dimensional printing technology. Using free, downloadable software, a 3-dimensional model of a bone reduction clamp for hand fractures was created based on the senior author's (M.V.M.) specific design, previous experience, and preferences for fracture fixation. Once deemed satisfactory, the computer files were sent to a 3-dimensional printing company for the production of the prototypes. Multiple plastic prototypes were made and adjusted, affording a fast, low-cost working model of the proposed clamp. Once a workable design was obtained, a printing company produced the surgical clamp prototype directly from the 3-dimensional model represented in the computer files. This prototype was used in the operating room, meeting the expectations of the surgeon. Three-dimensional printing is affordable and offers the benefits of reducing production time and nurturing innovations in hand surgery. This article presents a step-by-step description of our design process using online software programs and 3-dimensional printing services. As medical technology advances, it is important that hand surgeons remain aware of available resources, are knowledgeable about how the process works, and are able to take advantage of
3-dimensional imaging at nanometer resolutions
Werner, James H.; Goodwin, Peter M.; Shreve, Andrew P.
2010-03-09
An apparatus and method for enabling precise, 3-dimensional, photoactivation localization microscopy (PALM) using selective, two-photon activation of fluorophores in a single z-slice of a sample in cooperation with time-gated imaging for reducing the background radiation from other image planes to levels suitable for single-molecule detection and spatial location, are described.
NASA Astrophysics Data System (ADS)
Jackson, C.; Hartley, L. J.; Hoch, A.; Holton, D.; Hunter, F. M.; McCarthy, R.; Marsic, N.; Gylling, B.
2006-12-01
The Swedish Nuclear Fuel and Waste Management Company (SKB) is carrying out site investigations in two different areas in Sweden with the objective of describing the in-situ conditions for a deep rock repository for spent nuclear fuel. The two candidate areas are Forsmark and Laxemar, both located on the east coast of Sweden. An important aspect of site investigations is to develop and demonstrate an understanding of groundwater flow and solute transport. Since the geology in both candidate areas is comprised of hard crystalline rocks, the groundwater flow is predominantly contained within fractures, and therefore a discrete fracture network (DFN) concept has been applied to describe and model the hydrogeological situation at the sites. Much observed field data from several different disciplines (geology, rock mechanics, geophysics, hydrogeology and hydrochemistry) has been acquired from the sites, including from several deep cored boreholes, to inform an overall description. Many aspects of the site description are brought together in constructing a regional scale hydrogeological model to integrate the concepts and data interpretations, which are then tested against a range of field observations to build confidence that the models are representative. A methodology has been developed based on assembling a regional hydrogeological model from three main components: hydraulic conductor domains (HCD) that represent deterministic large scale deformation zones; hydraulic rock domains (HRD) that use a stochastic DFN model to represent the background rock between the deformation zones; and hydraulic soil domains (HSD) that represent near-surface Quaternary deposits. The HCD are interpreted from geophysical methods, drilling and single-hole hydraulic tests. For the HRD, borehole image- and core-logs, outcrop maps, and short-interval flow-logging are integrated to parameterise a DFN model for specific hydrogeological rock domains. Geological information, statistical analysis
Biochemical Applications Of 3-Dimensional Fluorescence Spectrometry
NASA Astrophysics Data System (ADS)
Leiner, Marc J.; Wolfbeis, Otto S.
1988-06-01
We investigated the 3-dimensional fluorescence of complex mixtures of bioloquids such as human serum, serum ultrafiltrate, human urine, and human plasma low density lipoproteins. The total fluorescence of human serum can be divided into a few peaks. When comparing fluorescence topograms of sera, from normal and cancerous subjects, we found significant differences in tryptophan fluorescence. Although the total fluorescence of human urine can be resolved into 3-5 distinct peaks, some of them. do not result from single fluorescent urinary metabolites, but rather from. several species having similar spectral properties. Human plasma, low density lipoproteins possess a native fluorescence that changes when submitted to in-vitro autoxidation. The 3-dimensional fluorescence demonstrated the presence of 7 fluorophores in the lipid domain, and 6 fluorophores in the protein. dovain- The above results demonstrated that 3-dimensional fluorescence can resolve the spectral properties of complex ,lxtures much better than other methods. Moreover, other parameters than excitation and emission wavelength and intensity (for instance fluorescence lifetime, polarization, or quenchability) may be exploited to give a multidl,ensio,a1 matrix, that is unique for each sample. Consequently, 3-dimensio:Hhal fluorescence as such, or in combination with separation techniques is therefore considered to have the potential of becoming a useful new H.ethod in clinical chemistry and analytical biochemistry.
Hydroelectric structures studies using 3-dimensional methods
Harrell, T.R.; Jones, G.V.; Toner, C.K. )
1989-01-01
Deterioration and degradation of aged, hydroelectric project structures can significantly affect the operation and safety of a project. In many cases, hydroelectric headworks (in particular) have complicated geometrical configurations, loading patterns and hence, stress conditions. An accurate study of such structures can be performed using 3-dimensional computer models. 3-D computer models can be used for both stability evaluation and for finite element stress analysis. Computer aided engineering processes facilitate the use of 3-D methods in both pre-processing and post-processing of data. Two actual project examples are used to emphasize the authors' points.
Theory of relativistic Brownian motion: the (1+3) -dimensional case.
Dunkel, Jörn; Hänggi, Peter
2005-09-01
A theory for (1+3) -dimensional relativistic Brownian motion under the influence of external force fields is put forward. Starting out from a set of relativistically covariant, but multiplicative Langevin equations we describe the relativistic stochastic dynamics of a forced Brownian particle. The corresponding Fokker-Planck equations are studied in the laboratory frame coordinates. In particular, the stochastic integration prescription--i.e., the discretization rule dilemma--is elucidated (prepoint discretization rule versus midpoint discretization rule versus postpoint discretization rule). Remarkably, within our relativistic scheme we find that the postpoint rule (or the transport form) yields the only Fokker-Planck dynamics from which the relativistic Maxwell-Boltzmann statistics is recovered as the stationary solution. The relativistic velocity effects become distinctly more pronounced by going from one to three spatial dimensions. Moreover, we present numerical results for the asymptotic mean-square displacement of a free relativistic Brownian particle moving in 1+3 dimensions. PMID:16241514
3-dimensional bioprinting for tissue engineering applications.
Gu, Bon Kang; Choi, Dong Jin; Park, Sang Jun; Kim, Min Sup; Kang, Chang Mo; Kim, Chun-Ho
2016-01-01
The 3-dimensional (3D) printing technologies, referred to as additive manufacturing (AM) or rapid prototyping (RP), have acquired reputation over the past few years for art, architectural modeling, lightweight machines, and tissue engineering applications. Among these applications, tissue engineering field using 3D printing has attracted the attention from many researchers. 3D bioprinting has an advantage in the manufacture of a scaffold for tissue engineering applications, because of rapid-fabrication, high-precision, and customized-production, etc. In this review, we will introduce the principles and the current state of the 3D bioprinting methods. Focusing on some of studies that are being current application for biomedical and tissue engineering fields using printed 3D scaffolds. PMID:27114828
Honigmann, Philipp; Thieringer, Florian; Steiger, Regula; Haefeli, Mathias; Schumacher, Ralf; Henning, Julia
2016-03-01
The reconstruction of malunited distal radius fractures is often challenging. Virtual planning techniques and guides for drilling and resection have been used for several years to achieve anatomic reconstruction. These guides have the advantage of leading to better operative results and faster surgery. Here, we describe a technique using a simple implant independent 3-dimensional printed drill guide and template to simplify the surgical reconstruction of a malunited distal radius fracture. PMID:26787406
Vehmeijer, Maarten; van Eijnatten, Maureen; Liberton, Niels; Wolff, Jan
2016-08-01
Fractures of the orbital floor are often a result of traffic accidents or interpersonal violence. To date, numerous materials and methods have been used to reconstruct the orbital floor. However, simple and cost-effective 3-dimensional (3D) printing technologies for the treatment of orbital floor fractures are still sought. This study describes a simple, precise, cost-effective method of treating orbital fractures using 3D printing technologies in combination with autologous bone. Enophthalmos and diplopia developed in a 64-year-old female patient with an orbital floor fracture. A virtual 3D model of the fracture site was generated from computed tomography images of the patient. The fracture was virtually closed using spline interpolation. Furthermore, a virtual individualized mold of the defect site was created, which was manufactured using an inkjet printer. The tangible mold was subsequently used during surgery to sculpture an individualized autologous orbital floor implant. Virtual reconstruction of the orbital floor and the resulting mold enhanced the overall accuracy and efficiency of the surgical procedure. The sculptured autologous orbital floor implant showed an excellent fit in vivo. The combination of virtual planning and 3D printing offers an accurate and cost-effective treatment method for orbital floor fractures. PMID:27137437
Cardiothoracic Applications of 3-dimensional Printing.
Giannopoulos, Andreas A; Steigner, Michael L; George, Elizabeth; Barile, Maria; Hunsaker, Andetta R; Rybicki, Frank J; Mitsouras, Dimitris
2016-09-01
Medical 3-dimensional (3D) printing is emerging as a clinically relevant imaging tool in directing preoperative and intraoperative planning in many surgical specialties and will therefore likely lead to interdisciplinary collaboration between engineers, radiologists, and surgeons. Data from standard imaging modalities such as computed tomography, magnetic resonance imaging, echocardiography, and rotational angiography can be used to fabricate life-sized models of human anatomy and pathology, as well as patient-specific implants and surgical guides. Cardiovascular 3D-printed models can improve diagnosis and allow for advanced preoperative planning. The majority of applications reported involve congenital heart diseases and valvular and great vessels pathologies. Printed models are suitable for planning both surgical and minimally invasive procedures. Added value has been reported toward improving outcomes, minimizing perioperative risk, and developing new procedures such as transcatheter mitral valve replacements. Similarly, thoracic surgeons are using 3D printing to assess invasion of vital structures by tumors and to assist in diagnosis and treatment of upper and lower airway diseases. Anatomic models enable surgeons to assimilate information more quickly than image review, choose the optimal surgical approach, and achieve surgery in a shorter time. Patient-specific 3D-printed implants are beginning to appear and may have significant impact on cosmetic and life-saving procedures in the future. In summary, cardiothoracic 3D printing is rapidly evolving and may be a potential game-changer for surgeons. The imager who is equipped with the tools to apply this new imaging science to cardiothoracic care is thus ideally positioned to innovate in this new emerging imaging modality. PMID:27149367
NASA Astrophysics Data System (ADS)
Wilson, C. E.; Aydin, A.; Durlofsky, L.; Karimi-Fard, M.; Brownlow, D. T.
2008-12-01
An active quarry near Uvalde, TX which mines asphaltic limestone from the Anacacho Formation offers an ideal setting to study fluid-flow in fractured and faulted carbonate rocks. Semi-3D exposures of normal faults and fractures in addition to visual evidence of asphalt concentrations in the quarry help constrain relationships between geologic structures and the flow and transport of hydrocarbons. Furthermore, a subsurface dataset which includes thin sections and measured asphalt concentration from the surrounding region provides a basis to estimate asphalt concentrations and constrain the depositional architecture of both the previously mined portions of the quarry and the un-mined surrounding rock volume. We characterized a series of normal faults and opening mode fractures at the quarry and documented a correlation between the intensity and distribution of these structures with increased concentrations of asphalt. The three-dimensional depositional architecture of the Anacacho Formation was characterized using the subsurface thin sections. Then outcrop exposures of faults, fractured beds, and stratigraphic contacts were mapped and their three-dimensional positions were recorded with differential gps devices. These two datasets were assimilated and a quarry-scale, geologically realistic, three-dimensional Discrete Feature Network (DFN) which represents the geometries and material properties of the matrix, normal faults, and fractures within the quarry was constructed. We then performed two-point flux, control-volume finite- difference fluid-flow simulations with the DFN to investigate the 3D flow and transport of fluids. The results were compared and contrasted with available asphalt concentration estimates from the mine and the aforementioned data from the surrounding drill cores.
Incorporating 3-dimensional models in online articles
Cevidanes, Lucia H. S.; Ruellasa, Antonio C. O.; Jomier, Julien; Nguyen, Tung; Pieper, Steve; Budin, Francois; Styner, Martin; Paniagua, Beatriz
2015-01-01
Introduction The aims of this article were to introduce the capability to view and interact with 3-dimensional (3D) surface models in online publications, and to describe how to prepare surface models for such online 3D visualizations. Methods Three-dimensional image analysis methods include image acquisition, construction of surface models, registration in a common coordinate system, visualization of overlays, and quantification of changes. Cone-beam computed tomography scans were acquired as volumetric images that can be visualized as 3D projected images or used to construct polygonal meshes or surfaces of specific anatomic structures of interest. The anatomic structures of interest in the scans can be labeled with color (3D volumetric label maps), and then the scans are registered in a common coordinate system using a target region as the reference. The registered 3D volumetric label maps can be saved in .obj, .ply, .stl, or .vtk file formats and used for overlays, quantification of differences in each of the 3 planes of space, or color-coded graphic displays of 3D surface distances. Results All registered 3D surface models in this study were saved in .vtk file format and loaded in the Elsevier 3D viewer. In this study, we describe possible ways to visualize the surface models constructed from cone-beam computed tomography images using 2D and 3D figures. The 3D surface models are available in the article’s online version for viewing and downloading using the reader’s software of choice. These 3D graphic displays are represented in the print version as 2D snapshots. Overlays and color-coded distance maps can be displayed using the reader’s software of choice, allowing graphic assessment of the location and direction of changes or morphologic differences relative to the structure of reference. The interpretation of 3D overlays and quantitative color-coded maps requires basic knowledge of 3D image analysis. Conclusions When submitting manuscripts, authors can
Fracture of the nose; Broken nose; Nasal fracture; Nasal bone fracture; Nasal septal fracture ... A fractured nose is the most common fracture of the face. It ... with other fractures of the face. Sometimes a blunt injury can ...
From 2-dimensional cephalograms to 3-dimensional computed tomography scans.
Halazonetis, Demetrios J
2005-05-01
Computed tomography is entering the orthodontic specialty as a mainstream diagnostic modality. Radiation exposure and cost have decreased significantly, and the diagnostic value is very high compared with traditional radiographic options. However, 3-dimensional data present new challenges and need a different approach from traditional viewing of static images to make the most of the available possibilities. Advances in computer hardware and software now enable interactive display of the data on personal computers, with the ability to selectively view soft or hard tissues from any angle. Transfer functions are used to apply transparency and color. Cephalometric measurements can be taken by digitizing points in 3-dimensional coordinates. Application of 3-dimensional data is expected to increase significantly soon and might eventually replace many conventional orthodontic records that are in use today. PMID:15877045
Fracture of the nose; Broken nose; Nasal fracture; Nasal bone fracture; Nasal septal fracture ... A fractured nose is the most common fracture of the face. It usually occurs after an injury and often occurs with ...
NASA Astrophysics Data System (ADS)
Vivaldi, Franco
2015-12-01
The concept of resonance has been instrumental to the study of Hamiltonian systems with divided phase space. One can also define such systems over discrete spaces, which have a finite or countable number of points, but in this new setting the notion of resonance must be re-considered from scratch. I review some recent developments in the area of arithmetic dynamics which outline some salient features of linear and nonlinear stable (elliptic) orbits over a discrete space, and also underline the difficulties that emerge in their analysis.
NASA Astrophysics Data System (ADS)
Vivaldi, Franco
The concept of resonance has been instrumental to the study of Hamiltonian systems with divided phase space. One can also define such systems over discrete spaces, which have a finite or countable number of points, but in this new setting the notion of resonance must be re-considered from scratch. I review some recent developments in the area of arithmetic dynamics which outline some salient features of linear and nonlinear stable (elliptic) orbits over a discrete space, and also underline the difficulties that emerge in their analysis.
Differential Cross Section Kinematics for 3-dimensional Transport Codes
NASA Technical Reports Server (NTRS)
Norbury, John W.; Dick, Frank
2008-01-01
In support of the development of 3-dimensional transport codes, this paper derives the relevant relativistic particle kinematic theory. Formulas are given for invariant, spectral and angular distributions in both the lab (spacecraft) and center of momentum frames, for collisions involving 2, 3 and n - body final states.
Controlled teleportation of a 3-dimensional bipartite quantum state
NASA Astrophysics Data System (ADS)
Cao, Hai-Jing; Chen, Zhong-Hua; Song, He-Shan
2008-07-01
A controlled teleportation scheme of an unknown 3-dimensional (3D) two-particle quantum state is proposed, where a 3D Bell state and 3D GHZ state function as the quantum channel. This teleportation scheme can be directly generalized to teleport an unknown d-dimensional bipartite quantum state.
Basilar skull fracture; Depressed skull fracture; Linear skull fracture ... Skull fractures may occur with head injuries . The skull provides good protection for the brain. However, a severe impact ...
RESEARCH PROGRAM ON FRACTURED PETROLEUM RESERVOIRS
Abbas Firoozabadi
2002-04-12
Numerical simulation of water injection in discrete fractured media with capillary pressure is a challenge. Dual-porosity models in view of their strength and simplicity can be mainly used for sugar-cube representation of fractured media. In such a representation, the transfer function between the fracture and the matrix block can be readily calculated for water-wet media. For a mixed-wet system, the evaluation of the transfer function becomes complicated due to the effect of gravity. In this work, they use a discrete-fracture model in which the fractures are discretized as one dimensional entities to account for fracture thickness by an integral form of the flow equations. This simple step greatly improves the numerical solution. Then the discrete-fracture model is implemented using a Galerkin finite element method. The robustness and the accuracy of the approach are shown through several examples. First they consider a single fracture in a rock matrix and compare the results of the discrete-fracture model with a single-porosity model. Then, they use the discrete-fracture model in more complex configurations. Numerical simulations are carried out in water-wet media as well as in mixed-wet media to study the effect of matrix and fracture capillary pressures.
Multimodality 3-Dimensional Image Integration for Congenital Cardiac Catheterization
2014-01-01
Cardiac catheterization procedures for patients with congenital and structural heart disease are becoming more complex. New imaging strategies involving integration of 3-dimensional images from rotational angiography, magnetic resonance imaging (MRI), computerized tomography (CT), and transesophageal echocardiography (TEE) are employed to facilitate these procedures. We discuss the current use of these new 3D imaging technologies and their advantages and challenges when used to guide complex diagnostic and interventional catheterization procedures in patients with congenital heart disease. PMID:25114757
Fractures in anisotropic media
NASA Astrophysics Data System (ADS)
Shao, Siyi
theory and experimental results in this report demonstrate that the presence of fractures in anisotropic material can be unambiguously interpreted if experimental measurements are made as a function of stress, which eliminates many fracture-generated discrete modes (e.g., interface waves, and leaky guided-modes). Orthogonal fracture networks that are often encountered in field exploration bring in additional challenges for seismic/acoustic data interpretation. An innovative wavefront imaging system with a bi-axial load frame was designed and implemented on orthogonally-fractured samples to determine the effect of fracture networks on elastic wave propagation. The effects of central wave guiding and extra time delays along a fracture intersection were observed in experiments and was analyzed. Interpreting data from media with intersecting fracture sets must account for fracture intersections and the non-uniformity of fracture properties caused by local tectonic conditions or other physical process such as non-uniform fluid distributions within a network and/or chemical alterations.
Fracture imaging with converted elastic waves
Nihei, K.T.; Nakagawa, S.; Myer, L.R.
2001-05-29
This paper examines the seismic signatures of discrete, finite-length fractures, and outlines an approach for elastic, prestack reverse-time imaging of discrete fractures. The results of this study highlight the importance of incorporating fracture-generated P-S converted waves into the imaging method, and presents an alternate imaging condition that can be used in elastic reverse-time imaging when a direct wave is recorded (e.g., for crosswell and VSP acquisition geometries).
The 3-dimensional cellular automata for HIV infection
NASA Astrophysics Data System (ADS)
Mo, Youbin; Ren, Bin; Yang, Wencao; Shuai, Jianwei
2014-04-01
The HIV infection dynamics is discussed in detail with a 3-dimensional cellular automata model in this paper. The model can reproduce the three-phase development, i.e., the acute period, the asymptotic period and the AIDS period, observed in the HIV-infected patients in a clinic. We show that the 3D HIV model performs a better robustness on the model parameters than the 2D cellular automata. Furthermore, we reveal that the occurrence of a perpetual source to successively generate infectious waves to spread to the whole system drives the model from the asymptotic state to the AIDS state.
NASA Technical Reports Server (NTRS)
Atluri, S. N.; Nakagaki, M.; Kathiresan, K.
1980-01-01
In this paper, efficient numerical methods for the analysis of crack-closure effects on fatigue-crack-growth-rates, in plane stress situations, and for the solution of stress-intensity factors for arbitrary shaped surface flaws in pressure vessels, are presented. For the former problem, an elastic-plastic finite element procedure valid for the case of finite deformation gradients is developed and crack growth is simulated by the translation of near-crack-tip elements with embedded plastic singularities. For the latter problem, an embedded-elastic-singularity hybrid finite element method, which leads to a direct evaluation of K-factors, is employed.
A numerical approach for pressure transient analysis of a vertical well with complex fractures
NASA Astrophysics Data System (ADS)
Wan, Yizhao; Liu, Yuewu; Liu, Wenchao; Han, Guofeng; Niu, Congcong
2016-05-01
A new well test model for a vertical fractured well is developed based on a discrete-fracture model in which the fractures are discretized as one dimensional (1-D) entities. The model overcomes the weakness of complex meshing, a large number of grids, and instability in conventional stripe-fracture models. Then, the discrete-fracture model is implemented using a hybrid element finite-element method. Triangular elements are used for matrix and line elements for the fractures. The finite element formulation is validated by comparing with the semi-analytical solution of a single vertical fractured well. The accuracy of the approach is shown through several examples with different fracture apertures, fracture conductivity, and fracture amount. Results from the discrete-fracture model agree reasonably well with the stripe-fracture model and the analytic solutions. The advantages of the discrete-fracture model are presented in mesh generation, computational improvement, and abilities to handle complex fractures like wedge-shaped fractures and fractures with branches. Analytical results show that the number of grids in the discrete-fracture model is 10 % less than stripe-fracture model, and computational efficiency increases by about 50 %. The more fractures there are, the more the computational efficiency increases.
3-Dimensional Imaging Modalities for Phenotyping Genetically Engineered Mice
Powell, K. A.; Wilson, D.
2013-01-01
A variety of 3-dimensional (3D) digital imaging modalities are available for whole-body assessment of genetically engineered mice: magnetic resonance microscopy (MRM), X-ray microcomputed tomography (microCT), optical projection tomography (OPT), episcopic and cryoimaging, and ultrasound biomicroscopy (UBM). Embryo and adult mouse phenotyping can be accomplished at microscopy or near microscopy spatial resolutions using these modalities. MRM and microCT are particularly well-suited for evaluating structural information at the organ level, whereas episcopic and OPT imaging provide structural and functional information from molecular fluorescence imaging at the cellular level. UBM can be used to monitor embryonic development longitudinally in utero. Specimens are not significantly altered during preparation, and structures can be viewed in their native orientations. Technologies for rapid automated data acquisition and high-throughput phenotyping have been developed and continually improve as this exciting field evolves. PMID:22146851
Protalign: a 3-dimensional protein alignment assessment tool.
Meads, D; Hansen, M D; Pang, A
1999-01-01
Protein fold recognition (sometimes called threading) is the prediction of a protein's 3-dimensional shape based on its similarity to a protein of known structure. Fold predictions are low resolution; that is, no effort is made to rotate the protein's component amino acid side chains into their correct spatial orientations. The goal is simply to recognize the protein family member that most closely resembles the target sequence of unknown structure and to create a sensible alignment of the target to the known structure (i.e., a structure-sequence alignment). To facilitate this type of structure prediction, we have designed a low resolution molecular graphics tool. ProtAlign introduces the ability to interact with and edit alignments directly in the 3-dimensional structure as well as in the usual 2-dimensional layout. It also contains several functions and features to help the user assess areas within the alignment. ProtAlign implements an open pipe architecture to allow other programs to access its molecular graphics capabilities. In addition, it is capable of "driving" other programs. Because amino acid side chain orientation is not relevant in fold recognition, we represent amino acid residues as abstract shapes or glyphs much like Lego (tm) blocks and we borrow techniques from comparative flow visualization using streamlines to provide clean depictions of the entire protein model. By creating a low resolution representation of protein structure, we are able to at least double the amount of information on the screen. At the same time, we create a view that is not as busy as the corresponding representations using traditional high resolution visualization methods which show detailed atomic structure. This eliminates distracting and possibly misleading visual clutter resulting from the mapping of protein alignment information onto a high resolution display of the known structure. This molecular graphics program is implemented in Open GL to facilitate porting to
Transport in fractured porous solids
NASA Astrophysics Data System (ADS)
Mattisson, Charlotte; Knackstedt, Mark A.; Senden, Tim J.
Laboratory measurements are made of the permeability and the resistivity of sintered porous media with disordered fractures over a wide range of matrix porosity. We discuss the preparation and characterisation of the samples. Approximating the topology of a rough fracture by a single discrete fracture can introduce large errors in the prediction of the permeability. We test the validity of empirical expressions relating permeability, resistivity and porosity for fractured samples. Resistivity correlations with porosity are independent of the presence of fractures. In contrast, permeability correlations show a strong dependence. Attempts to decouple the permeability of the fractured sample as a parallel sum of matrix and fracture permeability leads to large errors. The results indicate that transport in a medium with two distinct families of pathways cannot be described by a single-valued transport coefficient.
Hydraulically controlled discrete sampling from open boreholes
Harte, Philip T.
2013-01-01
Groundwater sampling from open boreholes in fractured-rock aquifers is particularly challenging because of mixing and dilution of fluid within the borehole from multiple fractures. This note presents an alternative to traditional sampling in open boreholes with packer assemblies. The alternative system called ZONFLO (zonal flow) is based on hydraulic control of borehole flow conditions. Fluid from discrete fractures zones are hydraulically isolated allowing for the collection of representative samples. In rough-faced open boreholes and formations with less competent rock, hydraulic containment may offer an attractive alternative to physical containment with packers. Preliminary test results indicate a discrete zone can be effectively hydraulically isolated from other zones within a borehole for the purpose of groundwater sampling using this new method.
Chromosome Conformation of Human Fibroblasts Grown in 3-Dimensional Spheroids
Chen, Haiming; Comment, Nicholas; Chen, Jie; Ronquist, Scott; Hero, Alfred; Ried, Thomas; Rajapakse, Indika
2015-01-01
In the study of interphase chromosome organization, genome-wide chromosome conformation capture (Hi-C) maps are often generated using 2-dimensional (2D) monolayer cultures. These 2D cells have morphological deviations from cells that exist in 3-dimensional (3D) tissues in vivo, and may not maintain the same chromosome conformation. We used Hi-C maps to test the extent of differences in chromosome conformation between human fibroblasts grown in 2D cultures and those grown in 3D spheroids. Significant differences in chromosome conformation were found between 2D cells and those grown in spheroids. Intra-chromosomal interactions were generally increased in spheroid cells, with a few exceptions, while inter-chromosomal interactions were generally decreased. Overall, chromosomes located closer to the nuclear periphery had increased intra-chromosomal contacts in spheroid cells, while those located more centrally had decreased interactions. This study highlights the necessity to conduct studies on the topography of the interphase nucleus under conditions that mimic an in vivo environment. PMID:25738643
Thermal crosstalk in 3-dimensional RRAM crossbar array
Sun, Pengxiao; Lu, Nianduan; Li, Ling; Li, Yingtao; Wang, Hong; Lv, Hangbing; Liu, Qi; Long, Shibing; Liu, Su; Liu, Ming
2015-01-01
High density 3-dimensional (3D) crossbar resistive random access memory (RRAM) is one of the major focus of the new age technologies. To compete with the ultra-high density NAND and NOR memories, understanding of reliability mechanisms and scaling potential of 3D RRAM crossbar array is needed. Thermal crosstalk is one of the most critical effects that should be considered in 3D crossbar array application. The Joule heat generated inside the RRAM device will determine the switching behavior itself, and for dense memory arrays, the temperature surrounding may lead to a consequent resistance degradation of neighboring devices. In this work, thermal crosstalk effect and scaling potential under thermal effect in 3D RRAM crossbar array are systematically investigated. It is revealed that the reset process is dominated by transient thermal effect in 3D RRAM array. More importantly, thermal crosstalk phenomena could deteriorate device retention performance and even lead to data storage state failure from LRS (low resistance state) to HRS (high resistance state) of the disturbed RRAM cell. In addition, the resistance state degradation will be more serious with continuously scaling down the feature size. Possible methods for alleviating thermal crosstalk effect while further advancing the scaling potential are also provided and verified by numerical simulation. PMID:26310537
Thermal crosstalk in 3-dimensional RRAM crossbar array.
Sun, Pengxiao; Lu, Nianduan; Li, Ling; Li, Yingtao; Wang, Hong; Lv, Hangbing; Liu, Qi; Long, Shibing; Liu, Su; Liu, Ming
2015-01-01
High density 3-dimensional (3D) crossbar resistive random access memory (RRAM) is one of the major focus of the new age technologies. To compete with the ultra-high density NAND and NOR memories, understanding of reliability mechanisms and scaling potential of 3D RRAM crossbar array is needed. Thermal crosstalk is one of the most critical effects that should be considered in 3D crossbar array application. The Joule heat generated inside the RRAM device will determine the switching behavior itself, and for dense memory arrays, the temperature surrounding may lead to a consequent resistance degradation of neighboring devices. In this work, thermal crosstalk effect and scaling potential under thermal effect in 3D RRAM crossbar array are systematically investigated. It is revealed that the reset process is dominated by transient thermal effect in 3D RRAM array. More importantly, thermal crosstalk phenomena could deteriorate device retention performance and even lead to data storage state failure from LRS (low resistance state) to HRS (high resistance state) of the disturbed RRAM cell. In addition, the resistance state degradation will be more serious with continuously scaling down the feature size. Possible methods for alleviating thermal crosstalk effect while further advancing the scaling potential are also provided and verified by numerical simulation. PMID:26310537
ERIC Educational Resources Information Center
Ghezzi, Patrick M.
2007-01-01
The advantages of emphasizing discrete trials "teaching" over discrete trials "training" are presented first, followed by a discussion of discrete trials as a method of teaching that emerged historically--and as a matter of necessity for difficult learners such as those with autism--from discrete trials as a method for laboratory research. The…
Video Based Sensor for Tracking 3-Dimensional Targets
NASA Technical Reports Server (NTRS)
Howard, R. T.; Book, Michael L.; Bryan, Thomas C.
2000-01-01
Video-Based Sensor for Tracking 3-Dimensional Targets The National Aeronautics and Space Administration's (NASAs) Marshall Space Flight Center (MSFC) has been developing and testing video-based sensors for automated spacecraft guidance for several years, and the next generation of video sensor will have tracking rates up to 100 Hz and will be able to track multiple reflectors and targets. The Video Guidance Sensor (VGS) developed over the past several years has performed well in testing and met the objective of being used as the terminal guidance sensor for an automated rendezvous and capture system. The first VGS was successfully tested in closed-loop 3-degree-of-freedom (3- DOF) tests in 1989 and then in 6-DOF open-loop tests in 1992 and closed-loop tests in 1993-4. Development and testing continued, and in 1995 approval was given to test the VGS in an experiment on the Space Shuttle. The VGS flew in 1997 and in 1998, performing well for both flights. During the development and testing before, during, and after the flight experiments, numerous areas for improvement were found. The VGS was developed with a sensor head and an electronics box, connected by cables. The VGS was used in conjunction with a target that had wavelength-filtered retro-reflectors in a specific pattern, The sensor head contained the laser diodes, video camera, and heaters and coolers. The electronics box contained a frame grabber, image processor, the electronics to control the components in the sensor head, the communications electronics, and the power supply. The system works by sequentially firing two different wavelengths of laser diodes at the target and processing the two images. Since the target only reflects one wavelength, it shows up well in one image and not at all in the other. Because the target's dimensions are known, the relative positions and attitudes of the target and the sensor can be computed from the spots reflected from the target. The system was designed to work from I
The 3-dimensional construction of the Rae craton, central Canada
NASA Astrophysics Data System (ADS)
Snyder, David B.; Craven, James A.; Pilkington, Mark; Hillier, Michael J.
2015-10-01
Reconstruction of the 3-dimensional tectonic assembly of early continents, first as Archean cratons and then Proterozoic shields, remains poorly understood. In this paper, all readily available geophysical and geochemical data are assembled in a 3-D model with the most accurate bedrock geology in order to understand better the geometry of major structures within the Rae craton of central Canada. Analysis of geophysical observations of gravity and seismic wave speed variations revealed several lithospheric-scale discontinuities in physical properties. Where these discontinuities project upward to correlate with mapped upper crustal geological structures, the discontinuities can be interpreted as shear zones. Radiometric dating of xenoliths provides estimates of rock types and ages at depth beneath sparse kimberlite occurrences. These ages can also be correlated to surface rocks. The 3.6-2.6 Ga Rae craton comprises at least three smaller continental terranes, which "cratonized" during a granitic bloom. Cratonization probably represents final differentiation of early crust into a relatively homogeneous, uniformly thin (35-42 km), tonalite-trondhjemite-granodiorite crust with pyroxenite layers near the Moho. The peak thermotectonic event at 1.86-1.7 Ga was associated with the Hudsonian orogeny that assembled several cratons and lesser continental blocks into the Canadian Shield using a number of southeast-dipping megathrusts. This orogeny metasomatized, mineralized, and recrystallized mantle and lower crustal rocks, apparently making them more conductive by introducing or concentrating sulfides or graphite. Little evidence exists of thin slabs similar to modern oceanic lithosphere in this Precambrian construction history whereas underthrusting and wedging of continental lithosphere is inferred from multiple dipping discontinuities.
A 3-Dimensional Anatomic Study of the Distal Biceps Tendon
Walton, Christine; Li, Zhi; Pennings, Amanda; Agur, Anne; Elmaraghy, Amr
2015-01-01
Background Complete rupture of the distal biceps tendon from its osseous attachment is most often treated with operative intervention. Knowledge of the overall tendon morphology as well as the orientation of the collagenous fibers throughout the musculotendinous junction are key to intraoperative decision making and surgical technique in both the acute and chronic setting. Unfortunately, there is little information available in the literature. Purpose To comprehensively describe the morphology of the distal biceps tendon. Study Design Descriptive laboratory study. Methods The distal biceps terminal musculature, musculotendinous junction, and tendon were digitized in 10 cadaveric specimens and data reconstructed using 3-dimensional modeling. Results The average length, width, and thickness of the external distal biceps tendon were found to be 63.0, 6.0, and 3.0 mm, respectively. A unique expansion of the tendon fibers within the distal muscle was characterized, creating a thick collagenous network along the central component between the long and short heads. Conclusion This study documents the morphologic parameters of the native distal biceps tendon. Reconstruction may be necessary, especially in chronic distal biceps tendon ruptures, if the remaining tendon morphology is significantly compromised compared with the native distal biceps tendon. Knowledge of normal anatomical distal biceps tendon parameters may also guide the selection of a substitute graft with similar morphological characteristics. Clinical Relevance A thorough description of distal biceps tendon morphology is important to guide intraoperative decision making between primary repair and reconstruction and to better select the most appropriate graft. The detailed description of the tendinous expansion into the muscle may provide insight into better graft-weaving and suture-grasping techniques to maximize proximal graft incorporation. PMID:26665092
Development and Validation of a 3-Dimensional CFB Furnace Model
NASA Astrophysics Data System (ADS)
Vepsäläinen, Arl; Myöhänen, Karl; Hyppäneni, Timo; Leino, Timo; Tourunen, Antti
At Foster Wheeler, a three-dimensional CFB furnace model is essential part of knowledge development of CFB furnace process regarding solid mixing, combustion, emission formation and heat transfer. Results of laboratory and pilot scale phenomenon research are utilized in development of sub-models. Analyses of field-test results in industrial-scale CFB boilers including furnace profile measurements are simultaneously carried out with development of 3-dimensional process modeling, which provides a chain of knowledge that is utilized as feedback for phenomenon research. Knowledge gathered by model validation studies and up-to-date parameter databases are utilized in performance prediction and design development of CFB boiler furnaces. This paper reports recent development steps related to modeling of combustion and formation of char and volatiles of various fuel types in CFB conditions. Also a new model for predicting the formation of nitrogen oxides is presented. Validation of mixing and combustion parameters for solids and gases are based on test balances at several large-scale CFB boilers combusting coal, peat and bio-fuels. Field-tests including lateral and vertical furnace profile measurements and characterization of solid materials provides a window for characterization of fuel specific mixing and combustion behavior in CFB furnace at different loads and operation conditions. Measured horizontal gas profiles are projection of balance between fuel mixing and reactions at lower part of furnace and are used together with both lateral temperature profiles at bed and upper parts of furnace for determination of solid mixing and combustion model parameters. Modeling of char and volatile based formation of NO profiles is followed by analysis of oxidizing and reducing regions formed due lower furnace design and mixing characteristics of fuel and combustion airs effecting to formation ofNO furnace profile by reduction and volatile-nitrogen reactions. This paper presents
NASA Astrophysics Data System (ADS)
Kettermann, Michael; von Hagke, Christoph; Virgo, Simon; Urai, Janos L.
2015-04-01
Brittle rocks are often affected by different generations of fractures that influence each other. We study pre-existing vertical joints followed by a faulting event. Understanding the effect of these interactions on fracture/fault geometries as well as the development of dilatancy and the formation of cavities as potential fluid pathways is crucial for reservoir quality prediction and production. Our approach combines scaled analogue and numerical modeling. Using cohesive hemihydrate powder allows us to create open fractures prior to faulting. The physical models are reproduced using the ESyS-Particle discrete element Modeling Software (DEM), and different parameters are investigated. Analogue models were carried out in a manually driven deformation box (30x28x20 cm) with a 60° dipping pre-defined basement fault and 4.5 cm of displacement. To produce open joints prior to faulting, sheets of paper were mounted in the box to a depth of 5 cm at a spacing of 2.5 cm. Powder was then sieved into the box, embedding the paper almost entirely (column height of 19 cm), and the paper was removed. We tested the influence of different angles between the strike of the basement fault and the joint set (0°, 4°, 8°, 12°, 16°, 20°, and 25°). During deformation we captured structural information by time-lapse photography that allows particle imaging velocimetry analyses (PIV) to detect localized deformation at every increment of displacement. Post-mortem photogrammetry preserves the final 3-dimensional structure of the fault zone. We observe that no faults or fractures occur parallel to basement-fault strike. Secondary fractures are mostly oriented normal to primary joints. At the final stage of the experiments we analyzed semi-quantitatively the number of connected joints, number of secondary fractures, degree of segmentation (i.e. number of joints accommodating strain), damage zone width, and the map-view area fraction of open gaps. Whereas the area fraction does not change
Unification of color postprocessing techniques for 3-dimensional computational mechanics
NASA Technical Reports Server (NTRS)
Bailey, Bruce Charles
1985-01-01
To facilitate the understanding of complex three-dimensional numerical models, advanced interactive color postprocessing techniques are introduced. These techniques are sufficiently flexible so that postprocessing difficulties arising from model size, geometric complexity, response variation, and analysis type can be adequately overcome. Finite element, finite difference, and boundary element models may be evaluated with the prototype postprocessor. Elements may be removed from parent models to be studied as independent subobjects. Discontinuous responses may be contoured including responses which become singular, and nonlinear color scales may be input by the user for the enhancement of the contouring operation. Hit testing can be performed to extract precise geometric, response, mesh, or material information from the database. In addition, stress intensity factors may be contoured along the crack front of a fracture model. Stepwise analyses can be studied, and the user can recontour responses repeatedly, as if he were paging through the response sets. As a system, these tools allow effective interpretation of complex analysis results.
Computer-assisted 3-dimensional anthropometry of the scaphoid.
Pichler, Wolfgang; Windisch, Gunther; Schaffler, Gottfried; Heidari, Nima; Dorr, Katrin; Grechenig, Wolfgang
2010-02-01
Scaphoid fracture fixation using a cannulated headless compression screw and the Matti-Russe procedure for the treatment of scaphoid nonunions are performed routinely. Surgeons performing these procedures need to be familiar with the anatomy of the scaphoid. A literature review reveals relatively few articles on this subject. The goal of this anatomical study was to measure the scaphoid using current technology and to discuss the findings with respect to the current, relevant literature.Computed tomography scans of 30 wrists were performed using a 64-slice SOMATOM Sensation CT system (resolution 0.6 mm) (Siemens Medical Solutions Inc, Malvern, Pennsylvania). Three-dimensional reconstructions from the raw data were generated by MIMICS software (Materialise, Leuven, Belgium). The scaphoid had a mean length of 26.0 mm (range, 22.3-30.7 mm), and men had a significantly longer (P<.001) scaphoid than women (27.861.6 mm vs 24.561.6 mm, respectively). The width and height were measured at 3 different levels for volume calculations, resulting in a mean volume of 3389.5 mm(3). Men had a significantly larger (P<.001) scaphoid volume than women (4057.86740.7 mm(3) vs 2846.56617.5 mm(3), respectively).We found considerable variation in the length and volume of the scaphoid in our cohort. We also demonstrated a clear correlation between scaphoid size and sex. Surgeons performing operative fixation of scaphoid fractures and corticocancellous bone grafting for nonunions need to be familiar with these anatomical variations. PMID:20192143
A 3-dimensional Analysis of the Cassiopeia A Supernova Remnant
NASA Astrophysics Data System (ADS)
Isensee, Karl
We present a multi-wavelength study of the nearby supernova remnant Cassiopeia A (Cas A). Easily resolvable supernova remnants such as Cas A provide a unique opportunity to test supernova explosion models. Additionally, we can observe key processes in the interstellar medium as the ejecta from the initial explosion encounter Cas A's powerful shocks. In order to accomplish these science goals, we used the Spitzer Space Telescope's Infrared Spectrograph to create a high resolution spectral map of select regions of Cas A, allowing us to make a Doppler reconstruction of its 3-dimensional structure structure. In the center of the remnant, we find relatively pristine ejecta that have not yet reached Cas A's reverse shock or interacted with the circumstellar environment. We observe O, Si, and S emission. These ejecta can form both sheet-like structures as well as filaments. Si and O, which come from different nucleosynthetic layers of the star, are observed to be coincident in some regions, and separated by >500 km s -1 in others. Observed ejecta traveling toward us are, on average, ˜800 km s -1 slower than the material traveling away from us. We compare our observations to recent supernova explosion models and find that no single model can simultaneously reproduce all the observed features. However, models of different supernova explosions can collectively produce the observed geometries and structures of the emission interior to Cas A's reverse shock. We use the results from the models to address the conditions during the supernova explosion, concentrating on asymmetries in the shock structure. We also predict that the back surface of Cassiopeia A will begin brightening in ∼30 years, and the front surface in ˜100 years. We then used similar observations from 3 regions on Cas A's reverse shock in order to create more 3-dimensional maps. In these regions, we observe supernova ejecta both immediately before and during the shock-ejecta interaction. We determine that the
Zhang, Boning; Herbold, Eric B.; Homel, Michael A.; Regueiro, Richard A.
2015-12-01
An adaptive particle fracture model in poly-ellipsoidal Discrete Element Method is developed. The poly-ellipsoidal particle will break into several sub-poly-ellipsoids by Hoek-Brown fracture criterion based on continuum stress and the maximum tensile stress in contacts. Also Weibull theory is introduced to consider the statistics and size effects on particle strength. Finally, high strain-rate split Hopkinson pressure bar experiment of silica sand is simulated using this newly developed model. Comparisons with experiments show that our particle fracture model can capture the mechanical behavior of this experiment very well, both in stress-strain response and particle size redistribution. The effects of density and packings o the samples are also studied in numerical examples.
3-Dimensional Marine CSEM Modeling by Employing TDFEM with Parallel Solvers
NASA Astrophysics Data System (ADS)
Wu, X.; Yang, T.
2013-12-01
In this paper, parallel fulfillment is developed for forward modeling of the 3-Dimensional controlled source electromagnetic (CSEM) by using time-domain finite element method (TDFEM). Recently, a greater attention rises on research of hydrocarbon (HC) reservoir detection mechanism in the seabed. Since China has vast ocean resources, seeking hydrocarbon reservoirs become significant in the national economy. However, traditional methods of seismic exploration shown a crucial obstacle to detect hydrocarbon reservoirs in the seabed with a complex structure, due to relatively high acquisition costs and high-risking exploration. In addition, the development of EM simulations typically requires both a deep knowledge of the computational electromagnetics (CEM) and a proper use of sophisticated techniques and tools from computer science. However, the complexity of large-scale EM simulations often requires large memory because of a large amount of data, or solution time to address problems concerning matrix solvers, function transforms, optimization, etc. The objective of this paper is to present parallelized implementation of the time-domain finite element method for analysis of three-dimensional (3D) marine controlled source electromagnetic problems. Firstly, we established a three-dimensional basic background model according to the seismic data, then electromagnetic simulation of marine CSEM was carried out by using time-domain finite element method, which works on a MPI (Message Passing Interface) platform with exact orientation to allow fast detecting of hydrocarbons targets in ocean environment. To speed up the calculation process, SuperLU of an MPI (Message Passing Interface) version called SuperLU_DIST is employed in this approach. Regarding the representation of three-dimension seabed terrain with sense of reality, the region is discretized into an unstructured mesh rather than a uniform one in order to reduce the number of unknowns. Moreover, high-order Whitney
Method and apparatus for imaging through 3-dimensional tracking of protons
NASA Technical Reports Server (NTRS)
Ryan, James M. (Inventor); Macri, John R. (Inventor); McConnell, Mark L. (Inventor)
2001-01-01
A method and apparatus for creating density images of an object through the 3-dimensional tracking of protons that have passed through the object are provided. More specifically, the 3-dimensional tracking of the protons is accomplished by gathering and analyzing images of the ionization tracks of the protons in a closely packed stack of scintillating fibers.
A compendium of fracture flow models, 1994
Diodato, D.M.
1994-11-01
The report is designed to be used as a decision-making aid for individuals who need to simulate fluid flow in fractured porous media. Fracture flow codes of varying capability in the public and private domain were identified in a survey of government, academia, and industry. The selection and use of an appropriate code requires conceptualization of the geology, physics, and chemistry (for transport) of the fracture flow problem to be solved. Conceptual models that have been invoked to describe fluid flow in fractured porous media include explicit discrete fracture, dual continuum (porosity and/or permeability), discrete fracture network, multiple interacting continua, multipermeability/multiporosity, and single equivalent continuum. The explicit discrete-fracture model is a ``near-field`` representation, the single equivalent continuum model is a ``far-field`` representation, and the dual-continuum model is intermediate to those end members. Of these, the dual-continuum model is the most widely employed. The concept of multiple interacting continua has been applied in a limited number of examples. Multipermeability/multiporosity provides a unified conceptual model. The ability to accurately describe fracture flow phenomena will continue to improve as a result of advances in fracture flow research and computing technology. This improvement will result in enhanced capability to protect the public environment, safety, and health.
Fracturing operations in a dry geothermal reservoir
Rowley, J.C.; Pettitt, R.A.; Hendron, R.H.; Sinclair, A.R.; Nicholson, R.W.
1983-01-01
Fracturing and completing deep wells in hot, nonporous crystalline basement rock challenges conventional equipment use, procedures, and techniques common in oil and gas and normal geothermal completions. Fracturing operations at the Fenton Hill, New Mexico, Hot Dry Rock (HDR) Geothermal Test Site initiated unique developments necessary to solve problems caused by an extremely harsh down-hole environment. Hydraulic fracturing experiments to connect the 2 wells have used openhole packers, hydraulic jet notching of the borehole wall, cemented-in isolation liners, and casing packers. Problems were encountered with hole drag, high fracture gradients, H/sub 2/S in vent back fluids, stress corrosion cracking of tubulars, and the complex nature of 3-dimensional fracture growth that requires large volumes of injected water. 20 references.
Principles of Discrete Time Mechanics
NASA Astrophysics Data System (ADS)
Jaroszkiewicz, George
2014-04-01
1. Introduction; 2. The physics of discreteness; 3. The road to calculus; 4. Temporal discretization; 5. Discrete time dynamics architecture; 6. Some models; 7. Classical cellular automata; 8. The action sum; 9. Worked examples; 10. Lee's approach to discrete time mechanics; 11. Elliptic billiards; 12. The construction of system functions; 13. The classical discrete time oscillator; 14. Type 2 temporal discretization; 15. Intermission; 16. Discrete time quantum mechanics; 17. The quantized discrete time oscillator; 18. Path integrals; 19. Quantum encoding; 20. Discrete time classical field equations; 21. The discrete time Schrodinger equation; 22. The discrete time Klein-Gordon equation; 23. The discrete time Dirac equation; 24. Discrete time Maxwell's equations; 25. The discrete time Skyrme model; 26. Discrete time quantum field theory; 27. Interacting discrete time scalar fields; 28. Space, time and gravitation; 29. Causality and observation; 30. Concluding remarks; Appendix A. Coherent states; Appendix B. The time-dependent oscillator; Appendix C. Quaternions; Appendix D. Quantum registers; References; Index.
... and held together with pins and wires or plates and screws. Fractures of the distal humerus (see ... doctor. These fractures usually require surgical repair with plates and/or screw, unless they are stable. SIGNS ...
Brolin, Tyler J; Throckmorton, Thomas
2015-11-01
Olecranon fractures are common upper extremity injuries, with all but nondisplaced fractures treated surgically. There has been a recent shift in the surgical management of these fractures from tension band wiring to locking plate fixation and intramedullary nailing; however, this comes with increased implant cost. Although most patients can expect good outcomes after these various techniques, there is little information to guide a surgeon's treatment plan. This article reviews the epidemiology, classification, treatment, and outcomes of olecranon fractures. PMID:26498547
Conventional 2.0 mm miniplates versus 3-D plates in mandibular fractures
Sadhwani, Bipin S.; Anchlia, Sonal
2013-01-01
Aim: To compare and evaluate the treatment outcome and postoperative complications in mandibular fractures using 2- and 3-dimensional miniplates. Materials and Methods: This study consisted of a sample of 28 patients (40 fracture sites) divided randomly but equally (single-blind control trial study) into two groups. Each group contains 14 patients (20 similar fracture sites in each group). Group 1 was treated with open reduction and internal fixation using 3-dimensional (3-D) miniplates. Group II was treated using 2-dimensional (2-D) 2-mm miniplates. Results: Out of 14 patients treated by conventional 2-mm miniplates, 2 patients developed occlusal discrepancy, another 2 had postoperative mobility at fracture site, and 1 developed plate failure and subsequent infection, which was treated by removal of the plate under antibiotic coverage. One patient treated by 3-dimensional plates had tooth damage. Statistical Analysis: Chi-square test. Conclusion: The results of this study suggested that the treatment of mandibular fractures (symphysis, parasymphysis, and angle) with 3-dimensional plates provided 3-dimensional stability and carried low morbidity and infection rates. The only probable limitations of 3-dimensional plates were excessive implant material due to the extra vertical bars incorporated for countering the torque forces. PMID:24205475
Flow Rate- and Fracture Property Dependence of Fracture-Matrix Ensemble Relative Permeability
NASA Astrophysics Data System (ADS)
Matthai, S. K.; Lang, P.; Bazrafkan, S.
2012-12-01
The grid-block scale ensemble relative permeability, kri of fractured porous rock with appreciable matrix permeability is of decisive interest to reservoir simulation and the prediction of production, injector-producer water breakthrough, and ultimate recovery. While the dynamic behaviour of naturally fractured reservoirs (NFR) already provides many clues about (pseudo) kri on the inter-well length scale, such data are difficult to interpret because, in the subsurface, the exact fracture geometry is unknown. Here we present numerical simulation results from discrete fracture and matrix (DFM) unstructured grid hybrid FEM-FVM simulation models, predicting the shape of fracture-matrix kri curves. In contrast to our earlier work, we also simulate capillary fracture matrix transfer (CFMT) and without relying the frequently made simplifying assumption that fracture saturation reflects fracture-matrix capillary pressure equilibrium. We also employ a novel discretization of saturation which permits jump discontinuities to develop across the fracture-matrix interface. This increased physical realism permits - for the first time - to test our earlier semi-analytical model of the flow rate dependence of relative permeability, ensuing from CFMT. The sensitivity analysis presented here constrains CMFT-related flow rate dependence of kri and illustrates how it manifests itself in two geometries of layer-restricted well-developed fracture patterns mapped in the field. We have also investigated the dependence of kri on fracture aperture as computed using discrete element analysis for plausible states of in situ stress. Our results indicate that fracture-matrix ensemble relative permeability can be matched with a new semi-analytic model taking into account the fracture-matrix flux ratio, the wetted fracture-matrix interface area as a function of saturation and the breakthrough saturation. However, we also detect a scale dependence of kri requiring a more elaborate treatment.
DeCoster, T. A.; Stevens, M. A.; Albright, J. P.
1994-01-01
Fractures occur in athletes and dramatically influence performance during competitive and recreational activities. Fractures occur in athletes as the result of repetitive stress, acute sports-related trauma and trauma outside of athletics. The literature provides general guidelines for treatment as well as a variety of statistics on the epidemiology of fractures by sport and level of participation. Athletes are healthy and motivated patients, and have high expectations regarding their level of function. These qualities make them good surgical candidates. Although closed treatment methods are appropriate for most sports fractures, an aggressive approach to more complicated fractures employing current techniques may optimize their subsequent performance. PMID:7719781
NASA Astrophysics Data System (ADS)
Fitzenz, D. D.; Miller, S. A.
2001-12-01
We present preliminary results from a 3-dimensional fault interaction model, with the fault system specified by the geometry and tectonics of the San Andreas Fault (SAF) system. We use the forward model for earthquake generation on interacting faults of Fitzenz and Miller [2001] that incorporates the analytical solutions of Okada [85,92], GPS-constrained tectonic loading, creep compaction and frictional dilatancy [Sleep and Blanpied, 1994, Sleep, 1995], and undrained poro-elasticity. The model fault system is centered at the Big Bend, and includes three large strike-slip faults (each discretized into multiple subfaults); 1) a 300km, right-lateral segment of the SAF to the North, 2) a 200km-long left-lateral segment of the Garlock fault to the East, and 3) a 100km-long right-lateral segment of the SAF to the South. In the initial configuration, three shallow-dipping faults are also included that correspond to the thrust belt sub-parallel to the SAF. Tectonic loading is decomposed into basal shear drag parallel to the plate boundary with a 35mm yr-1 plate velocity, and East-West compression approximated by a vertical dislocation surface applied at the far-field boundary resulting in fault-normal compression rates in the model space about 4mm yr-1. Our aim is to study the long-term seismicity characteristics, tectonic evolution, and fault interaction of this system. We find that overpressured faults through creep compaction are a necessary consequence of the tectonic loading, specifically where high normal stress acts on long straight fault segments. The optimal orientation of thrust faults is a function of the strike-slip behavior, and therefore results in a complex stress state in the elastic body. This stress state is then used to generate new fault surfaces, and preliminary results of dynamically generated faults will also be presented. Our long-term aim is to target measurable properties in or around fault zones, (e.g. pore pressures, hydrofractures, seismicity
Discrete Feature Approach for Heterogeneous Reservoir Production Enhancement
Dershowitz, William S.; Cladouhos, Trenton
2001-09-06
This progress report describes activities during the period January 1, 1999 to June 30, 1999. Work was carried out on 21 tasks. The major activity during the reporting period was the development and preliminary application of discrete fracture network (DFN) models for Stoney Point, South Oregon Basin, and North Oregon Basins project study sites. In addition, research was carried out on analysis algorithms for discrete future orientation.
Hydrodynamics of a vertical hydraulic fracture
Narasimhan, T.N.
1987-03-24
We have developed a numerical algorithm, HUBBERT, to simulate the hydrodynamics of a propagating vertical, rectangular fracture in an elastic porous medium. Based on the IFD method, this algorithm assumes fracture geometry to be prescribed. The breakdown and the creation of the incipient fracture is carried out according to the Hubbert-Willis theory. The propagation of the fracture is based on the criterion provided by Griffith, based on energy considerations. The deformation properties of the open fracture are based on simple elasticity solutions. The fracture is assumed to have an elliptical shape to a distance equal to the fracture height, beyond which the shape is assumed to be parallel plate. A consequence of Griffith's criterion is that the fracture must propagate in discrete steps. The parametric studies carried out suggest that for a clear understanding of the hydrodynamics of the hydraulic fracture many hitherto unrecognized parameters must be better understood. Among these parameters one might mention, efficiency, aperture of the newly formed fracture, stiffness of the newly formed fracture, relation between fracture aperture and permeability, and well bore compliance. The results of the studies indicate that the patterns of pressure transients and the magnitudes of fracture length appear to conform to field observations. In particular, the discrete nature of fracture propagation as well as the relevant time scales of interest inferred from the present work seem to be corroborated by seismic monitoring in the field. The results suggest that the estimation of least principal stress can be reliably made either with shut in data or with reinjection data provided that injection rates are very small.
NASA Astrophysics Data System (ADS)
Hooker, J. N.; Eichhubl, P.; Xu, G.; Ahn, H.; Fall, A.; Hargrove, P.; Laubach, S.; Ukar, E.
2011-12-01
The Cambrian Eriboll Formation quartzarenites contain abundant fractures with varying degrees of quartz cement infill. Fractures exist that are entirely sealed; are locally sealed by bridging cements but preserve pore space among bridges; are mostly open but lined with veneers of cement; or are devoid of cement. Fracture propagation in the Eriboll Formation is highly sensitive to the presence of pre-existing fractures. Fracture reactivation occurs in opening mode as individual fractures repeatedly open and are filled or bridged by syn-kinematic cements. As well, reactivation occurs in shear as opening of one fracture orientation coincides with shear displacement along pre-existing fractures of different orientations. The tendency for pre-existing fractures to slip varies in part by the extent of cement infill, yet we observe shear and opening-mode reactivation even among sealed fractures. Paleotemperature analysis of fluid inclusions within fracture cements suggests some fractures now in outcrop formed deep in the subsurface. Fractures within the Eriboll Formation may therefore affect later fracture propagation throughout geologic time. With progressive strain, fault zones develop within fracture networks by a sequence of opening-mode fracture formation, fracture reactivation and linkage, fragmentation, cataclasis, and the formation of slip surfaces. Cataclasite within fault zones is commonly more thoroughly cemented than fractures in the damage zone or outside the fault zone. This variance of cement abundance is likely the result of (1) continued exposure of freshly broken quartz surfaces within cataclasite, promoting quartz precipitation, and (2) possibly more interconnected pathways for mass transfer within the fault zone. Enhanced cementation of cataclasite results in strengthening or diagenetic strain hardening of the evolving fault zone. Further slip is accommodated by shear localization along discrete slip surfaces. With further linkage of fault segments
Inter-trochanteric fracture repair - discharge; Subtrochanteric fracture repair - discharge; Femoral neck fracture repair - discharge; Trochanteric fracture repair - discharge; Hip pinning surgery - discharge
Morris, J; Johnson, S
2007-12-03
The Distinct Element Method (also frequently referred to as the Discrete Element Method) (DEM) is a Lagrangian numerical technique where the computational domain consists of discrete solid elements which interact via compliant contacts. This can be contrasted with Finite Element Methods where the computational domain is assumed to represent a continuum (although many modern implementations of the FEM can accommodate some Distinct Element capabilities). Often the terms Discrete Element Method and Distinct Element Method are used interchangeably in the literature, although Cundall and Hart (1992) suggested that Discrete Element Methods should be a more inclusive term covering Distinct Element Methods, Displacement Discontinuity Analysis and Modal Methods. In this work, DEM specifically refers to the Distinct Element Method, where the discrete elements interact via compliant contacts, in contrast with Displacement Discontinuity Analysis where the contacts are rigid and all compliance is taken up by the adjacent intact material.
Synchronous Discrete Harmonic Oscillator
Antippa, Adel F.; Dubois, Daniel M.
2008-10-17
We introduce the synchronous discrete harmonic oscillator, and present an analytical, numerical and graphical study of its characteristics. The oscillator is synchronous when the time T for one revolution covering an angle of 2{pi} in phase space, is an integral multiple N of the discrete time step {delta}t. It is fully synchronous when N is even. It is pseudo-synchronous when T/{delta}t is rational. In the energy conserving hyperincursive representation, the phase space trajectories are perfectly stable at all time scales, and in both synchronous and pseudo-synchronous modes they cycle through a finite number of phase space points. Consequently, both the synchronous and the pseudo-synchronous hyperincursive modes of time-discretization provide a physically realistic and mathematically coherent, procedure for dynamic, background independent, discretization of spacetime. The procedure is applicable to any stable periodic dynamical system, and provokes an intrinsic correlation between space and time, whereby space-discretization is a direct consequence of background-independent time-discretization. Hence, synchronous discretization moves the formalism of classical mechanics towards that of special relativity. The frequency of the hyperincursive discrete harmonic oscillator is ''blue shifted'' relative to its continuum counterpart. The frequency shift has the precise value needed to make the speed of the system point in phase space independent of the discretizing time interval {delta}t. That is the speed of the system point is the same on the polygonal (in the discrete case) and the circular (in the continuum case) phase space trajectories.
Fracture-permeability behavior of shale
Carey, J. William; Lei, Zhou; Rougier, Esteban; Mori, Hiroko; Viswanathan, Hari
2015-05-08
The fracture-permeability behavior of Utica shale, an important play for shale gas and oil, was investigated using a triaxial coreflood device and X-ray tomography in combination with finite-discrete element modeling (FDEM). Fractures generated in both compression and in a direct-shear configuration allowed permeability to be measured across the faces of cylindrical core. Shale with bedding planes perpendicular to direct-shear loading developed complex fracture networks and peak permeability of 30 mD that fell to 5 mD under hydrostatic conditions. Shale with bedding planes parallel to shear loading developed simple fractures with peak permeability as high as 900 mD. In addition to the large anisotropy in fracture permeability, the amount of deformation required to initiate fractures was greater for perpendicular layering (about 1% versus 0.4%), and in both cases activation of existing fractures are more likely sources of permeability in shale gas plays or damaged caprock in CO₂ sequestration because of the significant deformation required to form new fracture networks. FDEM numerical simulations were able to replicate the main features of the fracturing processes while showing the importance of fluid penetration into fractures as well as layering in determining fracture patterns.
Fracture-permeability behavior of shale
Carey, J. William; Lei, Zhou; Rougier, Esteban; Mori, Hiroko; Viswanathan, Hari
2015-05-08
The fracture-permeability behavior of Utica shale, an important play for shale gas and oil, was investigated using a triaxial coreflood device and X-ray tomography in combination with finite-discrete element modeling (FDEM). Fractures generated in both compression and in a direct-shear configuration allowed permeability to be measured across the faces of cylindrical core. Shale with bedding planes perpendicular to direct-shear loading developed complex fracture networks and peak permeability of 30 mD that fell to 5 mD under hydrostatic conditions. Shale with bedding planes parallel to shear loading developed simple fractures with peak permeability as high as 900 mD. In addition tomore » the large anisotropy in fracture permeability, the amount of deformation required to initiate fractures was greater for perpendicular layering (about 1% versus 0.4%), and in both cases activation of existing fractures are more likely sources of permeability in shale gas plays or damaged caprock in CO₂ sequestration because of the significant deformation required to form new fracture networks. FDEM numerical simulations were able to replicate the main features of the fracturing processes while showing the importance of fluid penetration into fractures as well as layering in determining fracture patterns.« less
Where Does Water Go During Hydraulic Fracturing?
O'Malley, D; Karra, S; Currier, R P; Makedonska, N; Hyman, J D; Viswanathan, H S
2016-07-01
During hydraulic fracturing millions of gallons of water are typically injected at high pressure into deep shale formations. This water can be housed in fractures, within the shale matrix, and can potentially migrate beyond the shale formation via fractures and/or faults raising environmental concerns. We describe a generic framework for producing estimates of the volume available in fractures and undamaged shale matrix where water injected into a representative shale site could reside during hydraulic fracturing, and apply it to a representative site that incorporates available field data. The amount of water that can be stored in the fractures is estimated by calculating the volume of all the fractures associated with a discrete fracture network (DFN) based on real data and using probability theory to estimate the volume of smaller fractures that are below the lower cutoff for the fracture radius in the DFN. The amount of water stored in the matrix is estimated utilizing two distinct methods-one using a two-phase model at the pore-scale and the other using a single-phase model at the continuum scale. Based on these calculations, it appears that most of the water resides in the matrix with a lesser amount in the fractures. PMID:26469857
Carlsten, B.E.; Haynes, W.B.
1996-08-01
The authors theoretically and numerically investigate the operation and behavior of the discrete monotron oscillator, a novel high-power microwave source. The discrete monotron differs from conventional monotrons and transit time oscillators by shielding the electron beam from the monotron cavity`s RF fields except at two distinct locations. This makes the discrete monotron act more like a klystron than a distributed traveling wave device. As a result, the oscillator has higher efficiency and can operate with higher beam powers than other single cavity oscillators and has more stable operation without requiring a seed input signal than mildly relativistic, intense-beam klystron oscillators.
Statistical physics models for nacre fracture simulation
NASA Astrophysics Data System (ADS)
Nukala, Phani Kumar V. V.; Šimunović, Srđan
2005-10-01
Natural biological materials such as nacre (or mother-of-pearl), exhibit phenomenal fracture strength and toughness properties despite the brittle nature of their constituents. For example, nacre’s work of fracture is three orders of magnitude greater than that of a single crystal of its constituent mineral. This study investigates the fracture properties of nacre using a simple discrete lattice model based on continuous damage random thresholds fuse network. The discrete lattice topology of the proposed model is based on nacre’s unique brick and mortar microarchitecture, and the mechanical behavior of each of the bonds in the discrete lattice model is governed by the characteristic modular damage evolution of the organic matrix that includes the mineral bridges between the aragonite platelets. The analysis indicates that the excellent fracture properties of nacre are a result of their unique microarchitecture, repeated unfolding of protein molecules (modular damage evolution) in the organic polymer, and the presence of fiber bundle of mineral bridges between the aragonite platelets. The numerical results obtained using this simple discrete lattice model are in excellent agreement with the previously obtained experimental results, such as nacre’s stiffness, tensile strength, and work of fracture.
Statistical physics models for nacre fracture simulation.
Nukala, Phani Kumar V V; Simunović, Srdan
2005-10-01
Natural biological materials such as nacre (or mother-of-pearl), exhibit phenomenal fracture strength and toughness properties despite the brittle nature of their constituents. For example, nacre's work of fracture is three orders of magnitude greater than that of a single crystal of its constituent mineral. This study investigates the fracture properties of nacre using a simple discrete lattice model based on continuous damage random thresholds fuse network. The discrete lattice topology of the proposed model is based on nacre's unique brick and mortar microarchitecture, and the mechanical behavior of each of the bonds in the discrete lattice model is governed by the characteristic modular damage evolution of the organic matrix that includes the mineral bridges between the aragonite platelets. The analysis indicates that the excellent fracture properties of nacre are a result of their unique microarchitecture, repeated unfolding of protein molecules (modular damage evolution) in the organic polymer, and the presence of fiber bundle of mineral bridges between the aragonite platelets. The numerical results obtained using this simple discrete lattice model are in excellent agreement with the previously obtained experimental results, such as nacre's stiffness, tensile strength, and work of fracture. PMID:16383432
ERIC Educational Resources Information Center
Peters, James V.
2004-01-01
Using the methods of finite difference equations the discrete analogue of the parabolic and catenary cable are analysed. The fibonacci numbers and the golden ratio arise in the treatment of the catenary.
Discretizations of axisymmetric systems
NASA Astrophysics Data System (ADS)
Frauendiener, Jörg
2002-11-01
In this paper we discuss stability properties of various discretizations for axisymmetric systems including the so-called cartoon method which was proposed by Alcubierre et al. for the simulation of such systems on Cartesian grids. We show that within the context of the method of lines such discretizations tend to be unstable unless one takes care in the way individual singular terms are treated. Examples are given for the linear axisymmetric wave equation in flat space.
Construction of 3-Dimensional Printed Ultrasound Phantoms With Wall-less Vessels.
Nikitichev, Daniil I; Barburas, Anamaria; McPherson, Kirstie; Mari, Jean-Martial; West, Simeon J; Desjardins, Adrien E
2016-06-01
Ultrasound phantoms are invaluable as training tools for vascular access procedures. We developed ultrasound phantoms with wall-less vessels using 3-dimensional printed chambers. Agar was used as a soft tissue-mimicking material, and the wall-less vessels were created with rods that were retracted after the agar was set. The chambers had integrated luer connectors to allow for fluid injections with clinical syringes. Several variations on this design are presented, which include branched and stenotic vessels. The results show that 3-dimensional printing can be well suited to the construction of wall-less ultrasound phantoms, with designs that can be readily customized and shared electronically. PMID:27162278
3-Dimensional Terraced NAND (3D TNAND) Flash Memory-Stacked Version of Folded NAND Array
NASA Astrophysics Data System (ADS)
Kim, Yoon; Cho, Seongjae; Lee, Gil Sung; Park, Il Han; Lee, Jong Duk; Shin, Hyungcheol; Park, Byung-Gook
We propose a 3-dimensional terraced NAND flash memory. It has a vertical channel so it is possible to make a long enough channel in 1F2 size. And it has 3-dimensional structure whose channel is connected vertically along with two stairs. So we can obtain high density as in the stacked array structure, without silicon stacking process. We can make NAND flash memory with 3F2 cell size. Using SILVACO ATLAS simulation, we study terraced NAND flash memory characteristics such as program, erase, and read. Also, its fabrication method is proposed.
Magnetic topologies of coronal mass ejection events: Effects of 3-dimensional reconnection
Gosling, J.T.
1995-09-01
New magnetic loops formed in the corona following coronal mass ejection, CME, liftoffs provide strong evidence that magnetic reconnection commonly occurs within the magnetic ``legs`` of the departing CMEs. Such reconnection is inherently 3-dimensional and naturally produces CMEs having magnetic flux rope topologies. Sustained reconnection behind CMEs can produce a mixture of open and disconnected field lines threading the CMES. In contrast to the results of 2-dimensional reconnection. the disconnected field lines are attached to the outer heliosphere at both ends. A variety of solar and solar wind observations are consistent with the concept of sustained 3-dimensional reconnection within the magnetic legs of CMEs close to the Sun.
Goggin, L M; Descovich, M; McGuinness, C; Shiao, S; Pouliot, J; Park, C
2016-06-01
Accelerated partial breast irradiation is an attractive alternative to conventional whole breast radiotherapy for selected patients. Recently, CyberKnife has emerged as a possible alternative to conventional techniques for accelerated partial breast irradiation. In this retrospective study, we present a dosimetric comparison between 3-dimensional conformal radiotherapy plans and CyberKnife plans using circular (Iris) and multi-leaf collimators. Nine patients who had undergone breast-conserving surgery followed by whole breast radiation were included in this retrospective study. The CyberKnife planning target volume (PTV) was defined as the lumpectomy cavity + 10 mm + 2 mm with prescription dose of 30 Gy in 5 fractions. Two sets of 3-dimensional conformal radiotherapy plans were created, one used the same definitions as described for CyberKnife and the second used the RTOG-0413 definition of the PTV: lumpectomy cavity + 15 mm + 10 mm with prescription dose of 38.5 Gy in 10 fractions. Using both PTV definitions allowed us to compare the dose delivery capabilities of each technology and to evaluate the advantage of CyberKnife tracking. For the dosimetric comparison using the same PTV margins, CyberKnife and 3-dimensional plans resulted in similar tumor coverage and dose to critical structures, with the exception of the lung V5%, which was significantly smaller for 3-dimensional conformal radiotherapy, 6.2% when compared to 39.4% for CyberKnife-Iris and 17.9% for CyberKnife-multi-leaf collimator. When the inability of 3-dimensional conformal radiotherapy to track motion is considered, the result increased to 25.6%. Both CyberKnife-Iris and CyberKnife-multi-leaf collimator plans demonstrated significantly lower average ipsilateral breast V50% (25.5% and 24.2%, respectively) than 3-dimensional conformal radiotherapy (56.2%). The CyberKnife plans were more conformal but less homogeneous than the 3-dimensional conformal radiotherapy plans. Approximately 50% shorter
Dual permeability modeling of flow in a fractured geothermal reservoir
Miller, J.D.; Allman, D.W.
1986-01-01
A three dimensional fracture system synthesis and flow simulation has been developed to correlate drawdown characteristics measured in a geothermal well and to provide the basis for an analysis of tracer tests. A new dual permeability approach was developed which incorporates simulations at two levels to better represent a discrete fracture system within computer limitations. The first incorporates a discrete simulation of the largest fractures in the system plus distributed or representative element stimulation of the smaller fractures. The second determines the representative element properties by discrete simulation of the smaller fractures. The fracture system was synthesized from acoustic televiewer data on the orientation and separation of three distinct fracture sets, together with additional data from the literature. Lognormal and exponential distributions of fracture spacing and radius were studied with the exponential distribution providing more reasonable results. Hydraulic apertures were estimated as a function of distance from the model boundary to a constant head boundary. Mean values of 6.7, 101 and 46 ..mu..m were chosen as the most representative values for the three fracture sets. Recommendations are given for the additional fracture characterization needed to reduce the uncertainties in the model. 20 refs., 6 figs.
Simulating conservative tracers in fractured till under realistic timescales.
Helmke, M F; Simpkins, W W; Horton, R
2005-01-01
Discrete-fracture and dual-porosity models are infrequently used to simulate solute transport through fractured unconsolidated deposits, despite their more common application in fractured rock where distinct flow regimes are hypothesized. In this study, we apply four fracture transport models--the mobile-immobile model (MIM), parallel-plate discrete-fracture model (PDFM), and stochastic and deterministic discrete-fracture models (DFMs)--to demonstrate their utility for simulating solute transport through fractured till. Model results were compared to breakthrough curves (BTCs) for the conservative tracers potassium bromide (KBr), pentafluorobenzoic acid (PFBA), and 1,4-piperazinediethanesulfonic acid (PIPES) in a large-diameter column of fractured till. Input parameters were determined from independent field and laboratory methods. Predictions of Br BTCs were not significantly different among models; however, the stochastic and deterministic DFMs were more accurate than the MIM or PDFM when predicting PFBA and PIPES BTCs. DFMs may be more applicable than the MIM for tracers with small effective diffusion coefficients (De) or for short timescales due to differences in how these models simulate diffusion or incorporate heterogeneities by their fracture networks. At large scales of investigation, the more computationally efficient MIM and PDFM may be more practical to implement than the three-dimensional DFMs, or a combination of model approaches could be employed. Regardless of the modeling approach used, fractures should be incorporated routinely into solute transport models in glaciated terrain. PMID:16324009
Dual Permeability Modeling of Flow in a Fractured Geothermal Reservoir
Miller, John D.; Allman, David W.
1986-01-21
A three dimensional fracture system synthesis and flow simulation has been developed to correlate drawdown characteristics measured in a geothermal well and to provide the basis for an analysis of tracer tests. A new dual permeability approach was developed which incorporates simulations at two levels to better represent a discrete fracture system within computer limitations. The first incorporates a discrete simulation of the largest fractures in the system plus distributed or representative element simulation of the smaller fractures. the second determines the representative element properties by discrete simulation of the smaller fractures. The fracture system was synthesized from acoustic televiewer data on the orientation and separation of three distinct fracture sets, together with additional data from the literature. Lognormal and exponential distributions of fracture spacing and radius were studied with the exponential distribution providing more reasonable results. Hydraulic apertures were estimated as a function of distance from the model boundary to a constant head boundary. Mean values of 6.7, 101 and 46 {micro}m were chosen as the most representative values for the three fracture sets. Recommendations are given for the additional fracture characterization needed to reduce the uncertainties in the model.
Fracture Characterization in the Astor Pass Geothermal Field, Nevada
NASA Astrophysics Data System (ADS)
Walsh, D. C.; Reeves, D. M.; Pohll, G.; Lyles, B. F.; Cooper, C. A.
2011-12-01
The Astor Pass geothermal field, near Pyramid Lake, NV, is under study as a site of potential geothermal energy production. Three wells have been completed in the graben of this typical Basin and Range geologic setting. Lithologies include a layer of unconsolidated sediment (basin fill) underlain by various tertiary volcanic units and granodiorite and metavolcanic basement rock. Characterization of fractures within the relatively impermeable rock matrix is being conducted for the three wells. Statistical analysis of fracture orientation, densities, and spacing obtained from borehole imaging logs is used to determine stress orientation and to generate a statistically equivalent Discrete Fracture Network (DFN) model. Fractures at depth are compared to fracture data collected in nearby outcrops of the same lithologic stratigraphy. Fracture geometry and density is correlated to mechanically discrete layers within the stratigraphy to test whether variations in fracturing can be attributed to variations in Young's modulus. Correlation of fracture geometry and densities with spinner flowmeter logs and distributed temperature sensor records are made in an effort to identify potential flowing fracture zones intersecting the borehole. Mean fracture aperture is obtained from open fracture counts and reservoir-scale transmissivity values (computed from a 30 day pump test) in the absence of readily available aperture data. The goal of this thorough fracture characterization is to create a physically relevant model which may be coupled with a multipurpose fluid flow and thermal simulator for investigation of geothermal reservoir behavior, particularly at the borehole scale.
3-dimensional root phenotyping with a novel imaging and software platform
Technology Transfer Automated Retrieval System (TEKTRAN)
A novel imaging and software platform was developed for the high-throughput phenotyping of 3-dimensional root traits during seedling development. To demonstrate the platform’s capacity, plants of two rice (Oryza sativa) genotypes, Azucena and IR64, were grown in a transparent gellan gum system and ...
3-DIMENSIONAL MEASURED AND SIMULATED FLOW FOR SCOUR NEAR SPUR DIKES
Technology Transfer Automated Retrieval System (TEKTRAN)
To improve understanding of the flow and scour processes associated with spur dikes more fully, 3-dimensional flow velocities were measured using an acoustic Doppler velocimeter at a closely spaced grid over a fixed flat bed with a submerged spur dike. Some 2592 three-dimensional velocities around a...
3-dimensional orthodontics visualization system with dental study models and orthopantomograms
NASA Astrophysics Data System (ADS)
Zhang, Hua; Ong, S. H.; Foong, K. W. C.; Dhar, T.
2005-04-01
The aim of this study is to develop a system that provides 3-dimensional visualization of orthodontic treatments. Dental plaster models and corresponding orthopantomogram (dental panoramic tomogram) are first digitized and fed into the system. A semi-auto segmentation technique is applied to the plaster models to detect the dental arches, tooth interstices and gum margins, which are used to extract individual crown models. 3-dimensional representation of roots, generated by deforming generic tooth models with orthopantomogram using radial basis functions, is attached to corresponding crowns to enable visualization of complete teeth. An optional algorithm to close the gaps between deformed roots and actual crowns by using multi-quadratic radial basis functions is also presented, which is capable of generating smooth mesh representation of complete 3-dimensional teeth. User interface is carefully designed to achieve a flexible system with as much user friendliness as possible. Manual calibration and correction is possible throughout the data processing steps to compensate occasional misbehaviors of automatic procedures. By allowing the users to move and re-arrange individual teeth (with their roots) on a full dentition, this orthodontic visualization system provides an easy and accurate way of simulation and planning of orthodontic treatment. Its capability of presenting 3-dimensional root information with only study models and orthopantomogram is especially useful for patients who do not undergo CT scanning, which is not a routine procedure in most orthodontic cases.
Characterization of Students' Reasoning and Proof Abilities in 3-Dimensional Geometry
ERIC Educational Resources Information Center
Gutierrez, Angel; Pegg, John; Lawrie, Christine
2004-01-01
In this paper we report on a research aimed to identify and characterize secondary school students' reasoning and proof abilities when working with 3-dimensional geometric solids. We analyze students' answers to two problems asking them to prove certain properties of prisms. As results of this analysis, we get, on the one side, a characterization…
ERIC Educational Resources Information Center
Ince, Elif; Kirbaslar, Fatma Gulay; Yolcu, Ergun; Aslan, Ayse Esra; Kayacan, Zeynep Cigdem; Alkan Olsson, Johanna; Akbasli, Ayse Ceylan; Aytekin, Mesut; Bauer, Thomas; Charalambis, Dimitris; Gunes, Zeliha Ozsoy; Kandemir, Ceyhan; Sari, Umit; Turkoglu, Suleyman; Yaman, Yavuz; Yolcu, Ozgu
2014-01-01
The purpose of this study is to develop a 3-dimensional interactive multi-user and multi-admin IUVIRLAB featuring active learning methods and techniques for university students and to introduce the Virtual Laboratory of Istanbul University and to show effects of IUVIRLAB on students' attitudes on communication skills and IUVIRLAB. Although…
A discrete fractional random transform
NASA Astrophysics Data System (ADS)
Liu, Zhengjun; Zhao, Haifa; Liu, Shutian
2005-11-01
We propose a discrete fractional random transform based on a generalization of the discrete fractional Fourier transform with an intrinsic randomness. Such discrete fractional random transform inheres excellent mathematical properties of the fractional Fourier transform along with some fantastic features of its own. As a primary application, the discrete fractional random transform has been used for image encryption and decryption.
... to hold the fracture in the correct position. • Fiberglass casting is lighter and stronger and the exterior ... with your physician if this occurs. • When a fiberglass cast is used in conjunction with a GORE- ...
Discrete Newtonian cosmology: perturbations
NASA Astrophysics Data System (ADS)
Ellis, George F. R.; Gibbons, Gary W.
2015-03-01
In a previous paper (Gibbons and Ellis 2014 Discrete Newtonian cosmology Class. Quantum Grav. 31 025003), we showed how a finite system of discrete particles interacting with each other via Newtonian gravitational attraction would lead to precisely the same dynamical equations for homothetic motion as in the case of the pressure-free Friedmann-Lemaître-Robertson-Walker cosmological models of general relativity theory, provided the distribution of particles obeys the central configuration equation. In this paper we show that one can obtain perturbed such Newtonian solutions that give the same linearized structure growth equations as in the general relativity case. We also obtain the Dmitriev-Zel’dovich equations for subsystems in this discrete gravitational model, and show how it leads to the conclusion that voids have an apparent negative mass.
Flow upscaling in propped fracture
NASA Astrophysics Data System (ADS)
Jasinski, Lukasz; Dabrowski, Marcin
2016-04-01
Proppants in combination with hydraulic fracturing are widely used to maintain the production of oil or gas from low permeability formations (i.e. shale rocks). There are also examples of proppants use in geothermal reservoirs. Flow patterns in propped fracture control transport processes and give information about fracture/matrix exchange surface. Our main motivation is to understand flow behavior in such structures using direct numerical simulations and to find a good upscaling technique to be able to investigate models on reservoir scale. We study fracture made of two parallel plane walls, where void space between them is filled with partial monolayer of proppant. As the fracture is affected by closing pressure, the proppant grains are squeezed between two opposite fracture walls which can change the grain shapes or embed the grains into impermeable rock matrix. To take this effect into account and simplify the geometry, the grains are approximated as cylinders. Imposed macroscopic pressure gradient invokes flow in such medium. As the flow is considered in the low Reynolds number regime, a stationary velocity flow field is obtained by solving the Stokes equations in 3D by means of finite element method. Void space between the grains is accurately discretized by using tetrahedral mesh. To reduce computational effort, the Stokes equation is reduced over the fracture aperture to 2D Stokes-Brinkman equation, which is further numerically solved and compared against numerical solution in 3D. Systematic flow calculations using 2D Stokes-Brinkman equation are performed for periodic domain and no slip boundary condition on the grain surface. Results are discussed in terms of effective properties as a function of geometrical parameters of the medium, such as proppant packing fraction and proppant grain diameter to fracture aperture ratio.
Discrete breathers in crystals
NASA Astrophysics Data System (ADS)
Dmitriev, S. V.; Korznikova, E. A.; Baimova, Yu A.; Velarde, M. G.
2016-05-01
It is well known that periodic discrete defect-containing systems, in addition to traveling waves, support vibrational defect-localized modes. It turned out that if a periodic discrete system is nonlinear, it can support spatially localized vibrational modes as exact solutions even in the absence of defects. Since the nodes of the system are all on equal footing, it is only through the special choice of initial conditions that a group of nodes can be found on which such a mode, called a discrete breather (DB), will be excited. The DB frequency must be outside the frequency range of the small-amplitude traveling waves. Not resonating with and expending no energy on the excitation of traveling waves, a DB can theoretically conserve its vibrational energy forever provided no thermal vibrations or other perturbations are present. Crystals are nonlinear discrete systems, and the discovery in them of DBs was only a matter of time. It is well known that periodic discrete defect-containing systems support both traveling waves and vibrational defect-localized modes. It turns out that if a periodic discrete system is nonlinear, it can support spatially localized vibrational modes as exact solutions even in the absence of defects. Because the nodes of the system are all on equal footing, only a special choice of the initial conditions allows selecting a group of nodes on which such a mode, called a discrete breather (DB), can be excited. The DB frequency must be outside the frequency range of small-amplitude traveling waves. Not resonating with and expending no energy on the excitation of traveling waves, a DB can theoretically preserve its vibrational energy forever if no thermal vibrations or other perturbations are present. Crystals are nonlinear discrete systems, and the discovery of DBs in them was only a matter of time. Experimental studies of DBs encounter major technical difficulties, leaving atomistic computer simulations as the primary investigation tool. Despite
NASA Astrophysics Data System (ADS)
Arzano, Michele; Kowalski-Glikman, Jerzy
2016-09-01
We construct discrete symmetry transformations for deformed relativistic kinematics based on group valued momenta. We focus on the specific example of κ-deformations of the Poincaré algebra with associated momenta living on (a sub-manifold of) de Sitter space. Our approach relies on the description of quantum states constructed from deformed kinematics and the observable charges associated with them. The results we present provide the first step towards the analysis of experimental bounds on the deformation parameter κ to be derived via precision measurements of discrete symmetries and CPT.
ERIC Educational Resources Information Center
Sharp, Karen Tobey
This paper cites information received from a number of sources, e.g., mathematics teachers in two-year colleges, publishers, and convention speakers, about the nature of discrete mathematics and about what topics a course in this subject should contain. Note is taken of the book edited by Ralston and Young which discusses the future of college…
Method for hydraulic fracturing cased wellbores
Schmidt, J.H.
1991-12-24
This patent describes a method of hydraulically fracturing a cased wellbore in an earth formation. It comprises determining the angle with respect to the wellbore axis and a reference point on the circumference of the wellbore which will provide for initiation of a hydraulic fracture in the formation which will turn with the largest radius of curvature into a fracture plane normal to the minimum in situ stress in the formation; perforating the wellbore casing at the angle with respect to the reference point; initiating a hydraulic fracture in the formation by pumping a liquid through the perforation and into the formation to force the initiation of a fracture in the formation at a point which develops the highest tensile stress in the formation in relation to increasing the hydraulic pressure in the wellbore; extending the fracture by pumping a relatively proppant-free quantities of proppant per unit volume of pumped fluid and in successive discrete stages of increasing proppant density to provide a propped portion of increasing proppant density to provide a propped portion of the fracture in the near wellbore region of the fracture which will prevent reclosing of the fracture in the near wellbore region.
Well test analysis in fractured media
Karasaki, K.
1986-04-01
In this study the behavior of fracture systems under well test conditions and methods for analyzing well test data from fractured media are investigated. Several analytical models are developed to be used for analyzing well test data from fractured media. Numerical tools that may be used to simulate fluid flow in fractured media are also presented. Three types of composite models for constant flux tests are investigated. Several slug test models with different geometric conditions that may be present in fractured media are also investigated. A finite element model that can simulate transient fluid flow in fracture networks is used to study the behavior of various two-dimensional fracture systems under well test conditions. A mesh generator that can be used to model mass and heat flow in a fractured-porous media is presented. This model develops an explicit solution in the porous matrix as well as in the discrete fractures. Because the model does not require the assumptions of the conventional double porosity approach, it may be used to simulate cases where double porosity models fail.
Creating 3-dimensional Models of the Photosphere using the SIR Code
NASA Astrophysics Data System (ADS)
Thonhofer, S.; Utz, D.; Jurčák, J.; Pauritsch, J.; Hanslmeier, A.; Lemmerer, B.
A high-resolution 3-dimensional model of the photospheric magnetic field is essential for the investigation of magnetic features such as sunspots, pores or smaller elements like single flux tubes seen as magnetic bright points. The SIR code is an advanced inversion code that retrieves physical quantities, e.g. magnetic field, from Stokes profiles. Based on this code, we developed a program for automated inversion of Hinode SOT/SP data and for storing these results in 3-dimensional data cubes in the form of fits files. We obtained models of the temperature, magnetic field strength, magnetic field angles and LOS-velocity in a region of the quiet sun. We will give a first discussion of those parameters in regards of small scale magnetic fields and what we can obtain and learn in the future.
Collawn, Sherry S; Mobley, James A; Banerjee, N Sanjib; Chow, Louise T
2016-04-01
Wound healing involves a number of factors that results in the production of a "closed" wound. Studies have shown, in animal models, acceleration of wound healing with the addition of adipose-derived stromal cells (ADSC). The cause for the positive effect which these cells have on wound healing has not been elucidated. We have previously shown that addition of ADSC to the dermal equivalent in 3-dimensional skin cultures accelerates reepithelialization. We now demonstrate that conditioned media (CM) from cultured ADSC produced a similar rate of healing. This result suggests that a feedback from the 3-dimensional epithelial cultures to ADSC was not necessary to effect the accelerated reepithelialization. Mass spectrometry of CM from ADSC and primary human fibroblasts revealed differences in secretomes, some of which might have roles in the accelerating wound healing. Thus, the use of CM has provided some preliminary information on a possible mode of action. PMID:26954733
Towards non-AdS holography in 3-dimensional higher spin gravity
NASA Astrophysics Data System (ADS)
Gary, Michael; Grumiller, Daniel; Rashkov, Radoslav
2012-03-01
We take the first steps towards non-AdS holography in higher spin gravity. Namely, we propose a variational principle for generic 3-dimensional higher spin gravity that accommodates asymptotic backgrounds beyond AdS, like asymptotically Schrödinger, Lifshitz or warped AdS spacetimes. As examples we study in some detail the four sl(2) embeddings of spin-4 gravity and provide associated geometries, including an asymptotic Lifshitz black hole.
Energy Sources of the Dominant Frequency Dependent 3-dimensional Atmospheric Modes
NASA Technical Reports Server (NTRS)
Schubert, S.
1985-01-01
The energy sources and sinks associated with the zonally asymmetric winter mean flow are investigated as part of an on-going study of atmospheric variability. Distinctly different horizontal structures for the long, intermediate and short time scale atmospheric variations were noted. In previous observations, the 3-dimensional structure of the fluctuations is investigated and the relative roles of barotropic and baroclinic terms are assessed.
Pediatric Thighbone (Femur) Fracture
... fractures in infants under 1 year old is child abuse. Child abuse is also a leading cause of thighbone fracture ... contact sports • Being in a motor vehicle accident • Child abuse Types of Femur Fractures (Classification) Femur fractures vary ...
Not Available
1990-01-01
This book contains papers on hydraulic fracturing. Topics covered include: An overview of recent advances in hydraulic fracturing technology; Containment of massive hydraulic fracture; and Fracturing with a high-strength proppant.
NASA Astrophysics Data System (ADS)
Zou, Yushi; Zhang, Shicheng; Ma, Xinfang; Zhou, Tong; Zeng, Bo
2016-03-01
Hydraulic fracture network (HFN) propagation in naturally fractured shale formations is investigated numerically using a 3D complex fracturing model based on the discrete element method. To account for the plastic deformation behavior of shales, the Drucker-Prager plasticity model is incorporated into the fracturing model. Parametric studies are then conducted for different Young's moduli, horizontal differential stresses, natural fracture (NF) properties, injection rates, and number and spacing of perforation clusters. Numerical results show that horizontal differential stress primarily determines the generation of a complex HFN. The plastic deformation of shale can reduce the stimulated reservoir volume; this is more obvious with Young's modulus of less than 20 GPa. In addition, a higher injection rate could largely increase the fracture complexity index (FCI). Moreover, increasing perforation cluster numbers per fracturing stage is beneficial for increasing the FCI, but it also increases the potential merging of neighboring fractures, which may lead to non-uniform development of HFN in far-wellbore regions. To achieve uniform development of HFN within a fracturing stage, the distribution of NFs should be fully considered. The results presented here may provide improved understanding of HFN generation and are favorable for optimizing fracturing treatment designs for shale formations.
Fracture process zone in granite
Zang, A.; Wagner, F.C.; Stanchits, S.; Janssen, C.; Dresen, G.
2000-01-01
In uniaxial compression tests performed on Aue granite cores (diameter 50 mm, length 100 mm), a steel loading plate was used to induce the formation of a discrete shear fracture. A zone of distributed microcracks surrounds the tip of the propagating fracture. This process zone is imaged by locating acoustic emission events using 12 piezoceramic sensors attached to the samples. Propagation velocity of the process zone is varied by using the rate of acoustic emissions to control the applied axial force. The resulting velocities range from 2 mm/s in displacement-controlled tests to 2 ??m/s in tests controlled by acoustic emission rate. Wave velocities and amplitudes are monitored during fault formation. P waves transmitted through the approaching process zone show a drop in amplitude of 26 dB, and ultrasonic velocities are reduced by 10%. The width of the process zone is ???9 times the grain diameter inferred from acoustic data but is only 2 times the grain size from optical crack inspection. The process zone of fast propagating fractures is wider than for slow ones. The density of microcracks and acoustic emissions increases approaching the main fracture. Shear displacement scales linearly with fracture length. Fault plane solutions from acoustic events show similar orientation of nodal planes on both sides of the shear fracture. The ratio of the process zone width to the fault length in Aue granite ranges from 0.01 to 0.1 inferred from crack data and acoustic emissions, respectively. The fracture surface energy is estimated from microstructure analysis to be ???2 J. A lower bound estimate for the energy dissipated by acoustic events is 0.1 J. Copyright 2000 by the American Geophysical Union.
Discreteness induced extinction
NASA Astrophysics Data System (ADS)
dos Santos, Renato Vieira; da Silva, Linaena Méricy
2015-11-01
Two simple models based on ecological problems are discussed from the point of view of non-equilibrium statistical mechanics. It is shown how discrepant may be the results of the models that include spatial distribution with discrete interactions when compared with the continuous analogous models. In the continuous case we have, under certain circumstances, the population explosion. When we take into account the finiteness of the population, we get the opposite result, extinction. We will analyze how these results depend on the dimension d of the space and describe the phenomenon of the "Discreteness Inducing Extinction" (DIE). The results are interpreted in the context of the "paradox of sex", an old problem of evolutionary biology.
Su, Yuliang; Ren, Long; Meng, Fankun; Xu, Chen; Wang, Wendong
2015-01-01
Stimulated reservoir volume (SRV) fracturing in tight oil reservoirs often induces complex fracture-network growth, which has a fundamentally different formation mechanism from traditional planar bi-winged fracturing. To reveal the mechanism of fracture network propagation, this paper employs a modified displacement discontinuity method (DDM), mechanical mechanism analysis and initiation and propagation criteria for the theoretical model of fracture network propagation and its derivation. A reasonable solution of the theoretical model for a tight oil reservoir is obtained and verified by a numerical discrete method. Through theoretical calculation and computer programming, the variation rules of formation stress fields, hydraulic fracture propagation patterns (FPP) and branch fracture propagation angles and pressures are analyzed. The results show that during the process of fracture propagation, the initial orientation of the principal stress deflects, and the stress fields at the fracture tips change dramatically in the region surrounding the fracture. Whether the ideal fracture network can be produced depends on the geological conditions and on the engineering treatments. This study has both theoretical significance and practical application value by contributing to a better understanding of fracture network propagation mechanisms in unconventional oil/gas reservoirs and to the improvement of the science and design efficiency of reservoir fracturing. PMID:25966285
Su, Yuliang; Ren, Long; Meng, Fankun; Xu, Chen; Wang, Wendong
2015-01-01
Stimulated reservoir volume (SRV) fracturing in tight oil reservoirs often induces complex fracture-network growth, which has a fundamentally different formation mechanism from traditional planar bi-winged fracturing. To reveal the mechanism of fracture network propagation, this paper employs a modified displacement discontinuity method (DDM), mechanical mechanism analysis and initiation and propagation criteria for the theoretical model of fracture network propagation and its derivation. A reasonable solution of the theoretical model for a tight oil reservoir is obtained and verified by a numerical discrete method. Through theoretical calculation and computer programming, the variation rules of formation stress fields, hydraulic fracture propagation patterns (FPP) and branch fracture propagation angles and pressures are analyzed. The results show that during the process of fracture propagation, the initial orientation of the principal stress deflects, and the stress fields at the fracture tips change dramatically in the region surrounding the fracture. Whether the ideal fracture network can be produced depends on the geological conditions and on the engineering treatments. This study has both theoretical significance and practical application value by contributing to a better understanding of fracture network propagation mechanisms in unconventional oil/gas reservoirs and to the improvement of the science and design efficiency of reservoir fracturing. PMID:25966285
A paradigm for discrete physics
Noyes, H.P.; McGoveran, D.; Etter, T.; Manthey, M.J.; Gefwert, C.
1987-01-01
An example is outlined for constructing a discrete physics using as a starting point the insight from quantum physics that events are discrete, indivisible and non-local. Initial postulates are finiteness, discreteness, finite computability, absolute nonuniqueness (i.e., homogeneity in the absence of specific cause) and additivity.
Efficient and robust compositional two-phase reservoir simulation in fractured media
NASA Astrophysics Data System (ADS)
Zidane, A.; Firoozabadi, A.
2015-12-01
Compositional and compressible two-phase flow in fractured media has wide applications including CO2 injection. Accurate simulations are currently based on the discrete fracture approach using the cross-flow equilibrium model. In this approach the fractures and a small part of the matrix blocks are combined to form a grid cell. The major drawback is low computational efficiency. In this work we use the discrete-fracture approach to model the fractures where the fracture entities are described explicitly in the computational domain. We use the concept of cross-flow equilibrium in the fractures (FCFE). This allows using large matrix elements in the neighborhood of the fractures. We solve the fracture transport equations implicitly to overcome the Courant-Freidricks-Levy (CFL) condition in the small fracture elements. Our implicit approach is based on calculation of the derivative of the molar concentration of component i in phase (cαi ) with respect to the total molar concentration (ci ) at constant volume V and temperature T. This contributes to significant speed up of the code. The hybrid mixed finite element method (MFE) is used to solve for the velocity in both the matrix and the fractures coupled with the discontinuous Galerkin (DG) method to solve the species transport equations in the matrix, and a finite volume (FV) discretization in the fractures. In large scale problems the proposed approach is orders of magnitude faster than the existing models.
Fracture sealing utilising microbially induced carbonate precipitation
NASA Astrophysics Data System (ADS)
Moir, H.; El Mountassir, G.; Lunn, R. J.; Gilfillan, L. J.
2011-12-01
is focussed along a gradually decreasing number of discrete tortuous channels (within the fracture plane). Experiments also show that above a certain threshold velocity, precipitation can no longer occur. As a consequence, once the number and aperture of channels within a fracture plane is sufficiently small, the fluid velocity exceeds this threshold and a steady state is achieved i.e. no further precipitation occurs. These experiments have important consequences for understanding fluid flow and contaminant migration within the crust. They imply that, for bio-mineral processes at least, after each mechanical failure of the crust, regional fracture networks may gradually reduce to form a small number of discrete stable pathways within which the velocity exceeds a given threshold value that prohibits further mineralisation.
Chua, Julianto; Mathews, Nripan; Jennings, James R; Yang, Guangwu; Wang, Qing; Mhaisalkar, Subodh G
2011-11-21
We describe the application of 3-dimensional metal grid electrodes (3D-MGEs) as electron collectors in dye-sensitized solar cells (DSCs) as a replacement for fluorinated tin oxide (FTO) electrodes. Requirements, structure, advantages, and limitations of the metal grid electrodes are discussed. Solar conversion efficiencies of 6.2% have been achieved in 3D-MGE based solar cells, comparable to that fabricated on FTO (7.1%). The charge transport properties and collection efficiencies in these novel solar cells have been studied using electrochemical impedance spectroscopy. PMID:21989708
Incorporating a 3-dimensional printer into the management of early-stage cervical cancer.
Baek, Min-Hyun; Kim, Dae-Yeon; Kim, Namkug; Rhim, Chae Chun; Kim, Jong-Hyeok; Nam, Joo-Hyun
2016-08-01
We used a 3-dimensional (3D) printer to create anatomical replicas of real lesions and tested its application in cervical cancer. Our study patient decided to undergo radical hysterectomy after seeing her 3D model which was then used to plan and simulate this surgery. Using 3D printers to create patient-specific 3D tumor models may aid cervical cancer patients make treatment decisions. This technology will lead to better surgical and oncological outcomes for cervical cancer patients. J. Surg. Oncol. 2016;114:150-152. © 2016 Wiley Periodicals, Inc. PMID:27222318
MT3D: a 3 dimensional magnetotelluric modeling program (user's guide and documentation for Rev. 1)
Nutter, C.; Wannamaker, P.E.
1980-11-01
MT3D.REV1 is a non-interactive computer program written in FORTRAN to do 3-dimensional magnetotelluric modeling. A 3-D volume integral equation has been adapted to simulate the MT response of a 3D body in the earth. An integro-difference scheme has been incorporated to increase the accuracy. This is a user's guide for MT3D.REV1 on the University of Utah Research Institute's (UURI) PRIME 400 computer operating under PRIMOS IV, Rev. 17.
International "Intercomparison of 3-Dimensional (3D) Radiation Codes" (13RC)
NASA Technical Reports Server (NTRS)
Cahalan, Robert F.; Einaudi, Franco (Technical Monitor)
2000-01-01
An international "Intercomparison of 3-dimensional (3D) Radiation Codes" 13RC) has been initiated. It is endorsed by the GEWEX Radiation Panel, and funded jointly by the United States Department of Energy ARM program, and by the National Aeronautics and Space Administration Radiation Sciences program. It is a 3-phase effort that has as its goals to: (1) understand the errors and limits of 3D methods; (2) provide 'baseline' cases for future 3D code development; (3) promote sharing of 3D tools; (4) derive guidelines for 3D tool selection; and (5) improve atmospheric science education in 3D radiation.
High-speed 3-dimensional imaging in robot-assisted thoracic surgical procedures.
Kajiwara, Naohiro; Akata, Soichi; Hagiwara, Masaru; Yoshida, Koichi; Kato, Yasufumi; Kakihana, Masatoshi; Ohira, Tatsuo; Kawate, Norihiko; Ikeda, Norihiko
2014-06-01
We used a high-speed 3-dimensional (3D) image analysis system (SYNAPSE VINCENT, Fujifilm Corp, Tokyo, Japan) to determine the best positioning of robotic arms and instruments preoperatively. The da Vinci S (Intuitive Surgical Inc, Sunnyvale, CA) was easily set up accurately and rapidly for this operation. Preoperative simulation and intraoperative navigation using the SYNAPSE VINCENT for robot-assisted thoracic operations enabled efficient planning of the operation settings. The SYNAPSE VINCENT can detect the tumor location and depict surrounding tissues quickly, accurately, and safely. This system is also excellent for navigational and educational use. PMID:24882302
NASA Technical Reports Server (NTRS)
Kubrynski, Krzysztof
1991-01-01
A subcritical panel method applied to flow analysis and aerodynamic design of complex aircraft configurations is presented. The analysis method is based on linearized, compressible, subsonic flow equations and indirect Dirichlet boundary conditions. Quadratic dipol and linear source distribution on flat panels are applied. In the case of aerodynamic design, the geometry which minimizes differences between design and actual pressure distribution is found iteratively, using numerical optimization technique. Geometry modifications are modeled by surface transpiration concept. Constraints in respect to resulting geometry can be specified. A number of complex 3-dimensional design examples are presented. The software is adopted to personal computers, and as result an unexpected low cost of computations is obtained.
A 3-dimensional finite-difference method for calculating the dynamic coefficients of seals
NASA Technical Reports Server (NTRS)
Dietzen, F. J.; Nordmann, R.
1989-01-01
A method to calculate the dynamic coefficients of seals with arbitrary geometry is presented. The Navier-Stokes equations are used in conjunction with the k-e turbulence model to describe the turbulent flow. These equations are solved by a full 3-dimensional finite-difference procedure instead of the normally used perturbation analysis. The time dependence of the equations is introduced by working with a coordinate system rotating with the precession frequency of the shaft. The results of this theory are compared with coefficients calculated by a perturbation analysis and with experimental results.
NASA Astrophysics Data System (ADS)
Kotulski, Zbigniew; Szczepaski, Janusz
In the paper we propose a new method of constructing cryptosystems utilising a nonpredictability property of discrete chaotic systems. We formulate the requirements for such systems to assure their safety. We also give examples of practical realisation of chaotic cryptosystems, using a generalisation of the method presented in [7]. The proposed algorithm of encryption and decryption is based on multiple iteration of a certain dynamical chaotic system. We assume that some part of the initial condition is a plain message. As the secret key we assume the system parameter(s) and additionally another part of the initial condition.
NASA Astrophysics Data System (ADS)
Wuensche, Andrew
DDLab is interactive graphics software for creating, visualizing, and analyzing many aspects of Cellular Automata, Random Boolean Networks, and Discrete Dynamical Networks in general and studying their behavior, both from the time-series perspective — space-time patterns, and from the state-space perspective — attractor basins. DDLab is relevant to research, applications, and education in the fields of complexity, self-organization, emergent phenomena, chaos, collision-based computing, neural networks, content addressable memory, genetic regulatory networks, dynamical encryption, generative art and music, and the study of the abstract mathematical/physical/dynamical phenomena in their own right.
Breakdown of Continuum Fracture Mechanics at the Nanoscale
Shimada, Takahiro; Ouchi, Kenji; Chihara, Yuu; Kitamura, Takayuki
2015-01-01
Materials fail by the nucleation and propagation of a crack, the critical condition of which is quantitatively described by fracture mechanics that uses an intensity of singular stress field characteristically formed near the crack-tip. However, the continuum assumption basing fracture mechanics obscures the prediction of failure of materials at the nanoscale due to discreteness of atoms. Here, we demonstrate the ultimate dimensional limit of fracture mechanics at the nanoscale, where only a small number of atoms are included in a singular field of continuum stress formed near a crack tip. Surprisingly, a singular stress field of only several nanometers still governs fracture as successfully as that at the macroscale, whereas both the stress intensity factor and the energy release rate fail to describe fracture below a critically confined singular field of 2–3 nm, i.e., breakdown of fracture mechanics within the framework of the continuum theory. We further propose an energy-based theory that explicitly accounts for the discrete nature of atoms, and demonstrate that our theory not only successfully describes fracture even below the critical size but also seamlessly connects the atomic to macroscales. It thus provides a more universal fracture criterion, and novel atomistic insights into fracture. PMID:25716684
Breakdown of Continuum Fracture Mechanics at the Nanoscale
NASA Astrophysics Data System (ADS)
Shimada, Takahiro; Ouchi, Kenji; Chihara, Yuu; Kitamura, Takayuki
2015-02-01
Materials fail by the nucleation and propagation of a crack, the critical condition of which is quantitatively described by fracture mechanics that uses an intensity of singular stress field characteristically formed near the crack-tip. However, the continuum assumption basing fracture mechanics obscures the prediction of failure of materials at the nanoscale due to discreteness of atoms. Here, we demonstrate the ultimate dimensional limit of fracture mechanics at the nanoscale, where only a small number of atoms are included in a singular field of continuum stress formed near a crack tip. Surprisingly, a singular stress field of only several nanometers still governs fracture as successfully as that at the macroscale, whereas both the stress intensity factor and the energy release rate fail to describe fracture below a critically confined singular field of 2-3 nm, i.e., breakdown of fracture mechanics within the framework of the continuum theory. We further propose an energy-based theory that explicitly accounts for the discrete nature of atoms, and demonstrate that our theory not only successfully describes fracture even below the critical size but also seamlessly connects the atomic to macroscales. It thus provides a more universal fracture criterion, and novel atomistic insights into fracture.
NASA Astrophysics Data System (ADS)
Calogero, Francesco
2011-08-01
The original continuous-time ''goldfish'' dynamical system is characterized by two neat formulas, the first of which provides the N Newtonian equations of motion of this dynamical system, while the second provides the solution of the corresponding initial-value problem. Several other, more general, solvable dynamical systems ''of goldfish type'' have been identified over time, featuring, in the right-hand (''forces'') side of their Newtonian equations of motion, in addition to other contributions, a velocity-dependent term such as that appearing in the right-hand side of the first formula mentioned above. The solvable character of these models allows detailed analyses of their behavior, which in some cases is quite remarkable (for instance isochronous or asymptotically isochronous). In this paper we introduce and discuss various discrete-time dynamical systems, which are as well solvable, which also display interesting behaviors (including isochrony and asymptotic isochrony) and which reduce to dynamical systems of goldfish type in the limit when the discrete-time independent variable l=0,1,2,... becomes the standard continuous-time independent variable t, 0≤t<∞.
NASA Astrophysics Data System (ADS)
Rahn, Helene; Alexiou, Christoph; Trahms, Lutz; Odenbach, Stefan
2014-06-01
X-ray computed tomography is nowadays used for a wide range of applications in medicine, science and technology. X-ray microcomputed tomography (XµCT) follows the same principles used for conventional medical CT scanners, but improves the spatial resolution to a few micrometers. We present an example of an application of X-ray microtomography, a study of 3-dimensional biodistribution, as along with the quantification of nanoparticle content in tumoral tissue after minimally invasive cancer therapy. One of these minimal invasive cancer treatments is magnetic drug targeting, where the magnetic nanoparticles are used as controllable drug carriers. The quantification is based on a calibration of the XµCT-equipment. The developed calibration procedure of the X-ray-µCT-equipment is based on a phantom system which allows the discrimination between the various gray values of the data set. These phantoms consist of a biological tissue substitute and magnetic nanoparticles. The phantoms have been studied with XµCT and have been examined magnetically. The obtained gray values and nanoparticle concentration lead to a calibration curve. This curve can be applied to tomographic data sets. Accordingly, this calibration enables a voxel-wise assignment of gray values in the digital tomographic data set to nanoparticle content. Thus, the calibration procedure enables a 3-dimensional study of nanoparticle distribution as well as concentration.
Feng, Guodong; Zhao, Yang; Tian, Xu; Gao, Zhiqiang
2014-01-01
This study aimed to establish a 3-dimensional dynamic quantitative facial motion analysis system, and then determine its accuracy and test-retest reliability. The system could automatically reconstruct the motion of the observational points. Standardized T-shaped rod and L-shaped rods were used to evaluate the static and dynamic accuracy of the system. Nineteen healthy volunteers were recruited to test the reliability of the system. The average static distance error measurement was 0.19 mm, and the average angular error was 0.29°. The measuring results decreased with the increase of distance between the cameras and objects, 80 cm of which was considered to be optimal. It took only 58 seconds to perform the full facial measurement process. The average intra-class correlation coefficient for distance measurement and angular measurement was 0.973 and 0.794 respectively. The results demonstrated that we successfully established a practical 3-dimensional dynamic quantitative analysis system that is accurate and reliable enough to meet both clinical and research needs. PMID:25390881
A customizable 3-dimensional digital atlas of the canary brain in multiple modalities.
Vellema, Michiel; Verschueren, Jacob; Van Meir, Vincent; Van der Linden, Annemie
2011-07-15
Songbirds are well known for their ability to learn their vocalizations by imitating conspecific adults. This uncommon skill has led to many studies examining the behavioral and neurobiological processes involved in vocal learning. Canaries display a variable, seasonally dependent, vocal behavior throughout their lives. This trait makes this bird species particularly valuable to study the functional relationship between the continued plasticity in the singing behavior and alterations in the anatomy and physiology of the brain. In order to optimally interpret these types of studies, a detailed understanding of the brain anatomy is essential. Because traditional 2-dimensional brain atlases are limited in the information they can provide about the anatomy of the brain, here we present a 3-dimensional MRI-based atlas of the canary brain. Using multiple imaging protocols we were able to maximize the number of detectable brain regions, including most of the areas involved in song perception, learning, and production. The brain atlas can readily be used to determine the stereotactic location of delineated brain areas at any desirable head angle. Alternatively the brain data can be used to determine the ideal orientation of the brain for stereotactic injections, electrophysiological recordings, and brain sectioning. The 3-dimensional canary brain atlas presented here is freely available and is easily adaptable to support many types of neurobiological studies, including anatomical, electrophysiological, histological, explant, and tracer studies. PMID:21565273
Realization of integral 3-dimensional image using fabricated tunable liquid lens array
NASA Astrophysics Data System (ADS)
Lee, Muyoung; Kim, Junoh; Kim, Cheol Joong; Lee, Jin Su; Won, Yong Hyub
2015-03-01
Electrowetting has been widely studied for various optical applications such as optical switch, sensor, prism, and display. In this study, vari-focal liquid lens array is developed using electrowetting principle to construct integral 3-dimensional imaging. The electrowetting principle that changes the surface tension by applying voltage has several advantages to realize active optical device such as fast response time, low electrical consumption, and no mechanical moving parts. Two immiscible liquids that are water and oil are used for forming lens. By applying a voltage to the water, the focal length of the lens could be tuned as changing contact angle of water. The fabricated electrowetting vari-focal liquid lens array has 1mm diameter spherical lens shape that has 1.6mm distance between each lens. The number of lenses on the panel is 23x23 and the focal length of the lens array is simultaneously tuned from -125 to 110 diopters depending on the applied voltage. The fabricated lens array is implemented to integral 3-dimensional imaging. A 3D object is reconstructed by fabricated liquid lens array with 23x23 elemental images that are generated by 3D max tools. When liquid lens array is tuned as convex state. From vari-focal liquid lens array implemented integral imaging system, we expect that depth enhanced integral imaging can be realized in the near future.
Morphological Control of Cells on 3-Dimensional Multi-Layer Nanotopographic Structures.
Jeong, Heon-Ho; Noh, Young-Mu; Song, Hwan-Moon; Lee, Sang-Ho; Park, Jin-Sung; Lee, Chang-Soo
2015-05-01
The extracellular matrix (ECM) environment is known to play an important role in the process of various cell regulatory mechanisms. We have investigated the ability of 3-dimensional ECM geometries to induce morphological changes in cells. Bi-layer polymeric structures with submicron scale stripe patterns were fabricated using a two-step nano-imprinting technique, and the orientation angle (θ(α)) of the upper layer was controlled by changing its alignment with respect to the orientation of the bottom layer. When cells were grown on the mono-layer stripe structure with a single orientation, they elongated along the direction of the stripe pattern. On bi-layer polymer structures, the cell morphologies gradually changed and became rounded, with an increase of θα up to 90 degrees, but the polarities of these cells were still aligned along the orientation of the upper layer. As a result, we show that the polarity and the roundness of cells can be independently regulated by adjusting the orientation of 3-dimensional hierarchical ECM topography. PMID:26505024
Particle trajectory computation on a 3-dimensional engine inlet. Final Report Ph.D. Thesis
NASA Technical Reports Server (NTRS)
Kim, J. J.
1986-01-01
A 3-dimensional particle trajectory computer code was developed to compute the distribution of water droplet impingement efficiency on a 3-dimensional engine inlet. The computed results provide the essential droplet impingement data required for the engine inlet anti-icing system design and analysis. The droplet trajectories are obtained by solving the trajectory equation using the fourth order Runge-Kutta and Adams predictor-corrector schemes. A compressible 3-D full potential flow code is employed to obtain a cylindrical grid definition of the flowfield on and about the engine inlet. The inlet surface is defined mathematically through a system of bi-cubic parametric patches in order to compute the droplet impingement points accurately. Analysis results of the 3-D trajectory code obtained for an axisymmetric droplet impingement problem are in good agreement with NACA experimental data. Experimental data are not yet available for the engine inlet impingement problem analyzed. Applicability of the method to solid particle impingement problems, such as engine sand ingestion, is also demonstrated.
Galland-Girodet, Sigolène; Pashtan, Itai; MacDonald, Shannon M.; Ancukiewicz, Marek; Hirsch, Ariel E.; Kachnic, Lisa A.; Specht, Michelle; Gadd, Michele; Smith, Barbara L.; Powell, Simon N.; Recht, Abram; Taghian, Alphonse G.
2014-11-01
Purpose: To present long-term outcomes of a prospective feasibility trial using either protons or 3-dimensional conformal photon-based (accelerated partial-breast irradiation [APBI]) techniques. Methods and Materials: From October 2003 to April 2006, 98 evaluable patients with stage I breast cancer were treated with APBI (32 Gy in 8 fractions given twice daily) on a prospective clinical trial: 19 with proton beam therapy (PBT) and 79 with photons or mixed photons/electrons. Median follow-up was 82.5 months (range, 2-104 months). Toxicity and patient satisfaction evaluations were performed at each visit. Results: At 7 years, the physician rating of overall cosmesis was good or excellent for 62% of PBT patients, compared with 94% for photon patients (P=.03). Skin toxicities were more common for the PBT group: telangiectasia, 69% and 16% (P=.0013); pigmentation changes, 54% and 22% (P=.02); and other late skin toxicities, 62% and 18% (P=.029) for PBT and photons, respectively. There were no significant differences between the groups in the incidences of breast pain, edema, fibrosis, fat necrosis, skin desquamation, and rib pain or fracture. Patient-reported cosmetic outcomes at 7 years were good or excellent for 92% and 96% of PBT and photon patients, respectively (P=.95). Overall patient satisfaction was 93% for the entire cohort. The 7-year local failure rate for all patients was 6%, with 3 local recurrences in the PBT group (7-year rate, 11%) and 2 in photon-treated patients (4%) (P=.22). Conclusions: Local failure rates of 3-dimensional APBI and PBT were similar in this study. However, PBT, as delivered in this study, led to higher rates of long-term telangiectasia, skin color changes, and skin toxicities. We recommend the use of multiple fields and treatment of all fields per treatment session or the use of scanning techniques to minimize skin toxicity.
Noyes, H.P. ); Starson, S. )
1991-03-01
Discrete physics, because it replaces time evolution generated by the energy operator with a global bit-string generator (program universe) and replaces fields'' with the relativistic Wheeler-Feynman action at a distance,'' allows the consistent formulation of the concept of signed gravitational charge for massive particles. The resulting prediction made by this version of the theory is that free anti-particles near the surface of the earth will fall'' up with the same acceleration that the corresponding particles fall down. So far as we can see, no current experimental information is in conflict with this prediction of our theory. The experiment crusis will be one of the anti-proton or anti-hydrogen experiments at CERN. Our prediction should be much easier to test than the small effects which those experiments are currently designed to detect or bound. 23 refs.
Discrete Sibson interpolation.
Park, Sung W; Linsen, Lars; Kreylos, Oliver; Owens, John D; Hamann, Bernd
2006-01-01
Natural-neighbor interpolation methods, such as Sibson's method, are well-known schemes for multivariate data fitting and reconstruction. Despite its many desirable properties, Sibson's method is computationally expensive and difficult to implement, especially when applied to higher-dimensional data. The main reason for both problems is the method's implementation based on a Voronoi diagram of all data points. We describe a discrete approach to evaluating Sibson's interpolant on a regular grid, based solely on finding nearest neighbors and rendering and blending d-dimensional spheres. Our approach does not require us to construct an explicit Voronoi diagram, is easily implemented using commodity three-dimensional graphics hardware, leads to a significant speed increase compared to traditional approaches, and generalizes easily to higher dimensions. For large scattered data sets, we achieve two-dimensional (2D) interpolation at interactive rates and 3D interpolation (3D) with computation times of a few seconds. PMID:16509383
Discrete Pearson distributions
Bowman, K.O.; Shenton, L.R.; Kastenbaum, M.A.
1991-11-01
These distributions are generated by a first order recursive scheme which equates the ratio of successive probabilities to the ratio of two corresponding quadratics. The use of a linearized form of this model will produce equations in the unknowns matched by an appropriate set of moments (assumed to exist). Given the moments we may find valid solutions. These are two cases; (1) distributions defined on the non-negative integers (finite or infinite) and (2) distributions defined on negative integers as well. For (1), given the first four moments, it is possible to set this up as equations of finite or infinite degree in the probability of a zero occurrence, the sth component being a product of s ratios of linear forms in this probability in general. For (2) the equation for the zero probability is purely linear but may involve slowly converging series; here a particular case is the discrete normal. Regions of validity are being studied. 11 refs.
Immigration and Prosecutorial Discretion
Apollonio, Dorie; Lochner, Todd; Heddens, Myriah
2015-01-01
Immigration has become an increasingly salient national issue in the US, and the Department of Justice recently increased federal efforts to prosecute immigration offenses. This shift, however, relies on the cooperation of US attorneys and their assistants. Traditionally federal prosecutors have enjoyed enormous discretion and have been responsive to local concerns. To consider how the centralized goal of immigration enforcement may have influenced federal prosecutors in regional offices, we review their prosecution of immigration offenses in California using over a decade's worth of data. Our findings suggest that although centralizing forces influence immigration prosecutions, individual US attorneys' offices retain distinct characteristics. Local factors influence federal prosecutors' behavior in different ways depending on the office. Contrary to expectations, unemployment rates did not affect prosecutors' willingness to pursue immigration offenses, nor did local popular opinion about illegal immigration. PMID:26146530
Discrete stability in stochastic programming
Lepp, R.
1994-12-31
In this lecture we study stability properties of stochastic programs with recourse where the probability measure is approximated by a sequence of weakly convergent discrete measures. Such discrete approximation approach gives us a possibility to analyze explicitly the behavior of the second stage correction function. The approach is based on modern functional analytical methods of an approximation of extremum problems in function spaces, especially on the notion of the discrete convergence of vectors to an essentially bounded measurable function.
Estimating the hydraulic conductivity of two-dimensional fracture networks
NASA Astrophysics Data System (ADS)
Leung, C. T.; Zimmerman, R. W.
2010-12-01
Most oil and gas reservoirs, as well as most potential sites for nuclear waste disposal, are naturally fractured. In these sites, the network of fractures will provide the main path for fluid to flow through the rock mass. In many cases, the fracture density is so high as to make it impractical to model it with a discrete fracture network (DFN) approach. For such rock masses, it would be useful to have recourse to analytical, or semi-analytical, methods to estimate the macroscopic hydraulic conductivity of the fracture network. We have investigated single-phase fluid flow through stochastically generated two-dimensional fracture networks. The centres and orientations of the fractures are uniformly distributed, whereas their lengths follow either a lognormal distribution or a power law distribution. We have considered the case where the fractures in the network each have the same aperture, as well as the case where the aperture of each fracture is directly proportional to the fracture length. The discrete fracture network flow and transport simulator NAPSAC, developed by Serco (Didcot, UK), is used to establish the “true” macroscopic hydraulic conductivity of the network. We then attempt to match this conductivity using a simple estimation method that does not require extensive computation. For our calculations, fracture networks are represented as networks composed of conducting segments (bonds) between nodes. Each bond represents the region of a single fracture between two adjacent intersections with other fractures. We assume that the bonds are arranged on a kagome lattice, with some fraction of the bonds randomly missing. The conductance of each bond is then replaced with some effective conductance, Ceff, which we take to be the arithmetic mean of the individual conductances, averaged over each bond, rather than over each fracture. This is in contrast to the usual approximation used in effective medium theories, wherein the geometric mean is used. Our
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
Radial head fracture - aftercare
Elbow fracture - radial head - aftercare ... the radius bone, just below your elbow. A fracture is a break in your bone. The most common cause of a radial head fracture is falling with an outstretched arm.
... this page: //medlineplus.gov/ency/patientinstructions/000552.htm Hand fracture - aftercare To use the sharing features on ... need to be repaired with surgery. Types of hand fractures Your fracture may be in one of ...
NASA Astrophysics Data System (ADS)
Roubinet, D.; Linde, N.; Jougnot, D.; Irving, J.
2016-05-01
Numerous field experiments suggest that the self-potential (SP) geophysical method may allow for the detection of hydraulically active fractures and provide information about fracture properties. However, a lack of suitable numerical tools for modeling streaming potentials in fractured media prevents quantitative interpretation and limits our understanding of how the SP method can be used in this regard. To address this issue, we present a highly efficient two-dimensional discrete-dual-porosity approach for solving the fluid flow and associated self-potential problems in fractured rock. Our approach is specifically designed for complex fracture networks that cannot be investigated using standard numerical methods. We then simulate SP signals associated with pumping conditions for a number of examples to show that (i) accounting for matrix fluid flow is essential for accurate SP modeling and (ii) the sensitivity of SP to hydraulically active fractures is intimately linked with fracture-matrix fluid interactions. This implies that fractures associated with strong SP amplitudes are likely to be hydraulically conductive, attracting fluid flow from the surrounding matrix.
NASA Astrophysics Data System (ADS)
Jia, Pin; Cheng, Linsong; Huang, Shijun; Wu, Yonghui
2016-06-01
This paper presents a semi-analytical model for the flow behavior of naturally fractured formations with multi-scale fracture networks. The model dynamically couples an analytical dual-porosity model with a numerical discrete fracture model. The small-scale fractures with the matrix are idealized as a dual-porosity continuum and an analytical flow solution is derived based on source functions in Laplace domain. The large-scale fractures are represented explicitly as the major fluid conduits and the flow is numerically modeled, also in Laplace domain. This approach allows us to include finer details of the fracture network characteristics while keeping the computational work manageable. For example, the large-scale fracture network may have complex geometry and varying conductivity, and the computations can be done at predetermined, discrete times, without any grids in the dual-porosity continuum. The validation of the semi-analytical model is demonstrated in comparison to the solution of ECLIPSE reservoir simulator. The simulation is fast, gridless and enables rapid model setup. On the basis of the model, we provide detailed analysis of the flow behavior of a horizontal production well in fractured reservoir with multi-scale fracture networks. The study has shown that the system may exhibit six flow regimes: large-scale fracture network linear flow, bilinear flow, small-scale fracture network linear flow, pseudosteady-state flow, interporosity flow and pseudoradial flow. During the first four flow periods, the large-scale fracture network behaves as if it only drains in the small-scale fracture network; that is, the effect of the matrix is negligibly small. The characteristics of the bilinear flow and the small-scale fracture network linear flow are predominantly determined by the dimensionless large-scale fracture conductivity. And low dimensionless fracture conductivity will generate large pressure drops in the large-scale fractures surrounding the wellbore. With
NASA Technical Reports Server (NTRS)
2005-01-01
[figure removed for brevity, see original site] Context image for PIA03084 Fractured Surface
These fractures and graben are part of Gordii Fossae, a large region that has undergone stresses which have cracked the surface.
Image information: VIS instrument. Latitude 16.6S, Longitude 234.3E. 18 meter/pixel resolution.
Note: this THEMIS visual image has not been radiometrically nor geometrically calibrated for this preliminary release. An empirical correction has been performed to remove instrumental effects. A linear shift has been applied in the cross-track and down-track direction to approximate spacecraft and planetary motion. Fully calibrated and geometrically projected images will be released through the Planetary Data System in accordance with Project policies at a later time.
NASA's Jet Propulsion Laboratory manages the 2001 Mars Odyssey mission for NASA's Office of Space Science, Washington, D.C. The Thermal Emission Imaging System (THEMIS) was developed by Arizona State University, Tempe, in collaboration with Raytheon Santa Barbara Remote Sensing. The THEMIS investigation is led by Dr. Philip Christensen at Arizona State University. Lockheed Martin Astronautics, Denver, is the prime contractor for the Odyssey project, and developed and built the orbiter. Mission operations are conducted jointly from Lockheed Martin and from JPL, a division of the California Institute of Technology in Pasadena.
3-dimensional interactive space (3DIS); A new dimension in asset protection
Veitch, S. ); Veitch, J. ); West, S.J. )
1991-01-01
This paper reports on the 3DIS security system which uses standard CCTV cameras to create 3-Dimensional detection zones around valuable assets within protected areas. An intrusion into a zone changes light values and triggers an alarm that is annunciated, while images from multiple cameras are recorded. 3DIS lowers nuisance alarm rates and provides superior automated surveillance capability. Performance is improved over 2-D systems because activity around, above or below the zone does to cause an alarm. Invisible 3-D zones protect assets as small as a pin or as large as a 747 jetliner. Detection zones are created by excising subspaces from the overlapping fields of view of two or more video cameras. Hundred of zones may co-exist, operating simultaneously. Intrusion into any 3-D zone will cause a coincidental change in light values, triggering an alarm specific to that space.
Tracking Error analysis of Concentrator Photovoltaic Module Using Total 3-Dimensional Simulator
NASA Astrophysics Data System (ADS)
Ota, Yasuyuki; Nishioka, Kensuke
2011-12-01
A 3-dimensional (3D) operating simulator for concentrator photovoltaic (CPV) module using triple-junction solar cell was developed. By connecting 3D equivalent circuit simulation for triple-junction solar cell and ray-trace simulation for optics model, the operating characteristics of CPV module were calculated. A typical flat Fresnel lens and homogenizer were adapted to the optics model. The influence of tracking error on the performance of CPV module was calculated. There was the correlation between the optical efficiency and Isc. However, Pm was not correlated with these values, and was strongly dependent on FF. We can use this total simulator for the evaluation and optimization from the light incidence to operating characteristic of CPV modules.
NASA Technical Reports Server (NTRS)
UijtdeHaag, Maarten; Thomas, Robert; Rankin, James R.
2004-01-01
The report discusses the architecture and the flight test results of a 3-Dimensional Cockpit Display of Traffic and terrain Information (3D-CDTI). The presented 3D-CDTI is a perspective display format that combines existing Synthetic Vision System (SVS) research and Automatic Dependent Surveillance-Broadcast (ADS-B) technology to improve the pilot's situational awareness. The goal of the 3D-CDTI is to contribute to the development of new display concepts for NASA's Small Aircraft Transportation System research program. Papers were presented at the PLANS 2002 meeting and the ION-GPS 2002 meeting. The contents of this report are derived from the results discussed in those papers.
NASA Astrophysics Data System (ADS)
Kushner, Mark J.; Grapperhaus, Michael J.
1996-10-01
Inductively Coupled Plasma (ICP) reactors have the potential for scaling to large area substrates while maintaining azimuthal symmetry or side-to-side uniformity across the wafer. Asymmetric etch properties in these devices have been attributed to transmission line properties of the coil, internal structures (such as wafer clamps) and non-uniform gas injection or pumping. To investigate the origins of asymmetric etch properties, a 3-dimensional hybrid model has been developed. The hybrid model contains electromagnetic, electric circuit, electron energy equation, and fluid modules. Continuity and momentum equations are solved in the fluid module along with Poisson's equation. We will discuss results for ion and radical flux uniformity to the substrate while varying the transmission line characteristics of the coil, symmetry of gas inlets/pumping, and internal structures. Comparisons will be made to expermental measurements of etch rates. ^*Work supported by SRC, NSF, ARPA/AFOSR and LAM Research.
Pei, Jimin; Grishin, Nick V.
2015-01-01
SUMMARY Multiple sequence alignment (MSA) is an essential tool with many applications in bioinformatics and computational biology. Accurate MSA construction for divergent proteins remains a difficult computational task. The constantly increasing protein sequences and structures in public databases could be used to improve alignment quality. PROMALS3D is a tool for protein MSA construction enhanced with additional evolutionary and structural information from database searches. PROMALS3D automatically identifies homologs from sequence and structure databases for input proteins, derives structure-based constraints from alignments of 3-dimensional structures, and combines them with sequence-based constraints of profile-profile alignments in a consistency-based framework to construct high-quality multiple sequence alignments. PROMALS3D output is a consensus alignment enriched with sequence and structural information about input proteins and their homologs. PROMALS3D web server and package are available at http://prodata.swmed.edu/PROMALS3D. PMID:24170408
Using 3-dimensional printing to create presurgical models for endodontic surgery.
Bahcall, James K
2014-09-01
Advances in endodontic surgery--from both a technological and procedural perspective-have been significant over the last 18 years. Although these technologies and procedural enhancements have significantly improved endodontic surgical treatment outcomes, there is still an ongoing challenge of overcoming the limitations of interpreting preoperative 2-dimensional (2-D) radiographic representation of a 3-dimensional (3-D) in vivo surgical field. Cone-beam Computed Tomography (CBCT) has helped to address this issue by providing a 3-D enhancement of the 2-D radiograph. The next logical step to further improve a presurgical case 3-D assessment is to create a surgical model from the CBCT scan. The purpose of this article is to introduce 3-D printing of CBCT scans for creating presurgical models for endodontic surgery. PMID:25197746
The program FANS-3D (finite analytic numerical simulation 3-dimensional) and its applications
NASA Technical Reports Server (NTRS)
Bravo, Ramiro H.; Chen, Ching-Jen
1992-01-01
In this study, the program named FANS-3D (Finite Analytic Numerical Simulation-3 Dimensional) is presented. FANS-3D was designed to solve problems of incompressible fluid flow and combined modes of heat transfer. It solves problems with conduction and convection modes of heat transfer in laminar flow, with provisions for radiation and turbulent flows. It can solve singular or conjugate modes of heat transfer. It also solves problems in natural convection, using the Boussinesq approximation. FANS-3D was designed to solve heat transfer problems inside one, two and three dimensional geometries that can be represented by orthogonal planes in a Cartesian coordinate system. It can solve internal and external flows using appropriate boundary conditions such as symmetric, periodic and user specified.
Evaluation of the 3-Dimensional, Weight-bearing Orientation of the Normal Adult Knee
Nam, Denis; Shah, Ritesh R.; Nunley, Ryan M.; Barrack, Robert L.
2014-01-01
The purpose of this study was to use 3-dimensional, weight-bearing images corrected for rotation to establish normative data of limb alignment and joint line orientation in asymptomatic, adult knees. One hundred adults (200 lower extremities) were recruited to receive weight-bearing, simultaneous biplanar imaging of both lower extremities. Multiple radiographic parameters were measured from 3D images, corrected for limb rotation. 70.0% of knees were in neutral, 19.5% in varus, and 10.5% in valgus overall alignment. Only 31 % of knees possessed both a neutral mechanical axis and the absence of joint line obliquity. There was substantial agreement between the 2D and 3D images for overall mechanical alignment (κ=0.77), but only a moderate agreement for joint line obliquity (κ=0.58). A substantial portion of asymptomatic adults possess either a varus or valgus mechanical alignment and joint line obliquity, PMID:24315446
Zhou, Xichun; Turchi, Craig; Wang, Denong
2009-01-01
We reported here a novel, ready-to-use bioarray platform and methodology for construction of sensitive carbohydrate cluster microarrays. This technology utilizes a 3-dimensional (3-D) poly(amidoamine) starburst dendrimer monolayer assembled on glass surface, which is functionalized with terminal aminooxy and hydrazide groups for site-specific coupling of carbohydrates. A wide range of saccharides, including monosaccharides, oligosaccharides and polysaccharides of diverse structures, are applicable for the 3-D bioarray platform without prior chemical derivatization. The process of carbohydrate coupling is effectively accelerated by microwave radiation energy. The carbohydrate concentration required for microarray fabrication is substantially reduced using this technology. Importantly, this bioarray platform presents sugar chains in defined orientation and cluster configurations. It is, thus, uniquely useful for exploration of the structural and conformational diversities of glyco-epitope and their functional properties. PMID:19791771
DIEP Flap Breast Reconstruction Using 3-dimensional Surface Imaging and a Printed Mold.
Tomita, Koichi; Yano, Kenji; Hata, Yuki; Nishibayashi, Akimitsu; Hosokawa, Ko
2015-03-01
Recent advances in 3-dimensional (3D) surface imaging technologies allow for digital quantification of complex breast tissue. We performed 11 unilateral breast reconstructions with deep inferior epigastric artery perforator (DIEP) flaps (5 immediate, 6 delayed) using 3D surface imaging for easier surgery planning and 3D-printed molds for shaping the breast neoparenchyma. A single- or double-pedicle flap was preoperatively planned according to the estimated tissue volume required and estimated total flap volume. The DIEP flap was then intraoperatively shaped with a 3D-printed mold that was based on a horizontally inverted shape of the contralateral breast. Cosmetic outcomes were assessed as satisfactory, as confirmed by the postoperative 3D measurements of bilateral breasts. We believe that DIEP flap reconstruction assisted with 3D surface imaging and a 3D-printed mold is a simple and quick method for rebuilding a symmetric breast. PMID:25878927