In this work, the Extended Direct Algebraic Method (EDAM) is utilized to analyze and solve the fractional (2+1)-dimensional Conformable Broer-Kaup-Kupershmit System (CBKKS) and investigate different types of traveling wave solutions and study the soliton like-solutions. Using the suggested method, the fractional nonlinear partial differential equation (FNPDE) is primarily reduced to an integer-order nonlinear ordinary differential equation (NODE) under the traveling wave transformation, yielding an algebraic system of nonlinear equations. The ensuing algebraic systems are then solved to construct some families of soliton-like solutions and many other physical solutions. Some derived solutions are numerically analyzed using suitable values for the related parameters. The discovered soliton solutions grasp vital importance in fluid mechanics as they offer significant insight into the nonlinear behavior of the targeted model, opening the way for a deeper comprehension of complex physical phenomena and offering valuable applications in the associated areas.

The present study investigates the fractional Dullin-Gottwald-Holm equation by using the Riccati-Bernoulli sub-optimal differential equation method with the Bäcklund transformation. By employing a well-established criterion, the present study reveals novel cusp soliton solutions that resemble peakons and offers valuable insights into their dynamic behaviors and mysterious phenomena. The solution family encompasses various analytical solutions, such as peakons, periodic, and kink-wave solutions. Furthermore, the impact of both the time- and space-fractional parameters on all derived solutions' profiles is examined. This investigation's significance lies in its contribution to understanding intricate dynamics inside physical systems, offering valuable insights into various domains like fluid mechanics and nonlinear phenomena across different physical models. The computational technique's straightforward, effective, and concise nature is demonstrated through introduction of some graphical representations in two- and three-dimensional plots generated by adjusting the related parameters. The findings underscore the versatility of this methodology and demonstrate its applicability as a tool to solve more complicated nonlinear problems as well as its ability to explain many mysterious phenomena.