Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite Elements”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
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 OSTI2022

Reciprocal swimming at intermediate Reynolds number

In Stokes flow, Purcell's scallop theorem forbids objects with time-reversible (reciprocal) swimming strokes from moving. In the presence of inertia, this restriction is eased and reciprocally deforming bodies can swim. A number of recent works have investigated dimer models that swim reciprocally at intermediate Reynolds numbers Re≈1–1000. These show interesting results (e.g. switches of the swim direction as a function of inertia) but the results vary and seem to be case specific. Here, in this work, we introduce a general model and investigate the behaviour of an asymmetric spherical dimer of oscillating length for small-amplitude motion at intermediate Re. In our analysis we make the important distinction between particle and fluid inertia, both of which need to be considered separately. We asymptotically expand the Navier–Stokes equations in the small-amplitude limit to obtain a system of linear partial differential equations. Using a combination of numerical (finite element) and analytical (reciprocal theorem, method of reflections) methods we solve the system to obtain the dimer's swim speed and show that there are two mechanisms that give rise to motion: boundary conditions (an effective slip velocity) and Reynolds stresses. Each mechanism is driven by two classes of sphere–sphere interactions, between one sphere's motion and (1) the oscillating background flow induced by the other's motion, and (2) a geometric asymmetry induced by the other's presence. We can thus unify and explain behaviours observed in other works. Our results show how sensitive, counterintuitive and rich motility is in the parameter space of finite inertia of particles and fluid.

42 ENGINEERING↗
DOE OSTI2021

An FFT-based approach for Bloch wave analysis: application to polycrystals

A method based on the Fast Fourier Transform is proposed to obtain the dispersion relation of acoustic waves in heterogeneous periodic media with arbitrary microstructures. The microstructure is explicitly considered using a voxelized Representative Volume Element (RVE). The dispersion diagram is obtained solving an eigenvalue problem for Bloch waves in Fourier space. To this aim, two linear operators representing stiffness and mass are defined through the use of differential operators in Fourier space. The smallest eigenvalues are obtained using the implicitly restarted Lanczos and the subspace iteration methods, and the required inverse of the stiffness operator is done using the conjugate gradient with a preconditioner. The method is used to study the propagation of acoustic waves in elastic polycrystals, showing the strong effect of crystal anistropy and polycrystaline texture on the propagation. It is shown that the method combines the simplicity of classical Fourier series analysis with the versatility of Finite Elements to account for complex geometries proving an efficient and general approach which allows the use of large RVEs in 3D.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A heat flux sensor leveraging the transverse Seebeck effect in elemental antimony

Certain configurations of anisotropic single crystal materials can generate a thermoelectric voltage orthogonal to an induced temperature gradient. This phenomenon is known as the Transverse Seebeck Effect (TSE) and can be leveraged to fabricate simple and robust heat flux sensors. Only a small number of materials have been considered as TSE-based transducers and, among these, few have been developed into sensors with ruggedization against chemical and mechanical degradation. Here, we report on the fabrication and characterization of a rugged TSE-based heat flux sensor using prismatic antimony single crystals. The heat flux sensor was tested under static and dynamic heating scenarios. The sensor has a linear responsivity of 16.8 µV/(W/cm2) to heat fluxes spanning more than two orders of magnitude and a time constant of 4.4 s. The sensor’s response to localized heating, probed with a laser scanning technique, validated that the transduction mechanism is primarily the TSE by ruling out a sizable contribution from the conventional Seebeck effect. Finite element analysis corroborated that components used in the sensor package are the primary determinants of the time constant and the decrement of the responsivity from its theoretical maximum. Design principles that may be applied to elicit a faster transient response or higher responsivity are proposed. Furthermore, the results establish single crystal antimony as a promising transducer material for heat flux measurement systems and demonstrate potential effects of ruggedization on sensor performance.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Discrete Green’s functions and spectral graph theory for computationally efficient thermal modeling

Here, this work concerns solutions of the heat equation with the spectral graph method, for which the temperature is defined at discrete points in the domain and the spatial relationship among the points is described by a graph. The heat equation on the graph is solved using matrix techniques involving the eigenvectors and eigenvalues of the Laplacian matrix. The spectral graph approach precludes the computationally intensive meshing and numerous time-integration steps of the finite element method. In the present work, the spectral graph method is extended to include heat loss at the boundaries with a generalized boundary condition, and physics-based edge weights are introduced which simplify the calibration process. From this approach a discrete Green’s function is defined which allows for solutions under a variety of heating conditions including: space-varying initial conditions; time-and-space varying internal heating; and, time-and-space-varying heating at boundaries of type 1 (Dirichlet), type 2 (Neumann) and type 3 (Robin). Results are provided for benchmark heat transfer problems in one spatial dimension and in three spatial dimensions, and verification is provided by comparison with exact analytical solutions and finite difference solutions. The spectral graph method converges within 0.4% error of the analytical solution. The practical utility of the approach is demonstrated by thermal simulation of a multilayer additive manufacturing process. The spectral graph results are compared to experimentally-obtained temperature data for two metal parts, with error less than 5% of the experimental measurements, with computation time less than one minute on a desktop computer.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Smoothed boundary method for simulating incompressible flow in complex geometries

Simulating flow through porous media with explicit considerations of complex microstructures is very challenging using conventional sharp-interface methods because of the difficulties in generating meshes conformal to complex geometries. In this work, a diffuse interface embedded boundary method known as the Smoothed Boundary Method (SBM) is utilized to facilitate simulations of fluid dynamics involving complex geometries. In diffuse-interface methods, the geometry is described by a domain parameter. The SBM allows the straightforward reformulation of the time-dependent Navier-Stokes equations in terms of this domain parameter, using only algebraic identities. Thus, enforcing the appropriate boundary conditions at the irregular embedded boundary is greatly simplified. Adaptive mesh refinement is used to increase the accuracy of the diffuse interface method by allowing thinner interfacial thickness to be used in the domain parameter. Here, the SBM-formulated Navier-Stokes equations are solved with the Finite Difference Method on refined mesh systems. Sharp-interface Finite Element Method simulations using the commercial software COMSOL on body-conforming meshes are also provided for comparison. Favorable agreement between the two methods is observed. Since it is no longer necessary for the mesh to conform to the complex geometry, the grid system for the SBM simulations can be generated rapidly and without additional manual interventions, making the entire simulation process more expedient.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Characterization of W production during ICRF operations: experiments and modeling

For successfully heating plasma with waves in the ion cyclotron range of frequencies (ICRFs), mitigating impurity production is just as crucial as maximizing power coupling, especially in high-Z environments. ICRF can effectively deposit energy on ions, modify turbulence-driven transport, and enhance fusion reaction efficiency, but only when its power coupling has minimal impact on impurity production. To do so, one must rely on a toroidal array of at least three active elements excited with appropriate phasing and power ratio to reduce the currents induced on the antenna frame below levels critical for physical sputtering. In contrast to classic two-strap antennas, which are optimized for dipole phasing with equal power on both straps, three-strap antennas in ASDEX Upgrade (AUG)—but also four-strap antennas in JET, Alcator C-Mod, SPARC and ITER—offer the possibility to act also on the power ratio between the central and outer straps. With optimal settings, impurity production can be reduced substantially, making the ICRF compatible with the high-Z wall. This paper explores the characteristics of the AUG three-strap antennas in terms of impurity production, as well as the key role of plasma composition in this process. Numerical simulations were performed using SSWICH and Petra-M (finite element codes) to quantify impurity production and compare with experimental results. Energies of ions falling on antenna limiters (measured with probes) are well predicted by both codes. These tools are then used to further describe the source of the impurity, namely the gross erosion of tungsten from an ICRF antenna, for different plasma mixtures. Results are also compared to spectroscopy data. Ultimately, we show that deleterious effects of the ICRF on plasma surface interactions will be weaker in plasmas containing larger fractions of highly ionized heavier low-Z impurity, which is typically relevant for experiments relying on impurity seeding.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Challenges and Lessons Learned from fabrication, testing and analysis of eight MQXFA Low Beta Quadrupole magnets for HL-LHC

By the end of October 2022, the US HL-LHC Accelerator Upgrade Project (AUP) had completed fabrication of ten MQXFA magnets and tested eight of them. The MQXFA magnets are the low beta quadrupole magnets to be used in the Q1 and Q3 Inner Triplet elements of the High Luminosity LHC. This AUP effort is shared by BNL, Fermilab, and LBNL, with strand verification tests at NHMFL. An important step of the AUP QA plan is the testing of MQXFA magnets in a vertical cryostat at BNL. The acceptance criteria that could be tested at BNL were all met by the first four production magnets (MQXFA03-MQXFA06). Subsequently, two magnets (MQXFA07 and MQXFA08) did not meet some criteria and were disassembled. Lessons learned during the disassembly of MQXFA07 caused a revision to the assembly specifications that were used for MQXFA10 and subsequent magnets. In this paper, we present a summary of: 1) the fabrication and test data of all the MQXFA magnets; 2) the analysis of MQXFA07/A08 test results with characterization of the limiting mechanism; 3) the outcome of the investigation, including the lessons learned during MQXFA07 disassembly; and 4) the finite element analysis correlating observations with test performance.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Challenges and Lessons Learned From Fabrication, Testing, and Analysis of Eight MQXFA Low Beta Quadrupole Magnets for HL-LHC

By the end of October 2022, the US HL-LHC Accelerator Upgrade Project (AUP) had completed fabrication of ten MQXFA magnets and tested eight of them. The MQXFA magnets are the low-beta quadrupole magnets to be used in the Q1 and Q3 Inner Triplet elements of the High Luminosity LHC. This AUP effort is shared by BNL, Fermilab, and LBNL, with strand verification tests at NHMFL. An important step of the AUP QA plan is the testing of MQXFA magnets in a vertical cryostat at BNL. The acceptance criteria that could be tested at BNL were all met by the first four production magnets (MQXFA03-MQXFA06). Subsequently, two magnets (MQXFA07 and MQXFA08) did not meet some of the criteria and were disassembled. Furthermore, lessons learned during the disassembly of MQXFA07 caused a revision to the assembly specifications that were used for MQXFA10 and subsequent magnets. In this article, we present a summary of 1) the fabrication and test data for all the MQXFA magnets; 2) the analysis of MQXFA07/A08 test results with characterization of the limiting mechanism; 3) the outcome of the investigation, including the lessons learned during MQXFA07 disassembly; and 4) the finite element analysis correlating observations with test performance.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Mixing Cell Model: A One-Dimensional Numerical Model for Assessment of Water Flow and Contaminant Transport in the Unsaturated Zone

This report describes the Mixing Cell Model code, a one dimensional model for water flow and solute transport in the unsaturated zone under steady state or transient flow conditions. The model is based on the principles and assumptions underlying mixing-cell model formulations. The unsaturated zone is discretized into a series of independent mixing cells. Each cell may have unique hydrologic, lithologic, and sorptive properties. Ordinary differential equations describe the material (water and solute) balance within each cell. Water-flow equations are derived from the continuity equation, assuming that unit gradient conditions exist at all times in each cell. Pressure gradients are considered implicitly through model discretization. Unsaturated hydraulic conductivity and moisture contents are determined by the material specific moisture-characteristic curves. Solute-transport processes include explicit treatment of advective processes, first order chain decay, and linear sorption reactions. Dispersion is addressed through implicit and explicit dispersion. Implicit dispersion is an inherent feature of all mixing-cell models and originates from the formulation of the problem in terms of mass balance around fully mixed volume elements. Expressions are provided that relate implicit dispersion to the physical dispersion of the system. Two FORTRAN codes were developed to solve the water flow and solute-transport equations: (1) the Mixing Cell Model for Flow (MCMF) solves transient water-flow problems and (2) the Mixing Cell Model for Transport (MCMT) solves the solute-transport problem. The transient water-flow problem is typically solved first by estimating the water flux through each cell in the model domain as a function of time using the MCMF code. These data are stored in either ASCII or binary files that are later read by the solute transport code (MCMT). Code output includes solute pore water concentrations, water and solute inventories in each cell and at each specified output time, and water and solute fluxes through each cell and specified output time. Computer run times for coupled transient water flow and solute transport were typically several seconds on a 2 GHz Intel Pentium IV desktop computer. The model was benchmarked against analytical solutions and finite element approximations to the partial differential equations (PDE) describing unsaturated flow and transport. Differences between the maximum solute flux estimated by the mixing cell model and the PDE models were typically less than two percent. This revision includes an option for a fixed concentration lower boundary condition for diffusive fluxes for versions 020321 and later.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Enhancing the cooling performance of thermocouples: a power-constrained topology optimization procedure

Abstract Heat pumping through thermoelectric devices has many advantages over traditional cooling. However, their current efficiency is a limiting factor in their implementation. In this paper, we approach the non-convex topology optimization of thermoelectrical elements for cooling applications through the method of moving asymptotes (MMA) to improve their cooling capabilities per watt usage. The optimization problem is defined for a given power budget, aiming for the minimum temperature with a known heat pumping need. The introduction of power as a constraint justifies the introduction of the voltage gradient across the thermocouple as a design variable to maintain the thermoelectrical device in its optimum power-to-heat extraction ratio. To better understand the convergence of this non-convex problem, we present a two-variable analytical thermoelectric optimization model. This example provides information on how to select the penalty parameters used to scale the three material coefficients involved in the problem to obtain lower objective values and better convergence using MMA. The analytical model shows the non-convexity of the problem and provides the recommendation to use penalization coefficients of the form $$p_k=p_{\sigma }>p_{\alpha }=1$$ p k = p σ > p α = 1 for the thermal conductivity, electrical conductivity, and Seebeck coefficients. We tested these penalization coefficients through optimizations of a model based on the 1MC10-031 commercial thermoelectric-cooler (TEC) using the finite element method (FEM). These penalization coefficients provided local minima without the need for volume constraints. With this procedure, we found designs that provided temperatures close to 10 degrees lower using 60% less semiconductor material volume compared to the initial design.

Gutiérrez, G. Reales↗
DOE OSTI2025

Tailoring composition and deformation modes at the microstructural level for next generation low-cost high-strength austenitic stainless steels

The objective of this project is to enable deliberate development of cost-effective, hydrogen resistant alloys by establishing detailed relationships specific to the effects of alloy composition, short-range order (SRO), and microsegregation in the presence of hydrogen on the transition between homogeneous deformation and localized plasticity in shear bands. In collaboration with the International Institute for Carbon-Neutral Energy Research, I2CNER, at Kyushu University in Japan, we conceptualized, designed, and manufactured four austenitic alloys that maintain corrosion resistance and ensure lower cost relative to baseline commercial alloys. The mechanical properties and deformation modes of the novel alloys (KU alloys) were assessed in the presence of hydrogen (H). Correlations between composition and performance revealed that two of the KU alloys are suitable replacements for 316 steel, while another is a viable replacement for 304 steel at room temperature. We found that, in the presence of other austenite stabilizing elements namely Mn and N, replacing Ni with Cu does not lead to martensite formation as has been previously reported.1–3 Furthermore, we found that the addition of Cu leads to an earlier onset of multiple slip resulting in an relative earlier onset of a higher work hardening rate (WHR). Greater understanding of the relationships between alloy composition and SRO required the development of a novel advanced electron diffraction methodology to characterize SRO in complex FCC alloys. This innovative approach, which combines fluctuation and correlation analyses of diffuse-scattering signals, successfully differentiated between SRO and long-range ordering (LRO). Further investigations into annealed austenitic stainless steels could provide insights into manipulating SRO and its effects on material properties. Atomistic simulations provided understanding of SRO behavior that was difficult to capture experimentally. This project created the first spin cluster expansion model that is able to capture and describe SRO effects in Fe-Ni-Cr FCC alloys, accounting for the non-negligible effects of magnetism. An automated computational workflow was established to provide reliable predictions of SRO in Fe-Ni-Cr austenitic alloys, both with and without the presence of H atoms. Analysis of the propensity for SRO in Fe-Ni-Cr alloys revealed that H tends to cluster with specific, well-defined SRO domains. The computational framework is general purpose and can be extended to realistic stainless steels across diverse composition ranges. With confidence that SRO is possible in austenitic stainless steels, we developed a discrete dislocation finite element code to understand the interaction of dislocations with SRO in the presence of H. By incorporating H effects on the dislocation emission and SRO stress field we show that the critical stress for the dislocation pileup to breakthrough the SRO domain decreases in the presence of H, which directly contributes localized deformation at the macroscale. Through the simulation of a uniaxial tension test, we demonstrated that H-induced weakening of SRO stress field and H-enhanced dislocation emission can lead to the onset of shear localization at lower macroscopic strains. As a whole, this project identified three novel alloys that show improvements in performance and cost efficiency for H-facing applications by studying correlations between alloy chemistry and deformation behavior. We also made significant advancements to experimental and computational methodologies necessary to study the chemistry and distribution of SRO across a range of alloys, which in turn allowed us to demonstrate how deformation mechanisms change due to the contributions of SRO in austenitic alloys in the presence of H. The combined advancements in fundamental understanding with novel alloy development in this project has increased the viability of next generation H-technologies for the broader public through accessible low-cost alloys and accelerated development towards future H-infrastructure.

08 HYDROGEN↗
DOE OSTI2026

Mechanical Behavior of Neutron Irradiated Refractory Multi-Principal Element Alloys Processed via Spark Plasma Sintering

The search for advanced materials capable of withstanding the extreme conditions of Generation IV reactors is a critical area in materials science research. These reactors operate under severe environments, including high temperatures, corrosion, stress, and irradiation damage. Consequently, there is a need for innovative alloy systems to ensure the reliability and longevity of proposed Generation IV reactor components. Refractory multi-principal-element alloys (RMPEA) have emerged as a promising candidate due to their exceptional properties. These alloys, characterized by their composition of multiple principal elements in near-equiatomic ratios, exhibit superior resistance to irradiation damage, reduced void swelling, enhanced microstructural stability, and minimal irradiation-induced hardening. While initial studies on RMPEAs have shown promising results, most research has been limited to thin films, nanocrystalline microstructures, and ion irradiation, which do not accurately represent the behavior of bulk materials. To address this gap, our research focused on the neutron irradiation of bulk RMPEAs. We aim to conduct comprehensive post-irradiation examinations (PIE) of RMPEAs irradiated at the Advanced Test Reactor at Idaho National Laboratory. The RMPEAs were synthesized using spark plasma sintering (SPS) with mechanically alloyed metallurgical powder. The RMPEA specimens are a MoNbTi alloy system with additions of -Zr, and -ZrV. Furthermore, PIE consisted of mechanical testing and advanced materials characterization. The mechanical testing consisted of sub-sized tensile testing, micro- and nano- indentation. Microstructural characterization included scanning electron microscopy and transmission electron microscopy. Mechanical testing coupled with advanced microscopy techniques provides insight into phase morphology and its effects on the mechanical properties of the RMPEA specimens. The results indicate that both pristine and irradiated RMPEA specimens exhibited brittle behavior during tensile testing, which can be attributed to their heterogeneous microstructure. The SPS manufacturing process did not include any post treatment, which resulted in a heterogeneous microstructure. Energy-dispersive X-ray spectroscopy revealed the presence of intermetallic such as laves phases within the microstructure. Specifically, Ti-rich precipitates were observed in the MoNbTi specimen, while Mo-rich precipitates were found in the MoNbTiZrV specimen. Nano-hardness testing of pristine samples showed that the laves phases exhibited higher hardness values compared to the matrix phase, suggesting that precipitate hardening is likely the dominant hardening mechanism in these specimens. The results from this work will be used to build a finite element model to predict mechanical behavior of future MPEA compositions. Thus, enabling for a streamlined approach to developing novel MPEAs for the nuclear industry.

36 - MATERIALS SCIENCE↗
DOE OSTI2022

Radio-Frequency Rapid Thermal Processing Enabling Spatial Phase Transformation and Nanocrystallization of Soft Magnetic Amorphous Alloys

Thermal processing of soft magnetic amorphous and nanocrystalline alloys is explored under the influence of radio-frequency induction-heating techniques. Direct induction-heating concepts based on longitudinal and transverse flux heating are examined and the details of electromagnetic fields interaction with metallic strips are discussed by analytical calculations as well as finite element analysis. Initial experimental results confirming spatial control of phase transformations and nanocrystallization within a single strip of Finemet Fe-based amorphous ribbons are reported. Here, the degree to which primary and secondary crystallization temperature are achieved depends on the spacing between the ribbon relative to the induction coil as well as the coil design and configuration. For transverse coil configurations, the local temperature and therefore microstructural evolution is different across the lateral dimension of processed ribbons, with reduced gap sizes producing enhanced peak temperatures and larger temperature distributions with greater spatial variation in microstructure. In addition, indirect susceptor-based induction heating under tension is performed and the impact of microstructure is demonstrated. Herein, potential for exploiting spatially optimized phase transformations is illustrated through electromagnetic field–assisted processing in a scalable manufacturing process with amorphous and nanocrystalline soft magnetic alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Selective Heating Through Y-Junction Waveguide Designed by Acoustic Shape Optimization

Unlike phononic crystals or systems designed by topology optimization, waveguides designed by shape optimization do not have voids or internal defects, making the fabrication process more suitable for additive manufacturing. By designing a Y-junction waveguide through shape optimization, an ultrasonic wave can be controlled so that it propagates to a predetermined location just by adjusting its frequency. These demultiplexed ultrasonic waves can be used to transport signals or stimulate nearby materials. As an example, the ultrasonic wave is converted to heat at different locations, which can be applied to mechanisms that can take advantage of heating. First, shape optimization is performed on a cylindrical structure to selectively propagate ultrasonic waves of a particular frequency while attenuating others, which is analyzed through a finite element model. The numerical study results are compared with experimental measurements from samples fabricated through additive manufacturing methods. After verifying the concept, the Y-junction waveguide is fabricated to demultiplex the wave and selectively heat different locations. Here the results show that the method of combining shape optimization with additive manufacturing is exceptionally simple and capable of demultiplexing ultrasonic waves, which can replace complex electrical components with single-material waveguides.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Kirigami-Inspired Self-Assembly of 3D Structures

Self-assembly of 3D structures introduce an attractive and scalable route to realize reconfigurable and functionally capable mesoscale devices without human intervention. A common approach for achieving this is to utilize stimuli-responsive folding of hinged structures, which requires the integration of different materials and/or geometric arrangements along the hinges. It is also demonstrated that the inclusion of Kirigami cuts in planar, hingeless bilayer thin sheets can be used to produce complex 3D shapes in an on-demand manner. Nonlinear finite element models are developed to elucidate the mechanics of shape morphing in bilayer thin sheets and verify the predictions through swelling experiments of planar, millimeter-scaled PDMS (polydimethylsiloxane) bilayers in organic solvents. Building upon the mechanistic understandings, The transformation of Kirigami-cut simple bilayers into 3D shapes such as letters from the Roman alphabet (to make “ADVANCED FUNCTIONAL MATERIALS”) and open/closed polyhedral architectures is experimentally demonstrated. A possible application of the bilayers as tether-less optical metamaterials with dynamically tunable light transmission and reflection behaviors is also shown. As the proposed mechanistic design principles could be applied to a variety of materials, this research broadly contributes toward the development of smart, tetherless, and reconfigurable multifunctional systems.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Boosting Piezoelectricity by 3D Printing PVDF‐MoS 2 Composite as a Conformal and High‐Sensitivity Piezoelectric Sensor

Abstract Additively manufactured flexible and high‐performance piezoelectric devices are highly desirable for sensing and energy harvesting of 3D conformal structures. Herein, the study reports a significantly enhanced piezoelectricity in polyvinylidene fluoride (PVDF) achieved through the in situ dipole alignment of PVDF within PVDF‐2D molybdenum disulfide (2D MoS 2 ) composite by 3D printing. The shear stress‐induced dipole poling of PVDF and 2D MoS 2 alignment are harnessed during 3D printing to boost piezoelectricity without requiring a post‐poling process. The results show a remarkable, more than the eight‐fold increment in the piezoelectric coefficient ( d 33 ) for 3D printed PVDF‐8wt.% MoS 2 composite over cast neat PVDF. The underlying mechanism of piezoelectric property enhancement is attributed to the increased volume fraction of β phase in PVDF, filler fraction, heterogeneous strain distribution around PVDF‐MoS 2 interfaces, and strain transfer to the nanofillers as confirmed by microstructural analysis and finite element simulation. These results provide a promising route to design and fabricate high‐performance 3D piezoelectric devices via 3D printing for next‐generation sensors and mechanical–electronic conformal devices.

2D MoS2↗