Search NASA⌕ Search

SEARCH · Search NASA

Results for “standard problem”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2012

MARS-KS code validation activity through the atlas domestic standard problem

The 2 nd Domestic Standard Problem (DSP-02) exercise using the ATLAS integral effect test data was executed to transfer the integral effect test data to domestic nuclear industries and to contribute to improving the safety analysis methodology for PWRs. A small break loss of coolant accident of a 6-inch break at the cold leg was determined as a target scenario by considering its technical importance and by incorporating interests from participants. Ten calculation results using MARS-KS code were collected, major prediction results were described qualitatively and code prediction accuracy was assessed quantitatively using the FFTBM. In addition, special code assessment activities were carried out to find out the area where the model improvement is required in the MARS-KS code. The lessons from this DSP-02 and recommendations to code developers are described in this paper. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2010

Standard Problems for CFD Validation for NGNP - Status Report

The U.S. Department of Energy (DOE) is conducting research and development to support the resurgence of nuclear power in the United States for both electrical power generation and production of process heat required for industrial processes such as the manufacture of hydrogen for use as a fuel in automobiles. The project is called the Next Generation Nuclear Plant (NGNP) Project, which is based on a Generation IV reactor concept called the very high temperature reactor (VHTR). The VHTR will be of the prismatic or pebble bed type; the former is considered herein. The VHTR will use helium as the coolant at temperatures ranging from 250°C to perhaps 1000°C. While computational fluid dynamics (CFD) has not previously been used for the safety analysis of nuclear reactors in the United States, it is being considered for existing and future reactors. It is fully recognized that CFD simulation codes will have to be validated for flow physics reasonably close to actual fluid dynamic conditions expected in normal operational and accident situations. The “Standard Problem” is an experimental data set that represents an important physical phenomenon or phenomena, whose selection is based on a phenomena identification and ranking table (PIRT) for the reactor in question. It will be necessary to build a database that contains a number of standard problems for use to validate CFD and systems analysis codes for the many physical problems that will need to be analyzed. The first two standard problems that have been developed for CFD validation consider flow in the lower plenum of the VHTR and bypass flow in the prismatic core. Both involve scaled models built from quartz and designed to be installed in the INL’s matched index of refraction (MIR) test facility. The MIR facility employs mineral oil as the working fluid at a constant temperature. At this temperature, the index of refraction of the mineral oil is the same as that of the quartz. This provides an advantage to the optics used for data gathering. Particle image velocimetry (PIV) is used to take the data. The first standard problem represents several flow physics expected to be present in the lower plenum of the prismatic VHTR. In the lower plenum, heated helium coolant in the form of jets issues downward into the plenum and is then forced to turn ninety degrees and flow toward the exit duct. The lower plenum is filled with cylindrical graphite posts that hold up the core. Figure S-1 provides a plan view of the lower plenum. The red circles represent support posts holding up columns of heated blocks. Grey circles represent support posts under columns of reflector blocks. Helium enters the lower plenum at the junctions of the hexagonal blocks.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2011

Pre-test CFD Calculations for a Bypass Flow Standard Problem

The bypass flow in a prismatic high temperature gas-cooled reactor (HTGR) is the flow that occurs between adjacent graphite blocks. Gaps exist between blocks due to variances in their manufacture and installation and because of the expansion and shrinkage of the blocks from heating and irradiation. Although the temperature of fuel compacts and graphite is sensitive to the presence of bypass flow, there is great uncertainty in the level and effects of the bypass flow. The Next Generation Nuclear Plant (NGNP) program at the Idaho National Laboratory has undertaken to produce experimental data of isothermal bypass flow between three adjacent graphite blocks. These data are intended to provide validation for computational fluid dynamic (CFD) analyses of the bypass flow. Such validation data sets are called Standard Problems in the nuclear safety analysis field. Details of the experimental apparatus as well as several pre-test calculations of the bypass flow are provided. Pre-test calculations are useful in examining the nature of the flow and to see if there are any problems associated with the flow and its measurement. The apparatus is designed to be able to provide three different gap widths in the vertical direction (the direction of the normal coolant flow) and two gap widths in the horizontal direction. It is expected that the vertical bypass flow will range from laminar to transitional to turbulent flow for the different gap widths that will be available.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Leptogenesis in parity solutions to the strong CP problem and Standard Model parameters

We study the simplest theories with exact spacetime parity that solve the strong CP problem and successfully generate the cosmological baryon asymmetry via decays of right-handed neutrinos. Lower bounds are derived for the masses of the right-handed neutrinos and for the scale of spontaneous parity breaking, v R . For generic thermal leptogenesis, v R ≳ 10 12 GeV, unless the small observed neutrino masses arise from fine-tuning. We compute v R in terms of the top quark mass, the QCD coupling, and the Higgs boson mass and find this bound is consistent with current data at 1σ. Future precision measurements of these parameters may provide support for the theory or, if v R is determined to be below 10 12 GeV, force modifications. However, modified cosmologies do not easily allow reductions in v R — no reduction is possible if leptogenesis occurs in the collisions of domain walls formed at parity breaking, and at most a factor 10 reduction is possible with non-thermal leptogenesis. Standard Model parameters that yield low values for v R can only be accommodated by having a high degree of degeneracy among the right-handed neutrinos involved in leptogenesis. If future precision measurements determine v R to be above 10 12 GeV, it is likely that higher-dimensional operators of the theory will yield a neutron electric dipole moment accessible to ongoing experiments. This is especially true in a simple UV completion of the neutrino sector, involving gauge singlet fermions, where the bound from successful leptogenesis is strengthened to v R ≳ 10 13 GeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2013

RELAP5-3D Restart and Backup Verification Testing

Existing testing methodology for RELAP5-3D employs a set of test cases collected over two decades to test a variety of code features and run on a Linux or Windows platform. However, this set has numerous deficiencies in terms of code coverage, detail of comparison, running time, and testing fidelity of RELAP5-3D restart and backup capabilities. The test suite covers less than three quarters of the lines of code in the relap directory and just over half those in the environmental library. Even in terms of code features, many are not covered. Moreover, the test set runs many problems long past the point necessary to test the relevant features. It requires standard problems to run to completion. This is unnecessary for features can be tested in a short-running problem. For example, many trips and controls can be tested in the first few time steps, as can a number of fluid flow options. The testing system is also inaccurate. For the past decade, the diffem script has been the primary tool for checking that printouts from two different RELAP5-3D executables agree. This tool compares two output files to verify that all characters are the same except for those relating to date, time and a few other excluded items. The variable values printed on the output file are accurate to no more than eight decimal places. Therefore, calculations with errors in decimal places beyond those printed remain undetected. Finally, fidelity of restart is not tested except in the PVM sub-suite and backup is not specifically tested at all. When a restart is made from any midway point of the base-case transient, the restart must produce the same values. When a backup condition occurs, the code repeats advancements with the same time step. A perfect backup can be tested by forcing RELAP5 to perform a backup by falsely setting a backup condition flag at a user-specified-time. Comparison of the calculations of that run and those produced by the same input w/o the spurious condition should be identical. Backup testing is more difficult the other kinds of testing described above because it requires additional coding to implement. The testing system constructed and described in this document resolves all of these issues. A matrix of test features and short-running cases that exercise them is presented. A small information file that contains sufficient data to verify calculations to the last decimal place and bit is produced. This testing system is used to test base cases (called null testing) as well as restart and backup cases. The programming that implements these new capabilities is presented.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2018

Phase field benchmark problems for dendritic growth and linear elasticity

We present the second set of benchmark problems for phase field models that are being jointly developed by the Center for Hierarchical Materials Design (CHiMaD) and the National Institute of Standards and Technology (NIST) along with input from other members in the phase field community. As the integrated computational materials engineering (ICME) approach to materials design has gained traction, there is an increasing need for quantitative phase field results. New algorithms and numerical implementations increase computational capabilities, necessitating standard problems to evaluate their impact on simulated microstructure evolution as well as their computational performance. We propose one benchmark problem for solidifiication and dendritic growth in a single-component system, and one problem for linear elasticity via the shape evolution of an elastically constrained precipitate. We demonstrate the utility and sensitivity of the benchmark problems by comparing the results of 1) dendritic growth simulations performed with different time integrators and 2) elastically constrained precipitate simulations with different precipitate sizes, initial conditions, and elastic moduli. As a result, these numerical benchmark problems will provide a consistent basis for evaluating different algorithms, both existing and those to be developed in the future, for accuracy and computational efficiency when applied to simulate physics often incorporated in phase field models.

36 MATERIALS SCIENCE↗
DOE OSTI2018

A note on the S-matrix bootstrap for the 2d O(N) bosonic model

In this work we apply the S-matrix bootstrap maximization program to the 2d bosonic O ( N ) integrable model which has N species of scalar particles of mass m and no bound states. Since in previous studies theories were defined by maximizing the coupling between particles and their bound states, the main problem appears to be to find what other functional can be used to define this model. Instead, we argue that the defining property of this integrable model is that it resides at a vertex of the convex space determined by the unitarity and crossing constraints. Thus, the integrable model can be found by maximizing any linear functional whose gradient points in the general direction of the vertex, namely within a cone determined by the normals to the faces intersecting at the vertex. This is a standard problem in applied mathematics, related to semi-definite programming and solvable by fast available numerical algorithms. The information provided by the numerical solution is enough to reproduce the known analytical solution without using integrability, namely the Yang-Baxter equation. This situation seems quite generic so we expect that other theories without continuous parameters can also be found by maximizing linear functionals in the convex space of allowed S-matrices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2010

Dark matter decaying into a Fermi sea of neutrinos

We study the possible decay of a coherently oscillating scalar field, interpreted as dark matter, into light fermions. Specifically, we consider a scalar field with sub-eV mass decaying into a Fermi sea of neutrinos. We recognize the similarity between our scenario and inflationary preheating where a coherently oscillating scalar field decays into standard model particles. Like the case of fermionic preheating, we find that Pauli blocking controls the dark matter decay into the neutrino sea. The radius of the Fermi sphere depends on the expansion of the universe leading to a time varying equation of state of dark matter. This makes the scenario very rich and we show that the decay rate might be different at different cosmological epochs. We categorize this in two interesting regimes and then study the cosmological perturbations to find the impact on structure formation. We find that the decay may help in alleviating some of the standard problems related to cold dark matter.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

A precision test of averaging in AdS/CFT

We reconsider the role of wormholes in the AdS/CFT correspondence. We focus on Euclidean wormholes that connect two asymptotically AdS or hyperbolic regions with $\mathbb{S}$ 1 × $\mathbb{S}$ d-1 boundary. There is no solution to Einstein’s equations of this sort, as the wormholes possess a modulus that runs to infinity. To find on-shell wormholes we must stabilize this modulus, which we can do by fixing the total energy on the two boundaries. Such a wormhole gives the saddle point approximation to a non-standard problem in quantum gravity, where we fix two asymptotic boundaries and constrain the common energy. Crucially the dual quantity does not factorize even when the bulk is dual to a single CFT, on account of the fixed energy constraint. From this quantity we extract a smeared version of the microcanonical spectral form factor. For a chaotic theory this quantity is self-averaging, i.e. well-approximated by averaging over energy windows, or over coupling constants. We go on to give a precision test involving the microcanonical spectral form factor where the two replicas have slightly different coupling constants. In chaotic theories this form factor is known to smoothly decay at a rate universally predicted in terms of one replica physics, provided that there is an average either over a window or over couplings. We compute the expected decay rate for holographic theories, and the form factor from a wormhole, and the two exactly agree for a wide range of two-derivative effective field theories in AdS. This gives a precision test of averaging in AdS/CFT. Our results interpret a number of confusing facts about wormholes and factorization in AdS and suggest that we should regard gravitational effective field theory as a mesoscopic description, analogous to semiclassical mesoscopic descriptions of quantum chaotic systems.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Applications of the worldline Monte Carlo formalism in quantum mechanics

Highlights: Universal approach to the numerical computation of the Euclidean path integral. Inspired by recent work in relativistic quantum field theory. Here adapted to non-relativistic quantum mechanics. Worked out for the computation of propagators and ground-state energies. Special smoothing procedure for singular potentials. -- Abstract: In recent years efficient algorithms have been developed for the numerical computation of relativistic single-particle path integrals in quantum field theory. Here, we adapt this “worldline Monte Carlo” approach to the standard problem of the numerical approximation of the non-relativistic path integral, resulting in a formalism whose characteristic feature is the fast, non-recursive generation of an ensemble of trajectories that is independent of the potential, and is thus universally applicable. The numerical implementation discretises the trajectories with respect to their time parametrisation but maintains a continuous spatial domain. In the case of singular potentials, the discretised action gets adapted to the singularity through a “smoothing” procedure. We show for a variety of examples (the harmonic oscillator in various dimensions, the modified Pöschl–Teller potential, delta-function potentials, the Coulomb and Yukawa potentials) that the method allows one to obtain fast and reliable estimates for the Euclidean propagator and use them in a certain time window suitable for extracting the ground state energy. Usually Monte Carlo calculations in quantum mechanics are performed using algorithms that involve a more specific adaptation to the potential, and we briefly compare our method with the most widely used such algorithms, which are of Metropolis type. As an aside, we apply it for studying the classical limit where nearly classical trajectories are expected to dominate in the path integral. We expect the advances made here to be useful also in the relativistic case.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Benchmarking impact hydrocodes in the strength regime: Implications for modeling deflection by a kinetic impactor

The Double Asteroid Redirection Test (DART) is a NASA-sponsored mission that will be the first direct test of the kinetic impactor technique for planetary defense. Here, the DART spacecraft will impact into Didymos-B, the moon of the binary system 65803 Didymos and the resulting period change will be measured from Earth. Impact simulations will be used to predict the crater size and momentum enhancement expected from the DART impact. Because the specific material properties (strength, porosity, internal structure) of the Didymos-B target are unknown, a wide variety of numerical simulations must be performed to better understand possible impact outcomes. This simulation campaign will involve a large parameter space being simulated using multiple different shock physics hydrocodes. In order to understand better the behaviors and properties of numerical simulation codes applicable to the DART impact, a benchmarking and validation program using different numerical codes to solve a set of standard problems was designed and implemented.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Dynamics, stability, and statistics on lattices and networks

These lectures aim at surveying some dynamical models that have been widely explored in the recent scientific literature as case studies of complex dynamical evolution, emerging from the spatio-temporal organization of several coupled dynamical variables. The first message is that a suitable mathematical description of such models needs tools and concepts borrowed from the general theory of dynamical systems and from out-of-equilibrium statistical mechanics. The second message is that the overall scenario is definitely reacher than the standard problems in these fields. For instance, systems exhibiting complex unpredictable evolution do not necessarily exhibit deterministic chaotic behavior (i.e., Lyapunov chaos) as it happens for dynamical models made of a few degrees of freedom. In fact, a very large number of spatially organized dynamical variables may yield unpredictable evolution even in the absence of Lyapunov instability. Such a mechanism may emerge from the combination of spatial extension and nonlinearity. Moreover, spatial extension allows one to introduce naturally disorder, or heterogeneity of the interactions as important ingredients for complex evolution. It is worth to point out that the models discussed in these lectures share such features, despite they have been inspired by quite different physical and biological problems. Along these lectures we describe also some of the technical tools employed for the study of such models, e.g., Lyapunov stability analysis, unpredictability indicators for “stable chaos,” hydrodynamic description of transport in low spatial dimension, spectral decomposition of stochastic dynamics on directed networks, etc.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

CLOAKED GAMMA-RAY BURSTS

It is suggested that many gamma-ray bursts (GRBs) are cloaked by an ultrarelativistic baryonic shell that has high optical depth when the photons are manufactured. Such a shell would not fully block photons reflected or emitted from its inner surface, because the radial velocity of the photons can be less than that of the shell. This avoids the standard problem associated with GRBs that the thermal component should be produced where the flow is still obscured by high optical depth. The radiation that escapes high optical depth obeys the Amati relation. Observational implications may include (1) anomalously high ratios of afterglow to prompt emission, such as may have been the case in the recently discovered PTF 11agg, and (2) ultrahigh-energy neutrino pulses that are non-coincident with detectable GRB. It is suggested that GRB 090510, a short, very hard GRB with very little afterglow, was an exposed GRB, in contrast to those cloaked by baryonic shells.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Non-linear Model Predictive Control of Wave Energy Converters with Realistic Power Take-off Configurations and Loss Model

Model Predictive Control (MPC) has been recognized as a well-developed control strategy for optimizing the performance of wave energy converters (WECs). The standard problem involves maximizing the absorbed power from oncoming waves while respecting motion constraints of the device and force constraints of the power take-off (PTO). The generated electrical power is then calculated by assuming an ideal PTO with no power conversion losses or limitations to power flow. In this article, we remove the assumption of an ideal PTO and present four different control options that reflect unique combinations of PTO control and power flow constraints. We also recast the MPC problem to maximize the generated power instead of the absorbed power by accounting for the losses in the power train. Using a WEC with a hydraulic PTO, we develop a loss model suitable for controls optimization. Lastly, we illustrate the process of setting up the non-linear MPC problem and compare the performance of competing control strategies using a common benchmark based on annual energy captured at a reference site.

16 TIDAL AND WAVE POWER↗
DOE OSTI2010

Dynamics of a vertical flight in the stationary gravitational field of a celestial body: Post-newtonian corrections and gravitational redshift

The standard problem of a radial motion of test particles in the stationary gravitational field of a spherically symmetric celestial body is solved and is used to determine the time features of this motion. The problem is solved for the equations of motion of general relativity (GR), and the time features are obtained in the post-Newtonian approximation, with linear GR corrections proportional to r{sub g}/r and {beta}{sup 2} (in the solution being considered, they are of the same order of smallness) being taken rigorously into account. Total times obtained by integrating the time differentials along the trajectories of motion are considered as the time features in question. It is shown that, for any parameters of the motion, the proper time (which corresponds to watches comoving with a test particle) exceeds the time of watches at rest (watches at the surface of the celestial body being considered). The mass and the radius of the celestial body, as well as the initial velocity of the test particle, serve as arbitrary parameters of the motion. The time difference indicated above implies a leading role of the gravitational redshift, which decreases somewhat because of the opposite effect of the Doppler shift. The results are estimated quantitatively for the important (from the experimental point of view) case of vertical flights of rockets starting from the Earth's surface. In this case, the GR corrections, albeit being extremely small (a few microseconds for several hours of the flight), aremeasurable with atomic (quantum) watches.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Quantum signatures of chaos or quantum chaos?

A critical analysis of the present-day concept of chaos in quantum systems as nothing but a “quantum signature” of chaos in classical mechanics is given. In contrast to the existing semi-intuitive guesses, a definition of classical and quantum chaos is proposed on the basis of the Liouville–Arnold theorem: a quantum chaotic system featuring N degrees of freedom should have M < N independent first integrals of motion (good quantum numbers) specified by the symmetry of the Hamiltonian of the system. Quantitative measures of quantum chaos that, in the classical limit, go over to the Lyapunov exponent and the classical stability parameter are proposed. The proposed criteria of quantum chaos are applied to solving standard problems of modern dynamical chaos theory.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2018

A Nonlinear Algebraic Multigrid Framework for the Power Flow Equations

Multigrid is a highly scalable class of methods most often used for solving large linear systems. In this paper we develop a nonlinear algebraic multigrid framework for the power flow equations, a complex quadratic system of the form ${diag}({v})\overline{Y{v}}={s}$, where $Y$ is approximately a complex scalar rotation of a real graph Laplacian. This is a standard problem that needs to be solved repeatedly during power grid simulations. A key difference between our multigrid framework and typical multigrid approaches is the use of a novel multiplicative coarse-grid correction to enable a dynamic multigrid hierarchy. We also develop a new type of smoother that allows one to coarsen together the different types of nodes that appear in power grid simulations. In developing a specific multigrid method, one must make a number of choices that can significantly affect the method's performance, such as how to construct the restriction and interpolation operators, what smoother to use, and how aggressively to coarsen. In this paper, we make simple but reasonable choices that result in a scalable and robust power flow solver. Experiments demonstrate this scalability and show that it is significantly more robust to poor initial guesses than current state-of-the-art solvers.

97 MATHEMATICS AND COMPUTING↗