Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density matrix methods”

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 541 records · Page 30

High Velocity Impact Performance of a Dual Layer Thermal Protection System for the Mars Sample Return Earth Entry Vehicle

The Mars Sample Return (MSR) Earth Entry Vehicle (EEV) is currently planned on being released from its Micro-Meteorite/Orbital Debri (MM/OD) shielding housing about two days before the Earth entry phase. This leaves the EEV exposed to incoming MM/OD impacts, potentially damaging the heat shield and compromising its Entry, Decent, and Landing (EDL) integrity. Currently, two materials are proposed to comprise the MSR-EEV heat shield, a dual layer material Heat-shield for Extreme Entry Environment Technology (HEEET) and Phenolic Impregnated Carbon Ablator (PICA). PICA has been well characterized for OD class impacts, ~7km/s high mass impacts, from previous testing done in the Orion program, but hasn't under-gone extensive MM impact testing. HEEET is a rela-tively new material with minimal prior testing in re-gards to High Velocity Impacts (HVI). In order to in-form selectability of a material, it is crucial to under-stand the material performance when faced with an HVI, directly affecting mission success probability.The current measure of a Thermal Protection Sys-tem's (TPS) performance against an HVI is evaluating a thermally sized material against its derived Balistic Limit Equation (BLE). A BLE is generated empirically from multiple shots of HVI testing, and is used as a first order method in evaluating TPS's performance against the expected MM/OD environment. This method has proved useful for previous uniform densi-ty TPS materials, but has never been validated against a dual layer recession material such as HEEET. Testing in the MSR program for FY19 has a re-quirement to assess the effects of a dual layer TPS by testing various thicknesses of HEEET's Recession Layer, seen in Figure 1. From this data, BLE's will be derived for the individual thickness ratio samples, as well as the material as a whole to evaluate if a heritage form of the BLE can capture the complex physics associated with a dual layer system. Crater morpholo-gy will also be assessed with post-test Non-Destructive Evaluation (NDE) methods such as CT scanning to visualize if a BLE can well predict the associated pene-tration depths, since a BLE assumes full disinigration of the impacting particle and is generally only used to size a spherical crater ? disregarding any shrapenel effects from a high density impactor.To inform this analysis, expected MM environ-ments from the Meteoroid Engineering Model (MEM) and analytical equations for mass flux of incoming MM were evaluated against the notional MSR-EEV trajectory [1]. Using those dispersions, a monte-carlo was run to determine the most probable particle pa-rameters, as well as the riskiest in terms of full bondline penetration. From these probabilities, a test matrix was designed to test against bounding cases for the various parameters of the BLE: projectile density, projectile mass, projectile velocity, and the impact angle.This presentation will discuss the performance of the dual layer TPS material HEEET against a wide range of impact kinetic energies and densitites, as well as the comparison of HEEET to PICA in terms of MM/OD performance and selectability criteria.

Libben, B. J.↗

Probabilistic Modeling of a Three-Stage Human Landing System Architecture

Space Policy Directive-1 has led to NASA partnerships with commercial entities on procurement which includes the development of the Human Landing System (HLS) [1]. With the goal of delivering human crew to the lunar surface by 2024, system uncertainties become an important obstacle to the maturation of multiple new, driving technologies and mission concepts of the HLS program. As unmitigated uncertainties have previously led to failed development programs, these risks and their impacts must be understood and handled to ensure program success [2]. Sources of uncertainty include novel engine designs and configurations, increased reliance on cryogenic fluid management(CFM), and refueling technologies—which propagate as high-level performance metrics such as overall propellant mass and engine performance. Also, the occurrence of operational uncertainties—e.g. launch conditions or need to abort during the mission—can cause cascading effects on the rest of the mission that are difficult to definitively quantify, and are outside the scope of control. These concrete examples and other occurrences can be categorized as either epistemic or aleatory uncertainties.Epistemic uncertainty arises due to a lack of knowledge and can be alleviated with design and program maturation. Aleatory uncertainty is due to the inherent randomness of the system and cannot be directly reduced, unlike epistemic uncertainty. Robust design and probabilistic methods can compensate for aleatory effects. A taxonomy of uncertainty is referred to for this work [3]. In this paper, a probabilistic methodology to handle uncertainties has been demonstrated on a three-element HLS concept [1, 4], which allows tracking of current best estimates of the concept and assessment of concept design robustness against uncertainties. A sample case has been completed for this abstract, and an expansion on the methodology will be included in the final paper. This methodology has two key parts: first, the creation of a dynamic architecture model of a three-element HLS concept; and second, its use with surrogate modeling and range estimating techniques to capture and propagate uncertainties. This abstract will cover the basics of the approach used, and further details and justifications will be in the final paper.The mission profile associated with this three-element concept (Fig 1) was modeled as a set of mission events that facilitated mass changes, idles, or spacecraft maneuvers. The mission profile scope starts with each element’s NRHO orbit insertion and aggregation and ends at post-sortie rendezvous with Orion. More detail on the mission profile will be in the final paper. The DYnamic Rocket EQuation Tool (DYREQT), a space systems synthesis and sizing framework used by NASA, was used as the physics framework to model the HLS architecture for applying the probabilistic methodology [5, 6]. Specifically, a parametric representation of the lander, ascent, and transfer elements and the mission profile of each element was established, with vehicle and mission parameters available as inputs to allow for a dynamic model. Each vehicle stage was modeled with high-level performance metrics, using Isp and propellant mass fraction (PMF) to remain parametric. For the probabilistic analysis, uncertainties of interest within the HLS concept were enumerated and represented as parameters within the DYREQT model as inputs for vehicle stages or mission profile events. These parameters were frozen at their nominal values for the purposes of baselining architecture performance and sizing the vehicle appropriately based on reference documentation [1]. Range estimating—a probabilistic method that combines Monte Carlo sampling, focus on critical parameters, and heuristics to assess risk and opportunities—is traditionally used with Mass Equipment Lists (MELs), but has been adapted with operational parameters as well as vehicle parameters in theDYREQT model to capture mission uncertainty alongside vehicle uncertainty [7, 3]. This method was selected due to its application and insight on a system from a bottom-up perspective, independence from historical rules of thumb, and ability to generate sensitivities based on design decisions and uncertainties. As a sample case for the abstract, the boiloff rates of the vehicle elements and the loiter times during the mission (simulating launch time variations and changing window of opportunities) were used with range estimating to provide preliminary results. To perform the range estimation portion of this methodology (depicted in Fig. 3, further details in final paper), the DYREQT model was sampled using a Design of Experiments (DoE) to efficiently explore the architecture design space with respect to the sample set of uncertainty parameters; 5,000 cases via Latin Hypercube Sampling were computed on the DYREQT architecture model. Then, the results were used to create surrogate models, multivariate regressions that can visualize hypercube trends in the design space, of the architecture with respect to the uncertainty parameters. Range estimating was applied to the surrogates instead of the actual models, which saves computational expense due to the bulk of cases needed for the Monte Carlo simulation as part of range estimating. Uncertainty parameters were sampled independently from triangular distributions using the DoE ranges as ‘min’ and ‘max’, and the nominal value as ‘most likely’. Based engineering intuition, some uncertainty parameters are correlated—e.g. if the main propellant has a high boil-off rate, the oxidizer should follow suit as both are related to CFM technology.While a Monte Carlo simulation samples all inputs as independent, the results would show model correlations; thus, it is efficient to sample the inputs as correlated. Using a correlation matrix constructed for the uncertainty parameters, previously independent samples were transformed to perform a Correlated Monte Carlo. A table for the DoE ranges and probability distribution parameters is shown in Table 1, and more details on Correlated Monte Carlo Simulations will be discussed in the final paper. The model’s resulting DoE showed that multivariate polynomial equations fit via least squares method captured its behavior accurately for the sample case. For the Correlated Monte Carlo Simulation, a positive correlation between fuel and oxidizer boiloff rates was used as a demonstration. 10,000 cases were computed with the surrogates and the launched masses for each vehicle element was collated. The results can be displayed in a probability density function (PDF), showing the impact of the uncertainty parameters chosen. Integrating the PDFs will yield a cumulative distribution function (CDF) that shows the cumulative probability of a given value on the x-axis. For the sample case, the elements’ launch mass margin was calculated and represented in as CDFs, as a demonstrated representation of figures of merit for the HLS concept. For the lander and ascent elements, the NRHO mass insertion limit is 16t; the transfer element has a limit of 30t [1]. It can be seen with Figure 2 that this probabilistic methodology can provide insight into mass margin with respect to the uncertainties being modeled. Currently, the results show that the lander (descent) vehicle element has the most restrictive design space; it is the only element to show a 10% probability of negative margin. Further analysis on the Monte Carlo results will show sensitivities for driving constraints and parameters for architecture feasibility, which can lead to establishing potential mission rules.The combination of range estimating with a parametric architecture model for HLS demonstrated the capability of this probabilistic methodology in a sample case. As the HLS development progresses, this methodology has the potential for keeping current best estimates of architecture performance for awarded concepts due to the flexibility in DYREQT’s modeling framework and its parametric nature. Concept maturation and increased epistemic knowledge can be injected into the model probabilistic modeling, and thus continue to track probability of mission success.

Stephanie Y Zhu↗

Constraining Hamiltonians from chiral effective field theory with neutron-star data

Multi-messenger observations of neutron stars (NSs) and their mergers have placed strong constraints on the dense-matter equation of state (EOS). The EOS, in turn, depends on microscopic nuclear interactions that are described by nuclear Hamiltonians. These Hamiltonians are commonly derived within chiral effective field theory (EFT). Ideally, multi-messenger observations of NSs could be used to directly inform our understanding of EFT interactions, but such a direct inference necessitates millions of model evaluations. This is computationally prohibitive because each evaluation requires us to calculate the EOS from a Hamiltonian by solving the quantum many-body problem with methods such as auxiliary-field diffusion Monte Carlo (AFDMC), which provides very accurate and precise solutions but at a significant computational cost. Additionally, we need to solve the stellar structure equations for each EOS which further slows down each model evaluation by a few seconds. In this work, we combine emulators for AFDMC calculations of neutron matter, built using parametric matrix models, and for the stellar structure equations, built using multilayer perceptron neural networks, with the PyCBC data-analysis framework to enable a direct inference of coupling constants in an EFT Hamiltonian using multi-messenger observations of NSs. We find that astrophysical data can provide informative constraints on two-nucleon couplings despite the high densities probed in NS interiors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Optimal control of systems with capacity: Related noises

In the ordinary theory of optimal control (LQR and Kalman filter), the variances of the actuators and the sensors are assumed to be known (not related to the capacities of the devices). This assumption is not true in practice. Generally, a device with greater capacity to exert actuating forces and a sensor capable of sensing greater sensing range will generate noise of greater power spectral density. When the ordinary theory of optimal control is used to estimate the errors of the outputs in such cases it will lead to faulty results, because the capacities of such devices are unknown before the system is designed. The performance of the system designed by the ordinary theory will not be optimal as the variances of the sensors and the actuators are neither known nor constant. The interaction between the control system and structure could be serious because the ordinary method will lead to greater feedback (Kalman gain) matrices. Methods which can optimize the performance of systems when noises of the actuators and the sensors are related to their capacities are developed. These methods will result in smaller feedback (Kalman gain) matrix.

Ruan, Milfang↗

Theoretical Investigation of OCN(-) Charge Transfer Complexes in Condensed Phase Media: Spectroscopic Properties in Amorphous Ice

Density functional theory (DFT) calculations of cyanate (OCN(-)) charge-transfer complexes were performed to model the "XCN" feature observed in interstellar icy grain mantles. OCN(-) charge-transfer complexes were formed from precursor combinations of HNCO or HOCN with either NH3 or H2O. Three different solvation strategies for realistically modeling the ice matrix environment were explored, including (1) continuum solvation, (2) pure DFT cluster calculations, and (3) an ONIOM DFT/PM3 cluster calculation. The model complexes were evaluated by their ability to reproduce seven spectroscopic measurements associated with XCN: the band origin of the OCN(-) asymmetric stretching mode, shifts in that frequency due to isotopic substitutions of C, N, O, and H, plus two weak features. The continuum solvent field method produced results consistent with some of the experimental data but failed to account for other behavior due to its limited capacity to describe molecular interactions with solvent. DFT cluster calculations successfully reproduced the available spectroscopic measurements very well. In particular, the deuterium shift showed excellent agreement in complexes where OCN(-) was fully solvated. Detailed studies of representative complexes including from two to twelve water molecules allowed the exploration of various possible solvation structures and provided insights into solvation trends. Moreover, complexes arising from cyanic or isocyanic acid in pure water suggested an alternative mechanism for the formation of OCN(-) charge-transfer complexes without the need for a strong base such as NH3 to be present. An extended ONIOM (B3LYP/PM3) cluster calculation was also performed to assess the impact of a more realistic environment on HNCO dissociation in pure water.

Park, Jin-Young↗

Hadron spectroscopy of pseudoscalar and vector meson photoproduction using linearly polarised photons at CLAS

The field of hadron spectroscopy is informed by experiments which study excited states of the nucleon by detecting the decay products following excitation with a hadronic, leptonic or electromagnetic probe. Experiments using a polarised beam of photons provide an important contribution to the world data. This thesis describes the analysis of data from one such experiment, the CEBAF Large Acceptance Spectrometer at Jefferson Laboratory, Virginia, USA. In the experiment, an incident beam of electrons passes through a diamond radiator producing a linearly polarised photon beam of energy from 1.1 to 2.1 GeV on a liquid hydrogen target. The analysis method and results for two reaction channels are presented. In the first analysis, the polarisation observables {?, P, T, Ox, Oz} for the reaction ??? p ? K0 S ?+ are extracted simultaneously using likelihood sampling with Markov-Chain Monte-Carlo. The results are presented as angular distributions (cos ?K0S) for bins in photon energy E? for a total of 21 bins. The values extracted for T , Ox and Oz are a first measurement of these quantities for this reaction and the beam asymmetry (?) measurement supplements the one previous measurement for this quantity and extends the energy range over which it has been extracted. The measurement of the recoil polarisation, P, is consistent with previous measurements. In a second analysis, the spin density matrix elements (SDMEs) for ? p ? p ? are extracted using similar analysis techniques. The results are presented as angular distributions (cos ?? ) for bins in photon energy E? for a total of 24 bins. This work is the first to extract all nine SDMEs for the full angular range of the ? meson production. A comparison is made to the predictions of the Vector Meson Dominance (VMD) model and the results provide evidence for contributing processes other than VMD, particularly at more backward angles of the ? meson production. The data also indicate non-helicity conserving process in both the s-channel and the t-channel. The results from the first analysis have been passed to the Juelich-Bonn theory group and the effect of including the new data within their dynamical coupled channel model is described. Several nucleon resonances are affected and the impact on the pole positions, widths and photon decay amplitudes of the affected resonances is described.

Clark, Louise↗

Performance of multi level error correction in binary holographic memory

At the Optical Computing Lab in the Jet Propulsion Laboratory (JPL) a binary holographic data storage system was designed and tested with methods of recording and retrieving the binary information. Levels of error correction were introduced to the system including pixel averaging, thresholding, and parity checks. Errors were artificially introduced into the binary holographic data storage system and were monitored as a function of the defect area fraction, which showed a strong influence on data integrity. Average area fractions exceeding one quarter of the bit area caused unrecoverable errors. Efficient use of the available data density was discussed. .

data storage↗

A new color correlation method applied to XRF Al/Si ratios and other lunar remote sensing data

Orbital X-ray fluorescence Al/Si intensity ratios, corrected for variations in solar activity, are correlated with normal albedo, elevation measurements from laser altimetry data, and gamma ray data in the 2.75-8.60 MeV range. Each of these data sets is placed into a digital array consisting of 1/4 deg latitude by 1/4 deg longitude pixels. Information relative to the correlation of Al/Si ratios with each of the other data sets is presented in the following forms: (1) histograms are given for each data set to show the frequency distribution within the areas of common coverage; (2) density plots are produced from the plot of a two-dimensional array consisting of the Al/Si ratio vs the other parameter value for each pixel; and (3) color correlation maps are produced by placing the two-dimensional array into a 3 x 3 matrix consisting of nine equal subarrays containing an equal number of data points.

Clark, P. E.↗

Direct Feed High-Level Waste APPS Model Glass Testing (DFHLW APPS) Matrix

This report summarizes the data collected during the batching and melting of the Direct Feed High-Level Waste APPS Model Glass Matrix (DFHLW APPS) to serve as a quality-assured validation of the Aspen Process Performance Simulation (APPS) formulation method. Of 15 glasses tested, 12 satisfied all target property constraints. Two glasses, APPS-05 and -06, formed nepheline on canister centerline cooling heat-treatment and failed the Product Consistency Test response limits. Glass APPS-07-2 formed unacceptably high concentrations of crystals (primarily Na3Nd(PO4)2) when heat treated at 950 °C. All other glasses were found to be satisfactory. The measured property values were compared to predicted values from a set of current models. In many cases the current models were found to be inadequate for design of DFHLW glasses. These models are being adjusted to correct for mispredictions. Other models, e.g., density, toxicity characteristic leaching procedure, and sulfur solubility, are adequate for formulation of DFHLW glasses.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Efficient Langevin dynamics for “noisy” forces

Efficient Boltzmann-sampling using first-principles methods is challenging for extended systems due to the steep scaling of electronic structure methods with the system size. Additionally, stochastic approaches provide a gentler system-size dependency at the cost of introducing “noisy” forces, which could limit the efficiency of the sampling. When the forces are deterministic, the first-order Langevin dynamics (FOLD) offers efficient sampling by combining a well-chosen preconditioning matrix S with a time-step-bias-mitigating propagator [G. Mazzola and S. Sorella, Phys. Rev. Lett. 118, 015703 (2017)]. However, when forces are noisy, S is set equal to the force-covariance matrix, a procedure that severely limits the efficiency and the stability of the sampling. Here, we develop a new, general, optimal, and stable sampling approach for FOLD under noisy forces. We apply it for silicon nanocrystals treated with stochastic density functional theory and show efficiency improvements by an order-of-magnitude.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Too big, too small, or just right? A benchmark assessment of density functional theory for predicting the spatial extent of the electron density of small chemical systems

Multipole moments are the first-order responses of the energy to spatial derivatives of the electric field strength. The quality of density functional theory prediction of molecular multipole moments thus characterizes errors in modeling the electron density itself, as well as the performance in describing molecules interacting with external electric fields. However, only the lowest non-zero moment is translationally invariant, making the higher-order moments origin-dependent. Therefore, instead of using the 3 × 3 quadrupole moment matrix, we utilize the translationally invariant 3 × 3 matrix of second cumulants (or spatial variances) of the electron density as the quantity of interest (denoted by K). The principal components of K are the square of the spatial extent of the electron density along each axis. A benchmark dataset of the principal components of K for 100 small molecules at the coupled cluster singles and doubles with perturbative triples at the complete basis set limit is developed, resulting in 213 independent K components. The performance of 47 popular and recent density functionals is assessed against this Var213 dataset. Several functionals, especially double hybrids, and also SCAN and SCAN0 predict reliable second cumulants, although some modern, empirically parameterized functionals yield more disappointing performance. The H, Li, and Be atoms, in particular, are challenging for nearly all methods, indicating that future functional development could benefit from the inclusion of their density information in training or testing protocols.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An approximate Kalman filter for ocean data assimilation: An example with an idealized Gulf Stream model

A practical method of data assimilation for use with large, nonlinear, ocean general circulation models is explored. A Kalman filter based on approximation of the state error covariance matrix is presented, employing a reduction of the effective model dimension, the error's asymptotic steady state limit, and a time-invariant linearization of the dynamic model for the error integration. The approximations lead to dramatic computational savings in applying estimation theory to large complex systems. We examine the utility of the approximate filter in assimilating different measurement types using a twin experiment of an idealized Gulf Stream. A nonlinear primitive equation model of an unstable east-west jet is studied with a state dimension exceeding 170,000 elements. Assimilation of various pseudomeasurements are examined, including velocity, density, and volume transport at localized arrays and realistic distributions of satellite altimetry and acoustic tomography observations. Results are compared in terms of their effects on the accuracies of the estimation. The approximate filter is shown to outperform an empirical nudging scheme used in a previous study. The examples demonstrate that useful approximate estimation errors can be computed in a practical manner for general circulation models.

Fukumori, Ichiro↗

DMRG on Top of Plane-Wave Kohn–Sham Orbitals: A Case Study of Defected Boron Nitride

In this paper, we analyze the numerical aspects of the inherent multi reference density matrix renormalization group (DMRG) calculations on top of the periodic Kohn-Sham density functional theory using the complete active space approach. The potential of the framework is illustrated by studying hexagonal boron nitride nanoflakes embedding a charged single boron vacancy point defect by revealing a vertical energy spectrum with a prominent multireference character. We investigate the consistency of the DMRG energy spectrum from the perspective of sample size, basis size, and active space selection protocol. Results obtained from standard quantum chemical atom-centered basis calculations and plane-wave based counterparts show excellent agreement. Furthermore, we also discuss the spectrum of the periodic sheet which is in good agreement with extrapolated data of finite clusters. These results pave the way toward applying the DMRG method in extended correlated solid-state systems, such as point defect qubit in wide band gap semiconductors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Methodology for Generating Covariance Data of Thermal Neutron Scattering Cross Sections

This paper details and implements a framework for evaluating thermal neutron scattering cross sections that provide S(α,β) data and covariance data for hydrogen in light water. This methodology involves perturbing model parameters of molecular dynamics potentials and fitting the simulation results to experimental data. The framework is general and can be applied to any material or simulation method. The fit is made using the Unified Monte Carlo method to experimentally measure double-differential scattering cross sections of light water at the Spallation Neutron Source at Oak Ridge National Laboratory. Mean values and covariance data were generated for model parameters, phonon density of states, double-differential cross sections, and total scattering cross sections. These posterior parameter values were very similar to their prior values with a maximum relative error of 0.54%. This falls within in the Unified Monte Carlo–calculated uncertainties on the order of 2.7%. Additionally, posterior double-differential cross sections agree favorably with ENDF/B-VIII.0 cross sections. The new thermal scattering law was tested by comparing it against benchmarks from the International Criticality Safety Benchmark Evaluation Project Handbook, which showed a slight improvement over the ENDF/B-VIII.0 library. Additionally, the covariance matrix of the phonon density of states was validated to confirm that the spread of k eff from the density of states used to generate the covariance matrix was similar to the spread of k eff from the density of states of the sampled covariance matrix.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Chlorophylls as primary electron acceptors in reaction centers

Photosynthetic Reaction Centers (RCs) can be considered blueprints for highly efficient energy transfer. Embedded with an array of cofactors, including (bacterio)chlorophyll ((B)Chl) and (B)pheophytin ((B)Pheo) molecules, RCs function with a high quantum yield that spans a wide spectral range. Understanding the principles that underlie their function can influence the design of the next generation of artificial photosynthetic devices. We are particularly interested in the factors that influence the early stages of light-driven charge separation in RCs. With the recent publication of several highly anticipated RC structures and advanced computational methods available, it is possible to probe both the geometric and electronic structures of an array of RCs. In this chapter, we review the electronic and geometric structures of the (B)Chl and (B)Pheo primary electron acceptors from fiveRCs, comprising both Type I and Type II RCs and representing both heterodimeric and homodimeric systems. We showcase the dimeric A0•- state of Type I RCs, whereby the unpaired electron is delocalized, to various extents, over two (B)Chl molecules (B)Chl2 and (B)Chl3. This delocalization is controlled by several factors, including the structure of the (B)Chls, interactions with the surrounding protein matrix, and the orientation and distances of the cofactors themselves. In contrast, the primary acceptors of Type II RCs are entirely monomeric, with electron density residing solely on the (B)Pheo. We compare the natural design of the primary acceptors of the Type I and Type II RCs from both an evolutionary and application based perspective.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Bio-inspired redox mediated electrolyte for high performance flexible supercapacitor applications over broad temperature domain

For practical supercapacitors, recent advances in electrolytes have diminished long-standing problems such chemical stability, low ionic conductivity and limited electrochemical performance. Herewith, we communicate a facile method to produce anhydrous redox-active Bio-inspired electrolytes for flexible supercapacitors with application over a broad temperature range. The glycerol (Gly) gel electrolyte consists of various fractions of potassium hydroxide (KOH) and redox active unit, i.e., ammonium molybdate, (NH 4 ) 2 MoO 4 (abbreviated as Mo), in the anhydrous state. The doped glycerol matrix acts as a framework that provides a stable gel electrolyte, forming hierarchical pathway for ion diffusion. The developed gel network with Gly-KOH-MoX (X = 3, 5, 10, and 15%, w/w) provides an effective route to carry the ions onto the pores of the electrodes. The supercapacitor with Gly-KOH5-Mo10 gel provides a specific capacitance of 328 F g –1 , and the same device achieved an energy density of 45.6 Wh kg –1 at a power density of 497 W kg –1 . In addition, the supercapacitor maintained extraordinary cyclic robustness and wide temperature tolerance. Thus, the low cost and environmentally benign, cyclic durability of Gly-KOH-Mo makes these gel electrolytes a potential candidate for use in various energy storage devices.

25 ENERGY STORAGE↗

An Approximate Kalman Filter for Ocean Data Assimilation; An Example with an Idealized Gulf Stream Model

A practical method of data assimilation for use with large, nonlinear, ocean general circulation models is explored. A Kalman filter based on approximations of the state error covariance matrix is presented, employing a reduction of the effective model dimension, the error's asymptotic steady-state limit, and a time-invariant linearization of the dynamic model for the error integration. The approximations lead to dramatic computational savings in applying estimation theory to large complex systems. We examine the utility of the approximate filter in assimilating different measurement types using a twin experiment of an idealized Gulf Stream. A nonlinear primitive equation model of an unstable east-west jet is studied with a state dimension exceeding 170,000 elements. Assimilation of various pseudo measurements is examined, including velocity, density, and volume transport at localized arrays, and realistic distributions of satellite altimetry and acoustic tomography observations. Results are compared in terms of their effects on the accuracies of the estimation. The approximate filter is shown to outperform a previous study that used an empirical nudging scheme. The examples demonstrate that useful approximate estimation errors can be computed in a practical manner for general circulation models.

['Kalman↗

Evaluation of a Geotechnical Compaction Model to Estimate Sediment Deformation

When surface ruts from vehicle activity occur above archeological deposits, there is a potential for an impact to the archeological record to occur. While impacts to the soil surface are relatively straightforward to describe, sub-surface disturbances are more difficult to quantify. In particular, disturbances to buried archeological materials can be difficult to characterize. Materials can be damaged by surface disturbances, but they can also be impacted even if they stay intact. Interpretation of intact materials relocated within the sediment column can be problematic. Traditional characterization methods (e.g. excavation) can also disturb archaeological materials. Here, we present a geotechnical compaction model that can be used to assist with sub-surface characterization under vehicle ruts formed during off-road travel. The model is based on a Boussinesq analysis of stress increase below a surface load. The change in bulk density resulting from compaction is assumed to be proportional to the stress increase. The model results were compared against pilot scale tests conducted with homogeneous sediment and field scale in-situ tests. The results indicate that the model generally agreed with measurements to within 20%, although only a single soil type was evaluated in these studies, with a narrow range of initial bulk densities and rut depths. While application of the model across a broader matrix of soil types, bulk densities and rut depths would increase confidence in the model, it appears that the model may be suitable for use in archeological field investigations. Potential uses include both compliance and research based field applications.

Fritz, Brad G.↗