• Nuclear Engineering and Technology
• /
• v.53 no.7
• /
• pp.2084-2094
• /
• 2021

Delta-form-based methods for solving high order spatial discretization schemes are introduced into the reactor S_N transport equation. Due to the nature of the delta-form, the final numerical accuracy only depends on the residuals on the right side of the discrete equations and have nothing to do with the parts on the left side. Therefore, various high order spatial discretization methods can be easily adopted for only the transport term on the right side of the discrete equations. Then the simplest step or other robust schemes can be adopted to discretize the increment on the left hand side to ensure the good iterative convergence. The delta-form framework makes the sweeping and iterative strategies of various high order spatial discretization methods be completely the same with those of the traditional S_N codes, only by adding the residuals into the source terms. In this paper, the flux limiter method and weighted essentially non-oscillatory scheme are used for the verification purpose to only show the advantages of the introduction of delta-form-based solving methods and other high order spatial discretization methods can be also easily extended to solve the S_N transport equations. Numerical solutions indicate the correctness and effectiveness of delta-form-based solving method.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1994.04a
• /
• pp.1-8
• /
• 1994

An automated procedure which allows adaptation of spatial and time discretization simultaneously in finite element analysis of linear elastodynamic problems is presented. For dynamic problems having responses dominated by high frequency modes, such as those with impact, explosive, traveling and earthquake loads high gradient stress regions change their locations from time to time. And the time step size may need to vary in order to deal with whole process ranging from transient phase to steady state phase. As the sizes of elements in space vary in different regions, the procedure also permits different time stepping. In such a way, the best performance attainable by the finite element method can be achieved. In this study, we estimate both of the kinetic energy error and stran energy error induced by spatial and time discretization in a consistent manner. Numerical examples are used to demonstrate the performance of the procedure.

• Nuclear Engineering and Technology
• /
• v.53 no.12
• /
• pp.3861-3878
• /
• 2021

The steady-state simplified spherical harmonics equations (SP_N equations) are a higher order approximation to the neutron transport equations than the neutron diffusion equation that also have reasonable computational demands. This work extends these results for the analysis of transients by comparing of two formulations of time-dependent SP_N equations considering different treatments for the time derivatives of the field moments. The first is the full system of equations and the second is a diffusive approximation of these equations that neglects the time derivatives of the odd moments. The spatial discretization of these methodologies is made by using a high order finite element method. For the time discretization, a semi-implicit Euler method is used. Numerical results show that the diffusive formulation for the time-dependent simplified spherical harmonics equations does not present a relevant loss of accuracy while being more computationally efficient than the full system.

• Nuclear Engineering and Technology
• /
• v.54 no.7
• /
• pp.2625-2634
• /
• 2022

Discretization errors are extremely challenging conundrums of discrete ordinates calculations for radiation transport problems with void regions. In previous work, we have presented a multi-collision source method (MCS) to overcome discretization errors, but the efficiency needs to be improved. This paper proposes a goal-oriented algorithm for the MCS method to adaptively determine the partitioning of the geometry and dynamically change the angular quadrature in remaining iterations. The importance factor based on the adjoint transport calculation obtains the response function to get a problem-dependent, goal-oriented spatial decomposition. The difference in the scalar fluxes from one high-order quadrature set to a lower one provides the error estimation as a driving force behind the dynamic quadrature. The goal-oriented algorithm allows optimizing by using ray-tracing technology or high-order quadrature sets in the first few iterations and arranging the integration order of the remaining iterations from high to low. The algorithm has been implemented in the 3D transport code ARES and was tested on the Kobayashi benchmarks. The numerical results show a reduction in computation time on these problems for the same desired level of accuracy as compared to the standard ARES code, and it has clear advantages over the traditional MCS method in solving radiation transport problems with reflective boundary conditions.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.22 no.6
• /
• pp.72-83
• /
• 2018

In this study, a series of CFD analyses were carried out for a hydrogen rocket combustor with a single shear coaxial injector. A hybrid RANS/LES approach was used for the turbulent combustion analysis with a two-dimensional axisymmetric configuration. Three reaction mechanisms, three spatial discretization methods, and three levels of grid resolution were compared to determine an appropriate CFD approach. The performance of the CFD prediction were investigated by comparing the wall heat flux with experimental data. Investigation of the flow field results provides an insight into the characteristics of the turbulent reacting flow of a rocket combustor with a shear coaxial injector.

• The Transactions of the Korean Institute of Electrical Engineers D
• /
• v.50 no.10
• /
• pp.464-472
• /
• 2001

This study uses distributed parameter systems as the spatial discretization technique, modelling in lumped parameter systems, and applies fast WALSH transform and the Picard's iteration method to high order partial differential equations and matrix partial differential equations. This thesis presents a new algorithm which usefully exercises the optimal control in the distributed parameter systems. In exercising optimal control of distributed parameter systems, excellent consequences are found without using the existing decentralized control or hierarchical control method. This study will help apply to linear time-varying systems and non-linear systems. Further research on algorithm will be required to solve the problems of convergence in case of numerous applicable intervals.

• Proceedings of the KIEE Conference
• /
• 2001.07d
• /
• pp.2084-2086
• /
• 2001

This study uses distributed parameter systems as the spatial discretization technique, modelling in lumped parameter systems, and applies orthogonal transform and the Picard's iteration method to high order partial differential equations and matrix partial differential equations. In exercising optimal control of distributed parameter systems, excellent consequences are found without using the existing decentralized control or hierarchical control method.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.386-391
• /
• 2001

A finite element code based on P2P1 tetra element has been developed for the large eddy simulation (LES) of turbulent flows around a complex geometry. Fractional 4-step algorithm is employed to obtain time accurate solution since it is less expensive than the integrated formulation, in which the velocity and pressure fields are solved at the same time. Crank-Nicolson method is used for second order temporal discretization and Galerkin method is adopted for spatial discretization. For very high Reynolds number flows, which would require a formidable number of nodes to resolve the flow field, SUPG (Streamline Upwind Petrov-Galerkin) method is applied to the quadratic interpolation function for velocity variables, Noting that the calculation of intrinsic time scale is very complicated when using SUPG for quadratic tetra element of velocity variables, the present study uses a unique intrinsic time scale proposed by Codina et al. since it makes the present three-dimensional unstructured code much simpler in terms of implementing SUPG. In order to see the effect of numerical diffusion caused by using an upwind scheme (SUPG), those obtained from P2P1 Galerkin method and P2P1 Petrov-Galerkin approach are compared for the flow around a sphere at some Reynolds number. Smagorinsky model is adopted as subgrid scale models in the context of P2P1 finite element method. As a benchmark problem for code validation, turbulent flows around a sphere and a MIRA model have been studied at various Reynolds numbers.