Re: [deal.II] Applying non-homogeneous Neumann BC in step-18
Dear Wolfgang
Sorry for late response. Due to internet censorship in Iran, I didn't have
access to google group for few days.
I consider a specific region on triangulation for Neumann BC (a force
vector in 3D space).
There must be some mistakes in my program, but I couldn't figure out.
I tried to follow your instruction on:
https://groups.google.com/g/dealii/c/_4eeywhhxk4/m/myLiTWCqBgAJ .
So, I used this:
for (unsigned int i = 0; i < dofs_per_cell; ++i)
{
for (unsigned int q_point = 0; q_point < n_face_q_points ; ++q_point)
{
const Tensor<1,dim> n = fe_face_values.normal_vector(q_point);
cell_rhs (i)+=( (n*traction_component*fe_face_values.shape_value(i,
q_point))*
fe_face_values.JxW(q_point));
}//for q_point
}//for dof_per_cell
But in computation of cell rhs i get this error:
/home/amir/eclipse-workspace/mech1/mech.cc:983:22: error: no match for
‘operator+=’ (operand types are ‘double’ and ‘dealii::Tensor<1, 3>’)
983 | cell_rhs (i)+=(
(n*traction_component*fe_face_values.shape_value(i, q_point))*
|
^~
984 | fe_face_values.JxW(q_point));
As error says, the normal vector with 3 component cannot be multiplied by
double values in other terms.
Additionally, my traction component (for example in x-direction) is
preserved for dofs of other directions.
I know these mistakes, but I cannot figure out proper solution for them.
I have read steps 7, 18, 20, 21 and also reviewed lectures 21.~21.65. But I
couldn't find similar codes.
Thanks for your consideration.
Best regards
Amir
On Thursday, July 13, 2023 at 4:43:37 AM UTC+3:30 Wolfgang Bangerth wrote:
> On 7/11/23 05:32, Mohammad Amir Kiani Fordoei wrote:
> > Screenshot from 2023-07-11 14-01-30.png
> > As I understood, I just need to specify a boundary_id on triangulation
> and
> > apply traction force on that with code like this:
> > { - loop over cells
> > -loop over faces per cell && controlling being at_boundary
> > -loop over quadrature points of face
> > -loop over dof per cell && controlling intended boundary_id_NBC
> > -finally adding this term to rhs as
> > cell_rhs (i)+=
> > (fe_face_values.shape_value (i,q_point)
> >
> > *traction_component
> >
> > *fe_face_values.JxW(q_point));}
> > But, I didn't see noticble variation of displacement or norm of stress
> in my
> > solution, even for large traction force.
> > is it necessary to multiply traction_component to "* present_timestep *
> > velocity " (like Dirichlet BC displacement)? (I tried this too, but
> similar to
> > previous one no significant effect of load was observed.)
>
> Is this a boundary where you are still *also* applying a prescribed
> displacement? You can only do one or the other.
>
> Best
> W.
>
> --
>
> Wolfgang Bangerth email: bang...@colostate.edu
> www: http://www.math.colostate.edu/~bangerth/
>
>
>
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/7f105f4d-a4eb-473d-8396-9a4190c2ab14n%40googlegroups.com.
[deal.II] Applying non-homogeneous Neumann BC in step-18
Dear Deal.II group
I am trying to add boundary force (b) to step-18 as below
[image: Screenshot from 2023-07-11 14-01-30.png]
As I understood, I just need to specify a boundary_id on triangulation and
apply traction force on that with code like this:
{ - loop over cells
-loop over faces per cell && controlling being at_boundary
-loop over quadrature points of face
-loop over dof per cell && controlling intended boundary_id_NBC
-finally adding this term to rhs as
cell_rhs (i)+=
(fe_face_values.shape_value (i,q_point)
*traction_component
*fe_face_values.JxW(q_point));}
But, I didn't see noticble variation of displacement or norm of stress in
my solution, even for large traction force.
is it necessary to multiply traction_component to "* present_timestep *
velocity " (like Dirichlet BC displacement)? (I tried this too, but similar
to previous one no significant effect of load was observed.)
instead of procedure described in step-18, can I use functions "
create_right_hand_side" and "create_boundary_right_hand_side" of
VectorTools for applying body and traction force, respectively or they just
belongs to Laplace equation?
Thanks for considering this!
Amir
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/c5c74fbf-7bf5-47e9-906c-6e7db0d23a76n%40googlegroups.com.
Re: [deal.II] Component mask
Dear Wolfgang Thanks for your answer Finally, I succeed in replicating step-18 example. Everything looks fine. I really appreciate you and your colleagues on providing such rich and accurate FEM library. I will continue my study by adding boundary traction force, plasticity, etc, to reach a fully coupled THM model. Maybe one day I will be honored to be a useful member of deal.II community. Best regards Amir On Friday, July 7, 2023 at 1:05:25 AM UTC+3:30 Wolfgang Bangerth wrote: > On 7/6/23 05:57, Mohammad Amir Kiani Fordoei wrote: > > I think I find the problem. I didn't specify map of boundary_values. > > Then you have already found what the problem is. The compiler tells you in > the > error message that you are trying to call the function with the following > argument list: > > interpolate_boundary_values(dealii::DoFHandler<3, 3>&, > unsigned int&, > Step18::IncrementalBoundaryValues<3>, > dealii::ComponentMask)’ > > That is because no such function exists. The function I suspect you *want* > to > call is > > interpolate_boundary_values( > const DoFHandler & dof, > const types::boundary_id boundary_indicator, > const Function & boundary_function, > std::map &boundary_values, > const ComponentMask &component_mask = ComponentMask()); > > which has the 'boundary_values' argument in position 4. > > As a general rule for how you could have spotted the mistake: You are > calling > a non-member-function with only *input* arguments. Since it's not a member > function, it cannot modify the object it is a part of (it is not part of > any > object). Because it only has input arguments, that would mean that it > isn't > actually computing anything -- because none of its arguments is an output > argument to put anything in. So something is conceptually not right with > the call. > > > > I still worried about correctness of my method. > > Well, have you tried what happens? > > Best > W. > > -- > > Wolfgang Bangerth email: bang...@colostate.edu > www: http://www.math.colostate.edu/~bangerth/ > > > -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/f2e022d1-c91f-4d31-ba2f-6e01dd6c31c1n%40googlegroups.com.
Re: [deal.II] Component mask
I think I find the problem. I didn't specify map of boundary_values.
I still worried about correctness of my method.
On Thursday, July 6, 2023 at 2:15:38 PM UTC+3:30 Mohammad Amir Kiani
Fordoei wrote:
> Dear Wolfgang
>
> *The error is:*
>
> /home/amir/eclipse-workspace/mech1/mech.cc:894:45: error: no matching
> function for call to ‘interpolate_boundary_values(dealii::DoFHandler<3,
> 3>&, unsigned int&, Step18::IncrementalBoundaryValues<3>,
> dealii::ComponentMask)’
>
> 894 | VectorTools::interpolate_boundary_values(
>
> | ^
>
> 895 | dof_handler
>
> | ~~~
>
> 896 | , boundary_id_tag
>
> | ~
>
> 897 | , IncrementalBoundaryValues (present_time,
> present_timestep,Vector_of_BC_Value)
>
> |
>
>
> 898 | , fe.component_mask(Vector_of_BC_status) );
>
>
> *That part of code:*
>
>
> prm.enter_subsection("Geometry_info");
>
> unsigned int Number_of_BC=prm.get_integer("Number_of_boundary_id");
>
> prm.leave_subsection();
>
>
> prm.enter_subsection("Boundary_condition_info");
>
>
> std::string temp_BC_value_info
> =prm.get("Assign_boundary_value_info");//getting
> BC value
>
> std::vector vetor_of_BC_value_info= Utilities::
> split_string_list (temp_BC_value_info, ' ');//splitting BC values to
> vector
>
>
> std::string temp_BC_status_info =prm.get("Assign_status_of_component_info"
> );//getting status of BC value
>
> std::vector BC_comp_status_info= Utilities::split_string_list
> (temp_BC_status_info, ' ');//splitting status of BC value
>
>
> unsigned int Number_of_dof_for_each_BC_value_components =prm.get_integer(
> "Number_of_boundary_dof");
>
> unsigned int
> Length_of_BC_value_and_component=Number_of_dof_for_each_BC_value_components+1;//
>
> length of information for BC (value or component status)
>
> prm.leave_subsection();
>
>
> for (unsigned int i=0;i
> {
>
> unsigned int boundary_id_tag = std::stoi(vetor_of_BC_value_info
> [Length_of_BC_value_and_component*i]);//getting tag of BC
>
> std::vector Vector_of_BC_Value;
>
> std::vector Vector_of_BC_status;
>
>
> for (unsigned int j=1;j dof of BC
>
> {
>
> double boundary_id_value= std::stod(vetor_of_BC_value_info
> [Length_of_BC_value_and_component*i+j]);
>
> Vector_of_BC_Value.push_back(boundary_id_value);//assigning value of each
> component
>
> bool boundary_id_comp_status=false;
>
> if (BC_comp_status_info [Length_of_BC_value_and_component*i+j]=="true")
>
> {
>
> boundary_id_comp_status=true;
>
> }
>
> else
>
> {
>
> boundary_id_comp_status=false;
>
> }
>
> Vector_of_BC_status.push_back(boundary_id_comp_status);
>
> }
>
>
> VectorTools::interpolate_boundary_values(
>
> dof_handler
>
> , boundary_id_tag
>
> , IncrementalBoundaryValues (present_time, present_timestep
> ,Vector_of_BC_Value)
>
> , fe.component_mask(Vector_of_BC_status) );
>
> //, Vector_of_BC_status );
>
>
> temp_BC_status_info.clear();
>
> temp_BC_value_info.clear();
>
> Vector_of_BC_status.clear();
>
> Vector_of_BC_Value.clear();
>
> }
>
>
> I also added my doe and its prm
> I really appreciate your kindness.
> On Wednesday, July 5, 2023 at 7:23:27 PM UTC+3:30 Wolfgang Bangerth wrote:
>
>> On 7/4/23 09:28, Mohammad Amir Kiani Fordoei wrote:
>> > So, I defined two extra vector which indicate presence or
>> absence
>> > (being free dof) of displacement vector as below
>> > a) [true,true,true]
>> > b) [true, false,true]
>> > But, I couldn't use them as input in interpolate_boundary_values
>> (error: no
>> > matching function for call to ‘interpolate_boundary_values)
>> > Also, I couldn't find a matching function FEValuesExtractors::Vector
>> for these
>> > cases.
>>
>> Amir,
>> can you show us the code and what the complete error message is?
>>
>> Best
>> W.
>>
>> --
>>
>> Wolfgang Bangerth email: bang...@colostate.edu
>> www: http://www.math.colostate.edu/~bangerth/
>>
>>
>>
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/24b37c05-d454-48d1-974c-7a2f4cba21c7n%40googlegroups.com.
Re: [deal.II] Component mask
Dear Wolfgang
*The error is:*
/home/amir/eclipse-workspace/mech1/mech.cc:894:45: error: no matching
function for call to ‘interpolate_boundary_values(dealii::DoFHandler<3,
3>&, unsigned int&, Step18::IncrementalBoundaryValues<3>,
dealii::ComponentMask)’
894 | VectorTools::interpolate_boundary_values(
| ^
895 | dof_handler
| ~~~
896 | , boundary_id_tag
| ~
897 | , IncrementalBoundaryValues (present_time,
present_timestep,Vector_of_BC_Value)
|
898 | , fe.component_mask(Vector_of_BC_status) );
*That part of code:*
prm.enter_subsection("Geometry_info");
unsigned int Number_of_BC=prm.get_integer("Number_of_boundary_id");
prm.leave_subsection();
prm.enter_subsection("Boundary_condition_info");
std::string temp_BC_value_info =prm.get("Assign_boundary_value_info");//getting
BC value
std::vector vetor_of_BC_value_info= Utilities::
split_string_list (temp_BC_value_info, ' ');//splitting BC values to vector
std::string temp_BC_status_info
=prm.get("Assign_status_of_component_info");//getting
status of BC value
std::vector BC_comp_status_info= Utilities::split_string_list
(temp_BC_status_info, ' ');//splitting status of BC value
unsigned int Number_of_dof_for_each_BC_value_components =prm.get_integer(
"Number_of_boundary_dof");
unsigned int
Length_of_BC_value_and_component=Number_of_dof_for_each_BC_value_components+1;//
length of information for BC (value or component status)
prm.leave_subsection();
for (unsigned int i=0;i Vector_of_BC_Value;
std::vector Vector_of_BC_status;
for (unsigned int j=1;j(present_time, present_timestep
,Vector_of_BC_Value)
, fe.component_mask(Vector_of_BC_status) );
//, Vector_of_BC_status );
temp_BC_status_info.clear();
temp_BC_value_info.clear();
Vector_of_BC_status.clear();
Vector_of_BC_Value.clear();
}
I also added my doe and its prm
I really appreciate your kindness.
On Wednesday, July 5, 2023 at 7:23:27 PM UTC+3:30 Wolfgang Bangerth wrote:
> On 7/4/23 09:28, Mohammad Amir Kiani Fordoei wrote:
> > So, I defined two extra vector which indicate presence or absence
> > (being free dof) of displacement vector as below
> > a) [true,true,true]
> > b) [true, false,true]
> > But, I couldn't use them as input in interpolate_boundary_values (error:
> no
> > matching function for call to ‘interpolate_boundary_values)
> > Also, I couldn't find a matching function FEValuesExtractors::Vector for
> these
> > cases.
>
> Amir,
> can you show us the code and what the complete error message is?
>
> Best
> W.
>
> --
>
> Wolfgang Bangerth email: bang...@colostate.edu
> www: http://www.math.colostate.edu/~bangerth/
>
>
>
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/754b4e2e-77fa-4622-8be2-c4b3d4453e7bn%40googlegroups.com.
/* -
*
* Copyright (C) 2000 - 2022 by the deal.II authors
*
* This file is part of the deal.II library.
*
* The deal.II library is free software; you can use it, redistribute
* it, and/or modify it under the terms of the GNU Lesser General
* Public License as published by the Free Software Foundation; either
* version 2.1 of the License, or (at your option) any later version.
* The full text of the license can be found in the file LICENSE.md at
* the top level directory of deal.II.
*
* -
*
* Author: Wolfgang Bangerth, University of Texas at Austin, 2000, 2004, 2005,
* Timo Heister, 2013
*/
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
namespace Step18
{
using namespace dealii;
class ParameterReader : public Subscriptor
{
public:
ParameterReader(ParameterHandler &);
void read_parameters(const std::string &);
private:
void declare_parameters();
ParameterHandle
[deal.II] Component mask
Dear deal.ii developers,
I want to model time-dependent small elastic deformations (step-18).
Inputs (like triangulation, B.C., etc.) are handled by ParameterHandler and
prm file.
In applying multiple Dirichlet B.C. by looping over
"VectorTools::interpolate_boundary_values"
for a 3D problem, I will apply a displacement vector with two
scenarios:
1) displacement {1,0,2} applies on direction 1 to 3 for boundary id=10;
2)displacement {1,0,2} applies only on direction 1 and 3 for boundary id=20;
As mentioned in documentation, I tried to make component mask for each
scenario by hand.
So, I defined two extra vector which indicate presence or absence
(being free dof) of displacement vector as below
a) [true,true,true]
b) [true, false,true]
But, I couldn't use them as input in interpolate_boundary_values (error: no
matching function for call to ‘interpolate_boundary_values)
Also, I couldn't find a matching function FEValuesExtractors::Vector for
these cases.
Thanks in advance
Amir
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/f917b3d2-f5fb-4af2-b5ba-05a07733ac31n%40googlegroups.com.
[deal.II] Applying multiple in-homogeneous Dirichlet boundary condition
Hello
I am trying to extend step-18 to a runtime program.
To this end, I generate a triangulation made of several
generate_from_name_and_arguments and merged them.
then I assigned locus of various boundary id to that merged triangulation
(obviously only at-boundary faces).
Now I would like to apply non-homogeneous Dirichlet boundary values to each
boundary id through looping over boundary ids and using
VectorTools::interpolate_boundary_values.
As follow:
for (unsigned int i=0;i<=Number_of_DBC;i++)
{
VectorTools::interpolate_boundary_values(
dof_handler,
i,
IncrementalBoundaryValues(present_time, present_timestep),
boundary_values);
}
And modifying IncrementalBoundaryValues::vector_value to this:
template
void
IncrementalBoundaryValues::vector_value(const Point & /*p*/,
Vector &values) const
{
AssertDimension(values.size(), dim);
std::vector Displacement_Magnitude_Factor (3,0);
for (int k=0;k<3;k++)
{
double disp_mag_factor;
std::cout<< "Please Enter component "
[deal.II] Applying multiple in-homogeneous Dirichlet boundary condition
Hello
I am trying to extend step-18 to extend it to a runtime program.
To this end, I generate a triangulation made of several
generate_from_name_and_arguments and merged them.
then I assigned locus of various boundary id to that merged triangulation (
obviously only at-boundary faces).
Now I would like to apply non-homogeneous Dirichlet boundary values to each
boundary id through looping over boundary ids and using
VectorTools::interpolate_boundary_values.
As follow:
for (unsigned int i=0;i<=Number_of_DBC;i++)
{
VectorTools::interpolate_boundary_values(
dof_handler,
i,
IncrementalBoundaryValues(present_time, present_timestep),
boundary_values);
}
And modifying IncrementalBoundaryValues::vector_value to this:
template
void
IncrementalBoundaryValues::vector_value(const Point & /*p*/,
Vector &values) const
{
AssertDimension(values.size(), dim);
std::vector Displacement_Magnitude_Factor (3,0);
for (int k=0;k<3;k++)
{
double disp_mag_factor;
std::cout<< "Please Enter component "
Re: [deal.II] Using Dealii to developing aa finite element software on Windows
Dear Timo I appreciate your answer. We tried updating our compilers. But a lot of errors were reported during the compilation of our software. On Mon, Dec 28, 2020, 20:05 Timo Heister wrote: > Hi Amir, > > deal.II 9.1 requires c++11 > deal.II 9.2 requires c++11 > deal.II 9.3 will require c++14 > > Your only choice would be to use an older version of deal.II but I > would strongly suggest updating your compilers instead. > > On Sun, Dec 27, 2020 at 12:09 AM amir kiani > wrote: > > > > Hello. > > Recently I started to learn finite element programming. > > By surfing the net, I came to C++ and Deal.ii. > > One of my concerns in using Deal.ii and C++ is about any possible future > limitations. > > By reading tutorials on the website of Deal.ii, I realized that using > Deal.ii directly on windows is experimentally supported. > > Now my main question: > > As presented on the Deal.ii Website, If I build Deal.ii via each > approach (and in any operating system like Linux), Can I use classes and > libraries of Deal.ii in other C++ programs? Consider I want to use Deal.ii > features in C++03 with Visual Studio 2008 on Windows 7. > > Because I am working on developing an existing software, I am restricted > to use this standard of C++ and this version of Visual Studio. > > I appreciate all answers in advance > > > > -- > > The deal.II project is located at http://www.dealii.org/ > > For mailing list/forum options, see > https://groups.google.com/d/forum/dealii?hl=en > > --- > > You received this message because you are subscribed to the Google > Groups "deal.II User Group" group. > > To unsubscribe from this group and stop receiving emails from it, send > an email to dealii+unsubscr...@googlegroups.com. > > To view this discussion on the web visit > https://groups.google.com/d/msgid/dealii/a4df4788-b1c4-44e6-97c1-b06b6912554dn%40googlegroups.com > . > > > > -- > Timo Heister > http://www.math.clemson.edu/~heister/ > > -- > The deal.II project is located at http://www.dealii.org/ > For mailing list/forum options, see > https://groups.google.com/d/forum/dealii?hl=en > --- > You received this message because you are subscribed to the Google Groups > "deal.II User Group" group. > To unsubscribe from this group and stop receiving emails from it, send an > email to dealii+unsubscr...@googlegroups.com. > To view this discussion on the web visit > https://groups.google.com/d/msgid/dealii/CAMRj59H%2B8vXqn0WJb4aAQ_MyfWkfAz7CogB9kp55pycfe6vfgg%40mail.gmail.com > . > -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/CALNzT0coBMX0fHN36YmFGQeG%2B5HnOX3rYPXkxwW89BN3gOkcpg%40mail.gmail.com.
Re: [deal.II] Using Dealii to developing aa finite element software on Windows
Thanks dear Wolfgang The features of Deall.ii is very good. I will try to use these wonderfull features. You and your colleagues did a great job. I wish best for you. On Mon, Dec 28, 2020, 22:38 Wolfgang Bangerth wrote: > On 12/26/20 10:08 PM, amir kiani wrote: > > By reading tutorials on the website of Deal.ii, I realized that using > > Deal.ii directly on windows is experimentally supported. > > That may no longer be an accurate description. deal.II 9.2 is used in a > number > of industrial applications that run on Windows. It definitely works on > Windows. > > Best > W. > > -- > > Wolfgang Bangerth email: bange...@colostate.edu > www: http://www.math.colostate.edu/~bangerth/ > > -- > The deal.II project is located at http://www.dealii.org/ > For mailing list/forum options, see > https://groups.google.com/d/forum/dealii?hl=en > --- > You received this message because you are subscribed to the Google Groups > "deal.II User Group" group. > To unsubscribe from this group and stop receiving emails from it, send an > email to dealii+unsubscr...@googlegroups.com. > To view this discussion on the web visit > https://groups.google.com/d/msgid/dealii/75e252ef-0e7e-c5dc-6703-57cc514ee84c%40colostate.edu > . > -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/CALNzT0d%2BU8DgVEaWWrN4cXAhV-Viua2VWfeZZUmapAetZ-FqWg%40mail.gmail.com.
[deal.II] Using Dealii to developing aa finite element software on Windows
Hello. Recently I started to learn finite element programming. By surfing the net, I came to C++ and Deal.ii. One of my concerns in using Deal.ii and C++ is about any possible future limitations. By reading tutorials on the website of Deal.ii, I realized that using Deal.ii directly on windows is experimentally supported. Now my main question: As presented on the Deal.ii Website, If I build Deal.ii via each approach (and in any operating system like Linux), Can I use classes and libraries of Deal.ii in other C++ programs? Consider I want to use Deal.ii features in C++03 with Visual Studio 2008 on Windows 7. Because I am working on developing an existing software, I am restricted to use this standard of C++ and this version of Visual Studio. I appreciate all answers in advance -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/a4df4788-b1c4-44e6-97c1-b06b6912554dn%40googlegroups.com.
[deal.II] using Deal.ii in programming a new finite element software
Hello. Recently I started to learn finite element programming. By surfing the net, I came to C++ and Deal.ii. One of my concerns in using Deal.ii and C++ is about any possible future limitations. By reading tutorials on the website of Deal.ii, I realized that using Deal.ii directly on windows is not fully developed yet. Now my main question: As presented on the Deal.ii Website, If I build Deal.ii via each approach (and in any operating system like Linux), Can I use classes and libraries of Deal.ii in other C++ programs? Consider I want to use Deal.ii features in C++ with Visual Studio 2008 on Windows 7. I appreciate all answers in advance -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/e72055db-2f5c-40f3-9977-1a3b608f30a5n%40googlegroups.com.
Re: [deal.II] Thermohydromechanical finite element modeling
Dear Wolfgang Thank you for your support. Good luck and good wishes. Amir On Friday, January 24, 2020 at 10:09:42 PM UTC+3:30, Wolfgang Bangerth wrote: > > On 1/24/20 9:45 AM, amir kiani wrote: > > I am trying to find an FEM library that helps me to model a fully > > coupled THM modeling. > > I want to implement a THM formulation exactly the same as the attached > file. > > Can I use deal.ii for this purpose? > > Yes, but there is no code freely available to the best of my knowledge > that already does what you want to do. But people have implemented the > kind of coupled problem you are looking at many times over. For example, > I have had at least two students in my classes who as a class project > implemented a poroelastic solver -- essentially your hydromechanical > component. Adding a temperature solver is not going to make the problem > conceptually more difficult, it's just a bit of extra work. > > So yes, you can use deal.II! > > Best > W. > > -- > > Wolfgang Bangerth email: bang...@colostate.edu > > www: http://www.math.colostate.edu/~bangerth/ > -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/5a3fe118-07dc-49ad-9618-8d4887668e23%40googlegroups.com.
