• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.23 no.1
• /
• pp.1-18
• /
• 2019

In contrast to the well-developed convex splitting schemes for gradient flows of two-component system, there were few efforts on applying the convex splitting idea to gradient flows of multi-component system, such as the vector-valued Cahn-Hilliard (vCH) equation. In the case of the vCH equation, one need to consider not only the convex splitting idea but also a specific method to manage the partition of unity constraint to design an unconditionally energy stable scheme. In this paper, we propose a constrained Convex Splitting (cCS) scheme for the vCH equation, which is based on a convex splitting of the energy functional for the vCH equation under the constraint. We show analytically that the cCS scheme is mass conserving and unconditionally uniquely solvable. And it satisfies the constraint at the next time level for any time step thus is unconditionally energy stable. Numerical experiments are presented demonstrating the accuracy, energy stability, and efficiency of the proposed cCS scheme.

• Bulletin of the Korean Mathematical Society
• /
• v.56 no.1
• /
• pp.265-276
• /
• 2019

For a simple implementation, a linear convex splitting scheme was coupled with the Fourier spectral method for the Cahn-Hilliard equation with a logarithmic free energy. However, an inappropriate value of the splitting parameter of the linear scheme may lead to incorrect morphologies in the phase separation process. In order to overcome this problem, we present a nonlinear convex splitting Fourier spectral scheme for the Cahn-Hilliard equation with a logarithmic free energy, which is an appropriate extension of Eyre's idea of convex-concave decomposition of the energy functional. Using the nonlinear scheme, we derive a useful formula for the relation between the gradient energy coefficient and the thickness of the interfacial layer. And we present numerical simulations showing the different evolution of the solution using the linear and nonlinear schemes. The numerical results demonstrate that the nonlinear scheme is more accurate than the linear one.

• Bulletin of the Korean Mathematical Society
• /
• v.61 no.4
• /
• pp.969-998
• /
• 2024

We present numerical schemes to deal with nonconservative terms in the two-dimensional shallow water equations with variable topography. Relying on the dimensional splitting technique, we construct Godunov-type schemes. Such schemes can be categorized into two classes, namely the partly and fully splitting ones, depending on how deeply the scheme employs the splitting method. An upwind scheme technique is employed for the evolution of the velocity component for the partly splitting scheme. These schemes are shown to possess interesting properties: They can preserve the positivity of the water height, and they are well-balanced.

• KSII Transactions on Internet and Information Systems (TIIS)
• /
• v.14 no.1
• /
• pp.1-19
• /
• 2020

In an underlay cognitive simultaneous wireless information and power transfer (SWIPT) network, communication from secondary user (SU) to secondary destination (SD) is accomplished with decode-and-forward (DF) relays. Multiple energy-constrained relays are assumed to harvest energy from SU via power splitting (PS) protocol and complete SU secure information transmission with beamforming. Hence, physical layer security (PLS) is investigated in cognitive SWIPT network. In order to interfere with eavesdropper and improve relay's energy efficiency, a destination-assisted jamming scheme is proposed. Namely, SD transmits artificial noise (AN) to interfere with eavesdropping, while jamming signal can also provide harvested energy to relays. Beamforming vector and power splitting ratio are jointly optimized with the objective of SU secrecy capacity maximization. We solve this non-convex optimization problem via a general two-stage procedure. Firstly, we obtain the optimal beamforming vector through semi-definite relaxation (SDR) method with a fixed power splitting ratio. Secondly, the best power splitting ratio can be obtained by one-dimensional search. We provide simulation results to verify the proposed solution. Simulation results show that the scheme achieves the maximum SD secrecy rate with appropriate selection of power splitting ratio, and the proposed scheme guarantees security in cognitive SWIPT networks.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.17 no.3
• /
• pp.197-207
• /
• 2013

The Cahn-Hilliard equation was proposed as a phenomenological model for describing the process of phase separation of a binary alloy. The equation has been applied to many physical applications such as amorphological instability caused by elastic non-equilibrium, image inpainting, two- and three-phase fluid flow, phase separation, flow visualization and the formation of the quantum dots. To solve the Cahn-Hillard equation, many numerical methods have been proposed such as the explicit Euler's, the implicit Euler's, the Crank-Nicolson, the semi-implicit Euler's, the linearly stabilized splitting and the non-linearly stabilized splitting schemes. In this paper, we investigate each scheme in finite-difference schemes by comparing their performances, especially stability and efficiency. Except the explicit Euler's method, we use the fast solver which is called a multigrid method. Our numerical investigation shows that the linearly stabilized stabilized splitting scheme is not unconditionally gradient stable in time unlike the known result. And the Crank-Nicolson scheme is accurate but unstable in time, whereas the non-linearly stabilized splitting scheme has advantage over other schemes on the time step restriction.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.32 no.10C
• /
• pp.942-949
• /
• 2007

In this paper, we propose a new rate-control scheme, called splining, to construct low-rate codes from high-rate codes by splitting rows of the parity-check matrices of LDPC codes, which can construct rate-compatible LDPC codes having good initial transmission performance. Good low-rate codes can be constructed by making the number of distinct check node degrees as small as possible after splitting. The proposed scheme achieves good cycle property, low decoding complexity, and fast convergence speed, especially compared to the puncturing. Especially, rate-compatible repeat accumulate-type LDPC (RA-Type LDPC) code is constructed using splitting, which covers the code rates from 1/3 to 4/5. Through simulation it is shown that this code outperforms other rate-compatible RA-Type LDPC codes for all rates and can be decoded conveniently and efficiently.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.12 no.2
• /
• pp.70-80
• /
• 2000

Princeton Ocean Model (POM) is modified to construct a three-dimensional, semi-implicit hydro¬dynamic model with a wetting-and-drying scheme. The model employs semi-implicit treatment of the barotropic pressure gradient terms and the vertical mixing terms in the momentum equations, and the velocity divergence term in the vertically-integrated continuity equation. Such treatment removes the external mode and thus the mode splitting scheme in POM, allowing the semi-implicit model to use a larger time step. Applied to hypothetical systems, both the semi-implicit model and POM give nearly the same results. The semi-implicit model, however, runs approximately 4.4 times faster than POM showing its improved computational efficiency. Applied to a hypothetical system with intertidal flats, POM employing the mode splitting scheme produces noises at the intertidal flats, that propagate into the main channel resulting in unstable current velocities. Despite its larger time step, the semi-implicit model gives stable current velocities both at the intertidal flats and main channel. The semi-implicit model when applied to Kyeonggi Bay gives a good reproduction of the observed tides and tidal currents throughout the modeling domain, demonstrating its prototype applicability.

• Journal of the Society of Naval Architects of Korea
• /
• v.47 no.1
• /
• pp.11-19
• /
• 2010

A code is developed to simulate incompressible free surface flows using the Roe's flux-difference splitting scheme. An interface of two fluids is considered as a moving contact discontinuity. The continuities of pressure and normal velocity across the interface are enforced by the conservation law in the integral sense. The fluxes are computed using the Roe's flux-difference splitting scheme for two incompressible fluids. The interface can be identified based on the computed density distribution. However, no additional treatment is required along the interface during the whole computations. Complicated time evolution of the interface including topological change can be captured without any difficulties. The developed code is applied to simulate the Rayleigh-Taylor instability of two incompressible fluids in the density ratio of 7.2:1 and the broken dam problem of water-air. The present results are compared with other available results and good agreements are achieved for the both cases.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.189-190
• /
• 2003

A New kinetic scheme for the compressible Navier-Stokes equations is developed. While the conventional approach, such as KFVS scheme, employs the splitting algorithm and computes the numerical flux on the basis of the collisionless equation, the present approach employs the splitting algorithm in the evaluation of the numerical flux, where the collision effect is explicitly taken into account. However, the initial condition employed in the computation is slightly different from the conventional Chapman-Enskog NS distribution function. The present study also reveals the background of the existing kinetic schemes. such as the KFVS scheme and Gas-Kinetic BGK scheme.

• Journal of computational fluids engineering
• /
• v.8 no.4
• /
• pp.1-5
• /
• 2003

A Thin-layer Wavier-Stokes equations are applied for the hypersonic flow over blunt bodies with applications to laminar as well as turbulent flows. The equations are expressed in the forms of flux-vector splitting and explicit algorithm. The upwind schemes of Steger-Warming and Van Leer are investigated to predict accurately the heating loads along the surface of the body. A mixed scheme has been presented for the differencing the convective terms and the mixed scheme is found to be less dissipative producing accurate solutions.