Temperature Oscillations into a Couette-Poiseuille Flow

Following previous work that discussed temperature fluctuations without flowing media a physical model of temperature oscillations into a Couette-Poiseuille flow was built. The temperature distribution into the flow was calculated according to oscillations constraints on the upper and lower plates, and heat dissipation due to shear stresses into the fluid. The physical model deals with different temperature amplitudes and different frequencies constraints on the upper and the lower plates. A physical superposition and complex numbers were used. It was shown that when the constraint frequency increases, its penetration capacity is reduced. Increasing gap width between plates leads to increased fluid temperature values due to enlarged fluid velocity. Increasing thermal diffusivity, increases constrains temperatures penetration intensity.


Introduction
The study of fluid flow in a Couette flow encompasses a few subjects as magneto hydrodynamics, porous media, particles suspension or two phases flow, and heat transfer in micro and nanoscale.These subjects applications relates for example to heat exchangers, heat pipes, electronic cooling, geophysics, biomedical instrumentation and gas flow at microscales.Several works relevant to these topics are cited below.
By using the Network Simulation Method (NSM), Beg et al. [1] examined the unsteady Hartmann-Coeutte flow and heat transfer in a parallel plate channel system containing a Darcian porous medium with lateral wall mass flux.The fluid was assumed as laminar, viscous and incompressible.Pressure gradient S. Sadik DOI: 10.4236/epe.2018.109026415 Energy and Power Engineering along the plates was assumed to be constant.Beg et al. work deals with some effects as Hall current, ion slip, viscous dissipation and Joule heating under a strong uniform transverse magnetic field.Constants values were determined for the temperatures boundary conditions on the upper and the lower plates.It was revealed that fluid temperature increased with increasing Darcy number (Da).Some conclusions: Increasing lateral wall mass flux-suction at the upper plate and injection at the lower plate decreased fluid principles velocities and temperature values.Increasing Hartmann number (Ha), which means stronger magnetic field, reduced primary velocity along the plates and increased secondary velocity due to Hall current, temperature value was reduced.Kuznetsov [2] performed an analytical solution to Couette flow in a composite channel, partially filled with a porous medium and partially with a clear fluid.The lower plate is moving and transfer constant heat flux, while the upper plate is fixed and adiabatic.It was shown that temperature value increases from the lower plate to the upper plate while the temperature gradient is decreases.Buonomo et al. [3] investigated fully developed, steady state forced convection in parallel-plate micro channels filled with a porous medium saturated with rarefied gases at high temperatures in local thermal non-equilibrium (LTNE) condition for the first-order slip-flow regime ( ) . A uniform heat flux was applied at both walls of the parallel-plate channel.The result showed that the internal dissipation increases as the velocity slip increases.Xu [4] investigated theoretically the thermal performance of multi-layered micro heat exchangers with porous media.It was shown that counter-flow heat exchanger is much more efficient than parallel-flow heat exchanger, and the micro heat exchanger with porous medium is preferable than the heat exchanger without a porous medium, especially for counter-flow heat exchanger.Makinde and Onyejekwe [5] investigated the steady flow and heat transfer of an electrically conducting fluid with variable viscosity and electrical conductivity between two parallel plates in the presence of a transverse magnetic field.It was assumed that couple of effects drove the flow: pressure gradient along the plates and upper plate constant velocity (the lower plate was fixed).The plates were kept at constant but different temperatures.The work showed that increasing the viscosity dependence at temperature caused the heat transfer rate to be reduced through the standing plate and increased through the moving plate.Hatami et al. [6] showed a study of magneto-hydrodynamic (MHD) Couette flows between two parallel infinite plates.The fluid between the plates was characterized as electrically conducting with particles suspension.The two phases: fluid phase and particle phase were tested.The upper plate was fixed and lower plate was time dependent.It was shown that when the magnetic field was fixed to the moving plate, increasing the intensity of the magnetic field (increasing Hartmann number) increases velocities of both phases, while the magnetic field was fixed to the standing plate, an inverse trend was received.Attia [7] investigated the MHD flow and heat transfer of a dusty and electrically conducting fluid in the presence of a uniform magnetic fluid while the variations of the viscosity and the electric conductivity of the flu-S.Sadik DOI: 10.4236/epe.2018.109026416 Energy and Power Engineering id with temperature was taken under consideration.The governing equations were solved numerically using finite difference method.It was shown that increasing the fluid viscosity, increases the temperatures and velocities of the fluid and the dust particles.A similar investigation was performed by Eguia et al. [8] using network simulation method (NSM).Same results were received.Abd-Alla et al. [9] investigated the effect of both rotation and magnetic field of a micro polar fluid through a porous medium induced by sinusoidal peristaltic waves travelling down the channel walls.It was shown that the pressure gradient increases with increasing the rotation while it decreases with increasing the Hartmann number.It may be concluded that the rotation oppose the fluid flow while the magnetic field support the flow.Makinde and Chinyoka [10]  Reese study it did not cover all Knudsen possible range.It was shown that the temperature distribution in slip flow is higher than that in non-slip flow.
The above works do not deal with the effect of temperature fluctuations on the two surfaces that surround the fluid in a Couette-Poisseuille flow on the fluid convection ability or on its ability to remove heat.The current work deals with oscillating temperature constraints on both plates.The temperature oscillations may be differ by amplitude and frequency.Heat dissipation due to shear stresses was taken under consideration.The current physical model was introduced as a function of dimensional parameters since it was wished to focus the results on truly parameters.

Couette-Poisseuille Flow
Couette-Poiseuille flow is a steady, one-dimensional flow between two plates with constant gap; the flow is along the plates or along the x direction.Poiseuille flow structure and schematic velocity distribution is shown in Figure 1.
Poiseuille velocity distribution is parabolic and given as:  is the piezo metric pressure gradient, for horizontal plates the piezo metric pressure gradient is equal to the pure pressure gradient.More information about Couette-Poisseuille flow may be found for example in Fox et al. [13].

The Energy Equation
The energy equation for a fluid is given as: where ρ is the fluid density, v c is the fluid specific capacity in a constant pressure, DT Dt is the temperature material derivative, k is the fluid conductivity, 2 T ∇ is the temperature Laplacian, µ is the dynamic viscosity and v φ is the dissipation function.

The Energy Model Equation
For Poiseuille flow, the dissipation value is received as: The energy equation for two-dimensional coordinates is written as: ( ) In Couette-Poiseuille flow there is no velocity component in the y direction, 0 y v = and there is no changes in the x direction, 0 T x ∂ ∂ = .The energy equa- tion will be received as: or as: By marking: and: the energy equation will be received as: ( ) or as: ( ) by marking: ( ) it is received that: ( ) ( ) With θ the energy equation will be received as: According to the above definitions (Equation ( 12)), the following relations are received: ( ) or: Depending on θ only, the energy equation is received as:  where the following relations are received: Finally, the energy equation is received as: The last equation for oscillating solution is received (Sadik [14]) as:

The Current Model-Basic Structure
The basic model structure is shown in Figure 2.

The Current Model-Stage 1
The first stage structure of the current model is shown in Figure 3.The first boundary condition leads to the following relation: or to: The second boundary condition leads to the following connection: or after dividing Equation (33) by it leads to: It is received that: By inserting 2 c value from Equation (32) into Equation (35), the following relations are received: Full solution for φ is θ is , for T the full solution is received as ( ) For T, the initial condition is

The Current Model-Stage 2
The second stage structure of the current model is shown in Figure 4.
The boundary conditions for this stage may written as: The first boundary condition leads to the following connections: The second boundary condition leads to the following relations: Equation ( 43) may be divided by  In a parallel way to the first stage, Full solution for φ is and for T the full solution is received as For T, the initial condition is

The Current Model-Stage 3
Stage 3 is a simple addition of the results received at stage 1 and 2.

Results and Discussion
,           fluctuations were obtained without a fluid media (Sadik [14] and [15]).This result can be applied, for example, if it is necessary to remove heat by convection from an electronic chip that leads current with high frequency fluctuations or with high frequency of Ohmic power which causes to high temperature oscillations, and there is a limit on the convection liquid thickness, the thickness of the liquid needed may be reduced.

Conclusion
Main conclusion referred to the reduced penetration intensity of the enlarged frequency constraint.When the constraint frequency increases, its penetration

Declare
This work is dedicated to the memory of my mother and mother-in-law.

Figure 2 .
Figure 2. The current model-basic structure.

Figure 3 .
Figure 3.The current model-first structure.

,
the first boundary condition may be written with its real value as and related to T the first boundary condition is received as ( )

Figure 4 .
Figure 4.The current model second structure.

Figure 5
Figure 5 shows the functions ( ) y φ

Figure 5 (
b) shows change of temperature over time.It is shown that fluid temperature increases with time progression compare to Figure 5(a) plots.Without heat dissipation the phi function would represent the temperature distribution as is shown in Figure 5(a).The temperature increasing in Figure 5(b) over the temperature area in Figure 5(a) is due to the dissipation function.All the plot lines in Figure 5 are symmetric related to the ŷ axis, according to identity of phi amplitudes and frequencies.What was written above can be shown clearly by following the black line in Figure 5(b) ( ) 13T t .The temperatures, which the black line represent, are highest according to maximum heat dissipation time.Its parallel line in the phi plot (Figure 5(a)) is the red line ( ) can be considered as a reference line for the black line in Figure 5(b).Following the black line in Figure 5(b) shows maximum temperatures values near the plates where maximum shear stresses are produced and smaller temperature in the center according to extremum point area and smaller shear stresses.

Figure 7
Figure 7 shows ( ) y φ Figure 7 shows referring to both Figure 7(a) and Figure 7(b) the larger constraint amplitude penetration compare to the smaller amplitude penetration.By following every plot line, it may be concluded and seen that fluid temperatures follow the constraint temperature.By following for example the light blue line in Figure 7(a) it is shown that larger absolute temperatures values are received near the lower plate where the amplitude constraint is larger.By following the black line in Figure 7(b) where heat dissipation is also considered, the same diagnosis is shown, the temperature values near the lower plate are larger than the temperatures values near the upper plate.

Figure 11
Figure 11 shows ( ) y φ capacity is reduced.Other conclusion relates to gap width between plates; larger gap width increases fluid temperatures due to fluid increased velocity.Referring to fluid thermal parameters, increasing thermal fluid conductivity or decreasing fluid density or decreasing heat specific capacity, increases constrains temperatures penetration intensity; this is due to enlarged thermal fluid diffusivity.Some expected results are shown: increasing gradient pressure leads to enlarged fluid velocity and more heat dissipation resulted in larger temperatures values into the fluid, increasing temperature constraint amplitude, and increasing heat penetration that leads to enlarged temperature into the fluid.