• Korea-Australia Rheology Journal
• /
• v.21 no.1
• /
• pp.47-58
• /
• 2009

We have investigated the transient behavior of 1D fully developed Poiseuille viscoelastic flow under finite pressure gradient described by the Oldroyd-B and Leonov constitutive equations. For analysis we employ a simple $2^{nd}$ order discretization scheme such as central difference for space and the Crank-Nicolson for time approximation. For the analysis of the Oldroyd-B model, we also apply the analytical solution, which is obtained again in this work in terms of elementary solution procedure simpler than the previous one (Waters and King, 1970). Both models demonstrate qualitatively similar solutions, but their eventual steady flowrate exhibits noticeable difference due to the absence or presence of shear thinning behavior. In the inertialess flow, the flowrate instantaneously attains a large value corresponding to the Newtonian creeping flow and then decreases to its steady value when the applied pressure gradient is low. However with finite liquid density the flow field shows severe fluctuation even accompanying reversals of flow directions. As the assigned pressure gradient increases, the flowrate achieves its steady value significantly higher than its value during oscillations after quite long period of time. We have also illustrated comparison between 1D and 2D results and possible mechanism of complex 2D flow rearrangement employing a previous solution of [mite element computation. In addition, we discuss some mathematical points regarding missing boundary conditions in 2D modeling due to the change of the type of differential equations when varying from inertialess to inertial flow.

• Wind and Structures
• /
• v.22 no.3
• /
• pp.349-368
• /
• 2016

A preconditioning technique is presented for a simultaneous solution to wind-membrane interaction. In the simultaneous equations, a linear elastic model was employed to deal with the fluid-structure data transfer at the interface. A Lagrange multiplier was introduced to impose the specified boundary conditions at the interface and strongly coupled simultaneous equations are derived after space and time discretization. An initial linear elastic model preconditioner and modified one were derived by treating the linearized elastic model equation as a saddle point problem, respectively. Accordingly, initial and modified fluid-structure interaction (FSI) preconditioner for the simultaneous equations were derived based on the initial and modified linear elastic model preconditioners, respectively. Wind-membrane interaction analysis by the proposed preconditioners, for two and three dimensional membranous structures respectively, was performed. Comparison was made between the performance of initial and modified preconditioners by comparing parameters such as iteration numbers, relative residuals and convergence in FSI computation. The results show that the proposed preconditioning technique greatly improves calculation accuracy and efficiency. The priority of the modified FSI preconditioner is verified. The proposed preconditioning technique provides an efficient solution procedure and paves the way for practical application of simultaneous solution for wind-structure interaction computation.

• Journal of KSNVE
• /
• v.8 no.4
• /
• pp.632-642
• /
• 1998

A method based on finite element discretization is developed for optimizing the polarization profile of PVDF film to create the modal transducer for specific modes. Using this concept, one can design the modal transducer in two-dimensional structure having arbitrary geometry and boundary conditions. As a practical means for implementing this polarization profile without repoling the PVDF film the polarization profile is approximated by optimizing electrode patterns, lamination angles, and poling directions of the multi-layered PVDF transducer. This corresponds to the approximation of a continuous function using discrete values. The electrode pattern of each PVDF layer is optimized by deciding the electrode of each finite element to be used or not. Genetic algorithm, suitable for discrete problems, is used as an optimization scheme. For the optimization of each layers lamination angle, the continuous lamination angle is encoded into discrete value using binary 5 bit string. For the experimental demonstration, a modal sensor for first and second modes of cantilevered composite plate is designed using two layers of PVDF films. The actuator is designed based on the criterion of minimizing the system energy in the control modes under a given initial condition. Experimental results show that the signals from residual modes are successfully reduced using the optimized multi-layered PVDF sensor. Using discrete LQG control law, the modal peaks of first and second modes are reduced in the amount of 12 dB and 4 dB, resepctively.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1992.10a
• /
• pp.117-123
• /
• 1992

In recent years, the finite element method has become one of the most popular numerical technique for obtaining solutions of engineering science problems. However, there exist various uncertainties in modeling the problems, such as the dimensions(geometry shape), the material properties, boundary conditions, etc. The consideration for the uncertainties inherent in the problems can be made by understanding the influences of uncertain parameters[1]. Determining the influences of uncertainties as statistical quantities using the standard finite element method requires enormous computing time, while the probabilistic finite element method is realized as an efficient scheme[2,3] yielding statistical solution with just a few direct computations. In this paper, a formulation of the probabilistic fluid-structure interaction problem accounting for the first order perturbation of geometric shape is derived, and especially probabilistical acoustic pressure scattering from the structure with surrounding fluid is focused on. In Section 2, governing equations for the fluid-structure problems are given. In Section 3, a finite element formulation, based on the functional, is presented. First order perturbation of geometric shape with randomness is incorporated into the finite element formulation in conjunction with discretization of the random fields in Section 4 and 5. Finally, the proposed formulation is applied to a acoustic pressure scattering problem from an infinitely long cylindrical shell structure with randomness of radial perturbation.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.10
• /
• pp.1384-1391
• /
• 2001

Three-dimensional flow analysis is implemented to investigate the flow through transonic axial-flow compressor rotor(NASA R67) and to evaluate the performances of Abid's low-Reynolds-number k-$\omega$ and Baldwin-Lomax turbulence models. A finite volume method is used fur spatial discretization. The equations are solved implicitly in time by the use of approximate factorization. The upwind difference scheme is used for inviscid terms and viscous terms are approximated with central difference. The flux-difference-splitting method of Roe is used to obtain fluxes at the cell faces. Numerical analysis is performed near peak efficiency and near stall. The results are compared with the experimental data for NASA R67 rotor. Blade-to-Blade Mach number distributions are compared to confirm the accuracy of the code. From the results, it is concluded that Abid'k-$\omega$ model is better for the calculation of flow rate and efficiency than Baldwin-Lomax model. But, the predictions for Mach number and shock structure are almost the same.

• Proceedings of the KIEE Conference
• /
• 1994.07a
• /
• pp.208-210
• /
• 1994

In this paper, the design sensitivity evaluation based on the 2-dimensional finite clement discretization is presented for the voltage source and eddy current problem. And it is applied to the two shape design problems of the rotor slot of 3-phase squirrel cage induction motors. The first is to increase the starting torque while keeping the rated torque fixed. The other is only to increase the torque at the rated speed while keeping the starting torque fixed. As an optimization method, the Gradient Projection method is used to control casily the torques for various speeds of rotor. One fourth of rotor is analyzed by using a semi-periodic boundary condition. Because the shape of rotor slot has much influence on the slip torque characteristic, the 10 design parameters are taken on the interface between rotor core and rotor bar. The initial shape of rotor slot is the trapezoidal typo with rounding corners.

• Nuclear Engineering and Technology
• /
• v.39 no.5
• /
• pp.627-636
• /
• 2007

The present work is aimed at further discussing the effectiveness of dimensionless parameters recently proposed for the analysis of flow stability in heated channels with supercritical fluids. In this purpose, after presenting the main motivations for the introduction of these parameters in place of previously proposed ones, additional information on the theoretical bases and on the consequences of this development is provided. Stability maps, generated by an in-house program adapted from a previous application to boiling channels, are also shown for different combinations of the operating parameters. The maps are obtained as contour plots of an amplification parameter obtained from numerical discretization and subsequent linearization of governing equations; as such, they provide a quantitatively clear perspective of the effect of different boundary conditions on the stability of heated channels with supercritical fluids. In order to assess the validity of the assumptions at the basis of the in-house model, supporting calculations have been performed making use of the RELAP5/MOD3.3 computer code, detecting the values of the dimensionless parameters at the threshold for the occurrence of instability for a heated channel representative of SCWR proposed core configurations. The obtained results show reasonable agreement with the maps, supporting the applicability of the proposed scaling parameters for describing the dynamic behaviour of heated channels with supercritical fluids.

• Wind and Structures
• /
• v.17 no.1
• /
• pp.1-19
• /
• 2013

A multiscale finite element method is applied to the Spalart-Allmaras turbulence model based detached-eddy simulation (DES). The multiscale arises from a decomposition of the scalar field into coarse (resolved) and fine (unresolved) scales. It corrects the lack of stability of the standard Galerkin formulation by modeling the scales that cannot be resolved by a given spatial discretization. The stabilization terms appear naturally and the resulting formulation provides effective stabilization in turbulent computations, where reaction-dominated effects strongly influence near-wall predictions. The multiscale DES is applied in the context of high-Reynolds flow over the Commonwealth Advisory Aeronautical Council (CAARC) standard tall building model, for both uniform and turbulent inflows. Time-averaged pressure coefficients on the exterior walls are compared with experiments and it is demonstrated that DES is able to resolve the turbulent features of the flow and accurately predict the surface pressure distributions under atmospheric boundary layer flows.

• Structural Engineering and Mechanics
• /
• v.40 no.2
• /
• pp.257-277
• /
• 2011

In this paper the geometrically nonlinear continuum plate finite element model, hitherto not reported in the literature, is developed using the total Lagrange formulation. With the layerwise displacement field of Reddy, nonlinear Green-Lagrange small strain large displacements relations (in the von Karman sense) and linear elastic orthotropic material properties for each lamina, the 3D elasticity equations are reduced to 2D problem and the nonlinear equilibrium integral form is obtained. By performing the linearization on nonlinear integral form and then the discretization on linearized integral form, tangent stiffness matrix is obtained with less manipulation and in more consistent form, compared to the one obtained using laminated element approach. Symmetric tangent stiffness matrixes, together with internal force vector are then utilized in Newton Raphson's method for the numerical solution of nonlinear incremental finite element equilibrium equations. Despite of its complex layer dependent numerical nature, the present model has no shear locking problems, compared to ESL (Equivalent Single Layer) models, or aspect ratio problems, as the 3D finite element may have when analyzing thin plate behavior. The originally coded MATLAB computer program for the finite element solution is used to verify the accuracy of the numerical model, by calculating nonlinear response of plates with different mechanical properties, which are isotropic, orthotropic and anisotropic (cross ply and angle ply), different plate thickness, different boundary conditions and different load direction (unloading/loading). The obtained results are compared with available results from the literature and the linear solutions from the author's previous papers.

• Journal of the Society of Naval Architects of Korea
• /
• v.31 no.3
• /
• pp.87-99
• /
• 1994

Reynolds Averaged Navier-Stokes equations are solved numerically for the computation of turbulent flow around a Wigley double model. A second order finite difference method is applied for the spatial discretization on the nonstaggered grid system and 4-stage Runge-Kutta scheme for the numerical integration in time. In order to increase the time step, residual averaging scheme of Jameson is adopted. Pressure field is obtained by solving the pressure-Poisson equation with the appropriate Neumann boundary condition. For the turbulence closure, 0-equation turbulence model of Baldwin-Lomax is used. Numerical computation is carried out for the Reynolds number of 4.5 million. Comparisons of the computed results with the available experimental data show good agreements for the velocity and pressure distributions.