A compact difference scheme for numerical solutions of second order dual-phase-lagging models of microscale heat transfer

Dual-phase-lagging (DPL) models constitute a family of non-Fourier models of heat conduction that allow for the presence of time lags in the heat flux and the temperature gradient. These lags may need to be considered when modeling microscale heat transfer, and thus DPL models have found application in the last years in a wide range of theoretical and technical heat transfer problems. Consequently, analytical solutions and methods for computing numerical approximations have been proposed for particular DPL models in different settings. In this work, a compact difference scheme for second order DPL models is developed, providing higher order precision than a previously proposed method. The scheme is shown to be unconditionally stable and convergent, and its accuracy is illustrated with numerical examples.

A second order control-volume finite-element least-squares strategy for simulating diffusion in strongly anisotropic media

An unstructured mesh �nite volume discretisation method for simulating di�usion in anisotropic media in two-dimensional space is discussed. This technique is considered as an extension of the fully implicit hybrid control-volume �nite-element method and it retains the local continuity of the ux at the control volume faces. A least squares function recon- struction technique together with a new ux decomposition strategy is used to obtain an accurate ux approximation at the control volume face, ensuring that the overall accuracy of the spatial discretisation maintains second order. This paper highlights that the new technique coincides with the traditional shape function technique when the correction term is neglected and that it signi�cantly increases the accuracy of the previous linear scheme on coarse meshes when applied to media that exhibit very strong to extreme anisotropy ratios. It is concluded that the method can be used on both regular and irregular meshes, and appears independent of the mesh quality.

Numerical methods for second-order stochastic differential equations

We seek numerical methods for second‐order stochastic differential equations that reproduce the stationary density accurately for all values of damping. A complete analysis is possible for scalar linear second‐order equations (damped harmonic oscillators with additive noise), where the statistics are Gaussian and can be calculated exactly in the continuous‐time and discrete‐time cases. A matrix equation is given for the stationary variances and correlation for methods using one Gaussian random variable per timestep. The only Runge–Kutta method with a nonsingular tableau matrix that gives the exact steady state density for all values of damping is the implicit midpoint rule. Numerical experiments, comparing the implicit midpoint rule with Heun and leapfrog methods on nonlinear equations with additive or multiplicative noise, produce behavior similar to the linear case.

Numerical analysis for a variable-order nonlinear cable equation

In this paper, a variable-order nonlinear cable equation is considered. A numerical method with first-order temporal accuracy and fourth-order spatial accuracy is proposed. The convergence and stability of the numerical method are analyzed by Fourier analysis. We also propose an improved numerical method with second-order temporal accuracy and fourth-order spatial accuracy. Finally, the results of a numerical example support the theoretical analysis.

Analysis of steel structures for first- and second-order element (displacement and force) solutions

Finite element method (FEM) relies on an approximate function to fit into a governing equation and minimizes the residual error in the integral sense in order to generate solutions for the boundary value problems (nodal solutions). Because of this FEM does not show simultaneous capacities for accurate displacement and force solutions at node and along an element, especially when under the element loads, which is of much ubiquity. If the displacement and force solutions are strictly confined to an element’s or member’s ends (nodal response), the structural safety along an element (member) is inevitably ignored, which can definitely hinder the design of a structure for both serviceability and ultimate limit states. Although the continuous element deflection and force solutions can be transformed into the discrete nodal solutions by mesh refinement of an element (member), this setback can also hinder the effective and efficient structural assessment as well as the whole-domain accuracy for structural safety of a structure. To this end, this paper presents an effective, robust, applicable and innovative approach to generate accurate nodal and element solutions in both fields of displacement and force, in which the salient and unique features embodies its versatility in applications for the structures to account for the accurate linear and second-order elastic displacement and force solutions along an element continuously as well as at its nodes. The significance of this paper is on shifting the nodal responses (robust global system analysis) into both nodal and element responses (sophisticated element formulation).

Numerical analysis of a new space-time variable fractional order advection-dispersion equation

Many physical processes appear to exhibit fractional order behavior that may vary with time and/or space. The continuum of order in the fractional calculus allows the order of the fractional operator to be considered as a variable. In this paper, we consider a new space–time variable fractional order advection–dispersion equation on a finite domain. The equation is obtained from the standard advection–dispersion equation by replacing the first-order time derivative by Coimbra’s variable fractional derivative of order α(x)∈(0,1]α(x)∈(0,1], and the first-order and second-order space derivatives by the Riemann–Liouville derivatives of order γ(x,t)∈(0,1]γ(x,t)∈(0,1] and β(x,t)∈(1,2]β(x,t)∈(1,2], respectively. We propose an implicit Euler approximation for the equation and investigate the stability and convergence of the approximation. Finally, numerical examples are provided to show that the implicit Euler approximation is computationally efficient.

Compressible second-order boundary layers for three-dimensional stagnation point flows with mass transfer

All the second-order boundary-layer effects have been studied for the steady laminar compressible 3-dimensional stagnation-point flows with variable properties and mass transfer for both saddle and nodal point regions. The governing equations have been solved numerically using an implicit finite-difference scheme. Results for the heat transfer and skin friction have been obtained for several values of the mass-transfer rate, wall temperature, and also for several values of parameters characterizing the nature of stagnation point and variable gas properties. The second-order effects on the heat transfer and skin friction at the wall are found to be significant and at large injection rates, they dominate over the results of the first-order boundary layer, but the effect of large suction is just the opposite. In general, the second-order effects are more pronounced in the saddle-point region than in the nodal-point region. The overall heat-transfer rate for the 3-dimensional flows is found to be more than that of the 2-dimensional flows.

Unsteady laminar incompressible second-order boundary-layer flow at a three-dimensional stagnation point

The unsteadely laminar incompressible second-order boundary-layer flow at the stagnation point of a three-dimensional body has been studied for both nodal and saddle point regions. The effects of mass transfer and Prandtl number have been taken into account. The equations governing the flow have been solved numerically using an implicit finite-difference scheme. It has been found that the parameter characterizing the unsteadiness in the velocity of the free stream, the nature of the stagnation point, the mass transfer and Prandtl number strongly affect the second-order skin friction and heat transfer. The overall skin friction becomes less due to second-order effects but the heat transfer has the opposite behaviour. For large injection, the second-order skin-friction and heat-transfer results prevail over the first-order boundary layer results whereas for the case of large suction the behaviour is just the opposite.

Unsteady compressible second-order boundary layers at the stagnation point of two-dimensional and axisymmetric bodies

The unsteady laminar compressible boundary-layer flow over two-dimensional and axisymmetric bodies at the stagnation point with mass transfer has been studied for all second-order boundary layer effects when the basic potential flow admits selfsimilarity. The solutions for the governing equations are obtained by using an implicit finite-difference scheme. Computations have been carried out for different values of the parameters characterizing the unsteadiness in the free stream velocity, wall temperature, mass transfer rate and variable gas properties. The results are found to be strongly affected by the unsteadiness in the free stream velocity. For large injection rates the second-orderboundary layer effects may prevail over the first-order boundary layer, but reverse is true for suction. The wall temperature and the variation of the density-viscosity product across the boundary layer appreciably change the skin-friction and heat-transfer rates due to second-order boundary-layer effects.

Semi-similar solutions of the unsteady compressible second-order boundary layer flow at the stagnation point

Semi-similar solutions of the unsteady compressible laminar boundary layer flow over two-dimensional and axisymmetric bodies at the stagnation point with mass transfer are studied for all the second-order boundary layer effects when the free stream velocity varies arbitrarily with time. The set of partial differential equations governing the unsteady compressible second-order boundary layers representing all the effects are derived for the first time. These partial differential equations are solved numerically using an implicit finite-difference scheme. The results are obtained for two particular unsteady free stream velocity distributions: (a) an accelerating stream and (b) a fluctuating stream. It is observed that the total skin friction and heat transfer are strongly affected by the surface mass transfer and wall temperature. However, their variation with time is significant only for large times. The second-order boundary layer effects are found to be more pronounced in the case of no mass transfer or injection as compared to that for suction. Résumé Des solutions semi-similaires d'écoulement variable compressible de couche limite sur des corps bi-dimensionnels thermique, sont étudiées pour tous les effets de couche limite du second ordre, lorsque la vitesse de l'écoulement libre varie arbitrairement avec le temps. Le systéme d'équations aux dérivées partielles représentant tous les effets est écrit pour la premiére fois. On le résout numériquement á l'aide d'un schéma implicite aux différences finies. Les résultats sont obtenus pour deux cas de vitesse variable d'écoulement libre: (a) un écoulement accéléré et (b) un écoulement fluctuant. On observe que le frottement pariétal total et le transfert de chaleur sont fortement affectés par le transfert de masse et la température pariétaux. Néanmoins, leur variation avec le temps est sensible seulement pour des grandes durées. Les effets sont trouvés plus prononcés dans le cas de l'absence du transfert de masse ou de l'injection par rapport au cas de l'aspiration.

Second-order boundary-layer effects for the unsteady laminar incompressible three-dimensional stagnation-point flow

All the second-order boundary-layer effects on the unsteady laminar incompressible flow at the stagnation-point of a three-dimensional body for both nodal and saddle point regions have been studied. It has been assumed that the free-stream velocity, wall temperature and mass transfer vary arbitrarily with time. The effect of the Prandtl number has been taken into account. The partial differential equations governing the flow have been derived for the first time and then solved numerically unsteady free-stream velocity distributions, the nature of the using an implicit finite-difference scheme. It is found that the stagnation point and the mass transfer strongly affect the skin friction and heat transfer whereas the effects of the Prandtl number and the variation of the wall temperature with time are only on the heat transfer. The skin friction due to the combined effects of first- and second-order boundary layers is less than the skin friction due to, the first-order boundary layers whereas the heat transfer has the opposite behaviour. Suction increases the skin friction and heat transfer but injection does the opposite

Solution of Time Dependent Joule Heat Equation for a Graphene Sheet Under Thomson Effect

We address a physics-based solution of joule heating phenomenon in a single-layer graphene (SLG) sheet under the presence of Thomson effect. We demonstrate that the temperature in an isotopically pure (containing only C-12) SLG sheet attains its saturation level quicker than when doped with its isotopes (C-13). From the solution of the joule heating equation, we find that the thermal time constant of the SLG sheet is in the order of tenths of a nanosecond for SLG dimensions of a few micrometers. These results have been formulated using the electron interactions with the inplane and flexural phonons to demonstrate a field-dependent Landauer transmission coefficient. We further develop an analytical model of the SLG specific heat using the quadratic (out of plane) phonon band structure over the room temperature. Additionally, we show that a cooling effect in the SLG sheet can be substantially enhanced with the addition of C-13. The methodologies as discussed in this paper can be put forward to analyze the graphene heat spreader theory.

Composite second-order performance improvement in optical fibre CATV transmission system using chirped fibre grating

Theoretically, we analyse the dispersion compensation characteristics of the chirped fibre grating (CFG) in an optical fibre cable television (CATV) system and obtain the analytic expression of the composite second-order (CSO) distortion using the time-domain form of the field envelope wave equation. The obtained result is in good agreement with the numerical simulation result. Experimentally, we verify the result by making use of the tunable characteristics of CFG to change the dispersion compensation amount and obtain an optimal CSO performance in a 125km fibre transmission link. Both the theoretical and experimental results show that the CSO performance can be improved by properly choosing the dispersion compensation amount for a certain fibre transmission link.

Semi-implicit second order finite volume jet stream computations

Semi-implicit, second order temporal and spatial finite volume computations of the flow in a differentially heated rotating annulus are presented. For the regime considered, three cyclones and anticyclones separated by a relatively fast moving jet of fluid or "jet stream" are predicted. Two second order methods are compared with, first order spatial predictions, and experimental measurements. Velocity vector plots are used to illustrate the predicted flow structure. Computations made using second order central differences are shown to agree best with experimental measurements, and to be stable for integrations over long time periods (> 1000s). No periodic smoothing is required to prevent divergence.