• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.9 no.5
• /
• pp.678-686
• /
• 1998

This paper describes an anisotropic perfectly matched layer (APML) for mesh truncation of the Finite Difference Time Domain(FDTD) method. The APML method can classified by a split type and an unsplit type, in case of the split type be made up 12 equations or 8 equations largely, and in case of the unsplit type be made of 6 equations. Therefore the latter is more simple as compare with the former. For the APML method presented in this paper is the latter, is directly derived from the time domain equations of Maxwell and extend to the three dimensional problem for the method suggested by Chen. Especially, in the edge and corner parts except the planes, the APML method effectively treated as compound with the Chen's method and Gedney's method newly. The results of the numerical method in this paper show the radiation patterns and the time responses of electromagnetic fields of the wire antennas according to wavelengths and the APML results are compared with Mur's first order absorbing boundary condition and Kraus's analytical results. Eventually, Mur's first order condition have many errors at the edge and corner. On the other hand, in comparison with Kraus's analytical results, it is quite good agreement, and the validity of present method is confirmed.

• Advances in Computational Design
• /
• v.5 no.4
• /
• pp.397-416
• /
• 2020

Functionally graded material (FGM) illustrates a novel class of composites that consists of a graded pattern of material composition. FGM is engineered to have a continuously varying spatial composition profile. Current work focused on buckling analysis of beam made of stepwise linear and quadratic graded material in axial direction subjected to axial span-load with piecewise function and rested on shear layer based on classical beam theory. The various boundary and natural conditions including simply supported (S-S), pinned - clamped (P-C), axial hinge - pinned (AH-P), axial hinge - clamped (AH-C), pinned - shear hinge (P-SHH), pinned - shear force released (P-SHR), axial hinge - shear force released (AH-SHR) and axial hinge - shear hinge (AH-SHH) are considered. To the best of the author's knowledge, buckling behavior of this kind of Euler-Bernoulli beams has not been studied yet. The equilibrium differential equation is derived by minimizing total potential energy via variational calculus and solved analytically. The boundary conditions, natural conditions and deformation continuity at concentrated load insertion point are expressed in matrix form and nontrivial solution is employed to calculate first buckling loads and corresponding mode shapes. By increasing truncation order, the relative error reduction and convergence of solution are observed. Fast convergence and good compatibility with various conditions are advantages of the proposed method. A MATLAB code is provided in appendix to employ the numerical procedure based on proposed method.

• Coupled systems mechanics
• /
• v.8 no.5
• /
• pp.393-414
• /
• 2019

This paper presents seismic analysis of concrete gravity dams considering soil-structure-fluid interaction. Displacement based plane strain finite element formulation is considered for the dam and foundation domain whereas pressure based finite element formulation is considered for the reservoir domain. A direct coupling method has been adopted to obtain the interaction effects among the dam, foundation and reservoir domain to obtain the dynamic responses of the dam. An efficient absorbing boundary condition has been implemented at the truncation surfaces of the foundation and reservoir domains. A parametric study has been carried out considering each domain separately and collectively based on natural frequencies, crest displacement and stress at the neck level of the dam body. The combined frequency of the entire coupled system is very less than that of the each individual sub-system. The crest displacement and neck level stresses of the dam shows prominent enhancement when coupling effect is taken into consideration. These outcomes suggest that a complete coupled analysis is necessary to obtain the actual responses of the concrete gravity dam. The developed methodology can easily be implemented in finite element code for analyzing the coupled problem to obtain the desired responses of the individual subdomains.

• Kyungpook Mathematical Journal
• /
• v.62 no.4
• /
• pp.699-713
• /
• 2022

This paper proposes a hybrid successive over-relaxation iterative method for the numerical solution of a nonlinear diffusion in a one-dimensional porous medium. The considered mathematical model is discretized using a computational complexity reduction scheme called half-sweep finite differences. The local truncation error and the analysis of the stability of the scheme are discussed. The proposed iterative method, which uses explicit group technique and modified successive over-relaxation, is formulated systematically. This method improves the efficiency of obtaining the solution in terms of total iterations and program elapsed time. The accuracy of the proposed method, which is measured using the magnitude of absolute errors, is promising. Numerical convergence tests of the proposed method are also provided. Some numerical experiments are delivered using initial-boundary value problems to show the superiority of the proposed method against some existing numerical methods.

• Economic and Environmental Geology
• /
• v.24 no.4
• /
• pp.435-446
• /
• 1991

Ductile shear zones developed in Jurassic granites in the Yonggwang area show NE trend at the eastern part and nearly EW trend at the western part, respectively. Judged from shear sense indicators, they have resulted from dextral strike-slip movement. The intersection of both trends is thought to be due to the truncation and offset of NE shear zone Chonju Shear zone by the brittle Yonggwang fault which runs in near EW direction with sinistral movement sense. The simple shear deformation was predominate through the deformation in this ductile shear zone. Based on this deformation mechanism, the shear strain (${\gamma}$) estimated in domain 1 increases from 0.14 at the shear zone margin to 9.41 toward the center of shear zone. Total displacement obtained only from this measured section(JK 59 to JK14) appecars to be 1434.5 meters. The sequential development of microstructures can be divided into three stages; weakly-foliated, well-foliated and banded-foliated stages. In the weakly-foliated stage dislocation glide mechanism might be predominant. In the well-foliated stage grain boundary migration and progressive misorientation of subgrains was remarkable during dynamic recovery and recrystallization. In the banded-foliated stage grain boundary sliding and microfracturing mechanisms accompanied with crushing and cracking were marked. According to strain analysis from quartzites of the metasedimentary rocks, strain intensity (${\gamma}$) of the samples within the ductile shear zone ranges from 2.7 to 5.7, while that of the samples out of the ductile shear zone appears to be about 1.7.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.50 no.5
• /
• pp.297-308
• /
• 2022

Most air-breathing aircraft operated in the hypersonic region are equipped with a scramjet engine. In a scramjet engine, a shock wave generated at an inlet acts as a compressor for a general gas turbine engine instead, so total pressure loss caused by the shock wave is considered very important. In this study, to minimize total pressure loss, a method of designing an external compression inlet using isentropic compression surface was proposed, and an external compression inlet with 3-deflection angles and Busemann inlet were designed under the same conditions. After that, through computational analysis, the performance characteristics at off-design conditions were compared. Each inlet shape was truncated according to the length of the 3-ramp external compression inlet, and the boundary layer correction was performed. The isentropic external compression inlet showed superior performance at the design point, but under the off-design conditions, its performance was degraded compared to the 3-ramp external compression inlet.

• Journal of applied mathematics & informatics
• /
• v.15 no.1_2
• /
• pp.29-51
• /
• 2004

In this paper, we propose a new mixed finite element method, called the characteristics-mixed method, for approximating the solution to Burgers' equation. This method is based upon a space-time variational form of Burgers' equation. The hyperbolic part of the equation is approximated along the characteristics in time and the diffusion part is approximated by a mixed finite element method of lowest order. The scheme is locally conservative since fluid is transported along the approximate characteristics on the discrete level and the test function can be piecewise constant. Our analysis show the new method approximate the scalar unknown and the vector flux optimally and simultaneously. We also show this scheme has much smaller time-truncation errors than those of standard methods. Numerical example is presented to show that the new scheme is easily implemented, shocks and boundary layers are handled with almost no oscillations. One of the contributions of the paper is to show how the optimal error estimates in $L^2(\Omega)$ are obtained which are much more difficult than in the standard finite element methods. These results seem to be new in the literature of finite element methods.

• Structural Engineering and Mechanics
• /
• v.44 no.1
• /
• pp.85-107
• /
• 2012

The building frame and its foundation along with the soil on which it rests, together constitute a complete structural system. In the conventional analysis, a structure is analysed as an independent frame assuming unyielding supports and the interactive response of soil-foundation is disregarded. This kind of analysis does not provide realistic behaviour and sometimes may cause failure of the structure. Also, the conventional analysis considers infill wall as non-structural elements and ignores its interaction with the bounding frame. In fact, the infill wall provides lateral stiffness and thus plays vital role in resisting the seismic forces. Thus, it is essential to consider its effect especially in case of high rise buildings. In the present research work the building frame, infill wall, isolated column footings (open foundation) and soil mass are considered to act as a single integral compatible structural unit to predict the nonlinear interaction behaviour of the composite system under seismic forces. The coupled isoparametric finite-infinite elements have been used for modelling of the interaction system. The material of the frame, infill and column footings has been assumed to follow perfectly linear elastic relationship whereas the well known hyperbolic soil model is used to account for the nonlinearity of the soil mass.

• Coupled systems mechanics
• /
• v.3 no.3
• /
• pp.267-289
• /
• 2014

A proper physical modeling of infilled building frame-foundation beam-soil mass interaction system is needed to predict more realistic and accurate structural behavior under static vertical loading. This is achieved via finite element method considering the superstructure, foundation and soil mass as a single integral compatible structural unit. The physical modelling is achieved via use of finite element method, which requires the use of variety of isoparametric elements with different degrees of freedom. The unbounded domain of the soil mass has been discretized with coupled finite-infinite elements to achieve computational economy. The nonlinearity of soil mass plays an important role in the redistribution of forces in the superstructure. The nonlinear behaviour of the soil mass is modeled using hyperbolic model. The incremental-iterative nonlinear solution algorithm has been adopted for carrying out the nonlinear elastic interaction analysis of a two-bay two-storey infilled building frame. The frame and the infill have been considered to behave in linear elastic manner, whereas the subsoil in nonlinear elastic manner. In this paper, the computational methodology adopted for nonlinear soil-structure interaction analysis of infilled frame-foundation-soil system has been presented.

• Journal of the Society of Naval Architects of Korea
• /
• v.48 no.3
• /
• pp.233-237
• /
• 2011

In the present experiments, vortical structures behind the hydrofoil trailing edge are visualized and analyzed as an elementary study for propeller singing phenomena. Two sorts of hydrofoil are selected for the measurement of shedding vortices. One was KH45 hydrofoil section and the other is KH45 with the truncated trailing edge that is positioned at X/C = 0.9523(C=chord length). Assuming the Strouhal number of 0.23, the shedding frequencies of vortices are extracted by analyzing the boundary layer thickness and the flow speed. The frequency distribution of shedding vortices is obtained with the variation of angle-of-attack while the flow speed is fixed to 8m/s. The truncation of the trailing edge makes the frequency of shedding vortices about 120Hz lower than that of original trailing edge and makes the vorticity value higher than the original trailing edge.