Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Non-physical results with geometric transformations (TransFile) #245

Open
izod2 opened this issue Oct 4, 2023 · 0 comments
Open

Non-physical results with geometric transformations (TransFile) #245

izod2 opened this issue Oct 4, 2023 · 0 comments

Comments

@izod2
Copy link

izod2 commented Oct 4, 2023

The concept behind modifying geometric orientations is documented here, but it does not appear to work.

Here is a MWE example for a scuff-scatter calculation involving two side-by-side spheres. Relevant files are copied below.

Case 1: No TransFile is used (it is commented out in Args):
Output of the PFT file makes sense (as expected, both spheres have the same absorbed/scattered powers, etc.

# omega surface-label absorbed power (watts) scattered power (watts) ...
1.1 Sphere1 8.085983e-05 9.520283e-03 ...
1.1 Sphere2 8.081348e-05 9.501876e-03 ...

Case 2: A TransFile displaces both spheres along y:

Spheres.trans

TRANS Default
TRANS 1.0 OBJECT Sphere1 DISP 0 1 0 OBJECT Sphere2 DISP 0 1 0

Output of the PFT file should be the same as Case 1 - this is physically the same configuration. Instead, reported powers are not only different for the two spheres, but negative values show up for Sphere2 (?!)

# omega transform surface-label absorbed power (watts) scattered power (watts) ...
1.1 Default Sphere1 7.719171e-05 9.213488e-03 ...
1.1 Default Sphere2 -5.444921e-04 8.524173e-03 ...
1.1 1.0 Sphere1 7.719171e-05 9.213488e-03 ...
1.1 1.0 Sphere2 -5.444921e-04 8.524173e-03 ...

Does anyone know what is going on here? I tried several different configurations/frequencies (for both regular and abbreviated .trans file syntax) but the results never seem to make sense.

Relevant files:
Args

ARGS="${ARGS} --geometry    Spheres_327.scuffgeo"
ARGS="${ARGS} --Omega       1.1"
ARGS="${ARGS} --TransFile   Spheres.trans"   # comment out 
ARGS="${ARGS} --PFTFile     Spheres_327.pft"
ARGS="${ARGS} --pwDirection     0 0 1"
ARGS="${ARGS} --pwPolarization  1 0 0"

Spheres_327.scuffgeo

MATERIAL Gold
    wp = 1.37e16;
    gamma = 5.32e13;
    Eps(w) = 1 - wp^2 / (w * (w + i*gamma));
ENDMATERIAL

OBJECT Sphere1
	MESHFILE Sphere_327.msh
	MATERIAL Gold
ENDOBJECT

OBJECT Sphere2
	MESHFILE Sphere_327.msh
	MATERIAL Gold
        DISPLACED 2.5 0 0 
ENDOBJECT

(Sphere_327.msh is from the solid-sphere built-in example in scuff-em)

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant