• Coupled systems mechanics
• /
• v.6 no.4
• /
• pp.439-464
• /
• 2017

The hydro-elastic system consisting of a pre-stretched highly elastic plate, compressible Newtonian viscous fluid, and the rigid wall is considered and it is assumed that on the plate a lineal-located time-harmonic force acts. It is required to investigate the dynamic behavior of this system and determine how the problem parameters and especially the pre-straining of the plate acts on this behavior. The elasticity relations of the plate are described through the harmonic potential and linearized (with respect to perturbations caused by external time-harmonic force) form of these relations is used in the present investigation. The plane-strain state in the plate is considered and the motion of that is described within the scope of the three-dimensional linearized equations of elastic waves in elastic bodies with initial stresses. The motion of the fluid is described by the linearized Navier-Stokes equations and it is considered the plane-parallel flow of this fluid. The Fourier transform with respect to the space coordinate is applied for a solution to the corresponding boundary-value problem. Numerical results on the frequency response of the interface normal stress and normal velocity and the influence of the initial stretching of the plate on this response are presented and discussed. In particular, it is established that the initial stretching of the plate can decrease significantly the absolute values of the aforementioned quantities.

• Journal of Wetlands Research
• /
• v.14 no.1
• /
• pp.75-87
• /
• 2012

Hyporheic zone is a region beneath and alongside a stream, river, or lake bed, where there is mixing of shallow groundwater and surfacewater. Hyporheic exchange controls a variety of physical, biogeochemical and thermal processes, and provides unique ecotones in a aquatic ecosystem. Field and experimental observations, and modeling studies indicate that hyporheic exchange is mainly in response to pressure gradients driven by the geomorphological features of stream beds. In the reach scale of a stream, pool-riffle structures dominate the exchange patterns. Flow over a pool-riffle sequence develops recirculation zones and stagnation points, and this flow structures make irregular pressure gradient which is driving force of the hyporheic exchange. In this study, 3 D hydro-dynamic model solves the Reynolds-averaged Navier-Stokes equations for the surface water and Darcy's Law and the continuity equation for ground water. The two sets of equations are coupled via the pressure distribution along the interface. Simulation results show that recirculation zones and stagnation points in the pool-riffle structures dominantly control the upwelling and downwelling patterns. With decrease of recirculation zones, length of donwelling zone formed in front of riffles is reduced and position of maximum downwelling point moves downward. The numerical simulation could successfully predict the behavior of hyporheic exchange and contribute the field study, river management and restoration.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.36 no.1
• /
• pp.72-78
• /
• 2008

In this study, advanced (fluid-structure interaction (FSI)) analysis system has been developed in order to predict turbine cascade performance with blade deformation effect due to aerodynamic loads. Intereference effects due to the relative movement of the rotor cascade with respect to the stator cascade are also considered. Reynolds-averaged Navier-Stokes equations with one equation Spalart-Allmaras and two-equation k-ω SST turbulence models are solved to accurately predict fluid dynamic loads considering flow separation effects. A fully implicit time marching scheme based on the (coupled Newmark time-integration method) with high artificial damping is efficiently used to compute the complex fluid-structure interaction problem. Predicted aerodynamic performance considering structural deformation effect of the blade shows somewhat different results compared to the case of rigid blade model. Cascade performance evaluations for different elastic axis positions are importantly presented and its aeroelastic effects are investigated.

• International Journal of Aeronautical and Space Sciences
• /
• v.18 no.2
• /
• pp.186-196
• /
• 2017

In order to resolve the trimming difficulty in rotor CFD calculations, a high-efficiency and improved "delta trim method" is established to compute the blade control settings that are necessary to identify the blade motion. In this method, a simplified model which combines the blade element theory and different inflow models is employed to calculate the control settings according to the target aerodynamic forces, then it is coupled into a CFD solver with unsteady Navier-Stokes equations by the delta methodology, which makes the control settings and aerodynamics calculated and updated in the meantime at every trim cycle. Different from the previous work, the current research combines the inflow model based on prescribed wake theory. Using the method established, the control settings and aerodynamic characteristics of Helishape 7A, AH-1G and Caradonna-Tung rotors are calculated. The influence of different inflow models on trimming calculations is analyzed and the computational efficiency of the current "delta trim method" is compared with that of the "CFD-based trim method". Furthermore, for the sake of improving the calculation efficiency, a novel acceleration factor method is introduced to accelerate the trimming process. From the numerical cases, it is demonstrated that the current "delta trim method" has higher computational efficiency than "CFD-based trim method" in both hover and forward flight, and up to 70% of the amount of calculation can be saved by current "delta trim method" which turns out to be satisfactory for engineering applications. In addition, the proposed acceleration factor shows a good ability to accelerate the trim procedure, and the prescribed wake inflow model is always of better stability than other simple inflow models whether the acceleration factor is utilized in trimming calculations.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.7
• /
• pp.1-8
• /
• 2005

For study on the unsteady wall interference effect, flows around a circular cylinder in closed test-section wind tunnels have been numerically investigated by solving compressible Navier-Stokes equations. The numerical scheme is based on a node-based finite-volume method with the Roe's flux-difference splitting and an implicit time-integration method coupled with dual time-step sub-iteration. The computed results showed that the unsteady pressure gradient over the cylinder is enhanced by the wall interference, and as a result the fluctuations of lift and drag are augmented. The drag is further increased because of the lower base pressure. The vortex shedding frequency is also increased by the wall interference. The pressure on the test section wall shows the harmonics having the shedding frequency contained in the wall effect.

• Ocean Systems Engineering
• /
• v.7 no.4
• /
• pp.435-460
• /
• 2017

The main objective of this study is to investigate the turning and zig-zag maneuvering performance of the well-known naval surface combatant DTMB (David Taylor Model Basin) 5415 hull with URANS (Unsteady Reynolds-averaged Navier-Stokes) method. Numerical simulations of static drift tests have been performed by a commercial RANS solver based on a finite volume method (FVM) in an unsteady manner. The fluid flow is considered as 3-D, incompressible and fully turbulent. Hydrodynamic analyses have been carried out for a fixed Froude number 0.28. During the analyses, the free surface effects have been taken into account using VOF (Volume of Fluid) method and the hull is considered as fixed. First, the code has been validated with the available experimental data in literature. After validation, static drift, static rudder and drift and rudder tests have been simulated. The forces and moments acting on the hull have been computed with URANS approach. Numerical results have been applied to determine the hydrodynamic maneuvering coefficients, such as, velocity terms and rudder terms. The acceleration, angular velocity and cross-coupled terms have been taken from the available experimental data. A computer program has been developed to apply a fast maneuvering simulation technique. Abkowitz's non-linear mathematical model has been used to calculate the forces and moment acting on the hull during the maneuvering motion. Euler method on the other hand has been applied to solve the simultaneous differential equations. Turning and zig-zag maneuvering simulations have been carried out and the maneuvering characteristics have been determined and the numerical simulation results have been compared with the available data in literature. In addition, viscous effects have been investigated using Eulerian approach for several static drift cases.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.35 no.10
• /
• pp.1053-1060
• /
• 2011

This study is performed to determine the drug concentration profiles of drug-eluting stents (DES) for an ideal circular ring stent and intertwined stent models for various Reynolds numbers (Re = 200, 400, and 800). The Navier.Stokes equations coupled with the advection-diffusion equation are solved numerically in order to determine how the flow patterns and drug deposition are affected in the in-stent and post-stent regions where flow separation and recirculation occur. The presence of DES within the arterial segment affects the local drug distribution in the flow field. As a result, the drug concentration for the intertwined stent is higher over the in-stent region in comparison with the ideal stents. For a given stent geometry, the local drug concentration in the in-stent region decreases with Reynolds number, while for a given Reynolds number, the local drug concentration is relatively insensitive to the stent geometry. The results show that drug concentration along the arterial wall is significantly higher within the in-stent and post-stent regions for the intertwined stent geometry than for the ideal stent geometries.

• Current Optics and Photonics
• /
• v.2 no.3
• /
• pp.261-268
• /
• 2018

Stimulated polariton scattering (SPS) from the $A_1$ transverse optical (TO) modes of $BaTiO_3$ bulk crystal generating a terahertz (THz) wave with the noncollinear phase-matching (NPM) condition is theoretically investigated. To our best knowledge, this is the first report on THz wave generation from $BaTiO_3$ bulk crystal via SPS. Phase-matching (PM) characteristics in the NPM configuration are analyzed. Effective parametric gain lengths for the Stokes and THz waves in the NPM configuration are calculated. The effective parametric gain coefficient and absorption coefficient of the THz wave in $BaTiO_3$ are theoretically simulated. The THz phonon flux densities generated via SPS in $BaTiO_3$ are theoretically calculated by solving the coupled wave equations under the NPM condition. The PM characteristics and THz-wave parametric gain characteristics in $BaTiO_3$ are compared to those in $MgO:LiNbO_3$. The results of the analysis indicate that $BaTiO_3$ is an attractive optical crystal for efficient THz wave generation via SPS.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.11 no.2
• /
• pp.262-270
• /
• 1987

A numerical method is presented which can solve the unsteady momentum and thermal boundary layers, coupled through the agency of buoyancy force, over a heated circular cylinder impulsively started from rest. By linearizing the nonlinear finite difference equations without sacrificing accuracy, numerical solutions are obtained at each time step without iteration. To get rid of the requirement of excessive number of grid points in the region of reversed flow, special form of transformed variables are used, by which the computational boundary layer thickness is maintained almost constant. These numerical properties enable the method to easily handle the region of reversed flow and how the singularity develops in the interior of the boundary layer. In order to investigated the thermal effects on the skin friction, heat flux, displacement thickness and on the separation, we have successfully solved three different cases of the buoyancy parameter .alpha.(Gr/Re$^{2}$).

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.18 no.3
• /
• pp.189-197
• /
• 2006

The influence of leodo Ocean Research Station structure to surrounding atmospheric flow is carefully investigated using CFD techniques. Moreover, the validation works of computational results are performed by the comparison with the observed data of leodo Ocean Research station. In this paper, we performed 3-dimensional CAD modelling of the station, generated the grid system for numerical analysis and carried out flow analyses using Navier-Stokes equations coupled with two-equation turbulence model. For suitable free stream conditions of wind speed and direction, the interference of the research station structure on the flow field is predicted. Beside, the computational results are benchmarked by observed data to confirm the accuracy of measured date and reliable data range of each measuring position according to the wind direction. Through the results of this research, now the quantitative evaluation of the error range of interfered gauge data is possible, which is expected to be applied to provide base data of accurate sea surface wind around research stations.