Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “accelerator simulation”

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 973 records · Page 54

Plateau–Rayleigh instability with a grain boundary twist

We demonstrate using theoretical, computational, and experimental studies a morphological instability, in which a polycrystalline nanorod breaks up at grain boundaries (GBs) into an array of isolated domains. Our theoretical model is used to establish a neutral stability surface demarcating stable and unstable perturbations. It is shown that GBs play a destabilizing role in which the critical wavelength for the instability decreases with the increase in the GB energy. We carry out phase field simulations, which reveal accelerated pinch-off kinetics with the increase in the GB energy and predict temporal evolution of interfacial profiles in quantitative agreement with experimental observations.

Hussein, Omar↗

MHD-induced beta limits in the Large Helical Device

Using the extended-magnetohydrodynamics code, M3D-C1, we perform a systematic numerical study of the effect of externally applied heating on the achievable plasma beta in a ten field-period heliotron. Heat sources of varying intensity are applied to a vacuum magnetic field that is representative of the standard configuration of the Large Helical Device, with R 0 = 3.66 m, where R 0 is the radial position of the magnetic axis in vacuum. As the system is driven to a state that is unstable to low-n magnetohydrodynamic (MHD) modes, nonlinear mode interactions lead to the formation of chaotic magnetic fields. With sufficiently strong heating, a collapse of the electron temperature profile is observed. This demonstrates the necessity of simulating the self-consistent evolution of plasma profiles, without imposing assumptions on the structure of the magnetic field, to accurately determine transport properties in stellarator plasmas. It also highlights the value of these advanced simulation capabilities for accelerating the development of high-performance stellarator operating scenarios.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Acceleration of the matrix multiplication of Radiance three phase daylighting simulations with parallel computing on heterogeneous hardware of personal computer

Building designers are increasingly relying on complex fenestration systems to reduce energy consumed for lighting and HVAC in low energy buildings. Radiance, a lighting simulation program, has been used to conduct daylighting simulations for complex fenestration systems. Depending on the configurations, the simulation can take hours or even days using a personal computer. This paper describes how to accelerate the matrix multiplication portion of a Radiance three-phase daylight simulation by conducting parallel computing on heterogeneous hardware of a personal computer. The algorithm was optimized and the computational part was implemented in parallel using OpenCL. The speed of new approach was evaluated using various daylighting simulation cases on a multicore central processing unit and a graphics processing unit. Based on the measurements and analysis of the time usage for the Radiance daylighting simulation, further speedups can be achieved by using fast I/O devices and storing the data in a binary format.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

SPIKE PENETRATION IN BLAST-WAVE-DRIVEN INSTABILITIES

The problem of interest is the unstable growth of structure at density transitions affected by blast waves, which arise in natural environments such as core-collapse supernovae and in laboratory experiments. The resulting spikes of dense material, which penetrate the less dense material, develop broadened tips, but the degree of broadening varies substantially across both experiments and simulations. The variable broadening presumably produces variations in the drag experienced by the spike tips as they penetrate the less dense material. The present work has used semianalytic theory to address the question of how the variation in drag might affect the spike penetration, for cases in which the post-shock interface deceleration can be described by a power law in a normalized time variable. It did so by following the evolution of structure on the interface through the initial shock passage, the subsequent small-amplitude phase of Rayleigh-Taylor instability growth, and the later phase in which the spike growth involves the competition of buoyancy and drag. In all phases, the expansion of the system during its evolution was accounted for and was important. The calculated spike length is strongly affected by the drag attributed to spike tip broadening. One finds from such a calculation that it is not unreasonable for narrow spikes to keep up with the shock front of the blast wave. The implication is that the accuracy of prediction of spike penetration and consequent structure by simulations very likely depends on how accurately they treat the broadening of the spike tips and the associated drag. Experimental validation of spike morphology in simulations would be useful.

79 ASTRONOMY AND ASTROPHYSICS↗

Supporting shared resource usage for a diverse user community: The OSG experience and lessons learned

The Open Science Grid (OSG) supports a diverse community of new and existing users in adopting and making effective use of the Distributed High Throughput Computing (DHTC) model. The LHC user community has deep local support within the experiments. For other smaller communities and individual users the OSG provides consulting and technical services through the User Support area. We describe these sometimes successful and sometimes not so successful experiences and analyze lessons learned that are helping us improve our services. The services offered include forums to enable shared learning and mutual support, tutorials and documentation for new technology, and troubleshooting of problematic or systemic failure modes. For new communities and users, we bootstrap their use of the distributed high throughput computing technologies and resources available on the OSG by following a phased approach. We first adapt the application and run a small production campaign on a subset of friendly sites. Only then do we move the user to run full production campaigns across the many remote sites on the OSG, adding to the community resources up to hundreds of thousands of CPU hours per day. This scaling up generates new challenges - like no determinism in the time to job completion, and diverse errors due to the heterogeneity of the configurations and environments - so some attention is needed to get good results. We cover recent experiences with image simulation for the Large Synoptic Survey Telescope (LSST), small-file large volume data movement for the Dark Energy Survey (DES), civil engineering simulation with the Network for Earthquake Engineering Simulation (NEES), and accelerator modeling with the Electron Ion Collider group at BNL. We will categorize and analyze the use cases and describe how our processes are evolving based on lessons learned.

Garzoglio, Gabriele↗

Mechanistic studies of small molecule ligands selective to RNA single G bulges

Abstract Small-molecule RNA binders have emerged as an important pharmacological modality. A profound understanding of the ligand selectivity, binding mode, and influential factors governing ligand engagement with RNA targets is the foundation for rational ligand design. Here, we report a novel class of coumarin derivatives exhibiting selective binding affinity towards single G RNA bulges. Harnessing the computational power of all-atom Gaussian accelerated molecular dynamics simulations, we unveiled a rare minor groove binding mode of the ligand with a key interaction between the coumarin moiety and the G bulge. This predicted binding mode is consistent with results obtained from structure-activity relationship studies and transverse relaxation measurements by nuclear magnetic resonance spectroscopy. We further generated 444 molecular descriptors from 69 coumarin derivatives and identified key contributors to the binding events, such as charge state and planarity, by lasso (least absolute shrinkage and selection operator) regression. Our work deepened the understanding of RNA-small molecule interactions and integrated a new framework for the rational design of selective small-molecule RNA binders.

Biochemistry & Molecular Biology↗

Extracting symplectic maps for space-charge dominated beams

Symplecticity of transfer maps is important for reliable evaluation of space-charge dominated beams in accelerators. Unfortunately, most simulation codes that include collective effects, such as space charge, do not use canonical phase-space variables and therefore are not symplectic in the presence of electromagnetic fields. In this paper, we present a numerical method to extract local linear symplectic transfer maps using particle tracking simulation code for space-charge dominated beams. We demonstrate this method for the photoinjector (113 MHz SRF gun) section of the Coherent electron Cooling (CeC) Proof of Principle (POP) experiment.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Ion acceleration from microstructured targets irradiated by high-intensity picosecond laser pulses

Structures on the front surface of thin foil targets for laser-driven ion acceleration have been proposed to increase the ion source maximum energy and conversion efficiency. While structures have been shown to significantly boost the proton acceleration from pulses of moderate-energy fluence, their performance on tightly focused and high-energy lasers remains unclear. Here, we report the results of laser-driven three-dimensional (3D)-printed microtube targets, focusing on their efficacy for ion acceleration. Using the high-contrast (~10 12) PHELIX laser (150 J, 10 21 W / cm 2 ), we studied the acceleration of ions from 1-μm-thick foils covered with micropillars or microtubes, which we compared with flat foils. The front-surface structures significantly increased the conversion efficiency from laser to light ions, with up to a factor of 5 higher proton number with respect to a flat target, albeit without an increase of the cutoff energy. An optimum diameter was found for the microtube targets. Our findings in this work are supported by a systematic particle-in-cell modeling investigation of ion acceleration using 2D simulations with various structure dimensions. Simulations reproduce the experimental data with good agreement, including the observation of the optimum tube diameter, and reveal that the laser is shuttered by the plasma filling the tubes, explaining why the ion cutoff energy was not increased in this regime.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Exascale workflow applications and middleware: An ExaWorks retrospective

Exascale computers offer transformative capabilities to combine data-driven and learning-based approaches with traditional simulation applications to accelerate scientific discovery and insight. However, these software combinations and integrations are difficult to achieve due to the challenges of coordinating and deploying heterogeneous software components on diverse and massive platforms. Here, we present the ExaWorks project, which addresses many of these challenges. We developed a workflow Software Development Toolkit (SDK), a curated collection of workflow technologies that can be composed and interoperated through a common interface, engineered following current best practices, and specifically designed to work on HPC platforms. ExaWorks also developed PSI/J, a job management abstraction API, to simplify the construction of portable software components and applications that can be used over various HPC schedulers. The PSI/J API is a minimal interface for submitting and monitoring jobs and their execution state across multiple and commonly used HPC schedulers. We also describe several leading and innovative workflow examples of ExaWorks tools used on DOE leadership platforms. Furthermore, we discuss how our project is working with the workflow community, large computing facilities, and HPC platform vendors to address the requirements of workflows sustainably at the exascale.

97 MATHEMATICS AND COMPUTING↗

Virtual Observation System for Earth System Model: An Application to ACME Land Model Simulations

Investigating and evaluating physical-chemical-biological processes within an Earth system model (EMS) can be very challenging due to the complexity of both model design and software implementation. A virtual observation system (VOS) is presented to enable interactive observation of these processes during system simulation. Based on advance computing technologies, such as compiler-based software analysis, automatic code instrumentation, and high-performance data transport, the VOS provides run-time observation capability, in-situ data analytics for Earth system model simulation, model behavior adjustment opportunities through simulation steering. A VOS for a terrestrial land model simulation within the Accelerated Climate Modeling for Energy model is also presented to demonstrate the implementation details and system innovations.

58 GEOSCIENCES↗

Reference Model MHK Turbine Array Optimization Study within a Generic River System.

Increasing interest in marine hydrokinetic (MHK) energy has spurred to significant research on optimal placement of emerging technologies to maximize energy conversion and minimize potential effects on the environment. However, these devices will be deployed as an array in order to reduce the cost of energy and little work has been done to understand the impact these arrays will have on the flow dynamics, sediment-bed transport and benthic habitats and how best to optimize these arrays for both performance and environmental considerations. An "MHK-friendly" routine has been developed and implemented by Sandia National Laboratories (SNL) into the flow, sediment dynamics and water-quality code, SNL-EFDC. This routine has been verified and validated against three separate sets of experimental data. With SNL-EFDC, water quality and array optimization studies can be carried out to optimize an MHK array in a resource and study its effects on the environment. The present study examines the effect streamwise and spanwise spacing has on the array performance. Various hypothetical MHK array configurations are simulated within a trapezoidal river channel. Results show a non-linear increase in array-power efficiency as turbine spacing is increased in each direction, which matches the trends seen experimentally. While the sediment transport routines were not used in these simulations, the flow acceleration seen around the MHK arrays has the potential to significantly affect the sediment transport characteristics and benthic habitat of a resource. Evaluation Only. Created with Aspose.Pdf.Kit. Copyright 2002-2011 Aspose Pty Ltd Evaluation Only. Created with Aspose.Pdf.Kit. Copyright 2002-2011 Aspose Pty Ltd

54 ENVIRONMENTAL SCIENCES↗

Crystal growth and nucleation at metal-organic framework / solution interfaces

Metal-organic frameworks (MOFs) represent an incredibly diverse class of materials with equally diverse applications in catalysis, separations, and gas storage. Despite the prevalence of MOF structures in the literature, in many cases the synthesis of pure, crystalline MOF materials remains a significant challenge, with failed syntheses often yielding mixed crystalline or non-porous amorphous products. Our recent work aims to enable targeted MOF synthesis by developing a comprehensive and predictive model of the fundamental processes occurring at the MOF-solution interface under synthetically-relevant conditions. Our recent accomplishments include: (i) development and implementation of a graphics processor-accelerated approach to simulation the nucleation of weak electrolytes, and the extension of these methods to large/complex solutes such as MOFs; and (ii) studies of the structure of ZIF interfaces under conditions relevant to crystal growth; and (iii) addressing important open questions in the mechanism of MOF solvent-assisted linker exchange, and important synthetic pathway for MOFs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fast multipole methods for multiparticle simulations

The fast multipole method (FMM) reduces the computation cost of the pairwise non-oscillating interaction between N particles from O(N^2) to O(N). In the FMM, the contribution from a source particle is represented as a multipole expansion, while the contributions from multiple faraway sources can be combined into a local expansion around an objective particle. Without the dependence on a grid covering the whole domain under study, the FMM treats any charge distribution and geometry in a natural way. It avoids artificial smoothing due to the grid size and redundant computation on the free space grids. We will introduce the concept of the FMM using the Coulomb interaction as an example and then explain how the FMM can be extended to arbitrary non-oscillating interactions. Examples and discussions on how the FMM can be used in scientific simulations, especially in accelerator physics will also be provided.

Zhang, He↗

E3SM: Improved Climate Prediction with Exascale Capability

The Energy Exascale Earth System Model (E3SM) project is an ongoing, state-of-the-science earth system modeling, simulation, and prediction effort that optimizes Department of Energy (DOE) computing resources to meet the science needs of the nation and the agency’s mission objectives. Climate simulation has become a proven tool for identifying and quantifying the impacts of climate change, but even greater accuracy is required at all levels to improve forecast precision. Understanding the impact of climate change on global and regional water cycles is one of the highest priorities and most difficult challenges in climate change prediction. As part of a subproject of DOE’s Exascale Computing Project, a multidisciplinary team including geophysical and computational scientists developed a multiscale modeling framework (MMF) to refine cloud representation in E3SM climate simulation on GPU accelerated supercomputers, making higher resolution, more computationally efficient predictions possible.

54 ENVIRONMENTAL SCIENCES↗

Web-Based Tools for Data-Informed Remedy Optimization: Software Theory and User Guide

This report documents the development and application of two web-based decision-support tools for pump-and-treat (P&T) groundwater remediation systems: PTOLEMY (Pump-and-Treat Optimized Location Evaluation to Maximize Yields) and OPTIMA (Optimization for Pump-and-Treat Implementation, Management, & Assessment). These tools enhance remedy design and management by leveraging advanced computational methods – specifically deep learning and multi-objective optimization – within a user-friendly platform. By integrating data-driven models with established hydrogeological knowledge, PTOLEMY and OPTIMA enable more efficient evaluation of well placement and operational strategies, helping site managers balance multiple remediation objectives under complex conditions. Both tools are implemented as modules within the SOCRATES (Suite Of Comprehensive Rapid Analysis Tools for Environmental Sites) web platform, which provides data access, visualization, and analytics to support remedy optimization across sites in the U.S. Department of Energy Office of Environmental Management complex. PTOLEMY is a rapid screening module designed to identify promising locations for new extraction wells. It employs a multi-channel three-dimensional convolutional neural network (MC3D-CNN) trained on high-fidelity simulation data to predict the relative performance (in terms of contaminant mass recovery) of potential well sites. Through an interactive web interface, PTOLEMY visualizes the probability of high performance across a site, highlighting areas where an extraction well is likely to yield above-threshold contaminant removal over a multi-year period. PTOLEMY’s map-based displays and exportable results support transparent communication of screening analyses. By focusing attention on the most favorable candidate locations, the tool augments traditional engineering judgment and physics-based modeling, providing a data informed basis for subsequent detailed evaluations. OPTIMA is a multi objective optimization module designed to find wellfield layouts and operating schedules that meet various cleanup goals. It quickly evaluates thousands of candidate setups – combinations of well locations, timing, and rates – and returns a small set of best trade-off options for comparison. At its core, OPTIMA uses a U-Net-based surrogate model – a deep-learning emulator of a groundwater flow and transport simulator – to dramatically accelerate scenario evaluations. Coupling this fast surrogate with the NSGA-II (Non-dominated Sorting Genetic Algorithm II) evolutionary algorithm, OPTIMA explores a wide decision space of well locations and schedules to identify Pareto-optimal solutions that trade off key objectives (e.g., minimizing cleanup time, maximizing contaminant mass removal, and minimizing plume extent). The tool outputs a family of optimal configurations and visualizes their trade-offs (Pareto frontiers of cleanup metrics and maps of optimized well placements). Site managers can use these results to understand the range of viable strategies and to select candidate designs for more detailed verification. OPTIMA is currently under active development and not yet fully released; this guide provides early documentation to support planning and gather user feedback.

54 ENVIRONMENTAL SCIENCES↗

Localized Exponential Time Differencing Method for Shallow Water Equations: Algorithms and Numerical Study

Here, we explore the performance of the exponential time differencing (ETD) method applied to the rotating shallow water equations. Comparing with explicit time stepping of the same order accuracy in time, the ETD algorithms could reduce the computational time in many cases by allowing the use of large time step sizes while still maintaining numerical stability. To accelerate the ETD simulations, we propose a localized approach that synthesizes the ETD method and overlapping domain decomposition. By dividing the original problem into many subdomain problems of smaller sizes and solving them locally, the proposed approach could speed up the calculation of matrix exponential vector products. Several standard test cases for shallow water equations of one or multiple layers are considered. The results show great potential of the localized ETD method for high-performance computing because each subdomain problem can be naturally solved in parallel at every time step.

58 GEOSCIENCES↗