Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “modeling simulation”

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,477 records · Page 82

Changing seasonality of moderate and extreme precipitation events in the Alps

The intensity of precipitation events is expected to increase in the future.The rate of increase depends on the strength or rarity of the events; very strong and rare events tend to follow the Clausius–Clapeyron relation,whereas weaker events or precipitation averages increase at a smaller rate than expected from the Clausius–Clapeyron relation. An often overlooked aspect is seasonal occurrence of such events, which might change in the future. To address the impact of seasonality, we use a large ensemble of regional and global climate model simulations, comprising tens of thousands of model years of daily temperature and precipitation for the past, present,and future. In order to make the data comparable, they are quantile mapped to observation-based time series representative of the Aare catchment in Switzerland. Model simulations show no increase in annual maximum 1-day precipitation events (Rx1day) over the last 400 years and an increase of 10 %–20 % until the end of the century for a strong (RCP8.5) forcing scenario. This fits with a Clausius–Clapeyron scaling of temperature at the event day, which increases less than annual mean temperature. An important reason for this is a shift in seasonality. Rx1day events become less frequent in late summer and more frequent in early summer and early autumn, when it is cooler. The seasonality shift is shown to be related to summer drying. Models with decreasing annual mean or summer mean precipitation show this behaviour more strongly. The highest Rx1day per decade, in contrast, shows no change in seasonality in the future. This discrepancy implies that decadal-scale extremes are thermodynamically limited; conditions conducive to strong events still occur during the hottest time of the year on a decadal scale. In contrast, Rx1day events are also limited by other factors. Conducive conditions are not reached every summer in the present, and even less so in the future. Finally, results suggest that changes in the seasonal cycle need to be accounted for when preparing for moderately extreme precipitation events and assessing their socio-economic impacts.

54 ENVIRONMENTAL SCIENCES↗

Final Report for DOE contract DE-FG02-98ER62609

Several publications document, analyze, and map the paleoclimatic data available for the last glacial maximum (LGM) 21,000 years ago. Data coverage is fairly extensive for North America, Eurasia, and Africa. The data show much lower temperatures for that time in earth history. These data were used to check climate model simulations for the LGM. Simulations forced by the CLIMAP SST for the LGM, where the west/east SST gradient across the Pacific is reduced compared to present, produce a reduction in the strength of the trade winds and a decrease in the west/east slope of the equatorial thermocline that is incompatible with thermocline depths newly inferred from foraminiferal assemblages. Stronger than-present trade winds, and a more realistic simulation of the thermocline slope, are produced when eastern Pacific SSTs are 2°C cooler than western Pacific SSTs. Sensitivity experiments for the model simulations of the LGM were also used to show which aspects of the climate boundary conditions produced various changes in model simulations for the LGM. An experiment which has the LGM ice sheets for the Northern Hemisphere but no Eurasian ice sheet shows that the downstream cooling in Asia is largely dependent on sea-ice growth in the North Atlantic, which is a product of the Laurentide ice sheet. Weaker Asian monsoons develop as a result of this downstream cooling, which is largest for low Laurentide ice sheets of large areal extent.

54 ENVIRONMENTAL SCIENCES↗

Modeling and simulation of consumer response to dynamic pricing.

Assessing the impacts of dynamic-pricing under the smart grid concept is becoming extremely important for deciding its full deployment. In this paper, we develop a model that represents the response of consumers to dynamic pricing. In the model, consumers use forecasted day-ahead prices to shift daily energy consumption from hours when the price is expected to be high to hours when the price is expected to be low while maintaining the total energy consumption as unchanged. We integrate the consumer response model into the Electricity Market Complex Adaptive System (EMCAS). EMCAS is an agent-based model that simulates restructured electricity markets. We explore the impacts of dynamic-pricing on price spikes, peak demand, consumer energy bills, power supplier profits, and congestion costs. A simulation of an 11-node test network that includes eight generation companies and five aggregated consumers is performed for a period of 1 month. In addition, we simulate the Korean power system.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Energy savings and life cycle cost analysis of advanced double skin facade system applied to old apartments in South Korea

South Korea is one of the most densely populated countries in the world, and the population density in urban areas is much higher among others. Apartments are the most common form of residential buildings due to their high population density. More than 60% of domestic residential buildings are apartments, with more than 10 million households. The high population density is a problem that has long plagued South Korea because more than 45% of apartments are old apartments that are more than 20 years old. Old apartment buildings have low thermal performance that results in lots of energy consumption. Balconies, which worked as a thermal buffer, are also being removed by residents to expand indoor space. According to the green remodeling project supported by the Korean government, the typical retrofitting method is replacing external windows with the high-efficiency window in old apartments in South Korea. Here, this paper proposes a retrofitting method using an extended double-skin façade (DSF) system that replaces exterior windows and acts as a thermal buffer. The simulation model is developed with the EnergyPlus simulation program to conduct energy performance analyses. The simulation model has been calibrated using actual experimental data collected between October 1, 2019, and January 15, 2020. Results show that the cooling, heating, and lighting energy can be reduced up to 44.1% by fully utilizing electric energy generated by solar PV panels integrated with the DSF system. The payback period is about 15 years based on the energy price provided by the utility company. Although the payback period is long, it has great potential of energy savings and CO 2 emission reductions. The DSF system should be considered as a way of renovation, considering other factors such as increased real estate values and energy cost increases in the future.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Roll-up of validation results to a target application.

Suites of experiments are preformed over a validation hierarchy to test computational simulation models for complex applications. Experiments within the hierarchy can be performed at different conditions and configurations than those for an intended application, with each experiment testing only part of the physics relevant for the application. The purpose of the present work is to develop methodology to roll-up validation results to an application, and to assess the impact the validation hierarchy design has on the roll-up results. The roll-up is accomplished through the development of a meta-model that relates validation measurements throughout a hierarchy to the desired response quantities for the target application. The meta-model is developed using the computation simulation models for the experiments and the application. The meta-model approach is applied to a series of example transport problems that represent complete and incomplete coverage of the physics of the target application by the validation experiments.

Hills, Richard Guy↗

Blueprints for Training Information Bottlenecks for Collider Analyses

Dimensionality reduction is a crucial aspect of data analysis in high energy physics, even if accompanied by information loss. Several methods, including histogram- and kernel-based analyses, are only computationally feasible for low-dimensional data. Furthermore, simulation models used in HEP can often only be validated for low-dimensional data. We provide several blueprints for using machine learning to create low-dimensional data representations (continuous event variables and discrete classification labels) for use in signal discovery and parameter estimation tasks. We also describe how to design the learned representation to facilitate a) searches with unknown model parameters and b) validation of simulation models in data control regions.

43 PARTICLE ACCELERATORS↗

A Eulerian-Lagrangian Model to Simulate Two-Phase/Particulate Flows

Figure 1 shows a snapshot of liquid fuel spray coming out of an injector nozzle in a realistic gas-turbine combustor. Here the spray atomization was simulated using a stochastic secondary breakup model (Apte et al. 2003a) with point-particle approximation for the droplets. Very close to the injector, it is observed that the spray density is large and the droplets cannot be treated as point-particles. The volume displaced by the liquid in this region is significant and can alter the gas-phase ow and spray evolution. In order to address this issue, one can compute the dense spray regime by an Eulerian-Lagrangian technique using advanced interface tracking/level-set methods (Sussman et al. 1994; Tryggvason et al. 2001; Herrmann 2003). This, however, is computationally intensive and may not be viable in realistic complex configurations. We therefore plan to develop a methodology based on Eulerian-Lagrangian technique which will allow us to capture the essential features of primary atomization using models to capture interactions between the fluid and droplets and which can be directly applied to the standard atomization models used in practice. The numerical scheme for unstructured grids developed by Mahesh et al. (2003) for incompressible flows is modified to take into account the droplet volume fraction. The numerical framework is directly applicable to realistic combustor geometries. Our main objectives in this work are: Develop a numerical formulation based on Eulerian-Lagrangian techniques with models for interaction terms between the fluid and particles to capture the Kelvin- Helmholtz type instabilities observed during primary atomization. Validate this technique for various two-phase and particulate flows. Assess its applicability to capture primary atomization of liquid jets in conjunction with secondary atomization models.

Apte, S. V.↗

Rapid Prototyping of an Aircraft Model in an Object-Oriented Simulation

A team was created to participate in the Mars Scout Opportunity. Trade studies determined that an aircraft provided the best opportunity to complete the science objectives of the team. A high fidelity six degree of freedom flight simulation was required to provide credible evidence that the aircraft design fulfilled mission objectives and to support the aircraft design process by providing performance evaluations. The team created the simulation using the Langley Standard Real-Time Simulation in C++ (LaSRS++) application framework. A rapid prototyping approach was necessary because the team had only three months to both develop the aircraft simulation model and evaluate aircraft performance as the design and mission parameters matured. The design of LaSRS++ enabled rapid-prototyping in several ways. First, the framework allowed component models to be designed, implemented, unit-tested, and integrated quickly. Next, the framework provides a highly reusable infrastructure that allowed developers to maximize code reuse while concentrating on aircraft and mission specific features. Finally, the framework reduces risk by providing reusable components that allow developers to build a quality product with a compressed testing cycle that relies heavily on unit testing of new components.

Kenney, P. Sean↗

Ground Contact Modeling for the Morpheus Test Vehicle Simulation

The Morpheus vertical test vehicle is an autonomous robotic lander being developed at Johnson Space Center (JSC) to test hazard detection technology. Because the initial ground contact simulation model was not very realistic, it was decided to improve the model without making it too computationally expensive. The first development cycle added capability to define vehicle attachment points (AP) and to keep track of their states in the lander reference frame (LFRAME). These states are used with a spring damper model to compute an AP contact force. The lateral force is then overwritten, if necessary, by the Coulomb static or kinetic friction force. The second development cycle added capability to use the PolySurface class as the contact surface. The class can load CAD data in STL (Stereo Lithography) format, and use the data to compute line of sight (LOS) intercepts. A polygon frame (PFRAME) is computed from the facet intercept normal and used to convert the AP state to PFRAME. Three flat plane tests validate the transitions from kinetic to static, static to kinetic, and vertical impact. The hazardous terrain test will be used to test for visual reasonableness. The improved model is numerically inexpensive, robust, and produces results that are reasonable.

Cordova, Luis↗

Ground Contact Modeling for the Morpheus Test Vehicle Simulation

The Morpheus vertical test vehicle is an autonomous robotic lander being developed at Johnson Space Center (JSC) to test hazard detection technology. Because the initial ground contact simulation model was not very realistic, it was decided to improve the model without making it too computationally expensive. The first development cycle added capability to define vehicle attachment points (AP) and to keep track of their states in the lander reference frame (LFRAME). These states are used with a spring damper model to compute an AP contact force. The lateral force is then overwritten, if necessary, by the Coulomb static or kinetic friction force. The second development cycle added capability to use the PolySurface class as the contact surface. The class can load CAD data in STL (Stereo Lithography) format, and use the data to compute line of sight (LOS) intercepts. A polygon frame (PFRAME) is computed from the facet intercept normal and used to convert the AP state to PFRAME. Three flat plane tests validate the transitions from kinetic to static, static to kinetic, and vertical impact. The hazardous terrain test will be used to test for visual reasonableness. The improved model is numerically inexpensive, robust, and produces results that are reasonable.

Cordova, Luis↗

Changes in Anthropogenic Precursor Emissions Drive Shifts in the Ozone Seasonal Cycle Throughout the Northern Midlatitude Troposphere

Simulations by six Coupled Model Intercomparison Project Phase 6 (CMIP6) Earth system models indicate that the seasonal cycle of baseline tropospheric ozone at northern midlatitudes has been shifting since the mid-20th century. Beginning in ∼ 1940, the magnitude of the seasonal cycle increased by ∼10 ppb (measured from seasonal minimum to maximum), and the seasonal maximum shifted to later in the year by about 3 weeks. This shift maximized in the mid-1980s, followed by a reversal – the seasonal cycle decreased in amplitude and the maximum shifted back to earlier in the year. Similar changes are seen in measurements collected from the 1970s to the present. The timing of the seasonal cycle changes is generally concurrent with the rise and fall of anthropogenic emissions that followed industrialization and the subsequent implementation of air quality emission controls. A quantitative comparison of the temporal changes in the ozone seasonal cycle at sites in both Europe and North America with the temporal changes in ozone precursor emissions across the northern midlatitudes found a high degree of similarity between these two temporal patterns. We hypothesize that changing precursor emissions are responsible for the shift in the ozone seasonal cycle; this is supported by the absence of such seasonal shifts in southern midlatitudes where anthropogenic emissions are much smaller. We also suggest a mechanism by which changing emissions drive the changing seasonal cycle: increasing emissions of NOx allow summertime photochemical production of ozone to become more important than ozone transported from the stratosphere, and increasing volatile organic compounds (VOCs) lead to progressively greater photochemical ozone production in the summer months, thereby increasing the amplitude of the seasonal ozone cycle. Decreasing emissions of both precursor classes then reverse these changes. The quantitative parameter values that characterize the seasonal shifts provide useful benchmarks for evaluating model simulations, both against observations and between models.

Ozone↗

Interannual to Decadal Variability in Subarctic Atlantic Ocean Primary Productivity Historical Variability, Dynamics, and Response to Warming

Phytoplankton primary productivity in the Subarctic Atlantic Ocean sustains valuable fisheries and exports anthropogenic atmospheric carbon dioxide, but multiple Earth system models show that this productivity declines substantially in 21st century global warming scenarios. Prior studies have investigated the spatial and temporal dynamics of the seasonal cycle of the phytoplankton, which exhibit a prominent but spatially variable spring bloom that is caused by the closely related seasonal cycles in solar radiation and vertical mixing that modulate both phytoplankton growth and loss rates. However, recent work indicates that multi-decade declines in this spring bloom with global warming in multiple Earth system models are driven by reductions in circulation of nutrients from the subtropical and tropical thermocline, rather than local changes in vertical mixing or shortwave radiation in the Subarctic. It is unclear if current models of how phytoplankton productivity varies seasonally or over a century in global warming can be interpolated to explain the interannual to decadal variability. In this poster, I will quantify the interannual to decadal variability of phytoplankton chlorophyll using historical satellite ocean color observations and model simulations with the Community Earth System Model. I will then use the observations and model simulations to elucidate the mechanisms of interannual to decadal variability and how that variability is restructured in 21st century global warming scenarios.

nterannual↗

Impacts of aerosols on seasonal precipitation and snowpack in California based on convection-permitting WRF-Chem simulations

Here, a version of the WRF-Chem model with fully coupled aerosol–meteorology–snowpack is employed to investigate the impacts of various aerosol sources on precipitation and snowpack in California. In particular, the impacts of locally emitted anthropogenic and dust aerosols, and aerosols transported from outside California are studied. We differentiate three pathways of aerosol effects: aerosol–radiation interaction (ARI), aerosol–snow interaction (ASI), and aerosol–cloud interaction (ACI). The convection-permitting model simulations show that precipitation, snow water equivalent (SWE), and surface air temperature averaged over the whole domain (34–42° N, 117–124° W, not including ocean points) are reduced when aerosols are included, therefore reducing large biases in these variables due to the absence of aerosol effects in the model. Aerosols affect California water resources through the warming of mountaintops and the reduction of precipitation; however, different aerosol sources play different roles in changing surface temperature, precipitation, and snowpack in California by means of various weights of the three pathways. ARI by all aerosols mainly cools the surface, leading to slightly increased SWE over the mountains. Locally emitted dust aerosols warm the surface of mountaintops through ASI, in which the reduced snow albedo associated with dusty snow leads to more surface absorption of solar radiation and reduced SWE. Transported aerosols and local anthropogenic aerosols play a dominant role in increasing nonprecipitating clouds but reducing precipitation through ACI, leading to reduced SWE and runoff on the Sierra Nevada, as well as the warming of mountaintops associated with decreased SWE and hence lower surface albedo. The average changes in surface temperature from October 2012 to June 2013 are about –0.19 and 0.22 K for the whole domain and over mountaintops, respectively. Overall, the averaged reduction during October to June is about 7 % for precipitation, 3 % for SWE, and 7 % for surface runoff for the whole domain, while the corresponding numbers are 12, 10, and 10 % for the mountaintops. The reduction in SWE is more significant in a dry year, with 9 % for the whole domain and 16 % for the mountaintops. The maximum reduction of ~20 % in precipitation occurs in May and is associated with the maximum aerosol loading, leading to the largest decrease in SWE and surface runoff over that period. It is also found that dust aerosols can cause early snowmelt on the mountaintops and reduced surface runoff after April.

54 ENVIRONMENTAL SCIENCES↗

Modeling and simulation to investigate the electrification potential of medium- and heavy-duty vehicle fleets

This project involves developing and integrating new modeling tools to simulate the dynamics of electric medium- and heavy-duty fleet vehicle adoption. A technical and economic modeling tool, combining a data-driven hardware cost model with a cost-optimal charging strategy microsimulation, enables tailored analysis of the costs and benefits of electrifying individual fleets. Next, a novel text synthesis process, applied to a curated corpus of literature, quantifies trade-offs between technical, economic, and other factors in the fleet vehicle procurement decision. The outcomes of these tasks combine with knowledge from recent literature on fleet decision processes to specify the vehicle procurement model used by fleets in an agent-based model of the medium- and heavy-duty electric vehicle market. This model embodies an especially disaggregated approach to adoption modeling, internalizing factors and dynamics that conventional adoption models externalize. In particular, explicitly modeling the formation and diffusion of opinions among agents enables experiments that conventional models cannot support. Demonstrations show, for example, that increasing the extent of interactions between populations with different proclivities to electric vehicles has an asymmetrical outcome. High-proclivity electric vehicle adoption is generally unaffected as interactions increase, but low-proclivity adoption is accelerated. By representing individual fleets' requirements and costs at a high level of detail, incorporating an adoption decision model informed by a wide body of empirical research, and broadening the array of variables and dynamics available for experimentation, this integrated model offers a new way to understand the urgent challenge of eliminating emissions from the most emissions-intensive transportation sectors.

Trinko, David A.↗

Compressible homogeneous shear: Simulation and modeling

Compressibility effects were studied on turbulence by direct numerical simulation of homogeneous shear flow. A primary observation is that the growth of the turbulent kinetic energy decreases with increasing turbulent Mach number. The sinks provided by compressible dissipation and the pressure dilatation, along with reduced Reynolds shear stress, are shown to contribute to the reduced growth of kinetic energy. Models are proposed for these dilatational terms and verified by direct comparison with the simulations. The differences between the incompressible and compressible fields are brought out by the examination of spectra, statistical moments, and structure of the rate of strain tensor.

Sarkar, S.↗