Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “accelerator modelling and simulations”

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 469 records · Page 26

Dark matter, dark energy and fundamental acceleration

We discuss the existence of an acceleration scale in galaxies and galaxy clusters and its relevance for the nature of dark matter. The presence of the same acceleration scale found at very different length scales, and in very different astrophysical objects, strongly supports the existence of a fundamental acceleration scale governing the observed gravitational physics. We comment on the implications of such a fundamental acceleration scale for constraining cold dark matter models as well as its relevance for structure formation to be explored in future numerical simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Preliminary Design of a LAMP DTL (Update)

Present report is a continuation of the previously published report LA-UR-23-26600 from 19-Jun-2023 and incorporates the mentioned report as its part. All the additional information is contained in the part III of this report, the DTL Transverse Focusing Scheme. A preliminary Drift Tube Linac (DTL) layout was designed to create an algorithm for developing the conceptual design of the proposed LANSCE Accelerator Modernization Project (LAMP) final section: the proton linear accelerator from 3 MeV to 100 MeV. Preceding reports describe the proposed layouts of the LEBT, the RFQ, and the MEBT subsections of the linac. Initial estimates of the needed RF power for the DTL are also presented in the report. Present report describes the initial layout of the DTL and the first longitudinal beam dynamics (BD) results of the simulations in simplified models. The part III of this report includes the transverse focusing scheme and transverse BD in details.

43 PARTICLE ACCELERATORS↗

Modeling and simulation of multiphase flows

This presentation provides an overview of the National Energy Technology Laboratory’s (NETL) multiphase computational fluid dynamics codes. The highly successful Multiphase Flows with Interphase eXchanges (MFIX) suite has been used to model a wide range of applications including post-combustion carbon capture, bioreactor optimization, and bio-FCC regeneration. MFIX-Exa, a state-of-the-art CFD code, developed under DOE’s Exascale Computing Project, is built on the AMReX software framework (https://amrex-codes.github.io/) and is designed to leverage modern accelerator-based compute architectures. This presentation further reviews the underlying physical models of both MFIX and MFIX-Exa and contrasts their similarities and differences. Examples of past and present CFD simulations will illustrate how scientific computing at NETL is being used not only for scientific exploration but also for design, optimization and scale-up of multiphase flow devices.

Musser, Jordan [NETL]↗

Initial Steps Towards a Clinical FLASH Radiotherapy System: Pediatric Whole Brain Irradiation with 40 MeV Electrons at FLASH Dose Rates

In this work, we investigated the delivery of a clinically acceptable pediatric whole brain radiotherapy plan at FLASH dose rates using two lateral opposing 40-MeV electron beams produced by a practically realizable linear accelerator system. The EGSnrc Monte Carlo software modules, BEAMnrc and DOSXYZnrc, were used to generate whole brain radiotherapy plans for a pediatric patient using two lateral opposing 40-MeV electron beams. Electron beam phase space files were simulated using a model of a diverging beam with a diameter of 10 cm at 50 cm SAD (defined at brain midline). The electron beams were collimated using a 10-cm-thick block composed of 5 cm of aluminum oxide and 5 cm of tungsten. For comparison, a 6-MV photon plan was calculated with the Varian AAA algorithm. Electron beam parameters were based on a novel linear accelerator designed for the PHASER system and powered by a commercial 6-MW klystron. Calculations of the linear accelerator's performance indicated an average beam current of at least 6.25 µA, providing a dose rate of 115 Gy/s at isocenter, high enough for cognition-sparing FLASH effects. The electron plan was less homogenous with a homogeneity index of 0.133 compared to the photon plan's index of 0.087. Overall, the dosimetric characteristics of the 40-MeV electron plan were suitable for treatment. In conclusion, Monte Carlo simulations performed in this work indicate that two lateral opposing 40-MeV electron beams can be used for pediatric whole brain irradiation at FLASH dose rates of >115 Gy/s and serve as motivation for a practical clinical FLASH radiotherapy system, which can be implemented in the near future.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Particle-in-cell Simulations of Relativistic Magnetic Reconnection with Advanced Maxwell Solver Algorithms

Abstract Relativistic magnetic reconnection is a nonideal plasma process that is a source of nonthermal particle acceleration in many high-energy astrophysical systems. Particle-in-cell (PIC) methods are commonly used for simulating reconnection from first principles. While much progress has been made in understanding the physics of reconnection, especially in 2D, the adoption of advanced algorithms and numerical techniques for efficiently modeling such systems has been limited. With the GPU-accelerated PIC code WarpX, we explore the accuracy and potential performance benefits of two advanced Maxwell solver algorithms: a nonstandard finite-difference scheme (CKC) and an ultrahigh-order pseudo-spectral method (PSATD). We find that, for the relativistic reconnection problem, CKC and PSATD qualitatively and quantitatively match the standard Yee-grid finite-difference method. CKC and PSATD both admit a time step that is 40% longer than that of Yee, resulting in a ∼40% faster time to solution for CKC, but no performance benefit for PSATD when using a current deposition scheme that satisfies Gauss’s law. Relaxing this constraint maintains accuracy and yields a 30% speedup. Unlike Yee and CKC, PSATD is numerically stable at any time step, allowing for a larger time step than with the finite-difference methods. We found that increasing the time step 2.4–3 times over the standard Yee step still yields accurate results, but it only translates to modest performance improvements over CKC, due to the current deposition scheme used with PSATD. Further optimization of this scheme will likely improve the effective performance of PSATD.

79 ASTRONOMY AND ASTROPHYSICS↗

Modelling the development of capillary pressure in freshly 3D-printed concrete elements

3D concrete printing is a promising technology recently developed to automate construction. Since no formwork is used in this technology to support and protect fresh concrete, there are two aspects which considerably accelerate the development of capillary pressure in 3D-printed concrete in comparison to conventionally placed concrete: i) high stiffness of 3D-printed needed to provide sufficient buildability, and ii) very early and fast evaporation of pore water. Accelerated development of capillary pressure may lead to severe plastic shrinkage cracking in 3D-printed elements and, hence, need to be mitigated. This investigation aims at providing a poromechanical model for capillary pressure development in 3D-printed elements. To simulate the development of capillary pressure and plastic shrinkage, environmental factors, material properties, and element geometry need to be considered as a whole. The model inputs – coefficient of permeability, static bulk modulus, air entry pressure and chemical shrinkage rate – were determined experimentally. The model was validated for two fine-grained concretes. Both 3D-printed materials yielded faster capillary pressure increase in comparison to cast concrete, while partial substitution of cement with silica fume further accelerated the capillary pressure development. Furthermore, due to the lower permeability of the mixture containing silica fume, the gradient of capillary pressure between 3D-printed layers increased, as did the gradient of plastic shrinkage.

36 MATERIALS SCIENCE↗

A Feasibility Study on Use of an Accelerometer to Measure the Dynamic Forces in Turning

Metal cutting is a highly dynamic process that generates continuously varying forces. Measurement of these forces is essential to characterizing a cutting process and to fully utilize the machine tool. A force dynamometer is often used to measure these varying forces; however, it is limited by the natural frequency of the sensor and can sometimes be hard to set up on a machine. Therefore, this study aims to estimate the dynamic component of the cutting force using an accelerometer placed directly below the cutting edge. The placement of the sensor helps achieving better signal to noise ratio (S/N). The cutting tool is modelled as a single degree of freedom (SDOF) system relative to the workpiece and dynamic loads are experimentally simulated on it. The measured acceleration is numerically integrated to obtain velocity and deflection, which along with the natural frequency (ω n ) and damping ratio (ζ), is used to predict the dynamic load using the system equation. The dynamic forces tested over a wide range of frequencies, show good agreement with the data from simultaneous dynamometer measurement and the force estimated using the proposed method. The study shows the feasibility of this method in a real cutting scenario and the capability to apply it to a multi DOF system. One major limitation is the inability to capture the static or quasi-static component of the cutting force.

Cutting vibrations↗

The Role of Topography in Controlling Evapotranspiration Age

Abstract Evapotranspiration (ET) age is a key metric of water sustainability but a major unknown partly due to the extreme difficulty in modeling it. Groundwater is found to be important in ET age variations in small‐scale studies, yet our understanding is insufficient because groundwater systems are nested across scales. Here, we conducted GPU‐accelerated particle tracking with integrated hydrologic modeling to quantify the variations in ET age at a regional scale of ∼0.4 M km 2 . Simulation results reveal topography‐driven flow paths shaping the spatial and temporal patterns of ET age variations. On ridges, where root zone decoupling with deep subsurface storage, ET age is generally young, with seasonal variations dominated by meteorological conditions. In the valley bottom, ET age is generally old, with significant subseasonal variations caused by the convergence of subsurface flow paths. On hillslopes with water table depths ranging from 1 to 10 m, ET age shows strong seasonal variations caused by the connections with lateral groundwater regulated by ET demand. Our modeling approach provides insights into the basic linkages between ET age and topography at large scale. Our work highlights the perspective of multiscale studies of ET age, suggesting new field experiments to test these process connections and to determine if such linkages warrant inclusion in Earth System Models.

54 ENVIRONMENTAL SCIENCES↗

Modeling Interactive Effects of Manganese Bioavailability, Nitrogen Deposition, and Warming on Soil Carbon Storage

Manganese (Mn) is a redox–active micronutrient that has been shown to accelerate plant litter decomposition; however, the effect of Mn–promoted decomposition on soil C storage is unclear. We present a novel biogeochemical model simulating how Mn bioavailability influences soil organic C (SOC) stocks in a soil profile (<50 cm) within a temperate forest. In our model, foliar Mn increased in response to increasing soluble Mn released through Mn–oxide (birnessite) dissolution in mineral soil layers. The ensuing Mn enrichment in leaf litter redistributed Mn to the surface forest floor layer, promoted enzymatic oxidation of lignin, and decreased SOC stocks. Total SOC loss was partially mitigated by accumulation of lignin–oxidation products as mineral–associated organic C. We also explored how Mn–driven changes to C storage interacted with effects of N deposition and warming. Nitrogen enrichment inhibited Mn–dependent lignin degradation, increasing SOC stocks and weakening their dependence on Mn bioavailability. Warming stimulated decomposition and reduced C stocks but was less effective at low Mn bioavailability. Our model results suggest that SOC stocks are sensitive to Mn bioavailability because increased plant uptake redistributes Mn to surface soils where it can enhance litter decomposition. Finally, based on our simulations, we predict that Mn becomes limiting to litter decomposition where Mn is poorly soluble. Depletion of bioavailable Mn or other cofactors that are critical to decomposition could limit the response of organic C stocks to warming over time, but quantitative projections of the role of Mn bioavailability in regulating decomposition requires additional measurements to constrain model uncertainties.

58 GEOSCIENCES↗

Data assimilation empowered neural network parametrizations for subgrid processes in geophysical flows

In the past couple of years, there has been a proliferation in the use of machine learning approaches to represent subgrid-scale processes in geophysical flows with an aim to improve the forecasting capability and to accelerate numerical simulations of these flows. Despite its success for different types of flow, the online deployment of a data-driven closure model can cause instabilities and biases in modeling the overall effect of subgrid-scale processes, which in turn leads to inaccurate prediction. To tackle this issue, we exploit the data assimilation technique to correct the physics-based model coupled with the neural network as a surrogate for unresolved flow dynamics in multiscale systems. In particular, we use a set of neural network architectures to learn the correlation between resolved flow variables and the parametrizations of unresolved flow dynamics and formulate a data assimilation approach to correct the hybrid model during their online deployment. We illustrate our framework in a set of applications of the multiscale Lorenz 96 system for which the parametrization model for unresolved scales is exactly known, and the two-dimensional Kraichnan turbulence system for which the parametrization model for unresolved scales is not known a priori. Our analysis, therefore, comprises a predictive dynamical core empowered by (i) a data-driven closure model for subgrid-scale processes, (ii) a data assimilation approach for forecast error correction, and (iii) both data-driven closure and data assimilation procedures. We show significant improvement in the long-term prediction of the underlying chaotic dynamics with our framework compared to using only neural network parametrizations for future prediction. Moreover, we demonstrate that these data-driven parametrization models can handle the non-Gaussian statistics of subgrid-scale processes, and effectively improve the accuracy of outer data assimilation workflow loops in a modular nonintrusive way.

42 ENGINEERING↗

A simulation of the Fermilab Main Injector dual power amplifier cavities

The Fermilab Main Injector accelerating cavities have sparking issues when they are run at voltages higher than those required by the PIP-II project. This is a problem Fermilab is working on as planning begins for the next upgrade to the accelerator complex. One of the methods being used to address the issue is the development of a CST Microwave Studio simulation to accurately model the PIP-II dual power amplifier cavities and identify which part(s) of the cavity is causing sparking to develop. The model will also be used to determine if changes to the cavity geometry may allow the cavity to be used at higher voltages before sparking occurs.

Stevenson, Susanna [Fermilab] (ORCID:0009000891261↗

A Simulation of the Fermilab Main Injector Cavities

The Fermilab Main Injector accelerating cavities have sparking issues when they are run at voltages higher than those required by the PIP-II project. This is a problem Fermilab is working on as planning begins for the next upgrade to the accelerator complex. One of the methods being used to address the issue is the development of a CST Microwave Studio simulation to accurately model the PIP-II dual power amplifier cavities and identify which part(s) of the cavity is causing sparking to develop. The model will also be used to determine if changes to the cavity geometry may allow the cavity to be used at higher voltages before sparking occurs

Stevenson, Susanna M. [Fermilab] (ORCID:0009000891↗

Multi-scale models of whole cells: progress and challenges

Whole-cell modeling is “the ultimate goal” of computational systems biology and “a grand challenge for 21st century”. These complex, highly detailed models account for the activity of every molecule in a cell and serve as comprehensive knowledgebases for the modeled system. Their scope and utility far surpass those of other systems models. In fact, whole-cell models (WCMs) are an amalgam of several types of “system” models. The models are simulated using a hybrid modeling method where the appropriate mathematical methods for each biological process are used to simulate their behavior. Given the complexity of the models, the process of developing and curating these models is labor-intensive and to date only a handful of these models have been developed. While whole-cell models provide valuable and novel biological insights, and to date have identified some novel biological phenomena, their most important contribution has been to highlight the discrepancy between available data and observations that are used for the parametrization and validation of complex biological models. Another realization has been that current whole-cell modeling simulators are slow and to run models that mimic more complex (e.g., multi-cellular) biosystems, those need to be executed in an accelerated fashion on high-performance computing platforms. In this manuscript, we review the progress of whole-cell modeling to date and discuss some of the ways that they can be improved.

59 BASIC BIOLOGICAL SCIENCES↗

Molecular Dynamics Simulation of Complex Reactivity with the Rapid Approach for Proton Transport and Other Reactions (RAPTOR) Software Package

Simulating chemically reactive phenomena such as proton transport on nanosecond to microsecond and beyond time scales is a challenging task. Ab initio methods are unable to currently access these time scales routinely, and traditional molecular dynamics methods feature fixed bonding arrangements that cannot account for changes in the system’s bonding topology. The Multiscale Reactive Molecular Dynamics (MS-RMD) method, as implemented in the Rapid Approach for Proton Transport and Other Reactions (RAPTOR) software package for the LAMMPS molecular dynamics code, offers a method to routinely sample longer time scale reactive simulation data with statistical precision. RAPTOR may also be interfaced with enhanced sampling methods to drive simulations toward the analysis of reactive rare events, and a number of collective variables (CVs) have been developed to facilitate this. Key advances to this methodology, including GPU acceleration efforts and novel CVs to model water wire formation are reviewed, along with recent applications of the method which demonstrate its versatility and robustness.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Scalable Deep-Learning-Accelerated Topology Optimization for Additively Manufactured Materials

Topology optimization (TO) is a popular and powerful computational approach for designing novel structures, materials, and devices. Two computational challenges have limited the applicability of TO to a variety of industrial applications. First, a TO problem often involves a large number of design variables to guarantee sufficient expressive power. Second, many TO problems require a large number of expensive physical model simulations, and those simulations cannot be parallelized. To address these issues, we propose a general scalable deep-learning (DL) based TO framework, referred to as SDL-TO, which utilizes parallel schemes in high performance computing (HPC) to accelerate the TO process for designing additively manufactured (AM) materials. Unlike the existing studies of DL for TO, our framework accelerates TO by learning the iterative history data and simultaneously training on the mapping between the given design and its gradient. The surrogate gradient is learned by utilizing parallel computing on multiple CPUs incorporated with a distributed DL training on multiple GPUs. The learned TO gradient enables a fast online update scheme instead of an expensive update based on the physical simulator or solver. Using a local sampling strategy, we achieve to reduce the intrinsic high dimensionality of the design space and improve the training accuracy and the scalability of the SDL-TO framework. The method is demonstrated by benchmark examples and AM materials design for heat conduction. The proposed SDL-TO framework shows competitive performance compared to the baseline methods but significantly reduces the computational cost by a speed up of around 8.6x over the standard TO implementation.

Bi, Sirui↗

Performance of Michigan sMDT prototype chambers for the HL-LHC ATLAS muon detector upgrade

A new small-diameter Monitored Drift Tube (sMDT) chamber has been developed for the muon spectrometer of the ATLAS experiment to handle the higher collision rates expected at the CERN High Luminosity Large Hadron Collider (HL-LHC). This paper presents measurements of the tracking resolution and hit efficiency of two prototype sMDT chambers constructed at the University of Michigan. Using cosmic-ray muons the sMDT tracking resolution of 103.7 ± 8.1 μm was measured for one chamber and 101.8 ± 7.8 μm for the other, compared with a design resolution of 106 μm. A further tracking resolution improvement to 83.4 ± 7.8 μm was obtained by using new high-gain readout electronics which will be added for HL-LHC. An average tracking efficiency of (98.5 ± 0.2)% was found for both chambers. The methodology used to determine the detector tracking resolution and efficiency, including reconstruction of sMDT data and a Geant4 simulation of the test chamber, is presented in detail.

43 PARTICLE ACCELERATORS↗

Predicting beam transmission using 2-dimensional phase space projections of hadron accelerators

We present a method to compress the 2D transverse phase space projections from a hadron accelerator and use that information to predict the beam transmission. This method assumes that obtaining at least three projections of the 4D transverse phase space is possible and that an accurate simulation model is available for the beamline. Using a simulated model, we show that—a computer can train a convolutional autoencoder to reduce phase-space information which can later be used to predict the beam transmission. Finally, we argue that although using projections from a realistic nonlinear distribution produces less accurate results, the method still generalizes well.

43 PARTICLE ACCELERATORS↗