Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Building Energy Modeling”

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,801 records · Page 100

Renewable Thermal Hybridization Framework for Industrial Process Heat Applications

Solar industrial process heat (SIPH) technologies, such as concentrating solar power collectors, could economically replace the steam or heat needs at many industrial sites by providing high-temperature heat transfer fluids (HTFs) such as pressurized water, synthetic-oil, or direct steam. Renewable thermal energy systems (RTES) could be hybridized with different renewable options e.g., flat plate collectors with parabolic trough collectors, or combined with existing heat supplies (e.g., fossil fuels), to give options for targeted SIPH applications, industrial decarbonization and the reduction of fuel consumption. Hybrid solutions and thermal energy storage will be important for the dispatch of heat at optimal times needed by the demand side of the buildings and industrial applications. At present, there is no integrated modeling tool for hybrid RTES, and this paper highlights the development of a renewable thermal hybridization framework for IPH use that is built from existing tools like System Advisor Model. The long-term vision for the framework (through significant further research) is to develop a coupled hybrid energy generation and cost analysis tool, where the tool could help the user in determining the most suitable and cost-effective technologies for their applications. Ongoing work will look to add costs for RTES options and further refinement on the selection of suitable technologies. This future tool could calculate the levelized cost of heat of various RTES hybrid options, by taking the user's solar resource, fuel costs, industrial heat demand profile, available land, and other factors into account to determine the applicability into their process.

concentrating solar power↗

Modeling the First Hydrogen Direct Reduction Pilot Reactor for Ironmaking in the USA Using Finite Element Analysis and Its Validation Using Pilot Plant Trial Data

Direct reduction of hematite pellets with hydrogen (H 2 ) was used to produce directly reduced iron (DRI) in a pilot scale reactor at a pellet feed rate of 21.4 kg/h. At a steady state, operational parameters of the pilot plant (gas recycling rate and inlet temperature) along with key reactor output parameters, the pellet metallization, and the internal temperature profile of the reactor were reported for two scenarios with high recycle and low recycle rate of H 2 . Scenario 1, with a high recycle rate of 400 L/min H 2 along with external heating of 870 °C, gave an average metallization of 91.8%, while Scenario 2, with low recycle rate of 100 L/min H 2 and external heating of 850 °C gave a metallization of 67.8% due to the higher moles of H 2 available for reduction and the external energy required for the endothermic reduction reaction in Scenario 1 as compared with Scenario 2. Finite element analysis was used to build a model of the shaft reactor, which was validated against the metallization and internal temperature profile data. The average metallization values predicted by the model were very close to the metallization values obtained from the pilot plant samples, with 90.9% average metallization for Scenario 1 and 65.6% average metallization for Scenario 2. The internal temperature profiles in the lower region of the reactor obtained from the model were very close to these pilot plant data, with a maximum difference of 52.7 °C and 67.6 °C for Scenarios 1 and 2, respectively. The pilot plant reactor model was used extensively in the commissioning of the pilot plant and to predict the startup outcomes for a given set of operating parameters.

08 HYDROGEN↗

Occupancy-Based Controls for an All-Electric Residential Community in a Cold Climate

In residential buildings, rapid improvements in sensors, communication, and information technology have enabled occupancy-based building controls. These controls utilize occupancy information and modify the operation of the heating, ventilation, and air-conditioning (HVAC) system to minimize excess HVAC energy use, especially when the building is unoccupied. This reduces the total building energy consumption and utility bills while maintaining thermal comfort. In this paper, we present two novel occupancy-driven controls - reactive control and predictive control - and compare their performance. We model an all-electric residential community based on a 27- home community in Basalt, Colorado, in the United States. We simulated various scenarios, considering different temperature setback and control algorithms, to analyze the community-scale impact of these occupancy-based controls. The results show that total HVAC energy savings in a building ranges from 1%- 20% compared to the baseline scenario without occupancy-based controls. The energy-saving potential is highly correlated with the occupancy pattern and temperature setback in the building.

ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATION↗

Occupancy-Based Controls for an All-Electric Residential Community in a Cold Climate: Preprint

In residential buildings, rapid improvements in sensors, communication, and information technology have enabled occupancy-based building controls. These controls utilize occupancy information and modify the operation of the heating, ventilation, and air-conditioning (HVAC) system to minimize excess HVAC energy use, especially when the building is unoccupied. This reduces the total building energy consumption and utility bills while maintaining thermal comfort. In this paper, we present two novel occupancy-driven controls - reactive control and predictive control - and compare their performance. We model an all-electric residential community based on a 27- home community in Basalt, Colorado, in the United States. We simulated various scenarios, considering different temperature setback and control algorithms, to analyze the community-scale impact of these occupancy-based controls. The results show that total HVAC energy savings in a building ranges from 1%- 20% compared to the baseline scenario without occupancy-based controls. The energy-saving potential is highly correlated with the occupancy pattern and temperature setback in the building.

ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATION↗

The Distribution System Operator with Transactive (DSO+T) Study

The Distribution System Operator with Transactive (DSO+T) study investigates the engineering and economic performance of a transactive energy retail market coordinating a high penetration of customer-side flexible energy assets. The study seeks to answer whether such an implementation is cost effective for customers, recovers sufficient revenue for DSOs, and is equally applicable and beneficial to a range of flexible asset types, renewable generation scenarios, and market assumptions. Using a highly interdisciplinary co-simulation and valuation framework, this assessment encompasses the entire electrical delivery system from bulk system generation and transmission, through the distribution system, to the modeling of individual customer buildings and flexible assets (including heating, ventilation, and air conditioning [HVAC] units, water heaters, batteries, and electric vehicles). The study exercises a transactive energy retail market coordination scheme designed to integrate with an existing day-ahead and real-time competitive wholesale electricity market. Software decision-making agents are designed for the retail market operator as well as various price-responsive flexible assets. The engineering and economic performance of the transactive energy scheme is studied for two separate flexible asset deployments: flexible loads (HVAC units and residential water heaters) and behind-the-meter batteries. The results of each transactive case are compared to a business-as-usual case. These cases are subject to two different renewable generation scenarios, a moderate renewable generation scenario, representative of current levels of renewable generation deployment, and a future high renewables scenario, including the increased deployment of rooftop solar photovoltaic and electric vehicles. The transactive coordination scheme is shown to produce effective and stable control and decrease peak loads 9–15%. The resulting annual demand flexibility provides net economic savings of $3.3–5.0B per year for a region the size of Texas. Detailed analysis shows that net benefits were seen for a range of distribution system operator, customer, and flexible asset types. Both participating customer (with transactive flexible assets) and nonparticipating customers (with nonflexible assets) see reductions in annual utility bills and net annual energy expenses in the range of 10–16%.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Distribution System Operator with Transactive (DSO+T) Study: Volume 1 (Main Report)

The Distribution System Operator with Transactive (DSO+T) study investigates the engineering and economic performance of a transactive energy retail market coordinating a high penetration of customer-side flexible energy assets. The study seeks to answer whether such an implementation is cost effective for customers, recovers sufficient revenue for DSOs, and is equally applicable and beneficial to a range of flexible asset types, renewable generation scenarios, and market assumptions. Using a highly interdisciplinary co-simulation and valuation framework, this assessment encompasses the entire electrical delivery system from bulk system generation and transmission, through the distribution system, to the modeling of individual customer buildings and flexible assets (including heating, ventilation, and air conditioning [HVAC] units, water heaters, batteries, and electric vehicles). The study exercises a transactive energy retail market coordination scheme designed to integrate with an existing day-ahead and real-time competitive wholesale electricity market. Software decision-making agents are designed for the retail market operator as well as various price-responsive flexible assets. The engineering and economic performance of the transactive energy scheme is studied for two separate flexible asset deployments: flexible loads (HVAC units and residential water heaters) and behind-the-meter batteries. The results of each transactive case are compared to a business-as-usual case. These cases are subject to two different renewable generation scenarios, a moderate renewable generation scenario, representative of current levels of renewable generation deployment, and a future high renewables scenario, including the increased deployment of rooftop solar photovoltaic and electric vehicles. The transactive coordination scheme is shown to produce effective and stable control and decrease peak loads 9–15%. The resulting annual demand flexibility provides net economic savings of $3.3–5.0B per year for a region the size of Texas. Detailed analysis shows that net benefits were seen for a range of distribution system operator, customer, and flexible asset types. Both participating customer (with transactive flexible assets) and nonparticipating customers (with nonflexible assets) see reductions in annual utility bills and net annual energy expenses in the range of 10–16%.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Building a Cyberinfrastructure for Environmental System Science: Modeling Frameworks, Data Management, and Scientific Workflows, Workshop Report, Potomac, Maryland, April 30-May 1, 2015

The Environmental System Science (ESS) activity within the Department of Energy’s (DOE) Office of Biological and Environmental Research (BER) seeks to advance a robust, predictive understanding of terrestrial environments, extending from bedrock to the top of the vegetative canopy and from molecular to global scales, through an iterative cycle of model-driven experimentation and observation dubbed MODEX. Considerable progress has been made toward achieving this overarching goal, but widely recognized is the fragmentation across projects and disciplines of the relevant modeling and simulation capabilities, observational and experimental data, analysis algorithms, and workflow tools. This fragmentation creates significant challenges for studying impacts and feedbacks in these complex multiscale systems. These challenges are further exacerbated by ongoing disruptive changes in high-performance computational architectures and the exponential growth in the types and volume of data that render obsolete the conventional approaches to software development and data management. Overcoming these challenges will require the development of a BER Climate and Environmental Sciences Division (CESD)–wide enabling cyberinfrastructure to support data management, cross-domain modeling, data analysis, and collaborative research. To explore the potential for working groups to initiate and guide a more integrated and community-based cyberinfrastructure, BER held the ESS Workshop on Model-Data Integration: Modeling Frameworks, Data Management, and Scientific Workflows on April 30-May 1, 2015, following the ESS Principal Investigator Meeting in Potomac, Maryland. Participants included model developers, software engineers, and data management specialists from eight national laboratories, which represented a wide range of projects and programs from CESD as well as cross-cutting projects from DOE’s Office of Advanced Scientific Computing Research (ASCR; both BER and ASCR are operated from within DOE’s Office of Science). A series of plenary talks provided background information and clarified the workshop’s three objectives: (1) develop requirements for this community-based cyberinfrastructure to ensure enhanced scientific productivity of the community as a whole; (2) identify challenges associated with developing this new cyberinfrastructure using a phased approach guided by project-driven use cases; and, given these requirements and challenges, (3) chart a path forward for ESS working groups to lead the phased development of the new cyberinfrastructure.

54 ENVIRONMENTAL SCIENCES↗

Time-of-Use and Demand Charge Battery Controller Using Stochastic Model Predictive Control: Preprint

Stationary batteries in residential and commercial buildings are often used to smooth customer load profiles and to lower customer electricity bills. Controllers for these battery systems should account for customer energy consumption, rate structures, and high internal battery temperatures, which can lead to reduced performance over the battery lifetime. It is important to consider the uncertainty in forecasting energy consumption and temperature, especially for customers with highly variable and uncertain loads. We propose a novel battery controller using stochastic model predictive control that accounts for these uncertainties and can handle complex rate structures, including demand charges. We show that the controller performs better than standard model predictive control when there is significant uncertainty in the forecast. We also show improvements in the performance with more accurate forecasts and with a more aggressive control strategy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Grid-Interactive Multi-Zone Building Control Using Reinforcement Learning with Global-Local Policy Search

In this paper, we develop a grid-interactive multi-zone building controller based on a deep reinforcement learning (RL) approach. The controller is designed to facilitate building operation during normal conditions and demand response events, while ensuring occupants comfort and energy efficiency. We leverage a continuous action space RL formulation, and devise a two-stage global-local RL training framework. In the first stage, a global fast policy search is performed using a gradient-free RL algorithm. In the second stage, a local fine-tuning is conducted using a policy gradient method. In contrast to the state-of-the-art model predictive control (MPC) approach, the proposed RL controller does not require complex computation during real-time operation and can adapt to nonlinear building models. We illustrate the controller performance numerically using a five-zone commercial building.

30 DIRECT ENERGY CONVERSION↗

Adsorption of CO on Low-Energy, Low-Symmetry Pt Nanoparticles: Energy Decomposition Analysis and Prediction via Machine-Learning Models

We present a systematic analysis of CO adsorption on Pt nanoclusters in the 0.2–1.5 nm size range with the aim of unraveling size-dependent trends and developing predictive models for site-specific adsorption behavior. Using an empirical-potential-based genetic algorithm and density functional theory (DFT) modeling, we show that there exists a size window (40–70 atoms) over which Pt nanoclusters bind CO weakly, the binding energies being comparable to those on (111) or (100) facets. The size-dependent adsorption energy trends are, however, distinctly nonmonotonic and are not readily captured using traditional descriptors such as d-band energies or (generalized) coordination numbers of the Pt binding sites. Instead, by applying machine-learning algorithms, we show that multiple descriptors, broadly categorized as structural and electronic descriptors, are essential for qualitatively capturing the CO adsorption trends. Nevertheless, attaining quantitative accuracy requires further refinement, and we propose the use of an additional descriptors—the fully frozen adsorption energy—that is a computationally inexpensive probe of CO–Pt bond formation. With these three categories of descriptors, we achieve an absolute mean error in CO adsorption energy prediction of 0.12 eV, which is similar to the underlying error of DFT adsorption calculations. Our approach allows for building quantitatively predictive models of site-specific adsorbate binding on realistic, low-symmetry nanostructures, which is an important step in modeling reaction networks as well as for rational catalyst design in general.

36 MATERIALS SCIENCE↗

National Modeling of Geothermal District Energy Systems with Ambient-Temperature Loops Using dGeo

Geothermal district energy systems (DES) with ambient-temperature loops, also known as thermal energy networks, are one option for decarbonizing space heating and cooling loads. Geothermal fifth-generation DES include an "ambient" temperature thermal loop that connects heat pumps at each building with thermal balancing sources such as geothermal borehole fields. Heating and cooling are provided via a water-source heat pump at each end-user. This project seeks to analyze the nationwide potential for ambient-temperature loop districts by creating a new module within the Distributed Geothermal Market Demand Model (dGeo). dGeo is an agent-based modeling tool for distributed geothermal resources; it can investigate potential on a nationwide or statewide scale using geospatial data for all 50 states and thermal demands for existing buildings. This process allows for high-level estimates of technical and economic potential for ambient-temperature loop districts across the United States. A lookup table was created using GHEDesigner to size borehole fields for different thermal loads and ground conditions experienced across the country. A cost and financing structure, along with incentives, were applied. Cost estimates include costs for the distribution network, borehole field installation and operation, and circulation pump operation, while savings are calculated based on energy bills for building owners (agents). This newly developed module can be used for assessing which areas of the country have the highest potential for agent benefits from ambient-temperature loop installation and assess the impact of future cost and price scenarios. Initial results for statewide analysis (for Vermont) and nationwide (for United States) are provided. Future work includes expanding the module to consider mixed residential and commercial districts as well as evaluating multiple cost scenarios.

ambient-temperature loop↗

Recovery Act: Advanced Load Identification and Management for Buildings

In response to the U.S. Department of Energy (DoE)’s goal of achieving market ready, net-zero energy residential and commercial buildings by 2020 and 2025, Eaton partnered with the Department of Energy’s National Renewable Energy Laboratory (NREL) and Georgia Institute of Technology to develop an intelligent load identification and management technology enabled by a novel “smart power strip” to provide critical intelligence and information to improve the capability and functionality of building load analysis and building power management systems. Buildings account for 41% of the energy consumption in the United States, significantly more than either transportation or industrial. Within the building sector, plug loads account for a significant portion of energy consumption. Plug load consumes 15-20% of building energy on average. As building managers implement aggressive energy conservation measures, the proportion of plug load energy can increase to as much as 50% of building energy leaving plug loads as the largest remaining single source of energy consumption. This project focused on addressing plug-in load control and management to further improve building energy efficiency accomplished through effective load identification. The execution of the project falls into the following three major aspects; An intelligent load modeling, identification and prediction technology was developed to automatically determine the type, energy consumption, power quality, operation status and performance status of plug-in loads, using electric waveforms at a power outlet level. This project demonstrated the effectiveness of the developed technology through a large set of plug-in loads measurements and testing; A novel “Smart Power Strip (SPS) / Receptacle” prototype was developed to act as a vehicle to demonstrate the feasibility of load identification technology as a low-cost, embedded solution; and Market environment for plug-in load control and management solutions, in particular, advanced power strips (APSs) was studied. The project evaluated the market potential for Smart Power Strips (SPSs) with load identification and the likely impact of a load identification feature on APS adoption and effectiveness. The project also identified other success factors required for widespread APS adoption and market acceptance. Even though the developed technology is applicable for both residential and commercial buildings, this project is focused on effective plug-in load control and management for commercial buildings, accomplished through effective load identification. The project has completed Smart Receptacle (SR) prototype development with integration of Load ID, Control/Management, WiFi communication, and Web Service. Twenty SR units were built, tested, and demonstrated in the Eaton lab; eight SR units were tested in the National Renewable Energy Lab (NREL) for one-month of field testing. Load ID algorithm testing for extended load sets was conducted within the Eaton facility and at local university campuses. This report is to summarize the major achievements, activities, and outcomes under the execution of the project.

24 POWER TRANSMISSION AND DISTRIBUTION↗

LA100 Equity Strategies. Chapter 8: Equitable Rooftop Solar Access and Benefits

The LA100 Equity Strategies project integrates community guidance with robust research, modeling, and analysis to identify strategy options that can increase equitable outcomes in Los Angeles' clean energy transition. This chapter focuses on analysis of customer-sited rooftop solar and storage as a means to reduce electricity bills for low- and moderate-income (LMI) households, multifamily building residents, and renters, who traditionally lack access to bill savings from rooftop solar. Specifically, NREL modeled customer-sited solar and storage adoption using the Distributed Generation Market Demand (dGen™)1 model through 2035 and developed scenarios to identify programs or policies that could support equitable access to bill savings from rooftop solar or solar-plus-storage. Scenarios tested include a direct-install program for LMI customers, net metering for LMI customers, and equitable distribution of benefits from installing solar between owners and renters of renter-occupied buildings. Research was guided by input from the community engagement process, and equity strategies are presented in alignment with that guidance.

14 SOLAR ENERGY↗

Grid-Interactive Multi-Zone Building Control Using Reinforcement Learning with Global-Local Policy Search: Preprint

In this paper, we develop a grid-interactive multi-zone building controller based on a deep reinforcement learning (RL) approach. The controller is designed to facilitate building operation during normal conditions and demand response events, while ensuring occupants comfort and energy efficiency. We leverage a continuous action space RL formulation, and devise a two-stage global-local RL training framework. In the first stage, a global fast policy search is performed using a gradient-free RL algorithm. In the second stage, a local fine-tuning is conducted using a policy gradient method. In contrast to the state-of-the-art model predictive control (MPC) approach, the proposed RL controller does not require complex computation during real-time operation and can adapt to non-linear building models. We illustrate the controller performance numerically using a five-zone commercial building.

30 DIRECT ENERGY CONVERSION↗

Improving operational flexibility of integrated energy system with uncertain renewable generations considering thermal inertia of buildings

Insufficient flexibility in system operation caused by traditional "heat-set" operating modes of combined heat and power (CHP) units in winter heating periods is a key issue that limits renewable energy consumption. In order to reduce the curtailment of renewable energy resources through improving the operational flexibility, a novel optimal scheduling model based on chance-constrained programming (CCP), aiming at minimizing the lowest generation cost, is proposed for a small-scale integrated energy system (IES) with CHP units, thermal power units, renewable generations and representative auxiliary equipments. In this model, due to the uncertainties of renewable generations including wind turbines and photovoltaic units, the probabilistic spinning reserves are supplied in the form of chance-constrained; from the perspective of user experience, a heating load model is built with consideration of heat comfort and inertia in buildings. To solve the model, a solution approach based on sequence operation theory (SOT) is developed, where the original CCP-based scheduling model is tackled into a solvable mixed-integer linear programming (MILP) formulation by converting a chance constraint into its deterministic equivalence class, and thereby is solved via the CPLEX solver. We report the simulation results on the modified IEEE 30-bus system demonstrate that the presented method manages to improve operational flexibility of the IES with uncertain renewable generations by comprehensively leveraging thermal inertia of buildings and different kinds of auxiliary equipments, which provides a fundamental way for promoting renewable energy consumption.

30 DIRECT ENERGY CONVERSION↗

A simulation-based assessment of technologies to reduce heat emissions from buildings

Heat emissions from buildings are part of anthropogenic heat leading to urban overheating. This paper aims to assess how technologies (i.e., energy conservation measures - ECMs), used to decrease energy use, may also reduce heat emissions from buildings. This study employs the physics-based engine EnergyPlus to simulate the main components of heat emissions from buildings to ambient air: envelope, zone, and systems. Hourly simulations are run for IECC single- and multi-family reference models with three representative climates: Miami, Baltimore, and Chicago. The results show that the performance of ECMs varies among weather, seasons, and residential typologies. Particularly, some ECMs (i.e., cool coatings, heat pumps, additional insulation, energy-awareness occupants) show a strong decrease in heat emissions, yet they are not always correlated with proportional decreases in energy use. When all ECMs are combined, the reductions are larger on heat emissions (89%) than on site energy (65%) from the base cases. During summer in Miami, the combination of ECMs shows a decrease in heat emissions from the building surface component of 80% during daytime, 92% for the HVAC component and a counterbalanced increase in the zone component of 88%, bringing to a daily decrease in total heat emissions. The main contributions of this study are quantifying how typical ECMs influence residential building heat emissions using EnergyPlus simulations and informing urban planners and stakeholders on prioritizing measures for mitigating urban overheating problems.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Assessing the Energy Resilience of Office Buildings: Development and Testing of a Simplified Metric for Real Estate Stakeholders

Increasing concern over higher frequency extreme weather events is driving a push towards a more resilient built environment. In recent years there has been growing interest in understanding how to evaluate, measure, and improve building energy resilience, i.e., the ability of a building to provide energy-related services in the event of a local or regional power outage. In addition to human health and safety, many stakeholders are keenly interested in the ability of a building to allow continuity of operations and minimize business disruption. Office buildings are subject to significant economic losses when building operations are disrupted due to a power outage. We propose “occupant hours lost” (OHL) as a means to measure the business productivity lost as the result of a power outage in office buildings. OHL is determined based on indoor conditions in each space for each hour during a power outage, and then aggregated spatially and temporally to determine the whole building OHL. We used quasi-Monte Carlo parametric energy simulations to demonstrate how the OHL metric varies due to different building characteristics across different climate zones and seasons. The simulation dataset was then used to develop simple regression models for assessing the impact of ten key building characteristics on OHL. The most impactful were window-to-wall ratio and window characteristics. The regression models show promise as a simple means to assess and screen for resilience using basic building characteristics, especially for non-critical facilities where it may not be viable to conduct detailed engineering analysis.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Building blocks for a digital twin of additive manufacturing

Properties and serviceability of additively manufactured components are affected by their geometry, microstructure and defects. These important attributes are now optimized by trial and error because the essential process variables cannot currently be selected from scientific principles. A recourse is to build and rigorously validate a digital twin of the additive manufacturing process that can provide accurate predictions of the spatial and temporal variations of metallurgical parameters that affect the structure and properties of components. Key building blocks of a computationally efficient first-generation digital twin of laser-based directed energy deposition additive manufacturing utilize a transient, three-dimensional model that calculates temperature and velocity fields, cooling rates, solidification parameters and deposit geometry. The measured profiles of stainless steel 316L and Alloy 800H deposits as well as the secondary dendrite arm spacing (SDAS) and Vickers hardness measurements are used to validate the proposed digital twin. The predicted cooling rates, temperature gradients, solidification rates, SDAS and micro-hardness values are shown to be more accurate than those obtained from a commonly used heat conduction calculation. Finally, these metallurgical building blocks serve as a phenomenological framework for the development of a digital twin that will make the expanding knowledge base of additive manufacturing usable in a practical way for all scientists and engineers.

36 MATERIALS SCIENCE↗

Refine your search to explore more results.