• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.28 no.2
• /
• pp.123-130
• /
• 2015

This paper presents a new formulation for transient simulations of microwave propagation in heterogeneous unbounded domains. In particular, perfectly-matched-layers(PMLs) are introduced to allow for wave absorption at artificial boundaries used to truncate the infinite extent of the physical domains. The development of the electromagnetic PML targets the application to engineering mechanics problems such as structural health monitoring and inverse medium problems. To formulate the PML for plane electromagnetic waves, a complex coordinate transformation is introduced to Maxwell's equations in the frequency-domain. Then the PML-endowed partial differential equations(PDEs) for transient electromagnetic waves are recovered by the application of the inverse Fourier transform to the frequency-domain equations. A mixed finite element method is employed to solve the time-domain PDEs for electric and magnetic fields in the PML-truncated domain. Numerical results are presented for plane microwaves propagating through concrete structures, and the accuracy of solutions is investigated by a series of error analyses.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.9
• /
• pp.56-65
• /
• 2005

A dynamic analysis method that freezes a time domain by discretization and solves the spatial propagation equation has a unique feature that provides a degree of freedom on spatial domain compared with the space discretization or space-time discretization finite element method. Using this feature, the time finite element analysis can be effectively applied to optimize the spatial characteristics of distributed type actuators. In this research, the time domain finite element method was used to discretize the model. A state variable vector was used in the discretization to include arbitrary initial conditions. A performance index was proposed on spatial domain to consider both potential and vibrational energy, so that the resulting shape of the distributed actuator was optimized for dynamic control of the structure. It is assumed that the structure satisfies the final rest condition using the realizable control scheme although the initial disturbance can affect the system response. Both equations on states and costates were derived based on the selected performance index and structural model. Ricatti matrix differential equations on state and costate variables were derived by the reconfiguration of the sub-matrices and application of time/space boundary conditions, and finally optimal actuator distribution was obtained. Numerical simulation results validated the proposed actuator shape optimization scheme.

• Communications of the Korean Mathematical Society
• /
• v.18 no.3
• /
• pp.569-580
• /
• 2003

The Lotka-McKendrick equation which describes the evolution of a single population under the phenomenological conditions is developed from the well-known Malthus’model. In this paper, we introduce the Lotka-McKendrick equation for the description of the dynamics of a population. We apply a discontinuous Galerkin finite element method in age-time domain to approximate the solution of the system. We provide some numerical results. It is experimentally shown that, when the mortality function is bounded, the scheme converges at the rate of $h^2$ in the case of piecewise linear polynomial space. It is also shown that the scheme converges at the rate of $h^{3/2}$ when the mortality function is unbounded.

• Journal of applied mathematics & informatics
• /
• v.35 no.3_4
• /
• pp.277-302
• /
• 2017

In this paper, we have proposed a numerical method for Singularly Perturbed Boundary Value Problems (SPBVPs) of reaction-diffusion type of third order Ordinary Differential Equations (ODEs). The SPBVP is reduced into a weakly coupled system of one first order and one second order ODEs, one without the parameter and the other with the parameter ${\varepsilon}$ multiplying the highest derivative subject to suitable initial and boundary conditions, respectively. The numerical method combines boundary value technique, asymptotic expansion approximation, shooting method and finite difference scheme. The weakly coupled system is decoupled by replacing one of the unknowns by its zero-order asymptotic expansion. Finally the present numerical method is applied to the decoupled system. In order to get a numerical solution for the derivative of the solution, the domain is divided into three regions namely two inner regions and one outer region. The Shooting method is applied to two inner regions whereas for the outer region, standard finite difference (FD) scheme is applied. Necessary error estimates are derived for the method. Computational efficiency and accuracy are verified through numerical examples. The method is easy to implement and suitable for parallel computing. The main advantage of this method is that due to decoupling the system, the computation time is very much reduced.

• Steel and Composite Structures
• /
• v.19 no.6
• /
• pp.1421-1447
• /
• 2015

In this paper, the effect of material-temperature dependent on the wave propagation of a cantilever beam composed of functionally graded material (FGM) under the effect of an impact force is investigated. The beam is excited by a transverse triangular force impulse modulated by a harmonic motion. Material properties of the beam are temperature-dependent and change in the thickness direction. The Kelvin-Voigt model for the material of the beam is used. The considered problem is investigated within the Euler-Bernoulli beam theory by using energy based finite element method. The system of equations of motion is derived by using Lagrange's equations. The obtained system of linear differential equations is reduced to a linear algebraic equation system and solved in the time domain and frequency domain by using Newmark average acceleration method. In order to establish the accuracy of the present formulation and results, the comparison study is performed with the published results available in the literature. Good agreement is observed. In the study, the effects of material distributions and temperature rising on the wave propagation of the FGM beam are investigated in detail.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.13 no.1
• /
• pp.37-49
• /
• 2000

This paper is for the first essential study on the development of unified finite element formulations for solving problems related to the dynamics/thermoelastics behavior of solids. In the first part of formulations, the finite element method is based on the introduction of a new quantity defined as heat displacement, which allows the heat conduction equations to be written in a form equivalent to the equation of motion, and the equations of coupled thermoelasticity to be written in a unified form. The equations obtained are used to express a variational formulation which, together with the concept of generalized coordinates, yields a set of differential equations with the time as an independent variable. Using the Laplace transform, the resulting finite element equations are described in the transform domain. In the second, the Laplace transform is applied to both the equation of heat conduction derived in the first part and the equations of motions and their corresponding boundary conditions, which is referred to the transformed equation. Selections of interpolation functions dependent on only the space variable and an application of the weighted residual method to the coupled equation result in the necessary finite element matrices in the transformed domain. Finally, to prove the validity of two approaches, a comparison with one finite element equation and the other is made term by term.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.25 no.2
• /
• pp.54-65
• /
• 2021

In this study, a two-dimensional thermoelastic problem under hyperbolic heat conduction theory with an internal heat source is considered. The general solution for the temperature field, stress components and displacement field are obtained using the reduced differential transform method. The stress and displacement components are obtained using the thermal stress function in the reduced differential transform domain. All the solutions are obtained in the form of power series. The special case with a time-dependent laser heat source has been considered. The problem is considered for homogeneous material with finite rectangular cross-section heated with a non-Gaussian temporal profile. The effect of the heat source on all the characteristics of a material is discussed numerically and graphically for magnesium material taking a pulse duration of 0.2 ps. This study provides a powerful tool for finding the solution to the thermoelastic problem with less computational work as compared to other methods. The result obtained in the study may be useful for the investigation of thermal characteristics in engineering and industrial applications.

• The Transactions of the Korean Institute of Electrical Engineers
• /
• v.37 no.12
• /
• pp.839-846
• /
• 1988

If we use the 2-dimensional analyzing method with the magnetic vector potential for the analysis of eddy current distribution in electric machinery, we can obtain the magnitude of eddy current but can't have the characteristic of eddy current distribution. For the settlement of this problem, we have induced the governing equation with the current vector potential and attemptted 2-dimensional analysis of eddy current distribution by finite element method. And the time domain weighted residual method is used in treatment of time differential term and we have developed the algorithm by it. And then, we analyze eddy current distributions of analytic model and aluminium disk in singlephase watt hour meter. Consequently we have verified the propriety and utility of above mentioned method.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.26 no.4
• /
• pp.223-231
• /
• 2013

This paper presents a time-domain Gauss-Newton full-waveform inversion method for the material profile reconstruction in heterogeneous semi-infinite solid media. To implement the inverse problem in a finite computational domain, perfectly-matchedlayers( PMLs) are introduced as wave-absorbing boundaries within which the domain's wave velocity profile is to be reconstructed. The inverse problem is formulated in a partial-differential-equations(PDE)-constrained optimization framework, where a least-squares misfit between measured and calculated surface responses is minimized under the constraint of PML-endowed wave equations. A Gauss-Newton-Krylov optimization algorithm is utilized to iteratively update the unknown wave velocity profile with the aid of a specialized regularization scheme. Through a series of one-dimensional examples, the solution of the Gauss-Newton inversion was close enough to the target profile, and showed superior convergence behavior with reduced wall-clock time of implementation compared to a conventional inversion using Fletcher-Reeves optimization algorithm.

• Journal of the Institute of Electronics Engineers of Korea SD
• /
• v.41 no.3
• /
• pp.45-54
• /
• 2004

Microwave characteristics of ridge type CPW traveling-wave(TW) electroabsorption modulator and photodetector are affected by the thickness of intrinsic layer, width of guiding layer, and the separation of signal and ground electrodes. These factors are determined effective index of microwave and characteristic impedance due to changing of capacitance(C) and inductance(L) of device. However, conventional equivalent circuit of TW-structure is approximated to microstrip and CPW transmission line by distribution of electric and magnetic fields, respectively. In this paper, we analyzed microwave characteristics of TW-structure and found more accurate value of C and L by using finite difference time domain (FDTD) method. These values are adopted circuit element of equivalent circuit. Microwave characteristics obtained by the FDTD and equivalent circuit model show good agreement.