Search NASA⌕ Search

SEARCH · Search NASA

Results for “Parallelization Strategy”

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 955 records · Page 53

Exact signed distance fields using parallel Fast Sweeping Method

Signed distance fields are often used in multiphysics simulations to track material interfaces. We present a simple methodology based on the fast sweeping method to generate the exact signed distance from triangular meshes and linear paths on Cartesian grids. The methodology propagates the closest primitive to the boundary to the rest of the domain following the characteristics. A local upwind criterion is used to decide between the new and existing closest primitive at each grid point while capturing the correct sign of the global function. The methodology has optimal computational complexity and runs efficiently in distributed-memory architectures. We include 2D and 3D test cases along with a resolution study up to 0.512 trillion zones and 1,000 computer cores. The solution strategy can also be applied to other types of meshes or collections of primitives.

97 MATHEMATICS AND COMPUTING↗

Delayed leaf greening involves a major shift in the expression of cytosolic and mitochondrial ribosomes to plastid ribosomes in the highly phosphorus-use-efficient Hakea prostrata (Proteaceae)

Abstract Background and aims Hakea prostrata (Proteaceae) is a highly phosphorus-use-efficient plant native to southwest Australia. It maintains a high photosynthetic rate at low leaf phosphorus (P) and exhibits delayed leaf greening, a convergent adaptation that increases nutrient-use efficiency. This study aimed to provide broad physiological and gene expression profiles across leaf development, uncovering pathways leading from young leaves as nutrient sinks to mature leaves as low-nutrient, energy-transducing sources. Methods To explore gene expression underlying delayed greening, we analysed a de novo transcriptome for H. prostrata across five stages of leaf development. Photosynthesis and respiration rates, and foliar pigment, P and nitrogen (N) concentrations were determined, including the division of P into five biochemical fractions. Key results Transcripts encoding functions associated with leaf structure generally decreased in abundance across leaf development, concomitant with decreases in foliar concentrations of 85% for anthocyanins, 90% for P and 70% for N. The expression of genes associated with photosynthetic function increased during or after leaf expansion, in parallel with increases in photosynthetic pigments and activity, much later in leaf development than in species that do not have delayed greening. As leaves developed, transcript abundance for cytosolic and mitochondrial ribosomal proteins generally declined, whilst transcripts for chloroplast ribosomal proteins increased. Conclusions There was a much longer temporal separation of leaf cell growth from chloroplast development in H. prostrata than is found in species that lack delayed greening. Transcriptome-guided analysis of leaf development in H. prostrata provided insight into delayed greening as a nutrient-saving strategy in severely phosphorus-impoverished landscapes.

59 BASIC BIOLOGICAL SCIENCES↗

Twist-assisted all-antiferromagnetic tunnel junction in the atomic limit

Abstract Antiferromagnetic spintronics 1,2 shows great potential for high-density and ultrafast information devices. Magnetic tunnel junctions (MTJs), a key spintronic memory component that are typically formed from ferromagnetic materials, have seen rapid developments very recently using antiferromagnetic materials 3,4 . Here we demonstrate a twisting strategy for constructing all-antiferromagnetic tunnel junctions down to the atomic limit. By twisting two bilayers of CrSBr, a 2D antiferromagnet (AFM), a more than 700% nonvolatile tunnelling magnetoresistance (TMR) ratio is shown at zero field (ZF) with the entire twisted stack acting as the tunnel barrier. This is determined by twisting two CrSBr monolayers for which the TMR is shown to be derived from accumulative coherent tunnelling across the individual CrSBr monolayers. The dependence of the TMR on the twist angle is calculated from the electron-parallel momentum-dependent decay across the twisted monolayers. This is in excellent agreement with our experiments that consider twist angles that vary from 0° to 90°. Moreover, we also find that the temperature dependence of the TMR is, surprisingly, much weaker for the twisted as compared with the untwisted junctions, making the twisted junctions even more attractive for applications. Our work shows that it is possible to push nonvolatile magnetic information storage to the atomically thin limit.

Science & Technology - Other Topics↗

Grid Cyber-Security Strategy in an Attacker-Defender Model

The progression of cyber-attacks on the cyber-physical system is analyzed by the Probabilistic, Learning Attacker, and Dynamic Defender (PLADD) model. Although our research does apply to all cyber-physical systems, we focus on power grid infrastructure. The PLADD model evaluates the effectiveness of moving target defense (MTD) techniques. We consider the power grid attack scenarios in the AND configurations and OR configurations. In addition, we consider, for the first time ever, power grid attack scenarios involving both AND configurations and OR configurations simultaneously. Cyber-security managers can use the strategy introduced in this manuscript to optimize their defense strategies. Specifically, our research provides insight into when to reset access controls (such as passwords, internet protocol addresses, and session keys), to minimize the probability of a successful attack. Our mathematical proof for the OR configuration of multiple PLADD games shows that it is best if all access controls are reset simultaneously. For the AND configuration, our mathematical proof shows that it is best (in terms of minimizing the attacker's average probability of success) that the resets are equally spaced apart. We introduce a novel concept called hierarchical parallel PLADD system to cover additional attack scenarios that require combinations of AND and OR configurations.

97 MATHEMATICS AND COMPUTING↗

Mind the Gap: Building Simulation in the Architectural Design Studio

Building modelling and simulation approaches are increasingly being utilized in architectural design studios to guide and inform the design process and offer evidence-based feedback on proposed building performance. The development of intuitive and simplified simulation interfaces has greatly contributed to achieving this integration. One aspect that is often overlooked is the workflow that governs and regulates integrated design, which can have significant impacts on final design outcomes. Currently there are numerous software packages available for building performance simulations. This makes it challenging to select an appropriate tool that provides accurate results yet allows a designer to make informed architectural decisions with a designer-friendly interface. Furthermore, workflows to incorporate simulations into the design process proved to highly impact student’s project design integration. Yet, it is not clear what type of workflows are successful to achieve this goal, under what conditions, and/or for which building and site typologies. This paper addresses these issues by first reviewing three different workflows for integrating building performance simulation processes and highlighting their strengths and weaknesses. Second, a comparative case study approach was employed to test three of the most common workflows in three different integrated design architectural studios at the senior and vertical studio levels as well as in two courses that run parallel and complementary to the design studios. The workflows, processes, and the resultant student projects were further analyzed based on criteria for better integrated design and architectural excellence as outlined in the American Institute of Architects Committee on the Environment (AIA COTE) Top 10 Student’s Competition. Third, to situate the pedagogical case studies’ results within a larger context, a survey of the AIA COTE Top 10 student competition award recipients over the last five years was conducted. The results are summarized in a pedagogical framework that outlines best strategies of the type of workflows, software, design process used, methods to achieve desired interaction between design process and analytical feedback, and metrics for educators to evaluate the success of this integration and their learning outcomes in the design studio. The goal is to help bridge the gap between the building design and simulation within the design studio’s creative process for more integrated design outcomes.

Elzeyadi, Ihab↗

Synthetic Nuances to Maximize n-Type Organic Electrochemical Transistor and Thermoelectric Performance in Fused Lactam Polymers

A series of fully fused n-type mixed conduction lactam polymers p(g 7 NC n N), systematically increasing the alkyl side chain content, are synthesized via an inexpensive, nontoxic, precious-metal-free aldol polycondensation. Employing these polymers as channel materials in organic electrochemical transistors (OECTs) affords state-of-the-art n-type performance with p(g 7 NC 10 N) recording an OECT electron mobility of 1.20 x 10 -2 cm 2 V -1 s -1 and a μC* figure of merit of 1.83 F cm -1 V -1 s -1. In parallel to high OECT performance, upon solution doping with (4-(1,3-dimethyl-2,3-dihydro-1H-benzoimidazol-2-yl)phenyl)dimethylamine (N-DMBI), the highest thermoelectric performance is observed for p(g 7 NC 4 N), with a maximum electrical conductivity of 7.67 S cm -1 and a power factor of 10.4 μWm -1 K -2 . These results are among the highest reported for n-type polymers. Importantly, while this series of fused polylactam organic mixed ionic-electronic conductors (OMIECs) highlights that synthetic molecular design strategies to bolster OECT performance can be translated to also achieve high organic thermoelectric (OTE) performance, a nuanced synthetic approach must be used to optimize performance. Herein, we outline the performance metrics and provide new insights into the molecular design guidelines for the next generation of high-performance n-type materials for mixed conduction applications, presenting for the first time the results of a single polymer series within both OECT and OTE applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Fast Parallel Simulation Code for Interaction between Proto-Planetary Disk and Embedded Proto-Planets: Implementation for 3D Code

We develop a 3D simulation code for interaction between the proto-planetary disk and embedded proto-planets. The protoplanetary disk is treated as a three-dimensional (3D), self-gravitating gas whose motion is described by the locally isothermal Navier-Stokes equations in a spherical coordinate centered on the star. The differential equations for the disk are similar to those given in Kley et al. (2009) with a different gravitational potential that is defined in Nelson et al. (2000). The equations are solved by directional split Godunov method for the inviscid Euler equations plus operator-split method for the viscous source terms. We use a sub-cycling technique for the azimuthal sweep to alleviate the time step restriction. We also extend the FARGO scheme of Masset (2000) and modified in Li et al. (2001) to our 3D code to accelerate the transport in the azimuthal direction. Furthermore, we have implemented a reduced 2D (r, {theta}) and a fully 3D self-gravity solver on our uniform disk grid, which extends our 2D method (Li, Buoni, & Li 2008) to 3D. This solver uses a mode cut-off strategy and combines FFT in the azimuthal direction and direct summation in the radial and meridional direction. An initial axis-symmetric equilibrium disk is generated via iteration between the disk density profile and the 2D disk-self-gravity. We do not need any softening in the disk self-gravity calculation as we have used a shifted grid method (Li et al. 2008) to calculate the potential. The motion of the planet is limited on the mid-plane and the equations are the same as given in D'Angelo et al. (2005), which we adapted to the polar coordinates with a fourth-order Runge-Kutta solver. The disk gravitational force on the planet is assumed to evolve linearly with time between two hydrodynamics time steps. The Planetary potential acting on the disk is calculated accurately with a small softening given by a cubic-spline form (Kley et al. 2009). Since the torque is extremely sensitive to the position of the planet, we adopt the corotating frame that allows the planet moving only in radial direction if only one planet is present. This code has been extensively tested on a number of problems. For the earthmass planet with constant aspect ratio h = 0.05, the torque calculated using our code matches quite well with the the 3D linear theory results by Tanaka et al. (2002). The code is fully parallelized via message-passing interface (MPI) and has very high parallel efficiency. Several numerical examples for both fixed planet and moving planet are provided to demonstrate the efficacy of the numerical method and code.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Extending OpenSHMEM with Aggregation Support for Improved Message Rate Performance

OpenSHMEM is a highly efficient one-sided communication API that implements the PGAS parallel programming model, and is known for its low latency communication operations that can be mapped efficiently to RDMA capabilities of network interconnects. However, applications that use OpenSHMEM can be sensitive to point-to-point message rates, as many-to-many communication patterns can generate large amounts of small messages which tend to overwhelm network hardware that has predominantly been optimised for bandwidth over message rate. Additionally, many important emerging classes of problems such as data analytics are similarly troublesome for the irregular access patterns they employ. Message aggregation strategies have been proven to significantly enhance network performance, but their implementation often involves complex restructuring of user code, making them unwieldy. This paper shows how to combine the best qualities of message aggregation within the communication model of OpenSHMEM such that applications with small and irregular access patterns can improve network performance while maintaining their algorithmic simplicity. We do this by providing a path to a message aggregation framework called conveyors through a minimally intrusive OpenSHMEM extension introducing aggregation contexts that fit more naturally to the OpenSHMEM atomics, gets, and puts model. We test these extensions using four of the bale 3.0 applications which contain essential many-to-many access patterns to show how they can produce performance improvements of up to 65×.

Welch, Aaron↗

Performance upgrades to the MCNP6 burnup capability for large scale depletion calculations

The first MCNP based inline Monte Carlo depletion capability was officially released from the Radiation Safety Information and Computational Center as MCNPX 2.6.0. With the merger of MCNPX and MCNP5, MCNP6 combined the capability of both simulation tools, as well as providing new advanced technology, in a single radiation transport code. The new MCNP6 depletion capability was first showcased at the International Congress for Advancements in Nuclear Power Plants (ICAPP) meeting in 2012. At that conference the new capabilities addressed included the combined distributive and shared memory parallel architecture for the burnup capability, improved memory management, physics enhancements, and new predictability as compared to the H.B Robinson Benchmark. At Los Alamos National Laboratory, a special purpose cluster named “tebow,” was constructed such to maximize available RAM per CPU, as well as leveraging swap space with solid state hard drives, to allow larger scale depletion calculations (allowing for significantly more burnable regions than previously examined). As the MCNP6 burnup capability was scaled to larger numbers of burnable regions, a noticeable slowdown was realized.This paper details two specific computational performance strategies for improving calculation speedup: (1) retrieving cross sections during transport; and (2) tallying mechanisms specific to burnup in MCNP. To combat this slowdown new performance upgrades were developed and integrated into MCNP6 1.2.

42 ENGINEERING↗

Neglecting Model Parametric Uncertainty Can Drastically Underestimate Flood Risks

Abstract Floods drive dynamic and deeply uncertain risks for people and infrastructures. Uncertainty characterization is a crucial step in improving the predictive understanding of multi‐sector dynamics and the design of risk‐management strategies. Current approaches to estimate flood hazards often sample only a relatively small subset of the known unknowns, for example, the uncertainties surrounding the model parameters. This approach neglects the impacts of key uncertainties on hazards and system dynamics. Here we mainstream a recently developed method for Bayesian inference to calibrate a computationally expensive distributed hydrologic model. We compare three different calibration approaches: (a) stepwise line search, (b) precalibration or screening, and (c) the Fast Model Calibrations (FaMoS) approach. FaMoS deploys a particle‐based approach that takes advantage of the massive parallelization afforded by modern high‐performance computing systems. We quantify how neglecting parametric uncertainty and data discrepancy can drastically underestimate extreme flood events and risks. Precalibration improves prediction skill score over a stepwise line search. The Bayesian calibration improves the uncertainty characterization of model parameters and flood risk projections.

54 ENVIRONMENTAL SCIENCES↗

Revisiting Topographic Horizons in the Era of Big Data and Parallel Computing

Widely used to calculate illumination geometry forestimates of solar and emitted longwave radiation, and forcorrecting remotely sensed data for topographic effects, digitalelevation models (DEMs) are now extensive globally at 10–30-mspatial resolution and locally at spatial resolutions down to afew centimeters. Globally, regionally, or locally, elevation datasetshave many grid points. Many software packages calculate gradi-ents over every grid cell or point, but in the mountains, shadingby nearby terrain must also be assessed. Terrain may obscure aslope that would otherwise face the Sun. Four decadesago, a fastmethod to calculate topographic horizons at every point in anelevation grid required computations related only linearly to thesize of the grid, but grids now have so many points that parallelcomputing still provides an advantage. Exploiting parallelismover terrain grids can use alternative strategies: among columnsof a rotated grid, or simultaneously at multiple rotation angles,or on different tiles of a grid. On a multi-processor machine, theimprovement in computing time approaches 2/3 the number ofprocessors deployed,

Jeff Dozier↗

Demystifying asynchronous I/O Interference in HPC applications

With increasing complexity of HPC workflows, data management services need to perform expensive I/O operations asynchronously in the background, aiming to overlap the I/O with the application runtime. However, this may cause interference due to competition for resources: CPU, memory/network bandwidth. The advent of multi-core architectures has exacerbated this problem, as many I/O operations are issued concurrently, thereby competing not only with the application but also among themselves. Furthermore, the interference patterns can dynamically change as a response to variations in application behavior and I/O subsystems (e.g. multiple users sharing a parallel file system). Without a thorough understanding, I/O operations may perform suboptimally, potentially even worse than in the blocking case. To fill this gap, here we investigate the causes and consequences of interference due to asynchronous I/O on HPC systems. Specifically, we focus on multi-core CPUs and memory bandwidth, isolating the interference due to each resource. Then, we perform an in-depth study to explain the interplay and contention in a variety of resource sharing scenarios such as varying priority and number of background I/O threads and different I/O strategies: sendfile, read/write, mmap/write underlining trade-offs. The insights from this study are important both to enable guided optimizations of existing background I/O, as well as to open new opportunities to design advanced asynchronous I/O strategies.

97 MATHEMATICS AND COMPUTING↗

Neural Networks for Nuclear Reactions in MAESTROeX

We demonstrate the use of neural networks to accelerate the reaction steps in the MAESTROeX stellar hydrodynamics code. A traditional MAESTROeX simulation uses a stiff ODE integrator for the reactions; here, we employ a ResNet architecture and describe details relating to the architecture, training, and validation of our networks. Our customized approach includes options for the form of the loss functions, a demonstration that the use of parallel neural networks leads to increased accuracy, and a description of a perturbational approach in the training step that robustifies the model. We test our approach on millimeter-scale flames using a single-step, 3-isotope network describing the first stages of carbon fusion occurring in Type Ia supernovae. We train the neural networks using simulation data from a standard MAESTROeX simulation, and show that the resulting model can be effectively applied to different flame configurations. This work lays the groundwork for more complex networks, and iterative time-integration strategies that can leverage the efficiency of the neural networks.

79 ASTRONOMY AND ASTROPHYSICS↗

Career Development through Discreet Simulation Modeling - 19009

The cleanup mission of the Department of Energy (DOE) Office of River Protection (ORP) at the Hanford Site is to immobilize 212,000 m3 (56 million gallons) of nuclear waste, currently stored in 177 underground tanks, at the Waste Treatment and Immobilization Plant (WTP) and associated supplemental treatment facilities. Waste transfers to and from the WTP utilize the intricate Hanford Waste Feed Delivery (WFD) system, a network composed of: double-shell tanks (DSTs), single shell tanks (SSTs), the effluent treatment facility (ETF), a large number of pumps, valves, and leak detectors, and thousands of feet of underground transfer lines. This equipment must be prepared to safely and reliably transfer waste to the WTP. This study used discrete event simulation modeling to explore the behavior of the waste transfer system to identify areas of improvement. Two models and their subsequent results are presented: the Effluent Treatment Facility and the next-to-be-retrieved single-shell tanks (A/AX Farms). This study focused on identifying throughput bottlenecks, operational constraints, choke-point analysis, and redundancy in equipment, process design, and operational strategy. Analysis performed on the ETF OR Model indicated that the existing facility operation strategy will not result in a sufficient throughput to prevent WTP melter idling. Furthermore, in the absence of equipment failures, the facility was still shown to be unable to reach the required throughput to prevent melter idling. It is concluded that ETF system maintenance events are currently the main bottlenecks in the ETF facility. Analysis performed on the SST (A/AX) OR model indicate that the main drivers of the retrieval schedule are the hard retrieval start dates for each SST tank. In addition, it was concluded that operating the ventilation trains in parallel results in 15 times less retrieval interruptions and shortens the retrieval durations. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Towards the Structure Determination of a Modulated Protein Crystal: The Semicrystalline State of Profilin:Actin

One of the remaining challenges to structural biology is the solution of modulated structures. While small molecule crystallographers have championed this type of structure, to date, no modulated macromolecular structures have been determined. Modulation of the molecular structures within the crystal can produce satellite reflections or a superlattice of reflections in reciprocal space. We have developed the data collection methods and strategies that are needed to collect and analyze these data. If the macromolecule's crystal lattice is composed of physiologically relevant packing contacts, structural changes induced under physiological conditions can cause distortion relevant to the function and biophysical processes of the molecule making up the crystal. By careful measurement of the distortion, and the corresponding three-dimensional structure of the distorted molecule, we will visualize the motion and mechanism of the biological macromolecule(s). We have measured the modulated diffraction pattern produced by the semicrystalline state of profilin:actin crystals using highly parallel and highly monochromatic synchrotron radiation coupled with fine phi slicing (0.001-0.010 degrees) for structure determination. These crystals present these crystals present a unique opportunity to address an important question in structural biology. The modulation is believed to be due to the formation of actin helical filaments from the actin beta ribbon upon the pH-induced dissociation of profilin. To date, the filamentous state of actin has resisted crystallization and no detailed structures are available. The semicrystalline state profilin:actin crystals provides a unique opportunity to understand the many conformational states of actin. This knowledge is essential for understanding the dynamics underlying shape changes and motility of eukaryotic cells. Many essential processes, such as cytokinesis, phagocytosis, and cellular migration depend upon the capacity of the actin microfilament system to be restructured in a controlled manner via polymerization, depolymerization, severing, cross-linking, and anchorage. The structure the semicrystalline state of profilin:actin will challenge and validate current models of muscle contraction and cell motility. The methodology and theory under development will be easily extendable to other systems.

Borgstahl, G.↗

On Improving Efficiency of Differential Evolution for Aerodynamic Shape Optimization Applications

Differential Evolution (DE) is a simple and robust evolutionary strategy that has been proven effective in determining the global optimum for several difficult optimization problems. Although DE offers several advantages over traditional optimization approaches, its use in applications such as aerodynamic shape optimization where the objective function evaluations are computationally expensive is limited by the large number of function evaluations often required. In this paper various approaches for improving the efficiency of DE are reviewed and discussed. These approaches are implemented in a DE-based aerodynamic shape optimization method that uses a Navier-Stokes solver for the objective function evaluations. Parallelization techniques on distributed computers are used to reduce turnaround times. Results are presented for the inverse design of a turbine airfoil. The efficiency improvements achieved by the different approaches are evaluated and compared.

Madavan, Nateri K.↗

Direct Functionalization of Established 3D-Printed Aza-Michael Liquid Crystal Elastomers with Donor–Acceptor Stenhouse Adducts

Extrusion 3D printing has advanced the manufacturing of complex liquid crystal elastomer (LCE) architectures. In parallel, donor–acceptor Stenhouse adducts (DASAs), a class of white-light-responsive photoswitches, have enabled both photochemical and photothermal LCE actuation. However, DASA–LCEs have yet to be extruded and 3D-printed. Two key challenges exist: DASA’s inherent sensitivity to heat and radicals can lead to degradation during ink preparation and printing, and small changes in the concentration of the added DASA component impact the properties of the extrudable ink, requiring reoptimization of well-established 3D-printing protocols. To overcome these challenges, we present a post-printing functionalization strategy that circumvents these limitations. Residual secondary amines, inherent to inks synthesized via standard aza-Michael addition, serve as active sites for covalent attachment of DASA photoresponsive groups following printing and cross-linking. Our method means that DASAs can be directly grafted onto 3D-printed aza-Michael LCEs without modifying the ink formulation or processing. The resulting DASA–LCEs exhibit wavelength tunability within the visible range and a variety of photothermal and photochemical responses. The post-functionalization can occur within 2 min and enables spatial control of the DASA concentration, producing films with tunable color gradients and locally varied photothermal and photochemical responses under visible light. In conclusion, this approach enables the rapid fabrication of DASA-based light-responsive LCEs using established ink formulations with the potential for the design of complex 3D architectures.

3D printing↗

Switching Modes for Reduction of Peak Voltage Transients in GaN-Based Three Level ANPC Inverter

The switching pole voltage can transition between P,0, and N states in a three-level active neutral point clamped (ANPC) inverter. State transitions between 0-P and 0-N can be realized with different switching modes. Three switching modes- Short 1 , Short 2 , and Full are studied considering the effect of capacitive current paths. The role of clamping and inner switch on/off conditions is determined using simplified equivalent circuit models. Short 2 with clamping switch off is beneficial for turn-on overvoltage suppression but results in high turn-off overvoltage. Full mode with a parallel current sharing path is effective for turn-off overvoltage suppression but leads to high turn-on overvoltage. Hence, a Modified Full mode is proposed to achieve the simultaneous objective of overvoltage suppression at both turn-on and turn-off transient during high load currents. This benefits 3L ANPC operation at a low power factor with a high load current at fundamental voltage zero crossing points. Here, a 650V GaN-based three-level ANPC inverter prototype is designed and used to evaluate the switching modes and overvoltage suppression strategies through double pulse tests and continuous operation.

42 ENGINEERING↗