Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parallel beam dynamics 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 37 records · Page 2

RLGBS: Reinforcement Learning-Guided Beam Search for process optimization in a paper machine dryer section

Paper drying is responsible for over two-thirds of energy consumption in the U.S. pulp and paper industry, presenting significant potential for energy savings through optimization of process parameters. Current approaches often assume fixed operating conditions, neglecting dynamic ambient and process variations that limit achievable savings and real-world applicability. To this end, we develop a physics-based simulation environment for a paper machine dryer section and propose a reinforcement learning (RL) framework to minimize overall energy consumption by optimizing drying process parameters under diverse operating conditions. To mitigate overdrying and numerical instabilities caused by suboptimal local RL actions, we introduce Reinforcement Learning-Guided Beam Search (RLGBS), which explores multiple action sequences in parallel using beam search. Instead of making step-by-step decisions, RLGBS prioritizes solutions based on cumulative probability, reducing the impact of individual suboptimal actions. Experiments demonstrate that RLGBS achieves consistent energy savings under unseen operating conditions not encountered during training, outperforming conventional RL methods. While validated in drying optimization, this framework is broadly applicable to other RL-based industrial process control problems.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

UQpy: A general purpose Python package and development environment for uncertainty quantification

In this paper, we present the UQpy software toolbox, an open-source Python package for general uncertainty quantification (UQ) in mathematical and physical systems. The software serves as both a user-ready toolbox that includes many of the latest methods for UQ in computational modeling and a convenient development environment for Python programmers advancing the field of UQ. The paper presents an introduction to the software's architecture and existing capabilities, divided in the code in a set of modules centered around different UQ tasks such as sampling methods, generation of random processes and random fields, probabilistic inverse modeling, reliability analysis, surrogate modeling, and active learning. The paper also highlights the importance of the RunModel module, which is used to drive simulations in the uncertainty analyses performed in UQpy. This module conveniently allows the user to define computational models directly in Python, or to run simulations from a third-party software in serial or in parallel. To illustrate the various capabilities, two examples are tracked throughout the paper and analyzed repeatedly for various UQ tasks. The first is a Python model solving a nonlinear structural dynamics problem, used to illustrate UQpy's capabilities in sampling and forward propagation of high dimensional random vectors (stochastic processes), and probabilistic inference. The second model is a third-party Abaqus finite element model solving the thermomechanical response of a beam structure. This example is used to illustrate UQpy's capabilities in variance reduction sampling techniques, reliability analysis, surrogate modeling and active learning techniques.

97 MATHEMATICS AND COMPUTING↗

Dynamics of a split torque helicopter transmission

A high reduction ratio split torque gear train has been proposed as an alternative to a planetary configuration for the final stage of a helicopter transmission. A split torque design allows a high ratio of power-to-weight for the transmission. The design studied in this work includes a pivoting beam that acts to balance thrust loads produced by the helical gear meshes in each of two parallel power paths. When the thrust loads are balanced, the torque is split evenly. A mathematical model was developed to study the dynamics of the system. The effects of time varying gear mesh stiffness, static transmission errors, and flexible bearing supports are included in the model. The model was demonstrated with a test case. Results show that although the gearbox has a symmetric configuration, the simulated dynamic behavior of the first and second compound gears are not the same. Also, results show that shaft location and mesh stiffness tuning are significant design parameters that influence the motions of the system.

Rashidi, Majid↗

The effects of multispecies Hohlraum walls on stimulated Brillouin scattering, Hohlraum dynamics, and beam propagation

We report experiments and simulations have been conducted to investigate the efficacy of Ta 2 O 5 -lined Hohlraum walls at reducing stimulated Brillouin backscattering (SBS) as well as any subsequent effects on the Hohlraum dynamics and capsule implosions in indirect drive experiments at the National Ignition Facility. Using a 1.1 MJ 400 TW, 351 nm, shaped laser pulse, we measure a 5× reduction in SBS power in the peak of the pulse from the wall on the outer 50° cone beams. The SBS spectrum indicates a reduction in the high-Z spectral signature when using multispecies wall materials. Detailed hydrodynamic simulations were performed using different heat conduction models with flux limiters. Additional simulations were performed on the plasma maps using the 3D parallel paraxial code pF3D to compare backscatter powers between the pure Au and Ta 2 O 5 -lined Hohlraums. Further analysis, using hydrodynamically equivalent plasmas, shows that the SBS reduction is clearly a result of the added ion Landau damping caused by the oxygen ions and not from differences in plasma conditions. The experimental and simulation results also show an increase in the wall plasma expansion when using the Ta 2 O 5 liner leading to a 70% more oblate implosion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Steerable terahertz beams using surface waves on an active metasurface

The development of dynamic components for controlling wave fronts in the sub-terahertz region of the electromagnetic spectrum has emerged as a frontier research topic for many applications in sensing and communications. One approach which has attracted much attention involves the use of active metasurfaces, tiled arrays of sub-wavelength elements with properties that can be reconfigured via external actuation. In nearly all cases, these metasurfaces are employed as either transmissive or reflective elements, taking advantage of their strong and tunable interaction with free-space electromagnetic waves. These interactions can be significantly enhanced through the use of surface waves propagating parallel to the metasurface array, although very few studies have exploited this option. Here, we integrate a metasurface into the interior of a parallel-plate waveguide in a configuration explicitly designed to exploit this surface-wave geometry. We show that varying the electrical properties of the active metasurface changes the wave vector of the guided mode, and thereby alters the emission angle of radiation out-coupled through a leaky-wave slot aperture. These results, which are consistent with numerical simulations, represent a new approach to broadband beam steering suitable for the sub-terahertz spectral range.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

O(+) acceleration due to resistive momentum transfer in the auroral field line plasma

An analytical model is defined to demonstrate that parallel acceleration of an O(+) ion beam in the ionosphere can occur naturally due to the presence of a quasi-static parallel electric field. Momentum equations are defined for friction between hydrogen ions and electrons, which produces a quasi-static electric field. The field can accelerate ions, e.g., the O(+) ions, which do not participate in the frictional momentum exchange. The conditions are shown to be applicable to the auroral field line plasma if a current is present along the magnetic field. A simulation performed with the equations shows that the field line plasma exhibits dynamic behavior after a field-aligned current appears. The resulting momentum gain by O(+) ions can be sufficient for causing a potential drop of several kilovolts along the field line.

Mitchell, H. G., Jr.↗

LaRC computational dynamics overview

Present research centers on the development of advanced computational methods for transient simulation analyses. Aircraft, launch vehicles and space structure components are potential applications, but primary focus is presently on large space structures. There are both in-house and out-of-house activities. The in-house activity centers around the development of a multibody simulation tool for truss-like structures called LATDYN for Large Angle Transient DYNamics. Multibody analysis involves articulation of structural components as well as robotic maneuvers. These items are necessary for construction (erection or deployment) of large space structures in orbit and the carrying out of certain operations on board the space station. Thus, part of the in-house activity involves the development of methods which treat the changing mass, stiffness and constraints associated with articulating systems. The out-of-house activity involves subcycling, development of large deformation/motion beam formulation, constraint stabilization and direct time integration transient algorithms in parallel computing.

Husner, J. M.↗

Electron acceleration by Alfven waves in the magnetosphere

The self-consistent electron kinetics of Alfven waves on the electron inertial scale is studied using a 2D hybrid-kinetic description. The ions follow a fluid description for Alfven waves at frequencies below the ion cyclotron frequency. The parallel electron dynamics is treated kinetically using particle-in-cell techniques. In this model, the electron plasma mode is eliminated, and only the physics of the Alfven waves is retained. At sufficiently large amplitudes, it is found that oblique Alfven waves break due to finite electron inertia in a cold plasma. The consequence of wave breaking is the formation of an electron beam which can be unstable to the beam-plasma instability. The electrons supporting the parallel current thermalize into a non-Maxwellian distribution with an energetic tail up to several keV, assuming a reasonable magnetospheric Alfven speed. In hot plasma simulations, electron trapping is the principal mechanism of electron acceleration. It is proposed that wave breaking or electron trapping of oblique Alfven waves at 1 R(E) can result in electron acceleration and may explain some observed auroral phenomena.

Hui, C.-H.↗

Real-time dynamics and control strategies for space operations of flexible structures

This project (NAG9-574) was meant to be a three-year research project. However, due to NASA's reorganizations during 1992, the project was funded only for one year. Accordingly, every effort was made to make the present final report as if the project was meant to be for one-year duration. Originally, during the first year we were planning to accomplish the following: we were to start with a three dimensional flexible manipulator beam with articulated joints and with a linear control-based controller applied at the joints; using this simple example, we were to design the software systems requirements for real-time processing, introduce the streamlining of various computational algorithms, perform the necessary reorganization of the partitioned simulation procedures, and assess the potential speed-up realization of the solution process by parallel computations. The three reports included as part of the final report address: the streamlining of various computational algorithms; the necessary reorganization of the partitioned simulation procedures, in particular the observer models; and an initial attempt of reconfiguring the flexible space structures.

Park, K. C.↗

Toward fully automated UED operation using two-stage machine learning model

To demonstrate the feasibility of automating UED operation and diagnosing the machine performance in real time, a two-stage machine learning (ML) model based on self-consistent start-to-end simulations has been implemented. This model will not only provide the machine parameters with adequate precision, toward the full automation of the UED instrument, but also make real-time electron beam information available as single-shot nondestructive diagnostics. Furthermore, based on a deep understanding of the root connection between the electron beam properties and the features of Bragg-diffraction patterns, we have applied the hidden symmetry as model constraints, successfully improving the accuracy of energy spread prediction by a factor of five and making the beam divergence prediction two times faster. The capability enabled by the global optimization via ML provides us with better opportunities for discoveries using near-parallel, bright, and ultrafast electron beams for single-shot imaging. It also enables directly visualizing the dynamics of defects and nanostructured materials, which is impossible using present electron-beam technologies.

36 MATERIALS SCIENCE↗

Examination of stiff ion temperature gradient mode physics in simulations of DIII-D H-mode transport

A systematic evaluation of gyrokinetic and gyrofluid model predictions of ion temperature gradient (ITG) stability and transport using parameters from DIII-D high confinement mode (H-mode) plasmas has been performed. The nonlinear CGYRO code is used to make the gyrokinetic predictions, and the quasilinear TGLF model for the corresponding gyrofluid predictions. The assessments are made at three radii (normalized toroidal flux ρtor = 0.4, 0.55, and 0.7) in three different plasma scenarios with varying levels of neutral beam heating and torque. For each of the nine cases (3 radii × 3 scenarios) considered, ITG turbulence is found to be the dominant long-wavelength instability and transport mechanism. The inclusions of both transverse magnetic fluctuations and dynamic fast beam ions are stabilizing for all cases considered, with strongest effects seen at ρor = 0.4 where the fast ion population and normalized plasma pressure β = 2μ0nT/B2 are highest. The further inclusion of parallel magnetic fluctuations does not have a meaningful impact on the ITG turbulence in these scenarios, but does destabilize (in combination with fast ions) new high-frequency instabilities at ρtor = 0.4 in the high power scenarios. In each case the linear and nonlinear ITG critical gradients are predicted to be lower than the measured ITG scale lengths and their associated uncertainties. Inclusion of equilibrium flow shear in the transport predictions generally leads to an upshift in effective critical gradient rather than a qualitative change in the predicted stiffness, with stronger responses typically seen in the gyrokinetic predictions than in the gyrofluid results. However, in most cases these upshifted gradients still remain below the measured values and their uncertainties. Although the predicted critical gradients are below the measured gradients, both models predicted flux-matching gradients consistent with measured values in six of the nine cases considered, with no clear systematic over- or underprediction. Thus, while the experimental ion temperature profiles do not appear to be closely pinned to the ITG critical gradient, both gyrokinetic and gyrofluid models are able to accurately match the measured gradients reasonably well in most cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Detector and Beamline Simulation for Next-Generation High Energy Physics Experiments

The success of high energy physics programs relies heavily on accurate detector simulations and beam interaction modeling. The increasingly complex detector geometries and beam dynamics require sophisticated techniques in order to meet the demands of current and future experiments. Common software tools used today are unable to fully utilize modern computational resources, while data-recording rates are often orders of magnitude larger than what can be produced via simulation. In this paper, we describe the state, current and future needs of high energy physics detector and beamline simulations and related challenges, and we propose a number of possible ways to address them.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Numerical Simulation of Flow Field Within Parallel Plate Plastometer

Parallel Plate Plastometer (PPP) is a device commonly used for measuring the viscosity of high polymers at low rates of shear in the range 10 4 to 10 9 poises. This device is being validated for use in measuring the viscosity of liquid glasses at high temperatures having similar ranges for the viscosity values. PPP instrument consists of two similar parallel plates, both in the range of 1 inch in diameter with the upper plate being movable while the lower one is kept stationary. Load is applied to the upper plate by means of a beam connected to shaft attached to the upper plate. The viscosity of the fluid is deduced from measuring the variation of the plate separation, h, as a function of time when a specified fixed load is applied on the beam. Operating plate speeds measured with the PPP is usually in the range of 10 -3 cm/s or lower. The flow field within the PPP can be simulated using the equations of motion of fluid flow for this configuration. With flow speeds in the range quoted above the flow field between the two plates is certainly incompressible and laminar. Such flows can be easily simulated using numerical modeling with computational fluid dynamics (CFD) codes. We present below the mathematical model used to simulate this flow field and also the solutions obtained for the flow using a commercially available finite element CFD code.

Basil N. Antar↗

Dynamic Analysis of a Six-Cable Parallel Robot for Automated Panelized Building Retrofits

Cable-Driven Parallel Robots (CDPRs) are highly suitable for automated panelized building retrofits, thanks to their compact footprint and high payload-to-weight ratio. A common CDPR configuration featuring eight cables, where the anchors form a rectangular prism in front of the building facade, offers a large wrench feasible workspace and good control versatility. However, installing upper anchor points requires additional support structures, such as towers or beams, increasing setup complexity and posing logistical challenges in construction settings. To mitigate these challenges, we propose a six-cable CDPR model specifically designed for automated panelized building retrofits. Although the feasible workspace is limited, our analysis shows that the proposed CDPR adequately covers the critical areas required for panel installation. To validate that the six-cable system can effectively transport the end effector to the desired installation pose, we calculated optimal trajectories based on a constrained dynamic model. The simulation results of the six-cable CDPR demonstrate promising potential for automated panelized building retrofits, effectively balancing simplicity, cost-effectiveness, and functionality.

Liu, Yifang [ORNL] (ORCID:0000000190817417)↗

Heating and scattering of ring-beam distributions by turbulence

Pickup ions in the solar wind are initially are born in ring-beam distributions, i.e. f(v) varies as delta(v(sub perpendicular) - V(sub sw)sin(Theta)) delta(v(sub parallel) - V(sub sw)cos(Theta)), where Theta is the angle between the solar wind velocity and the IMF(Interplanetary Magnetic Field), and V(sub sw) is the solar wind speed. Often the distribution has been presumed to relax to a distribution that is isotropic in Theta and essentially mono-energetic, a shell or a 'bi-spherical distribution.' However solar wind turbulence is capable of heating the ring distribution on the timescale of a few tens of gyroperiods, a timescale not greatly distinct from that required for pitch angle scattering to a shell. To describe this effect, we have performed test-particle studies of the heating/scattering of the ring beam distribution by MHD turbulence, adopting various models for the MHD fluctuations, including slab and fully dynamic 2D and 3D incompressible turbulence. Furthermore, a system composed of a cold ion ring and a background plasma is unstable to several kinetic plasma instabilities. We carried out kinetic simulations of the ring beam distribution, showing that plasma instabilities also rapidly energize and scatter particles. Results will be presented comparing relaxation and heating rates of the ring-beam distribution by the various mechanisms.

Gray, P. C.↗

On Lidar Sounding of the Atmosphere to Estimate Static and Dynamic Characteristics of Aerosol Inhomogeneities

1. A possible application of intensity fluctuations of a pulse light signal reflected by atmospheric aerosols is analyzed by the correlation method to evaluate static (medium sizes, shape) and dynamic (speed and direction of movement, lifetime) characteristics of aerosol inhomogeneities. The aerosol inhomogeneities are assumed to be expanded, pressed, disintegrated and originated constantly in accordance with random laws, the set of inhomogeneities as a whole traveling together with air masses and having predominant movement in wind direction. It is shown that the characteristics of aerosol inhomogeneities considered can be expressed by the coefficients of the correlation function expansion of the reflected signal fluctuation intensity in Tailor series. 2/ Correlation systems for evaluating static and dynamic characteristics of driving objects can be divided into two types according to the kind and quantity of used information: the systems with coordinates of the information removal "points" to be fixed in space, and the systems with a parallel simultaneous information removal at discrete moments of time. The systems for determination of wind direction considered in are the examples of the first type system. However, the operating information removal for two points is insufficient to estimate completely static and dynamic characteristics of inhomogenities, their quantity ought to be increased up to three of them for two-dimensional problem and up to four of them for three-dimensional problem as it is usually done in the ionospheric studies. The second type systems are used for the investigation of a medium shape and speed of the clouds according to photographs made from satellites. These systems are also used for solution of navigation problems. The use of optical quantum generators with a scanning beam is seen to increase greatly the working information removal in comparison with the first type systems. Nevertheless, scanning rate is not sufficient sometimes in order· to consider a general picture of aerosol inhomogeneities to be stationary. In this connection the use of the systems of second type treatment becomes a matter of essential difficulty. 3. Aerosol inhomogeneities simulation has been carried out on the basis of the digital computer experiments with the aim of estimating static and dynamic characteristics of inhomogeneities by an optical beam in the atmosphere at different scanning procedures. The dependence of determination accuracy of these characteristics on the type of chosen laws of aerosol particle distribution in the atmosphere, the parameters of inhomogeneities geometry, their speed and the law of scanning have been obtained.

Zuev, V. E.↗

Distributed Stochastic Optimization of a Neural Representation Network for Time-Space Tomography Reconstruction

4D time-space reconstruction of dynamic events or deforming objects using X-ray computed tomography (CT) is an important inverse problem in non-destructive evaluation. Conventional back-projection based reconstruction methods assume that the object remains static for the duration of several tens or hundreds of X-ray projection measurement images (reconstruction of consecutive limited-angle CT scans). However, this is an unrealistic assumption for many in-situ experiments that causes spurious artifacts and inaccurate morphological reconstructions of the object. To solve this problem, we propose to perform a 4D time-space reconstruction using a distributed implicit neural representation (DINR) network that is trained using a novel distributed stochastic training algorithm. Our DINR network learns to reconstruct the object at its output by iterative optimization of its network parameters such that the measured projection images best match the output of the CT forward measurement model. Here, we use a forward measurement model that is a function of the DINR outputs at a sparsely sampled set of continuous valued 4D object coordinates. Unlike previous neural representation architectures that forward and back propagate through dense voxel grids that sample the object's entire time-space coordinates, we only propagate through the DINR at a small subset of object coordinates in each iteration resulting in an order-of-magnitude reduction in memory and compute for training. DINR leverages distributed computation across several compute nodes and GPUs to produce high-fidelity 4D time-space reconstructions. We use both simulated parallel-beam and experimental cone-beam X-ray CT datasets to demonstrate the superior performance of our approach.

36 MATERIALS SCIENCE↗

Status and future plans for C 3 R&D

C 3 is an opportunity to realize an e + e - collider for the study of the Higgs boson at √s = 250 GeV, with a well defined upgrade path to 550 GeV while staying on the same short facility footprint. C 3 is based on a fundamentally new approach to normal conducting linear accelerators that achieves both high gradient and high efficiency at relatively low cost. Given the advanced state of linear collider designs, the key system that requires technical maturation for C 3 is the main linac. This paper presents the staged approach towards a facility to demonstrate C 3 technology with both Direct (source and main linac) and Parallel (beam delivery, damping ring, ancillary component) R&D. The primary goal of the C 3 Demonstration R&D Plan is to reduce technical and cost risk by building and operating the key components of C 3 at an adequate scale. This R&D plan starts with the engineering design, and demonstration of one cryomodule and will culminate in the construction of a 3 cryomodule linac with pre-production prototypes. This R&D program would also demonstrate the linac rf fundamentals including achievable gradient and gradient stability over a full electron bunch train and breakdown rates. It will also investigate beam dynamics including energy spread, wakefields, and emittance growth. This work will be critical to confirm the suitability of the C 3 beam parameters for the physics reach and detector performance in preparation for a Conceptual Design Report (CDR), as well as for follow-on technology development and industrialization. The C 3 Demonstration R&D Plan will open up significant new scientific and technical opportunities based on development of high-gradient and high-efficiency accelerator technology. It will push this technology to operate both at the GeV scale and mature the technology to be reliable and provide high-brightness electron beams. The timeline for progressing with C 3 technology development will be governed by practical limitations on both the technical progress and resource availability. It consists of four stages: Stage 0) Ongoing fundamental R&D on structure prototypes, damping and vibrations. Stage 1) Advancing the engineering maturity of the design and developing start-to-end simulations including space-charge and wakefield effects. This stage will include testing of strucutres operating at cryogenic temperatures. Beam tests would be performed with high beam current to test full beam loading. Stage 2) Production and testing of the first cryomodule at cryogenic temperatures. This would provide sufficient experimental data to compile a CDR and it is anticipated for Stage 2 to last 3 years and to culminate with the transport of photo-electrons through the first cryomodule. Stage 3) Updates to the engineering design of the cryomodules, production of the second and third cryomodule and their installation. Lower charge and lower emittance beams will be used to investigate emittance growth. The successful full demonstration of the 3 cryomodules to deliver up to a 3 GeV beam and achieve the C 3 five gradient will allow a comprehensive and robust evaluation of the technical design of C 3 as well as mitigate technical, schedule, and cost risks required to proceed with a Technical Design Report (TDR).

radiation hardened magnets↗