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,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

SSMI Wind Speed Climatology of the Time of Monsoon Wind Offset in the Western Arabian Sea

Forecasting the time of onset of monsoon wind in the western Arabian Sea, which is believed to precede the onset of rainfall along the west coast of India, is an important unsolved problem. Prior to measurements of the surface wind field by satellite, there was an absence of suitable surface wind observations. NASA scatterometer (NSCAT) surface wind vectors revealed that the time of the 1997 onset of 12 m/s southwest monsoon wind speeds in the western Arabian Sea preceded the onset of monsoon rainfall in Goa, India, by 3 - 4 days. Wind speed and direction data were necessary to establish a dynamical mechanism between times of onset of 12 m/s wind speed off Somalia and rainfall in Goa. Except for NSCAT, no satellite scatterometer wind product recorded adequately sampled 2-day 1deg x 1deg averaged wind vectors, which are the required space and time scales, to examine the wind-rain relationship in other years. However, the greater-than-95% steadiness of summer monsoon winds allows an opportunity to use satellite measurements of surface wind speed. The Special Sensor Microwave Imager (SSMI) recorded surface wind speed with adequate sampling to produce a 1-day, 1deg x 1deg data product during 1988 - 1998. SSMI data had been uniformly processed throughout the period. Times of onset of 12 m/s wind speed off Somalia determined with the SSMI data set were 21 May 1988, 24 May 1989, 17 May 1990, 28 May 1991, 8 June 1992, 28 May 1993, 30 May 1994, 7 June 1995, 29 May 1996, 12 June 1997, and 15 May 1998. Uncertainty of the 1992 and 1996 times of onset were increased because of the absence of SSMI data on 6 and 7 June 1992 and on 30 May 1996. Correlations of timing of monsoon wind onset with El Nino will be described. Variability of the time difference between times of onset of 12 m/s wind speed and Goa rainfall will be discussed. At the time of submission of the abstract, the Goa rainfall data have not arrived from the India Meteorological Department.

Halpern, David↗
NASA NTRS2003

On the Computation of Space-Time Correlations by Large-Eddy Simulation

Numerical comparisons in decaying isotropic turbulence suggest that there exist discrepancies in time correlations evaluated by DNS and LES using eddy-viscosity-type SGS models. This is consistent with the previous observations in forced isotropic turbulence. Therefore, forcing is not the main cause of the discrepancies. Comparisons among different SGS models in the LES also indicate that the model choice affects the time correlations in the LES. The multi-scale LES method using the dynamic Smagorinsky model on the small scale equation is the most accurate of the all models, the classic Smagorinsky model is the least accurate and the dynamic Smagorinsky model and spectral eddy viscosity model give intermediate results with small differences. The generalized sweeping hypothesis implies that time correlations in decaying isotropic turbulence are mainly determined by the instantaneous energy spectra and sweeping velocities. The analysis based on the sweeping hypothesis explains the discrepancies in our numerical simulations: the LES overpredicts the decorrelation time scales because the sweeping velocities are smaller than the DNS values, and underpredicts the magnitudes of time correlations because the energy spectrum levels are lower than the DNS ones. Since the sweeping velocity is determined by the energy spectra, one concludes that an accurate prediction of the instantaneous energy spectra guarantees the accuracy of time correlations. An analytical expression of sound power spectra based on Lighthill's theory and the quasi-normal closure assumption suggests that the sound power spectra are sensitive to errors in time correlations. Small errors in time correlations can cause significant errors in the sound power spectra, which exhibit a sizable drop at moderate to high frequencies accompanied by a shift of the peaks to lower frequencies. Based on the above analysis, two possible ways to improve the acoustic power spectrum predictions can be considered. The first is to construct better SGS models to improve the LES accuracy for time correlations. The second is to remedy the temporal statistics of the Lighthill stress tensor in order to "recover" the contribution from the unresolved scales in LES to time correlations.

He, Guo-Wei↗
NASA NTRS2011

Cycle Time Reduction in Trapped Mercury Ion Atomic Frequency Standards

The use of the mercury ion isotope (201)Hg(+) was examined for an atomic clock. Taking advantage of the faster optical pumping time in (201)Hg(+) reduces both the state preparation and the state readout times, thereby decreasing the overall cycle time of the clock and reducing the impact of medium-term LO noise on the performance of the frequency standard. The spectral overlap between the plasma discharge lamp used for (201)Hg(+) state preparation and readout is much larger than that of the lamp used for the more conventional (199)Hg(+). There has been little study of (201)Hg(+) for clock applications (in fact, all trapped ion clock work in mercury has been with (199)Hg(+); however, recently the optical pumping time in (201)Hg(+) has been measured and found to be 0.45 second, or about three times faster than in (199)Hg(+) due largely to the better spectral overlap. This can be used to reduce the overall clock cycle time by over 2 seconds, or up to a factor of 2 improvement. The use of the (201)Hg(+) for an atomic clock is totally new. Most attempts to reduce the impact of LO noise have focused on reducing the interrogation time. In the trapped ion frequency standards built so far at JPL, the optical pumping time is already at its minimum so that no enhancement can be had by shortening it. However, by using (201)Hg(+), this is no longer the case. Furthermore, integrity monitoring, the mechanism that determines whether the clock is functioning normally, cannot happen faster than the clock cycle time. Therefore, a shorter cycle time will enable quicker detection of failure modes and recovery from them.

Burt, Eric A.↗
NASA NTRS2017

Time and Measurement Days

Questions in data analysis involving the concepts of time and measurement are often pushed into the background or reserved for a philosophical discussion. Some examples are: a) Is causality a consequence of the laws of physics, or can the arrow of time be reversed? b) Can we determine the arrow of time of an event? c) Do we need the continuum hypothesis for the underlying function in any measurement process? d) Can we say anything about the analyticity of the underlying process of an event? e) Would it be valid to model a non-analytical process as function of time? f) What are the implications of all these questions for classical Fourier techniques? However, in the age of big data gathered either from space missions supplying ultra-precise long time series, or e.g. LIGO data from the ground, the moment to bring these questions to the foreground seems arrived. The limitations of our understanding of some fundamental processes is emphasized by the lack of solution for problems open for more than 2 decades, such as the non-detection of solar g-modes, or the modal identification of main sequence stellar pulsators like delta Scuti stars. Flicker noise or 1/f noise, for example, attributed in solar-like stars to granulation, is analyzed mostly only to apply noise reduction techniques, neither considering the classical problem of 1/f noise that was introduced a 100 years ago, nor taking into account ergodic or non-ergodic solutions that make inapplicable spectral analysis techniques in practice. This topic was discussed by Nicholas W. Watkins during the ITISE meeting held in Granada in 2016. There he presented preliminary results of his research on Mandelbrot's related work. We reproduce here his quotation of Mandelbrot (1999) "There is a sharp contrast between a highly anomalous ("non-white") noise that proceeds in ordinary clock time and a noise whose principal anomaly is that it is restricted to fractal time", suggesting a connection with the above proposed topics that could be phrased as the following additional questions:a) Is self-organized criticality (SOC) frequent in astrophysical phenomena? b) Could all fractals in nature be considered stochastic? c) Could we establish mathematical/physical relationships between chaotic and fractal behaviors in time series? d) Could the differences between fractals and chaos in terms of analyticity be used to understand the residuals of the fitting of stellar light curves? In this meeting we would like to approximate these problems from a holistic and multidisciplinary perspective, taking into account not only technical issues but also the deeper implications. In particular the concept of connectivity (introduced in Pascual-Granado et al. A&A, 2015) could be used to implement, within the framework of ARMA processes, an "arrow of time" (see attached document), and so studying the possible implications in the concept of time as envisaged by Watkins.

data analysis↗
NASA NTRS2022

Investigation of the Differences in Onset Times for Magnetically Conjugate Magnetometers

We have identified nearly 1,000 onsets using two pairs of hemispheric conjugate ground magnetometers where the onset is defined based on a sharp decline in the H component of the magnetic field at a ground magnetometer station. Specifically, we used the pair of stations at West Antarctica Ice Sheet Divide and Sanikiluaq, Canada; Syowa, Antarctica; and Tjörnes, Iceland. While the onset time in the southern hemisphere is identified by eye, the value of the differences in the onset time between the northern and southern hemispheres is determined using cross covariance. We observe differences in the onset time between the two hemispheres as large as several minutes, but 53% of the events show no difference in the onset time. Using statistics, we show that the largest differences in onset time are associated with the summer and winter seasons and when the IMF By value is limited between 0.5 and 2.5 nT, which is the IMF By range when the local time difference between the northern and southern hemisphere foot points is the smallest. The results indicate that ionospheric conductivity associated with solar illumination plays a role in the differences in onset time between the northern and southern hemisphere when only non-zero differences in onset time are considered. We validate these results with two other less robust methods. The median value of the differences in onset time indicates that the onsets occur ∼23 s earlier in the winter hemisphere than that in the summer hemisphere. It has been reported that the time difference between the start of the substorm in the magnetotail and the observed auroral break up (substorm auroral onset) in the ionosphere is 30 s to 2 min in the current disruption model and the near earth neutral line model, respectively. Our results may be of interest to those two models.

James M. Weygand↗
DOE OSTI2017

Estimation of Core Recovery Time with TRACE

Following the severe accident in Fukushima Daiichi Nuclear Power Plant in Japan, safety related concerns about the zirconium based claddings used in current light water reactors have increased. High corrosion resistance and low neutron absorption cross section make zirconium a suitable cladding material for normal nuclear plant operational conditions. However, for nuclear accident scenarios when safety functions of the plant are lost, Zr claddings undergo exothermic reaction with steam and exacerbate the accident conditions. The exothermic reaction increases the temperature of the cladding and accelerates core melting in severe accident progression. And the generation of explosive hydrogen gas accumulates and challenges the structural integrity. The accident tolerant fuel (ATF) concept is the fuel and cladding system that gives a longer time margin, compared to standard Zr claddings, to cope with the accident when core cooling systems are lost in a nuclear power plant. The ATF cladding should have slower oxidation kinetics and less hydrogen generation than Zr to alleviate consequences of cladding interaction with steam at high temperatures. The maximum time available to mitigate progression of severe accident is an important measure to evaluate benefits of ATF claddings. For this work, the maximum time available to mitigate progression of severe accident is selected as the time at which cladding melting reaction starts. Melting reaction is selected as an upper limit because of the degradation of coolable geometry. Maintaining the coolable geometry is a regulation criteria for both loss of coolant (LOCA) analyses and outside the LOCA analyses. Once the melting of the cladding starts, it exacerbates the deformation of the fuel assembly and hottest rods within the core cannot be cooled by water reflood and results in fuel melting progression. The aim of this work is to evaluate time margin gain of ATF cladding in a severe accident scenario of LOCA without safety injections and station blackout considering oxidation kinetics only. Also it is assumed that safety injections are recovered before melting of the claddings start and the sensitivity of safety injection actuation time on core recovery is analyzed. The core recovery time is the time when makeup water from safety injections quench the fuel rods and halt the progression of the severe accident. This time is estimated by initiating safety injections at different times and observing the hot rod behavior which represents 1/5. of the core in this work. It is considered that core is not recoverable when 20% of the core is damaged.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Interface flux recovery framework for constructing partitioned heterogeneous time-integration methods

A common approach for the development of partitioned schemes employing different time integrators on different subdomains is to lag the coupling terms in time. This can lead to accuracy issues, especially in multistage methods. Here, in this article, we present a novel framework for partitioned heterogeneous time-integration methods, which allows the coupling of arbitrary multistage and multistep methods without reducing their order of accuracy. At the core of our approach are accurate estimates of the interface flux obtained from the Schur complement of an auxiliary monolithic system. We use these estimates to construct a polynomial-in-time approximation of the interface flux over the current time coupling window. This approximation provides the interface boundary conditions necessary to decouple the subdomain problems at any point within the coupling window. In so doing our framework enables a flexible choice of time-integrators for the individual subproblems without compromising the time-accuracy at the coupled problem level. This feature is the main distinction between our framework and other approaches. To demonstrate the framework, we construct a family of partitioned heterogeneous time-integration methods, combining multistage and multistep methods, for a simplified tracer transport component of the coupled air-sea system in Earth system models. We report numerical tests evaluating accuracy and flux conservation for different pairs of time-integrators from the explicit Runge-Kutta and Adams-Moulton families.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

A non-intrusive parallel-in-time approach for simultaneous optimization with unsteady PDEs

This paper presents a non-intrusive framework for integrating existing unsteady partial differential equation (PDE) solvers into a parallel-in-time simultaneous optimization algorithm. The time-parallelization is provided by the non-intrusive software library XBraid [Parallel multigrid in time, software available at http://llnl.gov/casc/xbraid], which applies an iterative multigrid reduction technique to the time domain of existing time-marching schemes for solving unsteady PDEs. Its general user-interface has been extended in [S. Günther, N.R. Gauger, and J.B. Schroder, A non-intrusive parallel-in-time adjoint solver with the XBraid library, Comput. Vis. Sci. 19 (2018), pp. 85–95. Available at arXiv, math.OC/1705.00663] for computing adjoint sensitivities such that gradients of output quantities with respect to design changes can be computed parallel-in-time alongside with the primal PDE solution. In this paper, the primal and adjoint XBraid iterations are embedded into a simultaneous optimization framework, namely the One-shot method. In this method, design updates towards optimality are employed after each state and adjoint update such that optimality and feasibility of the design and the PDE solution are reached simultaneously. Finally, the time-parallel optimization method is validated on an advection-dominated flow control problem which shows significant speedup over a classical time-serial optimization algorithm.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Two-Level Convergence Theory for Multigrid Reduction in Time (MGRIT)

Here we develop a two-grid convergence theory for the parallel-in-time scheme known as multigrid reduction in time (MGRIT), as it is implemented in the open-source package [XBraid: Parallel Multigrid in Time, http://llnl.gov/casc/xbraid]. MGRIT is a scalable and multilevel approach to parallel-in-time simulations that nonintrusively uses existing time-stepping schemes, and in a specific two-level setting it is equivalent to the widely known parareal algorithm. The goal of this paper is twofold. First, we present a two-level MGRIT convergence analysis for linear problems where the spatial discretization matrix can be diagonalized, and then apply this analysis to our two basic model problems, the heat equation and the advection equation. One important assumption is that the coarse and fine time-grid propagators can be diagaonalized by the same set of eigenvectors, which is often the case when the same spatial discretization operator is used on the coarse and fine time grids. In many cases, the MGRIT algorithm is guaranteed to converge, and we demonstrate numerically that the theoretically predicted convergence rates are sharp in practice for our model problems. Second, we explore how the convergence of MGRIT compares to the stability of the chosen time-stepping scheme. In particular, we demonstrate that a stable time-stepping scheme does not necessarily imply convergence of MGRIT, although MGRIT with FCF-relaxation always converges for the diffusion dominated problems considered here.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Multigrid Reduction in Time for Chaotic Dynamical Systems

As CPU clock speeds have stagnated and high performance computers continue to have ever higher core counts, increased parallelism is needed to take advantage of these new architectures. Traditional serial time-marching schemes can be a significant bottleneck, as many types of simulations require large numbers of time-steps which must be computed sequentially. Parallel-in-time schemes, such as the Multigrid Reduction in Time (MGRIT) method, remedy this by parallelizing across time-steps and have shown promising results for parabolic problems. However, chaotic problems have proved more difficult, since chaotic initial value problems (IVPs) are inherently ill-conditioned. MGRIT relies on a hierarchy of successively coarser time-grids to iteratively correct the solution on the finest time-grid, but due to the nature of chaotic systems, small inaccuracies on the coarser levels can be greatly magnified and lead to poor coarse-grid corrections. Here we introduce a modified MGRIT algorithm based on an existing quadratically converging nonlinear extension to the multigrid Full Approximation Scheme (FAS), as well as a novel time-coarsening scheme. Together, these approaches better capture long-term chaotic behavior on coarse-grids and greatly improve convergence of MGRIT for chaotic IVPs. Further, we introduce a novel low-memory variant of the algorithm for solving chaotic PDEs with MGRIT which not only solves the IVP, but also provides estimates for the unstable Lyapunov vectors of the system. Finally, we provide supporting numerical results for the Lorenz system and demonstrate parallel speedup for the chaotic Kuramoto–Sivashinsky PDE over a significantly longer time-domain than in previous works.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A prototype scintillator real‐time beam monitor for ultra‐high dose rate radiotherapy

Background: FLASH Radiotherapy (RT) is an emergent cancer RT modality where an entire therapeutic dose is delivered at more than 1000 times higher dose rate than conventional RT. For clinical trials to be conducted safely, a precise and fast beam monitor that can generate out-of-tolerance beam interrupts is required. This paper describes the overall concept and provides results from a prototype ultra-fast, scintillator-based beam monitor for both proton and electron beam FLASH applications. Purpose: A FLASH Beam Scintillator Monitor (FBSM) is being developed that employs a novel proprietary scintillator material. The FBSM has capabilities that conventional RT detector technologies are unable to simultaneously provide: (1) large area coverage; (2) a low mass profile; (3) a linear response over a broad dynamic range; (4) radiation hardness; (5) real-time analysis to provide an IEC-compliant fast beam-interrupt signal based on true two-dimensional beam imaging, radiation dosimetry and excellent spatial resolution. Methods: The FBSM uses a proprietary low mass, less than 0.5 mm water equivalent, non-hygroscopic, radiation tolerant scintillator material (designated HM: hybrid material) that is viewed by high frame rate CMOS cameras. Folded optics using mirrors enable a thin monitor profile of ∼10 cm. A field programmable gate array (FPGA) data acquisition system generates real-time analysis on a time scale appropriate to the FLASH RT beam modality: 100–1000 Hz for pulsed electrons and 10–20 kHz for quasi-continuous scanning proton pencil beams. An ion beam monitor served as the initial development platform for this work and was tested in low energy heavy-ion beams ( 86 Kr +26 and protons). A prototype FBSM was fabricated and then tested in various radiation beams that included FLASH level dose per pulse electron beams, and a hospital RT clinic with electron beams. Results: Results presented in this report include image quality, response linearity, radiation hardness, spatial resolution, and real-time data processing. Furthermore, the HM scintillator was found to be highly radiation damage resistant. It exhibited a small 0.025%/kGy signal decrease from a 216 kGy cumulative dose resulting from continuous exposure for 15 min at a FLASH compatible dose rate of 237 Gy/s. Measurements of the signal amplitude versus beam fluence demonstrate linear response of the FBSM at FLASH compatible dose rates of >40 Gy/s. Comparison with commercial Gafchromic film indicates that the FBSM produces a high resolution 2D beam image and can reproduce a nearly identical beam profile, including primary beam tails. The spatial resolution was measured at 35–40 µm. Tests of the firmware beta version show successful operation at 20 000 Hz frame rate or 50 µs/frame, where the real-time analysis of the beam parameters is achieved in less than 1 µs. Conclusions: The FBSM is designed to provide real-time beam profile monitoring over a large active area without significantly degrading the beam quality. A prototype device has been staged in particle beams at currents of single particles up to FLASH level dose rates, using both continuous ion beams and pulsed electron beams. Using a novel scintillator, beam profiling has been demonstrated for currents extending from single particles to 10 nA currents. Radiation damage is minimal and even under FLASH conditions would require ≥50 kGy of accumulated exposure in a single spot to result in a 1% decrease in signal output. Beam imaging is comparable to radiochromic films, and provides immediate images without hours of processing. Real-time data processing, taking less than 50 µs (combined data transfer and analysis times), has been implemented in firmware for 20 kHz frame rates for continuous proton beams.

2D beam imaging↗
DOE OSTI2017

Measurements and interpretation of shock tube ignition delay times in highly CO 2 diluted mixtures using multiple diagnostics

Common definitions for ignition delay time are often hard to determine due to the issue of bifurcation and other non-idealities that result from high levels of CO 2 addition. Using high-speed camera imagery in comparison with more standard methods (e.g., pressure, emission, and laser absorption spectroscopy) to measure the ignition delay time, the effect of bifurcation has been examined in this study. Experiments were performed at pressures between 0.6 and 1.2 atm for temperatures between 1650 and 2040 K. The equivalence ratio for all experiments was kept at a constant value of 1 with methane as the fuel. The CO 2 mole fraction was varied between a value of X CO2 = 0.00 to 0.895. The ignition delay time was determined from three different measurements at the sidewall: broadband chemiluminescent emission captured via a photodetector, CH 4 concentrations determined using a distributed feedback interband cascade laser centered at 3403.4 nm, and pressure recorded via a dynamic Kistler type transducer. All methods for the ignition delay time were compared to high-speed camera images taken of the axial cross-section during combustion. Methane time-histories and the methane decay times were also measured using the laser. It was determined that the flame could be correlated to the ignition delay time measured at the side wall but that the flame as captured by the camera was not homogeneous as assumed in typical shock tube experiments. The bifurcation of the shock wave resulted in smaller flames with large boundary layers and that the flame could be as small as 30% of the cross-sectional area of the shock tube at the highest levels of CO 2 dilution. Here, comparisons between the camera images and the different ignition delay time methods show that care must be taken in interpreting traditional ignition delay data for experiments with large bifurcation effects as different methods in measuring the ignition delay time could result in different interpretations of kinetic mechanisms and impede the development of future mechanisms.

03 NATURAL GAS↗
DOE OSTI2020

An error-controlled adaptive time-stepping method for particle advancement in coupled CFD-DEM simulations

Coupled Computational-Fluid-Dynamics (CFD) and Discrete-Element-Method (DEM) models provide an accurate description of multiphase physical systems where a solid granular particle phase exists in an underlying gaseous continuous medium. The time integration of the granular phase in these simulations is typically handled using an explicit scheme with a constant time-step among all particles that is invariant in time to resolve inter-particle collisions. A locally third-order accurate adaptive time integration technique for particles that employs an embedded locally second-order scheme for error determination is presented in this work. The particle time-step size is dynamically adapted based on solution error, thus leading to significant savings in computational time. The efficacy of our scheme is quantified using four test cases of varying complexity (binary collision, homogeneous cooling system, fluidized bed and hopper discharge). The adaptive time-stepping method exhibits improved performance (~ 2–3 times in most of the cases studied) compared to three commonly used non-adaptive time-step methods (first-order Euler-explicit, second-order Adams-Bashforth and third-order Runge-Kutta schemes), while maintaining the same level of accuracy and parallel scalability.

42 ENGINEERING↗
DOE OSTI2014

Spin-dependent delay time in ferromagnet/insulator/ferromagnet heterostructures

We study theoretically spin-dependent group delay and dwell time in ferromagnet/insulator/ferromagnet (FM/I/FM) heterostructure. The results indicate that, when the electrons with different spin orientations tunnel through the FM/I/FM junction, the spin-up process and the spin-down process are separated on the time scales. As the self-interference delay has the spin-dependent features, the variations of spin-dependent dwell-time and spin-dependent group-delay time with the structure parameters appear different features, especially, in low incident energy range. These different features show up as that the group delay times for the spin-up electrons are always longer than those for spin-down electrons when the barrier height or incident energy increase. In contrast, the dwell times for the spin-up electrons are longer (shorter) than those for spin-down electrons when the barrier heights (the incident energy) are under a certain value. When the barrier heights (the incident energy) exceed a certain value, the dwell times for the spin-up electrons turn out to be shorter (longer) than those for spin-down electrons. In addition, the group delay time and the dwell time for spin-up and down electrons also relies on the comparative direction of magnetization in two FM layers and tends to saturation with the thickness of the barrier.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2013

TIME-DEPENDENT PHOTOIONIZATION OF GASEOUS NEBULAE: THE PURE HYDROGEN CASE

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.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

THE NANOGRAV NINE-YEAR DATA SET: OBSERVATIONS, ARRIVAL TIME MEASUREMENTS, AND ANALYSIS OF 37 MILLISECOND PULSARS

We present high-precision timing observations spanning up to nine years for 37 millisecond pulsars monitored with the Green Bank and Arecibo radio telescopes as part of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) project. We describe the observational and instrumental setups used to collect the data, and methodology applied for calculating pulse times of arrival; these include novel methods for measuring instrumental offsets and characterizing low signal-to-noise ratio timing results. The time of arrival data are fit to a physical timing model for each source, including terms that characterize time-variable dispersion measure and frequency-dependent pulse shape evolution. In conjunction with the timing model fit, we have performed a Bayesian analysis of a parameterized timing noise model for each source, and detect evidence for excess low-frequency, or “red,” timing noise in 10 of the pulsars. For 5 of these cases this is likely due to interstellar medium propagation effects rather than intrisic spin variations. Subsequent papers in this series will present further analysis of this data set aimed at detecting or limiting the presence of nanohertz-frequency gravitational wave signals.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Assessing the sequestration time scales of some ocean-based carbon dioxide reduction strategies

Abstract Ocean-based carbon dioxide (CO 2 ) removal (CDR) strategies are an important part of the portfolio of approaches needed to achieve negative greenhouse gas emissions. Many ocean-based CDR strategies rely on injecting CO 2 or organic carbon (that will eventually become CO 2 ) into the ocean interior, or enhancing the ocean’s biological pump. These approaches will not result in permanent sequestration, because ocean currents will eventually return the injected CO 2 back to the surface, where it will be brought into equilibrium with the atmosphere. Here, a model of steady state global ocean circulation and mixing is used to assess the time scales over which CO 2 injected in the ocean interior remains sequestered from the atmosphere. There will be a distribution of sequestration times for any single discharge location due to the infinite number of pathways connecting a location at depth with the sea surface. The resulting probability distribution is highly skewed with a long tail of very long transit times, making mean sequestration times much longer than typical time scales. Deeper discharge locations will sequester purposefully injected CO 2 much longer than shallower ones and median sequestration times are typically decades to centuries, and approach 1000 years in the deep North Pacific. Large differences in sequestration times occur both within and between the major ocean basins, with the Pacific and Indian basins generally having longer sequestration times than the Atlantic and Southern Oceans. Assessments made over a 50 year time horizon illustrates that most of the injected carbon will be retained for injection depths greater than 1000 m, with several geographic exceptions such as the Western North Atlantic. Ocean CDR strategies that increase upper ocean ecosystem productivity with the goal of exporting more carbon to depth will have mainly a short-term influence on atmospheric CO 2 levels because ∼70% will be transported back to the surface ocean within 50 years. The results presented here will help plan appropriate ocean CDR strategies that can help limit climate damage caused by fossil fuel CO 2 emissions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Development of Real-Time High-Density Pulsar Data Transmission and Processing for Grid Synchronization

Taking advantage of the extreme stability of the pulsar period, it can serve as the timing source for grid synchronization to compensate for the timing drift instigated by the loss of GPS signal. Nevertheless, the real-time transmission and processing of the pulsar data suffer from its high-frequency data rate, varying from megahertz to gigahertz, resulting in reduced computing speed and increased time delay. To mitigate this issue, the hardware and software frameworks are implemented for the high-density pulsar data transmission and processing for grid synchronization in this research. Initially, the high-density pulsar data is transferred using open-source software. The complementary duty cycle timing module is designed to coordinate the operation of the dual-channel high-speed interface and software. Subsequently, the multiple-threading is applied to the receiving, parsing, and splicing pulsar data. Next, the pulsar signal extraction method is implemented based on the polyphase filterbank and time of arrival estimation. Ultimately, real-time performance verification experiments are carried out for different components under two hardware platforms. Finally, the results demonstrate that only 0.482 s is required for processing 4 Gigabyte data through multiple-threading, which is 3.8 times faster than the single thread. The pulsar signal extraction can also be executed within 707 ms for 4.8 seconds of data, thereby indicating that real-time requirements can be met.

24 POWER TRANSMISSION AND DISTRIBUTION↗