1D space

[festung]

1D space

[reuterbal]

Hi Dieter,

it is definitely possible. Have a look at core/generateGridData1D.m, which allows to generate a "1D grid". 1D modal basis functions and their gradients are in core/phi1D.m and core/gradPhi1D.m, respectively.
Most of the assembly routines are specific to 2D problems (since there is a lot more to be done than in 1D) but you might find some for 1D (which always have 1D in their name).

You can also look at our latest preprint which uses a 1D grid to discretize the free surface elevation and a 2D grid for the velocity unknowns for a free surface flow problem in a vertical slice. This might give you some more insight into the 1D capabilities of FESTUNG.

Feel free to come back with any questions.

Best regards,
Balthasar

[dietergrass]

Hi Balthasar, thanks a lot. My idea is to integrate the algorithm in a continuation framework for age specific optimal control problems. So I will definitely come back to you with questions :) Best Dieter

…

________________________________ From: Balthasar Reuter <[email protected]> Sent: Friday, October 11, 2019 9:56 AM To: FESTUNG/FESTUNG Cc: Grass, Dieter; Author Subject: Re: [FESTUNG/FESTUNG] 1D space (#5) Hi Dieter, it is definitely possible. Have a look at core/generateGridData1D.m<https://github.com/FESTUNG/FESTUNG/blob/master/core/generateGridData1D.m>, which allows to generate a "1D grid". 1D modal basis functions and their gradients are in core/phi1D.m<https://github.com/FESTUNG/FESTUNG/blob/master/core/phi1D.m> and core/gradPhi1D.m<https://github.com/FESTUNG/FESTUNG/blob/master/core/gradPhi1D.m>, respectively. Most of the assembly routines are specific to 2D problems (since there is a lot more to be done than in 1D) but you might find some for 1D (which always have 1D in their name). You can also look at our latest preprint<https://arxiv.org/abs/1806.03908> which uses a 1D grid to discretize the free surface elevation and a 2D grid for the velocity unknowns for a free surface flow problem in a vertical slice. This might give you some more insight into the 1D capabilities of FESTUNG. Feel free to come back with any questions. Best regards, Balthasar - You are receiving this because you authored the thread. Reply to this email directly, view it on GitHub<#5?email_source=notifications&email_token=ANOQLUFJTAMGLJ2Q7YOJX5DQOAWTFA5CNFSM4I7V2XXKYY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOEA7F2JQ#issuecomment-540957990>, or unsubscribe<https://github.com/notifications/unsubscribe-auth/ANOQLUHTSQXKAR2HLCU25S3QOAWTFANCNFSM4I7V2XXA>.

[dietergrass]

Hi Balthasar, is it possible to adapt the code in FESTUNG considering nonlocal (integral) boundary conditions and for a higher dimensional state vector? Thanks Dieter

…

________________________________ From: Balthasar Reuter <[email protected]> Sent: Friday, October 11, 2019 9:56 AM To: FESTUNG/FESTUNG Cc: Grass, Dieter; Author Subject: Re: [FESTUNG/FESTUNG] 1D space (#5) Hi Dieter, it is definitely possible. Have a look at core/generateGridData1D.m<https://github.com/FESTUNG/FESTUNG/blob/master/core/generateGridData1D.m>, which allows to generate a "1D grid". 1D modal basis functions and their gradients are in core/phi1D.m<https://github.com/FESTUNG/FESTUNG/blob/master/core/phi1D.m> and core/gradPhi1D.m<https://github.com/FESTUNG/FESTUNG/blob/master/core/gradPhi1D.m>, respectively. Most of the assembly routines are specific to 2D problems (since there is a lot more to be done than in 1D) but you might find some for 1D (which always have 1D in their name). You can also look at our latest preprint<https://arxiv.org/abs/1806.03908> which uses a 1D grid to discretize the free surface elevation and a 2D grid for the velocity unknowns for a free surface flow problem in a vertical slice. This might give you some more insight into the 1D capabilities of FESTUNG. Feel free to come back with any questions. Best regards, Balthasar - You are receiving this because you authored the thread. Reply to this email directly, view it on GitHub<#5?email_source=notifications&email_token=ANOQLUFJTAMGLJ2Q7YOJX5DQOAWTFA5CNFSM4I7V2XXKYY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOEA7F2JQ#issuecomment-540957990>, or unsubscribe<https://github.com/notifications/unsubscribe-auth/ANOQLUHTSQXKAR2HLCU25S3QOAWTFANCNFSM4I7V2XXA>.

[reuterbal]

Hi,

adding integral boundary conditions should be fairly easy to do in 1D. You will have to take care of the evaluation yourself (e.g., using the provided quadrature rules) and then put it in the right place of your linear system.

What do you mean by higher dimensional state vector? Basis functions in 1D are available up to local polynomial degree 5, which is already reasonable higher-order.

HTH,
Balthasar

[dietergrass]

Hi, can u be chosen in R^m, e.g. for d/dt u+d/dxf(u)=g(u) I think the answer is yes, what I read, but I wanted to be sure. Dieter

…

________________________________ From: Balthasar Reuter <[email protected]> Sent: Tuesday, October 15, 2019 10:03 AM To: FESTUNG/FESTUNG Cc: Grass, Dieter; Author Subject: Re: [FESTUNG/FESTUNG] 1D space (#5) Hi, adding integral boundary conditions should be fairly easy to do in 1D. You will have to take care of the evaluation yourself (e.g., using the provided quadrature rules) and then put it in the right place of your linear system. What do you mean by higher dimensional state vector? Basis functions in 1D are available up to local polynomial degree 5, which is already reasonable higher-order. HTH, Balthasar - You are receiving this because you authored the thread. Reply to this email directly, view it on GitHub<#5?email_source=notifications&email_token=ANOQLUCJGSIWH4JOQSLLLPTQOV2O7A5CNFSM4I7V2XXKYY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOEBH2NZI#issuecomment-542090981>, or unsubscribe<https://github.com/notifications/unsubscribe-auth/ANOQLUG7DVGHSRRF6I23SL3QOV2O7ANCNFSM4I7V2XXA>.

[reuterbal]

Meaning, you have a system of 1D equations? That is definitely possible.

[dietergrass]

yes, thanks

…

________________________________ From: Balthasar Reuter <[email protected]> Sent: Tuesday, October 15, 2019 10:24 AM To: FESTUNG/FESTUNG Cc: Grass, Dieter; Author Subject: Re: [FESTUNG/FESTUNG] 1D space (#5) Meaning, you have a system of 1D equations? That is definitely possible. - You are receiving this because you authored the thread. Reply to this email directly, view it on GitHub<#5?email_source=notifications&email_token=ANOQLUANAPH43R7S45DPSPLQOV423A5CNFSM4I7V2XXKYY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOEBH4IMI#issuecomment-542098481>, or unsubscribe<https://github.com/notifications/unsubscribe-auth/ANOQLUBTVQCMU4CSATJHUG3QOV423ANCNFSM4I7V2XXA>.

[dietergrass]

Hi, again and sorry for bothering you. But could you help me / give me an advice for setting up a test example, like the telegrapher equations in the attached file. Thanks a lot Dieter PS: we could also communicate in German.

…

________________________________ From: Balthasar Reuter <[email protected]> Sent: Tuesday, October 15, 2019 10:24 AM To: FESTUNG/FESTUNG Cc: Grass, Dieter; Author Subject: Re: [FESTUNG/FESTUNG] 1D space (#5) Meaning, you have a system of 1D equations? That is definitely possible. - You are receiving this because you authored the thread. Reply to this email directly, view it on GitHub<#5?email_source=notifications&email_token=ANOQLUANAPH43R7S45DPSPLQOV423A5CNFSM4I7V2XXKYY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOEBH4IMI#issuecomment-542098481>, or unsubscribe<https://github.com/notifications/unsubscribe-auth/ANOQLUBTVQCMU4CSATJHUG3QOV423ANCNFSM4I7V2XXA>.