• Nuclear Engineering and Technology
• /
• v.52 no.2
• /
• pp.230-237
• /
• 2020

In this part, an implicit time dependent solution is presented for the Boltzmann transport equation discretized by the analytic coarse mesh finite difference method (ACMFD) over the spatial domain as well as the simplified P₃ (SP₃) for the angular variable. In the first part of this work we proposed a SP₃-ACMFD approach to solve the static eigenvalue equations which provide the initial conditions for temp dependent equations. Having solved the 3D multi-group SP₃-ACMFD static equations, an implicit approach is resorted to ensure stability of time steps. An exponential behavior is assumed in transverse integrated equations to establish a relationship between flux moments and currents. Also, analytic integration is benefited for the time-dependent solution of precursor concentration equations. Finally, a multi-channel one-phase thermal hydraulic model is coupled to the proposed methodology. Transient equations are then solved at each step using the GMRES technique. To show the sufficiency of proposed transient SP₃-ACMFD approximation for a full core analysis, a comparison is made using transport peers as the reference. To further demonstrate superiority, results are compared with a 3D multi-group transient diffusion solver developed as a byproduct of this work. Outcomes confirm that the idea can be considered as an economic interim approach which is superior to the diffusion approximation, and comparable with transport in results.

• Journal of Korean Society of Steel Construction
• /
• v.18 no.5
• /
• pp.575-584
• /
• 2006

In this paper, a shear-flexible finite element model is developed for the buckling analysis of axially loaded, thin-walled composite I-beams. Based on an orthogonal Cartesian coordinate system, the displacement fields are defined using the first-order shear-deformable beam theory. The derived element takes into account flexural shear deformation and torsional warping deformation. Three different types of beam elements, namely, the two-noded, three-noded, and four-noded beam elements, were developed to solve the governing equations. An inverse iteration with shift eigenvalue solution was used to solve the resulting linearized buckling problem. A parametric study was conducted to show the importance of shear flexibility and fiber orientation on the buckling behavior of thin-walled composite beams. A good agreement was obtained among the proposed shear-flexible model, other results available in literature, and the finite element solution.

• Journal of Korean Society of Steel Construction
• /
• v.14 no.1
• /
• pp.23-29
• /
• 2002

Many of the bridge systems, including the girders and cross-beams, and concrete decks behave as the special orthotropic plates. Such systems with boundary conditions other than Navier or Levy solution types, or with irregular cross sections, analytical solution is very difficult to obtain. Numerical method for eigenvalue problems are also very much involved in seeking such a solution. A method of calculating the natural frequency corresponding to the first mode of vibration of beam and tower structures with irregular cross-sections was developed and reported by Kim in 1974. Recently, this method was extended to two dimensional problems including composite laminates, and has been applied to composite plates with shear deformation effects. In this paper, application of this method to the specially orthotropic laminated plates with various boundary condition is accomplished and the result of analysis is presented.

• The Bulletin of The Korean Astronomical Society
• /
• v.43 no.1
• /
• pp.53.1-53.1
• /
• 2018

During a major solar eruption, the erupting plasma is possibly out of the equilibrium ionization state because of its rapid heating or cooling. The non-equilibrium ionization process is important in a rapidly evolving system where the thermodynamical time scale is shorter than the ionization or recombination time scales. We investigate the effects of non-equilibrium ionization on EUV and X-ray observations by the Atmospheric Imaging Assembly (AIA) on board Solar Dynamic Observatory and X-ray Telescope (XRT) on board Hinode. For the investigation, first, we find the emissivities for all the lines of ions of elements using CHIANTI 8.07, and then we find the temperature responses multiplying the emissivities by the effective area for each AIA and XRT passband. Second, we obtain the ion fractions using a time-dependent ionization model (Shen et al. 2015), which uses an eigenvalue method, for all the lines of ion, as a function of temperature, and a characteristic time scale, $n_et$, where $n_e$ and t are density and time, respectively. Lastly, the ion fractions are multiplied to the temperature response for each passband, which results in a 2D grid for each combination of temperature and the characteristic time scale. This is the set of passband responses for plasma that is rapidly ionized in a current sheet or a shock. We investigate an observed event which has a relatively large uncertainty in an analysis using a differential emission measure method assuming equilibrium ionization state. We verify whether the observed coronal plasmas are in non-equilibrium or equilibrium ionization state using the passband responses.

• Journal of Ocean Engineering and Technology
• /
• v.23 no.1
• /
• pp.74-80
• /
• 2009

In this study, the static and modal analyses to find the characteristic of eigenvalues for a towed cable were with a free boundary condition at the bottom end carried out with numerical study. The resulting numerical code with finite element method was used to study sample problems for a cable with towing speeds. After tracing the equilibrium state with a towing speed through the static analysis, modal analysis on the basis of static results was performed. The static top tension for a critical towing speed is nearly 50 percent of what it was for a free hanging pipe. From static analyses, it is found that towing speed has a noticeable effect on top tension of a towed pipe. At a high towing speed, differences between the first and second periods become larger. Compared to the fundamental period for a free hanging pipe, that for a towed pipe with a critical towing speed is approximately 1.4 times larger. This result is very important point in that the lock in condition and tension of the towed cable system with top excitation can be predicted. The corrected close form solution to solve natural periods for a towed cable was presented in this study. The code is validated by comparison of the results of theoretical and numerical studies. Two results were in very good agreement. This study can contribute to predicting the lock-in condition and tension for a towed cable or pipe with top excitation.

• Earthquakes and Structures
• /
• v.14 no.2
• /
• pp.103-115
• /
• 2018

In this work, a simple but accurate hyperbolic plate theory for the free vibration analysis of functionally graded material (FGM) sandwich plates is developed. The significant feature of this formulation is that, in addition to including the shear deformation effect, it deals with only 3 unknowns as the classical plate theory (CPT), instead of 5 as in the well-known first shear deformation theory (FSDT) and higher-order shear deformation theory (HSDT). A shear correction factor is, therefore, not required. Two common types of FGM sandwich plates are considered, namely, the sandwich with the FGM face sheet and the homogeneous core and the sandwich with the homogeneous face sheet and the FGM core. The equation of motion for the FGM sandwich plates is obtained based on Hamilton's principle. The closed form solutions are obtained by using the Navier technique. The fundamental frequencies are found by solving the eigenvalue problems. Numerical results of the present theory are compared with the CPT, FSDT, order shear deformation theories (HSDTs), and 3D solutions. Verification studies show that the proposed theory is not only accurate and simple in solving the free vibration behaviour of FGM sandwich plates, but also comparable with the higher-order shear deformation theories which contain more number of unknowns.

• Journal of the Korean Mathematical Society
• /
• v.15 no.1
• /
• pp.39-57
• /
• 1978

It is given in the Lecture Note [1] of Berger, Gauduchon and Mazet that, if ${\pi}$n: (${\tilde{M}}$, ${\tilde{g}}$)${\rightarrow}$(${\tilde{M}}$, ${\tilde{g}}$) is a Riemannian submersion with totally geodesic fibers, ${\Delta}$ and ${\tilde{\Delta}}$ are Laplace operators on (${\tilde{M}}$, ${\tilde{g}}$) and (M, g) respectively and f is an eigenfunction of ${\Delta}$, then its lift $f^L$ is also an eigenfunction of ${\tilde{\Delta}}$ with the common eigenvalue. But such a simple relation does not hold for an eigenform of the Laplace-Beltrami operator ${\Delta}=d{\delta}+{\delta}d$. If ${\omega}$ is an eigenform of ${\Delta}$ and ${\omega}^L$ is the horizontal lift of ${\omega}$, ${\omega}^L$ is not in genera an eigenform of the Laplace-Beltrami operator ${\tilde{\Delta}}$ of (${\tilde{M}}$, ${\tilde{g}}$). The present author has obtained a set of formulas which gives the relation between ${\tilde{\Delta}}{\omega}^L$ and ${\Delta}{\omega}$ in another paper [7]. In the present paper a Sasakian submersion is treated. A Sasakian manifold (${\tilde{M}}$, ${\tilde{g}}$, ${\tilde{\xi}}$) considered in this paper is such a one which admits a Riemannian submersion where the base manifold is a Kaehler manifold (M, g, J) and the fibers are geodesics generated by the unit Killing vector field ${\tilde{\xi}}$. Then the submersion is called a Sasakian submersion. If ${\omega}$ is a eigenform of ${\Delta}$ on (M, g, J) and its lift ${\omega}^L$ is an eigenform of ${\tilde{\Delta}}$ on (${\tilde{M}}$, ${\tilde{g}}$, ${\tilde{\xi}}$), then ${\omega}$ is called an eigenform of the first kind. We define a relative eigenform of ${\tilde{\Delta}}$. If the lift ${\omega}^L$ of an eigenform ${\omega}$ of ${\Delta}$ is a relative eigenform of ${\tilde{\Delta}}$ we call ${\omega}$ an eigenform of the second kind. Such objects are studied.

• Journal of Advanced Navigation Technology
• /
• v.22 no.3
• /
• pp.228-232
• /
• 2018

Dealing with quntization error in a control system properly becomes much more important as many devices are connected through network and controlled. Thus, in this paper, we study a data rate condition on quantizer to achieve practical stability in a discrete time linear time invariant system with state feedback control. First, required data rate is shown to depend on eigenvalue of the closed loop system, the size of the initial state vector, the magnitude of initial quantization error, and control gain in the absence of process noise. It additionally depends on the maximum magnitude of process noise when noise is not zero. Asymptotic analysis shows that a new design method may be needed to reduce the date rate for a networked control in the presence of quantization error and noise.. We provide a simple numerical evaluation of uniform quantizer and logarithmic qunatizer to assess their characteristics of practical stability depending on data rate in the presence of noise.

• Nuclear Engineering and Technology
• /
• v.46 no.6
• /
• pp.783-796
• /
• 2014

The coarse mesh finite difference (CMFD) method is applied to the discontinuous finite element method based discrete ordinate calculation for source convergence acceleration. The three-dimensional (3-D) DFEM-Sn code FEDONA is developed for general geometry applications as a framework for the CMFD implementation. Detailed methods for applying the CMFD acceleration are established, such as the method to acquire the coarse mesh flux and current by combining unstructured tetrahedron elements to rectangular coarse mesh geometry, and the alternating calculation method to exchange the updated flux information between the CMFD and DFEM-Sn. The partial current based CMFD (p-CMFD) is also implemented for comparison of the acceleration performance. The modified p-CMFD method is proposed to correct the weakness of the original p-CMFD formulation. The performance of CMFD acceleration is examined first for simple two-dimensional multigroup problems to investigate the effect of the problem and coarse mesh sizes. It is shown that smaller coarse meshes are more effective in the CMFD acceleration and the modified p-CMFD has similar effectiveness as the standard CMFD. The effectiveness of CMFD acceleration is then assessed for three-dimensional benchmark problems such as the IAEA (International Atomic Energy Agency) and C5G7MOX problems. It is demonstrated that a sufficiently converged solution is obtained within 7 outer iterations which would require 175 iterations with the normal DFEM-Sn calculations for the IAEA problem. It is claimed that the CMFD accelerated DFEM-Sn method can be effectively used in the practical eigenvalue calculations involving general geometries.

• The Journal of the Institute of Internet, Broadcasting and Communication
• /
• v.20 no.5
• /
• pp.151-156
• /
• 2020

To analyze the diffraction properties of optical signals and the significant impacts of tapered sidewall profile at periodic trapezoidal 2D diffraction gratings, Toeplitz dielectric tensor is first defined and formulated by 2D spatial Fourier expansions associated with trapezoidal profile. The characteristic modes in each layer is then based on eigenvalue problem, and the complete solution is found rigorously in terms of modal transmission-line theory (MTLT) to address the pertinent boundary-value problems. Based on those one, the numerical analysis is performed on how the tapered side profile of grating structures with trapezoidal refractive index distribution affects the design of a sub-wavelength grating reflector. The numerical results reveal that this tapered sidewall profile plays a critical role in determining the reflection bandwidth, the average reflectance, and the band edge.