MFIX-Exa: CFD-DEM Simulations of Thermodynamics and Chemical Reactions in Multiphase Flows
For presentation at the AIChE Annual Meeting, Orlando, FL, November 5-10, 2023.
Porcu, Roberto↗
SEARCH · Engineering Papers
Results for “time integrator”
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.
PETSc/TAO Users Manual (Rev. 3.19)
This manual describes the use of the Portable, Extensible Toolkit for Scientific Computation (PETSc) and the Toolkit for Advanced Optimization (TAO) for the numerical solution of partial differential equations and related problems on high-performance computers. PETSc/TAO is a suite of data structures and routines that provide the building blocks for the implementation of large-scale application codes on parallel (and serial) computers. PETSc uses the MPI standard for all distributed memory communication. PETSc/TAO includes a large suite of parallel linear solvers, nonlinear solvers, time integrators, and opti mization that may be used in application codes written in Fortran, C, C++, and Python (via petsc4py; see Getting Started). PETSc provides many of the mechanisms needed within parallel application codes, such as parallel matrix and vector assembly routines. The library is organized hierarchically, enabling users to employ the level of abstraction that is most appropriate for a particular problem. By using techniques of object-oriented programming, PETSc provides enormous flexibility for users. PETSc is a sophisticated set of software tools; as such, for some users it initially has a much steeper learning curve than packages such as MATLAB or a simple subroutine library. In particular, for individuals without some computer science background, experience programming in C, C++, python, or Fortran and experience using a debugger such as gdb or lldb, it may require a significant amount of time to take full advantage of the features that enable efficient software use. However, the power of the PETSc design and the algorithms it incorporates may make the efficient implementation of many application codes simpler than “rolling them” yourself. For many tasks a package such as MATLAB is often the best tool; PETSc is not intended for the classes of problems for which effective MATLAB code can be written. There are several packages, built on PETSc, that may satisfy your needs without requiring directly using PETSc. We recommend reviewing these packages functionality before starting to code directly with PETSc. PETSc can be used to provide a “MPI parallel linear solver” in an otherwise sequential, or OpenMP parallel code. This approach cannot provide extremely large improvements in the application time by utilizing large numbers of MPI processes but can still improve the performance. Certainly all parts of a previously sequential code need not be parallelized but the matrix generation portion must be parallelized to expect true scalability to large numbers of MPI processes. See PCMPI for details on how to utilize the PETSc MPI linear solver server. Since PETSc is under continued development, small changes in usage and calling sequences of routines will occur. PETSc has been supported for twenty-five years; see mailing list information on our website for information on contacting support.
97 MATHEMATICS AND COMPUTING↗
Advancements of On-Line Monitoring and Real-Time Characterization of Actual Low-Level Hanford Tank-Waste Samples
Waste and process stream characterization is essential to the safe, efficient, and cost-effective processing of Hanford Site materials. The demand for materials characterization including chemical, radionuclide, and physical attributes during waste management, transfer, and staging operations in the Hanford tank farms is expected to significantly increase as the U.S. Department of Energy River Protection Project progresses toward Hanford waste treatment and immobilization. Our past work has applied Raman on-line monitoring and chemometric modeling to the supernate of Hanford tank 241-AP-105 and has demonstrated the quantitative measurement of nine analytes within this waste. To demonstrate broader applicability, this method has now been expanded to waste from two additional Hanford tanks, AW-102 and AP-107, and the offgas condensate from lab-scale melter runs of actual-tank waste feeds. To enhance the performance of the Raman method, instruments with various excitation wavelengths were compared, since it is known that the Raman response is enhanced by shorter wavelength excitation. The laser excitations for the three systems tested were 404 nm, 532 nm, and 671 nm. In addition to laser wavelength selection, the laser power, measurement integration times, and signal averaging techniques were also investigated to determine their effect on detection limits for oxy-anion analytes within tank wastes. These advancements in Raman capability are compared with past standards and will be presented.
Bryan, Samuel A.↗
An improved Method for Determining Transient Fuel Dilution of Oil in an Internal-Combustion Engine Using Laser-Induced Florescence and Multivariate Least Square Calibration
An optical diagnostic, based on laser-induced fluorescence (LIF), has been developed for on-engine measurements of real-time fuel dilution of engine oil or fuel in oil (FIO). Fuel dilution of oil is broadly relevant to advancing engine technology including durability, calibration, and catalyst-system management, and believed to promote destructive stochastic pre-ignition (SPI) during high-load engine operations. While standard (e.g., ASTM D3524-90) methods are not capable of real-time transient measurements, the LIF technique resolves transient dilution on the minutes time scale. We have expanded on our original FIO instrument development by introducing an improved analysis based on multivariate least square chemometrics analysis. The measurement uses a fuel dye (180–1300 parts per million, by mass) and monitors for its presence in the oil using 532 nm excitation and LIF. While the original FIO instrument utilized a two-color ratio method for analysis, the improved chemometric analysis uses the fully resolved LIF dye spectra to provide better predictive FIO accuracy (>92%) over a wide FIO range (1.5–14%) typical of engine application. We also investigate the effect of oil temperature on the LIF signal. Limited engine applications for demonstrating and validating the improved FIO instrument are shown, and the related data used to quantify practical detection limit and sensitivity. The improved analysis is insensitive to laser power fluctuation and change in detector integration time, providing an excellent FIO sensitivity (1–2%) and detection limit (0.01 %FIO) over a wide range of loads and injection timings, illustrating this updated approach to be a promising tool for advancing engine technology.
42 ENGINEERING↗
Upgrades and redesign of the National Ignition Facility’s soft x-ray opacity spectrometer (OpSpec)
The soft x-ray Opacity Spectrometer (OpSpec) used on the National Ignition Facility (NIF) has recently incorporated an elliptically shaped crystal. The original OpSpec used two convex cylindrical crystals for time-integrated measurements of point-projection spectra from 540 to 2100 eV. However, with the convex geometry, the low-energy portion of the spectrum suffered from high backgrounds due to scattered x-rays as well as reflections from alternate crystal planes. An elliptically shaped crystal allows an acceptance aperture at the crossover focus between the crystal and the detector, which reduces background and eliminates nearly all reflections from alternate crystal planes. The current elliptical design is an improvement from the convex cylindrical design but has a usable energy range from 900 to 2100 eV. In addition, OpSpec is currently used on 18 NIF shots/year, in which both crystals are typically damaged beyond reuse, so efficient production of 36 crystals/year is required. Design efforts to improve the existing system focus on mounting reliability, reducing crystal strain to increase survivability between mounting and shot time, and extending the energy range of the instrument down to 520 eV. The elliptical design, results, and future options are presented.
Wallace, M. S. (ORCID:0000000310871216)↗
All-sky Neutrino Point-source Search with IceCube Combined Track and Cascade Data
Despite extensive efforts, discovery of high-energy astrophysical neutrino sources remains elusive. We present an event-level simultaneous maximum likelihood analysis of tracks and cascades using IceCube data collected from 2008 April 6 to 2022 May 23 to search the whole sky for neutrino sources, and using a source catalog, for coincidence of neutrino emission with gamma-ray emission. This is the first time a simultaneous fit of different detection channels is used to conduct a time-integrated all-sky scan with IceCube. Combining all-sky tracks, with superior pointing power and sensitivity in the northern sky, with all-sky cascades, with good energy resolution and sensitivity in the southern sky, we have developed the most sensitive point-source search to date by IceCube that targets the entire sky. The most significant point in the northern sky aligns with NGC 1068, a Seyfert II galaxy, which, from the catalog search, shows a 3.5σ excess over background after accounting for trials. The most significant point in the southern sky does not align with any source in the catalog and is not significant after accounting for trials. A search for the single most significant Gaussian flare at the locations of NGC 1068, PKS 1424+240, and the southern highest-significance point shows results consistent with expectations for steady emission. Notably, this is the first time that a flare shorter than four years has been excluded as being responsible for NGC 1068’s emergence as a neutrino source. Our results show that combining tracks and cascades when conducting neutrino source searches improves sensitivity and can lead to new discoveries.
Abbasi, R. [Loyola University, Chicago, IL (United↗
Calculation of electron-ion temperature equilibration rates and friction coefficients in plasmas and liquid metals using quantum molecular dynamics
In this work, we discuss a method to calculate with quantum molecular dynamics simulations the rate of energy exchanges between electrons and ions in two-temperature plasmas, liquid metals, and hot solids. Promising results from this method were recently reported for various materials and physical conditions. Like other ab initio calculations, the approach offers a very useful comparison with the experimental measurements and permits an extension into conditions not covered by the experiments. The energy relaxation rate is related to the friction coefficients felt by individual ions due to their nonadiabatic interactions with electrons. Each coefficient satisfies a Kubo relation given by the time integral of the autocorrelation function of the interaction force between an ion and the electrons. These Kubo relations are evaluated using the output of quantum molecular dynamics calculations in which electrons are treated in the framework of finite-temperature density functional theory. The calculation presents difficulties that are unlike those encountered with the Kubo formulas for the electrical and thermal conductivities. In particular, the widely used Kubo-Greenwood approximation is inapplicable here. Indeed, the friction coefficients and the energy relaxation rate diverge in this approximation since it does not properly account for the electronic screening of electron-ion interactions. The inclusion of screening effects considerably complicates the calculations. We discuss the physically motivated approximations we applied to deal with these complications in order to investigate a widest range of materials and physical conditions. Unlike the standard method used for the electronic conductivities, the Kubo formulas are evaluated directly in the time domain and not in the energy domain, which spares one from needing to introduce an extraneous undetermined numerical parameter to account for the discrete character of the numerical density of states. We highlight interesting properties of the energy relaxation rate not shared by other electronic properties, in particular its self-averaging character. We then present a detailed parametric and convergence study with the numerical parameters, including the system size, the number of bands and k-points, and the physical approximations for the dielectric function and the exchange-correlation energy.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Review of Radio Frequency Interference and Potential Impacts on the CMB-S4 Cosmic Microwave Background Survey
CMB-S4 will map the cosmic microwave background to unprecedented precision, while simultaneously surveying the millimeter-wave time-domain sky, in order to advance our understanding of cosmology and the universe. CMB-S4 will observe from two sites, the South Pole and the Atacama Desert of Chile. A combination of small- and large-aperture telescopes with hundreds of thousands of polarization-sensitive detectors will observe in several frequency bands from 20-300 GHz, surveying more than 50% of the sky to arcminute resolution with unprecedented sensitivity. CMB-54 seeks to make a dramatic leap in sensitivity while observing across a broad range of largely unprotected spectrum which is increasingly being utilized for terrestrial and satellite transmissions. Fundamental aspects of CMB instrument technology leave them vulnerable to radio frequency interference (RFI) across a wide range of frequencies, including frequencies outside of their observing bands. Ground-based CMB instruments achieve their extraordinary sensitivities by deploying large focal planes of superconducting bolometers to extremely dry, high-altitude sites, with large fractional bandwidths, wide fields of view, and years of integration time. Suitable observing sites have historically offered significant protection from RFI, both naturally through their extremely remote locations as well as through restrictions on local emissions. Since the coupling mechanisms are complex, "safe" levels or frequencies of emission that would not interfere with CMB measurements cannot always be determined through straightforward calculations. We discuss models of interference for various types of RFI relevant to CMB-54, mitigation strategies, and the potential impacts on survey sensitivity.
Cosmic Microwave Background↗
Contrast and Temperature Dependence of Multi-epoch High-resolution Cross-correlation Exoplanet Spectroscopy
While high-resolution cross-correlation spectroscopy (HRCCS) techniques have proven effective at characterizing the atmospheres of transiting and nontransiting hot Jupiters, the limitations of these techniques are not well understood. We present a series of simulations of one HRCCS technique, which combines the cross-correlation functions from multiple epochs, to place temperature and contrast limits on the accessible exoplanet population for the first time. We find that planets approximately Saturn-sized and larger within ∼0.2 au of a Sun-like star are likely to be detectable with current instrumentation in the L band, a significant expansion compared with the previously studied population. Cooler (T {sub eq} ≤ 1000 K) exoplanets are more detectable than suggested by their photometric contrast alone as a result of chemical changes that increase spectroscopic contrast. The L-band CH{sub 4} spectrum of cooler exoplanets enables robust constraints on the atmospheric C/O ratio at T {sub eq} ∼ 900 K, which have proven difficult to obtain for hot Jupiters. These results suggest that the multi-epoch approach to HRCCS can detect and characterize exoplanet atmospheres throughout the inner regions of Sun-like systems with existing high-resolution spectrographs. We find that many epochs of modest signal-to-noise ratio (S/N{sub epoch} ∼ 1500) yield the clearest detections and constraints on C/O, emphasizing the need for high-precision near-infrared telluric correction with short integration times.
46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
A large deformation multiphase continuum mechanics model for shock loading of soft porous materials
A large deformation, coupled finite-element (FE) model is developed to simulate the multiphase response of soft porous materials subjected to high strain-rate loading. The approach is based on the theory of porous media (TPM) at large deformations. Simplifications to the one-dimensional regime studied in the numerical simulations follow. An overview of several different time integration schemes is presented for the purpose of solving the nonlinear dynamic coupled balance of momenta (mixture and fluid) and balance of mass of the mixture equations. Numerical examples are presented for (i) verification against closed-form analytical solutions assuming small loads, (ii) demonstrating large deformation effects at high strain-rate, and (iii) showing differences in deformations between a single-phase elastodynamics model with occluded compressible pore fluid and a multiphase poroelastodynamics model at high strain-rate. The multiphase model shows that the relative motion of the pore fluid significantly dampens the deformation response of the solid skeleton as compared to the single-phase model, and makes it possible to extract quantitative values for the stresses of the different constituents, thereby allowing one to form preliminary conclusions about the onset of damage in the solid skeleton. The novelty of the current work is developing a multiphase, large deformation, mixture theory numerical model for high strain-rate loading of soft porous materials. It was discovered that explicit, adaptive time-stepping Runge–Kutta schemes offer high accuracy at relatively low cost when compared to traditional implicit or explicit central difference time-stepping schemes for shock-like loadings. Here, shock viscosity is added to the mixture momentum balance equation to regularize the shock front, and a stabilization term is added to the mixture mass balance equation to stabilize equal order interpolation finite elements for the coupled finite element solution of multiphase materials.
Engineering↗
Ultrafast time-resolved 2D imaging of laser-driven fast electron transport in solid density matter using an x-ray free electron laser
High-power, short-pulse laser-driven fast electrons can rapidly heat and ionize a high-density target before it hydrodynamically expands. The transport of such electrons within a solid target has been studied using two-dimensional (2D) imaging of electron-induced Kα radiation. However, it is currently limited to no or picosecond scale temporal resolutions. Here, we demonstrate femtosecond time-resolved 2D imaging of fast electron transport in a solid copper foil using the SACLA x-ray free electron laser (XFEL). An unfocused collimated x-ray beam produced transmission images with sub-micron and ~10 fs resolutions. The XFEL beam, tuned to its photon energy slightly above the Cu K-edge, enabled 2D imaging of transmission changes induced by electron isochoric heating. Time-resolved measurements obtained by varying the time delay between the x-ray probe and the optical laser show that the signature of the electron-heated region expands at ~25% of the speed of light in a picosecond duration. Time-integrated Cu Kα images support the electron energy and propagation distance observed with the transmission imaging. The x-ray near-edge transmission imaging with a tunable XFEL beam could be broadly applicable for imaging isochorically heated targets by laser-driven relativistic electrons, energetic protons, or an intense x-ray beam.
47 OTHER INSTRUMENTATION↗
A parallel p ‐adaptive discontinuous Galerkin method for the Euler equations with dynamic load‐balancing on tetrahedral grids
Abstract A novel p ‐adaptive discontinuous Galerkin (DG) method has been developed to solve the Euler equations on three‐dimensional tetrahedral grids. Hierarchical orthogonal basis functions are adopted for the DG spatial discretization while a third order TVD Runge‐Kutta method is used for the time integration. A vertex‐based limiter is applied to the numerical solution in order to eliminate oscillations in the high order method. An error indicator constructed from the solution of order and is used to adapt degrees of freedom in each computational element, which remarkably reduces the computational cost while still maintaining an accurate solution. The developed method is implemented with under the Charm++ parallel computing framework. Charm++ is a parallel computing framework that includes various load‐balancing strategies. Implementing the numerical solver under Charm++ system provides us with access to a suite of dynamic load balancing strategies. This can be efficiently used to alleviate the load imbalances created by p ‐adaptation. A number of numerical experiments are performed to demonstrate both the numerical accuracy and parallel performance of the developed p ‐adaptive DG method. It is observed that the unbalanced load distribution caused by the parallel p ‐adaptive DG method can be alleviated by the dynamic load balancing from Charm++ system. Due to this, high performance gain can be achieved. For the testcases studied in the current work, the parallel performance gain ranged from 1.5× to 3.7×. Therefore, the developed p ‐adaptive DG method can significantly reduce the total simulation time in comparison to the standard DG method without p ‐adaptation.
97 MATHEMATICS AND COMPUTING↗
Assessment of Potential Dose and Environmental Impacts from Proposed Testing at the INL Radiological Response Training Range
This assessment uses screening-level models to calculate potential environmental impacts from proposed tests at the Idaho National Laboratory (INL) Radiological Response Training Range (RRTR) site. Proposed tests could be conducted using 11 different radioactive material types that include K2O, LaBr3, KBr, Cu, Zr, F, Ga, Ga2O3, NaNO2, Ga-68, and Tc-99m. The tests could potentially release radioactive material to the atmosphere and radionuclides and other contaminants to the soil, which could leach into the unsaturated zone and migrate to the aquifer. Atmospheric transport of radionuclides to potential human receptors and time-integrated air concentrations were calculated with a Gaussian plume model and three years of hourly meteorological data. Potential surface-soil impacts were calculated with the computer program mixing-cell model (MCM). Groundwater impacts were calculated with the computer programs MCM and GWSCREEN. Radiological doses from potential atmospheric releases were calculated for public receptors off the INL Site and for workers at nearby INL facilities. Results were compared to regulatory dose limits. Maximum potential groundwater concentrations were estimated in the aquifer below the NSTR site and compared to drinking water standards or risk-based screening levels for resident tap water. Soil concentrations were calculated and compared to risk-based screening levels for workers and potential future residents. All impacts were estimated assuming 12 tests are conducted annually using all 11 material types for a period of 15 years. This document provides the resources to enable a subject matter expert in the field of environmental assessments to replicate the modeling and calculations. The methodology and parameters are presented in the text. All electronic files, including computer-code input, output, executable files, batch files, scripts, and spreadsheet files are contained in a zip file that can be accessed by selecting “Additional Information” (select Native File) in the INL Electronic Document Management System (EDMS). It is highly unlikely the test scenarios evaluated in this ECAR will adversely impact human health based on comparisons of calculated dose and concentration against regulatory standards and risk-based screening levels. Conservative estimates of dose to workers and the public from atmospheric transport of possible radionuclide releases are far below federal radiation protection standards. Conservative estimates of potential contaminant concentrations in groundwater are less than federal drinking water standards or screening levels. Predicted radionuclide concentrations in surface soils are below risk-based screening levels, except for Ge-68 (material Ga-68) for the worker. The Ge-68 soil concentration can be made less than the worker PRG, if the number of annual tests using Ga-68 is reduced from 12 to 6. However, the sum of ratios still exceeds one because of the high K-40 ratio. If the EF of the worker (number of days the worker is in the contaminated testing area) is reduced from 225 days/yr (default value for full time worker) to 112 days/yr, the Ge-68 ratio is less than one and the sum of ratios is less than one. Actual radiation doses and groundwater and surface-soil concentrations are likely to be much less than those calculated because of the conservative assumptions and parameters employed in the modeling. For example, atmospheric-transport calculations assume the entire inventory of each material type is readily released to the atmosphere and no plume deposition, depletion, or radioactive decay occurs during transport. The calculations also assume the same meteorological conditions (e.g., wind velocity, wind direction, stability class) that produce the maximum 95th percentile concentration (i.e., concentration representing the 95th percentile of a distribution of concentrations derived from 3 years of hourly meteorological data) at each receptor location are the same for all 12 tests during the year, and each receptor is assumed to be present during all 12 tests. The surface-soil assessment assumes the entire inventory of each test is deposited in the top 5 cm of soil. No atmospheric dispersal is assumed, and the radionuclides are subject only to leaching and radioactive decay. The groundwater-pathway modeling is conservative in that it is one-dimensional in the unsaturated zone (no lateral spreading/dilution) and assumes the entire inventory of contaminants infiltrates into the ground at the same location for every test. This is especially conservative for particulate radionuclides because they would have to dissolve or corrode first and some would be dispersed into the atmosphere. The groundwater receptor is also assumed to consume water directly from a hypothetical well positioned in the location of maximum concentration. In addition, conservative degradation rates were used, and volatilization was not considered for the nonradioactive chemical
54 ENVIRONMENTAL SCIENCES↗
Assessment of Potential Dose and Environmental Impacts from Proposed Testing at the INL National Security Test Range
This assessment uses screening-level models to calculate potential environmental impacts from proposed tests at two locations at the Idaho National Laboratory (INL) National Security Test Range (NSTR) site. Proposed tests could be conducted using 11 different radioactive material types that include K 2 O, LaBr 3 , KBr, Cu, Zr, F, Ga, Ga 2 O 3 , NaNO 2 , Ga-68, and Tc-99m. The tests could potentially release radioactive material to the atmosphere and radionuclides and other contaminants to the soil, which could leach into the unsaturated zone and migrate to the aquifer. Atmospheric transport of radionuclides to potential human receptors and time-integrated air concentrations were calculated with a Gaussian plume model and three years of hourly meteorological data. Potential surface-soil impacts were calculated with the computer program mixing-cell model (MCM). Groundwater impacts were calculated with the computer programs MCM and GWSCREEN. Radiological doses from potential atmospheric releases were calculated for public receptors off the INL Site and for workers at nearby INL facilities. Results were compared to regulatory dose limits. Maximum potential groundwater concentrations were estimated in the aquifer below the NSTR site and compared to drinking water standards or risk-based screening levels for resident tap water. Soil concentrations were calculated and compared to risk-based screening levels for workers and potential future residents. All impacts were estimated based on the assumption that 12 tests are conducted annually using all 11 material types for 15 years. This document provides the resources to enable a subject matter expert in the field of environmental assessments to replicate the modeling and calculations. The methodology and parameters are presented in the text. All electronic files, including computer-code input, output, executable files, batch files, scripts, and spreadsheet files, are contained in a zip file that can be accessed by selecting “Additional Information” (select Native File) in the INL Electronic Document Management System (EDMS). It is highly unlikely the test scenarios evaluated in this ECAR will adversely impact human health based on comparisons of calculated dose and concentration against regulatory standards and risk-based screening levels. Conservative estimates of dose to workers and the public from atmospheric transport of possible radionuclide releases are far below federal radiation protection standards. Conservative estimates of potential contaminant concentrations in groundwater are less than federal drinking water standards or screening levels. Predicted radionuclide concentrations in surface soils are below risk-based screening levels, except for Ge-68 (material Ga-68) for the worker. The Ge-68 soil concentration can be made less than the worker PRG, if the number of annual tests using Ga-68 is reduced from 12 to 6. However, the sum of ratios still exceeds one because of the high K-40 ratio. If the EF of the worker (number of days the worker is in the contaminated testing area) is reduced from 225 days/yr (default value for full time worker) to 112 days/yr, the Ge-68 ratio is less than one and the sum of ratios is less than one. Actual radiation doses and groundwater and surface-soil concentrations are likely to be much less than those calculated because of the conservative assumptions and parameters employed in the modeling. For example, atmospheric-transport calculations assume the entire inventory of each material type is readily released to the atmosphere and no plume deposition, depletion, or radioactive decay occurs during transport. The calculations also assume the same meteorological conditions (e.g., wind velocity, wind direction, stability class) that produce the maximum 95th percentile concentration (i.e., concentration representing the 95th percentile of a distribution of concentrations derived from 3 years of hourly meteorological data) at each receptor location are the same for all 12 tests during the year, and each receptor is assumed to be present during all 12 tests. The surface-soil assessment assumes the entire inventory of each test is deposited in the top 5 cm of soil. No atmospheric dispersal is assumed, and the radionuclides are subject only to leaching and radioactive decay. The groundwater-pathway modeling is conservative in that it is one-dimensional in the unsaturated zone (no lateral spreading/dilution) and assumes the entire inventory of contaminants infiltrates into the ground at the same location for every test. This is especially conservative for particulate radionuclides because they would have to dissolve or corrode first and some would be dispersed into the atmosphere. The groundwater receptor is also assumed to consume water directly from a hypothetical well positioned in the location of maximum concentration. In addition, conservative degradation rates were used, and volatilization was not considered for the nonradradioactive chemicals modeled. And finally, the calculations assume all 12 tests will be performed at the same place at both locations, and all 11 radioactive material types will be used for each test. This is conservative because it is anticipated that no more than two material types will be used per test.
99 GENERAL AND MISCELLANEOUS↗
Asteroids seen by JWST-MIRI: Radiometric size, distance, and orbit constraints
Infrared measurements of asteroids are crucial for the determination of physical and thermal properties of individual objects, and for understanding the small-body populations in the solar system as a whole. However, standard radiometric methods can only be applied if the orbit of an object is known, hence its position at the time of the observation. With JWST-MIRI observations the situation will change and many unknown, often very small, solar system objects will be detected. Later orbit determinations are difficult due to the faintness of the objects and the lack of dedicated follow-up concepts. We present MIRI observations of the outer-belt asteroid (10920) 1998 BC1 and an unknown object, detected in all nine MIRI bands in close apparent proximity to (10920). We developed a new method called STM-ORBIT to interpret the multi-band measurements without knowing the object’s true location. The power of the new technique is that it determines the most-likely heliocentric and observer-centric distance and phase angle ranges, allowing us to make a radiometric size estimate. The application to the MIRI fluxes of (10920) was used to validate the method. It leads to a confirmation of the known radiometric size-albedo solution, and puts constraints on the asteroid’s location and orbit in agreement with its true orbit. To back up the validation of the method, we obtained additional ground-based light curve observations of (10920), combined with Gaia data, which indicate a very elongated object (a/b ≥ 1.5), with a spin-pole at (λ, β) ecl = (178°, +81°), with an estimated error of about 20°, and a rotation period of 4.861191 ± 0.000015 h. A thermophysical study of all available JWST-MIRI and WISE measurements leads to a size of 14.5–16.5 km (diameter of an equal-volume sphere), a geometric albedo pV between 0.05 and 0.10, and a thermal inertia in the range 9–35 (best value 15) J m -2 s -0.5 K -1 . For the newly discovered MIRI object, the STM-ORBIT method revealed a size of 100–230 m. The new asteroid must be on a low-inclination orbit (0.7° < i < 2.0°) and it was located in the inner main-belt region during JWST observations. A beaming parameter η larger than 1.0 would push the size even below 100 meters, a main-belt regime that has escaped IR detections so far. This kind of MIRI observations can therefore contribute to formation and evolution studies via classical size-frequency studies, which are currently limited to objects larger than about one kilometer in size. We estimate that MIRI frames with pointings close to the ecliptic and short integration times of only a few seconds will always include a few asteroids; most of them will be unknown objects.
79 ASTRONOMY AND ASTROPHYSICS↗
Room‐Temperature Mid‐Infrared Detection Using Metasurface‐Absorber‐Integrated Phononic Crystal Oscillator
Mid-infrared (MIR) detectors find extensive applications in chemical sensing, spectroscopy, communications, biomedical diagnosis, and space exploration. Alternative to semiconductor MIR photodiodes and bolometers, mechanical-resonator-based MIR detectors show advantages in higher sensitivity and lower noise at room temperature, especially toward longer wavelength infrared. Here, uncooled room-temperature MIR detectors based on lithium niobate surface acoustic wave phononic crystal (PnC) resonators integrated with wavelength-and-polarization-selective metasurface absorber arrays are demonstrated. The detection is based on the resonant frequency shift induced by the local temperature change due to MIR absorptions. The PnC resonator is configured in an oscillating mode, enabling active readout and low-frequency noise. The 1-GHz oscillator-based MIR detector shows a relative frequency deviation of 5.24 × 10 −10 Hz −1/2 at an integration time of 50 µs, leading to an incident noise equivalent power of 197 pW Hz −1/2 when input 6-µm MIR light is modulated at 1.8 kHz, and a large dynamic range of 10 7 in incident MIR power. The device architecture is compatible with the scalable manufacturing process and can be readily extended to a broader spectral range by tailoring the absorbing wavelengths of metasurface absorbers.
Xi, Zichen [Virginia Polytechnic Inst. and State U↗
Spectral deferred correction methods for high-order accuracy in poroelastic problems
In this work, we investigate high-order accuracy in time integration by examining two operator splitting methods for poroelastic problems: the two-pass and the spectral deferred correction (SDC) methods. To enhance the order of accuracy, the two-pass method partitions a coupled operator symmetrically, whereas the SDC method corrects truncation errors by establishing an error equation. These high-order methods are applied to underlying solution strategies, i.e., monolithic, fixed-stress sequential, and undrained sequential methods. We observe that semi-discretized systems from spatial discretization have forms similar to those of index-1 differential algebraic equations (DAEs), causing order reduction against the two-pass method when it is used in conjunction with either the monolithic or sequential method. On the other hand, the SDC in conjunction with the monolithic method exhibits the desired second-order accuracy in poroelastic problems while increasing the order of accuracy for index-1 DAEs. However, the SDC in conjunction with either of the two sequential methods does not achieve the desired order of accuracy, and maintains first order because the flow equation for poroelasticity has an additional approximation associated with the volumetric strain rate term, which does not yield exactly the same forms as those of conventional DAEs. Thus, the monolithic SDC method can achieve higher-order accuracy, but may require higher computational costs because it involves solving matrix systems larger than those for the sequential methods.
02 PETROLEUM↗
Arbitrary Order Virtual Element Methods for High‐Order Phase‐Field Modeling of Dynamic Fracture
ABSTRACT Accurate modeling of fracture nucleation and propagation in brittle and ductile materials subjected to dynamic loading is important in predicting material damage and failure under extreme conditions. Phase‐field fracture models have garnered a lot of attention in recent years due to their success in representing damage and fracture processes in a wide class of materials and under a variety of loading conditions. Second‐order phase‐field fracture models are by far the most popular among researchers (and increasingly, among practitioners), but fourth‐order models have started to gain broader acceptance since their more recent introduction. The exact solution corresponding to these high‐order phase‐field fracture models has higher regularity. Thus, numerical solutions of the model equations can achieve improved accuracy and higher spatial convergence rates. In this work, we develop a virtual element framework for the high‐order phase‐field model of dynamic fracture. The virtual element method (VEM) can be regarded as a generalization of the classical finite element method. In addition to many other desirable characteristics, the VEM allows computing on polytopal meshes. Here, we use ‐conforming virtual elements and the generalized‐ time integration method for the momentum balance equation, and adopt ‐conforming virtual elements for the high‐order phase‐field equation. We verify our virtual element framework using classical quasi‐static benchmark problems and demonstrate its capabilities with the aid of numerical simulations of dynamic fracture in brittle materials.
42 ENGINEERING↗