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

Duplicate variable name

$
0
0
Hello all,

I got some duplicate variable name error when computing my model(see the attachment), which including the piezoelectric device, laminar flow and deformed geometry. Should I use the physics such as FSI to get less variables ? And I notice that the duplicate variables are all concerning the coordinates.

BTW, when I built the model use laminar flow Physics in version 3.5a, I found there is normally the dependent variable mod1.u_lm and mod1.v_lm (N/m2) in Study , but when I try to build the same model in 4.2a, there are only dependent variables pre-defined by Physics( such as mod1.u for velocity). Is there any chance we can modify the dependent variables based on our needs?

Thanks for any help!

Best regards.

co-efficient of thermal expansion

$
0
0
Hello evrybody,

Can anyone help me to know if I can input my co-efficient of thermal expansion as a function of Stress and Temperature or not. And if possible how to do that.

Thank you.

Output pressure of the transducer adding Vac + Vdc

$
0
0
hello,
I 'm working on a capacitive transducer, using the electromechanics(emi) model. Now I want to observe the output pressure of the transducer in the frequency domain and I add a voltage V=Vac+50 volt.

First, in the Model 1—Definitions,I define Waveform and Analytic Functions,and plot the figure well. The function Analytic's name is an2 ,and Expression is wv2(x)+50.
Then I add the fixed voltage on the terminal of the top plate,I wrote an2(x) in the Electric potential V0.
Last, run it, there is an error which tells me that "constraint found for variables in diferent Segregated groups.Try to merge these groups. —Feature: Stationary Solver(sol1/s1) —constraint found for variables in diferent Segregated groups.Try to merge these groups"

So I beg your suggestions and introductions:
1. Am I right to add the Vac like this?
2. How can I solve the error?
3. Could you give me some directions about how to solve the output pressure please?
Thanks for any help!

Best regards

error in statonary analysis of a biased resonator

$
0
0
Dear comsol community,
I'm working on comsol4.2a , I would be glad if i can get a solution for my following two queries:

1) I have gone through the modelling instructions given in the 'Stationary Analysis of a Biased Resonator' doc. and after computation i found a following error which is difficult to interpret ,
'FAILED TO FIND A SOLUTION,
MAXIMUM nUMBER OF nEWTON ITERATIONS REACHEAD.
RETURN SOLUTION IS NOT CONVERGED.
-feature: stationary solver1(sol1/s1)
-error: failed to find a solution.


2) Is it possible to excite the cantilever beam and find the frequency modes and plot the graph of frequency vs. DC Bias
as in mems module doc. it is the fixed-fixed beam which is considered.

plot

$
0
0
Hi
I simulate the mode profile of the rib waveguide.
Now i want to plot the normalised intensity(arb unit) in the x and y directions.
Do you know how can i do this?what is the arb unit?
how can i find the normalised intensity?

Thanks in advance

Resonance frequency of beam (calculation+plot)

$
0
0
HI

I want to calculate and plot resonance frequency of cantilever that I designed in MEMS module. Based on the previous discussions here I understand that I should do the study in "Eigenfrequency" and then i can get different modes of resonance frequency, is that right?
I did what I have mentioned above, and I got 6 eigenfrequency. However there are two problems.

1- I can not understand what equation comsol is using for calculating this resonance frequency. For mu theoretical calculation of first mode of resonance frequency I used this formula "f=0.162(t/L^2)squrt(E/Rho)
by using this equation I get resonance frequency of 95 K, but the simulations give me first eigenfrequency of 117K. Is this differences normal?

2-I can not plot find out how I can plot frequency vs. amplitude to have a simple resonance frequency plot with a pick on resonance frequency.

I would be grateful if someone can help me with this.

Regards
Bahar

Out of memory during LU factorization

$
0
0
Hi there,

I am a new member in the COMSOL community. I am trying to solve a simple 3D problem having diffusion and chemical reaction within the subdomain. However, I am getting this message "out of memory during LU factorization." Is my computer running out of memory to solve this particular problem or should I use a better solver to solve this problem? Could anyone share his/her experiences with COMSOL?

Many thanks,

Mohammed

Piezo Loading

$
0
0
I have a 3D piezo model and the output is currently a floating potential. This tells me the theoretical unloaded voltage gain. I'm working on possibly using a global equation of a simple ohms law relation between voltage resistance and power. Since the resistive load can be purely numerical. But I haven't been able to figure out how to make this work. In the actual simulation. Since the voltage output is a boundary condition and the relationship need to be global.

Any feedback would be greatly appreciated, or any links to pages with similar information.

Placing a finite element model into an electrical circuit in COMSOL 4.1

$
0
0
Hello,

I have recently made the switch from COMSOL 3.5a to COMSOL 4.1, and I am currently unable to do things that I used to be able to do very easily, namely placing my 3-D model (drawn and defined in COMSOL) into an electrical circuit simulation. I easily achieved this in COMSOL 3.5a using the “Conductive Media DC (emdc)” physics module which includes the SPICE Circuit Editor. In the circuit editor, I defined my circuit with the following code:

V1 0 1 9
R1 1 2 1000
C1 2 0 0.00000000003
X1 2 0 resistor
.SUBCKT resistor a b COMSOL: *
.ENDS

This corresponds to a voltage source (9V) across nodes 0 and 1, a 1k resistor between nodes 1 and 2, and the parallel connection of a 30 pF capacitor and my 3-D silicon nanowire model (X1) across nodes 2 and 0. Nodes 2 and 0 correspond with circuit terminals a and b, respectively, which are defined in the Boundary Settings of Conductive Media DC (emdc). This was very easy and convenient to set up.

I have been trying to do this same thing with COMSOL 4.1, trying to learn this from one of the models in the AC/DC Module Model Gallery on the website called “Inductor in an Amplifier Circuit”:
www.comsol.com/showroom/documentation/model/990/

This tutorial could not teach me how to place a 3-D model into the Electric Circuit (cir) module; the explanation in given in the PDF document is quite unclear to me. Their SPICE Netlist appears as follows:

* BJT Amplifier circuit
.OPTIONS TNOM=27
.TEMP 27
Vin 1 0 sin(0 1 10kHz)
Vcc 4 0 15
Rg 1 2 100
Cin 2 3 10u
R1 4 3 47k
R2 3 0 10k
X1 4 5 inductor
RE 7 0 1k
Cout 5 6 10u
Rl 6 0 10k
Q1 5 3 7 BJT
.MODEL BJT NPN(Is=15f Ise=15f Isc=0 Bf=260 Br=6.1
+ Ikf=.3 Xtb=1.5 Ne=1.3 Nc=2 Rc=1 Rb=10 Eg=1.11
+ Cjc=7.5p Mjc=.35 Vjc=.75 Fc=.5 Cje=20p Mje=0.4 Vje=0.75
+ Vaf=75 Xtf=3 Xti=3)
.SUBCKT inductor V_coil I_coil COMSOL: *
.ENDS
.END

First off, I don’t see any way to create/edit the SPICE netlist in this version of COMSOL, the only option is to import and it is asking for a file with a .cir extension. Can I write the SPICE netlist in a .txt document and import that? Or do I make a .txt file and then change file type to .cir and then import that?

Second, why are there two lines “.ENDS” and “.END” here and in version 3.5a I used only one line “.ENDS”? Is this because of the BJT model? If I only have basic parts and one 3-D model (V1, R1, C1, X1) do I use both lines .ENDS and .END? Or just .ENDS?

Third, the line of code which defines X1 as the 3-D model now has variable names "V_coil" and "I_coil" instead of terminal names "a" and "b". The PDF document states:

“The variable V_coil must give the voltage over the device, defined in the global scope. I_coil must be a global variable used in the model as a current through the device.”

However, there are not even any global variables defined in the simulation file provided, there are just global parameters, none of which are V_coil or I_coil. I could not find one instance of the text “V_coil” or “I_coil” anywhere in this simulation file. Also they say V_coil is defined in the global scope. What is the global scope? How can I view the global scope? They don’t mean the Global Definitions tab do they? Because there are no variables defined in this file’s Global Definitions.

Fourth, I cannot tell from the boundary conditions which boundaries are corresponding to the terminals of the device which will be connected to nodes in the circuit. It is not clear which edges/surfaces are corresponding to the two terminals of the device. For example, I have a 3-D rectangular prism and want to make two opposite faces be the two terminals of the resistor. What kind of boundary conditions do I use for these faces with the physics module “Joule heating”? Also, how would I define variable V_wire in the “global scope” and I_wire as a global variable for this simple resistor model?

Fifth, it seems that their 3-D model appears as a component “External I Vs. U1” in the Electrical Circuit (cir) section of the model builder, as its node names are 4 and 5. Is it possible to define a circuit including a 3-D model using just the model builder (i.e. no SPICE netlist) since all components seem to be available in the Electric Circuit (cir) model builder?


I used to implement my 3-D models in the SPICE circuit so easily in COMSOL 3.5a, but to do the same exact thing in COMSOL 4.1 has been extremely difficult for me and I cannot figure it out even after reading the tutorial and looking at the example simulation file. I would very much appreciate some help.

Thank you,
Adam

Mark my experimental values to the curve point in the plot window of COMSOL

$
0
0
Hi Everybody

The purpose of my question is to compare the theoretical curve obtained on COMSOL with my experimental values on the plot window.

I have a curve point on the plot window. I have a few experimental values and I'd like to mark them to the curve point precisely.

Is this possible and how? if not is there another solution to my problem.
Thank you.

Error message "Cannot import a 2D geometry into a 3D sequence." in RE

$
0
0
Hi!

I have a 2D Reaction Engineering modell and want to add some more physics through the Synchronize with Comsol Multiphysics. I have imported my 2D geometry to the RE model.

When I then import the same geometry to the Model 2 that Comsol adds I get the error message above even though there is no 3D models or geometries involved.


Anyone have a clue whats going on ?


regards
Nils Eivind

How to Convert a 2D model Into 3d model

$
0
0

I want to know that how i can convert a 2D model into 3D model for simulation.
--
Vineet Bansal
Project Associate, NPMASS
Electronic Science Department
Kurukshetra University, Haryana

3D Poisson equation problem coupled with Nernst-Planck equation

$
0
0
I used coupled Electrostatics physics and Transport of Diluted Species physics to simulate this problem.
The source term in Poisson equation is the difference between cations and anions.
Usually, the difference should be zero due to the electroneutrality. In some conditions, it would be not.

My problem is that, when I used 2D model, the simulation converged even with a coarse mesh. When it comes to 3D situation, it no longer converged, no matter how fine the mesh is.(as far as I think, the mesh is fine enough)

The error is as below:

Failed to find consistent initial values.
Segregated group 1

Undefined value found.
Undefined value found in the stiffness matrix..
For mesh-case 1 there are 111061 equations giving NaN/Inf in the matrix rows for the variable mod1.H.
and similarly for the degrees of freedom, NaN/Inf in the matrix columns.
Last time step is not converged.

I'm not clear about the numerical methods. It seems that the solving processes of 2D and 3D are quite different. 3D used segregated groups, which is not shown in 2D process.

I need your help.
Thanks in advance.
Mingjie Jia

Structural Mechanics Equation System

$
0
0
I am trying to understand how my structural mechanics model is computing its solution. The model solves and gives me a reasonable solution. However, I don't understand how it is arriving at its answer. I am looking at the Subdomain Settings - Equation System panel where it shows the PDE which it is solving. I am viewing it in coefficient form, and all of the coefficients are set to zero. How am I arriving at a solution when the PDE which is being solved is zero?

Ryan

Sloshing Tank with FSI

$
0
0
Hello All

I calculate the waves of a water tank and then damp with a structure the waves. At the beginning I start with a single excitation an then wait and just let the gravity do the job. My model is similar to the sloshing tank example but rebuild with the FSI. I added then a bar which is made of rubber with a damping loss factor.

Now the problem: When the young modulus is 10e9 the model works fine, but the bar moves only for um. When I reduce the young modulus to 1e9, the model didn't converge. Does anybody know how I can make this model more robust? I would like to see the bar swinging within the waves and set the damping factor an the size of the bar to optimally damp the waves.

My model is attached.

2D geometry but illuminated from below

$
0
0
Hi,

I want to create a 2D geometry but illuminate it from below, or out of plane, with a plane wave. Whenever I try and set the Background Field as exp(-i*emw.k0*z) (in x direction) I get errors as it fails to evaluate the variable z.

Any help is greatly appreciated!

Daniel

problem with port boundary condition

$
0
0
Hi,

I am trying to simulate a coupler so I need to set a port as inport and I have to select "wave excitation at this port " check box but it's inactive.
How can I activate it?

Thanks,

sara

Flux units in surface reaction interface module, and modeling of surface diffusion in axi-symmetric system

$
0
0
Hi,
I am trying to model a bulk - surface reaction-diffusion system (diffusion in the solution, diffusion on the surface, and adsorption/desorption which couples the two together).
The system has a cylindrical symmetry so I want to solve it in an axi-symmetric configuration, and I have two questions regarding the implementation of such system:
1. Initially I constructed the system similarly to the example of Transport and Adsorption case study (Model Library > COMSOL Multiphysics > Chemical Engineering>Transport and Adsorption) which is quite similar. In this case the surface is modeled using a General Form Boundary PDE (gb). I found out that the surface concentration was not solved correctly, and I suspect it is related to the axial symmetry.
Did anyone encountered such a problem?
Does general form PDE work with 1D axisymmetric configuration?

2. I rebuilt the model with surface reaction interface instead of general PDE, and this indeed gives a correct solution for the surface concentration.
I want to check a flux on a point along the surface using the total flux magnitude in the post processing (chsr.tfluxMag_). However, it reports it with the weird units [ mol^2/(m^3*s) ] ,which are definitely not flux units.
The value is also not in agreement to the predicted flux (in fact it is the correct flux*constant).
What is the "flux" which is given there?
Is there another way to calculate the exact flux in this interface?

Thanks!
Henry

Error: Failed to compute elastoplastic strain variables

$
0
0
I am using solid mechanics and Darcy's law modules. In the solid mechanics, I add soil plasticity model and use D-P model, but calculation is ending by this error. Could you tell me what is the problem of this?

solver settings

$
0
0
Hi
I'm using two physics piezoelectric and moving mesh. I want to change the thickness of the PZT layer using moving mesh and want to calculate the voltage value for each thickness. how to do it in comsol 4.2a. I can do this in comsol 3.5a. Moving mesh is solved using stationary solver and PZT is solved using frequency solver. I used parametric segregated solver in 3.5a to compute the above said operations.
Viewing all 26527 articles
Browse latest View live