On Thu, Aug 12, 2021 at 3:59 PM Renato Poli wrote:
> Thanks,
>
> As I understood, it just works without any major concern in residual and
> jacobian calculations, right?
> I must add the DirichletBoundary to the system dof_map and LibMesh handles
> the matrix and vectors so Petsc will receive the
Thanks,
As I understood, it just works without any major concern in residual and
jacobian calculations, right?
I must add the DirichletBoundary to the system dof_map and LibMesh handles
the matrix and vectors so Petsc will receive the constrained system.
Is that correct?
Thanks,
Renato
On Thu, A
Hello Renato,
I think example 8 in systems_of_equations should be
useful in this context.
Vikram Garg
On Thu, Aug 12, 2021 at 8:48 AM Renato Poli wrote:
> Hi all,
>
> I had no reply on this one, I am struggling with this design.
> I could not make this work, should I go f
Hi all,
I had no reply on this one, I am struggling with this design.
I could not make this work, should I go for the Penalty method?
Renato
On Tue, Aug 10, 2021 at 10:49 PM Renato Poli wrote:
> Hi,
>
> In linear systems, I was used to adding dirichlet constraints with
> dogmap.add_dirichlet_
Hi,
In linear systems, I was used to adding dirichlet constraints with
dogmap.add_dirichlet_constraints(). I needed to use
heterogeneously_constrain_matrix_and_vector after assembly.
What would be the way to do that in NonlinearImplicitSystem? I understand
that NR solvers solve for the increment,
> Are you using the homogeneous or the heterogeneous constraint
> functions when constraining your system of equations? Often we solve
> for a delta between two solutions (Newton steps, time steps with fixed
> BC values, etc) rather than directly solving for the constrained
> solution, so we have
On Sun, 7 Jan 2018, Renato Poli wrote:
I am struggling with a Dirichlet boundary condition setting the
boundary to zero, instead of the value I impose (3.e7). Previously, I
was using the penalty method, and it was working fine.
Is there any known corner condition where the DirichletBoundary wo
Hi
I am struggling with a Dirichlet boundary condition setting the
boundary to zero, instead of the value I impose (3.e7). Previously, I
was using the penalty method, and it was working fine.
Is there any known corner condition where the DirichletBoundary would
drive the variable to zero? Any ide
Roy, thanks for the detailed reply. I look forward to browsing through the
linked PDF.
Alex
> On Jan 3, 2017, at 8:51 AM, Roy Stogner wrote:
>
>
>> On Mon, 2 Jan 2017, Alexander Lindsay wrote:
>>
>> While thinking about possibly implementing H(curl) projection sometime
>> in the future, I'v
On Mon, 2 Jan 2017, Alexander Lindsay wrote:
> While thinking about possibly implementing H(curl) projection sometime
> in the future, I've been trying to familiarize myself with libmesh's
> current projection implementations. I have a basic question: in
> dof_map_constraints.C how come projectio
Hi all,
While thinking about possibly implementing H(curl) projection sometime
in the future, I've been trying to familiarize myself with libmesh's
current projection implementations. I have a basic question: in
dof_map_constraints.C how come projections for side degrees of freedom
(roughly li
Thanks to David and John for your hints. I need to further digest it.
--Junchao Zhang
On Wed, Jan 13, 2016 at 4:33 PM, John Peterson wrote:
>
>
> On Wed, Jan 13, 2016 at 3:23 PM, Junchao Zhang
> wrote:
>
>> Hello,
>> I tried to follow the textbook to write weak forms of PDEs with Dirichlet
>>
On Wed, Jan 13, 2016 at 3:23 PM, Junchao Zhang
wrote:
> Hello,
> I tried to follow the textbook to write weak forms of PDEs with Dirichlet
> boundary conditions like some libmesh examples. The textbook usually lets u
> = 0 on boundary, making it easy to understand why the boundary term
> vanishes
Hello,
I tried to follow the textbook to write weak forms of PDEs with Dirichlet
boundary conditions like some libmesh examples. The textbook usually lets u
= 0 on boundary, making it easy to understand why the boundary term
vanishes in weak form. But when u is not zero on boundary, I found libmesh
Hi,
I'm trying to impose tw types of Dirichlet condition, one is homogeneous
and in both variables and the other is nonhomogeneus and only in the second
variable. It seems to work out in the initial projected solution but when I
solve iteratively the problem the last solution does not satisfy the
12:43:21 -0500
> From: [email protected]
> To: [email protected]
> Subject: RE: [Libmesh-users] dirichlet boundaries not being applied.
>
>
> On Wed, 28 May 2014, subramanya sadasiva wrote:
>
> > I can't find print constraints, but print_info from the dof
On Wed, 28 May 2014, subramanya sadasiva wrote:
> I am trying to use dirichlet boundaries to apply a dirichlet
> boundary on an internal boundary. This doesn't seem to be working
> and I seem to keep getting 0 . I have heterogenously constrained the
> element matrix and vector, but the dirichlet
I am trying to use dirichlet boundaries to apply a dirichlet boundary on an
internal boundary. This doesn't seem to be working and I seem to keep getting 0
. I have heterogenously constrained the element matrix and vector, but the
dirichlet boundary conditions just don't seem to be called. Any
boundary side (y=y_max) is supposed to deform evenly
>along both directions of the x-axis. However, in both figures, the
>deformation occurs mainly in one direction.
> Thanks very much for your help. Cheers, Dafang Date: Tue, 4 Mar 2014
> 18:47:51 -0500 From: David Knezevic
&g
along both directions of the x-axis. However, in both figures, the
deformation occurs mainly in one direction.
Thanks very much for your help. Cheers, Dafang Date: Tue, 4 Mar 2014
18:47:51 -0500 From: David Knezevic
Subject: Re: [Libmesh-users] Dirichlet Boundary Condition on Part of the
Hi Dafang,
I've done this a lot for elasticity problems, and haven't had any
problems with it.
If you can send through an example that illustrates what is happening
(e.g. a modified version of systems_of_equations_ex6), I'd be happy to
have a look.
David
On 03/04/2014 06:08 PM, Dafang Wang
Hi,
The Libmesh Examples tell how to use the DirichletBoundary class to
constrain all the variables in a system, but I am not sure how to apply
zero Dirichlet BC to PART of the field variables in a system. I am
solving a simple 3D elastic problem and I want to fix a boundary in the
y direction
Thanks, David.
I will be able to do that. I also have a patch coming for some work on
sensitivity analysis.
I figured that the enforce_constraints_exactly method is what provides the
capability for other systems for application of Dirichlet bcs. I should be
able to get the dof information from t
Hi Manav,
It would be great to update CondensedEigenSystem in the way you propose.
I originally implemented that class since I needed it for some of the
Reduced Basis code, and it pre-dates the DirichletBoundary stuff.
Current the user has to initialize the list of Dirichlet dofs in
CondensedE
Hi,
I am trying to understand the application of DirichletBoundary
constraints in the system solution.
For a system, say LinearImplicitSystem or NonlinearImplicitSystem, once
the user calls solve, at what stage (and which routine) in this function
call are the Dirichlet boundary condition
On Wed, 15 May 2013, Ataollah Mesgarnejad wrote:
> That's exactly is what I was wondering.
>
> As always many thanks for the help,
You're welcome! Glad I finally started to make sense.
One of these days I need to just get in the habit of writing actual
thorough documentation rather than expect
On May 15, 2013, at 5:44 PM, Roy Stogner wrote:
>
> On Wed, 15 May 2013, Ataollah Mesgarnejad wrote:
>
>> Unless I'm missing something, I don't think that would solve the problem.
>>
>> Maybe I'm not explaining my problem clear enough: I have parameters in
>> my_equation_system.parameters li
On Wed, 15 May 2013, Ataollah Mesgarnejad wrote:
> Unless I'm missing something, I don't think that would solve the problem.
>
> Maybe I'm not explaining my problem clear enough: I have parameters in
> my_equation_system.parameters like "boundary_type_01" with a value like
> "ramp" and "nounda
Roy,
Unless I'm missing something, I don't think that would solve the problem.
Maybe I'm not explaining my problem clear enough: I have parameters in
my_equation_system.parameters like "boundary_type_01" with a value like "ramp"
and "noundary_amp_01" with a value like 10.0 that is my way to im
On Wed, 15 May 2013, Ataollah Mesgarnejad wrote:
> I want to replace my dirichlet BCs implementation using penalty
> method with the DrichletBoundary class implementation.
>
> My problem is that I want different inhomogenous dirichlet BCs on
> different boundary_ids, but as far as I can see none
Dear all,
I want to replace my dirichlet BCs implementation using penalty method with the
DrichletBoundary class implementation.
My problem is that I want different inhomogenous dirichlet BCs on different
boundary_ids, but as far as I can see none of derived classes from FunctionBase
supports
imple states,independent of their specific history."
Herbert Callen . Thermodynamics and an Introduction to Thermostatics.
- Original Message -
From: "Paul T. Bauman"
To: "Subramanya Gautam Sadasiva"
Cc: "libmesh-users"
Sent: Tuesday, March 19, 2013
On Tue, Mar 19, 2013 at 9:37 PM, Subramanya Gautam Sadasiva <
[email protected]> wrote:
> Hi,
> I have a small question about applying boundary conditions using the
> Dirichlet Boundary Conditions while solving a problem with
> lagrangian_vector elements. How do I constrain only one of the compo
Hi,
I have a small question about applying boundary conditions using the Dirichlet
Boundary Conditions while solving a problem with lagrangian_vector elements.
How do I constrain only one of the components?
Subramanya Sadasiva
"But memory eventually fades. Turbulences damp out, internal stra
34 matches
Mail list logo
