Hi Martin, So, that means effectivePressure[eIdx] should be multiplied by alpha in model.hh:408 to get the correct total and principal stresses in the output files. Correct?
best regards, Edscott -----Mensaje original----- De: Martin Beck [mailto:[email protected]] Enviado el: lunes, 11 de marzo de 2019 04:10 a. m. Para: DuMuX User Forum; Ed Scott Wilson Garcia CC: Manuel Coronado Gallardo Asunto: Re: [DuMuX] Using Biot's Modulus and coefficient in el2p model Hi Edscott, yes, the model assumes alpha = 1. If you want to adapt this, you have to change it in the elementvolumevariables.hh, but also for the assembly of the momentum balance equation in localoperator.hh:346 (- alpha * peff instead of - peff). The post-processing evaluation of the total stresses and the principal stresses in model.hh also contains this assumption. Cheers, Martin On 3/8/19 9:02 PM, Ed Scott Wilson Garcia wrote: > We are now introducing Biot's modulus into the el2p problem and using this > to calculate Skempton's coefficient instead of the other way around (Biot's > modulus = 1/c0_). > In order to obtain the undrained bulk modulus we also need Biot's coefficient > (alpha), so we also provide that in the input file. > > My question is as follows. > To make the use of Biot's coefficient complete, would we also have to modify > the el2p model? > Specifically, multiply the divergence of u with alpha in > elementvolumevariables.hh:259? > > best regards, > > Edscott > _______________________________________________ > Dumux mailing list > [email protected] > https://listserv.uni-stuttgart.de/mailman/listinfo/dumux -- ******************************************************* Martin Beck Institute for Modelling Hydraulic and Environmental Systems Department of Hydromechanics and Modelling of Hydrosystems Pfaffenwaldring 61 D-70569 Stuttgart Tel.: (+49) 711/ 685-64899 [email protected] ******************************************************* _______________________________________________ Dumux mailing list [email protected] https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
