Scribd
Upload
135 views2 pages

LPG Leak Simulation Analysis

This document contains a simulation of the dispersion of an LPG gas leak from a source point at PT Mandom Indonesia. The simulation defines variables and equations to model the diffusion and flow of the gas (C) and pressure (P) over time. Boundary conditions with no flux are applied on the edges of the simulation region. Concentration (C) is initialized at 0 and injected at 1 from the source point. The simulation will run from 0 to 308,000 seconds in 100 second intervals, outputting contour plots and velocity vectors of the dispersing gas plume over 5 cycles.

Uploaded by

Muhammad Hanief
Copyright
© All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
135 views2 pages

LPG Leak Simulation Analysis

This document contains a simulation of the dispersion of an LPG gas leak from a source point at PT Mandom Indonesia. The simulation defines variables and equations to model the diffusion and flow of the gas (C) and pressure (P) over time. Boundary conditions with no flux are applied on the edges of the simulation region. Concentration (C) is initialized at 0 and injected at 1 from the source point. The simulation will run from 0 to 308,000 seconds in 100 second intervals, outputting contour plots and velocity vectors of the dispersing gas plume over 5 cycles.

Uploaded by

Muhammad Hanief
Copyright
© All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
  • Script Introduction & Setup
  • Script Continuation & Results

Aloha Script :

title
'Simulasi dispersi Kebocoran Gas LPG di PT Mandom
Indonesia'
select
nominmax
contours= 50
PAINTREGIONS
variables
C(threshold=0.0000000001)
P(threshold=0.000000000001)

definitions
A= 0.01
B = -1
epsvisc = 1.e-6*10
muw = 0.01861598
muo = 0.017668201
M = 10 !C^2/muw + (1-C)^2/muo
D = 0.001
! D >= 1 for explosion, D <= 1e-3 for leak
WF = 1
{wind factor, WF = 1 if there's wind and 0
if no wind}
K =6e-11
K2 = 6e-11
krw = C^2/muw
Pin = 1.01e5 !pascal
Pexhaust = 1496 !pascal
Pout = 1496
phi =1
hour = 60*60
day = hour * 24

initial values
C = 0
P = Pin

equations
C: phi*dt(C) - div(K*krw*grad(p)) -
epsvisc*div(grad(C)) = 0
p: div(K*M*grad(p)) = 0

boundaries
region 1
start(0,0)

natural(C) = 0 line to (3,0)


natural(C) = 0 line to (3,0.7)
natural(C) = 0 Natural(p) = 0
line to (4.3,0.7) line to (4.3,1.8)
natural(C) = 0 Natural(p) = 0
line to (0,1.8)

line to close

start (3.3,0.7)
natural(C) = 0 Natural (p) = 0
line to (3.35,0.7) line to (3.35,0.9) line to (3.3,0.9) line
to close

start (3,1.05)
natural(C) = 0 Natural (p) = 0
line to (3.05,1.05) line to (3.05,1.8) line to (3,1.8) line to
close

start (3.3,1.05)
natural(C) = 0 Natural (p) = 0
line to (3.35,1.05) line to (3.35,1.8) line to (3.3,1.8) line
to close

start (3.4, 0.95)


line to (3.41,0.95) line to (3.41,0.96) line to (3.4,0.96)
value (C) = 1 line to (3.4,0.95) line to close
! SUMBER KEBOCORAN LPG

TIME 0 to 308000 by 100


monitors
for cycle = 5
vector(-K*M*grad(p)) norm as "Flow Velocity"

plots
for cycle = 5

contour(min(max(c,0),1)) as "Chemical Dispersion Contour


Profile non painted" range(-0.26,1)
contour(min(max(c,0),1)) as "Chemical Dispersion Contour
Profile Painted" range(-0.26,1) painted nominmax
surface(min(max(c,0),1)) as "Chemical Dispersion Contour
Profile 3D" range(-0.26,1) VIEWPOINT(3, 0.7, 30) !fixed range
(0,10)

histories
end

Aloha Script : title 'Simulasi dispersi Kebocoran Gas LPG di PT Mandom Indonesia' select nominmax contours= 50 PAINTREGION
line to (4.3,0.7) line to (4.3,1.8) natural(C) = 0 Natural(p) = 0 line to (0,1.8) line to close start (3.3,0.7) natural(C) =

You might also like