Latest Discussions - COMSOL Forums

How to replace a conical "Release from Grid" ray with a Gaussian Beam?

Wed, 18 Jun 2025 14:19:20 +0000

I have a simulation of an off-axis parabolic mirror that colimates a diverging beam from a point source. I have designed. (added below)

At the moment I am doing ray tracing from a "Release from grid" point with a conical shaped ray and a specific cone angle. It works great, but I would like to replace this ray with a gaussian beam. How can I do this in comsol?

I have been looking at the documentation "gaussian beam", but it seems I cannot find a way to get it to work in my comsol... https://doc.comsol.com/6.2/doc/com.comsol.help.roptics/roptics_ug_optics.6.23.html#:~:text=Use%20the%20Gaussian%20Beam%20feature,Value%20of%20Auxiliary%20Dependent%20Variables.&text=Enter%20the%20Beam%20waist%20position,positive%20x%20direction%20in%202D.

Best Approach for Modeling Strain in a Thin Lamella Due to Differential Thermal Expansion

Wed, 18 Jun 2025 14:06:56 +0000

I'm looking for the optimal approach to model the strain induced in a thin lamella on a substrate caused by differential thermal expansion. The system undergoes a slow temperature change from 300K to sub-Kelvin levels. Are there specific modules, boundary conditions or studies in COMSOL that would be best suited for this scenario? Any insights on material properties or thermal stress analysis techniques would be greatly appreciated.

Thank you in advance for your guidance!

Li-ion batteries aging model behavior's issues

Wed, 18 Jun 2025 08:46:08 +0000

Hello, I have developed an electrochemical aging model for an LFP cell (with graphite on the negative electrode). The assumptions on which the model is based are: 1. Aging is solely due to SEI growth and lithium plating at the negative electrode. 2. Lithium plating occurs only when the overpotential at the negative electrode drops below 0 V. However, I am encountering three issues with the model's behavior: 1. The cell capacity drops more quickly at lower C-rates, whereas the drop should slow down when the current rate decreases. 2. During cycling, when the overpotential at the negative electrode falls below 0 V, the capacity fade slows down instead of accelerating when cycling at et voltage range between 3.2V and 3.65V. 3. Temperature variation does not appear to influence the rate of capacity degradation.

I kindly ask for your support in resolving this issue, which has been a significant obstacle to my progress for several weeks. To reduce the file size, I removed the stored simulation results. Please remember to run the simulation after opening the file.

Request for Trial code

Wed, 18 Jun 2025 07:02:31 +0000

I am working on a multiphysics model involving fluid flow and heat transfer as part of my B.Tech project

Electron Heating in COMSOL

Tue, 17 Jun 2025 10:19:01 +0000

I am trying to simulate heating through electron particle impingment in COMSOL. Can anyone help me? Ideally the electron particles are 1000, they propagate through a vacuum and collide with a tungsten block placed at an angle to direct xrays. I am open to suggestions

Setting up Parametric Eigenfrequency Sweep with Piecewise Functions

Mon, 16 Jun 2025 21:03:13 +0000

Hi!

I'm trying to run an eigenfrequency study of a lattice (a GKMG sweep). To do this, I want to parametric sweep over the parameter e. To define the path that e follows, I created piecewise parameters kx and ky, which are defined as functions of e.

How can I define e such that the parametric sweep follows the path set out by kx and ky?

Thank you so much! I've been stuck on this for a while so I'd really appreciate any help.

Query regarding Multiscale Modelling in COMSOL

Mon, 16 Jun 2025 12:58:37 +0000

Hello,

I have looked into the Multiscale Model example provided by COMSOL "https://www.comsol.com/model/micromechanics-of-failure-multiscale-analysis-of-a-composite-structure-107951". The following model allows a 1-way communication where homogenized properties are transferred to the Macro-scale structure. Is it possible to expand this implementation to allow a 2-way communication where strains from macro-scale are transferred to the RVE to compute the stresses and upscaled to the Macro-structure?

Moving Mesh Simulation – Fatal Error During Assembly/Calculation

Fri, 13 Jun 2025 07:39:29 +0000

Dear colleagues,

I would like to seek your advice on a moving mesh simulation issue.

Current Model: I am simulating "the motion of a current-carrying conductor rod under the influence of a coil's magnetic field."

Physics Interfaces Used: Moving Mesh Magnetic Fields (AC/DC Module) Global ODEs and DAEs (for equation-based modeling)

Implementation Details: In the Magnetic Fields interface, I added a Force Calculation condition to compute the electromagnetic (Lorentz) force. The calculated force is passed to a global differential equation: F=ma=mxtt The displacement is then fed into the Moving Mesh interface to update the conductor rod's position. Since the electromagnetic force changes as the rod moves, I use Transient Analysis in the study step.

Issue: Despite multiple adjustments, the simulation consistently fails with the error: "A fatal error occurred during assembly or calculation."

Request: Has anyone encountered similar issues or have suggestions for troubleshooting? Any insights on: Correct coupling between Magnetic Fields, ODE, and Moving Mesh? Solver settings for such coupled electromechanical problems? Common pitfalls in force-driven moving mesh simulations? Thank you in advance for your expertise! Best regards,

Liu

Setting of anisotropic conductivity of carbon fiber braided composites

Fri, 13 Jun 2025 05:03:12 +0000

I studied the cases in the official case library, trying to define the anisotropic conductivity of carbon fiber with curvilinear coordinates（cc）, but I don't know how to correctly quote the curve coordinates in the current (ec) and magnetic field (mf), do you have a partner who has done similar research?!!!

Modelling of B-Field with AC/DC Module

Thu, 12 Jun 2025 21:48:35 +0000

I am currently working on building a model to simulate the B-Field of an electro-magnet. This is my first project on COMSOL and I am self taught so I am know i am making errors in my simulation approach for sure. I've attached an outline of my approach showing my Boundary conditions, Mesh and results. The results are not too far of the measured data (in a small location - I still need to check the rest of the field area) but i don't think my model is very stable.

RLC circuit not giving correct resonance

Thu, 12 Jun 2025 21:05:29 +0000

Hi, I am trying to test a circuit with one lumped port and 3 lumped elements, one for R, L and C respectively. I am using 1 mohm, 1 nF and 1 nH as my values, so I should get a resonance at 159.15 MHz, with a Q of 1000 giving a width of 0.2 MHz. But when I sweep the frequency around this, I do not see any resonance. Why might this be? Thanks.

Possible to Obtain Current/Current Density in Electrical Currents, EC, from different Physics?

Thu, 12 Jun 2025 18:15:05 +0000

Hello,

Is it possible to impose a current/current density in a Electrical Current domain, from a different physics? In my model, i have a coil, hts superconducting coil, using Magnetic Field Formulation, mfh, where i have a current density in the hts coil domain. How can i have the same current density, seen in the mfh domain, in the Eletrical Currents side? How can i impose it/connect it?

My goal is to then have a terminal so i can connect it (hts coil) to an Electrical Circuit. Is there any other/possible way to connect this superconducting coil to an electrical circuit?

Thanks

Incorrect S11 resonant frequency

Thu, 12 Jun 2025 16:24:23 +0000

So based on previous discussion, I have managed to obtain better S11 values in negative range. However, I am unable to get the correct resonant frequency.

I am not sure if it because of improper modelling, or that of wrong lumped port definition.

Based on the image below https://gusbertianalog.com/content/images/2024/11/MEAS_ROUND_CC.png

from the author of the article https://gusbertianalog.com/loop-antenna-design/

I am supposed to get the resonant frequency at 434 MHz, but instead i get it at 460 MHz.

I have uploaded the comsol file.

(P.S. I tried to upload this question, but i got a message saying i cant view this question in the forum, so i am uploading the question again. )

Activated carbon fiber packed bed adsorption

Thu, 12 Jun 2025 08:39:32 +0000

How to couple chemical adsorption and physical adsorption when both exist?

Model amorphous silicon crystallization with COMSOL

Wed, 11 Jun 2025 20:28:13 +0000

Hello, I am interested in modeling laser induced crystallization of amorphous silicon (a-Si). The chain of phase transformations I am interested in is: a-Si -> Liquid Si -> polycrystalline Si. Is there a good way to model this in COMSOL? The main difficulty I have with Phase Change Material subnode of Heat Transfer in Solids module is that upon cooling the material drops back into a-Si phase which is not what should happen.

Thank you!

CFD - Phase Field in Fluids - Solver Settings

Wed, 11 Jun 2025 14:48:39 +0000

Hey Comsol Peeps,

I struggle with setting up a CFD simulation in which a small amount of liquid (micro liter regime) should be rotated within a hollow shell of 2mm radius. I always run into convergence errors after less than 1000 time steps no matter my settings.

The simulation is 2D and so far I only try to simulate the fluid resting inside the shell (no rotation). The used physics are the lamina flow interface with gravity included (the gravity will ramp up over time to make the initiation less demanding) and the Two-Phase Flow, Phase field interface. As phase field settings I chose for the interface thickness the default: pf.ep_default as mobility tuning parameter 0.5 m*s/kg as external free energy 0 J/m^3 (these values are rather guesses)

When computing I run into the error "Time adaptation interrupted." This happens at the time 0.057244s Looking at the velocity surface plot one can see two velocity hotspots at the fluid interfaces. However, checking the mesh at these hotspot , the quality doesn't seem to be too bad here. (see pictures) Does someone have an idea how to fix my simulation?

Fluid solution mapping error in Fluid-Acoustic Simulation

Wed, 11 Jun 2025 14:40:13 +0000

Hello everyone,

I am currently doing an Acoustic-fluid simulation with stationnary turbulent flow, convected wave equation and background fluid flow coupling modules. I have only one geometry for the whole study, in the first place i do the fluid simulation, then I map the fluid fields onto my acoustic mesh and finally I do my acoustic simulation with the background flow set up as a variable not solved for. My problem is that i already did many simulation using a simple geometry for my conduct (3D cylinder) which worked totally fine.

Now that I am using a more complex one based on CT-scans, I face issues during the mapping of my fluid fields onto the acoustic mesh. The error message is the following :

" Extrusion coupling operator interpolation failed for at least one point: no source point found. - Geometry: geom1 - Domain: 1 - Destination coordinates: [-0.00146266445760077; -0.0021864553895055991; 0.00060947754522279061], [-0.0014856891031490996; -0.0021811204460443394; 0.0006247936128450256], [-0.001450351221716216; -0.002194843619005469; 0.0006247936128450256], [-0.0014270113859349211; -0.0021981395225136669; 0.00060828294733559498], [-0.0014966022294479824; -0.00217111449601754; 0.00060828294733559509], ... Failed to evaluate expression. - Expression: withsol('sol1',comp1.spf.pA) Failed to evaluate expression. - Expression: (comp1.bffc1.p_map-withsol('sol1',comp1.spf.pA))dvol_spatial Failed to evaluate Jacobian of expression. - Expression: ((comp1.bffc1.p_map-withsol('sol1',comp1.spf.pA))test(comp1.bffc1.p_map))*(dvol_spatial) "

I understand that COMSOL struggles to find source point in my fluid domain but i can't understand why since my CFD simulation works fine and i have good results. Even with a mesh refinement of the source domain i keep having this problem, do you have any idea why and how I could solve this ?

I sincerely thank you for your time and your answers !

Acoustic streaming due to bubble oscillations

Tue, 10 Jun 2025 21:23:04 +0000

Hi,

I'm trying to simulate the acoustic streaming due to the bubble oscillations in a cylindrical cavity under acoustic fields from literature. I'm unable to get the expected streaming. It would be of great help if someone could help me figure out what is wrong in my setup. I expected the flow away from the bubble as in the reference but mine is reverse.

Thanks,

Learning Center Simulation: Modeling a Coil Connected to a Circuit

Mon, 09 Jun 2025 23:36:29 +0000

Hello,

I am following along to the COMSOL Learning Center video titled Modeling of a Coil Connected to a Circuit and have come across an issue. I am able to follow along seamlessly with the tutorial until 41:25, where the instructor changes the core to instead be a hollow solid instead of just being off to the side. When it's done in the video, COMSOL immediately understands what all the domains should be assigned to and computes the solution. But in my file (COMSOL 6.2.0.339) domain 5, which is the air inside the hollow core, can not have any domains applied to it as it is "Not applicable." Magnetic insulation is applied to both the coils and the core by default and I get a trivial solution. Is there something that I'm doing wrong when trying to follow along?

I attached both a .mph file from just before the hollow core is added, and a .mph file with the trivial solution (this one has Issue appended to the title).

surface impedance and sound absorption coefficient

Mon, 09 Jun 2025 11:51:30 +0000

How to calculate surface impedance and sound absorption coefficient for a two cavity parallel MPP system.