Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Optimal Distributed Control of Spatially Invariant Time Varying Linear Systems

In this paper, the optimal distributed linear time varying H2 control problem for cone causal spatially invariant linear time varying (SILTV) systems is considered. This class of systems was initially defined in the work of Voulgaris et al. for cone causal spatially invariant (SI) linear time invariant (LTI) systems. First, the optimal time varying H2 problem is posed using a version of the Youla parameterization. This allows to transform the problem into an affine form which results in a convex but infinite dimensional distributed optimization. Next, the problem is solved by computing an operator projection on a class of time varying Youla parameters with a cone causal distributed structure. The cone causal spatially invariant systems are viewed as input-output multiplicative operators with a mixed causal time varying structure with respect to time, and spatially invariant with respect to the space dimension.

Djouadi, Seddik M.↗
DOE OSTI2022

An Alternative Timing and Synchronization Approach for Situational Awareness and Predictive Analytics

Accurate and synchronized timing information is required by power system operators for controlling the grid infrastructure (relays, Phasor Measurement Units (PMUs), etc.) and determining asset positions. Satellite-based global positioning system (GPS) is the primary source of timing information. However, GPS disruptions today (both intentional and unintentional) can significantly compromise the reliability and security of our electric grids. A robust alternate source for accurate timing is critical to serve both as a deterrent against malicious attacks and as a redundant system in enhancing the resilience against extreme events that could disrupt the GPS network. To achieve this, we rely on the highly accurate, terrestrial atomic clock-based network for alternative timing and synchronization. In this paper, we discuss an experimental setup for an alternative timing approach. The data obtained from this experimental setup is continuously monitored and analyzed using various time deviation metrics. We also use these metrics to compute deviations of our clock with respect to the National Institute of Standards and Technologys (NIST) GPS data. The results obtained from these metric computations are elaborately discussed. Finally, we discuss the integration of the procedures involved, like real-time data ingestion, metric computation, and result visualization, in a novel microservices-based architecture for situational awareness.

Chinthavali, Supriya↗
DOE OSTI2021

Precise Timing Based on Pulsar Observation for Grid Synchronization

While Global Positioning System signals are widely used as the synchronized timing Sources for wide-area measurement systems for power grid monitoring and control, they are vulnerable to failures and malicious attacks. To address this Problem, an alternative kind of timing sources, the millisecond pulsars, were proposed to serve the power grid as they are more precise and physically indestructible. Millisecond pulsars are a kind of neutron stars that emit radio pulse signals that have extremely stable periods of milliseconds. However, the observed pulse Signals from pulsars cannot be directly used for timing because of their low signal-to-noise ratio and the dispersion effects of interstellar medium. In this paper, the precise timing method based on pulsar observation data is proposed to interpret the pulsar observation data and therefrom to generate precise timing signal for the power grid synchronization. In our proposed method,. incoherent de-dispersion is used on the observation data to mitigate the effect of interstellar medium. dispersion. Then spectrum analysis and segment folding are used to precisely estimate the local pulsar period. Finally the standard pulse per second signal is generated from. the estimated pulsar period and its precision. is analyzed. Experiments with the observation data of the pulsar based timing signal can be more accurate compared to GPS timing signals.

Luo, Xiqian↗
DOE OSTI2012

Toward Real Time Data Analysis for Smart Grids

This paper describes the architecture and design of a novel system for supporting large-scale real-time data analysis for future power grid systems. The widespread deployment of renewable generation, smart grid controls, energy storage, plug-in hybrids, and new conducting materials will require fundamental changes in the operational concepts and principal components of the grid. As a result, the whole system becomes highly dynamic and requires constant adjusting based on real time data. Even though millions of sensors such as phase measurement units (PMU) and smart meters are being widely deployed, a data layer that can analyze this amount of data in real time is needed. Unlike the data fabric in other cloud services, the data layer for smart grids has some unique design requirements. First, this layer must provide real time guarantees. Second, this layer must be scalable to allow a large number of applications to access the data from millions of sensors in real time. Third, reliability is critical and this layer must be able to continue to provide service in face of failures. Fourth, this layer must be secure. We address these challenges though a scalable system architecture that integrates the I/O and data processing capability in a devise set of devices. Data process operations can be placed anywhere from sensors, data storage devices, to control centers. We further employ compression to improve performance. We design a lightweight compression customized for power grid data. Our system can reduce end-to-end response time by reduce I/O overhead through compression and overlap compression operations with I/O. The initial prototype of our system was demonstrated with several use cases from PNNL’s FPGI and show that our system can provide real time guarantees to a diverse set of applications.

Data Analytics, Power Grid, Big Data↗
DOE OSTI2016

Multi-time Scale Coordination of Distributed Energy Resources in Isolated Power Systems

In isolated power systems, including microgrids, distributed assets, such as renewable energy resources (e.g. wind, solar) and energy storage, can be actively coordinated to reduce dependency on fossil fuel generation. The key challenge of such coordination arises from significant uncertainty and variability occurring at small time scales associated with increased penetration of renewables. Specifically, the problem is with ensuring economic and efficient utilization of DERs, while also meeting operational objectives such as adequate frequency performance. One possible solution is to reduce the time step at which tertiary controls are implemented and to ensure feedback and look-ahead capability are incorporated to handle variability and uncertainty. However, reducing the time step of tertiary controls necessitates investigating time-scale coupling with primary controls so as not to exacerbate system stability issues. In this paper, an optimal coordination (OC) strategy, which considers multiple time-scales, is proposed for isolated microgrid systems with a mix of DERs. This coordination strategy is based on an online moving horizon optimization approach. The effectiveness of the strategy was evaluated in terms of economics, technical performance, and computation time by varying key parameters that significantly impact performance. The illustrative example with realistic scenarios on a simulated isolated microgrid test system suggests that the proposed approach is generalizable towards designing multi-time scale optimal coordination strategies for isolated power systems.

Mayhorn, Ebony↗
DOE OSTI2014

Effect of UV curing time on physical and electrical properties and reliability of low dielectric constant materials

This study comprehensively investigates the effect of ultraviolet (UV) curing time on the physical, electrical, and reliability characteristics of porous low-k materials. Following UV irradiation for various periods, the depth profiles of the chemical composition in the low-k dielectrics were homogeneous. Initially, the UV curing process preferentially removed porogen-related CH{sub x} groups and then modified Si-CH{sub 3} and cage Si-O bonds to form network Si-O bonds. The lowest dielectric constant (k value) was thus obtained at a UV curing time of 300 s. Additionally, UV irradiation made porogen-based low-k materials hydrophobic and to an extent that increased with UV curing time. With a short curing time (<300 s), porogen was not completely removed and the residues degraded reliability performance. A long curing time (>300 s) was associated with improved mechanical strength, electrical performance, and reliability of the low-k materials, but none of these increased linearly with UV curing time. Therefore, UV curing is necessary, but the process time must be optimized for porous low-k materials on back-end of line integration in 45 nm or below technology nodes.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Real-time soft tissue motion estimation for lung tumors during radiotherapy delivery

Purpose: To provide real-time lung tumor motion estimation during radiotherapy treatment delivery without the need for implanted fiducial markers or additional imaging dose to the patient.Methods: 2D radiographs from the therapy beam's-eye-view (BEV) perspective are captured at a frame rate of 12.8 Hz with a frame grabber allowing direct RAM access to the image buffer. An in-house developed real-time soft tissue localization algorithm is utilized to calculate soft tissue displacement from these images in real-time. The system is tested with a Varian TX linear accelerator and an AS-1000 amorphous silicon electronic portal imaging device operating at a resolution of 512 × 384 pixels. The accuracy of the motion estimation is verified with a dynamic motion phantom. Clinical accuracy was tested on lung SBRT images acquired at 2 fps.Results: Real-time lung tumor motion estimation from BEV images without fiducial markers is successfully demonstrated. For the phantom study, a mean tracking error <1.0 mm [root mean square (rms) error of 0.3 mm] was observed. The tracking rms accuracy on BEV images from a lung SBRT patient (≈20 mm tumor motion range) is 1.0 mm.Conclusions: The authors demonstrate for the first time real-time markerless lung tumor motion estimation from BEV images alone. The described system can operate at a frame rate of 12.8 Hz and does not require prior knowledge to establish traceable landmarks for tracking on the fly. The authors show that the geometric accuracy is similar to (or better than) previously published markerless algorithms not operating in real-time.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

SU-D-18C-05: Variable Bolus Arterial Spin Labeling MRI for Accurate Cerebral Blood Flow and Arterial Transit Time Mapping

Purpose: Arterial spin labeling (ASL) is an MRI perfusion imaging method from which quantitative cerebral blood flow (CBF) maps can be calculated. Acquisition with variable post-labeling delays (PLD) and variable TRs allows for arterial transit time (ATT) mapping and leads to more accurate CBF quantification with a scan time saving of 48%. In addition, T1 and M0 maps can be obtained without a separate scan. In order to accurately estimate ATT and T1 of brain tissue from the ASL data, variable labeling durations were invented, entitled variable-bolus ASL. Methods: All images were collected on a healthy subject with a 3T Siemens Skyra scanner. Variable-bolus Psuedo-continuous ASL (PCASL) images were collected with 7 TI times ranging 100-4300ms in increments of 700ms with TR ranging 1000-5200ms. All boluses were 1600ms when the TI allowed, otherwise the bolus duration was 100ms shorter than the TI. All TI times were interleaved to reduce sensitivity to motion. Voxel-wise T1 and M0 maps were estimated using a linear least squares fitting routine from the average singal from each TI time. Then pairwise subtraction of each label/control pair and averaging for each TI time was performed. CBF and ATT maps were created using the standard model by Buxton et al. with a nonlinear fitting routine using the T1 tissue map. Results: CBF maps insensitive to ATT were produced along with ATT maps. Both maps show patterns and averages consistent with literature. The T1 map also shows typical T1 contrast. Conclusion: It has been demonstrated that variablebolus ASL produces CBF maps free from the errors due to ATT and tissue T1 variations and provides M0, T1, and ATT maps which have potential utility. This is accomplished with a single scan in a feasible scan time (under 6 minutes) with low sensivity to motion.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-T-39: Analyses of Plan Quality Using Different Gantry Rotation Time On Volumetric Modulated Arc Therapy in Rectal Carcinoma

Purpose: To compare plan quality of volumetric modulated arc radiotherapy (VMAT) using different gantry rotation times. Methods: In ten rectal carcinoma patients, plans with different gantry rotation times were designed by Oncentra4.3 planning system, each of them was optimized by experienced planners using almost similar parameters and dose constraints except the different gantry rotation time. The gantry rotation time was increased from 30s to 150s by step of 30(30s, 60s, 90s, 120s, and 150s). Plans (VMAT30s, VMAT60s, VMAT90s, VMAT120s and VMAT150s) were normalized so that at least 95% of PTV would receive the prescription dose 50 Gray in 25 fractions. Dose-volume histograms and isodose distributions were evaluated. Conformity indices (CI), homogeneous indices (HI), monitor units (MUs), delivery time were also compared. Results: As the increasing of gantry rotation time, VMAT150s provided comparable organs at risk sparing and better target coverage and conformity than VMAT30s. In the normal tissue such as bladder and femurs, the VMAT plans show almost the same V30 and Dmean for them compared with IMRT plans, meanwhile the conformity indices (CI) was better than IMRT plans(p=0.021). However, there was no significant different among the VMAT60s, VMAT90s, VMAT120s and VMAT150s (p=0.054). VMAT150s increased delivery time and monitor units in plans without improving plan quality compared to VMAT60s, VMAT90s and VMAT120s. Conclusion: VMAT60s, VMAT90s, VMAT120s and VMAT150s achieved comparable treatment plans in rectal carcinoma. Slower gantry movement in VMAT150s only improves slight improvement than the others.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

WE-D-17A-03: Improvement of Accuracy of Spot-Scanning Proton Beam Delivery for Liver Tumor by Real-Time Tumor-Monitoring and Gating System: A Simulation Study

Purpose: To improve the accuracy of spot-scanning proton beam delivery for target in motion, a real-time tumor-monitoring and gating system using fluoroscopy images was developed. This study investigates the efficacy of this method for treatment of liver tumors using simulation. Methods: Three-dimensional position of a fiducial marker inserted close to the tumor is calculated in real time and proton beam is gated according to the marker's distance from the planned position (Shirato, 2012). The efficient beam delivery is realized even for the irregular and sporadic motion signals, by employing the multiple-gated irradiations per operation cycle (Umezawa, 2012). For each of two breath-hold CTs (CTV=14.6cc, 63.1cc), dose distributions were calculated with internal margins corresponding to freebreathing (FB) and real-time gating (RG) with a 2-mm gating window. We applied 8 trajectories of liver tumor recorded during the treatment of RTRT in X-ray therapy and 6 initial timings. Dmax/Dmin in CTV, mean liver dose (MLD), and irradiation time to administer 3 Gy (RBE) dose were estimated assuming rigid motion of targets by using in-house simulation tools and VQA treatment planning system (Hitachi, Ltd., Tokyo). Results: Dmax/Dmin was degraded by less than 5% compared to the prescribed dose with all motion parameters for smaller CTV and less than 7% for larger CTV with one exception. Irradiation time showed only a modest increase if RG was used instead of FB; the average value over motion parameters was 113 (FB) and 138 s (RG) for smaller CTV and 120 (FB) and 207 s (RG) for larger CTV. In RG, it was within 5 min for all but one trajectory. MLD was markedly decreased by 14% and 5–6% for smaller and larger CTVs respectively, if RG was applied. Conclusions: Spot-scanning proton beam was shown to be delivered successfully to liver tumor without much lengthening of treatment time. This research was supported by the Cabinet Office, Government of Japan and the Japan Society for the Promotion of Science (JSPS) through the Funding Program for World-Leading Innovative R and D on Science and Technology (FIRST Program), initiated by the Council for Science and Technology Policy (CSTP)

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-P-27: Efficient Process for AccuBoost Planning and Treatment Delivery to Minimize Patient Compression Time

Purpose: To implement an efficient and robust process for AccuBoost planning and treatment delivery that can be safely performed by a single Physicist while minimizing patient’s total session time. Methods: Following a thorough commissioning and validation process, templates were created in the brachytherapy planning system for each AccuBoost applicator. Tables of individual and total nominal dwell times for each applicator as a function of separation were generated to streamline planning while an Excel-based nomogram provided by the vendor functions as a secondary verification of the treatment parameters. Tables of surface dose as a function of separation and applicator, along with concise guidance documents for applicator selection, are readily available during the planning process. The entire process is described in a set of detailed Standard Operating Procedures which, in addition to the items described above, include a verbal time-out between the primary planner and the individual performing the secondary verification as well as direct visual confirmation of applicator placement using an articulated mirror. Prior to treatment initiation, a final time-out is conducted with the Radiation Oncologist. Chart documentation is finalized after the patient is released from compression following completion of the treatment. Results: With the aforementioned procedures, it has been possible to consistently limit the time required to prepare each treatment such that the patient is typically under compression for less than 10 minutes per orientation prior to the initiation of the treatment, which is particularly important for APBI cases. This process can be overseen by a single physicist assisted by a dosimetrist and has been optimized during the past 16 months, with 180 treatment sessions safely completed to date. Conclusion: This work demonstrates the implementation of an efficient and robust process for real-time-planned AccuBoost treatments that effectively minimizes the time that the patient is under compression.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Fast leaf-fitting with generalized underdose/overdose constraints for real-time MLC tracking

Purpose: Real-time multileaf collimator (MLC) tracking is a promising approach to the management of intrafractional tumor motion during thoracic and abdominal radiotherapy. MLC tracking is typically performed in two steps: transforming a planned MLC aperture in response to patient motion and refitting the leaves to the newly generated aperture. One of the challenges of this approach is the inability to faithfully reproduce the desired motion-adapted aperture. This work presents an optimization-based framework with which to solve this leaf-fitting problem in real-time. Methods: This optimization framework is designed to facilitate the determination of leaf positions in real-time while accounting for the trade-off between coverage of the PTV and avoidance of organs at risk (OARs). Derived within this framework, an algorithm is presented that can account for general linear transformations of the planned MLC aperture, particularly 3D translations and in-plane rotations. This algorithm, together with algorithms presented in Sawant et al. [“Management of three-dimensional intrafraction motion through real-time DMLC tracking,” Med. Phys. 35, 2050–2061 (2008)] and Ruan and Keall [Presented at the 2011 IEEE Power Engineering and Automation Conference (PEAM) (2011) (unpublished)], was applied to apertures derived from eight lung intensity modulated radiotherapy plans subjected to six-degree-of-freedom motion traces acquired from lung cancer patients using the kilovoltage intrafraction monitoring system developed at the University of Sydney. A quality-of-fit metric was defined, and each algorithm was evaluated in terms of quality-of-fit and computation time. Results: This algorithm is shown to perform leaf-fittings of apertures, each with 80 leaf pairs, in 0.226 ms on average as compared to 0.082 and 64.2 ms for the algorithms of Sawant et al., Ruan, and Keall, respectively. The algorithm shows approximately 12% improvement in quality-of-fit over the Sawant et al. approach, while performing comparably to Ruan and Keall. Conclusions: This work improves upon the quality of the Sawant et al. approach, but does so without sacrificing run-time performance. In addition, using this framework allows for complex leaf-fitting strategies that can be used to account for PTV/OAR trade-off during real-time MLC tracking.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2016

SU-F-J-54: Towards Real-Time Volumetric Imaging Using the Treatment Beam and KV Beam

Purpose: Existing real-time imaging uses dual (orthogonal) kV beam fluoroscopies and may result in significant amount of extra radiation to patients, especially for prolonged treatment cases. In addition, kV projections only provide 2D information, which is insufficient for in vivo dose reconstruction. We propose real-time volumetric imaging using prior knowledge of pre-treatment 4D images and real-time 2D transit data of treatment beam and kV beam. Methods: The pre-treatment multi-snapshot volumetric images are used to simulate 2D projections of both the treatment beam and kV beam, respectively, for each treatment field defined by the control point. During radiation delivery, the transit signals acquired by the electronic portal image device (EPID) are processed for every projection and compared with pre-calculation by cross-correlation for phase matching and thus 3D snapshot identification or real-time volumetric imaging. The data processing involves taking logarithmic ratios of EPID signals with respect to the air scan to reduce modeling uncertainties in head scatter fluence and EPID response. Simulated 2D projections are also used to pre-calculate confidence levels in phase matching. Treatment beam projections that have a low confidence level either in pre-calculation or real-time acquisition will trigger kV beams so that complementary information can be exploited. In case both the treatment beam and kV beam return low confidence in phase matching, a predicted phase based on linear regression will be generated. Results: Simulation studies indicated treatment beams provide sufficient confidence in phase matching for most cases. At times of low confidence from treatment beams, kV imaging provides sufficient confidence in phase matching due to its complementary configuration. Conclusion: The proposed real-time volumetric imaging utilizes the treatment beam and triggers kV beams for complementary information when the treatment beam along does not provide sufficient confidence for phase matching. This strategy minimizes the use of extra radiation to patients. This project is partially supported by a Varian MRA grant.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

MO-DE-207A-08: Four-Dimensional Cone-Beam CT Iterative Reconstruction with Time-Ordered Chain Graph Model for Non-Periodic Organ Motion and Deformation

Purpose: The purpose of this study is to propose a new concept of four-dimensional (4D) cone-beam CT (CBCT) reconstruction for non-periodic organ motion using the Time-ordered Chain Graph Model (TCGM), and to compare the reconstructed results with the previously proposed methods, the total variation-based compressed sensing (TVCS) and prior-image constrained compressed sensing (PICCS). Methods: CBCT reconstruction method introduced in this study consisted of maximum a posteriori (MAP) iterative reconstruction combined with a regularization term derived from a concept of TCGM, which includes a constraint coming from the images of neighbouring time-phases. The time-ordered image series were concurrently reconstructed in the MAP iterative reconstruction framework. Angular range of projections for each time-phase was 90 degrees for TCGM and PICCS, and 200 degrees for TVCS. Two kinds of projection data, an elliptic-cylindrical digital phantom data and two clinical patients’ data, were used for reconstruction. The digital phantom contained an air sphere moving 3 cm along longitudinal axis, and temporal resolution of each method was evaluated by measuring the penumbral width of reconstructed moving air sphere. The clinical feasibility of non-periodic time-ordered 4D CBCT reconstruction was also examined using projection data of prostate cancer patients. Results: The results of reconstructed digital phantom shows that the penumbral widths of TCGM yielded the narrowest result; PICCS and TCGM were 10.6% and 17.4% narrower than that of TVCS, respectively. This suggests that the TCGM has the better temporal resolution than the others. Patients’ CBCT projection data were also reconstructed and all three reconstructed results showed motion of rectal gas and stool. The result of TCGM provided visually clearer and less blurring images. Conclusion: The present study demonstrates that the new concept for 4D CBCT reconstruction, TCGM, combined with MAP iterative reconstruction framework enables time-ordered image reconstruction with narrower time-window.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

TH-CD-209-11: Simulation Study of Real-Time-Image Gating On Spot Scanning Proton Therapy for Lung Tumors

Purpose: To study the impact of a real-time-image gating on spot scanning proton therapy for lung tumors and to examine the suitable size of the gating window (GW). Methods: We investigated a real-time-image gated proton therapy (RGPT), in which two fluoroscopic units monitor a gold sphere fiducial in real-time, and the proton beam is irradiated only when the marker enters within the pre-assigned GW. We designed 5 treatment plans for 7 lung cancer patients: RGPT with a GW of ±1, 2, 5, and 8 mm and free-breathing proton therapy (FBPT) using the end-exhale and average images of 4-dimensional (4D) CT, respectively. 70 Gy(RBE)/10fr was prescribed to 99% of the targets. The time-series data of the three-dimensional marker positions (RTRT data) were grouped into 10 phases to associate with the phases of 4DCT. The 4D dose distributions were calculated using the plan information, RTRT Data, 4DCT, and modeled accelerator pattern. The dose distribution in each respiratory phase was deformed into the end-exhale CT. The D99 and D5-95 of CTV (with a criteria of D99>95% and D5-95<5%), V20 of Lung-GTV, and treatment times were evaluated. Results: GWs ≤ ±2 mm satisfied the criteria of CTV in all cases, whereas GWs ≥ ±5 mm did not satisfy the criteria in some cases. The V20 was reduced by more than 18.9% (relative to FBPT) for GW ≤ ±2 mm, but equaled or even surpassed the FBPT for GWs ≥ ±5 mm. The irradiation times for the ±1, 2, 5, and 8 mm GWs and FBPT were 372.4±208.3, 215.2±51.5, 180.9±31.6, 178.4±21.2, and 140.1±15.2 s, respectively. The GW of ±1 mm caused large variation in irradiation time among the patients. Conclusion: In RGPT for lung cancer, the most suitable GW, in terms of good dose preservation without prolonging the therapeutic beam delivery, is ±2 mm.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

Study of the Formation Time of a Self-Sustained Subnanosecond Discharge at High and Ultrahigh Gas Pressures

The formation times of self-sustained subnanosecond discharges in nitrogen at pressures of 1‒40 atm and in hydrogen at pressures of 1–60 atm are analyzed in terms of the avalanche model. In experiments, a subnanosecond voltage pulse with an amplitude of 102 ± 2 kV was applied to a 0.5-mm-long discharge gap with a uniformly distributed electric field (the curvature radii of both the cathode and anode ends were 1 cm). The rise time of the voltage pulse from 0.1 to 0.9 of its amplitude value was about 250 ps. Breakdown occurred at the leading edge of the pulse. The discharge formation time was measured at different gas pressures with a step of 5–10 atm. Analysis of the experimental results shows that, in nitrogen at pressures of 10–40 atm and in hydrogen at pressures of 20–50 atm, breakdown occurs earlier than the electron avalanche reaches its critical length and that the critical avalanche length lies in the range of (2–8) × 10{sup –2} mm, which is one order of magnitude shorter than the discharge gap length. This means that the avalanche–streamer model is inapplicable in this case. The fast formation of a conducting channel under these conditions can be explained by ionization of gas by runaway electrons. In this case, the conducting column develops as a result of simultaneous development of a large number of electron avalanches in the gas volume. An increase in the hydrogen pressure from 50 to 60 atm leads to an abrupt increase in the discharge formation time by about 50%. As a result, the growth time of the electron avalanche to its critical length becomes shorter than the discharge formation time. In this case, the electrons cease to pass into the runaway regime and the discharge is initiated from the cathode due to field emission from microinhomogeneities on its surface. Under these conditions, the discharge formation time is well described by the avalanche–streamer model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Fast Solution of Fully Implicit Runge--Kutta and Discontinuous Galerkin in Time for Numerical PDEs, Part I: the Linear Setting

Fully implicit Runge--Kutta (IRK) methods have many desirable properties as time integration schemes in terms of accuracy and stability, but high-order IRK methods are not commonly used in practice with numerical PDEs due to the difficulty of solving the stage equations. This paper introduces a theoretical and algorithmic preconditioning framework for solving the systems of equations that arise from IRK methods applied to linear numerical PDEs (without algebraic constraints). Additionally, this framework also naturally applies to discontinuous Galerkin discretizations in time. Under quite general assumptions on the spatial discretization that yield stable time integration, the preconditioned operator is proven to have condition number bounded by a small, order-one constant, independent of the spatial mesh and time-step size, and with only weak dependence on number of stages/polynomial order; for example, the preconditioned operator for 10th-order Gauss IRK has condition number less than two, independent of the spatial discretization and time step. The new method can be used with arbitrary existing preconditioners for backward Euler-type time-stepping schemes and is amenable to the use of three-term recursion Krylov methods when the underlying spatial discretization is symmetric. The new method is demonstrated to be effective on various high-order finite-difference and finite element discretizations of linear parabolic and hyperbolic problems, demonstrating fast, scalable solution of up to 10th-order accuracy. The new method consistently outperforms existing block preconditioning approaches, and in several cases, the new method can achieve 4th-order accuracy using Gauss integration with roughly half the number of preconditioner applications and wallclock time as required using standard diagonally IRK methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Code for Calculating Regional Seismic Travel Time

The RSTT software computes predictions of the travel time of seismic energy traveling from a source to a receiver through 2.5D models of the seismic velocity distribution within the Earth. The two primary applications for the RSTT library are tomographic inversion studies and seismic event location calculations. In tomographic inversions studies, a seismologist begins with number of source-receiver travel time observations and an initial starting model of the velocity distribution within the Earth. A forward travel time calculator, such as the RSTT library, is used to compute predictions of each observed travel time and all of the residuals (observed minus predicted travel time) are calculated. The Earth model is then modified in some systematic way with the goal of minimizing the residuals. The Earth model obtained in this way is assumed to be a better model than the starting model if it has lower residuals. The other major application for the RSTT library is seismic event location. Given an Earth model, an initial estimate of the location of a seismic event, and some number of observations of seismic travel time thought to have originated from that event, location codes systematically modify the estimate of the location of the event with the goal of minimizing the difference between the observed and predicted travel times. The second application, seismic event location, is routinely implemented by the military as part of its effort to monitor the Earth for nuclear tests conducted by foreign countries.

BALLARD, SANFORD↗