The Influence of Basalt Layers on Seismic Wave Propagation
Hanssen, Peter
Large areas of the Earth are covered by intrusive and extrusive lavas, forming high velocity basalts, which are often opaque to cinventional seismicreflection surveys. On the North East Atlantic margin these basalts cover ...
Basaltic island sand provenance
Marsaglia, K.M. . Dept. of Geological Sciences)
1992-01-01
The Hawaiian Islands are an ideal location to study basaltic sand provenance in that they are a series of progressively older basaltic shield volcanoes with arid to humid microclimates. Sixty-two sand samples were collected from beaches on the islands of Hawaii, Maui, Oahu and Kauai and petrographically analyzed. The major sand components are calcareous bioclasts, volcanic lithic fragments, and monomineralic grains of dense minerals and plagioclase. Proportions of these components vary from island to island, with bioclastic end members being more prevalent on older islands exhibiting well-developed fringing reef systems and volcanic end members more prevalent on younger, volcanically active islands. Climatic variations across the island of Hawaii are reflected in the percentage of weathered detritus, which is greater on the wetter, northern side of the island. The groundmass of glassy, basaltic lithics is predominantly black tachylite, with lesser brown sideromelane; microlitic and lathwork textures are more common than holohyaline vitric textures. Other common basaltic volcanic lithic fragments are holocrystalline aggregates of silt-sized pyroxene or olivine, opaque minerals and plagioclase. Sands derived from alkalic lavas are texturally and compositionally indistinguishable from sands derived from tholeiitic lavas. Although Hawaiian basaltic sands overlap in composition with magmatic arc-derived sands in terms of their relative QFL, QmPK and LmLvLs percentages, they are dissimilar in that they lack felsic components and are more enriched in lathwork volcanic lithic fragments, holocrystalline volcanic lithic fragments, and dense minerals.
BasaltBasalt Rock units Cased interval Well completion Slotted casing Submersible pump Water level bls. Submersible pump intake near 262 ft bls; depth to water is 245.24 taken October 2, 2014. #12;GAM. Submersible pump intake near 604 ft bls. Depth to water is 600.32 taken October 9, 2014. #12;GAM(NAT) 0 75API
Bennecke, W.M.
1996-10-01
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.
Laske, Gabi
fracture in a rock Fissure eruption extrusion of lava along a fissure Flood basalt an extensive flow and steam from which minerals precipitate onto surrounding surfaces; the temperature in a fumarole typically
Lithium isotopes in global mid-ocean ridge basalts Paul B. Tomascak a,*, Charles H. Langmuir b
Langmuir, Charles H.
Lithium isotopes in global mid-ocean ridge basalts Paul B. Tomascak a,*, Charles H. Langmuir b January 2008 Abstract The lithium isotope compositions of 30 well-characterized samples of glassy lavas
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
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 QA:QAsource History View NewTexas:Montezuma, Arizona: Energy ResourcesProject |433014Â°, -82.9931607Â° ShowLava Dome
Injection and Monitoring at the Wallula Basalt Pilot Project
McGrail, B. Peter; Spane, Frank A.; Amonette, James E.; Thompson, Christopher J.; Brown, Christopher F.
2014-01-01
Continental flood basalts represent one of the largest geologic structures on earth but have received comparatively little attention for geologic storage of CO2. Flood basalt lava flows have flow tops that are porous, permeable, and have large potential capacity for storage of CO2. In appropriate geologic settings, interbedded sediment layers and dense low-permeability basalt rock flow interior sections may act as effective seals allowing time for mineralization reactions to occur. Previous laboratory experiments showed the relatively rapid chemical reaction of CO2-saturated pore water with basalts to form stable carbonate minerals. However, recent laboratory tests with water-saturated supercritical CO2 show that mineralization reactions occur in this phase as well, providing a second and potentially more important mineralization pathway than was previously understood. Field testing of these concepts is proceeding with drilling of the world’s first supercritical CO2 injection well in flood basalt being completed in May 2009 near the township of Wallula in Washington State and corresponding CO2 injection permit granted by the State of Washington in March 2011. Injection of a nominal 1000 MT of CO2 was completed in August 2013 and site monitoring is in progress. Well logging conducted immediately after injection termination confirmed the presence of CO2 predominantly within the upper flow top region, and showed no evidence of vertical CO2 migration outside the well casing. Shallow soil gas samples collected around the injection well show no evidence of leakage and fluid and gas samples collected from the injection zone show strongly elevated concentrations of Ca, Mg, Mn, and Fe and 13C/18O isotopic shifts that are consistent with basalt-water chemical reactions. If proven viable by this field test and others that are in progress or being planned, major flood basalts in the U.S., India, and perhaps Australia would provide significant additional CO2 storage capacity and additional geologic sequestration options in regions of these countries where conventional storage options are limited.
Injection and Monitoring at the Wallula Basalt Pilot Project
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
McGrail, B. Peter; Spane, Frank A.; Amonette, James E.; Thompson, Christopher J.; Brown, Christopher F.
2014-01-01
Continental flood basalts represent one of the largest geologic structures on earth but have received comparatively little attention for geologic storage of CO2. Flood basalt lava flows have flow tops that are porous, permeable, and have large potential capacity for storage of CO2. In appropriate geologic settings, interbedded sediment layers and dense low-permeability basalt rock flow interior sections may act as effective seals allowing time for mineralization reactions to occur. Previous laboratory experiments showed the relatively rapid chemical reaction of CO2-saturated pore water with basalts to form stable carbonate minerals. However, recent laboratory tests with water-saturated supercritical CO2 show thatmore »mineralization reactions occur in this phase as well, providing a second and potentially more important mineralization pathway than was previously understood. Field testing of these concepts is proceeding with drilling of the world’s first supercritical CO2 injection well in flood basalt being completed in May 2009 near the township of Wallula in Washington State and corresponding CO2 injection permit granted by the State of Washington in March 2011. Injection of a nominal 1000 MT of CO2 was completed in August 2013 and site monitoring is in progress. Well logging conducted immediately after injection termination confirmed the presence of CO2 predominantly within the upper flow top region, and showed no evidence of vertical CO2 migration outside the well casing. Shallow soil gas samples collected around the injection well show no evidence of leakage and fluid and gas samples collected from the injection zone show strongly elevated concentrations of Ca, Mg, Mn, and Fe and 13C/18O isotopic shifts that are consistent with basalt-water chemical reactions. If proven viable by this field test and others that are in progress or being planned, major flood basalts in the U.S., India, and perhaps Australia would provide significant additional CO2 storage capacity and additional geologic sequestration options in regions of these countries where conventional storage options are limited.« less
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: Alternative Fuels Data Center Homesum_a_epg0_fpd_mmcf_m.xls" ,"Available from WebQuantityBonneville Power Administration wouldMass map shinesSolarNewsusceptometer under pressureNavyNumerical simulations of current
Analytical-Numerical Modeling Of Komatiite Lava Emplacement And...
erosion becomes more effective for substrates that are increasingly unconsolidated and water rich. For thermo-mechanical erosion to be responsible for the formation of the...
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 QA:QA J-E-1 SECTION JEnvironmental Jump to:EAandAmminex A S Jump to: navigation,Inof Ground Source Heat| Open1986)
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 QA:QA J-E-1 SECTION J APPENDIXsource History ViewInformationWinds Jump to:Laredo Ridge WindHill JumpLava Creek
C. Doc Richardson; Nancy W. Hinman; Lindsay J. McHenry; J. Michelle Kotler; Jill R. Scott
2012-05-01
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.
Seismic Structure And Seismicity Of The Cooling Lava Lake Of...
Seismic Structure And Seismicity Of The Cooling Lava Lake Of Kilauea Iki, Hawaii Jump to: navigation, search OpenEI Reference LibraryAdd to library Journal Article: Seismic...
Aciego, S.M.; Jourdan, F.; DePaolo, D.J.; Kennedy, B.M.; Renne, P.R.; Sims, K.W.W.
2009-10-01
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.
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
Ecology of microbe/basaltic glass interactions : mechanisms and diversity
Sudek, Lisa A.
2011-01-01
nutrient and energy sources to support microbial growth.primary energy source for some basalt-associated microbialmicrobial communities inhabiting basalt surfaces to a greater extent than nutrients and energy
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-01
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.
Repository site definition in basalt: Pasco Basin, Washington
Guzowski, R.V.; Nimick, F.B.; Muller, A.B.
1982-03-01
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.
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 QA:QAsource History View NewTexas:Montezuma, Arizona: Energy ResourcesProject |433014Â°, -82.9931607Â° ShowLavaEnergy
A Primitive Alkali Basaltic Stratovolcano And Associated Eruptive...
Significance Jump to: navigation, search OpenEI Reference LibraryAdd to library Journal Article: A Primitive Alkali Basaltic Stratovolcano And Associated Eruptive Centres,...
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
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
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
THERMAL PROPERTIES OF GABLE MOUNTAIN BASALT CORES HANFORD NUCLEAR RESERVATION
Martinez-Baez, L.F.
2011-01-01
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
Insulation from basaltic stamp sand. Final technical report
Williams, F. D.
1981-04-01
A Midwest Appropriate Technology Grant was awarded to determine the technical and economic feasibility of producing mineral-fiber insulation directly from extensive deposits of basaltic sand produced during former mining and milling operations in the Keweenaw Peninsula region of Michigan's Upper Peninsula. The amounts of local basaltic sands available and representative chemical compositions were determined. The variation of viscosity with temperature and chemical composition was estimated. Samples were melted and either pulled or blown into fiber. In all cases fiber could be made with a reasonable tensile strength to ensure usefulness. It was concluded that it was technically feasible to produce fibers from basaltic stamp sands of the Upper Peninsula of Michigan. A technical feasibility study using published data, a cost and design analysis of a basalt fiber production plant, a market survey of fiber needs, and an economic analysis for investing in a basalt fiber venture was undertaken. These studies concluded that the local production of basaltic insulation was both feasible and economically reasonable. It was suggested that the plant be located in a region of greater population density with lower utility costs. A representative one-third of these studies is included as appendices A, B, C, and D.
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-01
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-01
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.
LAVA/CIS Version 2. 0: A software system for vulnerability and risk assessment
Smith, S.T.
1990-01-01
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.
Pedogenic Thresholds and Soil Process Domains in Basalt-Derived Soils
Vitousek, PM; Chadwick, OA
2013-01-01
rejuvenation of weathering-derived nutri- ent supply in anProcess Domains in Basalt-Derived Soils Peter M. Vitousekand domains in basalt-derived soils on two rainfall
Lava channel formation during the 2001 eruption on Mount Etna: evidence for mechanical erosion
Carmelo Ferlito; Jens Siewert
2005-12-14
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.
-destructive method using quantitative petrography of lunar basaltic samples Clive R. Neal a,b,c, , Patrick Donohue a
The Hadley-Apennine KREEP basalt igneous province G. Jeffrey TAYLOR1*
Spudis, Paul D.
-examine the geological setting of KREEP basalts at the Apollo 15 site, describe their petrography and mineral chemistry
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
O'Hara, Mhairi
2012-08-08
are in extremely good agreement with the real flow fields, with an average accuracy of 94.42%. The Lava Flow model is also used to produce hazard maps for both Hekla and Eyjafjallajökull, by identifying zones on the volcanoes flanks that will confine the lava...
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-14
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.
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
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
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
Suckale, Jenny
2011-01-01
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 ...
Spane, F.A. Jr.; Raymond, R.G.
1993-09-01
This report presents the first comprehensive Hanford Site-wide potentiometric map for the upper-basalt confined aquifer system (i.e., the upper Saddle Mountains Basalt). In constructing the potentiometric map, over forty on-site and off-site monitoring wells and boreholes were used. The potentiometric map developed for the upper-basalt confined aquifer is consistent with the areal head pattern indicated for the Mabton interbed, which is a deeper and more areally extensive confined aquifer underlying the Hanford Site. Salient features for the upper-basalt confined aquifer system potentiometric map are described.
Final Reclamation Report: Basalt Waste Isolation Project exploratory shaft site
Brandt, C.A.; Rickard, W.H. Jr.
1990-06-01
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.
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
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
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
Correlation of Cycles in Lava Lake Motion and Degassing at Erebus Volcano, Antarctica
Peters, Nial; Oppenheimer, Clive; Killingsworth, Drea Rae; Frechette, Jed; Kyle, Philip
2014-08-19
pulsatory behaviour either through experiments or computer simulations.111 D R A F T July 23, 2014, 11:09am D R A F T PETERS ET AL.: CYCLIC BEHAVIOUR OF THE EREBUS LAVA LAKE X - 7 2. Summary of Activity Although in general terms the inter-annual variability... to minimise the distance above the lake that211 we are measuring.212 The weather conditions at the summit of Erebus on 14 December 2013 were clear, with213 no clouds within the fields of view of the spectrometers and light winds of between 1–214 3 m s?1...
Lava Hot Springs Pool & Spa 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 QA:QAsource History View NewTexas:Montezuma, Arizona: Energy ResourcesProject |433014Â°, -82.9931607Â° ShowLava
Nakata, Dorene Samantha
2010-01-01
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 ...
Jackson, Matthew G. (Matthew Gerald)
2008-01-01
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 ...
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
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
Schuler, Juerg; Christie, Philip A. F.; White, Robert S.
2015-04-20
large uncertainties were associated with estimating constant phase shifts of the seismic arrivals within the topmost 250 m of the basalt sequences, where complex scattering occurred. Within this interval of the Brugdan I well, a phase-only compensation...
Basalt petrogenesis beneath slow- and ultraslow-spreading Arctic mid-ocean ridges
Elkins, Lynne J
2009-01-01
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 ...
Hammer, Julia Eve
, and the kinetics of silicate mineral nucleation and growth. Oxide abundance and magnetic intensity are most typical terrestrial basalt [Connerney et al., 1999]. [3] Most rock magnetic studies have focused
Cervantes, Pablo
1999-01-01
and their H?O and CO? concentrations. The compositions of the inclusions range from basalt to andesite and include both alkaline and calc-alkaline varieties. Most inclusions have low H?O concentrations (average values 0.4[]0.2 wt.% H?O for Xitle, 0.6[]0.5 wt...
Reidel, Steve P.; Spane, Frank A.; Johnson, Vernon G.
2002-08-08
This report provides the technical background and a guide to characterizing a site for storing natural gas in the Columbia River Basalt
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
Aciego, S.M.
2010-01-01
E. (1977) Subcommission on Geochronology - Convention on Useof Decay Constants in Geochronology and Cosmochronology.U-Th/He and 40 Ar/ 39 Ar Geochronology of Post-shield Lavas
Spane, F.A. Jr.; Webber, W.D.
1995-09-01
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-01
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
stratigraphy of plagioclase Hejiu Hui , Jocelyn G. Oshrin, Clive R. Neal Department of Civil Engineering-Al basaltic melts was studied using crystal stratigraphy, which involves textural (crystal size distributions the importance of crystal stratigraphy in revealing the intricacies of lunar basalt petrogenesis. Ó 2011 Elsevier
A lithospheric instability origin for Columbia River flood basalts and Wallowa Mountains uplift in
Humphreys, Eugene
A lithospheric instability origin for Columbia River flood basalts and Wallowa Mountains uplift-eruptive sub- sidence took place in the Wallowa Mountains, followed by syn- eruptive uplift of several hundred metres and a long-term uplift of about 2 km. The mapped surface uplift mimics regional topo- graphy
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-01
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...
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
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
Brandt, C.A.; Rickard, W.H. Jr.; Biehert, R.W.; Newell, R.L.; Page, T.L.
1989-01-01
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-02
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-01
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.
Johnson, Phillip Ray, II
2007-04-25
This thesis presents a material-centered characterization of 120 geologic samples from four fine-grained basalt quarries on the Samoan Island of Tutuila. Previous unsuccessful attempts at definitive Tutuilan quarry ...
Intermediate-scale sodium-concrete reaction tests with basalt and limestone concrete
Hassberger, J.A.; Muhlestein, L.D.
1981-01-01
Ten tests were performed to investigate the chemical reactions and rate and extent of attack between sodium and basalt and limestone concretes. Test temperatures ranged from 510 to 870/sup 0/C (950 to 1600/sup 0/F) and test times from 2 to 24 hours. Sodium hydroxide was added to some of the tests to assess the impact of a sodium hydroxide-aided reaction on the overall penetration characteristics. Data suggest that the sodium penetration of concrete surfaces is limited. Penetration of basalt concrete in the presence of sodium hydroxide is shown to be less severe than attack by the metallic sodium alone. Presence of sodium hydroxide changes the characteristics of sodium penetration of limestone concrete, but no major differences in bulk penetration were observed as compared to penetration by metallic sodium.
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
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-01
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-01
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.
Reidel, S.P.; Ledgerwood, R.K.; Myers, C.W.; Jones, M.G.; Landon, R.D.
1980-03-01
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.
Spane, Frank A.
2013-04-29
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
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
Non-basaltic asteroidal magmatism during the earliest stages of solar system evolution: A view from Sciences, Univ. of New Mexico, Albuquerque, NM 87131, USA b Dept. of Civil Eng. and Geo. Sci., Univ. of Notre Dame, IN 46556, USA c Dept. of Earth and Planetary Sciences, Univ. of New Mexico, Albuquerque, NM
Kaminski, Edouard
Estimation of ash injection in the atmosphere by basaltic volcanic plumes: The case explosive eruptions, volcanic plumes inject ash into the atmosphere and may severely affect air traffic, as illustrated by the 2010 Eyjafjallajökull eruption. Quantitative estimates of ash injection can be deduced from
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
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-15
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.
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
Santelli, Cara M
2007-01-01
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, ...
Teague, L.S.
2010-01-01
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
Constable, Steve
-basalt hydrocarbon exploration with marine controlled source electromagnetic and magnetotelluric data G. Michael accepted January 2015 ABSTRACT The recent use of marine electromagnetic technology for exploration from 84 marine controlled source electromagnetic and magnetotelluric stations for imaging volcanic
M. N. Pace; R. C. Bartholomay (USGS); J. J. Rosentreter (ISU)
1999-07-01
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.
Site characterization report for the basalt waste isolation project. Volume II
1982-11-01
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.
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-01
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.
Souza, Paulo Victor Santos
2015-01-01
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.
Aldrich, Jeffrey Brooks
1983-01-01
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...
Kelley, Jill Lynn
2002-01-01
. These rocks, when combined with 0-7 Ma altered dredge basalts from the same region, provide an opportunity to study long-term alteration and to determine the most important intluencing factorts). Mgo content, loss on ignition, and amount of secondary minerals...
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 ...
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-24
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-24
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.
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
Busby, Cathy
spe438-20 page 1 1 Garrison, N.J., Busby, C.J., Gans, P.B., Putirka, K., and Wagner, D.L., 2008-Miocene Lovejoy flood basalt, northern California Noah J. Garrison Cathy J. Busby Phillip B. Gans Department basalt. #12;2 Garrison et al. spe438-20 page 2 INTRODUCTION Mid-Miocene volcanism in the northern Sierra
C. Doc Richardson; Nancy W. Hinman; Jill R. Scott
2013-10-01
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.
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 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
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
Disruptive Innovation in Numerical Hydrodynamics
Waltz, Jacob I.
2012-09-06
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 Simulations of Bouncing Jets
Bonito, Andrea; Lee, Sanghyun
2015-01-01
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.
Class Generation for Numerical Wind Atlases
Class Generation for Numerical Wind Atlases Risø National Laboratory Wind Energy Department Constructing a Numerical Wind Atlas 5 2.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . 13 2.4 The existing procedure at Risø . . . . . . . . . . . . . . . . . . . 13 3 Representing a Wind
An implicit numerical algorithm general relativistic hydrodynamics
A. Hujeirat
2008-01-09
An implicit numerical algorithm general relativistic hydrodynamics This article has been replaced by arXiv:0801.1017
Petroleum Engineering 301 Petroleum Engineering Numerical Methods
10 Petroleum Engineering 301 Petroleum Engineering Numerical Methods Credit 3: (2-3) Required for Juniors Catalog Description: Use of numerical methods in a variety of petroleum engineering problems methods. Prerequisites(s): MATH 308 Textbook Required: Numerical Methods for Engineers, 6th Edition
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-10
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.
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
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
Numerical Simulation of Free Standing Hybrid Risers
Hou, Tiancong
2014-08-13
and offloading (FPSO). The numerical scheme is mainly based on an existing in-house numerical code, known as COUPLE. In using COUPLE, URA is modeled as a rigid body. The URA connects a vertical steel riser and a flexible jumper, both of which are modeled by beam...
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
Boal, David
. (1992) for further reading. There are several approaches to the integration of differential equationsCHAPTER 1 - NUMERICAL INTEGRATION There are many techniques for performing numerical integrations or numerical calculation of derivatives. Another approach, which we show for completeness in Sec. 1.1, is based
Okawa, Hirotada
2013-01-01
Numerical relativity became a powerful tool to investigate the dynamics of binary problems with black holes or neutron stars as well as the very structure of General Relativity. Although public numerical relativity codes are available to evolve such systems, a proper understanding of the methods involved is quite important. Here we focus on the numerical solution of elliptic partial differential equations. Such equations arise when preparing initial data for numerical relativity, but also for monitoring the evolution of black holes. Because such elliptic equations play an important role in many branches of physics, we give an overview of the topic, and show how to numerically solve them with simple examples and sample codes written in C++ and Fortran90 for beginners in numerical relativity or other fields requiring numerical expertise.
Hirotada Okawa
2013-08-15
Numerical relativity became a powerful tool to investigate the dynamics of binary problems with black holes or neutron stars as well as the very structure of General Relativity. Although public numerical relativity codes are available to evolve such systems, a proper understanding of the methods involved is quite important. Here we focus on the numerical solution of elliptic partial differential equations. Such equations arise when preparing initial data for numerical relativity, but also for monitoring the evolution of black holes. Because such elliptic equations play an important role in many branches of physics, we give an overview of the topic, and show how to numerically solve them with simple examples and sample codes written in C++ and Fortran90 for beginners in numerical relativity or other fields requiring numerical expertise.
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-06
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.
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids (2012)
Jameson, Antony
2012-01-01
dynamics, transitional or massively detached flows, turbulent mixing and combustion, or aerodynamic noise. CASTONGUAY AND A. JAMESON numerical methods available within the framework of industrial CFD applications
INDEX TO VOLUME 134 This index provides coverage for both the Initial Reports and Scientific Results
:388390 lava, A:153154 petrography, B:338, 342 petrology, B:353, 364367 andesite, basaltic, origin, A:16 andesite, porphyritic, blebs, A:338340 andesite, two-pyroxenes, petrography, A:338341 anions
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
Numerical simulations of strong incompressible magnetohydrodynamic turbulence
Mason, J.; Cattaneo, F.; Perez, J. C.; Boldyrev, S.
2012-05-15
Magnetised plasma turbulence pervades the universe and is likely to play an important role in a variety of astrophysical settings. Magnetohydrodynamics (MHD) provides the simplest theoretical framework in which phenomenological models for the turbulent dynamics can be built. Numerical simulations of MHD turbulence are widely used to guide and test the theoretical predictions; however, simulating MHD turbulence and accurately measuring its scaling properties is far from straightforward. Computational power limits the calculations to moderate Reynolds numbers and often simplifying assumptions are made in order that a wider range of scales can be accessed. After describing the theoretical predictions and the numerical approaches that are often employed in studying strong incompressible MHD turbulence, we present the findings of a series of high-resolution direct numerical simulations. We discuss the effects that insufficiencies in the computational approach can have on the solution and its physical interpretation.
Mechanical diode: Comparing numerical and experimental characterizations
Sagartz, M.J.; Segalman, D.; Simmermacher, T.
1998-02-01
In this introductory work, joint compliance is studied in both a numerical and experimental setting. A simple bolted interface is used as the test article and compliance is measured for the joint in both compression and in tension. This simple interface is shown to exhibit a strong non-linearity near the transition from compression to tension (or vice-versa). Modeling issues pertaining to numerically solving for the compliance are addressed. It is shown that the model predictions, in spite of convergence being very sensitive to numerical artifacts of the interface model, are in good agreement with experimentally measured strains and joint compliances. The joint behavior is a mechanical analogy to a diode, i.e., in compression, the joint is very stiff, acting almost as a rigid link, while in tension the joint is relatively soft, acting as a spring.
High performance computing and numerical modelling
,
2014-01-01
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-02
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.
Numerical simulation of an electroweak oscillon
Graham, N. [Department of Physics, Middlebury College, Middlebury, Vermont 05753 (United States)
2007-10-15
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.
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng (2010)
Li, Shaofan
2010-01-01
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
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 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 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
Torres-Verdín, Carlos
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2006; 65: Electrostatics D. Pardo1,,, L. Demkowicz1, C. Torres-Verdín2 and L. Tabarovsky3 1Institute for Computational. The refinement strategy is an extension of a fully automatic, energy-norm based, hp-adaptive algorithm. We
Numerical Verification of Bounce Harmonic Resonances in Neoclassical...
Office of Scientific and Technical Information (OSTI)
Numerical Verification of Bounce Harmonic Resonances in Neoclassical Toroidal Viscosity for Tokamaks Citation Details In-Document Search Title: Numerical Verification of Bounce...
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
13.024 Numerical Marine Hydrodynamics, Spring 2003
Milgram, Jerome H.
Introduction to numerical methods: interpolation, differentiation, integration, systems of linear equations. Solution of differential equations by numerical integration, partial differential equations of inviscid hydrodynamics: ...
Development of Numerical Simulation Capabilities for In Situ...
Office of Scientific and Technical Information (OSTI)
Development of Numerical Simulation Capabilities for In Situ Heating of Oil Shale Citation Details In-Document Search Title: Development of Numerical Simulation Capabilities for In...
Numerical simulations for low energy nuclear reactions including...
Office of Scientific and Technical Information (OSTI)
Numerical simulations for low energy nuclear reactions including direct channels to validate statistical models Citation Details In-Document Search Title: Numerical simulations for...
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
Numerical Verification of Bounce Harmonic Resonances in Neoclassical
Office of Scientific and Technical Information (OSTI)
for Tokamaks Kimin Kim, Jong-Kyu Park and Allen H. Boozer 70 PLASMA PHYSICS AND FUSION TECHNOLOGY Tokamaks, Numerical Verification Tokamaks, Numerical Verification This...
Stottlemyre, J.A.; Petrie, G.M.; Benson, G.L.; Zellmer, J.T.
1981-01-01
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.
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.
Two physical characteristics of numerical apparent horizons
Ivan Booth
2007-12-20
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.
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
Lateral boundary errors in regional numerical weather
?umer, Slobodan
Lateral boundary errors in regional numerical weather prediction models Author: Ana Car Advisor, they describe evolution of atmospher - weather forecast. Every NWP model solves the same system of equations (1: assoc. prof. dr. Nedjeljka Zagar January 5, 2015 Abstract Regional models are used in many national
Numerical likelihood analysis of cosmic ray anisotropies
Carlos Hojvat et al.
2003-07-02
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
Computer Algebra Applications for Numerical Relativity
Sascha Husa; Christiane Lechner
2003-01-21
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.
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 Model Construction with Closed Observables
Felix Dietrich; Gerta Köster; Hans-Joachim Bungartz
2015-10-18
Performing analysis, optimization and control using simulations of many-particle systems is computationally demanding when no macroscopic model for the dynamics of the variables of interest is available. In case observations on the macroscopic scale can only be produced via legacy simulator code or live experiments, finding a model for these macroscopic variables is challenging. In this paper, we employ time-lagged embedding theory to construct macroscopic numerical models from output data of a black box, such as a simulator or live experiments. Since the state space variables of the constructed, coarse model are dynamically closed and observable by an observation function, we call these variables closed observables. The approach is an online-offline procedure, as model construction from observation data is performed offline and the new model can then be used in an online phase, independent of the original. We illustrate the theoretical findings with numerical models constructed from time series of a two-dimensional ordinary differential equation system, and from the density evolution of a transport-diffusion system. Applicability is demonstrated in a real-world example, where passengers leave a train and the macroscopic model for the density flow onto the platform is constructed with our approach. If only the macroscopic variables are of interest, simulation runtimes with the numerical model are three orders of magnitude lower compared to simulations with the original fine scale model. We conclude with a brief discussion of possibilities of numerical model construction in systematic upscaling, network optimization and uncertainty quantification.
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
Observational and Numerical of the Adriatic Bora
Gohm, Alexander
and radiosoundings) 2. Mesoscale numerical model RAMS: · 5 or 6 nested grids · Horizontal mesh size as low as 800 m RAMS domain 5 #12;The event 28 March 2002: · Deep bora: deep NE flow, no critical level field: · The bora flow weakens and separates during daytime. RAMS vertical transects across Vratnik Pass
Volcanic Processes, the Resulting Landforms, and What These Mean To the Residents Nearby
Polly, David
amounts of silica. Silica molecules tend to form long polymer chains that aren't easily broken. High of silica the less polymer chains that form on a molecular level and the less dangerous the lava will be. Low viscosity magmas/lavas would form flood basalts and shield volcanos. Intermediate magmas
Fast numerical test of hyperbolic chaos
Pavel V. Kuptsov
2011-11-21
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-01
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.
A numerical study of pseudoscalar inflation
Cheng, Shu-Lin; Ng, Kin-Wang
2015-01-01
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.
A numerical study of pseudoscalar inflation
Shu-Lin Cheng; Wolung Lee; Kin-Wang Ng
2015-08-02
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.
Lattice Boltzmann Model for Numerical Relativity
Ilseven, E
2015-01-01
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.
Numerical Simulations Unravel the Cosmic Web
C. -A. Faucher-Giguere; A. Lidz; L. Hernquist
2008-03-03
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.
Beyond Numerics: The Existence of Pure Filtrations
Eisenbud, David; Schreyer, Frank-Olaf
2010-01-01
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.
Category:Numerical Modeling | 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 QA:QA J-E-1 SECTION JEnvironmentalBowerbank,CammackFLIR Jump to: navigation,GroundNear Infrared SurveysNumerical
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
Numerical Results for the Blue Phases
G. P. Alexander; J. M. Yeomans
2009-06-04
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.
RELAP-7 Numerical Stabilization: Entropy Viscosity Method
R. A. Berry; M. O. Delchini; J. Ragusa
2014-06-01
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.
GRChombo : Numerical Relativity with Adaptive Mesh Refinement
Clough, Katy; Finkel, Hal; Kunesch, Markus; Lim, Eugene A; Tunyasuvunakool, Saran
2015-01-01
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 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
Numerical Modeling of Eastern Connecticut's Visual Resources1
Numerical Modeling of Eastern Connecticut's Visual Resources1 Daniel L. Civco 2/ l Submitted Conservation, The University of Connecticut, Storrs, Connecticut. / Abstract: A numerical model capable of accurately predicting the preference for landscape photographs of selected points in eastern Connecticut
Testing and improving the numerical accuracy of the NLO predictions
R. Pittau
2010-09-24
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 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
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
NUMERICAL APPROXIMATION OF A MACROSCOPIC MODEL OF PEDESTRIAN FLOWS
Chalons, Christophe
NUMERICAL APPROXIMATION OF A MACROSCOPIC MODEL OF PEDESTRIAN FLOWS CHRISTOPHE CHALONS Abstract for the description of the flow of pedestrians. Solutions of the associated Riemann problem are known to be possibly in obtaining an efficient numerical scheme. Numerical evidences are proposed. Key words. macroscopic pedestrian
Numerical computation of elastoacoustic vibrations with interface damping
Rodríguez, Rodolfo
is shown to be well posed and optimal error estimates are proved for the numerical method. Numerical Babuska[1] for indefinite variational problems. Finally numerical results are given for some real test problem Let us consider an elastic (two-dimensional) vessel enclosing an inviscid com- pressible
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
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl. (in press)
Pasciak, Joseph
and is in the public domain in the U.S.A. Contract/grant sponsor: U.S. Department of Energy by the University. The latter assumption is important since it allows for easy construction of nested multilevel spaces on regular auxiliary meshes. Numerical experiments in both two and three space dimensions illustrate
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-01
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-01
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.
Numerical Analysis of Black Hole Evaporation
Tsvi Piran; Andrew Strominger
1993-04-28
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.
NUMERICAL MODELING OF CATHODE CONTACT MATERIAL DENSIFICATION
Koeppel, Brian J.; Liu, Wenning N.; Stephens, Elizabeth V.; Khaleel, Mohammad A.
2011-11-01
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.
Direct numerical simulation of turbulent reacting flows
Chen, J.H. [Sandia National Laboratories, Livermore, CA (United States)
1993-12-01
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.
Smith, D.K.
1990-01-01
Mineralogical, textural and compositional data accompanying greenschist facies metamorphism (to 300{degrees}C) of basalts of the East Rift Zone (ERZ), Kilauea, Hawaii may be evaluated relative to published and experimental results for the surface corrosion of borosilicate glass. The ERZ alteration sequence is dominated by intermittent palagonite, interlayered smectite-chlorite, chlorite, and actinolite-epidote-anhydrite. Alteration is best developed in fractures and vesicles where surface reaction layers root on the glass matrix forming rinds in excess of 100 microns thick. Fractures control fluid circulation and the alteration sequence. Proximal to the glass surface, palagonite, Fe-Ti oxides and clays replace fresh glass as the surface reaction layer migrates inwards; away from the surface, amphibole, anhydrite, quartz and calcite crystallize from hydrothermal fluids in contact with the glass. The texture and composition of basaltic glass surfaces are similar to those of a SRL-165 glass leached statically for sixty days at 150 {degrees}C. While the ERZ reservoir is a complex open system, conservative comparisons between the alteration of ERZ and synthetic borosilicate glass are warranted. 31 refs., 2 figs.
Numerical studies of Phi^2-Oscillatons
Miguel Alcubierre; Ricardo Becerril; F. Siddhartha Guzman; Tonatiuh Matos; Dario Nunez; L. Arturo Urena-Lopez
2003-01-24
We present an exhaustive analysis of the numerical evolution of the Einstein-Klein-Gordon equations for the case of a real scalar field endowed with a quadratic self-interaction potential. The self-gravitating equilibrium configurations are called oscillatons and are close relatives of boson stars, their complex counterparts. Unlike boson stars, for which the oscillations of the two components of the complex scalar field are such that the spacetime geometry remains static, oscillatons give rise to a geometry that is time-dependent and oscillatory in nature. However, they can still be classified into stable (S-branch) and unstable (U-branch) cases. We have found that S-oscillatons are indeed stable configurations under small perturbations and typically migrate to other S-profiles when perturbed strongly. On the other hand, U-oscillatons are intrinsically unstable: they migrate to the S-branch if their mass is decreased and collapse to black holes if their mass is increased even by a small amount. The S-oscillatons can also be made to collapse to black holes if enough mass is added to them, but such collapse can be efficiently prevented by the gravitational cooling mechanism in the case of diluted oscillatons.
Numerical treatment of disorder in PHC slabs
Minkov, Momchil
2015-01-01
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 ...
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-15
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.
Direct numerical simulations of aeolian sand ripples
Orencio Duran; Philippe Claudin; Bruno Andreotti
2014-11-07
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.
Self-similar radiation from numerical Rosenau-Hyman compactons
Rus, Francisco Villatoro, Francisco R.
2007-11-10
The numerical simulation of compactons, solitary waves with compact support, is characterized by the presence of spurious phenomena, as numerically induced radiation, which is illustrated here using four numerical methods applied to the Rosenau-Hyman K(p, p) equation. Both forward and backward radiations are emitted from the compacton presenting a self-similar shape which has been illustrated graphically by the proper scaling. A grid refinement study shows that the amplitude of the radiations decreases as the grid size does, confirming its numerical origin. The front velocity and the amplitude of both radiations have been studied as a function of both the compacton and the numerical parameters. The amplitude of the radiations decreases exponentially in time, being characterized by a nearly constant scaling exponent. An ansatz for both the backward and forward radiations corresponding to a self-similar function characterized by the scaling exponent is suggested by the present numerical results.
Numerical simulation of three-dimensional electrical flow through geomaterials
Akhtar, Anwar Saeed
1998-01-01
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-26
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.
A Numerical Comparison of Rule Ensemble Methods and Support Vector...
Office of Scientific and Technical Information (OSTI)
Title: A Numerical Comparison of Rule Ensemble Methods and Support Vector Machines Machine or statistical learning is a growing field that encompasses many scientific problems...
A FAST AND ACCURATE NUMERICAL SCHEME FOR THE ...
1999-10-21
[13] M. Jarraud and A. P. M. Baede, The use of spectral techniques in numerical weather prediction, in Large Scale Computations in Fluid Mechanics, B. Enquist,
Numerical simulations for width fluctuations in compound elastic...
Office of Scientific and Technical Information (OSTI)
simulations for width fluctuations in compound elastic and inelastic scattering at low energies Citation Details In-Document Search Title: Numerical simulations for width...
A Hydro-Thermo-Mechanical Numerical Model For Hdr Geothermal...
A Hydro-Thermo-Mechanical Numerical Model For Hdr Geothermal Reservoir Evaluation Jump to: navigation, search OpenEI Reference LibraryAdd to library Journal Article: A...
Numerical Modeling of the Nucleation Conditions of Petal-Centerline...
Numerical Modeling of the Nucleation Conditions of Petal-Centerline Fractures below a Borehole Floor, A Sensitivity Study and Application to the Coso Geothermal Field Jump to:...
Toward portable programming of numerical linear algebra on manycore...
Office of Scientific and Technical Information (OSTI)
Resource Relation: Conference: Toward petaflop numerical simulation on parallel hybrid architectures held June 6-10, 2011 in Sophia Antipolis, France.; Related Information:...
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
Numerical Modeling of PCCI Combustion | Department of Energy
Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site
PCCI Combustion Numerical Modeling of PCCI Combustion 2004 Diesel Engine Emissions Reduction (DEER) Conference Presentation: Lawrence Livermore National LaboratoryUniversity of...
Development of Numerical Simulation Capabilities for In Situ...
Office of Scientific and Technical Information (OSTI)
Numerical Simulation Capabilities for In Situ Heating of Oil Shale Hoda, Nazish ExxonMobil Upstream Research Company, Houston, TX, USA; Fang, Chen ExxonMobil Upstream Research...
An Iterated, Multipoint Differential Transform Method for Numerically...
Office of Scientific and Technical Information (OSTI)
Differential Transform Method for Numerically Evolving Partial Differential Equation Initial-Value Problems Citation Details In-Document Search Title: An Iterated,...
Numerical simulation of the environmental impact of hydraulic...
Office of Scientific and Technical Information (OSTI)
Numerical simulation of the environmental impact of hydraulic fracturing of tightshale gas reservoirs on near-surface groundwater: Background, base cases, shallow reservoirs,...
Optimization Online - On mutual impact of numerical linear algebra ...
Marco D'Apuzzo
2008-03-12
Mar 12, 2008 ... On mutual impact of numerical linear algebra and large-scale optimization with focus on interior point methods.
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 ...
Direct Numerical Simulations and Robust Predictions of Cloud...
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
cloud. Credit: Computational Science and Engineering Laboratory, ETH Zurich, Switzerland Direct Numerical Simulations and Robust Predictions of Cloud Cavitation Collapse PI Name:...
Numerical Data Advisory Board assembly of mathematial and physical sciences
Not Available
1980-07-30
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.
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
Numerical Study of Acoustic Modes in Ducted Shear Flow
Rienstra, Sjoerd W.
Numerical Study of Acoustic Modes in Ducted Shear Flow Gregory G. Vilenski & Sjoerd W. Rienstra mean flow inside a duct is studied numerically. For isentropic flow in a circular duct with zero swirl 26, 2007 #12;1 Introduction Normal mode analysis of small-amplitude disturbances in an annular duct
Validation of Eddy-renewal model by numerical simulation
Garbe, Christoph S.
Validation of Eddy-renewal model by numerical simulation Li-Ping Hung1,2 , Christoph S. Garbe1, E-mail: wttsai@ncu.edu.tw Abstract. The eddy-renewal model proposes that the dominant vortical flows as the diffusivity effect in the thin diffusive sublayer beneath surface. Key Words: eddy-renewal model, numerical
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
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
High-performance numerical optimization on multicore clusters
Dimakopoulos, Vassilios
High-performance numerical optimization on multicore clusters Panagiotis E. Hadjidoukas1 Department of Computer Science {phadjido,voglis,dimako,lagaris}@cs.uoi.gr 2 Department of Materials Science presents a software infrastructure for high per- formance numerical optimization on clusters of multicore
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
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
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
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
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
NUMERICAL MODELING OF LOW FREQUENCY HYDRO-ACOUSTIC WAVES
Kirby, James T.
NUMERICAL MODELING OF LOW FREQUENCY HYDRO-ACOUSTIC WAVES GENERATED BY SUBMARINE TSUNAMIGENIC#al to increase the reliability of the system Â· Can we use precursors of tsunami? Hydro numerical models applicable on an oceanic scale #12;Index Â· Introduc#on on hydro
Astrophysical jets: Observations, numerical simulations, and laboratory experiments
Bellan, Paul M.
Astrophysical jets: Observations, numerical simulations, and laboratory experiments P. M. Bellan,1; published online 22 April 2009 This paper provides summaries of ten talks on astrophysical jets given of observation, numerical modeling, and laboratory experiment. One essential feature of jets, namely
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
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 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
Mathiesen, Patrick James
2013-01-01
of numerical weather prediction solar irradiance forecastsof numerical weather prediction solar irradiance forecastsnumerical weather prediction model for solar irradiance
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
Numerical Simulations of MHD Turbulence in Accretion Disks
Steven A. Balbus; John F. Hawley
2002-03-20
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.}
Grid-independent Issue in Numerical Heat Transfer
Yao Wei; Wang Jian; Liao Guangxuan
2006-09-26
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.
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-09
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 simulation of micro-fluidic passive and active mixers
Kumar, Saurabh
2002-01-01
the flow rate. Hence, output channels exhibit predetermined concentration values, which allow concentration dependent chemistry experiments in each output channel. Convective diffusive transport in this micro mixer is studied numerically and theoretically...
Numerical models of phosphate esters in the Chattahoochee River
Haffey, Samuel Fraad, 1973-
2004-01-01
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 ...
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
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
Numerical Analysis of a Cold Air Distribution System
Zhu, L.; Li, R.; Yuan, D.
2006-01-01
non-uniform so that residents usually feel uncomfortable. The distribution of indoor airflow by cold air distribution is researched in this paper. We study indoor air distribution under different low temperature air supply conditions by numerical...
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
Assessment of the MUSTA approach for numerical relativistic dynamics
Blakely, P. M.; Nikiforakis, N.; Henshaw, W. D.
2015-03-04
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 Approximations of Stochastic Optimal Stopping and Control Problems
Siska, David
2007-01-01
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 and Experimental Investigation of Tidal Current Energy Extraction
Sun, Xiaojing
2008-01-01
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 Study of Flexible Flapping Wing Propulsion , Mingjun Wei
Wei, Mingjun
Numerical Study of Flexible Flapping Wing Propulsion Tao Yang , Mingjun Wei New Mexico State in 1912.4 The Knoller-Betz effect was verified experimentally in a wind tunnel by Katzmayr in 1922
Experimental and numerical study of primary consolidation of soft clay
Korchaiyapruk, Attasit, 1977-
2007-01-01
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. ...
An experimental and numerical study of wind turbine seismic behavior
Prowell, I.
2011-01-01
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
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 ...
Fully kinetic numerical modeling of a plasma thruster
Szabo, James Joseph, 1969-
2001-01-01
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 ...
Experimental and numerical analysis of a deepwater mini-TLP
Guichard, Aurelien
2001-01-01
As the quest for oil and gas resources drives the industry to ever deeper waters, model testing still represents an essential step after numerical modeling when designing offshore platforms in these hostile environments. ...
A numerical model simulation of longshore transport for Galveston Island
Gilbreath, Stephen Alexander
1995-01-01
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...
Numerical Chaos in a Fractional Order Logistic Map
Joakim Munkhammar
2010-11-10
In this paper we investigate a fractional order logistic map and its discrete time dynamics. We show some basic properties of the fractional logistic map and numerically study its period-doubling route to chaos.
Direct Numerical Simulation of the Flow in a Pebble Bed
Ward, Paul
2014-06-24
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 APPROXIMATIONS OF ALLEN-CAHN AND CAHN ...
2010-03-18
schemes that satisfy a corresponding discrete energy law, or in other words, energy ... Our first objective is to design stabilized semi-implicit schemes that satisfy an energy ...... Numerical analysis of a continuum model of phase transition.
Numerical Simulations of Gamma-Ray Burst Explosions
Lazzati, Davide; López-Cámara, Diego
2015-01-01
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...
Numerical simulation of flow separation control by oscillatory fluid injection
Resendiz Rosas, Celerino
2005-08-29
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...
Numerical Study of Abutment Scour in Cohesive Soils
Chen, Xingnian
2010-01-16
. Numerical simulations of overtopping flow in straight rectangular channel, straight compound channel and channel bend have been conducted. The bridge deck is found to be able to change the flow distribution and the bed shear stress will increase...
Numerical solutions of differential equations on FPGA-enhanced computers
He, Chuan
2009-05-15
, or even coding with assembly languages. With the emergence of Field Programmable Gate Array (FPGA) based Reconfigurable Computing (RC) technology, numerical scientists and engineers now have another option using FPGA devices as core components to address...
Numerical Weather Forecasting at the Savannah River Site
Buckley, R.L. [Westinghouse Savannah River Company, AIKEN, SC (United States)
1998-08-01
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.
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
Experimental and Numerical Study of Polymer Scratch Behavior
Jiang, Han
2010-10-12
As part of a larger effort to understand the fundamental knowledge of polymer scratch behavior, this dissertation is focused on both experimental study and numerical analysis of scratch deformation of a broad range of polymers, with an emphasis...
Numerical Simulation Study on Transpired Solar Air Collector
Wang, C.; Guan, Z.; Zhao, X.; Wang, D.
2006-01-01
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 analysis of vapor flow in a micro heat pipe
Liu, Xiaoqin
1996-01-01
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...
Notes 09. Numerical evaluation of natural modes and frequencies
San Andres, Luis
2008-01-01
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...
Nagurney, Anna
Contributions Motivation Literature Models Synergy Measure Numerical Examples Conclusions Supply Motivation Literature Models Synergy Measure Numerical Examples Conclusions 1 Contributions 2 Motivation 3 University of Massachusetts Amherst #12;Contributions Motivation Literature Models Synergy Measure Numerical
Kim, Yong Jung
Features and numerical Tools Laplace Biharmonic Helmholtz Maxwell Solving Fredholm second kind and numerical Tools Laplace Biharmonic Helmholtz Maxwell Outline Recursively Compressed Inverse Preconditioning. There are, of course, other methods. #12;Features and numerical Tools Laplace Biharmonic Helmholtz Maxwell
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
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
Tokamak elongation: how much is too much? II Numerical results
Lee, Jungpyo; Freidberg, Jeffrey P
2015-01-01
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...
Numerical Methods for a Nonlinear BVP Arising in Physical Oceanography
Riccardo Fazio; Alessandra Jannelli
2013-10-08
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-01
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...
Propagation of numerical noise in particle-in-cell tracking
Kesting, Frederik
2015-01-01
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 SIMULATION OF NATURAL GAS-SWIRL BURNER
Ala Qubbaj
2005-03-01
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.
Numerical simulation of water flow around a rigid fishing net
Roger Lewandowski; Géraldine Pichot
2006-12-20
This paper is devoted to the simulation of the flow around and inside a rigid axisymmetric net. We describe first how experimental data have been obtained. We show in detail the modelization. The model is based on a Reynolds Averaged Navier-Stokes turbulence model penalized by a term based on the Brinkman law. At the out-boundary of the computational box, we have used a "ghost" boundary condition. We show that the corresponding variational problem has a solution. Then the numerical scheme is given and the paper finishes with numerical simulations compared with the experimental data.
Numerical calculation of wave refraction by digital computer
Orr, Terry Edwin
1969-01-01
+~ 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 ...
Mathiesen, Patrick; Collier, Craig; Kleissl, Jan
2013-01-01
of numerical weather prediction solar irradiance forecasts numerical weather prediction model for solar irradiance weather prediction for intra?day solar forecasting in the
EXPERIMENTAL CHARACTERIZATION AND NUMERICAL MODELING OF THERMAL AND ELECTROCHEMISTRY EFFECTS IN
Kassegne, Samuel Kinde
EXPERIMENTAL CHARACTERIZATION AND NUMERICAL MODELING OF THERMAL AND ELECTROCHEMISTRY EFFECTS IN 3D and Numerical Modeling of Thermal and Electrochemistry Effects in 3D Bionanoelectronics Platform by Neha
Numerical Modeling At Neal Hot Springs Geothermal Area (U.S....
Area Exploration Technique Numerical Modeling Activity Date 2011 - 2011 Usefulness useful DOE-funding Unknown Exploration Basis A numerical reservoir model was created to...
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
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
Numerical Methods for the Inverse Nonlinear Fourier Transform
Civelli, Stella; Secondini, Marco
2015-01-01
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
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
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
ARTICULATORY SYNTHESIS: NUMERICAL SOLUTION OF A HYPERBOLIC DIFFERENTIAL EQUATION
into the context of finite-difference approximations to a differential equation describing acoustic wave-that of solving the differ- ential equation describing acoustic (small amplitude), one-dimensional propagationARTICULATORY SYNTHESIS: NUMERICAL SOLUTION OF A HYPERBOLIC DIFFERENTIAL EQUATION Richard S. Mc
NUMERICAL STUDY OF A NONLINEAR HEAT EQUATION FOR PLASMA PHYSICS
Filbet, Francis
, in particular the plasma-wall interactions as well as the occurring turbulence have an important impact of a tokamak. The nonlinearity implies some numerical difficulties, in particular long time behavior, when (SOL) of a tokamak, is nowadays one of the main problems for fusion generated energy production (ITER
NUMERICAL STUDY OF A NONLINEAR HEAT EQUATION FOR PLASMA PHYSICS
Negulescu, Claudia
of the tokamak, in particular the plasma-wall interactions as well as the occurring turbulence have an important of a tokamak. The nonlinearity implies some numerical difficulties, in particular for the long time behavior of the transport of charged particles confined in a tokamak is nowadays one of the main problems for fusion
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
Why we Couldn't Use Numerical Libraries for PETSc
Gropp, Bill
, Interoperability, Composition 1 INTRODUCTION This paper discusses some of the reasons for the crisis in numerical of Computational and Technology Research, U.S. Department of Energy, under Contract W31109Eng38. \\Lambda. A typical method may need to solve a linear system on the entire domain. Since the PDE discretization
Some Books Relevant to Computational Mathematics Numerical Differential Equations
Fasshauer, Greg
, 1997. 2. Briggs, W., A Multigrid Tutorial, SIAM, 1987. 3. Gear, C. W., Numerical Initial Value Problems and Surfaces for Computer Aided Geometric Design (3rd ed.), Academic Press, 1993. 12. Golub, G. H. and Ortega., Fundamentals of Computer Aided Geometric Design, A K Peters, 1993. 17. Isaacson, E. and Keller, H. B., Analysis
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
NUMERICAL SIMULATION OF AIR POLLUTION DYNAMICS DUE TO
Olszewski Jr., Edward A.
spots. Global weather patterns effected. CO2 and other pollutant emissions contributing to globalNUMERICAL SIMULATION OF AIR POLLUTION DYNAMICS DUE TO POINT SOURCE EMISSIONS FROM AN INDUSTRIAL, and ultimately effects the global climate balance. · About 60% of emissions from point sources · Major pollutants
Thermodynamics cycle analysis and numerical modeling of thermoelastic cooling systems
Rubloff, Gary W.
Thermodynamics cycle analysis and numerical modeling of thermoelastic cooling systems Suxin Qian level. However, a thermoelastic cooling system integrated with heat transfer fluid loops have not been;2012) (a.k.a. elastocaloric cooling). These solid-state cooling systems offer us alternatives to eliminate
NUMERICAL INVESTIGATION OF AN ANISOTHERMAL TURBULENT FLOW WITH EFFUSION
Nicoud, Franck
NUMERICAL INVESTIGATION OF AN ANISOTHERMAL TURBULENT FLOW WITH EFFUSION Simon Mendez CFD Team the plate. CONTEXT In gas turbines, the turbine blades and the liner of the combustion chamber are submitted. For economical reasons, computational fluid dynamics is now widely used as a design tool by combustion cham- ber
Flow Forcing Techniques for Numerical Simulation of Combustion Instabilities
Nicoud, Franck
Flow Forcing Techniques for Numerical Simulation of Combustion Instabilities A. KAUFMANN* and F of combustion instabilities in gas turbine combustors require the knowledge of flame transfer functions. Those flame) and for one case where a CFD code is necessary (a laminar Bunsen-type flame). © 2002
Numerical simulation of tsunami waves generated by deformable submarine landslides
Kirby, James T.
wave model Tsunami wave Numerical modeling a b s t r a c t This paper presents a new submarine energy is mostly concentrated on a narrow band of the dominant slide direction for the waves generated-up along the coast. For example, submarine mass failure is considered as one of the major sources
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
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, Inc., Cranston, RI, USA ABSTRACT We present the development and application of small buoy systems climates, of two new types of buoy systems equipped with an embedded Linear Electric Generator (LEG; made
On the numerical solution of some problems of environmental pollution
Ehrhardt, Matthias
On the numerical solution of some problems of environmental pollution Quang A Dang1 , Matthias world and their importance will even increase in the future. High pollution of air, water and soil may to the mathematical formula- tion of air pollution models. Let G be a cylindrical domain in the three
Simple intrinsic defects in InAs : numerical predictions.
Schultz, Peter Andrew
2013-03-01
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.
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
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
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
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
Black hole free energy during charged collapse: a numerical study
Hugues Beauchesne; Ariel Edery
2012-05-19
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.
INCORPORATION OF GROUNDWATER FLOW INTO NUMERICAL MODELS AND DESIGN MODELS
-coupled, ground-source heat pumps, groundwater, heat pump, heat exchanger, heat transfer, numerical models-loop ground-coupled heat exchangers. Green and Perry (1961) demonstrated that the value of effective thermal on the design and performance of vertical closed-loop ground heat exchangers. Based on the investigation results
Numerical Simulation in Applied Geophysics. From the Mesoscale to the
Santos, Juan
Seismic wave propagation is a common technique used in hydrocarbon exploration geophysics, mining's crust and induce attenuation, dispersion and anisotropy of the seismic waves observed at the macroscale process. Numerical Simulation in Applied Geophysics. From the Mesoscale to the Macroscale p. #12
Numerical Algorithms manuscript No. (will be inserted by the editor)
Gander, Martin J.
-time is performed through classical or optimized transmission conditions (TCs). We analyze in this pa- per with Robin TCs applied to the heat equation. After a general convergence analysis using energy estimates, we with a numerical experiment. Keywords Schwarz Waveform Relaxation Â· Optimized Robin Transmission Conditions Â· Short
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
Numerical Modeling of Periodic Composite Media for Electromagnetic Shielding Application
Koledintseva, Marina Y.
on a conventional mixing theory, have served as the fundamentals for these techniques. In these formulationsNumerical Modeling of Periodic Composite Media for Electromagnetic Shielding Application Dagang Wu-difference time-domain (FDTD) method. The results are compared with conventional mixing theories and 3D Fourier
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
Numerical Model of a Tensioner System and Flex Joint
Huang, Han
2013-07-27
D will be integrated into a numerical code, known as COUPLE, for the simulation of the dynamic interaction among the hull of a floating structure, such as SPAR or TLP, its mooring system and riser system under the impact of wind, current and waves...
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 well stimulation Fracturing Fluid Proppant #12;5 Quarries #12;6 Magma flow Tarkastad #12;7 Model EQ 1
Hydraulic Fractures: multiscale phenomena, asymptotic and numerical solutions
Peirce, Anthony
Hydraulic Fractures: multiscale phenomena, asymptotic and numerical solutions CSIRO CSS TCP Detournay (UMN) Eduard Siebrits (SLB) #12;2 Outline · Examples of hydraulic fractures · Governing equations well stimulation Fracturing Fluid Proppant #12;5 Quarries #12;6 Magma flow Tarkastad #12;7 Model EQ 1
Stress and diffusion induced interface motion: Modelling and numerical simulations
Styles, Vanessa
Stress and diffusion induced interface motion: Modelling and numerical simulations Harald Garcke of Mathematics, University of Sussex, Brighton, BN1 9QH, U.K. Abstract We propose a phase field model for stress stress effects. In this paper we will demonstrate that the model can also be used to describe other
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
Numerical Study of Flexible Flapping Wing Propulsion and Mingjun Wei
Wei, Mingjun
Numerical Study of Flexible Flapping Wing Propulsion Tao Yang and Mingjun Wei New Mexico State highly flexible flapping wings interacting with fluid flows. Here, the fluid motion, solid motion necessary flapping mechanism, control cells are implemented in solid area (i.e., the wing) as "skeleton
Direct Numerical Simulation of Solid Deformation During Dendritic Solidification
Beckermann, Christoph
solidification is a common phe- nomenon in metal casting and can lead to defects such as hot tears, macro is caused by external forces, for example through the rolls in continuous casting, mold wall movement of the mush can lead to numerous defects in a solidified casting, including hot tears, macrosegregation
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
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
Numerical Study of Noise Shielding by Airframe Structures Changzheng Huang*
Papamoschou, Dimitri
Numerical Study of Noise Shielding by Airframe Structures Changzheng Huang* and Dimitri Papamoschou-quiet advent aircraft that use jet noise shielding by the airframe. Current methods to predict shielding from predictive tools for jet noise shielding therefore requires a different approach. In this study we use
Landolt-Bornstein Numerical Data and Functional Relationships
Landolt-B¨ornstein Numerical Data and Functional Relationships NEW SERIES 4.1.1.4: Solar, the temperature structure of theoretical models is obtained by solving an energy equation. In the solar atmosphere Potsdam, Germany msteffen@aip.de October 29, 2009 #12;#12;Ref. p. 44] 4.1.1.4 Solar photosphere
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
NUMERICAL SOLUTION OF THE POLYMER SYSTEM BY FRONT TRACKING
, models twophase, threecomponent polymer flooding in an oil reservoir. Two examples are presented, where. Key words and phrases. Front tracking, nonstrictly hyperbolic systems, polymer flooding, dimensionalNUMERICAL SOLUTION OF THE POLYMER SYSTEM BY FRONT TRACKING V. HAUGSE, K. H. KARLSEN, K.--A. LIE
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 Simulations Concerning the Propagation of Protostellar Jets
Numerical Simulations Concerning the Propagation of Protostellar Jets A thesis submitted, Miruna, Natalie, Olena, Pat, Prakash, Ram, Rhona, Shane, Simon, Sri, Timur, Toby, Tolis, Tony, Tom;#12;Summary A protostellar jet is a highly supersonic stream of material which heralds the birth of a star
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
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
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 Code for LHCD Simulations with Self-consistent Treatment of Alpha Particles in Tokamak Geometry
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
Carmignani, B
2005-01-01
Numerical simulations of welds of thick steel pieces of interest for the thermonuclear fusion ITER machine
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
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-10
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 Regularization of Electromagnetic Quantum Fluctuations in Inhomogeneous Dielectric Media
Shin-itiro Goto; Alison C. Hale; Robin W. Tucker; Timothy J. Walton
2012-01-05
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 study of thermoacoustic convection in a cavity
Fusegi, Toru; Farouk, B.; Oran, E.S.
1995-12-31
Thermoacoustic convection in a two-dimensional cavity is numerically studied. Part of a compressible fluid (Helium) near the center line of the cavity is suddenly energized to generate pressure waves. Numerical solutions are secured by employing a highly accurate explicit method termed LCPFCT algorithm for the convection terms of the full Navier-Stokes equations. Thermoacoustic waves, which decay in large time due to the viscosity of fluid, are of the oscillatory nature. Much higher heat transfer rate can be achieved in an initial stage of transient processes, compared to that due to conduction. When a partial length of the cavity center line is heated, resulting thermoacoustic waves exhibit remarkable two-dimensional patterns.
Fractional Calculus in Hydrologic Modeling: A Numerical Perspective
David A. Benson; Mark M. Meerschaert; Jordan Revielle
2012-01-01
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.
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-01
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.
Numerical Estimation of Frictional Torques with Rate and State Friction
Arun K. Singh; T. N. Singh
2015-01-20
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.
On Numerical Considerations for Modeling Reactive Astrophysical Shocks
Papatheodore, Thomas L; Messer, Bronson
2014-01-01
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.
Non-Gaussian numerical errors versus mass hierarchy
Y. Meurice; M. B. Oktay
2000-05-12
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.
Numerical Simulation of Groundwater Withdrawal at the Nevada Test Site
Carroll, Rosemary; Giroux, Brian; Pohll, Greg; Hershey, Ronald; Russell, Charles; Howcroft, William
2004-01-28
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.
Numerical simulations of the decay of primordial magnetic turbulence
Kahniashvili, Tina [McWilliams Center for Cosmology and Department of Physics, Carnegie Mellon University, 5000 Forbes Ave, Pittsburgh, Pennsylvania 15213 (United States); Department of Physics, Laurentian University, Ramsey Lake Road, Sudbury, ON P3E 2C (Canada); Abastumani Astrophysical Observatory, Ilia State University, 2A Kazbegi Ave, Tbilisi, GE-0160 (Georgia); Brandenburg, Axel [Nordita, AlbaNova University Center, Roslagstullsbacken 23, 10691 Stockholm (Sweden); Department of Astronomy, Stockholm University, SE 10691 Stockholm (Sweden); Tevzadze, Alexander G. [Abastumani Astrophysical Observatory, Ilia State University, 2A Kazbegi Ave, Tbilisi, GE-0160 (Georgia); Faculty of Exact and Natural Sciences, Tbilisi State University, 1 Chavchavadze Avenue Tbilisi, GE-0128 (Georgia); Ratra, Bharat [Department of Physics, Kansas State University, 116 Cardwell Hall, Manhattan, Kansas 66506 (United States)
2010-06-15
We perform direct numerical simulations of forced and freely decaying 3D magnetohydrodynamic turbulence in order to model magnetic field evolution during cosmological phase transitions in the early Universe. Our approach assumes the existence of a magnetic field generated either by a process during inflation or shortly thereafter, or by bubble collisions during a phase transition. We show that the final configuration of the magnetic field depends on the initial conditions, while the velocity field is nearly independent of initial conditions.
Efficiency improvements for the numerical computation of NLO corrections
Sebastian Becker; Christian Reuschle; Stefan Weinzierl
2012-06-25
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 studies of galaxy formation using special purpose hardware
Matthias Steinmetz
2002-01-25
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.
Transient productivity index for numerical well test simulations
Blanc, G.; Ding, D.Y.; Ene, A.
1997-08-01
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.
Spectral Methods for Numerical Relativity. The Initial Data Problem
Lawrence E. Kidder; Lee Samuel Finn
1999-11-03
Numerical relativity has traditionally been pursued via finite differencing. Here we explore pseudospectral collocation (PSC) as an alternative to finite differencing, focusing particularly on the solution of the Hamiltonian constraint (an elliptic partial differential equation) for a black hole spacetime with angular momentum and for a black hole spacetime superposed with gravitational radiation. In PSC, an approximate solution, generally expressed as a sum over a set of orthogonal basis functions (e.g., Chebyshev polynomials), is substituted into the exact system of equations and the residual minimized. For systems with analytic solutions the approximate solutions converge upon the exact solution exponentially as the number of basis functions is increased. Consequently, PSC has a high computational efficiency: for solutions of even modest accuracy we find that PSC is substantially more efficient, as measured by either execution time or memory required, than finite differencing; furthermore, these savings increase rapidly with increasing accuracy. The solution provided by PSC is an analytic function given everywhere; consequently, no interpolation operators need to be defined to determine the function values at intermediate points and no special arrangements need to be made to evaluate the solution or its derivatives on the boundaries. Since the practice of numerical relativity by finite differencing has been, and continues to be, hampered by both high computational resource demands and the difficulty of formulating acceptable finite difference alternatives to the analytic boundary conditions, PSC should be further pursued as an alternative way of formulating the computational problem of finding numerical solutions to the field equations of general relativity.
Numerical implication of Riemann problem theory for fluid dynamics
Menikoff, R.
1988-01-01
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.
Numerical simulation of the hydrodynamical combustion to strange quark matter
Niebergal, Brian; Ouyed, Rachid; Jaikumar, Prashanth
2010-12-15
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 one dimension with neutrino emission from weak equilibrating reactions, and strange quark diffusion across the burning front. We also include entropy change from 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 {approx_equal}2 times saturation density). In a two-dimensional setting, such rapid speeds and neutrino cooling may allow for a flame wrinkle instability to develop, possibly leading to detonation.
Towards a novel wave-extraction method for numerical relativity
Andrea Nerozzi; Marco Bruni; Lior M. Burko; Virginia Re
2006-07-17
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.
Numerical simulation model for vertical flow in geothermal wells
Tachimori, M.
1982-01-01
A numerical simulation model for vertical flow in geothermal wells is presented. The model consists of equations for the conservation of mass, momentum, and energy, for thermodynamic state of water, for friction losses, for slip velocity relations, and of the criteria for various flow regimes. A new set of correlations and criteria is presented for two-phase flow to improve the accuracy of predictions; bubbly flow - Griffith and Wallis correlation, slug flow - Nicklin et al. one, annular-mist flow - Inoue and Aoki and modified by the author. The simulation method was verified by data from actual wells.
Advanced numerical methods in mesh generation and mesh adaptation
Lipnikov, Konstantine; Danilov, A; Vassilevski, Y; Agonzal, A
2010-01-01
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-15
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.
Connecting Numerical Relativity and Data Analysis of Gravitational Wave Detectors
Deirdre Shoemaker; Karan Jani; Lionel London; Larne Pekowsky
2015-03-09
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.
High numerical aperture projection system for extreme ultraviolet projection lithography
Hudyma, Russell M. (San Ramon, CA)
2000-01-01
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.
Numerical routines for predicting ignition in pyrotechnic devices
Pierce, K.G.
1986-06-01
Two numerical models of the thermal processes leading to ignition in a pyrotechnic device have been developed. These models are based on finite difference approximations to the heat diffusion equation, with temperature-dependent thermal properties, in a single spatial coordinate. The derivation of the finite difference equations is discussed and the methods employed at boundaries and interfaces are given. The sources of the thermal-properties data are identified and how these data are used is explained. The program structure is explained and example runs of the programs are given.
Diffusive mesh relaxation in ALE finite element numerical simulations
Dube, E.I.
1996-06-01
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.
Worldline Numerics for Energy-Momentum Tensors in Casimir Geometries
Schafer, Marco; Gies, Holger
2015-01-01
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.
Numerical model of Electron Cyclotron Resonance Ion Source
Mironov, V; Bondarchenko, A; Efremov, A; Loginov, V
2015-01-01
Important features of Electron Cyclotron Resonance Ion Source (ECRIS) operation are accurately reproduced with a numerical code. The code uses the particle-in-cell technique to model a dynamics of ions in ECRIS plasma. It is shown that gas dynamical ion confinement mechanism is sufficient to provide the ion production rates in ECRIS close to the experimentally observed values. Extracted ion currents are calculated and compared to the experiment for few sources. Changes in the extracted ion currents are obtained with varying the gas flow into the source chamber and the microwave power. Empirical scaling laws for ECRIS design are studied and the underlying physical effects are discussed.
Numerical simulation of carbon arc discharge for nanoparticle synthesis
Kundrapu, M.; Keidar, M.
2012-07-15
Arc discharge with catalyst-filled carbon anode in helium background was used for the synthesis of carbon nanoparticles. In this paper, we present the results of numerical simulation of carbon arc discharges with arc current varying from 10 A to 100 A in a background gas pressure of 68 kPa. Anode sublimation rate and current voltage characteristics are compared with experiments. Distribution of temperature and species density, which is important for the estimation of the growth of nanoparticles, is obtained. The probable location of nanoparticle growth region is identified based on the temperature range for the formation of catalyst clusters.
Numerical analysis of impact processes of granular jets
Tomohiko G. Sano; Hisao Hayakawa
2013-03-20
The rheology of a three-dimensional granular jet during an impact is investigated numerically. The cone-like scattering pattern and the sheet-like pattern observed in an experiment [X. Cheng, et al. Phys. Rev. Lett. 99, 188001 (2007)] can be reproduced through our calculation. We discuss the constitutive equation for granular jet impact in terms of our simulation. From the analysis of an effective friction constant, which is the ratio between the shear stress and the pressure the assumption of the zero yield stress would be natural in our setup and the shear visocity is not small in contrast to the suggestion by the experiment.
Property:Iso3166Numeric | Open Energy Information
AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on QA:QAsource History ViewMayo,AltFuelVehicle2 Jump to: navigation, searchContDiv JumpTechDsc Jump to:consIsDOELabNumeric
Sandia Energy - Numerical Manufacturing And Design Tool (NuMAD)
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: Alternative Fuels Data Center Homesum_a_epg0_fpd_mmcf_m.xls" ,"Available from WebQuantity ofkandz-cm11 Outreach Home RoomPreservation of Fe(II)Geothermal Energy &WaterNew CREW DatabaseNuclear FuelsNuclearNumerical
Numerical Model for Conduction-Cooled Current Lead Heat Loads
White, M.J.; Wang, X.L.; Brueck, H.D.; /DESY
2011-06-10
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).
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-15
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 simulation of multi-layered textile composite reinforcement forming
Wang, P.; Hamila, N.; Boisse, P.
2011-05-04
One important perspective in aeronautics is to produce large, thick or/and complex structural composite parts. The forming stage presents an important role during the whole manufacturing process, especially for LCM processes (Liquid Composites Moulding) or CFRTP (Continuous Fibre Reinforcements and Thermoplastic resin). Numerical simulations corresponding to multi-layered composite forming allow the prediction for a successful process to produce the thick parts, and importantly, the positions of the fibres after forming to be known. This paper details a set of simulation examples carried out by using a semi-discrete shell finite element made up of unit woven cells. The internal virtual work is applied on all woven cells of the element taking into account tensions, in-plane shear and bending effects. As one key problem, the contact behaviours of tool/ply and ply/ply are described in the numerical model. The simulation results not only improve our understanding of the multi-layered composite forming process but also point out the importance of the fibre orientation and inter-ply friction during formability.
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-01
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.
New insights on numerical error in symplectic integration
Hugo Jiménez-Pérez; Jean-Pierre Vilotte; Barbara Romanowicz
2015-08-13
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.
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 QA:QA J-E-1 SECTION JEnvironmental Jump to:EAandAmminexInformation Bartholomew County,Creek, Texas:Bartow
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-01
transport and numerical weather modeling. J. Applied cross correlations. Weather and Forecasting, 8:4, 401?of radiation for numerical weather prediction and climate
Pearce, Fred
2003-01-01
We use a 3-D finite difference numerical model to generate synthetic seismograms from a simple fractured reservoir
PHYSICAL AND NUMERICAL SIMULATIONS OF SUBSIDENCE ABOVE HIGH EXTRACTION COAL MINES
Sutherland, Herbert; Heckes, Albert; Taylor, Lee
1984-01-01
and the formation of a subsidence trough. The data fromNUMERICAL SIMULATIONS OF SUBSIDENCE ABOVE HIGH EXTRACTION
Shervais, John W.
drill core show that the parent magmas of these lavas evolved by crystal fractionation at shal- low to intermediate crustal depths, punctuated by episodic recharge with more primitive compositions and assimilation a 1136-m-thick strati- graphic section of basalt sampled by scientific drilling that preserves
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
Numerical Methods for the Bogoliubov-Tolmachev-Shirkov model in superconductivity theory
Zhihao Ge; Ruihua Li
2015-03-08
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.
Functional Imaging of Numerical Processing in Adults and 4-y-Old Children
Makous, Walter
Functional Imaging of Numerical Processing in Adults and 4-y-Old Children Jessica F. Cantlon1 to numerical deviants similarly in 4-y-old children and adults. To our knowledge, this is the first evidence of numerical processing in adults and 4-y-old children. PLoS Biol 4(5): e125. DOI: 10.1371/ journal
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 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
NUMERICAL SOLUTION OF RESERVOIR FLOW MODELS BASED ON LARGE TIME STEP OPERATOR SPLITTING ALGORITHMS
NUMERICAL SOLUTION OF RESERVOIR FLOW MODELS BASED ON LARGE TIME STEP OPERATOR SPLITTING ALGORITHMS. Special focus is posed on the numerical solution algorithms for the saturation equation, which. The general background for the reservoir ow model is reviewed, and the main features of the numerical
NUMERICAL SOLUTION OF RESERVOIR FLOW MODELS BASED ON LARGE TIME STEP OPERATOR SPLITTING ALGORITHMS
NUMERICAL SOLUTION OF RESERVOIR FLOW MODELS BASED ON LARGE TIME STEP OPERATOR SPLITTING ALGORITHMS focus is posed on the numerical solution algorithms for the saturation equation, which is a convectionÂ eral background for the reservoir flow model is reviewed, and the main features of the numerical
A MultiLevel Preconditioner with Applications to the Numerical Simulation of Coating Problems
Zhang, Jun
A MultiLevel Preconditioner with Applications to the Numerical Simulation of Coating Problems of unstructured sparse linear systems arising from the numerical simulation of coating problems. The coef ficient unstructured sparse linear systems from the numerical simulation of coating problems. Coating is a delicate
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
Mathematical and Numerical Analyses of Peridynamics for Multiscale Materials Modeling
Du, Qiang
2014-11-12
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 studyof 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.
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-01
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 determination of the material properties of porous dust cakes
D. Paszun; C. Dominik
2008-02-13
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.
Numerical code SELFAS-3 and electrodynamic aggregation of magnetized nanodust
Kukushkin, A B
2010-01-01
The principles of the parallel numerical code SELFAS-3 are presented. The code modifies previous version of the code to enable parallel computations of electrodynamic aggregation in a many-body system of basic blocks which are taken as strongly magnetized thin rods (i.e., one-dimensional static magnetic dipoles), with electric conductivity and static electric charge, screened with its own static plasma sheath. The aggregation modelling includes the electric current dynamics in a complicated evolving network to describe the processes of external and internal electric short-circuiting. The code enables the continuous modelling of a transition between the following states: randomly situated ensemble of solitary basic blocks; electric current-carrying filamentary system; restructured filamentary network with a trend towards a fractal skeletal structuring. The latter trend is illustrated with generation of a bigger magnetic dipole in (i) homogeneous random ensemble between the biased electrodes in the presence of ...
Laser spark plug numerical design process with experimental validation
McIntyre, D.; Woodruff, S.
2011-01-01
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.
Design and numerical simulation of thermionic electron gun
Hosseinzadeh, M
2015-01-01
This paper reports the simulation of an electron gun. The effects of some parameters on the beam quality were studied and optimal choices were identified. It gives numerical beam qualities in common electrostatic triode gun, and the dependences on design parameters such as electrode geometries and bias voltages to these electrodes are shown. An electron beam of diameter 5 mm with energy of five kilo electron volt was assumed for simulation process. Some design parameters were identified as variable parameters in the presence of space charge. These parameters are the inclination angle of emission electrode, the applied voltage to focusing electrode, the gap width between the emission electrode and the focusing electrode and the diameter of the focusing electrode. The triode extraction system is designed and optimized by using CST software (for Particle Beam Simulations). The physical design of the extraction system is given in this paper. From the simulation results, it is concluded that the inclination angle ...
A numerical model of aerosol scavenging: Part 1, Microphysics parameterization
Molenkamp, C.R.; Bradley, M.M.
1991-09-01
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.
MHD Remote Numerical Simulations: Evolution of Coronal Mass Ejections
L. Hernandez-Cervantes; A. Santillan; A. R. Gonzalez-Ponce
2008-12-22
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.
Electromagnetic Pulse Propagation over Nonuniform Earth Surface: Numerical Simulation
Alexei V. Popov; Vladimir V. Kopeikin
2007-04-14
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.
A fast direct numerical simulation method for characterising hydraulic roughness
Chung, Daniel; MacDonald, Michael; Hutchins, Nicholas; Ooi, Andrew
2015-01-01
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...
Progress report on LBL's numerical modeling studies on Cerro Prieto
Halfman-Dooley, S.E.; Lippman, M.J.; Bodvarsson, G.S.
1989-04-01
An exploitation model of the Cerro Prieto geothermal system is needed to assess the energy capacity of the field, estimate its productive lifetime and develop an optimal reservoir management plan. The model must consider the natural state (i.e., pre-exploitation) conditions of the system and be able to predict changes in the reservoir thermodynamic conditions (and fluid chemistry) in response to fluid production (and injection). This paper discusses the results of a three-dimensional numerical simulation of the natural state conditions of the Cerro Prieto field and compares computed and observed pressure and temperature/enthalpy changes for the 1973--1987 production period. 16 refs., 24 figs., 2 tabs.
Three dimensional numerical simulations of the UPS-292-SC engine
O'Rourke, P.J.; Amsden, A.A.
1987-01-01
We present and analyze three-dimensional calculations of the spray, mixing and combustion in the UPS-292 stratified charge engine for three different operating conditions, corresponding to overall air-fuel ratios between 22.4 and 61.0. The numerical calculations are performed with KIVA, a multidimensional arbitrary-mesh, finite-difference hydrodynamics program for internal combustion engine applications. The calculations use a mesh of 10,000 computational cells, which conform to the shape of the piston bowl and cylinder and move to follow piston motion. Each operating condition is calculated from intake valve closure at 118/sup 0/ BTDC to 90/sup 0/ ATDC and requires approximately three hours of CRAY-XMP computer time.
Numerical Weather Forecasting at the Savannah River Site
Buckley, R.L. [Westinghouse Savannah River Company, AIKEN, SC (United States)
1998-11-01
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.
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-01
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.
Inertial range scaling in numerical turbulence with hyperviscosity
Nils Erland L. Haugen; Axel Brandenburg
2004-10-26
Numerical turbulence with hyperviscosity is studied and compared with direct simulations using ordinary viscosity and data from wind tunnel experiments. It is shown that the inertial range scaling is similar in all three cases. Furthermore, the bottleneck effect is approximately equally broad (about one order of magnitude) in these cases and only its height is increased in the hyperviscous case--presumably as a consequence of the steeper decent of the spectrum in the hyperviscous subrange. The mean normalized dissipation rate is found to be in agreement with both wind tunnel experiments and direct simulations. The structure function exponents agree with the She-Leveque model. Decaying turbulence with hyperviscosity still gives the usual t^{-1.25} decay law for the kinetic energy, and also the bottleneck effect is still present and about equally strong.
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-15
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.
Monitoring Erebus volcano's active lava lake; tools, techniques and observations
Peters, Nial John
2015-03-03
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... and wind generators. Instead, power is generated 0.5 km down-slope at the Nausea Knob (NKB) seismic station site, where a 1000 Ah battery bank is charged using a ? 0.5 kW array of photovoltaic panels, and a 100 W wind turbine. A Schaefer AEP-1500 inverter...
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...
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
Wang, Zhi Jian "ZJ"
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2002; 39 consuming if `non-water-tight' geometries are given. Signiÿcant user time can be taken to repair or clean cleaning or repairing, dramatically reducing grid generation time. An I2B adaptive Cartesian grid
Guzina, Bojan
INTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS Int. J. Numer. Anal, University of Minnesota, 500 Pillsbury Drive S.E., Minneapolis, MN 55455, U.S.A. SUMMARY A mathematical of Minnesota, 500 Pillsbury Drive S.E., Minneapolis, MN 55455, U.S.A. } Graduate Research Assistant. Contract
Numerical simulation of linear fiction welding (LFW) processes
Fratini, L.; La Spisa, D. [University of Palermo-Dept. of Industrial engineering (Italy)
2011-05-04
Solid state welding processes are becoming increasingly important due to a large number of advantages related to joining ''unweldable'' materials and in particular light weight alloys. Linear friction welding (LFW) has been used successfully to bond non-axisymmetric components of a range of materials including titanium alloys, steels, aluminum alloys, nickel, copper, and also dissimilar material combinations. The technique is useful in the research of quality of the joints and in reducing costs of components and parts of the aeronautic and automotive industries.LFW involves parts to be welded through the relative reciprocating motion of two components under an axial force. In such process the heat source is given by the frictional forces work decaying into heat determining a local softening of the material and proper bonding conditions due to both the temperature increase and the local pressure of the two edges to be welded. This paper is a comparative test between the numerical model in two dimensions, i.e. in plane strain conditions, and in three dimensions of a LFW process of AISI1045 steel specimens. It must be observed that the 3D model assures a faithful simulation of the actual threedimensional material flow, even if the two-dimensional simulation computational times are very short, a few hours instead of several ones as the 3D model. The obtained results were compared with experimental values found out in the scientific literature.
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-01
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
Probing Modified Gravity with Atom-Interferometry: a Numerical Approach
Schlogel, Sandrine; Fuzfa, Andre
2015-01-01
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...
Numerical tests of dynamical friction in gravitational inhomogeneous systems
A. Del Popolo
2003-05-05
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 modeling of the wind flow over a transverse dune
Ascânio D. Araújo; Eric J. R. Parteli; Thorsten Poeschel; José S. Andrade Jr.; Hans J. Herrmann
2013-09-30
Transverse dunes, which form under unidirectional winds and have fixed profile in the direction perpendicular to the wind, occur on all celestial objects of our solar system where dunes have been detected. Here we perform a numerical study of the average turbulent wind flow over a transverse dune by means of computational fluid dynamics simulations. We find that the length of the zone of recirculating flow at the dune lee --- the {\\em{separation bubble}} --- displays a surprisingly strong dependence on the wind shear velocity, $u_{\\ast}$: it is nearly independent of $u_{\\ast}$ for shear velocities within the range between $0.2\\,$m$$s and $0.8\\,$m$$s but increases linearly with $u_{\\ast}$ for larger shear velocities. Our calculations show that transport in the direction opposite to dune migration within the separation bubble can be sustained if $u_{\\ast}$ is larger than approximately $0.39\\,$m$$s, whereas a larger value of $u_{\\ast}$ (about $0.49\\,$m$$s) is required to initiate this reverse transport.
Numerical modeling of the wind flow over a transverse dune
Araújo, Ascânio D; Poeschel, Thorsten; Andrade, José S; Herrmann, Hans J
2013-01-01
Transverse dunes, which form under unidirectional winds and have fixed profile in the direction perpendicular to the wind, occur on all celestial objects of our solar system where dunes have been detected. Here we perform a numerical study of the average turbulent wind flow over a transverse dune by means of computational fluid dynamics simulations. We find that the length of the zone of recirculating flow at the dune lee --- the {\\em{separation bubble}} --- displays a surprisingly strong dependence on the wind shear velocity, $u_{\\ast}$: it is nearly independent of $u_{\\ast}$ for shear velocities within the range between $0.2\\,$m$$s and $0.8\\,$m$$s but increases linearly with $u_{\\ast}$ for larger shear velocities. Our calculations show that transport in the direction opposite to dune migration within the separation bubble can be sustained if $u_{\\ast}$ is larger than approximately $0.39\\,$m$$s, whereas a larger value of $u_{\\ast}$ (about $0.49\\,$m$$s) is required to initiate this reverse transport.
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-01
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.
Numeric modeling of fire suppression by organophosphorous inhibitors
Makhviladze, G M; Zykov, A P
2008-01-01
Numerical calculations of the effect of organophosphorous inhibitor (CF3CH2O)3P and its mixtures with carbon dioxide on propane flames are carried out using the three dimensional Reynolds-averaged Navier-Stokes (RANS) equations in the low Mach number approximation. The k-e model of turbulence, the EDC combustion model and the weighted-sum-of-gray-gases model of radiation are used. The Westbrook global-kinetic scheme with fractional order of reaction was used for the calculation of chemical reaction rate of propane combustion. The empirical expression for the correction factor for the chemical reaction rate was used to model the effect of organophosphorous inhibitor no the reaction. Two series of test calculations for different values of the correction factor are carried out. Dependences of the minimum extinguishing concentration of the inhibitor per carbon dioxide volume concentration in the extinguishing mixtures were obtained. The results of test calculations are shown to agree reasonably with the experimen...
Comparison of numerical models of a pyrotechnic device
Pierce, K.G.
1986-01-01
The predictions of two numerical models of a hot-wire initiated pyrotechnic device are compared to each other and to experimental results. Both models employ finite difference approximations to the heat diffusion equation in cylindrical coordinates. The temperature dependence of the thermal properties of the pyrotechnic materials and of the bridgewire are modeled. An Arrhenius' model is used to describe the exothermic reaction in the powder. One model employs a single radial coordinate and predicts the radial temperature distribution in the bridgewire and surrounding powder mass. In addition to the radial coordinate, the other model also employs a longitudinal coordinate to predict the temperature distribution parallel to the axis of the bridgewire. The predictions of the two-dimensional model concerning the energy requirements for ignition and the energy losses from the ends of the bridgewire are presented. A comparison of the predictions of the two models and the development of thermal gradients are employed to define the regime where the assumption, in the one-dimensional model, of negligible heat transfer axial to the bridgewire does not lead to significant error. The general problems associated with predicting ignition from a diffusion model are also discussed.
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-01
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.
Testing numerical relativity with the shifted gauge wave
Maria C. Babiuc; Bela Szilagyi; Jeffrey Winicour
2006-02-17
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.
Numerical study of energy diffusion in King models
Tom Theuns
1995-11-07
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-31
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.
Numerical simulation of plasma heating of a composite powder particle
Demetriou, M.D.; Lavine, A.S.; Ghoniem, N.M.
1999-07-01
The use of fine composite powder particles (composed of a ceramic core and a metallic coating) in plasma spraying processes is desirable in developing thin film coatings that possess high abrasion as well as high fracture resistance. Quantitative knowledge of the thermal behavior of a composite particle in a plasma beam is essential in optimizing the process variables to achieve uniform melting of the coating material. In this work, a numerical model is developed to analyze the in-flight thermal behavior of a spherically symmetric WC-Co composite particle travelling in an argon arc-jet DC plasma under strongly unsteady plasma conditions. The model gives quantitative as well as qualitative information about the thermal response of the heated particle. The important features that are addressed are the temperature response of the particle; the history of the location of the melting and vaporization fronts; and the physical state of the particle at the end of its flight. For the conditions investigated, it was determined that the internal conduction resistance is negligible as compared to the net external resistance. However, the presence of the ceramic base was found to affect the transient heating process since its content in the particle composition determines the time constant of the process. Another interesting observation is that proper selection of the particle injection speed and injection location can be effective means for optimizing the heating process and achieving uniform melting of the coating material.
Numerical Relativity in Spherical Polar Coordinates: Off-center Simulations
Thomas W. Baumgarte; Pedro J. Montero; Ewald Müller
2015-06-03
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.
Probing Modified Gravity with Atom-Interferometry: a Numerical Approach
Sandrine Schlogel; Sebastien Clesse; Andre Fuzfa
2015-07-11
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.
Numerical investigation of spontaneous flame propagation under RCCI conditions
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Bhagatwala, Ankit V; Sankaran, Ramanan; Kokjohn, Sage; Chen, Jacqueline H
2015-06-30
This paper presents results from one and two-dimensional direct numerical simulations under Reactivity Controlled Compression Ignition (RCCI) conditions of a primary reference fuel (PRF) mixture consisting of n-heptane and iso-octane. RCCI uses in-cylinder blending of two fuels with different autoignition characteristics to control combustion phasing and the rate of heat release. These simulations employ an improved model of compression heating through mass source/sink terms developed in a previous work by Bhagatwala et al. (2014), which incorporates feedback from the flow to follow a predetermined experimental pressure trace. Two-dimensional simulations explored parametric variations with respect to temperature stratification, pressure profiles andmore »n-heptane concentration. Furthermore, statistics derived from analysis of diffusion/reaction balances locally normal to the flame surface were used to elucidate combustion characteristics for the different cases. Both deflagration and spontaneous ignition fronts were observed to co-exist, however it was found that higher n-heptane concentration provided a greater degree of flame propagation, whereas lower n-heptane concentration (higher fraction of iso-octane) resulted in more spontaneous ignition fronts. A significant finding was that simulations initialized with a uniform initial temperature and a stratified n-heptane concentration field, resulted in a large fraction of combustion occurring through flame propagation. The proportion of spontaneous ignition fronts increased at higher pressures due to shorter ignition delay when other factors were held constant. For the same pressure and fuel concentration, the contribution of flame propagation to the overall combustion was found to depend on the level of thermal stratification, with higher initial temperature gradients resulting in more deflagration and lower gradients generating more ignition fronts. Statistics of ignition delay are computed to assess the Zel’dovich (1980) theory for the mode of combustion propagation based on ignition delay gradients.« less
Mathiesen, Patrick; Kleissl, Jan
2011-01-01
Numerical Weather Prediction (NWP), Solar Forecasting 1. to more accurate prediction of solar irradiance, given a to create daily solar electricity predictions accurate to
X. CORRELATIONS BETWEEN DIFFUSE AND GLOBAL IRRADIANCE For numerous applications, particularly those
Oregon, University of
33 X. CORRELATIONS BETWEEN DIFFUSE AND GLOBAL IRRADIANCE For numerous applications, particularly these components, values for direct beam and diffuse irradiance are usu- ally unavailable. However, considerable in
Mass (gas and water) transfers within concrete are an important issue dealing with numerous appli, the latter aims at extending its nuclear powerplants lifetime which requires accurate quantifications
Geotextiles and Geomembranes 11 (1992) 115-131 Numerical Investigation of Controlled Yielding of
Horvath, John S.
1992-01-01
Geotextiles and Geomembranes 11 (1992) 115-131 Numerical Investigation of Controlled Yielding- constructed to active conditions is dependent on the cohesion and 115 Geotextiles and Geomembranes 0266
Mathiesen, Patrick James
2013-01-01
Evaluation of numerical weather prediction solar irradiancecycle: The RUC. Monthly Weather Review, 132 (2), 495-518.representations. Monthly Weather Review, 139 (6), 1972-1995.
NUMERICAL STUDY OF A TURBULENT HYDRAULIC JUMP Qun Zhao 1 Shubhra K. Misra1
Zhao, Qun
. Hydraulic jumps are commonly used as energy dissipators and they have been studied intensively by hydraulicNUMERICAL STUDY OF A TURBULENT HYDRAULIC JUMP Qun Zhao 1 Shubhra K. Misra1 Ib A. Svendsen 1 (Member of a turbulent hydraulic jump. The numerical model is based on RIPPLE (Kothe et al., 1994) with two turbulence
The VOLNA code for the numerical modelling of tsunami waves: generation, propagation and
of the finite volume scheme implemented in the code. We explain the numerical treatment of the wet is decided on the base of inundation maps which are produced with this type of numerical tools. Finally we and the perspectives for future research presented. Key words: tsunami waves, shallow water equations, tsunami
Numerical aperture influence on 3-D multi-layer optical data storage systems , Edwin P. Walkera
Esener, Sadik C.
Numerical aperture influence on 3-D multi-layer optical data storage systems Yi Zhanga* , Edwin P storage system is analyzed. Keywords: NA, multi-layer data storage, two-photon recording, capacity) 550-0596, Fax: (858) 550-0917 #12;Numerical aperture influence on 3-D multi-layer optical data storage
Fast Fourier Analysis for SL 2 over a Finite Field and Related Numerical Experiments
Fast Fourier Analysis for SL 2 over a Finite Field and Related Numerical Experiments John D of Fourier Analysis 4. Implementing the Computation 5. Discussion of Numerical Results 6. Speculations and Open Problems 7. Appendix: Fourier Inversion and Convolution for SL 2 Acknowledgements References AMS
Fast Fourier Analysis for SL2 over a Finite Field and Related Numerical Experiments
Fast Fourier Analysis for SL2 over a Finite Field and Related Numerical Experiments John D of Fourier Analysis 4. Implementing the Computation 5. Discussion of Numerical Results 6. Speculations and Open Problems 7. Appendix: Fourier Inversion and Convolution for SL2 Acknowledgements References AMS
Adequate Numerical Solution of Air Pollution Problems by positive Difference Schemes on
Ehrhardt, Matthias
Adequate Numerical Solution of Air Pollution Problems by positive Difference Schemes on unbounded In this work we deal with the numerical solution of some problems of air pollution. Since the problems) pollution. The efficiency and accuracy of our method is illustrated by an example. Key words: air pollution
Adequate Numerical Solution of Air Pollution Problems by positive Di#erence Schemes on
Ehrhardt, Matthias
Adequate Numerical Solution of Air Pollution Problems by positive Di#erence Schemes on unbounded In this work we deal with the numerical solution of some problems of air pollution. Since the problems) pollution. The e#ciency and accuracy of our method is illustrated by an example. Key words: air pollution
On the numerical relaxation of single-slip plasticity in finite strains
Moroz, Vitaly
;[2] Conti S., Relaxation of single-slip single-crystal plasticity with linear self-hardening, In: Proceed., Orlando A., Mixed analytical-numerical relaxation in single-slip crystal plasticity, AcceptedOn the numerical relaxation of single-slip plasticity in finite strains Antonio Orlando (Swansea
Merging of coronal and heliospheric numerical two-dimensional D. Odstrcil,1,5
California at Berkeley, University of
spatial and temporal scales between the Sun and Earth. Specialized physically based numerical models have been developed to address particular aspects of the entire system. However, an integrated modeling physically based numerical models have been developed to address particular aspects of the entire system
Efficient numerical methods for the instationary solution of laminar reacting gas
Vuik, Kees
Efficient numerical methods for the instationary solution of laminar reacting gas flow problems Universiteit Delft Efficient numerical methods for the instationary solution of laminar reacting gas flow of laminar reacting gas flow problems Sander van Veldhuizen Production processes of high-purity, high
The effective viscosity of rocksalt: implementation of steady-state creep laws in numerical
van Keken, Peter
of salt tectonics have received renewed attention in the last few years. This has been stimulated in numerical models of salt diapirism P.E. van Keken, C.J. Spiers, A.P. van den Berg and E.J. Muyzert and fluid-enhanced grain-boundary diffusion creep, has been used in numerical models of salt diapirism
Numerical Simulation of Pulse-Tube Refrigerators: 1D model I.A. Lyulina1
Eindhoven, Technische Universiteit
of a piston, an aftercooler (AC), a regenerator, a cold heat exchanger (CHX), a tube, a hot heat exchanger numerical model has been introduced to study steady oscillatory heat and mass transfer in the tube section, numerical simulation, high resolution scheme 1 Introduction The pulse tube is a relatively new type
A Flux-Limited Numerical Method for the MHD Equations to Simulate Propulsive Plasma Flows
Choueiri, Edgar
to be effective tools in plasma propulsion research, a higher order accu- rate solver that captures MHD shocks approach, numerical simulations are valuable tools in plasma thruster research. More- over, simulations can Simula- tions The importance of numerical simulation in advancing plasma thruster research was realized
Numerical tools applied to power reactor noise analysis Christophe Demazie`re*, Imre Pazsit
Demazière, Christophe
Review Numerical tools applied to power reactor noise analysis Christophe Demazie`re*, Imre Pa and application of the numerical tools employed. The code that was developed yields the space and non-critical systems with an external source. Some appli- cations of these tools to power reactor
Young, William R.
of mechanical energy required to mix the global ocean (Garrett and St. Laurent 2002). Tidal forces perform aboutNumerical and Analytical Estimates of M2 Tidal Conversion at Steep Oceanic Ridges EMANUELE DI) ABSTRACT Numerical calculations of the rate at which energy is converted from the external to internal
Numerical study on optimal Stirling engine regenerator matrix designs taking into account matrix design that improves the efficiency of a Stirling engine has been developed in a numerical study of the existing SM5 Stirling engine. A new, detailed, one-dimensional Stirling engine model that delivers results
Numerical and experimental investigations of the effective thermal conductivity of snow
Numerical and experimental investigations of the effective thermal conductivity of snow N. Calonne numerical simulations of the conductiv- ity of snow using microtomographic images. The full tensor of the effective thermal conductivity (keff) was computed from 30 threedimensional images of the snow microstruc
A regional numerical ocean model of the circulation in the Bay of Biscay
Drijfhout, Sybren
A regional numerical ocean model of the circulation in the Bay of Biscay Y. Friocourt,1,2,3 B Peninsula and in the Bay of Biscay is investigated by means of a regional ocean model. In particular numerical ocean model of the circulation in the Bay of Biscay, J. Geophys. Res., 112, C09008, doi:10
ENINALS Experimental and Numerical Investigations of Nonlinearity in soils using Advanced of Nonlinearity in Soils Using Advanced Laboratory-scaled Models Concluding Workshop Ispra, May 28-30 2013 IFSTTAR-geological constraints - centrifuge modelling - numerical modelling MILESTONES - Experiencing techniques for realising
van Leeuwen, Peter Jan
Lagrangian validation of numerical drifter trajectories using drifting buoys: Application validation KolmogorovSmirnov statistics Drifting buoys a b s t r a c t The skill of numerical Lagrangian to the trajectories of drifting buoys in the real ocean. The skill assess- ment is performed using the two
Numerical Modeling of CIGS Solar Cells: Definition of the Baseline and
Sites, James R.
Thesis Numerical Modeling of CIGS Solar Cells: Definition of the Baseline and Explanation our supervision by Markus Gloeckler entitled "Numerical Modeling of CIGS Solar Cells: Definition. A three-layer structure, simulating a Cu(InGa)Se2 (CIGS) heterojunction solar cell, was set up using
NUMERICAL MODELING OF CIGS AND CdTe SOLAR CELLS: SETTING THE BASELINE
Sites, James R.
NUMERICAL MODELING OF CIGS AND CdTe SOLAR CELLS: SETTING THE BASELINE M. Gloeckler, A important complications that are often found in experimental CIGS and CdTe solar cells. 1. INTRODUCTION Numerical modeling of polycrystalline thin-film solar cells is an important strategy to test the viability
Physical mechanisms involved in grooved flat heat pipes: experimental and numerical analyses
Paris-Sud XI, Université de
1 Physical mechanisms involved in grooved flat heat pipes: experimental and numerical analyses.lefevre@insa-lyon.fr) Abstract An experimental database, obtained with flat plate heat pipes (FPHP) with longitudinal grooves to results of numerical models. Physical mechanisms involved in grooved heat pipes are discussed, including
Using the Context-enhanced Additive Heuristic for Temporal and Numeric Planning
Vetter, Thomas
, the main dif- ference lies in the heuristic: Sapa utilizes a temporal version of planning graphs firstUsing the Context-enhanced Additive Heuristic for Temporal and Numeric Planning Patrick Eyerich,mattmuel,roeger}@informatik.uni-freiburg.de Abstract Planning systems for real-world applications need the abil- ity to handle concurrency and numeric
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
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
Circular Harmonic Decomposition Approach for Numerical Inversion of Circular Radon Transforms
Louis, Alfred K.
Circular Harmonic Decomposition Approach for Numerical Inversion of Circular Radon Transforms Gaël@num.uni-sb.de ABSTRACT Numerical inversions via circular harmonic decomposition for two classes of circular Radon transforms are established. The first class deals with the Radon transform (RT) de- fined on circular arcs
Numerical modeling of extreme rogue waves generated by directional energy focusing
Grilli, Stéphan T.
Numerical modeling of extreme rogue waves generated by directional energy focusing Christophe that contributes to the generation of extreme waves, also known as rogue waves, in the ocean. To simulate and analyze this phenomenon, we generate extreme waves in a 3D numerical wave tank (NWT), by specifying
Numerical modeling of alongshore sediment transport and shoreline change along the Galveston coast
Sitanggang, Khairil Irfan
2005-02-17
coast for the period of 2001-2011. The first assessment shows that GENESIS has a numerical error which tends to be large for low energy wave (small breaking wave height) and large breaking wave angle. This numerical inaccuracy cannot be neglected...
A Simulation and Decision Framework for Selection of Numerical Solvers in Scientific Computing
Burns, Peter
A Simulation and Decision Framework for Selection of Numerical Solvers in Scientific Computing by the journal ACM Transactions On Mathematical Software (ACM 2004 [1]) and the Numerical Recipes published Peter Bunus Department of Computer and Information Science Linköping University, Sweden petbu
ON THE NUMERICAL RANGE OF A MATRIX TRANSLATED FROM THE GERMAN BY PAUL F. ZACHLIN
Eindhoven, Technische Universiteit
-7058, USA paul.zachlin@case.edu. Department of Mathematics and Computing Science, Eindhoven UniversityON THE NUMERICAL RANGE OF A MATRIX TRANSLATED FROM THE GERMAN BY PAUL F. ZACHLIN AND MICHIEL E by mathematicians who work in the area of numerical ranges, thus it is hoped that this translation may be useful
Formal Numerical Program Analysis William M. Farmer and F. Javier Thayer
Farmer, William M.
to produce a machine-checked proof that the program computes an abstractly specified mathematical function computations involving numerical datatypes--for example, programs for computing solutions of all sorts, a precise mathematical analysis of even a simple numerical program is very difficult for the following
Huerta, Antonio
NUMERICAL MODELLING OF IMPACT NOISE Cristina Díaz Cereceda | Laboratori de Càlcul Numèric Modelling of Impact Noise. In Congreso de Métodos Numéricos en Ingeniería 2009. [2] J. Hetherington, J. Poblet-Puig, A. Rodríguez-Ferran. Numerical Modelling of Impact Noise in Lightweight Floors: a Modal
Stutzmann, Eléonore
Detection of microseismic compressional (P) body waves aided by numerical modeling of oceanic noise), Detection of microseismic compressional (P) body waves aided by numerical modeling of oceanic noise sources in seismic noise, body waves present appealing properties but are still challenging to extract. Here we first
Numerical simulation of flow and mixing behavior of solids on a moving grate combustion system
Columbia University
Numerical simulation of flow and mixing behavior of solids on a moving grate combustion system by #12;ii Numerical simulation of flow and mixing behavior of solids on a moving grate combustion system, and to a large extent influences the combustion process. Municipal solid waste (MSW) is not a uniform fuel
A Personal Perspective on the History of the Numerical Analysis of Fredholm Integral
Atkinson, Kendall
A Personal Perspective on the History of the Numerical Analysis of Fredholm Integral Equations of the Second Kind Kendall Atkinson The University of Iowa July 25, 2008 Abstract This is a personal perspective on the development of numerical meth- ods for solving Fredholm integral equations of the second kind, discussing work
Numerical Analysis of a one dimensional Diffusion Equation for a single chamber Microbial Fuel Cell generation within our society. Microbial fuel cells (MFCs) represent a new form of renewable energy using a Linked Simulation Optimization (LSO) technique E521: Advanced Numerical Methods Eric A. Zielke
CONVERGENCE OF NUMERICAL APPROXIMATIONS TO A PHASE FIELD BENDING ELASTICITY MODEL OF
CONVERGENCE OF NUMERICAL APPROXIMATIONS TO A PHASE FIELD BENDING ELASTICITY MODEL OF MEMBRANE, Number 1, Pages 1Â18 CONVERGENCE OF NUMERICAL APPROXIMATIONS TO A PHASE FIELD BENDING ELASTICITY MODEL OF MEMBRANE DEFORMATIONS QIANG DU AND XIAOQIANG WANG This paper is dedicated to Prof.Max Gunzburger
Zhu, Chen
of radioactive waste repositories. A case study of a regional aquifer in northeastern Arizona shows the rechargeEstimate of recharge from radiocarbon dating of groundwater and numerical flow and transport ages using a linked numerical 14 C transport and flow model while hydraulic conductivity values
OBJECT ORIENTED PROGRAMMING TECHNIQUES AND FAC METHOD IN NUMERICAL RESERVOIR SIMULATION \\Lambda
the reservoir. Combining the need for an accurate approximation of these moving features with the needOBJECT ORIENTED PROGRAMMING TECHNIQUES AND FAC METHOD IN NUMERICAL RESERVOIR SIMULATION \\Lambda in numerical simulation of flow through hydrocarbon reservoirs within limitations in computing time and memory
NUMERICAL NULL CONTROLLABILITY OF THE HEAT EQUATION THROUGH A LEAST SQUARES AND VARIATIONAL APPROACH
NUMERICAL NULL CONTROLLABILITY OF THE HEAT EQUATION THROUGH A LEAST SQUARES AND VARIATIONAL of null controls for the heat equation. The goal is to compute an approximation of controls that drives-space dimensional case. Key Words. Heat equation, Null controllability, Numerical approximation, Variational ap
Numerical null controllability of the 1D heat equation: dual methods
Paris-Sud XI, Université de
Numerical null controllability of the 1D heat equation: dual methods Enrique Fern for the 1D heat equation. The goal is to compute a control that drives (a numerical approximation of and efficiency. Keywords: one-dimensional heat equation, null controllability, finite element methods, dual meth
Numerical exact controllability of the 1D heat equation: dual algorithms
Sart, Remi
Numerical exact controllability of the 1D heat equation: dual algorithms Enrique Fern, which is closely related to the regularization effect of the heat equation. There, the controllability´andez-Cara and Arnaud M¨unch Abstract This work addresses the numerical approximation of distributed null controls (with
Numerical null controllability of semi-linear 1-D heat equations: fixed point, least squares and
Sart, Remi
Numerical null controllability of semi-linear 1-D heat equations: fixed point, least squares with the numerical computation of distributed null controls for semi- linear 1D heat equations, in the sublinear in [Fernandez-Cara & Zuazua, Null and approximate controllability for weakly blowing up semi-linear heat
NUMERICAL NULL CONTROLLABILITY OF THE HEAT EQUATION THROUGH A VARIATIONAL APPROACH
Sart, Remi
NUMERICAL NULL CONTROLLABILITY OF THE HEAT EQUATION THROUGH A VARIATIONAL APPROACH ARNAUD M for the heat equation. The goal is to compute a control that drives (an approximation of) the solution from¨UNCH AND PABLO PEDREGAL Abstract. This work is concerned with the numerical computation of null controls
NUMERICAL SIMULATIONS OF LONG TERM UNSATURATED FLOW AND ACID MINE DRAINAGE AT WASTE ROCK PILES
Aubertin, Michel
NUMERICAL SIMULATIONS OF LONG TERM UNSATURATED FLOW AND ACID MINE DRAINAGE AT WASTE ROCK PILES Omar present a numerical modeling study of unsaturated water flow and acid mine drainage in idealized (but of oxygen diffusion and acid mine drainage through the waste rock piles showed that oxygen is generally
Numerical evaluation of one-loop diagrams near exceptional momentum configurations
Walter T Giele; Giulia Zanderighi; E.W.N. Glover
2004-07-06
One problem which plagues the numerical evaluation of one-loop Feynman diagrams using recursive integration by part relations is a numerical instability near exceptional momentum configurations. In this contribution we will discuss a generic solution to this problem. As an example we consider the case of forward light-by-light scattering.
Experimental and numerical investigation of an optimized airfoil for vertical axis wind turbines
Experimental and numerical investigation of an optimized airfoil for vertical axis wind turbines and numerical verification of the per- formance of a new airfoil design for lift driven vertical-axis wind-turbines-driven vertical-axis wind-turbines VAWTs, with particular attention to floating installations (see Akimoto et al
Impact of Turbulence Closures and Numerical Errors for the Optimization of Flow Control Devices
Paris-Sud XI, Université de
Impact of Turbulence Closures and Numerical Errors for the Optimization of Flow Control Devices J the use of a Kriging-based global optimization method to determine optimal control parameters conduct an optimization process and measure the impact of numerical and modeling errors on the optimal
Fusion Engineering and Design 81 (2006) 14511458 Numerical study of MHD effect on liquid metal free
Abdou, Mohamed
2006-01-01
Fusion Engineering and Design 81 (2006) 14511458 Numerical study of MHD effect on liquid metal In this paper, we present numerical studies of liquid metal free surface jet characteristic behavior under with experimental observations. A 3D liquid metal MHD code based on an induced magnetic field formulation
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
Attribute Transformations on Numerical Applications to Stock Market and Economic Data
Lin, Tsau Young
Attribute Transformations on Numerical Databases Applications to Stock Market and Economic Data T of attribu- tion transformations, linear transformations, is applied to stock market and economic data. Some@cs.berkely.edu Abstract. The e#11;ects of attribute transformations on numerical data mining are investigated. Theoretical
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
Numerical Zoom for Multiscale Problems with an Application to Nuclear Waste Disposal
Numerical Zoom for Multiscale Problems with an Application to Nuclear Waste Disposal Jean of a nuclear waste repository site. Key words: Multiscale, Finite Element, Domain Decomposition, Chimera, Numerical Zoom, Nuclear Waste. PACS: 02.30.Jr, 47.11.Fg, 28.41.Kw, 47.55.P- 1 Introduction The present paper
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