• Multiscale and Multiphysics Mechanics
• /
• v.1 no.2
• /
• pp.171-188
• /
• 2016

Gradient enhanced theories of crystal plasticity enjoy great research interest. The focus of this work is on thermodynamically consistent modeling of grain size dependent hardening effects. In this contribution, we develop a model framework for damage coupled to gradient enhanced crystal thermoplasticity. The damage initiation is directly linked to the accumulated plastic slip. The theoretical setting is that of finite strains. Numerical results on single-crystalline metal showing the development of damage conclude the paper.

• Coupled systems mechanics
• /
• v.3 no.1
• /
• pp.27-52
• /
• 2014

For landfill monitoring and aftercare, long-term prognoses of emission and deformation behaviour are required. Landfills may be considered as heterogeneous porous soil-like structures, in which flow and transport processes of gases and liquids interact with local material degradation and mechanical deformation of the solid skeleton. Therefore, in the framework of continuous porous media mechanics a model is developed that permits the investigation of coupled mechanical, hydraulical and biochemical processes in municipal solid waste landfills.

• Proceedings of the Korean Institute of Electrical and Electronic Material Engineers Conference
• /
• 2009.11a
• /
• pp.111-111
• /
• 2009

In this work, we discuss simulation of surface acoustic wave device using Comsol Multiphysics. The structure SAW device based on piezoelectric thin film aluminum-nitride (AlN) on silicon was simulated. Some parameters of SAW device such as surface velocity, displacement of piezoelectric thin film were evaluated by software. Many modes and shapes of wave are also discussed in this paper. For evaluation physical parameters of AlN piezoelectric layer, the SAW resonator was modeled and simulation results were also compared with experiment results. we simulated arid evaluated the surface Rayleigh wave of AlN thin film on silicon substrate. Results simulation and experiment showed the surface velocity of AlN thin film was about 5200 m/s and shape of surface wave was also displayed. This paper has also proposed as method to study SAW characteristic of piezoelectric thin film and found out measurement values accurately of film such as stiffness matrix, piezoelectric matrix. These values are very important in calculation and design SAW device or MEMS device based on AlN piezoelectric layer.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03b
• /
• pp.479-482
• /
• 2008

In order to more closely examine the vacuum arc phenomena, it is necessary to predict the magnetohydrodynamic (MHD) characteristics by the multidisciplinary numerical modeling, which is coupled with the electromagnetic and hydrodynamic fields, simultaneously. In this study, the thermal-fluid characteristics of high current vacuum arcs were calculated by a commercial multiphysics package, ANSYS, in order to obtain Joule heat, Lorentz force and the interactions with flow variables. We assumed the diffused-mode arc within an AMF vacuum interrupter. It was found with four different currents that the temperature distributions on the anode surface are diffused uniformly without concentration in 7kA for both types (cup and coil-type). But the arc plasma transition and an increase of thermal flux density for increasing the applied current have caused the change of temperature distributions on the anode surface. We should need further studies on the two-way coupling method and radiation model for arc plasmas in order to accomplish the advanced analysis method for multiphysics.

• Nuclear Engineering and Technology
• /
• v.53 no.10
• /
• pp.3133-3150
• /
• 2021

This work deals with the validation of a high-fidelity multiphysics system coupling the Serpent 2 Monte Carlo neutron transport code with SUBCHANFLOW, a subchannel thermalhydraulics code, and TRANSURANUS, a fuel-performance analysis code. The results for a full-core pin-by-pin burnup calculation for the ninth operating cycle of the Temelín II VVER-1000 plant, which starts from a fresh core, are presented and assessed using experimental data. A good agreement is found comparing the critical boron concentration and a set of pin-level neutron flux profiles against measurements. In addition, the calculated axial and radial power distributions match closely the values reported by the core monitoring system. To demonstrate the modeling capabilities of the three-code coupling, pin-level neutronic, thermalhydraulic and thermomechanic results are shown as well. These studies are encompassed in the final phase of the EU Horizon 2020 McSAFE project, during which the Serpent-SUBCHANFLOW-TRANSURANUS system was developed.

• The Journal of Engineering Geology
• /
• v.32 no.4
• /
• pp.643-659
• /
• 2022

Geologic sequestration technologies such as CCS (carbon capture and storage), EGS (enhanced geothermal systems), and EOR (enhanced oil recovery) have been widely implemented in recent years, prompting evaluation of the mechanical stability of storage sites. As fluid injection can stimulate mechanical instability in storage layers by perturbing the stress state and pore pressure, poroelastic models considering various injection scenarios are required. In this study, we calculate the pore pressure, stress distribution, and vertical displacement along a surface using commercial finite element software (COMSOL); fault slips are subsequently simulated using PyLith, an open-source finite element software. The displacement fields, are obtained from PyLith is transferred back to COMSOL to determine changes in coseismic stresses and surface displacements. Our sequential use of COMSOL-PyLith-COMSOL for poroelastic modeling of fluid-injection and induced-earthquakes reveals large variations of pore pressure, vertical displacement, and Coulomb failure stress change during injection periods. On the other hand, the residual stress diffuses into the remote field after injection stops. This flow pattern suggests the necessity of numerical modeling and long-term monitoring, even after injection has stopped. We found that the time at which the Coulomb failure stress reaches the critical point greatly varies with the hydraulic and poroelastic properties (e.g., permeability and Biot-Willis coefficient) of the fault and injection layer. We suggest that an understanding of the detailed physical properties of the surrounding layer is important in selecting the injection site. Our numerical results showing the surface displacement and deviatoric stress distribution with different amounts of fault slip highlight the need to test more variable fault slip scenarios.

• Journal of Nuclear Fuel Cycle and Waste Technology(JNFCWT)
• /
• v.17 no.2
• /
• pp.227-234
• /
• 2019

A user-friendly modeling interface is developed for a process-based total system performance assessment framework (APro) specialized for a generic geological disposal system for high-level radioactive waste. The APro modeling interface is constructed using MATLAB, and the operator splitting scheme is used to combine COMSOL for simulation of multiphysics and PHREEQC for the calculation of geochemical reactions. As APro limits the modeling domain to the generic disposal system, the degree of freedom of the model is low. In contrast, the user-friendliness of the model is improved. Thermal, hydraulic, mechanical and chemical processes considered in the disposal system are modularized, and users can select one of multiple modules: "Default process" and multi "Alternative process". APro mainly consists of an input data part and calculation execution part. The input data are prepared in a single EXCEL file with a given format, and the calculation part is coded using MATLAB. The final results of the calculation are created as an independent COMSOL file for further analysis.

• Journal of Nuclear Fuel Cycle and Waste Technology(JNFCWT)
• /
• v.15 no.4
• /
• pp.321-332
• /
• 2017

In order to build a general model of a high-throughput uranium electrorefining process according to the electrode configuration, numerical analysis was conducted using the COMSOL Multiphysics V5.3 electrodeposition module with Ordinary Differential Equation (ODE) interfaces. The generated model was validated by comparing a current density-potential curve according to the distance between the anode and cathode and the electrode array, using a lab-scale (1kg U/day) multi-electrode electrorefiner made by the Korea Atomic Energy Research Institute (KAERI). The operating temperature was $500^{\circ}C$ and LiCl-KCl eutectic with 3.5wt% $UCl_3$ was used for molten salt. The efficiency of the uranium electrorefining apparatus was improved by lowering the cell potential as the distance between the electrodes decreased and the anode/cathode area ratio increased. This approach will be useful for constructing database for safety design of high throughput spent nuclear fuel electrorefiners.

• Progress in Superconductivity
• /
• v.13 no.3
• /
• pp.151-157
• /
• 2012

In this paper, we used E-J constitutive law and H-formulation to calculate magnetic field profile, current density, and magnetization of high temperature superconductor (HTS) placed in time varying applied magnetic field. Finite element method (FEM)-based software, Comsol Multiphysics 3.5a, was employed to simulate 2-dimensional model of a superconducting thin strip. The numerical results based on Kim's critical state model were compared with the case of strip in a perpendicular field in the Brandt's paper as well as experimental data observed by Scanning Hall Probe and SQUID.

• Multiscale and Multiphysics Mechanics
• /
• v.1 no.1
• /
• pp.65-76
• /
• 2016

In this study, an influence of local variation of nanoparticulate volume fraction on the homogenized elastic properties is investigated. It is well known that interface effect is dependent on the radius and volume fraction of reinforced nanofillers. However, there is no study on the multiscale modeling and analysis of polymer nanocomposites including polydispersed nanoparticles with consideration of interphase zone, which is dependent on the volume fraction of corresponding nanoparticles. As results of numerical examples, it is confirmed that an influence of local variation of nanoparticulate volume fraction should be considered for non-dilute system such as cluster of nanoparticles. Therefore representative volume element analysis is conducted by considering local variation of nanoparticle volume fraction in order to analyze the practical size of cell including hundreds of nanoparticles. It is expected that this study could be extended to the multiparticulate nanocomposite systems including polydispersed nanoparticles.