• Journal of computational fluids engineering
• /
• v.19 no.2
• /
• pp.58-65
• /
• 2014

The adaptive wavelet method is studied for the enhancement of computational efficiency of three-dimensional flows. For implementation of the method for three-dimensional Euler equation, wavelet decomposition process is introduced based on the previous two-dimensional adaptive wavelet method. The order of numerical accuracy of an original solver is preserved by applying modified thresholding value. In order to assess the efficiency of the proposed algorithm, the method is applied to the computation of flow field around ONERA-M6 wing in transonic regime with 4th and 6th order interpolating polynomial respectively. Through the application, it is confirmed that the three-dimensional adaptive wavelet method can reduce the computational time while conserving the numerical accuracy of an original solver.

• Journal of computational fluids engineering
• /
• v.2 no.1
• /
• pp.8-12
• /
• 1997

The three-dimensional Euler equations are solved numerically for the analysis of contraction flows in wind or water tunnels. A second-order finite difference method is used for the spatial discretization on the nonstaggered grid system and the 4-stage Runge-Kutta scheme for the numerical integration in time. In order to speed up the convergence, the local time stepping and the implicit residual-averaging schemes are introduced. The pressure field is obtained by solving the pressure-Poisson equation with the Neumann boundary condition. For the evaluation of the present Euler solver, numerical computations are carried out for three contraction geometries, one of which was adopted in the Large Cavitation Channel for the U.S. Navy. The comparison of the computational results with the available experimental data shows good agreement.

• 한국전산유체공학회:학술대회논문집
• /
• 1995.10a
• /
• pp.175-181
• /
• 1995

Three-Dimensional Euler equations are solved numerically for the analysis of contraction flows in wind or water tunnels. A second-order finite difference method is used for the spatial discretization on the nonstaggered grid system and the 4-stage Runge-Kutta scheme for the numerical integration in time. In order to speed up the convergence, the local time stepping and the implicit residual-averaging schemes are introduced. The pressure field is obtained by solving the pressure-Poisson equation with the Neumann boundary condition. For the evaluation of the present Euler solver, numerical computations are carried out for the various contraction geometries, one of which was adopted in the Large Cavitation Channel for the U.S. Navy. The comparison of the computational results with the available experimental data shows good agreements.

• Journal of computational fluids engineering
• /
• v.17 no.4
• /
• pp.103-109
• /
• 2012

Sound scattering in 3D complex geometry is difficult to model with body-fitted grid. Thus Brinkman Penalization method is used to compute sound scattering in 3D complex geometry. Sound propagation of monitor/TV is studied. The sound field for monitor/TV is simulated by applying Brinkman Penalization method to Linearized Euler Equation. Solid Structure and ambient air are represented as penalty terms in Linearized Euler Equation.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1997.04a
• /
• pp.627-635
• /
• 1997

In general, the booming noise intensity at tunnel exit is strongly related to the gradient of the compression wave front created by high speed train entering the tunnel. This paper presents some results in relation with the compression wave front produced when the high speed train enters a tunnel. Four kinds of tunnel entrance shape with real dimensions were studied to investigate the formation of compression wave front inside tunnel by train entering tunnel. Computations were carried out using three-dimensional compressible Euler equation with vanishing viscosity and conductivity of fluid. According to the reslts, the flow disturbance occured at tunnel entrance were eliminated by tunnel hood with same cross sectional area. The compression wave front is formed completely at 30-40m from tunnel entrance. The maximum pressure gradient of compression wave front is reduced by 29.8% for the inclined tunnel hood and reduced by 21.5% for the tunnel hood with holes at the top face with tunnel without hood. The length of the inclined hood is 15m and the length of the hood with holes is 20m.

• Journal of KSNVE
• /
• v.7 no.6
• /
• pp.959-966
• /
• 1997

In general, the booming noise intensity at tunnel exit is strongly related to the gradient of the compression wave front created by high speed train entering the tunnel. This paper presents some results in relation with the compression wave front produced when the high speed train enters a tunnel. Four kinds of tunnel entrance shape with real dimensions were studied to investigate the formation of compression wave front inside tunnel by train entering tunnel. Computations were carried out using three-dimensional compressible Euler equation with vanishing viscosity and conductivity of fluid. According to the results, the flow disturbances occured at tunnel entrance were eliminated by tunnel hood with same cross sectional area. The compression wave front is formed completely at 30-40m from tunnel entrance. The maximum pressure gradient of compression wave front is reduced by 29.8% for the inclined tunnel hood and reduced by 21.5% for the tunnel hood with holes at the top face with tunnel without hood. The length of the inclined hood is 15m and the length of the hood with holes is 20m.

• Journal of applied mathematics & informatics
• /
• v.42 no.1
• /
• pp.85-92
• /
• 2024

In this work we give the necessary and sufficient conditions for f-biharmonic curves in three-dimensional generalized symmetric spaces.

• Journal of Advanced Marine Engineering and Technology
• /
• v.39 no.10
• /
• pp.1023-1030
• /
• 2015

The purpose of this paper is to present the basic mathematical modeling of a hexacopter, which could be used to develop proper methods for stabilization and trajectory control. A hexacopter consists of six rotors with three pairs of counter-rotating fixed-pitch blades. This mechanism is an under-actuated, dynamically unstable, six-degrees-of-freedom system. The whole motion of this object consists of translational and rotational motion in three dimensions, where the translational motion is created by changing the direction and magnitude of the upward propeller thrust. The hexacopter is controlled by adjusting the angular velocities of the rotors, which are spun by electric motors. It is assumed to be a rigid body; thus, the differential equation of the hexacopter dynamics can be derived from the Newton-Euler equation. The Euler-angle parametrization of the three-dimensional rotations contains singular points in the coordinate space that can cause failure of both the dynamical model and control. In order to avoid singularities, the rotations of the hexacopter are parametrized in terms of quaternions. This choice has been made considering the linearity of the quaternion formulation and their stability and efficiency. Further, control simulation of a hexacopter applying cascaded-PID control is also presented in this paper.

• Journal of the Chungcheong Mathematical Society
• /
• v.31 no.4
• /
• pp.363-368
• /
• 2018

The aim of this work is to derive a partial differential equation that explains the movement of vortex lines on a vortex trajectory surface in a three dimensional incompressible inviscid flow.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.10 no.1
• /
• pp.102-109
• /
• 1986

Based on the Hartenberg-Denavit symbolic equation, which is one of equations for the kinematic analysis of three dimensional (3-D) linkage, a generalized kinematic motion equation is derived utilizing Euler angles and employing the coordinates transformation. The derived equation can feasibly be used for the motion analysis of any type of 3-D linkages as well as 2-D ones. In order to simulate the general motion of 3-D linkgages on digital computer, the generalized equation is programmed through the process of numerical analysis after converting the equation to the type of Newton-Raphson formula and denoting it in matrix form. The feasibility of theoretically derived equation is experimentally proved by comparing the results from the computer with those from experimental setup of three differrent but generally empolyed 3-D linkages.