Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reduced dimensions forms”

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.

66 records · Page 4

Determination of electron-hole pair creation energy in Cd 0.9 Zn 0.1 Te 0.98 Se 0.02 quaternary semiconductor for room-temperature gamma-ray detection

We report the first-time measurement of the electron-hole pair (ehp) creation energy (W ehp ) in novel Cd 0.9 Zn 0.1 Te 0.98 Se 0.02 (CZTS) quaternary semiconductor. CZTS in single crystalline form is poised to be the future of large-volume room-temperature gamma-ray detectors due to its excellent compositional homogeneity with highly reduced defects, high-Z (atomic number) constituents, wide bandgap (1.6 eV), and superior charge transport properties. Despite a great deal of study of the material and device properties since its inception, the W ehp in CZTS has not been measured experimentally. Accurate determination of W ehp is essential for calibration of the spectrometer and other theoretical calculations. In this study we have used an absolute calibration approach, which is based on an iterative approach that yields the Wehp as the best-fit parameter. Using a 241 Am alpha emitting radioisotope and a planar CZTS detector, the Wehp in CZTS was calculated to be 4.47 eV. The obtained value has been validated by accurately predicting the peak energy for gamma rays emitted by a 137 Cs source and read by a CZTS detector with different dimensions. The dependences of the calculated W ehp value on the detector dimensions, type of interaction, and effect of charge trapping are also discussed.

42 ENGINEERING↗

Extending Parsimonious Bayesian Inference

Parsimonious Bayesian inference is a theoretical framework for efficient data assimilation that seeks to balance increased consistency between predictions and training data against corresponding increases in model complexity. Within this framework, over-training is understood as optimization that encodes excessive information within model parameters while only achieving small improvements between predictions and training data. This project aims to develop practical methods of limiting excess model information during optimization. One key observation is that practical heuristics for parsimonious learning in high-dimensions must balance expressivity, i.e. the ability of the model to capture diverse predictions with only a few non-zero parameters, against discoverability, i.e. the ability to train the model with gradient-based optimization and drive parameters to low information states. As such, we developed logical activation functions that are able to adaptively approximate arbitrary truth tables that define Boolean logic operations within a probabilistic framework. These functions have demonstrated the ability to learn exclusive disjunction (XOR) and conditioned disjunction (if [condition] then [result_if_true] else [result_if_false]) within a single layer of a neural network. To efficiently exploit these activation functions to drive parsimonious learning required several other advances within the domain of variational inference. The most efficient form of complexity suppression is structured sparsification, driving most model parameters to zero while achieving the structural coherence among nonzeros needed for bandwidth reduction. Such models are not only far more efficient at suppressing information-theoretic complexity, they also reduce the other forms of complexity (computations, communication, storage, and the number of dependencies needed to evaluate predictions). Aiming to support enhanced sparsification, this project examined new approaches to high-dimensional variational inference that allow us to calibrate and control parameter uncertainty during optimization. By identifying which parameters can sustain sparsifying perturbations with little impact on prediction quality, we can develop better pruning strategies by framing them as approximate Bayesian inference. These advances also open paths to mitigate concerns with deploying advanced learning methods in resource-constrained environments, such as running models on power-limited or communication-limited devices.

97 MATHEMATICS AND COMPUTING↗

Noise and error analysis and optimization in particle-based kinetic plasma simulations

In this paper we analyze the noise in macro-particle methods used in plasma physics and fluid dynamics, leading to approaches for minimizing the total error, focusing on electrostatic models in one dimension. We begin by describing kernel density estimation for continuous values of the spatial variable x, expressing the kernel in a form in which its shape and width are represented separately. The covariance matrix of the noise in the density is computed, first for uniform true density. The bandwidth of the covariance matrix C(x,y) is related to the width of the kernel. A feature that stands out is the presence of constant negative terms in the elements of the covariance matrix both on and off-diagonal. These negative correlations are related to the fact that the total number of particles is fixed at each time step; they also lead to the property ∫C(x,y)dy = 0. We investigate the effect of these negative correlations on the electric field computed by Gauss's law, finding that the noise in the electric field is related to a process called the Ornstein-Uhlenbeck bridge, leading to a covariance matrix of the electric field with variance significantly reduced relative to that of a Brownian process. For non-constant density, p(x), still with continuous x, we analyze the total error in the density estimation and discuss it in terms of bias-variance optimization (BVO). For some characteristic length l, determined by the density and its second derivative, and kernel width h, having too few particles within h leads to too much variance; for h that is large relative to l, there is too much smoothing of the density. The optimum between these two limits is found by BVO. For kernels of the same width, it is shown that this optimum (minimum) is weakly sensitive to the kernel shape. Next, we repeat the analysis for x discretized on a grid. In this case the charge deposition rule is determined by a particle shape. An important property to be respected in the discrete system is the exact preservation of total charge on the grid; this property is necessary to ensure that the electric field is equal at both ends, consistent with periodic boundary conditions. We find that if the particle shapes satisfy a partition of unity property, the particle charge deposited on the grid is conserved exactly. Further, if the particle shape is expressed as the convolution of a kernel with another kernel that satisfies the partition of unity, then the particle shape obeys the partition of unity. This property holds for kernels of arbitrary width, including widths that are not integer multiples of the grid spacing. Furthermore, we show results relaxing the approximations used to do BVO optimization analytically, by doing numerical computations of the total error as a function of the kernel width, on a grid in x. The comparison between numerical and analytical results shows good agreement over a range of particle shapes. We discuss the practical implications of our results, including the criteria for design and implementation of computationally efficient particle shapes that take advantage of the developed theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Efficient computation of N -point correlation functions in D dimensions

We present efficient algorithms for computing the N-point correlation functions (NPCFs) of random fields in arbitrary D-dimensional homogeneous and isotropic spaces. Such statistics appear throughout the physical sciences and provide a natural tool to describe stochastic processes. Typically, algorithms for computing the NPCF components have $\mathscr O$(n N ) complexity (for a dataset containing n particles); their application is thus computationally infeasible unless N is small. By projecting the statistic onto a suitably defined angular basis, we show that the estimators can be written in a separable form, with complexity $\mathscr O$(n 2 ) or $\mathscr O$(n g log n g ) if evaluated using a Fast Fourier Transform on a grid of size n g . Our decomposition is built upon the D-dimensional hyperspherical harmonics; these form a complete basis on the (D – 1) sphere and are intrinsically related to angular momentum operators. Concatenation of (N – 1) such harmonics gives states of definite combined angular momentum, forming a natural separable basis for the NPCF. As N and D grow, the number of basis components quickly becomes large, providing a practical limitation to this (and all other) approaches: However, the dimensionality is greatly reduced in the presence of symmetries; for example, isotropic correlation functions require only states of zero combined angular momentum. We provide a Julia package implementing our estimators and show how they can be applied to a variety of scenarios within cosmology and fluid dynamics. The efficiency of such estimators will allow higher-order correlators to become a standard tool in the analysis of random fields.

97 MATHEMATICS AND COMPUTING↗

Bayesian Poroelastic Aquifer Characterization From InSAR Surface Deformation Data. 2. Quantifying the Uncertainty

Uncertainty quantification of groundwater (GW) aquifer parameters is critical for efficient management and sustainable extraction of GW resources. These uncertainties are introduced by the data, model, and prior information on the parameters. We develop a Bayesian inversion framework that uses Interferometric Synthetic Aperture Radar (InSAR) surface deformation data to infer the laterally heterogeneous permeability of a transient linear poroelastic model of a confined GW aquifer. The Bayesian solution of this inverse problem takes the form of a posterior probability density of the permeability. Exploring this posterior using classical Markov chain Monte Carlo (MCMC) methods is computationally prohibitive due to the large dimension of the discretized permeability field and the expense of solving the poroelastic forward problem. However, in many partial differential equation (PDE)-based Bayesian inversion problems, the data are only informative in a few directions in parameter space. For the poroelasticity problem, we prove this property theoretically for a one-dimensional problem and demonstrate it numerically for a three-dimensional aquifer model. Here we design a generalized preconditioned Crank-Nicolson (gpCN) MCMC method that exploits this intrinsic low dimensionality by using a low-rank-based Laplace approximation of the posterior as a proposal, which we build scalably. The feasibility of our approach is demonstrated through a real GW aquifer test in Nevada. The inherently two-dimensional nature of InSAR surface deformation data informs a sufficient number of modes of the permeability field to allow detection of major structures within the aquifer, significantly reducing the uncertainty in the pressure and the displacement quantities of interest.

54 ENVIRONMENTAL SCIENCES↗

Selective posttranslational inhibition of Ca V β 1 -associated voltage-dependent calcium channels with a functionalized nanobody

Ca 2+ influx through high-voltage-activated calcium channels (HVACCs) controls diverse cellular functions. A critical feature enabling a singular signal, Ca 2+ influx, to mediate disparate functions is diversity of HVACC pore-forming α1 and auxiliary Ca V β 1 –Ca V β 4 subunits. Selective Ca V α 1 blockers have enabled deciphering their unique physiological roles. By contrast, the capacity to post-translationally inhibit HVACCs based on Ca V β isoform is non-existent. Conventional gene knockout/shRNA approaches do not adequately address this deficit owing to subunit reshuffling and partially overlapping functions of Ca V β isoforms. Here, we identify a nanobody (nb.E8) that selectively binds Ca V β 1 SH3 domain and inhibits Ca V β 1 -associated HVACCs by reducing channel surface density, decreasing open probability, and speeding inactivation. Functionalizing nb.E8 with Nedd4L HECT domain yielded Chisel-1 which eliminated current through Ca V β 1 -reconstituted Ca V 1/Ca V 2 and native Ca V 1.1 channels in skeletal muscle, strongly suppressed depolarization-evoked Ca 2+ influx and excitation-transcription coupling in hippocampal neurons, but was inert against Ca V β 2 -associated Ca V 1.2 in cardiomyocytes. The results introduce an original method for probing distinctive functions of ion channel auxiliary subunit isoforms, reveal additional dimensions of Ca V β 1 signaling in neurons, and describe a genetically-encoded HVACC inhibitor with unique properties.

60 APPLIED LIFE SCIENCES↗

Criticality Analysis of FSV Spent Nuclear Fuel in the DOE Standard Canister

The U.S. Department of Energy (DOE) is responsible for managing over 300 types of spent nuclear fuel (SNF). To manage this large variety of fuel types, DOE plans to employ standardized canisters for the transportation, long-term storage, and eventual disposal of SNF. Idaho National Laboratory is currently supporting DOE’s SNF Packaging Demonstration Project, in which Fort Saint Vrain (FSV) fuel assemblies will be loaded into a DOE Standard Canister. This paper presents criticality calculations demonstrating that all four or five FSV fuel assemblies loaded into the DOE Standard Canister will remain subcritical in any expected normal or credible abnormal conditions. Previous criticality analyses were performed for one FSV fuel assembly and 12 Peach Bottom Core 2 fuel elements loaded into a DOE Standard Canister. This paper covers the criticality analysis performed for loading both four and five FSV fuel assemblies into a DOE Standard Canister. Various intact and degraded mode configurations were modeled in conducting the criticality calculations. This analysis encompassed three different configurations: (1) a single DOE Standard Canister loaded into a concrete storage overpack, (2) seven DOE Standard Canisters loaded into a concrete storage overpack, and (3) nine DOE Standard Canisters loaded into a concrete storage overpack. The overpack dimensions were varied for each of the three configurations, and transport, storage, and disposal scenarios were analyzed for each configuration. For transport scenarios, a pair of degradation cases was analyzed. In the first case, the fuel compacts became degraded and were removed from the fuel block, then deposited at the bottom of a horizontally placed canister, thereby simulating a drop event. The canister was considered to remain intact. In the second case, the spacing between horizontally placed canisters in a nine-canister overpack was reduced such that the canisters were piled on top of each other, simulating a drop event. For this case, no degradation of the canister internals or fuel was considered. For storage scenarios, the water moderator location in the system was varied to enable identification of the most reactive configurations. Dry and wet conditions were analyzed for the fuel materials, canister, and overpack. For disposal scenarios, two degradation cases were analyzed. In the first, the stainless-steel internals of the canister degraded to either hematite or goethite under both dry and wet conditions. In the second case, degraded FSV fuel formed a uranium-water slurry that filled the coolant/void holes. None of the cases presented exceeded the application specific upper subcritical limit.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Anatomy of Protein Electrospray Ionization Mass Spectra by Superconducting Tunnel Junction Mass and Energy Spectrometry

Cryogenic superconducting tunnel junction (STJ) detectors have the advantage of single-particle sensitivity, high quantum efficiency, low noise, and the ability to detect the time and relative impact energy of deposited ions. This makes them attractive for use in mass spectrometry (MS) and as a form of energy spectrometry. STJ cryodetectors have been coupled to time-offlight (TOF) mass spectrometers equipped with a matrix-assisted laser desorption ionization (MALDI) source and to an electrospray ionization (ESI) TOF mass spectrometer. Here, a lab-made linear ion trap (LIT) mass spectrometer system was coupled to an ESI source and a 16-channel Nb-STJ array with improved readout electronics. The goal was to investigate fundamentals of ESI-generated protein ions by further exploiting the advantage of resolving these ions in a third dimension of the relative energy deposited into the STJs. The proteins equine cytochrome c, bovine carbonic anhydrase, bovine serum albumin, and murine immunoglobulin G were studied using this ESI-LIT-STJ-MS instrument. Multiply charged monomers, multimers, and fragments from metastable ions were resolved from monomer peaks by differences in ion deposition energy even when these ions have the same mass-to-charge ratio as the corresponding monomer. The determination of fragment mass from metastable decomposition is accomplished without knowing the charge state of the fragment. The average charge state of the multimers is reduced with each addition of a protein which is a direct reflection of the surface area available for charging. Multiply charged in-source fragments have also been observed and distinguished in the mass spectrum of carbonic anhydrase by using the differences in the energy deposited in the STJs.

42 ENGINEERING↗

Adaptive experimental design for multi-fidelity surrogate modeling of multi-disciplinary systems

We present an adaptive algorithm for constructing surrogate models of multi-disciplinary systems composed of a set of coupled components. With this goal we introduce “coupling” variables with a priori unknown distributions that allow surrogates of each component to be built independently. Once built, the surrogates of the components are combined to form an integrated-surrogate that can be used to predict system-level quantities of interest at a fraction of the cost of the original model. The error in the integrated-surrogate is greedily minimized using an experimental design procedure that allocates the amount of training data, used to construct each component-surrogate, based on the contribution of those surrogates to the error of the integrated-surrogate. Here, the multi-fidelity procedure presented is a generalization of multi-index stochastic collocation that can leverage ensembles of models of varying cost and accuracy, for one or more components, to reduce the computational cost of constructing the integrated-surrogate. Extensive numerical results demonstrate that, for a fixed computational budget, our algorithm is able to produce surrogates that are orders of magnitude more accurate than methods that treat the integrated system as a black-box.

42 ENGINEERING↗

Atomically Precise Membranes for the Separation of Gases

Industrial separations require enormous amounts of energy, accounting for approximately half the industrial energy use and 10–15% of the total energy consumption. Distillation alone accounts for about half the energy demand for industrial separations. If these processes could be replaced by an energy-efficient membrane separation process, this energy demand could be reduced by 90%. However, although some membrane processes have made inroads into thermal distillation, for membrane-based separations to replace the distillation process to a far more significant and practical level, new membranes with higher robustness, selectivity, and flux still need be developed. In this membrane development program, we leveraged Temple University’s prior work in spiroligomers to develop robust membrane structures. These structures can be formed into atomically precise pores by controlling the chemical synthesis and the oligomer building blocks to first build precisely controlled macrocycles and then subsequently crosslinking these macrocycles to produce a membrane. Furthermore, through the highly controllable chemistry of our molecular building blocks, not only can we integrate pores with highly controllable and reproducible size and morphology, but we can also target internal functionalization. By using a range of scalable membrane synthesis approaches, combined with control of each and every pores internal chemistry and molecular conformation we can achieve membrane structures that can display the ultimate in high selectivity and permeance. The Mainstream-Temple University membranes can be designed to achieve ultra-selective separations based on the key factors of molecular size, shape, and functionality. In our approach to fabricate scalable atomically precise membranes, we used our molecular Lego nanostructures. We demonstrated an approach to create atomically precise pores within the membrane with pores that are the dimensions of the molecule we are trying to separate. Moreover, in addition to controlling the morphology of the membrane structure by controlling the size of every pore, we can also decorate every pore with precisely targeted and placed functional groups. These accurately placed functional groups can provide selective binding to molecules and provide enhanced selectivity via a facilitated transport mechanism. In Phase II, the Mainstream-Temple University team demonstrated the scalable synthesis of oligomers and the ability to control the pore internal and external functionality, or chemistry, to allow the fabrication of a thin-film membrane. The initial steps of the spiroligomer synthesis to produce the two key bis-amino enantiomer building blocks were scaled from the gram scale to the kilogram scale, obtaining 25 kg of the materials by transferring the procedure to a toll manufacturer. Finally, we successfully optimized the fabrication of these macrocycles into robust membranes. We successfully transitioned from a laboratory based, hard-to-scale Langmuir trough synthesis to a highly scalable, roll-to-roll applicable, interfacial polymerization process. During this Phase II program, we established a platform of atomically precise membranes where our highly controllable, atomically precise macrocycles served as a scaffold with precisely and uniformly controlled pores. Furthermore, this layer can be tailored to accommodate a diverse range of functional groups, which both further controlled the pore size to enhance the sieving effect as well as imparting precisely controlled targeted selectivity through biomimetic molecular interactions. In this Phase II, we established and scaled up a platform approach to both tune the pore size and chemistry as well as scale it to produce membranes that can be applied to a wide range of industries. In future development, the project team expects to scale up both the macrocycle building block production and interfacial polymerization process to produce the atomically precise membrane with targeted pore sizes and pore chemistries.

36 MATERIALS SCIENCE↗

Tandem particle-slurry batch reactors for solar water splitting (Final Scientific/Technical Report)

Economically, particle slurry reactors are projected to be one of the most promising technologies for solar photoelectrochemical hydrogen production, according to a 2009 techno-economic analysis commissioned by the US DOE and performed by Directed Technologies, Inc. The Fuel Cell Technologies Office’s Multi-Year Research, Development and Demonstration (MYRD&D) goals and targets are to reduce the cost of H 2 produced from renewable sources at the plant gate (i.e. not including delivery, dispensing, or storage) to < $2.00/gge, equivalent to ~$2.00/kg H 2 . Research results from our techno-economic modeling research suggest that this target could be met using particle slurry reactors assuming STH efficiencies in the range of 5 – 10%, materials lifetimes of < 1 year, and nanoparticles that cost up to 20 times more than projected costs of TiO 2 -coated Fe 2 O 3 nanoparticles. Although most large worldwide research efforts directed at solar photoelectrochemical hydrogen production focus on wafer-based designs, the projected lower cost for a particle slurry reactor at these disparate projected STH efficiencies clearly suggests that particle slurry reactors could be a scalable and deployable technology, assuming several challenges are overcome. These major technological challenges include the demonstration of a vertically-stacked-vessel architecture that is capable of operating sustainably while mostly relying on diffusion and natural convection to mix the redox shuttles between the vessels, and the demonstration that photocatalyst particles can operate at an overall 1% STH efficiency or larger when incorporated into this two-vessel design. Our research adds to the understanding of photocatalytic reactors for solar water splitting through numerical modeling results and experimental results. Numerical models were developed to simulate relevant device physics including particle and reactor dimensions which affect optical, transport, and rheological properties, electrocatalytic and photovoltaic properties of particles at various temperatures, and properties of redox shuttles and separators. Moreover, theoretical maximum solar-to-hydrogen efficiencies for ensembles of particles like in photocatalyst reactors were modeled and simulated and shown to equal or exceed those of photoelectrochemical designs under most scenarios. These results help determine constraints on the reactor that will enable more optimal designs for future prototypes. In parallel, experiments were performed to identify the most effective redox shuttles and to empirically validate the numerical models and simulations. Toward the latter, state-of-the-art light-absorbing particles and electrocatalysts were synthesized and characterized physically and photoelectrochemically for water electrolysis and redox chemistry with redox shuttles in the form factor of mesoporous electrodes and free-floating particles. The most promising materials candidates were used in a suspension reactor to evaluate performance toward photocatalytic H 2 production and results from the two measurements were compared. Predominantly, state-of-the-art cocatalyst-modified Rh-doped SrTiO 3 and BiVO 4 particles were further characterized to assess for their ability to perform visible-light-driven H 2 and O 2 evolution, respectively, and results were similar to those reported for the state-of-the-art in the peer-reviewed literature. Outcomes from this work inform the public of the effectiveness and promise of solar photocatalytic water splitting for clean and renewable hydrogen production. This work may also help increase research interest and funding for photocatalysis projects, which will accelerate development of a technology that will benefit the public by generating fuel while emitting few greenhouse gases and pollutants.

08 HYDROGEN↗

Development, Implementation, and Verification of Partially-Ionized Collisional Multifluid Plasma Models in Drekar

We consider the development of multifluid models for partially ionized multispecies plasmas. The models are composed of a standard set of five-moment fluid equations for each species plus a description of electromagnetics. The most general model considered utilizes a full set of fluid equations for each charge state of each atomic species, plus a set of fluid equations for electrons. The fluid equations are coupled through source terms describing electromagnetic coupling, ionization, recombination, charge exchange, and elastic scattering collisions in the low-density coronal limit. The form of each of these source terms is described in detail, and references for required rate coefficients are identified for a diverse range of atomic species. Initial efforts have been made to extend these models to incorporate some higher-density collisional effects, including ionization potential depression and three- body recombination. Some reductions of the general multifluid model are considered. First, a reduced multifluid model is derived which averages over all of the charge states (including neutrals) of each atomic species in the general multifluid model. The resulting model maintains full consistency with the general multifluid model from which it is derived by leveraging a quasi-steady-state collisional ionization equilibrium assumption to recover the ionization fractions required to make use of the general collision models. Further reductions are briefly considered to derive certain components of a single-fluid magnetohydrodynamics (MHD) model. In this case, a generalized Ohm's law is obtained, and the standard MHD resistivity is expressed in terms of the collisional models used in the general multifluid model. A number of numerical considerations required to obtain robust implementations of these multifluid models are discussed. First, an algebraic flux correction (AFC) stabilization approach for a continuous Galerkin finite element discretization of the multifluid system is described in which the characteristic speeds used in the stabilization of the fluid systems are synchronized across all species in the model. It is demonstrated that this synchronization is crucial in order to obtain a robust discretization of the multifluid system. Additionally, several different formulations are considered for describing the electromagnetics portion of the multifluid system using nodal continuous Galerkin finite element discretizations. The formulations considered include a parabolic divergence cleaning method and an implicit projection method for the traditional curl formulation of Maxwell's equations, a purely- hyperbolic potential-based formulation of Maxwell's equations, and a mixed hyperbolic-elliptic potential-based formulation of Maxwell's equations. Some advantages and disadvantages of each formulation are explored to compare solution robustness and the ease of use of each formulation. Numerical results are presented to demonstrate the accuracy and robustness of various components of our implementation. Analytic solutions for a spatially homogeneous damped plasma oscillation are derived in order to verify the implementation of the source terms for electromagnetic coupling and elastic collisions between fluid species. Ionization balance as a function of electron temperature is evaluated for several atomic species of interest by comparing to steady-state calculations using various sets of ionization and recombination rate coefficients. Several test problems in one and two spatial dimensions are used to demonstrate the accuracy and robustness of the discretization and stabilization approach for the fluid components of the multifluid system. This includes standard test problems for electrostatic and electromagnetic shock tubes in the two-fluid and ideal shock-MHD limits, a cylindrical diocotron instability, and the GEM challenge magnetic reconnection problem. A one-dimensional simplified prototype of an argon gas puff configuration as deployed on Sandia's Z-machine is used as a demonstration to exercise the full range of capabilities associated with the general multifluid model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗