News:

alphaMELTS 2.3 standalone & for MATLAB/Python is now open source and available on GitHub (https://github.com/magmasource/alphaMELTS).
alphaMELTS 1.9 is available at the legacy download and information site.
For news of all MELTS software see the MELTS Facebook page.

Main Menu

Magmatic mixing, isenthalpic process, operation failed.

Started by GIGN, March 28, 2025, 06:52:14 AM

Previous topic - Next topic

GIGN

Hello, I want to simulate the enrichment process of a trace element in minerals after magma mixing using Alphamelts 1.9 on Windows. The scenario involves mixing a 1200°C magma with a 1180°C magma in a 0.3:0.7 ratio. Both magmas have identical compositions, including H2O and trace elements. However, during the mixing calculation, I consistently encounter errors related to the reference enthalpy settings. Even after adjusting the enthalpy values, the simulation fails. I've tried debugging parameters like temperature, pressure, and oxygen fugacity (fO2) but the problem persists. I'm at my wit's end and urgently need assistance. Below are the steps I followed:
env
ALPHAMELTS_VERSION MELTS
ALPHAMELTS_MODE isenthalpic
ALPHAMELTS_CELSIUS_OUTPUT true
ALPHAMELTS_FLUX_MELTING true
ALPHAMELTS ASSIMILATE=0.3
ALPHAMELTS_ASSIMILATE_FILE=Enrichment.melts
ALPHAMELTS_FRACTIONATE_SOLIDS true
ALPHAMELTS_DO_TRACE true
ALPHAMELTS_DO_TRACE H2O true
ALPHAMELTS_TRACE_INPUT_FILE new_trace_data.txt
ALPHAMELTS_SAVE_ALL true


Initial.melts:
Initial Composition: SiO2 45.4076
Initial Composition: TiO2 4.10505
Initial Composition: Al2O3 15.9333
Initial Composition: Cr2O3 0.00
Initial Composition: FeO 10.6709
Initial Composition: Fe2O3 2.14495
Initial Composition: MnO 0.2072
Initial Composition: MgO 5.5787
Initial Composition: CaO 9.62144
Initial Composition: Na2O 3.59729
Initial Composition: K2O 1.26834
Initial Composition: P2O5 0.918457
Initial Composition: H2O 0.5467
Initial Temperature: 1130.00
Initial Pressure: 1000.00
Initial Trace: Sc 26.6834
Initial Trace: Ni 1.72881
log fo2 Path: FMQ

Enrichment.melts
Initial Composition: SiO2 41.8471
Initial Composition: TiO2 5.23912
Initial Composition: Al2O3 14.5134
Initial Composition: Cr2O3 0.008769
Initial Composition: FeO 14.8942
Initial Composition: Fe2O3 2.1862
Initial Composition: MnO 0.172914
Initial Composition: MgO 6.81381
Initial Composition: CaO 9.53249
Initial Composition: Na2O 3.0196
Initial Composition: K2O 1.0099
Initial Composition: P2O5 0.762578
Initial Composition: H2O 0.0
Initial Temperature: 1200.00
Initial Pressure: 1000.0
Initial Trace: Sc 24.7462
Initial Trace: Ni 3.52352
log fo2 Path: FMQ

-f-----1------3-----1 Superliquidus ------12----1  Mix traces and majors----Mass proportion 0.3

Old reference Enthalpy (H0): not set yet
Old total Enthalpy: -1.18037e+006
and over
..............................

I have tried manually inputting the enthalpy values and other steps, but all attempts have failed,I really need your help.Looking forward to your suggestions!

Paula

Check the ALPHAMELTS_ASSIMILATE entry in the documentation.
  • Note that ALPHAMELTS ASSIMILATE [sic] is missing an underscore in your file and should be to a value like 'true', not the proportion you want.
  • '=' is not a valid character in the environment file. 
  • There is no 'ALPHAMELTS_ASSIMILATE_FILE' variable. You will be asked for the file name when you call option 4.
  • If you are doing isenthalpic calculations you need to mix majors and traces. This means having ALPHAMELTS_ASSIMILATE set and leave ALPHAMELTS_FLUX_MELTING unset (with a '!' in front).
  • Menu option 12 is the Source Mixer for one off mixing. You want the Remixer, which will be triggered when you call menu option 4 if you have ALPHAMELTS_ASSIMILATE set correctly (and ALPHAMELTS_FLUX_MELTING unset).

GIGN

Thank you very much for your prompt assistance. I sincerely appreciate it.
I'm conducting an isenthalpic calculation: adding a small amount of new magma (new_magma.melts) to a shallow magma chamber (Initial.melts).(A small amount of magma from the deep magma chamber has replenished the shallow magma chamber.)
Following your advice, I gave it a try:
env.txt
ALPHAMELTS_VERSION MELTS
ALPHAMELTS_MODE isenthalpic
ALPHAMELTS_CELSIUS_OUTPUT true
ALPHAMELTS_ASSIMILATE true
ALPHAMELTS_FRACTIONATE_SOLIDS true
ALPHAMELTS_DO_TRACE true
ALPHAMELTS_DO_TRACE H2O true
ALPHAMELTS_TRACE_INPUT_FILE new_trace_data.txt
ALPHAMELTS_SAVE_ALL true


Steps:
-f -----env.txt
1------Initial.melts
3-----1--Superliquidus
12----new_magma.melts-------1--Mix traces and majors
                       20----Mass ...... added     ###Is that correct?
          Tip:  Old reference Enthalpy (H0): not set yet
                   Old total Enthalpy: -1.18037e+006
4------1--Superliquidus
0------Number of assimilant MELTS files (or '0' for binary input file)
1------Mass of enriching agent to be added in grams per cycle1
           Assimilant filename        ###I don't understand here.

       ###Do I need files for many different minerals here? I saw the discussion on "isenthalpic AFC" in the forum(2012). I don't have these files. What should I do?
After completing the isenthalpic mixing, I want to continue by allowing the mixed magma to undergo fractional crystallization under isobaric conditions. Do you have any suggestions?  Unfortunately, I'm stuck at the stage of isenthalpic mixing.My idea is to use the output from the successful mixing as the initial composition for a fractional crystallization simulation under isobaric conditions.

I have two more questions:
(1) Can Alphamelts simulate magnetite crystallization? The results I obtained were rhm-oxide. Is this ilmenite? If Alphamelts can simulate magnetite, how should I set it up?

(2) Can only integer oxygen fugacities be simulated? For example, +1FMQ, +2FMQ, but not +0.5FMQ? When I used non-integer oxygen fugacities, only liquid_0 and olivine_0 appeared, and no other minerals.

Thank you so much for your time and effort. I extend my best wishes to you!


Paula

The Assimilant filename is Enrichment.melts.

As described in the ALPHAMELTS_ASSIMILATE documentation you can have a single file or multiple files. That forum thread was about adding the option to have multiple files, but you can still have just 1 file.

See: https://magmasource.caltech.edu/alphamelts/1/manual_text.html#0.1__Toc192218456  (that link is from an older version but it's not changed significantly in the v1.9 pdf).

Magnetite is part of the spinel solid solution, and ilmenite is part of the run-oxide solid solution (see: https://magmasource.caltech.edu/forum/index.php/topic,85.0.html). If you are not getting spinel crystallizing then try suppressing rhm-oxide or change the fO2.

Check the MELTS file section of the documentation (https://magmasource.caltech.edu/alphamelts/1/alphamelts_manual.pdf) for the fO2 format. In alphaMELTS 1.9 you can specify non-integer offsets with "Log fO2 Delta:". Note that in alphaMELTS 2 that line is replaced by "Log fO2 Offset:" for more consistency other MELTS software.