text:
Office of Scientific and Technical Information
Effect of molybdenum plus chromium on the corrosion of iron-, nickel-, and cobalt-base alloys in basaltic lava and simulated magmatic gas at 1150/sup 0/C
SciTech Connect (OSTI)
The compatibility of several binary and ternary alloys in a magma environment was studied. Binary alloys containing molybdenum and ternary alloys containing chromium and molybdenum were exposed to basaltic lava at 1150/sup 0/C for periods of 24 and 96 hours. A cover gas was used to produce oxygen and sulfur fugacities corresponding to those of the gases dissolved in basaltic melts. Three base metals were used. These included iron, nickel, and cobalt. The primary reactions in binary alloys were found to be sulfidation. Oxide scales with a spinel layer formed on ternary alloys. The synergistic effect of molybdenum and chromium additions in ternary alloys exhibited superior corrosion resistance to binary alloys which formed base-metal sulfides down grain-boundaries. Extensive analyses of the reaction products by scanning electron microscopy, X-ray energy dispersive analysis, electron microprobe analysis, and metallography are presented for each alloys. The products formed are discussed with reference to thermodynamic stability diagrams, and the reaction path concept is used to explain some of the corrosion.
Ehrlich, S.A.; Douglass, D.L.
Basalt features observed in outcrops, cores, borehole video imagery and geophysical logs, and basalt hydrogeologic study at the Idaho National Engineering Laboratory, Eastern Idaho
A study was undertaken to examine permeable zones identified in boreholes open to the underlying basalt and to describe the vertical cross flows present in the boreholes. To understand the permeable zones in the boreholes detailed descriptions and measurements of three outcrops in the Snake River Plain, three cores at the Idaho Chemical Processing Plant (ICPP) at the INEL, and over fifty borehole TV logs from the INEL were carried out. Based on the observations made on the three outcrops an idealized basalt lava flow model was generated that used a set of nomenclature that would be standard for the basalt lava flows studied. An upper vesicular zone, a sometimes absent columnar zone, central zone, and lower vesicular zone make up the basalt lava flow model. The overall distinction between the different zones are based on the vesicle shape size, vesicularity, and fractures present. The results of the studies also indicated that the basalt lava flows at the INEL are distal to medial facies pahoehoe lava flows with close fitting contacts. The most permeable zones identified in these basalts are fractured vesiculated portions of the top of the lava flow, the columnar areas, and basalt-flow contacts in order of importance. This was determined from impeller flowmeter logging at the INEL. Having this information a detailed stratigraphy of individual basalt lava flows and the corresponding permeable units were generated. From this it was concluded that groundwater flow at the ICPP prefers to travel along thin basalt lava flows or flow-units. Flow direction and velocity of intrawell flows detected by flowmeter is controlled by a nearby pumping well.
Bennecke, W.M.
JOURNAL OF PETROLOGY VOLUME 39 NUMBER 1 PAGES 125154 1998 Chemical and Isotopic Composition of Lavas
E-Print Network [OSTI]
of the Mariana Trough extension axis. The source, reflecting the arc-like nature of the source region and isotopic compositions of igneous rocks the mantle. `True' back-arc basin basalts are generated upwelling in mature ex-basalts in the south which become increasingly similar to arc lavas tensional
Stern, Robert J.
Lithium isotopes in global mid-ocean ridge basalts Paul B. Tomascak a,*, Charles H. Langmuir b
Lithium isotopes in global mid-ocean ridge basalts Paul B. Tomascak a,*, Charles H. Langmuir b January 2008 Abstract The lithium isotope compositions of 30 well-characterized samples of glassy lavas
Langmuir, Charles H.
Pressure grouting of fractured basalt flows
This report describes a field trial of pressure grouting in basalt and the results of subsequent coring and permeability measurement activities. The objective was to show that the hydraulic conductivity of fractured basalt bedrock can be significantly reduced by pressure injection of cementitious materials. The effectiveness of the pressure grout procedure was evaluated by measuring the change in the hydraulic conductivity of the bedrock. The extent of grout penetration was established by analyzing postgrout injection drilling chips for the presence of a tracer in the grout and also by examining cores of the treated basalt. Downhole radar mapping was used to establish major lava flow patterns and follow water movement during a surface infiltration test. A site called Box Canyon, which is located northwest of the INEL, was chosen for this study due to the similarity of this surface outcrop geology to that of the underlying bedrock fracture system found at the Radioactive Waste Management Complex. This study showed that hydraulic conductivity of basalt can be reduced through pressure grouting of cementitious material.
Shaw, P.; Weidner, J.; Phillips, S.; Alexander, J.
Eruption and emplacement of flood basalt. An example from the large-volume Teepee Butte Member, Columbia River Basalt Group
Flows of the Teepee Butte Member, Grande Ronde Basalt, issued from a vent system in southeastern Washington, northeastern Oregon, and western Idaho. Three distinct basalt flows were erupted: the Limekiln Rapids flow, the Joseph Creek flow, and the Pruitt Draw flow. Together these mappable flows cover more than 52,000 km[sup 2] and have a volume exceeding 5,000 km[sup 3]. A portion of the vent system for the Joseph Creek flow is exposed in cross section in Joseph Canyon, Washington; it is one of the best preserved Columbia River Basalt Group vent complexes known. The vent complex is about 1 km in cross section, 30 m high, and composed of deposits characteristic of Hawaiian-type volcanism. The vent is asymmetrical; the eastern rampart consists of intercalated pyroclastic deposits and thin pahoehoe flows; the western rampart is composed wholly of pahoehoe flows. Flows of the Teepee Butte Member are compositionally homogeneous and were emplaced as sheet flows, each having several local flow units. Our study supports the importance of linear vent systems and the westward Palouse Slope, along with the large-volume lava flows, in controlling the distribution of Columbia River Basalt Group flows. Other factors, including the number of active fissure segments and topography, modified the shape of the flows and the number of flow units. 45 refs., 19 figs., 2 tabs.
Reidel, S.P. (Washington State Univ., Pullman (United States)); Tolan, T.L. (Portland State Univ., OR (United States))
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE:1 First Use of Energy for All Purposes (Fuel and Nonfuel),Feet) Year Jan Feb Mar Apr MayAtmospheric Optical Depth7-1D: Vegetation ProposedUsingFun withconfinementEtching. | EMSLtheIndustryMitch204 Latitude:Barbara Lava
AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE:1 First Use of Energy for All Purposes (Fuel and Nonfuel),Feet) Year Jan Feb Mar Apr MayAtmosphericNuclear Security Administration the ContributionsArms Control R&DNuclear fuel recycling in 4vs. Secondary9 0Numerical
Analytical-Numerical Modeling Of Komatiite Lava Emplacement And...
Open Energy Info (EERE)
substrates with various degrees of consolidation and water contents. Thermo-mechanical erosion becomes more effective for substrates that are increasingly unconsolidated...
Analytical-Numerical Modeling Of Komatiite Lava Emplacement And Thermal
AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home5b9fcbce19 NoPublic Utilities Address: 160 East 300Algoil JumpAltergyExperiments | OpenThe Tomoves Active| OpenErosion At
LAVA (Los Alamos Vulnerability and Risk Assessment Methodology): A conceptual framework for automated risk analysis
At Los Alamos National Laboratory, we have developed an original methodology for performing risk analyses on subject systems characterized by a general set of asset categories, a general spectrum of threats, a definable system-specific set of safeguards protecting the assets from the threats, and a general set of outcomes resulting from threats exploiting weaknesses in the safeguards system. The Los Alamos Vulnerability and Risk Assessment Methodology (LAVA) models complex systems having large amounts of ''soft'' information about both the system itself and occurrences related to the system. Its structure lends itself well to automation on a portable computer, making it possible to analyze numerous similar but geographically separated installations consistently and in as much depth as the subject system warrants. LAVA is based on hierarchical systems theory, event trees, fuzzy sets, natural-language processing, decision theory, and utility theory. LAVA's framework is a hierarchical set of fuzzy event trees that relate the results of several embedded (or sub-) analyses: a vulnerability assessment providing information about the presence and efficacy of system safeguards, a threat analysis providing information about static (background) and dynamic (changing) threat components coupled with an analysis of asset ''attractiveness'' to the dynamic threat, and a consequence analysis providing information about the outcome spectrum's severity measures and impact values. By using LAVA, we have modeled our widely used computer security application as well as LAVA/CS systems for physical protection, transborder data flow, contract awards, and property management. It is presently being applied for modeling risk management in embedded systems, survivability systems, and weapons systems security. LAVA is especially effective in modeling subject systems that include a large human component.
Smith, S.T.; Lim, J.J.; Phillips, J.R.; Tisinger, R.M.; Brown, D.C.; FitzGerald, P.D.
Cooling of Kilauea Iki lava lake
In 1959 Kilauea Iki erupted leaving a 110 to 120 m lake of molten lava in its crater. The resulting lava lake has provided a unique opportunity to study the cooling dynamics of a molten body and its associated hydrothermal system. Field measurements taken at Kilauea Iki indicate that the hydrothermal system above the cooling magma body goes through several stages, some of which are well modeled analytically. Field measurements also indicate that during most of the solidification period of the lake, cooling from above is controlled by 2-phase convection while conduction dominates the cooling of the lake from below. A summary of the field work related to the study of the cooling dynamics of Kilauea Iki is presented. Quantitative and qualitative cooling models for the lake are discussed.
Hills, R.G.
Hydrologic testing methodology and results from deep basalt boreholes
The objective of the hydrologic field-testing program is to provide data for characterization of the groundwater systems wihin the Pasco Basin that are significant to understanding waste isolation. The effort is directed toward characterizing the areal and vertical distributions of hydraulic head, hydraulic properties, and hydrochemistry. Data obtained from these studies provide input for numerical modeling of groundwater flow and solute transport. These models are then used for evaluating potential waste migration as a function of space and time. The groundwater system beneath the Hanford Site and surrounding area consists of a thick, accordantly layered sequence of basalt flows and associated sedimentary interbed that primarily occur in the upper part of the Columbia River basalt. Permeable horizons of the sequence are associated with the interbeds and the interflow zones within the basalt. The columnar interiors of a flow act as low-permeability aquitards, separating the more-permeable interflows or interbeds. This paper discusses the hydrologic field-gathering activities, specifically, field-testing methodology and test results from deep basalt boreholes.
Strait, S R; Spane, F A; Jackson, R L; Pidcoe, W W
Subduction Controls of Hf and Nd Isotopes in Lavas of the Aleutian Island Arc
The Hf and Nd isotopic compositions of 71 Quaternary lavas collected from locations along the full length of the Aleutian island arc are used to constrain the sources of Aleutian magmas and to provide insight into the geochemical behavior of Nd and Hf and related elements in the Aleutian subduction-magmatic system. Isotopic compositions of Aleutian lavas fall approximately at the center of, and form a trend parallel to, the terrestrial Hf-Nd isotopic array with {var_epsilon}{sub Hf} of +12.0 to +15.5 and {var_epsilon}{sub Nd} of +6.5 to +10.5. Basalts, andesites, and dacites within volcanic centers or in nearby volcanoes generally all have similar isotopic compositions, indicating that there is little measurable effect of crustal or other lithospheric assimilation within the volcanic plumbing systems of Aleutian volcanoes. Hafnium isotopic compositions have a clear pattern of along-arc increase that is continuous from the eastern-most locations near Cold Bay to Piip Seamount in the western-most part of the arc. This pattern is interpreted to reflect a westward decrease in the subducted sediment component present in Aleutian lavas, reflecting progressively lower rates of subduction westward as well as decreasing availability of trench sediment. Binary bulk mixing models (sediment + peridotite) demonstrate that 1-2% of the Hf in Aleutian lavas is derived from subducted sediment, indicating that Hf is mobilized out of the subducted sediment with an efficiency that is similar to that of Sr, Pb and Nd. Low published solubility for Hf and Nd in aqueous subduction fluids lead us to conclude that these elements are mobilized out of the subducted component and transferred to the mantle wedge as bulk sediment or as a silicate melt. Neodymium isotopes also generally increase from east to west, but the pattern is absent in the eastern third of the arc, where the sediment flux is high and increases from east to west, due to the presence of abundant terrigenous sediment in the trench east of the Amlia Fracture Zone, which is being subducting beneath the arc at Seguam Island. Mixing trends between mantle wedge and sediment end members become flatter in Hf-Nd isotope space at locations further west along the arc, indicating that the sediment end member in the west has either higher Nd/Hf or is more radiogenic in Hf compared to Nd. This pattern is interpreted to reflect an increase in pelagic clay relative to the terrigenous subducted sedimentary component westward along the arc. Results of this study imply that Hf does not behave as a conservative element in the Aleutian subduction system, as has been proposed for some other arcs.
Yogodzinski, Gene; Vervoort, Jeffery; Brown, Shaun Tyler; Gerseny, Megan
REFERENCES Ames, L. L. , Hanford basalt flow mineralogy, inProgress report for the Hanford Waste Isolation Project,ST-137, Atlantic Richfield Hanford Company, Richland, Wa. ,
Benson, L.V.
Seismic Structure And Seismicity Of The Cooling Lava Lake Of...
Jump to: navigation, search OpenEI Reference LibraryAdd to library Journal Article: Seismic Structure And Seismicity Of The Cooling Lava Lake Of Kilauea Iki, Hawaii Abstract The...
Combined U-Th/He and 40Ar/39Ar geochronology of post-shield lavas from the Mauna Kea and Kohala volcanoes, Hawaii
Late Quaternary, post-shield lavas from the Mauna Kea and Kohala volcanoes on the Big Island of Hawaii have been dated using the {sup 40}Ar/{sup 39}Ar and U-Th/He methods. The objective of the study is to compare the recently demonstrated U-Th/He age method, which uses basaltic olivine phenocrysts, with {sup 40}Ar/{sup 39}Ar ages measured on groundmass from the same samples. As a corollary, the age data also increase the precision of the chronology of volcanism on the Big Island. For the U-Th/He ages, U, Th and He concentrations and isotopes were measured to account for U-series disequilibrium and initial He. Single analyses U-Th/He ages for Hamakua lavas from Mauna Kea are 87 {+-} 40 ka to 119 {+-} 23 ka (2{sigma} uncertainties), which are in general equal to or younger than {sup 40}Ar/{sup 39}Ar ages. Basalt from the Polulu sequence on Kohala gives a U-Th/He age of 354 {+-} 54 ka and a {sup 40}Ar/{sup 39}Ar age of 450 {+-} 40 ka. All of the U-Th/He ages, and all but one spurious {sup 40}Ar/{sup 39}Ar ages conform to the previously proposed stratigraphy and published {sup 14}C and K-Ar ages. The ages also compare favorably to U-Th whole rock-olivine ages calculated from {sup 238}U - {sup 230}Th disequilibria. The U-Th/He and {sup 40}Ar/{sup 39}Ar results agree best where there is a relatively large amount of radiogenic {sup 40}Ar (>10%), and where the {sup 40}Ar/{sup 36}Ar intercept calculated from the Ar isochron diagram is close to the atmospheric value. In two cases, it is not clear why U-Th/He and {sup 40}Ar/{sup 39}Ar ages do not agree within uncertainty. U-Th/He and {sup 40}Ar/{sup 39}Ar results diverge the most on a low-K transitional tholeiitic basalt with abundant olivine. For the most alkalic basalts with negligible olivine phenocrysts, U-Th/He ages were unattainable while {sup 40}Ar/{sup 39}Ar results provide good precision even on ages as low as 19 {+-} 4 ka. Hence, the strengths and weaknesses of the U-Th/He and {sup 40}Ar/{sup 39}Ar methods are complimentary for basalts with ages of order 100-500 ka.
Aciego, S.M.; Jourdan, F.; DePaolo, D.J.; Kennedy, B.M.; Renne, P.R.; Sims, K.W.W.
The Influence of Basalt Layers on Seismic Wave Propagation�
are to examine the effects of basalts on seismic wave propagation and the concequent implications for imaging sedimentary structures beneath them. From studies of basalt propertiesand borehole data in connection with foreward modelling and real data, I show...
Hanssen, Peter
1974. 7. Atlantic Richfield Hanford Company, Research andGABLE MOUNTAIN BASALT CORES HANFORD NUCLEAR RESERVATION L.
Martinez-Baez, L.F.
Basics of lava-lamp convection Balzs Gyre1
marketed only since the sixties 1 . The most essential ingredients are two immis- cible fluids of densities function, the lava lamp has been used also as visual aid in geoscience courses 2 demonstrat- ing phenomena such as the transformation of energy, force and motion, adiabatic circulations in the atmosphere and oceans, or magma
J�nosi, Imre M.
Physical Constraints on Geologic CO2 Sequestration in Low-Volume Basalt Formations
Deep basalt formations within large igneous provinces have been proposed as target reservoirs for carbon capture and sequestration on the basis of favorable CO2-water-rock reaction kinetics that suggest carbonate mineralization rates on the order of 102�103 d. Although these results are encouraging, there exists much uncertainty surrounding the influence of fracture-controlled reservoir heterogeneity on commercial-scale CO2 injections in basalt formations. This work investigates the physical response of a low-volume basalt reservoir to commercial-scale CO2 injections using a Monte Carlo numerical modeling experiment such that model variability is solely a function of spatially distributed reservoir heterogeneity. Fifty equally probable reservoirs are simulated using properties inferred from the deep eastern Snake River Plain aquifer in southeast Idaho, and CO2 injections are modeled within each reservoir for 20 yr at a constant mass rate of 21.6 kg s�1. Results from this work suggest that (1) formation injectivity is generally favorable, although injection pressures in excess of the fracture gradient were observed in 4% of the simulations; (2) for an extensional stress regime (as exists within the eastern Snake River Plain), shear failure is theoretically possible for optimally oriented fractures if Sh is less than or equal to 0.70SV; and (3) low-volume basalt reservoirs exhibit sufficient CO2 confinement potential over a 20 yr injection program to accommodate mineral trapping rates suggested in the literature.
Ryan M. Pollyea; Jerry P. Fairley; Robert K. Podgorney; Travis L. McLing
Pulsatory magma supply to a phonolite lava lake Clive Oppenheimer a,
Pulsatory magma supply to a phonolite lava lake Clive Oppenheimer a, , Alexandra S. Lomakina b processes (Tazieff, 1997; Oppenheimer et al., 2004). The decadal persistence of some lava lakes, with no net, though it may well have existed for a century or more (Oppenheimer and Kyle, 2008). Since 2001, the lake
Kingsbury, Nick
Plio-Pleistocene time-averaged field in southern Patagonia recorded in lava flows
Plio-Pleistocene time-averaged field in southern Patagonia recorded in lava flows V. Mejia and N. D of 53 lava flows from southern Patagonia (latitudes 49.5°52.1°S) that include the Pali-Aike volcanic, 7 figures, 3 tables. Keywords: paleomagnetic secular variation; Patagonia; Pali-Aike Volcanic Field
Singer, Bradley S.
Initial advance of long lava flows in open channels Daisuke Takagi, Herbert E. Huppert
cooling-limited when it reaches a maximum length attainable for a given supply of lava from the vent to cooling become important. Theoretical results are applied to the study of long lava flows that descended the dynamic properties of the bulk flow. In contrast, effects due to cooling played a major role throughout
Huppert, Herbert
Hands On Science with NOAA TITLE: Plate Tectonics and Lava Lamps
and geological events?" Answer = heat from radioactive decay · Compare the light bulb in lava to Earth's internal in the lava lamp. · Plot recent earthquakes and volcanic activity on a world map and compare to the location.education.noaa.gov/Ocean_and_Coasts/Ocean_Floor_Features.html USGS Earthquake Hazards - http
Preliminary assessment of regional dispersivity of the Hanford basalts
A HANFORD SITE Yakima ) BOUNDARY REFERENCE m I REPOSITORY LOCATION. . . . COLD CREEK . 4 SYNCLINE "+-. . """. . . ~ . Richland Pasco Prosser f Kennewick WALLULA GAP Figure 1. Hanford Site location. tiirahluke Syncline 1FDRD McGEE ~ ~ Ehi... MEMBER & SADDLE MOUNTA IN 5 BASALT PRIEST RAPIDS MEMBER ROCKY COULEE FLO -2000 -3000 -4000 -750 -1000 -1250 GRANDE RONDE BASALT COHASSETT FLOW McCQY ~ CANYON UMTANUM FLOW TD = 1053 m (3396 ff) TO =1034 m (3335 N) VE =25X FRENCHMAN...
LaVenue, Arthur Marsh
Preliminary Hydrogeologic Characterization Results from the Wallula Basalt Pilot Study
The DOE's Big Sky Regional Carbon Sequestration Partnership has completed drilling the first continental flood basalt sequestration pilot borehole to a total depth (TD) of 4,110 feet on the Boise White Paper Mill property at Wallula, Washington. Site suitability was assessed prior to drilling by the 2007-2008 acquisition, processing and analysis of a four-mile, five-line three component seismic swath, which was processed as a single data-dense line. Analysis of the seismic survey data indicated a composite basalt formation thickness of {approx}8,000 feet and absence of major geologic structures (i.e., faults) along the line imaged by the seismic swath. Drilling of Wallula pilot borehole was initiated on January 13, 2009 and reached TD on April 6, 2009. Based on characterization results obtained during drilling, three basalt breccia zones were identified between the depth interval of 2,716 and 2,910 feet, as being suitable injection reservoir for a subsequent CO2 injection pilot study. The targeted injection reservoir lies stratigraphically below the massive Umtanum Member of the Grande Ronde Basalt, whose flow-interior section possesses regionally recognized low-permeability characteristics. The identified composite injection zone reservoir provides a unique and attractive opportunity to scientifically study the reservoir behavior of three inter-connected reservoir intervals below primary and secondary caprock confining zones. Drill cuttings, wireline geophysical logs, and 31one-inch diameter rotary sidewall cores provided geologic data for characterization of rock properties. XRF analyses of selected rock samples provided geochemical characterizations of the rocks and stratigraphic control for the basalt flows encountered by the Wallula pilot borehole. Based on the geochemical results, the pilot borehole was terminated in the Wapshilla Ridge 1 flow of the Grande Ronde Basalt Formation. Detailed hydrologic test characterizations of 12 basalt interflow reservoir zones and 3 flow-interior/caprock intervals were performed during drilling and immediately following reaching the final borehole drilling depth (i.e., 4,110 ft). In addition, six of the 12 basalt interflow zones were selected for detailed hydrochemical characterization. Results from the detailed hydrologic test characterization program provided the primary information on basalt interflow zone transmissivity/injectivity, and caprock permeability characteristics.
B.P. McGrail; E. C. Sullivan; F. A. Spane; D. H. Bacon; G. Hund; P. D. Thorne; C. J. Thompson; S. P. Reidel; F. S. Colwell
Simulating infiltration tests in fractured basalt at the Box Canyon Site, Idaho
The results of a series of ponded infiltration tests in variably saturated fractured basalt at Box Canyon, Idaho, were used to build confidence in conceptual and numerical modeling approaches used to simulate infiltration in fractured rock. Specifically, we constructed a dual-permeability model using TOUGH2 to represent both the matrix and fracture continua of the upper basalt flow at the Box Canyon site. A consistent set of hydrogeological parameters was obtained by calibrating the model to infiltration front arrival times in the fracture continuum as inferred from bromide samples collected from fracture/borehole intersections observed during the infiltrating tests. These parameters included the permeability of the fracture and matrix continua, the interfacial area between the fracture and matrix continua, and the porosity of the fracture continuum. To calibrate the model, we multiplied the fracture-matrix interfacial area by a factor between 0.1 and 0.01 to reduce imbibition of water from the fracture continuum into the matrix continuum during the infiltration tests. Furthermore, the porosity of the fracture continuum, as calculated using the fracture aperture inferred from pneumatic-test permeabilities, was increased by a factor of 50 yielding porosity values for the upper basalt flow in the range of 0.01 to 0.02. The fracture-continuum porosity was a highly sensitive parameter controlling the arrival times of the simulated infiltration fronts. Porosity values are consistent with those determined during the Large-Scale Aquifer Pumping and Infiltration Test at the Idaho National Engineering and Environmental Laboratory.
Unger, Andre J.A.; Faybishenko, Boris; Bodvarsson, Gudmundur S.; Simmons, Ardyth M.
TESTING MODELS FOR BASALTIC VOLCANISM: IMPLICATIONS FOR YUCCA MOUNTAIN, NEVADA Eugene Smith 1 The determination of volcanic risk to the proposed high- level nuclear waste repository at Yucca Mountain requires, then volcanism in the future may not be a significant threat to Yucca Mountain. On the other hand, if melting
Conrad, Clint
Tuesday, March 14, 2006 POSTER SESSION I: LUNAR BASALTIC VOLCANISM
Africa 003-A is a 124 g new lunar meteorite (low-Ti mare basalt) found in northern Libya in the wadi Zam lithology of a new, 124 g lunar meteorite Northeast Africa 003 found in northern Libya in the wadi Zam Zam
Rathbun, Julie A.
Volatile contents in subduction-related basaltic magmas from Central Mexico: considerations on mantle enrichment processes and low pressure magma degassing during basaltic eruptions
Melt inclusions in olivine phenocrysts from primitive lavas from five different volcanoes in the Sierra Chichinautzin (Central Mexico) were analyzed by infrared spectroscopy and electron microprobe to determine their major element chemistry...
Cervantes, Pablo
Interim reclamation report, Basalt Waste Isolation project: Boreholes, 1989
In 1968, a program was started to assess the feasibility of storing Hanford Site defense waste in deep caverns constructed in basalt. This program was expanded in 1976 to include investigations of the Hanford Site as a potential location for a mined commercial nuclear waste repository. An extensive site characterization program was begun to determine the feasibility of using the basalts beneath the Hanford Site for the repository. Site research focused primarily on determining the direction and speed of groundwater movement, the uniformity of basalt layers, and tectonic stability. Some 98 boreholes were sited, drilled, deepened, or modified by BWIP between 1977 and 1988 to test the geologic properties of the Site. On December 22, 1987, President Reagan signed into law the Nuclear Waste Policy Amendments Act of 1987, which effectively stopped all repository-related activities except reclamation of disturbed lands at the Hanford Site. This report describes the development of the reclamation program for the BWIP boreholes, its implementation, and preliminary estimates of its success. The goal of the reclamation program is to return sites disturbed by the repository program as nearly as practicable to their original conditions using native plant species. 48 refs., 28 figs., 14 tabs.
Brandt, C.A.; Rickard, W.H. Jr.; Hefty, M.G.
Hf isotope compositions of northern Luzon arc lavas suggest involvement of pelagic sediments in their source
1 Hf isotope compositions of northern Luzon arc lavas suggest involvement of pelagic sediments passing through the data intersect fields for depleted mantle and pelagic sediments suggesting. The relationship between Nd and Hf isotopic compositions in the Luzon volcanics show that the type of sediment
Boyer, Edmond
Direct numerical simulations of multiphase flow with applications to basaltic volcanism and planetary evolution
Multiphase flows are an essential component of natural systems: They affect the explosivity of volcanic eruptions, shape the landscape of terrestrial planets, and govern subsurface flow in hydrocarbon reservoirs. Advancing ...
Suckale, Jenny
Isotopic and trace element geochemistry of lavas from the northern Mariana and southern volcano arcs
Samples from submarine volcanoes and islands were analyzed for concentrations of K, Rb, Sr, Ba, REE, {sup 87}Sr/{sup 86}Sr and some selected samples for {sup 143}Nd/{sup 144}Nd. These data show strong variations along the arc, being relatively depleted in the tholeiitic and low-K calc-alkaline volcanoes of the Volcano Arc (VA) and the Mariana Central Island Province (CIP). All of the Mariana Northern Seamount Province (NSP) and Volcano arc Iwo Jima (IJ) are enriched in LIL and LREE, particularly in the northern half, where the lavas have strong shoshonitic affinities. Chemical characteristics of these lavas suggest source- or melt-mixing, with the NSP shoshonites being derived from a LIL- and LREE-enriched OIB-like source or melt, while Mariana CIP and Volcano Arc melts are derived from a depleted MORB-like mangle that has been recharged with K, Rb, Sr and Ba by hydrous fluids. Neodymium and strontium isotopic data reveal {var epsilon}{sub Nd} values ranging from +2.4 to +9.5 and {sup 87}Sr/{sup 86}Sr from 0.70320 to 0.70405. Anomalous trends of {sup 87}Sr/{sup 86}Sr and Ba/La found in some S-NSP lavas suggest that the addition of a sedimentary component may be superimposed on the two component mixing. The lavas from the Mariana and Volcano arcs, therefore, are interpreted as resulting from mixing of at least three components. The bulk of the lavas derive from an OIB-like mantle source (or melt) mixing with various proportions of a metasomatized depleted mantle source (or melt). These hybrid sources may be contaminated with minor amounts of subducted sediment and fluxed by multistage-fractionated metasomatic fluid which is derived from subducted sediment and slab after the mixing of the first two components.
Lin Pingnan.
Engineered barrier development for a nuclear waste repository in basalt: an integration of current knowledge
This document represents a compilation of data and interpretive studies conducted as part of the engineered barriers program of the Basalt Waste Isolation Project. The overall objective of these studies is to provide information on barrier system designs, emplacement and isolation techniques, and chemical reactions expected in a nuclear waste repository located in the basalts underlying the Hanford Site within the state of Washington. Backfills, waste-basalt interactions, sorption, borehole plugging, etc., are among the topics discussed.
Smith, M.J.
doi:10.1016/j.gca.2005.02.002 Sulfur diffusion in basaltic melts
doi:10.1016/j.gca.2005.02.002 Sulfur diffusion in basaltic melts CARMELA FREDA,1, * DON R. BAKER,1,2 February 3, 2005) Abstract--We measured the diffusion coefficients of sulfur in two different basaltic for sulfur diffusion in anhydrous basalts: D 2.19 10 4 exp 226.3 58.3 RT where D is the diffusion coefficient
Long, Bernard
Final Reclamation Report: Basalt Waste Isolation Project exploratory shaft site
The restoration of areas disturbed by activities of the Basalt Waste Isolation Project (BWIP) constitutes a unique operation at the US Department of Energy's (DOE) Hanford Site, both from the standpoint of restoration objectives and the time frame for accomplishing these objectives. The BWIP reclamation program comprises three separate projects: borehole reclamation, Near Surface Test Facility (NSTF) reclamation, and Exploratory Shaft Facility (ESF) reclamation. The main focus of this report is on determining the success of the revegetation effort 1 year after work was completed. This report also provides a brief overview of the ESF reclamation program. 21 refs., 7 figs., 14 tabs.
Brandt, C.A.; Rickard, W.H. Jr.
Conceptual model for regional radionuclide transport from a basalt repository site. Final draft, technical memorandum
This technical memorandum was prepared to: (1) describe a typical basalt radionuclide repository site, (2) describe geologic and hydrologic processes associated with regional radionuclide transport in basalts, (3) define the parameters required to model regional radionuclide transport from a basalt repository site, and (4) develop a ''conceptual model'' of radionuclide transport from a basalt repository site. In a general hydrological sense, basalts may be described as layered sequences of aquifers and aquitards. The Columbia River Basalt, centered near the semi-arid Pasco Basin, is considered by many to be typical basalt repository host rock. Detailed description of the flow system including flow velocities with high-low hydraulic conductivity sequences are not possible with existing data. However, according to theory, waste-transport routes are ultimately towards the Columbia River and the lengths of flow paths from the repository to the biosphere may be relatively short. There are many physical, chemical, thermal, and nuclear processes with associated parameters that together determine the possible pattern of radionuclide migration in basalts and surrounding formations. Brief process descriptions and associated parameter lists are provided. Emphasis has been placed on the use of the distribution coefficient in simulating ion exchange. The use of the distribution coefficient approach is limited because it takes into account only relatively fast mass transfer processes. In general, knowledge of hydrogeochemical processes is primitive.
Walton, W.C.; Voorhees, M.L.; Prickett, T.A.
Flood and Shield Basalts from Ethiopia: Magmas from the African Superswell
Flood and Shield Basalts from Ethiopia: Magmas from the African Superswell BRUNO KIEFFER1, ETHIOPIA 4 DEEPARTEMENT DES SCIENCES DE LA TERRE ET DE L'ENVIRONNEMENT, UNIVERSITEE LIBRE DE BRUXELLES 50 the shield volcanoes. KEY WORDS: Ethiopia; flood basalts; shield volcanism; superswell INTRODUCTION According
Demouchy, Sylvie
Effects of Biochar and Basalt Additions on Carbon Sequestration and Fluxes of Greenhouse Gases in Soils
Effects of Biochar and Basalt Additions on Carbon Sequestration and Fluxes of Greenhouse Gases Emissions--Carbon Dioxide Emissions--Sequestration and Storage--Biochar--Basalt--Organic Fertilizers, this investigation focuses on the range of potential of different soil additives to enhance sequestration and storage
Vallino, Joseph J.
Interim reclamation report: Basalt Waste Isolation Project exploration shaft site
In 1968, a program was started to assess the feasibility of storing Hanford Site defense waste in deep caverns constructed in basalt. This program was expanded in 1976 to include investigations of the Hanford Site as a potential location for a mined commercial nuclear waste repository. Extensive studies of the geotechnical aspects of the site were undertaken, including preparations for drilling a large diameter Exploratory Shaft. This report describes the development of the reclamation program for the Exploratory Shaft Facility, its implementation, and preliminary estimates of its success. The goal of the reclamation program is to return sites disturbed by the repository program as nearly as practicable to their original conditions using native plant species. 43 refs., 19 figs., 9 tabs.
Elevation of surficial sediment/basalt contact in the Subsurface Disposal Area, Idaho National Engineering Laboratory
The elevation of the surficial sediment/basalt contact at the Subsurface Disposal Area (SDA), within the Radioactive Waste Management Complex (RWMC) is presented to provide a data base for future remedial actions at this site. About 1,300 elevation data from published and unpublished reports, maps, and surveyors notes were compiled to generate maps and cross-sections of the surficial sediment/basalt contact. In general, an east to west trending depression exists in the south central portion of the SDA with basalt closer to land surface on the northern and southern boundaries of the SDA. The lowest elevation of the surficial sediment/basalt contact is 4,979 ft and the greatest is land surface at 5,012 ft. The median elevation of the sediment/basalt interface is 4,994 ft. The median depth to basalt in the SDA is 16 ft if land surface elevation is assumed to be 5,010 ft. The depth from land surface to the sediment/basalt interface ranges from 24 ft in the southeast corner of the SDA to less than 3 ft at the north-central boundary of the SDA.
Hubbell, J.M.
Dismantling the deep earth : geochemical constraints from hotspot lavas for the origin and lengthscales of mantle heterogeneity
Chapter 1 presents the first published measurements of Sr-isotope variability in olivine-hosted melt inclusions. Melt inclusions in just two Samoan basalt hand samples exhibit most of the total Sr-isotope variability ...
Jackson, Matthew G. (Matthew Gerald)
FREQUENT BASALTIC PLINIAN ERUPTIONS IN THE HISTORY OF CHIKURACHKI VOLCANO, KURILE ISLANDS, RUSSIA, Petropavlovsk- Kamchatsky, Russia; email@example.com. (2) - GEOMAR, Kiel, Germany. Chikurachki is a young, active
Belousov, Alexander
Basalt petrogenesis beneath slow- and ultraslow-spreading Arctic mid-ocean ridges
To explore the ability of melting mafic lithologies to produce alkaline ocean-island basalts (OIB), an experimental study was carried out measuring clinopyroxene (Cpx)melt and garnet (Gt)-melt partition coefficients during ...
Elkins, Lynne J
An experimental investigation of high-temperature interactions between seawater and rhyolite, andesite, basalt and peridotite
AN EXPERIMENTAL INVESTIGATION OF HIGH-TEMPERATURE INTERACTIONS BETWEEN SEAWATER AND RHYOLITE, ANDESITE, BASALT AND PERIDOTITE A Thesis by GARY WAYNE CHANDLER Submitted to the Graduate College of Texas A&M University in Partial fulfillment... of the requirement for the degree of MASTER OF SCIENCE August 1979 Major Subject: Geology AN EXPERIMENTAL INVESTIGATION OF HIGH-TEMPERATURE INTERACTIONS BETWEEN SEAWATER AND RHYOLITE, ANDESITE, BASALT AND PERIDOTITE A Thesis by GARY WAYNE CHANDLER Approved...
Chandler, Gary Wayne
Journal of Volcanology and Geothermal Research69 (1995) 105-l 16 Mount St. Helens and Santiaguito lava domes: The effect of short-
of silicic lava flows, we studied surface characteristics and obtained water content and hydrogen isotopic for these patterns is most clearly preserved in lavas erupted during early, rapid stages of dome growth. Petrologists to sample flows early in their emplacement while paying attention to surface texture, position relative
Rose, William I.
Evidence for serpentinite fluid in convergent margin systems: The example of El Salvador (Central America) arc lavas
Evidence for serpentinite fluid in convergent margin systems: The example of El Salvador (Central isotopes, has been carried out on El Salvador subduction-related lavas. The rocks have arc), Evidence for serpentinite fluid in convergent margin systems: The example of El Salvador (Central America
Doglioni, Carlo
Development and evaluation of a thermodynamic dataset for phases of interest in CO2 mineral sequestration in basaltic rocks
evaluation of a thermodynamic dataset for phases of interestKeywords: Thermodynamic dataset CO2�water� basaltABSTRACT A thermodynamic dataset describing 36 mineral
Aradottir, E.S.P.
Natural Gas Storage in Basalt Aquifers of the Columbia Basin, Pacific Northwest USA: A Guide to Site Characterization
This report provides the technical background and a guide to characterizing a site for storing natural gas in the Columbia River Basalt
Reidel, Steve P.; Spane, Frank A.; Johnson, Vernon G.
Preconceptual systems and equipment for plugging of man-made accesses to a repository in basalt
This report presents results of a study leading to preconceptual designs for plugging boreholes, shafts, and tunnels to a nuclear waste repository in basalt. Beginning design criteria include a list of preferred plug materials and plugging machines that were selected to suit the environmental conditions, and depths, diameters, and orientations of the accesses to a nuclear waste repository in the Columbia River basalts located in eastern Washington State. The environmental conditions are described. The fiscal year 1979-1980 Task II work is presented in two parts: preliminary testing of materials for plugging of man-made accesses to a repository in basalt (described in a separate report); and preconceptual systems and equipment for plugging of man-made accesses to a repository in basalt (described in this report). To fulfill the scope of the Task II work, Woodward-Clyde Consultants (WCC) was requested to: provide preconceptual systems for plugging boreholes, tunnels, and shafts in basalt; describe preconceptual borehole plugging equipment for placing the selected materials in man-made accesses; utilize the quality assurance program, program plan and schedule, and work plans previously developed for Task II; and prepare a preliminary report.
Taylor, C.L.; O'Rourke, J.E.; Allirot, D.; O'Connor, K.
Modeling of the Sedimentary Interbedded Basalt Stratigraphy for the Idaho National Laboratory Probabilistic Seismic Hazard Analysis
This report summarizes how the effects of the sedimentary interbedded basalt stratigraphy were modeled in the probabilistic seismic hazard analysis (PSHA) of the Idaho National Laboratory (INL). Drill holes indicate the bedrock beneath INL facilities is composed of about 1.1 km of alternating layers of basalt rock and loosely consolidated sediments. Alternating layers of hard rock and �soft� loose sediments tend to attenuate seismic energy greater than uniform rock due to scattering and damping. The INL PSHA incorporated the effects of the sedimentary interbedded basalt stratigraphy by developing site-specific shear (S) wave velocity profiles. The profiles were used in the PSHA to model the near-surface site response by developing site-specific stochastic attenuation relationships.
Suzette Payne
Hydrochemistry and hydrogeologic conditions within the Hanford Site upper basalt confined aquifer system
As part of the Hanford Site Ground-Water Surveillance Project, Flow System Characterization Task. Pacific Northwest Laboratory examines the potential for offsite migration of contamination within the upper basalt confined aquifer system for the US Department of Energy (DOE). As part of this activity, groundwater samples were collected over the past 2 years from selected wells completed in the upper Saddle Mountains Basalt. The hydrochemical and isotopic information obtained from these groundwater samples provides hydrologic information concerning the aquifer-flow system. Ideally, when combined with other hydrologic property information, hydrochemical and isotopic data can be used to evaluate the origin and source of groundwater, areal groundwater-flow patterns, residence and groundwater travel time, rock/groundwater reactions, and aquifer intercommunication for the upper basalt confined aquifer system. This report presents the first comprehensive Hanford Site-wide summary of hydrochemical properties for the upper basalt confined aquifer system. This report provides the hydrogeologic characteristics (Section 2.0) and hydrochemical properties (Section 3.0) for groundwater within this system. A detailed description of the range of the identified hydrochemical parameter subgroups for groundwater in the upper basalt confined aquifer system is also presented in Section 3.0. Evidence that is indicative of aquifer contamination/aquifer intercommunication and an assessment of the potential for offsite migration of contaminants in groundwater within the upper basalt aquifer is provided in Section 4.0. The references cited throughout the report are given in Section 5.0. Tables that summarize groundwater sample analysis results for individual test interval/well sites are included in the Appendix.
Spane, F.A. Jr.; Webber, W.D.
Basalt waste isolation project. Quarterly report, October 1, 1980-December 31, 1980
In September 1977, the National Waste Terminal Storage Program was restructured to support investigations of two US DOE sites - Hanford and Nevada. The Basalt Waste Isolation Project within Rockwell Hanford Operations has been chartered with the responsibility of conducting these investigations. The overall Basalt Waste Isolation Project is divided into the following principal work areas: systems integration, geosciences, hydrology, engineered barriers, near-surface test facility, engineering testing, and repository studies. Summaries of major accomplishments for each of these areas are reported in this document.
Deju, R.A.
Palaeomagnetism of flood basalts in the Pilbara Craton, Western Australia: Late Archaean continental drift and the oldest known
Palaeomagnetism of flood basalts in the Pilbara Craton, Western Australia: Late Archaean in the Nullagine Synclinorium (and Meentheena Centrocline) of the East Pilbara Basin, Western Australia, has been. Langereis, Palaeomagnetism of flood basalts in the Pilbara Craton, Western Australia: Late Archaean
Utrecht, Universiteit
Ancient sub-seafloor alteration of basaltic andesites of the Ongeluk Formation, South Africa: implications for the chemistry of
Ancient sub-seafloor alteration of basaltic andesites of the Ongeluk Formation, South Africa, Auckland Park 2006, South Africa b School of Earth Sciences, University of Leeds, Woodhouse Lane, Leeds LS2¨ttingen, Germany e Medkim 2, Du Toitspan Road, Kimberley 8301, South Africa Accepted 25 June 2003 Abstract Basaltic
Banks, David
Restoration of areas disturbed by site studies for a mined commercial radioactive waste repository: The Basalt Waste Isolation Project (BWIP)
The Basalt Waste Isolation Project (BWIP) was undertaken to environmentally characterize a portion of the US Department of Energy's Hanford Site in Washington State as a potential host for the nation's first mined commercial nuclear waste repository. Studies were terminated by Congress in 1987. Between 1976 and 1987, 72 areas located across the Hanford Site were disturbed by the BWIP. These areas include borehole pads, a large Exploratory Shaft Facility, and the Near Surface Test Facility. Most boreholes were cleared of vegetation, leveled, and stabilized with a thick layer of compacted pit-run gravel and sand. The Near Surface Test Facility consists of three mined adits, a rock-spoils bench, and numerous support facilities. Restoration began in 1988 with the objective of returning sites to pre-existing conditions using native species. The Hanford Site retains some of the last remnants of the shrub-steppe ecosystem in Washington. The primary constraints to restoring native vegetation at Hanford are low precipitation and the presence of cheatgrass, an extremely capable alien competitor. 5 figs.
Brandt, C.A.; Rickard, W.H. Jr.; Biehert, R.W.; Newell, R.L.; Page, T.L.
Upper Basalt-Confined Aquifer System in the Southern Hanford Site
The 1990 DOE Tiger Team Finding GW/CF-202 found that the hydrogeologic regime at the Hanford Site was inadequately characterized. This finding also identified the need for completing a study of the confined aquifer in the central and southern portions of the Hanford Site. The southern portion of the site is of particular interest because hydraulic-head patterns in the upper basalt-confined aquifer system indicate that groundwater from the Hanford central plateau area, where contaminants have been found in the aquifer, flows southeast toward the southern site boundary. This results in a potential for offsite migration of contaminants through the upper basalt-confined aquifer system. Based on the review presented in this report, available hydrogeologic characterization information for the upper basalt-confined aquifer system in this area is considered adequate to close the action item. Recently drilled offsite wells have provided additional information on the structure of the aquifer system in and near the southern part of the Hanford Site. Information on hydraulic properties, hydrochemistry, hydraulic heads and flow directions for the upper basalt-confined aquifer system has been re-examined and compiled in recent reports including Spane and Raymond (1993), Spane and Vermeul ( 1994), and Spane and Webber (1995).
Thorne, P.
Generation of CO2-rich melts during basalt magma ascent and degassing
Generation of CO2-rich melts during basalt magma ascent and degassing Michel Pichavant . Ida Di magma degassing, continuous decompressions of volatile-bearing (2.7-3.8 wt% H2O, 600-1300 ppm CO2 to solubilities. In contrast, the rate of vesiculation controls the final melt CO2 concentration. High
Are oxidized shergottite-like basalts an alternative to ``andesite'' on Victoria E. Hamilton
Are oxidized shergottite-like basalts an alternative to ``andesite'' on Mars? Victoria E. Hamilton; 6225 Planetology: Solar System Objects: Mars. Citation: Hamilton, V. E., and M. E. Minitti meteorites [Hamilton et al., 1997]. The mineralogies of the Martian TIR spectra were inter- preted
Hamilton, Victoria E.
Use of low frequencies for sub-basalt imaging Anton Ziolkowski,1
, Murchison House, West Mains Road, Edinburgh EH9 3LA, 3 Veritas DGC Limited, Crompton Way, Manor Royal Estate14 4HA, UK Received June 2002, revision accepted August 2002 ABSTRACT Many prospective passive ocean of passive ocean margins are covered by basalts, which are often opaque to conventional seismic reflection
Edinburgh, University of
E-Print Network 3.0 - andesitic terrains dominica Sample Search...
and sequences of Miocene basaltic andesite lavas (Tpb). Stone... depositsMCoarse, unconsolidated debris composed of local bedrock, particularly basaltic ... Source: Faulds, James...
Mineralization of Basalts in the CO2-H2O-SO2-O2 System. | EMSL
SO2-O2 System. Mineralization of Basalts in the CO2-H2O-SO2-O2 System. Abstract: Sequestering carbon dioxide (CO2) containing minor amounts of co-contaminants in geologic...
Instrumental neutron activation analysis (INAA) characterization of pre-contact basalt quarries on the American Samoan Island of Tutuila
This thesis presents a material-centered characterization of 120 geologic samples from four fine-grained basalt quarries on the Samoan Island of Tutuila. Previous unsuccessful attempts at definitive Tutuilan quarry differentiation have utilized x...
Johnson, Phillip Ray, II
Scenarios constructed for basaltic igneous activity at Yucca Mountain and vicinity; Yucca Mountain Site Characterization Project
Basaltic volcanism has been identified as a possible future event initiating a release of radionuclides from a potential repository at the proposed Yucca Mountain high-level waste repository site. The performance assessment method set forth in the Site Characterization Plan (DOE, 1988) requires that a set of scenarios encompassing all significant radionuclide release paths to the accessible environment be described. This report attempts to catalogue the details of the interactions between the features and processes produced by basaltic volcanism in the presence of the presumed groundwater flow system and a repository structure, the engineered barrier system (EBS), and waste. This catalogue is developed in the form of scenarios. We define a scenario as a well-posed problem, starting from an initiating event or process and proceeding through a logically connected and physically possible combination or sequence of features, events, and processes (FEPs) to the release of contaminants.
Barr, G.E.; Dunn, E.; Dockery, H.; Barnard, R. [Sandia National Labs., Albuquerque, NM (United States)] [Sandia National Labs., Albuquerque, NM (United States); Valentine, G.; Crowe, B. [Los Alamos National Lab., NM (United States)] [Los Alamos National Lab., NM (United States)
Summary and evaluation of hydraulic property data available for the Hanford Site upper basalt confined aquifer system
Pacific Northwest Laboratory, as part of the Hanford Site Ground-Water Surveillance Project, examines the potential for offsite migration of contamination within the upper basalt confined aquifer system. For the past 40 years, hydrologic testing of the upper basalt confined aquifer has been conducted by a number of Hanford Site programs. Hydraulic property estimates are important for evaluating aquifer flow characteristics (i.e., ground-water flow patterns, flow velocity, transport travel time). Presented are the first comprehensive Hanford Site-wide summary of hydraulic properties for the upper basalt confined aquifer system (i.e., the upper Saddle Mountains Basalt). Available hydrologic test data were reevaluated using recently developed diagnostic test analysis methods. A comparison of calculated transmissivity estimates indicates that, for most test results, a general correspondence within a factor of two between reanalysis and previously reported test values was obtained. For a majority of the tests, previously reported values are greater than reanalysis estimates. This overestimation is attributed to a number of factors, including, in many cases, a misapplication of nonleaky confined aquifer analysis methods in previous analysis reports to tests that exhibit leaky confined aquifer response behavior. Results of the test analyses indicate a similar range for transmissivity values for the various hydro-geologic units making up the upper basalt confined aquifer. Approximately 90% of the calculated transmissivity values for upper basalt confined aquifer hydrogeologic units occur within the range of 10{sup 0} to 10{sup 2} m{sup 2}/d, with 65% of the calculated estimate values occurring between 10{sup 1} to 10{sup 2} m{sup 2}d. These summary findings are consistent with the general range of values previously reported for basalt interflow contact zones and sedimentary interbeds within the Saddle Mountains Basalt.
Spane, F.A. Jr.; Vermeul, V.R.
Geology of the Wallula Gap Area, Washington. [Grande Ronde, Wanapum, and Saddle Mountains basalts
This study focuses on the structure and stratigraphy of an 80-km/sup 2/ area at the southern margin of the Pasco Basin in Wallula Gap. Field stratigraphy, petrography, natural remanent magnetism, and major-element chemistry indicate that the tholeiitic basalt flows of the Wallula Gap area correlate with units of the Grande Ronde, Wanapum, and Saddle Mountains Formations of the Yakima Subgroup of the Columbia River Basalt Group. Flows of the Frenchman Springs, Umatilla, Pomona, Elephant Mountain, and Ice Harbor Mmebers are present in the area. The Frenchman Springs Member exposed in the Wallula Gap is more than 185 m thick and consists of eight to nine flows. Its thickness and possible contemporaneous structural deformation apparently prevented emplacement of both the Roza and Priest Rapids Members at this locality. Structural uplift of the Horse Heaven Hills began prior to extrusion of the Pomona flow. Both the Pomona and Elephant Mountain Members thin and pinch out over the crest of the uplift near Mound Pond. The Ice Harbor flow was apparently confined to the basin north of the Horse Heaven uplift, but an exposure at Mound Pond suggests it flowed through Wallula Gap as an intracanyon flow. The Wallula Gap fault zone trends N65/sup 0/W and can be traced for at least 11 km along the north flank of the Horse Heaven Hills uplift. Where the fault intersects the Olympic-Wallowa Lineament at Van Sycle Canyon 8 km east of Wallula Gap, it is a broad zone of normal faulting, 300 m wide, with as much as 310 m of displacement of the basalt stratigraphy. Two faults occur in the northern portion of Van Sycle Canyon and define a graben trending N45/sup 0/W. A third fault, roughly parallel to the Wallula Gap fault, transects the The Nub and offsets 14 m of Ice Harbor basalt.
Gardner, J.N.; Snow, M.G.; Fecht, K.R.
Uranium-lead isotope systematics of Mars inferred from the basaltic shergottite QUE 94201
Uranium-lead ratios (commonly represented as {sup 238}U/{sup 204}Pb = {mu}) calculated for the sources of martian basalts preserve a record of petrogenetic processes that operated during early planetary differentiation and formation of martian geochemical reservoirs. To better define the range of {mu} values represented by the source regions of martian basalts, we completed U-Pb elemental and isotopic analyses on whole rock, mineral and leachate fractions from the martian meteorite Queen Alexandra Range 94201 (QUE 94201). The whole rock and silicate mineral fractions have unradiogenic Pb isotopic compositions that define a narrow range ({sup 206}Pb/{sup 204}Pb = 11.16-11.61). In contrast, the Pb isotopic compositions of weak HCl leachates are more variable and radiogenic. The intersection of the QUE 94201 data array with terrestrial Pb in {sup 206}Pb/{sup 204}Pb-{sup 207}Pb/{sup 204}Pb-{sup 208}Pb/{sup 204}Pb compositional space is consistent with varying amounts of terrestrial contamination in these fractions. We calculate that only 1-7% contamination is present in the purified silicate mineral and whole rock fractions, whereas the HCl leachates contain up to 86% terrestrial contamination. Despite the contamination, we are able to use the U-Pb data to determine the initial {sup 206}Pb/{sup 204}Pb of QUE 94201 (11.086 {+-} 0.008) and calculate the {mu} value of the QUE 94201 mantle source to be 1.823 {+-} 0.008. This is the lowest {mu} value calculated for any martian basalt source, and, when compared to the highest values determined for martian basalt sources, indicates that {mu} values in martian source reservoirs vary by at least 100%. The range of source {mu} values further indicates that the {mu} value of bulk silicate Mars is approximately three. The amount of variation in the {mu} values of the mantle sources ({mu} {approx} 2-4) is greater than can be explained by igneous processes involving silicate phases alone. We suggest the possibility that a small amount of sulfide crystallization may generate large extents of U-Pb fractionation during formation of the mantle sources of martian basalts.
Gaffney, A M; Borg, L E; Connelly, J N
Rate of deformation in the Pasco Basin during the Miocene as determined by distribution of Columbia River basalt flows
Detailed mapping of over 8000 square kilometers and logs from 20 core holes were used to determine the distribution and thickness of basalt flows and interbeds in the Pasco Basin. The data indicate the high-MgO Grande Ronde Basalt and Wanapum Basalt thicken from the northeast to the southwest. Deformation began in late Frenchman Springs time in the Saddle Mountains along a northwest-southeast trend and in Roza time along an east-west trend. By late Wanapum time, basalt flows were more restricted on the east side. Saddle Mountains Basalt flows spread out in the basin from narrow channels to the east. The Umatilla Member entered from the southeast and is confined to the south-central basin, while the Wilbur Creek, Asotin, Esquatzel, Pomona, and Elephant Mountain Members entered from the east and northeast. The distribution of these members is controlled by flow volume, boundaries of other flows, and developing ridges. The Wilbur Creek, Asotin, and Esquatzel flows exited from the basin in a channel along the northern margin of the Umatilla flow, while the Pomona and Elephant Mountain flows exited between Umtanum Ridge and Wallula Gap. The thickness of sedimentary interbeds and basalt flows indicated subsidence and/or uplift began in post-Grande Ronde time (14.5 million years before present) and continued through Saddle Mountains time (10.5 million years before present). Maximum subsidence occurred 40 kilometers (24 miles) north of Richland, Washington with an approximate rate of 25 meters (81 feet) per million years during the eruption of the basalt. Maximum uplift along the developing ridges was 70 meters (230 feet) per million years.
Reidel, S.P.; Ledgerwood, R.K.; Myers, C.W.; Jones, M.G.; Landon, R.D.
SECONDARY MINERALS FOUND IN CORES DC2 Al AND DC2 TAKEN FROM THE GRANDE RONDE BASALT FORMATION, PASCO BASIN, WASHINGTON
Fiscal Year 1980 for Rockwell Hanford Operations of Rockwell1. Ames, L. L. , 1976. Hanford basalt flow mineralogy.report to the Rockwell Hanford Operations Office of North
Teague, L.S.
Preliminary Analysis of Grande Ronde Basalt Formation Flow Top Transmissivity as it Relates to Assessment and Site Selection Applications for Fluid/Energy Storage and Sequestration Projects
Preliminary Analysis of Grande Ronde Basalt Formation Flow Top Transmissivity as it Relates to Assessment and Site Selection Applications for Fluid/Energy Storage and Sequestration Projects
Spane, Frank A.
Partitioning of light lithophile elements during basalt eruptions on Earth and application to Martian shergottites
/La and pressure derived from H2O-CO2 saturation (not shown); this is due to both melt 315 heterogeneity inducing variable relationships between degassing and crystallization; and to 316 the effects of post-entrapment crystallization on CO2 sequestration into a... result in a much higher Dv m than 367 observed for terrestrial basalts. The anomalies are unlikely to be due to CO2 fluxing, as 368 observed in carbon-rich systems on Earth (M�trich and Wallace, 2008), as unreasonable 369 fluxes of CO2 would...
Edmonds, Marie
NEW MEXICO BUREAU OF GEOLOGY & MINERAL RESOURCES, BULLETIN 160, 2004 Cosmogenic 36Cl ages of lava flows in the ZuniBandera
309 NEW MEXICO BUREAU OF GEOLOGY & MINERAL RESOURCES, BULLETIN 160, 2004 Cosmogenic 36Cl ages of lava flows in the ZuniBandera volcanic field, northcentral New Mexico, U.S.A. Nelia W. Dunbar1 and Fred M. Phillips2 1New Mexico Bureau of Geology and Mineral Resources, New Mexico Institute of Mining
Dunbar, Nelia W.
LAB #8 Numerical Methods
Required tools: Matlab routine dfield ; numerical routines eul, rk2, rk4; m-files. Discussion ... You should get a piecewise linear graph made up .... in the matrix y
Hydrogen Chemistry of Basalt Aquifers --Treiman et al. 282 (5397): 21... http://www.sciencemag.org/cgi/content/full/282/5397/2194e?maxtosh... 1 of 2 2/19/2008 1:26 PM
Hydrogen Chemistry of Basalt Aquifers -- Treiman et al. 282 (5397): 21... http. 2194 DOI: 10.1126/science.282.5397.2194e LETTERS Hydrogen Chemistry of Basalt Aquifers In their report "Evidence against hydrogen-based microbial ecosystems in basalt aquifers" (14 Aug., p. 976), Robert T
Lovley, Derek
Numerical Analysis Gordon K. Smyth
Numerical Analysis Gordon K. Smyth May 1997 Numerical analysis is concerned with the accurate discipline of numer ical analysis is almost entirely a product of the period since 1950 during which biostatisticians can benefit from familiarity with numerical analysis. An understanding of the numerical methods
Smyth, Gordon K.
Simulating Geologic Co-sequestration of Carbon Dioxide and Hydrogen Sulfide in a Basalt Formation
Co-sequestered CO2 with H2S impurities could affect geologic storage, causing changes in pH and oxidation state that affect mineral dissolution and precipitation reactions and the mobility of metals present in the reservoir rocks. We have developed a variable component, non-isothermal simulator, STOMP-COMP (Water, Multiple Components, Salt and Energy), which simulates multiphase flow gas mixtures in deep saline reservoirs, and the resulting reactions with reservoir minerals. We use this simulator to model the co-injection of CO2 and H2S into brecciated basalt flow top. A 1000 metric ton injection of these supercritical fluids, with 99% CO2 and 1% H2S, is sequestered rapidly by solubility and mineral trapping. CO2 is trapped mainly as calcite within a few decades and H2S is trapped as pyrite within several years.
Bacon, Diana H.; Ramanathan, Ramya; Schaef, Herbert T.; McGrail, B. Peter
Oxygen buffering of Kilauea volcanic gases and the oxygen fugacity of Kilauea basalt
Volcanic gases collected during episode 1 of the Puu Oo eruption along the east rift zone of Kilauea Volcano, Hawaii, have uniform C-O-H-S-Cl-F compositions that are sharply depleted in CO[sub 2]. The CO[sub 2]-poor gases are typical of Type II volcanic gases (GERLACH and GRAEBER, 1985) and were emitted from evolved magma stored for a prolonged period of time in the east rift zone after releasing CO[sub 2]-rich gases during an earlier period of temporary residence in the summit magma chamber. The samples are remarkably free of contamination by atmospheric gases and meteoric water. Thermodynamic evaluation of the analytical data shows that the episode 1 gases have equilibrium compositions appropriate for temperatures between 935 and 1032[degrees]C. Open- and closed-system equilibrium models of species distributions for the episode 1 gases show unequivocally that coexisting lavas buffered the gas oxygen fugacities during cooling. These models indicate that the F[sub o[sub 2
Gerlach, T.M. (Geological Survey, Vancouver, WA (United States))
Geomicrobiology of the ocean crust : the phylogenetic diversity, abundance, and distribution of microbial communities inhabiting basalt and implications for rock alteration processes
Basaltic ocean crust has the potential to host one of the largest endolithic communities on Earth. This portion of the biosphere, however, remains largely unexplored. In this study, we utilize molecular biological, ...
Santelli, Cara M
AEGIS technology demonstration for a nuclear waste repository in basalt. Assessment of effectiveness of geologic isolation systems
A technology demonstration of current performance assessment techniques as applied to a nuclear waste repository in the Columbia Plateau Basalts was conducted. Hypothetical repository coordinates were selected for an actual geographical setting on the Hanford Reservation in the state of Washington. Published hydrologic and geologic data used in the analyses were gathered in 1979 or earlier. The following report documents the technology demonstration in basalt. Available information has been used to establish the data base and initial hydrologic and geologic interpretations for this site-specific application. A simplified diagram of the AEGIS analyses is shown. Because an understanding of the dynamics of ground-water flow is essential to the development of release scenarios and consequence analyses, a key step in the demonstration is the systems characterization contained in the conceptual model. Regional and local ground-water movement patterns have been defined with the aid of hydrologic computer models. Hypothetical release scenarios have been developed and evaluated by a process involving expert opinion and a Geologic Simulation Model for basalt. (The Geologic Simulation Model can also be used to forecast future boundary conditions for the hydrologic simulation.) Chemical reactivity of the basalt with ground water will influence the leaching and transport of radionuclides; solubility equilibria based on available data are estimated with geochemical models. After the radionuclide concentrations are mathematically introduced into the ground-water movement patterns, waste movement patterns are outlined over elapsed time. Contaminant transport results are summarized for significant radionuclides that are hypothetically released to the accessible environment and to the biosphere.
Dove, F.H.; Cole, C.R.; Foley, M.G.
Site characterization report for the basalt waste isolation project. Volume II
The reference location for a repository in basalt for the terminal storage of nuclear wastes on the Hanford Site and the candidate horizons within this reference repository location have been identified and the preliminary characterization work in support of the site screening process has been completed. Fifteen technical questions regarding the qualification of the site were identified to be addressed during the detailed site characterization phase of the US Department of Energy-National Waste Terminal Storage Program site selection process. Resolution of these questions will be provided in the final site characterization progress report, currently planned to be issued in 1987, and in the safety analysis report to be submitted with the License Application. The additional information needed to resolve these questions and the plans for obtaining the information have been identified. This Site Characterization Report documents the results of the site screening process, the preliminary site characterization data, the technical issues that need to be addressed, and the plans for resolving these issues. Volume 2 contains chapters 6 through 12: geochemistry; surface hydrology; climatology, meteorology, and air quality; environmental, land-use, and socioeconomic characteristics; repository design; waste package; and performance assessment.
Numerical control user experience
AUTOCON is an acronym for Automatic Contouring, a program which generates a contour toolpath for the external profile of printed writing boards (PWB) using ICEM DDN. AUTOCON originates from the Computer Aided Process Planning (CAPP) Procedure developed by Electrical Manufacturing at Allied Bendix. The CAPP program uses electronically transmitted graphics files received from a design agency to generate automatic panel layouts and travellers for Process Engineers. The objective of AUTOCON is to utilize the information from the CAPP program to generate Numerical Control tapes. If the CAPP program can be utilized to its full potential, it is estimated that a good portion of the 60 to 90 PWB tapes shipped each month from our Numerical Control Department could be eliminated from conventional N/C Programming. The purpose of AUTOCON is to automatically generate a toolpath around the external profile of a PWB with no user intervention. The Process Engineer can generate a profile contour N/C tape after he has completed the panel layout and traveller through CAPPS for a part. After the layout and traveller illustrations have been completed, the same geometry can then be utilized again for the N/C tape. An N/C Analysist then reviews the listing for a final checkout of the job. 13 figs.
Butler, J.
Dynamical Spacetimes from Numerical Hydrodynamics
We numerically construct dynamical asymptotically-AdS$_4$ metrics by evaluating the fluid/gravity metric on numerical solutions of dissipative hydrodynamics in (2+1) dimensions. The resulting numerical metrics satisfy Einstein's equations in (3+1) dimensions to high accuracy.
Allan Adams; Nathan Benjamin; Arvin Moghaddam; Wojciech Musial
Numerical Analysis Gordon K. Smyth in Encyclopedia of Biostatistics (ISBN 0471 975761) Edited by Peter Armitage and Theodore Colton John Wiley & Sons, Ltd, Chichester, 1998 #12;Numerical Analysis Numerical analysis is concerned with the accurate and efficient evalua- tion of mathematical expressions
Summary of INEL research on the iron-enriched basalt waste form
This report summarizes the knowledge base on the iron-enriched basalt (IEB) waste form developed at the Idaho National Engineering Laboratory (INEL) during 1979--1982. The results presented discuss the applicability of IEB in converting retrieved transuranic (TRU) waste from INEL`s Radioactive Waste Management Complex (RWMC) into a vitreous/ceramic (glassy/rock) stable waste form suitable for permanent disposal in an appropriate repository, such as the Waste Isolation Pilot Plant (WIPP) in New Mexico. Borosilicate glass (BSG), the approved high-level waste form, appears unsuited for this application. Melting the average waste-soil mix from the RWMC produces the IEB composition and attempting to convert IEB to the BSG composition would require additions of substantial B{sub 2}0{sub 3}, Na, and SiO{sub 2} (glass frit). IEB requires processing temperatures of 1400 to 1600{degrees}C, depending upon the waste composition. Production of the IEB waste form, using Joule heated melters, has proved difficult in the past because of electrode and refractory corrosion problems associated with the high temperature melts. Higher temperature electric melters (arc and plasma) are available to produce this final waste form. Past research focused on extensive slag property measurements, waste form leachability tests, mechanical, composition, and microstructure evaluations, as well as a host of experiments to improve production of the waste form. Past INEL studies indicated that the IEB glass-ceramic is a material that will accommodate and stabilize a wide range of heterogeneous waste materials, including long lived radionuclides and scrap metals, while maintaining a superior level of chemical and physical performance characteristics. Controlled cooling of the molten IEB and subsequent heat treatment will produce a glass-ceramic waste form with superior leach resistance.
Reimann, G.A.; Grandy, J.D.; Eddy, T.L.; Anderson, G.L.
This report summarizes the knowledge base on the iron-enriched basalt (IEB) waste form developed at the Idaho National Engineering Laboratory (INEL) during 1979--1982. The results presented discuss the applicability of IEB in converting retrieved transuranic (TRU) waste from INEL's Radioactive Waste Management Complex (RWMC) into a vitreous/ceramic (glassy/rock) stable waste form suitable for permanent disposal in an appropriate repository, such as the Waste Isolation Pilot Plant (WIPP) in New Mexico. Borosilicate glass (BSG), the approved high-level waste form, appears unsuited for this application. Melting the average waste-soil mix from the RWMC produces the IEB composition and attempting to convert IEB to the BSG composition would require additions of substantial B{sub 2}0{sub 3}, Na, and SiO{sub 2} (glass frit). IEB requires processing temperatures of 1400 to 1600{degrees}C, depending upon the waste composition. Production of the IEB waste form, using Joule heated melters, has proved difficult in the past because of electrode and refractory corrosion problems associated with the high temperature melts. Higher temperature electric melters (arc and plasma) are available to produce this final waste form. Past research focused on extensive slag property measurements, waste form leachability tests, mechanical, composition, and microstructure evaluations, as well as a host of experiments to improve production of the waste form. Past INEL studies indicated that the IEB glass-ceramic is a material that will accommodate and stabilize a wide range of heterogeneous waste materials, including long lived radionuclides and scrap metals, while maintaining a superior level of chemical and physical performance characteristics. Controlled cooling of the molten IEB and subsequent heat treatment will produce a glass-ceramic waste form with superior leach resistance.
Numerical Modeling of HCCI Combustion
Broader source: Energy.gov (indexed) [DOE]
Numerical Modeling of HCCI Combustion Salvador M. Aceves, Daniel L. Flowers, J. Ray Smith, Joel Martinez-Frias, Francisco Espinosa-Loza, Tim Ross, Bruce Buchholz, Nick...
Depths and temperatures of <10.5?Ma mantle melting and the lithosphere-asthenosphere boundary below southern Oregon and northern California
Plagioclase and spinel lherzolite thermometry and barometry are applied to an extensive geochemical dataset of young (<10.5?Ma) primitive basaltic lavas from across Oregon's High Lava Plains, California's Modoc Plateau, ...
Till, Christy B.
Numerical Simulations in Cosmology I
The purpose of these lectures is to give a short introduction into a very vast field of numerical simulations for cosmological applications. I focus on major features of the simulations: the equations, main numerical techniques, effects of resolution, and methods of halo identification.
A. Klypin
Numerical Relativity at the Frontier
Numerical relativity is an essential tool for solving Einstein's equations of general relativity for dynamical systems characterized by high velocities and strong gravitational fields. The implementation of new algorithms that can solve these nonlinear equations in 3+1 dimensions has enabled us to tackle many long-standing problems of astrophysical interest, leading to an explosion of important new results. Numerical relativity has been used to simulate the evolution of a diverse array of physical systems, including coalescing black hole and neutron star binaries, rotating and collapsing compact objects (stars, collisionless clusters, and scalar fields), and magnetic and viscous stars, to name a few. Numerical relativity has been exploited to address fundamental points of principle, including critical phenomena and cosmic censorship. It holds great promise as a guide for interpreting observations of gravitational waves and gamma-ray bursts and identifying the sources of such radiation. Highlights of a few recent developments in numerical relativity are sketched in this brief overview.
Stuart L. Shapiro
A laser-ablation ICP-MS study of Apollo 15 low-titanium olivine-normative and quartz-normative mare basalts
differences in the mantle source to explain the $25% difference in whole-rock TiO2, and corresponding-normative basalt types, based on their different SiO2, FeO, and TiO2 whole-rock compositions. Previous studies have as originating from different lunar mantle source regions. To provide new information on the compositions
Perfect, Ed
Reinforcement of glulam beams with basalt-and glass fibre sheets Wood is one of the oldest building material's men have utilized. The strength of the timber
Reinforcement of glulam beams with basalt- and glass fibre sheets Wood is one of the oldest beams (glued laminated timber beams). By strengthening glulam beams with material who has more tensile strength than the glulam beams themselves, it distributes the load from defected areas to areas with higher
Karlsson, Brynjar
Quantitative 3-D Elemental Mapping by LA-ICP-MS of a Basaltic Clast from the Hanford 300 Area, Washington, USA
Quantitative 3-D Elemental Mapping by LA-ICP-MS of a Basaltic Clast from the Hanford 300 Area collected from the Hanford 300 Area in south-central Washington State, United States. A calibration method and riparian quality in many locations, most notably at the Hanford, Savannah River, Oak Ridge, and Nevada Test
Hu, Qinhong "Max"
Potential for Natural Gas Storage in Deep Basalt Formations at Canoe Ridge, Washington State: A Hydrogeologic Assessment
Between 1999 and 2002, Pacific Gas Transmission Company (PGT) (now TransCanada Pipeline Company) and AVISTA Corporation, together with technical support provided by the Pacific Northwest National Laboratory and the U.S. Department of Energy (DOE) examined the feasibility of developing a subsurface, natural gas-storage facility in deep, underlying Columbia River basalt in south-central Washington state. As part of this project, the 100 Circles #1 well was drilled and characterized in addition to surface studies. This report provides data and interpretations of the geology and hydrology collected specific to the Canoe Ridge site as part of the U.S. DOE funding to the Pacific Northwest National Laboratory in support of this project.
New developments in measurement technology relevant to the studies of deep geological repositories in domed salt and basalt
This report briefly describes recent geophysical and geotechnical instrumentation developments relevant to the studies of deep geologic repositories. Special emphasis has been placed on techniques that appear to minimize measurement problems associated with repositories constructed in basalt or domed salt. Included in the listing are existing measurement capabilities and deficiencies that have been identified by a few authors and instrumentation workshops that have assessed the capabilities of existing instrumentation with respect to repository applications. These deficiencies have been compared with the reported advantages and limitations of the new developments described. Based on these comparisons, areas that merit further research and development have been identified. The report is based on a thorough literature review and on discussions with several instrumentation specialists involved in instrumentation development.
Ramirez, A.L.; Mao, N.H.
Status of iron-enriched basalt as a medium for nuclear waste immobilization: a report by an independent peer review panel
The purpose of the Peer Review Panel was to provide an independent review by experts in nuclear waste processing and materials on the adequacy of the existing data base for the iron-enriched basalt waste form developed by EG and G Idaho, and to evaluate the broad range of proposed applications for this waste form. It was not the purpose of this review to specifically rank iron-enriched basalt against other nuclear waste forms. It was the concensus of the Peer Review Panel that the concept, experimental research, and identification of potential applications of the iron-enriched basalt waste form were of high quality. Iron-enriched basalt is a primarily ceramic waste form with a residual glass phase. It has a broad range of composition, permitting the incorporation of a wide variety of nuclear wastes. The product has good mechanical strength and produces very low quantities of respirable particles under impact conditions. Matrix dissolution rates under neutral pH conditions are comparable to or lower than those of borosilicate glass. In the area of waste form characterization, the Panel recommended additional static and dynamic leaching tests as a function of pH and CO/sub 2/ in solution, and in brine solutions of varying composition. The panel also recommended that unprocessed transuranic (TRU) wastes be subjected to leach tests. Large-scale iron-enriched basalt castings in which the grain growth was uncontrolled have been observed to be less durable than controlled-grain-growth laboratory-scale castings. Therefore, the Panel also recommended leaching tests as a function of microstructure to determine ranges of acceptable microstructure. In the area of the IEB production process, the Panel recommended a variety of laboratory-scale and pilot plant-scale research.
Palmour, H. III; Dosch, R.G.; Macedo, P.B.; Machiels, A.J.; Owen, D.E.
Physical and Numerical Space Running Head: Biases in Physical and Numerical Space
Physical and Numerical Space 1 Running Head: Biases in Physical and Numerical Space Elementary school children's attentional biases in physical and numerical space Tilbe Göksun: April 30, 2012 #12; Physical and Numerical Space 2 Abstract Numbers
Chatterjee, Anjan
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng. 2000; 47:15691603
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng. 2000; 47:15691603
Thompson, Lonny L.
Numerical methods in heat transfer
This book contains nine papers. Some of the titles are: Numerical calculation of bubble growth in nucleate boiling from inception through departure; An evaluation of a translator for finite element data to resistor/capacitor data for the heat diffusion; Thermophoretic deposition due to jet impingement on an inclined plane; and A three-dimensional boundary-fitted coordinate system.
Emery, A.F.; Douglass, R.W.
Numerical Modelling of Interaction between
plasma in arc furnace used in toxic waste destruction Plasma etching semiconductor High intensity arc lamp Electrode temperature after 1ms of arcing with power density of 3x109 W.m-2 Electrode temperatureNumerical Modelling of Interaction between the Electric Arc and Electrodes Principal researcher: W
Sóbester, András
THE ORBITAL PHASES AND SECONDARY TRANSITS OF KEPLER-10b. A PHYSICAL INTERPRETATION BASED ON THE LAVA-OCEAN PLANET MODEL
The Kepler mission has made an important observation: the first detection of photons from a terrestrial planet by observing its phase curve (Kepler-10b). This opens a new field in exoplanet science: the possibility of obtaining information about the atmosphere and surface of rocky planets, objects of prime interest. In this Letter, we apply the Lava-ocean model to interpret the observed phase curve. The model, a planet without atmosphere and a surface partially made of molten rocks, has been proposed for planets of the class of CoRoT-7b, i.e., rocky planets very close to their star (at a few stellar radii). Kepler-10b is a typical member of this family. It predicts that the light from the planet has an important emission component in addition to the reflected one, even in the Kepler spectral band. Assuming an isotropical reflection of light by the planetary surface (Lambertian-like approximation), we find that a Bond albedo of {approx}50% can account for the observed amplitude of the phase curve, as opposed to a first attempt where an unusually high value was found. We propose a physical process to explain this still large value of the albedo. The overall interpretation can be tested in the future with instruments such as the James Webb Space Telescope or the Exoplanet Characterization Observatory. Our model predicts a spectral dependence that is clearly distinguishable from that of purely reflected light and from that of a planet at a uniform temperature.
Rouan, D. [LESIA, UMR 8109 CNRS, Observatoire de Paris, UVSQ, Universite Paris-Diderot, 5 pl. J. Janssen, 92195 Meudon (France); Deeg, H. J. [Instituto de Astrofisica de Canarias, E-38205 La Laguna, Tenerife (Spain); Demangeon, O.; Samuel, B.; Cavarroc, C.; Leger, A. [Institut d'Astrophysique Spatiale, Universite Paris-Sud, bat 121, F-91405 Orsay (France); Fegley, B., E-mail: firstname.lastname@example.com [Planetary Chemistry Laboratory, McDonnell Center for the Space Sciences, Department of Earth and Planetary Sciences, Washington University, St. Louis, MO (United States)
spe438-20 page 1 Garrison, N.J., Busby, C.J., Gans, P.B., Putirka, K., and Wagner, D.L., 2008, A mantle plume beneath California? The mid-Miocene Lovejoy flood basalt, northern
-Miocene Lovejoy flood basalt, northern California Noah J. Garrison Cathy J. Busby Phillip B. Gans Department the eastern Snake River Plain toward the Yellowstone caldera (Armstrong et al., 1975; Rodgers et al., 1990
Busby, Cathy
Numerical integration of variational equations
We present and compare different numerical schemes for the integration of the variational equations of autonomous Hamiltonian systems whose kinetic energy is quadratic in the generalized momenta and whose potential is a function of the generalized positions. We apply these techniques to Hamiltonian systems of various degrees of freedom, and investigate their efficiency in accurately reproducing well-known properties of chaos indicators like the Lyapunov Characteristic Exponents (LCEs) and the Generalized Alignment Indices (GALIs). We find that the best numerical performance is exhibited by the \\textit{`tangent map (TM) method'}, a scheme based on symplectic integration techniques which proves to be optimal in speed and accuracy. According to this method, a symplectic integrator is used to approximate the solution of the Hamilton's equations of motion by the repeated action of a symplectic map $S$, while the corresponding tangent map $TS$, is used for the integration of the variational equations. A simple and systematic technique to construct $TS$ is also presented.
Ch. Skokos; E. Gerlach
Transient Inverse Calibration of Site-Wide Groundwater Model to Hanford Operational Impacts from 1943 to 1996--Alternative Conceptual Model Considering Interaction with Uppermost Basalt Confined Aquifer
The baseline three-dimensional transient inverse model for the estimation of site-wide scale flow parameters, including their uncertainties, using data on the transient behavior of the unconfined aquifer system over the entire historical period of Hanford operations, has been modified to account for the effects of basalt intercommunication between the Hanford unconfined aquifer and the underlying upper basalt confined aquifer. Both the baseline and alternative conceptual models (ACM-1) considered only the groundwater flow component and corresponding observational data in the 3-Dl transient inverse calibration efforts. Subsequent efforts will examine both groundwater flow and transport. Comparisons of goodness of fit measures and parameter estimation results for the ACM-1 transient inverse calibrated model with those from previous site-wide groundwater modeling efforts illustrate that the new 3-D transient inverse model approach will strengthen the technical defensibility of the final model(s) and provide the ability to incorporate uncertainty in predictions related to both conceptual model and parameter uncertainty. These results, however, indicate that additional improvements are required to the conceptual model framework. An investigation was initiated at the end of this basalt inverse modeling effort to determine whether facies-based zonation would improve specific yield parameter estimation results (ACM-2). A description of the justification and methodology to develop this zonation is discussed.
Vermeul, Vincent R.; Cole, Charles R.; Bergeron, Marcel P.; Thorne, Paul D.; Wurstner, Signe K.
Seismic Characterization of Basalt Topography at Two Candidate Sites for the INL Remote-Handled Low-Level Waste Disposal Project
This report presents the seismic refraction results from the depth to bed rock surveys for two areas being considered for the Remote-Handled Low-Level Waste (RH-LLW) disposal facility at the Idaho National Laboratory. The first area (Site 5) surveyed is located southwest of the Advanced Test Reactor Complex and the second (Site 34) is located west of Lincoln Boulevard near the southwest corner of the Idaho Nuclear Technology and Engineering Center (INTEC). At Site 5, large area and smaller-scale detailed surveys were performed. At Site 34, a large area survey was performed. The purpose of the surveys was to define the topography of the interface between the surficial alluvium and underlying basalt. Seismic data were first collected and processed using seismic refraction tomographic inversion. Three-dimensional images for both sites were rendered from the data to image the depth and velocities of the subsurface layers. Based on the interpreted top of basalt data at Site 5, a more detailed survey was conducted to refine depth to basalt. This report briefly covers relevant issues in the collection, processing and inversion of the seismic refraction data and in the imaging process. Included are the parameters for inversion and result rendering and visualization such as the inclusion of physical features. Results from the processing effort presented in this report include fence diagrams of the earth model, for the large area surveys and iso-velocity surfaces and cross sections from the detailed survey.
Jeff Sondrup; Gail Heath; Trent Armstrong; Annette Shafer; Jesse Bennett; Clark Scott
Uranium potential of southwestern New Mexico (southern Hidalgo County), including observations on crystallization history of lavas and ash tuffs and the release of uranium from them. Final report
Geological environments present in southwestern New Mexico include thick sequences of sedimentary rock including limestone, conglomerates, sandstone, and shale: igneous intrusions with associated metal deposits; caldera centers, margins, and outflow facies; and basins with marginal faults and thick late Cenozoic sedimentary fillings. Predominant rock types are Paleozoic carbonates, Mesozoic terrigeneous rocks and carbonates, and Cenozoic volcanic rocks and basin-filling terrigeneous rocks. Consideration of information available in Preliminary Reconnaissance Reports and in Hydrogeochemical and Stream Reconnaissance Reports together with 347 new whole rock chemical analyses points to three areas of anomalous uranium abundance in Hidalgo County, New Mexico. The area has experienced three major periods of igneous activity in Phanerozoic time: one associated with the Laramide cycle of the Late Cretaceous and early Tertiary, mid-Tertiary cycle of silicic volcanism with abundant calderas, and a late Tertiary cycle of mafic volcanism. Silicic volcanic rocks are the most common exposed rock type in the area, and the most enriched in uranium (range, 0.4 to 19 ppM). The most likely source for any uranium ore-forming solutions lies with this cycle of volcanism. Solutions might have been introduced during volcanism or formed later by groundwater leaching of cooled volcanic rocks. Results indicate that groundwater leaching of cooled volcanic rocks was not an effective means of mobilizing uranium in the area. Study of several rhyolite lava flows indicates that they were emplaced in supercooled condition and may have crystallized completely at temperatures well below their liquids, or they may have warmed as crystallization released latent heat. Statistical comparison of the uranium concentration revealed no differences between vitrophyres and associated felsites.
Walton, A.W.; Salter, T.L.; Zetterlund, D.
Disruptive Innovation in Numerical Hydrodynamics
We propose the research and development of a high-fidelity hydrodynamic algorithm for tetrahedral meshes that will lead to a disruptive innovation in the numerical modeling of Laboratory problems. Our proposed innovation has the potential to reduce turnaround time by orders of magnitude relative to Advanced Simulation and Computing (ASC) codes; reduce simulation setup costs by millions of dollars per year; and effectively leverage Graphics Processing Unit (GPU) and future Exascale computing hardware. If successful, this work will lead to a dramatic leap forward in the Laboratory's quest for a predictive simulation capability.
Waltz, Jacob I. [Los Alamos National Laboratory
Numerical Simulations of Bouncing Jets
Bouncing jets are fascinating phenomenons occurring under certain conditions when a jet impinges on a free surface. This effect is observed when the fluid is Newtonian and the jet falls in a bath undergoing a solid motion. It occurs also for non-Newtonian fluids when the jets falls in a vessel at rest containing the same fluid. We investigate numerically the impact of the experimental setting and the rheological properties of the fluid on the onset of the bouncing phenomenon. Our investigations show that the occurrence of a thin lubricating layer of air separating the jet and the rest of the liquid is a key factor for the bouncing of the jet to happen. The numerical technique that is used consists of a projection method for the Navier-Stokes system coupled with a level set formulation for the representation of the interface. The space approximation is done with adaptive finite elements. Adaptive refinement is shown to be very important to capture the thin layer of air that is responsible for the bouncing.
Bonito, Andrea; Lee, Sanghyun
Numerical simulation of electrokinetically driven micro flows
Spectral element based numerical solvers are developed to simulate electrokinetically driven flows for micro-fluidic applications. Based on these numerical solvers, basic phenomena and devices for electrokinetic applications in micro and nano flows...
Hahm, Jungyoon
An implicit numerical algorithm general relativistic hydrodynamics
An implicit numerical algorithm general relativistic hydrodynamics This article has been replaced by arXiv:0801.1017
A. Hujeirat
Satellite Application Facility for Numerical Weather Prediction
into Numerical Weather Prediction (NWP) models. However, the impact of such observations often critNWP SAF Satellite Application Facility for Numerical Weather Prediction Document NWPSAF-KN-VS-002 Stoffelen KNMI #12;NWP SAF Satellite Application Facility for Numerical Weather Prediction Report
Stoffelen, Ad
Amending Numerical Weather Prediction forecasts using GPS
to validate the amounts of humidity in Numerical Weather Prediction (NWP) model forecasts. This paper presents. Satellite images and Numerical Weather Prediction (NWP) models are used together with the synoptic surface. In this paper, a case is presented for which the operational Numerical Weather Prediction Model (NWP) HIRLAM
Lava Dome | Open Energy Information
AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data CenterFranconia, Virginia: Energy Resources Jump to:46 - 429Lacey, Washington:Lakeville, MN)Lauderhill, Florida: EnergyLaurelDome Jump to:
Numerical Simulations of Thermobaric Explosions
A Model of the energy evolution in thermobaric explosions is presented. It is based on the two-phase formulation: conservation laws for the gas and particle phases along with inter-phase interaction terms. It incorporates a Combustion Model based on the mass conservation laws for fuel, air and products; source/sink terms are treated in the fast-chemistry limit appropriate for such gas dynamic fields. The Model takes into account both the afterburning of the detonation products of the booster with air, and the combustion of the fuel (Al or TNT detonation products) with air. Numerical simulations were performed for 1.5-g thermobaric explosions in five different chambers (volumes ranging from 6.6 to 40 liters and length-to-diameter ratios from 1 to 12.5). Computed pressure waveforms were very similar to measured waveforms in all cases - thereby proving that the Model correctly predicts the energy evolution in such explosions. The computed global fuel consumption {mu}(t) behaved as an exponential life function. Its derivative {dot {mu}}(t) represents the global rate of fuel consumption. It depends on the rate of turbulent mixing which controls the rate of energy release in thermobaric explosions.
Kuhl, A L; Bell, J B; Beckner, V E; Khasainov, B
Numerical Investigations of Magnetohydrodynamic Hypersonic Flows.
??Numerical simulations of magnetohydrodynamic (MHD) hypersonic flow are presented for both laminar and turbulent flow over a cylinder and flow entering a scramjet inlet. ANSYS� (more)
Guarendi, Andrew N
Adjoint Sensitivity Analysis for Numerical Weather Prediction
Sep 2, 2011 ... Adjoint Sensitivity Analysis for Numerical Weather Prediction: Applications to Power Grid Optimization. Alexandru Cioaca(alexgc ***at*** vt.edu)
Alexandru Cioaca
February 25, 2005 NUMERICAL ENERGY CONSERVATION FOR
February 25, 2005 NUMERICAL ENERGY CONSERVATION FOR MULTIFREQUENCY OSCILLATORY DIFFERENTIAL systems, modulated Fourier expansion, energy conservation, oscillatory solutions. 1 Introduction
Tübingen, Universität
February 6, 2004 NUMERICAL ENERGY CONSERVATION FOR
February 6, 2004 NUMERICAL ENERGY CONSERVATION FOR MULTIFREQUENCY OSCILLATORY DIFFERENTIAL, Hamiltonian systems, modulated Fourier expansion, energy conservation, oscillatory solutions. 1 Introduction
Hairer, Ernst
Using MATLAB to solve differential equations numerically
Using MATLAB to solve differential equations numerically Morten Brøns Department of Mathematics of the programming language MATLAB. We will focus on practical matters and readers interested in numerical analysis as a mathematical subject should look elsewhere. In the G-databar at DTU, MATLAB can be accessed either by typing �?
Klein, David
Numerical Study of a Turbulent Hydraulic Jump
Numerical Study of a Turbulent Hydraulic Jump Qun Zhao, Shubhra Misra, Ib. A. Svendsen and James T of a Turbulent Hydraulic Jump � p.1/14 #12;Objective Our ultimate goal is to study the breaking waves. Numerical Study of a Turbulent Hydraulic Jump � p.2/14 #12;A moving bore Qiantang Bore China (Courtesy of Dr J
Zhao, Qun
Toward Practical Aerodynamic Design Through Numerical Optimization
Toward Practical Aerodynamic Design Through Numerical Optimization David W. Zingg, and Laura6, Canada A Newton-Krylov algorithm for aerodynamic optimization is applied to the multipoint design aerodynamic design. I. Introduction Beginning with the work of Hicks et al.1 and Hicks and Henne,2 numerical
Zingg, David W.
The Joint Essential Numerical Range of operators
The Joint Essential Numerical Range of operators: Convexity and Related Results Chi-Kwong Li Classification 47A12, 47A13, 47A55. Keywords Joint essential numerical range, self-adjoint operator, Hilbert the joint behavior of several operators A1, . . . , Am. One may see [1, 5, 12, 14, 15, 16, 19, 23, 28, 31
Li, Chi-Kwong
Numerical Simulation of Cooling Gas Injection Using
Numerical Simulation of Cooling Gas Injection Using Adaptive Multiscale Techniques Wolfgang Dahmen: finite volume method, film cooling, cooling gas injection, multiscale techniques, grid adaptation firstname.lastname@example.com (Thomas Gotzen) #12;Numerical simulation of cooling gas injection using adaptive multiscale techniques
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2012; 68:377402
of its extension to operational 3-D numerical weather prediction models. Copyright 2011 John Wiley & Sons complicated Numerical Weather Prediction (NWP) models. Therefore, derivation and testing of various algorithms by Chen et al. (Int. J. Numer. Meth. Fluids 2009), we attempt to obtain a reduced order model of the above
Navon, Michael
Sandia National Laboratories: Numerical Manufacturing And Design...
NuMAD (Numerical Manufacturing And Design) is an open-source software tool written in Matlab which simplifies the process of creating a three-dimensional model of a wind turbine...
High performance computing and numerical modelling
Numerical methods play an ever more important role in astrophysics. This is especially true in theoretical works, but of course, even in purely observational projects, data analysis without massive use of computational methods has become unthinkable. The key utility of computer simulations comes from their ability to solve complex systems of equations that are either intractable with analytic techniques or only amenable to highly approximative treatments. Simulations are best viewed as a powerful complement to analytic reasoning, and as the method of choice to model systems that feature enormous physical complexity such as star formation in evolving galaxies, the topic of this 43rd Saas Fee Advanced Course. The organizers asked me to lecture about high performance computing and numerical modelling in this winter school, and to specifically cover the basics of numerically treating gravity and hydrodynamics in the context of galaxy evolution. This is still a vast field, and I necessarily had to select a subset ...
13.024 Numerical Marine Hydrodynamics, Spring 2003
Introduction to numerical methods: interpolation, differentiation, integration, systems of linear equations. Solution of differential equations by numerical integration, partial differential equations of inviscid hydrodynamics: ...
Milgram, Jerome H.
Developments in Petroleum Science, 6 FUNDAMENTALS OF NUMERICAL
Developments in Petroleum Science, 6 FUNDAMENTALS OF NUMERICAL RESERVOIR SIMULATION DONALD WCongressCatalogingin PublicationData Peaceman, Donald W Fundamentals of numerical reservoir simulation. (develrpents in petroleum
Santos, Juan
Effective velocities in fractured media: a numerical study using the ...
and precise numerical study of effective velocities in fractured structures. ... In this paper, we ..... A final result is that our numerical simulations of P-, SV- and.
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2011; 00:123
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2011; 00 for the treatment of discontinuous pressures in multifluid flows Roberto F. Ausas1 , Gustavo C. Buscaglia1 WORDS: Multifluids, Twophase flows, Embedded interfaces, Finite element method, Surface tension
Buscaglia, Gustavo C.
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 31: 345358 (1999)
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 31: 345358 (1999 AND S. ROYc,3 a Engineering Science Program, MAES Department, 316A Perkins Hall, Uni6ersity of Tennessee Corporation, Burr Ridge, IL, USA SUMMARY The quest continues for accurate and efficient computational fluid
Roy, Subrata
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING Int. J. Numer. Meth. Biomed. Engng. 2014; 30:15971613
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING Int. J. Numer. Meth. Biomed BIOMEDICAL ENGINEERING 2013 A finite strain nonlinear human mitral valve model with fluid. © 2014 The Authors. International Journal for Numerical Methods in Biomedical Engineering published
Luo, Xiaoyu
Symplectic numerical integration of Hamiltonian systems
This paper describes some general techniques available for symplectic or Lie-Poisson integration and illustrate the results with some numerical computations. In this spirit, I also discuss reversible integration, equivariant integration, integration of volume preserving flows, and symplectic cellular automata. My intention is not to be exhaustive but to give a representative review. 76 refs., 5 figs.
Scovel, C.
Numerical likelihood analysis of cosmic ray anisotropies
A numerical likelihood approach to the determination of cosmic ray anisotropies is presented which offers many advantages over other approaches. It allows a wide range of statistically meaningful hypotheses to be compared even when full sky coverage is unavailable, can be readily extended in order to include measurement errors, and makes maximum unbiased use of all available information.
Carlos Hojvat et al.
Numerical techniques for lattice gauge theories
The motivation for formulating gauge theories on a lattice is reviewed. Monte Carlo simulation techniques are then discussed for these systems. Finally, the Monte Carlo methods are combined with renormalization group analysis to give strong numerical evidence for confinement of quarks by non-Abelian gauge fields.
Creutz, M.
Numerical Valuation of Discrete Barrier Options with
Numerical Valuation of Discrete Barrier Options with the Adaptive Mesh Model and Other Competing for discrete barrier options such that many methods have been suggested and declared to price discrete barrier options fast and accurately but no one can tell exactly that what method is the best. We also make
Chu, Hao-hua
GPU Acceleration of Numerical Weather John Michalakes
GPU Acceleration of Numerical Weather Prediction John Michalakes National Center for Atmospheric parallelism will prove ineffective for many scenarios. We present an alternative method of scaling model Exponentially increasing processor power has fueled fifty years of continuous improvement in weather and climate
Colorado at Boulder, University of
BIT Numerical Mathematics 2005, to appear
. AKRIVIS AND S. LARSSON This system models the electric heating of a conductive body, with u beingDEPENDENT JOULE HEATING PROBLEM # G. AKRIVIS 1+ and S. LARSSON 2 June 27, 2005 1 Computer Science firstname.lastname@example.com Abstract. Completely discrete numerical methods for a nonlinear ellipticparabolic system, the time
Larsson, Stig
the electric heating of a conductive body, with u being the temperature, the electric potential-DEPENDENT JOULE HEATING PROBLEM G. AKRIVIS1 and S. LARSSON2 June 27, 2005 1 Computer Science email@example.com Abstract. Completely discrete numerical methods for a nonlinear elliptic-parabolic system, the time
On Aerodynamic Design Through Multipoint Numerical Optimization
On Aerodynamic Design Through Multipoint Numerical Optimization Howard P. Buckley, and David WT6, Canada A multipoint optimization approach is used to solve aerodynamic design problems en integral as an objective function is demonstrated to improve aerodynamic performance over a range of on
Comparing Aerodynamic Models for Numerical Simulation of
Comparing Aerodynamic Models for Numerical Simulation of Dynamics and Control of Aircraft and simulation of aircraft, yet other aerodynamics models exist that can provide more accurate results for certain simulations without a large increase in computational time. In this paper, sev- eral aerodynamics
Peraire, Jaime
by heating Ceramic block Glass workpiece before heating Figure 1: Thermal Replication (after Smith et al. [14NUMERICAL DESIGN TOOLS FOR THERMAL REPLICATION OF OPTICALQUALITY SURFACES Y.M. Stokes 1 Department. #12; Thermal replication of optical surfaces 1 1 Introduction Thermal replication is an industrial
Stokes, Yvonne
Numerical Analysis by L. Ridgway Scott
-dimensional spaces, i.e. square matrices. This chapter is mainly about matrix norms, their relationshipReview of Numerical Analysis by L. Ridgway Scott Princeton University Press, 2011 David S. Watkins as a mathematical subject. Scott wishes to reverse this trend, so he has written a text in which the mathematics
Random Matrix Theory, Numerical Computation and Applications
a precursor to a book on Random Matrix Theory that will be forthcoming. We reserve the right to reuseRandom Matrix Theory, Numerical Computation and Applications Alan Edelman, Brian D. Sutton new approaches to theory. We illustrate by describing such random matrix techniques as the stochastic
Edelman, Alan
Parallel Worldline Numerics: Implementation and Error Analysis
We give an overview of the worldline numerics technique, and discuss the parallel CUDA implementation of a worldline numerics algorithm. In the worldline numerics technique, we wish to generate an ensemble of representative closed-loop particle trajectories, and use these to compute an approximate average value for Wilson loops. We show how this can be done with a specific emphasis on cylindrically symmetric magnetic fields. The fine-grained, massive parallelism provided by the GPU architecture results in considerable speedup in computing Wilson loop averages. Furthermore, we give a brief overview of uncertainty analysis in the worldline numerics method. There are uncertainties from discretizing each loop, and from using a statistical ensemble of representative loops. The former can be minimized so that the latter dominates. However, determining the statistical uncertainties is complicated by two subtleties. Firstly, the distributions generated by the worldline ensembles are highly non-Gaussian, and so the standard error in the mean is not a good measure of the statistical uncertainty. Secondly, because the same ensemble of worldlines is used to compute the Wilson loops at different values of $T$ and $x_\\mathrm{ cm}$, the uncertainties associated with each computed value of the integrand are strongly correlated. We recommend a form of jackknife analysis which deals with both of these problems.
Dan Mazur; Jeremy S. Heyl
Numerical Linear Algebra and Optimization on Facebook
Numerical Linear Algebra and Optimization on Facebook: "In a relationship" or "It's complicated indicate their "relationship status" on Facebook (!!!). (I don't need to define Facebook, right?) #12;Facebook is close to ubiquitous (1.2 billion users in March 2014). 556 million people log on to Facebook
Sidorov, Nikita
NUMERICAL APPROACHES TO THERMALLY COUPLED PERFECT PLASTICITY S�OREN BARTELS AND TOM�AS ROUB at small strains exhibiting also stress-driven Prandtl-Reuss perfect plasticity are considered limit passage. Keywords: Thermodynamics, Prandtl-Reuss plasticity, Kelvin-Voigt rheology, thermal
Bartels, Soeren
ENERGY CONSERVATION BY ST ¨ ORMERTYPE NUMERICAL INTEGRATORS ERNST HAIRER \\Lambda AND CHRISTIAN the longtime energy conservation by a class of methods which integrate the linear part of oscillatory. HAIRER AND CH. LUBICH 2. Energy conservation by symmetric methods used with small step size
Plasma astrophysics in numerical simulations A. Nordlund
in periodic domains are discussed. Modifications of the solutions by numerically motivated alterations-excited dynamos were invented by the Danish inventor Søren Hjorth, who received the patent for this discovery industrial importance in producing the most powerful generators at the time, for which he, in turn, received
Brandenburg, Axel
Lecture Notes Numerical Methods for Incompressible Flow
1.2 The Gas Turbine Combustion Chamber . . . . . . . . . . . . . . . . . . 7 2 Derivation of models to fluid mechanics, including the derivation and discussion of the most important models and equations. The numerical methods discussed in the subsequent part are ordered due to the model they solve, i.e., we start
Soatto, Stefano
Microsoft Word - ICPRP-Appendices_corrmade.doc
hot-spot volcano Basalt; hawaiite at beginning of eruption Eruption along a 2-km fis- sure that produced a cone 100 m high, 0.3 km 3 lava, 0.25 km 3 tephra, max. Eruption...
computer languages. INIS descriptors; DATA ACQUISITION; DATA PROCESSING; G CODES; NUMERICAL DATA UDC 681
Applications of Numerical Optimal Control to Nonlinear Hybrid Systems
, Eaton Corporation, Eden Prairie, Minnesota Abstract This paper develops a technique for numerically
Zefran, Milo�?
PLUTO: a Numerical Code for Computational Astrophysics
We present a new numerical code, PLUTO, for the solution of hypersonic flows in 1, 2 and 3 spatial dimensions and different systems of coordinates. The code provides a multi-physics, multi-algorithm modular environment particularly oriented towards the treatment of astrophysical flows in presence of discontinuities. Different hydrodynamic modules and algorithms may be independently selected to properly describe Newtonian, relativistic, MHD or relativistic MHD fluids. The modular structure exploits a general framework for integrating a system of conservation laws, built on modern Godunov-type shock-capturing schemes. Although a plethora of numerical methods has been successfully developed over the past two decades, the vast majority shares a common discretization recipe, involving three general steps: a piecewise polynomial reconstruction followed by the solution of Riemann problems at zone interfaces and a final evolution stage. We have checked and validated the code against several benchmarks available in literature. Test problems in 1, 2 and 3 dimensions are discussed.
A. Mignone; G. Bodo; S. Massaglia; T. Matsakos; O. Tesileanu; C. Zanni; A. Ferrari
Advanced Numerical Weather Prediction Techniques for Solar Irradiance Forecasting : : Statistical, Data-Assimilation, and Ensemble Forecasting
Multiscale numerical weather prediction model. Progress inassimilating numerical weather prediction model for solarwith numerical weather prediction models. In: Solar Energy
Mathiesen, Patrick James
A numerical approach to beam deflections
to Beam Deflections. (December 1974) Joel Ernest Lovell, B. S. , Texas ASM University Chairman of Advisory Committee: Prof. E. S. Holdr edge A numerical method for determining deflections of canti- lever and simply supported beams directly from their M/EI... de- flections due to shear is solvable analytically for many types of loadings and cross sectional variations. However, the solution of the equation is cumbersome for any situations where M/EI cannot be described as a continuous function of length...
Lovell, Joel Ernest
Numerical Simulations of Hyperfine Transitions of Antihydrogen
One of the ASACUSA (Atomic Spectroscopy And Collisions Using Slow Antiprotons) collaboration's goals is the measurement of the ground state hyperfine transition frequency in antihydrogen, the antimatter counterpart of one of the best known systems in physics. This high precision experiment yields a sensitive test of the fundamental symmetry of CPT. Numerical simulations of hyperfine transitions of antihydrogen atoms have been performed providing information on the required antihydrogen events and the achievable precision.
Kolbinger, B; Diermaier, M; Lehner, S; Malbrunot, C; Massiczek, O; Sauerzopf, C; Simon, M C; Widmann, E
Exploring New Physics Frontiers Through Numerical Relativity
The demand to obtain answers to highly complex problems within strong-field gravity has been met with significant progress in the numerical solution of Einstein's equations - along with some spectacular results - in various setups. We review techniques for solving Einstein's equations in generic spacetimes, focusing on fully nonlinear evolutions but also on how to benchmark those results with perturbative approaches. The results address problems in high-energy physics, holography, mathematical physics, fundamental physics, astrophysics and cosmology.
Vitor Cardoso; Leonardo Gualtieri; Carlos Herdeiro; Ulrich Sperhake
The numerical evaluation of many dimensional integrals
assistance, and to the members of my committee, Drs. H. A. Luther and R, R. Hocking. TABLE OF CONTENTS Chapter Page I. INTRODUCTION Il. EXISTING METHODS OF NUMERICAL INTEGRATION . 3 Classical Methods Mo nte Carlo Methods . III. AN ERROR ANALYSIS . 10... Derivation of a Criterion. Error Approximation Using the Classical Method . . Error Estimate Using the Crude Monte Carlo Method 10 12 16 IV, EXPERIMENTS WITH CRUDE MONTE CARLO AND THE CLASSICAL METHOD. 18 A Method for Finding Exact Solutions...
Rhodes, Duane Lewis
Numerical Analysis of some Generalized Casimir Pistons
The Casimir force due to a scalar field on a piston in a cylinder of radius $r$ with a spherical cap of radius $R>r$ is computed numerically in the world-line approach. A geometrical subtraction scheme gives the finite interaction energy that determines the Casimir force. The spectral function of convex domains is obtained from a probability measure on convex surfaces that is induced by the Wiener measure on Brownian bridges the convex surfaces are the hulls of. The vacuum force on the piston by a scalar field satisfying Dirichlet boundary conditions is attractive in these geometries, but the strength and short-distance behavior of the force depends crucially on the shape of the piston casing. For a cylindrical casing with a hemispherical head, the force for $a/R\\sim 0$ does not depend on the dimension of the casing and numerically approaches $\\sim - 0.00326(4)\\hbar c/a^2$. Semiclassically this asymptotic force is due to short, closed and non-periodic trajectories that reflect once off the piston near its periphery. The semiclassical estimate $-\\hbar c/(96\\pi a^2)(1+2\\sqrt{R^2-r^2}/a)$ for the force when $a/r\\ll r/R\\leq 1$ reproduces the numerical results within statistical errors.
Martin Schaden
Numerical Tests of the Improved Fermilab Action
Recently, the Fermilab heavy-quark action was extended to include dimension-six and -seven operators in order to reduce the discretization errors. In this talk, we present results of the first numerical simulations with this action (the OK action), where we study the masses of the quarkonium and heavy-light systems. We calculate combinations of masses designed to test improvement and compare results obtained with the OK action to their counterparts obtained with the clover action. Our preliminary results show a clear improvement.
C. Detar; A. S. Kronfeld; M. B. Oktay
Issues in Numerical Simulation of Fire Suppression
This paper outlines general physical and computational issues associated with performing numerical simulation of fire suppression. Fire suppression encompasses a broad range of chemistry and physics over a large range of time and length scales. The authors discuss the dominant physical/chemical processes important to fire suppression that must be captured by a fire suppression model to be of engineering usefulness. First-principles solutions are not possible due to computational limitations, even with the new generation of tera-flop computers. A basic strategy combining computational fluid dynamics (CFD) simulation techniques with sub-grid model approximations for processes that have length scales unresolvable by gridding is presented.
Tieszen, S.R.; Lopez, A.R.
Numerical Simulations Unravel the Cosmic Web
The universe is permeated by a network of filaments, sheets, and knots collectively forming a "cosmic web.'' The discovery of the cosmic web, especially through its signature of absorption of light from distant sources by neutral hydrogen in the intergalactic medium, exemplifies the interplay between theory and experiment that drives science, and is one of the great examples in which numerical simulations have played a key and decisive role. We recount the milestones in our understanding of cosmic structure, summarize its impact on astronomy, cosmology, and physics, and look ahead by outlining the challenges faced as we prepare to probe the cosmic web at new wavelengths.
C. -A. Faucher-Giguere; A. Lidz; L. Hernquist
Template:IsNumeric | Open Energy Information
AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data CenterFranconia, Virginia:FAQ < RAPID Jump to:Seadov Pty LtdSteen,Ltd Jump to:TaosISGANAttribution Jump to: navigation, searchIsNumeric
GRChombo : Numerical Relativity with Adaptive Mesh Refinement
Numerical relativity has undergone a revolution in the past decade. With a well-understood mathematical formalism, and full control over the gauge modes, it is now entering an era in which the science can be properly explored. In this work, we introduce GRChombo, a new numerical relativity code written to take full advantage of modern parallel computing techniques. GRChombo's features include full adaptive mesh refinement with block structured Berger-Rigoutsos grid generation which supports non-trivial "many-boxes-in-many-boxes" meshing hierarchies, and massive parallelism through the Message Passing Interface (MPI). GRChombo evolves the Einstein equation with the standard BSSN formalism, with an option to turn on CCZ4 constraint damping if required. We show that GRChombo passes all the standard "Apples-to-Apples" code comparison tests. We also show that it can stably and accurately evolve vacuum black hole spacetimes such as binary black hole mergers, and non-vacuum spacetimes such as scalar collapses into black holes. As an illustration of its AMR capability, we demonstrate the evolution of triple black hole merger, which can be set up trivially in GRChombo.
Katy Clough; Pau Figueras; Hal Finkel; Markus Kunesch; Eugene A. Lim; Saran Tunyasuvunakool
Numerical relativity has undergone a revolution in the past decade. With a well-understood mathematical formalism, and full control over the gauge modes, it is now entering an era in which the science can be properly explored. In this work, we introduce GRChombo, a new numerical relativity code written to take full advantage of modern parallel computing techniques. GRChombo's features include full adaptive mesh refinement with block structured Berger-Rigoutsos grid generation which supports non-trivial "many-boxes-in-many-boxes" meshing hierarchies, and massive parallelism through the Message Passing Interface (MPI). GRChombo evolves the Einstein equation with the standard BSSN formalism, with an option to turn on CCZ4 constraint damping if required. We show that GRChombo passes all the standard "Apples-to-Apples" code comparison tests. We also show that it can stably and accurately evolve vacuum black hole spacetimes such as binary black hole mergers, and non-vacuum spacetimes such as scalar collapses into b...
Clough, Katy; Finkel, Hal; Kunesch, Markus; Lim, Eugene A; Tunyasuvunakool, Saran
Numerical simulation of detonation failure in nitromethane
Detonation failure in the homogeneous liquid explosive nitromethane has been observed experimentally in a wide variety of confining geometries. However, numerical simulation of these failure situations with a wave propagation code has been essentially non-existent due to the large differences between the critical diameter and the length of the reaction zone - characteristic dimensions which differ by about two orders of magnitude. This inability to spatially resolve both the reaction zone and geometries of significant size has led us to propose a new numerical technique, based on the stability criterion for rate-type material models, in which only temporal resolution of the reaction zone is required. Using an improved model for nitromethane, we have carried out a series of two-dimensional calculations which illustrate the utility of the present approach in predicting a wide range of experimental observations. Of particular computational significance is the removal of the difficulty requiring spatial resolution of the reaction zone, so that problems of practical size can be analyzed with existing computer capabilities.
Kipp, M.E.; Nunziato, J.W.
RELAP-7 Numerical Stabilization: Entropy Viscosity Method
The RELAP-7 code is the next generation nuclear reactor system safety analysis code being developed at the Idaho National Laboratory (INL). The code is based on the INL's modern scientific software development framework, MOOSE (Multi-Physics Object Oriented Simulation Environment). The overall design goal of RELAP-7 is to take advantage of the previous thirty years of advancements in computer architecture, software design, numerical integration methods, and physical models. The end result will be a reactor systems analysis capability that retains and improves upon RELAP5's capability and extends the analysis capability for all reactor system simulation scenarios. RELAP-7 utilizes a single phase and a novel seven-equation two-phase flow models as described in the RELAP-7 Theory Manual (INL/EXT-14-31366). The basic equation systems are hyperbolic, which generally require some type of stabilization (or artificial viscosity) to capture nonlinear discontinuities and to suppress advection-caused oscillations. This report documents one of the available options for this stabilization in RELAP-7 -- a new and novel approach known as the entropy viscosity method. Because the code is an ongoing development effort in which the physical sub models, numerics, and coding are evolving, so too must the specific details of the entropy viscosity stabilization method. Here the fundamentals of the method in their current state are presented.
R. A. Berry; M. O. Delchini; J. Ragusa
Numerical Results for the Blue Phases
We review recent numerical work investigating the equilibrium phase diagram, and the dynamics, of the cholesteric blue phases. In equilibrium numerical results confirm the predictions of the classic analytical theories, and extend them to incorporate different values of the elastic constants, or the effects of an applied electric field. There is a striking increase in the stability of blue phase I in systems where the cholesteric undergoes helical sense inversion, and the anomalous electrostriction observed in this phase is reproduced. Solving the equations of motion allows us to present results for the phase transition kinetics of blue phase I under dielectric or flexoelectric coupling to an applied electric field. We also present simulations of the blue phases in a flow field, showing how the disclination network acts to oppose the flow. The results are based on the Landau-de Gennes exapnsion of the liquid crystal free energy: that such a simple and elegant theory can predict such complex and subtle physical behaviour is remarkable.
G. P. Alexander; J. M. Yeomans
NEW NUMERICAL TECHNOLOGIES FOR THE SIMULATION OF ARC WELDING PROCESSES Michel Bellet 1 , Makhlouf Antipolis, France; firstname.lastname@example.com Keywords: welding, finite elements, material deposit, adaptive for arc welding simulation and analysis. The new numerical technologies essentially consist first
Paris-Sud XI, Université de
Numerical simulations for nodal domains and spectral minimal partitions
unpublished results of [HHO2] with efficient numerical computations. This is the main goal of this paper
Vial, Gr�gory
Heft 87 Olaf Cirpka CONTRACT: A Numerical Tool
Heft 87 Olaf Cirpka CONTRACT: A Numerical Tool for Contaminant Transport and Chemical Transformations #12;Heft 87 CONTRACT: A Numerical Tool for Contaminant Transport and Chemical Transformations der Universität Stuttgart #12;CONTRACT: A Numerical Tool for Contaminant Transport and Chemical
Cirpka, Olaf Arie
Data Assimilation for Idealised Mathematical Models of Numerical Weather Prediction
Data Assimilation for Idealised Mathematical Models of Numerical Weather Prediction Supervisors). Background: Numerical Weather Prediction (NWP) has seen significant gains in accuracy in recent years due is directed at achieving real-world impact in numerical weather prediction by addressing fundamental issues
Wirosoetisno, Djoko
Forecasting wave height probabilities with numerical weather prediction models
Forecasting wave height probabilities with numerical weather prediction models Mark S. Roulstona; Numerical weather prediction 1. Introduction Wave forecasting is now an integral part of operational weather methods for generating such forecasts from numerical model output from the European Centre for Medium
Stevenson, Paul
On the Numerical Simulation of Waterflooding of Heterogeneous
On the Numerical Simulation of Waterflooding of Heterogeneous Petroleum Reservoirs Jim Douglas, Jr displacement in petroleum reservoirs. A very detailed description of the numerical method is presented. Follow, 22290 Rio de Janeiro, RJ, Brazil #12; On the Numerical Simulation of Waterflooding of Heterogeneous
Douglas Jr., Jim
Numerical Simulation of Unsteady Three-Dimensional Sheet Cavitation
Numerical Simulation of Unsteady Three-Dimensional Sheet Cavitation A.H. Koop #12;Numerical Simulation of Unsteady Three-Dimensional Sheet Cavitation A.H. Koop Thesis University of Twente, Enschede, the Netherlands #12;NUMERICAL SIMULATION OF UNSTEADY THREE-DIMENSIONAL SHEET CAVITATION PROEFSCHRIFT ter
Twente, Universiteit
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2011; 86:15441557
. Frey1,2 1UPMC Univ Paris 06, UMR 7598, Laboratoire J.L. Lions, F-75005 Paris, France 2Universidad de present a three-dimensional numerical simulation on a mechanical device. Copyright 2011 John Wiley & Sons.L. Lions, F-75005 Paris, France. E-mail: firstname.lastname@example.com Copyright 2011 John Wiley & Sons, Ltd. #12
Frey, Pascal
Numerical modeling was used to simulate the constrained sintering process of the cathode contact layer during assembly of solid oxide fuel cells (SOFCs). A finite element model based on the continuum theory for sintering of porous bodies was developed and used to investigate candidate low-temperature cathode contact materials. Constitutive parameters for various contact materials under investigation were estimated from dilatometry screening tests, and the influence of processing time, processing temperature, initial grain size, and applied compressive stress on the free sintering response was predicted for selected candidate materials. The densification behavior and generated stresses within a 5-cell planar SOFC stack during sintering, high temperature operation, and room temperature shutdown were predicted. Insufficient constrained densification was observed in the stack at the proposed heat treatment, but beneficial effects of reduced grain size, compressive stack preload, and reduced thermal expansion coefficient on the contact layer densification and stresses were observed.
Koeppel, Brian J.; Liu, Wenning N.; Stephens, Elizabeth V.; Khaleel, Mohammad A.
Numerical Investigation of Josephson Junction Structures
Multilayered long Josephson Junction Structures form an interesting physical system where both nonlinearity and interaction between subsystems play an important role. Such systems allow to study physical effects that do not occur in single Josephson junction.The Sakai-Bodin-Pedersen model--a system of perturbed sine-Gordon equations--is used to study the dynamic states of stacks of inductively coupled long Josephson Junctions (LJJs). The corresponding static problem is numerically investigated as well. In order to study the stability of possible static solutions a Sturm-Liouville problem is generated and solved.The transitions from static to dynamic state and the scenario of these transitions are analyzed depending on the model parameters. Different physical characteristics--current-voltage characteristics, individual instant voltages and internal magnetic fields, are calculated and interpreted.
Hristov, I.; Dimova, S.; Boyadjiev, T. [Faculty of Mathematics and Informatics, Sofia University 5 James Bourchier Blvd., 1164 Sofia (Bulgaria)
Direct numerical simulation of turbulent reacting flows
The development of turbulent combustion models that reflect some of the most important characteristics of turbulent reacting flows requires knowledge about the behavior of key quantities in well defined combustion regimes. In turbulent flames, the coupling between the turbulence and the chemistry is so strong in certain regimes that is is very difficult to isolate the role played by one individual phenomenon. Direct numerical simulation (DNS) is an extremely useful tool to study in detail the turbulence-chemistry interactions in certain well defined regimes. Globally, non-premixed flames are controlled by two limiting cases: the fast chemistry limit, where the turbulent fluctuations. In between these two limits, finite-rate chemical effects are important and the turbulence interacts strongly with the chemical processes. This regime is important because industrial burners operate in regimes in which, locally the flame undergoes extinction, or is at least in some nonequilibrium condition. Furthermore, these nonequilibrium conditions strongly influence the production of pollutants. To quantify the finite-rate chemistry effect, direct numerical simulations are performed to study the interaction between an initially laminar non-premixed flame and a three-dimensional field of homogeneous isotropic decaying turbulence. Emphasis is placed on the dynamics of extinction and on transient effects on the fine scale mixing process. Differential molecular diffusion among species is also examined with this approach, both for nonreacting and reacting situations. To address the problem of large-scale mixing and to examine the effects of mean shear, efforts are underway to perform large eddy simulations of round three-dimensional jets.
Chen, J.H. [Sandia National Laboratories, Livermore, CA (United States)
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING Int. J. Numer. Meth. Biomed. Engng. 2014; 30:232248
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING Int. J. Numer. Meth. Biomed. International Journal for Numerical Methods in Biomedical Engineering published by John Wiley & Sons, Ltd in Biomedical Engineering published by John Wiley & Sons, Ltd. #12;SEGMENTATION OF BIOMEDICAL IMAGES 233
Grant, P. W.
Collisionless microinstabilities in stellarators. II. Numerical simulations
Microinstabilities exhibit a rich variety of behavior in stellarators due to the many degrees of freedom in the magnetic geometry. It has recently been found that certain stellarators (quasi-isodynamic ones with maximum-J geometry) are partly resilient to trapped-particle instabilities, because fast-bouncing particles tend to extract energy from these modes near marginal stability. In reality, stellarators are never perfectly quasi-isodynamic, and the question thus arises whether they still benefit from enhanced stability. Here, the stability properties of Wendelstein 7-X and a more quasi-isodynamic configuration, QIPC, are investigated numerically and compared with the National Compact Stellarator Experiment and the DIII-D tokamak. In gyrokinetic simulations, performed with the gyrokinetic code GENE in the electrostatic and collisionless approximation, ion-temperature-gradient modes, trapped-electron modes, and mixed-type instabilities are studied. Wendelstein 7-X and QIPC exhibit significantly reduced growth rates for all simulations that include kinetic electrons, and the latter are indeed found to be stabilizing in the energy budget. These results suggest that imperfectly optimized stellarators can retain most of the stabilizing properties predicted for perfect maximum-J configurations.
Proll, J. H. E.; Xanthopoulos, P.; Helander, P. [Max-Planck-Institut f�r Plasmaphysik, EURATOM Association, Teilinstitut Greifswald, Wendelsteinstra�e 1, 17491 Greifswald, Germany and Max-Planck/Princeton Research Center for Plasma Physics, 17491 Greifswald (Germany)] [Max-Planck-Institut f�r Plasmaphysik, EURATOM Association, Teilinstitut Greifswald, Wendelsteinstra�e 1, 17491 Greifswald, Germany and Max-Planck/Princeton Research Center for Plasma Physics, 17491 Greifswald (Germany)
Lava Creek Geothermal Area | Open Energy Information
AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home5b9fcbce19 No revision hasInformation Earth's Heat JumpInc Place:Keystone CleanLaton, California: Energy
Evaluation of numerical weather prediction for intra-day solar forecasting in the continental United States
for�numerical�weather�prediction�and�climate�models. ��Abstract:� Numerical�weather�prediction�(NWP)�models�are�Model�output�statistics�(MOS),�Numerical�Weather�Prediction�(
Mathiesen, Patrick; Kleissl, Jan
A high-resolution, cloud-assimilating numerical weather prediction model for solar irradiance forecasting
Multiscale�Numerical�Weather�Prediction�Model. ��Progress�assimilating�numerical�weather�prediction�model�for�solar�customizable� numerical�weather�prediction�model�that�is�
Mathiesen, Patrick; Collier, Craig; Kleissl, Jan
Notes 09. Numerical evaluation of natural modes and frequencies
. From Eq. (3) ?=K? M? (6) MEEN 617 ? HD#9. Numerical methods for finding eigenvalues & eigenvectors L. San Andr?s ? 2008 4 One can define the product 1? B = M K (6a) and write Eq...MEEN 617 ? HD#9. Numerical methods for finding eigenvalues & eigenvectors L. San Andr?s ? 2008 1 ME617 - Handout 9 Solving the eigenvalue problem - Numerical Evaluation of Natural Modes and Frequencies in MDOF systems The standard...
San Andres, Luis
Status of NINJA: the Numerical INJection Analysis project
The 2008 NRDA conference introduced the Numerical INJection Analysis project (NINJA), a new collaborative effort between the numerical relativity community and the data analysis community. NINJA focuses on modeling and searching for gravitational wave signatures from the coalescence of binary system of compact objects. We review the scope of this collaboration and the components of the first NINJA project, where numerical relativity groups shared waveforms and data analysis teams applied various techniques to detect them when embedded in colored Gaussian noise.
Benjamin Aylott; John G. Baker; William D. Boggs; Michael Boyle; Patrick R. Brady; Duncan A. Brown; Bernd Br�gmann; Luisa T. Buchman; Alessandra Buonanno; Laura Cadonati; Jordan Camp; Manuela Campanelli; Joan Centrella; Shourov Chatterjis; Nelson Christensen; Tony Chu; Peter Diener; Nils Dorband; Zachariah B. Etienne; Joshua Faber; Stephen Fairhurst; Benjamin Farr; Sebastian Fischetti; Gianluca Guidi; Lisa M. Goggin; Mark Hannam; Frank Herrmann; Ian Hinder; Sascha Husa; Vicky Kalogera; Drew Keppel; Lawrence E. Kidder; Bernard J. Kelly; Badri Krishnan; Pablo Laguna; Carlos O. Lousto; Ilya Mandel; Pedro Marronetti; Richard Matzner; Sean T. McWilliams; Keith D. Matthews; R. Adam Mercer; Satyanarayan R. P. Mohapatra; Abdul H. Mrou�; Hiroyuki Nakano; Evan Ochsner; Yi Pan; Larne Pekowsky; Harald P. Pfeiffer; Denis Pollney; Frans Pretorius; Vivien Raymond; Christian Reisswig; Luciano Rezzolla; Oliver Rinne; Craig Robinson; Christian R�ver; Luc�a Santamar�a; Bangalore Sathyaprakash; Mark A. Scheel; Erik Schnetter; Jennifer Seiler; Stuart L. Shapiro; Deirdre Shoemaker; Ulrich Sperhake; Alexander Stroeer; Riccardo Sturani; Wolfgang Tichy; Yuk Tung Liu; Marc van der Sluys; James R. van Meter; Ruslan Vaulin; Alberto Vecchio; John Veitch; Andrea Vicer�; John T. Whelan; Yosef Zlochower
A Numerical Upscaling Procedure to Estimate Effective Plane Wave ...
Jan 12, 2008 ... significant energy losses and velocity dispersion. To perform numerical ..... saturated transient flow in fractal soils. Set ? = ?L ? ?B ? ?R ? ?T ,�...
Numerical Modeling At Dixie Valley Geothermal Area (McKenna ...
McKenna & Blackwell, 2003) Jump to: navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Numerical Modeling At Dixie Valley Geothermal Area (McKenna &...
Energy stable schemes and numerical simulations of two phase ...
Title: Energy stable schemes and numerical simulations of two phase complex fluids by the phase-field method Abstact: We present an energetic variational�...
Numerical Modeling At Dixie Valley Geothermal Area (Iovenitti...
Iovenitti, Et Al., 2013) Jump to: navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Numerical Modeling At Dixie Valley Geothermal Area (Iovenitti, Et Al.,...
Direct Numerical Simulations and Robust Predictions of Cloud...
the center of the bubble cloud. Credit: Computational Science and Engineering Laboratory, ETH Zurich, Switzerland Direct Numerical Simulations and Robust Predictions of Cloud...
The Numerical Approximation of Solutions of ... - Purdue University
numerical simulation of waterflooding of heterogeneous petroleum reservoirs. Computational Geosciences, 1 (1997). 155�190. J. Douglas, Jr., C.-S. Huang, and�...
Numerical Data Advisory Board assembly of mathematial and physical sciences
The Numerical Data Advisory Board (NDAB) is an advisory body that provides expert overview, on a broad basis, of data needs and data programs as required for the advancement of science and technology. Board members, representing various disciplines, concern themselves with the quality, reliability, availability, accessibility, and dissemination of numerical data in physical, chemical, engineering, and interdisciplinary subjects as well as numeric and non-numeric data that arise in biology and geology. Topics of concern are addressed by the NDAB membership, or by specific, carefully chosen committees and panels established by NDAB in order to include experts appropriate to the subject at hand.
Not Available
Direct Numerical Simulations in Engine-like Geometries | Argonne...
Direct Numerical Simulations in Engine-like Geometries Event Sponsor: Mathematics and Computing ScienceArgonne Leadership Computing Facility Seminar Start Date: Nov 14 2014 -...
An Integrated Experimental and Numerical Study: Developing a...
and natural geothermal systems. The simulator can also be used to evaluate long-term CO2 sequestration potentials. State Minnesota Objectives Modify a numerical simulator (TOUGH2)...
Preliminary LSG Numerical Reliability DATt 5 January 1971
of the component application analysis and the FMEA. I . f.: I #12;� .... ATM-937 HY.MO. Preliminary l.SG Numerical
A numerical scheme for ionizing shock waves
A two-dimensional (2D) visual computer code to solve the steady state (SS) or transient shock problems including partially ionizing plasma is presented. Since the flows considered are hypersonic and the resulting temperatures are high, the plasma is partially ionized. Hence the plasma constituents are electrons, ions and neutral atoms. It is assumed that all the above species are in thermal equilibrium, namely, that they all have the same temperature. The ionization degree is calculated from Saha equation as a function of electron density and pressure by means of a nonlinear Newton type root finding algorithms. The code utilizes a wave model and numerical fluctuation distribution (FD) scheme that runs on structured or unstructured triangular meshes. This scheme is based on evaluating the mesh averaged fluctuations arising from a number of waves and distributing them to the nodes of these meshes in an upwind manner. The physical properties (directions, strengths, etc.) of these wave patterns are obtained by a new wave model: ION-A developed from the eigen-system of the flux Jacobian matrices. Since the equation of state (EOS) which is used to close up the conservation laws includes electronic effects, it is a nonlinear function and it must be inverted by iterations to determine the ionization degree as a function of density and temperature. For the time advancement, the scheme utilizes a multi-stage Runge-Kutta (RK) algorithm with time steps carefully evaluated from the maximum possible propagation speed in the solution domain. The code runs interactively with the user and allows to create different meshes to use different initial and boundary conditions and to see changes of desired physical quantities in the form of color and vector graphics. The details of the visual properties of the code has been published before (see [N. Aslan, A visual fluctuation splitting scheme for magneto-hydrodynamics with a new sonic fix and Euler limit, J. Comput. Phys. 197 (2004) 1-27]). The two-dimensional nature of ION-A was presented by a planar shock wave propagating over a circular obstacle. It was demonstrated that including the effects of ionization in calculating complex flows is important, even when they appear initially negligible. This code can be used to accurately simulate the nonlinear time dependent evolution of neutral or ionized plasma flows from supersonic to hypersonic regimes.
Aslan, Necdet [Yeditepe University, Physics Department, Kayisda g-circumflex i, 34755 Istanbul (Turkey)]. E-mail: email@example.com; Mond, Michael [Ben Gurion University, Mechanical Engineering Department, Beer Sheva (Israel)
Numerical simulation of transpiration cooling through porous , T. Gotzen1
55, 52056 Aachen SUMMARY Transpiration cooling using ceramic matrix composite (CMC) materials to facilitate such numerical simulations for a carbon/carbon material mounted in the side wall of a hot gasNumerical simulation of transpiration cooling through porous material W. Dahmen1 , T. Gotzen1 and S
Richard # 9501 EXPERIMENTAL AND NUMERICAL ANALYSIS OF CORROSION-
Richard # 9501 Page 1/9 EXPERIMENTAL AND NUMERICAL ANALYSIS OF CORROSION- INDUCED COVER CRACKING to present both experimental and numerical studies of the corrosion induced cracking pattern evolution of a reinforced concrete sample subjected to accelerated corrosion. The sample was a concrete beam (1000x250x100
Numerical Simulation of Wave Loads on Static Offshore Structures
Numerical Simulation of Wave Loads on Static Offshore Structures Hrvoje Jasak, Inno Gatin, Vuko Workshop, Cambridge, 30 July 2014 Numerical Simulation of Wave Loads on Static Offshore Structures p. #12 of Wave Loads on Static Offshore Structures p. #12;VOF Free Surface Flow Model Modelling of Free Surface
A new shock-capturing numerical scheme for ideal hydrodynamics
We present a new algorithm for solving ideal relativistic hydrodynamics based on Godunov method with an exact solution of Riemann problem for an arbitrary equation of state. Standard numerical tests are executed, such as the sound wave propagation and the shock tube problem. Low numerical viscosity and high precision are attained with proper discretization.
Zuzana Feckova; Boris Tomasik
A Numerical quadrature for the Schwarz-Chimera Method
A Numerical quadrature for the Schwarz-Chimera Method J. -B. Apoung Kamga1 and Olivier Pironneau2 1 Curie email@example.com Summary. Chimera [9] happens to be a version of Schwarz' method quadrature for the Schwarz-Chimera Method 3 3 Numerical Quadrature As such, the scheme is too costly
A Numerical Investigation of Wettability Alteration during Immiscible CO2
A Numerical Investigation of Wettability Alteration during Immiscible CO2 Flooding Process, April 2012 #12;2 Table of Contest Abstract 3 Introduction 3 Literature Review 5 CO2 Flooding 7 New alteration during CO2 flooding. However, limited research on numerical and/or analytical modeling
Hossain, M. Enamul
A Mobile Computing Architecture for Numerical Simulation Cyril Dumont
A Mobile Computing Architecture for Numerical Simulation Cyril Dumont Paris 12 University LACL of a numeric computation have to react when their context changes. This paper offers a new architecture, a mobile computing architecture, based on mobile agents and JavaSpace. At the end of this paper, we apply
CSU ATS703 Fall 2012 Numerical Weather Prediction
CSU ATS703 Fall 2012 Numerical Weather Prediction ATS703 is based on the course notes and papers method. A crucial element of accurate weather prediction is initialization, which is briefly discussed in Chapter 11. In the next decade, numerical weather prediction will expe- rience a revolution in model
KH Computational Physics-2006 Basic Numerical Algorithms Integration
KH Computational Physics- 2006 Basic Numerical Algorithms Integration Numerical integration competing factors one needs to consider � speed - number of function evaluations or grid points � precision or "smart" meshes with lower order routines Kristjan Haule, 2006 �1� #12;KH Computational Physics- 2006
Glashausser, Charles
Feckova, Zuzana
NUMERICAL INVESTIGATION OF CAVITATION IN MULTI-DIMENSIONAL COMPRESSIBLE FLOWS KRISTEN J. DEVAULT, cavitation, pseudospectral AMS subject classifications. 35Q30, 65M70, 76N99 1. Introduction. A long to first make a careful numerical study of the simplest possible scenario where one could expect cavitation
Numerical modelling of ice sheets, streams, and shelves
Numerical modelling of ice sheets, streams, and shelves Ed Bueler Karthaus, September 2012 models can actually help you, the reader, understand the behavior of partic- ular ice flows. Their most numerical codes produce numbers, but we want numbers that actually come from our continuum model. We analyse
Bueler, Ed
Numerical Simulation of a Natural Circulation Steam Generator
Numerical Simulation of a Natural Circulation Steam Generator W. Linzer \\Lambda , K. Ponweiser circulation steam generator. We focus on a model with a simple geometry consisting of two vertical pipes properties of water and steam. We present a numerical algorithm based on an explicit upwind discretization
Weinmüller, Ewa B.
Simple intrinsic defects in GaAs : numerical supplement.
This Report presents numerical tables summarizing properties of intrinsic defects in gallium arsenide, GaAs, as computed by density functional theory. This Report serves as a numerical supplement to the results published in: P.A. Schultz and O.A. von Lilienfeld, 'Simple intrinsic defects in GaAs', Modelling Simul. Mater. Sci Eng., Vol. 17, 084007 (2009), and intended for use as reference tables for a defect physics package in device models. The numerical results for density functional theory calculations of properties of simple intrinsic defects in gallium arsenide are presented.
Schultz, Peter Andrew
Grid-independent Issue in Numerical Heat Transfer
Grid independent is associated with the accuracy or even rationality of numerical results. This paper takes two-dimensional steady heat transfer for example to reveal the effect of grid resolution on numerical results. The law of grid dependence is obtained and a simple mathematical formula is presented. The production acquired here can be used as the guidance in choosing grid density in numerical simulation and get exact grid independent value without using infinite fine grid. Through analyzing grid independent, we can find the minimum number of grid cells that is needed to get grid-independent results. Such strategy can save computational resource while ensure a rational computational result.
Yao Wei; Wang Jian; Liao Guangxuan
Testing gravitational-wave searches with numerical relativity waveforms: Results from the first Numerical INJection Analysis (NINJA) project
The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational-wave data analysis communities. The purpose of NINJA is to study the sensitivity of existing gravitational-wave search algorithms using numerically generated waveforms and to foster closer collaboration between the numerical relativity and data analysis communities. We describe the results of the first NINJA analysis which focused on gravitational waveforms from binary black hole coalescence. Ten numerical relativity groups contributed numerical data which were used to generate a set of gravitational-wave signals. These signals were injected into a simulated data set, designed to mimic the response of the Initial LIGO and Virgo gravitational-wave detectors. Nine groups analysed this data using search and parameter-estimation pipelines. Matched filter algorithms, un-modelled-burst searches and Bayesian parameter-estimation and model-selection algorithms were applied to the data. We report the efficiency of these search methods in detecting the numerical waveforms and measuring their parameters. We describe preliminary comparisons between the different search methods and suggest improvements for future NINJA analyses.
Benjamin Aylott; John G. Baker; William D. Boggs; Michael Boyle; Patrick R. Brady; Duncan A. Brown; Bernd Br�gmann; Luisa T. Buchman; Alessandra Buonanno; Laura Cadonati; Jordan Camp; Manuela Campanelli; Joan Centrella; Shourov Chatterji; Nelson Christensen; Tony Chu; Peter Diener; Nils Dorband; Zachariah B. Etienne; Joshua Faber; Stephen Fairhurst; Benjamin Farr; Sebastian Fischetti; Gianluca Guidi; Lisa M. Goggin; Mark Hannam; Frank Herrmann; Ian Hinder; Sascha Husa; Vicky Kalogera; Drew Keppel; Lawrence E. Kidder; Bernard J. Kelly; Badri Krishnan; Pablo Laguna; Carlos O. Lousto; Ilya Mandel; Pedro Marronetti; Richard Matzner; Sean T. McWilliams; Keith D. Matthews; R. Adam Mercer; Satyanarayan R. P. Mohapatra; Abdul H. Mrou�; Hiroyuki Nakano; Evan Ochsner; Yi Pan; Larne Pekowsky; Harald P. Pfeiffer; Denis Pollney; Frans Pretorius; Vivien Raymond; Christian Reisswig; Luciano Rezzolla; Oliver Rinne; Craig Robinson; Christian R�ver; Luc�a Santamar�a; Bangalore Sathyaprakash; Mark A. Scheel; Erik Schnetter; Jennifer Seiler; Stuart L. Shapiro; Deirdre Shoemaker; Ulrich Sperhake; Alexander Stroeer; Riccardo Sturani; Wolfgang Tichy; Yuk Tung Liu; Marc van der Sluys; James R. van Meter; Ruslan Vaulin; Alberto Vecchio; John Veitch; Andrea Vicer�; John T. Whelan; Yosef Zlochower
Numerical Test of Born-Oppenheimer Approximation in Chaotic Systems
We study the validity of the Born-Oppenheimer approximation in chaotic dynamics. Using numerical solutions of autonomous Fermi accelerators, we show that the general adiabatic conditions can be interpreted as the narrowness of the chaotic region in phase space.
Jeong-Bo Shim; Mahir S. Hussein; Martina Hentschel
Numerical simulation of ceramic breeder pebble bed thermal creep behavior
Numerical simulation of ceramic breeder pebble bed thermal creep behavior Alice Ying *, Hulin Huang-310 825 1715/2599. E-mail address: firstname.lastname@example.com (A. Ying). 0022-3115/02/$ - see front matter �?
Abdou, Mohamed
Numerical Simulation Study on Transpired Solar Air Collector
The unglazed transpired solar air collector is now a well-recognized solar air heater for heating outside air directly. In this article, researchers introduced numerical simulation tools into the solar air collector research area, analyzed...
Wang, C.; Guan, Z.; Zhao, X.; Wang, D.
2.086 Numerical Computation for Mechanical Engineers, Spring 2012
This class introduces elementary programming concepts including variable types, data structures, and flow control. After an introduction to linear algebra and probability, it covers numerical methods relevant to mechanical ...
Patera, Anthony
Householder Symposium on Numerical Linear Algebra June 1721, 2002
for discussions. This year's symposium is held at Peebles Hotel Hydro in the small town of Peebles (populationHouseholder Symposium on Numerical Linear Algebra June 1721, 2002 Peebles Hotel Hydro, Scotland
Higham, Nicholas J.
problem in robotics, computer modeling and graphics, computational geometry and numerical optimization. We) the Research Council KU Leuven: (a) OT/10/038, (b) PF/10/002 Optimization in Engineering (OPTEC), (2
problem in robotics, computer modeling and graphics, computational geometry and numerical optimization. We) the Research Council KU Leuven: (a) project OT/10/038, (b) CoE EF/05/006 Optimization in Engineering (OPTEC
LONGTIME ENERGY CONSERVATION OF NUMERICAL METHODS FOR OSCILLATORY DIFFERENTIAL EQUATIONS ERNSTtime energy conservation, secondorder symmetric methods, frequency expansion, backward error analysis, Fermi
An investigation of the numerical treatment of condensation
The simulation of complete condensation continues to challenge the numerical methods currently used for multi-phase flow modeling; especially at low pressures, the change of phase process from a two-phase mixture to liquid ...
Sasson, Joseph
Experimental and numerical analysis of a deepwater mini-TLP
As the quest for oil and gas resources drives the industry to ever deeper waters, model testing still represents an essential step after numerical modeling when designing offshore platforms in these hostile environments. In an attempt to better...
Guichard, Aurelien
Numerical simulations of mine penetration in soft sediments
This study is concerned with developing a predictive model for depth of penetration of mines in soft sediments. The primary objective of this numerical study is to model soil resistance as a function of penetration depth using the finite element...
Yao, Zhigang
Numerical and Experimental Investigation of Internal Short Circuit...
Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site
in a Li-ion Cell Numerical and Experimental Investigation of Internal Short Circuit in a Li-ion Cell 2011 DOE Hydrogen and Fuel Cells Program, and Vehicle Technologies Program...
A numerical investigation of collisionality and turbulent transport
An investigation of collisionality's role in turbulent transport in magnetized plasma using the GS2 gyrokinetic simulation software is presented. The investigation consists of three parts, conducted by way of numerical ...
Lindsey, Martin L., S.B. Massachusetts Institute of Technology
Numerical Simulations of Gamma-Ray Burst Explosions
Gamma-ray bursts are a complex, non-linear system that evolves very rapidly through stages of vastly different conditions. They evolve from scales of few hundred kilometers where they are very dense and hot to cold and tenuous on scales of parsecs. As such, our understanding of such a phenomenon can truly increase by combining theoretical and numerical studies adopting different numerical techniques to face different problems and deal with diverse conditions. In this review, we will describe the tremendous advancement in our comprehension of the bursts phenomenology through numerical modeling. Though we will discuss studies mainly based on jet dynamics across the progenitor star and the interstellar medium, we will also touch upon other problems such as the jet launching, its acceleration, and the radiation mechanisms. Finally, we will describe how combining numerical results with observations from Swift and other instruments resulted in true understanding of the bursts phenomenon and the challenges still lyi...
Lazzati, Davide; L�pez-C�mara, Diego
Higher-Rank Numerical Ranges of Unitary and Normal Matrices
We verify a conjecture on the structure of higher-rank numerical ranges for a wide class of unitary and normal matrices. Using analytic and geometric techniques, we show precisely how the higher-rank numerical ranges for a generic unitary matrix are given by complex polygons determined by the spectral structure of the matrix. We discuss applications of the results to quantum error correction, specifically to the problem of identification and construction of codes for binary unitary noise models.
Man-Duen Choi; John A. Holbrook; David W. Kribs; Karol Zyczkowski
Numerical simulation of imbibition process in fractured cores
of the equations describing the imbibition of pure and carbonated water and the countercurrent flow of oil in porous rocks. Comparison between numerical solutions and experimental observations indicates that the models do simulate water imbibition and the COp.... ~~ +Blair 1964 presented a numerical solution for equations describing oil displacement by countercurrent water imbibition. Calculations were made for water imbibition into both linear and radial systems. Blair studied the effects on imbibition rates as a...
Jiang, Chunhuan
Numerical simulation of flow separation control by oscillatory fluid injection
NUMERICAL SIMULATION OF FLOW SEPARATION CONTROL BY OSCILLATORY FLUID INJECTION A Dissertation by CELERINO RESENDIZ ROSAS Submitted to the O?ce of Graduate Studies of Texas A&M University in partial ful?llment of the requirements for the degree... of DOCTOR OF PHILOSOPHY May 2005 Major Subject: Aerospace Engineering NUMERICAL SIMULATION OF FLOW SEPARATION CONTROL BY OSCILLATORY FLUID INJECTION A Dissertation by CELERINO RESENDIZ ROSAS Submitted to Texas A&M University in partial ful...
Resendiz Rosas, Celerino
Numerical simulation is vital to analyse complex chemi-cal or physical processes. The Institute for Numerical
, says: "The power supplies and the cooling systems could not cope with the constant heavy loads of data computer networks supplied by specialist organisations have performed the institute's numerical simulation
Griebel, Michael
Numerical investigation of transient hydrothermal processes around intrusions: heat-transfer and fluid-
Numerical investigation of transient hydrothermal processes around intrusions: heat the intrusion. Keywords: Hydrothermal processes, numerical modelling, magmatic intrusion, permeability- depth around magmatic intrusions have been obtained through coupled hydrothermal numerical modelling that takes
Multigrid methods for improving the variational data assimilation in numerical weather prediction
conditions are needed to solve numerical weather prediction models: initial condition and boundary conditionMultigrid methods for improving the variational data assimilation in numerical weather prediction: numerical weather prediction, variational data assimilation, minimization procedure, multigrid methods, cell
Kwak, Do Young
Chimera: A hybrid approach to numerical loop quantum cosmology
The existence of a quantum bounce in isotropic spacetimes is a key result in loop quantum cosmology (LQC), which has been demonstrated to arise in all the models studied so far. In most of the models, the bounce has been studied using numerical simulations involving states which are sharply peaked and which bounce at volumes much larger than the Planck volume. An important issue is to confirm the existence of the bounce for states which have a wide spread, or which bounce closer to the Planck volume. Numerical simulations with such states demand large computational domains, making them very expensive and practically infeasible with the techniques which have been implemented so far. To overcome these difficulties, we present an efficient hybrid numerical scheme using the property that at the small spacetime curvature, the quantum Hamiltonian constraint in LQC, which is a difference equation with uniform discretization in volume, can be approximated by a Wheeler-DeWitt differential equation. By carefully choosi...
Diener, Peter; Singh, Parampreet
Accelerating numerical solution of Stochastic Differential Equations with CUDA
Numerical integration of stochastic differential equations is commonly used in many branches of science. In this paper we present how to accelerate this kind of numerical calculations with popular NVIDIA Graphics Processing Units using the CUDA programming environment. We address general aspects of numerical programming on stream processors and illustrate them by two examples: the noisy phase dynamics in a Josephson junction and the noisy Kuramoto model. In presented cases the measured speedup can be as high as 675x compared to a standard CPU, which corresponds to sev eral billion integration steps per second. This means that calculations which took weeks can now be completed in less than one hour. This brings stochastic simulation to a completely new level, opening for research a whole new range of problems which can now be solved interactively.
Januszewski, M
Numerical simulation of water flow around a rigid fishing net
This paper is devoted to the simulation of the flow around and inside a rigid axisymmetric net. We describe first how experimental data have been obtained. We show in detail the modelization. The model is based on a Reynolds Averaged Navier-Stokes turbulence model penalized by a term based on the Brinkman law. At the out-boundary of the computational box, we have used a "ghost" boundary condition. We show that the corresponding variational problem has a solution. Then the numerical scheme is given and the paper finishes with numerical simulations compared with the experimental data.
Roger Lewandowski; G�raldine Pichot
Numerical Analysis of Spherically Convergent Rayleigh-Taylor Experiments
In the frame of a CEA/US DOE collaboration, radiation driven spherically convergent experiments were performed on the Nova laser in order m measure the Rayleigh-Taylor growth at the ablation front. Numerical simulations using the 2D Lagrangian code FCI2 have correctly reproduced experiments in moderate convergent geometry. [C. Cherfils et al., PRL 83, 5507 (1999)]. Experiments have addressed convergence ratios up to 4 by considering larger capsules, larger hohlraum and longer laser pulses [S.G. Glendinning et al., to be published in Physics of Plasmas]. Numerical analysis of these high convergence implosions is presented, and the effect of convergence on the Rayleigh-Taylor growth is investigated.
Galmiche, D.; Cherfils, C.; Glendinning, S.G.; Remington, B.A.; Richard, A.
NCPART: management of ICEMDDN output for numerical control users
NCPART is a procedure developed by the Numerical Control Department at Bendix Kansas City Division to handle the entry to and exit from ICEMDDN, and process all of the local files output by ICEMDDN. The NCPART procedure is menu driven, and provides automatic access to ICEMDDN and any files necessary to process information with ICEM for numerical Control users. Basically, the procedure handles all of the ICEMDDN operations that involve operating system commands, and frees the NC programmer to concentrate on his/her work as a programmer.
Rossini, B.F.
Geodesic Motion on Closed Spaces: Two Numerical Examples
The geodesic structure is very closely related to the trace of the Laplace operator, involved in the calculation of the expectation value of the energy momentum tensor in Universes with non trivial topology. The purpose of this work is to provide concrete numerical examples of geodesic flows. Two manifolds with genus $g=0$ are given. In one the chaotic regions, form sets of negligible or zero measure. In the second example the geodesic flow, shows the presence of measurable chaotic regions. The approach is "experimental", numerical, and there is no attempt to an analytical calculation.
Daniel Muller
Numerical Methodology to Evaluate Fast Reactor Sodium Combustion
In the present study, a numerical methodology for sodium combustion has been developed for the safety evaluation of a liquid-metal-cooled fast reactor. The methodology includes a fast-running zone model computer program for safety evaluation, a field model program for multidimensional thermal hydraulics, and a chemical reaction analysis program based on chemical equilibrium theory. Two recently performed experiments have been analyzed using the computer programs, and the numerical results are in good agreement with the experiments. Although sodium combustion is a complex phenomenon, use of these computer programs gives better understanding of the coupled thermal hydraulics and chemical reaction.
Yamaguchi, Akira; Takata, Takashi; Okano, Yasushi [Japan Nuclear Cycle Development Institute (Japan)
Numerical calculation of wave refraction by digital computer
+~ with time. 20 22 22 25 28 WiL), H, H s o vs. F(i ), DL L 0 43 10 Refraction pattern for uniform slope 10 second p~riod 60 approach; numerical analysis solution 0 Refraction pattern for uniform slope 10 second period 120 approach; numerical... of varying radii (3) as shown in Figure l. It can be seen that the crest deforms and turns toward the segment of lowest celerity. Let M and N be two adjacent points along the crest separated by a distance An at time t (Figure 2a). The corresponding wave...
Orr, Terry Edwin
10.34 Numerical Methods Applied to Chemical Engineering, Fall 2001
Numerical methods for solving problems arising in heat and mass transfer, fluid mechanics, chemical reaction engineering, and molecular simulation. Topics: numerical linear algebra, solution of nonlinear algebraic equations ...
Beers, Kenneth J.
Probabilistic methods applied to 2D electromagnetic numerical dosimetry
the number of realizations. When the modeling of the 3D interaction between the human body and a cellphoneProbabilistic methods applied to 2D electromagnetic numerical dosimetry D. Voyer F. Musy L. Nicolas dosimetry problems in order to take into account the variability of the input parameters. Methodology
Numerical Prediction of High-Impact Local Weather: A
Chapter 6 Numerical Prediction of High-Impact Local Weather: A Driver for Petascale Computing Ming winds, lightning, hurricanes and winter storms, cause hundreds of deaths and average annual economic of mitigating the impacts of such events on the economy and society is obvious, our ability to do so
Xue, Ming
Department of Numerical Analysis Modeling the Austenite Ferrite
Department of Numerical Analysis Modeling the Austenite Ferrite Transformation by Cellular Ferrite Transformation by Cellular Automaton Improving Interface Stability Master of Science Thesis. Computational Materials Science 48.3 (2010): 692-699] for the austenite to ferrite transformation in low
Vuik, Kees
110112 NatSciLab -Numerical Software Introduction to MATLAB
Outline 110112 NatSciLab - Numerical Software Introduction to MATLAB Onur Oktay Jacobs University Bremen Spring 2010 #12;MATLAB Desktop Environment The Command line Operators MATLAB data types Outline 1 MATLAB Desktop Environment 2 The Command line A quick start Indexing 3 Operators Arithmetic (Array
Vasylkevych, Sergiy
Simple intrinsic defects in InAs : numerical predictions.
This Report presents numerical tables summarizing properties of intrinsic defects in indium arsenide, InAs, as computed by density functional theory using semi-local density functionals, intended for use as reference tables for a defect physics package in device models.
Inertia ratchets: A numerical study versus theory B. Lindner,1
Inertia ratchets: A numerical study versus theory B. Lindner,1 L. Schimansky-Geier,1 P. Reimann,2 P in a periodic but asymmetric potential ratchet is investigated. Our special focus is on the influence of inertia cur- rent in periodic structures that lack reflection symmetry ratchets , e.g. see the reviews 1
Lindner, Benjamin
On the numerical solution of some problems of environmental pollution
On the numerical solution of some problems of environmental pollution Quang A Dang1 , Matthias world and their importance will even increase in the future. High pollution of air, water and soil may to the mathematical formula- tion of air pollution models. Let G be a cylindrical domain in the three
Ehrhardt, Matthias
Numerical Modeling of Brain Dynamics in Traumatic Situations -Impulsive Translations
.S.A. Abstract We numerically model the brain dy- namics during and after impulsive head translations using brain injuries appear among boxers and shaken babies despite minimal rotations of their heads. Modeling head translations also helps understand the brain dynamics during head rotations about an arbitrary
Burtscher, Martin
Numerical Methods for the Valuation of Synthetic Collateralized Debt Obligations
Numerical Methods for the Valuation of Synthetic Collateralized Debt Obligations by Xiaofang Ma Methods for the Valuation of Synthetic Collateralized Debt Obligations Xiaofang Ma Doctor of Philosophy methods for synthetic CDO valuation are presented. iii #12;Acknowledgements Although many people have
Toronto, University of
Numerical modeling of magnetohydrodynamic activity in the Swarthmore Spheromak Experiment
Numerical modeling of magnetohydrodynamic activity in the Swarthmore Spheromak Experiment V. S resistive magnetohydrodynamic MHD simulation are compared to experimental data from the Swarthmore Spheromak is shown to reproduce global equilibrium magnetic field profiles of the spheromaks as well as much
Brown, Michael R.
RisR1448(EN) Numerical Study of the
is the better understanding of the physics of the aeroelastic motion of wind turbine blades in order to improve of wind turbine blades in order to improve the numerical simulation of such dynamical systems. In previous works, both aerodynamic damping cal- culations and fully-coupled aeroelastic simulations of wind turbine
Theoretical and Numerical Analysis of Polarization for Time Dependent Radiative
transport equation with respect to the polariza- tion parameters solve the matrix-valued radiative transferTheoretical and Numerical Analysis of Polarization for Time Dependent Radiative Transfer firstname.lastname@example.com Abstract We consider the matrix-valued radiative transfer equations for the Stokes param- eters
Bal, Guillaume
Numerical Aerodynamic Optimization Incorporating Laminar-Turbulent Transition Prediction
Numerical Aerodynamic Optimization Incorporating Laminar-Turbulent Transition Prediction J. Driver-dimensional NewtonKrylov aerodynamic shape optimization algorithm is applied to several optimization problems a striking demonstration of the capability of the Newton Krylov aerodynamic optimization algorithm to design
Numerical simulation of the truss spar 'Horn Mountain' using COUPLE
using an integrated marine monitoring system attached to the truss spar. In this study, dynamic analysis of the truss spar interacting with its mooring and riser system was performed using a time-domain numerical code, known as �?COUPLE�?. The simulated...
Theckum Purath, Basil
A Generic Framework for Interprocedural Analysis of Numerical Properties
A Generic Framework for Interprocedural Analysis of Numerical Properties Markus M¨uller-Olm1, Germany firstname.lastname@example.com 2 TU M¨unchen, Institut f¨ur Informatik, I2 80333 M of the abstracted collecting semantics is also used by M¨uller-Olm and Seidl in [12] where (in absence of equality
Müller-Olm, Markus
A Unified Framework for Numerically Inverting Laplace Joseph Abate
A Unified Framework for Numerically Inverting Laplace Transforms Joseph Abate 900 Hammond Road been many applications; e.g., see the survey by Abate et al. (1999) and the textbook treatment by Kao transforms; e.g., see the surveys in Abate and Whitt (1992) and Chapter 19 of Davies (2002), the extensive
Whitt, Ward
The existence of a quantum bounce in isotropic spacetimes is a key result in loop quantum cosmology (LQC), which has been demonstrated to arise in all the models studied so far. In most of the models, the bounce has been studied using numerical simulations involving states which are sharply peaked and which bounce at volumes much larger than the Planck volume. An important issue is to confirm the existence of the bounce for states which have a wide spread, or which bounce closer to the Planck volume. Numerical simulations with such states demand large computational domains, making them very expensive and practically infeasible with the techniques which have been implemented so far. To overcome these difficulties, we present an efficient hybrid numerical scheme using the property that at the small spacetime curvature, the quantum Hamiltonian constraint in LQC, which is a difference equation with uniform discretization in volume, can be approximated by a Wheeler-DeWitt differential equation. By carefully choosing a hybrid spatial grid allowing the use of partial differential equations at large volumes, and with a simple change of geometrical coordinate, we obtain a surprising reduction in the computational cost. This scheme enables us to explore regimes which were so far unachievable for the isotropic model in LQC. Our approach also promises to significantly reduce the computational cost for numerical simulations in anisotropic LQC using high performance computing.
Peter Diener; Brajesh Gupt; Parampreet Singh
A probabilistic numerical method for optimal multiple switching problem
performing the numerical resolution, the storage of the Monte Carlo simulation paths is not needed. Then, we at the minimum cost under a given constraint on the loss of load probability or on the level of energy non to thank Thomas Vareschi, Xavier Warin and the participants of the FiME seminar and the Energy Finance
Numerical simulation of fluid flow in porous/fractured media
Theoretical models of fluid flow in porous/fractured media can help in the design of in situ fossil energy and mineral extraction technologies. Because of the complexity of these processes, numerical solutions are usually required. Sample calculations illustrate the capabilities of present day computer models.
Travis B.J.; Cook, T.L.
Computing Flowpipe of Nonlinear Hybrid Systems with Numerical Methods
(which embed information on the plant dynamics) are more and more used: for such systems, the controller simulation used nowadays are very powerful and efficient to approximate the behavior of complex dynamical guaranteed generalization of classical efficient numerical simulation methods, including with variable
Cytosolic Ca2+ ions regulate numerous aspects of cellular
Cytosolic Ca2+ ions regulate numerous aspects of cellular activity in virtually all cell types (Berridge et al. 2000). This versatility is made possible by the diverse mechanisms by which Ca2+ signals; Marchant & Parker, 2000). Endogenous Ca2+ binding proteins play a key role in determining the magnitude
Parker, Ian
Numerical Optimization Methodology for the Design of Power Equipment
health effects from overhead power line conductors [6,7]. Even though no conclusive evidence has been of utility companies have been involved in a redesign of their overhead transmission lines to avoid potentialNumerical Optimization Methodology for the Design of Power Equipment Gaurav Tewari* My M. Hua
Mamishev, Alexander
Friction versus dilation revisited: Insights from theoretical and numerical models
Friction versus dilation revisited: Insights from theoretical and numerical models N. Makedonska,1 controlled by the frictional strength of the fault gouge, a granular layer that accumulates between the fault friction coefficient) of such granular layers is the systems resistance to dilation, a byprocess
Einat, Aharonov
NUMERICAL PREDICTION OF TEMPERATURE DISTRIBUTION IN PEM FUEL CELLS S. Shimpalee and S. Dutta distribution inside a straight channel proton exchange membrane ( PEM) fuel cell and the effect of heat is called the proton exchange membrane (PEM) fuel cell that operates at a signi¢cantly lower temperature
Van Zee, John W.
Numerical Model of a Tensioner System and Flex Joint
Top Tensioned Riser (TTR) and Steel Catenary Riser (SCR) are often used in a floating oil/gas production system deployed in deep water for oil transport. This study focuses on the improvements to the existing numerical code, known as CABLE3D...
Huang, Han
Numerical simulation of turbulent jet primary breakup in Diesel engines
Numerical simulation of turbulent jet primary breakup in Diesel engines Peng Zeng1 Marcus Herrmann" IRMA Strasbourg, 23.Jan.2008 #12;Introduction DNS of Primary Breakup in Diesel Injection Phase Transition Modeling Turbulence Modeling Summary Outline 1 Introduction 2 DNS of Primary Breakup in Diesel
Helluy, Philippe
Numerical algorithms for the computation of the Smith normal
Numerical algorithms for the computation of the Smith normal form of integral matrices C of the Smith normal form of integral matrices are described. More specifically, the com pound matrix method of the algorithms. AMS Subject Classification: Primary 65F30, Secondary 15A21, 15A36. Key words and phrases: Smith
Seberry, Jennifer
On the Numerical Simulation of Waterflooding of Heterogeneous Petroleum Reservoirs Jim Douglas, Jr Simulation of Waterflooding of Heterogeneous Petroleum Reservoirs · A two-stage operator-splitting allows consider as a model problem the two-phase immiscible displacement in petroleum reservoirs. A very detailed
Furtado, Fred
Cavitation level-acoustic intensity hysteresis: experimental and numerical characterization
Cavitation level-acoustic intensity hysteresis: experimental and numerical characterization P such as sonoporation, inertial cavitation is commonly considered as the main candidate inducing membrane poration. Thus, characterizing inertial cavitation, as related to bubble size distribution and medium history, is of great
Numerical studies of the metamodel fitting and validation processes
algorithms and application to a nuclear safety computer code show the relevance of this new sequential this problem consists in replacing cpu time expensive computer models by cpu inexpensive mathematical functions to fit the metamodel) has to provide adequate space filling properties. We adopt a numerical approach
Efficient Numerical Simulation for Long Range Wave Propagation 1
Efficient Numerical Simulation for Long Range Wave Propagation 1 Kai Huang 2 George Papanicolaou 3 for simulating wave propagation over long dis- tances with both weak and strong scatterers. In domains with weak heterogeneities the wave field is decomposed into forward propagating and back scattered modes using two coupled
Solna, Knut
Mercury's thermo-chemical evolution from numerical models constrained
Mercury's thermo-chemical evolution from numerical models constrained by MESSENGER observations Globe de Paris, France #12;Basics facts about Mercury · Semi-major axis: 0.39 AU · 3:2 spin Earth!) · Black body temperature: 440 K #12;Exploration of Mercury Mariner10 ·First spacecraft to use
Cerveny, Vlastislav
Numerical modeling on behavior under loads of self compacting concrete floors with hollow space 134 15, 2011 1 Table of Contents 1. David Procházka, Jirí Brozovský, Klára Krízová High-Strength Concrete-Ionut Stefancu Dynamic nonlinear analysis of structural reinforced concrete walls energy dissipators with shear
Fernandez, Thomas
155: Numerical Models of Groundwater Flow and Transport
155: Numerical Models of Groundwater Flow and Transport EKKEHARD HOLZBECHER1 AND SHAUL SOREK2 1. #12;2402 GROUNDWATER Calibration as a task cannot be separated from the other tasks. Inverse modeling of the Negev, J. Blaustein Institutes for Desert Research, Sede Boker, Israel The article gives an introduction
Sorek, Shaul
RisR1252(EN) The Numerical Wind Atlas
potential from high quality wind measurements. It estimates the local influences on the wind by small hillsRisøR1252(EN) The Numerical Wind Atlas -- the KAMM/WAsP Method Helmut P. Frank, Ole Rathmann The method of combining the Karlsruhe Atmospheric Mesoscale Model, KAMM, with the Wind Atlas Analysis
Numerical Analysis of Non-Uniformities and Anisotropy in High-Temperature Superconducting Coils
, critical current density (superconductivity), high-temperature superconductors, numerical analysis, superconducting coils, transport ac loss....
Hu, Di; Ainslie, Mark D.; Zou, Jin; Cardwell, David A.
Numerical simulations of welds of thick steel pieces of interest for the thermonuclear fusion ITER machine
Numerical simulations of welds of thick steel pieces of interest for the thermonuclear fusion ITER machine
Carmignani, B
On Numerical Considerations for Modeling Reactive Astrophysical Shocks
Simulating detonations in astrophysical environments is often complicated by numerical approximations to shock structure. A common prescription to ensure correct detonation speeds (and associated quantities) is to prohibit burning inside the numerically broadened shock (Fryxell et al. 1989). We have performed a series of simulations to verify the efficacy of this approximation and to understand how resolution and dimensionality might affect its use. Our results show that, in one dimension, prohibiting burning in the shock is important wherever the carbon burning length is not resolved, in keeping with the results of Fryxell et al. (1989). In two dimensions, we find that the prohibition of shock burning effectively inhibits the development of cellular structure for all but the most highly-resolved cases. We discuss the possible impacts this outcome may have on sub-grid models and detonation propagation in Type Ia supernovae.
Papatheodore, Thomas L [ORNL] [ORNL; Messer, Bronson [ORNL] [ORNL
Hydrodynamics of Hypersonic Jets: Experiments and Numerical Simulations
Stars form in regions of the galaxy that are denser and cooler than the mean interstellar medium. These regions are called Giant Molecular Clouds. At the beginning of their life, up to $10^5-10^6$ years, stars accrete matter from their rich surrounding environment and are origin of a peculiar phenomenon that is the jet emission. Jets from Young Stellar Objects (YSOs) are intensively studied by the astrophysical community by observations at different wavelengths, analytical and numerical modeling and laboratory experiments. Indications about the jet propagation and its resulting morphologies are here obtained by means of a combined study of hypersonic jets carried out both in the laboratory and by numerical simulations.
Belan, Marco; Tordella, Daniela; Massaglia, Silvano; Ferrari, Attilio; Mignone, Andrea; Bodenschatz, Eberhard
Numerical Relativity as a tool for studying the Early Universe
Numerical simulations are becoming a more effective tool for conducting detailed investigations into the evolution of our universe. In this article, we show how the framework of numerical relativity can be used for studying cosmological models. The author is working to develop a large-scale simulation of the dynamical processes in the early universe. These take into account interactions of dark matter, scalar perturbations, gravitational waves, magnetic fields and a turbulent plasma. The code described in this report is a GRMHD code based on the Cactus framework and is structured to utilize one of several different differencing methods chosen at run-time. It is being developed and tested on the University of Houston's Maxwell cluster.
David Garrison
Numerical Modeling of Charged Black Holes with Massive Dilaton
In this paper the static, spherically symmetric and electrically charged black hole solutions in Einstein-Born-Infeld gravity with massive dilaton are investigated numerically. The Continuous Analog of Newton Method (CANM) is used to solve the corresponding nonlinear multipoint boundary value problems (BVPs). The linearized BVPs are solved numerically by means of collocation scheme of fourth order. A special class of solutions are the extremal ones. We show that the extremal horizons within the framework of the model satisfy some nonlinear system of algebraic equations. Depending on the charge $q$ and dilaton mass $\\gamma$, the black holes can have no more than three horizons. This allows us to construct some Hermite polynomial of third order. Its real roots describe the number, the type and other characteristics of the horizons.
T. L. Boyadjiev; P. P. Fiziev
Fractional Calculus in Hydrologic Modeling: A Numerical Perspective
Fractional derivatives can be viewed either as a handy extension of classical calculus or, more deeply, as mathematical operators defined by natural phenomena. This follows the view that the diffusion equation is defined as the governing equation of a Brownian motion. In this paper, we emphasize that fractional derivatives come from the governing equations of stable Levy motion, and that fractional integration is the corresponding inverse operator. Fractional integration, and its multi-dimensional extensions derived in this way, are intimately tied to fractional Brownian (and Levy) motions and noises. By following these general principles, we discuss the Eulerian and Lagrangian numerical solutions to fractional partial differential equations, and Eulerian methods for stochastic integrals. These numerical approximations illuminate the essential nature of the fractional calculus.
David A. Benson; Mark M. Meerschaert; Jordan Revielle
Numerical investigation of recirculation in the UTSI MHD combustor
Numerical studies were carried out to investigate the gross structure of flow in cylindrical combustors. The combustor configurations studied are variations of a working design used at the University of Tennessee Space Institute to burn pulverized coal at temperatures in excess of 3000K for generation of a plasma feeding a magnetohydrodynamic channel. The numerical studies were conducted for an isothermal fluid; the main objective of the calculations was to study the effect of the oxidant injection pattern on the gross structure of recirculating flows within the combustor. The calculations illustrate the basic features of the flow in combustors of this type and suggest implications for the injection of coal and oxidizer in this type of combustor.
Schulz, R.J.; Lee, J.J.; Giel, T.V. Jr.
Accurate complex scaling of three dimensional numerical potentials
The complex scaling method, which consists in continuing spatial coordinates into the complex plane, is a well-established method that allows to compute resonant eigenfunctions of the time-independent Schroedinger operator. Whenever it is desirable to apply the complex scaling to investigate resonances in physical systems defined on numerical discrete grids, the most direct approach relies on the application of a similarity transformation to the original, unscaled Hamiltonian. We show that such an approach can be conveniently implemented in the Daubechies wavelet basis set, featuring a very promising level of generality, high accuracy, and no need for artificial convergence parameters. Complex scaling of three dimensional numerical potentials can be efficiently and accurately performed. By carrying out an illustrative resonant state computation in the case of a one-dimensional model potential, we then show that our wavelet-based approach may disclose new exciting opportunities in the field of computational non-Hermitian quantum mechanics.
Cerioni, Alessandro [European Synchrotron Radiation Facility, 6 rue Horowitz, BP220 38043 Grenoble Cedex 9 (France); Genovese, Luigi; Duchemin, Ivan; Deutsch, Thierry [Laboratoire de simulation atomistique (L-Sim), SP2M, UMR-E CEA/UJF-Grenoble 1, INAC, Grenoble F-38054 (France)
Numerical analysis of nanostructures for enhanced light extraction from OLEDs
Nanostructures, like periodic arrays of scatters or low-index gratings, are used to improve the light outcoupling from organic light-emitting diodes (OLED). In order to optimize geometrical and material properties of such structures, simulations of the outcoupling process are very helpful. The finite element method is best suited for an accurate discretization of the geometry and the singular-like field profile within the structured layer and the emitting layer. However, a finite element simulation of the overall OLED stack is often beyond available computer resources. The main focus of this paper is the simulation of a single dipole source embedded into a twofold infinitely periodic OLED structure. To overcome the numerical burden we apply the Floquet transform, so that the computational domain reduces to the unit cell. The relevant outcoupling data are then gained by inverse Flouqet transforming. This step requires a careful numerical treatment as reported in this paper.
Zschiedrich, L; Burger, S; Schmidt, F; 10.1117/12.2001132
Regularization of spherical and axisymmetric evolution codes in numerical relativity
Several interesting astrophysical phenomena are symmetric with respect to the rotation axis, like the head-on collision of compact bodies, the collapse and/or accretion of fields with a large variety of geometries, or some forms of gravitational waves. Most current numerical relativity codes, however, can not take advantage of these symmetries due to the fact that singularities in the adapted coordinates, either at the origin or at the axis of symmetry, rapidly cause the simulation to crash. Because of this regularity problem it has become common practice to use full-blown Cartesian three-dimensional codes to simulate axi-symmetric systems. In this work we follow a recent idea idea of Rinne and Stewart and present a simple procedure to regularize the equations both in spherical and axi-symmetric spaces. We explicitly show the regularity of the evolution equations, describe the corresponding numerical code, and present several examples clearly showing the regularity of our evolutions.
Milton Ruiz; Miguel Alcubierre; Dario Nunez
Numerical method for shear bands in ductile metal with inclusions
A numerical method for mesoscale simulation of high strain-rate loading of ductile metal containing inclusions is described. Because of small-scale inhomogeneities, such a composite material is prone to localized shear deformation (adiabatic shear bands). The modeling framework is the Generalized Method of Cells of Paley and Aboudi [Mech. Materials, vol. 14, pp. /27-139, 1992], which ensures that the micromechanical response of the material is reflected in the behavior of the composite at the mesoscale. To calculate the effective plastic strain rate when shear bands are present, the analytic and numerical analysis of shear bands by Glimm, Plohr, and Sharp [Mech. Materials, vol. 24, pp. 31-41, 1996] is adapted and extended.
Plohr, Jee Yeon N [Los Alamos National Laboratory; Plohr, Bradley J [Los Alamos National Laboratory
Numerical simulations of gravitational collapse in Einstein-aether theory
We study gravitational collapse of a spherically symmetric scalar field in Einstein-aether theory (general relativity coupled to a dynamical unit timelike vector field). The initial value formulation is developed, and numerical simulations are performed. The collapse produces regular, stationary black holes, as long as the aether coupling constants are not too large. For larger couplings a finite area singularity occurs. These results are shown to be consistent with the stationary solutions found previously.
David Garfinkle; Christopher Eling; Ted Jacobson
A Generic Framework for Interprocedural Analysis of Numerical Properties Markus M�?ullerOlm 1, Germany firstname.lastname@example.com 2 TU M�?unchen, Institut f�?ur Informatik, I2 80333 M of the abstracted collecting semantics is also used by M�?ullerOlm and Seidl in [12] where (in absence of equality
Seidl, Helmut
Comparative study of variational chaos indicators and ODEs' numerical integrators
The reader can find in the literature a lot of different techniques to study the dynamics of a given system and also, many suitable numerical integrators to compute them. Notwithstanding the recent work of Maffione et al. (2011a) for mappings, a detailed comparison among the widespread indicators of chaos in a general system is still lacking. Such a comparison could lead to select the most efficient algorithms given a certain dynamical problem. Furthermore, in order to choose the appropriate numerical integrators to compute them, more comparative studies among numerical integrators are also needed. This work deals with both problems. We first extend the work of Maffione et al. (2011) for mappings to the 2D H\\'enon & Heiles (1964) potential, and compare several variational indicators of chaos: the Lyapunov Indicator (LI); the Mean Exponential Growth Factor of Nearby Orbits (MEGNO); the Smaller Alignment Index (SALI) and its generalized version, the Generalized Alignment Index (GALI); the Fast Lyapunov Indicator (FLI) and its variant, the Orthogonal Fast Lyapunov Indicator (OFLI); the Spectral Distance (D) and the Dynamical Spectras of Stretching Numbers (SSNs). We also include in the record the Relative Lyapunov Indicator (RLI), which is not a variational indicator as the others. Then, we test a numerical technique to integrate Ordinary Differential Equations (ODEs) based on the Taylor method implemented by Jorba & Zou (2005) (called taylor), and we compare its performance with other two well-known efficient integrators: the Prince & Dormand (1981) implementation of a Runge-Kutta of order 7-8 (DOPRI8) and a Bulirsch-St\\"oer implementation. These tests are run under two very different systems from the complexity of their equations point of view: a triaxial galactic potential model and a perturbed 3D quartic oscillator.
Luciano A. Darriba; Nicol�s P. Maffione; Pablo M. Cincotta; Claudia M. Giordano
Numerical procedure for calculating temperature profiles in LMFBR coolant channels
A new numerical procedure (which makes use of a weighted residuals procedure in space and a fully-implicit finite difference procedure in time), for calculating temperatures in an LMFBR coolant channel has been developed and incorporated into the Super System Code (SSC). This procedure is highly accurate on a nodal basis and has greatly increased computational efficiency as compared to the method formerly in SSC.
Horak, W.C.; Kennett, R.J.; Guppy, J.G.
Segmented vs conventional numerals: legibility and long term retention
the environmental chamber. Sub jects The subjects were thirty male students between the ages of 1g and 27 from the Industrial Engineering department at Texas A&M University. Subjects were divided into three groups of 10. Procedure Exposure time and number... December 1971 Ma]or Subject: Industrial Engineering SEGMENTED VS CONVENTIONAL NUMERALS: LEGIBILITY AND LONG TERM RETENTION A Thesis STEVE EDGAR HILL Approved as to style and content by: Elias Chairman of Committee) r. A. W. ortham (Head...
Hill, Steve Edgar
Numerical heat conduction in hydrodynamical models of colliding hypersonic flows
Hydrodynamical models of colliding hypersonic flows are presented which explore the dependence of the resulting dynamics and the characteristics of the derived X-ray emission on numerical conduction and viscosity. For the purpose of our investigation we present models of colliding flow with plane-parallel and cylindrical divergence. Numerical conduction causes erroneous heating of gas across the contact discontinuity which has implications for the rate at which the gas cools. We find that the dynamics of the shocked gas and the resulting X-ray emission are strongly dependent on the contrast in the density and temperature either side of the contact discontinuity, these effects being strongest where the postshock gas of one flow behaves quasi-adiabatically while the postshock gas of the other flow is strongly radiative. Introducing additional numerical viscosity into the simulations has the effect of damping the growth of instabilities, which in some cases act to increase the volume of shocked gas and can re-he...
Parkin, E R
Numerical Methods of Computational Electromagnetics for Complex Inhomogeneous Systems
Understanding electromagnetic phenomena is the key in many scientific investigation and engineering designs such as solar cell designs, studying biological ion channels for diseases, and creating clean fusion energies, among other things. The objectives of the project are to develop high order numerical methods to simulate evanescent electromagnetic waves occurring in plasmon solar cells and biological ion-channels, where local field enhancement within random media in the former and long range electrostatic interactions in the latter are of major challenges for accurate and efficient numerical computations. We have accomplished these objectives by developing high order numerical methods for solving Maxwell equations such as high order finite element basis for discontinuous Galerkin methods, well-conditioned Nedelec edge element method, divergence free finite element basis for MHD, and fast integral equation methods for layered media. These methods can be used to model the complex local field enhancement in plasmon solar cells. On the other hand, to treat long range electrostatic interaction in ion channels, we have developed image charge based method for a hybrid model in combining atomistic electrostatics and continuum Poisson-Boltzmann electrostatics. Such a hybrid model will speed up the molecular dynamics simulation of transport in biological ion-channels.
Cai, Wei
Transient productivity index for numerical well test simulations
The most difficult aspect of numerical simulation of well tests is the treatment of the Bottom Hole Flowing (BHF) Pressure. In full field simulations, this pressure is derived from the Well-block Pressure (WBP) using a numerical productivity index which accounts for the grid size and permeability, and for the well completion. This productivity index is calculated assuming a pseudo-steady state flow regime in the vicinity of the well and is therefore constant during the well production period. Such a pseudo-steady state assumption is no longer valid for the early time of a well test simulation as long as the pressure perturbation has not reached several grid-blocks around the well. This paper offers two different solutions to this problem: (1) The first one is based on the derivation of a Numerical Transient Productivity Index (NTPI) to be applied to Cartesian grids; (2) The second one is based on the use of a Corrected Transmissibility and Accumulation Term (CTAT) in the flow equation. The representation of the pressure behavior given by both solutions is far more accurate than the conventional one as shown by several validation examples which are presented in the following pages.
Blanc, G.; Ding, D.Y.; Ene, A. [Institut Francais du Petrole, Pau (France)] [and others
Numerical implication of Riemann problem theory for fluid dynamics
The Riemann problem plays an important role in understanding the wave structure of fluid flow. It is also crucial step in some numerical algorithms for accurately and efficiently computing fluid flow; Godunov method, random choice method, and from tracking method. The standard wave structure consists of shock and rarefaction waves. Due to physical effects such as phase transitions, which often are indistinguishable from numerical errors in an equation of state, anomalkous waves may occur, ''rarefaction shocks'', split waves, and composites. The anomalous waves may appear in numerical calculations as waves smeared out by either too much artificial viscosity or insufficient resolution. In addition, the equation of state may lead to instabilities of fluid flow. Since these anomalous effects due to the equation of state occur for the continuum equations, they can be expected to occur for all computational algorithms. The equation of state may be characterized by three dimensionless variables: the adiabatic exponent ..gamma.., the Grueneisen coefficient GAMMA, and the fundamental derivative G. The fluid flow anomalies occur when inequalities relating these variables are violated. 18 refs.
Menikoff, R.
Kinetic theory and numerical simulations of two-species coagulation
In this work we study the stochastic process of two-species coagulation. This process consists in the aggregation dynamics taking place in a ring. Particles and clusters of particles are set in this ring and they can move either clockwise or counterclockwise. They have a probability to aggregate forming larger clusters when they collide with another particle or cluster. We study the stochastic process both analytically and numerically. Analytically, we derive a kinetic theory which approximately describes the process dynamics. One of our strongest assumptions in this respect is the so called well-stirred limit, that allows neglecting the appearance of spatial coordinates in the theory, so this becomes effectively reduced to a zeroth dimensional model. We determine the long time behavior of such a model, making emphasis in one special case in which it displays self-similar solutions. In particular these calculations answer the question of how the system gets ordered, with all particles and clusters moving in the same direction, in the long time. We compare our analytical results with direct numerical simulations of the stochastic process and both corroborate its predictions and check its limitations. In particular, we numerically confirm the ordering dynamics predicted by the kinetic theory and explore properties of the realizations of the stochastic process which are not accessible to our theoretical approach.
Carlos Escudero; Fabricio Macia; Raul Toral; Juan J. L. Velazquez
Advanced numerical methods in mesh generation and mesh adaptation
Numerical solution of partial differential equations requires appropriate meshes, efficient solvers and robust and reliable error estimates. Generation of high-quality meshes for complex engineering models is a non-trivial task. This task is made more difficult when the mesh has to be adapted to a problem solution. This article is focused on a synergistic approach to the mesh generation and mesh adaptation, where best properties of various mesh generation methods are combined to build efficiently simplicial meshes. First, the advancing front technique (AFT) is combined with the incremental Delaunay triangulation (DT) to build an initial mesh. Second, the metric-based mesh adaptation (MBA) method is employed to improve quality of the generated mesh and/or to adapt it to a problem solution. We demonstrate with numerical experiments that combination of all three methods is required for robust meshing of complex engineering models. The key to successful mesh generation is the high-quality of the triangles in the initial front. We use a black-box technique to improve surface meshes exported from an unattainable CAD system. The initial surface mesh is refined into a shape-regular triangulation which approximates the boundary with the same accuracy as the CAD mesh. The DT method adds robustness to the AFT. The resulting mesh is topologically correct but may contain a few slivers. The MBA uses seven local operations to modify the mesh topology. It improves significantly the mesh quality. The MBA method is also used to adapt the mesh to a problem solution to minimize computational resources required for solving the problem. The MBA has a solid theoretical background. In the first two experiments, we consider the convection-diffusion and elasticity problems. We demonstrate the optimal reduction rate of the discretization error on a sequence of adaptive strongly anisotropic meshes. The key element of the MBA method is construction of a tensor metric from hierarchical edge-based error estimates. We conclude that the quasi-optimal mesh must be quasi-uniform in this metric. All numerical experiments are based on the publicly available Ani3D package, the collection of advanced numerical instruments.
Lipnikov, Konstantine [Los Alamos National Laboratory; Danilov, A [MOSCOW, RUSSIA; Vassilevski, Y [MOSCOW, RUSSIA; Agonzal, A [UNIV OF LYON
Numerical investigation of closed-loop control for Hall accelerators
Low frequency discharge current oscillations in Hall accelerators are conventionally damped with external inductor-capacitor (LC) or resistor-inductor-capacitor (RLC) networks. The role of such network in the stabilization of the plasma discharge is investigated with a numerical model and the potential advantages of proportional-integral-derivative (PID) closed-loop control over RLC networks are subsequently assessed using either discharge voltage or magnetic field modulation. Simulations confirm the reduction of current oscillations in the presence of a RLC network, but suggest that PID control could ensure nearly oscillation-free operation with little sensitivity toward the PID settings.
Barral, S.; Miedzik, J. [Institute of Plasma Physics and Laser Microfusion, 01497 Warsaw (Poland)
Experimental and Numerical Investigation of Flows in Expanding Channels
We present an experimental realization of the classical Jeffery-Hamel flows inside a wedge-shaped channel. We compare the measured velocity fields with the predictions of Jeffery-Hamel theory. A detailed experimental study of bifurcation diagrams for the solutions reveals the absolute stability of the pure outflow solution and an interesting hysteretic structure for bifurcations. We also observe a multiple vortex flow regime predicted earlier numerically and analytically. Experimental studies of the stability of the flow to perturbations at the channel exit are also conducted.
Vorobieff, Peter; Putkaradze, Vakhtang
During the year 2008, the United States National Weather Service (NWS) completed an eight fold increase in sampling capability for weather radars to 250 m resolution. This increase is expected to improve warning lead times by detecting small scale features sooner with increased reliability; however, current NWS operational model domains utilize grid spacing an order of magnitude larger than the radar data resolution, and therefore the added resolution of radar data is not fully exploited. The assimilation of radar reflectivity and velocity data into high resolution numerical weather model forecasts where grid spacing is comparable to the radar data resolution was investigated under a Laboratory Directed Research and Development (LDRD) 'quick hit' grant to determine the impact of improved data resolution on model predictions with specific initial proof of concept application to daily Savannah River Site operations and emergency response. Development of software to process NWS radar reflectivity and radial velocity data was undertaken for assimilation of observations into numerical models. Data values within the radar data volume undergo automated quality control (QC) analysis routines developed in support of this project to eliminate empty/missing data points, decrease anomalous propagation values, and determine error thresholds by utilizing the calculated variances among data values. The Weather Research and Forecasting model (WRF) three dimensional variational data assimilation package (WRF-3DVAR) was used to incorporate the QC'ed radar data into input and boundary conditions. The lack of observational data in the vicinity of SRS available to NWS operational models signifies an important data void where radar observations can provide significant input. These observations greatly enhance the knowledge of storm structures and the environmental conditions which influence their development. As the increase in computational power and availability has made higher resolution real-time model simulations possible, the need to obtain observations to both initialize numerical models and verify their output has become increasingly important. The assimilation of high resolution radar observations therefore provides a vital component in the development and utility of numerical model forecasts for both weather forecasting and contaminant transport, including future opportunities to improve wet deposition computations explicitly.
Chiswell, S.; Buckley, R.
Connecting Numerical Relativity and Data Analysis of Gravitational Wave Detectors
Gravitational waves deliver information in exquisite detail about astrophysical phenomena, among them the collision of two black holes, a system completely invisible to the eyes of electromagnetic telescopes. Models that predict gravitational wave signals from likely sources are crucial for the success of this endeavor. Modeling binary black hole sources of gravitational radiation requires solving the Eintein equations of General Relativity using powerful computer hardware and sophisticated numerical algorithms. This proceeding presents where we are in understanding ground-based gravitational waves resulting from the merger of black holes and the implications of these sources for the advent of gravitational-wave astronomy.
Shoemaker, Deirdre; London, Lionel; Pekowsky, Larne
Numerical studies of a simple Coulomb blockade model
simple model of the Coulomb blockade is studied. In this model, two interacting electrons tunnel in a one-dimensional structure with two barriers in series. The two-particle, time-dependent Schrodinger equation is solved numerically. It is found... tunneling by the other, and the inside electron tunnels out only after the outside electron tunnels a second time. The distribution of the charge density for the two electrons is shown in a series of three-dimensional figures as a function of time. Also...
Shao, Jianfeng
Diffusive mesh relaxation in ALE finite element numerical simulations
The theory for a diffusive mesh relaxation algorithm is developed for use in three-dimensional Arbitary Lagrange/Eulerian (ALE) finite element simulation techniques. This mesh relaxer is derived by a variational principle for an unstructured 3D grid using finite elements, and incorporates hourglass controls in the numerical implementation. The diffusive coefficients are based on the geometric properties of the existing mesh, and are chosen so as to allow for a smooth grid that retains the general shape of the original mesh. The diffusive mesh relaxation algorithm is then applied to an ALE code system, and results from several test cases are discussed.
Dube, E.I.
Property:Iso3166Numeric | Open Energy Information
AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data CenterFranconia, Virginia: Energy ResourcesLoadingPenobscot County,ContAddr2 Jump to: navigation,PVYearsDisplay/Graphics JumpNumeric Jump
Management of numerical control data from the ICEM DDN products
NCPART is a procedure developed by the Numerical Control Department at Bendix Kansas City Division to handle the entry to and exit from ICEM DDN, and process all of the local files output by ICEM DDN. The NCPART procedure is menu driven, and provides automatic access to ICEM DDN and any files necessary to process information with ICEM for Numerical Control users. Basically, the procedure handles all of the ICEM DDN operations that involve operating system commands, and frees the NC Programmer to concentrate on his/her work as a programmer. Most NC Programmers at BKC are well versed in the use of operating system commands for attaching and cataloging files, but the use of procedures makes the programmer's job much simpler, especially in the case of ICEM, where multiple files must be present in order to take advantage of the available features. Also, much is dependent on what an individual plans to do once in ICEM, because different applications require different files, and all of these files may not be needed at any given user session.
Comparing Post-Newtonian and Numerical-Relativity Precession Dynamics
Binary black-hole systems are expected to be important sources of gravitational waves for upcoming gravitational-wave detectors. If the spins are not colinear with each other or with the orbital angular momentum, these systems exhibit complicated precession dynamics that are imprinted on the gravitational waveform. We develop a new procedure to match the precession dynamics computed by post-Newtonian (PN) theory to those of numerical binary black-hole simulations in full general relativity. For numerical relativity NR) simulations lasting approximately two precession cycles, we find that the PN and NR predictions for the directions of the orbital angular momentum and the spins agree to better than $\\sim 1^{\\circ}$ with NR during the inspiral, increasing to $5^{\\circ}$ near merger. Nutation of the orbital plane on the orbital time-scale agrees well between NR and PN, whereas nutation of the spin direction shows qualitatively different behavior in PN and NR. We also examine how the PN equations for precession a...
Ossokine, Serguei; Kidder, Lawrence E; Pfeiffer, Harald P; Scheel, Mark A; Szil�gyi, B�la
Binary black-hole systems are expected to be important sources of gravitational waves for upcoming gravitational-wave detectors. If the spins are not colinear with each other or with the orbital angular momentum, these systems exhibit complicated precession dynamics that are imprinted on the gravitational waveform. We develop a new procedure to match the precession dynamics computed by post-Newtonian (PN) theory to those of numerical binary black-hole simulations in full general relativity. For numerical relativity NR) simulations lasting approximately two precession cycles, we find that the PN and NR predictions for the directions of the orbital angular momentum and the spins agree to better than $\\sim 1^{\\circ}$ with NR during the inspiral, increasing to $5^{\\circ}$ near merger. Nutation of the orbital plane on the orbital time-scale agrees well between NR and PN, whereas nutation of the spin direction shows qualitatively different behavior in PN and NR. We also examine how the PN equations for precession and orbital-phase evolution converge with PN order, and we quantify the impact of various choices for handling partially known PN terms.
Serguei Ossokine; Michael Boyle; Lawrence E. Kidder; Harald P. Pfeiffer; Mark A. Scheel; B�la Szil�gyi
Free and constrained symplectic integrators for numerical general relativity
We consider symplectic time integrators in numerical General Relativity and discuss both free and constrained evolution schemes. For free evolution of ADM-like equations we propose the use of the Stoermer-Verlet method, a standard symplectic integrator which here is explicit in the computationally expensive curvature terms. For the constrained evolution we give a formulation of the evolution equations that enforces the momentum constraints in a holonomically constrained Hamiltonian system and turns the Hamilton constraint function from a weak to a strong invariant of the system. This formulation permits the use of the constraint-preserving symplectic RATTLE integrator, a constrained version of the Stoermer-Verlet method. The behavior of the methods is illustrated on two effectively 1+1-dimensional versions of Einstein's equations, that allow to investigate a perturbed Minkowski problem and the Schwarzschild space-time. We compare symplectic and non-symplectic integrators for free evolution, showing very different numerical behavior for nearly-conserved quantities in the perturbed Minkowski problem. Further we compare free and constrained evolution, demonstrating in our examples that enforcing the momentum constraints can turn an unstable free evolution into a stable constrained evolution. This is demonstrated in the stabilization of a perturbed Minkowski problem with Dirac gauge, and in the suppression of the propagation of boundary instabilities into the interior of the domain in Schwarzschild space-time.
Ronny Richter; Christian Lubich
Numerical and experimental investigation of vortical flow-flame interaction
A massively parallel coupled Eulerian-Lagrangian low Mach number reacting flow code is developed and used to study the structure and dynamics of a forced planar buoyant jet flame in two dimensions. The numerical construction uses a finite difference scheme with adaptive mesh refinement for solving the scalar conservation equations, and the vortex method for the momentum equations, with the necessary coupling terms. The numerical model construction is presented, along with computational issues regarding the parallel implementation. An experimental acoustically forced planar jet burner apparatus is also developed and used to study the velocity and scalar fields in this flow, and to provide useful data for validation of the computed jet. Burner design and laser diagnostic details are discussed, along with the measured laboratory jet flame dynamics. The computed reacting jet flow is also presented, with focus on both large-scale outer buoyant structures and the lifted flame stabilization dynamics. A triple flame structure is observed at the flame base in the computed flow, as is theoretically expected, but was not observable with present diagnostic techniques in the laboratory flame. Computed and experimental results are compared, along with implications for model improvements.
Najm, H.N.; Schefer, R.W.; Milne, R.B.; Mueller, C.J. [Sandia National Labs., Livermore, CA (United States); Devine, K.D.; Kempka, S.N. [Sandia National Labs., Albuquerque, NM (United States)
Radiative feedback and cosmic molecular gas: numerical method
We present results from self-consistent 3D numerical simulations of cosmic structure formation with a multi-frequency radiative transfer scheme and non-equilibrium molecular chemistry of 13 primordial species (e-, H, H+, H-, He, He+, He++, H2, H2+, D, D+, HD, HeH+), performed by using the simulation code GADGET. We describe our implementation and show tests for ionized sphere expansion in a static and dynamic density field around a central radiative source, and for cosmological abundance evolution coupled with the cosmic microwave background radiation. As a demonstrative application of radiative feedback on molecular gas, we run also cosmological simulations of early structure formation in a ~1Mpc size box. Our tests agree well with analytical and numerical expectations. Consistently with other works, we find that ionization fronts from central sources can boost H2 fractions in shock-compressed gas. The tight dependence on H2 lead to a corresponding boost of HD fractions, as well. We see a strong lowering of the the typical molecular abundances up to several orders of magnitudes which partially hinders further gas collapse of pristine neutral gas, and clearly suggests the need of re-ionized gas or metal cooling for the formation of the following generation of structures.
Margarita Petkova; Umberto Maio
The relative strength between forward and reverse shock emission in early gamma-ray burst (GRB) afterglow reflects that of magnetic energy densities in the two shock regions. We numerically show that with the current standard treatment, the fireball magnetization is underestimated by up to two orders of magnitude. This discrepancy is especially large in the sub-relativistic reverse shock regime (i.e., the thin shell and intermediate regime), where most optical flashes were detected. We provide new analytic estimates of the reverse shock emission based on a better shock approximation, which well describe numerical results in the intermediate regime. We show that the reverse shock temperature at the onset of afterglow is constant, ( {Gamma}-bar{sub d}-1){approx}8 Multiplication-Sign 10{sup -2}, when the dimensionless parameter {xi}{sub 0} is more than several. Our approach is applied to case studies of GRB 990123 and 090102, and we find that magnetic fields in the fireballs are even stronger than previously believed. However, these events are still likely to be due to a baryonic jet with {sigma} {approx} 10{sup -3} for GRB 990123 and {approx}3 Multiplication-Sign 10{sup -4} to 3 for GRB 090102.
Harrison, Richard; Kobayashi, Shiho, E-mail: firstname.lastname@example.com [Astrophysics Research Institute, Liverpool John Moores University, Liverpool, L3 5RF (United Kingdom)
Applied Numerical Mathematics 58 (2008) 871883 www.elsevier.com/locate/apnum
Applied Numerical Mathematics 58 (2008) 871883 www.elsevier.com/locate/apnum Updating the QR / Applied Numerical Mathematics 58 (2008) 871883 preliminary computation of QR decompositions
Gutknecht, Martin H.
This dissertation mainly studied on numerical approximation methods as a solution of the integrability problem and the measure of welfare changes, and demonstrated how numerical algorithms can be applied in empirical studies ...
Lim, Sung Jin
E-Print Network 3.0 - asymptotic numerical method Sample Search...
numerical ... Source: Ecole Polytechnique, Centre de mathmatiques Collection: Mathematics 2 Introduction Uniform Estimates for Transmission Problems 3D Multiscaled...
E-Print Network 3.0 - adaptive numerical methods Sample Search...
- Fachbereich Mathematik und Informatik, Philipps-Universitt Marburg Collection: Mathematics 2 Introduction Adaptive wavelet numerical simulation Turbulence modelling Summary...
Numerical study of high heat ux pool boiling heat transfer Ying He a,*, Masahiro Shoji b
Numerical study of high heat ¯ux pool boiling heat transfer Ying He a,*, Masahiro Shoji b , Shigeo simulation model of boiling heat transfer is proposed based on a numerical macrolayer model [S. Maruyama, M. Shoji, S. Shimizu, A numerical simulation of transition boiling heat transfer, in: Proceedings
Maruyama, Shigeo
The Numerical Modelling Research and Development Division is responsible for research into and develop-
into and develop- ment of numerical weather prediction models and other meteorological applications, that are opera in the field of numerical weather prediction: atmospheric and oceanographic modelling, physical and statistical132 The Numerical Modelling Research and Development Division is responsible for research
Haak, Hein
24 More Years of Numerical Weather Prediction: A Model Performance Model
24 More Years of Numerical Weather Prediction: A Model Performance Model Gerard Cats May 26, 2008 Abstract For two formulations of currently usual numerical weather prediction models the evolution in such a model is much 1 #12;24 More Years of Numerical Weather Prediction Gerard Cats higher than in a sis
Evaluating the ability of a numerical weather prediction model to forecast tracer concentrations during ETEX 2
Evaluating the ability of a numerical weather prediction model to forecast tracer concentrations an operational numerical weather prediction model to forecast air quality are also investigated. These potential a numerical weather prediction (NWP) model independently of the CTM. The NWP output is typically archived
Dacre, Helen
Temporal Changes in Wind as Objects for Evaluating Mesoscale Numerical Weather Prediction
a method of evaluating numerical weather prediction models by comparing the characteristics of temporal for biases in features forecast by the model. 1. Introduction Verification of numerical weather predictionTemporal Changes in Wind as Objects for Evaluating Mesoscale Numerical Weather Prediction DARAN L
Knievel, Jason Clark
Numerical Methods for the Bogoliubov-Tolmachev-Shirkov model in superconductivity theory
In the work, the numerical methods are designed for the Bogoliubov-Tolmachev-Shirkov model in superconductivity theory. The numerical methods are novel and effective to determine the critical transition temperature and approximate to the energy gap function of the above model. Finally, a numerical example confirming the theoretical results is presented.
Zhihao Ge; Ruihua Li
Numerical Age Computation of the Antarctic Ice Sheet for Dating Deep Ice Cores
Numerical Age Computation of the Antarctic Ice Sheet for Dating Deep Ice Cores Bernd M¨ugge1 for the computation of the age of ice is discussed within the frame of numerical ice sheet modelling. The first method of a numerical diffusion term to stabilize the solution and therefore produces arbitrary results in a near
Calov, Reinhard
Adaptive higher order numerical simulation of heat and mass transfer in fluidized beds
Adaptive higher order numerical simulation of heat and mass transfer in fluidized beds Ch. Nagaiah1 adaptive numerical results of heat and mass transfer in fluidized beds using higher order time stepping injection. The numerical results are tested with different time stepping methods for different spatial grid
Magdeburg, Universität
A fast direct numerical simulation method for characterising hydraulic roughness
We describe a fast direct numerical simulation (DNS) method that promises to directly characterise the hydraulic roughness of any given rough surface, from the hydraulically smooth to the fully rough regime. The method circumvents the unfavourable computational cost associated with simulating high-Reynolds-number flows by employing minimal-span channels (Jimenez & Moin 1991). Proof-of-concept simulations demonstrate that flows in minimal-span channels are sufficient for capturing the downward velocity shift, that is, the Hama roughness function, predicted by flows in full-span channels. We consider two sets of simulations, first with modelled roughness imposed by body forces, and second with explicit roughness described by roughness-conforming grids. Owing to the minimal cost, we are able to conduct DNSs with increasing roughness Reynolds numbers while maintaining a fixed blockage ratio, as is typical in full-scale applications. The present method promises a practical, fast and accurate tool for character...
Chung, Daniel; MacDonald, Michael; Hutchins, Nicholas; Ooi, Andrew
Numerically Tracking Contact Discontinuities with an Introduction for GPU Programming
We review some of the classic numerical techniques used to analyze contact discontinuities and compare their effectiveness. Several finite difference methods (the Lax-Wendroff method, a Multidimensional Positive Definite Advection Transport Algorithm (MPDATA) method and a Monotone Upstream Scheme for Conservation Laws (MUSCL) scheme with an Artificial Compression Method (ACM)) as well as the finite element Streamlined Upwind Petrov-Galerkin (SUPG) method were considered. These methods were applied to solve the 2D advection equation. Based on our results we concluded that the MUSCL scheme produces the sharpest interfaces but can inappropriately steepen the solution. The SUPG method seems to represent a good balance between stability and interface sharpness without any inappropriate steepening. However, for solutions with discontinuities, the MUSCL scheme is superior. In addition, a preliminary implementation in a GPU program is discussed.
Davis, Sean L [Los Alamos National Laboratory
Progress report on LBL's numerical modeling studies on Cerro Prieto
An exploitation model of the Cerro Prieto geothermal system is needed to assess the energy capacity of the field, estimate its productive lifetime and develop an optimal reservoir management plan. The model must consider the natural state (i.e., pre-exploitation) conditions of the system and be able to predict changes in the reservoir thermodynamic conditions (and fluid chemistry) in response to fluid production (and injection). This paper discusses the results of a three-dimensional numerical simulation of the natural state conditions of the Cerro Prieto field and compares computed and observed pressure and temperature/enthalpy changes for the 1973--1987 production period. 16 refs., 24 figs., 2 tabs.
Halfman-Dooley, S.E.; Lippman, M.J.; Bodvarsson, G.S.
Numerical simulation of the Tayler instability in liquid metals
The electrical current through an incompressible, viscous and resistive liquid conductor produces an azimuthal magnetic field that becomes unstable when the corresponding Hartmann number exceeds a critical value in the order of 20. This Tayler instability, which is not only discussed as a key ingredient of a non-linear stellar dynamo model (Tayler-Spruit dynamo), but also as a limiting factor for the maximum size of large liquid metal batteries, was recently observed experimentally in a column of a liquid metal (Seilmayer et al., Phys. Rev. Lett. 108, 244501, 2012}. On the basis of an integro-differential equation approach, we have developed a fully three-dimensional numerical code, and have utilized it for the simulation of the Tayler instability at typical viscosities and resistivities of liquid metals. The resulting growth rates are in good agreement with the experimental data. We illustrate the capabilities of the code for the detailed simulation of liquid metal battery problems in realistic geometries.
Weber, Norbert; Stefani, Frank; Weier, Tom; Wondrak, Thomas
Spherically symmetric cosmological spacetimes with dust and radiation � numerical implementation
We present new numerical cosmological solutions of the Einstein Field Equations. The spacetime is spherically symmetric with a source of dust and radiation approximated as a perfect fluid. The dust and radiation are necessarily non-comoving due to the inhomogeneity of the spacetime. Such a model can be used to investigate non-linear general relativistic effects present during decoupling or big-bang nucleosynthesis, as well as for investigating void models of dark energy with isocurvature degrees of freedom. We describe the full evolution of the spacetime as well as the redshift and luminosity distance for a central observer. After demonstrating accuracy of the code, we consider a few example models, and demonstrate the sensitivity of the late time model to the degree of inhomogeneity of the initial radiation contrast.
Lim, Woei Chet [Department of Mathematics, University of Waikato, Private Bag 3105, Hamilton 3240 (New Zealand); Regis, Marco [Dipartimento di Fisica, Universit� di Torino and INFN, Torino (Italy); Clarkson, Chris, E-mail: firstname.lastname@example.com, E-mail: firstname.lastname@example.com, E-mail: email@example.com [Astrophysics, Cosmology and Gravity Centre, and Department of Mathematics and Applied Mathematics, University of Cape Town, Rondebosch 7701, Cape Town (South Africa)
Quantum Cosmology will need to become a Numerical Subject
The inhomogeneous fluctuations that underlie structure formation - galaxies and CMB hotspots - might have been seeded by quantum cosmological fluctuations, as magnified by some inflationary mechanism. The Halliwell-Hawking model for these, as a lower-energy semiclassical limit, is expected to be shared by many theories. E.g. an O((H/m_pl)^2) suppression of power at large scales results from this. This model contains/suppresses very many terms; we want a qualitative understanding of the meaning of these terms and of different regimes resulting from different combinations of them. I study this with toy models that have tractable mathematics: minisuperspace and, especially, relational particle mechanics. In the present Seminar, I consider in particular averaged terms with some lessons from Hartree-Fock approach to Atomic and Molecular Physics. One needs to anchor this on variational principles; treating the subsequent equations is a numerical venture.
Anderson, Edward
Edward Anderson
An acoustic numerical analysis of single rotation torpedo propeller
&vers&ty C'hatrman of Adv&sorv C'ummittee Dr Kenneth D Ivnrkan 5 numerical study of the acoust&c field asscsciated ivith a t?rped& pr &!&cll&-r hav&ng fnur blades i&as accomplished thr&?&gh the application nf Su& ci method The Succi acnustic analvsi ivas m... I-d&g&t se?es air!'n&1 data ba. e to calculate the lift and &&rag f ir each blade segment, These values vere then converted ti thrust and torque radial distr&butinns. Alsn a rav tracing techniq?e ivas inc&tip&i&at&. c! int& tlic Succi numencal...
Kim, Soo Yong
Numerical tools applied to power reactor noise analysis
In this paper, the development of numerical tools allowing the determination of the neutron noise in power reactors is reported. These tools give the space-dependence of the fluctuations of the neutron flux induced by fluctuating properties of the medium in the 2-group diffusion approximation and in a 2-dimensional representation of heterogeneous systems. Some applications of these tools to power reactor noise analysis are then described. These applications include the unfolding of the noise source from the resulting neutron noise measured at a few discrete locations throughout the core, the study of the space-dependence of the Decay Ratio in Boiling Water Reactors, the noise-based estimation of the Moderator Temperature Coefficient of reactivity in Pressurized Water Reactors, the modeling of shell-mode core barrel vibrations in Pressurized Water Reactors, and the investigation of the validity of the point-kinetic approximation in subcritical systems. (authors)
Demaziere, C.; Pazsit, I. [Chalmers Univ. of Technology, Dept. of Nuclear Engineering, SE-412 96 Goeteborg (Sweden)
Flow focusing in unsaturated fracture networks: A numerical investigation
A numerical modeling study is presented to investigate flow-focusing phenomena in a large-scale fracture network, constructed using field data collected from the unsaturated zone of Yucca Mountain, Nevada, the proposed repository site for high-level nuclear waste. The two-dimensional fracture network for an area of 100 m x 150 m contains more than 20,000 fractures. Steady-state unsaturated flow in the fracture network is investigated for different boundary conditions and rock properties. Simulation results indicate that flow paths are generally vertical, and that horizontal fractures mainly provide pathways between neighboring vertical paths. In addition to fracture properties, flow-focusing phenomena are also affected by rock-matrix permeability, with lower matrix permeability leading to a high degree of flow focusing. The simulation results further indicate that the average spacing between flow paths in a layered system tends to increase and flow tends to becomes more focused, with depth.
Zhang, Keni; Wu, Yu-Shu; Bodvarsson, G.S.; Liu, Hui-Hai
Basalt, Colorado: Energy Resources | Open Energy Information
AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home5b9fcbce19 NoPublic Utilities Address: 160 EastMaine:Barbers Point Housing,Illinois: Energy Resources Jump8150781°
Topological insulators and C*-algebras: Theory and numerical practice
Research Highlights: > We classify topological insulators using C* algebras. > We present new K-theory invariants. > We develop efficient numerical algorithms based on this technique. > We observe unexpected quantum phase transitions using our algorithm. - Abstract: We apply ideas from C*-algebra to the study of disordered topological insulators. We extract certain almost commuting matrices from the free Fermi Hamiltonian, describing band projected coordinate matrices. By considering topological obstructions to approximating these matrices by exactly commuting matrices, we are able to compute invariants quantifying different topological phases. We generalize previous two dimensional results to higher dimensions; we give a general expression for the topological invariants for arbitrary dimension and several symmetry classes, including chiral symmetry classes, and we present a detailed K-theory treatment of this expression for time reversal invariant three dimensional systems. We can use these results to show non-existence of localized Wannier functions for these systems. We use this approach to calculate the index for time-reversal invariant systems with spin-orbit scattering in three dimensions, on sizes up to 12{sup 3}, averaging over a large number of samples. The results show an interesting separation between the localization transition and the point at which the average index (which can be viewed as an 'order parameter' for the topological insulator) begins to fluctuate from sample to sample, implying the existence of an unsuspected quantum phase transition separating two different delocalized phases in this system. One of the particular advantages of the C*-algebraic technique that we present is that it is significantly faster in practice than other methods of computing the index, allowing the study of larger systems. In this paper, we present a detailed discussion of numerical implementation of our method.
Hastings, Matthew B., E-mail: firstname.lastname@example.com [Department of Physics, Duke University, Durham, NC 27708 (United States); Microsoft Research, Station Q, Elings Hall, University of California, Santa Barbara, CA 93106 (United States); Department of Mathematics and Statistics, University of New Mexico, Albuquerque, NM 87131 (United States); Loring, Terry A. [Department of Physics, Duke University, Durham, NC 27708 (United States); Microsoft Research, Station Q, Elings Hall, University of California, Santa Barbara, CA 93106 (United States); Department of Mathematics and Statistics, University of New Mexico, Albuquerque, NM 87131 (United States)
Geology of the Waste Treatment Plant Seismic Boreholes
In 2006, DOE-ORP initiated the Seismic Boreholes Project (SBP) to emplace boreholes at the Waste Treatment Plant (WTP) site in order to obtain direct Vs measurements and other physical property measurements in Columbia River basalt and interbedded sediments of the Ellensburg Formation. The goal was to reduce the uncertainty in the response spectra and seismic design basis, and potentially recover design margin for the WTP. The characterization effort within the deep boreholes included 1) downhole measurements of the velocity properties of the suprabasalt, basalt, and sedimentary interbed sequences, 2) downhole measurements of the density of the subsurface basalt and sediments, and 3) confirmation of the geometry of the contact between the various basalt and interbedded sediments through examination of retrieved core from the corehole and data collected through geophysical logging of each borehole. This report describes the results of the geologic studies from three mud-rotary boreholes and one cored borehole at the WTP. All four boreholes penetrated the entire Saddle Mountains Basalt and the upper part of the Wanapum Basalt where thick sedimentary interbeds occur between the lava flows. The basalt flows penetrated in Saddle Mountains Basalt included the Umatilla Member, Esquatzel Member, Pomona Member and the Elephant Mountain Member. The underlying Priest Rapids Member of the Wanapum Basalt was also penetrated. The Ellensburg Formation sediments consist of the Mabton Interbed, the Cold Creek Interbed, the Selah Interbed and the Rattlesnake Ridge Interbed; the Byron Interbed occurs between two flows of the Priest Rapids Member. The Mabton Interbed marks the contact between the Wanapum and Saddle Mountains Basalts. The thicknesses of the basalts and interbedded sediments were within expected limits. However, a small reverse fault was found in the Pomona Member flow top. This fault has three periods of movement and less than 15 feet of repeated section. Most of the movement on the fault appears to have occurred before the youngest lava flow, the 10.5 million year old Elephant Mountain Member was emplaced above the Pomona Member.
Barnett, D. BRENT; Bjornstad, Bruce N.; Fecht, Karl R.; Lanigan, David C.; Reidel, Steve; Rust, Colleen F.
The Vishniac instability is thought to explain the complex structure of radiative supernova remnants in their Pressure-Driven Thin Shell (PDTS) phase after a blast wave (BW) has propagated from a central explosion. In this paper, the propagation of the BW and the evolution of the PDTS stage are studied numerically with the two-dimensional (2D) code HYDRO-MUSCL for a finite-thickness shell expanding in the interstellar medium (ISM). Special attention is paid to the adiabatic index, {gamma}, and three distinct values are taken for the cavity ({gamma}{sub 1}), the shell ({gamma}{sub 2}), and the ISM ({gamma}{sub 3}) with the condition {gamma}{sub 2} < {gamma}{sub 1}, {gamma}{sub 3}. This low value of {gamma}{sub 2} accounts for the high density in the shell achieved by a strong radiative cooling. Once the spherical background flow is obtained, the evolution of a 2D-axisymmetric perturbation is computed from the linear to the nonlinear regime. The overstable mechanism, previously demonstrated theoretically by E. T. Vishniac in 1983, is recovered numerically in the linear stage and is expected to produce and enhance anisotropies and clumps on the shock front, leading to the disruption of the shell in the nonlinear phase. The period of the increasing oscillations and the growth rate of the instability are derived from several points of view (the position of the perturbed shock front, mass fluxes along the shell, and density maps), and the most unstable mode differing from the value given by Vishniac is computed. In addition, the influence of several parameters (the Mach number, amplitude and wavelength of the perturbation, and adiabatic index) is examined and for wavelengths that are large enough compared to the shell thickness, the same conclusion arises: in the late stage of the evolution of the radiative supernova remnant, the instability is dampened and the angular initial deformation of the shock front is smoothed while the mass density becomes uniform with the angle. As a result, our model shows that the supernova remnant returns to a stable evolution and the Vishniac instability does not lead to the fragmentation of the shock as predicted by the theory.
Michaut, C.; Cavet, C.; Bouquet, S. E.; Roy, F.; Nguyen, H. C., E-mail: email@example.com [LUTH, Observatoire de Paris, CNRS, Universite Paris-Diderot, F-92190 Meudon (France)
Numerical simulations of super-luminous supernovae of type IIn
We present numerical simulations that include 1-D Eulerian multi-group radiation-hydrodynamics, 1-D non-LTE radiative transfer, and 2-D polarised radiative transfer for super-luminous interacting supernovae (SNe). Our reference model is a ~10Msun inner shell with 10^51erg ramming into a ~3Msun cold outer shell (the circumstellar-medium, or CSM) that extends from 10^15cm to 2x10^16cm and moves at 100km/s. We discuss the light curve evolution, which cannot be captured adequately with a grey approach. In these interactions, the shock-crossing time through the optically-thick CSM is much longer than the photon diffusion time. Radiation is thus continuously leaking from the shock through the CSM, in disagreement with the shell-shocked model that is often invoked. Our spectra redden with time, with a peak distribution in the near-UV during the first month gradually shifting to the optical range over the following year. Initially Balmer lines exhibit a narrow line core and the broad line wings that are characteristi...
Dessart, Luc; Hillier, D John
Topological Insulators and C^*-Algebras: Theory and Numerical Practice
We apply ideas from $C^*$-algebra to the study of disordered topological insulators. We extract certain almost commuting matrices from the free Fermi Hamiltonian, describing band projected coordinate matrices. By considering topological obstructions to approximating these matrices by exactly commuting matrices, we are able to compute invariants quantifying different topological phases. We generalize previous two dimensional results to higher dimensions; we give a general expression for the topological invariants for arbitrary dimension and several symmetry classes, including chiral symmetry classes, and we present a detailed $K$-theory treatment of this expression for time reversal invariant three dimensional systems. We can use these results to show non-existence of localized Wannier functions for these systems. We use this approach to calculate the index for time-reversal invariant systems with spin-orbit scattering in three dimensions, on sizes up to $12^3$, averaging over a large number of samples. The results show an interesting separation between the localization transition and the point at which the average index (which can be viewed as an "order parameter" for the topological insulator) begins to fluctuate from sample too sample, implying the existence of an unsuspected quantum phase transition separating two different delocalized phases in this system. One of the particular advantages of the $C^*$-algebraic technique that we present is that it is significantly faster in practice than other methods of computing the index, allowing the study of larger systems. In this paper, we present a detailed discussion of numerical implementation of our method.
M. B. Hastings; T. A. Loring
Numerical simulation of laminar flow in a curved duct
This paper describes numerical simulations that were performed to study laminar flow through a square duct with a 900 bend. The purpose of this work was two fold. First, an improved understanding was desired of the flow physics involved in the generation of secondary vortical flows in three-dimensions. Second, adaptive gridding techniques for structured grids in three- dimensions were investigated for the purpose of determining their utility in low Reynolds number, incompressible flows. It was also of interest to validate the commercial computer code CFD-ACE. Velocity predictions for both non-adaptive and adaptive grids are compared with experimental data. Flow visualization was used to examine the characteristics of the flow though the curved duct in order to better understand the viscous flow physics of this problem. Generally, moderate agreement with the experimental data was found but shortcomings in the experiment were demonstrated. The adaptive grids did not produce the same level of accuracy as the non-adaptive grid with a factor of four more grid points.
Lopez, A.R.; Oberkampf, W.L.
Determination of hydrogen cluster velocities and comparison with numerical calculations
The use of powerful hydrogen cluster jet targets in storage ring experiments led to the need of precise data on the mean cluster velocity as function of the stagnation temperature and pressure for the determination of the volume density of the target beams. For this purpose a large data set of hydrogen cluster velocity distributions and mean velocities was measured at a high density hydrogen cluster jet target using a trumpet shaped nozzle. The measurements have been performed at pressures above and below the critical pressure and for a broad range of temperatures relevant for target operation, e.g., at storage ring experiments. The used experimental method is described which allows for the velocity measurement of single clusters using a time-of-flight technique. Since this method is rather time-consuming and these measurements are typically interfering negatively with storage ring experiments, a method for a precise calculation of these mean velocities was needed. For this, the determined mean cluster velocities are compared with model calculations based on an isentropic one-dimensional van der Waals gas. Based on the obtained data and the presented numerical calculations, a new method has been developed which allows to predict the mean cluster velocities with an accuracy of about 5%. For this two cut-off parameters defining positions inside the nozzle are introduced, which can be determined for a given nozzle by only two velocity measurements.
T�schner, A.; K�hler, E.; Ortjohann, H.-W.; Khoukaz, A. [Institut f�r Kernphysik, Westf�lische Wilhelms-Universit�t M�nster, D-48149 M�nster (Germany)] [Institut f�r Kernphysik, Westf�lische Wilhelms-Universit�t M�nster, D-48149 M�nster (Germany)
Numerical simulation of linear fiction welding (LFW) processes
Solid state welding processes are becoming increasingly important due to a large number of advantages related to joining ''unweldable'' materials and in particular light weight alloys. Linear friction welding (LFW) has been used successfully to bond non-axisymmetric components of a range of materials including titanium alloys, steels, aluminum alloys, nickel, copper, and also dissimilar material combinations. The technique is useful in the research of quality of the joints and in reducing costs of components and parts of the aeronautic and automotive industries.LFW involves parts to be welded through the relative reciprocating motion of two components under an axial force. In such process the heat source is given by the frictional forces work decaying into heat determining a local softening of the material and proper bonding conditions due to both the temperature increase and the local pressure of the two edges to be welded. This paper is a comparative test between the numerical model in two dimensions, i.e. in plane strain conditions, and in three dimensions of a LFW process of AISI1045 steel specimens. It must be observed that the 3D model assures a faithful simulation of the actual threedimensional material flow, even if the two-dimensional simulation computational times are very short, a few hours instead of several ones as the 3D model. The obtained results were compared with experimental values found out in the scientific literature.
Fratini, L.; La Spisa, D. [University of Palermo-Dept. of Industrial engineering (Italy)
Evolutions in 3D numerical relativity using fixed mesh refinement
We present results of 3D numerical simulations using a finite difference code featuring fixed mesh refinement (FMR), in which a subset of the computational domain is refined in space and time. We apply this code to a series of test cases including a robust stability test, a nonlinear gauge wave and an excised Schwarzschild black hole in an evolving gauge. We find that the mesh refinement results are comparable in accuracy, stability and convergence to unigrid simulations with the same effective resolution. At the same time, the use of FMR reduces the computational resources needed to obtain a given accuracy. Particular care must be taken at the interfaces between coarse and fine grids to avoid a loss of convergence at higher resolutions, and we introduce the use of "buffer zones" as one resolution of this issue. We also introduce a new method for initial data generation, which enables higher-order interpolation in time even from the initial time slice. This FMR system, "Carpet", is a driver module in the freely available Cactus computational infrastructure, and is able to endow generic existing Cactus simulation modules ("thorns") with FMR with little or no extra effort.
Erik Schnetter; Scott H. Hawley; Ian Hawke
Measurement and numerical calculation of Rubidium Rydberg Stark spectra
We report on the measurement of Stark shifted energy levels of $^{87}$Rb Rydberg atoms in static electric fields by means of electromagnetically induced transparency (EIT). Electric field strengths of up to 500V/cm, ranging beyond the classical ionisation threshold, were applied using electrodes inside a glass cell with rubidium vapour. Stark maps for principal quantum numbers $n=35$ and $n=70$ have been obtained with high signal-to-noise ratio for comparison with results from ab initio calculations following the method described in [M. L. Zimmerman et al., Phys. Rev. A 20, 2251 (1979)], which was originally only verified for states around $n=15$. We also calculate the dipole matrix elements between low-lying states and Stark shifted Rydberg states to give a theoretical estimate of the relative strength of the EIT signal. The present work significantly extends the experimental verification of this numerical method in the range of both high principal quantum numbers and high electric fields with an accuracy of...
Grimmel, Jens; Karlewski, Florian; Jessen, Florian; Reinschmidt, Malte; S�ndor, N�ra; Fort�gh, J�zsef
The Dynamics of Deterministic Chaos in Numerical Weather Prediction Models
Atmospheric weather systems are coherent structures consisting of discrete cloud cells forming patterns of rows/streets, mesoscale clusters and spiral bands which maintain their identity for the duration of their appreciable life times in the turbulent shear flow of the planetary Atmospheric Boundary Layer. The existence of coherent structures (seemingly systematic motion) in turbulent flows has been well established during the last 20 years of research in turbulence. Numerical weather prediction models based on the inherently non-linear Navier-Stokes equations do not give realistic forecasts because of the following inherent limitations: (1) the non-linear governing equations for atmospheric flows do not have exact analytic solutions and being sensitive to initial conditions give chaotic solutions characteristic of deterministic chaos (2) the governing equations do not incorporate the dynamical interactions and co-existence of the complete spectrum of turbulent fluctuations which form an integral part of the large coherent weather systems (3) limitations of available computer capacity necessitates severe truncation of the governing equations, thereby generating errors of approximations (4) the computer precision related roundoff errors magnify the earlier mentioned uncertainties exponentially with time and the model predictions become unrealistic. The accurate modelling of weather phenomena therefore requires alternative concepts and computational techniques. In this paper a universal theory of deterministic chaos applicable to the formation of coherent weather structures in the ABL is presented.
A. Mary Selvam
Exploring inequality violations by classical hidden variables numerically
There are increasingly suggestions for computer simulations of quantum statistics which try to violate Bell type inequalities via classical, common cause correlations. The Clauser�Horne�Shimony�Holt (CHSH) inequality is very robust. However, we argue that with the Einstein�Podolsky�Rosen setup, the CHSH is inferior to the Bell inequality, although and because the latter must assume anti-correlation of entangled photon singlet states. We simulate how often quantum behavior violates both inequalities, depending on the number of photons. Violating Bell 99% of the time is argued to be an ideal benchmark. We present hidden variables that violate the Bell and CHSH inequalities with 50% probability, and ones which violate Bell 85% of the time when missing 13% anti-correlation. We discuss how to present the quantum correlations to a wide audience and conclude that, when defending against claims of hidden classicality, one should demand numerical simulations and insist on anti-correlation and the full amount of Bell violation. -- Highlights: �The widely assumed superiority of the CHSH fails in the EPR problem. �We simulate Bell type inequalities behavior depending on the number of photons. �The core of Bell�s theorem in the EPR setup is introduced in a simple way understandable to a wide audience. �We present hidden variables that violate both inequalities with 50% probability. �Algorithms have been supplied in form of Mathematica programs.
Vongehr, Sascha, E-mail: email@example.com
Numerical Relativity in Spherical Polar Coordinates: Off-center Simulations
We have recently presented a new approach for numerical relativity simulations in spherical polar coordinates, both for vacuum and for relativistic hydrodynamics. Our approach is based on a reference-metric formulation of the BSSN equations, a factoring of all tensor components, as well as a partially implicit Runge-Kutta method, and does not rely on a regularization of the equations, nor does it make any assumptions about the symmetry across the origin. In order to demonstrate this feature we present here several off-centered simulations, including simulations of single black holes and neutron stars whose center is placed away from the origin of the coordinate system, as well as the asymmetric head-on collision of two black holes. We also revisit our implementation of relativistic hydrodynamics and demonstrate that a reference-metric formulation of hydrodynamics together with a factoring of all tensor components avoids problems related to the coordinate singularities at the origin and on the axes. As a particularly demanding test we present results for a shock wave propagating through the origin of the spherical polar coordinate system.
Thomas W. Baumgarte; Pedro J. Montero; Ewald M�ller
We have recently presented a new approach for numerical relativity simulations in spherical polar coordinates, both for vacuum and for relativistic hydrodynamics. Our approach is based on a reference-metric formulation of the BSSN equations, a factoring of all tensor components, as well as a partially implicit Runge-Kutta method, and does not rely on a regularization of the equations, nor does it make any assumptions about the symmetry across the origin. In order to demonstrate this feature we present here several off-centered simulations, including simulations of single black holes and neutron stars whose center is placed away from the origin of the coordinate system, as well as the asymmetric head-on collision of two black holes. We also revisit our implementation of relativistic hydrodynamics and demonstrate that a reference-metric formulation of hydrodynamics together with a factoring of all tensor components avoids problems related to the coordinate singularities at the origin and on the axes. As a parti...
Baumgarte, Thomas W; M�ller, Ewald
Numerical evaluation of multi-gluon amplitudes for High Energy Factorization
We present a program to evaluate tree-level multi-gluon amplitudes with up to two of them off-shell. Furthermore, it evaluates squared amplitudes summed over colors and helicities for up to six external gluons. It employs both analytic expressions, obtained via BCFW recursion, and numerical BCFW recursion. It has been validated numerically with the help of an independent program employing numerical Dyson-Schwinger recursion.
Bury, M
Numerical Methods for Forward and Inverse Problems in Discontinuous Media
The research emphasis under this grant�??s funding is in the area of algebraic multigrid methods. The research has two main branches: 1) exploring interdisciplinary applications in which algebraic multigrid can make an impact and 2) extending the scope of algebraic multigrid methods with algorithmic improvements that are based in strong analysis.The work in interdisciplinary applications falls primarily in the field of biomedical imaging. Work under this grant demonstrated the effectiveness and robustness of multigrid for solving linear systems that result from highly heterogeneous finite element method models of the human head. The results in this work also give promise to medical advances possible with software that may be developed. Research to extend the scope of algebraic multigrid has been focused in several areas. In collaboration with researchers at the University of Colorado, Lawrence Livermore National Laboratory, and Los Alamos National Laboratory, the PI developed an adaptive multigrid with subcycling via complementary grids. This method has very cheap computing costs per iterate and is showing promise as a preconditioner for conjugate gradient. Recent work with Los Alamos National Laboratory concentrates on developing algorithms that take advantage of the recent advances in adaptive multigrid research. The results of the various efforts in this research could ultimately have direct use and impact to researchers for a wide variety of applications, including, astrophysics, neuroscience, contaminant transport in porous media, bi-domain heart modeling, modeling of tumor growth, and flow in heterogeneous porous media. This work has already led to basic advances in computational mathematics and numerical linear algebra and will continue to do so into the future.
Timothy P. Chartier
Numerical controlled polishing, continued force wear and part correction experiments
This abstract reports the near completion of the first phase of this program. It is the aim of this program to provide the operator of a N/C diamond turning machine or N/C grinding machine (jig grinder) with the wear characteristics necessary to achieve uniform material removal. The second phase of this program addresses a different problem, although solving this problem is highly dependent on the results of the first phase. Diamond turned, or any lathe turned surface, exhibits regular tool marks due to the tool passing over the surface being cut. Changes in depth of cut, feed rate and work rpm will change the character of these groves, but will not eliminate them. Optical surfaces produced by this process exhibit increased scattering as the light wavelength decreases limiting their use; at least for optical purposes, to IR and some visible applications. Utilizing wear information gathered in the first part of this program we will attempt to reduce these residual tool marks by polishing. The polishing of diamond turned surfaces is not new. Diamond turned metal surfaces, especially in electroless nickel and high phosphorus nickel electroplate have been polished to improve their scatter characteristics. What we believe is unique is the use of a spherical wheel, rotating on axis and being moved over the part in a prescribed manner by numerical control. Over the past year we have made some major changes in our polishing methods and procedures. We have listed below these changes, as a refresher for the reader as to our previous procedures. These changes will be addressed in the body of the text.
Hannah, P.R.; Day, R.D.; Hatch, D.J.; McClure, E.R.
Representing Cloud Processing of Aerosol in Numerical Models
The satellite imagery in Figure 1 provides dramatic examples of how aerosol influences the cloud field. Aerosol from ship exhaust can serve as nucleation centers in otherwise cloud-free regions, forming ship tracks (top image), or can enhance the reflectance/albedo in already cloudy regions. This image is a demonstration of the first indirect effect, in which changes in aerosol modulate cloud droplet radius and concentration, which influences albedo. It is thought that, through the effects it has on precipitation (drizzle), aerosol can also affect the structure and persistence of planetary boundary layer (PBL) clouds. Regions of cellular convection, or open pockets of cloudiness (bottom image) are thought to be remnants of strongly drizzling PBL clouds. Pockets of Open Cloudiness (POCs) (Stevens et al. 2005) or Albrecht's ''rifts'' are low cloud fraction regions characterized by anomalously low aerosol concentrations, implying they result from precipitation. These features may in fact be a demonstration of the second indirect effect. To accurately represent these clouds in numerical models, we have to treat the coupled cloud-aerosol system. We present the following series of mesoscale and large eddy simulation (LES) experiments to evaluate the important aspects of treating the coupled cloud-aerosol problem. 1. Drizzling and nondrizzling simulations demonstrate the effect of drizzle on a mesoscale forecast off the California coast. 2. LES experiments with explicit (bin) microphysics gauge the relative importance of the shape of the aerosol spectrum on the 3D dynamics and cloud structure. 3. Idealized mesoscale model simulations evaluate the relative roles of various processes, sources, and sinks.
Mechem, D.B.; Kogan, Y.L.
MATLAB muliplatform research license: Numeric computation, advanced graphics and visualization, and a high-level programming
MATLAB muliplatform research license: Numeric computation, advanced graphics and visualization license allows research and includes the following: Feature # of users MATLAB 7 Simulink 6 Control System
Dawson, Jeff W.
Numerical Study of Freestream Waves Receptivity and Nonlinear Breakdown in Hypersonic Boundary Layer
Numerical Methods for Hypersonic Boundary Layer Stability.of Instability in a Hypersonic Boundary Layer. TheoreticalA. P. , Receptivity of Hypersonic Boundary Layer to Wall
Lei, Jia
High Resolution Backside Imaging and Thermography using a Numerical Aperture Increasing Lens
High Resolution Backside Imaging and Thermography using a Numerical Aperture Increasing Lens M. S, representing the highest resolution subsurface thermography to date. Keywords: thermal imaging, high
Gravity effects on partially premixed flames: an experimental-numerical investigation
Gravity effects on partially premixed flames: an experimental-numerical investigation Andrew J and interactions between the various reaction zones are strongly influenced by gravity. The flames widen
Aggarwal, Suresh K.
Numerical tools for musical instruments acoustics: analysing nonlinear physical models using continuation
self-sustained oscillations. Several physical models (clarinet and saxophone) are formulated the dynamical properties of self-sustained musical instruments using tra- ditional numerical techniques
Numerical Modeling of Failure in Magnesium Alloys under Axial Compression and Bending for Crashworthiness Applications.
??Numerical modeling of failure was performed for magnesium alloys with circular and square cross-sections under axial compression. The failure criterion was employed using material model� (more)
Ali, Usman
13.002J Introduction to Numeric Analysis for Engineering, Spring 2003
An introduction to the formulation, methodology, and techniques for numerical solution of engineering problems. Fundamental principles of digital computing and the implications for algorithm accuracy and stability. Error ...
Leonard, John J.
E-Print Network 3.0 - acoustic sources numerical Sample Search...
analytical or numerical so- lutions. The natural physical and acoustical data are limited in accuracy... 400 150 250 350 SOURCE RECEIVER 1 RECEIVER 2 Simulation Experiment ......
E-Print Network 3.0 - adaptive numerical dissipative Sample Search...
rate. Through numerical simulations, we find... description (i) Dissipative, vortex-vortex forces ... Source: Grzybowski, Bartosz A. - Departments of Chemistry & Chemical...
Monitoring Erebus volcano's active lava lake; tools, techniques and observations
serves the crater rim was replaced. The old system was approaching ten years of service, and had become overly com- plex, and unreliable. A combination of extreme winds and corrosive gases makes the crater rim of Erebus an unsuitable site for solar panels... of the lake. I classify such events as being distinct from the far more frequently occurring �small�, metre and sub-metre scale bubbles (e.g. panel (d) in Fig. 4.9) which arrive at the surface of the lake, but do not result in explosive activity. A study...
Peters, Nial John
Lava Hot Springs Space Heating Low Temperature Geothermal Facility | Open
AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data CenterFranconia, Virginia: Energy Resources Jump to:46 - 429Lacey, Washington:Lakeville, MN)Lauderhill, Florida: EnergyLaurelDome
Direct Numerical Simulation of a Hypersonic Turbulent Boundary Layer on a Large Domain
Direct Numerical Simulation of a Hypersonic Turbulent Boundary Layer on a Large Domain Stephan Priebe , M. Pino Mart´in The direct numerical simulation (DNS) of a spatially-developing hypersonic There are few studies of hypersonic flows at Mach number greater than 5 and few involve the measurement of mean
Martín, Pino
Numerical evaluation of one-loop diagrams near exceptional momentum configurations
One problem which plagues the numerical evaluation of one-loop Feynman diagrams using recursive integration by part relations is a numerical instability near exceptional momentum configurations. In this contribution we will discuss a generic solution to this problem. As an example we consider the case of forward light-by-light scattering.
Walter T Giele; Giulia Zanderighi; E.W.N. Glover
Integration of Mathematics/Numeric Analysis with Chemistry/Chemical Engineering
1 Integration of Mathematics/Numeric Analysis with Chemistry/Chemical Engineering Claes Niklasson chemical engineering programs, such as fundamental chemistry, chemical reaction engineering and bioprocess/numerical analysis and chemistry/chemical engineering will be described and the implications of this change
Numerical Analysis of Non-constant Discounting with an Application to Renewable Resource Management
Numerical Analysis of Non-constant Discounting with an Application to Renewable Resource Management illustrate the approach by studying welfare and observational equivalence for a particular renewable resource man- agement problem. Keywords: Non-constant discounting, numerical methods, non-renewable resources
Karp, Larry S.
Numerical tools applied to power reactor noise analysis Christophe Demazie`re*, Imre Pazsit
Review Numerical tools applied to power reactor noise analysis Christophe Demazie`re*, Imre Pa and application of the numerical tools employed. The code that was developed yields the space and non-critical systems with an external source. Some appli- cations of these tools to power reactor
Pázsit, Imre
Numerical simulation of the heat transfer in amorphous silicon nitride membrane-based microcalorimeters
Numerical simulation of the heat transfer in amorphous silicon nitride membrane July 2003 Numerical simulations of the two-dimensional 2D heat flow in a membrane-based microcalorimeter have been performed. The steady-state isotherms and time-dependent heat flow have been calculated
Hellman, Frances
1 NUMERICAL SIMULATION OF POOL BOILING FOR STEADY STATE AND TRANSIENT HEATING Ying He, Masahiro role in nucleate and transition boiling heat transfer at high heat flux. Many experiments have been in the numerical simulation of boiling heat transfer. In this study, based on the macrolayer evaporation model
A Numerical quadrature for the Schwarz-Chimera J. -B. Apoung Kamga1
A Numerical quadrature for the Schwarz-Chimera Method J. -B. Apoung Kamga1 and Olivier Pironneau2 1 Curie email@example.com Abstract Chimera [10] happens to be a version of Schwarz' method fvih, (8) #12;A Numerical quadrature for the Schwarz-Chimera Method 3 Theorem 2. (see Hecht et al. [4
Pironneau, Olivier
NUMERICAL STUDIES OF COSMIC-RAY INJECTION AND ACCELERATION Hyesung Kang
; accepted 2002 July 9 ABSTRACT A numerical scheme that incorporates a thermal leakage injection model in the precursor region, especially in front of strong, highly modified shocks. The `` thermal leakage '' injectionNUMERICAL STUDIES OF COSMIC-RAY INJECTION AND ACCELERATION Hyesung Kang Department of Earth
Gieseler, Udo D. J.
NUMERICAL INVESTIGATION OF TEMPERATURE EFFECTS DURING THE INJECTION OF CARBON DIOXIDE INTO BRINE for the simulation of carbon dioxide injection into geological formations is currently an intensive field of research for the balance of thermal energy, we can investigate numerically the effects of temperature variations during
423A HIGH-PERFORMANCE COMPUTING/NUMERICAL The International Journal of High Performance Computing Applications,
423A HIGH-PERFORMANCE COMPUTING/NUMERICAL The International Journal of High Performance Computing and barriers in the development of high-performance computing (HPC) algorithms and software. The activity has computing, numerical analy- sis, roadmap, applications and algorithms, software 1 The High-performance
EVALUATION OF NUMERICAL WEATHER PREDICTION IN MODELING CLOUD- RADIATION INTERACTIONS OVER.bnl.gov ABSTRACT Numerical weather prediction (NWP) is the basis for present-day weather forecasts, and NWP for Medium-Range Weather Forecasts, the US North American Model, and the US Global Forecast System. Attempts
Johnson, Peter D.
Time Step Sensitivity of Nonlinear Atmospheric Models: Numerical Convergence, Truncation Error Growth, and Ensemble Design
1973; Oran and Boris 1987; Murray 1989; Gershenfeld 1999). Weather and climate prediction models, which to the initial conditions, which is a major source of uncertainty in Numerical Weather Prediction (NWP; eTime Step Sensitivity of Nonlinear Atmospheric Models: Numerical Convergence, Truncation Error
Judd, Kevin
Numerical Weather Prediction (NWP) and hybrid ARMA/ANN model to predict global radiation
Numerical Weather Prediction (NWP) and hybrid ARMA/ANN model to predict global radiation Cyril a hybrid ARMA/ANN model and data issued from a numerical weather prediction model (ALADIN). We email@example.com #12;Abstract. We propose in this paper an original technique to predict global radiation using
Hamiltonian-based numerical methods for forced-dissipative climate prediction
Hamiltonian-based numerical methods for forced-dissipative climate prediction Bob Peeters1 , Onno long-term weather forecast models fail at this point. But the question remains, however: Question: Is it advantageous to use numerical schemes with a Hamil- tonian core for realistic climate modeling? The primitive
Al Hanbali, Ahmad
Numerical study on optimal Stirling engine regenerator matrix designs taking into account the effects of matrix
Numerical study on optimal Stirling engine regenerator matrix designs taking into account matrix design that improves the efficiency of a Stirling engine has been developed in a numerical study of the existing SM5 Stirling engine. A new, detailed, one-dimensional Stirling engine model that delivers results
Numerical analysis of a model for Nickel-Iron alloy electrodeposition on rotating disk
Numerical analysis of a model for Nickel-Iron alloy electrodeposition on rotating disk electrode N the nickel-iron electrodeposition process, we have developed one-dimensional numerical model. This model ad can predict characteristic features of the nickel-iron sys- tem. this work was supported
Paris-Sud XI, Universit� de
Separation phenomena in the tritium source and numerical simulations of turbo-molecular pumps
Separation phenomena in the tritium source and numerical simulations of turbo-molecular pumps Felix In the previous works [1, 2], the results of numerical calculations of tritium flow from the buffer vessel up to the first vacuum system were reported. Two values of the tritium source temperature were considered, i.e. 27
Sharipov, Felix
Effects of Shear Rate on Propagation of Blood Clotting Determined Using Microfluidics and Numerical Simulations
Effects of Shear Rate on Propagation of Blood Clotting Determined Using Microfluidics and email@example.com Abstract: This paper describes microfluidic experiments with human blood plasma and numerical simulations removed. In addition, these results demonstrate the utility of simplified mechanisms and microfluidics
Ismagilov, Rustem F.
Comparison of Numerical Techniques Used for Simulating Variable-Density Flow and Transport Experiments
Comparison of Numerical Techniques Used for Simulating Variable-Density Flow and Transport code by simulating two new variable-density-flow and transport experimental data sets. The experiments transport of a sinking groundwater plume and a rising groundwater plume. The numerical techniques used
Clement, Prabhakar
Maps preserving the joint numerical radius distance of operators Chi-Kwong Li
Maps preserving the joint numerical radius distance of operators Chi-Kwong Li and Edward Poon Dedicated to Professor David Lutzer on the occasion of his retirement. Abstract Denote the joint numerical radius of an m-tuple of bounded operators A = (A1, . . . , Am) by w(A). We give a complete description
1. Department, Course Number, Title ORE 766 Numerical Methods in Ocean Engineering
1. Department, Course Number, Title ORE 766 Numerical Methods in Ocean Engineering 2. Designation of numerical methods for simulating and solving ocean engineering problems. Topics include: Mathematical, & engineering Program Outcome 5: Use of latest tools in ocean engineering Program Outcome 6: Problem formulation
Frandsen, Jannette B.
Numerical assessment of stability criteria from disturbance energies in gaseous combustion
Numerical assessment of stability criteria from disturbance energies in gaseous combustion A, which corresponds to a ducted, laminar premixed propane-air flame, is used to assess the different terms a contribution from the unsteady heat flux-pressure correlation, allows a better agreement with the numerical
Nicoud, Franck
1 ON ACCURACY OF NUMERICAL EMC/EMI MODELING OVER A WIDE FREQUENCY RANGE Sergey Loyka EMC Lab: firstname.lastname@example.com Abstract - Numerical EMC/EMI modeling over a wide frequency range requires computational efficiency is proposed. I. INTRODUCTION Almost all the EMC problems are wide frequency range ones
Loyka, Sergey
OBJECT ORIENTED PROGRAMMING TECHNIQUES AND FAC METHOD IN NUMERICAL RESERVOIR SIMULATION \\Lambda
OBJECT ORIENTED PROGRAMMING TECHNIQUES AND FAC METHOD IN NUMERICAL RESERVOIR SIMULATION \\Lambda in numerical simulation of flow through hydrocarbon reservoirs within limitations in computing time and memory. These consist of solution of the conservation equations whichs govern the motion of fluid through the reservoir
A Numerical Study of Multiphase Phenomena in Rayleigh-Taylor Mixing D. Sendersky and D. Saltz
of the Euler equations in two space dimensions. The data set, generated by the front- tracking code FronTier problem in a finite domain [9]. The numerical solution is accomplished using the hydrodynamics code FronTier. The FronTier code utilizes the front tracking algorithm, a numerical method in which surfaces
New York at Stoney Brook, State University of
Intelligent Splitting for Disjunctive Numerical Thomas Douillard, Christophe Jermann, and Frederic Benhamou
- ever, they focus essentially on the pruning operation. In this paper, we present experimental evidences First-order numerical CSPs A numerical CSP (NCSP) is composed of variables whose domains are subsets by the compositional structure of the mixer: M = (A (E I S)) (V (C D True)). Formalisms derived from CSP have
Course MA59800: Numerical Simulation in Applied Geophysics. From the Mesoscale to the Macroscale.
Description Wave propagation is a common technique used in hydrocarbon exploration geophysics, mining and reservoir characterization and production, among other fields. Local variations in the fluid and solid gradients via a slow-wave diffusion process that can be analyzed using numerical experiments. Numerical rock
t -software package for numerical simulations of radioactive contaminant transport in groundwater
r3 t - software package for numerical simulations of radioactive contaminant transport equations that arise from the modelling of radioactive contaminant transport in porous media. It can solve, see [6]) can help to numerically simulate the spreading of radioactive contaminants in flowing ground
Frolkovic, Peter
Numerical modeling of wave propagation in random anisotropic heterogeneous elastic media
Numerical modeling of wave propagation in random anisotropic heterogeneous elastic media Q.-A. Ta numerical experiments that were performed on wave propagation in a randomly generated anisotropic used for the propagation of waves in geophysical media are not compatible with the surface recordings
Two-dimensional numerical methods in electromagnetic hypersonics including fully coupled Maxwell equations
the magnetic field source (solenoid or permanent magnet). To test the technique, we show the results obtainedTwo-dimensional numerical methods in electromagnetic hypersonics including fully coupled Maxwell Keplerlaan 1, 2200 AG Noordwijk, The Netherlands We describe a numerical technique for solving the coupled
D'Ambrosio, Domenic
Numerical Zoom for Multiscale Problems with an Application to Nuclear Waste Disposal
Numerical Zoom for Multiscale Problems with an Application to Nuclear Waste Disposal Jean of a nuclear waste repository site. Key words: Multiscale, Finite Element, Domain Decomposition, Chimera, Numerical Zoom, Nuclear Waste. PACS: 02.30.Jr, 47.11.Fg, 28.41.Kw, 47.55.P- 1 Introduction The present paper
An efficient numerical terrestrial scheme (ENTS) for fast earth system modelling
An efficient numerical terrestrial scheme (ENTS) for fast earth system modelling Mark Williamson Working Paper 83 #12;An efficient numerical terrestrial scheme (ENTS) for fast earth system modelling Mark for long time period simulations and large ensemble studies in Earth system models of intermediate
Williamson, Mark
Penetration of solar radiation in the upper ocean: A numerical model for oceanic and coastal waters
Penetration of solar radiation in the upper ocean: A numerical model for oceanic and coastal waters in the upper ocean, the vertical distribution of solar radiation (ESR) in the shortwave domain plays (2005), Penetration of solar radiation in the upper ocean: A numerical model for oceanic and coastal
Lee, Zhongping
USING LEARNING MACHINES TO CREATE SOLAR RADIATION MAPS FROM NUMERICAL WEATHER PREDICTION MODELS to develop a methodology to generate solar radiation maps using information from different sources. First with conclusions and next works in the last section. Keywords: Solar Radiation maps, Numerical Weather Predictions
Experimental characterization and numerical simulations of a syntactic-foam/glass-bre composite sandwich
Experimental characterization and numerical simulations of a syntactic-foam/glass-®bre composite core. Such core consists of a syntactic foam made by hollow glass microspheres embedded in an epoxy. Keywords: A. Glass ®bre; Composite sandwich; Syntactic foam; Mechanical tests; Numerical simulations (FE) 1
Corigliano, Alberto
Parallel Public-domain Numerical Libraries CSCS Tutorial Series, 11.05.2010
;Overview Introduction and library overview Typical =numerical and semi-numerical problems Survey) Preconditioning of large systems (PETSc) Eigenvalue / singular value decompositions of sparse/dense matrices: LU factorization: Forward/back-substitution: Eigenvalue decomposition: Generalized eigen
Numerical methods of integration applied in the nonlinear dynamic analysis of shells of revolution
structural bel&avior a&d therefore cannot be used for nonlinear analyses. lhe second nðod of dynamic analysis uses numerical integra- tion of the equations of notion wit'bout necessitating a transforma- tion of coordinates. Tne numerical integration...!riinear beam analysis {first vibratory mode only) with the most prc!! ising methods b ing programmed for shell analy! es. Rea- listic test problems are used to evaluate three of the numerical integration procedures as applied in shell of revolution analyses...
Tillerson, Joe Richard
Determination of SNe explosions frequency distribution function.Method and numerical simulations
The method for determination of the Supernovae (SNe) explosions frequency distribution function based on the assumption of explosions independence are offered. The method is based on assumption that the sequence of SNe explosions in an individual galaxy is a Poisson sequence. The essence of the method is in the determination of statistical moments of the frequency of the SNe explosions and subsequent determination of distribution function . The program of numerical simulation has been developed for testing the efficiency of the method. Numerical simulations show that even for a small mean number of registered SNe explosions, method allows restoring initial distribution function. The results of numerical simulations are given.
A. A. Akopian
Notes 06. Numerical solution of EOM for a single degree of freedom (SDOF) system
with fidelity and accuracy. The numerical integration relies in representing time derivatives of a function with an algebraic approximation, for example ()() 0 1 lim ~ t ii ii dx x x dt t xt t xt xx x tt ?? + ? == ? ? +? ? ? ?? null null null (13) Eq... values at a prior time step, i.e. ( ) 1 , ii ii xxtf xt + =+? , (15) MEEN 617 HD 6 Numerical Integration for Time Response: SDOF system L. San Andr?s ? 2008 5 while in an implicit numerical scheme ( ) 11 , ii ii x xtf xt ++ =+? (16...
Numerical construction of initial data for Einstein's equations with static extension to space-like infinity
We describe a numerical method to construct Cauchy data extending to space-like infinity based on Corvino's (2000) gluing method. Adopting the setting of Giulini and Holzegel (2005), we restrict ourselves here to vacuum axisymmetric spacetimes and glue a Schwarzschildean end to Brill-Lindquist data describing two non-rotating black holes. Our numerical implementation is based on pseudo-spectral methods, and we carry out extensive convergence tests to check the validity of our numerical results. We also investigate the dependence of the total ADM mass on the details of the gluing construction.
Georgios Doulis; Oliver Rinne
Neutrino Factory / Muon Collider Target Meeting Numerical Simulations for Jet-Proton
breakup observed in simulations Mercury is able to sustain very large tension Jet oscillates after on the surface of a hydrofoil Pressure contour in mercury target. #12;13 The Bubble Insertion Model Numerical
McDonald, Kirk
Direct numerical simulation of electroconvective instability and hysteretic current-voltage response of a permselective membrane
We present a systematic, multiscale, fully detailed numerical modeling for dynamics of fluid flow and ion transport covering Ohmic, limiting, and overlimiting current regimes in conductance of ion-selective membrane. By ...
Pham, Van Sang
Numerical simulations of ion transport membrane oxy-fuel reactors for CO? capture applications
Numerical simulations were performed to investigate the key features of oxygen permeation and hydrocarbon conversion in ion transport membrane (ITM) reactors. ITM reactors have been suggested as a novel technology to enable ...
Hong, Jongsup
Numerical studies of hypersonic binary gas-mixture flows near a sphere
Numerical studies of hypersonic binary gas-mixture flows near a sphere V.V. Riabov 1 Diffusive] to study the flow. In the present study, diffusive effects in hypersonic flows of binary gas-mixtures near
Riabov, Vladimir V.
Numerical simulation of three-dimensional combined convective radiative heat transfer in rectangular channels
This dissertation presents a numerical simulation of three-dimensional flow and heat transfer in a channel with a backward-facing step. Flow was considered to be steady, incompressible, and laminar. The flow medium was treated to be radiatively...
Ko, Min Seok
Simulation of surface waves with porous boundaries in a 2-D numerical wave tank
with the BEM solutions. Good agreement is observed between the two independent solutions. After verifying the numerical methods, we studied the interaction of water waves with a porous or rigid bottom-mounted half cylinder. The reflection and transmission...
Koo, Weoncheol
Numerical analysis and thermographic investigation of induction heating Matej Kranjc, Anze Zupanic *, Damijan Miklavcic, Tomaz Jarm
Thermography a b s t r a c t Induction heating process was investigated numerically and experimentally, a non-contact measurement method based on thermography can be used. Thermography is a form of infrared
Ljubljana, University of
Comprehensive Modeling and Numerical Investigation of Entrained-Flow Coal Gasifiers.
??Numerical simulations of coal gasification process inside a generic 2-stage entrainedflow gasifier are carried out using the commercial CFD solver ANSYS/FLUENT. The 3-D Navier-Stokes equations� (more)
Silaen, Armin
A numerical analysis of condenser performance of a seawater desalination system
This thesis presents the numerical analysis of three type condensers for desalination of seawater system. The condensers that were analyzed were a finned tube condenser that was built in Malaysia desalination plant, a ...
Mohamed, Hassan, S.B. Massachusetts Institute of Technology
Numerical and experimental analyses of resin infusion manufacturing processes of composite materials
Liquid resin infusion (LRI) processes are promising manufacturing routes to produce large, thick, or complex structural parts. They are based on the resin flow induced, across its thickness, by a pressure applied onto a preform/resin stacking. However, both thickness and fiber volume fraction of the final piece are not well controlled since they result from complex mechanisms which drive the transient mechanical equilibrium leading to the final geometrical configuration. In order to optimize both design and manufacturing parameters, but also to monitor the LRI process, an isothermal numerical model has been developed which describes the mechanical interaction between the deformations of the porous medium and the resin flow during infusion.1, 2 With this numerical model, it is possible to investigate the LRI process of classical industrial part shapes. To validate the numerical model, first in 2D, and to improve the knowledge of the LRI process, this study details a comparison between numerical simulations and...
Wang, Peng; Molimard, J�r�me; Vautrin, Alain; Minni, Jean-Christophe; 10.1177/0021998311421990
Thermoelectric Generators for Automotive Waste Heat Recovery Systems Part I: Numerical Modeling
Thermoelectric Generators for Automotive Waste Heat Recovery Systems Part I: Numerical Modeling telluride TEMs. Key words: Thermoelectric generators, waste heat recovery, automotive exhaust, skutterudites bismuth telluride are considered for thermoelectric modules (TEMs) for conversion of waste heat from
Xu, Xianfan
Analytical and numerical Gubser solutions of the second-order hydrodynamics
Evolution of quark-gluon plasma (QGP) near equilibrium can be described by the second-order relativistic viscous hydrodynamic equations. Consistent and analytically verifiable numerical solutions are critical for phenomenological studies of the collective behavior of QGP in high-energy heavy-ion collisions. A novel analytical solution based on the conformal Gubser flow which is a boost-invariant solution with transverse fluid velocity is presented. It is used to verify with high precision the numerical solution with a newly developed $(3+1)$-dimensional second-order viscous hydro code (CLVisc). The perfect agreement between the analytical and numerical solutions demonstrates the reliability of the numerical simulations with the second-order viscous corrections. This lays the foundation for future phenomenological studies that allow one to gain access to the second-order transport coefficients.
Long-Gang Pang; Yoshitaka Hatta; Xin-Nian Wang; Bo-Wen Xiao
Some aspects of the numerical treatment of the conformal field equations
Some issues in the numerical treatment of the conformal field equations are discussed. In particular, the problem of obtaining smooth initial data for the hyperboloidal initial value problem is described and solution methods are presented.
J. Frauendiener
Numerical Analysis of Heat and Moisture Transfer in Underground Air-conditioning Systems
In view of the influence of humidity of room air on room heat load, indoor environment and building energy consumption in underground intermittent air-conditioning systems, numerical simulation was used to dynamically analyze the coupling condition...
Wang, Q.; Miao, X.; Cheng, B.; Fan, L.
Comparison of Numerical and Analytical Solutions of Heat Equation in R1n
Example Analytical Calculation of Heat Equation in 1-D 4 4 2-D Matrix Methods for Numerical Approximations Bibliography 11 8 Appendices 11 8.1 Graphs . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11 8.2 Matlab Code
Marzuola, Jeremy
2D-Modelling of pellet injection in the poloidal plane: results of numerical tests
2D-Modelling of pellet injection in the poloidal plane: results of numerical tests P. Lalousis developed for computing the expansion of pellet-produced clouds in the poloidal plane. The expansion
Experimental and Numerical Investigation of Pool Boiling Heat Transfer on Engineered Nano-Finned Surfaces
Pool boiling experiments for nanocoatings (or nanostructured surfaces) show that despite the lower thermal conductivity values than carbon, silicon yielded higher values of CHF (critical heat flux). Subsequently numerical studies showed...
Yang, Hongjoo
Multi-scale problems, high performance computing and hybrid numerical methods
Multi-scale problems, high performance computing and hybrid numerical methods G. Balarac, G of High Performance Computing G. Balarac LEGI, CNRS and Universit´e de Grenoble, BP 53, 38041 Grenoble
Cottet, Georges-Henri
Multi-scale problems, high performance computing and hybrid numerical methods G. Balarac, G of High Performance Computing (HPC) is not anymore restricted to academia and scientific grand challenges
CAVERN SOULEY Mountaka1 , MERCERAT Diego2 , DRIAD-LEBEAU Lynda1 , BERNARD Pascal2 1 INERIS, Ecole des collapse). KEYWORDS: cavern, numerical modelling, continuum-discrete, overburden, damage. R�?SUM�?: Dans l
On the numerical simulation of the instationary free fall of a solid in a fluid.
from meteorology, sedimentology and aerospace engineering to biology. We present a new numerical method], sedimentology [1], aerospace engineering [21] and biology (e.g. models for animal flights [2]). Recently
A RANS/DES Numerical Procedure for Axisymmetric Flows with and without Strong Rotation
A RANS/DES numerical procedure with an extended Lax-Wendroff control-volume scheme and turbulence model is described for the accurate simulation of internal/external axisymmetric flow with and without strong rotation. This new procedure is an extension, from Cartesian to cylindrical coordinates, of (1) a second order accurate multi-grid, control-volume integration scheme, and (2) a k-{omega} turbulence model. This paper outlines both the axisymmetric corrections to the mentioned numerical schemes and the developments of techniques pertaining to numerical dissipation, multi-block connectivity, parallelization, etc. Furthermore, analytical and experimental case studies are presented to demonstrate accuracy and computational efficiency. Notes are also made toward numerical stability of highly rotational flows.
Andrade, A J
Numerical Spherically Symmetric Static Solution of the RTG Equations Outside the Matter
There was obtained a numerical external solution for the exact system of the RTG equations with some natural boundary conditions in the static spherically symmetric case. The properties of the solution are discussed.
A. Godizov
Numerical simulation of the impeller tip clearance effect on centrifugal compressor performance
This thesis presents the numerical simulation of flow in centrifugal compressors. A three-dimensional Navier-Stokes solver was employed to simulate flow through two centrifugal compressors. The first compressor simulated was the NASA low speed...
Hoenninger, Corbett Reed
American Institute of Aeronautics and Astronautics Numerical Simulation of a Gas Turbine Combustor Using
1 American Institute of Aeronautics and Astronautics Numerical Simulation of a Gas Turbine of combustion by using nanosecond pulsed plasma actuators for a gas turbine combustor. Moreau [2] and Corke et
of electrical machines is the accurate computation of the power losses in the ferromagnetic components numerical results which allow us to confirm the theoretical estimates and to assess the performance al.,22,23 which consists in computing the electr
Rodríguez, Rodolfo
A theoretical and numerical procedure for predicting sailing yacht lift and drag
In this thesis, a theoretical and numerical procedure for predicting the effects of viscosity on the hydrodynamic forces developed by a sailing yacht hull is presented. A simultaneous viscous/inviscid algorithm is developed ...
Cairoli, Claudio, 1975-
Numerical study of flow and heat transfer in 3D serpentine channels using colocated grids
and average Nusselt number. The numerical code developed was validated by solving for fully developed flow and heat transfer in a square straight channel. Grid-independent solution was established for a reference case of serpentine channel with the highest...
Chintada, Sailesh Raju
Laboratory and numerical investigation of transport processes occurring above and within a saltwater wedge
Laboratory and numerical investigation of transport processes occurring above and within recent studies have investigated saltwater transport in coastal aquifers and the associated flow wedge itself or on studying contaminant transport processes occurring above the wedge. As per our
Acceleration of heavy and light particles in turbulence: comparison between experiments and direct numerical simulations
We compare experimental data and numerical simulations for the dynamics of inertial particles with finite density in turbulence. In the experiment, bubbles and solid particles are optically tracked in a turbulent flow of water using an Extended Laser Doppler Velocimetry technique. The probability density functions (PDF) of particle accelerations and their auto-correlation in time are computed. Numerical results are obtained from a direct numerical simulation in which a suspension of passive pointwise particles is tracked, with the same finite density and the same response time as in the experiment. We observe a good agreement for both the variance of acceleration and the autocorrelation timescale of the dynamics; small discrepancies on the shape of the acceleration PDF are observed. We discuss the effects induced by the finite size of the particles, not taken into account in the present numerical simulations.
R. Volk; E. Calzavarini; G. Verhille; D. Lohse; N. Mordant; J. -F. Pinton; F. Toschi
Modeling and numerics for two partial differential equation systems arising from nanoscale physics
), we study the operation of the device in several specific asymptotic regimes. Furthermore, we simulate such devices using a customized 2D hybrid discontinuous Galerkin finite element scheme and compare the numerical results to our asymptotics. Next, we...
Brinkman, Daniel
Numerical Simulation of Enhanced Mixing in Scramjet Combustor Using Ramp, Tabs and Suction Collar
Numerical simulations of the scramjet combustor by using the commercial CFD code Fluent with the coupled implicit method with second-order accurate discretization have been obtained for the reacting flows with the parallel fuel injection (ramp...
Hwang, Seung-Jae
Three-dimensional numerical manifold method simulations for blocky rock analysis
After decades of development, people realize that there are wider and more various applications of numerical modeling and analysis. However, current feasible software tools cannot satisfy engineering and commercial goals. ...
Shentu, Longfei
Numerical and analytical studies of single and multiphase starting jets and plumes
Multiphase starting jets and plumes are widely observed in nature and engineering systems. An environmental engineering example is open-water disposal of sediments. The present study numerically simulates such starting ...
Wang, Ruo-Qian
Consolidation theories for saturated-unsaturated soils and numerical simulation of residential buildings on expansive soils
and numerical techniques: different simulation methods for different boundary conditions such as tree, grass, and bare soils, coupled hydro-mechanical stress analysis to describe deformation of saturated-unsaturated soils, jointed elements simulation of soil...
Zhang, Xiong
Non-smooth Dynamics Using Differential-algebraic Equations Perspective: Modeling and Numerical Solutions
mathematical tools. On the other hand, the approach based on differential-algebraic equations gives more insight into the constitutive assumptions of a chosen model and easier to obtain numerical solutions. Bingham-type models in which the force cannot...
Gotika, Priyanka
Multiscale numerical methods for partial differential equations using limited global information and their applications
. These methods employ limited single or multiple global information. We apply these numerical methods to partial differential equations (elliptic, parabolic and wave equations) with continuum scales. To compute the solution of partial differential equations on a...
Jiang, Lijian
, vibration/noise control, and fluid mechanics, simultaneously. Higher product quality, less developing time and lower manufacturing cost will be achieved through a balanced and organic MDO method. In this paper, numerical stress analysis, optimization method...
Oh, Bong T.
Numerical Investigation of Temperature Distribution on a High Pressure Gas Turbine Blade
A numerical code is developed to calculate the temperature distributions on the surface of a gas turbine blade. This code is a tool for quick prediction of the temperatures by knowing the boundary conditions and the flow conditions, and doesn...
Zirakzadeh, Hootan
Numerical simulation of the response of sandy soils treated with pre-fabricated vertical drains
This research is part of the ongoing effort of the Seismic Risk Mitigation for Port Systems Grand Challenge. It addresses the problem of numerically simulating the response of sandy soils treated with earthquake drains, ...
Vytiniotis, Antonios
An investigation of analytical and numerical sucker rod pumping mathematical models
AN INVESTIGATION OF ANALYTICAL AND NUMERICAL SUCKER ROD PUMPING MATHEMATICAL MODELS A Thesis by DONALD JOSEPH SCHAFER Submitted to the Graduate College of Texas A&M University in partial fulfillment of the requirements for the degree... of MASTER OF SCIENCE May 1987 Major Subject: Petroleum Engineering AN INVESTIGATION OF ANALYTICAL AND NUMERICAL SUCKER ROD PUMPING MATHEMATICAL MODELS A Thesis by DONALD JOSEPH SCHAFER Approved as to style and content by: 7d JW. J ni (Chai...
Schafer, Donald Joseph
Numerical analysis of turbulent heat transfer in a nuclear reactor coolant channel
NUMERICAL ANALYSIS OF TURBULENT HEAT TRANSFER IN A NUCLEAR REACTOR COOLANT CHANNEL A Thesis Clarence William Garrard, Jr. Submitted to the Graduate College of the Texas A&M University in partial fulfillment of' the requirements for the degree... of' MASTER OF SC1ENCE May, 1965 Ma)or Subject Nuclear Engineering NUMERICAL ANALYSIS OF TURBULENT HEAT TRANSFER 1N A NUCLEAR REACTOR COOLANT CHANNEL A Thesis By Clarence William Garrard, Jr. Approved as to style and content by; Head...
Garrard, Clarence William
Numerical Simulation of Cold Pressing of Armstrong CP-Ti Powders
Numerical simulation results for the cold pressing of Armstrong CP-Ti Powders are presented. The computational model was implemented in the commercial finite element program ABAQUSTM. Several simulation cases were conducted for cylindrical samples with different friction coefficients and different compaction pressures, under both single-action and dual-action uniaxial pressing. Numerical simulation results for the density distribution are compared against experimental data in order to validate the computational model.
Sabau, Adrian S [ORNL] [ORNL; Gorti, Sarma B [ORNL] [ORNL; Peter, William H [ORNL] [ORNL; Chen, Wei [ORNL] [ORNL; Yamamoto, Yukinori [ORNL] [ORNL
Comparative time studies and error analyses of certain numerical methods of solving differential equations
, to whom the writer is deeply' indebted for her help and patient understanding. TABLE OF CONTENTS I. INTRODUCTION II. DEVELOPMENT OF CERTAIN NUMERICAL METHODS OF SOLVING DIFFERENTIAL EQUATIONS OF THE TYPE y' = f(x, y) Definitions Abbreviations Runga... 39 41 INTRODUCTION The purpose of this thesis is to study the aocuracy and calculation time of various numerical methods of solving first order differential equations. The ob]ective is to gain from actual computations an insight into the method...
Steger, Michael John
Numerical solutions of a system of partial differential equations by a linear programming technique
NUMERICAL SOLUTIONS OF A SYSTEM OF PARTIAL DIFFERENTIAL EQUATIONS BY A LINEAR PROGRAMMING TECHNIQUE A Thesis By PATRICK LOUIS MUSMAKER Submitted to the Graduate College of the Texas A&M University in partial fulfillment of the requirements... for the degree of MASTER OF SCIENCE May 1967 Major Subject: Mathematics NUMERICAL SOLUTIONS OF A SYSTEM OF PARTIAL DIFFERENTIAL EQUATIONS BY A LINEAR PROGRAMMING TECHNIQUE A Thesis By PATRICK LOUIS MUSMAKER Approved as to sty1e and content by: (Chairman...
Musmaker, Patrick Louis
A robotic crawler exploiting directional frictional interactions: experiments, numerics, and derivation of a reduced model
We present experimental and numerical results for a model crawler which is able to extract net positional changes from reciprocal shape changes, i.e. 'breathing-like' deformations, thanks to directional, frictional interactions with a textured solid substrate, mediated by flexible inclined feet. We also present a simple reduced model that captures the essential features of the kinematics and energetics of the gait, and compare its predictions with the results from experiments and from numerical simulations.
Giovanni Noselli; Antonio DeSimone
Numerical simulation of comminution in granular materials with an application to fault gouge evolution
NUMERICAL SIMULATION OF COMMINUTION IN GRANULAR MATERIALS WITH AN APPLICATION TO FAULT GOUGE EVOLUTION A Thesis by RICHARD ANTHONY LANG Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment... of the requirements for the degree of MASTER OF SCIENCE May 2002 Major Subject: Geophysics NUMERICAL SIMULATION OF COMMINUTION IN GRANULAR MATERIALS WITH AN APPLICATION TO FAULT GOUGE EVOLUTION A Thesis by RICHARD ANTHONY LANG Submitted...
Lang, Richard Anthony
Numerical simulation of scour process around bridge piers in cohesive soil
of streambed materials is also necessary. Consequently, progresses are very limitted in this area. Numerical computation of bridge pier scour in noncohesive soils can only be found in Olsen and Melaaen (1993) and Olsen (1996), and no reference exists... in numerical simulations of scour process. Furthermore, all the references are confined to noncohesive soils. Olsen and Melaaen (1993) calculated the shape of a scour hole around a cylinder mounted in a noncohesive streambed. They used a three...
Wei, Gengsheng
A numerical study of the baroclinic circulation in the Gulf of Maine
A NUMERICAL STUDY OF THE BAROCLINIC CIRCULATION IN THE GULF OF MAINE A Thesis JOSEPH CLAYTON BOGUSKY Submitted to the Graduate College of Texas A&M University in partial fulfillment of the requirements for the degree of MASTER OF SCIENCE... AUGUST 1985 Major Subject: Oceanography A NUMERICAL STUDY OF THE BAROCLINIC CIRCULATION IN THE GULF OF MAINE A Thesis by JOSEPH CLAYTON BOGUSXY Approved as to style and content by: David A. Brooks (Chairman) David C. Smith IV (Member) John M...
Bogusky, Joseph Clayton
A numerical study of seasonal circulation in the Gulf of Cariaco
A NUMERICAL STUDY OF SEASONAL CIRCULATION ZN THE GULF OF CARIACO A Thesis ANNICX JEAN LOPEZ-GARCIA Submitted to the Graduate College of Texas ARM University in partial fulfillment of the requirements for the degree of MASTER OF SCIENCE... December 1984 Major Subject: Oceanography A NUMERICAL STUDY OF SEASONAL CIRCULATION IN THE GULF OF CARIACO A Thesis by ANNICK JEAN LOPEZ-GARCIA approved as to style and content by: JOHN M. KLINCK (Chairman of Committee) ANDREW C. VA ANO (Member...
Lopez-Garcia, Annick Jean
Experimental and numerical investigation of one-dimensional waterflood in porous reservoir
Experimental and numerical investigation of relative permeability and oil recovery from the porous reservoir are described for short and long core samples. The relative permeability ratios, which are function of water saturation, obtained from laboratory core flooding experiments have been used for prediction of oil recovery through numerical simulation of non-dimensional Buckley-Leverett equation. The simulation results for oil recovery compared well with recovery results obtained from core flooding experiments. (author)
Hadia, N.; Chaudhari, L.; Mitra, Sushanta K. [IITB-ONGC Joint Research Centre, and Department of Mechanical Engineering, Indian Institute of Technology Bombay, Mumbai 400 076 (India); Aggarwal, A.; Vinjamur, M. [IITB-ONGC Joint Research Centre, and Department of Chemical Engineering, Indian Institute of Technology Bombay, Mumbai 400 076 (India); Singh, R. [IITB-ONGC Joint Research Centre, and Institute of Reservoir Studies, ONGC, Ahmedabad (India)
First principles predictions of intrinsic defects in aluminum arsenide, AlAs : numerical supplement.
This Report presents numerical tables summarizing properties of intrinsic defects in aluminum arsenide, AlAs, as computed by density functional theory. This Report serves as a numerical supplement to the results published in: P.A. Schultz, 'First principles predictions of intrinsic defects in Aluminum Arsenide, AlAs', Materials Research Society Symposia Proceedings 1370 (2011; SAND2011-2436C), and intended for use as reference tables for a defect physics package in device models.
Numerical simulation of large amplitude liquid sloshing in a rigid rectangular tank
NUMERICAL SIMULATION OF LARGE AMPLITUDE LIQUID SLOSHING IN A RIGID RECTANGULAR TANK A Thesis by THOMAS JACKSON BRIDGES Submitted to the Graduate College of Texas A8M University in partial fulfillment of the requirement for the degree... of MASTER OF SCIE'ICE December 1981 Major Subject: Ocean Engineering NUMERICAL SIMULATION OF LARGE AMPLITUDE LIQUID SLOSHING IN A RIGID RECTANGULAR TANK A Thesis by THOMAS JACKSON BRIDGES Approved as to sty1e and content by: (Chairman of Commit ee...
Bridges, Thomas J.
A numerical model of convective heat transfer in a three dimensional channel with baffles
to minimize the effects of the temperature. Chandrupatla and Sastri (1977) used a finite difference method to study laminar heat transfer and fluid flow for Non-Newtonian fluids. In the limiting case of Newtonian fluids their results showed good agreement... and Sastri (1977) for developing laminar flow in a rectangular smooth channel. Figure 4. 1 shows the numerically predicted centerline axial velocity compared with the experimental data of Goldstein and Kreid (1967). The agreement between the numerical...
Lopez Buso, Jorge Ricardo
Numerical modeling of geothermal systems with applications to Krafla, Iceland and Olkaria, Kenya
The use of numerical models for the evaluation of the generating potential of high temperature geothermal fields has increased rapidly in recent years. In the present paper a unified numerical approach to the modeling of geothermal systems is discussed and the results of recent modeling of the Krafla geothermal field in Iceland and the Olkaria, Kenya, are described. Emphasis is placed on describing the methodology using examples from the two geothermal fields.
Bodvarsson, G.S.
Force-controlled absorption in a fully-nonlinear numerical wave tank
An active control methodology for the absorption of water waves in a numerical wave tank is introduced. This methodology is based upon a force-feedback technique which has previously been shown to be very effective in physical wave tanks. Unlike other methods, an a-priori knowledge of the wave conditions in the tank is not required; the absorption controller being designed to automatically respond to a wide range of wave conditions. In comparison to numerical sponge layers, effective wave absorption is achieved on the boundary, thereby minimising the spatial extent of the numerical wave tank. In contrast to the imposition of radiation conditions, the scheme is inherently capable of absorbing irregular waves. Most importantly, simultaneous generation and absorption can be achieved. This is an important advance when considering inclusion of reflective bodies within the numerical wave tank. In designing the absorption controller, an infinite impulse response filter is adopted, thereby eliminating the problem of non-causality in the controller optimisation. Two alternative controllers are considered, both implemented in a fully-nonlinear wave tank based on a multiple-flux boundary element scheme. To simplify the problem under consideration, the present analysis is limited to water waves propagating in a two-dimensional domain. The paper presents an extensive numerical validation which demonstrates the success of the method for a wide range of wave conditions including regular, focused and random waves. The numerical investigation also highlights some of the limitations of the method, particularly in simultaneously generating and absorbing large amplitude or highly-nonlinear waves. The findings of the present numerical study are directly applicable to related fields where optimum absorption is sought; these include physical wavemaking, wave power absorption and a wide range of numerical wave tank schemes.
Spinneken, Johannes, E-mail: email@example.com; Christou, Marios; Swan, Chris
Numeric simulation of thyristor power converters taking into consideration the firing-circuits and the closed-loop control
Numeric simulation of thyristor power converters taking into consideration the firing-circuits and the closed-loop control
Bordry, Frederick; Proudlock, Paul
A numerical model of hydrothermal cooling and crustal accretion at a fast spreading mid-ocean ridge
: 7360 words, 11 figures. Keywords: Mid-ocean ridge processes; hydrothermal cooling; numerical modelA numerical model of hydrothermal cooling and crustal accretion at a fast spreading mid-ocean ridge present a steady state numerical representation of the sill model that explicitly includes hydrothermal
Toomey, Doug
ATS 680 A6: Applied Numerical Weather Prediction MW, 1:00-1:50 PM, ACRC Room 212B
experiments using a state-of-the-art numerical weather prediction model · Discuss the strengths and weaknesses, Parameterization Schemes: Keys to Understanding Numerical Weather Prediction Models, Cambridge University PressATS 680 A6: Applied Numerical Weather Prediction Fall 2013 MW, 1:00-1:50 PM, ACRC Room 212B Course
Estimation of the mean depth of boreal lakes for use in numerical weather prediction and climate modelling
in the numerical weather prediction (NWP) and climate models through parameterisation. For parameterisation, data. The effect of lakes should be parameterised in numerical weather prediction (NWP) and climate modellingEstimation of the mean depth of boreal lakes for use in numerical weather prediction and climate
Metr 6803: Numerical Weather Prediction Syllabus: Spring 2013 M / W 10:00-11:15, NWC 5930
Metr 6803: Numerical Weather Prediction Syllabus: Spring 2013 M / W 10:00-11:15, NWC 5930 weather analysis (NWA) and numerical weather prediction (NWP)? - why are they important? - how "good of numerics of GFS, RUC/RAP, CAPS, MM5, WRF models 11. Atmospheric Predictability - basic concepts
Droegemeier, Kelvin K.
Energy conservation in the numerical solution of Hamiltonian boundary value Pierluigi Amodio, Luigi Brugnano, and Felice Iavernaro
Energy conservation in the numerical solution of Hamiltonian boundary value problems Pierluigi://scitationnew.aip.org/termsconditions. Downloaded to IP: 22.214.171.124 On: Mon, 21 Oct 2013 16:20:27 #12;Energy Conservation in the Numerical Firenze, Italy Abstract. We consider the issue of energy conservation in the numerical solution
Brugnano, Luigi
Context DDFV Scheme Numerical analysis ECG Simulation 2D/3D DDFV scheme for anisotropic-heterogeneous
Context DDFV Scheme Numerical analysis ECG Simulation 2D/3D DDFV scheme for anisotropic Numerical analysis ECG Simulation Authors : Laboratoire de math´ematiques Jean Leray, Universit´e de Nantes Pays de l'Adour : · Charles Pierre #12;Context DDFV Scheme Numerical analysis ECG Simulation Outline
Pierre, Charles
Study on Applicability of Numerical Simulation to Evaluation of Gas Entrainment From Free Surface
An onset condition of gas entrainment (GE) due to free surface vortex has been studied to establish a design of fast breeder reactor with higher coolant velocity than conventional designs, because the GE might cause the reactor operation instability and therefore should be avoided. The onset condition of the GE has been investigated experimentally and theoretically, however, dependency of the vortex type GE on local geometry configuration of each experimental system and local velocity distribution has prevented researchers from formulating the universal onset condition of the vortex type GE. A real scale test is considered as an accurate method to evaluate the occurrence of the vortex type GE, but the real scale test is generally expensive and not useful in the design study of large and complicated FBR systems, because frequent displacement of inner equipments accompanied by the design change is difficult in the real scale test. Numerical simulation seems to be promising method as an alternative to the real scale test. In this research, to evaluate the applicability of the numerical simulation to the design work, numerical simulations were conducted on the basic experimental system of the vortex type GE. This basic experiment consisted of rectangular flow channel and two important equipments for vortex type GE in the channel, i.e. vortex generation and suction equipments. Generated vortex grew rapidly interacting with the suction flow and the grown vortex formed a free surface dent (gas core). When the tip of the gas core or the bubbles detached from the tip of the gas core reached the suction mouth, the gas was entrained to the suction tube. The results of numerical simulation under the experimental conditions were compared to the experiment in terms of velocity distributions and free surface shape. As a result, the numerical simulation showed qualitatively good agreement with experimental data. The numerical simulation results were similar to the experimental results in terms of the shape of free surface dent and the velocity distribution around the vortex, although the GE itself was not completely reproduced due to a lack of enough mesh partition. After confirming the applicability of the numerical simulation to the GE evaluation, several parameters, such as suction velocity and a configuration around suction mouth, were numerically examined to evaluate their influence on the GE. The tendencies of the GE occurrence enhanced by larger suction velocity or suction mouth on bottom surface enhance occurrence of the GE were obtained from the simulation results. These simulation results implied that the numerical simulation has enough potential to be used for the design work. (authors)
Kei Ito; Takaaki Sakai; Hiroyuki Ohshima [Japan Atomic Energy Agency (Japan)
Subtask 2.2 - Creating A Numerical Technique for Microseismic Data Inversion
Geomechanical and geophysical monitoring are the techniques which can complement each other and provide enhancement in the solutions of many problems of geotechnical engineering. One of the most promising geophysical techniques is passive seismic monitoring. The essence of the technique is recording the acoustic signals produced in the subsurface, either naturally or in response to human activity. The acoustic signals are produced by mechanical displacements on the contacts of structural elements (e.g., faults, boundaries of rock blocks, natural and induced fractures). The process can be modeled by modern numerical techniques developed in geomechanics. The report discusses a study that was aimed at the unification of the passive seismic monitoring and numerical modeling for the monitoring of the hydraulic fracture propagation. The approach adopted in the study consisted of numerical modeling of the seismicity accompanying hydraulic fracture propagation and defining seismic attributes and patterns characterizing the process and fracture parameters. Numerical experiments indicated that the spatial distribution of seismic events is correlated to geometrical parameters of hydrofracture. Namely, the highest density of the events is observed along fracture contour, and projection of the events to the fracture plane makes this effect most pronounced. The numerical experiments also showed that dividing the totality of the events into groups corresponding to the steps of fracture propagation allows for reconstructing the geometry of the resulting fracture more accurately than has been done in the majority of commercial applications.
Anastasia Dobroskok; Yevhen Holubnyak; James Sorensen
The RELAP-7 code is the next generation of nuclear reactor system safety analysis code being developed at the Idaho National Laboratory (INL). All the physics in RELAP-7 are fully coupled and the errors resulted from the traditional operator-splitting approach are eliminated. By using 2nd order methods in both time and space and eliminating operator-splitting errors, the numerical error of RELAP-7 can be minimized. Numerical verification is the process to verify the orders of numerical methods. It is an important part of modern verification and validation process. The core channel component in RELAP-7 is designed to simulate coolant flow as well as the conjugated heat transfer between coolant flow and the fuel rod. A special treatment at fuel centerline to avoid numerical singularity for the cylindrical heat conduction in the continuous finite element mesh is discussed. One steady state test case and one fast power up transient test case are utilized for the verification of the core channel model with single-phase flow. Analytical solution for the fuel pin temperature and figures of merit such as peak clad temperature and peak fuel temperature are used to define numerical errors. These cases prove that the mass and energy are well conserved and 2nd order convergence rates for both time and space are achieved in the core channel model.
Haihua Zhao; Ling Zou; Hongbin Zhang; Richard Martineau
Numerical and experimental investigation of melting with internal heat generation within cylindrical enclosures
There have been significant efforts by the heat transfer community to investigate the melting phenomenon of materials. These efforts have included the analytical development of equations to represent melting, numerical development of computer codes to assist in modeling the phenomena, and collection of experimental data. The understanding of the melting phenomenon has application in several areas of interest, for example, the melting of a Phase Change Material (PCM) used as a thermal storage medium as well as the melting of the fuel bundle in a nuclear power plant during an accident scenario. The objective of this research is two-fold. First a numerical investigation, using computational fluid dynamics (CFD), of melting with internal heat generation for a vertical cylindrical geometry is presented. Second, to the best of authors knowledge, there are very limited number of engineering experimental results available for the case of melting with Internal Heat Generation (IHG). An experiment was performed to produce such data using resistive, or Joule, heating as the IHG mechanism. The numerical results are compared against the experimental results and showed favorable correlation. Uncertainties in the numerical and experimental analysis are discussed. Based on the numerical and experimental analysis, recommendations are made for future work.
Amber Shrivastava; Brian Williams; Ali S. Siahpush; Bruce Savage; John Crepeau
Numerical study on the turbulence structures in closely spaced rod bundle subchannels
Fully developed turbulent flow through simulated rod bundle subchannels formed by a rod-trapezoidal duct was numerically studied. With a simple coordinate system transformation from an orthogonal cylindrical system to a nonorthogonal curvilinear system, the highly irregular flow passage of a rod-trapezoidal duct was converted to that of a regular rectangle. An empirical anisotropic eddy viscosity distribution based on existing experimental data was used in conjunction with the algebraic stress model to address the influence of coherent large-scale cross-gap eddy motion, whose existence in closely spaced rod bundle subchannels has bene substantiated by extensive hot-wire measurements. Results of the calculation are compared with experimental data, with emphasis on secondary flow and turbulence kinetic energy. The credibility of this numerical scheme was establishment through a series of numerical tests on simple geometry flows.
Wu, X. (Univ. of Manitoba, Winnipeg (Canada). Dept. of Mechanical Engineering)
Effect of virtual mass on the characteristics and the numerical stability in two-phase flows
It is known that the typical six equation two-fluid model of the two-phase flow possesses complex characteristics, exhibits unbounded instabilities in the short-wavelength limit and constitutes an ill-posed initial value problem. Among the suggestions to overcome these difficulties, one model for the virtual mass force terms were studied here, because the virtual mass represents real physical effects to accomplish the dissipation for numerical stability. It was found that the virtual mass has a profound effect upon the mathematical characteristic and numerical stability. Here a quantitative bound on the coefficient of the virtual mass terms was suggested for mathematical hyperbolicity and numerical stability. It was concluded that the finite difference scheme with the virtual mass model is restricted only by the convective stability conditions with the above suggested value.
No, H.C.; Kazimi, M.S.
Foundation heat exchangers for residential ground source heat pump systems Numerical modeling and experimental validation
A new type of ground heat exchanger that utilizes the excavation often made for basements or foundations has been proposed as an alternative to conventional ground heat exchangers. This article describes a numerical model that can be used to size these foundation heat exchanger (FHX) systems. The numerical model is a two-dimensional finite-volume model that considers a wide variety of factors, such as soil freezing and evapotranspiration. The FHX numerical model is validated with one year of experimental data collected at an experimental house located near Oak Ridge, Tennessee. The model shows good agreement with the experimental data-heat pump entering fluid temperatures typically within 1 C (1.8 F) - with minor discrepancies due to approximations, such as constant moisture content throughout the year, uniform evapotranspiration over the seasons, and lack of ground shading in the model.
Xing, Lu [Oklahoma State University; Cullin, James [Oklahoma State University; Spitler, Jeffery [Oklahoma State University; Im, Piljae [ORNL; Fisher, Daniel [Oklahoma State University
Focused particle transport in a nonuniform large-scale magnetic field is investigated numerically in the case of isotropic pitch-angle scattering. Evolving particle density profiles and distribution moments are computed from solutions of a system of stochastic differential equations, equivalent to the original Fokker-Planck equation for the particle distribution. Conflicting analytical predictions for the transport coefficients in the diffusion limit, independently calculated by Beeck and Wibberenz and Shalchi, are compared with the numerical results. The reasons for the discrepancies among the analytical and numerical treatments, as well as the general limitations of the diffusion model, are discussed. The telegraph equation, derived in a higher-order expansion of the particle distribution function, is shown to describe the particle transport much more accurately than the diffusion model, especially ahead of a moving density pulse.
Litvinenko, Yuri E. [Department of Mathematics, University of Waikato, P. B. 3105, Hamilton (New Zealand)] [Department of Mathematics, University of Waikato, P. B. 3105, Hamilton (New Zealand); Noble, P. L. [Sydney Institute for Astronomy, School of Physics, The University of Sydney, NSW 2006 (Australia)] [Sydney Institute for Astronomy, School of Physics, The University of Sydney, NSW 2006 (Australia)
Time-resolved temperature measurement and numerical simulation of millisecond laser irradiated silicon
Thermal process of 1064 nm millisecond pulsed Nd:YAG laser irradiated silicon was time-resolved temperature measured by an infrared radiation pyrometer, temperature evolutions of the spot center for wide range of laser energy densities were presented. The waveforms of temperature evolution curves contained much information about phase change, melting, solidification and vaporization. An axisymmetric numerical model was established for millisecond laser heating silicon. The transient temperature fields were obtained by using the finite element method. The numerical results of temperature evolutions of the spot center are in good agreement with the experimental results. Furthermore, the axial temperature distributions of the numerical results give a better understanding of the waveforms in the experimental results. The melting threshold, vaporizing threshold, melting duration, and melting depth were better identified by analyzing two kinds of results.
Li Zewen; Zhang Hongchao; Shen Zhonghua; Ni Xiaowu [School of Science, Nanjing University of Science and Technology, Nanjing 210094 (China)
Analytical formula for numerical evaluations of the Wigner rotation matrices at high spins
The Wigner d function, which is the essential part of an irreducible representation of SU(2) and SO(3) parameterized with Euler angles, has been know to suffer from a serious numerical errors at high spins, if it is calculated by means of the Wigner formula as a polynomial of cos and sin of half of the second Euler angle. This paper shows a way to avoid this problem by expressing the d functions as the Fourier series of the half angle. A precise numerical table of the coefficients of the series is obtainable from a web site.
Naoki Tajima
Numerical Solutions of the Einstein-Yang-Mills System with Cosmological Constant
Numerical evidence for a cosmological version of the Bartnik-McKinnon family of particle-like solutions of the Einstein-Yang-Mills system is presented. Our solutions are also static, but space has the topology of a three-sphere. By adjusting the cosmological constant we found numerically a spherically symmetric solution which can be regarded as an excitation of the unique SO(4)-invariant solution. We expect that for each node number there exists such a solution without a cosmological horizon.
P. Molnar
On the efficient numerical solution of lattice systems with low-order couplings
We apply the Quasi Monte Carlo (QMC) and recursive numerical integration methods to evaluate the Euclidean, discretized time path-integral for the quantum mechanical anharmonic oscillator and a topological quantum mechanical rotor model. For the anharmonic oscillator both methods outperform standard Markov Chain Monte Carlo methods and show a significantly improved error scaling. For the quantum mechanical rotor we could, however, not find a successful way employing QMC. On the other hand, the recursive numerical integration method works extremely well for this model and shows an at least exponentially fast error scaling.
A. Ammon; A. Genz; T. Hartung; K. Jansen; H. Le�vey; J. Volmer
On numerical simulation of flow, heat transfer and combustion processes in tangentially-fired furnace
In this work, an Eulerian/Lagrangian approach has been employed to investigate numerically flow characteristics, heat transfer and combustion processes inside corner-fired power plant boiler furnace. To avoid pseudo-diffusion that is significant in modeling tangentially-fired furnaces, some attempts have been made at improving the finite-difference scheme. Comparisons have been made between standard {kappa}-{epsilon} model and RNG {kappa}-{epsilon} model. Some new developments on turbulent diffusion of particles are taken into account in an attempt to improve computational accuracy. Finally, temperature deviation is studied numerically so as to gain deeper insight into tangentially fired furnace.
Sun, P.; Fan, J.; Cen, K.
Tajima, Naoki
Heat transfer coefficients in two-dimensional Yukawa systems (numerical simulations)
New data on heat transfer in two-dimensional Yukawa systems have been obtained. The results of a numerical study of the thermal conductivity for equilibrium systems with parameters close to the conditions of laboratory experiments in dusty plasma are presented. The Green-Kubo relations are used to calculate the heat transfer coefficients. The influence of dissipation (internal friction) on the heat transfer processes in nonideal systems is studied. New approximations are proposed for the thermal conductivity and diffusivity for nonideal dissipative systems. The results obtained are compared with the existing experimental and numerical data.
Khrustalyov, Yu. V., E-mail: firstname.lastname@example.com; Vaulina, O. S. [Russian Academy of Sciences, Joint Institute for High Temperatures (Russian Federation)
Numerical modeling and experimental measurements of water spray impact and transport over a cylinder.
This study compares experimental measurements and numerical simulations of liquid droplets over heated (to a near surface temperature of 423 K) and unheated cylinders. The numerical model is based on an unsteady Reynolds-averaged Navier-Stokes (RANS) formulation using a stochastic separated flow (SSF) approach for the droplets that includes submodels for droplet dispersion, heat and mass transfer, and impact on a solid surface. The details of the droplet impact model are presented and the model is used to simulate water spray impingement on a cylinder. Computational results are compared with experimental measurements using phase Doppler interferometry (PDI).
Avedisian, C. T. (Cornell University, Ithaca, NY); Presser, Cary (National Institute of Standard & Technology, Gaithersburg, MD); DesJardin, Paul Edward (University at Buffalo, New York, NY); Hewson, John C.; Yoon, Sam Sukgoo
Notes 13. Numerical methods for the dynamic response of MDOF damped systems
+DU+KU =F #0;#5;#0;#5; #0;#5; (1) where andU,U, U #0;#5;#0;#5;#0;#5; are the vectors of generalized displacement, velocity and acceleration, respectively; and ()t F is the vector of generalized (external forces) acting on the system. M,D,K represent... the matrices of inertia, viscous damping and stiffness coefficients, respectively 1 . The solution of Eq. (1) is unique with specified initial conditions , oo UU #0;#5; . Numerical solution methods ? direct numerical integration ? modal analysis (mode...
Massively-Parallel Direct Numerical Simulation of Gas Turbine Endwall Film-Cooling Conjugate Heat Transfer
MASSIVELY-PARALLEL DIRECT NUMERICAL SIMULATION OF GAS TURBINE ENDWALL FILM-COOLING CONJUGATE HEAT TRANSFER A Thesis by CHARLES MICHAEL MEADOR Submitted to the O ce of Graduate Studies of Texas A&M University in partial ful llment of the requirements... for the degree of MASTER OF SCIENCE December 2010 Major Subject: Mechanical Engineering MASSIVELY-PARALLEL DIRECT NUMERICAL SIMULATION OF GAS TURBINE ENDWALL FILM-COOLING CONJUGATE HEAT TRANSFER A Thesis by CHARLES MICHAEL MEADOR Submitted to the O ce of Graduate...
Meador, Charles Michael
The numerical solution of a nickel-cadmium battery cell model using the method of lines
THE NUMERICAL SOLUTION OF A NICKEL-CADMIUM BATTERY CELL MODEL USING THE METHOD OF LINES A Thesis by TESHOME HAILU Submitted to the Office of Graduate Studies Texas Adi:M University in partial fulfillment of the requirements for the degree... of MASTER OF SCIENCE December 1990 Major Subject: Chemical Engineering THE NUMERICAL SOLUTION OF A NICKEL-CADMIUM BATTERY CELL MODEL USING THE METHOD OF LINES A Thesis by TESHOME HAILU Approved as to style and content by: Ralph E. White (Chairman...
Hailu, Teshome
Ammon, A; Hartung, T; Jansen, K; Le�vey, H; Volmer, J
A study of some methods of numerical integration using floating point arithmetic
IN NUMERICAL INTEGRATION e ~ ~ e ~ ~ ~ ~ ~ e EVALUATION OF A FRESNEL INTEGRAL AS A POWER SERIES EXPANSION a ~ ~ o ~ ~ ~ ~ ~ DISCUSSIOM OF TWO MAJOR SOURCES OF ERROR WHEN INTEGRATING USING FLOATIMG POINT ARITHNEPICa e ~ ~ a ~ ~ e ~ ~ ~ ~ e ~ ~ ~ ~ 13 V...(x) ~ 10 for 0 c x c 10 f(x) = -10x + 110 for 10 c x c 11 ))0 c x c, 85 One of the Fresnel integrals was chosen because the Fresnel integrals are interesting to physicists as well as mathematicians. It was used only in the study of numerical...
Young, John Mackay
Arnold tongues for a resonant injection-locked frequency divider: analytical and numerical results
In this paper we consider a resonant injection-locked frequency divider which is of interest in electronics, and we investigate the frequency locking phenomenon when varying the amplitude and frequency of the injected signal. We study both analytically and numerically the structure of the Arnold tongues in the frequency-amplitude plane. In particular, we provide exact analytical formulae for the widths of the tongues, which correspond to the plateaux of the devil's staircase picture. The results account for numerical and experimental findings presented in the literature for special driving terms and, additionally, extend the analysis to a more general setting.
Bartuccelli, Michele V; Gentile, Guido; Schilder, Frank
Loop Current Mixed Layer Energy Response to Hurricane Lili (2002). Part II: Idealized Numerical Simulations
Loop Current Mixed Layer Energy Response to Hurricane Lili (2002). Part II: Idealized Numerical horizontal pressure gradient, wind energy transfer to the mixed layer can be more efficient in such a regime as compared to the case of an initially horizontally homogeneous ocean. However, nearly all energy is removed
Miami, University of
First numerical evidence of global Arnold diffusion in quasi--integrable systems
We provide numerical evidence of global diffusion occurring in slightly perturbed integrable Hamiltonian systems and symplectic maps. We show that even if a system is sufficiently close to be integrable, global diffusion occurs on a set with peculiar topology, the so--called Arnold web, and is qualitatively different from Chirikov diffusion, occurring in more perturbed systems.
Massimiliano Guzzo; Elena Lega; Claude Froeschle'
A numerical study of horizontal dispersion in a macro tidal basin
in Cobscook Bay, a macro tidal basin, is simulated using the three-dimensional, nonlinear, finite element regimes in different parts of Cobscook Bay. It is found that the effective Lagrangian dispersion ocean model, QUODDY_dry. Numerical particles are released from various transects in the bay at different
Maine, University of
American Solar Energy Society Proc. ASES Annual Conference, Raleigh, NC, EVALUATION OF NUMERICAL WEATHER PREDICTION
© American Solar Energy Society Proc. ASES Annual Conference, Raleigh, NC, EVALUATION;© American Solar Energy Society Proc. ASES Annual Conference, Raleigh, NC, irradiance forecasts over OF NUMERICAL WEATHER PREDICTION SOLAR IRRADIANCE FORECASTS IN THE US Richard Perez ASRC, Albany, NY, Perez
Perez, Richard R.
Analytical and Numerical Study of Joule Heating Effects on Electrokinetically Pumped Continuous Flow PCR Chips
Analytical and Numerical Study of Joule Heating Effects on Electrokinetically Pumped Continuous, and the potential for integration.1-3 Joule heating is inevitable when electrokinetic pumping is used Form: December 8, 2007 Joule heating is an inevitable phenomenon for microfluidic chips involving
Le Roy, Robert J.
Some links on this page may take you to non-federal websites. Their policies may differ from this site.
Website Policies and Important Links
Search capabilities provided by DOE/OSTI