Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Parameter Identifiability”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,009 records · Page 56

Microwave-assisted conversion of methane over H-(Fe)-ZSM-5: Evidence for formation of hot metal sites

We report microwave-assisted catalysis offers great promise as an “intensified” technology for chemical processing. The present study investigates microwave-assisted direct conversion of methane with focus on the design and evaluation of microwave-sensitive H-(Fe)ZSM-5 catalysts for the existence of hotspot formation. Isomorphous substituted H-(Fe)ZSM-5 catalysts are tuned to identify key design parameters that control their microwave sensitivity. Increasing the amount of Al and Fe substitution in the zeolite lattice is found to result in higher microwave sensitivity due to improved dielectric properties and hence facilitated microwave heating of the catalyst bed. Comparison between the performance of these catalysts in a conventional thermally-heated (CH) fixed-bed and a microwave (MW) reactor allows identification of the effect of microwave irradiation on the catalyst activity. Methane conversion is drastically enhanced at MW conditions (from 3% at CH conditions to 40% at MW conditions), indicating accelerated methane activation over the metal site and hence suggesting a hot metal site despite much lower catalyst bulk temperature measured in the microwave reactor. The existence of metal hotspot formation is further supported by the reaction product distribution and spent catalyst analysis. At MW conditions, C 2 (ethane and ethylene) and coke (mainly carbon nanotubes and nanofibers) selectivity are much higher at the cost of aromatics make compared to CH conditions. This can be explained by reduced aromatization activity at the Brønsted acid sites of the zeolite due to its lower relative temperature of the zeolite. The much-enhanced metal aggregation observed in the spent catalyst from the MW reactor along with the distribution of coke species further confirms the existence of selective heating (i.e. hotspot formation) occurring at the active metal site in the catalyst at MW conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Predicting the rheology of limestone calcined clay cements (LC 3 ): Linking composition and hydration kinetics to yield stress through Machine Learning

The physicochemical characteristics of calcined clay influence yield stress of limestone calcined clay cements (LC 3 ), but the independent influences the clay's physical and chemical characteristics as well as the effect of other variables on LC 3 rheology are less well-understood. Further, a relationship between LC 3 hydration kinetics and yield stress – important for informing mixture design – has not yet been established. Here, rheological properties were determined in pastes with varying water-to-solid ratio (w/s), constituent mass ratios (PC:metakaolin:limestone), limestone particle size and gypsum content. From these data, an ML model developed allowed the independent examination of the different mechanisms by which metakaolin fraction influences yield stress of LC 3 , identifying four predictors – packing index, Al 2 O 3 /SO 3 , total particle density and metakaolin fraction relative to limestone (MK/LS) – most significant for predicting LC 3 yield stress. A methodology based on kernel smoothing also identified hydration kinetics parameters best correlated with yield stress.

36 MATERIALS SCIENCE↗

Data-driven surrogates for high dimensional models using Gaussian process regression on the Grassmann manifold

This paper introduces a surrogate modeling scheme based on Grassmannian manifold learning to be used for cost-efficient predictions of high-dimensional stochastic systems. The method exploits subspace-structured features of each solution by projecting it onto a Grassmann manifold. This point-wise linear dimensionality reduction harnesses the structural information to assess the similarity between solutions at different points in the input parameter space. The method utilizes a solution clustering approach in order to identify regions of the parameter space over which solutions are sufficiently similarly such that they can be interpolated on the Grassmannian. In this clustering, the reduced-order solutions are partitioned into disjoint clusters on the Grassmann manifold using the eigen-structure of properly defined Grassmannian kernels and, the Karcher mean of each cluster is estimated. Then, the points in each cluster are projected onto the tangent space with origin at the corresponding Karcher mean using the exponential mapping. For each cluster, a Gaussian process regression model is trained that maps the input parameters of the system to the reduced solution points of the corresponding cluster projected onto the tangent space. Using this Gaussian process model, the full-field solution can be efficiently predicted at any new point in the parameter space. In certain cases, the solution clusters will span disjoint regions of the parameter space. In such cases, for each of the solution clusters we utilize a second, density-based spatial clustering to group their corresponding input parameter points in the Euclidean space. The proposed method is applied to two numerical examples. Here, the first is a nonlinear stochastic ordinary differential equation with uncertain initial conditions where the surrogate is used to predict the time history solution. The second involves modeling of plastic deformation in a model amorphous solid using the Shear Transformation Zone theory of plasticity, where the proposed surrogate is used to predict the full strain field of a material specimen under large shear strains.

42 ENGINEERING↗

Bridging interfacial properties and cell performance: A multiscale model for proton-exchange-membrane fuel cells

Here, to elucidate the impact of local interfaces on mass-transport resistance and overall cell performance of low-loaded proton-exchange-membrane fuel cells (PEMFCs), we present a multiscale modeling framework incorporating a novel modified agglomerate model. The model considers three distinct Pt-electrolyte interfaces: Pt on the carbon surface covered by either ionomer or water film and Pt inside carbon nanopores. Detailed mass-transport voltage-loss breakdowns reveal that coupled agglomerate-interface-scale mass transport dominates the mass-transport loss. The ionomer poisons the exterior-Pt surface through suppressing O 2 adsorption and intrinsic ORR activity, leading to low current-density performance. Conversely, interior-Pt interface enhances the kinetic performance but limits high current-density performance due to its low interfacial permeability. The exterior-Pt/water interface demonstrates superior kinetic performance and mass transport, though its practical implementation requires ensuring proton transport. By coupling the multiscale CL properties with ink parameters, the model identifies an optimal I to C ratio of approximately 0.5, a moderate value where the ionomer content is sufficient to guarantee proton transport without fully covering the Pt surface and forming large agglomeration, thus allowing the utilization of the Pt-water interface and avoiding high mass-transport loss. Overall, the model helps unravel limiting phenomena across different operating regimes and provides routes for optimizing performance.

Cell diagnostic↗

Design and performance analysis of a Swirl Pintle injector for a 1 MW th pressurized oxy-coal combustor

Pressurized Oxy-coal combustion has the potential to achieve high thermal efficiency by recovering the latent heat of steam from flue gases and providing options for up to 95% carbon capture. Here, this study proposes a novel burner technology for pressurized oxy-coal combustion and investigates its operability. Three burner versions are designed and tested; no swirl, swirl numbers 0.9 and 1.2 to identify optimum design parameters and performance. A commercially available Dantec Dynamics software is used to conduct the shadowgraph experiments to observe burner performance. Coal water slurry is used as fuel with three different slurry ratios where the coal percentage by weights were 30, 40 and 50. Droplet size, distribution, spray angle and jet breakup lengths are measured and analyzed to quantify the performance. The study found that increasing the swirl number enhanced atomization. Implementing swirlers allowed the jet to spread through faster primary jet separation, thus improving atomization. In addition, higher slurry ratios improved the generation of smaller particles. However, it was also found that excessive increases in oxidizer momentums negatively impact the atomization process by jet shrinkage. The study concluded that operating a pintle burner at 1.2 swirl conditions with a 40/60 slurry ratio produces optimum performance.

42 ENGINEERING↗

Developing custom fire behavior fuel models from ecologically complex fuel structures for upper Atlantic Coastal Plain forests.

Currently geospatial fire behavior analyses are performed with an array of fire behavior modeling systems such as FARSITE, FlamMap, and the Large Fire Simulation System. These systems currently require standard or customized surface fire behavior fuel models as inputs that are often assigned through remote sensing information. The ability to handle hundreds or thousands of measured surface fuelbeds representing the fine scale variation in fire behavior on the landscape is constrained in terms of creating compatible custom fire behavior fuel models. In this study, we demonstrate an objective method for taking ecologically complex fuelbeds from inventory observations and converting those into a set of custom fuel models that can be mapped to the original landscape. We use an original set of 629 fuel inventory plots measured on an 80,000 ha contiguous landscape in the upper Atlantic Coastal Plain of the southeastern United States. From models linking stand conditions to component fuel loads, we impute fuelbeds for over 6000 stands. These imputed fuelbeds were then converted to fire behavior parameters under extreme fuel moisture and wind conditions (97th percentile) using the fuel characteristic classification system (FCCS) to estimate surface fire rate of spread, surface fire flame length, shrub layer reaction intensity (heat load), non-woody layer reaction intensity, woody layer reaction intensity, and litter-lichen-moss layer reaction intensity. We performed hierarchical cluster analysis of the stands based on the values of the fire behavior parameters. The resulting 7 clusters were the basis for the development of 7 custom fire behavior fuel models from the cluster centroids that were calibrated against the FCCS point data for wind and fuel moisture. The latter process resulted in calibration against flame length as it was difficult to obtain a simultaneous calibration against both rate of spread and flame length. The clusters based on FCCS fire behavior parameters represent reasonably identifiable stand conditions, being: (1) pine dominated stands with more litter and down woody debriscomponents than other stands, (2) hardwood and pine stands with no shrubs, (3) hardwood dominated stands with low shrub and high non-woody biomass and high down woody debris, (4) stands with high grass and forb (i.e., non-woody) biomass as well as substantial shrub biomass, (5) stands with both high shrub and litter biomass, (6) pine-mixed hardwood stands with moderate litter biomass and low shrub biomass, and (7) baldcypress-tupelo stands. Models representing these stand clusters generated flame lengths from 0.6 to 2.3 musing a 30 km h{sub 1} wind speed and fireline intensities of 100-1500 kW m{sub 1} that are typical within the range of experience on this landscape. The fuel models ranked 1 < 2 < 7 < 5 < 4 < 3 < 6 in terms of both flame length and fireline intensity. The method allows for ecologically complex data to be utilized in order to create a landscape representative of measured fuel conditions and to create models that interface with geospatial fire models.

09 BIOMASS FUELS↗

Heat and mass transfer in hygroscopic hydrogels

Sorption and desorption with hygroscopic hydrogels hold significant promise for thermal management, passive cooling, thermal energy storage, and atmospheric water harvesting. However, a comprehensive understanding of the energy and mass transport mechanisms in hygroscopic hydrogels remains missing, impeding accurate modeling and optimization. In this paper, we develop a model for the simultaneous vapor, water, and heat transfer in hygroscopic hydrogels during sorption and desorption processes. We show that by considering vapor diffusion in the hydrogel micropores, water diffusion in the polymer mesh, and heat transfer in the porous hydrogel, we can accurately capture experimentally observed thermally-driven desorption rates in these hydrogels. Furthermore, we consider three typical operating configurations of hydrogels and elucidate the differences in the transport mechanisms depending on the configuration. Finally, for each of these configurations, we identify key design parameters, including hydrogel thickness, hydrogel shear modulus, heat transfer coefficient, and thermal conductivity, and we parametrically show that by varying these parameters, a hygroscopic hydrogel can desorb up to 128.5%, 14.9%, 69.7%, and 9.6% more water, respectively, relative to the initial water content. This work provides a generic framework to model sorption and desorption processes in hygroscopic hydrogels which can guide the design and optimization in applications of thermal management, passive cooling, thermal energy storage, and atmospheric water harvesting with hydrogels.

36 MATERIALS SCIENCE↗

An iterative particle filter approach for coupled hydro-geophysical inversion of a controlled infiltration experiment

The modeling of unsaturated groundwater flow is affected by a high degree of uncertainty related to both measurement and model errors. Geophysical methods such as Electrical Resistivity Tomography (ERT) can provide useful indirect information on the hydrological processes occurring in the vadose zone. In this paper, we propose and test an iterated particle filter method to solve the coupled hydrogeophysical inverse problem. We focus on an infiltration test monitored by time-lapse ERT and modeled using Richards equation. The goal is to identify hydrological model parameters from ERT electrical potential measurements. Traditional uncoupled inversion relies on the solution of two sequential inverse problems, the first one applied to the ERT measurements, the second one to Richards equation. This approach does not ensure an accurate quantitative description of the physical state, typically violating mass balance. To avoid one of these two inversions and incorporate in the process more physical simulation constraints, we cast the problem within the framework of a SIR (Sequential Importance Resampling) data assimilation approach that uses a Richards equation solver to model the hydrological dynamics and a forward ERT simulator combined with Archie's law to serve as measurement model. ERT observations are then used to update the state of the system as well as to estimate the model parameters and their posterior distribution. The limitations of the traditional sequential Bayesian approach are investigated and an innovative iterative approach is proposed to estimate the model parameters with high accuracy. The numerical properties of the developed algorithm are verified on both homogeneous and heterogeneous synthetic test cases based on a real-world field experiment.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Adaptive sparse polynomial dimensional decomposition for derivative-based sensitivity

In this report applications, such as the plasma-coupled combustion system we consider, derivative-based sensitivity indices (DSI) are known to have several advantages over Sobol's total sensitivity indices, especially for small sample sizes. Several properties of derivative-based sensitivity measures are leveraged to develop a new and efficient numerical approach to estimate them. It is based on computing the DSI measures by effectively cost-free Monte Carlo sampling of an adaptively constructed orthogonal polynomial surrogate with uncertain input parameters that can have arbitrary probability distributions. The adaptivity reduces the number of necessary model evaluations, which is demonstrated both in a constructed example (the Moon function) and in two plasma-combustion systems with up to 55 uncertain parameters. Unimportant parameters are successfully identified and neglected with a low number of model evaluations, which makes it an attractive non-intrusive approach when adjoint solutions are unavailable to provide sensitivity information.

42 ENGINEERING↗

A composite neural network that learns from multi-fidelity data: Application to function approximation and inverse PDE problems

Presently the training of neural networks relies on data of comparable accuracy but in real applications only a very small set of high-fidelity data is available while inexpensive lower fidelity data may be plentiful. We propose a new composite neural network (NN) that can be trained based on multi-fidelity data. It is comprised of three NNs, with the first NN trained using the low-fidelity data and coupled to two high-fidelity NNs, one with activation functions and another one without, in order to discover and exploit nonlinear and linear correlations, respectively, between the low-fidelity and the high-fidelity data. We first demonstrate the accuracy of the new multi-fidelity NN for approximating some standard benchmark functions but also a 20-dimensional function that is not easy to approximate with other methods, e.g. Gaussian process regression. Subsequently, we extend the recently developed physics-informed neural networks (PINNs) to be trained with multi-fidelity data sets (MPINNs). MPINNs contain four fully-connected neural networks, where the first one approximates the low-fidelity data, while the second and third construct the correlation between the low- and high-fidelity data and produce the multi-fidelity approximation, which is then used in the last NN that encodes the partial differential equations (PDEs). Specifically, by decomposing the correlation into a linear and nonlinear part, the present model is capable of learning both the linear and complex nonlinear correlations between the low- and high-fidelity data adaptively. By training the MPINNs, we can: (1) obtain the correlation between the low- and high-fidelity data, (2) infer the quantities of interest based on a few scattered data, and (3) identify the unknown parameters in the PDEs. In particular, we employ the MPINNs to learn the hydraulic conductivity field for unsaturated flows as well as the reactive models for reactive transport. The results demonstrate that MPINNs can achieve relatively high accuracy based on a very small set of high-fidelity data. Despite the relatively low dimension and limited number of fidelities (two-fidelity levels) for the benchmark problems in the present study, the proposed model can be readily extended to very high-dimensional regression and classification problems involving multi-fidelity data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Alteration of natural 37 Ar activity concentration in the subsurface by gas transport and water infiltration

High 37 Ar activity concentration in soil gas is proposed as a key evidence for the detection of underground nuclear explosion by the Comprehensive Nuclear Test-Ban Treaty. However, such a detection is challenged by the natural background of 37 Ar in the subsurface, mainly due to Ca activation by cosmic rays. A better understanding and improved capability to predict 37 Ar activity concentration in the subsurface and its spatial and temporal variability is thus required. A numerical model integrating 37 Ar production and transport in the subsurface is developed, including variable soil water content and water infiltration at the surface. A parameterized equation for 37 Ar production in the first 15 m below the surface is studied, taking into account the major production reactions and the moderation effect of soil water content. Using sensitivity analysis and uncertainty quantification, a realistic and comprehensive probability distribution of natural 37 Ar activity concentrations in soil gas is proposed, including the effects of water infiltration. Site location and soil composition are identified as the parameters allowing for a most effective reduction of the possible range of 37 Ar activity concentrations. The influence of soil water content on 37Ar production is shown to be negligible to first order, while 37 Ar activity concentration in soil gas and its temporal variability appear to be strongly influenced by transient water infiltration events. Overall, these results will be used as a basis for practical CTBTO concepts of operation during an OSI.

07 ISOTOPE AND RADIATION SOURCES↗

Mechanical behavior of bimetallic stainless steel and gray cast iron repairs via directed energy deposition additive manufacturing

The utility of gray cast iron in engine components remains tied to the mechanical performance and cost. Repair and remanufacturing of castings offer economical and sustainable benefits; however, high thermal input from traditional fusion-based welding is unable to restore the original mechanical quality owing to brittle microstructures and porosity formed in situ. Directed energy deposition (DED) is an additive manufacturing method that has received considerable interest for repairs owing to the highly controllable nature of the process. Despite this, few works have connected the effect of DED parameters on actual interfacial strength. Consequently, distinct DED parameter combinations were identified to maximize the strength and fatigue life of the repaired cast iron. Further, high speed melt pool imaging and residual stress measurements are provided to aid in the understanding of the metallurgical quality and strength seen in these structures. In general, higher scanning speeds and lower thermal gradients promoted comparable tensile strength to that of the original gray cast iron. The results presented here provide a foundation to tune in the DED process to generate the required mechanical quality as a starting point for future process advancements.

36 MATERIALS SCIENCE↗

A finite viscoelastic constitutive model for low to high strain rate response of elastomers with application of strain rate-induced glass transition

Amorphous elastomers exhibit significant rate-stiffening and unique viscous flow characteristics across a wide range of strain rates, often undergoing glass transition above a strain rate threshold. We have developed a thermodynamically-consistent and micromechanically-inspired constitutive model for soft elastomeric materials to capture the rate-dependent stress-strain behavior and hysteresis when subjected to low to high strain rates. Here, our proposed constitutive model encapsulates the viscous flow of materials through molecular motion at low strain rates and local rearrangement and alignment of the molecules trying to overcome the intermolecular resistance at high strain rates, essentially covering the glass transition. We applied our constitutive model to uniaxial compression experiments performed at low and high strain rates for polyborosiloxane (PBS) to identify the material parameters, and subsequently, performed numerical simulations of single and multi-cycle compression, stress relaxation, and small amplitude oscillatory tension-compression. Our analyses indicate that the model predicts higher total energy dissipation with increasing strain rate; however, dissipation associated with molecular relaxation decreases (forming a cusp) because, beyond a crossover strain rate, intermolecular rearrangement and alignment become dominant, which is consistent with the onset of the glass transition. For cyclic loading-unloading, we observed that dissipation over a cycle remains constant at low strain rates but decreases non-monotonically at high strain rates before becoming constant, with the peak stress over the cycle becoming higher, which can be interpreted as more loading being carried elastically by the polymer network as the intermolecular rearrangement process occurs. Additionally, our model was able to predict the qualitative nature of the storage modulus and loss modulus in the limit of small strain over a wide range of frequency sweeps.

36 MATERIALS SCIENCE↗

Fracture mechanics approach to TRISO fuel particle failure analysis

Weibull stress-based methods for failure probability assessment have been developed and analyzed to assess the integrity of tristructural isotropic (TRISO) fuel particles during fuel life cycles and accident operating conditions. While simple, these methods entail a number of drawbacks when stress concentrates near crack tips, including finite element mesh size dependency when the Weibull stress is averaged over the finite element domain. Fracture mechanics approaches involving the use of interaction integrals eliminate this lack of mesh convergence and produce consistent fracture predictions. In this work, we use an interaction integral approach to computing stress intensity factors for a crack in the inner pyrolytic carbon layer perpendicular to the silicon carbide layer, which is simplified representation of a failure mode in TRISO particles. Further, the interface between these two TRISO layers has been shown to become porous, which we simulate by considering a transition of mechanical properties over such porous length, i.e. the layers are modeled as a functionally graded material. Aspects such as porosity and thermal and irradiation eigenstrains are considered in computing the stress intensity factor from a fracture mechanics approach and compared with the known Weibull stress failure approach. The methodology introduced in this paper enables a more general fracture probability assessment in TRISO particles and eliminates the need to identify best suited parameters when using local or averaged stress-based failure criterion. Finally, the numerical sensitivity studies show how parameters such as the porous transition zone length, the material stiffness, and creep affect the probability of TRISO fuel particle failure.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Analysis of electrolyte imbibition through lithium-ion battery electrodes

A quantitative measurement of wettability between the porous electrode and the electrolyte in lithium-ion batteries can greatly improve our understanding of wetting behavior. Although the wetting balance method is widely used to measure the electrolyte transport rate in the porous electrodes, it suffers from several drawbacks and has limited accuracy. We here presented a combined experimental and theoretical investigation of the dynamics of electrolyte imbibition through electrodes. We suggested a novel method to accurately measure the electrolyte imbibition rate. Excellent agreement between the experimental data and the developed analytical model is obtained. The coefficient of penetrance (COP) and the solid permeability coefficient (SPC) are identified as important parameters, i.e., the electrolyte with higher COP value wets faster into an electrode, whereas for an electrolyte, the electrode with higher SPC value is more amenable to be impregnated. The effect of electrolyte salt concentration and electrolyte solvent has been studied in detail. The result implies that increasing salt concentration adversely influences electrolyte wetting rate, whereas switching from EC-DEC system to EC-EMC system improves electrolyte wetting rate. Moreover, for the electrolytes tested in this study, the imbibition into the uncalendered graphite anode is much faster than that into the uncalendered NMC532 cathode.

25 ENERGY STORAGE↗

Water limiting current measurements in anion exchange membrane fuel cells (AEMFCs); part 1: Water limiting current method development

We report water management is critical in optimizing anion exchange membrane fuel cell (AEMFCs) performance and durability. The role of membranes on determining cell water balance in AEMFCs has received limited study. A more thorough understanding of membrane water transport parameters is necessary for thorough cell optimization. In this study, a novel method for probing limiting current based on water flux across an AEMFC is presented for the first time. The water limiting current has been investigated as a function of relative humidity of anode and cathode streams, cathode backpressure and cell operating temperature. From these studies, it is clearly shown that water diffusion through the membrane from the anode is the most critical mechanism for achieving high current density performance and durability. The developed method can act as a novel screening technique that focuses on the impacts of water transport parameters and helps identify promising membrane materials based on their water transport properties. A critical but underappreciated facet of AEMFC performance and durability.

30 DIRECT ENERGY CONVERSION↗

Computational investigation of arranged scintillating particle composites for fast neutron detection

Composite materials have the potential to play an important role in enabling high-performance and cost-effective neutron detectors for fundamental science, global security, and dosimetry applications. This work presents neutron and gamma transport simulation results on composites consisting of 6 Li-containing GS20 scintillator glass particles arranged in an organic matrix. These composites achieve high sensitivity for neutrons while suppressing signals from gamma rays. The simulations successfully identified optimal composite parameters such as particle size, geometry, and inter-particle pitch. This information is key to the subsequent fabrication and scale-up of the desired GS20-in-acrylic composite. Finally, the simulation results presented in this paper demonstrate that arranged composites can provide high neutron detection efficiency with gamma misidentification values of <1 in 10 6 gamma rays, assuming irradiation with a 252 Cf spontaneous fission source.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗