Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “joint state estimation”

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 127 records · Page 7

Distributed Optimization for Nonrigid Nano-Tomography

Resolution level and reconstruction quality in nano-computed tomography (nano-CT) are in part limited by the stability of microscopes, because the magnitude of mechanical vibrations during scanning becomes comparable to the imaging resolution, and the ability of the samples to resist radiation induced deformations during data acquisition. In such cases, there is no incentive in recovering the sample state at different time steps like in time-resolved reconstruction methods, but instead the goal is to retrieve a single reconstruction at the highest possible spatial resolution and without any imaging artifacts. Here we propose a distributed optimization solver for tomographic imaging of samples at the nanoscale. Our approach solves the tomography problem jointly with projection data alignment, nonrigid sample deformation correction, and regularization. Projection data consistency is regulated by dense optical flow estimated by Farneback's algorithm, leading to sharp sample reconstructions with less artifacts. Synthetic data tests show robustness of the method to Poisson and low-frequency background noise. We accelerated the solver on multi-GPU systems and validated the method on three nano-imaging experimental data sets.

97 MATHEMATICS AND COMPUTING↗

Feedstock-efficient conversion through hydrogen and formate-driven metabolism in Escherichia coli

Product yields for biomanufacturing processes are often constrained by the tight coupling of cellular energy generation and carbon metabolism in sugar-based fermentation systems. To overcome this limitation, we engineered Escherichia coli to utilize hydrogen gas (H 2 ) and formate (HCOO - ) as alternative sources of energy and reducing equivalents, thereby decoupling energy generation from carbon metabolism. This approach enabled precise suppression of decarboxylative oxidation during acetate growth, with 86.6 ± 1.6 % of electrons from hydrogen gas (via soluble hydrogenase from Cupriavidus necator H16) and 98.4 ± 3.6 % of electrons from formate (via formate dehydrogenase from Pseudomonas sp. 101) offsetting acetate oxidation. Hydrogen gas supplementation led to a titratable and stoichiometric reduction in CO 2 evolution in acetate-fed cultures. Metabolomic analysis suggests that this metabolic decoupling redirects carbon flux through the glyoxylate shunt, partially bypassing two decarboxylative steps in the TCA cycle. Here, we demonstrated the utility of this strategy by applying it to mevalonate biosynthesis, where formate supplementation during glucose fermentation increased titers by 57.6 % in our best-performing strain. Flux balance analysis further estimated that 99.0 ± 2.8 % of electrons from formate were used to enhance mevalonate production. These findings highlight a broadly applicable strategy for enhancing biomanufacturing efficiency by leveraging external reducing power to optimize feedstock and energy use.

Biomanufacturing↗

Response of joint dominated space structures

An efficient linearization method is presented for calculating the transient response of nonlinear systems due to initial disturbances. The method is an extension of the describing function approach in which the steady state response of the system is calculated by representing the nonlinear element, typically joints in the case of space structures, by impedances which are functions of the amplitude of response. Thus, the problem of solving the differential equation for the steady state response becomes one of solving a set of nonlinear algebraic equations involving the steady state amplitudes and phases of the system. It is shown that for the transient case the steady state impedances can be averaged over the range of responses in order to provide equivalent values of stiffness and damping that, for a given set of initial displacements, may be treated as being constant for purposes of calculating system response. Single degree of freedom system are used to demonstrate the method and to develop an approach for optimizing the joint's characteristics so as to minimize transient response times. The use of this method for response estimation and optimization in multiple degree of freedom systems is investigated.

Source record↗

The Determination of the Rotational State and Interior Structure of Venus with VERITAS

Understanding the processes that led Venus to its current state and will drive its future evolution is a major objective of the next generation of orbiters. In this work we analyze the retrieval of the spin vector, the tidal response and the moment of inertia of Venus with VERITAS, a NASA Discovery-class mission. By simulating a systematic joint analysis of Doppler tracking data and tie points provided by the onboard synthetic aperture radar we show that VERITAS will provide accuracies (3σ) in the estimates of the tidal Love number k_2 to 4.6×10^(-4), its tidal phase lag to 0.05°, and the moment of inertia factor to 9.8×10^(-4) (0.3% of the expected value). Applying these results to recent models of the Venus interior, we show that VERITAS will provide much improved constraints on the interior structure of the planet.

G Cascioli↗

Elucidating microstructural evolution and hardness variation across friction self-piercing riveted Al-7055 using synchrotron X-ray scattering and advanced microscopy techniques

Friction self-piercing riveting (FSPR) is a unique hybrid joining technique that combines the advantages of mechanical interlocking, frictional heat, and solid-state joining (if metallurgically compatible) to produce crack free joints in high strength and/or low-ductility alloys at room temperature. Here, in the current study, Al-7055 sheets were joined using FSPR for lightweight automotive applications and significant microhardness variations were observed across the joint cross-section. A detailed microstructural characterization at multiple length scales was carried out using advanced electron microscopy and X-ray scattering techniques to provide a fundamental understanding of the process-structure-property relationships. The relative contributions of microstructural characteristics at various length scales (i.e., grain size, dislocation density, solute concentration, precipitate nature) to strengthening were estimated using existent formulations (i.e., Hall-Petch, Taylor, precipitate bypass/shear equations) and correlated to the observed microhardness values across different regions. Small-angle X-ray scattering and scanning transmission electron microscopy revealed significant changes in the size and volume fraction of precipitate species, i.e., GP-I Zones, η′, and Mg/Zn solute co-clusters, depending on the process region. It was observed that the dissolution of the small η′/GP-I zones (T ∼ 150–200 °C) in the heat-affected zone were the key reason for the hardness drop. Further, it was shown that solid-solution, dislocation, grain size and solute co-cluster strengthening played a key role in the thermo-mechanically affected zone and grain-refined zone (GRZ). Finally, these observations were leveraged along with the Zener-Holloman relationship and grain size in the GRZ to estimate the peak joining temperature of the GRZ (∼ 350 °C) near the steel rivet.

aluminum 7xxx alloy↗

Efficient Estimation of Pauli Observables by Derandomization

In this work, we consider the problem of jointly estimating expectation values of many Pauli observables, a crucial subroutine in variational quantum algorithms. Starting with randomized measurements, we propose an efficient derandomization procedure that iteratively replaces random single-qubit measurements by fixed Pauli measurements; the resulting deterministic measurement procedure is guaranteed to perform at least as well as the randomized one. In particular, for estimating any L low-weight Pauli observables, a deterministic measurement on only of order log(L) copies of a quantum state suffices. In some cases, for example, when some of the Pauli observables have high weight, the derandomized procedure is substantially better than the randomized one. Specifically, numerical experiments highlight the advantages of our derandomized protocol over various previous methods for estimating the ground-state energies of small molecules.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Three Dimensional Numerical Simulation and Characterization of Crack Growth in the Weld Region of a Friction Stir Welded Structure

Friction stir welding (FSW) fabrication technology is being adopted in aerospace applications. The use of this technology can reduce production cost, lead-times, reduce structural weight and need for fasteners and lap joints, which are typically the primary locations of crack initiation and multi-site fatigue damage in aerospace structures. FSW is a solid state welding process that is well-suited for joining aluminum alloy components; however, the process introduces residual stresses (both tensile and compressive) in joined components. The propagation of fatigue cracks in a residual stress field and the resulting redistribution of the residual stress field and its effect on crack closure have to be estimated. To insure the safe insertion of complex integral structures, an accurate understanding of the fatigue crack growth behavior and the complex crack path process must be understood. A life prediction methodology for fatigue crack growth through the weld under the influence of residual stresses in aluminum alloy structures fabricated using FSW will be detailed. The effects and significance of the magnitude of residual stress at a crack tip on the estimated crack tip driving force are highlighted. The location of the crack tip relative to the FSW and the effect of microstructure on fatigue crack growth are considered. A damage tolerant life prediction methodology accounting for microstructural variation in the weld zone and residual stress field will lead to the design of lighter and more reliable aerospace structures

Seshadri, Banavara R.↗

HLPW-4/GMGW-3: Mesh Adaptation for RANS Technology Focus Group Workshop Summary

The AIAA workshop series helps guide the Computational Fluid Dynamics (CFD) community through publicly available test cases and impartial evaluation to establish the state-of-the-art in applied CFD. A summary of the Mesh Adaptation Technical Focus Group submissions for the joint Fourth CFD High Lift Prediction and Third Geometry & Mesh Generation Workshop is provided. The Common Research Model High-Lift variant (CRM-HL) is the target of this investigation. A 2D high-lift airfoil solution verification case was extracted from the 3D CRM-HL. Low variation is shown between 2D solutions with mesh adaptation to control estimated solution interpolation and output errors. As expected, variation for these mesh adaptation techniques is larger for the 3D CRM-HL than the 2D airfoil. More importantly, the 3D mesh adapted results have a lower variation than expert-crafted meshes of equal or larger size. Suspected multiple solutions to Reynolds-averaged Navier-Stokes with the Spalart-Allmaras turbulence model appear to resolve into a common solution with mesh resolution.

HLPW-4↗

Comparison of Raindrop Size Distribution between NASA’s S-Band Polarimetric Radar and Two-Dimensional Video Disdrometers

The ground-based-radar-derived raindrop size distribution (DSD) parameters—mass-weighted drop diameter Dmass and normalized intercept parameter NW—are the sole resource for direct validation of the National Aeronautics and Space Administration (NASA) Global Precipitation Measurement (GPM) mission Core Observatory satellite-based retrieved DSD. Both Dmass and NW are obtained from radar-measured reflectivity ZH and differential reflectivity ZDR through empirical relationships. This study uses existing relationships that were determined for the GPM ground validation (GV) program and directly compares the NASA S-band polarimetric radar (NPOL) observables of ZH and ZDR and derivedDmass and NW with those calculated by two-dimensional video disdrometer (2DVD). The joint NPOL and 2DVD datasets were acquired during three GPM GV field campaigns conducted in eastern Iowa, southern Appalachia, and western Washington State. The comparative study quantifies the level of agreement for ZH, ZDR, Dmass, and log(NW) at an optimum distance (15–40 km) from the radar as well as at distances greater than 60 km from radar and over mountainous terrain. Interestingly, roughly 10%–15% of the NPOL ZH–ZDR pairs were well outside the envelope of 2DVD-estimated ZH–ZDR pairs. The exclusion of these pairs improved the comparisons noticeably.

Ali Tokay↗

Development of Digital Twin Predictive Model for PWR Components: Updates on Multi Times Series Temperature Prediction Using Recurrent Neural Network, DMW Fatigue Tests, System Level Thermal-Mechanical-Stress Analysis

The long-term operation (LTO) of nuclear power plant (NPP) beyond their original design life of 40 years, can lead to more material damage associated with cyclic fatigue under thermal-mechanical loading cycles and associated long-term exposure of reactor material to the deleterious reactor-coolant environments. However, under this LTO condition the reactor components can still safely operate but may require more frequent Nondestructive Evaluation (NDE) of reactor components. Frequent NDE requirement may lead to frequent shutdown of the NPP. This in turn can lead to power outage and additional NDE-inspection-cost related economic loss. The economic loss can be minimized by reducing uncertainty in life estimation of safety-critical pressure boundary components and by implementing more digital approach such as by using upcoming digital-twin (DT) technology for predicting the structural states (e.g., time and location dependent inside/outside thickness temperature, stress, strain, plastic deformation, etc.) and associated fatigue life of a component in real time. Towards this goal Argonne National Laboratory (ANL) with the sponsorship of DOE Light Water Reactor Sustainability (LWRS) program is working on the development of a DT framework that can be used for real time environmental fatigue prediction of reactor components. The DT framework is based on limited experiment-data, Artificial-intelligence (AI) – Machine-Learning (ML) - Deep-Learning (DL) based techniques and Multiphysics-computational-mechanics such as finite element (FE) based modeling tools. Towards this overall goal, following are some of the major contributions made during the FY21: 1) Multiple 82/182 dissimilar metal weld (DMW) specimens (both solid-weld and joint-weld representing the actual reactor multi-metal nozzles) were fatigue tested. The resulting fatigue lives were compared to the NUREG-6909 based best-fit and design fatigue curves. Additionally, the results of 52/152 DMW fatigue specimens (which were recently tested at Republic of Korea under the sponsorship of International Nuclear Energy Research Initiative - INERI program) were compared to the NUREG-6909 based best-fit and design fatigue curves. From the comparison of 82/182 and 52/152 DMW test data with NUREG-6909 best-fit curve, most of the reported test data fall way away from the NUREG-6909 suggested best-fit or mean curve. The NUREG-6909 suggested best-fit curve is the best-fit curve of austenitic stainless steel and due to lack of enough data on Nickel-based welds, this is currently being used for predicting the life of Nickel-alloy-based welded components. However, the above observation may require higher scaling factor (e.g., ASME suggested factor of 20 on cycles rather than the current NUREG-6909 suggested factor of 12 on cycles) for scaling the austenitic-stainless-steel best-fit-curve for estimating the design or safe-life of a welded component. Accordingly, for example, if a DMW component experience a strain amplitude of 0.6% the PWR-water life of the component would be 52 cycles instead of 85 cycles. However, more DMW tests are required to further ascertain the above-mentioned observations. 2) A system level CAD and finite element model were developed which consists of reactor pressure vessel (RPV), part of steam generator (SG), part of pressurizer (PRZ), hot leg (HL), and surge line (SL). This is with detailed nozzle geometry and thermal-mechanical material properties of different metals to simulate realistic thermal-mechanical stress under connected system global thermal-mechanical boundary conditions. 3) Different system level heat transfer analyses were performed with estimation of relevant heat transfer coefficients. The resulting data were used in subsequent system level thermal-mechanical stress analysis and for generating spatial-temporal training and validation data for a system level digital-twin based temperature predictor. Transient heat transfer analyses were performed considering thermal boundary condition under design-basis (DB) loading and EDF (Électricité de France) data-based grid-load-following (EDF-GLF) loading cycles. 4) System level thermal-mechanical stress analysis was performed for identifying damage-prone hotspots and for future extension of the model for cyclic state prediction. From the system-level model simulation under DB loading cycle it is found that HL and the SL nozzle that connect to the HL can experience significant stress and strain and could be one of the weakest links in the overall reactor coolant system (RCS). 5) An AI/ML based DT model was developed for multi-time-series temperature prediction at any inside/outside thickness locations of PWR pressure boundary components. This is by using Recurrent-neural-network (RNN) and keras machine learning libraries. The RNN model was validated against two laboratory test-based data sets with one obtained through ANL’s in-air fatigue test system and other through PWR-water test loop. The experimentally validated DT model further validated against FE model results to predict thermal scarification related spatialtemporal temperatures at random locations of a component. The well validated DT model was then used for demonstrating spatial-temporal temperature prediction under 100+ years of reactor operation subjected to combined DB, EDF-GLF and randomized grid-load-following (RANDOMGLF) loading Cycles. The expert-elicitation DT model framework was developed assuming field/input/process measurements can be available from a few existing plant sensors and can readily be used by the NPP operators. The above temperature prediction model will feed to the next-step stress analysis model based on which the life of a component can be predicted in realtime, which is one of our future works.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Latent Heating from TRMM Satellite Measurements

Rainfall production is a fundamental process within the Earth;s hydrological cycle because it represents both a principal forcing term in surface water budgets, and its energetics corollary, latent heating, is the principal source of atmospheric diabatic heating. Latent heat release itself is a consequence of phase changes between the vapor, liquid, and frozen states of water. The properties of the vertical distribution of latent heat release modulate large-scale meridional and zonal circulations with the Tropics - as well as modify the energetic efficiencies of mid-latitude weather systems. This paper highlights the retrieval of observatory, which was launched in November 1997 as a joint American-Japanese space endeavor. Since then, TRMM measurements have been providing an accurate four-dimensional amount of rainfall over the global Tropics and sub-tropics - information which can be used to estimate the spacetime structure of latent heating across the Earth's low latitudes. A set of algorithm methodologies has and continues to be developed to estimate latent heating based on rain rate profile retrievals obtained from TRMM measurements. These algorithms are briefly described followed by a discussion of the foremost latent heating products that can be generate from them. The investigation then provides an overview of how TRMM-derived latent heating information is currently being used in conjunction with global weather and climate models, concluding with remarks intended to stimulate further research on latent heating retrieval from satellites.

Tao, W.-K.↗

EVA Task and 3D Pose Recognition from Video

Extravehicular Activity (EVA) has been known to involve potential risks of biomechanical stresses and injuries to crewmembers. Gathering of EVA motion patterns is necessary for risk analysis and mitigation. However, many existing techniques, such as motion capture systems, are not only cost-prohibitive but are impractical for retrospective analysis of past missions. In this work, a software tool was developed, which can estimate the 3D poses of a spacesuit from photographs or videos, without using special sensors or equipment. The tool is based on the state-of-the-art artificial intelligence and machine learning (AI/ML) system, which was trained by studying and capturing motion patterns of past and current spacesuit test data. The AI/ML tool was further enhanced using synthetically generated data, in which the suit postures, backgrounds, camera angles and illumination conditions were parametrically adjusted and rendered for training. The tool, incorporated the methodologies of Convolutional Neural Network (CNN), was trained, and tested in the cloud computing environment. The trained model was then applied on new imagery and video to extract estimated joint positions and suit outlines. The joint positions were further processed to capture activity (“digging”), pose labels (“bending”), and other useful downstream information. The model performance on new imagery and video was successfully assessed for accuracy and reliability. This AI/ML based posture recognition tool thus allows for the quantification of injury risk and task performance characterization for both current and past missions and training, which can immensely help to improve EVA task and suit design.

Kyung Han Kim↗

Precise cosmological constraints from BOSS galaxy clustering with a simulation-based emulator of the wavelet scattering transform

For this study, we perform a reanalysis of the BOSS CMASS DR12 galaxy dataset using a simulation-based emulator for the wavelet scattering transform (WST) coefficients. Moving beyond our previous works, which laid the foundation for the first galaxy clustering application of this estimator, we construct a neural net-based emulator for the cosmological dependence of the WST coefficients and the 2-point correlation function multipoles, trained from the state-of-the-art suite of abacussummit simulations combined with a flexible halo occupation distribution (HOD) galaxy model. In order to confirm the accuracy of our pipeline, we subject it to a series of thorough internal and external mock parameter recovery tests, before applying it to reanalyze the CMASS observations in the redshift range 0.46 < z < 0.57. We find that a joint WST+2-point correlation function likelihood analysis allows us to obtain marginalized 1⁢σ errors on the Λ⁢ CDM parameters that are tighter by a factor of 2.5–6, compared to the 2-point correlation function, and by a factor of 1.4–2.5 compared to the WST-only results. This corresponds to a competitive 0.9%, 2.3% and 1% level of determination for parameters ω c , ⁢σ 8 &n s , respectively, and also to a 0.7% and 2.5% constraint on derived parameters h and ƒ⁡(z)⁢⁢σ 8 ⁡(z), in agreement with the Planck 2018 results. Our results reaffirm the constraining power of the WST and highlight the exciting prospect of employing higher-order statistics in order to fully exploit the power of upcoming stage-IV spectroscopic observations.

79 ASTRONOMY AND ASTROPHYSICS↗

A two-component rain model for the prediction of attenuation and diversity improvement

A new model was developed to predict attenuation statistics for a single Earth-satellite or terrestrial propagation path. The model was extended to provide predictions of the joint occurrences of specified or higher attenuation values on two closely spaced Earth-satellite paths. The joint statistics provide the information required to obtain diversity gain or diversity advantage estimates. The new model is meteorologically based. It was tested against available Earth-satellite beacon observations and terrestrial path measurements. The model employs the rain climate region descriptions of the Global rain model. The rms deviation between the predicted and observed attenuation values for the terrestrial path data was 35 percent, a result consistent with the expectations of the Global model when the rain rate distribution for the path is not used in the calculation. Within the United States the rms deviation between measurement and prediction was 36 percent but worldwide it was 79 percent.

Crane, R. K.↗

Short arc orbit determination for altimeter calibration and validation on TOPEX/POSEIDON

TOPEX/POSEIDON (T/P) is a joint mission of United States' National Aeronautics and Space Administration (NASA) and French Centre National d'Etudes Spatiales (CNES) design launched August 10, 1992. It carries two radar altimeters which alternately share a common antenna. There are two project designated verification sites, a NASA site off the coast at Pt. Conception, CA and a CNES site near Lampedusa Island in the Mediterranean Sea. Altimeter calibration and validation for T/P is performed over these highly instrumented sites by comparing the spacecraft's altimeter radar range to computed range based on in situ measurements which include the estimated orbit position. This paper presents selected results of orbit determination over each of these sites to support altimeter verification. A short arc orbit determination technique is used to estimate a locally accurate position determination of T/P from less than one revolution of satellite laser ranging (SLR) data. This technique is relatively insensitive to gravitational and non-gravitational force modeling errors and is demonstrated by covariance analysis and by comparison to orbits determined from longer arcs of data and other tracking data types, such as Doppler Orbitography and Radiopositioning Integrated by Satellite (DORIS) and Global Positioning System Demonstration Receiver (GPSDR) data.

Williams, B. G.↗

Comparing Anisotropic Output-Based Grid Adaptation Methods by Decomposition

Anisotropic grid adaptation is examined by decomposing the steps of flow solution, ad- joint solution, error estimation, metric construction, and simplex grid adaptation. Multiple implementations of each of these steps are evaluated by comparison to each other and expected analytic results when available. For example, grids are adapted to analytic metric fields and grid measures are computed to illustrate the properties of multiple independent implementations of grid adaptation mechanics. Different implementations of each step in the adaptation process can be evaluated in a system where the other components of the adaptive cycle are fixed. Detailed examination of these properties allows comparison of different methods to identify the current state of the art and where further development should be targeted.

Park, Michael A.↗

Investigation of Site Amplifications Using Ambient-Noise-Derived Shallow Velocity Structures Under a Dense Array in Oklahoma

The shear-wave velocity (V S ) structure plays an important role in characterizing site amplification. The Large-n Seismic Survey in Oklahoma (LASSO; 1820 stations) revealed large vertical ground-motion variability in a 25 km × 32 km area in northern Oklahoma. The LASSO array has a relatively simple and flat topography, typical in a sedimentary basin environment in the central United States. In this study, we use the dense array to investigate the velocity structure under the LASSO array and how vertical ground motions relate to the shallow-to-deep structures. We extract the fundamental-mode Rayleigh wave by cross-correlating one month of ambient noise (0.7–5 Hz). We use double-beamforming to measure the group and phase velocities and anisotropy. By jointly inverting the group and phase velocities, we obtain the V S structure. Here, we observe correlations between V S at depths of 0.1–1.5 km and vertical ground motions using sites on the stiffer Permian formations. The shallow Quaternary alluvium and terrace deposits can amplify vertical ground motions by a factor of 2–4.5 between 2 and 25 Hz and attenuate signals above 25 Hz. We use 1D V S profiles to simulate the SV-wave transfer functions. An average V S of 250 m/s in the upper 20–40 m may cause the observed amplification between 2 and 40 Hz. V S estimated by topographic slopes cannot predict the relative amplification. Our results highlight the large variability of site-dependent ground motion in a small local region and the importance of characterizing shallow structures to estimate seismic hazards. Small thickness variations of the shallow formation can significantly change the resonance amplitude and frequency, which likely reduces the coherency of the Rayleigh waves extracted from ambient noise.

Chang, Hilary [Massachusetts Inst. of Technology (↗

Adaptive Online Multivariate Signal Extraction With Locally Weighted Robust Polynomial Regression

High-frequency, multivariate data collected in real-time and used to control or make decisions regarding a process’ operation often contain some noise and outliers. Thus, a method to extract the signal is needed in order to reduce the number and magnitude of control-based adjustments that are implemented. Such a method must be (i) online, depending only on past and current observations; (ii) fast, producing a smooth value more quickly than the measurement frequency; (iii) robust, ignoring brief bursts of erroneously measured values; (iv) multivariate, ignoring observations that are jointly unusual; (v) adaptive, adjusting to periods of rapid fluctuation in the signal versus periods of stability; and (vi) purely data-driven, not incorporating any information about the process from which the data are collected. Most existing methods are only able to address a subset of these six features. Furthermore, we also require the method to be nonlinear, providing a local nonlinear estimate of the signal. In this work, we propose a novel, real-time signal extraction method based on a local, robust polynomial fit. We demonstrate the performance of our method compared to a state-of-the-art competitor through simulation. For illustration, the methodology is applied to data collected from a reverse osmosis water treatment process.

97 MATHEMATICS AND COMPUTING↗