• Journal of Advanced Marine Engineering and Technology
• /
• v.29 no.4
• /
• pp.391-398
• /
• 2005

In this investigation, turbulent premixed combustion and flame front propagation in a gas turbine combustion chamber is studied. Direct numerical simulation of turbulent reacting flows demands extremely high computational resources, especially in more complicated geometry. The alternative choice may be left for Large Eddy Simulation (LES) by which only large scales are solved directly. In combustion problems, capturing the large scales' behavior without solving the details of small scales is a difficult task. Using a transport equation for description of the flame front propagation and therefore avoiding the calculation of inner flame structure is the basic idea of this study. For this purpose. the so-called G-equation has been used by which any iso-level of the G variable provides the flame location. A comparison with the experiment indicates that the present method can predict a turbulent velocity field and also capture a instantaneous 3-dimensional flame structure.

• Computers and Concrete
• /
• v.2 no.5
• /
• pp.345-366
• /
• 2005

This work presents an assessment of the computational performance of a vector-parallel implementation of probabilistic model for concrete cracking in 3D. This paper shows the continuing efforts towards code optimization as reported in earlier works Paz, et al. (2002a,b and 2003). The probabilistic crack approach is based on the direct Monte Carlo method. Cracking is accounted by means of 3D interface elements. This approach considers that all nonlinearities are restricted to interface elements modeling cracks. The heterogeneity governs the overall cracking behavior and related size effects on concrete fracture. Computational kernels in the implementation are the inexact Newton iterative driver to solve the non-linear problem and a preconditioned conjugate gradient (PCG) driver to solve linearized equations, using an element by element (EBE) strategy to compute matrix-vector products. In particular the paper analyzes code behavior using OpenMP directives in parallel vector processors (PVP), such as the CRAY SV1 and CRAY T94. The impact of the memory architecture on code performance, and also some strategies devised to circumvent this issue are addressed by numerical experiment.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.4 no.2
• /
• pp.123-131
• /
• 2012

When the object is traveling in the water at tremendously high speeds, the cavity forms and grows up at a fore part of the object called cavitator, and the object is eventually enveloped by vaporized water, supercavitation. As a result, the only part of the object in direct contact with the water is the cavitator, so skin-friction drag is significantly reduced. This is why recently supercavitating objects have been interested in many applicable fields. In this study we are focused out attention on supercavitating flows around various shapes of two and three dimensional cavitators. First, general features of supercavitation are examined by analyzing results obtained by the previously developed numerical method. Second, experimental observations are carried out at a cavitation tunnel at the Chungnam National University (CNU CT), and supercavity dimensions are scrutinized.

• Elastomers and Composites
• /
• v.39 no.3
• /
• pp.228-233
• /
• 2004

A bushing is a device used in automotive suspension systems to reduce the load transmitted from the wheel to the frame of the vehicle. A bushing is a hollow cylinder, which is bonded to a solid steel shaft at its inner surface and a steel sleeve at its outer surface. The relation between the force applied to the shaft and the relative deformation of a bushing is nonlinear and exhibits features of viscoelasticity. Since a force-displacement relation for bushings is important for multibody dynamics numerical simulations, the relation is expressed in terms of a force relaxation function and a method of determination by experiments on bushings has been developed. For the nonlinear viscoelastic axial response, Pipkin-Rogers model, the direct relation of force and displacement, has been derived from experiment. It is shown that the predictions by the proposed force-displacement relation are in very good agreement with the experimental results.

• Nuclear Engineering and Technology
• /
• v.53 no.3
• /
• pp.732-740
• /
• 2021

Scattering resonance of medium mass nuclides leads complex spectrum in the fast reactor, which requires thousands of energy groups in the spectrum calculation. When the broad-group cross sections are collapsed, reaction rate cannot be completely conserved. To eliminate the error from energy collapsing, the Super-homogenization method in energy collapsing (ESPH) was employed in the fast reactor code SARAX. An ESPH factor was derived based on the ESPH-corrected SN transport equation. By applying the factor in problems with reflective boundary condition, both the effective multiplication factor and reaction rate were conserved. The fixed-source iteration was used to ensure the stability of ESPH iteration. However, in the energy collapsing process of SARAX, the vacuum boundary condition was adopted, which was necessary for fast reactors with strong heterogeneity. To further reduce the error caused by leakage, an additional conservation factor was proposed to correct the neutron current in energy collapsing. To evaluate the performance of ESPH with conservation factor, numerical benchmarks of fast reactors were calculated. The results of broad-group calculation agreed well with the direct full-core Monte-Carlo calculation, including the effective multiplication factor, radial power distribution, total control rod worth and sodium void worth.

• Geomechanics and Engineering
• /
• v.31 no.3
• /
• pp.305-318
• /
• 2022

Permeable reactive barriers used for groundwater treatment require proper estimation of the reactive material behavior regarding the emplacement method. This study evaluates the dry emplacement of zeolite (clinoptilolite) to be used as a reactive material in the barrier by carrying out several geotechnical laboratory tests. Dry zeolite samples, exhibited higher wetting-induced compression strains at the higher vertical stresses, up to 12% at 400 kN/m². The swelling potential was observed to be limited with a 3.5 swell index and less than 1% free swelling strain. Direct shear tests revealed that inundation reduces the shear strength of a dry zeolite column by a maximum of 10%. Falling head permeability tests indicate decreasing permeability values with increasing the vertical effective stress. Regarding self-loading and inundation, the porosity along the zeolite column was calculated using a proposed 1D numerical model to predict the permeability with depth considering the laboratory tests. The calculated discharge efficiency was significantly decreased with depth and less than 2% relative to the top for barrier depths deeper than 20 m. Finally, the importance of directional dependence in the permeability of the zeolite medium for calibrating 2D finite element flow analysis was highlighted by bench-scale tests performed under 2D flow conditions.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.1 no.1
• /
• pp.1-10
• /
• 1997

Current seismic design philosophy for reinforced concrete (RC) structures on energy dissipation through large inelastic defomations. A nonlinear dynamic analysis which is used to represent this behavior is time consuming and expensive, particularly if the computations have to be repeated many times. Therefore, the selection of an efficient yet accurate alogorithm becomes important. The main objective of the present study is to propose a new technique herein called the prediction approach with siffness measure (PASM) method in the convetional direct integration methods, the triangular decomposition of matrix is required for solving equations of motion in every time step or every iteration. The PASM method uses a limited number of predetermined decomposed effective matrices obtained from stiffness states of the structure when it is deformed into the nonlinear range by statically applied cyclic loading. The method to be developed herein will reduce the overall numerical effort when compared to approaches which recompute the stiffness in each time step or iteration.

• Steel and Composite Structures
• /
• v.18 no.4
• /
• pp.985-1000
• /
• 2015

The method of cross-section analysis for different sections in a structural frame has been widely investigated since the 1960s for determination of sectional capacities of beam-columns. Many hand-calculated equations and design graphs were proposed for the specific shape and type of sections in pre-computer age decades ago. In design of many practical sections, these equations may be uneconomical and inapplicable for sections with irregular shapes, leading to the high construction cost or inadequate safety. This paper not only proposes a versatile numerical procedure for sectional analysis of beam-columns, but also suggests a method to account for residual stress and geometric imperfections separately and the approach is applied to design of high strength steels requiring axial force-moment interaction for advanced analysis or direct analysis. A cross-section analysis technique that provides interaction curves of arbitrary welded sections with consideration of the effects of residual stress by meshing the entire section into small triangular fibers is formulated. In this study, two doubly symmetric sections (box-section and H-section) fabricated by high-strength steel is utilized to validate the accuracy and efficiency of the proposed method against a hand-calculation procedure. The effects of residual stress are mostly not considered explicitly in previous works and they are considered in an explicit manner in this paper which further discusses the basis of the yield surface theory for design of structures made of high strength steels.

• The Journal of The Korea Institute of Intelligent Transport Systems
• /
• v.20 no.3
• /
• pp.74-85
• /
• 2021

Global Navigation Satellite System (GNSS) receivers are intrinsically vulnerable to radio frequency jamming signals due to the fundamental property of radio navigation systems. A GNSS jamming monitoring system that is capable of jamming detection, classification and localization is essential for infrastructure for autonomous driving systems. For these 3 functionalities, a GNSS jamming monitoring network consisting of a multiple of low-cost GNSS receivers distributed in a certain area is needed, and the precise time synchronizaion between multiple independent GNSS receivers in the network is an essential element. This paper presents a precise time synchronization method based on the direct use of Time Difference of Arrival (TDOA) technique in signal domain. A block interpolation method is additionally incorporated into the method in order to maintain the precision of time synchronization even with the relatively low sampling rate of the received signals for computational efficiency. The feasibility of the proposed approach is verified in the numerical simualtions.

• Journal of the Korean GEO-environmental Society
• /
• v.13 no.2
• /
• pp.27-39
• /
• 2012

This research is to propose the reinforcing method and design code for the lateral behaviors of the abutment displacement induced from the rainfall infiltration on high landfill slope. First, to make the proper numerical analysis, in-situ soil (weathered granite soil) was taken, and the variance of strength parameters according to water content variance was examined by undrained direct shear test, furthermore, other soil parameters were calculated from the standard penetration test such as elastic modulus and Poisson's ratio etc,. Those parameters were used to calculate the lateral behavior of abutment by finite element method and the member force of pile in high landfill slope according to rainfall infiltration . From the results, the shoe displacement on abutment was calculated as 8.98cm, which is 3 times bigger than the allowable displacement, 3cm. To reinforce it, several reinforcing methods were selected and analyzed such as reinforced retaining wall, soil surcharge, pile reinforcing (5m enlargement, 3-line arrangement, 5m enlargement and 3-line arrangement). In case of 5m enlarged and 3-line arrangement piles, the lateral behavior of shoe showed lower value(2.26 cm) than allowable displacement.