DensityUpdate.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html

#pragma once

#include "ElementUserObject.h"
#include "MooseTypes.h"

class DensityUpdate : public ElementUserObject
{
public:
  static InputParameters validParams();

  DensityUpdate(const InputParameters & parameters);

  virtual void initialize() override{};
  virtual void timestepSetup() override;
  virtual void execute() override;
  virtual void finalize() override{};
  virtual void threadJoin(const UserObject &) override{};

protected:
  const MooseMesh & _mesh;
  const VariableName _design_density_name;
  const VariableName _density_sensitivity_name;
  MooseWritableVariable * _design_density;
  const MooseWritableVariable * _density_sensitivity;
  const Real _volume_fraction;

private:
  struct ElementData
  {
    Real old_density;
    Real sensitivity;
    Real volume;
    Real new_density;
    ElementData() = default;
    ElementData(Real dens, Real sens, Real vol, Real filt_dens)
      : old_density(dens), sensitivity(sens), volume(vol), new_density(filt_dens)
    {
    }
  };

  void gatherElementData();

  void performOptimCritLoop();

  Real computeUpdatedDensity(Real current_density, Real dc, Real lmid);

  Real _total_allowable_volume;

  std::map<dof_id_type, ElementData> _elem_data_map;

  const Real _lower_bound;
  const Real _upper_bound;
};