Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamical System Scaling”

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,135 records · Page 63

Initial Study of An Effective Fast-Time Simulation Platform for Unmanned Aircraft System Traffic Management

Small Unmanned Aerial Vehicles (sUAVs), typically 55 lbs and below, are envisioned to play a major role in surveilling critical assets, collecting important information, and delivering goods. Large scale small UAV operations are expected to happen in low altitude airspace in the near future. Many static and dynamic constraints exist in low altitude airspace because of manned aircraft or helicopter activities, various wind conditions, restricted airspace, terrain and man-made buildings, and conflict-avoidance among sUAVs. High sensitivity and high maneuverability are unique characteristics of sUAVs that bring challenges to effective system evaluations and mandate such a simulation platform different from existing simulations that were built for manned air traffic system and large unmanned fixed aircraft. NASA's Unmanned aircraft system Traffic Management (UTM) research initiative focuses on enabling safe and efficient sUAV operations in the future. In order to help define requirements and policies for a safe and efficient UTM system to accommodate a large amount of sUAV operations, it is necessary to develop a fast-time simulation platform that can effectively evaluate requirements, policies, and concepts in a close-to-reality environment. This work analyzed the impacts of some key factors including aforementioned sUAV's characteristics and demonstrated the importance of these factors in a successful UTM fast-time simulation platform.

UTM↗

Plasma Surface Interactions: Predicting the Performance and Impact of Dynamic PFC Surfaces

The objective of this project is to develop, and integrate, high-performance simulation tools capable of predicting plasma-facing component (PFC) operating lifetime and the impact of the evolving surface morphology of tungsten-based PFCs on plasma contamination, including the dynamic recycling of fuel species and tritium retention, in future magnetic fusion devices. Establishing a fundamental physical understanding and developing predictive capabilities of plasma-surface interactions (PSI) requires simultaneously addressing complex and diverse physics occurring over a wide range of length (Angstroms to meters) and time (femtoseconds to years) scales, as well as integrating extensive physical processes across the plasma–surface interface. This requires development of not only detailed physics models and computational strategies at each scale, but also algorithms and methods to couple them effectively in a way that can be robustly validated. Deploying these tools requires the continued development and coupling of leadership-scale computational codes to describe the boundary plasma and the evolving PFC surface, as well as a host of simulations that bridge disparate scales to address complex physical and computational issues at the plasma–surface interface in multi-component materials systems for magnetic fusion energy development beyond ITER.

36 MATERIALS SCIENCE↗

Role of the Terrestrial Biosphere in Atmospheric Chemistry and Climate

The terrestrial biosphere–atmosphere interface provides a key chemical, biological, and physical lower boundary for the atmosphere. The presence of vegetation itself modifies the physical boundary, or the biogeophysical aspects of the system, by controlling important climate drivers such as soil moisture, light environment, and temperature. The leaf surface area of the terrestrial biosphere provides additional surface area for emissions, and it can be up to 55% of the total Earth’s surface area during the boreal summer. Vegetation also influences the biogeochemical aspects of the system by emitting a broad suite of reactive trace gases such as biogenic volatile organic compound (BVOC) emissions and climate-relevant primary biological aerosol particles (PBAP). Many of these emissions are a function of meteorological and climatological conditions at the surface, including temperature, light environment, soil moisture, and winds. Once emitted, they can be processed in the troposphere through a suite of chemical reactions. BVOC can contribute to the formation of ozone and secondary organic aerosols (SOA), and PBAP can rupture to form smaller particles with climatic relevance. These emissions and subsequent aerosol products can influence atmospheric processes that affect the surface climate, such as the attenuation of radiation, the formation of greenhouse gases such as ozone that can feedback to surface air temperature, and the alteration of clouds and subsequent precipitation. These atmospheric changes can then feedback to the land surface and emissions themselves, creating positive or negative feedback loops that can dampen or amplify the emission response. For the dominant BVOC isoprene, the feedback response to temperature can be positive or negative depending on ambient temperatures that drive isoprene emissions. The feedback response to soil moisture and precipitation can be positive, negative, or uncoupled depending on the moisture content of the soil and the total atmospheric aerosol loading. For light, the isoprene response can be positive or negative depending on the role of diffuse light. Overall, these feedbacks highlight the dynamical response of the biosphere to changing atmospheric conditions across a range of time scales, from minutes for trace gases and aerosols, to months for phenological changes, to years for land cover and land use change. Here, the dynamic aspect of this system requires us to understand, simulate, and predict the complex feedbacks between the biosphere and atmosphere and understand their role in the simulation and understanding of climate and global change. From the observational perspective, these feedbacks are challenging to identify in observations, and predictive modeling tools provide a crucial link for understanding how these feedbacks will change under warming climate scenarios.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Renewable energy integration and system operation challenge: control and optimization of millions of devices

The electric power infrastructure, originally designed and built on large-scale power plants, is evolving into a more resilient power generation and delivery system in which millions of smaller units of distributed energy generation resources units will be installed in sub-transmission and distribution networks. In order to control, manage and optimize the future grid, a hierarchical design is presented in this chapter which enables the distributed control on grid edge while inheriting the existing centralized control structure. This layered design of large-scale power system operation and control uses the following principle: reactive power control is treated as a primary control for voltage stability, and the real power control is primarily a grid-level control but can also be a supplementary control for voltage support in the case of insufficient reactive power control capacity. For the purpose of active control and operation at the distribution level, a recursive power network model is derived from nodal injection and branch power flow models. Based on the model, the proposed algorithms of hierarchical control, grid-edge inference and dynamic hosting allowance are developed and presented for multi-level controlled operation. And, a co-simulation architecture of integrated T&D system is presented to validate and demonstrate the feasibility and scalability of proposed algorithms.

Xu, Ying↗

Molecular-gas-dynamics simulations of turbulent Couette flow over a mean-free-path-scale permeable substrate

Here we report flow statistics and visualizations from molecular-gas-dynamics simulations using the direct simulation Monte Carlo (DSMC) method for turbulent Couette flow in a minimal domain where the lower wall is replaced by an idealized permeable fibrous substrate representative of thermal-protection-system materials for which the Knudsen number is O(10 -1 ). Comparisons are made with smooth-wall DSMC simulations and smooth-wall direct numerical simulations (DNS) of the Navier-Stokes equations for the same conditions. Roughness, permeability, and noncontinuum effects are assessed. In the range of Reynolds numbers considered herein, the scalings of the skin friction on the permeable substrate and of the mean flow within the substrate suggest that they are dominated by viscous effects. While the regenerative cycle characteristic of smooth-wall turbulence remains intact for all cases considered, we observe that the near-wall velocity fluctuations are modulated by the permeable substrate with a wavelength equal to the pore spacing. Additionally, the flow within the substrate shows significant rarefaction effects, resulting in an apparent permeability that is 13% larger than the intrinsic permeability. In contrast, the smooth-wall DSMC and DNS simulations exhibit remarkably good agreement for the statistics examined, despite the Knudsen number based on the viscous length scale being as large as O(10 -1 ). This latter result is at variance with classical estimates for the breakdown of the continuum assumption and calls for further investigations into the interaction of noncontinuum effects and turbulence.

42 ENGINEERING↗

Development and application of hybrid AIMD/cDFT simulations for atomic-to-mesoscale chemistry

Many important chemical processes involve reactivity and dynamics in complex solutions. Gaining a fundamental understanding of these reaction mechanisms is a challenging goal that requires advanced computational and experimental approaches. However, important techniques such as molecular simulation have limitations in terms of scales of time, length, and system complexity. Furthermore, among the currently available solvation models, there are very few designed to describe the interaction between the molecular scale and the mesoscale. To help address this challenge, here, we establish a novel hybrid approach that couples first-principles plane-wave density functional theory with classical density functional theory (cDFT). In this approach, a region of interest described by ab initio molecular dynamics (AIMD) interacts with the surrounding medium described using cDFT to arrive at a self-consistent ground state. cDFT is a robust but efficient mesoscopic approach to accurate thermodynamics of bulk electrolyte solutions over a wide concentration range (up to 2M concentrations). Benchmarking against commonly used continuum models of solvation, such as SMD, as well as experiments, demonstrates that our hybrid AIMD–cDFT method is able to produce reasonable solvation energies for a variety of molecules and ions. With this model, we also examined the solvent effects on a prototype S N 2 reaction of the nucleophilic attack of a chloride ion on methyl chloride in the solution. The resulting reaction pathway profile and the solution phase barrier agree well with experiment, showing that our AIMD/cDFT hybrid approach can provide insight into the specific role of the solvent on the reaction coordinate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Rapid data-driven model reduction of nonlinear dynamical systems including chemical reaction networks using ℓ 1 -regularization

We develop a new data-driven paradigm for efficient model reduction of a broad class of nonlinear dynamical systems. Our model reduction method directly enables the interpretation of key components of the dynamical system, unlike traditional projection-based model reduction methods that focus on reducing computational complexity more than interpretability. Our method is not application specific and is simple to implement on nonlinear dynamical systems arising from a variety of different fields. It requires minimal parameterization using a single parameter to trade-off between model complexity and estimation error. We use a data-driven paradigm to formulate model reduction as an efficient convex optimization problem that scales polynomially in the original size of the complex system, enabling systems with as many as thousands of components to be reduced in a matter of minutes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Neutrons reveal the dynamics of leaf thylakoids in living plants

The study is the first known exploration of photosynthetic membranes dynamics in living plants by high resolution quasielastic neutron scattering spectroscopy. We investigated the mobility and flexibility of thylakoid membranes in common duckweed (Landoltia punctata) and identified dynamics across various length scales corresponding to individual membranes and membranes stack. We employed classical models typically used to study lipid bilayers to characterize the undulation modes and rigidity of the membranes and reveal how structural variations influence the observed complex dynamics. Our findings show that the stacks of thylakoids in duckweed behave as rigid systems, exhibiting an effective bending coefficient in the lower range associated with surfactant membranes. In contrast, the single thylakoid leaflets display greater apparent flexibility and are well situated within the bi-continuous surfactant phase dynamics. While our observations enhance the understanding of the intricate architecture and mobility of photosynthetic cellular machinery, they also highlight the limitations of applying ideal lipid membranes models to describe complex biological systems. This work opens more questions and the need for further investigations across extended length and time scales, as well as the importance of rigorous sample preparation and experimental control.

Applied physics↗

Development Efforts and Status of Plume-Surface Interaction Capabilities for Propulsive Landing Systems

Rocket plume-surface interaction is a multi-phase problem characterized by plume flow physics, erosion physics, and ejecta dynamics. All propulsive landers will experience plume-surface interaction. The risks posed by such effects can vary greatly as a function of the local environment, lander concept of operations, configuration, and physical scale. Prior Lunar and Martian landers have overcome challenges posed by these environments, on the basis of scaled ground testing and subsequent flight experience, but the landing systems for present and future missions are planning to operate increasingly outside of NASA's current experience with plume-surface interaction. This presentation discusses the current efforts and status of activities within NASA to progress understanding of fundamental plume-surface interaction physics, the capability to predict resulting environments and effects, and the definition of implications for current and future Lunar and planetary landing systems.

Landing↗

Atmosphere-Biosphere Interaction Study (Field Campaign Report)

The objective of this campaign was to collect hyperspectral reflectance measurements and complementary optical and infrared data of land surface (i.e., vegetation and soil) throughout the growing season (i.e., from March to September) at high temporal frequency (e.g., < 30 minutes). To collect the measurements, the EcoSpec system, a tower-based sensor system housing seven optical and infrared sensors that was developed at Argonne National Laboratory, was deployed to the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility Southern Great Plains (SGP) atmospheric observatory (within the crop, south of the eddy correlation flux measurement system [ECOR] tower). The EcoSpec system consists of a spectroradiometer having 2,151 channels (350–2,500 nm), red-greenblue (RGB) camera, diffuse radiometer, thermal infrared (TIR) sensors, albedometers, Photochemical Reflectance Index (PRI) sensors, Normalized Difference Vegetation Index (NDVI) sensors, and wetness sensor. The system is designed to collect an array of land surface measurements, including hyperspectral reflectance of plants and soils, RGB photos capturing contextual information of land surface, radiant temperature of sky and land surface, incoming and reflected shortwave radiation, and incoming radiation components (direct and diffuse) as frequent as every minute from up to 12 discrete positions around the tower. This would allow us to capture inherent heterogeneity of the terrestrial ecosystem and to obtain representative spectral signatures of the plant and soil surfaces via averaging those measurements. The SGP site provides advanced atmospheric and other physical measurements that are essential for our investigation of ecosystem functions using hyperspectral reflectance measurements of land surfaces. Located in the mid-latitudes, the SGP site exhibits dynamic seasonality and a wide range of meteorological variability and flux properties that make it ideal for studying interactions between near-surface atmosphere and biosphere across a range of environmental conditions. The atmosphere, plants, and soil control terrestrial carbon and water cycles. To make more accurate climate forecasts, researchers need to understand ecosystem dynamics at the biosphere–atmosphere interface. Even though our ability to forecast ecosystem dynamics and climate at both the regional and global scales has advanced significantly, we still do not know how local phenomena affecting water and carbon fluxes—such as diurnal variations—interrelate with large-scale atmospheric/climate phenomena, and vice versa. The lack of a more thorough understanding of ecosystem dynamics and climate interactions is a key knowledge gap for the improvement of earth system models (ESMs).

54 ENVIRONMENTAL SCIENCES↗

Multibody Based Digital Astronaut Dynamics Simulation

BACKGROUND: This study provides the Software, Robotics, & Simulation Division at the NASA Johnson Space Center with a verification tool for multibody dynamics simulation requiring human motion. The motivation stems from current studies of several Vibration Isolation & Stabilization(VIS)system designs that attenuate the moments and forces which would be transmitted to a spacecraft during an exercise. A multibody dynamics model for a proposed VIS was available previously[1], therefore modeling of the VIS was not needed for this work. The interest here is in creating the multibody dynamics model of an astronaut in motion which may be utilized independently or while attached to a mechanism. An existing simulation [2] that utilizes OpenSim [3,4] and an in-house multibody dynamics package (MBDyn) [5] is used in order to verify the astronaut model. The main advantage this model will have over the existing simulation is that everything will be processed in one tool. METHODS AND RESULTS: Creating the simulation required; estimation of Body Segment Inertial Parameters (BSIP),a multibody model of the human-VIS system, joint acceleration profiles, and input files for MBDyn, which is used for this analysis. The scaling factors provided by Dumas et al. [6] are utilized in estimating the BSIP. Anthropometric data are used for estimating these parameters, the Anthropometric Survey of US Army Personnel (ANSUR II) [7] was the source. The astronaut model consists of 15 bodies, 14 joints and 32 degrees of freedom, with the dynamics topology generated using the center of mass locations and anthropometric data. MBDyn has an option for prescribed joint motion (PJM), which requires joint acceleration data as input. The joint angle data is first obtained from a motion capture system and then processed through code that has been created to generate approximate joint acceleration profiles. The topology tree for the astronaut model begins at the right foot up to the pelvis where there is one branch for going down the left leg and another for the torso. The torso branch leads to branches for the arms and a leaf body for the head/neck segment. For attachment to the VIS, the heel of the right foot is connected to the VIS platform through a fixed joint, resembling a foot restraint. The left foot does not attach to the platform in order to prevent a system with a closed loop. Topology and symmetry of the astronaut model were verified through kinematic analysis. Further verification of the forces and moments transmitted to the VIS were verified against the existing simulation. There was a satisfactory level of agreement when testing a simple motion, for example, rocking back and forth. Full exercise motions are to be tested soon. The main outcome has been a novel application of MBDyn for biomechanics modeling that is now available for dynamic simulations involving human motion. The estimation of BSIP was another useful result of this study, requiring only 15 inputs for generating mass properties of a theoretical astronaut model. Expansion on this work is possible by going through an alternative route in obtaining the joint motion data. Instead of high-tech and often expensive motion capture systems, an individual may watch videos with high focus and at a slow motion for each individual segment in order to determine the initial and final time and angle for that specific degree of freedom. Synthetic trajectories may also be created if there is no video reference available.

F N Matari↗

A combined ensemble-volume average homogenization method for lattice structures with defects under dynamic and static loading

In the study of lattices structures, both experiments and numerical simulations are often conducted with small samples. Using combined ensemble and volume averaging, this work introduces a method to extract a macroscopic constitutive response of a lattice material from numerical simulations performed in periodic domains. The domain size needed to obtain statistically accurate results is investigated. Similar to molecular dynamics, the concept of the virial stress is introduced after homogenized equations are derived using the ensemble averaging method. Under static conditions, the virial stress is shown to agree with the volume averaged solid stress. Using the homogenization method, constitutive relations for this stress can be obtained from systems with uniform strains. Application of such obtained constitutive relations to more general cases results in an error proportional to the square of the ratio between the lattice length scale and the macroscopic length scale. Taking advantage of this property, numerical simulations are performed in systems with a uniform gradient of the average velocity. The volume average method is then used to accelerate convergence when studying lattices with defects. To avoid the artificial numerical time scale from the size of a representative volume element divided by the wave speed, a numerical scheme is developed to enforce a spatially uniform velocity gradient within the computational domain while allowing fluctuations of the velocity or displacement to develop naturally. To account for probability distribution of lattice defects, the stress is calculated as the ensemble-volume averaged value. For dynamic systems, energy dissipation properties are also studied.

36 MATERIALS SCIENCE↗

Adaptive variational simulation for open quantum systems

Emerging quantum hardware provides new possibilities for quantum simulation. While much of the research has focused on simulating closed quantum systems, the real-world quantum systems are mostly open. Therefore, it is essential to develop quantum algorithms that can effectively simulate open quantum systems. Here we present an adaptive variational quantum algorithm for simulating open quantum system dynamics described by the Lindblad equation. The algorithm is designed to build resource-efficient ansatze through the dynamical addition of operators by maintaining the simulation accuracy. We validate the effectiveness of our algorithm on both noiseless simulators and IBM quantum processors and observe good quantitative and qualitative agreement with the exact solution. We also investigate the scaling of the required resources with system size and accuracy and find polynomial behavior. Our results demonstrate that near-future quantum processors are capable of simulating open quantum systems.

97 MATHEMATICS AND COMPUTING↗

Improved Understanding of the Modeled QBO Using MLS Observations and MERRA Reanalysis

The Quasi-Biennial Oscillation (QBO) dominates the variability of the tropical stratosphere on interannual time scales. The QBO has been shown to extend its influence into the chemical composition of this region through dynamical mechanisms. We have started our analysis using the realistic QBO internally generated by the Goddard Earth Observing System Version 5 (GEOS-5) general circulation model coupled to a comprehensive stratospheric and tropospheric chemical mechanism forced with observed sea surface temperatures over the past 33 years. We will show targeted comparisons with observations from NASAs Aura satellite Microwave Limb Sounder (MLS) and the Modern Era Retrospective-Analysis for Research and Applications (MERRA) reanalysis to provide insight into the simulation of the primary and secondary circulations associated with the QBO. Using frequency spectrum analysis and multiple linear regression we can illuminate the resulting circulations and deduce the strengths and weaknesses in their modeled representation. Inclusion of the QBO in our simulation improves the representation of the subtropical barriers and overall tropical variability. The QBO impact on tropical upwelling is important to quantify when calculating trends in sub-decadal scale datasets.

MERRA↗

Conversion-Integration of MSFC Nonlinear Signal Diagnostic Analysis Algorithms for Realtime Execution of MSFC's MPP Prototype System

NASA's advanced propulsion system Small Scale Magnetic Disturbances/Advanced Technology Development (SSME/ATD) has been undergoing extensive flight certification and developmental testing, which involves large numbers of health monitoring measurements. To enhance engine safety and reliability, detailed analysis and evaluation of the measurement signals are mandatory to assess its dynamic characteristics and operational condition. Efficient and reliable signal detection techniques will reduce the risk of catastrophic system failures and expedite the evaluation of both flight and ground test data, and thereby reduce launch turn-around time. During the development of SSME, ASRI participated in the research and development of several advanced non- linear signal diagnostic methods for health monitoring and failure prediction in turbomachinery components. However, due to the intensive computational requirement associated with such advanced analysis tasks, current SSME dynamic data analysis and diagnostic evaluation is performed off-line following flight or ground test with a typical diagnostic turnaround time of one to two days. The objective of MSFC's MPP Prototype System is to eliminate such 'diagnostic lag time' by achieving signal processing and analysis in real-time. Such an on-line diagnostic system can provide sufficient lead time to initiate corrective action and also to enable efficient scheduling of inspection, maintenance and repair activities. The major objective of this project was to convert and implement a number of advanced nonlinear diagnostic DSP algorithms in a format consistent with that required for integration into the Vanderbilt Multigraph Architecture (MGA) Model Based Programming environment. This effort will allow the real-time execution of these algorithms using the MSFC MPP Prototype System. ASRI has completed the software conversion and integration of a sequence of nonlinear signal analysis techniques specified in the SOW for real-time execution on MSFC's MPP Prototype. This report documents and summarizes the results of the contract tasks; provides the complete computer source code; including all FORTRAN/C Utilities; and all other utilities/supporting software libraries that are required for operation.

Jong, Jen-Yi↗

Sliding Mode Transient Scaling Controller for Gas Turbine Engine Emulation on an Electric Machine

The testing of hybrid-electric aircraft propulsion system components and control strategies has dictated the need for a low-cost gas turbine engine (GTE) emulation testbed. Within that testbed, electric machines (EMs) are used to emulate the dynamics of a GTE. GTE mathematical models are used to provide control torque commands to EMs. Most GTEs operate at much higher power levels than commercially available EMs. Therefore, the transient shaft dynamics of the GTE must be scaled appropriately. The sliding mode controller being presented here aims to deliver a control torque to an electric machine with a fixed inertia that will allow it to mimic the response of a GTE with a different inertia. This solution provides an easily scalable and tunable transient emulation method without having to scale the GTE mathematical model and/or add inertia to the physical EM shaft. This controller has been proven in simulation to produce nearly identical transient rotational speed response to a reference GTE shaft model which has a large difference in rotational inertia compared to the motor. Also, the inherent robustness of the sliding mode control scheme results in negligible changes in performance for large amounts of parameter uncertainty while minimizing high frequency fluctuations in the control signal.

sliding mode↗

The Aerosol-Monsoon Climate System of Asia

In Asian monsoon countries such as China and India, human health and safety problems caused by air-pollution are worsening due to the increased loading of atmospheric pollutants stemming from rising energy demand associated with the rapid pace of industrialization and modernization. Meanwhile, uneven distribution of monsoon rain associated with flash flood or prolonged drought, has caused major loss of human lives, and damages in crop and properties with devastating societal impacts on Asian countries. Historically, air-pollution and monsoon research are treated as separate problems. However a growing number of recent studies have suggested that the two problems may be intrinsically intertwined and need to be studied jointly. Because of complexity of the dynamics of the monsoon systems, aerosol impacts on monsoons and vice versa must be studied and understood in the context of aerosol forcing in relationship to changes in fundamental driving forces of the monsoon climate system (e.g. sea surface temperature, land-sea contrast etc.) on time scales from intraseasonal variability (~weeks) to climate change (~ multi-decades). Indeed, because of the large contributions of aerosols to the global and regional energy balance of the atmosphere and earth surface, and possible effects of the microphysics of clouds and precipitation, a better understanding of the response to climate change in Asian monsoon regions requires that aerosols be considered as an integral component of a fully coupled aerosol-monsoon system on all time scales. In this paper, using observations and results from climate modeling, we will discuss the coherent variability of the coupled aerosol-monsoon climate system in South Asia and East Asia, including aerosol distribution and types, with respect to rainfall, moisture, winds, land-sea thermal contrast, heat sources and sink distributions in the atmosphere in seasonal, interannual to climate change time scales. We will show examples of how elevated absorbing aerosols (dust and black carbon) may interact with monsoon dynamics to produce feedback effects on the atmospheric water cycle, leading to in accelerated melting of snowpacks over the Himalayas and Tibetan Plateau, and subsequent changes in evolution of the pre-monsoon and peak monsoon rainfall, moisture and wind distributions in South Asia and East Asia.

Lau, William K. M.↗

Scaling dynamics in low-salt-rejection reverse osmosis for high-salinity produced water desalination: Mechanistic modeling and membrane autopsy

Membrane scaling remains a critical barrier to the reliable operation of desalination systems, particularly for hypersaline produced water (PW) treatment. This study fills the knowledge gap of autopsy-based model validation for PW desalination by elucidating scaling mechanisms in a Low-Salt-Rejection Reverse Osmosis (LSRRO) system through the integration of pilot-scale experimentation and complementary modeling approaches. A semi-empirical modeling framework was developed and applied to a multistage pilot LSRRO system equipped with nanofiltration and RO membranes treating high-salinity PW from the Permian Basin. Water quality analysis showed that total dissolved solids decreased from ~130,000 mg/L to ~1900 mg/L in the permeate, then further reduced to ~300 mg/L by a second-pass RO. Two different thermodynamic modeling approaches were evaluated: the first extends the LSRRO framework by incorporating system complexity and scaling phenomena, whereas the second method explicitly captures concentration polarization in localized supersaturation. Both methods illustrate the tendency for carbonate and sulfate scaling throughout the stages. Membrane autopsies revealed a silica-dominated deposit matrix, localized CaSO 4 at Stage 2, and minor barite/celestite despite their prominence in model predictions. Quantum-chemical calculations indicated silica scaling can be rationalized by favorable adsorption of H 4 SiO 4 on Fe-oxide surfaces (ΔG ≈ −44 kJ/mol), providing a kinetic pathway for interfacial inorganic polymerization even when bulk equilibrium predictions are conservative. Overall, the thermodynamic scaling modeling and membrane autopsy revealed heterogeneous, localized deposits with limited impact on LSRRO performance, while quantum analysis rationalized the thermodynamically unfavorable precipitation formation under bulk equilibrium, reconciling model–autopsy discrepancies. These insights support targeted pretreatment and silica-specific antiscalants to extend membrane lifetime and increase recovery, providing a transferable framework for hypersaline water desalination systems. The combined experimental–computational approach provides new mechanistic insight into scaling in hypersaline membrane systems and establishes a transferable framework for predicting and mitigating scaling in next-generation desalination technologies.

Low-salt-rejection reverse osmosis↗