- zero-sized faMeshSubset and fvMeshSubset had READ_IF_PRESENT instead
  of simply copying the schemes/solution setting from the baseMesh

- as part of the style change, include guards are increasingly
  prefixed with 'Foam_'

STYLE: replace stray 'Web:' entries with 'Website:'

sergio authored Jun 14, 2022 and Kutalmış Berçin committed Jun 21, 2022

ENH: saturationModels: new library for saturation models

sergio authored Apr 26, 2022 and Kutalmış Berçin committed Jun 21, 2022

TUT: icoReactingMultiphaseInterFoam: update thermophysicalProperties.liquid

sergio authored Jun 14, 2022 and Kutalmış Berçin committed Jun 21, 2022

phaseSystemModels function objects are relocated within
functionObjects in order to enable broader usage.

ENH: multiphaseInterHtcModel: new heatTransferCoeff function object model
COMP: createExternalCoupledPatchGeometry: add new dependencies
COMP: alphaContactAngle: avoid duplicate entries between multiphaseEuler and reactingEuler
TUT: damBreak4Phase: rename alphaContactAngle as multiphaseEuler::alphaContactAngle

sergio authored Jun 14, 2022 and Kutalmış Berçin committed Jun 21, 2022

Co-authored-by: Kutalmis Bercin <kutalmis.bercin@esi-group.com>

sergio authored Jun 14, 2022 and Kutalmış Berçin committed Jun 21, 2022

thermoTools is a relocation of various existing tools:

- src/TurbulenceModels/compressible/turbulentFluidThermoModels/derivedFvPatchFields/
- src/semiPermeableBaffle/derivedFvPatchFields/
- src/thermophysicalModels/thermophysicalPropertiesFvPatchFields/liquidProperties/

ENH: Allwmake: reordering various compilation steps
Co-authored-by: Kutalmis Bercin <kutalmis.bercin@esi-group.com>

This is on
- incompressible/pimpleFoam/laminar/mixerVesselAMI2D/mixerVesselAMI2D-topologyChange
- redistributePar -reconstruct
where the fvMesh::updateMesh does an early trigger of
mesh.phi() calculation

Specific to the VOF-to-lagrangian FO is to generate particles
which potentially do not relate to the mesh. So here they
are preserved instead of trying to locate them on the
reconstructed mesh. Note: this has the same effect
of actually copying the file...

sergio authored Jan 31, 2022 and Andrew Heather committed Jun 15, 2022

ENH: Adding report after matrix smoothing

BUG: inconsistent "U" name on re-read (buoyancyEnergy)

STYLE: fix spelling inconsistencies

sergio authored Sep 20, 2021 and Kutalmış Berçin committed Jun 14, 2022

speciesSorption is a zeroGradient BC which absorbs mass given by a first
order time derivative, absoprtion rate and an equilibrium value
calculated based on internal species values next to the wall.

patchCellsSource is a source fvOption which applies to the corresponding
species and apply the source calculated on the speciesSorption BC.

A new abstract virtual class was created to group BC's which
don't introduce a source to the matrix (i.e zeroGradient) but calculate
a mass sink/source which should be introduced into the matrix. This
is  done through the fvOption patchCellsSource.

Mark OLESEN authored Jun 08, 2022 and Andrew Heather committed Jun 13, 2022

- this allows the "relocation" of sampled surfaces. For example,
  to reposition into a different coordinate system for importing
  into CAD.

- incorporate output scaling for all surface writer types.

  This was previously done on an adhoc basis for different writers,
  but with now included in the base-level so that all writers
  can automatically use scale + transform.

  Example:

  formatOptions
  {
      vtk
      {
          scale 1000;  // m -> mm
          transform
          {
              origin  (0.05 0 0);
              rotation axisAngle;
              axis    (0 0 1);
              angle   -45;
          }
      }
  }

Vaggelis Papoutsis authored Jun 06, 2022 and Andrew Heather committed Jun 10, 2022

Necessary when grad(Ua) is cached

Vaggelis Papoutsis authored Jun 01, 2022 and Andrew Heather committed Jun 10, 2022

Partially fixes #2502 (see discussion there)

Vaggelis Papoutsis authored May 18, 2022 and Andrew Heather committed Jun 10, 2022

in RASModelVariables were doing this by checking whether the
corresponding pointer was allocated. In some cases, however, even if the
field does not exist, the pointer is not null, leading to the wrong
output. Made the correspding functions virtual and overwritten their
return values in the derived classes. Kept the initial implementation in
base to facilitate the clone function.

Vaggelis Papoutsis authored Apr 11, 2022 and Andrew Heather committed Jun 10, 2022

in cases with more than one primal or adjoint solvers

TUT: removed all occurances of useSolverNameForFields
from the optimisation tutorials since it is now set
automatically.

Vaggelis Papoutsis authored Dec 07, 2021 and Andrew Heather committed Jun 10, 2022

volBSplinesBase const, since the no private/protected
variable is changed there.

Vaggelis Papoutsis authored Dec 03, 2021 and Andrew Heather committed Jun 10, 2022

in the sensitivity patches, symmetry::evaluate() needs access to the
internalField which does exist, leading to wrong memory access.

Fixed by specifying a calculated type fvPatchField for all patches when
creating a boundaryField<Type>

Using a symmetry(Plane) as a sensitivity patch is quite rare and
borderline wrong, but this provides a fix nonetheless.

Vaggelis Papoutsis authored Dec 09, 2021 and Andrew Heather committed Jun 10, 2022

The multiplier of grad(dxdb) is a volTensorField which, by itself, is
memory consuming. The function computing it though was sloppy in terms
of memory management, constituting the peak memory consumption during an
adjoint optimisation. Initial changes to remedy the problem include the
deallocation of some of the volTensorFields included in the computation
of grad(dxdb) once unneeded, the utilisation of volSymmTensorFields
instead of volTensorFields where possible and avoiding allocating some
unnecessary intermediate fields.

Actions to further reduce memory consumption:

- For historical reasons, the code computes/stores the transpose of
  grad(dxdb), which is then transposed when used in the computation of
  the FI or the ESI sensitivity derivatives.  This redundant
  transposition can be avoid, saving the allocation of an additional
  volTensorField, but the changes need to permeate a number of places in
  the code that contribute to grad(dxdb) (e.g.  ATC, adjoint turbulence
  models, adjoint MRF, etc).
- Allocation of unnecessary pointers in the objective class should be
  avoided.