Search NASA⌕ Search

SEARCH · Search NASA

Results for “data intensive computing”

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 757 records · Page 42

Convective flow effects on protein crystal growth

The long-term stability of the interferometric setup for the monitoring of protein morphologies has been improved. Growth or dissolution of a crystal on a 100 A scale can now be clearly distinguished from dimensional changes occurring within the optical path of the interferometer. This capability of simultaneously monitoring the local interfacial displacement at several widely-spaced positions on the crystal surface with high local depth resolution, has already yielded novel results. We found with lysozyme that (1) the normal growth rate is oscillatory, and (2) the mean growth step density is greater at the periphery of a facet than in its center. The repartitioning of Na(+) and Cl(-) ions between lysozyme solutions and crystals was studied for a wide range of crystallization conditions. A nucleation-growth-repartitioning model was developed to interpret the large body of data in a unified way. The results strongly suggests that (1) the ion to lysozyme ratio in the crystal depends mostly on kinetic rather than crystallographic parameters, and (2) lysozyme crystals possess a salt-rich core with a diameter on the order of 10 microns. The computational model for diffusive-convective transport in protein crystallization (see the First Report) has been applied to a realistic growth cell geometry, taking into account the findings of the above repartitioning studies. These results show that some elements of a moving boundary problem must be incorporated into the model in order to obtain a more realistic description. Our experimental setup for light scattering investigations of aggregation and nucleation in protein solutions has been extensively tested. Scattering intensity measurements with a true Rayleigh scatterer produced systematically increased forward scattering, indicating problems with glare. These have been resolved. Preliminary measurements with supersaturated lysozyme solutions revealed that the scatterers grow with time. Work has begun on a computer program for the unified evaluation of simultaneously obtained, multi-angle static and dynamic light scattering data.

Rosenberger, Franz↗

Measurements of Neutron Importance Distribution in {sup 235}U Fast Critical Reactor

The reactor physics related integral experiments are significant to criticality safety analysis and the design of the next generation nuclear system. Measurements of parameters such as prompt neutron decay constant α, prompt neutron generation time Λ, effective delayed neutron fraction β{sub eff}, neutron importance distribution, etc., can provide data for the verification of calculation methods, computer codes and cross-section data in criticality safety analysis and test the design of new reactors. Among these parameters, neutron importance distribution has its unique value. Beside verification and validation of nuclear data, computer codes, and simulation methods, the experimental results of neutron importance distribution are also useful to the correcting of the spatial effects in Rossi-α measurement, which is based on point-reactor neutron kinetics. The relative neutron importance distribution experiments were performed on the {sup 235}U fast critical reactor of China Academy of Engineering Physics (CAEP). The key experimental objective was to obtain the importance distribution data for nuclear data validation and criticality safety analysis. We demonstrate in this paper the experiment and the data analysis. {sup 252}Cf neutron importance spatial distributions in the diametral hole of a depleted uranium reflected {sup 235}U fast critical reactor were measured. The evaluated uncertainties of the experiment data were varying with the source location but within 1.5% ranges in all. The experiment results can provide useful data for applications like nuclear data validation and the determination of rectification factor R of Rossi-a measurements. However, the experimental uncertainties are not satisfying and can be further reduced by using a {sup 252}Cf source with much higher intensity, recording more precisely the absolute position of the neutron source, and carrying out the experiment in more shallow subcritical states in which the system multiplication is larger. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Simple Method to Estimate Mean Heart Dose From Hodgkin Lymphoma Radiation Therapy According to Simulation X-Rays

Purpose: To describe a new method to estimate the mean heart dose for Hodgkin lymphoma patients treated several decades ago, using delineation of the heart on radiation therapy simulation X-rays. Mean heart dose is an important predictor for late cardiovascular complications after Hodgkin lymphoma (HL) treatment. For patients treated before the era of computed tomography (CT)-based radiotherapy planning, retrospective estimation of radiation dose to the heart can be labor intensive. Methods and Materials: Patients for whom cardiac radiation doses had previously been estimated by reconstruction of individual treatments on representative CT data sets were selected at random from a case–control study of 5-year Hodgkin lymphoma survivors (n=289). For 42 patients, cardiac contours were outlined on each patient's simulation X-ray by 4 different raters, and the mean heart dose was estimated as the percentage of the cardiac contour within the radiation field multiplied by the prescribed mediastinal dose and divided by a correction factor obtained by comparison with individual CT-based dosimetry. Results: According to the simulation X-ray method, the medians of the mean heart doses obtained from the cardiac contours outlined by the 4 raters were 30 Gy, 30 Gy, 31 Gy, and 31 Gy, respectively, following prescribed mediastinal doses of 25-42 Gy. The absolute-agreement intraclass correlation coefficient was 0.93 (95% confidence interval 0.85-0.97), indicating excellent agreement. Mean heart dose was 30.4 Gy with the simulation X-ray method, versus 30.2 Gy with the representative CT-based dosimetry, and the between-method absolute-agreement intraclass correlation coefficient was 0.87 (95% confidence interval 0.80-0.95), indicating good agreement between the two methods. Conclusion: Estimating mean heart dose from radiation therapy simulation X-rays is reproducible and fast, takes individual anatomy into account, and yields results comparable to the labor-intensive representative CT-based method. This simpler method may produce a meaningful measure of mean heart dose for use in studies of late cardiac complications.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Three-Dimensional Dosimetric Validation of a Magnetic Resonance Guided Intensity Modulated Radiation Therapy System

Purpose: To validate the dosimetric accuracy of a commercially available magnetic resonance guided intensity modulated radiation therapy (MRgIMRT) system using a hybrid approach: 3-dimensional (3D) measurements and Monte Carlo calculations. Methods and Materials: We used PRESAGE radiochromic plastic dosimeters with remote optical computed tomography readout to perform 3D high-resolution measurements, following a novel remote dosimetry protocol. We followed the intensity modulated radiation therapy commissioning recommendations of American Association of Physicists in Medicine Task Group 119, adapted to incorporate 3D data. Preliminary tests (“AP” and “3D-Bands”) were delivered to 9.5-cm usable diameter cylindrical PRESAGE dosimeters to validate the treatment planning system (TPS) for nonmodulated deliveries; assess the sensitivity, uniformity, and rotational symmetry of the PRESAGE dosimeters; and test the robustness of the remote dosimetry protocol. Following this, 4 clinical MRgIMRT plans (“MultiTarget,” “Prostate,” “Head/Neck,” and “C-Shape”) were measured using 13-cm usable diameter PRESAGE dosimeters. For all plans, 3D-γ (3% or 3 mm global, 10% threshold) passing rates were calculated and 3D-γ maps were examined. Point doses were measured with an IBA-CC01 ionization chamber for validation of absolute dose. Finally, by use of an in-house-developed, GPU-accelerated Monte Carlo algorithm (gPENELOPE), we independently calculated dose for all 6 Task Group 119 plans and compared against the TPS. Results: For PRESAGE measurements, 3D-γ analysis yielded passing rates of 98.7%, 99.2%, 98.5%, 98.0%, 99.2%, and 90.7% for AP, 3D-Bands, MultiTarget, Prostate, Head/Neck, and C-Shape, respectively. Ion chamber measurements were within an average of 0.5% (±1.1%) from the TPS dose. Monte Carlo calculations demonstrated good agreement with the TPS, with a mean 3D-γ passing rate of 98.5% ± 1.9% using a stricter 2%/2-mm criterion. Conclusions: We have validated the dosimetric accuracy of a commercial MRgIMRT system using high-resolution 3D techniques. We have demonstrated for the first time that hybrid 3D remote dosimetry is a comprehensive and feasible approach to commissioning MRgIMRT. This may provide better sensitivity in error detection compared with standard 2-dimensional measurements and could be used when implementing complex new magnetic resonance guided radiation therapy technologies.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Optimization of tomographic reconstruction workflows on geographically distributed resources

New technological advancements in synchrotron light sources enable data acquisitions at unprecedented levels. This emergent trend affects not only the size of the generated data but also the need for larger computational resources. Although beamline scientists and users have access to local computational resources, these are typically limited and can result in extended execution times. Applications that are based on iterative processing as in tomographic reconstruction methods require high-performance compute clusters for timely analysis of data. Here, time-sensitive analysis and processing of Advanced Photon Source data on geographically distributed resources are focused on. Two main challenges are considered: (i) modeling of the performance of tomographic reconstruction workflows and (ii) transparent execution of these workflows on distributed resources. For the former, three main stages are considered: (i) data transfer between storage and computational resources, (i) wait/queue time of reconstruction jobs at compute resources, and (iii) computation of reconstruction tasks. These performance models allow evaluation and estimation of the execution time of any given iterative tomographic reconstruction workflow that runs on geographically distributed resources. For the latter challenge, a workflow management system is built, which can automate the execution of workflows and minimize the user interaction with the underlying infrastructure. The system utilizes Globus to perform secure and efficient data transfer operations. The proposed models and the workflow management system are evaluated by using three high-performance computing and two storage resources, all of which are geographically distributed. Workflows were created with different computational requirements using two compute-intensive tomographic reconstruction algorithms. Experimental evaluation shows that the proposed models and system can be used for selecting the optimum resources, which in turn can provide up to 3.13× speedup (on experimented resources). Furthermore, the error rates of the models range between 2.1 and 23.3% (considering workflow execution times), where the accuracy of the model estimations increases with higher computational demands in reconstruction tasks.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Structure, Elastic Constants and XRD Spectra of Extended Solids under High Pressure

We present results of evolutionary simulations based on density functional calculations of a potentially new type of energetic materials called extended solids: P-N and N-H. High-density structures with covalent bonds generated using variable and fixed concentration methods were analysed in terms of thermo-dynamical stability and agreement with experimental X-ray diffraction (XRD) spectra. X-ray diffraction spectra were calculated using a virtual diffraction algorithm that computes kinematic diffraction intensity in three-dimensional reciprocal space before being reduced to a two-theta line profile. Calculated XRD patterns were used to search for the structure of extended solids present at experimental pressures by optimizing data according to experimental XRD peak position, peak intensity and theoretically calculated enthalpy. Elastic constants has been calculated for thermodynamically stable structures of P-N system.

Batyrev, I. G.↗

The measurement of boundary layers on a compressor blade in cascade. Volume 2: Data tables

Measurements were made of the boundary layers and wakes about a highly loaded, double-circular-arc compressor blade in cascade. These laser Doppler velocimetry measurements have yielded a very detailed and precise data base with which to test the application of viscous computational codes to turbomachinery. In order to test the computational codes at off-design conditions, the data have been acquired at a chord Reynolds number of 500,000 and at three incidence angles. Average values and 95 percent confidence bands were tabularized for the velocity, local turbulence intensity, skewness, kurtosis, and percent backflow. Tables also exist for the blade static-pressure distributions and boundary layer velocity profiles reconstructed to account for the normal pressure gradient.

Zierke, William C.↗

Global Bathymetric Prediction For Ocean Modeling and Marine Geophysics

We proposed to construct a complete bathymetric map of the oceans at a 3-10 km resolution by combining all of the available depth soundings collected over the past 30 years with high resolution marine gravity information provided by the Geosat, ERS-1/2, and Topex/Poseidon altimeters. Detailed bathymetry is essential for understanding physical oceanography and marine geophysics. Currents and tides are controlled by the overall shapes of the ocean basins as well as the smaller sharp ocean ridges and seamounts. Because erosion rates are low in the deep oceans, detailed bathymetry reveals the mantle convection patterns, the plate boundaries, the cooling/subsidence of the oceanic lithosphere, the oceanic plateaus, and the distribution of off-ridge volcanoes. We proposed to: (1) Accumulate all available depth soundings collected over the past 30 years; (2) Use the short wavelength (< 160 km) satellite gravity information to interpolate between sparse ship soundings; (3) Improve the resolution of the marine gravity field using enhanced estimates along repeat altimeter profiles together with the dense altimeter measurements; (4) Refine/improve bathymetric predictions using the improved resolution gravity field and also by investigating computer-intensive methods for bathymetric prediction such as inverse theory; and (5) Produce a 'Globe of the Earth' similar to the globe of Venus prepared by the NASA Magellan investigation. This will also include the best available digital land data.

Sandwell, David T.↗

Highly accurate potential energy surface, dipole moment surface, rovibrational energy levels, and infrared line list for {sup 32}S{sup 16}O{sub 2} up to 8000 cm{sup −1}

A purely ab initio potential energy surface (PES) was refined with selected {sup 32}S{sup 16}O{sub 2} HITRAN data. Compared to HITRAN, the root-mean-squares error (σ{sub RMS}) for all J = 0–80 rovibrational energy levels computed on the refined PES (denoted Ames-1) is 0.013 cm{sup −1}. Combined with a CCSD(T)/aug-cc-pV(Q+d)Z dipole moment surface (DMS), an infrared (IR) line list (denoted Ames-296K) has been computed at 296 K and covers up to 8000 cm{sup −1}. Compared to the HITRAN and CDMS databases, the intensity agreement for most vibrational bands is better than 85%–90%. Our predictions for {sup 34}S{sup 16}O{sub 2} band origins, higher energy {sup 32}S{sup 16}O{sub 2} band origins and missing {sup 32}S{sup 16}O{sub 2} IR bands have been verified by most recent experiments and available HITRAN data. We conclude that the Ames-1 PES is able to predict {sup 32/34}S{sup 16}O{sub 2} band origins below 5500 cm{sup −1} with 0.01–0.03 cm{sup −1} uncertainties, and the Ames-296K line list provides continuous, reliable and accurate IR simulations. The K{sub a}-dependence of both line position and line intensity errors is discussed. The line list will greatly facilitate SO{sub 2} IR spectral experimental analysis, as well as elimination of SO{sub 2} lines in high-resolution astronomical observations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Einstein@Home All-Sky Search for Periodic Gravitational Waves in LIGO S5 Data

This paper presents results of an all-sky search for periodic gravitational waves in the frequency range [50, 1 190] Hzand and with frequency derivative range of approx. [-20,1.1] x 10(exp-10) H(sub 0)z s(exp-1) for the fifth LIGO science run (S5). The search uses a noncoherent Hough-transform method to combine the information from coherent searches on time scales of about one day. Because these searches are very computationally intensive, they have been carried out with the Einstein@Home volunteer distributed computing project. Postprocessing identifies eight candidate signals; deeper follow-up studies rule them out. Hence, since no gravitational wave signals have been found, we report upper limits on the intrinsic gravitational wave strain amplitude ho. For example, in the 0.5 Hz-wide band at 152.5 Hz, we can exclude the presence of signals with h(exp 0) greater than 7.6 x 10(exp-25) at a 90% confidence level. This search is about a factor 3 more sensitive than the previous Einstein@Home search of early S5 LIGO data.

J Aasi↗

Insights of Boundary Layer Turbulence Over the Complex Terrain of Central Himalaya from GVAX Field Campaign

Limited observations hinder understanding of turbulent characteristics in mountainous terrain resulting from heating or cooling of slopes, wind, vertical motions, and heat or moisture advection, which disperse aerosols and other pollutants over the region. In this study, the 1290 MHz radar wind profiler data are utilized to compute the boundary layer height (BLH), the refractive index structure constant (C n 2 ), and the energy dissipation rate (ε) over the central Himalayan site for the period of November 2011 to March 2012, from the intense Ganges Valley Aerosol Experiment (GVAX) field measurements. The radar wind profiler (RWP) based estimation of BLH and ε is validated against the radiosonde, representing the effectiveness of the datasets for further investigation. The strong seasonal variation of log C n 2 and log ε, with average values of ≈ -12 m -2/3 and -2 m 2 s -3 , respectively, is associated with the mountain-induced local circulations and stability in the atmospheric boundary layer. The weak stratification during weak flow is found to be responsible for deep mixing, particularly in the nocturnal boundary layer in spring. Furthermore, the level of cloud cover significantly impacts the strength of turbulence, with the highest cloud cover resulting in a substantial increase in log C n 2 (approximately -11 m -2/3 ) due to intense updraft and downdraft motions compared to clear skies. Additionally, the distribution of aerosol loading across the site, coupled with the behavior of BLH, atmospheric stability, and orographic-induced circulations, implies distinctive seasonal mechanisms for transporting aerosols toward the mountains. This study offers valuable insights into the diurnal and seasonal patterns of turbulent mixing and the mechanisms behind the transport of pollutants through boundary layer processes over the region.

54 ENVIRONMENTAL SCIENCES↗

An Image-Guided Study of Setup Reproducibility of Postmastectomy Breast Cancer Patients Treated With Inverse-Planned Intensity Modulated Radiation Therapy

Purpose: To calculate planning target volume (PTV) margins for chest wall and regional nodal targets using daily orthogonal kilovolt (kV) imaging and to study residual setup error after kV alignment using volumetric cone-beam computed tomography (CBCT). Methods and Materials: Twenty-one postmastectomy patients were treated with intensity modulated radiation therapy with 7-mm PTV margins. Population-based PTV margins were calculated from translational shifts after daily kV positioning and/or weekly CBCT data for each of 8 patients, whose surgical clips were used as surrogates for target volumes. Errors from kV and CBCT data were mathematically combined to generate PTV margins for 3 simulated alignment workflows: (1) skin marks alone; (2) weekly kV imaging; and (3) daily kV imaging. Results: The kV data from 613 treatment fractions indicated that a 7-mm uniform margin would account for 95% of daily shifts if patients were positioned using only skin marks. Total setup errors incorporating both kV and CBCT data were larger than those from kV alone, yielding PTV expansions of 7 mm anterior–posterior, 9 mm left–right, and 9 mm superior–inferior. Required PTV margins after weekly kV imaging were similar in magnitude as alignment to skin marks, but rotational adjustments of patients were required in 32% ± 17% of treatments. These rotations would have remained uncorrected without the use of daily kV imaging. Despite the use of daily kV imaging, CBCT data taken at the treatment position indicate that an anisotropic PTV margin of 6 mm anterior–posterior, 4 mm left–right, and 8 mm superior–inferior must be retained to account for residual errors. Conclusions: Cone-beam CT provides additional information on 3-dimensional reproducibility of treatment setup for chest wall targets. Three-dimensional data indicate that a uniform 7-mm PTV margin is insufficient in the absence of daily IGRT. Interfraction movement is greater than suggested by 2-dimensional imaging, thus a margin of at least 4 to 8 mm must be retained despite the use of daily IGRT.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Atmospheric transmittance model for photosynthetically active radiation

A parametric model of the atmospheric transmittance in the PAR band is presented. The model can be straightforwardly applied for calculating the beam, diffuse and global components of the PAR solar irradiance. The required inputs are: air pressure, ozone, water vapor and nitrogen dioxide column content, Ångström's turbidity coefficient and single scattering albedo. Comparison with other models and ground measured data shows a reasonable level of accuracy for this model, making it suitable for practical applications. From the computational point of view the calculus is condensed into simple algebra which is a noticeable advantage. For users interested in speed-intensive computation of the effective PAR solar irradiance, a PC program based on the parametric equations along with a user guide are available online at http://solar.physics.uvt.ro/srms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Uncertainty Quantification of Fatigue Behavior of Rough AM Surfaces and Microstructures to Enable Hydrogen Gas Turbine

Modifying fossil-fueled industrial gas turbines to utilize low or zero-carbon fuels, such as hydrogen or hydrogen-natural gas blends, is a complex endeavor. The successful implementation of this technology hinges on three key design criteria: (1) developing new fuel injectors capable of efficiently burning alternative fuels, (2) ensuring manufacturability to meet cost and time-to-market goals, and (3) achieving component durability in the demanding environment of an operating gas turbine. Additive manufacturing (AM) accelerates product development, yet concerns persist regarding the durability of parts with rough AM surfaces. A fully experimental approach to quantify the fatigue performance of rough AM microstructures is both costly and labor-intensive. To address this, ORNL and Solar Turbines Incorporated (Solar) employed a crystal plasticity finite element (CPFE) model to identify the factors influencing AM surface fatigue behavior. These CPFE findings, combined with targeted experimental data, were used to develop a computationally efficient surrogate model suitable for assessing the lifespan of gas turbine engine components.

08 HYDROGEN↗

Acoustic analysis using numerical solutions of the Navier-Stokes equations

Numerical solutions of the time-dependent compressible Navier-Stokes equations are employed to analyze the unsteady features of jet flows which contribute to jet noise. The turbulent-mass flux spectra are examined in detail in the numerical analysis, and MacCormack's explicit finite difference algorithm is used to solve the governing equations. The computational grid is based on an axisymmetric jet with simple exponential stretching used to optimize the resolution of critical unsteady flow features. The time-varying numerical data are analyzed, and the time averages of the velocity and turbulence-intensity profiles are shown to agree with experimental data. Grid refinement and longer run times improved the prediction of the experimental mass-flux spectra. The numerical procedure can be extended to 3D jet-flow configurations to study the related 3D unsteady flow features.

Scott, James N.↗

The emulsion chamber technology experiment

Photographic emulsion has the unique property of recording tracks of ionizing particles with a spatial precision of 1 micron, while also being capable of deployment over detector areas of square meters or 10's of square meters. Detectors are passive, their cost to fly in Space is a fraction of that of instruments of similar collecting. A major problem in their continued use has been the labor intensiveness of data retrieval by traditional microscope methods. Two factors changing the acceptability of emulsion technology in space are the astronomical costs of flying large electronic instruments such as ionization calorimeters in Space, and the power and low cost of computers, a small revolution in the laboratory microscope data-taking. Our group at UAH made measurements of the high energy composition and spectra of cosmic rays. The Marshall group has also specialized in space radiation dosimetry. Ionization calorimeters, using alternating layers of lead and photographic emulsion, to measure particle energies up to 10(exp 15) eV were developed. Ten balloon flights were performed with them. No such calorimeters have ever flown in orbit. In the ECT program, a small emulsion chamber was developed and will be flown on the Shuttle mission OAST-2 to resolve the principal technological questions concerning space exposures. These include assessments of: (1) pre-flight and orbital exposure to background radiation, including both self-shielding and secondary particle generation; the practical limit to exposure time in space can then be determined; (2) dynamics of stack to optimize design for launch and weightlessness; and (3) thermal and vacuum constraints on emulsion performance. All these effects are cumulative and affect our ability to perform scientific measurements but cannot be adequately predicted by available methods.

Gregory, John C.↗

VLSI Processor For Vector Quantization

Pixel intensities in each kernel compared simultaneously with all code vectors. Prototype high-performance, low-power, very-large-scale integrated (VLSI) circuit designed to perform compression of image data by vector-quantization method. Contains relatively simple analog computational cells operating on direct or buffered outputs of photodetectors grouped into blocks in imaging array, yielding vector-quantization code word for each such block in sequence. Scheme exploits parallel-processing nature of vector-quantization architecture, with consequent increase in speed.

Tawel, Raoul↗

MIRO Computational Model

A computational model calculates the excitation of water rotational levels and emission-line spectra in a cometary coma with applications for the Micro-wave Instrument for Rosetta Orbiter (MIRO). MIRO is a millimeter-submillimeter spectrometer that will be used to study the nature of cometary nuclei, the physical processes of outgassing, and the formation of the head region of a comet (coma). The computational model is a means to interpret the data measured by MIRO. The model is based on the accelerated Monte Carlo method, which performs a random angular, spatial, and frequency sampling of the radiation field to calculate the local average intensity of the field. With the model, the water rotational level populations in the cometary coma and the line profiles for the emission from the water molecules as a function of cometary parameters (such as outgassing rate, gas temperature, and gas and electron density) and observation parameters (such as distance to the comet and beam width) are calculated.

Broderick, Daniel↗