# Density profile in latent heat model

**URL:** <https://community.geodynamics.org/t/density-profile-in-latent-heat-model/810>\
**Category:** ASPECT\
**Created:** [September 17, 2019, 3:55pm UTC](https://community.geodynamics.org/t/density-profile-in-latent-heat-model/810 "2019-09-17T15:55:50Z")\
**Posts on this page:** 5\
**Page:** 1

<div class="post-metadata">

**Author:** ![FELIPE](https://avatars.discourse-cdn.com/v4/letter/f/90db22/32.png) [@FELIPE](https://community.geodynamics.org/u/FELIPE)\
**Post date:** [September 17, 2019, 3:55pm UTC](https://community.geodynamics.org/t/density-profile-in-latent-heat-model/810/1 "2019-09-17T15:55:50Z")

</div>

Hi there,

I am testing some stuff with the ‘Latent heat’ material model. I have found some curious stuff that escapes a bit my understanding on what ASPECT is doing.

I have implemented a layered system with one single phase-change interface. My model is not compressible at least in the definition of the parameter file, no adiabatic heating, etc.

Upon plotting results in Paraview, carefully inspecting (plotting in different ways, and using ‘plot over line’) I realized that the density is increasing with depth (not solely jumping at the interface, but continuously inside the layers). Checking the manual, I realized that there is a default non-zero compressibility(~ +5e-12) I was unaware of.

```
I think there is nothing wrong with this in principle (physically), but I was trying to isolate everything and my maths did not give me my expected answer. I guess this compressibility is unavoidable in order to allow the software to perform the density jump across the phase change..

```

Can this compressibility set equal to 0 and have the system run correctly? I am trying to save some SUs on my account, and some explanations would be helpful.

cheers many,

Felipe

---

<div class="post-metadata">

**Author:** ![jbnaliboff](https://avatars.discourse-cdn.com/v4/letter/j/58956e/32.png) [@jbnaliboff](https://community.geodynamics.org/u/jbnaliboff)\
**Post date:** [September 17, 2019, 6:29pm UTC](https://community.geodynamics.org/t/density-profile-in-latent-heat-model/810/2 "2019-09-17T18:29:02Z")

</div>

Hi Felipe,

In the latent heat material model there are parameters that control the width of the phase transition (e.g., `Phase transition pressure widths` or `Phase transition widths`).

Can you let us know what you assigned for these values and also the other parameters related to the phase transitions?

Thanks!

Cheers,  
John

---

<div class="post-metadata">

**Author:** ![FELIPE](https://avatars.discourse-cdn.com/v4/letter/f/90db22/32.png) [@FELIPE](https://community.geodynamics.org/u/FELIPE)\
**Post date:** [September 18, 2019, 2:23pm UTC](https://community.geodynamics.org/t/density-profile-in-latent-heat-model/810/3 "2019-09-18T14:23:50Z")

</div>

Hi John,

Thanks for replying.

As I explained, although for my model I did not specify any compressibility nor adiabatic heating, the model is compressible (by default of Aspect settings), something I did not aimed for.

I have implemented 4 models of this kind. In all of them I set the same  
Clapeyron slope(-2.6 MPa/K) and the same density jump (350 kg/m3) (it’s  
like at the 660 km depth). In 2 of them I define the transition by pressure (2.904e10 Pa), with non-zero pressure widths (2.2e9 Pa); and in  
the 2 last ones by depth (660e3 m), with 0 width.

In all  
these models, the density jumps at the phase transition (by ~350, as expected), but it varies continuously inside the layers, i guess obviously in response to the default compressibility (5e-12).

cheers,

from Beijing,

Felipe

---

<div class="post-metadata">

**Author:** ![jbnaliboff](https://avatars.discourse-cdn.com/v4/letter/j/58956e/32.png) [@jbnaliboff](https://community.geodynamics.org/u/jbnaliboff)\
**Post date:** [September 18, 2019, 5:34pm UTC](https://community.geodynamics.org/t/density-profile-in-latent-heat-model/810/4 "2019-09-18T17:34:50Z")

</div>

Hi Felipe,

Thanks for the follow-up explanation. Just so we are on the same page, you are seeing density variations that would not be caused by temperature alone within each layer (e.g., assuming the thermal expansivity is not zero)?

If you are seeing density variations within the layers not explained by temperature variations alone, it is likely the default compressibility. Setting that value to zero will indeed turn off compressibility (e.g., material model returns that it is an incompressible system).

I suggest running a few low-resolution 2D runs (see test suite) to make sure the chosen parameters are giving exactly what you expect.

Please let us know if something is still amiss after running the tests.

Cheers,  
John

---

<div class="post-metadata">

**Author:** ![FELIPE](https://avatars.discourse-cdn.com/v4/letter/f/90db22/32.png) [@FELIPE](https://community.geodynamics.org/u/FELIPE)\
**Post date:** [September 19, 2019, 3:18pm UTC](https://community.geodynamics.org/t/density-profile-in-latent-heat-model/810/5 "2019-09-19T15:18:05Z")

</div>

Hi John,

Thanks a lot for replying.

I set up those models isothermic except for BCs. So, I expected the mass density, at the initial condition, to be constant inside layers…

So yes, it has to be the compressibility.

I will proceed with preparing a more complex model now that I understand how things are with the latent heat model.

cheers many,

from Beijing,

Felipe