Spatial behaviour of solutions of the dual‐phase‐lag heat equation

In this paper we study the spatial behaviour of solutions of some problems for the dual‐phase‐lag heat equation on a semi‐infinite cylinder. The theory of dual‐phase‐lag heat conduction leads to a hyperbolic partial differential equation with a third derivative with respect to time. First, we investigate the spatial evolution of solutions of an initial boundary‐value problem with zero boundary conditions on the lateral surface of the cylinder. Under a boundedness restriction on the initial data, an energy estimate is obtained. An upper bound for the amplitude term in this estimate in terms of the initial and boundary data is also established. For the case of zero initial conditions, a more explicit estimate is obtained which shows that solutions decay exponentially along certain spatial‐time lines. A class of non‐standard problems is also considered for which the temperature and its first two time derivatives at a fixed time T are assumed proportional to their initial values. Copyright © 2004 John Wiley & Sons, Ltd.     

The dynamics ofn weakly coupled identical oscillators

Summary We present a framework for analysing arbitrary networks of identical dissipative oscillators assuming weak coupling. Using the symmetry of the network, we find dynamically invariant regions in the phase space existing purely by virtue of their spatio-temporal symmetry (the temporal symmetry corresponds to phase shifts). We focus on arrays which are symmetric under all permutations of the oscillators (this arises with global coupling) and also on rings of oscillators with both directed and bidirectional coupling. For these examples, we classify all spatio-temporal symmetries, including limit cycle solutions such as in-phase oscillation and those involving phase shifts. We also show the existence of “submaximal” limit cycle solutions under generic conditions. The canonical invariant region of the phase space is defined and used to investigate the dynamics. We discuss how the limit cycles lose and gain stability, and how symmetry can give rise to structurally stable heteroclinic cycles, a phenomenon not generically found in systems without symmetry. We also investigate how certain types of coupling (including linear coupling between oscillators with symmetric waveforms) can give rise to degenerate behaviour, where the oscillators decouple into smaller groups.     

Similarity solutions for solute transport in fractal porous media using a time- and scale-dependent dispersivity

A specific form of the Fokker–Planck equation with a time- and scale-dependent dispersivity is presented for modelling solute transport in saturated heterogeneous porous media. By taking a dispersivity in the form of separable power-law dependence on both time and scale, we are able to show the existence of similarity solutions. Explicit closed-form solutions are then derived for an instantaneous point-source (Dirac delta function) input, and for constant concentration and constant flux boundary conditions on a semi-infinite domain. The solutions have realistic behaviour when compared to tracer breakthrough curves observed under both field and laboratory conditions. Direct comparison with the experimental laboratory data of Pang and Hunt [J. Contam. Hydrol. 53 (2001) 21] shows good agreement between the source solutions and the measured breakthrough curves.     

Modelling and remodelling of biological tissue in the framework of continuum biomechanics

A biological tissue in general is formed by cells, extracellular matrix (ECM) and fluids. Consequently, its overall material behaviour results from its components and their interaction among each other. Furthermore, in case of living tissues, the material properties do not remain constant but naturally change due to adaptation processes or diseases. In the context of the Theory of Porous Media (TPM), a continuum-mechanical model is introduced to describe the complex fluid-structure interaction in biological tissue on a macroscopic scale. The tissue is treated as an aggregate of two immiscible constituents, where the cells and the ECM are summarised to a solid phase, whereas the fluid phase represents the extracellular and interstitial liquids as well as necrotic debris and cell or matrix precursors in solution. The growth and remodelling processes are described by a distinct mass exchange between the fluid and solid phase, which also results in a change of the constituent material behaviour. To furthermore guarantee the compliance with the entropy principle, the growth energy is introduced as an additional quantity. It measures the average of chemical energy available for cell metabolism, and thus, controls the growth and remodelling processes. To set an example, the presented model is applied for the simulation of the early stages of avascular tumour growth in the framework of the finite element method (FEM). (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Analytical and numerical solutions to some one-dimensional relativistic heat equations

Analytical and numerical solutions to a family of one-dimensional (nonlinear) relativistic heat equations in finite domains are presented. The analytical solutions correspond to steady state conditions in the absence of source terms and have been obtained as functions of the (absolute) temperature at one of the boundaries, a characteristic exponent, a Péclet number based on the speed of light and the heat flux, while the numerical ones correspond to an initial Gaussian temperature distribution, adiabatic boundary conditions and different values of the Péclet number and a characteristic exponent. It is shown that, for steady conditions, the difference between the (nondimensional) temperature of the relativistic heat equation and that corresponding to Fourier law is very large for large values of both the coefficient and the exponent of the nonlinearity that characterize the relativistic contribution to the heat flux, small values of the temperature at one of the boundaries and large heat fluxes. Travelling-wave solutions of the wave-front type are reported for odd values of the nonlinearity exponent in infinite domains and in the absence of source terms. For an initial Gaussian distribution, it is shown that the relativistic contribution to heat transfer results in the formation of two triangular corner regions where the temperature is equal to the initial one, and the formation of two temperature fronts that propagate towards the domain’s boundaries. The amplitude and steepness of these fronts increase whereas their width and speed decrease as the Péclet number is decreased. It is also shown that the effects of the characteristic exponent are small provided that its value is greater than about two, and that, in the absence of source terms, the temperature becomes uniform in space and constant in time for adiabatic boundary conditions. In the presence of source terms and for adiabatic boundary conditions, it is shown that, soon after the temperature fronts hit the boundaries, the temperature becomes uniform in space but may either increase or decrease with time until it reaches a stable fixed point of the source term. For a cubic source term that exhibits bistability, it is shown that the temperature tends to the attractor of lowest temperature.     

A non-isentropic Euler–Poisson model for a collisionless plasma

We analyse transonic solutions of the one-dimensional Euler–Poisson model for a collisionless gas of charged particles in the non-isentropic steady-state case. The model consists of the conservation of mass, momentum and energy equations. The electric field is modelled self-consistently (Coulomb field). Boundary conditions on the particle density and particle temperature are imposed. The analysis is based on representing solutions piecewise as orbits in the particle-density-electric-field phase plane and connecting the orbit segments by the jump and entropy conditions. We characterize the set of all solutions of the Euler–Poisson problem. In particular, we show that, depending upon the length of the interval on which the boundary value problem is posed, fully subsonic, one-shock and (in certain cases) two-shock transonic and smooth transonic solutions exist. Also, numerical computations illustrating the structure of the solutions are reported.     

Spatial behaviour of solutions of the three-phase-lag heat equation

In this paper we study the spatial behaviour of solutions for the three-phase-lag heat equation on a semi-infinite cylinder. The theory of three-phase-lag heat conduction leads to a hyperbolic partial differential equation with a fourth-order derivative with respect to time. First, we investigate the spatial evolution of solutions of an initial boundary-value problem with zero boundary conditions on the lateral surface of the cylinder. Under a boundedness restriction on the initial data, an energy estimate is obtained. An upper bound for the amplitude term in this estimate in terms of the initial and boundary data is also established. For the case of zero initial conditions, a more explicit estimate is obtained which shows that solutions decay exponentially along certain spatial-time lines. A class of non-standard problems is also considered for which the temperature and its first two time derivatives at a fixed time T₀ are assumed proportional to their initial values. These results are relevant in the context of the Saint-Venant Principle for heat conduction problems.     

Numerical study of melting in an enclosure with discrete protruding heat sources

In order to explore the capability of a solid–liquid phase change material (PCM) for cooling electronic or heat storage applications, melting of a PCM in a vertical rectangular enclosure was studied. Three protruding generating heat sources are attached on one of the vertical walls of the enclosure, and generating heat at a constant and uniform volumetric rate. The horizontal walls are adiabatic. The power generated in heat sources is dissipated in PCM (n-eicosane with the melting temperature, T_m = 36 °C) that filled the rectangular enclosure. The advantage of using PCM is that it is able to absorb high amount of heat generated by heat sources due to its relatively high energy density. To investigate the thermal behaviour and thermal performance of the proposed system, a mathematical model based on the mass, momentum and energy conservation equations was developed. The governing equations are next discretised using a control volume approach in a staggered mesh and a pressure correction equation method is employed for the pressure–velocity coupling. The PCM energy equation is solved using the enthalpy method. The solid regions (wall and heat sources) are treated as fluid regions with infinite viscosity and the thermal coupling between solid and fluid regions is taken into account using the harmonic mean of the thermal conductivity method. The dimensionless independent parameters that govern the thermal behaviour of the system were next identified. After validating the proposed mathematical model against experimental data, a numerical investigation was next conducted in order to examine the thermal behaviour of the system by analyzing the flow structure and the heat transfer during the melting process, for a given values of governing parameters.     

Properties of nonadiabatic combustion waves in competitive exothermic reactions

We study the properties of travelling combustion waves in a diffusional thermal model with a two-step competitive exothermic reaction mechanism considering the system in one spatial dimension under nonadiabatic conditions. Based on the notion of the crossover temperature, the model is first examined analytically to predict the behaviour of travelling combustion waves in the limit of large activation energies. It is then studied numerically over a wide range of parameter values, such as those describing the ratios of the enthalpies, pre-exponential factors and activation energies. It is demonstrated that the nonadiabatic flame speed as a function of these parameters is a single-valued monotonic function with two flame regimes identified, each of which represents the region of the parameter values when one reaction dominates the other, while the flame speed as a function of the activation energy of the reaction R2 (the exothermicity parameter) is either a c- or m-shaped. The extinction conditions for these two flame regimes are investigated analytically.     

On resonant interactions of gravity‐capillary waves without energy exchange

We consider resonant triad interactions of gravity‐capillary waves and investigate in detail special resonant triads that exchange no energy during their interactions so that the wave amplitudes remain constant in time. After writing the resonance conditions in terms of two parameters (or two angles of wave propagation), we first identify a region in the two‐dimensional parameter space, where resonant triads can be always found, and then describe the variations of resonant wavenumbers and wave frequencies over the resonance region. Using the amplitude equations recovered from a Hamiltonian formulation for water waves, it is shown that any resonant triad inside the resonance region can interact without energy exchange if the initial wave amplitudes and relative phase satisfy the two conditions for fixed point solutions of the amplitude equations. Furthermore, it is shown that the symmetric resonant triad exchanging no energy forms a transversely modulated traveling wave field, which can be considered a two‐dimensional generalization of Wilton ripples.     

Numerical simulation of powder mixed electric discharge machining (PMEDM) using finite element method

In the present paper, an axisymmetric two-dimensional model for powder mixed electric discharge machining (PMEDM) has been developed using the finite element method (FEM). The model utilizes the several important aspects such as temperature-sensitive material properties, shape and size of heat source (Gaussian heat distribution), percentage distribution of heat among tool, workpiece and dielectric fluid, pulse on/off time, material ejection efficiency and phase change (enthalpy) etc. to predict the thermal behaviour and material removal mechanism in PMEDM process. The developed model first calculates the temperature distribution in the workpiece material using ANSYS (version 5.4) software and then material removal rate (MRR) is estimated from the temperature profiles. The effect of various process parameters on temperature distributions along the radius and depth of the workpiece has been reported. Finally, the model has been validated by comparing the theoretical MRR with the experimental one obtained from a newly designed experimental setup developed in the laboratory.     

Two-parameter extremum problems of boundary control for stationary thermal convection equations

Two-parameter extremum problems of boundary control are formulated for the stationary thermal convection equations with Dirichlet boundary conditions for velocity and with mixed boundary conditions for temperature. The cost functional is defined as the root mean square integral deviation of the desired velocity (vorticity, or pressure) field from one given in some part of the flow region. Controls are the boundary functions involved in the Dirichlet condition for velocity on the boundary of the flow region and in the Neumann condition for temperature on part of the boundary. The uniqueness of the extremum problems is analyzed, and the stability of solutions with respect to certain perturbations in the cost functional and one of the functional parameters of the original model is estimated. Numerical results for a control problem associated with the minimization of the vorticity norm aimed at drag reduction are discussed.     

运动热源产生的温度场的有限元分析

本文采用运动元法,考虑了热辐射、热交换和变热传导系数,在较宽的速度范围内研究了运动热源所产生的温度场.给出了在运动坐标和静止坐标系中,不同速度情况下温度场随时间变化的情形,给出了在运动坐标系中不同速度情况下温度场的(定常)等温线分布图.本文还讨论了动态裂纹扩展中裂纹端部过程区的塑性变形所引起的温度场,结果表明,结构钢中过程区的温度一般不会超过1000℃,或1832℉.     

Global existence and blow-up of solutions for a system of Petrovsky equations

The initial-boundary value problem for a system of Petrovsky equations with memory and nonlinear source terms in bounded domain is studied. The existence of global solutions for this problem is proved by constructing a stable set, and obtain the exponential decay estimate of global solutions. Meanwhile, under suitable conditions on relaxation functions and the positive initial energy as well as non-positive initial energy, it is proved that the solutions blow up in the finite time and the lifespan estimates of solutions are also given.     

A phase field approach to solidification and solute separation in water solutions

We propose a phase field model for the solid–liquid phase transition in a water-salt (sodium chloride) solution in the absence of macroscopic motion, under possibly non-isothermal conditions. A thermodynamic approach based on a free energy functional is assumed. The model consists of three evolution equations: a time-dependent Ginzburg–Landau equation for the solid–liquid phase change, a diffusion equation of the Cahn–Hilliard kind for the solute dynamics and the heat equation for the temperature change. The proposed system is aimed to contribute to the modelling of the brine channels formation in the ice of the polar seas.     

Studies of a mathematical model for temperature-modulated bioluminescence tomography

We introduce and study a mathematical model for temperature-modulated bioluminescence tomography (TBT). The model is capable of self-adjusting values of experimental parameters that are used in the formulation. Major theoretical results of this article include: Solution existence of the model, convergence of numerical solutions, an iterative scheme based on linearization, studies of the solution limiting behaviours when normalized total energy function and/or some or all the energy percentages in individual spectral bands are known exactly. Several numerical examples are included to illustrate the improvement of the accuracy of the reconstructed bioluminescent source distribution due to the employment of measurements from multiple temperature distributions.     

Gas injection into a moist porous medium with the formation of a gas hydrate

The model problem of the formation of a gas hydrate when a gas is injected into a porous medium, filled in the initial state with a gas and water, is considered in the one-dimensional approximation. A detailed pattern of the seepage flow with phase transitions for different modes of gas injection is obtained. Three seepage modes in a porous medium are possible, which differ qualitatively in the temperature and hydrate saturation fields. At low boundary pressures no hydrate is formed and the temperature distribution increases monotonically. As the boundary pressure increases, when the corresponding values of the pressure and temperature on the phase diagram lie in the region of gas-hydrate stability (below the equilibrium curve), a purely frontal pattern of hydrate formation is obtained with a monotonic temperature distribution. When the boundary pressure is increased further, an extended region of hydrate formation appears with a convex temperature profile, where, depending on the values of the boundary pressure, the hydrate saturation may be continuous (at high boundary pressures) or change abruptly at lower boundary pressures.     

Wärmewellen in einem leitenden und strahlenden Medium

Summary Plane thermal waves in a heat conducting and radiating (emitting and absorbing) medium that occupies the half-spacex>0 are investigated. The governing equations for a gray medium are linearized with regard to small perturbations of the radiative equilibrium. Solutions are given for the thermal wave that is due to harmonic oscillations of either the wall temperature or the radiative energy flux produced by an outer source. The behaviour of the thermal wave is then discussed for the asymptotic cases of weak, strong, optical thin and optical thick radiation, respectively, and also for the special case that the Bouguer numberBu and the radiation-conduction parameterK as defined in the text are equal to one. Then the equations and their solutions are generalized in order to apply to certain models of frequency-dependent absorption coefficients (non-gray media). Finally it is shown that nonlinear terms, although being small of higher order in the differential equations, cause the perturbation solution to be not uniformly valid as the distance from the boundary surface goes to infinity.     

Polynomial decay and blow up of solutions for variable coefficients viscoelastic wave equation with acoustic boundary conditions

A variable coefficient viscoelastic wave equation with acoustic boundary conditions and nonlinear source term is considered. Under suitable conditions on the initial data and the relaxation function g, we show the polynomial decay of the energy solution and the blow up of solutions by energy methods. The estimates for the lifespan of solutions are also given.