Lattice Boltzmann simulation of convective heat transfer of non-Newtonian fluids in impeller stirred tank

Lattice Boltzmann simulation of convective heat transfer of non-Newtonian fluids in an impeller stirred tank is performed. The curved and moving boundary methods combined with the unknown-index algorithm are used to solve the flow and thermal fields induced in a cold tank by an oscillating hot impeller. For a given maximum radius of the blades, the simulation results show that a rectangular impeller of large aspect ratio induces stronger heat transfer effect on the tank walls than the small aspect ratio. This is because the latter would cause worse field synergy than the former, i.e. the induced local velocities of fluid are mostly perpendicular to the temperature gradients. The convection effects on the tank walls are also improved as the oscillation amplitude of impeller increases until the swept areas of impeller are close to whole azimuth of the tank, i.e., oscillation amplitude of 90°. The maximum Nusselt number on the tank walls for power-law fluid flows of n = 0.7, 1 and 1.5 occurs at oscillation amplitude of 75°. Finally, it is found that the heat transfer effect on the tank walls is reduced as the power-law index of fluid increases.     

The influence of heat transfer on peristaltic transport of a Jeffrey fluid in a vertical porous stratum

The peristaltic flow of a Jeffrey fluid in a vertical porous stratum with heat transfer is studied under long wavelength and low Reynolds number assumptions. The nonlinear governing equations are solved using perturbation technique. The expressions for velocity, temperature and the pressure rise per one wave length are determined. The effects of different parameters on the velocity, the temperature and the pumping characteristics are discussed. It is observed that the effects of the Jeffrey number λ₁, the Grashof number Gr, the perturbation parameter N = EcPr, and the peristaltic wall deformation parameter ϕ are the strongest on the trapping bolus phenomenon. The results obtained for the flow and heat transfer characteristics reveal many interesting behaviors that warrant further study on the non-Newtonian fluid phenomena, especially the shear-thinning phenomena. Shear-thinning reduces the wall shear stress.     

Numerical simulation of magnetohydrodynamic buoyancy-induced flow in a non-isothermally heated square enclosure

Numerical simulation of magnetohydrodynamic (MHD) buoyancy-induced heat transfer and fluid flow has been analyzed in a non-isothermally heated square enclosure using finite volume method. The bottom wall of enclosure were heated and cooled with a sinusoidal function and top wall was cooled isothermally. Vertical walls of the enclosure were adiabatic. Effects of Rayleigh number (Ra = 10⁴, 10⁵ and 10⁶), Hartman number (Ha = 0, 50 and 100) and amplitude of sinusoidal function (n = 0.25, 0.5 and 1) on temperature and flow fields were analyzed. It was observed that heat transfer was decreased with increasing Hartmann number and decreasing value of amplitude of sinusoidal function.     

Effect of variable density on hydromagnetic mixed convection flow of a non-Newtonian fluid past a moving vertical plate

The effect of temperature-dependent density on MHD mixed convection flow of power-law fluid past a moving semi-infinite vertical plate for high temperature differences between the plate and the ambient fluid is studied. The fluid density is assumed to decrease exponentially with temperature. The usual Boussinesq approximations are not considered due to the large temperature differences. The surface temperature of the moving plate was assumed to vary according to a power-law form, that is, T_w(x) = T_∞ + Ax^γ. The fluid is permeated by a uniform magnetic field imposed perpendicularly to the plate on the assumption of small magnetic Reynolds number. A numerical shooting algorithm for two unknown initial conditions with fourth-order Runge–Kutta integration scheme has been used to solve the coupled non-linear boundary value problem. The effects of various parameters on the velocity and temperature profiles as well as the local skin-friction coefficient and the local Nusselt number are presented graphically and in the tabular form. The results show that application of Boussinesq approximations in a non-Newtonian fluid subjected to high temperature differences gives a significant error in the values of the skin-friction coefficient and the application of an external magnetic field reduces this error markedly in the case of shear-thickening fluid.     

Space–time symmetry violation,configuration mixing model and E-infinity theory

We have studied the time reversal symmetry violation on the bases of the configuration mixing model and E-infinity theory. With the use of the Cabibbo angle approximation, we have presented the transformation matrix in terms of the golden ratio (?), and shown that the time reversal symmetry violation is described by the configuration mixing of the unstable and stable manifolds (W^u, W^s). The magnitude of the mixing for the weak interaction field is given by the expression sin² θ_T(theor) ≈ sin⁴ θ_C(theor) ≈ (?)¹² = 3.105 × 10^?3, which is compared to the Kaon decay experiment ～2.3 × 10^?3. We have also discussed the space–time symmetry violation by using the CPT theorem.     

Transient nonlinear optically-thick radiative–convective double-diffusive boundary layers in a Darcian porous medium adjacent to an impulsively started surface: Network simulation solutions

A boundary-layer model is described for the two-dimensional nonlinear transient thermal convection heat and mass transfer in an optically-thick fluid in a Darcian porous medium adjacent to an impulsively started vertical surface, in the presence of significant thermal radiation and buoyancy forces in an (X¹,Y¹,t¹) coordinate system. An algebraic approximation is employed to simplify the integro-differential equation of radiative transfer for unidirectional flux normal to the plate into the boundary-layer regime, by incorporating this flux term in the energy conservation equation. The conservation equations are non-dimensionalized into an (X,Y,T) coordinate system and solved using the Network Simulation Method (NSM), a robust numerical technique which demonstrates high efficiency and accuracy. The transient variation of non-dimensional streamwise velocity component (u) and temperature (T) and concentration (C) functions is computed for various selected values of Stark number (radiation–conduction interaction parameter) and Darcy number. Transient velocity (u) and steady-state local skin friction (τ_X) are also studied for various thermal Grashof number (Gr), species Grashof number (Gm), Schmidt number (Sc) and Stark number (N) values. These computations for the infinite permeability case (Da → ∞) are compared with previous finite difference solutions [Prasad et al. Int J Therm Sci 2007;46(12):1251–8] and shown to be in excellent agreement. An increase in Darcy number is seen to accelerate the flow and boost velocity. A decrease in Stark number (corresponding to an increase in thermal radiation heat transfer contribution) is shown to increase the velocity values. Temperature function is observed to fall in value with a rise in Da and increase with decrease in N (corresponding to an increase in thermal radiation heat transfer contribution). Applications of the study include rocket combustion chambers, astrophysical flows, spacecraft thermal fluid dynamics in debris-laden environments (cosmic dust), heat transfer in forest fire spread, geochemical contamination and ceramic materials processing.     

Solution of the Burgers equation using an implicit linearizing transformation

The initial-boundary-value problem on the semi-infinite interval and on a finite interval for the Burgers equation u_t = u_xx + 2u_xu is solved using a stream function ? and a linearizing transformation w = e^?. The transformation reduces the equation to a heat equation with appropriate initial and homogeneous time-dependent linear boundary conditions. One advantage of this method is that we never need to find an explicit expression for ? in our computations.     

Simulación numérica del flujo turbulento de aire con gotas dispersas de agua a través de separadores de torres de refrigeración

This work presents a numerical study on the turbulent flow of air with dispersed water droplets in separators of mechanical cooling towers. The averaged Navier-Stokes equations are discretised through a finite volume method, using the Fluent and Phoenics codes, and alternatively employing the turbulence models k ? ?, k ? ω and the Reynolds stress model, with low-Re version and wall enhanced treatment refinements. The results obtained are compared with numerical and experimental results taken from the literature. The degree of accuracy obtained with each of the considered models of turbulence is stated. The influence of considering whether or not the simulation of the turbulent dispersion of droplets is analyzed, as well as the effects of other relevant parameters on the collection efficiency and the coefficient of pressure drop. Focusing on four specific eliminators (‘Belgian wave’, ‘H1-V’, ‘L-shaped’ and ‘Zig-zag’), the following ranges of parameters are outlined: 1 ≤ U_e ≤ 5 m/s for the entrance velocity, 2 ≤ D_p ≤ 50 μm for the droplet diameter, 650 ≤ Re ≤ 8.500 for Reynolds number, and 0.05 ≤ P_i ≤ 5 for the inertial parameter. Results reached alternately with Fluent and Phoenics codes are compared. The best results correspond to the simulations performed with Fluent, using the SST k ? ω turbulence model, with values of the dimensionless scaled distance to wall y⁺ in the range 0.2 to 0.5. Finally, correlations are presented to predict the conditions for maximum collection efficiency (100 %), depending on the geometric parameter of removal efficiency of each of the separators, which is introduced in this work.     

Unsteady magnetohydrodynamic Hartmann–Couette flow and heat transfer in a Darcian channel with Hall current,ionslip, viscous and Joule heating effects: Network numerical solutions

A theoretical study of unsteady magnetohydrodynamic viscous Hartmann–Couette laminar flow and heat transfer in a Darcian porous medium intercalated between parallel plates, under a constant pressure gradient is presented. Viscous dissipation, Joule heating, Hall current and ionslip current effects are included as is lateral mass flux at both plates. The dimensionless conservation equations for the primary (x¹-direction), secondary (z¹-direction) momentum and also energy conservation equation are derived and solved using a computational technique known as Network Simulation Methodology (NSM). Velocity distributions (u¹, w¹) and temperature distribution (T¹) at the channel centre (y¹ = 0) over time (t¹) are studied graphically for the effects of Darcy number (Da), Hartmann number (Ha), transpiration (N_t), Hall current parameter (Be), ionslip parameter (Bi), pressure gradient parameter (dP/dx¹) with Prandtl number prescribed at 7.0 (electrically conducting water), Eckert number held constant at 0.25 (heat convection from the plates to the fluid) and Reynolds number (Re) fixed at 5.0 (for Re < 10, Darcian model is generally valid). Increasing Darcy number causes an increase in temperature, T¹; values are however significantly reduced for the higher Hartmann number case (Ha = 10). For the case of low transpiration (i.e. N_t = 1 which corresponds to weak suction at the upper plate and weak injection at the lower plate), both primary velocity (u¹) and secondary velocity (w¹) are increased with a rise in Darcy number (owing to a simultaneous decrease in Darcian porous drag); temperature T¹ is also increased considerably with increasing Da. However, for stronger transpiration (N_t = 10), magnitudes of u¹, w¹ and T¹ are significantly reduced and also significant overshoots are detected prior to the establishment of steady state flow. With increasing Hall current parameter, Be, (for the purely fluid regime i.e. Da → ∞), primary velocity is considerably increased, whereas secondary velocity is reduced; temperatures are decreased in the early stages of flow but effectively increased in the steady state with increasing Be. With strong Darcian drag present (Da = 0.01 i.e. very low permeability), magnitudes of u¹, w¹ and T¹ are considerably reduced and temperatures are found to be reduced for all t¹, with increasing Hall current effect (Be). Increasing ionslip current parameter (Bi) increases primary velocity (u¹), decreases secondary velocity (w¹) and also temperature (T¹) for all time (t¹), in the infinite permeability case (Da → ∞). For weakly Darcian flow, ionslip parameter (Bi) has a much reduced effect on the velocity distributions. Temperature, T¹ is strongly increased with a rise in pressure gradient parameter, dP¹/dx¹, as is primary velocity (u¹); however, secondary velocity (w¹) is reduced. The present study has applications in hybrid magnetohydrodynamic (MHD) energy generators, materials processing, geophysical hydromagnetics, etc.     

Local sensitivity analysis for the heat flux-temperature integral relationship in the half-space

This paper presents insight into the heat flux-temperature (q″ ? T) integral relationship based on constant thermophysical properties. This relationship is often used in one-dimensional, transient heat transfer studies involving null-point calorimetry and heat flux investigations. This study focuses on a short transient studies where energy has not fully penetrated the body as the result of an imposed surface heating condition. A full nonlinear heat transfer model is developed involving a half-space planar region. Temperature results are then introduced into the constant property integral relationship and a newly derived Kirchoff integral relationship for retrieving the local heat flux. Good agreement is observed between the fully nonlinear results and locally linearized system. Additionally, a sensitivity study is presented which involves perturbing the average thermophysical properties of thermal conductivity and heat capacity.     

Heat transfer analysis of unsteady boundary layer flow by homotopy analysis method

This paper aims to present complete analytic solution to the unsteady heat transfer flow of an incompressible viscous fluid over a permeable plane wall. The flow is started due to an impulsively stretching porous plate. Homotopy analysis method (HAM) has been used to get accurate and complete analytic solution. The solution is uniformly valid for all time τ ∈ [0, ∞) throughout the spatial domain η ∈ [0, ∞). The accuracy of the present results is shown by giving a comparison between the present results and the results already present in the literature. This comparison proves the validity and accuracy of our present results. Finally, the effects of different parameters on temperature distribution are discussed through graphs.     

Vacuum polarization in light two-electron atoms and ions

A general theory of the vacuum polarization in light atomic and muon-atomic systems is considered. We derive the closed analytical expression for the Uehling potential and evaluate corrections on vacuum polarization for the 1¹S-state of the two-electron ³He and ⁴He atoms and for some two-electron ions, including the Li⁺, Be²⁺, B³⁺ and C⁴⁺ ions. The correction for vacuum polarization in two-electron He atoms has been evaluated as ΔE_ueh ≈ 7.253 ± 0.0025 × 10⁻⁷ a.u. The analogous corrections in the two-electron He-like ions rapidly increase with the nuclear charge Q (ΔE_ueh ≈ 2.7061 × 10⁻⁶ a.u. for the Li⁺ ion and ΔE_ueh ≈ 2.3495 × 10⁻⁵ a.u. for the C⁴⁺ ion). The corresponding corrections have also been evaluated for the electron–nucleus and electron–electron interactions.     

Integrals,Quantum Galois Extensions,and the Affineness Criterion for Quantum Yetter–Drinfel'd Modules

In this paper we shall generalize the notion of an integral on a Hopf algebra introduced by Sweedler, by defining the more general concept of an integral of a threetuple (H, A, C), where H is a Hopf algebra coacting on an algebra A and acting on a coalgebra C. We prove that there exists a total integral γ: C → Hom(C, A) of (H, A, C) if and only if any representation of (H, A, C) is injective in a functorial way, as a corepresentation of C. In particular, the quantum integrals associated to Yetter–Drinfel'd modules are defined. Let now A be an H-bicomodule algebra, ^H$\text{Y}$$\text{D}$_A the category of quantum Yetter–Drinfel'd modules, and B = {a ∈ A|∑S^− 1(a_〈1〉)a_{〈 − 1〉} ⊗ a_〈0〉 = 1_H ⊗ a}, the subalgebra of coinvariants of the Verma structure A ∈ ^H$\text{Y}$$\text{D}$_A. We shall prove the following affineness criterion: if there exists γ: H → Hom(H, A) a total quantum integral and the canonical map β: A ⊗ _B A → H ⊗ A, β(a ⊗ _B b) = ∑S^− 1(b_〈1〉)b_{〈 − 1〉} ⊗ ab_〈0〉 is surjective (i.e., A/B is a quantum homogeneous space), then the induction functor – ⊗ _B A: $\text{M}$_B → ^H$\text{Y}$$\text{D}$_A is an equivalence of categories. The affineness criteria proven by Cline, Parshall, and Scott, and independently by Oberst (for affine algebraic groups schemes) and Schneider (in the noncommutative case), are recovered as special cases.     

Magnetohydrodynamic stationary and oscillatory convective stability in a mushy layer during binary alloy solidification

For the case of solidification of a bottom cooled binary alloy, the magnetohydrodynamic stationary and oscillatory convective stability in the mushy layer is investigated analytically using normal mode linear stability analysis. In the limit of large Stefan number (S_t), a near–eutectic approximation with large far field temperature is considered in the present research. To ascertain the instability in the mushy layer, the strength of the superimposed magnetic field is so chosen that it corresponds to a given mush Hartmann number (Ha_m) of the problem. The results are presented for various values of mush Hartmann numbers in the range, 0 ≤ Ha_m ≤ 50. The critical Rayleigh number for stationary convection shows a linear relationship with increasing Ha_m. The magnetohydrodynamic effect imparts a stabilizing influence during stationary convection. In comparison to that of the stationary convective mode, the oscillatory mode appears to be critically susceptible at higher values of β (β = S_t/℘² ϒ², ℘ is the compositional ratio, ϒ = 1 + S_t/℘), and vice versa for lower β values. Analogous to the behavior for stationary convection, the magnetic field also offers a stabilizing effect in oscillatory convection and thus influences global stability of the mushy layer. Increasing magnetic strength shows reduction in the wavenumber and in the number of rolls formed in the mushy layer.     

Free convection heat and mass transfer in a power law fluid past an inclined surface with thermophoresis

The problem of heat and mass transfer in a power law, two-dimensional, laminar, boundary layer flow of a viscous incompressible fluid over an inclined plate with heat generation and thermophoresis is investigated by the characteristic function method. The governing non-linear partial differential equations describing the flow and heat transfer problem are transformed into a set of coupled non-linear ordinary differential equation which was solved using Runge–Kutta shooting method. Exact solutions for the dimensionless temperature and concentration profiles, are presented graphically for different physical parameters and for the different power law exponents 0 < n < 0.5 and for n > 0.5.     

Computational modelling of fluid structure interaction at nano-scale boundaries with modified Maxwellian velocity distribution

The soft collisions among fluid–fluid and fluid-wall molecules are modeled from first principles. In particular, the assumption of Maxwellian distribution of velocities for thermalized molecules, in both parallel and perpendicular directions to the wall, has been re-evaluated with supporting experimental and/or numerical evidence.It is proposed that the normal component of molecular velocity post collision is conserved for all fluid molecules. The slip effect at the wall boundary, introduced by the surface roughness, is accounted by an accommodation coefficient f. A moving least square method is used to calculate macroscopic velocity values. The influence of molecular interaction on the macroscopic velocity distribution is investigated at 40 MPa and 300 K for slit pore, inclined and stepped wall configurations. The accommodation coefficient values f = 0, 0.07, 0.257, 0.45, 0.681 and 1; and acceleration values ranging from zero to 1 × 10¹¹ m/s² and 250 × 10¹¹ m/s² are used for comparison.The distribution of macroscopic velocity parallel to the wall is studied to observe the effect of the slip behaviour. The detailed study of average of velocity values at various magnitudes of acceleration has shown an evidence of characteristic low and high speed of molecular flows that is considered as significant and a comparison is sought with an equivalent laminar and turbulent flow style behaviour. The two dimensional vector and contour plots of macroscopic velocity provide further insights in understanding Continuum velocity distributions resulting from molecular fluid-wall interaction at nanoscale. The research has highlighted the need to develop molecular dynamics simulation techniques for non-periodic boundary conditions.     

New exact solitary-wave special solutions for the nonlinear dispersive K(m, n) equations

In this paper, we study the nonlinear dispersive K(m, n) equations: u_t + (u^m)_x − (uⁿ)_xxx = 0 which exhibit solutions with solitary patterns. New exact solitary solutions are found. The two special cases, K(2, 2) and K(3, 3), are chosen to illustrate the concrete features of the decomposition method in K(m, n) equations. The nonlinear equations K(m, n) are studied for two different cases, namely when m = n being odd and even integers. General formulas for the solutions of K(m, n) equations are established.     

Statistical properties of nucleotides in human chromosomes 21 and 22

In this paper the statistical properties of nucleotides in human chromosomes 21 and 22 are investigated. The n-tuple Zipf analysis with n = 3, 4, 5, 6, and 7 is used in our investigation. It is found that the most common n-tuples are those which consist only of adenine (A) and thymine (T), and the rarest n-tuples are those in which GC or CG pattern appears twice. With the n-tuples become more and more frequent, the double GC or CG pattern becomes a single GC or CG pattern. The percentage of four nucleotides in the rarest ten and the most common ten n-tuples are also considered in human chromosomes 21 and 22, and different behaviors are found in the percentage of four nucleotides. Frequency of appearance of n-tuple f(r) as a function of rank r is also examined. We find the n-tuple Zipf plot shows a power-law behavior for r < 4ⁿ⁻¹ and a rapid decrease for r > 4ⁿ⁻¹. In order to explore the interior statistical properties of human chromosomes 21 and 22 in detail, we divide the chromosome sequence into some moving windows and we discuss the percentage of ξη (ξ, η = A, C, G, T) pair in those moving windows. In some particular regions, there are some obvious changes in the percentage of ξη pair, and there maybe exist functional differences. The normalized number of repeats N₀(l) can be described by a power law: N₀(l) ∼ l^−μ. The distance distributions P₀(S) between two nucleotides in human chromosomes 21 and 22 are also discussed. A two-order polynomial fit exists in those distance distributions: log P₀(S) = a + bS + cS², and it is quite different from the random sequence.     

Special Functions on the Sphere with Applications to Minimal Surfaces

A function which is homogeneous in x, y, z of degree n and satisfies V_xx + V_yy + V_zz = 0 is called a spherical harmonic. In polar coordinates, the spherical harmonics take the form rⁿf_n, where f_n is a spherical surface harmonic of degree n. On a sphere, f_n satisfies ▵ f_n + n(n + 1)f_n = 0, where ▵ is the spherical Laplacian. Bounded spherical surface harmonics are well studied, but in certain instances, unbounded spherical surface harmonics may be of interest. For example, if X is a parameterization of a minimal surface and n is the corresponding unit normal, it is known that the support function, w = X · n, satisfies ▵w + 2w = 0 on a branched covering of a sphere with some points removed. While simple in form, the boundary value problem for the support function has a very rich solution set. We illustrate this by using spherical harmonics of degree one to construct a number of classical genus-zero minimal surfaces such as the catenoid, the helicoid, Enneper's surface, and Hennenberg's surface, and Riemann's family of singly periodic genus-one minimal surfaces.