Search NASA⌕ Search

SEARCH · Search NASA

Results for “simulators”

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 19 records

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)↗

Quantum simulations of nuclei and nuclear pasta with the multiresolution adaptive numerical environment for scientific simulations

Neutron star and supernova matter at densities just below the nuclear matter saturation density is expected to form a lattice of exotic shapes. These so-called nuclear pasta phases are caused by Coulomb frustration. Their elastic and transport properties are believed to play an important role for thermal and magnetic field evolution, rotation, and oscillation of neutron stars. Furthermore, they can impact neutrino opacities in core-collapse supernovae. In this work, we present proof-of-principle three-dimensional (3D) Skyrme Hartree-Fock (SHF) simulations of nuclear pasta with the Multi-resolution ADaptive Numerical Environment for Scientific Simulations (MADNESS). We perform benchmark studies of 16 O, 208 Pb, and 238 U nuclear ground states and calculate binding energies via 3D SHF simulations. Results are compared with experimentally measured binding energies as well as with theoretically predicted values from an established SHF code. The nuclear pasta simulation is initialized in the so-called waffle geometry as obtained by the Indiana University Molecular Dynamics (IUMD) code. The size of the unit cell is 24 fm with an average density of about ρ = 0.05fm -3 , proton fraction of Yp = 0.3, and temperature of T = 0 MeV. Overall, our calculations reproduce the binding energies and shapes of light and heavy nuclei with different geometries. For the pasta simulation, we find that the final geometry is very similar to the initial waffle state. We compare calculations with and without spin-orbit forces. We find that while subtle differences are present, the pasta phase remains in the waffle geometry. Within the MADNESS framework, we can successfully perform calculations of inhomogeneous nuclear matter. By using pasta configurations from IUMD it is possible to explore different geometries and test the impact of self-consistent calculations on the latter.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Modeling and simulation challenges pursued by the Consortium for Advanced Simulation of Light Water Reactors (CASL)

The Consortium for the Advanced Simulation of Light Water Reactors (CASL), the first Energy Innovation Hub of the Department of Energy, was established in 2010 with the goal of providing modeling and simulation (M&S) capabilities that support and accelerate the improvement of nuclear energy's economic competitiveness and the reduction of spent nuclear fuel volume per unit energy, and all while assuring nuclear safety. To accomplish this requires advances in M&S capabilities in radiation transport, thermal-hydraulics, fuel performance and corrosion chemistry. To focus CASL's R&D, industry challenge problems have been defined, which equate with long standing issues of the nuclear power industry that M&S can assist in addressing. To date CASL has developed a multi-physics “core simulator” based upon pin-resolved radiation transport and subchannel (within fuel assembly) thermal-hydraulics, capitalizing on the capabilities of high performance computing. CASL's fuel performance M&S capability can also be optionally integrated into the core simulator, yielding a coupled multi-physics capability with untapped predictive potential. Material models have been developed to enhance predictive capabilities of fuel clad creep and growth, along with deeper understanding of zirconium alloy clad oxidation and hydrogen pickup. Understanding of corrosion chemistry (e.g., CRUD formation) has evolved at all scales: micro, meso and macro. CFD R&D has focused on improvement in closure models for subcooled boiling and bubbly flow, and the formulation of robust numerical solution algorithms. For multiphysics integration, several iterative acceleration methods have been assessed, illuminating areas where further research is needed. Finally, uncertainty quantification and data assimilation techniques, based upon sampling approaches, have been made more feasible for practicing nuclear engineers via R&D on dimensional reduction and biased sampling. Industry adoption of CASL's evolving M&S capabilities, which is in progress, will assist in addressing long-standing and future operational and safety challenges of the nuclear industry. - Highlights: • Complexity of physics based modeling of light water reactor cores being addressed. • Capability developed to help address problems that have challenged the nuclear power industry. • Simulation capabilities that take advantage of high performance computing developed.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

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↗

Hierarchical calibration and validation framework of bench-scale computational fluid dynamics simulations for solvent-based carbon capture. Part 2: Chemical absorption across a wetted wall column: Original Research Article: Hierarchical calibration and validation framework of bench-scale computational fluid dynamics simulations

Part 1 of this paper presents a numerical model for non-reactive physical mass transfer across a wetted wall column (WWC). In Part 2, we improved the existing computational fluid dynamics (CFD) model to simulate chemical absorption occurring in a WWC as a bench-scale study of solvent-based carbon dioxide (CO 2 ) capture. In this study, to generate data for WWC model validation, CO 2 mass transfer across a monoethanolamine (MEA) solvent was first measured on a WWC experimental apparatus. The numerical model developed in this work can account for both chemical absorption and desorption of CO 2 in MEA. In addition, the overall mass transfer coefficient predicted using traditional/empirical correlations is conducted and compared with CFD prediction results for both steady and wavy falling films. A Bayesian statistical calibration algorithm is adopted to calibrate the reaction rate constants in chemical absorption/desorption of CO 2 across a falling film of MEA. The posterior distributions of the two transport properties, i.e., Henry's constant and gas diffusivity in the non-reacting nitrous oxide (N 2 O)/MEA system obtained from Part 1 of this study, serves as priors for the calibration of CO 2 reaction rate constants after using the N 2 O/CO 2 analogy method. Finally, the calibrated model can be used to predict the CO 2 mass transfer in a WWC for a wider range of operating conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

Conversion and Validation of Distribution System Model from a QSTS-Based Tool to a Real-Time Dynamic Phasor Simulator

A software process is developed to convert distribution network models from a quasi-static time-series tool (OpenDSS) to a real-time dynamic phasor simulator (ePHASORSIM). The description of this process in this paper would be helpful for researchers who intend to perform similar conversions. The converter could be utilized directly by users of real-time simulators who intend to perform software-in-the-loop or hardware-in-the-loop tests on large distribution test feeders for a range of use cases, including testing functions of advanced distribution management systems against a simulated distribution system. In the future, the developers intend to release the conversion tool as open source to enable use by others.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Two-temperature model in molecular dynamics simulations of cascades in Ni-based alloys

In high-energy irradiation events, energy from the fast moving ion is transferred to the system via nuclear and electronic energy loss mechanisms. The nuclear energy loss results in the creation of point defects and clusters, while the energy transferred to the electrons results in the creation of high electronic temperatures, which can affect the damage evolution. We perform molecular dynamics simulations of 30 keV and 50 keV Ni ion cascades in nickel-based alloys without and with the electronic effects taken into account. We compare the results of classical molecular dynamics (MD) simulations, where the electronic effects are ignored, with results from simulations that include the electronic stopping only, as well as simulations where both the electronic stopping and the electron-phonon coupling are incorporated, as described by the two temperature model (2T-MD). Our results indicate that the 2T-MD leads to a smaller amount of damage, more isolated defects and smaller defect clusters.

Zarkadoula, Eva↗

The structure of liquid water up to 360 MPa from x-ray diffraction measurements using a high Q-range and from molecular simulation

X-ray diffraction measurements of liquid water are reported at pressures up to 360 MPa corresponding to a density of 0.0373 molecules per Å{sup 3}. The measurements were conducted at a spatial resolution corresponding to Q{sub max} = 16 Å{sup −1}. The method of data analysis and measurement in this study follows the earlier benchmark results reported for water under ambient conditions having a density of 0.0333 molecules per Å{sup 3} and Q{sub max} = 20 Å{sup −1} [J. Chem. Phys. 138, 074506 (2013)] and at 70 °C having a density of 0.0327 molecules per Å{sup 3} and Q{sub max} = 20 Å{sup −1} [J. Chem. Phys. 141, 214507 (2014)]. The structure of water is very different at these three different T and P state points and thus they provide the basis for evaluating the fidelity of molecular simulation. Measurements show that at 360 MPa, the 4 waters residing in the region between 2.3 and 3 Å are nearly unchanged: the peak position, shape, and coordination number are nearly identical to their values under ambient conditions. However, in the region above 3 Å, large structural changes occur with the collapse of the well-defined 2nd shell and shifting of higher shells to shorter distances. The measured structure is compared to simulated structure using intermolecular potentials described by both first-principles methods (revPBE-D3) and classical potentials (TIP4P/2005, MB-pol, and mW). The DFT-based, revPBE-D3, method and the many-body empirical potential model, MB-pol, provide the best overall representation of the ambient, high-temperature, and high-pressure data. The revPBE-D3, MB-pol, and the TIP4P/2005 models capture the densification mechanism, whereby the non-bonded 5th nearest neighbor molecule, which partially encroaches the 1st shell at ambient pressure, is pushed further into the local tetrahedral arrangement at higher pressures by the more distant molecules filling the void space in the network between the 1st and 2nd shells.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modulating factors of hydrologic exchanges in a large-scale river reach: Insights from three-dimensional computational fluid dynamics simulations

Hyporheic exchange is a critical mechanism shaping hydrological and biogeochemical processes along a river corridor. Recent studies on quantifying the hyporheic exchange were mostly limited to local scales due to field inaccessibility, computational demand, and complexity of geomorphology and subsurface geology. Surface flow conditions and subsurface physical properties are well known factors on modulating the hyporheic exchange, but quantitative understanding of their impacts on the strength and direction of hyporheic exchanges at reach scales is absent. In this study, a high resolution computational fluid dynamics (CFD) model that couples surface and subsurface flow and transport is employed to simulate hyporheic exchanges in a 7-km long reach along the main-stem of the Columbia River. Assuming that the hyporheic exchange does not affect surface water flow conditions due to its negligible magnitude compared to the volume and velocity of river water, we developed a one-way coupled surface and subsurface water flow model using the commercial CFD software STAR-CCM+. The model integrates the Reynolds-averaged Navier-Stokes (RANS) equation solver with a realizable k-e two-layer turbulence model, a two-layer all y^+ wall treatment, and the volume of fluid (VOF) method, and is used to simulate hyporheic exchanges by tracking the free water-air interface as well as flow in the river and the subsurface porous media. The model is validated against measurements from acoustic Doppler current profiler (ADCP) in the stream water and hyporheic fluxes derived from a set of temperature profilers installed across the riverbed. The validated model is then employed to systematically investigate how hyporheic exchanges are influenced by surface water fluid dynamics strongly regulated by upstream dam operations, as well as subsurface structures (e.g. thickness of riverbed and subsurface formation layers) and hydrogeological properties (e.g. permeability). The results suggest that the thickness of riverbed alluvium layer is the dominant factor for reach-scale hyporheic exchanges, followed by the alluvium permeability, the depth of the underlying impermeable layer, and the assumption of hydrostatic pressure.

Bao, Jie↗

High Level Architecture Distributed Space System Simulation for Simulation Interoperability Standards Organization Simulation Smackdown

Modeling and Simulation plays a very important role in mission design. It not only reduces design cost, but also prepares astronauts for their mission tasks. The SISO Smackdown is a simulation event that facilitates modeling and simulation in academia. The scenario of this year s Smackdown was to simulate a lunar base supply mission. The mission objective was to transfer Earth supply cargo to a lunar base supply depot and retrieve He-3 to take back to Earth. Federates for this scenario include the environment federate, Earth-Moon transfer vehicle, lunar shuttle, lunar rover, supply depot, mobile ISRU plant, exploratory hopper, and communication satellite. These federates were built by teams from all around the world, including teams from MIT, JSC, University of Alabama in Huntsville, University of Bordeaux from France, and University of Genoa from Italy. This paper focuses on the lunar shuttle federate, which was programmed by the USRP intern team from NASA JSC. The shuttle was responsible for provide transportation between lunar orbit and the lunar surface. The lunar shuttle federate was built using the NASA standard simulation package called Trick, and it was extended with HLA functions using TrickHLA. HLA functions of the lunar shuttle federate include sending and receiving interaction, publishing and subscribing attributes, and packing and unpacking fixed record data. The dynamics model of the lunar shuttle was modeled with three degrees of freedom, and the state propagation was obeying the law of two body dynamics. The descending trajectory of the lunar shuttle was designed by first defining a unique descending orbit in 2D space, and then defining a unique orbit in 3D space with the assumption of a non-rotating moon. Finally this assumption was taken away to define the initial position of the lunar shuttle so that it will start descending a second after it joins the execution. VPN software from SonicWall was used to connect federates with RTI during testing and the Smackdown event. HLA software from Pitch Technology and MAK Technology were used to edit and extend FOM and provide HLA services for federation execution. The SISO Smackdown event for 2011 was held in Boston, Massachusetts. The federation execution lasted for one hour, and the event was very successful in catching the attention of university students and faculties.

Li, Zuqun↗