Asymmetric heat conduction through a weak link

We study the heat conduction of two nonlinear lattices joined by a weak harmonic link. When the system reaches a steady state, the heat conduction of the system is decided by the tunneling heat flow through the weak link. We present an analytical analysis by the combination of the self-consistent phonon theory and the heat tunneling transport formalism, and then the tunneling heat flow can be obtained. Moreover, the nonequilibrium molecular dynamics simulations are performed and the simulations results are consistent with the analytical predictions.

Dimensional crossover of heat conduction in low dimensions

The dimensional crossover phenomena of heat conduction is studied by a two-dimensional (2D) Fermi-Pasta-Ulam lattice. The 2D divergence law of the thermal conductivity is confirmed by the simulations results. The divergence law of the thermal conductivity will change from the 2D class to 1D class as delta=N-y/N-x decreases, here N-y is the size in transverse direction and N-x in longitude direction. The simulation's results suggest that the dimensional crossover happens in delta(*)-> 0 as N-x ->infinity.

Asymmetric heat conduction in nonlinear lattices

In this Letter, we conduct an extensive study of the two-segment Frenkel-Kontorova model. We show that the rectification effect of the heat flux reported in recent literature is possible only in the weak interfacial coupling limit. The rectification effect will be reversed when the properties of the interface and the system size change. These two types of asymmetric heat conduction are governed by different mechanisms though both are induced by nonlinearity. An intuitive physical picture is proposed to interpret the reversal of the rectification effect. Since asymmetric heat conduction depends critically on the properties of the interface and the system size, it is probably not an easy task to fabricate a thermal rectifier or thermal diode in practice.

Performance evaluation of a distributed algorithm for an inverse heat conduction problem

Numerical solutions of realistic 2-D and 3-D inverse problems may require a very large amount of computation. A two-level concept on parallelism is often used to solve such problems. The primary level uses the problem partitioning concept which is a decomposition based on the mathematical/physical problem. The secondary level utilizes the widely used data partitioning concept. A theoretical performance model is built based on the two-level parallelism. The observed performance results obtained from a network of general purpose Sun Sparc stations are compared with the theoretical values. Restrictions of the theoretical model are also discussed.

Parallel genetic algorithms for the solution of inverse heat conduction problems

A parallel genetic algorithm (PGA) is proposed for the solution of two-dimensional inverse heat conduction problems involving unknown thermophysical material properties. Experimental results show that the proposed PGA is a feasible and effective optimization tool for inverse heat conduction problems

Transient Heat Conduction of Variable Stiffness Composite Laminate

Recent research on Variable Stiffness (VS) laminates, which are constructed by steering the fiber orientation as a spatial function of location, have shown to improve laminate performance under mechanical loads. Two distinct cases of stiffness variation can be achieved either by variation of the fiber orientation in the direction of the global x-axis, or perpendicular to it. In the present paper, thermal analysis of a VS laminate is performed to study the effect of steering fibers on transient heat conduction under uniform heat flux using finite element method. The goal of the present paper is a parametric study of the effect of variable stiffness properties on transient response including time to reach steady state and temperature profile. Also, stress resultants and maximum stress location are investigated under different boundary conditions. A FEM algorithm is applied to exactly incorporate the boundary conditions for stress resultant analysis.

Transient Heat Conduction of Variable Stiffness Composite Laminate

Recent research on Variable Stiffness (VS) laminates, which are constructed by steering the fiber orientation as a spatial function of location, have shown to improve laminate performance under mechanical loads. Two distinct cases of stiffness variation can be achieved either by variation of the fiber orientation in the direction of the global x-axis, or perpendicular to it. In the present paper, thermal analysis of VS laminate is performed to study the effect of steering fibers on transient heat conduction under uniform heat flux using finite element method. The goal of the present paper is a parametric study of the
effect of variable stiffness properties on transient response including time to reach steady state and temperature profile. Also, stress resultants and maximum stress location are investigated under different boundary conditions. A FEM algorithm is applied to exactly incorporate the boundary conditions.

Anisotropic adaptive resolution of boundary layers for heat conduction problems

We deal with the numerical solution of heat conduction problems featuring steep gradients. In order to solve the associated partial differential equation a finite volume technique is used and unstructured grids are employed. A discrete maximum principle for triangulations of a Delaunay type is developed. To capture thin boundary layers incorporating steep gradients an anisotropic mesh adaptation technique is implemented. Computational tests are performed for an academic problem where the exact solution is known as well as for a real world problem of a computer simulation of the thermoregulation of premature infants.

Heat Conduction Equation Solution in the Presence of a Change of State in a Bounded Axisymmetric Cylindrical Domain

The heat conduction problem, in the presence of a change of state, was solved for the case of an indefinitely long cylindrical layer cavity. As boundary conditions, it is imposed that the internal surface of the cavity is maintained below the fusion temperature of the infilling substance and the external surface is kept above it. The solution, obtained in nondimensional variables, consists in two closed form heat conduction equation solutions for the solidified and liquid regions, which formally depend of the, at first, unknown position of the phase change front. The energy balance through the phase change front furnishes the equation for time dependence of the front position, which is numerically solved. Substitution of the front position for a particular instant in the heat conduction equation solutions gives the temperature distribution inside the cavity at that moment. The solution is illustrated with numerical examples. [DOI: 10.1115/1.4003542]

Non-newtonian heat transfer on a plate heat exchanger with generalized configurations

For the configuration optimization of plate heat exchangers (PHEs), the mathematical models for heat transfer and pressure drop must be valid for a wide range of operational conditions of all configurations of the exchanger or the design results may be compromised. In this investigation, the thermal model of a PHE is adjusted to fit experimental data obtained from non-Newtonian heat transfer for eight different configurations, using carboxymethylcellulose solutions (CMC) as test fluid. Although it is possible to successfully adjust the model parameters, Newtonian and non-Newtonian heat transfer cannot be represented by a single generalized correlation. In addition, the specific heat, thermal conductivity and power-law rheological parameters of CMC solutions were correlated with temperature, over a range compatible with a continuous pasteurization process.

Transient heat conduction problems using B.I.E.M

A method, using boundary elements, is presented as a solution to plane transient heat conduction. The proposed method considers the governing equation to be a Helmholtz's equation and solves the problem of time variation using step by step integration. A numerical procedure is developed and its effectiveness verified. Several examples are provided and their results compared with the theoretical ones.

FTCS finite difference scheme GPGPU parallel computing for the heat conduction equation = Programación en paralelo GPGPU del método en diferencias finitas FTCS para la ecuación del calor

En el presente artículo se muestran las ventajas de la programación en paralelo resolviendo numéricamente la ecuación del calor en dos dimensiones a través del método de diferencias finitas explícito centrado en el espacio FTCS. De las conclusiones de este trabajo se pone de manifiesto la importancia de la programación en paralelo para tratar problemas grandes, en los que se requiere un elevado número de cálculos, para los cuales la programación secuencial resulta impracticable por el elevado tiempo de ejecución. En la primera sección se describe brevemente los conceptos básicos de programación en paralelo. Seguidamente se resume el método de diferencias finitas explícito centrado en el espacio FTCS aplicado a la ecuación parabólica del calor. Seguidamente se describe el problema de condiciones de contorno y valores iniciales específico al que se va a aplicar el método de diferencias finitas FTCS, proporcionando pseudocódigos de una implementación secuencial y dos implementaciones en paralelo. Finalmente tras la discusión de los resultados se presentan algunas conclusiones. In this paper the advantages of parallel computing are shown by solving the heat conduction equation in two dimensions with the forward in time central in space (FTCS) finite difference method. Two different levels of parallelization are consider and compared with traditional serial procedures. We show in this work the importance of parallel computing when dealing with large problems that are impractical or impossible to solve them with a serial computing procedure. In the first section a summary of parallel computing approach is presented. Subsequently, the forward in time central in space (FTCS) finite difference method for the heat conduction equation is outline, describing how the heat flow equation is derived in two dimensions and the particularities of the finite difference numerical technique considered. Then, a specific initial boundary value problem is solved by the FTCS finite difference method and serial and parallel pseudo codes are provided. Finally after results are discussed some conclusions are presented.