Search NASA⌕ Search

SEARCH · Search NASA

Results for “Time parallelization”

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 487 records · Page 27

PREEMPT: Scalable Epidemic Interventions Using Submodular Optimization on Multi-GPU Systems

Preventing and slowing the spread of epidemics is achieved through techniques such as vaccination and social distancing. Given practical limitations on the number of vaccines and cost of administration, optimization becomes a necessity. Previous approaches using mathematical programming methods have shown to be effective but are limited by computational costs. In this work, we make several contributions: First, we present a new approach for intervention via maximizing the influence of vaccinated nodes on the network. We call this method \preempt. Next, we prove submodular properties associated with the objective function of our method so that it aids in construction of an efficient greedy approximation strategy. Consequently, we present a new parallel algorithm based on greedy hill climbing for \preempt, and present an efficient parallel implementation for distributed CPU-GPU heterogeneous platforms. Our results demonstrate that \preempt{} is able to achieve a significant reduction (up to 6.75$\times$) in the percentage of people infected on a city-scale network. We also show strong scaling results of \preempt{} on 128 nodes of the Summit supercomputer. Our parallel implementation is able to significantly reduce time to solution, from hours to minutes on large networks. This work represents a first-of-its-kind effort in parallelizing greedy hill climbing and applying it toward devising effective interventions for epidemics.

Minutoli, Marco↗

Real-time fault detection and transient identification using extended Kalman filters - 254

In this paper, real-time fault detection and transient identification are studied with adaptive estimation methods, in particular parallel processing methods by implementing a bank of extended Kalman filters. Adaptive estimation is combined with probabilistic risk assessment to allow for real-time continuous monitoring of initiating events. The approach provides a probability of an initiating event happening. Real-time monitoring of initiating events can reduce labor demand. Each extended Kalman filter is configured to operate at a single configuration, either at normal operation or a faulted state. Fault detection is based on calculating a likelihood function for each extended Kalman filter to determine the posteriori conditional probability that the system configuration is the true configuration. The transient identification is determined by performing a multivariate statistical analysis on each extended Kalman filter to provide a posterior probability assessment that the reactor states and its dynamics are within the allowable limits. The probability distribution being outside the allowable limits indicate a transient has occurred. The combination of parallel processing methods with the multivariate statistical analysis for transient identification allows for the operator to perform transient identification with having a high level of certainty of the current system configuration. We demonstrate the real-time fault detection and transient identification on a simple reactor example to determine the reactor power, reactor temperature, and coolant temperature in the presence of various faults. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A scalable exponential-DG approach for nonlinear conservation laws: With application to Burger and Euler equations

In this work, we propose an Exponential DG framework for partial differential equations. We decompose 7 governing equations into linear and nonlinear parts to which we apply the discontinuous Galerkin 8 (DG) spatial discretization. In particular, we construct the linear part using Jacobian that effectively 9 capture stiff characteristics in the system. The former is integrated analytically, whereas the latter 10 is approximated. This approach i) is stable with a large Courant number (Cr > 1); ii) supports 11 high-order solutions both in time and space; iii) is computationally favorable compared to IMEX 12 DG methods with no preconditioner; iv) becomes comparable to explicit RKDG methods on uniform 13 mesh and beneficial on non-uniform grid for Euler equations; v) is scalable in a modern massively 14 parallel computing architecture due to its explicit nature of exponential time integrators and com15 pact communication stencil of DG method. Numerical results demonstrate the performance of our 16 proposed methods through various examples. We also discuss the stability and convergence analysis 17 for our exponential DG scheme in the context of Burgers equation.

42 ENGINEERING↗

Large-Scale Welding Process Simulation by GPU Parallelized Computing

The computational design of industrially relevant welded structures is extremely time consuming due to coupled physics and high nonlinearity. Previously, most welding distortion and residual stress simulations have been limited to small coupons and reduced order (from three-dimensional [3D] to two-dimensional [2D]), or inherent strain approximations were used for large structures. In this current study, an explicit finite element code based on a graphics processing unit was utilized to perform 3D transient thermomechanical simulation of structural components during welding. Laser brazing of aluminum alloy panels as representative of automotive manufacturing scenarios was simulated to predict out-of-plane distortion under different clamping conditions. The predicted deformation pattern and magnitude were validated by laser scanning data of physical assemblies. In addition, the code was used to investigate residual stresses developed during multipass arc welding of a nuclear industry pressurizer surge nozzle and subsequent welding repair where a 3D simulation was necessary. Taking the experimental data as reference, the 3D model predicted better residual stress distribution than a typical 2D asymmetrical model. Stress evolution in welding repair was also presented and discussed in this study. Furthermore, the efficient numerical model made it feasible to use integrated computational welding engineering to simulate welding processes for large-scale structures.

97 MATHEMATICS AND COMPUTING↗

Parallel replica dynamics method for bistable stochastic reaction networks: Simulation and sensitivity analysis

Stochastic reaction networks that exhibit bistable behavior are common in systems biology, materials science, and catalysis. Sampling of stationary distributions is crucial for understanding and characterizing the long-time dynamics of bistable stochastic dynamical systems. However, simulations are often hindered by the insufficient sampling of rare transitions between the two metastable regions. In this paper, we apply the parallel replica method for a continuous time Markov chain in order to improve sampling of the stationary distribution in bistable stochastic reaction networks. The proposed method uses parallel computing to accelerate the sampling of rare transitions. Furthermore, it can be combined with the path-space information bounds for parametric sensitivity analysis. With the proposed methodology, we study three bistable biological networks: the Schlögl model, the genetic switch network, and the enzymatic futile cycle network. We demonstrate the algorithmic speedup achieved in these numerical benchmarks. More significant acceleration is expected when multi-core or graphics processing unit computer architectures and programming tools such as CUDA are employed

97 MATHEMATICS AND COMPUTING↗

Optimizing inference of segmentation on high-resolution images in MLExchange

MLExchange is a machine learning (ML) operations platform providing web user-interfaces (UIs) for data visualization and analysis pipelines at synchrotron facilities. Among these UIs is the segmentation app which helps synchrotron users utilize ML algorithms to automatically segment high-resolution scientific images with minimal manual annotation effort. In this work, we share code optimizations that significantly speed up the segmentation inference workflow of large data in short time. By optimizing the sequence of CPU-GPU data transfers and introducing CPU parallelization to key operations, we improve the per-device, per-image frame computational efficiency and observe close to 3×$$\times$$ speedup over the original segmentation inference workflow run time when utilizing a single GPU. Further adaptations enabling multi-GPU inference yield more than 40×$$\times$$ speedup with 100 GPUs compared to the optimized single GPU inference workflow. This acceleration of the segmentation inference workflow will provide MLExchange users with easy access to segmentation results with little wait time.

Lu, Shizhao↗

Core inductive electric field during sawtooth crashes on DIII-D

Sawtooth crashes on tokamak plasmas exhibit relaxation much faster than resistive time scales via a mechanism not fully understood. Using core magnetic measurements from the Radial Interferometer-Polarimeter (RIP) diagnostic on the DIII-D tokamak, Grad–Shafranov equilibria constrained by internal magnetic measurements that have high time resolution (<1μs) can be computed, allowing analysis of how equilibrium parameters such as safety factor q, current density J, and parallel electric field E ∥ , particularly on-axis, evolve. At the sawtooth crash, on-axis safety factor q 0 is observed to rise by 5% but remain below 1 throughout the cycle, and on-axis current density J 0 is observed to drop by 5%. On-axis parallel electric field E ∥ (0) is found to be balanced by ηJ 0 (resistivity times on-axis current density) except during the 200 µs crash period, where E ∥ (0) reaches 22 V m –1 , exceeding ηJ 0 by a factor of more than 2000. These first measurements in tokamak plasmas verify that generalized Ohm's law is not balanced during the crash by resistive effects alone; this is a finding expected due to the relaxation being much faster than resistive timescales. As a result, measurement of the electric field during the tokamak sawtooth serves to illuminate the physical mechanisms at work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

LP-BNN: Ultra-low-Latency BNN Inference with Layer Parallelism

High inference latency seriously limits the deployment of DNNs in real-time domains such as autonomous driving, robotic control, and many others. To address this emerging challenge, researchers have proposed approximate DNNs with reduced precision, e.g., Binarized Neural Networks (BNNs). While BNNs can be built to have little loss in accuracy, latency reduction still has much room for improvement. In this paper, we propose a single-FPGA-based BNN accelerator that achieves microsecond-level ultra-low-latency inference of ImageNet, LP- BNN. We obtain this performance via several design optimizations. First, we optimize the network structure by removing Batch Normalization (BN) functions which leads to significant latency in BNNs without any loss on accuracy. Second, we propose a parameterized architecture which is based on layer parallelism and supports nearly perfect load balancing for various types of BNNs. Third, we fuse all the convolution layers and the first fully connected layer. We process them in parallel through fine- grained inter-layer pipelining. With our proposed accelerator, the inference of binarized AlexNet, VGGNet, and ResNet are completed within 21.5µs, 355µs, and 67.8µs respectively, with no loss in accuracy as compared with other BNN implementations.

Geng, Tong↗

Comparison of Hybrid K-Edge Densitometer (HKED) Performance Operating with the Canberra Lynx MCA and the Canberra ICB-NIM Electronics

From the 1991 until 2008 the Canberra Hybrid K-Edge Densitometer systems were provided with ICB-NIM (Integrated Control Bus – Nuclear Instrument Module) acquisition electronics. Newer electronics modules, such as the Lynx, were not supported under the VMS based operating system. The LYNX module was provided as the standard acquisition electronics following the release of the Windows based CHKED software. This report compares the electronics dead-time, gain shifts, detector resolution and measurement performance of the HKED system operated with the two types of acquisition modules. The comparison was performed using measurements obtained with the ORNL HKED system. The original intent of this study was to take advantage of both the timing and energy outputs from the HPGE detector to acquire data with both sets of electronics in parallel. Although this approach has been applied successfully with other systems, in this case we found the timing output produced a significant amount of noise such that a comparison between the electronics would be invalid. So the comparative measurements were performed sequentially. The ICB-NIM data was acquired over the course of 12 months with 255 measurements while the LYNX data was acquired over a period of 10 months with 75 measurements. To simplify the comparison, all data used in this study was acquired using the Canberra CHKED (V1.0) software package. The performance analysis was based primarily on the peak locations, peak widths and concentration values reported by the CHKED software. The raw spectra from the XRF measurements were also examined to extract additional 109Cd peak location and width data for the hybrid measurements (the standard hybrid report template does not report these values).

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Numerical Prediction of Cyclic Variability in a Spark Ignition Engine Using a Parallel Large Eddy Simulation Approach

Cycle-to-cycle variability (CCV) is detrimental to IC engine operation and can lead to partial burn, misfire, and knock. Predicting CCV numerically is extremely challenging due to two key reasons. First, high-fidelity methods such as large eddy simulation (LES) are required to accurately resolve the in-cylinder turbulent flow field both spatially and temporally. Second, CCV is experienced over long timescales and hence the simulations need to be performed for hundreds of consecutive cycles. Ameen et al. (2017, "Parallel Methodology to Capture Cyclic Variability in Motored Engines,"Int. J. Engine Res., 18(4), pp. 366-377.) developed a parallel perturbation model (PPM) approach to dissociate this long time-scale problem into several shorter time-scale problems. This strategy was demonstrated for motored engine and it was shown that the mean and variance of the in-cylinder flow field was captured reasonably well by this approach. In the present study, this PPM approach is extended to simulate the CCV in a fired port-fuel injected (PFI) spark ignition (SI) engine. The predictions from this approach are also shown to be similar to the consecutive LES cycles. It is shown that the parallel approach is able to predict the coefficient of variation (COV) of the in-cylinder pressure and burn raterelated parameters with sufficient accuracy, and is also able to predict the qualitative trends in CCV with changing operating conditions. It is shown that this new approach is able to give accurate predictions of the CCV in fired engines in less than one-tenth of the time required for the conventional approach of simulating consecutive engine cycles.

Ameen, Muhsin M.↗

Defining Lipidomic Responses to Coronavirus Infection

Highly pathogenic human coronavirus infection can cause a severe atypical, rapid onset pneumonia with a mortality rate of up to 10% for severe acute respiratory syndrome coronavirus 1 (SARS-CoV 1), 35% for Middle East respiratory syndrome coronavirus (MERS-CoV), or 1% for severe acute respiratory syndrome coronavirus 2 (SARS-CoV 2 causative agent of COVID 19). While medical countermeasures successfully controlled the worldwide SARS-CoV epidemic, the MERS-CoV epidemic is still ongoing and continues to be a concern for travelers in the Middle East and the multi-year SARS-CoV 2 pandemic underscore the importance of defining biomarkers that are diagnostic and/or predictive of severe disease outcomes for respiratory viruses. Systems biology approaches provide global snapshots of infection induced changes in host cells/tissues and provide extremely rich datasets for understanding host pathogen interactions. Metabolites, especially lipids, are critical for viral replication but less is understood about infection induced changes to lipids due to limits in lipid species detection and identification. To characterize how individual lipid species and proteins with lipid associated functions contribute to highly pathogenic human coronavirus replication and disease severity, existing datasets were probed to determine cell type specific lipid responses to MERS-CoV infection and verification studies were performed to determine if modification of the host lipid signature (by inhibiting enzymatic functions that produce specific lipid species) can perturb CoV replication in human lung cells. MERS-CoV infects human lung epithelial, endothelial and fibroblast cells. All three cell types were infected with MERS-CoV and samples collected to analyze lipids, proteins, metabolites, and transcripts from 12 to 48 hours post infection. Time matched mock-infected cells were collected in parallel for each cell type. Following sample and statistical analysis, functional enrichment and bioinformatic analysis was performed to identify differentially expressed lipids and proteins. Two lipid species were found to be significantly upregulated following MERS-CoV infection, ceramides, and triacylglycerol both of which are indicative of cells undergoing apoptotic cell death. In contrast, sphingomyelins (lipid molecules that can serve as a precursor for one pathway for ceramide synthesis) had significantly decreased expression in MERS-CoV infected cells. An inhibitor of acid sphingomyelinase (that converts sphingomyelin to ceramides and phosphorylcholine) reduced MERS-CoV replication suggesting that production of ceramides is key for successful viral replication and transmission. Acyl-CoA-synthetase 3 (ACSL3), the only protein whose function is lipid associated and had increased differential expression in our dataset, regulates the synthesis of triacylglycerol (increased expression). Inhibitors that directly block ACSL3 (Triacsin C) but not steps later in the triacylglycerol synthesis pathway (Etomoxir) inhibit MERS-CoV replication suggesting that triacylglycerol production and/or ACSL3 activity is also key for viral replication. ACSL3 expression is also upregulated in lung cancer cells and is predicted to promote continued cell viability which would also enhance viral replication. The differentially expressed lipid and lipid-associated protein species identified in our studies suggest that MERS-CoV infection is activating cellular death pathways to limit the number of cells that become infected but also stimulating the production of lipid-associated enzymes that can prolong host cell viability and the amount of time progeny virions can be produced and released. As the inhibitors that worked against MERS-CoV infection were also efficacious against SARS-CoV 2 infection, countermeasures that target these host pathways may provide novel ways to block highly pathogenic human coronavirus infection and/or prevent severe disease outcomes.

59 BASIC BIOLOGICAL SCIENCES↗

A broadband wavelet implementation for rapid ultrasound pulse-echo time-of-flight measurements

A broadband wavelet approach to ultrasonic pulse-echo time-of-flight measurements is described. The broadband approach significantly reduces the time required for frequency-dependent pulse-echo measurements, enabling studies of dynamic systems ranging from biological systems to solid-state phase transitions. The described broadband approach is demonstrated in parallel with the more traditional frequency stepping approach to perform ultrasound time-of-flight measurements inside a large volume Paris-Edinburgh press in-situ at a synchrotron source. The broadband wavelet data acquisition process was found to be 1-2 orders of magnitude faster than the steppedfrequency approach, with no compromise in data quality or determined results.

47 OTHER INSTRUMENTATION↗

Acceleration and sensitivity analysis of lattice kinetic Monte Carlo simulations using parallel processing and rate constant rescaling

Kinetic Monte Carlo (KMC) simulation provides insights into catalytic reactions unobtainable with either experiments or mean-field microkinetic models. Sensitivity analysis of KMC models assesses the robustness of the predictions to parametric perturbations and identifies rate determining steps in a chemical reaction network. Stiffness in the chemical reaction network, a ubiquitous feature, demands lengthy run times for KMC models and renders efficient sensitivity analysis based on the likelihood ratio method unusable. We address the challenge of efficiently conducting KMC simulations and performing accurate sensitivity analysis in systems with unknown time scales by employing two acceleration techniques: rate constant rescaling and parallel processing. We develop statistical criteria that ensure sufficient sampling of non-equilibrium steady state conditions. Our approach provides the twofold benefit of accelerating the simulation itself and enabling likelihood ratio sensitivity analysis, which provides further speedup relative to finite difference sensitivity analysis. As a result, the likelihood ratio method can be applied to real chemistry. We apply our methodology to the water-gas shift reaction on Pt(111).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Experimental Analysis of File Transfer Rates over Wide-Area Dedicated Connections

File transfers over dedicated connections, supported by large parallel file systems, have become increasingly important in high-performance computing and big data workflows. It remains a challenge to achieve peak rates for such transfers due to the complexities of file I/O, host, and network transport subsystems, and equally importantly, their interactions. We present extensive measurements of disk-to-disk file transfers using Lustre and XFS file systems mounted on multi-core servers over a suite of 10 Gbps emulated connections with 0-366 ms round trip times. Our results indicate that large buffer sizes and many parallel flows do not always guarantee high transfer rates. Furthermore, large variations in the measured rates necessitate repeated measurements to ensure confidence in inferences based on them. We propose a new method to efficiently identify the optimal joint file I/O and network transport parameters using a small number of measurements. We show that for XFS and Lustre with direct I/O, this method identifies configurations achieving 97% of the peak transfer rate while probing only 12% of the parameter space.

Rao, Nageswara S.↗

Progressive Hedging Decomposition for Solutions of Large-Scale Process Family Design Problems

Rapid, wide-scale deployment of green process systems, such as carbon capture or water desalination systems, is essential for combatting climate change. Methods relying on traditional design or modularity fail to capture the benefits of both economies of numbers and economies of scale. We have proposed process family design, which designs a family of processes simultaneously exploiting opportunities for common elements. In previous work, we explored different optimization formulations to solve this problem. In this work, we develop a decomposition approach to tackle larger problems efficiently. We solve a water desalination case study, which is too large to solve within a reasonable timeframe with the discretization formulation. We exploit the block angular structure of the discretization problem to decompose and solve using Progressive Hedging (PH). We use the open-source Python package mpi-sppy to execute PH which allows us to leverage parallelization and a HPC cluster to further improve solution time.

Stinchfield, Georgia↗

Nur77 inhibits oxLDL induced apoptosis of macrophages via the p38 MAPK signaling pathway

The interaction between macrophages and oxLDL plays a crucial role in the initiation and progression of atherosclerosis. As a key initiator in a number of plaque promoting processes, oxLDL induces variable effects such as cell apoptosis or proliferation. Orphan nuclear receptor Nur77 is potently induced in macrophages by diverse stimuli, suggesting that it is of importance in vascular inflammation resulting in atherosclerosis, but whether Nur77 induction is detrimental or protective is unclear. In our study, we explore the role of Nur77 in the regulation of oxLDL-induced macrophage apoptosis and the signaling pathways that are involved. We found that oxLDL induced Nur77 expression in a dose and time dependent fashion, and cell viability was decreased in parallel. To determine whether Nur77 induction contributes to the loss of cell viability or is a protective mechanism, the effect of Nur77 overexpression was examined. Importantly, Nur77 overexpression inhibited the oxLDL-induced decrease of cell viability, inhibited the production of apoptotic bodies and restored DNA synthesis following oxLDL exposure. Furthermore, we found that Nur77 induction is mediated through the p38 MAPK signaling pathway. After pretreatment with SB203580, cell viability was decreased, the expression of CyclinA2 and PCNA was attenuated and the percentage of cell apoptosis was enhanced. Likewise, Nur77 overexpression increased the expression of the cell cycle genes PCNA and p21, and attenuated the increase in caspase-3. On the other hand, knockdown of Nur77 expression by specific siRNA resulted in the increased expression of caspase 3. The results demonstrate that Nur77 is induced by oxLDL via the p38 MAPK signaling pathway, which is involved in the regulation of cell survival. Nur77 enhanced cell survival via suppressing apoptosis, without affecting cell proliferation of activated macrophages, which may be beneficial in patients with atherosclerosis. - Highlights: • oxLDL could induce Nur77 expression. • Nur77 overexpression inhibited oxLDL-induced cell viability, production of apoptotic bodies and restored DNA synthesis. • Cell viability, CyclinA2 and PCNA expression and cell apoptosis were mediated through the p38 MAPK signaling pathway. • Nur77 overexpression mediated the expression of genes PCNA, p21, and caspase-3.

60 APPLIED LIFE SCIENCES↗

Observations of Subduction, Downward Heat Flux and Dense Filament Collapse in the Northern Gulf of Mexico

Submesoscale processes are important contributors to the global heat budget and generally support upward heat transport through restratification. However, in salinity‐stratified regions, such as the northern Gulf of Mexico with its influx of freshwater from the Mississippi‐Atchafalaya river system, temperature can act like a passive tracer and submesoscale processes can contribute to downward heat transport. Oceanic heat content is a factor in many environmental risks the region faces, for example, hurricane intensification, and marine heatwaves. During the 2022 field campaign of the Submesoscales Under Near‐Resonant Inertial Shear Experiment, a sampling plan was developed to study such submesoscale processes in high resolution. Over 31 hr, four assets (two research ships and two remotely controlled boats) drove in parallel across a dense filament, capturing its evolution in time and space. The observations show that surface waters, warmed by daytime solar radiation, were subducted and that the associated overturning circulation transported heat below the surface layer where it was later irreversibly mixed away. The estimated downward heat flux was as strong as the concurrent net air‐sea heat flux into the ocean. The filament was then observed to rapidly collapse which we attribute to boundary layer turbulence and the breakdown of geostrophic balance. The collapsing fronts display behaviors indicative of gravity currents. These observations highlight how in salinity‐stratified regions, frontal dynamics can be associated with downward heat flux and how the submesoscale can play an important role in the oceanic heat budget.

58 GEOSCIENCES↗

Plasmon resonant amplification of hot electron-driven photocatalysis

We discuss plasmon resonant excitation of hot electrons in a metal based photocatalyst in the oxygen evolution half reaction in aqueous solution. In this report, the photocatalyst consists of a 100-nm thick Au film deposited on a corrugated silicon substrate. Here, hot electrons photoexcited in the metal are injected into the solution, ultimately reversing the water oxidation reaction (O 2 +4H + +4e–⇌2H 2 O) and producing a photocurrent. In order to amplify this process, the gold electrode is patterned into a plasmon resonant grating structure with a pitch of 500 nm. The photocurrent (i.e., charge transfer rate) is measured as a function of incident angle using 633 nm wavelength light. We observe peaks in the photocurrent at incident angles of ±9° from normal when the light is polarized parallel to the incident plane (p-polarization) and perpendicular to the lines on the grating. Due to these peaks, we estimate an overall plasmonic gain (or amplification) factor of 2.1× in the charge transfer rate. At these same angles, we also observe sharp dips in the photoreflectance, corresponding to the condition when there is wavevector matching between the incident light and the plasmon mode in the grating. No angle dependence is observed in the photocurrent or photoreflectance when the incident light is polarized perpendicular to the incident plane (s-polarization) and parallel to the lines on the grating. Finite difference time domain simulations also predict sharp dips in the photoreflectance at ±9°, and the electric field intensity profiles show clear excitation of a plasmon-resonant mode when illuminated at those angles with p-polarized light.

42 ENGINEERING↗