• Journal of applied mathematics & informatics
• /
• v.34 no.1_2
• /
• pp.19-34
• /
• 2016

In this paper, we present a split least-squares characteristic mixed finite element method(MFEM) to get the approximate solutions of the convection dominated Sobolev equations. First, to manage both convection term and time derivative term efficiently, we apply a least-squares characteristic MFEM to get the system of equations in the primal unknown and the flux unknown. Then, we obtain a split least-squares characteristic MFEM to convert the coupled system in two unknowns derived from the least-squares characteristic MFEM into two uncoupled systems in the unknowns. We theoretically prove that the approximations constructed by the split least-squares characteristic MFEM converge with the optimal order in L² and H¹ normed spaces for the primal unknown and with the optimal order in L² normed space for the flux unknown. And we provide some numerical results to confirm the validity of our theoretical results.

• East Asian mathematical journal
• /
• v.15 no.1
• /
• pp.121-133
• /
• 1999

• Journal of the Korean Mathematical Society
• /
• v.46 no.5
• /
• pp.1007-1025
• /
• 2009

First-order least-squares method of a distributed optimal control problem for the incompressible Stokes equations is considered. An optimality system for the optimal solution are reformulated to the equivalent first-order system by introducing the vorticity and then the least-squares functional corresponding to the system is defined in terms of the sum of the squared $H^{-1}$ and $L^2$ norms of the residual equations of the system. Finite element approximations are studied and optimal error estimates are obtained. Resulting linear system of the optimality system is symmetric and positive definite. The V-cycle multigrid method is applied to the system to test computational efficiency.

• 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.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.10
• /
• pp.26-34
• /
• 2005

In this paper, a novel moving least squares interface welding method is proposed to carry out the coupled analysis of whole model composed of independently modeled finite element substructures with nodal mismatching interfaces. To verify the validity, and efficiency of the proposed interface welding method, various numerical examples are worked out including patch tests, convergence tests, and examples of coupled analyses of the structural systems with mismatching substructures. From the numerical tests, it is confirmed that one can efficiently carry out the coupled analysis of whole model composed of mismatching finite element substructures through the proposed method without any remeshing or any additional unknown.

• Communications of the Korean Mathematical Society
• /
• v.20 no.3
• /
• pp.563-578
• /
• 2005

In [Z. Cai and B. C. Shin, SIAM J. Numer. Anal. 40 (2002), 307-318], we developed the discrete first-order system least squares method for the second-order elliptic boundary value problem by directly approximating $H(div){\cap}H(curl)-type$ space based on the Helmholtz decomposition. Under general assumptions, error estimates were established in the $L^2\;and\;H^1$ norms for the vector and scalar variables, respectively. Such error estimates are optimal with respect to the required regularity of the solution. In this paper, we study solution methods for solving the system of linear equations arising from the discretization of variational formulation which possesses discrete biharmonic term and focus on numerical results including the performances of multigrid preconditioners and the finite element accuracy.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.31 no.4
• /
• pp.26-34
• /
• 2003

In this paper, a novel technique is proposed to arbitrarily activate the nodal points in finite element model through the meshless approximation methods such as MLS(moving least squares method), and theoretical investigations are carried out including the consistency and boundeness of numerical solution to prove the validity of the proposed method. By using the proposed node activation technique, one can activate and handle only the concerned nodes as unknown variables among the large number of nodal points in the finite element model. Therefore, the proposed technique has a great potential in design and reanalysis procedure.

• Proceedings of the KSME Conference
• /
• 2007.05a
• /
• pp.354-359
• /
• 2007

A new level set based topology optimization employing inner-front creation algorithm is presented. In the conventional level set based topology optimization, the optimum topology strongly depends on the initial level set distribution due to the incapability of inner-front creation during optimization process. In the present work, an inner-front creation algorithm is proposed, in which the sizes, positions, and number of new inner-fronts during the optimization process can be globally and consistently identified. To update the level set function during the optimization process, the least-squares finite element method is employed. As demonstrative examples for the flexibility and usefulness of the proposed method, the level set based topology optimization considering lightweight design of 3D shell structure is carried out.

• Journal of Ocean Engineering and Technology
• /
• v.24 no.1
• /
• pp.20-33
• /
• 2010

The paper deals with strength design of a riser support installed on floating production storage and offloading (FPSO) vessel under various loading conditions - operation, extreme, damaged, one line failure case (OLFC) and installation. The design problem is formulated such that thickness sizing variables are determined by minimizing the weight of a riser support structure subject to stresses constraints. The initial design model is generated based on an actual FPSO riser support specification. The finite element analysis (FEA) is conducted using MSC/NASTRAN, and optimal solutions are obtained via moving least squares method (MLSM) in the context of response surface based approximate optimization. For the meta-modeling of inequality constraint functions of stresses, a constraint-feasible moving least squares method (CF-MLSM) is used in the present study. The method of CF-MLSM, compared to a conventional MLSM, has been shown to ensure the constraint feasibility in a case where the approximate optimization process is employed. The optimization results present improved design performances under various riser operating conditions.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.14 no.5
• /
• pp.957-964
• /
• 2011

New time schemes based on the FEM were developed and their performances were tested with 2D wave equation. The least-squares and weighted residual methods are used to construct new time schemes based on traditional residual minimization method. To overcome some drawbacks that time schemes based on the least-squares and weighted residual methods have, ad-hoc method is considered to minimize residuals multiplied by others residuals as a new approach. And variational method is used to get necessary conditions of ad-hoc minimization. A-stability was chosen to check the stability of newly developed time schemes. Specific values of new time schemes are presented along with their numerical solutions which were compared with analytic solution.