• The Bulletin of The Korean Astronomical Society
• /
• v.46 no.1
• /
• pp.30.3-30.3
• /
• 2021

In order to generate the 3D structure of the 21-cm signal during the reionization, semi-numerical simulations based on Excursion set formalism are broadly used. However, semi-numerical simulations in the realization of the 3D structure are known to be the ionizing photon non-conserving by the structure of the Excursion set approach. Recently, explicit photon conserving algorithms for semi-numerical simulations introduced, but they are still too slow when forward modelling the 21-cm signal with high-dimensional parameter spaces. Here, we introduce a new method for approximately correcting photon non-conservation, which can be applied on-the-fly. This method is tailored towards the efficient simulation and Bayesian inference with high-dimensional parameter space. Then, we investigate how large an impact that photon non-conservation has on astrophysical parameter inference by performing an MCMC analysis. We find that the ionizing escape parameter is deviated from the fiducial value by 2 sigma when we infer astrophysical parameters without this correction.

• Nuclear Engineering and Technology
• /
• v.53 no.1
• /
• pp.30-43
• /
• 2021

In the 2D/1D method, a global adjoint CMFD based on the generalized equivalence theory is built to synthesize the 2D radial MOC adjoint and 1D axial NEM adjoint calculation and also to accelerate the iteration convergence of 3D whole-core adjoint transport calculation. Even more important, an advanced yet accurate two-level (TL) CMFD acceleration technique is proposed, in which an equivalent one-group adjoint CMFD is established to accelerate the multi-group adjoint CMFD and then to accelerate the 3D whole-core adjoint transport calculation efficiently. Based on these method, a new code is developed to perform 3D adjoint neutron flux calculation. Then a set of VERA and C5G7 benchmark problems are chosen to verify the capability of the 3D adjoint calculations and the effectiveness of TL CMFD acceleration. The numerical results demonstrate that acceptable accuracy of 2D/1D adjoint calculations and superior acceleration of TL CMFD are achievable.

• Wind and Structures
• /
• v.5 no.2_3_4
• /
• pp.379-392
• /
• 2002

Numerical flow computations around an aeroelastic 3D square cylinder immersed in the turbulent boundary layer are shown. Present computational code can be characterized by three numerical aspects which are 1) the method of artificial compressibility is adopted for the incompressible flow computations, 2) the domain decomposition technique is used to get better grid point distributions, and 3) to achieve the conservation law both in time and space when the flow is computed a with moving and transformed grid, the time derivatives of metrics are evaluated using the time-and-space volume. To provide time-dependant inflow boundary conditions satisfying prescribed time-averaged velocity profiles, a convenient way for generating inflow turbulence is proposed. The square cylinder is modeled as a 4-lumped-mass system and it vibrates with two-degree of freedom of heaving motion. Those blocks which surround the cylinder are deformed according to the cylinder's motion. Vigorous oscillations occur as the vortex shedding frequency approaches cylinder's natural frequencies.

• Journal of Power System Engineering
• /
• v.2 no.3
• /
• pp.21-26
• /
• 1998

Numerical simulations for various fluid flows require enormous computing time during iterations. In order to solve this problem, several techniques have been proposed. A SOR method is one of the effective methods for solving elliptic equations. However, it is very difficult to find the optimum relaxation factor, the value of this factor for practical problems used to be estimated on the basis of expertise. In this paper, the implication of the relaxation factor are translated into fuzzy control rules on the basis of the expertise of numerical analysers, and fuzzy controller incorporated into a numerical algorithm. From two cases of study, Poisson equation and cavity flow problem, we confirmed the possibility of computational acceleration with fuzzy logic and qualitative reasoning in numerical simulations. Numerical experiments with the fuzzy controller resulted in generating a good performance.

• Geomechanics and Engineering
• /
• v.28 no.5
• /
• pp.437-454
• /
• 2022

The hydro-mechanical (HM) two-way sequential coupling in the TOUGH2-FLAC3D linking algorithm is validated completely and successfully in both M to H and H to M directions, which are initiated by mechanical surface loading for geomechanical validation and hydrological groundwater pumping for hydrogeological validation, respectively. For such complete and successful validation, a TOUGH2-FLAC3D linked numerical model is developed first by adopting the TOUGH2-FLAC3D linking algorithm, which uses the two-way (fixed-stress split) sequential coupling scheme and the implicit backward time stepping method. Two geomechanical and two hydrogeological validation problems are then simulated using the linked numerical model together with basic validation strategies and prerequisites. The second geomechanical and second hydrogeological validation problems are also associated with the Mandel effect and the Noordbergum and Rhade effects, respectively, which are three phenomenally well-known but numerically challenging HM effects. Finally, sequentially coupled numerical solutions are compared with either analytical solutions (verification) or fully coupled numerical solutions (benchmarking). In all the four validation problems, they show almost perfect to extremely or very good agreement. In addition, the second geomechanical validation problem clearly displays the Mandel effect and suggests a proper or minimum geometrical ratio of the height to the width for the rectangular domain to maximize agreement between the numerical and analytical solutions. In the meantime, the second hydrogeological validation problem clearly displays the Noordbergum and Rhade effects and implies that the HM two-way sequential coupling scheme used in the linked numerical model is as rigorous as the HM two-way full coupling scheme used in a fully coupled numerical model.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.27-30
• /
• 2002

As computer capacity has been progressed continuously, the studies of the flow characteristics have been performing by the numerical methods actively. In this study, 3-dimensional unsteady incompressible Wavier-Stokes equation was solved by numerical method using the fractional step method with the fourth order compact pade' scheme to achieve high accuracy To validate the present code and algorithm, 3D flow-field around a cylinder was simulated. The drag coefficient and lift coefficient were computed and, then, compared with experiment. The present code will be tailored to LES simulation for more accurate turbulent flow analysis.

• 제어로봇시스템학회:학술대회논문집
• /
• 2005.06a
• /
• pp.2410-2413
• /
• 2005

Finite element method (FEM) has been widely used as a useful numerical method that can analyze complex engineering problems in electro-magnetics, mechanics, and others. CEMTool, which is similar to MATLAB, is a command style design and analyzing package for scientific and technological algorithm and a matrix based computation language. In this paper, we present new 3D FEM package in CEMTool environment. In contrast to the existing CEMTool 2D FEM package and MATLAB PDE (Partial Differential Equation) Toolbox, our proposed 3D FEM package can deal with complex 3D models, not a cross-section of 3D models. In the pre-processor of 3D FEM package, a new 3D mesh generating algorithm can make information on 3D Delaunay tetrahedral mesh elements for analyses of 3D FEM problems. The solver of the 3D FEM package offers three methods for solving the linear algebraic matrix equation, i.e., Gauss-Jordan elimination solver, Band solver, and Skyline solver. The post-processor visualizes the results for 3D FEM problems such as the deformed position and the stress. Consequently, with our new 3D FEM toolbox, we can analyze more diverse engineering problems which the existing CEMTool 2D FEM package or MATLAB PDE Toolbox can not solve.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.29 no.5
• /
• pp.471-482
• /
• 2016

In this study, we focused on a speed-up of KIESSI-3D program, which is based on FE-IE techniques, by introducing a p-version dynamic infinite element method. In order to evaluate performance of the KIESSI-3D, numerical analyses for eight real-scale SSI problems are carried out. We considered three types of KIESSI-3D numerical models whose radii of near-field soil region($r_0$)are 1.2, 1.5, and 3.0 times of basemat radius of structure(R). In addition, SSI analyses using the SASSI2010 program are carried out used for comparison of accuracy and runtime against those of the KIESSI-3D. Numerical results show that the KIESSI-3D model of $r_0=1.2R$ is enough to give accurate solution. In view of the computing speed, the new KIESSI-3D was up to 25 times faster than the old KIESSI-3D.

• The Transactions of the Korean Institute of Electrical Engineers B
• /
• v.51 no.8
• /
• pp.449-456
• /
• 2002

3D Equivalent magnetic Circuit Network Method (EMCNM) is comparatively the easy way that analyzes 3D models of Electric Machine by using permeance as a distributive magnetic circuit parameter under the existing magnetic equivalent circuit method and Numerical Method. The existing 3D EMCNM could not correctly describe the shape of an analysis target when using rectangular shape element or fan shape element, so it made errors when calculating permeance. Therefore, this paper proposes the trapezoidal element contained rectangular element, fan-shape element, and quadrilateral element to express a shape. The proposed method in this research was confirmed as a useful and an accurate method through comparing with the analysis result of SRM model that is sufficiently guaranteed by 2D-Analysis.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.29 no.2
• /
• pp.149-159
• /
• 2016

In this study, the solitary wave-induced tsunami force acting on an onshore bridges in coastal area was numerically modelled by means of TWOPM-3D based on Navier-Stokes solver and VOF method which can track free surface effectively. The validity of numerical analysis was verified by comparing the experimental tsunami bore force acting on vertical wall and column structure. In particular, the characteristics of tsunami force with the changing tsunami intensity were surveyed through numerical experiments. The availability of 3-dimensional numerical analysis was reviewed through the comparison between the existing numerical results and design criteria for each drag force coefficient by applying Morison equation considering only drag force. As reasonable and high-precision estimation method of tsunami force, it was suggested to apply the estimation method taking drag and inertial force into consideration at the same time.