Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Data Intensive Science”

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 937 records · Page 52

Direct-methods structure determination of a trypanosome RNA-editing substrate fragment with translational pseudosymmetry

Using direct methods starting from random phases, the crystal structure of a 32-base-pair RNA (675 non-H RNA atoms in the asymmetric unit) was determined using only the native diffraction data (resolution limit 1.05 Å) and the computer program SIR 2014. The almost three helical turns of the RNA in the asymmetric unit introduced partial or imperfect translational pseudosymmetry (TPS) that modulated the intensities when averaged by the l Miller indices but still escaped automated detection. Almost six times as many random phase sets had to be tested on average to reach a correct structure compared with a similar-sized RNA hairpin (27 nucleotides, 580 non-H RNA atoms) without TPS. Lastly, more sensitive methods are needed for the automated detection of partial TPS.

59 BASIC BIOLOGICAL SCIENCES↗

Untargeted metabolite data from a root surface in a rhizobox

Raw data is provided from samples analyzed using separate reverse phase chromatographic methods on a high performance liquid chromatograph with mass spectrometry. These porewater samples were collected from a microdialysis which generated samples along the surface of a growing A. sative root (all_hc). This data was used to answer questions connecting rhizosphere metabolite (putatively identified metabolites, hc_putative_norm) changes over time (root growth) with changes in the surrounding rhizosphere biogeochemistry (DOC, redox, pH). Rhizosphere biogeochemistry values are provided in the hc_putative_norm file as averages over their respective range of time that they were collected at. The hc_putative_norm file also contains all normalized values over only the intensity values collected for putatively identified metabolites.

54 ENVIRONMENTAL SCIENCES↗

Deep Mapping of Small Solar System Bodies with Galactic Cosmic Ray Secondary Particle Showers

We will investigate the use of galactic cosmic ray (GCR) secondary particles to probe the deep interiors of small solar system bodies (SSBs), including comets, asteroids, and geologic structures on the surfaces of airless bodies. Applications include solar system science, planetary defense, and resource utilization. Our Phase I study demonstrated that muons, the long-range charged component of GCR showers, can penetrate SSBs up to a km in diameter, providing information on their interior structure. Muons produced in Earth’s atmosphere have been applied to image the interior of large objects for science and engineering. In Phase I, we found that the production of muons in the solid surfaces of airless bodies is much smaller than in Earth’s atmosphere. Nevertheless, the flux of transmitted muons is sufficient to detect inclusions within an asteroid or comet in a reasonable amount of time, ranging from hours to weeks, depending on the size of the SSB and the density contrast, position and size of the inclusion. For asteroids and comets, large density variations (e.g., porous soil or ice versus solid rock) are relatively easy to detect. The intrinsic spatial resolution of muon radiography (“muography”) is on the scale of a few meters. The spatial resolution that can be achieved in practice depends on signal intensity and integration time (counting statistics), the angular resolution of the muon tracker (hodoscope) and details of data reduction and analysis methodology. Our Phase II project will assess remaining unknowns for the application of muography to determining the interior structure of SSBs, assess risks for implementation, and provide a roadmap for development of SSB muography beyond the NIAC program. To achieve our objectives, we will focus on four interrelated tasks: Task1) Signal and background characterization: Characterize the production and transmission of muons and secondary particle backgrounds made by cosmic ray showers in SSBs; and near-surface features from radiographic and tomographic data; Task2) Imaging studies: Develop methods to determine the density structure of SSB interiors and near-surface features from radiographic and tomographic data; Task3) Instrument design: Using simulations and bench-top laboratory experiments, investigate specific concepts for the design of compact hodoscopes and components; Task4) Synthesis: Combine the results of the first three tasks to determine the range of applicability of the method, identify the steps needed for maturation of the concept, and explore concepts for a pilot muography mission.

Prettyman, Thomas H.↗

Nuclear Graphite - Fracture Behavior and Modeling

Evidence for the graphite fracture mechanism is reviewed and discussed. The roles of certain microstructural features in the graphite fracture process are reported. The Burchell fracture model is described and its derivation reported. The successful application of the fracture model to uniaxial tensile data from several graphites with widely ranging structure and texture is reported. The extension of the model to multiaxial loading scenarios using two criteria is discussed. Initially, multiaxial strength data for H-451 graphite were modeled using the fracture model and the Principle of Independent Action. The predicted 4th stress quadrant failure envelope was satisfactory but the 1st quadrant predictions were not conservative and thus were unsatisfactory. Multiaxial strength data from the 1st and 4th stress quadrant for NBG-18 graphite are reported. To improve the conservatism of the predicted 1st quadrant failure envelope for NBG-18 the Shetty criterion has been applied to obtain the equivalent critical stress intensity factor, KIc (Equi), for each applied biaxial stress ratio. The equivalent KIc value is used in the Burchell fracture model to predict the failure envelope. The predicted 1st stress quadrant failure envelope is conservative and thus more satisfactory than achieved previously using the fracture model combined with the Principle of Independent Action.

36 MATERIALS SCIENCE↗

Response of flame thickness and propagation speed under intense turbulence in spatially developing lean premixed methane–air jet flames

Direct numerical simulations of three-dimensional spatially-developing turbulent Bunsen flames were performed at three different turbulence intensities. We performed these simulations using a reduced methane–air chemical mechanism which was specifically tailored for the lean premixed conditions simulated here. A planar-jet turbulent Bunsen flame configuration was used in which turbulent preheated methane–air mixture at 0.7 equivalence ratio issued through a central jet and was surrounded by a hot laminar coflow of burned products. The turbulence characteristics at the jet inflow were selected such that combustion occured in the thin reaction zones (TRZ) regime. At the lowest turbulence intensity, the conditions fall on the boundary between the TRZ regime and the corrugated flamelet regime, and progressively moved further into the TRZ regime by increasing the turbulent intensity. The data from the three simulations was analyzed to understand the effect of turbulent stirring on the flame structure and thickness. Furthermore, statistical analysis of the data showed that the thermal preheat layer of the flame was thickened due to the action of turbulence, but the reaction zone was not significantly affected. A global and local analysis of the burning velocity of the flame was performed to compare the different flames. Detailed statistical averages of the flame speed were also obtained to study the spatial dependence of displacement speed and its correlation to strain rate and curvature.

36 MATERIALS SCIENCE↗

Start-to-end simulation of single-particle imaging using ultra-short pulses at the European X-ray Free-Electron Laser

Single-particle imaging with X-ray free-electron lasers (XFELs) has the potential to provide structural information at atomic resolution for non-crystalline biomolecules. This potential exists because ultra-short intense pulses can produce interpretable diffraction data notwithstanding radiation damage. This paper explores the impact of pulse duration on the interpretability of diffraction data using comprehensive and realistic simulations of an imaging experiment at the European X-ray Free-Electron Laser. In conclusion, it is found that the optimal pulse duration for molecules with a few thousand atoms at 5 keV lies between 3 and 9 fs.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

A hybrid reconstruction algorithm for fast and accurate 4D cone-beam CT imaging

Purpose: 4D cone beam CT (4D-CBCT) has been utilized in radiation therapy to provide 4D image guidance in lung and upper abdomen area. However, clinical application of 4D-CBCT is currently limited due to the long scan time and low image quality. The purpose of this paper is to develop a new 4D-CBCT reconstruction method that restores volumetric images based on the 1-min scan data acquired with a standard 3D-CBCT protocol. Methods: The model optimizes a deformation vector field that deforms a patient-specific planning CT (p-CT), so that the calculated 4D-CBCT projections match measurements. A forward-backward splitting (FBS) method is invented to solve the optimization problem. It splits the original problem into two well-studied subproblems, i.e., image reconstruction and deformable image registration. By iteratively solving the two subproblems, FBS gradually yields correct deformation information, while maintaining high image quality. The whole workflow is implemented on a graphic-processing-unit to improve efficiency. Comprehensive evaluations have been conducted on a moving phantom and three real patient cases regarding the accuracy and quality of the reconstructed images, as well as the algorithm robustness and efficiency. Results: The proposed algorithm reconstructs 4D-CBCT images from highly under-sampled projection data acquired with 1-min scans. Regarding the anatomical structure location accuracy, 0.204 mm average differences and 0.484 mm maximum difference are found for the phantom case, and the maximum differences of 0.3–0.5 mm for patients 1–3 are observed. As for the image quality, intensity errors below 5 and 20 HU compared to the planning CT are achieved for the phantom and the patient cases, respectively. Signal-noise-ratio values are improved by 12.74 and 5.12 times compared to results from FDK algorithm using the 1-min data and 4-min data, respectively. The computation time of the algorithm on a NVIDIA GTX590 card is 1–1.5 min per phase. Conclusions: High-quality 4D-CBCT imaging based on the clinically standard 1-min 3D CBCT scanning protocol is feasible via the proposed hybrid reconstruction algorithm.

60 APPLIED LIFE SCIENCES↗

Data from: Soil carbon change in intensive agriculture after 25 years of conservation management

Changes in soil organic carbon (SOC) and nitrogen (SON) are strongly affected by land management, but few long-term comparative studies have surveyed changes throughout the whole soil profile. We quantified 25-year SOC and SON changes to 1 m in 10 replicate ecosystems at an Upper Midwest, USA site. We compared four annual cropping systems in maize ( Zea mays )-soybean ( Glycine max)-winter wheat ( Triticum aestivum ) rotations, each managed differently (Conventional, No-till, Reduced input, and Biologically based); in three managed perennial systems (hybrid Poplar ( Populus × euramericana ), Alfalfa ( Medicago sativa ), and Conifer ( Pinus spp.); and in three successional systems (Early, Mid- and Late succession undergoing a gradual change in species composition and structure over time).

agricultural sciences↗

Assessment of SRS ambient air monitoring network

Three methodologies have been used to assess the effectiveness of the existing ambient air monitoring system in place at the Savannah River Site in Aiken, SC. Effectiveness was measured using two metrics that have been utilized in previous quantification of air-monitoring network performance; frequency of detection (a measurement of how frequently a minimum number of samplers within the network detect an event), and network intensity (a measurement of how consistent each sampler within the network is at detecting events). In addition to determining the effectiveness of the current system, the objective of performing this assessment was to determine what, if any, changes could make the system more effective. Methodologies included 1) the Waite method of determining sampler distribution, 2) the CAP88- PC annual dose model, and 3) a puff/plume transport model used to predict air concentrations at sampler locations. Data collected from air samplers at SRS in 2015 compared with predicted data resulting from the methodologies determined that the frequency of detection for the current system is 79.2% with sampler efficiencies ranging from 5% to 45%, and a mean network intensity of 21.5%. One of the air monitoring stations had an efficiency of less than 10%, and detected releases during just one sampling period of the entire year, adding little to the overall network intensity. By moving or removing this sampler, the mean network intensity increased to about 23%. Further work in increasing the network intensity and simulating accident scenarios to further test the ambient air system at SRS is planned

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Beyond integration: modeling every pixel to obtain better structure factors from stills

Most crystallographic data processing methods use pixel integration. In serial femtosecond crystallography (SFX), the intricate interaction between the reciprocal lattice point and the Ewald sphere is integrated out by averaging symmetrically equivalent observations recorded across a large number (10 4 -10 6 ) of exposures. Although sufficient for generating biological insights, this approach converges slowly, and using it to accurately measure anomalous differences has proved difficult. This report presents a novel approach for increasing the accuracy of structure factors obtained from SFX data. A physical model describing all observed pixels is defined to a degree of complexity such that it can decouple the various contributions to the pixel intensities. Model dependencies include lattice orientation, unit-cell dimensions, mosaic structure, incident photon spectra and structure factor amplitudes. Maximum likelihood estimation is used to optimize all model parameters. The application of prior knowledge that structure factor amplitudes are positive quantities is included in the form of a reparameterization. The method is tested using a synthesized SFX dataset of ytterbium(III) lysozyme, where each X-ray laser pulse energy is centered at 9034 eV. This energy is 100 eV above the Yb 3+ L -III absorption edge, so the anomalous difference signal is stable at 10 electrons despite the inherent energy jitter of each femtosecond X-ray laser pulse. This work demonstrates that this approach allows the determination of anomalous structure factors with very high accuracy while requiring an order-of-magnitude fewer shots than conventional integration-based methods would require to achieve similar results.

59 BASIC BIOLOGICAL SCIENCES↗

Metrics for evaluating tropical cyclones in climate data

This manuscript describes a software suite that can be used for objective evaluation of tropical cyclones (TCs) in gridded climate data. Using cyclone trajectories derived from 6-hourly data, a comprehensive set of metrics is defined to systematically compare and contrast products to one another. In addition to annual TC climatologies, attention is paid to spatial and temporal patterns of storm occurrence and intensity. Assessment can be performed either on the global scale or regional domains. Simple to visualize ‘scorecards’ allow for rapid credibility assessment. We showcase three key findings enabled by this suite. First, we compare the representation of TCs in seven current-generation global reanalyses and conclude that higher resolution models and those with TC-specific assimilation contain more accurate storm climatologies. Second, using a free-running Earth system model (ESM) we find that full basin refinement is required in variable-resolution configurations to adequately simulate North Atlantic TC frequency. Furthermore, upstream refinement over northern Africa offers little benefit in simulating storm occurrence but spatial genesis patterns are improved. Finally, we show that TCs simulated by ESMs can be highly sensitive to individual parameterizations in climate models, with North Atlantic TC metrics varying greatly depending on version of the Morrison-Gettelman microphysics package.

54 ENVIRONMENTAL SCIENCES↗

Building Climate Resilience at NASA Ames

NASA Ames Research Center, located at the southern end of the San Francisco Bay (SFB) estuary has identified three primary vulnerabilities to changes in climate. The Ames Climate Adaptation Science Investigator (CASI) workgroup has studied each of these challenges to operations and the potential exposure of infrastructure and employees to an increased frequency of hazards. Sea level rise inundation scenarios for the SFB Area generally refer to projected scenarios in mean sea level rather than changes in extreme tides that could occur during future storm conditions. In the Summer of 2014, high resolution 3-D mapping of the low lying portion of Ames was performed. Those data are integrated with improved sea level inundation scenarios to identify the buildings, basements and drainage systems potentially affected. We will also identify the impacts of sea level and storm surge effects on transportation to and from the Center. This information will help Center Management develop future Master Plans. Climate change will also lead to changes in temperature, storm frequency and intensity. These changes have potential impacts on localized floods and ecosystems, as well as on electricity and water availability. Over the coming decades, these changes are going to be imposed on top of ongoing land use and land cover changes, especially those deriving from continued urbanization and increase in impervious surface areas. These coupled changes have the potential to create a series of cascading impacts on ecosystems, including changes in primary productivity and disturbance of hydrological properties and increased flood risk.The majority of the electricity used at Ames is supplied by hydroelectric dams, which will be influenced by reductions in precipitation or changes in the timing or phase of precipitation which reduces snow pack. Coupled with increased demand for summertime air conditioning and other cooling needs, NASA Ames is at risk for electricity shortfalls. To assess the anticipated energy usage as climate changes, the Ames CASI team is collecting historical energy usage data from Ames facilities, historical weather data, and projected future weather parameters from the CASI Climate subgroup. This data will be incorporated into the RETScreen model to predict how energy usage at Ames will change over the coming century.

climate↗

Responsivity calibration of the LoWEUS spectrometer

We performed an in situ calibration of the relative responsivity function of the Long-Wavelength Extreme Ultraviolet Spectrometer (LoWEUS), while operating on the Lithium Tokamak Experiment (LTX) at Princeton Plasma Physics Laboratory. The calibration was accomplished by measuring oxygen lines, which are typically present in LTX plasmas. The measured spectral line intensities of each oxygen charge state were then compared to the calculated emission strengths given in the CHIANTI atomic database. Normalizing the strongest line in each charge state to the CHIANTI predictions, we obtained the differences between the measured and predicted values for the relative strengths of the other lines of a given charge state. We find that a 3rd degree polynomial function provides a good fit to the data points. Our measurements show that the responsivity between about 120 and 300 Å varies by factor of ∼30.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

PyDREAM: high-dimensional parameter inference for biological models in python

Biological models contain many parameters whose values are difficult to measure directly via experimentation and therefore require calibration against experimental data. Markov chain Monte Carlo (MCMC) methods are suitable to estimate multivariate posterior model parameter distributions, but these methods may exhibit slow or premature convergence in high-dimensional search spaces. Here, we present PyDREAM, a Python implementation of the (Multiple-Try) Differential Evolution Adaptive Metropolis [DREAM(ZS)] algorithm developed by Vrugt and ter Braak (2008) and Laloy and Vrugt (2012). PyDREAM achieves excellent performance for complex, parameter-rich models and takes full advantage of distributed computing resources, facilitating parameter inference and uncertainty estimation of CPU-intensive biological models.

59 BASIC BIOLOGICAL SCIENCES↗

DoCeph: DPU-Offloaded Messaging in Ceph for Reduced Host CPU Utilization

Ceph is a widely used distributed object store, but its messenger layer imposes substantial CPU overhead on the host. To address this limitation, we propose DoCeph, a DPU-offloaded storage architecture for Ceph that disaggregates the system by offloading the communication-intensive messaging component to the DPU while retaining the storage backend on the host. The DPU efficiently manages communication, using lightweight RPC for metadata operations and DMA for data transfer. Moreover, DoCeph introduces a pipelining technique that overlaps data transmission with buffer preparation, mitigating hardware-imposed transfer size limitations. We implemented DoCeph on a Ceph cluster with NVIDIA BlueField-3 DPUs. Evaluation results indicate that DoCeph cuts host CPU usage by up to 92% while sustaining stable throughput and providing larger performance benefits for object writes over 1 MB.

Park, Kuri [Sogang University]↗

Modeling Low Intensity Fires: Lessons Learned from 2012 RxCADRE

Coupled fire-atmosphere models are increasingly being used to study low-intensity fires, such as those that are used in prescribed fire applications. Thus, the need arises to evaluate these models for their ability to accurately represent fire spread in marginal burning conditions. In this study, wind and fuel data collected during the Prescribed Fire Combustion and Atmospheric Dynamics Research Experiments (RxCADRE) fire campaign were used to generate initial and boundary conditions for coupled fire-atmosphere simulations. We present a novel method to obtain fuels representation at the model grid scale using a combination of imagery, machine learning, and field sampling. Several methods to generate wind input conditions for the model from eight different anemometer measurements are explored. We find a strong sensitivity of fire outcomes to wind inputs. This result highlights the critical need to include variable wind fields as inputs in modeling marginal fire conditions. This work highlights the complexities of comparing physics-based model results against observations, which are more acute in marginal burning conditions, where stronger sensitivities to local variability in wind and fuels drive fire outcomes.

54 ENVIRONMENTAL SCIENCES↗

Computational approaches to model X-ray photon correlation spectroscopy from molecular dynamics

X-ray photon correlation spectroscopy (XPCS) allows for the resolution of dynamic processes within a material across a wide range of length and time scales. X-ray speckle visibility spectroscopy (XSVS) is a related method that uses a single diffraction pattern to probe ultrafast dynamics. Interpretation of the XPCS and XSVS data in terms of underlying physical processes is necessary to establish the connection between the macroscopic responses and the microstructural dynamics. To aid the interpretation of the XPCS and XSVS data, we present a computational framework to model these experiments by computing the X-ray scattering intensity directly from the atomic positions obtained from molecular dynamics simulations. We compare the efficiency and accuracy of two alternative computational methods: the direct method computing the intensity at each diffraction vector separately, and a method based on fast Fourier transform that computes the intensities at all diffraction vectors at once. Furthermore, the computed X-ray speckle patterns capture the density fluctuations over a range of length and time scales and are shown to reproduce the known properties and relations of experimental XPCS and XSVS for liquids.

36 MATERIALS SCIENCE↗

Reconstructing large-scale structure with neutral hydrogen surveys

Upcoming 21-cm intensity surveys will use the hyperfine transition in emission to map out neutral hydrogen in large volumes of the universe. Unfortunately, large spatial scales are completely contaminated with spectrally smooth astrophysical foregrounds which are orders of magnitude brighter than the signal. This contamination also leaks into smaller radial and angular modes to form a foreground wedge, further limiting the usefulness of 21-cm observations for different science cases, especially cross-correlations with tracers that have wide kernels in the radial direction. Here in this paper, we investigate reconstructing these modes within a forward modeling framework. Starting with an initial density field, a suitable bias parameterization and non-linear dynamics to model the observed 21-cm field, our reconstruction proceeds by {combining} the likelihood of a forward simulation to match the observations (under given modeling error and a data noise model) {with the Gaussian prior on initial conditions and maximizing the obtained posterior}. For redshifts z=2 and 4, we are able to reconstruct 21cm field with cross correlation, r c > 0.8 on all scales for both our optimistic and pessimistic assumptions about foreground contamination and for different levels of thermal noise. The performance deteriorates slightly at z=6. The large-scale line-of-sight modes are reconstructed almost perfectly. We demonstrate how our method also provides a technique for density field reconstruction for baryon acoustic oscillations, outperforming standard methods on all scales. We also describe how our reconstructed field can provide superb clustering redshift estimation at high redshifts, where it is otherwise extremely difficult to obtain dense spectroscopic samples, as well as open up a wealth of cross-correlation opportunities with projected fields (e.g. lensing) which are restricted to modes transverse to the line of sight.

21 cm↗