Commit b554c55e authored by Benoit Urruty's avatar Benoit Urruty

treatment_inversion

parent 02772f06
This diff is collapsed.
......@@ -50,7 +50,7 @@ $$
$$
Mean Viscosity Derivative :
$$
\partial\eta=2\times prefactor \times\mu
\partial\eta=2\times prefactor \times\muth
$$
Friction Parameter :
$$
......@@ -58,9 +58,9 @@ $$
$$
Friction parameter derivative :
$$
\frac{\partial J}{\partial\alpha}=\frac{\partial J}{\partial\beta}\times 10 ^ {\alpha}\times ln 10
\frac{\partial }{\partial\alpha}=\frac{\partial J}{\partial\beta}\times 10 ^ {\alpha}\times ln 10
$$
Grounded area coefficient : give the value 1 if it's grounded or 0 if it the grounding line
Grounded area coefficient : give the value 1 if it's grounded or 0 if it’s at the grounding line
Flux:
$$
......@@ -72,11 +72,11 @@ $$
The SSA is computed with the first guess on the value $\eta$ and $\beta$. We solve the equation :
$$
\frac{\partial}{\partial x}(2\eta h (2 \frac{\partial u}{\partial x} + \eta \frac{\partial v}{\partial y})) + \frac{\partial}{\partial y}(\eta h (\frac{\partial u}{\partial y} + \eta \frac{\partial v}{\partial x})) - \beta^2u =\rho_igh\frac{\partial s}{\partial x}
\frac{\partial}{\partial x}(2\eta h (2 \frac{\partial u}{\partial x} + \frac{\partial v}{\partial y})) + \frac{\partial}{\partial y}(\eta h (\frac{\partial u}{\partial y} + \frac{\partial v}{\partial x})) - \beta^2u =\rho_igh\frac{\partial s}{\partial x}
$$
$$
\frac{\partial}{\partial y}(2\eta h (2 \frac{\partial v}{\partial y} + \eta \frac{\partial u}{\partial x})) + \frac{\partial}{\partial x}(\eta h (\frac{\partial v}{\partial x} + \eta \frac{\partial u}{\partial y})) - \beta^2v =\rho_igh\frac{\partial s}{\partial y}
\frac{\partial}{\partial y}(2\eta h (2 \frac{\partial v}{\partial y} + \frac{\partial u}{\partial x})) + \frac{\partial}{\partial x}(\eta h (\frac{\partial v}{\partial x} + \frac{\partial u}{\partial y})) - \beta^2v =\rho_igh\frac{\partial s}{\partial y}
$$
......@@ -112,19 +112,9 @@ $$
The solver is taking as input the variable $Velocityb$. This variable contains the sensitivity of the cost function $xb$.
### The gradient
We need to compute the gradient for each parameter we have to optimize. Theoretically, we have the right solution when the gradient is equal to zero but we never obtain this value. So we define a lower limit which defines the gradient is enough low to consider the result as the best we can obtain.
$$
\frac{\partial J}{\partial \beta}=\frac{J(i)-J(i-1)}{\beta(i)-\beta(i-1)}\\
\frac{\partial J}{\partial \eta}=\frac{J(i)-J(i-1)}{\eta(i)-\eta(i-1)}
$$
The output from this solver are the nodal derivative of the friction parameter and the mean viscosity.
......
146271 275809 16731
135233 255419 15045
2
202 16731
303 275809
202 15045
303 255419
......@@ -52,37 +52,27 @@ Body Force 1
End
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
Solver 1
Equation = Reader
Procedure = "GridDataReader" "GridDataReader"
Variable = -nooutput dumy
!---- NOTE: File is case sensitive, String is not!
Filename = File "$TOPOGRAPHY_DATA$"
X Dim Name = String "x"
Y Dim Name = String "y"
!---
Variable 1 = File "thickness"
Target Variable 1 = String "H"
Exported Variable 1 = H
End
Equation = "ScatteredInter"
Variable = -nooutput "dumy"
procedure = "Scattered2DDataInterpolator" "Scattered2DDataInterpolator"
Bounding Box dx = Real 20000.0
!!!!! NNI or linear (nn-c library)
Variable 1 = File "thickness"
Target Variable 1 = String "H"
Variable 1 data file = File "$TOPOGRAPHY_DATA"
Variable 1 method = String "li"
Variable 2 = File "norm"
Target Variable 2 = String "Vobs"
Variable 2 data file = File "$VELOCITY_DATA"
Variable 2 method = String "li"
End
Solver 2
Equation = Reader2
Procedure = "GridDataReader" "GridDataReader"
Variable = -nooutput dumy
!---- NOTE: File is case sensitive, String is not!
Filename = File "$VELOCITY_DATA$"
X Dim Name = String "x"
Y Dim Name = String "y"
!--- Interpolation variables
Variable 1 = File "norm"
Target Variable 1 = String "Vobs"
Exported Variable 1 = "Vobs"
End
Solver 3
Equation = "Distance"
Variable = distance
Procedure = "DistanceSolve" "DistanceSolver1"
......@@ -92,7 +82,7 @@ End
! Compute the metric associated with f1
! Compute the metric associated with f1
! 1- compute projected gradient of f1
Solver 4
Solver 3
Equation = "Nodal Gradient 1"
Variable = -nooutput "Gradient1"
Variable DOFs = 2
......@@ -103,7 +93,7 @@ End
! 2- compute:
! - the hessian matrix by solving a diffusion equation:
! - the metric tensor
Solver 5
Solver 4
Equation = "Metric1"
Procedure = "ElmerIce_MeshAdapt2D" "MMG2D_MetricAniso"
......@@ -117,14 +107,14 @@ Solver 5
Linear System Refactorize = False
End
! Compute the metric associated with f2
Solver 6
Solver 5
Equation = "Nodal Gradient 2"
Variable = -nooutput "Gradient2"
Variable DOFs = 2
Procedure = "ElmerIce_MeshAdapt2D" "Compute2DNodalGradient"
Variable Name = string "h"
End
Solver 7
Solver 6
Equation = "Metric2"
Procedure = "ElmerIce_MeshAdapt2D" "MMG2D_MetricAniso"
......@@ -138,7 +128,7 @@ Solver 7
Linear System Refactorize = False
End
!! do the intersection of M1 and M2
Solver 8
Solver 7
Equation = "Metric"
Procedure = "ElmerIce_MeshAdapt2D" "MMG2D_MetricIntersect"
......@@ -147,7 +137,7 @@ Solver 8
Metric 2 Variable Name = String "Mu"
End
Solver 9
Solver 8
Equation = SaveScalars
Procedure = "SaveData" "SaveScalars"
Filename = "f_$RUN$.dat"
......@@ -166,7 +156,7 @@ Solver 9
End
!! Anisotropic mesh adaptation using the MMG library
Solver 10
Solver 9
!! mandatory else Model % Mesh % Changed reset to .FALSE. in coupled simulations
Exec Solver = after timestep
......@@ -188,7 +178,7 @@ End
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
Equation 1
Active Solvers(10) = 1 2 3 4 5 6 7 8 9 10
Active Solvers(10) = 1 2 3 4 5 6 7 8 9
End
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
......
......@@ -2776,3 +2776,19 @@ index de butressing (travail de johannes)
Entered on [2020-05-19 mar. 10:28]
[[file:~/Documents/PhD_TiPACCS/TiPACCS_git/journal/journal.org]]
*** 2020-05-26 mardi
**** reunion
comité de thèse!!!!!!!
MESH
changer scattered2Ddatainterpolator
retrouver fichier julien
temps des methodes faire que 10 itérations GMRES, GCR, MUMPS
eqqayé sur 2 noeuds 48 partitions si je vais 2 fois plus vite
ne pas tout faire tourner jusque qq itérations
Entered on [2020-05-26 mar. 14:58]
[[file:~/Documents/PhD_TiPACCS/TiPACCS_git/journal/journal.org::*regularisation][regularisation]]
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment