Analysis of a model for adhesive contact with thermal effects - - PowerPoint PPT Presentation

analysis of a model for adhesive contact with thermal
SMART_READER_LITE
LIVE PREVIEW

Analysis of a model for adhesive contact with thermal effects - - PowerPoint PPT Presentation

Derivation of the model The PDE system Existence Long-time analysis Analysis of a model for adhesive contact with thermal effects Riccarda Rossi (Universit` a di Brescia) joint work with Elena Bonetti (Universit` a di Pavia), Giovanna


slide-1
SLIDE 1

Derivation of the model The PDE system Existence Long-time analysis

Analysis of a model for adhesive contact with thermal effects

Riccarda Rossi (Universit` a di Brescia) joint work with Elena Bonetti (Universit` a di Pavia), Giovanna Bonfanti (Universit` a di Brescia), Mathematical Models and Analytical Problems for Special Materials, Brescia, July 09–11 2009

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-2
SLIDE 2

Derivation of the model The PDE system Existence Long-time analysis

Fr´ emond’s modeling of adhesive contact

Setting

A viscoelastic body Ω ⊂ R3 in adhesive contact with a rigid support on a (flat) part ΓCont of its boundary ∂Ω = ΓDir ∪ ΓNeu ∪ ΓCont.

Related contributions

.... Andrews, Cang´ emi, Chau, Cocou, Eck, Fern´ andez, Figuereido, Han, Jaruˇ sek, Klarbring, Krbec, Kuttler, Martins, Mu˜ noz-Rivera, Point, Racke, Raous, Shi, Shillor, Sofonea, Telega, Trabucho, Wright.....

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-3
SLIDE 3

Derivation of the model The PDE system Existence Long-time analysis

Fr´ emond’s modeling of adhesive contact

Setting

A viscoelastic body Ω ⊂ R3 in adhesive contact with a rigid support on a (flat) part ΓCont of its boundary ∂Ω = ΓDir ∪ ΓNeu ∪ ΓCont.

Related contributions

.... Andrews, Cang´ emi, Chau, Cocou, Eck, Fern´ andez, Figuereido, Han, Jaruˇ sek, Klarbring, Krbec, Kuttler, Martins, Mu˜ noz-Rivera, Point, Racke, Raous, Shi, Shillor, Sofonea, Telega, Trabucho, Wright.....

Fr´ emond’s approach

Based on [M. Fr´

emond, Non-smooth Thermomechanics, 2002]

Account for microscopic motions in the macroscopic predictive theory

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-4
SLIDE 4

Derivation of the model The PDE system Existence Long-time analysis

Fr´ emond’s modeling of adhesive contact

Setting

A viscoelastic body Ω ⊂ R3 in adhesive contact with a rigid support on a (flat) part ΓCont of its boundary ∂Ω = ΓDir ∪ ΓNeu ∪ ΓCont.

Related contributions

.... Andrews, Cang´ emi, Chau, Cocou, Eck, Fern´ andez, Figuereido, Han, Jaruˇ sek, Klarbring, Krbec, Kuttler, Martins, Mu˜ noz-Rivera, Point, Racke, Raous, Shi, Shillor, Sofonea, Telega, Trabucho, Wright.....

Fr´ emond’s approach

Based on [M. Fr´

emond, Non-smooth Thermomechanics, 2002]

Account for microscopic motions in the macroscopic predictive theory

◮ microscopic bonds are responsible for the adhesion, microscopic motions

lead to rupture

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-5
SLIDE 5

Derivation of the model The PDE system Existence Long-time analysis

Fr´ emond’s modeling of adhesive contact

Setting

A viscoelastic body Ω ⊂ R3 in adhesive contact with a rigid support on a (flat) part ΓCont of its boundary ∂Ω = ΓDir ∪ ΓNeu ∪ ΓCont.

Related contributions

.... Andrews, Cang´ emi, Chau, Cocou, Eck, Fern´ andez, Figuereido, Han, Jaruˇ sek, Klarbring, Krbec, Kuttler, Martins, Mu˜ noz-Rivera, Point, Racke, Raous, Shi, Shillor, Sofonea, Telega, Trabucho, Wright.....

Fr´ emond’s approach

Based on [M. Fr´

emond, Non-smooth Thermomechanics, 2002]

Account for microscopic motions in the macroscopic predictive theory

◮ microscopic bonds are responsible for the adhesion, microscopic motions

lead to rupture

◮ account for the power of the microscopic motions in the power of the

interior forces

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-6
SLIDE 6

Derivation of the model The PDE system Existence Long-time analysis

State variables

In the isothermal case [Bonetti-Bonfanti-R. ’07,’08]

◮ in the volume domain Ω:

◮ small deformation (ε(u) symm. linear. strain tensor) (small perturbation

assumption)

thermal effects (θ absolute temperature)

◮ on the contact surface ΓCont:

◮ adhesion (χ “phase parameter” related to the active bonds of the adhesion

“damage parameter”)

◮ effects of displacement (u|ΓCont trace of the displacement) ◮

thermal effects (θs absolute temperature)

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-7
SLIDE 7

Derivation of the model The PDE system Existence Long-time analysis

State variables

To account for thermal effects: [Bonetti-Bonfanti-R. preprint’08]

◮ in the volume domain Ω:

◮ small deformation (ε(u) symm. linear. strain tensor) (small perturbation

assumption)

thermal effects (θ absolute temperature)

◮ on the contact surface ΓCont:

◮ adhesion (χ “phase parameter” related to the active bonds of the adhesion

“damage parameter”)

◮ effects of displacement (u|ΓCont trace of the displacement) ◮

thermal effects (θs absolute temperature)

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-8
SLIDE 8

Derivation of the model The PDE system Existence Long-time analysis

The equations for u and χ

From the principle of virtual power (interior & exterior forces, no acceleration forces)

◮ momentum balance:

8 > > < > > :

− div Σ = f in Ω × (0, T), 8 > < > : Σn = R in ΓCont × (0, T), u = 0 in ΓDir × (0, T), Σn = g in ΓNeu × (0, T),

8 < :

Σ stress tensor R reaction on the contact surface f volume force, g traction

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-9
SLIDE 9

Derivation of the model The PDE system Existence Long-time analysis

The equations for u and χ

From the principle of virtual power (interior & exterior forces, no acceleration forces)

◮ momentum balance:

8 > > < > > :

− div Σ = f in Ω × (0, T), 8 > < > : Σn = R in ΓCont × (0, T), u = 0 in ΓDir × (0, T), Σn = g in ΓNeu × (0, T),

8 < :

Σ stress tensor R reaction on the contact surface f volume force, g traction

◮ equation for the microscopic motions:

B − div H = 0 in ΓCont × (0, T), H · ns = 0 on ∂ΓCont × (0, T),

B interior microscopic work H microscopic work flux vector

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-10
SLIDE 10

Derivation of the model The PDE system Existence Long-time analysis

The equations for θ and θs

Entropy balance for θ and θs:

◮ for θ:

8 < :

st + div Q = h in Ω × (0, T), ( Q · n = F on ΓCont × (0, T), Q · n = 0 on ∂Ω \ ΓCont × (0, T),

8 > > < > > :

s internal entropy Q entropy flux vector h entropy source, F entropy flux through ΓCont

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-11
SLIDE 11

Derivation of the model The PDE system Existence Long-time analysis

The equations for θ and θs

Entropy balance for θ and θs:

◮ for θ:

8 < :

st + div Q = h in Ω × (0, T), ( Q · n = F on ΓCont × (0, T), Q · n = 0 on ∂Ω \ ΓCont × (0, T),

8 > > < > > :

s internal entropy Q entropy flux vector h entropy source, F entropy flux through ΓCont

◮ for θs :

∂tss + div Qs = F in Ω × (0, T), Qs · ns = 0 on ∂ΓCont × (0, T),

ss contact surface entropy Qs surface entropy flux vector

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-12
SLIDE 12

Derivation of the model The PDE system Existence Long-time analysis

The equations for θ and θs

Entropy balance for θ and θs:

◮ for θ:

8 < :

st + div Q = h in Ω × (0, T), ( Q · n = F on ΓCont × (0, T), Q · n = 0 on ∂Ω \ ΓCont × (0, T),

8 > > < > > :

s internal entropy Q entropy flux vector h entropy source, F entropy flux through ΓCont

◮ for θs :

∂tss + div Qs = F in Ω × (0, T), Qs · ns = 0 on ∂ΓCont × (0, T),

ss contact surface entropy Qs surface entropy flux vector

Entropy balance: obtained by rescaling the internal energy balance (under small perturbation assumpt.): see [Bonetti-Fr´

emond’03, Bonetti-Colli-Fr´ emond’03, Bonetti’06, Bonetti-Colli-Fabrizio-Gilardi’06,’07,08, Bonetti-Rocca-Fr´ emond’07]

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-13
SLIDE 13

Derivation of the model The PDE system Existence Long-time analysis

Constitutive laws

Constitutive relations for Σ, R, s, Q, F, B, H, ss, Qs derive from the volume & surface free energies ΨΩ = ΨΩ(u, θ), ΨΓCont = ΨΓCont(u|ΓCont , χ, θs) and the pseudo-potentials of dissipation ΦΩ = ΦΩ(∇θ, ε(ut)), ΦΓCont = ΦΓCont(∇θs, χt, θ|ΓCont − θs)

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-14
SLIDE 14

Derivation of the model The PDE system Existence Long-time analysis

Constitutive laws

Constitutive relations for Σ, R, s, Q, F, B, H, ss, Qs derive from the volume & surface free energies ΨΩ = ΨΩ(u, θ), ΨΓCont = ΨΓCont(u|ΓCont , χ, θs) and the pseudo-potentials of dissipation ΦΩ = ΦΩ(∇θ, ε(ut)), ΦΓCont = ΦΓCont(∇θs, χt, θ|ΓCont − θs) The energy potentials include constraints on the variables for physical consistency: nonsmooth (multivalued) operators in the constitutive eqns

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-15
SLIDE 15

Derivation of the model The PDE system Existence Long-time analysis

Constitutive laws

Constitutive relations for Σ, R, s, Q, F, B, H, ss, Qs derive from the volume & surface free energies ΨΩ = ΨΩ(u, θ), ΨΓCont = ΨΓCont(u|ΓCont , χ, θs) and the pseudo-potentials of dissipation ΦΩ = ΦΩ(∇θ, ε(ut)), ΦΓCont = ΦΓCont(∇θs, χt, θ|ΓCont − θs) The energy potentials include constraints on the variables for physical consistency: nonsmooth (multivalued) operators in the constitutive eqns

Constraints

◮ admissible values for χ and (possibly) irreversibility ◮ impenetrability condition between the body and the support ◮ positivity of the absolute temperatures θ and θs

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-16
SLIDE 16

Derivation of the model The PDE system Existence Long-time analysis

The adhesion phenomenon

The “damage parameter” χ denotes the fraction of active glue fibers at each point of the contact surface

◮ χ = 0 no adhesion (completely broken bonds) ◮ χ = 1 active adhesion (unbroken bonds)

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-17
SLIDE 17

Derivation of the model The PDE system Existence Long-time analysis

The adhesion phenomenon

The “damage parameter” χ denotes the fraction of active glue fibers at each point of the contact surface

◮ χ = 0 no adhesion (completely broken bonds) ◮ χ = 1 active adhesion (unbroken bonds)

Physical constraints

◮ χ ∈ [0, 1] ◮ χt ≤ 0 (irreversible phenomenon)

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-18
SLIDE 18

Derivation of the model The PDE system Existence Long-time analysis

The adhesion phenomenon

The “damage parameter” χ denotes the fraction of active glue fibers at each point of the contact surface

◮ χ = 0 no adhesion (completely broken bonds) ◮ χ = 1 active adhesion (unbroken bonds)

Physical constraints

◮ χ ∈ [0, 1]

As a first step, neglect irreversibility (∼ fresh, liquid glue..)

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-19
SLIDE 19

Derivation of the model The PDE system Existence Long-time analysis

Unilateral conditions

The impenetrability condition u|ΓCont · n ≤ 0

  • n ΓCont.

is ensured by the reaction on the contact surface

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-20
SLIDE 20

Derivation of the model The PDE system Existence Long-time analysis

Unilateral conditions

The impenetrability condition u|ΓCont · n ≤ 0

  • n ΓCont.

is ensured by the reaction on the contact surface If χ = 0 (no adhesion) the reaction on ΓCont is R = −χu − ∂I]−∞,0](u · n)n (omitting traces) where ∂I(−∞,0](u · n) = ( if u · n < 0 [0, +∞[ if u · n = 0

◮ R is normal to ΓCont ◮ R = 0 if u · n < 0 ◮ R = γn, γ ≤ 0 if u · n = 0

This is in agreement with the Signorini conditions

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-21
SLIDE 21

Derivation of the model The PDE system Existence Long-time analysis

Unilateral conditions

The impenetrability condition u|ΓCont · n ≤ 0

  • n ΓCont.

is ensured by the reaction on the contact surface In the case the adhesion is active χ > 0 R = −χu − ∂I]−∞,0](u · n)n i.e., there is a reaction (with rigidity ∼ χ) counteracting separation: R · n = −χ u · n > 0 if u · n < 0

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-22
SLIDE 22

Derivation of the model The PDE system Existence Long-time analysis

The PDE system: the momentum balance

Recall Ω ⊂ R3 smooth, bounded and ∂Ω = ΓDir ∪ ΓNeu ∪ ΓCont

◮ The momentum balance

−div (Kε(u) + Kvε(ut) + θ1) = f in Ω × (0, T) where: K elasticity tensor, Kv viscosity tensor, θ1 ↔ thermal deformation

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-23
SLIDE 23

Derivation of the model The PDE system Existence Long-time analysis

The PDE system: the momentum balance

Recall Ω ⊂ R3 smooth, bounded and ∂Ω = ΓDir ∪ ΓNeu ∪ ΓCont

◮ The momentum balance

−div (Kε(u) + Kvε(ut) + θ1) = f in Ω × (0, T) where: K elasticity tensor, Kv viscosity tensor, θ1 ↔ thermal deformation

◮ the boundary conditions

u = 0

  • n ΓDir,

(Kε(u) + Kvε(ut) + θ1)n = g

  • n ΓNeu × (0, T),

(Kε(u) + Kvε(ut) + θ1)n + χu + ∂I]−∞,0](u · n)n ∋ 0

  • n ΓCont × (0, T),

where −χu − ∂I]−∞,0](u · n)n on ΓCont is the reaction

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-24
SLIDE 24

Derivation of the model The PDE system Existence Long-time analysis

The PDE system: the evolution of the adhesion

We consider on ΓCont χt − ∆χ + ∂I[0,1](χ) ∋ −λ′(χ)(θs − θeq) − 1 2|u|2 in ΓCont × (0, T) ∂nχ = 0,

  • n ∂ΓCont × (0, T)

◮ ∂I[0,1](χ) ⇒ χ ∈ [0, 1]

(physical consistency) ∂I[0,1](χ) = 8 > < > : ] − ∞, 0] if χ = 0 if 0 < χ < 1 [0, +∞[ if χ = 1

◮ λ (quadratic) function, related to the latent heat ◮ θs temperature of the glue, θeq constant ◮ − 1 2|u|2 source of damage due to macroscopic movements

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-25
SLIDE 25

Derivation of the model The PDE system Existence Long-time analysis

The PDE system: the temperature equations

◮ The entropy equation (rescaled energy balance) for θ in Ω

∂t(log θ) − div ut − ∆θ = h in Ω × (0, T), ∂nθ = (

  • n ∂Ω \ ΓCont × (0, T)

−χ(θ|ΓCont − θs) on ΓCont × (0, T)

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-26
SLIDE 26

Derivation of the model The PDE system Existence Long-time analysis

The PDE system: the temperature equations

◮ The entropy equation (rescaled energy balance) for θ in Ω

∂t(log θ) − div ut − ∆θ = h in Ω × (0, T), ∂nθ = (

  • n ∂Ω \ ΓCont × (0, T)

−χ(θ|ΓCont − θs) on ΓCont × (0, T)

◮ The entropy equation for θs on ΓCont is

∂t(log θs) − λ′(χ)χt − ∆θs = χ(θ|ΓCont − θs) in ΓCont × (0, T) ∂nθs = 0

  • n ∂ΓCont × (0, T)

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-27
SLIDE 27

Derivation of the model The PDE system Existence Long-time analysis

The complete PDE system: difficulties

− div (Kε(u) + Kvε(ut) + θ1) = f in Ω × (0, T), u = 0

  • n ΓDir × (0, T),

(Kε(u) + Kvε(ut) + θ1)n = g

  • n ΓNeu × (0, T),

(Kε(u) + Kvε(ut) + θ1)n + χu + ∂I]−∞,0](u · n)n ∋ 0

  • n ΓCont × (0, T),

χt − ∆χ + ∂I[0,1](χ) ∋ −λ′(χ)(θs − θeq) − 1 2|u|2 in ΓCont × (0, T), ∂nχ = 0

  • n ∂ΓCont × (0, T)

∂t(log θ) − div ut − ∆θ = h in Ω × (0, T), ∂nθ = (

  • n ∂Ω \ ΓCont × (0, T),

−χ(θ − θs) on ΓCont × (0, T), ∂t(log θs) − λ′(χ)χt − ∆θs = χ(θ − θs) in ΓCont × (0, T), ∂nθs = 0

  • n ∂ΓCont × (0, T),

+Cauchy conditions → singular character of the θ, θs-equations (θ-equation is coupled with a third type boundary condition) → we deduce directly θ, θs > 0, crucial for thermodynamical consistency!

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-28
SLIDE 28

Derivation of the model The PDE system Existence Long-time analysis

The complete PDE system: difficulties

− div (Kε(u) + Kvε(ut) + θ1) = f in Ω × (0, T), u = 0 in ΓDir × (0, T), (Kε(u) + Kvε(ut) + θ1)n = g in ΓNeu × (0, T), (Kε(u) + Kvε(ut) + θ1)n + χu + ∂I]−∞,0](u · n)n ∋ 0 in ΓCont × (0, T), χt − ∆χ + ∂I[0,1](χ) ∋ −λ′(χ)(θs − θeq) − 1 2|u|2 in ΓCont × (0, T), ∂nχ = 0

  • n ∂ΓCont × (0, T)

∂t(log θ) − div ut − ∆θ = h in Ω × (0, T), ∂nθ = (

  • n ∂Ω \ ΓCont × (0, T),

−χ(θ − θs) on ΓCont × (0, T), ∂t(log θs) − λ′(χ)χt − ∆θs = χ(θ − θs) in ΓCont × (0, T) ∂nθs = 0

  • n ∂ΓCont × (0, T)

+ Cauchy conditions → (quadratic) coupling terms on the boundary → we need sufficient regularity on θ and u to control their traces.

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-29
SLIDE 29

Derivation of the model The PDE system Existence Long-time analysis

The complete PDE system: difficulties

− div (Kε(u) + Kvε(ut) + θ1) = f in Ω × (0, T), u = 0

  • n ΓDir × (0, T),

(Kε(u) + Kvε(ut) + θ1)n = g

  • n ΓNeu × (0, T),

(Kε(u) + Kvε(ut) + θ1)n + χu + ∂I]−∞,0](u · n)n ∋ 0

  • n ΓCont × (0, T),

χt − ∆χ + ∂I[0,1](χ) ∋ −λ′(χ)(θs − θeq) − 1 2|u|2 in ΓCont × (0, T), ∂nχ = 0

  • n ∂ΓCont × (0, T),

∂t(log θ) − div ut − ∆θ = h in Ω × (0, T), ∂nθ = (

  • n ∂Ω \ ΓCont × (0, T),

−χ(θ − θs) on ΓCont × (0, T), ∂t(log θs) − λ′(χ)χt − ∆θs = χ(θ − θs) in ΓCont × (0, T), ∂nθs = 0

  • n ∂ΓCont × (0, T),

+Cauchy conditions → non-smooth (multivalued) constraints on χ and u · n

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-30
SLIDE 30

Derivation of the model The PDE system Existence Long-time analysis

The existence theorem

Theorem. Given initial data ( θ0 ∈ L1(Ω), log(θ0) ∈ L2(Ω), θs,0 ∈ L1(ΓCont), log(θs,0) ∈ L2(ΓCont), u0 ∈ H1(Ω)3, χ0 ∈ H1(ΓCont) there exist (u, χ, θ, θs) solving the weak, variational formulation of the IBVP

− div (Kε(u) + Kv ε(ut ) + θ1) = f in Ω × (0, T), (Kε(u) + Kv ε(ut ) + θ1)n + χu + ∂I]−∞,0](u · n)n ∋ 0

  • n Γc × (0, T),

χt − ∆χ + ∂I[0,1](χ) ∋ −λ′(χ)(θs − θeq) − 1 2 |u|2 in Γc × (0, T), ∂t (log θ) − div ut − ∆θ = h in Ω × (0, T), ∂nθ = (

  • n ∂Ω \ Γc × (0, T),

−χ(θ − θs ) on Γc × (0, T), ∂t (log θs ) − λ′(χ)χt − ∆θs = χ(θ − θs ) in Γc × (0, T), Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-31
SLIDE 31

Derivation of the model The PDE system Existence Long-time analysis

The existence theorem

Theorem. Given initial data ( θ0 ∈ L1(Ω), log(θ0) ∈ L2(Ω), θs,0 ∈ L1(ΓCont), log(θs,0) ∈ L2(ΓCont), u0 ∈ H1(Ω)3, χ0 ∈ H1(ΓCont) there exist (u, χ, θ, θs) solving the weak, variational formulation of the IBVP

− div (Kε(u) + Kv ε(ut ) + θ1) = f in Ω × (0, T), (Kε(u) + Kv ε(ut ) + θ1)n + χu + ∂I]−∞,0](u · n)n ∋ 0

  • n Γc × (0, T),

χt − ∆χ + ∂I[0,1](χ) ∋ −λ′(χ)(θs − θeq) − 1 2 |u|2 in Γc × (0, T), ∂t (log θ) − div ut − ∆θ = h in Ω × (0, T), ∂nθ = (

  • n ∂Ω \ Γc × (0, T),

−χ(θ − θs ) on Γc × (0, T), ∂t (log θs ) − λ′(χ)χt − ∆θs = χ(θ − θs ) in Γc × (0, T),

with u ∈ H1(0, T; H1(Ω)3) χ ∈ L2(0, T; H2(ΓCont)) ∩ L∞(0, T; H1(ΓCont)) ∩ H1(0, T; L2(ΓCont)) θ ∈ L2(0, T; H1(Ω)) ∩ L∞(0, T; L1(Ω)), log θ ∈ L∞(0, T; L2(Ω)) ∩ H1(0, T; H1(Ω)′) θs ∈ L2(0, T; H1(ΓCont)) ∩ L∞(0, T; L1(ΓCont)), log θs ∈ L∞(0, T; L2(ΓCont)) ∩ H1(0, T; H1(ΓCont)′)

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-32
SLIDE 32

Derivation of the model The PDE system Existence Long-time analysis

The existence theorem

Theorem. Given initial data ( θ0 ∈ L1(Ω), log(θ0) ∈ L2(Ω), θs,0 ∈ L1(ΓCont), log(θs,0) ∈ L2(ΓCont), u0 ∈ H1(Ω)3, χ0 ∈ H1(ΓCont) there exist (u, χ, θ, θs) solving the weak, variational formulation of the IBVP

− div (Kε(u) + Kv ε(ut ) + θ1) = f in Ω × (0, T), (Kε(u) + Kv ε(ut ) + θ1)n + χu + ∂I]−∞,0](u · n)n ∋ 0

  • n Γc × (0, T),

χt − ∆χ + ∂I[0,1](χ) ∋ −λ′(χ)(θs − θeq) − 1 2 |u|2 in Γc × (0, T), ∂t (log θ) − div ut − ∆θ = h in Ω × (0, T), ∂nθ = (

  • n ∂Ω \ Γc × (0, T),

−χ(θ − θs ) on Γc × (0, T), ∂t (log θs ) − λ′(χ)χt − ∆θs = χ(θ − θs ) in Γc × (0, T),

with u ∈ H1(0, T; H1(Ω)3) χ ∈ L2(0, T; H2(ΓCont)) ∩ L∞(0, T; H1(ΓCont)) ∩ H1(0, T; L2(ΓCont)) θ ∈ L2(0, T; H1(Ω)) ∩ L∞(0, T; L1(Ω)), log θ ∈ L∞(0, T; L2(Ω)) ∩ H1(0, T; H1(Ω)′) θs ∈ L2(0, T; H1(ΓCont)) ∩ L∞(0, T; L1(ΓCont)), log θs ∈ L∞(0, T; L2(ΓCont)) ∩ H1(0, T; H1(ΓCont)′)

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-33
SLIDE 33

Derivation of the model The PDE system Existence Long-time analysis

First step of the proof: approximation

  • Main difficulty: lack of regularity for θ (presence of log θ and boundary

conditions)

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-34
SLIDE 34

Derivation of the model The PDE system Existence Long-time analysis

First step of the proof: approximation

  • Main difficulty: lack of regularity for θ (presence of log θ and boundary

conditions) ε∂tθ+∂t(log θ) − div ut − ∆θ−ε∆ (∂tθ) = h in Ω × (0, T), ∂nθ = (

  • n ∂Ω \ ΓCont × (0, T)

−χ(θ − θs) on ΓCont × (0, T)

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-35
SLIDE 35

Derivation of the model The PDE system Existence Long-time analysis

First step of the proof: approximation

  • Main difficulty: lack of regularity for θ (presence of log θ and boundary

conditions) ε∂tθ+∂t(logεθ) − div ut − ∆θ−ε∆ (∂tθ) = h in Ω × (0, T), ∂nθ = (

  • n ∂Ω \ ΓCont × (0, T)

−χ(θ − θs) on ΓCont × (0, T)

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-36
SLIDE 36

Derivation of the model The PDE system Existence Long-time analysis

First step of the proof: approximation

  • Main difficulty: lack of regularity for θ (presence of log θ and boundary

conditions) ε∂tθ+∂t(logεθ) − div ut − ∆θ−ε∆ (∂tθ) = h in Ω × (0, T), ∂nθ = (

  • n ∂Ω \ ΓCont × (0, T)

−χ(θ − θs) on ΓCont × (0, T) where logε denotes a Yosida-type approximation of log

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-37
SLIDE 37

Derivation of the model The PDE system Existence Long-time analysis

First step of the proof: approximation

  • Main difficulty: lack of regularity for θ (presence of log θ and boundary

conditions) ε∂tθ+∂t(logεθ) − div ut − ∆θ−ε∆ (∂tθ) = h in Ω × (0, T), ∂nθ = (

  • n ∂Ω \ ΓCont × (0, T)

−χ(θ − θs) on ΓCont × (0, T) where logε denotes a Yosida-type approximation of log

  • Analogously for θs:

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-38
SLIDE 38

Derivation of the model The PDE system Existence Long-time analysis

First step of the proof: approximation

  • Main difficulty: lack of regularity for θ (presence of log θ and boundary

conditions) ε∂tθ+∂t(logεθ) − div ut − ∆θ−ε∆ (∂tθ) = h in Ω × (0, T), ∂nθ = (

  • n ∂Ω \ ΓCont × (0, T)

−χ(θ − θs) on ΓCont × (0, T) where logε denotes a Yosida-type approximation of log

  • Analogously for θs:

ε∂tθs+∂t(log θs) − λ′(χ)χt − ∆θs−ε∆ (∂tθs) = χ(θ − θs) in ΓCont × (0, T), ∂nθs = 0

  • n ∂ΓCont × (0, T)

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-39
SLIDE 39

Derivation of the model The PDE system Existence Long-time analysis

First step of the proof: approximation

  • Main difficulty: lack of regularity for θ (presence of log θ and boundary

conditions) ε∂tθ+∂t(logεθ) − div ut − ∆θ−ε∆ (∂tθ) = h in Ω × (0, T), ∂nθ = (

  • n ∂Ω \ ΓCont × (0, T)

−χ(θ − θs) on ΓCont × (0, T) where logε denotes a Yosida-type approximation of log

  • Analogously for θs:

ε∂tθs+∂t(logεθs) − λ′(χ)χt − ∆θs−ε∆ (∂tθs) = χ(θ − θs) in ΓCont × (0, T), ∂nθs = 0

  • n ∂ΓCont × (0, T)

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-40
SLIDE 40

Derivation of the model The PDE system Existence Long-time analysis

First step of the proof: approximation

  • Main difficulty: lack of regularity for θ (presence of log θ and boundary

conditions) ε∂tθ+∂t(logεθ) − div ut − ∆θ−ε∆ (∂tθ) = h in Ω × (0, T), ∂nθ = (

  • n ∂Ω \ ΓCont × (0, T)

−χ(θ − θs) on ΓCont × (0, T) where logε denotes a Yosida-type approximation of log

  • Analogously for θs:

ε∂tθs+∂t(logεθs) − λ′(χ)χt − ∆θs−ε∆ (∂tθs) = χ(θ − θs) in ΓCont × (0, T), ∂nθs = 0

  • n ∂ΓCont × (0, T)
  • we do not approximate ∂I[0,1](χ), ∂I]−∞,0](u · n) in the eq.’s for u and χ

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-41
SLIDE 41

Derivation of the model The PDE system Existence Long-time analysis

Outline of the proof

Existence:

◮ approximating procedure

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-42
SLIDE 42

Derivation of the model The PDE system Existence Long-time analysis

Outline of the proof

Existence:

◮ approximating procedure ◮ Schauder theorem to solve the approximated problem

◮ theory of evolution equations with maximal monotone operators for the

equations for θ and θs

◮ use previous results of [Bonetti, Bonfanti, R. ’07, ’08] for the system in u

and χ

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-43
SLIDE 43

Derivation of the model The PDE system Existence Long-time analysis

Outline of the proof

Existence:

◮ approximating procedure ◮ Schauder theorem to solve the approximated problem

◮ theory of evolution equations with maximal monotone operators for the

equations for θ and θs

◮ use previous results of [Bonetti, Bonfanti, R. ’07, ’08] for the system in u

and χ

◮ uniform a priori estimates (... some cancellations work, positivity of χ...)

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-44
SLIDE 44

Derivation of the model The PDE system Existence Long-time analysis

Outline of the proof

Existence:

◮ approximating procedure ◮ Schauder theorem to solve the approximated problem

◮ theory of evolution equations with maximal monotone operators for the

equations for θ and θs

◮ use previous results of [Bonetti, Bonfanti, R. ’07, ’08] for the system in u

and χ

◮ uniform a priori estimates (... some cancellations work, positivity of χ...) ◮ passage to the limit by compactness and semicontinuity arguments

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-45
SLIDE 45

Derivation of the model The PDE system Existence Long-time analysis

Outline of the proof

Existence:

◮ approximating procedure ◮ Schauder theorem to solve the approximated problem

◮ theory of evolution equations with maximal monotone operators for the

equations for θ and θs

◮ use previous results of [Bonetti, Bonfanti, R. ’07, ’08] for the system in u

and χ

◮ uniform a priori estimates (... some cancellations work, positivity of χ...) ◮ passage to the limit by compactness and semicontinuity arguments

Uniqueness?

◮ NOT expected due to the nonlinear structure of the equations, the lack

  • f regularity of θ and the boundary conditions

◮ holds for the approximating problem

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-46
SLIDE 46

Derivation of the model The PDE system Existence Long-time analysis

Trajectories on (0, +∞)

There exists a solution in (0, T) for any T: we can suitably define the notion of solution on (0, +∞)

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-47
SLIDE 47

Derivation of the model The PDE system Existence Long-time analysis

Trajectories on (0, +∞)

There exists a solution in (0, T) for any T: we can suitably define the notion of solution on (0, +∞) How the trajectories (u(t), χ(t), θ(t), θs(t)) behave as t → +∞?

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-48
SLIDE 48

Derivation of the model The PDE system Existence Long-time analysis

Trajectories on (0, +∞)

There exists a solution in (0, T) for any T: we can suitably define the notion of solution on (0, +∞) How the trajectories (u(t), χ(t), θ(t), θs(t)) behave as t → +∞?

◮ need uniform estimates on the solutions independent of the final time T ◮ further requirement on the entropy flux through Γc

∂nθ = −(χ + c)(θ|Γc − θs) , c > 0 on Γc residual flux even if χ = 0 (crucial to obtain a L2-estimate for θ|Γc − θs)

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-49
SLIDE 49

Derivation of the model The PDE system Existence Long-time analysis

Long-time a priori estimates

Due to the dissipative character of the system |θ|L∞(1,+∞;H1(Ω)) + |∇θ|L2(0,+∞;L2(Ω)) ≤ C |θt|L2(1,+∞;Lp(Ω)) ≤ C

for p ≤ 12/7

|θs|L∞(1,+∞;H1(ΓCont)) + |∇θs|L2(0,+∞;L2(ΓCont)) ≤ C |∂tθs|L2(1,+∞;Lq(ΓCont)) ≤ C,

for q < 2

|θ|ΓCont − θs|L2(0,+∞;L2(ΓCont)) ≤ C |χ|L∞(1,+∞;H2(ΓCont)) + |χt|L2(1,+∞;H1(ΓCont)) ≤ C |u|L∞(0,+∞;H1(Ω)3) + |ut|L2(0,+∞;H1(Ω)3) ≤ C |∂t(log θ)|L2(0,+∞;H1(Ω)′) + |∂t(log θs)|L2(0,+∞;H1(ΓCont)′) ≤ C the “energy” and the “dissipation” are uniformly bounded the solutions trajectories converge in a suitable sense to some cluster points as t∞

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-50
SLIDE 50

Derivation of the model The PDE system Existence Long-time analysis

The ω-limit set

The set of the possible cluster points (u∞, χ∞, θ∞, θs∞) of the solutions trajectories

◮ is non-empty, connected and compact (w.r.t. to a suitable topology)

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-51
SLIDE 51

Derivation of the model The PDE system Existence Long-time analysis

The ω-limit set

The set of the possible cluster points (u∞, χ∞, θ∞, θs∞) of the solutions trajectories

◮ is non-empty, connected and compact (w.r.t. to a suitable topology) ◮ its elements (u∞, χ∞, θ∞, θs∞) solve the stationary problem

div(Kε(u∞) + θ∞1) = f∞ in Ω + boundary conditions − ∆χ∞ + ∂I[0,1](χ∞) ∋ −λ′(χ∞)(θs∞ − θeq) − 1 2|u∞|2 in ΓCont + b.c. ∃ ¯ θ∞ ≥ 0 : θ∞ ≡ ¯ θ∞ in Ω , θs∞ ≡ ¯ θ∞ in ΓCont

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects

slide-52
SLIDE 52

Derivation of the model The PDE system Existence Long-time analysis

The ω-limit set

The set of the possible cluster points (u∞, χ∞, θ∞, θs∞) of the solutions trajectories

◮ is non-empty, connected and compact (w.r.t. to a suitable topology) ◮ its elements (u∞, χ∞, θ∞, θs∞) solve the stationary problem

div(Kε(u∞) + θ∞1) = f∞ in Ω + boundary conditions − ∆χ∞ + ∂I[0,1](χ∞) ∋ −λ′(χ∞)(θs∞ − θeq) − 1 2|u∞|2 in ΓCont + b.c. ∃ ¯ θ∞ ≥ 0 : θ∞ ≡ ¯ θ∞ in Ω , θs∞ ≡ ¯ θ∞ in ΓCont

◮ thermomechanical equilibrium (no dissipation) in the limit t → ∞.

Riccarda Rossi Analysis of a model for adhesive contact with thermal effects