Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model based development”

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 181 records · Page 10

An analysis of parameters affecting ampacity in aircraft bipolar MVDC power cables via coupled electrical, thermal, and computational fluid dynamic modelling

Abstract The next generation of aircraft, including more electric aircraft and all‐electric aircraft (AEA), requires electric power systems with high power density and low system mass specifications. Increasing the voltage of the system to the range of a few kV, medium voltage (MV), is a reasonable approach to achieving high‐power‐density and low‐system‐mass EPSs for aircraft applications. Higher voltages, however, pose many challenges for aviation MV power cables such as arcs and arc tracking, partial discharges (PDs), and thermal management. In this regard, thermal management is more challenging since heat transfer by convection is greatly reduced at wide‐body aircraft's cruising altitudes due to the reduced air pressure. In this paper, a finite element method (FEM) model is developed in COMSOL Multiphysics for an aircraft bipolar MVDC (±5 kV) power cable. Using the model, the maximum permissible cable current at a low pressure of 18.8 kPa (at an altitude of 12.2 km from sea level, the usual cruising altitude for wide‐body aircraft) is calculated. Also, an analytical model is developed based on analytical and proven empirical correlations governing conductive, radiative, and convective heat transfers at the steady state to estimate the ampacity of the bipolar cable system at reduced pressure. It was shown that the proposed analytical model can be used for atmospheric pressure and systems with a larger number of poles, expanding its range of applications. The results of the FEM and analytical models correlate at wide ranges of parameters such as ambient temperature, duct size, distance between the positive and negative pole cables, and the overall diameter of the cables. The influence of horizontal and vertical arrangement of poles is included in the analytical model. The results of this study can be used to design bipolar MVDC power cable systems for the envisaged wide‐body AEA.

Azizi, Arian↗

A pencil beam dose calculation model for CyberKnife system

Purpose: CyberKnife system is initially equipped with fixed circular cones for stereotactic radiosurgery. Two dose calculation algorithms, Ray-Tracing and Monte Carlo, are available in the supplied treatment planning system. A multileaf collimator system was recently introduced in the latest generation of system, capable of arbitrarily shaped treatment field. The purpose of this study is to develop a model based dose calculation algorithm to better handle the lateral scatter in an irregularly shaped small field for the CyberKnife system. Methods: A pencil beam dose calculation algorithm widely used in linac based treatment planning system was modified. The kernel parameters and intensity profile were systematically determined by fitting to the commissioning data. The model was tuned using only a subset of measured data (4 out of 12 cones) and applied to all fixed circular cones for evaluation. The root mean square (RMS) of the difference between the measured and calculated tissue-phantom-ratios (TPRs) and off-center-ratio (OCR) was compared. Three cone size correction techniques were developed to better fit the OCRs at the penumbra region, which are further evaluated by the output factors (OFs). The pencil beam model was further validated against measurement data on the variable dodecagon-shaped Iris collimators and a half-beam blocked field. Comparison with Ray-Tracing and Monte Carlo methods was also performed on a lung SBRT case. Results: The RMS between the measured and calculated TPRs is 0.7% averaged for all cones, with the descending region at 0.5%. The RMSs of OCR at infield and outfield regions are both at 0.5%. The distance to agreement (DTA) at the OCR penumbra region is 0.2 mm. All three cone size correction models achieve the same improvement in OCR agreement, with the effective source shift model (SSM) preferred, due to their ability to predict more accurately the OF variations with the source to axis distance (SAD). In noncircular field validation, the pencil beam calculated results agreed well with the film measurement of both Iris collimators and the half-beam blocked field, fared much better than the Ray-Tracing calculation. Conclusions: The authors have developed a pencil beam dose calculation model for the CyberKnife system. The dose calculation accuracy is better than the standard linac based system because the model parameters were specifically tuned to the CyberKnife system and geometry correction factors. The model handles better the lateral scatter and has the potential to be used for the irregularly shaped fields. Comprehensive validations on MLC equipped system are necessary for its clinical implementation. It is reasonably fast enough to be used during plan optimization.

60 APPLIED LIFE SCIENCES↗

ROP model for PDC bits in Geothermal drilling

Geothermal energy is a renewable source of energy, where heat extraction is preferentially balanced with the reservoir's natural heat recharge rate. The objective of this paper is to present and validate a novel rate of penetration (ROP) model for drilling hard and abrasive formations including granite formations for polycrystalline diamond compact (PDC) bits. The ROP model was developed based on a derived relationship of a threshold weight on cutter (WOC) and its corresponding depth of cut (DOC) for a single cutter. Laboratory data was used to scale the derived single cutter relationship to a full-hole ROP model for PDC bits. The ROP model includes a non-linear correlation for Phase I (inefficient drilling due to low WOB values) and a linear Phase II (efficient drilling) ROP response to WOB. The ROP model was verified using measured drilling parameter data from Utah FORGE well# 58-32 and data from Chocolate Mountains well # 17-8 in Southern California. When compared to oil and gas well drilling, geothermal drilling in granitic formations can be more difficult and complicated due to rock hardness and high temperatures. PDC bits can increase ROP and optimize drilling for these types of hard formations. This paper provides novel insight into the ROP response of PDC bits to drilling operational parameters.

15 GEOTHERMAL ENERGY↗

A Predictive Transient Model of the TREAT-SIRIUS Experiments

To support experiments modeling and transient analysis of the NASA-sponsored SIRIUS experiments in the Transient Reactor Test Facility (TREAT), an innovative predictive transient model has been developed to simulate various transient experiments. The developed model is based on a new approach that utilizes steady-state Monte Carlo calculations along with a surrogate model based on polynomial regression to determine the core reactivity. The reactivity is supplied to an exact point kinetics model that determines the reactor power, and it is coupled to an adiabatic feedback model to determine the fuel average temperature. In this work, the newly developed model is introduced along with validation test results considering SIRIUS-1 experiments at different power levels. The initial test results show a very good agreement with experimental data at different power levels, which indicates that the model can be used to predict the power of other transient tests.

42 ENGINEERING↗

cabbi-bio/BECT-ABM

The emergence of the cellulosic bioeconomy requires not only adequate technological, economic, and policy advances, but also effective communication and coordination among the multiple stakeholders in the bioeconomy community. Aiming to facilitate stakeholder communication and collaboration for cellulosic bioeconomy growth, an agent-based model is developed as a community communication tool in this study. The tool simulates the behaviors and interactions of key stakeholders, as well as the responses and interactions between the stakeholders and biophysical, socioeconomic, and environmental conditions. The agent behavior rules are developed and validated based on behavior and economic theories and information attained from a stakeholder focus group study and a land use survey.

Yang, Pan↗

AI-Batt (Autonomous Identification of Battery Life Models) [SWR 21-36]

Autonomous Identification of Battery Life Models (AI-Batt) AI-Batt is a MATLAB code base for developing lifetime models for batteries from accelerated aging data. The code base provides many functions for processing, visualizing, and modeling battery aging data, making the data processing, exploration, and modeling workflow substantially faster. These tools are tailored for working with battery aging data sets, which usually consist of many separate time-series for each cell, with many test conditions and possible replicates at each condition, which makes it difficult to simply process or visualize the data set. Complex modeling tasks, such as cross-validation, sensitivity analysis, and uncertainty quantification have been implemented to enable thorough statistical investigation of model predictions. Additionally, several machine-learning algorithms are implemented to autonomously identify suitable models via symbolic regression. Data processing functions automatically cast data from the struct data type, which is commonly used to store experimental data, but is not an acceptable input for most algorithms, to the table data type, which can be easily used as input to any optimization algorithm. Also, the data can be separated into time-invariant and time-variant data tables, which is helpful for exploring the data set as well as developing separate models for time-variant and time-invariant aging mechanisms. For example, in aging tests with constant temperature, temperature is a time-invariant experimental condition. Visualization tools enable plotting of data, model fits, and model simulations possible with single-line function calls, empowering data exploration of complex data sets with both time-varying and time-invariant trends. Plots can be automatically generated for the whole data set, or separated by data group (groups of test replicates) or individual data series. Data points or data series can be automatically colored by the value of a variable with a variety of color maps, and model predictions can also be colored by the value of a fit statistic. Comparisons between data sets and the predictions/simulations of different models on the same data set can be easily plotted as well. Distributions of parameter values from bootstrap resampling can be plotted to visualize the reliability of parameter estimation, or determine any correlations between parameters. Modeling tools handle the complex task of creating and parsing symbolic equations for modeling battery lifetime. Equations are parsed to grab relevant data variables, parameter values, or specified sub-models for input into optimization, evaluation, or simulation functions. Models can be optimized locally (one set of parameters for each data series), bi-level (some parameters shared across the data set), or globally (single set of parameters for all data). Functions implementing symbolic regression algorithms help users to discover effective model equations, even in poorly sampled, high-dimensional data.

Smith, Kandler↗

Fluid structure interaction modelling for the vibration of tube bundles, part I: analysis of the fluid flow in a tube bundle

It is well known that a fluid may strongly influence the dynamic behaviour of a structure. Many different physical phenomena may take place, depending on the conditions: fluid flow, fluid at rest, little or high displacements of the structure. Inertial effects can take place, with lower vibration frequencies, dissipative effects also, with damping, instabilities due to the fluid flow (Fluid Induced Vibration). In this last case the structure is excited by the fluid. Tube bundles structures are very common in the nuclear industry. The reactor cores and the steam generators are both structures immersed in a fluid which may be submitted to a seismic excitation or an impact. In this case the structure moves under an external excitation, and the movement is influence by the fluid. The main point in such system is that the geometry is complex, and could lead to very huge sizes for a numerical analysis. Homogenization models have been developed based on the Euler equations for the fluid. Only inertial effects are taken into account. A next step in the modelling is to build models based on the homogenization of the Navier-Stokes equations. The papers presents results on an important step in the development of such model: the analysis of the fluid flow in a oscillating tube bundle. The analysis are made from the results of simulations based on the Navier-Stokes equations for the fluid. Comparisons are made with the case of the oscillations of a single tube, for which a lot of results are available in the literature. Different fluid flow pattern may be found, depending in the Reynolds number (related to the velocity of the bundle) and the Keulegan Carpenter number (related to the displacement of the bundle). A special attention is paid to the quantification of the inertial and dissipative effects, and to the forces exchanges between the bundle and the fluid. The results of such analysis will be used in the building of models based on the homogenization of the Navier-Stokes equations for the fluid. (authors)

42 ENGINEERING↗

Rate Theory Modeling and Simulation of Silicide Fuel at LWR Conditions

As a promising candidate for the accident tolerant fuel (ATF) used in light water reactors (LWRs), the fuel performance of uranium silicide (U 3 Si 2 ) at LWR conditions needs to be well understood. In this report, rate theory model was developed based on existing experimental data and density functional theory (DFT) calculations so as to predict the fission gas behavior in U 3 Si 2 at LWR conditions. The fission gas behavior of U 3 Si 2 can be divided into three temperature regimes. During steady-state operation, the majority of the fission gas stays in intragranular bubbles, whereas the dominance of intergranular bubbles and fission gas release only occurs beyond 1000 K. The steady-state rate theory model was also used as reference to establish a gaseous swelling correlation of U 3 Si 2 for the BISON code. Meanwhile, the overpressurized bubble model was also developed so that the fission gas behavior at LOCA can be simulated. LOCA simulation showed that intragranular bubbles are still dominant after a 70 second LOCA, resulting in a controllable gaseous swelling. The fission gas behavior of U 3 Si 2 at LWR conditions is benign according to the rate theory prediction at both steady-state and LOCA conditions, which provides important references to the qualification of U 3 Si 2 as a LWR fuel material with excellent fuel performance and enhanced accident tolerance.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Development of an Agent-based Model to Analyze Contemporary Helium Markets

Although U.S. helium demand has remained relatively flat since 2009, exports of helium have increased significantly since then, driven primarily by demand for electronic and semiconductor manufacturing in Asia. In the midst of this global demand shift, the Helium Act dictates a new procedure for pricing and distributing the gas through a reserve that historically functioned as a loose “oligarchy.” The new procedure requires prices to be determined by the open market through auctions and a survey of market prices, as opposed to increasing prices according to the consumer price index. Response to these changes has caused temporary shortages, price increases, and a significant increase in the development of the helium extraction technologies used to produce helium from formerly marginal sources. Technologies are being developed and refined to extract helium from formerly low-yielding natural gas fields containing much lower amounts of helium than the previously considered economic threshold of 0.3%. Combining these transformative policies with the potential for new and significant global supplies from Qatar, Algeria, and Russia could lead to new and unforeseen market behaviors and reactions from global helium markets. The objective of the project is to analyze the global helium markets.

99 GENERAL AND MISCELLANEOUS↗

Empirical Dynamic Modeling Reveals Complexity of Methane Fluxes in a Temperate Salt Marsh

Abstract Methane dynamics within salt marshes are complex because vegetation types, temperature, oscillating water levels, and changes in salinity and redox conditions influence CH 4 production, consumption, oxidation, and emissions. These non‐linear and complex interactions among variables affect the traditionally expected functional relationships and present challenges for interpreting and developing process‐based models. We employed empirical dynamic modeling (EDM) and convergent cross mapping (CCM) as a novel approach for characterizing seasonal/multiday and diurnal CH 4 dynamics by inferring causal variables, lags, and interconnections among multiple biophysical variables within a temperate salt marsh using 5 years of eddy covariance data. EDM/CCM is a nonparametric approach capable of quantifying the coupling between variables while determining time scales where variable interactions are the most relevant. We found that gross primary productivity, tidal creek dissolved oxygen, and temperature were important for seasonal/multiday dynamics (rho = 0.73–0.80), while water level was most important for diurnal dynamics during both the growing and dormancy phenoperiods (rho = 0.72 and 0.56, respectively). Lags for the top‐ranked variables (i.e., gross primary productivity, dissolved oxygen, temperature, water level) occurred between 1 and 5 weeks at the seasonal scale and 1–24 hr at the diurnal scale. The EDM had high prediction capabilities for intra‐/inter‐seasonal patterns and annual CH 4 sums but had limitations in representing large, infrequent fluxes. Results highlight the importance of non‐linearity, drivers, lag times, and interconnections among multiple biophysical variables that regulate CH 4 fluxes in tidal wetlands. This research introduces a novel approach to examining CH 4 fluxes, which will aid in evaluating current paradigms in wetlands and other ecosystems.

Environmental Sciences & Ecology↗

Model-based optimizations of packaged rooftop air conditioners using low global warming potential refrigerants

Based on laboratory investigations for R-22 and R-410A alternative low GWP refrigerants in two baseline rooftop air conditioners (RTU), the DOE/ORNL Heat Pump Design Model was used to model the two RTUs and the models were calibrated against the experimental data. We compared the compressor efficiencies and heat exchanger performances. An efficiency-based compressor mapping method was developed. Extensive model-based optimizations were conducted to provide a fair comparison between all the low GWP candidates by selecting optimal configurations. The results illustrate that all the R-22 low GWP refrigerants will lead to slightly lower COPs. ARM-20B appears to be the best R-22 replacement at normal conditions. At higher ambient temperatures, ARM-20A exhibits better performance. All R-410A low GWP candidates will result in similar or better efficiencies than R-410A. R-32 has the best COP while requiring the smallest compressor. Finally, R-452B uses the closest compressor displacement volume and achieves the same efficiency as R-410A.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Development, Verification, and Validation of Multiphase Models for Polydisperse Flows

This report describes in detail the technical findings of the DOE Award entitled 'Development, Verification, and Validation of Multiphase Models for Polydisperse Flows.' The focus was on high-velocity, gas-solid flows with a range of particle sizes. A complete mathematical model was developed based on first principles and incorporated into MFIX. The solid-phase description took two forms: the Kinetic Theory of Granular Flows (KTGF) and Discrete Quadrature Method of Moments (DQMOM). The gas-solid drag law for polydisperse flows was developed over a range of flow conditions using Discrete Numerical Simulations (DNS). These models were verified via examination of a range of limiting cases and comparison with Discrete Element Method (DEM) data. Validation took the form of comparison with both DEM and experimental data. Experiments were conducted in three separate circulating fluidized beds (CFB's), with emphasis on the riser section. Measurements included bulk quantities like pressure drop and elutriation, as well as axial and radial measurements of bubble characteristics, cluster characteristics, solids flux, and differential pressure drops (axial only). Monodisperse systems were compared to their binary and continuous particle size distribution (PSD) counterparts. The continuous distributions examined included Gaussian, lognormal, and NETL-provided data for a coal gasifier.

01 COAL, LIGNITE, AND PEAT↗

Reliability of emerging bonded interface materials for large-area attachments

In this study, conventional thermal interface materials (TIMs), such as greases, gels, and phase change materials, pose bottlenecks to heat removal and have long caused reliability issues in automotive power electronics packages. Bonded interface materials (BIMs) with superior thermal performance have the potential to be a replacement to the conventional TIMs. However, due to coefficient of thermal expansion mismatches between different components in a package and resultant thermomechanical stresses, fractures or delamination could occur, causing serious reliability concerns. These defects manifest themselves in increased thermal resistance in the package. In this paper, the results of reliability evaluation of emerging BIMs for large-area attachments in power electronics packaging are reported. Thermoplastic (polyamide) adhesive with embedded near-vertical-aligned carbon fibers, sintered silver, and conventional lead solder (Sn 63 Pb 37 ) materials were bonded between 50.8 mm x 50.8 mm cross-sectional footprint silicon nitride substrates and copper base plate samples, and were subjected to accelerated thermal cycling until failure or 2500 cycles. Damage in the BIMs was monitored every 100 cycles by scanning acoustic microscopy. Thermoplastic with embedded carbon fibers performed the best with no defects, whereas sintered silver and lead solder failed at 2300 and 1400 thermal cycles, respectively. Besides thermal cycling, additional lead solder samples were subjected to thermal shock and thermal cycling with extended dwell periods. A finite element method (FEM)-based model was developed to simulate the behavior of lead solder under thermomechanical loading. Strain energy density per cycle results were calculated from the FEM simulations. A predictive lifetime model was formulated for lead solder by correlating strain energy density results extracted from modeling with cycles-to-failure obtained from experimental accelerated tests. A power-law-based approach was used to formulate the - redictive lifetime model.

30 DIRECT ENERGY CONVERSION↗

Integrated Computational Materials Engineering (ICME) for Third Generation Advanced High-Strength Steel Development

This paper presents an overview of a four-year project focused on development of an integrated computational materials engineering (ICME) toolset for third generation advanced high-strength steels (3GAHSS). Following a brief look at ICME as an emerging discipline within the Materials Genome Initiative, technical tasks in the ICME project will be discussed. Specific aims of the individual tasks are multi-scale, microstructure-based material model development using state-of-the-art computational and experimental techniques, forming, toolset assembly, design optimization, integration and technical cost modeling. The integrated approach is initially illustrated using a 980 grade transformation induced plasticity (TRIP) steel, subject to a two-step quenching and partitioning (Q&P) heat treatment, as an example.

Savic, Vesna↗

Dynamic Differential Image Circle Diameter Measurement Precision Assessment: Application to Burning Droplets

Dynamic measurement precision assessment has been achieved for a differential circle measurement application. Differential circle diameter measurement, in image analysis, typically requires fitting a circle model that optimizes for image distortions, defects or occlusions. The differential task occurs when precise measurements of diameter change are required given object size variation with time. An automated system was designed to provide diameter measurements and associated measurement precision of images of a fuel droplet undergoing combustion in zero gravity for the FLEX-2 dataset. An image gradient-based, least-squares boundary point fitting method to a circle or ellipse model is used for diameter measurement. The presence of soot aggregates poses significant challenges for diameter measurements when it occludes part of the droplet boundary. The precision of the diameter measurements depends upon the image quality. Using synthetic image simulations that model the soot behavior, we developed a model based on image quality measures that assesses the measurement precision for each individual diameter measurement. Thus, diameter measurements with precision assessments were made available for follow-up scientific analysis. As a result, the algorithm's success rate for measurable runs was 98%. In cases of limited occlusion, a measurement precision of ±0.2 pixels for the FLEX-2 dataset was achieved.

42 ENGINEERING↗

Low–Cost Bio-Based Carbon Fiber for High-Temperature Processing

GrafTech International Holdings Inc. (GTI), worked with Oak Ridge National Laboratory (ORNL) under CRADA No. NFE-15-05807 to develop lignin-based carbon fiber (LBCF) technology and to demonstrate LBCF performance in high-temperature products and applications. This work was unique and different from other reported LBCF work in that this study was application-focused and scalability-focused. Accordingly, the executed work was based on meeting criteria based on technology development, cost, and application suitability. The focus of this work was to demonstrate lab-scale LBCF from at least 4 different precursor feedstock sources that could meet the estimated production cost of $5.00/pound and have ash level of less than 500 ppm in the carbonized insulation-grade fiber. Accordingly, a preliminary cost model was developed based on publicly available information. The team demonstrated that 4 lignin samples met the cost criteria, as highlighted in Table 1. In addition, the ash level for the 4 carbonized lignin samples were below 500 ppm. Processing asreceived lignin to produce a high purity lignin fiber was a significant accomplishment in that most industrial lignin, prior to purification, had greater than 4X the ash level needed for this project, and prior to this work there was not a clear path of how to achieve the purity target. The lab scale development of LBCF was performed with a specific functional application in mind, specifically for high temperature rigid insulation. GTI is currently a consumer of foreignsourced pitch and rayon based carbon fibers for use in its high temperature insulation products, and the motivation was that LBCF had potential to decrease costs and increase product competitiveness in the marketplace through lowered raw material costs, lowered energy costs, and decreased environmental footprint. At the end of this project, the Technology Readiness Level (TRL) remained at 5 for LBCF in high temperature insulation.

36 MATERIALS SCIENCE↗

Building America Case Study: Southern Energy Homes, First DOE Zero Energy Ready Manufactured Home, Russellville, Alabama

Three side-by-side lab houses were built, instrumented and monitored in an effort to determine through field testing and analysis the relative contributions of select technologies toward reducing energy use in new manufactured homes.The lab houses in Russellville, Alabama compared the performance of three homes built to varying levels of thermal integrity and HVAC equipment: a baseline HUD-code home equipped with an electric furnace and a split system air conditioner; an ENERGY STAR manufactured home with an enhanced thermal envelope and traditional split system heat pump; and a house designed to qualify for Zero Energy Ready Home designation with a ductless mini-split heat pump with transfer fan distribution system in place of the traditional duct system for distribution. Experiments were conducted in the lab houses to evaluate impact on energy and comfort of interior door position, window blind position and transfer fan operation. The report describes results of tracer gas and co-heating tests and presents calculation of the heat pump coefficient of performance for both the traditional heat pump and the ductless mini-split. A series of calibrated energy models was developed based on measured data and run in three locations in the Southeast to compare annual energy usage of the three homes.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Thermal creep model for CWSR zircaloy-4 cladding taking into account the annealing of the irradiation hardening

After irradiation and cooling in a pool, spent nuclear fuel assemblies are either transported for wet storage to a devoted site or loaded in casks for dry storage. During dry transportation or at the beginning of dry storage, the cladding is expected to be submitted to creep deformation under the hoop stress induced by the internal pressure of the fuel rod. The thermal creep is a potential mechanism that might lead to cladding failure. A new creep model was developed, based on a database of creep tests on as-received and irradiated cold-worked stress-relieved Zircaloy-4 cladding in a wide range of temperatures (310 degrees C to 470 degrees C) and hoop stress (80 to 260 MPa). Based on three laws-a flow law, a strain-hardening recovery law, and an annealing of irradiation hardening law this model allows the simulation of not only the transient creep and the steady-state creep, but also the early creep acceleration observed on irradiated samples tested in severe conditions, which was not taken into account in the previous models. The extrapolation of the creep model in the conditions of very long-term creep tests is reassuring, proving the robustness of the chosen formalism. The creep model has been assessed in progressively decreasing stress conditions, more representative of a transport. Set up to predict the cladding creep behavior under variable temperature and stress conditions, this model can easily be implemented into codes in order to simulate the thermomechanical behavior of spent fuel rods in various scenarios of postirradiation phases. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗