SEARCH · Search NASA
Results for “Simulation System”
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.
Utilization of Hydroturbines in Wastewater Treatment Plants
Wastewater treatment plants (WWTPs) are a significant energy consumer, yet there are several opportunities for implementing on-site power generation systems. Within the treatment process, the high flow rate of effluent is produced and discharged to a nearby water body by gravity. Thus, hydroturbines can be utilized to generate power in such an application due to a difference in elevation and high flow rate. This paper presents a case study of introducing a hydroturbine in a WWTP in Wisconsin and evaluating the power output in addition to determining the energy savings. The WWTP considered in this study has an effluent flow rate of 190 MGD (million gallons per day) and elevation difference of 3 m (10 ft) between the final stage of treatment and the discharge point. Based on the parameters above; hubless rim-drive Kaplan type hydroturbine (RDT) is the optimal choice to be used in such an application. The RDT is designed and optimized by using in-house code. A computational fluid dynamics (CFD) software is applied to evaluate the performance of the proposed model, and the system is simulated through homer software to validate the results generated by the CFD. The expected savings is estimated to be 1564 MWh/yr.
Situating the Vector Density Approach Among Contemporary Continuum Theories of Dislocation Dynamics
Abstract For the past century, dislocations have been understood to be the carriers of plastic deformation in crystalline solids. However, their collective behavior is still poorly understood. Progress in understanding the collective behavior of dislocations has primarily come in one of two modes: the simulation of systems of interacting discrete dislocations and the treatment of density measures of varying complexity that are considered as continuum fields. A summary of contemporary models of continuum dislocation dynamics is presented. These include, in order of complexity, the two-dimensional statistical theory of dislocations, the field dislocation mechanics treating the total Kröner–Nye tensor, vector density approaches that treat geometrically necessary dislocations on each slip system of a crystal, and high-order theories that examine the effect of dislocation curvature and distribution over orientation. Each of theories contain common themes, including statistical closure of the kinetic dislocation transport equations and treatment of dislocation reactions such as junction formation. An emphasis is placed on how these common themes rely on closure relations obtained by analysis of discrete dislocation dynamics experiments. The outlook of these various continuum theories of dislocation motion is then discussed.
Toward Polydisperse Flows With MFIX-EXA
In the presence of large size disparities, single-grid neighbor search algorithms lead to inflated neighbor lists that significantly degrade the performance of Lagrangian particle solvers. If Eulerian–Lagrangian (EL) frameworks are to remain performant when simulating realistic systems, improved neighbor detection approaches must be adopted. To this end, we consider the application of a multigrid neighbor search (MGNS) algorithm in the mfix-exa software package, an exascale EL solver built upon the AMReX library. Here, details regarding the implementation and verification of MGNS are provided along with speedup curves for a bidisperse mixing layer. MGNS is shown to yield up to 15$\times$ speedup on CPU and 6$\times$ speedup on GPU for the problems considered here. The mfix-exa software is then validated for a variety of polydisperse flows. Finally, a brief discussion is given for how dynamic MGNS may be completed, with application to spatially varying particle size distributions.
Energy Impact of Radiative Cooling Paints in Warehouses Under Various United States Climates
Although radiative cooling research is widely found in the literature, no comprehensive study has yet been conducted on the impact of novel radiant cooling (>0.91 reflectance) on the energy efficiency of warehouses. Here, in this work, we develop three building models based on a Department of Energy prototype warehouse model using trnsys, representing a typical warehouse with a black roof, a typical warehouse with a white roof, and a warehouse with novel radiative cooling (RC) paint on its roof. These models are run for 15 different cities, each representative of a different ASHRAE climate zone, to better understand the impact of RC in many different climates. It was found that an RC-coated roof in a warehouse could reduce the building's annual heating, ventilation, and air conditioning (HVAC) loads by up to 14.11 kWh/m 2 of the roof area compared to a black roof, resulting in a maximum reduction in energy costs of 0.55 $\$$/m 2 or $\$$2646/year for a large 4835 m 2 warehouse. Similarly, replacing the typical white roof coating with an RC coating could reduce the warehouse's energy consumption by up to 8.17 kWh/ m 2 of roof area, thus reducing energy costs by as much as 0.29 $\$$/m 2 or $\$$1386/year for a 4835 m 2 warehouse. In addition, applying RC paint to an unconditioned warehouse could reduce the building's ASHRAE Standard 55 indoor temperature exceedance by up to 1330 h/year compared to a black roof and up to 532 h/year compared to a white roof.
CONTROLLED END GAS AUTO IGNITION WITH EXHAUST GAS RECIRCULATION ON A STOICHIOMETRIC, SPARK IGNITED, NATURAL GAS ENGINE
The goal of this study is to address fundamental limitations to achieving diesel-like efficiencies in heavy duty on-highway natural gas (NG) engines. Engine knock and misfire are barriers to pathways leading to higher efficiency engines. This study explores enabling technologies for development of high efficiency stoichiometric, spark ignited, natural gas engines. These include design strategies for fast and stable combustion and higher dilution tolerance. Additionally, advanced control methodologies are implemented to maintain stable operation between knock and misfire limits. To implement controlled end- gas autoignition (C-EGAI) strategies a Combustion Intensity Metric (CIM) is used for ignition control with the use of a Woodward large engine control module (LECM). Tests were conducted using a single cylinder, variable compression ratio, cooperative fuel research (CFR) engine with baseline conditions of 900 RPM, engine load of 800 kPa indicated mean effective pressure (IMEP), and stoichiometric air/fuel ratio. Exhaust gas recirculation (EGR) tests were performed using a custom EGR system that simulates a high pressure EGR loop and can provide a range of EGR rates from 0 to 40%. The experimental measurements included the variance of EGR rate, compression ratio, engine speed, IMEP, and CIM. These five variables were optimized through a Modified Box- Benken design Surface Response Method (RSM), with brake efficiency as the merit function. A positive linear correlation between CIM and f-EGAI was identified. Consequently, CIM was used as the feedback control parameter for C-EGAI. As such, implementation of C-EGAI effectively allowed for the utilization of high EGR rates and CRs, controlling combustion between a narrower gap between knock and lean limits. The change from fixed to parametric ignition timing with CIM targeted select values of f-EGAI with an average coefficient of variance (COV) of peak pressure of 5.4. The RSM efficiency optimization concluded with operational conditions of 1080 RPM, 1150 kPa IMEP, 10.55:1 compression ratio, and 17.8% EGR rate with a brake efficiency of 21.3%. At this optimized point of peak performance, a f-EGAI for C-EGAI was observed at 34.1% heat release due to auto ignition, a knock onset crank angle value of 10.3° aTDC and ignition timing of - 24.7° aTDC. This work has demonstrated that combustion at a fixed f-EGAI can be maintained through advanced ignition control of CIM without experiencing heavy knocking events.
Bio-inspired Bistable Piezoelectric Energy Harvester for Powering Animal Telemetry Tags: Concept Design and Preliminary Experimental Validation
This paper presents the concept design, preliminary experimental validation, and performance evaluation of a novel bio-inspired bi-stable piezoelectric energy harvester for self-powered animal telemetry tags. The overall concept design, which includes a bio-inspired attachment and a bi-stable piezoelectric energy harvester, is introduced firstly with a specific application example of marine fish tracking. The self-powered telemetry tag can be externally deployed on fish (dorsal fin) to monitor fish habitats, population, and underwater environment. Inspired by the Venus flytrap’s rapid shape transition, a bistable piezoelectric energy harvester is developed to scavenge energy from fish maneuvering and the surrounding fluid flow for a sustainable power supply. The bistability of the harvester is characterized by the measured force-displacement curve and double potential wells. A bluff body is integrated to the free end of the bistable piezoelectric energy harvester to enhance the structure-fluid interaction for the large-amplitude snap-through vibrations and higher voltage output. Controlled laboratory experiments are conducted in a water tank on the bio-inspired bi-stable piezoelectric energy harvester using a servo motor system to simulate fish swing motion at various conditions to evaluate the power generation performance. The preliminary underwater experimental results demonstrated that the proposed bio-inspired bi-stable piezoelectric could effectively converter fish swing motions into electricity. 17.25 mJ energy was obtained within 130 s under the peak-to-peak swing angle of 30o at 1.5 Hz in the capacitor charging experiments.
WE-G-204-01: BEST IN PHYSICS (IMAGING): Effect of Image Processing Parameters On Nodule Detectability in Chest Radiography
Purpose: To investigate the effect of varying system image processing parameters on lung nodule detectability in digital radiography. Methods: An anthropomorphic chest phantom was imaged in the posterior-anterior position using a GE Discovery XR656 digital radiography system. To simulate lung nodules, a polystyrene board with 6.35mm diameter PMMA spheres was placed adjacent to the phantom (into the x-ray path). Due to magnification, the projected simulated nodules had a diameter in the radiographs of approximately 7.5 mm. The images were processed using one of GE’s default chest settings (Factory3) and reprocessed by varying the “Edge” and “Tissue Contrast” processing parameters, which were the two user-configurable parameters for a single edge and contrast enhancement algorithm. For each parameter setting, the nodule signals were calculated by subtracting the chest-only image from the image with simulated nodules. Twenty nodule signals were averaged, Gaussian filtered, and radially averaged in order to generate an approximately noiseless signal. For each processing parameter setting, this noise-free signal and 180 background samples from across the lung were used to estimate ideal observer performance in a signal-known-exactly detection task. Performance was estimated using a channelized Hotelling observer with 10 Laguerre-Gauss channel functions. Results: The “Edge” and “Tissue Contrast” parameters each had an effect on the detectability as calculated by the model observer. The CHO-estimated signal detectability ranged from 2.36 to 2.93 and was highest for “Edge” = 4 and “Tissue Contrast” = −0.15. In general, detectability tended to decrease as “Edge” was increased and as “Tissue Contrast” was increased. A human observer study should be performed to validate the relation to human detection performance. Conclusion: Image processing parameters can affect lung nodule detection performance in radiography. While validation with a human observer study is needed, model observer detectability for common tasks could provide a means for optimizing image processing parameters.
SU-F-J-05: The Effect of Air Pockets in the Urinary Bladder During Bladder Hyperthermia Treatment
Purpose: Hyperthermia combined with Mitomycin C is used for the treatment of non-muscle invasive bladder cancer (NMIBC), using a phased array system of microwave antennas for bladder heating. Often some air is present in the bladder, which effectively blocks the microwave radiation, potentially preventing proper treatment of that part of the bladder. Air can be a relevant fraction of the bladder content and large air pockets are expected to have a noticeable influence on achieved temperatures. Methods: We analysed 14 NMIBC patients treated at our institute with our AMC-4 hyperthermia device with four 70MHz antennas around the pelvis. A CT scan was made after treatment and a physician delineated the bladder on the CT scan. On the same scan, the amount of air present in the bladder was delineated. Using our in-house developed hyperthermia treatment planning system, we simulated the treatment using the clinically applied device settings. We did this once with the air pocket delineated on the CT scan, and once with the same volume filled with bladder tissue. Results: The patients had on average 4.2ml (range 0.8–10.1ml) air in the bladder. The bladder volume was delineated by the physician, that is including air pocket and bladder wall, was on average 253ml (range 93–452ml). The average volume in which changes exceeded 0.25°C was 22ml (range 0–108 ml), with the bladder being up to 2°C cooler when an air pocket was present. Except for extreme cases, there was no evident relation between the quantity of air and the difference in temperature. Conclusion: The effect of an air pocket in the bladder during bladder hyperthermia treatment varies strongly between patients. Generally, this leads to lower temperatures in the bladder, potentially affecting treatment quality, and suggesting that care need be taken to minimise the size of air pockets during hyperthermia treatments. The KWF Dutch Cancer Society financially supported this work, grant UVA 2012-5539.
TH-AB-207A-01: Contrast-Enhanced CT: Correlation of Radiation Dose and Biological Effect
Purpose: The potential risk from CT is generally characterized in terms of radiation dose. The presence of iodinated-contrast medium increases radiation dose. However, it is unclear how much of this increase is biologically relevant. The purpose of this study was to establish the contribution of dose increase from iodine to biological effect. Methods: Radiation organ dose was estimated in 58 human (XCAT) phantoms “undergoing” chest CT examination (120 kVp, 9 mGy CTDI) on a simulated CT system (Definition Flash, Siemens) with and without iodinated-contrast agent (62.5 mL of iodine per subject). The dose without and with the presence of iodine was compared to the increase in foci per cell (a surrogate of DNA damage) measured before and after similar CT exams without and with contrast agent (Piechowiak et al. 2015). The data were analyzed to ascertain how the enhancement in biological effect in contrast-enhanced CTs correlated with the increase in dose due to the presence of iodine. Results: The presence of iodinated-contrast in CT increased the organ doses by 2% to 50% on average. Typical values were heart (50%±7%), kidney (19%±7%), and liver (2%±3%). The corresponding increase in the average foci per cell was 107%±19%, indicating biological effect of iodine was greater than what would be anticipated from the iodine-initiated increase in radiation dose alone. Conclusion: Mean foci per cell and organ dose both increase in the presence of contrast agent. The former, however, is at least twice as large as the latter, indicating that iodine contributes to an increase in the probability of DNA damage not only as a consequence of increased x-ray energy deposition but also from other mechanisms. Hence iodine radiation dose, while relevant to be included in estimating the risk associated with contrast-enhanced CT, still can underestimate the biological effects.
Reducing opioid use disorder and overdose deaths in the United States: A dynamic modeling analysis
Opioid overdose deaths remain a major public health crisis. We used a system dynamics simulation model of the U.S. opioid-using population age 12 and older to explore the impacts of 11 strategies on the prevalence of opioid use disorder (OUD) and fatal opioid overdoses from 2022 to 2032. These strategies spanned opioid misuse and OUD prevention, buprenorphine capacity, recovery support, and overdose harm reduction. By 2032, three strategies saved the most lives: (i) reducing the risk of opioid overdose involving fentanyl use, which may be achieved through fentanyl-focused harm reduction services; (ii) increasing naloxone distribution to people who use opioids; and (iii) recovery support for people in remission, which reduced deaths by reducing OUD. Increasing buprenorphine providers’ capacity to treat more people decreased fatal overdose, but only in the short term. Our analysis provides insight into the kinds of multifaceted approaches needed to save lives.
Adaptation of the model of tunneling in a metal/CaF{sub 2}/Si(111) system for use in industrial simulators of MIS devices
An approach toward simplification of the model of the tunneling transport of electrons through a thin layer of crystalline calcium fluoride into a silicon (111) substrate with subsequent implementation in simulators of semiconductor devices is suggested. The validity of the approach is proven by comparing the results of modeling using simplified formulas with the results of precise calculations and experimental data. The approach can be applied to calculations of tunneling currents in structures with any crystalline insulators on Si (111)
The modern temperature-accelerated dynamics approach
Accelerated molecular dynamics (AMD) is a class of MD-based methods used to simulate atomistic systems in which the metastable state-to-state evolution is slow compared with thermal vibrations. Temperature-accelerated dynamics (TAD) is a particularly efficient AMD procedure in which the predicted evolution is hastened by elevating the temperature of the system and then recovering the correct state-to-state dynamics at the temperature of interest. TAD has been used to study various materials applications, often revealing surprising behavior beyond the reach of direct MD. This success has inspired several algorithmic performance enhancements, as well as the analysis of its mathematical framework. Recently, these enhancements have leveraged parallel programming techniques to enhance both the spatial and temporal scaling of the traditional approach. Here, we review the ongoing evolution of the modern TAD method and introduce the latest development: speculatively parallel TAD.
Multi-Fluid and Kinetic Models of Partially Ionized Magnetic Reconnection
Magnetic reconnection in partially ionized plasmas is a ubiquitous and important phenomena in both laboratory and astrophysical systems. Here, simulations of partially ionized magnetic reconnection with well-matched initial conditions are performed using both multi-fluid and fully-kinetic approaches. Despite similar initial conditions, the time-dependent evolution differs between the two models. In multi-fluid models, the reconnection rate locally obeys either a decoupled Sweet-Parker scaling, where neutrals are unimportant, or a fully coupled Sweet-Parker scaling, where neutrals and ions are strongly coupled, depending on the resistivity. In contrast, kinetic models show a faster reconnection rate that is proportional to the fully-coupled, bulk Alfv\'en speed, $v_A^\star$. These differences are interpreted as the result of operating in different collisional regimes. Multi-fluid simulations are found to maintain $\nu_{ni}L/v_A^\star \gtrsim 1$, where $\nu_{ni}$ is the neutral-ion collision frequency and $L$ is the time-dependent current sheet half-length. This strongly couples neutrals to the reconnection outflow, while kinetic simulations evolve to allow $\nu_{ni}L/v_A^\star < 1$, decoupling neutrals from the reconnection outflow. Differences in the way reconnection is triggered may explain these discrepancies.
Geographic Visualization of Power-Grid Dynamics
The visualization enables the simulation analyst to see changes in the frequency through time and space. With this technology, the analyst has a bird's eye view of the frequency at loads and generators as the simulated power system responds to the loss of a generator, spikes in load, and other contingencies. The significance of a contingency to the operation of an electrical power system depends critically on how the resulting tansients evolve in time and space. Consequently, these dynamic events can only be understood when seen in their proper geographic context. this understanding is indispensable to engineers working on the next generation of distributed sensing and control systems for the smart grid. By making possible a natural and intuitive presentation of dynamic behavior, our new visualization technology is a situational-awareness tool for power-system engineers.
Probabilistic Resource Adequacy Suite (PRAS)
The Probabilistic Resource Adequacy Suite (PRAS) provides an open-source, research-oriented collection of tools for analysing the resource adequacy of a bulk power system. The simulation methods offered support everything from classical convolution-based analytical techniques through to high-performance sequential Monte Carlo methods supporting multi-region composite reliability assessment, including simulation of energy-limited resources such as storage.
PathParser [SWR-19-54]
PathParser is a python-based computational tool for the thermodynamics and robustness analysis of native and designed metabolic pathways. The following functionalities are provided: 1 Max-min driving force (MDF) optimization. The Gibbs free energy change of the least favorable reaction will be maximized to evaluate thermodynamic feasibility of the entire pathway. 2 Protein cost estimation. The smallest enzyme investment will be assessed by calculating the minimal total enzyme cost supporting a given pathway flux. 3 Robustness analysis. Specifically, an ensemble of models is generated to simulate the system response to enzyme perturbation based on bifurcation theory and a continuation method. Probability of system failure and flux fold change against enzyme perturbation are estimated as well as a flux control index.
Molecules in Studio v. 1.0
A Powersim Studio implementation of the system dynamics’ ‘Molecules of Structure’. The original implementation was in Ventana’s Vensim language by James Hines. The molecules are fundamental constructs of the system dynamics simulation methodology.