A multi-scale analysis of water-cooled reactor thermalhydraulics can be used to take advantage of increased computer power and improved simulation tools, including Direct Numerical Simulation (DNS), Computational Fluid Dynamics (CFD) (in both open and porous mediums), and system thermalhydraulic codes. This paper presents a general strategy for this procedure for various thermalhydraulic scales. A short state of the art is given for each scale, and the role of the scale in the overall multi-scale analysis process is defined. System thermalhydraulic codes will remain a privileged tool for many investigations related to safety. CFD in porous medium is already being frequently used for core thermalhydraulics, either in 3D modules of system codes or in component codes. CFD in open medium allows zooming on some reactor components in specific situations, and may be coupled to the system and component scales. Various modeling approaches exist in the domain from DNS to CFD which may be used to improve the understanding of flow processes, and as a basis for developing more physically based models for macroscopic tools. A few examples are given to illustrate the multi-scale approach. Perspectives for the future are drawn from the present state of the art and directions for future research and development are given.

In this paper we describe current activities on the project Multi-Scale Modeling and Analysis of convective boiling (MSMA), conducted jointly by the Paul Scherrer Institute (PSI) and the Swiss Nuclear Utilities (Swissnuclear). The long-term aim of the MSMA project is to formulate improved closure laws for Computational Fluid Dynamics (CFD) simulations for prediction of convective boiling and eventually of the Critical Heat Flux (CHF). As boiling is controlled by the competition of numerous phenomena at various length and time scales, a multi-scale approach is employed to tackle the problem at different scales. In the MSMA project, the scales on which we focus range from the CFD scale (macro-scale), bubble size scale (meso-scale), liquid micro-layer and triple interline scale (micro-scale), and molecular scale (nano-scale). The current focus of the project is on micro- and meso-scales modeling. The numerical framework comprises a highly efficient, parallel DNS solver, the PSI-BOIL code. The code has incorporated an Immersed Boundary Method (IBM) to tackle complex geometries. For simulation of meso-scales (bubbles), we use the Constrained Interpolation Profile method: Conservative Semi-Lagrangian $2^{nd}$ order (CIP-CSL2). The phase change is described either by applying conventional jump conditions at the interface, or by using the Phase Field (PF) approach. In this work, we present selected results for flows in complex geometry using the IBM, selected bubbly flow simulations using the CIP-CSL2 method and results for phase change using the PF approach. In the subsequent stage of the project, the importance of effects of nano-scale processes on the global boiling heat transfer will be evaluated. To validate the models, more experimental information will be needed in the future, so it is expected that the MSMA project will become the seed for a long-term, combined theoretical and experimental program.

A thermal-hydraulic code, named CUPID, has been being developed for the realistic analysis of transient two-phase flows in nuclear reactor components. The CUPID code development was motivated from very practical needs, including the analyses of a downcomer boiling, a two-phase flow mixing in a pool, and a two-phase flow in a direct vessel injection system. The CUPID code adopts a two-fluid, three-field model for two-phase flows, and the governing equations are solved over unstructured grids with a semi-implicit two-step method. This paper presents an overview of the CUPID code development and assessment strategy. It also presents the code couplings with a system code, MARS, and, a three-dimensional reactor kinetics code, MASTER.

Lead-alloys are very attractive nuclear coolants due to their thermo-hydraulic, chemical, and neutronic properties. By utilizing the HELIOS (Heavy Eutectic liquid metal Loop for Integral test of Operability and Safety of PEACER$^2$) facility, a thermal hydraulic benchmarking study has been conducted for the prediction of pressure loss in lead-alloy cooled advanced nuclear energy systems (LACANES). The loop has several complex components that cannot be readily characterized with available pressure loss correlations. Among these components is the core, composed of a vessel, a barrel, heaters separated by complex spacers, and the plenum. Due to the complex shape of the core, its pressure loss is comparable to that of the rest of the loop. Detailed CFD simulations employing different CFD codes are used to determine the pressure loss, and it is found that the spacers contribute to nearly 90 percent of the total pressure loss. In the system codes, spacers are usually accounted for; however, due to the lack of correlations for the exact spacer geometry, the accuracy of models relies strongly on assumptions used for modeling spacers. CFD can be used to determine an appropriate correlation. However, application of CFD also requires careful choice of turbulence models and numerical meshes, which are selected based on extensive experience with liquid metal flow simulations for the KALLA lab. In this paper consistent results of CFX and Star-CD are obtained and compared to measured data. Measured data of the pressure loss of the core are obtained with a differential pressure transducer located between the core inlet and outlet at a flow rate of 13.57kg/s.

To develop the manufacturing methods for the blanket first wall (FW) of the International Thermonuclear Experimental Reactor (ITER) and to verify the integrity of the joint, Be/Cu mockups were fabricated and tested at the KoHLT-1 (Korea Heat Load Test facility), a graphite heater facility located at the Korea Atomic Energy Research Institute (KAERI). Since Be and Cu joining is the focus of the present study, the fabricated mockups had a CuCrZr heat sink joined with three Be tiles as an armor material, unlike the original ITER blanket FW, which has a stainless steel structure and coolant tubes. Hot isostatic pressing (HIP) was carried out at $580^{\circ}C$ and 100 MPa for 2 hours as the method for Be/Cu joining. Three interlayers, namely, $1{\mu}mCr/10{\mu}mCu$, $1{\mu}mTi/0.5{\mu}mCr/10{\mu}mCu$, and $5{\mu}mTi/10{\mu}mCu$ were applied as a coating to the Be tiles by a physical vapor deposition (PVD) method. A shear test was performed with the specimens, which were fabricated by the same methods as those used to fabricate the mockups. The average values were 125 MPa to 180 MPa, and the samples with the $1{\mu}mCr/10{\mu}mCu$ interlayer showed the lowest value. No defect or delamination was found in the joints of the mockups by the developed ultrasonic test using a flat-type probe with a 10 MHz frequency and a 0.25 inch diameter. High heat flux (HHF) tests were performed at $1.0\;MW/m^2$ heat flux for each mockup using the given conditions, and the results were analyzed by ANSYS-CFX code. For the test criteria, an expected fatigue lifetime about 1,000 cycles was obtained by analysis with ANSYS-mechanical code. Mockups using the interlayers of $1{\mu}mTi/0.5{\mu}mCr/10{\mu}mCu$ and $5{\mu}mTi/10{\mu}mCu$ survived up to 1,100 cycles over the required number of cycles. However, one of the Be tiles in the other two mockups using the $1{\mu}mCr/10{\mu}mCu$ interlayer was detached during the screening test, and others were detached by discharge after 862 cycles. The integrity of the joints using the proposed interlayers was proven by the HHF test, but the other interlayer requires more study before it can be used for the joining of Be to Cu. Moreover, it was confirmed that the measured temperatures agreed well with the analysis temperatures, which were used to estimate the lifetime and that the developed facility showed its capability of the long time operation.

In this paper, a new approach using a pixel-based correction method was developed to fix the non-uniform responses of flat-bed type scanners used for radiochromic film dosimetry. In order to validate the method's performance, two cases were tested: the first consisted of simple dose distributions delivered by a single port; the second was a complicated dose distribution composed of multiple beams. In the case of the simple individual dose condition, ten different doses, from 8.3 cGy to 307.1 cGy, were measured, horizontal profiles were analyzed using the pixel-based correcton method and compared with results measured by an ionization chamber and results corrected using the existing correction method. A complicated inverse pyramid dose distribution was made by piling up four different field shapes, which were measured with GAFCHROMIC$^{(R)}$EBT film and compared with the Monte Carlo calculation; as well as the dose distribution corrected using a conventional method. The results showed that a pixel-based correction method reduced dose difference from the reference measurement down to 1% in the flat dose distribution region or 2 mm in a steep dose gradient region compared to the reference data, which were ionization chamber measurement data for simple cases and the MC computed data for the complicated case, with an exception for very low doses of less than about 10 cGy in the simple case. Therefore, the pixel-based scanner correction method is expected to enhance the accuracy of GAFCHROMIC$^{(R)}$EBT film dosimetry, which is a widely used tool for two-dimensional dosimetry.

The crack-tip stress fields and fracture mechanics assessment parameters for a surface crack, such as the elastic stress intensity factor or the elastic-plastic J-integral, can be affected significantly by the adjacent cracks. Such a crack interaction effect due to multiple cracks can alter the fracture mechanics assessment parameters significantly. There are many factors to be considered, for instance the relative distance between adjacent cracks, the crack shape, and the loading condition, to quantify the crack interaction effect on the fracture mechanics assessment parameters. Thus, the current assessment codes on crack interaction effects (crack combination rules), including ASME Sec. XI, BS7910, British Energy R6 and API 579-1/ASME FFS-1, provide different rules for combining multiple surface cracks into a single surface crack. The present paper investigates crack interaction effects by evaluating the elastic stress intensity factor and the elastic-plastic J-integral of adjacent in-plane surface cracks in a plate through detailed 3-dimensional elastic and elastic-plastic finite element analyses. The effects on the fracture mechanics assessment parameters of the geometric parameters, the relative distance between two cracks, and the crack shape are investigated systematically. As for the loading condition, an axial tension is considered. Based on the finite element results, the acceptability of the crack combination rules provided in the existing guidance was investigated, and the relevant recommendations on a crack interaction for in-plane surface cracks are discussed. The present results can be used to develop more concrete guidance on crack interaction effects for crack shape characterization to evaluate the integrity of defective components.

Sampling based uncertainty analysis was carried out to quantify uncertainty in predictions of best estimate code RELAP5/MOD3.2 for a thermal hydraulic test (10% hot leg break LOCA) performed in the Large Scale Test Facility (LSTF) as a part of an IAEA coordinated research project. The nodalisation of the test facility was qualified for both steady state and transient level by systematically applying the procedures led by uncertainty methodology based on accuracy extrapolation (UMAE); uncertainty analysis was carried out using the Latin hypercube sampling (LHS) method to evaluate uncertainty for ten input parameters. Sixteen output parameters were selected for uncertainty evaluation and uncertainty band between $5^{th}$ and $95^{th}$ percentile of the output parameters were evaluated. It was observed that the uncertainty band for the primary pressure during two phase blowdown is larger than that of the remaining period. Similarly, a larger uncertainty band is observed relating to accumulator injection flow during reflood phase. Importance analysis was also carried out and standard rank regression coefficients were computed to quantify the effect of each individual input parameter on output parameters. It was observed that the break discharge coefficient is the most important uncertain parameter relating to the prediction of all the primary side parameters and that the steam generator (SG) relief pressure setting is the most important parameter in predicting the SG secondary pressure.

Modal testing requires measuring the vibration of many points, for which an accelerometer, a gab sensor and laser vibrometer are generally used. Conventional modal testing requires mounting of these sensors to all measurement points in order to acquire the signals. However, this can be disadvantageous because it requires considerable measurement time and effort when there are many measurement points. In this paper, we propose a method for modal testing using a camera image. A camera can measure the vibration of many points at the same time. However, this task requires that the measurement points be classified frame by frame. While it is possible to classify the measurement points one by one, this also requires much time. Therefore, we try to classify multiple points using pattern recognition. The feasibility of the proposed method is verified by a beam experiment. The experimental results demonstrate that we can obtain good results.