PICurv 0.1.0
A Parallel Particle-In-Cell Solver for Curvilinear LES
Loading...
Searching...
No Matches
Functions
implicitsolvers.h File Reference
#include "variables.h"
#include "logging.h"
#include "rhs.h"
#include "Boundaries.h"
Include dependency graph for implicitsolvers.h:
This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Functions

PetscErrorCode RungeKutta (UserCtx *user, IBMNodes *ibm, FSInfo *fsi)
 Advances the momentum equations using an explicit 4th-order Runge-Kutta scheme.
 
PetscErrorCode ImplicitMomentumSolver (UserCtx *user, IBMNodes *ibm, FSInfo *fsi)
 Advances the momentum equations using an implicit scheme.
 

Function Documentation

◆ RungeKutta()

PetscErrorCode RungeKutta ( UserCtx *  user,
IBMNodes *  ibm,
FSInfo *  fsi 
)
extern

Advances the momentum equations using an explicit 4th-order Runge-Kutta scheme.

Parameters
userArray of UserCtx structs for all blocks.
ibm(Optional) Pointer to IBM data. Pass NULL if disabled.
fsi(Optional) Pointer to FSI data. Pass NULL if disabled.
Returns
PetscErrorCode 0 on success.

Advances the momentum equations using an explicit 4th-order Runge-Kutta scheme.

This function computes an intermediate, non-divergence-free contravariant velocity field (Ucont) at time t_{n+1} for all computational blocks.

This is a minimally-edited version of the legacy solver. It retains its internal loop over all blocks and is intended to be called once per time step from the main FlowSolver orchestrator. All former global variables are now accessed via the SimCtx passed in through the first block's UserCtx.

Parameters
userThe array of UserCtx structs for all blocks.
ibm(Optional) Pointer to the full array of IBM data structures. Pass NULL if disabled.
fsi(Optional) Pointer to the full array of FSI data structures. Pass NULL if disabled.
Returns
PetscErrorCode 0 on success.

Definition at line 22 of file implicitsolvers.c.

23{
24 PetscErrorCode ierr;
25 // --- Context Acquisition ---
26 // Get the master simulation context from the first block's UserCtx.
27 // This is the bridge to access all former global variables.
28 SimCtx *simCtx = user[0].simCtx;
29 PetscReal dt = simCtx->dt;
30 PetscReal st = simCtx->st;
31 PetscInt istage;
32 PetscReal alfa[] = {0.25, 1.0/3.0, 0.5, 1.0};
33
34 PetscFunctionBeginUser;
35
36 PROFILE_FUNCTION_BEGIN;
37
38 LOG_ALLOW(GLOBAL, LOG_INFO, "Executing explicit momentum solver (Runge-Kutta) for %d block(s).\n",simCtx->block_number);
39
40 // --- 1. Pre-Loop Initialization (Legacy Logic) ---
41 // This block prepares boundary conditions and allocates the RHS vector for all blocks
42 // before the main RK loop begins. This logic is preserved from the original code.
43 LOG_ALLOW(GLOBAL, LOG_DEBUG, "Preparing all blocks for Runge-Kutta solve...\n");
44 for (PetscInt bi = 0; bi < simCtx->block_number; bi++) {
45 ierr = InflowFlux(&user[bi]); CHKERRQ(ierr);
46 ierr = OutflowFlux(&user[bi]); CHKERRQ(ierr);
47 ierr = FormBCS(&user[bi]); CHKERRQ(ierr);
48
49 /*
50 // Immersed boundary interpolation (if enabled)
51 if (simCtx->immersed) {
52 LOG_ALLOW(LOCAL, LOG_DEBUG, " Performing pre-RK IBM interpolation for block %d.\n", bi);
53 for (PetscInt ibi = 0; ibi < simCtx->NumberOfBodies; ibi++) {
54 // The 'ibm' and 'fsi' pointers are passed directly from FlowSolver
55 ierr = ibm_interpolation_advanced(&user[bi], &ibm[ibi], ibi, 1); CHKERRQ(ierr);
56 }
57 }
58 */
59
60 // Allocate the persistent RHS vector for this block's context
61 ierr = VecDuplicate(user[bi].Ucont, &user[bi].Rhs); CHKERRQ(ierr);
62 }
63 LOG_ALLOW(GLOBAL, LOG_DEBUG, "Pre-loop initialization complete.\n");
64
65 // --- 2. Main Runge-Kutta Loop ---
66 // The legacy code had an outer `pseudot` loop that only ran once. We preserve it.
67 for (PetscInt pseudot = 0; pseudot < 1; pseudot++) {
68 // Loop over each block to perform the RK stages
69 for (PetscInt bi = 0; bi < simCtx->block_number; bi++) {
70 for (istage = 0; istage < 4; istage++) {
71 LOG_ALLOW(LOCAL, LOG_DEBUG, " Block %d, RK Stage %d (alpha=%.4f)...\n", bi, istage, alfa[istage]);
72
73 // a. Calculate the Right-Hand Side (RHS) of the momentum equation.
74 ierr = ComputeRHS(&user[bi], user[bi].Rhs); CHKERRQ(ierr);
75
76 // b. Advance Ucont to the next intermediate stage using the RK coefficient.
77 // Ucont_new = Ucont_old + alpha * dt * RHS
78 ierr = VecWAXPY(user[bi].Ucont, alfa[istage] * dt * st, user[bi].Rhs, user[bi].Ucont_o); CHKERRQ(ierr);
79
80 // c. Update local ghost values for the new intermediate Ucont.
81 ierr = DMGlobalToLocalBegin(user[bi].fda, user[bi].Ucont, INSERT_VALUES, user[bi].lUcont); CHKERRQ(ierr);
82 ierr = DMGlobalToLocalEnd(user[bi].fda, user[bi].Ucont, INSERT_VALUES, user[bi].lUcont); CHKERRQ(ierr);
83
84 // d. Re-apply boundary conditions for the new intermediate velocity.
85 // This is crucial for the stability and accuracy of the multi-stage scheme.
86 ierr = InflowFlux(&user[bi]); CHKERRQ(ierr);
87 ierr = OutflowFlux(&user[bi]); CHKERRQ(ierr);
88 ierr = FormBCS(&user[bi]); CHKERRQ(ierr);
89
90 } // End of RK stages for one block
91
92 /*
93 // Final IBM Interpolation for the block (if enabled)
94 if (simCtx->immersed) {
95 LOG_ALLOW(LOCAL, LOG_DEBUG, " Performing post-RK IBM interpolation for block %d.\n", bi);
96 for (PetscInt ibi = 0; ibi < simCtx->NumberOfBodies; ibi++) {
97 ierr = ibm_interpolation_advanced(&user[bi], &ibm[ibi], ibi, 1); CHKERRQ(ierr);
98 }
99 }
100 */
101
102 } // End loop over blocks
103
104 // --- 3. Inter-Block Communication (Legacy Logic) ---
105 // This is called after all blocks have completed their RK stages.
106 if (simCtx->block_number > 1) {
107 // LOG_ALLOW(GLOBAL, LOG_DEBUG, "Updating multi-block interfaces after RK stages.\n");
108 // ierr = Block_Interface_U(user); CHKERRQ(ierr);
109 }
110
111 } // End of pseudo-time loop
112
113 // --- 4. Cleanup ---
114 // Destroy the RHS vectors that were created at the start of this function.
115 for (PetscInt bi = 0; bi < simCtx->block_number; bi++) {
116 ierr = VecDestroy(&user[bi].Rhs); CHKERRQ(ierr);
117 }
118
119 LOG_ALLOW(GLOBAL, LOG_INFO, "Runge-Kutta solve completed for all blocks.\n");
120
121 PROFILE_FUNCTION_END;
122
123 PetscFunctionReturn(0);
124}
OutflowFlux
PetscErrorCode OutflowFlux(UserCtx *user)
Calculates the total outgoing flux through all OUTLET faces for reporting.
Definition Boundaries.c:932
InflowFlux
PetscErrorCode InflowFlux(UserCtx *user)
Applies inlet boundary conditions based on the modern BC handling system.
Definition Boundaries.c:706
FormBCS
PetscErrorCode FormBCS(UserCtx *user)
Definition Boundaries.c:1022
LOCAL
#define LOCAL
Logging scope definitions for controlling message output.
Definition logging.h:46
GLOBAL
#define GLOBAL
Scope for global logging across all processes.
Definition logging.h:47
LOG_ALLOW
#define LOG_ALLOW(scope, level, fmt,...)
Logging macro that checks both the log level and whether the calling function is in the allowed-funct...
Definition logging.h:201
PROFILE_FUNCTION_END
#define PROFILE_FUNCTION_END
Marks the end of a profiled code block.
Definition logging.h:740
LOG_INFO
@ LOG_INFO
Informational messages about program execution.
Definition logging.h:32
LOG_DEBUG
@ LOG_DEBUG
Detailed debugging information.
Definition logging.h:33
PROFILE_FUNCTION_BEGIN
#define PROFILE_FUNCTION_BEGIN
Marks the beginning of a profiled code block (typically a function).
Definition logging.h:731
ComputeRHS
PetscErrorCode ComputeRHS(UserCtx *user, Vec Rhs)
Computes the Right-Hand Side (RHS) of the momentum equations.
Definition rhs.c:1159
SimCtx::block_number
PetscInt block_number
Definition variables.h:593
UserCtx::simCtx
SimCtx * simCtx
Back-pointer to the master simulation context.
Definition variables.h:664
SimCtx::st
PetscReal st
Definition variables.h:581
SimCtx::dt
PetscReal dt
Definition variables.h:552
SimCtx
The master context for the entire simulation.
Definition variables.h:538
Here is the call graph for this function:
Here is the caller graph for this function:

◆ ImplicitMomentumSolver()

PetscErrorCode ImplicitMomentumSolver ( UserCtx *  user,
IBMNodes *  ibm,
FSInfo *  fsi 
)
extern

Advances the momentum equations using an implicit scheme.

Parameters
userArray of UserCtx structs for all blocks.
ibm(Optional) Pointer to IBM data. Pass NULL if disabled.
fsi(Optional) Pointer to FSI data. Pass NULL if disabled.
Returns
PetscErrorCode 0 on success.