• Steel and Composite Structures
• /
• v.48 no.4
• /
• pp.443-459
• /
• 2023

This paper investigates the stability of a bi-directional functionally graded (BD-FG) cylindrical beam made of imperfect concrete, taking into account size-dependency and the effect of geometry on its stability behavior. Both buckling and dynamic behavior are analyzed using the modified coupled stress theory and the classical beam theory. The BD-FG structure is created by using porosity-dependent FG concrete, with changing porosity voids and material distributions along the pipe radius, as well as uniform and nonuniform radius functions that vary along the beam length. Energy principles are used to generate partial differential equations (PDE) for stability analysis, which are then solved numerically. This study sheds light on the complex behavior of BD-FG structures, and the results can be useful for the design of stable cylindrical microstructures.

• Journal of the Korean Crystal Growth and Crystal Technology
• /
• v.32 no.1
• /
• pp.16-24
• /
• 2022

A computational study of combined thermal and solutal convection (double diffusive convection) in a sealed crystal growth reactor is presented, based on a two-dimensional numerical analysis of the nonlinear and strongly coupled partial differential equations and their associated boundary conditions. The average Nusselt numbers for the source regions are greater than those at the crystal regions for 9.73 × 10³ ≤ Gr_t ≤ 6.22 × 10⁵. The average Nusselt numbers for the source regions varies linearly and increases directly with the thermal Grashof number form 9.73 × 10³ ≤ Gr_t ≤ 6.22 × 10⁵ for aspect ratio, Ar (transport length-to-width) = 1 and 2. Additionally, the average Nusselt numbers for the crystal regions at Ar = 1 are much greater than those at Ar = 2. Also, the occurrence of one unicellular flow structure is caused by both the thermal and solutal convection, which is inherent during the physical vapor transport of Hg₂Br₂. When the aspect ratio of the enclosure increases, the fluid movement is hindered and results in the decrease of thermal buoyancy force.

• Advances in nano research
• /
• v.12 no.1
• /
• pp.81-99
• /
• 2022

The aim of current work is to evaluate thermo-electrical characteristics of graphene nanoplatelets Reinforced Composite (GNPRC) coupled with magneto-electro-elastic (MEE) face sheet. In this regard, a cylindrical smart nanocomposite made of GNPRC with an external MEE layer is considered. The bonding between the layers are assumed to be perfect. Because of the layer nature of the structure, the material characteristics of the whole structure is regarded as graded. Both mechanical and thermal boundary conditions are applied to this structure. The main objective of this work is to determine critical temperature and critical voltage as a function of thermal condition, support type, GNP weight fraction, and MEE thickness. The governing equation of the multilayer nanocomposites cylindrical shell is derived. The generalized differential quadrature method (GDQM) is employed to numerically solve the differential equations. This method is integrated with Deep Learning Network (DNN) with ADADELTA optimizer to determine the critical conditions of the current sandwich structure. This the first time that effects of several conditions including surrounding temperature, MEE layer thickness, and pattern of the layers of the GNPRC is investigated on two main parameters critical temperature and critical voltage of the nanostructure. Furthermore, Maxwell equation is derived for modeling of the MEE. The outcome reveals that MEE layer, temperature change, GNP weight function, and GNP distribution patterns GNP weight function have significant influence on the critical temperature and voltage of cylindrical shell made from GNP nanocomposites core with MEE face sheet on outer of the shell.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.31 no.6
• /
• pp.36-44
• /
• 2003

In this paper, an analytic solution method is proposed to overcome the numerical problems when the slewing dynamics of hybrid coordinate systems is investigated via time finite element analysis. It is shown that the dynamics of the hybrid coordinate systems is governed by the coupled dual differential equations for both slewing and structural modes. Structural modes are transformed into the time-based modal coordinates and analytic spatial propagation equations are derived for each space-dependent time mode. Slew angle history is obtained analytically by appropriate applications of the boundary conditions and structural propagation is re-calculated using the slew angle. Numerical examples are demonstrated to validate the proposed analytic method in comparison to the existing state transition matrix method.

• Proceedings of the SAREK Conference
• /
• 2008.06a
• /
• pp.109-114
• /
• 2008

A dynamic model has been developed to simulate dynamic operation of a real double-effect absorption chiller. Dynamic behavior of working fluids in main components was modeled in first-order nonlinear differential equations based on heat and mass balances. Mass transport mechanisms among the main components were modeled by valve throttling, 'U' tube overflow and solution sub-cooling. The nonlinear dynamic equations coupled with the subroutines to calculate thermodynamic properties of working fluids were solved by a numerical method. The dynamic performance of the model was compared with the test data of a commercial medium chiller. The model showed a good agreement with the test data except for the first 5,000 seconds during which different flow rates of the weak solution caused some discrepancy. It was found that the chiller dynamics is governed by the inlet temperatures of the cooling water and the chilled water when the heat input to the chiller is relatively constant.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.20 no.12
• /
• pp.781-788
• /
• 2008

A dynamic model has been developed to simulate dynamic operation of a real double-effect absorption chiller. Dynamic behavior of working fluids in main components was modeled in first-order nonlinear differential equations based on heat and mass balances. Mass transport mechanisms among the main components were modeled by valve throttling, 'U' tube overflow and solution sub-cooling. The nonlinear dynamic equations coupled with the subroutines to calculate thermodynamic properties of working fluids were solved by a numerical method. The dynamic performance of the model was compared with the test data of a commercial medium chiller. The model showed a good agreement with the test data except for the first 5,000 seconds during which different flow rates of the weak solution caused some discrepancy. It was found that the chiller dynamics is governed by the inlet temperatures of the cooling water and the chilled water when the heat input to the chiller is relatively constant.

• Nuclear Engineering and Technology
• /
• v.14 no.2
• /
• pp.78-85
• /
• 1982

The steady-state rotation of plasma centrifuge is theoretically analyzed to understand the physics of rotating plasmas and its feasibility for isotope separation. The centrifuge system under consideration consists of an annular gap between coaxial cylindrical anode and cathode in the presence of an externally-applied axial magnetic field. A problem for coupled partial differential equations describing centrifuge fields is formulated on the basis of the magnetohydrodynamic equations. Two-dimensional solutions are found analytically in the form of Fourier-Bessel series. The current density and velocity distributions are discussed in terms of the Hartmann number and the geometrical parameter of the system. At typical conditions, rotational speeds of the plasma up to the order of 10$^4$m/sec are achievable, and increase either with increasing Hartmann number, or with increasing ratio of the axial length to the inner radius of the cylinder. In view of much higher speeds of rotation which can be achieved in plasma centrifuge, it is expected that its efficiency is superior to mechanically driven gas centrifuges.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.2
• /
• pp.113-125
• /
• 2007

Improved numerical method to obtain the exact stiffness matrices is newly proposed to perform the spatially coupled elastic and stability analyses of non-symmetric and open/closed thin-walled beam on elastic foundation. This method overcomes drawbacks of the previous method to evaluate the exact stiffness matrix for the spatially coupled stability analysis of thin-walled beam-column This numerical technique is accomplished via a generalized eigenproblem associated with 14 displacement parameters by transforming equilibrium equations to a set of first order simultaneous ordinary differential equations. Next polynomial expressions as trial solutions are assumed for displacement parameters corresponding to zero eigenvalues and the eigenmodes containing undetermined parameters equal to the number of zero eigenvalues are determined by invoking the identity condition. And then the exact displacement functions are constructed by combining eigensolutions and polynomial solutions corresponding to non-zero and zero eigenvalues, respectively. Consequently an exact stiffness matrix is evaluated by applying the member force-deformation relationships to these displacement functions. In order to illustrate the accuracy and the practical usefulness of this study, the numerical solutions are compared with results obtained from the thin-walled beam and shell elements.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.26 no.2
• /
• pp.135-147
• /
• 2015

With the miniaturization and digitalization of electronics industry, demand for Multi-Layer Ceramic Capacitor(MLCC) has increased steadily because of its various applications such as DC Blocking, Decoupling and Filtering etc. The modeling techniques of MLCC has been studied for a long time but most of these modeling method can only be applied after measurement and this has some losses of material, time in both production stage and measurement stage. This paper proposes the modeling method which can predict the frequency characteristics of MLCC from structure data and material data in design stage. The impedance of N-Layer Capacitor can be expressed in differential mathematical form based on coupled transmission line equations. By using this formula, we can predict the impedance of MLCC. As a result, proposed modeling is correspond with simulation, and it takes much less time to obtain the result than the simulation.

• Journal of Computational Design and Engineering
• /
• v.3 no.4
• /
• pp.349-362
• /
• 2016

The numerical solutions of unsteady hydromagnetic natural convection Couette flow of a viscous, incompressible and electrically conducting fluid between the two vertical parallel plates in the presence of thermal radiation, thermal diffusion and diffusion thermo are obtained here. The fundamental dimensionless governing coupled linear partial differential equations for impulsive movement and uniformly accelerated movement of the plate were solved by an efficient Finite Element Method. Computations were performed for a wide range of the governing flow parameters, viz., Thermal diffusion (Soret) and Diffusion thermo (Dufour) parameters, Magnetic field parameter, Prandtl number, Thermal radiation and Schmidt number. The effects of these flow parameters on the velocity (u), temperature (${\theta}$) and Concentration (${\phi}$) are shown graphically. Also the effects of these pertinent parameters on the skin-friction, the rate of heat and mass transfer are obtained and discussed numerically through tabular forms. These are in good agreement with earlier reported studies. Analysis indicates that the fluid velocity is an increasing function of Grashof numbers for heat and mass transfer, Soret and Dufour numbers whereas the Magnetic parameter, Thermal radiation parameter, Prandtl number and Schmidt number lead to reduction of the velocity profiles. Also, it is noticed that the rate of heat transfer coefficient and temperature profiles increase with decrease in the thermal radiation parameter and Prandtl number, whereas the reverse effect is observed with increase of Dufour number. Further, the concentration profiles increase with increase in the Soret number whereas reverse effect is seen by increasing the values of the Schmidt number.