• Title/Summary/Keyword: Spatial gradient

Search Result 330, Processing Time 0.025 seconds

High-Resolution Numerical Simulation of Respiration-Induced Dynamic B0 Shift in the Head in High-Field MRI

  • Lee, So-Hee;Barg, Ji-Seong;Yeo, Seok-Jin;Lee, Seung-Kyun
    • Investigative Magnetic Resonance Imaging
    • /
    • v.23 no.1
    • /
    • pp.38-45
    • /
    • 2019
  • Purpose: To demonstrate the high-resolution numerical simulation of the respiration-induced dynamic $B_0$ shift in the head using generalized susceptibility voxel convolution (gSVC). Materials and Methods: Previous dynamic $B_0$ simulation research has been limited to low-resolution numerical models due to the large computational demands of conventional Fourier-based $B_0$ calculation methods. Here, we show that a recently-proposed gSVC method can simulate dynamic $B_0$ maps from a realistic breathing human body model with high spatiotemporal resolution in a time-efficient manner. For a human body model, we used the Extended Cardiac And Torso (XCAT) phantom originally developed for computed tomography. The spatial resolution (voxel size) was kept isotropic and varied from 1 to 10 mm. We calculated $B_0$ maps in the brain of the model at 10 equally spaced points in a respiration cycle and analyzed the spatial gradients of each of them. The results were compared with experimental measurements in the literature. Results: The simulation predicted a maximum temporal variation of the $B_0$ shift in the brain of about 7 Hz at 7T. The magnitudes of the respiration-induced $B_0$ gradient in the x (right/left), y (anterior/posterior), and z (head/feet) directions determined by volumetric linear fitting, were < 0.01 Hz/cm, 0.18 Hz/cm, and 0.26 Hz/cm, respectively. These compared favorably with previous reports. We found that simulation voxel sizes greater than 5 mm can produce unreliable results. Conclusion: We have presented an efficient simulation framework for respiration-induced $B_0$ variation in the head. The method can be used to predict $B_0$ shifts with high spatiotemporal resolution under different breathing conditions and aid in the design of dynamic $B_0$ compensation strategies.

Calculating the Actual Surface Area for Gangneung Forest Fire Area Using Slope-Aspect Algorithm (Slope-Aspect 알고리즘을 활용한 강릉시 산불 피해지역 실표면적 산출 방법)

  • Jeong, JongChul
    • Journal of Cadastre & Land InformatiX
    • /
    • v.52 no.1
    • /
    • pp.95-104
    • /
    • 2022
  • This study aims to find the exact area of the forest fire in Okgye-myeon, Gangneung, April 4, 2019. Since there is a gradient in our country's forests, we should find a surface area that takes into account The 5th numerical clinical map provided by the DEM and the Korea Forest Service provided by the National Geographic Information Service was used. In DEM, the center point of each pixel was created and all points were connected. The length of the connecting line is determined by the spatial resolution of the pixel and the cosine value, and the surface area is obtained along with the height value, which is called the Slope-Aspect algorithm. The surface area and floor area of the forest were shown according to the tree species and types of forest, and their quantitative numerical differences proved the validity of this study.

Active Exchange of Water and Nutrients between Seawater and Shallow Pore Water in Intertidal Sandflats

  • Hwang, Dong-Woon;Kim, Gue-Buem;Yang, Han-Soeb
    • Ocean Science Journal
    • /
    • v.43 no.4
    • /
    • pp.223-232
    • /
    • 2008
  • In order to determine the temporal and spatial variations of nutrient profiles in the shallow pore water columns (upper 30 cm depth) of intertidal sandflats, we measured the salinity and nutrient concentrations in pore water and seawater at various coastal environments along the southern coast of Korea. In the intertidal zone, salinity and nutrient concentrations in pore water showed marked vertical changes with depth, owing to the active exchange between the pore water and overlying seawater, while they are temporally more stable and vertically constant in the sublittoral zone. In some cases, the advective flow of fresh groundwater caused strong vertical gradients of salinity and nutrients in the upper 10 cm depth of surface sediments, indicating the active mixing of the fresher groundwater with overlying seawater. Such upper pore water column profiles clearly signified the temporal fluctuation of lower-salinity and higher-Si seawater intrusion into pore water in an intertidal sandflat near the mouth of an estuary. We also observed a semimonthly fluctuation of pore water nutrients due to spring-neap tide associated recirculation of seawater through the upper sediments. Our study shows that the exchange of water and nutrients between shallow pore water and overlying seawater is most active in the upper 20 cm layer of intertidal sandflats, due to physical forces such as tides, wave set-up, and density-thermal gradient.

Numerical Study on the Vortex Evolution from a Sharp-Edged, Wall-Mounted Obstacle (장애물 주위의 와구조 형성과정에 관한 수치적 연구)

  • Hwang, Jong-Yeon;Yang, Kyung-Soo
    • Transactions of the Korean Society of Mechanical Engineers B
    • /
    • v.28 no.6
    • /
    • pp.672-681
    • /
    • 2004
  • Direct numerical simulation was carried out to study the vortical structures of the flow around a wall-mounted cube in a channel at Re=1,000 and Re=3,500 based on cubic height and bulk mean velocity. The cubic obstacle is situated in the entrance region of the channel flow where the boundary layers are developing. Upstream of the obstacle, steady and unsteady laminar horseshoe vortex systems are observed at Re=1,000 and Re=3,500, respectively; the near-wake flow is turbulent in both cases. The flow separates at each leading sharp edge of the cube, and subsequent vortex roll-up is noticed in the corresponding free-shear layer. The vortex shedding from the upper leading edge (upper vortices) and that from the two lateral leading edges (lateral vortices) are both quasi-periodic and their frequencies are computed. The upper and lateral vortices further develop into hairpin and Λ vortices, respectively. A series of instantaneous contours of the second invariant of velocity gradient tensor helps us identify spatial and temporal behaviors of the vortices in detail. The results indicate that the length and time scales of the vortical structures at Re=3,500 are much shorter than those at Re:1,000. Correlations between the upper and lateral vortices are also reported.

Global Optimization Using Kriging Metamodel and DE algorithm (크리깅 메타모델과 미분진화 알고리듬을 이용한 전역최적설계)

  • Lee, Chang-Jin;Jung, Jae-Jun;Lee, Kwang-Ki;Lee, Tae-Hee
    • Proceedings of the KSME Conference
    • /
    • 2001.06c
    • /
    • pp.537-542
    • /
    • 2001
  • In recent engineering, the designer has become more and more dependent on computer simulation. But defining exact model using computer simulation is too expensive and time consuming in the complicate systems. Thus, designers often use approximation models, which express the relation between design variables and response variables. These models are called metamodel. In this paper, we introduce one of the metamodel, named Kriging. This model employs an interpolation scheme and is developed in the fields of spatial statistics and geostatistics. This class of interpolating model has flexibility to model response data with multiple local extreme. By reason of this multi modality, we can't use any gradient-based optimization algorithm to find global extreme value of this model. Thus we have to introduce global optimization algorithm. To do this, we introduce DE(Differential Evolution). DE algorithm is developed by Ken Price and Rainer Storn, and it has recently proven to be an efficient method for optimizing real-valued multi-modal objective functions. This algorithm is similar to GA(Genetic Algorithm) in populating points, crossing over, and mutating. But it introduces vector concept in populating process. So it is very simple and easy to use. Finally, we show how we determine Kriging metamodel and find global extreme value through two mathematical examples.

  • PDF

Development of 3-D Flow Analysis Code Using Unstructured Grid System (I) - Numerical Method - (비정렬격자계를 사용하는 3차원 유동해석코드 개발 (I) - 수치해석방법 -)

  • Kim, Jong-Tae;Myong, Hyon-Kook
    • Transactions of the Korean Society of Mechanical Engineers B
    • /
    • v.29 no.9 s.240
    • /
    • pp.1049-1056
    • /
    • 2005
  • A conservative pressure-based finite-volume numerical method has been developed for computing flow and heat transfer by using an unstructured grid system. The method admits arbitrary convex polyhedra. Care is taken in the discretization and solution procedures to avoid formulations that are cell-shape-specific. A collocated variable arrangement formulation is developed, i.e. all dependent variables such as pressure and velocity are stored at cell centers. Gradients required for the evaluation of diffusion fluxes and for second-order-accurate convective operators are found by a novel second-order accurate spatial discretization. Momentum interpolation is used to prevent pressure checkerboarding and the SIMPLE algorithm is used for pressure-velocity coupling. The resulting set of coupled nonlinear algebraic equations is solved by employing a segregated approach, leading to a decoupled set of linear algebraic equations fer each dependent variable, with a sparse diagonally dominant coefficient matrix. These equations are solved by an iterative preconditioned conjugate gradient solver which retains the sparsity of the coefficient matrix, thus achieving a very efficient use of computer resources.

Numerical Wear Analysis of a Three-dimensional Rough Surface (수치적 방법을 이용한 3차원 거친 표면의 마모 해석)

  • Kim, Yunji;Suh, Junho;Kim, Bongjun;Yu, Yonghun
    • Tribology and Lubricants
    • /
    • v.36 no.4
    • /
    • pp.232-243
    • /
    • 2020
  • It is essential to predict the amount of wear and surface parameters for a surface where relative motion occurs. In the asperity-based model for wear prediction, only the average contact pressure can be obtained. Hence, the accuracy of wear analysis is poor. In this study, DC-FFT is used to obtain the pressure of each node, and wear analysis is performed by considering the effect of the pressure gradient. The numerical surface generation method is used to create Gaussian, negatively skewed, and positively skewed surfaces for wear analysis. The spatial and height distributions of each surface are analyzed to confirm the effectiveness of the generated surface. Furthermore, wear analysis is performed using DC-FFT and Archard's wear formula. After analysis, it is confirmed that all peaks are removed and only valleys remain on the surface. The RMS roughness and Sk continue to decrease and Ku increases as the cycle progresses. It is observed that the surface parameters are significantly affected by the radius of curvature of the asperity. This analysis method is more accurate than the existing average wear and truncation models because the change in asperity shape during the wear process is reflected in detail.

A Sensitivity Analysis on Numerical Grid Size of a Three-Dimensional Hydrodynamic and Water Quality Model (EFDC) for the Saemangeum Reservoir (새만금호 3차원 수리.수질모델(EFDC)의 수치격자 민감도 분석)

  • Jeon, Ji Hye;Chung, Se Woong
    • Journal of Korean Society on Water Environment
    • /
    • v.28 no.1
    • /
    • pp.26-37
    • /
    • 2012
  • Multi-dimensional hydrodynamic and water quality models are widely used to simulate the physical and biogeochemical processes in the surface water systems such as reservoirs and estuaries. Most of the models have adopted the Eulerian grid modeling framework, mainly because it can reasonably simulate physical dynamics and chemical species concentrations throughout the entire model domain. Determining the optimum grid cell size is important when using the Eulerian grid-based three-dimensional water quality models because the characteristics of species are assumed uniform in each of the grid cells and chemical species are represented by concentration (mass per volume). The objective of this study was to examine the effect of grid-size of a three dimensional hydrodynamic and water quality model (EFDC) on hydrodynamics and mass transport in the Saemangeum Reservoir. Three grid resolutions, respectively representing coarse (CG), medium (MG), and fine (FG) grid cell sizes, were used for a sensitivity analysis. The simulation results of numerical tracer showed that the grid resolution affects on the flow path, mass transport, and mixing zone of upstream inflow, and results in a bias of temporal and spatial distribution of the tracer. With the CG, in particular, the model overestimates diffusion in the mixing zone, and fails to identify the gradient of concentrations between the inflow and the ambient water.

Estimation of the air temperature over the sea using the satellite data

  • Kwon B. H.;Hong G. M.;Kim Y. S.
    • Proceedings of the KSRS Conference
    • /
    • 2005.10a
    • /
    • pp.392-393
    • /
    • 2005
  • Due to the temporal and spatial simultaneity and the high-frequency repetition, the data set retrieved from the satellite observation is considered to be the most desirable ones for the study of air-sea interaction. With rapidly developing sensor technology, satellite-retrieved data has experienced improvement in the accuracy and the number of parameters. Nevertheless, since it is still impossible to directly measure the heat fluxes between air and sea, the bulk method is an exclusive way for the evaluation of the heat fluxes at the sea surface. It was noted that the large deviation of air temperature in the winter season by the linear regression despite good correlation coefficients. We propose a new algorithm based on the Fourier series with which the SST and the air temperature. We found that the mean of air temperature is a function of the mean of SST with the monthly gradient of SST inferred from the latitudinal variation of SST and the spectral energy of air temperature is related linearly to that of SST. An algorithm to obtain the air temperature over the sea was completed with a proper analysis on the relation between of air temperature and of SST. This algorithm was examined by buoy data and therefore the air temperature over the sea can be retrieved based on just satellite data.

  • PDF

Performance Test of Asynchronous Process of OGC WPS 2.0: A Case Study for Geo-based Image Processing

  • Yoon, Gooseon;Lee, Kiwon
    • Korean Journal of Remote Sensing
    • /
    • v.33 no.4
    • /
    • pp.391-400
    • /
    • 2017
  • Geo-based application services linked with the Open Geospatial Consortium (OGC) Web Processing Service (WPS) protocol have been regarded as an important standardized framework for of digital earth building in the web environments. The WPS protocol provides interface standards for analysis functionalities within geo-spatial processing in web-based service systems. Despite its significance, there is few performance tests of WPS applications. The main motivation in this study is to perform the comparative performance test on WPS standards. Test system, which was composed of WPS servers, WPS framework, data management module, geo-based data processing module and client-sided system, was implemented by fully open source stack. In this system, two kinds of geo-based image processing functions such as cloud detection and gradient magnitude computation were applied. The performance test of different server environments of non-WPS, synchronous WPS 1.0 and asynchronous WPS 2.0 was carried out using 100 threads and 400 threads corresponds client users on a web-based application service. As the result, at 100 threads, performance of three environments was within an adjacent range in the average response time to complete the processing of each thread. At 400 threads, the application case of WPS 2.0 showed the distinguished characteristics for higher performance in the response time than the small threads cases. It is thought that WPS 2.0 contributes to settlement of without performance problems such as time delay or thread accumulation.