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

learn more… | top users | synonyms

0
votes
0answers
6 views

Treatment of Neumann (Traction) boundary conditions using projection methods

I am looking to solve the incompressible Navier-Stokes equations in 3D, using an inflow boundary condition specifying a velocity: $\mathbf{u} = \mathbf{g}_0 \,\, \forall \,\, \mathbf{x} \in \Gamma_u$ ...
3
votes
1answer
42 views

Boundary conditions Chebyshev differentiation

I was wondering if anyone has any experience dealing with boundaries when implementing chebyshev differentiation. I am currently trying to implement a no slip boundary condition to solve the ...
0
votes
1answer
55 views

jump conditions for Poisson/Darcy equation in primal form versus mixed form

Consider the Darcy equation, $$\mathbf{v} + \dfrac{k}{\mu_0}\nabla p = \mathbf{f} \\ \mathrm{div}\; \mathbf{v} = 0$$ If the coefficient $k$ is piecewise constant across an interface $\Gamma$ in the ...
2
votes
1answer
76 views

Jump condition for elliptic equation in standard finite element method

Consider the elliptic PDE $ -\mathrm{div}(k(x)\nabla u) = f(x) \in \Omega$ and $u = 0$ on $\partial \Omega$ If $k(x)$ is piecewise constant across an interface $\Gamma$ in the domain, then we ...
5
votes
1answer
82 views

Is there a jump condition for this PDE? ( Brinkman model , piecewise constant permeability)

The Brinkman equations for steady flow of an incompressible fluid through rigid porous solid are: $-\dfrac{\mu_0}{k}\mathbf{v} - \mathrm{grad}p + \mu_0 \Delta \mathbf{v} =0$ and ...
4
votes
1answer
71 views

Applying Dirichlet b.c. to the Eigenvalue-Problem

If you use a FEM (on the variational formulation), you can discretize some continuous eigenvalue problem, $$L u = \lambda u \ \ \text{on} \ \Omega,$$ into some discrete, generalized eigenvalue ...
0
votes
0answers
111 views

Finite differences and Neumann boundary conditions

I am dealing with a highly nonlinear system of two PDEs. I already have a code to solve the system in case of Dirichlet boundary conditions. The explicit system is: $$ \begin{eqnarray*} \partial_{t}u ...
4
votes
3answers
150 views

Discrete Poisson Equation with Pure Neumann Boundary Conditions

I'm trying to implement the Helmholtz-Hodge Decomposition in 2D, which states that a vector field is composed by a rotational free component, a divergence free component and a harmonic component. ...
5
votes
2answers
112 views

Advice on numerical solution for 2D hyperbolic PDE with zero flux boundary conditions

I would like to numerically solve a hyperbolic PDE of the form $\frac{\partial\theta_t}{\partial t}(x,y)+\frac{\partial\left[\theta_t \gamma_t^x\right]}{\partial x}(x,y)+\frac{\partial\left[\theta_t ...
2
votes
1answer
69 views

Method of lines for inhomogeneous Dirichlet conditions

I understand how to set up the boundary conditions for a steady state problem discretized by Galerkin method, for a time dependent PDE below, $$\frac{\partial}{\partial t} u = c\nabla^2 u + a\nabla u ...
4
votes
0answers
54 views

Elliptical problem with Robin BC

Working in finite differences, I am using a transformation on the temperature variable $\Theta = \int_{T0}^T \kappa(T)dT$ to linearize the steady-state heat equation into a Poisson equation ...
2
votes
2answers
102 views

Solving Stokes flow with walls using Oseen tensor

Introduction I've developed a code to solve for generalised, incompressible 2D Stokes flow $\eta \nabla^2 \mathbf{v} - \nabla p + \mathbf{S} = 0$ $\nabla . \mathbf{v} = 0$ where $\mathbf{S}$ can ...
0
votes
0answers
19 views

Finding a parametrized surface with normal constraints on the boundary

I have a disk-like domain $D$ in $\bf{R}^2$ (let's assume the unit disk) and would like to find a parametric surface $f: D \rightarrow \bf{R}^3$ that satisfies the following constraints on the ...
0
votes
1answer
55 views

GMRES: Making the matrix square without solving for boundaries

How do we define the matrix for GMRES, if we do not want to solve the boundary elements but only the interior ones. I am using pentagonal elements so in a row there are 6 elements (cell itself + 5 ...
1
vote
0answers
80 views

boundary conditions with non-constant coefficients in cell centred finite volume method

Suppose am solving the heat conduction equation in 1d with Dirichlet boundary conditions. The thermal conductivity $k$ is a non constant function. So $-(k(x)u'(x))' = f(x)$ The value of $k$ enters ...
2
votes
1answer
101 views

Neumann / natural BCs in FEA

I'm trying to work out an as-general-as-possible 2D Laplace example for Finite Element Analysis. Starting from $\Delta u = 0$ for an unknown $u(x,y)$, I multiply both sides (well, in practice only the ...
1
vote
1answer
141 views

Moving airfoil boundary conditions

I am trying to simulate a moving airfoil with constant speed (Mach=0.755, aoa=1.25). I solve Euler equations with Roe's method. I have two boundary conditions: Farfield Slip wall (airfoil) For all ...
2
votes
1answer
61 views

Imposing symmetry plane boundary condtition

I want to impose symmetry plane boundary condition for a solid mechanics problem. I googled around and found out that in many places people say to "forbid displacemnts out of symmentry plane and ...
2
votes
0answers
37 views

Interface Formulation at Finite Volume Boundaries when using the Dual Mesh

When using the dual mesh (vertex-centered) for finite volume methods, you end up with a cell center at the boundaries between materials. It is possible that the equations being solved in each ...
1
vote
1answer
176 views

Euler's equations 1d for pipe, Inlet boundary conditions

$\def\rmin{{\mathrm{in}}}$ $\def\l{\left}\def\r{\right}$ $\def\tagl#1{\tag{#1}\label{#1}}$ I am using the one-dimensional finite volume method to calculate the air flow in some tube. For subsonic ...
3
votes
0answers
41 views

Discretizing boundary conditions for vortex methods

I am working on a fluid simulation using vortex methods. For this I must compute the vortex sheet on my boundaries given as: $$ \gamma(\mathbf{x}) - \frac{1}{\pi}\int_S ...
1
vote
0answers
17 views

What type of boundary condition(s) do I need to define for a process with diffusion of tracer concentration?

I want to simulate the diffusion of initial dye concentration in a small estuary with a 2D depth averaged model. Would a radiation boundary condition for concentration of tracers alone be enough or do ...
0
votes
0answers
42 views

Boundary layer in a stream function

i've to solve numerically this equation: where: psi_y=u psi_x=-v that is the boundary layer equation in a stream function. i want to solve this equation with a direct method. which are the ...
4
votes
1answer
96 views

How to handle inflow and outflow boundaries for a non-linear convection-diffusion equation (DGFEM)

Following "A conservative DGM for Convection-Diffusion and Navier-Stokes Problems" (Oden and Baumann), if we have a linear convection-diffusion equation of the following form: $$ ...
4
votes
2answers
94 views

Solving system of differential equations with interconnected boundary conditions

I am trying to solve the following system of differential equations numerically over the domain $x=0$ to $x=D$. The main difficulty is that the boundary conditions are interconnected and depend on the ...
3
votes
0answers
121 views

Transport Equation in a Tube: Source Term on Boundary

I'm modeling mass transport in a flow reactor. The flow reactor is a tube, which allows me to use cylindrical symmetry in solving the Convection-Diffusion-Reaction (CDR) Equation, which governs the ...
3
votes
1answer
162 views

Problem in Discretizing Convection-Diffusion-Reaction equation

I'm trying to solve the Convection-Diffusion-Reaction (CDR) equation on a rectangular domain, using cylindrical coordinates and Finite Difference Methods (FDM) (this approximates a flow reactor). ...
2
votes
1answer
202 views

Implementation of gradient zero boundary conditon in advection-diffusion equation

My question is about Finite Element Method. I want to know how to implement "gradient zero" conditions to advection-diffusion equations in conservative form like, $\frac{\partial \rho}{\partial t} + ...
2
votes
1answer
80 views

Boundary conditions for second order PDE

For a second order PDE, for example heat conduction equation $\frac{\partial T}{\partial t} = \frac{\alpha}{C_p} \nabla^2 T$, is it possible to determine the steady-state (or even transient) solution ...
1
vote
0answers
41 views

how to find outward normal for robin codition [duplicate]

I wrote the code to fix this problem; now I should validate it using a test function to see the error in my code. I can't figure out how starting from u can derive the function g on the edge, in ...
-1
votes
1answer
162 views

Pde problem with robin boundary condition

I have my pde 2D problem with robin condition (form: du/dn +ku=g) to solve with matlab. i have the exact function u and I want to find the function g in robin condition. How can i do it? thanks for ...
2
votes
2answers
83 views

How to impose a constant constraint PDE

What is the best way to impose a "constant constraint" for a PDE? Specifically, I want to solve an eigenvalue problem $Au=\lambda u$ on the rectangle $(0,2\pi)\times(-\pi/2,\pi/2)$ with periodicity ...
3
votes
1answer
153 views

The Lax-Milgram Lemma in FEM with non-homogenous Dirichlet BC

How can show that the prerequisites for the Lax-Milgram Lemma holds if I have different test and trial spaces (which I think is the natural thing to have if at least part of the boundary is ...
4
votes
1answer
122 views

Multigrid stops converging when more grid levels are used

I'm having a problem with multigrid code I wrote. If I solve Laplace's equation in 2D and use more than 5 grid levels, the V-cycles stop converging after a few cycles (see below, convergence factor > ...
9
votes
3answers
372 views

How to deal with curved boundary condition when using finite difference method?

I'm trying to learn about numerically solving PDE by myself. I've been beginning with finite difference method(FDM) for some time because I heard that FDM is the fundament of numerous numerical ...
2
votes
0answers
55 views

Finite Difference for Hamilton Jacobi Belman

I have hjb equation where $V=V(x,t)$ and $u=u(x,t)$ $V_t + \sup(u) [A(x,u)V_x + B(x,u)V_{xx}]=0$ for $x$ in $[0,1]$ and $t$ in $[0,1]$ I have been able to successfuly resolve it numerically having ...
2
votes
1answer
195 views

How to implement boundary condition in this case?

Let's start off by NUMERICALLY solving a 1-D steady-state heat transport problem using IMPLICIT FDM. $DT_{xx}=0; ~T(x=0)=T_{BL}; ~T(x=l)=T_{BR}$ where $D$ is diffusion; $T$ is temperature; subscript ...
8
votes
2answers
451 views

Solid mechanics with finite differences: How to handle “corner nodes”?

I have a question concerning coding boundary conditions for solid mechanics (linear elasticity). In the special case I have to use finite differences (3D). I am very new to this topic, so perhaps some ...
3
votes
1answer
102 views

boundary oscillations with Robin boundary conditions

When solving Poisson's equation on the unit square $\Omega$ with homogeneous Dirichlet boundary conditions for $x=0$ and Robin-type conditions at the rest of the boundary, $$ \begin{cases} -\Delta u = ...
6
votes
3answers
546 views

Applying Dirichlet boundary conditions to the Poisson equation with finite volume method

I would like to know how Dirichlet conditions are normally applied when using the finite volume method on a cell-centered non-uniform grid, My current implementation simply imposes the boundary ...
4
votes
1answer
298 views

boundary conditions for electrostatic problem

I'm solving an electrostatic problem governed by Laplace equation $$-\nabla \cdot (\rho^{-1} \nabla u) = 0$$ in the following domain: a brick ($\Omega_1$) with a cylindrical inclusion ($\Omega_2$), ...
5
votes
1answer
130 views

What Linear Equation Solver should be used for a problem with many dirichlet conditions?

I am solving a laplace equation on a finite-element mesh (tetrahedral, triagonal) and have many say 99% dirichlet conditions compared to the number of unknowns. Is there an efficient way to solve this ...
0
votes
0answers
315 views

solving a set of coupled differential equations in matlab

plot h vs t for the below coupled equations: $a_1\frac{d^2p}{dx^2}-b_1\frac{d^2g}{dx^2} = c\frac{dh}{dt}$ $a_2\frac{d^2p}{dx^2}-b_2\frac{d^2g}{dx^2} = c\frac{dh}{dt}+d\cdot g$ where, a1,b1,a2,b2 ...
5
votes
2answers
315 views

Finite-volume method: can Dirichlet boundary conditions be applied to the integral form?

I would like to apply Dirichlet conditions to the advection-diffusion equation using the finite-volume method. This answer, "How should boundary conditions be applied when using finite-volume method?" ...
3
votes
1answer
726 views

Finite Difference Method Neumann Boundary Condition with Variable Coefficients

Disclaimer In the process of typing up this question, I determine its solution. Since I went through the trouble of typing up the question in its entirety, I will post its answer as well. It may ...
13
votes
1answer
2k views

How should boundary conditions be applied when using finite-volume method?

Following from my previous question I am trying to apply boundary conditions to this non-uniform finite volume mesh, I would like to apply a Robin type boundary condition to the l.h.s. of the ...
5
votes
1answer
213 views

Are we free to choose the position of ghost cells on a non-uniform finite-volume mesh?

Following Hundsdorfer approach the finite volume discretisation of the advection-diffusion equation (conservative form) on non-uniform cell centered grid can be written as, $$ w_j^{\prime} = ...
3
votes
1answer
201 views

Implementation of a contraction force in Fenics

This might be a bit of an odd question, and sorry about the poor wording im having trouble conceptualizing what I want in my head let alone get it out in a way others can understand. Is there any way ...
0
votes
0answers
28 views

Computation of potential flow using dirichlet conditions [duplicate]

I know I already posted this question and I thank you for your answers, but unfortunately I didn't find what I was looking for among them. Anyway now I'll rewrite the question more clearly, because ...
5
votes
1answer
194 views

Potential flow around a non-symmetric obstacle using stream functions

I've seen that there is a way to use the finite differences method, on a cartesian orthogonal grid, to perform calculations on potential flow about an obstacle without using the Neumann conditions, ...