• Nuclear Engineering and Technology
• /
• v.49 no.6
• /
• pp.1358-1367
• /
• 2017

The data integration with modeled predictions (DIMP) model is a promising inverse radiation transport method for solving the special nuclear material (SNM) holdup problem. Unlike previous methods, DIMP is a completely passive nondestructive assay technique that requires no initial assumptions regarding the source distribution or active measurement time. DIMP predicts the most probable source location and distribution through Bayesian inference and quasi-Newtonian optimization of predicted detector responses (using the adjoint transport solution) with measured responses. DIMP performs well with forward hemispherical collimation and unshielded measurements, but several considerations are required when using narrow-view collimated detectors. DIMP converged well to the correct source distribution as the number of synthetic responses increased. DIMP also performed well for the first experimental validation exercise after applying a collimation factor, and sufficiently reducing the source search volume's extent to prevent the optimizer from getting stuck in local minima. DIMP's simple point detector response function (DRF) is being improved to address coplanar false positive/negative responses, and an angular DRF is being considered for integration with the next version of DIMP to account for highly collimated responses. Overall, DIMP shows promise for solving the SNM holdup inverse problem, especially once an improved optimization algorithm is implemented.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.23 no.11
• /
• pp.943-950
• /
• 2013

This paper discusses a harmonic response estimation method on the L$\acute{e}$vy plate with two opposite edges simply supported and the other two edges having free boundary conditions. Since the equation of motion of the plate is not self-adjoint, the modes are not orthogonal to each other on the domain. Noting that the L$\acute{e}$vy plate can be expressed using one term sinusoidal function that is orthogonal to other sinusoidal functions, this paper suggested the calculation method that is equivalent to finding a least square error minimization solution of the finite number of algebraic equations. Example problems subjected to a distributed area loading and a distributed line loading are defined and their solutions are provided. The solutions are compared to those of the commercial code, ANSYS. According to the verification results, it is expected that the suggested method will be useful to predict the forced response on the L$\acute{e}$vy plate with the distributed area or line loading conditions.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.20 no.10
• /
• pp.915-922
• /
• 2010

This paper presents an efficient method for determining the forced response of a spinning flexible disk-spindle system supported by fluid dynamic bearings(FDBs) in a computer hard disk drive(HDD). The spinning flexible disk-spindle system is represented by the asymmetric finite element equations of motion originating from the asymmetric dynamic coefficients of the FDBs and the gyroscopic moment of a spinning disk-spindle system. The proposed method utilizes only the right eigenvectors of the eigenvalue problem to transform the large asymmetric finite element equations of motion into a small number of coupled equations, guaranteeing the accuracy of their numerical integration. The results are then back-substituted into the equations of motion to determine the forced response. The effectiveness of the proposed method was verified by comparing it with the responses from the classical methods of mode superposition with the general eigenvalue problems, and mode superposition with modal approximation. The proposed method was shown to be effective in determining the forced response represented by the asymmetric finite element equations of motion of a spinning flexible disk-spindle system supported by FDBs.

• Steel and Composite Structures
• /
• v.30 no.1
• /
• pp.57-67
• /
• 2019

This paper deals with the static analysis of axially functionally graded rectangular plates. It is assumed that the flexural rigidity of the plate varies exponentially along one of the plate's in-plane dimensions. Both an analytical approach and a numerical method are utilized to solve the problem. The analytical solution is obtained by using the Green's function method. To employ this approach, the adjoint boundary value problem is established. Then, exact solutions for deflection of the plate for different boundary conditions are found. In another way, a finite element formulation for the problem is developed. In order to demonstrate the validity of the Authors' formulation, the results obtained via both mentioned schemes are compared with each other for functionally graded plates and with results of previously published works for homogeneous plates. The effect of plate parameters on the response of the plate is also investigated. To remind the research background, a brief review on the application of Green's function method in plates' analysis and functionally graded plates is also presented.

• 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.

• Nuclear Engineering and Technology
• /
• v.53 no.5
• /
• pp.1391-1402
• /
• 2021

The aim of this work is to prepare a neutron noise calculator based on the second order of average current nodal expansion method (ACNEM). Generally, nodal methods have the ability to fulfill the neutronic analysis with adequate precision using coarse meshes as large as a fuel assembly size. But, for the zeroth order of ACNEM, the accuracy of neutronic simulations may not be sufficient when coarse meshes are employed in the reactor core modeling. In this work, the capability of second order ACNEM is extended for solving the neutron diffusion equation in the frequency domain using coarse meshes. For this purpose, two problems are modeled and checked including a slab reactor and 2D BIBLIS PWR. For validating of results, a semi-analytical solution is utilized for 1D test case, and for 2D problem, the results of both forward and adjoint neutron noise calculations are exploited. Numerical results indicate that by increasing the order of method, the errors of frequency dependent coarse mesh solutions are considerably decreased in comparison to the reference. Accordingly, the accuracy of second order ACNEM can be acceptable for the neutron noise calculations by using coarse meshes in the nuclear reactor core.

• Atmosphere
• /
• v.32 no.3
• /
• pp.191-206
• /
• 2022

Assimilation trials were performed using the Korea Institute of Atmospheric Prediction Systems (KIAPS) Korea Integrated Model (KIM) semi-operational forecast system to assess the impact of ground-based Global Navigation Satellite System (GNSS) Zenith Total Delay (ZTD) on forecast. To use the optimal observation in data assimilation of KIM forecast system, in this study, the ZTD observation were pre-processed. It involves the bias correction using long term background of KIM, the quality control based on background and the thinning of ZTD data. Also, to give the effect of observation directly to data assimilation, the observation operator which include non-linear model, tangent linear model, adjoint model, and jacobian code was developed and verified. As a result, impact of ZTD observation in both analysis and forecast was neutral or slightly positive on most meteorological variables, but positive on geopotential height. In addition, ZTD observations contributed to the improvement on precipitation of KIM forecast, specially over 5 mm/day precipitation intensity.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2021.06a
• /
• pp.433-433
• /
• 2021

본 연구는 개수로 흐름에서 조류발전단지의 터빈 최적 배열의 거시적 특성에 관한 연구를 수행하였다. 천수방정식을 통해 직사각형 개수로의 흐름장을 해석하였고, 상류와 하류단에 대해 각각 유입경계조건(inlet boundary condition)과 Flather 형식의 개방경계조건(open boundary condition)을 부여하여 일정 유량으로 흐르는 개수로 흐름을 구현하였다. 더불어, Strickler의 법칙을 확장한 반력공식을 연계하여, 개수로 흐름에 대한 조류 터빈의 영향을 반영하였다. 주어진 상류의 흐름 조건에 대해 조류발전량을 최대로 하는 최적 배열을 구하기 위해 터빈 반력모형을 연계한 천수방정식, 터빈간 최소간격, 그리고 발전단지영역을 제한조건으로 하는 발전량 최대화 문제를 구성하였다. 여기서 조류 터빈의 위치를 나타내는 벡터를 설계변수로 두었는데, 설계되는 터빈의 수가 증가함에 따라 최적화 문제의 계산량이 증가하지 않도록 수반법(adjoint method)을 경사도기반법(gradient-based method)에 연계한 방법이 이용되었다. 다수의 터빈초기배치로 상당한 수치실험이 수행되었고, 발전량 최대화를 이루도록 최적화된 터빈의 배치들이 큰 규모에서 고유한 형상으로 수렴함을 확인하였다. 이러한 특성은 발전단지의 너비와 터빈의 최소간격의 함수로 정의된 무차원수 E를 바탕으로 설명되었다. 구체적으로, E가 1보다 작을 때에는 선형배열이 최적배열로 나타났고, E가 1을 넘어 점차 커짐에 따라 하류에 오목한 형상을 보이다가 V-형태로 발전하는 양상을 보였다. 또한, 어느 임계 수 이상의 터빈이 배치되는 경우 일열 배열을 유지하지 못하고 이열 배열로 분리됨이 관찰되었다.

• Steel and Composite Structures
• /
• v.51 no.4
• /
• pp.363-375
• /
• 2024

Within the optimization field, addressing the intricate posed by fluidic pressure loads on functionally graded structures with frequency-related designs is a kind of complex design challenges. This paper thus introduces an innovative density-based topology optimization strategy for frequency-constraint functionally graded structures incorporating Darcy's law and a drainage term. It ensures consistent treatment of design-dependent fluidic pressure loads to frequency-related structures that dynamically adjust their direction and location throughout the design evolution. The porosity of each finite element, coupled with its drainage term, is intricately linked to its density variable through a Heaviside function, ensuring a seamless transition between solid and void phases. A design-specific pressure field is established by employing Darcy's law, and the associated partial differential equation is solved using finite element analysis. Subsequently, this pressure field is utilized to ascertain consistent nodal loads, enabling an efficient evaluation of load sensitivities through the adjoint-variable method. Moreover, this novel approach incorporates load-dependent structures, frequency constraints, functionally graded material models, and polygonal meshes, expanding its applicability and flexibility to a broader range of engineering scenarios. The proposed methodology's effectiveness and robustness are demonstrated through numerical examples, including fluidic pressure-loaded frequency-constraint structures undergoing small deformations, where compliance is minimized for structures optimized within specified resource constraints.

• Computational Structural Engineering
• /
• v.4 no.2
• /
• pp.87-94
• /
• 1991

In this paper, a mixed variational principle applicable to the linear elasticity of inhomogeneous anisotropic materials is presented. For derivation of the general variational principle, a systematic procedure for the variational formulation of linear coupled boundary value problems developed by Sandhu et al. is employed. Consistency condition of the field operators with the boundary operators results in explicit inclusion of boundary conditions in the governing functional. Extensions of admissible state function spaces and specialization to a certain relation in the general governing functional lead to the desired mixed variational principle. In the physical sense, the present variational principle is analogous to the Reissner's recent formulation obtained by applying Lagrange multiplier technique followed by partial Legendre transform to the classical minimum potential energy principle. However, the present one is more advantageous for the application to the general anisotropic materials since Reissner's principle contains an implicit function which is not easily converted to an explicit form.