Skip to main content

All Questions

Filter by
Sorted by
Tagged with
0 votes
0 answers
35 views

Help Understanding the Weak Equation for Fully Developed Flow (Inlet) in COMSOL

I'm working on a simulation in COMSOL Multiphysics and trying to understand the Fully Developed Flow (Inlet) boundary condition. The documentation mentions: The Fully Developed Flow boundary condition ...
PDEPadawan's user avatar
0 votes
0 answers
29 views

Combining Study Steps with Multiple Geometries in COMSOL

Currently, I am using COMSOL to first solve an eigenmode problem in a simple element and subsequently extracting the displacement from the surface of said element and the results from the eigenmode-...
Benvz's user avatar
  • 1
0 votes
0 answers
128 views

Heat Equation for fast source with FiPy

I'm trying to solve the following differential equation with FiPy, basically laser irradiation on a surface $$ \rho_{s}C_{p,s}\frac{\partial T}{\partial t} = k_{s}\frac{\partial^{2}T}{\partial x^{2}} +...
clope99's user avatar
  • 11
3 votes
0 answers
53 views

Datasets for inverse heat transfer problems

I was wondering if there is an available, real-life known inverse heat transfer problem dataset to benchmark oneselfs algorithm, as in MNIST for deep learning. Talking about... (well in this case I ...
Aner's user avatar
  • 191
1 vote
1 answer
256 views

Coupled Partial Differential Equations

I'm trying to solve the following system of coupled differential equations, the two-temperature model for $e$ = electrons and $l$ = lattice. $$ \rho_{e}C_{p,e}\frac{\partial T_{e}}{\partial t} = k_{e}\...
clope99's user avatar
  • 11
1 vote
0 answers
50 views

How can I apply a mixed boundary condition to a multi-material heat transfer problem using Crank-Nicolson?

I am working on a mixed material model for a melting material and need to enforce both a Dirichlet and Neumann type condition at the interface. Subject to an external surface heat flux at the top of ...
ZebraEagle's user avatar
0 votes
0 answers
38 views

Thermo Hydraulic Mechanical modeling of energy wall slab in Comsol multiphysics

I am currently working on a complex simulation project involving an energy wall slab, and I need assistance in accurately modeling and validating it using COMSOL Multiphysics. Here are the details of ...
Hizbullah's user avatar
0 votes
0 answers
41 views

Fermi level determination

Is there a semiconductor physics module in COMSOL which correctly describes the Fermi level position for any bandgap system ( especially under different doping concentration)?? If not what is solving ...
Mohsen's user avatar
  • 1
1 vote
2 answers
121 views

How to handle non bilinear weak form?

I solved the 2D heat equation using the finite element method. It all went well first with the adiabatic case, however problems occured when I introduced cooling with the enviroment. I modeled the ...
Boiler4562's user avatar
2 votes
1 answer
411 views

Where am I making a mistake in solving the heat equation using the spectral method (Chebyshev's differentiation matrix)?

I would like to numerically solve the following heat equation problem: $$ u_t = \Bigg(2{a \over l}\Bigg)^2 u_{xx} \tag 1$$ $$ x \in [ -1, 1 ] \tag 2$$ $$ u(x, 0) = 0 \tag 3$$ $$ u(1, t) = A \sin \Bigg(...
FriendlyNeighborhoodEngineer's user avatar
1 vote
0 answers
41 views

What is the physical quantities in the equations for linear elasitc materials in a frequency study in comsol

I am using comsol to model a linear elastic materials and perform a frequency domain study. I try to understand the equations it provided but all symbols are not explained. For example, How should or ...
lsdragon's user avatar
  • 111
1 vote
1 answer
339 views

Why does scipy Conjugate Gradient solver fail to converge for non-steady heat equation using Crank-Nicolson method

Could someone please explain why my implementation of the Crank-Nicolson method applied to the non-steady heat equation won't converge? There shouldn't be any nonlinear aspects to my implementation ...
n1ck94's user avatar
  • 13
2 votes
0 answers
114 views

Efficient heat diffusion implementation with varying coefficients

I have the following heat diffusion equation: \begin{alignat}{3} \partial_t u(t, \vec{x}) &= g(\vec{x})\Delta u(t,\vec{x}), &\quad& \vec{x} \in\Omega, \, t\in(0,\infty],\\ \partial_n u(t,\...
lightxbulb's user avatar
  • 2,892
5 votes
1 answer
107 views

Prediction of sphere (i.e. roast) core temperature heated in an oven

The real-life problem Assume I put a spherical roast with initially constant temperature of start_temp=25 (°C) into an oven with ...
Dieter Menne's user avatar
1 vote
1 answer
313 views

2D Heat equation solved with finite element method converges in skewed way

I tried to solve the 2D heat equation with the finite element method, using triangles as elements. Currently generated by a Delaunay triangulation. The base function I'm currently using is basically ...
Boiler4562's user avatar
0 votes
1 answer
143 views

How to get a normalized gradient with FreeFem++?

I am trying to use FreeFem++ to solve the heat geodesics algorithm. The algorithm is: solve $\dot u = \Delta u$ at a specific time $t$. compute $X = \frac{\nabla u_t}{|\nabla u_t|}$ solve $\Delta\phi ...
Makogan's user avatar
  • 379
3 votes
0 answers
100 views

Is there any book about fundamental FEM theory using similar terminology as Comsol MultiPhysics?

I am considering to solve complexed PDE systems, like in this post, using Comsol MultiPhysics. The PDEs are different from the General Form provided by Comsol. The Weak Form module may be worth a ...
Charles6's user avatar
1 vote
0 answers
115 views

Accuracy of the Crank-Nicolson method for non-linear, inhomogeneous heat equation

I am currently coding a solution to the following PDE: $\frac{\partial T }{\partial t} =\frac{\partial}{\partial \theta}(A(\theta ,\phi )\frac{\partial T }{\partial \theta}) +\frac{\partial }{\partial ...
mathbruh67's user avatar
2 votes
0 answers
113 views

Conceptual doubt regarding 2D conjugate heat transfer modelling (COMSOL and Mathemtica)

I have been dealing with some conceptual flaws in my understanding of modelling, which I will elaborate herein. I am modelling conjugate heat transfer of a reciprocating fluid, which flows with ...
Avrana's user avatar
  • 41
0 votes
1 answer
148 views

Problem with my Octave code (unsteady heat equation with FEM)

I want help with my Octave code regarding the unsteady heat equation. My geometry and mesh are generated with FreeFEM++, so there is no problem with that (I tried it with the steady problem with no ...
mibo27's user avatar
  • 3
1 vote
2 answers
388 views

What is the difference between non-linear elastic simulation and linear elastic simulation with plasticity?

I'm learning how to do Finite Element calculations using Comsol Multiphysics. In Comsol, Linear Elastic Material and Nonlinear Elastic Material are available as material models: Using Linear Elastic ...
S. Rotos's user avatar
  • 133
2 votes
0 answers
92 views

How to accelerate the computing of implicit finite difference method for heat conduction between two solids

Edit on May 3rd: I have found the problem. Because the difference of between $k_1$ and $k_2$ is huge, a very small time step need to be chosen so that the right green part can "feel" the ...
Kai Jiao's user avatar
2 votes
0 answers
294 views

Rosenthal equation for multi track

Rosenthal's equation lets one calculate the temperature profile of a moving point heat source analytically for thin and thick plates. For simplicity I use the equation for thick plates defined as: $$T-...
vydesaster's user avatar
2 votes
0 answers
150 views

switching boundary conditions to simulate a reciprocating flow in COMSOL

This is a re-posted question which was poorly defined earlier I have been trying to simulate reciprocating flow of a convective fluid through a heated channel in COMSOL-Multiphysics. The schematic of ...
Avrana's user avatar
  • 41
1 vote
1 answer
2k views

How to solve heat equation in spherical coordinates with finite differences?

I have a problem dealing with heat transfer which is spherically symmetrical. I was thinking it should be possible to solve this as a 1d problem in spherical coordinates using the radius only. ...
Alex I's user avatar
  • 111
1 vote
0 answers
111 views

Weird "oscillatory" modes appearing in FEM simulations

I am using COMSOL to solve a mathematical model involving thermoelectric hydrodynamic (TEMHD) flow. I am running a very large parameter sweep and using the solutions obtained to make some plots. ...
epsilonD3LT4's user avatar
4 votes
1 answer
139 views

Solving geodesics on triangular meshes gives negative distances

I have implemented the heat method for geodesics: https://www.cs.cmu.edu/~kmcrane/Projects/HeatMethod/paperCACM.pdf When I run it I am getting a solution that, visually, seems correct: In this image, ...
Makogan's user avatar
  • 379
0 votes
1 answer
214 views

Incorporating heat flux into Laplace Equation

I need to find the temperature distribution of a square plate using the Laplace equation by using FDM: $$ \frac{d^2T}{dx^2} + \frac{d^2T}{dy^2} = 0$$ But there is a heat flux entering from the top ...
justauser's user avatar
  • 145
2 votes
2 answers
468 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 ...
Max_89's user avatar
  • 61
2 votes
2 answers
219 views

Heat equation in non-dimensional form behaving differently than in usual format

Starting from $$ c_p \frac{\partial u }{\partial t} = k \nabla^2 u $$ in a one dimensional domain [0,1] where $c_p$ and $k$ are modeling two different materials: $$ k = \begin{cases} 1 ~\text{if} ~x &...
balborian's user avatar
  • 601
0 votes
0 answers
54 views

Why does COMSOL treat a well as a rectangular cube?

I'm trying to learn COMSOL for a graduate research project and am struggling through a flow and transport in 3D tutorial. As such I'm walking through my boundary conditions and trying to figure out if ...
ldorchester12's user avatar
1 vote
2 answers
812 views

How to use the Thomas-Algorithm to the Heat-diffusion-equation correctly

My post is structured in four parts: I give you some information about the context my principal questions refer to. I will tell you what I believe to know about the Thomas Algorithm. If I am wrong ...
tom terenius's user avatar
2 votes
3 answers
553 views

Flux sign and face normal confusion in finite volume method

I implemented a solver for the 2D steady-state heat equation (without heat generation and homogeneous material) $\nabla. (k\nabla T) = 0$, using finite volume method, however, I am having some ...
Algo's user avatar
  • 304
0 votes
1 answer
163 views

How to connect two cylinders to form a knee in Comsol Multiphysics?

I have this I want it to be single bended wire. How to accomplish?
Dims's user avatar
  • 111
1 vote
2 answers
384 views

Writing a single PDE from a gradient equation

I have a differential system like this, where $\Phi$ is a scalar valued unknown function: $$\nabla\Phi = \left(f_1(x, y), f_2(x,y)\right)^T$$ I'm trying to solve it in a FEM solver (COMSOL ...
VP06's user avatar
  • 11
0 votes
0 answers
336 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. ...
Antonio Ganfornina Andrades's user avatar
1 vote
0 answers
456 views

Incorporating radiation boundary condition at the edge in finite difference

I am trying to solve the 2-d heat equation on a rectangle using finite difference method. I am confused as to how to incorporate non linear radiation boundary condition at the edge. $-k\frac{\partial ...
Nitin's user avatar
  • 19
3 votes
0 answers
108 views

Correct approach for thermal finite element simulation of layered assembly

I would like to optimise the heat transfer on a PCB. Several dies are on the top and cooling air is going through the fins in heat sink on the bottom. The assembly consists of several layers like ...
Ken Grimes's user avatar
2 votes
1 answer
434 views

Solution method of nonlinear heat transfer analysis

The governing equation of transient heat transfer analysis is described as follows: $$C \frac{dT}{dt}+K T = Q$$ When using backward difference scheme for the discretization of the time we get the ...
vydesaster's user avatar
0 votes
1 answer
51 views

Produce vertex displacements from volumetric shrinkage data on unstructured meshes

I was wondering what would be an efficient way to produce compatible displacements for mesh nodes/vertices if the computed data is volume shrinkage of each element/cell in the unstructured mesh? ...
Johntra Volta's user avatar
2 votes
1 answer
280 views

Does mass balance hold in convective diffusion

I'm trying to understand how convection-diffusion equations are solved in pipe flow modules available in CFD solvers. $$ \frac{\partial C}{\partial t} + \nabla \cdot (\mathbf{v} C) = \nabla \cdot (D \...
Natasha's user avatar
  • 433
0 votes
1 answer
219 views

A lot of identical staff in Comsol material database?

I got a lot of elements in Material Browser of Comsol Multiphysics of Optics section. ...
Dims's user avatar
  • 111
1 vote
0 answers
46 views

Fitting a multivariate PDE (using Java)

I'm doing simulations of 2 coupled PDE's with Comsol Multiphysics. I want to fit some data (using the Application method, whose language is Java) to those simulations. In order to answer my question ...
J.A's user avatar
  • 171
2 votes
1 answer
311 views

Lumped matrices in thermal analysis using finite elements

The governing equation of the transient heat transfer problem is $$C \frac{dT}{dt}+K T = Q$$ $C$ is the heat capacity matrix. $K$ is the thermal conductivity matrix. $T$ is the temperature vector. $...
vydesaster's user avatar
0 votes
0 answers
112 views

Simulating flow in a branched pipe

I am trying to simulate 1D advection and convection of a solute in the following blood vessel segment. I would like to know if this system can be simulated in COMSOL or MATLAB. I have used pdepe ...
Natasha's user avatar
  • 433
1 vote
1 answer
548 views

For a determined (known) Space charge density, what are the conditions to obtain the Electric potential/field distribution? (COMSOL, MATLAB)

Theoretic part From the theory, in Electrostatics inside a real dielectric material between real conductors, in a simple 1D plane geometry between points $P1$ and $P2$, according to the current ...
Lucian Viorel's user avatar
2 votes
1 answer
448 views

Derivatives over a Finite Element mesh

I have a data extracted from Comsol on some node points and I know the coordinates of each node. Does anyone know how Comsol calculate the partial derivative from the values at each node and also ...
starhd's user avatar
  • 23
1 vote
2 answers
191 views

Simulating the heat equation with insulating material

My plan is to solve the heat equation in the right half portion of the domain, while having the left half completely isolated with constant temperature. To do so, I model the left half with a very low ...
balborian's user avatar
  • 601
0 votes
1 answer
737 views

How to set an initial guess for the iterative solver in Comsol?

How to set the initial guess for the iterative solver GMRES or FGMRES for linear problems (Helmholtz equation of RF module) in Comsol?
praveen kr's user avatar
3 votes
0 answers
135 views

Strain propagation from surface to bulk in COMSOL

I am trying to simulate strain propagation from the surface into the bulk. I have a rectangular semiconductor block (~2 μm thick) on top of which metal gates (~25 nm thick) are deposited as seen in ...
duncanidaho's user avatar