• Advances in aircraft and spacecraft science
• /
• v.7 no.3
• /
• pp.215-228
• /
• 2020

As is shown in the paper, the Koltunov-Rzhanitsyn singular kernel of heredity (when constructing mathematical models of the dynamics problem of the hereditary theory of viscoelasticity) adequately describes real mechanical processes, best approximates experimental data for a long period of time. A mathematical model of the problem of the flutter of viscoelastic plates moving in a gas with a high supersonic velocity is given. Using the Bubnov-Galerkin method, discrete models of the problem of the flatter of viscoelastic plates flowed over by supersonic gas flow are obtained. A numerical method is developed to solve nonlinear integro-differential equations (IDE) for the problem of the hereditary theory of viscoelasticity with weakly singular kernels. A general computational algorithm and a system of application programs have been developed, which allow one to investigate the nonlinear dynamic problems of the hereditary theory of viscoelasticity with weakly singular kernels. On the basis of the proposed numerical method and algorithm, nonlinear problems of the flutter of viscoelastic plates flowed over in a gas flow at an arbitrary angle are investigated. In a wide range of changes in various parameters of the plate, the critical velocity of the flutter is determined. It is shown that the singularity parameter α affects not only the oscillations of viscoelastic systems, but the critical velocity of the flutter as well.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.13 no.5
• /
• pp.1032-1043
• /
• 1989

The heat transfer characteristics of the drag reducing polymer solutions are investigated experimentally in the thermal entrance region of circular tube flows. Fluids used in experiments are the aqueous solutions of high molecular polymer, polyacrylamide Separan AP-273 and the range of polymer concentrations is from 20 to 1000 wppm. Two stainless steel tubes with inside diameter 8.5mm(L/D=712) and 10.3mm(L/D=1160) are used for the heat transfer flow loops. The flow loop is set up to measure friction factors and heat transfer coefficients of test sections in two different modes; the recirculating flow system and once-through flow system. The test tubes are heated directly by electricity to apply the constant heat flux boundary conditions to the wall. Three different types of adaptors are used to observe the effects of the upstream flow conditions of the heat transfer test sections. The viscosity and characteristic relaxation time of the test fluids circulating in the flow system are measured by the capillary tube viscometer and falling ball viscometer at regular time intervals. The installed adaptors exhibit slight effect on the entrance heat transfer of Newtonian fluid. However, no noticeable effects are observed for the entrance heat transfer of the drag reducing fluids. The order of magnitude of the thermal entrance lengths of the drag reducing fluids which follow the minimum friction asymptote is much longer than that of Newtonian fluids in turbulent flows. A new dimensionless parameter, the viscoelastic Graetz number, is defined and all the experimental data are recasted in terms of the viscoelastic Graetz number. The local Nusselt number of the viscoelastic fluids is represented as a function of flow behavior index n and the viscoelastic Graetz number. As degradation continues the viscosity and the characteristic relaxation time of the testing fluids decrease. Weissenberg number defined by the relaxation time and D/V appears to be a proper dimensionless parameter in describing degradation effects on heat transfer of the viscoelastic fluids.

• Korea-Australia Rheology Journal
• /
• v.14 no.4
• /
• pp.161-174
• /
• 2002

A new technique for numerical calculation of viscoelastic flow based on the combination of Neural Net-works (NN) and Brownian Dynamics simulation or Stochastic Simulation Technique (SST) is presented in this paper. This method uses a "universal approximator" based on neural network methodology in combination with the kinetic theory of polymeric liquid in which the stress is computed from the molecular configuration rather than from closed form constitutive equations. Thus the new method obviates not only the need for a rheological constitutive equation to describe the fluid (as in the original Calculation Of Non-Newtonian Flows: Finite Elements St Stochastic Simulation Techniques (CONNFFESSIT) idea) but also any kind of finite element-type discretisation of the domain and its boundary for numerical solution of the governing PDE's. As an illustration of the method, the time development of the planar Couette flow is studied for two molecular kinetic models with finite extensibility, namely the Finitely Extensible Nonlinear Elastic (FENE) and FENE-Peterlin (FENE-P) models.P) models.

• Journal of Mechanical Science and Technology
• /
• v.20 no.3
• /
• pp.382-390
• /
• 2006

The present study proposes a modified temperature-dependent non-Newtonian viscosity model and investigates the flow characteristics and heat transfer enhancement of the viscoelastic non-Newtonian fluid in a 2:1 rectangular duct. The combined effects of temperature dependent viscosity, buoyancy, and secondary flow caused by the second normal stress difference are considered. Calculated Nusselt numbers by the modified temperature-dependent viscosity model give good agreement with the experimental results. The heat transfer enhancement of viscoelastic fluid in a rectangular duct is highly dependent on the secondary flow caused by the magnitude of second normal stress difference.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.05a
• /
• pp.192-198
• /
• 1999

The present study proposes modified temperature-dependent non-Newtonian viscosity model and investigates flow characters and heat transfer enhancement of the viscoelastic non-Newtonian fluid in a 2:1 rectangular duct. The proposed modified temperature dependent viscosity model has non-zero value near the high temperature and high shear rate region while on the existing viscosity models have zero value. Two versions of thermal boundary conditions involving difference combination of heated walls and adiabatic walls are analyzed in this study. The combined effect of temperature dependent viscosity, buoyancy, and secondary flow caused by second normal stress difference are ail considered. The Reiner-Rivlin model is adopted as a viscoelastic fluid model to simulate the secondary flow caused by second normal stress difference. Calculated Nusselt numbers by the modified temperature-dependent viscosity model gives under prediction than the existing temperature-dependent viscosity model in the regions of thermally developed with same secondary normal stress difference coefficients with experimental results in the regions of thermally developed. The heat transfer enhancement of the viscoelastic fluid in a 2:1 rectangular duct is highly dependent on the secondary flow caused by the magnitude of second normal stress difference.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.22 no.9
• /
• pp.1208-1216
• /
• 1998

The present numerical study investigates the effect of a secondary flow on the heat transfer in order to delineate the mechanism of laminar heat transfer enhancement of a viscoelastic fluid in rectangular ducts. The second normal stress generating a secondary flow is modeled by adopting the Reiner-Rivlin constitutive equation and the calculated secondary flow showed good agreement with experiments. The primary velocity U as well as the pressure drop were not affected by the secondary flow in rectangular ducts, whose order of magnitude is less than 0.1% of the primary velocity. The small magnitude of the secondary flow, however, affect moderately the temperature fields. The calculated Nusselt numbers with secondary flow show 50% heat transfer enhancement over those of a purely viscous non-Newtonian fluid, which are considerably lower than the experimental values. Therefore, we conclude that there should be an additional heat transfer enhancement mechanism involved in the viscoelastic fluid such as temperature-dependence.

• Korea-Australia Rheology Journal
• /
• v.18 no.2
• /
• pp.67-81
• /
• 2006

Using a strain-controlled rheometer, the dynamic viscoelastic properties of aqueous xanthan gum solutions with different concentrations were measured over a wide range of strain amplitudes and then the linear viscoelastic behavior in small amplitude oscillatory shear flow fields was investigated over a broad range of angular frequencies. In this article, both the strain amplitude and concentration dependencies of dynamic viscoelastic behavior were reported at full length from the experimental data obtained from strain-sweep tests. In addition, the linear viscoelastic behavior was explained in detail and the effects of angular frequency and concentration on this behavior were discussed using the well-known power-law type equations. Finally, a fractional derivative model originally developed by Ma and Barbosa-Canovas (1996) was employed to make a quantitative description of a linear viscoelastic behavior and then the applicability of this model was examined with a brief comment on its limitations. Main findings obtained from this study can be summarized as follows: (1) At strain amplitude range larger than 10%, the storage modulus shows a nonlinear strain-thinning behavior, indicating a decrease in storage modulus as an increase in strain amplitude. (2) At strain amplitude range larger than 80%, the loss modulus exhibits an exceptional nonlinear strain-overshoot behavior, indicating that the loss modulus is first increased up to a certain strain amplitude(${\gamma}_0{\approx}150%$) beyond which followed by a decrease in loss modulus with an increase in strain amplitude. (3) At sufficiently large strain amplitude range (${\gamma}_0>200%$), a viscous behavior becomes superior to an elastic behavior. (4) An ability to flow without fracture at large strain amplitudes is one of the most important differences between typical strong gel systems and concentrated xanthan gum solutions. (5) The linear viscoelastic behavior of concentrated xanthan gum solutions is dominated by an elastic nature rather than a viscous nature and a gel-like structure is present in these systems. (6) As the polymer concentration is increased, xanthan gum solutions become more elastic and can be characterized by a slower relaxation mechanism. (7) Concentrated xanthan gum solutions do not form a chemically cross-linked stable (strong) gel but exhibit a weak gel-like behavior. (8) A fractional derivative model may be an attractive means for predicting a linear viscoelastic behavior of concentrated xanthan gum solutions but classified as a semi-empirical relationship because there exists no real physical meaning for the model parameters.

• Korea-Australia Rheology Journal
• /
• v.19 no.4
• /
• pp.191-199
• /
• 2007

Precision injection molding process is of great importance since precision optical products such as CD, DVD and various lens are manufactured by those process. In such products, birefringence affects the optical performance while residual stress that determines the geometric precision level. Therefore, it is needed to study residual stress and birefringence that affect deformation and optical quality, respectively in precision optical product. In the present study, we tried to predict residual stress, final shrinkage and birefringence in injection molded parts in a systematic way, and compared numerical results with the corresponding experimental data. Residual stress and birefringence can be divided into two parts, namely flow induced and thermally induced portions. Flow induced birefringence is dominant during the flow, whereas thermally induced stress is much higher than flow induced one when amorphous polymer undergoes rapid cooling across the glass transition region. A numerical system that is able to predict birefringence, residual stress and final shrinkage in injection molding process has been developed using hybrid finite element-difference method for a general three dimensional thin part geometry. The present modeling attempts to integrate the analysis of the entire process consistently by assuming polymeric materials as nonlinear viscoelastic fluids above a no-flow temperature and as linear viscoelastic solids below the no-flow temperature, while calculating residual stress, shrinkage and birefringence accordingly. Thus, for flow induced ones, the Leonov model and stress-optical law are adopted, while the linear viscoelastic model, photoviscoelastic model and free volume theory taking into account the density relaxation phenomena are employed to predict thermally induced ones. Special cares are taken of the modeling of the lateral boundary condition which can consider product geometry, histories of pressure and residual stress. Deformations at and after ejection have been considered using thin shell viscoelastic finite element method. There were good correspondences between numerical results and experimental data if final shrinkage, residual stress and birefringence were compared.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2006.05a
• /
• pp.99-104
• /
• 2006

A semi-Lagrangian finite element scheme with objective time stepping algorithm for solving viscoelastic flow problem is presented. The convection terms in the momentum and constitutive equations are treated using a quasi-monotone semi-Lagrangian scheme, in which characteristic feet on a regular grid are traced backwards over a single time-step. Concerned with the generalized midpoint rule type of algorithms formulated to exactly preserve objectivity, we use the geometric transformation such as pull-back, push-forward operation. The method is applied to the 4:1 planar contraction problem for an Oldroyd B fluid for both creeping and inertial flow conditions.

• Journal of the Korean Mathematical Society
• /
• v.60 no.2
• /
• pp.273-302
• /
• 2023

In this paper, a fully discrete numerical scheme for the viscoelastic Oldroyd flow is considered with an introduced auxiliary variable. Our scheme is based on the finite element approximation for the spatial discretization and the backward Euler scheme for the time discretization. The integral term is discretized by the right trapezoidal rule. Firstly, we present the corresponding equivalent form of the considered model, and show the relationship between the origin problem and its equivalent system in finite element discretization. Secondly, unconditional stability and optimal error estimates of fully discrete numerical solutions in various norms are established. Finally, some numerical results are provided to confirm the established theoretical analysis and show the performances of the considered numerical scheme.