• Journal of the Geological Society of Korea
• /
• v.54 no.6
• /
• pp.677-682
• /
• 2018

Magnetic anomalies are sensitive to magnetic properties present in deep Earth and near surface structures. Such geophysical characteristics often can be quantified by numerical analyses. In this study, we developed a finite element method (FEM) approach to compute magnetic anomalies using COMOL $Multiphysics^{(R)}$. This FEM approach was verified by comparing its numerical results with the previously known analytic solution for a uniformly magnetized sphere. Then, we used the method to compute magnetic reversal patterns near mid-ocean ridge with various faulting scenarios. This COMSOL-based approach can be incorporated into advanced multi-physical numerical models to understand the Earth.

• Journal of the Geological Society of Korea
• /
• v.54 no.6
• /
• pp.683-694
• /
• 2018

Subduction has been the focal point of numerical studies for decades because it plays an important role in the Earth's mass and energy circulations and generates earthquakes and arc volcanoes which are closely related to the human lives. Among the studies on subduction, numerical modeling has been broadly applied to the quantitative studies on the subducting slab in the mantle which cannot be directly observed. In this study, we benchmark the numerical subduction modeling using a finite element package, COMSOL $Multiphysics^{(R)}$ and the results are consistent with the previously reported benchmark results.

• Geophysics and Geophysical Exploration
• /
• v.23 no.1
• /
• pp.1-12
• /
• 2020

In the development of geodynamic structures such as subduction and rift zones, a weakening mechanism is essential for localized weak zone formation in the lithosphere. Shear heating, a weakening mechanism, generates short-wavelength temperature elevation in the lithosphere; the increased temperature can reduce lithospheric strength and promote its breakup. A two-dimensional elastoplastic extensional basin model was used to conduct benchmarking based on previous numerical simulation studies to quantitatively analyze shear heating. The amount of shear heating was investigated by controlling the yield strength, extensional velocity, and strain- and temperature-dependent weakening. In the absence of the weakening mechanism, the higher yield strength and extensional velocity led to more vigorous shear heating. The reference model with a 100-MPa yield strength and 2-cm/year extension showed a temperature increase of ~ 50 K when the bulk extension was 20 km (i.e., 0.025 strain). However, in the yield-strength weakening mechanism, depending on the plastic strain and temperature, more efficient weakening induced stronger shear heating, which indicates positive feedback between the weakening mechanism and the shear heating. The rate of shear heating rapidly increased at the initial stage of deformation, and the rate decreased by 80% as the lithosphere weakened. This suggests that shear heating with the weakening mechanism can significantly influence the strength of relatively undamaged lithosphere.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.17 no.4
• /
• pp.419-427
• /
• 2011

The characteristics for hydrogen production and the thermochemical properties of high temperature steam electrolysis(HTSE) device have been numerically analyzed in a two-dimension, steady-state with using the COMSOL $Multiphysics^{(R)}$. The main parameters for the calculation are applied voltage, ASR(Area-specific Resistance), temperature and pressure of the inlet gas flow. The results showed that thermal-neutral voltage was 1.2454 V and rather than the cell temperature increases or decreases with increasing applied voltage by thermal-neutral voltage starting this voltage the temperature in high voltage tended to rise and temperature in the low voltage tended to fall. And with, increasing the values of ASR, temperature inside the cell and the hydrogen production rate were decreased.

• Journal of the Microelectronics and Packaging Society
• /
• v.17 no.1
• /
• pp.47-53
• /
• 2010

In this study, analyses on the stencil printing using solder paste were carried out. The key design parameters in the stencil printing process are printing conditions, stencil design, and solder paste properties. Among these parameters, the effects of physical properties of solder paste such as viscosity, surface tension, and contact angle on the stencil printing process were investigated. The analyses were performed for simple geometry and boundary conditions. In the analysis, solder paste was pushed into a stencil hole by pressure instead of printer pad. Considering the geometry and computational efficiency, axisymmetric analyses were adopted. A commercial software (COMSOL), which is well known in the area of micro-fluids analysis, was used. From the results, it was shown that viscosity of solder paste had an effect on the filling speed, while surface tension and contact angle had an effect on the filling shape.