Home
last modified time | relevance | path

Searched refs:E (Results 1 – 25 of 31) sorted by relevance

12

/libCEED/examples/solids/qfunctions/
H A Dmanufactured-force.h21 CeedScalar E; // Young's Modulus member
37 const CeedScalar E = context->E; in SetupMMSForce() local
48 …force[i + 0 * Q] = (-(E * (cos(x * 2.0) * cos(y * 3.0) * exp(z * 4.0) * 4.0 - cos(z * 4.0) * sin(y… in SetupMMSForce()
50 …(E * (cos(z * 4.0) * sin(y * 3.0) * exp(x * 2.0) * (4.5) + sin(x * 2.0) * sin(z * 4.0) * exp(y * 3… in SetupMMSForce()
52 … (E * nu * cos(x * 2.0) * cos(y * 3.0) * exp(z * 4.0) * 8.0) / ((nu * 2.0 - 1.0) * (nu + 1.0)) - in SetupMMSForce()
53 … (E * nu * sin(x * 2.0) * sin(z * 4.0) * exp(y * 3.0) * 6.0) / ((nu * 2.0 - 1.0) * (nu + 1.0)) - in SetupMMSForce()
54 …(E * cos(z * 4.0) * sin(y * 3.0) * exp(x * 2.0) * (nu - 1.0) * 4.0) / ((nu * 2.0 - 1.0) * (nu + 1.… in SetupMMSForce()
58 …force[i + 1 * Q] = (-(E * (cos(y * 3.0) * cos(z * 4.0) * exp(x * 2.0) * 3.0 - cos(x * 2.0) * sin(z… in SetupMMSForce()
60 …(E * (cos(x * 2.0) * sin(z * 4.0) * exp(y * 3.0) * 8.0 + sin(x * 2.0) * sin(y * 3.0) * exp(z * 4.0… in SetupMMSForce()
62 … (E * nu * cos(y * 3.0) * cos(z * 4.0) * exp(x * 2.0) * 6.0) / ((nu * 2.0 - 1.0) * (nu + 1.0)) - in SetupMMSForce()
[all …]
H A Dlinear.h21 CeedScalar E; // Young's Modulus member
39 const CeedScalar E = context->E; in ElasResidual_Linear() local
98 const CeedScalar ss = E / ((1 + nu) * (1 - 2 * nu)); in ElasResidual_Linear()
136 const CeedScalar E = context->E; in ElasJacobian_Linear() local
193 const CeedScalar ss = E / ((1 + nu) * (1 - 2 * nu)); in ElasJacobian_Linear()
228 const CeedScalar E = context->E; in ElasEnergy_Linear() local
232 const CeedScalar TwoMu = E / (1 + nu); in ElasEnergy_Linear()
234 const CeedScalar Kbulk = E / (3 * (1 - 2 * nu)); // Bulk Modulus in ElasEnergy_Linear()
296 const CeedScalar E = context->E; in ElasDiagnostic_Linear() local
300 const CeedScalar TwoMu = E / (1 + nu); in ElasDiagnostic_Linear()
[all …]
H A Dfinite-strain-neo-hookean.h21 CeedScalar E; // Young's Modulus member
158 const CeedScalar E = context->E; in ElasFSResidual_NH() local
160 const CeedScalar TwoMu = E / (1 + nu); in ElasFSResidual_NH()
162 const CeedScalar Kbulk = E / (3 * (1 - 2 * nu)); // Bulk Modulus in ElasFSResidual_NH()
265 const CeedScalar E = context->E; in ElasFSJacobian_NH() local
269 const CeedScalar TwoMu = E / (1 + nu); in ElasFSJacobian_NH()
271 const CeedScalar Kbulk = E / (3 * (1 - 2 * nu)); // Bulk Modulus in ElasFSJacobian_NH()
407 const CeedScalar E = context->E; in ElasFSEnergy_NH() local
409 const CeedScalar TwoMu = E / (1 + nu); in ElasFSEnergy_NH()
411 const CeedScalar Kbulk = E / (3 * (1 - 2 * nu)); // Bulk Modulus in ElasFSEnergy_NH()
[all …]
H A Dconstant-force.h21 CeedScalar E; // Young's Modulus member
/libCEED/examples/fluids/
H A DSTGRand.dat2 1.0E0 0.0E0 0.0E0 1.4E0 0.0E0 7.071067811865475E-1 7.071067811865475E…
3 0.0E0 7.071067811865475E-1 7.071067811865475E-1 2.4E0 1.0E0 0.0E0 0.0E0
H A Dconv_plot.py49 E = data['rel_error']
50 H = amin(E) * (h / amin(h))**p
51 ax.loglog(h, E, 'o', color=colors[i])
H A DREADME.md48 - `1E-11`
302 …nsity $\rho$ and momentum density $\rho \bm u$ constant while advecting "total energy density" $E$.
321 - `1E-2`
326 - `1E-2`
331 - `1E-6`
441 - `1E-2`
446 - `1E-2`
485 - `1E-2`
490 - `1E-2`
928 - `5.9E-4`
H A Dindex.md23 \frac{\partial E}{\partial t} + \nabla \cdot \left( \frac{(E + P)\bm{U}}{\rho} -\bm{u} \cdot \bm{\s…
28 …bm{u}$, where $\bm{u}$ is the vector velocity field), $E$ the total energy density (defined as $E …
31 P = \left( {c_p}/{c_v} -1\right) \left( E - {\bm{U}\cdot\bm{U}}/{(2 \rho)} \right) \, ,
45 …rix} \rho \\ \bm{U} \equiv \rho \bm{ u }\\ E \equiv \rho e …
56 {(E + P)\bm{U}}/{\rho}
216 (E + P)\diff\bm U/\rho + (\diff E + \diff P)\bm U/\rho - (E + P) \bm U/\rho^2 \diff\rho
525 \frac{\partial E}{\partial t} + \nabla \cdot (\bm{u} E ) - \kappa \nabla E = 0 \, ,
534 …We have solved {eq}`eq-advection` applying zero energy density $E$, and no-flux for $\bm{u}$ on th…
546 …For the outflow boundary conditions, we have used the current values of $E$, following {cite}`papa…
550 …_N) \cdot \widehat{\bm{n}} \,dS = \int_{\partial \Omega_{outflow}} \bm v \, E \, \bm u \cdot \wide…
[all …]
/libCEED/examples/solids/
H A Dindex.md16 …nd strain, respectively, in the small strain regime, while $\bm S$ and $\bm E$ are their finite-st…
20 {\overbrace{\bm S(\bm E)}^{\text{Finite Strain Hyperelastic}}}
22 {\overbrace{\bm S = \mathsf C \bm E}^{\text{St. Venant-Kirchoff}}} \\
23 …@V{\text{geometric}}V{\begin{smallmatrix}\bm E \to \bm \epsilon \\ \bm S \to \bm \sigma \end{small…
24 …@V{\begin{smallmatrix}\bm E \to \bm \epsilon \\ \bm S \to \bm \sigma \end{smallmatrix}}V{\text{geo…
114 \begin{aligned} \lambda &= \frac{E \nu}{(1 + \nu)(1 - 2 \nu)} \\ \mu &= \frac{E}{2(1 + \nu)} \end{a…
262 \bm E = \frac 1 2 (\bm C - \bm I_3) = \frac 1 2 \Big( \nabla_X \bm u + (\nabla_X \bm u)^T + (\nabla…
266 …, appropriate for large deformations, express $\bm S$ as a function of $\bm E$, similar to the lin…
269 We will assume without loss of generality that $\bm E$ is diagonal and take its set of eigenvalues …
270 …alternate choices, such as $\operatorname{trace}(\bm E), \operatorname{trace}(\bm E^2), \lvert \bm…
[all …]
H A DREADME.md20 ./elasticity -mesh [.exo file] -degree [degree] -nu [nu] -E [E] [boundary options] -problem [proble…
39 * - `-E [real]`
40 - [Young's modulus](https://en.wikipedia.org/wiki/Young%27s_modulus), $E > 0$
65 ./elasticity -mesh [.exo file] -degree 4 -E 1e6 -nu 0.3 -bc_clamp 998,999 -bc_clamp_998_translate 0…
76 ./elasticity -problem FS-NH -E 1 -nu 0.3 -num_steps 40 -snes_linesearch_type cp -dm_plex_box_faces …
158 ./elasticity -mesh [mesh] -degree [degree] -nu [nu] -E [E] -forcing mms
212 * - Young's modulus, $E$
219 ….g., {code}`-units_meter 100` to measure displacement in centimeters), but $E$ and $\int \rho \bm …
225 …t in the $z$ direction, pressure, $\operatorname{trace} \bm{E}$, $\operatorname{trace} \bm{E}^2$, …
244 - $\operatorname{trace} \bm{E}$
[all …]
/libCEED/examples/solids/problems/
H A Dneo-hookean.c66 phys->E = 0; in ProcessPhysics_NH()
77 PetscCall(PetscOptionsScalar("-E", "Young's Modulus", NULL, phys->E, &phys->E, &Young_flag)); in ProcessPhysics_NH()
102 phys->E *= units->Pascal; in ProcessPhysics_NH()
H A Dneo-hookean.h19 CeedScalar E; // Young's Modulus member
/libCEED/examples/fluids/qfunctions/
H A Deulervortex.h87 CeedScalar rho, P, T, E, u[3] = {0.}; in Exact_Euler() local
110 E = 2.; in Exact_Euler()
117 q[4] = E; in Exact_Euler()
121 E = 2.; in Exact_Euler()
130 q[4] = E; in Exact_Euler()
182 …bian_Euler(CeedScalar dF[3][5][5], const CeedScalar rho, const CeedScalar u[3], const CeedScalar E, in ConvectiveFluxJacobian_Euler() argument
191 …dF[i][4][k + 1] = ((i == k) ? (E * gamma / rho - (gamma - 1.) * u_sq / 2.) : 0.) - (gamma - 1.… in ConvectiveFluxJacobian_Euler()
195 dF[i][4][0] = u[i] * ((gamma - 1.) * u_sq - E * gamma / rho); in ConvectiveFluxJacobian_Euler()
285 const CeedScalar E = q[4][i]; in Euler() local
311 …alar E_kinetic = 0.5 * rho * (u[0] * u[0] + u[1] * u[1] + u[2] * u[2]), E_internal = E - E_kinetic, in Euler()
[all …]
H A Dshocktube.h114 …bian_Euler(CeedScalar dF[3][5][5], const CeedScalar rho, const CeedScalar u[3], const CeedScalar E, in ConvectiveFluxJacobian_Euler() argument
123 …dF[i][4][k + 1] = ((i == k) ? (E * gamma / rho - (gamma - 1.) * u_sq / 2.) : 0.) - (gamma - 1.… in ConvectiveFluxJacobian_Euler()
127 dF[i][4][0] = u[i] * ((gamma - 1.) * u_sq - E * gamma / rho); in ConvectiveFluxJacobian_Euler()
240 const CeedScalar E = q[4][i]; in EulerShockTube() local
268 …alar E_kinetic = 0.5 * rho * (u[0] * u[0] + u[1] * u[1] + u[2] * u[2]), E_internal = E - E_kinetic, in EulerShockTube()
291 …for (CeedInt j = 0; j < 3; j++) dv[j][4][i] += wdetJ * (E + P) * (u[0] * dXdx[j][0] + u[1] * dXdx[… in EulerShockTube()
329 ConvectiveFluxJacobian_Euler(jacob_F_conv, rho, u, E, gamma); in EulerShockTube()
H A Dchannel.h124 const CeedScalar E = rho_in * e_internal + E_kinetic; in Channel_Inflow() local
140 v[4][i] -= wdetJb * u_normal * (E + P); in Channel_Inflow()
163 const CeedScalar E = q[4][i]; in Channel_Outflow() local
180 v[4][i] -= wdetJb * u_normal * (E + P); in Channel_Outflow()
/libCEED/examples/nek/boxes/
H A Db1e.rea16 -0.500000E-02 p012 DT
25 0.100000E-05 p021 DIVERGENCE
26 0.100000E-07 p022 HELMHOLTZ
28 0.100000E-09 p024 TOLREL
29 0.100000E-09 p025 TOLABS
48 0.00000 p044 0=E-based/1=A-based prec.
107 0.500000E-01 p103 weight of stabilizing filter (.01)
/libCEED/examples/petsc/src/
H A Dmatops.c176 Vec E; in ComputeL2Error() local
179 PetscCall(VecDuplicate(X, &E)); in ComputeL2Error()
180 PetscCall(ApplyLocal_Ceed(X, E, op_error_ctx)); in ComputeL2Error()
181 PetscCall(VecViewFromOptions(E, NULL, "-error_view")); in ComputeL2Error()
183 PetscCall(VecSum(E, &error_sq)); in ComputeL2Error()
185 PetscCall(VecDestroy(&E)); in ComputeL2Error()
/libCEED/doc/sphinx/source/
H A DlibCEEDapi.md51 …us ($H^1$) elements, where we use the notions **T-vector**, **L-vector**, **E-vector** and **Q-vec…
56 - Element restriction $\bm{\mathcal{E}}$
60 …t and trial space differ, they get their own versions of $\bm{P}$, $\bm{\mathcal{E}}$ and $\bm{B}$.
68 …adaptive mesh refinement (AMR), the restrictions $\bm{P}$ and $\bm{\mathcal{E}}$ will involve not …
93 - Per element decomposition, **E-vector**:
98 > ```{image} ../../img/E-vector.svg
105 …> - this way, an **E-vector** can be derived from an **L-vector** without any communications and w…
106 …> - in other words, an entry in an **E-vector** is obtained by copying an entry from the correspon…
117 …> - this is similar to **E-vector** where instead of DoFs, the vector represents values at quadrat…
122 > - here X can be any of the T, L, E, or Q categories
[all …]
H A Dreferences.bib53 author = {Hughes, Thomas J R and Scovazzi, Guglielmo and Tezduyar, Tayfun E},
149 author = {Whiting, Christian H and Jansen, Kenneth E and Dey, Saikat},
175 author = {Tezduyar, Tayfun E and Senga, Masayoshi},
254 author = {Prakash, Aviral and Jansen, Kenneth E. and Evans, John A.},
265 author = {Prakash, Aviral and Jansen, Kenneth E. and Evans, John A.},
298 author = {Van Driest, E. R.},
H A DlibCEEDdev.md38 …In order to use operator fallback, the parent backend and fallback backend must use compatible E-v…
55 The `/cpu/self/ref/blocked` backend updates the {ref}`CeedOperator` to use an E-vector and Q-vector…
73 …E-vectors, then {ref}`CeedBasis` kernels map from the E-vectors to Q-vectors, then the {ref}`CeedQ…
91 Ceed backends are free to use any E-vector and Q-vector data layout (including never fully forming …
92 There are several common layouts for L-vectors, E-vectors, and Q-vectors, detailed below:
101 - **E-vector** layouts
104 …If the backend uses a strided E-vector layout, then the data for node `i`, component `j`, element …
105 …- Backends may choose to use a non-strided E-vector layout; however, the `t2**` tests will not fun…
112 …_EVAL_NONE` is properly ordered when the backend uses different layouts for E-vectors and Q-vector…
/libCEED/doc/sphinx/source/api/
H A Dindex.rst24 O(CeedOperator) --> E(CeedElemRestriction)
29 E --> V
/libCEED/doc/bib/
H A Dreferences.bib48 …title = {{H}igh-performance operator evaluations with ease of use: lib{C}{E}{E}{D}'s {P}ython …
/libCEED/examples/python/
H A Dtutorial-2-elemrestriction.ipynb125 …lly used for data stored at quadrature points or for vectors stored in the [E-vector](https://libc…
200 …ked element restriction (from an L-vector to an E-vector) and its transpose (inverse operation, fr…
270 …ked element restriction (from an L-vector to an E-vector) and its transpose (inverse operation, fr…
/libCEED/julia/LibCEED.jl/src/
H A DElemRestriction.jl271 Use the [`ElemRestriction`](@ref) to convert from L-vector to an E-vector (or apply the
290 Use the [`ElemRestriction`](@ref) to convert from L-vector to an E-vector (or apply the
318 Return a new [`CeedVector`](@ref) E-vector.
340 Return an (L-vector, E-vector) pair.
/libCEED/
H A D.mailmap23 Kenneth E. Jansen <Kenneth.Jansen@colorado.edu> <kenneth.jansen@colorado.edu>

12