| /libCEED/examples/solids/include/ |
| H A D | matops.h | 37 …ComputeStrainEnergy(DM dm_energy, UserMult user, CeedOperator op_energy, Vec X, PetscReal *energy);
40 PetscErrorCode RegressionTests_solids(AppCtx app_ctx, PetscReal energy);
|
| H A D | misc.h | 41 PetscErrorCode RegressionTests_solids(AppCtx app_ctx, PetscReal energy);
|
| H A D | structs.h | 131 CeedQFunctionUser setup_geo, residual, jacobian, energy, diagnostic, true_soln; member
|
| /libCEED/examples/solids/tests-output/ |
| H A D | MR-strain1.csv | 1 increment,energy
|
| H A D | MR-strain.csv | 1 increment,energy
|
| H A D | NH-strain.csv | 1 increment,energy
|
| /libCEED/examples/solids/src/ |
| H A D | matops.c | 219 …ComputeStrainEnergy(DM dmEnergy, UserMult user, CeedOperator op_energy, Vec X, PetscReal *energy) { in ComputeStrainEnergy() argument
254 (*energy) = 0; in ComputeStrainEnergy()
255 for (CeedInt i = 0; i < length; i++) (*energy) += e[i]; in ComputeStrainEnergy()
259 PetscCall(MPI_Allreduce(MPI_IN_PLACE, energy, 1, MPIU_REAL, MPIU_SUM, user->comm)); in ComputeStrainEnergy()
|
| H A D | misc.c | 237 PetscErrorCode RegressionTests_solids(AppCtx app_ctx, PetscReal energy) { in RegressionTests_solids() argument
242 PetscReal error = PetscAbsReal(energy - energy_ref) / energy_ref; in RegressionTests_solids()
245 …WORLD, "Energy %e does not match expected energy %e: relative tolerance %e > %e\n", (double)energy, in RegressionTests_solids()
|
| H A D | setup-libceed.c | 426 CeedQFunctionCreateInterior(ceed, 1, problem_data.energy, problem_data.energy_loc, &qf_energy); in SetupLibceedFineLevel()
|
| /libCEED/examples/solids/ |
| H A D | elasticity.c | 638 CeedScalar energy; in main() local
639 … PetscCall(ComputeStrainEnergy(dm_energy, res_ctx, ceed_data[fine_level]->op_energy, U, &energy)); in main()
643 PetscCall(PetscPrintf(comm, " Strain Energy : %.12e\n", energy)); in main()
645 PetscCall(PetscViewerASCIIPrintf(app_ctx->energy_viewer, "%f,%e\n", load_increment, energy)); in main()
773 PetscReal energy; in main() local
774 PetscCall(ComputeStrainEnergy(dm_energy, res_ctx, ceed_data[fine_level]->op_energy, U, &energy)); in main()
777 PetscCall(PetscPrintf(comm, " Strain Energy : %.12e\n", energy)); in main()
779 PetscCall(RegressionTests_solids(app_ctx, energy)); in main()
|
| H A D | index.md | 83 Here, we define a strain energy density functional $\Phi(\bm \epsilon) \in \mathbb R$ and obtain th…
87 $$ (strain-energy-grad)
90 The strain energy density functional cannot be an arbitrary function $\Phi(\bm \epsilon)$; it can o…
99 For the linear elasticity model, the strain energy density is given by
144 However, the strain energy density differs and is given by
268 Recall that the strain energy density functional can only depend upon invariants.
280 $$ (neo-hookean-energy)
284 To evaluate {eq}`strain-energy-grad`, we make use of
291 Carrying through the differentiation {eq}`strain-energy-grad` for the model {eq}`neo-hookean-energy…
322 A coupled Mooney-Rivlin strain energy density (cf. Neo-Hookean {eq}`neo-hookean-energy`) is {cite}`…
[all …]
|
| H A D | README.md | 225 …name{trace} \bm{E}$, $\operatorname{trace} \bm{E}^2$, $\lvert J \rvert$, and strain energy density.
|
| /libCEED/examples/solids/problems/ |
| H A D | linear.c | 30 .energy = ElasEnergy_Linear,
|
| H A D | finite-strain-mooney-rivlin.c | 34 .energy = ElasFSEnergy_MR,
|
| H A D | finite-strain-neo-hookean.c | 34 .energy = ElasFSEnergy_NH,
|
| /libCEED/examples/solids/qfunctions/ |
| H A D | linear.h | 224 CeedScalar(*energy) = (CeedScalar(*))out[0]; in ElasEnergy_Linear()
275 energy[i] = in ElasEnergy_Linear()
|
| H A D | finite-strain-neo-hookean.h | 403 CeedScalar(*energy) = (CeedScalar(*))out[0]; in ElasFSEnergy_NH()
458 …energy[i] = (lambda * logJ * logJ / 2. - mu * logJ + mu * (E2[0][0] + E2[1][1] + E2[2][2]) / 2.) *… in ElasFSEnergy_NH()
|
| H A D | finite-strain-mooney-rivlin.h | 424 CeedScalar(*energy) = (CeedScalar(*))out[0]; in ElasFSEnergy_MR()
500 …energy[i] = (0.5 * lambda * (logJ) * (logJ) - (mu_1 + 2 * mu_2) * logJ + (mu_1 / 2.) * (I_1 - 3) +… in ElasFSEnergy_MR()
|
| /libCEED/examples/fluids/ |
| H A D | index.md | 28 …$E$ the total energy density (defined as $E = \rho e$, where $e$ is the total energy including the…
45 … \leftarrow\textrm{ momentum density}\\ \leftarrow\textrm{ energy density} \e…
259 3. energy stabilization $\tau_E$
318 {ref}`problem-advection`, the problem of the transport of energy in a uniform vector velocity field…
522 A simplified version of system {eq}`eq-ns`, only accounting for the transport of total energy, is g…
529 In this particular test case, a blob of total energy (defined by a characteristic radius $r_c$) is …
533 In this case, a uniform circular velocity field transports the blob of total energy.
534 …We have solved {eq}`eq-advection` applying zero energy density $E$, and no-flux for $\bm{u}$ on th…
538 …linear velocity field, enters the domain and transports the blob of total energy out of the domain.
671 …ip and non-penetration boundary conditions for $\bm{u}$, and no-flux for mass and energy densities.
[all …]
|
| H A D | README.md | 8 The state variables are mass density, momentum density, and energy density.
302 …nsity $\rho$ and momentum density $\rho \bm u$ constant while advecting "total energy density" $E$.
380 - Total energy of inflow wind when `-wind_type translation`
573 - Stabilization energy constant, $C_E$
|
| /libCEED/doc/papers/joss/ |
| H A D | paper.md | 149 …isson problem) or CPU and GPU systems of similar (purchase/operating and energy) cost. These tests…
161 ![Strain energy density in a twisted Neo-Hookean beam.\label{fig:Solids}](img/SolidTwistExample.jpe…
|
| /libCEED/doc/sphinx/source/ |
| H A D | releasenotes.md | 283 - {ref}`example-petsc-elasticity` example updated with strain energy computation and more flexible …
|