Quantcast
Channel: Latest Discussions - COMSOL Forums
Viewing all 26323 articles
Browse latest View live

Different solvers and their influence on the obtained results

$
0
0

Hello!

I'm currently using COMSOL for a 3D microfluidics problem and I am performing a mesh refinement study. As I am quite new on the subject, I am only using different element sizes for a default mesh, and going through some of the sizes. So I have some questions about this:

1) Are these default meshes good enough for a two phase flow problem? I've seen in literature that quadrilateral (cubic in 3D) are use in such studies if the mesh is mentioned in the first place.

2) I started by using the 'Coarse' mesh and went all the to the 'Finer' one. I notived that COMSOL used PARDISO to solve the problem when the mesh was 'Coarse', 'Normal' and 'Fine', and the results were getting closer, as expected. However, GMRES was used with the 'Finer' mesh, and the results were a bit different. I was monitoring the fluid velocity at different points and in two plots, the velocity was quite higher than the two preivous ones. So I was wondering if the solver has any effect on the result, and how can I make COMSOL use PARDISO so I can be coherent all the way.


Undefined variable , Failed to evaluate Jacobian of expression

$
0
0
  • Feature: Stationary Solver 1 (sol1/s1) Undefined variable.
  • Variable: comp1.w
  • Geometry: geom1
  • Domain: 0 Failed to evaluate Jacobian of expression.
  • Expression: comp1.w Failed to evaluate Jacobian of operator.
  • Operator: mean
  • Geometry: geom1
  • Boundaries: 183-186 Failed to evaluate Jacobian of variable.
  • Variable: comp1.w
  • Geometry: geom1
  • Boundaries: 183-186 Failed to evaluate Jacobian of expression.
  • Expression: d((comp1.tffs.hcomp1.tffs.rho(comp1.tffs.vavextest(comp1.pfilmTx)+comp1.tffs.vaveytest(comp1.pfilmTy)+comp1.tffs.vaveztest(comp1.pfilmTz)))(dvol_spatial),{test@5}) Failed to evaluate Jacobian of expression.
  • Expression: (comp1.tffs.hcomp1.tffs.rho(comp1.tffs.vavextest(comp1.pfilmTx)+comp1.tffs.vaveytest(comp1.pfilmTy)+comp1.tffs.vaveztest(comp1.pfilmTz)))(dvol_spatial)

Question about Geometrical Optics

$
0
0

Hi all,

I'm new to COMSOL MultyPhysics, and I'm writing this email to ask you some questions about my Geometrical Optics simulation.

Please refer to the mph file I've attached.


  1. When I try to use help function, the error below shows up.

*Failed to establish a connection to the documentation server:

"connect timed out"

Try the online help preference setting*

How can I solve this?


  1. Before I start, I'm avoiding using "2D Axisymmetry" geometry since the real simulation I'm trying to run is not axisymmetric. To give you a brief concept of my simple simulation, the light with a uniform densityparallel to the +z-axis from the larger surface of the cone passes through the material with a certain attenuation coefficient (or imaginary refractive index) and the intensity is detected at the lower surface of the flat cylinder. I'm trying to get a "light intensity map" from "3D Plot Group 2".

To accomplish this purpose, I'm using an accumulator with these settings Accumulator Type: Density Accumulation over: Rays in boundary elements Accumulated variable name: rpb Source: gop.Q [W] Units: Custom unit [W/m^2] Smoothing radius : 0.02 [m]

However, the tutorial I found ( https://www.comsol.com/video/intro-to-modeling-ray-tracing-with-comsol-multiphysics-nov-22-2018 ) use the following settings

Accumulator Type: Density Accumulation over: Rays in boundary elements Accumulated variable name: rpb Source: gop.logI [1] Units: Custom unit [1/m^2] Smoothing radius : 0 [m]

Is the tutorial correct? In my instinct, it seems like I have to put total power ( gop.Q ) for the source so that the accumulator can calculate the "density" of "power", which is intensity. However, the tutorial is calculating the "density" of the (log of ) intensity over the area [m^2] which is definitely not the intensity itself.

So is my settings appropriate to get the intensity? or the tutorial is right?


  1. For the same purpose, I'm using the 'Density' instead of the default 'projected plane grid' for 'Inlet 1 Settings'. Would this help to generate the light with uniform density?

  2. The result '3D Plot Group - Surface 1' is very 'grainy' and doesn't seem to represent the real result (which should be smooth). Other than making the mesh finer and increasing the number of rays, is there something else I can do to make my simulation actually representing the real physical result?

  3. How can I check the "automatically assigned variables"? Such as gop.I, gop.logI, gop.I0, and gop.Q? I had to find each of them at several tutorials and it would be easier for me if I know how to check them up myself.


  1. From other tutorials regarding thermodynamics, 'T' is used for an automatically assigned variable for the temperature at each point. Is this('T' for temperature) different from the variables mentioned above (gop.I and so on)? The reason for asking this is that when I put in 'gop.I' in the Expression for Surface 1 of 3D Plot Group 2, it gives me an error below.

Cannot Evaluate Expression. Undefined variable. - Variable: comp1.gop.I - Geometry: geom1 - Boundary: 9 Failed to evaluate expression. - Expression: comp1.gop.I - Plot: surf1 (Surface)

Also, when using Matlab Livelink function,

*mphinterp(model,'T','coord',[0;0;0]) *seems to give the temperature at [0;0;0] point, while

mphinterp(model,'gop.I','coord',[0;0;0]) gives me an error below.

*Error using mphinterp mphinterp (0): Java exception occurred: Exception:

com.comsol.util.exceptions.FlException: Failed to evaluate expression

Messages:

Cannot evaluate expression.

Undefined variable.

  • Variable: comp1.gop.I

  • Geometry: geom1

  • Domain: 1

Failed to evaluate expression.

  • Expression: comp1.gop.I*

How can I solve this? And where can I find these (potentially different types of ) automatically assigned variables ('T' for temperature, 'V' for voltage, and so on...)?


  1. I would like to use the 'mphinterp()' function for the variable 'rpb' I mentioned above. I've tried both 'rpb' and 'gop.wall1.bacc1.rpb', but both doesn't work... What should I put in instead?

  1. I would like to get an "intensity map" I've made in "3D Plot group 2" as a '2D Plot group' format. How can I do this?

9.I'm currently using 'Compute intensity and power' option for 'Geometrical Optics Settings'. Then, I can put in the settings for 'Inlet 1', both the initial light intensity (I_0) and the total source power (P_src) which are definitely correlated. Is the initial light intensity automatically ignored when I put both of them?


  1. From '3D Plot group 2 - Surface 1', I would like to get a (numeric) table of the intensity for each point. Maybe with X=1mm and Y=1mm steps. How can I do this?

Thank you for reading these long questions!

Best regards, Rovert(Hyun Sung) Cho

Simulating X-ray

$
0
0

Hi all, I am new to comsol. Is it possible to simulate X-ray using COMSOL Multiphysics.

Thanks.

"Undefined" variable won't allow my study to execute

$
0
0

Hi all,

First off this this is my first real project in comsol. It was handed to me to fix, and I've gotten it this far. Now my defined variable shows up as undefined, causing the entire program to halt. If I could get any help I'd appreciate it. Thank you!

Non-Newtonian Laminar Flow, user defined flow law

$
0
0

I am trying to model ice (non-Newtonian laminar) flow under martian conditions using a flow law used in the literature for non-terrestrial ice modeling. It is essentially a power law form but combining multiple modes of deformation that use different exponents in the power law. This does not fit into any of the pre-defined models available and it appears that in version 5.5 it is no longer possible to use user-defined functions for flow laws. Am I missing something or do I need to go back to 5.4 to do this?

variable names

$
0
0

I would like to define functions for some material properties and need to use variables like shear rate in the function. I found from other posts that, for example, the laminar flow shear rate is spf.sr, but I was wondering if there's a way to find variable names in general either through the interface or documentation?

general stress-strain relation for fluid

$
0
0

Is it possible to define a general stress-strain relation for non-Newtonian laminar flow? Laminar flow seems to require you to define a dynamic viscosity. However I am trying to implement a flow law with a complicated stress-strain relation that does not yield a simple analytical equation for viscosity (see attached). Is it possible to define fluid properties more generally (other than externally solving for viscosity as a function of strain rate and then inputting an interpolation function for the viscosity)?


About Langmuir probe in comsol

$
0
0
Hello,everyone!
I have a question in my thesis topic.I wonder if the Langmuir probe can work in the established comsol model, and if so, is it possible to define the probe portion directly in the established model?
I'm so confused and try to find some useful articles,but not.Please reply me if you have any ideas.Thank you very much.

Explanation for the irregularity in my model

$
0
0

Hello, I made a 2D model for a single slit diffraction using EM waves and frequency domain. There seems to be some irregularity in the planar wave before the slit. Is this due to some reflection that I am not aware of or is it interference from the many source points along the bottom boundary?

https://www.youtube.com/watch?v=XnT45kzu53c This is the video I was basing my work off of and they seem to have a similar problem.

Simulating a half filled eppendorf tube in a centrifuge

$
0
0

Dear all,

We try to model the flows in an eppendorf tube while it's rotating in a centrifuge (positioned under an angle).

We where succesfull in simulating a fully filled tube using the rotating machinary module in 3D. Now we want to expand the simulations to also include half filled tubes (water at the bottom and air at the top). Unfortunatly we find that the rotating machinary module only offers mixture models for multiphase flow.

Has anybody a good suggestion on which modules or combination of modules is best used to tackle this kind of problem ?

With Kind Regards Sander Stroobants and colleagues.

Change the label of reference system

$
0
0

Hi, is it possible to change the default label of the reference system {x, y, z} in a different triplet named {1, 2, 3}? It is just a matter of graphics especially when exporting an image of the model. Thanks in advanced, Regards.

Additive Manufacturing model - requesting feedback

$
0
0

Hello fellow Comsol users,

I am currently working on a layer-by-layer thermal simulation for additive manufacturing and by now, my results are not bad, though still far from perfect. In order to further improve the current model, I'd love to hear your feedback and maybe a few tips.

The meshing is a topic I'm most uncertain about right now, since I feel there's the highest potential for optimization.

Feel free to take a look. Also attached is experimental data as well as a version of my current simulation model. I solved a problem earlier where the first layers temperature spike was much too high. As you can see, I still have a slight problem with high temperatures. Any ideas why?

solid.pw

$
0
0

Hi all,

I have some doubts regarding the variable solid.pw ; In particular, which is its meaning and which is the difference between solid.pw and solid.pm?

Thanks a lot, Maria

2DEG In COMSOL

$
0
0

If anyone could share their insight on acheiving the 2DEG in a HEMT in COMSOL, I would be grateful.


How to apply comsol moving mesh in weak form pdes application mode?

$
0
0

Hi,

I am simulating a bubble moving in fluid using Weak form PDEs application mode. I am applying moving mesh capability. The problem is, my mesh moves but the bubble do not move. Please can someone upload an mph file where moving mesh capability is utalized with weak from PDEs?. I cannot upload my file because the space is not enough.

Mobility analysis in electron beam interaction

$
0
0

Hi everyone, I am trying to simulate a model to calculate the mobility of the charge carriers in a semiconductor on which some potential is applied and it interacts with a beam of incoming electrons and photons. How to proceed for this? I am new to COMSOL. Any leads will be highly appreciated.

How to set the Values of dependent variables

$
0
0

I want to solve a pair of equations which are based on other two study steps, but in the setting of the "values of variables not solved for", I can only choose on study( and i actuallly want to add two study). How do I pass the solved value of the varibles into the equations that I try to solve subsequently? For example, I first solve a pair of question to get the Value of U and V, then I solve another equation to get the value of p, the finally I want to solve w and k, whose equations involve with U, V and p. How should I pass the value of U, V and p.

Use interior fan to resemble submersible mixer

$
0
0

I am trying to make a 2D model of a reservoir filled with liquid, with submersible mixers in it. In 2D, the only option I see is to use the interior fan, but I am uncertain what values I should use as input for the static pressure at no flow and the free delivery flow rate. I know the stationary power of the motor, the stationary RPM and the diameter of the propeller. Can anyone advise me on how I would translate these known values to the ones asked for the interior fan?

Thanks in advance!

How to find the temperature of meniscus during capillary filling

$
0
0

Hello,            I am trying to simulate 3d rectangular capillary filling using Leve set and two-phase flow, but simultaneously heat or (put at T+Tref) the one wall of the capillary, everything including (meniscus position(ls.Vf1), velocity, etc.) I found easily. The problem is, I am not able to calculate/plot the temperature profile along with the meniscus at some instant (look at attached figure). Thanks in advance!

Viewing all 26323 articles
Browse latest View live