• 천문학회지
• /
• 제34권4호
• /
• pp.281-283
• /
• 2001

As a companion to an adiabatic version developed by Ryu and his coworkers, we have built an isothermal magnetohydrodynamic code for astrophysical flows. It is suited for the dynamical simulations of flows where cooling timescale is much shorter than dynamical timescale, as well as for turbulence and dynamo simulations in which detailed energetics are unimportant. Since a simple isothermal equation of state substitutes the energy conservation equation, the numerical schemes for isothermal flows are simpler (no contact discontinuity) than those for adiabatic flows and the resulting code is faster. Tests for shock tubes and Alfven wave decay have shown that our isothermal code has not only a good shock capturing ability, but also numerical dissipation smaller than its adiabatic analogue. As a real astrophysical application of the code, we have simulated the nonlinear three-dimensional evolution of the Parker instability. A factor of two enhancement in vertical column density has been achieved at most, and the main structures formed are sheet-like and aligned with the mean field direction. We conclude that the Parker instability alone is not a viable formation mechanism of the giant molecular clouds.

• Bulletin of the Korean Chemical Society
• /
• 제14권1호
• /
• pp.21-29
• /
• 1993

For the study of relaxation processes in complex spin system, a general master equation, which can be used to simulate a vast range of pulse experiments, has been formulated using the Liouville representation of quantum mechanics. The state of a nonequilibrium spin system in magnetic field is described by a density vector in Liouville space and the time evolution of the system is followed by the application of a linear master operator to the density vector in this Liouville space. In this master equation the nuclear spin relaxation due to intramolecular dipolar interaction or randomly fluctuating field interaction is explicitly implemented as a relaxation supermatrix for a strong coupled two-spin (1/2) system. The whole dynamic information inherent in the spin system is thus contained in the density vector and the master operator. The radiofrequency pulses are applied in the same space by corresponding unitary rotational supertransformations of the density vector. If the resulting FID is analytically Fourier transformed, it is possible to represent the final nonstationary spectrum using a frequency dependent spectral vector and intensity determining shape vector. The overall algorithm including relaxation interactions is then translated into an ANSIFORTRAN computer program, which can simulate a variety of two dimensional spectra. Furthermore a new strategy is tested by simulation of multiple quantum signals to differentiate the two relaxation interaction types.

• Bulletin of the Korean Chemical Society
• /
• 제33권3호
• /
• pp.933-940
• /
• 2012

Recently, it has been shown that quantum coherence appears in energy transfers of various photosynthetic lightharvesting complexes at from cryogenic to even room temperatures. Because the photosynthetic systems are inherently complex, these findings have subsequently interested many researchers in the field of both experiment and theory. From the theoretical part, simplified dynamics or semiclassical approaches have been widely used. In these approaches, the quantum-classical Liouville equation (QCLE) is the fundamental starting point. Toward the semiclassical scheme, approximations are needed to simplify the equations of motion of various degrees of freedom. Here, we have adopted the Poisson bracket mapping equation (PBME) as an approximate form of QCLE and applied it to find the time evolution of the excitation in a photosynthetic complex from marine algae. The benefit of using PBME is its similarity to conventional Hamiltonian dynamics. Through this, we confirmed the coherent population transfer behaviors in short time domain as previously reported with a more accurate but more time-consuming iterative linearized density matrix approach. However, we find that the site populations do not behave according to the Boltzmann law in the long time limit. We also test the effect of adding spurious high frequency vibrations to the spectral density of the bath, and find that their existence does not alter the dynamics to any significant extent as long as the associated reorganization energy is changed not too drastically. This suggests that adopting classical trajectory based ensembles in semiclassical simulations should not influence the coherence dynamics in any practical manner, even though the classical trajectories often yield spurious high frequency vibrational features in the spectral density.

• Geomechanics and Engineering
• /
• 제7권5호
• /
• pp.569-587
• /
• 2014

Triaxial compression creep tests were performed on salt rock samples using cyclic confining pressure with a static axial pressure. The test results show that, up to a certain time, changes in the confining pressure have little influence on creep properties of salt rock, and the axial creep curve is smooth. After this point, the axial creep curve clearly fluctuates with the confining pressure, and is approximately a straight line both when the confining pressure decreases and when it increases within one cycle period. The slope of these lines differs: it is greater when the confining pressure decreases than when it increases. In accordance with rheology model theory, axial creep equations were deduced for Maxwell and Kelvin models under cyclic loading. These were combined to establish an axial creep equation for the Burgers model. We supposed that damage evolution follows an exponential law during creep process and replaced the apparent stress in creep equation for the Burgers model with the effective stress, the axial creep damage equation for the Burgers model was obtained. The model suitability was verified using creep test results for salt rock. The fitting curves are in excellent agreement with the test curves, so the proposed model can well reflect the creep behavior of salt rock under low-frequency cyclic loading. In particular, it reflects the fluctuations in creep deformation and creep rate as the confining pressure increasing and decreasing under different cycle periods.

• KIEE International Transactions on Electrophysics and Applications
• /
• 제4C권4호
• /
• pp.155-159
• /
• 2004

This work reports the theoretical and experimental investigations of capillary bust valves to regulate liquid flow in microchannels. The theoretical analysis uses the Young-Laplace equation and geometrical considerations to predict the pressure at the edge of the valve opening. Numerical simulations are employed to calculate the meniscus shape evolution while the interface is pinned at the valve edge. Microchannels and valves are fabricated using soft lithography. A wafer-rotating system, which can adjust the driving pressure by rotational speed, induces a liquid flow. Experimentally measured valve-bursting pressure agrees with theoretical predictions.

• KSII Transactions on Internet and Information Systems (TIIS)
• /
• 제8권3호
• /
• pp.735-748
• /
• 2014

The fault-tolerance routing problem is one of the most important issues in the application of the Internet of Things, and has been attracting growing research interests. In order to maintain the communication paths from source sensors to the macronodes, we present a hybrid routing scheme and model, in which alternate paths are created once the previous routing is broken. Then, we propose an improved efficient and intelligent fault-tolerance algorithm (IEIFTA) to provide the fast routing recovery and reconstruct the network topology for path failure in the Internet of Things. In the IEIFTA, mutation direction of the particle is determined by multi-swarm evolution equation, and its diversity is improved by the immune mechanism, which can improve the ability of global search and improve the converging rate of the algorithm. The simulation results indicate that the IEIFTA-based fault-tolerance algorithm outperforms the EARQ algorithm and the SPSOA algorithm due to its ability of fast routing recovery mechanism and prolonging the lifetime of the Internet of Things.

• International Journal of Aeronautical and Space Sciences
• /
• 제16권2호
• /
• pp.165-176
• /
• 2015

A numerical analysis was conducted to investigate the inlet buzz and combustion oscillation in an axisymmetric ramjet engine with wedge-type flame holders. The physical model of concern includes the entire engine flow path, extending from the leading edge of the inlet center-body through the exhaust nozzle. The theoretical formulation is based on the Farve-averaged conservation equations of mass, momentum, energy, and species concentration, and accommodates finite-rate chemical kinetics and variable thermo-physical properties. Turbulence closure is achieved using a combined scheme comprising of a low-Reynolds number k-${\varepsilon}$ two-equation model and Sarkar's compressible turbulence model. Detailed flow phenomena such as inlet flow aerodynamics, flame evolution, and acoustic excitation as well as their interactions, are investigated. Mechanisms responsible for driving the inlet buzz are identified and quantified for the engine operating at subcritical conditions.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2009년 추계학술대회논문집
• /
• pp.109-114
• /
• 2009

As one of the promising model on the multiphase fluid mixtures, the Lattice-Boltzmann Method(LBM) is being developed to simulate flows containing two immisible components which are different mass values. The equilibrium function in the LBM can have a nonideal gas model for the equation of state and use the interfacial energy for the phase separation effect. An example on the phase separation has been carried out through the time evolution. The LBM based on the statistic mechanics is appropriate to solve very complicated flow problems and this model gives comparative merits rather than the continuum mechanics model.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2009년 추계학술대회논문집
• /
• pp.223-227
• /
• 2009

In the present study, a least square weighted residual method and Taylor-Galerkin method were formulated and tested for the discretization of the two hyperbolic type equations of level set method; advection and reinitialization equations. The two approaches were compared by solving a time reversed vortex flow and three-dimensional broken dam flow by employing a four-step splitting finite element method for the solution of the incompressible Navier-Stokes equations. From the numerical experiments, it was shown that the least square method is more accurate and conservative than Taylor-Galerkin method and both methods are approximately first order accurate when both advection and reinitialization phase are involved in the evolution of free surface.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2004년도 추계학술대회
• /
• pp.465-467
• /
• 2004

Finite element analysis of welding processes, which entail phase evolution, heat transfer and deformation, is considered in this paper. Attention focuses on numerical implementation of the thermo-elastic-plastic constitutive equation proposed by Leblond et al in consideration of the transformation plasticity. Based upon the multiplicative decomposition of deformation gradient, hyperelastic formulation is employed for efficient numerical integration, and the algorithmic consistent moduli for elastic-plastic deformations including transformation plasticity are obtained in the closed form. The convergence behavior of the present implementation is demonstrated via a couple of numerical examples. Several locking phenomena removed by Solid-shell element.