• Advances in nano research
• /
• v.16 no.4
• /
• pp.413-426
• /
• 2024

This paper investigates the dynamic behavior of a simply supported viscoelastic plate made of functionally graded carbon nanotube reinforced composite under dynamic loading. Carbon nanotubes are distributed in 5 different shapes: U, V, A, O and X, depending on the shape they form through the thickness of the plate. The displacement fields are derived in the Laplace domain using a higher-order shear deformation theory. Equations of motion are obtained through the application of the energy method and Hamilton's principle. The resulting equations of motion are solved using Navier's method. Transforming the Laplace domain displacements into the time domain involves Durbin's modified inverse Laplace transform. To validate the accuracy of the developed algorithm, a free vibration analysis is conducted for simply supported plate made of functionally graded carbon nanotube reinforced composite and compared against existing literature. Subsequently, a parametric forced vibration analysis considers the influence of various parameters: volume fractions of carbon nanotubes, their distributions, and ratios of instantaneous value to retardation time in the relaxation function, using a linear standard viscoelastic model. In the forced vibration analysis, the dynamic distributed load applied to functionally graded carbon nanotube reinforced composite viscoelastic plate is obtained in terms of double trigonometric series. The study culminates in an examination of maximum displacement, exploring the effects of different carbon nanotube distributions, volume fractions, and ratios of instantaneous value to retardation times in the relaxation function on the amplitudes of maximum displacements.

• 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.

• Wind and Structures
• /
• v.20 no.4
• /
• pp.549-564
• /
• 2015

The three-dimensional unsteady incompressible Reynolds-averaged Navier-Stokes equations and k-${\varepsilon}$ double equations turbulent model were used to investigate the effect on the measurements of anemometers due to a passing high-speed train. Sliding mesh technology in Fluent was utilized to treat the moving boundary problem. The high-speed train considered in this paper was with bogies and inter-carriage gaps. Combined with the results of the wind tunnel test in a published paper, the accuracy of the present numerical method was validated to be used for further study. In addition, the difference of slipstream between three-car and eight-car grouping models was analyzed, and a series of numerical simulations were carried out to study the influences of the anemometer heights, the train speeds, the crosswind speeds and the directions of the induced slipstream on the measurements of the anemometers. The results show that the influence factors of the train-induced slipstream are the passing head car and tail car. Using the three-car grouping model to analyze the train-induced flow is reasonable. The maxima of horizontal slipstream velocity tend to reduce as the height of the anemometer increases. With the train speed increasing, the relationship between $V_{train}$ and $V_{induced\;slipstream}$ can be expressed with linear increment. In the absence of natural wind conditions, from the head car arriving to the tail car leaving, the induced wind direction changes about $330^{\circ}$, while under the crosswind condition the wind direction fluctuates around $-90^{\circ}$. With the crosswind speed increasing, the peaks of $V_X,{\mid}V_{XY}-V_{wind}{\mid}$ of the head car and that of $V_X$ of the tail car tend to enlarge. Thus, when anemometers are installed along high-speed railways, it is important to study the effect on the measurements of anemometers due to the train-induced slipstream.

• 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.