Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fluid”

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 487 records · Page 27

Use of Micro-Gravity Sensors for External Fluid Level Monitoring in Waste and Nuclear Related Applications - 20315

There are a number of applications in nuclear energy and hazardous waste disposal that require monitoring of fluids under extreme environments, including high levels of temperature, pressure, toxicity and radioactivity. Many of these applications will benefit from a monitoring technique that is external and non-invasive. Currently the sensors used are invasive, must reside inside the pressurized vessels and must penetrate the vessel walls, which can create a weakness in the vessel. Additionally, instruments that are used inside such containers must be exceptionally hardened to the environment. Information Systems Laboratories (ISL) has developed an external mass (gravimetric) measuring technique for monitoring nuclear coolant in Small Modular Reactors (SMRs), which will also work for measuring fluid levels in waste tanks, that avoids the problems inherent in invasive sensors. It utilizes a COTS gravitational sensor of unprecedented accuracy, leveraged via proper sensor placement geometry, to detect fluid changes of small amplitude from an outside position, obviating the need to penetrate the vessel. The technique is called Gravisense{sup TM}. ISL has proven via simulation and experiment that this concept can be usefully applied to monitoring fluid levels in both nuclear reactors and large waste tanks. Numerical simulation algorithms were developed to calculate the gravity effect of small changes in water level, which were verified by experiments at the NIST Physical Simulator facility at the Oregon State University. The measured ultralow noise levels of the superconducting gravimeter type which utilizes a Niobium sphere suspended in a magnetic field to attain its phenomenal accuracy, demonstrated that fluid levels in SMRs can be measured at least to within 3 cm. Furthermore, the method can distinguish between a contained leak (from reactor to containment vessel) from an external leak (from reactor to outside of containment). Additionally, simulations of waste canisters that hold spent fuel rods show that the fluid level measuring accuracy can potentially do better than 1 cm accuracy by measuring from below the vessel, and judicious placement of sensors on top of large waste tanks can potentially achieve a very impressive 2 mm measurement accuracy. These encouraging results prove that the Gravisense{sup TM} technique for fluid determination can be very useful in nuclear energy generation, testing, and research, as well as in waste monitoring situations that are difficult to monitor via traditional sensing technology. We believe that the next step should be to test the technique on canisters of the type that are currently storing waste in various DOE locations. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Flow of “stress power-law” fluids between parallel rotating discs with distinct axes

The problem of flow between parallel rotating discs with distinct axes corresponds to the case of flow in an orthogonal rheometer and has been studied extensively for different fluids since the instrument's inception. All the prior studies presume a constitutive prescription of the fluid stress in terms of the kinematical variables. In this paper, we approach the problem from a different perspective, i.e., a constitutive specification of the symmetric part of the velocity gradient in terms of the Cauchy stress. Such an approach ensures that the boundary conditions can be incorporated in a manner quite faithful to real world experiments with the instrument. Interestingly, the choice of the boundary condition is critical to the solvability of the problem for the case of creeping/Stokes flow. Furthermore, when the no-slip condition is enforced at the boundaries, depending on the model parameters and axes offset, the fluid response can show non-uniqueness or unsolvability, features which are absent in a conventional constitutive specification. In case of creeping/Stokes flow with prescribed values of the stress, the fluid response is indeterminate. We also record the response of a particular case of the given “stress power-law” fluid; one that cannot be attained by the conventional power-law fluids.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Validation and parameterization of a novel physics-constrained neural dynamics model applied to turbulent fluid flow

We report, in fluid physics, data-driven models to enhance or accelerate time to solution are becoming increasingly popular for many application domains, such as alternatives to turbulence closures, system surrogates, or for new physics discovery. In the context of reduced order models of high-dimensional time-dependent fluid systems, machine learning methods grant the benefit of automated learning from data, but the burden of a model lies on its reduced-order representation of both the fluid state and physical dynamics. In this work, we build a physics-constrained, data-driven reduced order model for Navier–Stokes equations to approximate spatiotemporal fluid dynamics in the canonical case of isotropic turbulence in a triply periodic box. The model design choices mimic numerical and physical constraints by, for example, implicitly enforcing the incompressibility constraint and utilizing continuous neural ordinary differential equations for tracking the evolution of the governing differential equation. We demonstrate this technique on a three-dimensional, moderate Reynolds number turbulent fluid flow. In assessing the statistical quality and characteristics of the machine-learned model through rigorous diagnostic tests, we find that our model is capable of reconstructing the dynamics of the flow over large integral timescales, favoring accuracy at the larger length scales. More significantly, comprehensive diagnostics suggest that physically interpretable model parameters, corresponding to the representations of the fluid state and dynamics, have attributable and quantifiable impact on the quality of the model predictions and computational complexity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Impact of viscous fingering and permeability heterogeneity on fluid mixing in porous media

Fluid mixing plays a fundamental role in many natural and engineered processes, including groundwater flows in porous media, enhanced oil recovery, and microfluidic lab-on-a-chip systems. Recent developments have explored the effect of viscosity contrast on mixing, suggesting that the unstable displacement of fluids with different viscosities, or viscous fingering, provides a powerful mechanism to increase fluid-fluid interfacial area and enhance mixing. However, existing studies have not incorporated the effect of medium heterogeneity on the mixing rate. Here, we characterize the evolution of mixing between two fluids of different viscosity in heterogeneous porous media. We focus on a practical scenario of divergent-convergent flow in a quarter five spot geometry prototypical of well-driven groundwater flows. We study by means of numerical simulations the impact of permeability heterogeneity and viscosity contrast on the breakthrough curves and mixing efficiency, and we rationalize the nontrivial mixing behavior that emerges from the competition between the creation of fluid-fluid interfacial area and channeling.

58 GEOSCIENCES↗

On the dynamics and kinematics of two-fluid-phase flow in porous media

A model formulated in terms of both conservation and kinematic equations for phases and interfaces in two-fluid-phase flow in a porous medium system is summarized. Macroscale kinematic equations are derived as extensions of averaging theorems and do not rely on conservation principles. Models based on both conservation and kinematic equations can describe multiphase flow with varying fidelity.When only phase-based equations are considered, a model similar in form to the traditional model for two-fluid-phase flow results. When interface conservation and kinematic equations are also included, a novel formulation results that naturally includes evolution equations that express dynamic changes in fluid saturations, pressures, the capillary pressure, and the fluid-fluid interfacial area density in a two-fluid-system. This dynamic equation set is unique to this work, and the importance of the modeled physics is shown through both microfluidic experiments and high-resolution lattice Boltzmann simulations. The validation work shows that the relaxation of interface distribution and shape toward an equilibrium state is a slow process relativeto the time scale typically allowed for a system to approach an apparent equilibrium state based upon observations of fluid saturations and external pressure measurements. Consequently, most pressure-saturation data intended to denote an equilibrium state are likely a sampling from a dynamic system under-going changes of interfacial curvatures that are not typically monitored. The results confirm the importance of kinematic analysis in combination with conservation equations for faithful modeling of system physics.

54 ENVIRONMENTAL SCIENCES↗

Fault slip controlled by stress path and fluid pressurization rate

Abstract The practice of injecting fluids into the crust is linked to regional increases in seismicity. Increasing fluid pressure along preexisting faults is believed to enhance seismicity rates by reducing the shear stress required for slip, but the processes that cause faults to slip under conditions of fluid pressurization are poorly constrained. We use experimental rock deformation to investigate the controls of fluid pressurization and pressurization rates on fault slip style. We show that pore fluid pressurization is less effective that mechanical changes in fault normal stress at initiating accelerated slip events. Fluid pressurization enhances the total slip, slip velocity, and shear stress drop of events initiated by mechanical changes in normal stress, and these parameters are correlated with pressurization rate, but not the magnitude of fluid pressure. This result is consistent with field‐scale observations and indicates that processes active at the pore network scale affect induced seismicity.

French, Melodie E.↗

A comparison of segmented abdominopelvic fluid volumes with conventional CT signs of abdominal compartment syndrome in a trauma population

Purpose: To compare the utility of abdominopelvic fluid volume measurements with established computed tomography signs for refractory post-traumatic abdominal compartment syndrome. Methods: This retrospective observational cohort study included 64 consecutive adult trauma patients with preoperative CT and diagnosis of refractory abdominal compartment syndrome requiring decompressive laparotomy at a level I trauma referral center between 2004 and 2014. We hypothesized that abdominal fluid volume measurements would be more predictive of the need for early laparotomy than previously described conventional CT signs of refractory ACS. Abdominopelvic fluid volumes were determined quantitatively using semi-automated segmentation software. The following conventional imaging parameters were recorded: abdominal anteroposterior:transverse ratio (round belly sign); infrahepatic vena cava diameter; distal abdominal aortic diameter; largest single small bowel wall diameter; hydronephrosis, inguinal herniation; and mesenteric and body wall edema. For outcome analysis, patients were stratified into two groups: those who underwent early (< 24 h) and late (≥ 24 h) decompressive laparotomy following CT. Correlation analysis, comparison of means, and multivariate logistic regression were performed. Results: Abdominal fluid volumes (p = 0.001) and anteroposterior:transverse ratio (p = 0.009) were increased and inferior vena cava diameter (p = 0.009) was decreased in the early decompressive laparotomy group. Multivariate analysis including conventional CT variables, fluid volumes, and laboratory values revealed abdominal fluid volumes (p = 0.012; Δ in log odds of 1.002/mL) as the only independent predictor of early decompressive laparotomy. Conclusions: Segmented abdominopelvic free fluid volumes had greater predictive utility for decision to perform early decompressive laparotomy than previously described ACS-related CT signs in trauma patients who developed refractory abdominal compartment syndrome.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Hard X-ray fluorescence spectroscopy of high pressure cavitating fluids in aluminum nozzles

X-rays are frequently used to study the internal geometry of dense objects, and to measure the density of multiphase flows. However, quantitatively measuring fluid density inside a metallic object such as a high pressure spray nozzle is difficult. X-rays of sufficiently high energy to penetrate a metal object are not appreciably absorbed by the fluid inside. This requires the use of plastic or beryllium test sections, which are not suited to high pressure conditions. We present a high-energy X-ray fluorescence technique which can overcome this problem. The experiments were conducted at the 7-BM beamline of the Advanced Photon Source at Argonne National Laboratory. A hydrocarbon fluid was seeded with cerium nanoparticles. The fluid was pumped at high pressure through aluminum nozzles with inner diameters of 0.36-0.90 mm and wall thicknesses of 2-3 mm. A collimated, monochromatic 42.5 keV X-ray beam excited K-edge fluorescence from the cerium. The K-alpha emission lines at 34-35 keV were recorded by a cryogenic germanium detector. Changes in fluid density due to cavitation of the liquid inside the nozzle were measured by raster scanning the nozzle across the beam. A spatial resolution of 20 x 20 mu m(2) was achieved with a slitted beam, which was improved to 5 x 10 mu m(2) with X-ray focusing mirrors. The uncertainty in the path-integrated vapor fraction was 30-40 mu m at 95% confidence. A limitation of this approach is that for low vapor pressure fluids, the nanoparticles increase the vapor pressure of the fluid and act as additional nucleation sites. These experiments demonstrate a path forward for measurements of multiphase flows inside metal components under conditions that are not feasible in optically accessible materials.

47 OTHER INSTRUMENTATION↗

Evidence for metal sources, fluid-mixing processes, and S isotope recycling within the feeder zone of an Irish type Zn-Pb deposit

The origin and evolution of fluids in Irish-type Zn-Pb deposits remains debated, particularly regarding the mobility of metals such as Cu and Ni, sources of sulphur, and the role of fluid mixing and replacement. The Lisheen Zn-Pb deposit, Ireland, offers a well-defined natural laboratory to investigate these questions. While most studies have focused on the Waulsortian Limestone Formation, the primary sulphide host, less is known about mineralisation in underlying units, such as the Lisduff Oolite Member (LOM). The LOM displays enrichment in Cu and Ni and displays intense replacement textures compared to other hosts at Lisheen, making it an ideal target for studying metal mobility and sulphur recycling in carbonate-hosted systems. Through characterising and studying LOM-hosted sulphides, valuable insights into mineralisation processes, especially related to Cu-Ni metals, can be defined. This study integrates petrography, EMPA, and in situ sulphur isotope (δ 34 S) analysis to investigate sulphide paragenesis, mineral chemistry, and fluid evolution across LOM ore zones. Results reveal a multistage mineralising system involving extensive replacement of early pyrite (Py0, δ 34 S = −28.4 to −21.9 ‰) by sphalerite and galena, with zoned pyrite (Py1) enriched in As-Cu-Ni-Tl. The δ 34 S values and trace element trends indicate mixing between hydrothermal and bacteriogenic sulphur-rich fluids, with evidence for sulphur recycling during replacement. Pyrite textures and compositions capture this evolving fluid regime, with trace element enrichment linked to paragenetic stage. The steel ore region, adjacent to major fault intersections, records intense hydrothermal fluid interaction, hosting Ni- and As-rich phases such as nickeline, gersdorffite, and arsenopyrite. These findings highlight the importance of structural controls and fluid mixing in metal transport and deposition, positioning the LOM as a key stratigraphic unit for understanding ore-forming processes in Irish-type systems. These results have implications for targeting similar carbonate-hosted systems globally, especially where deeper or structurally complex ore zones remain underexplored.

58 - GEOSCIENCES↗

Experimental Investigation of Barium Sources and Fluid–Rock Interaction in Unconventional Marcellus Shale Wells Using Ba Isotopes

Produced waters from unconventional Marcellus Shale gas wells have anomalously high barium (Ba) concentrations and yield some of the isotopically heaviest Ba measured to date. Experiments were conducted to constrain the source of Ba in these fluids and the controls on barite (BaSO 4 ) precipitation and dissolution in oil and gas wells. Experiments simulating the acidizing stage evaluated the solubility of pure barite and drilling mud in 2 M HCl at 80 °C for periods of 2, 6, and 48 h and resulted in <0.01% barite dissolution with no appreciable change in δ 138 Ba ( 138 Ba/ 134 Ba normalized to NIST standard 3104a). Static autoclave experiments conducted at 66 °C and 20.7 MPa with combinations of ground Marcellus Shale solids and/or barite-bearing drilling mud with synthetic low-Ba fracturing fluid resulted in decreased Ba concentrations in the fluid, with the largest decrease in the shale-only run. Fluid δ 138 Ba values increased by about 0.5‰ as Ba concentrations decreased, consistent with closed-system Rayleigh fractionation. Flow-through experiments in Marcellus Shale core conducted for 28 days resulted in effluent Ba concentrations an order of magnitude lower than the influent, while sulfate concentrations increased over time. Effluent δ 138 Ba values increased over the first 12 days and plateaued at about 1‰ higher than the influent. Modeling suggests a combination of the release of labile shale Ba and barite precipitation. This work indicates that the processes of Ba release from fluid–shale interactions and barite precipitation in fractures and the well bore, while capable of producing high δ 138 Ba fluids, are unlikely to generate fluids with high-Ba concentrations and δ 138 Ba values like those in Marcellus-produced waters. As a result, we find that the release of sulfate from shale pyrite oxidation rapidly catalyzes barite precipitation and that dissolution of drilling mud barite or natural barite in the shale is unlikely to be the major source of Ba in Marcellus-produced waters.

54 ENVIRONMENTAL SCIENCES↗

Dynamic Mixing Behaviors of Ionically Tethered Polymer Canopy of Nanoscale Hybrid Materials in Fluids of Varying Physical and Chemical Properties

An emerging area of sustainable energy and environmental research is focused on the development of novel electrolytes that can increase the solubility of target species and improve subsequent reaction performance. Electrolytes with chemical and structural tunability have allowed for significant advancements in flow batteries and CO 2 conversion integrated with CO 2 capture. Liquid-like nanoparticle organic hybrid materials (NOHMs) are nanoscale fluids that are composed of inorganic nanocores and an ionically tethered polymeric canopy. NOHMs have been shown to exhibit enhanced conductivity making them promising for electrolyte applications, though they are often challenged by high viscosity in the neat state. In this study, a series of binary mixtures of NOHM-I-HPE with five different secondary fluids, water, chloroform, toluene, acetonitrile, and ethyl acetate, were prepared to reduce the fluid viscosity and investigate the effects of secondary fluid properties (e.g., hydrogen bonding ability, polarity, and molar volume) on their transport behaviors, including viscosity and diffusivity. We report our results revealed that the molecular ratio of secondary fluid to the ether groups of Jeffamine M2070 ($λ_{\text{SF}}$) was able to describe the effect that secondary fluid has on transport properties. Our findings also suggest that in solution, the Jeffamine M2070 molecules exist in different nanoscale environments, where some are more strongly associated with the nanoparticle surface than others, and the conformation of the polymer canopy was dependent on the secondary fluid. This understanding of the polymer conformation in NOHMs can allow for the better design of an electrolyte capable of capturing and releasing small gaseous or ionic species.

25 ENERGY STORAGE↗

Confinement-mediated phase behavior of hydrocarbon fluids: Insights from Monte Carlo simulations

The phase behavior of hydrocarbon fluids confined in porous media has been discussed to deviate significantly from that in the bulk environment due to the existence of sub-10nm pores. Though experiments and simulations have measured the bubble/dew points and sorption isotherms of hydrocarbons confined in both natural and synthetic nanopores, the confinement effects in terms of the strength of fluid-pore interactions tuned by surface wettability and chemistry have received comparably less discussion. More significantly, the underlying physics of confinement-induced phenomena remain obfuscated. In this work, we studied the phase behavior and capillary condensation of n-hexane to understand the effects of confinement at the molecular level. To systematically investigate the pore effects, we constructed two types of wall confinements; one is a structureless virtual wall described by the Steele potential and the other one is an all-atom amorphous silica structure with surface modified by hydroxyl groups. Our numerical results demonstrated the importance of fluid-pore interaction, pore size, and pore morphology effects in mediating the pressure-volume-temperature (PVT) properties of hydrocarbons. The most remarkable finding of this work was that the saturation pressure predicted from the van der Waals-type adsorption isothermal loop could be elevated or suppressed relative to the bulk phase, as illustrated in the graphical abstract. As the surface energy (i.e., fluid-pore interaction) decreased, the isothermal vapor pressure increased, indicating a greater preference for the fluid to exist in the vapor state. Sufficient reduction of the fluid-pore interactions could even elevate the vapor pressure above that of the bulk fluid.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Room-temperature polariton quantum fluids in halide perovskites

Abstract Quantum fluids exhibit quantum mechanical effects at the macroscopic level, which contrast strongly with classical fluids. Gain-dissipative solid-state exciton-polaritons systems are promising emulation platforms for complex quantum fluid studies at elevated temperatures. Recently, halide perovskite polariton systems have emerged as materials with distinctive advantages over other room-temperature systems for future studies of topological physics, non-Abelian gauge fields, and spin-orbit interactions. However, the demonstration of nonlinear quantum hydrodynamics, such as superfluidity and Čerenkov flow, which is a consequence of the renormalized elementary excitation spectrum, remains elusive in halide perovskites. Here, using homogenous halide perovskites single crystals, we report, in both one- and two-dimensional cases, the complete set of quantum fluid phase transitions from normal classical fluids to scatterless polariton superfluids and supersonic fluids—all at room temperature, clear consequences of the Landau criterion. Specifically, the supersonic Čerenkov wave pattern was observed at room temperature. The experimental results are also in quantitative agreement with theoretical predictions from the dissipative Gross-Pitaevskii equation. Our results set the stage for exploring the rich non-equilibrium quantum fluid many-body physics at room temperature and also pave the way for important polaritonic device applications.

42 ENGINEERING↗

Fracture-matrix fluid exchange in oil-bearing unconventional mudstones

The poromechanical properties of unconventional reservoir materials are in large part dictated by their mineralogy. Since these properties govern the response to stress experienced during hydraulic fracturing, fluid production, and fluid injection, they play a central role in the formation of microcracks or bedding delaminations which ultimately dominate mass transport. In this work we study access to the porosity of end member unconventional reservoir materials, where the end members are predominantly dictated by carbonate content. Access to the porosity is quantified using state of the art 3D x-ray computed tomography coupled with physics informed data analytics. Xenon gas, which attenuates x-rays, provides a spatiotemporal map of access to the porosity. The accessible porosity is quantified over a range of net confining stress relevant to the manmade disturbances listed above. These experiments demonstrate that heavily carbonated mudstones are nearly impermeable at the core (~ cm) scale, while carbonate free analogues afford better access to the microstructure. Consistent with previous qualitative 2D radiographs, access to the interior of the clastic mudstones is first observed along planar microcracks, followed by slow penetration into the surrounding matrix. Physics informed data analytics of the 3D tomography measurements presented here show that these microcracks do not permit uniform access to the adjacent rock matrix. In addition, variation of the effective pressure elucidates the mechanisms that govern fracture/matrix fluid exchange. Under conditions consistent with hydrocarbon production fluid accumulates in the immediate vicinity of the nearest microcrack. While there is clear evidence that, as intended, part of this accumulation is from the more distant matrix, fluid is also squeezed out of the microcrack. The fluid build-up at the microcrack indicates that migration out of the rock is hindered by the coupled poroelastic response of the microcrack and adjacent rock matrix. We show that these mechanisms ultimately account for the meager oil recovery factors realized in practice. These insights have implications for making reservoir scale predictions based on core scale observations, and provide a basis for devising new asset development techniques to access more porosity, and enhance fluid extraction. Finally, these findings shed light on key features and mechanisms that govern shale storage capacity, with relevance to other important industrial processes, such as geologic CO 2 storage.

58 GEOSCIENCES↗

The role of fluid-wall interactions on confined liquid diffusion using Mori theory

The dynamics of fluid confined in a nano-channel with smooth walls have been studied through velocity autocorrelation function within the memory function approach by incorporating the atomic level interactions of fluid with the confining wall. Expressions for the second and fourth sum rules of velocity autocorrelation have been derived for nano-channel which involves fluid-fluid and fluid-wall interactions. These expressions, in addition, involve pair correlation function and density profiles. The numerical contributions of fluid-wall interaction to sum rules are found to play a very significant role, specifically at smaller channel width. Results obtained for velocity autocorrelation and self-diffusion coefficient of a fluid confined to different widths of the nanochannel have been compared with the computer simulation results. The comparison shows a good agreement except when the width of the channel is of the order of two atomic diameters, where it becomes difficult to estimate sum rules involving the triplet correlation’s contribution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Forward and backward in time dispersion of fluid and inertial particles in isotropic turbulence

In this paper, we investigate both theoretically and numerically the Forward-In-Time (FIT) and Backward-In-Time (BIT) dispersion of fluid and inertial particle-pairs in isotropic turbulence. Fluid particles are known to separate faster BIT than FIT in three-dimensional turbulence, and we find that inertial particles do the same. However, we find that the irreversibility in the inertial particle dispersion is in general much stronger than that for fluid particles. For example, the ratio of the BIT to FIT mean-square separation can be up to an order of magnitude larger for the inertial particles than for the fluid particles. We also find that for both the inertial and fluid particles, the irreversibility becomes stronger as the scale of their separation decreases. Regarding the physical mechanism for the irreversibility, we argue that whereas the irreversibility of fluid particle-pair dispersion can be understood in terms of a directional bias arising from the energy transfer process in turbulence, inertial particles experience an additional source of irreversibility arising from the non-local contribution to their velocity dynamics, a contribution that vanishes in the limit St → 0, where St is the particle Stokes number. For each given initial (final, in the BIT case) separation, r 0 , there is an optimum value of St for which the dispersion irreversibility is strongest, as such particles are optimally affected by both sources of irreversibility. We derive analytical expressions for the BIT, mean-square separation of inertial particles and compare the predictions with numerical data obtained from a Re λ ≈ 582 (where Re λ is the Taylor Reynolds number) Direct Numerical Simulation (DNS) of particle-laden isotropic turbulent flow. The small-time theory, which in the dissipation range is valid for times ≤max[Stτ η , τ η ] (where τ η is the Kolmogorov time scale), is in excellent agreement with the DNS. The theory for long-times is in good agreement with the DNS provided that St is small enough so that the inertial particle motion at long-times may be considered as a perturbation about the fluid particle motion, a condition that would in fact be satisfied for arbitrary St at sufficiently long-times in the limit Re λ → ∞.

42 ENGINEERING↗

Continuum kinetic and multi-fluid simulations of classical sheaths

The kinetic study of plasma sheaths is critical, among other things, to understand the deposition of heat on walls, the effect of sputtering, and contamination of the plasma with detrimental impurities. The plasma sheath also provides a boundary condition and can often have a significant global impact on the bulk plasma. In this paper, kinetic studies of classical sheaths are performed with the continuum kinetic code, Gkeyll, which directly solves the Vlasov-Maxwell equations. The code uses a novel version of the finite-element discontinuous Galerkin scheme that conserves energy in the continuous-time limit. The fields are computed using Maxwell equations. Ionization and scattering collisions are included; however, surface effects are neglected. The aim of this work is to introduce the continuum kinetic method and compare its results with those obtained from an already established finite-volume multi-fluid model also implemented in Gkeyll. Novel boundary conditions on the fluids allow the sheath to form without specifying wall fluxes, so the fluids and fields adjust self-consistently at the wall. Our work demonstrates that the kinetic and fluid results are in agreement for the momentum flux, showing that in certain regimes, a multifluid model can be a useful approximation for simulating the plasma boundary. There are differences in the electrostatic potential between the fluid and kinetic results. Further, the direct solutions of the distribution function presented here highlight the non-Maxwellian distribution of electrons in the sheath, emphasizing the need for a kinetic model. The densities, velocities, and the potential show a good agreement between the kinetic and fluid results. But, kinetic physics is highlighted through higher moments such as parallel and perpendicular temperatures which provide significant differences from the fluid results in which the temperature is assumed to be isotropic. Besides decompression cooling, the heat flux is shown to play a role in the temperature differences that are observed, especially inside the collisionless sheath. Published by AIP Publishing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A molecular Debye-Huckel theory of solvation in polar fluids: An extension of the Born model

A dielectric response theory of solvation beyond the conventional Born model for polar fluids is presented. The dielectric response of a polar fluid is described by a Born response mode and a linear combination of Debye-Hückel-like response modes that capture the nonlocal response of polar fluids. The Born mode is characterized by a bulk dielectric constant, while a Debye-Hückel mode is characterized by its corresponding Debye screening length. Both the bulk dielectric constant and the Debye screening lengths are determined from the bulk dielectric function of the polar fluid. The linear combination coefficients of the response modes are evaluated in a self-consistent way and can be used to evaluate the electrostatic contribution to the thermodynamic properties of a polar fluid. In conclusion, our theory is applied to a dipolar hard sphere fluid as well as interaction site models of polar fluids such as water, where the electrostatic contribution to their thermodynamic properties can be obtained accurately.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗