• Nuclear Engineering and Technology
• /
• 제48권3호
• /
• pp.765-772
• /
• 2016

In analyzing residual radiation, researchers generally use a two-step Monte Carlo (MC) simulation. The first step (MC1) simulates neutron transport, and the second step (MC2) transports the decay photons emitted from the activated materials. In this process, the stochastic uncertainty estimated by the MC2 appears only as a final result, but it is underestimated because the stochastic error generated in MC1 cannot be directly included in MC2. Hence, estimating the true stochastic uncertainty requires quantifying the propagation degree of the stochastic error in MC1. The brute force technique is a straightforward method to estimate the true uncertainty. However, it is a costly method to obtain reliable results. Another method, called the adjoint-based method, can reduce the computational time needed to evaluate the true uncertainty; however, there are limitations. To address those limitations, we propose a new strategy to estimate uncertainty propagation without any additional calculations in two-step MC simulations. To verify the proposed method, we applied it to activation benchmark problems and compared the results with those of previous methods. The results show that the proposed method increases the applicability and user-friendliness preserving accuracy in quantifying uncertainty propagation. We expect that the proposed strategy will contribute to efficient and accurate two-step MC calculations.

• 대한기계학회논문집
• /
• 제17권1호
• /
• pp.141-152
• /
• 1993

본 연구에서는 축대칭 열전도 고체의 형상설계민감도 해석을 위하여 2차원 문 제를 다룬 Lee, Choi와 Kwak의 방법을 축대칭 문제로 확장하였다.축대칭 형태로 표 시된 직접 및 간접 경계적분방정식의 정식화에 기초하여 전미분방접과 보조변수방법으 로 형상최적화 문제에서 발생하는 일반적인 성능 범함수의 형상설계민감도 공식을 유 도하고, 온도 및 열속의 제한조건에 이를 응용하였다. 제시된 민감도해석방법의 정 확성을 검증하기 위하여 해석적인 해를 갖는 원통문제와 구문제를 다루었는데, 두 문 제에 대하여 민감도 공식을 이용하여 수치계산된 결과를 해석적인 민감도와 비교하였 다. 또한 복잡한 수치해로서 냉각핀(cooling fin)문제를 다루었으며, 민감도 공식에 의한 계산 결과를 유한차분(finite difference)으로 수치미분한 결과와 비교하였다.

• 대한기계학회논문집B
• /
• 제33권1호
• /
• pp.17-27
• /
• 2009

Water jet impingement cooling is used to remove heat from high-temperature surfaces such as hot steel plates in the steel manufacturing process (thermo-mechanical cooling process; TMCP). In those processes, uniform cooling is the most critical factor to ensure high strength steel and good quality. In this study, experiments are performed to measure the heat transfer coefficient together with the inverse heat conduction problem (IHCP) analysis for a plate cooled by planar water jet. In the inverse heat transfer analysis, spatial and temporal variations of heat transfer coefficient, with no information regarding its functional form, are determined by employing the conjugate gradient method with an adjoint problem. To estimate the two dimensional distribution of heat transfer coefficient and heat flux for planar waterjet cooling, eight thermo-couple are installed inside the plate. The results show that heat transfer coefficient is approximately uniform in the span-wise direction in the early stage of cooling. In the later stage where the forced-convection effect is important, the heat transfer coefficient becomes larger in the edge region. The surface temperature vs. heat flux characteristics are also investigated for the entire boiling regimes. In addition, the heat transfer rate for the two different plate geometries are compared at the same Reynolds number.

• Nuclear Engineering and Technology
• /
• 제49권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.

• Steel and Composite Structures
• /
• 제51권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.

• Nuclear Engineering and Technology
• /
• 제56권7호
• /
• pp.2748-2755
• /
• 2024

This paper presents two new methods for variance reduction for shielding calculation in Monte Carlo radiation transport. One method is CADIS-NEE, which combines Consistent Adjoint Driven Importance Sampling (CADIS) and next-event estimator (NEE) methods to increase the calculation efficiency of tallies at points. The other is CADIS-deterministic transport (DXTRAN), which combines CADIS and DXTRAN to obtain higher performance than using CADIS and DXTRAN separately. The combination processes are derived and implemented in the hybrid Monte-Carlo-Deterministic particle-transport code NECP-MCX. Various problems are tested to demonstrate the effectiveness of the two methods. According to the results, the two combination methods have higher efficiency than using CADIS, NEE or DXTRAN separately. In a long-distance photon-transport problem, CADIS-NEE converges faster than NEE and the figure of merit (FOM) of CADIS-NEE is 75.6 times of NEE. In a labyrinthine problem, CADIS-DXTRAN's FOM surpasses that of DXTRAN and CADIS by a factor of 45.3 and 17.7, respectively. Therefore, it is advisable to employ these two novel methods selectively in appropriate scenarios to reduce variance.

• 한국전산구조공학회논문집
• /
• 제25권6호
• /
• pp.549-558
• /
• 2012

본 논문에서는 레벨셋 방법을 이용하여, 소음을 차단하기 위한 음향 구조물의 형상 최적설계를 수행하였다. 형상 최적설계의 목적은 특정한 각도와 각속도로 입사되는 입사파에 대해서 음향 투과율(acoustic transmittance)이 최소가 되도록 음향 결정의 형상(inclusion shape)을 결정하는 것이다. 음향 결정 구조에서는 음향이 흩어져 있는 결정 구조에 의해서 굴절되기 때문에 결정 모양을 조정함으로써, 음향 거동을 제어할 수 있다. 본 연구에서는 음향 구조물로 결정이 수평방향으로는 주기적으로 무한히 분포하고 수직방향으로는 유한한 층간 구조를 가지고 있는 소음 방어벽(Noise barrier)을 고려한다. 주기적 구조물을 고려하기 때문에 결정의 좌와 우에 Bloch 이론을 적용해 주기적 경계조건을 부과하였고, 소음 방어벽 위와 아래에는 임피던스 행렬(impedance matrix)를 이용하여, 무한 균질 영역과 소음 방어벽 사이의 음파 투과를 모사하였다. 결정의 위상과 형상변화를 묘사하기 위해서 레벨셋 방법(level set method)을 사용하였다. 레벨셋 방법에서는 초기 영역을 고정시킨 상태에서, 레벨셋으로 표현되는 임시적 경계(implicit moving boundary)를 변화시킴으로써 복잡한 형상을 다룰 수 있다. 몇몇 수치적 예제를 통해, 제시된 방법의 적용성을 검증하였다.

• 한국전산구조공학회논문집
• /
• 제27권1호
• /
• pp.9-16
• /
• 2014

레벨셋 기법과 위상민감도를 이용하여 선형 탄성 구조물에 대하여, 초기 설계형상에 의존성이 없는 위상 및 형상 최적설계 기법을 개발하였다. 레벨셋 기법에서는 복잡한 위상 형상변화를 쉽게 다루기 위해 초기 영역은 고정한 채 레벨셋 함수로 표현되는 암시적 이동경계로 경계를 표현한다. 해밀턴-자코비(H-J) 방정식과 수치적으로 강건한 기법인 'up-wind scheme'은 컴플라이언스 목적함수를 최소화시키고 허용체적 제약조건을 만족시키면서, 초기 암시적 경계를 법선 속도장에 따라 최적의 형상으로 이끌어 낸다. 점근적인 정규화 개념에 근거하여, 구멍의 반지름을 0으로 접근시켜 형상 미분의 극한을 취한 위상민감도를 고려하였다. 최적조건으로부터 유도된 라그란지안의 감소 방향을 이용하여 H-J 방정식을 갱신하기 위한 속도장을 결정하였다. 개발한 방법에서는 위상민감도로부터 얻어지는 지표를 이용하여 구멍을 언제든지 어디에서나 생성가능하기 때문에 초기 구멍이 최적 형상을 얻기 위해 요구되지 않는다는 사실을 확인하였다. 또한 효율적인 최적화 과정을 위해서는 구멍 생성을 위한 조정변수의 적절한 선택이 중요함을 확인하였다.

• 한국전자통신학회논문지
• /
• 제12권1호
• /
• pp.189-194
• /
• 2017

양자논리회로의 가역성은 정보 가역적 및 에너지 가역적 회로라는 두 가지 가역 조건을 만족할 때 실현될 수 있다. 본 논문은 다치 양자논리 회로에서 원래상태로의 정보가역성 회복에 필요한 연산 사이클을 모델링하였다. 모델링을 위해 유니터리 스위치를 산술 멱승 스위치로 사용하는 함수 임베딩 방법을 사용하였다. 양자논리회로에서 수반게이트 쌍이 대칭이면 유니터리 스위치함수가 균형함수 특성을 보임으로써 원래상태의 정보 가역성 회복에 1 사이클 연산이 소요되었다. 반대로 비대칭 구조이면 상수 함수에 의해 2 사이클 연산이 소요되었다. 본 논문은 ternary M-S 게이트로 hybrid MCT 게이트를 실현할 경우의 비대칭 구조에 따른 2 사이클 복원 문제는 비대칭 구조의 수반게이트들을 대칭구조의 수반게이트로 등가 변환하여 해결할 수 있음을 밝혔다.

• 대기
• /
• 제27권1호
• /
• pp.93-104
• /
• 2017

The effect of Advanced Microwave Sounding Unit-A (AMSU-A) observations on the short-range forecast in East Asia (EA) was investigated for the Northern Hemispheric (NH) summer and winter months, using the Forecast Sensitivity to Observations (FSO) method. For both periods, the contribution of radiosonde (TEMP) to the EA forecast was largest, followed by AIRCRAFT, AMSU-A, Infrared Atmospheric Sounding Interferometer (IASI), and the atmospheric motion vector of Communication, Ocean and Meteorological Satellite (COMS) or Multi-functional Transport Satellite (MTSAT). The contribution of AMSU-A sensor was largely originated from the NOAA 19, NOAA 18, and MetOp-A (NOAA 19 and 18) satellites in the NH summer (winter). The contribution of AMSU-A sensor on the MetOp-A (NOAA 18 and 19) satellites was large at 00 and 12 UTC (06 and 18 UTC) analysis times, which was associated with the scanning track of four satellites. The MetOp-A provided the radiance data over the Korea Peninsula in the morning (08:00~11:30 LST), which was important to the morning forecast. In the NH summer, the channel 5 observations on MetOp-A, NOAA 18, 19 along the seaside (along the ridge of the subtropical high) increased (decreased) the forecast error slightly (largely). In the NH winter, the channel 8 observations on NOAA 18 (NOAA 15 and MetOp-A) over the Eastern China (Tibetan Plateau) decreased (increased) the forecast error. The FSO provides useful information on the effect of each AMSU-A sensor on the EA forecasts, which leads guidance to better use of AMSU-A observations for EA regional numerical weather prediction.