Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coupling constant”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Multiple frequency radar observations of high-latitude E region irregularities in the HF modified ionosphere

This paper discusses the results of the September 1983 observations of artificial field-aligned irregularities (AFAIs) in the Tromso, Norway region, made by backscatter radars operating at 46.9, 143.8, 21.4, and 140.0 MHz. Four classes of resonant instability processes at work in the E and F regions are examined in detail: (1) the coupling of parametric decay instability waves across geomagnetic field lines, (2) thermal parametric instability, (3) four-wave interaction thermal parametric instability, and (4) the resonance instability. The characteristics of the AFAI scatter are described, with special attention given to the growth and decay time constants, functional dependence on the heater power and polarization, and the scattering cross sections of the irregularities.

Noble, S. T.↗
NASA NTRS2004

Anisotropic Piezocomposite Actuator Incorporating Machined PMN-PT Single Crystal Fibers

The design, fabrication, and testing of a flexible, planar, anisotropic piezoelectric composite actuator utilizing machined PMN-32%PT single crystal fibers is presented. The device consists of a layer of rectangular single crystal piezoelectric fibers in an epoxy matrix, packaged between interdigitated electrode polyimide films. Quasistatic free-strain measurements of the single crystal device are compared with measurements from geometrically identical specimens incorporating polycrystalline PZT-5A and PZT-5H piezoceramic fibers. Free-strain actuation of the single crystal actuator at low bipolar electric fields (+/- 250 V/mm) is approximately 400% greater than that of the baseline PZT-5A piezoceramic device, and 200% greater than that of the PZT-5H device. Free-strain actuation under high unipolar electric fields (0-4kV/mm) is approximately 200% of the PZT-5A baseline device, and 150% of the PZT-5H alternate piezoceramic device. Performance increases at low field are qualitatively consistent with predicted increases based on scaling the low-field d(sub 33) piezoelectric constants of the respective piezoelectric materials. High-field increases are much less than scaled d(sub 33) estimates, but appear consistent with high-field freestrain measurements reported for similar bulk single-crystal and piezoceramic compositions. Measurements of single crystal actuator capacitance and coupling coefficient are also provided. These properties were poorly predicted using scaled bulk material dielectric and coupling coefficient data. Rules-of-mixtures calculations of the effective elastic properties of the single crystal device and estimated actuation work energy densities are also presented. Results indicate longitudinal stiffnesses significantly lower (50% less) than either piezoceramic device. This suggests that single-crystal piezocomposite actuators will be best suited to low induced-stress, high strain and deflection applications.

Wilkie, W. Keats↗
NASA NTRS2016

Design of a Two-Stage Light Gas Gun for Muzzle Velocities of 10 - 11 km/s

Space debris poses a major risk to spacecraft. In low earth orbit, impact velocities can be 10 - 11 km/s and as high as 15 km/s. For debris shield design, it would be desirable to be able to launch projectiles of known shape and mass to these velocities. The design of the proposed 10 - 11 km/sec gun uses, as a starting point, the Ames 1.28”/0.22” two stage gun, which has achieved muzzle velocities of 10 - 11.3 km/sec. That gun is scaled up to a 0.3125” launch tube diameter. The gun is then optimized with respect to maximum pressures by varying the pump tube length to diameter ratio (L/D), the piston mass and the hydrogen pressure. A pump tube L/D of 36.4 is selected giving the best overall performance. Piezometric ratios for the optimized guns are found to be ~2.3, much more favorable than for more traditional two stage light gas guns, which range from 4 to 6. (The piezometric ratio for a gun is defined as the maximum projectile base pressure divided by the constant projectile base pressure which, acting over the entire barrel length, would produce the same muzzle velocity.) The maximum powder chamber pressures are 20 to 30 ksi. To reduce maximum pressures, the desirable range of the included angle of the cone of the high pressure coupling is found to be 7.3 to 14.6 degrees. Lowering the break valve rupture pressure is found to lower the maximum projectile base pressure, but to raise the maximum gun pressure. For the optimized gun with a pump tube L/D of 36.4, increasing the muzzle velocity by decreasing the projectile mass and increasing the powder loads is studied. It appears that saboted spheres could be launched to 10.25 and possibly as high as 10.8 km/sec, and that disc-like plastic models could be launched to 11.05 km/s. The use of a tantalum liner to greatly reduce bore erosion and increase muzzle velocity is discussed. With a tantalum liner, CFD code calculations predict muzzle velocities as high as 12 to 13 km/s.

Bogdanoff, David W.↗
DOE OSTI2026

Integration of Online Cross-Section Generation Capability with Depletion and Transient Solvers in Griffin

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE)-based reactor multiphysics analysis application jointly developed by Argonne and Idaho National Laboratories under the DOENE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. In FY25, an online crosssection generation capability based on the Self-Shielding Application Programming Interface (SSAPI) was demonstrated for TRISO-fueled reactor problems under steady-state conditions. This fiscal year, that capability was extended to support depletion and transient multiphysics calculations, enabling high-fidelity analyses that generate self-shielded cross sections on the fly from the actual evolving composition and temperature states rather than from pre-tabulated libraries. For depletion, a two-way coupling was established in which SSAPI computes compact-averaged self-shielded cross sections that the depletion solver then uses to advance the Bateman equations, with the updated compositions returned to SSAPI at each step; the depletion module was refactored to support both library-based and SSAPI-based cross sections, and additional logic was added to track daughter isotopes and to exclude minor isotopes for efficiency. For transient analysis, the SSAPI multigroup library was extended with the kinetics data required for time-dependent calculations, the Improved Quasi-Static (IQS) scheme was coupled with SSAPI, and several supporting capabilities were implemented, including a self-shielding treatment that lets control rods and drums move within a self-shielded model, which had previously been impossible and had ruled out rod- and drum-movement transients with on-the-fly cross sections altogether, a new mixing scheme for delayed-neutron precursor decay constants, a checkpoint-based restart workflow, and performance improvements such as pointwise cross-section interpolation and the bypassing of unnecessary Dancoff factor calculations. The implemented capabilities were verified against Serpent Monte Carlo solutions. For depletion, a prismatic pin-cell problem based on a Next Generation Nuclear Plant (NGNP) Very High Temperature Reactor benchmark showed excellent agreement, with eigenvalue differences within 200 pcm over the entire burnup range (up to 140 MWD/kgU) and fission-product and actinide inventories agreeing to within 0.8% and 2.5%, respectively; a heat-pipe microreactor assembly problem with a much higher fuel loading confirmed the same behavior and quantified the bias introduced when the multigroup equivalence effect is neglected. For transient analysis, a pin-cell problem with a step reactivity insertion and temperature feedback reproduced the analytically expected asymptotic power and showed close agreement between the direct and IQS solutions, and a two-dimensional microreactor core problem with control-drum rotation exercised the new moving-drum self-shielding treatment and demonstrated successful coupling of the online crosssection generation with both the direct and IQS transient methods. The capability was further exercised on a full-core pebble-bed problem, in which Griffin was coupled with the System Analysis Module (SAM) to simulate load-following operation of the gPBR with the Doppler feedback resolved at the TRISO fuel kernel temperature. These developments in Griffin provide a convenient, high-fidelity approach to cross-section generation for advanced thermal reactors with geometrically complex and highly heterogeneous configurations, including TRISO-fueled prismatic and pebble-bed systems, and support steady-state, depletion, and transient multiphysics calculations. They also enable self-shielded cross sections to be evaluated directly at the actual coupled state of the system, thereby establishing a foundation for high-fidelity, fully coupled multiphysics analysis of advanced reactors

Park, H.↗
NASA NTRS1983

HCN formation on Jupiter - The coupled photochemistry of ammonia and Acetylene

HCN formation in the upper troposphere and lower stratosphere of Jupiter is presently modeled in terms of UV pyrolysis of the C2H5N isomer aziridine, which is a product of the NH2 and C2H3 radicals that originate from ammonia photolysis and the addition of H atoms to acetylene, respectively. The sensitivity of the HCN column density to the individual rate constants and the eddy diffusion coefficient profile is considered, along with the possibility that additional HCN-yielding pathways may exist. Both ammonia and phosphine are strongly depleted by photolysis.

Kaye, J. A.↗
NASA NTRS1998

Acquisition and Early Losses of Rare Gases from the Deep Earth

Direct observations show that the deep Earth contains rare gases of solar composition distinct from those in the atmosphere. We examine the implications of mantle rare gas characteristics on acquisition of rare gases from the solar nebula and subsequent losses due to a large impact. Deep mantle rare gas concentrations and isotopic compositions can be obtained from a model of transport and distribution of mantle rare gases. This model assumes the lower mantle closed early, while the upper mantle is open to subduction from the atmosphere and mass transfer from the lower mantle. Constraints are derived that can be incorporated into models for terrestrial volatile acquisition: (1) Calculated lower-mantle Xe-isotopic ratios indicate that the fraction of radiogenic Xe produced by I-129 and Pu-244 during the first about 10(exp 8) yr was lost, a conclusion also drawn for atmospheric Xe. Thus, either the Earth was made from materials that had lost >99% of rare gases about (0.7-2) x 10(exp 8) yr after the solar system formed, or gases were then lost from the fully formed Earth. (2) Concentrations of 3He and 20Ne in the lower mantle were established after these losses. (3) Neon-isotopic data indicates that mantle Ne has solar composition. The model allows for solar Ar/Ne and Xe/Ne in the lower mantle if a dominant fraction of upper mantle Ar and Xe are subduction-derived. If Earth formed in the presence of the solar nebula, it could have been melted by accretional energy and the blanketing effect of a massive, nebula-derived atmosphere. Gases from this atmosphere would have been sequestered within the molten Earth by dissolution at the surface and downward mixing. It was found that too much Ne would be dissolved in the Earth unless the atmosphere began to escape when the Earth was only partially assembled. Here we consider conditions required to initially dissolve sufficient rare gases to account for the present lower mantle concentrations after subsequent losses at 10(exp 8) yr. It is assumed that equilibration of the atmosphere with a thoroughly molten mantle was rapid, so that initial abundances of gases retained in any mantle layer reflected surface conditions when the layer solidified. For subsequent gas loss of 99.5% and typical solubility coefficients, a total pressure of 100 atm was required for an atmosphere of solar composition. Calculations of the pressure at the base of a primordial atmosphere indicate that this value might be exceeded by an order of magnitude or more for an atmosphere supported by accretional energy. Surface temperatures of about 4000 K would have been produced, probably high enough to melt the deep mantle. Initial distributions of retained rare gases would then be determined by the history of surface pressure and temperature during mantle cooling and solidification, i.e., the coupled cooling of Earth and atmosphere. The Earth's thermal state was determined by its surface temperature and the efficiency of convection in the molten mantle, estimated to be sufficient to maintain an adiabatic gradient. Because the melting curve is steeper than the adiabat, solidification of the mantle proceeded outward from the interior. Incorporation of atmospheric gases in the mantle therefore occurred over a range in surface temperature of a few thousand degrees Kelvin. The thermal state of the atmosphere was controlled by total luminosity of the Earth (energy) released by accreting planetesimals and the cooling Earth), nebular temperature and pressure, and atmospheric opacity. The energy released by accretion declined with time as did nebular pressure. Analytical solutions for an idealized (constant opacity radiative atmosphere show that declining energy sources under constant nebular conditions result in slowly diminishing surface temperature but dramatically increasing surface pressure. For such an atmosphere with declining nebular pressure but constant total luminosity, surface pressure decreases gradually with decreasing temperaure. A decline in accretion luminosity might be compensated by energy released as the mantle cools for about 10(exp 5) year, after which luminosity must decline. The total complement of dissolved rare gases will depend on the particular evolutionary path determined by the declining accretional luminosity, the Earth thermal history, removal of the nebula, and opacity variations of the atmosphere. Models for these coupled evolutionary histories for Earth's acquisition of nebular-derived noble gases are in progress. The later losses required at about 10(exp 8) yr (depleting the interior concentrations of the sequestered solar gases by a factor of > 100) were presumably related to the major impact in which the Moon formed.

Porcelli, D.↗
DOE OSTI2024

Bioelectrochemical reactor to manage anthropogenic sulfate pollution for freshwater ecosystems: Mathematical modeling and experimental validation

Anthropogenic sulfate loading into otherwise low-sulfate freshwater systems can cause significant ecological consequences as a biogeochemical stressor. To address this challenge, in situ bioremediation technologies have been developed to leverage naturally occurring microorganisms that transform sulfate into sulfide rather than implementing resource-intensive physio-chemical processes. However, bioremediation technologies often require the supply of electron donors to facilitate biological sulfate reduction. Bioelectrochemical systems (BES) can be an alternative approach for supplying molecular hydrogen as an electron donor for sulfate-reducing bacteria through water electrolysis. Although the fundamental mechanisms behind BESs have been studied, limited research has evaluated the design and operational parameters of treatment systems when developing BESs on a scale relevant to environmental systems. Further, this study aimed to develop an application-based mathematical model to evaluate the performance of BESs across a range of reactor configurations and operational modes. The model was based on sulfate transformation by hydrogenotrophic sulfate-reducing bacteria coupled with the recovery of solid iron sulfide species formed by the oxidative dissolution of dissolved ferrous iron from a stainless steel anode. Sulfate removal closely corresponded to the rate of electrolytic hydrogen production and hydraulic residence time but was less sensitive to specific microbial rate constants. The mathematical model results were compared to experimental data from a pilot-scale BES tested with nonacidic mine drainage as a case study. The close agreement between the mathematical model and the pilot-scale BES experiment highlights the efficacy of using a mathematical model as a tool to develop a conceptual design of a scaled-up treatment system.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Bimetallic Metal–Organic Framework Fe/Co-MIL-88(NH 2 ) Exhibiting High Peroxidase-like Activity and Its Application in Detection of Extracellular Vesicles

Metal–organic frameworks (MOFs) have many attractive features, including tunable composition, rigid structure, controllable pore size, and large specific surface area, and thus are highly applicable in molecular analysis. Depending on the MOF structure, a high number of un-saturated metal sites can be exposed to catalyze chemical reactions. In the present work, we report that by using both Co(II) and Fe(III) to prepare the MIL-88(NH 2 ) MOF, we can produce the bimetallic MOF that can catalyze the conversion of 3,3', 5,5"-tetramethylbenzidine (TMB) to a color product through reaction with H 2 O 2 at a higher reaction rate than the monometallic Fe-MIL-88(NH 2 ). The Michaelis constants (K m ) of the catalytic reaction for TMB and H 2 O 2 are 3-5 times smaller, and the catalytic constants (k cat ) are 5-10 times higher than those of the horse-radish peroxidase (HRP), supporting ultrahigh peroxidase-like activity. These values are also much more superior to those of the HRP-mimicking MOFs reported previously. Interestingly, the bimetallic MOF can be coupled with glucose oxidase (GOx) to trigger the cascade enzymatic reaction for highly sensitive detection of extracellular vesicles (EVs), a family of important biomarkers. Through conjugation to the aptamer that recognizes the marker protein on EV surface, the MOF can help isolate the EVs from biological matrices, which are subsequently labeled by GOx via antibody recognition. The cascade enzymatic reaction between MOF and GOx enables detection of EVs at a concentration as low as 7.8 × 10 4 particles/ml. Further, the assay can be applied to monitor EV secretion by cultured cells, and also can successfully detect the different EV quantities in the sera samples collected from cancer patients and healthy controls. Overall, we prove that the bimetallic Fe/Co-MIL-88(NH 2 ) MOF, with its high peroxidase-activity and high biocompatibility, is a valuable tool deployable in clinical assays to facility disease diagnosis and prognosis.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Anyon superconductivity from topological criticality in a Hofstadter–Hubbard model

We argue that the combination of strong repulsive interactions and high magnetic fields can generate electron pairing and superconductivity. Inspired by the large lattice constants of moiré materials, which make large flux per unit cell accessible at laboratory fields, we study the triangular lattice Hofstadter-Hubbard model at one-quarter flux quantum per plaquette, where previous literature has argued that a chiral spin liquid separates a weak-coupling integer quantum Hall phase and a strong-coupling topologically trivial antiferromagnetic insulator at a density of one electron per site. We argue that topological superconductivity emerges upon doping in the vicinity of the integer quantum Hall to chiral spin liquid transition. We employ exact diagonalization and density matrix renormalization group methods to examine this theoretical scenario and find that electronic pairing indeed occurs on both sides of criticality over a remarkably broad range of interaction strengths. On the chiral spin liquid side, our results provide a concrete model realization of the long-hypothesized mechanism of anyon superconductivity. Our study thus establishes a beyond-Bardeen-Cooper-Schrieffer route to electron pairing in a well-controlled limit, relying crucially on the interplay between electron correlations and band topology.

Hubbard model↗
DOE OSTI2020

Wireless Interrogation of High Temperature Surface Acoustic Wave Dynamic Strain Sensor

Dynamic strain sensing is necessary for high-temperature harsh-environment applications, including powerplants, oil wells, aerospace, and metal manufacturing. Monitoring dynamic strain is important for structural health monitoring and condition-based maintenance in order to guarantee safety, increase process efficiency, and reduce operation and maintenance costs. Sensing in high-temperature (HT), harsh-environments (HE) comes with challenges including mounting and packaging, sensor stability, and data acquisition and processing. Wireless sensor operation at HT is desirable because it reduces the complexity of the sensor connection, increases reliability, and reduces costs. Surface acoustic wave resonators (SAWRs) are compact, can operate wirelessly and battery-free, and have been shown to operate above 1000°C, making them a potential option for HT HE dynamic strain sensing. This paper presents wirelessly interrogated SAWR dynamic strain sensors operating around 288.8MHz at room temperature and tested up to 400°C. The SAWRs were calibrated with a high-temperature wired commercial strain gauge. The sensors were mounted onto a tapered-type Inconel constant stress beam and the assembly was tested inside a box furnace. The SAWR sensitivity to dynamic strain excitation at 25°C, 100°C, and 400°C was .439 μV/με, 0.363μV/με, and .136 μV/με, respectively. Furthermore, the experimental outcomes verified that inductive coupled wirelessly interrogated SAWRs can be successfully used for dynamic strain sensing up to 400°C.

47 OTHER INSTRUMENTATION↗
NASA NTRS1984

Effects of dynamical heat fluxes on model climate sensitivity

A coupled high and low latitude radiative-dynamical model of the annual mean northern hemisphere has been constructed in order to study the interactions of the vertical and meridional heat fluxes and their feedback effect on model climate sensitivity. The model's climate sensitivity to solar constant changes and CO2 increases is investigated, and the effect of feedback in the dynamical fluxes on model climate sensitivity is examined. Nonlinear interactions between heat fluxes and other feedbacks such as radiation-temperature, ice albedo, and humidity are also discussed.

Wang, W.-C.↗
NASA NTRS1987

Biot theory and acoustical properties of high porosity fibrous materials and plastic foams

Experimental values of acoustic wave propagation constant and characteristic impedance in fibrous materials, and normal absorption for two plastic foams, were compared to theoretical predictions obtained with Biot's theory. The best agreement was observed for fibrous materials between Biot's theory and Delany and Bazley experiments for a nearly zero mass coupling parameter. For foams, the lambda/4 structure resonance effect on absorption was calculated by using four-pole modelling of the medium. A significant mass coupling parameter is then necessary for obtaining agreement between the behavior of the measured absorption coefficients and the theoretical predictions. It is shown how the formalism used for predicting foams absorption coefficients may be used for studying the acoustic behavior of multi-layered media.

Allard, J.↗
NASA NTRS1987

Comparison of Hirs' equation of Moody's equation for determining rotordynamic coefficients of annular pressure seals

The rotordynamic coefficients of an incompressible-flow annular pressure seal were determined using a bulk-flow model in conjunction with two different friction factor relationships. The first, Hirs' equation, assumes the friction factor is a function of Reynolds number only. The second, Moody's equation, approximates Moody's diagram and assumes the friction factor is a function of both Reynolds number and relative roughness. For each value of relative roughness, Hirs' constants were determined so that both equations gave the same magnitude and slope of the friction factor. For smooth seals, both relationships give the same results. For rough seals (e/2 H sub 0 = 0.05) Moody's equation predicts 44% greater direct stiffness, 35% greater cross-coupled stiffness, 19% smaller cross-coupled damping, 59% smaller cross-coupled inertia, and nominally the same direct damping and direct inertia.

Nelson, Clayton C.↗
NASA NTRS1997

[Photodissociation of Acetylene and Acetone using Step-Scan Time-Resolved FTIR Emission Spectroscopy]

The photodissociation of acetylene and acetone was investigated as a function of added quenching gas pressures using step-scan time-resolved FTIR emission spectroscopy. Its main components consist of Bruker IFS88, step-scan Fourier Transform Infrared (FTIR) spectrometer coupled to a flow cell equipped with Welsh collection optics. Vibrationally excited C2H radicals were produced from the photodissociation of acetylene in the unfocused experiments. The infrared (IR) emission from these excited C2H radicals was investigated as a function of added argon pressure. Argon quenching rate constants for all C2H emission bands are of the order of 10(exp -13)cc/molecule.sec. Quenching of these radicals by acetylene is efficient, with a rate constant in the range of 10(exp -11) cc/molecule.sec. The relative intensity of the different C2H emission bands did not change with the increasing argon or acetylene pressure. However, the overall IR emission intensity decreased, for example, by more than 50% when the argon partial pressure was raised from 0.2 to 2 Torr at fixed precursor pressure of 160mTorr. These observations provide evidence for the formation of a metastable C2H2 species, which are collisionally quenched by argon or acetylene. Problems encountered in the course of the experimental work are also described.

McLaren, Ian A.↗
DOE OSTI2020

Consequence Analysis of Residual Water in a Storage Canister - 20484

Recent observations from projects sponsored under the DOE-Nuclear Energy Spent Fuel and Waste Disposition campaign show that residual free water, well above the 0.4 gm-moles that correspond to a 3 torr rebound pressure after a 30-minute hold following prototypic drying, may remain within SNF canisters used for long-term dry storage. Specifically, experimental observations in an Integrated Research Project showed subzero temperatures and a local pocket of frozen water in a mockup fuel assembly subjected to an industry practiced drying process. In gas sampling and analysis from the High Burnup (HBU) Demonstration project, using a fully loaded and dried SNF canister, indicated presence of 100 g of residual water had remained in the canister. The effects of residual water in terms of potential adverse impacts on the condition of the SNF and the canister internals during long-term dry storage has been rigorously evaluated to address inadvertent residual water contents. The impact of residual water on the fuel and canister internals, with a focus on the oxidation of the cladding and of postulated exposed (breached-clad) fuel, were evaluated using an integrated set of time-dependent environmental conditions with semi-empirical materials oxidation models. An 'integration model,' with coupled sub-models, developed and previously used by the Center for Nuclear Waste Regulatory Analysis (CNWRA), was applied and expanded with improved sub-models for water radiolysis, and for cladding and fuel oxidations. Cladding oxidation models for various alloys were developed using the rate constants for various cladding alloys in the dry storage temperature range. A burnup-dependent fuel oxidation model was developed; it accounted for inhibition effect of fission products on fuel oxidation rates. The potential for extension of the initial breach flaw (1 mm hole) was evaluated by two separate mechanics criteria. In the first case, a critical cladding strain criterion was used. In the second case, a critical stress intensity was used to evaluate stability of a planar axial flaw postulate. A failure condition exists when a sufficient amount of fuel oxidation (at the postulated breached-clad location) causes fuel pellet swelling to load the cladding to meet a mechanics criterion for flaw extension. The integration model's results provide the following insights. The radiolysis-generated oxygen is consumed with a partitioning between the cladding and the fuel. The extent of cladding oxidation is no more than 2 μm additional consumption of the cladding metal even with 10 moles of residual water, indicating that changes in cladding conditions due to water are expected to be negligible. The conditions for cladding failure (breach opening size extension) directly correlates with the extent of fuel oxidation. Fuel oxidation significant to extend the breach could occur even with 5.5 moles of residual water. The results indicate that fuel oxidation vis-a-vis cladding oxidation is more likely under a higher radiation field, i.e., when radiolytic decomposition of the residual water occurs in few years compared to several decades, and that the conditions to cause failure are increased with fuel and cladding temperatures approaching the peak storage temperature limit of 400 deg. C. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Thermodynamically consistent algorithms for models of incompressible multiphase polymer solutions with a variable mobility

Here we present a general strategy for developing structure and property preserving numerical algorithms for thermodynamically consistent models of incompressible multiphase polymer solutions with a variable mobility. We first present a formalism to derive thermodynamically consistent, incompressible, multiphase polymer models. Then, we develop the general strategy, known as the supplementary variable method, to devise thermodynamically consistent numerical approximations to the models. We illustrate the numerical strategy using newly developed models of incompressible diblock copolymer solutions coupled with an electric and a magnetic field, respectively. Mesh refinement is conducted to verify convergence rates of the developed schemes. Some numerical examples are given to exhibit underlying dynamics absent from and driven by the external fields, respectively, highlighting differences between models with the variable and constant mobilities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Pore connectivity influences mass transport in natural rocks: Pore structure, gas diffusion and batch sorption studies

For this work, six rocks (one granodiorite, one limestone, two chalks, one mudstone, and one dolostone) with different extents of heterogeneity at six different particle sizes (from 75 to 8000 μm) were studied to describe the effects of pore connectivity on mass transport. The methods applied were (i) porosity measurement of granular rocks, (ii) analyses of gas-phase diffusive transport in a bed of packed particles, along with a solid quartz method at these six particle sizes being developed to identify the contribution of intraparticle diffusion, and (iii) batch sorption tests of multiple ions (anions and cations) with subsequent analyses of inductively coupled plasma-mass spectrometry. Granular porosity measurement results reveal that with decreasing particle sizes, the effective porosities for the “heterogenous” group of rocks (Grimsel granodiorite and Edwards limestone) increase, whereas the porosities of another “homogeneous” group (two Israel chalk samples, Japan mudstone, and Wyoming dolostone) remain constant. Gas diffusion results show that the intraparticle gas diffusion coefficient among these two sample groups, varying in the magnitude of 10 -8 to 10 -6 m 2 /s, are not directly correlated to the porosity differences. Moreover, the batch sorption work displays a different affinity of rocks for various tracers. For Grimsel granodiorite, Japan mudstone, and Wyoming dolostone, the adsorption capacity of Sm 3+ and Eu 3+ increases as the particle size decreases. In general, this integrated research of grain size distribution, granular rock porosity, intraparticle diffusivity, and ionic sorption capacity gives insights into the pore connectivity effect on both physical and chemical transport behaviors for different lithologies and/or different particle sizes.

58 GEOSCIENCES↗