## You limited your search to:

**Partner:**UNT Libraries Government Documents Department

- 3-D electromagnetic modeling of wakefields in accelerator components
- We discuss the use of 3-D finite-difference time-domain (FDTD) electromagnetic codes for modeling accelerator components. Computational modeling of cylindrically symmetric structures such as induction accelerator cells has been very successful in predicting the wake potential and wake impedances of these structures, but full 3-D modeling of complex structures has been limited due to substantial computer resources required for a full 3-D model. New massively parallel 3-D time domain electromagnetic codes now under development using conforming unstructured meshes allow a substantial increase in the geometric fidelity of the structures being modeled. Development of these new codes are discussed in context of applicability to accelerator problems. Various 3-D structures are tested with an existing cubical cell FDTD code and wake impedances compared with simple analytic models for the structures; results will be used as benchmarks for testing the new time time domain codes. Structures under consideration include a stripline beam position monitor as well as circular and elliptical apertures in circular waveguides. Excellent agreement for monopole and dipole impedances with models were found for these structures below the cutoff frequency of the beam line. digital.library.unt.edu/ark:/67531/metadc685010/
- 3-D Experimental Fracture Analysis at High Temperature
- T*e, which is an elastic-plastic fracture parameter based on incremental theory of plasticity, was determined numerically and experimentally. The T*e integral of a tunneling crack in 2024-T3 aluminum, three point bend specimen was obtained through a hybrid analysis of moire interferometry and 3-D elastic-plastic finite element analysis. The results were verified by the good agreement between the experimentally and numerically determined T*e on the specimen surface. digital.library.unt.edu/ark:/67531/metadc741358/
- 3-D field computation: The near-triumph of commerical codes
- In recent years, more and more of those who design and analyze magnets and other devices are using commercial codes rather than developing their own. This paper considers the commercial codes and the features available with them. Other recent trends with 3-D field computation include parallel computation and visualization methods such as virtual reality systems. digital.library.unt.edu/ark:/67531/metadc794259/
- 3-D Finite Element Analyses of the Egan Cavern Field
- Three-dimensional finite element analyses were performed for the two gas-filled storage caverns at the Egan field, Jennings dome, Louisiana. The effects of cavern enlargement on surface subsidence, storage loss, and cavern stability were investigated. The finite element model simulated the leaching of caverns to 6 and 8 billion cubic feet (BCF) and examined their performance at various operating conditions. Operating pressures varied from 0.15 psi/ft to 0.9 psi/ft at the bottom of the lowest cemented casing. The analysis also examined the stability of the web or pillar of salt between the caverns under differential pressure loadings. The 50-year simulations were performed using JAC3D, a three dimensional finite element analysis code for nonlinear quasistatic solids. A damage criterion based on onset of dilatancy was used to evaluate cavern instability. Dilation results from the development of microfractures in salt and, hence, potential increases in permeability onset occurs well before large scale failure. The analyses predicted stable caverns throughout the 50-year period for the range of pressures investigated. Some localized salt damage was predicted near the bottom walls of the caverns if the caverns are operated at minimum pressure for long periods of time. Volumetric cavern closures over time due to creep were moderate to excessive depending on the salt creep properties and operating pressures. However, subsidence above the cavern field was small and should pose no problem, to surface facilities. digital.library.unt.edu/ark:/67531/metadc686870/
- 3-D Finite Element Analysis of Induction Logging in a Dipping Formation
- Electromagnetic induction by a magnetic dipole located above a dipping interface is of relevance to the petroleum well-logging industry. The problem is fully three-dimensional (3-D) when formulated as above, but reduces to an analytically tractable one-dimensional (1-D) problem when cast as a small tilted coil above a horizontal interface. The two problems are related by a simple coordinate rotation. An examination of the induced eddy currents and the electric charge accumulation at the interface help to explain the inductive and polarization effects commonly observed in induction logs from dipping geological formations. The equivalence between the 1-D and 3-D formulations of the problem enables the validation of a previously published finite element solver for 3-D controlled-source electromagnetic induction. digital.library.unt.edu/ark:/67531/metadc705973/
- 3-d finite element model development for biomechanics: a software demonstration
- Finite element analysis is becoming an increasingly important part of biomechanics and orthopedic research, as computational resources become more powerful, and data handling algorithms become more sophisticated. Until recently, tools with sufficient power did not exist or were not accessible to adequately model complicated, three-dimensional, nonlinear biomechanical systems. In the past, finite element analyses in biomechanics have often been limited to two-dimensional approaches, linear analyses, or simulations of single tissue types. Today, we have the resources to model fully three-dimensional, nonlinear, multi-tissue, and even multi-joint systems. The authors will present the process of developing these kinds of finite element models, using human hand and knee examples, and will demonstrate their software tools. digital.library.unt.edu/ark:/67531/metadc697585/
- 3-D Force-balanced Magnetospheric Configurations
- The knowledge of plasma pressure is essential for many physics applications in the magnetosphere, such as computing magnetospheric currents and deriving magnetosphere-ionosphere coupling. A thorough knowledge of the 3-D pressure distribution has however eluded the community, as most in-situ pressure observations are either in the ionosphere or the equatorial region of the magnetosphere. With the assumption of pressure isotropy there have been attempts to obtain the pressure at different locations by either (a) mapping observed data (e.g., in the ionosphere) along the field lines of an empirical magnetospheric field model or (b) computing a pressure profile in the equatorial plane (in 2-D) or along the Sun-Earth axis (in 1-D) that is in force balance with the magnetic stresses of an empirical model. However, the pressure distributions obtained through these methods are not in force balance with the empirical magnetic field at all locations. In order to find a global 3-D plasma pressure distribution in force balance with the magnetospheric magnetic field, we have developed the MAG-3D code, that solves the 3-D force balance equation J x B = (upside-down delta) P computationally. Our calculation is performed in a flux coordinate system in which the magnetic field is expressed in terms of Euler potentials as B = (upside-down delta) psi x (upside-down delta) alpha. The pressure distribution, P = P(psi,alpha), is prescribed in the equatorial plane and is based on satellite measurements. In addition, computational boundary conditions for y surfaces are imposed using empirical field models. Our results provide 3-D distributions of magnetic field and plasma pressure as well as parallel and transverse currents for both quiet-time and disturbed magnetospheric conditions. digital.library.unt.edu/ark:/67531/metadc738597/
- 3-D full waveform inversion of seismic data; Part I. Theory
- Full waveform inversion of seismic data is a challenging subject partly because of the lack of precise knowledge of the source. Since currently available approaches involve some form of approximations to the source, inversion results are subject to the quality and the choice of the source information used. A new full waveform inversion scheme has been introduced (Lee and Kim, 2003) using normalized wavefield for simple two-dimensional (2-D) scalar problems. The method does not require source information, so potential inversion errors due to source estimation may be eliminated. A gather of seismic traces is first Fourier-transformed into the frequency domain and a normalized wavefield is obtained for each trace in the frequency domain. Normalization is done with respect to the frequency response of a reference trace selected from the gather, so the complex-valued normalized wavefield is source-independent and dimensionless. The inversion algorithm minimizes misfits between measured normalized wavefield and numerically computed normalized wavefield. In this paper the full waveform inversion is extended to three-dimensional (3-D) problems. digital.library.unt.edu/ark:/67531/metadc735736/
- A 3-D hydrodynamic dispersion model for modeling tracer transport in Geothermal Reservoirs
- A 3-D hydrodynamic dispersion model for tracer transport is developed and implemented into the TOUGH2 EOS3 (T2R3D) module. The model formulation incorporates a full dispersion tensor, based on a 3-D velocity field with a 3-D, irregular grid in a heterogeneous geological system. Two different weighting schemes are proposed for spatial average of 3-D velocity fields and concentration gradients to evaluate the mass flux by dispersion and diffusion of a tracer or a radionuclide. This new module of the TOUGH2 code is designed to simulate processes of tracer/radionuclide transport using an irregular, 3-D integral finite difference grid in non-isothermal, three-dimensional, multiphase, porous/fractured subsurface systems. The numerical method for this transport module is based on the integral finite difference scheme, as in the TOUGH2 code. The major assumptions of the tracer transport module are: (a) a tracer or a radionuclide is present and transported only within the liquid phase, (b) transport mechanisms include molecular diffusion and hydrodynamic dispersion in the liquid phase in addition to advection, and (c) first order decay and linear adsorption on rock grains are taken into account. The tracer or radionuclide is introduced as an additional mass component into the standard TOUGH2 formulation, time is discretized fully implicitly, and non-linearities of the conservation equations are handled using the Newton/Raphson iteration. We have verified this transport module by comparison with results of a 2-D transport problem for which an analytical solution is available. In addition, a field application is described to demonstrate the use of the proposed model. digital.library.unt.edu/ark:/67531/metadc719687/
- 3-D HYDRODYNAMIC MODELING IN A GEOSPATIAL FRAMEWORK
- 3-D hydrodynamic models are used by the Savannah River National Laboratory (SRNL) to simulate the transport of thermal and radionuclide discharges in coastal estuary systems. Development of such models requires accurate bathymetry, coastline, and boundary condition data in conjunction with the ability to rapidly discretize model domains and interpolate the required geospatial data onto the domain. To facilitate rapid and accurate hydrodynamic model development, SRNL has developed a pre- and post-processor application in a geospatial framework to automate the creation of models using existing data. This automated capability allows development of very detailed models to maximize exploitation of available surface water radionuclide sample data and thermal imagery. digital.library.unt.edu/ark:/67531/metadc883870/
- 3-D Measurement of Deformation Microstructure of Al(0.2%)Mg Using Submicron Resolution White X-Ray Microbeams
- We have used submicron-resolution white x-ray microbeams on the MHATT-CAT beamline 7-ID at the Advanced Photon Source to develop techniques for three-dimensional investigation of the deformation microstructure in a 20% plane strain compressed Al(0.2%)Mg tri-crystal. Kirkpatrick-Baez mirrors were used to focus white radiation from an undulator to a 0.7 x 0.7 {micro}m{sup 2} beam that was scanned over bi- and tri-crystal regions near the triple-junction of the tri-crystal. Depth resolution along the x-ray microbeam of less than 5 microns was achieved by triangulation to the diffractibn source point using images taken at a series of CCD distances from the microbeam. Computer indexing of the deformation cell structure in the bi-crystal region provided orientations of individual subgrains to {approximately}0.01 degrees, making possible detailed measurements of the rotation axes between individual cells. digital.library.unt.edu/ark:/67531/metadc707283/
- 3-D Micro-Scale Machining Implementation
- The KCP investigated milling, turning, and wire EDMing features at a micro-scale using existing equipment to find the limitations and hindrances to producing future designs. digital.library.unt.edu/ark:/67531/metadc932577/
- 3-D Microprobe Metrology
- This report documents the results of a project undertaken to develop an ultra-high-accuracy measurement capability, which is necessary to address a rising trend toward miniaturized mechanical products exhibiting dramatically reduced product tolerances. A significant improvement in measurement capability is therefore required to insure that a 4:1 ratio can be maintained between product tolerances and measurement uncertainty. digital.library.unt.edu/ark:/67531/metadc926346/
- 3-D Model of Broadband Emission from Supernova Remnants Undergoing Non-linear Diffusive Shock Acceleration
- We present a 3-dimensional model of supernova remnants (SNRs) where the hydrodynamical evolution of the remnant is modeled consistently with nonlinear diffusive shock acceleration occurring at the outer blast wave. The model includes particle escape and diffusion outside of the forward shock, and particle interactions with arbitrary distributions of external ambient material, such as molecular clouds. We include synchrotron emission and cooling, bremsstrahlung radiation, neutral pion production, inverse-Compton (IC), and Coulomb energy-loss. Boardband spectra have been calculated for typical parameters including dense regions of gas external to a 1000 year old SNR. In this paper, we describe the details of our model but do not attempt a detailed fit to any specific remnant. We also do not include magnetic field amplification (MFA), even though this effect may be important in some young remnants. In this first presentation of the model we don't attempt a detailed fit to any specific remnant. Our aim is to develop a flexible platform, which can be generalized to include effects such as MFA, and which can be easily adapted to various SNR environments, including Type Ia SNRs, which explode in a constant density medium, and Type II SNRs, which explode in a pre-supernova wind. When applied to a specific SNR, our model will predict cosmic-ray spectra and multi-wavelength morphology in projected images for instruments with varying spatial and spectral resolutions. We show examples of these spectra and images and emphasize the importance of measurements in the hard X-ray, GeV, and TeV gamma-ray bands for investigating key ingredients in the acceleration mechanism, and for deducing whether or not TeV emission is produced by IC from electrons or pion-decay from protons. digital.library.unt.edu/ark:/67531/metadc895492/
- A 3-D model of superfluid helium suitable for numerical analysis
- The two-fluid description is a very successful phenomenological representation of the properties of Helium II. A 3-D model suitable for numerical analysis based on the Landau-Khalatnikov description of Helium II is proposed. In this paper we introduce a system of partial differential equations that is both complete and consistent as well as practical, to be used for a 3-D solution of the flow of Helium II. The development of a 3-D numerical model for Helium II is motivated by the need to validate experimental results obtained by observing the normal component velocity distribution in a Helium II thermal counter-flow using the Particle Image Velocimetry (PIV) technique. digital.library.unt.edu/ark:/67531/metadc840706/
- A 3-d modular gripper design tool
- Modular fixturing kits are precisely machined sets of components used for flexible, short-turnaround construction of fixtures for a variety of manufacturing purposes. A modular vise is a parallel-jaw vise, where each jaw is a modular fixture plate with a regular grid of precisely positioned holes. A modular vise can be used to locate and hold parts for machining, assembly, and inspection tasks. To fixture a part, one places pins in some of the holes so that when the vise is closed, the part is reliably located and completely constrained. The modular vise concept can be adapted easily to the design of modular parallel-jaw grippers for robots. By attaching a grid plate to each jaw of a parallel-jaw gripper, the authors gain the ability to easily construct high-quality grasps for a wide variety of parts from a standard set of hardware. Wallack and Canny developed a previous algorithm for planning planar grasp configurations for the modular vise. In this paper, the authors expand this work to produce a 3-d fixture/gripper design tool. They describe several analyses added to the planar algorithm to improve its utility, including a three-dimensional grasp quality metric based on geometric and force information, three-dimensional geometric loading analysis, and inter-gripper interference analysis to determine the compatibility of multiple grasps for handing the part from one gripper to another. Finally, the authors describe two applications which combine the utility of modular vise-style grasping with inter-gripper interference: The first is the design of a flexible part-handling subsystem for a part cleaning workcell under development at Sandia National Laboratories; the second is the automatic design of grippers that support the assembly of multiple products on a single assembly line. digital.library.unt.edu/ark:/67531/metadc675560/
- A 3-d modular gripper design tool
- Modular fixturing kits are sets of components used for flexible, rapid construction of fixtures. A modular vise is a parallel-jaw vise, each jaw of which is a modular fixture plate with a regular grid of precisely positioned holes. To fixture a part, one places pins in some of the holes so that when the vise is closed, the part is reliably located and completely constrained. The modular vise concept can be adapted easily to the design of modular parallel-jaw grippers for robots. By attaching a grid-plate to each jaw of a parallel-jaw gripper, one gains the ability to easily construct high-quality grasps for a wide variety of parts from a standard set of hardware. Wallack and Canny developed an algorithm for planning planar grasp configurations for the modular vise. In this paper, the authors expand this work to produce a 3-d fixture/gripper design tool. They describe several analyses they have added to the planar algorithm, including a 3-d grasp quality metric based on force information, 3-d geometric loading analysis, and inter-gripper interference analysis. Finally, the authors describe two applications of their code. One of these is an internal application at Sandia, while the other shows a potential use of the code for designing part of an agile assembly line. digital.library.unt.edu/ark:/67531/metadc674520/
- 3-D Numerical Modeling of a Complex Salt Structure
- Reliably processing, imaging, and interpreting seismic data from areas with complicated structures, such as sub-salt, requires a thorough understanding of elastic as well as acoustic wave propagation. Elastic numerical modeling is an essential tool to develop that understanding. While 2-D elastic modeling is in common use, 3-D elastic modeling has been too computationally intensive to be used routinely. Recent advances in computing hardware, including commodity-based hardware, have substantially reduced computing costs. These advances are making 3-D elastic numerical modeling more feasible. A series of example 3-D elastic calculations were performed using a complicated structure, the SEG/EAGE salt structure. The synthetic traces show that the effects of shear wave propagation can be important for imaging and interpretation of images, and also for AVO and other applications that rely on trace amplitudes. Additional calculations are needed to better identify and understand the complex wave propagation effects produced in complicated structures, such as the SEG/EAGE salt structure. digital.library.unt.edu/ark:/67531/metadc742372/
- A 3-D numerical study of pinhole diffraction to predict the accuracy of EUV point diffraction interferometry
- A 3-D electromagnetic field simulation is used to model the propagation of extreme ultraviolet (EUV), 13-nm, light through sub-1500 {Angstrom} dia pinholes in a highly absorptive medium. Deviations of the diffracted wavefront phase from an ideal sphere are studied within 0.1 numerical aperture, to predict the accuracy of EUV point diffraction interferometersused in at-wavelength testing of nearly diffraction-limited EUV optical systems. Aberration magnitudes are studied for various 3-D pinhole models, including cylindrical and conical pinhole bores. digital.library.unt.edu/ark:/67531/metadc687713/
- A 3-D ray-trace model for an AMR code on distributed processors
- Distributed memory AMR codes provide a special challenge for laser ray-trace modeling. For RAGE we follow the energy-carrying rays through each cell, checking for crossings which require collection a new cell index (1 of 9 in 2D). Density gradients for ray deflections can be calculated 'on the fly.' Substantial excursions must be made from the legacy PIC particle area-weighting approach, but this serves as a useful 1st step. Similarly, IDL now offers a quick graphical rendering, while ENSIGHT graphics beautifully captures the 3D light refraction and deposition. digital.library.unt.edu/ark:/67531/metadc927626/
- 3-D Ray-tracing and 2-D Fokker-Planck Simulations of Radiofrequency Application to Tokamak Plasmas
- A state of the art numerical tool has been developed to simulate the propagation and the absorption of coexisting different types of waves in a tokamak geometry. The code includes a numerical solution of the three-dimensional (R, Z, {Phi}) toroidal wave equation for the electric field of the different waves in the WKBJ approximation. At each step of integration, the two-dimensional (v{sub {parallel}}, v{sub {perpendicular}}) Fokker-Planck equation is solved in the presence of quasilinear diffusion coefficients. The electron Landau damping of the waves is modeled taking into account the interaction of the wave electric fields with the quasilinearly modified distribution function. Consistently, the code calculates the radial profiles of non-inductively generated current densities, the transmitted power traces and the total power damping curves. Synergistic effects among the different type of waves (e.g., lower hybrid and ion Bernstein waves) are studied through the separation of the contributions of the single wave from the effects due to their coexistence. digital.library.unt.edu/ark:/67531/metadc703436/
- 3-D Reservoir and Stochastic Fracture Network Modeling for Enhanced Oil Recovery, Circle Ridge Phosphoria/Tensleep Reservoir, and River Reservation, Arapaho and Shoshone Tribes, Wyoming
- The goal of this project is to improve the recovery of oil from the Circle Ridge Oilfield, located on the Wind River Reservation in Wyoming, through an innovative integration of matrix characterization, structural reconstruction, and the characterization of the fracturing in the reservoir through the use of discrete fracture network models. digital.library.unt.edu/ark:/67531/metadc743258/
- 3-D Reservoir and Stochastic Fracture Network Modeling for Enhanced Oil Recovery, Circle Ridge Phosphoria/Tensleep Reservoir, and River Reservation, Arapaho and Shoshone Tribes, Wyoming
- The goal of this project is to improve the recovery of oil from the Circle Ridge Oilfield, located on the Wind River Reservation in Wyoming, through an innovative integration of matrix characterization, structural reconstruction, and the characterization of the fracturing in the reservoir through the use of discrete fracture network models. digital.library.unt.edu/ark:/67531/metadc742517/
- 3-D Reservoir and Stochastic Fracture Network Modeling for Enhanced Oil Recovery, Circle Ridge Phosphoria/Tensleep Reservoir, and River Reservation, Arapaho and Shoshone Tribes, Wyoming
- The goal of this project is to improve the recovery of oil from the Circle Ridge Oilfield, located on the Wind River Reservation in Wyoming, through an innovative integration of matrix characterization, structural reconstruction, and the characterization of the fracturing in the reservoir through the use of discrete fracture network models. digital.library.unt.edu/ark:/67531/metadc735678/
- The goal of this project is to improve the recovery of oil from the Circle Ridge Oilfield, located on the Wind River Reservation in Wyoming, through an innovative integration of matrix characterization, structural reconstruction, and the characterization of the fracturing in the reservoir through the use of discrete fracture network models. digital.library.unt.edu/ark:/67531/metadc737473/
- 3-D RESERVOIR AND STOCHASTIC FRACTURE NETWORK MODELING FOR ENHANCED OIL RECOVERY, CIRCLE RIDGE PHOSPHORIA/TENSLEEP RESERVOIR, WIND RIVER RESERVATION, ARAPAHO AND SHOSHONE TRIBES, WYOMING
- This report describes the results made in fulfillment of contract DE-FG26-00BC15190, ''3-D Reservoir and Stochastic Fracture Network Modeling for Enhanced Oil Recovery, Circle Ridge Phosphoria/Tensleep Reservoir, Wind River Reservation, Arapaho and Shoshone Tribes, Wyoming''. The goal of this project is to improve the recovery of oil from the Tensleep and Phosphoria Formations in Circle Ridge Oilfield, located on the Wind River Reservation in Wyoming, through an innovative integration of matrix characterization, structural reconstruction, and the characterization of the fracturing in the reservoir through the use of discrete fracture network models. Fields in which natural fractures dominate reservoir permeability, such as the Circle Ridge Field, often experience sub-optimal recovery when recovery processes are designed and implemented that do not take advantage of the fracture systems. For example, a conventional waterflood in a main structural block of the Field was implemented and later suspended due to unattractive results. It is estimated that somewhere less than 20% of the OOIP in the Circle Ridge Field have been recovered after more than 50 years' production. Marathon Oil Company identified the Circle Ridge Field as an attractive candidate for several advanced IOR processes that explicitly take advantage of the natural fracture system. These processes require knowledge of the distribution of matrix porosity, permeability and oil saturations; and understanding of where fracturing is likely to be well-developed or poorly developed; how the fracturing may compartmentalize the reservoir; and how smaller, relatively untested subthrust fault blocks may be connected to the main overthrust block. For this reason, the project focused on improving knowledge of the matrix properties, the fault block architecture and to develop a model that could be used to predict fracture intensity, orientation and fluid flow/connectivity properties. Knowledge of matrix properties was greatly extended by calibrating wireline logs from 113 wells with incomplete or older-vintage logging suites to wells with a full suite of modern logs. The model for the fault block architecture was derived by 3D palinspastic reconstruction. This involved field work to construct three new cross-sections at key areas in the Field; creation of horizon and fault surface maps from well penetrations and tops; and numerical modeling to derive the geometry, chronology, fault movement and folding history of the Field through a 3D restoration of the reservoir units to their original undeformed state. The methodology for predicting fracture intensity and orientation variations throughout the Field was accomplished by gathering outcrop and subsurface image log fracture data, and comparing it to the strain field produced by the various folding and faulting events determined through the 3D palinspastic reconstruction. It was found that the strains produced during the initial folding of the Tensleep and Phosphoria Formations corresponded well without both the orientations and relative fracture intensity measured in outcrop and in the subsurface. The results have led to a 15% to 20% increase in estimated matrix pore volume, and to the plan to drill two horizontal drain holes located and oriented based on the modeling results. Marathon Oil is also evaluating alternative tertiary recovery processes based on the quantitative 3D integrated reservoir model. digital.library.unt.edu/ark:/67531/metadc738026/
- A 3-D SAR approach to IFSAR processing
- Interferometric SAR (IFSAR) can be shown to be a special case of 3-D SAR image formation. In fact, traditional IFSAR processing results in the equivalent of merely a super-resolved, under-sampled, 3-D SAR image. However, when approached as a 3-D SAR problem, a number of IFSAR properties and anomalies are easily explained. For example, IFSAR decorrelation with height is merely ordinary migration in 3-D SAR. Consequently, treating IFSAR as a 3-D SAR problem allows insight and development of proper motion compensation techniques and image formation operations to facilitate optimal height estimation. Furthermore, multiple antenna phase centers and baselines are easily incorporated into this formulation, providing essentially a sparse array in the elevation dimension. This paper shows the Polar Format image formation algorithm extended to 3 dimensions, and then proceeds to apply it to the IFSAR collection geometry. This suggests a more optimal reordering of the traditional IFSAR processing steps. digital.library.unt.edu/ark:/67531/metadc702258/
- 3-D sedimentological and geophysical studies of clastic reservoir analogs: Facies architecture, reservoir properties, and flow behavior within delta front facies elements of the Cretaceous Wall Creek Member, Frontier Formation, Wyoming
- This project examined the internal architecture of delta front sandstones at two locations within the Turonian-age Wall Creek Member of the Frontier Formation, in Wyoming. The project involved traditional outcrop field work integrated with core-data, and 2D and 3D ground penetrating radar (GPR) imaging from behind the outcrops. The fluid-flow engineering work, handled through a collaborative grant given to PI Chris White at LSU, focused on effects on fluid flow of late-stage calcite cement nodules in 3D. In addition to the extensive field component, the work funded 2 PhD students (Gani and Lee) and resulted in publication of 10 technical papers, 17 abstracts, and 4 internal field guides. PI Bhattacharya also funded an additional 3 PhD students that worked on the Wall Creek sandstone funded separately through an industrial consortium, two of whom graduated in the fall 2006 ((Sadeque and Vakarelov). These additional funds provided significant leverage to expand the work to include a regional stratigraphic synthesis of the Wall Creek Member of the Frontier Formation, in addition to the reservoir-scale studies that DOE directly funded. Awards given to PI Bhattacharya included the prestigious AAPG Distinguished Lecture Award, which involved a tour of about 25 Universities and Geological Societies in the US and Canada in the fall of 2005 and Spring of 2006. Bhattacharya gave two talks, one entitled “Applying Deltaic and Shallow Marine Outcrop Analogs to the Subsurface”, which highlighted the DOE sponsored work and the other titled “Martian River Deltas and the Origin of Life”. The outcrop analog talk was given at about 1/2 of the venues visited. digital.library.unt.edu/ark:/67531/metadc887167/
- 3-D sedimentological and geophysical studies of clastic reservoir analogs: Facies architecture, reservoir properties, and flow behavior within delta front facies elements of the Cretaceous Wall Creek Member, Frontier Formation, Wyoming
- No Description digital.library.unt.edu/ark:/67531/metadc886580/
- 3-D Seismic Experimentation and Advanced Processing/Inversion Development for Investigations of the Shallow Subsurface
- Under ER63662, 3-D Seismic Experimentation and Advanced Processing/Inversion Development for Investigations of the Shallow Subsurface, we have completed a number of subprojects associated with the Hill Air Force Base (HAFB) high resolution 3-D reflection/tomography dataset. digital.library.unt.edu/ark:/67531/metadc786536/
- 3-D Seismic Experimentation and Advanced Processing/Inversion Development for Investigations of the Shallow Subsurface
- Gian Fradelizio, a Rice Ph.D. student has completed reprocessing the 3D seismic reflection data acquired at Hill AFB through post-stack depth migration for comparison to the traveltime and waveform tomography results. Zelt, Levander, Fradelizio, and 5 others spent a week at Hill AFB in September 2005, acquiring an elastic wave data set along 2 profiles. We used 60 3-component Galperin mounted 40 Hz geophones recorded by 3 GEOMETRICS Stratavision systems. The seismic source employed was a sledgehammer used to generate transverse, and radial, and vertical point source data. Data processing has begun at Rice to generate S-wave reflection and refraction images. We also acquired surface wave and ground penetrating rada data to complement the elastic wave dataset. digital.library.unt.edu/ark:/67531/metadc884379/
- 3-D Seismic Exploration Project, Ute Indian Tribe, Uintah and Ouray Reservation, Uintah County, Utah
- The objectives of this North Hill Creek 3-D seismic survey were to: (1) cover as large an area as possible with available budget; (2) obtain high quality data throughout the depth range of the prospective geologic formations of 2,000' to 12,000' to image both gross structures and more subtle structural and stratigraphic elements; (3) overcome the challenges posed by a hard, reflective sandstone that cropped out or was buried just a few feet below the surface under most of the survey area; and (4) run a safe survey. digital.library.unt.edu/ark:/67531/metadc742671/
- 3-D Seismic Methods for Geothermal Reservoir Exploration and Assessment--Summary
- A wide variety of seismic methods covering the spectrum from DC to kilohertz have been employed at one time or the other in geothermal environments. The reasons have varied from exploration for a heat source to attempting to find individual fractures producing hot fluids. For the purposes here we will assume that overall objective of seismic imaging is for siting wells for successful location of permeable pathways (often fracture permeability) that are controlling flow and transport in naturally fractured reservoirs. The application could be for exploration of new resources or for in-fill/step-out drilling in existing fields. In most geothermal environments the challenge has been to separate the ''background'' natural complexity and heterogeneity of the matrix from the fracture/fault heterogeneity controlling the fluid flow. Ideally one not only wants to find the fractures, but the fractures that are controlling the flow of the fluids. Evaluated in this work is current state-of-the-art surface (seismic reflection) and borehole seismic methods (Vertical Seismic Profiling (VSP), Crosswell and Single Well) to locate and quantify geothermal reservoir characteristics. The focus is on active methods; the assumption being that accuracy is needed for successful well siting. Passive methods are useful for exploration and detailed monitoring for in-fill drilling, but in general the passive methods lack the precision and accuracy for well siting in new or step out areas. In addition, MEQ activity is usually associated with production, after the field has been taken to a mature state, thus in most cases it is assumed that there is not enough MEQ activity in unproduced areas to accurately find the permeable pathways. The premise of this review is that there may new developments in theory and modeling, as well as in data acquisition and processing, which could make it possible to image the subsurface in much more detail than 15 years ago. New understanding of the effect of fractures on seismic wave propagation are now being applied to image fractures in gas and oil environments. It now may be appropriate to apply these methods, with modifications, to geothermal applications. It is assumed that to implement the appropriate methods an industry coupled program tightly linked to actual field cases, iterating between development and application will be pursued. The goal of this work is to evaluate the most promising methods and approaches that may be used for improved geothermal exploration and reservoir assessment. It is not a comprehensive review of all seismic methods used to date in geothermal environments. This work was motivated by a need to assess current and developing seismic technology that if applied in geothermal cases may greatly improve the chances for locating new geothermal resources and/or improve assessment of current ones. digital.library.unt.edu/ark:/67531/metadc787122/
- 3-D Silicon Photonic Lattices- Cornerstone of an Emerging Photonics Revolution
- Three-dimensional photonic lattices are engineered materials which are the photonic analogues of semiconductors. These structures were first proposed and demonstrated in the mid-to-late 1980's. However, due to fabrication difficulties, lattices active in the infrared are only just emerging. Wide ranges of structures and fabrication approaches have been investigated. The most promising approach for many potential applications is a diamond-like structure fabricated using silicon microprocessing techniques. This approach has enabled the fabrication of 3-D silicon photonic lattices active in the infrared. The structures display band gaps centered from 12{micro} down to 1.55{micro}. digital.library.unt.edu/ark:/67531/metadc793231/
- 3-D SIMULATION FOR ASSESSMENT OF TRANSPARENT WEAPON DISASSEMBLY OPERATIONS
- No abstract prepared. digital.library.unt.edu/ark:/67531/metadc717008/
- 3-D Simulations of Plasma Wakefield Acceleration with Non-Idealized Plasmas and Beams
- 3-D Particle-in-cell OSIRIS simulations of the current E-162 Plasma Wakefield Accelerator Experiment are presented in which a number of non-ideal conditions are modeled simultaneously. These include tilts on the beam in both planes, asymmetric beam emittance, beam energy spread and plasma inhomogeneities both longitudinally and transverse to the beam axis. The relative importance of the non-ideal conditions is discussed and a worst case estimate of the effect of these on energy gain is obtained. The simulation output is then propagated through the downstream optics, drift spaces and apertures leading to the experimental diagnostics to provide insight into the differences between actual beam conditions and what is measured. The work represents a milestone in the level of detail of simulation comparisons to plasma experiments. digital.library.unt.edu/ark:/67531/metadc877163/
- 3-D Spectral Induced Polarization (IP) Imaging: Non-Invasive Characterization Of Contaminant Plumes
- The overall objective of this project is to develop the scientific basis for characterizing contaminant plumes in the earth's subsurface using field measurements of induced polarization (IP) effects. Three specific objectives towards this end are 1. 2. 3. Understanding IP at the laboratory level through measurements of complex resistivity as a function of frequency in rock and soil samples with varying pore geometries, pore fluid conductivities and saturations, and contaminant chemistries and concentrations. Developing effective data acquisition techniques for measuring the critical IP responses (time domain or frequency domain) in the field. Developing modeling and inversion algorithms that permit the interpretation of field IP data in terms of subsurface geology and contaminant plume properties. digital.library.unt.edu/ark:/67531/metadc785886/
- 3-D Spectral IP Imaging: Non-Invasive Characterization DE FG02 96ER 14714
- The Earth Resources Laboratory (ERL) performed a broad foundational study of spectral induced polarization (SIP) for site characterization. The project encompassed laboratory studies of microgeometry and chemistry effects on Induced Polarization (IP), an investigation of electromagnetic coupling (emc) noise, and development of 3D modeling and inversion codes. The major finding of the project is that emc noise presents a critical limitation for field implementation of SIP and conventional correction methods are inadequate. The project developed a frequency domain 3D complex resistivity modeling and inversion code Laboratory experiments were conducted to study the effects of solution chemistry and microgeometry on the SIP response of sandstone. Results indicate that changes in chemistry affect the magnitude of the spectral IP response and changes in microgeometry affect the shape of the spectral IP response. The developed physiochemical IP model can be used to invert spectral IP data for an apparent grain size distribution. Laboratory studies over the last twenty years have shown that SIP data must be acquired over several decades of frequency and include frequencies greater than 1kHz. A model of the components of emc noise has been developed and investigation with this model showed that inductive coupling is the most significant component. The study concluded that emc limits the frequency range of usable field data to approximately 100 Hz and below for typical site conditions. Several correction schemes have been developed based on treating emc as noise to be removed from the data, but our investigation has shown that these are not adequate for high frequencies, greater than 100Hz. Laboratory studies have demonstrated that the greatest response is the frequency range greater than 1KHz, hence the emc problem must be resolved for field implementation of SIP to advance. The ERL developed 2D/3D time domain codes that perform inversions for charge abilities based on schemes introduced by Siegel (1959). The ERL has also developed a 3D complex resistivity code for inversion of frequency domain IP data. The algorithm accommodates a general earth model with a complex electrical resistivity as a function of frequency and 3-D spatial position. The forward problem is solved by the complex biconjugate gradient method, while the regularized inverse problem is solved by the nonlinear conjugate gradient method. Time domain field data was acquired along a single survey line at the FS-12 plume, Massachusetts Military Reservation. The data was inverted with the developed 2D time domain code with the results having an excellent match to monitoring well data. The data was further analyzed by an innovative scheme where a ''gross spectral chargeability'' was determined from time domain data. The results show that there is valuable information in the spectra of the data. digital.library.unt.edu/ark:/67531/metadc788396/
- 3-D Spectral IP Imaging: Non-Invasive Characterization of Contaminant Plumes
- The overall objective of this project is to develop the scientific basis for characterizing contaminant plumes in the earth's subsurface using field measurements of induced polarization (IP) effects. Three specific objectives towards this end are: (1) Understanding IP at the laboratory level through measurements of complex resistivity as a function of frequency in rock and soil samples with varying pore geometries, pore fluid conductivities and saturations, and contaminant chemistries and concentrations. (2) Developing effective data acquisition techniques for measuring the critical IP responses (time domain or frequency domain) in the field. (3) Developing modeling and inversion algorithms that permit the interpretation of field IP data in terms of subsurface geology and contaminant plume properties. digital.library.unt.edu/ark:/67531/metadc782451/
- 3-D spectral IP imaging: Non-invasive characterization of contaminant plumes. 1998 annual progress report
- 'The overall objective of this project is to develop the scientific basis for characterizing contaminant plumes in the earth''s subsurface using field measurements of induced polarization (IP) effects. Three specific objectives towards this end are: (1) understanding IP at the laboratory level through measurements of complex resistivity as a function of frequency in rock and soil samples with varying pore geometries, pore fluid conductivities and saturations, and contaminant chemistries and concentrations; (2) developing effective data acquisition techniques for measuring the critical IP responses (time domain or frequency domain) in the field; (3) developing modeling and inversion algorithms that permit the interpretation of field IP data in terms of subsurface geology and contaminant plume properties. The authors laboratory experiments to date are described in Appendices A and B, which consist of two papers submitted to the annual SAGEEP conference (Frye et al., 1998; Sturrock et al., 1998). The experiments involved measurements of complex resistivity vs. frequency on a suite of brine saturated sandstone samples. In one set of experiments, the fluid chemistry (pH, ionic strength, and cation type) was varied. In a second set of experiments, the microgeometry of the rock matrix was varied. The experiments showed that spectral IP responses are sensitive to subtle variations in both the solution chemistry and rock microgeometry. The results demonstrate that spectral IP responses have the potential of being sensitive indicators of in-situ chemistry and microgeometry, the latter of which may be related to the hydraulic properties. Data Acquisition The authors have been looking in some detail at the effects of electromagnetic coupling and how to practically deal with it. In this area, the results to date are summarized in Vandiver (1998). The progress in the development of modeling and inversion algorithms for IP is described in Appendix C, a paper submitted to the annual SAGEEP conference (Shi et al., 1998). The authors have developed algorithms for forward modeling and inversion of spectral IP data in 3-D media. The algorithms accommodate a general earth model with a complex electrical conductivity as a function of frequency and 3-D spatial position. Using regularization and optimization techniques, the inversion algorithm obtains a 3-D image of resistivity amplitude and phase for each frequency contained in the data set. They have begun testing their algorithms on synthetic data generated from a simple model of a contaminant plume. The complex resistivity parameters of the background medium and plume are based on the laboratory results described above.' digital.library.unt.edu/ark:/67531/metadc625746/
- 3-D spectral IP imaging: Non-invasive characterization of contaminant plumes. Annual progress report, September 15, 1996--September 14, 1997
- 'The objective of this project is to develop the scientific basis for characterizing contaminant plumes in the earth''s subsurface using field measurements of induced polarization (IP) effects. The first-year accomplishments are (1) laboratory experiments on fluid-saturated sandstones quantifying the dependence of spectral IP responses on solution chemistry and rock micro-geometry; (2) library research on the current understanding of electromagnetic coupling effects on IP data acquired in the field: and (3) development of prototype forward modeling and inversion algorithms for interpreting IP data in terms of 3-D models of complex resistivity.' digital.library.unt.edu/ark:/67531/metadc620009/
- 3-D surface profile measurements of large x-ray synchrotron radiation mirrors using stitching interferometry.
- Stitching interferometry, using small-aperture, high-resolution, phase-measuring interferometry, has been proposed for quite some time now as a metrology technique to obtain 3-dimensional profiles of surfaces of oversized optical components and substrates. The aim of this work is to apply this method to the specific case of long grazing-incidence x-ray mirrors, such as those used in beamlines at synchrotron radiation facilities around the world. Both fabrication and characterization of these mirrors would greatly benefit from this technique because it offers the potential for providing measurements with accuracy and resolution better than those obtained using existing noncontact laser profilers, such as the long trace profiler (LTP). Measurement data can be used as feedback for computer-controlled fabrication processes to correct for possible topography errors. The data can also be used for simulating and predicting mirror performance under realistic conditions. A semiautomated stitching system was built and tested at the X-ray Optics Metrology Laboratory of the Advanced Photon Source at Argonne National Laboratory. The initial objective was to achieve a measurement sensitivity on the order of 1 {micro}rad rms. Preliminary tests on a 1 m-long x-ray mirror showed system repeatability of less than 0.6 {micro}rad rms. This value is comparable to that of a conventional LTP. The measurement accuracy was mostly affected by environmental perturbations and system calibration effects. With a fully automated and improved system (to be built in the near future), we expect to achieve measurement sensitivity on the order of 0.0 {micro}rad rms or better. In this paper, after a brief review of basic principles and general technical difficulties and challenges of the stitching technique, a detailed description of the measurement setup is given and preliminary results obtained with it are analyzed and discussed. digital.library.unt.edu/ark:/67531/metadc741012/
- 3-D Target Location from Stereoscopic SAR Images
- SAR range-Doppler images are inherently 2-dimensional. Targets with a height offset lay over onto offset range and azimuth locations. Just which image locations are laid upon depends on the imaging geometry, including depression angle, squint angle, and target bearing. This is the well known layover phenomenon. Images formed with different aperture geometries will exhibit different layover characteristics. These differences can be exploited to ascertain target height information, in a stereoscopic manner. Depending on the imaging geometries, height accuracy can be on the order of horizontal position accuracies, thereby rivaling the best IFSAR capabilities in fine resolution SAR images. All that is required for this to work are two distinct passes with suitably different geometries from any plain old SAR. digital.library.unt.edu/ark:/67531/metadc619742/
- 3-D THERMAL EVALUATIONS FOR a FUELED EXPERIMENT in the ADVANCED TEST REACTOR
- The DOE Advanced Fuel Cycle Initiative and Generation IV reactor programs are developing new fuel types for use in the current Light Water Reactors and future advanced reactor concepts. The Advanced Gas Reactor program is planning to test fuel to be used in the Next Generation Nuclear Plant (NGNP) nuclear reactor. Preliminary information for assessing performance of the fuel will be obtained from irradiations performed in the Advanced Test Reactor large ''B'' experimental facility. A test configuration has been identified for demonstrating fuel types typical of gas cooled reactors or fast reactors that may play a role in closing the fuel cycle or increasing efficiency via high temperature operation Plans are to have 6 capsules, each containing 12 compacts, for the test configuration. Each capsule will have its own temperature control system. Passing a helium-neon gas through the void regions between the fuel compacts and the graphite carrier and between the graphite carrier and the capsule wall will control temperature. This design with three compacts per axial level was evaluated for thermal performance to ascertain the temperature distributions in the capsule and test specimens with heating rates that encompass the range of initial heat generation rates. digital.library.unt.edu/ark:/67531/metadc787091/
- 3-D Thermal Evaluations for a Fueled Experiment in the Advanced Test Reactor
- The DOE Advanced Fuel Cycle Initiative and Generation IV reactor programs are developing new fuel types for use in the current Light Water Reactors and future advanced reactor concepts. The Advanced Gas Reactor program is planning to test fuel to be used in the Next Generation Nuclear Plant (NGNP) nuclear reactor. Preliminary information for assessing performance of the fuel will be obtained from irradiations performed in the Advanced Test Reactor large “B” experimental facility. A test configurations has been identified for demonstrating fuel types typical of gas cooled reactors or fast reactors that may play a role in closing the fuel cycle or increasing efficiency via high temperature operation Plans are to have 6 capsules, each containing 12 compacts, for the test configuration. Each capsule will have its own temperature control system. Passing a helium-neon gas through the void regions between the fuel compacts and the graphite carrier and between the graphite carrier and the capsule wall will control temperature. This design with three compacts per axial level was evaluated for thermal performance to ascertain the temperature distributions in the capsule and test specimens with heating rates that encompass the range of initial heat generation rates. digital.library.unt.edu/ark:/67531/metadc881353/
- 3-D UNSTRUCTURED HEXAHEDRAL-MESH Sn TRANSPORT METHODS
- This is the final report of a three-year, Laboratory-Directed Research and Development (LDRD) project at the Los Alamos National Laboratory (LANL). We have developed a method for solving the neutral-particle transport equation on 3-D unstructured hexahedral meshes using a S{sub n} discretization in angle in conjunction with a discontinuous finite-element discretization in space and a multigroup discretization in energy. Previous methods for solving this equation in 3-D have been limited to rectangular meshes. The unstructured-mesh method that we have developed is far more efficient for solving problems with complex 3-D geometric features than rectangular-mesh methods. In spite of having to make several compromises in our spatial discretization technique and our iterative solution technique, our method has been found to be both accurate and efficient for a broad class of problems. digital.library.unt.edu/ark:/67531/metadc715757/
- 3-D vertical seismic profiling at LLNL Site 300
- The initial goal of the 3-D Vertical Seismic Profiling (VSP) work at LLNL was to characterize seismic wave velocities and frequencies below the vadose zone to design the acquisition geometry for a 3-D shallow surface seismic reflection survey. VSPs are also used routinely to provide a link between surface seismic data and well logs. However, a test 2-D seismic line recorded at LLNL in the Spring of 1994 indicated that obtaining high quality reflection images below the vadose zone, yet shallower that 50 m, would require an expensive, very finely sampled survey ({lt} 1 m receiver spacing). Extensive image processing of the LLNL 2-D test line indicated that the only reliable reflection was from the top of the water table. Surprisingly, these results were very different than recent 3-D seismic work recorded at other sites, where high quality, high frequency surface (up to 300 Hz) reflection images were obtained as shallow as 20m. We believe that the differences are primarily due to the comparatively deep vadose zone at LLNL (15 to 30m) as compared to 0-5m at other sites. The thick vadose zone attenuates the reflection signals, particularly at the high frequencies (above 100 @). In addition, the vadose zone at LLNL creates a seismogram in which surface-propagating noise overlaps with the reflection signals for reflections above 50 m. By contrast, when the vadose zone is not thick, high frequencies can propagate and noise will not overlap with reflections as severely. Based on the results from the 2-D seismic line and the encouraging results from a VSP run concurrent with the 2-D seismic experiment, we modified the objectives of the research and expanded the scope of the VSP imaging at LLNL. We conducted two 3-D multi-offset VSP experiments at LLNL in the Summer and Fall of 1994. These VSP experiments were designed to characterize the seismic propagation characteristics at two different locations on the LLNL site: the first was a well with a relatively shallow water table (10m), while the second was a well with a relatively deep water table (about 25m). Other goals of the VSP experiments were to characterize the velocity structure in the vicinity of boreholes, and to attempt to image reflections away from the boreholes. We found that the usable frequency content appeared to vary with water table level. digital.library.unt.edu/ark:/67531/metadc693116/
- A 3-D Vortex Code for Parachute Flow Predictions: VIPAR Version 1.0
- This report describes a 3-D fluid mechanics code for predicting flow past bluff bodies whose surfaces can be assumed to be made up of shell elements that are simply connected. Version 1.0 of the VIPAR code (Vortex Inflation PARachute code) is described herein. This version contains several first order algorithms that we are in the process of replacing with higher order ones. These enhancements will appear in the next version of VIPAR. The present code contains a motion generator that can be used to produce a large class of rigid body motions. The present code has also been fully coupled to a structural dynamics code in which the geometry undergoes large time dependent deformations. Initial surface geometry is generated from triangular shell elements using a code such as Patran and is written into an ExodusII database file for subsequent input into VIPAR. Surface and wake variable information is output into two ExodusII files that can be post processed and viewed using software such as EnSight{trademark}. digital.library.unt.edu/ark:/67531/metadc740481/
- 3-D woven, mullite matrix, composite filter
- Westinghouse, with Techniweave as a major subcontractor, is conducting a three-phase program aimed at providing advanced candle filters for a 1996 pilot scale demonstration in one of the two hot gas filter systems at Southern Company Service`s Wilsonville PSD Facility. The Base Program (Phases I and II) objective is to develop and demonstrate the suitability of the Westinghouse/Techniweave next generation composite candle filter for use in Pressurized Fluidized Bed Combustion (PFBC) and/or Integrated Gasification Combined Cycle (IGCC) power generation systems. The Optional Task (Phase M, Task 5) objective is to fabricate, inspect and ship to Wilsonville Hot gas particulate filters are key components for the successful commercializaion of advanced coal-based power-generation systems such as Pressurized Fluidized-bed Combustion (PFBC), including second-generation PFBC, and Integrated Gasification Combined Cycles (IGCC). Current generation monolithic ceramic filters are subject to catastrophic failure because they have very low resistance to crack propagation. To overcome this problem, a damage-tolerant ceramic filter element is needed. digital.library.unt.edu/ark:/67531/metadc623027/
- A 3-Dimensional discrete fracture network generator to examine fracture-matrix interaction using TOUGH2
- 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. digital.library.unt.edu/ark:/67531/metadc739350/