• The Magazine of the Society of Air-Conditioning and Refrigerating Engineers of Korea
• /
• v.7 no.1
• /
• pp.4-12
• /
• 1978

At Grashof numbers $10^{10},\;5{\times}10^{10}$, and $10^{11}$ nonlinear partial differential equations for two dimensional thermal circulation of air in a rectangular enclosure heated from below are solved numerically by finite difference explicit methood in time-dependent form. Two vertical walls and ceiling are held at low temperature and floor at high temperture. Results are compared with From's numerical solutions at $10^9{\lesssim}\;N_{Gr}\;<10^{13}$. The effective draft temperature fields are also obtained to examine cold draft problem, there included a line of constant effective draft temperature $-1.667^{\circ}C$ which is essentially Houghten's $80\%$ comfort data.

• Journal of Korean Society of Steel Construction
• /
• v.12 no.3 s.46
• /
• pp.303-315
• /
• 2000

The main object of this study is to analyze mechanical behaviors as anisotropic three-dimensional body under various static loads. This paper presents the applicability of the finite difference method to three dimensional problem of anisotropic body. The finite difference method as applied here is generalized to anisotropic three-dimensional problem of elastic body where the governing differential equations of equilibrium of such bodies are expressed in terms of the displacement u, v, and w in the coordinates axes x, y and z, care being taken to modify the finite difference expressions to satisfy the appropriate boundary conditions. By adopting a new three dimensional finite difference modelling including elimination of pivotal difference points in the case of free boundary condition, the three dimensional problem of anisotropic body was successfully completed. Several numerical results show quick convergence and numerical validity of finite difference technique in three dimensional problem.

• Journal of the Korea institute for structural maintenance and inspection
• /
• v.9 no.4
• /
• pp.159-168
• /
• 2005

This study carries out a finite difference nonlinear analysis of anisotropic advanced composite plate structures with various layer sequences. In the numerical analysis of various mechanical problems involving complex partial differential equations, the finite difference method (FDM) developed in this study has an advantage over the finite element method in its ability to avoid mesh generation and numerical integration. Many studies in FDM have been made on clamped or simple boundary conditions using merely an energy approach. These approaches cannot be satisfied, however, with pivotal points along the free boundary. Therefore, this study addresses the nonlinear problem of anisotropic plates by adopting a finite difference modeling elimination of pivotal difference points in the case of a free boundary condition. Complex nonlinear behaviors of composite plate structures for various parameters, especially for layer sequences, are analyzed using the proposed approach.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.8
• /
• pp.945-956
• /
• 1999

This work reports a set of approximate analytical solutions describing the initial transient process of close-contact melting between a rectangular parallelepiped solid and a flat plate on which either constant temperature or constant heat flux is imposed. Not only relative motion of the solid block tangential to the heating plate, but also the density difference between the solid and liquid phase is incorporated in the model. The thin film approximation reduces the force balance between the solid weight and liquid pressure, and the energy balance at the melting front into a simultaneous ordinary differential equation system. The normalized model equations admit compactly expressed analytical solutions which include the already approved two-dimensional solutions as a subset. In particular, the normalized liquid film thickness is independent of all pertinent parameters, thereby facilitating to define the transition period of close-contact melting. A unique behavior of the solid descending velocity due to the density difference is also resolved by the present solution. A new geometric function which alone represents the three-dimensional effect is introduced, and its properties are clarified. One of the representative results is that heat transfer is at least enhanced at the expense of the increase in friction as the cross-sectional shape deviates from the square under the same contact area.

• Journal of Navigation and Port Research
• /
• v.27 no.5
• /
• pp.473-478
• /
• 2003

Numerical study was made on the flow characteristics around a circular pipeline between parallel walls. The incompressible Navier-Stokes equations were solved by using a third-order upwind differential scheme. When the distance near a wall is small enough, the vortex shedding is almost completely suppressed because of the interaction with the wall boundary layer separation. This study aims to clarify the characteristics of the vortex shedding regime as the body approaches a wall as Reynolds number varies. The feature of separated vorticity dynamics is analyzed at different conditions with particular attention to the interaction between the pipeline wake and the induced separation on the plane walls.

• Journal of KSNVE
• /
• v.10 no.4
• /
• pp.610-614
• /
• 2000

Cellular Automata(CA)s are used as a simple mathematical model to investigate self-organization in statistical mechanics, which are originally introduced by von Neumann and S. Ulam at the end of the 1940s. CAs provide a framework for a large class of discrete models with homogeneous interactions, which are characterized by the following fundamental properties: 1) CAs are dynamical systems in which space and time are discrete. 2) The systems consist of a regular grid of cells. 3) Each cell is characterized by a state taken from a finite set of states and updated synchronously in discrete time steps according to a local, identical interaction rule. 4) The state of a cell is determined by the previous states of a surrounding neighborhood of cells. A cellular automaton has been attracted wide interest in modeling physical phenomena, which are described generally, partial differential equations such as diffusion and wave propagation. This paper describes one and two-dimensional analysis of wave propagation phenomena modeled by CA, where the local interaction rules were derived referring to the Lattice Gas Model reported by Chen et al., and also including finite difference scheme. Modeling processes by using CA are discussed and the simulation results of wave propagation with one wave source are compared with that by finite difference method.

• Journal of Korean Society of Steel Construction
• /
• v.12 no.4 s.47
• /
• pp.417-428
• /
• 2000

In recent years, composite materials have been used in civil engineering as well as architecture, automobile, aerospace, shipping industries. Composite materials are composed of two or more different materials to produce desirable properties for structural strength. The shell structures have the advantage of more efficient load resistance due to its curved shape as compared to the plate structures. And the shell structures with composite materials have many advantages in strength, corrosion resistance, and weight reduction. The objective of this study is to analyze circular conical shells with shear deformation effects and to prove the advantage of composite materials. To solve differential equations of conical shells, this paper used finite difference method.

• Coupled systems mechanics
• /
• v.8 no.2
• /
• pp.129-146
• /
• 2019

The main goal of this work is to develop a numerical simulator to study an isothermal single-phase two-component flow in a naturally fractured oil reservoir, taking into account advection and diffusion effects. We use the Peng-Robinson equation of state with a volume translation to evaluate the properties of the components, and the discretization of the governing partial differential equations is carried out using the Finite Difference Method, along with implicit and first-order upwind schemes. This process leads to a coupled non-linear algebraic system for the unknowns pressure and molar fractions. After a linearization and the use of an operator splitting, the Conjugate Gradient and Bi-conjugated Gradient Stabilized methods are then used to solve two algebraic subsystems, one for the pressure and another for the molar fraction. We studied the effects of fractures in both the flow field and mass transport, as well as in computing time, and the results show that the fractures affect, as expected, the flow creating a thin preferential path for the mass transport.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.12 no.3
• /
• pp.165-172
• /
• 2009

A set of equations for description of transformation of harmonic waves is proposed here. Velocity potential function and separation of variables are introduced for the derivation. The continuity equation is in a vertical plane is integrated through the water so that a horizontal one-dimensional wave equation is produced. The new equation composed of the complex velocity potential function, further be modified into. A set up of equations composed of the wave amplitude and wave phase gradient. The horizontally one-dimensional equations on the wave amplitude and wave phase gradient are the first and second-order ordinary differential equations. They are solved in a one-way marching manner starting from a side where boundary values are supplied, i.e. the wave amplitude, the wave amplitude gradient, and the wave phase gradient. Simple spatially-centered finite difference schemes are adopted for the present set of equations. The equations set is applied to three test cases, Booij's inclined plane slope profile, Massel's smooth bed profile, and Bragg's wavy bed profile. The present equations set is satisfactorily verified against existing theories including Massel's modified mild-slope equation, Berkhoff's mild-slope equation, and the full linear equation.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.17 no.10
• /
• pp.300-309
• /
• 2016

Curved beams are increasingly used in buildings, vehicles, ships, and aircraft, which has resulted in considerable effort towards developing an accurate method for analyzing the dynamic behavior of such structures. The stability behavior of elastic curved beams has been the subject of many investigations. Solutions to the relevant differential equations have traditionally been obtained by the standard finite difference or finite element methods. However, these techniques require a great deal of computer time for a large number of discrete nodes with conditions of complex geometry and loading. One efficient procedure for the solution of partial differential equations is the differential quadrature method (DQM). This method has been applied to many cases to overcome the difficulties of complex algorithms and high storage requirements for complex geometry and loading conditions. Out-of-plane buckling of curved beams with rotatory inertia were analyzed using DQM under uniformly distributed radial loads. Critical loads were calculated for the member with various parameter ratios, boundary conditions, and opening angles. The results were compared with exact results from other methods for available cases. The DQM used only a limited number of grid points and shows very good agreement with the exact results (less than 0.3% error). New results according to diverse variation are also suggested, which show important roles in the buckling behavior of curved beams and can be used for comparisons with other numerical solutions or experimental test data.