Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Data Flow Reversal”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS1983

Interplanetary flow systems associated with cosmic ray modulation in 1977 - 1980

The hydromagnetic flow configurations associated with cosmic ray modulation in 1977 to 1980 were determined using solar wind plasma and magnetic field data from Voyagers 1 and 2 and Helios 1. The modulation was related to two types of large scale systems of flows: one containing a number of transients such as shocks, post shock flows and magnetic clouds; the other consisting primarily of a series of quasi-stationary flows following interaction regions containing a stream interface and often bounded by a forward reverse shock pair. Each of the three major episodes of cosmic ray modulation was characterized by the passage of the system of transient flows. Plateaus in the cosmic ray intensity time profile were associated with the passage of systems of corotating streams.

Burlaga, L. F.↗
NASA NTRS1990

Analytical aerodynamic model of a high alpha research vehicle wind-tunnel model

A 6 DOF analytical aerodynamic model of a high alpha research vehicle is derived. The derivation is based on wind-tunnel model data valid in the altitude-Mach flight envelope centered at 15,000 ft altitude and 0.6 Mach number with Mach range between 0.3 and 0.9. The analytical models of the aerodynamics coefficients are nonlinear functions of alpha with all control variable and other states fixed. Interpolation is required between the parameterized nonlinear functions. The lift and pitching moment coefficients have unsteady flow parts due to the time range of change of angle-of-attack (alpha dot). The analytical models are plotted and compared with their corresponding wind-tunnel data. Piloted simulated maneuvers of the wind-tunnel model are used to evaluate the analytical model. The maneuvers considered are pitch-ups, 360 degree loaded and unloaded rolls, turn reversals, split S's, and level turns. The evaluation finds that (1) the analytical model is a good representation at Mach 0.6, (2) the longitudinal part is good for the Mach range 0.3 to 0.9, and (3) the lateral part is good for Mach numbers between 0.6 and 0.9. The computer simulations show that the storage requirement of the analytical model is about one tenth that of the wind-tunnel model and it runs twice as fast.

Cao, Jichang↗
NASA NTRS1993

Simulations of diffusion-reaction equations with implications to turbulent combustion modeling

An enhanced diffusion-reaction reaction system (DRS) is proposed as a statistical model for the evolution of multiple scalars undergoing mixing and reaction in an isotropic turbulence field. The DRS model is close enough to the scalar equations in a reacting flow that other statistical models of turbulent mixing that decouple the velocity field from scalar mixing and reaction (e.g. mapping closure model, assumed-pdf models) cannot distinguish the model equations from the original equations. Numerical simulations of DRS are performed for three scalars evolving from non-premixed initial conditions. A simple one-step reversible reaction is considered. The data from the simulations are used (1) to study the effect of chemical conversion on the evolution of scalar statistics, and (2) to evaluate other models (mapping-closure model, assumed multivariate beta-pdf model).

Girimaji, Sharath S.↗
NASA NTRS1975

Calculation of three-dimensional compressible laminar and turbulent boundary layers. An implicit finite-difference procedure for solving the three-dimensional compressible laminar, transitional, and turbulent boundary-layer equations

An implicit finite-difference procedure is presented for solving the compressible three-dimensional boundary-layer equations. The method is second-order accurate, unconditionally stable (conditional stability for reverse cross flow), and efficient from the viewpoint of computer storage and processing time. The Reynolds stress terms are modeled by (1) a single-layer mixing length model and (2) a two-layer eddy viscosity model. These models, although simple in concept, accurately predicted the equilibrium turbulent flow for the conditions considered. Numerical results are compared with experimental wall and profile data for a cone at an angle of attack larger than the cone semiapex angle. These comparisons clearly indicate that the numerical procedure and turbulence models accurately predict the experimental data with as few as 21 nodal points in the plane normal to the wall boundary.

Harris, J. E.↗
NASA NTRS1995

Simultaneous prenoon and postnoon observations of three field-aligned current systems from Viking and DMSP-F7

The spatial structure of dayside large-scale field-aligned current (FAC) systems is examined by using Viking and Defense Meteorological Satellite Program-F7 (DMSP-F7) data. We focus on four events in which the satellites simultaneously observed postnoon and prenoon three FAC systems: the region 2, the region 1, and the mantle (referred to as midday region O) systems, from equatorward to poleward. These events provide the most solid evidence to date that the midday region O system is a separate and unique FAC system, and is not an extension of the region 1 system from other local times. The events are examined comprehensively by making use of a mulit-instrumental data set, which includes magnetic field, particle flux, electric field, auroral UV image data from the satellites, and the Sondrestrom convection data. The results are summarized as follows: (1) Region 2 currents flow mostly in the central plasma sheet (CPS) precipitation region, often overlapping with the boundary plasma sheet (BPD) at their poleward edge. (2) The region 1 system is located in the core part of the auroral oval and is confined in a relatively narrow range in latitude which includes the convection reversal. The low-latitude boundary layer, possibly including the outer part of the plasma sheet, and the external cusp are the major source regions of dayside region 1 currents. (2) Midday region O currents flow on open field lines and are collocated with the shear of antisunward convection flows with velocites decreasing poleward. On the basis of these results we support the view that both prenoon and postnoon current systems consist of the three-sheet structure when the disctortion ofthe convection pattern associated with interplanetary magnetic field (IMF) B(sub Y) is small and both morningside and eveningside convection cells are crescent-shaped. We also propose that the midday region O and a part of the region 1 systems are closely coupled to the same source.

Ohtani, S.↗
DOE OSTI2021

Identification of crystal plasticity model parameters by multi-objective optimization integrating microstructural evolution and mechanical data

Crystal plasticity models evolve a polycrystalline yield surface using meso-scale descriptions of deformation mechanisms. The activation of deformation mechanisms is governed by crystallography and a set of model parameters, which are typically calibrated through the fitting of mechanical data such as stress–strain curves and elastic lattice strains. Microstructural data such as phase fractions and texture evolution are used for verifying crystal plasticity parameters. In this study, we use a multi-objective genetic algorithm to identify hardening parameters from flow stress curves with an option to incorporate texture into the optimization approach. Robust, generalized objective functions are developed and used to identify sets of parameters pertaining to dislocation density-based hardening laws in visco-plastic and elasto-plastic self-consistent (VPSC and EPSC) homogenization models. First, the parameters are identified for pure Nb directly from texture using an objective function based on generalized spherical harmonics. Since texture evolution is driven by the relative contribution of active slip systems, the parameters governing the evolution of slip resistance ratios can be recovered from fitting discrete textures at a series of strains. Next, a comprehensive set of load reversal data for dual phase (DP) 780 steel is used to fit a hardening law and a back-stress law in EPSC. Finally, parameters pertaining to a complex hardening law for the evolution of slip and twinning in pure α-Ti are identified. Remarkably, using texture as an objective in combination with stress–strain objectives constrains the model of Ti to fully reproduce not only stress–strain and texture evolution but also hierarchical twinning measurements as a function of initial grain size and texture. Furthermore, given an appropriate model fit to representative experimental texture evolution, underlying twin volume fractions contributing to texture evolution can be predicted.

42 ENGINEERING↗
DOE OSTI2025

A DATA EFFICIENT SPARSE MODELING FRAMEWORK FOR POWER ESTIMATION IN WATER TREATMENT SENSING OPERATIONS

With increasing freshwater scarcity, advanced process design mechanisms such as Closed-Circuit Reverse Osmosis (CCRO) and Digital/Physical Twin systems are gaining traction in water treatment and reuse operations. While digital and physical twin models enable improved system insight and control, their development is often expensive and computationally intensive, requiring large volumes of synthetic or experimental data to characterize underlying process dynamics. This work introduces a sparse surrogate modeling framework to estimate power consumption from measured flow and pressure variables, along with their nonlinear polynomial and interaction expansions. To ensure model reliability and reduce overfitting, a two-stage pipeline is proposed. First, a dynamic data filtering algorithm is employed to remove uninformative observations and transient operational states. Second, a sparse penalized regression technique is applied to select a minimal set of parsimonious features. The proposed model achieves high sparsity, retaining only 7 out of 34 candidate features (≈79.41% sparsity) while delivering a root mean square error (RMSE) of 0.072 on the test dataset.

Mukherjee, Subrata [ORNL] (ORCID:0000000309930338)↗
NASA NTRS2018

Numerical Simulation of a Simplified High-Lift CRM Configuration Embedded with Fluidic Actuators

Numerical simulations have been performed for a simplified high-lift configuration that is representative of a modern transport airplane. This configuration includes a leading-edge slat, fuselage, wing, nacelle-pylon and a simple hinged flap. The suction surface of the flap is embedded with multiple rows of fluidic actuators to reduce the extent of reversed flow regions and improve the aerodynamic performance of the configuration with flap in a deployed state. In the current paper, a Lattice Boltzmann Method based high-fidelity computational fluid dynamics (CFD) code, known as PowerFLOW® is used to simulate the entire flow field associated with this configuration, including the flow inside the actuators. A fully compressible version of the PowerFLOW® code that has been validated for high speed flows is used for the present simulations to accurately represent the transonic flow regimes that are encountered in the flow field generated by the actuators operating at higher mass flow (momentum) rates required to mitigate reverse flow regions on the suction surfaces of the main wing and the flap. The numerical solutions predict the expected trends in aerodynamic forces as the actuation levels are increased. More efficient active flow control (AFC) systems and actuator arrangement for lift augmentation are emerging based on the parametric studies conducted here prior to wind tunnel tests. These numerical solutions will be compared with experimental data, once such data becomes available.

Vatsa, Veer N.↗
DOE OSTI2023

Engineering a modular 44 Ti/ 44 Sc generator: eluate evaluation in preclinical models and estimation of human radiation dosimetry

Background: 44 Sc/ 47 Sc is an attractive theranostic pair for targeted in vivo positron emission tomographic (PET) imaging and beta-particle treatment of cancer. The 44 Ti/ 44 Sc generator allows daily onsite production of this diagnostic isotope, which may provide an attractive alternative for PET facilities that lack in-house irradiation capabilities. Early animal and patient studies have demonstrated the utility of 44 Sc. In our current study, we built and evaluated a novel clinical-scale 44 Ti/ 44 Sc generator, explored the pharmacokinetic profiles of 44 ScCl 3 , [ 44 Sc]-citrate and [ 44 Sc]-NODAGA (1,4,7-triazacyclononane,1-glutaric acid-4,7-acetic acid) in naïve mice, and estimated the radiation burden of 44 ScCl 3 in humans. Methods: 44 Ti/ 44 Sc (101.2 MBq) in 6 M HCl solution was utilized to assemble a modular ZR resin containing generator. After assembly, 44 Sc was eluted with 0.05 M HCl for further PET imaging and biodistribution studies in female Swiss Webster mice. Based on the biodistribution data, absorbed doses of 44 / 47 ScCl 3 in human adults were calculated for 18 organs and tissues using the IDAC-Dose software. Results: 44 Ti in 6 M HCl was loaded onto the organic resin generator with a yield of 99.97%. After loading and initial stabilization, 44 ScCl 3 was eluted with 0.05 M HCl in typical yields of 82.9 ± 5.3% (N = 16), which was normalized to the estimated generator capacity. Estimated generator capacity was computed based on elution time interval and the total amount of 44 Ti loaded on the generator. Run in forward and reverse directions, the 44 Sc/ 44 Ti ratio from a primary column was significantly improved from 1038 ± 440 to 3557 ± 680 (Bq/Bq) when a secondary, replaceable, ZR resin cartridge was employed at the flow outlet. In vivo imaging and ex vivo distribution studies of the reversible modular generator for 44 ScCl 3 , [ 44 Sc]-citrate and [ 44 Sc]-NODAGA show that free 44 Sc remained in the circulation significantly longer than the chelated 44 Sc. The dose estimation of 44 ScCl 3 reveals that the radiation burden is 0.146 mSv/MBq for a 70 kg adult male and 0.179 mSv/MBq for a 57 kg adult female. Liver, spleen and heart wall will receive the highest absorbed dose: 0.524, 0.502, and 0.303 mGy/MBq, respectively, for the adult male. Conclusions: A clinical-scale 44 Ti/ 44 Sc generator system with a modular design was developed to supply 44 ScCl 3 in 0.05 M HCl, which is suitable for further radiolabeling and in vivo use. Our data demonstrated that free 44 ScCl 3 remained in the circulation for extended periods, which resulted in approximately 10 times greater radiation burden than stably chelated 44 Sc. Stable 44 Sc/ 47 Sc-complexation will be more favorable for in vivo use and for clinical utility.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2023

Changes within the central stalk of E. coli F 1 F o ATP synthase observed after addition of ATP

F 1 F o ATP synthase functions as a biological generator and makes a major contribution to cellular energy production. Proton flow generates rotation in the F o motor that is transferred to the F 1 motor to catalyze ATP production, with flexible F 1 /F o coupling required for efficient catalysis. F 1 F o ATP synthase can also operate in reverse, hydrolyzing ATP and pumping protons, and in bacteria this function can be regulated by an inhibitory ε subunit. Here we present cryo-EM data showing E. coli F 1 F o ATP synthase in different rotational and inhibited sub-states, observed following incubation with 10 mM MgATP. Our structures demonstrate how structural transitions within the inhibitory ε subunit induce torsional movement in the central stalk, thereby enabling its rotation within the F ο motor. This highlights the importance of the central rotor for flexible coupling of the F 1 and F o motors and provides further insight into the regulatory mechanism mediated by subunit ε.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1989

Low-latitude zonal and vertical ion drifts seen by DE 2

Horizontal and vertical ion drift data from the DE 2 spacecraft have been used to determine average zonal and vertical plasma flow (electric field) characteristics in the +/- 26-deg dip latitude region during a time of high solar activity. The 'average data' local time profile for an apex height bin centered at 400 km indicates westward plasma flow from 0600 to 1900 solar local time ((SLT) with a maximum westward velocity of 80 m/s in the early afternoon. There is a sharp change to eastward flow at approximately 1900 hours with an early evening peak of 170 m/s. A secondary nighttime maximum exists at 0430 SLT preceeding the reversal to westward flow. This profile is in good agreement with Jicamarca, Peru, radar measurements made under similar solar maximum conditions. Haramonic analysis indicates a net superrotation which is strongest at lower apex altitudes. The diurnal term is dominant, but higher order terms through the quatradiurnal are significant.

Coley, W. R.↗
NASA NTRS2000

Size Exclusion Chromatography Studies of the Initial Self-Association Steps of Chicken Egg White Lysozyme Nucleation

Nucleation is one of the least understood aspects of crystallogenesis. In the case of macromolecule nucleation, this understanding is further hampered by uncertainty over what precisely is being discussed. We define the process of solute self-association (aggregation, oligomerization, interaction, clustering, etc.) whereby n-mers (n > or = 2) having a crystallographic or nascent crystallographic arrangement leading to the critical nucleus reversibly form in the solution, to be part of the nucleation process. This reversible self-association process is a fundamental part of the nucleation process, and occurs as a function of the solute concentration. In the case of chicken egg white lysozyme, a considerable body of experimental evidence leads us to the conclusion that it also forms the crystal growth units. Size exclusion chromatography is a simple and direct method for determining the equilibrium constants for the self-association process. A Pharmacia FPLC system was used to provide accurate solution flow rates. The column, injection valve, and sample loop were all mounted within a temperature-controlled chamber. Chromatographically re-purified lysozyme was first dialyzed against the column equilibration buffer, with injection onto the column after several hours pre-incubation at the running temperature. Preliminary experiments, were carried out using a Toyopearl HW-50F column (1 x 50cm), equilibrated with 0.1 M sodium acetate, 5% sodium chloride, pH 4.6, at 15C. Protein concentrations from 0.1 to 4 mg/ml were employed (C(sub sat) = 1.2 mg/ml). The data from several different protein preparations consistently shows a progressively decreasing elution volume with increasing protein concentration, indicating that reversible self-association is occurring. The dotted line indicates the monomeric lysozyme elution volume. However, lysozyme interacts with the column matrix in these experiments, which complicates data analysis.Accordingly, we are testing silica-based HPLC columns in an effort to eliminate this problem and substantially reduce the column volume and experimental run time. The results and data analysis from these and subsequent experiments will be presented.

Ewing, Felecia↗
DOE OSTI2022

Interfacial Reactivity of Silicon Electrodes: Impact of Electrolyte Solvent and Presence of Conductive Carbon

Silicon (Si) is a promising high-capacity material for lithium-ion batteries; however, its limited reversibility hinders commercial adoption. Approaches such as particle and crystallite size reduction, introduction of conductive carbon, and use of different electrolyte solvents have been explored to overcome these electrochemical limitations. Herein, operando isothermal microcalorimetry (IMC) is used to probe the influence of silicon particle size, electrode composition, and electrolyte additives fluoroethylene carbonate and vinylene carbonate on the heat flow during silicon lithiation. In this work, the IMC data are complemented by X-ray photoelectron and Raman spectroscopies to elucidate differences in solid electrolyte interphase (SEI) composition. Nanosized (~50 nm, n-Si) and micrometer-sized (~4 μm, μ-Si) silicon electrodes are formulated with and without amorphous carbon and electrochemically lithiated in ethylene carbonate (EC), fluoroethylene carbonate (FEC), or vinylene carbonate (VC) based electrolytes. Notably, n-Si electrodes generate 53–61% more normalized heat relative to their μ-Si counterparts, consistent with increased surface area and electrode/electrolyte reactivity. Introduction of amorphous carbon significantly alters the heat flow profile where multiple exothermic peaks and increased normalized heat dissipation are observed for all electrolyte types. Notably, the VC-containing electrolyte demonstrates the greatest normalized heat dissipation of the electrode compositions tested showing as much as a 50% increase compared to the EC or FEC counterparts. The results are relevant to the understanding of silicon negative electrode function in the presence of electrolyte additives and provide insight relative to silicon containing cell reactivity and safety.

25 ENERGY STORAGE↗
DOE OSTI2024

Offshore low-level jet observations and model representation using lidar buoy data off the California coast

Abstract. Low-level jets (LLJs) occur under a variety of atmospheric conditions and influence the available wind resource for wind energy projects. In 2020, lidar-mounted buoys owned by the US Department of Energy (DOE) were deployed off the California coast in two wind energy lease areas administered by the Bureau of Ocean Energy Management: Humboldt and Morro Bay. The wind profile observations from the lidars and collocated near-surface meteorological stations (4–240 m) provide valuable year-long analyses of offshore LLJ characteristics at heights relevant to wind turbines. At Humboldt, LLJs were associated with flow reversals and north-northeasterly winds, directions that are more aligned with terrain influences than the predominant northerly flow. At Morro Bay, coastal LLJs were observed primarily during northerly flow as opposed to the predominant north-northwesterly flow. LLJs were observed more frequently in colder seasons within the lowest 250 m a.s.l. (above sea level), in contrast with the summertime occurrence of the higher-altitude California coastal jet influenced by the North Pacific High, which typically occurs at heights of 300–400 m. The lidar buoy observations also validate LLJ representation in atmospheric models that estimate potential energy yield of offshore wind farms. The European Centre for Medium-Range Weather Forecasts Reanalysis version 5 (ERA5) was unsuccessful at identifying all observed LLJs at both buoy locations within the lowest 200 m. An extension of the National Renewable Energy Laboratory (NREL) 20-year wind resource dataset for the Outer Continental Shelf off the coast of California (CA20-Ext) yielded marginally greater captures of observed LLJs using the Mellor–Yamada–Nakanishi–Niino (MYNN) planetary boundary layer (PBL) scheme than the 2023 National Offshore Wind dataset (NOW-23), which uses the Yonsei University (YSU) scheme. However, CA20-Ext also produced the most LLJ false alarms, which are instances when a model identified an LLJ but no LLJ was observed. CA20-Ext and NOW-23 exhibited a tendency to overestimate the duration of LLJ events and underestimate LLJ core heights.

17 WIND ENERGY↗
DOE OSTI2021

L-H transition trigger physics in ITER-similar plasmas with applied n = 3 magnetic perturbations

The L–H transition power threshold PLH is observed to increase with applied n = 3 resonant magnetic perturbations (RMP) in ITER-similar-shape plasmas with balanced neutral beam torque injection in DIII-D. The increase is most pronounced with added electron–cyclotron heating: PLH increases with decreasing edge plasma collisionality as PLH/PLH-08 ~ (ν*)−0.5, where PLH-08 is the 2008 ITPA multi-machine power threshold scaling. This result raises concerns for H-mode access at low edge collisionality in ITER, where RMP may have to be applied before the L–H transition to safely suppress the first edge-localized mode. Non-axisymmetric modifications with RMP include a simultaneous reduction of the radial electric field (Er) well depth and E × B shear. This can be attributed to increasing edge toroidal co-current rotation, and is consistent with substantially increased local long-wavelength turbulence (measured via beam emission spectroscopy). At high RMP perturbation strength the edge electric field Er reverses sign locally (becomes positive), with changes in dominant turbulence modes. Edge magnetic stochasticity provides an attractive explanation of the observed modifications, and the observed changes in toroidal rotation and Er are consistent with a simple fluid model describing radial electron current flow along stochastic fieldlines. The observed collisionality dependence of the L-mode edge electric field with applied RMP is also qualitatively consistent with this model. Reflectometry data indicate a significant reduction of the normalized L-mode radial density gradient a/Ln at high RMP field with simultaneous increase in radial particle flux and electron thermal flux from power balance analysis. We conjecture that the increase of PLH with RMP results from the combined effects of reduced E × B flow shear (increasing turbulent transport levels) and toroidal/poloidal flow modulation due to edge stochasticity. Initial experiments indicate that non-resonant n = 3 magnetic perturbations lead only to relatively small changes in Er, E × B shear and fluctuation characteristics, and have less impact on the L–H transition power threshold. This motivates further exploration of the RMP spectrum dependence of PLH for possible mitigation of the observed threshold increase.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2019

Comparative Study of Active Flow Control Strategies for Lift Enhancement of a Simplified High-Lift Configuration

Numerical simulations have been performed for a simplified high-lift (SHL) version of the Common Research Model (CRM) configuration, where the Fowler flaps of the conventional high-lift (CRM-HL) configuration are replaced by a set of simple hinged flaps. These hinged flaps are equipped with integrated modular active flow control (AFC) cartridges on the suction surface, and the resulting geometry is known as the CRM-SHL-AFC configuration. The main objective is to make use of AFC devices on the CRM-SHL-AFC configuration to recover the aerodynamic performance (lift) of the CRM-HL configuration. In the current paper, a Lattice Boltzmann method-based computational fluid dynamics (CFD) code, known as PowerFLOWQ® is used to simulate the entire flow field associated with the CRM-SHL-AFC configuration equipped with several different types of AFC devices. The transonic version of the PowerFLOWQ® code that has been validated for high speed flows is used to accurately simulate the flow field generated by the high-momentum actuators required to mitigate reversed flow regions on the suction surfaces of the main wing and the flap. The numerical solutions predict the expected trends in aerodynamic forces as the actuation levels are increased. More efficient AFC systems and actuator arrangements emerged based on the parametric studies performed prior to a Fall 2018 wind tunnel test. Preliminary comparisons of the numerical solutions for lift and surface pressures are presented here with the experimental data, demonstrating the usefulness of CFD for predicting the flow field and lift characteristics of AFC-enabled high-lift configurations.

Vatsa, Veer N.↗
DOE OSTI2025

Machine learning modeling and model predictive control of a closed-circuit reverse osmosis system

Closed-circuit reverse osmosis (CCRO) offers a flexible and energy-efficient alternative to conventional reverse osmosis by operating in a semi-batch mode that recycles brine, enabling higher recovery rates and reduced specific energy consumption (SEC). However, developing accurate, system-level dynamic models for CCRO remains challenging due to its nonlinear, multi-phase operation and sensitivity to variable feed water conditions. Traditional modeling approaches, such as NARMAX (nonlinear autoregressive moving average with exogenous inputs), often struggle to generalize across varying inlet feed concentrations, necessitating frequent parameter re-estimation and limiting their utility for real-time control applications. To address these limitations, we developed a long short-term memory (LSTM) neural network model trained on an extensive experimental data set from a CCRO pilot plant. The model accepts three inputs, feed flow rate, recirculation flow rate, and initial feed conductivity, and predicts three key outputs: reject conductivity, feed pump power draw, and recirculation pump power draw. We validated the LSTM model against experimental data, demonstrating its ability to distinguish between different feed conductivities and adapt to variable flow rates. Subsequently, we incorporated the LSTM model within a nonlinear model predictive control (MPC) scheme and conducted closed-loop simulations to optimize the integrated SEC (iSEC). In conclusion, the results project up to a 6% reduction in iSEC by using MPC to optimize performance over the entire experiment duration, without requiring any random excitation for data collection or parameter re-estimation.

Desalination↗
DOE OSTI2021

Transport Modeling of Kinetically Limited Microscale Extraction Systems: Droplet and Supported Liquid Membrane Separations

Microfluidic droplet extraction devices can be used to efficiently measure solute phase transfer kinetics in liquid–liquid extraction systems. In this work, a droplet extraction apparatus was constructed using off-the-shelf and custom 3D printed parts for less than 110 USD. The interfacial kinetics of the transfer of Pu(IV) between an aqueous nitric acid phase and a 5 v/v% tridodecylamine (TDA) in toluene phase were characterized by determining the extent of Pu(IV) transfer between phases at different flow rates. Forward interfacial rate constants calculated from this data using an analytical model are 5.6 × 10 –3 mm s –1 under extraction (aqueous-to-organic solute transfer) conditions and 3.3 × 10 –3 mm s –1 under stripping (organic-to-aqueous solute transfer) conditions. The analytical model is valid when the rate of solute transfer across the interface between the two liquid phases has pseudo-first-order reversible reaction kinetics and is slow relative to the rate of solute transport to the interface. The change in the outlet concentration with residence time determined using the analytical model was compared with the results of a numerical model accounting for advection/diffusion to the interface in addition to kinetically limited interfacial solute transfer. These calculations revealed that for systems with forward interfacial rate constants greater than approximately 10 –2 mm s –1 , the analytical model must be supplemented with a numerical transport model accounting for an additional advection/diffusion mass transfer resistance. This threshold corresponds to a Damköhler number of one. The measured interfacial kinetic rate constants for Pu(IV)/TDA were used in a supported liquid membrane (SLM) numerical model, which showed that the Pu(IV) phase transfer kinetics are fast enough for use in a practical separation.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗