Latest Discussions - COMSOL Forums

Demagnetization of N52 magnets in a Halbach configuration

Mon, 11 Aug 2025 16:45:54 +0000

Hello,

I'm currently working on modeling the demagnetization of an array of N52 magnets in a cylindrical Halbach configuration using the AC/DC module. Currently I'm using the remnant magnetization model for the N52 magnets using COMSOL's materials library, but I am worried that with the field strengths being produced (well over 2 Tesla) that this system is operating far outside the linear regime for these magnets. I've attempted to model the magnets using the nonlinear permanent magnet model, with a B-H curve experimentally measured for the magnets supplied by the manufacturer, but running this simulation produces no field at all. Any help in modeling the non-linear behavior would be appreciated.

Messages: Failed to initialize physics interface

Mon, 11 Aug 2025 11:19:55 +0000

When I try to load a COMSOL file into Python using "client = mph.start()" and "model = client.load("./test.mph")", I encounter the following error. Some files load without any issues, while others fail. I suspect it is not a license or version problem, but rather an issue caused by specific settings within the file itself. "Even a blank file with no components, only physics(heat transfer in solids) added causes an error when loading.

Traceback (most recent call last): File "java.lang.Thread.java", line 0, in java.lang.Thread.run

Exception: Java Exception

com.comsol.util.exceptions.com.comsol.util.exceptions.FlException: Exception: com.comsol.util.exceptions.FlException: Failed to initialize physics interface. Messages: Failed to initialize physics interface.

Transient ray heating simulation in COMSOL (Ray Tracing + Heat Transfer in Solids)

Mon, 11 Aug 2025 10:58:28 +0000

Hello everyone,

I am trying to simulate the heating of a glass plate using Ray Tracing and Heat Transfer in Solids in combination with the Multiphysics node “Ray Heating”. The rays from the ray tracing act as the heat source.

My goal is to run a time-dependent study where the glass plate is heated by a pulsed laser with 10 ms ON time and 90 ms OFF time (repeating), which is essentially a quasi-CW operation mode.

However, I have the following problem: Single study approach

If I activate all physics in one study (see attached Model Tree / Study screenshot), the time-dependent solver does not allow me to resolve the heating process in the millisecond range. In this case, the rays appear to stop at the glass plate and the heat transfer is not computed as expected.

Two-study approach

I also tried separating the simulation into two studies: Study 1: Ray tracing (stationary) Study 2: Heat transfer in solids (time-dependent), using the Values of Dependent Variables from Study 1 as initial data.

Unfortunately, no heating occurs in the second study, even though the ray tracing results contain nonzero gop.I and gop.Q values.alues.

Has anyone encountered this behavior or knows how to properly set up a transient simulation with pulsed ray heating in COMSOL? Any hints on how to correctly couple the ray tracing results to the heat transfer simulation over time would be greatly appreciated.

Thank you in advance for your help!

Releasing light from a moving particle in solution?

Fri, 08 Aug 2025 22:10:55 +0000

Hello, I'm attempting to model light emission from a phosphorescent particle in a fluid flow. The flow and the particle motion are no problem. Modeling the release of light however is proving to be difficult. My thought was to solve the fluid flow in a stationary manner, then use a segregated solver to first update the position of the particle within the fluid, followed by updating the position of the ray release point based on that new particle position. And that's where I'm stuck and am open to suggestions. Thanks.

Method for determining the mobility tuning parameter in phase-field modeling of two-phase flow

Fri, 08 Aug 2025 14:35:21 +0000

In the Technical Support Knowledge Base, it is stated that when simulating two-phase flow problems using the phase-field method, the mobility tuning parameter can be calculated based on the velocity (as shown in the figure), rather than on the fluid viscosity, interface thickness parameter, and characteristic length scale. Is there any published reference that demonstrates the reliability of this algorithm?

Simulating Residual Stress of Thin Films on Silicon Wafer by Thickness Changing

Fri, 08 Aug 2025 12:22:10 +0000

Hello!

I am trying to simulate the stress (or strain) that is obtained from

the fabrication process, and

the experimental process (temperature evaluation).

Initially, there is a silicon wafer. On top of it, there is a supporting SiO₂ layer. On this support layer, VO₂ is deposited using a sputtering system. I assume that all thin films inherently have an initial intrinsic stress from the deposition process. Based on previous studies, I have found that the magnitude of this stress depends on the difference in lattice parameters between the substrate and the film, as well as on the film thickness. The fabrication process is performed at room temperature.

Moreover, after deposition, the film is tested by measuring its resistance while the temperature is increased from 25 °C to 100 °C. I would like to observe how the stress/strain in the film changes as the temperature varies.

Currently, I am using COMSOL 5.3, and my procedure consists of the following steps:

Fabrication Process (Goal: To show the initial stress distribution across the structure/geometry)

I assume that the wafer has negligible (zero) stress initially. After the second film is deposited, I use the Solid Mechanics module with a Stationary study.

I add another layer (m1) on top of the silicon wafer in the geometry.

I set the initial stress and initial strain in the silicon to zero.

I set the initial stress in the deposited layer to the given intrinsic stress, and the initial strain in the deposited layer to the calculated value from the lattice parameter difference:

(lattice parameter difference, such as strain=(a_m1 - a_film1)/(a_film1) ).

I would like to simulate cases where the top VO₂ layer has thicknesses of 10 nm, 30 nm, 50 nm, 70 nm, and 100 nm, and to show both:

the average stress applied over the entire surface area, and

the local stress distribution.

There is no temperature change and no external force or strain applied to the film. So, this simulation should ideally show the stress distribution that develops in the geometry due to the deposition process — specifically, the residual stress along with the initial intrinsic stress.

In addition, I would like to visualize stress vectors — arrows indicating the magnitude and direction of stress. I am not sure how to do this. I tried using an arrow array with (d.stress1x, d.stress1y), but the arrow directions did not match the expected theoretical results.

My questions for this step are:

Can the von Mises stress values obtained from this simulation be considered valid evidence of the residual stress formed in the thin film?

Are there any additional assumptions or considerations that should be added to make this simulation more realistic?

Experimental Process (Temperature Increase) I use the Solid Mechanics and Thermal Stress modules with a Stationary study.

Assuming the same structure as in Step 1 (Si + SiO₂ + VO₂), I apply a temperature change from 25 °C to 100 °C to the thin-film system.

The boundary condition is set to room temperature.

The initial temperature is set to 100 °C.

For the same VO₂ thicknesses (10 nm, 30 nm, 50 nm, 70 nm, 100 nm), I want to observe the average stress over the surface area and the local stress distribution as the temperature changes.

I would like to know:

Is this the correct way to set the temperature change in COMSOL?

What exactly is the meaning of reference temperature and ambient temperature in COMSOL?

General Question I would like to know if this is the best approach to obtain the stresses in the layers after a fabrication process. Is it possible to automatically use the simulation results from one step as the initial conditions for the next process step? For example, after simulating the first fabrication step, could I then add a new layer to the geometry and assign it an initial intrinsic stress while preserving the stress state of the existing layers from the previous step?

Thank you!

If you’d like, I can also prepare a clean, concise version optimized for COMSOL’s discussion forum style so it’s both detailed and easy for others to respond with technical advice. That would make it much more likely to get expert replies.

Need Help with Suction-Based Climbing Robot Simulation

Thu, 07 Aug 2025 20:38:32 +0000

Hi everyone,

I'm working on simulating a suction-based climbing robot using COMSOL Multiphysics. My objective is to check whether the robot will stick to a vertical surface or slip based on applied suction.

Project Details:

Suction pressure: –15,750 kPa applied to each suction pad

Robot weight: 117 N

I want to simulate how the robot behaves under gravity and suction force—whether it stays adhered or begins to slide.

I'm confused about contact pairs and setting up proper interaction between the robot body and the wall.

Also, deciding between stationary and time-dependent studies is a bit unclear.

Pipe flow and pump

Thu, 07 Aug 2025 13:12:27 +0000

Hi there, can someone please share how to set a pump with inlet and outlet of a pipe flow model (gas, not liquid) with the goal of simulating the pressure losses due to elbows and pipe length, please. The flow rate and ultimate vacuum of the pump are set, trying to work out how to achieve this flow rate reduction / pressure difference drop when piping gets very long. thanks for you help. I found the liquid example in pipe, but not anything with pump and gas. any reference to some documentation is appreciated.

tcd.kappaxx-zz in the Tertiary Current Distribution governing Equations

Thu, 07 Aug 2025 12:43:27 +0000

Hi All,

I took a deep dive on porous electrode reaction node equation view and i saw electrolyte conductivity term tcd.kappaxx/xy/xz/yx/.../zz. I tried find and replace (ctrl + F) for this variable but i could not find what this variable do and how does this contribute to the overall tertiary current distribution physics.

does anyone know how COMSOL uses this term? or is it just a post-processing variable?

How to stimulate the electrode on the semiconductor module?

Wed, 06 Aug 2025 05:57:05 +0000

Hi guys,

I want to stimulate a ITO-ETL-PVSK-HTL-FTO Perovskite solar cells by using the semiconductor module of COMSOL to investigate quantities like electrical potential, concentration of holes and electrons, and etcs, and I want all layers are included. Now I encounter a problem, which is I couldn't include FTO in my model although I can include ITO.

The first thing came into my mind for this issue is that maybe it is due to the PN effects on HTL(P)-FTO(N), such that there exists an large energy level gap on their interface, so the holes in HTL could hardly move into FTO, which then fail to solve the continuity equation on that boundary. However, after some adjustments on the electric affinity of FTO and HTL, I was almost close that gap(left a <0.5eV energy gap), the stimulation failed again.

As for the other b.c.s and i.c.s, I set 0V on the top of ITO, and a scanned volt Vapp on the bottom of FTO, and doping ITO/FTO with 1e20[1/cm^3] in N-type, HTL with 1e18[1/cm^3] in P-type, and ETL with 1e18[1/cm^3] in N-type.

Thanks Pals

How to define the contact angle in a model

Mon, 04 Aug 2025 11:29:20 +0000

Hello there, I am trying to simulate droplet impact on solid stationary particle. I am finding it difficult to define contact angle. I used wetting wall BC to give values of contact angle but it doesn't matter whatever value i give in contact angle, I am not able to see it in output, it feels like the contact angle value which i give has no effect upon simulation. Please guide me where i am doing wrong. Thanks

PML Wavelength

Sun, 03 Aug 2025 22:02:27 +0000

Hello I'm doing a structural analysis in Comsol using PML. The analysis is in the frequency domain, and I'm unsure which wavelength to use in PML.

Another question is, when I select the Comsol option to use the wavelength based on the physics of the problem, is there any way I can access that wavelength?

Schur complement of the global stiffness matrix

Sun, 03 Aug 2025 07:46:19 +0000

I am trying to calculate the Schur complement of the global stiffness matrix, I have already exported the matrix through the livelink. To calculate the Schur complement of the matrix, I have to figure out which part of the matrix corresponds to which node, therefore, are there any way for me to get how comsol assemble the matrix? Besides, I know that comsol has a schur complement solver, are there any way for me to perform this on the global matrix? (e.g. I got two same cascaded units, and I only need to calculate one stiffness matrix and cascaded them through the schur complement)

Silica Precipitation

Sun, 03 Aug 2025 02:50:44 +0000

I'd like to know how to achieve discontinuous precipitation of silica. My idea is to use the dissolution and precipitation of silica to change the porosity and permeability of the fractures, thereby achieving a discontinuous precipitation effect. I'd like to ask if this can be achieved.

How to integrate over infinite domain or scaling coordinate system

Sat, 02 Aug 2025 19:03:24 +0000

Hello!

Every integration I tried didn't recognize or consider the stretching of an infinite or scaling domain. What is the right or appropriate function or way to integrate over stretched scales? Whether in Results or in Component.

Thank you.

P.S. Datasets->Emped 1D/2D seemed a good way or workaround, but the 'Output space expression' doesn't recognize used component functions or expression operators, and I don't want to use redundancy or a workaround.

Nice greetings

Non Hermitian Acoustic Model

Sat, 02 Aug 2025 00:09:48 +0000

I'm summarizing my whole work.

With COMSOL 6.2. Non-Hermitian Study with Balanced Loss and gain with only COMSOL Weak Form PDE. I have an Hamiltonian, H=[0 (v+gamma)+w.exp(-ik);(v-gamma)+w.exp(-ik) 0]. I tried weak form of equation.

Model Details: 1. Parameter: L=60[mm], kx=pi/L, h1=5[mm], r=5[mm] g=1[1/m], c_air=343[m/s], rho_air=1[kg/m^3], k=1 [1/m].

Geometry: Rectangle1 (Width L/2, Height h1, Position, x=0, y=0); Rectangle2 (Width L/2, Height h1, Position, x=L/2, y=0); Circle1(Radius r, Position, x=L/4 y=h1/2), Circle2(Radius r, Position, x=3L/4 y=h1/2). Then Union. Then delete Entities to delete internal boundaries inside the circle.
Materials: Now we have 5 domains. Domain 2 is the channel between between circles A (domain 4) and B (domain 5). We used user defined materials. Material Type: Non-Solid, Speed of sound, c=c_air and density, rho=rho_air. All the domains have same material properties.
In Component 1> Definitions>Identity Boundary Pair 1, I selected Source boundary as Left (1) boundary and Destination boundary as Right (10) Boundary.
Weak form of PDE: Applied to Domain all Domain. Weak Form PDE 1< Weak Expressions: - (test(px)px + test(py)py)+ ((2pifreq/c_air)^2)test(p)p.

Then took Weak Form PDE 2 and selected only Domain 2 and put expression: igtest(p)*px. This is non-reciprocity term.

Next, for Periodic Condition, I used Weak Form PDE (w) >Pair> Weak Contribution. Then at Pairs section selected the predefined "Identity Boundary Pair 1". And at Weak expression field, wrote expression as "test(p)(p - dest(p)exp(ikxL))".

Mesh: Free Triangular with Predefined size 'fine'.
Study 1: Parametric Sweep: g=0:0.01:1 and kx=0:0.1*pi/L:pi/L. Step1: Eigenvalue, Desired number of eigenvalues 6 and Search for eigenfrequencies around 1+1i[rad/s].

I run the simulation. There was no error. But, the periodic condition "test(p)(p - dest(p)exp(ikxL))" is not working. How to solve ?

Simulating Grain Evolution with Surface Heating in COMSOL

Fri, 01 Aug 2025 14:47:17 +0000

Hello,

I'm new to using COMSOL. I’ve already defined the geometry and assigned materials to a cubic sample. Now, I would like to simulate how the grain structure evolves when heat is applied to the outer surfaces of the cube.

Could you please guide me on how to set up this simulation?

Thank you. Nafish (PhD. Student)

NACA 0012 transient

Fri, 01 Aug 2025 12:55:32 +0000

hello dear community

i was trying to do time dependent study but unfortunately always show inconsistent initial values even though i followed the steps of the model that is provided by the comsol community i would like to ask for guide and what is the exact issue here

Heat transfer from housing to fluid

Fri, 01 Aug 2025 09:43:07 +0000

I've been struggling with what would seem a fairly simple heat transfer model.

I have a domain with a power loss. Above that is another domain which represents a housing, within which is a cooling channel. I've been trying to calculate the power dissipated to the fluid to correlate to another FEA package.

When I calculate dissipated power from temperature rise of the fluid, I get more power than what is introduced to the model from the power loss domain.

All walls are insulated, so all power should be going into the fluid minus the power to increase each component temperature.

I've increased mesh density but the overall result hasn't changed which leads me to believe the problem is more fundamental.

Can anyone share some insight?

What is the meaning of (nitf1.qwf_u)?

Thu, 31 Jul 2025 01:00:13 +0000

I just know this formula for calculate the power of exchange but i do not know the exact formula contained in it. I will appreciate it, please help me.