• Proceedings of the KIEE Conference
• /
• 2004.07d
• /
• pp.2209-2211
• /
• 2004

Unlike the wavelet neural network, since a mother wavelet layer of the self-recurrent wavelet neural network (SRWNN) is composed of self-feedback neurons, it has the ability to store past information of the wavelet. Therefore we propose the prediction method for the nonlinear chaotic time series model using a SRWNN. The SRWNN model is learned for the modeling of a function such that the inputs arc known values of the time series and the output is the value in the future. The parameters of the network are tuned to minimize the difference between the nonlinear mapping of the chaotic time series and the output of SRWNN using the gradient-descent method for the adaptive backpropagation algorithm. Through the computer simulations, we demonstrate the feasibility and the effectiveness of our method for the prediction of the logistic map and the Mackey-Glass delay-differential equation as a nonlinear chaotic time series.

• Computers and Concrete
• /
• v.15 no.1
• /
• pp.127-140
• /
• 2015

The process of chloride ingress in saturated concrete was presented by a previous study that used a mathematical model for the same as that concrete. This model is to be studied chloride ion diffusion which is considered as a chemical phenomenon and is to be represented the chloride diffusion process to be a nonlinear partial differential equation (PDE). In this paper, this nonlinear PDE is solved by the Kirchhoff transformation to render into a linear PDE. This linear PDE associated with initial and boundary conditions is also solved by the Laplace transformation to obtain an analytical solution. To verify the serviceability and reliability of this proposed method, the practical application should be supplied. The input parameters were cited from the previous study. The free chloride concentration profiles obtained by the analytical solution of mathematical model for saturated concretes after 24 and 120 hrs of exposure were compared with the previous study. The predicted results obtained from proposed method have a tendency with experimental results obtained by the previous study and trend toward numerical results approximated by finite difference technique.

• 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.

• Proceedings of the Korean Society for Agricultural Machinery Conference
• /
• 1993.10a
• /
• pp.443-452
• /
• 1993

Direct freezing tests of soybean seed by liquid nitrogen were carried out at various moisture contents and the following important conclusions were drawn from the results of temperature measurements of soybean seed and photographs of bubbles generated on its surface : 1) Assuming that the temperature gradient in a soybean seed is negligible because of its small seed size and the freezing ratio is followed the Heiss's formula, and a differential equation based on the heat energy balance was introduced . The equation was easily solved by the Runge-Kutta-Gill method and the predicted values of the temperature were in good agreement with the observed data. 2) The photographs of bubble generation during freezing showed the boiling mode was nucleate, and then the most suitable formula on the nucleate boiling heat transfer was introduced from many formulate proposed up to now by fitting the calculated values based on the formula to the observed data. The formula used for the predict on of the seed temperature was as follows: $\frac{{\partial}T_s}{\partial\theta}\;=\;-\frac{{\alpha}(T_s\;-\;T_L)^{3.3}}{W(C_s\;-\;\frac{{\delta}m(CT_s\;+\;{\sigma})}{T_s^2})}$ where C = difference of the specific heat between pure ice and water m=moisture content of soybean seed $T_s$ = seed temperature $T_L$ = Temperature of liquid nitrogen W = mass of soybean seed $\alpha$ = proportional constant $\delta$ = constant depends on variety or the type of seed $\theta$ = time $\sigma$ = latent heat of melting of pure ice This study will give important information in the hydro-freezing technique by liquid nitrogen, available as a new technique of processing agricultural products in the near future.

• Journal of Korean Society of Steel Construction
• /
• v.10 no.4 s.37
• /
• pp.615-627
• /
• 1998

In order to determine the maximum shear force, the maximum bending moment, the maximum pure torsion. the maximum warping torsion, and the maximum bimoment for the curved girder grid bridges, it is important to find the location of live load applied to the curved girder grid bridges, so that the influence line can be estimated. The fundamental differential equation concerning the behaviour with warping effects for the curved girder is developed by Vlasov. In this paper, the influence line of shear force, bending moment, pure torsion, warping torsion, and bimoment due to unit vertical load and unit torsional moment for curved I-girder grid bridges with variable and constant cross section are obtained by using the finite difference method and compared with respectively.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.2
• /
• pp.61-71
• /
• 1993

The stress waves generated by the mechanical energies by impact or the chemical energies by the explosions are transmitted through medium. The wave propagation process through medium is a very complicated procedure due to the reflections and refractions of the waves at the free surfaces and interfaces. In this study the pressure independent Von-Mises model is employed for the wave propagation analysis in the layered systems. Governing equations of this study are conservation equations of momentum and mass in Lagrangian coordinate system which is fixed to the material. Due to the shock-front which violates the continuity assumptions inherent in the differential equations numerical artificial viscosity is used to spread the shock front over several computational zones. These equations are solved by Finite Difference Method with discretized time and space coordinates. The associate normality flow rule as a plastic theory is implemented to find the plastic strains.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.12 no.2
• /
• pp.21-33
• /
• 1992

Under the intense-impulsive loading, structures are subjected to the wide range of pressures at an instantaneous time. The constitutive laws capable to describe the material behavior under the extreme pressure as well as the low pressure are necessary for the analysis of the structural behavior under the intense -impulsive loadings. In this study, two plastic models, the pressure independent Von-Mises model and the pressure dependent Drucker-Prager model, are employed for the wave propagation analysis. Governing equations of this study are conservation equations of momentum and mass in Lagrangian coordinate system which is fixed to the material. Due to the shock-front which violates the continuity assumptions inherent in the differential equations numerical artificial viscosity is used to spread the shock front over several computational zones. These equations are solved by Finite Difference Method with discretized time and space coordinates. The associate normality flow rule as a plastic theory is implemented to find the plastic strains.

• Journal of the Korean Geotechnical Society
• /
• v.23 no.8
• /
• pp.137-147
• /
• 2007

Isotropic (instant loading) and Ko (gradual increase loading) consolidation tests were conducted in triaxial test equipment using cylindrical sample (5.0 cm in diameter and 10.0 cm in height) on two marine clay deposits. The duration of primary consolidation was estimated by two curve fitting methods using measured strain. A differential equation of consolidation for drainage in the radial and vertical direction was solved by the implicit finite difference scheme. The results of two curve fitting methods were compared with the numerical solutions to evaluate the appropriate axial loading rate of Ko consolidation and the primary consolidation periods. In addition, primary consolidation periods of the samples with a diameter of 35 mm and a height of 70 mm were calculated. The relation of radial and vertical consolidation coefficients is also presented.

• Geophysics and Geophysical Exploration
• /
• v.1 no.2
• /
• pp.127-135
• /
• 1998

The Born approximation is widely used for solving the complex scattering problems in electromagnetics. Approximating total internal electric field by the background field is reasonable for small material contrasts as long as scatterer is not too large and the frequency is not too high. However in many geophysical applications, moderate and high conductivity contrasts cause both real and imaginary part of internal electric field to differ greatly from background. In the extended Born approximation, which can improve the accuracy of Born approximation dramatically, the total electric field in the integral over the scattering volume is approximated by the background electric field projected to a depolarization tensor. The finite difference and elements methods are usually used in EM scattering problems with a 2D model and a 3D source, due to their capability for simulating complex subsurface conductivity distributions. The price paid for a 3D source is that many wavenumber domain solutions and their inverse Fourier transform must be computed. In these differential equation methods, all the area including homogeneous region should be discretized, which increases the number of nodes and matrix size. Therefore, the differential equation methods need a lot of computing time and large memory. In this study, EM modeling program for a 2D model and a 3D source is developed, which is based on the extended Born approximation. The solution is very fast and stable. Using the program, crosshole EM responses with a vertical magnetic dipole source are obtained and the results are compared with those of 3D integral equation solutions. The agreement between the integral equation solution and extended Born approximation is remarkable within the entire frequency range, but degrades with the increase of conductivity contrast between anomalous body and background medium. The extended Born approximation is accurate in the case conductivity contrast is lower than 1:10. Therefore, the location and conductivity of the anomalous body can be estimated effectively by the extended Born approximation although the quantitative estimate of conductivity is difficult for the case conductivity contrast is too high.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.9 no.3 s.43
• /
• pp.33-42
• /
• 2005

Frequency domain Volterra model is applied to nonlinear parameter identification procedure for dynamic systems modeled by nonlinear function. The frequency domain Volterra kernels, which correspond io linear, quadratic, and cubic transfer functions in lime domain, are incorporated in nonlinear parametric identification procedure. The nonlinear transfer functions, which can be derived from the Volterra series representation of the nonlinear differential equation of the system by Schetzen's method(1980), are directly used for modeling input output relation. The error is defined by the difference between the observed output and the estimated output which is calculated by substituting the observed input to nonlinear frequency domain model. The system parameters are searched by minimizing the error. Volterra model guarantees enough accuracy and convergence and the estimated coefficients have a good agreement with their actual values not only in the linear frequency region but also in the legion where the $2^{nd}\;or\;3^{rd}$ order nonlinearity is dominant.