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.

SourceResearch connectionFrom yearThrough year

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

At least 973 records · Page 54

All sources NASA NTRS DOE OSTI

DOE OSTI2018

Effect of intermediate annealing time on microstructure, texture and mechanical properties of Al-Mg-Si-Cu alloy

Highlights: • Long intermediate annealing time could be beneficial to the improvement of deep drawability. • Long intermediate annealing time could give rise to more particles with a larger average size. • As intermediate annealing time increases, strength will change and it has no positive correlation with annealing time, • A particle dissolution model is proposed to illustrate the variation of the strengths with intermediate annealing time. - Abstract: The effect of intermediate annealing time on the microstructure, texture and mechanical properties of Al-Mg-Si-Cu alloys was studied by mechanical properties, microstructure and texture characterization in the present study. The results show that the intermediate annealing time has some influence on the microstructure, recrystallization texture, and final mechanical properties. After the pre-aging treatment and natural aging treatment (T4P treatment), in contrast to the sheet with short intermediate annealing time, the sheet with long intermediate annealing time possesses higher plastic strain ratio r value. Strength has no positive correlation with annealing time, and when applying the annealing time of 16 h, the strength is highest. Before solution treatment, the sheet with long intermediate annealing time has more particles with a larger average size, compared with the sheet with short intermediate annealing time. After solution treatment, they possess the similar recrystallization equiaxed grain structure and recrystallization texture. The sheet with short intermediate annealing time is comprised of Cube{sub ND} {001} <310>, Goss {110}<001> and P {011}<122> orientations, while the sheet with long intermediate annealing time possesses weak texture including Cube{sub ND} {001}<310>, Goss {110}<001> and {015}<751> orientations. Finally, a particle dissolution model during solution treatment is proposed to explain the variation of strength with annealing time.

36 MATERIALS SCIENCE↗

DOE OSTI2013

THE BENEFITS OF VLBI ASTROMETRY TO PULSAR TIMING ARRAY SEARCHES FOR GRAVITATIONAL RADIATION

Precision astrometry is an integral component of successful pulsar timing campaigns. Astrometric parameters are commonly derived by fitting them as parameters of a timing model to a series of pulse times of arrival (TOAs). TOAs measured to microsecond precision over spans of several years can yield position measurements with sub-milliarcsecond precision. However, timing-based astrometry can become biased if a pulsar displays any red spin noise or a red signal produced by the stochastic gravitational wave background. We investigate how noise of different spectral types is absorbed by timing models, leading to significant estimation biases in the astrometric parameters. We find that commonly used techniques for fitting timing models in the presence of red noise (Cholesky whitening) prevent the absorption of noise into the timing model remarkably well if the time baseline of observations exceeds several years, but are inadequate for dealing with shorter pulsar data sets. Independent of timing, pulsar-optimized very long baseline interferometry (VLBI) is capable of providing position estimates precise to the sub-milliarcsecond levels needed for high-precision timing. In order to make VLBI astrometric parameters useful in pulsar timing models, the transformation between the International Celestial Reference Frame (ICRF) and the dynamical solar system ephemeris used for pulsar timing must be constrained to within a few microarcseconds. We compute a transformation between the ICRF and pulsar timing frames and quantitatively discuss how the transformation will improve in coming years. We find that incorporating VLBI astrometry into the timing models of pulsars for which only a couple of years of timing data exist will lead to more realistic assessments of red spin noise and could enhance the amplitude of gravitational wave signatures in post-fit timing residuals by factors of 20 or more.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

SU-E-T-629: Prediction of the ViewRay Radiotherapy Treatment Time for Clinical Logistics

Purpose: An algorithm is developed in our clinic, given a new treatment plan, to predict treatment delivery time for radiation therapy (RT) treatments of patients on ViewRay magnetic resonance-image guided radiation therapy (MR-IGRT) delivery system. This algorithm is necessary for managing patient treatment appointments, and is useful as an indicator to assess the treatment plan complexity. Methods: A patient’s total treatment delivery time, not including time required for localization, may be described as the sum of four components: (1) the treatment initialization time; (2) the total beam-on time; (3) the gantry rotation time; and (4) the multileaf collimator (MLC) motion time. Each of the four components is predicted separately. The total beam-on time can be calculated using both the planned beam-on time and the decay-corrected delivery dose rate. To predict the remaining components, we quantitatively analyze the patient treatment delivery record files. The initialization time is demonstrated to be random since it depends on the final gantry angle and MLC leaf positions of the previous treatment. Based on modeling the relationships between the gantry rotation angles and the corresponding rotation time, and between the furthest MLC leaf moving distance and the corresponding MLC motion time, the total delivery time is predicted using linear regression. Results: The proposed algorithm has demonstrated the feasibility of predicting the ViewRay treatment delivery time for any treatment plan of any patient. The average prediction error is 0.89 minutes or 5.34%, and the maximal prediction error is 2.09 minutes or 13.87%. Conclusion: We have developed a treatment delivery time prediction algorithm based on the analysis of previous patients’ treatment delivery records. The accuracy of our prediction is sufficient for guiding and arranging patient treatment appointments on a daily basis. The predicted delivery time could also be used as an indicator to assess the treatment plan complexity. This work was supported by a research grant from Viewray Inc.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

NASA NTRS2012

Extension of the Time-Spectral Approach to Overset Solvers for Arbitrary Motion

Forced periodic flows arise in a broad range of aerodynamic applications such as rotorcraft, turbomachinery, and flapping wing configurations. Standard practice involves solving the unsteady flow equations forward in time until the initial transient exits the domain and a statistically stationary flow is achieved. It is often required to simulate through several periods to remove the initial transient making unsteady design optimization prohibitively expensive for most realistic problems. An effort to reduce the computational cost of these calculations led to the development of the Harmonic Balance method [1, 2] which capitalizes on the periodic nature of the solution. The approach exploits the fact that forced temporally periodic flow, while varying in the time domain, is invariant in the frequency domain. Expanding the temporal variation at each spatial node into a Fourier series transforms the unsteady governing equations into a steady set of equations in integer harmonics that can be tackled with the acceleration techniques afforded to steady-state flow solvers. Other similar approaches, such as the Nonlinear Frequency Domain [3,4,5], Reduced Frequency [6] and Time-Spectral [7, 8, 9] methods, were developed shortly thereafter. Additionally, adjoint-based optimization techniques can be applied [10, 11] as well as frequency-adaptive methods [12, 13, 14] to provide even more flexibility to the method. The Fourier temporal basis functions imply spectral convergence as the number of harmonic modes, and correspondingly number of time samples, N, is increased. Some elect to solve the equations in the frequency domain directly, while others choose to transform the equations back into the time domain to simplify the process of adding this capability to existing solvers, but each harnesses the underlying steady solution in the frequency domain. These temporal projection methods will herein be collectively referred to as Time-Spectral methods. Time-Spectral methods have demonstrated marked success in reducing the computational costs associated with simulating periodic forced flows, but have yet to be fully applied to overset or Cartesian solvers for arbitrary motion with dynamic hole-cutting. Overset and Cartesian grid methodologies are versatile techniques capable of handling complex geometry configurations in practical engineering applications, and the combination of the Time-Spectral approach with this general capability potentially provides an enabling new design and analysis tool. In an arbitrary moving-body scenario for these approaches, a Lagrangian body moves through a fixed Eulerian mesh and mesh points in the Eulerian mesh interior to the solid body are removed (cut or blanked), leaving a hole in the Eulerian mesh. During the dynamic motion some gridpoints in the domain are blanked and do not have a complete set of time-samples preventing a direct implementation of the Time-Spectral method. Murman[6] demonstrated the Time-Spectral approach for a Cartesian solver with a rigid domain motion, wherein the hole cutting remains constant. Similarly, Custer et al. [15, 16] used the NASA overset OVERFLOW solver and limited the amount of relative motion to ensure static hole-cutting and interpolation. Recently, Mavriplis and Mundis[17] demonstrated a qualitative method for applying the Time-Spectral approach to an unstructured overset solver for arbitrary motion. The goal of the current work is to develop a robust and general method for handling arbitrary motion with the Time-Spectral approach within an overset or Cartesian mesh method, while still approaching the spectral convergence rate of the original Time-Spectral approach. The viscous OVERFLOW solver will be augmented with the new Time-Spectral algorithm and the capability of the method for benchmark problems in rotorcraft and turbomachinery will be demonstrated. This abstract begins with a brief synopsis of the Time-Spectral approach for overset grids and provides details of e current approach to allow for arbitrary motion. Model problem results in one and two dimensions are included to demonstrate the viability of the method and the convergence properties. Section IV briefly outlines the implementation into the OVERFLOW solver, and the abstract closes with a description of the benchmark test cases which will be included in the final paper.

Leffell, Joshua Isaac↗

DOE OSTI2014

Electronics for a Picosecond Time-of-flight Measurement

TITLE: Electronics for a Picosecond Time-of-flight Measurement ABSTRACT: Time-of-flight (TOF) detectors have historically been used as part of the particle identification capability of multi-purpose particle physics detectors. An accurate time measurement, combined with a momentum measurement based on the curvature of the track in a magnetic field, is often sufficient to determine the particle's mass, and thus its identity. Such detectors typically have measured the particle flight time extremely precisely, with an uncertainty of one hundred trillionths of a second (also referred to as 100 picoseconds). To put this in perspective it would be like counting all the people on the Earth and getting it right within 1 person! Another use of TOFs is to measure the vertex of the event, which is the location along the beam line where the incoming particles (typically protons) collide. This vertex positon is a well measured quantity for events where the protons collide “head on” as the outgoing particles produced when you blast the proton apart can be used to trace back to a vertex point from which they originated. More frequently the protons just strike a glancing blow and remain intact—in this case they are nearly parallel to the beam and you cannot tell their vertex without this ability to precisely measure the time of flight of the protons. Occasionally both happen in the same event, that is, a central system and two protons are produced. But are they from the same collision, or just a boring background where more than one collision in the same bunch crossing conspire to fake the signal of interest? That’s where the timing of the protons comes into play. The main idea is to measure the time it takes for the two protons to reach TOF detectors positioned equidistant from the center of the main detector. If the vertex is displaced to one side than that detector will measure a shorter time while the other side detector will measure a correspondingly longer time. Taking into account the speed of the particles, which is very close to the speed of light, an accuracy of 100 ps gives a vertex measurement of a few cm or a little more than an inch. At the Large Hadron Collider, where there are up to a hundred billion protons per bunch, and the collision region is compressed to a few inches, that is just not good enough. A higher level of precision is needed to determine whether the vertex of the protons and that of the central system are the same. A factor 10 improvement in the timing measurement to the 10 trillionths of a second level, for example would give the requisite 10 times improvement in the vertex measurement. An accurate measurement of the flight time depends on three key elements: the radiator that produces light when the proton passes through it, the photo-sensor that converts the light to an electrical signal, and the electronics that convert this electrical signal into a time measurement with a compact recordable format. With recent improvements in detector design featuring a series of quartz radiators connected to a micro-channel plate photomultiplier tube, this superior measurement capability is within reach if the readout electronics are sufficiently performant. As a result of the funding of this proposal, we have achieved our primary goals. 1) We have developed from scratch or improved upon existing designs of the full chain of electronics that can maintain the performance of a TOF detector from the output of the photo-sensor to the recording of a compressed data packet containing the timing information. We have accomplished this with a cost effective modular approach such that some or all of the components in the chain could easily be adapted for use in diverse particle physics experiments or in other areas where precise timing is required, such as medical and homeland security devices.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2017

Geothermal Drilling: A Baseline Study of Nonproductive Time Related to Lost Circulation

A single production geothermal well costs up to seven million dollars and a typical geothermal field contains 10 to 20 wells for a development cost of potentially $140 million according to a Geothermal Resources Group study. Capital cost coupled with exploration and resource uncertainty lead to high upfront project risk. This study analyzed the time and cost of detailed geothermal drilling operations. We found the largest cause of nonproductive time in geothermal wells is advancing through lost circulation zones, where drilling fluid is lost from the well. Mitigation measures implemented during lost circulation events can increase drilling time by a significant fraction of the drilling days. There is opportunity to reduce the drilling costs by addressing logistical inefficiencies and ineffective mitigation measures to minimize nonproductive time. This study reviewed daily drilling data for 38 geothermal wells drilled since 2009 in the U.S. The analysis presented here summarizes the lost circulation events identified in the available drilling data by assessing time per event, costs incurred, materials used, success of mitigation strategies, and resolution time. This study found that wells averaged over 100 hours of unprogrammed nonproductive time due to lost circulation, adding an estimated $185,000 to each well in rig costs. A plan for future research aimed at reducing lost circulation nonproductive time and cost is presented. operations. We found the largest cause of nonproductive time in geothermal wells is advancing through lost circulation zones, where drilling fluid is lost from the well. Mitigation measures implemented during lost circulation events can increase drilling time by a significant fraction of the drilling days. There is opportunity to reduce the drilling costs by addressing logistical inefficiencies and ineffective mitigation measures to minimize nonproductive time. This study reviewed daily drilling data for 38 geothermal wells drilled since 2009 in the U.S. The analysis presented here summarizes the lost circulation events identified in the available drilling data by assessing time per event, costs incurred, materials used, success of mitigation strategies, and resolution time. This study found that wells averaged over 100 hours of unprogrammed nonproductive time due to lost circulation, adding an estimated $185,000 to each well in rig costs. A plan for future research aimed at reducing lost circulation nonproductive time and cost is presented.

15 GEOTHERMAL ENERGY↗

DOE OSTI2014

Artifacts in digital coincidence timing

Digital methods are becoming increasingly popular for measuring time differences, and are the de facto standard in PET cameras. These methods usually include a master system clock and a (digital) arrival time estimate for each detector that is obtained by comparing the detector output signal to some reference portion of this clock (such as the rising edge). Time differences between detector signals are then obtained by subtracting the digitized estimates from a detector pair. A number of different methods can be used to generate the digitized arrival time of the detector output, such as sending a discriminator output into a time to digital converter (TDC) or digitizing the waveform and applying a more sophisticated algorithm to extract a timing estimator.All measurement methods are subject to error, and one generally wants to minimize these errors and so optimize the timing resolution. A common method for optimizing timing methods is to measure the coincidence timing resolution between two timing signals whose time difference should be constant (such as detecting gammas from positron annihilation) and selecting the method that minimizes the width of the distribution (i.e. the timing resolution). Unfortunately, a common form of error (a nonlinear transfer function) leads to artifacts that artificially narrow this resolution, which can lead to erroneous selection of the 'optimal' method. In conclusion, the purpose of this note is to demonstrate the origin of this artifact and suggest that caution should be used when optimizing time digitization systems solely on timing resolution minimization.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Eisenbud–Wigner–Smith time delay in atom–laser interactions

We report recent development of ultrafast measurement techniques on the sub-femtosecond time scale has enabled us to see the atom–laser interaction in real time. One of the important observables to study the temporally resolved dynamics is the various time delay in the process. Time delays is an experimentally measurable quantity, where the time is not. Time delays in quantum collisions and in photoionization/photodetachment of atomic and molecular systems is reviewed. The formalisms of time delay in the context of quantum collisions by Eisenbud and Wigner that by Smith and their equivalence is discussed. The time-reversal symmetry between solutions with outgoing and ingoing wave boundary conditions allow us to interpret photoionization/photodetachment as half-scattering. Subsequently, the formalism of Eisenbud–Wigner–Smith time delay is extended in photoionization/photodetachment from collisions. We discuss that the time delay can be written as a self-adjoint quantum operator which depicts its measurability. On the experimental side, the measurement protocols are reviewed for attosecond streaking and RABBITT. Along with the EWS time delay, the measurement introduces an additional (which is negative) time delay. From the measurement of total delay and an analytical/numerical estimate of the measurement induced part, one gets the EWS time delay in the process. A few illustrative examples of studies on time delay are given to get a flavour of the outstanding advances made in this field in the last two decades.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2022

Layer Time Control for Large Scale Additive Manufacturing Using High Performance Computing

This work proposes to optimize an additive manufacturing AM process to reduce energy and printing cost. The polymer AM process is inherently dependent on the time-temperature history of each layer to maintain geometric tolerances and mechanical integrity. Our preliminary study shows that regression-based layer time control model using thermal images could result in up to 30% build time reduction for simple geometries. This proposed work would use high-performance computing (HPC) to couple the data-driven model with thermal simulation for better predicting layer temperature profiles, improving throughput of large-scale additive manufacturing, and reducing its energy cost. We have developed a method to optimize a layer deposition time (a.k.a. layer time) for large-scale AM via physics-based simulations. A long layer time leads to an over-cooled surface on which a new layer is deposited, and therefore, it may result in a weak bonding or debonding between layers, cracking, or warping. A short layer time leads to a high temperature of the structure due to insufficient cooling, and therefore, the structure may not be stiff enough and may collapse during manufacturing. Therefore, it is important to estimate the optimal layer time in additive manufacturing for a high-quality product. The temperature of a top layer right before deposition is recommended to be slightly higher than the glass temperature of the material. A temperature cooling was approximated to an exponential function of time, and the optimized layer time was obtained based on a target temperature while maintaining a minimal printing time. The material used is carbon fiber-reinforced polycarbonate (CF/PC), and the large-scale deposition system used is LSAM TM from Thermwood Corporation. Three different layer time cases were used for experiments, and a series of thermal images were obtained via an infra-red (IR) camera during the entire AM processes. AM process simulations were performed using a finite element method and the temperature profiles from the simulation were in good agreements with those from experiments. The layer time optimization was performed based on the temperature profiles from the simulations. A layer temperature with the optimal layer time was confirmed as the target temperature through simulation. In addition to the development of a layer time optimization method, we have developed a numerical framework for AM simulation with element activations in sync with toolpath, based on an open source finite element framework, DEAL.II. A major portion of this work was presented at SAMPE 2022 Conference and Exhibition on May 2022, and published in Proceedings of SAMPE 2022.

42 ENGINEERING↗

NASA NTRS2008

Using LDPC Code Constraints to Aid Recovery of Symbol Timing

A method of utilizing information available in the constraints imposed by a low-density parity-check (LDPC) code has been proposed as a means of aiding the recovery of symbol timing in the reception of a binary-phase-shift-keying (BPSK) signal representing such a code in the presence of noise, timing error, and/or Doppler shift between the transmitter and the receiver. This method and the receiver architecture in which it would be implemented belong to a class of timing-recovery methods and corresponding receiver architectures characterized as pilotless in that they do not require transmission and reception of pilot signals. Acquisition and tracking of a signal of the type described above have traditionally been performed upstream of, and independently of, decoding and have typically involved utilization of a phase-locked loop (PLL). However, the LDPC decoding process, which is iterative, provides information that can be fed back to the timing-recovery receiver circuits to improve performance significantly over that attainable in the absence of such feedback. Prior methods of coupling LDPC decoding with timing recovery had focused on the use of output code words produced as the iterations progress. In contrast, in the present method, one exploits the information available from the metrics computed for the constraint nodes of an LDPC code during the decoding process. In addition, the method involves the use of a waveform model that captures, better than do the waveform models of the prior methods, distortions introduced by receiver timing errors and transmitter/ receiver motions. An LDPC code is commonly represented by use of a bipartite graph containing two sets of nodes. In the graph corresponding to an (n,k) code, the n variable nodes correspond to the code word symbols and the n-k constraint nodes represent the constraints that the code places on the variable nodes in order for them to form a valid code word. The decoding procedure involves iterative computation of values associated with these nodes. A constraint node represents a parity-check equation using a set of variable nodes as inputs. A valid decoded code word is obtained if all parity-check equations are satisfied. After each iteration, the metrics associated with each constraint node can be evaluated to determine the status of the associated parity check. Heretofore, normally, these metrics would be utilized only within the LDPC decoding process to assess whether or not variable nodes had converged to a codeword. In the present method, it is recognized that these metrics can be used to determine accuracy of the timing estimates used in acquiring the sampled data that constitute the input to the LDPC decoder. In fact, the number of constraints that are satisfied exhibits a peak near the optimal timing estimate. Coarse timing estimation (or first-stage estimation as described below) is found via a parametric search for this peak. The present method calls for a two-stage receiver architecture illustrated in the figure. The first stage would correct large time delays and frequency offsets; the second stage would track random walks and correct residual time and frequency offsets. In the first stage, constraint-node feedback from the LDPC decoder would be employed in a search algorithm in which the searches would be performed in successively narrower windows to find the correct time delay and/or frequency offset. The second stage would include a conventional first-order PLL with a decision-aided timing-error detector that would utilize, as its decision aid, decoded symbols from the LDPC decoder. The method has been tested by means of computational simulations in cases involving various timing and frequency errors. The results of the simulations ined in the ideal case of perfect timing in the receiver.

Jones, Christopher↗

NASA NTRS2009

Local-in-Time Adjoint-Based Method for Optimal Control/Design Optimization of Unsteady Compressible Flows

.We study local-in-time adjoint-based methods for minimization of ow matching functionals subject to the 2-D unsteady compressible Euler equations. The key idea of the local-in-time method is to construct a very accurate approximation of the global-in-time adjoint equations and the corresponding sensitivity derivative by using only local information available on each time subinterval. In contrast to conventional time-dependent adjoint-based optimization methods which require backward-in-time integration of the adjoint equations over the entire time interval, the local-in-time method solves local adjoint equations sequentially over each time subinterval. Since each subinterval contains relatively few time steps, the storage cost of the local-in-time method is much lower than that of the global adjoint formulation, thus making the time-dependent optimization feasible for practical applications. The paper presents a detailed comparison of the local- and global-in-time adjoint-based methods for minimization of a tracking functional governed by the Euler equations describing the ow around a circular bump. Our numerical results show that the local-in-time method converges to the same optimal solution obtained with the global counterpart, while drastically reducing the memory cost as compared to the global-in-time adjoint formulation.

Yamaleev, N. K.↗

NASA NTRS2004

Improved Timing Scheme for Spaceborne Precipitation Radar

An improved timing scheme has been conceived for operation of a scanning satellite-borne rain-measuring radar system. The scheme allows a real-time-generated solution, which is required for auto targeting. The current timing scheme used in radar satellites involves pre-computing a solution that allows the instrument to catch all transmitted pulses without transmitting and receiving at the same time. Satellite altitude requires many pulses in flight at any time, and the timing solution to prevent transmit and receive operations from colliding is usually found iteratively. The proposed satellite has a large number of scanning beams each with a different range to target and few pulses per beam. Furthermore, the satellite will be self-targeting, so the selection of which beams are used will change from sweep to sweep. The proposed timing solution guarantees no echo collisions, can be generated using simple FPGA-based hardware in real time, and can be mathematically shown to deliver the maximum number of pulses per second, given the timing constraints. The timing solution is computed every sweep, and consists of three phases: (1) a build-up phase, (2) a feedback phase, and (3) a build-down phase. Before the build-up phase can begin, the beams to be transmitted are sorted in numerical order. The numerical order of the beams is also the order from shortest range to longest range. Sorting the list guarantees no pulse collisions. The build-up phase begins by transmitting the first pulse from the first beam on the list. Transmission of this pulse starts a delay counter, which stores the beam number and the time delay to the beginning of the receive window for that beam. The timing generator waits just long enough to complete the transmit pulse plus one receive window, then sends out the second pulse. The second pulse starts a second delay counter, which stores its beam number and time delay. This process continues until an output from the first timer indicates there is less than one transmit pulse width until the start of the next receive event. This blocks future transmit pulses in the build-up phase. The feedback phase begins with the first timer paying off and starting the first receive window. When the first receive window is complete, the timing generator transmits the next beam from the list. When the second timer pays off, the second receive event is started. Following the second receive event, the timing generator will transmit the next beam on the list and start an additional timer. The timers work in a circular buffer fashion so there only need to be enough to cover the maximum number of echoes in flight.

Berkun, Andrew↗

NASA NTRS2004

Computationally-Efficient Minimum-Time Aircraft Routes in the Presence of Winds

A computationally efficient algorithm for minimizing the flight time of an aircraft in a variable wind field has been invented. The algorithm, referred to as Neighboring Optimal Wind Routing (NOWR), is based upon neighboring-optimal-control (NOC) concepts and achieves minimum-time paths by adjusting aircraft heading according to wind conditions at an arbitrary number of wind measurement points along the flight route. The NOWR algorithm may either be used in a fast-time mode to compute minimum- time routes prior to flight, or may be used in a feedback mode to adjust aircraft heading in real-time. By traveling minimum-time routes instead of direct great-circle (direct) routes, flights across the United States can save an average of about 7 minutes, and as much as one hour of flight time during periods of strong jet-stream winds. The neighboring optimal routes computed via the NOWR technique have been shown to be within 1.5 percent of the absolute minimum-time routes for flights across the continental United States. On a typical 450-MHz Sun Ultra workstation, the NOWR algorithm produces complete minimum-time routes in less than 40 milliseconds. This corresponds to a rate of 25 optimal routes per second. The closest comparable optimization technique runs approximately 10 times slower. Airlines currently use various trial-and-error search techniques to determine which of a set of commonly traveled routes will minimize flight time. These algorithms are too computationally expensive for use in real-time systems, or in systems where many optimal routes need to be computed in a short amount of time. Instead of operating in real-time, airlines will typically plan a trajectory several hours in advance using wind forecasts. If winds change significantly from forecasts, the resulting flights will no longer be minimum-time. The need for a computationally efficient wind-optimal routing algorithm is even greater in the case of new air-traffic-control automation concepts. For air-traffic-control automation, thousands of wind-optimal routes may need to be computed and checked for conflicts in just a few minutes. These factors motivated the need for a more efficient wind-optimal routing algorithm.

Jardin, Matthew R.↗

NASA NTRS2014

Contributions of Astronauts Aerobic Exercise Intensity and Time on Change in VO2peak during Spaceflight

There is considerable variability among astronauts with respect to changes in maximal aerobic capacity (VO2peak) during International Space Station (ISS) missions, ranging from a 5% increase to 30% decline. Individual differences may be due to in-flight aerobic exercise time and intensity. PURPOSE: To evaluate the effects of in-flight aerobic exercise time and intensity on change in VO2peak during ISS missions. METHODS: Astronauts (N=11) performed peak cycle tests approx 60 days before flight (L-60), on flight day (FD) approx 14, and every approx 30 days thereafter. Metabolic gas analysis and heart rate (HR) were measured continuously during the test using the portable pulmonary function system. HR and duration of each in-flight cycle ergometer and treadmill (TM) session were recorded and averaged in time segments corresponding to each peak test. Mixed effects linear regression with exercise mode (TM or cycle) as a categorical variable was used to assess the contributions of exercise intensity (%time >70% peak HR or %time >90% peak HR) and time (min/wk), adjusted for body weight, on %change in VO2peak during the mission, and incorporating the repeated-measures experimental design. RESULTS: 110 observations were included in the model (4-6 peak cycle tests per astronaut, 2 exercise devices). VO2peak was reduced from preflight throughout the mission (FD14: 13+/-13% and FD 105: 8+/-10%). Exercise intensity (%peak HR: FD14=66+/-14; FD105=75+/-8) and time (min/wk: FD14=82+/-46; FD105=158+/-40) increased during flight. The models showed main effects for exercise time and intensity with no interactions between time, intensity, and device (70% peak HR: time [z-score=2.39; P=0.017], intensity [z-score=3.51; P=0.000]; 90% peak HR: time [zscore= 3.31; P=0.001], intensity [z-score=2.24; P=0.025]). CONCLUSION: Exercise time and intensity independently contribute to %change in VO2peak during ISS missions, indicating that there are minimal values for exercise time and intensity required to maintain VO2peak. As the FD105 average exercise intensity and time did not prevent a decline in VO2peak from preflight, astronauts' exercise prescriptions should target at least 160 min of weekly aerobic exercise at an average above 75% peak HR with increased time at intensities above 90% of peak HR starting early in the mission.

Downs, Meghan E.↗

NASA NTRSDate not supplied

Changes in Oculomotor Behavior and PVT Reaction Time During One Night of Sleep Deprivation

INTRODUCTION: The Psychomotor Vigilance Test (PVT) is a widely used objective measure of sustained attention, alertness, and fatigue. Previous studies consistently show that prolonged wakefulness and disrupted sleep patterns lead to increased lapses of attention and slower reaction times. These findings have critical implications for various professions, including healthcare, aviation, transportation safety, and spaceflight. In this study, we employed linear mixed models (LMM) to investigate hourly-measured 5-minute PVT reaction times during one night of total sleep deprivation. Our goal was to characterize the contribution of homeostatic sleep pressure (time awake) and circadian phase (salivary melatonin levels) on the dynamics of PVT reaction time. MATERIALS & METHODS: Data from twelve human participants were used in the analysis. Participants were healthy non-smokers, aged 18 and 40, with normal sleep habits defined as Pittsburg Sleep Quality Index scores < 5, and Morningness – Eveningness Questionnaire scores > 42 and < 58, respectively. Participants followed a constant-routine protocol and arrived at the NASA sleep laboratory approximately 1 to 2 hours after their habitual waking time. Tympanic temperature was measured every 30 min to provide a real-time estimate of circadian phase. Participants were required to stay awake throughout the laboratory experiment, which continued until their tympanic temperature had returned to baseline levels. This approach ensured that we could capture each participant's circadian trough and recovery, typically occurring between 24 and 26 hours after waking. The PVT and saliva melatonin were measured every hour, starting approximately from 3 hours after awakening. More information about the measurement protocol can be found in the publication by Stone et al. We used PVT reaction times between 100 and 10000 ms to calculate Reciprocal Reaction Times (RRTs). LMM models were employed to assess the contribution of time awake and saliva melatonin levels on changes in RRTs. Two models were created: the first model treated time awake as a fixed effect and the second model considered both the salivary melatonin level and time awake as fixed effects. Participants were included as a random effect in both models. We compared the models using an Analysis of Variance (ANOVA). The models were computed using the 'lmer' function from the 'lme4' package in R, and LMM p-values were determined using Satterthwaite's method. RESULTS: We found that time awake had a statistically significant effect on PVT RRTs during one night of sleep deprivation [χ2 = 2132.4, df = 1, p < 0.001]. In addition, the model including both saliva melatonin levels and time awake was statistically significantly better than one that uses just time awake [χ2 = 17.4, df = 1, p < 0.001]. DISCUSSION & CONCLUSION: As expected, these preliminary results suggest that the time awake and circadian phase each contribute to performance reductions as measured by the PVT. Future analyses will explore the contributions of other factors to the model.

sleep deprivation↗

DOE OSTI2019

Time reversal symmetry of generalized quantum measurements with past and future boundary conditions

Here we expand the time reversal symmetry arguments of quantum mechanics, originally proposed by Wigner in the context of unitary dynamics, to contain situations including generalized measurements for monitored quantum systems. We propose a scheme to derive the time-reversed measurement operators by considering the Schrödinger picture dynamics of a qubit coupled to a measuring device, and show that the time-reversed measurement operators form a Positive Operator Valued Measure (POVM) set. We present three particular examples to illustrate time reversal of measurement operators: (1) the Gaussian spin measurement, (2) a dichotomous POVM for spin, and (3) the measurement of qubit fluorescence. We then propose a general rule to unravel any rank two qubit measurement, and show that the backward dynamics obeys the retrodicted equations of the forward dynamics starting from the time reversed final state. We demonstrate the time reversal invariance of dynamical equations using the example of qubit fluorescence. We also generalize the discussion of a statistical arrow of time for continuous quantum measurements introduced by Dressel et al. (Phys Rev Lett 119:220507, 2017): we show that the backward probabilities can be computed from a process similar to retrodiction from the time-reversed final state, and extend the definition of an arrow of time to ensembles prepared with pre- and post-selections, where we obtain a non-vanishing arrow of time in general. We discuss sufficient conditions for when time’s arrow vanishes and show our method also captures the contributions to time’s arrow due to natural physical processes like relaxation of an atom to its ground state. As a special case, we recover the time reversibility of the weak value as its complex conjugate using our method, and discuss how our conclusions differ from the time-symmetry argument of Aharonov–Bergmann–Lebowitz rule.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2015

Monte Carlo calculations of PET coincidence timing: single and double-ended readout

We present Monte Carlo computational methods for estimating the coincidence resolving time (CRT) of scintillator detector pairs in positron emission tomography (PET) and present results for Lu 2 SiO 5 : Ce (LSO), LaBr 3 : Ce, and a hypothetical ultra-fast scintillator with a 1 ns decay time. The calculations were applied to both single-ended and double-ended photodetector readout with constant-fraction triggering. They explicitly include (1) the intrinsic scintillator properties (luminosity, rise time, decay time, and index of refraction), (2) the exponentially distributed depths of interaction, (3) the optical photon transport efficiency, delay, and time dispersion, (4) the photodetector properties (fill factor, quantum efficiency, transit time jitter, and single electron response), and (5) the determination of the constant fraction trigger level that minimizes the CRT. The calculations for single-ended readout include the delayed photons from the opposite reflective surface. The calculations for double-ended readout include (1) the simple average of the two photodetector trigger times, (2) more accurate estimators of the annihilation photon entrance time using the pulse height ratio to estimate the depth of interaction and correct for annihilation photon, optical photon, and trigger delays, and (3) the statistical lower bound for interactions at the center of the crystal. For time-of-flight (TOF) PET we combine stopping power and TOF information in a figure of merit equal to the sensitivity gain relative to whole-body non-TOF PET using LSO. For LSO crystals 3 mm × 3 mm × 30 mm, a decay time of 37 ns, a total photoelectron count of 4000, and a photodetector with 0.2 ns full-width at half-maximum (fwhm) timing jitter, single-ended readout has a CRT of 0.16 ns fwhm and double-ended readout has a CRT of 0.111 ns fwhm. For LaBr 3 : Ce crystals 3 mm × 3 mm × 30 mm, a rise time of 0.2 ns, a decay time of 18 ns, and a total of 7600 photoelectrons the CRT numbers are 0.14 ns and 0.072 ns fwhm, respectively. For a hypothetical ultra-fast scintillator 3 mm × 3 mm × 30 mm, a decay time of 1 ns, and a total of 4000 photoelectrons, the CRT numbers are 0.070 and 0.020 ns fwhm, respectively. Over a range of examples, values for double-ended readout are about 10% larger than the statistical lower bound.

Monte Carlo↗

DOE OSTI2014

Dependence of Brownian and Néel relaxation times on magnetic field strength

Purpose: In magnetic particle imaging (MPI) and magnetic particle spectroscopy (MPS) the relaxation time of the magnetization in response to externally applied magnetic fields is determined by the Brownian and Néel relaxation mechanisms. Here the authors investigate the dependence of the relaxation times on the magnetic field strength and the implications for MPI and MPS. Methods: The Fokker–Planck equation with Brownian relaxation and the Fokker–Planck equation with Néel relaxation are solved numerically for a time-varying externally applied magnetic field, including a step-function, a sinusoidally varying, and a linearly ramped magnetic field. For magnetic fields that are applied as a step function, an eigenvalue approach is used to directly calculate both the Brownian and Néel relaxation times for a range of magnetic field strengths. For Néel relaxation, the eigenvalue calculations are compared to Brown's high-barrier approximation formula. Results: The relaxation times due to the Brownian or Néel mechanisms depend on the magnitude of the applied magnetic field. In particular, the Néel relaxation time is sensitive to the magnetic field strength, and varies by many orders of magnitude for nanoparticle properties and magnetic field strengths relevant for MPI and MPS. Therefore, the well-known zero-field relaxation times underestimate the actual relaxation times and, in particular, can underestimate the Néel relaxation time by many orders of magnitude. When only Néel relaxation is present—if the particles are embedded in a solid for instance—the authors found that there can be a strong magnetization response to a sinusoidal driving field, even if the period is much less than the zero-field relaxation time. For a ferrofluid in which both Brownian and Néel relaxation are present, only one relaxation mechanism may dominate depending on the magnetic field strength, the driving frequency (or ramp time), and the phase of the magnetization relative to the applied magnetic field. Conclusions: A simple treatment of Néel relaxation using the common zero-field relaxation time overestimates the relaxation time of the magnetization in situations relevant for MPI and MPS. For sinusoidally driven (or ramped) systems, whether or not a particular relaxation mechanism dominates or is even relevant depends on the magnetic field strength, the frequency (or ramp time), and the phase of the magnetization relative to the applied magnetic field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗