DensityUpdateTwoConstraints.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 DensityUpdateTwoConstraints : public ElementUserObject
{
public:
  static InputParameters validParams();

  DensityUpdateTwoConstraints(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;
  const VariableName _cost_density_sensitivity_name;
  const VariableName _cost_name;

  MooseWritableVariable * _design_density;
  const MooseWritableVariable * _density_sensitivity;
  const MooseWritableVariable * _cost_density_sensitivity;
  const MooseVariable & _cost;
  const Real _volume_fraction;
  const Real _cost_fraction;
  // Relative tolerance used to stop the two-variable bisection method.
  const Real _relative_tolerance;
  // Weights to solve a dual physics problem, e.g. a thermomechanical one.
  std::vector<Real> _weight_values;

private:
  struct ElementData
  {
    Real old_density;
    Real sensitivity;
    Real cost_sensitivity;
    Real thermal_sensitivity;
    Real cost;
    Real volume;
    Real new_density;

    ElementData() = default;
    ElementData(
        Real dens, Real sens, Real cost_sens, Real thermal_sens, Real cst, Real vol, Real filt_dens)
      : old_density(dens),
        sensitivity(sens),
        cost_sensitivity(cost_sens),
        thermal_sensitivity(thermal_sens),
        cost(cst),
        volume(vol),
        new_density(filt_dens)
    {
    }
  };

  void gatherElementData();
  void performOptimCritLoop();

  Real computeUpdatedDensity(Real current_density,
                             Real dc,
                             Real cost_sensitivity,
                             Real thermal_sensitivity,
                             Real cost,
                             Real lmid,
                             Real cmid);

  Real _total_allowable_volume;

  Real _total_allowable_cost;

  std::map<dof_id_type, ElementData> _elem_data_map;

  const Real _lower_bound;
  const Real _upper_bound;
  const Real _bisection_move;
  const bool _adaptive_move;
  VariableName _thermal_sensitivity_name;
  MooseVariable * _thermal_sensitivity;
};