def METHOD_NAME(db):
    base_obj = grab_thermo_Liq_Sol_FpcTP_base(db)
    # Our components for this problem are as follows:
    comp_list = ["H2O", "H_+", "OH_-", "Ca[OH]2", "Ca_2+", "CaOH_+"]
    # In this example, both Sol and Liq phases will be considered valid
    #   when searching for reactions in the database.
    phase_list = ["Liq"]
    base_obj = get_components_and_add_to_idaes_config(db, base_obj, comp_list)
    # Create a reaction config
    react_base = grab_base_reaction_config(db)
    # Add reactions to the reaction base as 'equilibrium'
    react_base = add_equilibrium_reactions_by_phase_to_react_base(
        db, react_base, comp_list, phase_list
    )
    thermo_config = base_obj.idaes_config
    reaction_config = react_base.idaes_config
    model = build_equilibrium_model(thermo_config, reaction_config)
    return model