• Journal of computational fluids engineering
• /
• v.12 no.4
• /
• pp.14-19
• /
• 2007

We present a direct numerical simulation technique and some preliminary results of the capillary spreading of a droplet containing particles on the solid substrate. We used the level-set method with the continuous surface stress for description of droplet spreading with interfacial tension and employed the discontinuous Galerkin method for the stabilization of the interface advection equation. The distributed Lagrangian-multipliers method has been combined for the implicit treatment of rigid particles. We investigated the droplet spreading by the capillary force and discussed effects of the presence of particles on the spreading behavior. It has been observed that a particulate drop spreads less than the pure liquid drop. The amount of spread of a particulate drop has been found smaller than that of the liquid with effectively the same viscosity as the particulate drop.

• 한국전산유체공학회:학술대회논문집
• /
• 2005.10a
• /
• pp.11-19
• /
• 2005

First, methods of numerical analysis of gas-particle flows is classified into micro, meso and macro scale approaches based on the concept of multi-scale mechanics. Next, the explanation moves on to discrete particle simulation where motion of individual particles is calculated numerically using the Newtonian equations of motion. The author focuses on the cases where particle-to-particle interaction has significant effects on the phenomena. Concerning the particle-to-particle interaction, two cases are considered: the one is collision-dominated flows and the other is the contact-dominated flows. To treat this interaction mathematically, techniques named DEM(Distinct Element Method) or DSMC (Direct Simulation Monte Carlo) have been developed DEM, which has been developed in the field of soil mechanics, is useful for the contact -dominated flows and DSMC method, developed in molecular gas flows, is for the collision-dominated flows. Combining DEM or DSMC with CFD (computer fluid dynamics), the discrete particle simulation becomes a more practical tool for industrial flows because not only the particle-particle interaction but particle-fluid interaction can be handled. As examples of simulations, various results are shown, such as hopper flows, particle segregation phenomena, particle mixing in a rotating drum, dense phase pneumatic conveying, spouted bed, dense phase fluidized bed, fast circulating fluidized bed and so on.

• Structural Engineering and Mechanics
• /
• v.25 no.3
• /
• pp.347-363
• /
• 2007

A novel methodology, referred to as Auxiliary Domain Method (ADM), allowing for a very efficient solution of nonlinear reliability problems is presented. The target nonlinear failure domain is first populated by samples generated with the help of a Markov Chain. Based on these samples an auxiliary failure domain (AFD), corresponding to an auxiliary reliability problem, is introduced. The criteria for selecting the AFD are discussed. The emphasis in this paper is on the selection of the auxiliary linear failure domain in the case where the original nonlinear reliability problem involves multiple objectives rather than a single objective. Each reliability objective is assumed to correspond to a particular response quantity not exceeding a corresponding threshold. Once the AFD has been specified the method proceeds with a modified subset simulation procedure where the first step involves the direct simulation of samples in the AFD, rather than standard Monte Carlo simulation as required in standard subset simulation. While the method is applicable to general nonlinear reliability problems herein the focus is on the calculation of the probability of failure of nonlinear dynamical systems subjected to Gaussian random excitations. The method is demonstrated through such a numerical example involving two reliability objectives and a very large number of random variables. It is found that ADM is very efficient and offers drastic improvements over standard subset simulation, especially when one deals with low probability failure events.

• Computers and Concrete
• /
• v.24 no.2
• /
• pp.103-115
• /
• 2019

In this study, the strength and stiffness (EI) of wood-concrete composite beams for bridges with T-shaped cross section were evaluated. Two types of connectors were used: connectors bonded with epoxy adhesive and connectors attached to the wood just by pre-drilling (without adhesive). The connectors consisted of common steel bars with a diameter of 12.5 mm. Initially, the strength and stiffness (EI) of the beams were analyzed by bending tests with the load applied at the third point of the beam. Subsequently, the composite beams were evaluated by numerical simulation using ANSYS software with focus on the connection system. To make the composite beams, Eucalyptus citriodora wood and medium strength concrete were used. The slip modulus K and the ultimate strength values of each type of connector were obtained by direct shear tests performed on composite specimens. The results showed that the connector glued with epoxy adhesive resulted in better strength and stiffness (EI) for the composite beams when compared to the connector fixed by pre-drilling. The differences observed were up to 10%. The strength and stiffness (EI) values obtained analytically by $M{\ddot{o}}hler^{\prime}$ model were lower than the values obtained experimentally from the bending tests, and the differences were up to 25%. The numerical simulations allowed, with reasonable approximation, the evaluation of stress distributions in the composite beams tested experimentally.

• The Bulletin of The Korean Astronomical Society
• /
• v.40 no.2
• /
• pp.45.2-45.2
• /
• 2015

It is obvious that there are plentiful of dust particles in the interplanetary spaces of the Solar System (IDPs), based on micrometeor craters, zodiacal light and direct measurements on the spacecraft. Because of photon drag and planetary perturbations, these particles are continuously falling to the Sun or planets, therefore continuous source of the IDPs are required. We studied the fractional contribution of each type of solar system objects to the IDPs complex through the optical properties of the potential dust sources and the zodiacal light. We found that more than 90% of the IDPs are originated from cometary nuclei. This result is discussed through the comparison with the dynamic simulation, micrometeors mineralogy and near-infrared spectrum of the zodiacal light. In addition, we introduce our new project on the numerical simulation for the dust particles ejected from the cometary nuclei, to verify the conclusion of dominant cometary contribution and its detailed consequences.

• Proceedings of the KIEE Conference
• /
• 2006.04a
• /
• pp.252-254
• /
• 2006

The recent improvement in the performance of digital processor, the application of control technology, which used in the HVDC(High Voltage Direct Current) system with the digital processors, has increased. Having this research development as the basis, this paper presents an achievement of progression by tuning the parameter of PI controller based on Genetic Algorithms(GAs) and by controlling with PI controller with a developed simulator by applying the Matrix operating function, voltage source switching element, modified nodal analysis which can include transformer and the backward Euler which does not create the problem of numerical oscillation. As a result, I expect this development in the simulator HVDC System to bring more application in the field of control technology research with an expanded practicality.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.7
• /
• pp.984-994
• /
• 2003

The suitability of high-order accurate, centered and upwind-biased compact difference schemes for large eddy simulation is evaluated by a spectral, static error analysis. To investigate the effect of numerical dissipation on LES solutions, power spectra of discretization errors are evaluated for isotropic turbulence models in both continuous and discrete wavevector spaces. Contrary to the common belief, the aliasing errors from upwind-biased schemes are larger than those from comparable non-dissipative schemes. However, this result is the direct consequence of the definition of the power spectral density of the aliasing error, which poses the limitation of the static error analysis for upwind schemes.

• Proceedings of the KSME Conference
• /
• 2008.11a
• /
• pp.973-978
• /
• 2008

This study concerned on the critical speed due to hunting and snake motion train to ensure the stability. First, the critical speed was calculated by using a numerical model, and calculated the critical speed of the vehicle through the simulation with the use of ADAMS/RAII. Also, the snake motion was confirmed through a modal analysis and running simulation. The calculated results, show that the rail irregularity becomes the influential factors of the stability since it is the direct source of excitation of the vehicle.

• Architectural research
• /
• v.8 no.2
• /
• pp.27-36
• /
• 2006

Structural uncertainties are generally modeled using probabilistic approaches in order to quantify uncertainties in behaviors of structures. This uncertainty results from the uncertainties of structural parameters. Monte Carlo methods have been usually carried out for analyses of uncertainty problems where no analytical expression is available for the forward relationship between data and model parameters. In such cases any direct mathematical treatment is impossible, however the forward relation materializes itself as an algorithm allowing data to be calculated for any given model. This study addresses a new method which is utilized as a basis for the uncertainty estimates of structural responses. It applies double uniform random numbers (i.e. DURN technique) to conventional Monte Carlo algorithm. In DURN method, the scenarios of uncertainties are sequentially selected and executed in its simulation. Numerical examples demonstrate the beneficial effect that the technique can increase uncertainty degree of structural properties with maintaining structural stability and safety up to the limit point of a breakdown of structural systems.

• Nuclear Engineering and Technology
• /
• v.41 no.9
• /
• pp.1191-1204
• /
• 2009

In this study, a CFD methodology is employed to address the problem of the prediction of the flow in a T-junction. An Unsteady Reynolds Averaged Navier-Stokes (URANS) approach has been selected for its low computational cost. Moreover, Unsteady Reynolds Navier-Stokes methodologies do not need complex boundary formulations for the inlet and the outlet such as those required when using Large Eddy Simulation (LES) or Direct Numerical Simulation (DNS). The results are compared with experimental data and an LES calculation. In the past, URANS has been tried on T-junctions with mixed results. The biggest limit observed was the underestimation of the oscillatory behavior of the temperature. In the present work, we propose a comprehensive approach able to correctly reproduce the root mean square (RMS) of the temperature directly downstream of the T-junction for cases where buoyancy is not present.