Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Functional priors”

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 163 records · Page 9

Analysis and Optimization of Seismo-Acoustic Monitoring Networks with Bayesian Optimal Experimental Design

The Bayesian optimal experimental design (OED) problem seeks to identify data, sensor configurations, or experiments which can optimally reduce uncertainty. The goal of OED is to find an experiment that maximizes the expected information gain (EIG) about quantities of interest given prior knowledge about expected data. Therefore, within the context of seismic monitoring, we can use Bayesian OED to configure sensor networks by choosing sensor locations, types, and fidelity in order to improve our ability to identify and locate seismic sources. In this work, we develop the framework necessary to use Bayesian OED to optimize the ability to locate seismic events from arrival time data of detected seismic phases. In order to do utilize Bayesian OED we must develop four elements:1. A likelihood function that describes the uncertainty of detection and travel times; 2. A Bayesian solver that takes a prior and likelihood to identify the posterior; 3. An algorithm to compute EIG; and, 4. An optimizer that finds a sensor network which maximizes EIG. Once we have developed this framework, we can explore many relevant questions to monitoring such as: how and what multiphenomenology data can be used to optimally reduce uncertainty, how to trade off sensor fidelity and earth model uncertainty, and how sensor types, number, and locations influence uncertainty

47 OTHER INSTRUMENTATION↗

A mixed, unified forward/inverse framework for earthquake problems: fault implementation and coseismic slip estimate

SUMMARY We introduce a new finite-element (FE) based computational framework to solve forward and inverse elastic deformation problems for earthquake faulting via the adjoint method. Based on two advanced computational libraries, FEniCS and hIPPYlib for the forward and inverse problems, respectively, this framework is flexible, transparent and easily extensible. We represent a fault discontinuity through a mixed FE elasticity formulation, which approximates the stress with higher order accuracy and exposes the prescribed slip explicitly in the variational form without using conventional split node and decomposition discrete approaches. This also allows the first order optimality condition, that is the vanishing of the gradient, to be expressed in continuous form, which leads to consistent discretizations of all field variables, including the slip. We show comparisons with the standard, pure displacement formulation and a model containing an in-plane mode II crack, whose slip is prescribed via the split node technique. We demonstrate the potential of this new computational framework by performing a linear coseismic slip inversion through adjoint-based optimization methods, without requiring computation of elastic Green’s functions. Specifically, we consider a penalized least squares formulation, which in a Bayesian setting—under the assumption of Gaussian noise and prior—reflects the negative log of the posterior distribution. The comparison of the inversion results with a standard, linear inverse theory approach based on Okada’s solutions shows analogous results. Preliminary uncertainties are estimated via eigenvalue analysis of the Hessian of the penalized least squares objective function. Our implementation is fully open-source and Jupyter notebooks to reproduce our results are provided. The extension to a fully Bayesian framework for detailed uncertainty quantification and non-linear inversions, including for heterogeneous media earthquake problems, will be analysed in a forthcoming paper.

58 GEOSCIENCES↗

Analysis of the Trusted Inertial Terrain-Aided Navigation Measurement Function

The trusted inertial terrain-aided navigation (TITAN) algorithm leverages an airborne vertical synthetic aperture radar to measure the range to the closest ground points along several prescribed iso-Doppler contours. These TITAN minimum-range, prescribed-Doppler measurements are the result of a constrained nonlinear optimization problem whose optimization function and constraints both depend on the radar position and velocity. Owing to the complexity of this measurement definition, analysis of the TITAN algorithm is lacking in prior work. This publication offers such an analysis, making the following three contributions: (1) an analytical solution to the TITAN constrained optimization measurement problem, (2) a derivation of the TITAN measurement function Jacobian, and (3) a derivation of the Cramér-Rao lower bound on the estimated position and velocity error covariance. These three contributions are verified via Monte Carlo simulations over synthetic terrain, which further reveal two remarkable properties of the TITAN algorithm: (1) the along-track positioning errors tend to be smaller than the cross-track positioning errors, and (2) the cross-track positioning errors are independent of the terrain roughness.

TITAN↗

How to Obtain the Redshift Distribution from Probabilistic Redshift Estimates

Abstract A reliable estimate of the redshift distribution n ( z ) is crucial for using weak gravitational lensing and large-scale structures of galaxy catalogs to study cosmology. Spectroscopic redshifts for the dim and numerous galaxies of next-generation weak-lensing surveys are expected to be unavailable, making photometric redshift (photo- z ) probability density functions (PDFs) the next best alternative for comprehensively encapsulating the nontrivial systematics affecting photo- z point estimation. The established stacked estimator of n ( z ) avoids reducing photo- z PDFs to point estimates but yields a systematically biased estimate of n ( z ) that worsens with a decreasing signal-to-noise ratio, the very regime where photo- z PDFs are most necessary. We introduce Cosmological Hierarchical Inference with Probabilistic Photometric Redshifts ( CHIPPR ), a statistically rigorous probabilistic graphical model of redshift-dependent photometry that correctly propagates the redshift uncertainty information beyond the best-fit estimator of n ( z ) produced by traditional procedures and is provably the only self-consistent way to recover n ( z ) from photo- z PDFs. We present the chippr prototype code, noting that the mathematically justifiable approach incurs computational cost. The CHIPPR approach is applicable to any one-point statistic of any random variable, provided the prior probability density used to produce the posteriors is explicitly known; if the prior is implicit, as may be the case for popular photo- z techniques, then the resulting posterior PDFs cannot be used for scientific inference. We therefore recommend that the photo- z community focus on developing methodologies that enable the recovery of photo- z likelihoods with support over all redshifts, either directly or via a known prior probability density.

79 ASTRONOMY AND ASTROPHYSICS↗

NREM sleep as a novel protective cognitive reserve factor in the face of Alzheimer's disease pathology

Alzheimer’s disease (AD) pathology impairs cognitive function. Yet some individuals with high amounts of AD pathology suffer marked memory impairment, while others with the same degree of pathology burden show little impairment. Why is this? One proposed explanation is cognitive reserve i.e., factors that confer resilience against, or compensation for the effects of AD pathology. Deep NREM slow wave sleep (SWS) is recognized to enhance functions of learning and memory in healthy older adults. However, that the quality of NREM SWS (NREM slow wave activity, SWA) represents a novel cognitive reserve factor in older adults with AD pathology, thereby providing compensation against memory dysfunction otherwise caused by high AD pathology burden, remains unknown. Here, we tested this hypothesis in cognitively normal older adults (N = 62) by combining 11 C-PiB (Pittsburgh compound B) positron emission tomography (PET) scanning for the quantification of β-amyloid (Aβ) with sleep electroencephalography (EEG) recordings to quantify NREM SWA and a hippocampal-dependent face-name learning task. We demonstrated that NREM SWA significantly moderates the effect of Aβ status on memory function. Specifically, NREM SWA selectively supported superior memory function in individuals suffering high Aβ burden, i.e., those most in need of cognitive reserve (B = 2.694, p = 0.019). In contrast, those without significant Aβ pathological burden, and thus without the same need for cognitive reserve, did not similarly benefit from the presence of NREM SWA (B = -0.115, p = 0.876). This interaction between NREM SWA and Aβ status predicting memory function was significant after correcting for age, sex, Body Mass Index, gray matter atrophy, and previously identified cognitive reserve factors, such as education and physical activity (p = 0.042). These findings indicate that NREM SWA is a novel cognitive reserve factor providing resilience against the memory impairment otherwise caused by high AD pathology burden. Furthermore, this cognitive reserve function of NREM SWA remained significant when accounting both for covariates, and factors previously linked to resilience, suggesting that sleep might be an independent cognitive reserve resource. Beyond such mechanistic insights are potential therapeutic implications. Unlike many other cognitive reserve factors (e.g., years of education, prior job complexity), sleep is a modifiable factor. As such, it represents an intervention possibility that may aid the preservation of cognitive function in the face of AD pathology, both present moment and longitudinally.

60 APPLIED LIFE SCIENCES↗

Shedding Light on Primary Donors in Photosynthetic Reaction Centers

Chlorophylls (Chl)s exist in a variety of flavors and are ubiquitous in both the energy and electron transfer processes of photosynthesis. The functions they perform often occur on the ultrafast (fs–ns) time scale and until recently, these have been difficult to measure in real time. Further, the complexity of the binding pockets and the resulting protein-matrix effects that alter the respective electronic properties have rendered theoretical modeling of these states difficult. Recent advances in experimental methodology, computational modeling, and emergence of new reaction center (RC) structures have renewed interest in these processes and allowed researchers to elucidate previously ambiguous functions of Chls and related pheophytins. This is complemented by a wealth of experimental data obtained from decades of prior research. Studying the electronic properties of Chl molecules has advanced our understanding of both the nature of the primary charge separation and subsequent electron transfer processes of RCs. In this review, we examine the structures of primary electron donors in Type I and Type II RCs in relation to the vast body of spectroscopic research that has been performed on them to date. Further, we present density functional theory calculations on each oxidized primary donor to study both their electronic properties and our ability to model experimental spectroscopic data. This allows us to directly compare the electronic properties of hetero- and homodimeric RCs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Understanding process-microstructure-property relationships in laser powder bed fusion of non-spherical Ti-6Al-4V powder

Powder feedstock is a major cost driver in metal additive manufacturing (AM). Replacing the spherical powder with the cost-efficient non-spherical one can reduce the feedstock cost up to 50% and attract more interest to adopt AM in production and new alloy development. Here, for this paper, a comprehensive study was conducted to understand process-microstructure-property relationships in laser powder bed fusion of hydride-dehydride Ti-6Al-4V powder. We demonstrated that variation of laser scan speed had a significant impact on the grain structure, pore evolution and properties compared to laser power. Dynamic X-ray radiography showed that with decreasing scan speed at a constant laser power, a transition from conduction to keyhole mode laser processing occurred, in which a deeper melt pool at lower scan speed intensified texture. In other words, an increase in laser scan speed resulted in formation of the refined prior β grains with shape factor of ~5, lowering the anisotropy. Furthermore, the degree of variant selection was evaluated based on the analyzed texture as a function of laser power and scan speed. With increasing laser scan speed, the dominant α/α boundary type was altered from type 2 to 4 and the degree of variant selection was noticeably decreased. On the other hand, increasing laser power left the morphology of prior β grains, their size, and the dominant α/α boundary (type 4) unchanged, while the texture and anisotropy were intensified, and the degree of variant selection was slightly decreased. Finally, dependency of surface roughness and microhardness were discussed as a function of laser processing parameters.

variant selection↗

Protective Effect of Phloretin against Hydrogen Peroxide-Induced Oxidative Damage by Enhancing Autophagic Flux in DF-1 Cells

Phloretin (PHL) is a dihydrochalcone flavonoid isolated from the peel and root bark of apples, strawberries, and other plants with antioxidative characteristic. In this study, we aimed to investigate the protective effect and the potential mechanism of PHL on hydrogen peroxide (H2O2)-induced oxidative damage in DF-1 cells. The results showed that PHL exhibited no cytotoxic effect on DF-1 cells at concentration below 20 μM. PHL markedly increased H2O2-reduced cell viability, decreased H2O2-induced apoptosis, as evidenced by reduced apoptosis rate, the upregulation of gene and protein level of Bcl-2, and the downregulation of gene and protein level of Bax and Cleaved caspase3. In addition, PHL reduced H2O2-induced reactive oxygen species (ROS) production and restored antioxidant enzymes activities as well as mitochondrial membrane potential in a dose-dependent manner. Moreover, PHL prior to H2O2 further increased LC3-II level, promoted p62 turnover and improved lysosomal function. Importantly, autophagy inhibitor chloroquine (CQ) reversed the protective effect of PHL, and increased H2O2-induced apoptosis. Furthermore, PHL inhibited the phosphorylation levels of ERK, p38, and JNK. Collectively, these results indicate that PHL could attenuate H2O2-induced oxidative injury and apoptosis by maintaining lysosomal function and promoting autophagic flux, and MAPKs pathway may be involved in this process. Our study provides evidence that PHL could as a new strategy to against oxidative damage in poultry industry.

Song, Dan↗

Heat Loss Correction Factor for Fireball Yield Measurements

Fireball yield calculations are performed using Taylor’s radius equation. However, at early times, the x-ray diffusion process has a significant effect on the size of the fireball, and at late times the shock wave degenerates into an acoustic wave. Thus, earlier research determined that the range of validity for Taylor’s equation was in the scaled time regime of 0.004–0.008 s. So, when performing a yield calculation, the only values that are analyzed are those that lie within that range of validity. While Taylor’s equations offer an approximation of the scaled yield, they also rely on several assumptions. For example, the equations assume that the shock wave is adiabatic, when, in fact, the first light pulse on the nuclear test films, which occurs as the shock wave is forming, suggests that the shock wave is not adiabatic. Furthermore, the amount of heat loss from the shock wave that occurs prior to the scaled time regime of 0.004–0.008 s can change significantly as a function of air density. This may cause the shock wave radius to be smaller relative to the assumed adiabatic shock wave. Therefore, the objective of this project was to determine the correction factor for heat transfer, which should be applied in order to produce a more accurate approximation of the weapon’s yield.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Using Bayesian Methodology to Estimate Liquefied Natural Gas Leak Frequencies

This analysis provides estimates on the leak frequencies of nine components found in liquefied natural gas (LNG) facilities. Data was taken from a variety of sources, with 25 different data sets included in the analysis. A hierarchical Bayesian model was used that assumes that the log leak frequency follows a normal distribution and the logarithm of the mean of this normal distribution is a linear function of the logarithm of the fractional leak area. This type of model uses uninformed prior distributions that are updated with applicable data. Separate models are fit for each component listed. Five order-of-magnitude fractional leak areas are considered, based on the flow area of the component. Three types of supporting analyses were performed: sensitivity of the model to the data set used, sensitivity of the leak frequency estimates to differences in the model structure or prior distributions, and sufficiency of sample sized used for convergence. Recommended leak frequency distributions for all component types and leak sizes are given. These leak frequency predictions can be used for quantitative risk assessments in the future.

03 NATURAL GAS↗

Dynamically Learning Incentives for Load Control

As electrical generation becomes more distributed and volatile, and loads become more uncertain, controllability of distributed energy resources (DERs), regardless of their ownership status, will be necessary for grid reliability. Grid operators lack direct control over end-users' grid interactions, such as energy usage, but incentives can influence behavior -- for example, an end-user that receives a grid-driven incentive may adjust their consumption or expose relevant control variables in response. A key challenge in studying such incentives is the lack of data about human behavior, which usually motivates strong assumptions, such as distributional assumptions on compliance or rational utility-maximization. In this paper, we propose a general incentive mechanism in the form of a constrained optimization problem -- our approach is distinguished from prior work by modeling human behavior (e.g., reactions to an incentive) as an arbitrary unknown function. We propose feedback-based optimization algorithms to solve this problem that each leverage different amounts of information and/or measurements. We show that each converges to an asymptotically stable incentive with (near)-optimality guarantees given mild assumptions on the problem. Finally, we evaluate our proposed techniques in voltage regulation simulations on standard test beds. We test a variety of settings, including those that break assumptions required for theoretical convergence (e.g., convexity, smoothness) to capture realistic settings. In this evaluation, our proposed algorithms are able to find near-optimal incentives even when the reaction to an incentive is modeled by a theoretically difficult (yet realistic) function.

demand response↗

A new approach to the evaluation and solution of the relativistic kinetic dispersion relation and verification with continuum kinetic simulation

Here, the present work describes a new approach to evaluation and root finding for the kinetic dispersion relation of Langmuir waves, which is central to the analytical understanding of collisionless damping in plasmas. The plasma dispersion function is solved to machine precision using direct integration in the complex plane in combination with an analytic evaluation of the residue to account for the deformation along the Landau contour. To efficiently attain machine precision, the contour is displaced in the complex plane prior to integration, and numerical subtleties related to the placement of the contour are discussed. The approach is generic in that it applies to arbitrary distribution functions, with the present manuscript focused on relativistic cases. Detailed verification of results via direct kinetic simulation in a variety of configuration space dimensions is also presented. Finally, the technique is applied to the challenging case of highly relativistic (i.e. extremely hot) plasmas. Here we show both qualitative agreement with prior work, as well as the disappearance of the Landau root which would have significant implication for real-life observation or experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Spike-and-Slab Shrinkage Priors for Structurally Sparse Bayesian Neural Networks

Network complexity and computational efficiency have become increasingly significant aspects of deep learning. Sparse deep learning addresses these challenges by recovering a sparse representation of the underlying target function by reducing heavily overparameterized deep neural networks. Specifically, deep neural architectures compressed via structured sparsity (e.g., node sparsity) provide low-latency inference, higher data throughput, and reduced energy consumption. In this article, we explore two well-established shrinkage techniques, Lasso and Horseshoe, for model compression in Bayesian neural networks (BNNs). To this end, we propose structurally sparse BNNs, which systematically prune excessive nodes with the following: 1) spike-and-slab group Lasso (SS-GL) and 2) SS group Horseshoe (SS-GHS) priors, and develop computationally tractable variational inference, including continuous relaxation of Bernoulli variables. We establish the contraction rates of the variational posterior of our proposed models as a function of the network topology, layerwise node cardinalities, and bounds on the network weights. Furthermore, we empirically demonstrate the competitive performance of our models compared with the baseline models in prediction accuracy, model compression, and inference latency.

97 MATHEMATICS AND COMPUTING↗

Validation of the Cossee–Arlman mechanism for propylene oligomerization on Ni/UiO-66

Steady state rate expressions can be derived to distinguish the Cossee–Arlman and metallacycle mechanisms postulated for propylene oligomerization on nickel-based catalysts based on product selectivities, where product selectivities for the former are a function of olefin pressure because sequential coordination and insertion steps lead to independent mechanistic pathways for different hexene isomers. In contrast, the metallacycle mechanism presents pressure-independent product selectivities due to successive coordination prior to the kinetically relevant steps in each mechanism. In this work, steady state propylene oligomerization rates and selectivities were measured in the absence of an activator on nickel functionalized UiO-66 metal organic framework (MOF), Ni/UiO-66, to validate the Cossee–Arlman mechanism for light olefin oligomerization. In situ NO titrations reveal that ~5% of nickel sites were active during the reaction, and thus, not all nickel sites are relevant for catalysis. Propylene dimerization was first order in propylene pressure from 5 to 500 kPa with an apparent activation energy of ~20 kJ mol -1 from 453 to 493 K. Calculated apparent activation energies with density functional theory (DFT) calculations on cluster models of Ni/UiO-66 are in agreement with experiment to corroborate the Cossee–Arlman mechanism. Selectivities of hexene products and the ratio of hexene product selectivities on Ni/UiO-66 are in accordance with selectivity expressions derived from the Cossee–Arlman mechanism. In conclusion, analysis of product selectivities can be used more extensively to demarcate the Cossee–Arlman and metallacycle mechanisms for olefin oligomerization on metal-based catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

\( \mathcal{N} \) = 4 supersymmetric Yang-Mills thermodynamics to order λ 2

We calculate the resummed perturbative free energy of \( \mathcal{N} \) = 4 supersymmetric Yang-Mills in four spacetime dimensions (SYM 4 , 4 ) through second order in the ’t Hooft coupling λ at finite temperature and zero chemical potential. Our final result is ultraviolet finite and all infrared divergences generated at three-loop level are canceled by summing over SYM 4 , 4 ring diagrams. Non-analytic terms at \( \mathcal{O} \) ( λ 3 / 2 ) and \( \mathcal{O} \) ( λ 2 log λ ) are generated by dressing the A 0 and scalar propagators. The gauge-field Debye mass m D and the scalar thermal mass M D are determined from their corresponding finite-temperature self-energies. Based on this, we obtain the three-loop thermodynamic functions of SYM 4 , 4 to \( \mathcal{O} \) ( λ 2 ). We compare our final result with prior results obtained in the weak- and strong-coupling limits and construct a generalized Padé approximant that interpolates between the weak-coupling result and the large- N c strong-coupling result. Our results suggest that the \( \mathcal{O} \) ( λ 2 ) weak-coupling result for the scaled entropy density is a quantitatively reliable approximation to the scaled entropy density for 0 ≤ λ ≲ 2.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Constraining Galaxy-Halo connection using machine learning

We investigate the potential of machine learning (ML) methods to model small-scale galaxy clustering for constraining Halo Occupation Distribution (HOD) parameters. Our analysis reveals that while many ML algorithms report good statistical fits, they often yield likelihood contours that are significantly biased in both mean values and variances relative to the true model parameters. This highlights the importance of careful data processing and algorithm selection in ML applications for galaxy clustering, as even seemingly robust methods can lead to biased results if not applied correctly. ML tools offer a promising approach to exploring the HOD parameter space with significantly reduced computational costs compared to traditional brute-force methods if their robustness is established. Using our ANN-based pipeline, we successfully recreate some standard results from recent literature. Properly restricting the HOD parameter space, transforming the training data, and carefully selecting ML algorithms are essential for achieving unbiased and robust predictions. Among the methods tested, artificial neural networks (ANNs) outperform random forests (RF) and ridge regression in predicting clustering statistics, when the HOD prior space is appropriately restricted. We demonstrate these findings using the projected two-point correlation function (w p (r p )), angular multipoles of the correlation function (ξ ℓ (r)), and the void probability function (VPF) of Luminous Red Galaxies from Dark Energy Spectroscopic Instrument mocks. Our results show that while combining w p (r p ) and VPF improves parameter constraints, adding the multipoles ξ 0 , ξ 2 , and ξ 4 to w p (r p ) does not significantly improve the constraints.

cosmology↗

Anisotropic physics-regularized interpretable machine learning of microstructure evolution

Anisotropic Physics-Regularized Interpretable Machine Learning Microstructure Evolution (APRIMME) is a general-purpose machine learning solution for grain growth simulations. In prior work, PRIMME employed a deep neural network to predict site-specific migration as a function of its neighboring sites to model normal, isotropic, grain growth behavior. This work aims to extend this method by incorporating grain boundary misorientation-based grain growth behavior. APRIMME is trained on anisotropic simulations created using the Monte Carlo-Potts (MCP) model. Furthermore, the results of this work are compared statistically using grain radius, number of sides per grain, mean neighborhood misorientations, and the standard deviation of triple junction dihedral angles, and are found to match in most cases. The exceptions are small and seem to be related to two causes: (1) the deterministic model of APRIMME is learning from the stochastic simulations of MCP, which seems to accentuate triple junction behaviors; and, (2) a bias against very small grains is made evident in a quicker decrease in grains than expected at the beginning of an APRIMME simulation. APRIMME is also evaluated for its general ability to capture anisotropic grain growth behavior by first investigating different test case initial conditions, including a circle grain, three grain, and hexagonal grain microstructures.

36 MATERIALS SCIENCE↗

In Situ Monitoring the Nucleation and Growth of Nanoscale CaCO 3 at the Oil–Water Interface

Interfaces can actively control the nucleation kinetics, orientations, and polymorphs of calcium carbonate (CaCO 3 ). Prior studies have revealed that CaCO 3 formation can be affected by the interplay between chemical functional moieties on solid–liquid or air–liquid interfaces as well as CaCO 3 ’s precursors and facets. Yet little is known about the roles of a liquid–liquid interface, specifically an oil–liquid interface, in directing CaCO 3 mineralization which are common in natural and engineered systems. Here, in this study, by using in situ X-ray scattering techniques to locate a meniscus formed between water and a representative oil, isooctane, we successfully monitored CaCO 3 formation at the pliable isooctane–water interface and systematically investigated the pivotal roles of the interface in the formation of CaCO 3 (i.e., particle size, its spatial distribution with respect to the interface, and its mineral phase). Different from bulk solution, ∼5 nm CaCO 3 nanoparticles form at the isooctane–water interface. They stably exist for a long time (36 h), which can result from interface-stabilized dehydrated prenucleation clusters of CaCO 3 . There is a clear tendency for enhanced amounts and faster crystallization of CaCO 3 at locations closer to isooctane, which is attributed to a higher pH and an easier dehydration environment created by the interface and oil. Our study provides insights into CaCO 3 nucleation at an oil–water interface, which can deepen our understanding of pliable interfaces interacting with CaCO 3 and benefit mineral scaling control during energy-related subsurface operation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗