Next Article in Journal
Reconstructing the Unknown Source Function of a Fractional Parabolic Equation from the Final Data with the Conformable Derivative
Previous Article in Journal
A Novel Deep Learning Model for Sea State Classification Using Visual-Range Sea Images
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Solution of Radiative and Viscous Dissipative Fluid Flow along an Oscillating Vertical Plate

by
Rajaraman Rengasamy
1,* and
Rajamanickam Muthucumaraswamy
2
1
Department of Mathematics, R.M.K. Engineering College, Kavaraipettai 601206, India
2
Department of Applied Mathematics, Sri Venkateswara College of Engineering, Sriperumbudur 602117, India
*
Author to whom correspondence should be addressed.
Symmetry 2022, 14(7), 1488; https://doi.org/10.3390/sym14071488
Submission received: 18 June 2022 / Revised: 14 July 2022 / Accepted: 18 July 2022 / Published: 21 July 2022

Abstract

:
A numerical analysis of an unsteady radiative and viscous dissipative fluid flow across a semi-infinite oscillating vertical plate with a constant temperature and mass diffusion is examined here. The fluid under consideration is optically thin gray, emitting, heat-retaining, and non-dispersing. A finite-difference approach known as the Crank−Nicolson strategy was applied to find the solution to the limitless governing equations. Applications of such engineering problems can be found in fields such as aerospace, solar power, the cooling of nuclear reactors, and chemical and mechanical engineering. Based on the computational solutions, the impacts of distinct physical parameters, such as velocity, concentration, temperature, skin friction, Nusselt number, and Sherwood number profiles, are portrayed pictorially and imparted. The numerical solution of the velocity profile developed in this study fits extremely well with what was previously published and validated. Furthermore, we discovered that plate oscillation, radiation, and viscous dissipation parameters surprisingly influence the flow pattern.

1. Introduction

Free convection or natural convection is caused by density differences due to temperature variations in the fluid. Free convection flow has several practical applications in agriculture, oceanography, and industry. Due to its extensive range of applications, several articles have been published on the temperamental natural convection progression of a viscid incompressible fluid alongside a semi-boundless upward plate. Soundalgekar [1] has demonstrated that the free convection impacts on the stream over an upward swaying plate can be deduced. He also established that when the plate cools quickly, the velocity of the stream increases, and when the plate heats quickly, the velocity of the stream drops. The transient free convection stream with intermittent temperature variation of a boundless upward plate has been addressed by Das et al. [2].
Scientists have been fascinated by the aftermath of radiative heat and mass transfer due to their general applications in aerospace, solar power, cooling nuclear reactors, and synthetic materials. Using the Rosseland transition model, Hossain and Takhar [3] examined the radiation consequences of heterogeneous convection on an upward plate with a uniform surface temperature. Muthucumaraswamy [4] analyzed the radiative effects of an erratic shear convective stream over a boundless upward swaying plate with a varying temperature and mass dispersion. By the Laplace transform method, he obtained an exact solution to the non-dimensional governing equations. Muthucumaraswamy and Vijayalakshmi [5] investigated the effect of mass exchange and heat transfer across an oscillating upward plate with varying temperatures. The numerical solution provided by Raptis [6] describes the convectional oscillatory flow and mass exchange for a fluid with a thin optical layer lying above a porous plate within the sight of radiation. Prasad et al. [7] examined the impact of thermal radiation on a natural convective viscous incompressible stream travelling over a semi-boundless upward plate with constant mass and heat flux.
Fluid flows periodically illustrate symmetricity under certain conditions. However, such a symmetric flow is not always attained if circumstances change. Because of their importance and relevance in flow control and heat transfer enhancement, the breaks of flow symmetry have been explored theoretically, experimentally, and numerically in fluid dynamics and thermal engineering. The conversion of energy from one form to another is called dissipation in fluid dynamics. The amount of dissipation is high in regions with large gradients as well as in regions with high levels of turbulence. A temperature variation may affect the stream characteristics due to the liquid thickness based on the viscous dissipation’s liquid property. Viscous dissipation is a significant factor in geophysical streams and food processing. A semi-limitless vertical isothermal plate by Gebhart [8] was the first to demonstrate the significance of viscous dissipation in free convection streams. According to Mahajan and Gebhart [9], viscous heating dissipation impacts natural convection. Furthermore, they showed that viscous dissipation effects are negligible compared to pressure effects for all Prandtl numbers (Pr). Soundalgekar et al. [10] illustrated the effect of an erratic shear convection stream of a viscous dissipative fluid along a semi-boundless isothermal plate. Muthucumaraswamy and Saravanan [11] examined the impacts of radiation-induced waves on streams passing through a swaying semi-boundless upward plate with invariant mass and varying temperatures. Mondal et al. [12] developed a numerical study of radiated and chemical reactions associated with the fluid flow across a porous medium with viscous dissemination. The motion of an impenetrable viscous fluid over a permeable upward plate beneath the significance of radiation and viscous dissemination was studied by Das et al. [13]. Dharmendar Reddy and Srinivasa Rao [14] investigated the influence of viscous dissipation in the presence of a heat source/sink on an unsteady, two-dimensional laminar convective boundary layer flow of a viscous, incompressible, chemically reactive fluid along a semi-infinite vertical plate with suction. They applied the finite element technique (FEM) for solving the governing boundary layer equations. Hemamalini and Suresh Kumar [15] examined the impacts of viscous dispersion and warm radiation on a swaying upward plate in a penetrable medium.
Regardless, no expert has suggested that a numerical study be conducted of a transient free convective stream passing through a semi-limitless swaying upward plate in the existence of heat radiation and viscous dissemination. This study employs an implicit finite-difference approach to examine an upward oscillating semi-infinite plate with a stable mass and temperature distribution considering the thermal radiation and viscous dissipation effects.

2. Mathematical Analysis

We propose to inspect the unstable laminar convective flow of an incompressible viscid liquid along an oscillating semi-boundless upward plate beneath the prominence of warm radiation and viscous dissemination. In this case, the plate is considered to be vertically aligned along the x- a x i s , with the y- a x i s being normal to the plate as shown in Figure 1. The concept of symmetry is quite apparent in the boundary layer theory. Here, the boundary is treated as a wall boundary, and the x- a x i s is considered as a line of symmetry. Both the liquid and plate are at an identical concentration C and temperature T at the beginning. The motion begins at the moment t > 0 , and the plate sways in its distinguished notch with velocity u 0 c o s ω t while the concentration and temperature of the plate are simultaneously uplifted to C w and T w , respectively. For the energy equation, viscous dissipation is taken into account. According to the above assumptions, the boundary layer equations for continuity, momentum, energy, and mass diffusion are presented below with the customary Boussinesq’s approximation (Gebart and Pera [16]):
u x + u y = 0
u t + u u x + v u y = g β ( T T ) + g β ( C C ) + υ u y y
T t + u T x + v T y = k ρ C p T y y 1 ρ C p ( q r ) y + υ ρ C p ( u y ) 2
C t + u C x + v C y = D C y y
where u   and   v are velocity factors in the x   and   y path, respectively, t —time, g —acceleration due to gravity, β * and   β —coefficient of volumetric concentration and thermal expansion, respectively, k —thermal conductivity, υ —kinematic viscosity, ρ —fluid density, μ —coefficient of viscosity, C p —specific heat at constant pressure, C —species concentration in the fluid, T —temperature of the fluid near the plate, and C and T —concentration and temperature at the free stream.
The associated boundary premises are as follows:
t 0 :   u = 0 ,   v = 0 ,   T = T ,   C = C t > 0 :   u = u 0 cos ω t ,   v = 0 ,   T = T w ,   C = C w     a t     y = 0 u = 0 ,   T = T ,   C = C     a t     x = 0 u 0 ,   T T ,   C C     a s     y }
where u 0 —velocity of the plate, ω t —phase angle, and C w   and   T w —concentration and temperature near the wall.
By employing the Rosseland approximation [17], the radiative heat flux q r is given by the following:
( q r ) y = 4 a σ s ( T 4 T 4 )
where a * —absorption constant, and σ s —Stefan−Boltzmann constant.
It should be noted that by using the Rosseland approximation, the present analysis is limited to optically thick fluids. The refractive index of the gas medium is constant. The local thermal radiation intensity is due to radiation emanating from proximate locations in the vicinity of which emission and scattering are comparable to the location of interest. The energy transfer depends on conditions only in the area adjacent to the plate regime, i.e., the boundary layer regime. If the temperature differences within the flow are sufficiently small, then Equation (6) can be linearized by expanding T 4 in a Taylor series about T and ignoring the terms of higher order, which results in the following:
T 4 4 T 3 T 3 T 4
Applying Equations (6) and (7), in Equation (3), we arrive at the following:
T t + u T x + v T y = k ρ C p ( T y y ) + 16 a σ s T 3 ( T T ) ρ C p + υ ρ C p ( u y ) 2
The following boundless quantities are identified by using the boundary layer equations and conditions:
X = x u 0 υ ,   Y = y u 0 υ ,   U = u u 0 ,   V = v u 0 ,   t = t u 0 2 υ ,   ω = ω υ u 0 2 , T = T T T w T ,   C = C C C w C ,   G r = υ g β ( T w T ) u 0 3 ,   G c = υ g β ( C w C ) u 0 3 , Pr = υ α ,   S c = υ D ,   E c = u 0 2 C p ( T w T ) ,   R = 16 a υ 2 σ s T 3 k u 0 2 }
where “ X —boundless spatial coordinate along the plate”, “ Y —boundless spatial coordinate normal to the plate”, U ,   V —boundless velocity components in X ,   Y —direction respectively, t —boundless time, ω —frequency of oscillation, T —boundless temperature, C —boundless concentration, “ G r —thermal Grashof number”, “ G c —mass Grashof number”, “ P r —Prandtl number”, “ S c —Schmidt number”, “ E c —Eckert number”, “ R —radiation parameter”, α —thermal diffusivity, and D —mass diffusion coefficient.
After applying Equation (9) in Equations (1)–(4), the following boundless equations are obtained:
U X + V Y = 0
U t + U U X + V U Y = G r T + G c C + U Y Y
T t + U T X + V T Y = 1 Pr ( T Y Y ) R T Pr + E c ( U Y ) 2
C t + U C X + V C Y = 1 S c C Y Y
The associated boundless boundary premises are as follows:
t 0 :   U = 0 ,   V = 0 ,   T = 0 ,   C = 0 t > 0 :   U = cos ω t ,   V = 0 ,   T = 1 ,   C = 1     a t   Y = 0 U = 0 ,   T = 0 ,   C = 0     a t   X = 0 U 0 ,   T 0 ,   C 0 ,   a s     Y }
where ω t —boundless phase angle.
The boundless form of the local and average skin friction is shown here:
τ X = ( U Y ) Y = 0   and     τ ¯ = 0 1 ( U Y ) Y = 0 d X
The boundless form of the local and average Nusselt number is as follows:
N u X = X [ ( T Y ) Y = 0 / T Y = 0 ]   and     N u ¯ = 0 1 [ ( T Y ) Y = 0 / T Y = 0 ] d X
The boundless form of the local and average Sherwood number is as follows:
S h X = X [ ( C Y ) Y = 0 / C Y = 0 ]   and     S h ¯ = 0 1 [ ( C Y ) Y = 0 / C Y = 0 ] d X

3. Numerical Solution of the Problem

The coupled non-linear partial differential Equations (10)–(13) under Equation (14) are converted into finite-difference equations by the Crank−Nicolson method.
[ U i , j n + 1 U i 1 , j n + 1 + U i , j n U i 1 , j n + U i , j 1 n + 1 U i 1 , j 1 n + 1 + U i , j 1 n U i 1 , j 1 n ] 4 Δ X + [ V i , j n + 1 V i , j 1 n + 1 + V i , j n V i , j 1 n ] 2 Δ Y = 0
[ U i , j n + 1 U i , j n ] Δ t + U i , j n [ U i , j n + 1 U i 1 , j n + 1 + U i , j n U i 1 , j n ] 2 Δ X + V i , j n [ U i , j + 1 n + 1 U i , j 1 n + 1 + U i , j + 1 n U i , j 1 n ] 4 Δ Y = G r 2 [ T i , j n + 1 + T i , j n ] + G c 2 [ C i , j n + 1 + C i , j n ] + [ U i , j 1 n + 1 2 U i , j n + 1 + U i , j + 1 n + 1 + U i , j 1 n 2 U i , j n + U i , j + 1 n ] 2 ( Δ Y ) 2
[ T i , j n + 1 T i , j n ] Δ t + U i , j n [ T i , j n + 1 T i 1 , j n + 1 + T i , j n T i 1 , j n ] 2 Δ X + V i , j n [ T i , j + 1 n + 1 T i , j 1 n + 1 + T i , j + 1 n T i , j 1 n ] 4 Δ Y = 1 Pr [ T i , j 1 n + 1 2 T i , j n + 1 + T i , j + 1 n + 1 + T i , j 1 n 2 T i , j n + T i , j + 1 n ] 2 ( Δ Y ) 2 R 2 Pr [ T i , j n + 1 + T i , j n ] + E c [ U i , j + 1 n U i , j n Δ Y ] 2
[ C i , j n + 1 C i , j n ] Δ t + U i , j n [ C i , j n + 1 C i 1 , j n + 1 + C i , j n C i 1 , j n ] 2 Δ X + V i , j n [ C i , j + 1 n + 1 C i , j 1 n + 1 + C i , j + 1 n C i , j 1 n ] 4 Δ Y = 1 S c [ C i , j 1 n + 1 2 C i , j n + 1 + C i , j + 1 n + 1 + C i , j 1 n 2 C i , j n + C i , j + 1 n ] 2 ( Δ Y ) 2
The detailed description of the functions and interrelationship connecting Equations (18)–(21) is referred to by John D. Anderson, Jr. [18]. At each interior nodal point on a specific i-level, the outputs of U ,   V , T   and   C are obtained from the finite-difference Equations (18)–(21), forming a tri-diagonal system of equations that is solved by the Thomas algorithm as discussed by Carnahan et al. [19]. A rectangular mesh is assumed here with sides X m a x = 1   and   Y m a x = 10 where Y m a x corresponds to Y = or those that lie beyond the thermal and momentum boundary layers. The limit of Y was picked as 10 after some fundamental mathematical calculations to such an extent that the last two of the limit conditions of Equation (14) are fulfilled inside as far as possible at 10 5 . The mesh sizes were set to Δ X = 0.05 ,   Δ Y = 0.125 with a time period of “ Δ t = 0.01 .

4. Graphical Results and Discussion

We used the Windows version of FORTRAN (Force 2.0) to compute the numerical values of U ,   V , T   and   C using the iterative procedure, and the outcomes of the numerical results were plotted graphically using MATLAB software. The steady-state computational values of velocity, concentration, temperature, local and average skin friction, Nusselt number, and Sherwood number were computed using physical parameters such as the phase angle “t = π/6”, thermal radiation “R = 2”, Grashof numbers “Gr = 5, Gc = 5”, Prandtl number “Pr = 0.71 (air)”, Schmidt number “Sc = 0.6 (water vapour)”, and Eckert number “Ec = 0.5”. All figures are based on these values unless specified otherwise.
Figure 2 illustrates the discrimination of the present velocity profile against the exact solution unveiled by Muthucumaraswamy [4]. The parameters “R = 0.2, 5, 20”, “ωt = π/6”, “Gr = 5”, “Gc = 2”, “Pr = 0.71”, “Sc = 0.3”, and “t = 0.4” are used for comparison, and the results are compatible with the exact solution. Figure 3a depicts velocity profiles for various phase angle values. The velocity is found to decrease as the phase angle increases. It can also be seen that velocity diminishes in the presence of high wavering. Furthermore, it has been noted that the time required to attain the consistent state on a vertical plate “ωt = 0” is much greater than on a flat plate “ωt = π/2”.
A consistent velocity profile for thermal and mass Grashof numbers “Gr” and “Gc” can be seen in Figure 3b. Rising the thermal Grashof number (or mass Grashof number) will result in more incredible air velocity. In the context of mass buoyancy, inflation in thermal and mass Grashof numbers increases the effect of mass buoyancy. The velocity is revamped by essential plate cooling. The steady-state velocity pattern corresponding to a range of the radiation parameter “R” is exhibited in Figure 3c. It was noticed that increasing the radiation parameter resulted in a decrease in velocity.
Figure 3d depicts the impact of velocity patterns for various Schmidt number “Sc” values. It has been found that velocity diminishes while expanding the Schmidt number. Furthermore, the time it takes to attain a consistent state depends on the Schmidt number. In addition, the limit layer is formed along the major end of a semi-endless upward plate moving in a liquid and the plate’s direction. According to the results, maximum overall velocity is produced by mass dispersion and buoyancy force. The consistent state velocity pattern for the various range of the Eckert number “Ec” is represented in Figure 3e. The graph clearly illustrates that velocity increases as the Eckert number increases.
Figure 4 depicts the impact of various Schmidt numbers on the concentration profile. It has been noticed that the concentration of plates diminished as the Schmidt number value was increased. Generally, it is true that an ascent in ‘Sc’ firmly squashes the molecular diffusivity. The temperature profile of various radiation parameters is depicted in Figure 5a. It can be identified that extending the radiation boundary “R” brings about a lessening of the velocity just as it does the temperature inside the limit layer. The temperature drops dramatically from the plate peak value at “y = 0” to the free stream “y → ∞” as the radiation parameter “R” increases. The Eckert number “Ec” has a significant impact on temperature, as can be seen in Figure 5b. The temperature rises as the dissemination boundary is extended. When viscous dissipation takes place in the stream field, there is an increase in energy, causing a rise in the temperature of the liquid and, therefore, in buoyancy.
Figure 6 demonstrates the local skin friction for various values of the phase angle “ωt” as an element of the pivotal coordinate “X”. The shear stress local wall increments with diminishing upsides of stage point “ωt”. As a result, an upward plate “ωt = 0” will have more significant wall shear stress than a flat plate “ωt = π/2”. Figure 7 illustrates the impacts of various phase angles on the average skin friction values over the time “t” at “X = 1.0”. When the phase angle “ωt” is expanded, the value of average skin friction diminishes. The impact of the Eckert number “Ec” over the local Nusselt number is represented in Figure 8 as a component of the pivotal coordinate “X”. The heat transfer rate generally increases as the dissipative parameter “Ec” decreases. The average Nusselt number as a function of time “t” at “X = 1.0” is shown in Figure 9. As the dissipative parameter “Ec” diminishes, the average Nusselt number increases. A local Sherwood number is shown in Figure 10 for various upsides of the Schmidt number “Sc”. When the Schmidt number “Sc” increases, the rate of mass transfer increases. On the other hand, the concentration profile on Figure 4 reverses this pattern. Figure 11 demonstrates the relation between the average Sherwood number and the Schmidt number. As the Schmidt number rises, the average Sherwood number increases as well.

5. Conclusions

In the existence of warm radiation and viscous dissipation, an upward swaying semi-endless plate with a constant mass and temperature is investigated. Analyzing the results for several parameters is depicted in the graph. The current study reveals the following intriguing features:
  • The velocity decreases when the phase angle for air in the cooling of the plate is raised.
  • Both the velocity and temperature fall while the radiation parameter increases.
  • Temperature falls while expanding the radiation parameter but rises when increasing the Eckert number.
  • Local and average skin friction decrease as the phase angle increases.
  • Eckert numbers enhance velocity and temperature profiles but diminish local and average Nusselt numbers.
  • Local and average Sherwood numbers are enhanced by an expanding Schmidt number, but in concentration, the pattern reverses.
  • According to the results of this study, the number of steps required to achieve a consistent state is highly influenced by the radiation parameter and phase angle.

Author Contributions

R.M. designed the research and formulated the mathematical model; R.R. computed the numerical results and wrote the manuscript; R.R. and R.M. analyzed the results. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to express their gratitude to the anonymous reviewers for their valuable comments and suggestions to improve this article and make it more distinguished.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Soundalgekar, V.M. Free convection effects on the flow past a vertical oscillating plate. Astrophys. Space Sci. 1979, 64, 165–172. Available online: https://ui.adsabs.harvard.edu/link_gateway/1979Ap&SS..64..165S/ (accessed on 17 June 2022). [CrossRef]
  2. Das, U.N.; Deka, R.K.; Soundalgekar, V.M. Transient free convection flow past an infinite vertical plate with periodic temperature variation. ASME J. Heat Transf. 1999, 121, 1091–1094. [Google Scholar] [CrossRef]
  3. Hossain, M.A.; Takhar, H.S. Radiation effects on mixed convection along a vertical plate with uniform surface temperature. J. Heat Mass Transf. 1996, 31, 243–248. [Google Scholar] [CrossRef]
  4. Muthucumaraswamy, R. Radiative flow past an infinite vertical oscillating plate with variable temperature and mass diffusion. Int. J. Math. Anal. 2006, 4, 35–46. [Google Scholar]
  5. Muthucumaraswamy, R.; Vijayalakshmi, A. Effects of heat and mass transfer on flow past an oscillating vertical plate with variable temperature. Int. J. Appl. Math. Mech. 2008, 4, 59–65. [Google Scholar]
  6. Raptis, A. Free convective oscillatory flow and mass transfer past a porous plate in the presence of radiation for an optically thin fluid. Therm. Sci. 2011, 15, 849–857. [Google Scholar] [CrossRef]
  7. Prasad, V.R.; Vasu, B.; Bhaskar Reddy, N. Radiation and Mass transfer effects on transient free convection flow of a dissipative fluid past semi-infinite vertical plate with uniform heat and mass flux. J. Appl. Fluid Mech. 2011, 4, 15–26. [Google Scholar] [CrossRef]
  8. Gebhart, B. Effects of viscous dissipation in natural convection. J. Fluid Mech. 1962, 14, 225–232. [Google Scholar] [CrossRef]
  9. Mahajan, R.I.; Gebhart, B.B. Viscous dissipation effects in buoyancy induced flows. Int. J. Heat Mass Transf. 1989, 32, 1380–1382. [Google Scholar] [CrossRef]
  10. Soundalgekar, V.M.; Jaiswal, B.S.; Uplekar, A.G.; Takhar, H.S. Transient free convection flow of a viscous dissipative fluid past a semi-infinite vertical plate. Appl. Mech. Eng. 1999, 4, 203–218. [Google Scholar]
  11. Muthucumaraswamy, R.; Saravanan, B. Finite difference solution of unsteady flow past an oscillating semi-infinite vertical plate with variable surface temperature and uniform mass flux. Int. J. Appl. Mech. Eng. 2014, 19, 709–724. [Google Scholar] [CrossRef]
  12. Mondal, R.K.; Hossain, M.A.; Ahmmed, S.F. Radiation and Chemical reaction effects on free convection and mass transfer flow of dissipative fluid past an exponentially infinite vertical plate through a porous medium. Elixir Int. J. 2015, 84C, 33523–33530. [Google Scholar]
  13. Das, S.; Jana, M.; Jana, R.N. Effects of radiation and viscous dissipation on unsteady free convective flow past a moving vertical porous plate embedded in porous medium. Commun. Appl. Sci. 2013, 1, 59–80. [Google Scholar]
  14. Dharmendar Reddy, Y.; Srinivasa Rao, V. Finite Element Solution of Thermal Radiation and Mass Transfer Flow past Semi- infinite Moving Vertical Plate with Viscous Dissipation in Presence of heat Source/Sink. J. Appl. Math. Fluid Mech. 2016, 8, 15–31. [Google Scholar]
  15. Hemamalini, P.T.; Suresh Kumar, N. The effects of thermal radiation and Viscous dissipation on heat and mass diffusion flow past an oscillating infinite vertical plate with variable temperature embedded in a porous medium. Int. J. Creat. Res. Thoughts 2018, 6, 1165–1173. Available online: https://ijcrt.org/papers/IJCRT1892173.pdf (accessed on 17 June 2022).
  16. Gebhart, B.; Pera, L. The nature of vertical natural convection flows resulting from the combined buoyancy effects of thermal and mass diffusion. Int. J. Heat Mass Transf. 1971, 14, 2025–2050. [Google Scholar] [CrossRef]
  17. Brewster, M.Q. Thermal Radiative Transfer and Properties; John Wiley & Sons: New York, NY, USA, 1992. [Google Scholar]
  18. Anderson, J.D., Jr. Computational Fluid Dynamics the Basics with Applications; McGraw-Hill, Inc.: New York, NY, USA, 1995. [Google Scholar]
  19. Carnahan, B.; Luther, H.A.; Wilkes, J.O. Applied Numerical Methods; John Wiley & Sons: New York, NY, USA, 1969. [Google Scholar]
Figure 1. Physical model of the problem.
Figure 1. Physical model of the problem.
Symmetry 14 01488 g001
Figure 2. Comparison of velocity profiles.
Figure 2. Comparison of velocity profiles.
Symmetry 14 01488 g002
Figure 3. Velocity profiles for different values of (a) ωt, (b) Gr, Gc, (c) R, (d) Sc, and (e) Ec.
Figure 3. Velocity profiles for different values of (a) ωt, (b) Gr, Gc, (c) R, (d) Sc, and (e) Ec.
Symmetry 14 01488 g003aSymmetry 14 01488 g003b
Figure 4. Concentration profile for different values of Sc.
Figure 4. Concentration profile for different values of Sc.
Symmetry 14 01488 g004
Figure 5. Temperature profiles for different values of (a) R and (b) Ec.
Figure 5. Temperature profiles for different values of (a) R and (b) Ec.
Symmetry 14 01488 g005
Figure 6. Local skin friction for different values of ωt.
Figure 6. Local skin friction for different values of ωt.
Symmetry 14 01488 g006
Figure 7. Average skin friction for different values of ωt.
Figure 7. Average skin friction for different values of ωt.
Symmetry 14 01488 g007
Figure 8. Local Nusselt number for different values of Ec.
Figure 8. Local Nusselt number for different values of Ec.
Symmetry 14 01488 g008
Figure 9. Average Nusselt number for different values of Ec.
Figure 9. Average Nusselt number for different values of Ec.
Symmetry 14 01488 g009
Figure 10. Local Sherwood number for different values of Sc.
Figure 10. Local Sherwood number for different values of Sc.
Symmetry 14 01488 g010
Figure 11. Average Sherwood number for different values of Sc.
Figure 11. Average Sherwood number for different values of Sc.
Symmetry 14 01488 g011
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rengasamy, R.; Muthucumaraswamy, R. Numerical Solution of Radiative and Viscous Dissipative Fluid Flow along an Oscillating Vertical Plate. Symmetry 2022, 14, 1488. https://doi.org/10.3390/sym14071488

AMA Style

Rengasamy R, Muthucumaraswamy R. Numerical Solution of Radiative and Viscous Dissipative Fluid Flow along an Oscillating Vertical Plate. Symmetry. 2022; 14(7):1488. https://doi.org/10.3390/sym14071488

Chicago/Turabian Style

Rengasamy, Rajaraman, and Rajamanickam Muthucumaraswamy. 2022. "Numerical Solution of Radiative and Viscous Dissipative Fluid Flow along an Oscillating Vertical Plate" Symmetry 14, no. 7: 1488. https://doi.org/10.3390/sym14071488

APA Style

Rengasamy, R., & Muthucumaraswamy, R. (2022). Numerical Solution of Radiative and Viscous Dissipative Fluid Flow along an Oscillating Vertical Plate. Symmetry, 14(7), 1488. https://doi.org/10.3390/sym14071488

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop