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,243 records · Page 69

Applying queueing theory to evaluate wait-time-savings of triage algorithms

Abstract In the past decade, artificial intelligence (AI) algorithms have made promising impacts in many areas of healthcare. One application is AI-enabled prioritization software known as computer-aided triage and notification (CADt). This type of software as a medical device is intended to prioritize reviews of radiological images with time-sensitive findings, thus shortening the waiting time for patients with these findings. While many CADt devices have been deployed into clinical workflows and have been shown to improve patient treatment and clinical outcomes, quantitative methods to evaluate the wait-time-savings from their deployment are not yet available. In this paper, we apply queueing theory methods to evaluate the wait-time-savings of a CADt by calculating the average waiting time per patient image without and with a CADt device being deployed. We study two workflow models with one or multiple radiologists (servers) for a range of AI diagnostic performances, radiologist’s reading rates, and patient image (customer) arrival rates. To evaluate the time-saving performance of a CADt, we use the difference in the mean waiting time between the diseased patient images in the with-CADt scenario and that in the without-CADt scenario as our performance metric. As part of this effort, we have developed and also share a software tool to simulate the radiology workflow around medical image interpretation, to verify theoretical results, and to provide confidence intervals for the performance metric we defined. We show quantitatively that a CADt triage device is more effective in a busy, short-staffed reading setting, which is consistent with our clinical intuition and simulation results. Although this work is motivated by the need for evaluating CADt devices, the evaluation methodology presented in this paper can be applied to assess the time-saving performance of other types of algorithms that prioritize a subset of customers based on binary outputs.

Thompson, Yee Lam Elim (ORCID:0000000196537707)↗

Particle detection and non-detection in a quantum time of arrival measurement

The standard time-of-arrival distribution cannot reproduce both the temporal and the spatial profile of the modulus squared of the time-evolved wave function for an arbitrary initial state. In particular, the time-of-arrival distribution gives a non-vanishing probability even if the wave function is zero at a given point for all values of time. This poses a problem in the standard formulation of quantum mechanics where one quantizes a classical observable and uses its spectral resolution to calculate the corresponding distribution. In this work, we show that the modulus squared of the time-evolved wave function is in fact contained in one of the degenerate eigenfunctions of the quantized time-of-arrival operator. This generalizes our understanding of quantum arrival phenomenon where particle detection is not a necessary requirement, thereby providing a direct link between time-of-arrival quantization and the outcomes of the two-slit experiment. -- Highlights: •The time-evolved position density is contained in the standard TOA distribution. •Particle may quantum mechanically arrive at a given point without being detected. •The eigenstates of the standard TOA operator are linked to the two-slit experiment.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Quantitative Analysis of the Relationship Between Radiation Therapy Use and Travel Time

Purpose: To model and quantify the relationship between radiation therapy (RT) use and travel time to RT services. Methods and Materials: Population-based registries and databases were used to identify both incident cancer patient and patients receiving RT within 1 year of diagnosis (RT1y) in British Columbia, Canada, between 1992 and 2011. The effects of age, gender, diagnosis year, income, prevailing wait time, and travel duration for RT on RT1y were assessed. Significant factors from univariate analyses were included in a multivariable logistic regression model. The shape of the travel time–RT1y curve was represented by generalized additive and segmented regression models. Analyses were conducted for breast, lung, and genitourinary cancer separately and for all cancer sites combined. Results: After adjustment for age, gender, diagnosis year, income, and prevailing wait times, increasing travel time to the closest RT facility had a negative impact RT1y. The shape of the travel time–RT1y curve varied with cancer type. For breast cancer, the odds of RT1y were constant for the first 2 driving hours and decreased at 17% per hour thereafter. For lung cancer, the odds of RT1y decreased by 16% after 20 minutes and then decreased at 6% per hour. Genitourinary cancer RT1y was relatively independent of travel time. For all cancer sites combined, the odds of RT1y were constant within the first 2 driving hours and decreased at 7% per hour thereafter. Conclusions: Travel time to receive RT has a different impact on RT1y for different tumor sites. The results provide evidence-based insights for the configuration of catchment areas for new and existing cancer centers providing RT.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

An adaptive time-stepping strategy for solving the phase field crystal model

In this work, we will propose an adaptive time step method for simulating the dynamics of the phase field crystal (PFC) model. The numerical simulation of the PFC model needs long time to reach steady state, and then large time-stepping method is necessary. Unconditionally energy stable schemes are used to solve the PFC model. The time steps are adaptively determined based on the time derivative of the corresponding energy. It is found that the use of the proposed time step adaptivity cannot only resolve the steady state solution, but also the dynamical development of the solution efficiently and accurately. The numerical experiments demonstrate that the CPU time is significantly saved for long time simulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Long-time uncertainty propagation using generalized polynomial chaos and flow map composition

We present an efficient and accurate method for long-time uncertainty propagation in dynamical systems. Uncertain initial conditions and parameters are both addressed. The method approximates the intermediate short-time flow maps by spectral polynomial bases, as in the generalized polynomial chaos (gPC) method, and uses flow map composition to construct the long-time flow map. In contrast to the gPC method, this approach has spectral error convergence for both short and long integration times. The short-time flow map is characterized by small stretching and folding of the associated trajectories and hence can be well represented by a relatively low-degree basis. The composition of these low-degree polynomial bases then accurately describes the uncertainty behavior for long integration times. The key to the method is that the degree of the resulting polynomial approximation increases exponentially in the number of time intervals, while the number of polynomial coefficients either remains constant (for an autonomous system) or increases linearly in the number of time intervals (for a non-autonomous system). The findings are illustrated on several numerical examples including a nonlinear ordinary differential equation (ODE) with an uncertain initial condition, a linear ODE with an uncertain model parameter, and a two-dimensional, non-autonomous double gyre flow.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The effects of low- Z shielding on uranium isotope discrimination using the time-emission profiles of long-lived delayed neutrons

Characteristic radiation signals can often be obscured or distorted by relatively small amounts of shielding, as is the case for highly-enriched uranium (HEU). Methods that can detect and identify special nuclear material (SNM) even in the presence of shielding are therefore of high value in many nuclear security and nonproliferation applications. It has been previously shown that the buildup and decay time profiles of long-lived delayed neutron emission from induced fission are unique and can be used to successfully discriminate between uranium isotopes and infer enrichment without the need for a calibration standard. Because the long-lived delayed neutron precursors exhibit decay times on the order of tens of seconds, their time profiles are not sensitive to the delays associated with scattering or diffusion in shielding materials, which occur on much shorter time scales. The use of capture-based neutron detectors can also mitigate the effects of shielding, because unlike recoil-based detectors, they do not have a low-energy detection threshold, so that neutrons that lose their energy in the shielding are still detectable. In this work, the effects of low-Z shielding on the discrimination of uranium isotopes using delayed neutron time profile measurements are investigated for the first time. Furthermore, the delayed neutron buildup and decay time profiles of polyethylene-shielded HEU and depleted uranium objects are measured using two different capture-gated detector designs. The recorded profiles show good agreement with the shapes predicted from six-group delayed neutron models even with the shielding in place, and the enrichment of the HEU and DU objects was estimated accurately to within 5%. While this method shows potential for discriminating and characterizing extreme cases of enrichment when a moderate amount of shielding is present, simulation results suggest that significant limitations exist for intermediate enrichment levels due to distortions to the delayed neutron time profile caused by neutron reflection in the shielding.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Simulating Electron Dynamics of Complex Molecules with Time-Dependent Complete Active Space Configuration Interaction

In this paper, time-dependent electronic structure methods are growing in popularity as tools for modeling ultrafast and/or nonlinear processes, for computing spectra, and as the electronic structure component of mean-field molecular dynamics simulations. Time-dependent configuration interaction (TD-CI) offers several advantages over the widely used real-time time-dependent density functional theory: namely, that it correctly models Rabi oscillations; it offers a spin-pure description of open-shell systems; and a hierarchy of TD-CI methods can be defined that systematically approach the exact solution of the time-dependent Schrodinger equation (TDSE). In this work, we present a novel TD-CI approach that extends TD-CI to large complete active-space configuration expansions. Such extension is enabled by use of a direct configuration interaction approach that eliminates the need to explicitly build, store, or diagonalize the Hamiltonian matrix. Graphics processing unit (GPU) acceleration enables fast solution of the TDSE even for large active spaces—up to 12 electrons in 12 orbitals (853776 determinants) in this work. A symplectic split operator propagator yields long-time norm conservation. We demonstrate the applicability of our approach by computing the response of a large molecule with a strongly correlated ground state, decacene (C 42 H 24 ), to various pulses (δ-function, transform limited, chirped). Our simulations predict that chirped pulses can be used to induce dipole-forbidden transitions. Simulations of decacene using the 6-31G(d) basis set and a 12 electrons/12 orbitals active space took 20.1 h to propagate for 100 fs with a 1 attosecond time step on a single NVIDIA K40 GPU. Convergence with respect to time step is found to depend on the property being computed and the chosen active space.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Measuring Activation and Luminescence Time Scales of Upconverting NaYF 4 :Yb,Er Nanocrystals

Accurate determination of upconversion and luminescence lifetimes requires kinetic modeling of the complete time-resolved response for upconversion luminescence of NaYF4:Yb,Er nanocrystals. Prior investigations typically perform exponential fitting of the tail in the time profile or employ complex systems of differential equations to extract lifetimes. To simplify analysis while fitting the entire time-resolved response, this work introduces a set of simplified models that model the response as a convolution of upconversion and fluorescence processes. Models for two- and three-photon upconversion processes are developed and tested for NaYF4:Yb,Er nanocrystals excited by a 980 nm laser. The results are presented for the transitions 2 H 11/2 , 4 S 3/2 → 4 I 15/2 (530 nm, green emission), 4 F 9/2 → 4 I 15/2 (650 nm, red emission), and 2 H 9/2 → 4 I 15/2 (400 nm, blue emission). Even with the same number of fitting parameters, the two- and three-photon models resulted in better fitting than the simpler models of single and two convolved exponentials for 30 nm, 400 nm, and silica-coated 30 nm nanoparticles. We provide evidence that the longest time scales (0.5–4 ms) are due to the luminescence of the final state, that the energy transfer waiting times leading to upconversion are shorter (<0.1 ms), and that more than two waiting times are required to explain the entire time response. The long-time response deviates from single-exponential behavior more strongly in uncoated and silica-coated 30 nm particles than for 400 nm particles. This deviation may result from heterogeneity of the local environment of ions within the 30 nm particles.

36 MATERIALS SCIENCE↗

Dynamics from noisy data with extreme timing uncertainty

Imperfect knowledge of the times at which ‘snapshots’ of a system are recorded degrades our ability to recover dynamical information, and can scramble the sequence of events. In X-ray free-electron lasers, for example, the uncertainty—the so-called timing jitter—between the arrival of an optical trigger (‘pump’) pulse and a probing X-ray pulse can exceed the length of the X-ray pulse by up to two orders of magnitude, marring the otherwise precise time-resolution capabilities of this class of instruments. The widespread notion that little dynamical information is available on timescales shorter than the timing uncertainty has led to various hardware schemes to reduce timing uncertainty. These schemes are expensive, tend to be specific to one experimental approach and cannot be used when the record was created under ill-defined or uncontrolled conditions such as during geological events. Here we present a data-analytical approach, based on singular-value decomposition and nonlinear Laplacian spectral analysis, that can recover the history and dynamics of a system from a dense collection of noisy snapshots spanning a sufficiently large multiple of the timing uncertainty. In this work, the power of the algorithm is demonstrated by extracting the underlying dynamics on the few-femtosecond timescale from noisy experimental X-ray free-electron laser data recorded with 300-femtosecond timing uncertainty. Using a noisy dataset from a pump-probe experiment on the Coulomb explosion of nitrogen molecules, our analysis reveals vibrational wave-packets consisting of components with periods as short as 15 femtoseconds, as well as more rapid changes, which have yet to be fully explored. Our approach can potentially be applied whenever dynamical or historical information is tainted by timing uncertainty.

47 OTHER INSTRUMENTATION↗

Classical discrete time crystals

The spontaneous breaking of time-translation symmetry in periodically driven quantum systems leads to a new phase of matter: the discrete time crystal (DTC). This phase exhibits collective subharmonic oscillations that depend upon an interplay of non-equilibrium driving, many-body interactions and the breakdown of ergodicity. However, subharmonic responses are also a well-known feature of classical dynamical systems ranging from predator–prey models to Faraday waves and a.c.-driven charge density waves. This raises the question of whether these classical phenomena display the same rigidity characteristic of a quantum DTC. In this work, we explore this question in the context of periodically driven Hamiltonian dynamics coupled to a finite-temperature bath, which provides both friction and, crucially, noise. Focusing on one-dimensional chains, where in equilibrium any transition would be forbidden at finite temperature, we provide evidence that the combination of noise and interactions drives a sharp, first-order dynamical phase transition between a discrete time-translation invariant phase and an activated classical discrete time crystal (CDTC) in which time-translation symmetry is broken out to exponentially long timescales. Power-law correlations are present along a first-order line, which terminates at a critical point. We analyse the transition by mapping it to the locked-to-sliding transition of a d.c.-driven charge density wave. To conclude, building upon results from the field of probabilistic cellular automata, we conjecture the existence of classical time crystals with true long-range order, where time-translation symmetry is broken out to infinite times.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

TDCOSMO VIII. A key test of systematics in the hierarchical method of time-delay cosmography

The largest source of systematic errors in the time-delay cosmography method likely arises from the lens model mass distribution, where an inaccurate choice of model could in principle bias the value of $H$ 0 . A Bayesian hierarchical framework has been proposed which combines lens systems with kinematic data, constraining the mass profile shape at a population level. The framework has been previously validated using a small sample of lensing galaxies drawn from hydro-simulations. The goal of this work is to expand the validation to a more general set of lenses consistent with observed systems, as well as confirm the capacity of the method to combine two lens populations: one which has time delay information and one which lacks time delays and has systematically different image radii. For this purpose, we generated samples of analytic lens mass distributions made of baryons+dark matter and fit the subsequent mock images with standard power-law models. Corresponding kinematics data were also emulated. The hierarchical framework applied to an ensemble of time-delay lenses allowed us to correct the $H$ 0 bias associated with model choice to find $H$ 0 within 1.5σ of the fiducial value. We then combined this set with a sample of corresponding lens systems which have no time delays and have a source at lower $z$, resulting in a systematically smaller image radius relative to their effective radius. The hierarchical framework has successfully accounted for this effect, recovering a value of $H$ 0 which is both more precise (σ ~ 2%) and more accurate (0.7% median offset) than the time-delay set alone. This result confirms that non-time-delay lenses can nonetheless contribute valuable constraining power to the determination of $H$ 0 via their kinematic constraints, assuming they come from the same global population as the time-delay set.

79 ASTRONOMY AND ASTROPHYSICS↗

Automated Signal Timing Plan Reconstruction Using High-Resolution Event-Based Controller Data for Digital Twins

Transportation digital twins are essential tools for evaluating emerging technologies such as connected and automated vehicles, adaptive traffic signal control, and mobility optimization strategies. Realistic digital twins require accurate emulation of real-world signal controllers and detailed signal timing plans. However, signal timing plans are often unavailable or difficult to access, forcing researchers and modelers to rely on assumed fixed timings or halt their analysis. To overcome this challenge, we present a method that directly estimates signal timing plan parameters using high-resolution, event-based data from traffic signal controllers. The proposed method extracts key parameters, including cycle length, offset, phase sequence, coordinated phases, phase-specific minimum and maximum green durations, vehicle extensions, and splits under coordination. A rule-based deterministic signal timing reconstruction algorithm based on traffic signal operation rules, such as those outlined in the Signal Timing Manual, is developed and validated. We evaluate this method, which uses high-resolution controller event logs and verified signal timing plans, on 94 signalized intersections in Nashville, Tennessee, demonstrating their ability to generate accurate, simulation-ready signal timing plans for tools such as SUMO and Vissim.

Saroj, Abhilasha [ORNL] (ORCID:0000000191178063)↗

Seismic energy data analysis of Merapi volcano to test the eruption time prediction using materials failure forecast method (FFM)

The test of eruption time prediction is an effort to prepare volcanic disaster mitigation, especially in the volcano’s inhabited slope area, such as Merapi Volcano. The test can be conducted by observing the increase of volcanic activity, such as seismicity degree, deformation and SO2 gas emission. One of methods that can be used to predict the time of eruption is Materials Failure Forecast Method (FFM). Materials Failure Forecast Method (FFM) is a predictive method to determine the time of volcanic eruption which was introduced by Voight (1988). This method requires an increase in the rate of change, or acceleration of the observed volcanic activity parameters. The parameter used in this study is the seismic energy value of Merapi Volcano from 1990 – 2012. The data was plotted in form of graphs of seismic energy rate inverse versus time with FFM graphical technique approach uses simple linear regression. The data quality control used to increase the time precision employs the data correlation coefficient value of the seismic energy rate inverse versus time. From the results of graph analysis, the precision of prediction time toward the real time of eruption vary between −2.86 up to 5.49 days.

42 ENGINEERING↗

Influence of Dupree diffusivity on the occurrence scattering time advance in turbulent plasmas

The influence of Dupree diffusivity on the occurrence scattering time advance for the electron-ion collision is investigated in turbulent plasmas. The second-order eikonal method and the effective Dupree potential term associated with the plasma turbulence are employed to obtain the occurrence scattering time as a function of the diffusion coefficient, impact parameter, collision energy, thermal energy, and Debye length. The result shows that the occurrence scattering time advance decreases with an increase of the Dupree diffusivity. Hence, we have found that the influence of plasma turbulence diminishes the occurrence time advance in forward electron-ion collisions in thermal turbulent plasmas. The occurrence time advance shows that the propensity of the occurrence time advance increases with increasing scattering angle. It is also found that the effect of turbulence due to the Dupree diffusivity on the occurrence scattering time advance decreases with an increase of the thermal energy. In addition, the variation of the plasma turbulence on the occurrence scattering time advance due to the plasma parameters is also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Forward and backward in time dispersion of fluid and inertial particles in isotropic turbulence

In this paper, we investigate both theoretically and numerically the Forward-In-Time (FIT) and Backward-In-Time (BIT) dispersion of fluid and inertial particle-pairs in isotropic turbulence. Fluid particles are known to separate faster BIT than FIT in three-dimensional turbulence, and we find that inertial particles do the same. However, we find that the irreversibility in the inertial particle dispersion is in general much stronger than that for fluid particles. For example, the ratio of the BIT to FIT mean-square separation can be up to an order of magnitude larger for the inertial particles than for the fluid particles. We also find that for both the inertial and fluid particles, the irreversibility becomes stronger as the scale of their separation decreases. Regarding the physical mechanism for the irreversibility, we argue that whereas the irreversibility of fluid particle-pair dispersion can be understood in terms of a directional bias arising from the energy transfer process in turbulence, inertial particles experience an additional source of irreversibility arising from the non-local contribution to their velocity dynamics, a contribution that vanishes in the limit St → 0, where St is the particle Stokes number. For each given initial (final, in the BIT case) separation, r 0 , there is an optimum value of St for which the dispersion irreversibility is strongest, as such particles are optimally affected by both sources of irreversibility. We derive analytical expressions for the BIT, mean-square separation of inertial particles and compare the predictions with numerical data obtained from a Re λ ≈ 582 (where Re λ is the Taylor Reynolds number) Direct Numerical Simulation (DNS) of particle-laden isotropic turbulent flow. The small-time theory, which in the dissipation range is valid for times ≤max[Stτ η , τ η ] (where τ η is the Kolmogorov time scale), is in excellent agreement with the DNS. The theory for long-times is in good agreement with the DNS provided that St is small enough so that the inertial particle motion at long-times may be considered as a perturbation about the fluid particle motion, a condition that would in fact be satisfied for arbitrary St at sufficiently long-times in the limit Re λ → ∞.

42 ENGINEERING↗

Path-accelerated stochastic molecular dynamics: Parallel-in-time integration using path integrals

Massively parallel computer architectures create new opportunities for the performance of long-time scale molecular dynamics (MD) simulations. Here, we introduce the path-accelerated molecular dynamics method that takes advantage of distributed computing to reduce the wall-clock time of MD simulation via parallelization with respect to stochastic MD time steps. The marginal distribution for the time evolution of a system is expressed in terms of a path integral, enabling the use of path sampling techniques to numerically integrate MD trajectories. By parallelizing the evaluation of the path action with respect to time and by initializing the path configurations from a nonequilibrium distribution, the algorithm enables significant speedups in terms of the length of MD trajectories that can be integrated in a given amount of wall-clock time. The method is demonstrated for Brownian dynamics, although it is generalizable to other stochastic equations of motion including open systems. We apply the method to two simple systems, a harmonic oscillator and a Lennard-Jones liquid, and we show that in comparison to the conventional Euler integration scheme for Brownian dynamics, the new method can reduce the wall-clock time for integrating trajectories of a given length by more than three orders of magnitude in the former system and more than two in the latter. This new method for parallelizing MD in the dimension of time can be trivially combined with algorithms for parallelizing the MD force evaluation to achieve further speedup.

Rosa-Raíces, Jorge L.↗

Time dependent model for non-inductive ECH X-I current ramp-up for SHPD tokamak facility

Non-inductive (NI) plasma current start-up and ramp-up is an important research topic for spherical tokamak (ST) based reactors and fusion pilot plant (FPP). For a compact FPP, the OH flux availability is highly restricted due to its compact geometry. Efficient fundamental extraordinary mode (X-I) electron cyclotron heating (ECH) current start-up and ramp-up regime was identified for a reactor-like high toroidal magnetic field range which has more than a hundred times higher current drive efficiency compared to more conventional ECH methods for the relevant start-up temperature range. High current drive efficiency is possible due to the strong X-I fundamental ECH interaction only with unidirectional passing electrons constrained by the wave accessibility conditions. Here, we extend the X-I electron cyclotron current drive (ECCD) investigation to a time dependent model to simulate the non-inductive current ramp-up to 10 MA for the Sustained High-Power Density (SHPD) facility. As the X-I ECCD driven current I EC rises, due to the back EMF driven negative current, the net plasma current Ip rises more slowly with the current resistive time scale. For tokamak confinement time (both L-mode and H-Mode) which tends to rise with I p , a positive feedback results and even with constant applied ECH power, T e0 , I EC , and I p can continue to rise to very high values. However, in a realistic situation, the T e0 rise should saturate due to a number of factors such as enhanced core radiation and increased power loss at high temperature. To simulate this effect, we adopt a maximum T e0 model which would limit the temperature rise to certain T e0 . With this model, we investigated the current ramp-up for various ECH power levels and the maximum T e0 of 15, 20, and 25 keV. We find that while the power required is reduced with increasing T e0 limit due to increased current drive efficiency, the time to reach 10 MA tends to go up due to the reduced plasma resistivity for the higher T e0 limit. We also find for a given Te0 limit, the time to reach 10 MA tends to be reduced by increasing the applied ECH power by over driving the current ramp-up where I EC is driven at significantly higher level than 10 MA. While the current ramp-up time may not be an issue for the steady-state reactor systems, if it is desirable to minimize the current ramp-up time, it is prudent to have a sufficient ECH power for current over-drive and have some T e0 limiting tools such as impurity seeding for enhanced radiation. In conclusion, a well-controlled NI ECH start-up also has a potential of improving the tokamak start-up reliability and avoid run-away electrons while the NI off-axis current drive could enhance MHD stability and plasma performance improvements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Kepler eclipsing binary stars. IV. Precise eclipse times for close binaries and identification of candidate three-body systems

We present a catalog of precise eclipse times and analysis of third-body signals among 1279 close binaries in the latest Kepler Eclipsing Binary Catalog. For these short-period binaries, Kepler's 30 minute exposure time causes significant smearing of light curves. In addition, common astrophysical phenomena such as chromospheric activity, as well as imperfections in the light curve detrending process, can create systematic artifacts that may produce fictitious signals in the eclipse timings. We present a method to measure precise eclipse times in the presence of distorted light curves, such as in contact and near-contact binaries which exhibit continuously changing light levels in and out of eclipse. We identify 236 systems for which we find a timing variation signal compatible with the presence of a third body. These are modeled for the light travel time effect and the basic properties of the third body are derived. This study complements J. A. Orosz et al. (in preparation), which focuses on eclipse timing variations of longer period binaries with flat out-of-eclipse regions. Together, these two papers provide comprehensive eclipse timings for all binaries in the Kepler Eclipsing Binary Catalog, as an ongoing resource freely accessible online to the community.

79 ASTRONOMY AND ASTROPHYSICS↗