• Journal of the Korean Physical Society
• /
• 제73권11호
• /
• pp.1774-1786
• /
• 2018

Phenotypic variation among clones (individuals with identical genes, i.e. isogenic individuals) has been recognized both theoretically and experimentally. We investigate the effects of phenotypic variation on evolutionary dynamics of a population. In a population, the individuals are assumed to be haploid with two genotypes : one genotype shows phenotypic variation and the other does not. We use an individual-based Moran model in which the individuals reproduce according to their fitness values and die at random. The evolutionary dynamics of an individual-based model is formulated in terms of a master equation and is approximated as the Fokker-Planck equation (FPE) and the coupled non-linear stochastic differential equations (SDEs) with multiplicative noise. We first analyze the deterministic part of the SDEs to obtain the fixed points and determine the stability of each fixed point. We find that there is a discrete phase transition in the population distribution when the probability of reproducing the fitter individual is equal to the critical value determined by the stability of the fixed points. Next, we take demographic stochasticity into account and analyze the FPE by eliminating the fast variable to reduce the coupled two-variable FPE to the single-variable FPE. We derive a quasi-stationary distribution of the reduced FPE and predict the fixation probabilities and the mean fixation times to absorbing states. We also carry out numerical simulations in the form of the Gillespie algorithm and find that the results of simulations are consistent with the analytic predictions.

• Structural Engineering and Mechanics
• /
• 제22권2호
• /
• pp.169-184
• /
• 2006

A close form solution of the maximum deflection for cracked columns with rectangular cross-sections was developed and thus the elastic buckling behavior and ultimate bearing capacity were studied analytically. First, taking into account the effect of the crack in the potential energy of elastic systems, a trigonometric series solution for the elastic deflection equation of an arbitrary crack position was derived by use of the Rayleigh-Ritz energy method and an analytical expression of the maximum deflection was obtained. By comparison with the rotational spring model (Okamura et al. 1969) and the equivalent stiffness method (Sinha et al. 2002), the advantages of the present solution are that there are few assumed conditions and the effect of axial compression on crack closure was considered. Second, based on the above solutions, the equilibrium paths of the elastic buckling were analytically described for cracked columns subjected to both axial and eccentric compressive load. Finally, as examples, the influence of crack depth, load eccentricity and column slenderness on the elastic buckling behavior was investigated in the case of a rectangular column with a single-edge crack. The relationship of the load capacity of the column with respect to crack depth and eccentricity or slenderness was also illustrated. The analytical and numerical results from the examples show that there are three kinds of collapse mechanisms for the various states of cracking, eccentricity and slenderness. These are the bifurcation for axial compression, the limit point instability for the condition of the deeper crack and lighter eccentricity and the fracture for higher eccentricity. As a result, the conception of critical transition eccentricity $(e/h)_c$, from limit-point buckling to fracture failure, was proposed and the critical values of $(e/h)_c$ were numerically determined for various eccentricities, crack depths and slenderness.

• 대한수학회지
• /
• 제57권6호
• /
• pp.1347-1372
• /
• 2020

We study bifurcation for the following fractional Schrödinger equation $$\{\left.\begin{eqnarray}(-{\Delta})^su+V(x)u&=&{\lambda}f(u)&&{\text{in}}\;{\Omega}\\u&>&0&&{\text{in}}\;{\Omega}\\u&=&0&&{\hspace{32}}{\text{in}}\;{\mathbb{R}}^n{\backslash}{\Omega}\end{eqnarray}\right$$ where 0 < s < 1, n > 2s, Ω is a bounded smooth domain of ℝⁿ, (-∆)^s is the fractional Laplacian of order s, V is the potential energy satisfying suitable assumptions and λ is a positive real parameter. The nonlinear term f is a positive nondecreasing convex function, asymptotically linear that is $\lim_{t{\rightarrow}+{\infty}}\;{\frac{f(t)}{t}}=a{\in}(0,+{\infty})$. We discuss the existence, uniqueness and stability of a positive solution and we also prove the existence of critical value and the uniqueness of extremal solutions. We take into account the types of Bifurcation problem for a class of quasilinear fractional Schrödinger equations, we also establish the asymptotic behavior of the solution around the bifurcation point.

• 대한수학회보
• /
• 제32권2호
• /
• pp.303-308
• /
• 1995

We consider a torus T, that is, a compact surface with genus 1 and $\Omega = D^2 \times S^1$ topologically with $\partial\Omega = T$, where $D^2$ is the open unit disk and $S^1$ is the unit circle. Let $\omega = (x,y)$ denote the generic point on T. For a smooth immersion $u : T \to R^3$, we define the Dirichlet functional by $$ E(u) = \frac{2}{1} \int_{T} $\mid$\nabla u$\mid$^2 d\omega $$ and the volume functional by $$ V(u) = \frac{3}{1} \int_{T} u \cdot u_x \Lambda u_y d\omege $$.

• 대한수학회보
• /
• 제32권2호
• /
• pp.211-214
• /
• 1995

Let M and N be compact Riemannian manifolds and $f : M \to N$ be a smooth map. Following J. Eells, f is exponentially harmonic if it represents a critical point of the exponential energy integral $$ E(f) = \int_{M} exp(\left\$\mid$ df \right\$\mid$^2) dM $$ where $(\left\ df $\mid$\right\$\mid$^2$ is the energy density defined as $\sum_{i=1}^{m} \left\$\mid$ df(e_i) \right\$\mid$^2$, m = dimM, for orthonormal frame $e_i$ of M. The Euler- Lagrange equation of the exponential energy functional E can be written $$ exp(\left\$\mid$ df \right\$\mid$^2)(\tau(f) + df(\nabla\left\$\mid$ df \right\$\mid$^2)) = 0 $$ where $\tau(f)$ is the tension field along f. Hence, if the energy density is constant, every harmonic map is exponentially harmonic and vice versa.

• 한국공간구조학회:학술대회논문집
• /
• 한국공간구조학회 2008년도 춘계 학술발표회 논문집
• /
• pp.89-94
• /
• 2008

This paper deals with the method of self-equilibrium stress mode analysis of cable dome structures. From the point of view of analysis, cable dome structure is a kind of unstable truss structure which is stabilized by means of introduction of prestressing. The prestress must be introduced according to a specific proportion among different structural member and it is determined by an analysis called self-equilibrium stress mode analysis. The mathematical equation involved in the self-equilibrium stress mode analysis is a system of linear equations which can be solved numerically by adopting the concept of Moore-Penrose generalized inverse. The calculation of the generalized inverse is carried out by rank factorization method. This method involves a parameter called epsilon which plays a critical role in self-equilibrium stress mode analysis. It is thus of interest to investigate the range of epsilon which produces consistent solution during the analysis of self-equilibrium stress mode.

• 한국유체기계학회 논문집
• /
• 제7권4호
• /
• pp.32-39
• /
• 2004

This study deals with the Benard Flow of Magnetic Fluids in a rectangular cavity which the ratio between height and width is 1 : 4 and the base side or left side is a heating face while other sides are to be cooling faces. When Magnetic field was equally impressed, considering the internal rotation of the elementary ferromagnetic particle, we found the following result from the numerical analysis of the GSMAC algorithm applied to the equation of the magnetic fluid. Benard flow is controlled by intensity and direction of magnetic fields, and critical point appears when especially magnetic field with a heating base and side area near H=-7000 and H=-10000 is applied.

• 산업기술연구
• /
• 제18권
• /
• pp.61-67
• /
• 1998

이 논문에서는 Diruchlet 경계 조건을 갖는 비선형 타원형 방정식 $-{\Delta}u+g(u)=f(x)$의 해의 존재에 대한 연구를 하였다. 존재하는 해의 다중성을 증명하기 위하여 임계점 이론과 롤의 정리를 사용하였으며, 대응되는 범함수에 따라서 방정식의 해와 임계점이 동시에 나타난다는 정리를 이용하였다. 이 때 $g(u)=bu^+-au^-$으로 나타날 때 외력항 (방정식의 우변)의 상수로 주어지는 경우 적어도 두 개의 해가 존재한다는 것을 증명하였다. 만약 우변(외력항)의 상수가 음수이거나 0인 경우이 방정식의 해가 존재하지 않거나 자명한 해만 존재하기 때문에 상수는 양수인 것으로 가정하였다.

• Korean Chemical Engineering Research
• /
• 제54권2호
• /
• pp.206-212
• /
• 2016

Experimental data are reported on the phase equilibrium of propoxylated neopentyl glycol diacrylate in supercritical carbon dioxide. Phase equilibria data were measured in static method at a temperature of (313.2, 333.2, 353.2, 373.2 and 393.2) K and at pressures up to 27.82 MPa. At a constant pressure, the solubility of propoxylated neopentyl glycol diacrylate for the (carbon dioxide + propoxylated neopentyl glycol diacrylate) system increases as temperature increases. The (carbon dioxide + propoxylated neopentyl glycol diacrylate) system exhibits type-I phase behavior. The experimental result for the (carbon dioxide + propoxylated neopentyl glycol diacrylate) system is correlated with Peng-Robinson equation of state using mixing rule. The critical property of propoxylated neopentyl glycol diacrylate is predicted with Joback and Lyderson method.

• 한국소음진동공학회논문집
• /
• 제16권6호
• /
• pp.634-642
• /
• 2006

For the analysis of dynamic contact between a catenary and a pantograph of high-speed electrical train, the numerical solution of the equations of motion of the vehicle pantograph and the catenary system subjected to the contact condition is obtained. The whole equations of motion of the catenary and the pantograph are simultaneously time integrated with the strict application of the contact condition. For the stability of the numerical solution, with the cubic spline interpolation of the catenary displacement, the velocity and acceleration constraints as well as the displacement constraint are imposed on the contact point. Especially it is shown that the Coriolis and centripetal accelerations are critical for the accuracy and stability of the computation.