Search NASA⌕ Search

SEARCH · Search NASA

Results for “Simulation System”

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 1,369 records · Page 76

Grating Polarizers at 170 GHz for ECRH Systems: Low Power Tests and Simulations

The elliptical polarization produced by a pair of sinusoidally profiled miter-bend grating polarizers has been measured experimentally and found to be in excellent agreement with theory. The polarizers were designed to be used at 170 GHz in the miter bends of a 63.5 mm corrugated metallic waveguide for creating the arbitrary elliptically polarized microwave beam needed for electron cyclotron resonance heating of plasma. Using a vector network analyzer to generate a linearly polarized HE 11 incident mode, each polarizer was individually tested to measure the rotation α and ellipticity β of the elliptically polarized reflected microwave beam as a function of the grating rotation angle. The grating polarizers were then tested together to measure the elliptical polarization as a function of the angles of the combination of the two polarizers arranged sequentially on the transmission line. The map of the ellipticity versus grating angles agreed very well with numerical simulations using high-frequency structure simulator. Numerical simulations show that up to eight combinations of rotation angles of the two grating polarizers can provide the same ellipticity but with the varying ohmic loss, thus allowing the choice of settings to minimize the ohmic loss, a necessary consideration to minimize heating effects in megawatt power level systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Grating Polarizers at 170 GHz for ECRH Systems: Low Power Tests and Simulations

The elliptical polarization produced by a pair of sinusoidally profiled miter-bend grating polarizers has been measured experimentally and found to be in excellent agreement with theory. The polarizers were designed to be used at 170 GHz in the miter bends of a 63.5 mm corrugated metallic waveguide for creating the arbitrary elliptically polarized microwave beam needed for electron cyclotron resonance heating of plasma. Using a vector network analyzer to generate a linearly polarized HE11 incident mode, each polarizer was individually tested to measure the rotation α and ellipticity β of the elliptically polarized reflected microwave beam as a function of the grating rotation angle. The grating polarizers were then tested together to measure the elliptical polarization as a function of the angles of the combination of the two polarizers arranged sequentially on the transmission line. The map of the ellipticity versus grating angles agreed very well with numerical simulations using high-frequency structure simulator. Numerical simulations show that up to eight combinations of rotation angles of the two grating polarizers can provide the same ellipticity but with the varying ohmic loss, thus allowing the choice of settings to minimize the ohmic loss, a necessary consideration to minimize heating effects in megawatt power level systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NREL Innovations Help Drive Wind Industry Transformation [Fact Sheet]

For nearly 30 years, NREL has helped the wind turbine industry through design and research innovations. The comprehensive capabilities of the National Wind Technology Center (NWTC), ranging from specialized computer simulation tools to unique test facilities, has been used to design, develop, and deploy several generations of advanced wind energy technology.

17 WIND ENERGY↗

Robust Restoration From Cyber-Physical Attacks in Active Distribution Grids With Grid-Edge IBRs

The inverter-based resources (IBRs) have enabled the integration of renewable energy at the grid edge with enhanced control capabilities to support the reliable operation of power grids. Different control frameworks, such as hierarchical or distributed architecture, have been proposed with the expansion of cyber networks for real-time monitoring and control. This evolution of critical infrastructure into cyber-physical systems also brings more vulnerabilities for the broadened attack surfaces, and significantly increases the possibility of physical system failures or outages caused by cyberattacks. Among tremendous efforts in the defense-in-depth approach, it remains challenging to provide prompt detection and accurate location of attack entry points or paths. Therefore, the prevailing restoration framework may struggle to fully consider the cyber-physical interdependence, successfully isolate the compromised cyber and physical components, and safely recover the systems without the potential risks leading to secondary outages. This paper is motivated to develop a cyber-physical restoration framework for distribution grids to recover from cyber attacks by harnessing grid-edge IBRs. The framework is first built on the operational guidelines of IBRs considering the compromised cyber layer. Then, an ambiguity set is established to represent the uncertainty of attack scenarios and their possibility levels. Next, a distributionally robust optimization model is developed to provide the optimal load restoration strategy across all scenarios. The effectiveness of the proposed model is demonstrated through various use cases on the modified IEEE 13-node and 123-node test systems. Finally, simulation results demonstrate the effectiveness and advancement of developed post-attack restoration strategies.

Cybersecurity↗

OC6 project Phase IV: validation of numerical models for novel floating offshore wind support structures

Abstract. This paper provides a summary of the work done within Phase IV of the Offshore Code Comparison Collaboration, Continued with Correlation and unCertainty (OC6) project, under International Energy Agency Wind Technology Collaboration Programme Task 30. This phase focused on validating the loading on and motion of a novel floating offshore wind system. Numerical models of a 3.6 MW horizontal-axis wind turbine atop the TetraSpar floating support structure were compared using measurement data from a 1:43-Froude-scale test performed in the University of Maine's Alfond Wind–Wave (W2) Ocean Engineering Laboratory. Participants in the project ran a series of simulations, including system equilibrium, surge offsets, free-decay tests, wind-only conditions, wave-only conditions, and a combination of wind and wave conditions. Validation of the models was performed by comparing the aerodynamic loading, floating support structure motion, tower base loading, mooring line tensions, and keel line tensions. The results show a relatively good estimation of the aerodynamic loading and a reasonable estimation of the platform motion and tower base fore–aft bending moment. However, there is a significant dispersion in the dynamic loading for the upwind mooring line. Very good agreement was observed between most of the numerical models and the experiment for the keel line tensions.

17 WIND ENERGY↗

A low-rank control variate for multilevel Monte Carlo simulation of high-dimensional uncertain systems

Multilevel Monte Carlo (MLMC) is a recently proposed variation of Monte Carlo (MC) simulation that achieves variance reduction by simulating the governing equations on a series of spatial (or temporal) grids with increasing resolution. Instead of directly employing the fine grid solutions, MLMC estimates the expectation of the quantity of interest from the coarsest grid solutions as well as differences between each two consecutive grid solutions. When the differences corresponding to finer grids become smaller, hence less variable, fewer MC realizations of finer grid solutions are needed to compute the difference expectations, thus leading to a reduction in the overall work. Here this paper presents an extension of MLMC, referred to as multilevel control variates (MLCV), where a low-rank approximation to the solution on each grid, obtained primarily based on coarser grid solutions, is used as a control variate for estimating the expectations involved in MLMC. Cost estimates as well as numerical examples are presented to demonstrate the advantage of this new MLCV approach over the standard MLMC when the solution of interest admits a low-rank approximation and the cost of simulating finer grids grows fast.

97 MATHEMATICS AND COMPUTING↗

Improved scaling of temperature-accelerated dynamics using localization

While temperature-accelerated dynamics (TAD) is a powerful method for carrying out non-equilibrium simulations of systems over extended time scales, the computational cost of serial TAD increases approximately as N{sup 3} where N is the number of atoms. In addition, although a parallel TAD method based on domain decomposition [Y. Shim et al., Phys. Rev. B 76, 205439 (2007)] has been shown to provide significantly improved scaling, the dynamics in such an approach is only approximate while the size of activated events is limited by the spatial decomposition size. Accordingly, it is of interest to develop methods to improve the scaling of serial TAD. As a first step in understanding the factors which determine the scaling behavior, we first present results for the overall scaling of serial TAD and its components, which were obtained from simulations of Ag/Ag(100) growth and Ag/Ag(100) annealing, and compare with theoretical predictions. We then discuss two methods based on localization which may be used to address two of the primary “bottlenecks” to the scaling of serial TAD with system size. By implementing both of these methods, we find that for intermediate system-sizes, the scaling is improved by almost a factor of N{sup 1/2}. Some additional possible methods to improve the scaling of TAD are also discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Exploring strongly correlated quantum spin systems with quantum computers

At inception, quantum annealing leveraged quantum mechanics for classical optimization tasks. However, as machine coherence increases, the realization of quantum spin systems has emerged as an even more fruitful application of this computational model. Alternatives to the standard gate model deserve exploration, to achieve useful quantum advantage via analog quantum computing. The recent demonstration of so-called coherent quantum annealing with 1,000s of qubits, further expands the potential for this hardware to explore quantum system dynamics where the effects of quantum fluctuations can be carefully controlled and observed directly. Implementing existing celebrated spin models—or in fact new and dedicated ones—into quantum annealers will lead to observation and detection of quantum phenomena not yet observed, or not visualized directly, in experimental physics laboratories. Rather than computing or simulating quantum systems, these quantum computers allow one to simply build quantum systems, and experiment on them in an uniquely controlled way, with characterization down to the constitutive degree of freedom. This provides an unprecedented opportunity for understanding the physics of quantum spin systems. Keywords: Theoretical physics currently abounds of interesting theoretical spin models to explore frustration, strongly correlated spins, spin liquids, fractionalized excitations, and topological matter. However enticing, such models are generally only weak proxies for the properties of actual materials. In quantum annealers these models could be realized and experimented upon. Moreover, many more realistic models of such materials could be realized in quantum annealers. Employed in this way, quantum annealers provide an extraordinary versatile platform to explore quantum effects that are hard to find, detect, and characterize in natural materials.

36 MATERIALS SCIENCE↗

Hardware-in-the-Loop (HIL) Simulations for Smart Grid Impact Studies: Preprint

Hardware-in-the-loop (HIL) simulation is increasingly employed in power engineering as more inverter-based generation and smart appliances are connected to the electric grid. HIL techniques allow for co-simulation of analytical models with actual devices whose complex behavior is computationally inefficient or difficult to capture. It also allows for testing of the behavior of these devices under adverse conditions that only occur rarely in the field, but that are important to evaluate. An overview of how HIL simulation has been used to date is provided. This paper further proposes that HIL simulation should play an important role in evaluating new control strategies, especially at the distribution level, that are being proposed to allow for continued affordable and reliable operation of the electric grid with the introduction of these new technologies. It presents a new capability that was developed to enable evaluation of interactions between residential loads and the smart grid: smart home hardware-in-the-loop. The paper includes results from an HIL simulation that incorporates multiple technologies and controls.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Evaluation of phase change thermal storage in a cascade heat pump

Here this study proposes a low-global warming potential cascade heat pump system with integrated phase change material (PCM) for combined space heating and cooling. PCM is embedded in the intermediate heat exchanger that couples an outdoor R290 refrigerant circuit and an indoor CO 2 circuit. We used a clustering method to develop seasonally representative load and ambient temperature profiles for residential buildings in Minneapolis, Minnesota. We simulated the system peak shaving performance on the selected representative days, assuming three different PCM options including ice (0 °C), tetradecane (5.6 °C), and a commercially available PCM, SP9 (10.5 °C). The baseline was the same cascade system without thermal storage. We investigated the effect of both thermal energy storage capacity (kWh th ) and PCM transition temperature on system performance. The results show that higher thermal energy storage capacity shaves more electric energy until it reaches a maximum value set by the load characteristics of the building. SP9, with the highest transition temperature, shaves the most electric energy, while ice shaves the least. Moreover, ice has the highest charging cost, especially in the cooling season when the outdoor circuit has to operate at an evaporating temperature lower than 0 °C to make ice. Our results show that SP9 presents the lowest levelized cost of storage, closely followed by ice. A lower off-peak electricity rate benefits ice more than SP9 or tetradecane because ice has the largest charging cost penalty.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Preliminary Designs for Criticality Safety Benchmarks—Iron/Steel/Chromium Series

The nuclear criticality safety community uses critical benchmarks extensively when analyzing nuclear systems. Having benchmarks sensitive to important materials of interest helps to ensure safety and constrain results. Critical benchmarks are also invaluable to the nuclear data community, which uses these benchmarks as a check and measure of quality of nuclear data libraries. Additionally, integral benchmarks are often used to validate the performance of computational tools for simulating nuclear systems. Many integral benchmarks have been documented in the International Criticality Safety Benchmark Evaluation Project (ICSBEP) handbook, the International Reactor Physics Experiment Evaluation Project (IRPhEP) handbook, and the Shielding Integral Benchmark Archive Database (SINBAD). The National Criticality Experiments Research Center (NCERC), operated by Los Alamos National Laboratory, is the only general purpose critical experiments facility in the United States of America. NCERC regularly designs and executes critical experiments and other measurements useful to a wide variety of fields including nuclear criticality safety, commercial nuclear energy, nonproliferation, and nuclear data. NCERC has recently developed tools to visualize gaps and deficiencies in existing criticality benchmarks – this work is still ongoing. Through this work and other reviews, some high priority deficiencies have been identified. This paper aims to highlight three areas where additional critical benchmarks would be useful to the nuclear criticality safety community and various other groups. Additionally, NCERC is seeking feedback in regard to these proposed measurements to make sure they are as useful to as many communities as possible.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NEAMS-Multiphysics MOOSE End Year Framework Activities FY20

The Multiphysics Object Oriented Simulation Environment (MOOSE) is a general finite element package meant for high performance solution of multiphysics problems in science and engineering. In this report, we present additions and enhancements to MOOSE funded by the Nuclear Engineering Advanced Modeling and Simulation (NEAMS) program. NEAMS-funded MOOSE framework library improvements include: expansion of support for multi-level multi-application restart; enhancement of coupling between native MOOSE applications and external libraries; addition of a sparse automatic differentiation (AD) container enabling non-local degree of freedom coupling; block-specific quadrature rules; further development of an eigenvalue executioner; faster setup of periodic boundary conditions; and creation of a mesh meta-data system streamlining simulation startup. Besides developing the framework library, NEAMS funds were used to overhaul a swath of important MOOSE infrastructure in order to substantially improve user experience. These critical infrastructure changes included adapting MOOSE to python 3, improving the test harness to catch race conditions, and most importantly transitioning MOOSE to use Conda, a globally known packaging system, that greatly eases software adoption.

97 MATHEMATICS AND COMPUTING↗

Simulating signatures of two-dimensional electronic spectra of the Fenna-Matthews-Olson complex: By using a numerical path integral

A framework for simulating electronic spectra from photon-echo experiments is constructed by using a numerical path integral technique. This method is non-Markovian and nonperturbative and, more importantly, is not limited by a fixed form of the spectral density functions of the environment. Next, a two-dimensional (2D) third-order electronic spectrum of a dimer system is simulated. The spectrum is in agreement with the experimental and theoretical results previously reported [for example, M. Khalil, N. Demirdöven, and A. Tokmakoff, Phys. Rev. Lett. 90, 047401 (2003)]. Finally, a 2D third-order electronic spectrum of the Fenna-Matthews-Olson (FMO) complex is simulated by using the Debye, Ohmic, and Adolphs and Renger spectral density functions. It is shown that this method can clearly produce the spectral signatures of the FMO complex by using only the Adolphs and Renger spectral density function. Plots of the evolution of the diagonal and cross-peaks show that they are oscillating with the population time.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

VAST PLANES OF SATELLITES IN A HIGH-RESOLUTION SIMULATION OF THE LOCAL GROUP: COMPARISON TO ANDROMEDA

We search for vast planes of satellites (VPoS) in a high-resolution simulation of the Local Group performed by the CLUES project, which improves significantly the resolution of previous similar studies. We use a simple method for detecting planar configurations of satellites, and validate it on the known plane of M31. We implement a range of prescriptions for modeling the satellite populations, roughly reproducing the variety of recipes used in the literature, and investigate the occurrence and properties of planar structures in these populations. The structure of the simulated satellite systems is strongly non-random and contains planes of satellites, predominantly co-rotating, with, in some cases, sizes comparable to the plane observed in M31 by Ibata et al. However, the latter is slightly richer in satellites, slightly thinner, and has stronger co-rotation, which makes it stand out as overall more exceptional than the simulated planes, when compared to a random population. Although the simulated planes we find are generally dominated by one real structure forming its backbone, they are also partly fortuitous and are thus not kinematically coherent structures as a whole. Provided that the simulated and observed planes of satellites are indeed of the same nature, our results suggest that the VPoS of M31 is not a coherent disk and that one-third to one-half of its satellites must have large proper motions perpendicular to the plane.

79 ASTRONOMY AND ASTROPHYSICS↗

Minimum critical masses and their uncertainties for several fissile nuclides and reflectors - 14296

A systematic effort was undertaken to reevaluate the minimum critical mass for several fissile nuclides and reflectors. This was done using MCNP6 in combination with several modern nuclear data libraries, viz. ENDF/B-VII.1, JENDL-4.0 and JEFF-3.1.1. This combination was benchmarked extensively against criticality safety benchmarks. The procedure to extract the minimum critical mass from the MCNP6 calculations was carefully implemented, so that it would have a negligible impact on the accuracy of the end results. All parts of the calculation, i.e. nuclear data processing, MCNP6 models, and post-processing, were applied consistently for all cases, involving varying fissile nuclides and reflectors. Results for the minimum critical mass were generated for {sup 233,235}U, {sup 236,239,241}Pu, {sup 242m}Am, {sup 243,245,247}Cm, and {sup 249,251}Cf, for the case of water moderation and water reflection. For {sup 235}U, the minimum critical mass was also calculated for cases with other reflectors, such as beryllium, heavy water, graphite, lead, and concrete. The results for the water reflection cases are in good agreement with most of the literature values, but for the cases with reflectors such as graphite, lead, or concrete, there are sizable differences. The accuracy of the results was assessed in two ways. First, a comparison was made with criticality safety benchmark results, by selecting benchmarks that resemble the minimum mass systems the most. It was clear that for {sup 233,235}U and {sup 239}Pu with water reflection, the differences in minimum critical masses between the libraries are consistent with the differences obtained for benchmark cases. The difference between the benchmark calculations and the benchmark values is then an indication of the accuracy with which one can simulate comparable systems. Secondly, sensitivity runs were performed with the new MCNP6 sensitivity option, based on ENDF/B-VII.1 nuclear data. For each fissile nuclide, the sensitivity of k{sub eff} to cross sections such as capture, fission, elastic, inelastic was calculated for the minimum critical mass case. The same was done for the {sup 1}H and {sup 16}O cross sections in the moderator and reflector. The uncertainty in k{sub eff} was calculated by combining the sensitivity profiles with the appropriate ENDF/B-VII.1 covariance matrix, using the SUSD code. As a final step, the uncertainty in k{sub eff} can be translated to an uncertainty in the minimum critical mass. (author)

42 ENGINEERING↗

Computer Simulation of Critical Behavior of Two-Dimensional Weakly Diluted Antiferromagnetic Potts Model on a Triangular Lattice

The critical behavior of a two-dimensional weakly diluted antiferromagnetic Potts model with spin states q = 3 on a triangular lattice is computationally simulated. The systems with periodic boundary conditions were calculated, when spin concentrations p were 1.0 and 0.9. Systems with linear dimensions L × L = N and L = 9–144 were considered. Static critical exponents of heat capacity α, the susceptibility γ, the order parameter β, and critical exponent ν for a correlation radius are calculated based on a theory of finite-dimensional scaling.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Scaled melter testing of alternative reductants for Low-Activity Waste vitrification: Melter feed properties and processing behavior

This study examines the use of boron nitride (BN) and coke dust as alternative reductants for low-activity waste vitrification in scaled melter systems, intended to enhance operational flexibility at the Hanford Waste Treatment and Immobilization Plant. The baseline reductant, sucrose, while widely considered advantageous for glass production rates, produces challenging off gas emissions during processing. Early crucible-scale testing indicated that BN and coke effectively control foaming and glass redox, while significantly reducing acetonitrile emissions. Both real tank waste and simulant waste melter feeds were tested with BN and coke in this study. Rheological properties for melter feeds with both reductants were within operational limits, although the particle size of coke should be controlled and sufficient agitation is necessary to ensure smooth operations. While promising for handling and reduction of emissions, both simulant melter feeds demonstrated lower average glass production rates in the simulant melter system compared with sucrose. In the radioactive system, the rates were more comparable to previously measured low activity waste feeds on the order of 1285 kg m 2 day −1 for BN and 1671 kg m −2 day −1 for coke. Further scale-up testing is required to interpret the processing rates in the context of the full scale melter systems and evaluate the cost-benefit impact of these reductants on efficiency and sustainability.

Rigby, Jessica C. [Pacific Northwest National Labo↗

Scheduling strategies for the ESPRESSO follow-up of TESS targets

ABSTRACT Radial-velocity follow-up of stars harbouring transiting planets detected by TESS is expected to require very large amounts of expensive telescope time in the next few years. Therefore, scheduling strategies should be implemented to maximize the amount of information gathered about the target planetary systems. We consider myopic and non-myopic versions of a novel uniform-in-phase scheduler, as well as a random scheduler, and compare these scheduling strategies with respect to the bias, accuracy and precision achieved in recovering the mass and orbital parameters of transiting and non-transiting planets. This comparison is carried out based on realistic simulations of radial-velocity follow-up with ESPRESSO of a sample of 50 TESS target stars, with simulated planetary systems containing at least one transiting planet with a radius below 4R⊕. Radial-velocity data sets were generated under reasonable assumptions about their noise component, including that resulting from stellar activity, and analysed using a fully Bayesian methodology. We find the random scheduler leads to a more biased, less accurate, and less precise, estimation of the mass of the transiting exoplanets. No significant differences are found between the results of the myopic and non-myopic implementations of the uniform-in-phase scheduler. With only about 22 radial velocity measurements per data set, our novel uniform-in-phase scheduler enables an unbiased (at the level of 1 per cent) measurement of the masses of the transiting planets, while keeping the average relative accuracy and precision around 16 per cent and 23 per cent, respectively. The number of non-transiting planets detected is similar for all the scheduling strategies considered, as well as the bias, accuracy and precision with which their masses and orbital parameters are recovered.

Cabona, L.↗