Ferroelectric material possesses spontaneous polarization at room temperature, which can be switched by an external electric field. The diverse domain structures within ferroelectric materials, consisting of polarizations in various directions, often significantly affect their physical properties and practical applications. Numerical simulations can aid in comprehending and validating the complex domains observed in experiments. They can also provide guidance for controlling such structures. One popular method for finding dipole configurations is to create an energy model and employ it in Monte-Carlo simulations to find dipole configuration. However, since these simulations usually reaches the ground state of the system (the state with the lowest global energy), they often miss the dipole configurations of interest, such as topological domain structures, which are usually metastable.

Here, in order to simulate complex domain, we introduce Brown's equation, which is originally used for micromagnetic simulation, into the large-scale simulation of ferroelectric materials. Using the effective Hamiltonian as the energy model, we derive the Brown's equations with respect to the electric dipoles in ferroelectric materials, and invesitgate perovskites such as\rmBaTiO_3 bulk, \rmPbTiO_3 bulk, and \rmSrTiO_3/\rmPbTiO_3/\rmSrTiO_3 sandwiched structures. We demonstrate the reliability and feasibility of Brown's equation in ferroelectrics through the simulation of \rmBaTiO_3 bulk and \rmPbTiO_3 bulk, which are consistent with experiments. Then, using Brown's equation derived in our work, we obtain various domain structures in \rmSrTiO_3/\rmPbTiO_3/\rmSrTiO_3 sandwiched structures, including periodic stripe domains and vortex domains. The simulation results are compared with related exprimental results.