Search NASA⌕ Search

SEARCH · Search NASA

Results for “STABLE OSCILLATION”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Experimental and numerical study of the shielding effect of two tandem rough cylinders in flow-induce oscillation

The shielding effect of the downstream cylinder in flow induced oscillation (FIO) of two cylinders arranged in tandem is studied experimentally and numerically at Reynolds number 30,000 to 120,000. Both cylinders are in one degree-of-freedom, transverse-oscillations, and have turbulence stimulation in the form of selective surface roughness to expand FIO beyond vortex-induced vibration (VIV) into galloping. Shielding of the downstream cylinder has a negative effect on harnessing hydrokinetic energy. Further, to study its effect and mechanics, selective cases are studied both numerically and experimentally and discussed to demonstrate the shielding effect on the downstream cylinder and understand its cause. The main conclusions are: (1) The shielding effect for the downstream cylinder shows a strong relation to the damping ratio. As the damping ratio increases, the shielding effect is mitigated. Additionally, the oscillation of the rear cylinder becomes stable and shows stable frequency. (2) In the VIV region, as the stiffness and natural frequency increase, the shielding effect decreases substantially. (3) In the VIV region, the vorticity of the vortices shedding from both the upper and the lower sides of the downstream cylinder does not accumulate enough due to the attraction by the vortices shed from the upstream cylinder, thus resulting in partial suppression of the oscillation on the downstream cylinder. (4) In the galloping region, the shielding effect for the downstream cylinder depends on whether the vorticity near the downstream cylinder is strengthened by the vortices generated by the shear layers of the upstream cylinder or weakened.

42 ENGINEERING↗
DOE OSTI2021

Structure and dynamical properties of two-dimensional dusty plasmas on one-dimensional periodic substrates

In the presence of the substrate, various structural and dynamical properties of two-dimensional dusty plasma (2DDP) were investigated using Langevin dynamical simulations. This paper reviews a series of results of the structural and dynamical properties of 2DDP modified by one-dimensional periodic substrates (1DPSs) as follows. First, when the depth of the 1DPS increases gradually, it is found that the static structural order of 2DDP along each potential well of the 1DPS increases first, and then decreases gradually. When the width of the 1DPS increases gradually, the particle arrangement in each potential well of the 1DPS changes from one straight row gradually to two rows, including the stable zigzag structure. Second, when there are two rows in each potential well of the 1DPS, the phonon spectra are split into two branches, corresponding to the breathing and sloshing modes, respectively. In addition, due to the pure repulsion between dust particles, from the obtained phonon spectra, the sloshing wave propagates backward at small wave numbers. Third, the calculated mean-square displacement (MSD) shown that, at the intermediate timescale between the initial ballistic and final diffusive motions, the particle motion tends to be more subdiffusive while the depth of the 1DPS increases. While the width of the 1DPS increases gradually, the long-time diffusive motion first increases, then decreases, and finally increases again, exhibiting the oscillation-like diffusion, due to the stable zigzag structure. Finally, when an external direct-current driving force is applied on all particles of 2DDP, three different depinning phases are discovered, which are the pinned, disordered plastic flow, and moving ordered states, respectively, as the driving force increases from zero. In addition, the continuous/discontinuous property of the phase transition between these different depinning states is investigated, showing that the transition property is modulated by the depth of the applied 1DPS

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Pressure-stabilized fixed-stress iterative solutions of compositional poromechanics

We consider the numerical behavior of the fixed-stress splitting method for coupled poromechanics as undrained regimes are approached. We explain that pressure stability is related to the splitting error of the scheme, not the fact that the discrete saddle point matrix never appears in the fixed-stress approach. This observation reconciles previous results regarding the pressure stability of the splitting method. Using examples of compositional poromechanics with application to geological CO sequestration, we see that solutions obtained using the fixed-stress scheme with a low order finite element-finite volume discretization which is not inherently inf-sup stable can exhibit the same pressure oscillations obtained with the corresponding fully implicit scheme. Moreover, pressure jump stabilization can effectively remove these spurious oscillations in the fixed-stress setting, while also improving the efficiency of the scheme in terms of the number of iterations required at every time step to reach convergence.

42 ENGINEERING↗
DOE OSTI2024

Three-wave coupling observed between a shear Alfvén wave and a kink-unstable magnetic flux rope

Results from a laboratory experiment are presented in which, for the first time, a shear Alfvén wave is launched using an antenna in a current-carrying plasma column that is tailored to be either stable or unstable to the kink oscillation. As the plasma is driven kink unstable, the frequency power spectrum of the Alfvén wave evolves from a single peak to a peak with multiple sidebands separated by integer multiples of the kink frequency. The main sidebands (one on either side of the launched wave peak in the power spectrum) are analyzed using azimuthal wavenumber matching, perpendicular and parallel wavenumber decomposition, and bispectral time series analysis. The dispersion relation and three-wave matching conditions are satisfied, given each sideband is a propagating Alfvén wave that results from the interaction of the pump Alfvén wave and the co-propagating component of a half-wavelength, standing kink mode. The interaction is shown to generate smaller perpendicular wavelength Alfvén waves that drive energy transport to scales that will approach the dissipation scale of k⊥ρs=1, with k⊥ being the perpendicular wavenumber and ρs being the ion gyroradius at the electron temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

On the number of stable solutions in the Kuramoto model

We consider a system of n coupled oscillators described by the Kuramoto model with the dynamics given by θ˙=ω+Kf(θ). In this system, an equilibrium solution θ∗ is considered stable when ω+Kf(θ∗)=0, and the Jacobian matrix Df(θ∗) has a simple eigenvalue of zero, indicating the presence of a direction in which the oscillators can adjust their phases. Additionally, the remaining eigenvalues of Df(θ∗) are negative, indicating stability in orthogonal directions. A crucial constraint imposed on the equilibrium solution is that |Γ(θ∗)|≤π, where |Γ(θ∗)| represents the length of the shortest arc on the unit circle that contains the equilibrium solution θ∗. We provide a proof that there exists a unique solution satisfying the aforementioned stability criteria. This analysis enhances our understanding of the stability and uniqueness of these solutions, offering valuable insights into the dynamics of coupled oscillators in this system.

Arenas, Alex (ORCID:0000000309370334)↗
DOE OSTI2021

Improved Control Strategy of Grid Interface for EV High-Power Dynamic Wireless Charging

This paper proposes a control strategy for the grid interface converter in high-power dynamic wireless charging system (DWCS) to address two issues on distribution network integration. Due to the unique pulsating load profile of DWCS, load transient response capability is critical for the grid interface to maintain the dc-bus voltage stable. Besides, the inherent unbalanced situation of distribution network would lead to 2nd-order oscillations on the dc-bus voltage, which would further affect the stable operation of the entire system. In this paper, the DWCS model is developed, and the relationship between the dc-bus voltage and the input/output power is analyzed. Based on the developed model, a control strategy based on direct power control is presented. Both simulation results and hardware-in-the-loop (HIL) results demonstrate that the proposed control strategy not only improves load transient response capability, but also eliminates the 2nd-order oscillations on the dc-bus voltage under imbalanced distribution network conditions.

Zeng, Rong↗
DOE OSTI2026

Harnessing the Second-Order Metal−Insulator Transition for Neuromorphic Computing

Vanadium oxides are widely studied phase change materials for brain-inspired computing architectures. Systems like VO 2 and V 2 O 3 exhibit first-order metal−insulator transitions (MITs) with hysteresis and percolative switching, increasing stochasticity and device variability. Here, we focus on the less explored Magnéli phase V 4 O 7 , which undergoes a continuous, non-hysteretic, second-order MIT. This surprisingly enables highly reproducible volatile resistive switching in spiking-neuron-type devices. We synthesize V 4 O 7 films, characterize their structural and transport properties, and demonstrate voltage and current-driven threshold switching with electrothermal feedback. In a Pearson–Anson oscillator, V 4 O 7 devices produce stable, tunable spiking across 20–200 kHz, with consistent operation among multiple devices. We introduce a numerical analog leaky-integrate-and-fire (aLIF) model that captures waveform shapes and their dependence on load resistance, temperature, and voltage. Furthermore, these findings suggest that second-order MIT materials like V 4 O 7 are promising for deterministic, scalable spiking neuron arrays for neuromorphic computing.

V4O7↗
DOE OSTI2023

ESR Dipole Power Supply Current Ripple and Noise Specifications

This note presents key findings for the ESR main magnet dipole power supplies (PS), where we find the current ripple specification to be close to or beyond the state-of-the-art. These specifications originate from beam-beam considerations, with the requirement to limit the ripple-induced hadron emittance growth to below 10%/hour. Beam dynamics that drive this PS ripple specification arise from the beam motions at the Interaction Point (IP). The frequency of the motions can be separated into "low", compared to the betatron frequency, and "high", i.e. around the betatron frequency and harmonics. In terms of the driving frequency, "low" implies ƒ<<ƒ 0 v x,y and "fast" means ƒ≈{ƒ 0 v x,y , ƒ 0 (1- v x,y , etc.}, where ƒ 0 =1/ T 0 =78.2 kHz is the revolution frequency, and v x,y are the fractional parts of the betatron tunes. Frequencies higher than ƒ 0 /2 will be folded back due to the particles sampling the field once per turn. To provide flexibility for future lattice adjustments and working point variations, we do not consider the tunes as fixed. Instead, we assume a certain margin and allow them to potentially fall within the range of 0.1< v x,y <0.5. In other words, the high-frequency region spans approximately from 8kHz to 40 kHz. Consequently,, we definte the dipole PS ripple in two distinct frequency ranges: the low-frequency range of [1-8000] Hz and the high-frequency range of [8-40] kHz. For the physics effects we analyzed in this note, there is no distinction between the ripple (which can be approximately reproduced in the frequency domain) and the random noise if both have some power within the frequency bandwidth of interest. Therefore, while we will use the term "ripple" for short, it should always be understood that we are referring to "ripple plus noise." Except for the lower end of the low-frequency range, the impact of the rippling PS current on the beam will be considerably reduced due to the eddy currents induced in the walls of the vacuum chamber. We will account for this effect in the PS ripple specifications to follow. The remaining sections of this note are structure as follows: Section II outlines the beam-beam physics requirements for the positional stability of the beam at the IP. Section III describes the anticipated shielding effect of the eddy currents induced in the vacuum chamber. Section IV derives the ripple requirement for the low-frequency range by propagating the closed orbit ripple resulting from the rippling dipoles to the IP (relevant lattice simulation results are summarized in the Appendix). In Section V, we present the analytical criterion for the ripple in the high-frequency range by considering resonant oscillations of the electron beam around a stable closed orbit near the betatron frequency. Finally, Section VI provides a summary of our findings and discusses related work.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Confrontation of Observation and Theory: High-frequency QPOs in X-Ray Binaries, Tidal Disruption Events, and Active Galactic Nuclei

In this paper, we compile observations of high-frequency quasi-periodic oscillations (QPOs) around black holes, both stellar and supermassive, and compare their positions in the parameter space of black hole mass, spin, and oscillation frequency. We find that supermassive black holes occupy a separate region of parameter space than stellar, and further, that QPOs seen around tidal disruption events rather than Seyfert-type active galactic nuclei (AGN) occupy an entirely different space. We then compare these results to the orbital resonance, diskoseismic, warped disk, and disk–jet coupling theoretical models for the origin of high-frequency QPOs. We find that while oscillations around stellar mass black holes are generally consistent with the above models, supermassive black holes are decidedly not. Oscillations seen in tidal disruption events are consistent with oscillations near the frequency of the innermost stable circular orbit, while QPOs in AGN are not accounted for by any of the physical models in consideration. This indicates that despite the scale invariance of accretion processes implied by a decades-wide correlation between QPO frequency and black hole mass, any theory of high-frequency QPOs must relate the frequency to more than just the mass and spin.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Lasing Below 170 nm Using an Oscillator FEL

The short wavelength operation of free-electron laser (FEL) oscillators is limited by the availability of high-reflectivity, thermally stable, and radiation-resistant FEL mirrors in the VUV wavelength. We report our recent work to extend the shortest lasing wavelength of the oscillator FEL to 168.6 nm using a storage ring FEL. This progress has been made possible by developing a new FEL configuration with substantially reduced undulator harmonic radiation on the FEL mirror, a thermally stable FEL optical cavity, and a new type of high-reflectivity fluoride-based multilayer coating with a protective capping layer. Using these fluoride-based mirrors, we have demonstrated storage ring FEL lasing from 168.6 to 179.7 nm with excellent beam stability. Employing this VUV FEL in Compton scattering, we have produced the first 120 MeV gamma rays at the High Intensity Gamma-ray Source (HIGS). Furthermore, operating the HIGS in this new high-energy region will create many new opportunities for photonuclear physics research, in particular, the low-energy QCD research.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Synthesis and Single Crystals of Refractory Oxides of Lanthanides and Thorium

At the completion of this program, we can report that we developed a considerable degree of technical improvements in our ability to perform hydrothermal reactions at high temperatures and pressures. We can now routinely perform reactions at 700-750°C and 200 MPa. Currently we are in the process of exploiting this new technology synthesizing a range of exotic new materials investigating relatively poorly understood materials. Our initial efforts focused on the chemistry of rare earth oxides with tetravalent and pentavalent oxides. We recently published a study of the lanthanides with Nb 5+ and Ta 5+ ions, where we grew oxides such as RENdO 4 and RETaO 4 as high quality single crystals. These compounds were targeted as potential hosts for luminescent and scintillation materials, particularly given that they are among the densest oxide hosts and hence have good potential as absorbers for high energy radiation like X-rays and gamma rays. We also isolated a range of unusual new rare earth tantalates with very complex structures. indicating that the chemistry is very sensitive to conditions. We performed some fairly comprehensive examinations of the solid-state chemistry of rare earth ions with various tetravalent metal ions especially Si 4+ , Ge 4+ , Sn 4+ and Ti 4+ . Given the potential role of rare earth silicates in immobilizing radioactive waste elements in long-term storage, and the similarity of our hydrothermal fluids with known geological conditions, this chemistry continues to be relevant. We prepared an extensive series of new lanthanide germanates (e.g. RE 13 Ge 6 O 31 (OH), BaRE 10 (GeO 4 ) 4 O 8 ). and found that there is there is almost no overlap between the chemistry of the rare earth silicates. Stannic oxide (SnO 2 ) is much more refractory and requires higher temperatures and of mineralizer concentrations. One significant result is the growth of RE 2 Sn 2 O 7 pyrochlore single crystals. These are of interest because the rare earth stannate pyrochlores are known to display a wide range of magnetic frustration such as spin ice behavior. We grew high quality single crystals of rare earth germanate and stannate pyrochlores and this led to a collaboration with Professor Kate Ross at Colorado State. Preliminary measurements, indicate that the crystals contain no detectable defects or site disorder. Initial neutron diffraction on single crystals was performed at Oak Ridge, and more detailed experiments involving the Ross group are underway at both NIST and ORNL. This particular chemistry has turned out to be the most potentially significant work on this project and the collaborative effort with Prof. Ross is the topic of a DoE renewal project on quantum materials. Our initial foray into the hydrothermal chemistry of rare earth titanates has also been very promising and a range of cubic and polar ferroic phases of the light rare earths RE 2 Ti 2 O 7 (RE = La - Pr) in the P2 1 phase. We also discovered an interesting new phase Ce 2 Ti 4 O 11 that can have implications in heavy metal immobilization and storage. along with a series of new rare earth titanates (La 5 Ti 4 O 15 (OH) Sm 3 TiO 5 (OH) 3 and Lu 5 Ti 2 O 11 (OH) with exceptionally complex structures. One interesting sidelight has been high temperature hydrothermal chemistry terbium, including the growth of large crystals of TbO(OH). This is not a new compound but it is the first time it has been grown as large single crystals. The Tb atom density is almost as high as that in Tb 2 O 3 and has a very high Verdet constant (ca. 70), making it a very attractive candidate as a Faraday rotator. Unfortunately it is not in a cubic structure but he material is hard, stable, pure and inexpensive, so should still be an attractive Faraday oscillator. We recently received a patent on this material. We also synthesized K 2 Tb(Ge 2 O 7 ) containing stable octahedral Tb 4+ ions, which appears to be the first example of a well-characterized Tb 4+ complex. Given that Tb 4+ has been proposed as a benign surrogate for more treacherous tetravalent ions such as Cf 4+ and Bk 4+ , we think that Tb 4+ silicates can be a particularly useful study for actinide immobilization and related work. We also began reaction studies with rare earths and both ReO 2 and RuO 2 . These resulted in large single crystals of species like RE 5 Ru 2 O 12 , RE 4 Re 2 O 11 , REReO 4 and RE 2 ReO 5 . Several of these samples have already been sent to ORNL for magnetic and neutron diffraction studies.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Numerical simulation of terahertz radiation by laser-driven plasma dipole oscillation

Laser-plasma interactions have many theoretical and technological applications. One is the use of coherently accelerated electrons to provide novel sources of THz radiation. Our research focuses on simulating the cross/self-interactions between two high intensity, ultra-short, counter propagating and detuned laser pulses and an initial neutral target for controlled ionization. Unlike our previous studies of laser-matter interaction over preformed plasma, we explore the injection and collision of laser pulses to induce background plasma driven by the self-guided laser wakefield mechanism, which is then used to perturb the plasma resulting in induced dipole oscillations leading to transverse radiation. Inducing a cylindrical spatial plasma column within the laser beam radius regime provides a stable, spatially localized plasma channel. The emitted radiation from the plasma dipole oscillation (PDO) will not be affected by surrounding plasma absorption, resulting in effective radiation distribution. Finally, results include 3D EM-PIC simulations and a comparison of the self- ionizing plasma against the preformed plasma to assess the efficiency of the mechanisms.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Duration of O 2 Exposure Determines Dominance of Fe II vs CH 4 Production in Tropical Forest Soils

Temporal fluctuations in redox conditions influence the availability of Fe III and greenhouse gas emissions in humid upland soils. However, the impact of fluctuation duration on biogeochemical processes remains unclear. We hypothesized that rates of Fe III reduction and CH 4 production are sensitive to the duration of soil oxygenation. To test this, surface soil from the Luquillo Forest, Puerto Rico, was subjected to fluctuating redox conditions with an anoxic interval of 6 days followed by oxic intervals of either 8, 24, or 72 h. Shorter oxic intervals enhanced Fe reduction, while longer oxic intervals enhanced CH 4 emissions. As O 2 exposure decreased from 72 to 8 h, Fe reduction rates increased from 0.12 ± 0.02 to 0.26 ± 0.05 mmol kg –1 h –1 , whereas cumulative CH 4 decreased from 44.0 ± 4.7 to 12.7 ± 4.6 μmol kg –1 . 13 C-amino acid spikes were preferentially incorporated into the DNA of iron reducers (Anaeromyxobacter sp.) in the shorter oxic treatment (8 h vs 24 h), suggesting that Fe reducers are less inhibited by shorter periods of oxidation. Conversely, longer oxygen pulses appear to suppress Fe reducers more than methanogens, leading to increased CH 4 emissions. These findings highlight the role of the redox oscillation length in modulating biogeochemical processes and greenhouse gas emissions in soils.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Strain-insensitive intrinsically stretchable transistors and circuits

An all-elastomer strain engineering approach, which uses patterned elastomer layers with tunable stiffnesses, can be used to create intrinsically stretchable transistor arrays with a device density of 340 transistors cm -2 and strain insensitivity of less than 5% performance variation when stretched to 100% strain. Intrinsically stretchable electronics can form intimate interfaces with the human body, creating devices that could be used to monitor physiological signals without constraining movement. However, mechanical strain invariably leads to the degradation of the electronic properties of the devices. Here we show that strain-insensitive intrinsically stretchable transistor arrays can be created using an all-elastomer strain engineering approach, in which the patterned elastomer layers with tunable stiffnesses are incorporated into the transistor structure. By varying the cross-linking density of the elastomers, areas of increased local stiffness are introduced, reducing strain on the active regions of the devices. This approach can be readily incorporated into existing fabrication processes, and we use it to create arrays with a device density of 340 transistors cm -2 and a strain insensitivity of less than 5% performance variation when stretched to 100% strain. We also show that it can be used to fabricate strain-insensitive circuit elements, including NOR gates, ring oscillators and high-gain amplifiers for the stable monitoring of electrophysiological signals.

42 ENGINEERING↗
DOE OSTI2025

Valley-Mediated Singlet- and Triplet-Polaron Interactions and Quantum Dynamics in a Doped WSe 2 Monolayer

In doped transition metal dichalcogenides, optically created excitons (bound electron-hole pairs) can strongly interact with a Fermi sea of electrons to form Fermi polaron quasiparticles. When there are two distinct Fermi seas, as is the case in WSe 2 , there are two flavors of lowest-energy (attractive) polarons—singlet and triplet—where the exciton is coupled to the Fermi sea in the same or opposite valley, respectively. Using two-dimensional coherent electronic spectroscopy, we analyze how their quantum decoherence evolves with doping density and determine the condition under which stable Fermi polarons form. Because of the large oscillator strength associated with these resonances, intrinsic quantum dynamics of polarons as well as valley coherence between coupled singlet- and triplet polarons occur on subpicosecond timescales. Surprisingly, we find that a dark-to-bright state conversion process leads to a particularly long-lived singlet polaron valley polarization, persisting up to 200–800 ps. Valley coherence between the singlet- and triplet polaron is correlated with their energy fluctuations. In conclusion, our finding provides valuable guidance for the electrical and optical control of spin and valley indexes in atomically thin semiconductors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Stabilize High-IBR Power Systems with Grid-Forming Inverters

The increasing share of IBRs challenges the stable operation of our power systems. Leveraging an oscillation event in Kaua`i Island as an example, we showcase the root cause of this event and highlight the potential of grid-forming inverters to reduce the stability risks in high-IBR power systems.

grid-forming inverters↗
DOE OSTI2020

Oscillation-like diffusion of two-dimensional liquid dusty plasmas on one-dimensional periodic substrates with varied widths

The long-time diffusion of two-dimensional dusty plasmas on a one-dimensional periodic substrate with varying widths is investigated using Langevin dynamical simulations. When the substrate is narrow and the dust particles form a single row, the diffusion is the smallest in both directions. In this study, we find that as the substrate width gradually increases to twice its initial value, the long-time diffusion of the two-dimensional dusty plasmas first increases, then decreases, and finally increases again, giving an oscillation-like diffusion with varying substrate widths. When the width increases to a specific value, the dust particles within each potential well arrange themselves in a stable zigzag pattern, greatly reducing the diffusion and leading to the observed oscillation in the diffusion with the increase in the width. In addition, the long-time oscillation-like diffusion is consistent with the number of dust particles that are hopping across the potential wells of the substrate.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

FY22 Progress on Computational Modeling of the Water-Based NSTF

This report summarizes the system-level modeling effort by Argonne National Laboratory (Argonne) of the Natural convection Shutdown heat removal Test Facility (NSTF) in FY22. As an extension of the effort from FY21, this year’s work focuses primarily on the two-phase modeling of the NSTF using RELAP5-3D, particularly with the inclusion of the cavity model. The results from simulations were used to compare against experimental data for benchmarking purposes of the RELAP5 deck. Additionally, RELAP5 was used as a predictive tool to guide planned test operations and identify expected system behaviors. In the first part of this report, details are provided of the cavity omitted model where heat flux is applied directly as a boundary condition to the risers. The general trend predicted by the RELAP5 model matches that from the experimental data when a single-phase natural circulation flow is first established, followed by an oscillatory two-phase period and finally a stable two-phase flow. However, the onset of oscillations is predicted early by the model due to the smaller thermal mixing region in the tank. However, by expanding the simulated thermal mixing region in the tank, the onset of oscillations predicted by the model is able to match that from the experiment. These oscillations where studied in depth and are deduced to be flashing-induced instability. The model was then modified to simulate an accident scenario case where a representative heat load based on the full-scale Framatome’s 625 MW t SC-HTGR was applied directly to the riser channels. The simulated initial and boundary conditions were identical to those performed experimentally, facilitating direct comparisons between the predicted and experimental results. It was determined that the results showed some discrepancies remain, likely due to the overprediction of vapor generation rate by the computer model. In the second part of this report, the cavity model is re-introduced where it is observed that the RELAP5 prediction is now able to capture the major trends of the observed flow commonly observed during two-phase conditions. However, the onset of oscillations is once again predicted early by the model, possibly caused by the underprediction of heat loss from the heater and cavity. This is likely due to the omission of support structures in the cavity that can act as additional pathways for heat to escape to the environment. To overcome the underprediction of heat loss, part of the insulation surrounding the cavity side panels and the back of the heaters are removed to allow heat to escape directly to the environment, which then improves the RELAP5 prediction. Parametric studies are also performed to investigate the effects of heater power, tank inventory level, and tank gas space pressure on flow behaviors, also in direct comparison to conditions tested experimentally. User option-61 in the RELAP5-3D input deck, which changes the heat transfer coefficient correlations used for calculating the vapor generation, is also investigated where it is found that by enabling the option, the overall duration of oscillations is increased and matches that from the experiment better. The RELAP5 model is further benchmarked with a header inlet- throttling case where it is observed that the prediction from the model fails to capture some major features observed in the experiment. By using a modified loss coefficient curve for the valve, the accuracy of the prediction is improved where most of the major features observed in the experiment are predicted by the model. Lastly, the model is benchmarked with an inventory depletion scenario where it is observed that despite the modeling limitation of RELAP5, the prediction shows good agreement with the experimental data where major trends and features are captured by the model. Future work will see continued development of the current RELAP5-3D input deck of the NSTF to both improve the accuracy of the model’s predictive capability and continuing serving the experimental program. The mutually beneficial relationship between analysis and experimental efforts has become integral to the parent NSTF program, and the greater objective to fully understand and accurately predict the heat removal performance of a full scale RCCS concept.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗