Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Anisotropic”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Overview of a Versatile Loading System for Anisotropic Material Property Characterization

Additive manufacturing, cold rolling, and other thermomechanical treatments on metals, especially HCP, can cause texture and oriented grain structures resulting in anisotropic mechanical properties. This can lead to macroscopic response that significantly deviates from isotropic assumptions which motivates studying mechanical properties in multiple directions. In addition, multiaxial loading and complex stress states overlapping the above mentioned material anisotropy can lead to unexpected outcomes. Idaho National Laboratory developed an advanced mechanical testing system to study a range of uniaxial to multiaxial stress states while measuring the anisotropic response bringing insight into the overlap between material properties and stress states. High-temperature capability and stress- or strain-controlled loading is available to enable a variety of experiment types and conditions to measure elastic, plastic, and viscoplastic properties. This poster presents the design and capabilities of this system with preliminary results highlighting the benefits of multiaxial loading and anisotropic analysis in an integrated system.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

An Eulerian crystal plasticity framework for modeling large anisotropic deformations in energetic materials under shocks

Here, this paper demonstrates a novel Eulerian computational framework for modeling anisotropic elastoplastic deformations of organic crystalline energetic materials (EM) under shocks. While Eulerian formulations are advantageous for handling large deformations, constitutive laws in such formulations have been limited to isotropic elastoplastic models, which may not fully capture the shock response of crystalline EM. The present Eulerian framework for high-strain rates, large deformation material dynamics of EM incorporates anisotropic isochoric elasticity via a hypo-elastic constitutive law and visco-plastic single-crystal models. The calculations are validated against atomistic calculations and experimental data and benchmarked against Lagrangian (finite element) crystal plasticity computations for shock-propagation in a monoclinic organic crystal, octahydro-1,3,5,7-tetranitro-1,3,5,7 tetrazocine (β-HMX). The Cauchy stress components and the resolved shear stresses calculated using the present Eulerian approach are shown to be in good agreement with the Lagrangian computations for different crystal orientations. The Eulerian framework is then used for computations of shock-induced inert void collapse in β-HMX to study the effects of crystal orientations on hotspot formation under different loading intensities. The computations show that the hotspot temperature distributions and the collapse profiles are sensitive to the crystal orientations at lower impact velocities (viz., 500 m/s); when the impact velocity is increased to 1000 m/s, the collapse is predominantly hydrodynamic and the role of anisotropy is modest. The present methodology will be useful to simulate energy localization in shocked porous energetic material microstructures and other situations where large deformations of single and polycrystals govern the thermomechanical response.

42 ENGINEERING↗
DOE OSTI2022

Thermal stability of anisotropic bonded magnets prepared by additive manufacturing

Abstract In this research, anisotropic NdFeB + SmFeN hybrid and NdFeB bonded magnets are additively printed in a polyphenylene sulfide (PPS) polymer binder. Printed NdFeB + SmFeN PPS bonded magnets displayed excellent magnetic properties ( B r [remanence] = 6.9 kG [0.69 T], H cj [coercivity] = 8.3 kOe [660 kA/m], and BH max [energy product] = 9.9 MGOe [79 kJ/m 3 ]) with superior corrosion resistance and thermal stability. The anisotropic NdFeB bonded magnet shows a high coercivity of 14.6 kOe (1162 kA/m) with a BH max of 8.7 MGOe (69 kJ/m 3 ). The coercivity and remanence temperature coefficients for NdFeB + SmFeN hybrid bonded magnets are −0.10%/K and −0.46%/K, and for NdFeB bonded magnets are −0.14%/K and −0.53%/K in the range of 300–400 K, indicating that the hybrid bonded magnets are thermally stable. The average flux aging loss for hybrid magnets was also determined to be very stable over 2000 h at 448 K (175°C) in air with 2.04% compared to that of NdFeB magnets with 3.62%.

thermal stability↗
DOE OSTI2022

Hybrid large-eddy simulation with adaptive downstream anisotropic eddy viscosity model

Turbulence pervades most flows of engineering interest, and its prediction remains a challenge on both accuracy and cost. One promising predictive approach that reduces cost combines large-eddy simulation (LES) with simulation based on Reynolds averaged Navier–Stokes equations (RANS). This study presents a method to overcome stability and accuracy issues associated with these hybrid LES-RANS methods. The method developed involves extracting the Reynolds stresses from the upstream LES solution and then using that information to convert the downstream RANS model from a scalar eddy-viscosity model to an anisotropic nonlinear eddy-viscosity model. The method developed differs from the downstream tensorial eddy-viscosity model by being independent of the coordinate system. The method developed was evaluated by computing film cooling of a flat plate with the coolant injected through one row of circular holes. Results obtained show instabilities at the LES-to-RANS interface to be eliminated. Also, the method developed yielded solutions that compare reasonably well with those from LES, even though a significant portion of the flow is computed by the adapted anisotropic RANS model instead of LES, which significantly reduced the number of grid points and computational time needed. Since modification of the downstream RANS model is based on information extracted from the upstream LES solution, the method developed is adaptive to the problem being studied.

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 OSTI2020

Evaluation of measuring thermal conductivity of isotropic and anisotropic thermally insulating materials by transient plane source (Hot Disk) technique

The transient plane source (TPS) technique, also referred as the Hot Disk method, has been widely used due to its ability to measure the thermal properties of an extensive range of materials (solids, liquids, and powder). Recently, it has been recognized that typical Hot Disk sensors can influence TPS results of thermally insulating materials and lead to an overestimation of thermal conductivity. Although improvements have been proposed, they have not yet been implemented in the commercial TPS, leaving researchers with non-standardized modifications or options provided by a commercial Hot Disk apparatus. An empirical study of thermally insulating materials such as extruded polystyrene (XPS) and aerogel blanket is conducted in order to address the factors that affect the reliability of thermal conductivity k obtained using the commercial TPS apparatus. Sensor size, input power, duration of the measurements, applied pressure, and, in the case of anisotropic materials, heat capacity are investigated, and the results are compared with those using a Heat Flow Meter apparatus. The effect of sensor size on the k value is ascribed to heat loss through connecting leads and is more pronounced in smaller sensors and in materials with lower k values . In the case of XPS and aerogel, the effect becomes minimal for sensors with a radius r ≥ 6.4 mm. The low input power yields a high scattering of the results and should be avoided. Applied contact pressure and the tested region of the specimen play an important role in experiments with low-density fibrous materials due to the large percentage of heat being transferred by radiation and the heterogeneous nature of the samples, respectively. Additionally, the sensitivity of anisotropic measurements to the value of the material’s volumetric heat capacity ( ρC p ) is shown, emphasizing the need for the precise determination.

thermal conductivity↗
DOE OSTI2020

A novel framework for elastoplastic behaviour of anisotropic solids

A general framework for developing nonlinear hyperelastic/plastic constitutive laws for anisotropic solids experiencing large strains and strain rates has been developed. The proposed framework does not rely on the “a priori” known strain energy function, but instead introduces a physical decomposition of the material element into seven physically independent stress bearing mechanisms, each of which has a constitutive law in terms of internal moments described by a scalar function of a single variable. Furthermore, the model has been encoded into a combined finite-discrete element method and tested against static geomechanical test data. The numerical validation experiments show the model can reproduce plastic anisotropic behavior in both biaxial and uniaxial loading of a geomaterial.

58 GEOSCIENCES↗
DOE OSTI2021

Single molecule characterization of anomalous transport in a thin, anisotropic film

Here, the diffusion of small, charged molecules incorporated in an anisotropic polyelectrolyte multilayers (PEMs) was tracked in three dimensions by combining single-molecule fluorescence localization (to characterize lateral diffusion) with Förster resonance energy transfer (FRET) between diffusing molecules and the supporting surface (to measure diffusion in the surface-normal direction). Analysis of the surface-normal diffusion required model-based statistical analysis to account for the inherently noisy FRET signal. Combining these distinct single-molecule methods, which are inherently sensitive to different length scales, permitted simultaneous characterization of severely anisotropic diffusion, which was more than three orders of magnitude slower in the surface-normal direction. We hypothesize that the anomalously slow surface-normal diffusion was related to the periodic distribution of charge in the PEM, which created electrostatic barriers. The motion was strongly subdiffusive, with anomalous temporal scaling exponents in lateral and normal directions, suggesting a connection to the transient, random fractal conformation of polymer chains in the film’s matrix.

36 MATERIALS SCIENCE↗
DOE OSTI2020

FFT-based solver for higher-order and multi-phase-field fracture models applied to strongly anisotropic brittle materials

This paper presents the application of a fast Fourier transform (FFT) based method to solve two phase field models designed to simulate crack growth of strongly anisotropic materials in the brittle regime. By leveraging the ability of the FFT-based solver to generate solutions with higher-order and global continuities, we design two simple algorithms to capture the complex fracture patterns (e.g. sawtooth, and curved crack growth) common in materials with strongly anisotropic surface energy via the multi-phase-field and high-order phase-field frameworks. A staggered operator-split solver is used where both the balance of linear momentum and the phase field governing equations are formulated in the periodic domain. The unit phase field of the initial failure region is prescribed by the penalty method to alleviate the sharp material contrast between the initial failure region and the base material. The discrete frequency vectors are generalized to estimate the second and fourth order gradients such that the Gibbs effect near shape interfaces or jump conditions can be suppressed. Furthermore, a preconditioner is adopted to improve the convergence rate of the iterative linear solver. Three numerical experiments are used to systematically compare the performance of the FFT-based method in the multi-phase-field and high-order phase-field frameworks.

42 ENGINEERING↗
DOE OSTI2020

SO(3)-invariance of informed-graph-based deep neural network for anisotropic elastoplastic materials

This work examines the frame-invariance (and the lack thereof) exhibited in simulated anisotropic elasto-plastic responses generated from supervised machine learning of classical multi-layer and informed-graph-based neural networks, and proposes different remedies to fix this drawback. The inherent hierarchical relations among physical quantities and state variables in an elasto-plasticity model are first represented as informed, directed graphs, where three variations of the graph are tested. While feed-forward neural networks are used to train path-independent constitutive relations (e.g., elasticity), recurrent neural networks are used to replicate responses that depends on the deformation history, i.e. or path dependent. In dealing with the objectivity deficiency, we use the spectral form to represent tensors and, subsequently, three metrics, the Euclidean distance between the Euler Angles, the distance from the identity matrix, and geodesic on the unit sphere in Lie algebra, can be employed to constitute objective functions for the supervised machine learning. In this, the aim is to minimize the measured distance between the true and the predicted 3D rotation entities. Following this, we conduct numerical experiments on how these metrics, which are theoretically equivalent, may lead to differences in the efficiency of the supervised machine learning as well as the accuracy and robustness of the resultant models. Neural network models trained with tensors represented in component form for a given Cartesian coordinate system are used as a benchmark. Our numerical tests show that, even given the same amount of information and data, the quality of the anisotropic elasto-plasticity model is highly sensitive to the way tensors are represented and measured. The results reveal that using a loss function based on geodesic on the unit sphere in Lie algebra together with an informed, directed graph yield significantly more accurate rotation prediction than the other tested approaches.

42 ENGINEERING↗
DOE OSTI2022

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↗
DOE OSTI2025

A provably stable numerical method for the anisotropic diffusion equation in confined magnetic fields

We present a novel numerical method for solving the anisotropic diffusion equation in magnetic fields confined to a periodic box which is accurate and provably stable. We derive energy estimates of the solution of the continuous initial boundary value problem. A discrete formulation is presented using operator splitting in time with the summation by parts finite difference approximation of spatial derivatives for the perpendicular diffusion operator. Weak penalty procedures are derived for implementing both boundary conditions and parallel diffusion operator obtained by field line tracing. We prove that the fully-discrete approximation is unconditionally stable. Discrete energy estimates are shown to match the continuous energy estimate given the correct choice of penalty parameters. A nonlinear penalty parameter is shown to provide an effective method for tuning the parallel diffusion penalty and significantly minimises rounding errors. Several numerical experiments, using manufactured solutions, the “NIMROD benchmark” problem and a single island problem, are presented to verify numerical accuracy, convergence, and asymptotic preserving properties of the method. Finally, we present a magnetic field with chaotic regions and islands and show the contours of the anisotropic diffusion equation reproduce key features in the field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Anisotropic thermal-conductivity degradation in the phase-field method accounting for crack directionality

Dynamic loads applied to metals may lead to brittle or ductile fracture depending on the imposed strain rates, material properties, and specimen geometry. With an increase in the applied velocity, brittle to ductile failure mode transition may also be observed. Furthermore, at high strain rates, ductile fracture may be preceded by shear bands, which are narrow bands of intense plastic deformation, typically accompanied by a significant rise in temperature. Reliable models are needed to predict the response of metals subject to dynamic loads. In particular, capturing the interplay between heat conduction and crack propagation using the phase-field fracture method is still an open research field. To accurately capture the heat transfer physics across crack surfaces, damage models degrading thermal-conductivity are necessary. While isotropic thermal-conductivity degradation models were proposed, they may lead to errors as they are indifferent to crack directionality. In our current work, a new anisotropic approach is proposed in which thermal-conductivity, which depends on the phase-field gradient, is degraded solely across the crack. Here, it is shown that this approach improves the near-field approximation of temperature and heat flux compared with isotropic degradation, when taking the discontinuous crack solutions as reference. Additionally, the anisotropic degradation approach is implemented in a unified dynamic fracture model and a couple of examples demonstrate the viability of this approach.

42 ENGINEERING↗
DOE OSTI2026

Finite-element boundary-integral simulation of thin wires and inhomogeneous penetrable bodies in subsurface multilayered anisotropic media

With the prevailing presence of drilling wells near the subsurface in mature oil and gas fields, the application of electromagnetic methods can be particularly challenging where the electromagnetic field is affected by the steel casing. In the past decades, borehole-to-surface and crosswell electromagnetic methods have been utilized for monitoring of reservoir and underground CO 2 storage. This paper presents a unified finite-element boundary-integral (FEBI) method capable of simultaneously modeling the complex electromagnetic interactions between thin metallic wires (representing steel casings) with 3D trajectory and arbitrary 3D inhomogeneous penetrable bodies (such as CO 2 plumes or hydrocarbon reservoirs) within anisotropic multilayered subsurface environments. Unlike existing approaches that treat these components separately or require dense discretization, or are limited to vertical wells, our unified formulation preserves flexible electromagnetic coupling while delivering improved computational efficiency. Assuming the background formation is multilayered anisotropic media, the surface integral equation method is applied to model the thin wires and boundaries of the inhomogeneous bodies. Meanwhile, the finite element method is applied to model the volume of inhomogeneous bodies. Here, the performance of the proposed FEBI method is assessed through comparison with reference numerical results and its practical significance is demonstrated through CO 2 plume monitoring scenarios.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Kink band and shear band localization in anisotropic perfectly plastic solids

Shear-driven strain localization has been observed in a wide variety of materials and may take the form of shear bands or kink bands. Based on observations of kink bands in plastically anisotropic metallic nanolaminates and single crystal metals, here we posit that, for the specific case of isochoric deformation, the kinematics of kink band formation are indistinguishable from those of plane strain shear band formation. The only distinction between shear bands and kink bands in these systems would then be that kink bands ‘lock up’ at a particular value of material rotation while shear bands may progress to arbitrarily high strains. In order to investigate whether strong material anisotropy is sufficient to arrest shear localization at a geometry that matches the classic kink band geometry, we model the development of a band of simple shear within an anisotropic perfectly plastic material. The resulting analytical model provides the stress state needed to maintain the kinematics of simple shear as a function of material anisotropy, deformation band orientation, and shear strain (or equivalently, material rotation). It is found that plastic anisotropy can promote either kink band or shear band formation depending on the loading orientation. When the deviatoric stress is positive parallel to the plane of anisotropy, shear localization may progress without bound and a shear band is produced. When the deviatoric stress is negative parallel to the plane of anisotropy, shear localization is arrested after a certain material rotation, resulting in a kink band. Examination of the requisite applied stress state during kink band formation provides an explanation for the experimentally-observed ‘lock up’ geometry. Solutions for the band boundary inclination angle are obtained and used to provide bounds on permissible band angles for both shear bands and kink bands.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A review of the application of 2D isotropic-anisotropic correlation NMR spectroscopy in structural studies of chalcogenide glasses

The application of solid-state high-resolution NMR spectroscopy in the structural investigation of chalcogenide glasses in Ge/As/P/Si-X (X = S, Se,Te) systems has remained challenging even for the spin-1/2 nuclides ( 29 Si, 31 P, 77 Se, 125 Te), owing to their low natural abundance (except for 31 P), slow spin-lattice relaxation rate and large CSA and chemical shift distribution induced line broadening effects. However, most of these deleterious effects can be successfully overcome in two-dimensional (2D) isotropic-anisotropic correlation NMR experiments, especially when performed at high magnetic field and in conjunction with the Car-Purcell-Meiboom-Gill (CPMG) echo train acquisition. Here, we present a short introduction to the basic principles of such experiments and review their applications over the last decade in deciphering various short- and intermediate- range structural characteristics of chalcogenide glasses in S-Se, Se-Te, Ge-Se, As-Se and Si-Se systems as well as in investigating the molecular dynamics in a P-Se supercooled liquid. We anticipate possible future applications of these 2D isotropic-anisotropic correlation NMR experiments, particularly in conjunction with density functional theory-based calculations of NMR chemical shift tensor parameters and additional signal enhancement schemes, in addressing complex structural correlations and distributions in chalcogenide glasses.

2D NMR↗
DOE OSTI2021

Anisotropic Transient Disordering of Colloidal, Two-Dimensional CdSe Nanoplatelets upon Optical Excitation

Nanoplatelets (NPLs)-colloidally synthesized, spatially anisotropic, two-dimensional semiconductor quantum wells-are of intense interest owing to exceptionally narrow transition line widths, coupled with solution processability and bandgap tunability. However, given large surface areas and undercoordinated bonding at facet corners and edges, excitation under sufficient intensities may induce anisotropic structural instabilities that impact desired properties. We employ time-resolved X-ray diffraction to study the crystal structure of CdSe NPLs in response to optical excitation. Photoexcitation induces greater out-of-plane than in-plane disordering in 4 and 5 monolayer (ML) NPLs, while 3 ML NPLs display the opposite behavior. Recovery dynamics suggest that out-of-plane cooling slightly outpaces in-plane cooling in 5 ML NPLs with recrystallization occurring on indistinguishable time scales. In comparison, for zero-dimensional CdSe nanocrystals, disordering is isotropic and recovery is faster. These results favor the use of NPLs in optoelectronic applications, where they are likely to exhibit superior performance over traditional, zero-dimensional nanocrystals.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2026

Anisotropic Exciton Transport in a Lamellar CsPbBr 3 Nanocrystal Superlattice

Colloidal self-assembly is one strategy for engineering anisotropic properties into otherwise isotropic materials. In this work, we demonstrate anisotropic exciton transport in an A 2 B-type superlattice containing columns of 5.3 nm CsPbBr 3 nanocubes assembled into a hexagonal lattice around 6.5 nm LaF 3 nanodisks. Using transient photoluminescence microscopy, we determined that the exciton diffusivity along the fast axis of the superlattice is more than twice as large as that along the slow axis at T = 5 K, but that anisotropy is greatly suppressed at room temperature. Calculations of the diffusivity anisotropy ratio based on Förster theory overestimate the measured values, highlighting the limitations of this theory in completely describing exciton transport. Overall, our results demonstrate how self-assembly of colloidal nanocrystals can be used to engineer directional energy transport, and raise more questions about the microscopic nature of dipole coupling in CsPbBr 3 NC superlattices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗