Welcome to the forums. Please post in English or French.

You are not logged in.

#1 2019-11-20 08:21:48

Ashok
Member
Registered: 2019-11-20
Posts: 3

Thermal non line (Iso Mat) with material phase change ?

I am doing the Thermal non line Analysis (Isotropic Material) with material phase change. IN between analysis there is a temp, where my material state changes from POWDER to SOLID. I have attached its graph herewith. Can we able to do such Thermal non line Analysis (including Phase Shift) in code_aster. Your help will be highly appreciated.


Attachments:
LAMBDA.jpg, Size: 45.68 KiB, Downloads: 20

Offline

#2 2019-11-20 10:18:20

Christophe Durand
Nuclear Fuel Division
From: EDF
Registered: 2007-11-22
Posts: 145
Website

Re: Thermal non line (Iso Mat) with material phase change ?

Indeed, Code_Aster can deal with phase transition in thermics. You should take an example on the test case TTNL02 - Transitoire thermique avec changement de phase.

Offline

#3 2019-11-20 15:08:25

Ashok
Member
Registered: 2019-11-20
Posts: 3

Re: Thermal non line (Iso Mat) with material phase change ?

Thanks for the reply.

I have some specific requirement. Can I able to use attached python program in which depending on State of the POWDER or  SOLID I will change the conductivity?

I think the UMAT and MFRONT is not supported for Thermal non linear.

---------- OR ---------------------

I would like to define a State variable at node like NUT1 (= 0 initial value) and depending on the max temp (like above 1336 degree) my powder will change its phase and so NUT1 become 1.0 and depending on these NUT1 value I would like to change my conductivity values.

(Its like handling three values at node)
Temp  K          NUT1
25    9.5175    1
50    10.015    1
75    10.5125    1
100    11.01    1
125    11.5075    1


Temp  K           NUT1
25    0.95175    0
50    1.0015    0
75    1.05125    0
100    1.101    0
125    1.15075    0
150    1.2005    0

Is it possible in code_aster?
Any help is highly appreciated.
Thanks


Attachments:
umatht_1.py, Size: 2.22 KiB, Downloads: 8

Offline

#4 2019-11-21 15:22:06

Christophe Durand
Nuclear Fuel Division
From: EDF
Registered: 2007-11-22
Posts: 145
Website

Re: Thermal non line (Iso Mat) with material phase change ?

On my opinion, not feasable in the actual version of Code_Aster without a specific development. The state of the matter (powder or solid) should be a specific internal variable to be added.

Thank you for your feedback. If this development interests a motivated developer, this feature could integrate a future version.

Offline

#5 2019-12-02 06:53:06

Ashok
Member
Registered: 2019-11-20
Posts: 3

Re: Thermal non line (Iso Mat) with material phase change ?

Hello All

I still tried this case, I am getting the error.
Can any one has any solution for it?

TEMPMAX4 = CREA_CHAMP(  RESULTAT = SOL_HEAT,
                       OPERATION = 'EXTR',
               NOM_CHAM = 'TEMP',
               TYPE_CHAM= 'NOEU_TEMP_R',
               TYPE_MAXI = 'MAXI',
               TYPE_RESU ='VALE',
               );

              
INCLUDE (UNITE=90, INFO= 2)              
STATV = FORMULE(NOM_PARA=('TEMP'),VALE='''findSTV(TEMP)''',      # Create a function of SateV depending on Temp
                 );              
              

FLD_STV=CREA_CHAMP(TYPE_CHAM='NOEU_NEUT_F',    # Create a field of functions (STATV) per node
                   OPERATION='AFFE',
                   MAILLAGE=WELDMESH,
                   AFFE=_F(
                   TOUT='OUI',
                   NOM_CMP='X1',
                   VALE_F=STATV,),);
                  

CHYOUNG=CREA_CHAMP(TYPE_CHAM='NOEU_NEUT_R',    # Evaluate the field FLD_STV with parameter field TEMPMAX4 at each node
                   OPERATION='EVAL',
                   CHAM_F=FLD_STV,
                   CHAM_PARA=TEMPMAX4,);

# --------------------------------------------------------------------------------------                  
              
IMPR_RESU(RESU=_F(RESULTAT=SOL_HEAT, LIST_INST=OUT_TIME,), FORMAT='MED',); 

R_MxTEMP =CREA_RESU(OPERATION='AFFE',
                  TYPE_RESU='FOURIER_THER',
                   NOM_CHAM='TEMP',
                   AFFE=_F(CHAM_GD=TEMPMAX4,),
                           );

 
FLD_STVR =CREA_RESU(OPERATION='AFFE',
                  TYPE_RESU='EVOL_ELAS',
                   NOM_CHAM='NEUT',
                   AFFE=_F(CHAM_GD=FLD_STV, INST = 1),
                           );
                          
IMPR_RESU(FORMAT='MED',
          RESU=_F(RESULTAT=FLD_STVR,));

I think grabbing the X1 component form NEUT is the problem.
Your valuable guidance will be greatly appriciacted.

Ashok


Attachments:
Componenet_Dev_error.zip, Size: 14.78 KiB, Downloads: 5

Offline