• Journal of Pharmaceutical Investigation
• /
• v.29 no.4
• /
• pp.295-307
• /
• 1999

Using a Rheometries Fluids Spectrometer (RFS II), the dynamic viscoelastic properties of aqueous poly(ethylene oxide) (PEO) solutions in small amplitude oscillatory shear flow fields have been measured over a wide range of angular frequencies. The angular frequency dependence of the storage and loss moduli at various molecular weights and concentrations was reported in detail, and the result was interpreted using the concept of a Deborah number De. In addition, the experimentally determined critical angular frequency at which the storage and loss moduli become equivalent was compared with the calculated characteristic time (or its inverse value), and their physical significance in analyzing the dynamic viscoelastic behavior was discussed. Finally, the relationship between steady shear flow and dynamic viscoelstic properties was examined by evaluating the applicability of some proposed models that describe the correlations between steady flow viscosity and dynamic viscosity, dynamic fluidity, and complex viscosity. Main results obtained from this study can be summarized as follows: (1) At lower angular frequencies where De<1, the loss modulus is larger than the storage modulus. However, such a relation between the two moduli is reversed at higher angular frequencies where De>l, indicating that the elastic behavior becomes dominant to the viscous behavior at frequency range higher than a critical angular frequency. (2) A critical angular frequency is decreased as an increase in concentration and/or molecular weight. Both the viscous and elastic properties show a stronger dependence on the molecular weight than on the concentration. (3) A characteristic time is increased with increasing concentration and/or molecular weight. The power-law relationship holds between the inverse value of a characteristic time and a critical angular frequency. (4) Among the previously proposed models, the Cox-Merz rule implying the equivalence between the steady flow viscosity and the magnitude of the complex viscosity has the best validity. The Osaki relation can be regarded to some extent as a suitable model. However, the DeWitt, Pao and HusebyBlyler models are not applicable to describe the correlations between steady shear flow and dynamic viscoelastic properties.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.26 no.5
• /
• pp.393-399
• /
• 2013

In this paper, a novel approach to estimate cohesive laws for the mode I fracture of the graphene is presented by combining molecular dynamic simulations and an inverse algorithm based on field projection method and finite element method. The determination of crack-tip cohesive laws of the graphene based on continuum mechanics is a non-trivial inverse problem of finding unknown tractions and separations from atomic simulations. The displacements of molecular dynamic simulations in a region far away from the crack tip are transferred to finite element nodes by using moving least square approximation. Inverse analyses for extracting unknown cohesive tractions and separation behind the crack tip can be carried out by using conservation nature of the interaction J- and M-integrals with numerical auxiliary fields which are generated by systematically imposing uniform surface tractions element-by-element along the crack surfaces in finite element models. The preset method can be a very successful approach to extract crack-tip cohesive laws from molecular dynamic simulations as a scale bridging method.

• Steel and Composite Structures
• /
• v.34 no.2
• /
• pp.261-277
• /
• 2020

The main objective of this research paper is to consider vibration analysis of vacancy defected graphene sheet as a nonisotropic structure via molecular dynamic and continuum approaches. The influence of structural defects on the vibration of graphene sheets is considered by applying the mechanical properties of defected graphene sheets. Molecular dynamic simulations have been performed to estimate the mechanical properties of graphene as a nonisotropic structure with single- and double- vacancy defects using open source well-known software i.e., large-scale atomic/molecular massively parallel simulator (LAMMPS). The interactions between the carbon atoms are modelled using Adaptive Intermolecular Reactive Empirical Bond Order (AIREBO) potential. An isogeometric analysis (IGA) based upon non-uniform rational B-spline (NURBS) is employed for approximation of single-layered graphene sheets deflection field and the governing equations are derived using nonlocal elasticity theory. The dependence of small-scale effects, chirality and different defect types on vibrational characteristic of graphene sheets is investigated in this comprehensive research work. In addition, numerical results are validated and compared with those achieved using other analysis, where an excellent agreement is found. The interesting results indicate that increasing the number of missing atoms can lead to decrease the natural frequencies of graphene sheets. It is seen that the degree of the detrimental effects differ with defect type. The Young's and shear modulus of the graphene with SV defects are much smaller than graphene with DV defects. It is also observed that Single Vacancy (SV) clusters cause more reduction in the natural frequencies of SLGS than Double Vacancy (DV) clusters. The effectiveness and the accuracy of the present IGA approach have been demonstrated and it is shown that the IGA is efficient, robust and accurate in terms of nanoplate problems.

• Bulletin of the Korean Chemical Society
• /
• v.26 no.4
• /
• pp.515-522
• /
• 2005

Enzyme-metal combo catalysis is described as a useful methodology for the synthesis of optically active compounds. The key point of the method is the use of enzyme and metal in combination as the catalysts for the complete transformation of racemic substrates to single enantiomeric products through dynamic kinetic resolution (DKR). In this approach, enzyme acts as an enantioselective resolving catalyst and metal does as a racemizing catalyst for the efficient DKR. Three kinds of enzyme-metal combinations - lipase-ruthenium, subtilisin-ruthenium, and lipase-palladium –have been developed as the catalysts for the DKRs of racemic alcohols, esters, and amines. The scope of the combination catalysts can be extended to the asymmetric transformations of ketones, enol acetates, and ketoximes via the DKRs. In most cases studied, enzyme-metal combo catalysis provided enantiomerically-enriched products in high yields.

• Structural Engineering and Mechanics
• /
• v.30 no.2
• /
• pp.177-190
• /
• 2008

Material failure behavior is generally dependent on loading rate. Especially in brittle and quasi-brittle materials, rate dependent material behavior can be significant. Empirical formulations are often used to predict the rate dependency, but such methods depend on extensive experimental works and are limited by practical constraints of physical testing. Numerical simulation can be an effective means for extracting knowledge about rate dependent behavior and for complementing the results obtained by testing. In this paper, the failure behavior of a brittle material under different loading rates is simulated by molecular dynamics analysis. A notched specimen is modeled by sub-million particles with a normalization scheme. Lennard-Jones potential is used to describe the interparticle force. Numerical simulations are performed with six different loading rates in a direct tensile test, where the loading velocity is normalized to the ratio of the pseudo-sonic speed. As a consequence, dynamic features are achieved from the numerical experiments. Remarkable failure characteristics, such as crack surface interaction/crack arrest, branching, and void nucleation, vary in case of the six loading cases. These characteristics are interpreted by the energy concept approach. This study provides insight into the change in dynamic failure mechanism under different loading rates.

• Tribology and Lubricants
• /
• v.36 no.5
• /
• pp.290-296
• /
• 2020

When two layers of carbon nanotube (CNT) arrays are loaded to mate, the free ends of individual CNTs come into contact at the interface of the two layers. This leads to a higher contact resistance due to a smaller contact region. However, when the free CNT ends of one array penetrate into the mating array, the contact region increases, effectively lowering the contact resistance. To explore the penetration of mating CNTs, we perform molecular dynamic simulations of a simple unit cell model, incorporating four CNTs in the lower array layer coupled with a single moving CNT on the upper layer. The interaction with neighboring CNTs is modelled by long-range carbon bond order potential (LCBOP I). The model structure is optimized by energy minimization through the conjugate gradient method. A NVT ensemble is used for maintain a room temperature during simulation. The time integration is performed through the velocity-Verlet algorithm. A significant vibrational motion of CNTs is captured when penetration is not available, resulting in a specific vibration mode with a high frequency. Due to this vibrational behavior, the random behaviors of CNT motion for predicting the penetration are confirmed under the specific gap distances between CNTs. Thus, the probability of penetration is examined according to the gap distance between CNTs in the lower array and the aspect ratio of CNTs. The penetration is significantly affected by the vibration mode due to the van der Waals forces between CNTs.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2003.10a
• /
• pp.339-346
• /
• 2003

Molecular dynamic simulations of nano indentation on single-crystal silicon (100) surface were performed using diamond indentor. Silicon substrate and diamond indentor were modeled diamond structure with Tersoff potential model. Phase transformation of silicon, incipient plastic deformation, change of incident temperature distribution are investigated through the change of potential energy distribution, displacement-load diagram, the change of kinetic energy distribution and displacements of silicon atoms. Phase transformation is highly localized and consists of a high-density region surrounding the tip. Axial load linearly increased according to the indenting depth. Number of atoms with high kinetic energy increased at the interface between substrate and indentor tip.

• Bulletin of the Korean Chemical Society
• /
• v.32 no.9
• /
• pp.3425-3428
• /
• 2011

Apolipoprotein A-I (Apo A-I) is a major component for high density lipoproteins (HDL). A number of mimetic peptides of Apo A-I were screened from the phase-displayed random peptide library by utilizing monoclonal antibodies (A12). Mimetic peptide for A12 epitope against Apo A-I was selected as CPFARLPVEHHDVVGL (pA1). From the BLAST search, the mimetic peptide pA1 had 40% homology with Apo A-I. As a result of the structural determination of this mimotope using homo/hetero nuclear 2D-NMR techniques and NMR-based distance geometry (DG)/molecular dynamic (MD) computations, DG structure had low penalty value of 0.3-0.7 ${\AA}^2$ and the total RMSD was 0.6-1.6 ${\AA}$. The mimotope pA1 exhibited characteristic conformation including a ${\beta}$-turn from Pro[7] to His[11].

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.19 no.12
• /
• pp.883-888
• /
• 2007

As a starting point of investigating what molecular dynamic simulations can reveal about the nature of atomic level of heating and cooling process, argon described by the LJ potential is considered. Stepwise heating and cooling of constant rates are simulated in the NPT (constant number, pressure and temperature) ensemble. Hysteresis is found due to the superheating and supercooling. Drastic change of volume and energy is involved with phase change, but the melting point can not be obtained by simply observing the changes of these quantities. Since liquid and solid phases can co-exist at the same temperature, Gibbs free energy should be calculated to find the temperature where the Gibbs free energy of liquid is equal to that of the solid since the equilibrium state is the state of minimum Gibbs free energy. The obtained melting temperature, $T^*=0.685$, is close to that of the experiment with only 2% error.

• Journal of the Korean Ceramic Society
• /
• v.34 no.7
• /
• pp.713-721
• /
• 1997

Molecular dynamic (MD) simulations with new interatomic potential function including the covalent bond were performed on the phase transition of $\alpha$-quartz-type GeO2 under high pressure. The optimized crystal structure and the pressure dependence of the lattice constant showed higher reproducibility than the previous models and were in very good agreement with the experimental data. A phase transition of $\alpha$-quartz and $\alpha$-quartz-type GeO2 by simulation was found approximately 24 GPa and 6-7 GPa, respectively. This phase transition involved an abrupt volume shrinkage and showed 4-6 coordination mixed structure with the increasing in the coordination number of cation.