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.

Sign up for GitHub

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

Jump to bottom

FDS issue wioth main.f90 and geom.f90 #12796

Open
Dancwr2911 opened this issue Apr 13, 2024 · 7 comments
Open

FDS issue wioth main.f90 and geom.f90 #12796

Dancwr2911 opened this issue Apr 13, 2024 · 7 comments
Assignees
@marcosvanella

Comments

@Dancwr2911
Copy link

I am with uni and using FDS 6.7.5
The error I get here is:

forrtl: severe (157): Program Exception - access violation
Image PC Routine Line Source
fds.exe 00007FF6533E8449 COMPLEX_GEOMETRY_ 12943 geom.f90
fds.exe 00007FF6533F25F4 Unknown Unknown Unknown
fds.exe 00007FF65311062B MAIN__ 779 main.f90
fds.exe 00007FF6530CF09E Unknown Unknown Unknown
fds.exe 00007FF653A66664 Unknown Unknown Unknown
KERNEL32.DLL 00007FFAF595257D Unknown Unknown Unknown
ntdll.dll 00007FFAF738AA58 Unknown Unknown Unknown

It is not allowing me to attach so here is my code:

#Title
&HEAD CHID='assignment_fire', TITLE='assignment_fire' /

#Mesh
&MESH XB= 0.0, 70.0, 0.0, 50.0, 0.0, 9.0, IJK= 280, 200, 36 /

#Time
&TIME T_END=180. /
DT=0.1

#Boundary Conditions

#Fire properties
&REAC FUEL = 'PROPANE', SOOT_YIELD=0.02

&SURF ID='FIRE', RAMP_Q='tsquared', HRRPUA=168.8, COLOR='RED' /
&RAMP ID='tsquared', T= 0.0, F=0.0000 /
&RAMP ID='tsquared', T= 10.0, F=0.0278 /
&RAMP ID='tsquared', T= 20.0, F=0.1111 /
&RAMP ID='tsquared', T= 30.0, F=0.2500 /
&RAMP ID='tsquared', T= 40.0, F=0.4444 /
&RAMP ID='tsquared', T= 50.0, F=0.6944 /
&RAMP ID='tsquared', T= 60.0, F=1.0000 /

&VENT XB= 39.5,41,14.5,16,0,0, SURF_ID='FIRE', COLOR='RED'/

#Devices
&DEVC ID='CO-rest', XYZ=47,11.5,8.25, QUANTITY='FED', SPEC_ID='CARBON DIOXIDE' /
&DEVC ID='CO2-rest', XYZ=47,11.5,8.25 QUANTITY='FED', SPEC_ID='CARBON MONOXIDE' /
&DEVC ID='O-rest', XYZ=47,11.5,8.25, QUANTITY='FED', SPEC_ID='OXYGEN' /
&DEVC ID='visibility-rest', XYZ=47,11.5,8.25, QUANTITY='VISIBILITY' /
&DEVC ID='CO-lift', XYZ=50.5,23.5,8.25, QUANTITY='FED', SPEC_ID='CARBON DIOXIDE' /
&DEVC ID='CO2-lift', XYZ=50.5,23.5,8.25, QUANTITY='FED', SPEC_ID='CARBON MONOXIDE' /
&DEVC ID='O-lift', XYZ=50.5,23.5,8.25, QUANTITY='FED', SPEC_ID='OXYGEN' /
&DEVC ID='visibility-lift', XYZ=50.5,23.5,8.25, QUANTITY='VISIBILITY' /
&DEVC ID='CO-showroom', XYZ=21,40,8.25, QUANTITY='FED', SPEC_ID='CARBON DIOXIDE' /
&DEVC ID='CO2-showroom', XYZ=21,40,8.25, QUANTITY='FED', SPEC_ID='CARBON MONOXIDE' /
&DEVC ID='O-showroom', XYZ=21,40,8.25, QUANTITY='FED', SPEC_ID='OXYGEN' /
&DEVC ID='visibility-showroom', XYZ=21,40,8.25, QUANTITY='VISIBILITY' /

#Obstructions
&SURF ID='Triangle_Surface', RGB=0,255,0 /
&SURF ID='RAMP_SURFACE', RGB=0,0,255 /
&SURF ID='WALL', RGB=255,0,0 /

#Second floor

&OBST XB=31.0, 70.0, 0.0, 2.0, 3.5, 5.5, COLOR='GREEN'/

&GEOM ID='Triangle_Extrusion', SURF_ID='Triangle_Surface',
VERTS=31.0, 1.75, 3.5,
31.0, 22.0, 3.5,
70.0, 1.75, 3.5,
POLY=1, 2, 3,
EXTRUDE=-2.0
COLOR='GREEN' /

&GEOM ID='Outside_Wall', SURF_ID='WALL',
VERTS=24.0, 15.0, 0.0,
24.0, 15.0, 3.5,
54.0, 30.0, 3.5,
54.0, 30.0, 0.0
POLY=1, 2, 3, 4,
EXTRUDE=10.0
COLOR='BLACK' /

&OBST XB=25.0, 36.0, 22.0, 30.0, 3.5, 5.5, COLOR='GREEN'/
&OBST XB=35.0, 53.75, 27.0, 30.0, 3.5, 5.5, COLOR='GREEN'/
&OBST XB=49.0, 52.0, 23.0, 28.0, 3.5, 5.5, COLOR='GREEN'/
&OBST XB=25.0, 36.0, 0.0, 30.0, 3.5, 5.5, COLOR='GREEN'/
&OBST XB=0.0, 32.0, 22.0, 50.0, 3.5, 5.5, COLOR='GREEN'/
&OBST XB=31.0, 42.0, 30.0, 50.0, 3.5, 5.5, COLOR='GREEN'/

#Solid Walls

#Building Boundaries
&OBST XB=0.0, 25.0, 21.0, 22.0, 0.0, 9.0, COLOR='BLACK'/
&OBST XB=24.0, 25.0, 0.0, 22.0, 0.0, 9.0, COLOR='BLACK'/
&OBST XB=42.0, 43.0, 30.0, 50.0, 0.0, 9.0, COLOR='BLACK'/
&OBST XB=42.0, 54.0, 30.0, 31.0, 0.0, 9.0, COLOR='BLACK'/

#Slanted wall
&OBST XB=67.4, 70.0, 6.0, 10.0, 0.0, 9.0, COLOR='BLACK'/
&OBST XB=64.9, 67.5, 10.0, 14.0, 0.0, 9.0, COLOR='BLACK'/
&OBST XB=62.4, 65.0, 14.0, 18.0, 0.0, 9.0, COLOR='BLACK'/
&OBST XB=59.9, 62.5, 18.0, 22.0, 0.0, 9.0, COLOR='BLACK'/
&OBST XB=57.4, 60.0, 22.0, 26.0, 0.0, 9.0, COLOR='BLACK'/
&OBST XB=53.9, 57.5, 26.0, 30.0, 0.0, 9.0, COLOR='BLACK'/

#Stairs
&GEOM ID='Ramp', SURF_ID='RAMP_SURFACE',
VERTS=52.0, 26.0, 5.5,
54, 23, 4.25,
56, 23.5, 4.25,
54, 26.5, 5.5,
POLY=1, 2, 3, 4,
EXTRUDE=-0.25,
COLOR ='RED' /

&GEOM ID='Landing', SURF_ID='RAMP_SURFACE',
VERTS=54, 23, 4.25,
55, 21.5, 4.25,
57, 22, 4.25,
56, 23.5, 4.25,
POLY=1, 2, 3, 4,
EXTRUDE=-0.25,
COLOR ='RED' /

&GEOM ID='Ramp', SURF_ID='RAMP_SURFACE',
VERTS=55, 21.5, 4.25,
57, 18.5, 3.0,
59, 19, 3.0,
57, 22, 4.25,
POLY=1, 2, 3, 4,
EXTRUDE=-0.25,
COLOR ='RED' /

&GEOM ID='Landing', SURF_ID='RAMP_SURFACE',
VERTS=57, 18.5, 3.0,
58, 17, 3.0,
60, 17.5, 3.0,
59, 19, 3.0,
POLY=1, 2, 3, 4,
EXTRUDE=-0.25,
COLOR ='RED' /

&GEOM ID='Ramp', SURF_ID='RAMP_SURFACE',
VERTS=58, 17, 3.0,
60, 14, 1.75,
62, 14.5, 1.75,
60, 17.5, 3.0,
POLY=1, 2, 3, 4,
EXTRUDE=-0.25,
COLOR ='RED' /

&GEOM ID='Landing', SURF_ID='RAMP_SURFACE',
VERTS=60, 14, 1.75,
61, 12.5, 1.75,
63, 13, 1.75,
62, 14.5, 1.75,
POLY=1, 2, 3, 4,
EXTRUDE=-0.25,
COLOR ='RED' /

&GEOM ID='Ramp', SURF_ID='RAMP_SURFACE',
VERTS=61, 12.5, 1.75,
63, 9.5, 0.5,
65, 10, 0.5,
63, 13, 1.75,
POLY=1, 2, 3, 4,
EXTRUDE=-0.25,
COLOR ='RED' /

&GEOM ID='Landing', SURF_ID='RAMP_SURFACE',
VERTS=63, 10, 0.5,
66, 10, 0.5,
66, 5, 0.5,
63, 5, 0.5,
POLY=1, 2, 3, 4,
EXTRUDE=0.25,
COLOR ='RED' /

&GEOM ID='Ramp', SURF_ID='RAMP_SURFACE',
VERTS=63, 10, 0.5,
60, 10, 0.0,
60, 5, 0.0,
63, 5, 0.5,
POLY=1, 2, 3, 4,
EXTRUDE=-0.25,
COLOR ='RED' /

#Internal Walls

#Restaurant
&OBST XB=47.5, 48.5, 0.0, 11.0, 5.5, 9.0, COLOR='BLUE'/
&OBST XB=36.0, 42.0,14.5, 14, 5.5, 9.0, COLOR='BLUE'/

#DIAGONALS
&MULT ID='WALL', DX=-0.5, DY=0.25, N_LOWER=0, N_UPPER=8/
&OBST XB=45.5,46.5,12.0,12.25,5.5,9.0, MULT_ID='WALL', COLOR='BLUE'/
&MULT ID='DOOR', DX=-0.5, DY=0.25, N_LOWER=0, N_UPPER=4/
&OBST XB=47.5,48.5,11.0,11.25,8,9, MULT_ID='DOOR', COLOR='BLUE'/

#Studios behind Restaurant
&OBST XB=36, 35.5,5.0, 14.0, 5.5, 9.0, COLOR='BLUE'/
&OBST XB=36.0, 34.0,5.5, 5, 5.5, 9.0, COLOR='BLUE'/
&OBST XB=33.5, 34,5.0, 22.5, 5.5, 9.0, COLOR='BLUE'/
&OBST XB=34.0, 25,22, 22.5, 5.5, 9.0, COLOR='BLUE'/

#Gallery Section
&OBST XB=32, 42,30.0, 30.5, 5.5, 9.0, COLOR='BLUE'/
&OBST XB=32.0, 32.5,30, 39, 5.5, 9.0, COLOR='BLUE'/
&OBST XB=28, 32,39.0, 38.5, 5.5, 9.0, COLOR='BLUE'/
&OBST XB=20.0, 20.5,22, 50, 5.5, 9.0, COLOR='BLUE'/

#Second Enclosure
&MULT ID='WALL2', DX=-0.5, DY=-0.5, N_LOWER=0, N_UPPER=8/
&OBST XB=28.5,29,40.0,39.5,5.5,9.0, MULT_ID='WALL2', COLOR='BLUE'/
&OBST XB=24.5,25,48.0,47.5,5.5,9.0, MULT_ID='WALL2', COLOR='BLUE'/
&MULT ID='WALL3', DX=-0.25, DY=0.5, N_LOWER=0, N_UPPER=16/
&OBST XB=28.5,28.75,40.0,39.5,5.5,9, MULT_ID='WALL3', COLOR='BLUE'/
&OBST XB=24.5,25,36.0,35.5,5.5,9, MULT_ID='WALL3', COLOR='BLUE'/

&SLCF PBY = 15, QUANTITY='TEMPERATURE'/
&SLCF PBY = 15, QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON DIOXIDE' /
&SLCF PBY = 15, QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON MONOXIDE' /
&SLCF QUANTITY='VISIBILITY', PBY=15/

&TAIL /
@gforney
Copy link
Contributor

gforney commented Apr 13, 2024 via email

@Dancwr2911
Copy link
Author

#Title.txt

@drjfloyd
Copy link
Contributor

first try your case with the current release

@Dancwr2911
Copy link
Author

tried and no luck

@marcosvanella
Copy link
Contributor

marcosvanella commented Apr 15, 2024 via email

@marcosvanella marcosvanella self-assigned this Apr 15, 2024
@marcosvanella
Copy link
Contributor

@Dancwr2911 swap your single mesh by this one:

&MULT ID='MyMult', DX=17.5, DY=12.5, I_UPPER=3, J_UPPER=3 /
&MESH XB= 0.0, 17.5, 0.0, 12.5, 0.0, 9.0, IJK= 70, 50, 36, MULT_ID='MyMult' /

Note you will have to run the case with 16 cores. A couple of things we need to work on:

  • Figure out a way to alert the user about large memory use of Poisson factorization matrix -> Pardiso info.
  • Either enable interaction of zero thinkness OBSTs with GEOMs or set Thicken to all OBSTs.

@obscureed
Copy link

I am not associated with NIST or FDS developers, so you are welcome to prioritize the advice from the real experts. I took a quick look at your model and can offer a few comments:
• FDS 6.7.5 is old (released 2020). As @drjfloyd said, use a recent release, specifically 6.9.1. It's not a difficult installation process (in my experience).
• Having said that, beta features bring complications. If you can proceed without them, you are more likely to make progress.
• Your case runs (slowly) in FDS6.9.1 when I remove all the &GEOM commands. (Actually, to my surprise, it runs slowly with the &GEOM commands as well, needing only 10GB or so of RAM.)
• If you decompose the mesh into several chunks, this brings the opportunity to reduce the cell count. You have several enclosed spaces that are not connected to the fire compartment. (And they have no other influence from the fire compartment, because you have adiabatic walls and no conduction.) These are presumably not useful parts of your model, so you can assemble &MESH blocks that don't fill these parts of space.
• Where there are still irrelevant enclosed spaces in the model (for example, on the other side of a stair-stepped angled wall), I would generally try to fill them in. Solid blocks of &OBST take less computational effort than idle regions of fluid, and have no potential to cause instability. Actually, enclosed fluid with not even slight heating from conduction should have zero flow, so perhaps instability is not a risk in this case. (FDS is better than many CFD codes at letting stagnant fluid do nothing. Usually, the risk with enclosed fluid spaces is that they will get some heat input from conduction, the air will expand and the pressure will zoom up.)
• Your fire compartment is also fully enclosed -- walls all around. The fire is not large compared to the space, but you will find that pressure builds up inside this enclosed space, almost certainly unrealistically. It will help stability (and realism) if you add a small amount of leakage somewhere. (My preferred method for this would be a small patch of "&VENT SURF_ID='OPEN', XB=..." on the exterior boundary, probably at least 2 pixels * 2 pixels. No need to involve the HVAC solver.)
• See the FDS User Guide about the default thermal boundary conditions for walls -- search for 'INERT'. This is almost never what you want; FDS declines to offer you a more realistic option, because it varies so much from case to case. You can see concrete and plasterboard, for example in the Bittern_Sprinkler_Experiments cases in the Validation examples.
• For an exploratory model, I would add a lot more postprocessing. (For example: slices of several quantities, including velocity vectors, at several heights and through a few vertical slices. The fire is unlikely to get close to underventilation, but you should prove that to yourself with a slice of oxygen mass fraction through the fire bed. Whenever I see a doorway in a model, I want to know the flowrate through it -- not easy for an angled doorway, but you'll find a way.)
• I wouldn't use red color for geometry next to a red fire bed.
• I don't know what you're hoping for when you specify a species (SPEC_ID=...) in "&DEVC QUANTITY='FED' ...". I suspect this is an example where FDS lets you specify something but does not use that input at all.
• Talking of which, your text "DT=0.1" early in the input file does nothing -- it is silently interpreted as a comment. Anything outside of a known command "&XXXX ... /" is a comment. In my view, this is a dangerous syntax, but it's probably too well embedded to challenge. A good way to start a comment (including at the end of a line inside a multiline command) is with "!".

Good luck!
Ed

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Assignees

@marcosvanella marcosvanella

Labels
None yet
Projects
None yet
Milestone
No milestone
Development

No branches or pull requests
5 participants
@marcosvanella @gforney @drjfloyd @obscureed @Dancwr2911