How can we implement a fully implicit linear source term? -- CFD Online Discussion Forums

February 8, 2021, 21:15	How can we implement a fully implicit linear source term?	#1
Rango New Member Join Date: Jun 2020 Location: UK Posts: 22 Rep Power: 5	Hi, I would like to solve the following equation: $-\nabla.(\kappa\nabla\phi) + s = 0$ , to obtain $\phi$ field for a 3-layered simple mesh setup. Here is a linear source term which is zero for 'upper' and 'lower' layers. For the 'mid' layer, however, we have: $s_{mid-layer}=f(\phi_{upper-layer}, \phi_{lower-layer})$ , where is a linear function. Example: for block 3 we have: $s_3=a_0\phi_1+a_1\phi_2+a_2$ where are constant properties of block 3. Is there a way to solve this equation for the whole mesh while treating the mid-layer as a buffer? Please note that the role of 'mid-layer' is nothing but coupling $\phi$ values in 'upper-layer' and 'lower-layer' through . I would appreciate your help on this. Cheers Last edited by Rango; February 21, 2021 at 18:31.

Posting Rules
You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is Off Trackbacks are Off Pingbacks are On Refbacks are On Forum Rules

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Custom Thermophysical Properties	wsmith02	OpenFOAM	4	June 1, 2023 14:30
[OpenFOAM.com] swak4foam compiling issues on a cluster	saj216	OpenFOAM Installation	5	January 17, 2023 16:05
[OpenFOAM.org] Error creating ParaView-4.1.0 OpenFOAM 2.3.0	tlcoons	OpenFOAM Installation	13	April 20, 2016 17:34
[foam-extend.org] problem when installing foam-extend-1.6	Thomas pan	OpenFOAM Installation	7	September 9, 2015 21:53
[Other] Adding solvers from DensityBasedTurbo to foam-extend 3.0	Seroga	OpenFOAM Community Contributions	9	June 12, 2015 17:18

All times are GMT -4. The time now is 11:48.