Search NASA⌕ Search

SEARCH · Search NASA

Results for “control volume methods”

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 523 records · Page 29

Development of a semigraphitic sulfur-doped ordered mesoporous carbon material for electroanalytical applications

The modification of traditional electrodes with mesoporous carbons is a promising strategy to produce high performance electrodes for electrochemical sensing. The high surface area of mesoporous carbons provides a large number of electroactive sites for binding analytes. Controlling the pore size and structure of mesoporous carbons and modifying their electronic properties via doping offers additional benefits like maximizing transport and tuning the electrochemical processes associated with analyte detection. This work reports a facile method to produce sulfur-doped ordered mesoporous carbon materials (S-OMC) with uniform pore structure, large pore volume, high surface area and semigraphitic structure. The synthesis used thiophenol as a single source of carbon and sulfur, and iron as a catalyst for low temperature carbonization. The S-OMC material was deposited on a glassy carbon electrode and used as a sensor with high sensitivity (11.7 A L mol-1) and selectivity for chloramphenicol detection in presence of other antibiotics. As a proof-of-concept, the sensor was applied to the direct analysis of the drug in reconstituted powdered milk and in commercial eye drops.

36 MATERIALS SCIENCE↗

The ever-evolving active site: transformation of single atoms to extended structures during the Rh-catalyzed reverse water–gas shift reaction

At low temperatures (< 400°C), single atoms of Rh supported on rutile TiO 2 (rTiO 2 ) are responsible for the formation of CO during the reverse water gas shift (RWGS), while methane production is associated with the Rh-TiO 2 interface due to the correlation between methane formation rates and the volume-averaged Rh nanoparticle diameter. As the temperature is increased to >540°C, there is a notable increase in CO selectivity and methane production rates tend towards zero. The time to reach zero depends on the temperature but independent of the initial Rh structure (single atoms and/or nanoparticles), which is controlled by the catalyst preparation method (wetness impregnation versus colloidal nanoparticles). At 600°C and > 4 h time on stream, the catalytic behaviour becomes completely agnostic to the initial Rh structure as well as weight loading, and the catalysts are highly selective for the reverse water gas-shift reaction. Post-reaction HR-TEM image analysis confirms Rh nanoparticles crystallize/order during the reaction; at 400°C, most of the Rh particles are disordered, while at 600°C, they are more ordered (i.e., development of defined faceting). Infrared spectroscopy of CO adsorption on Rh nanoparticles confirms the appearance of defined facets after annealing in nitrogen at high temperatures. Annealing the Rh/rTiO 2 catalysts prior to the RWGS reaction demonstrates the structural transformation of Rh depends only on time and temperature and not on reactant or product fugacity. Sites responsible for stabilizing Rh single atoms are no longer competent at higher temperatures, enabling their integration into existent nanoparticles. As the reaction temperature is increased to temperatures >540°C, the dominant active site for CO production evolves from single atoms to extended Rh structures.

36 MATERIALS SCIENCE↗

Robustness of pairwise kinematic Sunyaev–Zel’dovich effect to optical-cluster-selection bias

The pairwise kinematic Sunyaev–Zel’dovich (kSZ) effect measures both the pairwise motion between galaxy groups and clusters and the amount of gas within them, providing a tracer for cosmic growth. To interpret the cosmological information in the kSZ measurements, it is crucial to understand the optical-cluster-selection bias on the kSZ observables. Line-of-sight structures that contribute to both the optical observable (e.g., richness) and the cosmological signal can induce a correlation between these two quantities at a fixed cluster mass. The selection bias arising from this correlation is a key systematic effect for cosmological analyses. For cosmological observables such as cluster abundance and weak lensing, controlling this selection bias may help explain the tension between the DES-Y1 results and the Planck constraints. In order to test for a kSZ effect equivalent of such a bias, we adopted an alternative mock richness based on galaxy counts within cylindrical volumes along the line of sight. We applied the cylindrical count method to hydrodynamical simulations across a wide range of galaxy-selection criteria, assigning richness consistent with DES-Y1 to the mock clusters. When comparing optically selected clusters to mass-selected halos, we find no significant bias on pairwise kSZ signals, pairwise velocities, or optical depth within our uncertainty limits of approximately 16, 10, and 8%, respectively.

79 ASTRONOMY AND ASTROPHYSICS↗

AGC-1 Irradiation Induced Property Changes Analysis Report: Dynamic Elastic Modulus

Recent analyses of some of the test property data from the AGC-1 creep capsule are reported. Specifically the analyses of data for the effects of compressive creep strain on the Dynamic elastic modulus (by the fundamental frequency and ultrasonic velocity/tof methods) are reported. Statistical “t” testing was used to establish the statistical significance of any differences observed between the Dynamic elastic modulus of the creep (irradiated with compressive stress) and the control (irradiated without stress) specimens. Compressive creep strain does have a small effect on Dynamic elastic modulus, reducing it at the higher doses as volume turnaround is approached.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Effect of perturbations and a meal on superior mesenteric artery flow in patients with orthostatic hypotension

Our aims were to evaluate to role of superior mesenteric blood flow in the pathophysiology of orthostatic hypotension in patients with generalized autonomic failure. METHODS: Twelve patients with symptomatic neurogenic orthostatic hypotension and 12 healthy controls underwent superior mesenteric artery flow measurements using Doppler ultrasonography during head-up tilt and tilt plus meal ingestion. Autonomic failure was assessed using standard tests of the function of the sympathetic adrenergic, cardiovagal and postganglionic sympathetic sudomotor function. RESULTS: Superior mesenteric flow volume and time-averaged velocity were similar in patients and controls at supine rest; however, responses to cold pressor test and upright tilt were attenuated (p < 0.05) in patients compared to controls. Head-up tilt after the meal evoked a profound fall of blood pressure and mesenteric blood flow in the patients; the reduction of mesenteric blood flow correlated (r = 0.89) with the fall of blood pressure in these patients, providing another manifestation of failed baroreflexes. We make the novel finding that the severity of postprandial orthostatic hypotension regressed negatively with the postprandial increase in mesenteric flow in patients with orthostatic hypotension. CONCLUSION: Mesenteric flow is under baroreflex control, which when defective, results in, or worsens orthostatic hypotension. Its large size and baroreflexivity renders it quantitatively important in the maintenance of postural normotension. The effects of orthostatic stress can be significantly attenuated by reducing the splanchnic-mesenteric volume increase in response to food. Evaluation of mesenteric flow in response to eating and head-up tilt provide important information on intra-abdominal sympathetic adrenergic function, and the ability of the patient to cope with orthostatic stress.

Non-NASA Center↗

Standardization of a Volumetric Displacement Measurement for Two-Body Abrasion Scratch Test Data Analysis

A limitation has been identified in the existing test standards used for making controlled, two-body abrasion scratch measurements based solely on the width of the resultant score on the surface of the material. A new, more robust method is proposed for analyzing a surface scratch that takes into account the full three-dimensional profile of the displaced material. To accomplish this, a set of four volume displacement metrics are systematically defined by normalizing the overall surface profile to statistically denote the area of relevance, termed the Zone of Interaction (ZOI). From this baseline, depth of the trough and height of the ploughed material are factored into the overall deformation assessment. Proof of concept data were collected and analyzed to demonstrate the performance of this proposed methodology. This technique takes advantage of advanced imaging capabilities that now allow resolution of the scratched surface to be quantified in greater detail than was previously achievable. A quantified understanding of fundamental particle-material interaction is critical to anticipating how well components can withstand prolonged use in highly abrasive environments, specifically for our intended applications on the surface of the Moon and other planets or asteroids, as well as in similarly demanding, harsh terrestrial settings

Kobrick, Ryan L.↗

Standardization of a Volumetric Displacement Measurement for Two-Body Abrasion Scratch Test Data Analysis

A limitation has been identified in the existing test standards used for making controlled, two-body abrasion scratch measurements based solely on the width of the resultant score on the surface of the material. A new, more robust method is proposed for analyzing a surface scratch that takes into account the full three-dimensional profile of the displaced material. To accomplish this, a set of four volume- displacement metrics was systematically defined by normalizing the overall surface profile to denote statistically the area of relevance, termed the Zone of Interaction. From this baseline, depth of the trough and height of the plowed material are factored into the overall deformation assessment. Proof-of-concept data were collected and analyzed to demonstrate the performance of this proposed methodology. This technique takes advantage of advanced imaging capabilities that allow resolution of the scratched surface to be quantified in greater detail than was previously achievable. When reviewing existing data analysis techniques for conducting two-body abrasive scratch tests, it was found that the ASTM International Standard G 171 specified a generic metric based only on visually determined scratch width as a way to compare abraded materials. A limitation to this method was identified in that the scratch width is based on optical surface measurements, manually defined by approximating the boundaries, but does not consider the three-dimensional volume of material that was displaced. With large, potentially irregular deformations occurring on softer materials, it becomes unclear where to systematically determine the scratch width. Specifically, surface scratches on different samples may look the same from a top view, resulting in an identical scratch width measurement, but may vary in actual penetration depth and/or plowing deformation. Therefore, two different scratch profiles would be measured as having identical abrasion properties, although they differ significantly.

Street, K. W. Jr.↗

WIPP Space Faces TRU Waste Disposal Limits; The Volume Reducing Solution - 18365

Transuranic (TRU) waste packaging for removal from a glovebox has for many years involved the use of bags as the primary containment. In many cases, many hydrogenous layers of bags, yards of tape, and multiple filters are added to the waste stream in order to get the TRU waste into the final 55 gallon drum package to be shipped and buried at the Waste Isolation Pilot Plant (WIPP). The WIPP is currently limited to a total finite disposed volume of Contact Handled (CH) and Remote Handled (RH) TRU waste. As of early 2014 over half of the allotted WIPP volume had been filled with mainly CH-TRU waste. Based on current bag operation packaging methods used, the Department of Energy's (DOE) remaining CH and RH waste will not fit in the remaining limited WIPP space. Bag operations, can be performed safely by good operators following tedious repetitive processes and many Administrative Safety Controls. Since the bag is the only confinement layer between the contamination and the operator, contamination incidents can occur. Bag-out operations require close proximity to the waste to perform the process. With some waste streams, this close to waste proximity increases operator radiation dose. Central Research Laboratories (CRL), based in Red Wing Minnesota, who for many decades has developed proven engineered material transfer solutions based around their Rapid Transfer Port (RTP) technology, has worked with DOE sites to develop a TRU Waste Management Solution. The solution is a single filtered large Polyethylene (PE) or Stainless Steel (SST) Rapid Transfer Port container in a standard Department of Transportation (DOT) 7A, Type A drum, which allows for loading of TRU waste directly from the operational glovebox without the use of bags, sleeves or tapes to provide leak tight containment. All operator interventions to connect and disconnect a drum to and from the glovebox are achieved in an upright ergonomic position by one (1) operator in less than three (3) minutes. Minimal As Low As Reasonably Achievable (ALARA) levels are achieved during connect and disconnect of the drum, due to the engineered safety controls of the RTP. The CRL RTP solution allows glovebox or re-packaged TRU waste to be safely transferred to the final 55 Gallon drum container without additional bags or tape. The high diffusion coefficient filter on the CRL RTP 55G container and elimination of added hydrogenous packaging, improves the Decay Heat Limit (DHL) per drum for most TRU waste types compared with bagged drums, and in most cases allows more TRU waste to be loaded per drum. Inorganic TRU waste materials can be packaged directly into the CRL RTP SST drum container to the maximum DHL and Fissile Gram Equivalent (FGE) allowable for a drum. This paper will show how highly efficient TRAMPAC approved TRU waste packaging of CH and RH TRU waste can be accomplished in a safer, cost effective manner without adding any additional packaging materials to the waste stream, thus increasing the chance of placing all DOE's legacy and to be generated TRU waste in the remaining WIPP space. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Accelerating Thermochemical Equilibrium Calculations for Nuclear Reactor Applications

Thermochemical properties play a key role in modeling and simulation of several key phenomena in nuclear reactors. There has been an increasing interest in incorporating CALPHAD-based formulations in multiphysics simulations including for Molten Salt Reactors where knowledge of phase evolution of the salt and the chemical potentials of various elements are of utmost importance in source term analyses and redox control. However, the size of such simulations is often limited by the high computational cost of full thermodynamic equilibrium calculations. This work discusses the current efforts aimed at accelerating thermochemical equilibrium calculations for multiphysics simulations performed using the open-source finite element / finite volume code Multiphysics Object Oriented Simulation Environment (MOOSE) [1]. While several methods have been proposed for accelerating phase equilibrium calculations [2], most focus on relatively small systems and often rely on a- priori knowledge of the state-space of the system. Nuclear materials, however, are often multi-component systems owing to the evolution of composition under irradiation and an approach based on a-priori mapping of phase diagram is often not enough. This work is aimed at demonstrating an on-the-fly surrogate modeling framework that uses active learning to reduce the number of full equilibrium calculations that must be performed. By combining with efficient coupling approaches, the surrogate framework helps in reducing the computational cost of thermodynamic equilibrium informed multiphysics simulations of nuclear materials. The performance is benchmarked against full coupling with the thermochemistry library Thermochimica [3]. This work uses a machine learning based approach for constructing surrogate models to predict the stable phases in a multicomponent system. The surrogates were constructed using neural networks and Gaussian process classification. In this work, we compare the relative performance of the two methods. We also demonstrate the use of caching previous calculations by interpolating the values from nearest neighbors. References [1] Lindsay, A.D., et al. "2.0 – MOOSE: Enabling massively parallel multiphysics simulation", SoftwareX, 20 (2022): 101202. [2] Roos, W.A. and Zietsman J.H. "Accelerating complex chemical equilibrium calculations – A Review", Calphad, 77 (2022): 102380. [3] Piro, M.H.A., et al. "The thermochemistry library Thermochimica", Computational Materials Science, 67 (2013): 266-272.

36 MATERIALS SCIENCE↗

Comparing Cetane Number Measurement Methods

Cetane number is one of the most important fuel performance metrics for mixing controlled compression-ignition “diesel” engines, quantifying a fuel’s propensity for autoignition when injected into end-of-compression-type temperature and pressure conditions. The historical default and referee method on a Cooperative Fuel Research (CFR) engine configured with indirect fuel injection and variable compression ratio is cetane number (CN) rating. A subject fuel is evaluated against primary reference fuel blends, with heptamethylnonane defining a low-reactivity endpoint of CN = 15 and hexadecane defining a high-reactivity endpoint of CN = 100. While the CN scale covers the range from zero (0) to 100, typical testing is in the range of 30 to 65 CN. Alternatively, several constant-volume combustion chamber (CVCC)-based cetane rating devices have been developed to rate fuels with an equivalent derived cetane number (DCN) or indicated cetane number (ICN). These devices measure ignition delay for fuel injected into a fixed volume of high-temperature and high-pressure air to simulate end-of-compression-type conditions. In this study, a range of novel fuel compounds are evaluated across three CVCC methods: the Ignition Quality Tester (IQT), Fuel Ignition Tester (FIT), and Advanced Fuel Ignition Delay Analyzer (AFIDA). Resulting DCNs and ICNs are compared for fuels within the normal diesel fuel range of reactivity, as well as very high (~100) and very low DCNs/ICNs (~5). Distinct differences between results from various devices are discussed. This is important to consider because some new, high-efficiency advanced compression-ignition (CI) engine combustion strategies operate with more kinetically controlled distributed combustion as opposed to mixing controlled diffusion flames. These advanced combustion strategies may benefit from new fuel chemistries, but current rating methods of CN, DCN, and ICN may not fully describe their performance. In addition, recent evidence suggests ignition delay in modern on-road diesel engines with high-pressure common rail fuel injection systems may no longer directly correlate to traditional CN fuel ratings. Simulated end-of-compression conditions are compared for CN, DCN, and ICN and discussed in the context of modern diesel engines to provide additional insight. Results highlight the potential need for revised and/or multiple fuel test conditions to measure fuel performance for advanced CI strategies.

47 OTHER INSTRUMENTATION↗

Novel estimation of aerosol processes with particle size distribution measurements: a case study with the TOMAS algorithm v1.0.0

Atmospheric aerosol microphysical processes are a significant source of uncertainty in predicting climate change. Specifically, aerosol nucleation, emissions, and growth rates, which are simulated in chemical transport models to predict the particle size distribution, are not understood well. However, long-term size distribution measurements made at several ground-based sites across Europe implicitly contain information about the processes that created those size distributions. This work aims to extract that information by developing and applying an inverse technique to constrain aerosol emissions as well as nucleation and growth rates based on hourly size distribution measurements. We developed an inverse method based upon process control theory into an online estimation technique to scale aerosol nucleation, emissions, and growth so that the model–measurement bias in three measured aerosol properties exponentially decays. The properties, which are calculated from the measured and predicted size distributions, used to constrain aerosol nucleation, emission, and growth rates are the number of particles with a diameter between 3 and 6 nm, the number with a diameter greater than 10 nm, and the total dry volume of aerosol ( N 3–6 , N 10 , V dry ), respectively. In this paper, we focus on developing and applying the estimation methodology in a zero-dimensional “box” model as a proof of concept before applying it to a three-dimensional simulation in subsequent work. The methodology is first tested on a dataset of synthetic and perfect measurements that span diverse environments in which the true particle emissions, growth, and nucleation rates are known. The inverse technique accurately estimates the aerosol microphysical process rates with an average and maximum error of 2 % and 13 %, respectively. Next, we investigate the effect that measurement noise has on the estimated rates. The method is robust to typical instrument noise in the aerosol properties as there is a negligible increase in the bias of the estimated process rates. Finally, the methodology is applied to long-term datasets of in situ size distribution measurements in western Europe from May 2006 through June 2007. At Melpitz, Germany, and Hyytiälä, Finland, the average diurnal profiles of estimated 3 nm particle formation rates are reasonable, having peaks near noon local time with average peak values of 1 and 0.15 cm –3 s –1 , respectively. The normalized absolute error in estimated N 3–6 , N 10 , and V dry at three European measurement sites is less than 15 %, showing that the estimation framework developed here has potential to decrease model–measurement bias while constraining uncertain aerosol microphysical processes.

54 ENVIRONMENTAL SCIENCES↗

System and method for computed axial lithography (CAL) for 3D additive manufacturing

In one aspect the present disclosure relates to a system for forming a three dimensional (3D) object from a volume of photo-responsive material contained within a material container. A plurality of optical projection subsystems are adapted to be arranged circumferentially around a material container containing the volume of photo-responsive material, and are controlled by a controller. The optical projection subsystems direct optical projections at a plurality of angles θ through the volume of photo-responsive material. The optical projections are further directed about a z axis extending through the volume of photo-responsive material. The plurality of optical subsystems are further controlled by the controller to provide each of the optical projections with a calculated three-dimensional intensity distribution acting over a fixed temporal exposure period, which is sufficient to at least one of cure or remove selected portions of the volume of photo-responsive material, and leave other portions unmodified, to form a desired 3D part.

Kelly, Brett↗

System and method for control of compression in internal combustion engine via compression ratio and elastic piston

The present disclosure relates to a system for controlling ignition of an air/fuel mixture intake charge directed into an internal combustion engine. The system may have a longitudinally movable inner cylinder liner configured to fit within a cylinder wall portion of an internal combustion engine, and able to receive a piston of the engine therein. A portion of the inner cylinder liner defines an internal volume forming a combustion chamber, and the internal volume controls a compression ratio of the cylinder. The system also has a cylinder head assembly operatively associated with the inner cylinder liner and able to move linearly to cause longitudinal displacement of the inner cylinder liner relative to the cylinder wall portion. This enables the volume of the combustion chamber to be further varied, to thus further vary the compression ratio.

Flowers, Daniel L.↗

A Second Law Based Unstructured Finite Volume Procedure for Generalized Flow Simulation

An unstructured finite volume procedure has been developed for steady and transient thermo-fluid dynamic analysis of fluid systems and components. The procedure is applicable for a flow network consisting of pipes and various fittings where flow is assumed to be one dimensional. It can also be used to simulate flow in a component by modeling a multi-dimensional flow using the same numerical scheme. The flow domain is discretized into a number of interconnected control volumes located arbitrarily in space. The conservation equations for each control volume account for the transport of mass, momentum and entropy from the neighboring control volumes. In addition, they also include the sources of each conserved variable and time dependent terms. The source term of entropy equation contains entropy generation due to heat transfer and fluid friction. Thermodynamic properties are computed from the equation of state of a real fluid. The system of equations is solved by a hybrid numerical method which is a combination of simultaneous Newton-Raphson and successive substitution schemes. The paper also describes the application and verification of the procedure by comparing its predictions with the analytical and numerical solution of several benchmark problems.

Majumdar, Alok↗

OPERATIONAL TOPOLOGY ESTIMATION OF THREE PHASE UNBALANCED POWER DISTRIBUTION SYSTEMS WITH OUTAGES

Distribution grid, in general, is constructionally loopy but a radially operated network. Opening and closing of suitable tie-switches allows transition from one radial topology to another. Tie-switches ensure reliability and provide restoration paths when faults (outages) occur. The estimation of the system's current operational topology is a precursor to implementing any optimal control actions (during normal operation) or restorative actions (during outage condition). The existing literature in topology detection suffers with several limitations including: (1) computationally intensive nonlinear formulations that are not scalable for large distribution systems, (2) methods requiring exhaustive search for the correct topology by enumerating the complete topological search space, and (3) approaches requiring large volumes of time-series measurements and are not applicable for estimating continuously changing operational topology. Furthermore, the existing literature on distribution grid topology estimation cannot be extended to estimate the operational topology during outage conditions. For operational topology estimation during outages, our objective is to estimate both energized radial feeder and the outaged sections. Starting with a radial topology, these methods identify outaged feeder sections and hence, are inapplicable when the current radial operational topology is not known. In this work, we present a mathematical programming approach to estimate the operational topology of a three-phase unbalanced power distribution system under both outage and normal operating conditions. The aforementioned limitations of existing literature were countered by modeling the problem as a mixed-integer linear programming (MILP) problem which minimizes the errors in available measurements while satisfying the power flow and structural constraints of the system. Specically, a minimum weighted least absolute value estimator is proposed that uses the instantaneous line flow measurements, historical/ forecasted load data, and ping measurements from a limited set of smart meters and solves an MILP to identify the operational topology based on a linearized three-phase unbalanced power flow model. The proposed approach is tested on IEEE 123-bus feeder and 1069 bus multi-feeder test system and the numerical results prove the applicability of the proposed method to large three phase unbalanced grids with high errors in measurements. It is observed that the approach performs well even at high percentages of measurement errors.

GANDLURU NAGA, ANANDINI↗

350 ps Ultrafast room-temperature scintillation realized on CsPbBr 3 -based single crystals via Br 2 over-doping

Ultrafast scintillators are essential for next-generation radiation detection, positron emission tomography, and high-speed medical imaging. All-inorganic CsPbBr 3 perovskites are attractive candidates because of their high stopping power, and excellent optical quality, yet their long carrier lifetimes result in slow scintillation responses on the order of hundreds of nanoseconds. Here, we demonstrate that controlled over-doping with Br 2 produces CsPbBr 3.03 single crystals with sub-nanosecond scintillation at room temperature while preserving crystal quality. Single crystals grown by the Bridgman method exhibit high transparency and maintain the orthorhombic perovskite structure. Br 2 over-doping induces a slight lattice expansion (about 0.42% increase in unit-cell volume) while maintaining the orthorhombic perovskite phase and high optical transparency. Optical absorption reveals a slight redshift of the absorption edge after Br 2 introduction, indicating a modified defect landscape. Time-resolved photoluminescence and radioluminescence measurements show that Br 2 doping creates dense and efficient recombination centers that reduce the scintillation decay time from more than 100 ns in undoped crystals to 350 ps under 5.486 MeV α-particle excitation, and the scintillation decay time decreases by two orders of magnitude. The doped crystals also achieve a spatial resolution of 12 lp mm −1 in X-ray imaging. These results reveal a defect-engineering route for achieving ultrafast scintillation in halide perovskites and highlight the potential of Br 2 -modified CsPbBr 3 for fast timing applications.

Li, Zongxiao [Chinese Academy of Sciences (CAS), N↗

Optimizing Error-Bounded Lossy Compression for Scientific Data With Diverse Constraints

Vast volumes of data are produced by today's scientific simulations and advanced instruments. These data cannot be stored and transferred efficiently because of limited I/O bandwidth, network speed, and storage capacity. Error-bounded lossy compression can be an effective method for addressing these issues: not only can it significantly reduce data size, but it can also control the data distortion based on user-defined error bounds. In practice, many scientific applications have specific requirements or constraints for lossy compression, in order to guarantee that the reconstructed data are valid for post hoc analysis. For example, some datasets contain irrelevant data that should be isolated in particular and users often have intuition regarding value ranges, geospatial regions, and other data subsets that are crucial for subsequent analysis. Existing state-of-the-art error-bounded lossy compressors, however, do not consider these constraints during compression, resulting in inferior compression ratios with respect to user's post hoc analysis, due to the fact that the data itself provides little or no value for post hoc analysis. In this work we address this issue by proposing an optimized framework that can preserve diverse constraints during the error-bounded lossy compression, e.g., cleaning the irrelevant data, efficiently preserving different precision for multiple value intervals, and allowing users to set diverse precision over both regular and irregular regions. We perform our evaluation on a supercomputer with up to 2,100 cores. Experiments with six real-world applications show that our proposed diverse constraints based error-bounded lossy compressor can obtain a higher visual quality or data fidelity on reconstructed data with the same or even higher compression ratios compared with the traditional state-of-the-art compressor SZ. Furthermore, our experiments also demonstrate very good scalability in compression performance compared with the I/O throughput of the parallel file system.

97 MATHEMATICS AND COMPUTING↗

TU-AB-303-03: CBCT-Based On-Site Simulation, Planning, and Delivery (OSPD) for Whole Brain Radiotherapy

Purpose: To demonstrate the feasibility of a CBCT-based on-site simulation, planning, and delivery (OSPD) for whole brain radiotherapy, in which all steps from imaging, planning to treatment delivery are performed at the treatment unit in one appointment time slot. This work serves as the proof of concept for future OSPD single fraction radiation therapy. Methods: An integrated on-site imaging, planning and delivery workflow was developed and tested for whole brain radiotherapy. An automated two-opposed-oblique-beam plan is created by utilizing the treatment planning system scripting and simple field-in-field IMRT. The IMRT plan is designed with maximum 8 control points to cover the target volume consisting of the brain to C1/C2 of the spinal cord, with dose homogeneity criteria from −5% to +7% of the prescription dose. Due to inaccuracy of reconstructed Hounsfield unit numbers in CBCT images, the dose distribution is calculated with non-heterogeneity correction introducing only clinically insignificant dose discrepancy. A coherent and synchronized workflow was designed for a team of attending physician, physicist, therapists, and dosimetrist to work closely with the ability to quickly modify, approve, and implement the treatment. Results: Thirty-one patients have been treated with this OSPD treatment, without compromising the plan quality compared to our regular clinically used parallel apposed 2D plans. The average time for these procedures are 48.02 ±11.55 minutes from the time patient entered the treatment room until s/he exited, and 35.09 ±10.35 minutes from starting CBCT until last beam delivered. This time duration is comparable to the net time when individual tasks are summed up during our regular CT- based whole brain planning and delivery. Conclusions: The OSPD whole brain treatment has been tested to be clinically feasible. The next step is to further improve the efficiency and to streamline the workflow. Other disease sites will be also tested with this new technology.

60 APPLIED LIFE SCIENCES↗