Bennecke, W.M.
1996-10-01T23:59:59.000Z
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.
Poreda, Robert J.
. It remains in the public domain. 91 Geochemical Discrimination of Five Pleistocene Lava-Dam Outburst- Flood@usgs.gov) A B S T R A C T Pleistocene basaltic lava dams and outburst-flood deposits in the western Grand Canyon outburst-flood deposit came from a common source, a lava dam. With these data, it is possible
Mantle source provinces beneath the Northwestern USA delimited by helium isotopes in young basalts
Graham, David W.
phenocrysts in basalts from the eastern Snake River Plain (SRP), the Owyhee Plateau (OP) and the Oregon High: Received 17 April 2008 Accepted 11 December 2008 Available online 31 December 2008 Keywords: Snake River Plain High Lava Plains Yellowstone helium isotopes mantle plume We report new He, Nd and Sr isotope
Wells, S. [Univ. of California, Riverside, CA (United States). Dept. of Earth Sciences; McFadden, L.; Perry, F. [Univ. of New Mexico, Albuquerque, NM (United States). Dept. of Geology; Forman, S. [Ohio State Univ., Columbus, OH (United States). Byrd Polar Research Center; Crowe, B.; Pothis, J.; Olinger, C. [Los Alamos National Lab., NM (United States)
1992-12-31T23:59:59.000Z
Evaluation of volcanic hazards near the proposed high-level radioactive waste repository at Yucca Mountain provides the impetus for a series of detailed field and geochronologic studies of selected small late Quaternary basaltic scoria cones and lava flows in Nevada and California. Two of the most significant results of these studies are: the presence of chronostratigraphic units which indicate multiple eruptions with significant intervals of no activity between events (polycyclic volcanism); and a marked difference between conventional, numerical ages derived from K-Ar and Ar-40/Ar-39 methods and numerical, calibrated, and relative ages derived from thermoluminescence, cosmogenic He-3, the degree of soil development, and geomorphology of these volcanic landforms. Soil-bounded unconformities and buried stone pavements define the boundaries of chronostratigraphic units within these small volume basaltic centers. Volcanic centers displaying this type of stratigraphy may appear morphological simple but cannot be considered mongenetic. Recent studies by Perry and Crowe demonstrate that geochemical variations within a single basaltic volcanic center in NV are consistent with several magma batches forming a complex polycyclic volcano. The K-Ar and Ar-40/Ar-39 ages are 1--2 orders of magnitude older than either TL or cosmogenic He-3 and appear to have insufficient precision to constrain the ages of chronostratigraphic units within polycyclic volcanoes. In contrast, preliminary data indicate the TL and cosmogenic He-3 dating methods have the ability to resolve the late Quaternary volcanic stratigraphy, and results from these dating methods are consistent with the degree of soil development and geomorphic modification of the volcanic units. K-Ar and Ar-40/Ar-39 dates from these small basaltic volcanic centers have been used to calibrate new Quaternary dating methods, e.g. rock varnish, which in turn have been used to interpret landscape evolution in the SW US.
Decreasing Magmatic Footprints of Individual Volcanos in a Waning Basaltic Field
G.A> Valentine; F.V. Perry
2006-06-06T23:59:59.000Z
The distribution and characteristics of individual basaltic volcanoes in the waning Southwestern Nevada Volcanic Field provide insight into the changing physical nature of magmatism and the controls on volcano location. During Pliocene-Pleistocene times the volumes of individual volcanoes have decreased by more than one order of magnitude, as have fissure lengths and inferred lava effusion rates. Eruptions evolved from Hawaiian-style eruptions with extensive lavas to eruptions characterized by small pulses of lava and Strombolian to violent Strombolian mechanisms. These trends indicate progressively decreasing partial melting and length scales, or magmatic footprints, of mantle source zones for individual volcanoes. The location of each volcano is determined by the location of its magmatic footprint at depth, and only by shallow structural and topographic features that are within that footprint. The locations of future volcanoes in a waning system are less likely to be determined by large-scale topography or structures than were older, larger volume volcanoes.
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), 2002; Level:5(Million Cubic Feet) Oregon (Including Vehicle Fuel) (Million Cubic Feet) Natural Gas DeliveredStents -(NEI) ScienceNumerical simulations of current
G.A. Valentine; F.V. Perry; D. Krier; G.N. Keating; R.E. Kelley; A.H. Cogbill
2006-04-04T23:59:59.000Z
Five Pleistocene basaltic volcanoes in Crater Flat (southern Nevada) demonstrate the complexity of eruption processes associated with small-volume basalts and the effects of initial emplacement characteristics on post-eruptive geomorphic evolution of the volcanic surfaces. The volcanoes record eruptive processes in their pyroclastic facies ranging from ''classical'' Strombolian mechanisms to, potentially, violent Strombolian mechanisms. Cone growth was accompanied, and sometimes disrupted, by effusion of lavas from the bases of cones. Pyroclastic cones were built upon a gently southward-sloping surface and were prone to failure of their down-slope (southern) flanks. Early lavas flowed primarily southward and, at Red and Black Cone volcanoes, carried abundant rafts of cone material on the tops of the flows. These resulting early lava fields eventually built platforms such that later flows erupted from the eastern (at Red Cone) and northern (at Black Cone) bases of the cones. Three major surface features--scoria cones, lava fields with abundant rafts of pyroclastic material, and lava fields with little or no pyroclastic material--experienced different post-eruptive surficial processes. Contrary to previous interpretations, we argue that the Pleistocene Crater Flat volcanoes are monogenetic, each having formed in a single eruptive episode lasting months to a few years, and with all eruptive products having emanated from the area of the volcanoes main cones rather than from scattered vents. Geochemical variations within the volcanoes must be interpreted within a monogenetic framework, which implies preservation of magma source heterogeneities through ascent and eruption of the magmas.
Analytical-Numerical Modeling Of Komatiite Lava Emplacement And...
thicknesses over felsic tuffaceous substrates with various degrees of consolidation and water contents. Thermo-mechanical erosion becomes more effective for substrates that are...
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 Home Page on DeliciousPlasmaP a gHigh4-FD-a <Alaska UndergroundInformation An EmpiricalGeothermal-Energy| Open1986) |
C. Doc Richardson; Nancy W. Hinman; Lindsay J. McHenry; J. Michelle Kotler; Jill R. Scott
2012-05-01T23:59:59.000Z
Secondary deposits associated with the basaltic caves of Craters of the Moon National Monument (COM) in southern Idaho were examined using X-ray powder diffraction, X-ray fluorescence spectrometry, Fourier transform infrared spectrometry, and Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS). The secondary mineral assemblages are dominated by Na-sulfate minerals (thenardite, mirabilite) with a small fraction of the deposits containing minor concentrations of Na-carbonate minerals. The assemblages are found as white, efflorescent deposits in small cavities along the cave walls and ceilings and as localized mounds on the cave floors. Formation of the deposits is likely due to direct and indirect physiochemical leaching of meteoritic water through the overlying basalts. Whole rock data from the overlying basaltic flows are characterized by their extremely high iron concentrations, making them good analogs for martian basalts. Understanding the physiochemical pathways leading to secondary mineralization at COM is also important because lava tubes and basaltic caves are present on Mars. The ability of FTICR-MS to consistently and accurately identify mineral species within these heterogeneous mineral assemblages proves its validity as a valuable technique for the direct fingerprinting of mineral species by deductive reasoning or by comparison with reference spectra.
Karlsson, Brynjar
. It is known that basalt fibers have better tensile strength than E-glass fibers, greater failure strain than mention Compar model o fibres re time in regions failure o the mat then tak FRP com The use type is s function which m is anoth was des by dam has been of beam In colum column where om environ ned as conti red
Numerical Integration Numerical Summation
Cohen, Henri
Numerical Integration Numerical Summation Numerical Extrapolation Numerical Recipes for Multiprecision Computations #12;Numerical Integration Numerical Summation Numerical Extrapolation Multiprecision, integration, summation, extrapolation, evaluation of continued fractions, Euler products and sums, complete
Aciego, S.M.; Jourdan, F.; DePaolo, D.J.; Kennedy, B.M.; Renne, P.R.; Sims, K.W.W.
2009-10-01T23:59:59.000Z
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.
Pederson, Neil
Site name: Khorgo Lava Pine, Mongolia Collectors: A. Hessl, N. Pederson Empire, and modern Mongolia. Proceedings of the National Academy of Science collected in 2010 and 2012 across the Khorgo lava field near Tariat, Mongolia
Physical Constraints on Geologic CO2 Sequestration in Low-Volume Basalt Formations
Ryan M. Pollyea; Jerry P. Fairley; Robert K. Podgorney; Travis L. McLing
2014-03-01T23:59:59.000Z
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.
Surface temperature and spectral measurements at Santiaguito lava dome, Guatemala
Rose, William I.
Surface temperature and spectral measurements at Santiaguito lava dome, Guatemala Steve T. M, Michigan Technological University, Houghton, Michigan, USA Otoniel Matias INSIVUMEH, Guatemala City, Guatemala Received 4 June 2004; revised 23 July 2004; accepted 20 September 2004; published 13 October 2004
Cosmogenic 3 He production rates from Holocene lava flows
Licciardi, Joseph M.
rights reserved. Keywords: cosmogenic nuclides; helium; production rates; Icelandic Low; surface exposureCosmogenic 3 He production rates from Holocene lava flows in Iceland J.M. Licciardi a,, M.D. Kurz b Available online 25 April 2006 Editor: K. Farley Abstract We measured cosmogenic 3 He production rates
Distinctive upper mantle anisotropy beneath the High Lava Plains and Eastern Snake River Plain,
Distinctive upper mantle anisotropy beneath the High Lava Plains and Eastern Snake River Plain and continuing with the still- ongoing volcanism in the High Lava Plains (HLP) and eastern Snake River Plain (SRP waves; shear wave splitting; high lava plains; Snake River Plain; Yellowstone. Index Terms: 8137
Graham, David W.
Discussion Helium isotopes in historical lavas from Mount Vesuvius Comment on `Noble gas isotopic. Introduction Helium isotope results recently published by Tedesco et al. [1] appear to show a decrease in 3 He. Results Helium isotope results from our laboratory are reported in Table 1. The 3 He/4 He ratio has been
Repository site definition in basalt: Pasco Basin, Washington
Guzowski, R.V.; Nimick, F.B.; Muller, A.B.
1982-03-01T23:59:59.000Z
Discussion of the regional setting, geology, hydrology, and geochemistry of the Pasco Basin are included in this report. Pasco basin is a structural and topographic basin of approximately 2000 mi/sup 2/ (5180 km/sup 2/) located within the Yakima Fold Belt Subprovince of the Columbia Plateau. The stratigraphic sequence within the basin consists of an undetermined thickness of lower Miocene and younger flood basalts with interbedded and overlying sedimentary units. This sequence rests upon a basement of probably diverse rock types that may range in age from precambrian through early Tertiary. Although a large amount of information is available on the hydrology of the unconfined aquifer system, ground-water flow within the basin is, in general, poorly understood. Recharge areas for the Mabton interbed and the Saddle Mountains Formation are the highlands surrounding the basin with the flow for these units toward Gable Butte - Gable Mountain and Lake Wallula. Gable Butte - Gable Mountain probably is a ground-water sink, although the vertical flow direction in this zone is uncertain. The amount of upward vertical leakage from the Saddle Mountains Formation into the overlying sediments or to the Columbia River is unknown. Units underlying the Mabton interbed may have a flow scheme similar to those higher units or a flow scheme dominated by interbasin flow. Upward vertical leakage either throughout the basin, dominantly to the Columbia River, or dominantly to Lake Wallula has been proposed for the discharge of the lower units. None of these proposals is verified. The lateral and vertical distribution of major and minor ions in solution, Eh and pH, and ion exchange between basalt and ground-water are not well defined for the basin. Changes in the redox potential from the level of the subsurface facility to the higher stratigraphic levels along with the numerous other factors influencing K/sub d/, result in a poor understanding of the retardation process.
Elastic laboratory measurements and modeling of saturated basalts Ludmila Adam,1,2
Boise State University
the elastic behavior of basalt is important to seismically monitor volcanoes, subsea basalts, and carbon. Monitor- ing fluid substitution in basalts with seismic waves can be extended to volcanoes and subsea
Paces, J.B. (Michigan Technological Univ., Houghton (USA)); Bell, K. (Carleton Univ., Ottawa, Ontario (Canada))
1989-08-01T23:59:59.000Z
Midcontinent Rift flood basalts represent a sample of the relatively shallow, sub-continental upper mantle beneath the Canadian Shield at 1.1 Ga. A thick sequence of olivine tholeiite lavas, including minor intermediate to rhyolitic lavas, from the Portage Lake Volcanics (PLV) in northern Michigan have initial Nd and Sr isotopic compositions which cluster near Bulk Earth values. The effects of assimilation of old LREE-enriched continental crust into mantle-derived fractionating liquids are isotopically discernible in evolved lavas as well as in olivine tholeiites from the lowest portion of the volcanic pile. However, the effects of crustal contamination decrease with stratigraphic height and are absent in more primitive lavas in the upper half of the section. The source for PLV tholeiites is substantially less depleted than previously reported mantle values from the Superior Province. An origin for the PLV source is compatible with either of several mantle evolution models. The PLV source may have been associated with upwelling of a LIL element-enriched, asthenospheric plume which emplaced non-depleted material from deeper sources into the shallow sub-continental mantle beneath the Midcontinent Rift during continental break-up. Alternatively, the PLV source may have originated by enrichment of refractory sub-continental lithospheric mantle which was previously depleted in incompatible trace elements during Archean-aged melt extraction and continental crust formation. Concurrent generation of carbonatite magmas in other areas beneath the Superior Province indicates the widespread presence of sub-continental mantle with substantially higher {epsilon}{sub Nd}(T) and lower {epsilon}{sub Sr}(T) than the PLV source.
Consequences of melt transport for uranium series disequilibrium in young lavas
Spiegelman, Marc W.
Consequences of melt transport for uranium series disequilibrium in young lavas Marc Spiegelman do not actu- ally include melt transport. Here we explore the be- haviour of short
Peak discharge of a Pleistocene lava-dam outburst flood in Grand Canyon, Arizona, USA
in channel geometry, dam height, and hydraulic characteristics, ranged from 2.3 to 5.3 × 105 m3 s-1 Canyon; Colorado river; Pleistocene floods; Lava dams; Hydraulic modeling; Paleoflood indicators; DamPeak discharge of a Pleistocene lava-dam outburst flood in Grand Canyon, Arizona, USA Cassandra R
Microwave palaeointensities from a recent Mexican lava ow, baked sediments and reheated pottery
Utrecht, Universiteit
Microwave palaeointensities from a recent Mexican lava £ow, baked sediments and reheated pottery H form 10 June 2003; accepted 10 June 2003 Abstract Microwave palaeointensity (PI) estimates have been to overestimate the correct values. In addition to lava samples, microwave PIs were also determined from pottery
Cooling rate of some active lavas determined using an orbital imaging spectrometer
Wright, Robert
Click Here for Full Article Cooling rate of some active lavas determined using an orbital imaging flow is an important physical property to measure. Through its influence on lava crystallinity, cooling modeling problem that will aid in the analysis of data acquired by future hyperspectral remote sensing
THERMAL PROPERTIES OF GABLE MOUNTAIN BASALT CORES HANFORD NUCLEAR RESERVATION
Martinez-Baez, L.F.
2011-01-01T23:59:59.000Z
70 THERMAL PROPERTIES OF GABLE MOUNTAIN BASALT CORES HANFORDft); we used the data giyen for Gable Mountain K1005 for oursamples of Gable Mountain DB-5 (521 ft and 524 ft); and we
Origin of basalts by hybridisation in andesite-dominated arcs
Cassidy, Michael; Edmonds, Marie; Watt, Sebastian F. L.; Palmer, Martin R.; Gernon, Thomas M.
2015-02-16T23:59:59.000Z
basaltic magmas, with their higher bulk density over andesites from neighbouring volcanoes, ultimately may have been eruptible owing to both the transtensional tectonics imposed by offshore grabens (related to the oblique subduction of the Lesser Antilles...
Preliminary assessment of regional dispersivity of the Hanford basalts
LaVenue, Arthur Marsh
1985-01-01T23:59:59.000Z
PRELIMINARY ASSESSMENT OF REGIONAL DISPERSIVITY OF THE HANFORD BASALTS A Thesis by ARTHUR MARSH LAVENUE Submitted to the Graduate College of Texas A&M University in partial fulfillment of the requirements for the degree of MASTER OF SCIENCE... December 1985 Major Subject: Geophysics PRELIMINARY ASSESSMENT OF REGIONAL DISPERSIVITY OF THE HANFORD BASALTS A Thesis by ARTHUR MARSH LAVENUE Approved as to style and content by: ar . Hos ins (Co-Chair of Committee) atric . omenico (Co...
Shervais, John W.
High-K alkali basalts of the Western Snake River Plain: Abrupt transition from tholeiitic to mildly alkaline plume-derived basalts, Western Snake River Plain, Idaho John W. Shervais a, , Scott K. Vetter b Snake River Plain SRP Yellowstone plume hotspots Basaltic volcanism in the western Snake River Plain
Breakthroughs in Seismic and Borehole Characterization of Basalt Sequestration Targets
Sullivan, E. C.; Hardage, Bob A.; McGrail, B. Peter; Davis, Klarissa N.
2011-04-01T23:59:59.000Z
Mafic continental flood basalts form a globally important, but under-characterized CO2 sequestration target. The Columbia River Basalt Group (CRBG) in the northwestern U.S. is up to 5 km thick and covers over 168,000 km2. In India, flood basalts are 3 km thick and cover greater than 500,000 km2. Laboratory experiments demonstrate that the CRBG and other basalts react with formation water and super critical (sc) CO2 to precipitate carbonates, thus adding a potential mineral trapping mechanism to the standard trapping mechanisms of most other types of CO2 sequestration reservoirs. Brecciated tops of individual basalt flows in the CRBG form regional aquifers that locally have greater than 30% porosity and three Darcies of permeability. Porous flow tops are potential sites for sequestration of gigatons of scCO2 in areas where the basalts contain unpotable water and are at depths greater than 800 m. In this paper we report on the U.S. DOE Big Sky Regional Carbon Sequestration Partnership surface seismic and borehole geophysical characterization that supports a field test of capacity, integrity, and geochemical reactivity of CRBG reservoirs in eastern Washington, U.S.A. Traditional surface seismic methods have had little success in imaging basalt features in on-shore areas where the basalt is thinly covered by sediment. Processing of the experimental 6.5 km, 5 line 3C seismic swath included constructing an elastic wavefield model, identifying and separating seismic wave modes, and processing the swath as a single 2D line. Important findings include: (1) a wide variety of shear wave energy modes swamp the P-wave seismic records; (2) except at very short geophone offsets, ground roll overprints P-wave signal; and (3) because of extreme velocity contrasts, P-wave events are refracted at incidence angles greater than 7-15 degrees. Subsequent removal of S-wave and other noise during processing resulted in tremendous improvement in image quality. The application of wireline logging to onshore basalts is underexploited. Full waveform sonic logs and resistivity-based image logs acquired in the 1250 m basalt pilot borehole provide powerful tools for evaluating geomechanics and lithofacies. The azimuth of the fast shear wave is parallel to SH and records the changes through time in basalt flow and tectonic stress tensors. Combined with image log data, azimuthal S-wave data provide a borehole technique for assessing basalt emplacement and cooling history that is related to the development of reservoirs and seals, as well as the orientation of tectonic stresses and fracture systems that could affect CO2 transport or containment. Reservoir and seal properties are controlled by basalt lithofacies, and rescaled P- and S- wave slowness curves, integrated with image logs, provide a tool for improved recognition of subsurface lithofacies.
Preliminary Hydrogeologic Characterization Results from the Wallula Basalt Pilot Study
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
2009-12-01T23:59:59.000Z
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.
Origin of enriched ocean ridge basalts and implications for mantle dynamics
Langmuir, Charles H.
Origin of enriched ocean ridge basalts and implications for mantle dynamics Kathleen E. Donnellya that a small proportion of basalt are enriched in incompatible elements (enriched mid-ocean ridge basalts, E-MORB) from the MARK area. It is apparent that enriched magma sources, not associated with hot spots
Continental ood basalts: episodic magmatism above long-lived hotspots
Johnston, Stephen T.
November 1999 Abstract The eruption of continental flood basalt (CFB) may reflect episodic magmatism above long-lived mantle plumes. The Iceland and Yellowstone hotspots have generated successive CFB provinces in subducting oceanic lithosphere led to subsequent breakthrough and eruption of CFB. Since both mantle plume
TESTING MODELS FOR BASALTIC VOLCANISM: IMPLICATIONS FOR YUCCA MOUNTAIN, NEVADA
Conrad, Clint
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
Two new basaltic asteroids in the Outer Main Belt?
R. Duffard; F. Roig
2008-01-14T23:59:59.000Z
The identification of basaltic asteroids in the asteroid Main Belt and the description of their surface mineralogy is necessary to understand the diversity in the collection of basaltic meteorites. Basaltic asteroids can be identified from their visible reflectance spectra and are classified as V-type in the usual taxonomies. In this work, we report visible spectroscopic observations of two candidate V-type asteroids, (7472) Kumakiri and (10537) 1991 RY16, located in the outer Main Belt (a > 2.85 UA). These candidate have been previously identified by Roig and Gil-Hutton (2006, Icarus 183, 411) using the Sloan Digital Sky Survey colors. The spectroscopic observations have been obtained at the Calar Alto Observatory, Spain, during observational runs in November and December 2006. The spectra of these two asteroids show the steep slope shortwards of 0.70 microns and the deep absorption feature longwards of 0.75 microns that are characteristic of V-type asteroids. However, the presence of a shallow but conspicuous absorption band around 0.65 microns opens some questions about the actual mineralogy of these two asteroids. Such band has never been observed before in basaltic asteroids with the intensity we detected it. We discuss the possibility for this shallow absorption feature to be caused by the presence of chromium on the asteroid surface. Our results indicate that, together with (1459) Magnya, asteroids (7472) Kumakiri and (10537) 1991 RY16 may be the only traces of basaltic material found up to now in the outer Main Belt.
LAVA/CIS Version 2. 0: A software system for vulnerability and risk assessment
Smith, S.T.
1990-01-01T23:59:59.000Z
LAVA (the Los Alamos Vulnerability/Risk Assessment system) is an original systematic approach to risk assessment developed at the Los Alamos National Laboratory. It is an alternative to existing quantitative methods, providing an approach that is both objective and subjective, and producing results that are both quantitative and qualitative. LAVA was developed as a tool to help satisfy federal requirements for periodic vulnerability and risk assessments of a variety of systems and to satisfy the resulting need for an inexpensive, reusable, automated risk assessment tool firmly rooted in science. LAVA is a three-part systematic approach to risk assessment that can be used to model a variety of application systems such as computer security systems, communications security systems, information security systems, and others. The first part of LAVA is the mathematical model based on classical risk assessment, hierarchical multilevel system theory, decision theory, fuzzy possibility theory, expert system theory, utility theory, and cognitive science. The second part is the implementation of the mathematical risk model as a general software engine executed on a large class of personal computers. The third part is the application data sets written for a specific application system. The user of a LAVA application is not required to have knowledge of formal risk assessment techniques. All the technical expertise and specialized knowledge are built into the software engine and the application system itself. 36 refs., 5 figs.
The origin of young mare basalts inferred from lunar meteorites Northwest Africa 4734, 032, and La accepted 21 October 2013) AbstractNorthwest Africa (NWA) 4734 is an unbrecciated basaltic lunar meteorite basalts. NWA 4734 is a low-Ti mare basalt with a low Mg* (36.5) and elevated abundances of incompatible
Lava channel formation during the 2001 eruption on Mount Etna: evidence for mechanical erosion
Carmelo Ferlito; Jens Siewert
2005-12-14T23:59:59.000Z
We report the direct observation of a peculiar lava channel that was formed near the base of a parasitic cone during the 2001 eruption on Mount Etna. Erosive processes by flowing lava are commonly attributed to thermal erosion. However, field evidence strongly suggests that models of thermal erosion cannot explain the formation of this channel. Here, we put forward the idea that the essential erosion mechanism was abrasive wear. By applying a simple model from tribology we demonstrate that the available data agree favorably with our hypothesis. Consequently, we propose that erosional processes resembling the wear phenomena in glacial erosion are possible in a volcanic environment.
Singer, Bradley S.
Ar/Ar ages from transitionally magnetized lavas on La Palma, Canary Islands, and the geomagnetic the north and south walls of Barranco de los Tilos on the island of La Palma, Canary Islands, reveals from transitionally magnetized lavas on La Palma, Canary Islands, and the geomagnetic instability
The evolution of an active silicic lava flow field: an ETM+ perspective
Rose, William I.
'i, 2525 Correa Road, Honolulu, HI 96822, USA b INSIVUMEH, 7a Av. 14-57, Zona 13, Guatemala City, Guatemala, Houghton, MI 49931, USA d Santa Maria Volcano Observatory, c/o INSIVUMEH, 7a Av. 14-57, Zona 13, Guatemala City, Guatemala Accepted 5 December 2003 Abstract An active dacitic lava flow began advancing down
Short communication Origins and implications of zigzag rift patterns on lava lakes
Manga, Michael
the zigzags. Similar patterns are observed in analog wax experiments in which molten wax acts as cooling and solidifying lava. We perform a series of these wax experiments to find the relationship between , A of spreading and solidification speeds; the amplitude A is limited by the thickness of the solid wax crust
Smith, M.J.
1980-05-01T23:59:59.000Z
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.
Wright, Robert
Using infrared satellite data to drive a thermo-rheological/stochastic lava flow emplacement model the downhill paths that a lava flow could possibly follow and a thermorheological model to terminate, and A. J. L. Harris (2008), Using infrared satellite data to drive a thermo-rheological/stochastic lava
Suckale, Jenny
2011-01-01T23:59:59.000Z
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 ...
Final Reclamation Report: Basalt Waste Isolation Project exploratory shaft site
Brandt, C.A.; Rickard, W.H. Jr.
1990-06-01T23:59:59.000Z
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.
Evolution of Quaternary Tholeiitic Basalt Eruptive Centers on the Eastern Snake
Wetmore, Paul H.
Evolution of Quaternary Tholeiitic Basalt Eruptive Centers on the Eastern Snake River Plain, Idaho tholeiitic basalt eruptive centers on the eastern Snake River Plain, Idaho, in Bill Bonnichsen, C.M. White, and Michael McCurry, eds., Tectonic and Magmatic Evolution of the Snake River Plain Volcanic Province: Idaho
A lithospheric instability origin for Columbia River flood basalts and Wallowa Mountains uplift in
Humphreys, Eugene
from Yellowstone to the McDermitt caldera. Furthermore, the Snake River plain is a wellA lithospheric instability origin for Columbia River flood basalts and Wallowa Mountains uplift during the initiation of hotspot magmatism. The Columbia River basalts (CRB) represent the largest volume
Mason, Olivia U.; Di Meo-Savoie, Carol A.; Van Nostrand, Joy D.; Zhou, Jizhong; Fisk, Martin R.; Giovannoni, Stephen J.
2008-09-30T23:59:59.000Z
We used molecular techniques to analyze basalts of varying ages that were collected from the East Pacific Rise, 9 oN, from the rift axis of the Juan de Fuca Ridge, and from neighboring seamounts. Cluster analysis of 16S rDNA Terminal Restriction Fragment Polymorphism data revealed that basalt endoliths are distinct from seawater and that communities clustered, to some degree, based on the age of the host rock. This age-based clustering suggests that alteration processes may affect community structure. Cloning and sequencing of bacterial and archaeal 16S rRNA genes revealed twelve different phyla and sub-phyla associated with basalts. These include the Gemmatimonadetes, Nitrospirae, the candidate phylum SBR1093 in the c, andin the Archaea Marine Benthic Group B, none of which have been previously reported in basalts. We delineated novel ocean crust clades in the gamma-Proteobacteria, Planctomycetes, and Actinobacteria that are composed entirely of basalt associated microflora, and may represent basalt ecotypes. Finally, microarray analysis of functional genes in basalt revealed that genes coding for previously unreported processes such as carbon fixation, methane-oxidation, methanogenesis, and nitrogen fixation are present, suggesting that basalts harbor previously unrecognized metabolic diversity. These novel processes could exert a profound influence on ocean chemistry.
Thursday, March 26, 2009 POSTER SESSION II: SOURCES AND ERUPTION OF LUNAR BASALTS
Rathbun, Julie A.
Ar-39 Ar Ages and Chemical Composition for Lunar Mare Basalts: NWA 4734 and NWA 4898 [#1045] 40 Ar-39 Ar ages are reported for lunar mare basalts NWA 4734 and NWA 4898. The age obtained for NWA 4734 is 2
Vallino, Joseph J.
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
Petrology of some oceanic island basalts: PRIMELT2.XLS software for primary magma calculation
Petrology of some oceanic island basalts: PRIMELT2.XLS software for primary magma calculation C potential temperature; primary magma. Index Terms: 3621 Mineralogy and Petrology: Mantle processes (1038); 3610 Mineralogy and Petrology: Geochemical modeling (1009, 8410); 3630 Mineralogy and Petrology
Siebel, Wolfgang
Petrogenesis of the Linxi granitoids, northern Inner Mongolia of China: Constraints on basaltic in the Linxi area, northern Inner Mongolia of China, divided into two suites, are mainly granodiorite growth; Granitoid petrogenesis; Northern Inner Mongolia 1. Introduction Phanerozoic granitoids
Basalt petrogenesis beneath slow- and ultraslow-spreading Arctic mid-ocean ridges
Elkins, Lynne J
2009-01-01T23:59:59.000Z
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 ...
Jackson, Matthew G. (Matthew Gerald)
2008-01-01T23:59:59.000Z
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 ...
Nakata, Dorene Samantha
2010-01-01T23:59:59.000Z
Mid-ocean ridge basalts (MORBs) exhibit a wide range of CO2 concentrations, reflecting saturation to supersaturation (and rarely undersaturation) relative to their emplacement depths. In this study, we explore the mechanisms ...
Paris-Sud XI, Université de
-829 MOSSBAUER STUDIES OF ICELANDIC LAVAS 0. HELGASON, S. STEINTHORSSON Science Institute, University of Iceland ktudiks par la spectroscopie Mossbauer. On a trouvk qu'une haute vitesse de solidification cause un plus. Abstract. - Two samples of lava from Iceland have been studied by use of Mossbauer spec- troscopy
Reidel, Steve P.; Spane, Frank A.; Johnson, Vernon G.
2002-08-08T23:59:59.000Z
This report provides the technical background and a guide to characterizing a site for storing natural gas in the Columbia River Basalt
Hausrath, E.M.
2008-01-01T23:59:59.000Z
Basalt weathering rates on Earth and the duration of liquidCarlisle, PA 17013 Earth Sciences Division, Lawrencetime to liquid water: on Earth, mineral persistence times
Bindeman, Ilya N.
and compositional zoning in selected grains, and subtle to severe D18 O (melt-olivine) and D18 O (plagioclase isotope laser fluorination analyses of 55 individual and bulk olivine crystals, coexisting individual different lavas, and variable core-to-rim oxygen isotopic zoning is present in selected olivine grains. Many
O'Hara, Mhairi
2012-08-08T23:59:59.000Z
and Eyjfjallajökull in 2010, at the Fimmvörðuháls site in Iceland are used to validate the model. Results show that the model accurately replicates the paths taken by the lava flows during both events, while the calculated flow lengths from the Hekla 2000 eruption...
Demouchy, Sylvie
, northern Iceland P. Camps a, , B.S. Singer b , C. Carvallo c , A. Goguitchaichvili d , G. Fanjat a , B Brunhes reversal paleointensity ArAr dating Iceland Tjörnes We report paleomagnetic results from four overlapping stratigraphic sections (55 flows) through the lava pile in Tjörnes peninsula (North Iceland
Elliott, John
at Santiaguito volcano, Guatemala S.K. Ebmeier a,n , J. Biggs b , T.A. Mather a , J.R. Elliott a , G. Wadge c , F. We apply this to Santiaguito volcano, Guatemala, and measure increases in lava thickness of up to 140
Spane, F.A. Jr.; Webber, W.D.
1995-09-01T23:59:59.000Z
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.
Basalt waste isolation project. Quarterly report, October 1, 1980-December 31, 1980
Deju, R.A.
1981-02-01T23:59:59.000Z
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.
Lee, Cin-Ty Aeolus
Primitive off-rift basalts from Iceland and Jan Mayen: Os-isotopic evidence for a mantle source Volcanological Center, Institute of Earth Sciences, University of Iceland, Reykjavik, Iceland d Natural History elements, are presented for basalts from the three volcanic flank zones in Iceland and from Jan Mayen
Use of low frequencies for sub-basalt imaging Anton Ziolkowski,1
Edinburgh, University of
, 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
LETTER doi:10.1038/nature10749 Origin of Columbia River flood basalt controlled by
Liu, Lijun
Yellowstone hotspot track along the eastern Snake River plain (Fig. 1). Recent models trying to explainLETTER doi:10.1038/nature10749 Origin of Columbia River flood basalt controlled by propagating rupture of the Farallon slab Lijun Liu1 & Dave R. Stegman1 The origin of the SteensColumbia River (SCR
A shallow seismic reflection survey in basalts of the Snake River Plain, Idaho
Miller, Richard D.; Steeples, Don W.
1990-06-01T23:59:59.000Z
recorded along a 500 m long line over a basalt, rhyolite, and sedimentary sequence in the Snake River Plain. Some shallow reflections at 40 to 50 ms on the field files are of exceptional quality with frequency exceeding 150 Hz. Reflections and refractions...
Temperatures in ambient mantle and plumes: Constraints from basalts, picrites, and komatiites
below oceanic ridges, islands, and plateaus, on the basis of the petrology and geochemistry of erupted plumes; ambient mantle temperatures; petrology; mid ocean ridge basalt. Index Terms: 3610 Mineralogy and Petrology: Geochemical modeling (1009, 8410); 3619 Mineralogy and Petrology: Magma genesis and partial
Flow of Volatile-Free and Volatile-Rich Basaltic Magmas through Elastic-Walled Dikes
Al Hanbali, Ahmad
averaged in "one-dimensional" flow tube. #12;Questions: · Can we obtain a simple 1D symmetric model) · Anisotropy 1 in horizontal and vertical rock stress #12;Elastic flow equations · Continuity equation t(m bFlow of Volatile-Free and Volatile-Rich Basaltic Magmas through Elastic-Walled Dikes O. Bokhove
Generation of CO2-rich melts during basalt magma ascent and degassing
Boyer, Edmond
to solubilities. In contrast, the rate of vesiculation controls the final melt CO2 concentration. HighGeneration 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
The carbon dioxide solubility in alkali basalts: an experimental PRISCILLE LESNE 1,*
Paris-Sud XI, Université de
1 The carbon dioxide solubility in alkali basalts: an experimental study PRISCILLE LESNE 1 in both fluid and melt is required since, because of its low solubility, carbon dioxide is usually a major in silicate melts dramatically influence the physical properties of magmas, such as density, viscosity
Mineralization of Basalts in the CO2-H2O-H2S System. | EMSL
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
HT, BP McGrail, AT Owen, and BW Arey.2013."Mineralization of Basalts in the CO2-H2O-H2S System."International Journal of Greenhouse Gas Control 16:187-196. doi:10.1016...
Brandt, C.A.; Rickard, W.H. Jr.; Biehert, R.W.; Newell, R.L.; Page, T.L.
1989-01-01T23:59:59.000Z
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.
Zakharova, Natalia V.; Goldberg, David S.; Sullivan, E. C.; Herron, Michael M.; Grau, Jim A.
2012-11-02T23:59:59.000Z
Abstract This study presents borehole geophysical data and sidewall core chemistry from the Wallula Pilot Sequestration Project in the Columbia River flood basalt. The wireline logging data were reprocessed, core-calibrated and interpreted in the framework of reservoir and seal characterization for carbon dioxide storage. Particular attention is paid to the capabilities and limitations of borehole spectroscopy for chemical characterization of basalt. Neutron capture spectroscopy logging is shown to provide accurate concentrations for up to 8 major and minor elements but has limited sensitivity to natural alteration in fresh-water basaltic reservoirs. The Wallula borehole intersected 26 flows from 7 members of the Grande Ronde formation. The logging data demonstrate a cyclic pattern of sequential basalt flows with alternating porous flow tops (potential reservoirs) and massive flow interiors (potential caprock). The log-derived apparent porosity is extremely high in the flow tops (20%-45%), and considerably overestimates effective porosity obtained from hydraulic testing. The flow interiors are characterized by low apparent porosity (0-8%) but appear pervasively fractured in borehole images. Electrical resistivity images show diverse volcanic textures and provide an excellent tool for fracture analysis, but neither fracture density nor log-derived porosity uniquely correlate with hydraulic properties of the Grande Ronde formation. While porous flow tops in these deep flood basalts may offer reservoirs with high mineralization rates, long leakage migration paths, and thick sections of caprock for CO2 storage, a more extensive multi- well characterization would be necessary to assess lateral variations and establish sequestration capacity in this reservoir.
Jung, Hun Bok; Kabilan, Senthil; Carson, James P.; Kuprat, Andrew P.; Um, Wooyong; Martin, Paul F.; Dahl, Michael E.; Kafentzis, Tyler A.; Varga, Tamas; Stephens, Sean A.; Arey, Bruce W.; Carroll, KC; Bonneville, Alain; Fernandez, Carlos A.
2014-08-01T23:59:59.000Z
Composite Portland cement-basalt caprock cores with fractures, as well as neat Portland cement columns, were prepared to understand the geochemical and geomechanical effects on the integrity of wellbores with defects during geologic carbon sequestration. The samples were reacted with CO2-saturated groundwater at 50 ºC and 10 MPa for 3 months under static conditions, while one cement-basalt core was subjected to mechanical stress at 2.7 MPa before the CO2 reaction. Micro-XRD and SEM-EDS data collected along the cement-basalt interface after 3-month reaction with CO2-saturated groundwater indicate that carbonation of cement matrix was extensive with the precipitation of calcite, aragonite, and vaterite, whereas the alteration of basalt caprock was minor. X-ray microtomography (XMT) provided three-dimensional (3-D) visualization of the opening and interconnection of cement fractures due to mechanical stress. Computational fluid dynamics (CFD) modeling further revealed that this stress led to the increase in fluid flow and hence permeability. After the CO2-reaction, XMT images displayed that calcium carbonate precipitation occurred extensively within the fractures in the cement matrix, but only partially along the fracture located at the cement-basalt interface. The 3-D visualization and CFD modeling also showed that the precipitation of calcium carbonate within the cement fractures after the CO2-reaction resulted in the disconnection of cement fractures and permeability decrease. The permeability calculated based on CFD modeling was in agreement with the experimentally determined permeability. This study demonstrates that XMT imaging coupled with CFD modeling represent a powerful tool to visualize and quantify fracture evolution and permeability change in geologic materials and to predict their behavior during geologic carbon sequestration or hydraulic fracturing for shale gas production and enhanced geothermal systems.
Aciego, S.M.
2010-01-01T23:59:59.000Z
to Mineralogy and Petrology 84(4), 390-405. Frey F. A. ,to Mineralogy and Petrology 100, 383-397. Williams A.J. ,series lavas. Journal of Petrology 38(7), 911-939. Dunai T.
Washington at Seattle, University of
were probably triggered by destabilization of the chamber roof as magma was withdrawn during the Pagosa effusions of lava or high-energy plinian columns in which fragmentation of magma is efficient
Tian, Liyan
2011-01-01T23:59:59.000Z
and D. R. Hilton (2011), Petrology and Sr-Nd-Pb-He isotopeinvestigation of the petrology and geochemistry of basaltsF. Lonsdale (2008), Petrology and geochemistry of abandoned
Volcanism, the atmosphere and climate through time anja schmidt and alan robock
Robock, Alan
, 2003). Past episodes of continental flood basalt (CFB) volcanism produced huge lava volumes discussed in Chapters 5 and 11, CFB volcanism is typified by numerous, recurring large-volume eruptive an entire CFB province are short-lived typically about 1 Ma or less. For the assessment
Reidel, S.P.; Ledgerwood, R.K.; Myers, C.W.; Jones, M.G.; Landon, R.D.
1980-03-01T23:59:59.000Z
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.
Kelley, Jill Lynn
2002-01-01T23:59:59.000Z
reported values 44 5 Sample runs and duplicate runs 45 6 Major element analyses of glass, Leg 187 46 7 Major element and LOI analyses for fresh basalt samples, Leg 187 47 8 Major element analyses of altered basalts from Leg 187 49 9 Major oxides... and LOI analyses of altered dredge basalts 51 10 Computed mass change terms after Gresens and Grant 55 11 Correlation coefficients for Leg 187 altered basalts 59 ix LIST OF FIGURES FIGURE Page I Southeast Indian Ocean regional map showing...
Spane, Frank A.
2013-04-29T23:59:59.000Z
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
Edmonds, Marie
2014-12-18T23:59:59.000Z
/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...
Geometric Numerical Methods for Numerical Weather Prediction
Langdon, Stephen
-Mesh (HPM) Method · Label space is discretised into N particles with coordinates on the momentum phase space and Sij = (1 - ^2xx)-1. Geometric Numerical Methods for Numerical Weather Prediction p. 8/28 #12;HPM Equations of shallow water motions · The canonical HPM equations of 1D shallow water motion on TS1 are P
Dunbar, Nelia W.
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
Gent, Universiteit
Nuclear transparencies in relativistic A(e; e 0 p) models P. Lava a , M.C. Mart#19;#16;nez b , J Nuclear, Universidad de Sevilla, Apdo. 1065, E-41080 Sevilla, Spain c Departamento de F#19;#16;sica At#19;omica, Molecular y Nuclear, Facultad de Ciencias F#19;#16;sicas, Universidad Complutense de Madrid, E
Izett, G.A.; Pierce, K.L.; Naeser, N.D. (U.S.G.S., Denver, CO (United States)); Jaworowski, C. (Univ. of Wyoming, Laramie, WY (United States). Dept. of Geology and Geophysics)
1992-01-01T23:59:59.000Z
Along the Wind River near Kinnear (Pavillon quadrangle), a meter-thick tephra layer occurs near the middle of a main-stem gravel deposit about 100 m above the river. On Muddy Ridge 25 km east of Kinnear, a Lava Creek B tephra layer occurs at the base of a terrace deposit about 100 m above Muddy Creek. Another Lava Creek B tephra site 67 km northwest and upstream from Kinnear occurs within main-stem gravels of a terrace deposit 145 m above the river. This upstream increase of 45 m of the tephra horizon raises the concern that the two tephra layers might not be of the same age. All three tephras contain the same assemblage of phenocrysts as that in the Lava Creek Tuff, Member B in Yellowstone National Park and the Lava Creek B volcanic ash bed of the Western U.S., and therefore they are arguably correlatives. The authors confirmed this petrographic correlation by isotopic dating of sanidine crystals recovered from cm-size pumice lapilli in the Kinnear tephra and from coarse-grained tephra at the Muddy Creek site. Laser total-fusion Ar-40-Ar-39 ages of sanidine from the two sites are coeval, 0.66[plus minus]0.01 Ma and 0.67[plus minus]0.01 Ma at Muddy Creek. Conventional K-Ar dating of sanidine from the tephra at the Cl453 site resulted in an age of 0.60[plus minus]0.02 Ma. Glass-mantled zircon crystals from the Cl453 site yielded a fission-track age of 0.67[plus minus]0.16 Ma. These isotopic ages are compatible with conventional K-Ar, Ar-40-Ar-39, and fission-track ages of the Lava Creek Tuff, Member B in Yellowstone National Park and other occurrences of Lava Creek B ash beds. The authors suggest that the terrace deposit that contains the Lava Creek B tephra rises from the Kinnear site northwest up the Wind River as a result of Quaternary uplift in the area of the Yellowstone hotspot.
Seaman, S.; Dyar, D; Marinkovic, N
2009-01-01T23:59:59.000Z
This study focuses on the origin of flow-banded rhyolites that consist of compositionally similar darker and lighter flow bands of contrasting texture and color. Infrared radiation was used to obtain Fourier transform infrared (FTIR) spectra from which water concentrations were calculated, and to map variations in water concentrations across zones of spherulites and glass from the 23 million year old Sycamore Canyon lava flow of southern Arizona. Lighter-colored, thicker flow bands consist of gray glass, fine-grained quartz, and large (1.0 to 1.5 mm) spherulites. Darker-colored, thinner flow bands consist of orange glass and smaller (0.1 to 0.2 mm) spherulites. The centers of both large and small spherulites are occupied by either (1) a quartz or sanidine crystal, (2) a granophyric intergrowth, or (3) a vesicle. Mapping of water concentration (dominantly OH- in glass and OH- and H2O in sanidine crystals) illustrates fluctuating water availability during quenching of the host melt. Textures of large spherulites in the lighter (gray) bands in some cases indicate complex quenching histories that suggest that local water concentration controlled the generation of glass versus crystals. Small spherulites in darker (orange) bands have only one generation of radiating crystal growth. Both the glass surrounding spherulites, and the crystals in the spherulites contain more water in the gray flow bands than in the orange flow bands. Flow banding in the Sycamore Canyon lava flow may have originated by the stretching of a magma that contained pre-existing zones (vesicles or proto-vesicles) of contrasting water concentration, as the magma flowed in the conduit and on the surface. Variation in the original water concentration in the alternating layers is interpreted to have resulted in differences in undercooling textures in spherulites in the lighter compared to the darker flow bands.
Lovley, Derek
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
Stern, Robert J.
@jamstec.go.jp) Bradley R. Hacker Department of Earth Science, University of California, Santa Barbara, California 93106-9630, USA (hacker@geol.ucsb.edu) Peter E. van Keken Department of Geological Sciences, University. Hacker, P. E. van Keken, H. Kawabata, T. Yoshida, and R. J. Stern (2009), Arc Basalt Simulator version 2
temperature alteration. They represent products of very early planetesimal melting (4565.9 ± 0.3 Ma fundamentally basaltic in nature. These sodic plagioclase-rich rocks represent a series of early asteroidal high-temperature processes: (stage 1) melting and partial extraction of a low-temperature FeNiS melt, (stage 2) small
Simulating Geologic Co-sequestration of Carbon Dioxide and Hydrogen Sulfide in a Basalt Formation
Bacon, Diana H.; Ramanathan, Ramya; Schaef, Herbert T.; McGrail, B. Peter
2014-01-15T23:59:59.000Z
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.
Nicholson, S.W.; Cannon, W.F. (Geological Survey, Reston, VA (United States)); Woodruff, L.G. (Geological Survey, St. Paul, MN (United States))
1992-01-01T23:59:59.000Z
Whole-rock oxygen (n = 32) and strontium (n = 12) isotopic compositions of basalts erupted during the 1.1 Ga Midcontinent rift event provide evidence for both regional hydrothermal alteration during subsidence and more local alteration related to native Cu mineralization. In the Porcupine Mountains area of northern Michigan, the fissure-fed rift basalts of the Portage Lake volcanics show delta O-18 ca 8--9[per thousand] and Sr-87/Sr-86 = 0.7029--0.7042. Overlying the Portage Lake volcanics are the remnants of a stratavolcano, the Porcupine volcanics. Basalts near the base of the Porcupine Volcanics have delta O-18 ca 8--9[per thousand] and Sr-87/Sr-86 = 0.7027-0.7049, similar to those of the Portage Lake Volcanics. In contrast, basalts near the top of the Porcupine volcanics show higher delta O-18. These values are similar to values for thin basalt flows (Lake Shore Traps) intercalated with conglomerates and sandstones overlying the stratavolcano. The Sr isotopic analyses for basalts from the Porcupine Volcanics, Portage Lake Volcanics and Lake Shore Traps yield an age of 1,060.8 [+-] 66.2 Ma, most probably reflecting the age of the regional hydrothermal alteration. Elevated delta O-18 values of basalts from the lower part of the section can be attributed to pervasive regional prehnite-pumpellyite alteration during burial shortly after eruption. Native Cu mineralization is minor or absent in these rocks. Heavier O-isotopic compositions of the Lake Shore Traps and of the basalts at the top of the Porcupine volcanics reflect localized alteration near deposits of native copper; thus, these upper basalts may have reacted with hydrothermal fluids channeled through adjacent thick conglomerates.
2005-10-20T23:59:59.000Z
Page 1. LAB #8. Numerical Methods. Goal: The purpose of this lab is to explain how computers numerically ... Also you will examine what .... (7) Now consider the differential equation ... 3-exp(2*y)+sqrt(t)/y; (Don't forget the “;” at the end.).
Teague, L.S.
2010-01-01T23:59:59.000Z
of secondary and primary mineral phases in basalts of theB. , 1975. Perched silia minerals on mordenite fibers. Jour.lt;. UC—70 I SECONDARY MINERALS FOUND IN CORES DC2 AI AND
Santelli, Cara M
2007-01-01T23:59:59.000Z
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, ...
M. N. Pace; R. C. Bartholomay (USGS); J. J. Rosentreter (ISU)
1999-07-01T23:59:59.000Z
The U.S. Geological Survey and Idaho State University, in cooperation with the U.S. Department of Energy, are conducting a study to determine and evaluate strontium distribution coefficients (Kds) of subsurface materials at the Idaho National Engineering and Environmental Laboratory (INEEL). The purpose of this study is to aid in assessing the variability of strontium Kds at the INEEL as part of an ongoing investigation of chemical transport of strontium-90 in the Snake River Plain aquifer. Batch experimental techniques were used to determine Kds of six basalt core samples, five samples of sediment infill of vesicles and fractures, and six standard material samples. Analyses of data from these experiments indicate that the Kds of the sediment infill samples are significantly larger than those of the basalt samples. Quantification of such information is essential of furthering the understanding of transport processes of strontium-90 in the Snake River Plain aquifer and in similar environments.
Mineralization of Basalts in the CO2-H2O-SO2-O2 System
Schaef, Herbert T.; Horner, Jacob A.; Owen, Antionette T.; Thompson, Christopher J.; Loring, John S.; McGrail, B. Peter
2014-05-01T23:59:59.000Z
Sequestering carbon dioxide (CO2) containing minor amounts of co-contaminants in geologic formations was investigated in the laboratory through the use of high pressure static experiments. Five different basalt samples were immersed in water equilibrated with supercritical CO2 containing 1wt% sulfur dioxide (SO2) and 1wt% oxygen (O2) at reservoir conditions (~100 bar, 90°C) for 49 and 98 days. Gypsum (CaSO4) was a common precipitate, occurred early as elongated blades with striations, and served as substrates for other mineral products. Bimodal pulses of water released during dehydroxylation were key indicators along with X-ray diffraction for verifying the presences of jarosite-alunite group minerals. Well-developed pseudocubic jarosite crystals formed surface coatings, and in some instances mixtures of natrojarosite and natroalunite aggregated into spherically shaped structures measuring 100 ?m in diameter. Reaction products were also characterized using infrared spectroscopy, which indicated OH and Fe-O stretching modes. The presences of jarosite-alunite group minerals were found in the lower wavenumber region from 700–400 cm-1. A strong preferential incorporation of Fe(III) into natrojarosite was attributed to the oxidation potential of O2. Evidence of CO2 was detected during thermal decomposition of precipitates, suggesting the onset of mineral carbonation.
Site characterization report for the basalt waste isolation project. Volume II
None
1982-11-01T23:59:59.000Z
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.
Introduction to exact numerical computation
Escardó, Martín
Introduction to exact numerical computation Notes 2 Floating-point computation 4 3 Exact numerical (called the operational semantics) for computing a syntactical representative of the mathematical entity
Numerical Integration Gordon K. Smyth
Smyth, Gordon K.
Numerical Integration 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 Integration Numerical integration is the study of how the numerical value of an integral can be found. Also called
Aldrich, Jeffrey Brooks
1983-01-01T23:59:59.000Z
and tholeiitic glass 14 days 200'-500'C 500-800 bars smectite, anhydrite, albite, tremolite-actinolite, chalcopyrite, pyrrhotite, hematite. fluid became acidic rock ? + Mg, Na, N 0 Ca, K, Fe, Mn, Si, Ca Bischoff and Dickson (1975) basalt glass 4752... to matrix ratio is seldom less than 2:1 and frequently as high 27 as 5:1. The matrix composition varies considerably from nannofossil chalk in Unit 10 to extremely altered basalt glass s'nards in Unit 44. Smectite clays are cosssonly the primary...
11. NUMERICAL TECHNIQUES 1 Numerical identification of effective multipole
Boyer, Edmond
11. NUMERICAL TECHNIQUES 1 Numerical identification of effective multipole moments of polarizable of the induced multipole moments. A general multipole theory is available in the literature, however, only linear multipole model is usually exploited when determining numerically these effective moments. Since this axial
Busby, Cathy
the eastern Snake River Plain toward the Yellowstone caldera (Armstrong et al., 1975; Rodgers et al., 1990, and its age, volume, and chemistry indicate a genetic affinity with the Columbia River Basalt Group and suggest that it is coeval with the 16.115.0 Ma Imnaha and Grande Ronde flows of the Columbia River Basalt
Souza, Paulo Victor Santos
2015-01-01T23:59:59.000Z
Flow of viscous fluids are not usually discussed in detail in general and basic courses of physics. This is due in part to the fact that the Navier-Stokes equation has analytical solution only for a few restricted cases, while more sophisticated problems can only be solved by numerical methods. In this text, we present a computer simulation of wind tunnel, i.e., we present a set of programs to solve the Navier-Stokes equation for an arbitrary object inserted in a wind tunnel. The tunnel enables us to visualize the formation of vortices behind object, the so-called von K\\'arm\\'an vortices, and calculate the drag force on the object. We believe that this numerical wind tunnel can support the teacher and allow a more elaborate discussion of viscous flow. The potential of the tunnel is exemplified by the study of the drag on a simplified model of wing whose angle of attack can be controlled. A link to download the programs that make up the tunnel appears at the end.
Numerical Integration Gordon K. Smyth
Smyth, Gordon K.
Numerical Integration Gordon K. Smyth May 1997 Numerical integration is the study of how the numerical value of an integral can be found. Also called quadrature, which refers to finding a square whose \\Lambda . Of central interest is the process of approximating a definite integral from values of the in
Spera, Frank J.
-volume alkali basalt petrogenesis: implications for volcanic hazards at the proposed Yucca Mountain nuclear. Yucca Mountain (YM) in Nevada, USA was identified in the early 1980s as a potential geologic repository for nuclear waste. Yucca Mountain is made up of silicic volcanic tuffs- rocks composed chiefly of pyroclastic
Cohen, Barbara Anne
processing. Rapid extraction of partial melts from a range of depths could leave intact the range of characterized 171 feldspathic clasts in these meteorites in terms of texture, mineralogy, and mineral compo to represent a common igneous (generally basaltic) lithology and whose mafic minerals show a normal igneous
Kinetics of CO2-Fluid-Rock Reactions in a Basalt Aquifer, Soda Springs, Idaho
Maskell, Alexandra; Kampman, Niko; Chapman, Hazel; Condon, Daniel J.; Bickle, Mike
2015-06-24T23:59:59.000Z
the local geothermal gradient and spring temperatures. Values for ?G for plagioclase (?Gplagioclase) were determined using Eqn. 22, where a Keq value was calculated using SUPCRT92, using a plagioclase solid solution composition of Ca0.5Na0.5AlSi3O8... to the north (Armstrong et al., 1975; Hutsinpiller and Parry, 1985). 3. Hydrogeology and Flow Paths Numerous thermal CO2-rich springs, along with several artesian and abandoned geothermal exploration wells occur in the Blackfoot Volcanic Field...
Reidel, Steve P.; Spane, Frank A.; Johnson, Vernon G.
2005-09-24T23:59:59.000Z
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.
Suckale, Jenny
This is the second paper of two that examine numerical simulations of buoyancy-driven flow in the presence of large viscosity contrasts. In the first paper, we demonstrated that a combination of three numerical tools, an ...
The Numerical Simulation of Turbulence
W. Schmidt
2007-12-06T23:59:59.000Z
In this contribution, I give an overview of the various approaches toward the numerical modelling of turbulence, particularly, in the interstellar medium. The discussion is placed in a physical context, i. e. computational problems are motivated from basic physical considerations. Presenting selected examples for solutions to these problems, I introduce the basic ideas of the most commonly used numerical methods.
NUMERICAL ANALYSIS KENDALL E. ATKINSON
Atkinson, Kendall
of mathematics and computer science that creates, analyzes, and implements algorithms for solving nu- merically mathematical models in science and engineering, and numerical analysis of increasing sophistication has been of numerical analysis varies from quite theoretical mathematical studies (e.g. see [5]) to computer science
C. Doc Richardson; Nancy W. Hinman; Jill R. Scott
2013-10-01T23:59:59.000Z
Evidence of microbial activity associated with mineralization of secondary Na-sulphate minerals (thenardite, mirabilite) in the basaltic subsurface of Craters of the Moon National Monument (COM), Idaho were examined by scanning electron microscopy, X-ray diffraction, laser desorption Fourier transform ion cyclotron resonance mass spectrometry (LD-FTICR-MS), Fourier transform infrared spectroscopy (FTIR) and isotope ratio mass spectrometry. Peaks suggestive of bio/organic compounds were observed in the secondary Na-sulphate deposits by LD-FTICR-MS. FTIR provided additional evidence for the presence of bio/organic compounds. Sulphur fractionation was explored to assist in determining if microbes may play a role in oxidizing sulphur. The presence of bio/organic compounds associated with Na-sulphate deposits, along with the necessity of oxidizing reduced sulphur to sulphate, suggests that biological activity may be involved in the formation of these secondary minerals. The secondary Na-sulphate minerals probably form from the overlying basalt through leached sodium ions and sulphate ions produced by bio-oxidation of Fe-sulphide minerals. Since the COM basalts are one of the most comparable terrestrial analogues for their Martian counterparts, the occurrence of biological activity in the formation of sulphate minerals at COM has direct implications for the search for life on Mars. In addition, the presence of caves on Mars suggests the importance of these environments as possible locations for growth and preservation of microbial activity. Therefore, understanding the physiochemical pathways of abiotic and biotic mineralization in the COM subsurface and similar basaltic settings has direct implications for the search for extinct or extant life on Mars.
Johnson, Phillip R
2013-06-05T23:59:59.000Z
ELEMENTAL AND TECHNOLOGICAL ANALYSES OF BASALT ADZE MANUFACTURE ON TUTUILA, AMERIKA SAMOA: ECONOMIC INTENSIFICATION AND SPECIALIZATION DURING THE MONUMENT BUILDING PERIOD A Dissertation by PHILLIP RAY JOHNSON II Submitted... to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the requirements for the degree of DOCTOR OF PHILOSOPHY Chair of Committee, Suzanne Eckert Committee Members, William D. James Michael R. Waters Ted...
Crisman, D.P.; Jacobs, G.K.
1982-01-01T23:59:59.000Z
Copper deposits in the Portage Lake Volcanics of northern Michigan have been evaluated as a natural analogue for canister material to be emplaced in a nuclear waste repository located in the basalts beneath the Hanford Site near Richland, Washington. Te native copper, which precipitated from high temperature (200/sup 0/C to 300/sup 0/C) hydrothermal solutions between 500 and 800 million yr ago, has remained relatively unaltered in the Portage lake basalt-groundwater system. The results of this study illustrate the stability of copper in solutions of moderate pH, low to moderate Eh, and low total dissolved solids. Chemical trends and geochemical modeling of te near-surface waters suggestthat similarities (moderate pH, low E, low total dissolved solids) exist between the groundwater-basalt system of the Keweenaw Peninsula and the groundwater-basalt system at the Hanford Site. These similarities and the stability of copper in the Keweenswan basalts imply that copper alloys are adequate materials for nuclear waste canisters emplaced in a repository located in basalt.
Numerical Calculation of Granular Entropy
Asenjo, Daniel; Paillusson, Fabien; Frenkel, Daan
2014-03-05T23:59:59.000Z
We present numerical simulations that allow us to compute the number of ways in which N particles can pack into a given volume V. Our technique modifies the method of Xu, Frenkel, and Liu [Phys. Rev. Lett. 106, 245502 (2011)] and outperforms...
Planning numerical approximations Richard Power
Williams, Sandra
Planning numerical approximations Richard Power Sandra Williams 21st September 2009 #12;Table proportions (e.g., more than a quarter, 25.9 per cent) Proportions are a convenient well-defined subproblem Common in factual discourse (e.g., newspaper articles) Important for generating from data (but neglected
Numerical Simulation of Transpiration Cooling
to facilitate such numerical simulations for a carbon/carbon material mounted in the side wall of a hot gas channel that are able to capture a spatially varying interplay between the hot gas flow and the coolant itself. This calls for new cooling technologies, for example, in the combustion chamber of the rocket
Shear wave splitting and the pattern of mantle flow beneath eastern Oregon Maureen D. Long a,
period of bimodal (basaltic and silicic) volcanism in both the High Lava Plains and Snake River Plain October 2009 Editor: Y. Ricard Keywords: intraplate volcanism High Lava Plains Blue Mountains Pacific Oregon includes the volcanically active High Lava Plains (HLP) province and the accreted terrains
Numerical semigroups Easy bounds on ng
Elizalde, Sergi
Numerical semigroups Easy bounds on ng Improved bounds on ng The number of numerical semigroups #12;Numerical semigroups Easy bounds on ng Improved bounds on ng Motivation Definitions The tree on ng Improved bounds on ng Motivation Definitions The tree T of numerical semigroups Succession rules
Graduate Survey of Numerical Methods Background material
Corless, Robert M.
Numerical Methods, Numerical Analysis, Scientific Computing, Com- putational Mathematics, Computational. Numerical Meth- ods are what one uses to solve a problem from continuous mathematics (vide Nick Trefethen), and in what sense. Scientific Computing is the use of numerical methods to solve problems of scientific
Numerical bifurcation analysis of piecewise smooth systems
( ) ( ) ( ) ( ), , outoutinin inin xffxff xhhxgg xxxx == == #12;Numerical bifurcation analysis of piecewise smooth systems INRIA
Numerical simulation of hydraulic fracturing
Warner, Joseph Barnes
1987-01-01T23:59:59.000Z
~ared that the results of such treatments were not always adequately described by the two-dimensional models. With recent advances in hydraulic fracturing and computing technology, attempts have been made to formulate more realistic fracture models. These three...NUMERICAL SIMULATION OF HYDRAULIC FRACTURING A Thesis by JOSEPH BARNES WARNER 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 Maj or Subj ect...
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 Center Home Page on DeliciousPlasmaP aCentrothermDepew,IndependentCounty,KentPlatformLamoilleLauderdaleGeorgia: Energy ResourcesPage Edit
Disruptive Innovation in Numerical Hydrodynamics
Waltz, Jacob I. [Los Alamos National Laboratory
2012-09-06T23:59:59.000Z
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.
Numerical methods in stochastic processes
Corbett, Patrick James
1969-01-01T23:59:59.000Z
) for a large class of stochastic processes. The numerical solution of this problem was approached in the following manner. First of a1 1 equation (1. 1) was replaced by a di fference equation. Then the diff'erence ecuation wa. , solved for P(x, t...'s which would be used for i = n-1 would be those corresponding to r = 5. QP Using (2. 9) as an approximation for ? ', equation (2. 22) results 3x' in the difference equation (2. 10) where Z. is the expression on the right side of ecuation (2. 9...
Numerical Simulations of Bouncing Jets
Bonito, Andrea; Lee, Sanghyun
2015-01-01T23:59:59.000Z
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.
Numerical simulation of electrokinetically driven micro flows
Hahm, Jungyoon
2005-11-01T23:59:59.000Z
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...
Class Generation for Numerical Wind Atlases
Class Generation for Numerical Wind Atlases RisÃ¸ National Laboratory Wind Energy Department and The Technical University of Denmark Informatics and Mathematical Modelling Department Nicholas J. Cutler s000144 Constructing a Numerical Wind Atlas 5 2.1 Introduction
An implicit numerical algorithm general relativistic hydrodynamics
A. Hujeirat
2008-01-09T23:59:59.000Z
An implicit numerical algorithm general relativistic hydrodynamics This article has been replaced by arXiv:0801.1017
From Numerical Analysis to Computational Science
Li, Tiejun
. Numerical computing has, of course, been part of mathematics for a very long time. Al- gorithms by the namesFrom Numerical Analysis to Computational Science Bj¨orn Engquist · Gene Golub 1. Introduction The modern development of numerical computing is driven by the rapid in- crease in computer performance
Simulating Reionization in Numerical Cosmology
Aaron Sokasian; Tom Abel; Lars E. Hernquist
2001-05-10T23:59:59.000Z
The incorporation of radiative transfer effects into cosmological hydrodynamical simulations is essential for understanding how the intergalactic medium (IGM) makes the transition from a neutral medium to one that is almost fully ionized. Here, we present an approximate numerical method designed to study in a statistical sense how a cosmological density field is ionized by a set of discrete point sources. A diffuse background radiation field is also computed self-consistently in our procedure. The method requires relatively few time steps and can be employed with simulations having high resolution. We describe the details of the algorithm and provide a description of how the method can be applied to the output from a pre-existing cosmological simulation to study the systematic reionization of a particular ionic species. As a first application, we compute the reionization of He II by quasars in the redshift range 3 to 6.
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl. 2005; 12:683
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl. 2005; 12:683 Published Numerical Linear Algebra and its Applications The fourth workshop of the ERCIM Working Group on `Matrix Computations and Statistics' and the First International workshop on `Numerical Linear Algebra and its
Area E Numerics and Scientific Computing Foundation in Numerical Mathematics F4E1
Sturm, Karl-Theodor
Area E Numerics and Scientific Computing Foundation in Numerical Mathematics F4E1: · Scientific Selected Topics in Scientific Computing · V5E3 Advanced Topics in Numerical Methods in Science matrices) · Computational Finance (e.g. option pricing, fast numerical methods) · Visualization Methods (e
NUMERICAL ANALYSIS: This refers to the analysis of mathematical problems by numerical means, es-
Atkinson, Kendall
NUMERICAL ANALYSIS: This refers to the analysis of mathematical problems by numerical means, es- pecially mathematical problems arising from models based on calculus. Effective numerical analysis requires such as rootfinding and numerical integration; but we will also look at the structure of computers and the impli
MATHEMATICAL THEORY AND NUMERICAL METHODS FOR ...
2015-02-05T23:59:59.000Z
Weizhu Bao. Department of Mathematics and Center for Computational Science and Engineering ... Numerical methods for computing dynamics of GPE. 36. 4.1.
Adjoint Sensitivity Analysis for Numerical Weather Prediction
Alexandru Cioaca
2011-09-02T23:59:59.000Z
Sep 2, 2011 ... Adjoint Sensitivity Analysis for Numerical Weather Prediction: ... regional generation cost in the state of Illinois with respect to wind speed and ...
Lateral boundary errors in regional numerical weather
?umer, Slobodan
Lateral boundary errors in regional numerical weather prediction models Author: Ana Car Advisor weather services for short- range forecasts. These models are covering smaller areas with higher reso Introduction Equations for numerical weather prediction (NWP) are mathematical represen- ation of physical
Nonclassicality thresholds for multiqubit states: Numerical analysis
Gruca, Jacek; Zukowski, Marek [Institute of Theoretical Physics and Astrophysics, University of Gdansk, PL-80-952 Gdansk (Poland); Laskowski, Wieslaw [Institute of Theoretical Physics and Astrophysics, University of Gdansk, PL-80-952 Gdansk (Poland); Fakultaet fuer Physik, Ludwig-Maximilians Universitaet Muenchen, D-80799 Muenchen (Germany); Max-Planck-Institut fuer Quantenoptik, D-85748 Garching (Germany); Kiesel, Nikolai [Faculty of Physics, University of Vienna, A-1090 Vienna (Austria); Wieczorek, Witlef; Weinfurter, Harald [Fakultaet fuer Physik, Ludwig-Maximilians Universitaet Muenchen, D-80799 Muenchen (Germany); Max-Planck-Institut fuer Quantenoptik, D-85748 Garching (Germany); Schmid, Christian [Fakultaet fuer Physik, Ludwig-Maximilians Universitaet Muenchen, D-80799 Muenchen (Germany); Max-Planck-Institut fuer Quantenoptik, D-85748 Garching (Germany); European Organisation for Astronomical Research in the Southern Hemisphere, D-85748 Garching (Germany)
2010-07-15T23:59:59.000Z
States that strongly violate Bell's inequalities are required in many quantum-informational protocols as, for example, in cryptography, secret sharing, and the reduction of communication complexity. We investigate families of such states with a numerical method which allows us to reveal nonclassicality even without direct knowledge of Bell's inequalities for the given problem. An extensive set of numerical results is presented and discussed.
High performance computing and numerical Volker Springel
Masci, Frank
High performance computing and numerical modelling Volker Springel Lecture Notes 43rd Saas-mail: volker.springel@h-its.org 1 arXiv:1412.5187v1[astro-ph.GA]16Dec2014 #12;#12;Contents High performance computing and numerical modelling . . . . . . . . . . . . . . . 1 Volker Springel 1 Preamble
Numerical Study of a Turbulent Hydraulic Jump
Zhao, Qun
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
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl. 2011; 18:961980
De Sterck, Hans
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl. 2011; 18:961980 Published aggregation method of [1], and the Markov chain algebraic multigrid (MCAMG) method of [3] using the OTF
Buscaglia, Gustavo C.
2012-01-01T23:59:59.000Z
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING Int. J. Numer. Meth. Biomed iterative coupling of dimen- sionally heterogeneous flow models in computational hemodynamics. INTRODUCTION In recent years, several efforts have been directed at integrating different mathematical models
West, Matthew
momentum and energy conservation characteristics, as several numerical examples demonstrate. Copyright 2005INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2005; 64Center for Advanced Computing Research, California Institute of Technology, Mail Code 158-79, Pasadena
Stottlemyre, J.A.; Petrie, G.M.; Benson, G.L.; Zellmer, J.T.
1981-01-01T23:59:59.000Z
This report is a status report for an evolving methodology for release scenario development for underground nuclear waste repositories. As such, it is intended for use as a reference point and a preliminary description of an evolving geoscience methodology. When completed this methodology will be used as a tool in developing disruptive release scenarios for analyzing the long-term safety of geological nuclear waste repositories. While a basalt environment is used as an example, this report is not intended to reflect an actual site safety assessment for a repository in a media. It is rather intended to present a methodology system framework and to provide discussions of the geological phenomena and parameters that must be addressed in order to develop a methodology for potential release scenarios. It is also important to note that the phenomena, their interrelationships, and their relative importance along with the overall current structure of the model will change as new geological information is gathered through additional peer review, geotechnical input, site specific field work, and related research efforts.
Multiscale Problems: Numerical Analysis and Scientific Computing
Wirosoetisno, Djoko
/Transport Equations Alternative (Wind, Solar, Wave, . . . ) Weather and Climate Prediction Wave Energy Fuel CellsMultiscale Problems: Numerical Analysis and Scientific Computing with Applications in Energy inverse problems, data assimilation, stochastic differential equations Networks, compressed sensing
Toward standard testbeds for numerical relativity
Miguel Alcubierre; Gabrielle Allen; Carles Bona; David Fiske; Tom Goodale; F. Siddharta Guzman; Ian Hawke; Scott H. Hawley; Sascha Husa; Michael Koppitz; Christiane Lechner; Denis Pollney; David Rideout; Marcelo Salgado; Erik Schnetter; Edward Seidel; Hisa-aki Shinkai; Bela Szilagyi; Deirdre Shoemaker; Ryoji Takahashi; Jeffrey Winicour
2003-05-06T23:59:59.000Z
In recent years, many different numerical evolution schemes for Einstein's equations have been proposed to address stability and accuracy problems that have plagued the numerical relativity community for decades. Some of these approaches have been tested on different spacetimes, and conclusions have been drawn based on these tests. However, differences in results originate from many sources, including not only formulations of the equations, but also gauges, boundary conditions, numerical methods, and so on. We propose to build up a suite of standardized testbeds for comparing approaches to the numerical evolution of Einstein's equations that are designed to both probe their strengths and weaknesses and to separate out different effects, and their causes, seen in the results. We discuss general design principles of suitable testbeds, and we present an initial round of simple tests with periodic boundary conditions. This is a pivotal first step toward building a suite of testbeds to serve the numerical relativists and researchers from related fields who wish to assess the capabilities of numerical relativity codes. We present some examples of how these tests can be quite effective in revealing various limitations of different approaches, and illustrating their differences. The tests are presently limited to vacuum spacetimes, can be run on modest computational resources, and can be used with many different approaches used in the relativity community.
Volume 13, Number 1 21 August 2012
the Steens/Columbia River flood basalts (CRB) and the striking spatiotemporal trends of the Yellowstone/Snake River Plain (Y/SRP) and High Lava Plains (HLP) regions. Several different models have been proposed
UNITED STATES DEPARTMENT OF THE INTERIOR
Rathbun, Julie A.
UNITED STATES DEPARTMENT OF THE INTERIOR GEOLOGICAL SURVEY INTERAGENCY REPORT: ASTROGEOLOGY 26 · · · · · · · 0 · · · · · 4b. Skeletonized exploration of geologic units · o · · · · · · · · · · · · · . . . iii by a sequence of congealed basaltic lava flows, partly covered by surficial debris. The bedrock units exposed
Jones, Alan G.
by the extrusion of basaltic lava in Eocene and Miocene times that forms a variably thick sheet, averaging 100 m metres of gas and a billion cubic metres of oil (Hannigan et al., 1994). The thick volcanic cover limits
Provably Efficient Algorithms for Numerical Tensor Edgar Solomonik
California at Berkeley, University of
Provably Efficient Algorithms for Numerical Tensor Algebra Edgar Solomonik Electrical Engineering for Numerical Tensor Algebra by Edgar Solomonik A dissertation submitted in partial satisfaction-Gordon Fall 2014 #12;Provably Efficient Algorithms for Numerical Tensor Algebra Copyright 2014 by Edgar
High performance computing and numerical modelling
,
2014-01-01T23:59:59.000Z
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 ...
Numerical Simulation on Laser Fusion in China
Zhu Shaoping; Pei Wenbing; Xu Yan; Gu Peijun; Lan Ke; Ye Wenhua; Wu Junfeng; Li Jinghong; Gao Yaoming; Zheng Chunyang; Li Shuanggui; Mo Zeyao; Yan Jun [Institute of Applied Physics and Computational Mathematics, Beijing 100088 (China); Zhang Weiyan [National High-Tech Inertial Confinement Fusion Committee of China, Beijing 100088 (China)
2009-05-02T23:59:59.000Z
Numerical simulation is a powerful tool to get insight into the physics of laser fusion. Much effort has been devoted to develop the numerical simulation code series named LARED in China. The code series LARED are composed of six parts and enable us to have the simulation capability for the key processes in laser fusion. In recent years, a number of numerical simulations using LARED have been carried out and the simulation is checked by experiments done at the laser facility SG-II and SG-III prototype. In the present talk, some details of LARED code series will be introduced, and some simulation results, especially recent work on the opacities, will be shown.
KWIC Index for Numerical Linear Algebra
Carpenter, J.A.
1983-07-01T23:59:59.000Z
This report is a sequel to ORNL/CSD-106 in the ongoing supplements to Professor A.S. Householder's KWIC Index for Numerical Algebra. Beginning with the previous supplement, the subject has been restricted to Numerical Linear Algebra, roughly characterized by the American Mathematical Society's classification sections 15 and 65F but with little coverage of infinite matrices, matrices over fields of characteristics other than zero, operator theory, optimization and those parts of matrix theory primarily combinatorial in nature. Some consideration is given to the uses of graph theory in Numerical Linear Algebra, particularly with respect to algorithms for sparse matrix computations. The period covered by this report is roughly the calendar year 1982 as measured by the appearance of the articles in the American Mathematical Society's Contents of Mathematical Publications lagging actual appearance dates by up to nearly half a year. The review citations are limited to the Mathematical Reviews (MR).
Numerical simulation of an electroweak oscillon
Graham, N. [Department of Physics, Middlebury College, Middlebury, Vermont 05753 (United States)
2007-10-15T23:59:59.000Z
Numerical simulations of the bosonic sector of the SU(2)xU(1) electroweak standard model in 3+1 dimensions have demonstrated the existence of an oscillon--an extremely long-lived, localized, oscillatory solution to the equations of motion--when the Higgs mass is equal to twice the W{sup {+-}} boson mass. It contains total energy roughly 30 TeV localized in a region of radius 0.05 fm. A detailed description of these numerical results is presented.
Developments in Petroleum Science, 6 FUNDAMENTALS OF NUMERICAL
Santos, Juan
Developments in Petroleum Science, 6 FUNDAMENTALS OF NUMERICAL RESERVOIR SIMULATION DONALD WCongressCatalogingin PublicationData Peaceman, Donald W Fundamentals of numerical reservoir simulation. (develrpents in petroleum
Software for Numerical Methods for Partial Differential Equations
Software for Numerical Methods for Partial Differential Equations. This software was developed for and by the students in CS 615, Numerical Methods for Partial
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl. 2008; 15:115139
De Sterck, Hans
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl. 2008; 15:115139 Published and Ulrike Meier Yang2,, 1Department of Applied Mathematics, University of Waterloo, Waterloo, Ont., Canada N2L 3G1 2Center for Applied Scientific Computing, Lawrence Livermore National Laboratory, P.O. Box
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng (2014)
Leyendecker, Sigrid
2014-01-01T23:59:59.000Z
to ensure exact energy conservation in conservative systems, and the right rate of energy decayINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng (2014Computational Mechanics Group, California Institute of Technology, Mail Code 105-50, Pasadena, CA 91125, USA 2
Ortiz, Michael
in order to ensure exact energy conservation in conservative systems, and the right rate of energy decayINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2014; 100,*, and M. Ortiz1 1Computational Mechanics Group, California Institute of Technology, Mail Code 105
Jiao, Xiangmin "Jim"
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2004; 61 Jiao and Michael T. Heath, Computational Science and Engineering, University of Illinois at Urbana, and combustion, in an efficient and accurate manner. The multidisciplinary nature of such systems suggests
Kirby, Mike
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2008; 76 of Energy; contract/grant number: W-7405-ENG-48 Copyright q 2008 John Wiley & Sons, Ltd. #12;ANALYSIS AND REDUCTION OF QUADRATURE ERRORS IN MPM 923 by solving Newton's laws of motion for the internal force due
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng (2010)
Li, Shaofan
2010-01-01T23:59:59.000Z
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng (2010 Engineering, University of California, Berkeley, CA 94720, U.S.A. 2Faculty of Management and Administration is a subject of increasing importance to energy conversion/transfer, biochemistry, cellular and molecular
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2012; 91:126
Komatitsch, Dimitri
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2012; 91 Application of an elastoplastic spectral-element method to 3D slope stability analysis Hom Nath Gharti1-element method for 3D time-independent elastoplastic problems in geomechanics. As a first application, we use
Numerical Analysis Department Robustness improvement of
Vuik, Kees
Numerical Analysis Department Robustness improvement of polyhedral mesh method for airbag of polyhedral mesh method for airbag deployment simulations MASTER OF SCIENCE THESIS For the degree of Master entitled ROBUSTNESS IMPROVEMENT OF POLYHEDRAL MESH METHOD FOR AIRBAG DEPLOYMENT SIMULATIONS by SANTIAGO
Two physical characteristics of numerical apparent horizons
Ivan Booth
2007-12-20T23:59:59.000Z
This article translates some recent results on quasilocal horizons into the language of $(3+1)$ general relativity so as to make them more useful to numerical relativists. In particular quantities are described which characterize how quickly an apparent horizon is evolving and how close it is to either equilibrium or extremality.
Numerical Linear Algebra and Optimization on Facebook
Sidorov, Nikita
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
Numerical likelihood analysis of cosmic ray anisotropies
Carlos Hojvat et al.
2003-07-02T23:59:59.000Z
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.
A Numerical Method for Extended Boussinesq
Walkley, Mark
physi- cal, mathematical and computational considerations. A new formulation of internal wave generationA Numerical Method for Extended Boussinesq Shallow-Water Wave Equations by Mark Andrew Walkley School of Computer Studies September 1999 The candidate confirms that the work submitted is his own
Fluid = particles A Conservative Numerical Scheme
Al Hanbali, Ahmad
of this continuum formulation will 'directly' lead to a discrete Hamiltonian scheme: HPM. · alternative is fully Lagrangian description, which can be discretized into SPH. · both HPM and SPH require some sort of smoothing to maintain accuracy over long-time integrations. Eulerian grid in HPM allows for efficient numerical
Threedimensional numerical simulation for various geometries
Herbin, Raphaèle
modelling and numerical simulation of natural gasfed solid oxide cells (Solid Oxide Fuel Cell, SOFC) at a stationary regime. The principle of a Solid Oxide Fuel Cell (SOFC) is based on the conversion of the chemical is taken into account in the present model. The SOFC systems seem to be of great interest for use
Computer Algebra Applications for Numerical Relativity
Sascha Husa; Christiane Lechner
2003-01-21T23:59:59.000Z
We discuss the application of computer algebra to problems commonly arising in numerical relativity, such as the derivation of 3+1-splits, manipulation of evolution equations and automatic code generation. Particular emphasis is put on working with abstract index tensor quantities as much as possible.
MINIMIZING COMMUNICATION IN NUMERICAL LINEAR ALGEBRA*
California at Berkeley, University of
MINIMIZING COMMUNICATION IN NUMERICAL LINEAR ALGEBRA* GREY BALLARD , JAMES DEMMEL , OLGA HOLTZ, i.e., essentially all direct methods of linear al- gebra. The proof works for dense or sparse our lower bound technique to compositions of linear algebra operations (like computing powers
Numerical Solution of Ordinary Di erential Equations
... algorithms for ODEs and the mathematical analysis of their behaviour, cov- ..... (14). and then applying the numerical integration rule. Zxn+1. xn. g(x)dx hg(x n) ; ..... that Euler's method is a member of this family of methods, corresponding to c.
A numerical study of pseudoscalar inflation
Cheng, Shu-Lin; Ng, Kin-Wang
2015-01-01T23:59:59.000Z
A numerical study of a pseudoscalar inflation having an axion-photon-like coupling is performed by solving numerically the coupled differential equations of motion for inflaton and photon mode functions from the onset of inflation to the end of reheating. The backreaction due to particle production is also included self-consistently. We find that this particular inflation model realizes the idea of a warm inflation in which a steady thermal bath is established by the particle production. In most cases this thermal bath exceeds the amount of radiation released in the reheating process. In the strong coupling regime, the transition from the inflationary to the radiation-dominated phase does not involve either a preheating or reheating process. In addition, energy density peaks produced near the end of inflation may lead to the formation of primordial black holes.
Fast numerical test of hyperbolic chaos
Pavel V. Kuptsov
2011-11-21T23:59:59.000Z
The effective numerical method is developed performing the test of the hyperbolicity of chaotic dynamics. The method employs ideas of algorithms for covariant Lyapunov vectors but avoids their explicit computation. The outcome is a distribution of a characteristic value which is bounded within the unit interval and whose zero indicate the presence of tangency between expanding and contracting subspaces. To perform the test one needs to solve several copies of equations for infinitesimal perturbations whose amount is equal to the sum of numbers of positive and zero Lyapunov exponents. Since for high-dimensional system this amount is normally much less then the full phase space dimension, this method provide the fast and memory saving way for numerical hyperbolicity test of such systems.
Efficient Numerical Evaluation of Feynman Integral
Li, Zhao; Yan, Qi-Shu; Zhao, Xiaoran
2015-01-01T23:59:59.000Z
Feynman loop integral is the key ingredient of high order radiation effect, which is responsible for reliable and accurate theoretical prediction. We improve the efficiency of numerical integration in sector decomposition by implementing quasi-Monte Carlo method associated with the technique of CUDA/GPU. For demonstration we present the results of several Feynman integrals up to two loops in both Euclidean and physical kinematic regions in comparison with those obtained from FIESTA3. It is shown that both planar and non-planar two-loop master integrals in physical kinematic region can be evaluated in less than half minute with $\\mathcal{O}(10^{-3})$ accuracy, which makes the direct numerical approach viable for the precise investigation on the high order effect in multi-loop processes, e.g. the next-to-leading order QCD effect in Higgs pair production via gluon fusion with finite top quark mass.
Lattice Boltzmann Model for Numerical Relativity
Ilseven, E
2015-01-01T23:59:59.000Z
In the Bona-Masso formulation, Einstein equations are written as a set of flux conservative first order hyperbolic equations that resemble fluid dynamics equations. Based on this formulation, we construct a lattice Boltzmann model for Numerical Relativity. Our model is validated with well-established tests, showing good agreement with analytical solutions. Furthermore, we show that by increasing the relaxation time, we gain stability at the cost of losing accuracy, and by decreasing the lattice spacings while keeping a constant numerical diffusivity, the accuracy and stability of our simulations improves. Finally, in order to show the potential of our approach a linear scaling law for parallelisation with respect to number of CPU cores is demonstrated. Our model represents the first step in using lattice kinetic theory to solve gravitational problems.
A numerical study of pseudoscalar inflation
Shu-Lin Cheng; Wolung Lee; Kin-Wang Ng
2015-08-02T23:59:59.000Z
A numerical study of a pseudoscalar inflation having an axion-photon-like coupling is performed by solving numerically the coupled differential equations of motion for inflaton and photon mode functions from the onset of inflation to the end of reheating. The backreaction due to particle production is also included self-consistently. We find that this particular inflation model realizes the idea of a warm inflation in which a steady thermal bath is established by the particle production. In most cases this thermal bath exceeds the amount of radiation released in the reheating process. In the strong coupling regime, the transition from the inflationary to the radiation-dominated phase does not involve either a preheating or reheating process. In addition, energy density peaks produced near the end of inflation may lead to the formation of primordial black holes.
Origin and evolution of magmas on the Ontong Java Plateau J. GODFREY FITTON 1 & MARGUERITE GODARD 2
Demouchy, Sylvie
DrillingProject (DSDP) and Ocean DrillingProgram (ODP) drill sites on the plateau and in the adjacent on the plateau (Deep Sea Drilling Project (DSDP) Site 289 and Ocean Drilling Program (ODP) Sites 803 and 807) had penetrated basaltic basement. Five more basement sites were drilled during ODP Leg 192; basaltic lava flows
A numerical model of perturbation gas chromatography
DeBarro, Marc Joseph
1985-01-01T23:59:59.000Z
for in the model presented by Glover and Lsu. These mass transfer resistances are accommodated through an overall gss side mass transfer coefficient. Numerically generated peaks are compared to experimental profiles to con- firm the validity to the equations... and the method. Various solvent/polybutadiene systems are explored and individual profiles of the species are examined to under- stand the interaction which occurs between components for multicomponent sys- tems. Both ternary and quaternary systems...
Exploring New Physics Frontiers Through Numerical Relativity
Cardoso, Vitor; Herdeiro, Carlos; Sperhake, Ulrich
2014-01-01T23:59:59.000Z
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.
Beyond Numerics: The Existence of Pure Filtrations
Eisenbud, David; Schreyer, Frank-Olaf
2010-01-01T23:59:59.000Z
A recent result of Boij-Soederberg and Eisenbud-Schreyer proves that the Betti diagram of any graded module decomposes as a positive rational linear combination of pure diagrams. We consider the follow-up question of whether this numerical decomposition ever corresponds to an actual filtration of the minimal free resolution itself. Our main result is an affirmative answer to this question in many surprising cases. As applications of our technique, we also obtain new results about the semigroup of Betti diagrams and about very singular spaces of matrices.
Numerical modeling of vertical cavity semiconductor lasers
Chow, W.W.; Hadley, G.R.
1996-08-01T23:59:59.000Z
A vertical cavity surface emitting laser (VCSEL) is a diode laser whose optical cavity is formed by growing or depositing DBR mirror stacks that sandwich an active gain region. The resulting short cavity supports lasing into a single longitudinal mode normal to the wafer, making these devices ideal for a multitude of applications, ranging from high-speed communication to high-power sources (from 2D arrays). This report describes the development of a numerical VCSEL model, whose goal is to both further their understanding of these complex devices and provide a tool for accurate design and data analysis.
Numerical Simulations Unravel the Cosmic Web
C. -A. Faucher-Giguere; A. Lidz; L. Hernquist
2008-03-03T23:59:59.000Z
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.
RELAP-7 Numerical Stabilization: Entropy Viscosity Method
R. A. Berry; M. O. Delchini; J. Ragusa
2014-06-01T23:59:59.000Z
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.
Numerical Results for the Blue Phases
G. P. Alexander; J. M. Yeomans
2009-06-04T23:59:59.000Z
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.
GRChombo : Numerical Relativity with Adaptive Mesh Refinement
Clough, Katy; Finkel, Hal; Kunesch, Markus; Lim, Eugene A; Tunyasuvunakool, Saran
2015-01-01T23:59:59.000Z
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...
Numerical uncertainty in computational engineering and physics
Hemez, Francois M [Los Alamos National Laboratory
2009-01-01T23:59:59.000Z
Obtaining a solution that approximates ordinary or partial differential equations on a computational mesh or grid does not necessarily mean that the solution is accurate or even 'correct'. Unfortunately assessing the quality of discrete solutions by questioning the role played by spatial and temporal discretizations generally comes as a distant third to test-analysis comparison and model calibration. This publication is contributed to raise awareness of the fact that discrete solutions introduce numerical uncertainty. This uncertainty may, in some cases, overwhelm in complexity and magnitude other sources of uncertainty that include experimental variability, parametric uncertainty and modeling assumptions. The concepts of consistency, convergence and truncation error are overviewed to explain the articulation between the exact solution of continuous equations, the solution of modified equations and discrete solutions computed by a code. The current state-of-the-practice of code and solution verification activities is discussed. An example in the discipline of hydro-dynamics illustrates the significant effect that meshing can have on the quality of code predictions. A simple method is proposed to derive bounds of solution uncertainty in cases where the exact solution of the continuous equations, or its modified equations, is unknown. It is argued that numerical uncertainty originating from mesh discretization should always be quantified and accounted for in the overall uncertainty 'budget' that supports decision-making for applications in computational physics and engineering.
Testing and improving the numerical accuracy of the NLO predictions
R. Pittau
2010-09-24T23:59:59.000Z
I present a new and reliable method to test the numerical accuracy of NLO calculations based on modern OPP/Generalized Unitarity techniques. A convenient solution to rescue most of the detected numerically inaccurate points is also proposed.
Numerical simulation of micro-fluidic passive and active mixers
Kumar, Saurabh
2002-01-01T23:59:59.000Z
Numerical simulations of mixing using passive and active techniques are performed. For passive mixing, numerical modeling of a micro-fluidic device, build by Holden and Cremer, was performed. The micro-fluidic device consists of a Y...
NEW NUMERICAL TECHNOLOGIES FOR THE SIMULATION OF ARC WELDING PROCESSES
Paris-Sud XI, Université de
NEW NUMERICAL TECHNOLOGIES FOR THE SIMULATION OF ARC WELDING PROCESSES Michel Bellet 1 , Makhlouf Antipolis, France; soudage@transvalor.com Keywords: welding, finite elements, material deposit, adaptive for arc welding simulation and analysis. The new numerical technologies essentially consist first
Numerical Simulation of Laminar Reacting Flows with Complex Chemistry
Bell, John B.
Numerical Simulation of Laminar Reacting Flows with Complex Chemistry M S Day and J B Bell Lawrence: 47.40.Fw, 82.40.Py Submitted to: Combust. Theory Modelling #12;Numerical Simulation of Laminar
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 Home Page on DeliciousPlasmaP a gHigh Plains Wind Farm JumpKahului, Hawaii: EnergyKonicaLakefront
Numerical simulations for nodal domains and spectral minimal partitions
Vial, Grégory
unpublished results of [HHO2] with efficient numerical computations. This is the main goal of this paper
Symbolic Numeric Index Analysis Algorithm for Differential Algebraic Equations
Linninger, Andreas A.
. For the numerical solution of BVPs, global methods such as orthogonal collocation,1 finite elements,2 or spectral
Generation and Optimization of Numerical Programs by Symbolic Mathematical Methods
Fateman, Richard J.
Generation and Optimization of Numerical Programs by Symbolic Mathematical Methods Richard J FOR NUMERICAL COMPUTING, May 17-19, 1993 1 #12;5 The Future 21 5.1 Symbolic tools available in some form or numeric . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3 1.2 Recent CAS
Using the IMSL MATH/LIBRARY in Numerical Methods Courses
Toronto, University of
users, not to clone more numerical analysts or mathematical software developers. Some knowledge computing practices. In previous courses, we used software packages distributed with the pop ular numericalUsing the IMSL MATH/LIBRARY in Numerical Methods Courses by K. R. Jackson T. E. Hull University
Edmund G. Brown, Jr. NUMERICAL INVESTIGATION OF AIR-
Edmund G. Brown, Jr. Governor NUMERICAL INVESTIGATION OF AIR- COOLED STEAM CONDENSER PERFORMANCE.T.F. and D.G. Kröger. 2011. Numerical Investigation of AirCooled Steam Condenser Performance Under Windy · Transportation Numerical Investigation of AirCooled Steam Condenser Performance Under Windy Conditions
Upcoming Numerical Linear Algebra Issues in Systems and Control Theory
Upcoming Numerical Linear Algebra Issues in Systems and Control Theory Paul Van Dooren University in the interdisciplinary area of numerical linear algebra and systems and control theory. Although we do not claimto, and normalized coprime factorizations in robust control. Key words. Numerical algorithms, linear algebra, sparse
Upcoming Numerical Linear Algebra Issues in Systems and Control Theory
Upcoming Numerical Linear Algebra Issues in Systems and Control Theory Paul Van Dooren University in the interdisciplinary area of numerical linear algebra and systems and control theory. Although we do not claim, and normalized coprime factorizations in robust control. Key words. Numerical algorithms, linear algebra, sparse
Paris-Sud XI, Université de
INTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS Int. J. Numer. Anal Methods in Geomechanics (2012) 19 pages" DOI : 10.1002/nag.1138 #12;2 B. NEDJAR AND R. LE ROY asymptotic
Peirce, Anthony
INTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS Int. J. Numer. Anal. In the petroleum industry, hydraulic fractures are deliberately created in oil and gas reservoirs to substantially
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2013; 71:175184
Navon, Michael
laser heat deposition on the flow field near shock crossing. The influence of energy impulse was imi.1002/fld.3652 Numerical control of two-dimensional shock waves in dual solution domain by instant reflections in the dual solution domain is addressed in an optimization statement. The gradient
High numerical aperture multilayer Laue lenses
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Morgan, Andrew J.; Prasciolu, Mauro; Andrejczuk, Andrzej; Krzywinski, Jacek; Meents, Alke; Pennicard, David; Graafsma, Heinz; Barty, Anton; Bean, Richard J.; Barthelmess, Miriam; et al
2015-06-01T23:59:59.000Z
The ever-increasing brightness of synchrotron radiation sources demands improved X-ray optics to utilise their capability for imaging and probing biological cells, nanodevices, and functional matter on the nanometer scale with chemical sensitivity. Here we demonstrate focusing a hard X-ray beam to an 8 nm focus using a volume zone plate (also referred to as a wedged multilayer Laue lens). This lens was constructed using a new deposition technique that enabled the independent control of the angle and thickness of diffracting layers to microradian and nanometer precision, respectively. This ensured that the Bragg condition is satisfied at each point along themore »lens, leading to a high numerical aperture that is limited only by its extent. We developed a phase-shifting interferometric method based on ptychography to characterise the lens focus. The precision of the fabrication and characterisation demonstrated here provides the path to efficient X-ray optics for imaging at 1 nm resolution.« less
Very high numerical aperture light transmitting device
Allison, Stephen W. (Knoxville, TN); Boatner, Lynn A. (Oak Ridge, TN); Sales, Brian C. (Knoxville, TN)
1998-01-01T23:59:59.000Z
A new light-transmitting device using a SCIN glass core and a novel calcium sodium cladding has been developed. The very high index of refraction, radiation hardness, similar solubility for rare earths and similar melt and viscosity characteristics of core and cladding materials makes them attractive for several applications such as high-numerical-aperture optical fibers and specialty lenses. Optical fibers up to 60 m in length have been drawn, and several simple lenses have been designed, ground, and polished. Preliminary results on the ability to directly cast optical components of lead-indium phosphate glass are also discussed as well as the suitability of these glasses as a host medium for rare-earth ion lasers and amplifiers.
A numerical approximation to distribution function
Tuttle, Keith Allan
1977-01-01T23:59:59.000Z
then is to approximate F(Y) numerically. We sub- divide the unit cube Q into N cubes eqch of length h I/N on a n n side. Within each individual subcube Q , we will define an affine approximation to f, W = g(x , . . . , x ), which requires the gradient of 1' ' 'n our... & h/2, k k i + 1/2 2 n k k = 1, 2, . . . , n]. Let f(x) f C [Q ] the space of twice continuously n differentiable functions on Q , and define the auxiliary function g(x) as the tangent to f at x. i + 1 2 for x f Q. . That is, if x I Q i (x) f (x1 2...
NUMERICAL MODELING OF CATHODE CONTACT MATERIAL DENSIFICATION
Koeppel, Brian J.; Liu, Wenning N.; Stephens, Elizabeth V.; Khaleel, Mohammad A.
2011-11-01T23:59:59.000Z
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.
Dynamics of Dry Friction: A Numerical Investigation
Y. F. Lim; Kan Chen
1998-03-18T23:59:59.000Z
We perform extended numerical simulation of the dynamics of dry friction, based on a model derived from the phenomenological description proposed by T. Baumberger et al.. In the case of small deviation from the steady sliding motion, the model is shown to be equivalent to the state- and rate-dependent friction law which was first introduced by Rice and Ruina on the basis of experiments on rocks. We obtain the dynamical phase diagram that agrees well with the experimental results on the paper-on-paper systems. In particular, the bifurcation between stick-slip and steady sliding are shown to change from a direct (supercritical) Hopf type to an inverted (subcritical) one as the driving velocity increases, in agreement with the experiments.
Numerical Analysis of Black Hole Evaporation
Tsvi Piran; Andrew Strominger
1993-04-28T23:59:59.000Z
Black hole formation/evaporation in two-dimensional dilaton gravity can be described, in the limit where the number $N$ of matter fields becomes large, by a set of second-order partial differential equations. In this paper we solve these equations numerically. It is shown that, contrary to some previous suggestions, black holes evaporate completely a finite time after formation. A boundary condition is required to evolve the system beyond the naked singularity at the evaporation endpoint. It is argued that this may be naturally chosen so as to restore the system to the vacuum. The analysis also applies to the low-energy scattering of $S$-wave fermions by four-dimensional extremal, magnetic, dilatonic black holes.
Direct numerical simulation of turbulent reacting flows
Chen, J.H. [Sandia National Laboratories, Livermore, CA (United States)
1993-12-01T23:59:59.000Z
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.
Direct numerical simulations of convective heat transfer
Pointel, G.; Acharya, S.; Sharma, C. [Louisiana State Univ., Baton Rouge, LA (United States). Mechanical Engineering Dept.
1996-11-01T23:59:59.000Z
This paper deals with the development of a direct numerical simulation (DNS) code for solving the incompressible Navier-Stokes equation using higher order finite difference schemes. The time dependent Navier Stokes equation has been discretized using semi-implicit second order time splitting scheme, which requires the solution of pressure Poisson equation. For this purpose a Galerkin Fourier transform in the spanwise direction and a matrix diagonalization technique is used. The convection terms are formulated in non-conservative form on a collocated grid. A fifth order upwind biased scheme is used for this purpose. Diffusion terms are differenced using a sixth order central difference scheme. The algorithm is implemented on the MasPar MP-1, a Single Instruction Multiple Data computer where efficient data parallelization is used to get DNS results. The code has been used to get results for smooth channel flow at Re{sub {tau}} = 180. Results are now being obtained for the energy equation and for flow in a periodic ribbed channel.
Numerical treatment of disorder in PHC slabs
Minkov, Momchil
2015-01-01T23:59:59.000Z
This work concentrates on numerical simulations of Photonic Crystal structures using basis-expansion methods, with a main focus on simulating disorder. The plane-wave and guided-mode expansions are outlined as tools to compute the Bloch modes of a structure, on the basis of which the Bloch-mode expansion formalism is outlined - the latter allowing for simulations of large structures in presence of disorder. As a first illustration of the method, we apply it to three gentle-confinement cavities, to obtain results for their quality factors similar to the theoretically predicted in the literature. Furthermore, we compute that random disorder can drive those factors down to the experimentally measured values. As a second application, we study the effect of irregular hole shapes in a PHC waveguide, and find that the correlation length of the irregularity (i.e. the typical scale of the roughness of the features) matters: for higher correlation lengths, the computed modes show both higher band broadening and higher ...
Force budget: I. Theory and numerical methods
van der Veen, Cornelis J.; Whillans, I. M.
1989-10-05T23:59:59.000Z
,y) and vertical (z). Balance of forces is expressed by: 8 h- J Rxxdz8x b ~ 8 h ~ Rxx(h)- + - J Rxydz - Rxy(h)-8 +8x 8y yb (11 ) 8C1XX 8C1xy 8C1xz + + 0 (1)8x 8y 8z 8C1xy 8C1yy 8C1yz 0 (2)+ +8x 8y 8z 8 h 8x J Rxydz b (12) 8C1XZ 8C1yz 8C1ZZ + + ---- - pg = 0 (3)8x 8... are small compared to pg. This is satisfied if 8Rxz/ax = aC1xz/aX * 0 and 8Ryz/ay = aC1yz/ay * O.When the balance equations are solved numerically, this approximation is not necessary, and the value of Rzz can be obtained from Equation (3) for vertical force...
Direct numerical simulations of aeolian sand ripples
Orencio Duran; Philippe Claudin; Bruno Andreotti
2014-11-07T23:59:59.000Z
Aeolian sand beds exhibit regular patterns of ripples resulting from the interaction between topography and sediment transport. Their characteristics have been so far related to reptation transport caused by the impacts on the ground of grains entrained by the wind into saltation. By means of direct numerical simulations of grains interacting with a wind flow, we show that the instability turns out to be driven by resonant grain trajectories, whose length is close to a ripple wavelength and whose splash leads to a mass displacement towards the ripple crests. The pattern selection results from a compromise between this destabilizing mechanism and a diffusive downslope transport which stabilizes small wavelengths. The initial wavelength is set by the ratio of the sediment flux and the erosion/deposition rate, a ratio which increases linearly with the wind velocity. We show that this scaling law, in agreement with experiments, originates from an interfacial layer separating the saltation zone from the static sand bed, where momentum transfers are dominated by mid-air collisions. Finally, we provide quantitative support for the use the propagation of these ripples as a proxy for remote measurements of sediment transport.
A numerical model of aerosol scavenging
Bradley, M.M.; Molenkamp, C.R.
1991-10-01T23:59:59.000Z
Using a three-dimensional numerical cloud/smoke-plume model, we have simulated the burning of a large, mid-latitude city following a nuclear exchange. The model includes 18 dynamic and microphysical equations that predict the fire-driven airflow, cloud processes, and smoke-cloud interactions. In the simulation, the intense heating from the burning city produces a firestorm with updraft velocities exceeding 60 m/s. Within 15 minutes of ignition, the smoke plume penetrates the tropopause. The updraft triggers a cumulonimbus cloud that produces significant quantities of ice, snow, and hail. These solid hydrometeors, as well as cloud droplets and rain, interact with the smoke particles from the fire. At the end of the one-hour simulation, over 20% of the smoke is in slowly falling snowflakes. If the snow reaches the ground before the flakes completely sublimate (or melt and then evaporate), then only approximately 50% of the smoke will survive the scavenging processes and remain in the atmosphere to affect the global climate.
Collisionless microinstabilities in stellarators. II. Numerical simulations
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)
2013-12-15T23:59:59.000Z
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.
Numerical simulation of three-dimensional electrical flow through geomaterials
Akhtar, Anwar Saeed
1998-01-01T23:59:59.000Z
95 99 V ELECTRICAL FLOW AROUND AN ELECTRICAL CONE PENETROMETER 104 5. 1 INTRODUCTION 5. 2 ANALYTICAL SOLUTION FOR ELECTRICAL FLOW AROUND AN ELECTRICAL CONE PENETROMETER 5. 3 NUMERICAL INVESTIGATION 5. 4 COMPARISON OF ANALYTICAL AND NUMERICAL... RESULTS 5. 5 CONCLUSION AND APPLICATION 5. 5. 1 Utilization of Numerical Results 104 106 110 113 115 116 VI EXPERIMENTAL EQUIPMENT DESIGN 121 6. 1 INTRODUCTION 6. 2 ELECTRICAL POWER SOURCE 6. 3 ELECTRICAL RESISTIVITY CONE PENETROMETER 6. 4...
Status of NINJA: the Numerical INJection Analysis project
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
2009-05-26T23:59:59.000Z
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.
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 ...
Optimization Online - On mutual impact of numerical linear algebra ...
Marco D'Apuzzo
2008-03-12T23:59:59.000Z
Mar 12, 2008 ... On mutual impact of numerical linear algebra and large-scale optimization with focus on interior point methods.
Efficient numerical methods for computing ground states and ...
2010-08-06T23:59:59.000Z
Jul 8, 2010 ... New efficient and accurate numerical methods are proposed to ... ating integrals with high singularity and thus they are more efficient and ...
Numerical Modeling of the Nucleation Conditions of Petal-Centerline...
Coso Geothermal Field Jump to: navigation, search OpenEI Reference LibraryAdd to library Conference Proceedings: Numerical Modeling of the Nucleation Conditions of...
Preliminary LSG Numerical Reliability DATt 5 January 1971
Rathbun, Julie A.
of the component application analysis and the FMEA. I . f.: I #12;· .... ATM-937 HY.MO. Preliminary l.SG Numerical
A Numerical Evaluation Of Electromagnetic Methods In Geothermal...
Published Journal International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts, 1996 DOI 10.1016S0148-9062(97)87449-9 Citation . 1996. A Numerical...
Numerical Modeling At Dixie Valley Geothermal Area (McKenna ...
models were made to analyze the transient behavior of extensional basin and range geothermal systems. Notes Numerical models were made utilizing a range of bulk rock...
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)...
Numerical study of reaction in porous catalysts under composition modulation
Hsiao, Hsu-Wen
2010-01-01T23:59:59.000Z
Numerical methods for chemical engineering: applications incycling of the feed." Chemical Engineering Science 54(20): [computer." Computers & Chemical Engineering 26(3): 439-
Numerical Data Advisory Board assembly of mathematial and physical sciences
Not Available
1980-07-30T23:59:59.000Z
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.
Evolution of the Bohemian Massif: Insights from numerical modeling
Cerveny, Vlastislav
Evolution of the Bohemian Massif: Insights from numerical modeling Petra Maierová Supervisor: Doc of Geophysics Faculty of Mathematics and Physics Charles University in Prague #12;February 4, 2013Evolution Conclusions Outline #12;February 4, 2013Evolution of the Bohemian Massif: Insights from numerical modeling 3
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
Numerical Analysis of Nonlinear Multiharmonic Eddy Current Problems
Schoeberl, Joachim
Numerical Analysis of Nonlinear Multiharmonic Eddy Current Problems F. Bachinger U. Langer J. Sch¨oberl April 2004 Abstract This work provides a complete analysis of eddy current problems, rang- ing from at providing a complete analysis of general nonlinear eddy current problems and their numerical treatment
Numerical Prediction of High-Impact Local Weather: A
Xue, Ming
Chapter 6 Numerical Prediction of High-Impact Local Weather: A Driver for Petascale Computing Ming .................................................. 116 6.4 Numerical Weather Prediction Challenges and Requirements .......... 119 6.5 Summary The so-called mesoscale and convective scale weather events, including floods, tornadoes, hail, strong
Data Assimilation for Idealised Mathematical Models of Numerical Weather Prediction
Wirosoetisno, Djoko
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 in weather dynamics, e.g., the asymptotic balance seen in high and low pressure systems. Aims of the project
Numerical Modeling of Eastern Connecticut's Visual Resources1
Conservation, The University of Connecticut, Storrs, Connecticut. / Abstract: A numerical model capableNumerical Modeling of Eastern Connecticut's Visual Resources1 Daniel L. Civco 2/ l Submitted of accurately predicting the preference for landscape photographs of selected points in eastern Connecticut
Practical Experience in the Numerical Dangers of Heterogeneous Computing
Dongarra, Jack
in writing reliable numerical library software for heterogeneous computing environments. Although a lot ways to avoid potential pitfalls, or if that is not possible, we recommend that the software;1. INTRODUCTION There are special challenges associated with writing reliable numerical software on networks
Applications of Numerical Optimal Control to Nonlinear Hybrid Systems
Zefran, Milo?
Applications of Numerical Optimal Control to Nonlinear Hybrid Systems Shangming Wei a , Kasemsak the system has no mem- ory, the embedding technique can be used to reduce the hybrid optimal control problem of the controllers. Key words: Hybrid systems, optimal control, numerical optimization, model predictive control
A NUMERICAL STUDY OF A FULLY CONSERVATIVE METHOD FOR HYPERELASTICVISCOPLASTIC
New York at Stoney Brook, State University of
based on a new approach to the simulation of materials undergoing large de formation. Our numericalA NUMERICAL STUDY OF A FULLY CONSERVATIVE METHOD FOR HYPERELASTICVISCOPLASTIC MATERIALS XIAO LIN algorithm for the simulation of the impact of hyperelasticviscoplastic materials in two dimensions
Computation and Visualisation in the NumLab Numerical Laboratory
Telea, Alexandru C.
for the underlying mathematics. Assembling numerical simulations from computational and visualisation blocks, as wellComputation and Visualisation in the NumLab Numerical Laboratory J.M.L. Maubach1 and A.C. Telea1 Eindhoven University of Technology, Department of Mathematics and Computer Science, Postbox 513, NL-5600 MB
The Method of Generated Solutions for Numerical Verification of ICE
for verifying numerical code. The exact solution to the set of equations is derived using mathematical methodsThe Method of Generated Solutions for Numerical Verification of ICE Code Amjidanutpan Ramanujam, Christopher Sikorski, Todd Harman* UUCS07006 School of Computing University of Utah Salt Lake City, UT 84112
The Method of Generated Solutions for Numerical Verification of ICE
for verifying numerical code. The exact solution to the set of equations is derived using mathematical methodsThe Method of Generated Solutions for Numerical Verification of ICE Code Amjidanutpan Ramanujam, Christopher Sikorski, Todd Harman* UUCS-07-006 School of Computing University of Utah Salt Lake City, UT 84112
INCORPORATION OF GROUNDWATER FLOW INTO NUMERICAL MODELS AND DESIGN MODELS
03/10/99 1 INCORPORATION OF GROUNDWATER FLOW INTO NUMERICAL MODELS AND DESIGN MODELS Jeffrey D-coupled, ground-source heat pumps, groundwater, heat pump, heat exchanger, heat transfer, numerical models transport of heat by moving groundwater may be an important factor in reducing the necessary size of closed
Hamilton-Jacobi-Bellman Equations Analysis and Numerical Analysis
Flynn, E. Victor
Hamilton-Jacobi-Bellman Equations Analysis and Numerical Analysis Iain Smears #12;My deepest thanks at Durham University. Abstract. This work treats Hamilton-Jacobi-Bellman equations. Their relation and the inaugural papers on mean-field games. Original research on numerical methods for Hamilton-Jacobi
MULTISCALE NUMERICAL STUDY OF TURBULENT FLOW AND BUBBLE ENTRAINMENT
Kirby, James T.
MULTISCALE NUMERICAL STUDY OF TURBULENT FLOW AND BUBBLE ENTRAINMENT IN THE SURF ZONE BY GANGFENG MA . . . . . . . . . . . . . . . . . . . . . 3 1.1.2 Numerical Investigations . . . . . . . . . . . . . . . . . . . . . 5 1.2 Bubble Entrainment . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15 2.4 Bubble Entrainment Model . . . . . . . . . . . . . . . . . . . . . . . 18 2.5 Bubble Breakup
Numerical integration and importance 2.1 Quadrature
Schofield, Jeremy
2 Numerical integration and importance sampling 2.1 Quadrature Consider the numerical evaluation of the integral I(a, b) = b a dx f(x) · Rectangle rule: on small interval, construct interpolating function and integrate over interval. Polynomial of degree 0 using mid-point of interval: (a+1)h ah dx f(x) h f ((ah
Numerical Simulation of a Natural Circulation Steam Generator
Weinmüller, Ewa B.
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
NUMERICAL SIMULATIONS OF TRANSVERSE COMPRESSION AND DENSIFICATION IN WOOD
Nairn, John A.
NUMERICAL SIMULATIONS OF TRANSVERSE COMPRESSION AND DENSIFICATION IN WOOD John A. Nairn1 Professor- terials is a useful tool for stress analysis and for failure modeling. Although FEA of wood as an anisotropic continuum is used, numerical modeling of realistic wood structures, including details of wood
Numerical simulation of tsunami waves generated by deformable submarine landslides
Kirby, James T.
Numerical simulation of tsunami waves generated by deformable submarine landslides Gangfeng Ma a wave model Tsunami wave Numerical modeling a b s t r a c t This paper presents a new submarine of landslide motion and associated tsunami wave generation on parameters including sediment settling velocity
Numerical Simulation in Applied Geophysics. From the Mesoscale to the
Santos, Juan
Numerical Simulation in Applied Geophysics. From the Mesoscale to the Macroscale Juan E. Santos Numerical Simulation in Applied Geophysics. From the Mesoscale to the Macroscale p. #12;Introduction. I layering, fractures and craks at the mesoscale (on the order of centimeters) are common in the earth
Buoyant melting instabilities beneath extending lithosphere: 1. Numerical models
Tackley, Paul J.
Buoyant melting instabilities beneath extending lithosphere: 1. Numerical models John W. Hernlund,1,2 Paul J. Tackley,1,3 and David J. Stevenson4 Received 18 November 2006; revised 18 October 2007 diffusely extending lithosphere is studied using numerical convection models covering a wide range
NUMERICAL SIMULATION OF AIR POLLUTION DYNAMICS DUE TO
Olszewski Jr., Edward A.
NUMERICAL SIMULATION OF AIR POLLUTION DYNAMICS DUE TO POINT SOURCE EMISSIONS FROM AN INDUSTRIAL · Works Cited #12;Statement of the Problem · Numerical simulation of air pollution from a point source threat of air pollution Major health hazard to humans and ecosystem Billions effected, mainly
DYNAMOACTION,BETWEEN NUMERICAL EXPERIMENTS AND LIQUID SODIUM DEVICES
Guermond, Jean-Luc
and time, with a zero mean. A real fluid dynamo cannot be based on a purely laminar (steady) flow nor numerical dynamos are based either on laminar or on low kinetic Reynolds number flows. Moreover on the flow geometry. Positive numerical results are known either for ``laminar dynamos'', where the flow
Starting laminar plumes: Comparison of laboratory and numerical modeling
van Keken, Peter
Starting laminar plumes: Comparison of laboratory and numerical modeling Judith Vatteville Institut, France (davaille@fast.u-psud.fr) [1] A detailed comparison of starting laminar plumes in viscous fluids. Davaille (2009), Starting laminar plumes: Comparison of laboratory and numerical modeling, Geochem. Geophys
KH Computational Physics-2015 Basic Numerical Algorithms Ordinary differential equations
Gustafsson, Torgny
KH Computational Physics- 2015 Basic Numerical Algorithms Ordinary differential equations The set(xl) at certain points xl. Kristjan Haule, 2015 1 #12;KH Computational Physics- 2015 Basic Numerical Algorithms purpose routine · Numerov's algorithm: ¨y = f(t)y(t) ( for Schroedinger equation) · Verlet algorithm: ¨y
Sidorov, Nikita
The development of Numerical Analysis in Manchester was stimulated by the pioneering work on automatic computers on the first computer at Cambridge) to promote research in mathematical computing in Manchester. Some which required numerical solutions could benefit directly from faster computation. This was particularly
Grid-independent Issue in Numerical Heat Transfer
Yao Wei; Wang Jian; Liao Guangxuan
2006-09-26T23:59:59.000Z
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.
Simple intrinsic defects in GaAs : numerical supplement.
Schultz, Peter Andrew
2012-04-01T23:59:59.000Z
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.
Numerical Simulations of MHD Turbulence in Accretion Disks
Steven A. Balbus; John F. Hawley
2002-03-20T23:59:59.000Z
We review numerical simulations of MHD turbulence. The last decade has witnessed fundamental advances both in the technical capabilities of direct numerical simulation, and in our understanding of key physical processes. Magnetic fields tap directly into the free energy sources in a sufficiently ionized gas. The result is that adverse angular velocity and adverse temperature gradients, not the classical angular momentum and entropy gradients, destabilize laminar and stratified flow. This has profound consequences for astrophysical accretion flows, and has opened the door to a new era of numerical simulation experiments.}
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
2009-07-09T23:59:59.000Z
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.
Numerical investigation of the heating process inside an industrial furnace
Wolper, Pierre
Numerical investigation of the heating process inside an industrial furnace Proposition: Combined furnace taking into account convective, conductive and radiative heat transfer. The model: Catalysis, Energy Materials, Performance Materials and Recycling. Each business area is divided into market
Numerical simulation of flow separation control by oscillatory fluid injection
Resendiz Rosas, Celerino
2005-08-29T23:59:59.000Z
In this work, numerical simulations of flow separation control are performed. The sep-aration control technique studied is called 'synthetic jet actuation'. The developed code employs a cell centered finite volume scheme which handles viscous...
GEOMETRIC NUMERICAL INTEGRATION OF INEQUALITY CONSTRAINED, NONSMOOTH HAMILTONIAN SYSTEMS
Columbia University
and obtain an associated framework for composing ge- ometric numerical integration methods for inequality-equalityÂconstrained geometric integration setting, we target long-term structure preservation. In the inequality-constrained
Multiscale numerical methods for some types of parabolic equations
Nam, Dukjin
2009-05-15T23:59:59.000Z
method. The goal of the second problem is to develop efficient multiscale numerical techniques for solving turbulent diffusion equations governed by celluar flows. The solution near the separatrices can be approximated by the solution of a system of one...
Numerical modeling of fluid flow and time-lapse seismograms ...
gabriela
1. Inst. del Gas y del Petr´oleo - FI, Univ. de Buenos Aires, ARGENTINA. 2 .... 1 + i??s ). (9) where ?e > ?s are relaxation times and Mr = KG,µm. Numerical ...
Numerical Weather Forecasting at the Savannah River Site
Buckley, R.L. [Westinghouse Savannah River Company, AIKEN, SC (United States)
1998-08-01T23:59:59.000Z
This paper discusses the use of an advanced three-dimensional prognostic numerical model to provide space and time-dependent meteorological data for use in the WIND System dispersion models.
Numerical simulations of mine penetration in soft sediments
Yao, Zhigang
2003-01-01T23:59:59.000Z
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...
Characterization of Lung's Emphysema Distribution: Numerical Assessment of Disease Development
Paris-Sud XI, Université de
Characterization of Lung's Emphysema Distribution: Numerical Assessment of Disease Development M, Egypt. Abstract--Chronic Obstructive Pulmonary Disease (COPD) refers to a group of lung diseases bronchitis. Pulmonary emphysema is defined as a lung disease characterized by "abnormal enlargement
A numerical model simulation of longshore transport for Galveston Island
Gilbreath, Stephen Alexander
1995-01-01T23:59:59.000Z
The shoreline changes, deposition patterns, and longshore transport rates were calculated for the coast of Galveston Island using a numerical model simulation. The model only simulated changes due to waves creating longshore currents. East Beach...
Geometric numerical integration of nonholonomic systems and optimal control problems
M. de Leon; D. Martin de Diego; A. Santamaria Merino
2002-12-02T23:59:59.000Z
A geometric derivation of numerical integrators for nonholonomic systems and optimal control problems is obtained. It is based in the classical technique of generating functions adapted to the special features of nonholonomic systems and optimal control problems.
Assessment of the MUSTA approach for numerical relativistic dynamics
Blakely, P. M.; Nikiforakis, N.; Henshaw, W. D.
2015-03-04T23:59:59.000Z
waves 1. Introduction The numerical solution of the relativistic hydrodynamical equa- tions is of importance to the simulation of astrophysical phe- nomena such as gamma-ray bursts, supernova core-collapse, and relativistic wind accretion. Although...
Numerical Simulation Study on Transpired Solar Air Collector
Wang, C.; Guan, Z.; Zhao, X.; Wang, D.
2006-01-01T23:59:59.000Z
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...
Numerical and Experimental Investigation of Tidal Current Energy Extraction
Sun, Xiaojing
2008-01-01T23:59:59.000Z
Numerical and experimental investigations of tidal current energy extraction have been conducted in this study. A laboratory-scale water flume was simulated using commercial computational fluid dynamics (CFD) code FLUENT. ...
Numerical Simulation Study on Transpired Solar Air Collector
Wang, C.; Guan, Z.; Zhao, X.; Wang, D.
2006-01-01T23:59:59.000Z
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...
Prediction of Thermoacoustic Instabilities: Numerical Study of Mach number Effects.
Nicoud, Franck
Prediction of Thermoacoustic Instabilities: Numerical Study of Mach number Effects. K. Wieczorek equations. I. Introduction In the calculation of thermoacoustic instabilities, an assumption lead to significant changes in the evaluation of the thermo-acoustic modes present in the combustion
Householder Symposium on Numerical Linear Algebra June 1721, 2002
Higham, Nicholas J.
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
Effective velocities in fractured media: a numerical study using the ...
2002-03-18T23:59:59.000Z
theories that predict the effective P- and S-wave velocities in fractured materials in ... can be treated only by numerical techniques because analyt- ... apply the rotated staggered grid (Saenger, Gold and Shapiro ..... (r ? 0:2) and for SH-
An investigation of the numerical treatment of condensation
Sasson, Joseph
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 ...
Numerical Simulations of Gamma-Ray Burst Explosions
Lazzati, Davide; López-Cámara, Diego
2015-01-01T23:59:59.000Z
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...
Direct Numerical Simulation of the Flow in a Pebble Bed
Ward, Paul
2014-06-24T23:59:59.000Z
at Argonne National Laboratory, to conduct both large eddy simulation (LES) and direct numerical simulation (DNS) of fluid flow through a single face-centered cubic sphere lattice with periodic boundary conditions. Multiple LES were conducted with varying...
Numerical modelling of hyperbolic conservation laws using bicharacteristics
Hanke-Bourgeois, Martin
volume methods -1- #12;Overview I. Hyperbolic Conservation Laws Theory of bicharacteristics and evolution, bicharacteritsics stability, accuracy, error analysis III. Numerical Experiments: Wave equation system, Euler eqs finite volume methods -2- #12;· airflow · hydraulic schock · meteorological flow Source: efluid
Wavelet methods For the numerical simulation of incompressible fluids
Starck, Jean-Luc
Wavelet methods For the numerical simulation of incompressible fluids Erwan Deriaz Erwan Numerik, Seminar February 16th 2006 0-0 #12; Wavelets for the Navier-Stokes equations homogeneous or non homogeneous) With a wavelet discretization: £ ¦ ¨ 2 ©$# £ % ¦ ¨ 2 © & '' ( 1
NUMERICAL MODELING FOR MULTIPHASE INCOMPRESSIBLE FLOW WITH PHASE CHANGE
Abdou, Mohamed
NUMERICAL MODELING FOR MULTIPHASE INCOMPRESSIBLE FLOW WITH PHASE CHANGE Xiao-Yong Luo, Ming-Jiu Ni for multiphase flows. A con- tinuum surface force (CSF) tension model is used in the present cases. Phase change
Direct Numerical Simulation of the Flow in a Pebble Bed
Ward, Paul
2014-06-24T23:59:59.000Z
bed reactors: dust generation and scaling, proceedings of ICAPP 2012, Chicago, June 24–28, 2012 [3] A. Shams, F. Roelofs, EMJ. Komen, E. Baglietto, 2013. “Quasi-direct numerical simulation of a pebble bed configuration. Part I: Flow (velocity...
Numerical Approximations of Stochastic Optimal Stopping and Control Problems
Siska, David
2007-01-01T23:59:59.000Z
We study numerical approximations for the payoff function of the stochastic optimal stopping and control problem. It is known that the payoff function of the optimal stopping and control problem corresponds to the solution ...
Numerical simulation of flow separation control by oscillatory fluid injection
Resendiz Rosas, Celerino
2005-08-29T23:59:59.000Z
and Lomax. The application of synthetic jet actuators is based in their ability to energize the boundary layer, thereby providing signifcant increase in the lift coefficient. This has been corroborated experimentally and it is corroborated numerically...
Numerical models of phosphate esters in the Chattahoochee River
Haffey, Samuel Fraad, 1973-
2004-01-01T23:59:59.000Z
A numerical model was constructed to assess the magnitude of organophosphoric acid triester sinks in the Chattahoochee River and to identify concentration patterns downstream of point source discharges. The model was built ...
Numerical Simulation of Laminar Reacting Flows with Complex Chemistry
Bell, John B.
Numerical Simulation of Laminar Reacting Flows with Complex Chemistry M S Day and J B Bell Lawrence Simulation of Laminar Reacting Flows 2 1. Introduction Detailed modelling of time-dependent reacting ows
Experimental and Numerical Investigation of Aerosol Scavenging by Sprays
Goldmann, Andrew S.
2011-02-22T23:59:59.000Z
the poten- tial to substantially decrease the amount of radioactive material that could reach a population zone in a postulated accident scenario. An experimental and numerical investigation of airborne particulate scavenging by water sprays was conducted...
Nondestructive Testing of Overhead Transmission Lines: Numerical and Experimental Investigation
Kulkarni, Salil Subhash
2011-02-22T23:59:59.000Z
line is the primary reason for absence of a theoretical solution to get the analytical dispersion curves. The numerical results are then verified with experimental tests using a non-contact and broadband laser detection technique. The methodology...
Experimental and numerical study of primary consolidation of soft clay
Korchaiyapruk, Attasit, 1977-
2007-01-01T23:59:59.000Z
This research investigates, experimentally and numerically, the three main factors that differentiate predictions based on Terzaghi's one-dimensional consolidation theory from those observed in the laboratory and field. ...
CSE/Math 555: Numerical Optimization Techniques Course Announcement
Shontz, Suzanne M.
towards graduate students, researchers and faculty in · computer science and engineering · mathematics portfolios Description: The course will emphasize the design and mathematical analysis of numerical op: unconstrained optimization methods, automatic differentiation, nonlinear equations, constrained optimization
Numerical analysis of vapor flow in a micro heat pipe
Liu, Xiaoqin
1996-01-01T23:59:59.000Z
The vapor flow in a flat plate micro heat pipe with both uniform and linear heat flux boundary conditions has been numerically analyzed. For both types of boundary conditions, the Navier-Stokes equations with steady incompressible two...
An experimental and numerical study of wind turbine seismic behavior
Prowell, I.
2011-01-01T23:59:59.000Z
3.2.1 Description of Test Wind Turbine . . . . . .Figure 1.2: Components of a modern wind turbine . . . . . .Experimental and Numerical Seismic Response of a 65-kW Wind
Fully kinetic numerical modeling of a plasma thruster
Szabo, James Joseph, 1969-
2001-01-01T23:59:59.000Z
A Hall effect plasma thruster with conductive acceleration channel walls was numerically modeled using 2D3V Particle-in-Cell (PIC) and Monte-Carlo Collision (MCC) methodolo- gies. Electron, ion, and neutral dynamics were ...
Notes 09. Numerical evaluation of natural modes and frequencies
San Andres, Luis
2008-01-01T23:59:59.000Z
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... will give you an appreciation of what goes on inside such canned routines. The knowledge below will help you to make an intelligent choice when using or selecting one of the methods detailed. The references listed at the end of this document...
Numerical studies of a simple Coulomb blockade model
Shao, Jianfeng
1991-01-01T23:59:59.000Z
NUMERICAL STUDIES OF A SIMPLE COULOMB BLOCKADE MODEL A Thesis by JIANFENG SHAO Submitted to the Office of Graduate Studies of Texas ASM University in partial fulfillment of the requirements for the degree of MASTER OF SCIENCE August 1991... Major Subject: Physics NUMERICAL STUDIES OF A SIMPLE COULOMB BLOCKADE MODEL A Thesis by JIANFENG SHAO Approved as to style and content by: Roland E, Allen (Chair of Committee) /, 1 r oseph H. R s ( Member) Chin B. Su (Member) Richard L...
Segmented vs conventional numerals: legibility and long term retention
Hill, Steve Edgar
1971-01-01T23:59:59.000Z
SEGMENTED VS CONVENTIONAL NUMERALS: LEGIBILITY AND LONG TERM RETENTION A Thesis by STEVE EDGAR HILL Submitted to the Graduate College of Texas A8cM University in partial fulfillment of the requirements for the degree of MASTER OF SCIENCE... 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...
Multigrid methods for improving the variational data assimilation in numerical weather prediction
Kwak, Do Young
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 condition: numerical weather prediction, variational data assimilation, minimization procedure, multigrid methods, cell
X. Parallel and Distributed Scientific A Numerical Linear Algebra Problem Solving Environment
Dongarra, Jack
X. Parallel and Distributed Scientific Computing A Numerical Linear Algebra Problem SolvingQuality, Reusable, Mathematical Software : : : : : : : : : : : : : : 467 3. Automatic Generation of Tuned Numerical : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : 456 2. Numerical Linear Algebra Libraries : : : : : : : : : : : : : : : : : : : : : : : : : : : : 459
Numerical Methods for a Nonlinear BVP Arising in Physical Oceanography
Riccardo Fazio; Alessandra Jannelli
2013-10-08T23:59:59.000Z
In this paper we report and compare the numerical results for an ocean circulation model obtained by the classical truncated boundary formulation, the free boundary approach and a quasi-uniform grid treatment of the problem. We apply a shooting method to the truncated boundary formulation and finite difference methods to both the free boundary approach and the quasi-uniform grid treatment. Using the shooting method, supplemented by the Newton's iterations, we show that the ocean circulation model cannot be considered as a simple test case. In fact, for this method we are forced to use as initial iterate a value close to the correct missing initial condition in order to be able to get a convergent numerical solution. The reported numerical results allow us to point out how the finite difference method with a quasi-uniform grid is the less demanding approach and that the free boundary approach provides a more reliable formulation than the classical truncated boundary formulation.
The instanton method and its numerical implementation in fluid mechanics
Grafke, Tobias; Schäfer, Tobias
2015-01-01T23:59:59.000Z
A precise characterization of structures occurring in turbulent fluid flows at high Reynolds numbers is one of the last open problems of classical physics. In this review we discuss recent developments related to the application of instanton methods to turbulence. Instantons are saddle point configurations of the underlying path integrals. They are equivalent to minimizers of the related Freidlin-Wentzell action and known to be able to characterize rare events in such systems. While there is an impressive body of work concerning their analytical description, this review focuses on the question on how to compute these minimizers numerically. In a short introduction we present the relevant mathematical and physical background before we discuss the stochastic Burgers equation in detail. We present algorithms to compute instantons numerically by an efficient solution of the corresponding Euler-Lagrange equations. A second focus is the discussion of a recently developed numerical filtering technique that allows to...
NUMERICAL SIMULATION OF NATURAL GAS-SWIRL BURNER
Ala Qubbaj
2005-03-01T23:59:59.000Z
A numerical simulation of a turbulent natural gas jet diffusion flame at a Reynolds number of 9000 in a swirling air stream is presented. The numerical computations were carried out using the commercially available software package CFDRC. The instantaneous chemistry model was used as the reaction model. The thermal, composition, flow (velocity), as well as stream function fields for both the baseline and air-swirling flames were numerically simulated in the near-burner region, where most of the mixing and reactions occur. The results were useful to interpret the effects of swirl in enhancing the mixing rates in the combustion zone as well as in stabilizing the flame. The results showed the generation of two recirculating regimes induced by the swirling air stream, which account for such effects. The present investigation will be used as a benchmark study of swirl flow combustion analysis as a step in developing an enhanced swirl-cascade burner technology.
Tokamak elongation: how much is too much? II Numerical results
Lee, Jungpyo; Freidberg, Jeffrey P
2015-01-01T23:59:59.000Z
The analytic theory presented in Paper I is converted into a form convenient for numerical analysis. A fast and accurate code has been written using this numerical formulation. The results are presented by first defining a reference set of physical parameters based on experimental data from high performance discharges. Numerically obtained scaling relations of maximum achievable elongation versus inverse aspect ratio are obtained for various values of poloidal beta, wall radius and feedback capability parameter in ranges near the reference values. It is also shown that each value of maximum elongation occurs at a corresponding value of optimized triangularity, whose scaling is also determined as a function of inverse aspect ratio. The results show that the theoretical predictions of maximum elongation are slightly higher than experimental observations for high performance discharges as measured by high average pressure. The theoretical optimized triangularity values are noticeably lower. We suggest that the e...
Propagation of numerical noise in particle-in-cell tracking
Kesting, Frederik
2015-01-01T23:59:59.000Z
The particle-in-cell (PIC) method is the most used algorithm to perform self-consistent tracking of intense particle beams. For short term tracking the PIC method gives reliable results. However, the limited accuracy of the algorithm becomes important for long term tracking. It acts on single particles as a numerical noise and therefore ultimately leads to an artificial diffusion of the tracked particle beam. In this work we analyze the effect of PIC induced noise on the beam evolution. A detailed analysis of the particle tracker conjoint with the PIC solver leads to the concepts of correlated and decorrelated noise. For decorrelated numerical noise we derive a useful scaling law on the simulation parameters, that can be applied to mitigate artificial noise effects. Further, a strategy is proposed to avoid correlations in the numerical noise and therefore to assure the reliability of simulation studies.
Numerical Modeling At Coso Geothermal Area (2010) | 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 Home Page on DeliciousPlasmaP a gHigh Plains WindInformationNVN-079666New look toNotes from the CallNuclearNumericalNumerical
Numerical calculation of wave refraction by digital computer
Orr, Terry Edwin
1969-01-01T23:59:59.000Z
+~ 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...
10.34 Numerical Methods Applied to Chemical Engineering, Fall 2001
Beers, Kenneth J.
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 ...
ATS 680 A6: Applied Numerical Weather Prediction MW, 1:00-1:50 PM, ACRC Room 212B
, Parameterization Schemes: Keys to Understanding Numerical Weather Prediction Models, Cambridge University Press in class. Numerical model The primary numerical model that will be u
CGMurphi: Automatic Synthesis of Numerical Controllers for Nonlinear Hybrid Systems
Tronci, Enrico
CGMurphi: Automatic Synthesis of Numerical Controllers for Nonlinear Hybrid Systems Giuseppe Della, thus much work is being done to create automatic controller synthesis tools. When dealing with critical systems, most of the times such controllers are required to be optimal and robust, i.e., they must achieve
Numerical simulation of ceramic breeder pebble bed thermal creep behavior
Abdou, Mohamed
Numerical simulation of ceramic breeder pebble bed thermal creep behavior Alice Ying *, Hulin Huang Abstract The evolution of ceramic breeder pebble bed thermal creep deformation subjected to an external of ceramic breeder pebble beds under thermomechanical loads is necessary to ensure that the integrity of beds
Direct Numerical Simulations and Modeling of Jets in Crossflow
Mahesh, Krishnan
Direct Numerical Simulations and Modeling of Jets in Crossflow A THESIS SUBMITTED TO THE FACULTY. i #12;To my parents and my grandparents, and to Ramnath ii #12;Abstract Jets in crossflow are used to study the different aspects of round jets in a crossflow. The first problem studies
Control of Jets in Crossflow using Direct Numerical Simulations
Mahesh, Krishnan
Control of Jets in Crossflow using Direct Numerical Simulations A THESIS SUBMITTED TO THE FACULTY in crossflow by axial pulsing. Our main idea is that pulsing generates vortex rings; the effect of pulsing on jets in crossflow can therefore be explained by studying the behavior of vortex rings in crossflow
Validation of Erosion Modeling: Physical and Numerical Mehrad Kamalzare1
Franklin, W. Randolph
-3590 ABSTRACT The overall intent of this research is to develop numerical models of erosion of levees, dams is necessary for emergency plans for levee or dam breaches. Griffis, 2007 addressed the overall design and hydraulic shear stress. #12;2 Xu and Zhang (2009) found that in addition to soil type, the degree
Numerical Methods for the Inverse Nonlinear Fourier Transform
Civelli, Stella; Secondini, Marco
2015-01-01T23:59:59.000Z
We introduce a new numerical method for the computation of the inverse nonlinear Fourier transform and compare its computational complexity and accuracy to those of other methods available in the literature. For a given accuracy, the proposed method requires the lowest number of operations
Particle acceleration in solar flares: observations versus numerical simulations
Particle acceleration in solar flares: observations versus numerical simulations A O Benz, P C processes such as isotropization and magnetic trapping are made. Keywords: Particle acceleration, hard X. As the electric field of reconnection with possible parallel component capable of particle acceleration is limited
Numerical Modeling of Brain Dynamics in Traumatic Situations -Impulsive Translations
Burtscher, Martin
.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
Capturing Incomplete Information in Resource Allocation Problems through Numerical
Powell, Warren B.
Capturing Incomplete Information in Resource Allocation Problems through Numerical Patterns Arun information in resource allocation models. We use a popular goodness-of-fit measure known as the Cramer-Von Mises metric based on the empirical distribution function (EDF) as the foundation of our approach. We
Communication Avoiding and Overlapping for Numerical Linear Algebra
Touriño, Juan
Communication Avoiding and Overlapping for Numerical Linear Algebra Evangelos Georganas1, Jorge exascale systems, communication cost must be avoided or overlapped. Communication-avoiding 2.5D algorithms. Communication overlap attempts to hide messaging latency by pipelining messages and overlapping
Numerical Simulation of Titanium Production in the Plasma Quench Reactor
Numerical Simulation of Titanium Production in the Plasma Quench Reactor Ray A. Beny and Randall A the nucleation of condensates in the steady-state supersonic nozzle flow generated in a plasma quench reactor reactions. The device has been termed the Plasma Quench Reactor or PQR. The PQR has demonstrated the ability
Numerical Approximation of Vortex Density Evolution in a Superconductor.
Styles, Vanessa
Numerical Approximation of Vortex Density Evolution in a Superconductor. C.M. Elliott & V. Styles Abstract A #12;nite volume/element approximation of a mean #12;eld model of superconducting vortices in one approximations of a two-dimensional version of the mean #12;eld model of superconducting vortices considered
Probabilistic methods applied to 2D electromagnetic numerical dosimetry
Paris-Sud XI, Université de
the determinist solution is widely recast [Ghanem and Kruger, 1996] and the so called non intrusive method on the non intrusive method because of the limitation of the numerical resources. Indeed, the intrusive using the polyno- mial chaos. There are two approaches: the so called intrusive method where
Breakdown of analyticity: From rigorous results to numerics
Doedel, Eusebius
work with Rafael de la Llave #12;Computing the breakdown Quasi-Periodic Solutions Outline QuasiBreakdown of analyticity: From rigorous results to numerics Renato Calleja Applied Mathematics-Periodic Solutions Twist Maps Models arising in Statistical Mechanics Computation of the Breakdown Small divisors
Black hole free energy during charged collapse: a numerical study
Hugues Beauchesne; Ariel Edery
2012-05-19T23:59:59.000Z
We perform a numerical investigation of the thermodynamics during the collapse of a charged (complex) scalar field to a Reissner-Nordstr\\"om (RN) black hole in isotropic coordinates. Numerical work on gravitational collapse in isotropic coordinates has recently shown that the negative of the total Lagrangian approaches the Helmholtz free energy F= E-TS of a Schwarzschild black hole at late times of the collapse (where E is the black hole mass, T the temperature and S the entropy). The relevant thermodynamic potential for the RN black hole is the Gibbs free energy G=E-TS-$\\Phi_H$ Q where Q is the charge and $\\Phi_H$ the electrostatic potential at the outer horizon. In charged collapse, there is a large outgoing matter wave which prevents the exterior from settling quickly to a static state. However, the interior region is not affected significantly by the wave. We find numerically that the interior contribution to the Gibbs free energy is entirely gravitational and accumulates in a thin shell just inside the horizon. The entropy is gravitational in origin and one observes dynamically that it resides on the horizon. We also compare the numerical value of the interior Lagrangian to the expected analytical value of the interior Gibbs free energy for different initial states and we find that they agree to within 10-13%. The two values are approaching each other so that their difference decreases with more evolution time.
Conference: Numerical Analysis and Computers50 years of progress,
Higham, Nicholas J.
Report Conference: Numerical Analysis and Computers50 years of progress, University of Manchester, June 1617, 1998 A. Makroglou, School of Computer Science and Mathematics, Division of Mathematics and Computers50 years of progress ``Conference. It was held at the Mathematics Tower building, Oxford Rd
Conference: Numerical Analysis and Computers50 years of progress,
Higham, Nicholas J.
Report Conference: Numerical Analysis and Computers50 years of progress, University of Manchester, June 1617, 1998 A. Makroglou, School of Computer Science and Mathematics, Division of Mathematics and Computers50 years of progress'' Conference. It was held at the Mathematics Tower building, Oxford Rd
Geometric Integration: Numerical Solution of Differential Equations on Manifolds
Scheichl, Robert
and the solar system. Conserved quantities of a Hamiltonian system, such as energy, linear and angular momentumGeometric Integration: Numerical Solution of Differential Equations on Manifolds C.J. Budd 1 & A riches. Psalms 104:24 Since their introduction by Sir Isaac Newton, differential equations have played
Pyrolysis in Porous Media: Numerical Analysis and Comparison to Experiments.
Paris-Sud XI, Université de
by the combustor [7,8]. Carbon- or oxide-based composite materials are preferred for thermal, mechanical components in aerospace vehicles such as combustors in rocket engines. The numerical simulation = activation energy (J.mol-1 ) H = enthalpy of reaction (J.mol-1 ) effk = effective thermal conductivity
Numerical Construction of Likelihood Distributions and the Propagation of Errors
J. Swain; L. Taylor
1997-12-12T23:59:59.000Z
The standard method for the propagation of errors, based on a Taylor series expansion, is approximate and frequently inadequate for realistic problems. A simple and generic technique is described in which the likelihood is constructed numerically, thereby greatly facilitating the propagation of errors.
Numerical Methods for Controlled Hamilton-Jacobi-Bellman
Labahn, George
Numerical Methods for Controlled Hamilton-Jacobi-Bellman PDEs in Finance P.A. Forsyth , G. Labahn problems, leading to Hamilton-Jacobi-Bellman (HJB) or Hamilton-Jacobi-Bellman-Isaacs (HJBI) equations. We Hamilton-Jacobi-Bellman (HJB) PDEs, partial integro differential equations (PIDEs) or Hamilton-Jacobi
Numerical Methods for Controlled Hamilton-Jacobi-Bellman
Forsyth, Peter A.
Numerical Methods for Controlled Hamilton-Jacobi-Bellman PDEs in Finance P.A. Forsyth , G. Labahn problems, lead- ing to Hamilton-Jacobi-Bellman (HJB) or Hamilton-Jacobi-Bellman-Isaacs (HJBI) equations. We Hamilton-Jacobi-Bellman (HJB) PDEs, partial integro differential equations (PIDEs) or Hamilton-Jacobi
On Numerical Methods for Hyperbolic Conservation Laws and Related Equations
Bürger, Raimund
A classical kinematical model of sedimentation of small equal-sized particles dispersed in a viscous fluidOn Numerical Methods for Hyperbolic Conservation Laws and Related Equations Modelling Sedimentation with nonlocal flux, systems of nonlinear conservation modelling the sedimentation of polydisperse suspensions
Simulation of the secondary settling process with reliable numerical methods
Bürger, Raimund
are placed in alphabetical order Abstract A consistent model for the settling-compression-dispersion process The sedimentation process in the secondary settling tank (SST) is still a challenge in modelling the full Sedimentation, clarifier, conservation law, partial differential equation, numerical flux INTRODUCTION
Numerical calculation of wave refraction by digital computer
Orr, Terry Edwin
1969-01-01T23:59:59.000Z
. Navy Hydrographic Office, Pub. No. 603, 1963. 10. Selby, S. M. and Girling, M. S, C. , Standard Mathematical Tables, The Chemical Rubber Co. , Cleveland, Ohio, 1963. 15. McCormick, J. M. , Sa) vadori, M. G. , Numer ical Methods in Fortran, Prentice...
Numerical Convergence Properties of Option Pricing PDEs with Uncertain Volatility
Forsyth, Peter A.
is then simply the cost of this hedge. The principal source of risk is the price of the underlying assetNumerical Convergence Properties of Option Pricing PDEs with Uncertain Volatility D. M. Pooley , P, 2001 Abstract The pricing equations derived from uncertain volatility models in finance are often cast
NUMERICAL SIMULATIONS OF ROTATING SUNSPOTS G. J. J. Botha1
Rucklidge, Alastair
, at a depth of 0 to 5 Mm there exist subsurface horizontal vortical flows, while there are also suggestions's convection zone. The solution forms a central flux tube in the cylindrical numerical domain, with convection, vortical flow. As a result, the azimuthal velocity and mag- netic field have their maximum values close
NUMERICAL APPROXIMATION OF P-SYSTEMS WITH COULOMB FRICTIONAL DAMPING
in Scientific Computation, North Carolina State University, Raleigh, NC 27695-8205, USA. Current address: Sandia Computation, North Carolina State University, Raleigh, NC 27695-8205, USA (gremaud@unity.ncsu.edu). Partially analyzed. Numerical results show the asymptotic stress profiles to be distinct but "close" to those
Numerical Model of a Tensioner System and Flex Joint
Huang, Han
2013-07-27T23:59:59.000Z
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...
Numerical modelling and analysis of a room temperature magnetic
Numerical modelling and analysis of a room temperature magnetic refrigeration system Thomas Frank and analysis of a room temperature magnetic refrigeration system Department: Fuel Cells and Solid State-dimensional mathematical model of an Active Magnetic Regenerator (AMR) system which is used for magnetic refrigeration
Numerical Experimentation: A Third Way to Study Nature
École Normale Supérieure
of the numerical approach and trace back the origins of the use of computers to carry out simulations approach goes back much further than the appearance of the first computers. In a paper submitted in 1822 [1, and applied it to find out the atmospheric situation on May 20th 1910. After six Frontiers of Computational
Translation and integration of numerical atomic orbitals in linear molecules
Heinäsmäki, Sami, E-mail: sami.heinasmaki@gmail.com [Department of Physics, University of Oulu, FIN-90014, Oulu (Finland)] [Department of Physics, University of Oulu, FIN-90014, Oulu (Finland)
2014-02-14T23:59:59.000Z
We present algorithms for translation and integration of atomic orbitals for LCAO calculations in linear molecules. The method applies to arbitrary radial functions given on a numerical mesh. The algorithms are based on pseudospectral differentiation matrices in two dimensions and the corresponding two-dimensional Gaussian quadratures. As a result, multicenter overlap and Coulomb integrals can be evaluated effectively.
Numerical integration of functions originating from quantum mechanics
Armiento, Rickard
Numerical integration of functions originating from quantum mechanics R. Armiento Department Applications in quantum physics commonly involve large batches of integrals of smooth but very oscillatory for evaluating such integrals. The routines studied include: two from the QUADPACK package based on Gauss
On the numerical solution of some problems of environmental pollution
Ehrhardt, Matthias
world and their importance will even increase in the future. High pollution of air, water and soil mayOn the numerical solution of some problems of environmental pollution Quang A Dang1 , Matthias to the mathematical formula- tion of air pollution models. Let G be a cylindrical domain in the three
university-logo Numerical stability analysis for thin film flow
Marzuola, Jeremy
university-logo Numerical stability analysis for thin film flow: toward rigorous verification Blake Barker Indiana University October 2, 2013 B. Rigorous verification #12;university-logo Viscous roll waves (Picture courtesy Neil Balmforth, UBC.) B. Rigorous verification #12;university-logo Viscous roll waves 0 2
Numerical Modeling of Human Effect on Indoor Propagation
Paris-Sud XI, Université de
several types of obstacles such a brick enclosure walls, one metal heater, two metal wardrobes, two metal antenna placed 85 cm above the floor level. The spatial step is /10, or is the wavelength, chosen by the numerical dispersions. The choice of the spatial step is a compromise between the minimization of inaccuracy
Numerical modelling of tsunami mitigation by mangroves Putu Harry Gunawan
Mancini, Simona
Numerical modelling of tsunami mitigation by mangroves Putu Harry Gunawan LAMA (Laboratoire d'Analyse et de Mathmatiques Appliques) UPEM putu-harry.gunawan@univ.paris-est.fr Abstract Figure 1: Mangrove-Tsunami Model. The role of mangroves (coastal forests) in the mitigation of tsunami impacts is a debated topic
Numerical Simulation of the December 26, 2004: Indian Ocean Tsunami
Kirby, James T.
Numerical Simulation of the December 26, 2004: Indian Ocean Tsunami J. Asavanant1, M. Ioualalen2, N. Kaewbanjak1, S. Grilli3, P. Watts4, and J. Kirby5 Abstract: The December 26, 2004 tsunami is one of the most devastating tsunami in recorded history. It was generated in the Indian Ocean off the western coast
Diffusion and Dispersion Characterization of a Numerical Tsunami Model
Tolkova, Elena
Diffusion and Dispersion Characterization of a Numerical Tsunami Model David Burwell, Elena Tolkova and Arun Chawla February 6, 2007 1 Introduction Method Of Splitting Tsunami (MOST) is a depth averaged long wave tsunami inundation model that was originally developed by Titov and Synolakis (1995) for 1D
Simple intrinsic defects in InAs : numerical predictions.
Schultz, Peter Andrew
2013-03-01T23:59:59.000Z
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.
Friction versus dilation revisited: Insights from theoretical and numerical models
Einat, Aharonov
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
Numerical simulation of transpiration cooling through porous , T. Gotzen1
to facilitate such numerical simulations for a carbon/carbon material mounted in the side wall of a hot gas channel that are able to capture a spatially varying interplay between the hot gas flow and the coolant itself. This calls for new cooling technologies, for example, in the combustion chamber of the rocket
IVCNZ 2002 SUBMISSION 1 Pose Estimation by Applied Numerical Techniques
McCane, Brendan
process of how to deform the model in order to bring it into agreement with the input image. An energy. Keywords--- Numerical Optimisation, Pose Estimation, AnalysisÂbyÂSynthesis, Optimisation Algorithms (HCI). Traditionally, researchers have split approaches into an appearanceÂbased approach and a model
Numerical Modeling of the 2011 Tohoku Earthquake Tsunami
Frandsen, Jannette B.
Evolution of Ocean WAVEs) to reconstruct the tsunami across the Pacific and its transformation aroundNumerical Modeling of the 2011 Tohoku Earthquake Tsunami Yoshiki Yamazaki Post-doctoral Research 2011 Tohoku earthquake (Mw 9.0) generated a massive tsunami devastated the entire Pacific coast
Experimental and Numerical Investigation of DNAPL Dissolution Processes
Clement, Prabhakar
Experimental and Numerical Investigation of DNAPL Dissolution Processes in a Laboratory Aquifer) are extensively used in various industrial operations. Improper disposal and/or accidental release. The objective of this research was to complete an exper- imental investigation to study the dissolution patterns
Introduction Animals exploit numerous sources of information while
Lohmann, Kenneth J.
geographic position. For human navigators, this need is now usually met by a global positioning system (GPS3697 Introduction Animals exploit numerous sources of information while migrating, homing with two potential types of information. The simplest of these is directional or compass information, which
Numerical Simulation of the Wave Bottom Boundary Layer
Slinn, Donald
boundary layer. Oscillatory boundary layers are examined using a high-resolution time-dependent threeNumerical Simulation of the Wave Bottom Boundary Layer Over a Smooth Surface. Part 1: Three for turbulent boundary layers that occur over a smooth bottom. Results indicate that turbulence levels
ROOF CHARACTERISATION RELATED TO FIRE PROPAGATION RISK BY A NUMERICAL
Paris-Sud XI, Université de
ROOF CHARACTERISATION RELATED TO FIRE PROPAGATION RISK BY A NUMERICAL APPROACH L. Fournier1 , A by thé roof: - one is thé use of intumescent strips on thé roof, - thé other consists of extending thé fire walls (typically 70 cm or 1 m) above thé roof levé1 in order to prevent thé flame from being blown
A Numerical Investigation of Wettability Alteration during Immiscible CO2
Hossain, M. Enamul
A Numerical Investigation of Wettability Alteration during Immiscible CO2 Flooding Process was devoted to laboratory experiments including core and micro-model flooding to investigate wettability that the inclusion of wettability alteration is very influential parameter in enhancing oil recovery and sweeping
142 Los Alamos Science Number 29 2005 Direct Numerical
Kurien, Susan
142 Los Alamos Science Number 29 2005 Direct Numerical Simulations of Turbulence Data Generation visualization tools. Los Alamos scientists have been able to simulate flows of Reynolds numbers up to 105 subdomain of the 20483 turbulence simulation performed on the ASC Q machine at Los Alamos. The ASC Q machine
Numerical Study of Flexible Flapping Wing Propulsion , Mingjun Wei
Wei, Mingjun
Numerical Study of Flexible Flapping Wing Propulsion Tao Yang , Mingjun Wei New Mexico State Air Vehicles (MAVs), there is renewed interest in the flapping-wing design. At low Reynolds number-wing design offers other desirable characteristics in maneuverability, flutter, and energy efficiency.2
Hydraulic Fractures: multiscale phenomena, asymptotic and numerical solutions
Peirce, Anthony
Hydraulic Fractures: multiscale phenomena, asymptotic and numerical solutions SANUM Conference (UMN) Eduard Siebrits (SLB) #12;2 Outline · Examples of hydraulic fractures · Governing equations Energy Breaking rock Leak-off Viscous energy loss #12;10 2-3D HF Equations · Elasticity · Lubrication
Electron beam kinetics: numerical results Discussion of the experiments
Zharkova, Valentina V.
Electron beam kinetics: numerical results Discussion of the experiments In all the experiments the first set of Figures presents the differential energy spectra dN/dE for electron beams at a given depth presents the beam's mean electron flux. For comparison all the results for Experiments 1-3 are presented
Small Buoys for Energy Harvesting : Experimental and Numerical Modeling Studies
Grilli, Stéphan T.
Small Buoys for Energy Harvesting : Experimental and Numerical Modeling Studies St´ephan T. Grilli for wave energy harvesting (free-floating or slackly moored), to produce about 1 KW per unit at full scale-contained (water tight) resonating multiple-spar buoy (or Starspar), in which a longer central spar houses the LEG
Numerical simulation of turbulent jet primary breakup in Diesel engines
Helluy, Philippe
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
RisR1252(EN) The Numerical Wind Atlas
Risø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 and Application Program, WAsP, to make local predictions of the wind resource is presented. It combines
EXPERIMENTAL AND NUMERICAL STUDY OF TRANSIENT ELECTRONIC CHIP COOLING BY
EXPERIMENTAL AND NUMERICAL STUDY OF TRANSIENT ELECTRONIC CHIP COOLING BY LIQUID FLOW State University, Columbus, Ohio, USA Cooling of electronic chips has become a critical aspect on the transient electronic chip cooling by liquid flow in microchannel heat sinks. The temperature change
Numerical and experimental investigation of nonlinear properties of rubber absorber
Paris-Sud XI, Université de
Numerical and experimental investigation of nonlinear properties of rubber absorber in rail, Shanghai, 201804 Abstract: A dynamic nonlinear model of rubber absorber in railway fastening is proposed in this paper based on a superposition principal demonstrating that the restoring force of rubber components
Carmignani, B
2005-01-01T23:59:59.000Z
Numerical simulations of welds of thick steel pieces of interest for the thermonuclear fusion ITER machine
Numerical Analysis of Non-Uniformities and Anisotropy in High-Temperature Superconducting Coils
Hu, Di; Ainslie, Mark D.; Zou, Jin; Cardwell, David A.
2014-12-12T23:59:59.000Z
, critical current density (superconductivity), high-temperature superconductors, numerical analysis, superconducting coils, transport ac loss....
Olshanskii, Maxim A.
and accurate numerical methods for computing flows with free surfaces and interfaces, see, e.g., [1, 2 is studied in a series of numerical experiments. Institute of Numerical Mathematics, Russian Academy@math.uh.edu Department of Energy Resources Engineering, Stanford University and Institute of Numerical Mathematics
Basalt-Atmosphere Interactions on Venus -
Treiman, Allan H.
Atmosphere? Â· CaAl2Si2O8 + SO3 CaSO4 + Al2SiO5 + SiO2 Â Anhydrite + andalusite + quartz !!! Â· Does this reaction proceed? Â Venus atmosphere est'd 0.2 - 0.3 ppt SO3. #12;Reaction Position Â· From this, very possible that SO3 is buffered! Â· But ... SO3 value is not known very well - who knows what is really going
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 Home Page on DeliciousPlasmaP a gHigh4-FD-a <AlaskaAppliedInformation PageInformation source History View New
Hamann, Bernd
INTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS Int. J. Numer. Anal geomechanics Boris Jeremi!cc1, *,y , Gerik Scheuermann2,3 , Jan Frey3 , Zhaohui Yang1 , Bernd Hamann2 , Kenneth is the visualization of stress tensors resulting from 3D numerical simulations in computational geomechanics
Numerical Regularization of Electromagnetic Quantum Fluctuations in Inhomogeneous Dielectric Media
Shin-itiro Goto; Alison C. Hale; Robin W. Tucker; Timothy J. Walton
2012-01-05T23:59:59.000Z
Electromagnetic Casimir stresses are of relevance to many technologies based on mesoscopic devices such as MEMS embedded in dielectric media, Casimir induced friction in nano-machinery, micro-fluidics and molecular electronics. Computation of such stresses based on cavity QED generally require numerical analysis based on a regularization process. A new scheme is described that has the potential for wide applicability to systems involving realistic inhomogeneous media. From a knowledge of the spectrum of the stationary modes of the electromagnetic field the scheme is illustrated by estimating numerically the Casimir stress on opposite faces of a pair of perfectly conducting planes separated by a vacuum and the change in this result when the region between the plates is filled with an incompressible inhomogeneous non-dispersive dielectric.
Numerical analysis of nanostructures for enhanced light extraction from OLEDs
Zschiedrich, L; Burger, S; Schmidt, F; 10.1117/12.2001132
2013-01-01T23:59:59.000Z
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.
Noncommutative Field Theory: Numerical Analysis with the Fuzzy Disc
Fedele Lizzi; Bernardino Spisso
2012-07-28T23:59:59.000Z
The fuzzy disc is a discretization of the algebra of functions on the two dimensional disc using finite matrices which preserves the action of the rotation group. We define a $\\varphi^4$ scalar field theory on it and analyze numerically for three different limits for the rank of the matrix going to infinity. The numerical simulations reveal three different phases: uniform and disordered phases already the present in the commutative scalar field theory and a nonuniform ordered phase as a noncommutative effects. We have computed the transition curves between phases and their scaling. This is in agreement with studies on the fuzzy sphere, although the speed of convergence for the disc seems to be better. We have performed also three the limits for the theory in the cases of the theory going to the commutative plane or commutative disc. In this case the theory behaves differently, showing the intimate relationship between the nonuniform phase and noncommutative geometry.
On Numerical Considerations for Modeling Reactive Astrophysical Shocks
Papatheodore, Thomas L [ORNL] [ORNL; Messer, Bronson [ORNL] [ORNL
2014-01-01T23:59:59.000Z
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.
Numerical Estimation of Frictional Torques with Rate and State Friction
Arun K. Singh; T. N. Singh
2015-01-20T23:59:59.000Z
In this paper, numerical estimation of frictional torques is carried out of a rotary elastic disc on a hard and rough surface under different rotating conditions. A one dimensional spring- mass rotary system is numerically solved under the quasistatic condition with the rate and state dependent friction model. It is established that torque of frictional strength as well as torque of steady dynamic stress increases with radius and found to be maximum at the periphery of the disc. Torque corresponding to frictional strength estimated using the analytical solution matches closely with the simulation only in the case of high stiffness of the connecting spring. In steady relaxation simulation, a steadily rotating disc is suddenly stopped and relaxational angular velocity and corresponding frictional torque decreases with both steady angular velocity and stiffness of the connecting spring in the velocity strengthening regime. In velocity weakening regime, in contrast, torque of relaxation stress deceases but relaxation velocity increases. The reason for the contradiction is explained.
Numerical method for shear bands in ductile metal with inclusions
Plohr, Jee Yeon N [Los Alamos National Laboratory; Plohr, Bradley J [Los Alamos National Laboratory
2010-01-01T23:59:59.000Z
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.
Accurate complex scaling of three dimensional numerical potentials
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)
2013-05-28T23:59:59.000Z
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.
Fractional Calculus in Hydrologic Modeling: A Numerical Perspective
David A. Benson; Mark M. Meerschaert; Jordan Revielle
2012-01-01T23:59:59.000Z
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.
Refined numerical models for multidimensional Type Ia supernova simulations
Reinecke, M; Niemeyer, J C
2002-01-01T23:59:59.000Z
Following up on earlier work on this topic (Reinecke et al. 1999, A&A 347, pp. 724 and 739), we present an improved set of numerical models for simulations of white dwarfs exploding as Type Ia supernovae (SNe Ia). Two-dimensional simulations were used to test the reliability and numerical robustness of these algorithms; the results indicate that integral quantities like the total energy release are insensitive to changes of the grid resolution (above a certain threshold), which was not the case for our former code. The models were further enhanced to allow fully three-dimensional simulations of SNe Ia. A direct comparison of a 2D and a 3D calculation with identical initial conditions shows that the explosion is considerably more energetic in three dimensions; this is most likely caused by the assumption of axisymmetry in 2D, which inhibits the growth of flame instabilities in the azimuthal direction and thereby decreases the flame surface.
Refined numerical models for multidimensional Type Ia supernova simulations
M. Reinecke; W. Hillebrandt; J. C. Niemeyer
2001-11-26T23:59:59.000Z
Following up on earlier work on this topic (Reinecke et al. 1999, A&A 347, pp. 724 and 739), we present an improved set of numerical models for simulations of white dwarfs exploding as Type Ia supernovae (SNe Ia). Two-dimensional simulations were used to test the reliability and numerical robustness of these algorithms; the results indicate that integral quantities like the total energy release are insensitive to changes of the grid resolution (above a certain threshold), which was not the case for our former code. The models were further enhanced to allow fully three-dimensional simulations of SNe Ia. A direct comparison of a 2D and a 3D calculation with identical initial conditions shows that the explosion is considerably more energetic in three dimensions; this is most likely caused by the assumption of axisymmetry in 2D, which inhibits the growth of flame instabilities in the azimuthal direction and thereby decreases the flame surface.
Numerical solution of plasma fluid equations using locally refined grids
Colella, P., LLNL
1997-01-26T23:59:59.000Z
This paper describes a numerical method for the solution of plasma fluid equations on block-structured, locally refined grids. The plasma under consideration is typical of those used for the processing of semiconductors. The governing equations consist of a drift-diffusion model of the electrons and an isothermal model of the ions coupled by Poisson's equation. A discretization of the equations is given for a uniform spatial grid, and a time-split integration scheme is developed. The algorithm is then extended to accommodate locally refined grids. This extension involves the advancement of the discrete system on a hierarchy of levels, each of which represents a degree of refinement, together with synchronization steps to ensure consistency across levels. A brief discussion of a software implementation is followed by a presentation of numerical results.
Hyperboloidal Slices and Artificial Cosmology for Numerical Relativity
Charles W. Misner
2006-03-11T23:59:59.000Z
This preliminary report proposes integrating the Maxwell equations in Minkowski spacetime using coordinates where the spacelike surfaces are hyperboloids asymptotic to null cones at spatial infinity. The space coordinates are chosen so that Scri+ occurs at a finite coordinate and a smooth extension beyond Scri+ is obtained. The question addressed is whether a Cauchy evolution numerical integration program can be easily modified to compute this evolution. In the spirit of the von Neumann and Richtmyer artificial viscosity which thickens a shock by many orders of magnitude to facilitate numerical simulation, I propose artificial cosmology to thicken null infinity Scri+ to approximate it by a de Sitter cosmological horizon where, in conformally compactified presentation, it provides a shell of purely outgoing null cones where asymptotic waves can be read off as data on a spacelike pure outflow outer boundary. This should be simpler than finding Scri+ as an isolated null boundary or imposing outgoing wave conditions at a timelike boundary at finite radius.
Non-Gaussian numerical errors versus mass hierarchy
Y. Meurice; M. B. Oktay
2000-05-12T23:59:59.000Z
We probe the numerical errors made in renormalization group calculations by varying slightly the rescaling factor of the fields and rescaling back in order to get the same (if there were no round-off errors) zero momentum 2-point function (magnetic susceptibility). The actual calculations were performed with Dyson's hierarchical model and a simplified version of it. We compare the distributions of numerical values obtained from a large sample of rescaling factors with the (Gaussian by design) distribution of a random number generator and find significant departures from the Gaussian behavior. In addition, the average value differ (robustly) from the exact answer by a quantity which is of the same order as the standard deviation. We provide a simple model in which the errors made at shorter distance have a larger weight than those made at larger distance. This model explains in part the non-Gaussian features and why the central-limit theorem does not apply.
On numerical considerations for modeling reactive astrophysical shocks
Papatheodore, Thomas L.; Messer, O. E. Bronson, E-mail: tpapathe@utk.edu, E-mail: bronson@ornl.gov [Department of Physics and Astronomy, University of Tennessee, Knoxville, TN 37996 (United States)
2014-02-10T23:59:59.000Z
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. 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. 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 models of Type Ia supernovae, including potential impacts on observables.
Numerical Verification of Bounce Harmonic Resonances in Neoclassical
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), 2002; Level:5(Million Cubic Feet) Oregon (Including Vehicle Fuel) (Million Cubic Feet) Natural Gas DeliveredStents -(NEI) ScienceNumerical simulations of
Numerical analysis of a microwave torch with axial gas injection
Gritsinin, S. I.; Davydov, A. M.; Kossyi, I. A., E-mail: kossyi@fpl.gpi.ru [Russian Academy of Sciences, Prokhorov General Physics Institute (Russian Federation); Kulumbaev, E. B. [National Research University Belgorod State University (Russian Federation)] [National Research University Belgorod State University (Russian Federation); Lelevkin, V. M. [Kyrgyz-Russian Slavic University (Kyrgyzstan)] [Kyrgyz-Russian Slavic University (Kyrgyzstan)
2013-07-15T23:59:59.000Z
The characteristics of a microwave discharge in an argon jet injected axially into a coaxial channel with a shortened inner electrode are numerically analyzed using a self-consistent equilibrium gas-dynamic model. The specific features of the excitation and maintenance of the microwave discharge are determined, and the dependences of the discharge characteristics on the supplied electromagnetic power and gas flow rate are obtained. The calculated results are compared with experimental data.
Numerical studies of galaxy formation using special purpose hardware
Matthias Steinmetz
2002-01-25T23:59:59.000Z
I review recent progress in numerically simulating the formation and evolution of galaxies in hierarchically clustering universes. Special emphasis is given to results based on high-resolution gas dynamical simulations using the N-body hardware integrator GRAPE. Applications address the origin of the spin of disk galaxies, the structure and kinematics of damped Lyman-alpha systems, and the origin of galaxy morphology and of galaxy scaling laws.
Efficiency improvements for the numerical computation of NLO corrections
Sebastian Becker; Christian Reuschle; Stefan Weinzierl
2012-06-25T23:59:59.000Z
In this paper we discuss techniques, which lead to a significant improvement of the efficiency of the Monte Carlo integration, when one-loop QCD amplitudes are calculated numerically with the help of the subtraction method and contour deformation. The techniques discussed are: holomorphic and non-holomorphic division into sub-channels, optimisation of the integration contour, improvement of the ultraviolet subtraction terms, importance sampling and antithetic variates in loop momentum space, recurrence relations.
Numerical Simulation of Groundwater Withdrawal at the Nevada Test Site
Carroll, Rosemary; Giroux, Brian; Pohll, Greg; Hershey, Ronald; Russell, Charles; Howcroft, William
2004-01-28T23:59:59.000Z
Alternative uses of the Nevada Test Site (NTS) may require large amounts of water to construct and/or operate. The only abundant source of water at the NTS is groundwater. This report describes preliminary modeling to quantify the amount of groundwater available for development from three hydrographic areas at the NTS. Modeling was conducted with a three-dimensional transient numerical groundwater flow model.
Transient productivity index for numerical well test simulations
Blanc, G.; Ding, D.Y.; Ene, A. [Institut Francais du Petrole, Pau (France)] [and others
1997-08-01T23:59:59.000Z
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.
Numerical Methods of Computational Electromagnetics for Complex Inhomogeneous Systems
Cai, Wei
2014-05-15T23:59:59.000Z
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.
Numerical Simulation of the Hydrodynamical Combustion to Strange Quark Matter
Brian Niebergal; Rachid Ouyed; Prashanth Jaikumar
2010-08-27T23:59:59.000Z
We present results from a numerical solution to the burning of neutron matter inside a cold neutron star into stable (u,d,s) quark matter. Our method solves hydrodynamical flow equations in 1D with neutrino emission from weak equilibrating reactions, and strange quark diffusion across the burning front. We also include entropy change due to heat released in forming the stable quark phase. Our numerical results suggest burning front laminar speeds of 0.002-0.04 times the speed of light, much faster than previous estimates derived using only a reactive-diffusive description. Analytic solutions to hydrodynamical jump conditions with a temperature dependent equation of state agree very well with our numerical findings for fluid velocities. The most important effect of neutrino cooling is that the conversion front stalls at lower density (below approximately 2 times saturation density). In a 2-dimensional setting, such rapid speeds and neutrino cooling may allow for a flame wrinkle instability to develop, possibly leading to detonation.
Notes on the integration of numerical relativity waveforms
Christian Reisswig; Denis Pollney
2011-09-22T23:59:59.000Z
A primary goal of numerical relativity is to provide estimates of the wave strain, $h$, from strong gravitational wave sources, to be used in detector templates. The simulations, however, typically measure waves in terms of the Weyl curvature component, $\\psi_4$. Assuming Bondi gauge, transforming to the strain $h$ reduces to integration of $\\psi_4$ twice in time. Integrations performed in either the time or frequency domain, however, lead to secular non-linear drifts in the resulting strain $h$. These non-linear drifts are not explained by the two unknown integration constants which can at most result in linear drifts. We identify a number of fundamental difficulties which can arise from integrating finite length, discretely sampled and noisy data streams. These issues are an artifact of post-processing data. They are independent of the characteristics of the original simulation, such as gauge or numerical method used. We suggest, however, a simple procedure for integrating numerical waveforms in the frequency domain, which is effective at strongly reducing spurious secular non-linear drifts in the resulting strain.
Numerical implication of Riemann problem theory for fluid dynamics
Menikoff, R.
1988-01-01T23:59:59.000Z
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.
Towards a novel wave-extraction method for numerical relativity
Andrea Nerozzi; Marco Bruni; Lior M. Burko; Virginia Re
2006-07-17T23:59:59.000Z
We present the recent results of a research project aimed at constructing a robust wave extraction technique for numerical relativity. Our procedure makes use of Weyl scalars to achieve wave extraction. It is well known that, with a correct choice of null tetrad, Weyl scalars are directly associated to physical properties of the space-time under analysis in some well understood way. In particular it is possible to associate $\\Psi_4$ with the outgoing gravitational radiation degrees of freedom, thus making it a promising tool for numerical wave--extraction. The right choice of the tetrad is, however, the problem to be addressed. We have made progress towards identifying a general procedure for choosing this tetrad, by looking at transverse tetrads where $\\Psi_1=\\Psi_3=0$. As a direct application of these concepts, we present a numerical study of the evolution of a non-linearly disturbed black hole described by the Bondi--Sachs metric. This particular scenario allows us to compare the results coming from Weyl scalars with the results coming from the news function which, in this particular case, is directly associated with the radiative degrees of freedom. We show that, if we did not take particular care in choosing the right tetrad, we would end up with incorrect results.
Stochastic algorithms for the analysis of numerical flame simulations
Bell, John B.; Day, Marcus S.; Grcar, Joseph F.; Lijewski, Michael J.
2004-04-26T23:59:59.000Z
Recent progress in simulation methodologies and high-performance parallel computers have made it is possible to perform detailed simulations of multidimensional reacting flow phenomena using comprehensive kinetics mechanisms. As simulations become larger and more complex, it becomes increasingly difficult to extract useful information from the numerical solution, particularly regarding the interactions of the chemical reaction and diffusion processes. In this paper we present a new diagnostic tool for analysis of numerical simulations of reacting flow. Our approach is based on recasting an Eulerian flow solution in a Lagrangian frame. Unlike a conventional Lagrangian view point that follows the evolution of a volume of the fluid, we instead follow specific chemical elements, e.g., carbon, nitrogen, etc., as they move through the system . From this perspective an ''atom'' is part of some molecule of a species that is transported through the domain by advection and diffusion. Reactions cause the atom to shift from one chemical host species to another and the subsequent transport of the atom is given by the movement of the new species. We represent these processes using a stochastic particle formulation that treats advection deterministically and models diffusion and chemistry as stochastic processes. In this paper, we discuss the numerical issues in detail and demonstrate that an ensemble of stochastic trajectories can accurately capture key features of the continuum solution. The capabilities of this diagnostic are then demonstrated by applications to study the modulation of carbon chemistry during a vortex-flame interaction, and the role of cyano chemistry in rm NO{sub x} production for a steady diffusion flame.
Numerical simulation of three-dimensional electrical flow through geomaterials
Akhtar, Anwar Saeed
1998-01-01T23:59:59.000Z
investigation results further validate this fact. DEDICATION To my parents, who desired it more for me than I did, and to my dearest wife Sofia and loving daughter Aneeka, who endured most of it. ACKNOWLEDGMENTS The author would like to express his... 46 46 50 50 50 51 52 53 54 57 57 58 59 60 62 64 65 67 70 CHAPTER IV ELECTRICAL FLOW THROUGH GEOMATERIALS WITH ELECTRODES AT GROUND SURFACE Page 73 4. 1 INTRODUCTION 4. 2 NUMERICAL INVESTIGATION OF ELECTRICAL FLOW THROUGH...
Advanced numerical methods in mesh generation and mesh adaptation
Lipnikov, Konstantine [Los Alamos National Laboratory; Danilov, A [MOSCOW, RUSSIA; Vassilevski, Y [MOSCOW, RUSSIA; Agonzal, A [UNIV OF LYON
2010-01-01T23:59:59.000Z
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.
ASSIMILATION OF DOPPLER RADAR DATA INTO NUMERICAL WEATHER MODELS
Chiswell, S.; Buckley, R.
2009-01-15T23:59:59.000Z
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.
High numerical aperture projection system for extreme ultraviolet projection lithography
Hudyma, Russell M. (San Ramon, CA)
2000-01-01T23:59:59.000Z
An optical system is described that is compatible with extreme ultraviolet radiation and comprises five reflective elements for projecting a mask image onto a substrate. The five optical elements are characterized in order from object to image as concave, convex, concave, convex, and concave mirrors. The optical system is particularly suited for ring field, step and scan lithography methods. The invention uses aspheric mirrors to minimize static distortion and balance the static distortion across the ring field width which effectively minimizes dynamic distortion. The present invention allows for higher device density because the optical system has improved resolution that results from the high numerical aperture, which is at least 0.14.
Worldline Numerics for Energy-Momentum Tensors in Casimir Geometries
Schafer, Marco; Gies, Holger
2015-01-01T23:59:59.000Z
We develop the worldline formalism for computations of composite operators such as the fluctuation induced energy-momentum tensor. As an example, we use a fluctuating real scalar field subject to Dirichlet boundary conditions. The resulting worldline representation can be evaluated by worldline Monte-Carlo methods in continuous spacetime. We benchmark this worldline numerical algorithm with the aid of analytically accessible single-plate and parallel-plate Casimir configurations, providing a detailed analysis of statistical and systematic errors. The method generalizes straightforwardly to arbitrary Casimir geometries and general background potentials.
Diffusive mesh relaxation in ALE finite element numerical simulations
Dube, E.I.
1996-06-01T23:59:59.000Z
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.
Numerical analysis of decoy state quantum key distribution protocols
Harrington, Jim W [Los Alamos National Laboratory; Rice, Patrick R [Los Alamos National Laboratory
2008-01-01T23:59:59.000Z
Decoy state protocols are a useful tool for many quantum key distribution systems implemented with weak coherent pulses, allowing significantly better secret bit rates and longer maximum distances. In this paper we present a method to numerically find optimal three-level protocols, and we examine how the secret bit rate and the optimized parameters are dependent on various system properties, such as session length, transmission loss, and visibility. Additionally, we show how to modify the decoy state analysis to handle partially distinguishable decoy states as well as uncertainty in the prepared intensities.
Direct numerical integration for multi-loop integrals
Sebastian Becker; Stefan Weinzierl
2013-03-18T23:59:59.000Z
We present a method to construct a suitable contour deformation in loop momentum space for multi-loop integrals. This contour deformation can be used to perform the integration for multi-loop integrals numerically. The integration can be performed directly in loop momentum space without the introduction of Feynman or Schwinger parameters. The method can be applied to finite multi-loop integrals and to divergent multi-loop integrals with suitable subtraction terms. The algorithm extends techniques from the one-loop case to the multi-loop case. Examples at two and three loops are discussed explicitly.
Numerical solution of the linear least squares problem for splines
Vonderhaar, Thomas Jerome
1975-01-01T23:59:59.000Z
) Ci?. . &: , i ~ '(8'??. (Member) (Member) May 1975 ABSTRACT Numerical Solution of the Linear Least Squares Problem for Splines. (May 1975) Thomas Jerome Vonderhaar B. A. , Lewis University Chairman of Advisory Committee: Dr. Philip W. Smith... in the preparation of this thesis. I would also like to thank Dr. C. K. Chui and Dr. R. D. Larsen for serving on my committee. My sincere gratitude goes to Myra McWatt for her understanding and moral support during the completion of this paper. TABLE OF CONTENTS...
Numerical Modeling At Coso Geothermal Area (1995) | 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 Home Page on DeliciousPlasmaPLawrence County,1980) | OpenAl., 2001) |NewNorthNorwalk,Novomer Jump to:NuPowerNumerical Modeling At
Numerical Modeling At Dixie Valley Geothermal Area (Benoit, 1999) | 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 Center Home Page on DeliciousPlasmaPLawrence County,1980) | OpenAl., 2001) |NewNorthNorwalk,Novomer Jump to:NuPowerNumerical Modeling
Numerical Modeling At Dixie Valley Geothermal Area (Iovenitti, Et Al.,
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 Home Page on DeliciousPlasmaPLawrence County,1980) | OpenAl., 2001) |NewNorthNorwalk,Novomer Jump to:NuPowerNumerical Modeling2013)
Numerical Modeling At Raft River Geothermal Area (1983) | Open Energy
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 Home Page on DeliciousPlasmaPLawrence County,1980) | OpenAl., 2001) |NewNorthNorwalk,Novomer Jump to:NuPowerNumerical1984)
Numerical Modeling of Transient Basin and Range Extensional Geothermal
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 Home Page on DeliciousPlasmaPLawrence County,1980) | OpenAl., 2001) |NewNorthNorwalk,Novomer Jump to:NuPowerNumerical1984)Systems |
Numerical Modelling of Geothermal Systems a Short Introduction | 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 Center Home Page on DeliciousPlasmaPLawrence County,1980) | OpenAl., 2001) |NewNorthNorwalk,Novomer Jump to:NuPowerNumerical1984)Systems
Numerical studies of the stochastic Korteweg-de Vries equation
Lin Guang [Division of Applied Mathematics, Brown University, 182 George Street, Box F, Providence, RI 02912 (United States); Grinberg, Leopold [Division of Applied Mathematics, Brown University, 182 George Street, Box F, Providence, RI 02912 (United States); Karniadakis, George Em [Division of Applied Mathematics, Brown University, 182 George Street, Box F, Providence, RI 02912 (United States)]. E-mail: gk@dam.brown.edu
2006-04-10T23:59:59.000Z
We present numerical solutions of the stochastic Korteweg-de Vries equation for three cases corresponding to additive time-dependent noise, multiplicative space-dependent noise and a combination of the two. We employ polynomial chaos for discretization in random space, and discontinuous Galerkin and finite difference for discretization in physical space. The accuracy of the stochastic solutions is investigated by comparing the first two moments against analytical and Monte Carlo simulation results. Of particular interest is the interplay of spatial discretization error with the stochastic approximation error, which is examined for different orders of spatial and stochastic approximation.
Numerical evaluation of propeller noise, including non-linear effects
White, Terence Alan
1985-01-01T23:59:59.000Z
of Hanson and Fink CHAPTER III. PREDICTION OF TRANSONIC FLOW FIELD 19 19 21 25 A. Pr opeller Geometr y B. Calculation of the Flow Field 25 27 1. Establishment of the Computational Domain. 27 2. Numerical Method 3. Boundary Conditions C... are given by K = (1-M sin u ) - M cosv 2 . 2 1/2 CD r (2. 23) M sinu Y= r ~&& M, )&/2 (2. 24) The method of Jou, unlike the other methods presented here, has not been used for propeller noise calculations. 7. The Method of Hanson and Fink 24 All...
Connecting Numerical Relativity and Data Analysis of Gravitational Wave Detectors
Deirdre Shoemaker; Karan Jani; Lionel London; Larne Pekowsky
2015-03-09T23:59:59.000Z
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.
Numerical Modeling At Coso Geothermal Area (2000) | 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 Home Page on DeliciousPlasmaP a gHigh Plains WindInformationNVN-079666New look toNotes from the CallNuclearNumerical Modeling At
Numerical Modeling At Coso Geothermal Area (2006) | 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 Home Page on DeliciousPlasmaP a gHigh Plains WindInformationNVN-079666New look toNotes from the CallNuclearNumerical Modeling
Numerical Modeling At Coso Geothermal Area (2007) | 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 Home Page on DeliciousPlasmaP a gHigh Plains WindInformationNVN-079666New look toNotes from the CallNuclearNumerical
Theory comparison and numerical benchmarking on neoclassical toroidal viscosity torque
Wang, Zhirui; Park, Jong-Kyu; Logan, Nikolas; Kim, Kimin; Menard, Jonathan E. [Princeton Plasma Physics Laboratory, Princeton, New Jersey 08543 (United States)] [Princeton Plasma Physics Laboratory, Princeton, New Jersey 08543 (United States); Liu, Yueqiang [Euratom/CCFE Association, Culham Science Centre, Abingdon OX14 3DB (United Kingdom)] [Euratom/CCFE Association, Culham Science Centre, Abingdon OX14 3DB (United Kingdom)
2014-04-15T23:59:59.000Z
Systematic comparison and numerical benchmarking have been successfully carried out among three different approaches of neoclassical toroidal viscosity (NTV) theory and the corresponding codes: IPEC-PENT is developed based on the combined NTV theory but without geometric simplifications [Park et al., Phys. Rev. Lett. 102, 065002 (2009)]; MARS-Q includes smoothly connected NTV formula [Shaing et al., Nucl. Fusion 50, 025022 (2010)] based on Shaing's analytic formulation in various collisionality regimes; MARS-K, originally computing the drift kinetic energy, is upgraded to compute the NTV torque based on the equivalence between drift kinetic energy and NTV torque [J.-K. Park, Phys. Plasma 18, 110702 (2011)]. The derivation and numerical results both indicate that the imaginary part of drift kinetic energy computed by MARS-K is equivalent to the NTV torque in IPEC-PENT. In the benchmark of precession resonance between MARS-Q and MARS-K/IPEC-PENT, the agreement and correlation between the connected NTV formula and the combined NTV theory in different collisionality regimes are shown for the first time. Additionally, both IPEC-PENT and MARS-K indicate the importance of the bounce harmonic resonance which can greatly enhance the NTV torque when E×B drift frequency reaches the bounce resonance condition.
Numerical treatment of interfaces for second-order wave equations
F. Parisi; M. Cécere; M. Iriondo; O. Reula
2014-06-12T23:59:59.000Z
In this article we develop a numerical scheme to deal with interfaces between touching numerical grids when solving the second-order wave equation. We show that it is possible to implement an interface scheme of "penalty" type for the second-order wave equation, similar to the ones used for first-order hyperbolic and parabolic equations, and the second-order scheme used by Mattsson et al (2008). These schemes, known as SAT schemes for finite difference approximations and penalties for spectral ones, and ours share similar properties but in our case one needs to pass at the interface a smaller amount of data than previously known schemes. This is important for multi-block parallelizations in several dimensions, for it implies that one obtains the same solution quality while sharing among different computational grids only a fraction of the data one would need for a comparable (in accuracy) SAT or Mattsson et al.'s scheme. The semi-discrete approximation used here preserves the norm and uses standard finite-difference operators satisfying summation by parts. For the time integrator we use a semi-implicit IMEX Runge-Kutta method. This is crucial, since the explicit Runge-Kutta method would be impractical given the severe restrictions that arise from the stiff parts of the equations.
New insights on numerical error in symplectic integration
Hugo Jiménez-Pérez; Jean-Pierre Vilotte; Barbara Romanowicz
2015-08-13T23:59:59.000Z
We implement and investigate the numerical properties of a new family of integrators connecting both variants of the symplectic Euler schemes, and including an alternative to the classical symplectic mid-point scheme, with some additional terms. This family is derived from a new method, introduced in a previous study, for generating symplectic integrators based on the concept of special symplectic manifold. The use of symplectic rotations and a particular type of projection keeps the whole procedure within the symplectic framework. We show that it is possible to define a set of parameters that control the additional terms providing a way of "tuning" these new symplectic schemes. We test the "tuned" symplectic integrators with the perturbed pendulum and we compare its behavior with an explicit scheme for perturbed systems. Remarkably, for the given examples, the error in the energy integral can be reduced considerably. There is a natural geometrical explanation, sketched at the end of this paper. This is the subject of a parallel article where a finer analysis is performed. Numerical results obtained in this paper open a new point of view on symplectic integrators and Hamiltonian error.
Fast numerical treatment of nonlinear wave equations by spectral methods
Skjaeraasen, Olaf [ProsTek, Institute for Energy Technology, P.O. Box 40, N-2027 Kjeller (Norway); Robinson, P. A. [School of Physics, University of Sydney, New South Wales 2006 (Australia); Newman, D. L. [Center for Integrated Plasma Studies, University of Colorado at Boulder, Boulder, Colorado 80309 (United States)
2011-02-15T23:59:59.000Z
A method is presented that accelerates spectral methods for numerical solution of a broad class of nonlinear partial differential wave equations that are first order in time and that arise in plasma wave theory. The approach involves exact analytical treatment of the linear part of the wave evolution including growth and damping as well as dispersion. After introducing the method for general scalar and vector equations, we discuss and illustrate it in more detail in the context of the coupling of high- and low-frequency plasma wave modes, as modeled by the electrostatic and electromagnetic Zakharov equations in multiple dimensions. For computational efficiency, the method uses eigenvector decomposition, which is particularly advantageous when the wave damping is mode-dependent and anisotropic in wavenumber space. In this context, it is shown that the method can significantly speed up numerical integration relative to standard spectral or finite difference methods by allowing much longer time steps, especially in the limit in which the nonlinear Schroedinger equation applies.
Numerical Model for Conduction-Cooled Current Lead Heat Loads
White, M.J.; Wang, X.L.; /Fermilab; Brueck, H.D.; /DESY
2011-06-10T23:59:59.000Z
Current leads are utilized to deliver electrical power from a room temperature junction mounted on the vacuum vessel to a superconducting magnet located within the vacuum space of a cryostat. There are many types of current leads used at laboratories throughout the world; however, conduction-cooled current leads are often chosen for their simplicity and reliability. Conduction-cooled leads have the advantage of using common materials, have no superconducting/normal state transition, and have no boil-off vapor to collect. This paper presents a numerical model for conduction-cooled current lead heat loads. This model takes into account varying material and fluid thermal properties, varying thicknesses along the length of the lead, heat transfer in the circumferential and longitudinal directions, electrical power dissipation, and the effect of thermal intercepts. The model is validated by comparing the numerical model results to ideal cases where analytical equations are valid. In addition, the XFEL (X-Ray Free Electron Laser) prototype current leads are modeled and compared to the experimental results from testing at DESY's XFEL Magnet Test Stand (XMTS) and Cryomodule Test Bench (CMTB).
MAGNETIZATION DEGREE OF GAMMA-RAY BURST FIREBALLS: NUMERICAL STUDY
Harrison, Richard; Kobayashi, Shiho, E-mail: R.M.Harrison@2006.ljmu.ac.uk [Astrophysics Research Institute, Liverpool John Moores University, Liverpool, L3 5RF (United Kingdom)
2013-08-01T23:59:59.000Z
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.
NUMERICAL ANALYSIS OF A FINITE ELEMENT SCHEME FOR THE APPROXIMATION OF HARMONIC MAPS INTO SURFACES
Bartels, Soeren
NUMERICAL ANALYSIS OF A FINITE ELEMENT SCHEME FOR THE APPROXIMATION OF HARMONIC MAPS INTO SURFACES Geometric partial differential equations and their analysis as well as numerical simulation have recently
Mathiesen, Patrick; Kleissl, Jan
2011-01-01T23:59:59.000Z
transport and numerical weather modeling. J. Applied cross correlations. Weather and Forecasting, 8:4, 401?of radiation for numerical weather prediction and climate
Direct Numerical Simulation of Autoignition in a Jet in a Cross...
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Direct Numerical Simulation of Autoignition in a Jet in a Cross-Flow Direct Numerical Simulation of Autoignition in a Jet in a Cross-Flow PI Name: Christos Frouzakis PI Email:...
Farquharson, Colin G.
in Numerical Analysis, Scientific Computing or Computational and Applied Geophysics Applications are invited University of Newfoundland Centre for Numerical Analysis and Scientific Computing, within the DepartmentMemorial University, Department of Mathematics and Statistics Postdoctoral Research Position
Numerical simulation of oxidation and its effect on the crack growth resistance of titanium alloys
Triharjanto, Robertus Heru
1997-01-01T23:59:59.000Z
of titanium alloys used as metal matrix in MMC is imperative. This research effort develops a numerical modeling procedure to simulate the oxidation propagating from the crack surfaces and evaluate the effect on the mechanical behavior. The numerical modeling...
A numerical procedure to model and monitor CO2 sequestration in
Santos, Juan
Oceanografia e di Geofisica Sperimentale (OGS), Trieste, ITALY IC-MSQUARE 2012, September 3-7, 2012A numerical
Pearce, Fred
2003-01-01T23:59:59.000Z
We use a 3-D finite difference numerical model to generate synthetic seismograms from a simple fractured reservoir
THESIS FOR THE DEGREE OF LICENTIATE OF ENGINEERING Numerical Folding of Airbags
Patriksson, Michael
THESIS FOR THE DEGREE OF LICENTIATE OF ENGINEERING Numerical Folding of Airbags Based of Technology and GÂ¨oteborg University GÂ¨oteborg, Sweden 2007 #12;Numerical Folding of Airbags BasedÂ¨oteborg, Sweden 2007 #12;Numerical Folding of Airbags Based on Optimization and Origami Christoffer Cromvik
Numerical Methods for the Bogoliubov-Tolmachev-Shirkov model in superconductivity theory
Zhihao Ge; Ruihua Li
2015-03-08T23:59:59.000Z
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.
Numerical exact controllability of the 1D heat equation: Carleman weights and duality
Numerical exact controllability of the 1D heat equation: Carleman weights and duality Enrique Fern for the 1D heat equation. The goal is to compute a control that drives (a numerical approximation of [11], where primal methods were considered. Keywords: Heat equation, null controllability, numerical
Developing numerical libraries in Java RONALD F. BOISVERT1, JACK J. DONGARRA2, ROLDAN POZO1,
Dongarra, Jack
A. REMINGTON1 AND G.W. STEWART1;3 1Mathematical and Computational Sciences Division, Information. Finally, we outline the Java Numerical Toolkit1 JNT, which is meant to provide a base of computational;Developing numerical libraries in Java 2 2. NUMERICAL COMPUTING IN JAVA Java is both a computer language
The NumLab Numerical Laboratory J. Maubach 1 and A. Telea 2
Eindhoven, Technische Universiteit
The NumLab Numerical Laboratory J. Maubach 1 and A. Telea 2 Department of Mathematics and Computer environments addresses numerical simulation, interac- tive visualisation and computational steering. Most- bling numerical simulations from computational and visualisation blocks, as well as building such blocks
An Interactive Course in Numerical Methods for the Earth Susan E. Allen
Stockie, John
are typically in physics, mathematics, chemistry, biology or geography, and exposure to numerical computing science are becoming more and more dependent on numerical simulation, modelling, and computationally courses in numerical methods using MatlabTM as the programming language, or a graduate-level mathematics
A LOCAL DIMENSION TEST FOR NUMERICALLY APPROXIMATED POINTS ON ALGEBRAIC SETS
Sommese, Andrew J.
of magnitude improvement in the computation of a numerical irreducible decomposition. Also, to compute the irreducible components of a fixed dimension, it is no longer necessary to compute the numerical irreducible components that contain a specified point p; 4. computing the numerical irreducible decomposition of V more e
Mat 684: Numerical Analysis II, Spring 2013 Instructor: Dr. Yuesheng Xu
McConnell, Terry
:30-11:30, Monday and Wednesday Text book: Numerical Analysis: Mathematics of Scientific Computing, by David KincaidMat 684: Numerical Analysis II, Spring 2013 Instructor: Dr. Yuesheng Xu Office: 206 E Carnegie Hall and Ward Cheney, 3rd Edition, Brooks/Cole. Credit: 3 hours. Topics Covered by the Course: · Numerical
Perturbing Numerical Calculations for Statistical Analysis of Floating-Point Program (In)Stability
Su, Zhendong
mathematical techniques for performing error and stability analysis of numerical algorithms. However to numerical computing tend to suffer from misconceptions that arise from thinking in terms of R, not itsPerturbing Numerical Calculations for Statistical Analysis of Floating-Point Program (In
NUMERICAL SOLUTION OF HYBRID OPTIMAL CONTROL PROBLEMS WITH APPLICATIONS IN ROBOTICS
Stryk, Oskar von
NUMERICAL SOLUTION OF HYBRID OPTIMAL CONTROL PROBLEMS WITH APPLICATIONS IN ROBOTICS Martin Buss.g. (Buss et al., 2000; Buss, 2001; Tomlin, 1999). The key to numerically solving HOCPs appears to lie (Buss et al., 2000; Stryk and Glocker, 2000). This paper presents numerical solution techniques
Numerical study of high heat ux pool boiling heat transfer Ying He a,*, Masahiro Shoji b
Maruyama, Shigeo
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
Field Survey and Numerical Simulations: A Review of the 1998 Papua New Guinea Tsunami
Lynett, Patrick
Field Survey and Numerical Simulations: A Review of the 1998 Papua New Guinea Tsunami PATRICK J (PNG) tsunami of 1998 is re-examined through a detailed review of the field survey as well as numerous of frequency dispersion on the landslide-generated tsunami. The numerical comparisons indicate that the NLSW
Regier, Terry
Numeral systems across languages support efficient communication: From approximate numerosity this spectrum. We find that these numeral systems all reflect a functional need for efficient communication differ qualitatively in their numeral systems. At one extreme, some languages have a small set of number
Mathematical and Numerical Analyses of Peridynamics for Multiscale Materials Modeling
Du, Qiang
2014-11-12T23:59:59.000Z
The rational design of materials, the development of accurate and efficient material simulation algorithms, and the determination of the response of materials to environments and loads occurring in practice all require an understanding of mechanics at disparate spatial and temporal scales. The project addresses mathematical and numerical analyses for material problems for which relevant scales range from those usually treated by molecular dynamics all the way up to those most often treated by classical elasticity. The prevalent approach towards developing a multiscale material model couples two or more well known models, e.g., molecular dynamics and classical elasticity, each of which is useful at a different scale, creating a multiscale multi-model. However, the challenges behind such a coupling are formidable and largely arise because the atomistic and continuum models employ nonlocal and local models of force, respectively. The project focuses on a multiscale analysis of the peridynamics materials model. Peridynamics can be used as a transition between molecular dynamics and classical elasticity so that the difficulties encountered when directly coupling those two models are mitigated. In addition, in some situations, peridynamics can be used all by itself as a material model that accurately and efficiently captures the behavior of materials over a wide range of spatial and temporal scales. Peridynamics is well suited to these purposes because it employs a nonlocal model of force, analogous to that of molecular dynamics; furthermore, at sufficiently large length scales and assuming smooth deformation, peridynamics can be approximated by classical elasticity. The project will extend the emerging mathematical and numerical analysis of peridynamics. One goal is to develop a peridynamics-enabled multiscale multi-model that potentially provides a new and more extensive mathematical basis for coupling classical elasticity and molecular dynamics, thus enabling next generation atomistic-to-continuum multiscale simulations. In addition, a rigorous study of nite element discretizations of peridynamics will be considered. Using the fact that peridynamics is spatially derivative free, we will also characterize the space of admissible peridynamic solutions and carry out systematic analyses of the models, in particular rigorously showing how peridynamics encompasses fracture and other failure phenomena. Additional aspects of the project include the mathematical and numerical analysis of peridynamics applied to stochastic peridynamics models. In summary, the project will make feasible mathematically consistent multiscale models for the analysis and design of advanced materials.
Numerical Study of Induced False Vacuum Decay at High Energies
A. N. Kuznetsov; P. G. Tinyakov
1995-11-02T23:59:59.000Z
We calculate numerically the probability $\\exp[ {1\\over\\lambda} F(E/E_{sph},N/N_{sph})]$ of the false vacuum decay in the massive four-dimensional $-\\lambda\\phi^4$ model from multiparticle initial states with fixed number of particles $N$ and energy $E$ greater than the height of the barrier $E_{sph}$. We find that at $E\\lsim 3E_{sph}$ and $N\\lsim 0.4N_{sph}$ the decay is classically forbidden and thus is exponentially suppressed. We argue that the classically forbidden region extends at small $N$ at least up to $E\\sim 10 E_{sph}$ and, most likely, to all energies. Our data suggest that the false vacuum decay induced by two-particle collisions is exponentially suppressed at all energies.
Numerical determination of the material properties of porous dust cakes
D. Paszun; C. Dominik
2008-02-13T23:59:59.000Z
The formation of planetesimals requires the growth of dust particles through collisions. Micron-sized particles must grow by many orders of magnitude in mass. In order to understand and model the processes during this growth, the mechanical properties, and the interaction cross sections of aggregates with surrounding gas must be well understood. Recent advances in experimental (laboratory) studies now provide the background for pushing numerical aggregate models onto a new level. We present the calibration of a previously tested model of aggregate dynamics. We use plastic deformation of surface asperities as the physical model to bring critical velocities for sticking into accordance with experimental results. The modified code is then used to compute compression strength and the velocity of sound in the aggregate at different densities. We compare these predictions with experimental results and conclude that the new code is capable of studying the properties of small aggregates.
MHD Remote Numerical Simulations: Evolution of Coronal Mass Ejections
L. Hernandez-Cervantes; A. Santillan; A. R. Gonzalez-Ponce
2008-12-22T23:59:59.000Z
Coronal mass ejections (CMEs) are solar eruptions into interplanetary space of as much as a few billion tons of plasma, with embedded magnetic fields from the Sun's corona. These perturbations play a very important role in solar--terrestrial relations, in particular in the spaceweather. In this work we present some preliminary results of the software development at the Universidad Nacional Autonoma de Mexico to perform Remote MHD Numerical Simulations. This is done to study the evolution of the CMEs in the interplanetary medium through a Web-based interface and the results are store into a database. The new astrophysical computational tool is called the Mexican Virtual Solar Observatory (MVSO) and is aimed to create theoretical models that may be helpful in the interpretation of observational solar data.
Numerical Weather Forecasting at the Savannah River Site
Buckley, R.L. [Westinghouse Savannah River Company, AIKEN, SC (United States)
1998-11-01T23:59:59.000Z
Weather forecasts at the Savannah River Site (SRS) are important for applications to emergency response. The fate of accidentally-released radiological materials and toxic chemicals can be determined by providing wind and turbulence input to atmospheric transport models. This operation has been routinely performed at SRS using the WIND System, a system of computer models and monitors which collect data from towers situated throughout the SRS. However, the information provided to these models is spatially homogeneous (in one or two dimensions) with an elementary forecasting capability. This paper discusses the use of an advanced three-dimensional prognostic numerical model to provide space and time-dependent meteorological data for use in the WIND System dispersion models. The extensive meteorological data collection at SRS serves as a ground truth for further model development as well as for use in other applications.
Laser spark plug numerical design process with experimental validation
McIntyre, D.; Woodruff, S.
2011-01-01T23:59:59.000Z
This work reports the numerical modeling design procedure for a miniaturized laser spark plug. In previous work both side pumped and end pumped laser spark plugs were empirically designed and tested. Experimental data from the previous laser spark plug development cycles is compared to the output predicted by a known set of rate equations. The rate equations are used to develop interrelated inter cavity time dependent waveforms that are then used to identify key variables. These variables are then input to a set of secondary equations for determining the output pulse energy, output power, and output pulse width of the simulated laser system. The physical meaning and the operation of the rate equations is explained in detail. This paper concentrates on the process and decision points needed to successfully design a solid state passively Q-switched laser system, either side pumped or end pumped, that produces the appropriate output needed for use as a laser spark plug for internal combustion engines.
Direct numerical simulation of pattern formation in subaqueous sediment
Kidanemariam, Aman G
2014-01-01T23:59:59.000Z
We present results of direct numerical simulation of incompressible fluid flow over a thick bed of mobile, spherically-shaped particles. The algorithm is based upon the immersed boundary technique for fluid-solid coupling and uses a soft-sphere model for the solid-solid contact. Two parameter points in the laminar flow regime are chosen, leading to the emergence of sediment patterns classified as `small dunes', while one case under turbulent flow conditions leads to `vortex dunes' with significant flow separation on the lee side. Wavelength, amplitude and propagation speed of the patterns extracted from the spanwise-averaged fluid-bed interface are found to be consistent with available experimental data. The particle transport rates are well represented by available empirical models for flow over a plane sediment bed in both the laminar and the turbulent regimes.
Numerical Study of Shock-Cylinder Banks Interactions
Wang, S.P.; Anderson, M.H.; Oakley, J.G.; Bonazza, R. [University of Wisconsin-Madison (United States)
2003-07-15T23:59:59.000Z
A numerical parametric study of shock-cylinder banks interactions is presented using a high resolution Euler solver. Staggered cylinder banks of five rows are chosen with the purpose of modeling IFE reactor cooling tube banks. The effect of the aspect ratio of the intercylinder pitch to the distance between successive cylinder rows on the vertical pressure forces acting on the cylinders with different geometries is investigated. Preliminary results show that the largest vertical force develops on the cylinders of the second or third row. This peak pressure force increases with decreasing values of the aspect ratio. It is shown that an increasing second force peak also appears on the successive rows, starting with the second one, with decreasing aspect ratio. It is also observed that the force on the last-row cylinders basically decreases to the level of that on the first row. The results are useful for the optimal design of the cooling tubes system of IFE reactors.
Numerical analysis of kinematic soil-pile interaction
Castelli, Francesco; Maugeri, Michele [Department of Civil and Environmental Engineering, University of Catania, Viale Andrea Doria no. 6, 95125, Catania (Italy); Mylonakis, George [Department of Civil Engineering, University of Patras, Rio GR-26500, Patras (Greece)
2008-07-08T23:59:59.000Z
In the present study, the response of singles pile to kinematic seismic loading is investigated using the computer program SAP2000. The objectives of the study are: (1) to develop a numerical model that can realistically simulate kinematic soil-structure interaction for piles accounting for discontinuity conditions at the pile-soil interface, energy dissipation and wave propagation; (2) to use the model for evaluating kinematic interaction effects on pile response as function of input ground motion; and (3) to present a case study in which theoretical predictions are compared with results obtained from other formulations. To evaluate the effects of kinematic loading, the responses of both the free-field soil (with no piles) and the pile were compared. Time history and static pushover analyses were conducted to estimate the displacement and kinematic pile bending under seismic loadings.
Numerical computations of the dynamics of fluidic membranes and vesicles
Barrett, John W; Nürnberg, Robert
2015-01-01T23:59:59.000Z
Vesicles and many biological membranes are made of two monolayers of lipid molecules and form closed lipid bilayers. The dynamical behaviour of vesicles is very complex and a variety of forms and shapes appear. Lipid bilayers can be considered as a surface fluid and hence the governing equations for the evolution include the surface (Navier--)Stokes equations, which in particular take the membrane viscosity into account. The evolution is driven by forces stemming from the curvature elasticity of the membrane. In addition, the surface fluid equations are coupled to bulk (Navier--)Stokes equations. We introduce a parametric finite element method to solve this complex free boundary problem, and present the first three dimensional numerical computations based on the full (Navier--)Stokes system for several different scenarios. For example, the effects of the membrane viscosity, spontaneous curvature and area difference elasticity (ADE) are studied. In particular, it turns out, that even in the case of no viscosit...
A numerical model of aerosol scavenging: Part 1, Microphysics parameterization
Molenkamp, C.R.; Bradley, M.M.
1991-09-01T23:59:59.000Z
We have developed a three-dimensional numerical model (OCTET) to simulate the dynamics and microphysics of clouds and the transport, diffusion and precipitation scavenging of aerosol particles. In this paper we describe the cloud microphysics and scavenging parameterizations. The representation of cloud microphysics is a bulk- water parameterization which includes water vapor and five types of hydrometeors (cloud droplets, rain drops, ice crystals, snow, and graupel). A parallel parameterization represents the scavenging interactions between pollutant particles and hydrometeors including collection of particles because of condensation nucleation, Brownian and phoretic attachment, and inertial capture, resuspension because of evaporation and sublimation; and transfer interactions where particles collected by one type of hydrometeor are transferred to another type of freezing, melting, accretion, riming and autoconversion.
Electromagnetic Pulse Propagation over Nonuniform Earth Surface: Numerical Simulation
Alexei V. Popov; Vladimir V. Kopeikin
2007-04-14T23:59:59.000Z
We simulate EM pulse propagation along the nonuniform earth surface using so called time-domain parabolic equation. To solve it by finite differences, we introduce a time-domain analog of the impedance boundary condition and a nonlocal BC of transparency reducing open computational domain to a strip of finite width. Numerical examples demonstrate influence of soil conductivity on the wide-band pulse waveform. For a high-frequency modulated EM pulse, we develop an asymptotic approach based on the ray structure of the monochromatic wave field at carrier frequency. This radically diminishes the computation costs and allows for pulsed wave field calculation in vast domains measured by tens of thousands wavelengths.
Spherically symmetric cosmological spacetimes with dust and radiation — numerical implementation
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: wclim@waikato.ac.nz, E-mail: regis@to.infn.it, E-mail: chris.clarkson@gmail.com [Astrophysics, Cosmology and Gravity Centre, and Department of Mathematics and Applied Mathematics, University of Cape Town, Rondebosch 7701, Cape Town (South Africa)
2013-10-01T23:59:59.000Z
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.
Numerical integration of thermal noise in relativistic hydrodynamics
Clint Young
2013-06-03T23:59:59.000Z
Thermal fluctuations affect the dynamics of systems near critical points, the evolution of the early universe, and two-particle correlations in heavy-ion collisions. For the latter, numerical simulations of nearly-ideal, relativistic fluids are necessary. The correlation functions of noise in relativistic fluids are calculated, stochastic integration of the noise in 3+1-dimensional viscous hydrodynamics is implemented, and the effect of noise on observables in heavy-ion collisions is discussed. Thermal fluctuations will cause significant variance in the event-by-event distributions of integrated v2 while changing average values even when using the same initial conditions, suggesting that including thermal noise will lead to refitting of the hydrodynamical parameters with implications for understanding the physics of hot QCD.
A fast direct numerical simulation method for characterising hydraulic roughness
Chung, Daniel; MacDonald, Michael; Hutchins, Nicholas; Ooi, Andrew
2015-01-01T23:59:59.000Z
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...
Numerical simulation of laminar reacting flows with complex chemistry
Day, Marcus S.; Bell, John B.
1999-12-01T23:59:59.000Z
We present an adaptive algorithm for low Mach number reacting flows with complex chemistry. Our approach uses a form of the low Mach number equations that discretely conserves both mass and energy. The discretization methodology is based on a robust projection formulation that accommodates large density contrasts. The algorithm uses an operator-split treatment of stiff reaction terms and includes effects of differential diffusion. The basic computational approach is embedded in an adaptive projection framework that uses structured hierarchical grids with subcycling in time that preserves the discrete conservation properties of the underlying single-grid algorithm. We present numerical examples illustrating the performance of the method on both premixed and non-premixed flames.
Integrating Numerical Computation into the Modeling Instruction Curriculum
Caballero, Marcos D; Aiken, John M; Douglas, Scott S; Scanlon, Erin M; Thoms, Brian; Schatz, Michael F
2012-01-01T23:59:59.000Z
We describe a way to introduce physics high school students with no background in programming to computational problem-solving experiences. Our approach builds on the great strides made by the Modeling Instruction reform curriculum. This approach emphasizes the practices of "Developing and using models" and "Computational thinking" highlighted by the NRC K-12 science standards framework. We taught 9th-grade students in a Modeling-Instruction-based physics course to construct computational models using the VPython programming environment. Numerical computation within the Modeling Instruction curriculum provides coherence among the curriculum's different force and motion models, links the various representations which the curriculum employs, and extends the curriculum to include real-world problems that are inaccessible to a purely analytic approach.
Numerical investigation of the threshold for primordial black hole formation
J. C. Niemeyer
1998-06-02T23:59:59.000Z
First results of a numerical investigation of primordial black hole formation in the radiation dominated phase of the Early Universe are presented. The simulations follow the gravitational collapse of three different families of high-amplitude density fluctuations imposed at the time of horizon crossing. The threshold for black hole formation, \\delta_{c} \\approx 0.7, is found to be nearly identical for all perturbation families if the control parameter, \\delta, is chosen as the total excess mass within the initial horizon volume. Furthermore, we demonstrate that the scaling of black hole mass with distance from the formation threshold, known to occur in near-critical gravitational collapse, applies to primordial black hole formation.
Explicit Gravitational Radiation in Hyperbolic Systems for Numerical Relativity
C. Bona; C. Palenzuela
2002-02-28T23:59:59.000Z
A method for studying the causal structure of space-time evolution systems is presented. This method, based on a generalization of the well known Riemann problem, provides intrinsic results which can be interpreted from the geometrical point of view. A one-parameter family of hyperbolic evolution systems is presented and the physical relevance of their characteristic speeds and eigenfields is discussed. The two degrees of freedom corresponding to gravitational radiation are identified in an intrinsic way, independent of the space coordinate system. A covariant interpretation of these degrees of freedom is provided in terms of the geometry of the wave fronts. The requirement of a consistent geometrical interpretation of the gravitational radiation degrees of freedom is used to solve the ordering ambiguity that arises when obtaining first order evolution systems from the second order field equations. This achievement provides a benchmark which can be used to check both the existing and future first order hyperbolic formalisms for Numerical Relativity.
Flynn, Morris R.
MecE 390 Numerical Methods (Fall 2014) Objectives: To provide a basic familiarity or by appoint- ment. Text: Numerical Methods for Engineers, Sixth edition by Steven C. Chapra and Raymond P of MEC. (For more information, please consult the detailed course schedule on the course homepage
Seismic Structure And Seismicity Of The Cooling Lava Lake Of...
results from various active and passive experiments conducted on Kilauea Iki in March 1976 suggest a rather thin (less than 10 m) zone of residual melt, the lateral extent of...
Harnessing Sun, Wind and Lava for Islands' Energy Needs | Department...
Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site
Nations (EDIN) project - this international collaboration between the United States, Iceland and New Zealand is aimed at helping islands adopt clean energy policies, technology...
Monitoring Erebus volcano's active lava lake; tools, techniques and observations
Peters, Nial John
2015-03-03T23:59:59.000Z
are those that persistently expose magma at the Earth’s surface. In this respect, they encompass many of the popular stereo- types of volcanic activity, exhibiting sustained passive degassing, which is often ac- companied sporadically by varying degrees... 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...
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 Center Home Page on DeliciousPlasmaP aCentrothermDepew,IndependentCounty,KentPlatformLamoilleLauderdaleGeorgia: Energy ResourcesPageEnergy
Testing numerical relativity with the shifted gauge wave
Maria C. Babiuc; Bela Szilagyi; Jeffrey Winicour
2006-02-17T23:59:59.000Z
Computational methods are essential to provide waveforms from coalescing black holes, which are expected to produce strong signals for the gravitational wave observatories being developed. Although partial simulations of the coalescence have been reported, scientifically useful waveforms have so far not been delivered. The goal of the AppleswithApples (AwA) Alliance is to design, coordinate and document standardized code tests for comparing numerical relativity codes. The first round of AwA tests have now being completed and the results are being analyzed. These initial tests are based upon periodic boundary conditions designed to isolate performance of the main evolution code. Here we describe and carry out an additional test with periodic boundary conditions which deals with an essential feature of the black hole excision problem, namely a non-vanishing shift. The test is a shifted version of the existing AwA gauge wave test. We show how a shift introduces an exponentially growing instability which violates the constraints of a standard harmonic formulation of Einstein's equations. We analyze the Cauchy problem in a harmonic gauge and discuss particular options for suppressing instabilities in the gauge wave tests. We implement these techniques in a finite difference evolution algorithm and present test results. Although our application here is limited to a model problem, the techniques should benefit the simulation of black holes using harmonic evolution codes.
Probing Modified Gravity with Atom-Interferometry: a Numerical Approach
Schlogel, Sandrine; Fuzfa, Andre
2015-01-01T23:59:59.000Z
Refined constraints on chameleon theories are calculated for atom-interferometry experiments, using a numerical approach consisting in solving for a four-region model the static and spherically symmetric Klein-Gordon equation for the chameleon field. By modeling not only the test mass and the vacuum chamber but also its walls and the exterior environment, the method allows to probe new effects on the scalar field profile and the induced acceleration of atoms. In the case of a weakly perturbing test mass, the effect of the wall is to enhance the field profile and to lower the acceleration inside the chamber by up to one order of magnitude. In the thin-shell regime, significant deviations from the analytical estimations are found, even when measurements are realized in the immediate vicinity of the test mass. Close to the vacuum chamber wall, the acceleration becomes negative and potentially measurable. This prediction could be used to discriminate between fifth-force effects and systematic experimental uncerta...
Effective viscosity of active suspensions: Three-dimensional numerical modeling
Levan Jibuti; Walter Zimmermann; Salima Rafaï; Philippe Peyla
2014-12-10T23:59:59.000Z
A three-dimensional model is proposed for Chlamydomonas Reinhardtii swimming with a breaststroke-like beating of its two flagella. The model reveals unusual angular orbits of the active swimmer under a linear shear flow. Namely, the swimmer sustains orientation transiently across the flow when flagella plane is perpendicular to the shear plane, and amplify the shear-induced rotation along the flow. Such behavior is a result of the interplay between shear-induced deformation and swimmer's periodic beating motion that exerts internal torques on the torque-free swimmer. This particular behavior has some significant consequences on the rheological properties of the suspension that tends to confirm previous experimental results [Phys. Rev. Lett. 104, 098102 (2010)]. We calculated the intrinsic viscosity of the suspension with such isolated modeled microswimmers (dilute case) in shear flow using numerical simulations based on Rotne-Prager approximation. The results show an increased intrinsic viscosity for active swimmer suspensions in comparison to non-active ones in accordance with previous experimental measurements. A major enhancement of the active swimmer viscosity occurs due to the effectively extended shape of the deformable swimming cells. We also recover the experimentally observed shear thinning behavior.
Numerical study of energy diffusion in King models
Tom Theuns
1995-11-07T23:59:59.000Z
The energy diffusion coefficients D_n(E) (n=1,2) for a system of equal mass particles moving self-consistently in an N-body realisation of a King model are computed from the probability per unit time, P(E, Delta E), that a star with initial energy E will undergo an energy change Delta E. In turn, P is computed from the number of times during the simulation that a particle in a state of given energy undergoes a transition to another state. These particle states are defined directly from the time evolution of E by identifying them with the event occuring between two local maxima in the E(t) curve. If one assumes next that energy changes are uncorrelated between different states, one can use diffusion theory to compute D_n(E). The simulations employ N=512, 2048,... , 32768 particles and are performed using an implementation of Aarseth's direct integrator N-body1 on a massively parallel computer. The more than seven million transitions measured in the largest N simulation provide excellent statistics. The numerically determined D(E)'s are compared against their theoretical counterparts which are computed from phase-space averaged rates of energy change due to independent binary encounters. The overall agreement between them is impressive over most of the energy range, notwithstanding the very different type of approximations involved, giving considerable support to the valid usage of these theoretical expressions to simulate dynamical evolution in Fokker-Planck type calculations.
Numerical weather forecasting at the Savannah River Site
Buckley, R.L. [Westinghouse Savannah River Site, Aiken, SC (United States)
1998-12-31T23:59:59.000Z
Weather forecasts at the Savannah River Site (SRS) are important for applications to emergency response. The fate of accidentally released radiological materials and toxic chemicals can be determined by providing wind and turbulence input to atmospheric transport models. This operation has been routinely performed at SRS using the WIND system, a system of computer models and monitors that collects data from towers situated throughout the SRS. However, the information provided to these models is spatially homogeneous (in one or two dimensions) with an elementary forecasting capability. This paper discusses the use of an advanced three-dimensional prognostic numerical model to provide space- and time-dependent meteorological data for use in the WIND system dispersion models. The extensive meteorological data collection at SRS serves as a ground truth for further model development as well as for use in other applications. A prognostic mesoscale model, the regional atmospheric modeling system (RAMS), is used to provide these forecasts. Use of RAMS allows for incorporation of mesoscale features such as the sea breeze, which has been shown to affect local weather conditions. This paper discusses the mesoscale model and its configuration for the operational simulation, as well as an application using a dispersion model at the SRS.
Advanced Numerical Methods and Software Approaches for Semiconductor Device Simulation
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Carey, Graham F.; Pardhanani, A. L.; Bova, S. W.
2000-01-01T23:59:59.000Z
In this article we concisely present several modern strategies that are applicable to driftdominated carrier transport in higher-order deterministic models such as the driftdiffusion, hydrodynamic, and quantum hydrodynamic systems. The approaches include extensions of “upwind” and artificial dissipation schemes, generalization of the traditional Scharfetter – Gummel approach, Petrov – Galerkin and streamline-upwind Petrov Galerkin (SUPG), “entropy” variables, transformations, least-squares mixed methods and other stabilized Galerkin schemes such as Galerkin least squares and discontinuous Galerkin schemes. The treatment is representative rather than an exhaustive review and several schemes are mentioned only briefly with appropriate reference to the literature. Some of themore »methods have been applied to the semiconductor device problem while others are still in the early stages of development for this class of applications. We have included numerical examples from our recent research tests with some of the methods. A second aspect of the work deals with algorithms that employ unstructured grids in conjunction with adaptive refinement strategies. The full benefits of such approaches have not yet been developed in this application area and we emphasize the need for further work on analysis, data structures and software to support adaptivity. Finally, we briefly consider some aspects of software frameworks. These include dial-an-operator approaches such as that used in the industrial simulator PROPHET, and object-oriented software support such as those in the SANDIA National Laboratory framework SIERRA.« less
Numeric spectral radiation hydrodynamic calculations of supernova shock breakouts
Sapir, Nir; Halbertal, Dorri [Department of Particle Physics and Astrophysics, Weizmann Institute of Science, Rehovot 76100 (Israel)
2014-12-01T23:59:59.000Z
We present here an efficient numerical scheme for solving the non-relativistic one-dimensional radiation-hydrodynamics equations including inelastic Compton scattering, which is not included in most codes and is crucial for solving problems such as shock breakout. The devised code is applied to the problems of a steady-state planar radiation mediated shock (RMS) and RMS breakout from a stellar envelope. The results are in agreement with those of a previous work on shock breakout, in which Compton equilibrium between matter and radiation was assumed and the 'effective photon' approximation was used to describe the radiation spectrum. In particular, we show that the luminosity and its temporal dependence, the peak temperature at breakout, and the universal shape of the spectral fluence derived in this earlier work are all accurate. Although there is a discrepancy between the spectral calculations and the effective photon approximation due to the inaccuracy of the effective photon approximation estimate of the effective photon production rate, which grows with lower densities and higher velocities, the difference in peak temperature reaches only 30% for the most discrepant cases of fast shocks in blue supergiants. The presented model is exemplified by calculations for supernova 1987A, showing the detailed evolution of the burst spectrum. The incompatibility of the stellar envelope shock breakout model results with observed properties of X-ray flashes (XRFs) and the discrepancy between the predicted and observed rates of XRFs remain unexplained.
Numerical models of pressure pulse generation by imploding metal liners
Humphries, S. Jr. [Acceleration Associates, Albuquerque, NM (United States)] [Acceleration Associates, Albuquerque, NM (United States); [Univ. of New Mexico, Albuquerque, NM (United States). Dept. of Electrical Engineering; Ekdahl, C.A. [Los Alamos National Lab., NM (United States)] [Los Alamos National Lab., NM (United States)
1996-12-01T23:59:59.000Z
The authors describe numerical calculations of pressure pulse generation using imploding liners. Liners are metal cylinders that are magnetically compressed by an intense axial current flow from a high-power pulse generator. The simulations cover the acceleration of the liner, collision with an internal diagnostic target, followed by compression and shock wave heating of the target. With the projected current waveform of the Atlas capacitor bank (in development at Los Alamos National Laboratory), initial results suggest that it may be possible to achieve pressures exceeding 3,000 Gpa (30 Mbar) in a 4 mm diameter sample over an interval of 100--200 ns. The simulations were carried out with Crunch, a new one-dimensional hydrodynamics package for advanced personal computers. The program uses finite-element techniques to solve the coupled problems of hydrodynamics and magnetic diffusion. Crunch fully supports loading and interpolating Sesame equation-of-state tables. The program exhibits excellent stability, even for collisions between material shells and shock convergence on axis. In contrast to previous work, the present studies follow the full process through solid target collision and compression. The work supports the High-Energy Density Physics Program of Los Alamos National Laboratory (LANL), a component of the US Department of Energy Stockpile Stewardship Program. The purpose of this program is maintenance of the nuclear weapons stockpile through improved computational ability and above-ground experiments. Imploding liners driven by conventional capacitor banks constitute a portion of the program to study matter at high pressure.
Numerical studies on the geomechanical stability ofhydrate-bearing sediments
Rutqvist, Jonny; Moridis, George J.
2007-05-01T23:59:59.000Z
The thermal and mechanical loading of oceanicHydrate-Bearing Sediments (HBS) can result in hydrate dissociation and asignificant pressure increase, with potentially adverse consequences onthe integrity and stability of the wellbore assembly, the HBS, and thebounding formations. The perception of HBS instability, coupled withinsufficient knowledge of their geomechanical behavior and the absence ofpredictive capabilities, have resulted in a strategy of avoidance of HBSwhen locating offshore production platforms, and can impede thedevelopment of hydrate deposits as gas resources.In this study weinvestigate in three cases of coupled hydraulic, thermodynamic andgeomechanical behavior of oceanic hydrate-bearing sediments. The firstinvolves hydrate heating as warm fluids from deeper conventionalreservoirs ascend to the ocean floor through uninsulated pipesintersecting the HBS. The second case describes system response duringgas production from a hydrate deposit, and the third involves mechanicalloading caused by the weight of structures placed on the ocean flooroverlying hydrate-bearing sediments.For the analysis of the geomechanicalstability of HBS, we developed and used a numerical model that integratesa commercial geomechanical code and a simulator describing the coupledprocesses of fluid flow, heat transport and thermodynamic behavior in theHBS. Our simulation results indicate that the stability of HBS in thevicinity of warm pipes may be significantly affected, especially if thesediments are unconsolidated and more compressible. Gas production fromoceanic deposits may also affect the geomechanical stability of HBS underthe conditions that are deemed desirablefor production. Conversely, theincreased pressure caused by the weight of structures on the ocean floorincreases the stability of underlying hydrates.
Numerical simulations of compressively driven interstellar turbulence: I. Isothermal gas
Schmidt, Wolfram; Hupp, Markus; Kern, Sebastian; Niemeyer, Jens C
2008-01-01T23:59:59.000Z
We performed numerical simulations of supersonic isothermal turbulence driven by mostly compressive large-scale forcing, using both a static grid and adaptive mesh refinement with an effective resolution N=768^3. After a transient phase dominated by shocks, turbulence evolves into a steady state with an RMS Mach number about 2.5, in which cloud-like structures of over-dense gas are surrounded by highly rarefied gas. The index of the turbulence energy spectrum function beta = 2.0 in the shock-dominated phase. As the flow approaches statistical equilibrium, the spectrum flattens, with beta = 1.9. For the scaling exponent of the root mean square velocity fluctuation, we obtain gamma = 0.43 from the velocity structure functions of second order. These results are well within the range of observed scaling properties for the velocity dispersion in molecular clouds. Calculating structure functions of order p=1,...,5, we find for all scaling exponents significant deviations from the Kolmogorov-Burgers model proposed b...
Numerical Relativity in Spherical Polar Coordinates: Off-center Simulations
Thomas W. Baumgarte; Pedro J. Montero; Ewald Müller
2015-06-03T23:59:59.000Z
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.
Numerical tests of dynamical friction in gravitational inhomogeneous systems
A. Del Popolo
2003-05-05T23:59:59.000Z
In this paper, I test by numerical simulations the results of Del Popolo & Gambera (1998),dealing with the extension of Chandrasekhar and von Neumann's analysis of the statistics of the gravitational field to systems in which particles (e.g., stars, galaxies) are inhomogeneously distributed. The paper is an extension of that of Ahmad & Cohen (1974), in which the authors tested some results of the stochastic theory of dynamical friction developed by Chandrasekhar & von Neumann (1943) in the case of homogeneous gravitational systems. It is also a continuation of the work developed in Del Popolo (1996a,b), which extended the results of Ahmad & Cohen (1973), (dealing with the study of the probability distribution of the stochastic force in homogeneous gravitational systems) to inhomogeneous gravitational systems. Similarly to what was done by Ahmad & Cohen (1974) in the case of homogeneous systems, I test, by means of the evolution of an inhomogeneous system of particles, that the theoretical rate of force fluctuation d F/dt describes correctly the experimental one, I find that the stochastic force distribution obtained for the evolved system is in good agreement with the Del Popolo & Gambera (1998) theory. Moreover, in an inhomogeneous background the friction force is actually enhanced relative to the homogeneous case.
Numerical simulations of super-luminous supernovae of type IIn
Dessart, Luc; Hillier, D John
2015-01-01T23:59:59.000Z
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...
DIPOLE COLLAPSE AND DYNAMO WAVES IN GLOBAL DIRECT NUMERICAL SIMULATIONS
Schrinner, Martin; Dormy, Emmanuel [MAG (ENS/IPGP), LRA, Ecole Normale Superieure, 24 Rue Lhomond, 75252 Paris Cedex 05 (France); Petitdemange, Ludovic, E-mail: martin@schrinner.eu [Previously at Max-Planck-Institut fuer Astronomie, Koenigstuhl 17, 69117 Heidelberg, Germany. (Germany)
2012-06-20T23:59:59.000Z
Magnetic fields of low-mass stars and planets are thought to originate from self-excited dynamo action in their convective interiors. Observations reveal a variety of field topologies ranging from large-scale, axial dipoles to more structured magnetic fields. In this article, we investigate more than 70 three-dimensional, self-consistent dynamo models in the Boussinesq approximation obtained by direct numerical simulations. The control parameters, the aspect ratio, and the mechanical boundary conditions have been varied to build up this sample of models. Both strongly dipolar and multipolar models have been obtained. We show that these dynamo regimes in general can be distinguished by the ratio of a typical convective length scale to the Rossby radius. Models with a predominantly dipolar magnetic field were obtained, if the convective length scale is at least an order of magnitude larger than the Rossby radius. Moreover, we highlight the role of the strong shear associated with the geostrophic zonal flow for models with stress-free boundary conditions. In this case the above transition disappears and is replaced by a region of bistability for which dipolar and multipolar dynamos coexist. We interpret our results in terms of dynamo eigenmodes using the so-called test-field method. We can thus show that models in the dipolar regime are characterized by an isolated 'single mode'. Competing overtones become significant as the boundary to multipolar dynamos is approached. We discuss how these findings relate to previous models and to observations.
Numerical simulations of X-rays Free Electron Lasers (XFEL)
Paolo Antonelli; Agissilaos Athanassoulis; Zhongyi Huang; Peter A. Markowich
2014-06-17T23:59:59.000Z
We study a nonlinear Schr\\"odinger equation which arises as an effective single particle model in X-ray Free Electron Lasers (XFEL). This equation appears as a first-principles model for the beam-matter interactions that would take place in an XFEL molecular imaging experiment in \\cite{frat1}. Since XFEL is more powerful by several orders of magnitude than more conventional lasers, the systematic investigation of many of the standard assumptions and approximations has attracted increased attention. In this model the electrons move under a rapidly oscillating electromagnetic field, and the convergence of the problem to an effective time-averaged one is examined. We use an operator splitting pseudo-spectral method to investigate numerically the behaviour of the model versus its time-averaged version in complex situations, namely the energy subcritical/mass supercritical case, and in the presence of a periodic lattice. We find the time averaged model to be an effective approximation, even close to blowup, for fast enough oscillations of the external field. This work extends previous analytical results for simpler cases \\cite{xfel1}.
Numerical evaluation of multi-gluon amplitudes for High Energy Factorization
Bury, M
2015-01-01T23:59:59.000Z
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.
Probing Modified Gravity with Atom-Interferometry: a Numerical Approach
Sandrine Schlogel; Sebastien Clesse; Andre Fuzfa
2015-07-11T23:59:59.000Z
Refined constraints on chameleon theories are calculated for atom-interferometry experiments, using a numerical approach consisting in solving for a four-region model the static and spherically symmetric Klein-Gordon equation for the chameleon field. By modeling not only the test mass and the vacuum chamber but also its walls and the exterior environment, the method allows to probe new effects on the scalar field profile and the induced acceleration of atoms. In the case of a weakly perturbing test mass, the effect of the wall is to enhance the field profile and to lower the acceleration inside the chamber by up to one order of magnitude. In the thin-shell regime, significant deviations from the analytical estimations are found, even when measurements are realized in the immediate vicinity of the test mass. Close to the vacuum chamber wall, the acceleration becomes negative and potentially measurable. This prediction could be used to discriminate between fifth-force effects and systematic experimental uncertainties, by doing the experiment at several key positions inside the vacuum chamber. The influence of the wall thickness and density is also studied. For the chameleon potential $V(\\phi) = \\Lambda^{4+\\alpha} / \\phi^\\alpha$ and a coupling function $A(\\phi) = \\exp(\\phi /M)$, one finds $M \\gtrsim 7 \\times 10^{16}$ GeV, independently of the power-law index. For $V(\\phi) = \\Lambda^4 (1+ \\Lambda/ \\phi)$ one finds $M \\gtrsim 4 \\times 10^{16}$ GeV. Future experiments able to measure an acceleration $a \\sim 10^{-11} \\mathrm{m/s^2}$ would probe the chameleon parameter space up to the Planck scale. Our method can easily be extended to constrain other models with a screening mechanism, such as symmetron, dilaton and f(R) theories.
Mathiesen, Patrick James
2013-01-01T23:59:59.000Z
Evaluation of numerical weather prediction solar irradiancecycle: The RUC. Monthly Weather Review, 132 (2), 495-518.representations. Monthly Weather Review, 139 (6), 1972-1995.
Ismagilov, Rustem F.
Effects of Shear Rate on Propagation of Blood Clotting Determined Using Microfluidics and Numerical-ismagilov@uchicago.edu 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
Energy and enstrophy transfer in numerical simulations of two-dimensional' turbulence
Vallis, Geoff
Energy and enstrophy transfer in numerical simulations of two-dimensional' turbulence Mathew E a significant fraction of the flow field,w and energy spectra from these simulations have slopes significantly October 1992; accepted 25 March 1993) Numerical simulations of statistically steady two-dimensional (2-D
NUMERICAL SIMULATION OF POOL BOILING FOR STEADY STATE AND TRANSIENT HEATING
Maruyama, Shigeo
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
Higham, Nicholas J.
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
Numerical analysis of a one-dimensional elastodynamic model of dry friction and unilateral contact
Renard, Yves - Pôle de Mathématiques, Institut National des Sciences Appliquées de Lyon
.e. hyperbolic) model with dry friction. Since we consider a Coulomb friction law with a slip velocity dependentNumerical analysis of a one-dimensional elastodynamic model of dry friction and unilateral contact in the numerical analysis of more elaborated dynamic purely elastic problems with dry friction. Ó 2001 Elsevier
Penetration of solar radiation in the upper ocean: A numerical model for oceanic and coastal waters
Lee, Zhongping
Penetration of solar radiation in the upper ocean: A numerical model for oceanic and coastal waters (2005), Penetration of solar radiation in the upper ocean: A numerical model for oceanic and coastal; Siegel et al., 1995] have demonstrated that the penetration of EVIS in the upper layer of the ocean plays
Kirby, James T.
Numerical simulation of the 2011 Tohoku tsunami: Comparison with field observations and sensitivity history, created a major tsunami that caused numerous deaths and enormous destruction on the nearby Hon- shu coast. Various tsunami sources were developed for this event, based on inverting seismic or GPS
NUMERICAL MODELING OF COASTAL TSUNAMI DISSIPATION AND Stphan T. Grilli1
Grilli, Stéphan T.
1 NUMERICAL MODELING OF COASTAL TSUNAMI DISSIPATION AND IMPACT Stéphan T. Grilli1 , Jeffrey C Tehranirad2 Recent observations of the coastal impact of large tsunamis (e.g., Indian Ocean 2004; Tohoku 2011) and related numerical and theoretical works have made it increasingly clear that tsunami waves arrive
NUMERICAL MODELING OF COASTAL TSUNAMI IMPACT DISSIPATION Stphan T. Grilli1
Kirby, James T.
1 NUMERICAL MODELING OF COASTAL TSUNAMI IMPACT DISSIPATION AND IMPACT Stéphan T. Grilli1 , Jeffrey Tehranirad2 Recent observations of the coastal impact of large tsunamis (e.g., Indian Ocean 2004; Tohoku 2011) and related numerical and theoretical works have made it increasingly clear that tsunami waves arrive
Numerical Investigation of a Transverse Jet in a Supersonic Crossflow using Large Eddy Simulation
Paris-Sud XI, Université de
Numerical Investigation of a Transverse Jet in a Supersonic Crossflow using Large Eddy Simulation injected fuel jets in hot supersonic crossflows. This paper describes the numerical algorithm being used into a supersonic crossflow computed on a coarse mesh. These results are discussed and similarity