Search found 214 matches

by mika
09 Aug 2022, 17:46
Forum: General
Topic: Derivation magnetic co-energy
Replies: 7
Views: 3117

Re: Derivation magnetic co-energy

As further information, is this a simulation done in the frequency domain?
by mika
07 Jul 2022, 15:02
Forum: General
Topic: thermal expansion
Replies: 74
Views: 2508

Re: thermal expansion

The sif files repeats the constraint Displacement 3 = 0 twice in the same BC section. You should at least add constraints which prevent both rigid body translations in the X-Y plane and rotations around the Z-axis to have a unique solution. The lack of uniqueness may cause convergence problems. In a...
by mika
06 Jul 2022, 17:15
Forum: ElmerSolver
Topic: Definition of Concentration Diffusivity
Replies: 3
Views: 188

Re: Definition of Concentration Diffusivity

It seems that the advection-diffusion solver reads the diffusion coefficient from the keyword command in the form my_solver_variable Diffusivity = Real ... where the part "my_solver_variable" should be the string which is used as the value of the solver variable. That is, the command Conce...
by mika
30 Jun 2022, 15:16
Forum: ElmerSolver
Topic: ShellSolver calculation of stress strain
Replies: 2
Views: 155

Re: ShellSolver calculation of stress strain

Calculate Stresses = True Calculate Strains = True These commands work in the case of 3-D models, but strain/stress computation has not yet been implemented into the solver module ShellSolver. The old facet shell solver module FacetShellSolve (which is undocumented and whose use in shell analysis m...
by mika
19 Jan 2022, 19:52
Forum: ElmerSolver
Topic: Problem with transient heat transfer
Replies: 9
Views: 897

Re: Problem with transient heat transfer

Do you maybe have some further ideas to better physically introduce the flux in FEM methods that I could try out? If not, do you then suggest to check out finite volumes for these types of problems? I believe the mixed FEM formulation presented in the section "Mixed Approximation of the Poisso...
by mika
14 Dec 2021, 10:27
Forum: ElmerSolver
Topic: Modelpde with no stabilization
Replies: 3
Views: 609

Re: Modelpde with no stabilization

In Elmer the basic element definition for using the p-version of FEM is just "p:k", where k defines the order of approximation. Elmer then assigns automatically the right numbers of DOFs associated with different geometric entities (edges, faces and element interiors). Mixing the definitio...
by mika
09 Dec 2021, 19:16
Forum: ElmerSolver
Topic: Modelpde with no stabilization
Replies: 3
Views: 609

Re: Modelpde with no stabilization

It should be possible to apply bubble augmentation by giving a p-element definition. For example

Element = "p:1 b:1"

creates one elementwise bubble function. This construct doesn't depend on the keyword Bubbles and needs the background mesh consisting of the lowest-order elements.

--Mika
by mika
08 Sep 2021, 16:02
Forum: ElmerGUI
Topic: Magnetization
Replies: 7
Views: 2545

Re: Magnetization

Rich_B wrote: 08 Sep 2021, 15:45 The WhitneyAVsolver is the only solver that uses edge elements, correct?
The modules EMWaveSolver and VectorHelmholtz for electromagnetic waves also use edge finite elements.

-- Mika
by mika
05 Sep 2021, 13:29
Forum: Software development
Topic: GOLF rheology in ParStokes Solver?
Replies: 8
Views: 3880

Re: GOLF rheology in ParStokes Solver?

I committed changes which should be quite close to what are needed here, see https://github.com/ElmerCSC/elmerfem/commit/8b78545b963eae289b74320b9fc91bb298d39eba and a small fix https://github.com/ElmerCSC/elmerfem/commit/9a9a5e3537b14a7546b7dd2f99a1d4678d1e4c05 I followed the convention that the st...
by mika
01 Sep 2021, 09:49
Forum: Software development
Topic: GOLF rheology in ParStokes Solver?
Replies: 8
Views: 3880

Re: GOLF rheology in ParStokes Solver?

I think the book mentioned in the paper doesn't really treat anisotropy. If you can create a code which enables to represent the material law as T + p I = C^{ijkl}D_{kl}, with the righ-hand side expressed by using the Voigt notation (that is, the tensor C is then expected to be a (6x6)-matrix), I co...