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,675 records · Page 93

Method of making a distributed optical fiber sensor having enhanced Rayleigh scattering and enhanced temperature stability, and monitoring systems employing same

A method of making an optical fiber sensor device for distributed sensing includes generating a laser beam comprising a plurality of ultrafast pulses, and focusing the laser beam into a core of an optical fiber to form a nanograting structure within the core, wherein the nanograting structure includes a plurality of spaced nanograting elements each extending substantially parallel to a longitudinal axis of optical fiber. Also, an optical fiber sensor device for distributed sensing includes an optical fiber having a longitudinal axis, a core, and a nanograting structure within the core, wherein the nanograting structure includes a plurality of spaced nanograting elements each extending substantially parallel to the longitudinal axis of the optical fiber. Also, a distributed sensing method and system and an energy production system that employs such an optical fiber sensor device.

Chen, Peng Kevin↗

Efficiency Improvements to the Displacement Based Multilevel Structural Optimization Algorithm

Multilevel Structural Optimization (MSO) continues to be an area of research interest in engineering optimization. In the present project, the weight optimization of beams and trusses using Displacement based Multilevel Structural Optimization (DMSO), a member of the MSO set of methodologies, is investigated. In the DMSO approach, the optimization task is subdivided into a single system and multiple subsystems level optimizations. The system level optimization minimizes the load unbalance resulting from the use of displacement functions to approximate the structural displacements. The function coefficients are then the design variables. Alternately, the system level optimization can be solved using the displacements themselves as design variables, as was shown in previous research. Both approaches ensure that the calculated loads match the applied loads. In the subsystems level, the weight of the structure is minimized using the element dimensions as design variables. The approach is expected to be very efficient for large structures, since parallel computing can be utilized in the different levels of the problem. In this paper, the method is applied to a one-dimensional beam and a large three-dimensional truss. The beam was tested to study possible simplifications to the system level optimization. In previous research, polynomials were used to approximate the global nodal displacements. The number of coefficients of the polynomials equally matched the number of degrees of freedom of the problem. Here it was desired to see if it is possible to only match a subset of the degrees of freedom in the system level. This would lead to a simplification of the system level, with a resulting increase in overall efficiency. However, the methods tested for this type of system level simplification did not yield positive results. The large truss was utilized to test further improvements in the efficiency of DMSO. In previous work, parallel processing was applied to the subsystems level, where the derivative verification feature of the optimizer NPSOL had been utilized in the optimizations. This resulted in large runtimes. In this paper, the optimizations were repeated without using the derivative verification, and the results are compared to those from the previous work. Also, the optimizations were run on both, a network of SUN workstations using the MPICH implementation of the Message Passing Interface (MPI) and on the faster Beowulf cluster at ICASE, NASA Langley Research Center, using the LAM implementation of UP]. The results on both systems were consistent and showed that it is not necessary to verify the derivatives and that this gives a large increase in efficiency of the DMSO algorithm.

Plunkett, C. L.↗

A Parallel Cut-Cell Algorithm for the Free-Boundary Grad--Shafranov Problem

A parallel cut-cell algorithm is described to solve the free-boundary problem of the Grad--Shafranov equation. The algorithm reformulates the free-boundary problem in an irregular bounded domain and its important aspects include a searching algorithm for the magnetic axis and separatrix, a surface integral along the irregular boundary to determine the boundary values, an approach to optimize the coil current based on a targeting plasma shape, Picard iterations with Aitken's acceleration for the resulting nonlinear problem, and a Cartesian grid embedded boundary method to handle the complex geometry. Here the algorithm is implemented in parallel using a standard domain-decomposition approach and a good parallel scaling is observed. Numerical results verify the accuracy and efficiency of the free-boundary Grad--Shafranov solver.

97 MATHEMATICS AND COMPUTING↗

New Computational Methods for the Prediction and Analysis of Helicopter Noise

This paper describes several new methods to predict and analyze rotorcraft noise. These methods are: 1) a combined computational fluid dynamics and Kirchhoff scheme for far-field noise predictions, 2) parallel computer implementation of the Kirchhoff integrations, 3) audio and visual rendering of the computed acoustic predictions over large far-field regions, and 4) acoustic tracebacks to the Kirchhoff surface to pinpoint the sources of the rotor noise. The paper describes each method and presents sample results for three test cases. The first case consists of in-plane high-speed impulsive noise and the other two cases show idealized parallel and oblique blade-vortex interactions. The computed results show good agreement with available experimental data but convey much more information about the far-field noise propagation. When taken together, these new analysis methods exploit the power of new computer technologies and offer the potential to significantly improve our prediction and understanding of rotorcraft noise.

Strawn, Roger C.↗

Performance Improvements for the Griffin Transport Solvers

Griffin is a Multiphysics Object-Oriented Simulation Environment based reactor multiphysics analysis application jointly developed by Idaho National Laboratory and Argonne National Laboratory. Griffin includes a variety of deterministic radiation transport solvers for fixed source, k-eigenvalue, adjoint, and subcritical multiplication, as well as transient solvers for point-kinetics, improved quasi-static, and spatial dynamics. A code assessment performed in FY-20 identified two significant issues with the transport solvers in Griffin: first, the primary heterogeneous SN (discrete ordinates) transport solver based on continuous finite element methods required significant mesh refinement and higher memory usage compared to solvers based on the method of characteristic for equivalent accuracy. Second, the homogeneous PN (spherical harmonics expansion) transport solver did not adequately support polynomial refinement, which is a feature usually required for problems with spatial homogenization and pronounced streaming, typical in fast or gas-cooled reactor systems. To address the first issue, the development effort focused on the more promising discontinuous finite element method (DFEM)-based SN transport solver in Griffin. The addition of an asynchronous parallel transport sweeper and coarse mesh finite difference (CMFD) acceleration have rendered a superior heterogeneous SN transport capability for multiphysics problems that requires far less computing resources in terms of both CPU time and memory usage. This is demonstrated with typical thermal- and fast-spectrum reactor benchmark problems, including 2D Transient Reactor Test, 3D Advanced Burner Test Reactor (ABTR), and 2D and 3D Empire microreactor. For the second issue, the development effort focused on a new transport solver based on the hybrid finite element PN method (HFEM-PN), equivalent to the variational nodal method, as well as a new diffusion solver based on HFEM-Diffusion. This solver is intended for homogenized domains with multiphysics coupling (i.e., supports mesh displacement, seamless temperature feedback, etc.). Initial calculations with the HFEM-Diffusion implementation show very good parallel efficiency for the residual evaluations with the 2D ABTR benchmark. A future development effort will be centered on further improvements to the CMFD, HFEM-PN, and DFEM diffusion solvers to ensure Griffin meets performance and software quality assurance requirements for advanced reactor design and analysis.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

TGE: Machine Learning Based Task Graph Embedding for Large-Scale Topology Mapping

Task mapping is an important problem in parallel and distributed computing. The goal in task mapping is to find an optimal layout of the processes of an application (or a task) onto a given network topology. We target this problem in the context of staging applications. A staging application consists of two or more parallel applications (also referred to as staging tasks) which run concurrently and exchange data over the course of computation. Task mapping becomes a more challenging problem in staging applications, because not only data is exchanged between the staging tasks, but also the processes of a staging task may exchange data with each other. We propose a novel method, called Task Graph Embedding (TGE), that harnesses the observable graph structures of parallel applications and network topologies. TGE employs a machine learning based algorithm to find the best representation of a graph, called an embedding, onto a space in which the task-to-processor mapping problem can be solved. We evaluate and demonstrate the effectiveness of TGE experimentally with the communication patterns extracted from runs of XGC, a large-scale fusion simulation code, on Titan.

Choi, Jong Youl↗

A parallelized elliptic solver for reacting flows

A modified Newton algorithm for the solution of nonlinear elliptic boundary value problems via finite discretization methods is presented. A serial implementation of this algorithm which has recently been applied successfully to the computation of an axisymmetric over-ventilated subsonic laminar methane-air jet diffusion flame is described. Parallel implementation issues and a complexity theory are presented. Included as well are actual performance data for model systems obtained on the Intel Hypercube and a discussion of its implications for modeling realistic systems.

Keyes, David E.↗

HPLC-Parallel accelerator and molecular mass spectrometry analysis of 14 C-labeled amino acids

Accelerator mass spectrometry (AMS) is the method of choice for quantitation of low amounts of 14 C-labeled biomolecules. Despite exquisite sensitivity, an important limitation of AMS is its inability to provide structural information about the analyte. This limitation is not critical when the labeled compounds are well-characterized prior to AMS analysis. However, analyte identity is important in other experiments where, for example, a compound is metabolized and the structures of its metabolites are not known. We previously described a moving wire interface that enables direct AMS measurement of liquid sample in the form of discrete drops or HPLC eluent without the need for individual fraction collection, termed liquid sample-AMS (LS-AMS). Here, we now report the coupling of LS-AMS with a molecular mass spectrometer, providing parallel accelerator and molecular mass spectrometry (PAMMS) detection of analytes separated by liquid chromatography. The repeatability of the method was examined by performing repeated injections of 14 C-labeled tryptophan, and relative standard deviations of the 14 C peak areas were ≤10.57% after applying a normalization factor based on a standard. Five 14 C-labeled amino acids were separated and detected to provide simultaneous quantitative AMS and structural MS data, and AMS results were compared with solid sample-AMS (SS-AMS) data using Bland-Altman plots. To demonstrate the utility of the workflow, yeast cells were grown in a medium with 14 C-labeled tryptophan. The cell extracts were analyzed by PAMMS, and 14 C was detected in tryptophan and its metabolite kynurenine.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Polybrominated diphenyl ether (PBDE) neurotoxicity: a systematic review and meta-analysis of animal evidence

Background: Polybrominated diphenyl ethers (PBDEs) are ubiquitous in the environment, with well-documented and widespread human exposure. A recent systematic review and meta-analysis of human studies found an association between prenatal serum PBDE concentrations and a decrease in the IQ of children. Objective: Our goal was to demonstrate the application of systematic review methods, including meta-analysis (if possible), to animal toxicology data to draw conclusions about hazard. We performed a systematic review of experimental animal studies on the neurotoxicity of PBDEs following developmental exposure. Outcomes assessed included measures related to learning, memory, and attention, which parallel the intelligence-related outcomes evaluated in the systematic review of human studies. Methods: We searched PubMed, Embase, and Toxline for experimental studies in non-human mammals published in English through August 15, 2016. Evaluation of risk of bias and the overall body of evidence followed guidance developed by the National Toxicology Program’s Office of Health Assessment and Translation (OHAT). Results: The animal data on PBDEs and learning, memory, and attention were available for BDEs 47, 99, 153, 203, 206, and 209 and the technical mixture DE-71. These data were diverse and complex, with studies using varying designs and outcomes. There were concerns about risk of bias for many studies, and study reporting of methods and results was often incomplete. A meta-analysis of six studies showed consistent evidence of an increase in latency (effect size of 25.8 [CI, 20.3 to 31.2]) in the last trial of the Morris water maze in PBDE-exposed animals, with low heterogeneity. However, for most endpoints, there were unexplained inconsistencies across studies. There was also no consistent evidence of a dose-response relationship. Conclusions: Our analyses illustrates the application of systematic review methods to a chemical hazard assessment. There is a moderate level of evidence that exposure to BDEs 47, 99, and 209 affects learning. For other PBDEs and other endpoints, the level of evidence was low or very low. The meta-analysis led to stronger conclusions than would have been drawn based on a qualitative review of the evidence. The systematic review also identified factors that increased concerns about risk of bias in individual studies that might be remedied by better study reporting.

Polybrominated Diphenyl Ether↗

Neutron-Induced Fission Cross Section Measurements for Uranium Isotopes and Other Actinides at LANSCE

A well established program of neutron-induced fission cross section measurement at Los Alamos Neutron Science Center (LANSCE) is supporting the Fuel Cycle Research program (FC R&D). The incident neutron energy range spans from sub-thermal up to 200 MeV by combining two LANSCE facilities, the Lujan Center and the Weapons Neutron Research center (WNR). The time-of-flight method is implemented to measure the incident neutron energy. A parallel-plate fission ionization chamber was used as a fission fragment detector. The event rate ratio between the investigated foil and a standard {sup 235}U foil is translated into a fission cross section ratio. Thin actinide targets with deposits of <200 {micro}g/cm{sup 2} on stainless steel backing were loaded into a fission chamber. In addition to previously measured data for {sup 237}Np, {sup 239-242}Pu, {sup 243}Am, new measurements include the recently completed {sup 233,238}U isotopes, {sup 236}U data which is being analyzed, and {sup 234}U data acquired in the 2011-2012 LANSCE run cycle. The new data complete the full suite of Uranium isotopes which were investigated with this experimental approach. When analysis of the new measured data is completed, data will be delivered to evaluators. Having data for multiple Uranium isotopes will support theoretical modeling capabilities and strengthens nuclear data evaluation.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

A patch-based pseudo-CT approach for MRI-only radiotherapy in the pelvis

Purpose: In radiotherapy based only on magnetic resonance imaging (MRI), knowledge about tissue electron densities must be derived from the MRI. This can be achieved by converting the MRI scan to the so-called pseudo-computed tomography (pCT). An obstacle is that the voxel intensities in conventional MRI scans are not uniquely related to electron density. The authors previously demonstrated that a patch-based method could produce accurate pCTs of the brain using conventional T{sub 1}-weighted MRI scans. The method was driven mainly by local patch similarities and relied on simple affine registrations between an atlas database of the co-registered MRI/CT scan pairs and the MRI scan to be converted. In this study, the authors investigate the applicability of the patch-based approach in the pelvis. This region is challenging for a method based on local similarities due to the greater inter-patient variation. The authors benchmark the method against a baseline pCT strategy where all voxels inside the body contour are assigned a water-equivalent bulk density. Furthermore, the authors implement a parallelized approximate patch search strategy to speed up the pCT generation time to a more clinically relevant level. Methods: The data consisted of CT and T{sub 1}-weighted MRI scans of 10 prostate patients. pCTs were generated using an approximate patch search algorithm in a leave-one-out fashion and compared with the CT using frequently described metrics such as the voxel-wise mean absolute error (MAE{sub vox}) and the deviation in water-equivalent path lengths. Furthermore, the dosimetric accuracy was tested for a volumetric modulated arc therapy plan using dose–volume histogram (DVH) point deviations and γ-index analysis. Results: The patch-based approach had an average MAE{sub vox} of 54 HU; median deviations of less than 0.4% in relevant DVH points and a γ-index pass rate of 0.97 using a 1%/1 mm criterion. The patch-based approach showed a significantly better performance than the baseline water pCT in almost all metrics. The approximate patch search strategy was 70x faster than a brute-force search, with an average prediction time of 20.8 min. Conclusions: The authors showed that a patch-based method based on affine registrations and T{sub 1}-weighted MRI could generate accurate pCTs of the pelvis. The main source of differences between pCT and CT was positional changes of air pockets and body outline.

60 APPLIED LIFE SCIENCES↗

Mach 4 and Mach 8 axisymmetric nozzles for a shock tunnel

The performance of two axisymmetric nozzles which were designed to produce uniform, parallel flow with nominal Mach numbers of 4 and 8 is examined. A free-piston-driven shock tube was used to supply the nozzle with high-temperature, high-pressure test gas. The inviscid design procedure treated the nozzle expansion in two stages. Close to the nozzle throat, the nozzle wall was specified as conical and the gas flow was treated as a quasi-one-dimensional chemically-reacting flow. At the end of the conical expansion, the gas was assumed to be calorically perfect, and a contoured wall was designed (using method of characteristics) to convert the source flow into a uniform and parallel flow at the end of the nozzle. Performance was assessed by measuring Pitot pressures across the exit plane of the nozzles and, over the range of operating conditions examined, the nozzles produced satisfactory test flows. However, there were flow disturbances in the Mach 8 nozzle flow that persisted for significant times after flow initiation.

Jacobs, P. A.↗

Mach 4 and Mach 8 axisymmetric nozzles for a high-enthalpy shock tunnel

The performance of two axisymmetric nozzles which were designed to produce uniform, parallel flow with nominal Mach numbers of 4 and 8 is examined. A free-piston-driven shock tube was used to supply the nozzle with high-temperature, high-pressure test gas. The inviscid design procedure treated the nozzle expansion in two stages. Close to the nozzle throat, the nozzle wall was specified as conical and the gas flow was treated as a quasi-one-dimensional chemically-reacting flow. At the end of the conical expansion, the gas was assumed to be calorically perfect, and a contoured wall was designed (using method of characteristics) to convert the source flow into a uniform and parallel flow at the end of the nozzle. Performance was assessed by measuring Pitot pressures across the exit plane of the nozzles and, over the range of operating conditions examined, the nozzles produced satisfactory test flows. However, there were flow disturbances in the Mach 8 nozzle flow that persisted for significant times after flow initiation.

Jacobs, P. A.↗

Concurrent implicit spectral deferred correction scheme for low-Mach number combustion with detailed chemistry

We present a parallel multi-implicit time integration scheme for the advection-diffusion-reaction systems arising from the equations governing low-Mach number combustion with complex chemistry. Here, our strategy employs parallelisation across the method to accelerate the serial Multi-Implicit Spectral Deferred Correction (MISDC) scheme used to couple the advection, diffusion, and reaction processes. In our approach, the diffusion solves and the reaction solves are performed concurrently by different processors. Our analysis shows that the proposed parallel scheme is stable for stiff problems and that the sweeps converge to the fixed-point solution at a faster rate than with serial MISDC. We present numerical examples to demonstrate that the new algorithm is high-order accurate in time, and achieves a parallel speedup compared to serial MISDC.

97 MATHEMATICS AND COMPUTING↗

A Conservative, Scalable, Space-Time Blade Element Rotor Model for Multi-Rotor Vehicles

The development of a parallel blade-element rotor model and its implementation into an adaptive Cartesian method is described. The unsteady version of the rotor model applies a body force to all cells contained in the swept space-time volume at each timestep and special care is taken to maintain axisymmetry on the Cartesian grid. Mesh convergence of rotor thrust and torque is obtained with around 10000 cells in the disk for the steady model. Parallelization is accomplished using OpenMP and the rotor force computation is distributed across all available nodes. Simulations of an isolated XV-15 rotor in hover show good correlation with experimental data and predictions of multi-rotor thrust variation closely match previous high fidelity simulations. The final paper will also include results from the unsteady rotor model and parallel scaling tests.

Conservative↗

A Generalized Stress Inversion Approach With Application to Residual Stress Estimation

In this work develop a generalized stress inversion technique (or the generalized inversion method) capable of recovering stresses in linear elastic bodies subjected to arbitrary cuts. Specifically, given a set of displacement measurements found experimentally from digital image correlation (DIC), we formulate a stress estimation inverse problem as a partial differential equation-constrained optimization problem. We use gradient-based optimization methods, and we accordingly derive the necessary gradient and Hessian information in a matrix-free form to allow for parallel, large-scale operations. By using a combination of finite elements, DIC, and a matrix-free optimization framework, the generalized inversion method can be used on any arbitrary geometry, provided that the DIC camera can view a sufficient part of the surface. We present numerical simulations and experiments, and we demonstrate that the generalized inversion method can be applied to estimate residual stress.

42 ENGINEERING↗

A modified T-value method for selection of strain gages for measuring loads on a low aspect ratio wing

A technique which may be useful for selecting strain gages for use in load equations is described. The technique is an adaptation of the previously used T-value method and is applied to a multispar structure. The technique, called the modified T-value method, is used to reduce the number of strain gages used in a load equation from twelve to two. A parallel reduction is made by calculating relative equation accuracies from three applied load distributions. The equations developed from the modified T-value method proved to be accurate more consistently than the T-value method.

Tang, M. H.↗

Net radiation method for enclosure systems involving partially transparent walls

The net radiation method is developed for analyzing radiation heat transfer in enclosure systems involving partially transparent walls. One such system is an enclosure with windows in it. The conventional net radiation method was developed for enclosures having opaque walls. If a partially transparent wall is present, it will permit radiation to enter and leave the enclosure. The net radiation equations are developed here for gray and semigray enclosures with one or more windows. Another system of interest, such as in a flat plate solar collector, consists of a series of parallel transparent layers. The transmission characteristics of such window systems are obtained by the net radiation method, and the technique appears to be more convenient than the ray tracing method which has been used in the past. Relations are developed for windows consisting of any number of parallel layers having differing absorption coefficients and differing surface reflectivities, and for systems composed of parallel transmitting layers and opaque plates.

Siegel, R.↗