• Bulletin of the Korean Mathematical Society
• /
• v.31 no.2
• /
• pp.309-318
• /
• 1994

The approach presented in this paper is based on the transformation of the Stefan problem in one space dimension to an initial-boundary value problem for the heat equation in a fixed domain. Of course, the problem is non-linear. The finite element approximation adopted here is the standared continuous Galerkin method in time. In this paper, only the regular case is discussed. This means the error analysis is based on the assumption that the solution is sufficiently smooth. The aim of this paper is the existence of the solution in a finite Galerkin system of ordinary equations.

• Journal of applied mathematics & informatics
• /
• v.28 no.3_4
• /
• pp.651-661
• /
• 2010

We consider the inverse problem of the identification of a piecewise-constant conductivity in a bar given the extra information of the heat flux through one end of the bar. Our theoretical results show that such an identification is unique. This approach utilizes a "layer peeling" argument. A computational algorithm based on this method is proposed and implemented. The advantage of this algorithm is that it requires only 3D minimizations irrespective of the number of the unknown discontinuities. Its numerical effectiveness is investigated for several conductivities.

• The Pure and Applied Mathematics
• /
• v.1 no.1
• /
• pp.29-33
• /
• 1994

There are various aspects of the solution of boundary-value problems for second-order linear elliptic equations in two independent variables. One useful method of solving such boundary-value problems for Laplace's equation is by means of suitable integral representations of solutions and these representations are obtained most directly in terms of particular singular solutions, termed Green's functions.(omitted)

• Geotechnical Engineering
• /
• v.4 no.4
• /
• pp.19-28
• /
• 1988

An underground explosion crests shock waves, which propagate to a buried structure through the이 ground. Due to the explosion, very high stresses and large deformation occur in the ground so that the shock waves decay gradually. In this study the numerical simulation of the ground shock attenuation has teen performed. One dimensional wave equation is presented and the finite difference method is applies. A Cap model is adopted to describe the stress-strain behavior of the ground. The results are expressed by the attenuation of the peak stress and the particle vrelocity by the time and the distance.

• Korean Journal of Radiology
• /
• v.24 no.7
• /
• pp.601-605
• /
• 2023

• Structural Engineering and Mechanics
• /
• v.20 no.2
• /
• pp.191-207
• /
• 2005

Time delay exists inevitably in active control, which may not only degrade the system performance but also render instability to the dynamic system. In this paper, a novel active controller is developed to solve the time delay problem in flexible structures. By using the independent modal space control method, the differential equation of the controlled mode with time delay is obtained from the time-delay system dynamics. Then it is discretized and changed into a first-order difference equation without any explicit time delay by augmenting the state variables. The modal controller is derived based on the augmented system using the discrete variable structure control method. The switching surface is determined by minimizing a discrete quadratic performance index. The modal coordinate is extracted from sensor measurements and the actuator control force is converted from the modal one. Since the time delay is explicitly included throughout the entire controller design without any approximation, the system performance and stability are guaranteed. Numerical simulations show that the proposed controller is feasible and effective in active vibration control of dynamic systems with time delay. If the time delay is not explicitly included in the controller design, instability may occur.

• 제어로봇시스템학회:학술대회논문집
• /
• 2003.10a
• /
• pp.1798-1803
• /
• 2003

In this work the dynamic heat transfer occurring in a cable penetration fire stop system built in the firewall of nuclear power plants is three-dimensionally investigated to develop a test-simulator that can be used to verify effectiveness of the sealants. The dynamic heat transfer can be described by a partial differential equation (PDE) and its initial and boundary conditions. For the shake of simplicity PDE is divided into two parts; one corresponding to the heat transfer in the axial direction and the other corresponding to the heat transfer on the vertical layers. Two numerical methods, SOR (Sequential Over-Relaxation) and FEM (Finite Element Method), are implemented to solve these equations respectively. The axial line is discretized, and SOR is applied. Similarly, all the layers are separated into finite elements, where the time and spatial functions are assumed to be of orthogonal collocation state at each element. The heat fluxes on the layers are calculated by FEM. It is shown that the penetration cable influences the temperature distribution of the fire stop system very significantly. The simulation results are shown in the three-dimensional graphics for the understanding of the transient temperature distribution in the fire stop system.

• Journal of the Korean Institute of Illuminating and Electrical Installation Engineers
• /
• v.19 no.6
• /
• pp.29-35
• /
• 2005

This paper presents the basic quench protection idea for the HTS(High-Temperature Superconducting) cable. In Korea power system, the transfer capability of transmission line is limited by the voltage stability, HTS cable could be one of the countermeasure to enhance the transfer limit with its higher current capacity and lower impedance[1]. However, the quench characteristic makes not only HTS cable to loss its superconductivity, but also change the impedance of the transmission line and power system operating condition dramatically. This pheonominum threats HTS cable safety as well as power system security, therefore a proper protection scheme and security control counterplan have to be established before HTS cable implementation. In this paper, the quench characteristics of HTS cable for the fault current based on heat balance equation was established and a proper protection method regarding conventional protection system was suggested.

• 한국전산유체공학회:학술대회논문집
• /
• 2001.05a
• /
• pp.189-194
• /
• 2001

In this paper, non-reacting and reacting flowfields were computed using a preconditioned Navier-Stokes solver. The preconditioning technique of Merkle et al. and TVD scheme or Chakravarthy and Osher was employed and the results obtained using developed code have a good agreement with the previous results and experimental data. The preconditioned Wavier-Stokes equation set with low Reynolds number $\kappa-\epsilon$ equation and species continuity equations, are discretized with strongly implicit manner and time integrated with LU-SSOR scheme. For the purpose of treating unsteady problem the duel-time stepping scheme was employed. For the validation of the code in incompressible flow regime, steady driven square cavity flow was considered and calculation result shows reasonably good agreement with the result of incompressible code. Shock wave/boundary layer interaction problem was considered to show the shock capturing performance of preconditioned-TVD scheme. To validate unsteady flow, acoustic oscillation problem was calculated, and supersonic premix flame of $H_2$-air reaction problem which is calculated with turbulence model, 9-species/18-reaction step reaction model, shows reasonable agreement with the previous results. As a result, the preconditioning method has an advantage to calculate incompressible and compressible flow through one code and preconditioned solver easily developed from standard compressible code with minor efforts. But additional computational time and computer memory is required due to preconditioning matrix.

• Journal of KSNVE
• /
• v.6 no.5
• /
• pp.607-616
• /
• 1996

A combined computational fluid dynamics(CFD)-Kirchhoff method is presented for predicting high-speed impulsive noise generated by a hovering blade. Two types of Kirchhoff integral formula are used; one for the classical linear Kirchhoff formulation and the other for the nonlinear Kirchhoff formulation. An Euler finite difference solver is solved first to obtain the flow field close to the blade, and then this flow field is used as an input to a Kirchhoff formulation to predict the acoustic far-field. These formulas are used at Mach numbers of 0.90 and 0.95 to investigate the effectiveness of the linear and nonlinear Kirchhoff formulas for delocalized flow. During these calculiations, the retarded time equation is also carefully examined, in particular, for the cases of the control surface located outside of the sonic cylinder, where multiple roots are obtained. Predicted results of acoustic far-field pressure with the linear Kirchhoff formulation agree well with experimental data when the control surface is at the certain location(R=1.46), but the correlation is getting worse before or after this specific location of the control surface due to the delocalized nonlinear aerodynamic flow field. Calculations based on the nonlinear Kirchhoff equation using a linear sonic cylinder as a control surface show a reasonable agreement with experimental data in negative amplitudes for both tip Mach numbers of 0.90 and 0.95, except some computational integration problems over a shock. This concliudes that a nonlinear formulation is necessary if the control surface is close to the blade and the flow is delocalized.