Large-eddy simulation of coal combustion in 15 MW pilot-scale boiler-simulation facility
Not provided.
Engineering↗
SEARCH · Search NASA
Results for “simulations”
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.
Improving Land Surface Temperature Simulation in CoLM Over the Tibetan Plateau through Fractional Vegetation Cover Derived from a Remotely Sensed Clumping Index and Model-Simulated Leaf Area Index
Parameterizations of fractional vegetation cover (FVC) in land surface models have important effects on simulations of surface energy budget, especially in arid and semiarid regions. This study uses a FVC scheme in which FVC is derived from leaf area index and a remotely sensed clumping index. The performance of the new scheme (SMFVC) is evaluated against Moderate Resolution Imaging Spectroradiometer (MODIS) land surface temperature (LST) and in situ soil temperature observations, together with two other FVC schemes, a general FVC scheme (CTL) based on land cover map and a climatology-based FVC scheme (RSFVC) that uses long-term remotely sensed Normalized Difference Vegetation Index of MODIS. The three FVC schemes were implemented in the Common Land Model (CoLM) and applied in the Tibetan Plateau using the same forcing data and default parameters. Our results demonstrate that FVC schemes have significant influence on the CoLM performance: (1) the RSFVC and SMFVC schemes significantly reduce the LST biases found in CTL, particularly in grassland and during summer; (2) soil temperature evaluation by in situ observations from three networks on the Tibetan Plateau corroborates the LST results; and (3) the improvements are mainly related to representing temporal (seasonal) variability and subgrid heterogeneity of FVC, which improves surface albedo and surface energy balance. In other words, by including more vegetation characteristics, such as using a clumping index, land surface models may better simulate surface vegetation condition and further better represent the land surface energy budget over the Tibetan Plateau.
Li, Chengwei↗
Dark Energy Survey Year 3 results: Simulation-based cosmological inference with wavelet harmonics, scattering transforms, and moments of weak lensing mass maps. Validation on simulations
Beyond-two-point statistics contain additional information on cosmological as well as astrophysical and observational (systematics) parameters. In this methodology paper we provide an end-to-end simulationbased analysis of a set of Gaussian and non-Gaussian weak lensing statistics using detailed mock catalogs of the Dark Energy Survey (DES). Here, we implement: 1) second and third moments; 2) wavelet phase harmonics (WPH); 3) the scattering transform (ST). Our analysis is fully based on simulations, it spans a space of seven $νw$CDM cosmological parameters, and it forward models the most relevant sources of systematics of the data (masks, noise variations, clustering of the sources, intrinsic alignments, and shear and redshift calibration). We implement a neural network compression of the summary statistics, and we estimate the parameter posteriors using a likelihood-free-inference approach. We validate the pipeline extensively, and we find that WPH exhibits the strongest performance when combined with second moments, followed by ST, and then by third moments. The combination of all the different statistics further enhances constraints with respect to second moments, up to 25 percent, 15 percent, and 90 percent for S 8 , Ω m , and the figure-of-merit FoM S8;Ωm , respectively. We further find that non-Gaussian statistics improve constraints on w and on the amplitude of intrinsic alignment with respect to second moments constraints. The methodological advances presented here are suitable for application to Stage IV surveys from Euclid, Rubin-LSST, and Roman with additional validation on mock catalogs for each survey. In a companion paper we present an application to DES Year 3 data.
79 ASTRONOMY AND ASTROPHYSICS↗
Accelerating Simulation for High-Fidelity PV Inverter System Reliability Assessment with High-Performance Computing
The overall cost of photovoltaic (PV) systems has shown a downward trend during the last decade; however, PV inverter failures account for the highest cost of operation and maintenance. To address this, reliability tools with powerful computation and better accuracy are required for the lifetime prediction and degradation evaluation of PV inverters. This paper proposes an event-driven parallel computing-based simulator. The proposed simulator applies high-performance computing techniques and other accessory optimization techniques-including cluster merging, adaptive model updates, and steady-state identification-to make reliability assessments for PV inverters under given input mission profiles and operating conditions with high efficiency and high fidelity. The main idea of the simulator and its workflow are introduced. Then, a demo PV inverter system simulator is implemented, and the speedup of the total simulations of the switching model reaches 123.03 times.
high-performance computing↗
Tropical Tree Crop Simulation with a Process-Based, Daily Timestep Simulation Model (ALMANAC): Description of Model Adaptation and Examples with Coffee and Cocoa Simulations
Coffee (Coffea species) and Cocoa (Theobroma cacao) are important cash crops grown in the tropics but traded globally. This study was conducted to apply the ALMANAC model to these crops for the first time, and to test its ability to simulate them under agroforestry management schemes and varying precipitation amounts. To create this simulation, coffee was grown on a site in Kaua’i, Hawai’i, USA, and cocoa was grown on a site in Sefwi Bekwai, Ghana. A stand-in for a tropical overstory tree was created for agroforestry simulations using altered parameters for carob, a common taller tropical tree for these regions. For both crops, ALMANAC was able to realistically simulate yields when compared to the collected total yield data. On Kaua’i, the mean simulated yield was 2% different from the mean measured yield, and in all three years, the simulated values were within 10% of the measured values. For cocoa, the mean simulated yield was 3% different from the mean measured yield and the simulated yield was within 10% of measured yields for all four available years. When precipitation patterns were altered, in Ghana, the wetter site showed lower percent changes in yield than the drier site in Hawai’i. When agroforestry-style management was simulated, a low Leaf Area Index (LAI) of the overstory showed positive or no effect on yields, but when LAI climbed too high, the simulation was able to show the detrimental effect this competition had on crop yields. These simulation results are supported by other literature documenting the effects of agroforestry on tropical crops. This research has applied ALMANAC to new crops and demonstrated its simulation of different management and environmental conditions. The results show promise for ALMANAC’s applicability to these scenarios as well as its potential to be further tested and utilized in new circumstances.
60 APPLIED LIFE SCIENCES↗
Comparison of Simulation of Dual-Active Bridge in Different Simulators and Using Different Simulation Methods
The accuracy of states measured from simulation of a dual-active bridge (DAB) converter with an open-loop controller has been used for comparing simulation software and algorithms. The accuracy of the measured states is sensitive to simulation settings and parameters chosen. In this paper, the simulation of the DAB converter with open-loop control is run in the MATLAB/Simscape Electrical, PLECS, and PSCAD software environments, with various simulation settings. The simulation results are also compared with a new simulation algorithm proposed for DAB simulations. Different simulation settings are evaluated, including different solver types, operating conditions, and the on-state resistance values of the MOSFET switches. The simulation results are compared with the theoretical calculations. It is shown that unexpected divergence and DC offset may occur in some cases. This paper provides guidelines to simulate the DAB converter with open-loop control.
Yin, Shiyuan↗
COMPARISON OF FIELD MEASUREMENTS AND LARGE EDDY SIMULATIONS OF THE SCALED WIND FARM TECHNOLOGY (SWIFT) SITE.
Abstract not provided.
Blaylock, Myra L.↗
Continuum stress intensity factors from atomistic fracture simulations
Stress intensity factors (SIFs) are utilized in continuum fracture mechanics to quantify the stress fields surrounding a crack in a homogeneous material in the linear elastic regime. Critical values of the SIFs define an intrinsic measure of the resistance of a material to propagate a crack. At atomic scales, yet, fracture occurs as a series of atomic bonds breaking, differing from the continuum description. As a consequence, a formal analog of the continuum SIFs calculated from atomistic simulations can have spatially localized, microstructural contributions that originate from varying bond configurations. The ability to characterize fracture at the atomic scale in terms of the SIFs offers both an opportunity to probe the effects of chemistry, as well as how the addition of a microstructural component affects the accuracy. We present a novel numerical method to determine SIFs from molecular dynamics (MD) simulations. The accuracy of this approach is first examined for a simple model, and then applied to atomistic simulations of fracture in amorphous silica. MD simulations provide time and spatially dependent SIFs, with results that are shown to be in good agreement with experimental values for fracture toughness in silica glass.
42 ENGINEERING↗
Dehalogenation of electrochemical processing salt simulants with ammonium phosphates and immobilization of salt cations in an iron phosphate glass waste form
This study demonstrates the conversion of simulated pyroprocessing salt wastes [i.e., ER(SF) and ERV2] into chemically durable iron phosphate glass waste forms through reactions with ammonium hydrogen phosphate precursors [i.e., NH4H2PO4 and (NH4)2HPO4] and added Fe2O3. The reaction of the salt simulants with the phosphate precursors resulted in dehalogenation of the salt through the production of NH4Cl and NH4I (along with water), of which the NH4Cl can be used to produce UCl3 oxidant for use in the electrorefiner during pyroprocessing. These iron phosphate glasses contain a notably higher waste loading than the baseline glass-bonded sodalite ceramic waste form, which requires about 70 vol% sodalite to accommodate Cl from the salt, because they are dehalogenated prior to vitrification. The high waste loadings coupled with higher bulk densities (~3×103 kg/m3) than glass-bonded sodalite (~2.3×103 kg/m3) mean that these glasses have a higher storage volume, i.e., the volume required to immobilize a given mass of waste. The iron phosphate glass waste forms were characterized using a variety of techniques including physical property measurements and chemical durability testing. Following initial testing, a proof-of-concept demonstration was conducted to show that the NH4Cl generated from this process could be used to produce UCl3 using depleted uranium dendrites as a feedstock.
phosphate glass, pyroprocessing, waste forms, salt↗
Large-eddy simulations of the Northeastern US coastal marine boundary layer
In this study, large eddy simulations (LES) of offshore boundary layers near the Nantucket coast are performed using Nalu-Wind. The marine boundary layer conditions are chosen to match the predominantly unstable and neutral conditions measured by the Cape Wind platform. The appropriate domain, resolution, and boundary condition settings required for the LES are established through this work. Differences between stable and unstable cases are found in the wind speed profiles, averaged statistics, and wind spectra, and explained in terms of stratification effects.
54 ENVIRONMENTAL SCIENCES↗
Investigation of Farm-to-Farm Interactions and Blockage Effects from AWAKEN using Large-Scale Numerical Simulations
A large-scale numerical computation of five wind farms was performed as a part of the American Wake Experiment (AWAKEN). This high-fidelity computation used the ExaWind/AMR-Wind LES solver to simulate a 100km × 100km domain containing 541 turbines under unstable atmospheric conditions matching previous measurements. The turbines were represented by Joukowski and OpenFAST coupled actuator disk models. Results of this qualitative comparison illustrate the interactions with wind farms with large scale ABL structures in the flow, as well as the extent of downstream wake penetration in the flow and blockage effects around wind farms.
Cheung, Lawrence C.↗
Fast Quasi-Static Time-Series Simulation for Accurate PV Inverter Semiconductor Fatigue Analysis with a Long-Term Solar Profile
Power system simulations with long-term data typically have large time steps varying from one second to a few minutes. However, for PV inverter semiconductors, the minimum thermal stress cycle occurs over the fundamental grid frequency (50 or 60 Hz). This requires the time step of the fatigue simulation to be around 100 µs. This small time step requires long computation times to process yearly power production profiles. This paper proposes a fast fatigue simulation for inverter semiconductors using the quasi-static time series (QSTS) simulation concept. The proposed simulation calculates the steady state of the semiconductor junction temperature by using a Fast Fourier Transform (FFT). The small thermal cycling during a switching period and even over the fundamental waveform is disregarded to further accelerate the simulation speed. The resulting time step of the fatigue simulation is 15 minutes, which is consistent with the solar dataset. The error of the proposed simulation is 0.16% compared to the fatigue simulation results using the complete thermal stress profile. A PV inverter that responds to a Transactive Energy System (TES) is simulated to demonstrate the use of the proposed fatigue simulation. The proposed simulation has the potential to co-simulate with system level simulation tools that also adopt the QSTS concept.
Liu, Yunting↗
Reactive Molecular Dynamics Simulations and Quantum Chemistry Calculations To Investigate Soot-Relevant Reaction Pathways for Hexylamine Isomers
Sooting tendencies of a series of nitrogen-containing hydrocarbons (NHCs) have been recently characterized experimentally using the yield sooting index (YSI) methodology. This work aims to identify soot-relevant reaction pathways for three selected C6H15N amines, namely, dipropylamine (DPA), diisopropylamine (DIPA), and 3,3-dimethylbutylamine (DMBA) using ReaxFF molecular dynamics (MD) simulations and quantum mechanical (QM) calculations and to interpret the experimentally observed trends. ReaxFF MD simulations are performed to determine the important intermediate species and radicals involved in the fuel decomposition and soot formation processes. QM calculations are employed to extensively search for chemical reactions involving these species and radicals based on the ReaxFF MD results and also to quantitatively characterize the potential energy surfaces. Specifically, ReaxFF simulations are carried out in the NVT ensemble at 1400, 1600, and 1800 K, where soot has been identified to form in the YSI experiment. These simulations account for the interactions among test fuel molecules and pre-existing radicals and intermediate species generated from rich methane combustion, using a recently proposed simulation framework. ReaxFF simulations predict that the reactivity of the amines decrease in the order DIPA > DPA > DMBA, independent of temperature. Both QM calculations and ReaxFF simulations predict that C2H4, C3H6, and C4H8 are the main nonaromatic soot precursors formed during the decomposition of DPA, DIPA, and DMBA, respectively, and the associated reaction pathways are identified for each amine. Both theoretical methods predict that sooting tendency increases in the order DPA, DIPA, and DMBA, consistent with the experimentally measured trend in YSI. This work demonstrates that sooting tendencies and soot-relevant reaction pathways of fuels with unknown chemical kinetics can be identified efficiently through combined ReaxFF and QM simulations. Overall, predictions from ReaxFF simulations and QM calculations are consistent, in terms of fuel reactivity, major intermediates, and major nonaromatic soot precursors.
fuel decomposition↗
Advanced Modeling and Simulation Methods for Evaluation of Thermal Neutron Scattering Materials
With the rise of interest in thermal neutron scattering data for advanced reactor, criticality safety, and shielding applications, new experimental data are required for evaluation of new materials or for re-evaluation (or validations) of previously evaluated materials. New experimental data are evaluated in a three-step process: (1) computing the phonon characteristics, (2) computing the dynamic structure factor (DSF) from the data, and (3) using the experimental setup to simulate the experimental data. All three steps have challenges, ranging from the need for a sufficiently general material simulation code—a processing code that can compute the corresponding DSF—to having a detailed layout of the instrument/beamline/facility where the data were measured. Whereas phonon characteristics of materials can be computed using various methods (molecular dynamics, density functional theory, etc.), a high-fidelity computation of the DSF and the simulation of the experiment based on the DSF is vital to the accuracy of the evaluation. The latter two steps can be achieved by using the two corresponding code systems developed by instrument scientists at the Spallation Neutron Source (SNS) at Oak Ridge National Laboratory: (1) OCLIMAX, a program that calculates the dynamic structure factor from DFT and MD simulation results, and (2) MCViNE, a Monte Carlo neutron ray-tracing program designed to simulate neutron scattering experiments. Recently, polyethylene and yttrium hydride were measured at the Wide Angular-Range Chopper (ARCS) and SEQUOIA instrument stations of the SNS. These experiments are simulated using the density functional theory code, the Cambridge Serial Total Energy Package (CASTEP), to compute its phonon characteristics (eigenvalues/vectors and PDOS), which is then processed using OCLIMAX to yield the DSF, and finally the data at each instrument station are simulated by the MCViNE for comparison to the measured data for evaluation. For comparison to conventional evaluation methods, the scattering data processed from OCLIMAX are compared against those processed from the LEAPR module of NJOY, and the results from MCViNE simulations are compared against previously used simplified beamline models implemented in the Monte Carlo N-Particle (MCNP) code.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Accuracy Enhancement of Nuclear Power Plant Simulators Utilizing High Accuracy Simulation Predictions
More recently, reactor core simulators for core designs associated with commercial nuclear power plants that utilize what is believed to be higher fidelity models have been developed. Features such as neutronics models that utilize transport equation solvers with fine spatial meshes and many energy-groups, thermal-hydraulic models that utilize sub-channel solvers with fine spatial mesh and capable of treating a wide range of fluid conditions, and fuel-coolant chemistry interaction models capable of treating CRUD deposition are to be found in these higher fidelity core simulators. These reactor core simulators require access to higher performance computers, characterized by many processors, cores and large memory. So associated with utilization of these simulators is access to high performance computers and ability to accommodate in one’s workflow longer execution times. By contrast, currently used core simulators by the nuclear industry can execute on engineering workstations and have execution times of seconds to minutes. The desirability for having short execution times is not only desired for support of time critical tasks but supports the mental process of decision making by engineers. The goal of the work reported upon here has the objective of retaining the fidelity of higher fidelity models while retaining the ability to utilize engineering workstations. Beyond the core simulator goal, additional goals of this work include incorporating the just described core simulator capability into a Nuclear Steam Supply System (NSSS) simulator, and to incorporate the resulting capability into an environment supportive of design and operational decision making associated with nuclear power stations. The model selected for the core neutronics model is the NESTLE code, for the core thermal-hydraulic model is the CTF code utilizing coarse mesh, and for the NSSS model is the RELAP5-3D code. WSC’s proprietary 3KEYMASTERTM platform is being used to provide software coupling, user interface, visualization, and reporting. The NESTLE core neutronics simulator was first integrated with the CTF core thermal-hydraulic simulator using CTF developed communication commands which are also used for CTF to communicate with RELAP5-3D under WSC’s proprietary 3KEYMASTERTM platform. To assure NESTLE prediction consistency with higher fidelity core neutronic simulators, buffer codes have been created to automatically generate from output files written by the VERA core simulator the NESTLE nodal neutronic parameter’ library, geometry, and pin-power reconstruction input files, thereby avoiding a number of challenges associated with utilizing lattice physics codes and providing consistency with VERA predictions. To treat absorber rod effects a multi-set library is utilized, where a set refers to a specific absorber rod fully inserted pattern. A coarse spatial mesh CTF model was developed with features added that support using CTF as envisioned in the engineering quality simulator. A hybrid meshing approach was implemented to allow for automated construction of models with mixed levels of refinement. Specifically, a core model could resolve some assemblies at a nodal level (4 subchannels per assembly) and others at a pin-resolution (one subchannel per coolant subchannel in the assembly). The intention is that this will allow for better resolution of limiting conditions such as DNBR and PCT, which are based on local rod and subchannel conditions. Further development was done of features that enhance the capabilities for the envisioned engineering quality simulator that has been developed, but now for RELAP-3D. The RELAP5-3D code development includes ability to model more than 999 components and the addition of the cross-channels turbulence mixing model and the void drift model that are implemented in CTF, aiming to achieve closer prediction agreement of the two codes for transient simulations, specifically, more accurate matches of the overall mass, momentum, and energy exchanges of both the liquid and gas phases between the neighboring core assemblies. Graphics were also developed for the Instructor Station for this project under WSC’s proprietary 3KEYMASTERTM platform to facilitate design and operational decision making.
42 ENGINEERING↗
Simulation of Rough Electrodes Coupled with Structural Dynamics
Explore the source record for details and available documents.
Dankesreiter, Benjamin Thomas↗
Finite Element Method Simulation Portion of eXtremeMAT:Modeling and Simulation Thrust-2: Surrogate Modeling Approach and Finite Element Simulations
Brief overview of the highlights of work done on the integration of a surrogate reduced order model into the FEM solver code MOOSE under the eXtremeMAT: Modeling and Simulation project
36 MATERIALS SCIENCE↗
Advanced EMT and Phasor-Domain Hybrid Simulation With Simulation Mode Switching Capability for Transmission and Distribution Systems
Conventional electromagnetic transient (EMT) and phasor do-main hybrid simulation approaches presently exist for transmission system level studies and their simulation efficiency is generally constrained by the EMT simulation. With an increasing number of distributed energy resources and non-conventional loads being installed in the distribution systems, it is imperative to ex-tend the hybrid simulation application to include distribution systems and integrated transmission and distribution systems. Meanwhile, it is equally important to improve the simulation efficiency as the modeling scope and complexity of the detailed system in the EMT simulation increases. To meet both requirements, this paper introduces an advanced electromagnetic transient (EMT) and phasor domain hybrid simulation approach. This approach has two outstanding features: 1) a comprehensive phasor domain modeling framework which supports positive-sequence, three-sequence, three-phase and mixed three-sequence/three-phase representations; 2) a robust and flexible simulation mode switching scheme which enables simulation switching from hybrid simulation mode back to pure phasor domain dynamic simulation mode to achieve significantly improved simulation efficiency. The proposed method has been tested on integrated transmission and distribution systems, and the results show that with the developed simulation switching feature, the total computational time is significantly reduced compared to running the hybrid simulation for the whole simulation period, while good accuracy is maintained.
Huang, Qiuhua (ORCID:0000000184577429)↗