You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
Copy file name to clipboardExpand all lines: src/darcy.jl
+1-1Lines changed: 1 addition & 1 deletion
Original file line number
Diff line number
Diff line change
@@ -78,7 +78,7 @@ uD = VectorValue(0.0,0.0)
78
78
U =TrialFESpace(V,uD)
79
79
P =TrialFESpace(Q)
80
80
81
-
# When the singe-field spaces have been designed, the multi-field test and trial spaces are expressed as arrays of single-field ones in a natural way.
81
+
# When the sinlge-field spaces have been designed, the multi-field test and trial spaces are expressed as arrays of single-field ones in a natural way.
Copy file name to clipboardExpand all lines: src/elasticity.jl
+14-2Lines changed: 14 additions & 2 deletions
Original file line number
Diff line number
Diff line change
@@ -177,9 +177,9 @@ end
177
177
178
178
a(u,v) =∫( ε(v) ⊙ (σ_bimat∘(ε(u),tags)) )*dΩ
179
179
180
-
# In previous line, pay attention in the usage of the new constitutive law `σ_bimat`. Note that we have passed the vector `tags` containing the material identifiers in the last argument of the function`.
180
+
# In previous line, pay attention in the usage of the new constitutive law `σ_bimat`. Note that we have passed the vector `tags` containing the material identifiers in the last argument of the function.
181
181
#
182
-
# At this point, we can build the FE problem again and solve it
182
+
# At this point, we can build the FE problem again and solve it:
# A material law depending on the model tags but not on the fields can be defined as follow:
195
+
196
+
tags_field =CellField(tags, Ω)
197
+
σ_from_tag(tag) = tag==alu_tag ?1.:0.
198
+
σ_bimat_cst = σ_from_tag ∘ tags_field
199
+
200
+
# `tags_field` is a field which value at $x$ is the tag of the cell containing $x$. `σ_bimat_cst` is used like a constant in (bi)linear form definition and solution export:
Copy file name to clipboardExpand all lines: src/fsi_tutorial.jl
+1-1Lines changed: 1 addition & 1 deletion
Original file line number
Diff line number
Diff line change
@@ -66,7 +66,7 @@
66
66
# <a name="geometry"></a>
67
67
# ```
68
68
# ### Geometry and Discrete model
69
-
# In this tutorial we solve the benchmark descrived in [1], consisting on a flow over an elastic flag after a cylinder.
69
+
# In this tutorial we solve the benchmark described in [1], consisting on a flow over an elastic flag after a cylinder.
70
70
# The computational domain is defined by a channel of size $\Omega \doteq (0,4.5)\times(0,0.41)$, with an embedded cylinder of radius $R=0.05$ and center at $C=(0.2,0.2)$. The associated FE triangulation is denoted by $\mathcal{T}$, the fluid and solid domain and their associated triangulations will be denoted by $\Omega_{\rm F}$, $\Omega_{\rm S}$, $\mathcal{T}_{\rm F}$ and $\mathcal{T}_{\rm S}$, respectively.
71
71
#
72
72
# In order to load the discrete model we first setup Gridap
# In our case, the mass term ($m(t, \cdot, \cdot)$) is constant in time. We can take advantage of that to save some computational effort, and indicate it to Gridap as follows
0 commit comments