Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation optimization”

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 1,027 records · Page 57

Simulating large one-dimensional neutral-atom quantum systems

While abstract models of quantum computation assume a closed system of two-level states, practical quantum devices inevitably couple to the environment in some way, creating sources of noise. Understanding the tolerance to noise of specific quantum algorithms run on specific devices is important for determining the feasibility of quantum computing in the current noisy intermediate-scale quantum era. Of particular interest is understanding the noise sensitivity of these devices as more qubits are added to the system. Classical simulations are a useful tool to understand the effects of this noise, but direct classical simulations of open quantum systems are burdened by an exponentially growing cost in the number of qubits and a large local Hilbert space dimension. For onedimensional, shallow circuits, using tensor networks can replace this exponential cost with a linear one and simulate far wider systems than what would normally be available. In this paper, we describe a tensor network simulation of a neutral atom quantum system under the presence of noise, while introducing a purity-preserving truncation technique that compromises between the simplicity of the matrix product state and the positivity of the matrix product density operator. We apply this simulation to a near-optimized iteration of the quantum approximate optimization algorithm on a transverse field Ising model in order to investigate the influence of large system sizes on the performance of the algorithm. We find that while circuits with a large number of qubits fail more often under noise that depletes the qubit population, their outputs on a successful measurement are just as robust under Rydberg atom dissipation or qubit dephasing as smaller systems. However, such circuits might not perform as well under coherent multiqubit errors such as Rydberg atom crosstalk. We also find that the optimized parameters are especially robust to noise, suggesting that a noisier quantum system can be used to find the optimal parameters before switching to a cleaner system for measurements of observables.

Allen, James↗

Adaptive PID Gain Scheduling Control for Hydropower Turbine Using Neural CDE and Stochastic Distribution Shaping

This paper introduces a gain-scheduling PID controller design strategy for hydroturbine frequency control mode. This scheme first uses real data to learn the nonlinear dynamics of the hydroturbine using neural controlled differential equations and then perturbs the obtained nonlinear system at different equilibrium points, based on which a static output feedback adaptive dynamic programming algorithm is then used to optimize the PID gains for each equilibrium point. Moreover, a continuous-time version of stochastic distribution control is proposed to further fine-tune the optimized PID gains. Finally, the controller is obtained by implementing linear interpolation between the optimized PID control gains. The simulation results show that the proposed gain-scheduling PID controller can control a larger range of operation points compared with the given fixed PID controller and the baseline method. Compared with the given fixed PID controller, the proposed gain-scheduling PID controller can regulate hydroturbine frequency against disturbances induced by power-load variation with over 50% less overshoot for some operation points.

13 HYDRO ENERGY↗

Thermal Design Consideration of Medium Voltage High Frequency Transformers

This paper describes a thermal design using epoxy potting for high frequency transformer (HFT) for medium voltage high power electronic applications, namely emerging solid-state transformers. To address HFT’s thermal stress problem, thermal performance of the HFT with air and epoxy potting cooling method are both simulated by magnetothermal-fluid coupling method (MTFCM). The effectiveness of the simulation analysis is verified in a 100 kW DC-DC resonant converter. Influencing factors including distance between the high frequency transformer and the metal enclosure, epoxy thermal conductivity and heat sink placements are simulated and analyzed. Finally, an optimal cooling structure is proposed for future high power HFTs.

Wang, Haoming↗

Maximizing the Performance of a 3D Printed Heat Sink by Accounting for Anisotropic Thermal Conductivity During Filament Deposition

Newly developed high thermal conductivity polymer composite 3D printing filaments are used to characterize the thermal properties as a function of print orientation. The thermal conductivity of a printed part is anisotropic and varies by 2-6 times depending on the print direction - demonstrating higher conductivity in the deposition direction (in-plane) than in the two directions perpendicular to the deposition direction (cross-plane and through-plane). Therefore, deposition path planning greatly affects the overall heat dissipation rate and the performance of the heat sink. Traditionally, 3D printing slicers generate deposition paths based solely on geometric constraints. This work investigates a new approach of deposition path planning assisted by computational predictions of the heat sink thermal performance. The proposed approach uses a thermal simulation of a 3D-printed part, accounting for the anisotropic thermal properties, and the orientation of the local material properties are assigned based on the deposition path in multiple print orientations. The performances predicted via the simulations are compared, and the optimal deposition path is determined. For the highest thermal conductivity 3D printing filament (~ 12 W/m- K in-plane), a heat sink printed with the print direction parallel to the fins z-axis had ~20% improved performance in comparison to a heat sink with print direction perpendicular to the fins z-axis. Moreover, a plastic 3D printed heat sink was able to perform within 7% of an extruded Aluminum heat sink with similar geometry under natural convection. The computational predictions show the same trend as experimental measurements using 3D printed heat sinks.

Smith, Matthew↗

An Accurate Fast Fluid Dynamics Model for Data Center Applications

Traditional CFD is broadly useful in the design and operation of reliable and efficient data centers. It is nevertheless computationally expensive, particularly when employed for design optimization, which usually requires multiple simulations. To speed-up calculations while retaining essential airflow physics, several researchers have turned to an alternative CFD methodology, namely, fast fluid dynamics (FFD). FFD has been reported to be much faster than traditional CFD, but at an assumed (acceptable) trade-off of reduced accuracy. However, a recent comparison of FFD and traditional CFD for data center plenum applications produced nearly indistinguishable results and suggested that previously-reported FFD/traditional-CFD differences were due primarily to inconsistent: 1) advection schemes, 2) computational grids, and 3) turbulence models. The present paper extends this FFD/traditional-CFD comparison to the data center whitespace and confirms the finding that a true like-for-like comparison produces nearly identical predictions. Our FFD implementation utilizes a first-order upwind finite-volume scheme for advection like traditional CFD, so alternative advection schemes (e.g., semi-Lagrangian) are not considered further here. Likewise, the effect of grid choice on traditional-CFD predictions is well known. However, the effect of turbulence model for data center applications has not been reported extensively so we do consider this topic further here. We compare the standard k-e and a simpler algebraic model to benchmark experimental data from a real data center. We find that the algebraic turbulence model predicts rack-inlet temperatures at a similar level-of-accuracy as the k-e model for our fairly-simple-airflow reference data center.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Joint Fleet Sizing and Charging System Planning for Autonomous Electric Vehicles

This paper studies the joint fleet sizing and charging system planning problem for a company operating a fleet of autonomous electric vehicles (AEVs) for passenger and goods transportation. Most of the relevant published papers focus on intracity scenarios and adopt heuristic approaches, e.g., agent based simulation, which do not guarantee optimality. In contrast, we propose a mixed integer linear programming model for intercity scenarios. This model incorporates comprehensive considerations of 1) limited AEV driving range; 2) optimal AEV routing and relocating operations; 3) time-varying origin-destination transport demands; and 4) differentiated operation cost structure of passenger and goods transportation. The proposed model can be computational expensive when the scale of the transportation network is large. We then exploit the structure of this program to expedite its solution. Numerical experiments are conducted to validate the proposed method. Our experimental results show that AEVs in passenger and goods transportation have remarkable planning and operation differences. We also demonstrate that intelligent routing and relocating operations, charging system and vehicle parameters, e.g., charging power, battery capacity, driving speed etc., can significantly affect the economic efficiency and the planning results of an AEV fleet.

33 ADVANCED PROPULSION SYSTEMS↗

Correlated Production and Analog Transport of Fission Neutrons and Photons Using Fission Models FREYA, FIFRELIN, and the Monte Carlo Code TRIPOLI-4

Fission modeling in general-purpose Monte Carlo transport codes often relies on average nuclear data provided by the international evaluation libraries. As such, only average fission multiplicities are available and correlations between fission neutrons and photons are missing. Whereas uncorrelated fission physics is usually sufficient for standard reactor core and radiation shielding calculations, correlated fission secondaries are required for specialized nuclear instrumentation and detector modeling. For coincidence counting detector optimization for instance, the precise simulation of fission neutrons and photons that remain correlated in time from birth to detection is essential. New developments were recently integrated into the Monte Carlo transport code TRIPOLI-4 to model fission physics more precisely, the purpose being to access event-by-event fission events from two different fission models: Fission Reaction Event Yield Algorithm (FREYA) and Fission Fragment Evaporation Leading to an Investigation of Nuclear Data (FIFRELIN). TRIPOLI-4 simulations can now be performed, either by connecting via an application programming interface to the Lawrence Livermore National Laboratory fission library including FREYA, or by reading the external fission event data files produced by FIFRELIN beforehand. These new capabilities enable us to easily compare results from the Monte Carlo transport calculations using the two fission models in a nuclear instrumentation application. In the first part of this paper, the broad underlying principles of the two fission models are recalled. We then present the experimental measurements of neutron angular correlations for 252 Cf(sf) and 240 Pu(sf). The correlations were measured for several neutron kinetic energy thresholds. In the latter part of this paper, simulation results are compared with the experimental data. Spontaneous fissions in 252 Cf and 240 Pu are modeled by FREYA or FIFRELIN. Emitted neutrons and photons are subsequently transported to an array of scintillators by TRIPOLI-4 in the analog mode to preserve their correlations. Angular correlations between fission neutrons obtained independently from these TRIPOLI-4 simulations, using either FREYA or FIFRELIN, are compared with the experimental results. In conclusion, for 240 Pu(sf), the measured correlations were used to tune one of the FREYA model parameters. Similarly, different sets of parameters have been tested in FIFRELIN to try to improve the agreement with the experimental data.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Improving Longitudinal Impedance of Inductive Inserts

In this article, we discuss methods for improving the longitudinal impedance of inductive inserts used for space charge compensation in the proton storage ring (PSR) at Los Alamos Neutron Science Center (LANSCE) Facility. The PSR relies on inductive inserts to counteract the capacitive impedance produced by the space charge of the proton beam, which can otherwise lead to beam debunching. Existing inductive inserts, which utilize Toshiba M4C21A ferrite, have been effective but require heating to reduce resistive components that can cause instabilities at higher frequencies. This article explores the potential of alternative ferrite materials, specifically National Magnetics C2050, to enhance the performance of these inserts by reducing the need for heating and minimizing resistive impedance. Through both analytical modeling and CST Particle Studio simulations, it is demonstrated that optimized ferrite materials and geometries can significantly improve space charge compensation while mitigating instability risks, ultimately enhancing the performance of the PSR.

43 PARTICLE ACCELERATORS↗

Transfer of autocollimator calibration for use with scanning gantry profilometers for accurate determination of surface slope and curvature of state-of-the-art x-ray mirrors

X-ray optics, desired for beamlines at free-electron-laser and diffraction-limited-storage-ring x-ray light sources, must have almost perfect surfaces, capable of delivering light to experiments without significant degradation of brightness and coherence. To accurately characterize such optics at an optical metrology lab, two basic types of surface slope profilometers are used: The long trace profilers (LTPs) and nanometer optical measuring (NOM) like angular deflectometers, based on electronic autocollimator (AC) ELCOMAT-3000. The inherent systematic errors of the instrument's optical sensors set the principle limit to their measuring performance. Where autocollimator of a NOM-like profiler may be calibrated at a unique dedicated facility, this is for a particular configuration of distance, aperture size, and angular range that does not always match the exact use in a scanning measurement with the profiler. Here we discuss the developed methodology, experimental set-up, and numerical methods of transferring the calibration of one reference AC to the scanning AC of the Optical Surface Measuring System (OSMS), recently brought to operation at the ALS Xray Optics Laboratory. We show that precision calibration of the OSMS performed in three steps, allows us to provide high confidence and accuracy low-spatial-frequency metrology and not â€print into' measurements the inherent systematic error of tool in use. Here, with the examples of the OSMS measurements with a state-of-The-Art x-ray aspherical mirror, available from one of the most advanced vendors of X-ray optics, we demonstrate the high efficacy of the developed calibration procedure. The results of our work are important for obtaining high reliability data, needed for sophisticated numerical simulations of beamline performance and optimization of beamline usage of the optics. This work was supported by the U. S. Department of Energy under contract number DE-AC02-05CH11231.

47 OTHER INSTRUMENTATION↗

Inverse design of a pyrochlore lattice of DNA origami through model-driven experiments

Sophisticated statistical mechanics approaches and human intuition have demonstrated the possibility of self-assembling complex lattices or finite-size constructs. However, attempts so far have mostly only been successful in silico and often fail in experiment because of unpredicted traps associated with kinetic slowing down (gelation, glass transition) and competing ordered structures. Theoretical predictions also face the difficulty of encoding the desired interparticle interaction potential with the experimentally available nano- and micrometer-sized particles. To overcome these issues, we combine SAT assembly (a patchy-particle interaction design algorithm based on constrained optimization) with coarse-grained simulations of DNA nanotechnology to experimentally realize trap-free self-assembly pathways. In this paper, we use this approach to assemble a pyrochlore three-dimensional lattice, coveted for its promise in the construction of optical metamaterials, and characterize it with small-angle x-ray scattering and scanning electron microscopy visualization.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Advances in global optimization of high brightness beams

High brightness electron beams play an important role in accelerator-based applications such as driving X-ray free electron laser (FEL) radiation.Here in this paper, we report on advances in global beam dynamics optimization of an accelerator design using start-to-end simulations and a new parallel multi-objective differential evolution optimization method. The global optimization results in significant improvement of the final electron beam brightness.

43 PARTICLE ACCELERATORS↗

CGYRO

CGYRO is a multi-species spectral gyrokinetic solver with sonic rotation capability. It is an Eulerian GK solver specifically designed and optimized for collisional, electromagnetic, multiscale simulation. A key algorithmic aspect of CGYRO is the radially spectral formulation used to reduce the complicated integral gyroaveraging kernel into a multiplication in wavenumber space, but retaining the ability to treat profile variation important for edge plasmas. This makes it suitable for next-generation computational systems that require high levels of parallel concurrency.

Belli, EmilyA↗

TEAL

TEAL is a financial performance calculator plugin for the RAVEN code, framework, resolving around the computation of Net Present Value and associated financial metrics. TEAL can make use of inflation rates, taxation, escalation factors, capital expenditure economy of scale scaling factors. The unique feature of TEAL is the capability to be linked with RAVEN external models and build corresponding cash flows using the variables computed by those external models. In addition to be able to use the capability to generate cash flows derived from complex physical models generated by RAVEN, another distinctive feature of TEAL is the capability to provide financial risk/probabilistic metrics that can empower RAVEN to perform optimization/analysis driven by financial risk augmentations. Optimization, robust optimization, parametric studies, large parallel simulations, sensitivity analysis, data mining, etc. are just some of the capabilities that can be leveraged.

Alfonsi, Andrea↗

Deeper Fluids

A framework for building and training neural network surrogates for CFD simulation, with a focus on quickly simulating volume fraction dynamics for carbon capture efficiency prediction. Augments the Deep Fluids (Kim et al., 2019) approach as described in our paper "Latent Space Simulation for Carbon Capture Design Optimization" (Bartoldson et al., 2021).

Nguyen, PhuocChanh N↗

MAM4xx

SAND2023-05377O The Modal Aerosol Model with 4 fixed modes (MAM4xx) is a performance-portable, C++ implementation for modeling the chemical reactions of aerosols in the atmosphere. MAM4xx is part of the Energy Exascale Earth System Model (E3SM). The E3SM project is an ongoing Earth system modeling, simulation, and prediction project that optimizes the use of DOE laboratory resources to meet the science needs of the nation and the mission needs of DOE. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Overfelt, James↗

DEVELOPMENT OF PLASMA SIMULATION TOOL FOR THE ACCEL-ERATING CAVITIES

Plasma processing of superconducting radio frequency (SRF) cavities has shown an improvement in accelerating gradient by reducing the radiation due to field emission and multipacting. Plasma processing is a common tech-nique where the free oxygen produced by the plasma breaks down and removes hydrocarbons from surfaces. This increases the work function and reduces the second-ary emission coefficient. The hydrocarbon fragments of H2, CO, CO2, and H2O are removed from the system with the process gas which is flowing through the system. Here, we present COMSOL for the first time to simulate the plasma processing of an SRF cavity. In this work, we use Jefferson Lab's C75 SRF cavities design as our case study. Using simulation, we predict the condition of plasma ignition inside the SRF cavity. The simulation provides information about the optimal rf coupling to the cavity, mode for plasma ignition, choice of gas concen-tration, power, and pressure.

Senevirathne, I.↗

Xyce Parallel Electronic Simulator Users' Guide Version 6.8

This manual describes the use of the Xyce Parallel Electronic Simulator. Xyce has been de- signed as a SPICE-compatible, high-performance analog circuit simulator, and has been written to support the simulation needs of the Sandia National Laboratories electrical designers. This development has focused on improving capability over the current state-of-the-art in the following areas: Capability to solve extremely large circuit problems by supporting large-scale parallel com- puting platforms (up to thousands of processors). This includes support for most popular parallel and serial computers. A differential-algebraic-equation (DAE) formulation, which better isolates the device model package from solver algorithms. This allows one to develop new types of analysis without requiring the implementation of analysis-specific device models. Device models that are specifically tailored to meet Sandia's needs, including some radiation- aware devices (for Sandia users only). Object-oriented code design and implementation using modern coding practices. Xyce is a parallel code in the most general sense of the phrase$-$ a message passing parallel implementation $-$ which allows it to run efficiently a wide range of computing platforms. These include serial, shared-memory and distributed-memory parallel platforms. Attention has been paid to the specific nature of circuit-simulation problems to ensure that optimal parallel efficiency is achieved as the number of processors grows.

97 MATHEMATICS AND COMPUTING↗

Percheron Power – Small Business Vouchers Pilot Program Technology Assistance Report

Battelle Memorial Institute, under U.S. Department of Energy Contract No. DE-AC05-76RL01830 to operate the Pacific Northwest National Laboratory (PNNL), agreed to provide technology assistance to Percheron Power in accordance with Small Business Vouchers Pilot Program Agreement project number 68798. This agreement w as designed to provide assistance to Percheron Power in performing modeling and simulations for the design and optimization of composite Archimedes hydrodynamic screw (AHS) turbines.

42 ENGINEERING↗