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
2
votes
1answer
76 views

Non-reflective boundary condition

I'm currently solving incompressible Navier-Stokes system of equations with periodic flow and high viscosity. Is there any outlet boundary types that avoids the reflection of flow from the outlet back ...
4
votes
1answer
102 views

Computing second derivatives with Neumann boundary condition

I am implementing a finite difference method for a PDE with a Neumann boundary condition. I will simplify my question to a single dimension. Suppose I have a PDE $$\frac{\partial u}{\partial t} = \...
9
votes
1answer
428 views

How to solve a second order differential equation (diffusion) with boundary conditions using Python

I am having trouble implementing a model from a publication. Huang, K-L.; Holsen, T.M.; Selman, J.R. Ind. Eng. Chem. Res. 2003, 42, 15, 3620–3625 scihub link: https://sci-hub.se/10.1021/ie030109q I ...
1
vote
1answer
62 views

Finite element method with two different Dirichlet boundary conditions

I have the problem like this $$ -\triangle u = f \ \ on\ \Omega \\ u = g_1 \ \ on \ \partial \Omega_1 \\ u = g_2 \ \ on \ \partial \Omega_2 $$ If we choose $$ V_1 = \{ \nu_1 \in H^1 : \nu_1 = 0 \ ...
0
votes
1answer
78 views

deal.ii - ParaView "warp by scalar" of my output is not continuous

During our finite element course, we've solved the linear elasticity problem in 2D on a square (GridGenerator::hyper_cube) with $Q_1$ bilinear finite elements in ...
0
votes
1answer
77 views

Displacement field not correct?

Consider the elastic equation $$- \operatorname{div}(C \nabla \mathbf{u}) = \mathbf{f}$$ as presented in step-8. Here $\mathbf{u}$ is the displacement vector, let's consider the 2d case. As you can ...
4
votes
1answer
133 views

Weak Formulation of Poisson's Equation for Electrostatics with Surface Charge

I am currently attempting to use FEnICS to solve an electrostatic problem with two materials of different permittivity $\varepsilon_1$ and $\varepsilon_2$ forming an interface: Consider a domain $\...
1
vote
0answers
55 views

Constrained optimization for non-linear equations in octaveGNU

I have installed Optim1.6.1 package. I would like to solve a system of equations in non linear finite element analysis using constraints as u=1 at certain nodes. u=0 at certain nodes. Typically I find ...
1
vote
0answers
119 views

Weak form of Elliptic problem with mixed Dirichlet & Neumann conditions

Let $\Omega \subset \mathbb{R}$ in a bounded polygon domain and $f:\Omega \to \mathbb{R}$ known function.We split the boundary into two parts $\partial \Omega_{1}$ and $\partial \Omega_{2}$ such that $...
0
votes
0answers
23 views

Reynolds boundary conditions

I came across this paper comparing various boundary conditions. I am particularly interested to understand how to obtain the Reynolds boundary conditions (refer to equation 28). $$\left( \frac{1}{c} \...
0
votes
0answers
35 views

Insert a boundary condition without removing periodicity assumption

I've to perform a multi-fluid internal flow simulation with a code which intrinsically assumes periodicity on a given direction (say z) on a cartesian grid. Nonetheless, the problem I'm trying to ...
3
votes
2answers
217 views

How to solve second order coupled non linear differential equations

For a project I am doing, I have to solve the following system of differential equations numerically using my own code: $$ x^2K'' = KH^2 + K(K^2-1) $$ and, $$ x^2H'' = 2K^2H + \alpha H(H^2-x^2) $$ ...
1
vote
0answers
47 views

Why is this Advection-convection model with insulating boundary losing mass?

I am trying to model a 1-d advection-convection numerically, using an upwind scheme. I'm using the following equation to calculate the value of internal cells: $$C_x^{t+1} = C_x^{t} + D\frac{\Delta t}{...
0
votes
0answers
89 views

Boundary conditions for an FEM approximation of the Laplace operator

Using FEM, I want to approximate the Laplacian $$u = \nabla \cdot \nabla h \, ,$$ where $h(x,y)$ is an FEM approximated scalar field on the same mesh, i.e. piecewise differentiable. I am using MOOSE ...
2
votes
0answers
29 views

Semi-analytical/empirical modelling of wall boundary conditions in advection-diffusion-reaction equation with distributed source

Let's suppose I need to numerically solve a 3D steady-state transport equation of the form $$ \nabla \cdot (\mathbf{u} c) = \nabla \cdot (D \nabla c) - \lambda c + S $$ where $c$ is the transported ...
0
votes
0answers
27 views

(FD WENO) Correct symmetry boundary condition for Euler equations

I'm trying to solve 2D Euler equations in axisymmetric formulation with finite-difference WENO scheme. I found some info on high-order boundary conditions for plane formulation (in this thesis, for ...
1
vote
0answers
57 views

Solving Laplace equation with constraint on boundary

I have found the following PDE problem in a paper: Essentially, we have a rectangular domain where there is a unknown interface $z=\xi(x)$ (liquid-air interface) separating the domain into two medium ...
2
votes
0answers
87 views

How are finite volume method boundary conditions implemented without using ghost-cells?

I'm currently trying to implement my own FVM code in cpp, but when I try to calculate the laplacian of a test function, given by \begin{align}\phi_0=\sin(2\pi x)\sin(2\pi y)\end{align}, I get ...
0
votes
0answers
126 views

Poisson equation, stiffness matrix positive definiteness, Dirichlet boundary conditions

I have a question regarding the positive definiteness of the stiffness matrix. Specifically, I believe that it should be positive definite only when at least one Dirichlet point is given, so I would ...
2
votes
0answers
43 views

How to solve this boundary value problem which has more unknown than equation on MATLAB

I need your helps about solving the problem below with MATLAB. I am trying to solve 2D Stress Wave Propagation problem by using FDTD(Finite difference time domain) method on the cylindrical coord. I ...
0
votes
0answers
47 views

Cauchy problem ill-posed?

Find the solution to the Cauchy problem consisting of the wave equation : $$u_{xx}-u_{yy}=0$$ together with initial conditions: $$ u(x,0)=0,$$ $$u_{y}(x,0)=g(x)$$ for some known initial datum $g$. Is ...
2
votes
2answers
153 views

Two-dimensional heat equation with Neumann boundary conditions: any hope to find an analytical solution?

I am looking for references showing how to analytically solve the heat equation with Neumann boundary conditions in two dimensions. So far, I have found the problem solved analytically in one ...
0
votes
1answer
98 views

Implementing Dirichlet BC for the Advection-Diffusion equation using a second-order Upwind Scheme finite difference discretization

i am implementing a Matlab code to solve the following equation numerically : $$ (\frac{\partial c}{\partial t} =-D_{e} \frac{\partial^2 c}{\partial z^2} +U_{z}\frac{\partial c}{\partial z}) $$ with ...
0
votes
0answers
155 views

"This DAE appears to be of index greater than 1" daeic12 (line76) error code

Hi I am trying to solve a set of pde converted into ODE and DAE using central finite difference method. I have used the MATLAB 'solve' command to determine the coefficients of fictitious nodes for ...
1
vote
0answers
48 views

Integrating a wavelike equation with absorbing boundary conditions

I am trying to numerically solve the following equation: $\frac{\partial^{2} \phi}{\partial t^{2}}-\frac{\partial^{2} \phi}{\partial x^{2}}+V(x) \phi(x, t)=0$ On some domain, with: $\phi(x, 0) = I(x)$ ...
6
votes
0answers
104 views

A Question About a Claim from 1991 Computational EM paper about the Cancellation of certain Boundary Terms

Please let me know if this is not the appropriate site for this question. I found questions regarding EFIE/MFIE/CFIE on this site, so I thought my question might fit. I am studying the paper by Putnam ...
1
vote
2answers
83 views

inclined/general Dirichlet boundary conditions

For simpilcity, consider a single quad linear elasticity finite element in 2D. The Dirichlet boundary conditions on node 1 and node 2 are easy to implement and can be handled in the standard way. ...
0
votes
1answer
70 views

Are spatial boundary conditions required for PDEs discretized with Method of Lines?

As far as I understand, you need to define boundary conditions in time and space to select a unique solution to a PDE and make it solvable. However, in ODEs I only need to specific the initial value ...
1
vote
1answer
63 views

Where could error terms that blow up in SWE come from?

I have been working on a solver for shallow water equations with reflective boundary conditions. I have found that it diverges very fast. As a workaround I noticed yesterday that if I smooth the ...
0
votes
0answers
96 views

How to use FEniCS to calculate the electric field of an isolated charged sphere

Initially I thought that this is the kind of question which ought to have already been answered in the form of an example online, but so far I haven't found one. I will admit that I am very new to ...
2
votes
1answer
128 views

How to apply Dirichlet boundary conditions to time-dependent PDEs?

Assume the time-dependent linear elasticity equation. Using a finite element discretization we obtain $$M\ddot{u}=Ku+F_\text{ext}$$ where $M$ is the mass matrix,$K$ is the stiffness matrix, and $F_\...
2
votes
1answer
175 views

Finite element (1D) for steady state non-linear problem

I need to solve with linear finite elements the equation $$\frac{\partial }{\partial x}\Bigl(\text{sgn}(x) u \Big) +\frac{\partial}{\partial x} \Bigl[ \sqrt{u} \frac{\partial u}{\partial x} \Bigr] =0$...
0
votes
2answers
365 views

No flux Neumann boundary condition for non-stationary PDE equivalent to Dirichlet boundary?

When using no flux Neumann boundary conditions (i.e. zero derivative to the normal on the boundary) in a non-stationary PDE, I don't seem to recognize the difference to using Dirichlet boundary ...
4
votes
1answer
224 views

Imposing no flux boundary conditions on variants of the Cahn-Hilliard equation using Finite Differences in Python

I have been looking into simulations of phase separation in variants of the Cahn-Hilliard system and have been running into issues with implementing no flux boundary conditions on certain variants. ...
1
vote
0answers
39 views

Discretization formula for a system of two differential equations. "Solution to one of these is the initial condition of the other". In which sense?

Consider the following stochastic differential equation \begin{equation} dy=\left(A-\left(A+B\right)y\right)dt+C\sqrt{y\left(1-y\right)}dW\tag{1} \end{equation} where $A$, $B$ and $C$ are parameters ...
1
vote
1answer
112 views

traction boundary conditions in elasticity

I have a question about implementing traction boundary conditions in 2D and 3D linear elasticity. Consider the picture above. I want to apply traction boundary conditions on the boundary in red. My ...
0
votes
1answer
115 views

Which scheme for inhomogeneous convection-diffusion equation with highly variable coefficients?

I have a 1D convection-diffusion equation $\sigma_t = a(x,t) \sigma_{xx}+b(x,t)\sigma_x+f(x,t)$ defined on the unit interval, with nonzero Neumann boundary conditions at both ends. It should be noted ...
1
vote
0answers
43 views

Euler Explict Method for solving the PDE for option prices under the Schwartz mean reverting model. Numerical Finance

I have to solve the following PDE for a Call option : $\partial_tV + \{ \alpha - (\mu - \lambda/ \alpha -log(S))\}S\partial_SV + 1/2 \sigma^{2}S^{2}\partial_{S}^{2}V - rV = 0$ Where V(S,t) is the ...
3
votes
2answers
310 views

Finite difference method having a discontinuity

I am trying to understand the FDM which is a widely used method solving differential equations by using approximation below. $$\dfrac{\partial u}{\partial x}=\dfrac{u(i+1)-u(i-1)}{2\Delta x}$$ How can ...
1
vote
0answers
41 views

Effect of reducing flux consistency order at boundary on convergence order

Consider the 1D nonstationary convection-diffusion PDE $$ \begin{alignat}{2} \partial_t u &= -a \partial_x u + D \partial_{xx}u, &\qquad x \in (0,1), t \in (0,T), \\ f(t) &= \left.\left( a ...
0
votes
0answers
30 views

At what l/d ratio will a frame element start to behave as a shell element?

I'm working in ETABS. There are few columns of dimension 300mmx1400mm. The height of building is 36.6 meter above ground level and the dimension of building is 26mx68m. I'm getting the time period of ...
2
votes
1answer
128 views

Methodology Suggestion for Wave-propagation Problem using Finite Elements

I want to simulate the propagation of a sinusoidal plane wave in a rectangular domain using Finite Elements Method. First, the wave should propagate through a fluid medium, then it will encounter a ...
0
votes
0answers
65 views

COMSOL Circularl polarization

I'm having some problems trying to implement circularly polarized light in COMSOL Muliphysics. For a isotropic homogenous media, I've obtained without problems the TE and TM reflectance curves. ...
4
votes
2answers
529 views

Kronecker product representation of the finite difference laplacian

The laplacian equation when discretized gives a system of linear equations that can then be solved. See the answer to this question: https://math.stackexchange.com/questions/3120948/discretization-...
2
votes
1answer
619 views

How can I evaluate more accurate energy eigenvalues from Schrodinger equation using shooting method?

I am trying to use the "shooting method" for solving Schrodinger's equation for a reasonably arbitrary potential in 1D. But the eigenvalues so evaluated in the case of potentials that do not have hard ...
2
votes
1answer
247 views

Boundary conditions in a finite element eigenvalue problem

I've been reading multiple papers and related posts for a while now, but I can't seem to find a specific answer to the issues I'm having so I hope someone can clarify things here. I'll provide some ...
0
votes
0answers
77 views

Applying boundary Conditions on FEM

I have a partial differential equatons as shown below. $$\dfrac{d}{dx}((1+x)\dfrac{du(x)}{dx})=0$$ With the following boundary conditions. $$u(0)=0, u(3)=10$$ To solve it using FEM, I multiplied the ...
0
votes
1answer
97 views

Applying weak form

I have two dimensional equation and I want to solve it using Finite Element Methods. $$ \nabla . (\alpha(x,y)\nabla u(x,y)) + \dfrac{\partial u(x,y)}{\partial x}+\dfrac{\partial u(x,y)}{\partial y}+u(...
3
votes
2answers
165 views

Solving a Boundary Value Question $\frac{d^2y}{dx^2}=y\cos(x)+\frac{\sin(x)}{x^2+2}$ using Python

I'm looking to solve this boundary value question using the shooting method! $$\frac{d^2y}{dx^2}=y\cos(x)+\frac{\sin(x)}{x^2+2}$$ given the initial values: $$y'(x=-1)=-1\\y'(x=5)=0$$ I'm aware of ...
0
votes
1answer
98 views

Determining Displacement Field on a Sphere

I need to find the displacement field for this sphere shape in terms of $\delta$. So far, by applying boundary conditions, I know $u_r = u_\theta = u_\phi = 0$ at $r = a$ and $u_r = \delta, u_\theta = ...

1
2 3 4 5
8