Computational Reverse-Engineering Analysis for Scattering Experiments (CREASE) with Machine Learning
Explore the source record for details and available documents.
Heil, Christian M.↗
SEARCH · Engineering Papers
Results for “reversible computing”
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.
Logical and Physical Reversibility of Conservative Skyrmion Logic
Magnetic skyrmions are nanoscale whirls of magnetism that can be propagated with electrical currents. The repulsion between skyrmions inspires their use for reversible computing based on the elastic billiard ball collisions proposed for conservative logic in 1982. In this letter, we evaluate the logical and physical reversibility of this skyrmion logic paradigm, as well as the limitations that must be addressed before dissipation-free computation can be realized.
42 ENGINEERING↗
Electrochemical Oxidation and Speciation of Lanthanides in Potassium Carbonate Solution
Increasing lanthanide demand to support clean energy goals drives the need to develop more efficient approaches to separate adjacent lanthanides. Most approaches for lanthanide separations are not very selective and are based on small differences in lanthanide ionic radii. Concentrated potassium carbonate media has shown some potential to enable oxidation of praseodymium (Pr) and terbium (Tb) to their tetravalent states, which could ultimately enable a separation based on differences in oxidation states, but very little is known regarding the system's chemistry. This work completes a detailed examination of cerium (Ce) redox chemistry in concentrated carbonate media to support the development of Pr and Tb oxidation studies. The half-wave potential (E 1/2 ) of the Ce(III)/(IV) redox couple is evaluated under various solution conditions and computational modeling of carbonate coordination environments is discussed. Cyclic voltammetry shows higher carbonate concentrations and temperatures can lower the potential required to oxidize Ce(III) by 54 mV (3.5 to 5.5 M) and 39 mV (from 10 °C to 70 °C). Chronoabsorptometry shows Ce(III) and Ce(IV) carbonate complexes are chemically stable and reversible. Computational modelling suggests the most likely coordination environment for the Ce(IV) complex is Ce(CO 3 ) 4 (OH) 5– which is less entropically favorable than the lowest energy Ce(III) complex, Ce(CO 3 ) 4 5– .
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Quantifying dispersity in size and shape of nanoparticles from small-angle scattering data using machine learning based CREASE
Here, we use machine learning (ML) enhanced computational reverse engineering analysis of scattering experiments (CREASE) to interpret small-angle X-ray scattering (SAXS) data obtained from a system of nanoparticles without a priori knowledge of their exact shapes (e.g. spheres or ellipsoids), sizes (0.5–50 nm) and distributions. The SAXS measurements yielded three categories of scattering profiles exhibiting 'strong', 'weak' and 'no' features. Diminishing features (e.g. broadening or disappearing peaks) in scattering profiles have always been attributed to the presence of significant dispersity in the system. Such featureless SAXS data are not suitable for traditional analysis using analytical models. If one were to fit a relevant analytical model (e.g. the lmfit analytical model for polydisperse spheres) to these 'weak' and 'no' SAXS profiles from our nanoparticle systems, one would obtain non-unique interpretations of the data. Relying on electron microscopy to identify the distributions of nanoparticle shapes and sizes is also unfeasible, especially in high-throughput synthesis and characterization loops. In such situations, to identify the distributions of particle sizes and shapes that could be present in the sample, one must rely on methods like ML-CREASE to interpret the data quickly and output all relevant interpretations about the structure present in the system. The ML-CREASE optimization loop takes the experimental scattering profile as input and outputs multiple candidate solutions whose computed scattering profiles match the SAXS profile input. The ML-CREASE method outputs distributions of relevant structural features, such as the volume fraction of the nanoparticles in the system and the mean and standard deviation of the particle size and aspect ratio, assuming a type of distribution (e.g. normal, log-normal) for size and aspect ratio. We find that, for the SAXS profiles analyzed here, accounting for the shape dispersity along with size dispersity of the nanoparticles using ML-CREASE improved the match between the computed scattering profiles and input experimental profiles.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Tutorial: Machine-Learning-Based CREASE-2D Analysis of 2D SAXS Profiles to Characterize Anisotropic Nanostructures in Soft Materials
We present a tutorial to guide users on how to extend the Computational Reverse Engineering Analysis of Scattering Experiments-2D (CREASE-2D) framework to interpret their experimental two-dimensional small-angle scattering (SAS) data from soft materials (e.g., polymers, peptide amphiphiles, biomolecular fibrils). Unlike most traditional SAS analysis approaches, which typically rely on azimuthally averaged onedimensional (1D) profiles, CREASE-2D utilizes the complete 2D scattering profile to reveal information about anisotropy in the structure. In past applications, CREASE has provided insights into complex structural features, including the cross-sectional shapes of assembled nanostructures and dispersity in these features, which are difficult to discern with existing analytical models. While (1D- ) CREASE has been applied to SANS and SAXS data, this tutorial shares the steps for implementing CREASE-2D using an example of a dipeptide solution system, for which we have SAXS data. We present details for these steps involved in using CREASE-2D to interpret SAXS profiles: how to preprocess SAXS data, define relevant structural features, generate three-dimensional real-space structures for specific values of these features, train a machine learning (ML) surrogate model to predict scattering profiles for given structural features, and optimize these features using genetic algorithms (GA). Then, we use these steps to interpret complex 2DSAXS data collected from dipeptide solutions that, in microscopy images, exhibit nanoscale structures that could be elliptical tubes/ flat tapes/cylinders or a combination of these cross sections. Open-source codes, computational hardware, and software requirements, as well as the strengths and limitations of this protocol, are also presented. We expect researchers working with (soft) biomaterials, peptide amphiphiles, amphiphilic polymer solutions, polymer nanocomposites, and blends of particles/polymers will find this CREASE-2D method and this tutorial of use.
CREASE↗
BARC Element Classifier (barc)
SAND2024-02120O This program enumerates and classifies all possible functional elements in several different categories within the Ballistic Asynchronous Reversible Computing (barc) model of computation. barc software is a research tool that is being used to help document the possible digital behaviors of primitive functional elements in the ABRC a.k.a. BARC model of computation. The method of operation of this software leverages elementary concepts and methods from discrete mathematics: combinatorics (permutations), partitions and equivalence classes, symmetry transformations, and the theory of finite groups. Implementation is in Python programming language. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.
Frank, Michael↗
Pressure-driven tearing and thermal transport in finite-beta reversed field pinch computations
In this work, nonlinear resistive-magnetohydrodynamics (MHD) computation with heating and anisotropic transport is applied to examine the interaction between thermal energy and magnetic fluctuations in inductively driven reversed-field pinches (RFPs). The magnetic fluctuations underlie magnetic field reversal through dynamo-like correlations, and they enhance thermal energy transport through fluctuations of parallel heat flux density. With the unfavorable magnetic curvature that exists across the RFP profile, thermal energy also affects the magnetic fluctuations. Computations with the NIMROD code [Sovinec et al., J. Comput. Phys. 195, 355–386 (2004)] integrate nonlinear MHD dynamics with energy transport and reproduce an RFP state with experimentally relevant values of plasma-β. Equilibria constructed from results of the 3D computations are analyzed to assess the sources of free energy in the saturated nonlinear state. Linear computations for these profiles show unstable modes of tearing parity. Their eigenfunctions are used to evaluate and compare stabilizing and destabilizing contributions to the kinetic energy integral. An assessment of the drives in the integral reveals that the pressure gradient drive is of comparable magnitude to the parallel current drive, and only the sum of the two surpasses the stabilizing contributions. Correlation of magnetic and parallel heat flux density fluctuations in the nonlinear computations shows that fluctuation-induced thermal conduction is the dominant mode of energy loss, as expected from experimental evidence. Decomposition of the fluctuating heat flux density shows that second-order correlations, alone, do not explain the total energy transport. Higher-order correlations are also important.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Computer Science Research Needs for Parallel Discrete Event Simulation (PDES)
Historically, scientific computing efforts have demonstrated the clear need for, and effective use of, supercomputing with traditional time-stepped simulations. Nevertheless, there are several areas in the mission spaces of the U.S. Department of Energy and other agencies waiting to tap advanced computing research using a different, discrete event style of modeling, simulation, and analysis. These span a wide spectrum of applications including energy grid resilience, urban planning and policy, transportation science, building technologies, emergency response and planning, environmental impact analysis, computational epidemiology, Internet communications, cyber security, and cyber-physical systems, to name only a few. Even within traditional scientific applications, the role of discrete event modes of execution is increasing in the form of new event-based mathematical solvers such as quantized state integration methods and discrete-continuous hybrid system solvers. Co-design of advanced supercomputing hardware systems is another area that exploits discrete event simulation at its core for effective analyses. Complex systems, entity behaviors and interconnections play a significant role in all these applications, which are mapped to large-scale models with discrete event formulations. To make advancements in all the aforementioned scientific areas, many technical aspects need to be more thoroughly studied and deeply understood in parallel discrete event simulation (PDES). The unique dynamics inherent in a discrete event modeling approach, by their very nature, intersect and influence the entire stack of the computing system, including (a) the unique nature of the instruction sets exercised in PDES workloads without a predominance of high-precision floating point operations, (b) virtual time-constrained multi-threaded execution of many logical processes per processor, (c) extremely variable and difficult to predict network traffic characteristics, (d) interfaces and inter-dependencies with machine learning and artificial intelligence codes at higher software layers, and (e) highly challenging load balancing needs, especially in effectively accounting for accelerated/extremely heterogeneous computing in current and future high-performance computing systems. Efficient and accurate parallel execution of PDES workloads is also dominated by challenges in dealing with their asynchronous concurrency fundamentally present at the model level. Conservative synchronization, optimistic/speculative synchronization, and their hybrid schemes open new questions in fundamental computer science with respect to reversibility of computation and prediction (lookahead) of behaviors inherent within model codes. On the implementation front, there are relatively few scalable, general-purpose parallel discrete event simulators in the world, and even fewer have been studied on emerging hardware platforms. To enable scientific advances using PDES, the research needs in computer science must also be pursued and met in the intersection of the algorithmic and hardware-aware aspects of scalable PDES engines. This report is aimed at capturing a computer science-oriented view of this important area of research in PDES, presenting a sample of important applications with their inherent discrete event technology elements. Needs are outlined in core areas of parallel discrete event research as well as cross-cutting directions in computer science research that positively impact scientific advancements across several important application areas. A selection of priority research opportunities in advanced computing for PDES is identified to serve as reference for key research topics and their order of importance for scientific advancements.
97 MATHEMATICS AND COMPUTING↗
Designing Energy-Efficient Quantum Computers Through Prediction and Reduction of Cooling Requirements for Cryogenic Electronics
Quantum computing has been identified as a “wild card” by the International Energy Agency in predicting future global data center energy usage. This is primarily because both uncertainty in the extent to which quantum computing will be adopted, and uncertainty in the power consumption of individual quantum data centers. Unlike the classical counterparts, quantum computers need to be maintained at near absolute zero, requiring energy-intensive cryogenic cooling systems. Therefore, as quantum computers scale up from existing 50 qubit technology demonstrations to the 10,000 to 100,000 qubit systems that will be able to solve complex problems, the energy consumption of both the electronics and the required cooling systems will also increase. To predict this scaling, this work analyzes the energy requirements for both computation and cooling of quantum hardware. We show that the energy requirements for cooling of quantum computers is determined by several computing system parameters, including the number and type of physical qubits, the operating temperature, the packaging efficiency of the system, and the split between circuits operating at cryogenic temperatures and those operating at room temperature. The energy requirements can then be found based on thermal system parameters such as cooling efficiency and cryostat heat transfer. Analysis of these parameters shows that the energy required for cooling is significantly larger than that required for computation, a reversal from energy usage patterns seen in conventional computing. The results and discussions provide a road-map for creating energy efficient quantum computers through the selection of computer architectures and cryogenic system configurations that minimize cooling requirements.
energy efficiency↗
Mechanism of structural colors in binary mixtures of nanoparticle-based supraballs
Inspired by structural colors in avian species, various synthetic strategies have been developed to produce noniridescent, saturated colors using nanoparticle assemblies. Nanoparticle mixtures varying in particle chemistry and size have additional emergent properties that affect the color produced. For complex multicomponent systems, understanding the assembled structure and a robust optical modeling tool can empower scientists to identify structure-color relationships and fabricate designer materials with tailored color. Here, we demonstrate how we can reconstruct the assembled structure from small-angle scattering measurements using the computational reverse-engineering analysis for scattering experiments method and use the reconstructed structure in finite-difference time-domain calculations to predict color. We successfully, quantitatively predict experimentally observed color in mixtures containing strongly absorbing nanoparticles and demonstrate the influence of a single layer of segregated nanoparticles on color produced. The versatile computational approach that we present is useful for engineering synthetic materials with desired colors without laborious trial-and-error experiments.
36 MATERIALS SCIENCE↗
Two Circuits for Directing and Controlling Ballistic Fluxons
Reversible logic schemes using flux solitons (fluxons) on long Josephson junctions (LJJs) have recently been proposed. The attraction of the fluxon is that it propagates ballistically along an LJJ until it encounters a change in the character of the LJJ, often a designed circuit element. Logic gates involve fluxons interacting with circuit elements and with other fluxons. However, testing of ballistic fluxon circuits requires other circuits outside the logic family to direct and control fluxon motion. Here we discuss two such non-reversible fluxon control circuits. First, the polarity filter gate is a simple non-reversible gate that allows one polarity of fluxon to pass, while reflecting the other polarity. In the off state both polarities reflect. Second, the polarity separator generalizes on the polarity filter concept and allows separation of the two fluxon polarities into different LJJs. We discuss simulations of these structures and possible applications.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Qubit Assignment Using Time Reversal
As quantum computers with large numbers of qubits become increasingly available, experiments executed on a given device may not utilize all available qubits. In this case, the outcome of executing a quantum program will depend on the ability to efficiently select a subset of high-performing physical qubits. For any given quantum program and device there are many ways to assign physical qubits for execution of the program, and assignments will differ in performance due to the variability in quality across qubits and entangling operations on a single device. Evaluating the performance of each assignment using fidelity estimation introduces significant experimental overhead and will be infeasible for many applications, while relying on standard device benchmarks provides incomplete information about the performance of any specific program. Furthermore, the number of possible assignments grows combinatorially in the number of qubits on the device and in the program, motivating the use of heuristic optimization techniques. We demonstrate a practical solution to the problem of qubit assignment by using simulated annealing with a cost function based on the Loschmidt echo, a diagnostic that measures the reversibility of a quantum process. We provide theoretical justification for this choice of cost function by demonstrating that the optimal qubit assignment coincides with the optimal qubit assignment based on state fidelity in the weak error limit, and we provide experimental justification using diagnostics performed on Google’s superconducting qubit devices. We then establish the performance of simulated annealing for qubit assignment using classical simulations of noisy devices as well as optimization experiments performed on a quantum processor. Our results demonstrate that the use of Loschmidt echoes and simulated annealing provides a scalable and flexible approach to optimizing qubit assignment on near-term hardware.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Generative Thermodynamic Computing
Here, we introduce a generative modeling framework for thermodynamic computing, in which structured data are synthesized from noise by the natural time evolution of a physical system governed by Langevin dynamics. While conventional diffusion models use neural networks to perform denoising, here the information needed to generate structure from noise is encoded by the dynamics of a thermodynamic system. Training proceeds by maximizing the probability with which the computer generates the reverse of a noising trajectory, which ensures that the computer generates data with minimal heat emission. We demonstrate this framework within a digital simulation of a thermodynamic computer. If realized in analog hardware, such a system would function as a generative model that produces structured samples without the need for artificially injected noise or active control of denoising.
Whitelam, Stephen [Lawrence Berkeley National Labo↗
Reversing inverse Hall-Petch and direct computation of Hall-Petch coefficients
A 2D bicrystal atomistic model of dislocation transmission through a Σ11<101>{131} symmetric-tilt grain boundary reveals details of Hall-Petch breakdown in the single dislocation regime in Al, Ni, and Cu. Partly based on a previous study, this research determines the stress required for a single dislocation to be transmitted through the Σ11 boundary and finds that single dislocation transmission typically deviates from Hall-Petch behavior because the leading partial dislocation becomes trapped in the boundary and emits glissile grain boundary disconnections that cause boundary sliding deformation at stresses well below boundary transmission stresses. Thus, mechanisms of inverse Hall-Petch, namely grain boundary shear and sliding, are shown to operate in lieu of grain boundary transmission in the single dislocation regime. However, by controlling the applied stresses, inverse Hall-Petch can be reversed and typical Hall- Petch grain boundary transmission regained. This, in turn, allows the direct computation of Hall-Petch coefficients. Further, a second dislocation on the same slip plane preserves typical Hall-Petch behavior for Al and Ni and does not lead to boundary sliding events, thus implying a critical grain size for Hall-Petch breakdown based on the pileup model.
36 MATERIALS SCIENCE↗
Fostering a Guiding Multiscale Model for the Development of Advanced MgB 2 Hydrogen Storage Materials (Final Technical Report)
Project Goal and Objective. The demand for energy and for an upgraded energy infrastructure has steadily grown, as have the needs for energy independence and alternatives to our reliance on petroleum. Hydrogen is considered the most viable fuels for wide-scale implementation in the near future as it is less-polluting, non-toxic, and has more stored energy than petroleum. It is envisioned that hydrogen can eventually become the prime energy carrier, integrating the transportation, grid, and chemical sectors in a way that improves resiliency, diversifies feedstocks, and affords new economic opportunities. A key remaining challenge is the development materials with enhanced gravimetric and volumetric hydrogen storage capacities that offer a higher performance than compressed gas. These materials would eliminate the need for large-scale compression, thereby dramatically reducing the footprint and cost of gas storage. The high gravimetric and volumetric hydrogen capacities of complex hydrides has prompted an intensive investigation of the potential of this class of materials as hydrogen storage media over the past 25 years. Among the many complex hydrides that have been explored, magnesium borohydride, Mg(BH 4 ) 2 , has been found to possess the best combination of practical thermodynamic properties. These include a gravimetric H 2 density of 14.9 wt% H 2 and thermodynamics for the dehydrogenation of Mg(BH 4 ) 2 to MgB 2 (equation 1) (ΔH° = 39 kJ/mol H 2 , ΔS = 112 J/K mol H 2 ) which lie in the narrow window required Mg(BH 4 ) 2 $\Leftrightarrow$ MgB 2 + 4 H 2 (1) for reversibility under moderate pressure and temperature. However, overcoming the extremely slow kinetics of the reversible release of hydrogen by this material in the solid state is a daunting challenge. At temperatures greater than 400 °C, the borohydride releases up to 14 wt% hydrogen giving MgB 2 . We discovered that the direct re-hydrogenation of MgB 2 to Mg(BH 4 ) 2 can be accomplished under 950 bar H 2 at 400 °C. While this demonstrated that complete reversibility can be achieved, the conditions employed are far too extreme for commercial hydrogen storage applications. More recently, we found through US DOE funded research projects (EERE HyMARC and HySCOR), that hydrogen cycling, can be accomplish at much milder conditions upon modification of the borohydride or boride. Guided by these discoveries these discoveries, the objective of this research project was to obtain key information that will enable the development of a model of reversible hydrogenation of MgB 2 to Mg(BH 4 ) 2 . The ultimate goal of our efforts is to attain a model of this transformation that can be utilized to accelerate development further advanced materials. This project directly follows on discoveries that were made over the course of a US DOE, EERE HyMARC project that was focused on improvement of the hydrogen cycling kinetics of modified MgB 2 . We found that that mechanical milling with graphene results the desired, pronounced kinetic enhancement. The dramatic lowering of the conditions required for the hydrogenation of MgB 2 is a significant step towards overcoming its chemical inertness allowing its development as a practical onboard hydrogen storage material. However, the exact nature of the modification(s) of MgB 2 that is responsible for its activation towards hydrogenation is completely unknown. This situation is not unique, as efforts to develop hydrogen storage materials typically have a narrow focus rather than a comprehensive approach that takes atomic level bonding and structure; molecular dynamics; long range, nano- and mesoscale-structure and their interconnection all into account. The goal of this project was the development of a comprehensive, multi-scale computational model of reversible hydrogenation of MgB 2 to Mg(BH 4 ) 2 that can be utilized for development of higher performance versions of the modified material. Development of the model requires determination of: 1) the bulk, nano-scale, and meso-scale structural changes occurring at elevated pressure following mechano-chemical modification of MgB 2 ; 2) the reaction pathway of the reversible hydrogenation of MgB 2 to Mg(BH 4 ) 2 ; 3) the effect of elevated pressure and mechano-chemical modification on the chemical reaction pathways; 4) the interactions at solid-gas interfaces; and particle surfaces; and 5) the kinetics and thermodynamic parameters associated with each step of the hydrogenation reaction pathway. This investigation required advanced techniques as preliminary, standard XRD, 11 B NMR, and FTIR analysis showed no signs of material modification. In order to gain this level of understanding of modified MgB 2 , required the teaming of a diverse group of experts and state-of-the art experimental capabilities at the University of Hawaii at Manoa (UHM) and collaborating National Laboratories: Craig Jensen , Department of Chemistry (PI and Project Director), solid state, solution, and high pressure NMR spectroscopy; solid-state synthesis; and high pressure hydrogenation (collaboration with SNL); Godwin Severa , Hawaii Natural Energy Institute (co-PI) calorimetry; infrared and Raman spectroscopy (collaboration with NREL); Dera , high pressure X-ray diffraction including in situ experiments (collaboration with ANL); Hope Ishii , Hawaii Institute of Geophysics electron microscopy investigations (collaboration with LBNL); and Joe Brown , Mechanical Engineering , material electronic structure and electric field effects.
08 HYDROGEN↗
Phase-space sensitivity (weight functions) of 3 MeV proton diagnostics
Abstract The weight function W ( X ) is a diagnostic sensitivity to phase-space variables X that relates the measured signal C to the distribution function F ( X ) through the equation C = ∫ W ( X ) F ( X ) d X . In the present work, an algorithm to calculate W for a diagnostic that measures 3 MeV protons produced in d(d,p)t fusion reactions between a fast ion and a thermal deuteron is developed. The emitted protons escape the tokamak on curved orbits and are detected. These curved orbits constitute effective diagnostic ‘sightlines’. The presented algorithm accounts for the complications associated with these curved sightlines. An initial calculation of time-reversed proton orbits computes effective solid angles and sightlines for the relevant range of incident proton velocity vectors. These precomputed orbits are then used within the framework of FIDASIM (2020 Plasma Phys. Cont. Fusion 62 105008) to calculate the reactivity averaged over the thermal distribution of the ‘target’ deuterons and the probability that a fast ion of specified energy and pitch has a gyroangle that is consistent with the kinematic equations along each of the sightlines. Comparisons with analytic formulas and with independent calculations for the Mega Amp Spherical Tokamak 3 MeV proton diagnostic verify the algorithm.
Physics↗
Shear reflectivity compensation in full-waveform inversion using least-squares reverse-time migration
SUMMARY The computational cost of elastic-waveform inversion is too high for inverting PP reflections, while using acoustic full-waveform inversion (FWI) is inaccurate because it does not depend on the shear modulus/velocity/impedance that affects elastic PP wavefield amplitudes. To solve this problem, we develop a waveform inversion method that uses acoustic least-squares reverse-time migration (LSRTM) to compensate the shear reflectivity for acoustic FWI. Our method is based on the quasi-elastic-wave equation developed by Chapman et al. (2014). The quasi-elastic-wave equation uses a linearized acoustic-wave equation with shear modulus μ as a virtual source to correct the acoustic PP wavefield amplitudes toward elastic ones. Our waveform inversion method inverts for elastic parameters by minimizing the L2 norm of the difference between recorded and predicted PP reflections modelled using the quasi-elastic-wave equation. Numerical tests on synthetic and field data show that our method can properly handle the amplitudes of elastic PP reflections and provides an accurate estimate of the P- and S-wave velocities/impedances and, in some cases, the density. The method does not need the computationally expensive numerical solution to the elastic-wave equation. It also gives a better estimate of elastic parameters than a pure LSRTM method for elastic PP reflections.
Feng, Zongcai↗
Impact of large-scale effects on mass transfer and concentration polarization in Reverse Osmosis membrane systems
Here, we present well-resolved computational fluid dynamics simulations of a large-scale reverse osmosis membrane-spacer configuration (1 m). Our computational model solves the flow and transport equations with variable solute-dependent properties. We utilize a high resolution computational mesh to resolve all relevant length scales associated with spacer-induced mixing and thin concentration boundary layers. An important contribution of this work is the development of a modified mass-transfer correlation that accounts for the development of the concentration boundary layer along the channel. A set of 2D axisymmetric simulations were performed for a spiral wound module layer with varying cross-flow conditions and spacer diameters which indicate a significant entrance length effect for concentration profile development at lower flow rates while mixing effects dominate at higher flow rates. The mass-transfer correlations at higher flow rates compare well with published correlations while a surrogate model for Sherwood number was obtained that depends on an additional similarity variable that accounted for entrance length effects at lower flow rates. Finally, a large-scale membrane-spacer design relevant to high-pressure reverse osmosis is studied with a non-uniform arrangement of spacers, which indicate a substantial saving in pressure drop (~40%) compared to traditional uniformly spaced pattern with minor variations (~2%) in concentration polarization, product water quality (~1%) and water recovery (~7%) compared to a uniform spacer pattern.
97 MATHEMATICS AND COMPUTING↗