• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.25 no.4
• /
• pp.315-321
• /
• 2012

In this study, a sequential multiscale homogenization method to characterize the effective thermal conductivity of nano particulate polymer nanocomposites is proposed through a molecular dynamics(MD) simulations and a finite element-based homogenization method. The thermal conductivity of the nanocomposites embedding different-sized nanoparticles at a fixed volume fraction of 5.8% are obtained from MD simulations. Due to the Kapitza thermal resistance, the thermal conductivity of the nanocomposites decreases as the size of the embedded nanoparticle decreases. In order to describe the nanoparticle size effect using the homogenization method with accuracy, the Kapitza interface in which the temperature discontinuity condition appears and the effective interphase zone formed by highly densified matrix polymer are modeled as independent phases that constitutes the nanocomposites microstructure, thus, the overall nanocomposites domain is modeled as a four-phase structure consists of the nanoparticle, Kapitza interface, effective interphase, and polymer matrix. The thermal conductivity of the effective interphase is inversely predicted from the thermal conductivity of the nanocomposites through the multiscale homogenization method, then, exponentially fitted to a function of the particle radius. Using the multiscale homogenization method, the thermal conductivities of the nanocomposites at various particle radii and volume fractions are obtained, and parametric studies are conducted to examine the effect of the effective interphase on the overall thermal conductivity of the nanocomposites.

• Korean Chemical Engineering Research
• /
• v.57 no.4
• /
• pp.475-483
• /
• 2019

Special chemical warfare agents are lethal gases that attack the human respiratory system. One of such gases are blood agents that react with the irons present in the electron transfer system of the human body. This reaction stops internal respiration and eventually causes death. The molecular sizes of these agents are smaller than the pores of an activated carbon, making chemical adsorption the only alternative method for removing them. In this study, we carried out a Computational Fluid Dynamics simulation by passing a blood agent: cyanogen chloride gas through an SG-1 gas mask canister developed by SG Safety Corporation. The adsorption bed consisted of a Silver-Zinc-Molybdenum-Triethylenediamine activated carbon impregnated with copper, silver, zinc and molybdenum ions. The kinetic analysis of the chemical adsorption was performed in accordance with the test procedure for the gas mask canister and was validated by the kinetic data obtained from experimental results. We predicted the dynamic behaviors of the main variables such as the pressure drop inside the canister and the amount of gas adsorbed by chemisorption. By using a granular packed bed instead of the Ergun equation that is used to model porous materials in Computational Fluid Dynamics, applicable results of the activated carbon were obtained. Dynamic simulations and flow analyses of the chemical adsorption with varying gas flow rates were also executed.

• Advances in environmental research
• /
• v.7 no.2
• /
• pp.139-159
• /
• 2018

Green Leaf Volatiles (GLVs) is a class of biogenically emitted oxygenated hydrocarbons that have been identified as a potential source of Secondary Organic Aerosols (SOA) via aqueous oxidation. The physico-chemical properties of GLVs are vital to understanding their fate and transport in the atmosphere via fog processing, but few experimental data are available. We studied the aqueous solubility, 1-octanol/water partition coefficient, and Henry's law constant ($K_H$) of five GLVs at $25^{\circ}C$: methyl jasmonate, methyl salicylate, 2-methyl-3-buten-2-ol, cis-3-hexen-1-ol, and cis-3-hexenyl acetate. Henry's law constant was also measured at temperatures and ionic strengths typical of fog. Experimental values are compared to scarcely-available literature values, as well as estimations using group and bond contribution methods, property-specific correlations and molecular dynamics simulations. From these values, the partition coefficients to the air-water interface were also calculated. The large Henry's law constant of methyl jasmonate ($8091{\pm}1121M{\cdot}atm^{-1}$) made it the most significant GLV for aqueous phase photochemistry. The HENRYWIN program's bond contribution method from the Estimation Programs Interface Suite (EPI Suite) produced the best estimate of the Henry's constant for GLVs. Estimations of 1-octanol/water partition coefficient and solubility are best when correlating an experimental value of one to find the other. Finally, the scavenging efficiency was calculated for each GLV indicating aqueous phase processing will be most important for methyl jasmonate.

• Applied Chemistry for Engineering
• /
• v.10 no.2
• /
• pp.268-274
• /
• 1999

The chitosan solution was prepared by dissolving chitosan into 2 wt % aqueous acetic acid solution and then chitosan beads were made by sol-gel method. The average molecular weight and the degree of deacetylation of the chitosan used here were determined to be $8.2{\times}10^5$ and 85%, respectively. chitosan beads were highly porous which was confirmed by SEM photography and BET. Adsorption equilibrium of $Cu^{2+}$ on porous chitosan beads could be represented by Sips equation. The diffusion of cupric ions in the chitosan beads could be explained by pore and surface diffusion mechanisms. Adsorption dynamics of $Cu^{2+}$ in fixed-bed could be simulated by linear driving force approximation (LDFA). It was proven that porous chitosan beads manufactured in this work are good adsorbents for the removal of $Cu^{2+}$.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.29 no.6
• /
• pp.577-582
• /
• 2016

In this study, a multiscale method for solving a thermoelasticity problem for interphase in the polymeric nanocomposites is developed. Molecular dynamics simulation and finite element analysis were numerically combined to describe the geometrical boundaries and the local mechanical response of the interfacial region where the polymer networks were highly interacted with the nanoparticle surface. Also, the micrmechanical thermoelasticity equations were applied to the obtained equivalent continuum unit to compute the growth of interphase thickness according to the size of nanoparticles, as well as the thermal phase transition behavior at a wide range of temperatures. Accordingly, the equivalent continuum model obtained from the multiscale analysis provides a meaningful description of the thermoelastic behavior of interphase as well as its nanoparticle size effect on thermoelasticity at both below and above the glass transition temperature.

• Proceedings of the Korean Vacuum Society Conference
• /
• 2011.02a
• /
• pp.407-408
• /
• 2011

Natural boron consists of two stable isotopes 10B and 11B with natural abundance of 18.8 atom percent of 10B and 81.2 atom percent of 11B. The thermal neutron absorption cross-section for 10B and 11B are 3837 barn and 0.005 barn respectively. 10B enriched specific compounds are used for control rods and as a reactor coolant additives. In this work 2 methods for boron enrichment were analysed: 1) Gas irradiation in static conditions. Dissociation occurs due to multiphoton absorption by specific isotopes in appropriately tuned laser field. IR shifted laser pulses are usually used in combination with increasing the laser intensity also improves selectivity up to some degree. In order to prevent recombination of dissociated molecules BCl3 is mixed with H2S 2) SILARC method. Advantages of this method: a) Gas cooling is helpful to split and shrink boron isotopes absorption bands. In order to achieve better selectivity BCl3 gas has to be substantially rarefied (~0.01%-5%) in mixture with carrier gas. b) Laser intensity is lower than in the first method. Some preliminary calculations of dissociation and recombination with carrier gas molecules energetics for both methods will be demonstrated Boron separation in SILARC method can be represented as multistage process: 1) Mixture of BCl3 with carrier gas is putted in reservoir 2) Gas overcooling due to expansion through Laval nozzle 3) IR multiphoton absorption by gas irradiated by specifically tuned laser field with subsequent gradual gas condensation in outlet chamber It is planned to develop software which includes these stages. This software will rely on the following available software based on quantum molecular dynamics in external quantized field: 1) WavePacket: Each particle is treated semiclassicaly based on Wigner transform method 2) Turbomole: It is based on local density methods like density of functional methods (DFT) and its improvement- coupled clusters approach (CC) to take into account quantum correlation. These models will be used to extract information concerning kinetic coefficients, and their dependence on applied external field. Information on radiative corrections to equation of state induced by laser field which take into account possible phase transition (or crossover?) can be also revealed. This mixed phase equation of state with quantum corrections will be further used in hydrodynamical simulations. Moreover results of these hydrodynamical simulations can be compared with results of CFD calculations. The first reasonable question to ask before starting the CFD simulations is whether turbulent effects are significant or not, and how to model turbulence? The questions of laser beam parameters and outlet chamber geometry which are most optimal to make all gas volume irradiated is also discussed. Relationship between enrichment factor and stagnation pressure and temperature based on experimental data is also reported.

• Journal of Ginseng Research
• /
• v.47 no.3
• /
• pp.408-419
• /
• 2023

Background: Ginsenoside compound K (CK), the main active metabolite in Panax ginseng, has shown good safety and bioavailability in clinical trials and exerts neuroprotective effects in cerebral ischemic stroke. However, its potential role in the prevention of cerebral ischemia/reperfusion (I/R) injury remains unclear. Our study aimed to investigate the molecular mechanism of ginsenoside CK against cerebral I/R injury. Methods: We used a combination of in vitro and in vivo models, including oxygen and glucose deprivation/reperfusion induced PC12 cell model and middle cerebral artery occlusion/reperfusion induced rat model, to mimic I/R injury. Intracellular oxygen consumption and extracellular acidification rate were analyzed by Seahorse multifunctional energy metabolism system; ATP production was detected by luciferase method. The number and size of mitochondria were analyzed by transmission electron microscopy and MitoTracker probe combined with confocal laser microscopy. The potential mechanisms of ginsenoside CK on mitochondrial dynamics and bioenergy were evaluated by RNA interference, pharmacological antagonism combined with co-immunoprecipitation analysis and phenotypic analysis. Results: Ginsenoside CK pretreatment could attenuate mitochondrial translocation of DRP1, mitophagy, mitochondrial apoptosis, and neuronal bioenergy imbalance against cerebral I/R injury in both in vitro and in vivo models. Our data also confirmed that ginsenoside CK administration could reduce the binding affinity of Mul1 and Mfn2 to inhibit the ubiquitination and degradation of Mfn2, thereby elevating the protein level of Mfn2 in cerebral I/R injury. Conclusion: These data provide evidence that ginsenoside CK may be a promising therapeutic agent against cerebral I/R injury via Mul1/Mfn2 mediated mitochondrial dynamics and bioenergy.

• Journal of the Korean Institute of Intelligent Systems
• /
• v.25 no.6
• /
• pp.529-535
• /
• 2015

Biological networks have been handled with the static concept. However, life phenomena in cells occur depending on the cellular state and the external environment, and only a few proteins and their interactions are selectively activated. Therefore, we should adopt the dynamic network concept that the structure of a biological network varies along the flow of time. This concept is effective to analyze the progressive transition of the disease. In this paper, we applied the proposed method to Alzheimer's disease to analyze the structural and functional characteristics of the disease network. Using gene expression data and protein-protein interaction data, we constructed the sub-networks in accordance with the progress of disease (normal, early, middle and late). Based on this, we analyzed structural properties of the network. Furthermore, we found module structures in the network to analyze the functional properties of the sub-networks using the gene ontology analysis (GO). As a result, it was shown that the functional characteristics of the dynamics network is well compatible with the stage of the disease which shows that it can be used to describe important biological events of the disease. Via the proposed approach, it is possible to observe the molecular network change involved in the disease progression which is not generally investigated, and to understand the pathogenesis and progression mechanism of the disease at a molecular level.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.3
• /
• pp.137-145
• /
• 2014

We present a design sensitivity analysis(DSA) method for multiscale problems based on bridging scale decomposition. In this paper, we utilize a bridging scale method for the coupled system analysis. Since the analysis of full MD systems requires huge amount of computational costs, a coupled system of MD-level and continuum-level simulation is usually preferred. The information exchange between the MD and continuum levels is taken place at the MD-continuum boundary. In the bridging scale method, a generalized Langevin equation(GLE) is introduced for the reduced MD system and the GLE force using a time history kernel is applied at the boundary atoms in the MD system. Therefore, we can separately analyze the MD and continuum level simulations, which can accelerate the computing process. Once the simulation of coupled problems is successful, the need for the DSA is naturally arising for the optimization of macro-scale design, where the macro scale performance of the system is maximized considering the micro scale effects. The finite difference sensitivity is impractical for the gradient based optimization of large scale problems due to the restriction of computing costs but the analytical sensitivity for the coupled system is always accurate. In this study, we derive the analytical design sensitivity to verify the accuracy and applicability to the design optimization of the coupled system.

• Bulletin of the Korean Chemical Society
• /
• v.7 no.5
• /
• pp.367-371
• /
• 1986

Dynamic and equilibrium properties of n-alkane chains immersed in solvent molecules have been investigated by a molecular dynamics method. The n-alkane chain is assumed to be a chain of elements (CH$_2$) interconnected by bonds having a fixed bond length and bond angle, but each bond of the chain is allowed to execute hindered internal rotation. We studied the effect of the number of the chain elements (N$_c$ = 10, 15 and 20) on the equilibrium properties of the system, e.g., the pair correlation functions between a chain element and solvent molecules, g$_{cs}$(r), and between the chain elements, g$_{cc}$(r), and the configurational properties such as the mean-square end-to-end distance < R$^2$ >, the mean-square radius of gyration < S$^2$ >, and the eigenvalues of the moment-of-inertia tensor < S$_i^2$ > / < S$^2$ > (i = 1, 2 and 3). We also studied the dynamic properties of the system, e.g., the autocorrelation function C(A;t) where A = R$^2$(t), = S$^2$(t), or = ${\vec{V}}(t)({\vec{V}}$ = velocity of the center of mass), and the diffusion coefficient D. The g$_{cs}$(r)'s are almost equal irrespective of the change of Nc while g$_{cc}$(r) becomes larger as N$_c$ increases; The MD computed configurational properties < R$^2$2 > and < S$^2$ > were found to be a little different from the values calculated from the statistical equations of < R$^2$ > and < S$^2$ >, it may be due to the fact that our model for the MD simulations includes a long-range volume effect. From the < S$_i^2$ > / < S$^2$ >, it is found that the chain molecule has a nearly spherical shape irrespective of the variation of N$_c$. For the dynamic properties we found that the C(R$^2$;t) and C(S$^2$;t) of lower N$_c$ decay faster than those of higher N$_c$, while the C($\vec V$;t) of the center of mass in the chain is weakly dependent on the N$_c$. The center of mass diffusion coefficient D$_c$ decreases as N$_c$ increases while the end point diffusion coefficient D$_e$ is nearly equal irrespective of the change of N$_c$.