Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “point contacts”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A direct contact cooler design for simultaneously recovering latent heat and capturing SO x and NO x from pressurized flue gas

Integrated SO x and NO x removal technology using a direct contact cooler is a promising alternative for cleaning coal combustion flue gas, especially for pressurized combustion systems. Past investigations of integrated SO x and NO x removal have been limited to low-temperature processes, where the flue gas has been cooled upstream of the removal unit, and the moisture in the flue gas has been condensed out. Here this work presents a novel design that recovers the flue gas heat while simultaneously removing SO x and NO x in a single reactive-absorption column – direct contact cooler. A validated model was used to evaluate different direct contact cooler designs. Modeling results suggested that a counter-current column with a single water inlet can scrub 83% of the NO x and 96.5% of the SO x in 115 s, while an optimized design with multiple water inlets enhances the scrubbing of NO x by 9 percentage points and SO x by 3.5 percentage points, due to a 34% increase in residence time. It is also observed that the liquid-phase reaction between absorbed SO x and NO x plays a significant role in the removal of SO 2 from the flue gas for high-temperature scrubbing compared with low-temperature scrubbing. The design reinforces the potential of high-temperature SO x -NO x removal and heat recovery and suggests a means of reducing capital costs for such columns.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Embedding space approach to Lorentzian CFT amplitudes and causal spherical functions

Conformal field theory in a Minkowski setting is discussed in an embedding space approach, paying special attention to causality constraints for four-point amplitudes. The physics of dilatation and Lorentz boost is emphasized in specifying the noncompact maximal Abelian subgroup of S O ( d , 2 ) . Reduction of a conformal field theory four-point amplitudes as functions of cross ratios is shown to be equivalent to enforcing H bi-invariance, i.e., F ( h g h ′ ) = F ( g ) , with g ∈ S O ( d , 2 ) and H an appropriate subgroup. Causality is imposed by introducing appropriate semigroups. Causal zonal spherical functions are constructed, making contact with Minkowski conformal blocks introduced previously. Published by the American Physical Society 2024

Agarwal, Pulkit (ORCID:0000000346581691)↗
DOE OSTI2025

Homotopy Solver

This software implements parallel versions of an interior-point solver, based on the publicly available ipopt solver. Here we have full control over the linear solver and our algorithm is fully parallel thus enabling scalability to large-scale optimization problems. This package also has a parallel implementation of a homotopy solver developed under the scalable methods for contact LDRD project 23-ERD-017. This solver is an mfem-based implementation of algorithm described in ``A filter trust-region Newton continuation method for nonlinear complementarity problems''. Cosmin G. Petra, Nai-Yuan Chiang, Jingyi Wang, Tucker Hartland, and Michael Puso (submitted), LLNL-JRNL-869761.

Hartland, Tucker [Lawrence Livermore National Labo↗
DOE OSTI2021

Microstructural Evaluation of Phase Instability in Large Bandgap Metal Halide Perovskites

The optoelectronic performance of organic–inorganic halide perovskite (OIHP)-based devices has been improved in recent years. Particularly, solar cells fabricated using mixed-cations and mixed-halides have outperformed their single-cation and single-halide counterparts. Yet, a systematic evaluation of the microstructural behavior of mixed perovskites is missing despite their known composition-dependent photoinstability. Here, we explore microstructural inhomogeneity in (FAPbI 3 ) x (MAPbBr 3 ) 1– x using advanced scanning probe microscopy techniques. In this work, contact potential difference (CPD) maps measured by Kelvin probe force microscopy show an increased fraction of grains exhibiting a low CPD with flat topography as MAPbBr 3 concentration is increased. The higher portion of low CPD contributes to asymmetric CPD distribution curves. Chemical analysis reveals these grains being rich in MA, Pb, and I. The composition-dependent phase segregation upon illumination, reflected on the emergence of a low-energy peak emission in the original photoluminescence spectra, arises from the formation of such grains with flat topology. Bias-dependent piezo-response force microscopy measurements, in these grains, further confirm vigorous ion migration and cause a hysteretic piezo-response. Our results, therefore, provide insights into the microstructural evaluation of phase segregation and ion migration in OIHPs pointing toward process optimization as a mean to further enhance their optoelectronic performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

A Simplified Approach Based on Cellular Automata for Describing Direct Reduced Iron Production in Different Reducing Conditions

A quick computation approach based on cellular automata is developed and implemented to describe the reduction of iron ore pellets by a mixture of reducing agents featured by different H 2 /CO ratios. The evolution of oxygen concentration inside the pellet is followed from the beginning to the end of contact between reducing agent and pellet. The variation of thermal state of pellets and gas mixture is computed based on their initial temperature, considering the heat involved and the convective heat exchange between pellet and gas mixture. The use of cellular automata and finite‐difference method to solve the diffusion equation point out the absence of any diffusion coefficient value, allowing to make the model fit the experimental trial, because the problem is that it is not ruled just by diffusion but also by the concentration variation of reducing agent inside the pellet due to porosity increasing during reduction. The updating of the reducing agents concentration implies a sharp decrease of oxygen concentration that the cellular automata model considers. The developed model is able to provide the in‐line control of reduction process and could be used to adjust the chemical concentration and temperature of injected reducing agents.

Metallurgy & Metallurgical Engineering↗
DOE OSTI2023

Spatial ecology of translocated raccoons

Raccoons (Procyon lotor) are routinely translocated both legally and illegally to mitigate conflicts with humans, which has contributed to the spread of rabies virus across eastern North America. The movement behavior of translocated raccoons has important ramifications for disease transmission yet remains understudied and poorly quantified. To examine the spatial ecology of raccoons following experimental translocation, we performed reciprocal 16 km-distance translocations of 30 raccoons between habitats of high and low raccoon density (bottomland hardwood and upland pine, respectively) across the Savannah River Site (SRS) in Aiken, South Carolina, USA (2018–2019). Translocation influenced patterns of raccoon space use, with translocated animals exhibiting a 13-fold increase in 95% utilization distributions (UDs) post- compared to pre-translocation (mean 95% UD 35.8 ± 36.1 km 2 vs 1.96 ± 1.17 km 2 ). Raccoons originating from upland pine habitats consistently had greater space use and larger nightly movement distances post-translocation compared to raccoons moved from bottomland hardwood habitats, whereas these differences were generally not observed prior to translocation. Estimated home ranges of male raccoons were twice the area as estimated for female raccoons, on average, and this pattern was not affected by translocation. After a transient period lasting on average 36.5 days (SD = 30.0, range = 3.25–92.8), raccoons often resumed pre-experiment movement behavior, with 95% UD sizes not different from those prior to translocation (mean = 2.27 ± 1.63km 2 ). Most animals established new home ranges after translocation, whereas three raccoons moved > 16 km from their release point back to the original capture location. Four animals crossed a 100-m wide river within the SRS post-translocation, but this behavior was not documented among collared raccoons prior to translocation. Large increases in space use combined with the crossing of geographic barriers such as rivers may lead to elevated contact rates with conspecifics, which can heighten disease transmission risks following translocation. These results provide additional insights regarding the potential impacts of raccoon translocation towards population level risks of rabies outbreaks and underscore the need to discourage mesocarnivore translocations to prevent further spread of wildlife rabies.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

The latent variable proximal point algorithm for variational problems with inequality constraints

The latent variable proximal point (LVPP) algorithm is a framework for solving infinite-dimensional variational problems with pointwise inequality constraints. The algorithm is a saddle point reformulation of the Bregman proximal point algorithm. At the continuous level, the two formulations are equivalent, but the saddle point formulation is more amenable to discretization because it introduces a structure-preserving transformation between a latent function space and the feasible set. Working in this latent space is much more convenient for enforcing inequality constraints than the feasible set, as discretizations can employ general linear combinations of suitable basis functions, and nonlinear solvers can involve general additive updates. LVPP yields numerical methods with observed mesh-independence for obstacle problems, contact, fracture, plasticity, and others besides; in many cases, for the first time. The framework also extends to more complex constraints, providing means to enforce convexity in the Monge–Ampère equation and handling quasi-variational inequalities, where the underlying constraint depends implicitly on the unknown solution. Here, in this paper, we describe the LVPP algorithm in a general form and apply it to ten problems from across mathematics.

Inequality constraints↗
DOE OSTI2026

Nanoscale control over single vortex motion in an unconventional superconductor

Precise control of superconducting vortices is crucial for studying vortex dynamics and vortex braiding. We propose a new method to pull vortex lines in the layered superconductor FeSe using a scanning tunneling microscope (STM) tip. Weak contact of the STM tip with the FeSe surface locally reduces the superconducting gap, thereby creating a tunable vortex pinning potential on the nanometer scale. This enables controlled vortex line deformation even in dense vortex lattices. Analytical modeling reveals that the deformation strength scales logarithmically with conductance and depends on tip geometry. Our findings point to local strain-induced gap suppression as the mechanism for STM-mediated vortex manipulation, providing a fundamental insight relevant to vortex behavior in the context of quantum information studies.

Song, Sang Yong [Oak Ridge National Laboratory (OR↗
DOE OSTI2023

New alignment method for remote protein sequences by the direct use of pairwise sequence correlations and substitutions

Understanding protein sequences and how they relate to the functions of proteins is extremely important. One of the most basic operations in bioinformatics is sequence alignment and usually the first things learned from these are which positions are the most conserved and often these are critical parts of the structure, such as enzyme active site residues. In addition, the contact pairs in a protein usually correspond closely to the correlations between residue positions in the multiple sequence alignment, and these usually change in a systematic and coordinated way, if one position changes then the other member of the pair also changes to compensate. In the present work, these correlated pairs are taken as anchor points for a new type of sequence alignment. The main advantage of the method here is its combining the remote homolog detection from our method PROST with pairwise sequence substitutions in the rigorous method from Kleinjung et al. We show a few examples of some resulting sequence alignments, and how they can lead to improvements in alignments for function, even for a disordered protein.

Jia, Kejue↗
DOE OSTI2023

Co-deposition of hole-selective contact and absorber for improving the processability of perovskite solar cells

Simplifying the manufacturing processes of renewable energy technologies is crucial to lowering the barriers to commercialization. In this context, to improve the manufacturability of perovskite solar cells (PSCs), here we have developed a one-step solution-coating procedure in which the hole-selective contact and perovskite light absorber spontaneously form, resulting in efficient inverted PSCs. We observed that phosphonic or carboxylic acids, incorporated into perovskite precursor solutions, self-assemble on the indium tin oxide substrate during perovskite film processing. They form a robust self-assembled monolayer as an excellent hole-selective contact while the perovskite crystallizes. Our approach solves wettability issues and simplifies device fabrication, advancing the manufacturability of PSCs. Our PSC devices with positive-intrinsic-negative (p-i-n) geometry show a power conversion efficiency of 24.5% and retain >90% of their initial efficiency after 1,200 h of operating at the maximum power point under continuous illumination. The approach shows good generality as it is compatible with different self-assembled monolayer molecular systems, perovskites, solvents and processing methods.

14 SOLAR ENERGY↗
DOE OSTI2020

Conductive Thin Films over Large Areas by Supramolecular Self-Assembly

We report a “one-step” method for preparing conductive thin films with cylindrical microdomains oriented normal to the surface over large areas using the supramolecular assembly of poly(styrene-block-4-vinylpyridine) (PS19-b-P4VP 5 ) and 5,10,15,20-tetrakis(4-hydroxyphenyl)-21H,23H-porphine (HOTPP). HOTPP interacts with the P4VP block by hydrogen bonding between the hydroxyl group of HOTPP and pyridine ring of PS 19 -b-P4VP 5 , forming cylindrical P4VP(HOTPP) domains having an average diameter of ~17 nm in a PS matrix. Here, dynamic light scattering, contact angle, and in situ grazing incidence small-angle X-ray scattering measurements show a morphological transition from spherical micelles in solution to cylindrical microdomains oriented normal to the substrate surface during the drying process. From the dependence of current on voltage, an average current of ~4.0 nA is found to pass through a single microdomain, pointing to a promising route for organic semiconductor device applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Assessing Critical Conditions for Scour Near Obstructions using Bed Shear, Particle Onset of Motion Balances, and CFD-DEM Modeling of Granular Beds

Computational Fluid Dynamics combined with a Discrete Element Method is one of the computational methods that can be used to model multiphase flows. In this method various phases, gas and liquid or solid, are present in the same computational domain. The local averaged Navier–Stokes equations determine the flow of the continuous phase fluid and are solved using the traditional CFD finite volume approach. DEM is based on a Lagrangian formulation, which solves the equations of motion, expressed in ordinary differential equations, for representative particles as they move in space and time. The interactions between the continuous fluid phase and discrete solid phase are modeled with the use of Newton’s laws of motion via drag force. The particles interact with each other and with the boundaries of the fluid continuum, and the resulting contact forces are included in the equations of motion. The properties of solid particles and boundaries are treated as elastic bodies, with specified density, elastic modulus, and Poisson’s ratio. Particle shapes may vary from single spherical particles to more complex-shaped composite particles. The particles may be introduced into the domain by random or structured injection at a point, surface, or volume, depending on the application. More details on the formulation can be found in the Simcenter STAR-CCM+ User’s Manual and OpenFOAM website.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

BCJ amplitude relations for Anti-de Sitter boundary correlators in embedding space

We generalize the color/kinematics duality of flat-space scattering amplitudes to the embedding space formulation of AdS boundary correlators. Kinematic numerators and propagators are replaced with differential operators acting on a scalar contact diagram that is the AdS generalization of the momentum conserving delta function of flat space scattering amplitudes. We show that color/kinematics duality implies differential relations among AdS boundary correlators that naturally generalize the flat space BCJ amplitude relations and verify them for the correlators of Yang-Mills theory and of the Nonlinear Sigma Model through four- and six-points, respectively. For the latter we also find representations of the four- and six-point correlator that manifest the duality. Possible double-copy procedures in AdS space are also discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Unusual Li 2 O sublimation promotes single-crystal growth and sintering

Li 2 O is rarely used for cathode material synthesis due to its high melting point (1,438 °C). Here we discover that Li 2 O can sublimate at 800-1,000 °C under ambient pressure, opening new possibilities for cathode synthesis. We propose a mechanism that enables synthesis of single crystals-such as LiNi 0.8 Mn 0.1 Co 0.1 O 2 (NMC811) or LiNi 0.9 Mn 0.05 Co 0.05 O 2 (NMC90)-without direct contact with Li 2 O salts. We show that Li 2 O vapour successfully converts spent polycrystalline NMC811 into segregated single crystals without milling or post-treatment. The Li 2 O vapour, derived from Li 2 O solids, diffuses rapidly and reacts with precursors, mimicking a molten-salt environment, which facilitates single-crystal growth. The chemical lithiation process continuously drives Li 2 O sublimation, sintering the crystals. Single crystals derived from Li 2 O and fresh precursors or spent polycrystals exhibit outstanding cycling after 1,000 cycles in full cells. The demonstrated Li 2 O sublimation and its universal role in promoting single-crystal growth provides an effective approach for single-crystal synthesis, scale-up and recycling.

25 ENERGY STORAGE↗
DOE OSTI2020

CRNT4SBML: a Python package for the detection of bistability in biochemical reaction networks

Motivation: Signaling pathways capable of switching between two states are ubiquitous within living organisms. They provide the cells with the means to produce reversible or irreversible decisions. Switchlike behavior of biological systems is realized through biochemical reaction networks capable of having two or more distinct steady states which are dependent on initial conditions. Investigation of whether a certain signaling pathway can confer bistability involves a substantial amount of hypothesis testing. The cost of direct experimental testing can be prohibitive. Therefore, constraining the hypothesis space is highly bene?cial. One such methodology is based on Chemical Reaction Network Theory which uses computational techniques to rule out pathways that are not capable of bistability regardless of kinetic constant values and molecule concentrations. Although useful, these methods are complicated from both pureandcomputationalmathematicsperspectives.Thus,theiradoptionisverylimitedamongstbiologists. Results: We brought Chemical Reaction Network Theory approaches closer to experimental biologists by automating all the necessary steps in CRNT4SMBL. The input is based on SBML format which is the community standard for biological pathway communication. The tool parses SBML and derives C-graph representations of the biological pathway with mass action kinetics. Next steps involve an ef?cient search for potential saddle-node bifurcation points using an optimization technique. This type of bifurcation is important as it has the potential of acting as a switching point between two steady states. Finally, if any bifurcation points are present, numerical continuation analysis extends the equilibria branches for generating the diagram. Presence of an S-shaped bifurcation diagram indicates that the pathway acts as a bistable switch for the given optimization parameters. Availability: CRNT4SBML is available via the Python Package Index. The documentation can be found at https://crnt4sbml.readthedocs.io. CRNT4SBML is licensed under the Apache Software License 2.0. Contact: vladislav.petyuk@pnnl.gov

Reyes, Brandon C.↗
DOE OSTI2023

CCN data from TAMU TRACER campaign in the Houston TX region from July to September 2022

During TRACER, the Texas A&M Rapid Onsite Atmospheric Measurements Van (ROAM-V) was deployed to capture airmasses behind (maritime) and ahead (continental) of the passage of the sea-breeze front through Houston. On select sampling days, ROAM-V sampled in the morning/mid-day on the coast and then transited to a second inland site for the afternoon/evening. The suite of instruments deployed on ROAM-V included a Condensation Particle Counter (CPC; GRIMM Model 5.403 CPC), Scanning Mobility Particle Sizer (SMPS; TSI 3750 detector, TSI 3082 classifier, TSI 3088 neutralizer, TSI 3081A Differential Mobility Analyzer), Cloud Condensation Nuclei counter (Droplet Measurement Technologies CCN Counter), micro pulse lidar (Droplet Measurement Technologies Micro Pulse LiDAR (miniMPL)), and a Davis Rotating Uniform size-cut Monitor (DRUM; DRUMAir 4-DRUM). Before sampling at each location, the latitude and longitude were recorded using the GPS on the phone application “My Altitude”. Onboard the ROAM-V, aerosol samples are drawn through a shared isokinetic inlet at a flow rate ranging from 3.5 to 7.0 LPM. A portion of this flow is directed through a cyclone impactor (Brechtel, Inc. Model SCC 0.732) and 0.5 LPM is directed to the CCN. To calculate particle losses, we used a two-step method. First, the measured SMPS size distributions were used to calculate particle loss through the inlet during sampling. Second, the corrected SMPS data was used to calculate the average of the total losses per scan down the CCN line. Then, the correction was applied to the CCN data. This calculation was done separately for each deployment location due to changes in the measured size distributions between locations. Particle loss from diffusion (based on Kesten, 1991 and Gormley, 1949), inertial impaction in 90-degree bends (based on Aerosol Measurement, 2011 and Crane, 1977), and cyclone impactor efficiency (based on Dirgo, 1985) were included in the loss calculation. When the SMPS was not sampling at a location (in the case of an instrument malfunction or operator error), the reported CPC data was corrected with an average of the total losses for the entire campaign at the specified deployment location (e.g., if we needed to correct Galveston data, then the average of all calculated losses at Galveston was taken). These flatline corrections were used for all data on 22/07/13, 22/07/20, 22/07/22, and the data from Galveston on 22/08/09. The supersaturation uncertainty is estimated conservatively at +/- 0.03%, where variation in the inlet temperature, pressure, and calibration technique prevents a more accurate measurement. Confidence in the reported supersaturation measurements is based on a pre-campaign calibration (following the methods from our previous work and Deng, 2014 based on Rose, 2008) in addition to inter-comparisons with the DOE for two days (22/08/18 and 22/09/01) where TAMU was co-located with AMF1. The inter-comparisons show good agreement between our instrument and the DOEs instrument on both days at all supersaturations. After the last inter-comparison on 22/09/01, there was no indication of a malfunction by our instrument through the rest of the campaign. Unfortunately, the instrument was dropped during demobilization. A post-campaign calibration was conducted, which showed a substantial departure from the pre-campaign calibration. The drop may have damaged the instrument’s ability to produce the desired supersaturations. Therefore, we do consider the data after 22/09/01 to be correct, but it should be used with caution. The CCN counter sampled for 3 minutes at each supersaturation setpoint (0.2, 0.4, 0.6, 0.8, 1.0, and 1.2%). At the end of a cycle, the instrument was set to 0.01% supersaturation for 5 minutes. The data is comprised of the last 60 seconds of each supersaturation set point (0.2, 0.4, 0.6, 0.8, 1.0, and 1.2%) to ensure the instrument stabilized and was able to reach thermal equilibrium. We removed the data during the periods where there were operational difficulties, setup, or maintenance. This data was collected for ARM Field Campaign AFC07055 and supported by DOE ASR grant DE-SC0021047. For any further questions, please feel free to contact the instrument PI, Sarah D. Brooks, sbrooks@tamu.edu. Rose et. al. Calibration and Measurement Uncertainties of a Continuous-Flow Cloud Condensation Nuclei Counter (DMT-CCNC): CCN Activation of Ammonium Sulfate and Sodium Chloride Aerosol Particles in Theory and Experiment. Atmos. Chem. Phys., 8, 1153-1179, 2008. Deng et. al. Using Raman Microspectroscopy to Determine Chemical Composition and Mixing State of Airborne Marine Aerosols over the Pacific Ocean. Aerosol Science and Technology, Vol 48, Issue 2, 2014. Kesten et. al. Calibration of a TSI Model 3025 Ultrafine Condensation Particle Counter. Aerosol Science and Technology, 15:2, 107-111, 1991. Gormley et. al. Diffusion from a Stream Flowing through a Cylindrical Tube. Proceedings of the Royal Irish Academy, Vol 52, 163-169, 1948. Aerosol Measurement: Principles, Techniques, and Applications, Third Edition. John Wiley & Sons, Inc, 2011. Crane et. al. Inertial Deposition of Particles in a Bent Pipe. Journal of Aerosol Science, Vol 8, 161-170, 1977. Dirgo et. al. Cyclone Collection Efficiency: Comparison of Experimental Results with Theoretical Predictions. Aerosol Science and Technology, 4:4, 401-415, 1985.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Gamma Spectrometry Code Rodeo for Uranium Enrichment—FY 2021 Report

In FY22, the participating Department of Energy (DOE) laboratories continued the work on the project “Gamma Spectrometry Code Rodeo for Uranium Enrichment”. Oak Ridge National Laboratory (ORNL) and Lawrence Livermore National Laboratory (LLNL) continued their efforts to acquire spectra using uranium sources. A collimator (with side and back shields) was custom designed and constructed at ORNL. Spectra were acquired at ORNL using the uranium enrichment standards supplied by New Brunswick Laboratory (NBL) and the collimated M400 detector. Spectra were also acquired using steel absorbers of different thickness placed between the source and the collimated M400 detector. The plan was to ship the collimator to LLNL so that spectra could be acquired using the LLNL detector in collimated geometries. This activity was suspended since problems were noticed with the performance of the LLNL M400 detector. Also, it was brought to light to the project team that some of the aspects of the M400 detector had been updated by the vendor H3D in response to requests by the IAEA. In February 2021, the ORNL and LLNL M400 detectors were returned to the H3D factory in order to resolve the problems with the LLNL detector and to update ORNL and LLNL detectors to the same configuration as the M400 supplied to the IAEA. H3D repaired the LLNL M400 and performed the IAEA updates on both the ORNL and LLNL detectors. ORNL and LLNL received the detectors in April 2021. Data acquisition was continued using the upgraded detectors. ORNL re-collected the spectra using point sources and the NBL standards in uncollimated and collimated geometries. The collimator was shipped from ORNL to LLNL. Spectra were acquired at LLNL through Q4 of FY22. When the ORNL spectra were examined closely, it became apparent that there were spurious artifacts present in some of the spectra. This rendered suspect the ORNL spectra collected in June/July 2021 time frame. ORNL contacted H3D and per H3D’s advice, the M400 detector was returned to H3D for repairs in mid-August 2021. H3D diagnosed the problem and established that the root cause was the extra tight packaging that created a strain on the CZT crystals, causing one of them to break down. The repaired M400 was returned to ORNL in early September 2021. ORNL re-started data acquisition for the third time. Data was continued to be acquired at ORNL and LLNL through September 2021. The ORNL and LLNL data will be shared with the principal investigators from the analyzing laboratories. The GADRAS code and FRAM have been developed to have the capability to analyze the M400 spectra. Good quality spectra from FY21 have been analyzed using modified versions of GADRAS and FRAM. The spectra have also been analyzed using the GEM code. The analysis results from GADRAS, FRAM and GEM codes are presented and discussed in this report. Some preliminary results from FRAM and GEM analysis were presented in February 2021 at the U-Pu Isotopics Workshop sponsored by the IAEA. CZTU code is being modified to analyze M400 spectra. Data collection, code development and analysis will be continued in FY22 as per the FY22 Project Work Plan (PWP).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Fluid – Fluid Interfacial Area and Its Impact on Relative Permeability - A Pore Network Modeling Study

Abstract Relative permeability (kr) is commonly modeled as an empirical function of phase saturation. Although current empirical models can provide a good match of one or two measured relative permeabilities using saturation alone, they are unable to predict relative permeabilities well when there is hysteresis or when physical properties such as wettability change. Further, current models often result in relative permeability discontinuities that can cause convergence and accuracy problems in simulation. To overcome these problems, recent research has modeled relative permeability as a state function of both saturation (S) and phase connectivity (X). Pore network modeling (PNM) data, however, shows small differences in relative permeability for the same S-X value when approached from a different flow direction. This paper examines the impact of one additional Minkowski parameter (Mecke and Arns, 2005), the fluid-fluid interfacial area, on relative permeability to identify if that satisfactorily explains this discrepancy. We calculate the total fluid-fluid interfacial areas (IA) during two-phase (oil/water) flow in porous media using pore network modeling. The area is calculated from PNM simulations using the areas associated with corners and throats in pore elements of different shapes. The pore network is modeled after a Bentheimer sandstone, using square, triangular, and circular pore shapes. Simulations were conducted for numerous primary drainage and imbibition cycles at a constant contact angle of 0° for the wetting phase. Simultaneous measurements of capillary pressure, relative permeability, saturation, and phase connectivity are made for each displacement. Fluid-fluid interfacial area is calculated from the PNM capillary pressure, the fluid location in the pore elements, and the pore element dimensional data. The results show that differences in the relative permeability at the same (S,X) point is explained well by differences in the fluid-fluid interfacial area (IA). That is, for a larger change in IA at these intersection points, the permeability difference is greater. That difference in relative permeability approaches zero as the difference in IA approaches zero. This confirms that relative permeability can be modeled better as a unique function of S, X, and IA. The results also show that an increase in IA restricts flow decreasing the nonwetting (oil) phase permeability. This decrease is caused by an increase in the throat area fraction compared to the corner area as the total area IA increases. The wetting phase relative permeability, however, shows the inverse trend, in that its relative permeability is greater when IA becomes larger owing to a greater fraction of the total area associated with the corners. The area IA, however, impacts the nonwetting phase relative permeability more than the wetting phase relative permeability. Corner flow improves the wetting phase relative permeability because the wetting phase is continuous there. Finally, a sensitivity analysis shows that relative permeability a is more sensitive to change in S than they are for IA for the case studied implying that if only two parameters are used to model relative permeability it is better to choose S and X.

Engineering↗