Search NASA⌕ Search

SEARCH · Search NASA

Results for “simulation 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.

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

At least 1,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Control System of Multi-Port Autonomous Reconfigurable Solar Power Plant (MARS) & HIL Platforms for Design

Multi-port autonomous reconfigurable solar power plant (MARS) provides an attractive alternative to connect photovoltaic (PV) and energy storage systems (ESSs) to high-voltage direct current (HVdc) links and high-voltage alternating current (ac) grids. In this paper, a unique hierarchical control system of MARS is proposed and evaluated. To evaluate the control system and associated algorithms in early-stage research of complex architectures like MARS, it is important to develop unique suitable hardware-in-the-loop (HIL) platforms. In this paper, the HIL platforms for MARS to evaluate the performance of the hierarchical control system and the control algorithms implemented are presented. Further, they help with the design process of control systems. The real-time simulation models and algorithms that are utilized for MARS in the HIL platforms are also discussed in the paper. The HIL experiments of the control system of MARS showcase the capability to provide continuity of operation under faults and frequency support to the power grid during loss of generation. They also showcase the stability of the proposed hierarchical control system of MARS.

14 SOLAR ENERGY↗
DOE OSTI2023

Response of Extreme Rainfall to Atmospheric Warming and Wetting: Implications for Hydrologic Designs Under a Changing Climate

Understanding the processes of rainfall extremes and their response to anthropogenic climate change is pivotal for improved adaptation of unprecedented flood hazards around the world. Here we take the record-breaking 20 July 2021 storm over central China as an example. We investigate the response of this particular storm to atmospheric warming (i.e., increase in air temperature) and wetting (i.e., increase in atmospheric moisture content) based on a series of convection-permitting model simulations. Our results show non-monotonic changes of the space-time rainfall variability to either increased temperature or atmospheric moisture content. The most extreme rain rate is produced when relative humidity is increased by 20~40 % or temperature is increased by less than 2°C. The non-monotonic rainfall response is more clearly revealed at fine spatial (100-1000 km 2 ) and temporal scales (less than 6 hours) rather than over the entire domain (~10 4 km 2 ) and aggregated over the storm duration (around two days). This is mainly attributable to the distinct feedbacks from atmospheric dynamics (i.e., moisture convergence and interaction with regional topography) rather than regulated by thermodynamic changes alone. Furthermore, atmospheric warming poses notable changes in the vertical structure of storm cells, contributing to reduced areal reduction factors at small spatial scales and short durations, while atmospheric wetting additionally modifies storm evolution properties. Our modeling analyses challenge the existing practices for hydrologic designs under a changing climate, highlighting particular vulnerability for cities or small basins to short-duration rainfall extremes and the resultant flash flood hazards.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

The LAKE model input dataset for three Arctic lakes

This dataset contains meteorological data collected for three Arctic lakes and compiled to satisfy input requirements of the LAKE 2.0 model. The dataset was generated to act as a benchmarking dataset for future model-data inter-comparisons. The LAKE 2.0 model simulates temperatures within the water later and the sedimentary layer of a lake. The LAKE2.0. is an open-source code and available to download via this weblike http://tesla.parallel.ru/Viktor/LAKE/-/wikis/LAKE-model (last visit July 14, 2021). The meteorological data are required to simulate the surface energy balance at the surface of a lake. This dataset includes a compilation of the meteorological data pulled from multiple data streams, including National Oceanic and Atmospheric Administration (NOAA) climate data, Circumarctic Lakes Observation Network (CALON) data, and the United States Geological Survey (USGS) data. The data were compiled for three Arctic lakes: FoxDen (66.55877, -164.45670), Atqasuk (70.452497, -156.951984), and Toolik (68.63150, -149.60740). Each meteorological data is in comma-delimited format (file extension ‘.dat’) and includes eight columns: Temperature [K], Pressure [Pa], longwave downward radiation [W/m2], shortwave downward radiation [W/m2], “U” wind speed [m/s], ”V” wind speed [m/s], humidity [kg/kg], precipitation [m/s]. In addition to the meteorological data file, we included setup and driver files. The Toolik lake is the deepest out of three lakes and has inflowing and outflowing groundwater data. InflowOutflowREADME.txt has more information about inflow and outflow flies. The other two lakes are much shallower and modeled as a closed system (i.e. no water inflow or outflow).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Development of Thermal Scattering Law and Cross Sections of Hydrogen in Paraffin Wax

Paraffin wax is frequently used as a neutron moderator and shielding material. The main component of paraffin wax is straight-chain alkanes (n-alkanes). The deposition of paraffin wax is primarily attributed to the crystallization of n-alkanes. It is important to gain a deeper understanding of the mechanisms underlying the behavior of paraffin wax, which would impact the thermal scattering Law (TSL) and cross sections and affect the analysis of neutronic and critical systems. In this work, a classical molecular dynamics (CMD) simulation model was used in LAMMPS to create the TSL and cross sections at room temperature and pressure. To generate the required velocity autocorrelation functions (VACF), previously published data were used to validate the approach and models by comparing them with key model parameters. The phonon density of state (DOS) was calculated using Fourier transformation of the normalized VACF. This DOS was used as the primary input to estimate the TSL (S($a$, $β$)) and cross sections of hydrogen in paraffin wax. The TSL and cross sections of hydrogen were estimated using the Full Law Analysis Scattering System Hub (FLASSH) code. The cross section of hydrogen in paraffin wax is consistent with other hydrocarbon materials such as polyethylene with deviations due to structure in the lowest energy region.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Anatomy of a Summertime Convective Event over the Arabian Region

This study investigates the structure and evolution of a summertime convective event that occurred on 14 July 2015 over the Arabian region. We use the WRF Model with 1-km horizontal grid spacing and test three PBL parameterizations: the Mellor–Yamada–Nakanishi–Niino (MYNN) scheme; the Asymmetrical Convective Model, version 2, (ACM2) scheme; and the quasi-normal scale-elimination (QNSE) scheme. Convection initiates near the Al Hajar Mountains of northern Oman at around 1100 local time (LT; 0700 UTC) and propagates northwestward. A nonorographic convective band along the west coast of the United Arab Emirates (UAE) develops after 1500 LT as a result of the convergence of cold pools with the sea breeze from the Arabian Gulf. The model simulation employing the QNSE scheme simulates the convection initiation and propagation well. Although the MYNN and ACM2 simulations show convective initiation near the Al Hajar Mountains, they fail to simulate the development of the convective band along the UAE west coast. The MYNN run simulates colder near-surface temperatures and a weaker sea breeze, whereas the ACM2 run simulates a stronger sea breeze but a drier lower troposphere. Sensitivity simulations using horizontal grid spacings of 9 and 3 km show that lower-resolution runs develop broader convective structures and weaker cold pools and horizontal wind divergence, affecting the development of convection along the west coast of the UAE. The 1-km run using the QNSE PBL scheme realistically captures the sequence of events that leads to the moist convection over the UAE and adjacent mountains.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Impact of Open Communication Networks on Load Frequency Control with Plug-In Electric Vehicles by Cyber-Physical Dynamic Co-Simulation

With the increasing electrification of the transportation sector to achieve the carbon neutrality objective, despite the challenges of charging electric vehicles (EV), there are also opportunities through smart charging EVs to improve system frequency stability; however, EV control technologies might require nontraditional communication support. This paper investigates the impacts of communication variations of EV on power system load frequency control through a cyber-physical dynamic system (CPDS) co-simulation. Here, the CPDS is built upon our previously developed transmission-and-distribution dynamic co-simulation model with the added communication variation functions (i.e., delay and packet loss). The case studies consider multiple communication variation scenarios when the system experiences an N-1 generation trip contingency. The scenarios include communication delays and packet loss using both homogeneous and heterogeneous assumptions. The outcomes of this work can help improve EV frequency regulation services and provide robust and effective tests for different load frequency control algorithms of the future power systems.

ADVANCED PROPULSION SYSTEMS,ENERGY PLANNING, POLIC↗
DOE OSTI2021

Agent-based Disease Surveillance and Transmission Model

The SARS-CoV-2 pandemic has driven a need for models that can predict transmission among populations and test the effectiveness of proposed interventions to control it such as testing, quarantining, vaccinations, and nonpharmaceutical interventions. Some questions of interest include:1. How rapidly will the virus spread if no action is taken? 2. How frequently should people be tested and quarantined if necessary? 3. Should unvaccinated persons be tested more frequently than the vaccinated? 4. If test results are delayed, how will that affect the spread? 5. How do test false negative rates affect the overall efficacy of testing programs? 6. How do constraints on testing resources impact the overall efficacy of testing programs? 7. How does the patient’s failure to recognize COVID symptoms and quarantine affect the spread? 8. What level of vaccination coverage is needed to suppress the epidemic? 9. What is the risk of a large outbreak (probability exceeding x infections)? 10. What is the cost effectiveness of interventions (dollars per avoided infection)? This report describes the structure of a discrete event simulation model that provides actionable information to address these questions. Parameter inputs and results of analyses have been reported in briefings provided to LLNL management and are not discussed here.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Mechanistic study of chemical looping reactions between solid carbon fuels and CuO

Copper (Cu) based chemical-looping combustion (CLC) is a promising process that utilizes solid carbon fuel such as coal and biomass. Understanding the reaction kinetics in this process can facilitate the industrial design of CLC units. Here, in this work, molecular dynamics (MD) simulations were performed to investigate the reaction kinetics of n-butane and two simplified solid carbon fuels (lignite and anthracite) with/without copper oxide (CuO) nanoparticle. In addition, experiments were conducted on the thermal characteristics, flammability, and flame speeds of CuO and solid carbon mixtures. For n-butane, oxidation on the CuO surface is the primary reaction since the activation energy of the surface reaction is much lower than that of oxygen (O 2 ) combustion (9.2 vs. 53.3 kcal/mol). On the other hand, for lignite, there is a smaller difference in the activation energy of O 2 combustion and surface reaction (23.04 vs. 6.34 kcal/mol). We hypothesize that some solid carbon fuels can have different reaction kinetics dependent upon temperature. This is proven by the case of a simplified anthracite coal: an Arrhenius plot shows that this solid carbon fuel has two different reaction kinetic regimes and the critical temperature for the change in kinetics is related to the oxygen uncoupling of the CuO nanoparticle. Like the modeling simulations, a change in activation energy is observed in the experimental results, where desorption of molecular oxygen from CuO becomes important.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Inter-well connectivity detection in CO 2 WAG projects using statistical recurrent unit models

Routine well-wise injection and production measurements contain significant information on subsurface structure and properties. Data-driven technology that interprets surface data into subsurface structure or properties can assist operators in making informed decisions by providing a better understanding of field assets. Our machine-learning framework is built on the statistical recurrent unit (SRU) model and interprets well-based injection/production data into inter-well connectivity without relying on a geologic model. We test it on synthetic and field-scale CO 2 EOR projects utilizing the water-alternating-gas (WAG) process. SRU is a special type of recurrent neural network (RNN) that allows for better characterization of temporal trends, by learning various statistics of the input at different time scales. In our application, the complete states (injection rate, pressure and cumulative injection) at injectors and pressure states at producers are fed to SRU as the input and the phase rates at producers are treated as the output. Once the SRU is trained and validated, it is then used to assess the connectivity of each injector to any producer using permutation variable importance method, wherein inputs corresponding to an injector are shuffled and the increase in prediction error at a given producer is recorded as the importance (connectivity metric) of the injector to the producer. This method is tested in both synthetic and field-scale cases. The validation of the proposed data-driven inter-well connectivity assessment is performed using synthetic data from simulation models where inter-well connectivity can be easily measured using the streamline-based flux allocation. The SRU model is shown to offer excellent prediction performance on the synthetic case. Despite significant measurement noise and frequent well shut-ins imposed in the field-scale case, the SRU model offers good prediction accuracy, the overall relative error of the phase production rates at most producers ranges from 10% to 30%. It is shown that the dominant connections identified by the data-driven method and streamline method are in close agreement. This significantly improves confidence in our data-driven procedure. The novelty of this work is that it is purely data-driven method and can directly interpret routine surface measurements to intuitive subsurface knowledge. Furthermore, the streamline-based validation procedure provides physics-based backing to the results obtained from data analytics. This study results in a reliable and efficient data analytics framework that is well-suited for large field applications.

42 ENGINEERING↗
DOE OSTI2026

Urban Land Surface Effects on Summertime Clouds and Moist Convection in Houston Under Different Synoptic Conditions

Urban landscapes modify cloud formation and convection through complex thermodynamic and aerodynamic processes; however, their influence under different synoptic regimes remains poorly understood. This study investigates the impact of the Houston metropolitan area on summertime cloud cover and convective cell characteristics using a combination of satellite observations, radar data, and high-resolution process-based modeling. We isolate urban effects by comparing model simulations with realistic urban land cover against hypothetical scenarios where all urban areas are replaced by rural vegetation. Results reveal that Houston's land cover consistently enhances cloud fraction and convective activity relative to surrounding rural areas, altered by large-scale meteorological forcing. Under weakly forced conditions, enhanced surface heat flux primarily contributes to driving low-level convergence and vertical ascent, leading to over 8% cloud fraction increase between 2 and 6 km over the city. Under strongly forced conditions, urban influences manifest differently depending on cloud type. For non-convective clouds, the city acts as a barrier that decelerates and lifts moist southerly inflow, increasing low cloud cover over the urban core, while decreasing it downwind the city. For convective clouds, both synoptic ascent and urbanization modulate moisture redistribution and cloud structure, producing modest cloud enhancement over the city and slight suppression over the downwind area. Urbanization exerts small changes in the intensity of the convective cells; however, it significantly decreases their duration and traveling distance. This work highlights the importance of accounting for land surface heterogeneity in modeling clouds and precipitation and demonstrates that urban impacts on clouds are highly regime-dependent.

Liu, Ye [Pacific Northwest National Laboratory (PN↗
DOE OSTI2025

Technical, Economic, Energetic, and Environmental Evaluation of Pretreatment Strategies for Scaling Control in Brackish Water Desalination Brine Treatment

Effective pretreatment is essential for achieving long-term stable operation and high water recovery during the desalination of alternative waters. This study developed a process modeling approach for technical, economic, energetic, and environmental assessments of pretreatment technologies to identify the impacts of each technology treating brackish water desalination brine with high scaling propensity. The model simulations evaluated individual pretreatment technologies, including chemical softening (CS), chemical coagulation (CC), electrocoagulation (EC), and ion exchange (IX). In addition, combinations of these pretreatment technologies aiming at the effective reduction of key scaling constituents such as hardness and silica were investigated. The three evaluation parameters in this assessment consist of levelized cost of water (LCOW, $/m 3 ), specific energy consumption and cumulative energy demand (SEC|CED, kWh/m 3 ), and carbon dioxide emissions (CO 2 , kg CO 2-eq /m 3 ). The case study evaluated in this work was the desalination brine from the Kay Bailey Hutchison Desalination Plant (KBHDP) with a total dissolved solids (TDS) concentration of 11,000 mg/L and rich in hardness and silica. The evaluation of individual pretreatment units from the highest to lowest LCOW, SEC|CED, and CO 2 emissions in the KBHDP brine was IX > CS > EC > CC, CS > IX > EC > CC, and CC > CS > EC > IX, respectively. In the case of pretreatment combinations for the KBHDP, the EC + IX treatment combination was shown to be the best in terms of the LCOW and CO 2 emissions. The modeling and evaluation of these pretreatment units provide valuable guidance on the selection of cost-effective, energy-efficient, and environmentally sustainable pretreatment technologies tailored to desalination brine applications for minimal- or zero-liquid discharge.

14 SOLAR ENERGY↗
DOE OSTI2024

Core-Shell Oxidative Aromatization Catalysts for Single Step Liquefaction of Distributed Shale Gas (Final Technical Report)

The objective of this project was to design and demonstrate a core-shell structured multifunctional catalyst to convert the light (dry) components of shale gas into liquid aromatic compounds (primarily benzene and toluene) in a single step. Operated in a modular oxidative aromatization system (OAS) under a cyclic redox scheme, the novel catalyst and process can significantly improve the value and transportability of distributed shale gas. Since the project started, each quarter addressed a different set of tasks related to the completion of the milestone detailed in the project award. The yearly summaries of these tasks are summarized below: Q1-Q4: • Conducted project planning and literature search. • Investigated a number of SHC redox catalysts using thermogravimetric analysis and fixed-bed reactor experiments. • Initiated process modeling towards generating two process models for the methane DHA base case and OAS process. • Developed DHA catalysts capable of producing >500 g/kg-cat-hr aromatics at 80% or greater aromatics selectivity at 700°C. Q5-Q8: • Developed alternative approaches with sequential bed configurations to enhance the aromatic yields based on OCM+DHA • Improved the zeolite synthesis efficiency by using the microwave-assisted technique and investigated the synthesis conditions on the zeolite yield, crystalline structure and morphology • Constructed a set of Aspen Plus process models with significant energy savings for OAS as compared to the base case non-oxidative DHA. • Adapted conventional hydrothermal method to be applicable to the microwave synthesizer unit for more efficient catalyst synthesis. • Studied the structure of the OCM catalyst and the dispersion of the carbonate in the redox reactions and in methane flow with Raman Spectroscopy. Q9-Q12: • Scaled up the catalyst synthesis with the microwave synthesis method. Based on its performance, procedural characterizations and catalytic performance testing were further conducted for the new microwave synthesized catalysts with the newly-developed product analysis procedure. • Developed the reaction system setup for the C2-DHA or OCM+DHA reaction product and achieved a better product collection-analysis method for the aromatic products with an improved carbon balance. The product from the OCM reaction exhibited complicated effects on the DHA catalyst. • Conducted additional OCM catalyst characterization using Near Ambient Pressure X-ray Photoelectron Spectroscopy and in situ Raman characterization • Validated the significant energy savings for OAS as compared to the base case non-oxidative DHA. Successfully set up the simulation model for the OCM+DHA+SHC reaction system based on the updated experimental results from NCSU. Q13-End of project: • Synthesized new zeolite catalysts by the microwave method, conducted characterizations (XRD, SEM, and TEM) and catalytic behavior testing. • Explored the “wet” C 2 H 6 and C 2 H 4 DHA reactions with using steam co-feed. A subsequent reduction as the regeneration step can regenerate the DHA catalyst and recover 99% activity of the fresh performance. • Achieved a 15.3% single-pass aromatic yield from methane by rationally combining the OCM and DHA at different temperatures. • Conducted a 105-hour stability test with an improved regeneration procedure, with an average aromatic yield of 13.8%. • Developed new catalyst and achieved a record-high 23.2% yield.

03 NATURAL GAS↗
DOE OSTI2022

Evidence of Quadrupole and Octupole Deformations in Zr 96 + Zr 96 and Ru 96 + Ru 96 Collisions at Ultrarelativistic Energies

In the hydrodynamic model description of heavy-ion collisions, the elliptic flow v 2 and triangular flow v 3 are sensitive to the quadrupole deformation β 2 and octupole deformation β 3 of the colliding nuclei. The relations between v n and β n have recently been clarified and were found to follow a simple parametric form. The STAR Collaboration has just published precision v n data from isobaric Ru 96 + Ru 96 and Zr 96 + Zr 96 collisions, where they observe large differences in central collisions v 2 , Ru > v 2 , Zr and v 3 , Ru < v 3 , Zr . Using a transport model simulation, we show that these orderings are a natural consequence of $\beta_{2,Ru} \gg \beta_{2,Zr}$ and $\beta_{3,Ru} \ll \beta_{3,Zr}$ We reproduce the centrality dependence of the v 2 ratio qualitatively and v 3 ratio quantitatively and extract values of β 2 and β 3 that are consistent with those measured at low-energy nuclear structure experiments. STAR data provide the first direct evidence of strong octupole correlations in the ground state of Zr 96 in heavy-ion collisions. Our analysis demonstrates that flow measurements in high-energy, heavy-ion collisions, especially using isobaric systems, are a new precision tool to study nuclear structure physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Reliability estimation of an advanced nuclear fuel using coupled active learning, multifidelity modeling, and subset simulation

Tristructural isotropic (TRISO)-coated particle fuel is a robust nuclear fuel and determining its reliability is critical for the success of advanced nuclear technologies. However, TRISO failure probabilities are small and the associated computational models are expensive. Here we used coupled active learning, multifidelity modeling, and subset simulation to estimate the failure probabilities of TRISO fuels using several 1D and 2D models. With multifidelity modeling, we replaced expensive high-fidelity (HF) model evaluations with information fusion from two low-fidelity (LF) models. For the 1D TRISO models, we considered three multifidelity modeling strategies: only Kriging, Kriging LF prediction plus Kriging correction, and deep neural network (DNN) LF prediction plus Kriging correction. While the results across these multifidelity modeling strategies compared satisfactorily, strategies employing information fusion from two LF models called the HF model least often. Next, for the 2D TRISO model, we considered two multifidelity modeling strategies: DNN LF prediction plus Kriging correction (data-driven) and 1D TRISO LF prediction plus Kriging correction (physics-based). The physics-based strategy, as expected, consistently required the fewest calls to the HF model. However, the data-driven strategy had a lower overall simulation time since the DNN predictions are instantaneous, and the 1D TRISO model requires a non-negligible simulation time.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Networked HIL Simulation System for Modeling Large-scale Power Systems

This paper presents a network hardware-in-the-loop (HIL) simulation system for modeling large-scale power systems. Researchers have developed many HIL test systems for power systems in recent years. Those test systems can model both microsecond-level dynamic responses of power electronic systems and millisecond-level transients of transmission and distribution grids. By integrating individual HIL test systems into a network of HIL test systems, we can create large-scale power grid digital twins with flexible structures at required modeling resolution that fits for a wide range of system operating conditions. This will not only significantly reduce the need for field tests when developing new technologies but also greatly shorten the model development cycle. In this paper, we present a networked OPAL-RT based HIL test system for developing transmission-distribution coordinative Volt-VAR regulation technologies as an example to illustrate system setups, communication requirements among different HIL test systems, and connection mechanisms. Impacts of communication delays, information exchange cycles, and computing delays are illustrated. Simulation results show that the performance of a networked HIL test system is satisfactory.

digital twin, distribution system, co-simulation, ↗
DOE OSTI2023

Bayesian inference of in-medium baryon-baryon scattering cross sections from HADES proton flow data

Within a Bayesian statistical framework using a Gaussian Process emulator for an isospin-dependent Boltzmann-Uehling-Uhlenbeck (IBUU) transport model simulator of heavy-ion reactions at intermediate energies, we infer from the HADES proton flow data the posterior probability distribution functions of in-medium baryon-baryon scattering cross section modification factor X with respect to free-space and the corresponding incompressibility K of nuclear matter as well as their correlation function. In conclusion, the mean value of X is found to be $X$ = $1.32_ {+0.28} ^ {-0.40} $ at 68% confidence level assuming the nuclear incompressibility K will not exceed 400 MeV, providing circumstantial evidence for enhanced baryon-braryon scattering cross sections in hot and dense nuclear matter.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Parametric study and speciation analysis of rare earth precipitation using oxalic acid in a chloride solution system

Oxalic acid precipitation is a common step in the purification of rare earth elements (REE) from a concentrated pregnant leach solution (PLS). However, the presence of contaminants such as Al, Fe, and Ca in given amounts decreases the REE precipitation efficiency and product purity while also increasing the amount of oxalic acid needed to maximize recovery. As such, a statistically designed test program was performed to identify the optimal conditions necessary for a relatively low REE content PLS containing elevated concentrations of contaminant ions. The performance objective was maximization of REE precipitation efficiency while minimizing the oxalic acid dosage. A central composite design was utilized to quantify performance impacts and identify the ultimate set of parameter values for oxalic acid dosage, Fe(III) contamination concentration, solution pH, and reaction temperature. The resultant model suggested that oxalic acid dosage and reaction pH are the most significant factors for the REE precipitation efficiency, followed by the interaction of oxalic dosage and Fe concentration. Test results indicate that increasing the oxalic acid concentration from 0 g/L to 80 g/L improved the REE precipitation efficiency from approximately 4.2% to 95.0%. Furthermore, raising the solution pH from 0.5 to 2.5 considerably enhanced the precipitation efficiency from 0.0% to 98.9%. A solution temperature elevation decreased REE recovery, which indicated an exothermic reaction between REEs and oxalate anions. Finally, a high level of Fe contamination adversely impacted REE precipitation efficiency. Here, to further the understanding of the REE-oxalate system, a fundamental solution chemistry study was performed using the equilibrium constants of the reactions. The study resulted in the development of oxalate speciation diagrams and provided an analysis of the REE precipitation characteristics at various oxalate anion concentrations and Fe(III) contamination levels using MINTEQ software. The dominant Fe(III) species in the solution system were found to be Fe-(C 2 O 4 ) 3 3- , Fe-(C 2 O 4 ) 2- , and Fe-(C 2 O 4 ) + , which consume the majority of the oxalate anions. The simulated model was found to be in agreement with the experimental findings and helped to explain the adverse impact of increased iron concentrations on REE precipitation efficiency.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗