• Journal of Mechanical Science and Technology
• /
• v.17 no.1
• /
• pp.129-138
• /
• 2003

A 1 : 20 laboratory scale test rig of a 200 MW tangentially fired boiler is built up with completely simulated structures such as platen heaters and burners. Iso-thermal turbulent flow in the boiler is mapped by 3-D PDA (Particle Dynamic Analyzer). The 3-D numerical models for the same case are proposed based on the solution of к-$\varepsilon$ model closed RANS (Reynolds time-Averaged Navier-Stokes) equations, which are written in the framework of general coordinates and discretized in the corresponding body-fitted meshes. Not only are the grid lines arranged to fit the inner/outer boundaries. but also to align with the streamlines to the best possibility in order to reduce the NDE (numerical diffusion errors). Extensive comparisons of profiles of mean velocities are carried out between experiment and calculation. Predicted velocities in burner region were quantitatively similar with measured ones, while those in other area have same tendency with experimental counterpart.

• Nuclear Engineering and Technology
• /
• v.48 no.4
• /
• pp.915-924
• /
• 2016

After a severe accident to the nuclear reactor, the in-vessel retention strategy is a key way to prevent the leakage of radioactive material. Nanofluid is a steady suspension used to improve heat-transfer characteristics of working fluids, formed by adding solid particles with diameters below 100nm to the base fluids, and its thermal physical properties and heat-transfer characteristics are much different from the conventional working fluids. Thus, nanofluids with appropriate nanoparticle type and volume concentration can enhance the heat-transfer process. In this study, the moving particle semi-implicit method-meshless advection using flow-directional local grid method is used to simulate the bubble growth, departure, and sliding on the downward-facing heating surface in pure water and nanofluid (1.0 vol.% $Al_2O_3/H_2O$) flow boiling processes; additionally, the bubble critical departure angle and sliding characteristics and their influence are also investigated. The results indicate that the bubble in nanofluid departs from the heating surface more easily and the critical departure inclined angle of nanofluid is greater than that of pure water. In addition, the influence of nanofluid on bubble sliding is not significant compared with pure water.

• Structural Engineering and Mechanics
• /
• v.41 no.1
• /
• pp.25-41
• /
• 2012

Kriging surrogate model provides explicit functions to represent the relationships between the inputs and outputs of a linear or nonlinear system, which is a desirable advantage for response estimation and parameter identification in structural design and model updating problem. However, little research has been carried out in applying Kriging model to crack identification. In this work, a scheme for crack identification based on a Kriging surrogate model is proposed. A modified rectangular grid (MRG) is introduced to move some sample points lying on the boundary into the internal design region, which will provide more useful information for the construction of Kriging model. The initial Kriging model is then constructed by samples of varying crack parameters (locations and sizes) and their corresponding modal frequencies. For identifying crack parameters, a robust stochastic particle swarm optimization (SPSO) algorithm is used to find the global optimal solution beyond the constructed Kriging model. To improve the accuracy of surrogate model, the finite element (FE) analysis soft ANSYS is employed to deal with the re-meshing problem during surrogate model updating. Specially, a simple method for crack number identification is proposed by finding the maximum probability factor. Finally, numerical simulations and experimental research are performed to assess the effectiveness and noise immunity of this proposed scheme.

• Nuclear Engineering and Technology
• /
• v.55 no.12
• /
• pp.4395-4407
• /
• 2023

Bubble columns are widely encountered in several industries, especially in the field of nuclear safety. The Eulerian-Eulerian and the Eulerian-Lagrangian methods are commonly used to investigate bubble columns. Eulerian approaches require additional tasks such as strict volume conservation at the interface and a predefined well-structured grid. In contrast, the Lagrangian approach can be easily implemented. Hence, we introduce a fully Lagrangian approach for the simulation of bubble columns using the discrete bubble model (DBM) and moving particle semi-implicit (MPS) methods. Additionally, we propose a rigorous method to estimate the volume fraction accurately, and verified it through experimental data and analytical results. The MPS method was compared with the experimental data of Dambreak. The DBM was verified by analyzing the terminal velocity of a single bubble for each bubble size. It agreed with the analytical results for each of the four drag correlations. Additionally, the improved method for calculating the volume fraction showed agreement with the Ergun equation for the pressure drop in a packed bed. The implemented MPS-DBM was used to simulate the bubble column, and the results were compared with the experimental results. We demonstrated that the MPS-DBM was in quantitative agreement with the experimental data.

• Structural Engineering and Mechanics
• /
• v.91 no.1
• /
• pp.103-121
• /
• 2024

The smoothed particle hydrodynamics (SPH) method is a numerical technique used in dynamic analysis to simulate the fluid-like behavior of materials under extreme conditions, such as those encountered in explosions or high velocity impacts. In SPH, fluid or solid materials are discretized into particles. These particles interact with each other based on certain smoothing kernels, allowing the simulation of fluid flows and predict the response of solid materials to shock waves, like deformation, cracking or failure. One of the main advantages of SPH is its ability to simulate these phenomena without a fixed grid, making it particularly suitable for analyzing complex geometries. In this study, the structural damage to a masonry arch bridge subjected to blast loading was investigated. A high-fidelity micro-model was created and the explosives were modeled using the SPH approach. The Johnson-Holmquist II damage model and the Mohr-Coulomb material model were considered to evaluate the masonry and backfill properties. Consistent with the principles of the JH-II model, the authors developed a VUMAT code. The explosive charges (50 kg, 168 kg, 425 kg and 1000 kg) were placed in close proximity to the deck and pier of a bridge. The results showed that the 50 kg charges, which could have been placed near the pier by a terrorist, had only a limited effect on the piers. Instead, this charge caused a vertical displacement of the deck due to the confinement effect. Conversely, a 1000 kg TNT charge placed 100 cm above the deck caused significant damage to the bridge.

• Korean Journal of Soil Science and Fertilizer
• /
• v.22 no.3
• /
• pp.163-172
• /
• 1989

In order to study spatial variability of soil physical properties accross the land-partitioned boundary on Hwadong silt clay loam soil (Fine clayey, mixed, mesic family of Aquic Hapludalfs), all measured data were analyzed by means of kriging, fractile diagram, smooth frequency distribution, and autocorrelation. Sampling for soil particle size distribution analysis was made at 225 intersections of $15{\times}15$ grid with 10m interval. Field capacity, bulk density and saturated hydraulic conductivity were measured in situ at 594 intersections of $33{\times}18$ grid with 2.5m interval in only $6,000m^2$ reselected from $22,500m^2$ of sampling area for particle size distribution analysis. Sampled or measured soil depths were 0 to 10cm 25 to 35cm and 50 to 60cm at each intersections. The results are summarized as follows: 1. The coefficient of variance (CV) of various physical properties ranges from 4.8 to 128.8%. Saturated hydraulic conductivity is classified into the high variation group with CV greater than 100%, while the low variation group with CV smaller than 10% consists of bulk density. Other properties belong to the medium variation group with CV between 10 and 100%. 2. The appropriate number of soil samples for the determination of various physical properties with error smaller than 10% are calculated as one for bulk density, six for field moisture capacity, 16 for silt, 19 for clay, 69 for sand and 686 for saturated hydraulic conductivity. 3. Smooth frequency distribution and fractile diagram show that saturated hydraulic conductivity is in lognormal distribution while other physical properties are in normal distribution. 4. Serial correlation analysis reveals that the soil physical properties have spatial dependence between two nearest neighbouring grid points. Autocorrelation analysis of physical properties measured between the serial grid points in the direction of south to north following section boundary shows that the zone of influence showing stationarity ranges from 7.5 to 40m. In the direction of east to west across section boundary, the autocorrelogram of many physical properties shows peaks with the periodic interval of 30m, which are similar to the partitioned land width. This reveals that the land-partitioned boundary causes soil variability.

• Journal of the Korean Society of Safety
• /
• v.31 no.3
• /
• pp.22-27
• /
• 2016

The use of electromagnetic wave has been interested in various energy industry because it enhances a flame stability and provides higher safety environments. However it might increase the pollutant emissions such as NOx and soot, and have harmful influence on human and environments. Therefore, it is very important to understand interaction mechanism between flame and electromagnetic wave from environmental point of view. In this study, an experiment was performed with jet diffusion flames induced by electromagnetic wave. Microwave was used as representative electromagnetic wave and a flickering flame was introduced to simulate the more similar combustion condition to industry. The results show that the induced microwave enhances the flame stability and blowout limit. The unstable lifted flickering flames under low fuel/oxidizer velocity is changed to stable attached flames or lift-off flames when microwave applied to the flames, which results from the abundance of radical pool. However, NOx emission was increased monotonically with increasing the microwave power as microwave power increased up to 1.0 kW. The effects might be attributed to the heating of combustion field and thermal NOx mechanism will be prevailed. Soot particle was examined at the post flame region by TEM grid. The morphology of soot particle sampled in the microwave induced flames was similar to the incipient soot that is not agglomerated and contain a lots of liquid phase hydrocarbon such as PAH, which soot particle formed near reaction zone is oxidized on the extended yellow flame region and hence only unburned young particles are emitted on the post flame region.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.26 no.5B
• /
• pp.557-562
• /
• 2006

In order to predict the dispersion of suspended sediment arising from dredging operation in port and navigation channel, a hybrid model for dispersive transport of turbidity plume was developed using Lee's(1998) hybrid method. Using hybrid modeling scheme advection-diffusion equation was solved by the forward particle-tracking method for advection process and by the fixed Eulerian grid method for diffusion process. To examine numerical model simulation in accuracy, the simulated results for 1-D, 2-D, and 3-D cases were compared with the analytical solutions including Kuo, et al's (1985) 3-D mathematical model. The model results were in a good agreement with the analytical solutions and mathematical model for the dispersion of turbidity plume.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.19 no.1
• /
• pp.81-96
• /
• 2007

Nonlinear shoaling characteristics over surf zone are numerically investigated based on spatially averaged NavierStokes equation. We also test the validity of gradient model for turbulent stresses due to wave breaking using the data acquainted during SUPERTANK LABORATORY DATA COLLECTION PROJECT(Krauss et al., 1992). It turns out that the characteristics length scale of breaking induced current is not negligible, which firmly stands against ever popular gradient model, ${\kappa}-{\varepsilon}$ model, but favors Large Eddy Simulation with finer grid. Based on these observations, we model the residual stress of spatially averaged NavierStokes equation after Lagrangian Dynamic Smagorinsky(Meneveau et al., 1996). We numerically integrate newly proposed wave equations using SPH with Gaussian kernel function. Severely deformed water surface profile, free falling water particle, queuing splash after landing of water particle on the free surface and wave finger due to structured vortex on rear side of wave crest(Narayanaswamy and Dalrymple, 2002) are successfully duplicated in the numerical simulation of wave propagation over uniform slope beach, which so far have been regarded very difficult features to mimic in the computational fluid mechanics.

• International Journal of Contents
• /
• v.7 no.4
• /
• pp.1-5
• /
• 2011

Physics based fluid animation schemes need large computation cost due to tremendous degree of freedom. Many researchers tried to reduce the cost for solving the large linear system that is involved in grid-based schemes. GPU based algorithms and advanced numerical analysis methods are used to efficiently solve the system. Other groups studied local operation methods such as SPH (Smoothed Particle Hydrodynamics) and LBM (Lattice Boltzmann Method) for enhancing the efficiency. Our method investigates this efficiency problem thoroughly, and suggests novel paradigm in fluid animation field. Rather than physics based simulation, we propose a robust boundary handling technique for procedural fluid animation. Our method can be applied to arbitrary shaped objects and potential fields. Since only local operations are involved in our method, parallel computing can be easily implemented.