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 1,243 records · Page 69

Alternate Endpoints for Deep Vadose Zone Environments: Challenges, Opportunities, and Progress - 13036

Current requirements for site remediation and closure are standards-based and are often overly conservative, costly, and in some cases, technically impractical to achieve. Use of risk-informed alternate endpoints provide a means to achieve remediation goals that are permitted by regulations and are protective of human health and the environment. Alternate endpoints enable establishing a path for cleanup that may include intermediate remedial milestones and transition points and/or regulatory alternatives to standards-based remediation. A framework is presented that is centered around developing and refining conceptual models in conjunction with assessing risks and potential endpoints as part of a system-based assessment that integrates site data with scientific understanding of processes that control the distribution and transport of contaminants in the subsurface and pathways to receptors. This system-based assessment and subsequent implementation of the remediation strategy with appropriate monitoring are targeted at providing a holistic approach to addressing risks to human health and the environment. This holistic approach also enables effective predictive analysis of contaminant behavior to provide defensible criteria and data for making long-term decisions. Developing and implementing an alternate endpoint-based approach for remediation and waste site closure presents a number of challenges and opportunities. Categories of these challenges include scientific and technical, regulatory, institutional, and budget and resource allocation issues. Opportunities exist for developing and implementing systems-based approaches with respect to supportive characterization, monitoring, predictive modeling, and remediation approaches. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Kinetics of Cation and Oxyanion Adsorption and Desorption on Ferrihydrite: Roles of Ferrihydrite Binding Sites and a Unified Model

Understanding the kinetics of toxic ion reactions with ferrihydrite is crucial for predicting the dynamic behavior of contaminants in soil environments. In this study, the kinetics of As(V), Cr(VI), Cu, and Pb adsorption and desorption on ferrihydrite were investigated with a combination of laboratory macroscopic experiments, microscopic investigation and mechanistic modeling. The rates of As(V), Cr(VI), Cu, and Pb adsorption and desorption on ferrihydrite, as systematically studied using a stirred-flow method, was highly dependent on the reaction pH and metal concentrations and varied significantly among four metals. Spherical aberration-corrected scanning transmission electron microscopy (Cs-STEM) showed, at sub-nano scales, all four metals were distributed within the ferrihydrite particle aggregates homogeneously after adsorption reactions, with no evidence of surface diffusion-controlled processes. Based on experimental results, we developed a unifying kinetics model for both cation and oxyanion adsorption/desorption on ferrihydrite based on the mechanistic-based equilibrium model CD-MUSIC. Overall, the model described the kinetic results well, and we quantitatively demonstrated how the equilibrium properties of the cation and oxyanion binding to various ferrihydrite sites affected the adsorption and desorption rates. Our results provided a unifying quantitative modeling method for the kinetics of both cation and oxyanion adsorption/desorption on iron minerals.

Environmental Molecular Sciences Laboratory↗

Developing and Testing a Novel Stochastic Ice Microphysics Parameterization for Cloud and Climate Models Using ARM Field Campaign Data (Final Progress Report)

The major goals of this project were: 1) to use field campaign measurements from DOE’s Atmospheric Radiation Measurement (ARM) program to characterize variability of important parameters describing properties of ice particles in the atmosphere; 2) based on this observational analysis, to develop a parameterization scheme for weather and climate models that stochastically varies these parameters, and implement the new scheme into a weather model called the Weather Research and Forecasting model (WRF); 3) to use WRF coupled with the new stochastic scheme to simulate ARM field campaign thunderstorm cases and analyze how accounting for this parameter variability affects the model simulations. This work was performed jointly between the National Center for Atmospheric Research, University of Oklahoma, and University of Utah. To accomplish these goals, we extended an approach previously developed to characterize the variability in the size distribution of ice particles to parameters that are explicitly represented in models (i.e., relationships between ice particle mass and size, and between particle fall velocity and size). Our project was, to our knowledge, the first to apply observationally-constrained estimates of this parameter variability describing mass-size and fall velocity-size in a modeling framework. Our results showed efficacy of the approach, evaluated using ARM observations. Similarly, to our knowledge, work in this project was the first to propose and evaluate in detail a stochastic approach for unresolved turbulent mixing in high-resolution model simulations against detailed, benchmark large eddy simulations and ARM observations. Results showed some promising behavior, particularly with increased mixing and dilution of air in thunderstorm cores with surrounding environmental air, bringing the stochastic simulations closer to the benchmark large eddy simulations; however, results were somewhat degraded using stochastic mixing compared to observations from the AMIE/DYNAMO field campaign. This project also further refined and applied a modeling methodology called “piggybacking” that can robustly separate dynamical and thermodynamic impacts of model changes, and comparison studies of different models based on cases developed from ARM observations. Finally, this project directly supported three graduate students who completed their PhDs as well as a postdoctoral research fellow.

54 ENVIRONMENTAL SCIENCES↗

Dynamics of simulated Atlantic upwelling annual cycle in CMIP5 models

As one of the prominent climatic component of tropical climate system, the annual cycle of the equatorial Atlantic cold tongue region has been a subject of significant research. In this work, we examine the fidelity of the key part of the annual cycles, namely, the upwelling annual cycles, from the simulations of the Coupled Model Intercomparison Project Phase 5 (CMIP5). Utilizing the simple framework we have recently developed based on the Zebiak–Cane (ZC) ocean model, we first tested the validity of the simple framework of upwelling as a combination of the local wind–driven Ekman upwelling and nonlocal wind–driven wave upwelling. It is demonstrated that the theoretically reconstructed upwelling annual cycles to large extent are in good agreement with the simulations in these coupled climatic models. Comparing with observations, we show that the semiannual component of the upwelling is systematically too strong. This significant bias is attributable to the dynamic amplification to the wind–driven wave upwelling originated from a less significant bias in the semiannual harmonic component in the central equatorial wind field. The latter is dynamically amplified in the wind–driven wave upwelling owing to its relatively fast time scale. Furthermore, the realistic simulations of equatorial annual cycles may sensitively depend on good simulations of the semiannual harmonic components in equatorial wind fields.

54 ENVIRONMENTAL SCIENCES↗

A multiphase flow model of water droplets dielectrophoretic-induced air dehumidification phenomena

Air humidity in indoor spaces plays a critical role in human comfort and health. Dehumidification systems are used for building humidity controls, but they can take significant energy consumption, especially in geographic locations with high outdoor humidity and warm climates. Consequently, there is a growing demand for innovative dehumidification processes that consume minimal energy. Dielectrophoretic air dehumidification represents one such promising approach. However, it has not garnered significant attention due to the absence of engineering models and simulation tools capable of evaluating its performance and limitations at large-scale airflows. A new numerical multiphase CFD model, which is also experimentally validated, is developed in a customized Reacting Foam solver based on OpenFOAM® version 9. The newly developed model seeks to decrease substantial energy consumption and lower costs by leveraging the dielectrophoretic phenomenon to regulate moisture levels in the air. The solver integrates a hybrid Eulerian-Lagrangian framework to track the droplet's trajectory and growth rate while solving the continuum equations for the moist air. An electrospray produces electrically charged droplets, which grow during their in-flight trajectories as water vapor condenses onto their surfaces. The role of electrostatic forces in promoting vapor condensation within a high-gradient electrical field is investigated, and the dielectrophoretic vapor nucleation process on charged water droplets is discussed. The CFD model was validated against results from the literature and from proof-of-concept experiments conducted by the authors, which showed a 2 % air dehumidification with a single electrospray and airflow rate of 5 cubic feet per minute. The simulation results indicated that augmenting the number of electrically charged spray droplets increased the dehumidification of the air to 25 %. The initial mean droplet diameter, the orientation of the injector and relative humidity significantly influence the assessment of dehumidification. As a result, scaling up this approach to larger airflow volumes is identified as a potential future research direction.

42 ENGINEERING↗

A voxel-based multiscale model to simulate the radiation response of hypoxic tumors

Purpose: In radiotherapy, it is important to predict the response of tumors to irradiation prior to the treatment. This is especially important for hypoxic tumors, which are known to be highly radioresistant. Mathematical modeling based on the dose distribution, biological parameters, and medical images may help to improve this prediction and to optimize the treatment plan. Methods: A voxel-based multiscale tumor response model for simulating the radiation response of hypoxic tumors was developed. It considers viable and dead tumor cells, capillary and normal cells, as well as the most relevant biological processes such as (i) proliferation of tumor cells, (ii) hypoxia-induced angiogenesis, (iii) spatial exchange of cells leading to tumor growth, (iv) oxygen-dependent cell survival after irradiation, (v) resorption of dead cells, and (vi) spatial exchange of cells leading to tumor shrinkage. Oxygenation is described on a microscopic scale using a previously published tumor oxygenation model, which calculates the oxygen distribution for each voxel using the vascular fraction as the most important input parameter. To demonstrate the capabilities of the model, the dependence of the oxygen distribution on tumor growth and radiation-induced shrinkage is investigated. In addition, the impact of three different reoxygenation processes is compared and tumor control probability (TCP) curves for a squamous cells carcinoma of the head and neck (HNSSC) are simulated under normoxic and hypoxic conditions. Results: The model describes the spatiotemporal behavior of the tumor on three different scales: (i) on the macroscopic scale, it describes tumor growth and shrinkage during radiation treatment, (ii) on a mesoscopic scale, it provides the cell density and vascular fraction for each voxel, and (iii) on the microscopic scale, the oxygen distribution may be obtained in terms of oxygen histograms. With increasing tumor size, the simulated tumors develop a hypoxic core. Within the model, tumor shrinkage was found to be significantly more important for reoxygenation than angiogenesis or decreased oxygen consumption due to an increased fraction of dead cells. In the studied HNSSC-case, the TCD{sub 50} values (dose at 50% TCP) decreased from 71.0 Gy under hypoxic to 53.6 Gy under the oxic condition. Conclusions: The results obtained with the developed multiscale model are in accordance with expectations based on radiobiological principles and clinical experience. As the model is voxel-based, radiological imaging methods may help to provide the required 3D-characterization of the tumor prior to irradiation. For clinical application, the model has to be further validated with experimental and clinical data. If this is achieved, the model may be used to optimize fractionation schedules and dose distributions for the treatment of hypoxic tumors.

60 APPLIED LIFE SCIENCES↗

Initial microstructural model for creep-fatigue damage in Grade 91 steel

This report describes an initial Crystal Plasticity Finite Element Method (CPFEM) model for cyclic plasticity and damage in Grade 91 steel. The objective of this work is to develop a framework for modeling creep-fatigue interaction in Grade 91 steel to better predict the onset of damage in high temperature microreactor components. Many microreactor concepts envision low operating pressures but relatively high thermal stresses. Under these conditions, creep-fatigue will likely be the dominant design failure mechanism. Physically based models, like the one under development here, could lead to a better understanding of creep-fatigue mechanisms and the effect of stress multiaxiality and hold time on creep-fatigue damage. In turn, this could lead to more efficient microreactor component designs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

ARM Data-Oriented Metrics and Diagnostics Package for Climate Model Evaluation

A Python-based metrics and diagnostics package is currently being developed by the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Infrastructure Team at Lawrence Livermore National Laboratory (LLNL) to facilitate the use of long-term, high-frequency measurements from the ARM Facility in evaluating the regional climate simulation of clouds, radiation, and precipitation. This metrics and diagnostics package computes climatological means of targeted climate model simulation and generates tables and plots for comparing the model simulation with ARM observational data. The Coupled Model Intercomparison Project (CMIP) model data sets are also included in the package to enable model intercomparison as demonstrated in Zhang et al. (2017). The mean of the CMIP model can serve as a reference for individual models. Basic performance metrics are computed to measure the accuracy of mean state and variability of climate models. The evaluated physical quantities include cloud fraction, temperature, relative humidity, cloud liquid water path, total column water vapor, precipitation, sensible and latent heat fluxes, and radiative fluxes, with plan to extend to more fields, such as aerosol and microphysics properties. Process-oriented diagnostics focusing on individual cloud- and precipitation-related phenomena are also being developed for the evaluation and development of specific model physical parameterizations. The version 1.0 package is designed based on data collected at ARM’s Southern Great Plains (SGP) Research Facility, with the plan to extend to other ARM sites. The metrics and diagnostics package is currently built upon standard Python libraries and additional Python packages developed by DOE (such as CDMS and CDAT). The ARM metrics and diagnostic package is available publicly with the hope that it can serve as an easy entry point for climate modelers to compare their models with ARM data. In this report, we first present the input data, which constitutes the core content of the metrics and diagnostics package in section 2, and a user's guide documenting the workflow/structure of the version 1.0 codes, and including step-by-step instruction for running the package in section 3.

54 ENVIRONMENTAL SCIENCES↗

Energy Performance Assessment of Radiant Cooling System through Modeling and Calibration at Component Level

The paper describes a case study of an information technology office building with a radiant cooling system and a conventional variable air volume (VAV) system installed side by side so that performancecan be compared. First, a 3D model of the building involving architecture, occupancy, and HVAC operation was developed in EnergyPlus, a simulation tool. Second, a different calibration methodology was applied to develop the base case for assessing the energy saving potential. This paper details the calibration of the whole building energy model to the component level, including lighting, equipment, and HVAC components such as chillers, pumps, cooling towers, fans, etc. Also a new methodology for the systematic selection of influence parameter has been developed for the calibration of a simulated model which requires large time for the execution. The error at the whole building level [measured in mean bias error (MBE)] is 0.2%, and the coefficient of variation of root mean square error (CvRMSE) is 3.2%. The total errors in HVAC at the hourly are MBE = 8.7% and CvRMSE = 23.9%, which meet the criteria of ASHRAE 14 (2002) for hourly calibration. Different suggestions have been pointed out to generalize the energy saving of radiant cooling system through the existing building system. So a base case model was developed by using the calibrated model for quantifying the energy saving potential of the radiant cooling system. It was found that a base case radiant cooling system integrated with DOAS can save 28% energy compared with the conventional VAV system.

Building Energy Simulation↗

Network-Wide Traffic Signal Control Using Bilinear System Modeling and Adaptive Optimization

This study proposes a new multi-input multi-output optimal bilinear signal control method in which a bilinear dynamic model approximation is used to capture the nonlinear dynamics of the urban traffic networks. With signal green time splits as the control input and traffic delay changes as the output for each intersections in the network, a bilinear system model was developed, which, on the basis of linear system modeling, takes interactions among traffic delays and signal timing splits into consideration. Based on the bilinear system modeling framework, we conducted two steps in each time interval to derive traffic control strategies: (1) we used the normalized least-squared algorithm to estimate system parameters; and (2) we solved an online optimization problem to obtain the updated traffic control inputs for the signal timing that minimizes future traffic delays. We evaluated the proposed method in a microscopic traffic simulation environment (VISSIM) with a 35-intersection network of Bellevue city in Washington. Two different traffic demand patterns: (1) normal traffic demands; and (2) time-varying traffic demands were simulated to compare the performance of different control strategies. Experimental results show that (1) the proposed bilinear system model can better describe traffic system dynamics than linear-model based methods, such as our previously developed linear-quadratic regulator control; and (2) the proposed method outperforms the state-of-the-art signal control strategies, namely the max-pressure and the self-organizing traffic light control methods. We have also shown that the proposed method is applicable to all other possible network layouts and signal controller phasing structures.

42 ENGINEERING↗

Vibrational Energy Harvesting Using a Cantilever Model

Vibrational energy harvesting (VEH) is a method of capturing incidental mechanical vibrational energy and converting it to electrical energy. This is enabled by two technologies: Electromagnetic induction via a cantilever or piezoelectric devices. When designing a VEH system, a fast forward model is desired for response determination and optimal parameter estimation. An ordinary different equation (ODE) system model is developed for the cantilever system based upon the derivations of [1] and [2], and compared with a full electromechanical COMSOL model.

42 ENGINEERING↗

Predicting cell-to-cell failure propagation and limits of propagation in lithium-ion cell stacks

Thermal runaway of lithium-ion batteries is a risk that is magnified when stacks of lithium-ion cells are used for large scale energy storage. When limits of propagation can be identified so that systems can be designed to prevent large scale cascading failure even if a failure does occur, these systems will be safer. This work addresses the prediction of cell-to-cell failure propagation and the propagation limits in lithium-ion cell stacks to better understand and identify safe designs. A thermal-runaway model is presented based on recent developments in thermochemical source terms. It is noted that propagating failure is characterized by temperatures above which calorimetry data is available. Results show high temperature propagating failure predictions are too rapid unless an intra-particle diffusion limit is included, introducing a Damköhler number limiter into the rate expression. This new model form is evaluated against cell-to-cell failure propagation where the end cell of a stack is forced into thermal runaway through a nail-induced short circuit. Limits of propagation for this configuration are identified. Results show cell-to-cell propagation predictions are consistent with measurements over a range of cell states of charge and with the introduction of metal plates between cells to add system heat capacity representative of structural members. This consistency extends from scenarios where propagation occurs through scenarios where propagation is prevented.

25 ENERGY STORAGE↗

Deep learning for drag force modelling in dilute, poly-dispersed particle-laden flows with irregular-shaped particles

Here, this study applies machine learning-based approaches to develop a drag force model for irregular-shaped particles in incompressible flows. The particle-laden flows are studied through an in-house particle-resolved direct numerical simulation (PR-DNS) at low-intermediate Reynolds numbers (Re). We utilize the PR-DNS to obtain drag force coefficients and flow fields of single particles. A variational auto-encoder model is used to obtain latent vectors to represent the geometrical features of the particles, and artificial neural networks (ANN) are developed to predict drag force coefficients and flow fields of the single particles. This study applies a pairwise interaction extended point-particle (PIEP) model to obtain the coefficients assuming the flow fields of neighboring particles can be linearly superposed. The PIEP method shows significant improvement on prediction for a few neighbored particles. In addition, the results reveal R2 scores of 0.56-0.62 and errors of 9.1-10.0 % for the dilute, polydispersed systems with a volume fraction of 0.5 %.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Assimilation of Wave Imaging Radar Observations for Real-time Wave-by-Wave Forecasting

This project addressed Topic 3: “Wave Measurement Instrumentation for Feed Forward Controls” under the FOA number DE-FOA-0000971. The overall goal of the program was to develop a phase-resolving wave forecasting technique for application to the active control of Wave Energy Conversion (WEC) devices. We have developed an approach that couples a wave imaging marine radar with a phase-resolving linear wave model for real-time wave field reconstruction and forward propagation of the wave field in space and time. The scope of the project was to develop and assess the performance of this novel forecasting system. Specific project goals were as follows: Develop and verify a fast, GPU-based (Graphical Processing Unit) wave propagation model suitable for phase-resolved computation of nearshore wave transformation over variable bathymetry; Compare the accuracy and speed of performance of the wave model against a deep water model in their ability to predict wave field transformation in the intermediate water depths (50 to 70 m) typical of planned WEC sites; Develop and implement a variational assimilation algorithm that can ingest wave imaging radar observations and estimate the time-varying wave conditions offshore of the domain of interest such that the observed wave field is best reconstructed throughout the domain and then use this to produce model forecasts for a given WEC location; Collect wave-resolving marine radar data, along with relevant in situ wave data, at a suitable wave energy test site, apply the algorithm to the field data, assess performance, and identify any necessary improvements; and Develop a production cost estimate that addresses the affordability of the wave forecasting technology and include in the Final Report. The developed forecasting algorithm (“Wavecast”) was evaluated for both speed and accuracy against a substantial synthetic dataset. Early in the project, performance tests definitively demonstrated that the system was capable of forecasting in real-time, as the GPU-based wave model backbone was very computationally efficient. The data assimilation algorithm was developed on a polar grid domain in order to match the sampling characteristics of the observation system (wave imaging marine radar). For verification purposes, a substantial set of synthetic wave data (i.e. forward runs of the wave model) were generated to be used as ground truth for comparison to the reconstructions and forecasts produced by Wavecast. For these synthetic cases, Wavecast demonstrated very good accuracy, for example, typical forecast correlation coefficients were between 0.84-0.95 when compared to the input data. Dependencies on shadowing, observational noise, and forecast horizon were also identified. During the second year of the project, a short field deployment was conducted in order to assess forecast accuracy under field conditions. For this, a radar was installed on a fishing vessel and observations were collected at the South Energy Test Site (SETS) off the coast of Newport, OR. At the SETS site, simultaneous in situ wave observations were also available owing to an ongoing field project funded separately. Unfortunately, the position and heading information that was available for the fishing vessel were not of sufficient accuracy in order to validate the forecast in a phase-resolving sense. Instead, a spectral comparison was made between the Wavecast forecast and the data from the in situ wave buoy. Although the wave and wind conditions during the field test were complex, the comparison showed a promising reconstruction of the wave spectral shape, where both peaks in the bimodal spectrum were represented. However, the total reconstructed spectral energy (across all directions and frequencies) was limited to 44% of the observed spectrum. Overall, wave-by-wave forecasting using a data assimilation approach based on wave imaging radar observations and a physics-based wave model shows promise for short-term phase-resolved predictions. Two recommendations for future work are as follows: first, we would recommend additional focused field campaigns for algorithm validation. The field campaign should be long enough to capture a range of wave conditions relevant to the target application and WEC site. In addition, it will be crucial to make sure the vessel of choice has high accuracy position and heading instrumentation (this instrumentation is commercially available but not standard on commercial fishing vessels). The second recommendation is to expand the model physics in the wave model backbone to include some nonlinear effects. Specifically, the third-order correction to the wave speed due to amplitude dispersion would be the next step in order to more accurately represent the phase speeds of large amplitude waves.

16 TIDAL AND WAVE POWER↗

Plasma shape and position control development for NSTX-U using the GSEvolve plasma simulator

NSTX-U discharges have been reproduced using the GSEvolve plasma simulator in order to improve their control and reduce oscillations during the plasma current (Ip) ramp-up phase, and with the ultimate goal of enabling reliable H-mode access with a double-null shape at high elongation, low internal inductance, and high plasma performance. The control objective is to eliminate undesired oscillations in the plasma shape and vertical position that start when the plasma diverts during the Ip ramp-up phase, and that happened recurrently during the first NSTX-U experimental campaign. Model-based control solutions have been developed, including alternative H-mode access recipes and a Multi-Input Multi-Output (MIMO) decoupling controller. These solutions have demonstrated their capability to reduce or eliminate the aforementioned oscillations in GSEvolve simulations, enabling reproducible H-mode access. The results of this work aim to inform NSTX-U experimental operations on potential control developments to maximize its scientific output, and highlight the powerful capabilities of high-fidelity control-simulation frameworks like GSEvolve.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗