Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parallel 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 1,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Preliminary investigation of a water-based method for fast integrating mobility spectrometry

A water-based condensational growth channel was developed for imaging mobility-separated particles within a parallel plate separation channel of the Fast Integrated Mobility Spectrometer (FIMS). Reported are initial tests of that system, in which the alcohol condenser of the FIMS was replaced by a water-based condensational growth channel. Tests with monodispersed sodium chloride aerosol verify that the water-condensational growth maintained the laminar flow, while providing sufficient growth for particle imaging. Particle positions mapped onto particle mobility, in accordance with theoretical expectations. Particles ranging in size from 12 nm to 100 nm were counted with the same efficiency as with a butanol-based ultrafine particle counter, once inlet and line losses were taken into account.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

Robust and efficient overset grid assembly for partitioned unstructured meshes

This paper presents a method to perform efficient and automated Overset Grid Assembly (OGA) on a system of overlapping unstructured meshes in a parallel computing environment where all meshes are partitioned into multiple mesh-blocks and processed on multiple cores. The main task of the overset grid assembler is to identify, in parallel, among all points in the overlapping mesh system, at which points the flow solution should be computed (field points), interpolated (receptor points), or ignored (hole points). Point containment search or donor search, an algorithm to efficiently determine the cell that contains a given point, is the core procedure necessary for accomplishing this task. Donor search is particularly challenging for partitioned unstructured meshes because of the complex irregular boundaries that are often created during partitioning. Another challenge arises because of the large variation in the type of mesh-block overlap and the resulting large load imbalance on multiple processors. Desirable traits for the grid assembly method are efficiency (requiring only a small fraction of the solver time), robustness (correct identification of all point types), and full automation (no user input required other than the mesh system). Additionally, the method should be scalable, which is an important challenge due to the inherent load imbalance. This paper describes a fully-automated grid assembly method, which can use two different donor search algorithms. One is based on the use of auxiliary grids and Exact Inverse Maps (EIM), and the other is based on the use of Alternating Digital Trees (ADT). The EIM method is demonstrated to be more efficient than the ADT method, while retaining robustness. An adaptive load re-balance algorithm is also designed and implemented, which considerably improves the scalability of the method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

SU-E-I-36: A KWIC and Dirty Look at Dose Savings and Perfusion Metrics in Simulated CT Neuro Perfusion Exams

Purpose: CT neuro perfusion scans are one of the highest dose exams. Methods to reduce dose include decreasing the number of projections acquired per gantry rotation, however conventional reconstruction of such scans leads to sampling artifacts. In this study we investigated a projection view-sharing reconstruction algorithm used in dynamic MRI – “K-space Weighted Image Contrast” (KWIC) – applied to simulated perfusion exams and evaluated dose savings and impacts on perfusion metrics. Methods: A FORBILD head phantom containing simulated time-varying objects was developed and a set of parallel-beam CT projection data was created. The simulated scans were 60 seconds long, 1152 projections per turn, with a rotation time of one second. No noise was simulated. 5mm, 10mm, and 50mm objects were modeled in the brain. A baseline, “full dose” simulation used all projections and reduced dose cases were simulated by downsampling the number of projections per turn from 1152 to 576 (50% dose), 288 (25% dose), and 144 (12.5% dose). KWIC was further evaluated at 72 projections per rotation (6.25%). One image per second was reconstructed using filtered backprojection (FBP) and KWIC. KWIC reconstructions utilized view cores of 36, 72, 144, and 288 views and 16, 8, 4, and 2 subapertures respectively. From the reconstructed images, time-to-peak (TTP), cerebral blood flow (CBF) and the FWHM of the perfusion curve were calculated and compared against reference values from the full-dose FBP data. Results: TTP, CBF, and the FWHM were unaffected by dose reduction (to 12.5%) and reconstruction method, however image quality was improved when using KWIC. Conclusion: This pilot study suggests that KWIC preserves image quality and perfusion metrics when under-sampling projections and that the unique contrast weighting of KWIC could provided substantial dose-savings for perfusion CT scans. Evaluation of KWIC in clinical CT data will be performed in the near future. R01 EB014922, NCI Grant U01 CA181156 (Quantitative Imaging Network), and Tobacco Related Disease Research Project grant 22RT-0131.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SPHERICAL HARMONIC ANALYSES OF INTENSITY MAPPING POWER SPECTRA

Intensity mapping is a promising technique for surveying the large-scale structure of our universe from z = 0 to z ∼ 150, using the brightness temperature field of spectral lines to directly observe previously unexplored portions of our cosmic timeline. Examples of targeted lines include the 21 cm hyperfine transition of neutral hydrogen, rotational lines of carbon monoxide, and fine-structure lines of singly ionized carbon. Recent efforts have focused on detections of the power spectrum of spatial fluctuations, but have been hindered by systematics such as foreground contamination. This has motivated the decomposition of data into Fourier modes perpendicular and parallel to the line of sight, which has been shown to be a particularly powerful way to diagnose systematics. However, such a method is well-defined only in the limit of a narrow-field, flat-sky approximation. This limits the sensitivity of intensity mapping experiments, as it means that wide surveys must be separately analyzed as a patchwork of smaller fields. In this paper, we develop a framework for analyzing intensity mapping data in a spherical Fourier–Bessel basis, which incorporates curved sky effects without difficulty. We use our framework to generalize a number of techniques in intensity mapping data analysis from the flat sky to the curved sky. These include visibility-based estimators for the power spectrum, treatments of interloper lines, and the “foreground wedge” signature of spectrally smooth foregrounds.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Inter-Disciplinary Collaboration in Support of the Post-Standby TREAT Mission

Although analysis methods have advanced significantly in the last two decades, high fidelity multi- physics methods for reactors systems have been under development for only a few years and are not presently mature nor deployed. Furthermore, very few methods provide the ability to simulate rapid transients in three dimensions. Data for validation of advanced time-dependent multi- physics is sparse; at TREAT, historical data were not collected for the purpose of validating three-dimensional methods, let alone multi-physics simulations. Existing data continues to be collected to attempt to simulate the behavior of experiments and calibration transients, but it will be insufficient for the complete validation of analysis methods used for TREAT transient simulations. Hence, a 2018 restart will most likely occur without the direct application of advanced modeling and simulation methods. At present, the current INL modeling and simulation team plans to work with TREAT operations staff in performing reactor simulations with MAMMOTH, in parallel with the software packages currently being used in preparation for core restart (e.g., MCNP5, RELAP5, ABAQUS). The TREAT team has also requested specific measurements to be performed during startup testing, currently scheduled to run from February to August of 2018. These startup measurements will be crucial in validating the new analysis methods in preparation for ultimate application for TREAT operations and experiment design. This document describes the collaboration between modeling and simulation staff and restart, operations, instrumentation and experiment development teams to be able to effectively interact and achieve successful validation work during restart testing.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

GRID-BASED EXPLORATION OF COSMOLOGICAL PARAMETER SPACE WITH SNAKE

We present a fully parallelized grid-based parameter estimation algorithm for investigating multidimensional likelihoods called Snake, and apply it to cosmological parameter estimation. The basic idea is to map out the likelihood grid-cell by grid-cell according to decreasing likelihood, and stop when a certain threshold has been reached. This approach improves vastly on the 'curse of dimensionality' problem plaguing standard grid-based parameter estimation simply by disregarding grid cells with negligible likelihood. The main advantages of this method compared to standard Metropolis-Hastings Markov Chain Monte Carlo methods include (1) trivial extraction of arbitrary conditional distributions; (2) direct access to Bayesian evidences; (3) better sampling of the tails of the distribution; and (4) nearly perfect parallelization scaling. The main disadvantage is, as in the case of brute-force grid-based evaluation, a dependency on the number of parameters, N{sub par}. One of the main goals of the present paper is to determine how large N{sub par} can be, while still maintaining reasonable computational efficiency; we find that N{sub par} = 12 is well within the capabilities of the method. The performance of the code is tested by comparing cosmological parameters estimated using Snake and the WMAP-7 data with those obtained using CosmoMC, the current standard code in the field. We find fully consistent results, with similar computational expenses, but shorter wall time due to the perfect parallelization scheme.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

High Throughput Source-less Plasma Deposition of Structured Silicon Anodes for Lithium-Ion Batteries

Amprius developed a manufacturing solution for silicon nanowire anode that relies on an inexpensive, high throughput, and high gas precursor utilization plasma deposition method that uses the anode foils as electrodes for plasma generation. The capacitively couple plasma (CCP) method is used in semiconductor and photovoltaic industry and Amprius modified existing high throughput equipment to use anode foils and to deposit amorphous silicon. The equipment was installed ahead of the program at Amprius site. The rest of the tasks included foil handling and process development. The equipment passed site acceptance tests (SAT) and the process parameter mapping was completed, indicating that the target process window limits produce output materials at the rate and with yield and specifications that meet the manufacturing target criteria. Amprius has hired supporting personnel to optimize processes and run the equipment. A parallel task verified the baseline performance of the silicon anode material, to be used as reference for the new manufacturing method.

25 ENERGY STORAGE↗
DOE OSTI2019

High Order Anchoring and Reinitialization of Level Set Function for Simulating Interface Motion

A second order interface anchoring method has been developed and used with fast sweeping algorithm for reinitialization of a level set function. The algebraic anchoring formulation ensures that the location of the actual interface is preserved, leading to better mass conservation property. It also provides high order accurate algebraic constraint for solving the Eikonal equation on a finite difference grid. Geometric properties of the interface such as surface normal and curvature are subsequently computed from the reinitialized distance function. Various analytical functions for modeling distortion in level set field are considered and accuracy of reinitialization is evaluated using first and second order anchoring schemes. It is also shown that accurate computation of interface curvature requires a high order anchor in addition to a high order fast sweeping method. Mass conservation property of reinitialization is also analyzed by considering test problems from literature including the classic Rider–Kothe single vortex problem. The formulation is suitable for efficient parallelization for both distributed memory and on-node shared memory parallel systems. In conclusion, scalability and performance of the reinitialization scheme on multiple architectures are demonstrated.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Variance reduction for Fokker–Planck based particle Monte Carlo schemes

Recently, Fokker–Planck based particle Monte Carlo schemes have been proposed and evaluated for simulations of rarefied gas flows [1–3]. In this paper, the variance reduction for particle Monte Carlo simulations based on the Fokker–Planck model is considered. First, deviational based schemes were derived and reviewed, and it is shown that these deviational methods are not appropriate for practical Fokker–Planck based rarefied gas flow simulations. This is due to the fact that the deviational schemes considered in this study lead either to instabilities in the case of two-weight methods or to large statistical errors if the direct sampling method is applied. Motivated by this conclusion, we developed a novel scheme based on correlated stochastic processes. The main idea here is to synthesize an additional stochastic process with a known solution, which is simultaneously solved together with the main one. By correlating the two processes, the statistical errors can dramatically be reduced; especially for low Mach numbers. To assess the methods, homogeneous relaxation, planar Couette and lid-driven cavity flows were considered. For these test cases, it could be demonstrated that variance reduction based on parallel processes is very robust and effective.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

EXAGRAPH: Graph and combinatorial methods for enabling exascale applications

Combinatorial algorithms in general and graph algorithms in particular play a critical enabling role in numerous scientific applications. However, the irregular memory access nature of these algorithms makes them one of the hardest algorithmic kernels to implement on parallel systems. With tens of billions of hardware threads and deep memory hierarchies, the exascale computing systems in particular pose extreme challenges in scaling graph algorithms. The codesign center on combinatorial algorithms, ExaGraph, was established to design and develop methods and techniques for efficient implementation of key combinatorial (graph) algorithms chosen from a diverse set of exascale applications. Algebraic and combinatorial methods have a complementary role in the advancement of computational science and engineering, including playing an enabling role on each other. In this paper, we survey the algorithmic and software development activities performed under the auspices of ExaGraph from both a combinatorial and an algebraic perspective. In particular, we detail our recent efforts in porting the algorithms to manycore accelerator (GPU) architectures. We also provide a brief survey of the applications that have benefited from the scalable implementations of different combinatorial algorithms to enable scientific discovery at scale. We believe that several applications will benefit from the algorithmic and software tools developed by the ExaGraph team.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

Structured building model reduction toward parallel simulation

Building energy model reduction exchanges accuracy for improved simulation speed by reducing the number of dynamical equations. Parallel computing aims to improve simulation times without loss of accuracy but is poorly utilized by contemporary simulators and is inherently limited by inter-processor communication. This paper bridges these disparate techniques to implement efficient parallel building thermal simulation. We begin with a survey of three structured reduction approaches that compares their performance to a leading unstructured method. We then use structured model reduction to find thermal clusters in the building energy model and allocate processing resources. Experimental results demonstrate faster simulation and low error without any interprocessor communication.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

TH-CD-201-03: A Real-Time Method to Simultaneously Measure Linear Energy Transfer and Dose for Proton Therapy Using Organic Scintillators

Purpose: The light generated in organic scintillators depends on both the radiation dose and the linear energy transfer (LET). The LET dependence leads to an under-response of the detector in the Bragg peak of proton beams. This phenomenon, called ionization quenching, must be corrected to obtain accurate dose measurements of proton beams. This work exploits the ionization quenching phenomenon to provide a method of measuring LET and auto correcting quenching. Methods: We exposed simultaneously four different organic scintillators (BCF-12, PMMA, PVT, and LSD; 1mm in diameter) and a plane parallel ionization chamber in passively scattered proton beams to doses between 32 and 43 cGy and fluence averaged LET values from 0.47 to 1.26 keV/µm. The LET values for each irradiation condition were determined using a validated Monte Carlo model of the beam line. We determined the quenching parameter in the Birk’s equation for scintillation in BCF-12 for dose measurements. One set of irradiation conditions was used to correlate the scintillation response ratio to the LET values and plot a scintillation response ratio versus LET calibration curve. Irradiation conditions independent from the calibration ones were used to validate this method. Comparisons to the expected values were made on both the basis of dose and LET. Results: Among all the scintillators investigated, the ratio of PMMA to BCF-12 provided the best correlation to LET values and was used as the LET calibration curve. The expected LET values in the validation set were within 2%±6%, which resulted in dose accuracy of 1.5%±5.8% for the range of LET values investigated in this work. Conclusion: We have demonstrated the feasibility of using the ratio between the light output of two organic scintillators to simultaneously measure LET and dose of therapeutic proton beams. Further studies are needed to verify the response in higher LET values.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Peridynamic Multiscale Finite Element Methods

The problem of computing quantum-accurate design-scale solutions to mechanics problems is rich with applications and serves as the background to modern multiscale science research. The prob- lem can be broken into component problems comprised of communicating across adjacent scales, which when strung together create a pipeline for information to travel from quantum scales to design scales. Traditionally, this involves connections between a) quantum electronic structure calculations and molecular dynamics and between b) molecular dynamics and local partial differ- ential equation models at the design scale. The second step, b), is particularly challenging since the appropriate scales of molecular dynamic and local partial differential equation models do not overlap. The peridynamic model for continuum mechanics provides an advantage in this endeavor, as the basic equations of peridynamics are valid at a wide range of scales limiting from the classical partial differential equation models valid at the design scale to the scale of molecular dynamics. In this work we focus on the development of multiscale finite element methods for the peridynamic model, in an effort to create a mathematically consistent channel for microscale information to travel from the upper limits of the molecular dynamics scale to the design scale. In particular, we first develop a Nonlocal Multiscale Finite Element Method which solves the peridynamic model at multiple scales to include microscale information at the coarse-scale. We then consider a method that solves a fine-scale peridynamic model to build element-support basis functions for a coarse- scale local partial differential equation model, called the Mixed Locality Multiscale Finite Element Method. Given decades of research and development into finite element codes for the local partial differential equation models of continuum mechanics there is a strong desire to couple local and nonlocal models to leverage the speed and state of the art of local models with the flexibility and accuracy of the nonlocal peridynamic model. In the mixed locality method this coupling occurs across scales, so that the nonlocal model can be used to communicate material heterogeneity at scales inappropriate to local partial differential equation models. Additionally, the computational burden of the weak form of the peridynamic model is reduced dramatically by only requiring that the model be solved on local patches of the simulation domain which may be computed in parallel, taking advantage of the heterogeneous nature of next generation computing platforms. Addition- ally, we present a novel Galerkin framework, the 'Ambulant Galerkin Method', which represents a first step towards a unified mathematical analysis of local and nonlocal multiscale finite element methods, and whose future extension will allow the analysis of multiscale finite element methods that mix models across scales under certain assumptions of the consistency of those models.

42 ENGINEERING↗
DOE OSTI2019

Inter-digital capacitive sensor for evaluating cable jacket and insulation aging

An inter-digital capacitive (TDC) sensor has previously been used to measure dielectric properties of cable insulation polymer material when placed in direct contact with the insulation. Often cable insulation is covered by a polymer jacket. The dielectric properties of many cable jacket and insulation polymers are known to change due to thermal and radiation exposure-related damage. These dielectric properties frequently track with other measures of cable aging, such as tensile elongation-at-break and indenter modulus that have been broadly established as cable insulation polymer assessment methods. The external jacket of a cable is likely to have a different permittivity from the underlying insulation, and frequently the jacket material exhibits more severe damage than the insulation material due to environmental exposure. Because the jacket serves primarily to guard the cable during installation, as long as the underlying insulation condition is acceptable, the jacket condition is relatively unimportant in service. As part of a continuing program to develop and evaluate nondestructive examination methods that may be applied to cable condition assessment, a set of tools has been developed including (1) a parallel-plate sensor to directly measure the permittivity spectrum of flat sheet material and (2) an TDC and fixture to measure the effect of cable polymer dielectric property change on the sensor response. The TDC consists of two fork-like electrodes facing each other with the fork tines interspersed and separated by a small gap. The electrodes are printed on one side of a flexible substrate that can be conformed to the surface of a cylindrical cable, with tines parallel to the cable axis. The electrodes are connected to a broad-frequency-spectral impedance meter that senses the capacitance between the narrowly gapped electrode tines. This capacitance is known to vary as a function of the permittivity of any material in close proximity to the electrodes. By finite element modeling (FEM) and experimentation, this study investigates the effect of tine spacing and other design parameters associated with the TDC on the voltage (potential) distribution and electric field depth of penetration. The TDC measurement of an unshielded ethylene-propylene rubber (EPR)-insulated cable is shown to track with the degree of aging and quantities obtained by established methods. For jacketed cable systems, the TDC response is dominated by the jacket but, by analyzing measurements from TDC sensors with different depths-of-field penetration into the cable under test, the influence of the chlorinated polyethylene (CPE) cable jacket material degradation can be separated from an assessment of the cable insulation thereby enabling assessment of the insulation beneath/through the jacket.

Cable NDE, Interdigital Capacitance, Cable test↗
DOE OSTI2014

Static electric dipole polarizabilities of An{sup 5+/6+} and AnO{sub 2}{sup +/2+} (An = U, Np, and Pu) ions

The parallel components of static electric dipole polarizabilities have been calculated for the lowest lying spin-orbit states of the penta- and hexavalent oxidation states of the actinides (An) U, Np, and Pu, in both their atomic and molecular diyl ion forms (An{sup 5+/6+} and AnO{sub 2}{sup +/2+}) using the numerical finite-field technique within a four-component relativistic framework. The four-component Dirac-Hartree-Fock method formed the reference for MP2 and CCSD(T) calculations, while multireference Fock space coupled-cluster (FSCC), intermediate Hamiltonian Fock space coupled-cluster (IH-FSCC) and Kramers restricted configuration interaction (KRCI) methods were used to incorporate additional electron correlation. It is observed that electron correlation has significant (∼5 a.u.{sup 3}) impact upon the parallel component of the polarizabilities of the diyls. To the best of our knowledge, these quantities have not been previously reported and they can serve as reference values in the determination of various electronic and response properties (for example intermolecular forces, optical properties, etc.) relevant to the nuclear fuel cycle and material science applications. The highest quality numbers for the parallel components (α{sub zz}) of the polarizability for the lowest Ω levels corresponding to the ground electronic states are (in a.u.{sup 3}) 44.15 and 41.17 for UO{sub 2}{sup +} and UO{sub 2}{sup 2+}, respectively, 45.64 and 41.42 for NpO{sub 2}{sup +} and NpO{sub 2}{sup 2+}, respectively, and 47.15 for the PuO{sub 2}{sup +} ion.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Efficient reconstruction and validation of heterogeneous microstructures for energy applications

The digital reconstruction of microstructures is necessary for simulations in fields ranging from geology to electrochemistry, but the state-of-the-art digital reconstruction techniques often compromise between resolution and field of view. It is challenging to retain detailed microstructure information in large-scale reconstructions. Here, this study investigates different aspects of the Yeong-Torquato algorithm based on correlation functions to make it more efficient. We achieve this goal by reducing the computational complexity of the chord-length distribution function and the two-point correlation function, applying the random sphere-packing method as the initial condition, and restricting potential voxel swaps to interfaces. In addition, a novel superposition parallel scheme is introduced to aid in searching for potential voxel swaps. The algorithm proposed is validated by comparing the pore-size distributions of reconstructed 3D custom battery electrodes from a sample dataset obtained from transmission X-ray microscopy. From a sample image with 200 x 200 pixels, the code can reconstruct a 300 x 300 x 300 structure in under 22 h and reconstruct a 400 x 400 x 400 structure in 43 h with eight cores.

42 ENGINEERING↗
DOE OSTI2025

Microscale Metal Additive Manufacturing by Solid‐State Impact Bonding of Shaped Thin Films

The deposition of device-grade inorganic materials is one key challenge toward the implementation of additive manufacturing (AM) in microfabrication, and to that end, a broad range of physico-chemical principles has been explored for 3D fabrication with micro- and nanoscale resolution. Yet, for metals, a process that achieves material quality rivalling that of established thin-film deposition methods, and at the same time, has the potential to combine high throughput production with a broad palette of processable materials, is still lacking. Here, the kinetic, solid-state bonding of metal thin films for the additive assembly of high-purity, high-density metals with micrometer-scale precision is introduced. Indirect laser ablation accelerates micrometer-thick gold films to hundreds of meters per second without their heating or ablation. Their subsequent impact on the substrate above a critical velocity forms a permanent, metallic bond in the solid state. Stacked layers are of high density (>99%). By defining thin-film layers with established lithographic methods prior to launch, a variable feature size (2–50 µm), arbitrary shape of bonded layers, and parallel transfer of up to 36 independent film units in a single shot, is demonstrated. Thus, the solid-state kinetic bonding principle as a viable and potentially versatile route for micro-scale AM of metals is established.

3D printing↗