• Coupled systems mechanics
• /
• v.6 no.3
• /
• pp.287-316
• /
• 2017

This paper studies the particularities of the forced vibration of the hydro-elastic system consisting of a moving elastic plate, compressible viscous fluid and rigid wall. This study is made by employing the discrete-analytical solution method proposed in the paper by the authors (Akbarov and Panakhli (2015)). It is assumed that in the initial state the fluid flow is caused by the axial movement of the plate and the additional lineally-located time-harmonic forces act on the plate and these forces cause additional flow field in the fluid and a stress-strain state in the plate. The stress-strain state in the plate is described by utilizing the exact equations and relations of the linear elastodynamics. However, the additional fluid flow field is described with linearized Navier-Stokes equations for a compressible viscous fluid. Numerical results related to the influence of the problem parameters on the frequency response of the normal stress acting on the plate fluid interface plane and fluid flow velocity on this plane are presented and discussed. In this discussion, attention is focused on the influence of the initial plate axial moving velocity on these responses. At the same, it is established that as a result of the plate moving a resonance type of phenomenon can take place under forced vibration of the system. Moreover, numerical results regarding the influence of the fluid compressibility on these responses are also presented and discussed.

• Journal of computational fluids engineering
• /
• v.22 no.1
• /
• pp.59-66
• /
• 2017

The guideline of selecting the number of snapshot dataset, $N_s$ in proper orthogonal decomposition(POD) was presented via the analysis of Eigen values based on the singular value decomposition(SVD). In POD, snapshot datasets from the solutions of Euler or Navier-Stokes equations are utilized to SVD and a reduced order model(ROM) is constructed as the combination of Eigen vectors. The ROM is subsequently applied to reconstruct the flowfield data with new set of flow conditions, thereby enhancing the computational efficiency. The overall computational efficiency and accuracy of POD is dependent on the number of snapshot dataset; however, there is no reliable guideline of determining $N_s$. In order to resolve this problem, the order of maximum to minimum Eigen value ratio, O(R) from SVD was analyzed and presented for the decision of $N_s$; in case of steady flow, $N_s$ should be determined to make O(R) be $10^9$. For unsteady flow, $N_s$ should be increased to make O(R) be $10^{11\sim12}$. This strategy of selecting the snapshot dataset was applied to two dimensional NACA0012 airfoil and vortex flow problems including steady and unsteady cases and the numerical accuracies according to $N_s$ and O(R) were discussed.

• The KSFM Journal of Fluid Machinery
• /
• v.13 no.5
• /
• pp.43-53
• /
• 2010

This study investigates a design optimization of a rotating two-pass rectangular cooling channel with staggered arrays of pin-fins. The radial basis neural network method is used as an optimization technique with Reynolds-averaged Navier-Stokes analysis of fluid flow and heat transfer with shear stress transport turbulent model. The ratio of the diameter to height of the pin-fins and the ratio of the streamwise spacing between the pin-fins to height of the pin-fin are selected as design variables. The optimization problem has been defined as a minimization of the objective function, which is defined as a linear combination of heat transfer related term and friction loss related term with a weighting factor. Results are presented for streamlines, velocity vector fields, and contours of Nusselt numbers, friction coefficients, and turbulent kinetic energy. These results show how fluid flow in a two-pass square cooling channel evolves a converted secondary flows due to Coriolis force, staggered arrays of pin-fins, and a $180^{\circ}$ turn region. These results describe how the fluid flow affects surface heat transfer. The Coriolis force induces heat transfer discrepancy between leading and trailing surfaces, having higher Nusselt number on the leading surface in the second pass while having lower Nusselt number on the trailing surface. Dean vortices generated in $180^{\circ}$ turn region augment heat transfer in the turning region and in the upstream region of the second pass. As the result of optimization, in comparison with the reference geometry, thermal performance of the optimum geometry shows the improvement by 30.5%. Through the optimization, the diameter of pin-fin increased by 14.9% and the streamwise distance between pin-fins increased by 32.1%. And, the value of objective function decreased by 18.1%.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.2 no.4
• /
• pp.208-222
• /
• 1990

The nonlinear hydrodynamic pressure force and moment on hinged wavemakers of variable-draft are presented. A closed-form solution (correct to second-order) for the nonlinear wavemaker boundary value problem has been obtained by employing the Stokes perturbation expansion scheme. The physical significance of the second-order contributions to the hydrodynamic pressure moment are examined in detail. Design curves are presented which demonstrate both the magnitude of the second-order nonlinearities and the effects of the variable-draft hinge height. The second-order contributions to the total hydrodynamic force and moment consist of a time-dependent and a steady part. The sum of the first and second-order pressure force and moment show a significant increase over those predicted by linear wavemaker theory. The second-order effects are shown to vary with both relative water depth and wave amplitude. The second-order dynamic effects are relatively more important for hinged wavemakers with shallower drafts.

• Journal of Ship and Ocean Technology
• /
• v.12 no.1
• /
• pp.16-27
• /
• 2008

The interaction between advancing ships and the waves generated by them plays important roles in wave resistances and ship motions. Wave breaking phenomena near the ship bow at different speeds are investigated both numerically and experimentally. Numerical simulations of free surface profiles near the fore bodies of ships are performed and visualized to grasp the general trend or the mechanism of wave breaking phenomena from moderate waves rather than concentrating on local chaotic irregularities as ship speeds increase. Navier-Stokes equations are differentiated based on the finite difference method. The Marker and Cell (MAC) Method and Marker-Density Method are employed, and they are compared for the description of free surface conditions associated with the governing equations. Extra effort has been directed toward the realization of extremely complex free surface conditions at wave breaking. For this purpose, the air-water interface is treated with marker density, which is used for two layer flows of fluids with different properties. Adaptation schemes and refinement of the numerical grid system are also used at local complex flows to improve the accuracy of the solutions. In addition to numerical simulations, various model tests are performed in a ship model towing tank. The results are compared with numerical calculations for verification and for realizing better, more efficient research performance. It is expected that the present research results regarding wave breaking and the geometry of the fore body of ship will facilitate better hull form design productivity at the preliminary ship design stage, especially in the case of small and fast ship design. Also, the obtained knowledge on the impact due to the interaction of breaking waves and an advancing hull surface is expected to be applicable to investigation of the ship bow slamming problem as a specific application.

• Journal of Ocean Engineering and Technology
• /
• v.31 no.4
• /
• pp.274-280
• /
• 2017

Generally, global performance analysis in offshore platforms is performed using potential-based numerical tools, which neglect hydrodynamic viscous effects. In comparison with the potential theory, computational fluid dynamics (CFD) methods can take into account the viscous effects by solving the Navier-Stokes equation using the finite-volume method. The open-source field operation and manipulation (OpenFOAM) C++ libraries are employed for a finite volume method (FVM) numerical analysis. In this study, in order to apply CFD to the global performance analysis of a hull-mooring coupled system, we developed a numerical wave basin to analyze the global performance problem of a floating body with a catenary mooring system under regular wave conditions. The mooring system was modeled using a catenary equation and solved in a quasi-static condition, which excluded the dynamics of the mooring lines such as the inertia and drag effects. To demonstrate the capability of the numerical basin, the global performance of a barge with four mooring lines was simulated under regular wave conditions. The simulation results were compared to the analysis results from a commercial mooring analysis program, Orcaflex. The comparison included the motion of the barge, catenary shape, and tension in the mooring lines. The study found good agreement between the results from the developed CFD-based numerical calculation and commercial software.

• Journal of Ocean Engineering and Technology
• /
• v.34 no.6
• /
• pp.406-418
• /
• 2020

In this study, a nonlinear wave simulation code is developed using a higher-order spectral (HOS) method. The HOS method is very efficient because it can determine the solution of the boundary value problem using fast Fourier transform (FFT) without matrix operation. Based on the HOS order, the vertical velocity of the free surface boundary was estimated and applied to the nonlinear free surface boundary condition. Time integration was carried out using the fourth order Runge-Kutta method, which is known to be stable for nonlinear free-surface problems. Numerical stability against the aliasing effect was guaranteed by using the zero-padding method. In addition to simulating the initial wave field distribution, a nonlinear adjusted region for wave generation and a damping region for wave absorption were introduced for wave generation simulation. To validate the developed simulation code, the adjusted simulation was carried out and its results were compared to the eighth order Stokes theory. Long-time simulations were carried out on the irregular wave field distribution, and nonlinear wave propagation characteristics were observed from the results of the simulations. Nonlinear adjusted and damping regions were introduced to implement a numerical wave tank that successfully generated nonlinear regular waves. According to the variation in the mean wave steepness, irregular wave simulations were carried out in the numerical wave tank. The simulation results indicated an increase in the nonlinear interaction between the wave components, which was numerically verified as the mean wave steepness. The results of this study demonstrate that the HOS method is an accurate and efficient method for predicting the nonlinear interaction between waves, which increases with wave steepness.

• Nuclear Engineering and Technology
• /
• v.54 no.1
• /
• pp.36-48
• /
• 2022

As a valid numerical method to obtain a high-resolution result of a flow field, computational fluid dynamics (CFD) have been widely used to study coolant flow and heat transfer characteristics in fuel rod bundles. However, the time-consuming, iterative calculation of Navier-Stokes equations makes CFD unsuitable for the scenarios that require efficient simulation such as sensitivity analysis and uncertainty quantification. To solve this problem, a reduced-order model (ROM) based on proper orthogonal decomposition (POD) and machine learning (ML) is proposed to simulate the flow field efficiently. Firstly, a validated CFD model to output the flow field data set of the rod bundle is established. Secondly, based on the POD method, the modes and corresponding coefficients of the flow field were extracted. Then, an deep feed-forward neural network, due to its efficiency in approximating arbitrary functions and its ability to handle high-dimensional and strong nonlinear problems, is selected to build a model that maps the non-linear relationship between the mode coefficients and the boundary conditions. A trained surrogate model for modes coefficients prediction is obtained after a certain number of training iterations. Finally, the flow field is reconstructed by combining the product of the POD basis and coefficients. Based on the test dataset, an evaluation of the ROM is carried out. The evaluation results show that the proposed POD-ROM accurately describe the flow status of the fluid field in rod bundles with high resolution in only a few milliseconds.

• Wind and Structures
• /
• v.34 no.1
• /
• pp.43-58
• /
• 2022

The BARC flow is studied via Direct Numerical Simulation at a relatively low turbulent Reynolds number, with focus on the geometrical representation of the leading-edge (LE) corners. The study contributes to further our understanding of the discrepancies between existing numerical and experimental BARC data. In a first part, rounded LE corners with small curvature radii are considered. Results show that a small amount of rounding does not lead to abrupt changes of the mean fields, but that the effects increase with the curvature radius. The shear layer separates from the rounded LE at a lower angle, which reduces the size of the main recirculating region over the cylinder side. In contrast, the longitudinal size of the recirculating region behind the trailing edge (TE) increases, as the TE shear layer is accelerated. The effect of the curvature radii on the turbulent kinetic energy and on its production, dissipation and transport are addressed. The present results should be contrasted with the recent work of Rocchio et al. (2020), who found via implicit Large-Eddy Simulations at larger Reynolds numbers that even a small curvature radius leads to significant changes of the mean flow. In a second part, the LE corners are fully sharp and the exact analytical solution of the Stokes problem in the neighbourhood of the corners is used to locally restore the solution accuracy degraded by the singularity. Changes in the mean flow reveal that the analytical correction leads to streamlines that better follow the corners. The flow separates from the LE with a lower angle, resulting in a slightly smaller recirculating region. The corner-correction approach is valuable in general, and is expected to help developing high-quality numerical simulations at the high Reynolds numbers typical of the experiments with reasonable meshing requirements.

• Nuclear Engineering and Technology
• /
• v.55 no.10
• /
• pp.3874-3897
• /
• 2023

A high-fidelity computational fluid dynamics (CFD) analysis was performed using the Large Eddy Simulation (LES) model for the lower plenum of the High-Temperature Test Facility (HTTF), a ¼ scale test facility of the modular high temperature gas-cooled reactor (MHTGR) managed by Oregon State University. In most next-generation nuclear reactors, thermal stress due to thermal striping is one of the risks to be curiously considered. This is also true for HTGRs, especially since the exhaust helium gas temperature is high. In order to evaluate these risks and performance, organizations in the United States led by the OECD NEA are conducting a thermal hydraulic code benchmark for HTGR, and the test facility used for this benchmark is HTTF. HTTF can perform experiments in both normal and accident situations and provide high-quality experimental data. However, it is difficult to provide sufficient data for benchmarking through experiments, and there is a problem with the reliability of CFD analysis results based on Reynolds-averaged Navier-Stokes to analyze thermal hydraulic behavior without verification. To solve this problem, high-fidelity 3-D CFD analysis was performed using the LES model for HTTF. It was also verified that the LES model can properly simulate this jet mixing phenomenon via a unit cell test that provides experimental information. As a result of CFD analysis, the lower the dependency of the sub-grid scale model, the closer to the actual analysis result. In the case of unit cell test CFD analysis and HTTF CFD analysis, the volume-averaged sub-grid scale model dependency was calculated to be 13.0% and 9.16%, respectively. As a result of HTTF analysis, quantitative data of the fluid inside the HTTF lower plenum was provided in this paper. As a result of qualitative analysis, the temperature was highest at the center of the lower plenum, while the temperature fluctuation was highest near the edge of the lower plenum wall. The power spectral density of temperature was analyzed via fast Fourier transform (FFT) for specific points on the center and side of the lower plenum. FFT results did not reveal specific frequency-dominant temperature fluctuations in the center part. It was confirmed that the temperature power spectral density (PSD) at the top increased from the center to the wake. The vortex was visualized using the well-known scalar Q-criterion, and as a result, the closer to the outlet duct, the greater the influence of the mainstream, so that the inflow jet vortex was dissipated and mixed at the top of the lower plenum. Additionally, FFT analysis was performed on the support structure near the corner of the lower plenum with large temperature fluctuations, and as a result, it was confirmed that the temperature fluctuation of the flow did not have a significant effect near the corner wall. In addition, the vortices generated from the lower plenum to the outlet duct were identified in this paper. It is considered that the quantitative and qualitative results presented in this paper will serve as reference data for the benchmark.