Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mathematical model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 577 records · Page 32

Quantification of Load Flexibility in Residential Buildings Using Home Energy Management Systems

With increasing penetration of renewable energy resources, the flexibility of operating behind-the-meter (BTM) resources plays a key role in enhancing grid reliability and resilience. Residential buildings with home energy management systems (HEMS) can provide desired flexibility for the distribution system operator (DSO) while considering customer comfort and preferences. This paper discusses a methodology to quantify the flexibility of BTM resources of residential buildings using HEMS. First, we propose a model predictive control framework to formulate the flexibility band comprising nominal, upper, and lower demand profiles. Second, the paper proposes a dispatch method for HEMS to compute the control signals for each BTM resource (e.g., air conditioner, water heater, home battery system) upon receiving a flexibility service request from the DSO. The case study provides insight into the flexibility provided at the whole-home level with different user preferences and seasons. The results demonstrate that HEMS is capable of providing flexibility service at the request of the DSO while delivering primary services to the building occupants.

41 EE - Solar Energy Technologies Office (EE-4S)↗

Inelastic deformation of diamond single crystals shock compressed to multimegabar stresses: Wave profile calculations

As the archetypal strong solid, the response of diamond shock compressed to multimegabar stresses is important for fundamental science and for numerical simulations of wave profiles for applications in high energy density physics experiments. Previous experiments and analysis have shown that the commonly used hydrodynamic assumption is invalid for diamond shock compressed to stresses below melt and an elastic–inelastic description is needed. Here, we present a phenomenological material model for calculating wave profiles in shock compressed diamond single crystals that incorporates this description. Also, to support the modeling effort, we carried out wave profile measurements on shock compressed diamond single crystals at the Sandia Z facility to augment previous measurements. Wave profiles for [100] and [111] diamond calculated using the material model provide a good match to the elastic–inelastic response (observed two-wave structure) measured at ∼325 and ∼360 GPa. Furthermore, the calculated peak stresses for single (overdriven) waves provide a good match to the measured Hugoniot states for stresses reaching ∼700 GPa, which is near melting conditions. The present results show that the diamond single crystal response at multimegabar shock stresses is characteristic of a brittle solid—pressure-dependent strength and strength loss due to inelastic deformation.

Deformation↗

DC Arc Incident Energy in Photovoltaic Systems: Methods for Evaluation

Renewable energy systems continue to be one of the fastest growing segments of the energy industry. This article focuses on the understanding of how photovoltaic (PV) technology behaves under dc arc conditions. Emphasis is placed on the electrical safety aspect of dc arc-flash incident energy (IE) evaluation. Because of the fast proliferation of PV systems and lack of formal equivalent calculation guidelines, such as IEEE 1584 for ac systems, it has been necessary to rely on different equations and models presented by various researchers over the last few years. This article discusses the behavior of PV systems under arc conditions and presents the results of the available methods to estimate the dc arc-flash IE. It provides a comparative analysis of a proposed arc-flash IE calculation method against different laboratory tests, including those performed for this article at the National Renewable Energy Laboratory (NREL). Detailed explanations are provided regarding the effect of the PV module current-voltage (I-V) and power-voltage (P-V) curves under arcing conditions. Examples of the application of the proposed calculation method to the test measurements are included.

arc discharges↗

Co-design optimization of combined heat and power-based microgrids

With the emergent need for clean and reliable energy resources, hybrid energy systems, such as the microgrid, are widely adopted in the United States. A microgrid can consist of various distributed energy resources, for instance, combined heat and power (CHP) systems. Here, the CHP module is a distributed cogeneration technology that produces electricity and recaptures heat generated as a by-product. It is an energy-efficient technology converting heat that would otherwise be wasted to valuable thermal energy. For an optimal system configuration, this study develops a novel co-design optimization framework for CHP-based cogeneration microgrids. The framework provides the stakeholder with a method to optimize investments and attain resilient operations. The proposed co-design framework has a mixed integer programming (MIP) model that outputs decisions for both plant designs and operating controls. The microgrid considered in this study contains six components: the CHP, boiler, heat recovery unit, thermal storage system, power storage system, and photovoltaic plant. After solving the MIP model, the optimal design parameters of each component can be found to minimize the total installation cost of all components in the microgrid. Furthermore, the online costs from energy production, operation, maintenance, machine startup, and disruption-induced unsatisfied loads are minimized by solving the optimal control decisions for operations. Case studies based on designing a CHP-based microgrid with empirical data are conducted. Moreover, we consider both nominal and disruptive operational scenarios to validate the performance of the proposed co-design framework in terms of a cost-effective, resilient system.

42 ENGINEERING↗

System Modeling and Robust Stability Region Analysis for Multi-inverters Based on VSG

As the number of inverter-interfaced distributed generators (IIDGs) increases, the decrease in the inertia of the microgrid affects the stability of microgrid voltage and current. The virtual synchronous generator (VSG) control method simulates the output characteristics of the synchronous machine to improve the inertia and damping of the microgrid. However, low-frequency oscillations can occur in systems with multi-VSGs connected in parallel, primarily when operating in island mode. Therefore, this article develops the small-signal state-space model for multi-inverters based on VSG and analyses the robust stability region (RSR) of the control parameters. Firstly, an extensible small-signal state-space model is established and simplified in the dq axis, which is a necessary preprocessing for the subsequent stability analysis. Secondly, the guardian map method based on bialternate product operation is constructed, and the RSR is established through translation mapping and rotation mapping. Therefore, the solution of RSR is transformed into the problem of identifying the Hurwitz matrix. Finally, the simulations and experiments verify that the calculated region can ensure stable operation during startup and meet the requirement to suppress low-frequency oscillations. In conclusion, the method of solving RSR proposed in this paper provides a basis for selecting controller parameters in practical engineering.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Genetic Algorithm-Based Commutation Angle Control for Torque Ripple Mitigation in Switched Reluctance Motor Drives

This work addresses the application of the Genetic Algorithm (GA) technique to optimize the commutation angles of a 2 kW 8/6 switched reluctance machine (SRM). The primary goal is to reduce the well-known drawback of SRMs: the torque ripple. Firstly, the machine was modeled in Matlab /Simulink ® using lookup tables obtained via finite element method (FEM) simulations. Subsequently, the model was used to perform the GA routine aiming to find the optimal phase commutation angles that minimize the torque ripple factor. Notably, the torque performance of the SRM was significantly affected by the commutation angles during the search for the optimal solution. Afterwards, the GA results for four different operation points were verified experimentally through a developed drive platform with digital signal processor-based (DSP) control and an asymmetric bridge converter. As showed by the experiments, the proposed approach was suitable to reduce the torque ripple by more than 50% for one of the evaluated operating points. Furthermore, it was confirmed that the torque ripple mitigation led to acoustic noise improvement.

42 ENGINEERING↗

Differential contagiousness of respiratory disease across the United States

The initial contagiousness of a communicable disease within a given population is quantified by the basic reproduction number, $R_0$. This number depends on both pathogen and population properties. On the basis of compartmental models that reproduce Coronavirus Disease 2019 (COVID-19) surveillance data, we used Bayesian inference and the next-generation matrix approach to estimate region-specific $R_0$ values for 280 of 384 metropolitan statistical areas (MSAs) in the United States (US), which account for 95% of the US population living in urban areas and 82% of the total population. We focused on MSA populations after finding that these populations were more uniformly impacted by COVID-19 than state populations. Our maximum a posteriori (MAP) estimates for $R_0$ range from 1.9 to 7.7 and quantify the relative susceptibilities of regional populations to spread of respiratory diseases.

60 APPLIED LIFE SCIENCES↗

Bayesian Inference of State-Level COVID-19 Basic Reproduction Numbers across the United States

Although many persons in the United States have acquired immunity to COVID-19, either through vaccination or infection with SARS-CoV-2, COVID-19 will pose an ongoing threat to non-immune persons so long as disease transmission continues. We can estimate when sustained disease transmission will end in a population by calculating the population-specific basic reproduction number ℛ0, the expected number of secondary cases generated by an infected person in the absence of any interventions. The value of ℛ0 relates to a herd immunity threshold (HIT), which is given by 1−1/ℛ0. When the immune fraction of a population exceeds this threshold, sustained disease transmission becomes exponentially unlikely (barring mutations allowing SARS-CoV-2 to escape immunity). Here, we report state-level ℛ0 estimates obtained using Bayesian inference. Maximum a posteriori estimates range from 7.1 for New Jersey to 2.3 for Wyoming, indicating that disease transmission varies considerably across states and that reaching herd immunity will be more difficult in some states than others. ℛ0 estimates were obtained from compartmental models via the next-generation matrix approach after each model was parameterized using regional daily confirmed case reports of COVID-19 from 21 January 2020 to 21 June 2020. Our ℛ0 estimates characterize the infectiousness of ancestral strains, but they can be used to determine HITs for a distinct, currently dominant circulating strain, such as SARS-CoV-2 variant Delta (lineage B.1.617.2), if the relative infectiousness of the strain can be ascertained. On the basis of Delta-adjusted HITs, vaccination data, and seroprevalence survey data, we found that no state had achieved herd immunity as of 20 September 2021.

60 APPLIED LIFE SCIENCES↗

Analysis of heat and mass transfer potential of a dew-point cooling tower in different climatic conditions

In this study, the performance of the Dew-Point Cooling Tower (DPCT) was analyzed for different factors in a variety of climate conditions. For this purpose, a dedicated numerical model describing heat and mass transfer processes was developed and validated. The results of the numerical simulations allowed to analyze the potential of utilizing the heat and mass transfer process with the dew-point phenomenon for water cooling. It was established that the operational parameters that have a high impact on the performance of the DPCT are: inlet water temperature and inlet air humidity ratio. It was also established that DPCT achieves the highest COP and Specific Cooling Capacity for cold subtropical highland climates and that it achieves highest Wet-bulb Effectiveness for monsoon-influenced humid subtropical climates. In conclusion, the regions where all the efficiency factors achieved above-average values included warm, arid, and desert climates.

42 ENGINEERING↗

Performance study of a novel dew point evaporative cooler in the climate of central Europe using building simulation tools

This paper presents the performance research of a hybrid air conditioning system (System II) equipped with novel dew point evaporative cooler (DPIEC) for a retail building located in temperate climatic conditions. The new cooler has a dedicated, adaptable structure that maximizes its performance in the climate conditions of central Europe. The unit can operate as a typical regenerative air cooler when outdoor air is drier than indoor air, it can operate as a counter-flow heat recovery unit when indoor conditions are more dry than outdoor and it can operate as a combination of both arrangements. The year-round, hourly-stepped building energy simulations were carried-out using a prognostic tool of energy consumption. The novel dew point evaporative cooler application potential was established using the black-box model based on regression equations was established based on experimental tests on a prototype of novel DPIEC. The hybrid system operation was compared to the typical air handling unit equipped with a standard energy wheel for heat recovery (System I). Presented research shows that DPIEC can cover about 95% of total cooling loads while the energy wheel covers only 9% of total cooling load. System II allows us to save 65% of seasonal electricity consumption comparing to System I. It comes from the fact that during most of the cooling season DPIEC is in operation and it covers most of the cooling loads without the need to switch the chilled water system.

Air conditioning↗

Modeling and simulation of gas separations with spiral–wound membranes

Models for gas separations with spiral–wound membranes are developed and found to exhibit good agreement with experiments performed on N 2 /O 2 mixtures. The two–dimensional (2D) model can be accurately approximated by a one–dimensional (1D) surrogate model when the spacer widths are chosen to make the channel pressure drops small. Subsequently, the separation of propane/propylene mixtures from the recycle purge stream of a polypropylene reactor is investigated. Assuming ideal gas is found to lead to significant overestimations in membrane stage cuts (sometimes more than 10%), an extent comparable to that associated with extrapolating constant olefin permeance from a low–pressure condition. While olefin permeance can change significantly with pressure, using a constant–permeance formulation can result in a small (< 2.5%) underprediction in stage cut if the value for the permeance is taken from the feed condition. Finally, membrane properties and costs necessary for a viable separation process are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mitigating the Impacts of Uncertain Geomagnetic Disturbances on Electric Grids: A Distributionally Robust Optimization Approach

Severe geomagnetic disturbances (GMDs) increase the magnitude of the electric field on the Earth’s surface (E-field) and drive geomagnetically-induced currents (GICs) along the transmission lines in electric grids. These additional currents can pose severe risks, such as current distortions, transformer saturation and increased reactive power losses, each of which can lead to system unreliability. Today several mitigation actions (e.g., changing grid topology) exist that can reduce the harmful GIC effects on the grids. Making such decisions can be challenging, however, because the magnitude and direction of the E-field are uncertain and non-stationary. In this paper, we model uncertain E-fields using the distributionally robust optimization (DRO) approach that determines optimal transmission grid operations such that the worst-case expectation of the system cost is minimized. We also capture the effect of GICs on the nonlinear AC power flow equations. For solution approaches, we develop an accelerated column-and-constraint generation (CCG) algorithm by exploiting a special structure of the support set of uncertain parameters representing the E-field. Extensive numerical experiments based on “epri-21” and “uiuc-150” systems, designed for GMD studies, demonstrate (i) the computational performance of the accelerated CCG algorithm, (ii) the superior performance of distributionally robust grid operations that satisfy nonlinear, nonconvex AC power flow equations and GIC constraints, in comparison with standard stochastic programming-based methods during the out-of-sample testing.

42 ENGINEERING↗

Accelerating Hanford Site Cleanup through Operations Research Modeling - 20238

The Hanford Site cleanup effort will require the integration of dozens of unique facilities and processes, many of which will be first-of-a-kind in implementation and design. Each facility will be governed by its own set of operating logic, configured with a unique array of unit operations, and subject to a set of constraints that will affect its behavior. The collection of facilities have multiple points of interface, making the operations of any one facility potentially significant to the operations of other up- or downstream processes. It is therefore highly desirable to accurately predict these operations, as it allows for Site officials to identify and preempt bottlenecks and vulnerabilities before they unexpectedly inhibit the cleanup mission. With the quantity and complexity of the processes that will be on Site, building a pen-and-paper or even a spreadsheet-assisted model of the cleanup mission quickly becomes overwhelming in scope and inaccurate in execution. The Engineering organization for the Site's Tank Operations Contract (TOC) has therefore implemented the use of operations research (OR) modeling to simulate and predict future operations of Site facilities. These models are created using a discrete event simulation tool that allows for the development of detailed, versatile, and robust models. Not only can these models account for complex logical behaviors, but they can also simulate process details down to the level of vessel sizing, labor utilization, equipment reliability, and resource availability. To date, the TOC has developed OR models for several facilities on Site, including for single-shell tank (SST) farms, double-shell tank (DST) farms, the Effluent Treatment Facility (ETF), and the waste transfer system. These models have focused on identifying bottlenecks and operational constraints, and have been used to quantify the effects of implementing process changes. This latter point is particularly valuable, as it allows for several alternatives to be studied in a virtual setting before committing resources to making a change in the field. The decision to develop OR models has gained tremendous support from the Site's stakeholders and the U.S. Department of Energy (DOE) management, and has prompted the use of the tool to support additional internal and external initiatives. Recently, an initiative was proposed to use the models to help identify and provide quantitative backing for risks and opportunities for the TOC. This application of OR could not only help inform how the TOC manages its risks (e.g. quantities and types of spare parts), but could also help drive process improvements whose benefits might otherwise be hard to quantify. The models have also been used to drive the TOC's cloud computing, artificial intelligence (AI), and machine learning (ML) initiatives. These initiatives will not only improve the ability of the TOC to more rapidly respond to the needs of its customers, but it will also aid in the ability of the TOC to analyze and improve the processes it studies. Partnership with two external software development and consulting companies (Lanner and Ynformed) has furthered not only the application of AI and ML within the TOC, but has also spurred the development of new/improved software tools and platforms used by the companies. These partnerships have proven to be mutually beneficial and productive, and have set a precedent for the types of gains that can be made by exploring such options. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A computational investigation of high-flux, plate-and-frame membrane modules for industrial carbon capture

In this work, we study the application of membrane-based separation systems for carbon capture, considering plate-and-frame membrane modules. The successful deployment of membrane CO 2 capture system relies on high-performing membranes as well as effective membrane modules that can fully exploit the developed membranes. A plate-and-frame membrane module is especially attractive for CO 2 capture from industrial flue gas due to its lower pressure drop compared to its counterparts such as spiral wound modules and hollow fiber modules. To design better plate-and-frame modules, we investigate their basic unit - a single membrane stack through a combination of computational modeling and experimental investigations. The modeling approach is based on Computational Fluid Dynamics (CFD) to represent a multiphysics problem, including the fluid flow and diffusion processes within a membrane module. We use experimental data collected under different operating conditions to validate the CFD model. Numerical results suggest a good agreement between experiments and model outputs for the CO 2 recovery, CO 2 mole fraction in the retentate and permeate, and stage-cut. The CFD model is able to predict accurately the flow behavior, providing valuable insights on the effects of fluid dynamics on mass transfer of CO 2 . We also carry out a sensitivity analysis to identify the effect of key parameters on the CO 2 recovery and the CO 2 purity of the outlet streams.

CFD simulation↗

Impacts of Vaccination and Severe Acute Respiratory Syndrome Coronavirus 2 Variants Alpha and Delta on Coronavirus Disease 2019 Transmission Dynamics in Four Metropolitan Areas of the United States

To characterize Coronavirus Disease 2019 (COVID-19) transmission dynamics in each of the metropolitan statistical areas (MSAs) surrounding Dallas, Houston, New York City, and Phoenix in 2020 and 2021, we extended a previously reported compartmental model accounting for effects of multiple distinct periods of non-pharmaceutical interventions by adding consideration of vaccination and Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) variants Alpha (lineage B.1.1.7) and Delta (lineage B.1.617.2). For each MSA, we found region-specific parameterizations of the model using daily reports of new COVID-19 cases available from January 21, 2020 to October 31, 2021. In the process, we obtained estimates of the relative infectiousness of Alpha and Delta as well as their takeoff times in each MSA (the times at which sustained transmission began). The estimated infectiousness of Alpha ranged from 1.1x to 1.4x that of viral strains circulating in 2020 and early 2021. The estimated relative infectiousness of Delta was higher in all cases, ranging from 1.6x to 2.1x. The estimated Alpha takeoff times ranged from February 1 to February 28, 2021. The estimated Delta takeoff times ranged from June 2 to June 26, 2021. In conclusion, estimated takeoff times are consistent with genomic surveillance data.

60 APPLIED LIFE SCIENCES↗

Phase field modeling of dislocations and obstacles in InSb

We present a phase-field dislocation dynamics (PFDD) model informed by first-principle calculations to elucidate the competitive dislocation nucleation and propagation between the glide and shuffle sets in InSb diamond cubic crystal. The calculations are directly informed with generalized stacking fault energy curves on the (111) slip plane for both the “glide set,” with the smaller interplanar spacing, and the “shuffle set,” with the larger interplanar spacing. The formulation also includes elastic anisotropy and the gradient term associated with the dislocation core. The PFDD calculations show that under no stress the equilibrium structure of screw glide set dislocations dissociates into Shockley partials, while those of the shuffle set dislocations do not dissociate, remaining compact. The calculated dislocation core widths of these InSb dislocations agree well with the measured values for other semiconductor materials, such as Si and GaN. We find that a shuffle set dislocation emits from a dislocation source at an applied stress about three times smaller than that needed to emit leading and trailing partials successively on the glide set plane. Once the partial dislocations in the glide set are emitted, they propagate faster than the shuffle set perfect dislocation at the same stress level.

36 MATERIALS SCIENCE↗

Staged progression epidemic models for the transmission of invasive nontyphoidal Salmonella (iNTS) with treatment

We develop and analyze a stage-progression compartmental model to study the emerging invasive nontyphoidal Salmonella (iNTS) epidemic in sub-Saharan Africa. iNTS bloodstream infections are often fatal, and the diverse and non-specific clinical features of iNTS make it difficult to diagnose. We focus our study on identifying approaches that can reduce the incidence of new infections. In sub-Saharan Africa, transmission and mortality are correlated with the ongoing HIV epidemic and severe malnutrition. We use our model to quantify the impact that increasing antiretroviral therapy (ART) for HIV infected adults and reducing malnutrition in children would have on mortality from iNTS in the population. We consider immunocompromised subpopulations in the region with major risk factors for mortality, such as malaria and malnutrition among children and HIV infection and ART coverage in both children and adults. We parameterize the progression rates between infection stages using the branching probabilities and estimated time spent at each stage. We interpret the basic reproduction number R 0 as the total contribution from an infinite infection loop produced by the asymptomatic carriers in the infection chain. The results indicate that the asymptomatic HIV+ adults without ART serve as the driving force of infection for the iNTS epidemic. We conclude that the worst disease outcome is among the pediatric population, which has the highest infection rates and death counts. Our sensitivity analysis indicates that the most effective strategies to reduce iNTS mortality in the studied population are to improve the ART coverage among high-risk HIV+ adults and reduce malnutrition among children.

59 BASIC BIOLOGICAL SCIENCES↗

Performance evaluation of an air conditioning system based on quasi isothermal dehumidifcation

This work presents a first study of novel air conditioning system based on quasi isothermal dehumidification combined with dew-point evaporative cooling through the Maisotsenko cycle (M-Cycle). The dehumidification process is realized in desiccant coated heat exchanger with heat rejection system. The key idea behind the solution is to maximize the cooling effectiveness by creating favourable conditions for dehumidification and evaporative cooling process. This is achieved by rejecting the heat of sorption in the desiccant system and providing cold and dry air for M-Cycle air cooler, which allows the system to be powered by using low-grade heat energy at 50 °C. The analysis is performed on the basis of experimentally validated numerical models of each component combined in a unified model of the system. The goal of the study is to determine the performance of the solution under different operational conditions. In addition two types of M-Cycle air coolers (i.e. cross-flow and regenerative unit) are compared to establish which unit is more suitable as the source of cooling energy for the system. It was established that proposed system is able to provide comfortable room parameters (supply air temperature equal 17 °C and humidity ratio equal 0.010 kg/kg-da) for outdoor air temperatures ranging from 28 to 36 °C and outdoor air humidity ranging from 0.010 to 0.018 kg/kg-da. This means that the proposed solution is able to provide comfortable conditions in temperate and humid climates. Results also indicate that the system consisting of the regenerative M-Cycle HMX allows better cooling load management, offers higher thermal Coefficient of Performance (COP th ) and uses less water compared to the system with the cross-flow HMX. However, the system with cross-flow heat and mass exchanger (HMX) has a relatively greater electrical Coefficient of Performance (COP el ) in most cases.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗