Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mathematical model”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Multiphase Reacting Flow Simulations and Optimization of Commercial-Scale Aerobic Bioreactors

Aerobic bioreactors are widely used in the synthesis of value-added products in pharmaceutical and biotechnology industries, wherein microbial bioreactions aid in the conversion of sugars to macromolecules. These reactors are also being actively investigated as a cost-effective pathway for the production of low-value commodities such as biofuels and animal feed. Reactor scale-up is one of the important challenges when designing these systems. Gas and liquid phase transport, mass-transfer, and mixing physics at large length scales can significantly affect microbial conversion rates, particularly when the microbial reaction requires a narrow set of conditions. These phenomena are difficult to study in small-scale bench-top reactors that are typically well-mixed. Predictive computational fluid dynamics (CFD) based simulations can therefore aid in the design and optimization of these reactors. This work presents multiphase Euler-Euler CFD simulations of commercial-scale (~ 500 m3) aerobic bioreactors. Our mathematical model treats the gas and liquid as interpenetrating phases. This approach reduces the computational complexity of tracking individual gas bubbles that are several orders of magnitude smaller than reactor dimensions. We solve the Reynolds averaged Navier-Stokes (RANS) multiphase equations that account for phase and chemical species transport, interphase mass and momentum transfer and uses a phenomenological model for oxygen uptake by microbes. We use a customized solver derived from open-source CFD toolbox, OpenFOAM, to perform these simulations, which has been validated against small-scale reactors in our previous work. This work examines the performance of three different reactor designs, viz. bubble column reactor, airlift reactor with an internal draft tube, and a stirred-tank reactor with Rushton impellers. Reactor oxygen mass-transfer coefficient, gas hold-up, and oxygen distribution are critically analyzed among reactors, and sensitivity studies pertaining to gas flow rates and reactor geometry will be presented.

09 BIOMASS FUELS↗
DOE OSTI2020

System and method for measuring fluid properties using magnetic field techniques via magnetic tracer

A method and apparatus for determining a property of a fluid in a vessel. The method uses magnetic tracer particles and an externally applied magnetic field which orients the particles. When the fluid moves, it changes the orientation of the tracer particles, thus changing the magnetic fields. These changes are detected by external magnetic field sensors. By using mathematical models, the property of the fluid in the vessel is determined from the detected magnetic field. In this manner, fluid vorticity, velocity, strain and stress may be estimated.

42 ENGINEERING↗
DOE OSTI2021

Kinetic Measurements in Heterogeneous Catalysis

This contribution is about the experimental determination of the rate of a heterogeneous catalytic reaction and the analysis of kinetic data. In this case, the reaction rate can be defined as the frequency at which the closed sequence of elementary steps transforming reactants into products is occurring. However, the rate of chemical reaction is not directly observed; rather, one records the rate of substance change. The rate of chemical reaction is calculated based on the rate of substance change and assumed stoichiometry of the reaction. In order to make this an intensive quantity (i.e., independent of the volume in which the reaction is performed, this frequency is divided by the latter or by a quantity proportional to the latter). As the focus of what follows is on heterogeneously catalyzed reactions this can be the catalyst mass, the catalyst surface area or the total number of active sites present in the reaction volume confined by the walls of a chemical reactor. For example, the rate may be reported as one of the following: Net rate of consumption/production of component i mol m-3 s-1 Net specific rate of consumption/production of component i mol kg-1 s-1 Rate of substance change (per unit volume of catalyst) mol m-3 s-1 Specific reaction rate (per unit mass of catalyst) mol kg-1 s-1 Net rate of consumption/production of component i mol m-2 s-1 The turnover frequency (TOF) is a characteristic originally introduced by Boudart [1] is expressed as TOF=R/G_tot , where R is the steady-state rate of reactant consumption or product generation and G_totis the total areal density of active sites [mol m-2]. The latter is typically taken from experimental chemisorption data obtained at low temperature. The advancement of a chemical reaction leads to changes in the amounts of reactants and products which, for a stoichiometric single reaction, are connected by the stoichiometric coefficients. In many cases, for the sake of simplicity, we will assume that the investigated reaction proceeds according to a single reaction pathway, so that the net production rate of a reaction component is directly proportional to the reaction rate and that the latter equals the net production rate of any of the involved components when divided by the appropriate stoichiometric coefficient. For complex stoichiometric reactions the net production rates of the involved components are linear combinations of the reaction rates. Different goals for kinetic measurements can be formulated: Catalyst testing, i.e., obtaining the kinetic dependences for the development and discrimination of efficient catalytic materials. Precise kinetic characterization of active materials via high-throughput screening is a part of such activity. Detailed kinetics, i.e., revealing the detailed mechanism of complex catalytic reaction via systematic kinetic studies, both steady-state and non-steady-state Industrial kinetics, i.e., obtaining data and relationships for describing and predicting the behaviour of catalytic reactors and processes at industrial scale Finally, mathematical modeling and analysis where kinetic measurements provide data for understanding complex kinetic phenomena, e.g., oscillations, non-linear self-organization, etc.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Gamma Spectrometry Examination of the AGR-3/4 Irradiation

The results from gamma spectrometry examination of the different components from the combined third and fourth U.S. Advanced Gas Reactor (AGR) TRISO-coated particle fuel irradiation tests (AGR-3/4) have been analyzed. This experiment was designed to provide information about in-pile fission product migration. In each of the 12 capsules, a single stack of four compacts with designed-to-fail (DTF) particles surrounded by inner and outer graphite and/or graphitic matrix rings and a graphite sink ring were irradiated in the Advanced Test Reactor (ATR) at Idaho National Laboratory (INL). Gamma spectrometry has been used to evaluate the gamma-emitting fission product inventory of compacts from the irradiation and evaluate the burnup of these compacts based on the activity of the radioactive cesium isotopes (i.e., Cs-134 and Cs-137) in the compacts. Burnup from gamma spectrometry compares well with predicted burnup from simulations. The inner and outer rings were also examined by gamma spectrometry to evaluate the total fission product inventory of the rings, and gamma emission computed tomography (GECT) was used to investigate the spatial distribution of gamma-emitting fission products within the rings. This report focuses on results obtained via the Precision Gamma Scanner (PGS) for the compacts and inner and outer rings. These non-destructive gamma measurements are currently being compared to recent destructive measurements of spatial fission product distributions in the rings and to predictions made via mathematical modeling. A separate report discusses the entire AGR-3/4 fission product mass balance for all irradiation capsule components including the inner and outer rings, graphite sink rings, capsule spacers, capsule foils, through tubes, and felts. This report is similar to the conference proceedings: Harp, J.M., Demkowicz, P.A., and Stempien, J.D., “Initial gamma spectrometry examination of the AGR-3/4 irradiation,” International Topical Meeting High Temperature Reactor Technology (HTR 2016), Las Vegas, NV, USA, November 2016, Paper HTR2016-18593, but it has been updated to reflect all the data that has been collected on AGR-3/4 using PGS. The Appendices contain extensive data from the collected gamma spectra.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Conservation of Helium while Maintaining High System Purity

Recent helium shortages and helium price increases have lead to an increased emphasis being placed on conserving helium. The need to conserve helium must be balanced with need to maintain the high levels of purity necessary to prevent operational problems caused by contamination. Helium losses and contamination control are especially important for test stands that have cryogenic distribution systems operating continuously with frequent changeover of cryogenic temperature components that are being tested. This paper describes a mathematical model to estimate the quantity of helium lost and the purity of the helium after the pump and backfill procedure is complete. The process to determine the optimal time during pump down to cut off pumping and start backfilling is described. There is a tradeoff between trying to achieve the lowest possible pressure during pumping and the quantity of air leaking into the volume while pumping is occurring. An additional benefit of careful selection of pump and backfill parameters in conjunction with real-time pressure monitoring can reduce the labor and time required to complete a successful pump and backfill procedure. This paper is intended to be a tool for engineers to review their pump and backfill procedures and measured data to optimize helium losses, system purity, and labor required.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Droning on to Delivery: Examining the Energy Impacts of Using Drones for Moving Goods

The demand for fast, localized delivery has grown significantly in recent years. Whether for cheesy snacks, prepared food, medical solutions, or business deliveries, fast and efficient delivery is increasingly a demand and differentiator. Drone delivery in the freight sector offers to revolutionize last-mile logistics and improve services. This research analyzes the impacts of drone energy for delivery operations, aiming to compare various types of drones (large and small, rotary and VTOL) and various types of business methods. The project executes novel open-air and laboratory-based testing to look at the impacts of weights, operations, temperatures and weather conditions. The study combines this real-world experimental data with fleet optimization mathematical models to assess energy consumption of different scenarios as well as examine the minimum fleet size and additional battery requirements. These fleet optimization models compare energy between types of deployments with existing delivery methods. The model was also extended to look at mixed fleets of aerial drones and ground vehicles to accommodate different restrictions on drones or when weather prohibits their use. The analysis provided insights into factors such as drone design, payload weight, flight distance, weather conditions and operational parameters and impacts of each. It showed how to combine different types of vehicles to reduce energy and improve services. And it showed that drone speed, routing restrictions, and unfavorable weather significantly influence energy consumption. Our open experiment data and optimization models can assist stakeholders and industry in understanding drone package delivery and offers keys to improving deployment.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Fundamental Path Optimization Strategies for Extrusion-based Additive Manufacturing

Extrusion-based additive manufacturing processes begin with a software program, called a slicer, that generates layer geometry and fits toolpaths to each layer to define where material is to be extruded or deposited. Before the toolpaths are output as g-code for the additive manufacturing system to execute, the toolpaths should be optimized. Many complex optimization approaches using graph theory, Chinese postman problem, and other complex mathematical models exist, but these approaches are rarely used in daily printing operations and are not available through common slicing programs such as Cura and PrusaSlicer. Instead, path planning and optimization typically revolves around simpler, fully automated approaches such as inside out and next closest. This paper will explore the fundamental optimization strategies for toolpath planning and document a new implementation, available via open-source slicing software, that allows for greater control of the path planning process.

Roschli, Alex [ORNL] (ORCID:0000000213084632)↗
DOE OSTI2025

Neural Network Analysis of Nuclear Magnetic Resonance and Infrared Spectra

Nuclear magnetic resonance (NMR) spectroscopy and infrared (IR) spectroscopy are powerful chemical characterization techniques with broad general usage. However, the manual evaluation of the resulting spectra is time-consuming and requires significant expertise, preventing insights from being used in real-time applications. With recent advances in computation and artificial intelligence (AI), new tools are available for automating spectral interpretation. In this work, machine learning (ML) algorithms using 1-dimensional convolutional neural networks (CNNs) were applied to identify common functional groups from spectral information. Raw spectra were collected virtually from the Human Metabolome Database (HMDB) and National Institute of Standards and Technology (NIST) Chemistry WebBook and processed into a suitable standard. Algorithm design was tailored to best fit the nature of the problem, with built-in flexibility to accommodate relevant parameters beyond the raw spectral input, specifically solvent identity and magnetic frequency for NMR. The predictive capability of the algorithm in identifying functional groups is displayed in several examples. This methodology has been compiled into a code repository and could easily be modified to adapt alternative data sources, including other spectrum types. To mitigate overfitting, a common problem in mathematical modeling where overfamiliarity with training data produces trends that are not representative of the general data, a novel metric was developed, referred to as Accufit. Accufit includes a parameter that penalizes substantial differences in the training accuracy and the accuracy of an independent validation set. Examples are presented showing the effectiveness of Accufit in maintaining the model’s predictive capability while controlling the overfitting when used as a custom metric for hyperparameter tuning.

Sturgill, James↗
DOE OSTI2024

Statistical Estimation of EV Driver Charging Behavior and Influential Factors

INL received data collected via telematics from battery electric vehicles (BEVs), and these vehicles were owned by retail customers who had entered into a telematics user agreement. The goal of analyzing these data was to develop mathematical models to characterize how different sets of BEV drivers use charging infrastructure at home and away from home (i.e., public charging) and quantify how various factors influence BEV drivers’ decision to charge and use available infrastructure. The data used in this analysis are unique because they provide real world BEV driving and charging behavior at the individual driving and parking event level. In this study we seek to leverage this data to quantify BEV charging and driving metrics to help inform models that predict quantities like the specific times when loads are imposed on the electrical grid due to BEV charging. Most models that have been developed to predict electrical grid load due to BEV charging, use simulations of BEV driving events and rely on assumptions such as every vehicle charges every night. Using a statistical modelling framework, we seek to investigate BEV charging behavior and quantitatively assess these common assumptions of BEV charging behavior.

33 - ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

Collaboratory on Mathematics for Mesoscopic Modeling of Materials (CM4)(Final Technical Report)

Our work on the CM4 grant focused on new mathematical and computational approaches for stochastic methods for mesoscopic systems. Our work included continuum SPDEs, molecular coarse-grained models, and hybrid particle-continuum approaches. This work was motivated by applications in soft materials, fluidic devices, and interfacial phenomena. Throughout the project we collaborated with several staff researchers, including from Pacific Northwest National Laboratories (PNNL) and Sandia National Laboratory (SNL) as well as other members of the CM4 project. Below we summarize briefly a few specific accomplishments from this work over the funding period. A full list of published papers and software releases supported by the grant are listed at the end.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Security-Constrained Unit Commitment for Electricity Market: Modeling, Solution Methods, and Future Challenges

This paper summarizes the technical activities of the IEEE Task Force on Solving Large Scale Optimization Problems in Electricity Market and Power System Applications. This Task Force was established by the IEEE Technology and Innovation Subcommittee to first review the state-of-the-art of the security-constrained unit commitment (SCUC) business model, its mathematical formulation, and solution techniques in solving electricity market clearing problems. The Task Force then investigated the emerging challenges of future market clearing problems and presented efforts in building benchmark mathematical and business models.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Mathematical analysis of laboratory microbial experiments demonstrating deterministic chaotic dynamics

Presented is a system of four ordinary differential equations and a mathematical analysis of microbiological experiments in a four-component chemostat—nutrient n, rods r, cocci c, and predators p. The analysis is consistent with the conclusion that previous experiments produced features of deterministic chaotic and classical dynamics depending on dilution rate. The surrogate model incorporates as much experimental detail as possible, but necessarily contains unmeasured parameters. The objective is to understand better the differences between model simulations and experimental results in complex microbial populations. The key methodology for simulation of chaotic dynamics, consistent with the measured dilution rate and microbial volume averages, was to cause the preference of p for r vs. c to vary with the r and c concentrations, to make r more competitive for nutrient than c, and to recycle some dying p biomass, leading to a modified version of the Monod kinetics model. Our mathematical model demonstrated that the occurrence of chaotic dynamics requires a predator, p, preference for r versus c to increase significantly with increases in r and c populations. Also included is a discussion of several generalizations of the existing model and a possible involvement of the minimum energy dissipation principle. This principle appears fundamental to thermodynamic systems including living systems. Several new experiments are suggested.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Sparsifying priors for Bayesian uncertainty quantification in model discovery

We propose a probabilistic model discovery method for identifying ordinary differential equations governing the dynamics of observed multivariate data. Our method is based on the sparse identification of nonlinear dynamics (SINDy) framework, where models are expressed as sparse linear combinations of pre-specified candidate functions. Promoting parsimony through sparsity leads to interpretable models that generalize to unknown data. Instead of targeting point estimates of the SINDy coefficients, we estimate these coefficients via sparse Bayesian inference. The resulting method, uncertainty quantification SINDy (UQ-SINDy), quantifies not only the uncertainty in the values of the SINDy coefficients due to observation errors and limited data, but also the probability of inclusion of each candidate function in the linear combination. UQ-SINDy promotes robustness against observation noise and limited data, interpretability (in terms of model selection and inclusion probabilities) and generalization capacity for out-of-sample forecast. Sparse inference for UQ-SINDy employs Markov chain Monte Carlo, and we explore two sparsifying priors: the spike and slab prior, and the regularized horseshoe prior. UQ-SINDy is shown to discover accurate models in the presence of noise and with orders-of-magnitude less data than current model discovery methods, thus providing a transformative method for real-world applications which have limited data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Verifying infectious disease scenario planning for geographically diverse populations

In the face of the COVID-19 pandemic, the literature saw a spike in publications for epidemic models, and a renewed interest in capturing contact networks and geographic movement of populations. There remains a general lack of consensus in the modeling community around best practices for spatiotemporal epi-modeling, specifically as it pertains to the infection rate formulation and the underlying contact or mixing model. We mathematically verify several common modeling assumptions in the literature, to prove when certain choices can provide consistent results across different geographic resolutions, population densities and patterns, and mixing assumptions. The most common infection rate formulation, a computationally low cost per capita infection rate assumption, fails the consistency tests for heterogeneous populations and gravity-weighting assumptions. Future modeling efforts in spatiotemporal disease modeling should be wary of this limitation, particularly when working with more heterogeneous or sparse populations. Our results provide guidance for testing that a model preserves desirable properties even when model inputs mask potential problems due to symmetry or homogeneity. We also provide a recipe for performing this type of verification, strengthening decision support tools.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Particle size influences decay rates of environmental DNA in aquatic systems

Abstract Environmental DNA (eDNA) analysis is a powerful tool for remote detection of target organisms. However, obtaining quantitative and longitudinal information from eDNA data is challenging, requiring a deep understanding of eDNA ecology. Notably, if the various size components of eDNA decay at different rates, and we can separate them within a sample, their changing proportions could be used to obtain longitudinal dynamics information on targets. To test this possibility, we conducted an aquatic mesocosm experiment in which we separated fish‐derived eDNA components using sequential filtration to evaluate the decay rate and changing proportion of various eDNA particle sizes over time. We then fit four alternative mathematical decay models to the data, building towards a predictive framework to interpret eDNA data from various particle sizes. We found that medium‐sized particles (1–10 μm) decayed more slowly than other size classes (i.e., <1 and > 10 μm), and thus made up an increasing proportion of eDNA particles over time. We also observed distinct eDNA particle size distribution (PSD) between our Common carp and Rainbow trout samples, suggesting that target‐specific assays are required to determine starting eDNA PSDs. Additionally, we found evidence that different sizes of eDNA particles do not decay independently, with particle size conversion replenishing smaller particles over time. Nonetheless, a parsimonious mathematical model where particle sizes decay independently best explained the data. Given these results, we suggest a framework to discern target distance and abundance with eDNA data by applying sequential filtration, which theoretically has both metabarcoding and single‐target applications.

Brandão‐Dias, Pedro F. P.↗
DOE OSTI2023

Relaxing the constant molar overflow assumption in distillation optimization

Abstract The constant molar overflow (CMO) framework, while useful for shortcut distillation models, assumes that all components have the same latent heats of vaporization. A simple transformation, from molar flows to latent‐heat flows, allows shortcut models to retain the mathematical simplicity of the CMO framework while accounting for different latent heats, resulting in the constant heat transport (CHT) framework for adiabatic distillation columns. Although several past works have already proposed this transformation in the literature, it has not been well utilized in recent times. In this article, we show the utility of this transformation in upgrading various applications such as identifying energy‐efficient multicomponent distillation configurations based on heat duty rather than surrogate vapor flow. The method transforms the diagram to a diagram. Furthermore, we derive new and insightful analytical results in distillation, such as cumulative latent‐heat stage fractions having monotonic profiles within a distillation column under the CHT framework.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Comparison of steady-state analytical wake models implemented in wind farm analysis software

A common set of mathematical wind turbine wake models are implemented in a few, well-adopted computational tools for wind farm wake modelling. Although the referenced mathematical formulations are common, implementation details may lead to differences in results. This study presents a systematic comparison of the implementation of mathematical wake models in open source, Python-based wind turbine wake modelling software, and a set of the models are directly compared. Despite aligning only the mathematical model parameters and retaining the default computational model parameters, good agreement is found across most of the model implementations, and additional agreement is expected upon further parameters alignment.

17 WIND ENERGY↗