• Journal of the Korean Mathematical Society
• /
• v.32 no.3
• /
• pp.541-551
• /
• 1995

The nonlinear Klein Gordon equation $$ (1) \frac{\partial t^2}{\partial^2 u} - \Delta u + V_u(u) = f $$ where $\Delta$ is the Laplacian operator in $R^d (d = 1, 2, 3), V_u(u)$ is the derivative of the "potential function" V, and f is a source term independent of the solution u, in various areas of mathematical physics.l physics.

• Journal of the Korean Mathematical Society
• /
• v.50 no.6
• /
• pp.1291-1310
• /
• 2013

This paper develops least-squares pseudo-spectral collocation methods for elliptic boundary value problems having interface conditions given by discontinuous coefficients and singular source term. From the discontinuities of coefficients and singular source term, we derive the interface conditions and then we impose such interface conditions to solution spaces. We define two types of discrete least-squares functionals summing discontinuous spectral norms of the residual equations over two sub-domains. In this paper, we show that the homogeneous least-squares functionals are equivalent to appropriate product norms and the proposed methods have the spectral convergence. Finally, we present some numerical results to provide evidences for analysis and spectral convergence of the proposed methods.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.22 no.3
• /
• pp.360-372
• /
• 1998

This study deals with the interaction between buoyancy-induced convection and externally imposed excitation in the form of harmonic rocking and the effect of the interaction upon heat transfer in low-Pr liquids. A wide array of system responses are discussed using the spectral collocation numerical technique. The superposition of buoyancy and Coriolis forces leads to complex fluid flow and heat transfer. The transition to chaotic convection is accelerated, and heat transfer rates are reduced as the enclosure is excited at the fundamental frequency of oscillation associated with the pure buoyancy-driven case. Average heat transfer rates are correlated for Pr=0.02 and 0.03. The heat transfer is affected more in the Pr=0.03 liquid than the case of Pr=0.02.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.24 no.1 s.173
• /
• pp.16-28
• /
• 2000

The steady rotating flows in eccentric annulus has been studied by a numerical method based on the spectral collocation method. The inner cylinder has a constant angular velocity while the outer on e is stationary. Flow between eccentric cylinders is of considerable technical importance as it occurs in journal bearings. In the present work, the governing equations for laminar flow are expressed as Navier-Stokes equations, including the non-linear convection terms. The solutions were utilized i, estimate the effects of the nonlinear terms on the load acting on the rotating cylinder. Based on the half and the full Sommerfeld methods, the load on the rotating cylinder is evaluated with eccentricity, by integrating the pressure and skin friction around the cylinder. The attitude angle and Sommerfeld reciprocal are calculated from the load. Also, the torque on the rotating inner cylinder was calculated. considering the skin friction. The attitude angle and Sommerfeld reciprocal are decreased with eccentricity. Viscous damping coefficient due to the skin friction becomes larger with decreasing the annular space. It is found the non-linear effects of the convection terms on the flow and the load are important. especially on the attitude angle, for relatively wide annular configurations however, the effects on those are minor for very narrow annular ones.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.25 no.4
• /
• pp.173-195
• /
• 2021

A new implicit discontinuous Galerkin spectral element method (DGSEM) based on the first order hyperbolic system(FOHS) is presented for solving elliptic type partial different equations, such as the Poisson problems. By utilizing the idea of hyperbolic formulation of Nishikawa[1], the original Poisson equation was reformulated in the first-order hyperbolic system. Such hyperbolic system is solved implicitly by the collocation type DGSEM. The steady state solution in pseudo-time, which is the solution of the original Poisson problem, was obtained by the implicit solution of the global linear system. The optimal polynomial orders of 𝒪(𝒽^𝑝+1)) are obtained for both the solution and gradient variables from the test cases in 1D and 2D regular grids. Spectral accuracy of the solution and gradient variables are confirmed from all test cases of using the uniform grids in 2D.

• Journal of KSNVE
• /
• v.11 no.3
• /
• pp.428-436
• /
• 2001

The rotating inner cylinder executes a periodic translational motion in concentric annulus while the outer one is stationary. In the study of flow-induced vibrations and relaxed instabilities, it is of interest to evaluate the fluid-dynamic forces acting on the rotating inner cylinder. In the present work, the governing equations for the confined flow are expressed as Navier-Stokes equations, including the steady and unsteady terms. The fluid parameters for steady flow generated by the rotating cylinder are determined analytically while the unsteady ones by the oscillatory motion are evaluates by a numerical method based on the spectral collocation method. In order to validate the numerical approach, the numerical results are compared wish the analytical ones given by existing theories, for simple cases where the both approaches are applicable. Good agreement was found between the results. It is found the effects of the Reynolds number, defined by rotating velocity, on the fluid-dynamic forces are important for the case of relatively low oscillatory Reynolds number, defined by oscillatory frequency : j.e., in case of $Re_\omega\gg Re_S$.

• Honam Mathematical Journal
• /
• v.37 no.3
• /
• pp.299-315
• /
• 2015

In this paper, we study the first-order system least-squares (FOSLS) spectral method for parabolic partial differential equations. There were lots of least-squares approaches to solve elliptic partial differential equations using finite element approximation. Also, some approaches using spectral methods have been studied in recent. In order to solve the parabolic partial differential equations in parallel, we consider a parallel numerical method based on a hybrid method of the frequency-domain method and first-order system least-squares method. First, we transform the parabolic problem in the space-time domain to the elliptic problems in the space-frequency domain. Second, we solve each elliptic problem in parallel for some frequencies using the first-order system least-squares method. And then we take the discrete inverse Fourier transforms in order to obtain the approximate solution in the space-time domain. We will introduce such a hybrid method and then present a numerical experiment.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.22 no.7
• /
• pp.1030-1040
• /
• 1998

One of the main unresolved issues in large-eddy simulation(LES) of wall-bounded turbulent flows is the requirement of high spatial resolution in the near-wall region, especially in the spanwise direction. Such high resolution required in the near-wall region is generally used throughout the computational domain, making simulations of high Reynolds number, complex-geometry flows prohibitive. A grid-embedding strategy using a nonconforming spectral domain-decomposition method is proposed to address this limitation. This method provides an efficient way of clustering grid points in the near-wall region with spectral accuracy. LES of transitional and turbulent channel flow has been performed to evaluate the proposed grid-embedding technique. The computational domain is divided into three subdomains to resolve the near-wall regions in the spanwise direction. Spectral patching collocation methods are used for the grid-embedding and appropriate conditions are suggested for the interface matching. Results of LES using the grid-embedding strategy are promising compared to LES of global spectral method and direct numerical simulation. Overall, the results show that the spectral domain-decomposition grid-embedding technique provides an efficient method for resolving the near-wall region in LES of complex flows of engineering interest, allowing significant savings in the computational CPU and memory.

• International Journal of Aeronautical and Space Sciences
• /
• v.16 no.2
• /
• pp.264-277
• /
• 2015

This article investigates various transcription techniques for the Legendre pseudospectral (PS) method to compare the pros and cons of each approach. Eight combinations from four different types of collocation points and two discretization methods for dynamic constraints, which differentiate Legendre PS transcription techniques, are implemented to solve a carefully selected test set of nonlinear optimal control problems (NOCPs). The convergence property and prediction accuracy are compared to provide a useful guideline for selecting the best combination. The tested NOCPs consist of the minimum time, minimum energy, and problems with state and control constraints. Therefore, the results drawn from this comparative study apply to the solution of similar types of NOCPs and can mitigate much debate about the best combinations. Additionally, important findings from this study can be used to improve the numerical efficiency of the Legendre PS methods. Three PS applications to the aerospace engineering problems are demonstrated to prove this point.

• International Journal of Aeronautical and Space Sciences
• /
• v.15 no.2
• /
• pp.183-189
• /
• 2014

Pseudo-spectral method typically converges at an exponential rate. However, it requires a special set of fixed collocation points (CPs) to get highly accurate formulas for partial integration and differentiation. In this study, computational nodes for defining the discrete variables of states and controls are built independently of the CPs. The state and control variables at each CP, which are required to transcribe an NOCP into the corresponding NLP, are interpolated, using those variables allocated at each node. Additionally, Lagrange interpolation and spline interpolation are investigated, to provide a guideline for selecting a favorable interpolation method. The proposed techniques are applied to the solution of an NOCP using equally spaced nodes, and the computed results are compared to those using the standard PS approach, to validate the usefulness of the proposed methods.