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
0
votes
2answers
605 views

Correctly setting boundary condition for periodic linear elasticity problem

From an old, wise engineering book Peterson's Stress Concentration Factors (http://eu.wiley.com/WileyCDA/WileyTitle/productCd-0470048247.html page 324) I've got the following problem: There is 2D ...
1
vote
0answers
381 views

Implement Robin boundary condition (finite volume)

I have a PDE equation with Robin Boundary condition in an annulus system and I should solve it by finite volume method: \begin{align} \frac{\partial T_f}{\partial t} - k \left(\frac{\...
5
votes
3answers
991 views

Computing accurate fluxes with FEM

I have solved Poisson equation on a 3d domain with neumann and dirichlet boundary condition. I get the potential, take the gradient for each element and integrate on a surface of an element, I do this ...
1
vote
0answers
579 views

How to implement outgoing wave boundary condition

I am solving the one dimensional wave equation: $0=\Box\phi = -\partial_t^2\phi + \partial_r^2\phi ,$ using a Crank-Nicolson finite difference scheme, in the domain $r\in[0,R]$. First, I define $\xi\...
7
votes
1answer
510 views

Discrete wave simulation - absorbing boundaries?

I wrote a simple 2D wave simulation using the following equations: $$\frac{\partial^2 u}{\partial t^2}=c^2\nabla^2u$$ Where $\nabla^2$ is the discrete laplace operator using a Von Neumann neighborhood ...
1
vote
0answers
97 views

Discontinuos Galerkin Method - inhomogeneous advection problem

I'm currently trying to get into this topic. I've learned that the basic scheme for the advection problem ($D_{x}u+a*D_{t}u=0$) can be solved in a scheme like $$ M^{k}\frac{d}{dt}u^{k}_{h}-(S^{k})^{T}...
0
votes
1answer
770 views

How to implement an integral condition when solving a BVP in MatLab

I am trying to solve a coupled system of ODE's using MATLAB's bvp4c function. I want to impose the condition that $$\int_{0}^{\pi} y_{1}(t) y_{1}(t) dt = 1,$$ where ...
2
votes
0answers
211 views

Outflow boundary condition - second derivative of velocity

Consider a fluid flow simulation in a pipe. At the outflow, instead of explicitly imposing a boundary condition, I linearly extrapolate information from the interior (for velocity components). This ...
1
vote
0answers
126 views

boundary condition in 2-D planar finite difference problem

I'm working on a 2-D planar finite difference code. My differencing scheme at the boundary nodes involves introducing ghost nodes in the computation. My code also involves a multi-dimensional ...
0
votes
1answer
235 views

Well-posedness of Elasticity Boundary Conditions

For geotechnical engineering problems, it is common to fix a single component of displacement along a boundary as a Dirichlet boundary condition (roller boundary condition). However, I'm having ...
1
vote
4answers
320 views

Mass transport in porous media

My goal is to numerically solve the convection-diffusion equation of the form: $$\frac{\partial C}{\partial t} = \nabla (D \nabla C) - \nabla (v C)$$ $C$ is concentration and $D$ is the diffusivity ...
1
vote
0answers
254 views

Breather solutions of Sine-Gordon Using Finite Differences

I'm attempting to simulate a standing breather of the form $$ u(x,t)=4\tan^{-1}\left(\sqrt{3}\cos\left(\frac{t}{2}\right)sech\left(\frac{\sqrt{3}x}{2}\right)\right)$$ for the Sine-Gordon equation $$u_{...
4
votes
1answer
426 views

Symmetric boundary condition

I am solving a flow physics problem, in which I encounter a symmetric boundary. So I set the boundary conditions to be $$\frac{\partial u}{\partial r} = \frac{\partial v}{\partial r}=\frac{\partial w}...
0
votes
1answer
156 views

Initial Condition in a Numerical Problem

In a initial value problem does the initial condition has to satisfy the boundary condition and the governing equation? For example: If a non-homogeneous Neumann boundary condition for the pressure ...
2
votes
1answer
252 views

How do I simulate an open end?

When simulating a partial differential equation describing a physical phenomenon like vibrations on a string, fluid flow in a chamber or quantum wave functions, the most straight-forward way is to ...
4
votes
0answers
1k views

How to implement boundary conditions on Finite Difference WENO5 scheme for the Euler equations

I'm implementing a Finite Difference WENO5 with Lax-Friedrich flux splitting on a uniform, structured grid to solve the 2D Euler equations of fluid dynamics on a rectangular domain in cartesian ...
1
vote
1answer
834 views

Flux boundary condition in solute transport

I have a pretty naive question, though important to me. Usually when solving the following PDE in solute transport: $\frac{{\partial C}}{\partial t } = \nabla. (D\nabla C -vC )=0,$ one can be asked to ...
0
votes
2answers
565 views

Applying boundary conditions in a simple one element rotated mesh

I'm testing a Finite Elements code. To do this, i have created a simple one element quadrilateral mesh to verify the solution: This is a plane stress elasticity problem, and the local stiffness is ...
3
votes
1answer
368 views

How to reformulate a 1/x^2 singular term to 1/x so that bvp4c can solve it?

I have a ''cosmetic'' problem with a singular term in my Matlab script. I am trying to solve the following system of differential equations: $$ \begin{aligned} y_1' &= y_2,\\ y_2' &= \frac{...
1
vote
1answer
2k views

OpenFoam Mapped Boundary condition [closed]

I am using OpenFoam to do an LES simulation. I am using mapped boundary condition on the inlet plane. I know it maps from a source plane to the target plane(here the inlet plane). In "Pitzdaily-mapped"...
2
votes
0answers
185 views

Mass conservation in atmospheric continuity equation numerical solution

My phd project is heavily related to numerical modeling of planetary atmospheres. In particular now I am dealing with a particular expression of the continuity equation, involving a thermodynamic flux....
0
votes
1answer
255 views

How can an engineering student become a computational scinece expert in a short time [closed]

How can a student with zero computing or programming language knowledge, few engineering mathematics knowledge, understand computational science especially Finite Element Modelling (FEM) from ...
2
votes
1answer
2k views

Solving the Poisson equation with Neumann Boundary Conditions - Finite Difference, BiCGSTAB

I am trying to solve the Poisson equation in a rectangular domain using a finite difference scheme with a rectangular mesh. I have happily generated the matrix system of equations Ax = b which is ...
0
votes
1answer
273 views

The system matrix and the right hand side for diffusion equation with staggered grid

In the following staggered grid setting, I want to solve diffusion equation as a linear system. $$\frac{\partial u}{\partial t} = \frac{\partial^2 u}{\partial x^2} + \frac{\partial^2 u}{\partial y^2}$...
0
votes
1answer
515 views

Outlet boundary condisions in lattice Boltzmann Method

Here is the flow past a square cylinder configuration. The problem is a viscous and incompressible with parabolic velocity profile using freestream velocity U across single cylinder. I use the single ...
1
vote
2answers
1k views

Boundary condition for Pressure in Navier-Stokes equation

I am reading the paper, http://math.mit.edu/~gs/cse/codes/mit18086_navierstokes.pdf I could not really understand the description. Could someone explain a little bit more? It says, "For the lid ...
2
votes
2answers
292 views

$O(h^2)$ convergence for Elliptic PDE

I am trying to solve an elliptic PDE in 2-D: $$-\nabla^{2} u = 20tanh(10x-5)(10-10tanh(10x-5)^2) = f$$ I know that the solution is $u = tanh(10x-5)$ but I am unable to get $O(h^2)$ solution with a ...
3
votes
2answers
1k views

Dirichlet boundary condition for sparse matrix - Solving Ax=b only for free nodes?

I am solving Biot equation with sparse matrix in MATLAB. I have no problem with global sparse matrices assembly, but when I assign Dirichlet boundary condition, it is so slow. From this topic, How to ...
4
votes
2answers
2k views

How to handle boundary conditions in Crank-Nicolson solution of IVP-BVP?

I'm trying to solve the PDE for $c(r,t)$ $$c_t=(1/r)(rJ)_r$$ using Crank-Nicolson, and I'm having difficulty with the boundary conditions. $J$ is the flux, the initial condition is $c(0,r)=c_{init}$, ...
1
vote
0answers
135 views

How to apply Dirichlet boundary condition for lowest order Nedelec element over tetrahedral domain? [closed]

The Dirichlet boundary condition is $n \times \vec{E} = n \times \vec{g}$. I don't know how to calculate the boundary value. Any kind of help would be appreciated.
3
votes
0answers
184 views

Multigrid for Robin boundary conditions

I have a question regarding the treatment of Robin boundary conditions in a multigrid solver. I am solving the Poisson equation in $\Omega=(0,1)^2$ with Robin boundary conditions on the boundary, $$- \...
3
votes
1answer
281 views

Pressure-Pressure BC

My FDM code simulates Backward Facing Step flow when I use conventional BCs such as defining velocity profile at inlet and fully developed condition at outlet. I have validated the results and it ...
2
votes
1answer
244 views

Imposing boundary conditions for PDE quadratic eigenvalue problem

I have a quadratic eigenvalue problem of the form: $$(A_2 s^2 + A_1 s + A_0)\hat{v} = 0$$ where $s$ is the eigenvalue. The matrices $A_i$ contain derivatives up to order six, and I have six boundary ...
1
vote
0answers
197 views

Spherical Advection Discretization (boundary nodes)

Consider the spherical advection problem: describing the conservation of a property $u$ in a closed spherical domain. $$ \frac{\partial u}{\partial t}+\frac{1}{r^2}\frac{\partial }{\partial r}\left(r^...
1
vote
1answer
295 views

FD implementation of Absorbing Boundary condition for acoustic wave

I am simulating acoustic wave equation in which absorbing boundary condition has to be applied. It is applied in two ways Ist method is as mentioned in this paper. Boundary condition at bottom is (...
1
vote
2answers
1k views

Boundary conditions in conjugate gradient method for poisson's equation

I want to use the conjugate gradient method to solve poisson's equation in an electrostatic setup: \begin{align} \rho=-\nabla^2\phi \end{align} I am however a little confused when it comes to the ...
4
votes
0answers
71 views

Differential eigenproblem with eigenvalue in boundary condition

Statement of the problem I need to (numerically) solve an eigenproblem of the type $$-\omega^2\mathcal{D}_1\vec{x}=\mathcal{D}_2\vec{x}$$ on the interval $[-1,1]$, where $\mathcal{D}_1$ and $\...
1
vote
1answer
86 views

Perfectly matched layer simulation with two vibrating sources

I'm doing PML simulation, and while I was looking for correct geometry for my layers of propagation, I got a question. My device has two vibration sources, and this vibration will propagate through ...
1
vote
2answers
112 views

Definition of inflow boundary in CFD

If $w$ is the vector of conservative variables, $f=f(w)$ the flux function, I think have read somewhere (I can not find it anymore) that the inflow boundary $\Sigma$_ is characterized by: $\Sigma_{-}...
5
votes
1answer
792 views

Best practice for dealing with Dirichlet boundary conditions in finite-difference schemes: add artificial unknowns?

I know of at least two ways of dealing with Dirichlet boundary conditions in finite-difference schemes (and, to a lesser extent, finite-element schemes). Here I'm thinking of solving Poisson's ...
2
votes
1answer
239 views

Coupling Boundary Condition of one PDE with source term of another PDE

We have a system of equations, wherein the BC of one PDE is coupled with the source term of another PDE. We have a regular 2D unit grid in x and y. There are two PDEs to be solved The first PDE (...
1
vote
1answer
250 views

imposing “measured data” to Dirichlet boundary conditions in fenics

I'm relatively new to fenics and I just looked through all questions related to Dirichlet boundary conditions. I don't seem to find a well-described question or answer about what I'm about to ask. I'...
2
votes
1answer
1k views

heat equation on bounded and unbounded domain

I have been reading about the heat equation and I am confused about uniqueness in the case when the domain is bounded and when it is not. In the book I am following, it is common to write the heat ...
1
vote
0answers
318 views

How to add a Ricker Wavelet (Mexican Hat) to a 2D/ 3D fem mesh?

I have a 2D square mesh and a 3D beam shaped mesh and I want to propagate a seismic wave in them. I am trying to simulate them using Open source FEM codes (fenics). I have left the top surface to be ...
2
votes
1answer
1k views

Pressure boundary condition in Navier-Stokes equations

I would like to solve 3D transient incompressible Navier-Stokes with FEM, Newton method, Schur-based preconditioner, Lagrangean P2/P1 elements (no stabilization), in a rigid pipe discretized with ...
3
votes
2answers
431 views

2D Laplace problem with mixed boundary conditions using Conjugate Gradients

I am being asked for one of my classes to solve 2D Laplace equations with mixed boundary conditions using the Conjugate Gradient method. The equations and conditions are given as: $$ \frac{\partial^...
1
vote
0answers
589 views

Issues with self-consistent Poisson-Schrodinger solver

I'm currently in the process of writing a self-consistent Schrodinger-Poisson solver for a device heterstructure (High Mobility Electron Transistor). The algorithm is based off of this journal(1). I ...
5
votes
2answers
3k views

Recovering pressure from velocity or streamfunction fields

I am interested in 2D channel flow of an incompressible Stokes fluid (Re << 1), with periodic boundary conditions in the x-direction and no-slip at the walls in the y-direction. I have existing ...
2
votes
1answer
147 views

Definitiones of solid, fluid, and boundary nodes in the context of LBM

In the family of Lattice Boltzmann (LB) methods, like many others, one deals with three types of node; namely, fluid node, solid node, and boundary node. (I know; the boundary node is a subtype of ...
3
votes
2answers
228 views

Determine numerical infinity for Schrodinger equation $−\psi''(x) + x^ 2 \psi(x) = E\psi(x)$

Consider the following Schrodinger equation for the harmonic oscillator with real $x$: $$ −ψ''(x) + x^ 2 ψ(x) = Eψ(x). $$ I solve the last equation using shooting method and implicit Runge-Kutta ...

1 2 3
4
5
8