Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model analyses”

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

Shielding requirements for the Space Station habitability modules

The design, analysis, development, and tests of the total meteoroid/debris protection system for the Space Station Freedom habitability modules, such as the habitation module, the laboratory module, and the node structures, are described. Design requirements are discussed along with development efforts, including a combination of hypervelocity testing and analyses. Computer hydrocode analysis of hypervelocity impact phenomena associated with Space Station habitability structures is covered and the use of optimization techniques, engineering models, and parametric analyses is assessed. Explosive rail gun development efforts and protective capability and damage tolerance of multilayer insulation due to meteoroid/debris impact are considered. It is concluded that anticipated changes in the debris environment definition and requirements will require rescoping the tests and analysis required to develop a protection system.

Avans, Sherman L.↗

ISAC - A tool for aeroservoelastic modeling and analysis

This paper discusses the capabilities of the Interaction of Structures, Aerodynamics, and Controls (ISAC) system of program modules. The major modeling, analysis, and data management components of ISAC are identified. Equations of motion are displayed for a Laplace-domain representation of the unsteady aerodynamic forces. Options for approximating a frequency-domain representation of unsteady aerodynamic forces with rational functions of the Laplace variable are shown. Linear time invariant state-space equations of motion that result are discussed. Model generation and analyses of stability and dynamic response characteristics are shown for an aeroelastic vehicle which illustrate some of the capabilities of ISAC as a modeling and analysis tool for aeroelastic applications.

Adams, William M., Jr.↗

ISAC: A tool for aeroservoelastic modeling and analysis

The capabilities of the Interaction of Structures, Aerodynamics, and Controls (ISAC) system of program modules is discussed. The major modeling, analysis, and data management components of ISAC are identified. Equations of motion are displayed for a Laplace-domain representation of the unsteady aerodynamic forces. Options for approximating a frequency-domain representation of unsteady aerodynamic forces with rational functions of the Laplace variable are shown. Linear time invariant state-space equations of motion that result are discussed. Model generation and analyses of stability and dynamic response characteristics are shown for an aeroelastic vehicle which illustrates some of the capabilities of ISAC as a modeling and analysis tool for aeroelastic applications.

Adams, William M., Jr.↗

Dynamic genome evolution in a model fern

The large size and complexity of most fern genomes have hampered efforts to elucidate fundamental aspects of fern biology and land plant evolution through genome-enabled research. Here we present a chromosomal genome assembly and associated methylome, transcriptome and metabolome analyses for the model fern species Ceratopteris richardii. The assembly reveals a history of remarkably dynamic genome evolution including rapid changes in genome content and structure following the most recent whole-genome duplication approximately 60 million years ago. These changes include massive gene loss, rampant tandem duplications and multiple horizontal gene transfers from bacteria, contributing to the diversification of defence-related gene families. The insertion of transposable elements into introns has led to the large size of the Ceratopteris genome and to exceptionally long genes relative to other plants. Gene family analyses indicate that genes directing seed development were co-opted from those controlling the development of fern sporangia, providing insights into seed plant evolution. Our findings and annotated genome assembly extend the utility of Ceratopteris as a model for investigating and teaching plant biology.

59 BASIC BIOLOGICAL SCIENCES↗

Generating Performance Models for Irregular Applications

Many applications have irregular behavior --- non-uniform input data, input-dependent solvers, irregular memory accesses, unbiased branches --- that cannot be captured using today's automated performance modeling techniques. We describe new hierarchical critical path analyses for the \Palm model generation tool. To create a model's structure, we capture tasks along representative MPI critical paths. We create a histogram of critical tasks with parameterized task arguments and instance counts. To model each task, we identify hot instruction-level sub-paths and model each sub-path based on data flow, instruction scheduling, and data locality. We describe application models that generate accurate predictions for strong scaling when varying CPU speed, cache speed, memory speed, and architecture. We present results for the Sweep3D neutron transport benchmark; Page Rank on multiple graphs; Support Vector Machine with pruning; and PFLOTRAN's reactive flow/transport solver with domain-induced load imbalance.

analytical performance models↗

Do Coupled Climate Models Correctly Simulate the Upward Branch of the Deep Ocean Global Conveyor?

The large-scale meridional overturning circulation (MOC) connects the deep ocean, a major reservoir of carbon, to the other components of the climate system and must therefore be accurately represented in Earth System Models. Our project aims to address the specific question of the pathways and mechanisms controlling the upwelling branch of the MOC, a subject of significant disagreement between models and observational syntheses, and among general circulation models. Observations of these pathways are limited, particularly in regions of complex hydrography such as the Southern Ocean. As such, we rely on models to examine theories of the overturning circulation, both physically and biogeochemically. This grant focused on a particular aspect of the meridional overturning circulation (MOC) where there is currently significant disagreement between models and observationally based analyses of the MOC, and amongst general circulation models. In particular, the research focused on addressing the following questions: 1. Where does the deep water that sinks in the polar regions rise to the surface? 2. What processes are responsible for this rise? 3. Do state-of-the-art coupled GCMs capture these processes? Our research had three key components: observational synthesis, model development and model analysis. In this final report we outline the key results from these areas of research for the 2007 to 2012 grant period. The research described here was carried out primarily by graduate student, Daniele Bianchi (now a Postdoc at McGill University, Canada), and Postdoc Stephanie Downes (now a Research Fellow at The Australian national University, Australia). Additional support was provided for programmers Jennifer Simeon as well as Rick Slater.

3He data Synthesis↗

Distribution System V1G PEV Charging Impacts Report

Plug-in Electric Vehicles (PEVs) are a clean form of transportation that will require a new control paradigm between the energy supply and the transportation sectors as PEV adoption continues to grow. Meeting the projected PEV energy demands will require implementing PEV / grid integration methods beyond charging PEVs at incentivized times. MATLAB was used to simulate the daily vehicle trip demands and charging processes of hundreds of PEVs. The National Household Travel Survey (NHTS) data was used as the source of vehicle travel data, and EV Project provided measured PEV charging data to calibrate the PEV population models. GridLAB-D analyses used a prototypical feeder and PEV population models to evaluate residential PEV V1G (unidirectional) charging in multiple locations, V1G PEV distribution feeder impacts, Time-Of-Use (TOU) rate effects on PEV charging, and charging control methods. The key findings were: optimized charging using California TOU rates could result in savings of over $20 per month for non-TOU program participants; non-TOU participants consume 30-40% of PEV charging energy; about 50% of PEV drivers have their transportation fuel bill reduced by 50 dollars per month; 6.6kW charging will cause distribution transformer power limits to be exceeded on prototypical feeders if the household PEV adoption rate reaches 0.75 PEVs per home; distribution feeder loading is proportional to PEV charging rate; a method was devised to remotely identify when residential transformer overloads occur using changes in PEV charging rate; and the low-cost TOU period is composed of a distribution of PEVs needing small, medium and larger energy needs that can be more effectively shifted as groups to mitigate the TOU peak power and enable a higher PEV adoption rate per household without overloading distribution transformers.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Boundary layer models for calving marine outlet glaciers

Abstract. We consider the flow of marine-terminating outlet glaciers that are laterally confined in a channel of prescribed width. In that case, the drag exerted by the channel side walls on a floating ice shelf can reduce extensional stress at the grounding line. If ice flux through the grounding line increases with both ice thickness and extensional stress, then a longer shelf can reduce ice flux by decreasing extensional stress. Consequently, calving has an effect on flux through the grounding line by regulating the length of the shelf. In the absence of a shelf, it plays a similar role by controlling the above-flotation height of the calving cliff. Using two calving laws, one due to Nick et al. (2010) based on a model for crevasse propagation due to hydrofracture and the other simply asserting that calving occurs where the glacier ice becomes afloat, we pose and analyse a flowline model for a marine-terminating glacier by two methods: direct numerical solution and matched asymptotic expansions. The latter leads to a boundary layer formulation that predicts flux through the grounding line as a function of depth to bedrock, channel width, basal drag coefficient, and a calving parameter. By contrast with unbuttressed marine ice sheets, we find that flux can decrease with increasing depth to bedrock at the grounding line, reversing the usual stability criterion for steady grounding line location. Stable steady states can then have grounding lines located on retrograde slopes. We show how this anomalous behaviour relates to the strength of lateral versus basal drag on the grounded portion of the glacier and to the specifics of the calving law used.

54 ENVIRONMENTAL SCIENCES↗

Mobility Data and Models Informing Smart Cities

DOE SMART Mobility Urban Science Efforts are helping use/develop key data sets, models, and roles for DOE in engaging across seven Smart City Finalists as well as for ensuring useful/useable insights. Austin, Texas, was targeted for analyses of opportunity (model and data maturity). Other Urban Science/broader SMART initiatives were also helped by this research.

33 ADVANCED PROPULSION SYSTEMS↗

Numerical coupling of aerosol emissions, dry removal, and turbulent mixing in the E3SM Atmosphere Model version 1 (EAMv1) – Part 1: Dust budget analyses and the impacts of a revised coupling scheme

Abstract. An earlier study evaluating dust life cycle in the Energy Exascale Earth System Model (E3SM) Atmosphere Model version 1 (EAMv1) has revealed that the simulated global mean dust lifetime is substantially shorter when higher vertical resolution is used, primarily due to significant strengthening of dust dry removal in source regions. This paper demonstrates that the sequential splitting of aerosol emissions, dry removal, and turbulent mixing in the model's time integration loop, especially the calculation of dry removal after surface emissions and before turbulent mixing, is the primary reason for the vertical resolution sensitivity reported in that earlier study. Based on this reasoning, we propose a revised numerical process coupling scheme that requires the least amount of code changes, in which the surface emissions are applied before turbulent mixing instead of before dry removal. The revised scheme allows newly emitted particles to be transported aloft by turbulence before being removed from the atmosphere, and hence better resembles the dust life cycle in the real world. Sensitivity experiments show that the revised process coupling substantially weakens dry removal and strengthens vertical mixing in dust source regions. It also strengthens the large-scale transport from source to non-source regions, strengthens dry removal outside the source regions, and strengthens wet removal and activation globally. In transient simulations of the years 2000–2009 conducted using 1∘ horizontal grid spacing, 72 vertical layers, and unchanged tuning parameters of emission strength, the revised process coupling leads to a 40 % increase in the global total dust burden and an increase of dust lifetime from 1.8 to 2.5 d in terms of 10-year averages. Weakened dry removal and increased mixing ratios are also seen for other aerosol species that have substantial surface emissions, although the changes in mixing ratio are considerably smaller for the submicron species than for dust and sea salt. Numerical experiments confirm that the revised coupling scheme significantly reduces the strong and non-physical sensitivities of model results to vertical resolution in the original EAMv1. This provides a motivation for adopting the revised scheme in EAM as well as for further improvements on the simple revision presented in this paper.

54 ENVIRONMENTAL SCIENCES↗

Simple, accurate surrogate models of the elastic response of three-dimensional open truss micro-architectures with applications to multiscale topology design

Elastic meta-materials are those whose unique properties come from their micro-architecture, rather than, e.g., from their chemistry. The introduction of such architecture, which is increasingly able to be fabricated due to advances in additive manufacturing, expands the design domain and enables improved design, from the most complex multi-physics design problems to the simple compliance design problem that is our focus. Unfortunately, concurrent design of both the micro-scale and the macroscale is computationally very expensive when the former can vary spatially, particularly in three dimensions. Instead, we provide simple, accurate surrogate models of the homogenized linear elastic response of the isotruss, the octet truss, and the ORC truss based on high-fidelity continuum finite element analyses. These surrogate models are relatively accurate over the full range of relative densities, in contrast to analytical models in the literature, which we show lose accuracy as relative density increases. The surrogate models are also simple to implement, which we demonstrate by modifying Sigmund’s 99-line code to solve a three-dimensional, multiscale compliance design problem with spatially varying relative density. Here, we use this code to generate examples in both two and three dimensions that illustrate the advantage of elastic meta-materials over structures with a single length scale, i.e., those without micro-architectures.

42 ENGINEERING↗

An extension of Thwaites’ method for turbulent boundary layers

Thwaites developed an approximate method for determining the evolution of laminar boundary layers. The approximation follows from an assumption that the growth of a laminar boundary layer in the presence of pressure gradients could be parameterized solely as a function of the Holstein–Bohlen flow parameter, thus reducing the von Kármán momentum integral to a first-order ordinary differential equation. This method is useful for the analysis of laminar flows, and in computational potential flow solvers to account for the viscous effects. In this work, an approximate method for determining the momentum thickness of a two-dimensional, turbulent boundary layer is proposed following Thwaites’ work. It is shown that the method provides good estimates of the momentum thickness for multiple boundary layers, including both favourable and adverse pressure gradient effects, up to the point of separation. In the limit of high Reynolds numbers, it is possible to derive a criterion for the onset of separation from the proposed model, which is shown to be in agreement with prior empirical observations. The sensitivity of the separation location with respect to upstream perturbations is also analysed through this model for the NASA/Boeing speed bump and the transonic Bachalo–Johnson bump.

17 WIND ENERGY↗

The scoping, design, and plasma physics optimization of the Eos neutron source stellarator

On the path to a fusion pilot plant, Thea Energy plans to build Eos, a sub-breakeven, deuterium-deuterium, beam-target fusion, stellarator neutron source facility for producing tritium and other valuable radioisotopes. In this paper, a set of 1D plasma physics models are coupled and used to design the operating point of the facility and predict performance. At this foundational stage of the design, analytic and approximate models are sufficient to capture the leading-order effects, and fast enough to run in the inner loop of an optimizer. Higher-fidelity analyses will follow. Models of 1D profile-dependent neutral beam stopping, ion beam slowing down, beam-target fusion, electron-ion classical heat transfer, energy confinement (ISS04), beam pressure, beam heating of ions and electrons, beam-beam fusion fraction, and neutral beam injection and gyrotron heating electrical efficiencies are included. A numerical optimizer is used to determine the minimum required facility electric power to generate tritium at a given rate. A potentially advantageous regime is described in which modern precisely-quasisymmetric stellarators, new high-temperature superconductors, ITER-derived neutral beam injection, and new high-frequency gyrotrons enable a suitible target plasma with hot electrons, cold ions, peaked density and temperature profiles, and high beam-injected ion density. It appears possible at this time for a facility with a medium-scale and medium-strength stellarator whose required facility electric power is less than 40 MW to produce $2.5\times 10^{17}$ neutrons s -1 for the production of radioisotopes. With the addition of a tritium breeding blanket, such a facility could produce 0.2 grams d -1 or 70 grams yr -1 of tritium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A macro-micromechanics analysis of a notched metal matrix composite

A macro-micromechanics analysis was formulated to determine the matrix and fiber behavior near the notch tip in a center-notched metal matrix composite. Results are presented for a boron/aluminum monolayer. The macro-level analysis models the entire notched specimen using a three dimensional finite element program which uses the vanishing-fiber-diameter model to model the elastic-plastic behavior of the matrix and the elastic behavior of the fiber. The micro-behavior is analyzed using a Discrete Fiber-Matrix (DFM) model containing one fabric and the surrounding matrix. The dimensions of the DFM model were determined by the ply thickness and the fiber volume fraction and corresponded to the size of the notch-tip element in the macro-level analysis. The boundary conditions applied to the DFM model were determined from the macro-level analysis. Stress components within the DFM model were calculated and stress distributions are presented along selected planes and surfaces within the DFM model, including the fiber-matrix interface. Yielding in the matrix was examined at the notch tip in both the macro- and micro-level analyses. The DFM model predicted higher stresses (24 percent) in the fiber compared to the global analysis. In the notch-tip element, the interface stresses indicated that a multi-axial criterion may be required to predict interfacial failure. The DFM analysis predicted yielding to initiate in the notch-tip element at a stress level 28 percent lower than predicted by the global analysis.

Bigelow, C. A.↗

Finite-Difference Solutions of the Alternate Turbopump Development High-Pressure Oxidizer Turbopump Pump-End Ball-Bearing Cavity Flows

These analyses were undertaken to aid in the understanding of flow phenomena in the Alternate Turbopump Development (ATD) High-pressure Oxidizer Turbopump (HPOTP) Pump-end ball bearing (PEBB) cavities and their roles in turbopump vibration initiation and bearing distress. This effort was being performed to provide timely support to the program in a decision as to whether or not the program should be continued. In the first case, it was determined that a change in bearing through flow had no significant effect on axial preload. This was a follow-on to a previous study which had resulted in a redesign of the bearing exit cavity which virtually eliminated bearing axial loading. In the second case, a three-dimensional analysis of the inner-race-guided cage configuration was performed so as to determine the pressure distribution on the outer race when the shaft is 0.0002 inches off-center. The results indicate that there is virtually no circumferential pressure difference caused by the offset to contribute to bearing tilt. In the third case, axisymmetric analyses were performed on an outer-race guided cage configuration to determine the magnitude of tangential flow entering the bearing. The removed-shoulder case was analyzed as was the static diverter case. A third analysis where the preload spring was shielded by a sheet of metal for the baseline case was also performed. It was determined that the swirl entering the bearing was acceptable and the project decided to use the outer-race-guided cage configuration. In the fourth case, more bearing configurations were analyzed. These analyses included thermal modeling so as to determine the added benefit of injecting colder fluid directly onto the bearing inner-race contact area. The results of these analyses contributed to a programmatic decision to include coolant injection in the design.

Benjamin, Theodore G.↗

A STRONGLY COUPLED REACTOR CORE ISOLATION COOLING SYSTEM MODEL FOR EXTENDED STATION BLACK-OUT ANALYSES

The reactor core isolation cooling (RCIC) system in a boiling water reactor (BWR) provides makeup cooling water to the reactor pressure vessel (RPV) when the main steam lines are isolated and the normal supply of water to the reactor vessel is lost. The RCIC system operates independently of AC power, service air, or external cooling water systems. The only required external energy source is from the battery to maintain the logic circuits to control the opening and/or closure of valves in the RCIC systems in order to control the RPV water level by shutting down the RCIC pump to avoid overfilling the RPV and flooding the steam line to the RCIC turbine. It is generally considered in almost all the existing station black-out accidents (SBO) analyses that loss of the DC power would result in overfilling the steam line and allowing liquid water to flow into the RCIC turbine, where it is assumed that the turbine would then be disabled. This behavior, however, was not observed in the Fukushima Daiichi accidents, where the Unit 2 RCIC functioned without DC power for nearly three days. Therefore, more detailed mechanistic models for RCIC system components are needed to understand the extended SBO for BWRs. As part of the effort to develop the next generation reactor system safety analysis code RELAP-7, we have developed a strongly coupled RCIC system model, which consists of a turbine model, a pump model, a check valve model, a wet well model, and their coupling models. Unlike the traditional SBO simulations where mass flow rates are typically given in the input file through time dependent functions, the real mass flow rates through the turbine and the pump loops in our model are dynamically calculated according to conservation laws and turbine/pump operation curves. A simplified SBO demonstration RELAP-7 model with this RCIC model has been successfully developed. The demonstration model includes the major components for the primary system of a BWR, as well as the safety system components such as the safety relief valve (SRV), the RCIC system, the wet well, and the dry well. The results show reasonable system behaviors while exhibiting rich dynamics such as variable flow rates through RCIC turbine and pump during the SBO transient. The model has the potential to resolve the Fukushima RCIC mystery after adding the off-design two-phase turbine operation model and other additional improvements.

RCIC↗

Effects of Omitting Non-confounding Predictors From General Relative-Risk Models for Binary Outcomes

The effects, in terms of bias and precision, of omitting non-confounding predictive covariates from generalized linear models have been well studied, and it is known that such omission results in attenuation bias but increased precision with logistic regression. However, many epidemiologic risk analyses utilize alternative models that are not based on a linear predictor, and the effect of omitting non-confounding predictive covariates from such models has not been characterized. We employed simulation to study the effects on risk estimation of omitting non-confounding predictive covariates from an excess relative risk (ERR) model and a general additive-multiplicative relative-risk mixture model for binary outcome data in a case-control setting. We also compared the results to the effects with ordinary logistic regression. For these commonly employed alternative relative-risk models, the bias was similar to that with logistic regression when the risk was small. More generally, the bias and standard error of the risk-parameter estimates demonstrated patterns that are similar to those with logistic regression, but with greater magnitude depending on the true value of the risk. The magnitude of bias and standard error had little relation to study size or underlying disease prevalence. Prior conclusions regarding omitted covariates in logistic regression models can be qualitatively applied to the ERR and the general additive-multiplicative relative-risk mixture model without substantial change. Quantitatively, however, these alternative models may have slightly greater omitted-covariate bias, depending on the magnitude of the true risk being estimated.

60 APPLIED LIFE SCIENCES↗

Thermal Extraction and Delivery from Advanced Reactor Systems

This document presents preliminary thermodynamic analysis of thermal extraction and delivery from advanced nuclear reactors. Light-water, high temperature gas, sodium fast, and molten salt reactors are templated and analyzed for thermal extraction from nominal steam production conditions. Uncontrolled, controlled, and partial admission steam extraction methods are described and analyzed to show the differences in system impact from the various methods. Thermal delivery modeling results show comparisons that have become possible through the computational capabilities of the thermal delivery model. Due to the significant degrees of freedom of the problem, specific cases are used to compare the thermodynamic costs of heat transport using various heat transfer fluids. All this work is done using Modelica models, allowing for the future application of models in dynamic analyses beyond the steady-state calculations presented in this report. Appendices describe additional capabilities and methods that can help users analyze complex problems using HYBRID models developed in Modelica.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗