• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.27 no.11
• /
• pp.1907-1916
• /
• 2003

During decades, there has been much progress in understanding of the inelastic behavior of the materials and numerous inelastic constitutive equations have been developed. The complexity of these constitutive equations generally requires a stable and accurate numerical method. To obtain the increment of state variable, its evolution laws are linearized by several approximation methods, such as general midpoint rule(GMR) or general trapezoidal rule(GTR). In this investigation, semi-implicit integration schemes using GTR and GMR were developed and implemented into ABAQUS by means of UMAT subroutine. The comparison of integration schemes was conducted on the simple tension case, and simple shear case and nonproportional loading case. The fully implicit integration(FI) was the most stable but amplified the truncation error when the nonlinearity of state variable is strong. The semi-implicit integration using GTR gave the most accurate results at tension and shear problem. The numerical solutions with refined time increment were always placed between results of GTR and those of FI. GTR integration with adjusting midpoint parameter can be recommended as the best integration method for viscoplastic equation considering nonlinear kinematic hardening.

• Journal of The Korean Astronomical Society
• /
• v.34 no.4
• /
• pp.281-283
• /
• 2001

As a companion to an adiabatic version developed by Ryu and his coworkers, we have built an isothermal magnetohydrodynamic code for astrophysical flows. It is suited for the dynamical simulations of flows where cooling timescale is much shorter than dynamical timescale, as well as for turbulence and dynamo simulations in which detailed energetics are unimportant. Since a simple isothermal equation of state substitutes the energy conservation equation, the numerical schemes for isothermal flows are simpler (no contact discontinuity) than those for adiabatic flows and the resulting code is faster. Tests for shock tubes and Alfven wave decay have shown that our isothermal code has not only a good shock capturing ability, but also numerical dissipation smaller than its adiabatic analogue. As a real astrophysical application of the code, we have simulated the nonlinear three-dimensional evolution of the Parker instability. A factor of two enhancement in vertical column density has been achieved at most, and the main structures formed are sheet-like and aligned with the mean field direction. We conclude that the Parker instability alone is not a viable formation mechanism of the giant molecular clouds.

• Fire Science and Engineering
• /
• v.10 no.2
• /
• pp.28-39
• /
• 1996

A recently developed electro-thermal simulation methodology is used to analyze the behavior of a PWM(Pulse-Width-Modulated) voltage source inverter which uses IGBT(Insulated Gate Bipolar Transistor) as the switching devices. In the electro-thermal network simulation methdology, the simulator solves for the temperature distribution within the power semiconductor devices(IGBT electro-thermal model), control logic circuitry, the IGBT gate drivers, the thermal network component models for the power silicon chips, package, and heat sinks as well as the current and voltage within the electrical network. The thermal network describes the flow of heat form the chip surface through the package and heat sink and thus determines the evolution of the chip surface temperature used by the power semiconductor device models. The thermal component model for the device silicon chip, packages, and heat sink are developed by discretizing the nonlinear heat diffusion equation and are represented in component from so that the thermal component models for various package and heat sink can be readily connected to on another to form the thermal network.

• Computational Structural Engineering
• /
• v.9 no.1
• /
• pp.65-76
• /
• 1996

The initiation and growth of microcracks or microvoids inside concrete results in the progressive degradation of concrete. This damage processing along processing along with plastic deformation is main cause of nonlinear behavior of concrete. In this study, a continuum damage model of concrete is developed for the analysis of the nonlinear behavior of concrete due to damage and elasto-plastic deformation. Anisotropic damage tensor is used to describe the anisotropy of concrete and hypothesis of equivalent elastic energy is used to define the effective elastic tensor. The damage model including the damage evolution law and constitutive equation is derived with damage variable and damage surface which is defined by damage energy release rate by using the Helmholtz free energy and dissipation potential based on the thermodynamic principles. By adopting a typical plasticity model of concrete, plasticity of concrete is included to this model. Afinite element analysis program implemented with this model was developed and finite element analysis was performed for the analyses of concrete subjected to uniaxial and biaxial loadings. Comparison of the results of analysis with those of experiments and other models shows that the model successfully predicts the nonlinear behavior of concrete.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.28 no.1
• /
• pp.161-167
• /
• 2022

Ocean waves passing through the underwater bar at a shallow depth experience a shoaling effect caused by decreasing water depth, a nonlinear interaction therein owing to steepening wave slope, and a wave dispersion effect as the water depth increases again. Because this problem includes many complicated phenomena, it is used as a good example of validating a theoretical development or a CFD method for ocean wave applications. Validation is performed mainly for regular waves by comparing the wave elevation patterns in the time domain with the experimental results. In this study, the spectral evolution of wave spectrum is investigated in the frequency domain when a CFD method such as OpenFOAM is applied for this problem. In particular, the effects of initial phase conditions as well as the nonlinear interaction among harmonic waves are studied.

• Korea-Australia Rheology Journal
• /
• v.20 no.2
• /
• pp.59-63
• /
• 2008

A theoretical model was developed to describe the flow behavior of a filled polymer in the packing stage of reaction injection molding and predict the residual stress distribution of thin injection-molded parts. The model predictions were compared with experiments performed for phenol-formaldehyde resin filled with wood dust and cured by urotropine. The packing stage of reaction injection molding process presents a typical example of complex non-isothermal flow combined with chemical reaction. It is shown that the time evolution of pressure distribution along the mold cavity that determines the residual stress in the final product can be described by a single 1D partial differential equation (PDE) if the rheological behavior of reacting liquid is simplistically described by the power-law approach with some approximations made for describing cure reaction and non-isothermality. In the formulation, the dimensionless time variable is defined in such a way that it includes all necessary information on the cure reaction history. Employing the routine separation of variables made possible to obtain the analytical solution for the nonlinear PDE under specific initial condition. It is shown that direct numerical solution of the PDE exactly coincides with the analytical solution. With the use of the power-law approximation that describes highly shear thinning behavior, the theoretical calculations significantly deviate from the experimental data. Bearing in mind that in the packing stage the flow is extremely slow, we employed in our theory the Newtonian law for flow of reacting liquid and described well enough the experimental data on evolution of pressure.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.39 no.7
• /
• pp.653-657
• /
• 2015

A numerical investigation to simulate the selective growth of silicon-germanium quantum-dots via local surface diffusivity enhancement is presented. A nonlinear equation for the waviness evolution of film surface is derived to consider the effects of spatially-varying diffusivity, influenced by a surface temperature profile. Results show that the morphology of the initially planar film shapes into an undulated surface upon perturbation, and a steady-state solution describes a fully grown quantum-dot. The present study points toward a fabrication technique that can obtain selectivity for self-assembly.

• Journal of the Korean Mathematical Society
• /
• v.49 no.3
• /
• pp.585-604
• /
• 2012

In this paper we consider the evolution of the rolling stone with a rotationally symmetric nonconvex compact initial surface ${\Sigma}_0$ under the Gauss curvature flow. Let $X:S^n{\times}[0,\;{\infty}){\rightarrow}\mathbb{R}^{n+1}$ be the embeddings of the sphere in $\mathbb{R}^{n+1}$ such that $\Sigma(t)=X(S^n,t)$ is the surface at time t and ${\Sigma}(0)={\Sigma}_0$. As a consequence the parabolic equation describing the motion of the hypersurface becomes degenerate on the interface separating the nonconvex part from the strictly convex side, since one of the curvature will be zero on the interface. By expressing the strictly convex part of the surface near the interface as a graph of a function $z=f(r,t)$ and the non-convex part of the surface near the interface as a graph of a function $z={\varphi}(r)$, we show that if at time $t=0$, $g=\frac{1}{n}f^{n-1}_{r}$ vanishes linearly at the interface, the $g(r,t)$ will become smooth up to the interface for long time before focusing.

• Journal of the Korean Mathematical Society
• /
• v.53 no.1
• /
• pp.161-185
• /
• 2016

We consider a mathematical model which describes the quasistatic frictional contact between a piezoelectric body and an electrically conductive obstacle, the so-called foundation. A nonlinear electro-viscoelastic constitutive law is used to model the piezoelectric material. Contact is described with Signorini's conditions and a version of Coulomb's law of dry friction in which the adhesion of contact surfaces is taken into account. The evolution of the bonding field is described by a first order differential equation. We derive a variational formulation for the model, in the form of a system for the displacements, the electric potential and the adhesion. Under a smallness assumption which involves only the electrical data of the problem, we prove the existence of a unique weak solution of the model. The proof is based on arguments of time-dependent quasi-variational inequalities, differential equations and Banach's fixed point theorem.

• Journal of The Korean Astronomical Society
• /
• v.34 no.4
• /
• pp.285-287
• /
• 2001

We have studied the nonlinear evolution of a magnetized disk of isothermal gas, which is sustained by its self-gravity. Our objective is to investigate how the Jeans, Parker, and convective instabilities compete with each other in structuring/de-structuring large scale condensations in such disk. The Poisson equation for the self-gravity has been solved with a fourth-order accurate Fourier method along with the Green function, and the MHD part has been handled by an isothermal TVD code. When large wavelength perturbations are applied, the combined action of the Jeans and Parker instabilities suppresses the development of the convection and forms a dense core of prolate shape in the mid-plane. Peripheral structures around it are filamentary. The low density filaments connect the dense core to the diffuse upper region. On the other hand, when small wavelength perturbations are applied, the disk develops into an equilibrium state which is reminiscent of the Mouschovias's 2-D non-linear equilibrium of the classical Parker instability under an externally given gravity.