• 한국자동차공학회논문집
• /
• 제6권2호
• /
• pp.12-20
• /
• 1998

Exter ballistics of a typical high-speed projectile is studied through a flow-visualization experiment and an unstructured grid Navier-Srokes computation. Experiment produced a schlieren photograph that adequately shows the characteristic features of this complex flow, namely two kinds of oblique cone shocks and turbulent wake developing into the downstream. A hybrid scheme of finite volume-element method is used to simulate the compressible Reynolds-Averaged Navier-Stok- es solution on unstructured grids. Osher's approximate Riemann solver is used to discretize the cinvection term. Higher-order spatial accuracy is obtained by MUSCL extension and van Albada ty- pe flux limiter is used to stabilize the numerical oscillation near the solution discontinuity. Accurate Gakerkin method is used to discretize the viscous term. Explict fourth-order Runge-Kutta method is used for the time-stepping, which simplifies the application of MUSCL extension. A two-layer k-$\varepsilon$ turbulence model is used to simulate the turbulent wakes accurately. Axisymmetric folw and two-dimensional flow with an angle of attack have been computed. Grid-dependency is also checked by carrying out the computation with doubled meshes. 2-D calculation shows that effect of angle of attack on the flow field is negligible. Axi-symmetric results of the computation agrees well with the flow visualization. Primary oblique shock is represented within 2-3 meshes in numerical results, and the varicose mode of the vortex shedding is clearly captured in the turbulent wake region.

• Structural Engineering and Mechanics
• /
• 제40권3호
• /
• pp.431-451
• /
• 2011

Time delay inevitably exists in active control systems, and it may cause the degradation of control efficiency or instability of the systems. So time delay needs to be compensated in control design in order to eliminate its negative effect on control efficiency. Today time delay in linear systems has been more studied and some treating methods had been worked out. However, there are few treating methods for time delay in nonlinear systems. In this paper, an active controller for a nonlinear and hysteretic building structure with time delay is studied. The nonlinear and hysteretic behavior of the system is illustrated by the Bouc-Wen model. By specific transformation and augmentation of state parameters, the motion equation of the system with explicit time delay is transformed into the standard state space representation without any explicit time delay. Then the fourth-order Runge-Kutta method and instantaneous optimal control method are applied to the controller design with time delay. Finally, numerical simulations and comparisons of an eight-story building using the proposed time-delay controller are carried out. Simulation results indicate that the control performance will deteriorate if time delay is not taken into account in the control design. The simulations also prove the proposed time delay controller in this paper can not only effectively compensate time delay to get better control effectiveness, but also work well with both small and large time delay problems.

• Advances in aircraft and spacecraft science
• /
• 제8권4호
• /
• pp.345-372
• /
• 2021

The analysis of nonlinear vibrations, buckling, post-buckling, flutter boundary determination and post-flutter behavior of a homogeneous curved plate assuming cylindrical bending is conducted in this article. Other assumptions include simply-supported boundary conditions, supersonic aerodynamic flow at the top of the plate, constant pressure conditions below the plate, non-viscous flow model (using first- and third-order piston theory), nonlinear structural model with large deformations, and application of mechanical and thermal loads on the curved plate. The analysis is performed with constant environmental indicators (flow density, heat, Reynolds number and Mach number). The material properties (i.e., coefficient of thermal expansion and modulus of elasticity) are temperature-dependent. The equations are derived using the principle of virtual displacement. Furthermore, based on the definitions of virtual work, the potential and kinetic energy of the final relations in the integral form, and the governing nonlinear differential equations are obtained after fractional integration. This problem is solved using two approaches. The frequency analysis and flutter are studied in the first approach by transferring the handle of ordinary differential equations to the state space, calculating the system Jacobin matrix and analyzing the eigenvalue to determine the instability conditions. The second approach discusses the nonlinear frequency analysis and nonlinear flutter using the semi-analytical solution of governing differential equations based on the weighted residual method. The partial differential equations are converted to ordinary differential equations, after which they are solved based on the Runge-Kutta fourth- and fifth-order methods. The comparison between the results of frequency and flutter analysis of curved plate is linearly and nonlinearly performed for the first time. The results show that the plate curvature has a profound impact on the instability boundary of the plate under supersonic aerodynamic loading. The flutter boundary decreases with growing thermal load and increases with growing curvature.

• 대한의용생체공학회:의공학회지
• /
• 제27권4호
• /
• pp.143-153
• /
• 2006

The purpose of this work was to assess and quantify the beneficial effects of gas exchange, while testingto the various frequencies of the sinusoidal wave that was excited by the PZT actuator, for patients suffering from acute respiratory distress syndrome (ARDS) or chronic respiratory problems. Also, this paper considered a simulator to design a hollow type artificial lung, and a mathematical model was used to predict a behavior of blood. This simulation was carried out according to the Montecarno's simulation method, anda fourth order Runge-Kutta method was used to solve the equation. The experimental design and procedure are then applied to the construction of a new device to assess the effectiveness of the membrane vibrations. As a result, the vibration method is very effective in the increase of gas transport. The gas exchange efficiency for the vibrating intravascular lung assist device can be increased by emphasizing the following design features: consistent and reproducible fiber geometry, and most importantly, an active means of enhancing convective mixing of water around the hollow fiber membranes. The experimental results showed the effective performance of the vibrating intravascular lung assist device. Also, we concluded that important design parameters were blood flow rates, fiber outer diameter and oxygen pressure drop. Based on the present results, it was believed that the optimal level of blood flow rates was 200$cm^3$/min.

• 한국해양공학회지
• /
• 제33권5호
• /
• pp.400-410
• /
• 2019

In this study, the numerical code for the 3D nonlinear dynamic analysis of an SLWR (Steel Lazy Wave Riser) was developed using the lumped mass line model in a FORTRAN environment. Because the lumped mass line model is an explicit method, there is no matrix operation. Thus, the numerical algorithm is simple and fast. In the lumped mass line model, the equations of motion for the riser were derived by applying the various forces acting on each node of the line. The applied forces at the node of the riser consisted of the tension, shear force due to the bending moment, gravitational force, buoyancy force, riser/ground contact force, and hydrodynamic force based on the Morison equation. Time integration was carried out using a Runge-Kutta fourth-order method, which is known to be stable and accurate. To validate the accuracy of the developed numerical code, simulations using the commercial software OrcaFlex were carried out simultaneously and compared with the results of the developed numerical code. To understand the nonlinear dynamic characteristics of an SLWR, dynamic simulations of SLWRs excited at the hang-off point and of SLWRs in regular waves were carried out. From the results of these dynamic simulations, the displacements at the maximum bending moments at important points of the design, like the hang-off point, sagging point, hogging points, and touch-down point, were observed and analyzed.

• 한국생물공학회:학술대회논문집
• /
• 한국생물공학회 2005년도 생물공학의 동향(XVI)
• /
• pp.65-66
• /
• 2005

Anaerobic digestion is one of the well-known methods for biological treatment handling of concentrated organic matter such as swine $wastewater.^{1)} The anaerobic digestion can reduce organic loading but also hydrolyze non-biodegradable organic $matter.^{2)}$ The feces from the scrapper-type barn are usually collected to make compost and the urine is discarded with swine-slurry wastewater by ocean-dumping or treated by biological methods. The lagoon, aerobic digestion, anaerobic digestion, SBR, $A^{2}/O$, and UCT have been applied for treating swine $wastewater.^{3)} In this study, as a result of the analysis of swine wastewater, the total and soluble chemical oxygen demand was 130g/L and 60g/L, respectively. And the volatile fatty acid as chemical oxygen demand equivalent was 45g/L, which was 75% of soluble chemical oxygen demand. Before everything else, ammonia nitrogen concentration was 6.5 g/L. From biochemical acidogenic potential test, it was concluded that the enhanced acidification process to manage swine waste should be operated in the ammonia nitrogen concentration of less than 1.2 g/L. In the result of seeding ratio experiments with artificial $wastewater^{4)}, the lag period of acidogens was taken the long time because of the inhibition by the $ammonia^{5)}$, however no difference of period by the seeding ratio was not shown. The Haldane-based biokinetics were also evaluated using a method of fourth order Runge-Kutta $approximation.^{6,7)}$ The nonlinear least squares (NLLS) method with a 95% confidence interval was also used. The ranges of maximum microbial growth rate, ${/mu_{max}}$, and half saturation coefficient, $K_{s}$, for acidogenesis of various seeding ratio with artificial wastewater were 6.1 ~ 12.6 $d^{-1}$ and 45,000 ~ 53,500 mg glucose/L, respectively. Also, the methanogenic microbial yield coefficient, Y, and microbial decay rate coefficient, $k_{d}$, and inhibition substrate concentration, $K_{si}$, for the reactors were determined to be 0.32 ~ 0.465 ${/mu}g$/mg glucose; 0.42 ~ 1.01 $d^{-1}$ and 51,500 ~ 55,600 mg glucose/L, respectively.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제10권6호
• /
• pp.692-710
• /
• 2018

The Three-dimensional (3-D) dynamical behaviors of a fluid-conveying pipe subjected to vortex-induced vibration are investigated with different internal flow velocity ${\nu}$. The values of the internal flow velocity are considered in both subcritical and supercritical regimes. During the study, the 3-D nonlinear equations are discretized by the Galerkin method and solved by a fourth-order Runge-Kutta method. The results indicate that for a constant internal flow velocity ${\nu}$ in the subcritical regime, the peak Cross-flow (CF) amplitude increases firstly and then decrease accompanied by amplitude jumps with the increase of the external reduced velocity. While two response bands are observed in the In-line (IL) direction. For the dynamics in the lock-in condition, 3-D periodic, quasi-periodic and chaotic vibrations are observed. A variety of CF and IL responses can be detected for different modes with the increase of ${\nu}$. For the cases studied in the supercritical regime, the dynamics shows a great diversity with that in the subcritical regime. Various dynamical responses, which include 3-D periodic, quasi-periodic as well as chaotic motions, are found while both CF and IL responses are coupled while ${\nu}$ is beyond the critical value. Besides, the responses corresponding to different couples of ${\mu}_1$ and ${\mu}_2$ are obviously distinct from each other.

• 한국해양공학회지
• /
• 제32권6호
• /
• pp.447-457
• /
• 2018

In this study, a two-dimensional fully nonlinear transient wave numerical tank was developed using a desingularized indirect boundary integral equation method. The desingularized indirect boundary integral equation method is simpler and faster than the conventional boundary element method because special treatment is not required to compute the boundary integral. Numerical simulations were carried out in the time domain using the fourth order Runge-Kutta method. A mixed Eulerian-Lagrangian approach was adapted to reconstruct the free surface at each time step. A numerical damping zone was used to minimize the reflective wave in the downstream region. The interpolating method of a Gaussian radial basis function-type artificial neural network was used to calculate the gradient of the free surface elevation without element connectivity. The desingularized indirect boundary integral equation using an isolated point source and radial basis function has no need for information about the element connectivity and is a meshless method that is numerically more flexible. In order to validate the accuracy of the numerical wave tank based on the desingularized indirect boundary integral equation method and meshless technique, several numerical simulations were carried out. First, a comparison with numerical results according to the type of desingularized source was carried out and confirmed that continuous line sources can be replaced by simply isolated sources. In addition, a propagation simulation of a $2^{nd}$-order Stokes wave was carried out and compared with an analytical solution. Finally, simulations of propagating waves in shallow water and propagating waves over a submerged bar were also carried and compared with published data.