• Title/Summary/Keyword: iterative method

Search Result 2,058, Processing Time 0.036 seconds

Precise System Models using Crystal Penetration Error Compensation for Iterative Image Reconstruction of Preclinical Quad-Head PET

  • Lee, Sooyoung;Bae, Seungbin;Lee, Hakjae;Kim, Kwangdon;Lee, Kisung;Kim, Kyeong-Min;Bae, Jaekeon
    • Journal of the Korean Physical Society
    • /
    • v.73 no.11
    • /
    • pp.1764-1773
    • /
    • 2018
  • A-PET is a quad-head PET scanner developed for use in small-animal imaging. The dimensions of its volumetric field of view (FOV) are $46.1{\times}46.1{\times}46.1mm^3$ and the gap between the detector modules has been minimized in order to provide a highly sensitive system. However, such a small FOV together with the quad-head geometry causes image quality degradation. The main factor related to image degradation for the quad-head PET is the mispositioning of events caused by the penetration effect in the detector. In this paper, we propose a precise method for modelling the system at the high spatial resolution of the A-PET using a LOR (line of response) based ML-EM (maximum likelihood expectation maximization) that allows for penetration effects. The proposed system model provides the detection probability of every possible ray-path via crystal sampling methods. For the ray-path sampling, the sub-LORs are defined by connecting the sampling points of the crystal pair. We incorporate the detection probability of each sub-LOR into the model by calculating the penetration effect. For comparison, we used a standard LOR-based model and a Monte Carlo-based modeling approach, and evaluated the reconstructed images using both the National Electrical Manufacturers Association NU 4-2008 standards and the Geant4 Application for Tomographic Emission simulation toolkit (GATE). An average full width at half maximum (FWHM) at different locations of 1.77 mm and 1.79 mm are obtained using the proposed system model and standard LOR system model, which does not include penetration effects, respectively. The standard deviation of the uniform region in the NEMA image quality phantom is 2.14% for the proposed method and 14.3% for the LOR system model, indicating that the proposed model out-performs the standard LOR-based model.

Application of POD reduced-order algorithm on data-driven modeling of rod bundle

  • Kang, Huilun;Tian, Zhaofei;Chen, Guangliang;Li, Lei;Wang, Tianyu
    • Nuclear Engineering and Technology
    • /
    • v.54 no.1
    • /
    • pp.36-48
    • /
    • 2022
  • As a valid numerical method to obtain a high-resolution result of a flow field, computational fluid dynamics (CFD) have been widely used to study coolant flow and heat transfer characteristics in fuel rod bundles. However, the time-consuming, iterative calculation of Navier-Stokes equations makes CFD unsuitable for the scenarios that require efficient simulation such as sensitivity analysis and uncertainty quantification. To solve this problem, a reduced-order model (ROM) based on proper orthogonal decomposition (POD) and machine learning (ML) is proposed to simulate the flow field efficiently. Firstly, a validated CFD model to output the flow field data set of the rod bundle is established. Secondly, based on the POD method, the modes and corresponding coefficients of the flow field were extracted. Then, an deep feed-forward neural network, due to its efficiency in approximating arbitrary functions and its ability to handle high-dimensional and strong nonlinear problems, is selected to build a model that maps the non-linear relationship between the mode coefficients and the boundary conditions. A trained surrogate model for modes coefficients prediction is obtained after a certain number of training iterations. Finally, the flow field is reconstructed by combining the product of the POD basis and coefficients. Based on the test dataset, an evaluation of the ROM is carried out. The evaluation results show that the proposed POD-ROM accurately describe the flow status of the fluid field in rod bundles with high resolution in only a few milliseconds.

An improved time-domain approach for the spectra-compatible seismic motion generation considering intrinsic non-stationary features

  • Feng Cheng;Jianbo Li;Zhixin Ding;Gao Lin
    • Nuclear Engineering and Technology
    • /
    • v.55 no.3
    • /
    • pp.968-980
    • /
    • 2023
  • The dynamic structural responses are sensitive to the time-frequency content of seismic waves, and seismic input motions in time-history analysis are usually required to be compatible with design response spectra according to nuclear codes. In order to generate spectra-compatible input motions while maintaining the intrinsic non-stationarity of seismic waves, an improved time-domain approach is proposed in this paper. To maintain the nonstationary characteristics of the given seismic waves, a new time-frequency envelope function is constructed using the Hilbert amplitude spectrum. Based on the intrinsic mode functions (IMFs) obtained from given seismic waves through variational mode decomposition, a new corrective time history is constructed to locally modify the given seismic waves. The proposed corrective time history and time-frequency envelope function are unique for each earthquake records as they are extracted from the given seismic waves. In addition, a dimension reduction iterative technique is presented herein to simultaneously superimpose corrective time histories of all the damping ratios at a specific frequency in the time domain according to optimal weights, which are found by the genetic algorithm (GA). Examples are presented to show the capability of the proposed approach in generating spectra-compatible time histories, especially in maintaining the nonstationary characteristics of seismic records. And numerical results reveal that the modified time histories generated by the proposed method can obtain similar dynamic behaviors of AP1000 nuclear power plant with the natural seismic records. Thus, the proposed method can be efficiently used in the design practices.

Analytical study on cable shape and its lateral and vertical sags for earth-anchored suspension bridges with spatial cables

  • Gen-min Tian;Wen-ming Zhang;Jia-qi Chang;Zhao Liu
    • Structural Engineering and Mechanics
    • /
    • v.87 no.3
    • /
    • pp.255-272
    • /
    • 2023
  • Spatial cable systems can provide more transverse stiffness and torsional stiffness without sacrificing the vertical bearing capacity compared with conventional vertical cable systems, which is quite lucrative for long-span earth-anchored suspension bridges' development. Higher economy highlights the importance of refined form-finding analysis. Meanwhile, the internal connection between the lateral and vertical sags has not yet been specified. Given this, an analytic algorithm of form-finding for the earth-anchored suspension bridge with spatial cables is proposed in this paper. Through the geometric compatibility condition and mechanical equilibrium condition, the expressions for cable segment, the recurrence relationship between catenary parameters and control equations of spatial cable are established. Additionally, the nonlinear general reduced gradient method is introduced into fast and high-precision numerical analysis. Furthermore, the analytic expression of the lateral and vertical sags is deduced and discussed. This is very significant for the space design above the bridge deck and the optimization of the sag-to-span ratio in the preliminary design stage of the bridge. Finally, the proposed method is verified with the aid of two examples, one being an operational self-anchored suspension bridge (with spatial cables and a 260 m main span), and the other being an earth-anchored suspension bridge under design (with spatial cables and a 500 m main span). The necessity of an iterative calculation for hanger tensions on earth-anchored suspension bridges is confirmed. It is further concluded that the main cable and their connected hangers are in very close inclined planes.

The Study of Reducing Radiation Exposure Dose and Comparing SUV According to Applied IRIS (Iterative Reconstruction in Image Space) for PET/CT (PET/CT 검사 시 IRIS (Iterative Reconstruction in Image Space) 적용에 따른 CT 피폭선량 감소와 PET SUV 비교 연구)

  • Do, Yong Ho;Song, Ho Jun;Lee, Hyung Jin;Lee, Hong Jae;Kim, Jin Eui
    • The Korean Journal of Nuclear Medicine Technology
    • /
    • v.16 no.2
    • /
    • pp.29-34
    • /
    • 2012
  • Purpose : Presently, hardwares and softwares for reducing radiation exposure are continually developed for PET/CT examination. Purpose of this study is to evaluate effectiveness of reducing radiation exposure dose of CT and SUV changes of PET when applied each kernel to ACCT (Attenuation Correction Computed Tomography) according to adopted IRIS (Iterative Reconstruction in Image Space) software. Materials and Methods : Biograph mCT (Siemens, Germany) was used as a PET/CT scanner. Using AAPM CT performance phantom, from standard (120 kVp, 100 mAs), 7 scans were conducted by reducing 15 mAs each. After image reconstruction by FBP (Filtered Back Projection) and IRIS, noise and spatial resolution were evaluated. The same method was applied to anthropomorphic chest phantom and acquired images were compared. NEMA IEC body phantom was used for SUV evaluation. Injected dose rate for hot sphere (hot) and background cylinder (BKG) were 1:8. CT dose condition (120 kVp, 50 mAs) was the same for each scan and PET scan durations were 1, 2, 3 and 4min. After scanning, each kernel of IRIS was applied to ACCT. And PET images were reconstructed by ACCT adopted IRIS for comparing SUV changes. Results : AAPM phantom test for noise evaluation, SD for FBP 100 mAs, IRIS 55 mAs were 8.8 and 8.9. FBP 85 mAs, IRIS 40 mAs were 9.5 and 9.7. FBP 70 mAs, IRIS 25 mAs were 11.9 and 11.1. Above mAs condition for FBP and IRIS, SD showed similar values. And for spatial resolution test, there was no significant difference. For chest phantom test, when applied the same mAs and kernel to both of FBP and IRIS, every applied kernels showed reduced noise. Lower mAs and higher kernel value showed higher noise reduction. There was no considerable difference only except for I70 very sharp kernel for SUV comparison using NEMA IEC body phantom. Conclusion : In this study, low mAs (55 mAs) applied IRIS and standard mAs (100 mAs) applied FBP showed similar noise. And only except for I70 kernel, there was no significant SUV changes. It is possible to reduce needless radiation exposure and acquire better image quality than FBP's through applying appropriate kernel of IRIS to PET/CT.

  • PDF

Reducing Dose in SPECT/CT Using Adaptive Statistical Iterative Reconstruction Technique (Adaptive Statistical Iterative Reconstruction 기법을 이용한 Bone SPECT/CT 검사에서 피폭량 감소 방안)

  • Choi, Jin-Wook;Choi, Hyeon-Jun;Park, Chan-Rok;Cho, Sung-Wook;Kim, Jin-Eui;Lee, Jae-Sung;Lee, Dong-Soo
    • The Korean Journal of Nuclear Medicine Technology
    • /
    • v.18 no.1
    • /
    • pp.134-139
    • /
    • 2014
  • Purpose: Adaptive statistical iterative reconstruction (ASIR) technique is a reconstruction method of CT image using statistical noise modeling which is known to reduce image noise and to preserve image quality despite reducing radiation dose. The aim of this study is to evaluate images using ASIR on bone SPECT/CT which is primarily performed in our hospital. Materials and Methods: We compared the images of applied ASIR (ASIR level: 20-80%) and none ASIR by changing the mA based on 120 kVp, 100 mA using Discovery NM/CT 670 (GE, U.S.A). First, we evaluated attenuation correction in SPECT image by changing the ASIR level using Anthropomorphic phantom. Second, we compared the contrast to noise ratio (CNR), image noise and spatial resolution in CT image using ACR phantom. Third, after selecting the ASIR level applicable patient using lower torso phantom, we examined 2 patients who followed up bone SPECT/CT and we performed blind test. Results: The degree of attenuation correction in SPECT image showed no significant difference between applied ASIR and none ASIR (P>0.05). When applied ASIR, the noise of CT image were reduced at least 17 up to 52% by changing the mA. The CNR of image with ASIR was maintained more than 0.8 at 40 mA (ASIR 60%) while those without ASIR showed 0.42 at standard 40 mA. In comparison of the high contrast object, we distinguished 12 line pairs/cm at 40 mA regardless of appling ASIR. Comparison of the patients image applied ASIR level 60% (40 mA) which found out by spine image of lower torso phantom showed no signigicant difference between applied ASIR and none ASIR in blind test. The CTDIvol and DLP for applied ASIR 60% showed decreased by 60%, 60% on average than using standard mA. Conclusion: The study show that the radiation dose in SPECT/CT using ASIR can be reduced despite degradation of SPECT and CT images. In addition, higher ASIR level could be possibly applied characteristics of SPECT/CT that region of interest is limited to bone.

  • PDF

FEM-based Seismic Reliability Analysis of Real Structural Systems (실제 구조계의 유한요소법에 기초한 지진 신뢰성해석)

  • Huh Jung-Won;Haldar Achintya
    • Journal of the Computational Structural Engineering Institute of Korea
    • /
    • v.19 no.2 s.72
    • /
    • pp.171-185
    • /
    • 2006
  • A sophisticated reliability analysis method is proposed to evaluate the reliability of real nonlinear complicated dynamic structural systems excited by short duration dynamic loadings like earthquake motions by intelligently integrating the response surface method, the finite element method, the first-order reliability method, and the iterative linear interpolation scheme. The method explicitly considers all major sources of nonlinearity and uncertainty in the load and resistance-related random variables. The unique feature of the technique is that the seismic loading is applied in the time domain, providing an alternative to the classical random vibration approach. The four-parameter Richard model is used to represent the flexibility of connections of real steel frames. Uncertainties in the Richard parameters are also incorporated in the algorithm. The laterally flexible steel frame is then reinforced with reinforced concrete shear walls. The stiffness degradation of shear walls after cracking is also considered. The applicability of the method to estimate the reliability of real structures is demonstrated by considering three examples; a laterally flexible steel frame with fully restrained connections, the same steel frame with partially restrained connections with different rigidities, and a steel frame reinforced with concrete shear walls.

Seismic Traveltime Tomography using Neural Network (신경망 이론을 이용한 탄성파 주시 토모그래피의 연구)

  • Kim, Tae-Yeon;Yoon, Wang-Jung
    • Geophysics and Geophysical Exploration
    • /
    • v.2 no.4
    • /
    • pp.167-173
    • /
    • 1999
  • Since the resolution of the 2-D hole-to-hole seismic traveltime tomography is affected by the limited ray transmission angle, various methods were used to improve the resolution. Linear traveltime interpolation(LTI) ray tracing method was chosen for forward-modeling method. Inversion results using the LTI method were compared with those using the other ray tracing methods. As an inversion algorithm, SIRT method was used. In the iterative non-linear inversion method, the cost of ray tracing is quite expensive. To reduce the cost, each raypath was stored and the inversion was performed from this information. Using the proposed method, fast convergence was achieved. Inversion results are likely to be affected by the initial velocity guess, especially when the ray transmission angle was limited. To provide a good initial guess for the inversion, generalized regression neural network(GRNN) method was used. When the transmitted raypath angle is not limited or the geological model is very complex, the inversion results are not affected by initial velocity model very much. Since the raypath angles, however, are limited in most geophysical tomographic problems, the enhancement of resolution in tomography can be achieved by providing a proper initial velocity model by another inversion algorithm such as GRNN.

  • PDF

Direct Inelastic Strut-Tie Model Using Secant Stiffness (할선강성을 이용한 직접 비탄성 스트럿-타이 모델)

  • Park Hong-Gun;Kim Yun-Gon;Eom Tae-Sung
    • Journal of the Korea Concrete Institute
    • /
    • v.17 no.2 s.86
    • /
    • pp.201-212
    • /
    • 2005
  • A new strut-tie model using secant stiffness, Direct Inelastic Strut-Tie Model, was developed. Since basically the proposed design method uses linear analysis, it is convenient and stable in numerical analysis. At the same time, the proposed design method can accurately estimate the inelastic strength and ductility demands of struts and ties because it can analyzes the inelastic behavior of structure using iterative calculations for secant stiffness. In the present study, the procedure of the proposed design method was established, and a computer program incorporating the proposed method was developed. Design examples using the proposed method were presented, and its advantages were highlighted by the comparison with the traditional strut-tie model. The Direct Inelastic Strut-Tie Model, as an integrated analysis/design method, can directly address the design strategy intended by the engineer to prevent development of macro-cracks and brittle failure of struts. Since the proposed model can analyze the inelastic deformation, indeterminate strut-tie model can be used. Also, since the proposed model controls the local deformations of struts and ties, it can be used as a performance-based design method for various design criteria.

Ship Collision Risk of Suspension Bridge and Design Vessel Load (현수교의 선박충돌 위험 및 설계박하중)

  • Lee, Seong Lo;Bae, Yong Gwi
    • KSCE Journal of Civil and Environmental Engineering Research
    • /
    • v.26 no.1A
    • /
    • pp.11-19
    • /
    • 2006
  • In this study ship collision risk analysis is performed to determine the design vessel for collision impact analysis of suspension bridge. Method II in AASHTO LRFD bridge design specifications which is a more complicated probability based analysis procedure is used to select the design vessel for collision impact. From the assessment of ship collision risk for each bridge pier exposed to ship collision, the design impact lateral strength of bridge pier is determined. The analysis procedure is an iterative process in which a trial impact resistance is selected for a bridge component and a computed annual frequency of collapse(AF) is compared to the acceptance criterion, and revisions to the analysis variables are made as necessary to achieve compliance. The acceptance criterion is allocated to each pier using allocation weights based on the previous predictions. This AF allocation method is compared to the pylon concentration allocation method to obtain safety and economy in results. This method seems to be more reasonable than the pylon concentration allocation method because AF allocation by weights takes the design parameter characteristics quantitatively into consideration although the pylon concentration allocation method brings more economical results when the overestimated design collision strength of piers compared to the strength of pylon is moderately modified. The design vessel for each pier corresponding with the design impact lateral strength obtained from the ship collision risk assessment is then selected. The design impact lateral strength can vary greatly among the components of the same bridge, depending upon the waterway geometry, available water depth, bridge geometry, and vessel traffic characteristics. Therefore more researches on the allocation model of AF and the selection of design vessel are required.