• Structural Engineering and Mechanics
• /
• v.69 no.1
• /
• pp.33-42
• /
• 2019

Thermal cracks are cracks that commonly form at early ages in mass concrete. During the concrete pouring process, the elastic modulus changes continuously. This requires the time domain to be divided into several steps in order to solve for the temperature, stress, and displacement of the concrete. Numerical simulations of thermal crack propagation in concrete are more difficult at early ages. To solve this problem, this study divides crack propagation in concrete at early ages into two cases: the case in which cracks do not propagate but the elastic modulus of the concrete changes and the case in which cracks propagate at a certain time. This paper provides computational models for these two cases by integrating the characteristics of the extended finite element algorithm, compiles the corresponding computational programs, and verifies the accuracy of the proposed model using numerical comparisons. The model presented in this paper has the advantages of high computational accuracy and stable results in resolving thermal cracking and its propagation in concrete at early ages.

• Journal of Communications and Networks
• /
• v.12 no.6
• /
• pp.605-615
• /
• 2010

Intercarrier interference (ICI) in orthogonal frequency division multiplexing (OFDM) systems, which causes substantial performance degradation in time-varying fading channels, is analyzed. An equivalent spreading code formulation is derived based on the analogy of OFDM and code division multiple access (CDMA) systems. Techniques as linear multiuser detection in CDMA systems are applied to suppress the ICI in OFDM systems. The performance of linear detection, measured using multiuser efficiency and asymptotic multiuser efficiency, is analyzed given the assumption of perfect channel state information (CSI), which serves as an upper bound for the performance of practical systems. For systems without CSI, time domain and frequency domain channel estimation based linear detectors are proposed. The performance gains and robustness of a linear minimum mean square error (LMMSE) filter over a traditional filter (TF) and matched filter (MF) in the high signal-to-noise ratio (SNR) regime are demonstrated with numerical simulation results.

• The Journal of the Institute of Internet, Broadcasting and Communication
• /
• v.9 no.5
• /
• pp.139-146
• /
• 2009

Parametric array for audio applications is analyzed by numerical modeling and analytical approximation. The nonlinear wave equations are used to provide design guidelines for the audio parametric array. A time domain finite difference code that accurately solves the KZK (Khokhlov-Zabolotskaya-Kuznetsov) nonlinear parabolic wave equation is used to predict the response of the parametric array. The time domain code relates the source size and the carrier frequency to the audible signal response including the output level and beamwidth to considering the implementation issues for audio applications of the parametric array, the emphasis is given to the frequency response and distortion. We use the time domain code to find out the optimal parameters that will help produce the parametric array with highest achievable output in terms of the average power within the demodulated signal. Parameters such as primary input frequency, audio source radius and the modulation method are given utmost importance. The output effect of those parameters are demonstrated through the numerical simulation.

• Journal of Mechanical Science and Technology
• /
• v.18 no.3
• /
• pp.460-470
• /
• 2004

A computational analysis of engineering problems with moving domain or/and boundary according to either Lagrangian or Eulerian approach may encounter inherent numerical difficulties, the extreme mesh distortion in the former and the material boundary indistinctness in the latter. In order to overcome such defects in classical numerical approaches, the ALE(arbitrary Lagrangian Eulerian) method is widely being adopted in which the finite element mesh moves with arbitrary velocity. This paper is concerned with the ALE finite element formulation, aiming at the dynamic response analysis of baffled fuel-storage container in yawing motion, for which the coupled time integration scheme, the remeshing and smoothing algorithm and the mesh velocity determination are addressed. Numerical simulation illustrating theoretical works is also presented.

• Journal of the Society of Naval Architects of Korea
• /
• v.41 no.1
• /
• pp.23-30
• /
• 2004

In the Lagrangian vortex particle method based on the vorticity-velocity formulation for solving the incompressible Navier-Stokes equations, a numerical scheme for calculating pressure fields is presented. Implementation of the numerical method is directly connected with the well-established surface panel methods, just by dealing with the dynamic coupling among vorticity field. Assuming the vorticity and the velocity fields are to be calculated in time domain analysis, the pressure calculation for a complete set of solution at present time step is performed in a similar way to the one used in the Eulerian description. For a validation of the present method, we illustrate the early development of the viscous flow about an impulsive started circular cylinder for Reynolds number 550. The comparative study with the Eulerian finite Volume method provides an extensive understanding and application of the mesh-free Lagrangian vortex methods for numerical simulation of viscous flows around arbitrary bodies of general shape.

• Journal of Mechanical Science and Technology
• /
• v.18 no.4
• /
• pp.718-730
• /
• 2004

The solution for the natural convection in internally finned horizontal annuli is obtained by using a numerical simulation of time-dependent and two-dimensional governing equations. The fins existing in annuli influence the flow pattern, temperature distribution and heat transfer rate. The variations of the On configuration suppress or accelerate the free convective effects compared to those of the smooth tubes. The effects of fin configuration, number of fins and ratio of annulus gap width to the inner cylinder radius on the fluid flow and heat transfer in annuli are demonstrated by the distribution of the velocity vector, isotherms and streamlines. The governing equations are solved efficiently by using a parallel implementation. The technique is adopted for reduction of the computation cost. The parallelization is performed with the domain decomposition technique and message passing between sub-domains on the basis of the MPI library. The results from parallel computation reveal in consistency with those of the sequential program. Moreover, the speed-up ratio shows linearity with the number of processor.

• International Journal of Aeronautical and Space Sciences
• /
• v.3 no.2
• /
• pp.76-81
• /
• 2002

A numerical simulation of wake behavior behind three-dimensional wings in ground effect is done using an indirect boundary element method (Panel Method). An integral equation is obtained by applying Green's 2nd Identity on all surfaces of the flow domain. The AIC is constructed by imposing the no penetration condition on solid surfaces, and the Kutta at the wing's trailing edge. The ground effect is included using an image method. At each time step, a row of wake panels from wings' trailing edge are convected downstream following the force-free condition. The roll-up of wake vortices behind wings in close proximity is simulated.

• Proceedings of the Optical Society of Korea Conference
• /
• 2003.02a
• /
• pp.184-185
• /
• 2003

Unlike the conventional waveguide such as optical fiber using total internal reflection, photonic crystal waveguide(PCW), a waveguide made of a line defect in a photonic crystal(PC) structure, does not admit an analytic approach due to its complexity but requires a direct numerical approach. Here, we present numerical results of computer simulation for PCW by using the three-dimensional(3D) Finite-Difference Time -Domain(FDTD) algorithm. (omitted)

• Progress in Superconductivity and Cryogenics
• /
• v.21 no.2
• /
• pp.45-49
• /
• 2019

This paper reports comparison between analytic and numerical simulation approaches for calculation of screening current and screening current induced field in a high temperature superconductor magnet. Bean slab model is adopted to calculate screening current and SCF analytically, while the finite element method numerically. A case study of screening current and SCF calculation are conducted with a magnet, a 7 T 68 mm cold-bore multi-width no-insulation GdBCO magnet built and tested by Massachusetts Institute of Technology Francis Bitter Magnet Laboratory. In this study, we assume the magnet is dunked in liquid nitrogen at 77 K. Furthermore, the simulation results are compared in terms of computation time and accuracy. Finally, discussion on the different methods together with the comparison between the calculations and experiment is provided.

• Earthquakes and Structures
• /
• v.10 no.5
• /
• pp.989-1012
• /
• 2016

The structural dynamic behavior and yield strength considering both ductility and strain rate effects are analyzed in this article. For the single-degree-of-freedom (SDOF) system, the relationship between the relative velocity and the strain rate response is deduced and the strain rate spectrum is presented. The ductility factor can be incorporated into the strain rate spectrum conveniently based on the constant-ductility velocity response spectrum. With the application of strain rate spectrum, it is convenient to consider the ductility and strain rate effects in engineering practice. The modal combination method, i.e., square root of the sum of the squares (SRSS) method, is employed to calculate the maximum strain rate of the elastoplastic multiple-degree-of-freedom (MDOF) system under uniform excitation. Considering the spatially varying ground motions, a new response spectrum method is developed by incorporating the ductility factor and strain rate into the conventional response spectrum method. In order to further analyze the effects of strain rate and ductility on structural dynamic behavior and yield strength, the cantilever beam (one-dimensional) and the triangular element (two-dimensional) are taken as numerical examples to calculate their seismic responses in time domain. Numerical results show that the permanent displacements with and without considering the strain rate effect are significantly different from each other. It is not only necessary in theory but also significant in engineering practice to take the ductility and strain rate effects into consideration.