Time-Dependent Seismic Tomography of the Coso Geothermal Area...
to: navigation, search OpenEI Reference LibraryAdd to library Conference Proceedings: Time-Dependent Seismic Tomography of the Coso Geothermal Area, 1996-2004 Abstract...
Time-dependent seismic tomography of the Coso geothermal area...
to: navigation, search OpenEI Reference LibraryAdd to library Conference Proceedings: Time-dependent seismic tomography of the Coso geothermal area, 1996-2004 Abstract...
Time-dependent seismic tomography and its application to the...
changes in Earth structure are commonly determined using local earthquake tomography computer programs that invert multiple seismic-wave arrival time data sets separately and...
Evaluation of Cross-Hole Seismic Tomography for Imaging Low Resistance...
Office of Scientific and Technical Information (OSTI)
Evaluation of Cross-Hole Seismic Tomography for Imaging Low Resistance Intervals and ... Title: Evaluation of Cross-Hole Seismic Tomography for Imaging Low Resistance Intervals ...
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Wang, Hui; Liu, Li; Lu, Yang; Pan, Pan; Hooker, Jacob M.; Fowler, Joanna S.; Tonge, Peter J.
2015-07-14
PT70 is a diaryl ether inhibitor of InhA, the enoyl-ACP reductase in the Mycobacterium tuberculosis fatty acid biosynthesis pathway. It has a residence time of 24 min on the target, and also shows antibacterial activity in a mouse model of tuberculosis infection. Due to the interest in studying target tissue pharmacokinetics of PT70, we developed a method to radiolabel PT70 with carbon-11 and have studied its pharmacokinetics in mice and baboons using positron emission tomography.
Nurhandoko, Bagus Endar B.; Wely, Woen; Setiadi, Herlan; Riyanto, Erwin
2015-04-16
It is already known that tomography has a great impact for analyzing and mapping unknown objects based on inversion, travel time as well as waveform inversion. Therefore, tomography has used in wide area, not only in medical but also in petroleum as well as mining. Recently, tomography method is being applied in several mining industries. A case study of tomography imaging has been carried out in DOZ ( Deep Ore Zone ) block caving mine, Tembagapura, Papua. Many researchers are undergoing to investigate the properties of DOZ cave not only outside but also inside which is unknown. Tomography takes a part for determining this objective.The sources are natural from the seismic events that caused by mining induced seismicity and rocks deformation activity, therefore it is called as passive seismic. These microseismic travel time data are processed by Simultaneous Iterative Reconstruction Technique (SIRT). The result of the inversion can be used for DOZ cave monitoring. These information must be used for identifying weak zone inside the cave. In addition, these results of tomography can be used to determine DOZ and cave information to support mine activity in PT. Freeport Indonesia.
Dynamic ray tracing and traveltime corrections for global seismic tomography
Tian Yue Hung, S.-H.; Nolet, Guust; Montelli, Raffaella; Dahlen, F.A.
2007-09-10
We present a dynamic ray tracing program for a spherically symmetric Earth that may be used to compute Frechet kernels for traveltime and amplitude anomalies at finite frequency. The program works for arbitrarily defined phases and background models. The numerical precisions of kinematic and dynamic ray tracing are optimized to produce traveltime errors under 0.1 s, which is well below the data uncertainty in global seismology. This tolerance level is obtained for an integration step size of about 20 km for the most common seismic phases. We also give software to compute ellipticity, crustal and topographic corrections and attenuation.
Seismic Surface-Wave Tomography of Waste Sites - Final Report
Long, Timothy L.
2000-09-14
The objective of this study was to develop analysis programs for surface-wave group-velocity tomography, and apply these to three test areas. We succeeded by obtaining data covering two square areas that were 30 meters on a side, and a third area that was 16 meters on a side, in addition to a collaborative effort wherein we processed data from the Oak Ridge National Laboratory site. At all sites, usable group velocities were obtained for frequencies from 16 to 50 Hz using a sledgehammer source. The resulting tomographic images and velocity anomalies were sufficient to delineate suspected burial trenches (one 4-meters deep) and anomalous velocity structure related to rocks and disturbed soil. The success was not uniform because in portions of one area the inversion for shear-wave structure became unstable. More research is needed to establish a more robust inversion technique.
Seismic Surface-Wave Tomography of Waste Sites
Leland Timothy Long
2002-12-17
Surface-wave group-velocity tomography is an efficient way to obtain images of the group velocity over a test area. Because Rayleigh-wave group velocity depends on frequency, there are separate images for each frequency. Thus, at each point in these images the group velocities define a dispersion curve, a curve that relates group velocity to frequency. The objective of this study has been to find an accurate and efficient way to find the shear-wave structure from these dispersion curves. The conventional inversion techniques match theoretical and observed dispersion curves to determine the structure. These conventional methods do not always succeed in correctly differentiating the fundamental and higher modes, and for some velocity structures can become unstable. In this research a perturbation technique was developed. The perturbation method allows the pre-computation of a global inversion matrix which improves efficiency in obtaining solutions for the structure. Perturbation methods are stable and mimic the averaging process in wave propagation; hence. leading to more accurate solutions. Finite difference techniques and synthetic trace generation techniques were developed to define the perturbations. A new differential trace technique was developed for slight variations in dispersion. The improvements in analysis speed and the accuracy of the solution could lead to real-time field analysis systems, making it possible to obtain immediate results or to monitor temporal change in structure, such as might develop in using fluids for soil remediation.
Setiadi, Herlan; Nurhandoko, Bagus Endar B.; Wely, Woen; Riyanto, Erwin
2015-04-16
Fracture prediction in a block cave of underground mine is very important to monitor the structure of the fracture that can be harmful to the mining activities. Many methods can be used to obtain such information, such as TDR (Time Domain Relectometry) and open hole. Both of them have limitations in range measurement. Passive seismic tomography is one of the subsurface imaging method. It has advantage in terms of measurements, cost, and rich of rock physical information. This passive seismic tomography studies using Fresnel zone to model the wavepath by using frequency parameter. Fresnel zone was developed by Nurhandoko in 2000. The result of this study is tomography of P and S wave velocity which can predict position of fracture. The study also attempted to use sum of the wavefronts to obtain position and time of seismic event occurence. Fresnel zone tomography and the summation wavefront can predict location of geological structure of mine area as well.
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Fisk, Mark D.; Pasyanos, Michael E.
2016-05-03
Characterizing regional seismic signals continues to be a difficult problem due to their variability. Calibration of these signals is very important to many aspects of monitoring underground nuclear explosions, including detecting seismic signals, discriminating explosions from earthquakes, and reliably estimating magnitude and yield. Amplitude tomography, which simultaneously inverts for source, propagation, and site effects, is a leading method of calibrating these signals. A major issue in amplitude tomography is the data quality of the input amplitude measurements. Pre-event and prephase signal-to-noise ratio (SNR) tests are typically used but can frequently include bad signals and exclude good signals. The deficiencies ofmore » SNR criteria, which are demonstrated here, lead to large calibration errors. To ameliorate these issues, we introduce a semi-automated approach to assess the bandwidth of a spectrum where it behaves physically. We determine the maximum frequency (denoted as Fmax) where it deviates from this behavior due to inflections at which noise or spurious signals start to bias the spectra away from the expected decay. We compare two amplitude tomography runs using the SNR and new Fmax criteria and show significant improvements to the stability and accuracy of the tomography output for frequency bands higher than 2 Hz by using our assessments of valid S-wave bandwidth. We compare Q estimates, P/S residuals, and some detailed results to explain the improvements. Lastly, for frequency bands higher than 4 Hz, needed for effective P/S discrimination of explosions from earthquakes, the new bandwidth criteria sufficiently fix the instabilities and errors so that the residuals and calibration terms are useful for application.« less
Time-Dependent Reliability Analysis
Energy Science and Technology Software Center (OSTI)
1999-10-27
FRANTIC-3 was developed to evaluate system unreliability using time-dependent techniques. The code provides two major options: to evaluate standby system unavailability or, in addition to the unavailability to calculate the total system failure probability by including both the unavailability of the system on demand as well as the probability that it will operate for an arbitrary time period following the demand. The FRANTIC-3 time dependent reliability models provide a large selection of repair and testingmore » policies applicable to standby or continously operating systems consisting of periodically tested, monitored, and non-repairable (non-testable) components. Time-dependent and test frequency dependent failures, as well as demand stress related failure, test-caused degradation and wear-out, test associated human errors, test deficiencies, test override, unscheduled and scheduled maintenance, component renewal and replacement policies, and test strategies can be prescribed. The conditional system unavailabilities associated with the downtimes of the user specified failed component are also evaluated. Optionally, the code can perform a sensitivity study for system unavailability or total failure probability to the failure characteristics of the standby components.« less
Brauchler, R.; Doetsch, J.; Dietrich, P.; Sauter, M.
2012-01-10
In this study, hydraulic and seismic tomographic measurements were used to derive a site-specific relationship between the geophysical parameter p-wave velocity and the hydraulic parameters, diffusivity and specific storage. Our field study includes diffusivity tomograms derived from hydraulic travel time tomography, specific storage tomograms, derived from hydraulic attenuation tomography, and p-wave velocity tomograms, derived from seismic tomography. The tomographic inversion was performed in all three cases with the SIRT (Simultaneous Iterative Reconstruction Technique) algorithm, using a ray tracing technique with curved trajectories. The experimental set-up was designed such that the p-wave velocity tomogram overlaps the hydraulic tomograms by half. The experiments were performed at a wellcharacterized sand and gravel aquifer, located in the Leine River valley near Gttingen, Germany. Access to the shallow subsurface was provided by direct-push technology. The high spatial resolution of hydraulic and seismic tomography was exploited to derive representative site-specific relationships between the hydraulic and geophysical parameters, based on the area where geophysical and hydraulic tests were performed. The transformation of the p-wave velocities into hydraulic properties was undertaken using a k-means cluster analysis. Results demonstrate that the combination of hydraulic and geophysical tomographic data is a promising approach to improve hydrogeophysical site characterization.
Reducing Disk Storage of Full-3D Seismic Waveform Tomography (F3DT) Through Lossy Online Compression
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Lindstrom, Peter; Chen, Po; Lee, En-Jui
2016-05-05
Full-3D seismic waveform tomography (F3DT) is the latest seismic tomography technique that can assimilate broadband, multi-component seismic waveform observations into high-resolution 3D subsurface seismic structure models. The main drawback in the current F3DT implementation, in particular the scattering-integral implementation (F3DT-SI), is the high disk storage cost and the associated I/O overhead of archiving the 4D space-time wavefields of the receiver- or source-side strain tensors. The strain tensor fields are needed for computing the data sensitivity kernels, which are used for constructing the Jacobian matrix in the Gauss-Newton optimization algorithm. In this study, we have successfully integrated a lossy compression algorithmmore » into our F3DT SI workflow to significantly reduce the disk space for storing the strain tensor fields. The compressor supports a user-specified tolerance for bounding the error, and can be integrated into our finite-difference wave-propagation simulation code used for computing the strain fields. The decompressor can be integrated into the kernel calculation code that reads the strain fields from the disk and compute the data sensitivity kernels. During the wave-propagation simulations, we compress the strain fields before writing them to the disk. To compute the data sensitivity kernels, we read the compressed strain fields from the disk and decompress them before using them in kernel calculations. Experiments using a realistic dataset in our California statewide F3DT project have shown that we can reduce the strain-field disk storage by at least an order of magnitude with acceptable loss, and also improve the overall I/O performance of the entire F3DT-SI workflow significantly. The integration of the lossy online compressor may potentially open up the possibilities of the wide adoption of F3DT-SI in routine seismic tomography practices in the near future.« less
Topic 5: Time-Dependent Behavior
Pfeiffer, P.A. [Argonne National Lab., IL (United States); Tanabe, Tada-aki [Nagoya Univ. (Japan). Dept. of Civil Engineering
1991-12-31
This chapter is a report of the material presented at the International Workshop on Finite Element Analysis of Reinforced Concrete, Session 4 -- Time Dependent Behavior, held at Columbia University, New York on June 3--6, 1991. Dr. P.A. Pfeiffer presented recent developments in time-dependent behavior of concrete and Professor T. Tanabe presented a review of research in Japan on time-dependent behavior of concrete. The chapter discusses the recent research of time-dependent behavior of concrete in the past few years in both the USA-European and Japanese communities. The author appreciates the valuable information provided by Zdenek P. Bazant in preparing the USA-European Research section.
Investigations of Low Temperature Time Dependent Cracking
Van der Sluys, W A; Robitz, E S; Young, B A; Bloom, J
2002-09-30
The objective of this project was to investigate metallurgical and mechanical phenomena associated with time dependent cracking of cold bent carbon steel piping at temperatures between 327 C and 360 C. Boiler piping failures have demonstrated that understanding the fundamental metallurgical and mechanical parameters controlling these failures is insufficient to eliminate it from the field. The results of the project consisted of the development of a testing methodology to reproduce low temperature time dependent cracking in laboratory specimens. This methodology was used to evaluate the cracking resistance of candidate heats in order to identify the factors that enhance cracking sensitivity. The resultant data was integrated into current available life prediction tools.
Syracuse, E. M.; Maceira, M.; Zhang, H.; Thurber, C. H.
2015-02-18
Joint inversions of seismic data recover models that simultaneously fit multiple constraints while playing upon the strengths of each data type. Here, we jointly invert 14 years of local earthquake body wave arrival times from the Alaska Volcano Observatory catalog and Rayleigh wave dispersion curves based upon ambient noise measurements for local V_{p}, V_{s}, and hypocentral locations at Akutan and Makushin Volcanoes using a new joint inversion algorithm.The velocity structure and relocated seismicity of both volcanoes are significantly more complex than many other volcanoes studied using similar techniques. Seismicity is distributed among several areas beneath or beyond the flanks of both volcanoes, illuminating a variety of volcanic and tectonic features. The velocity structures of the two volcanoes are exemplified by the presence of narrow high-V_{p} features in the near surface, indicating likely current or remnant pathways of magma to the surface. A single broad low-V_{p} region beneath each volcano is slightly offset from each summit and centered at approximately 7 km depth, indicating a potential magma chamber, where magma is stored over longer time periods. Differing recovery capabilities of the Vp and Vs datasets indicate that the results of these types of joint inversions must be interpreted carefully.
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Syracuse, E. M.; Maceira, M.; Zhang, H.; Thurber, C. H.
2015-02-18
Joint inversions of seismic data recover models that simultaneously fit multiple constraints while playing upon the strengths of each data type. Here, we jointly invert 14 years of local earthquake body wave arrival times from the Alaska Volcano Observatory catalog and Rayleigh wave dispersion curves based upon ambient noise measurements for local Vp, Vs, and hypocentral locations at Akutan and Makushin Volcanoes using a new joint inversion algorithm.The velocity structure and relocated seismicity of both volcanoes are significantly more complex than many other volcanoes studied using similar techniques. Seismicity is distributed among several areas beneath or beyond the flanks ofmore » both volcanoes, illuminating a variety of volcanic and tectonic features. The velocity structures of the two volcanoes are exemplified by the presence of narrow high-Vp features in the near surface, indicating likely current or remnant pathways of magma to the surface. A single broad low-Vp region beneath each volcano is slightly offset from each summit and centered at approximately 7 km depth, indicating a potential magma chamber, where magma is stored over longer time periods. Differing recovery capabilities of the Vp and Vs datasets indicate that the results of these types of joint inversions must be interpreted carefully.« less
Trichandi, Rahmantara; Yudistira, Tedi; Nugraha, Andri Dian; Zulhan, Zulfakriza; Saygin, Erdinc
2015-04-24
Ambient noise tomography is relatively a new method for imaging the shallow structure of the Earth subsurface. We presents the application of this method to produce a Rayleigh wave group velocity maps around the Merapi Volcano, Central Java. Rayleigh waves group velocity maps were reconstructed from the cross-correlation of ambient noise recorded by the DOMERAPI array which consists 43 broadband seismometers. In the processing stage, we first filtered the observation data to separatethe noise from the signal that dominated by the strong volcanic activities. Next, we cross-correlate the filtered data and stack to obtain the Green’s function for all possible station pairs. Then we carefully picked the peak of each Green’s function to estimate the dispersion trend and appliedMultiple Filter Technique to obtain the dispersion curve. Inter-station group velocity curvesare inverted to produceRayleigh wave group velocity maps for periods 1 to 10 s. The resulted Rayleigh group velocity maps show the interesting features around the Merapi Volcano which generally agree with the previous studies. Merapi-Lawu Anomaly (MLA) is emerged as a relatively low anomaly in our group velocity maps.
Time dependent particle emission from fission products (Conference...
Office of Scientific and Technical Information (OSTI)
Time dependent particle emission from fission products Citation Details In-Document Search Title: Time dependent particle emission from fission products Decay heating following ...
Elbring, G.J.; Narbutovskih, S.M.
1994-02-01
This integrated test plan describes the demonstration of the crosswell acoustic tomography technique as part of the Volatile Organic Compounds-Arid Integrated Demonstration (VOC-Arid ID). The purpose of this demonstration is to image the subsurface seismic velocity structure and to relate the resulting velocity model to lithology and saturation. In fiscal year (FY) 1994 an initial fielding will test three different downhole sources at two different sites at the Hanford US Department of Energy facility to identify which sources will provide the energy required to propagate between existing steel-cased wells at these two sites. Once this has been established, a second fielding will perform a full compressional and shear wave tomographic survey at the most favorable site. Data reduction, analysis, and interpretation of this full data set will be completed by the end of this fiscal year. Data collection for a second survey will be completed by the end of the fiscal year, and data reduction for this data set will be completed in FY 1995. The specific need is detailed subsurface characterization with minimum intrusion. This technique also has applications for long term vadose zone monitoring for both Resource Conservation and Recovery Act (RCRA) waste storage facilities and for remediation monitoring. Images produced are continuous between boreholes. This is a significant improvement over the single point data derived solely from core information. Saturation changes, either naturally occurring (e.g., perched water tables) or remediation induced (e.g., water table mounding from injection wells or during inwell air sparging) could be imaged. These crosswell data allow optimal borehole placement for groundwater remediation, associated monitoring wells and possibly evaluation of the effective influence of a particular remediation technique.
Effect of Cusps in Time-Dependent Quantum Mechanics (Journal...
Office of Scientific and Technical Information (OSTI)
Effect of Cusps in Time-Dependent Quantum Mechanics Title: Effect of Cusps in Time-Dependent Quantum Mechanics Authors: Yang, Zeng-hui ; Maitra, Neepa T. ; Burke, Kieron ...
Time dependence of Bragg forward scattering and self-seeding...
Office of Scientific and Technical Information (OSTI)
Journal Article: Time dependence of Bragg forward scattering and self-seeding of hard x-ray free-electron lasers Citation Details In-Document Search Title: Time dependence of Bragg ...
Efficient Real-Time Time-Dependent Density Functional Theory...
Office of Scientific and Technical Information (OSTI)
Efficient Real-Time Time-Dependent Density Functional Theory Method and its Application to a Collision of an Ion with a 2D Material Title: Efficient Real-Time Time-Dependent ...
Time dependent deformation of metals at room temperature. (Conference...
Office of Scientific and Technical Information (OSTI)
Title: Time dependent deformation of metals at room temperature. Abstract not provided. Authors: Deibler, Lisa Anne ; Boyce, Brad Lee ; Puskar, Joseph D. Publication Date: ...
A time-dependent formulation of multi-reference perturbation...
Office of Scientific and Technical Information (OSTI)
Title: A time-dependent formulation of multi-reference perturbation theory Authors: Sokolov, Alexander Yu. 1 ; Chan, Garnet Kin-Lic 1 + Show Author Affiliations Department of ...
Finite Lifetime Effects on the Polarizability Within Time-dependent...
Office of Scientific and Technical Information (OSTI)
Density-functional Theory Citation Details In-Document Search Title: Finite Lifetime Effects on the Polarizability Within Time-dependent Density-functional Theory The research ...
Effect of noise on time-dependent quantum chaos
Ott, E.; Antonsen T.M. Jr.; Hanson, J.D.
1984-12-03
The dynamics of a time-dependent quantum system can be qualitatively different from that of its classical counterpart when the latter is chaotic. It is shown that small noise can strongly alter this situation.
Fragment Yields Calculated in a Time-Dependent Microscopic Theory...
Office of Scientific and Technical Information (OSTI)
Fragment Yields Calculated in a Time-Dependent Microscopic Theory of Fission Citation ... Word Cloud More Like This Full Text preview image File size NAView Full Text View Full ...
Time dependent particle emission from fission products (Conference) |
Office of Scientific and Technical Information (OSTI)
SciTech Connect Time dependent particle emission from fission products Citation Details In-Document Search Title: Time dependent particle emission from fission products Decay heating following nuclear fission is an important factor in the design of nuclear facilities; impacting a variety of aspects ranging from cooling requirements to shielding design. Calculations of decay heat, often assumed to be a simple product of activity and average decay product energy, are complicated by the so
Two-stream instability with time-dependent drift velocity
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Qin, Hong; Davidson, Ronald C.
2014-06-26
The classical two-stream instability driven by a constant relative drift velocity between two plasma components is extended to the case with time-dependent drift velocity. A solution method is developed to rigorously define and calculate the instability growth rate for linear perturbations relative to the time-dependent unperturbed two-stream motions. The stability diagrams for the oscillating two-stream instability are presented over a large region of parameter space. It is shown that the growth rate for the classical two-stream instability can be significantly reduced by adding an oscillatory component to the relative drift velocity.
Bellomo, Bruno; De Pasquale, Antonella; Gualdi, Giulia; Marzolino, Ugo
2010-12-15
We propose a procedure to fully reconstruct the time-dependent coefficients of convolutionless non-Markovian dissipative generators via a finite number of experimental measurements. By combining a tomography-based approach with a proper data sampling, our proposal allows to relate the time-dependent coefficients governing the dissipative evolution of a quantum system to experimentally accessible quantities. The proposed scheme not only provides a way to retrieve the full information about potentially unknown dissipative coefficients, but also, most valuably, can be employed as a reliable consistency test for the approximations involved in the theoretical derivation of a given non-Markovian convolutionless master equation.
Full Reviews: Seismicity and Seismic
Broader source: Energy.gov [DOE]
Below are the project presentations and respective peer reviewer comments for Seismicity and Seismic.
Sublinear scaling for time-dependent stochastic density functional theory
Gao, Yi; Neuhauser, Daniel; Baer, Roi; Rabani, Eran
2015-01-21
A stochastic approach to time-dependent density functional theory is developed for computing the absorption cross section and the random phase approximation (RPA) correlation energy. The core idea of the approach involves time-propagation of a small set of stochastic orbitals which are first projected on the occupied space and then propagated in time according to the time-dependent Kohn-Sham equations. The evolving electron density is exactly represented when the number of random orbitals is infinite, but even a small number (≈16) of such orbitals is enough to obtain meaningful results for absorption spectrum and the RPA correlation energy per electron. We implement the approach for silicon nanocrystals using real-space grids and find that the overall scaling of the algorithm is sublinear with computational time and memory.
Time-dependent first-principles approaches to PV materials
Miyamoto, Yoshiyuki
2013-12-10
Computational scheme for designing photovoltaic (PV) materials is presented. First-principles electron dynamics of photo-excitation and subsequent electron-hole splitting is performed based on the time-dependent density functional theory. Photo-induced enhancement of dipole moment was observed in a polar crystal and a donor-acceptor molecular pair. These experiences will pave a way to design PV material from first-principles simulations.
Neutrino flavor instabilities in a time-dependent supernova model
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Abbar, Sajad; Duan, Huaiyu
2015-10-19
In this study, a dense neutrino medium such as that inside a core-collapse supernova can experience collective flavor conversion or oscillations because of the neutral-current weak interaction among the neutrinos. This phenomenon has been studied in a restricted, stationary supernova model which possesses the (spatial) spherical symmetry about the center of the supernova and the (directional) axial symmetry around the radial direction. Recently it has been shown that these spatial and directional symmetries can be broken spontaneously by collective neutrino oscillations. In this letter we analyze the neutrino flavor instabilities in a time-dependent supernova model. Our results show that collectivemore »neutrino oscillations start at approximately the same radius in both the stationary and time-dependent supernova models unless there exist very rapid variations in local physical conditions on timescales of a few microseconds or shorter. Our results also suggest that collective neutrino oscillations can vary rapidly with time in the regimes where they do occur which need to be studied in time-dependent supernova models.« less
Progress Report on Alloy 617 Time Dependent Allowables
Wright, Julie Knibloe
2015-06-01
Time dependent allowable stresses are required in the ASME Boiler and Pressure Vessel Code for design of components in the temperature range where time dependent deformation (i.e., creep) is expected to become significant. There are time dependent allowable stresses in Section IID of the Code for use in the non-nuclear construction codes, however, there are additional criteria that must be considered in developing time dependent allowables for nuclear components. These criteria are specified in Section III NH. St is defined as the lesser of three quantities: 100% of the average stress required to obtain a total (elastic, plastic, primary and secondary creep) strain of 1%; 67% of the minimum stress to cause rupture; and 80% of the minimum stress to cause the initiation of tertiary creep. The values are reported for a range of temperatures and for time increments up to 100,000 hours. These values are determined from uniaxial creep tests, which involve the elevated temperature application of a constant load which is relatively small, resulting in deformation over a long time period prior to rupture. The stress which is the minimum resulting from these criteria is the time dependent allowable stress St. In this report data from a large number of creep and creep-rupture tests on Alloy 617 are analyzed using the ASME Section III NH criteria. Data which are used in the analysis are from the ongoing DOE sponsored high temperature materials program, form Korea Atomic Energy Institute through the Generation IV VHTR Materials Program and historical data from previous HTR research and vendor data generated in developing the alloy. It is found that the tertiary creep criterion determines St at highest temperatures, while the stress to cause 1% total strain controls at low temperatures. The ASME Section III Working Group on Allowable Stress Criteria has recommended that the uncertainties associated with determining the onset of tertiary creep and the lack of significant
Translation invariant time-dependent solutions to massive gravity
Mourad, J.; Steer, D.A. E-mail: steer@apc.univ-paris7.fr
2013-12-01
Homogeneous time-dependent solutions of massive gravity generalise the plane wave solutions of the linearised Fierz-Pauli equations for a massive spin-two particle, as well as the Kasner solutions of General Relativity. We show that they also allow a clear counting of the degrees of freedom and represent a simplified framework to work out the constraints, the equations of motion and the initial value formulation. We work in the vielbein formulation of massive gravity, find the phase space resulting from the constraints and show that several disconnected sectors of solutions exist some of which are unstable. The initial values determine the sector to which a solution belongs. Classically, the theory is not pathological but quantum mechanically the theory may suffer from instabilities. The latter are not due to an extra ghost-like degree of freedom.
Translation invariant time-dependent solutions to massive gravity II
Mourad, J.; Steer, D.A. E-mail: steer@apc.univ-paris7.fr
2014-06-01
This paper is a sequel to JCAP 12 (2013) 004 and is also devoted to translation-invariant solutions of ghost-free massive gravity in its moving frame formulation. Here we consider a mass term which is linear in the vielbein (corresponding to a ?{sub 3} term in the 4D metric formulation) in addition to the cosmological constant. We determine explicitly the constraints, and from the initial value formulation show that the time-dependent solutions can have singularities at a finite time. Although the constraints give, as in the ?{sub 1} case, the correct number of degrees of freedom for a massive spin two field, we show that the lapse function can change sign at a finite time causing a singular time evolution. This is very different to the ?{sub 1} case where time evolution is always well defined. We conclude that the ?{sub 3} mass term can be pathological and should be treated with care.
Spin-orbit torque induced spike-timing dependent plasticity
Sengupta, Abhronil Al Azim, Zubair; Fong, Xuanyao; Roy, Kaushik
2015-03-02
Nanoelectronic devices that mimic the functionality of synapses are a crucial requirement for performing cortical simulations of the brain. In this work, we propose a ferromagnet-heavy metal heterostructure that employs spin-orbit torque to implement spike-timing dependent plasticity. The proposed device offers the advantage of decoupled spike transmission and programming current paths, thereby leading to reliable operation during online learning. Possible arrangement of such devices in a crosspoint architecture can pave the way for ultra-dense neural networks. Simulation studies indicate that the device has the potential of achieving pico-Joule level energy consumption (maximum 2 pJ per synaptic event) which is comparable to the energy consumption for synaptic events in biological synapses.
Time-Dependent, Parallel Neutral Particle Transport Code System.
Energy Science and Technology Software Center (OSTI)
2009-09-10
Version 00 PARTISN (PARallel, TIme-Dependent SN) is the evolutionary successor to CCC-547/DANTSYS. The PARTISN code package is a modular computer program package designed to solve the time-independent or dependent multigroup discrete ordinates form of the Boltzmann transport equation in several different geometries. The modular construction of the package separates the input processing, the transport equation solving, and the post processing (or edit) functions into distinct code modules: the Input Module, the Solver Module, and themore » Edit Module, respectively. PARTISN is the evolutionary successor to the DANTSYSTM code system package. The Input and Edit Modules in PARTISN are very similar to those in DANTSYS. However, unlike DANTSYS, the Solver Module in PARTISN contains one, two, and three-dimensional solvers in a single module. In addition to the diamond-differencing method, the Solver Module also has Adaptive Weighted Diamond-Differencing (AWDD), Linear Discontinuous (LD), and Exponential Discontinuous (ED) spatial differencing methods. The spatial mesh may consist of either a standard orthogonal mesh or a block adaptive orthogonal mesh. The Solver Module may be run in parallel for two and three dimensional problems. One can now run 1-D problems in parallel using Energy Domain Decomposition (triggered by Block 5 input keyword npeg>0). EDD can also be used in 2-D/3-D with or without our standard Spatial Domain Decomposition. Both the static (fixed source or eigenvalue) and time-dependent forms of the transport equation are solved in forward or adjoint mode. In addition, PARTISN now has a probabilistic mode for Probability of Initiation (static) and Probability of Survival (dynamic) calculations. Vacuum, reflective, periodic, white, or inhomogeneous boundary conditions are solved. General anisotropic scattering and inhomogeneous sources are permitted. PARTISN solves the transport equation on orthogonal (single level or block-structured AMR) grids in 1-D
Seismic Imaging of the Earth's Interior (LBNL Summer Lecture Series)
Romanowicz, Barbara
2011-04-28
Summer Lecture Series 2006: Earth scientist Barbara Romanowicz discusses how she explores the deep structure and dynamics of the Earth using seismic tomography.
Seismic waveform viewer, processor and calculator
Energy Science and Technology Software Center (OSTI)
2015-02-15
SWIFT is a computer code that is designed to do research level signal analysis on seismic waveforms, including visualization, filtering and measurement. LLNL is using this code, amplitude and global tomography efforts.
MaRIE: A facility for time-dependent materials science at the...
Office of Scientific and Technical Information (OSTI)
MaRIE: A facility for time-dependent materials science at the mesoscale Citation Details In-Document Search Title: MaRIE: A facility for time-dependent materials science at the ...
Problem-free time-dependent variational principle for open quantum...
Office of Scientific and Technical Information (OSTI)
formulation of TDVP applied to a non-stochastic open system Schrdinger equation L. ... QUANTUM SYSTEMS; SCHROEDINGER EQUATION; STOCHASTIC PROCESSES; TIME DEPENDENCE; ...
Boyarinov, V. F.; Kondrushin, A. E.; Fomichenko, P. A.
2012-07-01
Finite-difference time-dependent equations of Surface Harmonics method have been obtained for plane geometry. Verification of these equations has been carried out by calculations of tasks from 'Benchmark Problem Book ANL-7416'. The capacity and efficiency of the Surface Harmonics method have been demonstrated by solution of the time-dependent neutron transport equation in diffusion approximation. The results of studies showed that implementation of Surface Harmonics method for full-scale calculations will lead to a significant progress in the efficient solution of the time-dependent neutron transport problems in nuclear reactors. (authors)
Measurement of Time-Dependent CP Asymmetry in the Decay of a...
Office of Scientific and Technical Information (OSTI)
a Long-Lived Neutral Kaon at BaBar Citation Details In-Document Search Title: Measurement of Time-Dependent CP Asymmetry in the Decay of a Neutral B Meson to a JPsi and a Long-Li...
Fragment Yields Calculated in a Time-Dependent Microscopic Theory of
Office of Scientific and Technical Information (OSTI)
Fission (Technical Report) | SciTech Connect Fragment Yields Calculated in a Time-Dependent Microscopic Theory of Fission Citation Details In-Document Search Title: Fragment Yields Calculated in a Time-Dependent Microscopic Theory of Fission Authors: Younes, W ; Gogny, D Publication Date: 2012-09-28 OSTI Identifier: 1053671 Report Number(s): LLNL-TR-586678 DOE Contract Number: W-7405-ENG-48 Resource Type: Technical Report Research Org: Lawrence Livermore National Laboratory (LLNL),
The precise time-dependent solution of the Fokker–Planck equation with anomalous diffusion
Guo, Ran; Du, Jiulin
2015-08-15
We study the time behavior of the Fokker–Planck equation in Zwanzig’s rule (the backward-Ito’s rule) based on the Langevin equation of Brownian motion with an anomalous diffusion in a complex medium. The diffusion coefficient is a function in momentum space and follows a generalized fluctuation–dissipation relation. We obtain the precise time-dependent analytical solution of the Fokker–Planck equation and at long time the solution approaches to a stationary power-law distribution in nonextensive statistics. As a test, numerically we have demonstrated the accuracy and validity of the time-dependent solution. - Highlights: • The precise time-dependent solution of the Fokker–Planck equation with anomalous diffusion is found. • The anomalous diffusion satisfies a generalized fluctuation–dissipation relation. • At long time the time-dependent solution approaches to a power-law distribution in nonextensive statistics. • Numerically we have demonstrated the accuracy and validity of the time-dependent solution.
Office of Scientific and Technical Information (OSTI)
... event tree SMA seismic margins analysis SNF spent nuclear fuel SRTC site rail transfer cart SSC ... Safeguards and Security System Safeguards and Security Entire NA NA NA NA System ...
MaRIE: A facility for time-dependent materials science at the mesoscale
Office of Scientific and Technical Information (OSTI)
(Technical Report) | SciTech Connect Technical Report: MaRIE: A facility for time-dependent materials science at the mesoscale Citation Details In-Document Search Title: MaRIE: A facility for time-dependent materials science at the mesoscale Authors: Barnes, Cris William [1] ; Kippen, Karen Elizabeth [1] + Show Author Affiliations Los Alamos National Laboratory Publication Date: 2015-02-11 OSTI Identifier: 1170260 Report Number(s): LA-UR-15-20995 DOE Contract Number: AC52-06NA25396 Resource
Charged Particle in a Time-dependent Electric Field: A White Noise Functional Approach
Gravador, E. B.; Bornales, J. B.; Liwanag, M. J.
2008-06-18
The propagator for a charged particle in a time-dependent electric field is calculated following Hida and Streit's framework where the propagator is the T-transform of Feynman functional. However, we have to regard the potential V = -qE({tau})x{identical_to}{radical}((m/({Dirac_h}/2{pi}))){xi}x following C. C. Bernido and M. V. Carpio-Bernido's prescription of time-dependent potentials. The result agrees with the limiting form of Eqn. (16) of N. Morgenstern Horing, H. L. Cui, and G. Fiorenza, when the magnetic field is switched off, and Eqn. (17) of [3] when the electric field is constant in time.
R. Quittmeyer
2006-09-25
This technical work plan (TWP) describes the efforts to develop and confirm seismic ground motion inputs used for preclosure design and probabilistic safety 'analyses and to assess the postclosure performance of a repository at Yucca Mountain, Nevada. As part of the effort to develop seismic inputs, the TWP covers testing and analyses that provide the technical basis for inputs to the seismic ground-motion site-response model. The TWP also addresses preparation of a seismic methodology report for submission to the U.S. Nuclear Regulatory Commission (NRC). The activities discussed in this TWP are planned for fiscal years (FY) 2006 through 2008. Some of the work enhances the technical basis for previously developed seismic inputs and reduces uncertainties and conservatism used in previous analyses and modeling. These activities support the defense of a license application. Other activities provide new results that will support development of the preclosure, safety case; these results directly support and will be included in the license application. Table 1 indicates which activities support the license application and which support licensing defense. The activities are listed in Section 1.2; the methods and approaches used to implement them are discussed in more detail in Section 2.2. Technical and performance objectives of this work scope are: (1) For annual ground motion exceedance probabilities appropriate for preclosure design analyses, provide site-specific seismic design acceleration response spectra for a range of damping values; strain-compatible soil properties; peak motions, strains, and curvatures as a function of depth; and time histories (acceleration, velocity, and displacement). Provide seismic design inputs for the waste emplacement level and for surface sites. Results should be consistent with the probabilistic seismic hazard analysis (PSHA) for Yucca Mountain and reflect, as appropriate, available knowledge on the limits to extreme ground motion at
2-D Time-Dependent Fuel Element, Thermal Analysis Code System.
Energy Science and Technology Software Center (OSTI)
2001-09-24
Version 00 WREM-TOODEE2 is a two dimensional, time-dependent, fuel-element thermal analysis program. Its primary purpose is to evaluate fuel-element thermal response during post-LOCA refill and reflood in a pressurized water reactor (PWR). TOODEE2 calculations are carried out in a two-dimensional mesh region defined in slab or cylindrical geometry by orthogonal grid lines. Coordinates which form order pairs are labeled x-y in slab geometry, and those in cylindrical geometry are labeled r-z for the axisymmetric casemore » and r-theta for the polar case. Conduction and radiation are the only heat transfer mechanisms assumed within the boundaries of the mesh region. Convective and boiling heat transfer mechanisms are assumed at the boundaries. The program numerically solves the two-dimensional, time-dependent, heat conduction equation within the mesh region. KEYWORDS: FUEL MANAGEMENT; HEAT TRANSFER; LOCA; PWR« less
Simulation of biochemical reactions with time-dependent rates by the rejection-based algorithm
Thanh, Vo Hong; Priami, Corrado
2015-08-07
We address the problem of simulating biochemical reaction networks with time-dependent rates and propose a new algorithm based on our rejection-based stochastic simulation algorithm (RSSA) [Thanh et al., J. Chem. Phys. 141(13), 134116 (2014)]. The computation for selecting next reaction firings by our time-dependent RSSA (tRSSA) is computationally efficient. Furthermore, the generated trajectory is exact by exploiting the rejection-based mechanism. We benchmark tRSSA on different biological systems with varying forms of reaction rates to demonstrate its applicability and efficiency. We reveal that for nontrivial cases, the selection of reaction firings in existing algorithms introduces approximations because the integration of reaction rates is very computationally demanding and simplifying assumptions are introduced. The selection of the next reaction firing by our approach is easier while preserving the exactness.
Nascimento, Daniel R.; DePrince, A. Eugene
2015-12-07
We present a combined cavity quantum electrodynamics/ab initio electronic structure approach for simulating plasmon-molecule interactions in the time domain. The simple Jaynes-Cummings-type model Hamiltonian typically utilized in such simulations is replaced with one in which the molecular component of the coupled system is treated in a fully ab initio way, resulting in a computationally efficient description of general plasmon-molecule interactions. Mutual polarization effects are easily incorporated within a standard ground-state Hartree-Fock computation, and time-dependent simulations carry the same formal computational scaling as real-time time-dependent Hartree-Fock theory. As a proof of principle, we apply this generalized method to the emergence of a Fano-like resonance in coupled molecule-plasmon systems; this feature is quite sensitive to the nanoparticle-molecule separation and the orientation of the molecule relative to the polarization of the external electric field.
Time-dependent density functional theory quantum transport simulation in non-orthogonal basis
Kwok, Yan Ho; Xie, Hang; Yam, Chi Yung; Chen, Guan Hua; Zheng, Xiao
2013-12-14
Basing on the earlier works on the hierarchical equations of motion for quantum transport, we present in this paper a first principles scheme for time-dependent quantum transport by combining time-dependent density functional theory (TDDFT) and Keldysh's non-equilibrium Green's function formalism. This scheme is beyond the wide band limit approximation and is directly applicable to the case of non-orthogonal basis without the need of basis transformation. The overlap between the basis in the lead and the device region is treated properly by including it in the self-energy and it can be shown that this approach is equivalent to a lead-device orthogonalization. This scheme has been implemented at both TDDFT and density functional tight-binding level. Simulation results are presented to demonstrate our method and comparison with wide band limit approximation is made. Finally, the sparsity of the matrices and computational complexity of this method are analyzed.
An analysis of time-dependence for Chernobyl fallout in Italy
Salvadori, G.; Ratti, S.P.; Belli, G.
1997-01-01
We present here an analysis of the airborne radioactivity measured in Italy after the Chernobyl accident. We provide some quality assurance, isolate suspicious data, and devise a mathematical model to aid in interpreting time-dependent fallout data. The model consists of an interpolating function whose parameters can be related to (1) the arrival time of the radioactive cloud; (2) the time of the maximum radioactive concentration; and (3) the decay-rate of airborne radioactivity as the pollutant cloud passes. Multiple arrivals of the radioactive cloud in a given site can also be considered. The parametrization can be used to estimate concentrations of {sup 137}Cs using measurements of {sup 131}I, {sup 103}Ru, or {sup 132}Te. The interpolating function is fitted to the data collected in several Italian Provinces. We feel this model is an useful tool for interpreting time-dependent fallout data. 13 refs., 7 figs., 13 tabs.
Jing, Longfei; Yang, Dong; Li, Hang; Zhang, Lu; Lin, Zhiwei; Li, Liling; Kuang, Longyu; Jiang, Shaoen Ding, Yongkun; Huang, Yunbao
2015-02-15
The x-ray drive on a capsule in an inertial confinement fusion setup is crucial for ignition. Unfortunately, a direct measurement has not been possible so far. We propose an angular radiation temperature simulation to predict the time-dependent drive on the capsule. A simple model, based on the view-factor method for the simulation of the radiation temperature, is presented and compared with the experimental data obtained using the OMEGA laser facility and the simulation results acquired with VISRAD code. We found a good agreement between the time-dependent measurements and the simulation results obtained using this model. The validated model was then used to analyze the experimental results from the Shenguang-III prototype laser facility. More specifically, the variations of the peak radiation temperatures at different view angles with the albedo of the hohlraum, the motion of the laser spots, the closure of the laser entrance holes, and the deviation of the laser power were investigated. Furthermore, the time-dependent radiation temperature at different orientations and the drive history on the capsule were calculated. The results indicate that the radiation temperature from “U20W112” (named according to the diagnostic hole ID on the target chamber) can be used to approximately predict the drive temperature on the capsule. In addition, the influence of the capsule on the peak radiation temperature is also presented.
The time-dependent chemistry of cometary debris in the solar corona
Pesnell, W. D.; Bryans, P.
2014-04-10
Recent improvements in solar observations have greatly progressed the study of sungrazing comets. They can now be imaged along the entirety of their perihelion passage through the solar atmosphere, revealing details of their composition and structure not measurable through previous observations in the less volatile region of the orbit further from the solar surface. Such comets are also unique probes of the solar atmosphere. The debris deposited by sungrazers is rapidly ionized and subsequently influenced by the ambient magnetic field. Measuring the spectral signature of the deposited material highlights the topology of the magnetic field and can reveal plasma parameters such as the electron temperature and density. Recovering these variables from the observable data requires a model of the interaction of the cometary species with the atmosphere through which they pass. The present paper offers such a model by considering the time-dependent chemistry of sublimated cometary species as they interact with the solar radiation field and coronal plasma. We expand on a previous simplified model by considering the fully time-dependent solutions of the emitting species' densities. To compare with observations, we consider a spherically symmetric expansion of the sublimated material into the corona and convert the time-dependent ion densities to radial profiles. Using emissivities from the CHIANTI database and plasma parameters derived from a magnetohydrodynamic simulation leads to a spatially dependent emission spectrum that can be directly compared with observations. We find our simulated spectra to be consistent with observation.
Towards time-dependent current-density-functional theory in the non-linear regime
Escartn, J. M.; Vincendon, M.; Dinh, P. M.; Suraud, E.; Romaniello, P.; Reinhard, P.-G.
2015-02-28
Time-Dependent Density-Functional Theory (TDDFT) is a well-established theoretical approach to describe and understand irradiation processes in clusters and molecules. However, within the so-called adiabatic local density approximation (ALDA) to the exchange-correlation (xc) potential, TDDFT can show insufficiencies, particularly in violently dynamical processes. This is because within ALDA the xc potential is instantaneous and is a local functional of the density, which means that this approximation neglects memory effects and long-range effects. A way to go beyond ALDA is to use Time-Dependent Current-Density-Functional Theory (TDCDFT), in which the basic quantity is the current density rather than the density as in TDDFT. This has been shown to offer an adequate account of dissipation in the linear domain when the Vignale-Kohn (VK) functional is used. Here, we go beyond the linear regime and we explore this formulation in the time domain. In this case, the equations become very involved making the computation out of reach; we hence propose an approximation to the VK functional which allows us to calculate the dynamics in real time and at the same time to keep most of the physics described by the VK functional. We apply this formulation to the calculation of the time-dependent dipole moment of Ca, Mg and Na{sub 2}. Our results show trends similar to what was previously observed in model systems or within linear response. In the non-linear domain, our results show that relaxation times do not decrease with increasing deposited excitation energy, which sets some limitations to the practical use of TDCDFT in such a domain of excitations.
Seismic Monitoring - Hanford Site
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Public Safety and Resource Protection (PSRP) Seismic Monitoring Public Safety and Resource Protection (PSRP) Public Safety and Resource Protection Home Cultural Resource Program and Curation Services Ecological Monitoring Environmental Surveillance Meteorology and Climatology Services Seismic Monitoring Email Email Page | Print Print Page | Text Increase Font Size Decrease Font Size Seismic Monitoring Seismic Monitoring Hanford Site Seismic Monitoring provides an uninterrupted collection of
Chakrabarti, Nikhil; Maity, Chandan; Schamel, Hans
2011-04-08
Compressional waves in a magnetized plasma of arbitrary resistivity are treated with the Lagrangian fluid approach. An exact nonlinear solution with a nontrivial space and time dependence is obtained with boundary conditions as in Harris' current sheet. The solution shows competition among hydrodynamic convection, magnetic field diffusion, and dispersion. This results in a collapse of density and the magnetic field in the absence of dispersion. The dispersion effects arrest the collapse of density but not of the magnetic field. A possible application is in the early stage of magnetic star formation.
Optical absorption in B{sub 13} cluster: A time-dependent density functional approach
Shinde, Ravindra; Tayade, Meenakshi
2013-02-05
The linear optical absorption spectra of three isomers of planar boron cluster B{sub 13} are calculated using time-dependent spin-polarized density functional approach. The geometries of these cluster are optimized at the B3LYP/6-311+G* level of theory. Even though the isomers are almost degenerate, the calculated spectra are quite different, indicating a strong structure-property relationship. Therefore, these computed spectra can be used in the photo-absorption experiments to distinguish between different isomers of a cluster.
Solvent effects in time-dependent self-consistent field methods. I. Optical response calculations
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Bjorgaard, J. A.; Kuzmenko, V.; Velizhanin, K. A.; Tretiak, S.
2015-01-22
In this study, we implement and examine three excited state solvent models in time-dependent self-consistent field methods using a consistent formalism which unambiguously shows their relationship. These are the linear response, state specific, and vertical excitation solvent models. Their effects on energies calculated with the equivalent of COSMO/CIS/AM1 are given for a set of test molecules with varying excited state charge transfer character. The resulting solvent effects are explained qualitatively using a dipole approximation. It is shown that the fundamental differences between these solvent models are reflected by the character of the calculated excitations.
Quantum work statistics of charged Dirac particles in time-dependent fields
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Deffner, Sebastian; Saxena, Avadh
2015-09-28
The quantum Jarzynski equality is an important theorem of modern quantum thermodynamics. We show that the Jarzynski equality readily generalizes to relativistic quantum mechanics described by the Dirac equation. After establishing the conceptual framework we solve a pedagogical, yet experimentally relevant, system analytically. As a main result we obtain the exact quantum work distributions for charged particles traveling through a time-dependent vector potential evolving under Schrdinger as well as under Dirac dynamics, and for which the Jarzynski equality is verified. Thus, special emphasis is put on the conceptual and technical subtleties arising from relativistic quantum mechanics.
Momentum space iterative solution of the time-dependent Schrödinger equation
Kiss, G. Zs.; Borbély, S.; Nagy, L.
2013-11-13
We present a novel approach, the iterative solution of the time-dependent Schrödinger equation (iTDSE model), for the investigation of atomic systems interacting with external laser fields. This model is the extension of the momentum-space strong-field approximation (MSSFA) [1], in which the Coulomb potential was considered only as a first order perturbation. In the iTDSE approach higher order terms were gradually introduced until convergence was achieved. Benchmark calculations were done on the hydrogen atom, and the obtained results were compared to the direct numerical solution [2].
Fu, Jin; Wu, Sheng; Li, Hong; Petzold, Linda R.
2014-10-01
The inhomogeneous stochastic simulation algorithm (ISSA) is a fundamental method for spatial stochastic simulation. However, when diffusion events occur more frequently than reaction events, simulating the diffusion events by ISSA is quite costly. To reduce this cost, we propose to use the time dependent propensity function in each step. In this way we can avoid simulating individual diffusion events, and use the time interval between two adjacent reaction events as the simulation stepsize. We demonstrate that the new algorithm can achieve orders of magnitude efficiency gains over widely-used exact algorithms, scales well with increasing grid resolution, and maintains a high level of accuracy.
Quantum work statistics of charged Dirac particles in time-dependent fields
Deffner, Sebastian; Saxena, Avadh
2015-09-28
The quantum Jarzynski equality is an important theorem of modern quantum thermodynamics. We show that the Jarzynski equality readily generalizes to relativistic quantum mechanics described by the Dirac equation. After establishing the conceptual framework we solve a pedagogical, yet experimentally relevant, system analytically. As a main result we obtain the exact quantum work distributions for charged particles traveling through a time-dependent vector potential evolving under Schrödinger as well as under Dirac dynamics, and for which the Jarzynski equality is verified. Thus, special emphasis is put on the conceptual and technical subtleties arising from relativistic quantum mechanics.
The tunneling solutions of the time-dependent Schroedinger equation for a square-potential barrier
Elci, A.; Hjalmarson, H. P.
2009-10-15
The exact tunneling solutions of the time-dependent Schroedinger equation with a square-potential barrier are derived using the continuous symmetry group G{sub S} for the partial differential equation. The infinitesimal generators and the elements for G{sub S} are represented and derived in the jet space. There exist six classes of wave functions. The representative (canonical) wave functions for the classes are labeled by the eigenvalue sets, whose elements arise partially from the reducibility of a Lie subgroup G{sub LS} of G{sub S} and partially from the separation of variables. Each eigenvalue set provides two or more time scales for the wave function. The ratio of two time scales can act as the duration of an intrinsic clock for the particle motion. The exact solutions of the time-dependent Schroedinger equation presented here can produce tunneling currents that are orders of magnitude larger than those produced by the energy eigenfunctions. The exact solutions show that tunneling current can be quantized under appropriate boundary conditions and tunneling probability can be affected by a transverse acceleration.
2007 Time_Dependent Density-Functional Therory (July 15-20, 2007 Colby College, Maine)
Ullrich Carsten Nancy Ryan Gray
2008-09-19
Time-dependent density-functional theory (TDDFT) provides an efficient, elegant, and formally exact way of describing the dynamics of interacting many-body quantum systems, circumventing the need for solving the full time-dependent Schroedinger equation. In the 20 years since it was first rigorously established in 1984, the field of TDDFT has made rapid and significant advances both formally as well as in terms of successful applications in chemistry, physics and materials science. Today, TDDFT has become the method of choice for calculating excitation energies of complex molecules, and is becoming increasingly popular for describing optical and spectroscopic properties of a variety of materials such as bulk solids, clusters and nanostructures. Other growing areas of applications of TDDFT are nonlinear dynamics of strongly excited electronic systems and molecular electronics. The purpose and scope of this Gordon Research Conference is to provide a platform for discussing the current state of the art of the rapidly progressing, highly interdisciplinary field of TDDFT, to identify and debate open questions, and to point out new promising research directions. The conference will bring together experts with a diverse background in chemistry, physics, and materials science.
Current density partitioning in time-dependent current density functional theory
Mosquera, Martín A.; Wasserman, Adam; Department of Physics, Purdue University, West Lafayette, Indiana 47907
2014-05-14
We adapt time-dependent current density functional theory to allow for a fragment-based solution of the many-electron problem of molecules in the presence of time-dependent electric and magnetic fields. Regarding a molecule as a set of non-interacting subsystems that individually evolve under the influence of an auxiliary external electromagnetic vector-scalar potential pair, the partition 4-potential, we show that there are one-to-one mappings between this auxiliary potential, a sharply-defined set of fragment current densities, and the total current density of the system. The partition electromagnetic (EM) 4-potential is expressed in terms of the real EM 4-potential of the system and a gluing EM 4-potential that accounts for exchange-correlation effects and mutual interaction forces between fragments that are required to yield the correct electron dynamics. We prove the zero-force theorem for the fragmented system, establish a variational formulation in terms of action functionals, and provide a simple illustration for a charged particle in a ring.
Time-dependent models for dark matter at the galactic center
Bertone, Gianfranco
2005-11-15
The prospects for indirect detection of dark matter at the galactic center with {gamma}-ray experiments like the space telescope GLAST, and air Cherenkov telescopes like HESS, CANGAROO, MAGIC and VERITAS depend sensitively on the mass profile within the inner parsec. We calculate the distribution of dark matter on subparsec scales by integrating the time-dependent Fokker-Planck equation, including the effects of self-annihilations, scattering of dark matter particles by stars, and capture in the supermassive black hole. We consider a variety of initial dark matter distributions, including models with very high densities ('spikes') near the black hole, and models with 'adiabatic compression' of the baryons. The annihilation signal after 10{sup 10} yr is found to be substantially reduced from its initial value, but in dark matter models with an initial spike, order-of-magnitude enhancements can persist compared with the rate in spike-free models.
Time-dependent Protein-directed Growth of Gold Nanoparticles within a Single Crystal of Lysozyme
H Wei; Z Wang; J Zhang; S House; Y Gao; L Yang; H Robinson; L Tan; H Xing; C Hou
2011-12-31
Gold nanoparticles are useful in biomedical applications due to their distinct optical properties and high chemical stability. Reports of the biogenic formation of gold colloids from gold complexes has also led to an increased level of interest in the biomineralization of gold. However, the mechanism responsible for biomolecule-directed gold nanoparticle formation remains unclear due to the lack of structural information about biological systems and the fast kinetics of biomimetic chemical systems in solution. Here we show that intact single crystals of lysozyme can be used to study the time-dependent, protein-directed growth of gold nanoparticles. The protein crystals slow down the growth of the gold nanoparticles, allowing detailed kinetic studies to be carried out, and permit a three-dimensional structural characterization that would be difficult to achieve in solution. Furthermore, we show that additional chemical species can be used to fine-tune the growth rate of the gold nanoparticles.
Proposed SPAR Modeling Method for Quantifying Time Dependent Station Blackout Cut Sets
John A. Schroeder
2010-06-01
Abstract: The U.S. Nuclear Regulatory Commissions (USNRCs) Standardized Plant Analysis Risk (SPAR) models and industry risk models take similar approaches to analyzing the risk associated with loss of offsite power and station blackout (LOOP/SBO) events at nuclear reactor plants. In both SPAR models and industry models, core damage risk resulting from a LOOP/SBO event is analyzed using a combination of event trees and fault trees that produce cut sets that are, in turn, quantified to obtain a numerical estimate of the resulting core damage risk. A proposed SPAR method for quantifying the time-dependent cut sets is sometimes referred to as a convolution method. The SPAR method reflects assumptions about the timing of emergency diesel failures, the timing of subsequent attempts at emergency diesel repair, and the timing of core damage that may be different than those often used in industry models. This paper describes the proposed SPAR method.
Time-dependent simulation of prebunched one and two-beam free electron laser
Mirian, N. S.; Maraghechi, B.
2014-04-15
A numerical simulation in one-dimension is conducted to study the slippage effects on prebunched free electron laser. A technique for the simulation of time dependent free electron lasers (FEL) to model the slippage effects is introduced, and the slowly varying envelope approximation in both z and t is used to illustrate the temporal behaviour in the prebunched FEL. Slippage effect on prebunched two-beam FEL is compared with the one-beam modeling. The evaluation of the radiation pulse energy, thermal and phase distribution, and radiation pulse shape in one-beam and two-beam modeling is studied. It was shown that the performance is considerably undermined when the slippage time is comparable to the pulse duration. However, prebunching reduces the slippage. Prebunching also leads to the radiation pulse with a single smooth spike.
Staff, Jan E.; Niebergal, Brian P.; Ouyed, Rachid; Pudritz, Ralph E.; Cai, Kai
2010-10-20
We perform state-of-the-art, three-dimensional, time-dependent simulations of magnetized disk winds, carried out to simulation scales of 60 AU, in order to confront optical Hubble Space Telescope observations of protostellar jets. We 'observe' the optical forbidden line emission produced by shocks within our simulated jets and compare these with actual observations. Our simulations reproduce the rich structure of time-varying jets, including jet rotation far from the source, an inner (up to 400 km s{sup -1}) and outer (less than 100 km s{sup -1}) component of the jet, and jet widths of up to 20 AU in agreement with observed jets. These simulations when compared with the data are able to constrain disk wind models. In particular, models featuring a disk magnetic field with a modest radial spatial variation across the disk are favored.
Chemical dynamics in the gas phase: Time-dependent quantum mechanics of chemical reactions
Gray, S.K.
1993-12-01
A major goal of this research is to obtain an understanding of the molecular reaction dynamics of three and four atom chemical reactions using numerically accurate quantum dynamics. This work involves: (i) the development and/or improvement of accurate quantum mechanical methods for the calculation and analysis of the properties of chemical reactions (e.g., rate constants and product distributions), and (ii) the determination of accurate dynamical results for selected chemical systems, which allow one to compare directly with experiment, determine the reliability of the underlying potential energy surfaces, and test the validity of approximate theories. This research emphasizes the use of recently developed time-dependent quantum mechanical methods, i.e. wave packet methods.
Buehring, W. A.; Samsa, M. E.; Decision and Information Sciences
2008-03-28
The examination of which protective measures are the most appropriate to be implemented in order to prevent, protect against, respond to, and recover from attacks on critical infrastructures and key resources typically involves a comparison of the consequences that could occur when the protective measure is implemented to those that could occur when it is not. This report describes a framework for evaluation that provides some additional capabilities for comparing optional protective measures. It illustrates some potentially important time-dependent factors, such as the implementation rate, that affect the relative pros and cons associated with widespread implementation of protective measures. It presents example results from the use of protective measures, such as detectors and pretrained responders, for an illustrative biological incident. Results show that the choice of an alternative measure can depend on whether or not policy and financial support can be maintained for extended periods of time. Choice of a time horizon greatly influences the comparison of alternatives.
Grey transport acceleration method for time-dependent radiative transfer problems
Larsen, E.
1988-10-01
A new iterative method for solving hte time-dependent multifrequency radiative transfer equations is described. The method is applicable to semi-implicit time discretizations that generate a linear steady-state multifrequency transport problem with pseudo-scattering within each time step. The standard ''lambda'' iteration method is shown to often converge slowly for such problems, and the new grey transport acceleration (GTA) method, based on accelerating the lambda method by employing a grey, or frequency-independent transport equation, is developed. The GTA method is shown, theoretically by an iterative Fourier analysis, and experimentally by numerical calculations, to converge significantly faster than the lambda method. In addition, the GTA method is conceptually simple to implement for general differencing schemes, on either Eulerian or Lagrangian meshes. copyright 1988 Academic Press, Inc.
Talamo, Alberto
2013-05-01
This study presents three numerical algorithms to solve the time dependent neutron transport equation by the method of the characteristics. The algorithms have been developed taking into account delayed neutrons and they have been implemented into the novel MCART code, which solves the neutron transport equation for two-dimensional geometry and an arbitrary number of energy groups. The MCART code uses regular mesh for the representation of the spatial domain, it models up-scattering, and takes advantage of OPENMP and OPENGL algorithms for parallel computing and plotting, respectively. The code has been benchmarked with the multiplication factor results of a Boiling Water Reactor, with the analytical results for a prompt jump transient in an infinite medium, and with PARTISN and TDTORT results for cross section and source transients. The numerical simulations have shown that only two numerical algorithms are stable for small time steps.
Measurement of time-dependent CP violation in B 0 → η'K 0 decays
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Šantelj, L.; Yusa, Y.; Abdesselam, A.; Adachi, I.; Aihara, H.; Al Said, S.; Asner, D. M.; Aulchenko, V.; Aushev, T.; Ayad, R.; et al
2014-10-29
We present a measurement of the time-dependent CP violation parameters in B0 → η'K0 decays. The measurement is based on the full data sample containing 772×106 BB-bar pairs collected at the Υ(4S) resonance using the Belle detector at the KEKB asymmetric-energy e+e- collider. The measured values of the mixing-induced and direct CP violation parameters are: sin 2φ1eff = +0.68 ± 0.07 ± 0.03, Aη'K0 = +0.03 ± 0.05 ± 0.04, where the first uncertainty is statistical and the second is systematic. The values obtained are the most accurate to date. Furthermore, these results are consistent with our previous measurements andmore » with the world-average value of sin 2φ1 measured in B0 → J/ψK0 decays.« less
Studies of time dependence of fields in TEVATRON superconducting dipole magnets
Hanft, R.W.; Brown, B.C.; Herrup, D.A.; Lamm, M.J.; McInturff, A.D.; Syphers, M.J.
1988-08-22
The time variation in the magnetic field of a model Tevatron dipole magnet at constant excitation current has been studied. Variations in symmetry allowed harmonic components over long time ranges show a log t behavior indicative of ''flux creep.'' Both short time range and long time range behavior depend in a detailed way on the excitation history. Similar effects are seen in the remnant fields present in full-scale Tevatron dipoles following current ramping. Both magnitudes and time dependences are observed to depend on details for the ramps, such as ramp rate, flattop duration, and number of ramps. In a few magnets, variations are also seen in symmetry unallowed harmonics. 9 refs., 10 figs.
Two-component hybrid time-dependent density functional theory within the Tamm-Dancoff approximation
Khn, Michael; Weigend, Florian
2015-01-21
We report the implementation of a two-component variant of time-dependent density functional theory (TDDFT) for hybrid functionals that accounts for spin-orbit effects within the Tamm-Dancoff approximation (TDA) for closed-shell systems. The influence of the admixture of Hartree-Fock exchange on excitation energies is investigated for several atoms and diatomic molecules by comparison to numbers for pure density functionals obtained previously [M. Khn and F. Weigend, J. Chem. Theory Comput. 9, 5341 (2013)]. It is further related to changes upon switching to the local density approximation or using the full TDDFT formalism instead of TDA. Efficiency is demonstrated for a comparably large system, Ir(ppy){sub 3} (61 atoms, 1501 basis functions, lowest 10 excited states), which is a prototype molecule for organic light-emitting diodes, due to its spin-forbidden triplet-singlet transition.
Adaptive finite volume methods for time-dependent P.D.E.S.
Ware, J.; Berzins, M.
1995-12-31
The aim of adaptive methods for time-dependent p.d.e.s is to control the numerical error so that it is less than a user-specified tolerance. This error depends on the spatial discretization method, the spatial mesh, the method of time integration and the timestep. The spatial discretization method and positioning of the spatial mesh points should attempt to ensure that the spatial error is controlled to meet the user`s requirements. It is then desirable to integrate the o.d.e. system in time with sufficient accuracy so that the temporal error does not corrupt the spatial accuracy or the reliability of the spatial error estimates. This paper is concerned with the development of a prototype algorithm of this type, based on a cell-centered triangular finite volume scheme, for two space dimensional convection-dominated problems.
Zhu, Hong-Ming; Chen, Jin-Wang; Pan, Xiao-Yin; Sahni, Viraht
2014-01-14
We derive via the interaction “representation” the many-body wave function for harmonically confined electrons in the presence of a magnetostatic field and perturbed by a spatially homogeneous time-dependent electric field—the Generalized Kohn Theorem (GKT) wave function. In the absence of the harmonic confinement – the uniform electron gas – the GKT wave function reduces to the Kohn Theorem wave function. Without the magnetostatic field, the GKT wave function is the Harmonic Potential Theorem wave function. We further prove the validity of the connection between the GKT wave function derived and the system in an accelerated frame of reference. Finally, we provide examples of the application of the GKT wave function.
Electrical analysis of wideband and distributed windows using time-dependent field codes
Shang, C.C.; Caplan, M.; Nickel, H.U.; Thumm, M. |
1993-09-16
Windows, which provide the barrier to maintain the vacuum envelope in a microwave tube, are critical components in high-average-power microwave sources, especially at millimeter wavelengths. As RF power levels approach the 100`s of kWs to 1 MW range (CW), the window assembly experiences severe thermal and mechanical stresses. Depending on the source, the bandwidth of the window may be less than 1 GHz for gyrotron oscillators or up to {approximately}20 GHz for the FOM Institute`s fast-tunable, free-electron-maser. The bandwidth requirements give rise to a number of window configurations where the common goal is locally distributed heat dissipation. In order to better understand the transmission and RF properties of these microwave structures, the authors use detailed time-dependent field solvers.
Time-dependent behavior of a localized electron at a heterojunction boundary of graphene
Jang, Min S.; Kim, Hyungjun; Atwater, Harry A.; Goddard, William A.
2010-01-01
We develop a finite-difference time-domain(FDTD) method for simulating the dynamics of graphene electrons, denoted GraFDTD. We then use GraFDTD to study the temporal behavior of a single localized electron wave packet, showing that it exhibits optical-like dynamics including the Goos–Hänchen effect [F. Goos and H. Hänchen, Ann. Phys.436, 333 (1947)] at a heterojunction, but the behavior is quantitatively different than for electromagnetic waves. This suggests issues that must be addressed in designing graphene-based electronic devices analogous to optical devices. GraFDTD should be useful for studying such complex time-dependent behavior of a quasiparticle in graphene.
Measurement of time-dependent CP violation in B0 ? ?'K0 decays
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
antelj, L.; Yusa, Y.; Abdesselam, A.; Adachi, I.; Aihara, H.; Al Said, S.; Asner, D. M.; Aulchenko, V.; Aushev, T.; Ayad, R.; et al
2014-10-29
We present a measurement of the time-dependent CP violation parameters in B0 ? ?'K0 decays. The measurement is based on the full data sample containing 772106 BB pairs collected at the ?(4S) resonance using the Belle detector at the KEKB asymmetric-energy e+e- collider. The measured values of the mixing-induced and direct CP violation parameters are: sin 2?eff1 = +0.68 0.07 0.03, A?'K0 = +0.03 0.05 0.04, where the first uncertainty is statistical and the second is systematic. The values obtained are the most accurate to date. Furthermore, these results are consistent with our previous measurements andmorewith the world-average value of sin 2?1 measured in B0 ? J/?K0 decays.less
Multi-fluid transport code modeling of time-dependent recycling in ELMy H-mode
Pigarov, A. Yu.; Krasheninnikov, S. I.; Hollmann, E. M.; Rognlien, T. D.; Lasnier, C. J.; Unterberg, E.
2014-06-15
Simulations of a high-confinement-mode (H-mode) tokamak discharge with infrequent giant type-I ELMs are performed by the multi-fluid, multi-species, two-dimensional transport code UEDGE-MB, which incorporates the Macro-Blob approach for intermittent non-diffusive transport due to filamentary coherent structures observed during the Edge Localized Modes (ELMs) and simple time-dependent multi-parametric models for cross-field plasma transport coefficients and working gas inventory in material surfaces. Temporal evolutions of pedestal plasma profiles, divertor recycling, and wall inventory in a sequence of ELMs are studied and compared to the experimental time-dependent data. Short- and long-time-scale variations of the pedestal and divertor plasmas where the ELM is described as a sequence of macro-blobs are discussed. It is shown that the ELM recovery includes the phase of relatively dense and cold post-ELM divertor plasma evolving on a several ms scale, which is set by the transport properties of H-mode barrier. The global gas balance in the discharge is also analyzed. The calculated rates of working gas deposition during each ELM and wall outgassing between ELMs are compared to the ELM particle losses from the pedestal and neutral-beam-injection fueling rate, correspondingly. A sensitivity study of the pedestal and divertor plasmas to model assumptions for gas deposition and release on material surfaces is presented. The performed simulations show that the dynamics of pedestal particle inventory is dominated by the transient intense gas deposition into the wall during each ELM followed by continuous gas release between ELMs at roughly a constant rate.
Multi-fluid transport code modeling of time-dependent recycling in ELMy H-mode
Pigarov, A. Yu.; Krasheninnikov, S. I.; Rognlien, T. D.; Hollmann, E. M.; Lasnier, C. J.; Unterberg, Ezekial A
2014-01-01
Simulations of a high-confinement-mode (H-mode) tokamak discharge with infrequent giant type-I ELMs are performed by the multi-fluid, multi-species, two-dimensional transport code UEDGE-MB, which incorporates the Macro-Blob approach for intermittent non-diffusive transport due to filamentary coherent structures observed during the Edge Localized Modes (ELMs) and simple time-dependent multi-parametric models for cross-field plasma transport coefficients and working gas inventory in material surfaces. Temporal evolutions of pedestal plasma profiles, divertor recycling, and wall inventory in a sequence of ELMs are studied and compared to the experimental time-dependent data. Short- and long-time-scale variations of the pedestal and divertor plasmas where the ELM is described as a sequence of macro-blobs are discussed. It is shown that the ELM recovery includes the phase of relatively dense and cold post-ELM divertor plasma evolving on a several ms scale, which is set by the transport properties of H-mode barrier. The global gas balance in the discharge is also analyzed. The calculated rates of working gas deposition during each ELM and wall outgassing between ELMs are compared to the ELM particle losses from the pedestal and neutral-beam-injection fueling rate, correspondingly. A sensitivity study of the pedestal and divertor plasmas to model assumptions for gas deposition and release on material surfaces is presented. The performed simulations show that the dynamics of pedestal particle inventory is dominated by the transient intense gas deposition into the wall during each ELM followed by continuous gas release between ELMs at roughly a constant rate.
Larson, M.K.; Maleki, H.
1996-12-01
The U.S. Bureau of Mines and Cyprus Shoshone Coal Corp. conducted a study of deformation mechanisms in strata around a longwall gate road system at two sites in an underground coal mine near Hanna, WY. Of particular interest was time-dependent (on the order of several months) roof deformation. Strata above and below the coal seam are very weak, carbonaceous mudstones that have cohesionless planes of weakness oriented along bedding, which dips approximately 8{degrees} to 16{degrees} at N 49{degrees} E. An extensive array of instrumented roof bolts, roof extensometers, and biaxial stressmeters were installed during development mining at these two sites under variable depth of cover [183 to 335 in (600 to 1,100 ft)], seam dip (10{degrees} to 16{degrees}), and top coal thickness [estimated to be 0.30 to 0.61 in (1 to 2 ft) on the downdip side of the entry]. Bolt load and roof deformation histories during the entry development periods were compared. Bolt loads at site 2 were less than those at site 1, which is consistent with the amount of overburden. Deformation magnitudes at site 2 were similar to those at site 1, but deformation rates were approximately 44% greater and nearly constant over time, whereas the deformation rates at site 1 decreased exponentially over time. Also, most deformation occurred above the bolt horizon at site two, but extensometer measurements at site 1 showed that a significant amount of deformation occurred within the bolt horizon. The most significant factor that influenced bolt tension and time-dependent strata deformation was strength and Young`s modulus of the strata. Finite-difference models support this conclusion.
TIME-DEPENDENT PHOTOIONIZATION OF GASEOUS NEBULAE: THE PURE HYDROGEN CASE
Garcia, J.; Elhoussieny, E. E.; Bautista, M. A.; Kallman, T. R. E-mail: manuel.bautista@wmich.edu E-mail: timothy.r.kallman@nasa.gov
2013-09-20
We study the problem of time-dependent photoionization of low density gaseous nebulae subjected to sudden changes in the intensity of ionizing radiation. To this end, we write a computer code that solves the full time-dependent energy balance, ionization balance, and radiation transfer equations in a self-consistent fashion for a simplified pure hydrogen case. It is shown that changes in the ionizing radiation yield ionization/thermal fronts that propagate through the cloud, but the propagation times and response times to such fronts vary widely and nonlinearly from the illuminated face of the cloud to the ionization front (IF). IF/thermal fronts are often supersonic, and in slabs initially in pressure equilibrium such fronts yield large pressure imbalances that are likely to produce important dynamical effects in the cloud. Further, we studied the case of periodic variations in the ionizing flux. It is found that the physical conditions of the plasma have complex behaviors that differ from any steady-state solution. Moreover, even the time average of ionization and temperature is different from any steady-state case. This time average is characterized by overionization and a broader IF with respect to the steady-state solution for a mean value of the radiation flux. Around the time average of physical conditions there is a large dispersion in instantaneous conditions, particularly across the IF, which increases with the period of radiation flux variations. Moreover, the variations in physical conditions are asynchronous along the slab due to the combination of nonlinear propagation times for thermal fronts/IFs and equilibration times.
Representing the thermal state in time-dependent density functional theory
Modine, N. A.; Hatcher, R. M.
2015-05-28
Classical molecular dynamics (MD) provides a powerful and widely used approach to determining thermodynamic properties by integrating the classical equations of motion of a system of atoms. Time-Dependent Density Functional Theory (TDDFT) provides a powerful and increasingly useful approach to integrating the quantum equations of motion for a system of electrons. TDDFT efficiently captures the unitary evolution of a many-electron state by mapping the system into a fictitious non-interacting system. In analogy to MD, one could imagine obtaining the thermodynamic properties of an electronic system from a TDDFT simulation in which the electrons are excited from their ground state by a time-dependent potential and then allowed to evolve freely in time while statistical data are captured from periodic snapshots of the system. For a variety of systems (e.g., many metals), the electrons reach an effective state of internal equilibrium due to electron-electron interactions on a time scale that is short compared to electron-phonon equilibration. During the initial time-evolution of such systems following electronic excitation, electron-phonon interactions should be negligible, and therefore, TDDFT should successfully capture the internal thermalization of the electrons. However, it is unclear how TDDFT represents the resulting thermal state. In particular, the thermal state is usually represented in quantum statistical mechanics as a mixed state, while the occupations of the TDDFT wavefunctions are fixed by the initial state in TDDFT. We work to address this puzzle by (A) reformulating quantum statistical mechanics so that thermodynamic expectations can be obtained as an unweighted average over a set of many-body pure states and (B) constructing a family of non-interacting (single determinant) TDDFT states that approximate the required many-body states for the canonical ensemble.
Two-electron Rabi oscillations in real-time time-dependent density-functional theory
Habenicht, Bradley F.; Tani, Noriyuki P.; Provorse, Makenzie R.; Isborn, Christine M.
2014-11-14
We investigate the Rabi oscillations of electrons excited by an applied electric field in several simple molecular systems using time-dependent configuration interaction (TDCI) and real-time time-dependent density-functional theory (RT-TDDFT) dynamics. While the TDCI simulations exhibit the expected single-electron Rabi oscillations at a single resonant electric field frequency, Rabi oscillations in the RT-TDDFT simulations are a two-electron process. The existence of two-electron Rabi oscillations is determined both by full population inversion between field-free molecular orbitals and the behavior of the instantaneous dipole moment during the simulations. Furthermore, the Rabi oscillations in RT-TDDFT are subject to an intensity threshold of the electric field, below which Rabi oscillations do not occur and above which the two-electron Rabi oscillations occur at a broad range of frequencies. It is also shown that at field intensities near the threshold intensity, the field frequency predicted to induce Rabi oscillations by linear response TDDFT only produces detuned Rabi oscillations. Instead, the field frequency that yields the full two-electron population inversion and Rabi oscillation behavior is shown to be the average of single-electron transition frequencies from the ground S{sub 0} state and the doubly-excited S{sub 2} state. The behavior of the two-electron Rabi oscillations is rationalized via two possible models. The first model is a multi-photon process that results from the electric field interacting with the three level system such that three level Rabi oscillations may occur. The second model suggests that the mean-field nature of RT-TDDFT induces paired electron propagation.
Spatial resolution of diffraction tomography
Dickens, T.A.; Winbow, G.A.
1997-01-01
Diffraction tomography is an imaging technique applicable to crosshole seismic data and aimed at achieving optimal spatial resolution away from the borehole. In principle the method can form acoustic images equivalent to extending acoustic well logs away from the wellbore and into the formation with a spatial resolution less than one wavelength of the radiation employed to gather the crosshole data. This paper reports on the capability of diffraction tomography to produce high-resolution reconstructions of simple targets from limited-view-angle data. The goal is to quantify the resolution and velocity-reconstruction capability of diffraction tomography with realistic source{endash}receiver geometries. Simple targets (disks and low-contrast sequences of layers) are used for this study. The scattering from these targets can be calculated without approximation, making them ideal test cases for the algorithm. The resolution capability of diffraction tomography is determined to be on the order of one wavelength for several experimental geometries. It is shown that the image-formation characteristics of diffraction tomography, in terms of its ability to determine object boundaries and velocities, are closely related to the experimental geometry. Reflection and vertical seismic profiling (VSP) experiments tend to reproduce boundaries well, while crosshole experiments give the best overall reconstruction of both target boundaries and velocity. The quantitative accuracy of the velocity reconstruction depends upon the match between the spatial-frequency content of the object and the spatial-frequency response of the algorithm. For some targets, the velocity cannot be correctly reproduced from limited-view-angle data. {copyright} {ital 1997 Acoustical Society of America.}
Boyarinov, V. F.; Kondrushin, A. E.; Fomichenko, P. A.
2013-07-01
Time-dependent equations of the Surface Harmonics Method (SHM) have been derived from the time-dependent neutron transport equation with explicit representation of delayed neutrons for solving the two-dimensional time-dependent problems. These equations have been realized in the SUHAM-TD code. The TWIGL benchmark problem has been used for verification of the SUHAM-TD code. The results of the study showed that computational costs required to achieve necessary accuracy of the solution can be an order of magnitude less than with the use of the conventional finite difference method (FDM). (authors)
Green, M.A.; Cook, N.G.W.; McEvilly, T.V.; Majer, E.L.; Witherspoon, P.A.
1987-04-20
Apparatus is described for placement in a borehole in the earth, which enables the generation of closely controlled seismic waves from the borehole. Pure torsional shear waves are generated by an apparatus which includes a stator element fixed to the borehole walls and a rotor element which is electrically driven to rapidly oscillate on the stator element to cause reaction forces transmitted through the borehole walls to the surrounding earth. Longitudinal shear waves are generated by an armature that is driven to rapidly oscillate along the axis of the borehole, to cause reaction forces transmitted to the surrounding earth. Pressure waves are generated by electrically driving pistons that press against opposite ends of a hydraulic reservoir that fills the borehole. High power is generated by energizing the elements for more than about one minute. 9 figs.
Representing the thermal state in time-dependent density functional theory
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Modine, N. A.; Hatcher, R. M.
2015-05-28
Classical molecular dynamics (MD) provides a powerful and widely used approach to determining thermodynamic properties by integrating the classical equations of motion of a system of atoms. Time-Dependent Density Functional Theory (TDDFT) provides a powerful and increasingly useful approach to integrating the quantum equations of motion for a system of electrons. TDDFT efficiently captures the unitary evolution of a many-electron state by mapping the system into a fictitious non-interacting system. In analogy to MD, one could imagine obtaining the thermodynamic properties of an electronic system from a TDDFT simulation in which the electrons are excited from their ground state bymore » a time-dependent potential and then allowed to evolve freely in time while statistical data are captured from periodic snapshots of the system. For a variety of systems (e.g., many metals), the electrons reach an effective state of internal equilibrium due to electron-electron interactions on a time scale that is short compared to electron-phonon equilibration. During the initial time-evolution of such systems following electronic excitation, electron-phonon interactions should be negligible, and therefore, TDDFT should successfully capture the internal thermalization of the electrons. However, it is unclear how TDDFT represents the resulting thermal state. In particular, the thermal state is usually represented in quantum statistical mechanics as a mixed state, while the occupations of the TDDFT wave functions are fixed by the initial state in TDDFT. Two key questions involve (1) reformulating quantum statistical mechanics so that thermodynamic expectations can be obtained as an unweighted average over a set of many-body pure states and (2) constructing a family of non-interacting (single determinant) TDDFT states that approximate the required many-body states for the canonical ensemble. In Section II, we will address these questions by first demonstrating that thermodynamic
Representing the thermal state in time-dependent density functional theory
Modine, N. A.; Hatcher, R. M.
2015-05-28
Classical molecular dynamics (MD) provides a powerful and widely used approach to determining thermodynamic properties by integrating the classical equations of motion of a system of atoms. Time-Dependent Density Functional Theory (TDDFT) provides a powerful and increasingly useful approach to integrating the quantum equations of motion for a system of electrons. TDDFT efficiently captures the unitary evolution of a many-electron state by mapping the system into a fictitious non-interacting system. In analogy to MD, one could imagine obtaining the thermodynamic properties of an electronic system from a TDDFT simulation in which the electrons are excited from their ground state by a time-dependent potential and then allowed to evolve freely in time while statistical data are captured from periodic snapshots of the system. For a variety of systems (e.g., many metals), the electrons reach an effective state of internal equilibrium due to electron-electron interactions on a time scale that is short compared to electron-phonon equilibration. During the initial time-evolution of such systems following electronic excitation, electron-phonon interactions should be negligible, and therefore, TDDFT should successfully capture the internal thermalization of the electrons. However, it is unclear how TDDFT represents the resulting thermal state. In particular, the thermal state is usually represented in quantum statistical mechanics as a mixed state, while the occupations of the TDDFT wave functions are fixed by the initial state in TDDFT. Two key questions involve (1) reformulating quantum statistical mechanics so that thermodynamic expectations can be obtained as an unweighted average over a set of many-body pure states and (2) constructing a family of non-interacting (single determinant) TDDFT states that approximate the required many-body states for the canonical ensemble. In Section II, we will address these questions by first demonstrating that thermodynamic expectations
Crommentuijn, T.; Doodeman, C.J.A.M.; Doornekamp, A.; Pol, J.J.C. van der; Bedaux, J.J.M.; Gestel, C.A.M. van )
1994-11-01
Time-dependent toxicity in bioassays is usually explained in terms of uptake and elimination kinetics of the toxicant. By comparing different species with essentially different accumulation kinetics, a firm test of this concept may be made. This article compares the sensitivity of six soil arthropods, the collembolans Orchesella cincta and Tomocerus minor, the oribatid mite Platynothrus peltifer, the isopods Porcellio scaber and Oniscus asellus, and the diplopod Cylindroiulus britannicus, when exposed to cadmium in the food. Survival was determined at various time intervals; accumulation of cadmium in the animals was measured at one time interval. Kinetic-based toxicity models were fitted to the data, and estimates were obtained for lethal body concentration, uptake rate constant, elimination rate constant, and ultimate LC50. Two different accumulation patterns could be discerned; these were correlated with time-survival relationships. One, species that have the possibility to eliminate cadmium will reach an equilibrium for the internal concentration and also an ultimate LC50. Two, species that are unable to eliminate cadmium but store it in the body will have an ultimate LC50 equal to zero. For these species the time in which the lethal body concentration is reached is more important. Taxonomically related species appeared to have comparable accumulation patterns, but lethal body concentrations differed. It is concluded that knowledge of the accumulation pattern is indispensable for the evaluation of species' sensitivities to toxicants.
Measurement of Time-Dependent CP Asymmetries in B0 to D(*)D Decays
Aubert, B.; Barate, R.; Boutigny, D.; Couderc, F.; Karyotakis, Y.; Lees, J.P.; Poireau, V.; Tisserand, V.; Zghiche, A.; Grauges, E.; Palano, A.; Pappagallo, M.; Pompili, A.; Chen, J.C.; Qi, N.D.; Rong, G.; Wang, P.; Zhu, Y.S.; Eigen, G.; Ofte, I.; Stugu, B. /Bergen U. /LBL, Berkeley /UC, Berkeley /Birmingham U. /Ruhr U., Bochum /Bristol U. /British Columbia U. /Brunel U. /Novosibirsk, IYF /UC, Irvine /UCLA /UC, Riverside /UC, San Diego /UC, Santa Barbara /UC, Santa Cruz /Caltech /Cincinnati U. /Colorado U. /Colorado State U. /Dortmund U. /Dresden, Tech. U. /Ecole Polytechnique /Edinburgh U. /Ferrara U. /INFN, Ferrara /Frascati /Genoa U. /INFN, Genoa /Harvard U. /Heidelberg U. /Imperial Coll., London /Iowa U. /Iowa State U. /Orsay, LAL /LLNL, Livermore /Liverpool U. /Queen Mary, U. of London /Royal Holloway, U. of London /Louisville U. /Manchester U. /Maryland U. /Massachusetts U., Amherst /MIT, LNS /McGill U. /Milan U. /INFN, Milan /Mississippi U. /Montreal U. /Mt. Holyoke Coll. /Naples U. /INFN, Naples /NIKHEF, Amsterdam /Notre Dame U. /Ohio State U. /Oregon U. /Padua U. /INFN, Padua /Paris U., VI-VII /Pennsylvania U. /Perugia U. /INFN, Perugia /Pisa U. /INFN, Pisa /Prairie View A-M /Princeton U. /Rome U. /INFN, Rome /Rostock U. /Rutherford /DAPNIA, Saclay /South Carolina U. /SLAC /Oregon U. /SLAC /SLAC /Stanford U., Phys. Dept. /SUNY, Stony Brook /Tennessee U. /Texas U. /Texas U., Dallas /Turin U. /INFN, Turin /Trieste U. /INFN, Trieste /Valencia U., IFIC /Vanderbilt U. /Victoria U. /Warwick U. /Wisconsin U., Madison /Yale U.
2005-06-24
The authors present a first measurement of CP asymmetries in neutral B decays to D{sup +}D{sup -}, and updated CP asymmetry measurements in decays to D*{sup +}D{sup -} and D*{sup -}D{sup +}. They use fully-reconstructed decays collected in a data sample of (232 {+-} 3) x 10{sup 6} {Upsilon}(4S) {yields} B{bar B} events in the BABAR detector at the PEp-II asymmetric-energy B Factory at SLAC. they determine the time-dependent asymmetry parameters to be S{sub D*{sup +}D{sup -}} = -0.54 {+-} 0.35 {+-} 0.07, C{sub D*{sup +}D{sup -}} = 0.09 {+-} 0.25 {+-} 0.06, S{sub D*{sup -}D{sup +}} = -0.29 {+-} 0.33 {+-} 0.07, C{sub D*{sup -}D{sup +}} = 0.17 {+-} 0.24 {+-} 0.04, S{sub D{sup +}D{sup -}} = -0.29 {+-} 0.63 {+-} 0.06, and C{sub D{sup +}D{sup -}} = 0.11 {+-} 0.35 {+-} 0.06, where in each case the first error is statistical and the second error is systematic.
A pseudospectral matrix method for time-dependent tensor fields on a spherical shell
Brügmann, Bernd
2013-02-15
We construct a pseudospectral method for the solution of time-dependent, non-linear partial differential equations on a three-dimensional spherical shell. The problem we address is the treatment of tensor fields on the sphere. As a test case we consider the evolution of a single black hole in numerical general relativity. A natural strategy would be the expansion in tensor spherical harmonics in spherical coordinates. Instead, we consider the simpler and potentially more efficient possibility of a double Fourier expansion on the sphere for tensors in Cartesian coordinates. As usual for the double Fourier method, we employ a filter to address time-step limitations and certain stability issues. We find that a tensor filter based on spin-weighted spherical harmonics is successful, while two simplified, non-spin-weighted filters do not lead to stable evolutions. The derivatives and the filter are implemented by matrix multiplication for efficiency. A key technical point is the construction of a matrix multiplication method for the spin-weighted spherical harmonic filter. As example for the efficient parallelization of the double Fourier, spin-weighted filter method we discuss an implementation on a GPU, which achieves a speed-up of up to a factor of 20 compared to a single core CPU implementation.
Probability of loss of assured safety in systems with multiple time-dependent failure modes.
Helton, Jon Craig; Pilch, Martin M.; Sallaberry, Cedric M.
2012-09-01
Weak link (WL)/strong link (SL) systems are important parts of the overall operational design of high-consequence systems. In such designs, the SL system is very robust and is intended to permit operation of the entire system under, and only under, intended conditions. In contrast, the WL system is intended to fail in a predictable and irreversible manner under accident conditions and render the entire system inoperable before an accidental operation of the SL system. The likelihood that the WL system will fail to deactivate the entire system before the SL system fails (i.e., degrades into a configuration that could allow an accidental operation of the entire system) is referred to as probability of loss of assured safety (PLOAS). Representations for PLOAS for situations in which both link physical properties and link failure properties are time-dependent are derived and numerically evaluated for a variety of WL/SL configurations, including PLOAS defined by (i) failure of all SLs before failure of any WL, (ii) failure of any SL before failure of any WL, (iii) failure of all SLs before failure of all WLs, and (iv) failure of any SL before failure of all WLs. The effects of aleatory uncertainty and epistemic uncertainty in the definition and numerical evaluation of PLOAS are considered.
Time-dependent models for blazar emission with the second-order Fermi acceleration
Asano, Katsuaki; Takahara, Fumio; Toma, Kenji; Kusunose, Masaaki; Kakuwa, Jun
2014-01-01
The second-order Fermi acceleration (Fermi-II) driven by turbulence may be responsible for the electron acceleration in blazar jets. We test this model with time-dependent simulations. The hard electron spectrum predicted by the Fermi-II process agrees with the hard photon spectrum of 1ES 1101–232. For other blazars that show softer spectra, the Fermi-II model requires radial evolution of the electron injection rate and/or diffusion coefficient in the outflow. Such evolutions can yield a curved electron spectrum, which can reproduce the synchrotron spectrum of Mrk 421 from the radio to the X-ray regime. The photon spectrum in the GeV energy range of Mrk 421 is hard to fit with a synchrotron self-Compton model. However, if we introduce an external radio photon field with a luminosity of 4.9 × 10{sup 38} erg s{sup –1}, GeV photons are successfully produced via inverse Compton scattering. The temporal variability of the diffusion coefficient or injection rate causes flare emission. The observed synchronicity of X-ray and TeV flares implies a decrease of the magnetic field in the flaring source region.
Relativistic Coulomb excitation within the time dependent superfluid local density approximation
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Stetcu, I.; Bertulani, C. A.; Bulgac, A.; Magierski, P.; Roche, K. J.
2015-01-06
Within the framework of the unrestricted time-dependent density functional theory, we present for the first time an analysis of the relativistic Coulomb excitation of the heavy deformed open shell nucleus 238U. The approach is based on the superfluid local density approximation formulated on a spatial lattice that can take into account coupling to the continuum, enabling self-consistent studies of superfluid dynamics of any nuclear shape. We compute the energy deposited in the target nucleus as a function of the impact parameter, finding it to be significantly larger than the estimate using the Goldhaber-Teller model. The isovector giant dipole resonance, themore » dipole pygmy resonance, and giant quadrupole modes are excited during the process. As a result, the one-body dissipation of collective dipole modes is shown to lead a damping width Γ↓≈0.4 MeV and the number of preequilibrium neutrons emitted has been quantified.« less
Online Support Vector Regression with Varying Parameters for Time-Dependent Data
Omitaomu, Olufemi A; Jeong, Myong K; Badiru, Adedeji B
2011-01-01
Support vector regression (SVR) is a machine learning technique that continues to receive interest in several domains including manufacturing, engineering, and medicine. In order to extend its application to problems in which datasets arrive constantly and in which batch processing of the datasets is infeasible or expensive, an accurate online support vector regression (AOSVR) technique was proposed. The AOSVR technique efficiently updates a trained SVR function whenever a sample is added to or removed from the training set without retraining the entire training data. However, the AOSVR technique assumes that the new samples and the training samples are of the same characteristics; hence, the same value of SVR parameters is used for training and prediction. This assumption is not applicable to data samples that are inherently noisy and non-stationary such as sensor data. As a result, we propose Accurate On-line Support Vector Regression with Varying Parameters (AOSVR-VP) that uses varying SVR parameters rather than fixed SVR parameters, and hence accounts for the variability that may exist in the samples. To accomplish this objective, we also propose a generalized weight function to automatically update the weights of SVR parameters in on-line monitoring applications. The proposed function allows for lower and upper bounds for SVR parameters. We tested our proposed approach and compared results with the conventional AOSVR approach using two benchmark time series data and sensor data from nuclear power plant. The results show that using varying SVR parameters is more applicable to time dependent data.
Neutrino and cosmic-ray release from gamma-ray bursts: Time-dependent simulations
Asano, Katsuaki [Institute for Cosmic Ray Research, The University of Tokyo, 5-1-5 Kashiwanoha, Kashiwa, Chiba 277-8582 (Japan); Mszros, Peter, E-mail: asanok@icrr.u-tokyo.ac.jp, E-mail: nnp@psu.edu [Department of Astronomy and Astrophysics, Department of Physics, Center for Particle and Gravitational Astrophysics, Pennsylvania State University, University Park, PA 16802 (United States)
2014-04-10
We revisit the neutrino and ultra-high-energy cosmic-ray (UHECR) production from gamma-ray bursts (GRBs) with time-dependent simulations for the proton-induced cascades. This method can generate self-consistent photon, neutrino, and escaped neutron spectra. To obtain the integrated background spectra, we take into account the distributions of the burst luminosity and pulse duration timescale. A benchmark case with standard GRB luminosity function, a bulk Lorentz factor ? = 300, and a proton to gamma-ray luminosity fraction f{sub p} = 10 is consistent with both the neutrino upper limits and the observed UHECR intensity at ?10{sup 20} eV, while requiring a different type of UHECR source at the ankle. For the benchmark case, the GRBs in the bright end of the luminosity function, which contribute most of the neutrinos, have their photon spectrum substantially distorted by secondary photons. Such bright GRBs are few in number, and reducing their f{sub p} eliminates the distortion and reduces the neutrino production. Even if we neglect the contribution of the brightest GRBs, the UHECR production rate at energies corresponding to the Greisen-Zatsepin-Kuzmin limit is almost unchanged. These nominal GRB models, especially with L {sub iso} ? 10{sup 53} erg s{sup 1}, appear to meet the current constraints as far as being candidate UHECR sources above the ankle energy.
Photon and neutrino spectra of time-dependent photospheric models of gamma-ray bursts
Asano, K. [Institute for Cosmic Ray Research, The University of Tokyo, 5-1-5 Kashiwanoha, Kashiwa, Chiba 277-8582 (Japan); Mszros, P., E-mail: asanok@icrr.u-tokyo.ac.jp, E-mail: nnp@astro.psu.edu [Department of Astronomy and Astrophysics, Department of Physics, Center for Particle and Gravitational Astrophysics, Pennsylvania State University, University Park, PA 16802 (United States)
2013-09-01
Thermal photons from the photosphere may be the primary source of the observed prompt emission of gamma-ray bursts (GRBs). In order to produce the observed non-thermal spectra, some kind of dissipation mechanism near the photosphere is required. In this paper we numerically simulate the evolution of the photon spectrum in a relativistically expanding shell with a time-dependent numerical code. We consider two basic models. One is a leptonic model, where a dissipation mechanism heats the thermal electrons maintaining their high temperature. The other model involves a cascade process induced by pp(pn)-collisions which produce high-energy electrons, modify the thermal spectrum, and emit neutrinos. The qualitative properties of the photon spectra are mainly determined by the optical depth at which the dissipation mechanism sets in. Too large optical depths lead to a broad and curved spectrum contradicting the observations, while for optical depths smaller than unity the spectral hardness becomes softer than observed. A significant shift of the spectral peak energy to higher energies due to a large energy injection can lead to an overly broad spectral shape. We show ideal parameter ranges for which these models are able to reproduce the observed spectra. For the pn-collision model, the neutrino fluence in the 10100 GeV range is well above the atmospheric neutrino fluence, but its detection is challenging for presently available detectors.
Problem-free time-dependent variational principle for open quantum systems
Joubert-Doriol, Loc; Izmaylov, Artur F.
2015-04-07
Methods of quantum nuclear wave-function dynamics have become very efficient in simulating large isolated systems using the time-dependent variational principle (TDVP). However, a straightforward extension of the TDVP to the density matrix framework gives rise to methods that do not conserve the energy in the isolated system limit and the total system population for open systems where only energy exchange with environment is allowed. These problems arise when the system density is in a mixed state and is simulated using an incomplete basis. Thus, the basis set incompleteness, which is inevitable in practical calculations, creates artificial channels for energy and population dissipation. To overcome this unphysical behavior, we have introduced a constrained Lagrangian formulation of TDVP applied to a non-stochastic open system Schrdinger equation [L. Joubert-Doriol, I. G. Ryabinkin, and A. F. Izmaylov, J. Chem. Phys. 141, 234112 (2014)]. While our formulation can be applied to any variational ansatz for the system density matrix, derivation of working equations and numerical assessment is done within the variational multiconfiguration Gaussian approach for a two-dimensional linear vibronic coupling model system interacting with a harmonic bath.
Katsevich, Alexander J.; Ramm, Alexander G.
1996-01-01
Local tomographic data is used to determine the location and value of a discontinuity between a first internal density of an object and a second density of a region within the object. A beam of radiation is directed in a predetermined pattern through the region of the object containing the discontinuity. Relative attenuation data of the beam is determined within the predetermined pattern having a first data component that includes attenuation data through the region. The relative attenuation data is input to a pseudo-local tomography function, where the difference between the internal density and the pseudo-local tomography function is computed across the discontinuity. The pseudo-local tomography function outputs the location of the discontinuity and the difference in density between the first density and the second density.
Katsevich, A.J.; Ramm, A.G.
1996-07-23
Local tomographic data is used to determine the location and value of a discontinuity between a first internal density of an object and a second density of a region within the object. A beam of radiation is directed in a predetermined pattern through the region of the object containing the discontinuity. Relative attenuation data of the beam is determined within the predetermined pattern having a first data component that includes attenuation data through the region. The relative attenuation data is input to a pseudo-local tomography function, where the difference between the internal density and the pseudo-local tomography function is computed across the discontinuity. The pseudo-local tomography function outputs the location of the discontinuity and the difference in density between the first density and the second density. 7 figs.
Advanced Seismic While Drilling System
Robert Radtke; John Fontenot; David Glowka; Robert Stokes; Jeffery Sutherland; Ron Evans; Jim Musser
2008-06-30
. An APS Turbine Alternator powered the SeismicPULSER{trademark} to produce two Hz frequency peak signals repeated every 20 seconds. Since the ION Geophysical, Inc. (ION) seismic survey surface recording system was designed to detect a minimum downhole signal of three Hz, successful performance was confirmed with a 5.3 Hz recording with the pumps running. The two Hz signal generated by the sparker was modulated with the 3.3 Hz signal produced by the mud pumps to create an intense 5.3 Hz peak frequency signal. The low frequency sparker source is ultimately capable of generating selectable peak frequencies of 1 to 40 Hz with high-frequency spectra content to 10 kHz. The lower frequencies and, perhaps, low-frequency sweeps, are needed to achieve sufficient range and resolution for realtime imaging in deep (15,000 ft+), high-temperature (150 C) wells for (a) geosteering, (b) accurate seismic hole depth, (c) accurate pore pressure determinations ahead of the bit, (d) near wellbore diagnostics with a downhole receiver and wired drill pipe, and (e) reservoir model verification. Furthermore, the pressure of the sparker bubble will disintegrate rock resulting in an increased overall rates of penetration. Other applications for the SeismicPULSER{trademark} technology are to deploy a low-frequency source for greater range on a wireline for Reverse Vertical Seismic Profiling (RVSP) and Cross-Well Tomography. Commercialization of the technology is being undertaken by first contacting stakeholders to define the value proposition for rig site services utilizing SeismicPULSER{trademark} technologies. Stakeholders include national oil companies, independent oil companies, independents, service companies, and commercial investors. Service companies will introduce a new Drill Bit SWD service for deep HTHP wells. Collaboration will be encouraged between stakeholders in the form of joint industry projects to develop prototype tools and initial field trials. No barriers have been identified
Reep, J. W.; Bradshaw, S. J.; Klimchuk, J. A. E-mail: stephen.bradshaw@rice.edu
2013-02-20
The time dependence of heating in solar active regions can be studied by analyzing the slope of the emission measure distribution coolward of the peak. In a previous study we showed that low-frequency heating can account for 0% to 77% of active region core emission measures. We now turn our attention to heating by a finite succession of impulsive events for which the timescale between events on a single magnetic strand is shorter than the cooling timescale. We refer to this scenario as a 'nanoflare train' and explore a parameter space of heating and coronal loop properties with a hydrodynamic model. Our conclusions are (1) nanoflare trains are consistent with 86% to 100% of observed active region cores when uncertainties in the atomic data are properly accounted for; (2) steeper slopes are found for larger values of the ratio of the train duration {Delta} {sub H} to the post-train cooling and draining timescale {Delta} {sub C}, where {Delta} {sub H} depends on the number of heating events, the event duration and the time interval between successive events ({tau} {sub C}); (3) {tau} {sub C} may be diagnosed from the width of the hot component of the emission measure provided that the temperature bins are much smaller than 0.1 dex; (4) the slope of the emission measure alone is not sufficient to provide information about any timescale associated with heating-the length and density of the heated structure must be measured for {Delta} {sub H} to be uniquely extracted from the ratio {Delta} {sub H}/{Delta} {sub C}.
Green, Michael A.; Cook, Neville G. W.; McEvilly, Thomas V.; Majer, Ernest L.; Witherspoon, Paul A.
1992-01-01
Apparatus is described for placement in a borehole in the earth, which enables the generation of closely controlled seismic waves from the borehole. Pure torsional shear waves are generated by an apparatus which includes a stator element fixed to the borehole walls and a rotor element which is electrically driven to rapidly oscillate on the stator element to cause reaction forces transmitted through the borehole walls to the surrounding earth. Logitudinal shear waves are generated by an armature that is driven to rapidly oscillate along the axis of the borehole relative to a stator that is clamped to the borehole, to cause reaction forces transmitted to the surrounding earth. Pressure waves are generated by electrically driving pistons that press against opposite ends of a hydraulic reservoir that fills the borehole. High power is generated by energizing the elements at a power level that causes heating to over 150.degree. C. within one minute of operation, but energizing the elements for no more than about one minute.
Seismic Design Expectations Report
Broader source: Energy.gov [DOE]
The Seismic Design Expectations Report (SDER) is a tool that assists DOE federal project review teams in evaluating the technical sufficiency of the project seismic design activities prior to...
Seismic intrusion detector system
Hawk, Hervey L.; Hawley, James G.; Portlock, John M.; Scheibner, James E.
1976-01-01
A system for monitoring man-associated seismic movements within a control area including a geophone for generating an electrical signal in response to seismic movement, a bandpass amplifier and threshold detector for eliminating unwanted signals, pulse counting system for counting and storing the number of seismic movements within the area, and a monitoring system operable on command having a variable frequency oscillator generating an audio frequency signal proportional to the number of said seismic movements.
3-D seismic velocity and attenuation structures in the geothermal field
Nugraha, Andri Dian [Global Geophysics Research Group, Faculty of Mining and Petroleum Engineering, Institute of Technology Bandung, Jalan Ganesha No. 10 Bandung, 40132 (Indonesia)] [Global Geophysics Research Group, Faculty of Mining and Petroleum Engineering, Institute of Technology Bandung, Jalan Ganesha No. 10 Bandung, 40132 (Indonesia); Syahputra, Ahmad [Geophyisical Engineering, Faculty of Mining and Petroleum Engineering, Institute of Technology Bandung, Jalan Ganesha No. 10 Bandung, 40132 (Indonesia)] [Geophyisical Engineering, Faculty of Mining and Petroleum Engineering, Institute of Technology Bandung, Jalan Ganesha No. 10 Bandung, 40132 (Indonesia); Fatkhan,; Sule, Rachmat [Applied Geophysics Research Group, Faculty of Mining and Petroleum Engineering, Institute of Technology Bandung, Jalan Ganesha No. 10 Bandung, 40132 (Indonesia)] [Applied Geophysics Research Group, Faculty of Mining and Petroleum Engineering, Institute of Technology Bandung, Jalan Ganesha No. 10 Bandung, 40132 (Indonesia)
2013-09-09
We conducted delay time tomography to determine 3-D seismic velocity structures (Vp, Vs, and Vp/Vs ratio) using micro-seismic events in the geothermal field. The P-and S-wave arrival times of these micro-seismic events have been used as input for the tomographic inversion. Our preliminary seismic velocity results show that the subsurface condition of geothermal field can be fairly delineated the characteristic of reservoir. We then extended our understanding of the subsurface physical properties through determining of attenuation structures (Qp, Qs, and Qs/Qp ratio) using micro-seismic waveform. We combined seismic velocities and attenuation structures to get much better interpretation of the reservoir characteristic. Our preliminary attanuation structures results show reservoir characterization can be more clearly by using the 3-D attenuation model of Qp, Qs, and Qs/Qp ratio combined with 3-D seismic velocity model of Vp, Vs, and Vp/Vs ratio.
Gonalves, W. C.; Sardella, E.; UNESP-Universidade Estadual Paulista, IPMet-Instituto de Pesquisas Meteorolgicas, CEP 17048-699 Bauru, SP ; Becerra, V. F.; Miloevi?, M. V.; Peeters, F. M.; Departamento de Fsica, Universidade Federal do Cear, 60455-900 Fortaleza, Cear
2014-04-15
The time-dependent Ginzburg-Landau formalism for (d + s)-wave superconductors and their representation using auxiliary fields is investigated. By using the link variable method, we then develop suitable discretization of these equations. Numerical simulations are carried out for a mesoscopic superconductor in a homogeneous perpendicular magnetic field which revealed peculiar vortex states.
Idaho National Laboratory (INL) Seismic Initiative | Department...
Office of Environmental Management (EM)
Presentation from the May 2015 Seismic Lessons-Learned Panel Meeting. INL Seismic ... Seismic Soil Structure Interaction (SSI) Method for Developing Non-Linear Seismic SSI ...
Abdel-Khalek, S.; Berrada, K.; Eleuch, H.
2015-10-15
The dynamics of a superconducting (SC) qubit interacting with a field under decoherence with and without time-dependent coupling effect is analyzed. Quantum features like the collapse–revivals for the dynamics of population inversion, sudden birth and sudden death of entanglement, and statistical properties are investigated under the phase damping effect. Analytic results for certain parametric conditions are obtained. We analyze the influence of decoherence on the negativity and Wehrl entropy for different values of the physical parameters. We also explore an interesting relation between the SC-field entanglement and Wehrl entropy behavior during the time evolution. We show that the amount of SC-field entanglement can be enhanced as the field tends to be more classical. The studied model of SC-field system with the time-dependent coupling has high practical importance due to their experimental accessibility which may open new perspectives in different tasks of quantum formation processing.
Seismic Imaging and Monitoring
Huang, Lianjie
2012-07-09
I give an overview of LANL's capability in seismic imaging and monitoring. I present some seismic imaging and monitoring results, including imaging of complex structures, subsalt imaging of Gulf of Mexico, fault/fracture zone imaging for geothermal exploration at the Jemez pueblo, time-lapse imaging of a walkway vertical seismic profiling data for monitoring CO{sub 2} inject at SACROC, and microseismic event locations for monitoring CO{sub 2} injection at Aneth. These examples demonstrate LANL's high-resolution and high-fidelity seismic imaging and monitoring capabilities.
Teleseismic-Seismic Monitoring At Central Nevada Seismic Zone...
Teleseismic-Seismic Monitoring At Central Nevada Seismic Zone Region (Biasi, Et Al., 2008) Jump to: navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity:...
USING MICRO-SEISMICITY AND SEISMIC VELOCITIES TO MAP SUBSURFACE...
some cases, although a significant portion of seismicity remains diffuse and does not cluster into sharply defined structures. The seismic velocity structure reveals heterogeneous...
Using Micro-Seismicity and Seismic Velocities to Map Subsurface...
some cases, although a significant portion of seismicity remains diffuse and does not cluster into sharply defined structures. The seismic velocity structure reveals heterogeneous...
Seismicity Protocol | Department of Energy
Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site
Seismicity Protocol Seismicity Protocol Project objectives: Develop an updated protocolbest engineering practices to address public and industry issues associated with induced ...
Double-Difference Tomography for Sequestration MVA
Westman, Erik
2008-12-31
Analysis of synthetic data was performed to determine the most cost-effective tomographic monitoring system for a geologic carbon sequestration injection site. Double-difference tomographic inversion was performed on 125 synthetic data sets: five stages of CO2 plume growth, five seismic event regions, and five geophone arrays. Each resulting velocity model was compared quantitatively to its respective synthetic velocity model to determine an accuracy value. The results were examined to determine a relationship between cost and accuracy in monitoring, verification, and accounting applications using double-difference tomography. The geophone arrays with widely-varying geophone locations, both laterally and vertically, performed best. Additionally, double difference seismic tomography was performed using travel time data from a carbon sequestration site at the Aneth oil field in southeast Utah as part of a Department of Energy initiative on monitoring, verification, and accounting (MVA) of sequestered CO2. A total of 1,211 seismic events were recorded from a borehole array consisting of 22 geophones. Artificial velocity models were created to determine the ease with which different CO2 plume locations and sizes can be detected. Most likely because of the poor geophone arrangement, a low velocity zone in the Desert Creek reservoir can only be detected when regions of test site containing the highest ray path coverage are considered. MVA accuracy and precision may be improved through the use of a receiver array that provides more comprehensive ray path coverage.
Hong Qin and Ronald C. Davidson
2011-07-19
In a linear trap confining a one-component nonneutral plasma, the external focusing force is a linear function of the configuration coordinates and/or the velocity coordinates. Linear traps include the classical Paul trap and the Penning trap, as well as the newly proposed rotating-radio- frequency traps and the Mobius accelerator. This paper describes a class of self-similar nonlinear solutions of nonneutral plasma in general time-dependent linear focusing devices, with self-consistent electrostatic field. This class of nonlinear solutions includes many known solutions as special cases.
Prado, F. O.; Duzzioni, E. I. [Universidade Federal de Uberlandia, Caixa Postal 593, 38400-902 Uberlandia, Minas Geraisn (Brazil); Almeida, N. G. de [Instituto de Fisica, Universidade Federal de Goias, 74001-970, Goiania, Goias (Brazil); Moussa, M. H. Y. [Instituto de Fisica de Sao Carlos, Universidade de Sao Paulo, Caixa Postal 369, 13560-970 Sao Carlos, Sao Paulo (Brazil); Villas-Boas, C. J. [Departamento de Fisica, Universidade Federal de Sao Carlos, P.O. Box 676, 13565-905 Sao Carlos, Sao Paulo (Brazil)
2011-07-15
In this paper we detail some results advanced in a recent letter [Prado et al., Phys. Rev. Lett. 102, 073008 (2009).] showing how to engineer reservoirs for two-level systems at absolute zero by means of a time-dependent master equation leading to a nonstationary superposition equilibrium state. We also present a general recipe showing how to build nonadiabatic coherent evolutions of a fermionic system interacting with a bosonic mode and investigate the influence of thermal reservoirs at finite temperature on the fidelity of the protected superposition state. Our analytical results are supported by numerical analysis of the full Hamiltonian model.
Zhang, Yu; Biggs, Jason D.; Healion, Daniel; Govind, Niranjan; Mukamel, Shaul
2012-11-21
We report simulations of X-ray absorption near edge structure (XANES), resonant inelastic X-ray scattering (RIXS) and 1D stimulated X-ray Raman spectroscopy (SXRS) signals of cysteine at the oxygen, nitrogen and sulfur K and L2,3 edges. The simulated XANES signals from the restricted window time-dependent density functional theory (REW-TDDFT) and the static exchange (STEX) method are compared with experiments, showing that REW-TDDFT is more accurate and computationally less expensive than STEX. Simulated RIXS and 1D SXRS signals from REW-TDDFT give some insights on the correlation of different excitations in the molecule.
Regnier, D.; Verriere, M.; Dubray, N.; Schunck, N.
2015-11-30
In this study, we describe the software package FELIX that solves the equations of the time-dependent generator coordinate method (TDGCM) in NN-dimensions (N ≥ 1) under the Gaussian overlap approximation. The numerical resolution is based on the Galerkin finite element discretization of the collective space and the Crank–Nicolson scheme for time integration. The TDGCM solver is implemented entirely in C++. Several additional tools written in C++, Python or bash scripting language are also included for convenience. In this paper, the solver is tested with a series of benchmarks calculations. We also demonstrate the ability of our code to handle a realistic calculation of fission dynamics.
Levander, Alan Richard; Zelt, Colin A.
2015-03-17
The work plan for this project was to develop and apply advanced seismic reflection and wide-angle processing and inversion techniques to high resolution seismic data for the shallow subsurface to seismically characterize the shallow subsurface at hazardous waste sites as an aid to containment and cleanup activities. We proposed to continue work on seismic data that we had already acquired under a previous DoE grant, as well as to acquire additional new datasets for analysis. The project successfully developed and/or implemented the use of 3D reflection seismology algorithms, waveform tomography and finite-frequency tomography using compressional and shear waves for high resolution characterization of the shallow subsurface at two waste sites. These two sites have markedly different near-surface structures, groundwater flow patterns, and hazardous waste problems. This is documented in the list of refereed documents, conference proceedings, and Rice graduate theses, listed below.
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Chaplin, Vernon H.; Bellan, Paul M.
2015-12-28
A time-dependent two-fluid model has been developed to understand axial variations in the plasma parameters in a very high density (peak ne~ > 5x1019 m–3) argon inductively coupled discharge in a long 1.1 cm radius tube. The model equations are written in 1D, with radial losses to the tube walls accounted for by the inclusion of effective particle and energy sink terms. The ambipolar diffusion equation and electron energy equation are solved to find the electron density ne(z,t) and temperature Te(z,t), and the populations of the neutral argon 4s metastable, 4s resonant, and 4p excited state manifolds are calculated inmore » order to determine the stepwise ionization rate and calculate radiative energy losses. The model has been validated through comparisons with Langmuir probe ion saturation current measurements; close agreement between the simulated and measured axial plasma density profiles and the initial density rise rate at each location was obtained at pAr = 30-60 mTorr. Lastly, we present detailed results from calculations at 60 mTorr, including the time-dependent electron temperature, excited state populations, and energy budget within and downstream of the radiofrequency (RF) antenna.« less
Measurement of the Time-Dependent CP Asymmetry of Partially Reconstructed B0 to D*+D*- Decays
Lees, J. P.
2012-08-13
We present a new measurement of the time-dependent CP asymmetry of B{sup 0} {yields}D*{sup +}D*{sup -} decays using (471 {+-} 5) million B{bar B} pairs collected with the BABAR detector at the PEP-II B Factory at the SLAC National Accelerator Laboratory. Using the technique of partial reconstruction, we measure the time-dependent CP asymmetry parameters S = -0.34 {+-} 0.12 {+-} 0.05 and C = +0:15 {+-} 0.09 {+-} 0.04. Using the value for the CP-odd fraction R{perpendicular} = 0.158 {+-} 0.028 {+-} 0.006, previously measured by BABAR with fully reconstructed B{sup 0} {yields} D{sup *+}D{sup *-} events, we extract the CP-even components S{sub +} = -0.49 {+-} 0.18 {+-} 0.07 {+-} 0.04 and C{sub +} = +0.15 {+-} 0.09 {+-} 0.04. In each case, the first uncertainty is statistical and the second is systematic; the third uncertainty on S{sub +} is the contribution from the uncertainty on R{perpendicular}. The measured value of the CP-even component S{sub +} is consistent with the value of sin 2{beta} measured in b {yields} (c{bar c})s transitions, and with the Standard Model expectation of small penguin contributions.
Chaplin, Vernon H.; Bellan, Paul M.
2015-12-28
A time-dependent two-fluid model has been developed to understand axial variations in the plasma parameters in a very high density (peak n_{e~} > 5x10^{19} m^{–3}) argon inductively coupled discharge in a long 1.1 cm radius tube. The model equations are written in 1D, with radial losses to the tube walls accounted for by the inclusion of effective particle and energy sink terms. The ambipolar diffusion equation and electron energy equation are solved to find the electron density n_{e}(z,t) and temperature T_{e}(z,t), and the populations of the neutral argon 4s metastable, 4s resonant, and 4p excited state manifolds are calculated in order to determine the stepwise ionization rate and calculate radiative energy losses. The model has been validated through comparisons with Langmuir probe ion saturation current measurements; close agreement between the simulated and measured axial plasma density profiles and the initial density rise rate at each location was obtained at p_{Ar} = 30-60 mTorr. Lastly, we present detailed results from calculations at 60 mTorr, including the time-dependent electron temperature, excited state populations, and energy budget within and downstream of the radiofrequency (RF) antenna.
Seismicity and Reservoir Fracture Characterization
Broader source: Energy.gov [DOE]
Below are the project presentations and respective peer review results for Seismicity and Reservoir Fracture Characterization.
Develpment of a low Cost Method to Estimate the Seismic Signiture...
Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site
... data (database building, seismic interferometry, interpretation) Previous seismic ... data (database building, seismic interferometry, interpretation) Previous seismic ...
Method of migrating seismic records
Ober, Curtis C.; Romero, Louis A.; Ghiglia, Dennis C.
2000-01-01
The present invention provides a method of migrating seismic records that retains the information in the seismic records and allows migration with significant reductions in computing cost. The present invention comprises phase encoding seismic records and combining the encoded seismic records before migration. Phase encoding can minimize the effect of unwanted cross terms while still allowing significant reductions in the cost to migrate a number of seismic records.
Deepwater seismic acquisition technology
Caldwell, J.
1996-09-01
Although truly new technology is not required for successful acquisition of seismic data in deep Gulf of Mexico waters, it is helpful to review some basic aspects of these seismic surveys. Additionally, such surveys are likely to see early use of some emerging new technology which can improve data quality. Because such items as depth imaging, borehole seismic, 4-D and marine 3-component recording were mentioned in the May 1996 issue of World Oil, they are not discussed again here. However, these technologies will also play some role in the deepwater seismic activities. What is covered in this paper are some new considerations for: (1) longer data records needed in deeper water, (2) some pros and cons of very long steamer use, and (3) two new commercial systems for quantifying data quality.
BUILDING 341 Seismic Evaluation
Halle, J.
2015-06-15
The Seismic Evaluation of Building 341 located at Lawrence Livermore National Laboratory in Livermore, California has been completed. The subject building consists of a main building, Increment 1, and two smaller additions; Increments 2 and 3.
Category:Seismic Techniques | Open Energy Information
Seismic Techniques Jump to: navigation, search GEOTHERMAL ENERGYGeothermal Home Geothermalpower.jpg Looking for the Seismic Techniques page? For detailed information on Seismic...
Measurement of the Time-dependent CP Asymmetry inB to D(*)_CP h0 Decays
Aubert, B.
2007-03-14
The authors report a measurement of the time-dependent CP-asymmetry parameters S and C in color-suppressed B{sup 0} {yields} D{sup (*)0}h{sup 0} decays, where h{sup 0} is a {pi}{sup 0}, {eta}, or {omega} meson, and the D{sup 0} decays to one of the CP eigenstates K{sup +}K{sup -}, K{sub S}{sup 0}{pi}{sup 0}, or K{sub S}{sup 0}{omega}. The data sample consists of 383 x 10{sup 6} {Upsilon}(4S) {yields} B{bar B} decays collected with the BABAR detector at the PEP-II asymmetric-energy B-factory at SLAC. The results are S = -0.56 {+-} 0.23 {+-} 0.05 and C = -0.23 {+-} 0.16 {+-} 0.04, where the first error is statistical and the second is systematic.
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Regnier, D.; Verriere, M.; Dubray, N.; Schunck, N.
2015-11-30
In this study, we describe the software package FELIX that solves the equations of the time-dependent generator coordinate method (TDGCM) in NN-dimensions (N ≥ 1) under the Gaussian overlap approximation. The numerical resolution is based on the Galerkin finite element discretization of the collective space and the Crank–Nicolson scheme for time integration. The TDGCM solver is implemented entirely in C++. Several additional tools written in C++, Python or bash scripting language are also included for convenience. In this paper, the solver is tested with a series of benchmarks calculations. We also demonstrate the ability of our code to handle amore » realistic calculation of fission dynamics.« less
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Bahinipati, S.; Trabelsi, K.; Kinoshita, K.; Arinstein, K.; Aulchenko, V.; Aushev, T.; Bakich, A. M.; Balagura, V.; Barberio, E.; Belous, K.; et al
2011-07-05
We report results on time-dependent CP asymmetries in B→D*∓π± decays based on a data sample containing 657×10⁶ BB¯ pairs collected with the Belle detector at the KEKB asymmetric-energy e⁺e⁻ collider at the Υ(4S) resonance. We use a partial reconstruction technique, wherein signal B→D*∓π± events are identified using information only from the fast pion from the B decay and the slow pion from the subsequent decay of the D*∓, where the former (latter) corresponds to D*⁺(D*⁻) final states. We obtain CP violation parameters S⁺=+0.061±0.018 (stat)±0.012 (syst) and S⁻=+0.031±0.019 (stat)±0.015 (syst).
Bjorgaard, J. A.; Velizhanin, K. A.; Tretiak, S.
2015-08-06
This study describes variational energy expressions and analytical excited state energy gradients for time-dependent self-consistent field methods with polarizable solvent effects. Linear response, vertical excitation, and state-specific solventmodels are examined. Enforcing a variational ground stateenergy expression in the state-specific model is found to reduce it to the vertical excitation model. Variational excited state energy expressions are then provided for the linear response and vertical excitation models and analytical gradients are formulated. Using semiempiricalmodel chemistry, the variational expressions are verified by numerical and analytical differentiation with respect to a static external electric field. Lastly, analytical gradients are further tested by performing microcanonical excited state molecular dynamics with p-nitroaniline.
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Bjorgaard, J. A.; Velizhanin, K. A.; Tretiak, S.
2015-08-06
This study describes variational energy expressions and analytical excited state energy gradients for time-dependent self-consistent field methods with polarizable solvent effects. Linear response, vertical excitation, and state-specific solventmodels are examined. Enforcing a variational ground stateenergy expression in the state-specific model is found to reduce it to the vertical excitation model. Variational excited state energy expressions are then provided for the linear response and vertical excitation models and analytical gradients are formulated. Using semiempiricalmodel chemistry, the variational expressions are verified by numerical and analytical differentiation with respect to a static external electric field. Lastly, analytical gradients are further tested by performingmore » microcanonical excited state molecular dynamics with p-nitroaniline.« less
Interaction between O{sub 2} and ZnO films probed by time-dependent second-harmonic generation
Andersen, S. V.; Vandalon, V.; Bosch, R. H. E. C.; Loo, B. W. H. van de; Kessels, W. M. M.; Pedersen, K.
2014-02-03
The interaction between O{sub 2} and ZnO thin films prepared by atomic layer deposition has been investigated by time-dependent second-harmonic generation, by probing the electric field induced by adsorbed oxygen molecules on the surface. The second-harmonic generated signal decays upon laser exposure due to two-photon assisted desorption of O{sub 2}. Blocking and unblocking the laser beam for different time intervals reveals the adsorption rate of O{sub 2} onto ZnO. The results demonstrate that electric field induced second-harmonic generation provides a versatile non-contact probe of the adsorption kinetics of molecules on ZnO thin films.
Wang, RuLin; Zheng, Xiao; Kwok, YanHo; Xie, Hang; Chen, GuanHua; Yam, ChiYung
2015-04-14
Understanding electronic dynamics on material surfaces is fundamentally important for applications including nanoelectronics, inhomogeneous catalysis, and photovoltaics. Practical approaches based on time-dependent density functional theory for open systems have been developed to characterize the dissipative dynamics of electrons in bulk materials. The accuracy and reliability of such approaches depend critically on how the electronic structure and memory effects of surrounding material environment are accounted for. In this work, we develop a novel squared-Lorentzian decomposition scheme, which preserves the positive semi-definiteness of the environment spectral matrix. The resulting electronic dynamics is guaranteed to be both accurate and convergent even in the long-time limit. The long-time stability of electronic dynamics simulation is thus greatly improved within the current decomposition scheme. The validity and usefulness of our new approach are exemplified via two prototypical model systems: quasi-one-dimensional atomic chains and two-dimensional bilayer graphene.
Seismic Consequence Abstraction
M. Gross
2004-10-25
The primary purpose of this model report is to develop abstractions for the response of engineered barrier system (EBS) components to seismic hazards at a geologic repository at Yucca Mountain, Nevada, and to define the methodology for using these abstractions in a seismic scenario class for the Total System Performance Assessment - License Application (TSPA-LA). A secondary purpose of this model report is to provide information for criticality studies related to seismic hazards. The seismic hazards addressed herein are vibratory ground motion, fault displacement, and rockfall due to ground motion. The EBS components are the drip shield, the waste package, and the fuel cladding. The requirements for development of the abstractions and the associated algorithms for the seismic scenario class are defined in ''Technical Work Plan For: Regulatory Integration Modeling of Drift Degradation, Waste Package and Drip Shield Vibratory Motion and Seismic Consequences'' (BSC 2004 [DIRS 171520]). The development of these abstractions will provide a more complete representation of flow into and transport from the EBS under disruptive events. The results from this development will also address portions of integrated subissue ENG2, Mechanical Disruption of Engineered Barriers, including the acceptance criteria for this subissue defined in Section 2.2.1.3.2.3 of the ''Yucca Mountain Review Plan, Final Report'' (NRC 2003 [DIRS 163274]).
Cochrane, Alexander P.; Merrett, Craig G.; Hilton, Harry H.
2014-12-10
The advent of new structural concepts employing composites in primary load carrying aerospace structures in UAVs, MAVs, Boeing 787s, Airbus A380s, etc., necessitates the inclusion of flexibility as well as viscoelasticity in static structural and aero-viscoelastic analyses. Differences and similarities between aeroelasticity and aero-viscoelasticity have been investigated in [2]. An investigation is undertaken as to the dependence and sensitivity of aerodynamic and stability derivatives to elastic and viscoelastic structural flexibility and as to time dependent flight and maneuver velocities. Longitudinal, lateral and directional stabilities are investigated. It has been a well established fact that elastic lifting surfaces are subject to loss of control effectiveness and control reversal at certain flight speeds, which depend on aerodynamic, structural and material properties [5]. Such elastic analyses are extended to linear viscoelastic materials under quasi-static, dynamic, and sudden and gradual loading conditions. In elastic wings one of the critical static parameters is the velocity at which control reversal takes place (V{sub REV}{sup E}). Since elastic formulations constitute viscoelastic initial conditions, viscoelastic reversal may occur at speeds V{sub REV<}{sup ≧}V{sub REV}{sup E}, but furthermore does so in time at 0 < t{sub REV} ≤ ∞. The influence of the twin effects of viscoelastic and elastic materials and of variable flight velocities on longitudinal, lateral, directional and spin stabilities are also investigated. It has been a well established fact that elastic lifting surfaces are subject to loss of control effectiveness and control reversal at certain flight speeds, which depend on aerodynamic, structural and material properties [5]. Such elastic analyses are here extended to linear viscoelastic materials under quasi-static, dynamic, and sudden and gradual loading conditions. In elastic wings the critical parameter is the velocity at
HANFORD DOUBLE SHELL TANK THERMAL AND SEISMIC PROJECT SEISMIC...
Office of Scientific and Technical Information (OSTI)
Title: HANFORD DOUBLE SHELL TANK THERMAL AND SEISMIC PROJECT SEISMIC ANALYSIS OF HANFORD DOUBLE SHELL TANKS M&D Professional Services, Inc. (M&D) is under subcontract to Pacific ...
Time-dependent Dalitz-Plot Analysis of the Charmless Decay B^0 -> K^0S Pi Pi- at BABAR
Ilic, J
2009-10-17
A time-dependent amplitude analysis of B{sup 0} {yields} K{sub S}{sup 0}{pi}{sup +}{pi}{sup -} decays is performed in order to extract the CP violation parameters of f{sub 0}(980)K{sub S}{sup 0} and {rho}{sup 0}(770)K{sub S}{sup 0} and direct CP asymmetries of K*{sup +}(892){pi}{sup -}. The results are obtained from the final BABAR data sample of (465 {+-} 5)10{sup 6} B{bar B} decays, collected with the BABAR detector at the PEP-II asymmetric-energy B factory at SLAC. The time dependent CP asymmetry for f{sub 0}(980)K{sub S}{sup 0} and {rho}{sup 0}(770)K{sub S}{sup 0} are measured to be S(f{sub 0}(980)K{sub S}{sup 0}) = -0.97 {+-} 0.09 {+-} 0.01 {+-} 0.01, and S({rho}{sup 0}(770)K{sub S}{sup 0}) = 0.67 {+-} 0.20 {+-} 0.06 {+-} 0.04, respectively. In decays to K*{sup +}(892){pi}{sup -} the direct CP asymmetry is found to be A{sub CP}(K*{sup {+-}}(892){pi}{sup {-+}}) = -0.18 {+-} 0.10 {+-} 0.04 {+-} 0.00. The relative phases between B{sup 0} {yields} K*{sup +}(892){pi}{sup -} and {bar B}{sup 0} {yields} K*{sup -}(892){pi}{sup +}, relevant for the extraction of the unitarity triangle angle {gamma}, is measured to be {Delta}{phi}(K*(892){pi}) = (34.9 {+-} 23.1 {+-} 7.5 {+-} 4.7){sup o}, where uncertainties are statistical, systematic and model-dependent, respectively. Fit fractions, direct CP asymmetries and the relative phases of different other resonant modes have also been measured. A new method for extracting longitudinal shower development information from longitudinally unsegmented calorimeters is also presented. This method has been implemented as a part of the BABAR final particle identification algorithm. A significant improvement in low momenta muon identification at BABAR is obtained.
Seismic Fracture Characterization Methods for Enhanced Geothermal...
Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site
Seismic Fracture Characterization Methods for Enhanced Geothermal Systems; 2010 Geothermal Technology Program Peer Review Report Seismic Fracture Characterization Methods for ...
Dias Da Silva, Luis G; Heidrich-Meisner, Fabian; Feiguin, Adrian E; Busser, C. A.; Martins, G. B.; Anda, E. V.; Dagotto, Elbio R
2008-01-01
We apply the adaptive time-dependent density-matrix renormalization-group method tDMRG to the study of transport properties of quantum-dot systems connected to metallic leads. Finite-size effects make the usual tDMRG description of the Kondo regime a numerically demanding task. We show that such effects can be attenuated by describing the leads by Wilson chains, in which the hopping matrix elements decay exponentially away from the impurity tn n/2. For a given system size and in the linear-response regime, results for 1 show several improvements over the undamped =1 case: perfect conductance is obtained deeper in the strongly interacting regime and current plateaus remain well defined for longer time scales. Similar improvements were obtained in the finite-bias regime up to bias voltages of the order of the Kondo temperature. These results show that with the proposed modification, the tDMRG characterization of Kondo correlations in the transport properties can be substantially improved, while it turns out to be sufficient to work with much smaller system sizes. We discuss the numerical cost of this approach with respect to the necessary system sizes and the entanglement growth during the time evolution.
Measurement of time-dependent CP violation in B ^{0} ? ?'K ^{0 }decays
antelj, L.; Yusa, Y.; Abdesselam, A.; Adachi, I.; Aihara, H.; Al Said, S.; Asner, D. M.; Aulchenko, V.; Aushev, T.; Ayad, R.; Bahinipati, S.
2014-10-29
We present a measurement of the time-dependent CP violation parameters in B^{0} ? ?'K^{0} decays. The measurement is based on the full data sample containing 77210^{6} BB-bar pairs collected at the ?(4S) resonance using the Belle detector at the KEKB asymmetric-energy e^{+}e^{-} collider. The measured values of the mixing-induced and direct CP violation parameters are: sin 2?_{1}^{eff} = +0.68 0.07 0.03, A_{?'K0} = +0.03 0.05 0.04, where the first uncertainty is statistical and the second is systematic. The values obtained are the most accurate to date. Furthermore, these results are consistent with our previous measurements and with the world-average value of sin 2?_{1} measured in B^{0} ? J/?K^{0} decays.
Ho, Cong Son Tan, Seng Ghee; Jalil, Mansoor B. A.
2014-05-14
The generation of spin current and spin polarization in a two-dimensional electron gas structure is studied in the presence of Dresselhaus and Rashba spin-orbit couplings (SOC), the strength of the latter being modulated in time by an ac gate voltage. By means of the non-Abelian gauge field approach, we established the relation between the Lorentz spin force and the spin current in the SOC system, and showed that the longitudinal component of the spin force induces a transverse spin current. For a constant (time-invariant) Rashba system, we recover the universal spin Hall conductivity of e/(8π) , derived previously via the Berry phase and semi-classical methods. In the case of a time-dependent SOC system, the spin current is sustained even under strong impurity scattering. We evaluated the ac spin current generated by a time-modulated Rashba SOC in the absence of any dc electric field. The magnitude of the spin current reaches a maximum when the modulation frequency matches the Larmor frequency of the electrons.
Steady-state benchmarks of DK4D: A time-dependent, axisymmetric drift-kinetic equation solver
Lyons, B. C.; Jardin, S. C.; Ramos, J. J.
2015-05-15
The DK4D code has been written to solve a set of time-dependent, axisymmetric, finite-Larmor-radius drift-kinetic equations (DKEs) for the non-Maxwellian part of the electron and ion distribution functions using the full, linearized Fokker–Planck–Landau collision operator. The plasma is assumed to be in the low- to finite-collisionality regime, as is found in the cores of modern and future magnetic confinement fusion experiments. Each DKE is formulated such that the perturbed distribution function carries no net density, parallel momentum, or kinetic energy. Rather, these quantities are contained within the background Maxwellians and would be evolved by an appropriate set of extended magnetohydrodynamic (MHD) equations. This formulation allows for straight-forward coupling of DK4D to existing extended MHD time evolution codes. DK4D uses a mix of implicit and explicit temporal representations and finite element and spectral spatial representations. These, along with other computational methods used, are discussed extensively. Steady-state benchmarks are then presented comparing the results of DK4D to expected analytic results at low collisionality, qualitatively, and to the Sauter analytic fits for the neoclassical conductivity and bootstrap current, quantitatively. These benchmarks confirm that DK4D is capable of solving for the correct, gyroaveraged distribution function in stationary magnetic equilibria. Furthermore, the results presented demonstrate how the exact drift-kinetic solution varies with collisionality as a function of the magnetic moment and the poloidal angle.
Yu, Haitao; Guo, Xinmeng; Wang, Jiang Deng, Bin; Wei, Xile
2014-09-01
The phenomenon of stochastic resonance in Newman-Watts small-world neuronal networks is investigated when the strength of synaptic connections between neurons is adaptively adjusted by spike-time-dependent plasticity (STDP). It is shown that irrespective of the synaptic connectivity is fixed or adaptive, the phenomenon of stochastic resonance occurs. The efficiency of network stochastic resonance can be largely enhanced by STDP in the coupling process. Particularly, the resonance for adaptive coupling can reach a much larger value than that for fixed one when the noise intensity is small or intermediate. STDP with dominant depression and small temporal window ratio is more efficient for the transmission of weak external signal in small-world neuronal networks. In addition, we demonstrate that the effect of stochastic resonance can be further improved via fine-tuning of the average coupling strength of the adaptive network. Furthermore, the small-world topology can significantly affect stochastic resonance of excitable neuronal networks. It is found that there exists an optimal probability of adding links by which the noise-induced transmission of weak periodic signal peaks.
Gomez, Antonio; DeRego, Paul Jeffrey; Ferrel, Patrick Andrew; Thom, Robert Anthony; Trujillo, Joshua J.; Herridge, Brian
2014-08-19
An apparatus for generating seismic waves includes a housing, a strike surface within the housing, and a hammer movably disposed within the housing. An actuator induces a striking motion in the hammer such that the hammer impacts the strike surface as part of the striking motion. The actuator is selectively adjustable to change characteristics of the striking motion and characteristics of seismic waves generated by the impact. The hammer may be modified to change the physical characteristics of the hammer, thereby changing characteristics of seismic waves generated by the hammer. The hammer may be disposed within a removable shock cavity, and the apparatus may include two hammers and two shock cavities positioned symmetrically about a center of the apparatus.
Gomez, Antonio; DeRego, Paul Jeffrey; Ferrell, Patrick Andrew; Thom, Robert Anthony; Trujillo, Joshua J.; Herridge, Brian
2015-09-29
An apparatus for generating seismic waves includes a housing, a strike surface within the housing, and a hammer movably disposed within the housing. An actuator induces a striking motion in the hammer such that the hammer impacts the strike surface as part of the striking motion. The actuator is selectively adjustable to change characteristics of the striking motion and characteristics of seismic waves generated by the impact. The hammer may be modified to change the physical characteristics of the hammer, thereby changing characteristics of seismic waves generated by the hammer. The hammer may be disposed within a removable shock cavity, and the apparatus may include two hammers and two shock cavities positioned symmetrically about a center of the apparatus.
Turbocharging Quantum Tomography.
Blume-Kohout, Robin J; Gamble, John King,; Nielsen, Erik; Maunz, Peter Lukas Wilhelm; Scholten, Travis L.; Rudinger, Kenneth Michael
2015-01-01
Quantum tomography is used to characterize quantum operations implemented in quantum information processing (QIP) hardware. Traditionally, state tomography has been used to characterize the quantum state prepared in an initialization procedure, while quantum process tomography is used to characterize dynamical operations on a QIP system. As such, tomography is critical to the development of QIP hardware (since it is necessary both for debugging and validating as-built devices, and its results are used to influence the next generation of devices). But tomography su %7C ers from several critical drawbacks. In this report, we present new research that resolves several of these flaws. We describe a new form of tomography called gate set tomography (GST), which unifies state and process tomography, avoids prior methods critical reliance on precalibrated operations that are not generally available, and can achieve unprecedented accuracies. We report on theory and experimental development of adaptive tomography protocols that achieve far higher fidelity in state reconstruction than non-adaptive methods. Finally, we present a new theoretical and experimental analysis of process tomography on multispin systems, and demonstrate how to more e %7C ectively detect and characterize quantum noise using carefully tailored ensembles of input states.
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Syracuse, Ellen M.; Maceira, Monica; Prieto, German A.; Zhang, Haijiang; Ammon, Charles J.
2016-04-12
Subduction beneath the northernmost Andes in Colombia is complex. Based on seismicity distributions, multiple segments of slab appear to be subducting, and arc volcanism ceases north of 5° N. Here, we illuminate the subduction system through hypocentral relocations and Vp and Vs models resulting from the joint inversion of local body wave arrivals, surface wave dispersion measurements, and gravity data. The simultaneous use of multiple data types takes advantage of the differing sensitivities of each data type, resulting in velocity models that have improved resolution at both shallower and deeper depths than would result from traditional travel time tomography alone.more » The relocated earthquake dataset and velocity model clearly indicate a tear in the Nazca slab at 5° N, corresponding to a 250-km shift in slab seismicity and the termination of arc volcanism. North of this tear, the slab is flat, and it comprises slabs of two sources: the Nazca and Caribbean plates. The Bucaramanga nest, a small region of among the most intense intermediate-depth seismicity globally, is associated with the boundary between these two plates and possibly with a zone of melting or elevated water content, based on reduced Vp and increased Vp/Vs. As a result, we also use relocated seismicity to identify two new faults in the South American plate, one related to plate convergence and one highlighted by induced seismicity.« less
Meng, Qingyong, E-mail: mengqingyong@dicp.ac.cn [State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Zhongshan Road 457, 116023 Dalian (China); Meyer, Hans-Dieter, E-mail: hans-dieter.meyer@pci.uni-heidelberg.de [Theoretische Chemie, Physikalisch-Chemisches Institut, Ruprecht-Karls Universitt Heidelberg, Im Neuenheimer Feld 229, D-69120 Heidelberg (Germany)
2014-09-28
Employing the multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) method in conjunction with the multistate multimode vibronic coupling Hamiltonian (MMVCH) model, we perform a full dimensional (9D) quantum dynamical study on the simplest Criegee intermediate, formaldehyde oxide, in five lower-lying singlet electronic states. The ultraviolet (UV) spectrum is then simulated by a Fourier transform of the auto-correlation function. The MMVCH model is built based on extensive MRCI(8e,8o)/aug-cc-pVTZ calculations. To ensure a fast convergence of the final calculations, a large number of ML-MCTDH test calculations is performed to find an appropriate multilayer separations (ML-trees) of the ML-MCTDH nuclear wave functions, and the dynamical calculations are carefully checked to ensure that the calculations are well converged. To compare the computational efficiency, standard MCTDH simulations using the same Hamiltonian are also performed. A comparison of the MCTDH and ML-MCTDH calculations shows that even for the present not-too-large system (9D here) the ML-MCTDH calculations can save a considerable amount of computational resources while producing identical spectra as the MCTDH calculations. Furthermore, the present theoretical B{sup ~} {sup 1}A{sup ?}?X{sup ~} {sup 1}A{sup ?} UV spectral band and the corresponding experimental measurements [J. M. Beames, F. Liu, L. Lu, and M. I. Lester, J. Am. Chem. Soc. 134, 2004520048 (2012); L. Sheps, J. Phys. Chem. Lett. 4, 42014205 (2013); W.-L. Ting, Y.-H. Chen, W. Chao, M. C. Smith, and J. J.-M. Lin, Phys. Chem. Chem. Phys. 16, 1043810443 (2014)] are discussed. To the best of our knowledge, this is the first theoretical UV spectrum simulated for this molecule including nuclear motion beyond an adiabatic harmonic approximation.
Morzan, Uriel N.; Ramrez, Francisco F.; Scherlis, Damin A. E-mail: mcgl@qb.ffyb.uba.ar; Lebrero, Mariano C. Gonzlez E-mail: mcgl@qb.ffyb.uba.ar
2014-04-28
This article presents a time dependent density functional theory (TDDFT) implementation to propagate the Kohn-Sham equations in real time, including the effects of a molecular environment through a Quantum-Mechanics Molecular-Mechanics (QM-MM) hamiltonian. The code delivers an all-electron description employing Gaussian basis functions, and incorporates the Amber force-field in the QM-MM treatment. The most expensive parts of the computation, comprising the commutators between the hamiltonian and the density matrixrequired to propagate the electron dynamics, and the evaluation of the exchange-correlation energy, were migrated to the CUDA platform to run on graphics processing units, which remarkably accelerates the performance of the code. The method was validated by reproducing linear-response TDDFT results for the absorption spectra of several molecular species. Two different schemes were tested to propagate the quantum dynamics: (i) a leap-frog Verlet algorithm, and (ii) the Magnus expansion to first-order. These two approaches were confronted, to find that the Magnus scheme is more efficient by a factor of six in small molecules. Interestingly, the presence of iron was found to seriously limitate the length of the integration time step, due to the high frequencies associated with the core-electrons. This highlights the importance of pseudopotentials to alleviate the cost of the propagation of the inner states when heavy nuclei are present. Finally, the methodology was applied to investigate the shifts induced by the chemical environment on the most intense UV absorption bands of two model systems of general relevance: the formamide molecule in water solution, and the carboxy-heme group in Flavohemoglobin. In both cases, shifts of several nanometers are observed, consistently with the available experimental data.
Miniati, Francesco
2015-02-10
We use the Matryoshka run to study the time-dependent statistics of structure-formation-driven turbulence in the intracluster medium of a 10{sup 15} M {sub ?} galaxy cluster. We investigate the turbulent cascade in the inner megaparsec for both compressional and incompressible velocity components. The flow maintains approximate conditions of fully developed turbulence, with departures thereof settling in about an eddy-turnover time. Turbulent velocity dispersion remains above 700kms{sup 1} even at low mass accretion rate, with the fraction of compressional energy between 10% and 40%. The normalization and the slope of the compressional turbulence are susceptible to large variations on short timescales, unlike the incompressible counterpart. A major merger occurs around redshift z ? 0 and is accompanied by a long period of enhanced turbulence, ascribed to temporal clustering of mass accretion related to spatial clustering of matter. We test models of stochastic acceleration by compressional modes for the origin of diffuse radio emission in galaxy clusters. The turbulence simulation model constrains an important unknown of this complex problem and brings forth its dependence on the elusive microphysics of the intracluster plasma. In particular, the specifics of the plasma collisionality and the dissipation physics of weak shocks affect the cascade of compressional modes with strong impact on the acceleration rates. In this context radio halos emerge as complex phenomena in which a hierarchy of processes acting on progressively smaller scales are at work. Stochastic acceleration by compressional modes implies statistical correlation of radio power and spectral index with merging cores distance, both testable in principle with radio surveys.
Nonstructural seismic restraint guidelines
Butler, D.M.; Czapinski, R.H.; Firneno, M.J.; Feemster, H.C.; Fornaciari, N.R.; Hillaire, R.G.; Kinzel, R.L.; Kirk, D.; McMahon, T.T.
1993-08-01
The Nonstructural Seismic Restraint Guidelines provide general information about how to secure or restrain items (such as material, equipment, furniture, and tools) in order to prevent injury and property, environmental, or programmatic damage during or following an earthquake. All SNL sites may experience earthquakes of magnitude 6.0 or higher on the Richter scale. Therefore, these guidelines are written for all SNL sites.
Subsurface void detection using seismic tomographic imaging
Gritto, Roland
2003-06-26
Tomographic imaging has been widely used in scientific and medical fields to remotely image media in a nondestructive way. This paper introduces a spectrum of seismic imaging applications to detect and characterize voids in coal mines. The application of seismic waves to detect changes in coal relies on two types of waves: body waves refracted along the interface between coal and bedrock (i.e., refracted P-waves) and channel waves that propagate directly through the coal (dispersive wave trains of the Rayleigh or Love type). For example, a P-wave tomography study to find underlying old mine workings in a coal mine in England, produced velocity patterns that revealed increases in velocity where high stress concentrations occur in the rock, which are most likely connected to old pillars left in support of the old working areas. At the same time, low velocities were found in areas of low stress concentrations, which are related to roof collapses indicating the locations of mined areas below. The application of channel wave tomography to directly image the presence of gaseous CO{sub 2} in a low velocity oil reservoir showed that the injected CO{sub 2} followed an ancient flow channel in the reservoir migrating from the injector to the producer well. The study showed how channel waves are preferable over refracted P-waves, as the latter were only marginally affected by the presence of the gas in the low-velocity channel. Similar approaches show great promise for the detection of voids in coal mines. Finally, a newly developed technique, based on scattering theory, revealed that the location and the size of a subsurface cavity could be accurately determined even in the presence of strong correlated and uncorrelated noise.
Alonso, Juan J.; Iaccarino, Gianluca
2013-08-25
solution to the long-time integration problem of spectral chaos approaches; 4. A rigorous methodology to account for aleatory and epistemic uncertainties, to emphasize the most important variables via dimension reduction and dimension-adaptive refinement, and to support fusion with experimental data using Bayesian inference; 5. The application of novel methodologies to time-dependent reliability studies in wind turbine applications including a number of efforts relating to the uncertainty quantification in vertical-axis wind turbine applications. In this report, we summarize all accomplishments in the project (during the time period specified) focusing on advances in UQ algorithms and deployment efforts to the wind turbine application area. Detailed publications in each of these areas have also been completed and are available from the respective conference proceedings and journals as detailed in a later section.
Amplitude Analysis and Measurement of the Time-dependent CP Asymmetry of B0 to KsKsKs Decays
Lees, J.P.
2012-04-11
We present the first results on the Dalitz-plot structure and improved measurements of the time-dependent CP-violation parameters of the process B{sup 0} {yields} K{sub S}{sup 0}K{sub S}{sup 0}K{sub S}{sup 0} obtained using 468 x 10{sup 6} B{bar B} decays collected with the BABAR detector at the PEP-II asymmetric-energy B factory at SLAC. The Dalitz-plot structure is probed by a time-integrated amplitude analysis that does not distinguish between B{sup 0} and {bar B}{sup 0} decays. We measure the total inclusive branching fraction {Beta}(B{sup 0} {yields} K{sub S}{sup 0}K{sub S}{sup 0}K{sub S}{sup 0}) = (6.19 {+-} 0.48 {+-} 0.15 {+-} 0.12) x 10{sup -6}, where the first uncertainty is statistical, the second is systematic, and the third represents the Dalitz-plot signal model dependence. We also observe evidence for the intermediate resonant states f{sub 0}(980), f{sub 0}(1710), and f{sub 2}(2010). Their respective product branching fractions are measured to be (2.70{sub -1.19}{sup +1.25} {+-} 0.36 {+-} 1.17) x 10{sup -6}, (0.50{sub -0.24}{sup +0.46} {+-} 0.04 {+-} 0.10) x 10{sup -6}, and (0.54{sub -0.20}{sup +0.21} {+-} 0.03 {+-} 0.52) x 10{sup -6}. Additionally, we determine the mixing-induced CP-violation parameters to be S = -0.94{sub -0.21}{sup +0.24} {+-} 0.06 and C = -0.17 {+-} 0.18 {+-} 0.04, where the first uncertainty is statistical and the second is systematic. These values are in agreement with the standard model expectation.
Hanford Sitewide Probabilistic Seismic Hazard Analysis
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
B - PPRP Closure Letter Hanford Sitewide Probabilistic Seismic Hazard Analysis 2014 B.1 Appendix B PPRP Closure Letter 2014 Hanford Sitewide Probabilistic Seismic Hazard Analysis B.2 Hanford Sitewide Probabilistic Seismic Hazard Analysis 2014 B.3 2014 Hanford Sitewide Probabilistic Seismic Hazard Analysis B.4 Hanford Sitewide Probabilistic Seismic Hazard Analysis 2014 B.5
High Order Seismic Simulations
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Order Seismic Simulations on the Intel Xeon Phi Processor (Knights Landing) Alexander Heinecke 1 , Alexander Breuer 2 , Michael Bader 3 , and Pradeep Dubey 1 1 Intel Corporation, 2200 Mission College Blvd., Santa Clara 95054, CA, USA 2 University of California, San Diego, 9500 Gilman Dr., La Jolla 92093, CA, USA 3 Technische Universit¨ at M¨ unchen, Boltzmannstr. 3, D-85748 Garching, Germany Abstract. We present a holistic optimization of the ADER-DG finite element software SeisSol targeting
Seismic Structure And Seismicity Of The Cooling Lava Lake Of...
Of The Cooling Lava Lake Of Kilauea Iki, Hawaii Abstract The use of multiple methods is indispensable for the determination of the seismic properties of a complex body...
Induced Seismicity | Open Energy Information
of Project DOE Funding Total Project Cost Mapping Diffuse Seismicity for Geothermal Reservoir Management with Matched Field Processing California Lawrence Livermore National...
Hanford Sitewide Probabilistic Seismic Hazard Analysis
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
F - Seismicity Relocation Analyses Hanford Sitewide Probabilistic Seismic Hazard Analysis 2014 F.1 Appendix F Seismicity Relocation Analyses Final Report: High-Resolution Seismicity Study of the Yakima Fold and Thrust Belt Region, Washington Prepared by Clifford H. Thurber Department of Geoscience University of Wisconsin-Madison 1215 W. Dayton St. Madison, WI 53706 January 31, 2014 Final Report: Hanford Site-Wide Probabilistic Seismic Hazard Analysis (PSHA): High-Resolution Seismicity Analysis
Automating Shallow Seismic Imaging
Steeples, Don W.
2004-12-09
This seven-year, shallow-seismic reflection research project had the aim of improving geophysical imaging of possible contaminant flow paths. Thousands of chemically contaminated sites exist in the United States, including at least 3,700 at Department of Energy (DOE) facilities. Imaging technologies such as shallow seismic reflection (SSR) and ground-penetrating radar (GPR) sometimes are capable of identifying geologic conditions that might indicate preferential contaminant-flow paths. Historically, SSR has been used very little at depths shallower than 30 m, and even more rarely at depths of 10 m or less. Conversely, GPR is rarely useful at depths greater than 10 m, especially in areas where clay or other electrically conductive materials are present near the surface. Efforts to image the cone of depression around a pumping well using seismic methods were only partially successful (for complete references of all research results, see the full Final Technical Report, DOE/ER/14826-F), but peripheral results included development of SSR methods for depths shallower than one meter, a depth range that had not been achieved before. Imaging at such shallow depths, however, requires geophone intervals of the order of 10 cm or less, which makes such surveys very expensive in terms of human time and effort. We also showed that SSR and GPR could be used in a complementary fashion to image the same volume of earth at very shallow depths. The primary research focus of the second three-year period of funding was to develop and demonstrate an automated method of conducting two-dimensional (2D) shallow-seismic surveys with the goal of saving time, effort, and money. Tests involving the second generation of the hydraulic geophone-planting device dubbed the ''Autojuggie'' showed that large numbers of geophones can be placed quickly and automatically and can acquire high-quality data, although not under rough topographic conditions. In some easy-access environments, this device could
Seismic & Natural Phenomena Hazards | Department of Energy
Broader source: Energy.gov (indexed) [DOE]
designed to withstand the hazards. CNS maintains a panel of experts known as the Seismic Lessons-Learned Panel, which meets periodically to discuss seismic issues impacting DOE...
Induced Seismicity Impact | Open Energy Information
Seismicity Impact Jump to: navigation, search Retrieved from "http:en.openei.orgwindex.php?titleInducedSeismicityImpact&oldid612409" Feedback Contact needs updating...
Elastic-Wavefield Seismic Stratigraphy: A New Seismic Imaging Technology
Bob A. Hardage; Milo M. Backus; Michael V. DeAngelo; Sergey Fomel; Khaled Fouad; Robert J. Graebner; Paul E. Murray; Randy Remington; Diana Sava
2006-07-31
The purpose of our research has been to develop and demonstrate a seismic technology that will provide the oil and gas industry a better methodology for understanding reservoir and seal architectures and for improving interpretations of hydrocarbon systems. Our research goal was to expand the valuable science of seismic stratigraphy beyond the constraints of compressional (P-P) seismic data by using all modes (P-P, P-SV, SH-SH, SV-SV, SV-P) of a seismic elastic wavefield to define depositional sequences and facies. Our objective was to demonstrate that one or more modes of an elastic wavefield may image stratal surfaces across some stratigraphic intervals that are not seen by companion wave modes and thus provide different, but equally valid, information regarding depositional sequences and sedimentary facies within that interval. We use the term elastic wavefield stratigraphy to describe the methodology we use to integrate seismic sequences and seismic facies from all modes of an elastic wavefield into a seismic interpretation. We interpreted both onshore and marine multicomponent seismic surveys to select the data examples that we use to document the principles of elastic wavefield stratigraphy. We have also used examples from published papers that illustrate some concepts better than did the multicomponent seismic data that were available for our analysis. In each interpretation study, we used rock physics modeling to explain how and why certain geological conditions caused differences in P and S reflectivities that resulted in P-wave seismic sequences and facies being different from depth-equivalent S-wave sequences and facies across the targets we studied.
Seismic Data Gathering and Validation
Coleman, Justin
2015-02-01
Three recent earthquakes in the last seven years have exceeded their design basis earthquake values (so it is implied that damage to SSC’s should have occurred). These seismic events were recorded at North Anna (August 2011, detailed information provided in [Virginia Electric and Power Company Memo]), Fukushima Daichii and Daini (March 2011 [TEPCO 1]), and Kaswazaki-Kariwa (2007, [TEPCO 2]). However, seismic walk downs at some of these plants indicate that very little damage occurred to safety class systems and components due to the seismic motion. This report presents seismic data gathered for two of the three events mentioned above and recommends a path for using that data for two purposes. One purpose is to determine what margins exist in current industry standard seismic soil-structure interaction (SSI) tools. The second purpose is the use the data to validated seismic site response tools and SSI tools. The gathered data represents free field soil and in-structure acceleration time histories data. Gathered data also includes elastic and dynamic soil properties and structural drawings. Gathering data and comparing with existing models has potential to identify areas of uncertainty that should be removed from current seismic analysis and SPRA approaches. Removing uncertainty (to the extent possible) from SPRA’s will allow NPP owners to make decisions on where to reduce risk. Once a realistic understanding of seismic response is established for a nuclear power plant (NPP) then decisions on needed protective measures, such as SI, can be made.
Teleseismic-Seismic Monitoring At Central Nevada Seismic Zone...
DOE-funding Unknown References Glenn Biasi, Leiph Preston, Ileana Tibuleac (2009) Body Wave Tomography For Regional Scale Assessment Of Geothermal Indicators In The Western...
Seismic event classification system
Dowla, Farid U.; Jarpe, Stephen P.; Maurer, William
1994-01-01
In the computer interpretation of seismic data, the critical first step is to identify the general class of an unknown event. For example, the classification might be: teleseismic, regional, local, vehicular, or noise. Self-organizing neural networks (SONNs) can be used for classifying such events. Both Kohonen and Adaptive Resonance Theory (ART) SONNs are useful for this purpose. Given the detection of a seismic event and the corresponding signal, computation is made of: the time-frequency distribution, its binary representation, and finally a shift-invariant representation, which is the magnitude of the two-dimensional Fourier transform (2-D FFT) of the binary time-frequency distribution. This pre-processed input is fed into the SONNs. These neural networks are able to group events that look similar. The ART SONN has an advantage in classifying the event because the types of cluster groups do not need to be pre-defined. The results from the SONNs together with an expert seismologist's classification are then used to derive event classification probabilities.
Seismic event classification system
Dowla, F.U.; Jarpe, S.P.; Maurer, W.
1994-12-13
In the computer interpretation of seismic data, the critical first step is to identify the general class of an unknown event. For example, the classification might be: teleseismic, regional, local, vehicular, or noise. Self-organizing neural networks (SONNs) can be used for classifying such events. Both Kohonen and Adaptive Resonance Theory (ART) SONNs are useful for this purpose. Given the detection of a seismic event and the corresponding signal, computation is made of: the time-frequency distribution, its binary representation, and finally a shift-invariant representation, which is the magnitude of the two-dimensional Fourier transform (2-D FFT) of the binary time-frequency distribution. This pre-processed input is fed into the SONNs. These neural networks are able to group events that look similar. The ART SONN has an advantage in classifying the event because the types of cluster groups do not need to be pre-defined. The results from the SONNs together with an expert seismologist's classification are then used to derive event classification probabilities. 21 figures.
Measurement of Time-Dependent CP Asymmetries and the CP-Odd Fraction in the Decay B0->D*+D*-
Aubert, B.; Barate, R.; Boutigny, D.; Couderc, F.; Karyotakis, Y.; Lees, J.P.; Poireau, V.; Tisserand, V.; Zghiche, A.; Grauges, E.; Palano, A.; Pappagallo, M.; Pompili, A.; Chen, J.C.; Qi, N.D.; Rong, G.; Wang, P.; Zhu, Y.S.; Eigen, G.; Ofte, I.; Stugu, B. /Bergen U. /LBL, Berkeley /UC, Berkeley /Birmingham U. /Ruhr U., Bochum /Bristol U. /British Columbia U. /Brunel U. /Novosibirsk, IYF /UC, Irvine /UCLA /UC, Riverside /UC, San Diego /UC, Santa Barbara /UC, Santa Cruz /Caltech /Cincinnati U. /Colorado U. /Colorado State U. /Dortmund U. /Dresden, Tech. U. /Ecole Polytechnique /Edinburgh U. /Ferrara U. /INFN, Ferrara /Frascati /Genoa U. /INFN, Genoa /Harvard U. /Heidelberg U. /Imperial Coll., London /Iowa U. /Iowa State U. /Orsay, LAL /LLNL, Livermore /Liverpool U. /Queen Mary, U. of London /Royal Holloway, U. of London /Louisville U. /Manchester U. /Maryland U. /Massachusetts U., Amherst /MIT, LNS /McGill U. /Milan U. /INFN, Milan /Mississippi U. /Montreal U. /Mt. Holyoke Coll. /Naples U. /INFN, Naples /NIKHEF, Amsterdam /Notre Dame U. /Ohio State U. /Oregon U. /Padua U. /INFN, Padua /Paris U., VI-VII /Pennsylvania U. /Perugia U. /INFN, Perugia /Pisa U. /INFN, Pisa /Prairie View A-M /Princeton U. /Rome U. /INFN, Rome /Rostock U. /Rutherford /DAPNIA, Saclay /South Carolina U. /SLAC /Stanford U., Phys. Dept. /SUNY, Stony Brook /Tennessee U. /Texas U. /Texas U., Dallas /Turin U. /INFN, Turin /Trieste U. /INFN, Trieste /Valencia U., IFIC /Vanderbilt U. /Victoria U. /Warwick U. /Wisconsin U., Madison /Yale U.
2005-07-06
We present an updated measurement of time-dependent CP asymmetries and the CP-odd fraction in the decay B{sup 0} D*{sup +}D*{sup -} using 232 x 10{sup 6} B{bar B} pairs collected by the BABAR detector at the PEP-II B factory. We determine the CP-odd fraction to be 0.125 {+-} 0.044(stat) {+-} 0.007(syst). The time-dependent CP asymmetry parameters C{sub +} and S{sub +} are determined to be 0.06 {+-} 0.17(stat) {+-} 0.03(syst) and -0.75 {+-} 0.25(stat) {+-} 0.03(syst), respectively. The Standard Model predicts these parameters to be 0 and -sin2{beta}, respectively, in the absence of penguin amplitude contributions.
Prayitno, T. B.
2014-03-24
We have imposed the conditions in order to preserve the real-valued partition function in the case of onedimensional Gross-Pitaevskii equation coupled by time-dependent potential. In this case we have solved the Gross-Pitaevskii equation by means of the time-dependent perturbation theory by extending the previous work of Kivshar et al. [Phys. Lett A 278, 225–230 (2001)]. To use the method, we have treated the equation as the macroscopic quantum oscillator and found that the expression of the partition function explicitly has complex values. In fact, we have to choose not only the appropriate functions but also the suitable several values of the potential to keep the real-valued partition function.
Elastic-Wavefield Seismic Stratigraphy: A New Seismic Imaging Technology
Bob A. Hardage
2005-07-31
We have developed a numerical technique that will adjust 3-D S-wave seismic images so that they are depth equivalent to 3-D P-wave seismic images. The ability to make this type of P-SV to P-P depth registration is critical to our elastic wavefield seismic stratigraphy research because we now have higher confidence that depth-equivalent data windows are being used in the P-SV to P-P comparisons that we are making.
Quantifying the {sup 12}C+{sup 12}C sub-Coulomb fusion with the time-dependent wave-packet method
Diaz-Torres, Alexis; Wiescher, Michael
2012-10-20
This contribution provides a preliminary study of the {sup 12}C+{sup 12}C sub-Coulomb fusion reaction using the time-dependent wave-packet method within a nuclear molecular picture. The theoretical sub-Coulomb fusion resonances seem to correspond well with observations. The present method might be a more suitable tool for expanding the cross-section predictions towards lower energies than the commonly used potential-model approximation.
White, J.R.
1985-04-01
This report provides the background theory, user input, and sample problems required for the efficient application of the DEPTH-CHARGE system - a code black for both static and time-dependent perturbation theory and data sensitivity analyses. The DEPTH-CHARGE system is of modular construction and has been implemented within the VENTURE-BURNER computational system at Oak Ridge National Laboratory. The DEPTH module (coupled with VENTURE) solves for the three adjoint functions of Depletion Perturbation Theory and calculates the desired time-dependent derivatives of the response with respect to the nuclide concentrations and nuclear data utilized in the reference model. The CHARGE code is a collection of utility routines for general data manipulation and input preparation and considerably extends the usefulness of the system through the automatic generation of adjoint sources, estimated perturbed responses, and relative data sensitivity coefficients. Combined, the DEPTH-CHARGE system provides, for the first time, a complete generalized first-order perturbation/sensitivity theory capability for both static and time-dependent analyses of realistic multidimensional reactor models. This current documentation incorporates minor revisions to the original DEPTH-CHARGE documentation (ORNL/CSD-78) to reflect some new capabilities within the individual codes.
Downhole hydraulic seismic generator
Gregory, Danny L.; Hardee, Harry C.; Smallwood, David O.
1992-01-01
A downhole hydraulic seismic generator system for transmitting energy wave vibrations into earth strata surrounding a borehole. The system contains an elongated, unitary housing operably connected to a well head aboveground by support and electrical cabling, and contains clamping apparatus for selectively clamping the housing to the walls of the borehole. The system further comprises a hydraulic oscillator containing a double-actuating piston whose movement is controlled by an electro-servovalve regulating a high pressure hydraulic fluid flow into and out of upper and lower chambers surrounding the piston. The spent hydraulic fluid from the hydraulic oscillator is stored and pumped back into the system to provide high pressure fluid for conducting another run at the same, or a different location within the borehole.
SEISMIC ANALYSIS FOR PRECLOSURE SAFETY
E.N. Lindner
2004-12-03
The purpose of this seismic preclosure safety analysis is to identify the potential seismically-initiated event sequences associated with preclosure operations of the repository at Yucca Mountain and assign appropriate design bases to provide assurance of achieving the performance objectives specified in the Code of Federal Regulations (CFR) 10 CFR Part 63 for radiological consequences. This seismic preclosure safety analysis is performed in support of the License Application for the Yucca Mountain Project. In more detail, this analysis identifies the systems, structures, and components (SSCs) that are subject to seismic design bases. This analysis assigns one of two design basis ground motion (DBGM) levels, DBGM-1 or DBGM-2, to SSCs important to safety (ITS) that are credited in the prevention or mitigation of seismically-initiated event sequences. An application of seismic margins approach is also demonstrated for SSCs assigned to DBGM-2 by showing a high confidence of a low probability of failure at a higher ground acceleration value, termed a beyond-design basis ground motion (BDBGM) level. The objective of this analysis is to meet the performance requirements of 10 CFR 63.111(a) and 10 CFR 63.111(b) for offsite and worker doses. The results of this calculation are used as inputs to the following: (1) A classification analysis of SSCs ITS by identifying potential seismically-initiated failures (loss of safety function) that could lead to undesired consequences; (2) An assignment of either DBGM-1 or DBGM-2 to each SSC ITS credited in the prevention or mitigation of a seismically-initiated event sequence; and (3) A nuclear safety design basis report that will state the seismic design requirements that are credited in this analysis. The present analysis reflects the design information available as of October 2004 and is considered preliminary. The evolving design of the repository will be re-evaluated periodically to ensure that seismic hazards are properly
Seismicity and seismic stress in the Coso Range, Coso geothermal...
and seismic stress in the Coso Range, Coso geothermal field, and Indian Wells Valley region, Southeast-Central California Jump to: navigation, search OpenEI Reference LibraryAdd to...
Seismic Fragility Analysis of a Degraded Condensate Storage Tank
Nie, J.; Braverman, J.; Hofmayer, C.; Choun, Y-S.; Kim, M.K.; Choi, I-K.
2011-05-16
for five cases: (1) a baseline analysis where the design condition (undegraded) is assumed, (2) a scenario with degraded stainless steel tank shell, (3) a scenario with degraded anchor bolts, (4) a scenario with anchorage concrete cracking, and (5) a perfect correlation of the above three degradation scenarios. This paper will present the methodology for the time-dependent fragility calculation and discuss the insights drawn from this study. To achieve a better understanding of the effects of aging on the performance of structures and passive components (SPCs) in nuclear power plants (NPPs), the Korea Atomic Energy Research Institute (KAERI) and Brookhaven National Laboratory (BNL) are collaborating to develop seismic fragility analysis methods that consider age-related degradation of SPCs. The rate of age-related degradation of SPCs was not found to be significantly large, but increasing as the plants get older. The slow but increasing rate of degradation of SPCs can potentially affect the safety of the older plants and become an important factor in decision making in the current trend of extending the operating license period of the plants (e.g., in the U.S. from 40 years to 60 years, and even potentially to 80 years). In this paper, a condensate storage tank (CST) was analyzed to estimate its seismic fragility capacities under various postulated degradation scenarios. This paper will present the methodology for the time-dependent fragility calculation and discuss the insights drawn from this study.
Seismicity and Improved Velocity Structure in Kuwait
Gok, R M; Rodgers, A J; Al-Enezi, A
2006-01-26
The Kuwait National Seismic Network (KNSN) began operation in 1997 and consists of nine three-component stations (eight short-period and one broadband) and is operated by the Kuwait Institute for Scientific Research. Although the region is largely believed to be aseismic, considerable local seismicity is recorded by KNSN. Seismic events in Kuwait are clustered in two main groups, one in the south and another in the north. The KNSN station distribution is able to capture the southern cluster within the footprint of the network but the northern cluster is poorly covered. Events tend to occur at depths ranging from the free surface to about 20 km. Events in the northern cluster tend to be deeper than those in south, however this might be an artifact of the station coverage. We analyzed KNSN recordings of nearly 200 local events to improve understanding of seismic events and crustal structure in Kuwait, performing several analyses with increasing complexity. First, we obtained an optimized one-dimensional (1D) velocity model for the entire region using the reported KNSN arrival times and routine locations. The resulting model is consistent with a recently obtained model from the joint inversion of receiver functions and surface wave group velocities. Crustal structure is capped by the thick ({approx} 7 km) sedimentary rocks of the Arabian Platform underlain by normal velocities for stable continental crust. Our new model has a crustal thickness of 44 km, constrained by an independent study of receiver functions and surface wave group velocities by Pasyanos et al (2006). Locations and depths of events after relocation with the new model are broadly consistent with those reported by KISR, although a few events move more than a few kilometers. We then used a double-difference tomography technique (tomoDD) to jointly locate the events and estimate three-dimensional (3D) velocity structure. TomoDD is based on hypoDD relocation algorithm and it makes use of both absolute and
Seismic Techniques | Open Energy Information
impedance boundary7 References (Majer, n.d.) "3-D Seismic Methods For Geothermal Reservoir Exploration and Assessment- Summary" 2.0 2.1 2.2 (Dobrin and Savit, 1988)...
Non-Linear Seismic Soil Structure Interaction (SSI) Method for...
Office of Environmental Management (EM)
Non-Linear Seismic Soil Structure Interaction (SSI) Method for Developing Non-Linear Seismic SSI Analysis Techniques Non-Linear Seismic Soil Structure Interaction (SSI) Method for ...
Seismic Analysis of Facilities and Evaluation of Risk | Department...
Office of Environmental Management (EM)
Seismic Analysis of Facilities and Evaluation of Risk Seismic Analysis of Facilities and Evaluation of Risk Seismic Analysis of Facilities and Evaluation of Risk Michael Salmon,...
Central Nevada Seismic Zone Geothermal Region | Open Energy Informatio...
Central Nevada Seismic Zone Geothermal Region (Redirected from Central Nevada Seismic Zone) Jump to: navigation, search GEOTHERMAL ENERGYGeothermal Home Central Nevada Seismic Zone...
The Idaho National Laboratory (INL) Senior Seismic Hazard Analysis...
Office of Environmental Management (EM)
The Idaho National Laboratory (INL) Senior Seismic Hazard Analysis Committee (SSHAC) Level 1 Seismic Hazard Analysis Presentation from the May 2015 Seismic Lessons-Learned Panel ...
Category:Borehole Seismic Techniques | Open Energy Information
out of 2 total. S Single-Well And Cross-Well Seismic Imaging 1 pages V Vertical Seismic Profiling 1 pages Pages in category "Borehole Seismic...
Down hole periodic seismic generator
Hardee, Harry C.; Hills, Richard G.; Striker, Richard P.
1989-01-01
A down hole periodic seismic generator system for transmitting variable frequency, predominantly shear-wave vibration into earth strata surrounding a borehole. The system comprises a unitary housing operably connected to a well head by support and electrical cabling and contains clamping apparatus for selectively clamping the housing to the walls of the borehole. The system further comprises a variable speed pneumatic oscillator and a self-contained pneumatic reservoir for producing a frequency-swept seismic output over a discrete frequency range.
Advanced downhole periodic seismic generator
Hardee, Harry C.; Hills, Richard G.; Striker, Richard P.
1991-07-16
An advanced downhole periodic seismic generator system for transmitting variable frequency, predominantly shear-wave vibration into earth strata surrounding a borehole. The system comprises a unitary housing operably connected to a well head by support and electrical cabling and contains clamping apparatus for selectively clamping the housing to the walls of the borehole. The system further comprises a variable speed pneumatic oscillator and a self-contained pneumatic reservoir for producing a frequency-swept seismic output over a discrete frequency range.
Byrd, John; Corlett, John; Doolittle, Larry; Fawley, William; Lidia, Steven; Penn, Gregory; Ratti, Alex; Staples, John; Wilcox Russell; Wurtele, Jonathan; Zholents, Alexander
2005-10-01
This report constitutes the third deliverable of LBNLs contracted role in the FERMI {at} Elettra Technical Optimization study. It describes proposed R&D activities for the baseline design of the Technical Optimization Study, initial studies of the RF gun mode-coupling and potential effects on beam dynamics, steady-state studies of FEL-2 performance to 10 nm, preliminary studies of time-dependent FEL-1 performance using electron bunch distribution from the start-to-end studies, and a preliminary investigation of a configuration with FEL sinclined at a small angle from the line of the linac.
Tussupbayev, Samat; Govind, Niranjan; Lopata, Kenneth A.; Cramer, Christopher J.
2015-03-10
We assess the performance of real-time time-dependent density functional theory (RT-TDDFT) for the calculation of absorption spectra of 12 organic dye molecules relevant to photovoltaics and dye sensitized solar cells with 8 exchange-correlation functionals (3 traditional, 3 global hybrids, and 2 range-separated hybrids). We compare the calculations with traditional linear-response (LR) TDDFT. In addition, we demonstrate the efficacy of the RT-TDDFT approach to calculate wide absorption spectra of two large chromophores relevant to photovoltaics and molecular switches.
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Muñoz Burgos, J. M.; Barbui, T.; Schmitz, O.; Stutman, D.; Tritz, K.
2016-07-11
Helium line-ratios for electron temperature (Te) and density (ne) plasma diagnostic in the Scrape-Off-Layer (SOL) and Edge regions of tokamaks are widely used. Due to their intensities and proximity of wavelengths, the singlet 667.8 and 728.1 nm, and triplet 706.5 nm visible lines have been typically preferred. Time-dependency of the triplet line (706.5 nm) has been previously analyzed in detail by including transient effects on line-ratios during gas-puff diagnostic applications. In this work, several line-ratio combinations within each of the two spin systems are analyzed with the purpose of eliminating transient effects to extend the application of this powerful diagnosticmore » to high temporal resolution characterization of plasmas. The analysis is done using synthetic emission modeling and diagnostic for low electron density NSTX SOL plasma conditions by several visible lines. Quasi-static equilibrium, and time-dependent models are employed to evaluate transient effects of the atomic population levels that may affect the derived electron temperatures and densities as the helium gas-puff penetrates the plasma. Ultimately, the analysis of a wider range of spectral lines will help to extend this powerful diagnostic to experiments where the wavelength range of the measured spectra may be constrained either by limitations of the spectrometer, or by other conflicting lines from different ions.« less
Zhu, Rui Dai, Jiao-Hua; Guo, Yong
2015-04-28
Interference between different quantum paths can generate Fano resonance. One of the examples is transport through a quasibound state driven by a time-dependent scattering potential. Previously it is found that Fano resonance occurs as a result of energy matching in one-dimensional systems. In this work, we demonstrate that when transverse motion is present, Fano resonance occurs precisely at the wavevector matching situation. Using the Floquet scattering theory, we considered the transport properties of a nonadiabatic time-dependent well both in a two-dimensional electron gas and monolayer graphene structure. Dispersion of the quasibound state of a static quantum well is obtained with transverse motion present. We found that Fano resonance occurs when the wavevector in the transport direction of one of the Floquet sidebands is exactly identical to that of the quasibound state in the well at equilibrium and follows the dispersion pattern of the latter. To observe the Fano resonance phenomenon in the transmission spectrum, we also considered the pumped shot noise properties when time and spatial symmetry secures vanishing current in the considered configuration. Prominent Fano resonance is found in the differential pumped shot noise with respect to the reservoir Fermi energy.
Chmura, Bartosz; Rode, Michal F.; Sobolewski, Andrzej L.; Lan Zhenggang
2009-10-07
The photoinduced electron-driven proton-transfer dynamics of the water-dimer system has been investigated by time-dependent quantum wave-packet calculations. The main nuclear degrees of freedom driving the system from the Frank-Condon region to the S{sub 0}-S{sub 1} conical intersection are the distance between the oxygen atoms and the displacement of the hydrogen atom from the oxygen-oxygen bond center. Two important coupling modes have been investigated: Rotation of the H-donating water dangling proton and asymmetric stretching of the H-accepting water dangling protons'O{sub a}H bonds. Potential energy surfaces of the ground and lowest excited electronic states have been constructed on the basis of ab initio calculations. The time-dependent quantum wave-packet propagation has been employed within the (2 + 1)-dimensional systems for the description of the nonadiabatic dynamics of water dimer. The effects of the initial vibrational state of the system on the electronic population transfer and dissociation dynamics are presented. To approximate the photochemical behavior of water dimer in bulk water, we add a boundary condition into the (2 + 1)-dimensional systems to simulate the existence of water bulk. The results provide insight into the mechanisms of excited state deactivation of the water-dimer system in gas phase and in bulk water through the electron-driven proton-transfer process.
DOE R&D Accomplishments [OSTI]
Hansche, B. D.
1983-01-01
Computed tomography (CT) is a relatively new radiographic technique which has become widely used in the medical field, where it is better known as computerized axial tomographic (CAT) scanning. This technique is also being adopted by the industrial radiographic community, although the greater range of densities, variation in samples sizes, plus possible requirement for finer resolution make it difficult to duplicate the excellent results that the medical scanners have achieved.
Grubsky, Victor; Romanoov, Volodymyr; Shoemaker, Keith; Patton, Edward Matthew; Jannson, Tomasz
2016-02-02
A Compton tomography system comprises an x-ray source configured to produce a planar x-ray beam. The beam irradiates a slice of an object to be imaged, producing Compton-scattered x-rays. The Compton-scattered x-rays are imaged by an x-ray camera. Translation of the object with respect to the source and camera or vice versa allows three-dimensional object imaging.
Seismic assessment of buried pipelines
Al-Chaar, G.; Brady, P.; Fernandez, G.
1995-12-31
A structure and its lifelines are closely linked because the disruption of lifeline systems will obstruct emergency service functions that are vitally needed after an earthquake. As an example of the criticality of these systems, the Association of Bay Area Government (ABAG) recorded thousands of leaks in pipelines that resulted in more than twenty million gallons of hazardous materials being released in several recorded earthquakes. The cost of cleaning the spills from these materials was very high. This information supports the development of seismic protection of lifeline systems. The US Army Corps of Engineers Construction Engineering Research Laboratories (USACERL) has, among its missions, the responsibility to develop seismic vulnerability assessment procedures for military installations. Within this mission, a preliminary research program to assess the seismic vulnerability of buried pipeline systems on military installations was initiated. Phase 1 of this research project resulted in two major studies. In the first, evaluating current procedures to seismically design or evaluate existing lifeline systems, the authors found several significant aspects that deserve special consideration and need to be addressed in future research. The second was focused on identifying parameters related to buried pipeline system vulnerability and developing a generalized analytical method to relate these parameters to the seismic vulnerability assessment of existing pipeline systems.
Generalized local emission tomography
Katsevich, Alexander J.
1998-01-01
Emission tomography enables locations and values of internal isotope density distributions to be determined from radiation emitted from the whole object. In the method for locating the values of discontinuities, the intensities of radiation emitted from either the whole object or a region of the object containing the discontinuities are inputted to a local tomography function .function..sub..LAMBDA..sup.(.PHI.) to define the location S of the isotope density discontinuity. The asymptotic behavior of .function..sub..LAMBDA..sup.(.PHI.) is determined in a neighborhood of S, and the value for the discontinuity is estimated from the asymptotic behavior of .function..sub..LAMBDA..sup.(.PHI.) knowing pointwise values of the attenuation coefficient within the object. In the method for determining the location of the discontinuity, the intensities of radiation emitted from an object are inputted to a local tomography function .function..sub..LAMBDA..sup.(.PHI.) to define the location S of the density discontinuity and the location .GAMMA. of the attenuation coefficient discontinuity. Pointwise values of the attenuation coefficient within the object need not be known in this case.
Katsevich, Alexander J.; Ramm, Alexander G.
1996-01-01
Local tomography is enhanced to determine the location and value of a discontinuity between a first internal density of an object and a second density of a region within the object. A beam of radiation is directed in a predetermined pattern through the region of the object containing the discontinuity. Relative attenuation data of the beam is determined within the predetermined pattern having a first data component that includes attenuation data through the region. In a first method for evaluating the value of the discontinuity, the relative attenuation data is inputted to a local tomography function .function..sub..LAMBDA. to define the location S of the density discontinuity. The asymptotic behavior of .function..sub..LAMBDA. is determined in a neighborhood of S, and the value for the discontinuity is estimated from the asymptotic behavior of .function..sub..LAMBDA.. In a second method for evaluating the value of the discontinuity, a gradient value for a mollified local tomography function .gradient..function..sub..LAMBDA..epsilon. (x.sub.ij) is determined along the discontinuity; and the value of the jump of the density across the discontinuity curve (or surface) S is estimated from the gradient values.
Advanced Seismic Data Analysis Program- The "Hot Pot" Project
Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site
Seismic Data Analysis Program - The "Hot Pot" Project Principal Investigator : Frank ... (2) * Innovative aspects Seismic data processing generally follows procedure ...
The Idaho National Laboratory (INL) Senior Seismic Hazard Analysis
Office of Environmental Management (EM)
Committee (SSHAC) Level 1 Seismic Hazard Analysis | Department of Energy The Idaho National Laboratory (INL) Senior Seismic Hazard Analysis Committee (SSHAC) Level 1 Seismic Hazard Analysis The Idaho National Laboratory (INL) Senior Seismic Hazard Analysis Committee (SSHAC) Level 1 Seismic Hazard Analysis Presentation from the May 2015 Seismic Lessons-Learned Panel Meeting. Seismic Hazard Definition: SSHAC Level 1 PSHA at MFC (416.86 KB) More Documents & Publications The INL Seismic Risk
Morris, C. L.; Bourke, M.; Byler, D. D.; Chen, C. F.; Hogan, G.; Hunter, J. F.; Kwiatkowski, K.; Mariam, F. G.; McClellan, K. J.; Merrill, F.; Morley, D. J.; Saunders, A.
2013-02-15
We present an assessment of x-rays and proton tomography as tools for studying the time dependence of the development of damage in fuel rods. We also show data taken with existing facilities at Los Alamos National Laboratory that support this assessment. Data on surrogate fuel rods have been taken using the 800 MeV proton radiography (pRad) facility at the Los Alamos Neutron Science Center (LANSCE), and with a 450 keV bremsstrahlung X-ray tomography facility. The proton radiography pRad facility at LANSCE can provide good position resolution (<70 {mu}m has been demonstrate, 20 {mu}m seems feasible with minor changes) for tomography on activated fuel rods. Bremsstrahlung x-rays may be able to provide better than 100 {mu}m resolution but further development of sources, collimation, and detectors is necessary for x-rays to deal with the background radiation for tomography of activated fuel rods.
Morris, Christopher L.; Bourke, Mark A.; Byler, Darrin D.; Chen, Ching-Fong; Hogan, Gary E.; Hunter, James F.; Kwiatkowski, Kris K.; Mariam, Fesseha G.; McClellan, Kenneth J.; Merrill, Frank E.; Morley, Deborah J.; Saunders, Alexander
2013-02-11
We present an assessment of x-rays and proton tomography as tools for studying the time dependence of the development of damage in fuel rods. Also, we show data taken with existing facilities at Los Alamos National Laboratory that support this assessment. Data on surrogate fuel rods has been taken using the 800 MeV proton radiography (pRad) facility at the Los Alamos Neutron Science Center (LANSCE), and with a 450 keV bremsstrahlung X-ray tomography facility. The proton radiography pRad facility at LANSCE can provide good position resolution (<70 μm has been demonstrate, 20 μm seems feasible with minor changes) for tomography on activated fuel rods. Bremsstrahlung x-rays may be able to provide better than 100 μm resolution but further development of sources, collimation and detectors is necessary for x-rays to deal with the background radiation for tomography of activated fuel rods.
Seismic Lessons-Learned Panel Meetings
Broader source: Energy.gov [DOE]
The Chief of Nuclear Security (CNS) maintains a panel of experts known as the Seismic Lessons-Learned Panel, which meets periodically to discuss seismic issues impacting DOE facilities.
National Nuclear Security Administration (NNSA)
Seismic Data for Evaluation of Ground Motion Hazards in Las Vegas in Support of Test Site ... Seismic Data for Evaluation of Ground Motion Hazards in Las Vegas in Support of Test Site ...
Development of a HT Seismic Tool
Broader source: Energy.gov [DOE]
The program objective is to design; fabricate and field test two high temperature (HT) seismic tools in an EGS application.
Active Seismic Techniques | Open Energy Information
Technique: Seismic Techniques Information Provided by Technique Lithology: Rock unit density influences elastic wave velocities. StratigraphicStructural: Structural geology-...
Vertical Seismic Profiling | Open Energy Information
Borehole Seismic Techniques Information Provided by Technique Lithology: Rock unit density influences elastic wave velocities. StratigraphicStructural: Structural geology-...
Teleseismic-Seismic Monitoring | Open Energy Information
Passive Seismic Techniques Information Provided by Technique Lithology: Rock unit density influences elastic wave velocities. StratigraphicStructural: Map geothermal...
Roy, M.; Mitchell, J. F.; Schiffer, P.
2000-02-17
The ground state of La{sub 1{minus}x}Ca{sub x}MnO{sub 3} changes from a ferromagnetic metallic to an antiferromagnetic charge-ordered state as a function of Ca concentration at x {approximately} 0.50. The authors present evidence from transport measurements on a sample with x = 0.50 that the two phases can coexist, in agreement with other observations of phase separation in these materials. They also observe that, by applying and then removing a magnetic field to the mainly charge-ordered state at some temperatures, they can ``magnetically anneal'' the charge order, resulting in a higher zero-field resistivity. They also observe logarithmic time dependence in both resistivity and magnetization after a field sweep at low temperatures.
Han, Lu; Liang, WanZhen; Zhao, Yi; Zhong, Xinxin
2014-06-07
The time-dependent wavepacket diffusive method [X. Zhong and Y. Zhao, J. Chem. Phys. 138, 014111 (2013)] is extended to investigate the energy relaxation and separation of a hot electron-hole pair in organic aggregates with incorporation of Coulomb interaction and electron-phonon coupling. The pair initial condition generated by laser pulse is represented by a Gaussian wavepacket with a central momentum. The results reveal that the hot electron energy relaxation is very well described by two rate processes with the fast rate much larger than the slow one, consistent with experimental observations, and an efficient electron-hole separation is accomplished accompanying the fast energy relaxation. Furthermore, although the extra energy indeed helps the separation by overcoming the Coulomb interaction, the width of initial wavepacket is much sensitive to the separation efficiency and the narrower wavepacket generates the more separated charges. This behavior may be useful to understand the experimental controversy of the hot carrier effect on charge separation.
Tao, Liang; Vanroose, Wim; Reps, Brian; Rescigno, Thomas N.; McCurdy, C. William
2009-09-08
We demonstrate that exterior complex scaling (ECS) can be used to impose outgoing wave boundary conditions exactly on solutions of the time-dependent Schrodinger equation for atoms in intense electromagnetic pulses using finite grid methods. The procedure is formally exact when applied in the appropriate gauge and is demonstrated in a calculation of high harmonic generation in which multiphoton resonances are seen for long pulse durations. However, we also demonstrate that while the application of ECS in this way is formally exact, numerical error can appear for long time propagations that can only be controlled by extending the finite grid. A mathematical analysis of the origins of that numerical error, illustrated with an analytically solvable model, is also given.
Tao Liang; Rescigno, T. N.; Vanroose, W.; Reps, B.; McCurdy, C. W.
2009-12-15
We demonstrate that exterior complex scaling (ECS) can be used to impose outgoing wave boundary conditions exactly on solutions of the time-dependent Schroedinger equation for atoms in intense electromagnetic pulses using finite grid methods. The procedure is formally exact when applied in the appropriate gauge and is demonstrated in a calculation of high-harmonic generation in which multiphoton resonances are seen for long pulse durations. However, we also demonstrate that while the application of ECS in this way is formally exact, numerical error can appear for long-time propagations that can only be controlled by extending the finite grid. A mathematical analysis of the origins of that numerical error, illustrated with an analytically solvable model, is also given.
Zhang, Yu; Mukamel, Shaul; Khalil, Munira; Govind, Niranjan
2015-11-09
Valence-to-core (VtC) X-ray emission spectroscopy (XES) has emerged as a power- ful technique for the structural characterization of complex organometallic compounds in realistic environments. Since the spectrum represents electronic transitions from the ligand molecular orbitals to the core holes of the metal centers, the approach is more chemically sensitive to the metal-ligand bonding character compared with con- ventional X-ray absorption techniques. In this paper we study how linear-response time-dependent density functional theory (LR-TDDFT) can be harnessed to simulate K-edge VtC X-ray emission spectra reliably. LR-TDDFT allows one to go beyond the single-particle picture that has been extensively used to simulate VtC-XES. We con- sider seven low- and high-spin model complexes involving chromium, manganese and iron transition metal centers. Our results are in good agreement with experiment.
Seismic Imaging Processing and Migration
Energy Science and Technology Software Center (OSTI)
2000-06-26
Salvo is a 3D, finite difference, prestack, depth migration code for parallel computers. It is also capable of processing 2D and poststack data. The code requires as input a seismic dataset, a velocity model and a file of parameters that allows the user to select various options. The code uses this information to produce a seismic image. Some of the options available to the user include the application of various filters and imaging conditions. Themore » code also incorporates phase encoding (patent applied for) to process multiple shots simultaneously.« less
Demján, Tamás; Vörös, Márton; Palummo, Maurizia; Gali, Adam
2014-08-14
Diamondoids are small diamond nanoparticles (NPs) that are built up from diamond cages. Unlike usual semiconductor NPs, their atomic structure is exactly known, thus they are ideal test-beds for benchmarking quantum chemical calculations. Their usage in spintronics and bioimaging applications requires a detailed knowledge of their electronic structure and optical properties. In this paper, we apply density functional theory (DFT) based methods to understand the electronic and optical properties of a few selected pure and modified diamondoids for which accurate experimental data exist. In particular, we use many-body perturbation theory methods, in the G{sub 0}W{sub 0} and G{sub 0}W{sub 0}+BSE approximations, and time-dependent DFT in the adiabatic local density approximation. We find large quasiparticle gap corrections that can exceed thrice the DFT gap. The electron-hole binding energy can be as large as 4 eV but it is considerably smaller than the GW corrections and thus G{sub 0}W{sub 0}+BSE optical gaps are about 50% larger than the Kohn-Sham (KS) DFT gaps. We find significant differences between KS time-dependent DFT and GW+BSE optical spectra on the selected diamondoids. The calculated G{sub 0}W{sub 0} quasiparticle levels agree well with the corresponding experimental vertical ionization energies. We show that nuclei dynamics in the ionization process can be significant and its contribution may reach about 0.5 eV in the adiabatic ionization energies.
Ma, HuiLi; Department of Chemical Physics, University of Science and Technology of China, Hefei 230026 ; Zhao, Yi; Liang, WanZhen
2014-03-07
A time-dependent approach is presented to simulate the two-photon absorption (TPA) and resonance hyper-Raman scattering (RHRS) spectra including Duschinsky rotation (mode-mixing) and Herzberg-Teller (HT) vibronic coupling effects. The computational obstacles for the excited-state geometries, vibrational frequencies, and nuclear derivatives of transition dipole moments, which enter the expressions of TPA and RHRS cross sections, are further overcome by the recently developed analytical excited-state energy derivative approaches in the framework of time-dependent density functional theory. The excited-state potential curvatures are evaluated at different levels of approximation to inspect the effects of frequency differences, mode-mixing and HT on TPA and RHRS spectra. Two types of molecules, one with high symmetry (formaldehyde, p-difluorobenzene, and benzotrifluoride) and the other with non-centrosymmetry (cis-hydroxybenzylidene-2,3-dimethylimidazolinone in the deprotonated anion state (HDBI{sup −})), are used as test systems. The calculated results reveal that it is crucial to adopt the exact excited-state potential curvatures in the calculations of TPA and RHRS spectra even for the high-symmetric molecules, and that the vertical gradient approximation leads to a large deviation. Furthermore, it is found that the HT contribution is evident in the TPA and RHRS spectra of HDBI{sup −} although its one- and two-photon transitions are strongly allowed, and its effect results in an obvious blueshift of the TPA maximum with respect to the one-photon absorption maximum. With the HT and solvent effects getting involved, the simulated blueshift of 1291 cm{sup −1} agrees well with the experimental measurement.
Seismic, shock, and vibration isolation - 1988
Chung, H. ); Mostaghel, N. )
1988-01-01
This book contains papers presented at a conference on pressure vessels and piping. Topics covered include: Design of R-FBI bearings for seismic isolation; Benefits of vertical and horizontal seismic isolation for LMR nuclear reactor units; and Some remarks on the use and perspectives of seismic isolation for fast reactors.
Radial reflection diffraction tomography
Lehman, Sean K.
2012-12-18
A wave-based tomographic imaging method and apparatus based upon one or more rotating radially outward oriented transmitting and receiving elements have been developed for non-destructive evaluation. At successive angular locations at a fixed radius, a predetermined transmitting element can launch a primary field and one or more predetermined receiving elements can collect the backscattered field in a "pitch/catch" operation. A Hilbert space inverse wave (HSIW) algorithm can construct images of the received scattered energy waves using operating modes chosen for a particular application. Applications include, improved intravascular imaging, bore hole tomography, and non-destructive evaluation (NDE) of parts having existing access holes.
Positron Emission Tomography (PET)
DOE R&D Accomplishments [OSTI]
Welch, M. J.
1990-01-01
Positron emission tomography (PET) assesses biochemical processes in the living subject, producing images of function rather than form. Using PET, physicians are able to obtain not the anatomical information provided by other medical imaging techniques, but pictures of physiological activity. In metaphoric terms, traditional imaging methods supply a map of the body's roadways, its, anatomy; PET shows the traffic along those paths, its biochemistry. This document discusses the principles of PET, the radiopharmaceuticals in PET, PET research, clinical applications of PET, the cost of PET, training of individuals for PET, the role of the United States Department of Energy in PET, and the futures of PET.
Functional performance requirements for seismic network upgrade
Lee, R.C.
1991-08-18
The SRL seismic network, established in 1976, was developed to monitor site and regional seismic activity that may have any potential to impact the safety or reduce containment capability of existing and planned structures and systems at the SRS, report seismic activity that may be relevant to emergency preparedness, including rapid assessments of earthquake location and magnitude, and estimates of potential on-site and off-site damage to facilities and lifelines for mitigation measures. All of these tasks require SRL seismologists to provide rapid analysis of large amounts of seismic data. The current seismic network upgrade, the subject of this Functional Performance Requirements Document, is necessary to improve system reliability and resolution. The upgrade provides equipment for the analysis of the network seismic data and replacement of old out-dated equipment. The digital network upgrade is configured for field station and laboratory digital processing systems. The upgrade consists of the purchase and installation of seismic sensors,, data telemetry digital upgrades, a dedicated Seismic Data Processing (SDP) system (already in procurement stage), and a Seismic Signal Analysis (SSA) system. The field stations and telephone telemetry upgrades include equipment necessary for three remote station upgrades including seismic amplifiers, voltage controlled oscillators, pulse calibrators, weather protection (including lightning protection) systems, seismometers, seismic amplifiers, and miscellaneous other parts. The central receiving and recording station upgrades will include discriminators, helicopter amplifier, omega timing system, strong motion instruments, wide-band velocity sensors, and other miscellaneous equipment.
Waste inspection tomography (WIT)
Bernardi, R.T.
1995-10-01
Waste Inspection Tomography (WIT) provides mobile semi-trailer mounted nondestructive examination (NDE) and assay (NDA) for nuclear waste drum characterization. WIT uses various computed tomography (CT) methods for both NDE and NDA of nuclear waste drums. Low level waste (LLW), transuranic (TRU), and mixed radioactive waste can be inspected and characterized without opening the drums. With externally transmitted x-ray NDE techniques, WIT has the ability to identify high density waste materials like heavy metals, define drum contents in two- and three-dimensional space, quantify free liquid volumes through density and x-ray attenuation coefficient discrimination, and measure drum wall thickness. With waste emitting gamma-ray NDA techniques, WIT can locate gamma emitting radioactive sources in two- and three-dimensional space, identify gamma emitting, isotopic species, identify the external activity levels of emitting gamma-ray sources, correct for waste matrix attenuation, provide internal activity approximations, and provide the data needed for waste classification as LLW or TRU.
High energy neutron Computed Tomography developed
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
objects. May 9, 2014 Neutron tomography horizontal "slice" of a tungsten and polyethylene test object containing tungsten carbide BBs. Neutron tomography horizontal "slice"...
T'Jampens, Stephane; /Orsay
2006-09-18
This thesis presents the full-angular time-dependent analysis of the vector-vector channel B{sub d}{sup 0} {yields} J/{psi}(K{sub S}{sup 0}{pi}{sup 0})*{sup 0}. After a review of the CP violation in the B meson system, the phenomenology of the charmonium-K*(892) channels is exposed. The method for the measurement of the transversity amplitudes of the B {yields} J/{psi}K*(892), based on a pseudo-likelihood method, is then exposed. The results from a 81.9 fb{sup -1} of collected data by the BABAR detector at the {Upsilon}(4S) resonance peak are |A{sub 0}|{sup 2} = 0.565 {+-} 0.011 {+-} 0.004, |A{sub {parallel}}|{sup 2} = 0.206 {+-} 0.016 {+-} 0.007, |A{sub {perpendicular}}|{sup 2} = 0.228 {+-} 0.016 {+-} 0.007, {delta}{sub {parallel}} = -2.766 {+-} 0.105 {+-} 0.040 and {delta}{sub {perpendicular}} = 2.935 {+-} 0.067 {+-} 0.040. Note that ({delta}{sub {parallel}}, {delta}{sub {perpendicular}}) {yields} (-{delta}{sub {parallel}}, {pi} - {delta}{sub {perpendicular}}) is also a solution. The strong phases {delta}{sub {parallel}} and {delta}{sub {perpendicular}} are at {approx}> 3{sigma} from {+-}{pi}, signing the presence of final state interactions and the breakdown of the factorization hypothesis. The forward-backward analysis of the K{pi} mass spectrum revealed the presence of a coherent S-wave interfering with the K*(892). It is the first evidence of this wave in the K{pi} system coming from a B meson. The particularity of the B{sub d}{sup 0} {yields} J/{psi}(K{sub S}{sup 0}{pi}{sup 0})*{sup 0} channel is to have a time-dependent but also an angular distribution which allows to measure sin 2{beta} but also cos2{beta}. The results from an unbinned maximum likelihood fit are sin 2{beta} = -0.10 {+-} 0.57 {+-} 0.14 and cos 2{beta} = 3.32{sub -0.96}{sup +0.76} {+-} 0.27 with the transversity amplitudes fixed to the values given above. The other solution for the strong phases flips the sign of cos 2{beta}. Theoretical considerations based on the s-quark helicity
Aubert, Bernard; Bona, M.; Karyotakis, Y.; Lees, J.P.; Poireau, V.; Prencipe, E.; Prudent, X.; Tisserand, V.; Garra Tico, J.; Grauges, E.; Lopez, L.; Palano, Antimo; Pappagallo, M.; Eigen, G.; Stugu, Bjarne; Sun, L.; Abrams, G.S.; Battaglia, M.; Brown, D.N.; Cahn, Robert N.; Jacobsen, R.G.; /LBL, Berkeley /Birmingham U. /Ruhr U., Bochum /Bristol U. /British Columbia U. /Brunel U. /Novosibirsk, IYF /UC, Irvine /UCLA /UC, Riverside /UC, San Diego /UC, Santa Barbara /UC, Santa Cruz /Caltech /Cincinnati U. /Colorado U. /Colorado State U. /Dortmund U. /Dresden, Tech. U. /Ecole Polytechnique /Edinburgh U. /Ferrara U. /INFN, Ferrara /Frascati /Genoa U. /INFN, Genoa /Harvard U. /Heidelberg U. /Humboldt U., Berlin /Imperial Coll., London /Iowa U. /Iowa State U. /Johns Hopkins U. /Karlsruhe U. /Orsay, LAL /LLNL, Livermore /Liverpool U. /Queen Mary, U. of London /Royal Holloway, U. of London /Louisville U. /Manchester U. /Maryland U. /Massachusetts U., Amherst /MIT /McGill U. /Consorzio Milano Ricerche /INFN, Milan /Mississippi U. /Montreal U. /Mt. Holyoke Coll. /Napoli Seconda U. /INFN, Naples /NIKHEF, Amsterdam /Notre Dame U. /Ohio State U. /Oregon U. /Padua U. /INFN, Padua /Paris U., VI-VII /Pennsylvania U. /Perugia U. /INFN, Perugia /INFN, Pisa /Princeton U. /Banca di Roma /Frascati /Rostock U. /Rutherford /DSM, DAPNIA, Saclay /South Carolina U. /SLAC /Stanford U., Phys. Dept. /SUNY, Albany /Tennessee U. /Texas U. /Texas U., Dallas /Turin U. /INFN, Turin /Trieste U. /INFN, Trieste /Valencia U., IFIC /Victoria U. /Warwick U. /Wisconsin U., Madison
2008-08-04
The authors present evidence of D{sup 0}-{bar D}{sup 0} mixing using a time-dependent amplitude analysis of the decay D{sup 0} {yields} K{sup +}{pi}{sup -}{pi}{sup 0} in a data sample of 384 fb{sup -1} collected with the BABAR detector at the PEP-II e{sup +}e{sup -} collider at SLAC. Assuming CP conservation, they measure the mixing parameters x{prime}{sub K{pi}{pi}{sup 0}} = [2.61{sub -0.68}{sup +0.57}(stat.) {+-} 0.39(syst.)]%, y{prime}{sub K{pi}{pi}{sup 0}} = [-0.06{sub -0.64}{sup +0.55}(stat.) {+-} 0.34(syst.)]%. The confidence level for the data to be consistent with the no-mixing hypothesis is 0.1%, including systematic uncertainties. This result is inconsistent with the no-mixing hypothesis with a significance of 3.2 standard deviations. They find no evidence of CP violation in mixing.
Zhao, Bin; Zhang, Dong-H.; Sun, Zhigang; Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei, Anhui 230026 ; Lee, Soo-Y.
2014-04-28
A framework for quantum state-to-state integral and differential cross sections of triatomic reactive scattering using the Multi-Configuration Time-Dependent Hartree (MCTDH) method is introduced, where a modified version of the Heidelberg MCTDH package is applied. Parity of the system is adopted using only non-negative helicity quantum numbers, which reduces the basis set size of the single particle functions in angular degree of freedom almost by half. The initial wave packet is constructed in the space-fixed frame, which can accurately account for the centrifugal potential. By using the reactant-coordinate-based method, the product state-resolved information can be accurately extracted. Test calculations are presented for the H + H{sub 2} reactive scattering. This work demonstrates the capability of the MCTDH method for extracting accurate state-to-state integral and differential cross sections. As an efficient scheme for high-dimensional problems, the MCTDH method may be promising for the study of product state-resolved cross sections for polyatomic reactive systems.
Baiardi, Alberto; Barone, Vincenzo; Bloino, Julien
2014-09-21
We present a new formulation of the time-dependent theory of Resonance-Raman spectroscopy (TD-RR). Particular attention has been devoted to the generality of the framework and to the possibility of including different effects (Duschinsky mixing, Herzberg-Teller contributions). Furthermore, the effects of different harmonic models for the intermediate electronic state are also investigated. Thanks to the implementation of the TD-RR procedure within a general-purpose quantum-chemistry program, both solvation and leading anharmonicity effects have been included in an effective way. The reliability and stability of our TD-RR implementation are validated against our previously proposed and well-tested time-independent procedure. Practical applications are illustrated with some closed- and open-shell medium-size molecules (anthracene, phenoxyl radical, benzyl radical) and the simulated spectra are compared to the experimental results. More complex and larger systems, not limited to organic compounds, can be also studied, as shown for the case of Tris(bipyridine)ruthenium(II) chloride.
Head-Marsden, Kade; Mazziotti, David A.
2015-02-07
For an open, time-dependent quantum system, Lindblad derived the most general modification of the quantum Liouville equation in the Markovian approximation that models environmental effects while preserving the non-negativity of the system’s density matrix. While Lindblad’s modification is correct for N-electron density matrices, solution of the Liouville equation with a Lindblad operator causes the one-electron reduced density matrix (1-RDM) to violate the Pauli exclusion principle. Consequently, after a short time, the 1-RDM is not representable by an ensemble N-electron density matrix (not ensemble N-representable). In this communication, we derive the necessary and sufficient constraints on the Lindbladian matrix within the Lindblad operator to ensure that the 1-RDM remains N-representable for all time. The theory is illustrated by considering the relaxation of an excitation in several molecules F{sub 2}, N{sub 2}, CO, and BeH{sub 2} subject to environmental noise.
Pautz, Andreas; Birkhofer, Adolf
2003-11-15
We introduce a new coupled neutronics/thermal hydraulics code system for analyzing transients of nuclear power plants and research reactors, based on a neutron transport theory approach. For the neutron kinetics, we have developed the code DORT-TD, a time-dependent extension of the well-known discrete ordinates code DORT. DORT-TD uses a fully implicit time integration scheme and is coupled via a general interface to the thermal hydraulics system code ATHLET, a generally applicable code for the analyses of LWR accident scenarios. Feedback is accounted for by interpolating multigroup cross sections from precalculated libraries, which are generated in advance for user-specified, discrete sets of thermal hydraulic parameters, e.g., fuel and coolant temperature. The coupled code system is applied to the high-flux research reactor FRM-II (Germany). Several design basis accidents are considered, namely the unintended control rod withdrawal, the loss of offsite power, and the loss of the secondary heat sink as well as a hypothetical transient with large reactivity insertion.
Xie, Wenbo; Liu, Lan; Sun, Zhigang; Guo, Hua; Dawes, Richard
2015-02-14
The title isotope exchange reaction was studied by converged time-dependent wave packet calculations, where an efficient 4th order split operator was applied to propagate the initial wave packet. State-to-state differential and integral cross sections up to the collision energy of 0.35 eV were obtained with {sup 32}O{sub 2} in the hypothetical j{sub 0} = 0 state. It is discovered that the differential cross sections are largely forward biased in the studied collision energy range, due to the fact that there is a considerable part of the reaction occurring with large impact parameter and short lifetime relative to the rotational period of the intermediate complex. The oscillations of the forward scattering amplitude as a function of collision energy, which result from coherent contribution of adjacent resonances, may be a sensitive probe for examining the quality of the underlying potential energy surface. A good agreement between the theoretical and recent experimental integral and differential cross sections at collision energy of 7.3 kcal/mol is obtained. However, the theoretical results predict slightly too much forward scattering and colder rotational distributions than the experimental observations at collision energy of 5.7 kcal/mol.
Sato, Shunsuke A.; Taniguchi, Yasutaka; Shinohara, Yasushi; Yabana, Kazuhiro
2015-12-14
We develop methods to calculate electron dynamics in crystalline solids in real-time time-dependent density functional theory employing exchange-correlation potentials which reproduce band gap energies of dielectrics; a meta-generalized gradient approximation was proposed by Tran and Blaha [Phys. Rev. Lett. 102, 226401 (2009)] (TBm-BJ) and a hybrid functional was proposed by Heyd, Scuseria, and Ernzerhof [J. Chem. Phys. 118, 8207 (2003)] (HSE). In time evolution calculations employing the TB-mBJ potential, we have found it necessary to adopt the predictor-corrector step for a stable time evolution. We have developed a method to evaluate electronic excitation energy without referring to the energy functional which is unknown for the TB-mBJ potential. For the HSE functional, we have developed a method for the operation of the Fock-like term in Fourier space to facilitate efficient use of massive parallel computers equipped with graphic processing units. We compare electronic excitations in silicon and germanium induced by femtosecond laser pulses using the TB-mBJ, HSE, and a simple local density approximation (LDA). At low laser intensities, electronic excitations are found to be sensitive to the band gap energy: they are close to each other using TB-mBJ and HSE and are much smaller in LDA. At high laser intensities close to the damage threshold, electronic excitation energies do not differ much among the three cases.
Seismic Velocity Measurements at Expanded Seismic Network Sites
Woolery, Edward W; Wang, Zhenming
2005-01-01
Structures at the Paducah Gaseous Diffusion Plant (PGDP), as well as at other locations in the northern Jackson Purchase of western Kentucky may be subjected to large far-field earthquake ground motions from the New Madrid seismic zone, as well as those from small and moderate-sized local events. The resultant ground motion a particular structure is exposed from such event will be a consequence of the earthquake magnitude, the structures' proximity to the event, and the dynamic and geometrical characteristics of the thick soils upon which they are, of necessity, constructed. This investigation evaluated the latter. Downhole and surface (i.e., refraction and reflection) seismic velocity data were collected at the Kentucky Seismic and Strong-Motion Network expansion sites in the vicinity of the Paducah Gaseous Diffusion Plant (PGDP) to define the dynamic properties of the deep sediment overburden that can produce modifying effects on earthquake waves. These effects are manifested as modifications of the earthquake waves' amplitude, frequency, and duration. Each of these three ground motion manifestations is also fundamental to the assessment of secondary earthquake engineering hazards such as liquefaction.
Hanford quarterly seismic monitoring report 96C
Reidel, S.P.
1996-09-24
Seismic monitoring at the Hanford Site was established in 1969 by the United States Geological Survey (USGS) under a contract with the U.S. Atomic Energy Commission. In 1975 the University of Washington assumed responsibility for and expanded the network. In 1979 the Basalt Waste Isolation Program (BWIP) became responsible for collecting seismic data for the site as part of site characterization. Rockwell International Operations followed by Westinghouse Hanford Company (WHC), Geosciences Group, operated the local network and were the contract technical advisors for the Eastern Washington Regional Network operated by the University of Washington. Funding ended for BWIP in December 1988. Seismic Monitoring and the University of Washington contract was then transferred WHC`s Environmental Division. Seismic Monitoring is currently assigned to WHC`s Hanford Technical Services (HTS), part of the Environmental Division. The Seismic Monitoring Analysis and Repair Team (SMART) operates, maintains, and analyzes data from the Hanford Seismic Network (HSN), extending the site historical seismic database and fulfilling U.S. Department of Energy, Richland Operations Office requirements and orders. The Seismic Monitoring Analysis and Repair Team also maintains the Eastern Washington Regional Network (EWRN). The University of Washington uses the data from the EWRN and other seismic networks in the Northwest to provide the SMART with necessary regional input for the seismic hazards analysis at the Hanford Site.
Seismic Isolation Working Meeting Gap Analysis Report
Justin Coleman; Piyush Sabharwall
2014-09-01
The ultimate goal in nuclear facility and nuclear power plant operations is operating safety during normal operations and maintaining core cooling capabilities during off-normal events including external hazards. Understanding the impact external hazards, such as flooding and earthquakes, have on nuclear facilities and NPPs is critical to deciding how to manage these hazards to expectable levels of risk. From a seismic risk perspective the goal is to manage seismic risk. Seismic risk is determined by convolving the seismic hazard with seismic fragilities (capacity of systems, structures, and components (SSCs)). There are large uncertainties associated with evolving nature of the seismic hazard curves. Additionally there are requirements within DOE and potential requirements within NRC to reconsider updated seismic hazard curves every 10 years. Therefore opportunity exists for engineered solutions to manage this seismic uncertainty. One engineered solution is seismic isolation. Current seismic isolation (SI) designs (used in commercial industry) reduce horizontal earthquake loads and protect critical infrastructure from the potentially destructive effects of large earthquakes. The benefit of SI application in the nuclear industry is being recognized and SI systems have been proposed, in the American Society of Civil Engineers (ASCE) 4 standard, to be released in 2014, for Light Water Reactors (LWR) facilities using commercially available technology. However, there is a lack of industry application to the nuclear industry and uncertainty with implementing the procedures outlined in ASCE-4. Opportunity exists to determine barriers associated with implementation of current ASCE-4 standard language.
Nuclear component horizontal seismic restraint
Snyder, Glenn J.
1988-01-01
A nuclear component horizontal seismic restraint. Small gaps limit horizontal displacement of components during a seismic occurrence and therefore reduce dynamic loadings on the free lower end. The reactor vessel and reactor guard vessel use thicker section roll-forged rings welded between the vessel straight shell sections and the bottom hemispherical head sections. The inside of the reactor guard vessel ring forging contains local vertical dovetail slots and upper ledge pockets to mount and retain field fitted and installed blocks. As an option, the horizontal displacement of the reactor vessel core support cone can be limited by including shop fitted/installed local blocks in opposing alignment with the reactor vessel forged ring. Beams embedded in the wall of the reactor building protrude into apertures in the thermal insulation shell adjacent the reactor guard vessel ring and have motion limit blocks attached thereto to provide to a predetermined clearance between the blocks and reactor guard vessel ring.
Korneev, Valeri A [LaFayette, CA
2009-05-05
The detailed analysis of cross well seismic data for a gas reservoir in Texas revealed two newly detected seismic wave effects, recorded approximately 2000 feet above the reservoir. A tube-wave (150) is initiated in a source well (110) by a source (111), travels in the source well (110), is coupled to a geological feature (140), propagates (151) through the geological feature (140), is coupled back to a tube-wave (152) at a receiver well (120), and is and received by receiver(s) (121) in either the same (110) or a different receiving well (120). The tube-wave has been shown to be extremely sensitive to changes in reservoir characteristics. Tube-waves appear to couple most effectively to reservoirs where the well casing is perforated, allowing direct fluid contact from the interior of a well case to the reservoir.
Korneev, Valeri A.; Bakulin, Andrey
2009-10-13
The detailed analysis of cross well seismic data for a gas reservoir in Texas revealed two newly detected seismic wave effects, recorded approximately 2000 feet above the reservoir. A tube-wave (150) is initiated in a source well (110) by a source (111), travels in the source well (110), is coupled to a geological feature (140), propagates (151) through the geological feature (140), is coupled back to a tube-wave (152) at a receiver well (120), and is and received by receiver(s) (121) in either the same (110) or a different receiving well (120). The tube-wave has been shown to be extremely sensitive to changes in reservoir characteristics. Tube-waves appear to couple most effectively to reservoirs where the well casing is perforated, allowing direct fluid contact from the interior of a well case to the reservoir.
seismic | National Nuclear Security Administration
National Nuclear Security Administration (NNSA)
seismic What rock has the biggest impact on the San Andreas? The moon. Last year the subject of a Hollywood disaster movie, this year the San Andreas Fault is getting attention from geophysicists at NNSA's Los Alamos National Laboratory. U.S. Geological Survey and Los Alamos scientists have released a new study detailing how the gravitational tug of the sun and... NNSA administrator visits NNSS to meet team, see national security work Last month, Department of Energy Under Secretary for Nuclear
Kuiken, Benjamin E. Van; Valiev, Marat; Daifuku, Stephanie L.; Bannan, Caitlin; Strader, Matthew L.; Cho, Hana; Huse, Nils; Schoenlein, Robert W.; Govind, Niranjan; Khalil, Munira
2013-04-26
Ruthenium L3-edge X-ray absorption (XA) spectroscopy probes unoccupied 4d orbitals of the metal atom and is increasingly being used to investigate the local electronic structure in ground and excited electronic states of Ru complexes. The simultaneous development of computational tools for simulating Ru L3-edge spectra is crucial for interpreting the spectral features at a molecular level. This study demonstrates that time-dependent density functional theory (TDDFT) is a viable and predictive tool for simulating ruthenium L3-edge XA spectroscopy. We systematically investigate the effects of exchange correlation functional and implicit and explicit solvent interactions on a series of RuII and RuIII complexes in their ground and electronic excited states. The TDDFT simulations reproduce all of the experimentally observed features in Ru L3-edge XA spectra within the experimental resolution (0.4 eV). Our simulations identify ligand-specific charge transfer features in complicated Ru L3-edge spectra of [Ru(CN)6]4- and RuII polypyridyl complexes illustrating the advantage of using TDDFT in complex systems. We conclude that the B3LYP functional most accurately predicts the transition energies of charge transfer features in these systems. We use our TDDFT approach to simulate experimental Ru L3-edge XA spectra of transition metal mixed-valence dimers of the form [(NC)5MII-CN-RuIII(NH3)5] (where M = Fe or Ru) dissolved in water. Our study determines the spectral signatures of electron delocalization in Ru L3-edge XA spectra. We find that the inclusion of explicit solvent molecules is necessary for reproducing the spectral features and the experimentally determined valencies in these mixed-valence complexes. This study validates the use of TDDFT for simulating Ru 2p excitations using popular quantum chemistry codes and providing a powerful interpretive tool for equilibrium and ultrafast Ru L3-edge XA spectroscopy.