Pavel Tomin
Pavel Tomin
GEOS formulation allows saturation values to go above 1 (not sure about below zero). The computed saturation values are used as is to compute rel perms that might be not...
Stress initialization does not properly work in sequential mode when mechanics-only solver is used to perform initialization. The reason is zero bulk density. Initialization works ok if coupled (flow+mechanics) is...
Case: [coupled_benchmark (2).zip](https://github.com/GEOS-DEV/GEOS/files/12562800/coupled_benchmark.2.zip) Running with 2400 MPIs instead of 1200 MPIs leads to the following error: ``` Time: 0.00e+00 s, dt: 10000 s, Cycle: 0 Attempt: 0, ConfigurationIter: 0, NewtonIter:...
Case from here https://github.com/GEOS-DEV/MAELSTROM/tree/master/usecases/francois/syntheticMultiphasePoromechanics Modified maxDt = 1e8 as in [InputInitialized.txt](https://github.com/GEOS-DEV/GEOS/files/12562779/InputInitialized.txt) For some reason more partitioning along Z direction leads to increase in linear iterations: 
Case from @sohailwaziri, just documenting here to keep the track. Modified version of https://github.com/GEOS-DEV/MAELSTROM/tree/master/usecases/francois/HI24L/isothermalMultiphasePoromechanics Modifications: increase coupling strength by reducing K_dr of the reservoir solid, increased the time step size...
- [ ] use `ParallelVector` and remove re-implemented `dot`, `axpy`, `scale` - [ ] threading - [ ] MPI support - [ ] other acceleration methods (such as Quasi-Newton)
- [ ] check and if possible remove `serialPolicy` here https://github.com/GEOS-DEV/GEOS/blob/d814b2213c1d15a988ff56581ecd64a25bea9dd2/src/coreComponents/physicsSolvers/multiphysics/SinglePhasePoromechanicsConformingFractures.cpp#L483 - [ ] convert to c-array https://github.com/GEOS-DEV/GEOS/blob/d814b2213c1d15a988ff56581ecd64a25bea9dd2/src/coreComponents/physicsSolvers/multiphysics/SinglePhasePoromechanicsConformingFractures.cpp#L507 and check if all others should be using stackArrays
Example 1 (conforming): [PoroElastic_conformingFracture_2d_faultSlip](https://github.com/GEOS-DEV/GEOS/blob/develop/inputFiles/poromechanicsFractures/PoroElastic_conformingFracture_2d_faultSlip_benchmark.xml) switch to MGR: ``` ``` run and observe: ``` Time: 0.00e+00 s, dt: 100 s, Cycle: 0 Number of element for each fracture state: stick: 80...
Introduce configuration tolerance, the configuration is assumed converged if changed area is below the defined fraction of total area: `changedArea =1 mean effectively configuration check is disabled No default behavior...