Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mesoscale”

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 559 records · Page 31

Realistic large eddy and dispersion simulation experiments during project sagebrush phase 1

Typical dispersion simulations use mesoscale models at grid spacing greater than 1 km to simulate pollutant transport. Advances in computing power have allowed Large Eddy Simulations (LES), (grid spacing less than 100 m) to be run with boundary conditions in the grayscale regime (grid size of 1 km–100 m). We used the WRF mesoscale/LES model and the HYSPLIT particle dispersion model to simulate SF 6 transport from Project Sagebrush Phase 1, using a grayscale-aware scheme with 4 nests of 3 km, 1 km, 200 m, 67 m. The WRF/LES simulation was in satisfactory agreement with meteorological observations but the HYSPLIT SF6 simulations indicated insufficient sub-grid turbulence kinetic energy, especially in the lowest 100 m. This resulted in over predictions of the surface SF6 concentration and insufficient vertical transport. Best agreement with the Sagebrush observations was obtained with the WRF/LES large-scale and sub-grid turbulence, enhanced with a HYSPLIT parameterized turbulence profile.

54 ENVIRONMENTAL SCIENCES↗

Opportunities in multiscale modeling of mosquito-borne flaviviruses

Mosquito-borne flaviviruses, such as Zika, dengue, West Nile, and yellow fever virus, represent a growing public health concern due to their widespread distribution and the severe diseases they cause. These viruses are difficult to control as climate change and urbanization help mosquitoes expand into new areas, increasing the risk of outbreaks. Mathematical models play a key role in understanding their spread, providing insights at every level—from how the virus multiplies inside cells to how it circulates through entire populations. This review examines various approaches used in modeling arboviruses, including microscale models that focus on cellular and molecular dynamics, mesoscale models that address within-host processes, and macroscale models that capture population-level transmission. We briefly summarize the methodology used for models at each scale, which primarily consists of sets of differential equations with parameters that represent physical rates of change for different subprocesses. We particularly highlight how temperature affects virus transmission, which is key to understanding the impact of climate change. We also show how multiscale models can connect viral replication, immune response, and the spread of infection at a larger scale. This is essential for developing better vaccines and treatments, evaluating disease control measures, predicting the impact of climate change, and improving public health responses to outbreaks.

60 APPLIED LIFE SCIENCES↗

Integral boundary conditions in phase field models

Modeling the chemical, electric and thermal transport as well as phase transitions and the accompanying mesoscale microstructure evolution within a material in an electronic device setting involves the solution of partial differential equations often with integral boundary conditions. Employing the familiar Poisson equation describing the electric potential evolution in a material exhibiting insulator to metal transitions, we exploit a special property of such an integral boundary condition, and we properly formulate the variational problem and establish its well-posedness. Next, we compare our method with the commonly-used Lagrange multiplier method that can also handle such boundary conditions. Numerical experiments demonstrate that our new method achieves optimal convergence rate in contrast to the conventional Lagrange multiplier method. Furthermore, the linear system derived from our method is symmetric positive definite, and can be efficiently solved by Conjugate Gradient method with algebraic multigrid preconditioning.

97 MATHEMATICS AND COMPUTING↗

High-performance organic photodetectors enabled by a refined fibrillar multiphase morphology

Film morphology of donor:acceptor blend layers plays a critical role in photon-to-current efficiency and dark/noise current of organic photodetectors (OPDs). One effective approach to manipulate crystallization and mesoscale phase separation of such blend layers is sequential casting (SC). However, the guiding strategies to control the morphology and the impacts on OPD performance of the SC films remain elusive, as the film structural evolution during SC is different from conventional blend casting (BC). In this report a refined fibrillar multiphase morphology is demonstrated by SC in an NT812:IEICO-4F blend film, where the high-quality IEICO-4F crystallites infiltrate the robust NT812 fibril network from the surface to the bulk. Such a morphology leads to improved charge generation and collection, reduced trap states, and enhanced charge block capability of resultant device, enabling simultaneous achievement of high external quantum efficiency and low dark/noise current. A maximum special detectivity of 5 × 10 13 Jones is achieved at 860 nm under –0.1 V, which is among the highest detectivities for vis-to-NIR OPDs. The linear dynamic range and response speed are also improved. Such enhancements are parallelly observed from OPD devices based on other blend systems with similar fibrillar refinement, which provides guidelines for film structure manipulation towards high OPD performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The interpenetration polymer network in a cement paste–waterborne epoxy system

The formation of the interpenetrating polymer network (IPN) structure within the cement-polymer system has been revealed by experiments from 2 dimensional to 3-dimensional scales. However, the microstructure design and performance prediction of IPN as a function of specific C-S-H/polymer components or ingredients' parameters e.g. water to cement ratios (w/c), polymer to cement ratios (p/c), monomer components, degree of polymerization (DP), etc. is by far not available. Here we developed a mesoscale model for IPN visualization based on the Flory-Huggins interaction theory which was applied to cement science for the first time. All the ingredients in the micro-structure were considered as soluble beads with rational sizes based on their properties obtained by molecular dynamic methods. The interaction parameters of each bead were then determined based on their element ratios and chemical backbones. The model was validated with the waterborne epoxy-cement material (WECM) in which a novel waterborne epoxy resin (WEP) was prepared and impregnated. The verification results on the 2D-3D scale show that the developed model predicts the WECM's IPN structures by rule and line concerning DP, w/c, and p/c in the mixture. The C-S-H beads were progressively scaled in size which alters the C-S-H texture to have completely different dissolution characteristics. Beads of diameter ~6 Å are more unstable and soluble which enable them to form a continuous phase in water or a composite structure with WEP. In contrast, beads with diameters larger than 10 Å have different properties with stronger nucleation effects. The results also suggest the impregnation content of WEP in cement-based material should be limited to 10% vol. to prevent the polymer IPN from decomposing into discrete clusters. The application of Flory-Huggins theory in cement-based composites demonstrates great potential in performance prediction and microstructure design.

36 MATERIALS SCIENCE↗

Influence of the distribution of expansive sites in aggregates on microscopic damage caused by alkali-silica reaction: Insights into the mechanical origin of expansion

The origin of damage in concrete due to the alkali-silica reaction (ASR) is attributed to the expansion site in the aggregate. To investigate the cracking process of the aggregate during ASR and its consequences on concrete damage, the effect of the distribution of the expansive sites in the aggregate on ASR expansion and the crack patterns must be evaluated. Thus, in this study, a mesoscale discrete model was applied to ASR modeling to represent the propagation of cracking during ASR accurately. The distribution of the expansive sites in the aggregate was based on the gel pocket and reaction rim models, which are two ASR mechanisms reported in the literature. These two expansion models highlight the different crack patterns obtained based on the aggregate characteristics. Further, the expansion cracking processes determined based on the gel pocket and reaction rim models are consistent with the evolution of cracking with the expansion level.

36 MATERIALS SCIENCE↗

Chemical structure and complex growth modes of magnesium silicate hydrate: Nanoparticle orientation, aggregation, and fusion

The extent of utilization of magnesium silicate hydrate (MSH) in construction is limited partly because of insufficient data ascertaining the kinetics of its precipitation. In this report MSH grown homogeneously or heterogeneously in the range of magnesium to silicon solution molar concentration ratios, [Mg]/[Si] = 0.5-1.5, and temperature = 25-80 °C was analyzed for its chemical structure and morphology. Infrared spectroscopy and geochemical modeling show that increasing the [Mg]/[Si] ratio results in less silicate polymerization in MSH. Atomic force microscopy reveals the presence of nanoparticles that aggregate and fuse to form films. Oriented attachment of nanoparticles and features indicative of enhanced crystallinity were observed at higher temperatures and longer reaction times. These data provide direct evidence for the persistence of hierarchical structures (i.e., nanoparticles forming 3D microparticles and 2D film layers) at the nanoscale to mesoscale. These findings offer insights into the precise chemical synthesis of MSH and its widespread use as a binder for construction purposes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling and simulation of vascular tumors embedded in evolving capillary networks

Here in this work, we present a coupled 3D–1D model of solid tumor growth within a dynamically changing vascular network to facilitate realistic simulations of angiogenesis. Additionally, the model includes erosion of the extracellular matrix, interstitial flow, and coupled flow in blood vessels and tissue. We employ continuum mixture theory with stochastic Cahn–Hilliard type phase-field models of tumor growth. The interstitial flow is governed by a mesoscale version of Darcy’s law. The flow in the blood vessels is controlled by Poiseuille flow, and Starling’s law is applied to model the mass transfer in and out of blood vessels. The evolution of the network of blood vessels is orchestrated by the concentration of the tumor angiogenesis factors (TAFs); blood vessels grow towards the increasing TAFs concentrations. This process is not deterministic, allowing random growth of blood vessels and, therefore, due to the coupling of nutrients in tissue and vessels, makes the growth of tumors stochastic. We demonstrate the performance of the model by applying it to a variety of scenarios. Numerical experiments illustrate the flexibility of the model and its ability to generate satellite tumors. Simulations of the effects of angiogenesis on tumor growth are presented as well as sample-independent features of cancer.

3D–1D coupled blood flow models↗

Influence of Hydrogen Isotopes on Vacancy Formation and Antisite Defect Diffusion in Palladium and Vanadium Metals

Density functional theory calculations are performed to study the interaction between hydrogen isotopes, vacancy, and antisite defects in Pd and V. Various defect configurations and defect-defect distances are systematically explored. Binding energies and migration energy barriers are calculated and tabulated. The study provides atomistic data for subsequent mesoscale simulations of tritium, vacancy, and antisite diffusion. In Pd, a strong binding is found between a vacancy and tritium (0.16 eV). This tritium binding increases the thermal concentration of vacancies by a factor of ~10 at 500 ?C. The increase in vacancy concentration enhances V antisite diffusion in Pd by a factor of ~6 at 500 ?C. The influence of tritium is even stronger in V, with a tritium and vacancy binding energy of 0.38 eV. Such a strong binding increases the thermal concentration of vacancies by a factor of ~300 at 500 ?C. The increase in vacancy concentration enhances Pd antisite diffusion in V by a factor of ~640 at 500 ?C. Vanadium and Pd exhibit a strong driving force to intermix with a formation energy of -1.57 eV for V antisite in Pd and -1.05 eV for Pd antisite in V. The results suggest vanadium diffusion into Pd is energetically stronger than the reverse. Zero-point-energy corrections are taken into account and calculations for hydrogen and tritium are presented.

Setyawan, Wahyu↗

The impact of anisotropic thermal expansion on the isothermal annealing of polycrystalline α -uranium

Although grain growth impacts microstructural evolution in a wide variety of materials systems, the effect of anisotropic thermal expansion on grain boundary mobility and texture evolution has not been widely studied. Anisotropic thermal expansion occurs in multiple non-cubic metals, and the thermomechanical processing behavior of these materials can be better understood with further study into the impact of thermal expansion on grain boundary mobility and texture evolution. In this work, we develop a mesoscale phase field model of grain growth that includes the effect of anisotropic thermal expansion, which is applied to study polycrystalline α-uranium, a highly anisotropic metal. Three-dimensional simulations on polycrystalline a-uranium with and without thermal expansion eigenstrains are performed to study the grain boundary mobility and texture evolution as a function of temperature. A strain-free temperature of 933 K is selected, and the system is studied at within the range of 873 – 933 K at intervals of ten degrees, resulting in increasing thermal eigenstrain with decreasing temperature. We also estimate a grain boundary mobility prefactor and activation energy based on existing experimental data of isothermal annealing of α-uranium. The grain boundary mobility is found to display significant deviation from Arrhenius behavior with the inclusion of thermal expansion eigenstrain as the amount of thermal eigenstrain (and thus elastic strain energy within the system) increases. This result explains an experimentally observed grain boundary mobility deviation from Arrhenius behavior. Furthermore, the texture evolution is affected, such that the grain orientations become less random with increasing thermal eigenstrain, which could explain experimentally observed texture behavior. These results indicate that the effect of thermal expansion should be considered when predicting the thermomechanical processing behavior of α-uranium and other materials with anisotropic thermal expansion.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Monte-Carlo modeling of phonon thermal transport using DFT-based anisotropic dispersion relations over the full Brillouin zone

In this work, we present a Monte Carlo (MC) approach to solve the phonon Boltzmann transport equation (BTE) in which the anisotropic phonon dispersion relations over the full Brillouin zone (BZ) are used. In this approach, the discretization of the BZ used to compute the phonon relaxation time places constraints on the direction of scattered phonons in the real-space simulation domain. The phonon dispersion and phonon relaxation times are calculated using the density functional theory (DFT) approach. The modified MC approach is validated by a close examination of its ability to simulate phonon transport in both the ballistic and diffusive regimes for multiple materials including GaAs, InAs, ThO 2 , and α-U. In doing so, the phonon thermal conductivities from 100 K to 1000 K are calculated and compared with traditional non-transport solution of the phonon BTE. It is found that the phonon thermal conductivities of α-U and ThO 2 obtained from MC simulations using isotropic dispersion are larger than the values obtained using anisotropic phonon dispersion relations over the full BZ. The effect of phonon-defect scattering on the thermal conductivity of ThO 2 is also studied as an application of the current MC approach and found to agree with previously computed values in the literature. The MC solver developed here has been parallelized as a step to demonstrate its potential to solving computationally intensive phonon thermal transport problems at the mesoscale.

36 MATERIALS SCIENCE↗

A sensitivity analysis of twinning crystal plasticity finite element model using single crystal and poly crystal Zircaloy

The popularity of crystal plasticity finite element method (CPFEM) models is increasing due to their ability to predict the mechanical response of crystalline materials such as metals and metal alloys more accurately than traditional continuum mechanics models. This is since the crystal plasticity models consider the effect of atomic structure, microstructural morphology, and properties of individual grains. These CPFEM models use a large number of material parameters in order to capture the mesoscale physics which comes with the downside of the tedious calibration process. In this paper, a CPFEM code was developed to include the twinning induced grain reorientation and subsequent crystallographic slip for HPC material. The developed code is incorporated in a large-scale, parallelized nonlinear solver WARP3D. Further, a sensitivity analysis with respect to 22 material parameters was then conducted using single crystal and polycrystal representative volume element (RVE) of Zircaloy material. Loading was applied along five different crystallographic orientations for single crystal RVE and along three directions namely, rolling (RD), transverse (TD), and normal (ND) direction for polycrystal RVE. Results obtained from the sensitivity analysis were used for the calibration of material parameters for Zircaloy. Finally, developed code along with calibrated material parameters was used to investigate the effect of the hydride phase formation in Zircaloy which is a typical case observed for nuclear applications. It was found that the volume fraction of the hydride phase has a significant impact on the mechanical properties of Zircaloy.

36 MATERIALS SCIENCE↗

Computational modeling of grain boundary segregation: A review

Nearly all metals, alloys, ceramics, and their associated composites are polycrystalline in nature, with grain boundaries that separate well-defined crystalline regions that influence materials properties. In all but the most pure elemental systems, intentional solutes or impurities are present and can segregate to, or less commonly away from, the grain boundaries, in turn influencing boundary behavior, their stability, and associated materials properties. In some cases, grain-boundary segregation can also trigger “phase-like” structural transitions that dramatically alter the essential nature of the boundary. With the development of advanced electron microscopy techniques, researchers can directly observe grain-boundary structures and segregation with atomic precision. Despite such spatial resolution, the underlying mechanisms governing grain-boundary segregation remain difficult to characterize. As a result, computational modeling techniques such as density functional theory, molecular dynamics, mesoscale phase-field, continuum defect theory, and others are important complementary tools to experimental observations for studying grain-boundary segregation behavior. In conclusion, these computational methods offer the ability to explore the underlying formation mechanisms of grain-boundary segregation, elucidate complex segregation behavior, and provide insights into solutions to effectively controlling microstructure.

36 MATERIALS SCIENCE↗

Phase-field modeling of radiation-induced composition redistribution: An application to additively manufactured austenitic Fe–Cr–Ni

Multicomponent alloys undergoing irradiation damage develop radiation-induced composition redistribution at point defect sinks such as grain boundaries (GBs) and dislocations. Such redistribution results in undesired changes to their mechanical behavior and corrosion resistance. Additively manufactured alloys proposed for future nuclear applications are expected to demonstrate a distinct response to irradiation owing to their unique microstructure with as-solidified dislocation density and chemical microsegregation. To capture the composition redistribution in such systems, we develop a mesoscale model with coupled evolution of atomic and point defect components in the presence of dislocation density, dislocation heterogeneity, and thermodynamic interactions at the GB. The model is parameterized for an FCC Fe–Cr–Ni alloy as a representative system for austenitic stainless steels, and simulations are performed in 1D and 2D as a function of irradiation temperature, dose, dislocation density, and grain size. Radiation-induced segregation (RIS) characterized by Cr depletion and Ni enrichment is predicted at both the GB and the dislocation cell wall, with RIS being lower in magnitude but wider at the cell wall. Strongly biased absorption of self-interstitials by dislocations is found to suppress Ni enrichment but slightly enhance Cr depletion under certain conditions. Thermodynamic segregation at the GB is predicted to be narrower and opposite in sign to RIS for both Cr and Ni. Importantly, non-monotonic segregation is found to occur when both thermodynamic and RIS mechanisms are considered, providing a novel physical interpretation of experimental observations. The model is expected to serve as a key tool in accelerated qualification of irradiated materials.

additively manufactured microstructure↗

Machine-learning-assisted deciphering of microstructural effects on ionic transport in composite materials: A case study of Li 7 La 3 Zr 2 O 12 -LiCoO 2

The effective diffusivity of ionic species in multiphase materials is critical for the design and function of composite materials for electrochemical energy storage. In practice, effective diffusivity depends sensitively not only on the intrinsic diffusivities of constituting materials but also on their topological arrangement; nevertheless, these coupled contributions are oversimplified in most analytical models. Here, we combine atomistically informed mesoscale modeling and machine learning (ML) analysis to unravel how such features affect effective diffusivity in two-phase composites. Using the Li 7 La 3 Zr 2 O 12 -LiCoO 2 composite solid-state battery cathode as a model system, we compute effective diffusivity for 600 distinct dense polycrystalline microstructures with different topological configurations of grains, grain boundaries, and heterointerfaces. We verify that in addition to atomic-scale variabilities, microstructural feature diversity can significantly impact effective transport properties. Across the ensemble of test microstructures, this often results in bimodal distributions of effective diffusivity that encompass two qualitatively distinct operating mechanisms, which we identify via flux analysis. An ML approach reveals that the most critical determining factors for effective diffusivity are the connectivity of bulk phases and their heterointerfaces. The role of ionic mobility at the heterointerfaces is also discussed. These insights highlight the combined importance of microstructure and interface engineering in tuning the transport properties of ionic species in composite materials. In conclusion, our framework can also be extended for understanding generic microstructure-property relationships in other complex multiphase materials.

25 ENERGY STORAGE↗

Modeling and Experimental Validation of a Direct-Contact Counter-Flow Fluidized Bed Heat Exchanger for Thermal Energy Storage (TES) Applications

Particle-based thermal energy storage (TES) systems are an emerging energy storage technology. The technological advances have reduced costs, making TES more competitive and reliable in the marketplace but an efficient and reliable operation is heavily dependent on coherent heat transfer between air to particles or vice versa. The particle-based TES technologies provide an intermediate system that can store energy for short (0-10 h), long (10-200 h) and seasonal (> 200 h) timescales. The TES systems store energy by converting electricity to thermal energy; electricity can be directly sourced intermittent generation technologies and/or the grid, helping manage peak loads and other mismatches in supply and demand. The overall efficiency of the TES system depends on the performance of system components (particle storage silos and particle transfer mechanism etc.). The particle heat exchanger is one of the key system components that affects the system efficiency. The pressurized fluidized bed heat exchanger (PFB HX) performance is challenging to predict due to the chaotic behavior of particle and fluid interaction. This research presents a computational study of a novel direct-contact, counter-flow and air-to-particles PFB HX, that contributes in advancing the particle-based long-duration TES technologies. For the current analysis an unsteady Eulerian-Eulerian CFD model was developed and validated against experiments performed at the National Laboratory of the Rockies for two particle sizes (600 ..mu..m and 825 ..mu..m ). Following validation, parametric simulations were conducted to evaluate the effects of interphase drag models (Syamlal-O'Brien and Gidaspow), particle size, bed height and the influence of a frictional-viscosity term on hydrodynamics and heat transfer between the air & particles. The key findings from the analysis are: (1) for the studied operating window Syamlal-O'Brien provides superior agreement with measured gas temperatures (errors generally < 10%) while Gidaspow shows large deviations for the coarse particle case; (2) model predictions are most sensitive in the lower 0.2 m above the air distributor where bubble initiation and local mixing dominate interphase heat transfer; (3) representation of the distributor (number of inlet ports) materially affects predicted local mixing and temperature stratification; and (4) the Eulerian-Eulerian framework reproduces bulk thermal trends but shows regime dependent limitations for coarse particles, motivating mesoscale informed closures for scale-up analysis for future studies. These results provide validated guidance for drag selection and distributor design in particle-based thermal energy storage applications. Collectively, the validated model and parametric results quantify key drivers of PHB-HX performance and provide practical guidance for design and optimization. The results provide confidence in the model predictability and provide a step forward to improve on heat exchange performance. The demonstrated performance and modeling approach support the deployment and further development of this novel PHB-HX concept for robust, particle-based long-duration thermal energy storage systems.

25 ENERGY STORAGE↗

Microstructure-Sensitive modeling of surface roughness and notch effects on extreme value fatigue response

Fatigue life knockdown due to surface roughness induces local stress concentrations. Physical fatigue experiments are crucial to qualify materials for use but are limited in their capability to examine the mechanisms underlying the fatigue crack formation driving force due to these stress concentrations. Microstructure-sensitive computational models can aid in this regard and are employed in this work to examine the effects of surface roughness in Al 7075-T6. Digital microstructure models are subjected to crystal plasticity finite element method (CPFEM) simulations to examine the combined effects of microstructure and grain scale asperities on the driving force for fatigue crack formation. Following elastic–plastic shakedown, mesoscale volume-averaged fatigue indicator parameters (FIPs) are computed within fatigue damage process zones of grains. An increase in the intensity of realistic surface roughness profiles corresponds accordingly to larger FIPs but it is difficult to precisely interpret the mechanisms and effects of this intensification. We thus parametrically investigate surface asperities that couple with microstructure and lack of constraint on slip at the free surface to describe microstructure-sensitive surface roughness knockdown effects on fatigue resistance. The depth of a single valley (i.e., notch) due to surface roughness is found to be more detrimental to fatigue resistance than the radius (i.e., notch width, not to be confused with notch acuity). Finally, we also find that different crystallographic textures, grain morphologies, and number of sampled grains have a limited effect on the depth of influence of these notches.

36 MATERIALS SCIENCE↗

A predictive modeling tool for damage analysis and design of hydrogen storage composite pressure vessels

In this study, a predictive modeling tool is developed for damage analysis and design of hydrogen (H 2 ) storage composite pressure vessels. It integrates micromechanics of matrix cracking into a continuum damage mechanics (CDM) description for damage evolution, and three-dimensional (3D) finite element (FE) modeling of the vessel structural response. At the scale of the composite layer (mesoscale), the temperature-dependent stiffness reduction law in terms of the damage variable for transverse matrix cracking is computed using an Eshelby-Mori-Tanaka approach for the initial composite thermoelastic properties and a self-consistent model for the stiffness reduction as a function of the damage variable. While transverse matrix cracking obeying a damage evolution relation can progressively evolve from an initiation to a saturation state, fiber failure is predicted by a micromechanical fiber rupture criterion that accounts for the fiber strength and matrix stress. The implementation of this integrated multiscale modeling model into a 3D FE formulation enables damage analysis and design of H 2 storage composite pressure vessels. The developed tool is illustrated through 3D damage analyses of a cryogenically compressed H 2 storage vessel model subjected to thermomechanical loadings to investigate effects of the helical layer fiber orientation and loading scenario on damage development, vessel integrity and burst pressure.

08 HYDROGEN↗