ations, I think you
> want to read this paper:
> https://www.math.colostate.edu/~bangerth/publications/2017-boussinesq.pdf
> Specifically, section 3.3 talks about this situation, and it shows in
> figures
> 2 and 4 that you do get spikes in the pressure.
>
> Best
>
This is exactly what I was thinking off. Thanks.
On Friday, January 3, 2025 at 7:44:48 PM UTC+3 Wolfgang Bangerth wrote:
>
> Najwa:
>
> > I know that merge triangulation in deal.ii will not work properly if
> there are
> > hanging nodes ( unmatching of the meshes at the interface).
> >
> > My
Dear Abbas,
I appreciated your answer as I did not know about step 46. Thank you.
I already merged cells to remove the hanging node, I think there is no
other solution for it now. Some studies are out there about gluing
different faces that do not match. Not sure if that would help anyone
sear
Dear Abbas,
Thank you for your reply. I quickly looked at step 46 and did not notice
that they had the solution to my question. I will read it more carefully
tomorrow.
I forgot to say that I know that a trivial solution is to refine the
corser side. However, my question is, can this be avoid
Dear all,
I have a quick question here.
I know that merge triangulation in deal.ii will not work properly if there
are hanging nodes ( unmatching of the meshes at the interface).
My question is, what if I have to have this hanging node? Is there a way I
can work it out, or are there any other
s,
Najwa
On Wednesday, December 4, 2024 at 8:37:04 AM UTC+3 Najwa Alshehri wrote:
> Thank you for your response.
>
> For those who might face something similar, I found that the saved file
> does not have the solution saved and the last thing that was saved is the
> dof_handler (as
the benefit of
others if I find a solution to it.
Best,
Najwa
On Wednesday, December 4, 2024 at 7:57:59 AM UTC+3 Wolfgang Bangerth wrote:
> On 12/3/24 02:32, Najwa Alshehri wrote:
> >
> > The exact same steps I used for elliptic problem with FE_Q and solution
> that
>
as
> small a program as possible that still shows the issue. It doesn't have
> to do anything useful, for example you can replace assemble_system() and
> solve() by just setting the vector to ones. Can you come up with such a
> minimal program?
>
> Best
>
problems is to create as
> small a program as possible that still shows the issue. It doesn't have
> to do anything useful, for example you can replace assemble_system() and
> solve() by just setting the vector to ones. Can you come up with such a
> minimal program?
>
Dear All,
Thank you for your always support.
I am trying to serialize a Stokes solution in 2d ( dim for velocity +1 for
pressure) to use it as a reference solution. Please stop me here if there
is a better way to do this.
Unofortunatly, for some reason that I do not know, the saved_solution is
Thanks
On Monday, September 16, 2024 at 5:14:03 PM UTC+3 Wolfgang Bangerth wrote:
> On 9/16/24 05:38, Najwa Alshehri wrote:
> >
> > It could be an issue related to having one of the meshes circle
> approximated
> > by a polygon. I have a guess that the error o
ersection calculation.
If those notes would help anyhow. Thank you all for your support.
Best,
Najwa
On Saturday, September 14, 2024 at 5:55:10 PM UTC+3 Najwa Alshehri wrote:
> Dear developers and users,
>
> I am currently working on a problem involving the coupling of multiple
> physics.
Wolfgang:
Thank you for bringing this to my attention. I would definitely love to
share the code in the code gallery for the benefit of others. In this case,
I will create a code for a simpler problem, such as the Stokes problem or
mixed Laplace problem, to illustrate the concept of the inf-sup
you did for the mass matrix, but replacing the inverse
> operator done with umfpack with one done with CG.
>
> Luca
>
> Il giorno 2 giu 2024, alle ore 09:42, Najwa Alshehri
> ha scritto:
>
> Wolfgang and all,
>
> I have a positive definite matrix M, and its inverse
Wolfgang Bangerth wrote:
> On 6/1/24 14:49, Najwa Alshehri wrote:
> >
> > I decided to solve the exact problem directly, namely AA x = \lambda M
> x. To
> > achieve this, I computed the inverse of the matrix AA= Bt * A^inv * B
> using a
> > Conjugate Gradient (CG) s
cal yields different
results, I would greatly appreciate your input on this matter.
Best,
Najwa
On Friday, May 31, 2024 at 11:41:27 AM UTC+3 Najwa Alshehri wrote:
> Dear Wolfgang:
>
> Thank you for your answer. I have checked the matrices and I printed the
> case for a 2 by 2 mes
x = \lambda M x
>
>
> Best,
> Matthias
>
>
> On Thu, May 30, 2024, at 10:55 CDT, Najwa Alshehri
> wrote:
>
> > Dear all,
> >
> >
> > Thank you for your always support and help.
> >
> >
> > I am encountering an issue while solvi
Dear all,
Thank you for your always support and help.
I am encountering an issue while solving an eigenvalue problem related to
the computation of the discrete inf-sup constant for a saddle point
problem. Specifically, I am solving the following system:
Bt A^-1 B eigenvector = eigenvalue
Dear developers and users,
I have two meshes one is immersed in the other. I wanted to find the
intersection between the two meshes, so I used the following function.
NonMatching::compute_intersection(omega_grid_tools_cache,
omega2_grid_tools_ca
UTC+3 Najwa Alshehri wrote:
> Hello Bruno,
>
> I have the solution serialized with the old version. Do you think this is
> the issue? I will do everything with the new version today and update you
> if it works.
>
> Best,
> Najwa
>
> On Monday, November 6, 20
Hello Bruno,
I have the solution serialized with the old version. Do you think this is
the issue? I will do everything with the new version today and update you
if it works.
Best,
Najwa
On Monday, November 6, 2023 at 5:33:11 PM UTC+3 bruno.t...@gmail.com wrote:
> Najwa,
>
> Are you serializin
Dear developers and group members,
I was using an old version of Dealii where I was serializing some solutions
and deserializing them where needed. In the code where I deserialize the
solution I used to use the following headers:
#include
#include
#include
#include
#include
and I used
r. I will try to write a patch to the deal.II
>>>> sources that implement the missing function. This might require some time.
>>>> I will come back here if needed.
>>>>
>>>>
>>>> Thank you all for your quick answers,
>>>>
answers,
>> Najwa
>>
>> On Wednesday, July 26, 2023 at 11:41:15 PM UTC+3 Wolfgang Bangerth wrote:
>>
>>> On 7/26/23 12:17, Najwa Alshehri wrote:
>>> >
>>> > Regarding the function "VectorFucntionFromScalarFunctionObject," I
>>&g
Wolfgang Bangerth wrote:
> On 7/26/23 12:17, Najwa Alshehri wrote:
> >
> > Regarding the function "VectorFucntionFromScalarFunctionObject," I have
> > observed that it works for computing the L2 norm of the error. However,
> when I
> > attempted to compute the
"FeFeildFunction."
Could you kindly advise me on the simplest way to accomplish this?
On a side note, although "VectorFucntionFromScalarFunctionObject" works
with the L2 norm of the error, it is quite computationally intensive.
Best regards,
Najwa
On Sunday, July 23, 2023 at 4
y_fun([&](const auto &p){ return
> fe_field_func.value(p);}, 0, 2);
>
> Luca
>
> Il giorno 21 lug 2023, alle ore 22:05, Najwa Alshehri
> ha scritto:
>
> Dear Daniel, your answer makes sense, Finally it worked.
>
>
> Dear Luca, thank you for your answer
ion to allocate two components:
>
> ExactSolution2::ExactSolution2( … )
> : Functions::Function(2)
>
> Best,
> Luca.
>
>
> > On 20 Jul 2023, at 12:15, Najwa Alshehri wrote:
> >
> > Hello again,
> >
> > I have a follow-up question. Does th
means that the mask is not really picking up the
first component of the solution u_omega2. Do you have any suggestions to
fix this issue?
Appreciate your help,
Najwa
On Thursday, July 20, 2023 at 11:48:04 AM UTC+3 Najwa Alshehri wrote:
> Thank you Daniel for the clear quick answer. I wil
> for all vector components uniformly.
>
> Best,
> Daniel
>
> On Wed, Jul 19, 2023 at 8:01 AM Najwa Alshehri
> wrote:
>
>> Dear group members,
>>
>> I have one question,
>>
>> I am trying to calculate the L2 norm of the error of a solution.
Dear group members,
I have one question,
I am trying to calculate the L2 norm of the error of a solution. In
particular, I have a *vector of the solution* u_omega and a *FeSystem* with
two fes.
I am interested in computing the L2 norm of the error related to the first
fe using *integrate_dif
31 matches
Mail list logo
