Australian Journal of Crop Science

Article | https://doi.org/10.21475/ajcs.26.20.07.pne168

Submitted: 03 November 2025 | Revised: 19 February 2026 | Accepted: 11 May 2026

Pages 550-555

Model design of white copra mechanical extrusion process to produce virgin coconut oil

Budy Rahmat*1, Yogi Nirwanto2, Zulfikar Noormasyah1, Tedi Hartoyo3, Iqbal Akmalludin4

1Doctoral Program of Agricultural Science, Postgraduate Program, Siliwangi University, Indonesia

2Department of Agrotechnology, Agriculture Faculty, Siliwangi University, Indonesia

3Department of Agrobusiness, Agriculture Faculty, Siliwangi University, Indonesia

4Department of Agrotechnology, Postgraduate Program, Siliwangi University, Indonesia

Adress: Siliwangi Street 24th, Tasikmalaya City-46115, West Java, Indonesia

Corresponding author: Budy Rahmat budyrahmat@unsil.ac.id

ORCHID ID: 0000-0002-2646-6445

Abstract: This study aims to design a mechanical extrusion process model - namely the process of extracting oil from vegetable materials using a screw press - in the production of virgin coconut oil (VCO) from white copra. White copra is coconut meat that preserved by a proper drying process to avoid microbial infection during storage. The focus of this study is to determine the relationship between the physical process parameters of extrusion and the yield and quality of the resulting VCO. This study addresses the limitations of conventional VCO extraction methods such as fermentation, enzymatic, and centrifugation processes, which although produce high-quality VCO, but have shortcomings in time and cost efficiency for industrial scale. The experiment was arranged using a Completely Randomized Design (CRD) with a factorial pattern consisting of three levels of copra moisture content (C1 = 6%, C2 = 8%, and C3 = 10%), two extrusion temperatures (T1 = 60°C and T2 = 80°C), and two screw speeds (S1 = 30 rpm and S2 = 60 rpm). The combination of these three factors formed 12 treatments and was carried out three replications, thus requiring 36 experimental batches. 100 g of grated white copra was provided for each batch of extrusion experiments. The observed VCO yield parameters included: oil yield, free fatty acid content, and peroxide value. The results of the study showed that the treatment factors of water content (C), extrusion temperature (T), and screw speed (S), significantly affected the oil yield and quality. The highest VCO yield (64.47%) was obtained in the treatment of 6% water copra content, 80oC process temperature and 60 ppm screw speed (C1T2S2), with a regression model Y1 = 54.82 − 0.65C + 2.15T + 3.24S. In terms of VCO quality, it turned out that the lowest FFA content (0.254%) occurred in the C1T1S1 treatment, while the highest (0.581%) occurred in the C3T2S2, explained by the model Y2 = −0.050 + 0.035C + 0.143T + 0.119S. Peroxide values ​​increased with increasing temperature and screw rotation speed, following the model Y3 = 1.623 + 0.078C + 0.065T + 0.510S. Overall, increasing temperature and screw speed increased the oil extraction efficiency but reduced the chemical stability of VCO due to increased lipid oxidation and hydrolysis. The resulting regression model effectively describes the empirical relationship between extrusion process parameters and oil yield, allowing optimization of extrusion conditions to produce high-yield and high-quality VCO.

Keywords: free fatty acids, mechanical extrusion, process modeling, virgin coconut oil, white copra.

Abbreviations: VCO_Virgin coconut oil; CRD_Completely randomized design; FFA_Free fatty acid; SNI_Standar Nasional Indonesia; PV_Peroxide value; DMRT_Duncan's Multiple Differences Test.

Introduction

The centrifugation and freeze-thaw methods as initial mechanical approaches have demonstrated superior quality, with a moisture content of 0.12%, free fatty acid (FFA) of 0.058%, and an adequate yield of 31.42%. However, these methods require a processing time of more than 48 hours for freezing and operate in batch mode (Pulungan et al., 2020). Centrifugation demands a large capital investment for high-speed equipment (Jo et al., 2024). This method is less optimal in meeting industrial needs for scalability and cost efficiency.

Vegetable oil extraction using screw extruders has become the industry standard due to its continuous processing capability and high yield. In coconut oil pressing research, the highest yield achieved was 63.62%, attained with high thermal application between 100 °C and 125 °C (Zikri et al., 2020). Although coconut oil produced at 125 °C technically meets SNI 01-2902:1992 standards with a moisture content of 0.41% and free fatty acid (FFA) of 0.32%, applying heat at this level violates the virgin coconut oil (VCO) principle, which requires a non-thermal process (Indahsari and Ramadhani, 2024).

The technical gap above shifts the research focus from equipment application to the design and modeling stage of process engineering. In extrusion processes, mechanical energy applied through screw rotation speed and compression pressure inevitably converts into thermal energy via friction, which can raise the internal temperature of the material. The model must be able to predict the frictional heat generation rate based on screw speed and pressure.

Based on the above explanation, it is necessary to design and model a mechanical extrusion process that can link physical parameters with an optimal VCO oil yield, thus ensuring the production of high-quality VCO. The objective of this research is to develop innovative process design and modeling for VCO production using a mechanical extrusion system. The scientific contribution of this study is expected to fill the knowledge gap in extrusion modeling for thermally sensitive emulsion materials.

Results and Discussion

Oil yield

Analysis of variance shows that the combination of moisture content (C), temperature (T), and screw speed (S) significantly affects the oil yield obtained from the extrusion process of white copra. As shown in Table 1, the highest the VCO yield was achieved under treatment C1T2S2 (low moisture content, high temperature, high speed), at 64.473%. This is due to increased oil flow efficiency caused by decreased material viscosity at higher extrusion temperatures and speeds (Sulhatun et al., 2020). Meanwhile, treatment C3T1S1 produced the lowest yield of 57.213%.

Low moisture content facilitates easier oil release during copra pressing (Mizera et al., 2023). Similarly, higher temperature, both from preheating and screw friction, enhances oil flow (Zikri et al., 2020). Meanwhile, high screw speed applies maximum pressing and friction forces required to efficiently extract the prepared oil from copra pulp, resulting in the highest VCO yield (Abdilova et al., 2023).

Regression analysis for the independent variables—copra moisture content (C), extrusion temperature (T), and extrusion speed (S)—on the VCO yield (Y1) produced the following equation:

Y1 =54.82−0.65C+2.15T+3.24S

This regression model indicates oil yield increases with temperature and screw speed but slightly decreases with higher moisture content. The screw speed factor (S) has the greatest influence on increasing oil yield.

Simultaneously, the treatments of moisture content, temperature, and speed affect copra oil yield. Increasing temperature and speed improve oil yield (Zikri et al., 2020; Abdilova et al., 2023), whereas increasing moisture content reduces it (Mizera et al., 2023). The mechanical extrusion of white copra at low moisture, higher temperature, and high screw speed results in the highest VCO yield because these factors synergistically optimize oil separation from the copra solid material.

Free fatty acid (FFA) content

Based on Duncan's further test results (Tabel 1), the lowest oil FFA content was found in treatment C1T1S1 (0.254%), while the highest was observed in C3T2S2 (0.581%). Generally, increases in material moisture content (C), extrusion temperature (T), and screw speed (S) tend to raise the oil’s FFA content. There is a trend where combinations of high moisture, high temperature, and high screw speed produce the highest FFA compared to other combinations. FFA content increases with rising temperature (Ng et al., 2021), moisture (Divya et al., 2023), and extrusion speed (Feng et al., 2014). This confirms that extreme extrusion conditions accelerate oxidation and hydrolysis of the oil.

The combination of high copra moisture (reactant), high temperature (reaction accelerator), and high screw speed (heat generator and reactant releaser) creates optimal conditions for hydrolysis reactions, directly increasing FFA levels in VCO. High FFA content is the main cause of rancid odor and taste in the oil (Song et al., 2022).

The regression relationship between copra moisture content (C), extrusion temperature (T), screw speed (S), and FFA content in VCO (Y2) is expressed as:

Y2 =−0.050+0.035C+0.143T+0.119S

This regression model shows that FFA content increases with temperature and screw speed but slightly decreases as moisture content rises. The screw speed factor (S) has the largest effect on oil yield increase. The model indicates all factors tend to increase FFA, with temperature (T) having the most significant effect. Higher temperature and screw speed accelerate fat hydrolysis, thus increasing free fatty acid content. Simultaneously, moisture, temperature, and speed treatments influence copra oil FFA levels. Increasing moisture, temperature, and speed will raise the copra oil FFA content (Ng et al., 2021; Divya et al., 2023; Feng et al., 2014).

Increasing moisture, temperature, and extrusion speed during mechanical extrusion of white copra can increase the FFA content in the produced VCO. This FFA increase indicates a decline in oil quality because these three factors stimulate and accelerate hydrolysis and thermal degradation reactions.

Peroxide value

Analysis of variance results indicate that all three tested factors—material moisture content (C), extrusion temperature (T), and screw speed (S)—affect the peroxide value of the resulting oil. The screw speed factor has the greatest effect on increasing peroxide value, followed by extrusion temperature, while moisture content shows an effect but with lower intensity (Table 1).

Peroxide value is an early indicator of fat oxidation. An increasing value indicates the formation of peroxide compounds from the reaction between oxygen and unsaturated fatty acids. Higher temperature and extrusion speed accelerate oxidative reactions, causing peroxide values to rise. PV significantly increase at high moisture, temperature, and speed. The highest value of 3.361% was recorded for treatment C3T2S2, indicating enhanced fat oxidation due to high temperature and friction during extrusion. These three factors create conditions that maximize oxidation rate (high temperature and friction) and

Table 1. The effect of copra water content, temperature, and extrusion speed of white copra on the yield, FFA content, and peroxide value of oil.

Treatments

Y1 = Oil yield

(%)

Y2 = FFA content

(% b/b)

Y3 = Peroxide value

(meq O₂/kg)

C1T1S1 59.570 bc 0.254 a 2.199 a
C1T1S2 63.477 efg 0.393 bc 3.272 cd
C1T2S1 61.143 cd 0.398 bc 2.875 b
C1T2S2 64.473 g 0.473 cd 3.373 d
C2T1S1 58.473 ab 0.281 A 2.268 a
C2T1S2 62.817 def 0.402 bc 3.231 cd
C2T2S1 60.973 cd 0.385 bc 2.919 bc
C2T2S2 63.920 fg 0.513 d 3.332 d
C3T1S1 57.213 a 0.294 a 2.171 a
C3T1S2 61.637 de 0.436 cd 3.090 bc
C3T2S1 60.563 cd 0.509 d 2.937 b
C3T2S2 64.047 fg 0.254 e 3.361 d

Note: Numbers followed by the same letter in each column indicate no difference according to Duncan's Multiple Range Test (DMRT).

Table 2. Combination of factors and treatment levels formed.

Water content of copra Extrusion temperature Extrusion speed
S1 (30 rpm) S2 (60 rpm)
C1 (6%) T1 (60°C) C1T1S1 C1T1S2
T2 (80°C) C1T2S1 C1T2S2
C2 (8%) T1 (60°C) C2T1S1 C2T1S2
T2 (80°C) C2T2S1 C2T2S2
C3 (10%) T1 (60°C) C3T1S1 C3T1S2
T2 (80°C) C3T2S1 C3T2S2

Table 3. Analysis of variance of experimental data.

Source of Variance Sum of Square (SS) Degree of Freedom (df) Mean of Square (MS) = SS/df FCal FTable
C SSC a -1 MSC FC
T SST b -1 MST FT
S SSS c -1 MSS FS
Error SSE dfE MSE - -
Total SSTotal N -1 - - -

cause vulnerable substrates to oxidize due to high moisture and FFA levels, all leading to increased peroxide compounds and decreased VCO quality (Chandran et al., 2017; Agung and Rismaya, 2024; Harimurti et al., 2024). Increasing moisture, temperature, and extrusion speed during mechanical extrusion of white copra can raise the peroxide number in VCO. An increase in PV is a key indicator of oil quality degradation due to primary oxidation.

Regression analysis for the effects of copra moisture content (C), extrusion temperature (T), and extrusion speed (S) on PV in VCO (Y3) is expressed as:

Y3 =1.623+0.078C+0.065T+0.510S

The study shows that physical process factors (moisture content, temperature, and extrusion speed) significantly affect all quality parameters of oil extracted from white copra. Thus, there is a compromise between extraction efficiency (high yield) and chemical stability of the oil (low FFA and PV). Optimal processing conditions can be determined by balancing these two aspects. Oil yield increases at high temperature and screw speed because viscosity decreases and extrusion pressure rises. However, FFA and PV also increase under these conditions due to more intense oxidation and hydrolysis processes.

Materials and Methods

Materials and equipment

White copra is coconut meat that preserved by a proper drying process to avoid microbial infection during storage as main material, dried to the appropriate moisture content: 6%, 8%, and 10%. This white copra was made from ripe coconut flesh from a tall variety obtained from smallholder coconut plantations on Cikalong Beach, Tasikmalaya Regency.

Figure 1. VCO extrusion process from white copra to determine yield, FFA content, and PV.

The Ossel Type OS-D035 oil extrusion machine for laboratory scale is equipped with a 452-watt Ceramic Band United Heater with automatic temperature control, a titration unit, and a 200-watt SG coconut grater. In addition, supporting equipment is available, such as analytical balances, graduated pipettes, beakers, etc.

Experimental design

The experiment consisted of three treatment factors, each with its own level: C: copra moisture content (6%, 8%, and 10%); T: extrusion temperature (60°C and 80°C); and S: screw speed (30 rpm and 60 rpm). The treatments were arranged in a 3×2×2 Completely Randomized Design, resulting in 12 treatment combinations, as shown in Table 2.

Each treatment combination was repeated three times, so that 36 experimental units were provided.

Experimental procedure

Three levels of whole white copra moisture were each grated separately to facilitate the extrusion process to extract virgin coconut oil. Each experimental unit was prepared with 100 g of grated white copra for each treatment. Each unit was loaded into the extruder with the temperature and speed settings specified for the treatment combination mentioned above. The VCO from the copra extruded was then weighed to calculate its yield and placed in labeled bottles for each treatment. The VCO in these bottles was prepared for analysis of FFA content and PV (Figure 1).

Observation parameters

1) Oil yield

The VCO yield produced from every 100 g of white copra is calculated using the following formula (Anwar and Salima, 2016):

\[Oil\ yield\ \ (\%) = \frac{Mass\ of\ VCO}{Mass\ of\ copra} \times 100\%\]

2) Free fatty acids

Determination of the free fatty acid (FFA) content in VCO is performed using the acid-base titration method, which neutralizes the free fatty acids present in the VCO with a standard base solution of NaOH (Uçar et al., 2024). The general formula for calculating the free fatty acid content in percentage is:

\[FFA\ Content\ (\%) = \frac{V \times N \times MW \times 100}{S}\]

Where: V = volume of NaOH solution (standard) used (in mL); N = normality of NaOH solution; MW = molecular weight of fatty acid as standard fatty acid (e.g. oleic acid = 282 g/mol); S = weight of oil sample (in mg or g, adjusted units)

The FFA content indicates the percentage of free fatty acids in total oil and is an indicator of oil quality, where high FFA levels indicate triglyceride hydrolysis and decreased oil quality (Anwar and Salima, 2016; Tan et al., 2023)

3) Peroxide value

Peroxide value (meq O2/kg) of VCO is determined using the following formula (Oktaviana et al., 2025):

\[{Peroxide\ Value\ (meq\ O}_{2}/kg) = \frac{V \times N \times 1000}{W}\]

Where: V: Volume of sodium thiosulfate (Na2S2O3) solution used in the titration (in milliliters, mL); N: Normality of the sodium thiosulfate (usually 0.01 N or according to the titration standard used); 1000: Conversion factor to convert grams of sample to kilograms; W: Weight of the oil sample being titrated (in grams, g).

This formula calculates the number of equivalents of peroxide (active oxygen) in 1 kg of oil based on the volume of titrant required to react with the iodine formed from the peroxide oxidation process in the sample (Laga et al., 2023).

Data analysis

1) Analysis of variance

The design used was a completely randomized design (CRD) with three factors, an additive factorial pattern without interaction. Factor C had 3 levels, T had 2 levels, S had 2 levels, and each combination had 3 replications. The total number of observations was N = 3 × 2 × 2 × 3 = 36 experimental units.

Linear model of observation data after treatment

Yijkl=μ+Ci+Tj+Skijkl

with: μ = general mean; Ci = effect of factor C at level i (i = 1..a); Tj = effect of factor T at level j (j = 1..b); Sk = effect of factor S at level k (k = 1..c); and εijkm ~ IID N(0,σ2) (error).

Based on the linear model, the analysis of variance is compiled as shown in Table 3.

Hypothesis formulation:

Null hypothesis (H0): All treatments have the same mean.

Alternative hypothesis (H1): At least one treatment has a different mean.

Hypothesis testing:

Compare the calculated F value with the table F value at a specific significance level (e.g., 0.05).

If the calculated F value > table F value, reject the null hypothesis, indicating a significant difference between the treatments.

If a significant difference is found, conduct a Duncan's Multiple Differences Test (DMRT) to determine which treatment is different (Huzsvai et al., 2022). This method can be implemented using statistical software such as SPSS or Excel to facilitate calculations and significance testing (Assaad et al., 2015).

2) Regression analysis

Multiple Linear Regression. Used when more than one extrusion process factor simultaneously influences oil yield. The multiple linear regression model is as follows:

\[Y = a + b_{1}X_{1} + b_{2}X_{2} + \ldots + b_{n}X_{n}\]

Where X1, X2,…, Xn are the process factors being tested.

This analysis provides information on the contribution of each factor to oil yield and its statistical significance (Irvan et al., 2025).

Model Validation Analysis. Using this regression method, a mathematical model of the relationship between extrusion process factors and oil yield can be obtained, allowing for empirical optimization based on experimental data (Espínola et al., 2011; Mursalykova et al., 2023).

This reference provides a comprehensive explanation of various regression techniques frequently used in experimental data analysis, such as the influence of extrusion process factors on oil yield.

Conclusion

This regression model was designed to show that oil yield increases with increasing temperature and screw speed, but decreases slightly with increasing moisture content. The regression model also shows that FFA content increases with increasing temperature and screw speed, but decreases slightly with increasing moisture content. Factor V (screw speed) has the greatest influence on increasing oil yield. Increasing moisture content, temperature, and extrusion speed in the mechanical extrusion of white copra can increase the peroxide value (PV) in VCO. This increase is a primary indicator of oil quality deterioration due to primary oxidation. The results showed that the combination of moisture content (C), extrusion temperature (T), and screw speed (S) significantly affected the yield, free fatty acid (FFA) content, and PV of oil produced from the mechanical extrusion process of white copra.

Author Contributions

BR: project management, experimental design and manuscript revision. YN and IA: experiment conduction, data analysis, and manuscript drafting. ZN and TH: experimental design and data analysis.

Declaration of Competing Interest

The authors declare no conflicts of interest regarding the publication of this paper.

Acknowledgments

The authors extend their sincere gratitude to the Director of Postgraduate Studies and the Dean of the Faculty of Agriculture for their human resource support, as well as to the Head of the Research and Community Service Institute of Siliwangi University for the funding support of the Professor Research Scheme.

References

Abdilova G, Sergibayeva Z, Orynbekov D, Shamenov M, Zhumadilova G, Bakiyeva A, Mukashev N, Bayadilova A, Dukenbayev D (2023). Influence of grinding degree and screw rotation speed on sunflower oil pressing process. Appl. Sci. 13:9958. https://doi.org/10.3390/app13179958

Agung GS, Rismaya R (2024). Pengaruh suhu pemanasan terhadap karakteristik mutu minyak goreng bekas pakai pedagang gorengan. Agritekno J. Teknol. Pertan. 13:15–23. https://doi.org/10.30598/jagritekno.2024.13.1.15

Anwar C, Salima R (2016) Perubahan rendemen dan mutu virgin coconut oil (VCO) pada berbagai kecepatan putar dan lama waktu sentrifugasi. J. Teknotan. 10:51–60. https://doi.org/10.24198/jt.vol10n2.8

Assaad HI, Hou Y, Zhou L, Carroll RJ, Wu G (2015). Rapid publication-ready MS-Word tables for two-way ANOVA. SpringerPlus. 4:33. https://doi.org/10.1186/s40064-015-0795-z

Chandran J, Nayana N, Roshini N, Nisha P (2017) Oxidative stability, thermal stability and acceptability of coconut oil flavored with essential oils from black pepper and ginger. J. Food Sci. Technol. 54:144–152. https://doi.org/10.1007/s13197-016-2446-y

Divya PM, Roopa BS, Manusha C, Balannara P (2023). A concise review on oil extraction methods, nutritional and therapeutic role of coconut products. J. Food Sci. Technol. 60:441–452. https://doi.org/10.1007/s13197-022-05352-0

Espínola F, Moya M, Fernández DG, Castro E (2011). Modelling of virgin olive oil extraction using response surface methodology. Int. J. Food Sci. Technol. 46:2576–2583. https://doi.org/10.1111/j.1365-2621.2011.02786.x

Feng H, Sui X, Chang Y, Qi B, Zhang Y, Li Y, Jiang L (2014) Effect of extruding full-fat soy flakes on trans fat content. Sci. World J. 2014:427423. https://doi.org/10.1155/2014/427423

Harimurti S, Prayogo GK, Widada H, Orbayinah S, Sukamdi DP, Erviana R, Amid A, Tho DCMV (2024). Degradation study of bromelain-extracted virgin coconut oil based on peroxide and fatty acid formation under elevated temperatures. BIO Web Conf. 144:05003. https://doi.org/10.1051/bioconf/202414405003

Huzsvai L, Mohammed S, Harsányi E, Széles A (2022). Novel approach for statistical interpretation: A case study from long-term crop production experiments (Hungary). Horticulturae. 8:48. https://doi.org/10.3390/horticulturae8010048

Indahsari LEN, Ramadhani AN (2024). Quality of virgin coconut oil production using pineapple flesh and skin based on chemical parameters. Berk. Penelit. Hayati. 30:1–4. https://doi.org/10.23869/bphjbr.30.2.20246

Irvan I, Zahedi Z, Nusantara BC, Romantica KP, Siswadi S, Nasution MD (2025). Influence of world development indicators on human development: A multiple linear regression and PSI ranking approach. ZERO J. Sci. Math. Appl. 9:77. https://doi.org/10.30829/zero.v9i1.23925

Jo P, Dewi AAK, Dika IW (2024). Rancang bangun mesin pemurnian virgin coconut oil (VCO) dengan metode sentrifugal berbasis internet of things (IoT). J. Manaj. Teknol. Inf. 15:15–22. https://doi.org/10.59819/jmti.v15i2.3986

Laga A, Bastian F, Musdalifa M, Muhpidah M, Djalal M (2023). Quality evaluation of virgin coconut oil produced with enzymatic extraction using coated crude papain. Curr. Res. Nutr. Food Sci. J. 11:174–186. https://doi.org/10.12944/CRNFSJ.11.1.12

Mizera Č, Aleš Z, Herák D, Hrabě P, Kabutey A, Napitupulu RAM, Ungureanu N (2023). Mechanical pressing of coconut oil and evaluation of its lubricant properties. Processes. 11:3034. https://doi.org/10.3390/pr11103034

Mursalykova M, Kakimov M, Kassenov A, Iskakov B, Sergibayeva Z, Kaspakov E, Zhumadilova G, Shulenova A, Kokayeva G, Suychinov A (2023). Mathematical modeling of screw press configuration for processing safflower oil. Appl. Sci. 13:3057. https://doi.org/10.3390/app13053057

Ng YJ, Tham PE, Khoo KS, Cheng CK, Chew KW, Show PL (2021). A comprehensive review on the techniques for coconut oil extraction and its application. Bioprocess Biosyst. Eng. 44:1807–1818. https://doi.org/10.1007/s00449-021-02577-9

Oktaviana L, Nugroho WA, Al Riza DF (2025). Evaluation of peroxide value and free fatty acid content in coconut oil and palm oil under different heating temperature treatments using reflectance–fluorescence-based computer vision. J. Pangan Agroind. 13:101–111. https://doi.org/10.21776/ub.jpa.2025.013.02.4

Pulungan MH, Fadhlillah ON, Dewi IA (2020). Optimasi proses pembekuan minyak kelapa murni (virgin coconut oil) menggunakan metode RSM: Kajian suhu dan lama waktu pembekuan. J. Pangan Agroind. 8:1–12. https://doi.org/10.21776/ub.jpa.2020.008.01.1

Sulhatun S, Mutiawati M, Kurniawan E (2020). Pengaruh temperatur dan waktu pemasakan terhadap perolehan minyak kemiri dengan menggunakan cara basah. J. Teknol. Kim. Unimal. 9:54–60. https://doi.org/10.29103/jtku.v9i2.4400

Tan BA, Nair A, Zakaria MIS, Low JYS, Kua SF, Koo KL, Wong YC, Neoh BK, Lim CM, Appleton DR (2023). Free fatty acid formation points in palm oil processing and the impact on oil quality. Agriculture. 13:957. https://doi.org/10.3390/agriculture13050957

Uçar B, Gholami Z, Svobodová K, Hradecká I, Hönig V (2024). A comprehensive study for determination of free fatty acids in selected biological materials: A review. Foods. 13:1891. https://doi.org/10.3390/foods13121891

Zikri A, Aswan A, Erlinawati E, Fatria F, Pratama Y, Anggraini T, Cendikia MB (2020). Uji kinerja screw oil press machine ditinjau dari rendemen dan kualitas minyak kelapa yang dihasilkan. Fluida. 13:46–53. https://doi.org/10.35313/fluida.v13i2.2305