GEOS
ProppantTransportKernels.hpp
Go to the documentation of this file.
1 /*
2  * ------------------------------------------------------------------------------------------------------------
3  * SPDX-License-Identifier: LGPL-2.1-only
4  *
5  * Copyright (c) 2016-2024 Lawrence Livermore National Security LLC
6  * Copyright (c) 2018-2024 TotalEnergies
7  * Copyright (c) 2018-2024 The Board of Trustees of the Leland Stanford Junior University
8  * Copyright (c) 2023-2024 Chevron
9  * Copyright (c) 2019- GEOS/GEOSX Contributors
10  * All rights reserved
11  *
12  * See top level LICENSE, COPYRIGHT, CONTRIBUTORS, NOTICE, and ACKNOWLEDGEMENTS files for details.
13  * ------------------------------------------------------------------------------------------------------------
14  */
15 
20 #ifndef GEOS_PHYSICSSOLVERS_FLUIDFLOW_PROPPANTTRANSPORT_PROPPANTTRANSPORTKERNELS_HPP_
21 #define GEOS_PHYSICSSOLVERS_FLUIDFLOW_PROPPANTTRANSPORT_PROPPANTTRANSPORTKERNELS_HPP_
22 
23 #include "common/DataTypes.hpp"
24 #include "common/GEOS_RAJA_Interface.hpp"
25 #include "constitutive/fluid/singlefluid/SingleFluidFields.hpp"
26 #include "constitutive/fluid/singlefluid/ParticleFluidBase.hpp"
27 #include "constitutive/fluid/singlefluid/ParticleFluidFields.hpp"
28 #include "constitutive/fluid/singlefluid/SlurryFluidBase.hpp"
29 #include "constitutive/fluid/singlefluid/SlurryFluidFields.hpp"
30 #include "constitutive/permeability/PermeabilityBase.hpp"
31 #include "constitutive/permeability/PermeabilityFields.hpp"
32 #include "finiteVolume/FluxApproximationBase.hpp"
33 #include "physicsSolvers/fluidFlow/FlowSolverBaseFields.hpp"
34 #include "physicsSolvers/fluidFlow/proppantTransport/ProppantTransportFields.hpp"
35 #include "physicsSolvers/fluidFlow/StencilAccessors.hpp"
36 #include "physicsSolvers/PhysicsSolverBaseKernels.hpp"
37 
38 namespace geos
39 {
40 
41 namespace proppantTransportKernels
42 {
43 /******************************** FluidUpdateKernel ********************************/
44 
45 struct FluidUpdateKernel
46 {
47  template< typename FLUID_WRAPPER >
48  static void launch( FLUID_WRAPPER const & fluidWrapper,
49  arrayView1d< real64 const > const & pres,
50  arrayView2d< real64 const > const & componentConcentration )
51  {
52  forAll< parallelDevicePolicy<> >( fluidWrapper.numElems(), [=] GEOS_HOST_DEVICE ( localIndex const a )
53  {
54  localIndex const NC = fluidWrapper.numFluidComponents();
55  stackArray1d< real64, constitutive::SlurryFluidBase::MAX_NUM_COMPONENTS > compConc( NC );
56 
57  for( localIndex c = 0; c < NC; ++c )
58  {
59  compConc[c] = componentConcentration[a][c];
60  }
61 
62  for( localIndex q = 0; q < fluidWrapper.numGauss(); ++q )
63  {
64  fluidWrapper.updateFluidProperty( a, q,
65  pres[a],
66  compConc,
67  0.0 );
68  }
69  } );
70  }
71 };
72 
73 /******************************** ComponentDensityUpdateKernel ********************************/
74 
75 struct ComponentDensityUpdateKernel
76 {
77  template< typename FLUID_WRAPPER >
78  static void launch( FLUID_WRAPPER const & fluidWrapper,
79  arrayView1d< real64 const > const & pres,
80  arrayView2d< real64 const > const & componentConcentration )
81  {
82  forAll< parallelDevicePolicy<> >( fluidWrapper.numElems(), [=] GEOS_HOST_DEVICE ( localIndex const a )
83  {
84  localIndex const NC = fluidWrapper.numFluidComponents();
85  stackArray1d< real64, constitutive::SlurryFluidBase::MAX_NUM_COMPONENTS > compConc( NC );
86 
87  for( localIndex c = 0; c < NC; ++c )
88  {
89  compConc[c] = componentConcentration[a][c];
90  }
91 
92  for( localIndex q = 0; q < fluidWrapper.numGauss(); ++q )
93  {
94  fluidWrapper.updateComponentDensity( a, q,
95  pres[a],
96  compConc );
97  }
98  } );
99  }
100 };
101 
102 /******************************** ProppantUpdateKernel ********************************/
103 
104 struct ProppantUpdateKernel
105 {
106  template< typename PROPPANT_WRAPPER >
107  static void launch( PROPPANT_WRAPPER const & proppantWrapper,
108  arrayView1d< real64 const > const & proppantConc,
109  arrayView2d< real64 const, constitutive::singlefluid::USD_FLUID > const & fluidDens,
110  arrayView2d< real64 const > const & dFluidDens_dPres,
111  arrayView3d< real64 const > const & dFluidDens_dCompConc,
112  arrayView2d< real64 const > const & fluidVisc,
113  arrayView2d< real64 const > const & dFluidVisc_dPres,
114  arrayView3d< real64 const > const & dFluidVisc_dCompConc )
115  {
116 
117  forAll< parallelDevicePolicy<> >( proppantWrapper.numElems(), [=] GEOS_HOST_DEVICE ( localIndex const a )
118  {
119  proppantWrapper.update( a,
120  proppantConc[a],
121  fluidDens[a][0],
122  dFluidDens_dPres[a][0],
123  dFluidDens_dCompConc[a][0],
124  fluidVisc[a][0],
125  dFluidVisc_dPres[a][0],
126  dFluidVisc_dCompConc[a][0] );
127  } );
128  }
129 };
130 
131 /******************************** AccumulationKernel ********************************/
132 
133 struct AccumulationKernel
134 {
135  GEOS_HOST_DEVICE
136  inline
137  static
138  void
139  compute( localIndex const NC,
140  real64 const proppantConc_n,
141  real64 const proppantConcNew,
142  arraySlice1d< real64 const > const & componentDens_n,
143  arraySlice1d< real64 const > const & componentDensNew,
144  arraySlice1d< real64 const > const & GEOS_UNUSED_PARAM( dCompDens_dPres ),
145  arraySlice2d< real64 const > const & dCompDens_dCompConc,
146  real64 const volume,
147  real64 const packPoreVolume,
148  real64 const proppantLiftVolume,
149  arraySlice1d< real64 > const & localAccum,
150  arraySlice2d< real64 > const & localAccumJacobian );
151 
152  static void
153  launch( localIndex const size,
154  localIndex const NC,
155  localIndex const NDOF,
156  globalIndex const rankOffset,
157  arrayView1d< globalIndex const > const & dofNumber,
158  arrayView1d< integer const > const & elemGhostRank,
159  arrayView1d< real64 const > const & proppantConc_n,
160  arrayView1d< real64 const > const & proppantConc,
161  arrayView2d< real64 const > const & componentDens_n,
162  arrayView3d< real64 const > const & componentDens,
163  arrayView3d< real64 const > const & dCompDens_dPres,
164  arrayView4d< real64 const > const & dCompDens_dCompConc,
165  arrayView1d< real64 const > const & volume,
166  arrayView1d< real64 const > const & proppantPackVolFrac,
167  arrayView1d< real64 const > const & proppantLiftFlux,
168  real64 const dt,
169  real64 const maxProppantConcentration,
170  CRSMatrixView< real64, globalIndex const > const & localMatrix,
171  arrayView1d< real64 > const & localRhs );
172 };
173 
174 /******************************** FluxKernel ********************************/
175 
176 struct FluxKernel
177 {
178  using FlowAccessors =
179  StencilAccessors< fields::ghostRank,
180  fields::elementAperture,
181  fields::flow::pressure,
182  fields::flow::gravityCoefficient,
183  fields::proppant::proppantConcentration,
184  fields::proppant::isProppantMobile >;
185 
186  using CellBasedFluxFlowAccessors =
187  StencilAccessors< fields::flow::pressure,
188  fields::flow::gravityCoefficient,
189  fields::elementAperture >;
190 
191  using ParticleFluidAccessors =
192  StencilMaterialAccessors< constitutive::ParticleFluidBase,
193  fields::particlefluid::settlingFactor,
194  fields::particlefluid::dSettlingFactor_dPressure,
195  fields::particlefluid::dSettlingFactor_dProppantConcentration,
196  fields::particlefluid::dSettlingFactor_dComponentConcentration,
197  fields::particlefluid::collisionFactor,
198  fields::particlefluid::dCollisionFactor_dProppantConcentration >;
199 
200  using SlurryFluidAccessors =
201  StencilMaterialAccessors< constitutive::SlurryFluidBase,
202  fields::singlefluid::density,
203  fields::singlefluid::dDensity,
204  fields::slurryfluid::dDensity_dProppantConcentration,
205  fields::slurryfluid::dDensity_dComponentConcentration,
206  fields::singlefluid::viscosity,
207  fields::singlefluid::dViscosity,
208  fields::slurryfluid::dViscosity_dProppantConcentration,
209  fields::slurryfluid::dViscosity_dComponentConcentration,
210  fields::slurryfluid::componentDensity,
211  fields::slurryfluid::dComponentDensity_dPressure,
212  fields::slurryfluid::dComponentDensity_dComponentConcentration,
213  fields::slurryfluid::fluidDensity,
214  fields::slurryfluid::dFluidDensity_dPressure,
215  fields::slurryfluid::dFluidDensity_dComponentConcentration >;
216 
217  using CellBasedFluxSlurryFluidAccessors =
218  StencilMaterialAccessors< constitutive::SlurryFluidBase,
219  fields::singlefluid::density,
220  fields::singlefluid::viscosity >;
221 
222  using PermeabilityAccessors =
223  StencilMaterialAccessors< constitutive::PermeabilityBase,
224  fields::permeability::permeability,
225  fields::permeability::permeabilityMultiplier >;
226 
227  using DerivOffset = constitutive::singlefluid::DerivativeOffsetC< 0 >;
235  template< typename VIEWTYPE >
236  using ElementViewConst = typename ElementRegionManager::ElementViewConst< VIEWTYPE >;
237 
238  template< typename VIEWTYPE >
239  using ElementView = typename ElementRegionManager::ElementView< VIEWTYPE >;
240 
241  static void
242  launch( SurfaceElementStencilWrapper const & stencilWrapper,
243  localIndex const numDofPerCell,
244  real64 const dt,
245  globalIndex const rankOffset,
246  integer const updateProppantPacking,
247  R1Tensor const & unitGravityVector,
248  ElementViewConst< arrayView1d< globalIndex const > > const & dofNumber,
249  ElementViewConst< arrayView1d< integer const > > const & ghostRank,
250  ElementViewConst< arrayView1d< real64 const > > const & pres,
251  ElementViewConst< arrayView1d< real64 const > > const & proppantConc,
252  ElementViewConst< arrayView3d< real64 const > > const & componentDens,
253  ElementViewConst< arrayView3d< real64 const > > const & dComponentDens_dPres,
254  ElementViewConst< arrayView4d< real64 const > > const & dComponentDens_dComponentConc,
255  ElementViewConst< arrayView1d< real64 const > > const & gravDepth,
256  ElementViewConst< arrayView2d< real64 const, constitutive::singlefluid::USD_FLUID > > const & dens,
257  ElementViewConst< arrayView3d< real64 const, constitutive::singlefluid::USD_FLUID_DER > > const & dDens,
258  ElementViewConst< arrayView2d< real64 const > > const & dDens_dProppantConc,
259  ElementViewConst< arrayView3d< real64 const > > const & dDens_dComponentConc,
260  ElementViewConst< arrayView2d< real64 const, constitutive::singlefluid::USD_FLUID > > const & visc,
261  ElementViewConst< arrayView3d< real64 const, constitutive::singlefluid::USD_FLUID_DER > > const & dVisc,
262  ElementViewConst< arrayView2d< real64 const > > const & dVisc_dProppantConc,
263  ElementViewConst< arrayView3d< real64 const > > const & dVisc_dComponentConc,
264  ElementViewConst< arrayView2d< real64 const, constitutive::singlefluid::USD_FLUID > > const & fluidDensity,
265  ElementViewConst< arrayView2d< real64 const > > const & dFluidDens_dPres,
266  ElementViewConst< arrayView3d< real64 const > > const & dFluidDens_dComponentConc,
267  ElementViewConst< arrayView1d< real64 const > > const & settlingFactor,
268  ElementViewConst< arrayView1d< real64 const > > const & dSettlingFactor_dPres,
269  ElementViewConst< arrayView1d< real64 const > > const & dSettlingFactor_dProppantConc,
270  ElementViewConst< arrayView2d< real64 const > > const & dSettlingFactor_dComponentConc,
271  ElementViewConst< arrayView1d< real64 const > > const & collisionFactor,
272  ElementViewConst< arrayView1d< real64 const > > const & dCollisionFactor_dProppantConc,
273  ElementViewConst< arrayView1d< integer const > > const & isProppantMobile,
274  ElementViewConst< arrayView3d< real64 const > > const & permeability,
275  ElementViewConst< arrayView3d< real64 const > > const & permeabilityMultiplier,
276  ElementViewConst< arrayView1d< real64 const > > const & aperture,
277  CRSMatrixView< real64, globalIndex const > const & localMatrix,
278  arrayView1d< real64 > const & localRhs );
279 
280  static void
281  launchCellBasedFluxCalculation( SurfaceElementStencilWrapper const & stencilWrapper,
282  R1Tensor const & unitGravityVector,
283  ElementViewConst< arrayView1d< real64 const > > const & pres,
284  ElementViewConst< arrayView1d< real64 const > > const & gravDepth,
285  ElementViewConst< arrayView2d< real64 const, constitutive::singlefluid::USD_FLUID > > const & dens,
286  ElementViewConst< arrayView2d< real64 const, constitutive::singlefluid::USD_FLUID > > const & visc,
287  ElementViewConst< arrayView3d< real64 const > > const & permeability,
288  ElementViewConst< arrayView3d< real64 const > > const & permeabilityMultiplier,
289  ElementViewConst< arrayView1d< real64 const > > const & aperture,
290  ElementView< arrayView2d< real64 > > const & cellBasedFlux );
291 
298  template< localIndex MAX_NUM_FLUX_ELEMS >
299  GEOS_HOST_DEVICE
300  static void
301  computeJunction( localIndex const numElems,
302  localIndex const numDofPerCell,
303  arraySlice1d< localIndex const > const & stencilElementIndices,
304  arrayView1d< real64 const > const & pres,
305  arrayView1d< real64 const > const & proppantConc,
306  arrayView3d< real64 const > const & componentDens,
307  arrayView3d< real64 const > const & dComponentDens_dPres,
308  arrayView4d< real64 const > const & dComponentDens_dComponentConc,
309  arrayView1d< real64 const > const & gravDepth,
310  arrayView2d< real64 const, constitutive::singlefluid::USD_FLUID > const & dens,
311  arrayView3d< real64 const, constitutive::singlefluid::USD_FLUID_DER > const & dDens,
312  arrayView2d< real64 const > const & dDens_dProppantConc,
313  arrayView3d< real64 const > const & dDens_dComponentConc,
314  arrayView2d< real64 const, constitutive::singlefluid::USD_FLUID > const & visc,
315  arrayView3d< real64 const, constitutive::singlefluid::USD_FLUID_DER > const & dVisc,
316  arrayView2d< real64 const > const & dVisc_dProppantConc,
317  arrayView3d< real64 const > const & dVisc_dComponentConc,
318  arrayView2d< real64 const, constitutive::singlefluid::USD_FLUID > const & fluidDensity,
319  arrayView2d< real64 const > const & dFluidDens_dPres,
320  arrayView3d< real64 const > const & dFluidDens_dComponentConc,
321  arrayView1d< real64 const > const & settlingFactor,
322  arrayView1d< real64 const > const & dSettlingFactor_dPres,
323  arrayView1d< real64 const > const & dSettlingFactor_dProppantConc,
324  arrayView2d< real64 const > const & dSettlingFactor_dComponentConc,
325  arrayView1d< real64 const > const & collisionFactor,
326  arrayView1d< real64 const > const & dCollisionFactor_dProppantConc,
327  arrayView1d< integer const > const & isProppantMobile,
328  real64 const (&transmissibility)[MAX_NUM_FLUX_ELEMS],
329  real64 const (&apertureWeight)[MAX_NUM_FLUX_ELEMS],
330  real64 const (&geometricWeight)[MAX_NUM_FLUX_ELEMS],
331  real64 const dt,
332  arraySlice1d< real64 > const & localFlux,
333  arraySlice2d< real64 > const & localFluxJacobian );
334 
335 
336  template< localIndex MAX_NUM_FLUX_ELEMS >
337  GEOS_HOST_DEVICE
338  static void
339  computeCellBasedFlux( localIndex const numElems,
340  arraySlice1d< localIndex const > const & stencilElementIndices,
341  arrayView1d< real64 const > const & pres_n,
342  arrayView1d< real64 const > const & gravDepth,
343  arrayView2d< real64 const, constitutive::singlefluid::USD_FLUID > const & dens,
344  arrayView2d< real64 const, constitutive::singlefluid::USD_FLUID > const & visc,
345  arraySlice1d< R1Tensor const > const & cellCenterToEdgeCenters,
346  real64 const (&transmissibility)[MAX_NUM_FLUX_ELEMS],
347  real64 const (&apertureWeight)[MAX_NUM_FLUX_ELEMS],
348  real64 const (&geometricWeight)[MAX_NUM_FLUX_ELEMS],
349  arrayView2d< real64 > const & cellBasedFlux );
350 };
351 
352 struct ProppantPackVolumeKernel
353 {
354 
355  using FlowAccessors =
356  StencilAccessors< fields::ghostRank,
357  fields::elementAperture,
358  fields::elementVolume,
359  fields::proppant::cellBasedFlux,
360  fields::proppant::isProppantMobile,
361  fields::proppant::isProppantBoundary >;
362 
363  using ParticleFluidAccessors =
364  StencilMaterialAccessors< constitutive::ParticleFluidBase, fields::particlefluid::settlingFactor >;
365 
366  using SlurryFluidAccessors =
367  StencilMaterialAccessors< constitutive::SlurryFluidBase,
368  fields::singlefluid::density,
369  fields::slurryfluid::fluidDensity,
370  fields::slurryfluid::fluidViscosity >;
371 
372 
373  template< typename VIEWTYPE >
374  using ElementView = ElementRegionManager::ElementView< VIEWTYPE >;
375 
376  template< typename VIEWTYPE >
377  using ElementViewConst = ElementRegionManager::ElementViewConst< VIEWTYPE >;
378 
379  static void
380  launchProppantPackVolumeCalculation( SurfaceElementStencil const & stencil,
381  real64 const dt,
382  real64 const proppantDensity,
383  real64 const proppantDiameter,
384  real64 const maxProppantConcentration,
385  R1Tensor const & unitGravityVector,
386  real64 const criticalShieldsNumber,
387  real64 const fricitonCoefficient,
388  ElementView< arrayView1d< real64 const > > const & settlingFactor,
389  ElementView< arrayView2d< real64 const, constitutive::singlefluid::USD_FLUID > > const & density,
390  ElementView< arrayView2d< real64 const, constitutive::singlefluid::USD_FLUID > > const & fluidDensity,
391  ElementView< arrayView1d< integer const > > const & isProppantMobile,
392  ElementView< arrayView1d< integer const > > const & isProppantBoundaryElement,
393  ElementView< arrayView1d< real64 const > > const & aperture,
394  ElementView< arrayView1d< real64 const > > const & volume,
395  ElementView< arrayView1d< integer const > > const & elemGhostRank,
396  ElementView< arrayView2d< real64 const > > const & cellBasedFlux,
397  ElementView< arrayView1d< real64 > > const & conc,
398  ElementView< arrayView1d< real64 > > const & proppantPackVolFrac,
399  ElementView< arrayView1d< real64 > > const & proppantExcessPackVolume,
400  ElementView< arrayView1d< real64 > > const & proppantLiftFlux );
401 
402  static void
403  launchProppantPackVolumeUpdate( SurfaceElementStencil const & stencil,
404  R1Tensor const & unitGravityVector,
405  real64 const maxProppantConcentration,
406  ElementView< arrayView1d< integer const > > const & isProppantMobile,
407  ElementView< arrayView1d< real64 const > > const & proppantExcessPackVolume,
408  ElementView< arrayView1d< real64 > > const & conc,
409  ElementView< arrayView1d< real64 > > const & proppantPackVolFrac );
410 
411  GEOS_HOST_DEVICE
412  inline
413  static
414  void
415  computeProppantPackVolume( localIndex const numElems,
416  real64 const dt,
417  real64 const proppantDensity,
418  real64 const proppantDiameter,
419  real64 const maxProppantConcentration,
420  R1Tensor const & unitGravityVector,
421  real64 const criticalShieldsNumber,
422  real64 const frictionCoefficient,
423  arraySlice1d< localIndex const > const & stencilElementIndices,
424  arraySlice1d< real64 const > const & stencilWeights,
425  arraySlice1d< R1Tensor const > const & stencilCellCenterToEdgeCenters,
426  arrayView1d< real64 const > const & settlingFactor,
427  arrayView2d< real64 const, constitutive::singlefluid::USD_FLUID > const & density,
428  arrayView2d< real64 const, constitutive::singlefluid::USD_FLUID > const & fluidDensity,
429  arrayView1d< real64 const > const & volume,
430  arrayView1d< real64 const > const & aperture,
431  arrayView1d< integer const > const & elemGhostRank,
432  arrayView1d< integer const > const & isProppantBoundaryElement,
433  arrayView1d< integer const > const & isProppantMobile,
434  arrayView2d< real64 const > const & cellBasedFlux,
435  arrayView1d< real64 > const & conc,
436  arrayView1d< real64 > const & proppantPackVolFrac,
437  arrayView1d< real64 > const & proppantExcessPackVolume,
438  arrayView1d< real64 > const & proppantLiftFlux );
439 
440  GEOS_HOST_DEVICE
441  inline
442  static
443  void
444  updateProppantPackVolume( localIndex const numElems,
445  arraySlice1d< localIndex const > const & stencilElementIndices,
446  arraySlice1d< real64 const > const & stencilWeights,
447  arraySlice1d< R1Tensor const > const & stencilCellCenterToEdgeCenters,
448  R1Tensor const & unitGravityVector,
449  real64 const maxProppantConcentration,
450  arrayView1d< integer const > const & isProppantMobile,
451  arrayView1d< real64 const > const & proppantExcessPackVolume,
452  arrayView1d< real64 > const & conc,
453  arrayView1d< real64 > const & proppantPackVolFrac );
454 };
455 
456 /******************************** ResidualNormKernel ********************************/
457 
461 class ResidualNormKernel : public physicsSolverBaseKernels::ResidualNormKernelBase< 1 >
462 {
463 public:
464 
465  using Base = physicsSolverBaseKernels::ResidualNormKernelBase< 1 >;
466  using Base::m_minNormalizer;
467  using Base::m_rankOffset;
468  using Base::m_localResidual;
469  using Base::m_dofNumber;
470 
471  ResidualNormKernel( globalIndex const rankOffset,
472  arrayView1d< real64 const > const & localResidual,
473  arrayView1d< globalIndex const > const & dofNumber,
474  arrayView1d< localIndex const > const & ghostRank,
475  integer const numComp,
476  ElementSubRegionBase const & subRegion,
477  real64 const minNormalizer )
478  : Base( rankOffset,
479  localResidual,
480  dofNumber,
481  ghostRank,
482  minNormalizer ),
483  m_numComp( numComp ),
484  m_volume( subRegion.getElementVolume() )
485  {}
486 
487  GEOS_HOST_DEVICE
488  virtual void computeLinf( localIndex const ei,
489  LinfStackVariables & stack ) const override
490  {
491  real64 const normalizer = LvArray::math::max( m_minNormalizer, m_volume[ei] );
492  for( integer idof = 0; idof < m_numComp; ++idof )
493  {
494  real64 const valMass = LvArray::math::abs( m_localResidual[stack.localRow + idof] ) / normalizer;
495  if( valMass > stack.localValue[0] )
496  {
497  stack.localValue[0] = valMass;
498  }
499  }
500  }
501 
502  GEOS_HOST_DEVICE
503  virtual void computeL2( localIndex const ei,
504  L2StackVariables & stack ) const override
505  {
506  real64 const normalizer = LvArray::math::max( m_minNormalizer, m_volume[ei] );
507  for( integer idof = 0; idof < m_numComp; ++idof )
508  {
509  stack.localValue[0] += m_localResidual[stack.localRow + idof] * m_localResidual[stack.localRow + idof];
510  stack.localNormalizer[0] += normalizer;
511  }
512  }
513 
514 
515 protected:
516 
518  integer const m_numComp;
519 
521  arrayView1d< real64 const > const m_volume;
522 };
523 
527 class ResidualNormKernelFactory
528 {
529 public:
530 
543  template< typename POLICY >
544  static void
545  createAndLaunch( physicsSolverBaseKernels::NormType const normType,
546  integer const numComp,
547  globalIndex const rankOffset,
548  string const & dofKey,
549  arrayView1d< real64 const > const & localResidual,
550  ElementSubRegionBase const & subRegion,
551  real64 const minNormalizer,
552  real64 (& residualNorm)[1],
553  real64 (& residualNormalizer)[1] )
554  {
555  arrayView1d< globalIndex const > const dofNumber = subRegion.getReference< array1d< globalIndex > >( dofKey );
556  arrayView1d< integer const > const ghostRank = subRegion.ghostRank();
557 
558  ResidualNormKernel kernel( rankOffset, localResidual, dofNumber, ghostRank,
559  numComp, subRegion, minNormalizer );
560  if( normType == physicsSolverBaseKernels::NormType::Linf )
561  {
562  ResidualNormKernel::launchLinf< POLICY >( subRegion.size(), kernel, residualNorm );
563  }
564  else // L2 norm
565  {
566  ResidualNormKernel::launchL2< POLICY >( subRegion.size(), kernel, residualNorm, residualNormalizer );
567  }
568  }
569 
570 };
571 
572 
573 
574 } // namespace proppantTransportKernels
575 
576 } // namespace geos
577 
578 #endif //GEOS_PHYSICSSOLVERS_FLUIDFLOW_PROPPANTTRANSPORT_PROPPANTTRANSPORTKERNELS_HPP_
GEOS_HOST_DEVICE
#define GEOS_HOST_DEVICE
Marks a host-device function.
Definition: GeosxMacros.hpp:49
GEOS_UNUSED_PARAM
#define GEOS_UNUSED_PARAM(X)
Mark an unused argument and silence compiler warnings.
Definition: GeosxMacros.hpp:72
geos::physicsSolverBaseKernels::NormType
NormType
Type of norm used to check convergence TODO: find a way to put this inside the class.
Definition: PhysicsSolverBaseKernels.hpp:335
geos::ElementRegionManager::ElementView
typename ElementViewAccessor< VIEWTYPE >::NestedViewType ElementView
The ElementViewAccessor at the ElementRegionManager level is the type resulting from ElementViewAcces...
Definition: ElementRegionManager.hpp:73
geos::ElementRegionManager::ElementViewConst
typename ElementViewAccessor< VIEWTYPE >::NestedViewTypeConst ElementViewConst
The ElementViewAccessor at the ElementRegionManager level is the type resulting from ElementViewAcces...
Definition: ElementRegionManager.hpp:81
geos::ElementSubRegionBase::getElementVolume
arrayView1d< real64 const > getElementVolume() const
Get the volume of each element in this subregion.
Definition: ElementSubRegionBase.hpp:160
geos::ObjectManagerBase::ghostRank
array1d< integer > const & ghostRank()
Get the ghost information of each object.
Definition: ObjectManagerBase.hpp:854
geos::StencilAccessors
A struct to automatically construct and store element view accessors.
Definition: StencilAccessors.hpp:39
geos::StencilMaterialAccessors
A struct to automatically construct and store element view accessors.
Definition: StencilAccessors.hpp:98
geos::SurfaceElementStencil
Provides management of the interior stencil points for a face elements when using Two-Point flux appr...
Definition: SurfaceElementStencil.hpp:230
geos::SurfaceElementStencilWrapper
Provides access to the SurfaceElementStencil that may be called from a kernel function.
Definition: SurfaceElementStencil.hpp:40
geos::dataRepository::Group::getReference
GEOS_DECLTYPE_AUTO_RETURN getReference(LOOKUP_TYPE const &lookup) const
Look up a wrapper and get reference to wrapped object.
Definition: Group.hpp:1275
geos::dataRepository::Group::size
localIndex size() const
Get the "size" of the group, which determines the number of elements in resizable wrappers.
Definition: Group.hpp:1317
geos::physicsSolverBaseKernels::ResidualNormKernelBase
Define the base interface for the residual calculations.
Definition: PhysicsSolverBaseKernels.hpp:42
geos::physicsSolverBaseKernels::ResidualNormKernelBase< 1 >::m_minNormalizer
real64 const m_minNormalizer
Value used to make sure that normalizers are never zero.
Definition: PhysicsSolverBaseKernels.hpp:236
geos::physicsSolverBaseKernels::ResidualNormKernelBase< 1 >::m_dofNumber
arrayView1d< globalIndex const > const m_dofNumber
View on the dof numbers.
Definition: PhysicsSolverBaseKernels.hpp:230
geos::physicsSolverBaseKernels::ResidualNormKernelBase< 1 >::ghostRank
GEOS_HOST_DEVICE integer ghostRank(localIndex const i) const
Getter for the ghost rank.
Definition: PhysicsSolverBaseKernels.hpp:91
geos::physicsSolverBaseKernels::ResidualNormKernelBase< 1 >::m_localResidual
arrayView1d< real64 const > const m_localResidual
View on the local residual.
Definition: PhysicsSolverBaseKernels.hpp:227
geos::proppantTransportKernels::ResidualNormKernelFactory::createAndLaunch
static void createAndLaunch(physicsSolverBaseKernels::NormType const normType, integer const numComp, globalIndex const rankOffset, string const &dofKey, arrayView1d< real64 const > const &localResidual, ElementSubRegionBase const &subRegion, real64 const minNormalizer, real64(&residualNorm)[1], real64(&residualNormalizer)[1])
Create a new kernel and launch.
Definition: ProppantTransportKernels.hpp:545
geos::proppantTransportKernels::ResidualNormKernel::computeL2
virtual GEOS_HOST_DEVICE void computeL2(localIndex const ei, L2StackVariables &stack) const override
Compute the local values and normalizer for the L2 norm.
Definition: ProppantTransportKernels.hpp:503
geos::proppantTransportKernels::ResidualNormKernel::m_numComp
integer const m_numComp
Number of fluid components.
Definition: ProppantTransportKernels.hpp:518
geos::proppantTransportKernels::ResidualNormKernel::computeLinf
virtual GEOS_HOST_DEVICE void computeLinf(localIndex const ei, LinfStackVariables &stack) const override
Compute the local values for the Linf norm.
Definition: ProppantTransportKernels.hpp:488
geos::proppantTransportKernels::ResidualNormKernel::m_volume
arrayView1d< real64 const > const m_volume
View on the element volume.
Definition: ProppantTransportKernels.hpp:521
geos::arrayView1d
ArrayView< T, 1 > arrayView1d
Alias for 1D array view.
Definition: DataTypes.hpp:180
geos::CRSMatrixView
LvArray::CRSMatrixView< T, COL_INDEX, localIndex const, LvArray::ChaiBuffer > CRSMatrixView
Alias for CRS Matrix View.
Definition: DataTypes.hpp:310
geos::globalIndex
GEOS_GLOBALINDEX_TYPE globalIndex
Global index type (for indexing objects across MPI partitions).
Definition: DataTypes.hpp:88
geos::stackArray1d
StackArray< T, 1, MAXSIZE > stackArray1d
Alias for 1D stack array.
Definition: DataTypes.hpp:188
geos::arraySlice2d
ArraySlice< T, 2, USD > arraySlice2d
Alias for 2D array slice.
Definition: DataTypes.hpp:200
geos::real64
double real64
64-bit floating point type.
Definition: DataTypes.hpp:99
geos::localIndex
GEOS_LOCALINDEX_TYPE localIndex
Local index type (for indexing objects within an MPI partition).
Definition: DataTypes.hpp:85
geos::arraySlice1d
ArraySlice< T, 1, USD > arraySlice1d
Alias for 1D array slice.
Definition: DataTypes.hpp:184
geos::integer
std::int32_t integer
Signed integer type.
Definition: DataTypes.hpp:82
geos::arrayView4d
ArrayView< T, 4, USD > arrayView4d
Alias for 4D array view.
Definition: DataTypes.hpp:228
geos::arrayView2d
ArrayView< T, 2, USD > arrayView2d
Alias for 2D array view.
Definition: DataTypes.hpp:196
geos::array1d
Array< T, 1 > array1d
Alias for 1D array.
Definition: DataTypes.hpp:176
geos::arrayView3d
ArrayView< T, 3, USD > arrayView3d
Alias for 3D array view.
Definition: DataTypes.hpp:212
geos::fields::elementAperture
Trait struct for elementAperture data.
Definition: MeshFields.hpp:116
geos::fields::elementVolume
Trait struct for elementVolume data.
Definition: MeshFields.hpp:108
geos::fields::ghostRank
Trait struct for ghostRank data.
Definition: MeshFields.hpp:100
geos::proppantTransportKernels::FluxKernel::computeJunction
static GEOS_HOST_DEVICE void computeJunction(localIndex const numElems, localIndex const numDofPerCell, arraySlice1d< localIndex const > const &stencilElementIndices, arrayView1d< real64 const > const &pres, arrayView1d< real64 const > const &proppantConc, arrayView3d< real64 const > const &componentDens, arrayView3d< real64 const > const &dComponentDens_dPres, arrayView4d< real64 const > const &dComponentDens_dComponentConc, arrayView1d< real64 const > const &gravDepth, arrayView2d< real64 const, constitutive::singlefluid::USD_FLUID > const &dens, arrayView3d< real64 const, constitutive::singlefluid::USD_FLUID_DER > const &dDens, arrayView2d< real64 const > const &dDens_dProppantConc, arrayView3d< real64 const > const &dDens_dComponentConc, arrayView2d< real64 const, constitutive::singlefluid::USD_FLUID > const &visc, arrayView3d< real64 const, constitutive::singlefluid::USD_FLUID_DER > const &dVisc, arrayView2d< real64 const > const &dVisc_dProppantConc, arrayView3d< real64 const > const &dVisc_dComponentConc, arrayView2d< real64 const, constitutive::singlefluid::USD_FLUID > const &fluidDensity, arrayView2d< real64 const > const &dFluidDens_dPres, arrayView3d< real64 const > const &dFluidDens_dComponentConc, arrayView1d< real64 const > const &settlingFactor, arrayView1d< real64 const > const &dSettlingFactor_dPres, arrayView1d< real64 const > const &dSettlingFactor_dProppantConc, arrayView2d< real64 const > const &dSettlingFactor_dComponentConc, arrayView1d< real64 const > const &collisionFactor, arrayView1d< real64 const > const &dCollisionFactor_dProppantConc, arrayView1d< integer const > const &isProppantMobile, real64 const (&transmissibility)[MAX_NUM_FLUX_ELEMS], real64 const (&apertureWeight)[MAX_NUM_FLUX_ELEMS], real64 const (&geometricWeight)[MAX_NUM_FLUX_ELEMS], real64 const dt, arraySlice1d< real64 > const &localFlux, arraySlice2d< real64 > const &localFluxJacobian)
Compute flux and its derivatives for a given multi-element connector.
geos::proppantTransportKernels::FluxKernel::ElementViewConst
typename ElementRegionManager::ElementViewConst< VIEWTYPE > ElementViewConst
The type for element-based non-constitutive data parameters. Consists entirely of ArrayView's.
Definition: ProppantTransportKernels.hpp:236
Generated by doxygen 1.9.1