Home
last modified time | relevance | path

Searched refs:nu (Results 1 – 23 of 23) sorted by relevance

/libCEED/examples/solids/qfunctions/
H A Dmanufactured-force.h20 CeedScalar nu; // Poisson's ratio member
38 const CeedScalar nu = context->nu; in SetupMMSForce() local
48 …s(y * 3.0) * exp(z * 4.0) * 4.0 - cos(z * 4.0) * sin(y * 3.0) * exp(x * 2.0) * 8.0) * (nu - 0.5)) / in SetupMMSForce()
49 ((nu * 2.0 - 1.0) * (nu + 1.0)) + in SetupMMSForce()
50 …y * 3.0) * exp(x * 2.0) * (4.5) + sin(x * 2.0) * sin(z * 4.0) * exp(y * 3.0) * 3.0) * (nu - 0.5)) / in SetupMMSForce()
51 ((nu * 2.0 - 1.0) * (nu + 1.0)) + 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 …s(z * 4.0) * exp(x * 2.0) * 3.0 - cos(x * 2.0) * sin(z * 4.0) * exp(y * 3.0) * 2.0) * (nu - 0.5)) / in SetupMMSForce()
[all …]
H A Dlinear.h20 CeedScalar nu; // Poisson's ratio member
40 const CeedScalar nu = context->nu; in ElasResidual_Linear() local
98 const CeedScalar ss = E / ((1 + nu) * (1 - 2 * nu)); in ElasResidual_Linear()
99 const CeedScalar sigma00 = ss * ((1 - nu) * e[0][0] + nu * e[1][1] + nu * e[2][2]), in ElasResidual_Linear()
100 sigma11 = ss * (nu * e[0][0] + (1 - nu) * e[1][1] + nu * e[2][2]), in ElasResidual_Linear()
101 …sigma22 = ss * (nu * e[0][0] + nu * e[1][1] + (1 - nu) * e[2][2]), sigma12 = ss * (1 - 2 * nu) * e… in ElasResidual_Linear()
102 … sigma02 = ss * (1 - 2 * nu) * e[0][2] * 0.5, sigma01 = ss * (1 - 2 * nu) * e[0][1] * 0.5; in ElasResidual_Linear()
137 const CeedScalar nu = context->nu; in ElasJacobian_Linear() local
193 const CeedScalar ss = E / ((1 + nu) * (1 - 2 * nu)); in ElasJacobian_Linear()
194 const CeedScalar dsigma00 = ss * ((1 - nu) * de[0][0] + nu * de[1][1] + nu * de[2][2]), in ElasJacobian_Linear()
[all …]
H A Dfinite-strain-neo-hookean.h20 CeedScalar nu; // Poisson's ratio member
159 const CeedScalar nu = context->nu; 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()
266 const CeedScalar nu = context->nu; 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()
408 const CeedScalar nu = context->nu; 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.h20 CeedScalar nu; // Poisson's ratio member
/libCEED/tests/
H A Dt502-operator-f.f9022 integer nx,nu
24 parameter(nu=nelem*(p-1)+1)
30 real*8 hu(nu*2),hv(nu*2)
57 call ceedelemrestrictioncreate(ceed,nelem,p,2,1,2*nu,ceed_mem_host,&
105 call ceedvectorcreate(ceed,2*nu,u,err)
107 do i=1,nu
112 call ceedvectorcreate(ceed,2*nu,v,err)
118 do i=1,nu
H A Dt505-operator-f.f9022 integer nx,nu
24 parameter(nu=nelem*(p-1)+1)
30 real*8 hv(nu)
57 call ceedelemrestrictioncreate(ceed,nelem,p,1,1,nu,ceed_mem_host,&
105 call ceedvectorcreate(ceed,nu,u,err)
107 call ceedvectorcreate(ceed,nu,v,err)
114 do i=1,nu
128 total=-nu
129 do i=1,nu
H A Dt506-operator-f.f9024 integer nx,nu
26 parameter(nu=nelem*(p-1)+1)
32 real*8 hu(nu*2),hv(nu*2)
59 call ceedelemrestrictioncreate(ceed,nelem,p,2,1,2*nu,ceed_mem_host,&
141 call ceedvectorcreate(ceed,2*nu,u,err)
143 do i=1,nu
148 call ceedvectorcreate(ceed,2*nu,v,err)
159 do i=1,nu
184 do i=1,nu
H A Dt500-operator-f.f9022 integer nx,nu
24 parameter(nu=nelem*(p-1)+1)
30 real*8 hv(nu)
56 call ceedelemrestrictioncreate(ceed,nelem,p,1,1,nu,ceed_mem_host,&
104 call ceedvectorcreate(ceed,nu,u,err)
106 call ceedvectorcreate(ceed,nu,v,err)
110 do i=1,nu
H A Dt503-operator-f.f9022 integer nx,nu
24 parameter(nu=nelem*(p-1)+1)
30 real*8 hv(nu)
40 call ceedvectorcreate(ceed,nu,u,err)
43 call ceedvectorcreate(ceed,nu,v,err)
67 call ceedelemrestrictioncreate(ceed,nelem,p,1,1,nu,ceed_mem_host,&
112 do i=1,nu
H A Dt501-operator-f.f9022 integer nx,nu
24 parameter(nu=nelem*(p-1)+1)
30 real*8 hv(nu)
57 call ceedelemrestrictioncreate(ceed,nelem,p,1,1,nu,ceed_mem_host,&
105 call ceedvectorcreate(ceed,nu,u,err)
107 call ceedvectorcreate(ceed,nu,v,err)
112 do i=1,nu
H A Dt504-operator-f.f9022 integer nx,nu
24 parameter(nu=nelem*(p-1)+1)
51 call ceedelemrestrictioncreate(ceed,nelem,p,2,1,2*nu,ceed_mem_host,&
/libCEED/examples/solids/problems/
H A Dmooney-rivlin.c66 PetscReal nu = -1; in ProcessPhysics_MR() local
84 PetscCall(PetscOptionsScalar("-nu", "Poisson ratio", NULL, nu, &nu, NULL)); in ProcessPhysics_MR()
85 …if (nu < 0 || nu >= 0.5) SETERRQ(PETSC_COMM_SELF, PETSC_ERR_SUP, "Mooney-Rivlin model requires Poi… in ProcessPhysics_MR()
86 phys->lambda = 2 * (phys->mu_1 + phys->mu_2) * nu / (1 - 2 * nu); in ProcessPhysics_MR()
H A Dneo-hookean.c51 phys_smoother->nu = nu_smoother; in PhysicsSmootherContext_NH()
65 phys->nu = 0; in ProcessPhysics_NH()
75 PetscCall(PetscOptionsScalar("-nu", "Poisson's ratio", NULL, phys->nu, &phys->nu, &nu_flag)); in ProcessPhysics_NH()
H A Dneo-hookean.h18 CeedScalar nu; // Poisson's ratio member
/libCEED/examples/solids/
H A DMakefile74 tests-output/NH-strain.csv: args = -problem FSInitial-NH1 -E 2.8 -nu 0.4
75 tests-output/MR-strain.csv: args = -problem FSInitial-MR1 -mu_1 1 -mu_2 0.0 -nu .4
76 tests-output/MR-strain1.csv: args = -problem FSInitial-MR1 -mu_1 .5 -mu_2 0.5 -nu .4
H A DREADME.md20 ./elasticity -mesh [.exo file] -degree [degree] -nu [nu] -E [E] [boundary options] -problem [proble…
41 * - `-nu [real]`
42 - [Poisson's ratio](https://en.wikipedia.org/wiki/Poisson%27s_ratio), $\nu < 0.5$
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 …
123 - Poisson's ratio for multigrid smoothers, $\nu < 0.5$
158 ./elasticity -mesh [mesh] -degree [degree] -nu [nu] -E [E] -forcing mms
H A Dindex.md114 \begin{aligned} \lambda &= \frac{E \nu}{(1 + \nu)(1 - 2 \nu)} \\ \mu &= \frac{E}{2(1 + \nu)} \end{a…
137 Note that the incompressible limit $\nu \to \frac 1 2$ causes $\lambda \to \infty$, and thus $\math…
362 nh["parameters"] = "E=2.8, nu=0.4"
366 mr["parameters"] = "mu_1=1, mu_2=0, nu=.4"
370 mr1["parameters"] = "mu_1=0.5, mu_2=0.5, nu=.4"
/libCEED/examples/fluids/qfunctions/
H A Dstg_shur14.h108 … const CeedScalar nu, CeedScalar *hmax, CeedScalar *ke, CeedScalar *keta, CeedScalar *kcut) { in SpectrumConstants() argument
111 *keta = 2 * M_PI * pow(Cube(nu) / eps, -0.25); in SpectrumConstants()
129 … const CeedScalar nu, CeedScalar qn[], const StgShur14Context stg_ctx) { in CalcSpectrum() argument
134 SpectrumConstants(wall_dist, eps, lt, h_node_sep, nu, &hmax, &ke, &keta, &kcut); in CalcSpectrum()
200 … const CeedScalar eps, const CeedScalar lt, const CeedScalar nu, CeedScalar u[3], in StgShur14Calc_PrecompEktot() argument
208 SpectrumConstants(wall_dist, eps, lt, h_node_sep, nu, &hmax, &ke, &keta, &kcut); in StgShur14Calc_PrecompEktot()
247 const CeedScalar nu = mu / rho; in StgShur14Preprocess() local
266 SpectrumConstants(wall_dist, eps, lt, h_node_sep, nu, &hmax, &ke, &keta, &kcut); in StgShur14Preprocess()
294 const CeedScalar nu = gas->mu / rho; in ICsStg() local
307 CalcSpectrum(x_i[1], eps, lt, h_node_sep, nu, qn, stg_ctx); in ICsStg()
[all …]
H A Dblasius.h79 CeedScalar nu = mu / rho_infty; in BlasiusSolution() local
81 CeedScalar eta = x[1] * sqrt(U_infty / (nu * (x0 + x[0] - x_inflow))); in BlasiusSolution()
89 *t12 = mu * U_infty * f[2] * sqrt(U_infty / (nu * (x0 + x[0] - x_inflow))); in BlasiusSolution()
93 Y[2] = 0.5 * sqrt(nu * U_infty / (x0 + x[0] - x_inflow)) * (eta * f[1] - f[0]); in BlasiusSolution()
/libCEED/python/tests/
H A Dtest-5-operator.py48 nu = nelem * (p - 1) + 1
58 u = ceed.Vector(nu)
59 v = ceed.Vector(nu)
73 ru = ceed.ElemRestriction(nelem, p, 1, 1, nu, indu,
137 nu = nelem * (p - 1) + 1
147 u = ceed.Vector(nu)
148 v = ceed.Vector(nu)
162 ru = ceed.ElemRestriction(nelem, p, 1, 1, nu, indu,
210 for i in range(nu):
227 nu = nelem * (p - 1) + 1
[all …]
/libCEED/examples/fluids/
H A Dindex.md493 where $y^+$ is the wall-friction scaled wall-distance ($y^+ = y u_\tau / \nu = y/\delta_\nu$), $A^+…
494 For this implementation, we assume that $\delta_\nu$ is constant across the wall and is defined by …
682 …^2}{c_p T_w}$ is the modified Eckert number, and $Re_h = \frac{u_{\max}H}{\nu}$ is the Reynolds nu…
765 (\nu^3/\varepsilon)^{-1/4}$ with $\nu$ the kinematic viscosity and
H A DREADME.md668 - Friction length associated with the flow, $\delta_\nu$. Used in wall-damping functions
/libCEED/
H A DREADME.md315 $ ./elasticity -ceed /cpu/self -mesh [.exo file] -degree 2 -E 1 -nu 0.3 -problem Linear -forcing mms
316 $ ./elasticity -ceed /gpu/cuda -mesh [.exo file] -degree 2 -E 1 -nu 0.3 -problem Linear -forcing mms