• Nuclear Engineering and Technology
• /
• v.41 no.7
• /
• pp.885-892
• /
• 2009

This manuscript will discuss a numerical method where the six equations of two-phase flow, the solid heat conduction equations, and the two equations that describe neutron diffusion and precursor concentration are solved together in a tightly coupled, nonlinear fashion for a simplified model of a nuclear reactor core. This approach has two important advantages. The first advantage is a higher level of accuracy. Because the equations are solved together in a single nonlinear system, the solution is more accurate than the traditional "operator split" approach where the two-phase flow equations are solved first, the heat conduction is solved second and the neutron diffusion is solved third, limiting the temporal accuracy to $1^{st}$ order because the nonlinear coupling between the physics is handled explicitly. The second advantage of the method described in this manuscript is that the time step control in the fully implicit system can be based on the timescale of the solution rather than a stability-based time step restriction like the material Courant limit required of operator-split methods. In this work, a pilot code was used which employs this tightly coupled, fully implicit method to simulate a reactor core. Results are presented from a simulated control rod movement which show $2^{nd}$ order accuracy in time. Also described in this paper is a simulated rod ejection demonstrating how the fastest timescale of the problem can change between the state variables of neutronics, conduction and two-phase flow during the course of a transient.

• Journal of Hydrogen and New Energy
• /
• v.33 no.6
• /
• pp.692-698
• /
• 2022

With the expansion of domestic hydrogen fuelling station infrastructure, it is necessary to secure reliability among hydrogen traders, and for this, technology to accurately measure hydrogen is important. In this study, 4 types of hydrogen trading volume calculation models (model 1-4) were presented to improve the accuracy of the hydrogen trading volume. In order to obtain the reference value of model 4, and experiment was conducted using a flow rate measurement equipment, and the error rate of the calculated value for each model was compared and analyzed. As a result, model 1 had the lowest metering accuracy, model 2 had the second highest metering accuracy and model 3 had the highest metering accuracy until a certain point. But after the point, model 2 had the highest metering accuracy and model 3 had the second metering accuracy.

• Ocean Systems Engineering
• /
• v.9 no.1
• /
• pp.97-109
• /
• 2019

Accurate calculation of the mean drift forces and moments is necessary when studying the higher order excitations on the floater in waves. When taking the time average of the second order forces and moments, the second order potential and motion diminish with only the first order terms remained. However, in the results of the first order forces or motions, the irregular frequency effects are often observed in higher frequencies, which will affect the accuracy of the calculation of the second order forces and moments. Therefore, we need to pay close attention to the irregular frequency effects in the mean drift forces. This paper will discuss about the irregular frequency effects in the calculations of the mean drift forces and validate our in-house program MDL Multi DYN using some examples which are known to have irregular frequency effects. Finally, we prove that it is necessary to remove the effects and demonstrate that the effectiveness of the formula and methods adopted in the development of our program.

• Steel and Composite Structures
• /
• v.29 no.2
• /
• pp.175-186
• /
• 2018

The plastic hinge method and the plastic zone method are extensively adopted in displacement-based elements and force-based elements respectively for second-order inelastic analysis. The former enhances the computational efficiency with relatively less accurate results while the latter precisely predicts the structural behavior but generally requires more computer time. The displacement-based elements receive criticism mainly on plasticity dominated problems not only in accuracy but also in longer computer time to redistribute the forces due to formation of plastic hinges. The multi-element-per-member model relieves this problem to some extent but will induce a new problem in modeling of member initial imperfections required in design codes for direct analysis. On the contrary, a force-based element with several integration points is sufficient for material yielding. However, use of more integration points or elements associated with fiber section reduces computational efficiency. In this paper, a new force-based element equipped with stress-resultant plasticity model with minimal computational cost is proposed for second-order inelastic analysis. This element is able to take the member initial bowing into account such that one-element-per-member model is adequate and complied with the codified requirements of direct analysis. This innovative solution is new and practical for routine design. Finally, several examples demonstrate the validity and accuracy of the proposed method.

• KSII Transactions on Internet and Information Systems (TIIS)
• /
• v.2 no.1
• /
• pp.51-62
• /
• 2008

This paper proposes a second-order implicit constraint enforcement method which yields enhanced controllability compared to a first-order implicit constraints enforcement method. Although the proposed method requires solving a linear system twice, it yields superior accuracy from the constraints error perspective and guarantees the precise and natural movement of objects, in contrast to the first-order method. Thus, the proposed method is the most suitable for exact prediction simulations. This paper describes the numerical formulation of second-order implicit constraints enforcement. To prove its superiority, the proposed method is compared with the firstorder method using a simple two-link simulation. In this paper, there is a reasonable discussion about the comparison of constraints error and the analysis of dynamic behavior using kinetic energy and potential energy.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.08a
• /
• pp.31-36
• /
• 2003

A field reconstruction scheme for a cell centered finite volume method on unstructured meshes is developed. Regardless of mesh quality, this method is exact within a machine accuracy if the solution is linear, which means it has full second order accuracy. It does not have any limitation on cell shape except convexity of the cells and recovers standard discretization stencils at structured orthogonal grids. Accuracy comparisons with other popular reconstruction schemes are performed on a simple example.

• Management Science and Financial Engineering
• /
• v.21 no.2
• /
• pp.25-30
• /
• 2015

This short article compares different solution methods for a basic RBC model (Hansen, 1985). We solve and simulate the model using two main algorithms: the methods of perturbation and projection, respectively. One novelty is that we offer a type of the hybrid method: we compute easily a second-order approximation to decision rules and use that approximation as an initial guess for finding Chebyshev polynomials. We also find that the second-order perturbation method is most competitive in terms of accuracy for standard RBC model.

• Proceedings of the KIEE Conference
• /
• 2003.07e
• /
• pp.55-57
• /
• 2003

This paper presents a new method for finding the Block Pulse series coefficients and deriving the Block Pulse integration operational matrices which are necessary for the control fields using the Block Pulse functions. In this paper, the accuracy of the Block Pulse series coefficients derived by using the Lagrange second order interpolation polynomial is approved by the mathematical method.

• 한국전산유체공학회:학술대회논문집
• /
• 2000.10a
• /
• pp.129-134
• /
• 2000

An efficient numerical method to solve the unsteady incompressible Navier-Stokes equations is developed. A fully implicit time advancement is employed to avoid the CFL(Courant-Friedrichs-Lewy) restriction, where the Crank-Nicholson discretization is used for both the diffusion and convection terms. Based on a block LU decomposition, velocity-pressure decoupling is achieved in conjunction with the approximate factorization. Main emphasis is placed on the additional decoupling of the intermediate velocity components with only n th time step velocity The temporal second-order accuracy is Preserved with the approximate factorization without any modification of boundary conditions. Since the decoupled momentum equations are solved without iteration, the computational time is reduced significantly. The present decoupling method is validated by solving the turbulent minimal channel flow unit.

• Journal of computational fluids engineering
• /
• v.6 no.3
• /
• pp.19-26
• /
• 2001

An Eulerian-Lagrangian method, so called immersed boundary method, is used for analysing viscous flow around arbitrary bodies, where governing equations are discretized on a regular grid by using a finite volume method. To improve the accuracy of flow near body boundaries, a second-order accurate interpolation scheme is used and a level-set based grid deformation method is presented to construct the adaptive grids around body boundaries. The present scheme is used to simulate steady flow around a semicircular cylinder mounted on the bottom of flow domain and calculated results are validated by results of a body fitted grid method. Finally, present method is applied to a complex flow around multi body and the usefulness is checked by investigating calculated results.