• Wind and Structures
• /
• v.20 no.4
• /
• pp.579-607
• /
• 2015

A Combination Random Flow Generation (CRFG) technique for obtaining the fluctuating inflow boundary conditions for Large Eddy Simulation (LES) is proposed. The CRFG technique was developed by combining the typical RFG technique with a novel calculation of k and ${\varepsilon}$ to estimate the length- and time-scales (l, ${\tau}$) of the target fluctuating turbulence field used as the inflow boundary conditions. Through comparatively analyzing the CRFG technique and other existing numerical/experimental results, the CRFG technique was verified for the generation of turbulent wind velocity fields with prescribed turbulent statistics. Using the turbulent velocity fluctuations generated by the CRFG technique, a series of LESs were conducted to investigate the wind flow around S-, R-, L- and U-shaped building models. As the pressures of the models were also measured in wind tunnel tests, the validity of the LES, and the effectiveness of the inflow boundary generated by the CRFG techniques were evaluated through comparing the simulation results to the wind tunnel measurements. The comparison showed that the LES accurately and reliably simulates the wind-induced pressure distributions on the building surfaces, which indirectly validates the CRFG technique in generating realistic fluctuating wind velocities for use in the LES. In addition to the pressure distribution, the LES results were investigated in terms of wind velocity profiles around the building models to reveal the wind flow dynamics around bluff bodies. The LES results quantitatively showed the decay of the bluff body influence when the flow moves away from the building model.

• The Journal of the Acoustical Society of Korea
• /
• v.16 no.5
• /
• pp.43-49
• /
• 1997

Generally, echo cancelers of a TV conference system or a audio conference system are to carry out a real time processing in the case of the closed room having long reverberation time because the system requires much time to modify filter coefficients to environmental changes. Therefore this paper proposes a new subband adaptive filtering method using polyphase filter banks of MPEG(Moving Picture Experts Group) audio system to solve the problems. This method divides signal spectra of input and output into several frequency bands, and each band is adaptively filtered by using ES-NLMS (Exponential Step-Normalized Least Mean Square) algorithm. The optimal number of subband is determined by computational simulations. According to the results of simulation, ERLE of the subband model is 2dB smaller than general full band, calculation rate's of the subband model is decreased about 88%.

• Journal of the Society of Naval Architects of Korea
• /
• v.32 no.3
• /
• pp.72-82
• /
• 1995

In the present paper, flow characteristics and free surface waves generated by a submerged hydrofoil advancing with an uniform speed are calculated. Using a numerical method based on a MAC(Marker And Cell) method, the Navier-Stokes and the continuity equations are solved to simulate flow fields around the hydrofoil. Computations are carried out in a rectangular grid system in which grids are concentrated near the foil and the free surface to improve numerical accuracies. Viscous flow phenomenas including pressure distributions are computed. Moreover, the influences of submerged depths upon the generated wave profiles and the wave breaking phenomena are also investigated. Experiments are performed at the towing tank of Inha University to measure free surface wave elevations due to the advancing hydrofoil. The computational results are compared with the present and the other available experimental data to show the accuracy of the numerical method developed.

• Journal of the Society of Naval Architects of Korea
• /
• v.50 no.2
• /
• pp.79-87
• /
• 2013

In this paper, an Euler equation solver based on a Cartesian-grid method and non-uniform staggered grid system is applied to predict the ship motion response and added resistance in waves. Water, air, and solid domains are identified by a volume-fraction function for each phase and in each cell. For capturing the interface between air and water, the tangent of hyperbola for interface capturing (THINC) scheme is used with a weighed line interface calculation (WLIC) method. The volume fraction of solid body embedded in a Cartesian-grid system is calculated by a level-set based algorithm, and the body boundary condition is imposed by volume weighted formula. Added resistance is calculated by direct pressure integration on the ship surface. Numerical simulations for a Wigley III hull and an S175 containership in regular waves have been carried out to validate the newly developed code, and the ship motion responses and added resistances are compared with experimental data. For S175 containership, grid convergence test has been conducted to investigate the sensitivity of grid spacing on the motion responses and added resistances.

• Nuclear Engineering and Technology
• /
• v.47 no.4
• /
• pp.472-478
• /
• 2015

The patient dose incurred from diagnostic procedures during advanced radiotherapy has become an important issue. Many researchers in medical physics are using computational simulations to calculate complex parameters in experiments. However, extended computation times make it difficult for personal computers to run the conventional Monte Carlo method to simulate radiological images with high-flux photons such as images produced by computed tomography (CT). To minimize the computation time without degrading imaging quality, we applied a deterministic adaptation to the Monte Carlo calculation and verified its effectiveness by simulating CT image reconstruction for an image evaluation phantom (Catphan; Phantom Laboratory, New York NY, USA) and a human-like voxel phantom (KTMAN-2) (Los Alamos National Laboratory, Los Alamos, NM, USA). For the deterministic adaptation, the relationship between iteration numbers and the simulations was estimated and the option to simulate scattered radiation was evaluated. The processing times of simulations using the adaptive method were at least 500 times faster than those using a conventional statistical process. In addition, compared with the conventional statistical method, the adaptive method provided images that were more similar to the experimental images, which proved that the adaptive method was highly effective for a simulation that requires a large number of iterations-assuming no radiation scattering in the vicinity of detectors minimized artifacts in the reconstructed image.

• Nuclear Engineering and Technology
• /
• v.48 no.3
• /
• pp.765-772
• /
• 2016

In analyzing residual radiation, researchers generally use a two-step Monte Carlo (MC) simulation. The first step (MC1) simulates neutron transport, and the second step (MC2) transports the decay photons emitted from the activated materials. In this process, the stochastic uncertainty estimated by the MC2 appears only as a final result, but it is underestimated because the stochastic error generated in MC1 cannot be directly included in MC2. Hence, estimating the true stochastic uncertainty requires quantifying the propagation degree of the stochastic error in MC1. The brute force technique is a straightforward method to estimate the true uncertainty. However, it is a costly method to obtain reliable results. Another method, called the adjoint-based method, can reduce the computational time needed to evaluate the true uncertainty; however, there are limitations. To address those limitations, we propose a new strategy to estimate uncertainty propagation without any additional calculations in two-step MC simulations. To verify the proposed method, we applied it to activation benchmark problems and compared the results with those of previous methods. The results show that the proposed method increases the applicability and user-friendliness preserving accuracy in quantifying uncertainty propagation. We expect that the proposed strategy will contribute to efficient and accurate two-step MC calculations.

• Economic and Environmental Geology
• /
• v.39 no.4 s.179
• /
• pp.473-483
• /
• 2006

Although application of electrical methods in Korea began with observation of self potentials before World War II, the methods were developed slowly by the beginning of 1980's when a major burst of development activity took place. DC resistivity methods are applied in Korea more to geotechnical problems rather than to environmental ones unlike other developed countries. As with every other branch of technology, the evolving speed of the silicon chip and of streaming data to hard disk has revolutionized data collection and noise reduction processing. The last two decades saw major advances in data collection, processing, and interpretation of electrical data. Development of smooth-model two-dimensional (2D) resistivity inversion is one of the most visible changes to geophysical interpretation of the last 40 years and is now routinely applied to apparent resistivity data. The ability to represent resistivities in section rather than pseudosection view has revolutionized interpretation. Although calculation of sensitivities for general electromagnetic problems require numerous forward modelings, DC resistivity methods can enjoy computational efficiencies if sources and receivers occupy the same position, and previously intractable 3D inversion is now becoming available.

• Journal of the Korean Institute of Gas
• /
• v.19 no.1
• /
• pp.38-44
• /
• 2015

FDS (Fire Dynamics Simulator) is the most widely used computational fluid dynamics software in the fire safety engineering community, and it is applicable to various evaluations of fire growth and its effects. This study made use of a range of outputs from FDS simulation to predict FED (Fractional Effective Concentration) and FEC (Fractional Effective Concentration) levels which are often adopted to evaluate toxicity of fire smoke. As it is not possible to calculate these values directly from outputs of FDS, it was necessary to produce them by means of additional calculation procedures incorporating results of evacuation simulation. In this study, the latest version of FDS, which was recently updated in November 2013, was utilized for the purpose of quantitative comparison with the old version of FDS. As a result, it was found that they make about 10 percent difference on average in predicting FED and FEC levels for the cable fire case study.

• Advances in Computational Design
• /
• v.2 no.3
• /
• pp.123-141
• /
• 2017

Solving a system of linear or non-linear equations is required to analyze any kind of structures. There are many ways to solve a system of equations, and they can be classified as implicit and explicit techniques. The explicit methods eliminate round-off errors and use less memory. The dynamic relaxation method (DR) is one of the powerful and simple explicit processes. The important point is that the DR does not require to store the global stiffness matrix, for which it just uses the residual loads vector. In this paper, a new approach to the DR method is expressed. In this approach, the damping, mass and time steps are similar to those of the traditional method of dynamic relaxation. The difference of this proposed method is focused on the method of calculating the damping. The proposed method is expressed such that the time step is constant, damping is equal to zero except in steps with maximum energy and the concentrated damping can be applied to minimize the energy of system in this step. In this condition, the calculation of damping in all steps is not required. Then the volume of computation is reduced. The DR method for form-finding of membrane structures is employed in this paper. The form-finding of the three plans related to the membrane structures with different loading is considered to investigate the efficiency of the proposed method. The numerical results show that the convergence rate based on the proposed method increases in all cases than other methods.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03b
• /
• pp.256-259
• /
• 2008

Ejector system is a device to transport a low-pressure secondary flow by using a high-pressure primary flow. Ejector system is, in general, composed of a primary nozzle, a mixing section, a casing part for suction of secondary flow and a diffuser. It can induce the secondary flow or affect the secondary chamber pressure by both shear stress and pressure drop which are generated in the primary jet boundary. Ejector system is simple in construction and has no moving parts, so it can not only compress and transport a massive capacity of fluid without trouble, but also has little need for maintenance. Ejectors are widely used in a range of applications such as a turbine-based combined-cycle propulsion system and a high altitude test facility for rocket engine, pressure recovery system, desalination plant and ejector ramjet etc. The primary interest of this study is to set up an applicable model and operating conditions for an ejector in the condition of sonic and subsonic, which can be extended to the hydrogen fuel cell vehicle. Experimental and theoretical investigation on the sonic and subsonic ejectors with a converging-diverging diffuser was carried out. Optimization technique and numerical simulation was adopted for an optimal geometry design and satisfying the required performance at design point of ejector for hydrogen recirculation. Also, some ejectors with a various of nozzle throat and mixing chamber diameter were manufactured precisely and tested for the comparison with the calculation results.