Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulation analytics”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Generation of polarization entanglement via the quantum Zeno effect

The quantum Zeno effect reveals that continuous observation of a quantum system can significantly alter its evolution. Here, we present a method for establishing polarization entanglement between two initially unentangled photons in coupled waveguides via the quantum Zeno effect. We support our analytical investigation with numerical simulations of the underlying Schrodinger equation describing the system. Further, we extend our technique to three coupled waveguides in a planar configuration and determine the parameters required to generate three-qubit W-states. In contrast to existing schemes based on a vacuum and single-photon encoding, the polarization encoding in our approach is compatible with quantum information protocols that remove photon loss through post-selection. Our findings offer a powerful quantum state engineering approach for photonic quantum information technologies.

Nodurft, Ian C.↗
DOE OSTI2023

Refractive axicon for X-ray microscopy applications: design, optimization, and experiment

In a full-field transmission X-ray microscopy (TXM) setup, a condenser X-ray optical element is used to illuminate the sample by condensing the X-ray beam delivered by the synchrotron storage ring. On-going and future upgrades of synchrotron facilities to diffraction-limited storage rings will pose new challenges to these TXM setups, such as much smaller X-ray beams on the condenser. Here, we demonstrate that a refractive axicon can be used as an X-ray beam shaper to match the ring-shaped aperture of the condenser. Aiming at more efficient use of the incoming X-ray intensity, we explore several axicon designs both analytically and with numerical simulations. The axicons were produced by two-photon polymerization 3D printing on thin silicon nitride membrane substrates. The first characterization of the axicon was carried out at the TOMCAT beamline of the Swiss Light Source (Switzerland).

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Improvements to the New CTH Code Verification & Validation Test Suite (FY2020)

The CTH multiphysics hydrocode, which is used for a wide range of important calculations, has undertaken in recent years to overhaul its software quality and testing processes. A key part of this effort entailed building a new, robust V&V test suite made up of traditional hydrocode verification problems, such as those listed in the ASC Tri-Lab Test Suite and the Enhanced Tri-Lab Test Suite, as well as validation problems for some of CTHs most frequently used equations of state, materials models, and other key capabilities. Substantial progress towards this goal was made in FY19. In FY20, this test suite has been expanded to include verification and validation tests of the Sesame and JWL equation of state models as well as the Mader verification problem from the Tri-Lab Test Suite and the Blake verification problem - a linear elastic analog to the Hunter problem from the Enhanced Tri-Lab Test Suite. This report documents CTH performance on the new test suite problems. Verification test results are compared to analytic solutions and, for most tests, convergence results are presented. Validation test results are compared to experimental data and mesh refinement studies are included. CTH performs well overall on the new test problems. Convergence rates for the Blake and Mader problems are comparable to those for similar ASC codes. The JWL and Sesame verification tests show good agreement with analytic solutions. Likewise, CTH simulation results show good agreement with experimental validation data for the Sesame and JWL equations of state for the materials tested. Future V&V work will focus on adding tests for other key capabilities like fracture and high explosive models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Beam Stability Overview [Slides]

Radiographic spot size is affected by many factors, including beam motion due to instabilities. Instabilities and effects that can influence the radiation-source spot size are being evaluated for the Scorpius flash-radiography accelerator. The Scorpius LIA design has a direct influence on beam stability. Many features of the design were specifically incorporated to mitigate instability. The major beam dynamics concerns for high-current LIAs are being evaluated for the long, 102-cell Scorpius design. For these results, all tools available were used including experimental data, analytic theory, and physics simulation codes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

5G Energy FRAME: The Design and Implementation of Data, Model, and Use Case (Year 2 Report)

This report summarizes the Year 2 work of Pacific Northwest National Laboratory’s (PNNL’s) 5G Fabricated Resource and Asset Management Encompassment for energy infrastructure (Energy FRAME) project funded by the Department of Energy Office of Science’s Advanced Scientific Computing Research (ASCR) Program. In this report, newest 5G equipment testing results are presented, along with two 5G-enabled AI/ML examples for grid applications; in addition, the work flow of grid edge, cloud, and High Performance Computing (HPC) platform is introduced, to support and interface the cross-domain simulation for power system transmission, distribution, and communication networks. Last but not least, the outlook for Year 3 work and the overarching impact of 5G Energy FRAME work to a multitude of stakeholders are provided. Additional 5G performance data now is shared through the publicly available weblink, https://www.pnnl.gov/projects/5g-energy-frame/publications

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Radial drift of plasma blobs in a toroidal magnetic field with fully kinetic and reduced fluid models

In curved magnetic geometries, field-aligned regions of enhanced plasma pressure and density, termed ‘blobs,’ move as coherent filaments across the magnetic field lines. Coherent blobs account for a significant fraction of transport at the edges of magnetic fusion experiments and arise in naturally-occurring space plasmas. This work examines the dynamics of blobs with a fully kinetic electromagnetic particle-in-cell code and with a drift-reduced fluid code. In low-beta regimes with moderate blob speeds, good agreement is found in the maximum blob velocity between the two simulation schemes and simple analytical estimates. The fully kinetic code demonstrates that blob speeds saturate near the initial sound speed, which is a regime outside the validity of the reduced fluid model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

NIOLD Run 4 Report

In this report we discuss the motivation, methods, and analysis of the IOTA Run 4 studies of the Nonlinear Integrable Optics, Landau Damping (NIOLD) experiment. We introduce Landau Damping, an effect that damps collective instabilities in particle accelerators. We also introduce and discuss a new method to measure stability diagrams which quantify the strength of Landau Damping, that employs an antidamper. We also present the data collection process, the data analysis procedure, and the preliminary results. The first results qualitatively agree with the analytical predictions and the simulations. Next steps for the current data and goals for future data are also discussed.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Fast and Thermal Neutron Spectrum Dosimetry Measurements in the Advanced Test Reactor large-B and small-I Positions Following the Sixth Core Internals Change-out [Slides]

The Advanced Test Reactor (ATR) has a wide variety of irradiation positions that have had experiments that were developed by users from around the world. Most experiment irradiations rely on thoroughly benchmarked numerical models. However, some irradiation positions in ATR are not as well- characterized and are complicated by spectral perturbations from control cylinder orientation. The “small I” and “large-B” irradiation positions are located nearby control cylinders and suffer from these flux perturbations from control cylinder orientations that change during an irradiation cycle to maintain the desired core power distribution. Models predicted that these types of position would exhibit both spectral shifts and amplitude changes in neutron flux, but few measurements have been conducted to benchmark these predictions. Recently, requalification testing was performed to confirm the operational readiness of the ATR following the completion of the Core Internals Change-out (CIC). These tests provided a unique opportunity to validate the analytical methods used to simulate the ATR because nearly all components in the reactor were in a clean as-built state, significantly reducing modelling uncertainties. One subset of the testing included characterization of the fast and thermal neutron flux in the “small-I” and “large-B” positions using silver, cobalt, and nickel neutron dosimetry. In contrast to typical irradiation cycles, the control cylinders were held in position during the post-CIC nuclear testing. The specific activity of these dosimeter wires was measured following two low-power tests with different control cylinder positions.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2025

Positron Beam Loading and Acceleration in the Blowout Regime of a Plasma Wakefield Accelerator

Plasma wakefield acceleration in the nonlinear blowout regime has achieved marked milestones in electron beam acceleration, demonstrating high acceleration gradients and energy efficiency while preserving excellent beam quality. However, this regime is deemed unsuitable for achieving positron acceleration of comparable results, which is vital for future compact electron–positron colliders. In this article, we find that an intense positron beam loaded at the back of beam-driven blowout cavity can self-consistently induce the focusing field and flatten the longitudinal wakefield, leading to stable, high-efficiency, and high-quality positron acceleration. This is achieved through the formation of an on-axis electron filament induced by positron beam load, which shapes the plasma wakefield in a distinct way compared to electron beam load in the blowout regime. Via a nonlinear analytic model and numerical simulations, we explain the novel beam loading effects of the interaction between the on-axis filament and the blowout cavity. High-fidelity simulations show that a high-charge positron beam can be accelerated with >20% energy transfer efficiency, ~1% energy spread, and ~1 mm·mrad normalized emittance, while considerably depleting the energy of the drive beam. The concept can also be extended to simultaneous acceleration of electron and positron beams and high transformer ratio positron acceleration as well. This development offers a new route for the application of plasma wakefield acceleration into particle physics.

Science & Technology - Other Topics↗
DOE OSTI2020

Milky Way Satellite Census. I. The Observational Selection Function for Milky Way Satellites in DES Y3 and Pan-STARRS DR1

We report the results of a systematic search for ultra-faint Milky Way satellite galaxies using data from the Dark Energy Survey (DES) and Pan-STARRS1 (PS1). Together, DES and PS1 provide multi-band photometry in optical/near-infrared wavelengths over ~80% of the sky. Our search for satellite galaxies targets ~25,000 deg$^2$ of the high-Galactic-latitude sky reaching a 10$\sigma$ point-source depth of $\gtrsim$ 22.5 mag in the $g$ and $r$ bands. While satellite galaxy searches have been performed independently on DES and PS1 before, this is the first time that a self-consistent search is performed across both data sets. We do not detect any new high-significance satellite galaxy candidates, while recovering the majority of satellites previously detected in surveys of comparable depth. We characterize the sensitivity of our search using a large set of simulated satellites injected into the survey data. We use these simulations to derive both analytic and machine-learning models that accurately predict the detectability of Milky Way satellites as a function of their distance, size, luminosity, and location on the sky. To demonstrate the utility of this observational selection function, we calculate the luminosity function of Milky Way satellite galaxies, assuming that the known population of satellite galaxies is representative of the underlying distribution. We provide access to our observational selection function to facilitate comparisons with cosmological models of galaxy formation and evolution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Full Transport General Relativistic Radiation Magnetohydrodynamics for Nucleosynthesis in Collapsars

We model a compact black hole-accretion disk system in the collapsar scenario with full transport, frequency dependent, general relativistic radiation magnetohydrodynamics. We examine whether or not winds from a collapsar disk can undergo rapid neutron capture (r-process) nucleosynthesis and significantly contribute to solar r-process abundances. We find the inclusion of accurate transport has significant effects on outflows, raising the electron fraction above ${Y}_{{\rm{e}}}\sim 0.3$ and preventing third-peak r-process material from being synthesized. We analyze the time evolution of neutrino processes and electron fraction in the disk and present a simple one-dimensional model for the vertical structure that emerges. Finally, we compare our simulation to semi-analytic expectations and argue that accurate neutrino transport and realistic initial and boundary conditions are required to capture the dynamics and nucleosynthetic outcome of a collapsar.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

xSDK: Building an ecosystem of highly efficient math libraries for exascale

Current efforts to build increasingly powerful computer architectures are opening up new avenues for more complex and higher fidelity simulations coupled with data analytics and learning, leading to new scientific insights and deeper understanding. At one extreme, exascale computers will be much faster than previous computer generations (performing 10 18 operations per second—that is, 1,000 times faster than petascale). To achieve these performance improvements, computer architectures are becoming increasingly complex, with deep memory hierarchies, very high node and core counts, and heterogeneous features such as graphics processing units (GPUs). Such architectural changes impact the full breadth of computing scales, as heterogeneity pervades even current-generation laptops, workstations, and moderate-sized clusters. While emerging advanced architectures provide unprecedented opportunities, they also present significant challenges for developers of scientific applications, such as multiphysics and multiscale codes, who must adapt their software to handle disruptive changes in architectures and new programming models that have not yet stabilized. Developers must consider increasing concurrency while reducing communication and synchronization, and other complexities such as the potential for using mixed precision to leverage the compute power available in low-precision tensor cores. On one hand, developers must implement new scientific capabilities, which in turn increase code complexity. On the other hand, the codes must be ported to new architectures, requiring the inclusion of new programming models and the restructuring of code to achieve good performance. Addressing these issues is beyond the capability of any single person or team—leading to the need for collaboration among many teams, who encapsulate their expertise in reusable software and work together to create sustainable software ecosystems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Determining Compression Characteristics of Honeycomb Material - 19674

The objective is to determine compression test characteristics of stainless steel honeycomb material to be able to represent honeycomb structures accurately in analytical models used to simulate hypothetical accident scenarios of shipping packages. Honeycomb is a material used primary in the aerospace industry due to its high strength-to-weight ratio. Because it doesn't have a shelf life and can withstand high heats, it is an excellent candidate for a structural material in package designs. Honeycomb has orthotropic material properties. The honeycomb currently being evaluated is made of metal ribbons spot welded together to form hexagon pattern between metal plates. The hexagon ribbons are brazed to the metal plates. The complexity of the honeycomb significantly lengthens the simulation time, which is further complicated when the brazing, welding, and imperfections of the material is considered. This means to affectively represent the honeycomb in simulation programs, such as Abaqus, the material needs to be approximated as a uniform orthotropic material. This requires material properties in each of the three directions, T, W, and L. The T direction is defined as the 'strong' direction, perpendicular to the honeycomb sheet. The L direction is parallel the ribbon and the W direction is perpendicular to the ribbon. Honeycomb has 3 different compression stages. Stage 1 is the initial compression where the honeycomb maintains its structural integrity and does not permanently deform from loads from normal operations. Stage 2 is when the initial buckling causes irreversible damage to the honeycomb. Stage 2 is the one we are most interested in because it absorbs the most energy from a hypothetical accident scenario. During Stage 2, the honeycomb fails layer-by-layer, indicating that the more layers, the longer material crushing is sustained, and thus the more energy that is absorbed. Stage 3 is final compression, similar to compressing solid metal. During stage 3 the stress increases with a diminishing rate of elongation and the honeycomb is completely failed where the plates between the honeycomb core sandwich the crushed honeycomb ribbon. The stress-strain graph bellow shows compression test of 3 different samples in the T direction. The red lines divide the three stages. The far left is Stage 1, the middle stage 2, and the right stage 3. The data collected is displayed in the table below. The expressions in the table represent the regression formula that represents each stage of the compression on a stress-strain graph. Stage 1 is assumed to intersect with the origin; data for stage 3 for the crush in W and L directions where unable to be gathered due to the nature of failure for those directions. Stage 1 and 2 are linear regressions while stage 3 is a degree 2, to best match the curve. Linear regression is used on stage 2 because ultimately that represents the energy absorbed. With more data, a degree n x 2 regression would be a more appropriate regression for stage 2, where n is the number of layers. The table below is the approximation of each scenario and stage. Each stage starts and ends at the intersection of the next stage. There where two main objectives for these tests: firstly to understand how honeycomb performs under extreme compression, and secondly to be able to numerically represent the honeycomb structure. Both of these where completed to varying degrees. It is important understood how the honeycomb fails. If it is crushed in the T direction it retains its integrity even after being crushed, however when crushed in the L and W crush direction, if it fails, it disintegrates, and loses all integrity. Also, the more layers the more time the material spends in stage 2. Failure is started by buckling, thus if there is any imperfection, the stress will not spike but transition straight into stage 2. The data collected and aggregated can be used for initial simulation of honeycomb used in packages. The initial testing has set the ground work for more data to be collected in order to verify results and to allow more confidence in the simulation results. This is only the initial data collected. The next steps is to continue to collect more data to verify results and to test more variants of honeycomb, with different brazing, layers, and shape.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A317_Fast and Thermal Neutron Spectrum Dosimetry Measurements in the Advanced Test Reactor large-B and small-I Positions Following the Sixth Core Internals Change-out

The Advanced Test Reactor (ATR) has a wide variety of irradiation positions that have had experiments that were developed by users from around the world. Most experiment irradiations rely on thoroughly benchmarked numerical models. However, some irradiation positions in ATR are not as well- characterized and are complicated by spectral perturbations from control cylinder orientation. The “small I” and “large-B” irradiation positions are located nearby control cylinders and suffer from these flux perturbations from control cylinder orientations that change during an irradiation cycle to maintain the desired core power distribution. Models predicted that these types of position would exhibit both spectral shifts and amplitude changes in neutron flux, but few measurements have been conducted to benchmark these predictions. Recently, requalification testing was performed to confirm the operational readiness of the ATR following the completion of the Core Internals Change-out (CIC). These tests provided a unique opportunity to validate the analytical methods used to simulate the ATR because nearly all components in the reactor were in a clean as-built state, significantly reducing modelling uncertainties. One subset of the testing included characterization of the fast and thermal neutron flux in the “small-I” and “large-B” positions using silver, cobalt, and nickel neutron dosimetry. In contrast to typical irradiation cycles, the control cylinders were held in position during the post-CIC nuclear testing. The specific activity of these dosimeter wires was measured following two low-power tests with different control cylinder positions.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2026

Analytical small-signal stability analysis of low-inertia power system frequency response considering secondary frequency regulation

Modern power systems are increasingly vulnerable to frequency instability as inverter-based resources (IBRs) replace synchronous machines and reduce system rotational inertia. Existing small-signal frequency stability assessment methods are either computationally intensive, relying on simulation-driven approaches, or lack analytical stability regions that explicitly account for secondary frequency response (SFR). This paper introduces new analytical frameworks that enable evaluation small-signal frequency stability while explicitly incorporating tunable IBR and SFR parameters. Using Kharitonov’s theorem with an overbounding approach, explicit small-signal stability criteria are derived. In addition, based on Białas’ criterion and Hurwitz matrix, analytical stability regions are established to reveal feasible design spaces for SFR and IBR parameters tuning. Extensive Matlab/Simulink-based simulations validate the accuracy and computational efficiency of the proposed methods, demonstrating that coordinated tuning of SFR and IBR parameters can substantially enhance system resilience. By bridging analytical rigor with practical tunability, this work provides an analytical framework for assessing small-signal frequency stability in low-inertia grids, supporting the real-time, scalable, and resilient operation of sustainable power systems.

14 SOLAR ENERGY↗
DOE OSTI2021

Hybrid simulations of sub-cyclotron compressional and global Alfven Eigenmode stability in spherical tokamaks

A comprehensive numerical study has been conducted in order to investigate the stability of beam-driven, sub-cyclotron frequency compressional (CAE) and global (GAE) Alfven Eigenmodes in low aspect ratio plasmas for a wide range of beam parameters. The presence of CAEs and GAEs has previously been linked to anomalous electron temperature profile flattening at high beam power in NSTX experiments, prompting further examination of the conditions for their excitation. Linear simulations are performed with the hybrid MHD-kinetic initial value code HYM in order to capture the general Doppler-shifted cyclotron resonance that drives the modes. Three distinct types of modes are found in simulations -- co-CAEs, cntr-GAEs, and co-GAEs -- with differing spectral and stability properties. The simulations reveal that unstable GAEs are more ubiquitous than unstable CAEs, consistent with experimental observations, as they are excited at lower beam energies and generally have larger growth rates. Local analytic theory is used to explain key features of the simulation results, including the preferential excitation of different modes based on beam injection geometry and the growth rate dependence on the beam injection velocity, critical velocity, and degree of velocity space anisotropy. The background damping rate is inferred from simulations and estimated analytically for relevant sources not present in the simulation model, indicating that co-CAEs are closer to marginal stability than modes driven by the cyclotron resonances.

Compressional Alfven Eigenmodes↗
DOE OSTI2020

Analyzing the Validity of Brazilian Testing Using Digital Image Correlation and Numerical Simulation Techniques

Characterizing the mechanical behavior of rocks plays a crucial role to optimize the fracturing process in unconventional reservoirs. However, due to the intrinsic anisotropy and heterogeneity in unconventional resources, fracture process prediction remains the most significant challenge for sustainable and economic hydrocarbon production. During the deformation tracking under compression, deploying conventional methods (strain gauge, extensometer, etc.) is insufficient to measure the deformation since the physical attachment of the device is restricted to the size of the sample, monitoring limited point-wise deformation, producing difficulties in data retrieval, and a tendency to lose track in failure points, etc. Where conventional methods are limited, the application of digital image correlation (DIC) provides detailed and additional information of strain evolution and fracture patterns under loading. DIC is an image-based optical method that records an object with a camera and monitors the random contrast speckle pattern painted on the facing surface of the specimen. To overcome the existing limitations, this paper presents numerical modeling of Brazilian disc tests under quasi-static conditions to understand the full-field deformation behaviors and finally, it is validated by DIC. As the direct tensile test has limitations in sample preparation and test execution, the Brazilian testing principle is commonly used to evaluate indirectly the tensile strength of rocks. The two-dimensional numerical model was built to predict the stress distribution and full-field deformation on Brazilian disc under compression based on the assumptions of a homogenous, isotropic and linear elastic material. The uniaxial compression test was conducted using the DIC technique to determine the elastic properties of Spider Berea sandstone, which were used as inputs for the simulation model. The model was verified by the analytical solution and compared with the digital image correlation. The numerical simulation results showed that the solutions matched reasonably with the analytical solutions where the maximum deviation of stress distribution was obtained as 14.59%. The strain evolution (normal and shear strains) and displacements along the central horizontal and vertical planes were investigated in three distinguishable percentages of peak loads (20%, 40%, and 90%) to understand the deformation behaviors in rock. The simulation results demonstrated that the strain evolution contours consistently matched with DIC generated contours with a reasonable agreement. The changes in displacement along the central horizontal and vertical planes showed that numerical simulation and DIC generated experimental results were repeatable and matched closely. In terms of validation, Brazilian testing to measure the indirect tensile strength of rocks is still an issue of debate. The numerical model of fracture propagation supported by digital image correlation from this study can be used to explain the fracturing process in the homogeneous material and can be extended to non-homogeneous cases by incorporating heterogeneity, which is essential for rock mechanics field applications.

58 GEOSCIENCES↗
DOE OSTI2026

Multibody for Everybody (M4E): A Symbolic Dynamics Modeling Tool with Applications in Simulation, Control, and Optimization

Developing the analytical model of a multibody system is often the initial step in control and optimization. The analytical model (equations of motion) describes a system’s time evolution under specified forcing conditions. Although developing these equations is easy for simple systems, this process becomes more complex for systems composed of multiple bodies. Deriving equations of motion for complex multibody systems requires specialized expertise in multibody dynamics, is time-consuming, and is susceptible to error. To address this issue, this paper presents an open-source, easy-to-use, systematic framework to derive symbolic equations of motion in both Python and MATLAB using the joint coordinate formulation. This formulation results in a set of ordinary differential equations that use the minimum set of coordinates needed to model a system. The symbolic representation provides better insight into the influence of design parameters on system performance, facilitates sensitivity analysis and parameter studies, and supports direct implementation of control and optimization routines. The tool enables numerical simulation for specified parameter sets, is modular for straightforward integration with other tools and libraries, and allows incorporation of hydrodynamics, mooring, and other external forces. The result is a reproducible, extensible pipeline for modeling, simulation, and design of complex multibody systems. The proposed tool is versatile and can be applied to domains such as robotics, control, and design. In addition, we integrated external libraries that provide capabilities for modeling offshore systems such as underwater robots and marine energy converters.

16 TIDAL AND WAVE POWER↗