• Nuclear Engineering and Technology
• /
• 제37권6호
• /
• pp.511-524
• /
• 2005

Thermalhydraulic reactor simulation of tomorrow will require a new generation of codes combining at least three scales, the CFD scale in open medium, the component scale and the system scale. DNS will be used as a support for modelling more macroscopic models. NEPTUNE is such a new generation multi-scale platform developed jointly by CEA-DEN and EDF-R&D and also supported by IRSN and FRAMATOME-ANP. The major steps towards the next generation lie in new physical models and improved numerical methods. This paper presents the advances obtained so far in physical modelling for each scale. Macroscopic models of system and component scales include multi-field modelling, transport of interfacial area, and turbulence modelling. Two-phase CFD or CMFD was first applied to boiling bubbly flow for departure from nucleate boiling investigations and to stratified flow for pressurised thermal shock investigations. The main challenges of the project are presented, some selected results are shown for each scale, and the perspectives for future are also drawn. Direct Numerical Simulation tools with Interface Tracking Techniques are also developed for even smaller scale investigations leading to a better understanding of basic physical processes and allowing the development of closure relations for macroscopic and CFD models.

• Interaction and multiscale mechanics
• /
• 제6권2호
• /
• pp.157-172
• /
• 2013

This paper presented an investigation of macromolecular suspension in a grooved channel by using the dissipative particle dynamics (DPD) with finitely extensible non-linear elastic (FENE) bead spring chains model. Before studying the movement and evolution of macromolecules, the DPD method was first validated by modeling the simple fluid flow in the grooved channel. For both simple fluid flow and macromolecular suspension, the flow fields were analyzed in detail. It is found that the structure of the grooved channel with sudden contraction and expansion strongly affects the velocity distribution. As the width of the channel reduces, the horizontal velocity increases simultaneously. Vortices can also be found at the top and bottom corners behind the contraction section. For macromolecular suspension, the macromolecular chains influence velocity and density distribution rather than the temperature and pressure. Macromolecules tend to drag simple fluid particles, reducing the velocity with density and velocity fluctuations. Particle trajectories and evolution of macromolecular conformation were investigated. The structure of the grooved channel with sudden contraction and expansion significantly influence the evolution of macromolecular conformation, while macromolecules display adaptivity to adjust their own conformation and angle to suit the structure so as to pass the channel smoothly.

• Nuclear Engineering and Technology
• /
• 제47권5호
• /
• pp.523-533
• /
• 2015

Three-dimensional flow phenomena in a wire-wrapped 37-pin fuel assembly mock-up of a Japanese loop-type sodium-cooled fast reactor, Monju, were investigated with a numerical analysis using a general-purpose commercial computational fluid dynamics code, CFX. Complicated and vortical flow phenomena in the wire-wrapped 37-pin fuel assembly were captured by a Reynolds-averaged Navier-Stokes flow simulation using a shear stress transport turbulence model. The main purpose of the current study is to understand the three-dimensional complex flow phenomena in a wire-wrapped fuel assembly to support the license issue for the core design. Computational fluid dynamics results show good agreement with friction factor correlation models. The secondary flow in the corner and edge subchannels is much stronger than that in an interior subchannel. The axial velocity averaged in the corner and edge subchannels is higher than that averaged in the interior subchannels. Three-dimensional multiscale vortex structures start to be formed by an interaction between secondary flows around each wire-wrapped pin. Behavior of the large-scale vortex structures in the corner and edge subchannels is closely related to the relative position between the hexagonal duct wall and the helically wrapped wire spacer. The small-scale vortex is axially developed in the interior subchannels. Furthermore, a driving force on each wire spacer surface is closely related to the relative position between the hexagonal duct wall and the wire spacer.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2021년도 학술발표회
• /
• pp.147-147
• /
• 2021

A partition of water distribution network (WDN) into district metered areas (DMAs) brings the efficiency and efficacy for water network operation and management (O&M), especially in monitoring pressure and leakage. Traditionally, the DMA configurations (i.e., number, shape, and size of DMAs) are permanent and cannot be changed occasionally. This leads to changes in water quality and reduced network redundancy lowering network resilience against abnormal conditions such as water demand variability and mechanical failures. This study proposes a framework to automatically divide a WDN into dynamic DMA configurations, in which the DMA layouts can self-adapt in response to abnormal scenarios. To that aim, a complex graph theory is adopted to sectorize a WDN into multiscale DMA layouts. Then, different failure-based scenarios are investigated on the existing DMA layouts. Here, an optimization-based model is proposed to convert existing DMA layouts into dynamic layouts by considering existing valves and possibly placing new valves. The objective is to minimize the alteration of flow paths (i.e., flow direction and velocity in the pipes) while preserving the hydraulic performance of the network. The proposed method is tested on a real complex WDN for demonstration and validation of the approach.

• Interaction and multiscale mechanics
• /
• 제5권3호
• /
• pp.287-318
• /
• 2012

We apply a partitioned-solution (iterative-staggered) coupling method based on a fixed Eulerian mesh with the level set function to a large-deformation fluid-structure interaction (FSI) problem where a large-deformable thin structure moves in a high-speed flow field, as an airbag does during deployment. This method combines advanced fluid and structure solvers-specifically, the constrained interpolation profile finite element method (CIP-FEM) for fluid Eulerian mesh and large-deformable structural elements for Lagrangian structural mesh. We express the large-deformable interface as a zero isosurface by the level set function, and introduce virtual nodes with level sets and structural normal velocities to generate the level set function according to the large-deformable interfacial geometry and enforce the kinematic condition at the interface. The virtual nodes are located in the direction normal to the structural mesh. It is confirmed that application of the method to unfolded airbag deployment simulation shows the adequacy of the method.

• Interaction and multiscale mechanics
• /
• 제6권2호
• /
• pp.211-235
• /
• 2013

An incompressible smoothed particle hydrodynamics (ISPH) method based on the incremental pressure projection method is developed in this study. The Rayleigh-B$\acute{e}$nard convection in a square enclosure is used as a validation case and the results obtained by the proposed ISPH model are compared to the benchmark solutions. The comparison shows that the established ISPH method has a good performance in terms of accuracy. Subsequently, the proposed ISPH method is employed to simulate natural convection from a heated cylinder in a square enclosure. It shows that the predictions obtained by the ISPH method are in good agreements with the results obtained by previous studies using alternative numerical methods. A rotating and heated cylinder is also considered to study the effect of the rotation on the heat transfer process in the enclosure space. The numerical results show that for a square enclosure at, the addition of kinetic energy in the form of rotation does not enhance the heat transfer process. The method is also applied to simulate forced convection from a circular cylinder in an unbounded uniform flow. In terms of results, it turns out that the proposed ISPH model is capable to simulate heat transfer problems with the complex and moving boundaries.

• 한국입자에어로졸학회지
• /
• 제12권1호
• /
• pp.1-9
• /
• 2016

In this work, we evaluated the characteristics of flow field and uniformity of the nose-only exposure chambers for the inhalation toxicity test. Computational fluid dynamics (CFD) modeling was carried out to demonstrate uniformity of the nose-only exposure chambers. Because it is very important in the inhalation toxicity experiments that test materials are distributed uniformly to each holder of the chamber. The test was done with these 3 types of chamber with different form to develop inhalation toxicity evaluation system, easy-to-operate system among exposure chamber used for evaluating inhalation toxicity of environmental chemical mixtures. Through CFD interpretation, nose-only exposure chamber was made with the selection of the optimal conditions. For its evaluation, one type of fragrance was selected and measured particle size distribution of each port. The gene becoming luminous to green fluorescence was combined with GPT-SPE, a type of tGFP vector, to be inhaled to the mouse. Based on this, luminous intensity was checked. As a result, total particle number concentration of each port had average value of $3.17{\times}10^6{\sharp}/cm^3$ and range of the highest and lowest concentration value was approximately ${\pm}4.8%$. Autopsy of lung tissues of mouse showed that it had clearly better delivery of gene compared to the control group.

• Wind and Structures
• /
• 제37권1호
• /
• pp.25-38
• /
• 2023

Due to the complexity and difficulty in meeting the multiphase flow complexity, similarity, and multiscale characteristics, the mechanism of snow drift is so complicated that the snow deposition prediction is still inaccurate and needs to be far improved. Meanwhile, the validation of prediction methods is also limited due to a lack of field-measured data about snow deposition. To this end, a field measurement activity about snow deposition around a cube with time was carried out, and the snow accumulation process was measured under blowing snow conditions in northwest China. The maximum snow depth, snow profile, and variation in snow depth around the cube were discussed and analyzed. The measured results indicated three stages of snow accumulation around the cube. First, snow is deposited in windward, lateral and leeward regions, and then the snow depth in windward and lateral regions increases. Secondly, when the snow in the windward region reaches its maximum, the downwash flow erodes the snow against the front wall. Meanwhile, snow range and depth in lateral regions have a significant increase. Thirdly, a narrow road in the leeward region is formed with the increase in snow range and depth, which results in higher wind speed and reforming snow deposition there. The field measurement study in this paper not only furthers understanding of the snow accumulation process instead of final deposition under complex conditions but also provides an important benchmark for validating prediction methods.

• 한국대기환경학회지
• /
• 제34권2호
• /
• pp.294-305
• /
• 2018

Quantitative assessment on the impact from North Korean emissions to surface particulate matter(PM) concentration in the Seoul Metropolitan Area (SMA), South Korea is conducted using a 3-dimensional chemistry transport model. Transboundary transport of air pollutants and their precursors are important to understand regional air quality in East Asian countries. As North Korea locates in the middle of main transport pathways of Chinese pollutants, quantifiable estimation of its impact is essential for policy making in South Korean air quality management. In this study, the Community Multiscale Air Quality Modeling System is utilized to simulate regional air quality and its sensitivity, using the Comprehensive Regional Emissions inventory for Atmospheric Transport Experiment 2015 and the Clean Air Policy Support System 2013 emissions inventories for North and South Korea, respectively. Contributions were estimated by a brute force method, perturbing 50% of North and South Korean emissions. Simulations demonstrate that North Korean emissions contribute $3.89{\mu}g/m^3$ of annual surface PM concentrations in the SMA, which accounts 14.7% of the region's average. Impacts are dominant in nitrate and organic carbon (OC) concentrations, attributing almost 40% of SMA OC concentration during January and February. Clear seasonal variations are also found in North Korean emissions contribution to South Korea (and vice versa) due to seasonal characteristics of synoptic weather, especially by the change of seasonal flow patterns.

• Interaction and multiscale mechanics
• /
• 제2권1호
• /
• pp.45-68
• /
• 2009

An aim of the study is to develop an efficient numerical simulation technique that can handle the two-scale analysis of fluid permeation filters fabricated by the partial sintering technique of small spherical ceramics. A solid-fluid mixture homogenization method is introduced to predict the mechanical characters such as rigidity and permeability of the porous ceramic filters from the micro-scale geometry and configuration of partially-sintered particles. An extended finite element (X-FE) discretization technique based on the enriched interpolations of respective characteristic functions at fluid-solid interfaces is proposed for the non-interface-fitted mesh solution of the micro-scale analysis that needs non-slip condition at the interface between solid and fluid phases of the unit cell. The homogenization and localization performances of the proposed method are shown in a typical two-dimensional benchmark problem whose model has a hole in center. Three-dimensional applications to the body-centered cubic (BCC) and face-centered cubic (FCC) unit cell models are also shown in the paper. The 3D application is prepared toward the computer-aided optimal design of ceramic filters. The accuracy and stability of the X-FEM based method are comparable to those of the standard interface-fitted FEM, and are superior to those of the voxel type FEM that is often used in such complex micro geometry cases.