Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Constraint Reduction”

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

WE-AB-303-06: Combining DAO with MV + KV Optimization to Improve Skin Dose Sparing with Real-Time Fluoroscopy

Purpose: Real-time fluoroscopy may allow for improved patient positioning and tumor tracking, particularly in the treatment of lung tumors. In order to mitigate the effects of the imaging dose, previous studies have demonstrated the effect of including both imaging dose and imaging constraints into the inverse treatment planning object function. That method of combined MV+kV optimization may Result in plans with treatment beams chosen to allow for more gentle imaging beam-on times. Direct-aperture optimization (DAO) is also known to produce treatment plans with fluence maps more conducive to lower beam-on times. Therefore, in this work we demonstrate the feasibility of a combination of DAO and MV+kV optimization for further optimized real-time kV imaging. Methods: Therapeutic and imaging beams were modeled in the EGSnrc Monte Carlo environment, and applied to a patient model for a previously treated lung patient to provide dose influence matrices from DOSXYZnrc. An MV + kV IMRT DAO treatment planning system was developed to compare DAO treatment plans with and without MV+kV optimization. The objective function was optimized using simulated annealing. In order to allow for comparisons between different cases of the stochastically optimized plans, the optimization was repeated twenty times. Results: Across twenty optimizations, combined MV+kV IMRT resulted in an average of 12.8% reduction in peak skin dose. Both non-optimized and MV+kV optimized imaging beams delivered, on average, mean dose of approximately 1 cGy per fraction to the target, with peak doses to target of approximately 6 cGy per fraction. Conclusion: When using DAO, MV+kV optimization is shown to Result in improvements to plan quality in terms of skin dose, when compared to the case of MV optimization with non-optimized kV imaging. The combination of DAO and MV+kV optimization may allow for real-time imaging without excessive imaging dose. Financial support for the work has been provided in part by NIH Grant T32 EB002103, ACS RSG-13-313-01-CCE, and NIH S10 RR021039 and P30 CA14599 grants. The contents of this submission do not necessarily represent the official views of any of the supporting organizations.

60 APPLIED LIFE SCIENCES↗

A Highly Miniaturized Inertial Grade Gyroscope for Space Applications

The evolution of inertial grade gyroscopes for space applications represents well over 50 years of technology development and an investment of hundreds of millions of dollars. The workhorse product which represents the current state-of-the art for commercially available high performance devices is the Litton-Hemishperical Resonator Gyro (HRG) Inertial Measurement Unit (IMU). This product has a performance figure of merit of 0.003 deg/hr bias drift, a volume of 567 cubic inches, weighs 19 pounds, draws about 30 watts and costs over $1 million each. Clearly devices of this magnitude are not conducive to the minimized mass, volume, power, and cost constraints of outer planet missions. An approach to breaking these potential barriers is the use of Microelectromechanical Systems (MEMS) based inertial devices. Although substantially reduced in size, mass power and cost, this approach has produced devices in the tactical performance range of greater than 1 deg/hour bias drift. This level of performance satisfies the preponderance of high market volume requirements such as automotive and tactical munitions but does not meet the limited market quantity requirements for the high precision space based market. Because of the very limited size of the space based market, there is little economic incentive for commercial fabricators of tactical grade devices to address the necessary performance improvements. The Jet Propulsion Laboratory (JPL) in conjunction with Boeing Space Systems (BSS) is addressing this void to satisfy our mutual requirements in this area. The project objective to is to achieve 0.01 deg/hr performance in an IMU which is less than 10 cubic inches in volume, weighs less than 0.5 pounds, draws less than 1 watt and is available in volume production for less than $2500. Reductions of this magnitude will be mission enabling capabilities for a variety of anticipated outer planet mission attributes such as autonomous control and docking, formation flying and robotic outposts. The improved performance will be realized using improved relative precision fabrication, enhanced vibratory drive and sense designs, and statistical data analysis.

Wiberg, D. V.↗

Comparison of box-wing and conventional aircraft mission performance using multidisciplinary analysis and optimization

Box-wing aircraft designs have the potential to achieve significant reductions in fuel consumption. Closed non-planar wing designs have been shown to reduce induced drag and the statically indeterminate wing structure can lead to reduced wing weight. In addition, the streamwise separation of the two main wings can provide the moments necessary for static stability and control, eliminating the weight and aerodynamic drag of a horizontal tail. Proper assessment of the disciplinary interactions in box-wing designs is essential to determine any realistic performance benefits arising from the use of such a configuration. This study analyzes both box-wing and conventional aircraft designed for representative regional-jet missions. A preliminary parametric investigation shows a lift-to-drag ratio advantage for box-wing designs, while a more detailed multidisciplinary study indicates that the requirement to carry the mission fuel in the wings leads to an increase of between 5% and 1% in total fuel burn compared to conventional designs. Furthermore, the multidisciplinary study identified operating conditions where the box-wing can have superior performance to conventional aircraft despite the fuel volume constraint.

42 ENGINEERING↗

Towards real-time monitoring: data assimilated time-lapse full waveform inversion for seismic velocity and uncertainty estimation

SUMMARY Rapid development of time-lapse seismic monitoring instrumentations has made it possible to collect dense time-lapse data for tomographically retrieving time-lapse (even continuous) images of subsurface changes. While traditional time-lapse full waveform inversion (TLFWI) algorithms are designed for sparse time-lapse surveys, they lack of effective temporal constraint on time-lapse data, and, more importantly, lack of the uncertainty estimation of the TLFWI results that is critical for further interpretation. Here, we propose a new data assimilation TLFWI method, using hierarchical matrix powered extended Kalman filter (HiEKF) to quantify the image uncertainty. Compared to existing Kalman filter algorithms, HiEKF allows to store and update a data-sparse representation of the cross-covariance matrices and propagate model errors without expensive operations involving covariance matrices. Hence, HiEKF is computationally efficient and applicable to 3-D TLFWI problems. Then, we reformulate TLFWI in the framework of HiEKF (termed hereafter as TLFWI-HiEKF) to predict time-lapse images of subsurface spatiotemporal velocity changes and simultaneously quantify the uncertainty of the inverted velocity changes over time. We demonstrate the validity and applicability of TLFWI–HiEKF with two realistic CO2 monitoring models derived from Frio-II and Cranfield CO2 injection sites, respectively. In both 2-D and 3-D examples, the inverted high-resolution time-lapse velocity results clearly reveal a continuous velocity reduction due to the injection of CO2. Moreover, the accuracy of the model is increasing over time by assimilating more time-lapse data while the standard deviation is decreasing over lapsed time. We expect TLFWI-HiEKF to be equipped with real-time seismic monitoring systems for continuously imaging the distribution of subsurface gas and fluids in the future large-scale CO2 sequestration experiments and reservoir management.

58 GEOSCIENCES↗

Hardware Evaluation: Abstract Machine Models and Proxy Architectures for Exascale Computing

The strategic and economic impacts of high performance computing (HPC) systems, over the last several decades, have enabled dramatic improvements in manufacturing, design, development, and research, across almost every sector of the U.S. economy, including genetic analysis for agriculture and medicine, oil exploration, electronics, aircraft, and automotive design, and of course national defense technologies. HPC has repeatedly compounded reductions in time-to-solution and time-to-market in each of these areas. In recognition of these strategic and economic impacts, the Department of Energy (DOE) created the Exascale Computing Project (ECP) as a broad-based technology project focusing on simultaneously co-evolving high-performance computing architecture, system software, and application software. While many details of the exascale architectures are undefined, communicating the form and direction of this co-evolution of advances, in each of these three components, across the project, presents a need for a common language through which these advances can be shared. In this document we present a series of abstract representations designed to allow application developers to focus on the aspects of the machine that are important or relevant to performance and code structure. These abstract machine models (AMMs) describe the proposed architectures at the component and system level and are intended as communication aids between application developers and hardware architects during the co-design process. In addition to providing a common language for communication, the constraints of protecting the intellectual property of ECP’s vendor partners requires that the architectural advances be presented in an abstract sense, rather than exposing the specifics of an individual vendor’s designs, while still providing enough detail such that application developers will be able to reason about the performance trade-offs in the design space of their applications.

97 MATHEMATICS AND COMPUTING↗

Deep reinforcement learning to discover multi-fuel injection strategies for compression ignition engines

Over the past several decades, regulation of compression ignition engine emissions has become increasingly stringent as concern about the environmental and health implications of these emissions has grown. These changing constraints have led to a series of new, alternative fuel injection strategies that aim to maintain power output while reducing in-cylinder generated emissions by operating in the low-temperature combustion (LTC) regime. These advanced injection strategies are created and retuned for individual combinations of engine geometry, fuel, and emissions constraints. Deep reinforcement learning has been shown to be an effective alternative to traditional optimization approaches for highly combinatorial control problems, such as discovering the optimal injection schedules for compression ignition engines. In this study, we deploy a previously presented deep reinforcement learning framework to iteratively optimize a series of engine geometries over a range of increasingly strict NO x emissions regulations. We then examine the resulting injection schedules. We discuss the potential for using this deep reinforcement learning framework for fuel selection screening and for discovering unique injection strategies for different engine geometries and future emissions standards.

33 ADVANCED PROPULSION SYSTEMS↗

Exploration Medical System Demonstration Project

A near-Earth Asteroid (NEA) mission will present significant new challenges including hazards to crew health created by exploring a beyond low earth orbit destination, traversing the terrain of asteroid surfaces, and the effects of variable gravity environments. Limited communications with ground-based personnel for diagnosis and consultation of medical events require increased crew autonomy when diagnosing conditions, creating treatment plans, and executing procedures. Scope: The Exploration Medical System Demonstration (EMSD) project will be a test bed on the International Space Station (ISS) to show an end-to-end medical system assisting the Crew Medical Officers (CMO) in optimizing medical care delivery and medical data management during a mission. NEA medical care challenges include resource and resupply constraints limiting the extent to which medical conditions can be treated, inability to evacuate to Earth during many mission phases, and rendering of medical care by a non-clinician. The system demonstrates the integration of medical technologies and medical informatics tools for managing evidence and decision making. Project Objectives: The objectives of the EMSD project are to: a) Reduce and possibly eliminate the time required for a crewmember and ground personnel to manage medical data from one application to another. b) Demonstrate crewmember's ability to access medical data/information via a software solution to assist/aid in the treatment of a medical condition. c) Develop a common data management architecture that can be ubiquitously used to automate repetitive data collection, management, and communications tasks for all crew health and life sciences activities. d) Develop a common data management architecture that allows for scalability, extensibility, and interoperability of data sources and data users. e) Lower total cost of ownership for development and sustainment of peripheral hardware and software that use EMSD for data management f) Provide better crew health via the reduction in crew errors, crew time, and ground time.

Chin, D. A.↗

Constraints on the Physical Properties of GW190814 through Simulations Based on DECam Follow-up Observations by the Dark Energy Survey

On 2019 August 14, the LIGO and Virgo Collaborations detected gravitational waves from a black hole and a 2.6 solar mass compact object, possibly the first neutron star–black hole merger. In search of an optical counterpart, the Dark Energy Survey (DES) obtained deep imaging of the entire 90% confidence level localization area with Blanco/DECam 0, 1, 2, 3, 6, and 16 nights after the merger. Objects with varying brightness were detected by the DES Pipeline, and we systematically reduced the candidate counterparts through catalog matching, light-curve properties, host-galaxy photometric redshifts, Southern Astrophysical Research spectroscopic follow-up observations, and machine-learning-based photometric classification. All candidates were rejected as counterparts to the merger. To quantify the sensitivity of our search, we applied our selection criteria to full light-curve simulations of supernovae and kilonovae as they would appear in the DECam observations. Because the source class of the merger was uncertain, we utilized an agnostic, three-component kilonova model based on tidally disrupted neutron star (NS) ejecta properties to quantify our detection efficiency of a counterpart if the merger included an NS. We find that, if a kilonova occurred during this merger, configurations where the ejected matter is greater than 0.07 solar masses, has lanthanide abundance less than 10 -8.56 , and has a velocity between 0.18c and 0.21c are disfavored at the 2σ level. Furthermore, we estimate that our background reduction methods are capable of associating gravitational wave signals with a detected electromagnetic counterpart at the 4σ level in 95% of future follow-up observations.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Emissions of nitrogen oxides from US urban areas: estimation from Ozone Monitoring Instrument retrievals for 2005–2014

Satellite remote sensing of tropospheric nitrogen dioxide (NO 2 ) can provide valuable information for estimating surface nitrogen oxides (NO x ) emissions. Using an exponentially modified Gaussian (EMG) method and taking into account the effect of wind on observed NO 2 distributions, we estimate 3-year moving-average emissions of summertime NO x from 35 US (United States) urban areas directly from NO 2 retrievals of the Ozone Monitoring Instrument (OMI) during 2005–2014. Following conclusions of previous studies that the EMG method provides robust and accurate emission estimates under strong-wind conditions, we derive top-down NO x emissions from each urban area by applying the EMG method to OMI data with wind speeds greater than 3–5 m s -1 . Meanwhile, we find that OMI NO 2 observations under weak-wind conditions (i.e., < 3 m s &minus;1 ) are qualitatively better correlated to the surface NO x source strength in comparison to all-wind OMI maps; therefore, we use them to calculate the satellite-observed NO 2 burdens of urban areas and compare with NO x emission estimates. The EMG results show that OMI-derived NO x emissions are highly correlated ( R > 0.93) with weak-wind OMI NO 2 burdens as well as with bottom-up NO x emission estimates over 35 urban areas, implying a linear response of the OMI observations to surface emissions under weak-wind conditions. The simultaneous EMG-obtained effective NO 2 lifetimes (~ 3.5 ± 1.3 h), however, are biased low in comparison to the summertime NO 2 chemical lifetimes. In general, isolated urban areas with NO x emission intensities greater than ~ 2 Mg h -1 produce statistically significant weak-wind signals in 3-year average OMI data. From 2005 to 2014, we estimate that total OMI-derived NO x emissions over all selected US urban areas decreased by 49 %, consistent with reductions of 43, 47, 49, and 44 % in the total bottom-up NO x emissions, the sum of weak-wind OMI NO 2 columns, the total weak-wind OMI NO 2 burdens, and the averaged NO 2 concentrations, respectively, reflecting the success of NO x control programs for both mobile sources and power plants. The decrease rates of these NO x -related quantities are found to be faster (i.e., -6.8 to -9.3 % yr &minus;1 ) before 2010 and slower (i.e., -3.4 to -4.9 % yr &minus;1 ) after 2010. For individual urban areas, we calculate the R values of pair-wise trends among the OMI-derived and bottom-up NO x emissions, the weak-wind OMI NO 2 burdens, and ground-based NO 2 measurements, and high correlations are found for all urban areas (median R = 0.8), particularly large ones ( R up to 0.97). The results of the current work indicate that using the EMG method and considering the wind effect, the OMI data allow for the estimation of NO x emissions from urban areas and the direct constraint of emission trends with reasonable accuracy.

54 ENVIRONMENTAL SCIENCES↗

Emissions of nitrogen oxides from US urban areas: estimation from Ozone Monitoring Instrument retrievals for 2005-2014

Satellite remote sensing of tropospheric nitrogen dioxide (NO 2 ) can provide valuable information for estimating surface nitrogen oxides (NO x ) emissions. Using an exponentially-modified Gaussian (EMG) method and taking into account the effect of wind on observed NO 2 distributions, we estimate three-year moving-average emissions of summertime NO x from 35 US urban areas directly from NO 2 retrievals of the Ozone Monitoring Instrument (OMI) during 2005–2014. Following the conclusions of previous studies that the EMG method provides robust and accurate emission estimates under strong-wind conditions, we derive top-down NO x emissions from each urban area by applying the EMG method to OMI data with wind speeds greater than 3–5 m s -1 . Meanwhile, we find that OMI NO 2 observations under weak-wind conditions (i.e., < 3 m s -1 ) are qualitatively better correlated with the surface NO x source strength in comparison to all-wind OMI maps; and therefore we use them to calculate the satellite-observed NO 2 burdens of urban areas and compare with NO x emission estimates. The EMG results show that OMI-derived NO x emissions are highly correlated ( R > 0.93) with weak-wind OMI NO 2 burdens as well as bottom-up NO x emission estimates over 35 urban areas, implying a linear response of the OMI observations to surface emissions under weak-wind conditions. The simultaneous, EMG-obtained, effective NO 2 lifetimes (~3.5 ± 1.3 h), however, are biased low in comparison to the summertime NO 2 chemical lifetimes. In general, isolated urban areas with NO x emission intensities greater than ~ 2 Mg h -1 produce statistically significant weak-wind signals in three-year average OMI data. From 2005 to 2014, we estimate that total OMI-derived NO x emissions over all selected US urban areas decreased by 49%, consistent with reductions of 43, 47, 49, and 44% in the total bottom-up NO x emissions, the sum of weak-wind OMI NO 2 columns, the total weak-wind OMI NO 2 burdens, and the averaged NO 2 concentrations, respectively, reflecting the success of NO x control programs for both mobile sources and power plants. The decrease rates of these NO x -related quantities are found to be faster (i.e., -6.8 to -9.3% yr -1 ) before 2010 and slower (i.e., -3.4 to -4.9% yr -1 ) after 2010. For individual urban areas, we calculate the R values of pair-wise trends among the OMI-derived and bottom-up NO x emissions, the weak-wind OMI NO 2 burdens, and ground-based NO 2 measurements; and high correlations are found for all urban areas (median R = 0.8), particularly large ones ( R up to 0.97). The results of the current work indicate that using the EMG method and considering the wind effect, the OMI data allow for the estimation of NO x emissions from urban areas and the direct constraint of emission trends with reasonable accuracy.

54 ENVIRONMENTAL SCIENCES↗

The Importance of Addressing Disagreements Between Nominal and Effective Treatments During Bat Mortality Minimization Validation Studies

As the wind energy industry grows, so too does our need for effective and low-cost bat mortality minimization solutions. Despite knowledge gaps in our understanding of what drives bats to collide with spinning turbines, minimization solutions have shown success during validation studies. Curtailment has been consistently effective, reducing bat mortality from 33%-79%, depending on the curtailment scenario and species present. Several validation studies have demonstrated deterrent solutions, specifically ultrasonic deterrents, to be effective, but overall, more variable compared to curtailment solutions. Further, a study combining curtailment with ultrasonic deterrents produced encouraging results, finding that adding ultrasonic deterrents to turbines that were designated to curtail, significantly reduced mortality rates compared to curtailment only turbines. Despite positive results, the strength of inference achieved with validation studies is limited by methodological constraints associated with appropriately assigning fatalities to treatments. Because of the temporal separation between when a bat collides with a turbine and when it is discovered during a ground-based carcass survey, we must assign carcasses to a treatment that ran during the prior night(s). This process may introduce errors at multiple stages. First, for studies that rotate treatments among turbines, mortality surveyors must be confident that carcasses are 'fresh' such that mortalities are correctly assigned to the treatment from the previous night(s). Second, we must reconcile any misalignments between how we assign treatments and how we implement treatments (e.g., turbines assigned with a deterrent treatment may be implemented as another treatment when deterrent devices are not operating as designed). Finally, for curtailment solutions, it is critical to recognize that measured effects are a function of not just the treatment as implemented but the proportion of the night the treatment is realized (e.g., if wind speeds are greater than the curtailment treatment cut-in speed for the entire night, we cannot expect there to me any differences in mortality reduction relative to control turbines). Using a dataset collected between June and October 2017, that rotated 3 treatments (Deterrent Only, 5 m/s Curtail only, and 5 m/s Curtail & Deterrent) and 1 control condition across 16 turbines each night, we explored the importance and implications of accounting for potential errors in assigning fatalities to treatments. We present preliminary results comparing the mortality associated with error-corrected treatments and control conditions highlighting how the measured effect of a treatment (mortality) greatly depends on site specific implementation.

bats↗

MO-AB-BRA-08: A Modular Multi-Source X-Ray Tube for Novel Computed Tomography Applications

Purpose: To develop and build a practical implementation of an x-ray line source for the rapidly increasing number of multi-source imaging applications in CT. Methods: An innovative x-ray tube was designed using CST Particle Studio, ANSYS, and SolidWorks. A slowly varying magnetic field is synchronized with microsecond gating of multiple thermionic electron sources. Electrostatic simulations were run to optimize the geometry of the optics and prevent electrode arcing. Magnetostatic simulations were used for beam deflection studies and solenoid design. Particle beam trajectories were explored with an emphasis on focusing, acceleration, deflection, and space charge effects. Thermal constraints were analyzed for both transient and steady-state regimes. Electromagnetic simulations informed the design of a prototype unit under construction. Results: Particle tracking simulations for a benchtop system demonstrate that three 80 keV electron beams are able to be finely controlled and laterally swept a combined distance of 15 cm over a stationary target with an oscillating magnetic field in the hundreds of gauss. The beams are pulsed according to scanning sequences developed for implementation in a mock stationary CT scanner capable of a 30 ms temporal resolution. Beam spot diameters are approximately 1 mm for 30 mA beams and the stationary target stays well within thermal limits. The relevant hardware and control circuits were developed for incorporation into a physical prototype. Conclusion: A new multi-source x-ray tube was designed in a modular form factor to push the barriers of high-speed CT and spur growth in emerging imaging applications. This technology can be used as the basis for a stationary high-speed CT scanner, a system for generating a virtual fan-beam for dose reduction, or for reducing scatter radiation in cone-beam CT utilizing a tetrahedron beam CT geometry. A 2.4 kW benchtop system is currently being built to show proof of concept for the tube. Support for this research was provided by the University of Wisconsin Madison, Office of the Vice Chancellor for Research and Graduate Education with funding from the Wisconsin Alumni Research Foundation.

60 APPLIED LIFE SCIENCES↗

Reduction of Marine Magnetic Data for Modeling the Main Field of the Earth

The marine data set archived at the National Geophysical Data Center (NGDC) consists of shipborne surveys conducted by various institutes worldwide. This data set spans four decades (1953, 1958, 1960-1987), and contains almost 13 million total intensity observations. These are often less than 1 km apart. These typically measure seafloor spreading anomalies with amplitudes of several hundred nanotesla (nT) which, since they originate in the crust, interfere with main field modeling. The source for these short wavelength features are confined within the magnetic crust (i.e., sources above the Curie isotherm). The main field, on the other hand, is of much longer wavelengths and originates within the earth's core. It is desirable to extract the long wavelength information from the marine data set for use in modeling the main field. This can be accomplished by averaging the data along the track. In addition, those data which are measured during periods of magnetic disturbance can be identified and eliminated. Thus, it should be possible to create a data set which has worldwide data distribution, spans several decades, is not contaminated with short wavelengths of the crustal field or with magnetic storm noise, and which is limited enough in size to be manageable for the main field modeling. The along track filtering described above has proved to be an effective means of condensing large numbers of shipborne magnetic data into a manageable and meaningful data set for main field modeling. Its simplicity and ability to adequately handle varying spatial and sampling constraints has outweighed consideration of more sophisticated approaches. This filtering technique also provides the benefits of smoothing out short wavelength crustal anomalies, discarding data recorded during magnetically noisy periods, and assigning reasonable error estimates to be used in the least square modeling. A useful data set now exists which spans 1953-1987.

Baldwin, R. T.↗

Improved Diurnal Cycle of Precipitation in E3SM With a Revised Convective Triggering Function

Abstract We revise the convective triggering function in Department of Energy's Energy Exascale Earth System Model (E3SM) Atmosphere Model version 1 (EAMv1) by introducing a dynamic constraint on the initiation of convection that emulates the collective dynamical effects to prevent convection from being triggered too frequently and allowing air parcels to launch above the boundary layer to capture nocturnal elevated convection. The former is referred to as the dynamic Convective Available Potential Energy (dCAPE) trigger and the latter as the Unrestricted Launch Level (ULL) trigger. Compared to the original trigger in EAMv1 that initiates convection whenever CAPE is larger than a threshold, the revised trigger substantially improves the simulated diurnal cycle of precipitation over both midlatitude and tropical lands. The nocturnal peak of precipitation and the eastward propagation of convection downstream of the Rockies and over the adjacent Great Plains are much better captured than those in the default model. The overall impact on mean precipitation is minor with some notable improvements over the Indo‐Western Pacific, subtropical Pacific and Atlantic, and South America. In general, the dCAPE trigger helps to better capture late afternoon rainfall peak, while ULL is key to capturing nocturnal elevated convection and the eastward propagation of convection. The dCAPE trigger also primarily contributes to the considerable reduction of convective precipitation over subtropical regions and the frequency of light‐to‐moderate precipitation occurrence. However, no clear improvement is seen in intense convection and the amplitude of diurnal precipitation.

54 ENVIRONMENTAL SCIENCES↗

Implicit–explicit (IMEX) Runge–Kutta methods for non-hydrostatic atmospheric models

Abstract. The efficient simulation of non-hydrostatic atmospheric dynamics requires time integration methods capable of overcoming the explicit stability constraints on time step size arising from acoustic waves. In this work, we investigate various implicit–explicit (IMEX) additive Runge–Kutta (ARK) methods for evolving acoustic waves implicitly to enable larger time step sizes in a global non-hydrostatic atmospheric model. The IMEX formulations considered include horizontally explicit – vertically implicit (HEVI) approaches as well as splittings that treat some horizontal dynamics implicitly. In each case, the impact of solving nonlinear systems in each implicit ARK stage in a linearly implicit fashion is also explored. The accuracy and efficiency of the IMEX splittings, ARK methods, and solver options are evaluated on a gravity wave and baroclinic wave test case. HEVI splittings that treat some vertical dynamics explicitly do not show a benefit in solution quality or run time over the most implicit HEVI formulation. While splittings that implicitly evolve some horizontal dynamics increase the maximum stable step size of a method, the gains are insufficient to overcome the additional cost of solving a globally coupled system. Solving implicit stage systems in a linearly implicit manner limits the solver cost but this is offset by a reduction in step size to achieve the desired accuracy for some methods. Overall, the third-order ARS343 and ARK324 methods performed the best, followed by the second-order ARS232 and ARK232 methods.

54 ENVIRONMENTAL SCIENCES↗

An Evaluation of Controller and Pilot Performance, Workload and Acceptability under a NextGen Concept for Dynamic Weather Adapted Arrival Routing

In todays terminal operations, controller workload increases and throughput decreases when fixed standard terminal arrival routes (STARs) are impacted by storms. To circumvent this operational constraint, Prete, Krozel, Mitchell, Kim and Zou (2008) proposed to use automation to dynamically adapt arrival and departure routing based on weather predictions. The present study examined this proposal in the context of a NextGen trajectory-based operation concept, focusing on the acceptability and its effect on the controllers ability to manage traffic flows. Six controllers and twelve transport pilots participated in a human-in-the-loop simulation of arrival operations into Louisville International Airport with interval management requirements. Three types of routing structures were used: Static STARs (similar to current routing, which require the trajectories of individual aircraft to be modified to avoid the weather), Dynamic routing (automated adaptive routing around weather), and Dynamic Adjusted routing (automated adaptive routing around weather with aircraft entry time adjusted to account for differences in route length). Spacing Responsibility, whether responsibility for interval management resided with the controllers (as today), or resided with the pilot (who used a flight deck based automated spacing algorithm), was also manipulated. Dynamic routing as a whole was rated superior to static routing, especially by pilots, both in terms of workload reduction and flight path safety. A downside of using dynamic routing was that the paths flown in the dynamic conditions tended to be somewhat longer than the paths flown in the static condition.

Johnson, Walter W.↗

Low-dose high-resolution $\mathrm{TOF}$-$\mathrm{PET}$ using ionization-activated multi-state low-$\mathrm{Z}$ detector media

We propose PET scanners using low atomic number media that undergo a persistent local change of state along the paths of the Compton recoil electrons. Measurement of the individual scattering locations and angles, deposited energies, and recoil electron directions allows using the kinematical constraints of the 2-body Compton scattering process to perform a statistical time-ordering of the scatterings, with a high probability of precisely identifying where the gamma first interacted in the detector. In these cases the Line-of-Response is measured with high resolution, determined by the underlying physics processes and not the detector segmentation. There are multiple such media that act through different mechanisms. As an example in which the change of state is quantum-mechanical through a change in molecular configuration, rather than thermodynamic, as in a bubble chamber, we present simulations of a two-state photoswitchable organic dye, a ‘Switchillator’, that is activated to a fluorescent-capable state by the ionization of the recoil electrons. The activated state is persistent, and can be optically excited multiple times to image individual activated molecules. Energy resolution is provided by counting the activated molecules. Location along the LOR is implemented by large-area time-of-flight MCP-PMT photodetectors with single photon time resolution in the tens of ps and sub-mm spatial resolution. Simulations indicate a large reduction of dose.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Comparative compressibility of hydrous wadsleyite and ringwoodite: Effect of H 2 O and implications for detecting water in the transition zone: EFFECT OF WATER ON D K /D P OF WADSLEYITE

Review of recent mineral physics literature shows consistent trends for the influence of Fe and H 2 O on the bulk modulus (K0) of wadsleyite and ringwoodite, the major phases of Earth's mantle transition zone (410–660 km). However, there is little consensus on the first pressure derivative, K0' = (dK/dP)P=0, which ranges from about 4 to >5 across experimental studies and compositions. Here we demonstrate the importance of K0' in evaluating the bulk sound velocity of the transition zone in terms of water content and provide new constraints on the effect of H 2 O on K0' for wadsleyite and ringwoodite by conducting a comparative compressibility study. In the experiment, multiple crystals of hydrous Fo90 wadsleyite containing 2.0 and 0.25 wt % H 2 O were loaded into the same diamond anvil cell, along with hydrous ringwoodite containing 1.4 wt % H 2 O. By measuring their pressure-volume evolution simultaneously up to 32 GPa, we constrain the difference in K0' independent of the pressure scale, finding that H 2 O has no effect on K0', whereas the effect of H 2 O on K0 is significant. The fitted K0' values of hydrous wadsleyite (0.25 and 2.0 wt % H 2 O) and hydrous ringwoodite (1.4 wt % H 2 O) examined in this study were found to be identical within uncertainty, with K0' ~3.7(2). New secondary-ion mass spectrometry measurements of the H 2 O content of these and previously investigated wadsleyite samples shows the bulk modulus of wadsleyite is reduced by 7.0(5) GPa/wt % H 2 O, independent of Fe content for upper mantle compositions. Because K0' is unaffected by H 2 O, the reduction of bulk sound velocity in very hydrous regions of transition zone is expected to be on the order of 1.6%, which is potentially detectible in high-resolution, regional seismology studies.

58 GEOSCIENCES↗