Issues: Multicomponent model

ASPECT
1

Hello.
I’m a newcomer of “ASEPCT”.
The cookbook seems good enough to start “ASPECT” but it’s quite difficult for me.
Like the jelly fish problem in Cookbook, I drew a simple subduction model in 2D using Inkscape and converted the figure into ascii values (x, y, phase).

Then, loading the file, however, an error occurs in some conditions.
Following in condition A, no error occur and the result show like as I expect.
It says two thing : the external material index file has no problem and input is not wrong, I believe.

But, following condition B and C, they are not working well.

I can’t understand. Any one has solutions?

subduction_test.txt includes x, y, phase(0, 1, 2). Phase denotes for mantle, two continental plates.

Condition A:

subsection Compositional fields
set Number of fields = 1
end

subsection Initial composition model
set Model name = ascii data
subsection Ascii data model
set Data directory = /home/sungho/Geodynamics/subduction_test/
set Data file name = subduction.txt
end
end

subsection Material model
set Model name = multicomponent
subsection Multicomponent
set Densities = 3200
set Specific heats = 1000
set Thermal conductivities = 2.5
set Thermal expansivities = 3.e-5
set Viscosities = 1.0e21
end
end

Condtion B:

subsection Compositional fields
set Number of fields = 3
end

subsection Initial composition model
set Model name = ascii data
subsection Ascii data model
set Data directory = /home/sungho/Geodynamics/subduction_test/
set Data file name = subduction.txt
end
end

subsection Material model
set Model name = multicomponent
subsection Multicomponent
set Densities = 3200
set Specific heats = 1000
set Thermal conductivities = 2.5
set Thermal expansivities = 3.e-5
set Viscosities = 1.0e21
end
end

Condition C:

subsection Compositional fields
set Number of fields = 1
end

subsection Initial composition model
set Model name = ascii data
subsection Ascii data model
set Data directory = /home/sungho/Geodynamics/subduction_test/
set Data file name = subduction.txt
end
end

subsection Material model
set Model name = multicomponent
subsection Multicomponent
set Densities = 3200, 3200, 3200
set Specific heats = 1000, 1000, 1000
set Thermal conductivities = 2.5, 2.5, 2.5
set Thermal expansivities = 3.e-5, 3.e-5, 3.e-5
set Viscosities = 1.0e21, 1.0e21, 1.0e21
end
end

2

Sungho: Can you state what error you get precisely?
Best
W.

3

In advance, I appreciate your attention to mine.

It shows quite long lines. I don’t know which are specifically related to main the error.
Please understand what I have faced difficulties.
Plus, if you want to have a look at my input files, I can email you.

Thanks

Best regards

Sungho

this is that ASPECT shows me

Screenshot%20from%202019-04-15%2010-15-19
Screenshot from 2019-04-15 10-15-19.png1670×1028 147 KB

4

Are you specifying the field “Densities” in the material model? If yes, how? This is what the error message is saying.

5

I assign field names from the external file(subduction.txt).
Then, I specify field densities using multicomponent model like below way.

subsection Compositional fields
set Number of fields = 1
end

subsection Initial composition model
set Model name = ascii data
subsection Ascii data model
set Data directory = /home/sungho/Geodynamics/subduction_test/
set Data file name = subduction.txt
end
end

subsection Material model
set Model name = multicomponent
subsection Multicomponent
set Densities = 3200, 3200, 3200
set Specific heats = 1000, 1000, 1000
set Thermal conductivities = 2.5, 2.5, 2.5
set Thermal expansivities = 3.e-5, 3.e-5, 3.e-5
set Viscosities = 1.0e21, 1.0e21, 1.0e21
end
end

6

The documentation of the parameter says

"List of densities for background mantle and compositional fields,"
"for a total of N+1 values, where N is the number of compositional fields."
"If only one value is given, then all use the same value.  Units: $kg / m^3$");

As you have one compositional field, you need to provide 1 or 2 values (but not 3).