Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation models”

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,261 records · Page 70

Technical and economic feasibility of an integrated ethanol and anthocyanin coproduction process using purple corn stover

The coproduction of high-value anthocyanin extract in the cellulosic ethanol process would diversify the co-product market, increase revenue, and potentially improve the economics of the process. The high anthocyanin concentration in the cob and structural carbohydrates in residual stover make purple corn stover an attractive source for anthocyanin and ethanol coproduction. In this study, we aimed to develop simulation models for processes integrating ethanol production and anthocyanin extraction using purple corn stover, to evaluate their techno-economic feasibility, and to compare their performance with the conventional ethanol production process using corn stover. The annual ethanol production for plants processing 2000 MT dry feedstock / day was 148.6 million L/year for the integrated processes compared with 222.6 million L/year for the conventional process. Anthocyanin production in the modified processes using dilute acid-based and water-based anthocyanin extraction processes was 1779 and 1099 MT/year, respectively. Capital investments for the integrated processes ($\$448.1$ to $\$443.8$ million) were higher than the conventional process ($\$371.9$ million). Due to high revenue from anthocyanin extract, the ethanol production cost for the integrated process using acid-based anthocyanin extraction ($\$0.36$/L) was 34.5% lower than conventional ethanol production ($\$0.55$/L). The ethanol production cost for the integrated process using water-based anthocyanin extraction ($\$0.68$/L) was higher than conventional ethanol production due to low ethanol and anthocyanin yields. The minimum ethanol selling price for the integrated process using acid-based anthocyanin extraction ($\$0.65$/L) was also lower than the conventional process ($\$0.72$/L), indicating an improvement in economic performance.

09 BIOMASS FUELS↗

Boosting the Oxygen Reduction Performance via Tuning the Synergy between Metal Core and Oxide Shell of Metal−Organic Frameworks‐Derived Co@CoO x

Abstract Co core@Co oxide shell (Co@CoO x ) catalysts represent a large family with promising oxygen reduction reaction (ORR) catalytic activity. However, inadequate understanding of Co@CoO x synergy prohibits further pursuit of catalytic performance enhancement. Herein, a Co zeolitic−imidazolate framework was converted into metallic Co, followed by controlled air treatment to form Co@CoO x . The composition and structure evolution as a function of air treatment temperature were studied thoroughly through conventional and synchrotron (both ex‐situ and in‐situ) characterizations, confirming the coexistence of CoO and Co 3 O 4 in the shell. The optimal catalyst showed an ORR half‐wave potential of 0.87 V (vs. RHE) in an alkaline half‐cell and delivered high discharge capacity in an aprotic Li−O 2 battery (7,124 mAh g Cat+C −1 ) and an aqueous Zn−air battery (694 mAh g Zn −1 ) with good performance retention after durability test. Modeling simulation and density functional theory calculation confirmed the charge donation from metal core to oxide shell and shed light on new insights of how metal@metal oxide synergy impacted the ORR via tuning the charge conductivity, oxygen affinity and intermediate transfer pathway. This work opens up a venue to boost ORR catalytic activity from an interfacial synergy perspective.

Wang, Hao↗

Carbon budget of the Harvard Forest Long-Term Ecological Research site: pattern, process, and response to global change

How, where, and why carbon (C) moves into and out of an ecosystem through time are long-standing questions in biogeochemistry. Here, we bring together hundreds of thousands of C-cycle observations at the Harvard Forest in central Massachusetts, USA, a mid-latitude landscape dominated by 80–120-yr-old closed-canopy forests. These data answered four questions: (1) where and how much C is presently stored in dominant forest types; (2) what are current rates of C accrual and loss; (3) what biotic and abiotic factors contribute to variability in these rates; and (4) how has climate change affected the forest’s C cycle? Harvard Forest is an active C sink resulting from forest regrowth following land abandonment. Soil and tree biomass comprise nearly equal portions of existing C stocks. Net primary production (NPP) averaged 680–750 g C·m –2 ·yr –1 ; belowground NPP contributed 38–47% of the total, but with large uncertainty. Mineral soil C measured in the same inventory plots in 1992 and 2013 was too heterogeneous to detect change in soil-C pools; however, radiocarbon data suggest a small but persistent sink of 10–30 g C·m –2 ·yr –1 . Net ecosystem production (NEP) in hardwood stands averaged ~300 g C·m –2 ·yr –1 . NEP in hemlock-dominated forests averaged ~450 g C·m –2 ·yr –1 until infestation by the hemlock woolly adelgid turned these stands into a net C source. Since 2000, NPP has increased by 26%. For the period 1992–2015, NEP increased 93%. The increase in mean annual temperature and growing season length alone accounted for ~30% of the increase in productivity. Interannual variations in GPP and NEP were also correlated with increases in red oak biomass, forest leaf area, and canopy-scale light-use efficiency. Compared to long-term global change experiments at the Harvard Forest, the C sink in regrowing biomass equaled or exceeded C cycle modifications imposed by soil warming, N saturation, and hemlock removal. Furthermore, the results of this synthesis and comparison to simulation models suggest that forests across the region are likely to accrue C for decades to come but may be disrupted if the frequency or severity of biotic and abiotic disturbances increases.

54 ENVIRONMENTAL SCIENCES↗

Trade-offs of forest management scenarios on forest carbon exchange and threatened and endangered species habitat

Management of forest lands under climate warming poses challenges to managers, some of which are difficult to predict. Examining the trade-offs associated with forest stewardship choices is essential to avoid consequences with loss of natural capital. We utilized LANDIS-II process model simulations for three forested sites located in disparate parts of the United States with the purpose of understanding the trade-offs imposed by management choices under climate warming and associated wildfire. There were only small trade-offs that emerged from the simulations, between habitat area for threatened and endangered species (TES), net ecosystem exchange for CO2 (NEE), and risk of wildfire. Stand management in the form of thinning and prescribed burning typically increased NEE while simultaneously increasing habitat for TES, while reducing the risk of wildfire. These benefits were also observed under a climate warming scenario; however, the benefits were greatly outweighed by the negative impacts f warming on both TES habitat and NEE. Balancing these ecosystem services via thinning and burning treatments is a strategic approach to mitigate risks of wildfire both currently and under a warming future climate.

54 ENVIRONMENTAL SCIENCES↗

Machine Learning Techniques for Data Reduction of Climate Applications

Scientists conduct large-scale simulations to compute derived quantities-of-interest (QoI) from primary data. Often, QoI are linked to specific features, regions, or time intervals, such that data can be adaptively reduced without compromising the integrity of QoI. For many spatiotemporal applications, these QoI are binary in nature and represent presence or absence of a physical phenomenon. We present a pipelined compression approach that first uses neural-network-based techniques to derive regions where QoI are highly likely to be present. Then, we employ a Guaranteed Autoencoder (GAE) to compress data with differential error bounds. GAE uses QoI information to apply low-error compression to only these regions. This results in overall high compression ratios while still achieving downstream goals of simulation or data collections. Experimental results are presented for climate data generated from the E3SM Simulation model for downstream quantities such as tropical cyclone and atmospheric river detection and tracking. These results show that our approach is superior to comparable methods in the literature.

Li, Xiao [University of Florida]↗

Phosphorus recycling by mineral-catalyzed ribonucleotide cleavage on iron and manganese oxides

Abstract Phosphorus is an essential element influencing both food security via plant fertilization, and water pollution through excessive phosphorus use, yet the phosphorus cycle in ecosystems is poorly known. In particular, beyond adsorption, the role of iron and manganese oxides in catalyzing the abiotic dephosphorylation of biomolecules is debated. Here, we studied the reactions of ribonucleotides, containing different phosphate bonding, with goethite, hematite, and birnessite. We employed both high-resolution mass spectrometry of solution species and molecular modeling simulations of ribonucleotide-mineral complexes. Results disclose an up to fivefold preferential hydrolytic cleavage of a phosphoanhydride bond over a phosphoester bond, indicating that mineral-catalyzed reactions reflect the hierarchy reported for the activity of phosphatase enzymes. The fourfold higher catalytic reactivity of goethite and birnessite versus hematite is explained by mineral-specific binding rather than surface area differences. Corresponding simulated adsorbate conformations at the water–mineral interfaces are proposed. Overall, our findings provide new insights on the catalytic recycling of organic phosphorus species by mineral oxides.

58 GEOSCIENCES↗

Climate science to inform adaptation policy: Heat waves over India in the 1.5°C and 2°C warmer worlds

Developing a better scientific understanding of anthropogenic climate change and climate variability, especially the prediction/projection of climate futures with useful temporal and geographical resolution and quantified uncertainties, and using that knowledge to inform adaptation planning and action will become crucially important in the coming years. Generating such policy-relevant knowledge may be particularly important for developing countries such as India. It is with this backdrop that, in this paper, we analyze future heat waves in India by using observations and a large number of model simulations of historical, + 1.5 °C, and + 2.0 °C warmer worlds. In both the future scenarios, there is an increased probability of heat waves during June and July when the Indian monsoon is in full swing and humidity is high, which makes the heat events even more of a health risk. While the highest temperatures in heat waves may not increase much in future climates, the duration and areal extent of the heat waves will most likely increase, leading to the emergence of new heat wave-prone zones in India. The results indicate that the joint frequencies of the longest duration and large area events could be nearly threefold greater in the + 1.5 °C and fivefold greater in the + 2.0 °C future scenarios compared to historical simulations. Thus, overall, the study indicates a substantial increase in the risk of heat events that typically elicit warnings from forecasters. The likely widespread and persistent nature of heat wave events in the future, as revealed by this study, will require planning and adaptation measures beyond the short-term disaster planning frameworks currently in place. Exploring what these measures might look like is beyond the scope of this study, but it underlines the importance of developing climate knowledge with high temporal and geographical resolution capable of informing adaptation policy and planning.

+1.5 °C and +2.0 °C↗

Characterization of Fe-6Si Soft Magnetic Alloy Produced by Laser-Directed Energy Deposition Additive Manufacturing

Abstract Commercial electrical steels, Fe-Si alloys with < 4 wt.% Si, are inexpensive and efficient materials for electrical power conversion. Further efficiency improvements require increasing the silicon concentration to 6 wt.%, at which point the material becomes brittle and difficult to form by conventional rolling and sheet fabrication methods. Additive manufacturing stands to overcome challenges with commercial manufacturing techniques by leveraging near-net-shape fabrication. The wide array of process conditions provides additive manufacturing with increased flexibility, enabling control over the microstructure and mechanical properties. This work explores the microstructures and magnetic properties of ring-shaped Fe-Si alloys produced using concentric and cross-hatch tool paths on a laser-directed energy deposition additive manufacturing system. Concentric-built samples exhibit elongated grain structures while cross-hatch-built samples comprise lower aspect ratio grain structures. Thermal finite element analysis simulations model the stress conditions produced by the different scan path geometries. Microhardness measurements probe the mechanical properties as a function of anneal temperature, providing a qualitative understanding of the intergranular defect density. Soft magnetic properties measured under quasistatic and AC conditions show frequency- and microstructure-dependent coercivity and permeability. Finally, analysis of the core loss quantifies how the build strategies and thermal treatments influence efficiency in electrical power conversion applications. Understanding the influences of scan path geometry and thermal treatment provides a pathway towards application of additively manufactured soft magnetic materials.

Adamczyk, Jesse M.↗

Reassessing the Climate Change Narrative

Abstract We note that the atmosphere has distinct tropical and extratropical regimes. The tropical regime is significantly dependent on the greenhouse effect and is characterized by temperatures that are largely horizontally homogenized. The extratropical regime is dominated by large scale unstable convective eddies that transport heat between the tropics and the poles (leaving the poles warmer than they otherwise would be) and serve to determine the temperature difference between the tropics and the poles. Changes in tropical temperature and in the tropics-to-pole temperature difference both contribute to changes in global mean temperature. It turns out that changes in global mean temperature associated with major climate change (i.e., the last glacial maximum and the warm period of the Eocene about 50 million years ago) were associated primarily with changes in the tropics-to-pole temperature differences. By contrast, changes in global mean temperature over the past 150 years or so are almost entirely associated with changes in tropical temperature. Thus, there is no intrinsic amplification associated with a change in the tropics-to-pole temperature difference. However, model simulations of climate behave differently from both observations and from each other. In particular, they all show more significant contributions for the tropics-to-pole temperature difference – sometimes much more significant. They also show excessive tropical warming.

Meteorology & Atmospheric Sciences↗

Shortened Duration of Global Warming Slowdowns with Elevated Greenhouse Gas Emissions

Continuous emissions of anthropogenic greenhouse gases (GHGs) and aerosols in the last 160 years have resulted in an increasing trend of global mean surface temperatures (GMSTs). Due to interactions with natural variability, rates of the combined anthropogenically and naturally induced warming trends are characterized by significant slowdowns and speedups on decadal timescales. In this work, by analyzing observed and model-simulated data, we investigate how the duration of these episodes will change with different strengths of GHG and aerosol forcing. We found that the duration of warming slowdowns can be more than 30 yr with a slower rate of anthropogenic emissions but would shorten to about 5 yr with a higher one. This duration reduction depends on both the magnitude of the climate response to anthropogenic forcing and the strength of the internal variability. Moreover, the warming slowdowns can still occur even towards the end of this century under high emissions scenarios but with significantly shortened duration.

54 ENVIRONMENTAL SCIENCES↗

Sensitivity Analysis of an Electrospray Dehumidification System

Electrically enhanced condensation achieved by an electrospray system can enhance dehumidification in HVAC applications if successfully scaled up. The first step towards scale-up is identifying critical design and operating parameters that significantly impact dehumidification for the electrospray system. Sensitivity analysis is one approach for identifying critical parameters. This paper presents a methodology that can efficiently perform quantitative sensitivity analysis on a computationally expensive high fidelity simulation model, which is a computational fluid dynamics (CFD) model of the electrospray system. The methodology employs surrogate modeling and Sobol sensitivity analysis. We estimate the main and total effect Sobol sensitivity indices of four operating parameters of the electrospray system. The results reveal that the most important parameter is the size of the sprayed droplets, followed by the volumetric flow rate of the bulk air and injected water in the system. The water injection velocity of the spray was found to have little to no impact on dehumidification.

Young, David↗

Feature extraction for subtle anomaly detection using semi-supervised learning

The demand for automated and effective monitoring techniques has soared with the increased digitization of industrial monitoring systems. State-of-the-art machine learning methods are effectively detecting abrupt changes in system states. However, these methods lack comparable maturity in detecting subtle changes that may be signs of incipient faults. This manuscript argues that the current anomaly detection methods can be enhanced by exploring weak patterns to enable subtle variation detection. Specifically, the concept of semi-supervised learning is employed, with labels representing knowledge about some anomalous conditions of a system. The basic idea is to extract a candidate set of weak patterns discarded by state-of-the-art baselining algorithms. With few labeled anomalous data, the algorithm selects the weak patterns and allows for their possible fusion using the highest sensitivity to the labeled anomalies. Here, the method’s applicability is demonstrated using a representative pressurized water reactor (PWR) model simulated by Dymola.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Fuel performance analysis of Cr-coated Zircaloy-4 cladding during a prototypical LOCA event using BISON

Here, deformation and failure of chromium (Cr) coated Zircaloy-4 (Zry-4) were studied in loss-of-coolant accident (LOCA) conditions using the BISON fuel performance code. The BISON validation model simulating Halden research reactor experiments was extended to include Cr coatings and higher rod internal pressures to simulate high-burnup fuel. The transient simulations show Cr coatings help relieve stress in the Zry-4 substrate during the transient, delaying the onset of high-temperature creep, which leads to ballooning and bursting of the cladding. A nominal Cr coating thickness of 30 µm delays clad failure by 26 seconds and increases the clad burst temperature by 40 K. A parametric study showed that time to failure and burst temperature both increase with coating thickness, and Cr-coated cladding offers burst resistance under a wide range of rod internal pressures simulating high-burnup fuel. Results indicate that a thin Cr coating provides resistance against ballooning and bursting of cladding during LOCA events.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Analyzing at-scale distribution grid response to extreme temperatures

Threats against power grids continue to increase, as extreme weather conditions and natural disasters (extreme events) become more frequent. Hence, there is a need for the simulation and modeling of power grids to reflect realistic conditions during extreme events conditions, especially distribution systems. Herein, this paper presents a modeling and simulation platform for electric distribution grids which can estimate overall power demand during extreme weather conditions. The presented platform's efficacy is shown by demonstrating estimation of electrical demand for 1) Electricity Reliability Council of Texas (ERCOT) during winter storm Uri in 2021, and 2) alternative hypothetical scenarios of integrating Distributed Energy Resources (DERs), weatherization, and load electrification. In comparing to the actual demand served by ERCOT during the winter storm Uri of 2021, the proposed platform estimates approximately 34 GW of peak capacity deficit. These numbers are consistent with state-of-the-art prediction results published in the literature. For the case of the future electrification of heating loads, peak capacity of 78 GW (124% increase) is estimated, which would be reduced to 47 GW (38% increase) with the adoption of efficient heating appliances and improved thermal insulation. Integrating distributed solar PV and storage into the grid causes improvement in the local energy utilization and hence reduces the potential unmet energy by 31% and 40%, respectively.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Adaptive primal–dual control for distributed energy resource management

With the increased adoption of distributed energy resources (DERs) in distribution networks, their coordinated control with a DER management system (DERMS) that provides grid services (e.g., voltage regulation, virtual power plant) is becoming more necessary. One particular type of DERMS using primal–dual control has recently been found to be very effective at providing multiple grid services among an aggregation of DERs; however, the main parameter, the primal–dual step size, must be manually tuned for the DERMS to be effective, which can take a considerable amount of engineering time and labor. To this end, we design a simple method that self-tunes the step size(s) and adapts it to changing system conditions. Additionally, it gives the DER management operator the ability to prioritize among possibly competing grid services. Here we evaluate the automatic tuning method on a simulation model of a real-world feeder in Colorado with data obtained from an electric utility. Through a variety of scenarios, we demonstrate that the DERMS with automatically and adaptively tuned step sizes provides higher-quality grid services than a DERMS with a manually tuned step size.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Observational study of relationships between entrainment rate, homogeneity of mixing, and cloud droplet relative dispersion

The entrainment–mixing process is key to the cloud’s dynamical and microphysical processes. However, to the best of our knowledge, no study has quantitatively linked entrainment rate and homogeneity of mixing. Additionally, their combined effects on cloud droplet relative dispersion and physical mechanisms are even more unclear. This study analyzed the aircraft measurements of cumulus clouds collected in the Holistic Interactions of Shallow Clouds, Aerosols, and Land-Ecosystems field campaign to identify these relationships. The relationship between the entrainment rate and homogeneous mixing degree was found to be moderately positive, and the underlying mechanisms related to mixing and evaporation time scales and environmental relative humidity were explored. The relationship between relative dispersion and entrainment rate changes from positive to negative with increasing entrainment rate, depending on the dominance of evaporation of large or small droplets. The relative dispersion first increases and then decreases with increasing homogeneous mixing degree, validating the results of previous model simulations. Further, the results shed new light on understanding the interactions between cloud microphysics and cloud dynamics through turbulent entrainment and mixing processes.

54 ENVIRONMENTAL SCIENCES↗

Molecular dynamics study of carbon dioxide and nitrogen selectivity through poly[bis((methoxyethoxy)ethoxy)phosphazene] (MEEP) membrane

Here, a molecular dynamics simulation model was developed to comprehensively understand carbon dioxide over nitrogen (CO 2 /N 2 ) selectivity through polyphosphazene-based membrane comprising of poly[bis((methoxyethoxy)ethoxy)phosphazene] (MEEP) selective layer at the molecular level. The effects of temperature, pressure, and initial feed gas composition on the CO 2 transport on a polymer membrane were studied. The computed free energy and density profile of the permeating gas mixture exhibited that CO 2 molecules express higher interactions with the membrane than N 2 molecules, resulting in higher CO 2 /N 2 selectivity. Statistical analysis of gas molecules (CO 2 , water (H 2 O), and N 2 ) transportation suggested that hydro- and CO 2 -philic functional groups in the membrane significantly impact CO 2 permeability and CO 2 /N 2 selectivity. This study suggested that Lewis acid–base and hydrogen bonding combinations contribute to CO 2 permeation and CO 2 /N 2 selectivity. An equal CO 2 /N 2 selectivity was observed with and without water vapor in the feed gas suggesting that water does not hinder CO 2 transport through the membrane.

36 MATERIALS SCIENCE↗

Intrepid MCMC: Metropolis-Hastings with exploration

In engineering examples, one often encounters the need to sample from unnormalized distributions with complex shapes that may also be implicitly defined through a physical or numerical simulation model, making it computationally expensive to evaluate the associated density function. For such cases, MCMC has proven to be an invaluable tool. Random-walk Metropolis Methods (also known as Metropolis-Hastings (MH)), in particular, are highly popular for their simplicity, flexibility, and ease of implementation. However, most MH algorithms suffer from significant limitations when attempting to sample from distributions with multiple modes (particularly disconnected ones). Here, in this paper, we present Intrepid MCMC - a novel MH scheme that utilizes a simple coordinate transformation to significantly improve the mode-finding ability and convergence rate to the target distribution of random-walk Markov chains while retaining most of the simplicity of the vanilla MH paradigm. Through multiple examples, we showcase the improvement in the performance of Intrepid MCMC over vanilla MH for a wide variety of target distribution shapes. We also provide an analysis of the mixing behavior of the Intrepid Markov chain, as well as the efficiency of our algorithm for increasing dimensions. A thorough discussion is presented on the practical implementation of the Intrepid MCMC algorithm. Finally, its utility is highlighted through a Bayesian parameter inference problem for a two-degree-of-freedom oscillator under free vibration.

97 - MATHEMATICS AND COMPUTING↗