• Kyungpook Mathematical Journal
• /
• v.63 no.2
• /
• pp.141-154
• /
• 2023

We study the stability of traveling waves of a certain chemotaxis model. The traveling wave solution is a central object of study in a chemotaxis model. Kim et al. [8] introduced a model on a population and nutrient densities based on a nonlinear diffusion law. They proved the existence of traveling waves for the one dimensional Cauchy problem. Existence theory for traveling waves is typically followed by stability analysis because any traveling waves that are not robust against a small perturbation would have little physical significance. We conduct a numerical nonlinear stability for a few relevant instances of traveling waves shown to exist in [8]. Results against absolute additive noises and relative additive noises are presented.

• Bulletin of the Korean Mathematical Society
• /
• v.61 no.3
• /
• pp.637-669
• /
• 2024

The goal of this paper is to study decay properties of weak solutions to Cauchy problem of the non-stationary fractional Navier-Stokes equations. By using the Fourier splitting method, we give the time L²-decay rate of weak solutions, which reveals that L²-decay is generally determined by its linear generalized Stokes flow. In second part, we establish various decay results and the uniqueness of the two dimensional fractional Navier-Stokes flows. In the end of this article, as an appendix, the existence of global weak solutions is given by making use of Galerkin' method, weak and strong compact convergence theorems.

• Structural Engineering and Mechanics
• /
• v.67 no.3
• /
• pp.277-281
• /
• 2018

The analytical solution of two functionally graded layers with Volterra type screw dislocation is investigated under anti-plane shear impact loading. The energy dissipation of FGM layers is modeled by viscous damping and the properties of the materials are assumed to change exponentially along the thickness of the layers. In this study, the rate of gradual change ofshear moduli, mass density and damping constant are assumed to be same. At first, the stress fields in the interface of the FGM layers are derived by using a single dislocation. Then, by determining a distributed dislocation density on the crack surface and by using the Fourier and Laplace integral transforms, the problem are reduce to a system ofsingular integral equations with simple Cauchy kernel. The dynamic stress intensity factors are determined by numerical Laplace inversion and the distributed dislocation technique. Finally, various examples are provided to investigate the effects of the geometrical parameters, material properties, viscous damping and cracks configuration on the dynamic fracture behavior of the interacting cracks.

• Bulletin of the Korean Mathematical Society
• /
• v.57 no.3
• /
• pp.649-660
• /
• 2020

This paper proves a new regularity criterion for solutions to the Cauchy problem of the 3D Boussinesq equations via one directional derivative of the horizontal component of the velocity field (i.e., (∂_iu₁; ∂_ju₂; 0) where i, j ∈ {1, 2, 3}) in the framework of the anisotropic Lebesgue spaces. More precisely, for 0 < T < ∞, if $$\large{\normalsize\displaystyle\smashmargin{2}{\int\nolimits_o}^T}({\HUGE\left\|{\small{\parallel}{\partial}_iu_1(t){\parallel}_{L^{\alpha}_{x_i}}}\right\|}{\small^{\gamma}_{L^{\beta}_{x_{\hat{i}}x_{\bar{i}}}}+}{\HUGE\left\|{\small{\parallel}{\partial}_iu_2(t){\parallel}_{L^{\alpha}_{x_j}}}\right\|}{\small^{\gamma}_{L^{\beta}_{x_{\hat{i}}x_{\bar{i}}}}})dt<{{\infty}},$$ where ${\frac{2}{{\gamma}}}+{\frac{1}{{\alpha}}}+{\frac{2}{{\beta}}}=m{\in}[1,{\frac{3}{2}})$ and ${\frac{3}{m}}{\leq}{\alpha}{\leq}{\beta}<{\frac{1}{m-1}}$, then the corresponding solution (u, θ) to the 3D Boussinesq equations is regular on [0, T]. Here, (i, ${\hat{i}}$, ${\tilde{i}}$) and (j, ${\hat{j}}$, ${\tilde{j}}$) belong to the permutation group on the set 𝕊₃ := {1, 2, 3}. This result reveals that the horizontal component of the velocity field plays a dominant role in regularity theory of the Boussinesq equations.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.10 no.5
• /
• pp.787-797
• /
• 1986

An implicit approach is employed to obtain a general boundary integral formulation of axisymmetric elastic problems in terms of a pair of singular integral equations. The corresponding kernel functions from the solutions of Navier's equation are derived by applying a three dimensional integral and a direct axisymmetrical approach. A numerical discretization schem including the evaluation of Cauchy principal values of the singular integral is described. Finally the typical axisymmetric elastic models are analyzed, i.e. the hollow sphere, the constant thickness and the V-notched round bar.

• Structural Engineering and Mechanics
• /
• v.54 no.3
• /
• pp.433-451
• /
• 2015

This paper focuses on large deflection static behavior of edge cracked simple supported beams subjected to a non-follower transversal point load at the midpoint of the beam by using the total Lagrangian Timoshenko beam element approximation. The cross section of the beam is circular. The cracked beam is modeled as an assembly of two sub-beams connected through a massless elastic rotational spring. It is known that large deflection problems are geometrically nonlinear problems. The considered highly nonlinear problem is solved considering full geometric non-linearity by using incremental displacement-based finite element method in conjunction with Newton-Raphson iteration method. There is no restriction on the magnitudes of deflections and rotations in contradistinction to von-Karman strain displacement relations of the beam. The beams considered in numerical examples are made of Aluminum. In the study, the effects of the location of crack and the depth of the crack on the non-linear static response of the beam are investigated in detail. The relationships between deflections, end rotational angles, end constraint forces, deflection configuration, Cauchy stresses of the edge-cracked beams and load rising are illustrated in detail in nonlinear case. Also, the difference between the geometrically linear and nonlinear analysis of edge-cracked beam is investigated in detail.

• Structural Engineering and Mechanics
• /
• v.66 no.3
• /
• pp.329-341
• /
• 2018

This paper considers the smooth receding contact problem between a homogeneous half-plane and a composite laminate composed of an inhomogeneously coated elastic layer. The inhomogeneity of the elastic modulus of the coating is approximated by an exponential function along the thickness dimension. The three-component structure is pressed together by either a concentrated force or uniform pressures applied at the top surface of the composite laminate. Both semianalytical and finite element analysis are performed to solve for the extent of contact and the contact pressure. In the semianalytical formulation, Fourier integral transformation of governing equations and boundary conditions leads to a singular integral equation of Cauchy-type, which can be numerically integrated by Gauss-Chebyshev quadrature to a desired degree of accuracy. In the finite element modeling, the functionally graded coating is divided into homogeneous sublayers and the shear modulus of each sublayer is assigned at its lower boundary following the predefined exponential variation. In postprocessing, the stresses of any node belonging to sublayer interfaces are averaged over its surrounding elements. The results obtained from the semianalytical analysis are successfully validated against literature results and those of the finite element modeling. Extensive parametric studies suggest the practicability of optimizing the receding contact peak stress and the extent of contact in multilayered structures by the introduction of functionally graded coatings.

• Structural Engineering and Mechanics
• /
• v.73 no.5
• /
• pp.585-598
• /
• 2020

We develop design model within which nucleation and propagation of crack in a fibrous composite is described. It is assumed that under loading, crack initiation and fracture of material happens in the composite. The problem of equilibrium of a composite with embryonic crack is reduced to the solution of the system of nonlinear singular integral equations with the Cauchy type kernel. Normal and tangential forces in the crack nucleation zone are determined from the solution of this system of equations. The crack appearance conditions in the composite are formed with regard to criterion of ultimate stretching of the material's bonds. We study the case when near the fiber, the binder has several arbitrary arranged rectilinear prefracture zones and a crack with interfacial bonds. The proposed computational model allows one to obtain the size and location of the zones of damages (prefracture zones) depending on geometric and mechanical characteristics of the fibrous composite and applied external load. Based on the suggested design model that takes into account the existence of damages (the zones of weakened interparticle bonds of the material) and cracks with end zones in the composite, we worked out a method for calculating the parameters of the composite, at which crack nucleation and crack growth occurs.

• Structural Engineering and Mechanics
• /
• v.86 no.4
• /
• pp.535-546
• /
• 2023

In the present paper, multiple interface cracks between a functionally graded orthotropic coating and an orthotropic half-plane substrate under concentrated loading are considered by means of the distribution dislocation technique (DDT). With the use of integration of Fourier transform the problem is reduced to a system of Cauchy-type singular integral equations which are solved numerically to compute the dislocation density on the surfaces of the cracks. The distribution dislocation is a powerful method to calculate accurate solutions to plane crack problems, especially this method is very good to find SIFs for multiple unequal cracks located at the interface. Hence this technique allows considering any number of interface cracks. The primary objective of this paper is to investigate the effects of the interaction of multiple interface cracks, load location, material orthotropy, nonhomogeneity parameters and geometry parameters on the modes I and II SIFs. Numerical results show that modes I/II SIFs decrease with increasing the nonhomogeneity parameter and the highest magnitude of SIF occurs where distances between the load location and crack tips are minimal.

• Journal of the Society of Naval Architects of Korea
• /
• v.29 no.4
• /
• pp.114-131
• /
• 1992

In the present work, a two-dimensional boundary-value problem for a large amplitude motion is treated as an initial-value problem by satisfying the exact body-boundary and nonlinear free-surface boundary conditions. The present nonlinear numerical scheme is similar to that described by Vinje and Brevig(1981) who utilized the Cauchy's theorem and assumed the periodicity in the horizontal coordinate. In the present thesis, however, the periodicity in the horizontal coordinate is not assumed. Thus the present method can treat more realistic problems, which allow radiating waves to infinities. In the present method of solution, the original infinite fluid domain, is divided into two subdomains ; ie the inner and outer subdomains which are a local nonlinear subdomain and the truncated infinite linear subdomain, respectively. By imposing an appropriate matching condition, the computation is carried out only in the inner domain which includes the body. Here we adopt the nonlinear scheme of Vinje & Brevig only in the inner domain and respresent the solution in the truncated infinite subdomains by distributing the time-dependent Green function on the matching boundaries. The matching condition is that the velocity potential and stream function are required to be continuous across the matching boundary. In the computations we used, if necessary, a regriding algorithm on the free surface which could give converged stable solutions successfully even for the breaking waves. In harmonic oscillation problem, each harmonic component and time-mean force are obtained by the Fourier transform of the computed forces in the time domain. The numerical calculations are made for the following problems. $\cdot$ Forced harmonic large-amplitude oscillation(${\omega}{\neq}0,\;U=0$) $\cdot$ Translation with a uniform speed(${\omega}=0,\;U{\neq}0$) The computed results are compared with available experimental data and other analytical results.