• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.4 no.4
• /
• pp.219-224
• /
• 1992

A numerical method using the truncated Fourier series is presented to predict the wave potential and water surface profile for two dimensional nonlinear standing waves. The unknown coefficients of the series are to be determined through the Newton solution of nonlinear simultaneous equations given by the governing equation and boundary conditions of the problem. In order to prove the effectiveness of the present method. an existing Stokes-like perturbation method is considered together, and a hydraulic experiment for measuring water surface profile and wave pressure is performed as well. The results are such that the present method can generally give exact solutions even for relatively big wave stiffness regardless of the water depth condition. It also demonstrates its validity by showing double humps in the crest of temporal wave pressure profile which normally appear in strongly nonlinear standing waves.

• Structural Engineering and Mechanics
• /
• v.17 no.1
• /
• pp.1-14
• /
• 2004

Numerical solution to linear bending analysis of circular plates is obtained by the method of harmonic differential quadrature (HDQ). In the method of differential quadrature (DQ), partial space derivatives of a function appearing in a differential equation are approximated by means of a polynomial expressed as the weighted linear sum of the function values at a preselected grid of discrete points. The method of HDQ that was used in the paper proposes a very simple algebraic formula to determine the weighting coefficients required by differential quadrature approximation without restricting the choice of mesh grids. Applying this concept to the governing differential equation of circular plate gives a set of linear simultaneous equations. Bending moments, stresses values in radial and tangential directions and vertical deflections are found for two different types of load. In the present study, the axisymmetric bending behavior is considered. Both the clamped and the simply supported edges are considered as boundary conditions. The obtained results are compared with existing solutions available from analytical and other numerical results such as finite elements and finite differences methods. A comparison between the HDQ results and the finite difference solutions for one example plate problem is also made. The method presented gives accurate results and is computationally efficient.

• Journal of Information Processing Systems
• /
• v.13 no.5
• /
• pp.1183-1202
• /
• 2017

Image restoration has been carried out by texture synthesis mostly for large regions and inpainting algorithms for small cracks in images. In this paper, we propose a new approach that allows for the simultaneous fill-in of different structures and textures by processing in a wavelet domain. A combination of structure inpainting and patch-based texture synthesis is carried out, which is known as patch-based inpainting, for filling and updating the target region. The wavelet transform is used for its very good multiresolution capabilities. The proposed algorithm uses the wavelet domain subbands to resolve the structure and texture components in smooth approximation and high frequency structural details. The subbands are processed separately by the prioritized patch-based inpainting with isophote energy driven texture synthesis at the core. The algorithm automatically estimates the wavelet coefficients of the target regions of various subbands using optimized patches from the surrounding DWT coefficients. The suggested performance improvement drastically improves execution speed over the existing algorithm. The proposed patch optimization strategy improves the quality of the fill. The fill-in is done with higher priority to structures and isophotes arriving at target boundaries. The effectiveness of the algorithm is demonstrated with natural and textured images with varying textural complexions.

• Nuclear Engineering and Technology
• /
• v.53 no.2
• /
• pp.386-392
• /
• 2021

MACCS (MELCOR Accident Consequence Code System), WinMACCS, and MelMACCS now facilitate a multi-unit consequence analysis. MACCS evaluates the consequences of an atmospheric release of radioactive gases and aerosols into the atmosphere and is most commonly used to perform probabilistic safety assessments (PSAs) and related consequence analyses for nuclear power plants (NPPs). WinMACCS is a user-friendly preprocessor for MACCS. MelMACCS extracts source-term information from a MELCOR plot file. The current development can combine an arbitrary number of source terms, representing simultaneous releases from a multi-unit facility, into a single consequence analysis. The development supports different release signatures, fission product inventories, and accident initiation times for each unit. The treatment is completely general except that the model is currently limited to collocated units. A major practical consideration for performing a multi-unit PSA is that a comprehensive treatment for more than two units may involve an intractable number of combinations of source terms. This paper proposes and evaluates an approach for reducing the number of calculations to be tractable, even for sites with eight or ten units. The approximation error introduced by the approach is acceptable and is considerably less than other errors and uncertainties inherent in a Level 3 PSA.

• Geophysics and Geophysical Exploration
• /
• v.11 no.1
• /
• pp.68-77
• /
• 2008

We have developed an inversion algorithm for loop-loop electromagnetic (EM) data, based on the localised non-linear or extended Born approximation to the solution of the 2.5D integral equation describing an EM scattering problem. Source and receiver configuration may be horizontal co-planar (HCP) or vertical co-planar (VCP). Both multi-frequency and multi-separation data can be incorporated. Our inversion code runs on a PC platform without heavy computational load. For the sake of stable and high-resolution performance of the inversion, we implemented an algorithm determining an optimum spatially varying Lagrangian multiplier as a function of sensitivity distribution, through parameter resolution matrix and Backus-Gilbert spread function analysis. Considering that the different source-receiver orientation characteristics cause inconsistent sensitivities to the resistivity structure in simultaneous inversion of HCP and VCP data, which affects the stability and resolution of the inversion result, we adapted a weighting scheme based on the variances of misfits between the measured and calculated datasets. The accuracy of the modelling code that we have developed has been proven over the frequency, conductivity, and geometric ranges typically used in a loop-loop EM system through comparison with 2.5D finite-element modelling results. We first applied the inversion to synthetic data, from a model with resistive as well as conductive inhomogeneities embedded in a homogeneous half-space, to validate its performance. Applying the inversion to field data and comparing the result with that of dc resistivity data, we conclude that the newly developed algorithm provides a reasonable image of the subsurface.

• Tunnel and Underground Space
• /
• v.31 no.6
• /
• pp.480-493
• /
• 2021

Since the generalized Hoek-Brown criterion (GHB) provides an efficient way of identifying its strength parameter values with the consideration of in-situ rock mass condition via Geological Strength Index (GSI), this criterion is recognized as one of the standard rock mass failure criteria in rock mechanics community. However, the nonlinear form of the GHB criterion makes its mathematical treatment inconvenient and limits the scope of its application. As an effort to overcome this disadvantage of the GHB criterion, the explicit approximate analytical equations for the minimum principal stress, which is associated with the maximum principal stress at failure, are formulated based on the Taylor polynomial approximation of the original GHB criterion. The accuracy of the derived approximate formula for the minimum principal stress is verified by comparing the resulting approximate minimum principal stress with the numerically calculated exact values. To provide an application example of the approximate formulation, the equivalent friction angle and cohesion for the expected plastic zone around a circular tunnel in a GHB rock mass are calculated by incorporating the formula for the approximate minimum principal stress. It is found that the simultaneous consideration of the values of mi, GSI and far-field stress is important for the accurate calculation of equivalent Mohr-Coulomb parameter values of the plastic zone.

• Journal of Computational Design and Engineering
• /
• v.1 no.3
• /
• pp.187-193
• /
• 2014

One of the hottest challenges in automotive industry is related to weight reduction in sheet metal forming processes, in order to produce a high quality metal part with minimal material cost. Stamping is the most widely used sheet metal forming process; but its implementation comes with several fabrication flaws such as springback and failure. A global and simple approach to circumvent these unwanted process drawbacks consists in optimizing the initial blank shape with innovative methods. The aim of this paper is to introduce an efficient methodology to deal with complex, computationally expensive multicriteria optimization problems. Our approach is based on the combination of methods to capture the Pareto Front, approximate criteria (to save computational costs) and global optimizers. To illustrate the efficiency, we consider the stamping of an industrial workpiece as test-case. Our approach is applied to the springback and failure criteria. To optimize these two criteria, a global optimization algorithm was chosen. It is the Simulated Annealing algorithm hybridized with the Simultaneous Perturbation Stochastic Approximation in order to gain in time and in precision. The multicriteria problems amounts to the capture of the Pareto Front associated to the two criteria. Normal Boundary Intersection and Normalized Normal Constraint Method are considered for generating a set of Pareto-optimal solutions with the characteristic of uniform distribution of front points. The computational results are compared to those obtained with the well-known Non-dominated Sorting Genetic Algorithm II. The results show that our proposed approach is efficient to deal with the multicriteria shape optimization of highly non-linear mechanical systems.

• Ocean and Polar Research
• /
• v.26 no.3
• /
• pp.489-499
• /
• 2004

The climatological characteristics of the polar ionospheric currents obtained from the simultaneous observations of the ionospheric electric field and conductivity are examined. For this purpose, 43 and 109 days of measurements from the Chatanika and Sondrestrom incoherent scatter radars are utilized respectively. The ionospheric current density is compared with the corresponding ground magnetic disturbance. Several interesting characteristics about the polar ionosphere are apparent from this study: (1) The sun determines largely the conductance over the Sondrestrom radar, while the nighttime conductance distribution over the Chatanika radar is significantly affected by auroral precipitation. (2) The regions of the maximum N-S electric field over the Chatanika radar are located approximately at the dawn and dusk sectors, while they tend to shift towards dayside over the Sondrestrom radar. The N-S component over Son-drestrom is slightly stronger than Chatanika. However, the E-W component over Chatanika is negligible compared to that of Sondrestrom. (3) The E-W ionospheric current flows dominantly in the night hemisphere over Chatanika, while it flows in the sunlit hemisphere over Sondrestrom. The N-S current over Chatanika flows prominently in the dawn and dusk sectors, while a strong southward current flows in the prenoon sector over Sondrestrom. (4) The assumption of infinite sheet current approximation is far from realistic, underestimating the current density by a factor of 2 or more. It is particularly serious for the higher latitude region. (5) The correlation between ${\Delta}H\;and\;J_E$ is higher than the one between ${\Delta}D\;and\;J_N$, indicating that field-aligned current affects ${\Delta}D$significantly.

• Proceedings of the KOSOMBE Conference
• /
• v.1998 no.11
• /
• pp.249-250
• /
• 1998

In this study, we modeled moving-actuator type Total Artificial Heart (TAH) with cardiovascular system as a form of electric circuit. The bronchial circulation, important for the imbalance between the left cardiac output and the right one, was considered and added to the model. In the model, the relations of hemodynamic variables, just as blood pressures, volumes, or flow rates of each part of body, can be expressed as simultaneous first order ordinary differential equations. To solve the equations by the numerical analysis, Runge-Kutta forth order approximation method was adopted. The simulation software (SimTAH), implemented in C++ as a window-based application program, was developed to display the hemodynamic variables and to receive control inputs from users. SimTAH was evaluated by comparison of the simulation results with the results of mock-circulation tests, in vitro.

• Journal of the Society of Naval Architects of Korea
• /
• v.34 no.4
• /
• pp.31-41
• /
• 1997

This paper describes a potential-based panel method formulated for the prediction of the steady performance of a waterjet propulsor. The method employs normal dipoles and sources distributed on the solid surfaces such as the impeller/stator blades, hub and duct, and normal dipoles in the shed wakes trailing the impeller and stator to represent the potential flow around the waterjet propulsor. To define a closed boundary surface, the inlet and outlet open boundary surfaces are introduced where the sources and dipoles are distributed. The kinematic boundary condition on the solid boundary surface is satisfied by requiring that the normal component of the total velocity should vanish. On the inlet surface, the total inflow flux into the duct is specified, and on the outlet surface the conservation of mass principle is applied to evaluate the source strength. The solid surfaces are discretized into a set of quadrilateral panel elements and the strengths of sources and dipoles are assumed constant at each panel. Applying this approximation to the boundary conditions leads to a set of simultaneous equations. Systematic numerical tests show that the present numerical method is fast and stable. In order to validate the present method, sample computations are carried out first for the case of a conventional axial flow fan which has a similar geometry as the waterjet propulsor, and then for the case of a waterjet propulsor on which experiments are carried out at KRISO(Korea Research Institute of Ships and Ocean Engineering).