Large-eddy simulation of coal combustion in 15 MW pilot-scale boiler-simulation facility
Not provided.
Engineering↗
SEARCH · Search NASA
Results for “Simulation',”
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.
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↗
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↗
DISCOVR strain pipeline screening – Part II: Winter and summer season areal productivities and biomass compositional shifts in climate-simulation photobioreactor cultures
Assessing the seasonal biomass productivity and compositional shift dynamics under simulated outdoor culture conditions of the top 22 algae strains selected during Tier I flask screening is an important step in the further prioritization of strains with regard to outdoor pond cultivation. These top 22 strains were subjected to Tier II testing in the PNNL Laboratory Environmental Algae Pond Simulator (LEAPS) photobioreactors, simulating light and temperature conditions of 20 cm deep outdoor ponds during the Arizona winter and summer season. All strains were grown in two consecutive nutrient-replete batch culture experiments at their particular optimal medium salinity to quantify their respective seasonal linear-phase areal biomass productivities. To determine biomass compositional shifts in response to nutrient-depletion, the LEAPS cultures were allowed to enter a 9-day nutrient depletion phase at the end of the second batch run. The following strains were evaluated in winter-season climate-simulated cultures and are listed in the order from highest (7.9 g m -2 day -1 ) to lowest (2.3 g m -2 day -1 ) areal N-replete biomass productivity: Monoraphidium minutum 26B-AM, Tetraselmis striata LANL1001, Chlorella vulgaris LRB AZ-1201, Micractinium reisseri NREL14-F2, Monoraphidium sp. MONOR1, Chlorella vulgaris NREL4-C12, Scenedesmus obliquus UTEX393, Scenedesmus acutus LRB-AP-0401, Nannochloropsis oceanica CCAP849/10, and Stichococcus minutus CCALA727. The following strains were evaluated in summer-season climate-simulated cultures and are listed in the order from highest (31.8 g m -2 day -1 ) to lowest (14.2 g m -2 day -1 ) areal N-replete biomass productivity: Picochlorum renovo NREL39-A8, Scenedesmus obliquus UTEX393, Porphyridium cruentum CCMP675, Picochlorum celeri TG2-WT-CSM/EMRE, Chlorella sorokiniana DOE1116, Stichococcus minor CCMP819, Picochlorum oklahomensis CCMP2329, Chlorella sorokiniana DOE1412 (UTEXB3016), Scenedesmus rubescens NREL46B-D3, Picochlorum soloecismus DOE101, Tetraselmis striata LANL1001, Scenedesmus obliquus DOE 0152.z, and Agmenellum quadruplicatum UTEX2268. All LEAPS cultures experienced a significant reduction in areal biomass productivity in response to nutrient-depletion, from 7-16% in the winter season simulation and up to 1-60% in the summer season simulation. Finally, for 10 of the strains tested, the carbohydrate content more than doubled upon nutrient depletion, and for 9 strains, the lipid content increased by over 50% of the initial content.
09 BIOMASS FUELS↗
Suppressing simulation bias in multi-modal data using transfer learning
Abstract Many problems in science and engineering require making predictions based on few observations. To build a robust predictive model, these sparse data may need to be augmented with simulated data, especially when the design space is multi-dimensional. Simulations, however, often suffer from an inherent bias. Estimation of this bias may be poorly constrained not only because of data sparsity, but also because traditional predictive models fit only one type of observed outputs, such as scalars or images, instead of all available output data modalities, which might have been acquired and simulated at great cost. To break this limitation and open up the path for multi-modal calibration, we propose to combine a novel, transfer learning technique for suppressing the bias with recent developments in deep learning, which allow building predictive models with multi-modal outputs. First, we train an initial neural network model on simulated data to learn important correlations between different output modalities and between simulation inputs and outputs. Then, the model is partially retrained, or transfer learned, to fit the experiments; a method that has never been implemented in this type of architecture. Using fewer than 10 inertial confinement fusion experiments for training, transfer learning systematically improves the simulation predictions while a simple output calibration, which we design as a baseline, makes the predictions worse. We also offer extensive cross-validation with real and carefully designed synthetic data. The method described in this paper can be applied to a wide range of problems that require transferring knowledge from simulations to the domain of experiments.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
EMT-TS Hybrid Simulation for Large Power Grids Considering IBR-Driven Dynamics
The escalating integration of inverter-based resources (IBRs) poses new challenges to power systems by introducing fast dynamics with higher frequencies, which may need to be simulated by an electromagnetic transient (EMT) program. As an alternative to conducting EMT simulations for the entire system, which is typically time consuming, hybrid simulation between EMT and phasor-domain transient stability (TS) can greatly reduce the computational burden while preserving the detailed fast dynamics in the EMT zone. This paper establishes an EMT-TS hybrid simulation platform using open-source tools, specifically ParaEMT, GridPACK, and HELICS, which are the EMT simulator, TS simulator, and interface framework, respectively. Case studies on the 240-bus Western Electricity Coordinating Council (WECC) system demonstrate that the developed ParaEMT-HELICS-GridPACK hybrid simulator can accurately capture both slow electromechanical and fast IBR-driven dynamics with a 2.4x speedup.
electromagnetic transient↗
Comparison of equilibrium climate sensitivity estimates from slab ocean, 150-year, and longer simulations
We compare equilibrium climate sensitivity (ECS) estimates from pairs of long (= 800-year) control and abruptly quadrupled CO2 simulations with shorter (150, 300 year) coupled atmosphere-ocean simulations and Slab Ocean Models (SOM). Consistent with previous work, ECS estimates from shorter coupled simulations based on annual averages for years 1-150 underestimate those from SOM (-8% ± 13%) and long (-14% ± 8%) simulations. Analysis of only years 21-150 improved agreement with SOM (-2% ± 14%) and long (-8% ± 10%) estimates. Use of pentadal averages for years 51-150 results in improved agreement with long simulations (-4% ± 11%). While ECS estimates from current generation US models based on SOM and coupled annual averages of years 1-150 range from 2.6°C to 5.3°C, estimates based longer simulations of the same models range from 3.2°C to 7.0°C. Such variations between methods argues for caution in comparison and interpretation of ECS estimates across models.
Dunne, John P.↗
Hybrid RANS-LES of the Atmospheric Boundary Layer for Wind Farm Simulations: Preprint
Wind farm simulations often do not accurately represent wake-atmospheric boundary layer (ABL) interactions, blade boundary layer (BL) dynamics, and turbine-turbine interactions. In this work, we use Active Model Split (AMS), a new hybrid Reynolds-Averaged Navier Stokes (RANS)-large eddy simulation (LES) model, which is well suited to capture these effects because the model can (i) accurately simulate the ABL with the Coriolis effect, (ii) is accurate in adverse pressure gradients such as those near wind turbine blades, and (iii) has sufficiently low computational cost to simulate multiple turbines while resolving the blade BL. For simplicity and consistency we develop AMS to be used throughout the domain rather than in a zonal method. We implement our work in the massively parallel flow solver, Nalu-Wind, so that our model can access the compute resources needed for blade-resolved simulations of multiple wind turbines. To accomplish these aims, we modify the baseline AMS by changing the RANS contribution to SST k - omega with a length scale limiter, adding the Coriolis effect, and developing an appropriate wall treatment. We show that AMS of the ABL with the Coriolis effect matches LES reference results better than those obtained with RANS. We describe our plans to add buoyancy effects and wind turbines to our AMS simulations.
atmospheric boundary layer↗