Search NASA⌕ Search

SEARCH · Search NASA

Results for “parallel cluster”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRS2006

Forecasting of Storm-Surge Floods Using ADCIRC and Optimized DEMs

Increasing the accuracy of storm-surge flood forecasts is essential for improving preparedness for hurricanes and other severe storms and, in particular, for optimizing evacuation scenarios. An interactive database, developed by WorldWinds, Inc., contains atlases of storm-surge flood levels for the Louisiana/Mississippi gulf coast region. These atlases were developed to improve forecasting of flooding along the coastline and estuaries and in adjacent inland areas. Storm-surge heights depend on a complex interaction of several factors, including: storm size, central minimum pressure, forward speed of motion, bottom topography near the point of landfall, astronomical tides, and, most importantly, maximum wind speed. The information in the atlases was generated in over 100 computational simulations, partly by use of a parallel-processing version of the ADvanced CIRCulation (ADCIRC) model. ADCIRC is a nonlinear computational model of hydrodynamics, developed by the U.S. Army Corps of Engineers and the US Navy, as a family of two- and three-dimensional finite-element-based codes. It affords a capability for simulating tidal circulation and storm-surge propagation over very large computational domains, while simultaneously providing high-resolution output in areas of complex shoreline and bathymetry. The ADCIRC finite-element grid for this project covered the Gulf of Mexico and contiguous basins, extending into the deep Atlantic Ocean with progressively higher resolution approaching the study area. The advantage of using ADCIRC over other storm-surge models, such as SLOSH, is that input conditions can include all or part of wind stress, tides, wave stress, and river discharge, which serve to make the model output more accurate. To keep the computational load manageable, this work was conducted using only the wind stress, calculated by using historical data from Hurricane Camille, as the input condition for the model. Hurricane storm-surge simulations were performed on an eight-node Linux computer cluster. Each node contained dual 2-GHz processors, 2GB of memory, and a 40GB hard drive. The digital elevation model (DEM) for this region was specified using a combination of Navy data (over water), NOAA data (for the coastline), and optimized Interferometric Synthetic Aperture Radar data (over land). This high-resolution topographical data of the Mississippi coastal region provided the ADCIRC model with improved input with which to calculate improved storm-surge forecasts.

Valenti, Elizabeth↗
DOE OSTI2025

Discovery of a Ferromagnetic Nickel Chalcogenide Nanocluster Ni 3 S 3 H(PEt 3 ) 5

Atomically precise ligated nanoclusters (NC) are promising cluster-based materials with novel molecular architectures and tunable magnetic properties. Herein, the synthesis and characterization of a nickel sulfide NC Ni 3 S 3 H(PEt 3 ) 5 (PEt 3 = triethylphosphine) with distinct magnetic properties are reported. Magnetization measurements reveal its magnetic moment of 1.5 µ B in the solid phase, consistent with the existence of one unpaired electron predicted by density functional theory (DFT) calculations. Additionally, experimental measurements indicate the presence of ferromagnetic ordering within each Ni 3 S 3 H(PEt 3 ) 5 NC and strong coercivity at temperatures below 20 K. Ion mobility-mass spectrometry is employed in conjunction with DFT calculations and collision cross-section simulations to investigate the structure of the isolated Ni 3 S 3 H(PEt 3 ) 5 . Theoretical studies show that [Ni 3 S 3 H(PEt 3 ) 5 ] + has a planar Ni 3 S 3 core where three Ni atoms are arranged in a triangle with three bridging S atoms residing in the same plane. This structure is preserved in both solution and solid phases, which is confirmed by spectroscopic studies of Ni 3 S 3 H(PEt 3 ) 5 . Additionally, DFT calculations indicate that all spins at the Ni sites are aligned parallel, confirming the presence of ferromagnetic coupling. Overall, this study provides key insights into the structure and magnetic properties of Ni 3 S 3 H(PEt 3 ) 5 , which will facilitate the design of new NC-based magnetic materials.

Nickel sulfide nanocluster↗
DOE OSTI2020

Bridging the Transuranics with Uranium(IV) Sulfate Aqueous Species and Solid Phases

Isolating isomorphic compounds of tetravalent actinides (i.e., Th IV , U IV , Np IV , and Pu IV ) improve our understanding of the bonding behavior across the series, in addition to their relationship with tetravalent transition metals (Zr and Hf) and lanthanides (Ce). Similarities between these tetravalent metals are particularly illuminated in their hydrolysis and condensation behavior in aqueous systems, leading to polynuclear clusters typified by the hexamer [M IV 6 O 4 (OH) 4 ] 12+ building block. Prior studies have shown the predominance and coexistence of smaller species for Th IV (monomers, dimers, and hexamers) and larger species for U IV , Np IV , and Pu IV (including 38-mers and 70-mers). We show here that aqueous uranium(IV) sulfate also displays behavior similar to that of Th IV (and Zr IV ) in its isolated solidphase and solution speciation. Two single-crystal X-ray structures are described: a dihydroxide-bridged dimer (U 2 ) formulated as U 2 (OH) 2 (SO 4 ) 3 (H 2 O) 4 and a monomer-linked hexamer framework (U-U 6 ) as (U(H 2 O) 3.5 ) 2 U 6 O 4 (OH) 4 (SO 4 ) 10 (H 2 O) 9 . These structures are similar to those previously described for ThIV. Moreover, cocrystallization of monomer and dimer and of dimer and monomer–hexamer phases for both Th IV (prior) and U IV (current) indicates the coexistence of these species in solution. Because it was not possible to effectively study the sulfate-rich solutions via X-ray scattering from which U 2 and U-U 6 crystallized, we provide a parallel solution speciation study in low sulfate conditions, as a function of the pH. Raman spectroscopy, UV–vis spectroscopy, and small-angle X-ray scattering of these show decreasing sulfate binding, increased hydrolysis, increased species size, and increased complexity, with increasing pH. Furthermore, this study describes a bridge across the first half the actinide series, highlighting U IV similarities to Th IV , in addition to the previously known similarities to the transuranic elements.

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

MemXCT: Design, Optimization, Scaling, and Reproducibility of X-Ray Tomography Imaging

Here, this work extends our previous research entitled "MemXCT: Memory-centric X-ray CT Reconstruction with Massive Parallelization" that was originally published at SC19 conference (Hidayetoglu et al., 2019) with reproducibility of the computational imaging performance. X-ray computed tomography (XCT) is regularly used at synchrotron light sources to study the internal morphology of materials at high resolution. However, experimental constraints, such as radiation sensitivity, can result in noisy or undersampled measurements. Further, depending on the resolution, sample size and data acquisition rates, the resulting noisy dataset can be in the order of terabytes. Advanced iterative reconstruction techniques can produce high-quality images from noisy measurements, but their computational requirements have made their use an exception rather than the rule. We propose a novel memory-centric approach that avoids redundant computations at the expense of additional memory complexity. We develop a memory-centric iterative reconstruction system, MemXCT, that uses an optimized SpMV implementation with two-level pseudo-Hilbert ordering and multi-stage input buffering. We evaluate MemXCT on various supercomputer architectures involving KNL and GPU. MemXCT can reconstruct a large (11Kx11K) mouse brain tomogram in 10 seconds using 4096 KNL nodes (256K cores). The results presented in our original article at the SC19 were based on large-scale supercomputing resources. The MemXCT application was selected for the Student Cluster Competition (SCC) Reproducibility Challenge and evaluated on a variety of cloud computing resources by universities around the world in the SC20 conference. We summarize the results of the top-ranked SCC Reproducibility Challenge teams and identify the most pertinent measures for ensuring the reproducibility of our experiments in this article.

47 OTHER INSTRUMENTATION↗
NASA NTRS1996

Non-Linear Cosmological Power Spectra in Real and Redshift Space

We present an expression for the non-linear evolution of the cosmological power spectrum based on Lagrangian trajectories. This is simplified using the Zel'dovich approximation to trace particle displacements, assuming Gaussian initial conditions. The model is found to exhibit the transfer of power from large to small scales expected in self-gravitating fields. Some exact solutions are found for power-law initial spectra. We have extended this analysis into red-shift space and found a solution for the non-linear, anisotropic redshift-space power spectrum in the limit of plane-parallel redshift distortions. The quadrupole-to-monopole ratio is calculated for the case of power-law initial spectra. We find that the shape of this ratio depends on the shape of the initial spectrum, but when scaled to linear theory depends only weakly on the redshift-space distortion parameter, beta. The point of zero-crossing of the quadrupole, kappa(sub o), is found to obey a simple scaling relation and we calculate this scale in the Zel'dovich approximation. This model is found to be in good agreement with a series of N-body simulations on scales down to the zero-crossing of the quadrupole, although the wavenumber at zero-crossing is underestimated. These results are applied to the quadrupole-to-monopole ratio found in the merged QDOT plus 1.2-Jy-IRAS redshift survey. Using a likelihood technique we have estimated that the distortion parameter is constrained to be beta greater than 0.5 at the 95 percent level. Our results are fairly insensitive to the local primordial spectral slope, but the likelihood analysis suggests n = -2 un the translinear regime. The zero-crossing scale of the quadrupole is k(sub 0) = 0.5 +/- 0.1 h Mpc(exp -1) and from this we infer that the amplitude of clustering is sigma(sub 8) = 0.7 +/- 0.05. We suggest that the success of this model is due to non-linear redshift-space effects arising from infall on to caustic and is not dominated by virialized cluster cores. The latter should start to dominate on scales below the zero-crossing of the quadrupole, where our model breaks down.

Taylor, A. N.↗
NASA NTRS2012

Simulations of Magnetic Reconnection - Kinetic Mechanisms Underlying the Fluid Description of Ions

Because of its ability to transfer the energy stored in magnetic field together with the breaking of the flux freezing constraint, magnetic reconnection is considered as one of the most important phenomena in plasma physics. When it happens in a collision less environment such as the terrestrial magnetosphere, it should a priori be modelled with in the framework of kinetic physics. The evidence of kinetic features has incidentally for a long time, been shown by researchers with the help of both numerical simulations and satellite observations. However, most of our understanding of the process comes from the more intuitive fluid interpretation with simple closure hypothesis which do not include kinetic effects. To what extent are these two separate descriptions of the same phenomenon related? What is the role of kinetic effects in the averaged/fluid dynamics of reconnection? This thesis addresses these questions for the proton population in the particular case of anti parallel merging with the help of 2D Hybrid simulations. We show that one can not assume, as is usually done, that the acceleration of the proton flow is only due to the Laplace force. Our results show, for symmetric and asymmetric connection, the importance of the pressure force, opposed to the electric one on the separatrices, in the decoupling region. In the symmetric case, we emphasize the kinetic origin of this force by analyzing the proton distribution functions and explain their structure by studying the underlying particle dynamics. Protons, as individual particles, are shown to bounce in the electric potential well created by the Hall effect. The spatial divergence of this well results in a mixing in phase space responsible for the observed structure of the pressure tensor. A detailed energy budget analysis confirms the role of the pressure force for the acceleration; but, contrary to what is sometimes assumed, it also reveals that the major part of the incoming Poynting flux is transferred to the thermal energy flux rather than to the convective kinetic energy flux, although the latter is generally supposed dominant. In the symmetric case, we propose the pressure tensor to be an additional proxy of the ion decoupling region in satellite data and verify this suggestion by studying a reconnection event encountered by the Cluster spacecrafts. Finally, the last part of this thesis is devoted to the study of the kinetic structure of asymmetric tangential current sheets where connection can develop. This theoretical part consists in finding a steady state solution to the Vlasov-Maxwell system for the protons in such a configuration. We present the theory and its first confrontation to numerical tests.

Aunai, icolas↗
NASA NTRS2016

DeepSAT: A Deep Learning Approach to Tree-Cover Delineation in 1-m NAIP Imagery for the Continental United States

High resolution tree cover classification maps are needed to increase the accuracy of current land ecosystem and climate model outputs. Limited studies are in place that demonstrates the state-of-the-art in deriving very high resolution (VHR) tree cover products. In addition, most methods heavily rely on commercial softwares that are difficult to scale given the region of study (e.g. continents to globe). Complexities in present approaches relate to (a) scalability of the algorithm, (b) large image data processing (compute and memory intensive), (c) computational cost, (d) massively parallel architecture, and (e) machine learning automation. In addition, VHR satellite datasets are of the order of terabytes and features extracted from these datasets are of the order of petabytes. In our present study, we have acquired the National Agriculture Imagery Program (NAIP) dataset for the Continental United States at a spatial resolution of 1-m. This data comes as image tiles (a total of quarter million image scenes with ~60 million pixels) and has a total size of ~65 terabytes for a single acquisition. Features extracted from the entire dataset would amount to ~8-10 petabytes. In our proposed approach, we have implemented a novel semi-automated machine learning algorithm rooted on the principles of "deep learning" to delineate the percentage of tree cover. Using the NASA Earth Exchange (NEX) initiative, we have developed an end-to-end architecture by integrating a segmentation module based on Statistical Region Merging, a classification algorithm using Deep Belief Network and a structured prediction algorithm using Conditional Random Fields to integrate the results from the segmentation and classification modules to create per-pixel class labels. The training process is scaled up using the power of GPUs and the prediction is scaled to quarter million NAIP tiles spanning the whole of Continental United States using the NEX HPC supercomputing cluster. An initial pilot over the state of California spanning a total of 11,095 NAIP tiles covering a total geographical area of 163,696 sq. miles has produced true positive rates of around 88 percent for fragmented forests and 74 percent for urban tree cover areas, with false positive rates lower than 2 percent for both landscapes.

Imagery↗
DOE OSTI2022

Galaxy Alignments with Surrounding Structure in the Sloan Digital Sky Survey

Using data from the Sloan Digital Sky Survey Legacy Survey, we study the alignment of luminous galaxies with spectroscopic data with the surrounding larger-scale structure as defined by galaxies with only photometric data. We find that galaxies from the red sequence have a statistically significant tendency for their apparent long axes to align parallel to the projected surrounding structure. Red galaxies more luminous than the median of our sample (M r < –21.78) have a mean alignment angle $\langle$Φ$\rangle$ < 45°, indicating preferred parallel alignment, at a significance level >4.5σ on projected scales 0.1 Mpc < r p ≤ 7.5 Mpc. Fainter red galaxies have $\langle$Φ$\rangle$ < 45° at a significance level >4.3σ at scales 1 Mpc < r p < 3 Mpc. At a projected scale r p = 3.0 Mpc, the mean alignment angle decreases steadily with increasing luminosity for red galaxies with M r ≲ –22.5, reaching $\langle$Φ$\rangle$ = 40$^{°}_{.}$49 ± 0$^{°}_{.}$56 for the most luminous 1% (M r ~ –23.57). Galaxies from the blue sequence show no statistically significant tendency for their axes to align with larger-scale structure, regardless of galaxy luminosity. Galaxies in higher-density regions do not show a statistically significant difference in the mean alignment angle from galaxies in lower-density regions; this holds true for the faint blue, luminous blue, faint red, and luminous red subsets.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2001

Carbon Nanotubes in Water: MD Simulations of Internal and External Flow, Self Organization

We have developed computational tools, based on particle codes, for molecular dynamics (MD) simulation of carbon nanotubes (CNT) in aqueous environments. The interaction of CNTs with water is envisioned as a prototype for the design of engineering nano-devices, such as artificial sterocillia and molecular biosensors. Large scale simulations involving thousands of water molecules are possible due to our efficient parallel MD code that takes long range electrostatic interactions into account. Since CNTs can be considered as rolled up sheets of graphite, we expect the CNT-water interaction to be similar to the interaction of graphite with water. However, there are fundamental differences between considering graphite and CNTs, since the curvature of CNTs affects their chemical activity and also since capillary effects play an important role for both dynamic and static behaviour of materials inside CNTs. In recent studies Gordillo and Marti described the hydrogen bond structure as well as time dependent properties of water confined in CNTs. We are presenting results from the development of force fields describing the interaction of CNTs and water based on ab-initio quantum mechanical calculations. Furthermore, our results include both water flows external to CNTs and the behaviour of water nanodroplets inside heated CNTs. In the first case (external flows) the hydrophobic behaviour of CNTs is quantified and we analyze structural properties of water in the vicinity of CNTs with diagnostics such as hydrogen bond distribution, water dipole orientation and radial distribution functions. The presence of water leads to attractive forces between CNTs as a result of their hydrophobicity. Through extensive simulations we quantify these attractive forces in terms of the number and separation of the CNT. Results of our simulations involving arrays of CNTs indicate that these exhibit a hydrophobic behaviour that leads to self-organising structures capable of trapping water clusters. In the second case (internal flows) we study the behaviour of water droplets confined inside CNTs. Constant temperature simulations allow us to capture structural properties such as the contact angles and density profiles of the equilibrated drops. By heating and subsequently cooling of the CNT, we are able to measure the evaporation and the condensation rate of the entrapped water.

Jaffe, Richard L.↗
DOE OSTI2021

Life cycle of streaks in the buffer layer of wall-bounded turbulence

In this work, streaks in the buffer layer of wall-bounded turbulence are tracked in time to study their life cycle. Spatially and temporally resolved direct numerical simulation data are used to analyze the strong wall-parallel movements conditioned to low-speed streamwise flow. The analysis of the streaks shows that there is a clear distinction between wall-attached and detached streaks, and that the wall-attached streaks can be further categorized into streaks that are contained in the buffer layer and the ones that reach the outer region. The results reveal that streaks are born in the buffer layer, coalescing with each other to create larger streaks that are still attached to the wall. Once the streak becomes large enough, it starts to meander due to the large streamwise-to-wall-normal aspect ratio, and consequently the elongation in the streamwise direction, which makes it more difficult for the streak to be oriented strictly in the streamwise direction. While the continuous interaction of the streaks allows the superstructure to span extremely long temporal and length scales, individual streak components are relatively small and short-lived. Tall-attached streaks eventually split into wall-attached and wall-detached components. These wall-detached streaks have a strong wall-normal velocity away from the wall, similar to ejections or bursts observed in the literature. Conditionally averaging the flow fields to these split events show that the detached streak has not only a larger wall-normal velocity compared to the wall-attached counterpart, it also has a larger (less negative) streamwise velocity, similar to the velocity field at the tip of a vortex cluster.

42 ENGINEERING↗
NASA NTRS2019

Accelerated Simulation of Air Pollution Using NVIDIA RAPIDS

Atmospheric chemistry models are a central tool to study and forecast the impact of air pollution on the environment, vegetation, and human health. However, the numerical simulation of chemical kinetics is computationally expensive due to the stiffness of the system of ordinary differential equations that describes atmospheric chemistry. Here we present an alternative approach to the computation of atmospheric chemistry based on machine learning. Our training data set is produced using the NASA Goddard Earth Observing System (GEOS) model with GEOS-Chem chemistry, run on the NASA Center for Climate Simulation (NCCS) Discover supercomputing cluster on 384 Intel Xeon Haswell cores. This model spends more than 50% of total run time on solving atmospheric chemistry. The data set contains as input features the air pollution concentrations before solving the differential equations, together with some key physical parameters such as temperature and sun intensity. As target variables we define the air pollution concentrations after solving the differential equations. Using Dask-cuDF and Dask-XGBoost on the NVIDIA RAPIDS platform on 8 Tesla V100 GPUs, we generate from this training set gradient boosted decision tree models that can reproduce the simulation of chemical kinetics. We do this on the NCCS Advanced Data Analytics Platform (ADAPT) science cloud environment. Our application takes full advantage of recent advances in Dask-XGBoost, such as multi-node and multi-GPU scaling for distributed training with large data sets. The increase in training data size enabled by this is critical to capture the full range of chemical environments encountered across the globe and all annual seasons.The boosted tree models offer good predictability and show many of the features of the full chemistry reference simulation. Further improvements can be achieved through mass balance considerations and by accounting for error correlations. We incorporate the boosted tree models into the GEOS reference model using XGBoost's C API. This enables a seamless integration of the GPU trained models into GEOS-Chem, which is written in Fortran and optimized for use in a massively parallel CPU environment. We show the benefits of this approach and discuss the potential speedup of this machine learning accelerated atmospheric chemistry model.

Keller, Christoph A.↗
DOE OSTI2022

Machine learning–driven multiscale modeling reveals lipid-dependent dynamics of RAS signaling proteins

Significance Here we present an unprecedented multiscale simulation platform that enables modeling, hypothesis generation, and discovery across biologically relevant length and time scales to predict mechanisms that can be tested experimentally. We demonstrate that our predictive simulation-experimental validation loop generates accurate insights into RAS-membrane biology. Evaluating over 100,000 correlated simulations, we show that RAS–lipid interactions are dynamic and evolving, resulting in: 1) a reordering and selection of lipid domains in realistic eight-lipid bilayers, 2) clustering of RAS into multimers correlating with specific lipid fingerprints, 3) changes in the orientation of the RAS G-domain impacting its ability to interact with effectors, and 4) demonstration that RAS–RAS G-domain interfaces are nonspecific in these putative signaling domains.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Equatorial Electronic Structure in the Uranyl Ion: Cs 2 UO 2 Cl 4 and Cs 2 UO 2 Br 4

Electric field gradient (EFG) tensors in the equatorial plane of the linear UO 2 2+ ion have been measured by nuclear magnetic resonance (NMR) and nuclear quadrupole resonance (NQR) experiments and computed by relativistic Kohn–Sham methods with and without environment embedding for Cs 2 UO 2 Cl 4 and Cs 2 UO 2 Br 4 . This approach expands the possibilities for probing the electronic structure in uranyl complexes beyond the strongly covalent U–O bonds. The combined analyses find that one of the two largest principal EFG tensor components at the halogen sites points along the U–X bond (X = Cl, Br), and the second is parallel to the UO 2 2+ ion; in Cs 2 UO 2 Cl 4 , the components are nearly equal in magnitude, whereas in Cs 2 UO 2 Br 4 , due to short-range bromide–cesium interactions, the equatorial component is dominant for one pair of Br sites and the axial component is larger for the second pair. Additionally, the directions and relative magnitudes of the field gradient principal axes are found to be sensitive to the σ and π electron donation by the ligands and the model of the environment. Chlorine-35 NQR spectra of 235 U-depleted and 235 U-enriched Cs 2 UO 2 Cl 4 exhibited no uranium-isotope-dependent shift, but the resonance of the depleted sample displayed a 58% broader line width.

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

Identification of 2-Hydroxyacyl-CoA Synthases with High Acyloin Condensation Activity for Orthogonal One-Carbon Bioconversion

One-carbon (C1) compounds are emerging as cost-effective and potentially carbon-negative feedstocks for biomanufacturing, which require efficient, versatile metabolic platforms for the synthesis of value-added products. Synthetic formyl-CoA elongation (FORCE) pathways allow diverse product synthesis from C1 compounds via iterative C1 elongation, operating independently from the host metabolism with reduced engineering complexity and improved theoretical yields. However, a major bottleneck was identified as the suboptimal kinetics of the core C1–C1 condensation enzyme, 2-hydroxyacyl-CoA synthase (HACS), catalyzing the acyloin condensation reaction between formaldehyde and formyl-CoA. Furthermore, we used a combinatorial approach of bioprospecting and rational protein engineering to identify multiple HACS variants with significantly improved activities toward C1 substrates. Sequence and structure alignment of the active variants elucidated the key regions for the catalytic function, which were targeted for mutagenesis, leading to improved catalytic efficiency. In parallel, a consecutive round of bioprospecting for homologs with high similarity with active variants revealed a highly active HACS variant exhibiting up to 7-fold improvement in catalytic efficiency (k cat /K M ) and 14-fold improvement in the FORCE pathway flux in vivo compared to the previous reports. Upon further optimization of the downstream pathway, the orthogonal C1-to-product bioconversion system showed a metabolic flux of up to 700 μM glycolate OD –1 h –1 (2.1 mmol gDCW –1 h –1 ) and an industrially relevant glycolate titer, rate, and yield of 5.2 g L –1 (67.8 mM), 0.22 g L –1 h –1 , and 94% carbon yield, respectively.

2-hydroxyacyl-CoA synthase↗
NASA NTRS2017

Washboard Terrain on Pluto

Washboard texture or patterning consists of fields of parallel to sub-parallel ridges typically spaced ~1-2 km crest to crest and a few 100 m in amplitude (Fig. 4a in Moore et al., 2016, Science, 351, 1284-1293). For the most part, underlying topography can be easily discerned. We will refer to discrete, well-bounded patches of these landforms as Washboard Terrain (WT). WT is observed to occur along the rim, and just beyond the rim, of Sputnik basin from the West to NNW. Where it is seen in high-resolution data, it has clearly defined limits, beyond which it would be able to be seen if it were there. WT doesn't occur at very low latitudes or very high latitudes (ranging from 22degN to 62degN). WT seems to occur most conspicuously on relatively level, gently sloping terrain. It is restricted to elevations between approximately 2 km to less than +1.5 km (i.e. not at high elevations). The most noticeable regional aspect of the area in which WT occurs is the sinuous valley network, which is suspected to have been formed, or at least substantially modified, by glaciation. WT also appears to occur mainly on an intermediate-albedo reddish material, where seen in enhanced color data. Where it occurs in level terrain, WT tends to trend ENE - there doesn't seem to be a strong local control of its orientation in response to valley drainage directions. WT can display a greater range of orientations where it occurs in higher-relief (not higher elevation) settings such as spurs. WT appears superposed on very ancient landscapes, but is itself cratered locally by clusters of small (approximately 1-3 km) craters, which may be secondaries. This implies that WT may be intermediate in age. Of several working hypotheses, we currently provisionally favor that WT may be akin to terrestrial recessional moraines (or de Geer moraines) associated with the retreat of a higher stand of N2 glaciation that once overfilled Sputnik basin. These putative moraine features may owe their spacing to superseasonal retreat on Milankovitch timescales of approximately 1 Ma. If this hypothesis has validity, then perhaps the intermediate-albedo reddish material may be akin to ground moraine deposits.

Moore, Jeffrey M.↗
NASA NTRS2011

Magnetic Reconnection: Theoretical and Observational Perspectives: Preface

Magnetic reconnection is a fundamental plasma-physical process by which energy stored in a magnetic field is converted, often explosively, into heat and the kinetic energy of the charged particles that constitute the plasma. It occurs in a variety of astrophysical settings, ranging from the solar corona to pulsar magnetospheres and winds, as well as in laboratory fusion experiments, where it is responsible for sawtooth crashes. First proposed by R.G. Giovanelli in the late I 940s as the mechanism responsible for solar flares, magnetic reconnection was invoked at the beginning of the space age to explain not just solar flares but also the transfer of energy, mass, and momentum from the solar wind to Earth's magnetosphere and the subsequent storage and release of the transferred energy in the magnetotai\. During the half century or so that has followed the seminal theoretical works by J.W. Dungey, P.A. Sweet, E.N. Parker, and H.E. Petschek, in-situ measurements by Earth-orbiting satellites and remote-sensing observations of the solar corona have provided a growing body of evidence for the occurrence of reconnection at the Sun, in the solar wind, and in the near-Earth space environment. The last thirty years have also seen the development of laboratory reconnection experiments at a number of institutions. In parallel with the efforts of experimentalists in both space and laboratory plasma physics, theorists have investigated, analytically and with the help of increasingly powerful MHD, hybrid, and kinetic numerical simulations, the structure of the diffusion region, the factors controlling the rate, onset, and cessation of reconnection, and the detailed physics that enables the demagnetization of the ions and electrons and the topological reconfiguration of the magnetic field. Moreover, the scope of theoretical reconnection studies has been extended well beyond solar system and laboratory plasmas to include more exotic astrophysical plasma systems whose strong (10(exp 14)-10(exp 15) G) magnetic fields require that models of reconnection in these systems incorporate quantum electrodynamical, special relativistic, and radiative effects. The papers collected in this topical issue of Space Science Reviews cover different aspects of recent theoretical and observational work on magnetic reconnection in solar and space physics, astrophysics, and laboratory plasma physics. They derive from presentations given at a workshop on magnetic reconnection held in the Yosemite National Park, February 8-12,2010. The intent of the workshop was to stimulate, through a combination of tutorial talks, shorter focused talks, and extensive informal discussions, an interdisciplinary dialogue among members of the different research communities working on the problem of magnetic reconnection. One of the motivating considerations for holding the workshop was its relevance to NASA's Magnetospheric Multiscale (MMS) mission, scheduled for launch in 2014. The four identically instrumented MMS spacecraft are designed to study reconnect ion in Earth's magnetosphere and, specifically, to probe the electron diffusion region in order to determine the microphysical processes that enable the change in the topology of the magnetic field. Building on the achievements of the multi spacecraft Cluster and THEMIS missions, MMS will use the magnetosphere as an astrophysical plasma laboratory in which to test, through in-situ measurement of the plasma, energetic particles, and electric and magnetic fields, various models and theories that have emerged during the past twenty years, a period of extraordinarily productive theoretical and observational work.

Lewis, W. S.↗
NASA NTRS1997

A 2000 Solar Mass Rotating Molecular Disk Around NGC 6334A

We present millimeter and centimeter wave spectroscopic observations of the H II region NGC 6334A. We have mapped the source in several transitions of CO, CS, and NH3. The molecular emission shows a distinct flattened structure in the east-west direction. This structure is probably a thick molecular disk or torus (2.2 x 0.9 pc) responsible for the bipolarity of the near-infrared (NIR) and radio continuum emission which extends in two "lobes" to the north and south of the shell-like H II region. The molecular disk is rotating from west to east (omega approximately equals 2.4 km/s.pc) about an axis approximately parallel to the radio and NIR emission lobes. By assuming virial equilibrium, we find that the molecular disk contains approximately 2000 solar mass. Single-component gas excitation model calculations show that the molecular gas in the disk is warmer and denser (T(sub k) approximately equals 60 K, n approximately equals 3000/cc) than the gas to the north and south (T(sub k) approximately equals 50 K, n approximately equals 400/cc). High resolution (approximately 5 sec) NH3 (3, 3) images of NGC 6334A reveal several small (approximately 0.1 pc) clumps, one of which lies southwest of the radio continuum shell, and is spatially coincident with a near-infrared source, IRS 20. A second NH3 clump is coincident with an H2O maser and the center of a molecular outflow. The dense gas tracers, CS J = 5 approaches 4 and 7 approaches 6, peak near IRS 20 and the H2O maser, not at NGC 6334A. IRS 20 has a substantial far-infrared (FIR) luminosity L(sub FIR) approximately 10(exp 5) solar luminosity, which indicates the presence of an O 7.5 star but has no detected radio continuum (F(sub 6 cm) < 0.02 Jy). The combination of dense gas, a large FIR luminosity and a lack of radio continuum can best be explained if IRS 20 is a protostar. A third clump of NH3 emission lies to the west of IRS 20 but is not associated with any other molecular or continuum features. The star formation activity in the region has moved west of NGC 6334A to IRS 20 and the H2O maser position. We suggest that NGC 6334A, IRS 20, and the H2O maser spot are part of a "protocluster" of stars which is condensing from the massive molecular disk. The similarity between the structure around NGC 6334A and other large (r approximately 1 pc), massive (M approximately 10(exp 3) solar mass), rotating disks (K3-50A and G10.6-0.4) suggests that this may be a common mechanism by which open clusters form.

Kraemer, Kathleen E.↗
DOE OSTI2022

Avoiding excess computation in asynchronous evolutionary algorithms

Abstract Asynchronous evolutionary algorithms are becoming increasingly popular as a means of making full use of many processors while solving computationally expensive search and optimization problems. These algorithms excel at keeping large clusters fully utilized, but may sometimes inefficiently sample an excess of fast‐evaluating solutions at the expense of higher‐quality, slow‐evaluating ones. We have previously introduced a steady‐state parent selection strategy, SWEET (“Selection whilE EvaluaTing”), that sometimes selects individuals that are still being evaluated and allows them to reproduce early. We perform a takeover‐time analysis that confirms that this strategy gives slow‐evaluating individuals that have higher fitnesses an increased ability to multiply in the population. We also find that SWEET appears effective at improving optimization performance on problems in which solution quality is positively correlated with evaluation time. We evaluate our approach on six simulated real‐valued optimization problems and three real‐world applications: an autonomous vehicle controller problem that involves tuning a spiking neural network and two adversarial EA problems. We further evaluate SWEET versus a basic asynchronous process in a simulated setting. We present evidence that SWEET outperforms basic asynchronous processes in a use‐case in which performance is positively correlated with evaluation time, and performs comparably (and often better) than basic asynchronous processes in several use‐cases where performance is negatively correlated with evaluation time. That said, in the cases where performance and evaluation time are negatively correlated the variance of outcomes for SWEET is notably high.

97 MATHEMATICS AND COMPUTING↗