On formation of singularities in one-dimensional nonlinear thermoelasticity

Dedicated to Bernard D. Coleman on the occasion of his 60^th birthday     

Global smooth solutions and asymptotic stability in one-dimensional nonlinear thermoelasticity

The global existence of smooth solutions to the equations of nonlinear thermoelasticity is shown for a one-dimensional homogeneous reference configuration. Dirichlet boundary conditions are studied and the asymptotic behaviour of the solutions as t is described.Dedicated to Professor Rolf Leis and to Proffessor Mutsuhide Matsumura on the occasion of their sixtieth birthdays in 1991     

Numerical nonlinear inverse problem of determining wall heat flux

The inverse problem of determining time-variable surface heat flux in a plane wall, with constant or temperature dependent thermal properties, is numerically studied. Different kinds of incident heat flux, including rectangular waveform, are assumed. The solution is numerically solved as a function estimation problem, so that no a priori information for the functional waveforms of the unknown heat flux is needed. In all cases, a solution in the form of a piece-wise function is used to approach the incident flux. Transient temperature measurements at the boundary, from the solution of the direct problem, served as the simulated experimental data needed as input for the inverse analysis. Both direct and inverse heat conduction problems are solved using the network simulation method. The solution is obtained step-by-step by minimising the classical functional that compares the above input data with those obtained from the solution of the inverse problem. A straight line of variable slope and length is used for each one of the stretches of the desired solution. The influence of random error, number of functional terms and the effect of sensor location are studied. In all cases, the results closely agree with the solution.     

Analysis of dispersion relations of a coupled thermoelasticity problem with regard to heat flux relaxation

Dispersion relations for a coupled thermoelasticity problem including a hyperbolic heat conduction equation are derived, and their asymptotic analysis is performed. Dependences of the wave number and characteristics of the vibration damping rate on frequency are obtained and compared with similar diagrams in the classical model.     

Unsteady one-dimensional flow in a saturated porous medium with a volume heat source

The one-dimensional process of the heating of a saturated porous medium by a volume heat source as a result of the absorption of the energy of a high-frequency (frequency R~ 10¹–10³ MHz) electromagnetic wave is investigated. It is assumed that in the initial state the saturating (second) component is in the high-viscosity liquid or solid state. Under the action of the heat it is heated, melts, expands, becomes less viscous and under the pressure head created may flow relative to the stationary rock skeleton (first component). On the basis of the mathematical model proposed the basic laws of the process are analyzed and numerically investigated in the case of one-dimensional axisymmetric motion. It is shown that under actual conditions the dimensions of the thermal influence zone may be very considerable. Thus, by varying certain external factors it is possible to modify the dynamics of the process and the distributions of the temperature, pressure and phase velocity fields.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 115–124, July–August, 1991.     

Analytical solution of a double moving boundary problem for nonlinear flows in one-dimensional semi-infinite long porous media with low permeability简

Based on Huang＇s accurate tri-sectional nonlin- ear kinematic equation （1997）, a dimensionless simplified mathematical model for nonlinear flow in one-dimensional semi-infinite long porous media with low permeability is presented for the case of a constant flow rate on the inner boundary. This model contains double moving boundaries, including an internal moving boundary and an external mov- ing boundary, which are different from the classical Stefan problem in heat conduction： The velocity of the external moving boundary is proportional to the second derivative of the unknown pressure function with respect to the distance parameter on this boundary. Through a similarity transfor- mation, the nonlinear partial differential equation （PDE） sys- tem is transformed into a linear PDE system. Then an ana- lytical solution is obtained for the dimensionless simplified mathematical model. This solution can be used for strictly checking the validity of numerical methods in solving such nonlinear mathematical models for flows in low-permeable porous media for petroleum engineering applications. Finally, through plotted comparison curves from the exact an- alytical solution, the sensitive effects of three characteristic parameters are discussed. It is concluded that with a decrease in the dimensionless critical pressure gradient, the sensi- tive effects of the dimensionless variable on the dimension- less pressure distribution and dimensionless pressure gradi- ent distribution become more serious; with an increase in the dimensionless pseudo threshold pressure gradient, the sensi- tive effects of the dimensionless variable become more serious; the dimensionless threshold pressure gradient （TPG） has a great effect on the external moving boundary but has little effect on the internal moving boundary.     

Influence of boundary conditions on the solution of a hyperbolic thermoelasticity problem

We consider a series of problems with a short laser impact on a thin metal layer accounting various boundary conditions of the first and second kind. The behavior of the material is modeled by the hyperbolic thermoelasticity of Lord–Shulman type. We obtain analytical solutions of the problems in the semi-coupled formulation and numerical solutions in the coupled formulation. Numerical solutions are compared with the analytical ones. The analytical solutions of the semi-coupled problems and numerical solutions of the coupled problems show qualitative match. The solutions of hyperbolic thermoelasticity problems are compared with those obtained in the frame of the classical thermoelasticity. It was determined that the most prominent difference between the classical and hyperbolic solutions arises in the problem with fixed boundaries and constant temperature on them. The smallest differences were observed in the problem with unconstrained, thermally insulated edges. It was shown that a cooling zone is observed if the boundary conditions of the first kind are given for the temperature. Analytical expressions for the velocities of the quasiacoustic and quasithermal fronts as well as the critical value for the attenuation coefficient of the excitation impulse are verified numerically.     

Numerical solution of a pulsatile flow problem

Summary A numerical method, developed and used in previous work, is used to solve the Navier-Stokes equations for a two-dimensional incompressible flow in variable motion. A pulsatile flow through a pipe cross is studied to simulate blood flow in vascular aorto-renal junctions.

---

Sommario In questo studio si è utilizzato uno schema numerico, già studiato e utilizzato in precedenti lavori, per risolvere le equazioni bidimensionali di Navier-Stokes per un fluido viscoso incomprimibile in un problema di moto vario. Nelle ipotesi di poter descrivere il moto fasico del sangue con il moto di un fluido viscoso newtoniano di opportuna viscosità, si è studiato il moto pulsante di questo fluido in due condotti raccordati per simulare il moto del sangue nel raccordo aorta-renale durante un ciclo temporale completo.

---

---

Research carried out at Istituto di Idraulica del Politecnico di Milano, Milan, Italy.     

Strict solution of a nonlinear particle transport problem

Summary We show that an evolution problem arising from particle transport with quadratic nonlinearities has a unique positive strict solution in the Banach space UCB.

---

Sommario Facendo uso della teoria dei semigruppi si prova che un problema di trasporto di particelle con non-linearità quadratiche ha un'unica soluzione forte nello spazio di Banach UCB.

     

Transition radiation in a nonlinear and infinite one-dimensional structure: a comparison of solution methods

Transition zones in railway tracks are locations with a significant variation of track properties (i.e. foundation stiffness) encountered near structures such as bridges and tunnels. Due to strong amplification of the track’s response, transition zones are prone to rapid degradation. To investigate the degradation mechanisms in transition zones, researchers have developed a multitude of models, some of them being very complex. This study compares three solution methods, namely an integral-transform method, a time-domain method, and a hybrid method, with the goal of solving these systems efficiently. The methods are compared in terms of accuracy, computational efficiency, and feasibility of application to more complex systems. The model employed in this paper consists of an infinite, inhomogeneous, and piecewise-linear 1-D structure subjected to a moving constant load. Although the 1-D model is not particularly demanding computationally, it is used to make qualitative observations as to which method is most suitable for the 2-D and 3-D models, which could lead to significant gains. Results show that all three methods can reach similar accuracy levels, and in doing so, the time-domain method is most computationally efficient. The integral-transform method appears to be efficient in dealing with frequency-dependent parameters, while the time-domain and hybrid methods are efficient in dealing with a smooth nonlinearity. For multi-dimensional models, if nonlinearities and inhomogeneities are considered throughout the depth, the time-domain method is likely to be most efficient; however, if nonlinearities and inhomogeneities are limited to the surface layers, the integral-transform and hybrid methods have the potential to be more efficient than the time-domain one. Finally, although the 1-D model presented in this study is mainly used to assess the three methods, it can also be used for preliminary designs of transition zones in railway tracks.