Search NASA⌕ Search

SEARCH · Search NASA

Results for “macroscale”

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 487 records · Page 27

Knudsen heat capacity

We present a “Knudsen heat capacity” as a more appropriate and useful fluid property in micro/nanoscale gas systems than the constant pressure heat capacity. At these scales, different fluid processes come to the fore that are not normally observed at the macroscale. For thermodynamic analyses that include these Knudsen processes, using the Knudsen heat capacity can be more effective and physical. We calculate this heat capacity theoretically for non-ideal monatomic and diatomic gases, in particular, helium, nitrogen, and hydrogen. The quantum modification for para and ortho hydrogen is also considered. We numerically model the Knudsen heat capacity using molecular dynamics simulations for the considered gases, and compare these results with the theoretical ones.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Dark-to-arc transition in field emission dominated atmospheric microdischarges

We study the voltage-current characteristics of gas discharges driven by field emission of electrons at the microscale. Particle-in-cell with Monte Carlo collision calculations are first verified by comparison with breakdown voltage measurements and then used to investigate atmospheric discharges in nitrogen at gaps from 1 to 10 μm. The results indicate the absence of the classical glow discharge regime because field electron emission replaces secondary electron emission as the discharge sustaining mechanism. Additionally, the onset of arcing is significantly delayed due to rarefied effects in electron transport. While field emission reduces the breakdown voltage, the power required to sustain an arc of the same density in microgaps is as much as 30% higher than at macroscale.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Shock interaction with deformable particles using a constrained interface reinitialization scheme

In this paper, we present axisymmetric numerical simulations of shock propagation in nitromethane over an aluminum particle for post-shock pressures up to 10 GPa. We use the Mie-Gruneisen equation of state to describe both the medium and the particle. The numerical method is a finite-volume based solver on a Cartesian grid, that allows for multi-material interfaces and shocks, and uses a novel constrained reinitialization scheme to precisely preserve particle mass and volume. We compute the unsteady inviscid drag coefficient as a function of time, and show that when normalized by post-shock conditions, the maximum drag coefficient decreases with increasing post-shock pressure. We also compute the mass-averaged particle pressure and show that the observed oscillations inside the particle are on the particle-acoustic time scale. Finally, we present simplified point-particle models that can be used for macroscale simulations. In the Appendix, we extend the isothermal or isentropic assumption concerning the point-force models to non-ideal equations of state, thus justifying their use for the current problem.

Sridharan, P. (ORCID:0000000287902423)↗

Inter-phase heat transfer and energy coupling in turbulent dispersed multiphase flows

Here, the present paper addresses important fundamental issues of inter-phase heat transfer and energy coupling in turbulent dispersed multiphase flows through scaling analysis. In typical point-particle or two-fluid approaches, the fluid motion and convective heat transfer at the particle scale are not resolved and the momentum and energy coupling between fluid and particles are provided by proper closure models. By examining the kinetic energy transfer due to the coupling forces from the macroscale to microscale fluid motion, closure models are obtained for the contributions of the coupling forces to the energy coupling. Due to the inviscid origin of the added-mass force, its contribution to the microscale kinetic energy does not contribute to dissipative transfer to fluid internal energy as was done by the quasi-steady force. Time scale analysis shows that when the particle is larger than a critical diameter, the diffusive-unsteady kernel decays at a time scale that is smaller than the Kolmogorov time scale. As a result, the computationally costly Basset-like integral form of diffusive-unsteady heat transfer can be simplified to a non-integral form. Conventionally, the fluid-to-particle volumetric heat capacity ratio is used to evaluate the relative importance of the unsteady heat transfer to the energy balance of the particles. Therefore, for gas-particle flows, where the fluid-to-particle volumetric heat capacity ratio is small, unsteady heat transfer is usually ignored. However, the present scaling analysis shows that for small fluid-to-particle volumetric heat capacity ratio, the importance of the unsteady heat transfer actually depends on the ratio between the particle size and the Kolmogorov scale. Furthermore, the particle mass loading multiplied by the heat capacity ratio is usually used to estimate the importance of the thermal two-way coupling effect. Through scaling argument, improved estimates are established for the energy coupling parameters of each energy exchange mechanism between the phases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Atomistic modeling of dropwise condensation

The basic aim of the atomistic modeling of condensation of water is to determine the size of the stable cluster and connect phenomena occurring at atomic scale to the macroscale. In this paper, a population balance model is described in terms of the rate equations to obtain the number density distribution of the resulting clusters. The residence time is taken to be large enough so that sufficient time is available for all the adatoms existing in vapor-phase to loose their latent heat and get condensed. The simulation assumes clusters of a given size to be formed from clusters of smaller sizes, but not by the disintegration of the larger clusters. The largest stable cluster size in the number density distribution is taken to be representative of the minimum drop radius formed in a dropwise condensation process. A numerical confirmation of this result against predictions based on a thermodynamic model has been obtained. Results show that the number density distribution is sensitive to the surface diffusion coefficient and the rate of vapor flux impinging on the substrate. The minimum drop radius increases with the diffusion coefficient and the impinging vapor flux; however, the dependence is weak. The minimum drop radius predicted from thermodynamic considerations matches the prediction of the cluster model, though the former does not take into account the effect of the surface properties on the nucleation phenomena. For a chemically passive surface, the diffusion coefficient and the residence time are dependent on the surface texture via the coefficient of friction. Thus, physical texturing provides a means of changing, within limits, the minimum drop radius. The study reveals that surface texturing at the scale of the minimum drop radius does not provide controllability of the macro-scale dropwise condensation at large timescales when a dynamic steady-state is reached.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Statistical mechanics of transport processes in active fluids: Equations of hydrodynamics

We present the equations of hydrodynamics including mass, linear momentum, angular momentum, and energy are derived by coarse-graining the microscopic equations of motion for systems consisting of rotary dumbbells driven by internal torques. In deriving the balance of linear momentum, we find that the symmetry of the stress tensor is broken due to the presence of non-zero torques on individual particles. The broken symmetry of the stress tensor induces internal spin in the fluid and leads us to consider the balance of internal angular momentum in addition to the usual moment of momentum. In the absence of spin, the moment of momentum is the same as the total angular momentum. In deriving the form of the balance of total angular momentum, we find the microscopic expressions for the couple stress tensor that drives the spin field. We show that the couple stress contains contributions from both intermolecular interactions and the active forces. The presence of spin leads to the idea of balance of moment of inertia due to the constant exchange of particles in a small neighborhood around a macroscopic point. We derive the associated balance of moment of inertia at the macroscale and identify the moment of inertia flux that induces its transport. Lastly, we obtain the balances of total and internal energy of the active fluid and identify the sources of heat and heat fluxes in the system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Shock interactions with heterogeneous energetic materials

The complex physical phenomenon of shock wave interaction with material heterogeneities has significant importance and nevertheless remains little understood. In many materials, the observed macroscale response to shock loading is governed by characteristics of the microstructure. Yet the majority of computational studies aimed at predicting phenomena affected by these processes, such as initiation and propagation of detonation waves in explosives, or shock propagation in geological materials, employ continuum material and reactive burn model treatment. In an effort to highlight the grain-scale processes that underlie the observable effects in an energetic system, a grain-scale model for hexanitrostilbene (HNS) has been developed. Measured microstructures were used to produce synthetic computational representations of the pore structure, and a density functional theory molecular dynamics (DFT-MD) derived equation of state (EOS) was used for the fully dense HNS matrix. The explicit inclusion of microstructure along with a fully-dense EOS resulted in close agreement with historical shock compression experiments. More recent experiments on dynamic reaction threshold were also reproduced by inclusion of a global kinetics model. The complete model was shown to reproduce accurately the expected response of this heterogeneous material to shock loading. Mesoscale simulations were shown to provide clear insight into the nature of threshold behavior, and are a way to understand complex physical phenomena.

36 MATERIALS SCIENCE↗

Hydrodynamic dispersion in thin channels with micro-structured porous walls

Flowand transport within porous- and microtextured-walled channels is relevant to a number of natural and industrial processes. Designing and optimizing the topology of the pores and/or microstructure to achieve target performance at the system scale (or macroscale) is still an open question. In this work, we study whether hydrodynamic dispersion in microfluidic channels with walls structured by obstacles can be modeled by dispersion in channels with porous walls described as continuous porous media of zero or finite permeability. We perform single phase microfluidic non-reactive flow experiments in channels embedded in micropatterns with different topologies. Specifically, we focus on transverse riblets and arrays of pillars as examples of impermeable and permeable obstructions, respectively. We compare the experimental results with three models: 3D pore-scale simulations which resolve the micropattern geometry explicitly and two upscaled models which treat the micropattern as a continuum of zero or finite permeability. This study demonstrates that polydimethylsiloxane micromodels with appropriately patterned surfaces can be successfully employed to validate various continuum-scale modeling approximations in different physical regimes, identified by the order of magnitude of the Peclet number and the obstruction permeability.

Ling, Bowen↗

On the dynamic tensile strength of an FCC metal

The tensile response of ductile, polycrystalline metals is often accompanied by the formation of pores within the material, which coalesce and fail a plane within a metal. This large deformation process is broadly progressive, with a physical path consisting of nucleation, growth, coalescence, and failure, which occur individually over short periods of time. Thus distinct micro-mechanisms operate, each influenced by microstructure, loading path, and loading profile, which remains a significant challenge to represent and predict numerically at the macroscale. In a previous study, the influence of loading path on damage evolution in high-purity tantalum has been presented; in this paper, complimentary measurements are made on a pure, FCC copper. Samples were shock loaded to three different peak stresses using both symmetric impact, and two different composite flyer plate configurations, such that upon unloading, the three samples were subject to nearly identical tensile loading. The damage evolution in the soft-recovered copper samples was then quantified using optical metallography and electron-back-scatter diffraction. We shall compare metallurgical observations, velocimetry histories and simulations to discuss the dynamic failure mechanics observed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Investigating deformation and mesoscale void creation in HMX based composites using tomography based grain scale FEM

The microstructure of plastic bonded explosives (PBXs) significantly affects their macroscale mechanical characteristics. Imaging and modeling of the mesoscale constituents allows for a detailed examination of the deformation of mechanically loaded PBXs. In this study, a pair of explosive HMX crystals were bound with an HTPB based polymer binder and imaged using micro Computed Tomography ($µ$CT). Delamination of the polymer binder from the crystals was induced via macro-scale application of a tensile load. Cohesive parameters for simulation of the crystal/binder interface were determined by comparing numerical and experimental results of the delamination of this bi-crystal system. Simulations showed macro-scale plastic behavior is accommodated through both crystal-binder delamination and plastic deformation of the binder itself. Imaging and simulation of this bi-crystal system demonstrated the complex interplay of material constitutive behavior and grain-scale geometry. Successful parameterization of this relatively simple mechanical system is a step toward a better understanding void nucleation in similarly composed polycrystalline composites.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Implementation of a dislocation-density based single-crystal model into a continuum shock hydrodynamics code

The dynamic thermomechanical responses of polycrystalline materials under shock loading are often dominated by the interaction of defects and interfaces. For example, within metals, a prescribed deformation associated with a shock wave may be accommodated by crystallographic slip, provided a sufficient population of mobile dislocations is available. However, if the deformation rate is large enough, there may be an insufficient number of freely mobile dislocations. In this case, additional dislocations may be nucleated, or alternate mechanisms (e.g. twinning, damage) activated in order to accommodate the deformation. Furthermore, direct numerical simulation at the mesoscale offers insight into these physical processes that can be invaluable to the development of macroscale constitutive theories, if the mesoscale models adequately represent the anisotropic nonlinear thermomechanical response of individual crystals and their interfaces.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Ultrathin oxide layers for nanoscale integration of molecular light absorbers, catalysts, and complete artificial photosystems

Artificial photosynthesis is an attractive approach for the generation of renewable fuels because such systems will be suitable for deployment on highly abundant, non-arable land. Recently emerged methods of nanoscience to create conformal, ultrathin oxide layers enable the hierarchical integration of light absorbers, catalysts, and membranes into systems with far simpler synthetic approaches than available till now. This holds in particular for the coupling of molecular light absorbers and catalysts for sunlight to fuel conversion, providing photoelectrodes with greatly improved stability. Moreover, the use of ultrathin inert oxides as proton conducting, molecule impermeable membranes has opened up the integration of reduction and oxidation half reactions into complete photosynthetic systems on the shortest possible length scale-the nanometer scale. This capability affords minimization of energy-degrading resistance losses caused by ion transport over macroscale distances while separating the incompatible water oxidation and carbon dioxide reduction catalysis environments on the nanoscale. Understanding of charge transport between molecular components embedded in the oxide layers is critical for guiding synthetic design improvements of the light absorber-catalyst units to optimize performance and integrate them into complete artificial photosystems. Recent results and insights from transient optical, vibrational, and photoelectrochemical studies are presented, and future challenges and opportunities for engaging dynamic spectroscopies to accelerate the development of nanoscale integrated artificial photosystems are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Full stress tensor measurement using fluorescence spectroscopy

Photoluminescent spectral peak positions are known to shift as a function of mechanical stress state. This has been demonstrated at macroscales to determine mean stress and mesoscales to determine mean stress and a quantity related to shear stress. Here, we propose a method to utilize traction-free surface conditions and knowledge of material orientation to solve for two in-plane displacement fields given two measured spectral peak positions measured at a grid of points. It is then possible to calculate the full stress tensor at each measurement point. This is a significant advancement over the previous ability to measure one or two stress quantities. We validate the proposed method using a simple, two-grain geometry and show that it produces the same mean stress and shear stress measure as the existing direct method. Furthermore, we also demonstrate determination of the full stress field in a polycrystalline alumina specimen.

47 OTHER INSTRUMENTATION↗

Revealing nanoscale dynamics during an epoxy curing reaction with x-ray photon correlation spectroscopy

In this work the evolution of nanoscale properties is measured during the thermally triggered curing of an industrial epoxy adhesive. We use x-ray photon correlation spectroscopy (XPCS) to track the progression of the curing reaction through the local dynamics of filler particles that reflect the formation of a thermoset network. Out-of-equilibrium dynamics are resolved through identification and analysis of the intensity– intensity autocorrelation functions obtained from XPCS. The characteristic time scale and local velocity of the filler is calculated as functions of time and temperature. We find that the dynamics speed up when approaching the curing temperature (T cure ), and decay rapidly once T cure is reached. We compare the results from XPCS to conventional macroscale characterization by differential scanning calorimetry (DSC). The demonstration and implementation of nanoscale characterization of curing reactions by XPCS proves useful for future development and optimization of epoxy thermoset materials and other industrial adhesive systems.

36 MATERIALS SCIENCE↗

Nucleation and growth of voids in shock loaded copper bicrystals

Understanding the evolution of damage and deformation due to spall at grain boundaries can provide a basis for connecting micro- to macroscale failure behavior in metals under extreme conditions. Copper bicrystal samples were shock loaded using flyer-plate impacts in a light gas gun with shock stresses ranging from 3 to 6 GPa. Pulse duration as well as crystal orientation along the shock direction were varied for a fixed boundary misorientation to determine their effects on void nucleation and coalescence. Samples were soft recovered and cross-sectioned to characterize damage using electron backscattering diffraction and scanning electron microscopy to gather information on damage characteristics at and around the GB, with emphasis on growth of boundary and bulk voids. Chemistry and composition analysis were also performed on samples to determine if trace elements present in a sample affected the threshold for void nucleation. Results show that the kinetics of damage growth at the boundary are strongly affected by stress level and impurities. It was found that the boundary selected had a similar or even lower tendency to show damage than the bulk at low pulse durations and amplitudes. As pulse duration and amplitude increased damage localized at the boundary, which was found to consist of many small voids, indicating that the boundary experienced rapid void nucleation and coalescence. Furthermore, the presence of impurities correlated strongly with scatter on damage evolution.

36 MATERIALS SCIENCE↗

Machine learning approaches for structural and thermodynamic properties of a Lennard-Jones fluid

Predicting the functional properties of many molecular systems relies on understanding how atomistic interactions give rise to macroscale observables. However, current attempts to develop predictive models for the structural and thermodynamic properties of condensed-phase systems often rely on extensive parameter fitting to empirically selected functional forms whose effectiveness is limited to a narrow range of physical conditions. Here, we illustrate how these traditional fitting paradigms can be superseded using machine learning. Specifically, we use the results of molecular dynamics simulations to train machine learning protocols that are able to produce the radial distribution function, pressure, and internal energy of a Lennard-Jones fluid with increased accuracy in comparison to previous theoretical methods. The radial distribution function is determined using a variant of the segmented linear regression with the multivariate function decomposition approach developed by Craven et al. [J. Phys. Chem. Lett. 11, 4372 (2020)]. The pressure and internal energy are determined using expressions containing the learned radial distribution function and also a kernel ridge regression process that is trained directly on thermodynamic properties measured in simulation. The presented results suggest that the structural and thermodynamic properties of fluids may be determined more accurately through machine learning than through human-guided functional forms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Interface learning in fluid dynamics: Statistical inference of closures within micro–macro-coupling models

Many complex multiphysics systems in fluid dynamics involve using solvers with varied levels of approximations in different regions of the computational domain to resolve multiple spatiotemporal scales present in the flow. The accuracy of the solution is governed by how the information is exchanged between these solvers at the interface and several methods have been devised for such coupling problems. In this article, we construct a data-driven model by spatially coupling a microscale lattice Boltzmann method (LBM) solver and macroscale finite difference method (FDM) solver for reaction-diffusion systems. The coupling between the micro-macro solvers has one to many mapping at the interface leading to the interface closure problem, and we propose a statistical inference method based on neural networks to learn this closure relation. The performance of the proposed framework in a bifidelity setting partitioned between the FDM and LBM domain shows its promise for complex systems where analytical relations between micro-macro solvers are not available.

42 ENGINEERING↗

Probabilistic effects of porosity and chemical kinetics on the shock initiation of an octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) based explosive

In this work, we investigate the effects of porosity and chemical reactivity on the shock initiation of an HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine)-based polymer bonded explosive due to thin flyer impact using two-dimensional mesoscale simulations. Microstructural samples are computationally generated, and the calculated shock response is compared to experiments. The sample dimensions are 0.1×1.05 mm 2 , and each sample accounts for grain scale heterogeneities, including randomized grain morphology and pore spacing. Simulations are performed using the Lawrence Livermore National Laboratory (LLNL) arbitrary Lagrangian–Eulerian multi-physics code ALE3D. Chemistry and equation of states are handled with LLNL's thermochemical code Cheetah. We employ a modified Arrhenius-based chemical kinetics model that scales the temperature exponent fitting parameter with flyer shock conditions. This new model is calibrated to the mean response of experimental shock initiation thresholds for three separate Kapton flyer thicknesses. A simplified shock initiation criterion is proposed, which classifies the macroscale initiation response based on cross-sectional temperature profiles of the shocked sample. This criterion is used to further investigate the probabilistic response of the sample to variations in the mean porosity and chemical kinetics rate, and we find an inverse correlation between the shock pressure and the confidence interval of the experimental threshold velocity. Furthermore, the methodology presented in this work for quantifying the probabilistic nature of shock initiation is material agnostic. These mesoscale studies may be used to inform continuum scale models and expedite the development of future, tailored energetic formulations.

36 MATERIALS SCIENCE↗