1 // Copyright (c) 2017-2026, Lawrence Livermore National Security, LLC and other CEED contributors.
2 // All Rights Reserved. See the top-level LICENSE and NOTICE files for details.
3 //
4 // SPDX-License-Identifier: BSD-2-Clause
5 //
6 // This file is part of CEED: http://github.com/ceed
7
8 #include "../problems/neo-hookean.h"
9
10 #include <ceed.h>
11 #include <petscsys.h>
12
13 // Build libCEED context object
PhysicsContext_NH(MPI_Comm comm,Ceed ceed,Units * units,CeedQFunctionContext * ctx)14 PetscErrorCode PhysicsContext_NH(MPI_Comm comm, Ceed ceed, Units *units, CeedQFunctionContext *ctx) {
15 Physics_NH phys;
16
17 PetscFunctionBegin;
18
19 PetscCall(PetscMalloc1(1, units));
20 PetscCall(PetscMalloc1(1, &phys));
21 PetscCall(ProcessPhysics_NH(comm, phys, *units));
22 CeedQFunctionContextCreate(ceed, ctx);
23 CeedQFunctionContextSetData(*ctx, CEED_MEM_HOST, CEED_COPY_VALUES, sizeof(*phys), phys);
24 PetscCall(PetscFree(phys));
25
26 PetscFunctionReturn(PETSC_SUCCESS);
27 }
28
29 // Build libCEED smoother context object
PhysicsSmootherContext_NH(MPI_Comm comm,Ceed ceed,CeedQFunctionContext ctx,CeedQFunctionContext * ctx_smoother)30 PetscErrorCode PhysicsSmootherContext_NH(MPI_Comm comm, Ceed ceed, CeedQFunctionContext ctx, CeedQFunctionContext *ctx_smoother) {
31 PetscScalar nu_smoother = 0;
32 PetscBool nu_flag = PETSC_FALSE;
33 Physics_NH phys, phys_smoother;
34
35 PetscFunctionBegin;
36
37 PetscOptionsBegin(comm, NULL, "Neo-Hookean physical parameters for smoother", NULL);
38
39 PetscCall(PetscOptionsScalar("-nu_smoother", "Poisson's ratio for smoother", NULL, nu_smoother, &nu_smoother, &nu_flag));
40
41 PetscOptionsEnd(); // End of setting Physics
42
43 if (nu_flag) {
44 // Copy context
45 CeedQFunctionContextGetData(ctx, CEED_MEM_HOST, &phys);
46 PetscCall(PetscMalloc1(1, &phys_smoother));
47 PetscCall(PetscMemcpy(phys_smoother, phys, sizeof(*phys)));
48 CeedQFunctionContextRestoreData(ctx, &phys);
49 // Create smoother context
50 CeedQFunctionContextCreate(ceed, ctx_smoother);
51 phys_smoother->nu = nu_smoother;
52 CeedQFunctionContextSetData(*ctx_smoother, CEED_MEM_HOST, CEED_COPY_VALUES, sizeof(*phys_smoother), phys_smoother);
53 PetscCall(PetscFree(phys_smoother));
54 } else {
55 *ctx_smoother = NULL;
56 }
57
58 PetscFunctionReturn(PETSC_SUCCESS);
59 }
60
61 // Process physics options - Neo-Hookean
ProcessPhysics_NH(MPI_Comm comm,Physics_NH phys,Units units)62 PetscErrorCode ProcessPhysics_NH(MPI_Comm comm, Physics_NH phys, Units units) {
63 PetscBool nu_flag = PETSC_FALSE;
64 PetscBool Young_flag = PETSC_FALSE;
65 phys->nu = 0;
66 phys->E = 0;
67 units->meter = 1; // 1 meter in scaled length units
68 units->second = 1; // 1 second in scaled time units
69 units->kilogram = 1; // 1 kilogram in scaled mass units
70
71 PetscFunctionBeginUser;
72
73 PetscOptionsBegin(comm, NULL, "Neo-Hookean physical parameters", NULL);
74
75 PetscCall(PetscOptionsScalar("-nu", "Poisson's ratio", NULL, phys->nu, &phys->nu, &nu_flag));
76
77 PetscCall(PetscOptionsScalar("-E", "Young's Modulus", NULL, phys->E, &phys->E, &Young_flag));
78
79 PetscCall(PetscOptionsScalar("-units_meter", "1 meter in scaled length units", NULL, units->meter, &units->meter, NULL));
80 units->meter = fabs(units->meter);
81
82 PetscCall(PetscOptionsScalar("-units_second", "1 second in scaled time units", NULL, units->second, &units->second, NULL));
83 units->second = fabs(units->second);
84
85 PetscCall(PetscOptionsScalar("-units_kilogram", "1 kilogram in scaled mass units", NULL, units->kilogram, &units->kilogram, NULL));
86 units->kilogram = fabs(units->kilogram);
87
88 PetscOptionsEnd(); // End of setting Physics
89
90 // Check for all required options to be set
91 if (!nu_flag) {
92 SETERRQ(PETSC_COMM_SELF, PETSC_ERR_SUP, "-nu option needed");
93 }
94 if (!Young_flag) {
95 SETERRQ(PETSC_COMM_SELF, PETSC_ERR_SUP, "-E option needed");
96 }
97
98 // Define derived units
99 units->Pascal = units->kilogram / (units->meter * PetscSqr(units->second));
100
101 // Scale E to Pa
102 phys->E *= units->Pascal;
103
104 PetscFunctionReturn(PETSC_SUCCESS);
105 };