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 829 records · Page 46

Strongly nonlinear wave propagation in elasto-plastic metamaterials: Low-order dynamic modeling

Nonlinear elastic metamaterials are known to support a variety of dynamic phenomena that enhance our capacity to manipulate elastic waves. Since these properties stem from complex, subwavelength geometry, full-scale dynamic simulations are often prohibitively expensive at scales of interest. Prior studies have therefore utilized low-order effective medium models, such as discrete mass-spring lattices, to capture essential properties in the long-wavelength limit. While models of this type have been successfully implemented for a wide variety of nonlinear elastic systems, they have predominantly considered dynamics depending only on the instantaneous kinematics of the lattice, neglecting history-dependent effects, such as wear and plasticity. Here, to address this limitation, the present study develops a lattice-based modeling framework for nonlinear elastic metamaterials undergoing plastic deformation. Due to the history- and rate-dependent nature of plasticity, the framework generally yields a system of differential-algebraic equations whose computational cost is significantly greater than an elastic system of comparable size. We demonstrate the method using several models inspired by classical lattice dynamics and continuum plasticity theory and explore means to obtain empirical plasticity models for general geometries, thereby gaining insight into the influence of microstructural plasticity on effective material performance, which can be used to improve the design of nonlinear mechanical metamaterials.

Dynamic simulation↗

Development of Application Function Blocks for Power-Hardware-in-the-Loop Testing of Grid-Connected Inverters: Preprint

Stability and accuracy of the PHIL platform has been a common problem and challenging issue for a long time, thus continuous efforts and research have been devoted to developing interface algorithms to improve upon this. Some of existing works exhibit promising solutions to compensate the errors in the PHIL platform. However, there is a need to harmonize various solutions into a standard framework that provides a general solution for various PHIL test beds and that can be reused with minimum modifications to avoid extra effort. A good way to do this is by using application function blocks (AFBs), which are model-based and can be developed as a standard library component, thus achieving satisfactory reusability and reconfigurability for various platforms. Since PHIL simulation has the flexibility of the software simulation, proper functions can be easily implemented in the digital real-time simulator (DRTS) to preprocess the reference signal before sending it to the power amplifier. Therefore, this paper aims to develop standard AFBs in DRTS for the PHIL interface to test grid-connected inverters. The PHIL interface algorithm is separated into five AFBs that will be presented in detail in this paper.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A Causal Inference Model Based on Random Forests to Identify the Effect of Soil Moisture on Precipitation

Soil moisture influences precipitation mainly through its impact on land–atmosphere interactions. Understanding and correctly modeling soil moisture–precipitation (SM–P) coupling is crucial for improving weather forecasting and subseasonal to seasonal climate predictions, especially when predicting the persistence and magnitude of drought. However, the sign and spatial structure of SM–P feedback are still being debated in the climate research community, mainly due to the difficulty in establishing causal relationships and the high degree of nonlinearity in land–atmosphere processes. To this end, we developed a causal inference model based on the Granger causality analysis and a nonlinear machine learning model. This model includes three steps: nonlinear anomaly decomposition, nonlinear Granger causality analysis, and evaluation of the quality of SM–P feedback, which eliminates the nonlinear response of interannual and seasonal variability and the memory effects of climatic factors and isolates the causal relationship of local SM–P feedback. We applied this model by using National Climate Assessment–Land Data Assimilation System (NCA-LDAS) datasets over the United States. Here, the results highlight the importance of nonlinear atmosphere responses in land–atmosphere interactions. In addition, the strong feedback over the southwestern United States and the Great Plains both highlight the impacts of topographic factors rather than only the sensitivity of evapotranspiration to soil moisture. Furthermore, the SM–P index defined by our framework is used to benchmark Earth system models (ESMs), which provides a new metric for efficiently identifying potential model biases in modeling local land–atmosphere interactions and may help the development of ESMs in improving simulations of water cycle variability and extremes.

54 ENVIRONMENTAL SCIENCES↗

Modeling the effects of carbon payments and forest owner cooperatives on carbon storage and revenue in Pacific Northwest forestlands

The forestlands of the Pacific Northwest store a significant amount of carbon, provide wildlife habitat, generate revenue for forestland owners, and provide for sustainable generation of forest products. Promoting forest carbon storage is increasingly recognized as a key management objective to mitigate climate change. We developed an agent-based model to examine the effectiveness of current land management practices in sustaining forest carbon storage and revenue flows, and the utility of market-based carbon payments in improving carbon storage while also improving revenue generation in the Pacific Northwest. We find that while market based incentives can improve both forest carbon storage and revenue, they are limited in their benefit to small forest owners. Our results highlight the crucial need for additional policy interventions to reduce barriers to market access for small private owners and ensure sustainable forest management in this region. We further show that a forest owners’ cooperative, which aggregates small owners to reduce the costs of participating in a carbon market, can help overcome these barriers and result in both greater carbon storage and a more equitable benefit from carbon payments.

54 ENVIRONMENTAL SCIENCES↗

SU-E-T-578: MCEBRT, A Monte Carlo Code for External Beam Treatment Plan Verifications

Purpose: Present a new Monte Carlo code (MCEBRT) for patient-specific dose calculations in external beam radiotherapy. The code MLC model is benchmarked and real patient plans are re-calculated using MCEBRT and compared with commercial TPS. Methods: MCEBRT is based on the GEPTS system (Med. Phys. 29 (2002) 835–846). Phase space data generated for Varian linac photon beams (6 – 15 MV) are used as source term. MCEBRT uses a realistic MLC model (tongue and groove, rounded ends). Patient CT and DICOM RT files are used to generate a 3D patient phantom and simulate the treatment configuration (gantry, collimator and couch angles; jaw positions; MLC sequences; MUs). MCEBRT dose distributions and DVHs are compared with those from TPS in absolute way (Gy). Results: Calculations based on the developed MLC model closely matches transmission measurements (pin-point ionization chamber at selected positions and film for lateral dose profile). See Fig.1. Dose calculations for two clinical cases (whole brain irradiation with opposed beams and lung case with eight fields) are carried out and outcomes are compared with the Eclipse AAA algorithm. Good agreement is observed for the brain case (Figs 2-3) except at the surface where MCEBRT dose can be higher by 20%. This is due to better modeling of electron contamination by MCEBRT. For the lung case an overall good agreement (91% gamma index passing rate with 3%/3mm DTA criterion) is observed (Fig.4) but dose in lung can be over-estimated by up to 10% by AAA (Fig.5). CTV and PTV DVHs from TPS and MCEBRT are nevertheless close (Fig.6). Conclusion: A new Monte Carlo code is developed for plan verification. Contrary to phantombased QA measurements, MCEBRT simulate the exact patient geometry and tissue composition. MCEBRT can be used as extra verification layer for plans where surface dose and tissue heterogeneity are an issue.

60 APPLIED LIFE SCIENCES↗

WE-EF-207-03: Design and Optimization of a CBCT Head Scanner for Detection of Acute Intracranial Hemorrhage

Purpose: To design a dedicated x-ray cone-beam CT (CBCT) system suitable to deployment at the point-of-care and offering reliable detection of acute intracranial hemorrhage (ICH), traumatic brain injury (TBI), stroke, and other head and neck injuries. Methods: A comprehensive task-based image quality model was developed to guide system design and optimization of a prototype head scanner suitable to imaging of acute TBI and ICH. Previously reported models were expanded to include the effects of x-ray scatter correction necessary for detection of low contrast ICH and the contribution of bit depth (digitization noise) to imaging performance. Task-based detectablity index provided the objective function for optimization of system geometry, x-ray source, detector type, anti-scatter grid, and technique at 10–25 mGy dose. Optimal characteristics were experimentally validated using a custom head phantom with 50 HU contrast ICH inserts imaged on a CBCT imaging bench allowing variation of system geometry, focal spot size, detector, grid selection, and x-ray technique. Results: The model guided selection of system geometry with a nominal source-detector distance 1100 mm and optimal magnification of 1.50. Focal spot size ∼0.6 mm was sufficient for spatial resolution requirements in ICH detection. Imaging at 90 kVp yielded the best tradeoff between noise and contrast. The model provided quantitation of tradeoffs between flat-panel and CMOS detectors with respect to electronic noise, field of view, and readout speed required for imaging of ICH. An anti-scatter grid was shown to provide modest benefit in conjunction with post-acquisition scatter correction. Images of the head phantom demonstrate visualization of millimeter-scale simulated ICH. Conclusions: Performance consistent with acute TBI and ICH detection is feasible with model-based system design and robust artifact correction in a dedicated head CBCT system. Further improvements can be achieved with incorporation of model-based iterative reconstruction techniques also within the scope of the task-based optimization framework. David Foos and Xiaohui Wang are employees of Carestream Health.

60 APPLIED LIFE SCIENCES↗

A Monte Carlo calculation model of electronic portal imaging device for transit dosimetry through heterogeneous media

Purpose: To develop and evaluate a fast Monte Carlo (MC) dose calculation model of electronic portal imaging device (EPID) based on its effective atomic number modeling in the XVMC code. Methods: A previously developed EPID model, based on the XVMC code by density scaling of EPID structures, was modified by additionally considering effective atomic number (Z{sub eff}) of each structure and adopting a phase space file from the EGSnrc code. The model was tested under various homogeneous and heterogeneous phantoms and field sizes by comparing the calculations in the model with measurements in EPID. In order to better evaluate the model, the performance of the XVMC code was separately tested by comparing calculated dose to water with ion chamber (IC) array measurement in the plane of EPID. Results: In the EPID plane, calculated dose to water by the code showed agreement with IC measurements within 1.8%. The difference was averaged across the in-field regions of the acquired profiles for all field sizes and phantoms. The maximum point difference was 2.8%, affected by proximity of the maximum points to penumbra and MC noise. The EPID model showed agreement with measured EPID images within 1.3%. The maximum point difference was 1.9%. The difference dropped from the higher value of the code by employing the calibration that is dependent on field sizes and thicknesses for the conversion of calculated images to measured images. Thanks to the Z{sub eff} correction, the EPID model showed a linear trend of the calibration factors unlike those of the density-only-scaled model. The phase space file from the EGSnrc code sharpened penumbra profiles significantly, improving agreement of calculated profiles with measured profiles. Conclusions: Demonstrating high accuracy, the EPID model with the associated calibration system may be used for in vivo dosimetry of radiation therapy. Through this study, a MC model of EPID has been developed, and their performance has been rigorously investigated for transit dosimetry.

60 APPLIED LIFE SCIENCES↗

THREE-DIMENSIONAL RADIO AND X-RAY MODELING AND DATA ANALYSIS SOFTWARE: REVEALING FLARE COMPLEXITY

Many problems in solar physics require analysis of imaging data obtained in multiple wavelength domains with differing spatial resolution in a framework supplied by advanced three-dimensional (3D) physical models. To facilitate this goal, we have undertaken a major enhancement of our IDL-based simulation tools developed earlier for modeling microwave and X-ray emission. The enhanced software architecture allows the user to (1) import photospheric magnetic field maps and perform magnetic field extrapolations to generate 3D magnetic field models; (2) investigate the magnetic topology by interactively creating field lines and associated flux tubes; (3) populate the flux tubes with user-defined nonuniform thermal plasma and anisotropic, nonuniform, nonthermal electron distributions; (4) investigate the spatial and spectral properties of radio and X-ray emission calculated from the model; and (5) compare the model-derived images and spectra with observational data. The package integrates shared-object libraries containing fast gyrosynchrotron emission codes, IDL-based soft and hard X-ray codes, and potential and linear force-free field extrapolation routines. The package accepts user-defined radiation and magnetic field extrapolation plug-ins. We use this tool to analyze a relatively simple single-loop flare and use the model to constrain the magnetic 3D structure and spatial distribution of the fast electrons inside this loop. We iteratively compute multi-frequency microwave and multi-energy X-ray images from realistic magnetic flux tubes obtained from pre-flare extrapolations, and compare them with imaging data obtained by SDO, NoRH, and RHESSI. We use this event to illustrate the tool's use for the general interpretation of solar flares to address disparate problems in solar physics.

79 ASTRONOMY AND ASTROPHYSICS↗

Neutron diffraction measurements and micromechanical modelling of temperature-dependent variations in TATB lattice parameters

Triaminotrinitrobenzene (TATB) is a highly anisotropic molecular crystal used in several plastic-bonded explosive (PBX) formulations. TATB-based explosives exhibit irreversible volume expansion (“ratchet growth”) when thermally cycled. A theoretical understanding of the relationship between anisotropy of the crystal, crystal orientation distribution (texture) of polycrystalline aggregates, and the intergranular interactions leading to this irreversible growth is necessary to accurately develop physics-based predictive models for TATB-based PBXs under various thermal environments. In this work, TATB lattice parameters were measured using neutron diffraction during thermal cycling of loose powder and a pressed pellet. The measured lattice parameters help clarify conflicting reports in the literature as these new results are more consistent with one set of previous results than another. The lattice parameters of pressed TATB were also measured as a function of temperature, showing some differences from the powder. This data is used along with anisotropic single-crystal stiffness moduli reported in the literature to model the nominal stresses associated with intergranular constraints during thermal expansion. The texture of both specimens were characterized and the pressed pellet exhibits preferential orientation of (001) poles along the pressing direction, whereas no preferred orientation was found for the loose powder. Lastly, thermal strains for single-crystal TATB computed from lattice parameter data for the powder is input to a self-consistent micromechanical model, which predicts the lattice parameters of the constrained TATB crystals within the pellet. The agreement of these model results with the diffraction data obtained from the pellet is discussed along with future directions of research.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Physical model and experimental validation of a high temperature proton exchange membrane electrochemical hydrogen pump cell for efficient single-stage extraction of low concentration hydrogen gas

There is interest in valorization of existing natural gas infrastructure to facilitate the co-transportation of hydrogen via blending of hydrogen gas initially at limited concentrations of 1–20 vol% H2 and to subsequently extract hydrogen at fuel cell quality standards (SAE J2719/ISO14687-2). High temperature proton exchange membrane electrochemical hydrogen pump (HT-PEM EHP) based on phosphoric acid doped polybenzimidazole (PA-PBI) exhibits good performance at elevated temperatures (>120 °C), which provides desirable tolerance to non-methane natural gas constituents that are problematic for lower temperature based EHP. To better understand the suitability of the HT-PEM EHP for such gas separation processes, a two-dimensional model of EHP based on PA-PBI was developed. The model is validated for several relevant operating conditions and across cells with differing amounts of phosphoric acid content in the electrodes. Operando micro x-ray computed tomography (CT) imaging of an HT-PEM EHP was used to further validate physical parameters and assumptions of the model. The impacts of pressure, relative humidity of the anode feed, and concentration of feed gases on separation performance are investigated. This study shows that a specific energy of separation of 5.1 kWh/kg H2 at a hydrogen recovery factor (HRF) of 50 vol% can be achieved in a single stage with the EHP, producing fuel cell quality hydrogen purity of 99.99 vol% H2 from a 2 vol% H2/CH4 feed blend, while pressurizing the product H2 at a pressure ratio of 1.3 relative to feed pressure.

Stansberry, John M↗

The impact of low-Z impurities on x-ray conversion efficiency from laser-produced plasmas of low-density gold foam targets

It is an important approach to improve the x-ray conversion efficiency of laser-ablated high-Z plasmas by using low initial density materials for various applications. However, unavoidable low-Z impurities in the manufacture process of low-density high-Z foam targets will depress this effect. A general easy-to-use analytical model based on simulations was developed to evaluate the quantitative impact of impurities within the gold foam target on laser to x-ray conversion efficiency. In addition, the x-ray conversion efficiencies of 1 g/cm{sup 3} gold foams with two different initial contents of impurities were experimentally investigated. Good agreements have been achieved between the model results and experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Shattered pellet injection simulations with NIMROD

Optimal strategies for disruption mitigation benefit from the understanding of details both spatially and temporally. Beyond the assessment of the efficacy of a particular proposed Disruption Mitigation System (DMS), ITER's longevity will require accounting of both mitigated and unmitigated disruptions. Accurate models and validated simulations that detail multiple ITER scenarios with mitigated and unmitigated disruptions are essential for accurate estimates of load damage. The primary candidate for ITER's DMS is Shattered Pellet Injection (SPI); its efficacy must be evaluated within the next several years. To perform critical time dependent 3-D nonlinear simulations, we have developed a particle based SPI model in the NIMROD code coupled to its modified single fluid equations with impurity and radiation [Izzo, Nucl. Fusion 46(5), 541 (2006)]. SPI validation simulations of the thermal quench and comparisons to DIII-D impurity scan experiments [Shiraki et al., Phys. Plasmas 23(6), 062516 (2016)] are presented. We also present an initial ITER Q = 10 pure neon SPI simulation and compare it with the DIII-D SPI simulations. NIMROD SPI simulations demonstrate that the ablating fragment drives strong parallel flows that transport the impurities and governs the thermal quench. Analysis of SPI simulations shows that the mixed deuterium/neon SPI results in a more benign thermal quench due to the enhanced transport caused by the additional deuterium. These results suggest that an optimal pellet mixture exists for the SPI system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An integrated QSAR-PBK/D modelling approach for predicting detoxification and DNA adduct formation of 18 acyclic food-borne α,β-unsaturated aldehydes

Acyclic α,β-unsaturated aldehydes present in food raise a concern because the α,β-unsaturated aldehyde moiety is considered a structural alert for genotoxicity. However, controversy remains on whether in vivo at realistic dietary exposure DNA adduct formation is significant. The aim of the present study was to develop physiologically based kinetic/dynamic (PBK/D) models to examine dose-dependent detoxification and DNA adduct formation of a group of 18 food-borne acyclic α,β-unsaturated aldehydes without 2- or 3-alkylation, and with no more than one conjugated double bond. Parameters for the PBK/D models were obtained using quantitative structure–activity relationships (QSARs) defined with a training set of six selected aldehydes. Using the QSARs, PBK/D models for the other 12 aldehydes were defined. Results revealed that DNA adduct formation in the liver increases with decreasing bulkiness of the molecule especially due to less efficient detoxification. 2-Propenal (acrolein) was identified to induce the highest DNA adduct levels. At realistic dietary intake, the predicted DNA adduct levels for all aldehydes were two orders of magnitude lower than endogenous background levels observed in disease free human liver, suggesting that for all 18 aldehydes DNA adduct formation is negligible at the relevant levels of dietary intake. The present study provides a proof of principle for the use of QSAR-based PBK/D modelling to facilitate group evaluations and read-across in risk assessment. - Highlights: • Physiologically based in silico models were made for 18 α,β-unsaturated aldehydes. • Kinetic parameters were determined by in vitro incubations and a QSAR approach. • DNA adduct formation was negligible at levels relevant for dietary intake. • The use of QSAR-based PBK/D modelling facilitates group evaluations and read-across.

60 APPLIED LIFE SCIENCES↗

Dynamic Modeling of Aerobic Growth of Shewanella oneidensis. Predicting Triauxic Growth, Flux Distributions and Energy Requirement for Growth

A model-based analysis is conducted to investigate metabolism of Shewanella oneidensis MR-1 strain in aerobic batch culture, which exhibits an intriguing growth pattern by sequentially consuming substrate (i.e., lactate) and by-products (i.e., pyruvate and acetate). A general protocol is presented for developing a detailed network-based dynamic model for S. oneidensis based on the Lumped Hybrid Cybernetic Model (LHCM) framework. The L-HCM, although developed from only limited data, is shown to accurately reproduce exacting dynamic metabolic shifts, and provide reasonable estimates of energy requirement for growth. Flux distributions in S. oneidensis predicted by the L-HCM compare very favorably with 13C-metabolic flux analysis results reported in the literature. Predictive accuracy is enhanced by incorporating measurements of only a few intracellular fluxes, in addition to extracellular metabolites. The L-HCM developed here for S. oneidensis is consequently a promising tool for the analysis of intracellular flux distribution and metabolic engineering.

ATP requirement for growth↗

Non-equilibrium approach to prediction of microstructure evolution for metals undergoing severe plastic deformation

Highlights: • New consistent physically-based model of dynamic recrystallization. • The model is based on evolution of fraction of high-angle grain boundaries. • Possibility to predict the evolution of dislocation cells and grain boundaries. • The model is implemented for FEM simulations. • Results demonstrate inhomogeneous distribution of strain within the sample. - Abstract: Available models of dynamic recrystallization have a number of disadvantages that in most cases make them inapplicable for practical predictions of material microstructure evolution. Both the microstructural and the empirically based approaches do not reflect physical processes leading to evolution of material defect structure in the process of plastic deformation. This work presents an attempt to develop a consistent physically-based model of dynamic recrystallization. This model, accounting for physical nature of processes of material defect structure evolution, should provide a possibility to predict evolution of several different experimentally measurable parameters of material microstructure without introduction of big number of fitting parameters. It is suggested that such a model should be based on equation for evolution of fraction of high-angle grain boundaries (HAGBs) in the process of deformation. It is shown, that the new model gives a possibility to predict the evolution of dislocation cells and grain boundaries in copper-based alloys providing good coincidence with experimental observations. Full 3-dimensional numerical simulation of multidirectional forging of copper is performed utilizing the developed dynamic recrystallization model. The same 3D simulations demonstrate new noteworthy effects connected to inhomogeneous distribution of plastic strain within the bulk of the material and material strain hardening.

36 MATERIALS SCIENCE↗

A thermodynamic and kinetic model for paste–aggregate interactions and the alkali–silica reaction

A new conceptual model is developed for ASR formation based on geochemical principles tied to aqueous speciation, silica solubility, kinetically controlled mineral dissolution, and diffusion. ASR development is driven largely by pH and silica gradients that establish geochemical microenvironments between paste and aggregate, with gradients the strongest within the aggregate adjacent to the paste boundary (i.e., where ASR initially forms). Super-saturation of magadiite and okenite (crystalline ASR surrogates) occurs in the zone defined by gradients in pH, dissolved silica, Na{sup +}, and Ca{sup 2} {sup +}. This model provides a thermodynamic rather than kinetic explanation of why quartz generally behaves differently from amorphous silica: quartz solubility does not produce sufficiently high concentrations of H{sub 4}SiO{sub 4} to super-saturate magadiite, whereas amorphous silica does. The model also explains why pozzolans do not generate ASR: their fine-grained character precludes formation of chemical gradients. Finally, these gradients have interesting implications beyond the development of ASR, creating unique biogeochemical environments.

36 MATERIALS SCIENCE↗

Development of a General Form CO 2 and Brine Flux Input Model

The National Risk Assessment Partnership (NRAP) project is developing a science-based toolset for the quantitative analysis of the potential risks associated with changes in groundwater chemistry from CO 2 injection. In order to address uncertainty probabilistically, NRAP is developing efficient, reduced-order models (ROMs) as part of its approach. These ROMs are built from detailed, physics-based process models to provide confidence in the predictions over a range of conditions. The ROMs are designed to reproduce accurately the predictions from the computationally intensive process models at a fraction of the computational time, thereby allowing the utilization of Monte Carlo methods to probe variability in key parameters. This report presents the procedures used to develop a generalized model for CO 2 and brine leakage fluxes based on the output of a numerical wellbore simulation. The resulting generalized parameters and ranges reported here will be used for the development of third-generation groundwater ROMs.

58 GEOSCIENCES↗

WTFgenes: What's The Function of these genes? Static sites for model-based gene set analysis

A common technique for interpreting experimentally-identified lists of genes is to look for enrichment of genes associated with particular ontology terms. The most common test uses the hypergeometric distribution; more recently, a model-based test was proposed. These approaches must typically be run using downloaded software, or on a server. We develop a collapsed likelihood for model-based gene set analysis and present WTFgenes, an implementation of both hypergeometric and model-based approaches, that can be published as a static site with computation run in JavaScript on the user's web browser client. Apart from hosting files, zero server resources are required: the site can (for example) be served directly from Amazon S3 or GitHub Pages. A C++11 implementation yielding identical results runs roughly twice as fast as the JavaScript version. WTFgenes is available from https://github.com/evoldoers/wtfgenes under the BSD3 license. A demonstration for the Gene Ontology is usable at https://evoldoers.github.io/wtfgo.

Mungall, Christopher J.↗