Thermomechanical formulation of ductile damage coupled to nonlinear isotropic hardening and multiplicative viscoplasticity

In this paper, we present a thermomechanical framework which makes use of the internal variable theory of thermodynamics for damage-coupled finite viscoplasticity with nonlinear isotropic hardening. Damage evolution, being an irreversible process, generates heat. In addition to its direct effect on material's strength and stiffness, it causes deterioration of the heat conduction. The formulation, following the footsteps of Simó and Miehe (1992), introduces inelastic entropy as an additional state variable. Given a temperature dependent damage dissipation potential, we show that the evolution of inelastic entropy assumes a split form relating to plastic and damage parts, respectively. The solution of the thermomechanical problem is based on the so-called isothermal split. This allows the use of the model in 2D and 3D example problems involving geometrical imperfection triggered necking in an axisymmetric bar and thermally triggered necking of a 3D rectangular bar.     

A note on the two temperature theory with dual-phase-lag delay: Some exact solutions

In this note we present exact solutions of two initial-boundary value problems (IBVP)s in the setting of a recently-introduced theory of heat conduction, wherein the two temperature theory of the late 1960s is merged with Tzou’s dual-phase-lag flux relation. First, we solve a one-dimensional problem on a finite interval for a simple, parabolic initial condition. We then describe how to extend the analysis to the general three-dimensional case. In particular, it is demonstrated that the instability which generally arises in connection with the dual-phase-lag model can be avoided under this hybrid formulation.     

Numerical simulation of the interface molten metal air in the shot sleeve chambre and mold cavity of a die casting machine

The objective of this study relates to the numerical simulation of the free surface during the two-dimensional flow and solidification of aluminum in the horizontal cylinder and mold cavity of the high pressure die casting HPDC machine with cold chamber. The flow is governed by the Navier–Stokes equations (the mass and the momentum conservations) and solved in the two phase’s liquid aluminum and air. The tracking of the free surface is ensured by the VOF method. The equivalent specific heat method is used to solve the phase change heat transfer problem in the solidification process. Considering the displacement of the plunger, the geometry of the problem is variable and the numerical resolution uses a dynamic grid. The study examines the influence of the plunger speed on the evolution of the interface aluminum liquid–air profile, the mass of air imprisoned and the stream function contours versus time. Filling of a mold is an essential part of HPDC process and affects significantly the heat transfer and solidification of the melt. For this reason, accurate prediction of the temperature field in the system can be achieved only by including simulation of filling in the analysis.     

Efficient reconstruction of local heat fluxes in pool boiling experiments by goal-oriented adaptive mesh refinement

In this paper, we consider the efficient estimation of local boiling heat fluxes from transient temperature measurements in the heater close to the heater surface. For accurate prediction, heat flux estimation is formulated as a transient three-dimensional (3D) inverse heat conduction problem (IHCP). This inverse problem is ill-posed and cannot be treated straightforwardly by established numerical methods. In order to obtain a regularized stable solution, a large-scale time-dependent PDE-constrained optimization problem has to be solved and an appropriate stopping criterion for the termination of the iterative solution process has to be chosen. Since the boiling heat flux is non-uniformly distributed on the heater surface due to the strong local activity of the boiling process, the use of a fixed uniform spatial discretization is not efficient. Instead, an adaptive mesh refinement strategy can be used to obtain an appropriate discretization which significantly reduces the total computational effort. In this work, we present an automatic algorithm incorporating an adaptive mesh refinement via a heat flux-based a-posteriori error estimation technique. The suggested algorithm can cope with both spatially point-wise or highly resolved temperature observations efficiently. It is applied to real measurement data obtained from two different types of pool boiling experiments. The numerical results show that the computational effort can be reduced significantly for given estimation quality. This adaptive IHCP solution technique can be also viewed as an efficient soft sensor to deduce unmeasurable local boiling heat fluxes.     

A numerical study of natural convection from a localized heat source in the lower plenum of a fast breeder reactor under failed conditions

A comprehensive numerical study has been done to investigate two-dimensional, steady state, conjugate natural heat convection in the hemi spherical lower plenum of a fast breeder reactor under failed conditions. The continuity, momentum and energy equations are solved over the entire domain, using the corresponding properties for the solid and fluid regions. The control volume approach is employed in order to discretize the governing equations for their numerical solution. A parametric study has been done to study the variation of the velocity vectors and isotherms for different constant temperature of the heat source, simulating different heat generation rates. The actual problem in a nuclear reactor involves a volumetric heat generation in the debris falling over the heat shield plate under failed conditions of the reactor and heat is removed by a decay heat exchanger serving as a sink. In this study we have reduced this transient problem to a quasi-steady problem with a prescribed temperature on the heat shield plate. This makes the problem more tractable. The fluid flow pattern, variation of the temperature along the axis in and around the heat source are presented to show the overall heat transfer characteristics inside the plenum.     

Nulle contrôlabilité régionale pour des équations de la chaleur dégénéréesRegional null controllability for degenerate heat equations

We are interested in a null controllability problem for a class of strongly degenerate heat equations.First for all T>0, we prove a regional null controllability result at time T at least in the region where the equation is not degenerate. The proof is based on an adequate observability inequality for the homogeneous adjoint problem. This inequality is obtained by application of Carleman estimates combined with the introduction of cut-off functions.Then we improve this result: for all T′>T, we obtain a result of persistent regional null controllability during the time interval [T,T′]. Finally we give similar results for the (non degenerate) heat equation in unbounded domain. To cite this article: P. Cannarsa et al., C. R. Mecanique 330 (2002) 397–401.     

Two‐objective optimization strategies using the adjoint method and game theory in inverse natural convection problems

This paper considers the problem of estimating the strengths of two time‐varying heat sources simultaneously, from measurements of the temperature inside the square domain in a porous medium, when prior knowledge of the source functions is not available. This problem is an inverse natural convection problem. In order to circumvent this problem, we define several optimization criteria (objective functionals) that measure discrepancies between model and measured data, where objective functionals depend on two heat sources and use multi‐criteria optimization to identify Nash equilibria, which are solutions to the non‐cooperative game according to game theory. Two non‐cooperative game strategies are considered: competitive (Nash) game and hierarchical (modified Stackelberg) game. The methodology that we employ relies on a combination of mixed finite element space approximations, finite difference time discretizations, adjoint equation and sensitivity equation techniques, and nonlinear conjugate gradient algorithms for the solutions of estimating two heat sources. Applying the Sobolev gradient for the noise removal is investigated. The performance of the present technique of inverse analysis is evaluated, by means of several numerical experiments, and is found to be very accurate as well as efficient. Copyright © 2012 John Wiley & Sons, Ltd.     

Interaction between a rigid inclusion and a line crack under uniform heat flux

The thermoelastic displacement boundary value problem for a rigid inclusion interacting with a line crack in an infinite plane subjected to a uniform heat flux is studied, in which the rigid body rotation of the inclusion is considered. To solve the prescribed problem, we use the principle of superposition to decompose it into two groups of problems, which are further reduced to several basic subproblems including Green’s functions of edge dislocation and heat source couple, as well as the problem of a plane containing the inclusion under uniform heat flux and the problem of the inclusion subjected to a small rotation. The problems are solved using the complex variable method along with the rational mapping function technique. The variations of the stress intensity factors at the crack tips and the rigid body rotation angles with various crack lengths and heat flux angles are shown. The effects of the inclusion shape and size are also investigated.     

A note on the delayed heat equation: Instability with respect to initial data

Working in the context of a simple, one-dimensional, initial-boundary value problem involving homogeneous Dirichlet boundary data, we show that the time delayed heat equation can exhibit a type of instability with respect to the initial condition (IC); specifically, we show that a slight (in the L² sense) change in the IC can change a well-posed problem to an ill-posed one. We also establish that a physically realistic solution is possible only if the IC is of a (very) specific form. The main implication of this study is that the single and dual phase lag models, which have been put forward as possible alternatives to Fourier’s law, are not valid constitutive relations for the thermal flux vector.     

The numerical method for analysing the transient temperature field and transient phase transformation of metal materials in heat treating

In this paper,the Kirchhoffs transformation is popularized to the nonlinear heat conduction problem which the heat conductivity can be expressd as a multinomial of temperature firstly,the boundary condition of heat conduction problem is determined by analytics.Secondly,the incubation peroid superposition and the linear combination law is employed to simulate the transient phasses transformation in the process of heat treatment of materials.That the begin time of phase transformation,the type of phase transformation and the amount of phase constitution is determined simply.Finally,the three-dimension Dual Reciprocity Boundary Element Method is usedto analysis the total process of various heat treatment of component,the results of numerical calculation of examples show that the method provided in this paper is effectivce.     

Identification of a moving boundary for a heat conduction problem in a multilayer medium

In this paper, we give a uniqueness theorem for the moving boundary of a heat problem in a composite medium. Through solving the Cauchy problem of heat equation in each subdomain, we finally find an approximation to the moving boundary for one-dimensional heat conduction problem in a multilayer medium. The numerical scheme is based on the use of the method of fundamental solutions and a discrete Tikhonov regularization technique with the generalized cross-validation choice rule for a regularization parameter. Numerical experiments for five examples show that our proposed method is effective and stable.     

Adiabatic and diabatic two-phase venting flow in a microchannel

The growth and advection of the vapor phase in two-phase microchannel heat exchangers increase the system pressure and cause flow instabilities. One solution is to locally vent the vapor formed by capping the microchannels with a porous, hydrophobic membrane. In this paper we visualize this venting process in a single 124 μm by 98 μm copper microchannel with a 65 μm thick, 220 nm pore diameter hydrophobic Teflon membrane wall to determine the impact of varying flow conditions on the flow structures and venting process during adiabatic and diabatic operation. We characterize liquid velocities of 0.14, 0.36 and 0.65 m/s with superficial air velocities varying from 0.3 to 8 m/s. Wavy-stratified and stratified flow dominated low liquid velocities while annular type flows dominated at the higher velocities. Gas/vapor venting can be improved by increasing the venting area, increasing the trans-membrane pressure or using thinner, high permeability membranes. Diabatic experiments with mass flux velocities of 140 and 340 kg/s/m² and exit qualities up to 20% found that stratified type flows dominate at lower mass fluxes while churn-annular flow became more prevalent at the higher mass-flux and quality. The diabatic flow regimes are believed to significantly influence the pressure-drop and heat transfer coefficient in vapor venting heat exchangers.     

Numerical approximation of the heat transfer between domains separated by thin walls

In this paper, we analyse the numerical approximation of the heat transfer problem between two subdomains that we will consider filled with a fluid and separated by a thin solid wall. First of all, we state the problem in the whole domain with discontinuous physical properties. As an alternative and under certain assumptions on the separating walls, a classical Robin boundary condition between the fluid domains is obtained, thus eliminating the solid wall, and according to which the heat flux is proportional to the temperature difference between the two subdomains. Apart from discussing the relation between both approaches, we consider their numerical approximation, considering different alternatives for the first case, that is, the case in which temperatures are also computed in the solid wall. Copyright © 2006 John Wiley & Sons, Ltd.     

冻土水热力耦合作用的数学模型及数值模拟

将冻土体视为空间弹性体,提出了土体在冻结过程中水分场、温度场、应力场三场耦合的一般数学模型,并给出了相应的离散方程及其解法,最后给出了数值算例,并与实测值比较,证明了该模型和算法的正确性。     

Decoupled three-dimensional finite element computation of thermoelastic damping using Zener’s approximation

We consider three dimensional finite element computations of thermoelastic damping ratios of arbitrary bodies using Zener’s approach. In our small-damping formulation, unlike existing fully coupled formulations, the calculation is split into three smaller parts. Of these, the first sub-calculation involves routine undamped modal analysis using ANSYS. The second sub-calculation takes the mode shape, and solves on the same mesh a periodic heat conduction problem. Finally, the damping coefficient is a volume integral, evaluated elementwise. In the only other decoupled three dimensional computation of thermoelastic damping reported in the literature, the heat conduction problem is solved much less efficiently, using a modal expansion. We provide numerical examples using some beam-like geometries, for which Zener’s and similar formulas are valid. Among these we examine tapered beams, including the limiting case of a sharp tip. The latter’s higher-mode damping ratios dramatically exceed those of a comparable uniform beam.     

Fundamental problems in brittle fracture: unstable cracks and delayed breaking

This presents a short review on two problems where brittle fracture is involved. In the first one, a hot glass plate is subject to a local stress when drowned into a cold bath. In the region of transition between the cold and hot side, large stresses build up, that can be related accurately to the various coefficients typical of the plate material, heat dilation coefficient, heat conductivity, Young's modulus, etc. Thanks to thais, one can compare well the experimental results and the theoretical predictions based on the Griffith criterion for the propagation of a straight crack, and for its instability against an undulating mode. In the other problem, one looks at the delay before a bent 2D crystal breaks: this (very long) delay is interpreted as the time required for homogeneous nucleation of a critical Griffith nucleus in a region of the crystal under extension. Although it agrees fairly well with some experimental data, other experimental facts are required to complicate the model by considering a multistep nucleation process. To cite this article: Y. Pomeau, C. R. Mecanique 330 (2002) 249–257.     

Nonlinear oscillators with real valued powers: an analytic treatment

We give an analytic treatment of a second order ordinary differential equation describing a nonlinear oscillatory process with real valued power. This equation occurs in studying flows through porous media, heat conduction or plasma physics. After giving suitable conditions for the solubility of this equation in closed form, we tackle the related first boundary value problem and apply our analytic method to solve it. In the numerical section we discuss some models with fractional powers.     

采用包络-准则法优化考虑瞬态传热的薄板结构

刚度和强度是薄板结构的两个主要性能。在瞬态传热过程中,考虑热-力耦合,随时间和空间变化的非均匀温度场在结构中会引起热变形和热应力,温度场随时间变化的规律和空间分布依赖于板的厚度变化,进而影响板的刚度和强度。因此,考虑瞬态传热的薄板优化问题具有更强的非线性,更加难以求解。本文给出一种包络-准则方法处理这类结构优化问题。首先,针对外力荷载,进行一个结构柔顺性的优化设计;以这一设计为基础,通过瞬态热-力耦合分析及优化准则,计算多个时刻的优化设计变量并取其包络,对上述优化结果进行迭代修正,以消除瞬态温度场作用下较高的局部应力。优化算例表明,该方法对于考虑瞬态传热薄板优化问题有效。     

On the effect of the heat generated during a stress-induced thermoelastic phase transformation

In this paper we examine the stress-elongation response of a bar undergoing a thermoelastic phase transition. Attention is focussed on how this response is affected by the heat generated during the transformation. The analysis is based on a continuum model consisting of a two-well Helmholtz free-energy function, a kinetic relation and a nucleation criterion. The governing mathematical problem is related to one that describes a moving heat source, except that here, the strength and speed of the source are not knowna priori and the energy field equation involves coupling between thermal and mechanical effects. A finite difference solution of this moving boundary-value problem is carried out. The heat generated by the transformation is found to have a significant effect on the mechanical response whenever the prescribed elongation-rate is moderately large.     

Analytical solution to the unsteady one-dimensional conduction problem with two time-varying boundary conditions: Duhamel’s theorem and separation of variables

An analytical resolution of the time-dependent one-dimensional heat conduction problem with time-dependent boundary conditions using the method of separation of variables and Duhamel’s theorem is presented. The two boundary conditions used are a time-dependent heat flux at one end and a varying temperature at the other end of the one-dimensional domain. It is put forth because the author found that the prescribed resolution method using separation of variables and Duhamel’s theorem presented in heat conduction textbooks is not directly applicable to problems with more than one time-dependent boundary condition. The analytical method presented in this paper makes use of one of the property of the heat conduction equation: the apparent linearity of the solutions. For that reason, in order to solve a problem with two time-dependent boundary conditions, the author first separates the initial problem into two independent but complementary problems, each with only one time-dependent boundary condition. Doing that, both simpler problems can be solved independently using a prescribed method that is known to work and the final solution can be obtained by joining the two independent solutions from the simpler separated problems. Every step of the resolution method is presented in this paper, along with a numerical validation of the final solution of three test case problems.