Claim Missing Document
Check
Articles

Pemodelan Magnetotellurik 2D Menggunakan Metode Elemen Batas Imran Hilman Mohammad; Wahyu Srigutomo; Doddy Sutarno
Jurnal Matematika & Sains Vol 16, No 2 (2011)
Publisher : Institut Teknologi Bandung

Show Abstract | Download Original | Original Source | Check in Google Scholar

Abstract

Dalam permasalahan elektromagnetik untuk aplikasi geofisika, rangkaian persamaan Maxwell dapat disederhanakan menjadi persamaan Helmholtz, sehingga dapat dicari solusi permasalahan menggunakan berbagai skema numerik. Metode elemen batas merupakan metode numerik untuk memecahkan persamaan diferensial parsial yang telah dikembangkan dalam beberapa dekade ini untuk memecahkan berbagai permasalahan medan elektromagnet. Metode ini memiliki keunikan dibandingkan metode numerik lain untuk memecahkan persamaan diferensial parsial, karena hanya membutuhkan diskretisasi pada bidang-bidang batas domain pemodelan. Solusi pada batas domain dapat digunakan untuk mencari solusi pada seluruh  domain pemodelan, membuat metode ini memiliki algoritma numerik yang sangat efisien. Dalam makalah ini, metode elemen batas digunakan untuk menghitung respon magnetotellurik 2D dalam bentuk resistivitas semu dan fasa impedansi modus TE (transverse electric). Pemodelan elemen hingga digunakan sebagai pembanding hasil pemodelan elemen batas yang dikembangkan. Kata kunci : Metode elemen batas, Magnetotelurik 2D.   2D Magnetotelluric Modelling using Boundary Element Method Abstract In most formulations of electromagnetic methods in geophysics, the set of Maxwell’s equations can be simplified into Helmholtz’s equation which leads to possibility of solving electromagnetic problems using various numerical schemes. The boundary element method is a numerical method for solving partial differential problems which in the last several decades has been applied in electromagnetic problems. The method poses unique advantage in comparison to other methods; it requires discretization only on the boundaries of the modeling domain.  Solutions in the boundaries can be used to find solutions on the entire modeling domain, which makes the boundary element method a highly efficient numerical method. This paper will discuss 2D magnetotelluric (MT) modeling using boundary element method. The method is applied to calculate 2D MT responses, expressed in apparent resistivity and phase of impedance, for transverse electric (TE) mode.  The results obtained by the method are compared to those calculated by analytical solution and finite element modeling to show the accuracy of boundary element method. Keywords: Boundary-element method, 2D magnetotelluric.
KOREKSI STATIK DATA CSAMT (CONTROLLED SOURCE AUDIO-FREQUENCY MAGNETOTELLURIC) MENGGUNAKAN TEKNIK SPATIAL FILTERING Hamdi, Husnul; Srigutomo, Wahyu
JOURNAL ONLINE OF PHYSICS Vol 1, No 1 (2015): JoP
Publisher : Prodi Fisika FST UNJA

Show Abstract | Download Original | Original Source | Check in Google Scholar

Abstract

CSAMT adalah metoda elektromagnetik dengan domain frekuensi dengan menggunakan dipol yang ditanahkan sebagai sumber sinyal buatan. Pengukuran CSAMT dengan jarak sumber medan primer terhadap penerima yang berhingga menyebabkan gelombang yang ditangkap penerima bersifat kompleks. Pergeseran nilai medan listrik karena ada efek statik menyebabkan kurva resistivitas terangkat naik atau turun. sehingga mempengaruhi data hasil pengukuran. Tujuan penelitian ini untuk memperoleh data CSAMT yang sudah terkoreksi efek statik. Metode penelitian yang digunakan untuk menghilangkan efek statik digunakan teknik spatial filtering. Pada penelitian ini diperoleh data CSAMT yang sudah terkoreksi efek statik sehingga dapat diinterpretasikan dengan baik. Kata kunci : CSAMT; Efek Statik; Spatial Filtering.
Pemodelan Aliran Fluida 2-D Pada Kasus Aliran Permukaan Menggunakan Metode Beda Hingga Hapsoro, Cahyo Aji; Srigutomo, Wahyu
Jurnal Matematika dan Sains Vol 18 No 3 (2013)
Publisher : Institut Teknologi Bandung

Show Abstract | Download Original | Original Source | Check in Google Scholar

Abstract

Persamaan dasar aliran fluida yang disebut persamaan Navier-Stokes merupakan persamaan diferensial parsial non linier yang kompleks. Untuk menyelesaikan dan memodelkan aliran fluida perlu dilakukan pendekatan numerik, salah satunya dengan metode beda hingga. Penyelesaian persamaan Navier-Stokes dilakukan dengan meninjau beberapa asumsi penyederhanaan yaitu: fluida bersifat tak termampatkan, parameter aliran bergantung pada arah spasial x dan y, serta semua variabel dianggap sebagai fungsi periodik. Pemodelan numerik dilakukan untuk menghitung ketidakstabilan Kelvin-Helmholtz lapisan campuran, evolusi struktur vorteks dan dipol vorteks. Ketidakstabilan Kelvin-Helmholtz divariasikan dengan nilai panjang gelombang gangguan sebesar l = 0.5Lx  dan l = 0.25Lx. Bilangan Reynolds (Re) divariasikan dengan nilai 1000, 3000, dan 5000. Hasil pemodelan menunjukkan bahwa untuk ketiga Re tersebut, aliran fluida bersifat laminar, kritis, dan turbulen. Hal ini terindikasi oleh arah medan vortisitas serta distribusi massa jenis fluida. Semakin besar kecepatan aliran maka sifat aliran akan menjadi semakin acak atau turbulen. Kata kunci: Aliran fluida, Persamaan Navier-Stokes, Bilangan Reynolds, Vortisitas, Ketidakstabilan Kelvin-Helmholtz.   2-D Fluid Surface Flow Modeling using Finite-Difference Method Abstract Navier-Stokes equation is a complex non-linear second-order partial differential equation describing a fluid flow. Solving and modeling the fluid flow, a numerical method called finite difference method is frequently used. Several assumptions are incorporated in solving numerically Navier-Stokes equation: the fluid is incompressible, fluid flow parameters depend on its positions, and all variables are considered as periodic functions. In this paper numerical calculation has been carried out to model the Kelvin-Helmholtz instability of mixed layer, evolution of vortex structure and vortex dipole. The calculation is done by varying perturbation wavelength l = 0.5Lx  and l = 0.25Lx. Reynolds number (Re) is varied at 1000, 3000, and 5000. The results show that for the three values of Re, the properties of the flows are laminar, critical, and turbulent, respectively as indicated by the vorticity direction and distribution of fluid density. The larger value of fluid velocity, the more random and turbulent the fluid is. Keywords: Fluid flow, Navier-Stokes equation, Reynolds number, Vorticity, Kelvin-Helmholtz instability.
Analisis Parameter Anomali SP Menggunakan Metode Inversi Non linier Pendekatan Linier Fajriani, Fajriani; Hamzah, Hardi; Srigutomo, Wahyu; Pratomo, Prihamdhanu Mukti
Jurnal Saintifik Vol 5, No 1 (2019): Volume 5 Nomor 1
Publisher : Fakultas MIPA UNSULBAR

Show Abstract | Download Original | Original Source | Check in Google Scholar | DOI: 10.31605/saintifik.v5i1.193

Abstract

nterpretasi anomali self-potential (SP) model geometri sederhana (fixed geometry) sederhana, dilakukan untuk menghitung kedalaman, sudut polarisasi benda, moment dipole listrik, dan bentuk struktur benda dibawah permukaan yang mempengaruhi data SP yang diamati diatas permukaan. Dalam mengidentifikasi parameter anomali SP akibat dari model fixed geometry dibawah permukaan bumi dapat dilakukan dengan cara pemodelan inversi. Dalam tulisan ini pemodelan inversi diselesaikan menggunakan pendekatan linier dengan metode least-square terbobot dan metode Levenberg-Marquardt. Kedua metode ini diaplikasikan untuk menentukan nilai parameter dari anomali self-potential akibat struktur sesar, dengan menganggap struktur sesar menyerupai bentuk lempeng. Metode least-square terbobot dan Levenberg-Marquardt mampu mengikuti pola hasil data pengukuran di lapangan. Kedua metode inversi tersebut efisien karena proses perhitungannya yang relatif singkat serta nilai parameter model struktur lempeng hasil perhitungan kedua metode tersebut menunjukkan hasil yang sama hanya terdapat sedikit selisih perhitungan pada nilai kemiringan sudut lempeng.Kata kunci: anomali self-potential, fixed geometry, pendekatan linier
Perpindahan Panas pada Media Berpori Menggunakan Metode Elemen Hingga Ningsih, Surya; Srigutomo, Wahyu
Jurnal Saintifik Vol 4, No 2 (2018): volume 4 nomor 2 juli 2018
Publisher : Fakultas MIPA UNSULBAR

Show Abstract | Download Original | Original Source | Check in Google Scholar | DOI: 10.31605/saintifik.v4i2.182

Abstract

Masuknya aliran panas dari bawah reservoir ke dalam medium berpori di dalamnya mempengaruhi kondisi tekanan, entalpi, dan temperatur di dalam reservoir. Hal ini terjadi karena panas mengalir dari daerah bertemperatur tinggi ke daerah bertemperatur lebih rendah. Panas mengalir melalui daerah patahan (fracture zone) yang memiliki permeabilitas dan porositas lebih tinggi dari daerah sekitarnya di dalam reservoir untuk kemudian terjadi perpindahan panas secara konduksi menuju daerah sekitarnya. Tulisan ini mensimulasikan proses perpindahan panas dengan mengamati distribusi temperatur dan aliran fluida pada medium berpori. Metode numerik yang digunakan adalah metode elemen hingga. Hasil yang ditampilkan adalah distribusi tekanan, kecepatan aliran fluida, dan distribusi temperatur dan kemudian diinterpretasikan.Kata kunci: permeabilitas, porositas, kecepatan, distribusi temperatur, metode elemen hingga.
Three-dimensional DC Resistivity Modeling using Galerkin Finite Element Method Composed by Tetrahedral Elements Srigutomo, Wahyu; Anwar, Hairil; Agustine, Eleonora; Mahardika, Harry
Journal of Engineering and Technological Sciences Vol 51, No 4 (2019)
Publisher : ITB Journal Publisher, LPPM ITB

Show Abstract | Download Original | Original Source | Check in Google Scholar | Full PDF (768.912 KB) | DOI: 10.5614/j.eng.technol.sci.2019.51.4.5

Abstract

Successful interpretation of DC resistivity data depends on the availability of a proper forward modeling scheme. In this study, a three-dimensional DC resistivity forward modeling scheme was developed using the finite element method. The finite element equations were obtained using a weakened form of the weighted-residual method called the Galerkin method. Discretization of the modeling domain was carried out by dividing it into smaller three-dimensional blocks and subdividing each block into five tetrahedral elements. A linear interpolation function was employed and elemental linear equations were set up, followed by formation of global matrix systems of equation and incorporation of proper boundary conditions. The conjugate gradient method was applied to solve the global system of equations, which in this study was proven to be more efficient than a direct solver, contributing to a 67% time reduction. Using a Wenner array configuration, comparison with theoretical calculation of the electric potential for a homogeneous model yielded a relative error of 3.66%. To confirm the applicability of this forward modeling scheme, apparent resistivity profiles for several basic three-dimensional subsurface resistivity models were compared with the analytical profiles, yielding an acceptable level of fitting.
Three-dimensional DC Resistivity Modeling using Galerkin Finite Element Method Composed by Tetrahedral Elements Wahyu Srigutomo; Hairil Anwar; Eleonora Agustine; Harry Mahardika
Journal of Engineering and Technological Sciences Vol. 51 No. 4 (2019)
Publisher : Institute for Research and Community Services, Institut Teknologi Bandung

Show Abstract | Download Original | Original Source | Check in Google Scholar | DOI: 10.5614/j.eng.technol.sci.2019.51.4.5

Abstract

Successful interpretation of DC resistivity data depends on the availability of a proper forward modeling scheme. In this study, a three-dimensional DC resistivity forward modeling scheme was developed using the finite element method. The finite element equations were obtained using a weakened form of the weighted-residual method called the Galerkin method. Discretization of the modeling domain was carried out by dividing it into smaller three-dimensional blocks and subdividing each block into five tetrahedral elements. A linear interpolation function was employed and elemental linear equations were set up, followed by formation of global matrix systems of equation and incorporation of proper boundary conditions. The conjugate gradient method was applied to solve the global system of equations, which in this study was proven to be more efficient than a direct solver, contributing to a 67% time reduction. Using a Wenner array configuration, comparison with theoretical calculation of the electric potential for a homogeneous model yielded a relative error of 3.66%. To confirm the applicability of this forward modeling scheme, apparent resistivity profiles for several basic three-dimensional subsurface resistivity models were compared with the analytical profiles, yielding an acceptable level of fitting.
Global Inversion of Grounded Electric Source Time-domain Electromagnetic Data Using Particle Swarm Optimization Cahyo Aji Hapsoro; Wahyu Srigutomo; Acep Purqon; Warsa warsa; Doddy Sutarno; Tsuneomi Kagiyama
Journal of Engineering and Technological Sciences Vol. 53 No. 1 (2021)
Publisher : Institute for Research and Community Services, Institut Teknologi Bandung

Show Abstract | Download Original | Original Source | Check in Google Scholar | DOI: 10.5614/j.eng.technol.sci.2021.53.1.1

Abstract

Global optimization inversion of grounded wire time-domain electromagnetic (TDEM) data was implemented through application of the particle swarm optimization (PSO) algorithm. This probabilistic approach is an alternative to the widely used deterministic local-optimization approach. In the PSO algorithm, each particle that constitutes the swarm epitomizes a probable geophysical model comprised by subsurface resistivity values at several layers and layer thicknesses. The forward formulation of the TDEM problem for calculating the vertical component of the induced magnetic field is first expressed in the Laplace domain. Transformation of the magnetic field from the Laplace domain into the time domain is performed by applying the Gaver-Stehfest numerical method. The implementation of PSO inversion to the TDEM problem is straightforward. It only requires adjustment of a few inversion parameters such as inertia, acceleration coefficients and numbers of iteration and particles. The PSO inversion scheme was tested on synthetic noise-free data and noisy synthetic data as well as to field data recorded in a volcanic-geothermal area. The results suggest that the PSO inversion scheme can effectively solve the TDEM 1D stratified earth problem. 
Curie Point Depth Analysis of Lesugolo Area, East Nusa Tenggara, Indonesia Based on Ground Magnetic Data Alamta Singarimbun; Umar Said; Dini Andriani; R. B. Astro; Bakrun Bakrun; I G. P. F. Soerya Djaja; Eleonora Agustine; Pepen Supendi; Wahyu Srigutomo
Journal of Engineering and Technological Sciences Vol. 54 No. 1 (2022)
Publisher : Institute for Research and Community Services, Institut Teknologi Bandung

Show Abstract | Download Original | Original Source | Check in Google Scholar | DOI: 10.5614/j.eng.technol.sci.2022.54.1.11

Abstract

The Curie point depth, or magnetic basal depth, of the Lesugolo geothermal area in Ende, Flores Island, East Nusa Tenggara, Indonesia was estimated by performing spectral analysis on spatial magnetic data and transforming it into the frequency domain, resulting in a link between the 2D spectrum of magnetic anomalies and the depths of the top and centroid of the magnetic sources. Shallow Curie point depths of 16 to 18 km were found in the north-northeast to southeast areas of Lesugolo, while deeper depths of 24 to 26 km were found in the southwest. The tectonic setting beneath the central part of Flores Island governs the distribution of the Curie point depths in the area. Shallow Curie point depth zones are associated with high thermal gradients (30 to 34 °C/km) and heat flow (80 to 100 mW/m2). Deep depths, on the other hand, correspond to zones of low thermal gradient (21 to 26 °C/km) and low heat flow (65 to 80 mW/m2). Both the derived thermal gradient and the heat flow maps contribute to a better understanding of the Lesugolo geothermal system’s configuration. This study suggests that the Lesugolo geothermal area’s prospect zone is located in the center of the investigated area, where the Lesugolo normal fault forms its southeastern boundary.
Topography changes and thermal distribution at the Kelud crater after the 2014 Plinian eruption Wahyudi Wahyudi; Ari Setiawan; Heriansyah Putra; Herlan Darmawan; Imam Suyanto; Irwan Meilano; Irzaman irzaman; Maria Evita; Mitra Djamal; Moh Yasin; Nina Siti Aminah; Perdinan Perdinan; Retna Apsari; Wahyu Srigutomo; Wiwit Suryanto
Indonesian Journal of Geography Vol 52, No 3 (2020): Indonesian Journal of Geography
Publisher : Faculty of Geography, Universitas Gadjah Mada

Show Abstract | Download Original | Original Source | Check in Google Scholar | DOI: 10.22146/ijg.51986

Abstract

Topography of a volcano crater can change due to endogenic processes such as deformation or eruption, or surface processes. Erosion and deposition are surface processes that may occur and gradually change the slope of the inner volcano crater. Here, we investigated erosion and deposition processes that occurred in the Kelud crater after the 2014 plinian eruption. We used high-resolution Digital Elevation Models (DEMs) and orthomosaic images derived by drone photogrammetry that acquired between September 2018 and July 2019. We obtained hundreds of aerial images which were reconstructed to obtain 3D models of Kelud’s crater by using Structure from Motion (SfM) technique. Results show erosions at alluvial fans that dominantly located at the east valleys of Kelud crater. The erosion removed the volcanic materials up to -5 m which transported and deposited close to the vicinity of the Kelud crater. The deposition process causes the increase of the Kelud crater lake up to 3 m. Moreover, we also mapped the thermal distribution of the Kelud crater lake by using low cost thermal camera. Our thermal investigation is able to identify some hotspots at the vicinity of the Kelud crater lake with range temperature of 43.7°C – 55.3°C, while the average apparent temperature of the Kelud crater lake is ~ 29°C. This high temperature area may indicate underwater active fractures that continuously release volcanic gasses which leads to convection heat transfer through Kelud’s water lake.
Co-Authors Abdullah, Adi Rahmansyah Amir Acep Purqon Acep Purqon Adam Sukma Putra Adhitya S. Sunarya Adi R.A. Abdullah Alamta Singarimbun Alan Maulana Alan Maulana Alief Sadlie Kasman Anggie Susilawati Anwar, Hairil Ari Setiawan Arizal Akbar Zikri Asep Harja Asep Harja Asep Saefullah Azrul Azwar Bakrun Bakrun Candra E. Novana Deny Juanda Puradimaja Dini Andriani Dini Fitriani Doddy Sutarno Doddy Sutarno Doddy Sutarno Doddy Sutarno Doddy Sutarno Doddy Sutarno Eddy Supriyana Eleonora Agustine Eleonora Agustine Eleonora Agustine Eleonora Agustine Enjang Jaenal Mustopa Enjang Jaenal Mustopa Enjang Jaenal Mustopa Fajriani Fajriani Fajriani Fajriani Fariduzzaman Fariduzzaman Fourier Dzar Eljabbar Latief Frafti Rejeki Fungki Iqlima Nasyidiyah Ghany Hanifan Muslim Hairil Anwar Hairil Anwar Hairil Anwar Halawa Ernwati Hapsoro, Cahyo Aji Harapan Marpaung Hardi Hamzah Harry Mahardika Haryo Seno Hendra Gunawan Herlan Darmawan Herlin Tarigan Hisashi Utada Humolungo, Ismail Husnul Hamdi Husnul Hamdi, Husnul I G. P. F. Soerya Djaja I Gede Putu Fadjar Soerya Djaja Imam Suyanto Imran Hilman Mohammad Imran Hilman Mohammad Irwan Ary Darmawan irwan meilano Irzaman, Irzaman Ismail Humolungo Jayawarsa, A.A. Ketut Joko Sampurno Kartika Hajar Kirana Kumalasari, Ratih Kuncoro, Alvina Kusumadewi Linus Pasasa M. Habri Zen Mahardika, Harry Maria Evita Martin Liess Martin Liess Melany Febrina Mitra Djamal Moh Yasin Mohammad Heriyanto Mohammad Heriyanto Muhamad Hilmi Aufa Muhammad Abdilllah Budianto Nadhira Gunawan Ni Ketut Lasmia Nina Siti Aminah Pepen Supendi Pratomo, Prihamdhanu Mukti Prihamdhanu Mukti Pratomo Prihandhanu M. Pratomo Purnama, Dikdik Sidik Putra, Heriansyah Rahmi Elzufiah Rakotovao Lovanantenaina Omega Rasito Tursinah, Rasito RATIH KUMALASARI Ratih Kumalasari Retna Apsari Richardo Barry Astro Rifky Nauval Rostikwati, Diana Ayu Ruddy Kurnia Satria, Eko Seno, Haryo Sesri Santurima Setyanto Cahyo Pranoto Sidik, Adi Permana Sitti Balkis Sparisoma Viridi Sumintadireja, Prihadi Suprijadi Suprijadi Suprijadi Suprijadi Supriyadi Supriyadi Surya Ningsih Surya Ningsih Suryanto, Wiwit Tsuneomi Kagiyama Tsuneomi Kagiyama Tsuneomi Kagiyama Umar Fauzi UMAR FAUZI Umar Said Wahyu Hidayat Wahyudi Wahyudi Warsa warsa Wataru Kanda Yanti Yanti Yazid Bindar Yuant Tiandho Zulfahmi Zulfahmi Zulfahmi Zulfahmi Zulhendra Zulhendra