Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed computation”

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,423 records · Page 79
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Robust Importance Sampling for Bayesian Model Calibration with Spatio-Temporal Data

This paper addresses two challenges in Bayesian calibration: 1) computational speed of existing sampling algorithms, and 2) calibration with spatio-temporal responses. The commonly used Markov Chain Monte Carlo (MCMC) approaches require many sequential model evaluations making the computational expense prohibitive. This paper proposes an efficient sampling algorithm: iterative importance sampling with genetic algorithm (IISGA). While iterative importance sampling enables computational efficiency, the genetic algorithm enables robustness by preventing sample degeneration and avoids getting stuck in multimodal search spaces. An inflated likelihood further enables robustness in high-dimensional parameter spaces by enlarging the target distribution. Spatio-temporal data complicate both surrogate modeling, which is necessary for expensive computational models, and the likelihood estimation. In this work, singular value decomposition is investigated for reducing the high-dimensional field data to a lower-dimensional space prior to Bayesian calibration. Then the likelihood is formulated and Bayesian inference is performed in the lower-dimension, latent space. An illustrative example is provided to demonstrate IISGA relative to existing sampling methods, and then IISGA is employed to calibrate a thermal battery model with 26 uncertain calibration parameters and spatio-temporal response data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

ADVANCED THERMAL-HYDRAULIC MODEL OF HEAT RECOVERY STEAM GENERATORS

Steam-side oxidation at elevated temperatures in heat recovery steam generators (HRSGs) can lead to overheating of local tubes, potential oxide scale exfoliation, and damage to critical equipment. Predicting the oxide growth distribution requires detailed knowledge of the distribution of tube metal temperatures. To address these challenges, a computational fluid dynamics (CFD) model of a high-pressure superheater (HPSH) was developed within the Siemens’ STAR-CCM+ CFD commercial software. The HPSH section of the steam generator was selected for modeling because of its high oxidation rates.The model development included several steps. In the first step, an oxidation model based on a parabolic oxidation relationship was implemented. The oxidation growth predicted by the model was then verified for a full-load boiler operation case. In the second step, a three-dimensional CFD model was developed to account for the thermal and flow physics involved in HRSGs. The CFD model was then verified and validated on a simplified geometry that contained six row-tubes which mimicked the tube arrangements of an HPSH section. Numerical results compared well with published experimental results. In the third step, the oxidation model was coupled to the HRSG’s CFD model, and a comparison study was performed for two cases: (a) regular finless tubes and (b) finned tubes. Preliminary results showed that the oxide thickness predicted by the case with tube fins was higher compared to that of the finless tube case, thus demonstrating the need to model the effects of the tube fins on oxide thickness.

Panicker, Nithin↗
DOE OSTI2020

Xyce Parallel Electronic Simulator Users' Guide (V.7.1)

This manual describes the use of the Xyce Parallel Electronic Simulator. Xyce has been designed as a SPICE-compatible, high-performance analog circuit simulator, and has been written to support the simulation needs of the Sandia National Laboratories electrical designers. This development has focused on improving capability over the current state-of-the-art in the following areas: 1) Capability to solve extremely large circuit problems by supporting large-scale parallel computing platforms (up to thousands of processors). This includes support for most popular parallel and serial computers. 2) A differential-algebraic-equation (DAE) formulation, which better isolates the device model package from solver algorithms. This allows one to develop new types of analysis without requiring the implementation of analysis-specific device models. 3) Device models that are specifically tailored to meet Sandia's needs, including some radiation-aware devices (for Sandia users only). 4) Object-oriented code design and implementation using modern coding practices. Xyce is a parallel code in the most general sense of the phrase a message passing parallel implementation which allows it to run efficiently a wide range of computing platforms. These include serial, shared-memory and distributed-memory parallel platforms. Attention has been paid to the specific nature of circuit-simulation problems to ensure that optimal parallel efficiency is achieved as the number of processors grows.

42 ENGINEERING↗
DOE OSTI2021

Radiation Transport Simulation of the Sulfur Counting System in the Radiation Metrology Laboratory (RML)

At the request of staff members from the Radiation Metrology Laboratory (RML), a series of Monte Carlo radiation transport calculations were performed using two different models of the detector geometry of the RMLs sulfur counting system. The fraction of electrons from each β-decay of 32 P in the sulfur pellet that enter the window of the sulfur counting system was calculated with both MCNP and ITS. In addition, the differential energy distributions of the electrons entering the counting system window were computed. There was significant agreement between the integral and differential quantities calculated by the two transport codes. Summary tables are for the integral efficiency values are found in the body of the report.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Measuring the Density Fields around Bright Quasars at z ~ 6 with XQR-30 Spectra

Measuring the density of the intergalactic medium using quasar sight lines in the epoch of reionization is challenging due to the saturation of Ly α absorption. Near a luminous quasar, however, the enhanced radiation creates a proximity zone observable in the quasar spectra where the Ly α absorption is not saturated. In this study, we use 10 high-resolution ( R ≳ 10,000) z ~ 6 quasar spectra from the extended XQR-30 sample to measure the density field in the quasar proximity zones. We find a variety of environments within 3 pMpc distance from the quasars. We compare the observed density cumulative distribution function (CDF) with models from the Cosmic Reionization on Computers simulation and find a good agreement between 1.5 and 3 pMpc from the quasar. This region is far away from the quasar hosts and hence approaching the mean density of the universe, which allows us to use the CDF to set constraints on the cosmological parameter σ 8 = 0.6 ± 0.3. The uncertainty is mainly due to the limited number of high-quality quasar sight lines currently available. Utilizing the more than 200 known quasars at z ≳ 6, this method will allow us to tighten the constraint on σ 8 to the percent level in the future. In the region closer to the quasar within 1.5 pMpc, we find that the density is higher than predicted in the simulation by 1.23 ± 0.17, suggesting that the typical host dark matter halo mass of a bright quasar ( M 1450 < -26.5) at z ~ 6 is .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

An Investigation into the Effects of Swirl on the Performance and Emissions of an Opposed-Piston Two-Stroke Engine using Large Eddy Simulations

Opposed-piston two-stroke (OP-2S) engines have the potential to achieve higher thermal efficiency than a conventional four-stroke diesel engine. However, the uniflow scavenging process is difficult to control over a wider range of speed and loads due to its sensitivity to pressure dynamics, port timings, and port design. Specifically, the angle of the intake ports can be used to generate swirl which has implications for open and closed cycle effects. This study proposes an analysis of the effects of port angle on the in-cylinder flow distribution and combustion performance of an OP-2S using computational fluid dynamics engine. Large Eddy Simulation (LES) was used to model turbulence given its ability to predict in-cylinder mixing and cyclic variability. A three-cylinder model was validated to experimental data collected by Achates Power and the grid was verified using an LES quality approach from the literature. The model was used to simulate port angles from 12 to 29 degrees at constant pressure and temperature boundary conditions. Results indicated that the higher bulk swirl ratio generated by larger port angles tends to trap more internal residuals. This effect on the scavenging performance, combined with the larger trapped swirl ratio, also has a significant impact on the combustion performance in a two-stroke engine. It was concluded that there exists a tradeoff of efficiency and emissions that must be considered when increasing the port angle in a uniflow two-stroke engine. 1. Introduction Two-stroke engines were developed in the 1800s and are commonly used in practice today for lightweight power applications such as motorcycles or handheld power tools, as well as large-bore low speed engines for stationary power and marine applications. Two-stroke engines provide the distinct advantage of higher torque/power output but present challenges due to its coupled intake and exhaust process. Four-stroke engines have the benefit of using two full piston strokes to induct intake charge and to expel exhaust gases. Controlling the trapped charge composition in a twostroke engine via a scavenging process is crucial for reliable and efficient engine operation.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Neutron Absorber Plate Characterization Plan for Criticality Experiments Design

After being used in nuclear installations, depleted fuel can still be highly reactive and must be handled securely to prevent any radiological or criticality concerns. In particular, spent fuel from use in nuclear power reactors must be stored and transported in specifically designed containers using neutron absorber materials to prevent criticality. Various neutron absorber material types exist and are manufactured by various entities, as thoroughly described in the Handbook of Neutron Absorber Materials for Spent Nuclear Fuel Storage and Transportation Applications written by EPRI. Presently, one of the most modern and most widely used types of neutron absorber material contains particles of boron carbide, or B 4 C, embedded in aluminum matrix: Boralcan, manufactured by Rio Tinto. It is very important for the community to know as much as possible about such neutron absorber materials. Therefore, in the recent years, a US Department of Energy National Nuclear Security Administration–Nuclear Criticality Safety Program funded project initiated design of an experiment that places Boralcan neutron-absorbing plates in an established critical assembly using low-enriched uranium fuel at the Sandia Pulsed Reactor Facility/Critical Experiments (SPRF/CX) apparatus at Sandia National Laboratories. The goal of the experiment is to produce high-quality benchmark data to submit to the International Criticality Safety Benchmark Evaluation Project (ICSBEP), for use in validating calculational tools and nuclear data by criticality safety analysts. The project, named IER-554, is currently in its final design stage, following a successful preliminary design. In the work documented in the design study, ten critical configurations using Boralcan neutron absorber plates were designed, and the experiment was proven to be feasible, with a predicted low k eff uncertainty around 100 pcm. An overview of the modeled cutout of the critical assembly with a Boralcan plate is shown in Figure 1, representing one of the configurations planned for the critical experiments. Before the plates are inserted in the critical assembly, it is necessary to know more about their composition and uniformity. This summary focuses on the plate characterization plans. Each plate will undergo (1) neutron transmission measurements at different locations to determine the 10 B areal density and (2) an in-depth x-ray computed tomography (XCT) examination to obtain the exact Sizes and distribution of the B4C powder particles inside the plates. In parallel, plate modeling studies are performed with a goal to determine the validity of the currently used approximation of modeling the neutron absorber plates as a homogeneous mixture of Aluminum 1100 alloy and B4C— instead of explicitly modeling the B4C particles. By using the experimental 10 B areal density measurements, and the exact size and location of the B4C particles obtained by XCT, a plate model can theoretically be built that reproduces the plate with extremely high fidelity. The results of this modeling study could increase the confidence of the criticality safety community in its modeling methods when using this type of neutron absorber material, and the industry could use these validations to change the boron loading credit limits from the U.S. Nuclear Regulatory Commission standard review plan for dry cask storage of spent nuclear fuel. The modeling calculations are performed with SCALE 6.3.0 using the KENO V.a sequence for criticality calculations with the ENDF/B-VIII.0 continuous-energy cross section library.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Community Input on the Need, Scope, and Development Roadmap of a Proposed National Virtual Biosecurity for Bioenergy Crops Center (NVBBCC)

Brookhaven National Laboratory (BNL) was awarded a pilot project in FY22 under the U.S. Department of Energy (DOE) Office of Science Biopreparedness Research Virtual Environment (BRaVE) initiative, to define research priorities, needs, and requirements for a national virtual center devoted to the biosecurity of bioenergy crops. The mission of the proposed center, referred to as the National Virtual Biosecurity for Bioenergy Crop Center (NVBBCC), would be to provide the scientific basis and tools to detect, characterize, model, and mitigate biothreats to bioenergy crops. This function will be essential to ensure the projected increased US reliance over the next few decades on key plant-based energy products, such as biojet fuel. The NVBBCC is envisioned as a distributed, virtual center with multiple national laboratories at its core to maximize the use of existing unique facilities and expertise across the DOE complex. A major goal of the pilot project was to develop a roadmap for establishing NVBBCC through a series of meetings to gather community input. A total of about 150 individuals, drawn from DOE laboratories, the USDA, academia, NIH, DHS and the private sector participated in six planning meetings held in FY23. Four of the meetings were focused on specific research topics (disease detection, dispersion and disease propagation, biomolecular characterization of plant-pathogen interaction, and mitigation strategies). These four meetings were followed by a meeting that focused on computational needs to support collaborative, data-intensive research within a distributed center as well as workforce development. A final meeting focused on establishing and maintaining preparedness within NVBBCC to respond to an emerging disease within bioenergy crops and how it would collaborate and coordinate with USDA and DHS.

09 BIOMASS FUELS↗
DOE OSTI2020

A Lumped Stochastic Model of Coupled Neutronic Assemblies

The objective of the present work is to apply the Harris ansatz to coupled assemblies in a lumped setting, i.e., ignoring the neutron phase dependence, which hitherto has not been considered. A Master equation for the single-time joint distribution of neutron number in each assembly, accounting for transfer of neutrons from one assembly to another via leakage transition probabilities and intercept or view factors is derived from which equations for the mean and variance in each assembly are developed. Because of the neutronic coupling, it is necessary to additionally compute the correlation function between pairs of assemblies which gives an expanded set of moment equations to solve. The reactivity of each assembly is allowed to vary with time, although a formal mechanism that is responsible for reactivity coupling between assemblies is not considered or proposed. Following Harris, it is then hypothesized that the number distribution in each assembly is a gamma distribution with assembly-specific mean and variance computed exactly from the moment equations. Illustrative numerical results are presented for the moments of a system of four coupled assemblies under various reactivity and neutronic coupling scenarios. The gamma distribution is then used to compute the probability that the neutron population exceeds a threshold value in any assembly. The report concludes with a discussion of recommended further work on this problem.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Evaluating Symmetry in 3-dimensional Simulation Data

The increasing use of three-dimensional simulations has led to an increased need for effective verification and validation, or V&V. Unfortunately, the high computational cost associated with three-dimensional simulations makes many standard V&V techniques cost-prohibitive. Given these circumstances, we must pay close attention to V&V metrics that can be evaluated using only a single simulation. One of the simplest of these metrics is the evaluation of physical symmetries. This report details a method of computing the deviation from symmetry in a three-dimensional dataset with an arbitrarily distributed mesh. Examples are presented for spherical and axisymmetric systems using the Paraview post-processing analysis workflow.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

An approach to grid-to-rod fretting wear modeling using dynamic mortar contact

This article presents the development of a grid-to-rod fretting (GTRF) modeling methodology that utilizes dynamic mortar contact. Here, we leverage a recently developed computational framework for modeling thermomechanical contact in dynamic simulations from the nuclear fuel performance code Bison. Usage of mortar contact ensures a smooth distribution of normal and frictional contact forces, displacements, and velocities on the contact interface, thus facilitating an accurate computation of wear. The integration within an advanced nuclear fuel performance code enables analysis of complex interactions between the fuel, cladding, and spacer grid. Such interactions include swelling, creep, fuel fragmentation, burnup, plenum pressure, and fission gas release. Our methodology is demonstrated via the modeling of two spacer grid geometries and the influence of fuel-cladding mechanical contact on the generation of fretting wear.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

An ODE-Enabled Distributed Transient Stability Analysis for Networked Microgrids

Networked microgrid (NMG) exhibits noteworthy resiliency and flexibility benefits for the mutual support from neighboring microgrids. With high penetration of distributed energy resources (DERs) and the associated controls, the transient stability analysis of NMGs is of critical significance. To address the issues of computation burdens and privacy in the centralized transient analysis, this paper devises an ordinary differential equation (ODE)-enabled distributed transient stability (DTS) methodology for NMGs. First, an ODE-based microgrid model is established to capture the dynamics in the droop control of DERs as well as network and load. Further, a distributed DTS is devised for the ODE representation of an NMG, allowing a privacy-preserving transient analysis of each microgrid while accurately reconstructing the frequency dynamics under droop controls in all DERs. In conclusion, extensive tests are performed to verify the validity of the ODE-based microgrid model through both dynamic response and eigenvalue analysis, and the efficacy of the DTS algorithm in simulating the large signal responses and the frequent fluctuations in NMG.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2026

Protonic nickelate device networks for spatiotemporal neuromorphic computing

Computation in biological neural circuits arises from the interplay of nonlinear temporal responses and spatially distributed dynamic network interactions. Replicating this richness in hardware has remained challenging, as most neuromorphic devices emulate only isolated neuron- or synapse-like functions. Here we introduce an integrated neuromorphic computing platform in which both nonlinear spatiotemporal processing and programmable memory are realized within a single perovskite nickelate material system. By engineering symmetric and asymmetric hydrogenated NdNiO 3 junction devices on the same wafer, we combine ultrafast, proton-mediated transient dynamics with stable multilevel resistance states. Networks of symmetric NdNiO 3 junctions exhibit emergent spatial interactions mediated by proton redistribution, while each node simultaneously provides short-term temporal memory, enabling nanosecond-scale operation with an energy cost of ~0.2 nJ per input. When interfaced with asymmetric output units serving as reconfigurable long-term weights, these networks allow both feature transformation and linear classification in the same material system. Leveraging these emergent interactions, the platform enables real-time pattern recognition and achieves high accuracy in spoken digit classification and early seizure detection, outperforming temporal-only or uncoupled architectures. These results position protonic nickelates as a compact, energy-efficient, CMOS-compatible platform that integrates processing and memory for scalable intelligent hardware.

Electrical and electronic engineering↗
DOE OSTI2024

DASSH-F: Subchannel Based Thermal Analysis

The DASSH thermal analysis code is designed to rapidly allow a reactor design engineer to obtain flow rates requirements that satisfy peak temperature constraints in the domain. The advantage of using DASSH over a hand calculation is that it has a more rigorous treatment of the pin power distribution and coolant heat transfer within an assembly and between assemblies. The advantage of using DASSH over a conventional 3D subchannel code or a computational fluid dynamics code (CFD) is that it can obtain the desired solution in a matter of minutes in serial with minor computer memory needs. The DASSH methodology is virtually identical to SUPERENERGY-2 with additional functionalities taken from follow on work to SUPERENERGY-2 done at ANL in the 1980s. DASSH today is an integral component of the Argonne Fast Reactor analysis suite for reactor design work. DASSH obtains the power distribution from a coupled neutron-gamma heating calculation in GAMSOR (including DIF3D) at each time point of a companion fuel cycle analysis calculation with REBUS. The domain in DASSH assumes a hexagonal grid typical for fast reactors with much of the geometry information taken from the DIF3D model. DASSH assumes the assemblies that are loaded into each grid position are ducted to control the coolant flow. Considerable detail is given on the subchannel formulation of DASSH in this document. Much of the formulation and design of the code builds upon research done by previous authors with little new investigation. Thus the decisions made in developing the subchannel model used in DASSH have their origins over 50 years ago. Much of the heat transfer methodology in DASSH is built upon correlations for both the coolant mixing and heat transfer coefficients for pins and ducts. DASSH is thus not a rigorous treatment of a given problem, but a rapid assessment of the temperature field that has known limitations with respect to an experimental measurement or CFD calculation. The DASSH input and output are detailed along with usage of the software. The DASSH output provides tables of evaluated material properties and key coolant and pin temperature results. DASSH can create Python scripts that generate domain summary pictures. DASSH can also generate assembly temperature maps and VTK output files which allow the DASSH solution to be visualized. As the primary purpose of the DASSH software is to compute the coolant and fuel pin temperature distribution for a given model of a reactor, much of the output focus is giving the user quick summary tables needed to assess the performance of a given orifice flow specification. The present version of DASSH has a crude orifice search capability and an efficient orifice flow search capability. The flow search tries to meet user specified constraints for 1) peak 2-sigma clad midwall temperature, 2) peak coolant temperature, and 3) desired bulk outlet temperature. There are many development shortcomings in DASSH detailed in this document, but this version is functional for modern analysis needs. This document serves as the manual for the Fortran based DASSH software that was developed to replace the Python version of DASSH developed as part of the VTR program.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Novel Geometric Operations for Linear Programming

This report summarizes the work performed under the project "Linear Programming in Strongly Polynomial Time." Linear programming (LP) is a classic combinatorial optimization problem heavily used directly and as an enabling subroutine in integer programming (IP). Specifically IP is the same as LP except that some solution variables must take integer values (e.g. to represent yes/no decisions). Together LP and IP have many applications in resource allocation including general logistics, and infrastructure design and vulnerability analysis. The project was motivated by the PI's recent success developing methods to efficiently sample Voronoi vertices (essentially finding nearest neighbors in high-dimensional point sets) in arbitrary dimension. His method seems applicable to exploring the high-dimensional convex feasible space of an LP problem. Although the project did not provably find a strongly-polynomial algorithm, it explored multiple algorithm classes. The new medial simplex algorithms may still lead to solvers with improved provable complexity. We describe medial simplex algorithms and some relevant structural/complexity results. We also designed a novel parallel LP algorithm based on our geometric insights and implemented it in the Spoke-LP code. A major part of the computational step is many independent vector dot products. Our parallel algorithm distributes the problem constraints across processors. Current commercial and high-quality free LP solvers require all problem details to fit onto a single processor or multicore. Our new algorithm might enable the solution of problems too large for any current LP solvers. We describe our new algorithm, give preliminary proof-of-concept experiments, and describe a new generator for arbitrarily large LP instances.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

An Efficient Bayesian Approach to Learning Droplet Collision Kernels: Proof of Concept Using “Cloudy,” a New n -Moment Bulk Microphysics Scheme

The small-scale microphysical processes governing the formation of precipitation particles cannot be resolved explicitly by cloud resolving and climate models. Instead, they are represented by microphysics schemes that are based on a combination of theoretical knowledge, statistical assumptions, and fitting to data (“tuning”). Historically, tuning was done in an ad hoc fashion, leading to parameter choices that are not explainable or repeatable. Recent work has treated it as an inverse problem that can be solved by Bayesian inference. The posterior distribution of the parameters given the data—the solution of Bayesian inference—is found through computationally expensive sampling methods, which require over $\mathcal{O}$(10 5 ) evaluations of the forward model; this is prohibitive for many models. We present a proof of concept of Bayesian learning applied to a new bulk microphysics scheme named “Cloudy,” using the recently developed Calibrate-Emulate-Sample (CES) algorithm. Cloudy models collision-coalescence and collisional breakup of cloud droplets with an adjustable number of prognostic moments and with easily modifiable assumptions for the cloud droplet mass distribution and the collision kernel. The CES algorithm uses machine learning tools to accelerate Bayesian inference by reducing the number of forward evaluations needed to $\mathcal{O}$(10 2 ). It also exhibits a smoothing effect when forward evaluations are polluted by noise. In a suite of perfect-model experiments, we show that CES enables computationally efficient Bayesian inference of parameters in Cloudy from noisy observations of moments of the droplet mass distribution. In an additional imperfect-model experiment, a collision kernel parameter is successfully learned from output generated by a Lagrangian particle-based microphysics model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Kinetic Modeling and Experimental Studies of Capacitively Coupled Macro- and Microplasmas at Microwave Frequencies

The effect of driving frequency in the range of 13.56 MHz to 73 MHz on electron energy distribution and electron heating modes in a 50 mTorr capacitively coupled argon plasma discharge is studied using 1D-3V particle-in-cell simulations. Calculated electron energy probability functions exhibit three distinct temperatures for low-, mid-, and high-energy electrons at all the studied driving frequencies. When compared to published experimental data, the calculated probability functions show a reasonable agreement for the energy range resolved in the measurements (about 2 eV to 10 eV). Discrepancies due to limitations in experimental energy resolution outside this range lead to differences between computational and experimental values of the electron number density determined from the distribution functions, and the predicted effective electron temperature is within 25 % of experimental values. The impedance of the discharge is interpreted in terms of a homogeneous equivalent circuit model and the driving frequency dependence of the inferred combined sheath thickness is found to obey a known, theoretically-derived, power law. The average power transferred from the field to the electrons (electron heating) is computed, and a region of negative heating near the sheath edge, particularly at higher driving frequencies, is identified. Analysis of the electron momentum equation shows that electron inertia, which on temporal averaging would be zero in a linear regime, is responsible for negative values of power deposition near the sheath edge at high driving frequencies due to the highly nonlinear behavior of the discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

An improved specimen preparation method for marine shallow gas-bearing sand sediments and its validations

Instabilities of shallow gas-charged seabed are potential geological hazards in ocean engineering. In practice, the conventional field sampling techniques failed to obtain undisturbed gas-bearing sediments from the seabed for laboratory mechanical testing because of sensitive gas exsolution and escape from sediments. However, preparation of representative remoulded gas-charged specimens is a challenging issue, because it is rather difficult to quantitatively control the gas content and obtain uniform distribution of gas bubbles within the specimen. Given the above problems, this work proposes a reliable approach to reconstitute the high-saturation specimen of gas-charged sediments in the laboratory by an improved multifunction integrated triaxial apparatus (MITA). This apparatus is developed based on an advanced stress path triaxial system by introducing a temperature-controlled system and a wave-monitoring system. The temperature-controlled system is used to accurately mimic the in situ environments of sediments in the seabed. The wave-monitoring system is utilized to identify exsolution point of free gas and examine the disturbance of gas to specimens during gas exsolution. The detailed procedure of gassy specimen preparation is introduced. Then, the quality of prepared specimens using our improved apparatus is validated by the high-resolution micro-X-ray computed tomography (μCT) scanning test, from which bubble occurrence and size distribution within the gassy sand specimen can be obtained; and preliminary mechanical tests on gassy sand specimens with various initial saturation degrees are performed. The proposed specimen preparation procedure succeeds in proving the postulated occurrence state of gas bubbles in coarse-grained sediments and accurately controlling the gas content.

58 GEOSCIENCES↗