Search NASA⌕ Search

SEARCH · Search NASA

Results for “loading simulations”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A combined experimental and crystal plasticity study of grain size effects in magnesium alloys

This work presents a method to incorporate the micro Hall-Petch equation into the crystal plasticity finite element (CPFE) framework accounting for the microstructural features to understand the coupling between grain size, texture, and loading direction in magnesium alloys. The effect of grain size and texture is accounted for by modifying the slip resistances of individual basal and prismatic systems based on the micro Hall-Petch equation. The modification based on the micro Hall-Petch equation endows every slip system at each microstructural point with a slip system-level grain size and maximum compatibility factor, which are in turn used to modify the slip resistance. While the slip-system level grain size is a measure of the grain size, the maximum compatibility factor encodes the effect of the grain boundary on the slip system resistance modification and is computed based on the Luster-Morris factor. The model is calibrated using experimental stress-strain curves of Mg-4Al samples with three different grain sizes from which the Hall-Petch coefficients are extracted and compared with Hall-Petch coefficients predicted using original parameters from previous work. The predictability of the model is then evaluated for a Mg-4Al sample with different texture and three grain sizes subjected to loading in different directions. The calibrated parameters are then used for some parametric studies to investigate the variation of Hall-Petch slope for different degrees of simulated spread in basal texture, variation of Hall-Petch slope with loading direction relative to basal poles for a microstructure with strong basal texture, and variation of yield strength with change in grain morphology. The proposed approach to incorporate the micro Hall-Petch equation into the CPFE framework provides a foundation to quantitatively model more complicated scenarios of coupling between grain size, texture and loading direction in the plasticity of Mg alloys.

Crystal plasticity↗
DOE OSTI2021

Open-Loop Control of Adjustable Tuned Mass Dampers for Floating Wind Turbine Platforms

In this study, we demonstrate the ability of adjustable tuned mass dampers (TMDs) to reduce the platform motion of floating offshore wind turbines (FOWTs). The TMDs are located in the hull and provide control authority by varying the amount of water ballast and compressed air in a reservoir over time. The optimal TMD settings depend on the wind and wave conditions of the FOWT. We present an open-loop control scheme that changes the TMD natural frequency based on sea state and a simple method for estimating the peak wave period and significant wave height. The performance of this open-loop control is compared to ideal control of the TMDs (set knowing the exact sea state), a TMD with a constant natural frequency targeting the design load case with the greatest platform motion, and a baseline set of simulations of the platform without TMDs. The constant natural frequency case performs nearly as effectively as the actively controlled case, suggesting that the TMDs can be simply designed for extreme load conditions and provide consistent control over a range of environments; however, the methods used to find the optimal TMD parameters and estimate sea state statistics could be used in other aspects of FOWT control and design.

artificial intelligence↗
DOE OSTI2022

Recovery of MOF-5 from Extreme High-Pressure Conditions Facilitated by a Modern Pressure Transmitting Medium

Here mechanisms underlying the mechanically induced amorphization of metal-organic frameworks (MOFs) are of current interest, and both high-pressure experimentation and molecular dynamics simulations have been used to reveal the fundamentals of load bearing, deformation, and pressure-induced amorphization (PIA) in these highly porous materials. Unfortunately, MOFs are typically highly susceptibility to amorphization, which limits the conditions under which they can be processed and used. However, their flexible structures can be stabilized to high pressures by incorporating guest species into the framework matrix. In this study, a large molecule pressure transmitting medium (DAPHNE 7575 ) is used as a structure-fortifying guest species to stabilize the prototypical MOF -5 to ultra-high pressures ( >9 GPa) and enable the recovery of crystalline material upon decompression. Structural changes associated with the penetration of the pressure transmitting medium on compression are examined using a combination of high pressure synchrotron powder diffraction and molecular dynamics simulations. This work enhances the understanding of PIA in MOFs while showcasing a potential route for the stabilization of MOFs to surprisingly high-pressures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Xyce(™) Parallel Electronic Simulator v.7.5

The Xyce Parallel Electronic Simulator simulates electronic circuit behavior in DC, AC, HB, MPDE and transient mode using standard analog (DAE) and/or device (PDE) device models including several age and radiation aware devices. It supports a variety of computing platforms (both serial and parallel) computers. Lastly, it uses a variety of modern solution algorithms dynamic parallel load-balancing and iterative solvers.! ! Xyce is primarily used to simulate the voltage and current behavior of a circuit network (a network of electronic devices connected via a conductive network). As a tool, it is mainly used for the design and analysis of electronic circuits.! ! Kirchoff's conservation laws are enforced over a network using modified nodal analysis. This results in a set of differential algebraic equations (DAEs). The resulting nonlinear problem is solved iteratively using a fully coupled Newton method, which in turn results in a linear system that is solved by either a standard sparse-direct solver or iteratively using Trilinos linear solver packages, also developed at Sandia National Laboratories.

Source record↗
DOE OSTI2023

MultiLoad-GAN: A GAN-Based Synthetic Load Group Generation Method Considering Spatial-Temporal Correlations

This paper presents a deep-learning framework, Multi-load Generative Adversarial Network (MultiLoad-GAN), for generating a group of synthetic load profiles (SLPs) simultaneously. The main contribution of MultiLoad-GAN is the capture of spatial-temporal correlations among a group of loads that are served by the same distribution transformer. This enables the generation of a large amount of correlated SLPs required for microgrid and distribution system studies. Here, the novelty and uniqueness of the MultiLoad-GAN framework are three-fold. First, to the best of our knowledge, this is the first method for generating a group of load profiles bearing realistic spatial- temporal correlations simultaneously. Second, two complementary realisticness metrics for evaluating generated load profiles are developed: computing statistics based on domain knowledge and comparing high-level features via a deep-learning classifier. Third, to tackle data scarcity, a novel iterative data augmentation mechanism is developed to generate training samples for enhancing the training of both the classifier and the MultiLoad-GAN model. Simulation results show that MultiLoad- GAN can generate more realistic load profiles than existing approaches, especially in group level characteristics. With little finetuning, MultiLoad-GAN can be readily extended to generate a group of load or PV profiles for a feeder or a service area.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Comparison of three DWM-based wake models at above-rated wind speeds

In this study we investigate three mid-fidelity wind turbine wake models based on the dynamic wake meandering (DWM) model principle, and compare their performance with a reference dataset, produced with large-eddy simulations using the actuator line model. The models are compared with respect to flow field, power, and loads on a row of four 5MW reference turbines experiencing above-rated wind conditions. In general, the DWM models show fairly good agreement with large-eddy simulation for the time-averaged flow fields, blade forces and power, with increasing differences along the turbine row. Also when comparing fatigue loads of blade root moments, the differences between the models increase further into the row, with deviations up to 25 % of the reference case. However, while the development in blade root moment fatigue along the turbine row is predominantly driven by the energy content at the frequency corresponding to the turbine's rotational period (1 P ) for the DWM models, the large-eddy simulation results suggest that the key drivers for the blade root and tower loads are the increase in meandering and energy at higher frequencies (> 1 P ) deeper into the turbine row. For the tower loads, the DWM models highly underestimate the fatigue for the waked turbines. From these results, we suggest priorities for future model developments so that robust model implementations can be used in wind farm design and operation.

17 WIND ENERGY↗
DOE OSTI2020

Cesium Removal Performance Comparisons of Crystalline Silicotitanate Media Batches with Savannah River Site Waste Simulant

The Tank Closure Cesium Removal (TCCR) system uses ion exchange columns filled with Crystalline Silicotitanate (CST) media to process radioactive waste solutions for the removal of Cs137. TCCR currently focuses on dissolving Savannah River Site (SRS) Tank 10H waste (primarily sodium salt cake solids) within the tank followed by at-tank ion exchange column treatment. Plans are underway to prepare and install a second TCCR unit at SRS. Capacity and particle size differences exist between archived (IE-911) and more recently-prepared CST media batches (9120- B and 9140-B). Side-by-side comparison testing was performed to evaluate the cesium removal performance of each batch to aid in selecting the preferred CST batch and media characteristics to load into the second TCCR unit. Batch contact equilibrium and flow-through column tests have been conducted with three CST batches using an SRS Average Simulant. Simulant batch contact equilibrium cesium loading results for the three CST batches (including two lots of one batch) are provided in Table ES-1. The 35 °C data indicates that the archived IE-911 batch has a higher cesium capacity than recently-prepared CST media and that the minor TCCR CST 9120-B media lot (2099000035) has similar cesium removal performance to the major TCCR lot (2099000034). Tests conducted at 25 ºC for the archived IE-911 CST batch indicated lower cesium removal performance with this simulant batch than was observed recently with a different SRS Average Simulant batch. A dilution factor (DF) is typically utilized when modeling engineered CST media cesium loading performance to account for mass contributions from the binder material. In cases where CST performance is lower than expected, this factor includes corrections for the binder and for low performance. The ZAM (Zheng, Anthony, and Miller) Isotherm Model DF values are provided in Table ES-1 for each simulant batch contact result. DF values near 0.5 were determined for the 9120-B and 9140-B CST batch contact tests while DF values near 0.6 (20% higher) were determined for the tests with IE-911. These dilution factors are lower than recently observed with a different SRS Average Simulant batch (9140-B DF = 0.68; IE-911 DF = 1.0). Based on these results, it appears that some component in the simulant solution used for equilibrium testing may have resulted in reduced cesium loading on the CST media.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Dual-season thermal energy storage-integrated heat pump system for load shifting: experiments, modeling and nationwide analysis

Here, this paper proposes and validates a dual-season thermal energy storage-integrated heat pump (TES-HP) system that shifts building electric load in both cooling and heating modes using a single phase-change TES. The design uses the TES as a heat sink for the condenser in summer and a heat source for the evaporator in winter, thereby reducing on-peak electric demand while ensuring compatible with existing air-distribution systems. A 14-kW TES-HP prototype with six operating modes and novel refrigerant charge management strategies was experimentally tested under varied conditions. Subsequently, data-driven polynomial performance curves were developed and validated against experiments, then coupled with a phase change material (PCM) model and a calibrated building model within a rule-based controller. Single-day case studies demonstrated effective on-peak demand reduction while maintaining thermal comfort. During cooling, hourly power was lowered by up to 1.5 kW; in heating, average hourly power decreased from over 3.5 kW to 1.8 kW, compared to an HP-only system, achieving 64.5 % load shifting, and the need for resistance heating was eliminated. Seasonal simulations showed typical on-peak electricity savings of 25–35 % in cooling and 40–65 % in heating, with the largest benefits on the hottest and coldest days. Extended response-surface analysis and nationwide mapping quantified load shifting as functions of ambient and TES temperatures, revealing a seasonal trade-off in TES phase-change temperature. These results demonstrate the TES-HP system as a practical and scalable solution for grid-interactive buildings that reduces on-peak demand and electricity use while maintaining thermal comfort and enhancing grid flexibility and reliability.

Grid-interactive buildings↗
DOE OSTI2021

The Simons Observatory: modeling optical systematics in the Large Aperture Telescope

We present geometrical and physical optics simulation results for the Simons Observatory Large Aperture Telescope. This work was developed as part of the general design process for the telescope, allowing us to evaluate the impact of various design choices on performance metrics and potential systematic effects. The primary goal of the simulations was to evaluate the final design of the reflectors and the cold optics that are now being built. We describe nonsequential ray tracing used to inform the design of the cold optics, including absorbers internal to each optics tube. We discuss ray tracing simulations of the telescope structure that allow us to determine geometries that minimize detector loading and mitigate spurious near-field effects that have not been resolved by the internal baffling. We also describe physical optics simulations, performed over a range of frequencies and field locations, that produce estimates of monochromatic far-field beam patterns, which in turn are used to gauge general optical performance. Finally, we describe simulations that shed light on beam sidelobes from panel gap diffraction.

Gudmundsson, Jon E. (ORCID:0000000317600355)↗
DOE OSTI2023

Extrapolation of the Rainflow-Counted Load Ranges for Fatigue Assessment of the Wind Turbine's Blades

Wind turbine design standards recommend the use of statistical modeling coupled with extrapolation of the short-term load data to long-term periods for fatigue reliability assessment. However, statistical error and computational expense can limit the accuracy of such approaches. In the case of wind turbine blades, the errors are more significant because of the high material fatigue exponent that makes the damage estimations more sensitive to variations. In addition, due to different excitation sources, the flapwise load range histogram is not unimodal, and thus its statistical modeling is complex. In the present work, we provide three methods for statistical modeling of the flapwise bending moment ranges including a novel approach based on frequency-based separation of the modes. The first two methods are simplified approaches for modeling the most crucial load ranges using unimodal distributions and the third method involves multimodal distribution fitting. The research is based on 3600 10-minute aeroelastic simulations of DTU 10MW case study wind turbine from which a benchmark damage equivalent load (DEL) is calculated. The DEL calculated by each of the three proposed methods is compared to this reference. The results show that the conventional approach based on using 6 seeds as well as using mixture models fitted on the limited data lead to under-conservative results with errors up to 23%. On the other hand, the simplified unimodal approaches provided in this work can provide conservative estimations of the fatigue damage with mean values 5% and 12% higher than the benchmark. However, the variability of the DEL estimates is higher when using unimodal extrapolation of the load ranges, and the data can be conservative by 17.5%. The proposed unimodal fits suggested for modeling and extrapolation of the blade's load ranges provide less errors relatively and most importantly conservative DEL estimations while maintaining computational efficiency.

blade fatigue↗
DOE OSTI2020

The Simons Observatory: Modeling Optical Systematics in the Large Aperture Telescope

We present geometrical and physical optics simulation results for the Simons Observatory Large Aperture Telescope. This work was developed as part of the general design process for the telescope; allowing us to evaluate the impact of various design choices on performance metrics and potential systematic effects. The primary goal of the simulations was to evaluate the final design of the reflectors and the cold optics which are now being built. We describe non-sequential ray tracing used to inform the design of the cold optics, including absorbers internal to each optics tube. We discuss ray tracing simulations of the telescope structure that allow us to determine geometries that minimize detector loading and mitigate spurious near-field effects that have not been resolved by the internal baffling. We also describe physical optics simulations, performed over a range of frequencies and field locations, that produce estimates of monochromatic far field beam patterns which in turn are used to gauge general optical performance. Finally, we describe simulations that shed light on beam sidelobes from panel gap diffraction.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Probabilistic Resilience of DER Systems - A Simulation Assisted Optimization Approach

Energy systems resilience is becoming increasingly important as the frequency of major grid outages increases. In this work, we present a methodology to optimize a behind-the-meter distributed energy resource system to sustain a site's critical loads during a pre-defined outage period. With the fixed system design, we then propose an outage simulation approach to estimate the resilience potential of the DER system to sustain loads beyond the fixed outage period - a yearlong resilience performance analysis. We apply statistical analysis to assess the system's resilience performance over a broader parametric problem space on an hourly, monthly, and yearly basis. We demonstrate the impact of the pre-defined outage period on the resilience performance through a case study. Results show that the probability of surviving a random outage of a given duration changes from 20% to 95% when the outage is modeled for a weekday instead of a weekend for the given load-profile.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Structures of fluoride containing aluminosilicate low activity nuclear waste glasses: A molecular dynamics simulations study

Fluorine and other halides commonly exist in nuclear waste forms, and due to their volatile nature, halide retention poses an issue affecting waste loading during vitrification. The compositional effect on fluorine incorporation in aluminosilicate glasses is investigated through molecular dynamics simulations. Oxygen and fluorine coordination numbers around glass former and modifier cations, bond angle distributions, and medium range structure features such as Q n distributions, ring size distributions and neutron diffraction structure factors were calculated. It was found that fluorine has higher preference to bond to Ca 2+ than to Na + , both in the melt and the glass, and there is no Si-F bond formation in the glass but they do exist in the melt. Consequently, CaO for Na 2 O substitution can be an effective way to help fluorine retention without significantly changing the glass chemistry. Furthermore, these results thus provide insights on fluorine incorporation in the aluminosilicate nuclear waste glasses and the strategy on how to improve fluorine retention both in the glass and the melt.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Stabilizing the Grid and Reducing Utility Bills Through Price-Responsive Controls for Heat Pump Water Heaters

The electricity grid is facing increasing challenges in cost-effectively balancing supply and demand. These challenges are exacerbated by increased penetration of photovoltaics, which causes mid-day overproduction, and electrification of gas appliances, which increases peak-period electricity demand. Decarbonization requires shifting building loads from fossil-intensive high-cost times to renewable-intensive low-cost times while maintaining quality of service to occupants. Utilities and ISO’s are investigating new ways of incentivizing this load shifting. One promising method is the use of Highly Dynamic Prices (HDPs). HDPs feature continuously changing prices that reflect real-time grid generation and distribution costs and capacity constraints, and thus incentivize consumers to shift their loads. California’s CPUC CalFUSE proposal and Hawaii’s recent changes demonstrate that electricity tariffs are moving towards this model. For this to work however, loads must have the capability to respond to these prices. Heat pump water heaters (HPWHs) are an ideal device for this purpose because their storage tanks decouple delivery of domestic hot water from electricity consumption. The storage enables control strategies that consume midday solar power to increase the energy stored in the tank, then provide evening peak domestic hot water services using the stored energy. Berkeley Lab's CalFlexHub project is pioneering price-driven load flexibility by developing and deploying cost-minimizing controls for many flexible loads - including HPWHs - in response to HDP. Control development is based on simulations using the Flexible Heat Pump Water Heater Performance Predictor which captures the control decisions of the on-board controller in a residential, integrated HPWH . The price-responsive controls a) shift load in ways that consume additional midday solar power to help stabilize the grid and reduce overall emissions, b) ensure that occupants receive equal or better hot water delivery service, and c) minimize the operating cost for each home in the fleet. On the grid level, the resulting shift will reduce utility operating costs and emissions, and can avoid expensive system capacity expansions. The control approach is customized to each home based on typical hot water consumption patterns. HPWH controllers, whether on the device or remotely, will receive a schedule of CTA-2045-B signals or set temperature adjustments customized to the current HDP price schedule and home. Simulation results for a fleet of 148 HPWHs on a summer day in Berkeley, California show cost savings of 29% and high price electricity consumption reductions of 80%, while maintaining full quality of service.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Probabilistic Resilience of DER Systems - A Simulation Assisted Optimization Approach

Energy systems resilience is becoming increasingly important as the frequency of major grid outages increases. In this work, we present a methodology to optimize a behind-tlie-meter distributed energy resource system to sustain a site's critical loads during a pre-defined outage period. With the fixed system design, we then propose an outage simulation approach to estimate the resilience potential of the DER system to sustain loads beyond the fixed outage period - a yearlong resilience performance analysis. We apply statistical analysis to assess the system's resilience performance over a broader parametric problem space on an hourly, monthly, and yearly basis. We demonstrate the impact of the pre-defined outage period on the resilience performance through a case study. Results show that the probability of surviving a random outage of a given duration changes from 20% to 95% when the outage is modeled for a weekday instead of a weekend for the given load-profile.

28 EE - Advanced Manufacturing Office (EE-5A)↗
DOE OSTI2022

Enhancement of phase change material hysteresis model: A case study of modeling building envelope in EnergyPlus

Nowadays, buildings are expected to offer demand side services to the power grid to enhance the electrical load flexibility, which leads to the concepts of grid-interactive efficient buildings (GEBs). Phase change material (PCM)-based thermal energy storage has seen increasing attention in recent years for peak load shifting of grid-interactive efficient buildings (GEBs). Numerical models are critical tools for design and evaluation of PCM-integrated systems. Most industrial-grade PCMs are reported to melt/freeze over a temperature range instead of at a unique temperature. Such thermal hysteresis effect significantly affects the reliability of simulation results because not only the heat transfer process depends on melting and freezing temperatures, the PCM thermal properties change significantly during the phase change process as well. This study is aimed to develop a model for the PCMs used in the building envelope with the capability to accurately simulate hysteretic behaviors. Further, this model is based on a two-phase assumption and is implemented in a whole building energy performance simulation program (i.e., EnergyPlus). A comparison between numerical results and experimental data shows that during a complete phase transition, the two-phase model could achieve a good agreement with the experimental data. During a partial phase transition, the two-phase model could lead to significant improvements compared to other alternative PCM models, including the existing PCM model in EnergyPlus. Last, whole building simulations were performed to study this model's performance regarding heating/cooling loads and zone mean air temperature of a given building. The results show that the difference in hourly heating/cooling loads introduced by the models was less than 1% in design conditions, while significant changes were observed in both hourly heating/cooling loads and zone mean air temperature when the PCM envelope underwent partial phase transition processes.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

High-performance and high-fidelity Monte Carlo solutions to the BEAVRS benchmark

The BEAVRS (Benchmark for Evaluation and Validation of Reactor Simulation) benchmark is solved by PRAGMA, the GPU-based continuous energy Monte Carlo code. The resulting solutions are comprised of the detailed simulation results of two cycles, each of which consists of the zero power physics test (ZPPT) and the core depletion calculations. The ZPPT consists of characteristic parameters, such as critical boron concentration (CBC), control rod bank worth, isothermal temperature coefficients, and assembly-wise detector signal, which are compared with measured data provided by the benchmark administration. The core depletion calculations were performed for both the hot full power and the load follow modes, and the comparison was made with the measured or deduced CBCs and assembly-wise detector signals. In the load follow calculations, the operating power history was approximately applied to simulate the real operation as closely as possible. PRAGMA performed the various calculations with a tremendous number of histories ranging up to hundreds of millions per cycle, exploiting GPUs' massively parallel performance. The load follow run time was shorter than 16 hours on a single rack of computing nodes mounded with 24 gaming GPUs with a remarkable agreement with the measurements for most comparisons. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Simulated aging processes of black carbon and its impact during a severe winter haze event in the Beijing-Tianjin- Hebei region

Black carbon (BC) can mitigate or worsen air pollution through perturbing meteorological conditions. BC aging processes are important for the evolution of particle size, concentration, and optical properties of BC that determines its influence on the meteorology. Here we use the online coupled Weather Research and Forecasting-Chemistry (WRF-Chem) model to quantify the role of BC aging processes, including physical processes (PP) and absorption enhancements (AE), in exerting BC-induced meteorological changes and the associated feedback to PM2.5 (particulate matter less than 2.5 µm in diameter) and O3 concentrations during a severe haze event in Beijing-Tianjin-Hebei (BTH) region during 21-27 February 2014. Our results show that, compared to the simulation without the PP treatment, the simulated near-surface BC concentration and BC mass loading in BTH region is lowered by 6.6 % and 12.1 %, respectively, during the haze event with the PP included. PP increases the proportion of large-size BC (particle diameter greater than 0.312 µm) from 28 %-33 % to 59 %-64 % below 1000 m in BTH region. Both PP and AE enhance the “dome effect” of BC. With both PP and AE considered, a reduction in PBL height due to BC-PBL interaction is 116.3 m (20.7 %), compared to 75.7 m (13.5 %) without AE and 66.6 m (11.9 %) without both PP and AE. However, during this haze event, anomalous northeasterly winds are produced by the direct radiative effect of BC, which further affects the mixing and transport of aerosols. When PBL height is decreased (from 10:00 to 21:00), combining all the impacts on multiple meteorological factors, BC effect without PP and AE, without AE, and with PP and AE, respectively, increases surface concentrations of PM2.5 by 8.3 µg m-3 (6.1 % relative to the mean value), 6.1 µg m-3 (4.5 %) and 9.6 µg m-3 (7.0 %) but decreases surface O3 concentrations by 2.8 ppbv (7.4 %), 4.0 ppbv (9.0 %) and 5.0 ppbv (10.8 %) in BTH region averaged over 21-27 February 2014. Our results highlight the importance of the aging processes and absorption enhancements of BC in simulating weather and air quality.

Chen, Donglin↗