• Journal of the Earthquake Engineering Society of Korea
• /
• v.8 no.2
• /
• pp.19-25
• /
• 2004

The linear analysis for the balance of linear momentum of a structure is relatively easy to perform, but the error becomes large when the structure experiences large deformation. Therefore, the material and geometric nonlinearity need to be considered for the precise calculations in that case. The plastic flow of a ductile steel-like metal mainly transforms its dissipated mechanical energy into heat, which transfers under the first and second law of thermodynamics. This heat increases the temperature of the material and the strength of the material decreases accordingly, which affects mechanical behavior of the given structure. This paper presents a finite-strain thermo-elasto-plastic steel model. This model can handle large deformation and thermal load simultaneously, which is common during earthquake periods. Two 3-dimensional finite element analyses verify this formulation.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.41 no.5
• /
• pp.361-366
• /
• 2017

As the recent development of computing architecture and application software technology, real world simulation, which is the ultimate destination of computer simulation, is emerging as a practical issue in several research sectors. In this paper, metal plate motion in a square shock tube for small time interval was calculated using a supercomputing-based fluid-structure-combustion multi-physics simulation tool called Illinois Rocstar, developed in a US national R amp; D program at the University of Illinois. Afterwards, the simulation results were compared with those from experiments. The coupled solvers for unsteady compressible fluid dynamics and for structural analysis were based on the finite volume structured grid system and the large deformation linear elastic model, respectively. In addition, a strong correlation between calculation and experiment was shown, probably because of the predictor-corrector time-integration scheme framework. In the future, additional validation studies and code improvements for higher accuracy will be conducted to obtain a reliable open-source software research tool.

• Geophysics and Geophysical Exploration
• /
• v.4 no.2
• /
• pp.25-33
• /
• 2001

Loop-loop electromagnetic (EM) survey in frequency domain has been carried out in order to provide basic solution to geotechnical applications. Source and receiver configuration may be horizontal co-planar (HCP) and/or vertical co-planar (VCP). Three quadrature components of mutual impedance ratio for each configuration are used to construct the subsurface image. For the purpose of obtaining the model response and validating the reasonable performance of the inversion, we obtained each responses of two-layered and three-layered earth models and two-dimensional (2-D) isolated anomalous body. The response of 2-D isolated anomalous body has been calculated using extended Born approximation for the solution of 2.5-D integral equation describing EM scattering problem. As a result of the least-squares inversion with variable Lagrangian multiplier, we could construct more resolvable image from HCP data than VCP data. Furthermore, joint inversion of HCP and VCP data made better stability and resolution of the inversion. Resistivity values, however, did not exactly match the true ones. Loop-loop EM field data was obtained with EM34-3XL system manufactured by Geonics Ltd. (Canada). Electrical resistivity survey was conducted on the same line for the comparison in advance. Since the constructed image from loop-loop EM data by 2-D inversion algorithm showed almost similar resistivity distribution to that from electrical resistivity one, we expect the developed 2.5-D loop-loop EM inversion program can be applied for the reconnaissance site survey.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.29 no.6
• /
• pp.583-590
• /
• 2016

This paper presents a dynamic crack propagation algorithm with Rayleigh damping effect based on the MLS(Moving Least Squares) Difference Method. Dynamic equilibrium equation and constitutive equation are derived by considering Rayliegh damping and governing equations are discretized by the MLS derivative approximation; the proportional damping, which has not been properly treated in the conventional strong formulations, was implemented in both the equilibrium equation and constitutive equation. Dynamic equilibrium equation including time relevant terms is integrated by the Central Difference Method and the discrete equations are simplified by lagging the velocity one step behind. A geometrical feature of crack is modeled by imposing the traction-free condition onto the nodes placed at crack surfaces and the effect of movement and addition of the nodes at every time step due to crack growth is appropriately reflected on the construction of total system. The robustness of the proposed numerical algorithm was proved by simulating single and multiple crack growth problems and the effect of proportional damping on the dynamic crack propagation analysis was effectively demonstrated.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.32 no.4
• /
• pp.215-222
• /
• 2019

The heat conduction problem with discontinuous material coefficients generally consists of the conservative equation, boundary condition, and interface condition, which should be additionally satisfied in the solution procedure. This feature often makes the development of new numerical schemes difficult as it induces a layered singularity in the solution fields; thus, a special approximation is required to capture the singular behavior. In addition to the approximation, the construction of a total system of equations is challenging. In this study, a wedge function is devised for enriching the approximation, and the interface condition itself is embedded in the moving least squares(MLS) derivative approximation to consistently satisfy the interface condition. The heat conduction problem is then discretized in a strong form using the developed derivative approximation, which is named as the interface immersed MLS difference method. This method is able to efficiently provide a numerical solution for such interface problems avoiding both numerical quadrature as well as extra difference equations related to the interface condition enforcement. Numerical experiments proved that the developed numerical method was highly accurate and computationally efficient at solving the heat conduction problem with interfacial jump as well as the problem with a geometrically induced interfacial singularity.

• Atmosphere
• /
• v.29 no.1
• /
• pp.87-103
• /
• 2019

This study introduces the overall characteristics of LOVECLIM version 1.3, the earth system model of intermediate complexity (EMIC), including the installation and operation processes by conducting two kinds of past climate simulation. First climate simulation is the equilibrium experiment during the mid-Holocene (6,000 BP), when orbital parameters were different compared to those at present. The overall accuracy of simulated global atmospheric fields by LOVECLIM is relatively lower than that in Coupled Model Intercomparison Project phase 5 (CMIP5) and Paleoclimate modelling Intercomparison Project phase 3 (PMIP3) simulations. However, surface temperature over the globe, the 800 hPa meridional wind over the mid-latitude coastal region, and the 200 hPa zonal wind from LOVECLIM show similar spatial distribution to those multi-model mean of CMIP5/PMIP3 climate models. Second one is the transient climate experiment from mid-Holocene to present. LOVECLIM well captures the major differences in surface temperature between preindustrial and mid-Holocene simulations by CMIP5/PMIP3 multi-model mean, even though it was performed with short integration time (i.e., about four days in a single CPU environment). In this way, although the earth system model of intermediate complexity has a limit due to its relatively low accuracy, it can be a very useful tool in the specific research area such as paleoclimate.

• The Journal of the Acoustical Society of Korea
• /
• v.27 no.3
• /
• pp.140-148
• /
• 2008

In this study, a numerical method for a time-domain acoustic wave backscattering analysis is established based on a physical optics and a Fourier transform. The frequency responses of underwater targets are calculated based on physical optics derived from the Kirchhoff-Helmholtz integral equation by applying Kirchhoff approximation and the time-domain signals are simulated taking inverse fast Fourier transform to the obtained frequency responses. Particularly, the adaptive triangular beam method is introduced to calculate the areas impinged directly by acoustic incident wave and the virtual surface concept is adopted to consider the multiple reflection effect. The numerical analysis result for an acoustic plane wave field incident normally upon a square flat plate is coincident with the result by the analytic time-domain physical optics derived theoretically from a conventional physical optics. The numerical simulation result for a hemi-spherical end-capped cylinder model is compared with the measurement result, so that it is recognized that the presented method is valid when the specular reflection effect is predominant, but, for small targets, gives errors due to higher order scattering components. The numerical analysis of an idealized submarine shows that the established method is effectively applicable to large and complex-shaped underwater targets.

• Journal of the Korea Society of Computer and Information
• /
• v.16 no.10
• /
• pp.101-109
• /
• 2011

The Riemann's zeta function $\zeta(s)$ has been known as answer for a number of primes $\pi$(x) less than given number x. In prime number theorem, there are another approximation function $\frac{x}{lnx}$,Li(x), and R(x). The error about $\pi$(x) is R(x) < Li(x) < $\frac{x}{lnx}$. The logarithmic integral function is Li(x) = $\int_{2}^{x}\frac{1}{lnt}dt$ ~ $\frac{x}{lnx}\sum\limits_{k=0}^{\infty}\frac{k!}{(lnx)^k}=\frac{x}{lnx}(1+\frac{1!}{(lnx)^1}+\frac{2!}{(lnx)^2}+\cdots)$. This paper shows that the $\pi$(x) can be represent with finite Li(x), and presents generalized prime counting function $\sqrt{{\alpha}x}{\pm}{\beta}$. Firstly, the $\pi$(x) can be represent to $Li_3(x)=\frac{x}{lnx}(\sum\limits_{t=0}^{{\alpha}}\frac{k!}{(lnx)^k}{\pm}{\beta})$ and $Li_4(x)=\lfloor\frac{x}{lnx}(1+{\alpha}\frac{k!}{(lnx)^k}{\pm}{\beta})}k\geq2$ such that $0{\leq}t{\leq}2k$. Then, $Li_3$(x) is adjusted by $\pi(x){\simeq}Li_3(x)$ with ${\alpha}$ and error compensation value ${\beta}$. As a results, this paper get the $Li_3(x)=Li_4(x)=\pi(x)$ for $x=10^k$. Then, this paper suggests a generalized function $\pi(x)=\sqrt{{\alpha}x}{\pm}{\beta}$. The $\pi(x)=\sqrt{{\alpha}x}{\pm}{\beta}$ function superior than Riemann's zeta function in representation of prime counting.

• Korean Journal of Fisheries and Aquatic Sciences
• /
• v.34 no.6
• /
• pp.643-651
• /
• 2001

Paddlewheel-driven circulation in a culture pond has been simulated based on the depth integrated 2 dimensional hydrodynamic model. Acceleration by paddlewheel is expressed as shaft force divided by water mass discharged by paddlewheel blades. The model has been calibrated and applied to culture ponds as following steps:- i) The model predicted velocities at every 10 m along longitudinal direction from the paddlewheel. The model was calibrated comparing the results with the measured values at mass correction factor $\alpha$ and dimensionless eddy viscosity constant $\gamma$, respectively, in a range $15\~20$ and 6. ii) Wind shear stress was simulated under conditions of direction $0^{\circ}C,\;90^{\circ}C\;and\;180^{\circ}C$ and speed 0.0, 2.5, 5.0 and 7.5 m/s. Change rate of current speed was <$1\%$ at wind in parallel or opposite direction to the paddlewheel-driven jet flow, while $4\%$ at orthogonal angle. iii) The model was then applied to 2 culture ponds located at the Western coast of Korea. The measured and predicted currents for the ponds were compared using the regression analysis. Analysis of flow direction and speed showed correlation coefficients 0.8928 and 0.6782 in pond A, 0.8539 and 0.7071 in pond B, respectively. Hence, the model is concluded to accurately predict circulation driven by paddlewheel such that it can be a useful tool to provide pond management strategy relating to paddlewheel operation and water quality.