• Journal of the Society of Naval Architects of Korea
• /
• v.30 no.1
• /
• pp.104-114
• /
• 1993

In this paper, effects of free surface on viscous flow is investigated. Continuity equation coupled with Navier-Stokes equations are solved numerically by using an artificial compressibility method[1, 2]. The body-fitted generalized curvilinear coordinate system is employed to deal with arbitrary body shape. The IAF scheme with finite difference method is used to solve the equations, and a diagonal algorithm is applied to time-varying Jacobian matrices for the computational economics. Free surface shape is obtained by applying zero pressure condition to still water surface at each time step. A numerical test is made for larminar flow around a circular cylinder vertically piercing the free surface. Computed flow patterns are largely affected by the existance of free surface in low Reynolds number flows treated in this paper. Free surface causes viscous pressure drag to vary much in depth direction in accordance with the variations of flow pattern.

• Structural Engineering and Mechanics
• /
• v.13 no.2
• /
• pp.135-154
• /
• 2002

This paper presents the convected material frame approach to study the nonlinear behavior of inelastic frame structures. The convected material frame approach is a modification of the co-rotational approximation by incorporating an adaptive convected material frame in the basic definition of the displacement vector and strain tensor. In the formulation, each discrete element is associated with a local coordinate system that rotates and translates with the element. For each load increment, the corresponding strain-displacement and nodal force-stress relationships are defined in the updated local coordinates, and based on the updated element geometry. The rigid body motion and deformation displacements are decoupled for each increment. This modified approach incorporates the geometrical nonlinearities through the continuous updating of the material frame geometry. A generalized nonlinear function is used to derive the inelastic constitutive relation and the kinematic hardening is considered. The equation of motion is integrated by an explicit procedure and it involves only vector assemblage and vector storage in the analysis by assuming a lumped mass matrix of diagonal form. Several numerical examples are demonstrated in close agreement with the solutions obtained by the ANSYS code. Numerical studies show that the proposed approach is capable of investigating large deflection of inelastic planar structures and providing an excellent numerical performance.

• Computers and Concrete
• /
• v.31 no.3
• /
• pp.223-239
• /
• 2023

Geometric nonlinear performance simulation and analysis of complex modern buildings and industrial products require high-performance shell elements. Balancing multiple aspects of performance in the one geometric nonlinear analysis element remains challenging. We present a new shell element, flat shell DKMGQ-CR (Co-rotational Discrete Kirchhoff-Mindlin Generalized Conforming Quadrilateral), for linear and geometric nonlinear analysis of both thick and thin shells. The DKMGQ-CR shell element was developed by combining the advantages of high-performance membrane and plate elements in a unified coordinate system and introducing the co-rotational formulation to adapt to large deformation analysis. The effectiveness of linear and geometric nonlinear analysis by DKMGQ-CR is verified through the tests of several classical numerical benchmarks. The computational results show that the proposed new element adapts to mesh distortion and effectively alleviates shear and membrane locking problems in linear and geometric nonlinear analysis. Furthermore, the DKMGQ-CR demonstrates high performance in analyzing thick and thin shells. The proposed element DKMGQ-CR is expected to provide an accurate, efficient, and convenient tool for the geometric nonlinear analysis of shells.

• Proceedings of the Korean Society of Medical Physics Conference
• /
• 2002.09a
• /
• pp.53-60
• /
• 2002

Motion of lung tumors from respiration has been reported in the literature to be as large as of 1-2 cm. This motion requires an additional margin between the Clinical Target Volume (CTV) and the Planning Target Volume (PTV). While such a margin is necessary, it may not be sufficient to ensure proper delivery of Intensity Modulated Radiotherapy (IMRT) to the CTV during the simultaneous movement of the DMLC. Gated treatment has been proposed to improve normal tissues sparing as well as to ensure accurate dose coverage of the tumor volume. The following questions have not been addressed in the literature: a) what is the dose error to a target volume without gated IMRT treatment\ulcorner b) what is an acceptable gating window for such treatment. In this study, we address these questions by proposing a novel technique for calculating the 3D dose error that would result if a lung IMRT plan were delivered without gating. The method is also generalized for gated treatment with an arbitrary triggering window. IMRT plans for three patients with lung tumor were studied. The treatment plans were generated with HELIOS for delivery with 6 MV on a CL2100 Varian linear accelerator with a 26 pair MLC. A CTV to PTV margin of 1 cm was used. An IMRT planning system searches for an optimized fluence map ${\Phi}$ (x,y) for each port, which is then converted into a dynamic MLC file (DMLC). The DMLC file contains information about MLC subfield shapes and the fractional Monitor Units (MUs) to be delivered for each subfield. With a lung tumor, a CTV that executes a quasi periodic motion z(t) does not receive ${\Phi}$ (x,y), but rather an Effective Incident Fluence EIF(x,y). We numerically evaluate the EIF(x,y) from a given DMLC file by a coordinate transformation to the Target's Eye View (TEV). In the TEV coordinate system, the CTV itself is stationary, and the MLC is seen to execute a motion -z(t) that is superimposed on the DMLC motion. The resulting EIF(x,y)is inputted back into the dose calculation engine to estimate the 3D dose to a moving CTV. In this study, we model respiratory motion as a sinusoidal function with an amplitude of 10 mm in the superior-inferior direction, a period of 5 seconds, and an initial phase of zero.

• The Mathematical Education
• /
• v.26 no.1
• /
• pp.41-45
• /
• 1987

In this paper, we give several fixed point theorems in a complete metric space for two multi-valued mappings commuting with two single-valued mappings. In fact, our main theorems show the existence of solutions of functional equations f($\chi$)=g($\chi$)$\in$S$\chi$∩T$\chi$ and $\chi$=f($\chi$)=g($\chi$)$\in$S$\chi$∩T$\chi$ under certain conditions. We also answer an open question proposed by Rhoades-Singh-Kulsherestha. Throughout this paper, let (X, d) be a complete metric space. We shall follow the following notations : CL(X) = {A; A is a nonempty closed subset of X}, CB(X)={A; A is a nonempty closed and founded subset of X}, C(X)={A; A is a nonempty compact subset of X}, For each A, B$\in$CL(X) and $\varepsilon$>0, N($\varepsilon$, A) = {$\chi$$\in$X; d($\chi$, ${\alpha}$) < $\varepsilon$ for some ${\alpha}$$\in$A}, E$\sub$A, B/={$\varepsilon$ > 0; A⊂N($\varepsilon$ B) and B⊂N($\varepsilon$, A)}, and (equation omitted). Then H is called the generalized Hausdorff distance function fot CL(X) induced by a metric d and H defined CB(X) is said to be the Hausdorff metric induced by d. D($\chi$, A) will denote the ordinary distance between $\chi$$\in$X and a nonempty subset A of X. Let R$\^$+/ and II$\^$+/ denote the sets of nonnegative real numbers and positive integers, respectively, and G the family of functions ${\Phi}$ from (R$\^$+/)$\^$s/ into R$\^$+/ satisfying the following conditions: (1) ${\Phi}$ is nondecreasing and upper semicontinuous in each coordinate variable, and (2) for each t>0, $\psi$(t)=max{$\psi$(t, 0, 0, t, t), ${\Phi}$(t, t, t, 2t, 0), ${\Phi}$(0, t, 0, 0, t)} $\psi$: R$\^$+/ \longrightarrow R$\^$+/ is a nondecreasing upper semicontinuous function from the right. Before sating and proving our main theorems, we give the following lemmas:

• Proceedings of the Korea Water Resources Association Conference
• /
• 2019.05a
• /
• pp.267-267
• /
• 2019

Driftwood is one of serious problems in a river environment. In several countries, such as Indonesia, Japan, and Italy, the driftwood frequently appears in a river basin, and it can alter the channel bed, flow configuration by wood deposition and jam formation. Therefore, the studies related to driftwood have been actively conducted by many researchers to understand the mechanism of driftwood dynamics. In particular, wood motion by collision is one of the difficult issues in the numerical simulation because the calculation for wood collision requires significantly expensive calculation time due to small time step. Thus, this study conducted the numerical simulation in consideration of the wood motion by water flow and wood collision to understand the wood dynamics in terms of computation. We used the 2D (two-dimensional) depth-averaged velocity model, Nays2DH, which is a Eulerian model to calculate the water flow on the generalized coordinate. A Lagrangian type driftwood model, which expresses the driftwood as connected sphere shape particles, was employed to Nays2DH. In addition, the present study considered root wad effect by using larger diameter for a particle at a head of driftwood. An anisotropic bed friction was considered for the sliding motion dependent on stemwise, streamwise and motion directions. We particularly considered changeable draft at each particle and projection area by an angle between stemwise and flow directions to precisely reproduce the wood motions. The simulation results were compared with experimental results to verify the model. As a result, the simulation results showed good agreement with experimental results. Through this study, it would be expected that this model is a useful tool to predict the driftwood effect in the river flow.

• KIPS Transactions on Computer and Communication Systems
• /
• v.6 no.7
• /
• pp.297-306
• /
• 2017

Citizens want more accurate forecast information using Bus Information System. However, most bus information systems that use an average based short-term prediction algorithm include many errors because they do not consider the effects of the traffic flow, signal period, and halting time. In this paper, we try to improve the precision of forecast information by analyzing the influencing factors of the error, thereby making the convenience of the citizens. We analyzed the influence factors of the error using BIS data. It is shown in the analyzed data that the effects of the time characteristics and geographical conditions are mixed, and that effects on halting time and passes speed is different. Therefore, the halt time is constructed using Generalized Additive Model with explanatory variable such as hour, GPS coordinate and number of routes, and we used Hidden Markov Model to construct a pattern considering the influence of traffic flow on the unit section. As a result of the pattern construction, accurate real-time forecasting and long-term prediction of route travel time were possible. Finally, it is shown that this model is suitable for travel time prediction through statistical test between observed data and predicted data. As a result of this paper, we can provide more precise forecast information to the citizens, and we think that long-term forecasting can play an important role in decision making such as route scheduling.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.37 no.9
• /
• pp.1069-1075
• /
• 2013

An effective method is necessary for the real-time analysis of an unmanned military robot. To achieve highly efficient simulations, a subsystem synthesis method has been developed. The subsystem synthesis method separately generates equations of motion for the base body and for the subsystem. The equations of motion are expressed by DAE, which consist of differential equations and algebraic equations. To increase the accuracy and efficiency of solutions, DAE solvers such as the Direct, CS (Constraint Stabilization), and GCP (Generalized Coordinate Partitioning) method are employed. In this study, the subsystem synthesis method is applied for effective multi-body dynamics analysis of an unmanned military robot, and a comparative study of three different DAE solvers is carried out.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.26 no.2B
• /
• pp.209-216
• /
• 2006

This study is to analyse the characteristics of bar migration, which is important roles to bank erosion and meandering development, for the meandering channels with erodible bed by using a 2-D numerical model in the generalized coordinate systems. The results of the numerical experiments showed that the features of bar migration were affected by the meandering wavelength to the width ratio, and had a relatively good agreement with the criterion for bar migration through a bend suggested by Kinoshita and Miwa (1974). The bar migrated with speed in the channel with long wavelength and broad width, and the criterion for bar migration was increased. The bar celerity was decreased abruptly near the criterion.

• Journal of Korea Water Resources Association
• /
• v.37 no.5
• /
• pp.387-395
• /
• 2004

A computer model was proposed to simulate channel changes and bed material sorting of the meandering channels with different grain size in time and space simultaneously. The bed at the outside of the meandering channel with mixed sediments was scoured deeply and composed of coarser materials, and at the inside was aggradated and composed of finer materials. The sorting process started at the upstream inflection point and was finished at the downstream inflection point. At the natural with complicated boundaries and non-uniform grain sizes, the bed near the outside at the bend and narrow width was scoured deeper with coarse materials than in the channel with uniform grain sizes. The point bars showed lip at the inside near the bend and the bed materials were finer The bed at the outside near the bend and in the narrow width was scoured deeply with the coarser materials.