• Journal of computational fluids engineering
• /
• v.19 no.3
• /
• pp.91-97
• /
• 2014

A high resolution numerical method aimed at solving cavitating flow was proposed and applied to gas-liquid two-phase shock tube problem with arbitrary void fraction. The present method with compressibility effects employs a finite-difference 4th-order Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL TVD scheme. The Jacobian matrix from the inviscid flux of constitute equation is diagonalized analytically and the speed of sound for the two-phase media is derived by eigenvalues. So that the present method is appropriate for the extension of high order upwind schemes based on the characteristic theory. By this method, a Riemann problem for Euler equations of one dimensional shock tube was computed. Numerical results of high speed flow phenomena such as detailed observations of shock and expansion wave propagations through the gas-liquid two-phase media and some data related to computational efficiency are made. Comparisons of predicted results and solutions at isothermal condition are provided and discussed.

• The Journal of the Acoustical Society of Korea
• /
• v.33 no.5
• /
• pp.282-290
• /
• 2014

Though there have been advances in the numerical analysis of the acoustic scattering by smooth objects, numerical analysis of the acoustic scattering by the objects that are partially exposed on the interface are still rare. In determining the backscattering amplitude by a partially buried cylinder on a seabed, reverberation by the interface changes the feature of the scattering form function. Current study adopted the Method of moments (MoM) to provide the numerical analysis on the backscattering amplitude for a partially buried cylinder on a flat interface. Suggested numerical analysis showed the good agreements with the measurements and the analytic solution obtained by the Kirchhoff approximation. Numerical analysis described in the current study can be applied to the backscattering problem of any shape of the objects partially imbedded on a seabed by combining the reverberation from the seabed with the scattered wave from the objects.

• KSII Transactions on Internet and Information Systems (TIIS)
• /
• v.12 no.4
• /
• pp.1504-1526
• /
• 2018

For uplink multi-user massive MIMO systems, conventional minimum mean square error (MMSE) linear detection method achieves near-optimal performance when the number of antennas at base station is much larger than that of the single-antenna users. However, MMSE detection involves complicated matrix inversion, thus making it cumbersome to be implemented cost-effectively and rapidly. In this paper, we first summarize in detail the state-of-the-art simplified MMSE detection algorithms that circumvent the complicated matrix inversion and hence reduce the computation complexity from ${\mathcal{O}}(K^3)$ to ${\mathcal{O}}(K^2)$ or ${\mathcal{O}}(NK)$ with some certain performance sacrifice. Meanwhile, we divide the simplified algorithms into two categories, namely the matrix inversion approximation and the classical iterative linear equation solving methods, and make comparisons between them in terms of detection performance and computation complexity. In order to further optimize the detection performance of the existing detection algorithms, we propose more proper solutions to set the initial values and relaxation parameters, and present a new way of reconstructing the exact effective noise variance to accelerate the convergence speed. Analysis and simulation results verify that with the help of proper initial values and parameters, the simplified matrix inversion based detection algorithms can achieve detection performance quite close to that of the ideal matrix inversion based MMSE algorithm with only a small number of series expansions or iterations.

• KSII Transactions on Internet and Information Systems (TIIS)
• /
• v.12 no.6
• /
• pp.2595-2618
• /
• 2018

In this paper, we investigate secure communication with the presence of multiple eavesdroppers (Eves) in a two-tier downlink dense heterogeneous network, wherein there is a macrocell base station (MBS) and multiple femtocell base stations (FBSs). Each base station (BS) has multiple users. And Eves attempt to wiretap a macrocell user (MU). To keep Eves ignorant of the confidential message, we propose a physical-layer security scheme based on cross-layer cooperation to exploit interference in the considered network. Under the constraints on the quality of service (QoS) of other legitimate users and transmit power, the secrecy rate of system can be maximized through jointly optimizing the beamforming vectors of MBS and cooperative FBSs. We explore the problem of maximizing secrecy rate in both non-colluding and colluding Eves scenarios, respectively. Firstly, in non-colluding Eves scenario, we approximate the original non-convex problem into a few semi-definite programs (SDPs) by employing the semi-definite relaxation (SDR) technique and conservative convex approximation under perfect channel state information (CSI) case. Furthermore, we extend the frame to imperfect CSI case and use the Lagrangian dual theory to cope with uncertain constraints on CSI. Secondly, in colluding Eves scenario, we transform the original problem into a two-tier optimization problem equivalently. Among them, the outer layer problem is a single variable optimization problem and can be solved by one-dimensional linear search. While the inner-layer optimization problem is transformed into a convex SDP problem with SDR technique and Charnes-Cooper transformation. In the perfect CSI case of both non-colluding and colluding Eves scenarios, we prove that the relaxation of SDR is tight and analyze the complexity of proposed algorithms. Finally, simulation results validate the effectiveness and robustness of proposed scheme.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.22 no.10
• /
• pp.1012-1015
• /
• 2011

The decision-directed blind equalization algorithm is often used due to its simplicity and good convergence property when the eye pattern is open. However, in a channel where the eye pattern is closed, the decision-directed algorithm is not guaranteed to converge. Hence, a modified Bussgang-type algorithm using a hyperbolic tangent function for zero-memory nonlinear(ZNL) function has been proposed and applied to avoid this problem by Filho et al. But application of this algorithm includes the calculation of hyperbolic tangent function and its derivative or a look-up table which may need a large amount of memory due to channel variations. To reduce the computational and/or hardware complexity of Filho's algorithm, in this paper, an improved method for the decision-directed algorithm is proposed. In the proposed scheme, the ZNL function and its derivative are respectively set to be the original signum function and a narrow rectangular pulse which is an approximation of Dirac delta function. It is shown that the proposed scheme, when it is combined with decision-directed algorithm, reduces the computational complexity drastically while it retains the convergence and steady-state performance of the Filho's algorithm.

• Journal of the Korea Computer Graphics Society
• /
• v.5 no.1
• /
• pp.27-33
• /
• 1999

In this paper, we propose an efficient technique for the animation of flexible thin objects. Mass-spring model was employed to represent the flexible objects. Till now, many techniques have used the mass-spring model to generate plausible animation of soft objects. A straight-forward approach to the animation with mass-spring model is explicit Euler method, but the explicit Euler method has serious disadvantage that it suffers from 'instability problem'. The implicit integration method is a possible solution to overcome the instability problem. However, the most critical flaw of the implicit method is that it involves a large linear system. This paper presents a fast animation technique for mass-spring model with approximated implicit method. The proposed technique stably updates the state of n mass-points in O(n) time when the number of total springs are O(n). We also consider the interaction of the flexible object and air in order to generate plausible results.

• Structural Engineering and Mechanics
• /
• v.56 no.4
• /
• pp.549-567
• /
• 2015

A first-order moment method (FORM) reliability analysis is commonly used for structural stability analysis. It requires the values and partial derivatives of the performance to function with respect to the random variables for the design. These calculations can be cumbersome when the performance functions are implicit. A Gaussian process (GP)-based response surface is adopted in this study to approximate the limit state function. By using a trained GP model, a large number of values and partial derivatives of the performance functions can be obtained for conventional reliability analysis with a FORM, thereby reducing the number of stability analysis calculations. This dynamic renewed knowledge source can provide great assistance in improving the predictive capacity of GP during the iterative process, particularly from the view of machine learning. An iterative algorithm is therefore proposed to improve the precision of GP approximation around the design point by constantly adding new design points to the initial training set. Examples are provided to illustrate the GP-based response surface for both structural and non-structural reliability analyses. The results show that the proposed approach is applicable to structural reliability analyses that involve implicit performance functions and structural response evaluations that entail time-consuming finite element analyses.

• Proceedings of the KSRS Conference
• /
• v.1
• /
• pp.86-89
• /
• 2006

Ocean color remote sensing community currently uses the different solar irradiance spectra covering the visible and near-infrared in the calibration/validation and deriving products of ocean color instruments. These spectra derived from single and / or multiple measurements sets or models have significant discrepancies, primarily due to variation of the solar activity and uncertainties in the measurements from various instruments and their different calibration standards. Thus, it is prudent to examine model-to-model differences and select a standard reference spectrum that can be adopted in the future calibration and validation processes, particularly of the first Geostationary Ocean Color Imager (GOCI) onboard its Communication Ocean and Meterological Satellite (COMS) planned to be launched in 2008. From an exhaustive survey that reveals a variety of solar spectra in the literature, only eight spectra are considered here seeing as reference in many remote sensing applications. Several criteria are designed to define the reference spectrum: i.e., minimum spectral range of 350-1200nm, based completely or mostly on direct measurements, possible update of data and less errors. A careful analysis of these spectra reveals that the Thuillier 2004 spectrum seems to be very identical compared to other spectra, primarily because it represents very high spectral resolution and the current state of the art in solar irradiance spectra of exceptionally low uncertainty ${\sim}0.1%.$ This study also suggests use of the Gueymard 2004 spectrum as an alternative for applications of multispectral/multipurpose satellite sensors covering the terrestrial regions of interest, where it provides spectral converge beyond 2400nm of the Thuillier 2004 spectrum. Since the solar-activity induced spectral variation is about less than 0.1% and a large portion of this variability occurs particularly in the ultraviolet portion of the electromagnetic spectrum that is the region of less interest for the ocean color community, we disregard considering this variability in the analysis of solar irradiance spectra, although determine the solar constant 1366.1 $Wm^{-2}$ to be proposed for an improved approximation of the extraterrestrial solar spectrum in the visible and NIR region.

• Nuclear Engineering and Technology
• /
• v.24 no.1
• /
• pp.63-74
• /
• 1992

The resonant formation rates of muonic molecular ion dtr in the muon-catalyzed fusion are calculated in various fuel temperatures and densities. The elastic scattering cross sections between t$\mu$ and deuterons are obtained by making use of the partial wave method. The transition property of the excited compound molecule [(dt$\mu$)dee]＊ derived by the impulse approximation in the form of a bound-state form factor. The radiative, Auger, and collisional deexcitations are considered as the deexcitation mechanisms of the excited dt$\mu$, and each deexcitation width is calculated as well as back decay width. The resultant reaction widths are used to calculate the formation cross sections of resonant dt$\mu$. The resonant formation rates for dt$\mu$-d and dt$\mu$-t collisions are computed as functions of fuel temperature and density. The calculations show that the resonant formation rates increase with fuel densities and have the maximum values at the particular temperatures where the relative collision energies are equal to the resonant ones.

• Steel and Composite Structures
• /
• v.24 no.5
• /
• pp.549-559
• /
• 2017

As few multi-tower single-box multi-cell cable-stayed bridges with corrugated steel webs have been built, analysis is mostly achieved by combining single-girder model, beam grillage model and solid model in support of the design. However, such analysis methods usually suffer from major limitations in terms of the engineering applications: single-girder model fails to account for spatial effect such as shear lag effect of the box girder and the relevant effective girder width and eccentric load coefficient; owing to the approximation in the principle equivalence, the plane grillage model cannot accurately capture shear stress distribution and local stress state in both top and bottom flange of composite box girder; and solid model is difficult to be practically combined with the overall calculation. The usual effective width method fails to provide a uniform and accurate "effective length" (and the codes fail to provide a unified design approach at those circumstance) considering different shear lag effects resulting from dead load, prestress and cable tension in the construction. Therefore, a novel spatial grid model has been developed to account for shear lag effect. The theoretical principle of the proposed spatial grid model has been elaborated along with the relevant illustrations of modeling parameters of composite box girder with corrugated steel webs. Then typical transverse and longitudinal shear lag coefficient distribution pattern at the side-span and mid-span key cross sections have been analyzed and summarized to provide reference for similar bridges. The effectiveness and accuracy of spatial grid analysis methods has been finally validated through a practical cable-stayed bridge.