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,765 records · Page 98
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Broadband turbulent spectra in gamma-ray burst light curves

Broadband power density spectra offer a window to understanding turbulent behavior in the emission mechanism and, at the highest frequencies, in the putative inner engines powering long gamma-ray bursts (GRBs). We describe a chirp search method alongside Fourier analysis for signal detection in the Poisson noise-dominated, 2 kHz sampled, BeppoSAX light curves. An efficient numerical implementation is described in O(Nnlog n) operations, where N is the number of chirp templates and n is the length of the light-curve time series, suited for embarrassingly parallel processing. For the detection of individual chirps over a 1 s duration, the method is one order of magnitude more sensitive in signal-to-noise ratio than Fourier analysis. The Fourier-chirp spectra of GRB 010408 and GRB 970816 show a continuation of the spectral slope with up to 1 kHz of turbulence identified in low-frequency Fourier analysis. The same continuation is observed in an average spectrum of 42 bright, long GRBs. An outlook on a similar analysis of upcoming gravitational wave data is included.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Nonlinear Fourier analysis for discontinuous conductivities: Computational results

Two reconstruction methods of Electrical Impedance Tomography (EIT) are numerically compared for nonsmooth conductivities in the plane based on the use of complex geometrical optics (CGO) solutions to D-bar equations involving the global uniqueness proofs for Calderón problem exposed in Nachman (1996) [43] and Astala and Päivärinta (2006) [6]: the Astala–Päivärinta theory-based low-pass transport matrix method implemented in Astala et al. (2011) [3] and the shortcut method which considers ingredients of both theories. The latter method is formally similar to the Nachman theory-based regularized EIT reconstruction algorithm studied in Knudsen et al. (2009) [34] and several references from there. New numerical results are presented using parallel computation with size parameters larger than ever, leading mainly to two conclusions as follows. First, both methods can approximate piecewise constant conductivities better and better as the cutoff frequency increases, and there seems to be a Gibbs-like phenomenon producing ringing artifacts. Second, the transport matrix method loses accuracy away from a (freely chosen) pivot point located outside of the object to be studied, whereas the shortcut method produces reconstructions with more uniform quality.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2000

Parallel Processing in Combustion Analysis

The objective of this research is to demonstrate the application of the Flow-field Dependent Variation (FDV) method to a problem of current interest in supersonic chemical combustion. Due in part to the stiffness of the chemical reactions, the solution of such problems on unstructured three dimensional grids often dictates the use of parallel computers. Preliminary results for the injection of a supersonic hydrogen stream into vitiated air are presented.

Schunk, Richard Gregory↗
DOE OSTI2019

Energy-Screened Many-Body Expansion: A Practical Yet Accurate Fragmentation Method for Quantum Chemistry

We introduce an implementation of the truncated many-body expansion, MBE(n), in which the n-body corrections are screened using the effective fragment potential force field, and only those that exceed a specified energy threshold are computed at a quantum-mechanical level of theory. This energy-screened MBE(n) approach is tested at the n = 3 level for a sequence of water clusters, (H2O) N=6–34 . A threshold of 0.25 kJ/mol eliminates more than 80% of the subsystem electronic structure calculations and is even more efficacious in that respect than is distance-based screening. Even so, the energy-screened MBE(3) method is faithful to a full-system quantum chemistry calculation to within 1–2 kJ/mol/monomer, even in good quality basis sets such as aug-cc-pVTZ. These errors can be reduced by means of a two-layer approach that involves a Hartree–Fock calculation for the entire cluster. Such a correction proves to be necessary in order to obtain accurate relative energies for conformational isomers of (H2O)20, but the cost of a full-system Hartree–Fock calculation remains smaller than the cost of three-body subsystem calculations at correlated levels of theory. At the level of second-order Møller–Plesset perturbation theory (MP2), a screened MBE(3) calculation plus a full-system Hartree–Fock calculation is less expensive than a full-system MP2 calculation starting at N = 12 water molecules. This is true even if all MBE(3) subsystem calculations are performed on a single 40-core compute node, i.e., without significant parallelization. Energy-screened MBE(n) thus provides a fragment-based method that is accurate, stable in large basis sets, and low in cost, even when the latter is measured in aggregate computer time.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

TuckerMPI: A Parallel C++/MPI Software Package for Large-scale Data Compression via the Tucker Tensor Decomposition

With this study, our goal is compression of massive-scale grid-structured data, such as the multi-terabyte output of a high-fidelity computational simulation. For such data sets, we have developed a new software package called TuckerMPI, a parallel C++/MPI software package for compressing distributed data. The approach is based on treating the data as a tensor, i.e., a multidimensional array, and computing its truncated Tucker decomposition, a higher-order analogue to the truncated singular value decomposition of a matrix. The result is a low-rank approximation of the original tensor-structured data. Compression efficiency is achieved by detecting latent global structure within the data, which we contrast to most compression methods that are focused on local structure. In this work, we describe TuckerMPI, our implementation of the truncated Tucker decomposition, including details of the data distribution and in-memory layouts, the parallel and serial implementations of the key kernels, and analysis of the storage, communication, and computational costs. We test the software on 4.5 and 6.7 terabyte data sets distributed across 100 s of nodes (1,000 s of MPI processes), achieving compression ratios between 100 and 200,000×, which equates to 99--99.999% compression (depending on the desired accuracy) in substantially less time than it would take to even read the same dataset from a parallel file system. Moreover, we show that our method also allows for reconstruction of partial or down-sampled data on a single node, without a parallel computer so long as the reconstructed portion is small enough to fit on a single machine, e.g., in the instance of reconstructing/visualizing a single down-sampled time step or computing summary statistics. The code is available at https://gitlab.com/tensors/TuckerMPI.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Discontinuous Galerkin methods for plasma physics in the scrape-off layer of tokamaks

A new parallel discontinuous Galerkin solver, called ArcOn, is developed to describe the intermittent turbulent transport of filamentary blobs in the scrape-off layer (SOL) of fusion plasma. The model is comprised of an elliptic subsystem coupled to two convection-dominated reaction–diffusion–convection equations. Upwinding is used for a class of numerical fluxes developed to accommodate cross product driven convection, and the elliptic solver uses SIPG, NIPG, IIPG, Brezzi, and Bassi–Rebay fluxes to formulate the stiffness matrix. A novel entropy sensor is developed for this system, designed for a space–time varying artificial diffusion/viscosity regularization algorithm. Some numerical experiments are performed to show convergence order on manufactured solutions, regularization of blob/streamer dynamics in the SOL given unstable parameterizations, long-time stability of modon (or dipole drift vortex) solutions arising in simulations of drift-wave turbulence, and finally the formation of edge mode turbulence in the scrape-off layer under turbulent saturation conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Making many-body interactions nearly pairwise additive: The polarized many-body expansion approach

The Many-Body Expansion (MBE) is a useful tool to simulate condensed phase chemical systems, often avoiding the steep computational cost of usual electronic structure methods. However, it often requires higher than 2-body terms to achieve quantitative accuracy. In this work, we propose the Polarized MBE (PolBE) method where each MBE energy contribution is treated as an embedding problem. In each energy term, a smaller fragment is embedded into a larger, polarized environment and only a small region is treated at the high-level of theory using embedded mean-field theory. The role of polarized environment was found to be crucial in providing quantitative accuracy at the 2-body level. PolBE accurately predicts noncovalent interaction energies for a number of systems, including CO 2 , water, and hydrated ion clusters, with a variety of interaction mechanisms, from weak dispersion to strong electrostatics considered in this work. We further demonstrate that the PolBE interaction energy is predominantly pairwise unlike the usual vacuum MBE that requires higher-order terms to achieve similar accuracy. We numerically show that PolBE often performs better than other widely used embedded MBE methods such as the electrostatically embedded MBE. Finally,owing to the lack of expensive diagonalization of Fock matrices and its embarrassingly parallel nature, PolBE is a promising way to access condensed phase systems with hybrid density functionals that are difficult to treat with currently available methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Model Exploration of an Information-Based Healthcare Intervention Using Parallelization and Active Learning

This paper describes the application of a large-scale active learning method to characterize the parameter space of a computational agent-based model developed to investigate the impact of CommunityRx, a clinical information-based health intervention that provides patients with personalized information about local community resources to meet basic and self-care needs. Additionally, the diffusion of information about community resources and their use is modeled via networked interactions and their subsequent effect on agents' use of community resources across an urban population. A random forest model is iteratively fitted to model evaluations to characterize the model parameter space with respect to observed empirical data. We demonstrate the feasibility of using high-performance computing and active learning model exploration techniques to characterize large parameter spaces; by partitioning the parameter space into potentially viable and non-viable regions, we rule out regions of space where simulation output is implausible to observed empirical data. We argue that such methods are necessary to enable model exploration in complex computational models that incorporate increasingly available micro-level behavior data. We provide public access to the model and high-performance computing experimentation code.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

False Negative Rates of a Macrofoam-Swab Sampling Method with Low Surface Concentrations of Two Bacillus anthracis Surrogates via Real-Time PCR

Surface sampling for Bacillus anthracis spores has traditionally relied on detection via bacterial cultivation methods. Although effective, this approach does not provide the level of organism specificity that can be gained through molecular techniques. False negative rates (FNR) and limits of detection (LOD) were determined for two B. anthracis surrogates with modified rapid viability-polymerase chain reaction (mRV-PCR) following macrofoam-swab sampling. This study was conducted in parallel with a previously reported study that analyzed spores using a plate-culture method. B. anthracis Sterne (BAS) or B. atrophaeus Nakamura (BG) spores were deposited onto four surface materials (glass, stainless steel, vinyl tile, and plastic) at nine target concentrations (2 to 500 spores/coupon; 0.078 to 19.375 colony-forming units [CFU] per cm 2 ). Mean FNR values for mRV-PCR analysis ranged from 0 to 0.917 for BAS and 0 to 0.875 for BG and increased as spore concentration decreased (over the concentrations investigated) for each surface material. FNRs based on mRV-PCR data were not statistically different for BAS and BG, but were significantly lower for glass than for vinyl tile. FNRs also tended to be lower for the mRV-PCR method compared to the culture method. The mRV-PCR LOD 95 was lowest for glass (0.429 CFU/cm 2 with BAS and 0.341 CFU/cm 2 with BG) and highest for vinyl tile (0.919 CFU/cm 2 with BAS and 0.917 CFU/cm 2 with BG). These mRV-PCR LOD 95 values were lower than the culture values (BAS: 0.678 to 1.023 CFU/cm 2 and BG: 0.820 to 1.489 CFU/cm 2 ). The FNR and LOD₉₅ values reported in this work provide guidance for environmental sampling of Bacillus spores at low concentrations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

False Negative Rates of a Macrofoam-Swab Sampling Method with Low Surface Concentrations of Two Bacillus anthracis Surrogates via Real-Time PCR

Surface sampling for Bacillus anthracis spores has traditionally relied on detection via bacterial cultivation methods. Although effective, this approach does not provide the level of organism specificity that can be gained through molecular techniques. False negative rates (FNR) and limits of detection (LOD) were determined for two B. anthracis surrogates with modified rapid viability-polymerase chain reaction (mRV-PCR) following macrofoam-swab sampling. This study was conducted in parallel with a previously reported study that analyzed spores using a plate-culture method. B. anthracis Sterne (BAS) or B. atrophaeus Nakamura (BG) spores were deposited onto four surface materials (glass, stainless steel, vinyl tile, and plastic) at nine target concentrations (2 to 500 spores/coupon; 0.078 to 19.375 colony-forming units [CFU] per cm2). Mean FNR values for mRV-PCR analysis ranged from 0 to 0.917 for BAS and 0 to 0.875 for BG and increased as spore concentration decreased (over the concentrations investigated) for each surface material. FNRs based on mRV-PCR data were not statistically different for BAS and BG, but were significantly lower for glass than for vinyl tile. FNRs also tended to be lower for the mRV-PCR method compared to the culture method. The mRV-PCR LOD95 was lowest for glass (0.429 CFU/cm2 with BAS and 0.341 CFU/cm2 with BG) and highest for vinyl tile (0.919 CFU/cm2 with BAS and 0.917 CFU/cm2 with BG). These mRV-PCR LOD95 values were lower than the culture values (BAS: 0.678 to 1.023 CFU/cm2 and BG: 0.820 to 1.489 CFU/cm2). The FNR and LOD95 values reported in this work provide guidance for environmental sampling of Bacillus spores at low concentrations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

On the Convergence of Overlapping Schwarz Decomposition for Nonlinear Optimal Control

Here, we study the convergence properties of an overlapping Schwarz decomposition algorithm for solving nonlinear optimal control problems (OCPs). The algorithm decomposes the time domain into a set of overlapping subdomains, and solves all subproblems defined over subdomains in parallel. The convergence is attained by updating primal-dual information at the boundaries of overlapping subdomains. We show that the algorithm exhibits local linear convergence, and that the convergence rate improves exponentially with the overlap size. We also establish global convergence results for a general quadratic programming, which enables the application of the Schwarz scheme inside second-order optimization algorithms (e.g., sequential quadratic programming). The theoretical foundation of our convergence analysis is a sensitivity result of nonlinear OCPs, which we call "exponential decay of sensitivity" (EDS). Intuitively, EDS states that the impact of perturbations at domain boundaries (i.e., initial and terminal time) on the solution decays exponentially as one moves into the domain. Here, we expand a previous analysis available in the literature by showing that EDS holds for both primal and dual solutions of nonlinear OCPs, under uniform second-order sufficient condition, controllability condition, and boundedness condition. We conduct experiments with a quadrotor motion planning problem and a partial differential equations (PDE) control problem to validate our theory, and show that the approach is significantly more efficient than alternating direction method of multipliers and as efficient as the centralized interior-point solver.

42 ENGINEERING↗
NASA NTRS2012

Fiber Bragg Grating Sensor System for Monitoring Smart Composite Aerospace Structures

Lightweight, electromagnetic interference (EMI) immune, fiber-optic, sensor- based structural health monitoring (SHM) will play an increasing role in aerospace structures ranging from aircraft wings to jet engine vanes. Fiber Bragg Grating (FBG) sensors for SHM include advanced signal processing, system and damage identification, and location and quantification algorithms. Potentially, the solution could be developed into an autonomous onboard system to inspect and perform non-destructive evaluation and SHM. A novel method has been developed to massively multiplex FBG sensors, supported by a parallel processing interrogator, which enables high sampling rates combined with highly distributed sensing (up to 96 sensors per system). The interrogation system comprises several subsystems. A broadband optical source subsystem (BOSS) and routing and interface module (RIM) send light from the interrogation system to a composite embedded FBG sensor matrix, which returns measurand-dependent wavelengths back to the interrogation system for measurement with subpicometer resolution. In particular, the returned wavelengths are channeled by the RIM to a photonic signal processing subsystem based on powerful optical chips, then passed through an optoelectronic interface to an analog post-detection electronics subsystem, digital post-detection electronics subsystem, and finally via a data interface to a computer. A range of composite structures has been fabricated with FBGs embedded. Stress tensile, bending, and dynamic strain tests were performed. The experimental work proved that the FBG sensors have a good level of accuracy in measuring the static response of the tested composite coupons (down to submicrostrain levels), the capability to detect and monitor dynamic loads, and the ability to detect defects in composites by a variety of methods including monitoring the decay time under different dynamic loading conditions. In addition to quasi-static and dynamic load monitoring, the system can capture acoustic emission events that can be a prelude to structural failure, as well as piezoactuator-induced ultrasonic Lamb-waves-based techniques as a basis for damage detection.

Moslehi, Behzad↗
NASA NTRS2013

Time-Accurate Local Time Stepping and High-Order Time CESE Methods for Multi-Dimensional Flows Using Unstructured Meshes

With the wide availability of affordable multiple-core parallel supercomputers, next generation numerical simulations of flow physics are being focused on unsteady computations for problems involving multiple time scales and multiple physics. These simulations require higher solution accuracy than most algorithms and computational fluid dynamics codes currently available. This paper focuses on the developmental effort for high-fidelity multi-dimensional, unstructured-mesh flow solvers using the space-time conservation element, solution element (CESE) framework. Two approaches have been investigated in this research in order to provide high-accuracy, cross-cutting numerical simulations for a variety of flow regimes: 1) time-accurate local time stepping and 2) highorder CESE method. The first approach utilizes consistent numerical formulations in the space-time flux integration to preserve temporal conservation across the cells with different marching time steps. Such approach relieves the stringent time step constraint associated with the smallest time step in the computational domain while preserving temporal accuracy for all the cells. For flows involving multiple scales, both numerical accuracy and efficiency can be significantly enhanced. The second approach extends the current CESE solver to higher-order accuracy. Unlike other existing explicit high-order methods for unstructured meshes, the CESE framework maintains a CFL condition of one for arbitrarily high-order formulations while retaining the same compact stencil as its second-order counterpart. For large-scale unsteady computations, this feature substantially enhances numerical efficiency. Numerical formulations and validations using benchmark problems are discussed in this paper along with realistic examples.

Chang, Chau-Lyan↗
DOE OSTI2013

In utero and early life arsenic exposure in relation to long-term health and disease

Background: There is a growing body of evidence that prenatal and early childhood exposure to arsenic from drinking water can have serious long-term health implications. Objectives: Our goal was to understand the potential long-term health and disease risks associated with in utero and early life exposure to arsenic, as well as to examine parallels between findings from epidemiological studies with those from experimental animal models. Methods: We examined the current literature and identified relevant studies through PubMed by using combinations of the search terms “arsenic”, “in utero”, “transplacental”, “prenatal” and “fetal”. Discussion: Ecological studies have indicated associations between in utero and/or early life exposure to arsenic at high levels and increases in mortality from cancer, cardiovascular disease and respiratory disease. Additional data from epidemiologic studies suggest intermediate effects in early life that are related to risk of these and other outcomes in adulthood. Experimental animal studies largely support studies in humans, with strong evidence of transplacental carcinogenesis, atherosclerosis and respiratory disease, as well as insight into potential underlying mechanisms of arsenic's health effects. Conclusions: As millions worldwide are exposed to arsenic and evidence continues to support a role for in utero arsenic exposure in the development of a range of later life diseases, there is a need for more prospective studies examining arsenic's relation to early indicators of disease and at lower exposure levels. - Highlights: • We review in utero and early-life As exposure impacts on lifelong disease risks. • Evidence indicates that early-life As increases risks of lung disease, cancer and CVD. • Animal work largely parallels human studies and may lead to new research directions. • Prospective studies and individual exposure assessments with biomarkers are needed. • Assessing intermediary endpoints may aid early intervention and establish causality.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

Multi-scale calculation based on dual domain material point method combined with molecular dynamics

This dissertation combines the dual domain material point method (DDMP) with molecular dynamics (MD) in an attempt to create a multi-scale numerical method to simulate materials undergoing large deformations with high strain rates. In these types of problems, the material is often in a thermodynamically non-equilibrium state, and conventional constitutive relations are often not available. In this method, the closure quantities, such as stress, at each material point are calculated from a MD simulation of a group of atoms surrounding the material point. Rather than restricting the multi-scale simulation in a small spatial region, such as phase interfaces, or crack tips, this multi-scale method can be used to consider non-equilibrium thermodynamic e ects in a macroscopic domain. This method takes advantage that the material points only communicate with mesh nodes, not among themselves; therefore MD simulations for material points can be performed independently in parallel. First, using a one-dimensional shock problem as an example, the numerical properties of the original material point method (MPM), the generalized interpolation material point (GIMP) method, the convected particle domain interpolation (CPDI) method, and the DDMP method are investigated. Among these methods, only the DDMP method converges as the number of particles increases, but the large number of particles needed for convergence makes the method very expensive especially in our multi-scale method where we calculate stress in each material point using MD simulation. To improve DDMP, the sub-point method is introduced in this dissertation, which provides high quality numerical solutions with a very small number of particles. The multi-scale method based on DDMP with sub-points is successfully implemented for a one dimensional problem of shock wave propagation in a cerium crystal. The MD simulation to calculate stress in each material point is performed in GPU using CUDA to accelerate the computation. The numerical properties of the multiscale method are investigated as well as the results from this multi-scale calculation are compared of particles needed for convergence makes the method very expensive especially in our multi-scale method where we calculate stress in each material point using MD simulation. To improve DDMP, the sub-point method is introduced in this dissertation, which provides high quality numerical solutions with a very small number of particles. The multi-scale method based on DDMP with sub-points is successfully implemented for a one dimensional problem of shock wave propagation in a cerium crystal. The MD simulation to calculate stress in each material point is performed in GPU using CUDA to accelerate the computation. The numerical properties of the multiscale method are investigated as well as the results from this multi-scale calculation are compared with direct MD simulation results to demonstrate the feasibility of the method. Also, the multi-scale method is applied for a two dimensional problem of jet formation around copper notch under a strong impact.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Micro-lensed fiber optic plate and methods of forming thereof

A method of forming a micro-lensed fiber optic plate. The method includes obtaining a fiber optic plate that has a plurality of fibers arranged in parallel. A first end of a fiber of the plurality of fibers can have a first numerical aperture. The fiber of the plurality of fibers includes a core material running along a length of the fiber that is circumscribed by a layer of cladding material. The method also includes creating a lensed surface at a second end of the fiber of the plurality of fibers. The second end is opposite the first end. The lensed surface of the fiber of the plurality of fibers can have a second numerical aperture that is different from the first numerical aperture.

Bagwell, Brett↗
NASA NTRS1985

Research in computer science

Synopses are given for NASA supported work in computer science at the University of Virginia. Some areas of research include: error seeding as a testing method; knowledge representation for engineering design; analysis of faults in a multi-version software experiment; implementation of a parallel programming environment; two computer graphics systems for visualization of pressure distribution and convective density particles; task decomposition for multiple robot arms; vectorized incomplete conjugate gradient; and iterative methods for solving linear equations on the Flex/32.

Ortega, J. M.↗
NASA NTRS1982

Radiative transfer with partial frequency redistribution in inhomogeneous atmospheres - Application to the Jovian aurora

A direct finite difference numerical solution for the equation of radiative transfer by the Feautrier method is developed for use in planetary atmospheres. The procedure described here uses a plane-parallel atmosphere, and can treat partial frequency redistribution, inhomogeneity, external or internal sources, and various boundary conditions. Isotropic scattering is assumed, but in the case of no frequency redistribution, Rayleigh scattering can also be handled. A program utilizing this method is tested in a variety of situations against more powerful and elaborate methods. The case of the Lyman alpha aurora on Jupiter is then considered, where the effects of partial frequency redistribution are shown to be of great importance. New results for the detailed line profiles for Lyman alpha in the Jovian aurora are presented. The method is quite versatile, and should be especially useful in studying a wide range of problems related to auroral or dayglow emissions in planetary atmospheres.

Gladstone, G. R.↗