Some Comments on the La Primavera Geothermal Field, Mexico (open access)

Some Comments on the La Primavera Geothermal Field, Mexico

The La Primavera geothermal field is located about 20 km west of the city of Guadalajara, Jalisco, in the western part of the Mexican Neovolcanic Axis. Initial results of five deep exploration wells (down to 2000 m depth) were very promising; measured downhole temperatures exceed 300{degrees}C. During production, however, downhole temperatures dropped, and the chemistry of the fluids changed. The analysis of geologic, mineralogic, geochemical, and well completion data indicate that colder fluids flow down the wellbore from shallower aqifers cooling the upper zones of the gothermal reservoir. This problem is attributed to inadequate well completions. Doubts have arisen about continuing the exploration of the field because of the somewhat disappointing drilling results. However, a more thorough analysis of all available data indicates that a good geothermal prospect might exist below 3000 m, and that it could be successfully developed with appropriately located and completed wells.
Date: December 15, 1983
Creator: A., Bernardo Dominguez & Lippmann, Marcelo J.
System: The UNT Digital Library
Crustal Rock Fracture Mechanics for Design and Control of Artificial Subsurface Cracks in Geothermal Energy Extraction Engineering ({Gamma}-Project) (open access)

Crustal Rock Fracture Mechanics for Design and Control of Artificial Subsurface Cracks in Geothermal Energy Extraction Engineering ({Gamma}-Project)

Recently a significant role of artificial and/or natural cracks in the geothermal reservoir has been demonstrated in the literatures (Abe, H., et al., 1983, Nielson, D.L. and Hullen, J.B., 1983), where the cracks behave as fluid paths and/or heat exchanging surfaces. Until now, however, there are several problems such as a design procedure of hydraulic fracturing, and a quantitative estimate of fluid and heat transfer for reservoir design. In order to develop a design methodology of geothermal reservoir cracks, a special distinguished research project, named as ''{Lambda}-Project'', started at Tohoku University (5 years project, 1983-1988). In this project a basic fracture mechanics model of geothermal reservoir cracks is being demonstrated and its validation is being discussed both theoretically and experimentally. This paper descibes an outline of ''{Lambda}-Project''.
Date: December 15, 1983
Creator: Abe, Hiroyuki & Takahashi, Hideaki
System: The UNT Digital Library
Behaviors of Crack-Like Reservoirs by Means of Fracturing at Nigorikawa and Kakkonda Geothermal Fields (open access)

Behaviors of Crack-Like Reservoirs by Means of Fracturing at Nigorikawa and Kakkonda Geothermal Fields

A basic concept of the geothermal reservoir as a set of cracks is first presented. Extensions of subsurface cracks during well stimulation treatments at Nigorikawa(Mori) and closure operations of production well-head valves at Kakkonda are analysed and their behaviors are demonstrated based on results of long-distance AE Measurements.
Date: December 15, 1983
Creator: Abe, Hiroyuki; Takahashi, Hideaki; Nakatsuka, Katsuto; Niitsuma, Hiroaki & Takanohashi, Morihiko
System: The UNT Digital Library
Gamma-ray optical counterpart search experiment (GROCSE) (open access)

Gamma-ray optical counterpart search experiment (GROCSE)

The requirements of a gamma-ray burst optical counterpart detector are reviewed. By taking advantage of real-time notification of bursts, new instruments can make sensitive searches while the gamma-ray transient is still in progress. A wide field of view camera at Livermore National Laboratories has recently been adapted for detecting GRB optical counterparts to a limiting magnitude of 8. A more sensitive camera, capable of reaching m{sub upsilon} = 14, is under development.
Date: December 15, 1993
Creator: Akerlof, C.; Fatuzzo, M.; Lee, B.; Bionta, R.; Ledebuhr, A.; Park, H. S. et al.
System: The UNT Digital Library
Neutrino telescopes as a direct probe of supersymmetrybreaking (open access)

Neutrino telescopes as a direct probe of supersymmetrybreaking

We consider supersymmetric models where the scale of supersymmetry breaking lies between 5 x 10{sup 6} GeV and 5 x 10{sup 8} GeV. In this class of theories, which includes models of gauge mediated supersymmetry breaking, the lightest supersymmetric particle is the gravitino. The next to lightest supersymmetric particle is typically a long lived charged slepton with a lifetime between a microsecond and a second, depending on its mass. Collisions of high energy neutrinos with nucleons in the earth can result in the production of a pair of these sleptons. Their very high boost means they typically decay outside the earth. We investigate the production of these particles by the diffuse flux of high energy neutrinos, and the potential for their observation in large ice or water Cerenkov detectors. The relatively small cross-section for the production of supersymmetric particles is partially compensated for by the very long range of heavy particles. The signal in the detector consists of two parallel charged tracks emerging from the earth about 100 meters apart, with very little background. A detailed calculation using the Waxman-Bahcall limit on the neutrino flux and realistic spectra shows that km{sup 3} experiments could see as many as 4 events …
Date: December 15, 2003
Creator: Albuquerque, Ivone; Burdman, Gustavo & Chacko, Z.
System: The UNT Digital Library
Physics of arcing, and implications to sputter deposition (open access)

Physics of arcing, and implications to sputter deposition

Arcing is a well-known, unwanted discharge regime observed on the surface of sputtering targets. The discharge voltage breaks down to less than 50 V while the current jumps to elevated levels. Arcing is unwanted because it prevents uniform deposition and creates particulates. The issue of arcing has been dealt with by target surface conditioning and by using modern power supplies that have arc suppression incorporated. With increasing quality requirements in terms of uniformity of coatings, and absence of particulates, especially for electrochromic and other advanced coatings applications, the issue of arcing warrants a closer examination with the goal to find other, physics-based, and hopefully better approaches of arcing prevention. From a physics point of view, the onset of arcing is nothing else than the transition of the discharge to a cathodic arc mode, which is characterized by the ignition of non-stationary arc spots. Arc spots operate by a sequence of microexplosions, enabling explosive electron emission, as opposed to secondary electron emission. Arc spots and their fragments have a size distribution in the micrometer and sub-micrometer range, and a characteristic time distribution that has components shorter than microseconds. Understanding the ignition conditions of arc spots are of central physical interest. Spot …
Date: December 15, 2003
Creator: Anders, Andre
System: The UNT Digital Library
Heavy leptons at the SSC (open access)

Heavy leptons at the SSC

It is argued that detection of heavy leptons at the Superconducting Super Collider seems to be very difficult but perhaps not impossible. The feasibility is shown to depend critically upon the ability to identify events with W's decaying hadronically and missing transverse momentum. (LEW)
Date: December 15, 1987
Creator: Anderson, G. & Hinchliffe, I.
System: The UNT Digital Library
Regional Heat Sources and the Active and Break Phases of Boreal Summer Intraseasonal Variability (open access)

Regional Heat Sources and the Active and Break Phases of Boreal Summer Intraseasonal Variability

The boreal summer intraseasonal variability (BSISV) associated with the 30-50 day mode is represented by the co-existence of three components, poleward propagation of convection over the Indian and tropical west Pacific longitudes and eastward propagation along the equator. The hypothesis that the three components influence each other has been investigated using observed OLR, NCEP-NCAR reanalysis, and solutions from an idealized linear model. The null hypothesis is that the three components are mutually independent. Cyclostationary EOF (CsEOF) analysis is applied on filtered OLR to extract the life-cycle of the BSISV. The dominant mode of CsEOF is significantly tied to observed rainfall over the Indian subcontinent. The components of the heating patterns from CsEOF analysis serve as prescribed forcings for the linear model. This allows us to ascertain which heat sources and sinks are instrumental in driving the large-scale monsoon circulation during the BSISV life-cycle. We identify three new findings: (1) the circulation anomalies that develop as a Rossby wave response to suppressed convection over the equatorial Indian Ocean associated with the previous break phase of the BSISV precondition the ocean-atmosphere system in the western Indian Ocean and trigger the next active phase of the BSISV, (2) the development of convection over …
Date: December 15, 2003
Creator: Annamalai, H & Sperber, K R
System: The UNT Digital Library
Gas Chemistry in Geothermal Systems (open access)

Gas Chemistry in Geothermal Systems

Five new gas geothermometers are introduced. They are useful for predicting subsurface temperatures in water dominated geothermal systems. The geothermometers use data on CO{sub 2}, H{sub 2}S and H{sub 2} concentrations in fumarole steam as well as CO{sub 2}/H{sub 2} and H{sub 2}S/H{sub 2} ratios. It is demonstrated that the gas composition of fumarole steam may be used with or withour drillhole data to evaluate steam condensation in the upflow zones of geothermal systems. Uncertainty exists, however, in distinguishing between the effects of steam condensation and phase separation at elevated pressures. The gas content in steam from discharging wells and the solute content of the water phase can be used to evaluate which boiling processes lead to "excess steam" in the discharge and at which temperature this "excess steam" is added to the fluid moving through the aquifer and into the well. Examples, using field data, are given to demonstrate all the mentioned applications of geothermal chemistry.
Date: December 15, 1983
Creator: Arnorsson, Stefan & Gunnlaugsson, Einar
System: The UNT Digital Library
Analysis of Flow Cytometry DNA Damage Response Protein Activation Kinetics Following X-rays and High Energy Iron Nuclei Exposure (open access)

Analysis of Flow Cytometry DNA Damage Response Protein Activation Kinetics Following X-rays and High Energy Iron Nuclei Exposure

We developed a mathematical method to analyze flow cytometry data to describe the kinetics of {gamma}H2AX and pATF2 phosphorylations ensuing various qualities of low dose radiation in normal human fibroblast cells. Previously reported flow cytometry kinetic results for these DSB repair phospho-proteins revealed that distributions of intensity were highly skewed, severely limiting the detection of differences in the very low dose range. Distributional analysis reveals significant differences between control and low dose samples when distributions are compared using the Kolmogorov-Smirnov test. Radiation quality differences are found in the distribution shapes and when a nonlinear model is used to relate dose and time to the decay of the mean ratio of phosphoprotein intensities of irradiated samples to controls. We analyzed cell cycle phase and radiation quality dependent characteristic repair times and residual phospho-protein levels with these methods. Characteristic repair times for {gamma}H2AX were higher following Fe nuclei as compared to X-rays in G1 cells (4.5 {+-} 0.46 h vs 3.26 {+-} 0.76 h, respectively), and in S/G2 cells (5.51 {+-} 2.94 h vs 2.87 {+-} 0.45 h, respectively). The RBE in G1 cells for Fe nuclei relative to X-rays for {gamma}H2AX was 2.05 {+-} 0.61 and 5.02 {+-} 3.47, at 2 …
Date: December 15, 2010
Creator: Association, Universities Space Research; Chappell, Lori J.; Whalen, Mary K.; Gurai, Sheena; Ponomarev, Artem; Cucinotta, Francis A. et al.
System: The UNT Digital Library
Non-Salado flow and transport position paper. Revision 1 (open access)

Non-Salado flow and transport position paper. Revision 1

The US Department of Energy (DOE) is preparing to request the US Environmental Protection Agency (EPA) to certify compliance of the Waste Isolation Pilot Plant (WIPP) with long-term requirements of the environmental Radiation Protection Standards for Management and Standards for Management and Disposal of Spent Nuclear Fuel, High-Level and Transuranic Waste (40 CFR Part 191). The DOE must also demonstrate compliance with the long-term requirements of the Land Disposal Restrictions of the Resource Conservation and Recovery Act (RCRA) (40 CFR Part 268.6). Sandia National Laboratories (SNL) has ben conducting iterative performance assessments (PAs) for the the WIPP to provide guidance to the project on the technical activities required to determine long-term performance of the WIPP disposal system. The most recent PA was conducted in 1992. The objectives of this paper are to: (1) Identify and describe the relationship between non-Salado hydrology and the array of scenarios that might be relevant to the long-term performance of the repository. (2) Identify and describe the array of conceptual and mechanistic models that are required to evaluate the scenarios for the purpose of compliance. (3) Identify and describe the data/information that are required to support the conceptual and mechanistic models.
Date: December 15, 1994
Creator: Axness, C.; Beauheim, R. & Behl, Y.
System: The UNT Digital Library
A Final Focus Model for Heavy Ion Fusion Driver System Codes (open access)

A Final Focus Model for Heavy Ion Fusion Driver System Codes

The need to reach high temperatures in an inertial fusion energy (IFE) target (or a target for the study of High Energy Density Physics, HEDP) requires the ability to focus ion beams down to a small spot. System models indicate that within the accelerator, the beam radius will be of order centimeters, whereas at the final focal spot on the target, a beam radius of order millimeters is required, so radial compression factors of order ten are required. The IFE target gain (and hence the overall cost of electricity) and the HEDP target temperature are sensitive functions of the final spot radius on target. Because of this sensitivity, careful attention needs to be paid to the spot radius calculation. We review our current understanding of the elements that enter into a systems model (such as emittance growth from chromatic, geometric, and non-linear space charge forces) for the final focus based on a quadrupolar magnet system.
Date: December 15, 2004
Creator: Barnard, J. J.; Bangerter, R. O.; Henestroza, E.; Kaganovich, I. D.; Logan, B. G.; Meier, W. R. et al.
System: The UNT Digital Library
High Frequency Mechanical Pyroshock Simulations for Payload Systems (open access)

High Frequency Mechanical Pyroshock Simulations for Payload Systems

Sandia National Laboratories (SNL) designs mechanical systems with components that must survive high frequency shock environments including pyrotechnic shock. These environments have not been simulated very well in the past at the payload system level because of weight limitations of traditional pyroshock mechanical simulations using resonant beams and plates. A new concept utilizing tuned resonators attached to the payload system and driven with the impact of an airgun projectile allow these simulations to be performed in the laboratory with high precision and repeatability without the use of explosives. A tuned resonator has been designed and constructed for a particular payload system. Comparison of laboratory responses with measurements made at the component locations during actual pyrotechnic events show excellent agreement for a bandwidth of DC to 4 kHz. The bases of comparison are shock spectra. This simple concept applies the mechanical pyroshock simulation simultaneously to all components with the correct boundary conditions in the payload system and is a considerable improvement over previous experimental techniques and simulations.
Date: December 15, 1999
Creator: Bateman, Vesta I.; Brown, Frederick A.; Cap, Jerome S. & Nusser, Michael A.
System: The UNT Digital Library
Interpretation of Interference Data from the Klamath Falls, Oregon Geothermal Resource (open access)

Interpretation of Interference Data from the Klamath Falls, Oregon Geothermal Resource

Data from a seven week pressure interference test in the Klamath Falls, Oregon geothermal resource have been analyzed. The data indicate that productive wells are fed by a highly permeable fracture network and that the less permeable matrix blocks contribute significantly to the reservoir storage capacity. Detailed analysis of data from two wells is presented. Data from both of the wells yield a reservoir permeability-thickness (kh) of approximately 1.3x10{sup 6} md-ft and a storativity of 6.8x10{sup -3} ft/psi. The parameters ({lamda} and {omega}), which are determined by the distribution of permeability and storativity between the matrix and fractures, vary by more than an order of magnitude. A sensitivity study shows that for these wells, the pressure transients are not very sensitive to the distribution of permeability and storativity between the fractures and matrix blocks. No hydrologic boundaries were detected during the test. This indicates that the fault which supplies hot water to the shallow hydrothermal system does not behave according to the cassical model of either a barrier or constant potential boundary.
Date: December 15, 1983
Creator: Benson, S.M.
System: The UNT Digital Library
Development and implementation of sensitivity coefficient equations for heat conduction problems (open access)

Development and implementation of sensitivity coefficient equations for heat conduction problems

Three different methods are discussed for computing the sensitivity of the temperature field to changes in material properties and initial-boundary condition parameters for heat conduction problems. The most general method is to derive sensitivity equations by differentiating the energy equation with respect to the parameter of interest and numerically solving the resulting sensitivity equations. An example problem in which there are twelve parameters of interest is presented and the resulting sensitivity equations are derived. Numerical results are presented for thermal conductivity and volumetric heat capacity sensitivity coefficients for heat conduction in a 2-D orthotropic body. The numerical results are compared with the analytical solution to demonstrate that the numerical method is second order accurate as the mesh is refined spatially.
Date: December 15, 1997
Creator: Blackwell, B. F.; Cochran, R. J. & Dowding, K. J.
System: The UNT Digital Library
Features, Events and Processes for the Used Fuel Disposition Campaign (open access)

Features, Events and Processes for the Used Fuel Disposition Campaign

The Used Fuel Disposition (UFD) Campaign within DOE-NE is evaluating storage and disposal options for a range of waste forms and a range of geologic environments. To assess the potential performance of conceptual repository designs for the combinations of waste form and geologic environment, a master set of Features, Events, and Processes (FEPs) has been developed and evaluated. These FEPs are based on prior lists developed by the Yucca Mountain Project (YMP) and the international repository community. The objective of the UFD FEPs activity is to identify and categorize FEPs that are important to disposal system performance for a variety of disposal alternatives (i.e., combinations of waste forms, disposal concepts, and geologic environments). FEP analysis provides guidance for the identification of (1) important considerations in disposal system design, and (2) gaps in the technical bases. The UFD FEPs also support the development of performance assessment (PA) models to evaluate the long-term performance of waste forms in the engineered and geologic environments of candidate disposal system alternatives. For the UFD FEP development, five waste form groups and seven geologic settings are being considered. A total of 208 FEPs have been identified, categorized by the physical components of the waste disposal system …
Date: December 15, 2010
Creator: Blink, J. A.; Greenberg, H. R.; Caporuscio, F. A.; Houseworth, J. E.; Freeze, G. A.; Mariner, P. et al.
System: The UNT Digital Library
The Disposal Systems Evaluation Framework for DOE-NE (open access)

The Disposal Systems Evaluation Framework for DOE-NE

The Used Fuel Disposition (UFD) Campaign within DOE-NE is evaluating storage and disposal options for a range of waste forms and a range of geologic environments. For each waste form and geologic environment combination, there are multiple options for repository conceptual design. The Disposal Systems Evaluation Framework (DSEF) is being developed to formalize the development and documentation of options for each waste form and environment combination. The DSEF is being implemented in two parts. One part is an Excel workbook with multiple sheets. This workbook is designed to be user friendly, such that anyone within the UFD Campaign can use it as a guide to develop and document repository conceptual designs that respect thermal, geometric, and other constraints. The other part is an Access relational database file that will be centrally maintained to document the ensemble of conceptual designs developed with individual implementations of the Excel workbook. The DSEF Excel workbook includes sheets for waste form, environment, geometric constraints, engineered barrier system (EBS) design, thermal, performance assessment (PA), materials, cost, and fuel cycle system impacts. Each of these sheets guides the user through the process of developing internally consistent design options, and documenting the thought process. The sheets interact with …
Date: December 15, 2010
Creator: Blink, J. A.; Greenberg, H. R.; Halsey, W. G.; Jove-Colon, C.; Nutt, W. M. & Sutton, M.
System: The UNT Digital Library
A site scale model for modeling unsaturated zone processes atYucca Mountain, Nevada (open access)

A site scale model for modeling unsaturated zone processes atYucca Mountain, Nevada

None
Date: December 15, 1997
Creator: Bodvarsson, G.S.; Wu, Y.S.; Sonnenthal, E.L.; Bandurraga, T.M.; Ahlers, C.F.; Haukwa, C. et al.
System: The UNT Digital Library
Modeling Studies of Geothermal Systems with a Free Water Surface (open access)

Modeling Studies of Geothermal Systems with a Free Water Surface

Numerical simulators developed for geothermal reservoir engineering applications generally only consider systems which are saturated with liquid water and/or steam. However, most geothermal fields are in hydraulic communicatino with shallow ground water aquifers having free surface (water level), so that production or injection operations will cause movement of the surface, and of the air in the pore spaces above the water level. In some geothermal fields the water level is located hundreds of meters below the surface (e.g. Olkaria, Kenya; Bjornsson, 1978), so that an extensive so that an extensive unsaturated zone is present. In other the caprock may be very leaky or nonexistent [e.g., Klamath Falls, oregon (Sammel, 1976)]; Cerro Prieto, Mexico; (Grant et al., 1984) in which case ther eis good hydraulic communication between the geothermal reservoir and the shallow unconfined aquifers. Thus, there is a need to explore the effect of shallow free-surface aquifers on reservoir behavior during production or injection operations. In a free-surface aquifer the water table moves depending upon the rate of recharge or discharge. This results in a high overall storativity; typically two orders of magnitude higher than that of compressed liquid systems, but one or two orders of magnitude lower than that …
Date: December 15, 1983
Creator: Bodvarsson, Gudmundur S. & Pruess, K.
System: The UNT Digital Library
New Security Results on Encrypted Key Exchange (open access)

New Security Results on Encrypted Key Exchange

Schemes for encrypted key exchange are designed to provide two entities communicating over a public network, and sharing a (short) password only, with a session key to be used to achieve data integrity and/or message confidentiality. An example of a very efficient and ''elegant'' scheme for encrypted key exchange considered for standardization by the IEEE P1363 Standard working group is AuthA. This scheme was conjectured secure when the symmetric-encryption primitive is instantiated via either a cipher that closely behaves like an ''ideal cipher,'' or a mask generation function that is the product of the message with a hash of the password. While the security of this scheme in the former case has been recently proven, the latter case was still an open problem. For the first time we prove in this paper that this scheme is secure under the assumptions that the hash function closely behaves like a random oracle and that the computational Diffie-Hellman problem is difficult. Furthermore, since Denial-of-Service (DoS) attacks have become a common threat we enhance AuthA with a mechanism to protect against them.
Date: December 15, 2003
Creator: Bresson, Emmanuel; Chevassut, Olivier & Pointcheval, David
System: The UNT Digital Library
Soft Error Vulnerability of Iterative Linear Algebra Methods (open access)

Soft Error Vulnerability of Iterative Linear Algebra Methods

Devices become increasingly vulnerable to soft errors as their feature sizes shrink. Previously, soft errors primarily caused problems for space and high-atmospheric computing applications. Modern architectures now use features so small at sufficiently low voltages that soft errors are becoming significant even at terrestrial altitudes. The soft error vulnerability of iterative linear algebra methods, which many scientific applications use, is a critical aspect of the overall application vulnerability. These methods are often considered invulnerable to many soft errors because they converge from an imprecise solution to a precise one. However, we show that iterative methods can be vulnerable to soft errors, with a high rate of silent data corruptions. We quantify this vulnerability, with algorithms generating up to 8.5% erroneous results when subjected to a single bit-flip. Further, we show that detecting soft errors in an iterative method depends on its detailed convergence properties and requires more complex mechanisms than simply checking the residual. Finally, we explore inexpensive techniques to tolerate soft errors in these methods.
Date: December 15, 2007
Creator: Bronevetsky, G & de Supinski, B
System: The UNT Digital Library
Spent fuel shipping, reprocessing, and recycle fabrication in the HTGR fuel cycle (open access)

Spent fuel shipping, reprocessing, and recycle fabrication in the HTGR fuel cycle

From joint power generation conference; New Orleans, Louasiana, USA (16 Sep 1973). The High Temperature Gas Cooled Reactor (HTGR) fuel recycle operation is described. The description includes the HTGR spent fuel shipping system and the proposed method of reprocessing the spent fuel to recover the bred /sup 233/U and /sup 235/U. The process for refabricating the recovered fuel into recycle fuel is also discussed. (auth)
Date: December 15, 1972
Creator: Brooks, L. H.; Davis, C. R.; Peterman, D. D. & Spaeth, M. E.
System: The UNT Digital Library
The Cost Effectiveness of Fracture Stimulation in Increasing the Flow from Geothermal Wells (open access)

The Cost Effectiveness of Fracture Stimulation in Increasing the Flow from Geothermal Wells

The cost effectiveness of fracture stimulation at The Geysers, the Imperial Valley, and other geothermal resource areas in the United States vas studied using GEOCOM, a computer code for analyzing the impact of completion activities on the life-cycle costs of geothermal wells. Technologies for fracturing the reservoir near the wellbore involve the creation of a pressure pulse in the wellbore by means of either hydraulic or explosive force. The cost of a single fracture stimulation job can vary from $50,000 to over $500,000, with a typical cost of around $300,000. The code shows that additional flow achieved by fracture stimulation must exceed 10,000 pounds per hour for each $100,000 invested in stimulation in order for a fracture treatment to be cost effective. In some reservoirs, this additional flow must be as great as 30,000 pounds per hour. The cost effectiveness of fracturing has not yet been demonstrated in the field. The Geothermal Well Stimulation Program achieved an overall average of about 10,000 pounds per hour for each $100,000 invested.
Date: December 15, 1983
Creator: Brown, Gerald L.
System: The UNT Digital Library
Using the Schur Complement to Reduce Runtime in KULL's Magnetic Diffusion Package (open access)

Using the Schur Complement to Reduce Runtime in KULL's Magnetic Diffusion Package

Recently a Resistive Magnetohydrodynamics (MHD) package has been added to the KULL code. In order to be compatible with the underlying hydrodynamics algorithm, a new sub-zonal magnetics discretization was developed that supports arbitrary polygonal and polyhedral zones. This flexibility comes at the cost of many more unknowns per zone - approximately ten times more for a hexahedral mesh. We can eliminate some (or all, depending on the dimensionality) of the extra unknowns from the global matrix during assembly by using a Schur complement approach. This trades expensive global work for cache-friendly local work, while still allowing solution for the full system. Significant improvements in the solution time are observed for several test problems.
Date: December 15, 2010
Creator: Brunner, T A & Kolev, T V
System: The UNT Digital Library