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 1,189 records · Page 66

Linking plant hydraulics and beta diversity in tropical forests

In tropical forests, studies of xylem traits governing water transport through plants, or ‘hydraulic architecture’ (Tyree et al., 1991), and changes in species composition across environmental gradients, or ‘beta diversity’ (Gentry, 1988; Ackerly & Cornwell, 2007), have progressedmostly in parallel until recently (Hao et al., 2008; Bartlett et al., 2016). In this issue of New Phytologist, Cosme et al. (pp. 000–5 000) present a timely contribution to the intersection of plant hydraulic architecture (HA) with trait-based community ecology. Building on previous biogeographical work that demonstrated shifts in species composition (beta diversity) across a gradient from valleys to plateaus in central Amazonia (Schietti et al., 2014), Cosme et al. explore how variation in HA might underpin this sorting, sampling pairs of congeneric species restrictedmostly to either plateau or valley habitats. Valley species had significantly lower wood density and higher hydraulically-weighted vessel diameter and vessel area. By contrast, trees with some of the largest hydraulically-weighted vessel diameters existed in tall, deciduous plateau species, while the leaf: sapwood area ratio decreased with height in valley but not plateau species. These intriguing results suggest that species differentiation in water transport traits mediate edaphic filtering along the valley-toplateau gradient, in contrast to previous work where wood mechanical support mediated valley-to-plateau environmental filtering (Fortunel et al., 2014).

Christoffersen, Bradley↗

Attribution of North American Subseasonal Precipitation Prediction Skill

The skill of NOAA’s official monthly U.S. precipitation forecasts (issued in the middle of the prior month) has historically been low, having shown modest skill over the southern United States, but little or no skill over large portions of the central United States. The goal of this study is to explain the seasonal and regional variations of the North American subseasonal (weeks 3–6) precipitation skill, specifically the reasons for its successes and its limitations. The performances of multiple recent-generation model reforecasts over 1999–2015 in predicting precipitation are compared to uninitialized simulation skill using the atmospheric component of the forecast systems. This parallel analysis permits attribution of precipitation skill to two distinct sources: one due to slowly evolving ocean surface boundary states and the other to faster time-scale initial atmospheric weather states. A strong regionality and seasonality in precipitation forecast performance is shown to be analogous to skill patterns dictated by boundary forcing constraints alone. The correspondence is found to be especially high for the North American pattern of the maximum monthly skill that is achieved in the reforecast. The boundary forcing of most importance originates from tropical Pacific SST influences, especially those related to El Niño–Southern Oscillation. Furthermore, we discuss physical constraints that may limit monthly precipitation skill and interpret the performance of existing models in the context of plausible upper limits.

54 ENVIRONMENTAL SCIENCES↗

Basic Physics Data: Measurement of Neutron Multiplicity from Induced Fission

From October 1 to October 17 a team of researchers from UM visited the LANSCE facility for an experiment during beam-time allotted from October 4 to October 17. A total of 24 detectors were used at LANSCE including liquid organic scintillation detectors (EJ-309), NaI scintillation detectors, and Li-6 enriched glass detectors. It is a double time-offlight (TOF) measurement using spallation neutrons generated by a target bombarded with pulsed high-energy protons. The neutrons travel to an LLNL-manufactured parallel plate avalanche chamber (PPAC) loaded with thin U-235 foils in which fission events are induced. The generated fission neutrons and photons are then detected in a detector array designed and built at UM and shipped to LANSCE. Preparations were made at UM, where setup and proposed detectors were tested. The UM equipment was then shipped to LANSCE for use at the 15L beam of the weapons neutron research (WNR) facility.

43 PARTICLE ACCELERATORS↗

Unraveling the Determinant Mechanisms in Flow-Mediated Crystal Growth and Phase Behaviors

To uncover the critical mechanisms responsible for mesoscopic level development during flow-mediated crystal growth, we develop a semi-two-way hydrodynamic coupled structural phase-field crystal formalism (HXPFC-s2). The new formalism, inspired by previous attempts at coupling hydrodynamic and phase-field crystal (PFC) equations, allows for studying mesoscopic flow-mediated crystallization at diffusive timescales pertinent to industrial applications. Unlike previous efforts, the devised coupling to the structural PFC (XPFC) equations allows generalization to more complex crystal structures through explicit parameterization of the direct correlation function (DCF). Utilizing the HXPFC-s2 formalism, we seek to uncover the determinant physical mechanisms in crystallization under simple shear flows by comparing temperature-driven crystallization to flow-mediated crystallization under varying flow-strengths. Parallels and deviations of under-cooling and flow-strength effects on crystal growth are drawn using the crystal cluster-size and system ordering time evolutions. In doing so, we identify scaling behaviors with a Peclet-like number, Pe∼, a critical Peclet-like number, Pe∼*, and flow-field-crystal plane-dependent interactions. Our findings may be relevant for controlling crystal growth and phase behaviors in flow applications.

Willis, L. Connor (ORCID:0009000961321848)↗

Parallel Tensor Compression for Large-Scale Scientific Data.

As parallel computing trends towards the exascale, scientific data produced by high-fidelity simulations are growing increasingly massive. For instance, a simulation on a three-dimensional spatial grid with 512 points per dimension that tracks 64 variables per grid point for 128 time steps yields 8 TB of data. By viewing the data as a dense five way tensor, we can compute a Tucker decomposition to find inherent low-dimensional multilinear structure, achieving compression ratios of up to 10000 on real-world data sets with negligible loss in accuracy. So that we can operate on such massive data, we present the first-ever distributed memory parallel implementation for the Tucker decomposition, whose key computations correspond to parallel linear algebra operations, albeit with nonstandard data layouts. Our approach specifies a data distribution for tensors that avoids any tensor data redistribution, either locally or in parallel. We provide accompanying analysis of the computation and communication costs of the algorithms. To demonstrate the compression and accuracy of the method, we apply our approach to real-world data sets from combustion science simulations. We also provide detailed performance results, including parallel performance in both weak and strong scaling experiments.

59 BASIC BIOLOGICAL SCIENCES↗

Observation of Persistent Zero Modes and Superconducting Vortex Doublets in UTe 2

Superconducting vortices can reveal electron pairing details and nucleate topologically protected states. Yet, vortices of bulk spin-triplet superconductors have never been visualized at the atomic scale. Recently, UTe 2 has emerged as a prime spin-triplet superconductor, but its superconducting order parameter is elusive, and whether time-reversal symmetry (TRS) is broken remains unsettled. Here, in this work, we visualize vortices on the (011) surface of ultraclean UTe 2 single crystals (T c = 2.1 K) using scanning tunneling microscopy (STM). We introduce $\frac{d^{2}𝐼}{d𝑉^{2}}$ imaging as an effective technique for vortex visualization in superconductors with substantial residual zero-energy density of states (DOS), as in UTe 2 . Anisotropic single-flux-quantum vortices, with coherence lengths of ∼12 nm (4 nm) parallel (perpendicular) to the a -axis, form a triangular vortex lattice (VL) under a small out-of-plane magnetic field. The invariance of vortex structures and VL under changes of field polarity and cooling history strongly supports time-reversal invariant superconductivity under zero field. At vortex cores (VCs), nonsplit, spectrally sharp zero-bias conductance peaks (ZBPs) persist up 8 T that are consistent with symmetry-protected Majorana zero modes (MZMs) in a topological vortex line. Close examination of vortex structures reveals a mirror-asymmetric doublet─one with ZBPs and another with an enhanced apparent gap, possibly originating from a field-induced multicomponent order parameter.

UTe2↗

Nonlinear fluid simulation of particle and heat fluxes during burst of ELMs on DIII-D with BOUT++ code [Fluid Simulation of Particle and Heat Fluxes during Burst of ELMs on DIID with BOUT++ code]

In order to study the distribution and evolution of the transient particle and heat fluxes during edge-localized mode (ELM) bursts, a BOUT++ six-field two-fluid model based on the Braginskii equations with non-ideal physics effects is used to simulate pedestal collapse in divertor geometry. We used the profiles from the DIII-D H-mode discharge #144382 with fast target heat flux measurements as the initial conditions for the simulations. Moreover, a flux-limited parallel thermal conduction is used with three values of the flux-limiting coefficient ${{\alpha}_{j}}$ , free streaming model with ${{\alpha}_{j}}=1$ , sheath-limit with ${{\alpha}_{j}}=0.05$ , and one value in between. The studies show that a 20 times increase in ${{\alpha}_{j}}$ leads to ~6 times increase in the heat flux amplitude to both the inner and outer targets, and the widths of the fluxes are also expanded. The sheath-limit model of flux-limiting coefficient is found to be the most appropriate one, which shows ELM sizes close to the measurements. The evolution of the density profile during the burst of ELMs of DIII-D discharge #144382 is simulated, and the collapse in width and depth of ${{n}_{\text{e}}}$ are reproduced at different time steps. The growing process of the profiles for the heat flux at divertor targets during the burst of ELMs measured by IRTV (infrared television) is also reproduced by this model. The widths of heat fluxes towards targets are a little narrower, and the peak amplitudes are twice the measurements possibly due to the lack of a model of divertor radiation which can effectively reduce the heat fluxes. The magnetic flutter combined with parallel thermal conduction is found to be able to increase the total heat loss by around 33% since the magnetic flutter terms provide the additional conductive heat transport in the radial direction. Finally, the heat flux profile at both the inner and outer targets is obviously broadened by magnetic flutter. The lobe structures near the X-point at LFS are both broadened and elongated due to the magnetic flutter.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Accelerating Simulation for High-Fidelity PV Inverter System Reliability Assessment with High-Performance Computing

The overall cost of photovoltaic (PV) systems has shown a downward trend during the last decade; however, PV inverter failures account for the highest cost of operation and maintenance. To address this, reliability tools with powerful computation and better accuracy are required for the lifetime prediction and degradation evaluation of PV inverters. This paper proposes an event-driven parallel computing-based simulator. The proposed simulator applies high-performance computing techniques and other accessory optimization techniques-including cluster merging, adaptive model updates, and steady-state identification-to make reliability assessments for PV inverters under given input mission profiles and operating conditions with high efficiency and high fidelity. The main idea of the simulator and its workflow are introduced. Then, a demo PV inverter system simulator is implemented, and the speedup of the total simulations of the switching model reaches 123.03 times.

high-performance computing↗

TuckerCompressMPI v. 1.0

As parallel computing trends towards the exascale, scientific data produced by high-fidelity simulations are growing increasingly massive. For instance, a simulation on a three-dimensional spatial grid with 512 points per dimension that tracks 64 variables per grid point for 128 time steps yields 8 TB of data. By viewing the data as a dense five-way tensor, we can compute a Tucker decomposition to find inherent low-dimensional multilinear structure, achieving compression ratios of up to 10000 on real-world data sets with negligible loss in accuracy. So that we can operate on such massive data, we present the first-ever distributed-memory parallel implementation for the Tucker decomposition, whose key computations correspond to parallel linear algebra operations, albeit with nonstandard data layouts. Our approach specifies a data distribution for tensors that avoids any tensor data redistribution, either locally or in parallel. This software provides a method for compressing large-scale multiway data.

Austin, Woody↗

OpenACC acceleration of an unstructured CFD solver based on a reconstructed discontinuous Galerkin method for compressible flows

Here, an OpenACC directive-based graphics processing unit (GPU) parallel scheme is presented for solving the compressible Navier–Stokes equations on 3D hybrid unstructured grids with a third-order reconstructed discontinuous Galerkin method. The developed scheme requires the minimum code intrusion and algorithm alteration for upgrading a legacy solver with the GPU computing capability at very little extra effort in programming, which leads to a unified and portable code development strategy. A face coloring algorithm is adopted to eliminate the memory contention because of the threading of internal and boundary face integrals. A number of flow problems are presented to verify the implementation of the developed scheme. Timing measurements were obtained by running the resulting GPU code on one Nvidia Tesla K20c GPU card (Nvidia Corporation, Santa Clara, CA, USA) and compared with those obtained by running the equivalent Message Passing Interface (MPI) parallel CPU code on a compute node (consisting of two AMD Opteron 6128 eight-core CPUs (Advanced Micro Devices, Inc., Sunnyvale, CA, USA)). Speedup factors of up to 24× and 1.6× for the GPU code were achieved with respect to one and 16 CPU cores, respectively. The numerical results indicate that this OpenACC-based parallel scheme is an effective and extensible approach to port unstructured high-order CFD solvers to GPU computing.

97 MATHEMATICS AND COMPUTING↗

Extending the Roofline Model for Asynchronous Many-Task Runtimes

New architectural trends and extreme scale parallelism challenge the efficient mapping of applications to large scale and emerging systems. Fine-grained Asynchronous Manytask Runtimes (AMR) offer a unique flexibility to hide the varying latencies of contentious operations, and facilitate the effective exploitation of abundant resources. By decomposing work into smaller asynchronous chunks, the impact of unpredictable operations can be mitigated through a tighter overlapping of tasks. Unfortunately, the cost of this decomposition comes at a potentially steep price in the form of runtime overhead which can at times rival the cost of the computation itself. This leads to the question, how large should a task be? A common practice for application developers is to experimentally determine the granularity of a task after a code has been parallelized. Instead, we propose a new methodology based on an extended Roofline model to provide practical upper bounds on the throughput performance of an application. First, we extend the Roofline model to support not only latency hiding analysis, but also a multidimensional amortized analysis. By combining this new methodology with a serial application and an AMR implementation, we can predict the worst case runtime overhead attribution of individual runtime features prior to the development of parallel code. Thus, this runtimecentric methodology can provide a vehicle for application/runtime codesign by providing a comprehensive bottleneck analysis based on existing runtime features.

Suetterlein, Joshua D.↗

Steady-state benchmarks of DK4D: A time-dependent, axisymmetric drift-kinetic equation solver

The DK4D code has been written to solve a set of time-dependent, axisymmetric, finite-Larmor-radius drift-kinetic equations (DKEs) for the non-Maxwellian part of the electron and ion distribution functions using the full, linearized Fokker–Planck–Landau collision operator. The plasma is assumed to be in the low- to finite-collisionality regime, as is found in the cores of modern and future magnetic confinement fusion experiments. Each DKE is formulated such that the perturbed distribution function carries no net density, parallel momentum, or kinetic energy. Rather, these quantities are contained within the background Maxwellians and would be evolved by an appropriate set of extended magnetohydrodynamic (MHD) equations. This formulation allows for straight-forward coupling of DK4D to existing extended MHD time evolution codes. DK4D uses a mix of implicit and explicit temporal representations and finite element and spectral spatial representations. These, along with other computational methods used, are discussed extensively. Steady-state benchmarks are then presented comparing the results of DK4D to expected analytic results at low collisionality, qualitatively, and to the Sauter analytic fits for the neoclassical conductivity and bootstrap current, quantitatively. These benchmarks confirm that DK4D is capable of solving for the correct, gyroaveraged distribution function in stationary magnetic equilibria. Furthermore, the results presented demonstrate how the exact drift-kinetic solution varies with collisionality as a function of the magnetic moment and the poloidal angle.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Short-time polymerization and convolution kinetics in the analysis of the initiation step of isoprene polymerization over neodymium catalyst

The paper considers a convolution model for polymerization of 1,3-dienes in the presence of organometallic catalysts. It is shown that initial steps of polymerization can be represented as a set of shifted parallel reactions of the chain propagation, in which shift value are governed by the kinetics of the initiation step. In this case, polymer accumulation curve is described by an equation containing convolution of functions which allow one to solve an inverse problem of restoring the function of active site accumulation. Based on the example of isoprene polymerization over neodymium catalyst, it was demonstrated that the results of the inverse problem solution are in good agreement with the experimental data on concentration of the active sites. These data were obtained by a method of controlled introduction of inhibitor in short-time polymerization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Picosecond UV single photon detectors with lateral drift field: Concept and technologies

Group III–V semiconductor materials are being considered as a Si replacement for advanced logic devices for quite some time. Advances in III–V processing technologies, such as interface and surface passivation, large area deep submicron lithography with high-aspect ratio etching primarily driven by the metal-oxide-semiconductor field-effect transistor development can also be used for other applications. In this paper we will focus on photodetectors with the drift field parallel to the surface. We compare the proposed concept to the state-of-the-art Si-based technology and discuss requirements which need to be satisfied for such detectors to be used in a single photon counting mode in blue and ultraviolet spectral region with about 10 ps photon timing resolution essential for numerous applications ranging from high-energy physics to medical imaging.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Experimental and modeling studies of IPDI-based polyurea elastomers – The role of hard segment fraction

Segmented polyureas (PUa) are industrially important class of polymers widely used in coatings, sealant, and adhesive applications. Here, in this paper, we report synthesis, characterization, and modeling of Isophorone Diisocyanate-Diethyl-Toluene-Diamine-Polyether amine (IPDI-DETDA-PO PUa) with varied hard segment contents of 20, 30, and 40 weight percent. For each of the three materials, we study its structure and phase behavior using FTIR, DSC, and TEM, and clearly show the presence of microphase separation between the hard and soft nanodomains. We then measure the linear viscoelastic response of the PUa-s using DMA (frequency sweeps at multiple temperatures). The DMA data are shown to obey the time-temperature superposition. Finally, we develop a new micromechanical model describing the DMA results; the model describes a phase-separated PUa as two “Fractional-order Maxwell gels” branches, connected in parallel, with the first FMG branch representing the “percolated hard phase" and the second one modeling the “filled soft phase". In agreement with the earlier thermodynamic theories, the volume-fraction of the percolated hard phase is related to the hard segment weight-fraction (HSWF), defined as the combined mass of IPDI and DETDA normalized to the total mass of the polymer. The data and model are found to be in a good qualitative and quantitative agreement.

36 MATERIALS SCIENCE↗

Ultrarelativistic electron butterfly distributions created by parallel acceleration due to magnetosonic waves

The Van Allen Probe observations during the recovery phase of a large storm that occurred on 17 March 2015 showed that the ultrarelativistic electrons at the inner boundary of the outer radiation belt (L* = 2.6–3.7) exhibited butterfly pitch angle distributions, while the inner belt and the slot region also showed evidence of sub-MeV electron butterfly distributions. In this paper, strong magnetosonic waves were observed in the same regions and at the same time periods as these butterfly distributions. Moreover, when these magnetosonic waves extended to higher altitudes (L* = 4.1), the butterfly distributions also extended to the same region. Combining test particle calculations and Fokker-Planck diffusion simulations, we successfully reproduced the formation of the ultrarelativistic electron butterfly distributions, which primarily result from parallel acceleration caused by Landau resonance with magnetosonic waves. Finally, the coexistence of ultrarelativistic electron butterfly distributions with magnetosonic waves was also observed in the 24 June 2015 storm, providing further support that the magnetosonic waves play a key role in forming butterfly distributions.

79 ASTRONOMY AND ASTROPHYSICS↗

Magnetohydrodynamic pressure drop and flow balancing of liquid metal flow in a prototypic fusion blanket manifold

Understanding magnetohydrodynamic (MHD) phenomena associated with the flow of electrically conducting fluids in complex geometry ducts subject to a strong magnetic field is required to effectively design liquid metal (LM) blankets for fusion reactors. Specifically, accurately predicting the 3D MHD pressure drop and flow distribution is important. To explore these topics, we simulate a LM MHD flow through an electrically non-conducting prototypic manifold for a wide range of flow and geometry parameters using a 3D MHD solver, HyPerComp incompressible MHD solver for arbitrary geometry. The reference manifold geometry consists of a rectangular feeding duct which suddenly expands such that the duct thickness in the magnetic field direction abruptly increases by a factor r exp . Downstream of the sudden expansion, the LM is distributed into several parallel channels. As a first step in qualifying the flow, a magnitude of the curl of the induced Lorentz force was used to distinguish between inviscid, irrotational core flows and boundary and internal shear layers where inertia and/or viscous forces are important. Scaling laws have been obtained which characterize the 3D MHD pressure drop and flow balancing as a function of the flow parameters and the manifold geometry. Associated Hartmann and Reynolds numbers in the computations were ~10 3 and ~10 1 -10 3 , respectively, while r exp was varied from 4 to 12. An accurate model for the pressure drop was developed for the first time for inertial-electromagnetic and viscous-electromagnetic regimes based on 96 computed cases. Analysis reflects that flow balance can be improved by lengthening the distance between the manifold inlet and the entrances of the parallel channels by utilizing the effect of flow transitioning to a quasi-two-dimensional state in the expansion region of the manifold.

42 ENGINEERING↗

HPC for Optimizing Process Parameters to Control Material Evolution in Seamless Induction Hardening of Wind Turbine Main Shaft Bearings

Work proposed in this project focused on understanding the effect of martensitic transformation in the steel on the potential for cracking during seamless induction hardening (SIH) as a function of process conditions to allow the process to optimally scale up. Large-scale, three-dimensional phase-field simulations of martensitic transformation were performed using MEUMAPPS-SS (Microstructure Evolution Using Massively Parallel Phase-field Simulations – Solid State) code developed at Oak Ridge National Laboratory. The simulations were guided by location-specific thermal history generated by experimental measurements of time-temperature history generated at The Timken Company. The simulations were able to capture the morphological evolution of the martensite variants in an Fe-1.0C-1.5Cr steel based on the Nishiyama-Wasserman (NW) orientation relationship. The simulations were also able to quantify the stress-state at the interface between impinging martensite variants. The simulations indicated that the magnitude of the various stress and strain components were dependent on the sizes of the impinging plates with a reduction in these quantities with reduced plate size in agreement with experimental findings. The results obtained from the simulations will be used to guide the optimization of the alloy thermal conditions to eliminate quench cracking during SIH of bearing steels.

99 GENERAL AND MISCELLANEOUS↗