• Nuclear Engineering and Technology
• /
• 제52권2호
• /
• pp.223-229
• /
• 2020

The present work proposes a solution to the static Boltzmann transport equation approximated by the simplified P₃ (SP₃) on angular, and the analytic coarse mesh finite difference (ACMFD) for spatial variables. Multi-group SP₃-ACMFD equations in 3D rectangular geometry are solved using the GMRES solution technique. As the core time dependent analysis necessitates the solution of an eigenvalue problem for an initial condition, this work is hence devoted to development and verification of the proposed static SP₃-ACMFD solver. A 3D multi-group static diffusion solver is also developed as a byproduct of this work to assess the improvement achieved using the SP₃ technique. Static results are then compared against transport benchmarks to assess the proximity of SP₃-ACMFD solutions to their full transport peers. Results prove that the approach can be considered as an acceptable interim approximation with outputs superior to the diffusion method, close to the transport results, and with the computational costs less than the full transport approach. The work would be further generalized to time dependent solutions in Part II.

• Nuclear Engineering and Technology
• /
• 제55권6호
• /
• pp.2172-2194
• /
• 2023

A variational nodal method (VNM) with unstructured-mesh is presented for solving steady-state and dynamic neutron diffusion equations. Orthogonal polynomials are employed for spatial discretization, and the stiffness confinement method (SCM) is implemented for temporal discretization. Coordinate transformation relations are derived to map unstructured triangular nodes to a standard node. Methods for constructing triangular prism space trial functions and identifying unique nodes are elaborated. Additionally, the partitioned matrix (PM) and generalized partitioned matrix (GPM) methods are proposed to accelerate the within-group and power iterations. Neutron diffusion problems with different fuel assembly geometries validate the method. With less than 5 pcm eigenvalue (k_eff) error and 1% relative power error, the accuracy is comparable to reference methods. In addition, a test case based on the kilowatt heat pipe reactor, KRUSTY, is created, simulated, and evaluated to illustrate the method's precision and geometrical flexibility. The Dodds problem with a step transient perturbation proves that the SCM allows for sufficiently accurate power predictions even with a large time-step of approximately 0.1 s. In addition, combining the PM and GPM results in a speedup ratio of 2-3.

• Wind and Structures
• /
• 제37권1호
• /
• pp.57-78
• /
• 2023

In this study, a generalized three-degree-of-freedom (3-DoF) analytical model is formulated to predict linear aerodynamic instabilities of a prism under quasi-steady (QS) conditions. The prism is assumed to possess a generic cross-section exposed to turbulent wind flow. The 3-DoFs encompass two orthogonal horizontal directions and rotation about the prism body axis. Inertial coupling is considered to account for the non-coincidence of the mass center and the rotation center. The aerodynamic force coefficients-drag, lift, and moment-depend on the Reynolds number based on relative flow velocity, angle of attack, and the angle between the wind and the cable. Aerodynamic forces are linearized with respect to the static equilibrium configuration and mean wind velocity. Routh-Hurwitz and Liénard and Chipart criteria are used in the eigenvalue problem, yielding an analytical solution for instabilities in galloping and static divergence types. Additionally, the minimum structural damping and stiffness required to prevent these instabilities are numerically determined. The proposed 3-DoF instability model is subsequently applied to a conductor with ice accretion and a full-scale dry inclined cable. In comparison to existing models, the developed model demonstrates superior prediction accuracy for unstable regions compared with results in wind tunnel tests.

• Steel and Composite Structures
• /
• 제35권5호
• /
• pp.671-685
• /
• 2020

This manuscript presents impacts of gradation of material functions and axial load functions on critical buckling loads and mode shapes of functionally graded (FG) thin and thick beams by using higher order shear deformation theory, for the first time. Volume fractions of metal and ceramic materials are assumed to be distributed through a beam thickness by both sigmoid law and symmetric power functions. Ceramic-metal-ceramic (CMC) and metal-ceramic-metal (MCM) symmetric distributions are proposed relative to mid-plane of the beam structure. The axial compressive load is depicted by constant, linear, and parabolic continuous functions through the axial direction. The equilibrium governing equations are derived by using Hamilton's principles. Numerical differential quadrature method (DQM) is developed to discretize the spatial domain and covert the governing variable coefficients differential equations and boundary conditions to system of algebraic equations. Algebraic equations are formed as a generalized matrix eigenvalue problem, that will be solved to get eigenvalues (buckling loads) and eigenvectors (mode shapes). The proposed model is verified with respectable published work. Numerical results depict influences of gradation function, gradation parameter, axial load function, slenderness ratio and boundary conditions on critical buckling loads and mode-shapes of FG beam structure. It is found that gradation types have different effects on the critical buckling. The proposed model can be effective in analysis and design of structure beam element subject to distributed axial compressive load, such as, spacecraft, nuclear structure, and naval structure.

• 한국강구조학회 논문집
• /
• 제14권4호
• /
• pp.509-518
• /
• 2002

비대칭 단면을 갖는 박벽보의 3차원 휨-배틂 좌굴해석 및 정적해석을 위하여, 평형방정식과 힘-변위 관계식을 이용하여 엄밀한 정적요소강성행렬을 수치적으로 산정하는 개선된 기법을 제시한다. 먼저 14개의 변위피라미터를 도입하여 고차의 연립미분방정식을 1차 연립미분방정식으로 변환하고, 복소수 영역에서 선형고유치문제를 해를 구한다. 이 경우 동적강성행렬을 산정하는 경우와는 달리 복수개의 '영'의 고유치가 발생한다. 이에 대응하는 변위피라미터의 다항식을 항등식 조거능로부터 구하고, 이를 고유치와 결합하여 박벽보 요소의 엄밀한 처짐함수를 구한다. 이렇게 구한 엄밀한 처짐함수에 재단력-변위 관계식을 적용하여 세가지 초기단면력 조건에 대응하는 엄밀한 정적요소강성행렬을 산정한다. 본 방법의 타당성을 보이기 위하여 비대칭 박벽보의 좌굴하중과 처짐값을 계산하고 해석해나 ABAQUS 쉘요소를 이용한 해석결과 및 직선보요소를 사용한 유한요소해의 결과와 비교, 검증한다.

• 한국전산구조공학회논문집
• /
• 제15권3호
• /
• pp.463-469
• /
• 2002

탄성지반 위에 놓인 보-기둥 요소의 총포텐셜 에너지로부터 변분원리를 적용하여 지배방정식과 힘-변위 관계식을 유도하였다. 4계 상미분방정식 형태의 지배방정식을 4개의 변위 파라메타를 도입하여 1계 연립미분방정식 형태의 선형 고유치 문제로 전환하고, 힘-변위 관계식을 적용하여 엄밀한 정적, 동적 요소강성행렬을 유도하였다. 직접강성법을 이용하여 구조물 강성행렬을 구하고, 2차원 보-기둥구조의 엄밀한 좌굴하중과 고유진동수를 구하고, 결과를 유한요소해와 비교함으로써 본 연구의 타당성을 검증하였다. 이러한 엄밀한 해석방법은 Hermitian 다항식을 형상함수로 도입하여 요소의 강성행렬을 산정하는 유한요소법과 비교할 때, 요소의 수를 대폭 줄일 수 있는 장점이 있다.