Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SeN”

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 55 records · Page 3

Atomic-Level Structure of Mesoporous Hexagonal Boron Nitride Determined by High-Resolution Solid-State Multinuclear Magnetic Resonance Spectroscopy and Density Functional Theory Calculations

Mesoporous hexagonal boron nitride (p-BN) has received significant attention over the last decade as a promising candidate for water cleaning/pollutant removal and hydrogen storage applications. In this work, high-resolution solid-state NMR spectroscopy and plane-wave density-functional theory (DFT) calculations are used to obtain an atomic-level description of p-BN. 1 H– 15 N or 1 H– 14 N heteronuclear (HETCOR) correlation experiments recorded with either conventional NMR at room temperature or dynamic nuclear polarization surface-enhanced spectroscopy (DNP-SENS) at ca. 100 K reveal NB 2 H, NBH 2 , NBH 3 + species residing on the edges of BN sheets. Ultra-high field 35.2 T 11 B NMR spectroscopy was used to resolve 11 B NMR signals from BN 3 , BN 2 O x (OH) 1–x (x = 0–1), BNO x (OH) 2–x (x = 0–2), BO x (OH) 3–x (x = 0–3), and BO x (OH) 4–x – (x = 0–4). Importantly, 2D 11 B dipolar double-quantum–single-quantum homonuclear correlation spectra reveal that many pore/defect sites are composed of boron oxide/hydroxide clusters connected to the BN framework through BN 2 O units. 1D and 2D 11 B{ 15 N} HETCOR NMR experiments, in addition to plane-wave DFT calculations of nine different structural models, further confirm the assignment of all NMR signals. The detailed structure determination of the pore and edge/defect sites within p-BN should further enable the rational design and development of next-generation p-BN-based materials. In addition, the techniques outlined here should be applicable to determine structure within other porous and/or boron-based materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Further developments on observations of the Taylor instability in a dusty plasma

We previously reported the observation of a hydrodynamic-type interchange instability occurring sporadically at the curved interface between two dust clouds of different densities, in a three-dimensional dusty plasma [Pacha et al., Phys. Plasmas 19, 014501 (2012)]. The mechanism for the instability was not conclusively identified. Further analysis of the video images has been performed, showing that the instability is associated with the appearance of large amplitude dust acoustic waves or shocks. Evidence of a dust acoustic shock at the boundary, which could drive a Richtmyer–Meshkov-type instability, will be presented. We also discuss the possibility, suggested by Avinash and Sen [Phys. Plasmas 22, 083707 (2015)], that the observed instability is due to a Rayleigh–Taylor instability driven by a component of the gravitational acceleration, which is perpendicular to the curved boundary.

Pacha, K. A.↗

Supergravity black holes, Love numbers, and harmonic coordinates

A problem which has now come to a head in general relativity, primarily as a result of the many recent detections of gravitational waves from the coalescence of compact binaries, is that there are very few viable theories against which to “test” it. To perform realistic tests of theories of gravity, we need to be able to look beyond general relativity and evaluate the consistency of a parametrized, physically acceptable, family of black hole metric alternatives with observational data from, especially, gravitational wave detections using, for example, an agnostic Bayesian approach. In this paper we further examine properties of one class of such metrics, which in fact arise as solutions of ungauged supergravity. In particular, we examine the massless, neutral, minimally coupled scalar wave equation in a general stationary, axisymmetric background metric such as that of a charged rotating black hole, when the scalar field is either time independent or in the low-frequency, near-zone limit, with a view to calculating the Love numbers of tidal perturbations, and of obtaining harmonic coordinates for the background metric. For a four-parameter family of charged asymptotically flat rotating black hole solutions of ungauged supergravity theory known as STU black holes, which includes Kaluza-Klein black holes and the Kerr-Sen black hole as special cases, we find that all time-independent solutions, and hence the harmonic coordinates of the metrics, are identical to those of the Kerr solution. In the low-frequency limit we find the scalar fields exhibit the same SL(2,R) symmetry as holds in the case of the Kerr solution. We point out extensions of our results to a wider class of metrics, which includes solutions of Einstein-Maxwell-dilaton theory.

79 ASTRONOMY AND ASTROPHYSICS↗

Automated Vehicle Multi-Object Tracking at Scale with CAN

Millions of vehicles are on the road with RADAR sensors in use for adaptive cruise control (ACC), and RADAR sen- sors are not tracking all of the objects in the field of view. This work shows a work-in-progress tool to improve track- ing from RADAR and controller area network (CAN) which should be vitally useful for safety of transportation systems and automated vehicle development. The CAN data provides object detections, but there is a lingering data association problem. The contribution of this work in progress is the solution to the data association problem by posing the data association as a minimum cost network flow problem, and doing it at low cost with an eye toward scalable CPS research.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Optical O 2 Sensors Also Respond to Redox Active Molecules Commonly Secreted by Bacteria

From a metabolic perspective, molecular oxygen (O 2 ) is arguably the most significant constituent of Earth’s atmosphere. Nearly every facet of microbial physiology is sensitive to the presence and concentration of O 2 , which is the most favorable terminal electron acceptor used by organisms and also a dangerously reactive oxidant. As O 2 has such sweeping implications for physiology, researchers have developed diverse approaches to measure O 2 concentrations in natural and laboratory settings. Recent improvements to phosphorescent O 2 sensors piqued our interest due to the promise of optical measurement of spatiotemporal O 2 dynamics. However, we found that our preferred bacterial model, Pseudomonas aeruginosa PA14, secretes more than one molecule that quenches such sen sors, complicating O 2 measurements in PA14 cultures and biofilms. Assaying supernatants from cultures of 9 bacterial species demonstrated that this phenotype is common: all super natants quenched a soluble O 2 probe substantially. Phosphorescent O 2 probes are often embedded in solid support for protection, but an embedded probe called O 2 NS was quenched by most supernatants as well. Measurements using pure compounds indicated that quenching is due to interactions with redox-active small molecules, including phena zines and flavins. Uncharged and weakly polar molecules like pyocyanin were especially potent quenchers of O 2 NS. These findings underscore that optical O 2 measurements made in the presence of bacteria should be carefully controlled to ensure that O 2 , and not bacterial secretions, is measured, and motivate the design of custom O 2 probes for specific organisms to circumvent sensitivity to redox-active metabolites.

59 BASIC BIOLOGICAL SCIENCES↗

SMC 2021 : Analyzing Resource Utilization and User Behavior on Titan Supercomputer

Resource utilization statistics of submitted jobs on a supercomputer can help us understand how users from various scientific domains use HPC platforms and better design a job scheduler. We explore to generate insight regarding workload distribution and usage pattern domains from job scheduler trace, GPU failure information, and project-specific information collected from Titan supercomputer. Furthermore, we want to know how the scheduler performance varies over time and how the users' scheduling behavior changes following a system failure. These observations have the potential to provide valuable insight, which is helpful to prepare for system failures. These practices will help us develop and apply novel machine learning algorithms in understanding system behavior, requirement, and better scheduling of HPC systems. There are two datasets, RUR and GPU. RUR: This dataset is the job scheduler traces collected from the Titan supercomputerfrom 01/01/2015 to 07/31/2019 (2015.csv - 2019.csv). These were collected usingResource Utilization Report (RUR), a Cray-developed resource-usage data collectionand reporting system. It contains the usage information of its critical resources (CPU,Memory, GPU, and I/O) of each running job on Titan during that period [2]. ProjectAreas: Every job is associated with a project ID. TheProjectAreas.csvdatasetprovides a mapping of the project ID to its domain science. GPU: There have been some hardware-related issues in the GPUs in Titan that caused some GPUs to fail, sometimes irrecoverably during some job runs. This dataset provides information regarding these failures during the execution of the submitted jobs. GPUs on Titan are uniquely identified by a serial number (SN), and they are installed in a location. A GPU can be installed in a location, then removed from that location following a failure, and then re-installed in a different location after fixing the problem. If the failure can't be recovered, the GPU might be removed entirely from Titan. There are two prominent types of failures that resulted in the removal of GPUs from Titan: Double Bit Error (DBE) and Out of the Bus (OTB). The dataset (gc_full.csv) has the following fields: 1. SN : Serial number of a GPU 2. location : The location where it is installed 3. insert : The time when it was inserted into that location 4. remove : The time when it was removed from that location 5. duration : Amount of time the GPU spent in this location 6. out : If the device was taken out entirely w/o a re-installment into a new location. 7. event : If the GPU was taken out entirely, the reason for its removal. To learn more about this dataset, please refer to the git repositoryhttps://github.com/olcf/TitanGPULifeand the related publication [1]. References [1] George Ostrouchov, Don Maxwell, Rizwan A Ashraf, Christian Engelmann, MallikarjunShankar, and James H Rogers. Gpu lifetimes on titan supercomputer: Survival analysisand reliability. InSC20: International Conference for High Performance Computing,Networking, Storage and Analysis, pages 1-14. IEEE, 2020. [2] Feiyi Wang, Sarp Oral, Satyabrata Sen, and Neena Imam. Learning from five-yearresource-utilization data of titan system. In2019 IEEE International Conference onCluster Computing (CLUSTER), pages 1-6. IEEE, 2019.

42 ENGINEERING↗

Interplay of Topography, Fire History, and Climate on Interior Alaska Boreal Forest Vegetation Dynamics in the 21st Century: A Landsat Time-Series Analysis

This study investigates vegetation dynamics in boreal forests of Interior Alaska, focusing on topography, fire history, and climate influences. The study area includes Bonanza Creek Experimental Forest (BCEF) and surrounding region, categorized by topography (upland, floodplain, lowland) and fire history. Using Mann–Kendall trend and Theil–Sen slope analyses on Landsat-derived spectral metrics: Normalized Difference Vegetation Index (NDVI) and Normalized Burn Ratio (NBR), we observed a shift from browning to greening trends, particularly in historically burned areas. The photosynthetic activity in burned upland converged with unburned areas ~30 years post-fire, coincident with a shift towards deciduous dominance during post-fire succession. Normalized Difference Moisture Index (NDMI) trends revealed a significant increase in vegetation moisture content across all topographies. We introduce Effective Seasonal Precipitation Index (ESPI), which combines prior-year annual precipitation with current-year spring snow depth. Its positive correlation with NDMI highlights its potential for monitoring vegetation moisture dynamics at the landscape scale. Furthermore, by correlating dendrochronology-based climate indices, we found strong correlation between NDMI and normalized Supplemental Precipitation Index (nSPI), across topographies. Overall, this research provides critical insights into how climate and fire influence interior boreal vegetation, highlighting the effects of increased precipitation, and topography on shaping differential vegetation responses across the landscape.

Google Earth Engine↗

Liquid Core Detection and Strand Condition Monitoring in a Continuous Caster Using Optical Fiber

Real-time monitoring of the liquid core position during the continuous casting of steel has been demonstrated using low-cost distributed optical-fiber-based strain sensors. These sensors were installed on the containment roll support structures in the segments of a production continuous caster to detect the position of the solid–liquid interface and monitor the strand condition during the continuous casting. Distributed Fiber Bragg Grating sensors (FBGs) were used in this work to monitor strain at six roll positions in the caster. The sensor performance was first validated by comparing optical strain measurements with conventional strain gauge measurements in the lab. Next, optical strain measurements were performed on an isolated caster segment in a segment maintenance facility using hydraulic jacks to simulate the presence of a liquid core under the roll. Finally, the sensors were evaluated during caster operation. The sensors successfully detected the load increase associated with the presence of a liquid core under each instrumented roll location. Incidents of bulging and roll eccentricity were also detected using frequency analysis of the optical strain signal. The liquid core position measurements were compared using predictions from computer models (digital twins) in use at the production site. The measurements were in good agreement with the model predictions, with a few exceptions. Under certain transient caster operating conditions, such as spraying practice changes and SEN exchanges, the model predictions deviated slightly from the liquid core position determined from strain measurements.

47 OTHER INSTRUMENTATION↗

The reflection effect in Hz Hercules

The magnitude brightness variation due to binary star reflection effect is considered, along with tidal distortion effects. The flat region in the Milne-Sen Law is also reviewed, including the nonfit observations. Several studies are included in this binary star program investigation.

Wilson, R. E.↗

Fatigue crack growth at elevated temperature 316 stainless steel and H-13 steel

Crack growths were measured at elevated temperatures under four types of loading: pp, pc, cp, and cc. In H-13 steel, all these four types of loading gave nearly the same crack growth rates, and the length of hold time had negligible effects. In AISI 316 stainless steel, the hold time effects on crack growth rate were negligible if the loading was tension-tension type; however, these effects were significant in reversed bending load, and the crack growth rates under these four types of loading varied considerably. Both tensile and compressive hold times caused increased crack growth rate, but the compressive hold period was more deleterious than the tensile one. Metallographic examination showed that all the crack paths under different types of loading were largely transgranular for both CTS tension-tension specimens and SEN reversed cantilever bending specimens. In addition, an electric potential technique was used to monitor crack growth at elevated temperature.

Chen, W. C.↗

Time-dependent perturbation of a two-state quantum system by a sinusoidal field

Different methods for solving the 'two-level problem' are discussed, namely, the problem of what happens to a material system having only two nondegenerate energy levels when it is perturbed by an electromagnetic field that varies with time in a monochromatic sinusoidal fashion. The various methods discussed include: (1) the Sen Gupta technique using nondegenerate Rayleigh-Schroedinger perturbation theory, (2) the Salwen-Winter-Shirley partitioning perturbation technique, (3) the Shirley and series degenerate Rayleigh-Schroedinger expansion, (4) the degenerate Rayleigh-Schroedinger technique for considering high frequency fields, and (5) the singular perturbation expansion technique.

Dion, D. R.↗

Crack diffusion coefficient - A candidate fracture toughness parameter for short fiber composites

In brittle matrix composites, crack propagation occurs along random trajectories reflecting the heterogeneous nature of the strength field. Considering the crack trajectory as a diffusive process, the 'crack diffusion coefficient' is introduced. From fatigue crack propagation experiments on a set of identical SEN polyester composite specimens, the variance of the crack tip position along the loading axis is found to be a linear function of the effective 'time'. The latter is taken as the effective crack length. The coefficient of proportionality between variance of the crack trajectory and the effective crack length defines the crack diffusion coefficient D which is found in the present study to be 0.165 mm. This parameter reflects the ability of the composite to deviate the crack from the energetically most efficient path and thus links fracture toughness to the microstructure.

Mull, M. A.↗

Elevated temperature crack growth

Alloy 718 crack growth experiments were conducted to assess the ability of the selected path-independent (P-I) integrals to describe the elevated temperature crack growth behavior. These tests were performed on single edge notch (SEN) specimens under displacement control with multiple extensometers to monitor the specimen and crack mouth opening displacement (CMOD). The displacements in these tests were sufficiently high to induce bulk cyclic inelastic deformation of the specimen. Under these conditions, the linear elastic fracture mechanics (LEFM) parameter K does not correlate the crack growth data. The experimentally measured displacement gradients at the end of specimen gage length were used as the boundary conditions in elastic-plastic finite element method (FEM) analyses. These analyses were performed with a node release approach using CYANIDE, a GEAE FEM code, which included a gap element which is capable of efficiently simulating crack closure. Excellent correlation was obtained between the experimentally measured and predicted variation of stress and CMOD with crack length and the stress-CMOD loops for Alloy 718 tests conducted at 538 C. This confirmed the accuracy of the FEM crack growth simulation approach. The experimentally measured crack growth rate data correlated well the selected P-I integrals. These investigations have produced significant progress in developing P-I integrals as non-linear fracture mechanics parameters. The results suggest that this methodology has the potential of accurately describing elevated temperature crack growth behavior under the combined influence of thermal cycling and bulk elastic-inelastic deformation states.

Kim, K. S.↗

Intermediate Temperature Stress Rupture of Woven SiC Fiber, BN Interphase, SiC Matrix Composites in Air

Tensile stress-rupture experiments were performed on woven Hi-Nicalon reinforced SiC matrix composites with BN interphases in air. Modal acoustic emission (AE) was used to monitor the damage accumulation in the composites during the tests and microstructural analysis was performed to determine the amount of matrix cracking that occurred for each sample. Fiber fractograph), was also performed for individual fiber failures at the specimen fracture surface to determine the strengths at which fibers failed. The rupture strengths were significantly worse than what would have been expected front the inherent degradation of the fibers themselves when subjected to similar rupture conditions. At higher applied stresses the rate of rupture "?as larger than at lower applied stresses. It was observed that the change in rupture rate corresponded to the onset of through-thickness cracking in the composites themselves. The primary cause of the sen,ere degradation was the ease with which fibers would bond to one another at their closest separation distances, less than 100 nanometers, when exposed to the environment. The near fiber-to-fiber contact in the woven tows enabled premature fiber failure over large areas of matrix cracks due to the stress-concentrations created b), fibers bonded to one another after one or a few fibers fail. i.e. the loss of global load sharing. An@, improvement in fiber-to-fiber separation of this composite system should result in improved stress- rupture properties. A model was den,eloped in order to predict the rupture life-time for these composites based on the probabilistic nature of indin,idual fiber failure at temperature. the matrix cracking state during the rupture test, and the rate of oxidation into a matrix crack. Also incorporated into the model were estimates of the stress-concentration that would occur between the outer rim of fibers in a load-bearing bundle and the unbridged region of a matrix crack after Xia et al. For the lower stresses, this source of stress-concentration was the likely cause for initial fiber failure that would trigger catastrophic failure of the composite.

Morscher, Gregory N.↗

Thermal Characterization of Functionally Graded Materials: Design of Optimum Experiments

This paper is a study of optimal experiment design applied to the measure of thermal properties in functionally graded materials. As a first step, a material with linearly-varying thermal properties is analyzed, and several different tran- sient experimental designs are discussed. An optimality criterion, based on sen- sitivity coefficients, is used to identify the best experimental design. Simulated experimental results are analyzed to verify that the identified best experiment design has the smallest errors in the estimated parameters. This procedure is general and can be applied to design of experiments for a variety of materials.

Cole, Kevin D.↗

Triboluminescent Materials for Smart Optical Damage Sensors for Space Applications

Triboluminescence is light that is produced by pressure, friction or mechanical shock. New composite materials are constantly being reengi neered in an effort to make lightweight spacecrafts for various NASA missions. For these materials there is interest in monitoring the con dition of the composite in real time to detect any delamination or cr acking due to damage, fatigue or external forces. Methods of periodic inspection of composite structures for mechanical damage such as ult rasonic testing are rather mature. However, there is a need to develop a new technique of damage detection for composites, which could dete ct cracking or delamination from any desired location within a materi al structure in real time. This could provide a valuable tool in the confident use of composite materials for various space applications. Recently, triboluminnescent materials have been proposed as smart sen sors of structural damage. To sense the damage, these materials can b e epoxy bonded or coated in a polymer matrix or embedded in a composi te host structure. When the damage or fracture takes place in the hos t structure, it will lead to the fracture of triboluminescent crystal s resulting in a light emission. This will warn, in real time, that a structural damage has occurred. The triboluminescent emission of the candidate phosphor has to be sufficiently bright, so that the light signal reaching from the point of fracture to the detector through a fiber optic cable is sufficiently strong to be detected. There are a large number of triboluminescent materials, but few satisfy the above criterion. Authors have synthesized a Eu based organic material know n as Europium tetrakis (dibenzoylmethide) triethylammonium .(EuD(sub 4)TEA), one of the bright triboluminescent materials, which is a pote ntial candidate for application as a damage sensor and could be made into a wireless sensor with the addition of microchip, antenna and el ectronics. Preliminary results on the synthesis and characterization of this material shall be presented.

Aggarwal, Mohan D.↗

Climate Change, Nutrition and Food Security in Sub-Saharan Africa

Food security and nutrition in sub-Saharan Africa have long been affected by variations in the weather. Vulnerability to these hazards, along with economic shocks and an adverse political environment, is often uneven in a community. Some individuals and households are more susceptible to emergencies or crises than others, and thus determining who is most vulnerable are and how they are responding to a shock or crises is essential to understand the impact on food security. Daily, quantitative and global observations derived from satellite remote sensing instruments can contribute to understanding how food production has declined due to drought, flood or other weather-related hazard, but it can say nothing about the likelihood that the people living in that area are suffering food insecurity as a result. As Amartya Sen argued, a famine can occur even when there is an absolute surplus of food in a region. Thus organizations like the US Agency for International Development's Famine Early Warning Systems Network (FEWS NET) work to integrate biophysical and socio-economic indicators together with on-the ground assessments to estimate the food security consequences of a variety of events. Climate change is likely to restructure local, regional and global agricultural systems and commodity markets. Although remote sensing information has been used to identify seasonal production declines for the past two decades, new ways of using the data will need to be developed in order to understand, document and respond to the impact of climate change on food security as it is manifested in shorter term shocks. In this article, the contribution of remote sensing is explained, along with the other factors that affect food security

Brown, Molly E.↗

Cruise-Efficient Short Takeoff and Landing (CESTOL): Potential Impact on Air Traffic Operations

The National Aeronautics and Space Administration (NASA) is investigating technological and operational concepts for introducing Cruise-Efficient Short Takeoff and Landing (CESTOL) aircraft into a future US National Airspace System (NAS) civil aviation environment. CESTOL is an aircraft design concept for future use to increase capacity and reduce emissions. CESTOL provides very flexible takeoff, climb, descent and landing performance capabilities and a high-speed cruise capability. In support of NASA, this study is a preliminary examination of the potential operational impact of CESTOL on airport and airspace capacity and delay. The study examines operational impacts at a subject site, Newark Liberty Intemational Airport (KEWR), New Jersey. The study extends these KEWR results to estimate potential impacts on NAS-wide network traffic operations due to the introduction of CESTOL at selected major airports. These are the 34 domestic airports identified in the Federal Aviation Administration's Operational Evolution Plan (OEP). The analysis process uses two fast-time simulation tools to separately model local and NAS-wide air traffic operations using predicted flight schedules for a 24-hour study period in 2016. These tools are the Sen sis AvTerminal model and NASA's Airspace Concept Evaluation System (ACES). We use both to simulate conventional-aircraft-only and CESTOL-mixed-with-conventional-aircraft operations. Both tools apply 4-dimension trajectory modeling to simulate individual flight movement. The study applies AvTerminal to model traffic operations and procedures for en route and terminal arrival and departures to and from KEWR. These AvTerminal applications model existing arrival and departure routes and profiles and runway use configurations, with the assumption jet-powered, large-sized civil CESTOL aircraft use a short runway and standard turboprop arrival and departure procedures. With these rules, the conventional jet and CESTOL aircraft are procedurally separated from each other geographically and in altitude during tenninal airspace approach and departure operations, and each use a different arrival runway. AvTeminal implements its unique Focal-point Scheduling Process to sequence, space and delay aircraft to resolve spacing and overtake conflicts among flights in the airspace and airport system serving KEWR. This Process effectively models integrated arrival and departure operations. AvTerminal assesses acceptance rates and delay magnitude and causality at selected locations, including en route outer boundary fixes, tenninal airspace arrival and departure boundary fixes, terminal airspace arrival merge and departure diverge fixes, and runway landing and takeoff runways. The analysis compares the resulting capacity impacts, flight delays and delay sources between CESTOL and conventional KEWR operations. AvTerminal quantitative results showed that CESTOL has significant capability to increase airport arrival acceptance rates (35-40% at KEWR) by taking advantage of otherwise underused airspace and runways where available. The study extrapolates the AvTerminal-derived KEWR peak arrival and departure acceptance rates to estimate capacity parameter values for each of the OEP airports in the ACES modeling of traffic through the entire NAS network. The extrapolations of acceptance rates allow full, partial or no achievement of CESTOL capacity gains at an OEP airport as determined by assessments of the degree to which local procedures allow leveraging of CESTOL capabilities. These assessments consider each OEP airport's runway geometries, runway system configurations, airport and airspace operations, and potential CESTOL traffic loadings. The ACES modeling, simulates airport and airspace spacing constraints imposed by airport runway system, terminal and en route air traffic control and traffic flow management operations using airport acceptance rates representing conventional-aircraft-only and CESTOL-mixed operations. CEOL aircraft are assumed to have Mach 0.8, and alternatively Mach 0.7, cruise speeds to examine compatibility with conventional aircraft operations in common airspace. The ACES results provides estimates of CESTOL delay impact NAS-wide and at OEP airports due to changes in OEP airport acceptance rates and changes in en route airspace potential conflict rates. Preliminary results show meaningful nationwide delay reductions (20%) due to CESTOL operations at 34 major domestic airports.

Couluris, G. J.↗