Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computer Programming”

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 577 records · Page 32

TAMM: Tensor algebra for many-body methods

Tensor algebra operations such as contractions in computational chemistry consume a significant fraction of the computing time on large-scale computing platforms. The widespread use of tensor contractions between large multi-dimensional tensors in describing electronic structure theory has motivated the development of multiple tensor algebra frameworks targeting heterogeneous computing platforms. In this paper, we present Tensor Algebra for Many-body Methods (TAMM), a framework for productive and performance-portable development of scalable computational chemistry methods. TAMM decouples the specification of the computation from the execution of these operations on available high-performance computing systems. With this design choice, the scientific application developers (domain scientists) can focus on the algorithmic requirements using the tensor algebra interface provided by TAMM, whereas high-performance computing developers can direct their attention to various optimizations on the underlying constructs, such as efficient data distribution, optimized scheduling algorithms, and efficient use of intra-node resources (e.g., graphics processing units). The modular structure of TAMM allows it to support different hardware architectures and incorporate new algorithmic advances. We describe the TAMM framework and our approach to the sustainable development of scalable ground- and excited-state electronic structure methods. We present case studies highlighting the ease of use, including the performance and productivity gains compared to other frameworks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Measurement of early time outer laser beam reflection inside a cylindrical hohlraum at the National Ignition Facility

In a National Ignition Facility (NIF) cylindrical hohlraum, any specular reflection (“glint”) off the wall from the outer cone (incident angle of 50° and 44° relative to the surface of the wall) laser beam will irradiate the capsule poles. If the glint power is sufficiently large during the picket (early time) of the laser pulse, it may seed high-mode perturbations on the capsule surface that can grow during the implosion. To quantify the glint power on the capsule during the picket by the outer beams, we performed dedicated experiments on NIF using a flat witness foil as a surrogate for the capsule in a half-hohlraum target. We found that the measured glint power is approximately comparable to simulations using a nominal electron conduction flux limiter of f = 0.15, and over an order of magnitude lower than that predicted using f = 0.03 in the wall. Based on our current understanding, we conclude that the glint from the outer beams plays an insignificant role in capsule drive asymmetry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Design and analysis of dudded fuel experiments at the National Ignition Facility

Recent experiments conducted at the National Ignition Facility (NIF) within the past 2 years have achieved the burning plasma state and exceeded the Lawson criterion for the first time in the laboratory. Here, we report on a set of experiments where the deuterium and tritium (DT) ice layers were replaced with dudded tritium, hydrogen, and deuterium (THD) fuel mixtures to remove the influence of alpha-heating on hot spot dynamics. The hot spot compression and yield in the absence of alpha particle self-heating were measured to assess the proximity of NIF implosions toward the ignition cliff. We find that the “burn-off” Lawson parameters χnoα inferred from the THD experiments are in good agreement with the inferences from postshot simulations of the DT-layered implosions. The THD for burning plasma shot N210307 yielded χnoα≈0.88±0.03 while the THD for ignition shot N210808 yielded χnoα≈1.04±0.04. These results also provide important context for the observed variability in the repeat attempts of ignition shot N210808 since implosions on the ignition cliff are expected to exhibit very large variations in the fusion yield from small changes in the initial conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

MHD-induced beta limits in the Large Helical Device

Using the extended-magnetohydrodynamics code, M3D-C1, we perform a systematic numerical study of the effect of externally applied heating on the achievable plasma beta in a ten field-period heliotron. Heat sources of varying intensity are applied to a vacuum magnetic field that is representative of the standard configuration of the Large Helical Device, with R 0 = 3.66 m, where R 0 is the radial position of the magnetic axis in vacuum. As the system is driven to a state that is unstable to low-n magnetohydrodynamic (MHD) modes, nonlinear mode interactions lead to the formation of chaotic magnetic fields. With sufficiently strong heating, a collapse of the electron temperature profile is observed. This demonstrates the necessity of simulating the self-consistent evolution of plasma profiles, without imposing assumptions on the structure of the magnetic field, to accurately determine transport properties in stellarator plasmas. It also highlights the value of these advanced simulation capabilities for accelerating the development of high-performance stellarator operating scenarios.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Spatial description of dislocation nucleation in the shock response of single-crystal aluminum

Nonequilibrium molecular dynamics simulations of shock loaded single-crystal Al in the $\langle$100$\rangle$, $\langle$110$\rangle$, $\langle$111$\rangle$, and $\langle$123$\rangle$ orientations are conducted to study elastic and plastic shockwave formation and details associated with dislocation activity. A computer vision-based approach is implemented to capture the presence of dislocations and describe their spatial characteristics in the zone of nucleation behind the propagating shockwave. The methodology developed relies on the sequences of images extracted during shock loading that show dislocation activity within a cross section of the sample. Results reveal that the spacing between activated slip systems is orientation dependent and exhibits a modest reduction for the $\langle$100$\rangle$ and $\langle$111$\rangle$ orientations as shock pressure increases. Comparisons are made to existing theoretical models. Such relationships between shock pressure and dislocation activity, extracted from molecular dynamics simulations, can be used to inform higher length scale simulations or modeling of dislocation-based plasticity during shock.

36 MATERIALS SCIENCE↗

Evaluating nonlocal heat transport in directly driven chromium spheres using x-ray spectroscopy

We report on experiments investigating heat transport in laser-generated plasmas using directly driven chromium spheres. The spheres are fielded at the OMEGA laser facility and are driven with laser intensities of 5×10 14 Wcm −2 . Plasma conditions in the corona and scattered light are measured experimentally and compared against predictions from two-dimensional (2D) radiation-hydrodynamic simulations using different heat transport models. Spectroscopic analysis of x-ray self-emission is used as an additional diagnostic. X-ray emission is integrated over a large region of the plasma, probing regions that are not observed by localized optical Thomson scattering. In particular, x-ray emission peaks near the plasma critical density, so emission from optically thin lines provides information on plasma conditions where nonlocal transport is most likely to be significant. Three common heat transport models are considered: local transport with flux limiters f = 0.15 and f = 0.03, and the nonlocal Schurtz–Nicolai–Busquet (SNB) model. Consistent with previous work, both the high-flux (f = 0.15) and SNB models show good agreement with experimentally measured plasma conditions in the corona despite overpredicting laser absorption, whereas the low-flux (f = 0.03) model fails to match any experimental data. Conditions inferred from x-ray self-emission line ratios support this conclusion during the period of laser peak power, although synthetic spectra for all models fail to match the experiment during the transient portions of the pulse. For these reasons, the low-flux model is again rejected.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Parallel interior-point solver for block-structured nonlinear programs on SIMD/GPU architectures

Here, we investigate how to port the standard interior-point method to new exascale architectures for block-structured nonlinear programs with state equations. Computationally, we decompose the interior-point algorithm into two successive operations: the evaluation of the derivatives and the solution of the associated Karush-Kuhn-Tucker (KKT) linear system. Our method accelerates both operations using two levels of parallelism. First, we distribute the computations on multiple processes using coarse parallelism. Second, each process uses SIMD/GPU accelerators locally to accelerate the operations using fine-grained parallelism. The KKT system is reduced by eliminating the inequalities and the state variables from the corresponding equations. We demonstrate our method's capability on the supercomputer Polaris, a testbed for the future exascale Aurora system. Each node is equipped with four GPUs, a setup amenable to our two-level approach. Our experiments on the stochastic optimal power flow problem show that the reduction method is 50x faster than the sparse linear solver HSL MA57 running in serial on the CPU, and 6x faster than Pardiso running in parallel on CPU on the same number of processes.

97 MATHEMATICS AND COMPUTING↗

Enhanced Dual Confocal Measurement System

A measurement instrument utilizing dual, chromatic, confocal, distance sensors has been jointly developed by General Atomics (GA) and Sandia National Laboratories (SNL) for thickness and flatness measurement of target components used in dynamic materials properties (DMP) experiments on the SNL Z-Machine (Z). Compared to previous methods used in production of these types of targets, the tool saves time and yields a 4x reduction in thickness uncertainty which is one of the largest sources of error in equation of state (EOS) measurements critical to supporting the NNSA’s Stockpile Stewardship program and computer modeling of high energy density experiments. It has numerous differences from earlier instruments operating on the dual confocal sensor principle to accommodate DMP components including larger lateral travel, longer working distance, ability to measure flatness in addition to thickness, built-in thickness calibration standards for quickly checking calibration before and after each measurement, and streamlined operation. Thickness and flatness of 0.2mm-3.3mm thick sections of diamond machined copper and aluminum can be measured to “sub-micron” accuracy. Sections up to 6mm thick can be measured with as-yet undermined accuracy. Furthermore, samples must have one surface which is flat to within 300µm, lateral dimensions of no more than 50mm x 50mm, and height less than 40mm.

36 MATERIALS SCIENCE↗

Multiscale evolution of charmed particles in a nuclear medium

Parton energy-momentum exchange with the quark gluon plasma (QGP) is a multiscale problem. In this work, we calculate the interaction of charm quarks with the QGP within the higher twist formalism at high virtuality and high energy using the Modular All Twist Transverse-scattering Elastic-drag and Radiation (MATTER) model, while the low-virtuality and high-energy portion is treated via a linearized Boltzmann transport formalism. Coherence effect that reduces the medium-induced emission rate in the MATTER model is also taken into account through a virtuality-dependent qˆ, leaving the simultaneous dependence of qˆ on heavy quark mass and virtuality for future studies. The interplay between these two formalisms is studied phenomenologically and used to produce a first description of the D-meson and charged hadron nuclear modification factor R AA across multiple centralities. As a result, all calculations were carried out utilizing the Jet Energy-loss Tomography with a Statistically and Computationally Advanced Program Envelope framework.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

dxtbx : the diffraction experiment toolbox

Data formats for recording X-ray diffraction data continue to evolve rapidly to accommodate new detector technologies developed in response to more intense light sources. Processing the data from single-crystal X-ray diffraction experiments therefore requires the ability to read, and correctly interpret, image data and metadata from a variety of instruments employing different experimental representations. Tools that have previously been developed to address this problem have been limited either by a lack of extensibility or by inconsistent treatment of image metadata. The dxtbx software package provides a consistent interface to both image data and experimental models, while supporting a completely generic user-extensible approach to reading the data files. The library is written in a mixture of C++ and Python and is distributed as part of the cctbx under an open-source licence at http://cctbx.sourceforge.net.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

On the release of cppxfel for processing X-ray free-electron laser images

As serial femtosecond crystallography expands towards a variety of delivery methods, including chip-based methods, and smaller collected data sets, the requirement to optimize the data analysis to produce maximum structure quality is becoming increasingly pressing. Here cppxfel , a software package primarily written in C++, which showcases several data analysis techniques, is released. This software package presently indexes images using DIALS (diffraction integration for advanced light sources) and performs an initial orientation matrix refinement, followed by post-refinement of individual images against a reference data set. Cppxfel is released with the hope that the unique and useful elements of this package can be repurposed for existing software packages. However, as released, it produces high-quality crystal structures and is therefore likely to be also useful to experienced users of X-ray free-electron laser (XFEL) software who wish to maximize the information extracted from a limited number of XFEL images.

97 MATHEMATICS AND COMPUTING↗

Condor : a simulation tool for flash X-ray imaging

Flash X-ray imaging has the potential to determine structures down to molecular resolution without the need for crystallization. The ability to accurately predict the diffraction signal and to identify the optimal experimental configuration within the limits of the instrument is important for successful data collection. This article introduces Condor , an open-source simulation tool to predict X-ray far-field scattering amplitudes of isolated particles for customized experimental designs and samples, which the user defines by an atomic or a refractive index model. The software enables researchers to test whether their envisaged imaging experiment is feasible, and to optimize critical parameters for reaching the best possible result. It also aims to support researchers who intend to create or advance reconstruction algorithms by simulating realistic test data. Condor is designed to be easy to use and can be either installed as a Python package or used from its web interface (http://lmb.icm.uu.se/condor). X-ray free-electron lasers have high running costs and beam time at these facilities is precious. Data quality can be substantially improved by using simulations to guide the experimental design and simplify data analysis.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

WavePropaGator : interactive framework for X-ray free-electron laser optics design and simulations

This article describes the WavePropaGator (WPG) package, a new interactive software framework for coherent and partially coherent X-ray wavefront propagation simulations. The package has been developed at European XFEL for users at the existing and emerging free-electron laser (FEL) facilities, as well as at the third-generation synchrotron sources and future diffraction-limited storage rings. The WPG addresses the needs of beamline scientists and user groups to facilitate the design, optimization and improvement of X-ray optics to meet their experimental requirements. The package uses the Synchrotron Radiation Workshop (SRW) C/C++ library and its Python binding for numerical wavefront propagation simulations. The framework runs reliably under Linux, Microsoft Windows 7 and Apple Mac OS X and is distributed under an open-source license. The available tools allow for varying source parameters and optics layouts and visualizing the results interactively. The wavefront history structure can be used for tracking changes in every particular wavefront during propagation. The batch propagation mode enables processing of multiple wavefronts in workflow mode. The paper presents a general description of the package and gives some recent application examples, including modeling of full X-ray FEL beamlines and start-to-end simulation of experiments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

CONTIN XPCS : software for inverse transform analysis of X-ray photon correlation spectroscopy dynamics

X-ray photon correlation spectroscopy (XPCS) and dynamic light scattering (DLS) reveal materials dynamics using coherent scattering, with XPCS permitting the investigation of dynamics in a more diverse array of materials than DLS. Heterogeneous dynamics occur in many material systems. The authors' recent work has shown how classic tools employed in the DLS analysis of heterogeneous dynamics can be extended to XPCS, revealing additional information that conventional Kohlrausch exponential fitting obscures. The present work describes the software implementation of inverse transform analysis of XPCS data. This software, called CONTIN XPCS , is an extension of traditional CONTIN analysis and accommodates the various dynamics encountered in equilibrium XPCS measurements.

36 MATERIALS SCIENCE↗

TXM-Sandbox : an open-source software for transmission X-ray microscopy data analysis

A transmission X-ray microscope (TXM) can investigate morphological and chemical information of a tens to hundred micrometre-thick specimen on a length scale of tens to hundreds of nanometres. It has broad applications in material sciences and battery research. TXM data processing is composed of multiple steps. A workflow software has been developed that integrates all the tools required for general TXM data processing and visualization. The software is written in Python and has a graphic user interface in Jupyter Notebook . Users have access to the intermediate analysis results within Jupyter Notebook and have options to insert extra data processing steps in addition to those that are integrated in the software. The software seamlessly integrates ImageJ as its primary image viewer, providing rich image visualization and processing routines. As a guide for users, several TXM specific data analysis issues and examples are also presented.

36 MATERIALS SCIENCE↗

Experimental Characterization Of The Interaction Between Carbon Fiber Composite Prepregs During The Preforming Process

Carbon fiber composites have received growing attention because of their high performance. One economic method to manufacturing the composite parts is the sequence of forming followed by the compression molding process. In this sequence, the preforming procedure forms the prepreg, which is the composite with the uncured resin, to the product geometry while the molding process cures the resin. Slip between different prepreg layers is observed in the preforming step and this paper reports a method to characterize the properties of the interaction between different prepreg layers, which is critical to predictive modeling and design optimization. An experimental setup was established to evaluate the interactions at various industrial production conditions. The experimental results were analyzed for an in-depth understanding about how the temperature, the relative sliding speed, and the fiber orientation affect the tangential interaction between two prepreg layers. The interaction factors measured from these experiments will be implemented in the computational preforming program.

Composite, Preforming, Material characterization, ↗

SU-E-T-100: Designing a QA Tool for Enhance Dynamic Wedges Based On Dynalog Files

Purpose: A robust quality assurance (QA) program for computer controlled enhanced dynamic wedge (EDW) has been designed and tested. Calculations to perform such QA test is based upon the EDW dynamic log files generated during dose delivery. Methods: Varian record and verify system generates dynamic log (dynalog) files during dynamic dose delivery. The system generated dynalog files contain information such as date and time of treatment, energy, monitor units, wedge orientation, and type of treatment. It also contains the expected calculated segmented treatment tables (STT) and the actual delivered STT for the treatment delivery as a verification record. These files can be used to assess the integrity and precision of the treatment plan delivery. The plans were delivered with a 6 MV beam from a Varian linear accelerator. For available EDW angles (10°, 15°, 20°, 25°, 30°, 45°, and 60°) Varian STT values were used to manually calculate monitor units for each segment. It can also be used to calculate the EDW factors. Independent verification of fractional MUs per segment was performed against those generated from dynalog files. The EDW factors used to calculate MUs in TPS were dosimetrically verified in solid water phantom with semiflex chamber on central axis. Results: EDW factors were generated from the STT provided by Varian and verified against practical measurements. The measurements were in agreement of the order of 1 % to the calculated EDW data. Variation between the MUs per segment obtained from dynalog files and those manually calculated was found to be less than 2%. Conclusion: An efficient and easy tool to perform routine QA procedure of EDW is suggested. The method can be easily implemented in any institution without a need for expensive QA equipment. An error of the order of ≥2% can be easily detected.

60 APPLIED LIFE SCIENCES↗

SU-E-T-616: Plan Quality Assessment of Both Treatment Planning System Dose and Measurement-Based 3D Reconstructed Dose in the Patient

Purpose: Systematic radiotherapy plan quality assessment promotes quality improvement. Software tools can perform this analysis by applying site-specific structure dose metrics. The next step is to similarly evaluate the quality of the dose delivery. This study defines metrics for acceptable doses to targets and normal organs for a particular treatment site and scores each plan accordingly. The input can be the TPS or the measurement-based 3D patient dose. From this analysis, one can determine whether the delivered dose distribution to the patient receives a score which is comparable to the TPS plan score, otherwise replanning may be indicated. Methods: Eleven neuroblastoma patient plans were exported from Eclipse to the Quality Reports program. A scoring algorithm defined a score for each normal and target structure based on dose-volume parameters. Each plan was scored by this algorithm and the percentage of total possible points was obtained. Each plan also underwent IMRT QA measurements with a Mapcheck2 or ArcCheck. These measurements were input into the 3DVH program to compute the patient 3D dose distribution which was analyzed using the same scoring algorithm as the TPS plan. Results: The mean quality score for the TPS plans was 75.37% (std dev=14.15%) compared to 71.95% (std dev=13.45%) for the 3DVH dose distribution. For 3/11 plans, the 3DVH-based quality score was higher than the TPS score, by between 0.5 to 8.4 percentage points. Eight/11 plans scores decreased based on IMRT QA measurements by 1.2 to 18.6 points. Conclusion: Software was used to determine the degree to which the plan quality score differed between the TPS and measurement-based dose. Although the delivery score was generally in good agreement with the planned dose score, there were some that improved while there was one plan whose delivered dose quality was significantly less than planned. This methodology helps evaluate both planned and delivered dose quality. Sun Nuclear Corporation has provded a license for the software described.

62 RADIOLOGY AND NUCLEAR MEDICINE↗