Search NASA⌕ Search

SEARCH · Search NASA

Results for “ISOTROPISM”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

On the numerical accuracy in finite-volume methods to accurately capture turbulence in compressible flows

The goal of the present article is to understand the impact of numerical schemes for the reconstruction of data at cell faces in finite-volume methods, and to assess their interaction with the quadrature rule used to compute the average over the cell volume. Here, third-, fifth- and seventh-order WENO-Z schemes are investigated. On a problem with a smooth solution, the theoretical order of convergence rate for each method is retrieved, and changing the order of the reconstruction at cell faces does not impact the results, whereas for a shock-driven problem all the methods collapse to first-order. Here, study of the decay of compressible homogeneous isotropic turbulence reveals that using a high-order quadrature rule to compute the average over a finite-volume cell does not improve the spectral accuracy and that all methods present a second-order convergence rate. However the choice of the numerical method to reconstruct data at cell faces is found to be critical to correctly capture turbulent spectra. In the context of simulations with finite-volume methods of practical flows encountered in engineering applications, it becomes apparent that an efficient strategy is to perform the average integration with a low-order quadrature rule on a fine mesh resolution, whereas high-order schemes should be used to reconstruct data at cell faces.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Novel pore size-controlled, susceptibility matched, 3D-printed MRI phantoms

We report the design concept and fabrication of MRI phantoms, containing blocks of aligned microcapillaires that can be stacked into larger arrays to construct diameter distribution phantoms or fractured, to create a “powder-averaged” emulsion of randomly oriented blocks for vetting or calibrating advanced MRI methods, that is, diffusion tensor imaging, AxCaliber MRI, MAP-MRI, and multiple pulsed field gradient or double diffusion-encoded microstructure imaging methods. The goal was to create a susceptibility-matched microscopically anisotropic but macroscopically isotropic phantom with a ground truth diameter that could be used to vet advanced diffusion methods for diameter determination in fibrous tissues. Two-photon polymerization, a novel three-dimensional printing method is used to fabricate blocks of capillaries. Double diffusion encoding methods were employed and analyzed to estimate the expected MRI diameter. Susceptibility-matched microcapillary blocks or modules that can be assembled into large-scale MRI phantoms have been fabricated and measured using advanced diffusion methods, resulting in microscopic anisotropy and random orientation. This phantom can vet and calibrate various advanced MRI methods and multiple pulsed field gradient or diffusion-encoded microstructure imaging methods. We demonstrated that two double diffusion encoding methods underestimated the ground truth diameter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Evolution of anisotropy with saturation and its implications for the elastoplastic responses of clay rocks

Many clay rocks have distinct bedding planes. Experimental studies have shown that their mechanical properties evolve with the degree of saturation (DOS), often with higher stiffness and strength after drying. For transversely isotropic rocks, the effects of saturation can differ between the bed-normal (BN) and bed-parallel (BP) directions, which gives rise to saturation-dependent stiffness and strength anisotropy. Accurate prediction of the mechanical behavior of clay rocks under partially saturated conditions requires numerical models that can capture the evolving elastic and plastic anisotropy with DOS. In this study, we present an anisotropy framework for coupled solid deformation-fluid flow in unsaturated elastoplastic media. Here we incorporate saturation-dependent strength anisotropy into an anisotropic modified Cam-Clay (MCC) model and consider the evolving anisotropy in both the elastic and plastic responses. The model was calibrated using experimental data from triaxial tests to demonstrate its capability in capturing strength anisotropy at various levels of saturation. Through numerical simulations, we demonstrate the role of evolving stiffness and strength anisotropy in the mechanical behavior of clay rocks. Plane strain simulations of triaxial compression tests were also conducted to demonstrate the impacts of material anisotropy and DOS on the mechanical and fluid flow responses.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Sequential formulation of all‐way coupled finite strain thermoporomechanics for largely deformable gas hydrate deposits

We develop a numerically stable sequential formulation of thermoporomechanics for largely deformable gas hydrate deposits, extended from the fixed stress split of infinitesimal transformation. Constitutive equations are based on the total Lagrangian approach for both flow and geomechanics, including dynamic full tensor permeability and thermal conductivity updated from the deformation gradient. For space discretization, we take the cell-centered finite volume and node-based finite element method for flow and geomechanics, respectively. Then, we propose a sequential implicit method for all-way coupled thermoporomechanics, where the nonisothermal multiphase flow problem of gas hydrates is solved implicitly first and then the geomechanics problem is solved implicitly at the next step. During solution of the flow problem, we fix the rate of first Pioal total stress for numerical stability as well as apply porosity correction and entropy correction to account for geomechanical effects. We test numerical examples where flow and geomechanics parameters are based on deep oceanic gas hydrate deposits. When applying depressurization, even though the results between the infinitesimal transformation and finite strain geomechanics are similar in the early stages due to small deformation, we find differences between them in the late times as deformation becomes large. Accordingly, permeability and thermal conductivity tensors become nonisotropic full tensors although they are initially isotropic. Furthermore, we identify numerical stability of the developed sequential method from the test cases that exhibit the highly complex coupled gas hydrate systems with large deformation. Thus, the proposed sequential formulation can be applied in largely deformable gas hydrate systems.

42 ENGINEERING↗
DOE OSTI2025

A Poromechanical Framework for Internal Interactions Induced by Solid Inclusions

The framework of poromechanics is generalized to simulate the multiscale behavior of porous media subjected to internal loadings stemming from the growth of solid inclusions. This generalization is designed to enable the study of anisotropic internal stress generation from solid growth within the pores, while recovering isotropic fluid-induced loading as a particular case. For this purpose, a mathematical strategy to define constitutive tensors in a thermodynamically consistent form is proposed, thus offering new opportunities for determining the poromechanical properties of a porous solid through advanced experimentation or micromechanical models. The framework is specialized by means of established elastic solutions for single pore–matrix interaction, as well as through homogenization schemes considering the interaction among congruent pores. In particular, the second Eshelby solution and the Tanaka–Mori–Benveniste homogenization scheme are used to derive a microporoelastic model. At an elemental scale, the model is tested under mixed control conditions by replicating different scenarios of geomaterial testing. In addition, the model characteristics are outlined with reference to inelastic microscopic loadings replicating chemo-mechanical forcing, such as expansive crystal formation. Through a series of parametric analyses, it is shown that the microstructure of the pores significantly influences the properties of porous media. Most notably, it is shown that the effects of a solid forming within the pores depend in a highly nonlinear fashion on the constitutive characteristics of the inhomogeneities and can therefore not be readily quantified or predicted without models capturing the diverse multiscale interactions among pores, inhomogeneities, and matrix.

chemomechanics↗
DOE OSTI2021

Multiscale design of nonlinear materials using a Eulerian shape optimization scheme

Motivated by recent advances in manufacturing, the design of materials is the focal point of interest in the material research community. One of the critical challenges in this field is finding optimal material microstructure for a desired macroscopic response. This work presents a computational method for the mesoscale-level design of particulate composites for an optimal macroscale-level response. The method relies on a custom shape optimization scheme to find the extrema of a nonlinear cost function subject to a set of constraints. Three key “modules” constitute the method: multiscale modeling, sensitivity analysis, and optimization. Multiscale modeling relies on a classical homogenization method and a nonlinear NURBS-based generalized finite element scheme to efficiently and accurately compute the structural response of particulate composites using a nonconformal discretization. A three-parameter isotropic damage law is used to model microstructure-level failure. An analytical sensitivity method is developed to compute the derivatives of the cost/constraint functions with respect to the design variables that control the microstructure's geometry. The derivation uncovers subtle but essential new terms contributing to the sensitivity of finite element shape functions and their spatial derivatives. Several structural problems are solved to demonstrate the applicability, performance, and accuracy of the method for the design of particulate composites with a desired macroscopic nonlinear stress-strain response.

42 ENGINEERING↗
DOE OSTI2021

Plastic work constrained elastoplastic topology optimization

An elastoplastic topology optimization framework for limiting plastic work generation while maximizing stiffness is presented. The kinematics and constitutive model are based on finite strain linear isotropic hardening plasticity, and the balance laws are solved using a total Lagrangian finite element formulation. Aggregation of the specific plastic work combined with an adaptive normalization scheme efficiently constrains the maximum specific plastic work. The optimization problem is regularized using an augmented partial differential equation filter, and is solved by the method of moving asymptotes where path-dependent sensitivities are derived using the adjoint method. The numerical examples show a clear dependence on the optimized maximum stiffness structures for different levels of constrained specific plastic work. It is also shown that due to the history dependency of the plasticity, the load path significantly influences the structural performance and optimized topology.

42 ENGINEERING↗
DOE OSTI2024

A multiscale anisotropic polymer network model coupled with phase field fracture

The study of polymers has continued to gain substantial attention due to their expanding range of applications, spanning essential engineering fields to emerging domains like stretchable electronics, soft robotics, and implantable sensors. These materials exhibit remarkable properties, primarily stemming from their intricate polymer chain network, which, in turn, increases the complexity of precisely modeling their behavior. Especially for modeling elastomers and their fracture behavior, accurately accounting for the deformations of the polymer chains is vital for predicting the rupture in highly stretched chains. Despite the importance, many robust multiscale continuum frameworks for modeling elastomer fracture tend to simplify network deformations by assuming uniform behavior among chains in all directions. Recognizing this limitation, our study proposes a multiscale fracture model that accounts for the anisotropic nature of elastomer network responses. At the microscale, damage in the chains is assumed to be driven by both the chain's entropy and the internal energy due to molecular bond distortions. In order to bridge the stretching in the chains to the macroscale deformation, we employ the maximal advance path constraint network model, inherently accommodating anisotropic network responses. As a result, chains oriented differently can be predicted to exhibit varying stretch and, consequently, different damage levels. To drive macroscale fracture based on damages in these chains, we utilize the micromorphic regularization theory, which involves the introduction of dual local-global damage variables at the macroscale. The macroscale local damage variable is obtained through the homogenization of the chain damage values, resulting in the prediction of an isotropic material response. The macroscale global damage variable is subjected to nonlocal effects and boundary conditions in a thermodynamically consistent phase field continuum formulation. Moreover, the total dissipation in the system is considered to be mainly due to the breaking of the molecular bonds at the microscale. To validate our model, we employ the double-edge notched tensile test as a benchmark, comparing simulation predictions with existing experimental data. Additionally, to enhance our understanding of the fracturing process, we conduct uniaxial tensile experiments on a square film made up of polydimethylsiloxane (PDMS) rubber embedded with a hole and notches and then compare our simulation predictions with the experimental observations. Furthermore, we visualize the evolution of stretch and damage values in chains oriented along different directions to assess the predictive capacity of the model. In conclusion, the results are also compared with another existing model to evaluate the utility of our model in accurately simulating the fracture behavior of rubber-like materials.

42 ENGINEERING↗
DOE OSTI2022

Tailoring compatibilization potential of maleic anhydride–grafted polypropylene by sequential rheochemical processing of polypropylene and polyamide 66 blends

Reactive compatibilization of immiscible polymer blends using a preferred compatibilizer leaves the following question: How much loading of a compatibilizer is good enough and what maximum properties can be achieved? A good understanding of the process can help solve the reutilization of mixed waste plastics. Here, the reactive compatibilization of polypropylene-graft-maleic anhydride (PP-g-MAH) on polypropylene/polyamide 66 (PP/PA66) blends is quantitatively assessed using thermorheological, microoptical, spectroscopic, and x-ray scattering-based characterization tools. The overall compositions of all PP(60%)/PA66(30%)/PP-g-MAH(10%) blends are kept constant while systematically controlling the degree of the chemical reaction by varying the sequential addition of PP-g-MAH to the melt mixture. The first feeding of PP-g-MAH (x%, 0 ≤ x ≤ 10) is conducted at 280°C and the second feeding of (10–x)% is at 200°C. During the first step, a larger amount of chemical activity is observed up to 4%–6% addition of PP-g-MAH. Furthermore, the constant chemical composition allows a systematic comparison of the compatibilizer-dependent thermal, rheological, morphological, and mechanical properties. With an increased degree of interfacial reaction, the size of dispersed PA66 domains decreases to register improved interfacial adhesion with the PP matrix, yielding enhanced Young's modulus and absolute failure strength of the isotropic matrix.

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

On the cogent formulation of an elastomeric silicone ink material for direct ink write ( DIW ) 3D printing

Abstract Adhesives and sealants show fine rheology with good physical and mechanical properties as viscous pastes, a possible starting point for developing direct ink writing (DIW) 3D printing ink. However, many commercial adhesives and sealants take days or weeks to cure fully. DIW 3D‐printed parts made directly from these sealants are not designed for a scalable manufacturing process and high‐volume production. Moreover, most of these adhesives and sealants have volume shrinkage during cure. A systematic understanding of formulation methods and design principles for an elastomeric silicone DIW ink can overcome these issues. This study presents the cogent formulation development of a 3D printable thermoset elastomer silicone that gels and cures isotropically in minutes, reducing cycle time for rapid ink development with no shrinkage during cure. More specifically, we outline the principles of raw material selection of a formulation to achieve excellent rheology, printability, synchronized working, and gel time fitting requirements closer to scalable manufacturing. The reaction kinetics and their corresponding 3D‐printed structural properties are also described. Interest in future work is toward a rational DIW 3D printing ink material development protocol and use of machine learning (ML). Highlights Formulation method flexibility and design principle of DIW ink. Raw material selection principle to achieve optimal rheology for DIW printing. Ink gel kinetics for large‐scale DIW manufacturing. Hydrosilylation conversion over time at different ambient temperatures. Structural properties of DIW 3D printed parts.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Imaging a solvent-swollen polymer gel network by open liquid transmission electron microscopy

Visualizing the network of a solvent-swollen polymer gel remains problematic. To address this challenge, open transmission electron microscopy (TEM) was applied to thin gel films permeated by a nonvolatile ionic liquid. The targeted physical gels were prepared by cooling concentrated solutions of poly(ethylene glycol) in 1-ethyl-3-methyl imidazolium ethyl sulfate [EMIM][EtSO 4 ]. During the cooling, gelation occurred by a frustrated crystallization of the dissolved polymer, leading to a percolated, solvent-permeated semicrystalline network in which nanoscale polymer crystals acted as crosslinks. Crystalline features ranging from ~5 to ~200 nm were observed, with the visible network strands dominantly consisting of long curvilinear crystallites of ~15–20 nm diameter. Nascent spherulites irregularly decorated the network, creating a complex structural hierarchy that complicated analyses. Lacking diffraction contrast, TEM did not visualize the many disordered, fully solvated PEG chains present in the voids between crystals. Finally, recognizing that a network's three dimensionality is ambiguous when assessed through two-dimensional microscopy projections, a small gel region was studied by TEM tomography, revealing a nearly isotropic three-dimensional arrangement of the curvilinear crystallites, which displayed remarkably uniform cylindrical cross sections.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Bayesian D‐Optimal Designs for Gaussian Process Surrogate Models

Computer experiments often employ space-filling strategies to create surrogate models with strong predictive performance. The impact of model parameter estimation for Gaussian process surrogates, however, is often overlooked. Obtaining a better initial estimate of the covariance lengthscale parameter, θ, can greatly improve the resulting Gaussian process fit through more effective sequential acquisitions during active learning. In this work, we propose a novel initial design maximizing the Bayesian D-optimality criterion of the Gaussian process lengthscale parameter. Previously published results have shown the emphasis on lengthscale estimation to be promising, but relied on an empirically driven design creation process. Our Bayesian D-optimal designs are rooted in information theory and lead to more informative sequential acquisitions by improving lengthscale estimation. In many cases, these gains eventually result in better surrogates than those seeded with space-filling initial designs. Furthermore, Bayesian D-optimal designs can be tailored to either isotropic or anisotropic covariance structures, and the Bayesian framework enables the inclusion of prior knowledge in the design process, offering greater flexibility and adaptability. Through several simulation studies, we demonstrate the advantages of Bayesian D-optimal designs in terms of both lengthscale estimation accuracy and predictive performance during active learning.

Bayesian experimental design↗
DOE OSTI2023

Band gap dependence of exciton Lande factor in lead halide perovskites

The optical properties of lead halide perovskite semiconductors in vicinity of the bandgap are controlled by excitons, so that investigation of their fundamental properties is of critical importance. The exciton Landé or g-factor g X is the key parameter, determining the exciton Zeeman spin splitting in magnetic fields. The exciton, electron, and hole carrier g-factors provide information on the band structure, including its anisotropy, and the parameters contributing to the electron and hole effective masses. Here, g X is measured by reflectivity in magnetic fields up to 60 T for lead halide perovskite crystals. The materials band gap energies at a liquid helium temperature vary widely across the visible spectral range from 1.520 up to 3.213 eV in hybrid organic–inorganic and fully inorganic perovskites with different cations and halogens: FA 0.9 Cs 0.1 PbI 2.8 Br 0.2 , MAPbI 3 , FAPbBr 3 , CsPbBr 3 , and MAPb(Br 0.05 Cl 0.95 ) 3 . The exciton g-factors are found to be nearly constant, ranging from +2.3 to +2.7. Thus, the strong dependences of the electron and hole g-factors on the bandgap roughly compensate each other when combining to the exciton g-factor. The same is true for the anisotropies of the carrier g-factors, resulting in a nearly isotropic exciton g-factor. The experimental data are compared favorably with model calculation results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Evidence of Directional Structural Superlubricity and Lévy Flights in a van der Waals Heterostructure

Abstract Structural superlubricity is a special frictionless contact in which two crystals are in incommensurate arrangement such that relative in‐plane translation is associated with vanishing energy barrier crossing. So far, it has been realized in multilayer graphene and other van der Waals (2D crystals with hexagonal or triangular crystalline symmetries, leading to isotropic frictionless contacts. Directional structural superlubricity, to date unrealized in 2D systems, is possible when the reciprocal lattices of the two crystals coincide in one direction only. Here, directional structural superlubricity a α‐bismuthene/graphite van der Waals system is evidenced, manifested by spontaneous hopping of the islands over hundreds of nanometers at room temperature, resolved by low‐energy electron microscopy and supported by registry simulations. Statistical analysis of individual and collective α‐bismuthene islands populations reveal a heavy‐tailed distribution of the hopping lengths and sticking times indicative of Lévy flight dynamics, largely unobserved in condensed‐matter systems.

Le Ster, Maxime↗
DOE OSTI2020

An Improved Fracture Mechanics-Informed Multiscale Thermomechanical Damage Model for Ceramic Matrix Composites

This paper extends recent work done by the authors in modeling length scale-dependent damage behavior of ceramic matrix composites (CMCs) to include effects of local anisotropy introduced by matrix cracking. This model captures scale-dependent damage initiation and propagation behavior of the brittle matrix by employing internal state variable (ISV) theory within a multiscale modeling framework to obtain damaged matrix stress/strain constitutive relationships at each length scale. The damage ISV captures the effects of matrix cracking and growth by using fracture mechanics and the self-consistent scheme to determine the reduced stiffness of the cracked matrix. Matrix cracks, which activate when stress intensity factors near manufacturing induced cavities exceed the fracture toughness of the material, are assumed to be transversely isotropic in the plane of the crack, and matrix anisotropy occurs when the damaged stiffness tensor is rotated from the crack plane to the global axes. The crack progression and temporal evolution of the damage ISV are governed by fracture mechanics and crack growth kinetics. The model effectively captures first matrix cracking, which is the first significant deviation from linear elasticity. The nonlinear predictive capabilities of the material model are demonstrated for monolithic silicon carbide (SiC) and a 2D woven five-harness satin (5HS) carbon fiber SiC matrix (C/SiC) CMC.

Skinner, Travis↗
DOE OSTI2022

Holographic Kolmogorov-Sinai entropy and the quantum Lyapunov spectrum

In classical chaotic systems the entropy, averaged over initial phase space distributions, follows a universal behavior. While approaching thermal equilibrium it passes through a stage where it grows linearly, while the growth rate, the Kolmogorov-Sinai entropy (rate), is given by the sum over all positive Lyapunov exponents. A natural question is whether a similar relation is valid for quantum systems. We argue that the Maldacena-Shenker-Stanford bound on quantum Lyapunov exponents implies that the upper bound on the growth rate of the entropy, averaged over states in Hilbert space that evolve towards a thermal state with temperature T, should be given by πT times the thermal state’s von Neumann entropy. Strongly coupled, large N theories with black hole duals should saturate the bound. To test this we study a large number of isotropization processes of random, spatially homogeneous, far from equilibrium initial states in large N, $\mathcal{N}$ = 4 Super Yang Mills theory at strong coupling and compute the ensemble averaged growth rate of the dual black hole’s apparent horizon area. We find both an analogous behavior as in classical chaotic systems and numerical evidence that the conjectured bound on averaged entropy growth is saturated granted that the Lyapunov exponents are degenerate and given by λ i = ±2πT. This fits to the behavior of classical systems with plus/minus symmetric Lyapunov spectra, a symmetry which implies the validity of Liouville’s theorem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Cosmic Axion Background detection using resonant cavity arrays

The axion is a well-motivated and generic extension of the Standard Model. If produced in the early universe, axions may still be relativistic today, forming a Cosmic Axion Background (CaB) potentially detectable in direct detection experiments. Although CaB is expected to be broadband, which makes it challenging to be detected, a high-quality-factor microwave cavity acts as a narrowband filter with response peaked at its resonant frequency. We propose a new strategy using multi-cavity arrays to distinguish signal from background noise by exploiting spatial correlations of the axion-induced electric field which are set by the cavity quality factor. We compute the two-point correlation function for electric fields in spatially separated cavities sourced by an isotropic CaB. Analyzing various cavity geometries, we find that stacked, wide-base cavity arrays offer coherent enhancement of the axion signal. We apply our formalism to prospective upgrades of the ADMX experiment, including configurations with four and eighteen coupled cavities. Although these arrays do not achieve a coherent enhancement, optimizing the geometry could potentially yield an $\mathcal{O}(1)$ improvement in the sensitivity to the CaB.

Axions and ALPs↗
DOE OSTI2024

The neutrino force in neutrino backgrounds: Spin dependence and parity-violating effects

The neutrino force results from the exchange of a pair of neutrinos. A neutrino background can significantly influence this force. In this work, we present a comprehensive calculation of the neutrino force in various neutrino backgrounds with spin dependence taken into account. In particular, we calculate the spin-independent and spin-dependent parity-conserving neutrino forces, in addition to the spin-dependent parity-violating neutrino forces with and without the presence of a neutrino background for both isotropic and anisotropic backgrounds. Compared with the vacuum case, the neutrino background can effectively violate Lorentz invariance and lead to additional parity-violating terms that are not suppressed by the velocity of external particles. We estimate the magnitude of the effect of atomic parity-violation experiments, and it turns out to be well below the current experimental sensitivity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗