• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.41 no.6
• /
• pp.507-515
• /
• 2017

Flexible multibody simulations are widely used in the industry to design mechanical systems. In flexible multibody dynamics, deformation coordinates are described either relatively in the body reference frame that is floating in the space or in the inertial reference frame. Moreover, these deformation coordinates are generated based on the discretization of the body according to the finite element approach. Therefore, the formulation of the flexible multibody system always deals with a huge number of degrees of freedom and the numerical solution methods require a substantial amount of computational time. Parallel computational methods are a solution for efficient computation. However, most of the parallel computational methods are focused on the efficient solution of large-sized linear equations. For multibody analysis, we need to develop an efficient formulation that could be suitable for parallel computation. In this paper, we developed a subsystem synthesis method for a flexible multibody system and proposed efficient parallel computational schemes based on the OpenMP API in order to achieve efficient computation. Simulations of a rotating blade system, which consists of three identical blades, were carried out with two different parallel computational schemes. Actual CPU times were measured to investigate the efficiency of the proposed parallel schemes.

• Journal of Chest Surgery
• /
• v.44 no.2
• /
• pp.89-98
• /
• 2011

Background: Although considerable efforts have been made to improve the graft patency in coronary artery bypass surgery, the role of biomechanical factors remains underrecognized. The aim of this study is to investigate the influences of geometric configurations of the bypass graft on hemodynamic characteristics in relation to anastomosis. Materials and Methods: The Numerical analysis focuses on understanding the flow patterns for different values of inlet and distal diameters and graft angles. The Blood flow field is treated as a two-dimensional incompressible laminar flow. A finite volume method is adopted for discretization of the governing equations. The Carreau model is employed as a constitutive equation for blood. In an attempt to obtain the optimal aorto-coronary bypass conditions, the blood flow characteristics are analyzed using in vitro models of the end-to-side anastomotic angles of $45^{\circ}$, $60^{\circ}$ and $90^{\circ}$. To find the optimal graft configurations, the mass flow rates at the outlets of the four models are compared quantitatively. Results: This study finds that Model 3, whose bypass diameter is the same as the inlet diameter of the stenosed coronary artery, delivers the largest amount of blood and the least pressure drop along the arteries. Conclusion: Biomechanical factors are speculated to contribute to the graft patency in coronary artery bypass grafting.

• Journal of Navigation and Port Research
• /
• v.27 no.5
• /
• pp.539-546
• /
• 2003

The numerical model of the flow analysis by finite element technique is described. The Galerkin method is employed for spatial discretization Two step explicit finite element scheme is used to discretize the time function, which has advantage in problems treating large numbers of elements and unsteady state. Two dimensional hydrodynamic model considering moving boundary condition is developed. Also it applied flow model which develop on flow portion of ideal fluid in the model flume and verified, and the results of this study confirm the efficiency of moving boundary treatment in Jeju harbor. The computed results have shown the good adaptability of moving boundary condition From these studies, it can be concluded that the present method is a useful and effective tool in tidal flow analysis.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.34 no.9
• /
• pp.45-52
• /
• 2006

In this study, a nonlinear control by feedback linearization method, one of nonlinear control schemes based on the nonlinear model, is proposed to suppress the flutter of a supersonic composite panel using piezoelectric materials. Most of the previous panel flutter controllers are the LQR(Linear Quadratic Regulator) which is based on the linear model. A nonlinear feedback linearizing controller proposed in this study considers the nonlinear characteristics of the system model. We use the actuator implemented by piezoceramic PZT. Using the principle of virtual displacements and a finite element discretization with the conforming four-node rectangular element, we first derive the discretized dynamic equations of motion, which are transformed into a nonlinear coupled-modal equations of motion of state space form. The effectiveness of the proposed method is also compared with the LQR based on the linear model through numerical simulations in the time domain using the Newmark method.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.26 no.4
• /
• pp.775-786
• /
• 2002

A Mixed Volume and Boundary Integral Equation Method is applied for the effective analysis of elastic wave scattering problems and plane elastostatic problems in unbounded solids containing general anisotropic inclusions and voids or isotropic inclusions. It should be noted that this newly developed numerical method does not require the Green's function for anisotropic inclusions to solve this class of problems since only Green's function for the unbounded isotropic matrix is involved in their formulation for the analysis. This new method can also be applied to general two-dimensional elastodynamic and elastostatic problems with arbitrary shapes and number of anisotropic inclusions and voids or isotropic inclusions. In the formulation of this method, the continuity condition at each interface is automatically satisfied, and in contrast to finite element methods, where the full domain needs to be discretized, this method requires discretization of the inclusions only. Finally, this method takes full advantage of the pre- and post-processing capabilities developed in FEM and BIEM. Through the analysis of plane elastostatic problems in unbounded isotropic matrix with orthotropic inclusions and voids or isotropic inclusions, and the analysis of plane wave scattering problems in unbounded isotropic matrix with isotropic inclusions and voids, it will be established that this new method is very accurate and effective for solving plane wave scattering problems and plane elastic problems in unbounded solids containing general anisotropic inclusions and voids/cracks or isotropic inclusions.

• Advances in nano research
• /
• v.7 no.3
• /
• pp.209-221
• /
• 2019

Graphene, with impressive electronic properties, have high potential in the microelectronic field. However, graphene itself is a zero bandgap material which is not suitable for digital logic gates and its application. Thus, much focus is on graphene nanoribbons (GNRs) that are narrow strips of graphene. During GNRs fabrication process, the occurrence of defects that ultimately change electronic properties of graphene is difficult to avoid. The modelling of GNRs with defects is crucial to study the non-idealities effects. In this work, nearest-neighbor tight-binding (TB) model for GNRs is presented with three main simplifying assumptions. They are utilization of basis function, Hamiltonian operator discretization and plane wave approximation. Two major edges of GNRs, armchair-edged GNRs (AGNRs) and zigzag-edged GNRs (ZGNRs) are explored. With single vacancy (SV) defects, the components within the Hamiltonian operator are transformed due to the disappearance of tight-binding energies around the missing carbon atoms in GNRs. The size of the lattices namely width and length are varied and studied. Non-equilibrium Green's function (NEGF) formalism is employed to obtain the electronics structure namely band structure and density of states (DOS) and all simulation is implemented in MATLAB. The band structure and DOS plot are then compared between pristine and defected GNRs under varying length and width of GNRs. It is revealed that there are clear distinctions between band structure, numerical DOS and Green's function DOS of pristine and defective GNRs.

• Advances in nano research
• /
• v.10 no.4
• /
• pp.359-371
• /
• 2021

A tumor immune interaction is a main topic of interest in the last couple of decades because majority of human population suffered by tumor, formed by the abnormal growth of cells and is continuously interacted with the immune system. Because of its wide range of applications, many researchers have modeled this tumor immune interaction in the form of ordinary, delay and fractional order differential equations as the majority of biological models have a long range temporal memory. So in the present work, tumor immune interaction in fractional form provides an excellent tool for the description of memory and hereditary properties of inter and intra cells. So the interaction between effector-cells, tumor cells and interleukin-2 (IL-2) are modeled by using the definition of Caputo fractional order derivative that provides the system with long-time memory and gives extra degree of freedom. Moreover, in order to achieve more efficient computational results of fractional-order system, a discretization process is performed to obtain its discrete counterpart. Furthermore, existence and local stability of fixed points are investigated for discrete model. Moreover, it is proved that two types of bifurcations such as Neimark-Sacker and flip bifurcations are studied. Finally, numerical examples are presented to support our analytical results.

• Advances in nano research
• /
• v.12 no.2
• /
• pp.213-221
• /
• 2022

Silicon carbide (SiC) is a binary carbon-silicon compound. In its two-dimensional form, monolayer SiC is composed of a monolayer carbon and silicon atoms constructed as a honeycomb lattice. SiC has recently been receiving increasing attention from researchers owing to its intriguing electronic properties. In this present work, SiC nanoribbons (SiCNRs) are modelled and simulated to obtain accurate electronic properties, which can further guide fabrication processes, through bandgap engineering. The primary objective of this work is to obtain the electronic properties of monolayer SiCNRs by applying numerical computation methods using nearest-neighbour tight-binding models. Hamiltonian operator discretization and approximation of plane wave are assumed for the models and simulation by applying the basis function. The computed electronic properties include the band structures and density of states of monolayer SiCNRs of varying width. Furthermore, the properties are compared with those of graphene nanoribbons. The bandgap of ASiCNR as a function of width are also benchmarked with published DFT-GW and DFT-GGA data. Our nearest neighbour tight-binding (NNTB) model predicted data closer to the calculations based on the standard DFT-GGA and underestimated the bandgap values projected from DFT-GW, which takes in account the exchange-correlation energy of many-body effects.

• Advances in aircraft and spacecraft science
• /
• v.9 no.5
• /
• pp.449-462
• /
• 2022

Space-telescopes placed in the Sun-Earth second Lagrange point (L2) observe the sky following a scan strategy that is usually based on a spin-precession motion. Knowing which regions of the sky will be more observed by the instrument is important for the science operations and the instrument calibration. Computing sky observation parameters numerically (discretizing time and the sky) can consume large amounts of time and computational resources, especially when high resolution isrequired.This problem becomesmore critical if quantities are evaluated at detector level instead of considering the instrument entire Field of View (FoV). In previous studies, the authors have derived analytic solutions for quantities that characterize the observation of each point in the sky in terms of observation time according to the scan strategy parameters and the instrument FoV. Analytic solutions allow to obtain results faster than using numerical methods as well as capture detailed characteristics which can be overseen due to discretization limitations. The original approach is based on the analytic expression of the instrument trace over the sky. Such equations are implicit and thusrequiresthe use of numeric solversto compute the quantities.In this work, a new and simpler approach for computing one ofsuch quantities(mean observation time) is presented.The quantity is first computed for pure spin motion and then the effect of the spin axis precession is incorporated under the assumption that the precession motion is slow compared to the spin motion.In this sense, this new approach further simplifies the analytic approach, sparing the use of numeric solvers, which reduces the complexity of the implementation and the computing time.

• Nuclear Engineering and Technology
• /
• v.56 no.9
• /
• pp.3950-3960
• /
• 2024

As a solution to the problem of spent nuclear fuels (SNFs), the disposal of SNF has gained attention from nations using nuclear energy because of hazards posed to the ecosystem. Among many proposed solutions, the most promising method is to dispose of SNF in a deep geological repository (DGR) which utilizes the multi-barrier concept developed by Finland and Sweden. Here, a new fully-coupled Thermal-Hydraulic (TH) code HADES (High-level rAdionuclide Disposal Evaluation Simulator) is developed using the MOOSE framework. This new code suggests basic numerical tools, such as a non-linear solver and finite element discretization, to assess the safety performance of disposal systems. The new TH code considered various TH behavior using Richards' flow approach, assuming gas pressure is constant. The HADES showed promising results when it was compared to various TH codes validated from DECOVAELX-THMC projects. When the single-canister model was utilized to estimate the TH behavior of the Korean Reference disposal System, although it showed significant saturation reduction due to the evaporation of water, the temperature was maintained under the thermal criteria limit, which is 100 ℃. In addition, the new code estimated temperature and degree of saturation of the multi-canisters model, considering two or three canisters, it showed a slightly lower temperature, 5 ℃, than the single-canister model. From these results, the following are concluded: (1) the new TH code contribute to an additional integrity by estimating TH behavior of KRS; (2) however, due to limitations in single-canister simulation, it is recommended to use multi-canisters simulation to estimate TH behavior accurately. Therefore, this model is anticipated not only to help licensing applications and estimation of various multi-physics phenomena and multi-canister at the disposal site.