• Water for future
• /
• v.29 no.4
• /
• pp.187-198
• /
• 1996

A numerical model for unsteady dispersion of horizontal line source in turbulent shear flow is developed. A fractional step finite difference method is used which splits the unsteady two-dimensional advective diffusion equation into the longitudinal advection and the vertical diffusion equations, and solves them alternately for half time intervals by the Holly-Preissmann scheme and the Crank-Nicholson scheme, respectively. The developed numerical model is verified using a semi-analytic solution for steady dispersion in turbulent shear flow. Dispersion of an instantaneous plane source in turbulent shear flow is analyzed using the model. The degree of mixing at the same dimensionless time is almost the same regardless of the friction factor, and the travel distance required to reach a certain degree of mixing is inversely proportional to the square root of the friction factor.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.26 no.5B
• /
• pp.551-555
• /
• 2006

In this study, new dispersion-correction terms are added to a leap-frog finite difference scheme for the linear shallow-water equations with the purpose of considering dispersion effects of the linear Boussinesq equations for propagation of tsunamis. The numerical model developed in this study is tested to the problem that the initial free surface displacement is a Gaussian hump over a constant water depth, and the predicted numerical results are compared with analytical solutions. The results of the present numerical model are accurate in comparison with those of existing models.

• Journal of the Korean Chemical Society
• /
• v.47 no.1
• /
• pp.31-37
• /
• 2003

The expanded uncertainties calculated by the application of GUM -approximation and Monte-Carlo simulation were compared about the model equation of one-point calibration which is widely used for the measurements and chemical analysis. For the comparisons, we assumed a set of artificial data at the various level of concentration and dispersion of t or normal distribution. Mistakes of more then 50 % was revealed at the values calculated by GUM-approximation in comparison with those of Monte-Carlo simulation because of the excess dispersion from t-distribution and non-linearity by division in the equation. In contrary, the mistake of calculation due to non-linearity of the equation was not observed in the level of detection limits with the equation of one-point calibration, because of the relatively large values of uncertainty in response.

• Coupled systems mechanics
• /
• v.12 no.1
• /
• pp.41-68
• /
• 2023

The present work is concerned with the study of the influence of inhomogeneous initial stresses in a hollow cylinder containing a compressible inviscid fluid on the propagation of axisymmetric longitudinal waves propagating in this cylinder. The study is carried out using the so-called three-dimensional linearized theory of elastic waves in bodies with initial stresses to describe the motion of the cylinder and using the linearized Euler equations to describe the flow of the compressible inviscid fluid. It is assumed that the inhomogeneous initial stresses in the cylinder are caused by the internal pressure of the fluid. To solve the corresponding eigenvalue problem, the discrete-analytic solution method is applied and the corresponding dispersion equation is obtained, which is solved numerically, after which the corresponding dispersion curves are constructed and analyzed. To obtain these dispersion curves, parameters characterizing the magnitude of the internal pressure, the ratio of the sound velocities in the cylinder material and in the fluid, and the ratio of the material densities of the fluid and the cylinder are introduced. Based on these parameters, the influence of the inhomogeneous initial stresses in the cylinder on the dispersion of the above-mentioned waves in the considered hydro-elastic system is investigated. Moreover, based on these results, appropriate conclusions about this influence are drawn. In particular, it is found that the character of the influence depends on the wavelength. Accordingly, the inhomogeneous initial stresses before (after) a certain value of the wavelength lead to a decrease (increase) of the wave propagation velocity in the zeroth and first modes.

• Geomechanics and Engineering
• /
• v.15 no.1
• /
• pp.645-657
• /
• 2018

The paper aims to analyze the behaviour of torsional type surface waves propagating through fluid saturated inhomogeneous porous media clamped between two inhomogeneous anisotropic media. We considered three types of inhomogeneities in upper anisotropic layer which varies exponentially, quadratically and hyperbolically with depth. The anisotropic half space inhomogeneity varies linearly with depth and intermediate layer is taken as inhomogeneous fluid saturated porous media with sinusoidal variation. Following Biot, the dispersion equation has been derived in a closed form which contains Whittaker's function and its derivative, for approximate result that have been expanded asymptotically up to second term. Possible particular cases have been established which are in perfect agreement with standard results and observe that when one of the upper layer vanishes and other layer is homogeneous isotropic over a homogeneous half space, the velocity of torsional type surface waves coincides with that of classical Love type wave. Comparative study has been made to identify the effects of various dimensionless parameters viz. inhomogeneity parameters, anisotropy parameters, porosity parameter, and initial stress parameters on the torsional wave propagation by means of graphs using MATLAB. The study has its own relevance in connection with the propagation of seismic waves in the earth where fluid saturated poroelastic layer is present.

• Journal of the Korean Society of Safety
• /
• v.10 no.3
• /
• pp.68-80
• /
• 1995

A series of parametric calculations have been performed in order to investigate the short-term and short-range plume and puff behavior of toxic gaseous and solid pollutant dispersion in an open atmosphere. The simulation is made by the use of the computer program developed by this laboratory, in which a control-volume based finite-difference method is used together with the SIMPLEC algorithm for the resolution of the pressure-velocity coupling appeared In Wavier-Stokes equation. The Reynolds stresses are solved by the standard two-equation k-$\varepsilon$ model modified for buoyancy together with the RNG(Renormalization Group) k-$\varepsilon$ model. The major parameters considered in this calculation are pollutant gas density and temperature, the relative velocity of pollutants to that of the surrounding atmospheric air, and particulate size and density together with the height released. The flow field is typically characterized by the formation of a strong recirculation region for the case of the low density gases such as $CH_4$ and air due to the strong buoyancy, while the flow is simply declining pattern toward the downstream ground for the case of heavy molecule like the $CH_2C1_2$and $CCl_4$, even for the high temperature, $200^{\circ}C$. The effect of gas temperature and velocity on the flow field together with the particle trajectory are presented and discussed in detail. In general, the results are physically acceptable and consistent.

• Composites Research
• /
• v.32 no.3
• /
• pp.141-147
• /
• 2019

In this study, micro-sized and nano-sized pure aluminum (Al) powders were compressed by unidirectional pressure at room temperature. Although neither type of Al bulk was heated, they had a high relative density and improved mechanical properties. The microstructural analysis showed a difference in the process of densification according to particle size, and the mechanical properties were measured by the Vickers hardness test and the nano indentation test. The Vickers hardness of micro Al and nano Al fabricated in this study was five to eight times that of ordinary Al. The grain refinement effect was considered to be one of the strengthening factors, and the Hall-Petch equation was introduced to analyze the improved hardness caused by grain size reduction. In addition, the effect of particle size and dispersion of aluminum oxide in the bulk were additionally considered. Based on these results, the present study facilitates the examination of the effect of particle size on the mechanical properties of compacted bulk fabricated by the powder metallurgy method and suggests the possible way to improve the mechanical properties of nano-crystalline powders.

• Journal of Korean Port Research
• /
• v.12 no.1
• /
• pp.75-86
• /
• 1998

To design a coastal structure in the nearshore region, engineers must have means to estimate wave climate. Waves, approaching the surf zone from offshore, experience changes caused by combined effects of bathymetric variations, interference of man-made structure, and nonlinear interactions among wave trains. This paper has attempted to find out the effects of two of the more subtle phenomena involving nonlinear shallow water waves, amplitude dispersion and secondary wave generation. Boussinesq-type equations can be used to model the nonlinear transformation of surface waves in shallow water due to effect of shoaling, refraction, diffraction, and reflection. In this paper, generalized Boussinesq equations under the complex bottom condition is derived using the depth averaged velocity with the series expansion of the velocity potential as a product of powers of the depth of flow. A time stepping finite difference method is used to solve the derived equation. Numerical results are compared to hydraulic model results. The result with the non-linear dispersive wave equation can describe an interesting transformation a sinusoidal wave to one with a cnoidal aspect of a rapid degradation into modulated high frequency waves and transient secondary waves in an intermediate region. The amplitude dispersion of the primary wave crest results in a convex wave front after passing through the shoal and the secondary waves generated by the shoal diffracted in a radial manner into surrounding waters.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.33 no.2
• /
• pp.495-506
• /
• 2013

This study developed the Time-split Mixing Model (TMM) which can represent the pollutant mixing process on a three-dimensional open channel through constructing the conceptual model based on Taylor's assumption (1954) that the shear flow dispersion is the result of combination of shear advection and diffusion by turbulence. The developed model splits the 2-D mixing process into longitudinal mixing and transverse mixing, and it represents the 2-D advection-dispersion by the repetitive calculation of concentration separation by the vertical non-uniformity of flow velocity and then vertical mixing by turbulent diffusion sequentially. The simulation results indicated that the proposed model explains the effect of concentration overlapping by boundary walls, and the simulated concentration was in good agreement with the analytical solution of the 2-D advection-dispersion equation in Taylor period (Chatwin, 1970). The proposed model could explain the correlation between hydraulic factors and the dispersion coefficient to provide the physical insight about the dispersion behavior. The longitudinal dispersion coefficient calculated by the TMM varied with the mixing time unlike the constant value suggested by Elder (1959), whereas the transverse dispersion coefficient was similar with the coefficient evaluated by experiments of Sayre and Chang (1968), Fischer et al. (1979).

• Proceedings of the Korean Society of Soil and Groundwater Environment Conference
• /
• 2002.04a
• /
• pp.122-126
• /
• 2002

Contaminant transport in porous media is characterized by solving an advection-dispersion equation(ADE). The ADE can cover equilibrium phenomena of interest, which include sorption, decay, and chemical reactions. Among these phenomena, sorption mechanism is described by several types of sorption isotherm. If we assume the sorption isotherm as linear, the solution of ADE can be easily procured. However, if we consider the sorption isotherm as non-linear isotherm like a Dual Reactive Domain Model (DRDM), the resulting differential equation becomes non-linear. In this case, the solution of ADE cannot be easily acquired by an analytic method. In this paper, we present the numerical analysis of ADE using a DRDM. The results reveal that even if sorption data may be fitted well using linear or non-linear isotherm, the characteristics of contaminant transport of the two cases are different from each other. To be concrete, the retardation of linear isotherm has stronger effect than that of the DRDM. As the non-linearity of sorption isotherm increases, the difference of retardation effects of the two cases becomes larger. For a pulse source, the maximum concentration of the linear model is higher than that of the DRDM, but the plume of the DRDM moves faster than that of the linear model. Behaviors of contaminant transport using the DRDM are consistent with common features of a linear model. For instance, biodegradation effect becomes larger as time goes by The faster the seepage velocity is, the faster the plume of contaminant moves. The plume of the contaminant is distributed evenly over overall domain in the event of high dispersion coefficient.