• Bulletin of the Korean Chemical Society
• /
• v.35 no.3
• /
• pp.858-864
• /
• 2014

Recently, many researches have shown that even photosynthetic light-harvesting pigment-protein complexes can have quantum coherence in their excitonic energy transfer at cryogenic and physiological temperatures. Because the protein supplies such noisy environment around pigments that conventional wisdom expects very short lived quantum coherence, elucidating the mechanism and searching for an applicability of the coherence have become an interesting topic in both experiment and theory. We have previously studied the quantum coherence of a phycocyanin 645 complex in a marine algae harvesting light system, using Poisson mapping bracket equation (PBME). PBME is one of the applicable methods for solving quantum-classical Liouville equation, for following the dynamics of such pigment-protein complexes. However, it may suffer from many defects mostly from mapping quantum degrees of freedom into classical ones. To make improvements against such defects, benchmarking targets with more accurately described dynamics is highly needed. Here, we fall back to reduced hierarchical equation of motion (HEOM), for such a purpose. Even though HEOM is known to applicable only to simplified system that is coupled to a set of harmonic oscillators, it can provide ultimate accuracy within the regime of quantum-classical description, thus providing perfect benchmark targets for certain systems. We compare the evolution of the density matrix of pigment excited states by HEOM against the PBME results at physiological temperature, and observe more sophisticated changes of density matrix elements from HEOM. In PBME, the population of states with intermediate energies display only monotonically increasing behaviors. Most importantly, PBME suffers a serious issue of wrong population in the long time limit, likely generated by the zero-point energy leaking problem. Future prospects for developments are briefly discussed as a concluding remark.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.3
• /
• pp.380-391
• /
• 1997

A numerical study has been done on the evolution of particle size distribution in particle laden high temperature jet flows undergoing convection, diffusion, thermophoresis and coagulation. The dynamic behavior of these particles have been modelled by approximating the particle size distribution by a lognormal function throughout the process and the moments of the particle size distribution have been used to solve the general dynamic equation. The size distributions of spherical particles in the radial and axial direction have been obtained including the effect of buoyancy. Of particular interests are the variations of geometric mean diameter, number concentration and polydispersity. Results show that buoyancy significantly alters the size distribution in both axial and radial direction. One dimensional analysis for non-spherical particles has also been done and the results have been compared with the existing experimental data.

• Journal of Welding and Joining
• /
• v.25 no.1
• /
• pp.63-69
• /
• 2007

Once assessment of material failure characteristics is captured precisely in a unified way, it can bedirectly incorporated into the structural failure assessment under various loading environments, based on the theoretical backgrounds so called Local Approach to Fracture. The aim of this study is to develop a numerical fatigue test method by continuum damage mechanics applicable for the assessment of structural integrity throughout crack initiation and structural failure based on the Local Approach to Fracture. The generalized elasto-visco-plastic constitutive equation, which can consider the internal damage evolution behavior, is developed and employed in the 3-D FEA code in order to numerically evaluate the material and/or structural responses. Explicit information of the relationships between the mechanical properties and material constants, which are required for the mechanical constitutive and damage evolution equations for each material, are implemented in numerical fatigue test method. The material constants selected from constitutive equations are used directly in the failure assessment of material and/or structures. The performance of the developed system has been evaluated with assessing the S-N diagram of stainless steel materials.

• Interaction and multiscale mechanics
• /
• v.2 no.4
• /
• pp.375-396
• /
• 2009

In the present work, the elasto-viscoplastic behavior, interactions between grains, and the texture evolution in polycrystalline materials subjected to finite deformations are modeled using a multiscale analysis procedure within a finite element framework. Computational homogenization is used to relate the grain (meso) scale to the macroscale. Specifically, a polycrystal is modeled by a material representative volume element (RVE) consisting of an aggregate of grains, and a periodic distribution of such unit cells is considered to describe material behavior locally on the macroscale. The elastic behavior is defined by a hyperelastic potential, and the viscoplastic response is modeled by a simple power law complemented by a work hardening equation. The finite element framework is based on a Lagrangian formulation, where a kinematic split of the deformation gradient into volume preserving and volumetric parts together with a three-field form of the Hu-Washizu variational principle is adopted to create a stable finite element method. Examples involving simple deformations of an aluminum alloy are modeled to predict inhomogeneous fields on the grain scale, and the macroscopic effective stress-strain curve and texture evolution are compared to those obtained using both upper and lower bound models.

• Korea-Australia Rheology Journal
• /
• v.21 no.3
• /
• pp.193-200
• /
• 2009

We present a direct simulation technique for two-dimensional mold-filling simulations of fluids filled with a large number of circular disk-like rigid particles. It is a direct simulation in that the hydrodynamic interaction between particles and fluid is fully considered. We employ a pseudo-concentration method for the evolution of the flow front and the DLM (distributed Lagrangian multipliers)-like fictitious domain method for the implicit treatment of the hydrodynamic interaction. Both methods allow the use of a fixed regular discretization during the entire computation. The discontinuous Galerkin method has been used to solve the concentration evolution equation and the rigid-ring description has been introduced for freely suspended particles. A buffer zone, the gate region of a finite area subject to the uniform velocity profile, has been introduced to put discrete particles into the computational domain avoiding any artificial discontinuity. From example problems of 450 particles, we investigated the particle motion and effects of particles on the flow for both Newtonian and shear-thinning fluid media. We report the prolonged particle movement toward the wall in case of a shear-thinning fluid, which has been interpreted with the shear rate distribution.

• Transactions of Materials Processing
• /
• v.9 no.1
• /
• pp.72-79
• /
• 2000

Dynamic recrystallization (DRX), which may occur during hot deformation, is important for the microsturctural evolution of 304 stainless steel. Especially, the current interest in modelling hot rolling demands quantitative relationships among the thermomechanical process variables, such as strain, temperature, strain rate, and etc. Thus, this paper individually presents the relationships for flow stress and volume fraction of DRX as a function of processing variables using torsion tests. The hot torsion tests of 304 stainless steel were performed at the temperature range of 900~110$0^{\circ}C$ and the strain rate range of 5x10-2~5s-1 to study the high temperature softening behavior. For the exact prediction of flow stress, the equation was divided into two regions, the work hardening (WH) and dynamic recovery (DRV) region and the DRX region. Especially, The flow stress of DRX region could be expressed by using the volume fraction of DRX (XDRX). Since XDRX was consisted of the critical strain($\varepsilon$c) for initiation of dynamic recrystallization (DRX) and the strain for maximum softening rate ($\varepsilon$*), that were related with the evolution of microstructure. The calculated results predicted the flow stress and the microstructure of the alloy at any deformation conditions well.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.37 no.5A
• /
• pp.311-316
• /
• 2012

In this paper, a bit-mapping scheme is proposed for partial Chase combining in LDPC-coded systems. Contrary to the previously known bit mapping that assigns the information bits to more reliable channels, the proposed mapping assigns the codeword bits of irregular LDPC codes to distinct Gaussian channels by considering the characteristics of LDPC codes and channels. The recursion equation for partial Chase combining is derived by using the density evolution technique, based on it, the best bit mapping among the various bit-mapping schemes is derived, and the validity of them is confirmed through simulation.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.7 no.3
• /
• pp.595-609
• /
• 2015

The wake characteristics of Contra-Rotating Propeller (CRP) were investigated using numerical simulation and flow measurement. The numerical simulation was carried out with a commercial CFD code based on a Reynolds Averaged Navier-Stokes (RANS) equations solver, and the flow measurement was performed with Stereoscopic Particle Image Velocimetry (SPIV) system. The simulation results were validated through the comparison with the experiment results measured around the leading edge of rudder to investigate the effect of propeller operation under the conditions without propeller, with forward propeller alone, and with both forward and aft propellers. The evolution of CRP wake was analyzed through velocity and vorticity contours on three transverse planes and one longitudinal plane based on CFD results. The trajectories of propeller tip vortex core in the cases with and without aft propeller were also compared, and larger wake contraction with CRP was confirmed.

• Geomechanics and Engineering
• /
• v.16 no.6
• /
• pp.619-626
• /
• 2018

A simple prediction procedure was investigated for calculating the stresses and displacements of a circular opening. Unlike existed approaches, the proposed approach starts each step with a radius increment. The stress for each annulus could be obtained analytically, while strain increments for each step can be determinate numerically from the compatility equation by finite difference approximation, flow rule and Hooke's law. In the successive manner, the distributions of stresses and displacements could be found. It should be noted that the finial radial stress and displacement were equal to the internal supporting pressure and deformation at the tunnel wall, respectively. By assuming different plastic radii, GRC and the evolution curve of plastic radii and internal supporting pressures could be obtained conveniently. Then the real plastic radius can be calculated by using linear interpolation in the evolution curve. Some numerical and engineering examples were performed to demonstrate the accuracy and validity for the proposed procedure. The comparisons results show that the proposed procedure was faster than that in Lee and Pietrucszczak (2008). The influence of annulus number and dilation on the accuracy of solutions was also investigated. Results show that the larger the annulus number was, the more accurate the solutions were. Solutions in Park et al. (2008) were significantly influenced by dilation.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.11 no.1
• /
• pp.177-188
• /
• 1991

The littoral drift in the nearshore zone occurs mainly due to wave action and nearshore currents. Beach topotraphical evolution model presented in this study shows the possibility which can be applied to the prediction of beach deformation over short time interval, associated with the construction of coastal structures. The 3-dimensional beach evolution model, based on Watanabe's equation and Deguchi's flux model, is developed and consists of three submodels of wave transformation, rip-currents, and beach deformation. This model is applied to the several cases with different conditions and compared with the results of Watanabe's numerical model. In addition, the effects of parameters involved are discussed.