Skip to content
Snippets Groups Projects
Commit 1074b3de authored by Vinayak Gholap's avatar Vinayak Gholap
Browse files

Updates in weak and strong coupling test parameters for air domain

parent e7be8b54
No related branches found
No related tags found
1 merge request!386Major update: Computing average IPVar over time period, rewriting last time...
Show whitespace changes
Inline Side-by-side
dG3D/benchmarks/strongCouplingSMP/StrongCouplingEMTM.py
+ 12
− 10
View file @ 1074b3de
......@@ -12,7 +12,7 @@ import math
# material law: SMP
lawnum= 1
rho = 270. # material density (Kg/m^3)
tinitial = 273.+25. # (K)
tinitial = 273.+47. # (K)
G=299.e6 # Shear modulus (Pa)
Mu=0.29 # Poisson ratio
alpha = beta=gamma=0. # parameter of anisotropy
......@@ -58,12 +58,12 @@ centralAxis_z = 0. # Unit vector along central axis of inductor coil
# material law: vacuum/free space region
lawnumvac = 3
rhovac = 1.2
Gvac=156.e0 # Shear modulus
Gvac=156.e1 # Shear modulus
Muvac = 0.35 # Poisson ratio
Evac= 2.0 * Gvac * (1.0 + Muvac) #youngs modulus
alphavac = betavac = gammavac = 0. # parameter of anisotropy
cpvac= 1012.*rhovac
Kxvac=Kyvac=Kzvac= 26.e-3 # thermal conductivity tensor components
cpvac= 1.e-12 #1012.*rhovac
  • Vinayak Gholap @vinugholap ·
    Author Developer

    Note: Cannot use cp = 0. for air in strong coupling test since solver fails to converge due to stiffness matrix.

    Edited by Vinayak Gholap
  • Please register or sign in to reply
Kxvac=Kyvac=Kzvac= 0.0 #26.e-3 # thermal conductivity tensor components
lxvac=lyvac=lzvac= 1.e-12 # electrical conductivity
seebeckvac= 0.
alphaThermvac= 0.0 # thermal dilatation
......@@ -73,6 +73,7 @@ mag_mu_x_vac = mag_mu_y_vac = mag_mu_z_vac = mag_r_vac * mag_mu_0 # mag permeabi
magpot0_x_vac = magpot0_y_vac = magpot0_z_vac = 0.0 # initial magnetic potential
useFluxT=True
evaluateCurlField = True;
# geometry
meshfile="undeformedSMPFullDomain.msh" # name of mesh file
......@@ -81,10 +82,10 @@ meshfile="undeformedSMPFullDomain.msh" # name of mesh file
sol = 2 # Gmm=0 (default) Taucs=1 PETsc=2
soltype = 1 # StaticLinear=0 (default) StaticNonLinear=1
#EMncycles = 5; # num of sinusoidal cycles of EM frequency application
nstep = 8000 # 8*EMncycles # number of step (used only if soltype=1)
nstep = 80 # 8*EMncycles # number of step (used only if soltype=1)
ftime = 5.e-3 # EMncycles/freq; # according to characteristic time
#EMftime = EMncycles/freq # Final time (used only if soltype=1)
tol=1.e-3 # relative tolerance for NR scheme (used only if soltype=1)
tol=1.e-6 # relative tolerance for NR scheme (used only if soltype=1)
nstepArch=1 # Number of step between 2 archiving (used only if soltype=1)
fullDg = 0 #O = CG, 1 = DG
beta1 = 1000.
......@@ -96,7 +97,7 @@ eqRatio =1.e8
# material law for smp
lawsmp = LinearElecMagTherMechDG3DMaterialLaw(lawnum,rho,E_matrix,E_matrix,E_matrix,Mu,Mu,Mu,G,G,G,alpha,beta,gamma,
tinitial,Kx,Ky,Kz,alphaTherm,alphaTherm,alphaTherm,lx,ly,lz,seebeck,cp,v0, mag_mu_x, mag_mu_y, mag_mu_z, magpot0_x,
magpot0_y, magpot0_z, Irms, freq, nTurnsCoil, coilLx, coilLy, coilLz, coilW,useFluxT)
magpot0_y, magpot0_z, Irms, freq, nTurnsCoil, coilLx, coilLy, coilLz, coilW,useFluxT,evaluateCurlField)
#lawsmp.setUseBarF(True)
......@@ -104,14 +105,14 @@ magpot0_y, magpot0_z, Irms, freq, nTurnsCoil, coilLx, coilLy, coilLz, coilW,useF
lawind = LinearElecMagInductorDG3DMaterialLaw(lawnumind,rhoind,Eind,Eind,Eind,Muind,Muind,Muind,Gind,Gind,Gind,alphaind,betaind,gammaind,
tinitial,Kxind,Kyind,Kzind,alphaThermind,alphaThermind,alphaThermind,lxind,lyind,lzind,seebeckind,cpind,v0ind, mag_mu_x_ind, mag_mu_y_ind, mag_mu_z_ind,
magpot0_x_ind,magpot0_y_ind, magpot0_z_ind, Irms, freq, nTurnsCoil, coilLx, coilLy, coilLz, coilW,
Centroid_x,Centroid_y,Centroid_z,centralAxis_x,centralAxis_y,centralAxis_z,useFluxT)
Centroid_x,Centroid_y,Centroid_z,centralAxis_x,centralAxis_y,centralAxis_z,useFluxT,evaluateCurlField)
#lawind.setUseBarF(True)
#material law for free space (vacuum)
lawvac = LinearElecMagTherMechDG3DMaterialLaw(lawnumvac,rhovac,Evac,Evac,Evac,Muvac,Muvac,Muvac,Gvac,Gvac,Gvac,alphavac,betavac,gammavac,tinitial,
Kxvac,Kyvac,Kzvac,alphaThermvac,alphaThermvac,alphaThermvac,lxvac,lyvac,lzvac,seebeckvac,cpvac,v0vac, mag_mu_x_vac, mag_mu_y_vac, mag_mu_z_vac, magpot0_x_vac,
magpot0_y_vac, magpot0_z_vac, Irms, freq, nTurnsCoil, coilLx, coilLy, coilLz, coilW,useFluxT)
magpot0_y_vac, magpot0_z_vac, Irms, freq, nTurnsCoil, coilLx, coilLy, coilLz, coilW,useFluxT,evaluateCurlField)
#lawvac.setUseBarF(True)
......@@ -168,7 +169,7 @@ mysolver.snlData(nstep,ftime,tol,1.e-9)#forsmp
mysolver.options("-ksp_type preonly")
mysolver.options("-pc_type lu")
mysolver.stepBetweenArchiving(nstepArch)
mysolver.snlManageTimeStep(25, 3, 2, 10)
mysolver.snlManageTimeStep(25, 10, 2, 10)
# BC
finalDisp = 0.01;
......@@ -250,6 +251,7 @@ mysolver.internalPointBuildView("js0_z",IPField.INDUCTORSOURCEVECTORFIELD_Z, 1,
mysolver.archivingIPOnPhysicalGroup("Volume",SMPfield, IPField.TEMPERATURE,IPField.MEAN_VALUE,nstepArch);
mysolver.archivingIPOnPhysicalGroup("Volume",SMPfield, IPField.VOLTAGE,IPField.MEAN_VALUE,nstepArch);
mysolver.archivingIPOnPhysicalGroup("Volume",SMPfield, IPField.EMFIELDSOURCE,IPField.MEAN_VALUE,nstepArch);
mysolver.archivingForceOnPhysicalGroup("Face", SurfSmp, 3, 1);
mysolver.archivingForceOnPhysicalGroup("Face", SurfSmp, 4, 1);
......
dG3D/benchmarks/weakCouplingSMP/EM_weakCoupling.py
+ 5
− 4
View file @ 1074b3de
......@@ -15,7 +15,7 @@ import pickle
# material law: SMP
lawnum= 1
rho = 270. # material density (Kg/m^3)
tinitial = 273.+25. # (K)
tinitial = 273.+47. # (K)
G=299.e6 # Shear modulus (Pa)
Mu=0.29 # Poisson ratio
alpha = beta=gamma=0. # parameter of anisotropy
......@@ -65,8 +65,8 @@ Gvac=156.e1 # Shear modulus (1.e-5 less than SMP)
Muvac = 0.35 # Poisson ratio (same as SMP)
Evac= 2.0 * Gvac * (1.0 + Muvac) #youngs modulus
alphavac = betavac = gammavac = 0. # parameter of anisotropy
cpvac= 1012.*rhovac
Kxvac=Kyvac=Kzvac= 26.e-3 # thermal conductivity tensor components
cpvac= 0.0 #1012.*rhovac
  • Vinayak Gholap @vinugholap ·
    Author Developer

    Note: Here since in weak coupling -> EM problem we don't solve for T we can set cp = 0. for air.

  • Please register or sign in to reply
Kxvac=Kyvac=Kzvac= 0.0 #26.e-3 # thermal conductivity tensor components
lxvac=lyvac=lzvac= 1.e-12 # electrical conductivity
seebeckvac= 0.
alphaThermvac= 0.0 # thermal dilatation
......@@ -100,7 +100,7 @@ soltype = 1 # StaticLinear=0 (default) StaticNonLinear=1
EMftime =EMncycles/freq # Final time (used only if soltype=1)
tol=1.e-3 # relative tolerance for NR scheme (used only if soltype=1)
tol=1.e-6 # relative tolerance for NR scheme (used only if soltype=1)
nstepArch=1 # Number of step between 2 archiving (used only if soltype=1)
fullDg = 0 #O = CG, 1 = DG
beta1 = 1000.
......@@ -262,6 +262,7 @@ mysolver.internalPointBuildView("eps_xz",IPField.STRAIN_XZ, 1, 1)
mysolver.archivingIPOnPhysicalGroup("Volume",SMPfield, IPField.TEMPERATURE,IPField.MEAN_VALUE,nstepArch);
mysolver.archivingIPOnPhysicalGroup("Volume",SMPfield, IPField.VOLTAGE,IPField.MEAN_VALUE,nstepArch);
mysolver.archivingIPOnPhysicalGroup("Volume",SMPfield, IPField.EMFIELDSOURCE,IPField.MEAN_VALUE,nstepArch);
mysolver.archivingForceOnPhysicalGroup("Face", SurfSmp, 3, 1);
mysolver.archivingForceOnPhysicalGroup("Face", SurfSmp, 4, 1);
......
dG3D/benchmarks/weakCouplingSMP/TM_weakCoupling.py
+ 3
− 2
View file @ 1074b3de
......@@ -15,7 +15,7 @@ import pickle
# material law: SMP
lawnum= 1
rho = 270. # material density (Kg/m^3)
tinitial = 273.+25. # (K)
tinitial = 273.+47. # (K)
G=299.e6 # Shear modulus (Pa)
Mu=0.29 # Poisson ratio
alpha = beta=gamma=0. # parameter of anisotropy
......@@ -59,7 +59,7 @@ print('Nb coupling: ', nbCoupling)
print('TM solver mesh file: ', meshfile)
print('Final disp: ', finalDisp)
tol=1.e-3 # relative tolerance for NR scheme (used only if soltype=1)
tol=1.e-6 # relative tolerance for NR scheme (used only if soltype=1)
nstepArch=1 # Number of step between 2 archiving (used only if soltype=1)
fullDg = 0 #O = CG, 1 = DG
beta1 = 1000.
......@@ -190,6 +190,7 @@ mysolver.internalPointBuildView("jez",IPField.ELECTRICALFLUX_Z, 1, 1)
mysolver.archivingIPOnPhysicalGroup("Volume",SMPfield, IPField.TEMPERATURE,IPField.MEAN_VALUE,nstepArch);
mysolver.archivingIPOnPhysicalGroup("Volume",SMPfield, IPField.VOLTAGE,IPField.MEAN_VALUE,nstepArch);
mysolver.archivingIPOnPhysicalGroup("Volume",SMPfield, IPField.EMFIELDSOURCE,IPField.MEAN_VALUE,nstepArch);
mysolver.archivingForceOnPhysicalGroup("Face", SurfSmp, 3, 1);
mysolver.archivingForceOnPhysicalGroup("Face", SurfSmp, 4, 1);
......
dG3D/benchmarks/weakCouplingSMP/run.py
+ 2
− 2
View file @ 1074b3de
......@@ -29,11 +29,11 @@ deltaDisp = 0.01;
finalDisp = 0.01;
# num of sinusoidal cycles of EM frequency application
EMncycles = 10;
EMncycles = 1;
# num of steps used in TM and EM solver
TMnstep = 8;
EMnstep = 8*EMncycles;
EMnstep = 80; #8*EMncycles;
stepIter = 0;
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment