Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation Cases”

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 865 records · Page 48

A Sparse and Low Rank Penalized Signal Decomposition Model with Constraints: Anomaly Detection in PV Systems

Recently, robust PCA has seen its wide application in various industries for its ability to perform the task of anomaly detection. The essence of robust PCA approach is to break down the signal into a low rank component and sparse component. In many applications, a simple breakdown of the signal without accounting for the signs of low rank components and sparse components would violate the physical constraints of the decomposed signal. In addition, often times, the signals in the real world collected for a long duration has smooth changes within a day and between days. As an example, the power signals collected in a photovoltaic (PV) system are cyclostationary, exhibiting these characteristics. Neglecting the smoothness of signals would result in miss detection of anomalous signals which are smooth within a day but non-smooth between days and vice versa. In this paper, we developed a signal decomposition approach for the purpose of anomaly detection based on the idea of low rank and sparse decomposition taking into consideration the signs of the decomposed low rank and sparse components and the within-day and between-day smooth changes in the original signals. The proposed unsupervised approach for fault detection eliminates the need for faulty samples required by other machine learning methods. It does not require the full I-V characteristics to work. Furthermore, there is no need for complex modelling of PV systems as in the case of power loss analysis. Using Monte Carlo simulations, we demonstrate the ability of our proposed approach for detecting anomalies of different duration and severity in PV systems.

14 SOLAR ENERGY↗

Integrated Coordination of Voltage Regulators with Distributed Cooperative Inverter Control in Systems with High Penetration of DGs

Increased penetration of distributed generation (DG) in recent years has a significant impact on operation and control of power systems. One of the main issues with the increased penetration of DGs is the effect of voltage rise. This paper investigates integration of coordinated voltage regulator (VRs) operation with cooperative distributed inverter control to maintain voltages in appropriate limits in distribution systems with high penetration of DGs. The coordinated VR control is used to optimize the settings of multiple nested VRs while the cooperative distributed control optimally dispatches reactive power of the DGs. The integrated approach is compared with the case without coordinated VRs control (uncoordinated local control). Simulation results illustrate the effective performance of the proposed method to achieve unified voltage profile.

coordinated voltage regulator control↗

Guaranteeing a Physically Realizable Battery Dispatch Without Charge-Discharge Complementarity Constraints

The non-convex complementarity constraints present a fundamental computational challenge in energy constrained optimization problems. In this work, we present a new, linear, and robust battery optimization formulation that sidesteps the need for battery complementarity constraints and integers and prove analytically that the formulation guarantees that all energy constraints are satisfied which ensures that the optimized battery dispatch is physically realizable. In addition, we bound the worst-case model mismatch and discuss conservativeness. In conclusion, simulation results further illustrate the effectiveness of this approach.

25 ENERGY STORAGE↗

Estimation of the Water Table Position in Unconfined Aquifers with MODFLOW 6

Abstract The numerical estimation of the position of the water table in unconfined aquifers is important for many practical applications. Its determination through observations or analytical methods is restricted to a few cases. Therefore, it is often estimated through numerical simulations, which may be affected by numerical artifacts and/or poor stability. We use MODFLOW to estimate the position of the water table for a seemingly simple example problem and demonstrate difficulties that can be faced when performing this kind of numerical simulation. We explain the causes for the numerical challenges that originate from the properties of the mathematical equations that must be solved. Based on the results of more than 600 steady‐state simulations, we show how the stability of the numerical solution can be affected by the values of physical parameters that define the problem (e.g., recharge rate, anisotropy ratio, and other parameters that control the numerical algorithm such as settings of the linear and nonlinear solution methods). Finally, we comment on some best practices to apply numerical simulations to estimate the water table position.

Geology↗

Bichromatic Wave Selection for Validation of the Difference-Frequency Transfer Function for the OC6 Validation Campaign: Preprint

The focus of the Offshore Code Comparison Collaboration, Continuation, with Correlation and unCertainity (OC6) project, which operates under the International Energy Agency (IEA) Wind Task 30, is to refine the accuracy of engineering tools used to design offshore wind turbines. In support of this work, a new validation campaign is being developed that seeks to better understand the nonlinear wave loading that excites floating wind systems at their low-frequency, rigid-body modes in surge and pitch. The validation data will be used in a three-way validation between simplified engineering tools and higher-fidelity tools, such as computational fluid dynamics (CFD). Irregular wave spectrums, which are traditionally used to examine the nonlinear wave interaction with offshore structures, are too computationally expensive to be simulated in CFD tools, and so bichromatic wave cases will be used instead. This paper reviews the process used to choose the bichromatic wave pairs to be used in the campaign to validate the second-order difference-frequency quadratic and potential loads at the surge and pitch natural frequencies of a floating semisubmersible.

17 WIND ENERGY↗

Leveraging on-shell interference to search for FCNCs of the top quark and the Z boson

Abstract Flavour-changing-neutral currents (FCNCs) involving the top quark are highly suppressed within the Standard Model (SM). Hence, any signal in current or planned future collider experiments would constitute a clear manifestation of physics beyond the SM. We propose a novel, interference-based strategy to search for top-quark FCNCs involving the Z boson that has the potential to complement traditional search strategies due to a more favourable luminosity scaling. The strategy leverages on-shell interference between the FCNC and SM decay of the top quark into hadronic final states. We estimate the feasibility of the most promising case of anomalous tZc couplings using Monte Carlo simulations and a simplified detector simulation. We consider the main background processes and discriminate the signal from the background with a deep neural network that is parametrised in the value of the anomalous tZc coupling. We present sensitivity projections for the HL-LHC and the FCC-hh. We find an expected 95% CL upper limit of $${\mathcal {B}}_{\textrm{excl}}(t\rightarrow Zc) = 6.4 \times 10^{-5}$$ B excl ( t → Z c ) = 6.4 × 10 - 5 for the HL-LHC. In general, we conclude that the interference-based approach has the potential to provide both competitive and complementary constraints to traditional multi-lepton searches and other strategies that have been proposed to search for tZc FCNCs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Energy Efficiency and Performance Evaluation of Self-Contained, Medium Temperature, Reach-In Refrigerated Display Cases

This project was part of an effort by Commonwealth Edison Company (ComEd) to evaluate the energy and peak demand saving potential of emerging technologies in the Chicago area. This document focuses on the assessment of energy-efficient, medium-temperature, self-contained refrigerated display cases utilizing environmentally friendly refrigerants. The results of this evaluation will be considered by CLEAResult to develop a new energy efficiency rebate measure for ComEd's incentive programs. This rebate measure will become an addition to the Technical Reference Manual (TRM). In 2016, the United Nations passed the Kigali Montreal Protocol Amendment which placed restrictions on certain types of refrigerants. In compliance with this amendment, the US Environmental Protection Agency (EPA) placed a ban on the manufacture of refrigeration systems using hydrofluorocarbons including R134a starting in January 2020. Although the ban has halted manufacture, the EPA continues to allow the use of these refrigerants. Therefore, it is critical to provide incentives for replacing these refrigerants with other environmentally friendly and energy-efficient alternatives. The energy-efficient refrigerator cases evaluated here (referred to as EE Case A and B) contain environmentally-friendly refrigerants in compliance with the EPA hydrofluorocarbon ban. These consist of natural refrigerant propane (R290), and HFC drop-in hydrofluoroolefin R513a, respectively. These cases also contain other energy-efficient components including efficient lighting, more robust evaporator and condenser fans and different-sized heat exchangers. EE Case A is also built with materials that better insulate the case, which improves energy efficiency. To ascertain the energy efficiency contribution of these design components, the consumption of the evaporator and condenser fan motors, compressor, and lighting/controller were evaluated individually. The medium-temperature, self-contained reach-in refrigerated display case was selected due to its widespread use in convenience stores and small supermarkets. Self-contained refrigeration has also seen increased use in restaurants due to curbside pickup during the COVID-19 pandemic. For this study, the refrigerated display cases' performance was evaluated in a controlled environmental chamber at representative indoor dry-bulb and humidity conditions found in supermarkets within ComEd's service territory climate zone. The test method used in this project was foundationally inspired by the ANSI/ASHRAE 72-2018 method of testing. However, modifications to the ANSI/ASHRAE test methodology were implemented to better represent customer operation of the units. In addition to the indoor supermarket conditions, the cases were also evaluated at the "upper target," or environmental conditions used in ANSI/ASHRAE 72-2018. The case total power and case components were metered to obtain their daily power (kW) and energy consumption(kWh). The cases were filled with thermal filler mass to replicate thermal mass of product loading. Additionally, product simulators were used to provide product temperature information. Door actuators were mounted to each of the cases' three doors to replicate regular door openings and effects of shopper traffic.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Final Report on Development of Microwave Resonant Cavity Transducer for Fluid Flow Sensing: Development of Microwave Cavity Flow Meter for Advanced Reactor High Temperature Fluids

We are investigating a microwave resonant cavity transducer for flow sensing in the vessel of a high temperature fluid advanced reactor (AR), such as a molten salt cooled reactor (MSCR) or a sodium fast reactor (SFR). This transducer is a hollow metallic cylindrical cavity, with one of the flat walls of the cylinder flexible enough to undergo microscopic deflection due to dynamic fluid pressure. Membrane deflection leads to a shift in the microwave resonant frequency, which can be detected with a spectrum analyzer. We have developed a continuum electromechanics model of the microwave resonant cavity transducer performance. The deflection of the membrane is calculated using equations from the theory of plates and shells for deflection of radially constrained thin circular disk under uniform load. The microwave frequency shift is calculated using the equations from the microwave cavity perturbation theory. The transducer model was validated with proof-of-concept experimental data of water flow sensing with a Brass 360 right circular cylinder with 22.2mm diameter and 203µm thick wall. The cavity was excited through a subwavelength hole in TE 011 mode with resonant frequency f ≈ 17.8GHz. Subsequently, we have performed a preliminary proof-of-principle test of flow sensing in high temperature liquid sodium in environment. For this test, we have developed a cylindrical resonator with the same dimensions as for the Brass 360 cavity. The resonator for high temperature liquid sodium test was machined from stainless steel 316 and electroplated with silver on the interior surfaces. We have also developed and insertion probe consisting of a 50cm WR-42 brass waveguide enclosed in a protective SS316 tube. The cavity was excited through a subwavelength hole on the side of the wall of the cylinder in the TE 011 mode with resonant frequency f ≈ 17.8GHz. The liquid sodium setup consists of a cylindrical vessel with a center feed line, where a transducer inserted through the top cap of the vessel measures velocity of the impinging liquid jet. The sodium flow sensing study was performed in a liquid sodium vessel at 340°C temperature and ambient pressure. The flow rate was changed by varying sodium pump power, and a frequency shift of several hundred kHz was observed. The transducer was removed after 70 days of testing. No visual damage was observed, and the electromagnetic response remained the same. As a calibration experiment, we have assembled a water vessel with a center feed, and with dimensions similar to those of the liquid sodium setup. Water flow sensing was performed before and after liquid sodium test. The responses in both cases are similar and agree with COMSOL computer simulations. Because long-term (multi-year) experimental tests of transducer resilience to harsh environment are not practical, we have developed a probabilistic model of creep to estimate transducer resilience to the harsh environment. The probabilistic model considers diffusion creep under the condition of high temperature and low stress, where the stress and temperature are allowed to be random variables with Gaussian distributions. Using the probabilistic model, we estimate inelastic membrane deflections due to creep for several temperature ranges. We conclude that for temperatures less than 650°C, creep has negligible long-term effect on the transducer performance.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

An Approach for Simulating Transport of Nonspherical, Sublimating Particles in Rarefied Gas Flows via the Direct Simulation Monte Carlo Method

This paper extends the particle Green’s function approach of Gallis et al. [Phys. Fluids 13, 3482–3492 (2001)] to compute the force and heat transfer on nonspherical, sublimating particles in rarefied gas flows. Green’s functions are derived for rods and plates, which are representative of ice particle shapes, and a numerical procedure is developed for computing the Green’s functions for a general discretized convex shape. Additionally, a sublimation model is developed to simulate the evolution of ice particles in rarefied gas flows. These models are implemented in a direct simulation Monte Carlo code and demonstrated for the case of ice-particle-laden hypersonic flow over a ramp.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Two-DoF Model-Informed Controller Gain Tuning for Several Floating Wind Platforms

Much effort has recently been put into the development of collective blade pitch controllers for floating offshore wind turbines, with the aim of overcoming negative damping issues that arise with traditional control methods. One proposed approach to this challenge involves using a two-degree-of-freedom model to inform the gain schedule of a nacelle velocity feedback term in an otherwise conventional proportional-integral controller. The model uses tower-top fore-aft and rotor angular displacements, and is used to calculate a nacelle velocity feedback gain that results in a specified increase in platform pitch damping. Earlier performance evaluations of this tuning method were favorable, suggesting its potential as an easy way for researchers to obtain an adequate controller. This paper expands on those previous results by examining the performance of the tuning method relative to baseline controllers for several hull configurations, and for several prescribed increases in platform pitch damping. Simulations were run in OpenFAST for several load cases above rated wind speed and show results consistent with trends in the earlier study. The tuning method is thus shown to be adaptable to many different types of hulls, making it useful for the evaluation of prototype designs.

artificial intelligence↗

Considerations on wakefield effects in a VUV FELO driven by a superconducting tesla-type linac

The effects on beam dynamics from long-range and short-range wakefields from TESLA-type cavities are considered in regard to a proposed FEL oscillator (FE-LO) operating at 120 nm. This would be driven by the Fermilab Accelerator Science and Technology (FAST) linac at 300 MeV with a 3-MHz micropulse repetition rate. Our wakefield studies showed measurable effects on submacropulse centroid stability and on submicropulse head-tail kicks that can lead to emittance degradation. In the case of the former, we use MINERVA/OPC to simulate the ~100-$\mu$m centroid slew effects on the saturated output power levels of the FELO.

43 PARTICLE ACCELERATORS↗

Portable Programming Model Exploration for LArTPC Simulation in a Heterogeneous Computing Environment: OpenMP vs. SYCL

The evolution of the computing landscape has resulted in the proliferation of diverse hardware architectures, with different flavors of GPUs and other compute accelerators becoming more widely available. To facilitate the efficient use of these architectures in a heterogeneous computing environment, several programming models are available to enable portability and performance across different computing systems, such as Kokkos, SYCL, OpenMP and others. As part of the High Energy Physics Center for Computational Excellence (HEP-CCE) project, we investigate if and how these different programming models may be suitable for experimental HEP workflows through a few representative use cases. One of such use cases is the Liquid Argon Time Projection Chamber (LArTPC) simulation which is essential for LArTPC detector design, validation and data analysis. Following up on our previous investigations of using Kokkos to port LArTPC simulation in the Wire-Cell Toolkit (WCT) to GPUs, we have explored OpenMP and SYCL as potential portable programming models for WCT, with the goal to make diverse computing resources accessible to the LArTPC simulations. In this work, we describe how we utilize relevant features of OpenMP and SYCL for the LArTPC simulation module in WCT. We also show performance benchmark results on multi-core CPUs, NVIDIA and AMD GPUs for both the OpenMP and the SYCL implementations. Comparisons with different compilers will also be given where appropriate.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The accuracy of multi-group models for nonlocal electron transport in magnetized plasmas

In the extreme conditions of inertial confinement fusion experiments, heat flow plays a vital role, but local diffusive models frequently break down and overestimate the heat flow. The situation becomes more complicated again in the significant magnetic fields generated during laser–plasma interactions or in magnetized fusion schemes. Accurate non-local and magnetized heat flow computations can be carried out using Vlasov–Fokker–Planck (VFP) simulations, but these are computationally expensive. There is, therefore, significant interest in using faster multi-group models to accurately calculate the non-local heat flow in magnetized plasmas. We benchmark two such multi-group models for calculating the heat flow, M1 and hybrid-AWBS-BGK, against diffusive models and full VFP simulations, before applying the models to realistic example test cases, both magnetized and unmagnetized. We find that the multi-group models generally perform very well for moderate non-localities up to kλmfp∼0.01, but the computational cost increases dramatically. hybrid-AWBS-BGK performs more effectively than M1 at high non-localities, up to kλmfp∼1, due to its adaptive solver and robust P1 closure, but tends to fail in very strong magnetic fields. Both codes are much faster than VFP simulations but are still slow in steep temperature gradients.

Arran, C. (ORCID:0000000286448118)↗

Unraveling the implications of finite specimen size on the interpretation of dynamic experiments for polycrystalline aluminum through direct numerical simulations

Normal and Pressure-shear plate impact (NPI and PSPI) experiments are popular experimental techniques for studying the mean-field macroscopic behavior of polycrystalline metals under high-rate dynamic loading. However, since both configurations rely upon geometry for subjecting the specimen to high strain rates, these experiments often involve a limited specimen size. Moreover, because of the inherent heterogeneities present within polycrystalline metals, it is difficult to ascertain if the size of the specimen and/or regions where measurements are made are sufficiently large for making representative inferences about the mean-field macroscopic properties from single-point velocity measurements. In the present study, we quantify the expected measurement variability on observable point measurements in NPI and PSPI experiments by carrying out direct numerical simulations (DNS) of statistically representative polycrystalline microstructures subjected to dynamic compression and compression-shear loading. In particular, we consider the role of specific material heterogeneities (e.g. the grain-to-grain difference in size, crystallographic orientation) on dispersion in the normal and transverse particle velocity records and on local fluctuations in key state variables (e.g. velocity, accumulated plastic strain) by incorporating these effects directly into a representative synthetic microstructure geometry and crystalline description of pure polycrystalline aluminum. The form of the present study is a large parametric investigation, consisting of ten ensembles of one hundred simulations. Each of the thousand simulations reflects a randomly realized synthetic microstructure in one of five cases of decreasing average grain size for the two loading configurations. Our analysis of the DNS results demonstrates that for both of these experimental configurations, the grain size directly correlates with the coefficient of variation (CV) in simulated point measurements, showing a convergent decrease in CV to zero (i.e. particle velocity record approaches the mean-field value) with decreasing grain size. Remarkably, the magnitude of variations in the particle velocity record is shown to be largest where the deviatoric stresses are most significant. In the case of NPI, this occurs at the elastic and plastic wavefront, whereas, in the case of PSPI, the magnitude of fluctuations are approximately constant throughout the experimental window time. The reasoning for the scatter in particle velocity due to the heterogeneous microstructure is demonstrated to be dependent on the mechanisms for accommodating deformation and on the interaction of reflection waves generated at sites of heterogeneities occurring at the scale of grains. Lastly, we develop a power-law description for the magnitude of scattering versus characteristic length, which provides a statistical framework for assessing the required number of grains per characteristic specimen dimension for minimizing scatter within these two experimental configurations (NPI, PSPI).

36 MATERIALS SCIENCE↗

Reservoir Dynamics in Proposed Operational Scenarios Where CO 2 -EOR Fields are Transitioned From CO 2 -Flood Enhanced Oil Recovery to Dedicated Carbon Storage: A Field Case Study

This study models the transition of CO 2 -enhanced oil recovery (CO 2 -EOR) fields to dedicated carbon storage using a generalized reservoir model informed by real field data, SACROC. Through scenario-based simulations, the work evaluates how reservoir depletion levels, boundary conditions, fluid properties, and domain size affect pressure buildup, CO 2 plume migration, and long-term containment performance. The results inform practices for repurposing oil fields into reliable CO 2 storage sites.

02 PETROLEUM↗

Thermal and Well Flow Performance of Closed-Loop Geothermal in the Wattenberg Area: Preprint

Closed-loop geothermal systems provide an alternative to resource-constrained hydrothermal systems and stimulation-intensive enhanced geothermal systems. In this work, we apply the slender-body theory (SBT) model, to simulate the well flow and heat transfer performance of U-loop well designs drilled in the Wattenberg area of the Denver-Julesburg Basin. Three U-loop well patterns are investigated including a single, double, and multi-lateral design. The subsurface within area is characterized by deep, hot (> 200degreesC) igneous/metamorphic basement rock underlying multiple sedimentary formations. The lateral section(s) of the U-loop lies within a target depth of 6 km where temperatures are estimated to approach 300degreesC. As a base case, conduction-only heat transfer is investigated through simulations with the SBT model within U-loops with open-hole laterals that exchange heat directly with the hot dry rock using water as a working fluid. The utilization of supercritical CO2 as a heat transfer fluid is also considered. For each scenario, the system performance in terms of annual heat production and temperature profile over a 20-year project lifetime are assessed. Also, the levelized costs of heat and electricity (LCOH and LCOE) are determined using a top-down technoeconomic analysis model. The results show that the performance and cost optimized U-loop design is one having an injection-production well spacing of 1,000 meters with ten 50-meter-spaced laterals that traverse a subsurface system with a temperature gradient of 60degreesC/km. By injecting 20 degreesC-water at a rate of 60 kg/s through this loop, an average heat production of 19 MWth can be achieved, resulting in an LCOE and LCOH of $136/MWh and $1.53/GJ, respectively, over a 20-year project life.

closed-loop geothermal↗