Demand bidding vs. demand response for industrial electrical loads
Not Available
Tang, Xin (ORCID:0009000720910131)↗
SEARCH · Search NASA
Results for “test case”
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.
Adaptive mesh refinement for steady flows in Nek5000
Adaptive mesh refinement is performed in the framework of the spectral element method augmented by approaches to error estimation and control. The h-refinement technique is used for adapting the mesh, where selected grid elements are split by a quadtree (2D) or octree (3D) structure. Continuity between parent-child elements is enforced by high-order interpolation of the solution across the common faces. Parallel mesh partitioning and grid management respectively, are taken care of by the external libraries ParMETIS and p4est. Two methods are considered for estimating and controlling the error of the solution. The first error estimate is local and based on the spectral properties of the solution on each element. This method gives a local measure of the L-2-norm of the solution over the entire computational domain. The second error estimate uses the dual-weighted residuals method - it is based on and takes into account both the local properties of the solution and the global dependence of the error in the solution via an adjoint problem. The objective of this second approach is to optimize the computation of a given functional of physical interest. The simulations are performed by using the code Nek5000 and three steady-state test cases are studied: a two-dimensional lid-driven cavity at Re = 7, 500, a two-dimensional flow past a cylinder at Re = 40, and a three-dimensional lid-driven cavity at Re = 2,000 with a moving lid tilted by an angle of 30 degrees. Finally, the efficiency of both error estimators is compared in terms of refinement patterns and accuracy on the functional of interest. In the case of the adjoint error estimators, the trend on the error of the functional is shown to be correctly represented up to a multiplicative constant.
42 ENGINEERING↗
A reconstructed discontinuous Galerkin method for compressible flows in Lagrangian formulation
In this work, a high-order accurate reconstructed discontinuous Galerkin (rDG) method is developed for solving two-dimensional hydrodynamic problems in cell-centered updated Lagrangian formulation. This method is the Lagrangian limit of the unsplit rDG-ALE formulation, and is obtained by assuming the equality of the grid velocity to the fluid velocity only at cell boundaries. The conservative variables and the Taylor basis defined on the time-dependent moving mesh, provide the piece-wise polynomial expansion in the updated Lagrangian formulation. A multi-directional nodal Riemann solver is implemented for computing the grid velocity at the vertices and the numerical flux at the cell boundaries. A characteristic limiting procedure is extended from the primitive variable version to the conservative variable version, and its performance is compared with the limiter on physical variables. A number of benchmark test cases are conducted to assess the accuracy, robustness, and non-oscillatory property of the DG(P0), DG(P1) and rDG(P1P2) methods. The numerical experiments demonstrate that the developed rDG method is able to attain the designed order of accuracy and the characteristic limiting procedure outperforms the limiter on physical variables in terms of the monotonicity and symmetry preservation for shock problems.
42 ENGINEERING↗
Comparison of Surface Tension Generation Methods in Smoothed Particle Hydrodynamics for Dynamic Systems
Developing robust numerical models of dynamic surface tension dominated multiphase systems is an ongoing challenge, especially in scenarios with large density and viscosity ratios. This is critical to the design and understanding of various physical and engineering systems, such as fluidized beds, fuel injectors, and drug delivery schemes. Much of the computational work in surface tension dominated multiphase flows has employed the continuum surface force method (CSF) of Brackbill et al. [1], which recasts surface tension from a surface force to a volumetric force that can be imposed in the vicinity of an interface. The CSF method produces accurate results across a variety of systems, however it relies on the identication of surface normals, which can be unreliable under certain conditions. Alternative methods of simulating surface tension have been proposed. Here the advantages and disadvantages of the CSF method in comparison to a pairwise forces (PF) method proposed by Tartakovsky and Panchenko [2] are explored. The CSF and PF methods are used in a smoothed particle hydrodynamics (SPH) framework to model dynamic systems. Results are compared to existing test cases from the literature and to analytic solutions derived from fundamental normal mode behavior of bubbles and droplets. A physical system for which the PF method is more stable and physically appropriate than the CSF method is identied.
smoothed particle hydrodynamics, multiphase flow↗
3D Cell-centered hydrodynamics with subscale closure model and multi-material remap
We extend a higher-order finite volume cell-centered hydrodynamic (CCH) formulation to include an interface-aware subscale closure model and a multi-material remap for simulating 3D compressible hydrodynamic problems within an arbitrary Lagrangian-Eulerian (ALE) framework. This CCH formulation involves a multidirectional approximate Riemann solution using quadratic polynomial reconstructions of the stress tensor and the velocity. At the subscale level, we determine pair-wise material interactions by solving a distinct approximate Riemann problem at the common interface, using the volume of fluids (VOF) method to find the interface. Material interactions are constrained to ensure smooth pressure equilibration among materials. The accuracy and robustness of the ALE method is demonstrated by simulating a suite of 3D Cartesian multi-material problems covering both gas and solid dynamics, where each test case has two or more materials.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
A particle resolved simulation approach for studying shock interactions with moving, colliding solid particles
This work applies a new combination of techniques for the fully resolved simulation of compressible, gas–particle multiphase flows. Here, the adaptive wavelet collocation method is used to dynamically, and efficiently, adapt the computational grid to localized flow features and the particles. A characteristic-based volume penalization method that imposes arbitrary Dirichlet, Neumann, or Robin-type immersed boundary conditions, is used to enforce the no-slip condition at particle surfaces. A hard-sphere collision model is applied to capture the particle–particle collisions. Proof of concept test cases are presented, showcasing the dynamic grid adaptation and fully resolved two-way coupling between the phases that is possible with this approach. Results for a shock-driven single cylinder under viscous and inviscid conditions are presented along with a demonstration of a shock interacting with a cloud of randomly distributed cylinders and spheres.
42 ENGINEERING↗
Highly-scalable GPU-accelerated compressible reacting flow solver for modeling high-speed flows
Emerging supercomputing systems utilize a combination of central processing units (CPUs) and graphics processing units (GPUs) in an effort to reach exascale capabilities while minimizing the energy footprint of operating such systems. Such heterogeneous machines introduce new challenges for fluids solvers because the hardware architecture and operation of a GPU are fundamentally different from conventional CPUs. In this work, a general approach for efficient implementation of finite-volume based reacting flow solvers on such heterogeneous systems is presented. Three main challenges, namely, data access pattern, thread divergence, and thread safety, are addressed. Since compressible reacting flows require special methods to deal with chemical reactions, hyperbolic and nonlinear convection terms, and the presence of turbulence, specific algorithms that ensure GPU-based efficiency are developed. The approach is demonstrated on the widely available OpenFOAM open source software by modifying core algorithms for GPU accessibility. The scalability of the resulting solver, is demonstrated using practical test cases, including flow through a scramjet engine and the dynamics of a rotating detonation engine. Here, the solver provides near-ideal scaleup on a large number of GPUs (>3000), and extremely efficient use of the GPUs, with throughput nearly a constant even when processing a large number of control volumes.
42 ENGINEERING↗
3D mesh regularization within an ALE code using a weighted line sweeping method
The Lagrangian formalism is widely used to simulate hydrodynamic responses in complex engineering applications, particularly those involving strong shock waves. However, as the mesh moves with the fluid, it can become highly distorted, requiring a regularization step. This involves constructing a new grid and remapping conservative quantities onto it to restore mesh quality. This work introduces a regularization method for block-structured meshes within a 3D ALE (Arbitrary Lagrangian-Eulerian) code. The proposed approach prevents mesh tangling while preserving the anisotropic features of the initial Lagrangian mesh. This regularization technique incorporates aspect ratio-based weights to control mesh smoothing. Unlike uniform rezoning techniques, this weighted approach maintains proximity to the Lagrangian mesh while improving mesh quality. Here, the method effectively handles concave geometries by mitigating the grid attraction phenomenon, which typically leads to mesh concentration along concave edges. Numerical experiments demonstrate its efficiency in regularizing severely deformed meshes, and its integration within the ALE framework is validated on challenging hydrodynamic test cases, including the triple point problem.
42 ENGINEERING↗
Accelerating dissipative particle dynamics simulations on GPUs: Algorithms, numerics and applications
We present a scalable dissipative particle dynamics simulation code, fully implemented on the Graphics Processing Units (GPUs) using a hybrid CUDA/MPI programming model, which achieves 10-30 times speedup on a single GPU over 16 CPU cores and almost linear weak scaling across a thousand nodes. A unified framework is developed within which the efficient generation of the neighbor list and maintaining particle data locality are addressed. Our algorithm generates strictly ordered neighbor lists in parallel, while the construction is deterministic and makes no use of atomic operations or sorting. Such neighbor list leads to optimal data loading efficiency when combined with a two-level particle reordering scheme. A faster in situ generation scheme for Gaussian random numbers is proposed using precomputed binary signatures. We designed custom transcendental functions that are fast and accurate for evaluating the pairwise interaction. The correctness and accuracy of the code is verified through a set of test cases simulating Poiseuille flow and spontaneous vesicle formation. Computer benchmarks demonstrate the speedup of our implementation over the CPU implementation as well as strong and weak scalability. A large-scale simulation of spontaneous vesicle formation consisting of 128 million particles was conducted to further illustrate the practicality of our code in real-world applications.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Extension of Landau-fluid closure to weakly collisional plasma regime
Both collisionless and collisional Landau fluid (LF) closure have been developed and implemented in two-fluid plasma simulation framework BOUT++, taking full advantage of the fast non-Fourier method (Dimits, 2014). The spectral range of good fit can be conveniently extended to larger domain by adding more Lorentzians. The scaling for upper limit of wavenumber resolved in collisionless LF operator and the fitting coefficients for the collisional LF operator, the later is extended to weakly collisional regime (ν$^{*}_{e}$ ~ 0.01), are presented in this paper. The flux limited expression q$^{FL}_{∥e}$ is reduced to free streaming expression q$^{FS}_{∥e}$ in collisionless limit and gives classical Spitzer–Härm result q$^{SH}_{∥e}$ in collisional limit as we expected. The collisional LF closure recovers the classical Spitzer–Härm result q$^{SH}_{∥e}$ in collisional limit and is same as the collisionless LF closure in collisionless limit. The test cases also obviously reflect the nonlocal effects from LF closures and responding to different boundary conditions in open surface region. Our report will make it more applicable and possible to include the Landau damping kinetic effect in the fluid simulation of tokamak plasma.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
ECsim-CYL: Energy Conserving Semi-Implicit particle in cell simulation in axially symmetric cylindrical coordinates
Based on the previously developed Energy Conserving Semi Implicit Method (ECsim) code, we present its cylindrical implementation, called ECsim-CYL, to be used for axially symmetric problems. The main motivation for the development of the cylindrical version is to greatly improve the computational speed by utilizing cylindrical symmetry. The ECsim-CYL discretizes the field equations in two-dimensional cylindrical coordinates using the finite volume method. For the particle mover, it uses a modification of ECsim’s mover for cylindrical coordinates by keeping track of all three components of velocity vectors, while only keeping radial and axial coordinates of particle positions. In this paper, we describe the details of the algorithm used in the ECsim-CYL and present a series of tests to validate the accuracy of the code including a wave spectrum in a homogeneous plasma inside a cylindrical waveguide and free expansion of a spherical plasma ball in vacuum. Finally, the ECsim-CYL retains the stability properties of ECsim and conserves the energy within machine precision, while accurately describing the plasma behavior in the test cases.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
A hybrid nodal-staggered pseudo-spectral electromagnetic particle-in-cell method with finite-order centering
Electromagnetic particle-in-cell (PIC) codes are widely used to perform computer simulations of a variety of physical systems, including fusion plasmas, astrophysical plasmas, plasma wakefield particle accelerators, and secondary photon sources driven by ultra-intense lasers. In a PIC code, Maxwell's equations are solved on a grid with a numerical method of choice. This article focuses on pseudo-spectral analytical time-domain (PSATD) algorithms and presents a novel hybrid PSATD PIC scheme that combines the respective advantages of standard nodal and staggered methods. The novelty of the hybrid scheme consists in using finite-order centering of grid quantities between nodal and staggered grids, in order to combine the solution of Maxwell's equations on a staggered grid with the deposition of charges and currents and the gathering of electromagnetic forces on a nodal grid. The correctness and performance of the novel hybrid scheme are assessed by means of numerical tests that employ different classes of PSATD equations in a variety of physical scenarios, ranging from the modeling of electron-positron pair creation in vacuum to the simulation of laser-driven and particle beam-driven plasma wakefield acceleration. It is shown that the novel hybrid scheme offers significant numerical and computational advantages, compared to purely nodal or staggered methods, for all the test cases presented.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
An immersed interface method for microstructure-scale electrochemical battery models: numerical formulation and performance portable implementation
We present the numerical formulation, verification, and performance portable implementation of an immersed interface method for microstructure scale electrochemical modeling of batteries. The innovation in this approach is the resolution of chemical species and electrostatic potential discontinuities at active interfaces without the use of interface conforming unstructured grids. A unified formulation on Cartesian grids for all domains (electrodes and electrolyte) is used with interfacial flux conditions applied using volume fraction or “color” function gradients. We have developed one dimensional and two dimensional test cases with analytic solutions for electrochemical modeling using which we verified the consistency and accuracy of our scheme. Our solver is also validated against solutions from a macroscale model and an unstructured multi-subdomain solver for a full lithium ion cell. We then demonstrated the utility of our solver on an image-based complex battery electrode microstructure at high charging rate. Our technique also exhibits good scalability on distributed memory architectures using central processing units (CPU), with problem sizes up to 1.8 billion degrees of freedom and with 5400 ranks. Initial performance studies of our open-source performance portable solver showed about 70 times speed up using a graphics processing unit (GPU) compared to single compute core for a problem with 4 million cells.
25 ENERGY STORAGE↗
Putting social impact assessment to the test as a method for implementing responsible tourism practice
The discourse on the social impacts of tourism needs to shift from the current descriptive critique of tourism to considering what can be done in actual practice to embed the management of tourism's social impacts into the existing planning, product development and operational processes of tourism businesses. A pragmatic approach for designing research methodologies, social management systems and initial actions, which is shaped by the real world operational constraints and existing systems used in the tourism industry, is needed. Our pilot study with a small Bulgarian travel company put social impact assessment (SIA) to the test to see if it could provide this desired approach and assist in implementing responsible tourism development practice, especially in small tourism businesses. Our findings showed that our adapted SIA method has value as a practical method for embedding a responsible tourism approach. While there were some challenges, SIA proved to be effective in assisting the staff of our test case tourism business to better understand their social impacts on their local communities and to identify actions to take. - Highlights: • Pragmatic approach is needed for the responsible management of social impacts of tourism. • Our adapted Social impact Assessment (SIA) method has value as a practical method. • SIA can be embedded into tourism businesses existing ‘ways of doing things’. • We identified challenges and ways to improve our method to better suit small tourism business context.
54 ENVIRONMENTAL SCIENCES↗
A knowledge brokerage approach for assessing the impacts of the setting up young farmers policy measure in Greece
This study explores Knowledge Brokerage (KB) aspects of an ex-post Impact Assessment (IA) for the Rural Development Programme (RDP) measure of setting up young farmers, under the Common Agricultural Policy (CAP), at the regional level in Northern Greece. The measure supports the entry of young farmers in agriculture by moving land from older to younger farmers. The aim of the study was to test a set of KB tools for improving the interaction between researchers and policy makers. Our analysis mainly focused on a suite of IA Support Modules to guide practitioners, and on a technical tool kit, a web-based contextualisation platform, to support the IA of the specific test case. Offering a structured approach towards IA, both the Support Modules and LIAISE-KIT allow framing the context, organisation, scheduling and method selection in the light of KB objectives. The evaluation of how IA Support Modules influence the Science Policy Interface (SPI), in the case of the ex-post assessment, demonstrated the high relevance of KB activities for facilitating the interaction between researchers and regional policy makers. The assessment bridges the gap between knowledge producers developing scientific output to be applied in a specific context, and knowledge users, who want clear messages regarding the policy challenges they face. Other conclusions include the need for specific guidelines and training for knowledge users, especially with regard to the use of tools. According to our findings, a consequent application of KB activities is a crucial pre-condition for successfully implementing IAs in future RDP measures.
54 ENVIRONMENTAL SCIENCES↗
Statistical change detection of building energy consumption: Applications to savings estimation
The surge in interval meter data availability and associated activity in energy data analytics has inspired new interest in advanced methods for building efficiency savings estimation. Statistical and machine learning approaches are being explored to improve the energy baseline models used to measure and verify savings. One outstanding challenge is the ability to identify and account for operational changes that may confound savings estimates. In the measurement and verification (M&V) context, ‘non-routine events’ (NREs) cause changes in building energy use that are not attributable to installed efficiency measures, and not accounted for in the baseline model’s independent variables. In the M&V process NREs must be accounted for as ‘adjustments’ to appropriately attribute the estimated energy savings to the specific efficiency interventions that were implemented. Currently this is a manual and custom process, conducted using professional judgment and engineering expertise. As such it remains a barrier in scaling and standardizing meter-based savings estimation. Here in this work, a data driven methodology was developed to (partially) automate, and therefore streamline the process of detecting NREs in the post-retrofit period and making associated savings adjustments. The proposed NRE detection algorithm is based on a statistical change point detection method and a dissimilarity metric. The dissimilarity metric measures the proximity between the actual time series of the post-retrofit energy consumption and the projected baseline, which is generated using a statistical baseline model. The suggested approach for NRE adjustment involves the NRE detection algorithm, the M&V practitioner, and a regression modeling algorithm. The performance of the detection and adjustment algorithm was evaluated using a simulation-generated test data set, and two benchmark algorithms. Results show a high true positive detection rate (75%-100% across the test cases), higher than ideal false positive detection rates (20%-70%), and low errors in energy adjustment (<0.7%). These results indicate that the algorithm holds for helping M&V practitioners to streamline the process of handling NREs. Moreover, the change point algorithm and underlying statistical principles could prove valuable for other building analytics applications such as anomaly detection and fault diagnostics.
29 ENERGY PLANNING, POLICY, AND ECONOMY↗
Energy performance evaluation of the ASHRAE Guideline 36 control and reinforcement learning–based control using field measurements
This study evaluates the energy performance of ASHRAE Guideline 36–compliant control (ASHRAE 36 control) and reinforcement learning (RL)–based control through experimental field tests and a simulation study. Three field tests were conducted at Oak Ridge National Laboratory’s commercial building test facility in Oak Ridge, Tennessee: a baseline with a baseline conventional control, a test with ASHRAE 36 control, and a test with RL-based control. The selected ASHRAE 36 controls were trim and respond control, as well as variable air volume (VAV) box control. We compared the measured supply air temperature of the rooftop unit, VAV box supply air temperature, and VAV box supply airflow rate across the three test cases. The field data indicated that ASHRAE 36 controls operated as specified by ASHRAE Guideline 36. Based on these data, ASHRAE 36 control achieved a 45 % reduction in hourly averaged HVAC energy consumption compared with the baseline, and RL-based control achieved a 66 % reduction. These potential annual energy savings were confirmed using a calibrated whole-building energy model. Compared with the baseline, ASHRAE 36 control reduced HVAC energy consumption by 42 %, and RL-based control achieved a 54 % reduction. Furthermore, RL-based control reduced total HVAC energy consumption by 21 % more than ASHRAE 36 control.
32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
Active multi-mode data analysis to improve fault diagnosis in AHUs
Faults in heating, ventilation and air conditioning systems can lead to increased energy consumption, occupant comfort issues, and reduced equipment lifetime. Commercial fault detection and diagnosis (FDD) tools has been increasingly deployed in U.S. commercial buildings. While they are helping to achieve energy efficiency and operational reliability, there remain gaps in their fault diagnostic capabilities. The diagnostic results often contain multiple distinct candidate root causes (CRCs) or offer no insight into CRCs. This study developed a novel active rule-based multi-mode data analysis method to enhance diagnostic resolution by applying proven rule sets and additional new rules to data from multiple known operational modes. The proposed method was demonstrated using enhanced air handling unit performance assessment rule sets and validated with the simulated data of two air handling units. New metrics, namely, reduced number of CRCs and improvement ratio, were developed to quantify the improvement of fault diagnostic resolution. The validation results showed that the proposed method effectively reduced the number of CRCs in contrast to analyzing data solely for a single mode of operation. It achieved a median improvement ratio of 80% in 19 test cases.
32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗