Search NASA⌕ Search

SEARCH · Search NASA

Results for “Finite size effects”

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 271 records · Page 15

Resonance ultrasound prediction of residual stress within a hybrid layer for additively manufactured samples

Hybrid additive manufacturing (AM) involves secondary processes or energy sources to alter specified locations within the build volume. Each hybrid step can refine the grain size, increase dislocation density, or modify residual stresses. Typically, the changes in mechanical properties are not confined within a single layer but have a compounding effect on preceding layers. Existing methods of measuring AM residual stress are limited in terms of their sensitivity, or they are destructive measurements. We propose using resonant ultrasound spectroscopy (RUS) to measure the residual stress in hybrid-AM components noninvasively, based on changes to the resonances, compared to a stress-free component. In this paper, we use finite element models to simulate residual stress in hybrid-AM components and to examine the sensitivity of RUS measurements in terms of frequency shifts and mode shapes with respect to single hybrid layers. Then, the RUS results are used to predict stress for a layer at a known location with unknown stress. Here, the approach highlights the capabilities of RUS to address an AM characterization challenge.

36 MATERIALS SCIENCE↗

A high-frequency, low-power resonant radio-frequency neutron spin flipper for high-resolution spectroscopy

Here, we present a resonant-mode, transverse-field, radio-frequency (rf) neutron spin flipper design that uses high-temperature superconducting films to ensure sharp transitions between uniform magnetic field regions. Resonant mode allows for low-power, high-frequency operation but requires strict homogeneity of the magnetic fields inside the device. This design was found to efficiently flip neutrons at 96.6 ± 0.6% at an effective frequency of 4 MHz in bootstrap configuration with a beam size of 2.4 × 2.5 cm 2 and a wavelength of 0.4 nm. The high frequency and efficiency enable this device to perform high-resolution neutron spectroscopy with comparable performance with currently implemented rf flipper designs. The limitation of the maximum frequency was found due to the field homogeneity of the device. We numerically analyze the maximum possible efficiency of this design using a Bloch solver simulation with magnetic fields generated from finite-element simulations. We also discuss future improvements of the efficiency and frequency to the design based on the experimental and simulation results.

47 OTHER INSTRUMENTATION↗

Numerical-heating effects in atmospheric pressure streamer discharges simulated with a PIC code

Artificial heating in plasma simulations is a well-known phenomenon which occurs when, among other things, the Debye length is poorly resolved by the simulation mesh. Here, in this work, the degree to which numerical-heating occurs during a simulation of a nanosecond atmospheric pressure streamer discharge is examined. The streamer is simulated using a two-dimensional finite-element, particle-in-cell code Empire, which uses direct simulation Monte Carlo for binary particle interactions. Initially, an estimate of the numerical-heating rate applied to Empire is performed using a simple plasma model. Second, a positive atmospheric pressure streamer discharge simulation is performed to study the effects of numerical heating on plasma density, electron temperature, and streamer velocity. The nominal Debye length is approximately 1 μm and the amount of numerical heating introduced in the simulation is varied by using mesh sizes ranging from 2 μm to 20 μm. A measurable numerical heating quantity is proposed that can be used to estimate the appropriate element size and quantify the numerical-heating that can be expected over the simulation time for an atmospheric pressure streamer. In conclusion while Δx/λ D violations can be an issue it is not likely to be an issue with streamer discharges that are temporally short and occur in environments where collision frequencies are high. This result validates the rationale of grid size choices for a large amount of previously published works where Δx/λ D violation was not clearly addressed. Primary finding of this work is that numerical heating is of minor concern for plasma simulations where electron–neutral collisions are numerous such that multiple collisions can occur within a single plasma period.

Nikic, Dejan [University of New Mexico, Albuquerqu↗

Blockage effects in the chemotaxis of diffusiophoretic particles

Transport mechanisms at the micro- and nano-scale play an essential role in regulating intracellular organization. Recent work indicates that directed motion of constituents inside cells can emerge through diffusiophoretic transport, in which colloidal particles move under the influence of chemical gradients. Here, we examine how blockers—passive or actively consuming—reshape those gradients and thereby influence the motion of diffusiophoretic particles. By combining analytical solutions with finite element simulations, we first show that a single blocker can distort a background gradient enough to create or eliminate stagnation points, significantly modifying particle transport. We then introduce a second, explicitly sized blocker at one of these stagnation points and measure how its finite radius alters the diffusiophoretic velocity field for a test particle. Even moderate changes in the second blockers size can cause noticeable shifts in the substrate distribution, highlighting the importance of accounting for explicit particle radii under crowded or consumption-driven conditions. Our findings underscore that subtle geometric variations—such as the radii and positions of two or more blockers—can profoundly affect diffusiophoretic motion, providing a more complete picture of how blocking and crowding phenomena shape intracellular transport.

Song, Zehao [Northwestern Univ., Evanston, IL (Uni↗

Response of 11 B enriched ZrB 2 ultra-high temperature ceramic to neutron irradiation at elevated temperatures

ZrB 2 , an ultra-high temperature ceramic (UHTC) is being considered for use in fusion reactor first-wall structures, yet its response to irradiation remains poorly understood. This study employed scanning/transmission electron microscopy (S/TEM), synchrotron X-ray diffraction (XRD), finite element calculations, and thermal property measurements to thoroughly investigate the neutron-irradiation effects on 11 B-enriched ZrB 2 . Neutron irradiations were conducted at 220 °C and 620 °C, with a neutron fluence of 2.2 × 10 25 neutron/m 2 (energy > 0.1 MeV), resulting in 3.9 dpa and 4200 appm He. The study revealed the unusual prevalence of prism loops and a > c anisotropic lattice swelling, likely linked to the low c/a ratio of ZrB 2 , leading to grain boundary microcracking. Reducing the grain sizes was effective in reducing intergranular cracking and macroscopic swelling. The observation of cavities in ZrB 2 irradiated at 620 °C, as opposed to 220 °C, prompts questions about the temperature at which vacancies in ZrB 2 become mobile, and the role of neutron absorption by 10 B in elevating irradiation temperatures. Isotopic enrichment in 11 B proves to be a viable strategy for mitigating helium production in transition-metal diborides, which is a critical consideration for nuclear applications. Irradiation-induced defects reduce the thermal diffusivity and conductivity of ZrB 2 by a factor of 4–9, which has important implications for its role as a plasma-facing material in fusion reactors that drive high heat fluxes through first-wall materials. Here, this comprehensive study lays the foundation for understanding ZrB 2 behavior under neutron irradiation and highlights important phenomena to consider for various material applications.

36 MATERIALS SCIENCE↗

Contact pressure analysis to allow improved design for the clearing plate in a biomass comminution system

For biomass processing systems, feedstocks with small but equal sizes can improve flowability and thermochemical outputs. The Crumbler rotary shear system was designed to achieve such ideal feedstocks to reduce costs and energy throughout the biomass supply chain. Improving the wear resistance of major components in this system is important to decrease the cost. Accelerating the feedstock flow through the rotary shear machine plays an equivalent role of improving its efficiency. An analytical analysis and a finite-element analysis of the stationary clearing plate are presented in this work. Here, these tools are used to investigate the effect of contact surface curvature on the impact pressure. The optimized new clearing plate design is expected to improve the component’s lifetime and wood chips flowability. The field tests confirm that the trend of the new clearing plate design guided by FEA. The optimization method, model verification, and validation experience applied in this work can also be applied to other static components with similar simple contact wear problems.

09 BIOMASS FUELS↗

Fast Solution of Fully Implicit Runge--Kutta and Discontinuous Galerkin in Time for Numerical PDEs, Part I: the Linear Setting

Fully implicit Runge--Kutta (IRK) methods have many desirable properties as time integration schemes in terms of accuracy and stability, but high-order IRK methods are not commonly used in practice with numerical PDEs due to the difficulty of solving the stage equations. This paper introduces a theoretical and algorithmic preconditioning framework for solving the systems of equations that arise from IRK methods applied to linear numerical PDEs (without algebraic constraints). Additionally, this framework also naturally applies to discontinuous Galerkin discretizations in time. Under quite general assumptions on the spatial discretization that yield stable time integration, the preconditioned operator is proven to have condition number bounded by a small, order-one constant, independent of the spatial mesh and time-step size, and with only weak dependence on number of stages/polynomial order; for example, the preconditioned operator for 10th-order Gauss IRK has condition number less than two, independent of the spatial discretization and time step. The new method can be used with arbitrary existing preconditioners for backward Euler-type time-stepping schemes and is amenable to the use of three-term recursion Krylov methods when the underlying spatial discretization is symmetric. The new method is demonstrated to be effective on various high-order finite-difference and finite element discretizations of linear parabolic and hyperbolic problems, demonstrating fast, scalable solution of up to 10th-order accuracy. The new method consistently outperforms existing block preconditioning approaches, and in several cases, the new method can achieve 4th-order accuracy using Gauss integration with roughly half the number of preconditioner applications and wallclock time as required using standard diagonally IRK methods.

97 MATHEMATICS AND COMPUTING↗

Microstructure based modeling of friction stir welded joint of dissimilar metals using crystal plasticity

We report the friction stir welding (FSW) process shows promising results in joining dissimilar metals which are otherwise almost impossible to join using traditional welding techniques. Being a new technique, the deformation and the failure mechanism of the joints made by the FSW process needs to be investigated. In this work, a joint between AZ31 Mg alloy and DP590 steel is modeled using phenomenological crystal plasticity formulation on the mesoscale in the form of a representative volume element (RVE). The interface of the two materials is modeled using a cohesive zone model. A parametric study has been performed to understand the effect of grain size and interface fracture toughness as well as strength on the mechanical performance of the joint. It was found that the grain size of AZ31 Mg alloy, as well as DP590 steel, has little effect on the overall joint performance. On the other hand, interface fracture toughness and strength have a significant impact on the mechanical properties of the joint

36 MATERIALS SCIENCE↗

Measurement of Transport Properties of Woody Biomass Feedstock Particles Before and After Pyrolysis by Numerical Analysis of X-Ray Tomographic Reconstructions

Lignocellulosic biomass has a complex, species-specific microstructure that governs heat and mass transport during conversion processes. A quantitative understanding of the evolution of pore size and structure is critical to optimize conversion processes for biofuel and bio-based chemical production. Further, improving our understanding of the microstructure of biochar coproduct will accelerate development of its myriad applications. This work quantitatively compares the microstructural features and the anisotropic permeabilities of two woody feedstocks, red oak and Douglas fir, using X-ray computed tomography (XCT) before and after the feedstocks are subjected to pyrolysis. Quantitative analysis of the three-dimensional (3D) reconstructions allows for direct calculations of void fractions, pore size distributions and tortuosity factors. Next, 3D images are imported into an immersed boundary based finite volume solver to simulate gas flow through the porous structure and to directly calculate the principal permeabilities along longitudinal, radial, and tangential directions. The permeabilities of native biomass are seen to differ by three to four orders of magnitude in the different principal directions, but we find that this anisotropy is substantially reduced in the biochar formed during pyrolysis. The quantitative transport properties reported here enhance the ability of pyrolysis simulations to account for feedstock-specific effects and thereby provide a useful touchstone for the biorefining community.

09 BIOMASS FUELS↗

Effect of Part Size, Displacement Rate, and Aging on Compressive Properties of Elastomeric Parts of Different Unit Cell Topologies Formed by Vat Photopolymerization Additive Manufacturing

Due to its ability to achieve geometric complexity at high resolution and low length scales, additive manufacturing (AM) has increasingly been used for fabricating cellular structures (e.g., foams and lattices) for a variety of applications. Specifically, elastomeric cellular structures offer tunability of compliance as well as energy absorption and dissipation characteristics. However, there are limited data available on compression properties for printed elastomeric cellular structures of different designs and testing parameters. In this work, the authors evaluate how unit cell topology, part size, the rate of compression, and aging affect the compressive response of polyurethane-based simple cubic, body-centered, and gyroid structures formed by vat photopolymerization AM. Finite element simulations incorporating hyperelastic and viscoelastic models were used to describe the data, and the simulated results compared well with the experimental data. Of the designs tested, only the parts with the body-centered unit cell exhibited differences in stress–strain responses at different part sizes. Of the compression rates tested, the highest displacement rate (1000 mm/min) often caused stiffer compressive behavior, indicating deviation from the quasi-static assumption and approaching the intermediate rate response. The cellular structures did not change in compression properties across five weeks of aging time, which is desirable for cushioning applications. This work advances knowledge on the structure–property relationships of printed elastomeric cellular materials, which will enable more predictable compressive properties that can be traced to specific unit cell designs.

36 MATERIALS SCIENCE↗

Mechanical Comparison of Short Models of Nb$_3$ Sn Low-$\beta$ Quadrupole for the Hi-Lumi LHC

MQXF is the Nb 3 Sn Low-β quadrupole magnet that the HL-LHC project is planning to install in the LHC interaction regions in 2026 to increase the LHC integrated luminosity. The magnet will be fabricated in two different lengths: 4.2 m for MQXFA, built in the US by the Accelerator Upgrade Project (AUP), and 7.15 m for MQXFB, fabricated by CERN. In order to qualify the magnet design and characterize its performance with different conductors, cable geometries and pre-load configurations, five short model magnets, called MQXFS, were fabricated, assembled and tested. We compare the mechanical behavior of short model magnets using experimental data and new numerical models that take into account the measured coil sizes as a function of position.

43 PARTICLE ACCELERATORS↗

Effect of sample size on the maximum value distribution of fatigue driving forces in metals and alloys

An analytical framework is presented to predict the effects of sample size on the maximum value distribution (MVD) of the driving forces for fatigue crack formation in metals and alloys. The distribution of the maximum driving force for fatigue crack formation over the domain follows the generalized extreme value theory in the limit as the domain size increases to infinity. Here, a simulation-based analysis of microstructure influences on fatigue resistance for polycrystalline metals and alloys is very costly, and reaching those limits is intractable. This work models the MVD of Fatigue Indicator Parameters (FIPs), which serve as surrogate measures for the driving force for fatigue crack formation, at finite sample sizes prior to their convergence to a limiting extreme value distribution. Large-scale crystal plasticity finite element (CPFE) simulations of FCC Al 7075-T6 with microstructure realizations of various sizes are incorporated to calibrate and evaluate the developed framework, and a total of ∼6.5 million grains of Al 7075-T6 are examined. The calibrated analytical solution agrees well with the brute force Monte Carlo simulation results extracted from the CPFE simulations. Furthermore, the developed formulation can predict the MVD of FIPs for different sample sizes using a size-dependent parameter, and it is capable of accurately extrapolating the MVD of FIPs for much larger microstructure sample sizes than the size used for its calibration.

Crystal plasticity↗

Benchmarking Microscale Ductility Measurements (Final Report of the Project DE-NE0008799)

Conventional macroscale experimentation is generally considered to be straightforward with few limitations. Conversely, micro/nanoscale experimentation presents numerous challenges in loading device design, sample preparation and handling, as well as accurate understanding of grain size and local texture effects on recorded measurements. Despite these challenges, nanopillar compression, MEMs based micro-tension, and nanoindentation approaches have been able to provide fundamental contributions to the understanding of material behavior at small lengthscales. However, the overarching shortcoming of these micro/nanoscale experimentation approaches, is the inability to directly translate measurements evaluated at the nm and µm length scales (e.g., hardness) to macroscale tensile material behavior (i.e., elastic modulus, yield strength, and ductility). The objectives of the proposed study are, 1) to establish best practices for obtaining tensile microscale ductility measurements, and 2) to validate methodologies to for comparing microscale ductility measurements to macroscale ductility measurements. In order to achieve these objectives, a multi-lengthscale, multi-temperature testing protocol and simulation framework are executed first on copper as a model material to validate the following approach, and second on reactor grade Zircaloy-2. Experiments are conducted on specimens extracted from the same test piece to ensure nominally identical grain size and texture from specimens to specimen. Motivated by the need to isolate the contribution of size-effects on obtained mechanical property measurements, specimens are manufactured with thicknesses at the micro- (1-10 µm), meso- (10-100's µm), and macroscales (sub-sized ASTM E8). In-situ full-field deformation techniques (scanning electron microscopy (SEM) grid methods and optical DIC) are incorporated into testing at each specimen length-scale to capture plasticity localization and evolution. Experimental testing for all specimens is conducted at both room temperature and elevated temperatures to probe the role of thermal activation on plastic deformation accommodation processes. Simulation efforts focus on examining the mechanical behavior of microscale specimens using a finite element approach with explicitly resolved grain morphologies, and an embedded crystal plasticity model. The cost-efficient implementation method allows for the modeling of a statistically significant number of both real (i.e., digital twin) and generated microstructures to obtain an understanding of the interrelationships between specimen microstructure and geometric variables (grain size, texture, specimen geometry, etc.) on microscale mechanical behavior.

36 MATERIALS SCIENCE↗

Holographic measurement and bulk teleportation

Holography has taught us that spacetime is emergent and its properties depend on the entanglement structure of the dual theory. In this paper, we describe how changes in the entanglement due to a local projective measurement (LPM) on a subregion A of the boundary theory modify the bulk dual spacetime. We find that LPMs destroy portions of the bulk geometry, yielding post-measurement bulk spacetimes dual to the complementary unmeasured region A c that are cut off by end-of-the-world branes. Using a bulk calculation in AdS 3 and tensor network models of holography (in particular, the HaPPY code and random tensor networks), we show that the portions of the bulk geometry that are preserved after the measurement depend on the size of A and the state we project onto. The post-measurement bulk dual to A c includes regions that were originally part of the entanglement wedge of A prior to measurement. This suggests that LPMs performed on a boundary subregion A teleport part of the bulk information originally encoded in A into the complementary region A c . In semiclassical holography an arbitrary amount of bulk information can be teleported in this way, while in tensor network models the teleported information is upper-bounded by the amount of entanglement shared between A and A c due to finite-N effects. When A is the union of two disjoint subregions, the measurement triggers an entangled/disentangled phase transition between the remaining two unmeasured subregions, corresponding to a connected/disconnected phase transition in the bulk description. Our results shed new light on the effects of measurement on the entanglement structure of holographic theories and give insight on how bulk information can be manipulated from the boundary theory. They could also represent a first step towards a holographic description of measurement-induced phase transitions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Influence of chemical strains on the electrocaloric response, polarization morphology, tetragonality, and negative-capacitance effect of ferroelectric core-shell nanorods and nanowires

Using Landau-Ginzburg-Devonshire (LGD) approach, we proposed the analytical description of the influence of chemical strains on spontaneous polarization and the electrocaloric response in ferroelectric core-shell nanorods. We postulate that the nanorod core presents a defect-free single-crystalline ferroelectric material, and elastic defects are accumulated in the ultrathin shell, where they can induce tensile or compressive chemical strains. Finite-element modeling (FEM) based on the LGD approach reveals transitions of domain-structure morphology induced by chemical strains in the Ba Ti O 3 nanorods. Namely, tensile chemical strains induce and support the single-domain state in the central part of the nanorod, while the curled domain structures appear near the unscreened or partially screened ends of the rod. The vortexlike domains propagate toward the central part of the rod and fill it entirely, when the rod is covered by a shell with compressive chemical strains above some critical value. The critical value depends on the nanorod sizes, aspect ratio, and screening conditions at its ends. Both analytical theory and FEM predict that the tensile chemical strains in the shell increase the nanorod polarization, lattice tetragonality, and electrocaloric response well above the values corresponding to the bulk material. The physical reason for the increase is strong electrostriction coupling between the mismatch-type elastic strains induced in the core by chemical strains in the shell. Comparison with earlier XRD data confirmed an increase of the tetragonality ratio in tensile Ba Ti O 3 nanorods compared to the bulk material. Obtained analytical expressions, which are suitable for the description of strain-induced changes in a wide range of multiaxial ferroelectric core-shell nanorods and nanowires, can be useful for strain engineering of advanced ferroelectric nanomaterials for energy storage, harvesting, electrocaloric applications, and negative capacitance elements. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗

Deformation and Delamination of Dynamically Bulging Bilayered Films

We employ finite element simulations to study the dynamic bulging of films under shock tube conditions and compare these responses with quasi-static bulging. Both scenarios subject the film to time-varying pressure. Two film types are considered: uniform material films and bilayers with a relatively soft interior layer and a hard exterior layer. The interface between the materials in the bilayer cases is modeled as a cohesive zone that follows a bilinear normal traction-separation cohesive law. Additionally, bilayer films with preexisting cracks in the interior layer are analyzed. Parametric studies are also conducted to explore the effects of loading rate and specimen size. The mechanical response of all materials is assumed to be strain-rate independent, focusing solely on the inertial effects in the response of the films under dynamic loading, and the absence of such effects under quasi-static loading. The results indicate that quasi-static bulging exhibits a predictable plate bending like deformation. In contrast, depending on the loading rate and specimen size, dynamic bulging can give rise to elastic waves, and a sequence of deformation processes including initial uniform acceleration, followed by bending-like deformation, and finally through-thickness biaxial stretching. Dynamic loading also leads to larger bulge formation and greater deformation than quasi-static loading. A key finding is the role of preexisting cracks in driving delamination under dynamic loading, which is absent in quasi-static loading. This work enhances our understanding of dynamic bulging and motivates further research on characterizing delamination in layered films using dynamic bulge tests.

36 MATERIALS SCIENCE↗

The Third Dimension of Ferroelectric Domain Walls

Abstract Ferroelectric domain walls are quasi‐2D systems that show great promise for the development of nonvolatile memory, memristor technology, and electronic components with ultrasmall feature size. Electric fields, for example, can change the domain wall orientation relative to the spontaneous polarization and switch between resistive and conductive states, controlling the electrical current. Being embedded in a 3D material, however, the domain walls are not perfectly flat and can form networks, which leads to complex physical structures. In this work, the importance of the nanoscale structure for the emergent transport properties is demonstrated, studying electronic conduction in the 3D network of neutral and charged domain walls in ErMnO 3 . By combining tomographic microscopy techniques and finite element modeling, the contribution of domain walls within the bulk is clarified and the significance of curvature effects for the local conduction is shown down to the nanoscale. The findings provide insights into the propagation of electrical currents in domain wall networks, reveal additional degrees of freedom for their control, and provide quantitative guidelines for the design of domain‐wall‐based technology.

36 MATERIALS SCIENCE↗

Effect of time and spatial domains on monolayer 2D material interface thermal conductance measurement using ns ET-Raman

Raman-based approaches to characterizing thermal properties of nanoscale and microscale materials have proven to be powerful experimental methods. Here, this work explores the potential of energy-transport state-resolved Raman (ET-Raman) for the precise determination of interfacial thermal conductance (G) between monolayer 2D materials and their substrates. ET-Raman employs a continuous wave laser and nanosecond (ns) pulsed laser to induce both steady state and transient heating while simultaneously collecting 2D material temperature information from the Raman signal. A 3D numerical model built using the finite volume method simulates the two differentiated thermal transport processes. The simulated temperature rises are compared with experimental Raman temperature information to extract G. In this study, we investigate the temporal effects of the pulse frequency and pulse width of the ns laser on precise characterization of G. By running multiple simulations for a range of pulsed laser frequencies and pulse widths, we determine that ET-Raman is most sensitive to precise interfacial thermal conductance measurements with short pulse widths (~below 20 ns) and with pulsing frequencies that ensure sufficient cooling time between consecutive pulses (~below 500 kHz). We also conclude that substrate thickness and laser spot size play critical roles in heat distribution and subsequent sensitivity to meaningful G measurements. As both laser spot size and SiO 2 substrate thickness decrease, detection of interfacial thermal transport becomes more difficult. These results offer essential guidance in the design of future Raman-based approaches to thermal measurements using pulsed lasers.

2D materials↗