• Computers and Concrete
• /
• v.8 no.2
• /
• pp.163-176
• /
• 2011

By virtue of chord-length density function from the field of statistical physics, this paper introduced a quantitative approach to estimate the distribution of cement paste thickness between aggregates in concrete. Dynamics mixing method based on molecular dynamics was employed to generate one model structure, then image analysis algorithm was used to obtain the distribution of thickness of cement paste in model structure for the purpose of verification. By comparison of probability density curves and cumulative probability curves of the cement paste thickness among neighboring aggregates, it is found that the theoretical results are consistent with the simulation. Furthermore, for the model mortar and concrete mixtures with practical volume fraction of Fuller-type aggregate, this analytical formula was employed to predict the influence of aggregate volume fraction and aggregate fineness. And evolution of its mean values were also investigated with the variation of volume fraction of aggregate as well as the fineness of aggregates in model mortars and concretes.

• The Bulletin of The Korean Astronomical Society
• /
• v.36 no.2
• /
• pp.106.2-106.2
• /
• 2011

The two competing theories of star formation are based on turbulence and ambipoar diffusion. I will first briefly explain the two theories. There have been analytical (or semi-analytic) models, which estimate star formation rates in a turbulent cloud. Most of them are based on the log-normal density PDF (probability density function) of the turbulent cloud without self-gravity. I will first show that the core (star) formation rate can be increased significantly once self-gravity of a turbulence cloud is taken into account. I will then present the evolution of molecular line profiles of HCO+ and C18O toward a dense core that is forming inside a magnetized turbulent molecular cloud. Features of the profiles can be affected more significantly by coupled velocity and abundance structures in the outer region than those in the inner dense part of the core. During the evolution of the core, the asymmetry of line profiles easily changes from blue to red, and vice versa. Finally, I will introduce a method for incorporating ambipolar diffusion in the strong coupling approximation into a multidimensional magnetohydrodynamic code.

• Wind and Structures
• /
• v.35 no.4
• /
• pp.287-296
• /
• 2022

Analyzing the typhoon wind hazards is crucial to determine the extreme wind load on engineering structures in the typhoon prone region. In essence, the typhoon hazard analysis is a high-dimensional problem with randomness arising from the typhoon genesis, environmental variables and the boundary layer wind field. This study suggests a dimension reduction approach by decoupling the original typhoon hazard analysis into two stages. At the first stage, the randomness of the typhoon genesis and environmental variables are propagated through the typhoon track model and intensity model into the randomness of the key typhoon parameters. At the second stage, the probability distribution information of the key typhoon parameters, combined with the randomness of the boundary layer wind field, could be used to estimate the extreme wind hazard. The Chinese southeast coastline is taken as an example to demonstrate the adequacy and efficiency of the suggested decoupling approach.

• Journal of Mechanical Science and Technology
• /
• v.18 no.8
• /
• pp.1470-1478
• /
• 2004

In this paper, a new 'presumed' Probability Density Function (PDF) approach coupled with a Lagrangian tracking method is proposed for turbulent combustion modeling. The test and the investigation of the model are conducted by comparing the model results with DNS data for a premixed flame subjected in a decaying turbulent field. The newly constructed PDF, which incorporates the instantaneous chemical reaction term, demonstrates consistent improvement over conventional assumed PDF models. It has been found that the time evolution of the mean scalar, the variance and the mean reaction rate are strongly influenced by a parameter deduced by a Lagrangian equation which takes into account explicitly the local reaction rate. Tests have been performed for a moderate Damkohler number, and it is expected the model may cover a broader range of Damkohler number. The comparison with the DNS data demonstrates that the proposed model may be promising and affordable for implementation in a moment-equation solver.

• Journal of the Korean Society of Combustion
• /
• v.14 no.2
• /
• pp.32-41
• /
• 2009

The transported probability density function model combined with the consistent finite volume (FV) method has been applied to simulate the turbulent bluff-body reacting flows. To realistically account for the non-isotropic turbulence effects on the turbulent bluff-body reacting flows, the present PDF transport approach is based on the joint velocity- turbulent frequency-composition PDF formulation. The evolution of the fluctuating velocity of a particle is modeled by a simplified Langevin equation and the particle turbulence frequency is represented by the modified Jayesh - Pope model. Effects of molecular diffusion are represented by the interaction by exchange with the mean (IEM) mixing model. To validate this hybrid FV/PDF transport model, the numerical results are compared with experimental data for the turbulent bluff-body reacting flows.

• Honam Mathematical Journal
• /
• v.44 no.3
• /
• pp.402-418
• /
• 2022

Noise is a fundamental factor to increased validity and regularity of spike propagation and neuronal firing in the nervous system. In this paper, we examine the stochastic version of the Izhikevich-FitzHugh neuron dynamical model. This approach is based on techniques presented by Luo and Guo, which provide a general framework for the bifurcation and stability analysis of two dimensional stochastic dynamical system as an Itô averaging diffusion system. By using largest lyapunov exponent, local and global stability of the stochastic system at the equilibrium point are investigated. We focus on the two kinds of stochastic bifurcations: the P-bifurcation and the D-bifurcations. By use of polar coordinate, Taylor expansion and stochastic averaging method, it is shown that there exists choices of diffusion and drift parameters such that these bifurcations occurs. Finally, numerical simulations in various viewpoints, including phase portrait, evolution in time and probability density, are presented to show the effects of the diffusion and drift coefficients that illustrate our theoretical results.

• Nuclear Engineering and Technology
• /
• v.50 no.5
• /
• pp.673-682
• /
• 2018

One of the most widely used methods to estimate core damage during a nuclear power plant accident is containment radiation measurement. The evolution of severe accidents is extremely complex, leading to uncertainty in the containment dose rate (CDR). Therefore, it is difficult to accurately determine core damage. This study proposes to conduct uncertainty analysis of CDR for core damage assessment. First, based on source term estimation, the Monte Carlo (MC) and point-kernel integration methods were used to estimate the probability density function of the CDR under different extents of core damage in accident scenarios with late containment failure. Second, the results were verified by comparing the results of both methods. The point-kernel integration method results were more dispersed than the MC results, and the MC method was used for both quantitative and qualitative analyses. Quantitative analysis indicated a linear relationship, rather than the expected proportional relationship, between the CDR and core damage fraction. The CDR distribution obeyed a logarithmic normal distribution in accidents with a small break in containment, but not in accidents with a large break in containment. A possible application of our analysis is a real-time core damage estimation program based on the CDR.

• Coupled systems mechanics
• /
• v.11 no.2
• /
• pp.167-198
• /
• 2022

In this paper we deal with classical instability problems of heterogeneous Euler beam under conservative loading. It is chosen as the model problem to systematically present several possible solution methods from simplest deterministic to more complex stochastic approach, both of which that can handle more complex engineering problems. We first present classical analytic solution along with rigorous definition of the classical Euler buckling problem starting from homogeneous beam with either simplified linearized theory or the most general geometrically exact beam theory. We then present the numerical solution to this problem by using reduced model constructed by discrete approximation based upon the weak form of the instability problem featuring von Karman (virtual) strain combined with the finite element method. We explain how such numerical approach can easily be adapted to solving instability problems much more complex than classical Euler's beam and in particular for heterogeneous beam, where analytic solution is not readily available. We finally present the stochastic approach making use of the Duffing oscillator, as the corresponding reduced model for heterogeneous Euler's beam within the dynamics framework. We show that such an approach allows computing probability density function quantifying all possible solutions to this instability problem. We conclude that increased computational cost of the stochastic framework is more than compensated by its ability to take into account beam material heterogeneities described in terms of fast oscillating stochastic process, which is typical of time evolution of internal variables describing plasticity and damage.

• 한국신재생에너지학회:학술대회논문집
• /
• 2007.11a
• /
• pp.417-423
• /
• 2007

Oxy-gasification or oxygen-blown gasification, enables a clean and efficient use of coal and opens a promising way to CO2 capture. The coal gasification process of a slurry feed type, entrained-flow coal gasifier was numerically predicted in this paper. The purposes of this study are to develop an evaluation technique for design and performance optimization of coal gasifiers using a numerical simulation technique, and to confirm the validity of the model. By dividing the complicated coal gasification process into several simplified stages such as slurry evaporation, coal devolatilization, mixture fraction model and two-phase reactions coupled with turbulent flow and two-phase heat transfer, a comprehensive numerical model was constructed to simulate the coal gasification process. The influence of turbulence on the gas properties was taken into account by the PDF (Probability Density Function) model. A numerical simulation with the coal gasification model is performed on the Conoco-Philips type gasifier for IGCC plant. Gas temperature distribution and product gas composition are also presented. Numerical computations were performed to assess the effect of variation in oxygen to coal ratio and steam to coal ratio on reactive flow field. The concentration of major products, CO and H2 were calculated with varying oxygen to coal ratio (0.2-1.5) and steam to coal ratio(0.3-0.7). To verify the validity of predictions, predicted values of CO and H2 concentrations at the exit of the gasifier were compared with previous work of the same geometry and operating points. Predictions showed that the CO and H2 concentration increased gradually to its maximum value with increasing oxygen-coal and hydrogen-coal ratio and decreased. When the oxygen-coal ratio was between 0.8 and 1.2, and the steam-coal ratio was between 0.4 and 0.5, high values of CO and H2 were obtained. This study also deals with the comparison of CFD (Computational Flow Dynamics) and STATNJAN results which consider the objective gasifier as chemical equilibrium to know the effect of flow on objective gasifier compared to equilibrium. This study makes objective gasifier divided into a few ranges to study the evolution of the gasification locally. By this method, we can find that there are characteristics in the each scope divided.