Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamical balances”

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 325 records · Page 18

Amphiphile conformation impacts aggregate morphology and solution structure across multiple lengthscales

Although the self-assembly of amphiphiles is well-studied in aqueous solutions, much less is understood about the fundamental driving forces and structure property relationships in non-polar media. In recent work [Journal of Physical Chemistry B, 2020, 124, 10822.] the authors have studied a series of malonomide-based amphiphiles that are relevant to liquid-liquid extraction. That work demonstrated that aggregation is largely driven by local dipole-dipole interactions between molecules. Here, we build upon this observation to develop a more detailed understanding of how the balance of dipole-dipole interactions (controlled by conformation) and molecular architecture influences the morphology of the aggregates across lengthscales. Using constrained molecular dynamics about key degrees of freedom, we demonstrate that the conformation of N,N-dimethyl,N,N-dioctylhexylethoxy malonamide (DMDOHEMA) and N,N-dimethyl,N,N-dibutyltetradecyl malonamide (DMDBTDMA) has a significant impact upon self-association - where appropriate conformational sampling is essential. To quantify the aggregate morphology, several graph theoretic and persistent homology based properties are determined. The former examines the patterns of intermolecular interactions within clusters, while the latter examines the 3-dimensional spatial distribution. Based upon these analyses, we find that the morphology of aggregates, particularly at higher malonamide concentration, depends on a balance of dipole alignment and alkyl tail sterics. Dipole alignment encourages linear patterns of the intermolecular interactions within aggregates, while the alkyl tail steric interactions between the malonamide result in noticeably less linear aggregates for DMDOHEMA than DMDBTDMA. This is reflected in the spatial distribution, where more holes or voids exist between extractants within the DMDOHEMA that distribute within the solution in more of a "swiss cheese" arrangement as opposed to the more filamentous distribution of DMBDTDMA. This study links conformation and molecular structure to the morphology of amphiphile assemblies, and serves as a basis for ongoing study of multicomponent amphiphile solutions with polar and other solutes, and how these impact aggregation phenomena.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials: Part II Resveratrol Exemplar

This SAND report summarizes work supported by an Engineering Sciences Research Foundation (ESRF) Lab Directed Research and Development (LDRD) project entitled “Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials.” This SAND report is written in two parts with Part 1 discusses recrystallization of our explosive exemplar and Part 2 summarizing our work with recrystallization of resveratrol. We have studied resveratrol recrystallization with a multiscale approach combining experiments, modeling and simulation. At the single crystal scale, microscopy experiments illuminate crystal time-dependent growth rates using advanced image analysis. Bench scale experiments were carried out to look at growth of multiple particles in a small reactor creating thousands of particles and analyzing the results with microscopy and μCT. For the modeling we combine kinetic Monte Carlo (kMC) models with subscale information from density functional theory (DFT) or molecular dynamics. This work is discussed in Part 1 and can also be found in a paper from the project discussing a coarse-grained kMC model specifically developed for resveratrol. For well-mixed systems, we have population balance equations (PBE) linked with species mass conservation forming a set of ordinary differential equations that can be solved quickly. For more complicated geometries, such as the vat crystallization used throughout the complex, a coupled computational fluid dynamic (CFD)/PBE method was developed to account for gradients in temperature and concentration and differences in crystallization rates throughout the domain. These simulations are more complex and require high performance computing. We present results for two cases: 5% seed fast cool with parameters fit to the well-mixed case and 5% seed slow cool using the same parameters. We show reasonable agreement with experiments though are particles are significantly larger than the experiments.

36 MATERIALS SCIENCE↗

AI-Assisted Conceptual Development of a Pre-Geometric Cosmological Model - An Exercise in AI-Assisted Conceptual Framework Generation, Paper I: Foundations and Replication Dynamics

We develop a pre-geometric cosmological framework in which existence is identified with a finite amount of unstructured energy possessing vibration as its only intrinsic property. This vibrational substrate occupies an open, bounded spectral interval $(\omega_{\min},\omega_{\max})$, ensuring finiteness of total energy and excluding infinitely stable configurations. The substrate evolves under two fundamental and competing tendencies---excitation, which amplifies coherence, and randomization, which scrambles it. Their balance produces a metastable unstructured regime in which rare fluctuations may form long-lived self-consistent spectral configurations. Because the substrate is finite and subject to competing order--disorder dynamics, no coherent configuration can be perpetually stable. We show that the only mechanism capable of sustaining long-lived organization is a replication instability: a coherent unit may reproduce into multiple offspring according to a general $1\!\to n$ rule. Replication consumes energy from the finite substrate, breaks the metastable symmetry, and induces a discrete notion of event time through the replication tick $\Delta\tau$. Temporal succession is defined through correlation ordering of spectral microstates, producing an intrinsic pre-causal structure. The compactness of the spectral domain imposes minimal and maximal timescales, bounds the internal coherence of emergent units, and limits their proliferation. These spectral constraints serve as precursors for the emergence of geometry, adjacency, and a limiting propagation speed, developed in subsequent papers of this series. Paper I provides the foundational axioms (PG1--PG9) governing the spectral substrate, its metastable dynamics, the formation of coherent units, and the necessity of replication, establishing a fully pre-geometric stage from which causal and geometric structure naturally emerge.

79 ASTRONOMY AND ASTROPHYSICS↗

Pseudospectral convex optimization for on-ramp merging control of connected vehicles

It can be a daunting task for human drivers to merge into highways because of the intricate vehicle negotiations and potential risk within limited time and space. Connected vehicle (CV) technologies could be a solution to this problem and offer many benefits to the road safety, traffic mobility, and energy efficiency. However, real-time optimal control of CVs is still an open challenge, due to the nonlinear vehicle dynamics, non-convex fuel consumption model, and highly dynamic uncertain inter-vehicle interactions. To tackle these issues, a novel real-time optimal control approach that balances the computational efficiency and solution optimality is proposed for the purpose of onboard application. To this end, the pseudospectral collocation method is integrated with a sequential convex programming approach to develop two new optimization algorithms, which are implemented within a model predictive control (MPC) framework to allow for real-time generation of optimal merging speed profiles. One algorithm leverages the line search technique to improve convergence, and the other benefits from the trust region method for better computational efficiency. The optimality and convergence process of both proposed algorithms are investigated by comparing their solutions with a popular non-linear solver. Furthermore, simulation results show that the proposed methods outperform the benchmark in terms of computational cost, fuel consumption, and traffic efficiency. In particular, the proposed fuel-economy merging rule can save 57.1% fuel consumption on average on four different traffic volumes. Meanwhile, the proposed optimal control algorithms can reduce 2.2% travel time on average comparing to the “first-in-first-out” merging rule.

33 ADVANCED PROPULSION SYSTEMS↗

Evaluating the Transport Mechanism for Lignocellulosic Inhibitors in Zymomonas mobilis

The conversion of lignocellulosic biomass to biofuels and bioproducts is limited by small inhibitory molecules generated during the pretreatment process. The diverse compounds within the cell wall hydrolysate including aldehydes, carboxylic acids, phenolics, and alcohols inhibit microbial growth by disrupting microbial membranes. Although Zymomonas mobilis combines high ethanol productivity with hopanoid-rich membranes that confer solvent tolerance, hydrolysate compounds still inhibit microbial growth. One proposed mitigation strategy is to pump these molecules out through active transport processes, reducing direct toxicity. However, if passive permeation rates are high, active transport would create a futile cycle where the exported molecules diffuse back in, creating a net drag on fitness. To investigate this hypothesis, we measure the passive permeation for thirty-three lignocellulose derived compounds using atomistic molecular dynamics simulations with enhanced sampling. Permeability coefficients determined spanned more than eight orders of magnitude driven largely by the balance between hydrophobicity and polarity. Hydrophobic inhibitors such as aromatic acids, aldehydes, and short chain alcohols readily partitioned into the membrane and showed high permeability consistent with rapid passive influx and strong toxicity. Most polar molecules have lower but still appreciable permeation rates, with only the most polar sugars having low permeabilities. By integrating these predictions together with hydrolysate composition information, we find that most of these molecules are so permeable that futile cycles where passive permeation overwhelms efflux is highly likely. These findings provide a predictive framework for strain engineering, highlighting that, for highly permeable inhibitors, reducing extracellular inhibitor concentrations through hydrolysate detoxification, pretreatment optimization, or dilution may be necessary to mitigate futile cycling and improve microbial productivity in industrial biorefineries.

Zymomonas mobilis↗

Observations of Subduction, Downward Heat Flux and Dense Filament Collapse in the Northern Gulf of Mexico

Submesoscale processes are important contributors to the global heat budget and generally support upward heat transport through restratification. However, in salinity‐stratified regions, such as the northern Gulf of Mexico with its influx of freshwater from the Mississippi‐Atchafalaya river system, temperature can act like a passive tracer and submesoscale processes can contribute to downward heat transport. Oceanic heat content is a factor in many environmental risks the region faces, for example, hurricane intensification, and marine heatwaves. During the 2022 field campaign of the Submesoscales Under Near‐Resonant Inertial Shear Experiment, a sampling plan was developed to study such submesoscale processes in high resolution. Over 31 hr, four assets (two research ships and two remotely controlled boats) drove in parallel across a dense filament, capturing its evolution in time and space. The observations show that surface waters, warmed by daytime solar radiation, were subducted and that the associated overturning circulation transported heat below the surface layer where it was later irreversibly mixed away. The estimated downward heat flux was as strong as the concurrent net air‐sea heat flux into the ocean. The filament was then observed to rapidly collapse which we attribute to boundary layer turbulence and the breakdown of geostrophic balance. The collapsing fronts display behaviors indicative of gravity currents. These observations highlight how in salinity‐stratified regions, frontal dynamics can be associated with downward heat flux and how the submesoscale can play an important role in the oceanic heat budget.

58 GEOSCIENCES↗

Estimating the carbon footprint of computational fluid dynamics

In this study, computational resources have grown exponentially in the past few decades. These machines make possible research and design in fields as diverse as medicine, astronomy, and engineering. Despite ever-increasing computational capabilities, direct simulation of complex systems has remained challenging owing to the degrees of freedom involved. At the cusp of exascale computing, high-resolution simulation of practical problems with minimal model assumptions may soon experience a renaissance. However, growing reliance on modern computers comes at the cost of a growing carbon footprint. To illustrate this, we examine historic computations in fluid dynamics where larger computers have afforded the opportunity to simulate flows at increasingly relevant Reynolds numbers. Under a variety of flow configurations, the carbon footprint of such simulations is found to scale roughly with the fourth power of Reynolds number. This is primarily explained by the computation cost in core-hours, which is also described by similar scaling, though regional differences in renewable energy use also play a role. Using the established correlation, we examine a large database of simulations to develop estimates for the carbon footprint of computational fluid dynamics in a given year. Collectively, the analysis provides an additional benchmark for new computations where, in addition to balancing considerations of model fidelity, carbon footprint should also be considered.

42 ENGINEERING↗

Resilient Inverter-Driven Black Start with Collective Parallel Grid-Forming Operation

As modern power systems are experiencing exceptional changes with increasing penetrations of inverter-based resources (IBRs), system restoration using IBRs has received attention. Using local grid-forming (GFM) assets near consumers, engineered to establish grid voltages in the absence of a stiff grid, i.e., bottom-up restoration, a distribution system could obtain high system resilience by not relying on the bulk power system restoration, which requires significant human intervention and procedure. This paper studies the technical feasibility of the novel approach with detailed electromagnetic transient (EMT) simulations. To thoroughly evaluate the potential of GFM inverters and the technical challenges in IBR-driven black start, a detailed three-phase inverter model is developed, including negative-sequence control for voltage balance and a phase-by- phase current limiter to sustain momentary overloading during the black start. To examine dynamic aspects of the black-start process, the EMT simulation also models transformer and motor dynamics to emulate their inrush and startup behaviors as well as network dynamics. In addition, active involvement of grid- following distributed energy resources is also studied to facilitate the black-start process. By allowing multiple GFM inverters to collectively black start without leader-follower coordination, we demonstrate that a system can achieve high resilience even with a fraction of assets lost. Two test cases of inverter-driven black start, using two and one GFM inverters, respectively, for a heavily unbalanced 2-MVA distribution feeder are demonstrated. Takeaways for further study and field deployment are provided.

grid-forming inverter↗

Metabolic Feedback Inhibition Influences Metabolite Secretion by the Human Gut Symbiont Bacteroides thetaiotaomicron

ABSTRACT Microbial metabolism and trophic interactions between microbes give rise to complex multispecies communities in microbe-host systems. Bacteroides thetaiotaomicron ( B. theta ) is a human gut symbiont thought to play an important role in maintaining host health. Untargeted nuclear magnetic resonance metabolomics revealed B. theta secretes specific organic acids and amino acids in defined minimal medium. Physiological concentrations of acetate and formate found in the human intestinal tract were shown to cause dose-dependent changes in secretion of metabolites known to play roles in host nutrition and pathogenesis. While secretion fluxes varied, biomass yield was unchanged, suggesting feedback inhibition does not affect metabolic bioenergetics but instead redirects carbon and energy to CO 2 and H 2 . Flux balance analysis modeling showed increased flux through CO 2 -producing reactions under glucose-limiting growth conditions. The metabolic dynamics observed for B. theta , a keystone symbiont organism, underscores the need for metabolic modeling to complement genomic predictions of microbial metabolism to infer mechanisms of microbe-microbe and microbe-host interactions. IMPORTANCE Bacteroides is a highly abundant taxon in the human gut, and Bacteroides thetaiotaomicron ( B. theta ) is a ubiquitous human symbiont that colonizes the host early in development and persists throughout its life span. The phenotypic plasticity of keystone organisms such as B. theta is important to understand in order to predict phenotype(s) and metabolic interactions under changing nutrient conditions such as those that occur in complex gut communities. Our study shows B. theta prioritizes energy conservation and suppresses secretion of “overflow metabolites” such as organic acids and amino acids when concentrations of acetate are high. Secreted metabolites, especially amino acids, can be a source of nutrients or signals for the host or other microbes in the community. Our study suggests that when metabolically stressed by acetate, B. theta stops sharing with its ecological partners.

59 BASIC BIOLOGICAL SCIENCES↗

Resilient Inverter-Driven Black Start with Collective Parallel Grid-Forming Operation: Preprint

As the modern power systems are experiencing exceptional changes with increasing penetration of inverter-based resources (IBRs), system restoration using IBRs has received attention. Using local grid-forming (GFM) assets near consumers, engineered to establish grid voltages in the absence of a stiff grid, i.e., bottom-up restoration, a distribution system would obtain high system resilience, by not relying on the bulk power system restoration requiring significant human intervention and procedure to restore. This paper studies the technical feasibility of the novel approach with detailed electromagnetic transient (EMT) simulations. To thoroughly evaluate the potential of GFM inverters and technical challenges in IBR-driven black start, a detailed three-phase inverter model is developed, including negative-sequence control for voltage balance and a phase-by-phase current limiter to sustain momentary overloading during the black start. To examine dynamic aspects of the black start process, the EMT simulation also models transformer and motor dynamics to emulate their inrush and start-up behaviors as well as network dynamics. In addition, active involvement of grid-following distributed energy resources (DER) is also studied to facilitate the black start process. It is shown that, by allowing multiple GFM inverters to collectively black start without master-slave coordination, a system can achieve high resilience even with a fraction of assets lost. Two test cases of inverter-driven black start using two and one GFM inverters, respectively, for a heavily unbalanced 2-MVA distribution feeder are demonstrated. Takeaways for further study and field deployment are provided.

grid-forming inverter↗

Regional inertia dynamics of U.S. interconnections: An event-based measurement approach

Power grid inertia plays a vital role in frequency stability following large disturbances, yet its distribution across the U.S. grid is highly uneven. While interconnection-wide inertia benchmarks are useful, they can mask regional variability driven by resource mix, network coupling, and geographic separation. This paper extends event-driven inertia estimation to the regional scale using field measurements from the Frequency Monitoring Network (FNET/GridEye). Starting from balancing authority and independent system operator footprints, candidate regions are refined using a composite coherency score that combines frequency-trajectory shape similarity, timing spread, and lead/lag behavior to ensure dynamic consistency. A filtered sliding difference method (FSDM) is then used to construct regional frequency trajectories, detect disturbance onset, and compute robust regional rate-of-change of frequency (RoCoF). Regional, local, and interconnection inertia are estimated by combining RoCoF with event power imbalance, and additional indicators (regional-to-system inertia ratio and inertial-support arrival time) quantify regional-to-interconnection coupling and relative regional contributions. The method is demonstrated on eleven regions across the Eastern Interconnection (EI) and the Western Electricity Coordinating Council (WECC), with the Electric Reliability Council of Texas (ERCOT) used for validation. In ERCOT, estimates compared against energy management system (EMS) values achieve a mean absolute percentage error of 17.94%. WECC exhibits consistently shorter inertial-support arrival times (0.15–0.3 s) than EI (0.7–1.1 s), highlighting contrasting coupling and disturbance-propagation behavior. Overall, the results reveal pronounced spatial heterogeneity in inertia and coupling, underscoring the value of regional monitoring for both operational decision-making and long-term system planning.

Disturbance events↗

A dynamic 2D Borehole Thermal Energy Storage (BTES) model for enhanced computational efficiency

Progressing toward a future increasingly reliant on renewable energy sources, the development of effective, durable energy storage solutions becomes essential to balance supply and demand fluctuations. Borehole Thermal Energy Storage (BTES) is a long-duration thermal energy storage technology that captures excess heat generated from renewable energy sources and stores it underground for later use, enabling the efficient utilization of sustainable energy. This approach is particularly valuable in district energy networks when integrated with Ground Source Heat Pumps (GSHP) to provide stable heating and cooling. However, traditional three-dimensional (3D) numerical models of BTES systems demand extensive computational resources, limiting their practicality for real-time and large-scale applications. This study introduces a novel two-dimensional (2D) modeling approach that reduces computational costs while maintaining high accuracy. By employing a radial ring-based discretization method, the model simulates heat injection, retention, and retrieval dynamics over seasonal cycles. A new thermal-mass weighted-average temperature parameter is introduced to evaluate the performance of BTES systems. Model validation against FEFLOW simulations demonstrates a 17-fold improvement in computational speed compared to traditional Computational Fluid Dynamics (CFD) models while achieving a mean absolute percentage error (MAPE) of 2 % during charging and 4 % during discharging. Additionally, a trade-off analysis between computational efficiency and accuracy is conducted, ensuring the model's applicability for real-world scenarios. The findings of this research contribute to the development of computationally efficient BTES models, facilitating better optimization, control, and integration into renewable energy systems. This work provides a foundation for further studies in techno-economic analysis, multi-year performance evaluation, and real-time operational strategies for BTES applications, supporting a more sustainable energy future.

2D modeling↗

New Particle Formation and Growth Dynamics for α-Pinene Ozonolysis in a Smog Chamber and Implications for Ambient Environments

Secondary organic aerosol (SOA) plays an important role in the Earth’s energy balance and air quality. Understanding the formation of SOA accompanying the oxidation of volatile organic compounds (VOCs) remains incomplete because of the complex reactions and partitioning mechanisms. In this study, we investigated new particle formation (NPF) of SOA derived from α-pinene ozonolysis in the Harvard Environmental Chamber using the Model for Simulating Aerosol Interactions and Chemistry (MOSAIC) with the condensing species as two products, one of lower volatility (LVOC, 0.01 μg m –3 ) and one of higher volatility (SVOC, 1.0 μg m –3 ). A module for the classical nucleation theory was included. The surface tension (σ), bulk diffusivity (D b ), and mass accommodation coefficient (α) were varied until convergence between the simulated and observed particle number size distributions. α of LVOC was constrained as 0.03–0.1 with the reported σ of α-pinene-ozonolysis-derived SOA and Tolman surface correction (~23.0–27.5 mN m –1 ) and Db from 10 –15.5 to >10 –11 cm 2 s –1 . The observed particle number size distribution of α-pinene ozonolysis was well-simulated to derive the physical parameters, including σ, Db, and α, suggesting the possible application to group the reported complex oxidation products and further application for other atmospheric-related oxidations for quantifying NPF.

54 ENVIRONMENTAL SCIENCES↗

Laser heating and evaporation of a single droplet

The laser technology is being abundantly studied for controlled energy deposition for a range of applications in aerodynamic flow control, material processing, ignition, and combustion. The absorption of laser radiation by liquid droplets affects further propagation of laser in the atmosphere and causes bleaching of suspended droplets while the ignition and combustion characteristics in combustors are influenced by the evaporation rate of the sprayed fuel. In this work, we present a multi-dimensional mathematical model built on OpenFOAM for laser heating and evaporation of a single droplet in the diffusion dominated regime taking into account absorption of the laser radiation, evaporation process and vapor flow dynamics. The developed solver is validated against available experimental and numerical data for the ethanol and water droplet heating and evaporation. For continuous heating the peak temperature is established by the balance of cooling, evaporation and heating and results in high temperature for larger droplets. It has been shown that for heating by a single laser pulse the maximum temperature of droplets depends only on the peak intensity of the laser radiation. Furthermore, for the peak irradiance close to the transition to the boiling regime, temporal dynamics of the droplet temperature is independent of the droplet size. With proper normalization of time, the dynamics of the droplet shrinkage and cooling is shown to be independent of droplet sizes and peak laser intensities. The influence of cooling and evaporation processes on droplet heating was found to be controlled by the pulse repetition rate for repeated pulse operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dynamic Analysis of a Six-Cable Parallel Robot for Automated Panelized Building Retrofits

Cable-Driven Parallel Robots (CDPRs) are highly suitable for automated panelized building retrofits, thanks to their compact footprint and high payload-to-weight ratio. A common CDPR configuration featuring eight cables, where the anchors form a rectangular prism in front of the building facade, offers a large wrench feasible workspace and good control versatility. However, installing upper anchor points requires additional support structures, such as towers or beams, increasing setup complexity and posing logistical challenges in construction settings. To mitigate these challenges, we propose a six-cable CDPR model specifically designed for automated panelized building retrofits. Although the feasible workspace is limited, our analysis shows that the proposed CDPR adequately covers the critical areas required for panel installation. To validate that the six-cable system can effectively transport the end effector to the desired installation pose, we calculated optimal trajectories based on a constrained dynamic model. The simulation results of the six-cable CDPR demonstrate promising potential for automated panelized building retrofits, effectively balancing simplicity, cost-effectiveness, and functionality.

Liu, Yifang [ORNL] (ORCID:0000000190817417)↗

IoT-based retrofit information diffusion in future smart communities

Community-scale building retrofits are not merely scaled-up versions of single-building retrofits. They involve complex challenges, such as reconciling individual interests with collective goals and managing the dynamic interplay between buildings through mechanisms like power grids and social connections. Internet of Things (IoT) connectivity holds the potential to leverage these interplays to balance individual and collective interests effectively in smart communities. One critical aspect of this interplay is information diffusion, which shapes how retrofit decisions spread among neighbors, influencing individual choices and ultimately impacting community-level retrofit outcomes. In other words, IoT-based smart devices automatically push tailored retrofit notifications to homeowners, which completely changes the format of information diffusion in the future. To investigate this influence by such information diffusion, the study used CityBES to simulate energy performance for different retrofits and applied an information diffusion model to analyze how decisions spread in a networked community of 192 buildings. The diffusion process was modeled on a weighted, directed network, capturing the dynamics of information flow and decision-making across 16 scenarios. Individual retrofit benefits were evaluated through payback years, while community-level retrofit outcomes were assessed using greenhouse gas (GHG) emission reductions. The results demonstrate that easier information diffusion among neighbors encourages households to prioritize retrofit measures that align with the majority’s optimal choices, even at the expense of individual financial benefits. In this case, such collective prioritization enhanced community-level retrofit performance, increasing GHG emission reductions by up to 29.4 %. However, this improvement came with trade-offs, as the average payback period for households extended by approximately 1.74 years. These findings highlight the potential of IoT-based information diffusion in future smart communities to coordinate individual interests with collective goals, ultimately accelerating community-level building retrofits.

Shu, Lei↗

Arbitrary Order Virtual Element Methods for High‐Order Phase‐Field Modeling of Dynamic Fracture

ABSTRACT Accurate modeling of fracture nucleation and propagation in brittle and ductile materials subjected to dynamic loading is important in predicting material damage and failure under extreme conditions. Phase‐field fracture models have garnered a lot of attention in recent years due to their success in representing damage and fracture processes in a wide class of materials and under a variety of loading conditions. Second‐order phase‐field fracture models are by far the most popular among researchers (and increasingly, among practitioners), but fourth‐order models have started to gain broader acceptance since their more recent introduction. The exact solution corresponding to these high‐order phase‐field fracture models has higher regularity. Thus, numerical solutions of the model equations can achieve improved accuracy and higher spatial convergence rates. In this work, we develop a virtual element framework for the high‐order phase‐field model of dynamic fracture. The virtual element method (VEM) can be regarded as a generalization of the classical finite element method. In addition to many other desirable characteristics, the VEM allows computing on polytopal meshes. Here, we use ‐conforming virtual elements and the generalized‐ time integration method for the momentum balance equation, and adopt ‐conforming virtual elements for the high‐order phase‐field equation. We verify our virtual element framework using classical quasi‐static benchmark problems and demonstrate its capabilities with the aid of numerical simulations of dynamic fracture in brittle materials.

42 ENGINEERING↗

The mechanisms and seasonal differences of the impact of aerosols on daytime surface urban heat island effect

The urban heat island intensity (UHII) is the temperature difference between urban areas and their rural surroundings. It is commonly attributed to changes in the underlying surface structure caused by urbanization. Air pollution caused by aerosol particles can affect the UHII through changing (1) the surface energy balance by the aerosol radiative effect (ARE) and (2) planetary-boundary-layer (PBL) stability and airflow intensity by modifying thermodynamic structure, which is referred to as the aerosol dynamic effect (ADE). By analyzing satellite data and ground-based observations collected from 2001 to 2010 at 35 cities in China and using the WRF-Chem model, we find that the impact of aerosols on UHII differs considerably: reducing the UHII in summer but increasing the UHII in winter. This seasonal contrast is proposed to be caused by the different strengths of the ARE and ADE between summer and winter. In summer, the ARE on UHII is dominant over the ADE, cooling down surface temperature more strongly in urban areas than in rural areas because of much higher aerosol loading, and offsets the urban heating, therefore weakening UHII. In winter, however, the ADE is more dominant, because aerosols stabilize the PBL more in the polluted condition, weakening the near-surface heat transport over urban areas in both vertical and horizontal directions. This means that the heat accumulated in urban areas is dispersed less effectively, and thus the UHII is enhanced. These findings shed new light on the impact of the interaction between urbanization-induced surface changes and air pollution on urban climate.urban areas and their rural surroundings. It is commonly attributed to changes in the underlying surface structure caused by urbanization. Air pollution caused by aerosol particles can affect the UHII through changing (1) the surface energy balance by the aerosol radiative effect (ARE) and (2) planetary-boundary-layer (PBL) stability and airflow intensity by modifying thermodynamic structure, which is referred to as the aerosol dynamic effect (ADE). By analyzing satellite data and ground-based observations collected from 2001 to 2010 at 35 cities in China and using the WRF-Chem model, we find that the impact of aerosols on UHII differs considerably: reducing the UHII in summer but increasing the UHII in winter. This seasonal contrast is proposed to be caused by the different strengths of the ARE and ADE between summer and winter. In summer, the ARE on UHII is dominant over the ADE, cooling down surface temperature more strongly in urban areas than in rural areas because of much higher aerosol loading, and offsets the urban heating, therefore weakening UHII. In winter, however, the ADE is more dominant, because aerosols stabilize the PBL more in the polluted condition, weakening the near-surface heat transport over urban areas in both vertical and horizontal directions. This means that the heat accumulated in urban areas is dispersed less effectively, and thus the UHII is enhanced. These findings shed new light on the impact of the interaction between urbanization-induced surface changes and air pollution on urban climate.

54 ENVIRONMENTAL SCIENCES↗