Re: [petsc-users] TS_SSP implementation for co-dependent variables

2019-10-23 Thread Jed Brown via petsc-users
Manuel Valera writes: > Yes, all of that sounds correct to me, > > No I haven't tried embedding the column integral into the RHS, right now I > am unable to think how to do this without the solution of the previous > intermediate stage. Any ideas are welcome, Do you have some technical notes on

Re: [petsc-users] TS_SSP implementation for co-dependent variables

2019-10-09 Thread Manuel Valera via petsc-users
Yes, all of that sounds correct to me, No I haven't tried embedding the column integral into the RHS, right now I am unable to think how to do this without the solution of the previous intermediate stage. Any ideas are welcome, Thanks, On Wed, Oct 9, 2019 at 4:18 PM Jed Brown wrote: > Manuel V

Re: [petsc-users] TS_SSP implementation for co-dependent variables

2019-10-09 Thread Jed Brown via petsc-users
Manuel Valera writes: > Sorry, I don't follow this last email, my spatial discretization is fixed, > the problem is caused by the choice of vertical coordinate, in this case > sigma, that calls for an integration of the hydrostatic pressure to correct > for the right velocities. Ah, fine. To ph

Re: [petsc-users] TS_SSP implementation for co-dependent variables

2019-10-09 Thread Manuel Valera via petsc-users
Sorry, I don't follow this last email, my spatial discretization is fixed, the problem is caused by the choice of vertical coordinate, in this case sigma, that calls for an integration of the hydrostatic pressure to correct for the right velocities. I had RK3 working before and SSP is much more s

Re: [petsc-users] TS_SSP implementation for co-dependent variables

2019-10-09 Thread Jed Brown via petsc-users
Is it a problem with the spatial discretization or with the time discretization that you've been using thus far? (This sort of problem can occur for either reason.) Note that an SSP method is merely "preserving" -- the spatial discretization needs to be strongly stable for an SSP method to preser

Re: [petsc-users] TS_SSP implementation for co-dependent variables

2019-10-09 Thread Manuel Valera via petsc-users
To correct for the deformation of the sigma coordinate grid... without this correction the velocity become unphysical in the zones of high slope of the grid. This is very specific of our model and probably will be solved by updating the equations transformation, but that's not nearly close to happe

Re: [petsc-users] TS_SSP implementation for co-dependent variables

2019-10-09 Thread Manuel Valera via petsc-users
Thanks, My time integration schemes are all explicit, sorry if this a very atypical setup. This is similar to the barotropic splitting but not exactly, we don't have free surface in the model, this is only to correct for sigma coordinates deformations in the velocity field. >From how i see it thi

Re: [petsc-users] TS_SSP implementation for co-dependent variables

2019-10-09 Thread Jed Brown via petsc-users
Manuel Valera writes: > Thanks, > > My time integration schemes are all explicit, sorry if this a very atypical > setup. This is similar to the barotropic splitting but not exactly, we > don't have free surface in the model, this is only to correct for sigma > coordinates deformations in the velo

Re: [petsc-users] TS_SSP implementation for co-dependent variables

2019-10-09 Thread Jed Brown via petsc-users
Manuel Valera writes: > Thanks for the answer, I will read the mentioned example, but to clarify > for Barry I will schematize the process: > > At time n, the program need to do all of these at once: > >1. Solve T as a function of u,v,w >2. Solve S as a function of u,v,w >3. Solve rho

Re: [petsc-users] TS_SSP implementation for co-dependent variables

2019-10-09 Thread Manuel Valera via petsc-users
Thanks for the answer, I will read the mentioned example, but to clarify for Barry I will schematize the process: At time n, the program need to do all of these at once: 1. Solve T as a function of u,v,w 2. Solve S as a function of u,v,w 3. Solve rho density as a function of T,S 4. De

Re: [petsc-users] TS_SSP implementation for co-dependent variables

2019-10-09 Thread Jed Brown via petsc-users
Manuel Valera via petsc-users writes: > Hello, > > I have a set of equations which are co-dependent when integrating in time, > this means the velocities u,v,w need a component from the Temperature and > Salinity integration at the same intermediate step. Same for Temperature > and Salinity, whic

Re: [petsc-users] TS_SSP implementation for co-dependent variables

2019-10-09 Thread Smith, Barry F. via petsc-users
How are you integrating the two sets of variables now? Would it be possible to just integrate all five variables together and not have any of these difficulties? Barry > On Oct 9, 2019, at 6:03 PM, Manuel Valera via petsc-users > wrote: > > Hello, > > I have a set of equation

[petsc-users] TS_SSP implementation for co-dependent variables

2019-10-09 Thread Manuel Valera via petsc-users
Hello, I have a set of equations which are co-dependent when integrating in time, this means the velocities u,v,w need a component from the Temperature and Salinity integration at the same intermediate step. Same for Temperature and Salinity, which need the current velocities (at the intermediate