• Title/Summary/Keyword: Numerical simulations

검색결과 4,005건 처리시간 0.031초

NUMERICAL SIMULATIONS FOR THE CONTRACTION FLOW USING GRID GENERATION

  • Salem, S.A.
    • Journal of applied mathematics & informatics
    • /
    • 제16권1_2호
    • /
    • pp.383-405
    • /
    • 2004
  • We study the incomprssible Navier Stokes equations for the flow inside contraction geometry. The governing equations are expressed in the vorticity-stream function formulations. A rectangular computational domain is arised by elliptic grid generation technique. The numerical solution is based on a technique of automatic numerical generation of acurvilinear coordinate system by transforming the governing equation into computational plane. The transformed equations are approximated using central differences and solved simultaneously by successive over relaxation iteration. The time dependent of the vorticity equation solved by using explicit marching procedure. We will apply the technique on several irregular-shapes.

COMPARISON OF NUMERICAL METHODS FOR OPTION PRICING UNDER THE CGMY MODEL

  • Lee, Ahram;Lee, Younhee
    • 충청수학회지
    • /
    • 제29권3호
    • /
    • pp.503-508
    • /
    • 2016
  • We propose a number of finite difference methods for the prices of a European option under the CGMY model. These numerical methods to solve a partial integro-differential equation (PIDE) are based on three time levels in order to avoid fixed point iterations arising from an integral operator. Numerical simulations are carried out to compare these methods with each other for pricing the European option under the CGMY model.

경사식 방파제의 전면에 설치된 수중방파제의 영향에 관한 연구 (Influence of Submerged Breakwater in front of Rubble Mound Breakwater)

  • 민현성;조용식
    • 한국방재학회:학술대회논문집
    • /
    • 한국방재학회 2008년도 정기총회 및 학술발표대회
    • /
    • pp.217-220
    • /
    • 2008
  • The reflection coefficients and the run-up heights affected by submerged structures are studied by using the numerical and the laboratory experimental methods. The three-point method is chosen to calculate the reflection coefficients in both the experimental and the numerical methods. The results of numerical simulations are shown a good agreement with laboratory measurements. The reflection coefficients increase and the run-up heights decrease when the rubble mound breakwater is defended by low-crested structures.

  • PDF

Development of a Plasma Training Lab kart: System Setup and Numerical Simulation

  • Joo, Junghoon
    • Applied Science and Convergence Technology
    • /
    • 제26권6호
    • /
    • pp.195-200
    • /
    • 2017
  • A mobile lab kart for plasma training is developed with a high vacuum pumping system, vacuum gauges and a glass discharge tube powered by a high voltage transformer connected to a household 60 Hz line. A numerical model is developed by using a commercial multiphysics software package, CFD-ACE+ to analyze the experimental data. Simulations for argon and nitrogen were carried out to provide fundamental discharge characteristics. Variations of the kart configuration were demonstrated: a glass tube with three electric probes, optical emission spectrometer attachment and infra red thermal imaging system to give more detailed analysis of the discharge characteristics.

비특이화 간접경계적분방정식방법을 이용한 2차원 수치수조 개발 및 적용 (Development and Application of Two-Dimensional Numerical Tank using Desingularized Indirect Boundary Integral Equation Method)

  • 오승훈;조석규;정동호;성홍근
    • 한국해양공학회지
    • /
    • 제32권6호
    • /
    • pp.447-457
    • /
    • 2018
  • In this study, a two-dimensional fully nonlinear transient wave numerical tank was developed using a desingularized indirect boundary integral equation method. The desingularized indirect boundary integral equation method is simpler and faster than the conventional boundary element method because special treatment is not required to compute the boundary integral. Numerical simulations were carried out in the time domain using the fourth order Runge-Kutta method. A mixed Eulerian-Lagrangian approach was adapted to reconstruct the free surface at each time step. A numerical damping zone was used to minimize the reflective wave in the downstream region. The interpolating method of a Gaussian radial basis function-type artificial neural network was used to calculate the gradient of the free surface elevation without element connectivity. The desingularized indirect boundary integral equation using an isolated point source and radial basis function has no need for information about the element connectivity and is a meshless method that is numerically more flexible. In order to validate the accuracy of the numerical wave tank based on the desingularized indirect boundary integral equation method and meshless technique, several numerical simulations were carried out. First, a comparison with numerical results according to the type of desingularized source was carried out and confirmed that continuous line sources can be replaced by simply isolated sources. In addition, a propagation simulation of a $2^{nd}$-order Stokes wave was carried out and compared with an analytical solution. Finally, simulations of propagating waves in shallow water and propagating waves over a submerged bar were also carried and compared with published data.

Numerical discrepancy between serial and MPI parallel computations

  • Lee, Sang Bong
    • International Journal of Naval Architecture and Ocean Engineering
    • /
    • 제8권5호
    • /
    • pp.434-441
    • /
    • 2016
  • Numerical simulations of 1D Burgers equation and 2D sloshing problem were carried out to study numerical discrepancy between serial and parallel computations. The numerical domain was decomposed into 2 and 4 subdomains for parallel computations with message passing interface. The numerical solution of Burgers equation disclosed that fully explicit boundary conditions used on subdomains of parallel computation was responsible for the numerical discrepancy of transient solution between serial and parallel computations. Two dimensional sloshing problems in a rectangular domain were solved using OpenFOAM. After a lapse of initial transient time sloshing patterns of water were significantly different in serial and parallel computations although the same numerical conditions were given. Based on the histograms of pressure measured at two points near the wall the statistical characteristics of numerical solution was not affected by the number of subdomains as much as the transient solution was dependent on the number of subdomains.

Numerical study on the influence of mesomechanical properties on macroscopic fracture of concrete

  • Zhu, W.C.;Tang, C.A.;Wang, S.Y.
    • Structural Engineering and Mechanics
    • /
    • 제19권5호
    • /
    • pp.519-533
    • /
    • 2005
  • The numerical simulations on the influence of mesoscopic structures on the macroscopic strength and fracture characteristics are carried out based on that the concrete is assumed to be a three-phase composite composed of matrix (mortar), aggregate and bond between them by using a numerical code named MFPA. The finite element program is employed as the basic stress analysis tool when the elastic damage mechanics is used to describe the constitutive law of meso-level element and the maximum tensile strain criterion and Mohr-Coulomb criterion are utilized as damage threshold. It can be found from the numerical results that the bond between matrix and aggregate has a significant effect on the macroscopic mechanical performance of concrete.

Numerical Simulation of the Effects of Moisture on the Reinforcement of a Tropopause Fold

  • Lee, Hong-Ran;Kim, Kyung-Eak;Lee, Yong-Hee
    • 한국지구과학회지
    • /
    • 제30권5호
    • /
    • pp.630-645
    • /
    • 2009
  • The tropopause fold event that took place on January 1, 1997 over mid-region on the Korean Peninsula is examined by means of a numerical simulation based on a Mesoscale Model (MM5). The purpose of this paper is to investigate the effects of moisture in reinforcing a tropopause fold linked to an explosive cyclone. Two types of simulations were carried out; 1) simulations for moist conditions in which full physical and dynamic processes are considered and 2) simulations for dry conditions in which cumulus parameterization and cloud microphysics process are excluded. The results of the moist condition simulations demonstrate that the intensity of the central pressure of the cyclone was overestimated compared with the observed values and that the location of the center and the pressure deepening rates (-17 hPa/12 hr) complied with the observed values. The potential vorticity (PV) anomaly on the isentropic surface at 305 K continued to move in a southeast direction on January 1, 1997 and thus created a single tube of tropopause fold covering the northern and the middle area of the Korean Peninsula and reaching the ground surface at 0300 UTC and 0600 UTC. The results of the dry condition simulations show that the tropopause descended to 500 and 670 hPa in 0300 and 0600 UTC, respectively at the same location for the moist condition simulation; however, there was no deep tropopause fold observed. A comparison of the simulated data between the moist and the dry conditions suggests that a deep tropopause fold should happen when there is sufficient moist in the atmosphere and significantly large PV in the lower atmosphere pulls down the upper atmosphere rather than when the tropopause descends itself due to dynamic causes. Thus, it is estimated that moisture in the atmosphere should have played a crucial role in a deep tropopause fold process.

주기 발열 파형을 이용한 열식 질량 유량계의 특성에 관한 수치적 연구 (A numerical study on the characteristics of a thermal mass air flow sensor with periodic heating pulses)

  • 전홍규;오동욱;박병규;이준식
    • 대한기계학회:학술대회논문집
    • /
    • 대한기계학회 2007년도 춘계학술대회B
    • /
    • pp.2482-2487
    • /
    • 2007
  • Numerical simulations are conducted for the analysis of a thermal mass air flow sensor with periodic heating pulses on silicon-nitride ($Si_3N_4$) thin membrane structure. This study aims to find the locations of temperature sensors on the thin membrane and the heating pulse conditions, that the higher sensitivity can be achieved, for the development of a MEMS fabricated mass air flow sensor which is driven in periodic heating pulse. The simulations, thus, focus on the membrane temperature profile according to variation of the flow velocity, heating duration time and imposed power. The flow velocity of the simulations is ranging from 3 m/s to 35 m/s, heating duration time from 1 ms to 3 ms and imposed power from 50 mW to 90 mW. The corresponding Reynolds numbers vary from 1000 to 10000.

  • PDF

Tracing the trajectory of pelagic Sargassum using satellite monitoring and Lagrangian transport simulations in the East China Sea and Yellow Sea

  • Kwon, Kyungman;Choi, Byoung-Ju;Kim, Kwang Young;Kim, Keunyong
    • ALGAE
    • /
    • 제34권4호
    • /
    • pp.315-326
    • /
    • 2019
  • Northeastward drifts of massive Sargassum patches were observed in the East China Sea (ECS) and Yellow Sea (YS) by the Geostationary Ocean Color Imager (GOCI) in May 2017. Coverage of the brown macroalgae patches was the largest ever recorded in the ECS and YS. Three-dimensional circulation modeling and Lagrangian particle tracking simulations were conducted to reproduce drifting trajectories of the macroalgae patches. The trajectories of the macroalgae patches were controlled by winds as well as surface currents. A windage (leeway) factor of 1% was chosen based on sensitivity simulations. Southerly winds in May 2017 contributed to farther northward intrusion of the brown macroalgae into the YS. Although satellite observation and numerical modeling have their own limitations and associated uncertainties, the two methods can be combined to find the best estimate of Sargassum patch trajectories. When satellites were unable to capture all patches because of clouds and sea fog in the ECS and YS, the Lagrangian particle tracking model helped to track and restore the missing patches in satellite images. This study suggests that satellite monitoring and numerical modeling are complementary to ensure accurate tracking of macroalgae patches in the ECS and YS.