Diffusive Mixing (1D <--> CFD) #50
Conversation
The resolution for the calculation of diffusion using an analytical solution including a fourier series is renamed to numFourierTerms and set as a private variable in the DiffusionMixingModel. It is set to ```int numFourierTerms = 100;``` It can be adapted using ```resetNumFourierTerms(int newNumFourierTerms)``` This solves Issue #43
MariaEmmerich
added
feature
|
@micheltakken I have added the connection between the abstract diffusive mixing to and from the CFD Modules.
For both cases we adhere to the law of mass conservation. The connection "on the CFD side" is still missing. Once we have included the missing parts we can start to test. |
src/simulation/MixingModels.hh
Outdated
| // auto& section = outflowDistributions.at(channel->getId()).front(); // there is only one channel connected to a node at a CFD hybrid node connection | ||
| // // std::vector<T> inflowSectionIntoCfd = section.concentrationField; | ||
| // } | ||
| int resolutionIntoCfd = 20; // TODO this should be a parameter |
There was a problem hiding this comment.
Here the lattice resolution is hard coded rather than passed from the simulator.
There was a problem hiding this comment.
The resolution values are stored in olbContinuous.h or essContinuous.h and can technically be obtained from a cfdSimulator with cfdSimulator->getResolution(). This is not yet implemented. I will take care of this.
TODO
There was a problem hiding this comment.
I added the getResolution(nodeId) function and a TODO in the olbMixing.hh
…ntration profiles.
This only works if the resolution, i.e., the number of lattice points at the node, does not change during the simulation (at different timesteps)
This PR updates
The connection of the CFD Module and 1D diffusive mixing, concentrations are passed as concentration fields to the CFD Module and concentration fields are translated as linear segments as input for the 1D diffusive mixing.
PR Checklist