You limited your search to:

  Partner: UNT Libraries Government Documents Department
2-D and 3-D Elastic Modeling with Shared Seismic Models
Several elastic models, both 2-D and 3-D, are being built for use in calculating synthetic elastic seismic data. The models will be made available to the research community, along with the synthetic data that are being calculated from them. These shared models have been proposed or contributed by participants in a collaborative industry, national laboratory, and university research project. The purpose of the modeling is to provide synthetic data to better understand elastic wave propagation and the effects of structural and stratigraphic complexities. The 2-D models are easier to design and change and synthetic calculations can be run relatively quickly in them. It will be possible to alter their layer properties and calculate time-lapse data sets from them. Field data will be available to accompany many of the 2-D models. 3-D models are more realistic, but more difficult to design and change. They also require considerably more computing resources to calculate synthetic data from them. A new 3-D model is being designed, and will be used for computing synthetic elastic data.
2-D electric fields and drifts near the magnetic separatrix in divertor tokamaks
A 2-D calculation is presented for the transport of plasma in the edge region of a divertor tokamak solving continuity, momentum, and energy balance fluid equations. The model uses anomalous radial diffusion, including perpendicular ion momentum, and classical cross-field drifts transport. Parallel and perpendicular currents yield a self-consistent electrostatic potential on both sides of the magnetic separatrix. Outside the separatrix, the simulation extends to material divertor plates where the incident plasma is recycled as neutral gas and where the plate sheath and parallel currents dominate the potential structure. Inside the separatrix, various radial current terms - from viscosity, charge-exchange and poloidal damping, inertia, and {triangledown}B - contribute to the determining the potential. The model rigorously enforces cancellation of gyro-viscous and magnetization terms from the transport equations. The results emphasize the importance of E x B particle flow under the X-point which depends on the sign of the toroidal magnetic field. Radial electric field (E{sub y}) profiles at the outer midplane are small with weak shear when high L-mode diffusion coefficients are used and are large with strong shear when smaller H-mode diffusion coefficients are used. The magnitude and shear of the electric field (E{sub y}) is larger both when the core toroidal rotation is co-moving with the inductive plasma current and when the ion {triangledown}B-drift is towards the single-null X-point.
2-D Finite Element Cable and Box IEMP Analysis
A 2-D finite element code has been developed for the solution of arbitrary geometry cable SGEMP and box IEMP problems. The quasi- static electric field equations with radiation- induced charge deposition and radiation-induced conductivity y are numerically solved on a triangular mesh. Multiple regions of different dielectric materials and multiple conductors are permitted.
2-D image segmentation using minimum spanning trees
This paper presents a new algorithm for partitioning a gray-level image into connected homogeneous regions. The novelty of this algorithm lies in the fact that by constructing a minimum spanning tree representation of a gray-level image, it reduces a region partitioning problem to a minimum spanning tree partitioning problem, and hence reduces the computational complexity of the region partitioning problem. The tree-partitioning algorithm, in essence, partitions a minimum spanning tree into subtrees, representing different homogeneous regions, by minimizing the sum of variations of gray levels over all subtrees under the constraints that each subtree should have at least a specified number of nodes, and two adjacent subtrees should have significantly different average gray-levels. Two (faster) heuristic implementations are also given for large-scale region partitioning problems. Test results have shown that the segmentation results are satisfactory and insensitive to noise.
2-D Imaging of Electron Temperature in Tokamak Plasmas
By taking advantage of recent developments in millimeter wave imaging technology, an Electron Cyclotron Emission Imaging (ECEI) instrument, capable of simultaneously measuring 128 channels of localized electron temperature over a 2-D map in the poloidal plane, has been developed for the TEXTOR tokamak. Data from the new instrument, detailing the MHD activity associated with a sawtooth crash, is presented.
2-D linear motion system. Innovative technology summary report
The US Department of Energy's (DOE's) nuclear facility decontamination and decommissioning (D and D) program requires buildings to be decontaminated, decommissioned, and surveyed for radiological contamination in an expeditious and cost-effective manner. Simultaneously, the health and safety of personnel involved in the D and D activities is of primary concern. D and D workers must perform duties high off the ground, requiring the use of manlifts or scaffolding, often, in radiologically or chemically contaminated areas or in areas with limited access. Survey and decontamination instruments that are used are sometimes heavy or awkward to use, particularly when the worker is operating from a manlift or scaffolding. Finding alternative methods of performing such work on manlifts or scaffolding is important. The 2-D Linear Motion System (2-D LMS), also known as the Wall Walker{trademark}, is designed to remotely position tools and instruments on walls for use in such activities as radiation surveys, decontamination, and painting. Traditional (baseline) methods for operating equipment for these tasks require workers to perform duties on elevated platforms, sometimes several meters above the ground surface and near potential sources of contamination. The Wall Walker 2-D LMS significantly improves health and safety conditions by facilitating remote operation of equipment. The Wall Walker 2-D LMS performed well in a demonstration of its precision, accuracy, maneuverability, payload capacity, and ease of use. Thus, this innovative technology is demonstrated to be a viable alternative to standard methods of performing work on large, high walls, especially those that have potential contamination concerns. The Wall Walker was used to perform a final release radiological survey on over 167 m{sup 2} of walls. In this application, surveying using a traditional (baseline) method that employs an aerial lift for manual access was 64% of the total cost of the improved technology. However, for areas over approximately 600 m{sup 2}, the Wall Walker would cost less than the baseline. Using the Wall Walker 2-D LMS, ALARA exposure and worker safety is improved, and there is potential for increased productivity. This innovative technology performed better than the baseline by providing real-time monitoring of the tool or instrument position. Also, the Wall Walker 2-D LMS can traverse any two-dimensional path at constant speeds of up to 18.3 linear meters per minute (60 linear feet per minute). The survey production rate for the innovative technology was about 0.6 m{sup 2}/min (6 ft{sup 2}/min); the baseline production rate was approximately 0.3 m{sup 2}/min (3 ft{sup 2}/min), using the same surveying instrument and maximum scanning rate.
A 2-D Pore-Network Model of the Drying of Single-Component Liquids in Porous Media
The drying of liquid-saturated porous media is typically approaching using macroscopic continuum models involving phenomenological coefficients. Insight on these coefficients can be obtained by a more fundamental study at the pore- and pore-network levels. In this report, a model based on pore-network representation of porous media that accounts for various process at the pore-scale is presented. These include mass transfer by advection and diffusion in the gas phase, viscous flow in liquid and gas phases and capillary effects at the gas-liquid menisci in the pore throats.
A 2-D Self-Consistent DSMC Model for Chemically Reacting Low Pressure Plasma Reactors
This paper will focus on the methodology of using a 2D plasma Direct Simulation Monte Carlo technique to simulate the species transport in an inductively coupled, low pressure, chemically reacting plasma system. The pressure in these systems is typically less than 20 mtorr with plasma densities of approximately 10{sup 17} {number_sign}/m{sup 3} and an ionization level of only 0.1%. This low ionization level tightly couples the neutral, ion, and electron chemistries and interactions in a system where the flow is subsonic. We present our strategy and compare simulation results to experimental data for Cl{sub 2} in a Gaseous Electronics Conference (GEC) reference cell modified with an inductive coil.
2-D Simulations of Drainage Winds and Diffusion Compared to Observations
A vertically integrated dynamical drainage flow model is developed from conservation equations for momentum and mass in a terrain-following coordinate system. Wind fields from the dynamical model drive a Monte Carlo transport and diffusion model. The model needs only topographic data, an Eulerian or Lagrangian time scale and a surface drag coefficient for input data, and can be started with a motionless atmosphere. Model wind and diffusion predictions are compared to observations from the rugged Geysers CA area. Model winds generally agree with observed surface winds, and in some cases may give better estimates of area-averaged flow than point observations. Tracer gas concentration contours agree qualitatively with observed contours, and point predictions of maximum concentrations were correctly predicted to within factors of 2 to 10. Standard statistical tests of model skill showed that the accuracy of the predictions varied significantly from canyon to canyon in the Geysers are a. Model wind predictions are also compared to observations from the Savannah River Plant of SC which has gently rolling terrain. The model correctly simulated the slower development of drainage winds and slower deepening of the drainage layer in the Savannah River Valley, relative to the Geysers CA simulations. The SC simulations and observations suggest that drainage winds are more frequent in the southeast United States than is generally recognized. They may be responsible for some of the errors in air pollution concentration predictions made by Gaussian models which assume homogeneous winds and turbulence.
2-Fold G.M. Coincidence Unit
Report containing a diagram of a coincidence unit, outlining resistor values and other technical specifications.
No abstract prepared.
2 Letters from concerned citizens in response to the recommendation regarding the Defense Office of Hearings and Appeals in Fort Meade
Community Correspondence - 2 Letters from concerned citizens in response to the recommendation regarding the Defense Office of Hearings and Appeals in Fort Meade
2-MV electrostatic quadrupole injector for heavy-ion fusion
High current and low emittance are principal requirements for heavy-ion injection into a linac driver for inertial fusion energy. An electrostatic quadrupole (ESQ) injector is capable of providing these high charge density and low emittance beams. We have modified the existing 2-MV Injector to reduce beam emittance and to double the pulse length. We characterize the beam delivered by the modified injector to the High Current Transport Experiment (HCX) and the effects of finite rise time of the extraction voltage pulse in the diode on the beam head. We demonstrate techniques for mitigating aberrations and reducing beam emittance growth in the injector.
2 MW upgrade of the Fermilab Main Injector
In January 2002, the Fermilab Director initiated a design study for a high average power, modest energy proton facility. An intensity upgrade to Fermilab's 120-GeV Main Injector (MI) represents an attractive concept for such a facility, which would leverage existing beam lines and experimental areas and would greatly enhance physics opportunities at Fermilab and in the U.S. With a Proton Driver replacing the present Booster, the beam intensity of the MI is expected to be increased by a factor of five. Accompanied by a shorter cycle, the beam power would reach 2 MW. This would make the MI a more powerful machine than the SNS or the J-PARC. Moreover, the high beam energy (120 GeV) and tunable energy range (8-120 GeV) would make it a unique high power proton facility. The upgrade study has been completed and published. This paper gives a summary report.
2-Page Summary for Neptunium solubility in the Near-field Environment of A Proposed Yucca Mountain Repository
The total system performance assessment (TSPA) for the proposed repository at Yucca Mountain, NV, includes a wide variety of processes to evaluate the potential release of radionuclides from the Engineered Barrier System into the unsaturated zone of the geosphere. The principal processes controlling radionuclide release and mobilization from the waste forms are captured in the model to assess the dissolved concentrations of radionuclides in the source-term. The TSPA model of the source-term incorporates the far-from-equilibrium dissolution of, for example, spent nuclear fuel (SNF) to capture bounding rates of radionuclide availability as the SNF degrades. In addition, for individual radionuclides, the source-term model evaluates solubility constraints that are more indicative of longer-term, equilibrium processes that can limit the potential mass transport from the source term in those cases. These solubility limits represent phase saturation and precipitation processes that can occur either at the waste form as it alters, or at other locations in the near-field environment (e.g., within the invert) if chemical conditions are different. Identification and selection of applicable constraints for solubility-limited radionuclide concentrations is a primary focus in formulating the source-term model for the TSPA. Neptunium is a long-lived radionuclide that becomes a larger fraction of the potential dose as radioactive decay of other radionuclides proceeds. To delineate appropriate long-term source-term controls on dissolved neptunium concentrations, a number of alternative models have been defined. The models are based on data both collected within the Yucca Mountain Project and taken from published literature, and have been evaluated against independent data sets to assess their applicability. The alternative models encompass ones based on precipitation of neptunium within its own separate oxide phases (i.e., ''pure'' Np-phases), and those where neptunium is incorporated into the secondary (tertiary, quaternary, etc.) uranyl phases forming as alteration products of SNF (secondary phases). The constraints on dissolved neptunium concentrations from ''pure'' Np-phase models provide more elevated bounds compared to the values derived from models of trace incorporation of Np into secondary uranyl phases. Such secondary phase models depend on a larger set of processes and variables compared to more idealized solubility models for Np-oxides, and therefore, secondary phase models generally will have more complex bases to define adequately the expected behavior.
2-pi Photoproduction from CLAS and CB-ELSA - The Search for Missing Resonances
2-pi-photoproduction is one of the promising reactions to search for baryon resonances that have been predicted but have not yet been observed. The gamma-rho --> rho-pi{sup 0}-pi{sup 0}(CB-ELSA) and the gamma-rho --> rho-pi{sup +}-pi{sup -} (CLAS) data show interesting resonance structures. A partial wave analysis (PWA) has to be done to determine which baryon resonances contribute what their quantum numbers and their relative couplings to the different accessible rho-2-pi-channels and to the photon are. First preliminary PWA-results on the lowest energy rho-pi{sup 0}-pi{sup 0} data (sq rt s<1.8 GeV)look very promising. From an extension of this analysis to higher energies combining the rho-pi{sup 0}-pi{sup 0} and the rho-pi{sup +}-pi{sup -}-data, one can expect; interesting results on resonances decaying into Delta-pi, N-rho, N(pi-pi)s, N*-pi, and Delta*-pi.
2 Questions pertaining to DON-0133 and the non-BRAC Scenario (Portsmouth Naval Shipyard)
2 Questions pertaining to DON-0133 and the non-BRAC Senario (Portsmouth Naval Shipyard. Department of Defense Clearinghouse Response: DoD Clearinghouse reply to a letter from the BRAC Commission regarding 2 Questions pertaining to DON-0133 and the non-BRAC scenario (Portsmouth Naval Shipyard)
N = 2 string amplitudes
In physics, solvable models have played very important roles. Understanding a simple model in detail teaches us a lot about more complicated models in generic situations. Five years ago, C. Vafa and I found that the closed N = 2 string theory, that is a string theory with the N = 2 local supersymmetry on the worldsheet, is classically equivalent to the self-dual Einstein gravity in four spacetime dimensions. Thus this string theory is solvable at the classical level. More recently, we have examined the N = 2 string partition function for spacial compactifications, and computed it to all order in the string perturbation expansion. The fact that such computation is possible at all suggests that the N = 2 string theory is solvable even quantum mechanically.
2 {times} 2 TeV {mu}{sup +}{mu}{sup {minus}} collider: Lattice and accelerator-detector interface study
The design for a high-luminosity {mu}{sup +}{mu}{sup {minus}} superconducting storage ring is presented based on first-pass calculations. Special attention is paid to two Iowa interaction regions (IR) whose optics are literally interlaced with the collider detectors. Various sources of backgrounds in IR are explored via realistic Monte Carlo simulations. An improved design of the collider lattice in the neighborhood of the interaction points (EP) is determined by the need to reduce significantly background levels in the detectors.
A 2 to 4 nm high power FEL on the SLAC linac
We report the results of preliminary studies of a 2 to 4 nm SASE FEL, using a photoinjector to produce the electron beam, and the SLAC linac to accelerate it to an energy up to 10 GeV. Longitudinal bunch compression is used to increases ten fold the peak current to 2.5 kA, while reducing the bunch length to the subpicosecond range. The saturated output power is in the multi-gigawatt range, producing about 10{sup 14} coherent photons within a bandwidth of about 0.2% rms, in a pulse of several millijoules. At 120Hz repetition rate the average power is about 1 W. The system is optimized for x-ray microscopy in the water window around 2 to 4 nm, and will permit imaging a biological sample in a single subpicosecond pulse.
2 x 2 Polyethylene Reflected and Moderated Highly Enriched Uranium System with Rhenium
The 2 × 2 array HEU-Re experiment was performed on the Planet universal critical assembly machine on November 4th, 2003 at the Los Alamos Critical Experiments Facility (LACEF) at Los Alamos National Laboratory (LANL). For this experiment, there were 10 ½ units, each full unit containing four HEU foils and two rhenium foils. The top unit contained only two HEU foils and two rhenium foils. A total of 42 HEU foils were used for this experiment. Rhenium is a desirable cladding material for space nuclear power applications. This experiment consisted of HEU foils interleaved with rhenium foils and is moderated and reflected by polyethylene plates. A unit consisted of a polyethylene plate, which has a recess for rhenium foils, and four HEU foils in a single layer in the top recess of each polyethylene plate. The Planet universal criticality assembly machine has been previously used in experiments containing HEU foils interspersed with SiO2 (HEU-MET-THERM-001), Al (HEU-MET-THERM-008), MgO (HEU-MET-THERM-009), Gd foils (HEU-MET-THERM-010), 2 × 2 × 26 Al (HEU-MET-THERM-012), Fe (HEU-MET-THERM-013 and HEU-MET-THERM-015), 2 × 2 × 23 SiO2 (HEU-MET-THERM-014), 2 × 2 × 11 hastalloy plates (HEU-MET-THERM-016), and concrete (HEU-MET-THERM-018). The 2 × 2 array of HEU-Re is considered acceptable for use as a benchmark critical experiment.
2 x 2 TeV mu(superscript +) mu (superscript) collider
The scenarios for high-luminosity 2 x 2 TeV and 250 x 250 GeV {mu}{sup +}{mu}{sup -} colliders are presented. Having a high physics potential, such a machine has specific physics and technical advantages and disadvantages when compared with an e{sup +}e{sup -} collider. Parameters for the candidate designs and the basic components - proton source, pion production and decay channel, cooling, acceleration and collider storage ring - are considered. Attention is paid to the areas mostly affecting the collider performance: targetry, energy spread, superconducting magnet survival, detector backgrounds, polarization, environmental issues. 13 refs., 9 figs., 4 tabs.
A 3.3 MJ, Rb + 1 driver design based on an integrated systems analysis
A computer model for systems analysis of heavy ion drivers has been developed and used to evaluate driver designs for inertial fusion energy (IFE). The present work examines a driver for a close-coupled target design that requires less total beam energy but also smaller beam spots sizes than previous target designs. Design parameters and a cost estimate for a 160 beam, 3.3 MJ driver using rubidium ions (A = 85) are reported, and the sensitivity of the results to variations in selected design parameters is given.
A 3.3 MJ, Rb{sup +1} Driver Design Based on an Integrated Systems Analysis
A computer model for systems analysis of heavy ion drivers has been developed and used to evaluate driver designs for inertial fusion energy (IFE). The present work examines a driver for a close-coupled target design that requires less total beam energy but also smaller beam spots sizes than previous target designs. Design parameters and a cost estimate for a 160 beam, 3.3 MJ driver using rubidium ions (A = 85) are reported, and the sensitivity of the results to variations in selected design parameters is given.
3,4,3-LI(1,2-HOPO): In Vitro Formation of Highly Stable Lanthanide Complexes Translates into Efficacious In Vivo Europium Decorporation
The spermine-based hydroxypyridonate octadentate chelator 3,4,3-LI(1,2-HOPO) was investigated for its ability to act as an antennae that sensitizes the emission of Sm{sup III}, Eu{sup III}, and Tb{sup III} in the Visible range (Φ{sub tot} = 0.2 - 7%) and the emission of Pr{sup III}, Nd{sup III}, Sm{sup III}, and Yb{sup III} in the Near Infra-Red range, with decay times varying from 1.78 μs to 805 μs at room temperature. The particular luminescence spectroscopic properties of these lanthanide complexes formed with 3,4,3-LI(1,2-HOPO) were used to characterize their respective solution thermodynamic stabilities as well as those of the corresponding La{sup III}, Gd{sup III}, Dy{sup III}, Ho{sup III}, Er{sup III}, Tm{sup III}, and Lu{sup III} complexes. The remarkably high affinity of 3,4,3-LI(1,2-HOPO) for lanthanide metal ions and the resulting high complex stabilities (pM values ranging from 17.2 for La{sup III} to 23.1 for Yb{sup III}) constitute a necessary but not sufficient criteria to consider this octadentate ligand an optimal candidate for in vivo metal decorporation. The in vivo lanthanide complex stability and decorporation capacity of the ligand were assessed, using the radioactive isotope {sup 152}Eu as a tracer in a rodent model, which provided a direct comparison with the in vitro thermodynamic results and demonstrated the great potential of 3,4,3-LI(1,2-HOPO) as a therapeutic metal chelating agent.
The 3.8% Medicare Contribution Tax on Unearned Income, Including Real Estate Transactions
No Description
A 3-axis force balanced accelerometer using a single proof-mass
This paper presents a new method for wideband force balancing a proof-mass in multiple axes simultaneously. Capacitive position sense and force feedback are accomplished using the same air-gap capacitors through time multiplexing. Proof of concept is experimentally demonstrated with a single-mass monolithic surface micromachined 3-axis accelerometer.
3-D Cavern Enlargement Analyses
Three-dimensional finite element analyses simulate the mechanical response of enlarging existing caverns at the Strategic Petroleum Reserve (SPR). The caverns are located in Gulf Coast salt domes and are enlarged by leaching during oil drawdowns as fresh water is injected to displace the crude oil from the caverns. The current criteria adopted by the SPR limits cavern usage to 5 drawdowns (leaches). As a base case, 5 leaches were modeled over a 25 year period to roughly double the volume of a 19 cavern field. Thirteen additional leaches where then simulated until caverns approached coalescence. The cavern field approximated the geometries and geologic properties found at the West Hackberry site. This enabled comparisons are data collected over nearly 20 years to analysis predictions. The analyses closely predicted the measured surface subsidence and cavern closure rates as inferred from historic well head pressures. This provided the necessary assurance that the model displacements, strains, and stresses are accurate. However, the cavern field has not yet experienced the large scale drawdowns being simulated. Should they occur in the future, code predictions should be validated with actual field behavior at that time. The simulations were performed using JAS3D, a three dimensional finite element analysis code for nonlinear quasi-static solids. The results examine the impacts of leaching and cavern workovers, where internal cavern pressures are reduced, on surface subsidence, well integrity, and cavern stability. The results suggest that the current limit of 5 oil drawdowns may be extended with some mitigative action required on the wells and later on to surface structure due to subsidence strains. The predicted stress state in the salt shows damage to start occurring after 15 drawdowns with significant failure occurring at the 16th drawdown, well beyond the current limit of 5 drawdowns.
3-D computer simulations of EM field sin the APS vacuum chamber. Part 2: Time-domain analysis
Our simulations suggest that the strong peak around 4 GHz in the narrow gap observed in the measurements is generated by TE modes. Therefore, one should not worry about this peak insofar as the coupling impedance is concerned. On the other hand, some discrepancies between our simulations and the measurements are noticed and remain to be resolved.
3-D computer simulations of EM fields in the APS vacuum chamber: Part 1, Frequency-domain analysis
The vacuum chamber proposed for the storage ring of the 7-GeV Advanced Photon Source (APS) basically consists of two parts: the beam chamber and the antechamber, connected to each other by a narrow gap. A sector of 1-meter-long chamber with dosed end plates, to which are attached the 1-inch-diameter beampipes centered at the beam chamber, has been built for experimental purposes. The 3-D code MAFIA has been used to simulate the frequency-domain behaviors of EM fields in this setup. The results are summarized in this note and are compared with that previously obtained from 2-D simulations and that from network analyzer measurements. They are in general agreement. A parallel analysis in the time-domain is reported in a separate note. The method of our simulations can be briefly described as follows. The 1-inch diameter beampipes are terminated by conducting walls at a length of 2 cm. The whole geometry can thus be considered as a cavity. The lowest RF modes of this geometry are computed using MAFIA. The eigenfrequencies of these modes are a direct output of the eigenvalue solver E3, whereas the type of each mode is determined by employing the postprocessor P3. The mesh sizes are chosen such that they are small enough for computations in the frequency region in which we are interested (the sampling theorem), while the total number of mesh points is still well within the range that our computer system can cope with.
No abstract prepared.
3-D elastic wave scattering by a layer containing vertical periodic fractures
No abstract prepared.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.