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,171 records · Page 65

First Cloud-to-Ground Lightning Timing Study

NASA's LSP, GSDO and other programs use the probability of cloud-to-ground (CG) lightning occurrence issued by the 45th Weather Squadron (45 WS) in their daily and weekly lightning probability forecasts. These organizations use this information when planning potentially hazardous outdoor activities, such as working with fuels, or rolling a vehicle to a launch pad, or whenever personnel will work outside and would be at-risk from lightning. These organizations would benefit greatly if the 45 WS could provide more accurate timing of the first CG lightning strike of the day. The Applied Meteorology Unit (AMU) has made significant improvements in forecasting the probability of lightning for the day, but forecasting the time of the first CG lightning with confidence has remained a challenge. To address this issue, the 45 WS requested the AMU to determine if flow regimes, wind speed categories, or a combination of the two could be used to forecast the timing of the first strike of the day in the Kennedy Space Center (KSC)/Cape Canaveral Air Force Station (CCAFS) lightning warning circles. The data was stratified by various sea breeze flow regimes and speed categories in the surface to 5,000-ft layer. The surface to 5,000-ft layer was selected since that is the layer the 45 WS uses to predict the behavior of sea breeze fronts, which are the dominant influence on the occurrence of first lightning in Florida during the warm season. Due to small data sample sizes after stratification, the AMU could not determine a statistical relationship between flow regimes or speed categories and the time of the first CG strike.. As expected, although the amount and timing of lightning activity varies by time of day based on the flow regimes and speed categories, there are extended tails of low lightning activity making it difficult to specify times when the threat of the first lightning flash can be avoided. However, the AMU developed a graphical user interface with input from the 45 WS that allows forecasters to visualize the climatological frequencies of the timing of the first lightning strike. This tool should contribute directly to the 45 WS goal of improving lightning timing capability for its NASA, US Air Force and commercial customers.

statistical analysis↗

Monitoring Land Surface Albedo and Vegetation Dynamics Using High Spatial and Temporal Resolution Synthetic Time Series from Landsat and the MODIS BRDF/NBAR/Albedo Product

Seasonal vegetation phenology can significantly alter surface albedo which in turn affects the global energy balance and the albedo warmingcooling feedbacks that impact climate change. To monitor and quantify the surface dynamics of heterogeneous landscapes, high temporal and spatial resolution synthetic time series of albedo and the enhanced vegetation index (EVI) were generated from the 500-meter Moderate Resolution Imaging Spectroradiometer (MODIS) operational Collection V006 daily BRDF (Bidirectional Reflectance Distribution Function) / NBAR (Nadir BRDF-Adjusted Reflectance) / albedo products and 30-meter Landsat 5 albedo and near-nadir reflectance data through the use of the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM). The traditional Landsat Albedo (Shuai et al., 2011) makes use of the MODIS BRDFAlbedo products (MCD43) by assigning appropriate BRDFs from coincident MODIS products to each Landsat image to generate a 30-meter Landsat albedo product for that acquisition date. The available cloud free Landsat 5 albedos (due to clouds, generated every 16 days at best) were used in conjunction with the daily MODIS albedos to determine the appropriate 30-meter albedos for the intervening daily time steps in this study. These enhanced daily 30-meter spatial resolution synthetic time series were then used to track albedo and vegetation phenology dynamics over three Ameriflux tower sites (Harvard Forest in 2007, Santa Rita in 2011 and Walker Branch in 2005). These Ameriflux sites were chosen as they are all quite nearby new towers coming on line for the National Ecological Observatory Network (NEON), and thus represent locations which will be served by spatially paired albedo measures in the near future. The availability of data from the NEON towers will greatly expand the sources of tower albedometer data available for evaluation of satellite products. At these three Ameriflux tower sites the synthetic time series of broadband shortwave albedos were evaluated using the tower albedo measurements with a Root Mean Square Error (RMSE) less than 0.013 and a bias within the range of 0.006. These synthetic time series provide much greater spatial detail than the 500 meter gridded MODIS data, especially over more heterogeneous surfaces, which improves the efforts to characterize and monitor the spatial variation across species and communities. The mean of the difference between maximum and minimum synthetic time series of albedo within the MODIS pixels over a subset of satellite data of Harvard Forest (16 kilometers by 14 kilometers) was as high as 0.2 during the snow-covered period and reduced to around 0.1 during the snow-free period. Similarly, we have used STARFM to also couple MODIS Nadir BRDF-Adjusted Reflectances (NBAR) values with Landsat 5 reflectances to generate daily synthetic times series of NBAR and thus Enhanced Vegetation Index (NBAR-EVI) at a 30-meter resolution. While normally STARFM is used with directional reflectances, the use of the view angle corrected daily MODIS NBAR values will provide more consistent time series. These synthetic times series of EVI are shown to capture seasonal vegetation dynamics with finer spatial and temporal details, especially over heterogeneous land surfaces.

Vegetation Index↗

Flight Tests of a Remaining Flying Time Prediction System for Small Electric Aircraft in the Presence of Faults

This paper addresses the problem of building trust in the online prediction of a battery powered aircraft's remaining flying time. A series of flight tests is described that make use of a small electric powered unmanned aerial vehicle (eUAV) to verify the performance of the remaining flying time prediction algorithm. The estimate of remaining flying time is used to activate an alarm when the predicted remaining time is two minutes. This notifies the pilot to transition to the landing phase of the flight. A second alarm is activated when the battery charge falls below a specified limit threshold. This threshold is the point at which the battery energy reserve would no longer safely support two repeated aborted landing attempts. During the test series, the motor system is operated with the same predefined timed airspeed profile for each test. To test the robustness of the prediction, half of the tests were performed with, and half were performed without, a simulated powertrain fault. The pilot remotely engages a resistor bank at a specified time during the test flight to simulate a partial powertrain fault. The flying time prediction system is agnostic of the pilot's activation of the fault and must adapt to the vehicle's state. The time at which the limit threshold on battery charge is reached is then used to measure the accuracy of the remaining flying time predictions. Accuracy requirements for the alarms are considered and the results discussed.

Hogge, Edward F.↗

Remaining Flying Time Prediction Implementing Battery Prognostics Framework for Electric UAV's

In this paper the problem of building trust in the online safety prediction of an fixed wing small electric unmanned aerial vehicles (e-UAV) for remaining flying time is addressed. A series of flight tests are described to verify the performance of the remaining flying time prediction algorithm. The estimate of remaining flying time is used to activate an alarm when the predicted remaining time falls below a threshold of two minutes. This updates the pilot to transition to the landing sequence of the flight profile. A second alarm is activated when the battery state of charge (SOC) falls below a specified safety limit threshold. This SOC threshold is the point at which the battery energy reserve would no longer safely support enough aborted landing attempts. During the test flights, the motor system is operated with the same predefined timed airspeed profile for each test. To test the robustness of the developed prediction algorithm, partial tests were performed with and remaining were performed without a simulated power train fault. To simulate a partial power train fault in the e-UAV the pilot engages a resistor bank at a specified time during the test flight. The flying time prediction system is agnostic of the pilot's activation of the fault and must adapt to the vehicle's state. The time at which the limit threshold on battery SOC is reached, it is then used to measure the accuracy of the remaining flying time predictions. This is demonstrated through comparing results from two battery models being developed. Accuracy requirements for the alarms are considered and the results discussed.

Kulkarni, Chetan↗

Thermospheric Heating and Cooling Times During Geomagnetic Storms, Including Extreme Events

We present the first quantitative calculations of thermospheric heating and cooling times for geomagnetic storms of different intensity, including extreme events. We utilize the neutral mass density database of the CHAllenging Mini‐satellite Payload and Gravity Recovery And Climate Experiment missions to produce thermospheric global system response to geomagnetic storms caused by coronal mass ejections via superposed epoch analysis during May 2001 to December 2015. Storm events are grouped in five different categories based on the minimum value of the SYM‐H index. We calculate the time from storm onset for the thermosphere to reach maximum intensification (heating time) and the time from onset for the thermosphere to recover (cooling time). We find that heating and cooling times decrease as storm intensity increases and the effect is more pronounced for the cooling times. For extreme storms, the thermospheric heating time is 9.5 hr, while the cooling time is 22 hr.

Zesta, Eftyhia↗

High-precision x-ray timing of three millisecond pulsars with NICER: Stability estimates and comparison with radio

The Neutron Star Interior Composition Explorer (NICER) is an X-ray astrophysics payload on the International Space Station. It enables unprecedented high-precision timing of millisecond pulsars (MSPs) without the pulse broadening and delays due to dispersion and scattering within the interstellar medium that plague radio timing. We present initial timing results from a year of data on the MSPs PSR B1937+21 and PSR J0218+4232, and nine months of data on PSR B1821−24. NICER time-of-arrival uncertainties for the three pulsars are consistent with theoretical lower bounds and simulations based on their pulse shape templates and average source and background photon count rates. To estimate timing stability, we use the σz measure, which is based on the average of the cubic coefficients of polynomial fits to subsets of timing residuals. So far we are achieving timing stabilities σz≈3×10−14 for PSR B1937+21 and on the order of 10−12 for PSRs B1821−24 and J0218+4232. Within the span of our NICER data we do not yet see the characteristic break point in the slope of σz; detection of such a break would indicate that further improvement in the cumulative root-mean-square timing residual is limited by timing noise. We see this break point in our comparison radio data sets for PSR B1821−24 and PSR B1937+21 on timescales of >2 yr.

Pulsars↗

VALVE TIMING AND AREA PROFILE SELECTION FOR HYDRAULIC PUMPS AND MOTORS

A significant amount of research has been conducted to select valve timing and area profiles that create efficient and quiet hydraulic pumps and motors. Numerous active valve architectures have been modeled and optimized, but the rationale for the final solution is often unclear. The solution is usually highly dependent on the modeled valve geometry constraints and the duty cycle of the pump or motor for which the valve was optimized. This paper first presents an ideal valve area profile and what compromises are made to design for a working pressure and speed range. The ideal valve area profile is resolved using a piston-cylinder model with a pressure dependent bulk modulus and variable valve openings. The piston-cylinder is first modeled as a pump with passive (check) valves and a specified pressure drop across the valves to determine the ideal area profile. A generalization of the ideal timing is demonstrated for positive pressure differential and positive rotation direction, also known as the first quadrant. The effects on valve timing when varying pressure are shown for the first quadrant, but the trends can be extrapolated for all quadrants of operation. The ideal valve profile is implemented as fixed valve timing, meaning the valve area is only a function of the timing angle, in a pump-motor. The pump-motor is simulated in one rotation direction through a pressure range. Performance is high in pumping operation, but when the pressure differential is reversed, cylinder pressure spikes ensue. Two strategies to modify an ideal valve area profile are presented: timing grooves and a pressure shifted valve timing. Timing grooves reduce pressure spikes and cavitation in the cylinder but generally increase throttling losses. A pressure shifted valve timing has lower throttling energy losses, making it the favored solution.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Adaptive time scaling for accelerating explicit finite element analysis

A method for accelerating an explicit finite element analysis (FEA) simulation of a modeled system or process includes performing an initial iteration of the FEA simulation according to a baseline time interval via an FEA computing network, and calculating a criteria ratio of a predetermined set of scaling criteria for the modeled system or process. The method includes determining a time-scaling factor using the criteria ratio via the FEA computing network as a function of the criteria ratio, and then applying the time-scaling factor to the baseline time interval to generate a scaled time interval. The scaled time interval accelerates simulation time of the FEA simulation. The method includes performing a subsequent iteration of the explicit FEA simulation at the scaled time interval using the FEA computing network. The process continues for subsequent iterations, with the time-scaling factor adapting with each iteration.

Chen, Jian↗

Thermal Neutron Die-Way-Time Studies for P and DGNAA of Radioactive Waste Drums at the MEDINA Facility

In Germany, radioactive waste with negligible heat production has to pass through a process of quality checking in order to check its conformance with national regulations prior to its transport, intermediate storage and final disposal. Additionally to its radioactive components, the waste may contain non-radioactive chemically toxic substances that can adversely affect human health and pollute the environment, especially the ground water. After an adequate decay time, the waste radioactivity will become harmless but the non-radioactive substances will persist over time. In principle, these hazardous substances may be quantified from traceability and quality controls performed during the production of the waste packages. As a consequence, a research and development program was initiated in 2007 with the aim to develop a nondestructive analytical technique for radioactive waste packages based on prompt and delayed gamma neutron activation analysis (P and DGNAA) employing a DT-neutron generator in pulsed mode. In a preliminary study it was experimentally demonstrated that P and DGNAA is suitable to determine the chemical composition of large samples. In 2010 a facility called MEDINA (Multi Element Detection based on Instrumental Neutron Activation) was developed for the qualitative and quantitative determination of nonradioactive, toxic elements and substances in 200-l steel drums. The determination of hazardous substances and elements is generally achieved measuring the prompt gamma-rays induced by thermal neutrons. Additional information about the composition of the waste matrix could be derived measuring the delayed gamma-rays from short life activation products. However a sensitive detection of these delayed gamma-rays requires that thermal neutrons have almost vanished. Therefore, the thermal neutron die-away-time has to be known in order to achieve an optimal discrimination between prompt and delayed gamma-ray spectra acquisition. Measurements Thermal neutron die-away times have been determined for the following cases: a) the empty chamber, b ) an empty 200-l steel drum, for a 200-l steel drum filled c) with concrete d) with polyethylene and e) with a mixture of polyethylene and concrete by measuring the prompt-gamma ray count rate of relevant isotopes like of {sup 1}H, {sup 10}B, {sup 12}C, {sup 28}Si, {sup 35}Cl, {sup 40}Ca and {sup 56}Fe which are emitted from different parts of the facility and the sample. Additionally, the average die-away-time was determined from the total detector count rate. The neutron generator was operated with a neutron emission of 8x10{sup 7} n.s{sup -1}, a neutron pulse with a length of 250 μs and a repetition time of 5 ms. The spectra were acquired between the neutron pulses over t{sub c}=500 μs after a pre-defined waiting time t{sub D} (multiple of 500 μs). The thermal neutron die-away time was ranging between 0.9 ms and 5 ms according to the sample composition. As an example the measured thermal neutron die-away-time Λ [μs] of a drum filled with concrete is presented. Detailed results of this study will be presented and discussed. (authors)

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A non-intrusive parallel-in-time adjoint solver with the XBraid library

In this paper, an adjoint solver for the multigrid-in-time software library XBraid is presented. XBraid provides a non-intrusive approach for simulating unsteady dynamics on multiple processors while parallelizing not only in space but also in the time domain (XBraid: Parallel multigrid in time, http://llnl.gov/casc/xbraid). It applies an iterative multigrid reduction in time algorithm to existing spatially parallel classical time propagators and computes the unsteady solution parallel in time. Techniques from Automatic Differentiation are used to develop a consistent discrete adjoint solver which provides sensitivity information of output quantities with respect to design parameter changes. The adjoint code runs backwards through the primal XBraid actions and accumulates gradient information parallel in time. It is highly non-intrusive as existing adjoint time propagators can easily be integrated through the adjoint interface. The adjoint code is validated on advection-dominated flow with periodic upstream boundary condition. We report it provides similar strong scaling results as the primal XBraid solver and offers great potential for speeding up the overall computational costs for sensitivity analysis using multiple processors.

97 MATHEMATICS AND COMPUTING↗

A time-parallel method for scalable heat transfer simulations of additive manufacturing

Here, a major challenge in simulating the thermal behavior in additive manufacturing processes is the disparate length and time scales between transport phenomena occurring in the melt pool and the component. A common simulation approach relies on spatial decomposition for parallel computing, but due to the nature of heat transfer in AM, where most of the computational expenditure is localized near the melt pool, the computational speedup from spatial parallelization saturates quickly. Therefore, additional parallelism by means of time-domain decomposition is needed to fully take advantage of high-performance computing (HPC) resources. This work introduces a time-parallel method to improve the computational scalability of additive manufacturing simulations on HPC systems, while maintaining high temporal resolution of heat transfer near the melt pool. The method, inspired by the nonlinear paraexp formalism, performs an iterative superposition of nonlinear solutions to the initial value problem, integrating the heat equation across overlapping time-parallel intervals. For a single layer of the NIST AMB2018–01 L7 benchmark problem, the method achieves a 38.51x speedup in wall-clock time with a maximum error in the global temperature solution of 0.99%. This reduces the total solution time from 196.72 min to 5.11 min on 128 nodes of the ORNL Frontier supercomputer. The tradeoff between accuracy and total wall-clock time is investigated and recommendations for time-parallel deployment for AM problems are made.

Additive manufacturing↗

An experimental comparison of a space-time multigrid method with PFASST for a reaction-diffusion problem

We consider two parallel-in-time approaches applied to a (reaction) diffusion problem, possibly non-linear. In particular, we consider PFASST (Parallel Full Approximation Scheme in Space and Time) and space-time multigrid strategies. For both approaches, we start from an integral formulation of the continuous time dependent problem. Then, a collocation form for PFASST and a discontinuous Galerkin discretization in time for the space-time multi-grid are employed, resulting in the same discrete solution at the time nodes. Strong and weak scaling of both multilevel strategies are compared for varying orders of the temporal discretization. Moreover, we investigate the respective convergence behavior for non-linear problems and highlight quantitative differences in execution times

97 MATHEMATICS AND COMPUTING↗

Quantifying groundwater travel time near managed recharge operations using 35S as an intrinsic tracer

By identifying groundwater retention times near managed aquifer recharge (MAR) facilities is a high priority for managing water quality, especially for operations that incorporate recycled wastewater. In order to protect public health, California guidelines for Groundwater Replenishment Reuse Projects require a minimum 2–6 month subsurface retention time for recycled water depending on the level of disinfection, which highlights the importance of quantifying groundwater travel times on short time scales. This study developed and evaluated a new intrinsic tracer method using the naturally occurring radioisotope sulfur-35 (35S). The 87.5 day half-life of 35S is ideal for investigating groundwater travel times on the <1 year timescale of interest to MAR managers. Natural concentrations of 35S found in water as dissolved sulfate (35SO4) were measured in source waters and groundwater at the Rio Hondo Spreading Grounds in Los Angeles County, CA, and Orange County Groundwater Recharge Facilities in Orange County, CA. 35SO4 travel times are comparable to travel times determined by well-established deliberate tracer studies. The study also revealed that 35SO4 in MAR source water can vary seasonally and therefore careful characterization of 35SO4 is needed to accurately quantify groundwater travel time. But, more data is needed to fully assess whether or not this tracer could become a valuable tool for managers.

54 ENVIRONMENTAL SCIENCES↗

Parallel-In-Time Multigrid with Adaptive Spatial Coarsening for The Linear Advection and Inviscid Burgers Equations

We apply a multigrid reduction-in-time (MGRIT) algorithm to hyperbolic partial differential equations in one spatial dimension. This study is motivated by the observation that sequential time-stepping is a computational bottleneck when attempting to implement highly concurrent algorithms; thus parallel-in-time methods are desirable. MGRIT adds parallelism by using a hierarchy of successively coarser temporal levels to accelerate the solution on the finest level. In the case of explicit time-stepping, spatial coarsening is a suitable approach to ensure that stability conditions are satisfied on all levels, and it may be useful for implicit time-stepping by producing cheaper multigrid cycles. Unfortunately, uniform spatial coarsening results in extremely slow convergence when the wave speed is near zero, even if only locally. We present an adaptive spatial coarsening strategy that addresses this issue for the variable coefficient linear advection equation and the inviscid Burgers equation using first-order explicit or implicit time-stepping methods. Serial numerical results show this method offers significant improvements over uniform coarsening and is convergent for the inviscid Burgers equation with and without shocks. Parallel scaling tests on up to 128K cores indicate that run-time improvements over serial time-stepping strategies are possible when spatial parallelism alone saturates, and that scalability is robust for oscillatory solutions which change on the scale of the grid spacing.

97 MATHEMATICS AND COMPUTING↗

Benchmarking Real-Time Control Platforms Using a Matlab/Simulink Coder with Applications in the Control of DC/AC Switched Power Converters

In the control of DC/AC switched power converters (SPC), one of the most important aspects to be considered is the selection of the real-time control platform. The real-time control platform must be able to meet the high performance efficiency and regulation requirements of the DC/AC SPC, as these typically operate at switching frequencies in the order of kHz to MHz. For this reason, the hardware characteristics of the ADC and PWM, and the processing capacity of the real-time control platform are of vital importance when implementing advanced digital controllers that maintain voltage and current levels within regulatory standards. In this context, we aimed to perform a comparative study of the computation times of different real-time control platforms when implementing different control strategies for DC/AC switched power converters. We also analyzed the impact of the real-time control platforms on the THD of the voltages generated by the DC/AC switched power converters. With the help of this paper, researchers and developers will have criteria to select which real-time control platform to use in real-time control for DC/AC SPC applications.

42 ENGINEERING↗

Kinetic Simulation of Collisional Magnetized Plasmas with Semi-implicit Time Integration

Plasmas with varying collisionalities occur in many applications, such as tokamak edge regions, where the flows are characterized by significant variations in density and temperature. While a kinetic model is necessary for weakly-collisional high-temperature plasmas, high collisionality in colder regions render the equations numerically stiff due to disparate time scales. In this study, we propose an implicit–explicit algorithm for such cases, where the collisional term is integrated implicitly in time, while the advective term is integrated explicitly in time, thus allowing time step sizes that are comparable to the advective time scales. This partitioning results in a more efficient algorithm than those using explicit time integrators, where the time step sizes are constrained by the stiff collisional time scales. Finally, we implement semi-implicit additive Runge–Kutta methods in COGENT, a high-order finite-volume gyrokinetic code and test the accuracy, convergence, and computational cost of these semi-implicit methods for test cases with highly-collisional plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Real-Time Analysis and Signal Optimization for Charge Detection Mass Spectrometry

Charge detection mass spectrometry (CDMS) is an important tool for measuring mass distributions for high mass samples and heterogeneous mixtures. In CDMS, single ions are trapped and their m/z and charge are measured simultaneously. As a single particle technique, the average signal must be optimized to maximize the number of single ion trapping events. If the average signal is too small, most of the trapping events will be empty, and if the average signal is too large, most of the trapping events will contain multiple ions. In recent embodiments, the time domain signal from the trapped ion is analyzed by fast Fourier transforms. The analysis time is much longer that the data collection time which precludes real-time optimization of the experimental conditions. In this paper, we describe the implementation of CDMS with real-time analysis. Processing the data in real time allows the average signal intensities to be dynamically optimized to maximize the number of single ion trapping events. Real-time analysis also allows the experimental settings to be optimized in a timely manner to target specific mass regimes to maximize the useful information content of the measurements. .

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

The precise time-dependent solution of the Fokker–Planck equation with anomalous diffusion

We study the time behavior of the Fokker–Planck equation in Zwanzig’s rule (the backward-Ito’s rule) based on the Langevin equation of Brownian motion with an anomalous diffusion in a complex medium. The diffusion coefficient is a function in momentum space and follows a generalized fluctuation–dissipation relation. We obtain the precise time-dependent analytical solution of the Fokker–Planck equation and at long time the solution approaches to a stationary power-law distribution in nonextensive statistics. As a test, numerically we have demonstrated the accuracy and validity of the time-dependent solution. - Highlights: • The precise time-dependent solution of the Fokker–Planck equation with anomalous diffusion is found. • The anomalous diffusion satisfies a generalized fluctuation–dissipation relation. • At long time the time-dependent solution approaches to a power-law distribution in nonextensive statistics. • Numerically we have demonstrated the accuracy and validity of the time-dependent solution.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗