def METHOD_NAME():
    m = ConcreteModel()
    m.fs = FlowsheetBlock(dynamic=False)
    # Set up thermo props and reaction props
    m.fs.solid_properties = SolidPhaseParameterBlock()
    m.fs.solid_state_block = m.fs.solid_properties.build_state_block(
        [0], parameters=m.fs.solid_properties, defined_state=True
    )
    m.fs.gas_properties = GasPhaseParameterBlock()
    m.fs.gas_state_block = m.fs.gas_properties.build_state_block(
        [0], parameters=m.fs.gas_properties, defined_state=True
    )
    m.fs.reactions = HeteroReactionParameterBlock(
        solid_property_package=m.fs.solid_properties,
        gas_property_package=m.fs.gas_properties,
    )
    m.fs.unit = m.fs.reactions.reaction_block_class(
        [0],
        parameters=m.fs.reactions,
        solid_state_block=m.fs.solid_state_block,
        gas_state_block=m.fs.gas_state_block,
        has_equilibrium=False,
    )
    # Fix required variables to make reaction model square
    # (gas mixture and component densities,
    # solid particle porosity, density and component fractions)
    m.fs.gas_state_block[0].dens_mol.fix(10)
    m.fs.gas_state_block[0].dens_mol_comp.fix(10)
    m.fs.solid_state_block[0].particle_porosity.fix(0.27)
    m.fs.solid_state_block[0].mass_frac_comp["Fe2O3"].fix(0.45)
    m.fs.solid_state_block[0].mass_frac_comp["Fe3O4"].fix(1e-9)
    m.fs.solid_state_block[0].mass_frac_comp["Al2O3"].fix(0.55)
    m.fs.solid_state_block[0].dens_mass_skeletal.fix(1)
    m.fs.solid_state_block[0].temperature.fix(1183.15)  # K
    return m