High-resolution Numerical Data Research Articles

The European GO-VIKING project on flow-induced vibrations: overview and current status

Abstract The interaction between cooling fluid and solid structures (rods, tubes) in nuclear power plants may lead to flow-induced vibrations (FIV), causing material fatigue, fretting wear, and eventually loss of component integrity. This can cause further safety issues as well as substantial standstill costs due to longer or unplanned outages. With the growing computational power, the application of modern 3D numerical simulation tools for the accurate prediction of FIV phenomena is rapidly increasing. In 2022, the GO-VIKING (Gathering expertise On Vibration ImpaKt In Nuclear power Generation) project received a grant within the Horizon Europe research and innovation funding program. Sixteen European and two US partners started their collaboration in the field of FIV experiments and analysis. Over four years, the GO-VIKING project investigates FIV phenomena occurring in nuclear reactor cores and steam generators under single- and two-phase flow conditions. The project’s main objectives are to expand the expertise in the field of FIV through generation of new experimental and high-resolution numerical data; development, improvement, and validation of fluid-structure interaction (FSI) methods for FIV evaluation; training stakeholders in the application of these methods; and synthesizing guidelines for the prediction and assessment of FIV phenomena in nuclear reactors. This paper provides an overview of the GO-VIKING objectives, scientific program, as well as of the main scientific achievements in the first project year.

Nucleate boiling simulation using interface tracking method

The development and validation of 3D multiphase computational fluid dynamics (M-CFD) models and physics-informed data-driven modeling require data of high-quality and high-resolution. Considering the difficulties in acquiring the corresponding experimental data in prototypical conditions, two-phase boiling simulations by Interface Tracking Method (ITM) based models can be used to generate high-resolution numerical data in a consistent and relatively economical manner. A boiling model is developed in one of the ITM-based multiphase-flow solvers, named PHASTA, to investigate the nucleate boiling phenomenon. The interaction between bubbles forming at adjacent nucleation sites is investigated with this ITM boiling model. Nucleate pool boiling simulations with multiple nucleation sites are presented in this paper. Theinfluences of site distance, neighboring bubble size and contact angle effect are investigated. The presented boiling simulations are conducted on 3D unstructured computational meshes. These simulation results improve our understanding of the physical mechanisms of the nucleate boiling phenomenon and provide high-resolution numerical data for M-CFD validation and advanced boiling model development.

One-armed spiral instability in neutron star mergers and its detectability in gravitational waves

We study the development and saturation of the $m=1$ one-armed spiral instability in remnants of binary neutron star mergers by means of high-resolution long-term numerical relativity simulations. Our results suggest that this instability is a generic outcome of neutron stars mergers in astrophysically relevant configurations; including both "stiff" and "soft" nuclear equations of state. We find that, once seeded at merger, the $m=1$ mode saturates within $\sim 10\ \mathrm{ms}$ and persists over secular timescales. Gravitational waves emitted by the $m=1$ instability have a peak frequency around $1-2\ \mathrm{kHz}$ and, if detected, could be used to constrain the equation of state of neutron stars. We construct hybrid waveforms spanning the entire Advanced LIGO band by combining our high-resolution numerical data with state-of-the-art effective-one-body waveforms including tidal effects. We use the complete hybrid waveforms to study the detectability of the one-armed spiral instability for both Advanced LIGO and the Einstein Telescope. We conclude that the one-armed spiral instability is not an efficient gravitational wave emitter. Its observation by current generation detectors is unlikely and will require third-generation interferometers.

Evolution of a double-front Rayleigh-Taylor system using a graphics-processing-unit-based high-resolution thermal lattice-Boltzmann model.

We study the turbulent evolution originated from a system subjected to a Rayleigh-Taylor instability with a double density at high resolution in a two-dimensional geometry using a highly optimized thermal lattice-Boltzmann code for GPUs. Our investigation's initial condition, given by the superposition of three layers with three different densities, leads to the development of two Rayleigh-Taylor fronts that expand upward and downward and collide in the middle of the cell. By using high-resolution numerical data we highlight the effects induced by the collision of the two turbulent fronts in the long-time asymptotic regime. We also provide details on the optimized lattice-Boltzmann code that we have run on a cluster of GPUs.

Characterization of flow and mixing in an SMX static mixer

AbstractLaminar flow and mixing of a Newtonian fluid are characterized in a four‐element Koch‐Glitsch SMX static mixer. The computational analysis on a fine, unstructured mesh containing more than 3.5 million tetrahedral elements led to high‐resolution numerical data for velocity and pressure. Computed pressure drops agreed excellently with those in the literature. The flow in this static mixer is essentially independent of the flow rate up to Re = 1, thereby causing the flow characteristics to deviate substantially from those observed when Re < 1. Furthermore, mixing behavior is examined using Lagrangian particle tracking simulations, as well as statistical methods to facilitate comparisons with experimental data. The mixing rate calculated from the simulations agreed closely with experimentally measured values. Both exponential decrease in the variation coefficient as a function of downstream length and the evolution of a self‐similar mixing structure are both evidence of chaotic mixing in the SMX mixer.

Estimating Ω from Galaxy Redshifts: Linear Flow Distortions and Nonlinear Clustering

We propose a method to determine the cosmic mass density Ω from redshift-space distortions induced by large-scale flows in the presence of nonlinear clustering. Nonlinear structures in redshift space, such as fingers of God, can contaminate distortions from linear flows on scales as large as several times the small-scale pairwise velocity dispersion σv. Following Peacock & Dodds, we work in the Fourier domain and propose a model to describe the anisotropy in the redshift-space power spectrum; tests with high-resolution numerical data demonstrate that the model is robust for both mass and biased galaxy halos on translinear scales and above. On the basis of this model, we propose an estimator of the linear growth parameter β = Ω0.6/b, where b measures bias, derived from sampling functions that are tuned to eliminate distortions from nonlinear clustering. The measure is tested on the numerical data and found to recover the true value of β to within ~10%. An analysis of IRAS 1.2 Jy galaxies yields β=0.8+0.4-0.3 at a scale of 1000 km s-1, which is close to optimal given the shot noise and finite size of the survey. This measurement is consistent with dynamical estimates of β derived from both real-space and redshift-space information. The importance of the method presented here is that nonlinear clustering effects are removed to enable linear correlation anisotropy measurements on scales approaching the translinear regime. We discuss implications for analyses of forthcoming optical redshift surveys in which the dispersion is more than a factor of 2 greater than in the IRAS data.

Contour Advection with Surgery: A Technique for Investigating Finescale Structure in Tracer Transport

We present a trajectory technique, contour advection with surgery (CAS), for tracing the evolution of material contours in a specified (including observed) evolving flow. CAS uses the algorithms developed by Dritschel for contour dynamics/surgery to trace the evolution of specified contours. The contours are represented by a series of particles, which are advected by a specified, gridded, wind distribution. The resolution of the contours is preserved by continually adjusting the number of particles, and finescale features are produced that are not present in the input data (and cannot easily be generated using standard trajectory techniques). The reliability, and dependence on the spatial and temporal resolution of the wind field, of the CAS procedure is examined by comparisons with high-resolution numerical data (from contour dynamics calculations and from a general circulation model), and with routine stratospheric analyses. These comparisons show that the large-scale motions dominate the deformation field and that CAS can accurately reproduce small scales from low-resolution wind fields. The CAS technique therefore enables examination of atmospheric tracer transport at previously unattainable resolution.