Main #274

JefferssonAgudelo · 2022-07-08T19:50:12Z

No description provided.

Not pretty, this really needs some consolidation, but it'll do for now.

…harris and turbulent configuration

…ble in up and zero in bottom

…e x direction. It seems the velocity is not properly setup. Also, the conducting and reflectic Bndc generate gauss and continuity errors

…if mime!=0. There is an issue with the plasma. There are no outflows but inflowsgit add .

… case now

…hat shouldnt be there. It might be due to the boundary conditions. Also there are some outputs in which the fields and moments files are return with an adittional 2 and I dont know why is this happening

…le which force the program to write the file name including an extra number. On the other hand the harris file has now two separate configurations for xz and yz geometries. There is still an issue with the Dim_xyz. It is not running because there is something worng with the lhs and rhs definition somewhere.

…ementation of langevin antena (still alot of work to do). However, it is necessary to get it as a usable function that can be called by the fields, to compute the current from the external magnetic field since it is the external current the one that is actully included into the pushing of the fields. Basically, to let the code do the selfconsistent calculation of the curl of Bext so it gets the correct fc or cc quantities. INCLUDE THE TIME STEPS AND THE ITERATIVE VALUES

…configuration for a pair plasma

…s no drifting speed. So it might not be necessary to change the way to inject particles in the setup file

…he issue now is to include the crd and the iteration

no actual changes

we permuted x, y, z, so that the xy perturbations show up in the 2d yz plane

Psc harris turb

since that's needed when calculating the curl

results look reasonable, at least

updated notebook

…x there is a segmentation fault when pushing particles

…I put them by hand

…he same mass ratio for direct comaprison. The temperatures might be different tho

… remove the 0. on the external current calculation

germasch and others added 30 commits May 15, 2022 21:47

turb: start out as copy from flatfoil_yz, get rid of CASEs

cf4f6df

turb: set background density to 1, don't inject initially

e346a83

turb: get rid of heating

ac3a4bb

turb: get rid of injection

64a0ded

turb: get rid of ELECTRON_HE

4161302

turb: now we have a minimal case

aaeab80

turb: switch output to WriterMRC (XDMF/HDF5)

8ea58f2

turb: setup intermediate mflds_alfven first

5e4f6ea

turb: use 6-component field for Alfven perturbation

ef4b5d7

turb: looks like initial H and V pert works

134ed56

turb: avoid hacky cast to double

39b366c

turb: add jupyter notebook

65d06f0

turb: fix compile with cuda

2d52c8a

Not pretty, this really needs some consolidation, but it'll do for now.

adding psc_turb_harris_xz.cxx This is the file that will include the …

c50b1ed

…harris and turbulent configuration

Including the up and bo cases doesn't really work. The density is dou…

688e588

…ble in up and zero in bottom

The harris set up is almost done. However, there is a net drift on th…

170ff4f

…e x direction. It seems the velocity is not properly setup. Also, the conducting and reflectic Bndc generate gauss and continuity errors

the drifting is because ions and electrons have different velocities …

2da3a41

…if mime!=0. There is an issue with the plasma. There are no outflows but inflowsgit add .

The psc_harris_xz.cxx seems to work for the simple pair plasma harris…

05f5573

… case now

It is still not properly working. It seems there is an hy component t…

a1a3d4e

…hat shouldnt be there. It might be due to the boundary conditions. Also there are some outputs in which the fields and moments files are return with an adittional 2 and I dont know why is this happening

The file turb_harris_yz.cxx is now ready as a standar harris current …

bef0561

…configuration for a pair plasma

The psc_turb_harris_yz.cxx works now for mime greater than 1. There i…

3ac0274

…s no drifting speed. So it might not be necessary to change the way to inject particles in the setup file

the file now have the correct calculation of Bext_x_r and Bext_y_r. T…

fcd2e04

…he issue now is to include the crd and the iteration

extra comments

bc8731c

auto-format only

cfc253d

no actual changes

evolve bn_k over time, set initial magnetic field

8624d6a

we permuted x, y, z, so that the xy perturbations show up in the 2d yz plane

initialize in 3d, via vector potential, div B == 0 works

59f2fba

Merge pull request #1 from germasch/psc_harris_turb

713e16d

Psc harris turb

updated notebook

cec882e

JefferssonAgudelo and others added 26 commits June 22, 2022 14:52

the div B works but not the curl of B. There is an issue with the size

7cec5c9

Merge remote-tracking branch 'jeffersson1/main' into psc_harris_turb

485b0a9

clang-format

f348539

calculate H including one ghost point on the higher side

3cec8ee

since that's needed when calculating the curl

second curl needs to be done going to -1 to get staggering right

665fbad

initialize random numbers on rank 0 and broadcast

d20c0f0

put random ext B evolution back in

bdc0556

calculate Az from random modes again, rather than test func

7dddb3c

use MY_{ELECTRON,ION} constants rather than hardcoded numbers

0c6c0a3

run 10 steps, output every step

ea3b0a7

put initial condition back to Harris-only, with no pert/random B

ca5a793

update jupyter notebook

41c7ea6

put initial perturbation back to turbulent modes, make box all periodic

c0bc660

use separate functions for calculating curl Az, H

a59bbcc

add ExtCurrent to psc as a hook to be able to modify fields

2a39bf5

run 100 steps and remove particles (setting plasma density to zero)

2c1babc

add external currents

cbef35c

results look reasonable, at least

update notebook

4d44dca

Merge pull request #2 from germasch/psc_harris_turb

251c961

updated notebook

Including the turb_alfven

7ed5092

When including the perturbation for the velocity field in turb_alf.cx…

e77ef10

…x there is a segmentation fault when pushing particles

the turb_alf now works. There is an issue with the random numbers so …

fcb33de

…I put them by hand

The turb_alf and the turb_langevin are set on the same geometry and t…

2e3ac46

…he same mass ratio for direct comaprison. The temperatures might be different tho

Including the correct background and harris populations

efceeb7

testing where the commit goes to

668682f

The drift velocity were too low. Also, to run the driving remember to…

716753c

… remove the 0. on the external current calculation

Provide feedback

Saved searches

Use saved searches to filter your results more quickly

Main #274

Main #274

JefferssonAgudelo commented Jul 8, 2022

Main #274

Are you sure you want to change the base?

Main #274

Conversation

JefferssonAgudelo commented Jul 8, 2022