Questions tagged [boundary-conditions]

For questions regarding the choice and/or appropriateness of conditions necessary to model a particular phenomenon with a partial differential equations.

Filter by
Sorted by
Tagged with
1
vote
0answers
52 views

Decrease in slope during convergence analysis

I am using the method of manufactured solutions to perform the order of accuracy testing. I am using a cube for the testing. The cube is size 1m on all sides. I used 5 refinements: $dx = dy = dz = ...
1
vote
1answer
569 views

Well-posedness of linear elasticity boundary conditions

I have several questions regarding suitable boundary conditions for linear elasticity. I have read that in order for modeling linear elasticity to be well-posed, the entire boundary cannot be ...
1
vote
0answers
122 views

Method of manufactured solutions - choice of type of boundary conditions

I am attempting to use the method of manufactured solutions (MMS) for code verification for linear elasticity. However, this is more of a general question regarding the general use of MMS. In MMS, ...
2
votes
1answer
117 views

FEM-Laplace with Dirichlet in only a few points: Nonsingular operator?

Let's consider the FEM discretization of the Laplace operator without boundary conditions, i.e., $$ a(u,v) = \int_\Omega \nabla u \cdot \nabla v - \int_\Gamma (n\cdot \nabla u) v. $$ For one-...
1
vote
1answer
282 views

Are mixed boundary conditions possible in structural mechanics?

For structural mechanics, such as linear elasticity, I am aware of BCs such as a prescribed displacement (Dirichlet) or a prescribed traction (Neumann). Is it possible that a boundary can have a ...
0
votes
1answer
305 views

Proving solution existence and uniqueness of the Helmholtz equation with Robin boundary conditions with complex coefficients

I am trying to solve the Helmholtz equation with Robin boundary conditions with complex coefficients and the weak formulation $$ \iint_\limits\Omega\nabla p_0(x,y)\nabla\left(\overline{v(x,y)}\right)...
2
votes
1answer
697 views

Traction -> stress; stress->displacement gradient

If given a displacement gradient tensor, we can easily obtain the stress tensor (using Hooke's law and the strain-displacement relationship), as well as the traction vector. If given a traction, and ...
2
votes
1answer
236 views

Closed boundary conditions in finite difference method for diffusive-advective equation

I am implementing a finite difference method in solving the diffusive-advective equation: $$ u_t + v \cdot u_x = D\cdot u_{xx} $$ (v, D are constants). Planning to use the operator splitting method (...
1
vote
0answers
214 views

Neumann boundary condition FD implementation for instationnary diffusion equation

I am trying to solve this diffusion equation : $\dfrac{\partial D\dfrac{\partial f}{\partial x}}{\partial x}+S = \dfrac{\partial f}{\partial t}$ ($D$ is not constant and varies according to $x$) with ...
3
votes
2answers
1k views

Implementing no-flux boundary condition reaction-diffusion PDE

I'm having trouble figuring out how to implement boundary conditions for this problem: \begin{align} \frac{\partial n}{\partial t} &= D_n\nabla^2n - \nabla\cdot\left(\frac{\chi}{1+\alpha c}n\nabla ...
2
votes
0answers
87 views

Finite differenced eigenvalue prob. of inhomogeneous boundary conditions?

I am basically asking about eigenvalue problems of differential equations using some finite difference method (FDM). Usually the system is subject to some boundary conditions (BC), e.g., Dirichlet or ...
3
votes
1answer
2k views

Applying Neumann boundaries to Crank-Nicolson solution in python

Consider the heat equation $$u_t = \kappa u_{xx}$$ with boundary conditions of $$u(x,0)=0\\ u(0,t)=100\\ u(l,t)=0$$ Numerical analysis by pyton can be done with ...
2
votes
0answers
39 views

Can I expect reasonable results when expanding in scaled Legendre polynomials?

Imagine I want to compute the eigenvalues of an operator $\hat O$ defined on $ L^2(\mathbb{R})$, however using a properly scaled N-dimensional polynomial basis of $ L^2([-a,a]) $ which fulfills ...
1
vote
0answers
70 views

Radially symmetric system of PDEs in deal.II

I am trying to solve the radially symmetric polar form of the PDE with homogeneous Neumann BC in deal.II on a unit circle: $$ u_t = \Delta u - \nabla \cdot (u \nabla h) $$ $$ h_t = \Delta h $$ I am ...
3
votes
0answers
157 views

Non-reflecting boundary conditions for compressible Navier-Stokes equations

I have some questions about the implementation of non-reflecting OUTFLOW boundary condition for Navier Stokes equations. Following Poinsot, Lele "Boundary Conditions for Direct Simulations of ...
0
votes
2answers
125 views

Second derivative in coordinate invariant form

To solve stationary, incompressible, inviscid and irrotational flow around a circular cylinder, I am using general coordinates. Since the flow is symmetrical, we only consider the upper half of the ...
2
votes
2answers
186 views

Find classical solution of transport equation with FDM

We know the classical solution of transport equation is determined by one initial (boundary?) condition, for example, the solution of $$\frac{\partial u(t,x)}{\partial t}+\frac{\partial u(t,x)}{\...
0
votes
2answers
74 views

Finite difference for 2nd order ode $y'^2+y y''+\frac{2}{x} y y' -0.1 y^2=0$ with $y'(1)=0$ and $y(1)=1$

How to solve second order non-linear ODE $$y'^2+y y''+\frac{2}{x} y y' -0.1 y^2=0$$ subject to $y'(1)=0$ and $y(1)=1$ over the interval $0 < x \le 1$. I turned the equation to a PDE $y'^2+y y''+\...
3
votes
1answer
446 views

Periodic boundary condition in solid

I want to solve a small deformation solid structure problem applying periodic boundary conditions in FEM. The geometry is a square and the equations are: $$ \text{div} \, \sigma = 0 \\ \sigma = f(\...
1
vote
1answer
785 views

Proper boundary conditions for potential flow around cylinder

I am computing the stationary, incompressible, inviscid and irrotational flow around a circular cylinder using a discretization in general coordinates. I derived a PDE and proper boundary conditions ...
6
votes
2answers
223 views

Why naively chopped finite difference matrix works for different ODE boundary conditions

We know finite difference method (FDM) can replace $y''(x)$ as $\frac{1}{h^2}[y(x+h)+y(x-h)-2y(x)]$ or so. One naive way to write down the matrix of the differential operator is like the following, ...
3
votes
1answer
722 views

Radiation boundary condition (heat transfer)

I am looking for reference on how to implement nonlinear boundary conditions. Specifically, I am interested in implementing a radiation boundary condition for heat transfer with the FEM: $-k \frac {\...
4
votes
0answers
174 views

How to apply an integrated constrain condition in FEM?

I'm running some simulation using FEM. In my model I need to apply a constraint condition to the governing equation. My governing equation similar to the diffusion equation as below: $$\frac{\...
1
vote
1answer
2k views

Structural mechanics - traction free = Zero displacement gradient?

This is a follow up question to my question yesterday Structural mechanics traction boundary condition question Does a traction free boundary condition also mean that the displacement gradient is 0 ...
1
vote
0answers
388 views

Structural mechanics traction boundary condition question

In structural mechanics, are the boundary conditions "free surface," "Traction free", "stress free" all equivalent Neumann boundary conditions?
6
votes
0answers
614 views

Understanding Boundary Condition in FEM

I am trying to understand Dirichlet and Neumann boundary conditions in FEM and I wanted to know if my inference is correct. To articulate my understanding, lets consider a simple case of TE and TM ...
0
votes
2answers
132 views

How do you apply boundary conditions in a time-stepping problem?

It looks to me like a very common problem, yet I haven't been able to find any practical guide on the subject despite many hours searching. Here is a clearer statement of my question: I have a ...
1
vote
1answer
146 views

Parallel calculation in finite elements

I am trying to solve a 1 Dimensional eigenvalue of poisson problem: $$\nabla \phi ^2 +\nabla \phi = k\phi$$ with the boundary condition: $\phi (0)=0 , \nabla \phi(1) = 0 $. I could solve this ...
3
votes
1answer
257 views

FEM current toy problem

I am solving the Dirichlet problem $$ \begin{cases} \Delta u = 0, \\ u|_{\partial D} = f, \end{cases} $$ in a $2d$ domain $D$ using the finite element method. What I want to get is the ...
4
votes
1answer
1k views

How to apply non zero Dirichlet boundary condition in finite elements?

I am writing a code for steady state heat transfer on a rectangular domain. I am specifying temperature on the edges - nonzero Dirichlet boundary condition. The equations can be written in form of $$...
2
votes
2answers
619 views

Neumann boundary conditions diffusion equations methods of lines

I want to solve the diffusion equation using the method of lines with Neumann boundary conditions $$ \frac{\partial p}{\partial t}=\frac{\partial^2p}{\partial x^2}\\ \frac{\partial p}{\partial x}(x=0)=...
1
vote
1answer
271 views

Raviart Thomas Mixed Finite Element with Mixed boundary conditions reference request

This is perhaps a more focused version of this question. Using standard notation, I have a code that works and solves the following PDE using a Raviart Thomas mixed method. $$\begin{align} 0 &= ...
6
votes
0answers
137 views

Do practice and theory differ substantially when implementing Neumann Boundary Conditions using a Mixed Method?

I have implemented a pretty straightforward finite element solver for the following Poisson equation. For the purposes of this question we can assume the source term and the Dirichlet data both ...
1
vote
0answers
233 views

the augmented global stiffness matrix is not positive semi-definite using Lagrange Multipliers method within FEM

The augmented global stiffness matrix is not positive semi-definite when using Lagrange Multipliers method to enforce boundary constraints on a simple square domain of integral form: I am considering ...
3
votes
0answers
502 views

Neumann-Neumann boundary intersection

I am trying to solve a Vertex Centered, Finite Difference, Poisson equation $-\nabla^{2}u=f$ with Dirichlet boundaries on the left and bottom and Neumann boundaries on the top and right using ...
1
vote
0answers
53 views

Appropriately handling boundary conditions in a PDE eigenvalue problem

Suppose I have a nonlinear second order Cauchy PDE $\dfrac{\partial p(x,t)}{\partial t}=N(p(x,t))$, where $N:L^2(\mathbb{R})\rightarrow L^2(\mathbb{R})$, and a known fixed point $u(x)$. Mathematically,...
2
votes
0answers
32 views

Minimum box size in order for an object in the center to not 'see' duplications of itself under periodic boundary conditions?

I am working on developing a Dissipative Particle Dynamics (DPD) model of a colloid in a bulk fluid. I essentially want to ensure that the perturbations of my meshed colloid are completely uncoupled ...
1
vote
0answers
160 views

shallow water equation maccormack method

I am trying to make a code for 1D shallow water equation (nonlinear without source terms) using the MacCormack method for sinusoidal wave propagation. My issue is that the wave fluctuates and does not ...
1
vote
0answers
53 views

Extrapolating to non-fluid cells (for Shallow Water Equations), for a shore/beach?

The water height $h$ and 2d velocity field $(u,w)$ are "extrapolated to non-fluid-cells, i.e., setting $h$ equal to the value in the nearest fluid cell." [Bridson] I'm using finite differences. ...
1
vote
0answers
56 views

Discrepancy in estimating boundary stencil for finite difference method

I am trying to estimate the FD stencil for boundary as mentioned in this paper (section 4.1) using MATLAB. The stencil order (6th) is higher than the one mentioned in paper (4th). $$ f_1' +\alpha f_2'...
0
votes
3answers
2k views

Neumann boundary conditions for the upwind scheme applied to the advection equation (Python)

I'm trying to solve the linear advection equation $$u_{t} = cu_{x}, \\ x \in [x_{0}, x_{e}], \quad t \in (0, T], \quad c \in \mathbb{R} \\ u(x,0) = f(x)$$ Note that for $c > 0$, the solution is ...
1
vote
0answers
107 views

How would you specify mixed boundary conditions for a 2D PDE in the matrix used for finite differences

I have the following PDE in 2D: $U_{x} + U_{xx} + U_y + U_{yy} + U_{xy} = f$ where $f$ is a constant. And I'm trying to create a matrix $A$ to solve the PDE through finite differences: $AU = f$. I ...
2
votes
1answer
1k views

Line integral along the edge of an isoparametrically mapped triangle

I need to integrate the following function on the line segment from $P_{1} = \begin{bmatrix} -2\\-1 \end{bmatrix}$ to $P_{2} = \begin{bmatrix} 1\\2 \end{bmatrix}$: $$\int_{P_{1}}^{P_{2}} 4x + y \ ds$$...
1
vote
0answers
39 views

Simulating periodic boundary condition with long_range interaction

Consider 2 particles in a 2D plane. There is a long-range interaction between two particles: $$f (r_i-r_j)= \frac{\vec{r_i} - \vec{r_j}}{|r_i-r_j|^3}$$ $$\vec{r_i} =(x_i,y_i)$$ $r_i$ determines the ...
1
vote
0answers
75 views

Numerically solving generalized eigenproblem with Neumann conditions

I am interested in finding the eigenvalues/eigenfunctions of problems such as $$ \partial_{xx} u = \lambda \partial_{yy} u, $$ which can be solved as the generalised eigenvalue problem $$ \mathbf{A}...
2
votes
1answer
153 views

How to determine if Kelvin-Voigt elements are dissipating stress correctly?

I am using a program that creates viscous/absorbing boundaries by implementing Kelvin-Voigt elements. The theory behind 1-D Kelvin Voigt elements is given in this Wikipedia page. In my case I am ...
6
votes
3answers
402 views

How to impose boundary condition with mixed derivatives?

I have the biharmonic equation on a 2D rectangular domain $\Omega$ with the following boundary conditions: $\Delta^2 u = f$ on $\Omega$ $\nabla u \bullet \mathbf{n}=0$ on $\partial \Omega$ (1) $u_{...
8
votes
0answers
615 views

Numerical implementation of the Dirichlet-to-Neumann map

I am solving the Dirichlet problem $$ \begin{cases} \Delta u = 0, \\ u|_{\partial D} = f, \end{cases} $$ in a $2d$ domain $D$ using the finite element method. What I want to get is the ...
1
vote
0answers
561 views

Implementing Neumann boundary condition for elasticity problem using the finite element method

I am solving the following mixed boundary elasticity problem with the finite element method on the unit circle, using a triangle mesh with 3 nodes on each element. In my problem I know $u_1$ and $t_2$...
5
votes
2answers
482 views

Solving Poisson equation with current BC using FEM

I have a 3-dimensional domain D, with 3 types of BC which I am trying to solve the Poisson equation on. $$-\nabla \cdot (\sigma(x,y,z)\nabla u)=0$$ Insulator (Neumann BC) Electrode set at some ...

1 2
3
4 5
8