• Steel and Composite Structures
• /
• v.46 no.6
• /
• pp.839-856
• /
• 2023

This work presents a simple four-unknown refined integral plate theory for mechanical and thermal buckling behaviors of functionally graded (FG) plates resting on Visco-Pasternak foundations. The proposed refined high order shear deformation theory has a new displacement field which includes indeterminate integral variables and contains only four unknowns in which any shear correction factor not used, with even less than the conventional theory of first shear strain (FSDT). Governing equations are deduced from the principle of minimum total potential energy and a Navier type analytical solution is adopted for simply supported FG plates. The Visco-Pasternak foundations is considered by adding the impact of damping to the usual foundation model which characterized by the linear Winkler's modulus and Pasternak's foundation modulus. The accuracy of the present model is demonstrated by comparing the computed results with those available in the literature. Some numerical results are presented to show the impact of material index, elastic foundation type, and damping coefficient of the foundation, on the mechanical and thermal buckling behaviors of FG plates.

• Steel and Composite Structures
• /
• v.45 no.1
• /
• pp.67-82
• /
• 2022

In this paper, buckling and free vibration of imperfect, functionally graded beams, including porosities, are investigated, using a higher order shear strain theory. Due to defects during the manufacturing process, micro porosities may appear in the material, hence the appearance of this imperfection in the structure. The material properties of the beams are assumed to vary regularly, with power and sigmoid law, in the direction of thickness. A novel porosity distribution affecting the functionally graded volume fraction is presented. For the compact formulation used for cementite-based materials and already used in P-FGM, we have adapted it for the distribution of S-FGM. The equations of motion in the FG beam are derived using Hamilton's principle. The boundary conditions for beam FG are assumed to be simply supported. Navier's solution is used to obtain the closed form solutions of the FG beam. The numerical results of this work are compared with those of other published research to verify accuracy and reliability. The comparisons of different shear shape functions, the influence of porosity, thickness and inhomogeneity parameters on buckling and free vibration of the FG beam are all discussed. It is established that the present work is more precise than certain theories developed previously.

• Structural Engineering and Mechanics
• /
• v.91 no.1
• /
• pp.1-23
• /
• 2024

This paper formulates a new integral shear deformation shell theory to investigate the free vibration response of carbon nanotube (CNT) reinforced structures with only four independent variables, unlike existing shell theories, which invariably and implicitly induce a host of unknowns. This approach guarantees traction-free boundary conditions without shear correction factors, using a non-polynomial hyperbolic warping function for transverse shear deformation and stress. By introducing undetermined integral terms, it will be possible to derive the motion equations with a low order of differentiation, which can facilitate a closed-form solution in conjunction with Navier's procedure. The mechanical properties of the CNT reinforcements are modeled to vary smoothly and gradually through the thickness coordinate, exhibiting different distribution patterns. A comparison study is performed to prove the efficacy of the formulated shell theory via obtained results from existing literature. Further numerical investigations are current and comprehensive in detailing the effects of CNT distribution patterns, volume fractions, and geometrical configurations on the fundamental frequencies of CNT-reinforced nanocomposite shells present here. The current shell theory is assumed to serve as a potent conceptual framework for designing reinforced structures and assessing their mechanical behavior.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.7
• /
• pp.2386-2396
• /
• 1996

A numerical simulation has been performed for isothermal and reacting flows in an exisymmetric, bluff-body research combustor. The present formulation is based on the density-weighted averaged Navier-Stokes equations together with a k-epsilon. turbulence model and a modified eddy-breakup combustion model. The PISO algorithm is employed for solution of thel Navier-Stokes system. Comparison between measurements and predictions are made for a centerline axial velocities, location of stagnation points, strength of recirculation zone, and temperature profile. Even though the numerical simulation gives acceptable agreement with experimental data in many respects, the present model is defictient in predicting the recoveryt rate of a central near-wake region, the non-isotropic turbulence effects, and variation of turbulent Schmidt number. Several possible explanations for these discrepancies have been discussed.

• Journal of Korea Water Resources Association
• /
• v.53 no.1
• /
• pp.1-13
• /
• 2020

The most common approach for computing engineering flow problems at high Reynolds number is still the Reynolds-averaged Navier-Stokes (RANS) computations based on turbulence models with wall functions. The recently developed generalized wall functions blending between the wall-limiting viscous and the outer logarithmic relations ensure a smooth transition of flow quantities across two regions. The performances and convergence properties of widely used turbulence models with wall functions that are applicable for turbulence kinetic energy (TKE), turbulent and specific dissipation rates, and eddy viscosity are presented through a series of near wall flow simulations. The present results show that RNG k-𝜖 model should be carefully applied with small tolerance to get the stable solution when the first grid lies in the buffer layer. The standard k-𝜖 and RNG k-𝜖 models are not sensitive to the selection of wall functions for both TKE and eddy viscosity, while the k-ω SST model should be applied together with kL-wall function for TKE and nutUB-wall functions for eddy viscosity to ensure accurate and stable boundary conditions. The applications to a backward-facing step flow at Re=155,000 reveal that the reattachment length is reasonably well predicted on appropriately refined mesh by all turbulence models, except the standard k-𝜖 model which about 13% underestimates the reattachment length regardless of the grid refinement.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.48 no.5
• /
• pp.323-334
• /
• 2020

Since boundary layer transition has a significant impact on the aero-thermodynamic performance of hypersonic flight vehicles, capability of accurate prediction of transition location is essential for design and performance analysis. In this study, γ-Re_θt model is improved to predict transition of hypersonic boundary layers and validated. A coefficient in the production term of the intermittency transport equation that affects the transition onset location is constructed and applied as a function of Mach number, wall temperature, and freestream stagnation temperature based on the similarity numerical solution of compressible boundary layer. To take into account a Mach number dependency of transition onset momentum thickness Reynolds number and transition length, additional correlation equations are determined as function of Mach number and applied to Re_θc and F_length correlations of the baseline model. The suggested model is implemented to a commercial CFD code in consideration of practical use. Analysis of hypersonic flat plate and circular cone boundary layers is carried out by using the model for validation purpose. An improvement of prediction capability with respect to variation of Mach number and unit Reynolds number is identified from the comparison with experimental data.

• Solar Energy
• /
• v.14 no.2
• /
• pp.75-90
• /
• 1994

Thermal energy transport in a two-dimensional horizontal and vertical channel with an isothermal rectangular beam attached to one adiabatic wall is investigated from the numerical solution of Navier-Stokes and energy equations. The solutions have been obtained for dimensionless aspect equations. The solutions have been obtained for dimensionless aspect ratios of beam, H/B=$0.25{sim}4$, Reynolds numbers, Re=$50{\sim}500$ and Grashof numbers, Gr=$0{\sim}5{\times}10^4$. The mean Nusselt number, $\overline{Nu}$ for horizontal and vertical channels shows same value at Gr=0 and increases as Gr increases and decreases as H/B increases at Re=100. $\overline{Nu}$ of vertical channel shows higher in $0.25{\leq}H/B<1.1$ and lower in $1.1{\leq}H/B{\leq}4.0$ than that of horizontal channel at $Gr=10^4$, Re=100. $\overline{Nu}$ of vertical channel shows higher in $0.25{\leq}H/B<1.1$ and lower in $1.1{\leq}H/B=1.0$ than that of horizontal channel at Re=100, $0<Gr{\leq}5{\times}10^4$. A comparison between the experimental and numerical results shows good agreement.

• The Sea:JOURNAL OF THE KOREAN SOCIETY OF OCEANOGRAPHY
• /
• v.24 no.2
• /
• pp.187-207
• /
• 2019

In order to develop a high performance ocean model, we used Julia, a Just-In-Time compile language, and to obtain the solution of the momentum equation, we made the code to solve the Poisson equation by the Successive Over-Relaxation method. And then we made two models to test Julia calculation codes. First, a simple channel form is modeled to test constant source/sink conditions. Second, the simplified Yellow Sea was modeled to test tidal forcing, Coriolis forces, and the effect of vertical eddy diffusivity coefficients. The model has been tested with a total of eight cases in the two scenarios. As a result of the test, the depth-averaged current speed of the three cases in Scenario 1 converged perfectly to the theoretical value, and that showed well a vertical flow velocity gradient due to the bottom friction. Also, the result of Scenario 2 represented well the amphidromic points of Yellow Sea and the tidal characteristics of mid-western and southwestern coast of Korea. Therefore, it is considered that the ocean model using Julia language has developed successfully, this suggests that the ocean model has come to the stage of successful transition from a classical compile language to a Just-In-Time compile language.