Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuum scale modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11

Integrated Modeling of Aerosol, Cloud, Precipitation and Land Processes at Satellite-Resolved Scales

In this talk, I will present recent results from a project led at NASA/GSFC, in collaboration with NASA/MSFC and JHU, focused on the development and application of an observation-driven integrated modeling system that represents aerosol, cloud, precipitation and land processes at satellite-resolved scales. The project, known as the NASA Unified WRF (NU-WRF), is funded by NASA's Modeling and Analysis Program, and leverages prior investments from the Air Force Weather Agency and NASA's Earth Science Technology Office (ESTO). We define "satellite-resolved" scales as being within a typical mesoscale atmospheric modeling grid (roughly 1-25 km), although this work is designed to bridge the continuum between local (microscale), regional (mesoscale) and global (synoptic) processes. NU-WRF is a superset of the standard NCAR Advanced Research WRF model, achieved by fully integrating the GSFC Land Information System (LIS, already coupled to WRF), the WRF/Chem enabled version of the Goddard Chemistry Aerosols Radiation Transport (GOCART) model, the Goddard Satellite Data Simulation Unit (SDSU), and boundary/initial condition preprocessors for MERRA and GEOS-5 into a single software release (with source code available by agreement with NASA/GSFC). I will show examples where the full coupling between aerosol, cloud, precipitation and land processes is critical for predicting local, regional, and global water and energy cycles, including some high-impact phenomena such as floods, hurricanes, mesoscale convective systems, droughts, and monsoons.

Peters-Lidard, Christa↗

A Micromorphic Length-Scale Coupling Framework for the Determination of Higher-Order Constitutive Models and the Multi-Scale Simulation of Heterogeneous Materials [Thesis]

Heterogeneous materials and materials with complex microstructures pose a unique challenge in the development of accurate models of their response to external stimuli. These difficulties principally arise due to the difficulty in characterizing and modeling the constituents and their interactions. It is usually possible, though non-trivial, to construct an explicit representation of the microstructure (a direct numeric simulation or DNS) but the method by which the complex modes of deformation and other processes can be homogenized to a reduced-order approximation is frequently unclear. Many homogenization approaches are ad-hoc and lack a strong justification beyond ease of computation. Furthermore, the homogenization approach can, in some cases, not utilize the full breadth of information available in the computation of the macro-scale stresses and deformation measures. It is also noteworthy that homogenization, by its very nature, will tend to obfuscate details of processes occurring at the lower length-scale. It is therefore of interest to include as much information as possible in the construction of the reduced-order model so as to be predictive in a variety of loading environments. We here present a length-scale bridging technique based upon the micromorphic continuum mechanics of Eringen which incorporates volume and surface area averages as a part of its construction. This approach enforces the balance equations at the micro-scale and then studies the effect of the spatially varying nature on the macro-scale. This leads, naturally, to further balance equations which are solved at the macro-scale. This work details a homogenization framework which arises naturally from the micromorphic construction of Eringen attempting to introduce no definitions beyond which are justifiable from micro-structural considerations. One of the results of this effort is the so-called “micromorphic filter” which has been developed and applied to several DNS to demonstrate its effectively. In order to determine the macroscopic degrees of freedom we utilize the special case of an overlap coupling technique where the macro-scale is fully constrained to the micro-scale. This enables us to study the resulting material properties as expected but also allows us to further study the boundary conditions on the additional degrees of freedom at the macro-scale.

42 ENGINEERING↗

From Pseudostreamer Jets to Coronal Mass Ejections:: Observations of the Breakout Continuum

The magnetic breakout model, in which reconnection in the corona leads to destabilization of a filament channel, explains numerous features of eruptive solar events, from small-scale jets to global-scale coronal mass ejections (CMEs). The underlying multipolar topology, pre-eruption activities, and sequence of magnetic-reconnection onsets (first breakout, then flare) of many observed fast CMEs/eruptive flares are fully consistent with the model. Recently, we have demonstrated that most observed coronal-hole jets in fan/spine topologies also are induced by breakout reconnection at the null point above a filament channel (with or without a filament).For these two types of eruptions occurring in similar topologies, the key question is, why do some events generate jets while others form CMEs? We focused on the initiation of eruptions in large bright points/small active regions that were located in coronal holes and clearly exhibited null-point (fan/spine) topologies: such configurations are referred to as pseudostreamers. We analyzed and compared SDO/AIA,SOHO/LASCO, and RHESSI observations of three events. Our analysis of the events revealed two new observable signatures of breakout reconnection prior to the explosive jet/CME outflows and flare onset: coronal dimming and the opening-up of field lines above the breakout current sheet. Most key properties were similar among the selected erupting structures, thereby eliminating region size, photospheric field strength, magnetic configuration, and pre-eruptive evolution as discriminating factors between jets and CMEs. We consider the factors that contribute to the different types of dynamic behavior, and conclude that the determining factor is the ratio of the magnetic free energy associated with the filament channel compared to the energy associated with the overlying flux inside and outside the pseudostreamer dome. Our results imply that a pseudostreamer eruption must be fast in order to produce an interplanetary CME.

Sun↗

Magnetism in Moiré Materials

This project focused on the theory of the electronic properties of magnetic moire materials. Moire materials can be formed by overlaying two or more two-dimensional crystals that are either semiconductors or semimetals and establishing a moire pattern. Low-energy electronic properties are then described by a continuum model with crystalline periodicity. Moire materials are therefore artificial two-dimensional crys- tals whose periodicity is at the tens of nanometers scale, increasing the unit cell area by a factor of about 10,000 compared to real crystals. They are important because the main tuning knob of chemistry, the number of electrons per effective atom or molecule (the analog of the number of electrons per unit cell in real crystals), can be shifted by more than one simply by using electrical gates. The research covered two diffeerent classes of moire materials, graphene multilayers twisted to a rotation angle at which correlations are strong, and transition metal dichalcogenide semiconductors. The proposal identifed one project in each area both related very broadly to magnetism: i) Advancing understanding of broken spin-valley and sublattice symmetries in Magic Angle Twisted Bilayer Graphene and ii) Developing the theory of itinerant electron magnetism in Gamma-valley and K-valley transition metal dichalcogenide moire materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Models of the Cartwheel ring galaxy: Spokes and starbursts

Recent observations of this famous ring galaxy, including optical and near-infrared CCD surface photometry, and VLA radio continuum and 21 cm line mapping (Higdon 1992b, in prep.), have inspired a renewed modeling effort. Toomre's (1978, in The Large-scale Structure of the Universe, eds. Longair and Einasto) series of restricted three-body simulations demonstrated how the multiple rings could be produced in a nearly head-on galaxy collision. New models with a halo-dominated potential based on the 21 cm rotation curve are able to reproduce such details as the spacing between rings, ring widths, offset of the nucleus, and several kinematical features, thus providing strong support for the collisional theory. The new observations have shown there are little or no old stars in Cartwheel; it may consist almost entirely of gas and stars produced as a result of compression in the ring wave. To model this process Smooth Particle Hydrodynamics (SPH) simulations of the Cartwheel disk have been performed. Fixed gravitational potentials were used to represent the Cartwheel and a roughly 30 percent mass collision partner. The interaction dynamics was treated as in the usual restricted three-body approximation, and the effects of local self-gravity between disk particles were calculated. We are particularly interested in testing the theory that enhanced star formation in waves is the result of gravitational instability in the compressed region (see e.g. Kennicutt 1989, ApJ 344, 685). The gas surface density in a number of simulations was initialized to a value slightly below the threshold for local gravitational instability throughout most of the disk. The first ring wave produces relatively modest compressions (a factor of order a few), triggering instability in a narrow range of wavelengths. Self-gravity in the disk is calculated over a comparable range of scales. Simulations were run with isothermal, adiabatic, and adiabatic with radiative cooling characterized by a relatively short timescale. The isothermal approximation is good except in the vicinity of the strong second (inner) ring, and several snapshots from one case are shown in the figure below. Flocculent spiral segments are present before the collision, and these are compressed into dense knots in the ring wave. These knots are likely to be sites of vigorous star formation. In the strong rarefaction behind the outer ring most of the knots are radially stretched and sheared, giving rise to spoke-like features. A few dense knots are evidently very tightly bound, because they retain their coherence and are stretched relatively little through the rarefaction. This is in accord with evidence for continuing star formation in some spokes (Marcum, Appleton and Higdon 1992). The number and spacing of spokes is a direct function of the scale of the gravitational instability in the disk. Thus, the gravitational instability theory, together with the hypothesis that massive stars are only formed in dense knots of gas, can account for most of the distinct morphology of the Cartwheel.

Struck-Marcell, Curtis↗

Evaluation of Nuclear Spent Fuel Disposal in Clay-Bearing Rock - Process Model Development and Experimental Studies

The DOE R&D program under the Spent Fuel Waste Science Technology (SFWST) campaign has made key progress in modeling and experimental approaches towards the characterization of chemical and physical phenomena that could impact the long-term safety assessment of heat-generating nuclear waste disposition in deep clay/shale/argillaceous rock. International collaboration activities such as heater tests and postmortem analysis of samples recovered from these have elucidated key information regarding changes in the engineered barrier system (EBS) material exposed to years of thermal loads. Chemical and structural analyses of sampled bentonite material from such tests has as well as experiments conducted on these are key to the characterization of thermal effects affecting bentonite clay barrier performance and the extent of sacrificial zones in the EBS during the thermal period. Thermal, hydrologic, and chemical data collected from heater tests and laboratory experiments has been used in the development, validation, and calibration of THMC simulators to model near-field coupled processes. This information leads to the development of simulation approaches (e.g., continuum vs. discrete) to tackle issues related to flow and transport at various scales of the host-rock and EBS design concept. Consideration of direct disposal of large capacity dual-purpose canisters (DPCs) as part of the back-end SNF waste disposition strategy has generated interest in improving our understanding of the effects of elevated temperatures on the EBS design. This is particularly important for backfilled repository concepts where temperature plays a key role in the EBS behavior and long-term performance. This report describes multiple R&D efforts on disposal in argillaceous geologic media through development and application of coupled THMC process models, experimental studies on clay/metal/cement barrier and host-rock (argillite) material interactions, molecular dynamic (MD) simulations of water transport during (swelling) clay dehydration, first-principles studies of metaschoepite (UO 2 corrosion product) stability, and advances in thermodynamic plus surface complexation database development. Drift-scale URL experiments provides key data for testing hydrological-chemical (HC) model involving strong couplings of fluid mixing and barrier material chemical interactions. The THM modeling focuses on heater test experiments in argillite rock and gas migration in bentonite as part of international collaboration activities at underground research laboratories (URLs). In addition, field testing at an URL involves in situ analysis of fault slip behavior and fault permeability. Pore-scale modeling of gas bubble migration is also being investigated within the gas migration modeling effort. Interaction experiments on bentonite samples from heater test under ambient and elevated temperatures permit the evaluation of ion exchange, phase stability, and mineral transformation changes that could impact clay swelling. Advances in the development, testing, and implementation of a spent nuclear fuel (SNF) degradation model coupled with canister corrosion focus on the effects of hydrogen gas generation and its integration with Geologic Disposal Safety Assessment (GDSA). GDSA integration activities includes evaluation of groundwater chemistries in shale formations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

The temperature scale of hot DA white dwarfs - Temperatures from far-ultraviolet continuum fluxes

Observed far-UV to visible flux ratios of all DA white dwarfs hotter than about 25,000 K are compared with ratios predicted by models in order to derive effective temperatures for the stars. All suitable IUE spectra of seven hot DA white dwarfs for which accurate temperatures based on hydrogen line profiles were available are used to derive a flux correction to the 1980 IUE absolute calibration. Application of the corrections to the IUE spectra results in derived temperature which are consistent with the line profile temperatures.

Finley, David S.↗

The Spectral Properties of Gamma-ray Bursts: a Review of Recent Development

Developments in the spectroscopy of gamma ray bursts (GRB) are reviewed. The general question of the validity of the spectral results, particularly with regard to features in the spectrum, is discussed. Confirmations of these spectral features are summarized. Results from the KONUS experiments on Venera 13 and 14 are reviewed. The status of models of the continuum spectrum is summarized. A number of different radiation mechanisms appear capable of fitting the data. These include thermal bremsstrahlung, thermal synchrotron and inverse Compton. Rapid variability of the spectra shape on time scales 76] 0.25 sec. was reported. The characteristic energy of the spectrum was observed to vary over nearly an order of magnitude during individual events. A strong correlation between spectral hardness and luminosity was found. Low-energy (50 keV) absorption features and high-energy (400 keV) emission features continue to appear in GRB spectra. Understanding the origin of these lines in the context of the existing continuum models remains a difficult problem.

Teegarden, B. J.↗

A stochastic solver based on the residence time algorithm for crystal plasticity models

Abstract The deformation of crystalline materials by dislocation motion takes place in discrete amounts determined by the Burgers vector. Dislocations may move individually or in bundles, potentially giving rise to intermittent slip. This confers plastic deformation with a certain degree of variability that can be interpreted as being caused by stochastic fluctuations in dislocation behavior. However, crystal plasticity (CP) models are almost always formulated in a continuum sense, assuming that fluctuations average out over large material volumes and/or cancel out due to multi-slip contributions. Nevertheless, plastic fluctuations are known to be important in confined volumes at or below the micron scale, at high temperatures, and under low strain rate/stress deformation conditions. Here, we develop a stochastic solver for CP models based on the residence-time algorithm that naturally captures plastic fluctuations by sampling among the set of active slip systems in the crystal. The method solves the evolution equations of explicit CP formulations, which are recast as stochastic ordinary differential equations and integrated discretely in time. The stochastic CP model is numerically stable by design and naturally breaks the symmetry of plastic slip by sampling among the active plastic shear rates with the correct probability. This can lead to phenomena such as intermittent slip or plastic localization without adding external symmetry-breaking operations to the model. The method is applied to body-centered cubic tungsten single crystals under a variety of temperatures, loading orientations, and imposed strain rates.

42 ENGINEERING↗

Updates on the Predictive Materials Modeling Software Tools

Updates on NASA‘s efforts to build a Predictive Material Modeling (PMM) framework from the micro-scale to the macro-scale are presented in this abstract. The PMM effort is part of the Entry Systems Modeling (ESM) project under NASA’s Game Changing Development (GCD) program. To reduce the need for extensive testing and accelerate the design cycle process, ESM is developing simulation and modeling tools that enable the characterization of the properties of thermal protection materials and their response to extremely hot plasma. The Porous Microstructure Analysis (PuMA) software has been developed to compute effective material properties and perform material response simulations on digitized microstructures of porous media. PuMA is able to import three-dimensional digital images obtained from X-ray microtomography or to generate artificial microstructures that mimic real materials. PuMA also provides a module for interactive 3D visualizations. Version 3, which was recently released as open-source, includes modules to compute simple morphological properties such as porosity, volume fractions, pore diameter, and specific surface area. Additional capabilities include the determination of effective thermal and electrical conductivity (both radiative and solid conduction - including the ability to simulate local anisotropy for the latter); effective diffusivity and tortuosity from the continuum to the rarefied regime; techniques to determine the local material orientation, as well as mechanical properties (elasticity coefficients), and permeability. Computed properties are then used to inform a macro-scale material response model, such as those implemented in the Porous material Analysis Toolbox based on OpenFOAM (PATO) software developed within ESM. The computational model in PATO is a generic heat and mass transfer model for porous reactive materials containing several solid phases and a single gas phase. The detailed chemical interactions occurring between the solid phases and the gas phase are modeled at the pore scale, assuming Local Thermal Equilibrium. Recent efforts include the development of a mechanical erosion model as well as a unified model allowing an intrinsic coupling between fluid and material. Comparison to flight data (Mars Science Laboratory [MSL] Entry Descent and Landing Instrument [MEDLI] and Mars 2020 MEDLI2) is critical in order to validate these computational tools. Examples of ablative material response using the code will be presented, including 3D simulations of the full-scale heatshield of the MSL capsule. The simulations demonstrated the ability of the modern material response code, PATO, to handle the material response of geometrically complex and large domains through the use of massively parallel computations.

material modeling↗

Asymptotic Freedom at the Berezinskii-Kosterlitz-Thouless Transition without Fine-Tuning Using a Qubit Regularization

We propose a two-dimensional hard-core loop-gas model as a way to regularize the asymptotically free massive continuum quantum field theory that emerges at the Berezinskii-Kosterlitz-Thouless transition. Without fine-tuning, our model can reproduce the universal step-scaling function of the classical lattice X Y model in the massive phase as we approach the phase transition. This is achieved by lowering the fugacity of Fock-vacuum sites in the loop-gas configuration space to zero in the thermodynamic limit. Some of the universal quantities at the Berezinskii-Kosterlitz-Thouless transition show smaller finite size effects in our model as compared to the traditional X Y model. Our model is a prime example of qubit regularization of an asymptotically free massive quantum field theory in Euclidean space-time and helps understand how asymptotic freedom can arise as a relevant perturbation at a decoupled fixed point without fine-tuning. Published by the American Physical Society 2024

Physics↗

Evaluating Basin-Scale Forest Adaptation Scenarios: Wildfire, Streamflow, Biomass, and Economic Recovery Synergies and Trade-Offs

Active forest management is applied in many parts of the western United States to reduce wildfire severity, mitigate vulnerability to drought and bark beetle mortality, and more recently, to increase snow retention and late-season streamflow. A rapidly warming climate accelerates the need for these restorative treatments, but the treatment priority among forest patches varies considerably. We simulated four treatment scenarios across the 3,450 km 2 Wenatchee River basin in eastern Washington, United States. We used a decision support tool (DST) to assess trade-offs and synergies within and among treatments on wildfire risk and smoke emissions, water yield and snow retention, biomass production, and economic return. Treatment scenarios emphasized prescribed burning ( BurnOnly ), biomass production ( MaxBiomass ), gap-based thinning to optimize water yield ( IdealWater ), and a principle-based restoration scenario ( RA1 ). Fire hazard, smoke emissions, and biomass production metrics were evaluated across scenarios using the Forest Vegetation Simulator, and water yields were modeled using the Distributed Hydrology Soil Vegetation Model. Simulations were summarized to both patch- (10 1 –10 2 ha) and subwatershed- (10 3 –10 4 ha) scales, and treatment effects were evaluated against an untreated baseline landscape. We used logic models to rank effect sizes by scenario across metrics along a continuum between −1 (no or weak effect) to +1 (large effect). All treatments produced benefits across one or more ecosystem services and led to synergistic benefits to water yield and wildfire hazard reduction. Tradeoffs among resource benefits were clear in wilderness where reliance on prescribed burning without mechanical treatment increased costs and eliminated the potential for biomass recovery. The BurnOnly scenario improved fire risk metrics and streamflow, but effect sizes were lower compared to other treatments. IdealWater showed the strongest benefits overall, demonstrating the ability to capture multiple resource benefits through spatially explicit thinning. Our study provides a framework for integrating strategic and tactical models that evaluate tradeoffs and synergies gained through varied management approaches. We demonstrate the utility of decision support modeling to enhance management synergies across large landscapes.

54 ENVIRONMENTAL SCIENCES↗

Computational Modeling as a Design Tool in Microelectronics Manufacturing

Plans to introduce pilot lines or fabs for 300 mm processing are in progress. The IC technology is simultaneously moving towards 0.25/0.18 micron. The convergence of these two trends places unprecedented stringent demands on processes and equipments. More than ever, computational modeling is called upon to play a complementary role in equipment and process design. The pace in hardware/process development needs a matching pace in software development: an aggressive move towards developing "virtual reactors" is desirable and essential to reduce design cycle and costs. This goal has three elements: reactor scale model, feature level model, and database of physical/chemical properties. With these elements coupled, the complete model should function as a design aid in a CAD environment. This talk would aim at the description of various elements. At the reactor level, continuum, DSMC(or particle) and hybrid models will be discussed and compared using examples of plasma and thermal process simulations. In microtopography evolution, approaches such as level set methods compete with conventional geometric models. Regardless of the approach, the reliance on empricism is to be eliminated through coupling to reactor model and computational surface science. This coupling poses challenging issues of orders of magnitude variation in length and time scales. Finally, database development has fallen behind; current situation is rapidly aggravated by the ever newer chemistries emerging to meet process metrics. The virtual reactor would be a useless concept without an accompanying reliable database that consists of: thermal reaction pathways and rate constants, electron-molecule cross sections, thermochemical properties, transport properties, and finally, surface data on the interaction of radicals, atoms and ions with various surfaces. Large scale computational chemistry efforts are critical as experiments alone cannot meet database needs due to the difficulties associated with such controlled experiments and costs.

Meyyappan, Meyya↗

Multi-scale Rule-of-Mixtures Model of Carbon Nanotube/Carbon Fiber/Epoxy Lamina

A unidirectional carbon fiber/epoxy lamina in which the carbon fibers are coated with single-walled carbon nanotubes is modeled with a multi-scale method, the atomistically informed rule-of-mixtures. This multi-scale model is designed to include the effect of the carbon nanotubes on the constitutive properties of the lamina. It included concepts from the molecular dynamics/equivalent continuum methods, micromechanics, and the strength of materials. Within the model both the nanotube volume fraction and nanotube distribution were varied. It was found that for a lamina with 60% carbon fiber volume fraction, the Young's modulus in the fiber direction varied with changes in the nanotube distribution, from 138.8 to 140 GPa with nanotube volume fractions ranging from 0.0001 to 0.0125. The presence of nanotube near the surface of the carbon fiber is therefore expected to have a small, but positive, effect on the constitutive properties of the lamina.

Frankland, Sarah-Jane V.↗

Characterizing leaf-scale fluorescence with spectral invariants

Sun-induced chlorophyll fluorescence (SIF) is increasingly recognized as a non-destructive probe for tracking terrestrial photosynthesis. Emerging developments in spectral invariants theory provide an innovative and efficient approach for representing SIF radiative transfer processes at the canopy scale. However, modeling leaf-scale fluorescence based on the spectral invariants properties (SIP) remains underexplored. In this study, the spectral invariants theory is employed for the first time to model the leaf-scale total, backward and forward fluorescence (leaf-SIP SIF). The leaf-SIP SIF model separates the leaf-scale radiative transfer process into two distinct components: the wavelength-dependent one associated with leaf biochemical properties, and the wavelength-independent component linked to leaf structural characteristics. The leaf structure-related effects are characterized by two spectrally invariant parameters: the photon recollision probability (p) and the scattering asymmetry parameter (q), which are parameterized using the directly measurable leaf dry matter. Evaluation against field measurements shows that the proposed leaf-SIP SIF model has a good performance, with coefficient of determination (R 2 ) of 0.89, 0.89, 0.90 and root mean squared errors (RMSE) of 1.28, 0.69, 0.74 Wm -2 µm -1 sr -1 , respectively for the total, backward, and forward fluorescence (660–800 nm). The leaf-SIP SIF model with a more concise formulation demonstrates comparable performance with the widely used Fluspect model. Further, the leaf-SIP SIF model provides a simple and efficient approach for simulating leaf-scale fluorescence, with the potential to be integrated into a unified SIP-based model framework for simulating the radiative transfer processes across the soil-leaf-canopy-atmosphere continuum.

59 BASIC BIOLOGICAL SCIENCES↗

Microemulsions in the driven Widom-Rowlinson lattice gas

An investigation of the two-dimensional Widom-Rowlinson lattice gas under an applied drive uncovered a remarkable nonequilibrium steady state in which uniform stripes (reminiscent of an equilibrium lamellar phase) form perpendicular to the drive direction [R. Dickman and R. K. P. Zia, Phys. Rev. E 97, 062126 (2018)]. Here we study this model at low particle densities in two and three dimensions, where we find a disordered phase with a characteristic length scale (a “microemulsion”) along the drive direction. Further we develop a continuum theory of this disordered phase to derive a coarse-grained field-theoretic action for the nonequilibrium dynamics. The action has the form of two coupled driven diffusive systems with different characteristic velocities, generated by an interplay between the particle repulsion and the drive. We then show how fluctuation corrections in the field theory may generate the characteristic features of the microemulsion phase, including a peak in the static structure factor corresponding to the characteristic length scale. This work lays the foundation for understanding the stripe phenomenon more generally.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Molecular scale insight of pore morphology relation with mechanical properties of amorphous silica using ReaxFF

Porous materials are typically heterogeneous and they contain large variations of micro-/nano-pore structures, causing complicated behaviors. In continuum models, most mechanical properties of porous materials are estimated based on porosity, while the variations of micro/nano structures are ignored. That could be problematic as the microscopic heterogeneity may affect the mechanical response of porous materials. Thus, understanding micro/nano heterogeneity impact has been the focus in many scientific and engineering subjects. In the present study, we investigated the effect of nanopore structure (including pore shape and orientation) as well as porosity on mechanical properties of amorphous silica (a-SiO 2 ). The pore sizes in our simulations are comparable to the corresponding ones observed in a-SiO 2 based materials. We found that the existing of nanopores strongly influences Young’s modulus (E) and critical energy release rate (G IC ). These properties decrease with increasing porosity. Importantly, the impact of nanopores was characterized by structural parameters of porous materials. In addition to dependency on porosity, Young’s modulus also was found to vary as a function of potential energy per atom, which highly depends on nanopore shape. Furthermore, critical energy release rate was found to increase with increasing ligament length (also known as pore wall thickness). Finally, the results highlighted the importance of nanopore structures, which must be taken into account when studying fracture mechanisms in porous materials. Based on our findings, it was proposed that mechanical properties of porous materials can be controlled by nano-engineering pore structures.

36 MATERIALS SCIENCE↗

Multiscale modeling of microwave-heated multiphase systems

Modeling of microwave heating of multiphase systems requires resolving the electromagnetic field from the single particle to the entire cavity. We introduce a multiscale methodology for computationally affordable simulations of microwave heating of multiphase systems consisting of one phase dispersed in a continuum phase. The methodology homogenizes the original multiphase system, by taking advantage of the large separation of length scales, to calculate the effective permittivity, effective thermal conductivity, and volumetric power absorbed in each phase. The methodology is rigorously assessed against particle-resolved detailed numerical simulations of model systems. We demonstrate a significant reduction in computational cost while retaining the accuracy of detailed simulations. Lastly, the present methodology enables high throughput exploration of novel multiphase system designs utilizing microwave heating.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗