Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Validated Models”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 847 records · Page 47

Consistent Evaluation of ACOS-GOSAT, BESD-SCIAMACHY, CarbonTracker, and MACC Through Comparisons to TCCON

Consistent validation of satellite CO2 estimates is a prerequisite for using multiple satellite CO2 measurements for joint flux inversion, and for establishing an accurate long-term atmospheric CO2 data record. Harmonizing satellite CO2 measurements is particularly important since the differences in instruments, observing geometries, sampling strategies, etc. imbue different measurement characteristics in the various satellite CO2 data products. We focus on validating model and satellite observation attributes that impact flux estimates and CO2 assimilation, including accurate error estimates, correlated and random errors, overall biases, biases by season and latitude, the impact of coincidence criteria, validation of seasonal cycle phase and amplitude, yearly growth, and daily variability. We evaluate dry-air mole fraction (X(sub CO2)) for Greenhouse gases Observing SATellite (GOSAT) (Atmospheric CO2 Observations from Space, ACOS b3.5) and SCanning Imaging Absorption spectroMeter for Atmospheric CHartographY (SCIAMACHY) (Bremen Optimal Estimation DOAS, BESD v2.00.08) as well as the CarbonTracker (CT2013b) simulated CO2 mole fraction fields and the Monitoring Atmospheric Composition and Climate (MACC) CO2 inversion system (v13.1) and compare these to Total Carbon Column Observing Network (TCCON) observations (GGG2012/2014). We find standard deviations of 0.9, 0.9, 1.7, and 2.1 parts per million vs. TCCON for CT2013b, MACC, GOSAT, and SCIAMACHY, respectively, with the single observation errors 1.9 and 0.9 times the predicted errors for GOSAT and SCIAMACHY, respectively. We quantify how satellite error drops with data averaging by interpreting according to (error(sup 2) equals a(sup 2) plus b(sup 2) divided by n (with n being the number of observations averaged, a the systematic (correlated) errors, and b the random (uncorrelated) errors). a and b are estimated by satellites, coincidence criteria, and hemisphere. Biases at individual stations have year-to-year variability of 0.3 parts per million, with biases larger than the TCCON predicted bias uncertainty of 0.4 parts per million at many stations. We find that GOSAT and CT2013b under-predict the seasonal cycle amplitude in the Northern Hemisphere (NH) between 46 and 53 degrees North latitude, MACC over-predicts between 26 and 37 degrees North latitude, and CT2013b under-predicts the seasonal cycle amplitude in the Southern Hemisphere (SH). The seasonal cycle phase indicates whether a data set or model lags another data set in time. We find that the GOSAT measurements improve the seasonal cycle phase substantially over the prior while SCIAMACHY measurements improve the phase significantly for just two of seven sites. The models reproduce the measured seasonal cycle phase well except for at Lauder_125HR (CT2013b) and Darwin (MACC). We compare the variability within 1 day between TCCON and models in June-July-August; there is correlation between 0.2 and 0.8 in the NH, with models showing 10-50 percent the variability of TCCON at different stations and CT2013b showing more variability than MACC. This paper highlights findings that provide inputs to estimate flux errors in model assimilations, and places where models and satellites need further investigation, e.g., the SH for models and 45-67 degrees North latitude for GOSAT and CT2013b.

ACOS-GOSAT↗

The Impact of High-Resolution Sea Surface Temperatures on the Simulated Nocturnal Florida Marine Boundary Layer

High- and low-resolution sea surface temperature (SST) analysis products are used to initialize the Weather Research and Forecasting (WRF) Model for May 2004 for short-term forecasts over Florida and surrounding waters. Initial and boundary conditions for the simulations were provided by a combination of observations, large-scale model output, and analysis products. The impact of using a 1-km Moderate Resolution Imaging Spectroradiometer (MODIS) SST composite on subsequent evolution of the marine atmospheric boundary layer (MABL) is assessed through simulation comparisons and limited validation. Model results are presented for individual simulations, as well as for aggregates of easterly- and westerly-dominated low-level flows. The simulation comparisons show that the use of MODIS SST composites results in enhanced convergence zones. earlier and more intense horizontal convective rolls. and an increase in precipitation as well as a change in precipitation location. Validation of 10-m winds with buoys shows a slight improvement in wind speed. The most significant results of this study are that 1) vertical wind stress divergence and pressure gradient accelerations across the Florida Current region vary in importance as a function of flow direction and stability and 2) the warmer Florida Current in the MODIS product transports heat vertically and downwind of this heat source, modifying the thermal structure and the MABL wind field primarily through pressure gradient adjustments.

LaCasse, Katherine M.↗

Learning continuous models for continuous physics

Abstract Dynamical systems that evolve continuously over time are ubiquitous throughout science and engineering. Machine learning (ML) provides data-driven approaches to model and predict the dynamics of such systems. A core issue with this approach is that ML models are typically trained on discrete data, using ML methodologies that are not aware of underlying continuity properties. This results in models that often do not capture any underlying continuous dynamics—either of the system of interest, or indeed of any related system. To address this challenge, we develop a convergence test based on numerical analysis theory. Our test verifies whether a model has learned a function that accurately approximates an underlying continuous dynamics. Models that fail this test fail to capture relevant dynamics, rendering them of limited utility for many scientific prediction tasks; while models that pass this test enable both better interpolation and better extrapolation in multiple ways. Our results illustrate how principled numerical analysis methods can be coupled with existing ML training/testing methodologies to validate models for science and engineering applications.

97 MATHEMATICS AND COMPUTING↗

Investigation of arterial gas occlusions

The effect of noncondensable gases on high-performance arterial heat pipes was investigated both analytically and experimentally. Models have been generated which characterize the dissolution of gases in condensate, and the diffusional loss of dissolved gases from condensate in arterial flow. These processes, and others, were used to postulate stability criteria for arterial heat pipes under isothermal and non-isothermal condensate flow conditions. A rigorous second-order gas-loaded heat pipe model, incorporating axial conduction and one-dimensional vapor transport, was produced and used for thermal and gas studies. A Freon-22 (CHCIF2) heat pipe was used with helium and xenon to validate modeling. With helium, experimental data compared well with theory. Unusual gas-control effects with xenon were attributed to high solubility.

Saaski, E. W.↗

A thermal-capillary mechanism for a growth rate limit in edge-defined film-fed growth of silicon sheets

Capillarity, acting to set the shape of the melt/gas interfaces, and heat transfer can interact to cause limits to steady-state growth of thin silicon sheets by the Edge-Defined Film-Fed Growth (EFG) method. A finite-element/Newton solution method for a two-dimensional thermal-capillary model of EFG is used to show that limiting values of pull rate exist beyond which steady-state growth is impossible. The pull rate limit is also predicted by a one-dimensional heat transfer model valid when the die sides and menisci are almost parallel and when the thermal conductivities of melt, crystal, and die are all equal. Both the one- and two-dimensional heat transfer models show that heat loss from the melt is dominated by conduction into the crystal and slow heat release to the ambient along the length of the ribbon. The limiting pull rate results from the reduced efficiency of conduction through the melt caused by the curvature of the meniscus which increases height of the die top above the level of the melt. Thermal-capillary limits are predicted for both positive and negative pressure differences across the meniscus.

Thomas, P. D.↗

The Contribution of CEOP Data to the Understanding and Modeling of Monsoon Systems

CEOP has contributed and will continue to provide integrated data sets from diverse platforms for better understanding of the water and energy cycles, and for validaintg models. In this talk, I will show examples of how CEOP has contributed to the formulation of a strategy for the study of the monsoon as a system. The CEOP data concept has led to the development of the CEOP Inter-Monsoon Studies (CIMS), which focuses on the identification of model bias, and improvement of model physics such as the diurnal and annual cycles. A multi-model validation project focusing on diurnal variability of the East Asian monsoon, and using CEOP reference site data, as well as CEOP integrated satellite data is now ongoing. Preliminary studies show that climate models have difficulties in simulating the diurnal signals of total rainfall, rainfall intensity and frequency of occurrence, which have different peak hours, depending on locations. Further more model diurnal cycle of rainfall in monsoon regions tend to lead the observed by about 2-3 hours. These model bias offer insight into lack of, or poor representation of, key components of the convective and stratiform rainfall. The CEOP data also stimulated studies to compare and contrasts monsoon variability in different parts of the world. It was found that seasonal wind reversal, orographic effects, monsoon depressions, meso-scale convective complexes, SST and land surface land influences are common features in all monsoon regions. Strong intraseasonal variability is present in all monsoon regions. While there is a clear demarcation of onset, breaks and withdrawal in the Asian and Australian monsoon region associated with climatological intraseasonal variabillity, it is less clear in the American and Africa monsoon regions. The examination of satellite and reference site data in monsoon has led to preliminary model experiments to study the impact of aerosol on monsoon variability. I will show examples of how the study of the dynamics of aerosol-water cycle interactions in the monsoon region, can be best achieved using the CEOP data and modeling strategy.

Lau, William K. M.↗

A novel machine learning based identification of potential adopter of rooftop solar photovoltaics

With the proliferation of rooftop solar photovoltaic installations, there is a need to proactively predict consumer potential for solar photovoltaic adoption, for improved electric utility planning and operation. Traditional analytical modeling approaches are limited to a few survey features and a larger part of the survey would remain untouched by the decision model. This article presents a novel, data-driven modeling approach that strategically prunes a large set of consumer profile features using a machine learning framework to train a model for predicting potential solar adoption. The approach utilizes the Gradient Boosting Decision Tree model through a Light Gradient Boosting framework that improves significantly over the poor prediction accuracy of the existing approaches. Model training using focal-loss based supervision is used to overcome the difficulty in identifying the potential adopters that is inherent in conventional data-driven models. In addition, to overcome possible data sparsity in a limited survey sample, a Generative Adversarial Network is presented to create synthetic user samples and its effectiveness on model performance is assessed. A Bayesian optimization approach is used to systematically arrive at the hyperparameters of the proposed model. Validation of the presented approach on a survey data collected by the National Rural Electric Cooperative Association in Virginia in 2018 demonstrates the excellent predictive capability of the machine learning based approach to modeling solar adoption reliably.

14 SOLAR ENERGY↗

Effects of various parameters of different porous transport layers in proton exchange membrane water electrolysis

Porous transport layers (PTLs) play an important role in proton exchange membrane water electrolysis (PEMWE) cells. The PTL facilitates water and gas transport, as well as thermal and electrical conduction, and is required to sustain good contact with adjacent components. It is expected that using PTLs with variations in material properties such as structure, composition, thickness and wettability results in performance changes of the PEMWE. Here, a general mathematical PEMWE model is developed that separates and analyzes the contributions of ohmic, activation, diffusion and Nernst potentials. For model validation, three inherently different anode PTL structures (carbon paper, sintered titanium particles, and titanium felt) are operated over a range of conditions. Additionally, the effects of PTL wettability were used to verify the model using Polytetrafluoroethylene (PTFE) treated Toray papers with PTFE loading ranging from 0% to 20%. The modeling results of both PTFE treated and untreated materials show good agreement with the experimental data. Mass transport or diffusion loss is the primary reason for performance differences between PTFE treated and untreated PTLs. Sintered titanium PTLs with thicknesses above 1 mm suffer from up to 33% increased ohmic losses without indicating any obvious changes in activation and diffusion losses when compared to untreated PTLs. The losses of the cell increase when using PTFE treated Toray paper. Individual contributions are quantified and assigned to increased ohmic, activation, and diffusion losses. In conclusion, the proposed model offers insights into the overpotential contributions of a PEMWE. It is a useful tool for predicting performance of various PTL materials and can be applied for PTL development and optimization efforts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Some Basic GN&C Modeling-and-Design Gap Areas for Low-Gee Slosh: Comparisons to Powered Flight Slosh

Potential gaps between powered-flight slosh and near-zero-gee slosh modeling intuition ... - Extremely-basic items like inertia (which is a strong function of frequency and liquid depth) - Visualization of response – ability to form a physical picture via mechanical analog - Availability of frequency-domain design insight from open-loop plant - Sense for edge-of model-validity – how does it behave “at the edges” And some suggestions on how NASA can help

Guidance Navigation and Control Modeling and Desig↗

Data-driven emulation of modal aerosol microphysics via neural operator-based modeling

The complexity and the small characteristic scales of aerosol microphysical processes pose a big challenge for accurate and efficient Earth system simulations at regional and global scales. In this work, we construct and evaluate a surrogate model: the aerosol deep operator network (ADON), a physics-inspired dual-net architecture for emulating the aerosol microphysics parameterization suite in the version 2 of the Energy Earth System Model (E3SMv2). The current version of the surrogate model is trained on a dataset comprising 9.8 million samples obtained from a global E3SMv2 simulation with the horizontal resolution of about one degree under cloud-free conditions. Incorporating domain spatial and temporal coordinates, as well as principle components extracted from training data, the dual-net surrogate model effectively captures the intricate representations of aerosol and the relationship with atmospheric state variables, achieving an R-squared score over $$95.7\%$$ for all the lognormal aerosol modes in the extrapolated regime. The validated model provides feature importance of input variables and their impact on the predictive capacity of the surrogate model in relation to the E3SM. The computational cost of online inference time deployed on CPUs and GPUs with lower precisions highlights ADON’s efficiency and potential in robust predictive modeling for large-scale Earth system computations.

Bai, Zhe↗

Computational aspects of real-time simulation of rotary-wing aircraft

A study was conducted to determine the effects of degrading a rotating blade element rotor mathematical model suitable for real-time simulation of rotorcraft. Three methods of degradation were studied, reduction of number of blades, reduction of number of blade segments, and increasing the integration interval, which has the corresponding effect of increasing blade azimuthal advance angle. The three degradation methods were studied through static trim comparisons, total rotor force and moment comparisons, single blade force and moment comparisons over one complete revolution, and total vehicle dynamic response comparisons. Recommendations are made concerning model degradation which should serve as a guide for future users of this mathematical model, and in general, they are in order of minimum impact on model validity: (1) reduction of number of blade segments; (2) reduction of number of blades; and (3) increase of integration interval and azimuthal advance angle. Extreme limits are specified beyond which a different rotor mathematical model should be used.

Houck, J. A.↗

Multiple Scattering of Laser Pulses in Snow Over Ice: Modeling the Potential Bias in ICESat Altimetry

The primary goal of NASA's current ICESat and future ICESat2 missions is to map the altitude of the Earth's land ice with high accuracy using laser altimetry technology, and to measure sea ice freeboard. Ice however is a highly transparent optical medium with variable scattering and absorption properties. Moreover, it is often covered by a layer of snow with varying depth and optical properties largely dependent on its age. We describe a modeling framework for estimating the potential altimetry bias caused by multiple scattering in the layered medium. We use both a Monte Carlo technique and an analytical diffusion model valid for optically thick media. Our preliminary numerical results are consistent with estimates of the multiple scattering delay from laboratory measurements using snow harvested in Greenland, namely, a few cm. Planned refinements of the models are described.

Davis, A. B.↗

Human Mars Entry, Descent, and Landing Architecture Study: Phase 3 Summary

Over the past four years, NASA has directed the Entry, Descent and Landing Architecture Study (EDLAS) team to evaluate candidate technologies to deliver human-scale vehicles (carrying 20t payloads) to a precise location on the surface of Mars. The study, which initially considered four candidate vehicles, narrowed the design space to focus on two vehicles in Phase 3, one low and one mid lift-to-drag vehicle. Key design challenges exist for both, and the purpose of the Phase 3 analysis was to identify specific technology investment areas and opportunities to mature the vehicle designs beyond simulations to include ground and flight tests. This paper summarizes the detailed analyses performed on the two vehicle configurations, including aerodynamic and propulsive interference effects during the powered flight phase, vehicle packaging, as well as outer mold line and parametric mass model upgrades. The analyses were used to update models in the vehicle performance simulations. The simulation results showing the impact of the Phase 3 analyses on vehicle performance are also presented. Finally, a summary of the technology investment recommendations, including opportunities to validate models using wind tunnel tests and evaluate technologies at the moon, are presented. This paper offers a systems level overview of the more detailed analysis that will be presented in this special session.

Alicia Dwyer Cianciolo↗

Continuum Model to Define the Chemistry and Mass Transfer in a Bicarbonate Electrolyzer

Bicarbonate electrolyzers are devices designed to convert CO 2 captured from point sources or the atmosphere into chemicals and fuels without needing to first isolate pure CO 2 gas. In this work, we report here an experimentally validated model that quantifies the reaction chemistry and mass transfer processes within the catalyst layer and cation exchange membrane layer of a bicarbonate electrolyzer. Our results demonstrate that two distinct chemical microenvironments are key to forming CO at high rates: an acidic membrane layer that promotes in situ CO 2 formation and a basic catalyst layer that suppresses the hydrogen evolution reaction. We show that the rate of CO product formation can be increased by modulating the catalyst and membrane layer properties to increase the rate of in situ CO 2 generation and transport to the cathode. These insights serve to inform the design of bicarbonate and BPM-based CO 2 electrolyzers while demonstrating the value of modeling for resolving rate-determining processes in electrochemical systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

SCALE 6.2.4 Validation: Nuclear Criticality Safety

The computational bias of criticality safety computer codes must be established through the validation of the codes to critical experiments. A large collection of suitable experiments has been vetted by the International Criticality Safety Benchmark Evaluation Project (ICSBEP) and made available in the International Handbook of Evaluated Criticality Safety Benchmark Experiments (ICSBEP Handbook). More than 600 cases from this handbook have been prepared and reviewed within the Verified, Archived Library of Inputs and Data (VALID), which is maintained by the Reactor and Nuclear Systems Division at Oak Ridge National Laboratory. The performance of the KENO V.a and KENO-VI Monte Carlo codes within the SCALE 6.2.4 code system is assessed using the VALID models of benchmark experiments. A range of nuclear cross section libraries based on Evaluated Nuclear Data File (ENDF)/B-VII.1 in both multigroup (MG) and continuous energy (CE) formats is considered. The critical experiments available to validate the KENO V.a code cover 15 broad categories of systems. These systems use a range of fissile materials, including a range of uranium enrichments, various plutonium isotopic vectors, and some mixed uranium/plutonium oxides. The physical forms of the fissile material also vary and are represented as metal, solutions, or arrays of rods or plates in a water moderator. The neutron energy spectra of the systems also vary and cover fast, intermediate, mixed, and thermal spectra. Over 550 of the total cases use the KENO V.a code for the four nuclear data libraries considered in this report.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Thermodynamics-guided machine learning model for predicting convective boundary layer height and its multi-site applicability

Accurate estimation of convective boundary layer height (CBLH) is vital for weather, climate, and air quality modeling. Machine learning (ML) shows promise in CBLH prediction, but input parameter selection often lacks physical grounding, limiting generalizability. This study introduces a novel ML framework for CBLH prediction, integrating thermodynamic constraints and the diurnal CBLH cycle as an implicit physical guide. Boundary layer growth is modeled as driven by surface heat fluxes and atmospheric heat absorption represented with the low tropospheric stability, using the diurnal cycle as input and output. TPOT and AutoKeras are employed to select optimal models, validated against Doppler lidar-derived CBLH data, achieving an R 2 of 0.84 across untrained years. Comparisons of eddy covariance (ECOR) and energy balance Bowen ratio (EBBR) flux measurements show the same prediction capability. Models trained on the ARM SGP C1 site with ECOR data and tested at E37 and E39 yield R 2 values of 0.79 and 0.81, respectively, demonstrating their adaptability. The ML model trained with all sites' data slightly enhances the performance compared with ML models trained over single-site data. The interquartile range for predicted CBLH is consistently narrower than that for DL-derived CBLH, reflecting lower variability in predicted CBLH compared to DL-derived CBLH, which is influenced by additional factors, which are not well represented with the model inputs. The model's generalizability across multiple sites at the ARM SGP site demonstrates its potential for transfer to greater distances, offering a scalable approach for enhancing boundary layer parameterization in atmospheric models.

Chu, Yufei [Stony Brook Univ., NY (United States)]↗

Flux Comparison of Master-8 and Ordem 3.1 Modelled Space Debris Population

With ESA’s Meteoroid And Space debris Terrestrial Environment Reference (MASTER-8) model and NASA’s Orbital Debris Engineering Model(ORDEM) 3.1, the two premier orbital debris engineering models have been officially released. The two models come with significant enhancements and now represent state-of-the-art orbital debris modelling for their respective agencies. Both models provide the community with estimates of the space debris environment from low Earth orbit (LEO) up to at least geostationary altitude. The MASTER population is an event-based simulation of all known events that generate debris and objects that are part of the U.S. Space Surveillance Network (SSN) catalog, which provides coverage of objects with diameters down to approximately 10cm in LEO and 1 m in geosynchronous Earth orbit (GEO). Different models are used to simulate the artificial objects and their orbit evolution over time. These models are called “sources” since they assign an origin to each individual object and consist of fragments, solid rocket motor (SRM) remainders, sodium-potassium(NaK) droplets, paintflakes, ejecta, and multi-layer insulation (MLI) fragments. The objects from each source are characterized by having individual release mechanisms, as well as orbital distributions, material composition, size, and mass distributions. Dedicated radar and telescope observation data is used to calibrate the model for objects larger than 1 cm in LEO and larger than 10 cm in GEO. For calibrating the small-sized objects, below 1 cm, impact data from returned surfaces are analyzed. Because the >1cm object population is dominated by fragments, the fragmentation event database was updated to include new events, as well as re-evaluate past events. Special attention was drawn to re-evaluating theFengyun-1C anti-satellite test from 2007 and Cosmos-Iridium collision event from 2009 since these events shape the fragment population because of their severity. After 2009, the two largest fragmentations in terms of number of tracked debris are the Briz-M explosion in 2012 and the NOAA-16explosion in 2015. In total, there are 261 confirmed fragmentations in the database up to November 2016.The baseline population for ORDEM 3.1 is based on the U.S. SSN catalog, and observational datasets from radar, in situ, and optical sources provide a foundation from which the model populations are statistically extrapolated to smaller size regions. These regions are not well-covered by the SSN catalog yet may pose the greatest threat to operational spacecraft. The NASA Standard Satellite Breakup Model is used to generate fragments greater than 1mm from collisions and explosions, and these fragment populations are scaled using ground-based radar data. Specific major debris-producing events, including the Fengyun-1C, Iridium 33, and Cosmos2251 debris clouds, and unique populations, such as NaK droplets, were re-examined, modelled, and added to the ORDEM environment separately. Optical measurement data is used to model the GEO population down to 10 cm. The debris environment is propagated using NASA’s LEO-to-GEO Environment Debris model, and future explosions of intact objects and collisions involving objects greater than 10 cm are assessed statistically. The environment from a few millimetres down to 10 𝜇m is modelled using a special degradation model where small particles are generated from intact spacecraft and rocket bodies, then the populations are scaled to fit in situ cratering data from Space Shuttle returned surfaces. Fragments smaller than 10 cm are differentiated based on material density categories, i.e., high-, medium-, and low-density, to better characterize the potential debris risk posed to upper stages and spacecraft. This paper will discuss the MASTER and ORDEM approaches for modelling populations and compare fluxes for specific orbits, including sun-synchronous, ISS-altitude, geosynchronous transfer, and GEO. In the end, a conclusion is drawn towards the importance of having multiple fundamentally different, yet validated, models to estimate the space debris population.

Andre Horstmann↗

Revealing the critical role of radical-involved pathways in high temperature cyclopentanone pyrolysis

Cyclopentanone (CPO) is a promising biofuel for spark-ignition engines due to its ring strain and high auto-ignition resistance. Understanding CPO decomposition is crucial for building a high-temperature combustion model. Here we present a comprehensive kinetic model for high-temperature pyrolysis of CPO with verified results from high-pressure shock tube (HPST) measurements. The time- histories of carbon monoxide (CO), ethylene (C 2 H 4 ), and CPO absorbances over the temperature range of 1156-1416 K and pressure range of 8.53-10.06 atm were measured during current experiments. A corresponding detailed kinetic model was generated using the Reaction Mechanism Generator (RMG) with dominant unimolecular/radical-involved decomposition pathways from either previous studies or quantum calculations within the current work. The obtained model containing 821 species and 79,859 reactions exhibited a good agreement with the experimental results. In this study, the absorbance ratio between C 2 H 4 and CO was used as an important factor to validate models and to prove that radical-involved bimolecular pathways were as significant as unimolecular decomposition of CPO. The rate of production (ROP) analysis showed H radicals play a major role in the decomposition, and the whole decomposition process could be divided into three stages based on the H radical concentration. Finally, the insights from present work can be used to generate a better CPO combustion model and help evaluate CPO as an advanced biofuel.

09 BIOMASS FUELS↗