• 천문학회지
• /
• 제40권4호
• /
• pp.91-97
• /
• 2007

The Fokker-Planck (FP) model is one of the commonly used methods for studies of the dynamical evolution of dense spherical stellar systems such as globular clusters and galactic nuclei. The FP model is numerically stable in most cases, but we find that it encounters numerical difficulties rather often when the effects of tidal shocks are included in two-dimensional (energy and angular momentum space) version of the FP model or when the initial condition is extreme (e.g., a very large cluster mass and a small cluster radius). To avoid such a problem, we have developed a new integration scheme for a two-dimensional FP equation by adopting an Alternating Direction Implicit (ADI) method given in the Douglas-Rachford split form. We find that our ADI method reduces the computing time by a factor of ${\sim}2$ compared to the fully implicit method, and resolves problems of numerical instability.

• Nuclear Engineering and Technology
• /
• 제56권8호
• /
• pp.2974-2989
• /
• 2024

The one-dimensional Two-Fluid concept uses an area-average approach to simplify the time and phase-averaged Two-Fluid conservation equations, making it more suitable for addressing difficulties at an industrial scale. Nevertheless, the mathematical framework has inherent weaknesses due to the loss of details throughout the averaging procedures. This limitation makes the conventional model inappropriate for some flow regimes, where short-wavelength perturbations experience uncontrolled amplification, leading to solutions that need to be physically accurate. The critical factor in resolving this problem is the integration of closure relations. These relationships play a crucial function in reintroducing essential physical characteristics, thus correcting the loss that occurs during averaging and guaranteeing the stability of the model. To improve the accuracy of predictions, it is essential to assess the stability and grid dependence of one-dimensional formulations, which are particularly affected by closure relations and numerical schemes. The current research presented in the text focuses on improving the well-posedness of the TFM, specifically within the TRACE code, which is widely utilized for nuclear reactor safety assessments. Incorporating a bubble collision model in the momentum equations is demonstrated to enhance the TFM's resilience, especially in scenarios with high void fractions where conventional TFMs may face challenges. The analysis presents a linear stability analysis performed for the transient one-dimensional Two-Fluid Model of system code TRACE within the framework of vertically dispersed flows. The main emphasis is on evaluating the stability characteristics of the model while also acknowledging its susceptibility to closure relations and numerical techniques.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2006년도 PARALLEL CFD 2006
• /
• pp.311-314
• /
• 2006

The supersonic flow around tandem cavities was investigated by three- dimensional numerical simulations using the Reynolds-Averaged Navier-Stokes(RANS) equation with the $\kappa-\omega$ thrbulence model. The flow around a cavity is characterized as unsteady flow because of the formation and dissipation of vortices due to the interaction between the freestream shear layer and cavity internal flow, the generation of shock and expansion waves, and the acoustic effect transmitted from wake flow to upstream. The upwind TVD scheme based on the flux vector split using van Leer's limiter was used as the numerical method. Numerical calculations were performed by the parallel processing with time discretizations carried out by the 4th-order Runge-Kutta method. The aspect ratio of cavities are 3 for the first cavity and 1 for the second cavity. The ratio of cavity interval to depth is 1. The ratio of cavity width to depth is 1 in the case of three dimensional flow. The Mach number and the Reynolds number were 1.5 and $4.5{\times}10^5$, respectively. The characteristics of the dominant frequency between two-dimensional and three-dimensional flows were compared, and the characteristics of the second cavity flow due to the fire cavity flow cavity flow was analyzed. Both two dimensional and three dimensional flow oscillations were in the 'shear layer mode', which is based on the feedback mechanism of Rossiter's formula. However, three dimensional flow was much less turbulent than two dimensional flow, depending on whether it could inflow and outflow laterally. The dominant frequencies of the two dimensional flow and three dimensional flows coincided with Rossiter's 2nd mode frequency. The another dominant frequency of the three dimensional flow corresponded to Rossiter's 1st mode frequency.

• Journal of Mechanical Science and Technology
• /
• 제20권8호
• /
• pp.1256-1265
• /
• 2006

The supersonic flows around tandem cavities were investigated by two-dimensional and three-dimensional numerical simulations using the Reynolds-Averaged Navier-Stokes (RANS) equation with the k- ω turbulence model. The flow around a cavity is characterized as unsteady flow because of the formation and dissipation of vortices due to the interaction between the freestream shear layer and cavity internal flow, the generation of shock and expansion waves, and the acoustic effect transmitted from wake flow to upstream. The upwind TVD scheme based on the flux vector split with van Leer's limiter was used as the numerical method. Numerical calculations were performed by the parallel processing with time discretizations carried out by the 4th-order Runge- Kutta method. The aspect ratios of cavities are 3 for the first cavity and 1 for the second cavity. The ratio of cavity interval to depth is 1. The ratio of cavity width to depth is 1 in the case of three dimensional flow. The Mach number and the Reynolds number were 1.5 and $4.5{\times}10^5$, respectively. The characteristics of the dominant frequency between two- dimensional and three-dimensional flows were compared, and the characteristics of the second cavity flow due to the first cavity flow was analyzed. Both two dimensional and three dimensional flow oscillations were in the 'shear layer mode', which is based on the feedback mechanism of Rossiter's formula. However, three dimensional flow was much less turbulent than two dimensional flow, depending on whether it could inflow and outflow laterally. The dominant frequencies of the two dimensional flow and three dimensional flows coincided with Rossiter's 2nd mode frequency. The another dominant frequency of the three dimensional flow corresponded to Rossiter's 1st mode frequency.

• 한국농공학회:학술대회논문집
• /
• 한국농공학회 2000년도 학술발표회 발표논문집
• /
• pp.162-168
• /
• 2000

A two-dimensional numerical model based on a finite volume method was formulated to solve the shallow water equations and applied for evaluating irrigation and drainage characteristics at large-sized paddy fields. Manning roughness coefficient was calibrated using the observed inundating depths at drainage tests and used for validating the model with the results from irrigation and drainage test. The simulated results were in good agreement with the observed inundating depths.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2015년도 학술발표회
• /
• pp.240-240
• /
• 2015

As Simpson et al. (1990) emphasized the importance of the straining process in the stratification and mixing in the estuarine circulation process, various researches have investigated on the relative contribution of each process to the overall potential energy anomaly dynamics. However, many numerical works have done only for two dimensional modeling along channel or the short distance cross sectional three dimensional simulations as Burchard et al. (2008) and the estuarine channel was not simulated so far. But, in the study on the physics of shallow coastal seas, spatial dimension in the three dimensional way affects significantly on results of a particular numerical model. Therefore, the comparison of two and three dimensional models is important to understand the real physics of mixing and stratification in an estuary. Also, as Geyer and MacCready (2013) pointed out that the lateral process seems to be important in determining the periodic stratifications, to study such process the three dimensional modeling must be required. The present study uses a numerical model to show the signification roles of each term of the time-dependent dynamic equation for the potential energy anomaly (PEA) in controlling along and lateral channel flows and different stratification structures. Moreover, we present the relationships between the ${\Phi}$-advection, the depth mean straining, vertical mixing and vertical advection can explain well how water level, salinity distribution and across velocity 2D model are slightly different from 3D.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2002년도 학술대회지
• /
• pp.709-710
• /
• 2002

Cancer angiogenesis is simulated using a two dimensional model. Governing equation of angiogenesis is a TAE (Tumor angiogenesis factor) conservation equation in time and space. A stochastic process model is utilized to simulate vessel formation, proliferation, and migration to a cancer pellet. Numerical results are presented especially in case of growing cancer.

• Advances in environmental research
• /
• 제5권4호
• /
• pp.263-282
• /
• 2016

Contaminant transport in groundwater induces major threat and harmful effect on the environment; hence, the fate of the contaminant migration in groundwater is seeking a lot of attention. In this paper a two dimensional numerical flow and transport model through saturated layered soil is developed. Groundwater flow and solute transport has been simulated numerically using proposed model. The model implements the finite volume time splitting method to discretize the main equations. The performance, accuracy and efficiency of the out coming numerical models have been successfully examined by two test cases. The verification test cases consist of two-dimensional, groundwater flow and solute transport. The final purpose of this paper is to discuss and compare the shape of contaminant plume in homogeneous and heterogeneous media with different soil properties and control of solute transport using a zone for minimizing the potential of groundwater contamination; furthermore, this model leads to select the effective and optimum remedial strategies for cleaning the contaminated aquifers.

• Water Engineering Research
• /
• 제2권3호
• /
• pp.151-160
• /
• 2001

A two-dimensional flow model is newly developed. Two-dimensional shallow-water equations are discretized by the finite volume method. A nonorthogonal coordinate system is then employed. The developed model is applied to simulations of flows in a 180 degree curved bend flow. Numerical prediction are compared to available laboratory measurement. A good agreement is observed.

• 한국연소학회:학술대회논문집
• /
• 한국연소학회 2006년도 제32회 KOSCO SYMPOSIUM 논문집
• /
• pp.89-95
• /
• 2006

This study has numerically modelled the combustion processes of the turbulent swirling premixed lifted flames in the low-swirl burner (LSB). In these turbulent swirling premixed flames, the four tangentially-injected air jets induce the turbulent swirling flow which plays the crucial role to stabilize the turbulent lifted flame. In the present approach, the turbulence-chemistry interaction is represented by the level-set based flamelet model. Two-dimensional and three-dimensional computations are made for the various swirl numbers and nozzle length. In terms of the centerline velocity profiles and flame liftoff heights, numerical results are compared with experimental data The three-dimensional approach yields the much better conformity with agreements with measurements without any analytic assumptions on the inlet swirl profiles, compared to the two-dimensional approach. Numerical clearly results indicate that the present level-set based flamelet approach has realistically simulated the structure and stabilization mechanism of the turbulent swirling stoichiometric and lean-premixed lifted flames in the low-swirl burner.