• Journal of Hydrogen and New Energy
• /
• v.33 no.1
• /
• pp.47-54
• /
• 2022

In this study, we performed numerical analysis for 1 kWe SOFC stack of internal manifold types according to the different manifold sizes to verify the influence of the flow uniformity into each cell. To simulate the flow phenomena in the stack, the continuity and momentum conservation equations including the standard k-𝜺 turbulent model for the steady-state conditions were applied. From the calculation results, we verified that the pressure drop from inlet pipes to outlet pipes decreased to a log scale as the manifold size increased in the internal manifold types. Also, we found that the flow uniformity increased on an exponential scale as the manifold size increased. In addition, the calculation results showed that the flow uniformity gradually improved as the fuel and oxygen utilization increased.

• Steel and Composite Structures
• /
• v.47 no.6
• /
• pp.743-757
• /
• 2023

In this paper, the dynamic response of a simply-supported composite sandwich plate with a viscoelastic core based on the Golla-Hughes-McTavish (GHM) viscoelastic model is investigated analytically. The formulation is developed using the three-layered sandwich panel theory. Hamilton's principle has been employed to derive the equations of motion. Since classical models, like kelvin-voigt and Maxwell models, cannot express a comprehensive description of the dynamic behavior of viscoelastic material, the GHM method is used to model the viscoelastic core of the plate in this research. The main advantage of the GHM model in comparison with classical models is the consideration of the frequency-dependent characteristic of viscoelastic material. Identification of the material parameters of GHM mini-oscillator terms is an essential procedure in applying the GHM model. In this study, the focus of viscoelastic modeling is on the development of GHM parameters identification. For this purpose, a new method is proposed to find these constants which express frequency-dependent behavior characterization of viscoelastic material. Natural frequencies and loss factors of the sandwich panel based on ESL and three-layered theories in different geometrics are described at 30℃ and 90℃; also, the comparisons show that obtained natural frequencies are grossly overestimated by ESL theory. The argumentations of differences in natural frequencies are also illustrated in detail. The obtained results show that the GHM model presents a more accurate description of the plate's dynamic response by considering the frequency dependency behavior of the viscoelastic core.

• Structural Engineering and Mechanics
• /
• v.84 no.3
• /
• pp.295-310
• /
• 2022

This paper proposes a modified bar simulation method for analyzing the shear lag effect of variable sectional box girder with multiple cells. This theoretical method formulates the equivalent area of stiffening bars and the allocation proportion of shear flows in webs, and re-derives the governing differential equations of bar simulation method. The feasibility of the proposed method is verified by the model test and finite element (FE) analysis of a simply supported multi-cell box girder with constant depth. Subsequently, parametric analysis is conducted to explore the mechanism of shear lag effect of varied sectional cantilever box girder with multiple cells. Results show that the shear lag behavior of variable box-section cantilever box girder is weaker than that of box girder with constant section. It is recommended to make the gradient of shear flow in the web with respect to span length vary as smoothly as possible for eliminating the shear lag effect of box girder. An effective countermeasure for diminishing shear lag effect is to increase the number of box chambers or change the variation manner of bridge depth. The shear lag effect of varied sectional cantilever box girder will get more server when the length of central flanges is shorter than 0.26 or longer than 0.36 times of total width of top flange, as well as the cantilever length exceeds 0.29 times of total length of box's flange. Therefore, the distance between central webs can adjust the shear lag effect of box girder. Especially, the width ratio of cantilever plate with respect to total length of top flange is proposed to be no more 1/3.

• Steel and Composite Structures
• /
• v.44 no.5
• /
• pp.707-720
• /
• 2022

In the present work, bending and free vibration analyses of multilayered functionally graded (FG) graphene platelet (GPL) and fiber-reinforced hybrid composite beams are carried out using the parabolic function based shear deformation theory. Parabolic variation of transverse shear stress across the thickness of beam and transverse shear stress-free conditions at top and bottom surfaces of the beam are considered, and the proposed formulation incorporates a transverse displacement field. The present theory works only with four unknowns and is computationally efficient. Hamilton's principle has been employed for deriving the governing equations. Analytical solutions are obtained for both the bending and free vibration problems in the present work considering different variations of GPLs and fibers distribution, namely, FG-X, FG-U, FG-Λ, and FG-O for beams having simply-supported boundary condition. First, the matrix is assumed to be strengthened using GPLs, and then the fibers are embedded. Multiscale modeling for material properties of functionally graded graphene platelet/fiber hybrid composites (FG-GPL/FHRC) is performed using Halpin-Tsai micromechanical model. The study reveals that the distributions of GPLs and fibers have significant impacts on the stresses, deflections, and natural frequencies of the beam. The number of layers and shape factors widely affect the behavior of FG-GPL-FHRC beams. The multilayered FG-GPL-FHRC beams turn out to be a good approximation to the FG beams without exhibiting the stress-channeling effects.

• The Journal of the Acoustical Society of Korea
• /
• v.27 no.3
• /
• pp.129-139
• /
• 2008

The purpose of the present work is to derive the reflection and transmission coefficients of $S_0\;and\;A_0$ mode Lamb waves in relation to the geometry of a rectangular notch when the waves propagate across the notch in an elastic plate. Firstly, the excitable modes of the Lamb wave were analyzed with respect to the plate thickness. The scattering phenomena were divided into three independent processes according to the boundary shape of the notch and the direction of the wave propagation. Linear equations for each process were derived with corresponding free or continuous boundary conditions to analyze the scattered waves. By the rule of linear superposition, the waves scattered at each process were summed for each mode. Then the steady-state reflection and transmission coefficients of the scattered waves were determined so that the difference of energy flux between the incident and the scattered waves would remain within 4%.

• Advances in nano research
• /
• v.16 no.2
• /
• pp.111-125
• /
• 2024

Recently, a lot of research has been done on the analysis of axial vibrations of homogeneous and FG nanotubes (nanorods) with various aspects of vibrations that have been fully mentioned in history. However, there is a lack of investigation of the dynamic internal resonances of FG nanotubes (nanorods) between them. This is one of the essential or substantial characteristics of nonlinear vibration systems that have many applications in various fields of engineering (making actuators, sensors, etc.) and medicine (improving the course of diseases such as cancers, etc.). For this reason, in this study, for the first time, the dynamic internal resonances of FG nanorods in the simultaneous presence of large-amplitude size dependent behaviour, inertial and shear effects are investigated for general state in detail. Such theoretical patterns permit as to carry out various numerical experiments, which is the key point in the expansion of advanced nano-devices in different sciences. This research presents an AFG novel nano resonator model based on the axial vibration of the elastic nanorod system in terms of derivation from large-amplitude size dependent internal modals interactions. The Hamilton's Principle is applied to achieve the basic equations in movement and boundary conditions, and a harmonic deferential quadrature method, and a multiple scale solution technique are employed to determine a semi-analytical solution. The interest of the current solution is seen in its specific procedure that useful for deriving general relationships of internal resonances of FG nanorods. The numerical results predicted by the presented formulation are compared with results already published in the literature to indicate the precision and efficiency of the used theory and method. The influences of gradient index, aspect ratio of FG nanorod, mode number, nonlinear effects, and nonlocal effects variations on the mechanical behavior of FG nanorods are examined and discussed in detail. Also, the inertial and shear traces on the formations of internal resonances of FG nanorods are studied, simultaneously. The obtained valid results of this research can be useful and practical as input data of experimental works and construction of devices related to axial vibrations of FG nanorods.

• Geomechanics and Engineering
• /
• v.37 no.1
• /
• pp.49-64
• /
• 2024

Rocks undergoing repeated loading and unloading over an extended period, such as due to earthquakes, human excavation, and blasting, may result in the gradual accumulation of stress and deformation within the rock mass, eventually reaching an unstable state. In this study, a CNN-CCM is proposed to address the mechanical behavior. The structure and hyperparameters of CNN-CCM include Conv2D layers × 5; Max pooling2D layers × 4; Dense layers × 4; learning rate=0.001; Epoch=50; Batch size=64; Dropout=0.5. Training and validation data for deep learning include 71 rock samples and 122,152 data points. The AI Rock Constitutive Model learned by CNN-CCM can predict strain values(ε₁) using Mass (M), Axial stress (σ₁), Density (ρ), Cyclic number (N), Confining pressure (σ₃), and Young's modulus (E). Five evaluation indicators R², MAPE, RMSE, MSE, and MAE yield respective values of 0.929, 16.44%, 0.954, 0.913, and 0.542, illustrating good predictive performance and generalization ability of model. Finally, interpreting the AI Rock Constitutive Model using the SHAP explaining method reveals that feature importance follows the order N > M > σ₁ > E > ρ > σ₃.Positive SHAP values indicate positive effects on predicting strain ε₁ for N, M, σ₁, and σ₃, while negative SHAP values have negative effects. For E, a positive value has a negative effect on predicting strain ε₁, consistent with the influence patterns of conventional physical rock constitutive equations. The present study offers a novel approach to the investigation of the mechanical constitutive model of rocks under cyclic loading and unloading conditions.

• Earthquakes and Structures
• /
• v.27 no.2
• /
• pp.127-142
• /
• 2024

The amount of wave propagation through a rubber concrete construction is the subject of the current investigation. Rubber tire waste was used to make two different types of cement mixtures. One type contains sand substitute in amounts ranging from 15% to 60% of the total volume, while the other has gravel with diameters of 3/8 and 8/15 and 15% sand in the same mixture. A wide variety of concrete forms and compositions were created, and their viscous and solid state characteristics were assessed, along with their short-, medium-, and long-term strengths. Diffusion, density, mechanical strength resistance to compressive force, and ultrasound wave propagation were also assessed. The water-to-cement ratio and plasticizer were used in this investigation. In the second part of the study, an analytical model is presented that simulates the experimental model in predicting the speed of waves and the frequencies accompanying them for this type of mixture. Higher order shear deformation beam theory for wave propagation in the rubberized concrete beam is developed, considering the bidirectional distribution, which is primarily expressed by the density, the Poisson coefficient, and Young's modulus. Hamilton's concept is used to determine the governing equations of the wave propagation in the rubberized concrete beam structure. When the analytical and experimental results for rubber concrete beams were compared, the outcomes were very comparable. The addition of rubber gravel and sandy rubber to the mixture both resulted in a discernible drop in velocities and frequencies, according to the data.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2003.05a
• /
• pp.91-93
• /
• 2003

A comprehensive numerical study is carried out to investigate for the understanding of the flow evolution and flame development in a supersonic combustor with normal injection of ncumally injecting hydrogen in airsupersonic flows. The formulation treats the complete conservation equations of mass, momentum, energy, and species concentration for a multi-component chemically reacting system. For the numerical simulation of supersonic combustion, multi-species Navier-Stokes equations and detailed chemistry of H2-Air is considered. It also accommodates a finite-rate chemical kinetics mechanism of hydrogen-air combustion GRI-Mech. 2.11[1], which consists of nine species and twenty-five reaction steps. Turbulence closure is achieved by means of a k-two-equation model (2). The governing equations are spatially discretized using a finite-volume approach, and temporally integrated by means of a second-order accurate implicit scheme (3-5).The supersonic combustor consists of a flat channel of 10 cm height and a fuel-injection slit of 0.1 cm width located at 10 cm downstream of the inlet. A cavity of 5 cm height and 20 cm width is installed at 15 cm downstream of the injection slit. A total of 936160 grids are used for the main-combustor flow passage, and 159161 grids for the cavity. The grids are clustered in the flow direction near the fuel injector and cavity, as well as in the vertical direction near the bottom wall. The no-slip and adiabatic conditions are assumed throughout the entire wall boundary. As a specific example, the inflow Mach number is assumed to be 3, and the temperature and pressure are 600 K and 0.1 MPa, respectively. Gaseous hydrogen at a temperature of 151.5 K is injected normal to the wall from a choked injector.A series of calculations were carried out by varying the fuel injection pressure from 0.5 to 1.5MPa. This amounts to changing the fuel mass flow rate or the overall equivalence ratio for different operating regimes. Figure 1 shows the instantaneous temperature fields in the supersonic combustor at four different conditions. The dark blue region represents the hot burned gases. At the fuel injection pressure of 0.5 MPa, the flame is stably anchored, but the flow field exhibits a high-amplitude oscillation. At the fuel injection pressure of 1.0 MPa, the Mach reflection occurs ahead of the injector. The interaction between the incoming air and the injection flow becomes much more complex, and the fuel/air mixing is strongly enhanced. The Mach reflection oscillates and results in a strong fluctuation in the combustor wall pressure. At the fuel injection pressure of 1.5MPa, the flow inside the combustor becomes nearly choked and the Mach reflection is displaced forward. The leading shock wave moves slowly toward the inlet, and eventually causes the combustor-upstart due to the thermal choking. The cavity appears to play a secondary role in driving the flow unsteadiness, in spite of its influence on the fuel/air mixing and flame evolution. Further investigation is necessary on this issue. The present study features detailed resolution of the flow and flame dynamics in the combustor, which was not typically available in most of the previous works. In particular, the oscillatory flow characteristics are captured at a scale sufficient to identify the underlying physical mechanisms. Much of the flow unsteadiness is not related to the cavity, but rather to the intrinsic unsteadiness in the flowfield, as also shown experimentally by Ben-Yakar et al. [6], The interactions between the unsteady flow and flame evolution may cause a large excursion of flow oscillation. The work appears to be the first of its kind in the numerical study of combustion oscillations in a supersonic combustor, although a similar phenomenon was previously reported experimentally. A more comprehensive discussion will be given in the final paper presented at the colloquium.

• Anatomy & Biological Anthropology
• /
• v.31 no.4
• /
• pp.121-131
• /
• 2018

Stature is one of the unique biological properties of a person, which can be used for identification of the individual. In this regard, statures are estimated for the unknown victims from crimes and disasters. However, the accuracy of estimates may be compromised by inappropriate methodologies and/or practices of stature estimation. Discussed in this study are the methodological issues related to the current practices of forensic anthropological stature estimation in Korea, followed by suggestions to enhance the accuracy of the stature estimates. Summaries of forensic anthropological examinations for 560 skeletal remains, which were conducted at the National Forensic Service (NFS), were reviewed. Mr. Yoo Byung-eun's case is utilized as an example of the NFS's practices. To estimate Mr. Yoo's stature, Trotter's (1970) femur equation was applied even though the fibula equation of a lower standard error was available. In his case report, the standard error associated with the equation (${\pm}3.8cm$) was interpreted as an 'error range', which gave a hasty impression that the prediction interval is that narrow. Also, stature shrinkage by aging was not considered, so the estimated stature in Mr. Yoo's case report should be regarded as his maximum living stature, rather than his stature-at-death. Lastly, applying Trotter's (1970) White female equations to Korean female remains is likely to underestimate their statures. The anatomical method will enhance the accuracy of stature estimates. However, in cases that the anatomical method is not feasible, the mathematical method based on Korean samples should be considered. Since 1980's, effort has been made to generate stature estimation equations using Korean samples. Applying the equations based on Korean samples to Korean skeletal remains will enhance the accuracy of the stature estimates, which will eventually increase the likelihood of successful identification of the unknowns.