Latest Discussions - COMSOL Forums

Inverted mesh element - phase change interface deforming geometry

Tue, 01 Apr 2025 18:47:22 +0000

Hi,

I have been getting this error and this warning:

The following feature has encountered a problem: - Feature: Automatic Remeshing 1 (sol1/t1/ar1) Failed to find consistent initial values. Last time step is not converged.

Inverted mesh element near material frame coordinates (0.0499996, 87998.2). In Material Frame Variables: There was an error message from the linear solver. The relative error (1) is greater than the relative tolerance.

I am trying to create a moving boundary problem so my geometry is deforming, I am assuming this is what may be causing inverted mesh elements.

I have tried adding an automatic re-meshing node and refining my mesh around the deforming boundary.

I am new to COMSOL so any help would be appreciated, thanks.

Transient magnetic field: mimic a wire

Tue, 01 Apr 2025 18:46:45 +0000

Hello, I have built a solid mechanics and magnetic fields multiphysics simulation. At a high level, the simulation includes a current-carrying wire (that generates a magnetic field), and a small, soft-magnet device at some distance from the wire. In the presence of the magnetic field, I study the deformation of the soft magnetic structure.

I now want to do this simulation in a transient domain, however I would like to eliminate the need to simulate the field from the wire. Instead, I would like to specify or easily generate a field with the same exact profile as the wire. The reason for this is that COMSOL is giving me trouble when I attempt to study the coupled physics involving the current carrying wire while in a transient domain. If I can instead specify a transient magnetic field (not coming from a wire, just a "magic" field, if you'd like), then I can avoid this problem.

I have made this field already by using the MFNC module and specifying the Hx Hy and Hz. It runs fine. However, I am wondering if there is/are better ways to do this?

Thank you.

prism based SPR sensors simulation

Tue, 01 Apr 2025 09:02:20 +0000

Hello, I am doing the SPR sensor design on COMSOL, i have done my structure(Upside down the prism, and other layers onto the prism), assigned the materials, and properties also. Now I have to give the ports and Scattering Boundary Conditions. Can you help?

Problem With ODE Time Dependent Study

Mon, 31 Mar 2025 21:03:02 +0000

I have in my simulation some variables. They are converging correctly but when I add a ODE that is not coupled with those other variables, my solver doesn't converge anymore. All variables is not converging. Someone could explain to me what is happening? I'm using segregated solver, so, in theory this variable doesn't affect the other because they are uncoupled.

Running Batch Sweep with Batch commands

Mon, 31 Mar 2025 13:50:01 +0000

Hi All,

Is there a way to run a Batch Sweep using batch commands? I'm trying to execute a huge parametric sweep with Batch Sweep to speed up the overall simulation. I'd like to initiate the simulation with a batch command, but it seems to only perform a parametric simulation. (The simulation took the same amount of time as a normal parameter sweep, and the individual sub-simulations do not appear in the batch folder.)

Persistent error with undefined variable in time dependent solver and stationary solver

Mon, 31 Mar 2025 12:16:06 +0000

Hello Comsol,

I am facing an issue in my COMSOL model that I have been unable to resolve for a few days now, even after multiple attempts and fixes. I am working on a time-dependent simulation of blood flow in the aorta (normal aorta and with coarctation), and I have added Global ODEs (RCR model) to model the pressure changes in the arteries. My Boundary conditions have been set to velocity inlet and outlets with pressure. My main issue is that when I try to run the simulation, I keep encountering the error: "Undefined variable - Variable: comp1.p" This happens both in time-dependent and stationary solvers. How do I fix this issue?

(C_B*comp1.pt + comp1.p/R2_B - Q_B) this is the equation I have used in the global equation in global ODEs and DAEs. where: C_B is the compliance comp1.pt is dp/dt R2_B is distal resistance Q_B is the flow rate. I have defined all these parameters.

Electrical and Thermal Characteristics

Mon, 31 Mar 2025 09:32:08 +0000

Hello everyone! I'm looking for a big guy in Beijing who is proficient in electric heat simulation, and I can pay a consulting fee if I have a question!

keep all plot features in matlab livelink

Mon, 31 Mar 2025 05:08:55 +0000

There's a matlab file (function out = model) that creates a model, solves and postprocesses. All plot features including titles, legends are written within this model m function. Then use mphplot(model,'pg') to plot 'pg'. But this loses all plot features. How to use mphplot while keeping all plot features without rewriting all features?

solvers: solving the same model with various equally suitable modules

Mon, 31 Mar 2025 03:55:06 +0000

The software automatically specifies solvers according to the specific modules used. When solving the same model with various equally suitable modules in one batch within the same mph file, by default only one solver is used for all the modules applied. This solver may be suitable for only one of the modules but not for the rest of the modules. In this case how to make the software use different solvers when solving different modules?

Issue in modeling natural convection to a pot

Mon, 31 Mar 2025 01:54:16 +0000

I am trying to model the heat transfer to a pot in natural convection mode. This is an axisymmetric model, where a stove setup is modeled. The temperature at a bottom boundary (represents the burner) is set to a constant value. There is a pot placed on the stove. There is water (modeled as a solid object with a high thermal conductivity) inside the pot. We are studying the heat transfer to the pot and to the water through conduction and natural convection in this case.

I am using the heat transfer in solids and fluids, turbulent flow, and nonisothermal flow multiphysics coupling with a fully coupled solver to get the solution. The issue is that the bottom surface of the pot is getting heated, however, there is negligible temperature gain at the interior surfaces of the pot and also the water in the pot. What could be be possible issue? I am attaching the model .mph file for your reference.

Thanks!

Electrical and Thermal Characteristics

Mon, 31 Mar 2025 00:39:40 +0000

Hello everyone! Using a porous medium model to simulate Joule heat, it is known that the bed resistance decreases with increasing temperature, and the resistance and temperature data are obtained experimentally, fitted by an equation and then defined inside the material conductivity of comsol, why is the bed temperature change small?

Question about COMSOL tutorial provided

Sat, 29 Mar 2025 20:03:14 +0000

Hello I am trying to understand this tutorial https://www.comsol.com/model/dc-corona-discharge-in-air-in-a-point-to-plane-configuration-100451 I want to change the wall boundary condition to just absorb the electrons and i want to plot it on the wall. I am not able to do so. can you please tell me why. My objective is to determine how much charge got deposited on the wall during the simulation

different time stepping under same 1D plot group

Sat, 29 Mar 2025 04:45:18 +0000

I want to compute the same thing with 2 different time stepping and plot them under the same 1D plot group. The 2 different time stepping are plotted in line graph 1 (sol1) and line graph 2 (sol2) under 1D plot group.

Do I have to use different studies for different time stepping (yielding sol1 and sol2)? Different steps under the same study don't seem to allow different time stepping.

Is there any way of sweeping time steps in the parent node of 1D plot group? sol1 and sol2 are from different studies.

2D Simulation of PDMS Membrane with PEG-PCL-DA Hydrogel Reservoir

Fri, 28 Mar 2025 19:12:13 +0000

Hi,

I am trying to model the stress experienced within a PEG-PCL-DA hydrogel reservoir which will house cells. I am using a 2D, solid mechanics, time-dependent study to do so. I have created a rectangular membrane made of polydimethylsiloxane (PDMS) with a stationary end and a moving end, this end will be stretched at the strain rates 5%, 10%, and 12.5%. Within the rectangular membrane is a rectangular reservoir with the assigned material, PEG-PCL-DA hydrogel. The model is time dependent and is stretched for 0.2 seconds. I want to set up a cyclic time dependent study to include a retraction phase for 0.4 seconds after reaching a distance of 2mm, but am unsure about how to set up a piecewise or interpolation function. Any thoughts or advice? I have attached my .mph file for reference.

Gradient Coil Simulation Data Analysis

Fri, 28 Mar 2025 17:04:53 +0000

I have a simulation where I have created a set of gradient coils along the z-axis, and the coils are to create a gradient magnetic field across the field of view. I want to plot the change in magnetic flux density, and the magnetic field norm in the x-direction, as a function of distance along the z-axis. I have tried both line integration, and plotting of the magnetic flux density w.r.t. distance, and both return 0, which is satisfactory for line integration, but the magnetic flux density always returns a 0 line value. This is incorrect because when I place points along the z-axis, measuring the magnetic flux density in the x-direction, I receive values. How would I go about creating a 1-D plot that gives me the magnetic flux density as a function of distance?

Evaluate how much charge deposited on wall

Fri, 28 Mar 2025 02:30:50 +0000

Hello I see in the documentation of Comsol 6.3 that

This text bellow is from documentation Electron Density Wall Settings The Use wall for electron density checkbox is selected by default and deactivates the wall boundary condition for the electron density only. When activated, the flux of electrons lost to the wall is self-consistently computed based on the electron temperature and electron density adjacent to the wall. When this is deactivated, the flux term is no longer applied for the corresponding equation, and the software will impose a zero flux boundary condition for that dependent variable. Select or enter a value for the Secondary emission flux (SI unit: 1/m2·s)). The secondary emission flux can either be a user-defined expression or come from another physics interface, typically The Heavy Species Transport Interface. When using this boundary condition in the Plasma interface, this setting is not editable, and the secondary emission flux comes from the surface reactions.

My question is this Please tell me if this understanding is correct. Suppose i do a electrostatic discharge and give wall boundary condition i can estimate how much electrons got deposited on the wall? If so how can i visualize it?

Defining port in 2D (not axisymmetric) with microwaves launched from open-ended circular waveguide

Thu, 27 Mar 2025 20:28:33 +0000

Hi COMSOL folks,

Here is my problem: I am trying to focus microwaves launched from an open-ended circular waveguide at 2.45 GHz to the inside of a small diverging nozzle (essentially a cone with metal walls). I am working inside of a vacuum chamber, so I am limited in space, therefore I am trying to implement a custom off-axis parabolic (OAP) mirror, assuming a mostly collimated incident beam which will be reflected and focused to a point, at which location I'd place the nozzle.

In the image attached, I have the OAP, section of circular waveguide, and the conical nozzle modeled in the domain. I would like to put this setup into a 2D simulation, but I am having a hard time knowing how to correctly define the port as there is no circular port option. I don't have the strongest EM background, but I believe I would need the e-field distribution for TE11 mode with a circular waveguide.

I would really appreciate any advice, comments, or guidance on how to approach this simulation or improve the setup as a whole. Thanks in advance!

Transient Induction Heating Simulation with BH Curve Material

Thu, 27 Mar 2025 15:18:58 +0000

Dear respected experts,

I am encountering difficulties in simulating an induction heating process and would appreciate your guidance. The configuration consists of a U-shaped single-turn coil carrying AC current to heat a soft magnetic material plate via induction.

The key technical issue is: When defining the soft magnetic material's constitutive relation using a BH curve, frequency-domain analysis becomes unavailable. However, when switching to transient calculation, I cannot apply current excitation to the single-turn coil. Attempting voltage excitation results in severe convergence issues.

Core requirements:

Maintain BH curve definition (necessary for accurate hysteresis loss calculation)

Perform transient analysis with AC current in U-shaped coil

Constraints observed:

Relative permeability models cannot replicate measured hysteresis losses

Experimental data confirms hysteresis loss dominates over eddy current loss in this heating process

Has anyone successfully implemented such a combination (BH curve material + transient AC current excitation for single-turn coils)? Any suggestions on solver settings/excitation implementation would be highly valuable.

singular matrix error

Thu, 27 Mar 2025 10:48:00 +0000

Hello, I am having trouble computing the multiphysics of multiphase flow in porous media for a water electrolyzer. The geometry is configured as a half-cell. Upon reviewing the conditions, I didn't find any issues, but I keep encountering a singular matrix error.

Could this be due to boundary conditions or solver settings?

I would appreciate any help.

Mobility plot of example: mobility model

Thu, 27 Mar 2025 07:38:34 +0000

Hello, I tried to do semiconductor MOSFET mobility model. When I plotted 2D mobility value, I found there is specific depth layer where mobility decreased. I can't understand why there is situation occurs.