synthetik-technologies / blastfoam

A CFD solver for multi-component compressible flow with application to high-explosive detonation, explosive safety and air blast
Other
215 stars 96 forks source link

How can I find the equations in blastFSIFoam #91

Closed zyy212121 closed 3 days ago

zyy212121 commented 1 month ago

Thank you for providing the algorithm; it has been very helpful. I am learning the process of coding to add some source terms to the scalar equations in the solver for a specific purpose. I am primarily using the blastFSIFoam solver. However, I found that fluid.solve.H does not directly include the equation assembly and solving. I discovered the timeIntegrators[i].integrate(); function, and in timeIntegrator.C, I found that this is a process for solving sub-time steps, updating fields such as pressure, velocity, displacement, etc., at each sub-time step through system_.solve(). Is this how it works? Have I made any mistakes in my understanding so far?

Furthermore, I have not found the definition of the solve() method, neither in the derived class nor in the base class timeIntegrationSystem. If I need to modify the scalar equation for a scalar fluid, for example, adding a temperature source term related to pressure, where should I make these modifications? Your response would be greatly appreciated, as it will be very helpful to me.