Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation Cases”

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

Modeling the Strategic Behavior of an Active Distribution Network in the ISO Markets

With increasing integration of distributed energy resources (DERs), active distribution networks (ADNs) can actively participate in the electricity markets by dispatching their DERs, which can change the existing electricity market paradigm. It is essential to investigate the strategic behaviors of ADNs and their DER dispatch when they participate in the wholesale market as price-makers. This paper proposes a bi-level optimization model to study the strategic behavior of an ADN in both energy and reserve markets. The optimal scheduling of DERs in the ADN is modeled as the upper level problem and the joint energy and reserve market-clearing of the ISO is modeled as the lower-level problem. The two-level optimization models exchange bidding information and energy/reserve prices with each other. The proposed bi-Ievel optimization problem is converted to a mathematical programming with equilibrium constraints (MPEC) by using Karush-Kuhn Tucker (KKT) conditions and strong duality theory. Further, the MPEC problem is reformulated as a computationally-solvable mixed integer second order cone programming (MISOCP) model. The simulation results on an illustrative case demonstrate the impact of the strategic bidding of the ADN on the day-ahead energy and reserve market prices.

Xue, Yaosuo↗

Active Learning for Metamaterial Optimization on HPC and QC Integrated Systems

Active learning algorithms, integrating machine learning, quantum computing and optics simulation in an iterative loop, offer a promising approach to optimizing metamaterials. However, these algorithms can face difficulties in optimizing highly complex structures due to computational limitations. High-performance computing (HPC) and quantum computing (QC) integrated systems can address these issues by enabling parallel computing. In this study, we develop an active learning algorithm working on HPC-QC integrated systems. We evaluate the performance of optimization processes within active learning (i.e., training a machine learning model, problem-solving with quantum computing, and evaluating optical properties through wave-optics simulation) for highly complex metamaterial cases. Our results showcase that utilizing multiple cores on the integrated system can significantly reduce computational time, thereby enhancing the efficiency of optimization processes. Therefore, we expect that leveraging HPC-QC integrated systems helps effectively tackle large-scale optimization challenges in general.

Kim, Seongmin↗

Process Modeling of Aerosol‐Cloud Interaction in Summertime Precipitating Shallow Cumulus Over the Western North Atlantic

Abstract Process modeling of Aerosol‐cloud interaction (ACI) is essential to bridging gaps between observational analysis and climate modeling of aerosol effects in the Earth system and eventually reducing climate projection uncertainties. In this study, we examine ACI in summertime precipitating shallow cumuli observed during the Aerosol Cloud meTeorology Interactions oVer the western ATlantic Experiment (ACTIVATE). Aerosols and precipitating shallow cumuli were extensively observed with in‐situ and remote‐sensing instruments during two research flight cases on 02 June and 07 June, respectively, during the ACTIVATE summer 2021 deployment phase. We perform observational analysis and large‐eddy simulation (LES) of aerosol effect on precipitating cumulus in these two cases. Given the measured aerosol size distributions and meteorological conditions, LES is able to reproduce the observed cloud properties by aircraft such as liquid water content (LWC), cloud droplet number concentration ( N c ) and effective radius r eff . However, it produces smaller liquid water path (LWP) and larger N c compared to the satellite retrievals. Both 02 and 07 June cases are over warm waters of the Gulf Stream and have a cloud top height over 3 km, but the 07 June case is more polluted and has larger LWC. We find that the N a ‐induced LWP adjustment is dominated by precipitation feedback for the 2 June precipitating case and there is no clear entrainment feedback in both cases. An increase of cloud fraction due to a decrease of aerosol number concentration is also shown in the simulations for the 02 June case.

54 ENVIRONMENTAL SCIENCES↗

Comparison of N2O Emissions from Soils at Three Temperate Agricultural Sites

Nitrous oxide (N2O) flux simulations by four models were compared with year-round field measurements from five temperate agricultural sites in three countries. The field sites included an unfertilized, semi-arid rangeland with low N2O fluxes in eastern Colorado, USA; two fertilizer treatments (urea and nitrate) on a fertilized grass ley cut for silage in Scotland; and two fertilized, cultivated crop fields in Germany where N2O loss during the winter was quite high. The models used were daily trace gas versions of the CENTURY model, DNDC, ExpertN, and the NASA-Ames version of the CASA model. These models included similar components (soil physics, decomposition, plant growth, and nitrogen transformations), but in some cases used very different algorithms for these processes. All models generated similar results for the general cycling of nitrogen through the agro-ecosystems, but simulated nitrogen trace gas fluxes were quite different. In most cases the simulated N20 fluxes were within a factor of about 2 of the observed annual fluxes, but even when models produced similar N2O fluxes they often produced very different estimates of gaseous N loss as nitric oxide (NO), dinitrogen (N2), and ammonia (NH3). Accurate simulation of soil moisture appears to be a key requirement for reliable simulation of N2O emissions. All models simulated the general pattern of low background fluxes with high fluxes following fertilization at the Scottish sites, but they could not (or were not designed to) accurately capture the observed effects of different fertilizer types on N2O flux. None of the models were able to reliably generate large pulses of N2O during brief winter thaws that were observed at the two German sites. All models except DNDC simulated very low N2O fluxes for the dry site in Colorado. The US Trace Gas Network (TRAGNET) has provided a mechanism for this model and site intercomparison. Additional intercomparisons are needed with these and other models and additional data sets; these should include both tropical agro-ecosystems and new agricultural management techniques designed for sustainability.

Frolking, S. E.↗

Reactor Pressure Vessel Fracture Mechanics Development and Concrete Application Testing for Grizzly

The Grizzly code is being developed to address degradation issues in nuclear reactor structures and components. For light-water reactors, Grizzly currently has capabilities to simulate degradation processes and their effects on structural integrity in two key areas: reactor pressure vessels (RPVs) and reinforced concrete structures. This report documents improvements made to Grizzly’s ability to address both of these structural systems. For RPVs, the reduced-order models (ROMs) used in fracture mechanics calculations have been expanded to allow their application over a broader range of the parameter space than was permitted by the previous models. The ROMs currently used in Grizzly for the evaluation of flaws that are fully embedded within the RPV (as opposed to surface-breaking flaws) are based on a model that is known to be conservative, indicating higher stress intensity factors than would be obtained from direct simulations. A more accurate model that eliminates these excess conservatisms has been recently included in the American Society of Mechanical Engineers Boiler and Pressure Vessel Code but was not applicable for flaws near the RPV surface, which is where the most critical flaws are usually located. That model has recently been extended for increased applicability in this near-surface region. The ROMs for embedded flaws in the Grizzly code have been expanded to include these recent extensions, which permit their use in a much broader set of cases than previously possible. Direct 3D simulations have been used to check these ROMs and have shown good agreement in most cases, although there are still some cases that need further investigation. There are considerable benefits to using these these more accurate and less conservative ROMs for embedded flaws. On a benchmark probabilistic fracture mechanics problem tested here, the conditional probability of fracture initiation computed for a population of flaws in a single plate in an RPV decreased by over a factor of 3. To address aging in reinforced concrete structures, a capability to simulate multiple degradation mechanisms, including alkali-silica reaction and radiation-induced volumetric expansion has been developed in Grizzly over the past several years. This had previously been demonstrated on laboratory-scale specimens but not on full-scale nuclear concrete structures with reinforcement. To demonstrate the applicability of Grizzly to the analysis of large-scale structures of interest, a full 3D model of a representative reinforced concrete structure, including a complex arrangement of reinforcing bars, was developed and demonstrated in Grizzly.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Evaluation of the Fitch Wind-Farm Wake Parameterization with Large-Eddy Simulations of Wakes Using the Weather Research and Forecasting Model

Abstract Wind-farm parameterizations in weather models can be used to predict both the power output and farm effects on the flow; however, their correctness has not been thoroughly assessed. We evaluate the wind-farm parameterization of the Weather Research and Forecasting Model with large-eddy simulations (LES) of the wake performed with the same model. We study the impact on the velocity and turbulence kinetic energy (TKE) of inflow velocity, roughness, resolution, number of turbines (one or two), and inversion height and strength. We compare the mesoscale with the LES by spatially averaging the LES within areas correspondent to the mesoscale horizontal spacing: one covering the turbine area and two downwind. We find an excellent agreement of the velocity within the turbine area between the two types of simulations. However, within the same area, we find the largest TKE discrepancies because in mesoscale simulations, the turbine-added TKE has to be highest at the turbine position to be advected downwind. Within the downwind areas, differences between velocities increase as the wake recovers faster in the LES, whereas for the TKE both types of simulations show similar levels. From the various configurations, the impact of inversion height and strength is small for these heights and inversion levels. The highest impact for the one-turbine simulations appears under the low-speed case due to the higher thrust, whereas the impact of resolution is low for the large-eddy simulations but high for the mesoscale simulations. Our findings demonstrate that higher-fidelity simulations are needed to validate wind-farm parameterizations.

17 WIND ENERGY↗

A numerical model simulation of the regional air pollution meteorology of the greater Chesapeake Bay area - Summer day case study

The mesoscale numerical model of the University of Virginia (UVMM), has been applied to the greater Chesapeake Bay area in order to provide a detailed description of the air pollution meteorology during a typical summer day. This model provides state of the art simulations for land-sea thermally induced circulations. The model-predicted results agree favorably with available observed data. The effects of synoptic flow and sea breeze coupling on air pollution meteorological characteristics in this region, are demonstrated by a spatial and temporal presentation of various model predicted fields. A transport analysis based on predicted wind velocities indicated possible recirculation of pollutants back onto the Atlantic coast due to the sea breeze circulation.

Segal, M.↗

High Ice Water Concentrations in the 19 August 2015 Coastal Mesoconvective System

During August 2015, NASA's DC-8 research aircraft was flown into High Ice Water Content (HIWC) events as part of a three-week campaign to collect airborne radar data and to obtain measurements from microphysical probes. Goals for this flight campaign included improved characterization of HIWC events, especially from an airborne radar perspective. This paper focuses on one of the flight days, in which a coastal mesoscale convective system (MCS) was investigated for HIWC conditions. The system appears to have been maintained by bands of convection flowing in from the Gulf of Mexico. These convective bands were capped by a large cloud canopy, which masks the underlying structure if viewed from an infrared sensing satellite. The DC-8 was equipped with an IsoKinetic Probe that measured ice concentrations of up to 2.3 g m(exp -3) within the cloud canopy of this system. Sustained measurements of ice crystals with concentrations exceeding 1 g m(exp -3) were encountered for up to ten minutes of flight time. Airborne Radar reflectivity factors were found to be weak within these regions of high ice water concentrations, suggesting that Radar detection of HIWC would be a challenging endeavor. This case is then investigated using a three-dimensional numerical cloud model. Profiles of ice water concentrations and radar reflectivity factor demonstrate similar magnitudes and scales between the flight measurements and model simulation. Also discussed are recent modifications to the numerical model's ice-microphysics that are based on measurements during the flight campaign. The numerical model and its updated ice-microphysics are further validated with a simulation of a well-known case of a supercell hailstorm measured during the Cooperative Convective Precipitation Experiment. Differences in HIWC between the continental supercell and the coastal MCS are discussed.

Proctor, Fred H.↗

Overcoming timestep limitations in boosted-frame particle-in-cell simulations of plasma-based acceleration

Explicit electromagnetic Particle-In-Cell (PIC) codes are typically limited by the Courant- Friedrichs-Lewy (CFL) condition, which implies that the timestep multiplied by the speed of light must be smaller than the smallest cell size. In the case of boosted-frame PIC simulations of plasma-based acceleration, this limitation can be a major hinderance as the cells are often very elongated along the longitudinal direction and the timestep is thus limited by the small, transverse cell size. This entails many small-timestep PIC iterations, and can limit the potential speed-up of the boosted-frame technique. Here, by using a CFL-free analytical spectral solver, and by mitigating additional numerical instabilities that arise at large timestep, we show that it is possible to overcome traditional limitations on the timestep and thereby realize the full potential of the boosted-frame technique over a much wider range of parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Experiments with the Mesoscale Atmospheric Simulation System (MASS) using the synthetic relative humidity

This study is intended to examine the impact of the synthetic relative humidity on the model simulation of mesoscale convective storm environment. The synthetic relative humidity is derived from the National Weather Services surface observations, and non-conventional sources including aircraft, radar, and satellite observations. The latter sources provide the mesoscale data of very high spatial and temporal resolution. The synthetic humidity data is used to complement the National Weather Services rawinsonde observations. It is believed that a realistic representation of initial moisture field in a mesoscale model is critical for the model simulation of thunderstorm development, and the formation of non-convective clouds as well as their effects on the surface energy budget. The impact will be investigated based on a real-data case study using the mesoscale atmospheric simulation system developed by Mesoscale Environmental Simulations Operations, Inc. The mesoscale atmospheric simulation system consists of objective analysis and initialization codes, and the coarse-mesh and fine-mesh dynamic prediction models. Both models are a three dimensional, primitive equation model containing the essential moist physics for simulating and forecasting mesoscale convective processes in the atmosphere. The modeling system is currently implemented at the Applied Meteorology Unit, Kennedy Space Center. Two procedures involving the synthetic relative humidity to define the model initial moisture fields are considered. It is proposed to perform several short-range (approximately 6 hours) comparative coarse-mesh simulation experiments with and without the synthetic data. They are aimed at revealing the model sensitivities should allow us both to refine the specification of the observational requirements, and to develop more accurate and efficient objective analysis schemes. The goal is to advance the MASS (Mesoscal Atmospheric Simulation System) modeling expertise so that the model output can provide reliable guidance for thunderstorm forecasting.

Chang, Chia-Bo↗

Simulation scaling studies of reactor core two-phase flow using direct numerical simulation

Tremendous growth in supercomputing power in recent years has resulted in the emergence of high-resolution flow analysis methods as an advanced research tool to evaluate single and two-phase flow behavior. In particular, unstructured mesh-based methods have been applied to analyze flows in complex reactor core geometries, including those of light water reactors (LWR). The finite-element based code, PHASTA, is utilized to perform large-scale simulations of two-phase bubbly flows in LWR geometries. Given the large computational cost of direct numerical simulation (DNS) coupled with interface tracking methods (ITM), typical domains encompass a portion of a single subchannel. In the presented research, the state-of-the-art analysis of turbulent two-phase flows in complex LWR subchannel geometries are demonstrated at both prototypical reactor parameters as well as scaled low pressure conditions. Three different cases are studied, a high-pressure simulation in prototypical reactor subchannel geometry, a low-pressure case in prototypical geometry and a final low-pressure case in a geometry scaled up to conserve the ratio between the bubble size and the domain pitch. Utilizing advanced statistical processing tools, these simulation conditions are compared to shed light on the relevancy of two-phase flow characteristics given the significant differences between LWR and low-pressure conditions. These findings can lead to the generation of useful guiding principles when researchers need to scale the two-phase flow behavior captured at low pressure and temperature conditions to those at reactor operating conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Regolith Simulant Report Cards: A User Friendly Simulant Certification System

As the nation prepares to return to the Moon, there is an increasing need for testing tools, instruments, and equipment in simulated lunar regolith to ensure successful operations during lunar missions. Because of this, it is important to use accurate lunar regolith simulants to prepare for these missions. To address this, we are building an objective certification and reporting system for regolith simulants to ensure suitable simulants are being used for the appropriate use cases. Previous systems like the Figures of Merit (FoM) [1] and the Lunar Regolith Simulant User’s Guide fit-to-use matrices [2] for lunar simulants were not particularly user friendly. Some of the Figures of Merit were not clearly defined and did not include recommended baseline values of lunar regolith properties to compare simulants to, and the fit-to-use matrices are out of date and include simulants that are no longer in production and do not include simulants that have been developed since.

S. R. Deitrick↗

OCEANFILMS (Organic Compounds from Ecosystems to Aerosols: Natural Films and Interfaces via Langmuir Molecular Surfactants) sea spray organic aerosol emissions – implementation in a global climate model and impacts on clouds

Sea spray aerosol is one of the major sources of atmospheric particulate matter globally. It has increasingly been recognized that organic matter derived from ocean biological precursors contributes significantly to the composition of submicron sea spray and may modify sea spray aerosol impacts on clouds and climate. This paper describes the implementation of the OCEANFILMS (Organic Compounds from Ecosystems to Aerosols: Natural Films and Interfaces via Langmuir Molecular Surfactants) parameterization for sea spray organic aerosol emissions in a global Earth system model, the Energy Exascale Earth System Model (E3SM). OCEANFILMS is a physically based model that links sea spray chemistry with ocean biogeochemistry using a Langmuir partitioning approach. We describe the implementation details of OCEANFILMS within E3SM, compare simulated aerosol fields with observations, and investigate impacts on simulated clouds and climate. Four sensitivity cases are tested, in which organic emissions either strictly add to or strictly replace sea salt emissions (in mass and number) and are either fully internally or fully externally mixed with sea salt. The simulation with internally mixed, added organics agrees reasonably well with observed seasonal cycles of organic matter in marine aerosol and has been selected as the default configuration of the E3SM. In this configuration, sea spray organic aerosol contributes an additional source of cloud condensation nuclei, adding up to 30 cm -3 to Southern Ocean boundary-layer cloud condensation nuclei concentrations (supersaturation = 0.1 %). The addition of this new aerosol source strengthens shortwave radiative cooling by clouds by -0.36 W m -2 in the global annual mean and contributes more than -3.5 W m -2 to summertime zonal mean cloud forcing in the Southern Ocean, with maximum zonal mean impacts of about -4 W m -2 around 50–60°S. This is consistent with a previous top-down, satellite-based empirical estimate of the radiative forcing by sea spray organic aerosol over the Southern Ocean. Through its mechanistic approach, OCEANFILMS offers a path towards improved understanding of the feedbacks between ocean biology, sea spray organic matter, and climate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Advanced Distributed Wind Turbine Controls Series: Part 2-Wind Energy in Isolated Grids – Microgrids, Infrastructure Resilience, and Advanced Controls Launchpad (MIRACL)

In an isolated grid, wind turbines are typically deployed to provide energy to maximize energy production, reduce diesel-fuel consumption, reduce carbon emissions, and reduce costs for energy and fuel transportation. However, in addition to solely providing energy to the power system, wind turbines contain rotating masses and inverter-based controls that can enable various reliability and resilience services through advance controls. As part of the Microgrids, Infrastructure Resilience, and Advanced Controls Launchpad (MIRACL) this paper demonstrates, through desktop simulations, advanced wind turbine controls that can be employed to support higher contributions of wind in isolated grids, and to demonstrate ways that wind can play a role in supporting stability of an isolated grid. This isolated grid used in these desktop simulations is comprised of a wind turbine (600 kW), solar PV (430 kW), battery energy storage system (1 MW/1MWh), a simulated diesel generator (2 MW) and various types of loads (critical, dynamic). We developed a model of the subsystems in MATLAB/Simulink and validated them with available data from real-world components on NREL's Flatirons Campus. These validated models are then configured for various case studies. We compare the output of the desktop simulation with a baseline case with the diesel generator. Active and reactive power control of the wind turbine can help improve frequency and voltage responses in the isolated grid, respectively. By utilizing a small integrated battery energy storage system in the DC-link of the wind turbine, we also demonstrate that wind turbines can help blackstart a critical load comparable to its rated power and support other renewables (e.g. solar PV) come online and pick up an additional load. This report illustrates some of these reliability and resilience services a wind turbine can provide in an isolated grid.

17 WIND ENERGY↗

Identifying Meteorological Drivers for Errors in Modeled Winds along the Northern California Coast

Abstract An accurate wind resource dataset is required for assessing the potential energy yield of floating offshore wind farms that are expected along the California outer continental shelf. The National Renewable Energy Laboratory has developed and disseminated an updated wind resource dataset offshore of California, using the Weather Research and Forecasting Model, referred to as the CA20 dataset. As compared to buoy lidar measurements that have become available recently, the CA20 dataset showed significant positive biases for 100-m wind speeds along Northern California wind energy lease areas. To investigate the meteorological drivers for the model errors, we first consider two 1-yr simulations run with two different planetary boundary layer (PBL) parameterizations: the Mellor–Yamada–Nakanishi–Niino (MYNN) PBL scheme (the chosen configuration in the CA20 dataset) and the Yonsei University PBL scheme (which significantly reduces the bias in modeled winds). By comparing the 1-yr simulations to the concurrent lidar buoy observations, we find that errors are larger with the MYNN PBL scheme in warm seasons. We then dive deeper into the analysis by running simulations for short-term (3-day) case studies to evaluate the sensitivity of initial/boundary condition forcings on model results. By analyzing the short-term simulations, we find that during synoptic-scale northerly flows driven by the North Pacific high and inland thermal low, a coastal warm bias in the MYNN simulation is mainly responsible for the modeled wind speed bias by altering the boundary layer thermodynamics. The results of our analysis will help guide the creation of an updated version of the CA20 dataset.

17 WIND ENERGY↗

Space-flight simulations of calcium metabolism using a mathematical model of calcium regulation

The results of a series of simulation studies of calcium matabolic changes which have been recorded during human exposure to bed rest and space flight are presented. Space flight and bed rest data demonstrate losses of total body calcium during exposure to hypogravic environments. These losses are evidenced by higher than normal rates of urine calcium excretion and by negative calcium balances. In addition, intestinal absorption rates and bone mineral content are assumed to decrease. The bed rest and space flight simulations were executed on a mathematical model of the calcium metabolic system. The purpose of the simulations is to theoretically test hypotheses and predict system responses which are occurring during given experimental stresses. In this case, hypogravity occurs through the comparison of simulation and experimental data and through the analysis of model structure and system responses. The model reliably simulates the responses of selected bed rest and space flight parameters. When experimental data are available, the simulated skeletal responses and regulatory factors involved in the responses agree with space flight data collected on rodents. In addition, areas within the model that need improvement are identified.

Brand, S. N.↗

Aircraft microwave observations and simulations of deep convection from 18 to 183 GHz. II - Model results

In this June 29, 1986 case study, a radiative transfer model is used to simulate the aircraft multichannel microwave brightness temperatures presented in the Adler et al. (1990) paper and to study the convective storm structure. Ground-based radar data are used to derive hydrometeor profiles of the storm, based on which the microwave upwelling brightness temperatures are calculated. Various vertical hydrometeor phase profiles and the Marshall and Palmer (M-P, 1948) and Sekhon and Srivastava (S-S, 1970) ice particle size distributions are experimented in the model. The results are compared with the aircraft radiometric data. The comparison reveals that the M-P distribution well represents the ice particle size distribution, especially in the upper tropospheric portion of the cloud; the S-S distribution appears to better simulate the ice particle size at the lower portion of the cloud, which has a greater effect on the low-frequency microwave upwelling brightness temperatures; and that, in deep convective regions, significant supercooled liquid water (about 0.5 g/cu m) may be present up to the -30 C layer, while in less convective areas, frozen hydrometeors are predominant above -10 C level.

Yeh, Hwa-Young M.↗

Transactive Energy Rationing in Islanded Electric Power Systems

Extreme weather events, many of which are climate change related, are occurring with increasing frequency and intensity and causing catastrophic outages. One of the major modern-day concerns of utilities is dealing with such extreme outages and consequently, its repercussions on the lives of people in society and social aspects. Traditionally, operators have been using rolling blackouts as a contingency plan to serve critical loads during such events when electricity supply is scarce. However, such blackouts practices are executed on a last-minute mandatory basis, depriving customers having low-capacity high priority loads (i.e. refrigeration, water, telecommunication, etc) that should also be serviced if at all possible. Historically, policymakers have often adopted quota-based regulatory actions or rationing for other commodities (such as gasoline, butter, sugar) to handle scarcity situations. Motivated from such quota based systems, our contribution presented in this work is an alternative transactive rationing mechanism that would provide some minimum level of service to all of the customers and serve the critical loads using market-based control during such scarcity-based contingencies. This is in contrast to the state-of-art TE mechanisms that allocate resources to customers solely based on their willingness-to-pay. The effectiveness of the proposed mechanism is demonstrated through simulation-based evaluation on two real-life use cases having feeder-level and microgrid-level configurations respectively. The simulation results clearly demonstrate the capability of the rationing scheme to serve customers’ high-priority loads through prolonged outages even during extreme scarcity scenarios.

42 ENGINEERING↗