SEARCH · Engineering Papers

Results for “Mixture modeling”

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,513 records · Page 84

All sources NASA NTRS DOE OSTI

NASA NTRS1990

Thermal/structural analyses of several hydrogen-cooled leading-edge concepts for hypersonic flight vehicles

The aerodynamic heating at high flight Mach numbers, when shock interference heating is included, can be extremely high and can exceed the capability of most conventional metallic and potential ceramic materials available. Numerical analyses of the heat transfer and thermal stresses are performed on three actively cooled leading-edge geometries (models) made of three different materials to address the issue of survivability in a hostile environment. These analyses show a mixture of results from one configuration to the next. Results for each configuration are presented and discussed. Combinations of enhanced internal film coefficients and high material thermal conductivity of copper and tungsten are predicted to maintain the maximum wall temperature for each concept within acceptable operating limits. The exception is the TD nickel material which is predicted to melt for most cases. The wide range of internal impingement film coefficients (based on correlations) for these conditions can lead to a significant uncertainty in expected leading-edge wall temperatures. The equivalent plastic strain, inherent in each configuration which results from the high thermal gradients, indicates a need for further cyclic analysis to determine component life.

Gladden, Herbert J.↗

NASA NTRS1991

Calculations and curve fits of thermodynamic and transport properties for equilibrium air to 30000 K

A self-consistent set of equilibrium air values were computed for enthalpy, total specific heat at constant pressure, compressibility factor, viscosity, total thermal conductivity, and total Prandtl number from 500 to 30,000 K over a range of 10(exp -4) atm to 10(exp 2) atm. The mixture values are calculated from the transport and thermodynamic properties of the individual species provided in a recent study by the authors. The concentrations of the individual species, required in the mixture relations, are obtained from a free energy minimization calculation procedure. Present calculations are based on an 11-species air model. For pressures less than 10(exp -2) atm and temperatures of about 15,000 K and greater, the concentrations of N(++) and O(++) become important, and consequently, they are included in the calculations determining the various properties. The computed properties are curve fitted as a function of temperature at a constant value of pressure. These curve fits reproduce the computed values within 5 percent for the entire temperature range considered here at specific pressures and provide an efficient means for computing the flowfield properties of equilibrium air, provided the elemental composition remains constant at 0.24 for oxygen and 0.76 for nitrogen by mass.

Gupta, Roop N.↗

NASA NTRS1993

Models of thermal relaxation mechanics for particle simulation methods

A comparative study is conducted of the simulation of molecular vibrational relaxation in rarefied flows by means of existing models and a novel phenomenological microscopic model. This improved technique iterates between translation-rotation and rotation-vibration exchanges; equilibrium is in this way promoted if the model for the latter process is compatible with quantized oscillators. The iteration-equipartition model retains computational simplicity and promotes thermal equilibrium, even in the case of nondegenerate anharmonic quantized oscillators and multispecies gas mixtures.

Haas, Brian L.↗

NASA NTRS1997

Force Balance and Substorm Effects in the Magnetotail

A model of the quiet time middle magnetotail is developed using a consistent orbit tracing technique. The momentum equation is used to calculate geocentric solar magnetospheric components of the particle and electromagnetic forces throughout the current sheet. Ions generate the dominant x and z force components. Electron and ion forces almost cancel in the y direction because the two species drift earthward at comparable speeds. The force viewpoint is applied to a study of some substorm processes. Generation of the rapid flows seen during substorm injection and bursty bulk flow events implies substantial force imbalances. The formation of a substorm diversion loop is one cause of changes in the magnetic field and therefore in the electromagnetic force. It is found that larger forces are produced when the cross-tail current is diverted to the ionosphere than would be produced if the entire tail current system simply decreased. Plasma is accelerated while the forces are unbalanced resulting in field lines within a diversion loop becoming more dipolar. Field lines become more stretched and the plasma sheet becomes thinner outside a diversion loop. Mechanisms that require thin current sheets to produce current disruption then can create additional diversion loops in the newly thinned regions. This process may be important during multiple expansion substorms and in differentiating pseudoexpansions from full substorms. It is found that the tail field model used here can be generated by a variety of particle distribution functions. However, for a given energy distribution the mixture of particle mirror or reflection points is constrained by the consistency requirement. The study of uniqueness also leads to the development of a technique to select guiding center electrons that will produce charge neutrality all along a flux tube containing nonguiding center ions without the imposition of a parallel electric field.

Kaufmann, Richard L.↗

NASA NTRS1995

Stealth Plumes on Io

We suggest that Io's eruptive activity may include a class of previously undetected SO2 geysers. The thermodynamic models for the eruptive plumes discovered by Voyager 'involve low to moderate entropy SO2 eruptions. The resulting plumes are a mixture of solid and gas which emerge from the vent and follow essentially ballistic trajectories. We show that intrusion of silicate magma into buried SO2 deposits can create the required conditions for high entropy eruptions which proceed entirely in the vapor phase. These purely gaseous plumes would have been invisible to Voyager's instruments. Hence, we call them "stealth" plumes. Such eruptions could explain the "patchy" SO2 atmosphere inferred from recent UV and micro-wave spectral observations. The magma intrusion rate required to support the required gas production for these plumes is a negligible fraction of estimated global magma intrusion rates.

Johnson, T. V.↗

NASA NTRS2000

Characterizing the Use of Ultrasonic Energy in Promoting Uniform Microstructural Dispersions in Immiscible Mixtures

Gravity driven separation and preferential wetting precludes uniform microstructural distributions during solidification processing of immiscible, liquid-liquid mixtures. Historically, it is, however, established that liquid/liquid suspensions can be established and maintained by utilizing ultrasound. Following a brief introduction the results of experiments on immiscible mixtures subjected to ultrasonic energy during solidification processing will be compared and evaluated in view of a recently developed mathematical model. The presentation continues by discussion of scaling the model to commercial viability and concludes with the implications of such processing in a microgravity environment.

Grugel, R. N.↗

NASA NTRS2001

Characterizing the Use of Ultrasonic Energy in Promoting Uniform Microstructural Dispersions in Immiscible Mixtures

Grugel, R. N.↗

NASA NTRS2004

Numerical Simulation of Nanostructure Growth

Nanoscale structures, such as nanowires and carbon nanotubes (CNTs), are often grown in gaseous or plasma environments. Successful growth of these structures is defined by achieving a specified crystallinity or chirality, size or diameter, alignment, etc., which in turn depend on gas mixture ratios. pressure, flow rate, substrate temperature, and other operating conditions. To date, there has not been a rigorous growth model that addresses the specific concerns of crystalline nanowire growth, while demonstrating the correct trends of the processing conditions on growth rates. Most crystal growth models are based on the Burton, Cabrera, and Frank (BCF) method, where adatoms are incorporated into a growing crystal at surface steps or spirals. When the supersaturation of the vapor is high, islands nucleate to form steps, and these steps subsequently spread (grow). The overall bulk growth rate is determined by solving for the evolving motion of the steps. Our approach is to use a phase field model to simulate the growth of finite sized nanowire crystals, linking the free energy equation with the diffusion equation of the adatoms. The phase field method solves for an order parameter that defines the evolving steps in a concentration field. This eliminates the need for explicit front tracking/location, or complicated shadowing routines, both of which can be computationally expensive, particularly in higher dimensions. We will present results demonstrating the effect of process conditions, such as substrate temperature, vapor supersaturation, etc. on the evolving morphologies and overall growth rates of the nanostructures.

Hwang, Helen H.↗

NASA NTRS2009

Insights into the Galactic Cosmic-ray Source from the TIGER Experiment

We report results from 50 days of data accumulated in two Antarctic flights of the Trans-Iron Galactic Element Recorder (TIGER). With a detector system composed of scintillators, Cherenkov detectors, and scintillating optical fibers, TIGER has a geometrical acceptance of 1.7 sq m sr and a charge resolution of 0.23 cu at Iron. TIGER has obtained abundance measurements of some of the rare galactic cosmic rays heavier than iron, including Zn, Ga, Ge, Se, and Sr, as well as the more abundant lighter elements (down to Si). The heavy elements have long been recognized as important probes of the nature of the galactic cosmic-ray source and accelerator. After accounting for fragmentation of cosmic-ray nuclei as they propagate through the Galaxy and the atmosphere above the detector system, the TIGER source abundances are consistent with a source that is a mixture of about 20% ejecta from massive stars and 80% interstellar medium with solar system composition. This result supports a model of cosmic-ray origin in OB associations previously inferred from ACE-CRIS data of more abundant lighter elements. These TIGER data also support a cosmic-ray acceleration model in which elements present in interstellar grains are accelerated preferentially compared with those found in interstellar gas.

Link, Jason T.↗

NASA NTRS2019

Constraining the Metallicities, Ages, Star Formation Histories, and Ionizing Continua of Extragalactic Massive Star Populations

We infer the properties of massive star populations using the far-ultraviolet stellar continua of 61 star-forming galaxies: 42 at low redshift observed with the Hubble Space Telescope and 19 at z ~ 2 from the MegaSaura sample. We fit each stellar continuum with a linear combination of up to 50 single-age and single-metallicity starburst99 models. From these fits, we derive light-weighted ages and metallicities, which agree with stellar wind and photospheric spectral features, and infer the spectral shapes and strengths of the ionizing continua. Inferred light-weighted stellar metallicities span 0.05–1.5 Z(sub ⊙) and are similar to the measured nebular metallicities. We quantify the ionizing continua using the ratio of the ionizing flux at 900 Å to the non-ionizing flux at 1500 Å and demonstrate the evolution of this ratio with stellar age and metallicity using theoretical single-burst models. These single-burst models only match the inferred ionizing continua of half of the sample, while the other half are described by a mixture of stellar ages. Mixed-age populations produce stronger and harder ionizing spectra than continuous star formation histories, but, contrary to previous studies that assume constant star formation, have similar stellar and nebular metallicities. Stellar population age and metallicity affect the far-UV continua in different and distinguishable ways; assuming a constant star formation history diminishes the diagnostic power. Finally, we provide simple prescriptions to determine the ionizing photon production efficiency (ξ(sub ion)) from the stellar population properties. The ξ(sub ion) inferred from the observed star-forming galaxies has a range of log(ξ(sub ion)) = 24.4–25.7 Hz erg(exp −1) that depends on the stellar population age, metallicity, star formation history, and contributions from binary star evolution. These stellar population properties must be observationally determined to accurately determine the number of ionizing photons generated by massive stars.

Chisholm, J.↗

NASA NTRS2019

Modeling and simulation challenges for softgoods in space deployable structures

This paper examines modeling and simulation challenges for deployable spacecraft structures with softgood components. Previous studies highlighted the numerical challenges due to the mixture of high stiffness and low stiffness contact in this class of structures. This paper examines the use of the commercial LS-DYNA nonlinear finite element code on a set of relevant benchmark problems.

Mobrem, Mehran↗

DOE OSTI2021

Experimental and modeling studies of PR and ND oxalate solubility to high ionic strengths: Insight into actinide(III) oxalates

Actinide oxalates are chemical compounds important to nuclear industry, ranging from actinide separation in waste reprocessing, to production of specialty actinides, and to disposal of high level nuclear waste (HLW) and spent nuclear fuel (SNF). In this study, the solubility constants for Pr 2 (C 2 O 4 ) 3 •10H 2 O and Nd 2 (C 2 O 4 ) 3 •10H 2 O by performing solubility experiments in HNO 3 and mixtures of HNO 3 and H 2 C 2 O 4 at 23.0 ± 0.2 °C have been determined. The targeted starting materials, Pr 2 (C 2 O 4 ) 3 •10H 2 O and Nd 2 (C 2 O 4 ) 3 •10H 2 O, were successfully synthesized at room temperature using PrCl 3 , NdCl 3 and oxalic acid as the source metrials. Then, we utilized the targeted solubility-controlling phases to conduct solubility measurements. There was no phase change over the entire periods of experiments, demonstrating that Pr 2 (C 2 O 4 ) 3 •10H 2 O and Nd 2 (C 2 O 4 ) 3 •10H 2 O were the solubility-controlling phases in our respective experiments. Based on our experimental data, we have developed a thermodynamic model for Pr 2 (C 2 O 4 ) 3 •10H 2 O and Nd 2 (C 2 O 4 ) 3 •10H 2 O in the mixtures of HNO 3 and H 2 C 2 O 4 to high ionic strengths. The model for Pr 2 (C 2 O 4 ) 3 •10H 2 O reproduces well the reported experimental data for Pu 2 (C 2 O 4 ) 3 •10H 2 O, which are not utilized for the model development, demonstrating that Pr(III) is an excellent analog for Pu(III). Similarly, the model for Nd 2 (C 2 O 4 ) 3 •10H 2 O reproduces the solubility of Am 2 (C 2 O 4 ) 3 •10H 2 O and Cm 2 (C 2 O 4 ) 3 •10H 2 O. The Pitzer model was used for the calculation of activity coefficients. Based on the published, well established model for dissociation constants for oxalic acid and stability constants for actinide-oxalate complexes [i.e., AmC 2 O 4 + , and Am(C 2 O 4 ) 2 - ] to high ionic strengths, we have obtained the solubility constants (log 10 K 0 ) for the following reactions at 25 °C, Pr 2 (C 2 O 4 ) 3 •10H 2 O ⇌ 2Pr 3+ + 3C 2 O 4 2- + 10H 2 O(l). Nd 2 (C 2 O 4 ) 3 •10H 2 O ⇌ 2Nd 3+ + 3C 2 O 4 2- + 10H 2 O(l). to be -30.82 ± 0.30 (2σ), and - 31.14 ± 0.35 (2σ), respectively. These values for can be directly applied to Pu 2 (C 2 O 4 ) 3 •10H 2 O, Am 2 (C 2 O 4 ) 3 •10H 2 O and Cm 2 (C 2 O 4 ) 3 •10H 2 O. The model established for actinide oxalates by this study provides the needed knowledge with regard to solubilities of actinide/REE oxalates at various ionic strengths, and is expected to find applications in many fields, including the geological disposal of nuclear waste and the mobility of REE under the surface conditions, as Pr 2 (C 2 O 4 ) 3 •10H 2 O and Nd 2 (C 2 O 4 ) 3 •10H 2 O can be regarded as the pure Pr and Nd end-members of deveroite, a recently discovered natural REE oxalate with the following stoichiometry, (Ce 1.01 Nd 0.33 La 0.32 Pr 0.11 Y 0.11 Sm 0.01 Pb 0.04 U 0.03 Th 0.01 Ca 0.04 ) 2.01 (C 2 O 4 ) 2.99 •9.99H 2 O. Regarding its importance in the geological disposal of nuclear waste, Am 2 (C 2 O 4 ) 3 •10H 2 O/Pu 2 (C 2 O 4 ) 3 •10H 2 O/Cm 2 (C 2 O 4 ) 3 •10H 2 O can be the source-term phase for actinides, as demonstrated by the instance in the disposal in clay/shale formations. This is exemplified by the stability of Am 2 (C 2 O 4 ) 3 •10H 2 O in comparison with Am(OH) 3 (am), Am(OH) 3 (s) and AmCO 3 (OH)(s) under the relevant geological repository conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Assessment of turbulent combustion models for simulating pre-chamber ignition in a natural gas engine

Lean combustion in an internal combustion engine is a promising strategy to increase thermal efficiency by leveraging a more favorable specific heat ratio of the fresh mixture and simultaneously suppressing the heat losses to the cylinder wall. However, unstable ignition events and slow flame propagation at fuel-lean condition lead to high cycle-to-cycle variability and hence limit the high-efficiency engine operating range. Pre-chamber ignition is considered an effective concept to extend the lean operating limit, by providing spatially distributed ignition with multiple turbulent flame-jets and enabling faster combustion rate compared to the conventional spark ignition approach. From a numerical modeling perspective, to date, still the science base and available simulation tools are inadequate for understanding and predicting the combustion processes in pre-chamber ignited engines. In this paper, conceptually different RANS combustion models widely adopted in the engine modeling community were used to simulate the ignition and combustion processes in a medium-duty natural gas engine with a pre-chamber spark-ignition system. A flamelet-based turbulent combustion model, i.e., G-equation, and a multi-zone well-stirred reactor model were employed for the multi-dimensional study. Simulation results were compared with experimental data in terms of in-cylinder pressure and heat release rate. Finally, the analysis of the performance of the two models is carried out to highlight the strengths and limitations of the two formulations respectively.

combustion model evaluation↗

DOE OSTI2022

Short- to Intermediate-Range Structure, Transport, and Thermophysical Properties of LiF–NaF–ZrF4 Molten Salts

LiF–NaF–ZrF 4 multicomponent molten salts are identified as promising candidates for coolant salts in molten salt reactors and advanced high-temperature reactors. This study focused on low-melting point salt compositions of interest: 38LiF–51NaF–11ZrF 4 , 42LiF–29NaF–29ZrF 4 , and 26LiF–37NaF–37ZrF 4 . Ab-initio molecular dynamics (AIMD) calculations were performed and compared with available experimental data to assess the ability of rigid ion models (RIM) to reproduce short- to intermediate-range structure, transport, and thermophysical properties of the LiF–NaF–ZrF 4 salt mixtures. It is found that as ZrF 4 mol% increases, the average cation–anion coordination number (CN) of monovalent cations (Li + , Na + ) obtained from RIM calculations decreases, while multivalent Zr 4+ CN varied from 15% to 19% in comparison to corresponding AIMD values. In addition, RIM is found to predict the existence of 7, 8, and 9 coordinated fluorozirconate complexes, while AIMD and the available experimental data showed an occurrence of 6, 7, and 8 coordinated complexes in the melt. The intermediate-range structure analysis revealed that while the RIM parameters are able to reproduce a local structure for lower ZrF 4 mol% salts such as in 38LiF–51NaF–11ZrF 4 , an extensive fluorozirconate network formation is observed in RIM simulations for higher ZrF 4 mol% compositions. The network generated by RIM parameters is found to be mainly connected by “corner-sharing” fluorozirconate complexes as opposed to both “edge-sharing” and “corner-sharing” connectively portrayed by AIMD. It is found that a close agreement between AIMD and the RIM salt structure for the 11-mol% ZrF 4 salt resulted in good agreement in the calculated Zr diffusivities and the viscosity values. However, due to the inaccurate short- to intermediate-range structure prediction by RIM for higher ZrF 4 mol% compositions, thermophysical properties such as densities and heat capacity differ by up to 26% and 27%, respectively, upon comparison with AIMD and experimental values. Also, the network-dominated properties such as diffusion coefficients and viscosities differed by up to two and three orders of magnitude, respectively. This study signifies the importance of accurate salt structure generation for an accurate prediction of transport and thermophysical properties of multicomponent molten salts.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

An experimental and chemical kinetic modeling study of 4-butoxyheptane combustion

Here, the combustion kinetics of a novel oxygenated bioblendstock for diesel, 4-butoxyheptane (4-BH), was investigated experimentally using a flow reactor and a heated, high-pressure shock tube. The flow reactor experiments employed oxygen as the oxidizer and helium as the diluent with oxidation conducted at atmospheric pressure and 10 bar for temperatures from 400 to 1000 K at 20-K intervals. The fuel, oxidizer, and diluent flow rates were varied at different temperatures to maintain a constant initial fuel mole fraction of 1000 ppm, with stoichiometric equivalence ratio, and a residence time of 2.0 s. The reacted gas was fed to two separate GC systems that could qualitatively and quantitatively detect product species. Additionally, real fuel-air ignition delay time (IDT) data were collected using a heated, high-pressure shock-tube facility. Fuel lean (φ = 0.5) and stoichiometric (φ = 1.0) mixtures were investigated at 10 atm as well as at 30 atm for the fuel lean case for temperatures between 847 and 1259 K. A detailed chemical kinetics mechanism was developed to model the product distribution from the flow reactor and IDTs from the shock tube. The proposed model was able to predict the double NTC behavior in flow reactor experiments reasonably well. Model predictions at low temperatures were observed to be highly sensitive to the rate constants of ketohydroperoxide (KHP) decomposition in the case of the OOH group in α position which were modeled based on existing literature studies on ethers. It was noted that in the absence of theoretical or experimental studies, the rate constants for KHP decomposition used in the literature were empirically set. Additional studies are required to address the gap in model prediction obtained in this study and to reduce the uncertainty in kinetics models for ether oxidation. Predicted product concentrations and IDTs showed some quantitative agreement with experimental data, but the overall reactivity of the IDTs is underpredicted. Additionally, significant deviation is observed for the IDT results at 10 atm for the stoichiometric case with minor deviations for the other cases. The reaction pathways to the missing products were then further analyzed theoretically through quantum-mechanical calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

iNovo479: Metabolic Modeling Provides a Roadmap to Optimize Bioproduct Yield from Deconstructed Lignin Aromatics by Novosphingobium aromaticivorans

Lignin is an abundant renewable source of aromatics and precursors for the production of other organic chemicals. However, lignin is a heterogeneous polymer, so the mixture of aromatics released during its depolymerization can make its conversion to chemicals challenging. Microbes are a potential solution to this challenge, as some can catabolize multiple aromatic substrates into one product. Novosphingobium aromaticivorans has this ability, and its use as a bacterial chassis for lignin valorization could be improved by the ability to predict product yields based on thermodynamic and metabolic inputs. In this work, we built a genome-scale metabolic model of N. aromaticivorans, iNovo479, to guide the engineering of strains for aromatic conversion into products. iNovo479 predicted product yields from single or multiple aromatics, and the impact of combinations of aromatic and non-aromatic substrates on product yields. We show that enzyme reactions from other organisms can be added to iNovo479 to predict the feasibility and profitability of producing additional products by engineered strains. Thus, we conclude that iNovo479 can help guide the design of bacteria to convert lignin aromatics into valuable chemicals.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

A multiplexed nanostructure-initiator mass spectrometry (NIMS) assay for simultaneously detecting glycosyl hydrolase and lignin modifying enzyme activities

Lignocellulosic biomass is composed of three major biopolymers: cellulose, hemicellulose and lignin. Analytical tools capable of quickly detecting both glycan and lignin deconstruction are needed to support the development and characterization of efficient enzymes/enzyme cocktails. Previously we have described nanostructure-initiator mass spectrometry-based assays for the analysis of glycosyl hydrolase and most recently an assay for lignin modifying enzymes. Here we integrate these two assays into a single multiplexed assay against both classes of enzymes and use it to characterize crude commercial enzyme mixtures. Application of our multiplexed platform based on nanostructure-initiator mass spectrometry enabled us to characterize crude mixtures of laccase enzymes from fungi Agaricus bisporus (Ab) and Myceliopthora thermophila (Mt) revealing activity on both carbohydrate and aromatic substrates. Using time-series analysis we determined that crude laccase from Ab has the higher GH activity and that laccase from Mt has the higher activity against our lignin model compound. Inhibitor studies showed a significant reduction in Mt GH activity under low oxygen conditions and increased activities in the presence of vanillin (common GH inhibitor). Ultimately, this assay can help to discover mixtures of enzymes that could be incorporated into biomass pretreatments to deconstruct diverse components of lignocellulosic biomass.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Analysis of Waste Material Feedstocks Using Laser-Induced Breakdown Spectroscopy and Machine Learning

Predicting properties such as heating value, ash fusion temperature, and mineral ash composition from Laser-Induced Breakdown Spectroscopy (LIBS) data can make gasifiers more flexible to different feedstocks. Understanding these feedstock properties in-situ improves feedstock conversion modelling methods that allow for consistent operation, higher carbon conversion, and reduced fouling and erosion rates. The purpose of this study is to demonstrate methods for model creation that take LIBS data as predictor features and estimate higher order material properties as a function of feedstock material properties. Six samples were chosen to represent a mixture of abundant and carbon rich waste materials. LIBS measurements were performed on these samples for elemental wavelengths and intensity values. Laboratory analytical results were obtained for each sample’s heating value, proximate and ultimate analysis, mineral ash composition, ash fusion temperatures, and viscosity temperatures. Thermal conductivity was measured using a HotDisk TPS 2500S. LIBS measurements were processed and used as predictor features for machine learning (ML) models to predict the sample’s material properties. Predictor feature selection algorithms, particularly minimum redundancy maximum relevance (mRMR), reduced the dimensionality of ML models. Many modelling methods such as Gaussian process regression (GPR), regression tree, neural networks (NN), and support vector machines (SVM) were demonstrated to be effective at predicting higher order properties; however, mRMR with GPR stood out as a clear winning combination.

01 COAL, LIGNITE, AND PEAT↗