Search NASA⌕ Search

SEARCH · Search NASA

Results for “Multiscale simulation”

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 991 records · Page 55

Ion-Scale Secondary Flux Ropes Generated by Magnetopause Reconnection as Resolved by MMS

New Magnetospheric Multiscale (MMS) observations of small-scale (approx. 7 ion inertial length radius) flux transfer events (FTEs) at the dayside magnetopause are reported. The 1O km MMS tetrahedron size enables their structure and properties to be calculated using a variety of multispacecraft techniques, allowing them to be identified as flux ropes, whose flux content is small (approx. 22 kWb).The current density, calculated using plasma and magnetic field measurements independently, is found to be filamentary. lntercomparison of the plasma moments with electric and magnetic field measurements reveals structured non-frozen-in ion behavior. The data are further compared with a particle-in-cell simulation. It is concluded that these small-scale flux ropes, which are not seen to be growing, represent a distinct class of FTE which is generated on the magnetopause by secondary reconnection.

Eastwood, J. P.↗

Advancing Organized Convection Representation in the Unified Model: Implementing and Enhancing Multiscale Coherent Structure Parameterization

To address the effect of stratiform latent heating on meso- to large-scale circulations, an enhanced implementation of the Multiscale Coherent Structure Parameterization (MCSP) is developed for the Met Office Unified Model. MCSP represents the top-heavy stratiform latent heating from under-resolved organized convection in general circulation models. We couple the MCSP with a mass-flux convection scheme (CoMorph-A) to improve storm lifecycle continuity. The improved MCSP trigger is specifically designed for mixed-phase deep convective cloud, combined with a background vertical wind shear, both known to be crucial for stratiform development. We also test a cloud top temperature dependent convective-stratiform heating partitioning, in contrast to the earlier fixed partitioning. Assessments from ensemble weather forecasts and decadal simulations demonstrate that MCSP directly reduces cloud deepening and precipitation areas by moderating mesoscale circulations. Indirectly, it amends tropical precipitation biases, notably correcting dry and wet biases over India and the Indian Ocean, respectively. Remarkably, the scheme outperforms a climate model ensemble by improving seasonal precipitation cycle predictions in these regions. The scheme also improves Madden-Julian Oscillation (MJO) spectra, achieving better alignment with observational and reanalysis data by intensifying the simulated MJO over the Indian Ocean during phases 4 to 5. However, the scheme increases precipitation overestimation over the Western Pacific. Shifting from fixed to temperature-dependent convective-stratiform partitioning reduces the Pacific precipitation overestimation and further improves the seasonal cycle in India. Spatially correlated biases highlight the necessity for advances beyond deterministic approaches to align MCSP with environmental conditions.

54 ENVIRONMENTAL SCIENCES↗

MAS: A versatile Landau-fluid eigenvalue code for plasma stability analysis in general geometry

We have developed a new global eigenvalue code, Multiscale Analysis for plasma Stabilities (MAS), for studying plasma problems with wave toroidal mode number (n) and frequency (ω) in a broad range of interest in general tokamak geometry, based on a five-field Landau-fluid description of thermal plasmas. Beyond keeping the necessary plasma fluid response, we further retain the important kinetic effects including diamagnetic drift, ion finite Larmor radius, finite parallel electric field (E||), ion and electron Landau resonances in a self-consistent and non-perturbative manner without sacrificing the attractive efficiency in computation. The physical capabilities of the code are evaluated and examined in the aspects of both theory and simulation. In theory, the comprehensive Landau-fluid model implemented in MAS can be reduced to the well-known ideal MHD model, electrostatic ion-fluid model, and drift-kinetic model in various limits, which clearly delineates the physics validity regime. In simulation, MAS has been well benchmarked with theory and other gyrokinetic and kinetic-MHD hybrid codes in a manner of adopting the unified physical and numerical framework, which covers the kinetic Alfv\'en wave (KAW), ion sound wave (ISW), low-n kink, high-n ion temperature gradient mode (ITG) and kinetic ballooning mode (KBM). Moreover, MAS is successfully applied to model the Alfv\'en eigenmode (AE) activities in DIII-D discharge #159243, which faithfully captures the frequency sweeping of reversed shear Alfv\'en eigenmode (RSAE), the tunneling damping of toroidal Alfv\'en eigenmode (TAE), as well as the polarization characteristics of kinetic beta-induced Alfv\'en eigenmode (KBAE) and beta-induced Alfv\'en-acoustic eigenmode (BAAE) being consistent with former gyrokinetic theory and simulation. With respect to the key progress contributed to the community, MAS has the advantage of combining rich physics ingredients, realistic global geometry and high computation efficiency together for plasma stability analysis in linear regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

MMS Extended Mission Design: Evaluation of a Lunar Gravity Assist Option

This paper will describe a study that was carried out on the design of a set of maneuvers that were considered for the later stages of an extended mission of the Magnetospheric Multiscale (MMS) mission. The goal of these maneuvers was to put MMS into a significantly different orbit from those flown heretofore, so allowing science collection in a different region of the magnetosphere. This study was made feasible by the fact that the rate at which fuel is being consumed to maintain small formations on the MMS high-apogee orbit is less than expected pre-flight: the current consumption rate is only about 2 kg/yr/spacecraft. In addition, the spacecraft finished the prime mission with a significant amount of fuel remaining: this was about 1-sigma above the mean when compared with pre-launch Monte Carlo simulations. The resulting situation is similar to that of a libration orbit mission, where station-keeping requires so little fuel that any margin at all will lead to an extensive mission lifetime. In the case of MMS, the spacecraft could, if desired, perform formation flying in the current orbit for several decades. Alternatively, the spacecraft could use a significant fraction of the remaining fuel to perform major orbit modifications, while still leaving enough to conduct formation flying for on the order of a decade. The extended mission maneuvers studied here are further apogee-raises, with the goal of setting up one or more lunar gravity assists. Geometry dictates that a lunar encounter is only achievable when the MMS apogee vector lies approximately in the lunar orbit plane: this limits the possible dates to mid-2021 or early 2027.

Design↗

Mechanism of Reconnection on Kinetic Scales Based on Magnetospheric Multiscale Mission Observations

We examine the role that ions and electrons play in reconnection using observations from the Magnetospheric Multiscale (MMS) mission on kinetic ion and electron scales, which are much shorter than magnetohydrodynamic scales. This study reports observations with unprecedented high resolution that MMS provides for magnetic field (7.8 ms) and plasma (30 ms for electrons and 150 ms for ions). We analyze and compare approaches to the magnetopause in 2016 November, to the electron diffusion region in the magnetotail in 2017 July followed by a current sheet crossing in 2018 July. Besides magnetic field reversals, changes in the direction of the flow velocity, and ion and electron heating, MMS observed large fluctuations in the electron flow speeds in the magnetotail. As expected from numerical simulations, we have verified that when the field lines and plasma become decoupled a large reconnecting electric field related to the Hall current (1–10 mV/m) is responsible for fast reconnection in the ion diffusion region. Although inertial accelerating forces remain moderate (1–2 mV/m), the electric fields resulting from the divergence of the full electron pressure tensor provide the main contribution to the generalized Ohm’s law at the neutral sheet (as large as 200 mV/m). In our view, this illustrates that when ions decouple electron physics dominates. The results obtained on kinetic scales may be useful for better understanding the physical mechanisms governing reconnection processes in various magnetized laboratory and space plasmas.

W. M. Macek↗

Quantum-Accurate Multiscale Modeling of Shock Hugoniots, Ramp Compression Paths, Structural and Magnetic Phase Transitions, and Transport Properties in Highly Compressed Metals

Fully characterizing high energy density (HED) phenomena using pulsed power facilities (Z machine) and coherent light sources is possible only with complementary numerical modeling for design, diagnostic development, and data interpretation. The exercise of creating numerical tests, that match experimental conditions, builds critical insight that is crucial for the development of a strong fundamental understanding of the physics behind HED phenomena and for the design of next generation pulsed power facilities. The persistence of electron correlation in HED materials arising from Coulomb interactions and the Pauli exclusion principle is one of the greatest challenges for accurate numerical modeling and has hitherto impeded our ability to model HED phenomena across multiple length and time scales at sufficient accuracy. An exemplar is a ferromagnetic material like iron, while familiar and widely used, we lack a simulation capability to characterize the interplay of structure and magnetic effects that govern material strength, kinetics of phase transitions and other transport properties. Herein we construct and demonstrate the Molecular-Spin Dynamics (MSD) simulation capability for iron from ambient to earth core conditions, all software advances are open source and presently available for broad usage. These methods are multi-scale in nature, direct comparisons between high fidelity density functional theory (DFT) and linear-scaling MSD simulations is done throughout this work, with advancements made to MSD allowing for electronic structure changes being reflected in classical dynamics. Main takeaways for the project include insight into the role of magnetic spins on mechanical properties and thermal conductivity, development of accurate interatomic potentials paired with spin Hamiltonians, and characterization of the high pressure melt boundary that is of critical importance to planetary modeling efforts.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Long-Range Hierarchical Nanocrystal Assembly Driven by Molecular Structural Transformation

According to the article, the hierarchical control in biogenic minerals, from precise nanomorphology control to subsequent macroscopic assembly, remains a formidable challenge in artificial synthesis. Studies in biomineralization, however, are largely limited to atomic andmolecular scale crystallization, devoting little attention to biomolecular higher-order structures (HOSs) which critically impact long-range assembly of biominerals. Here we demonstrate a biomimetic route and quantitative simulations that explore peptide HOSs on guiding nanocrystal formation and anisotropic assembly into hierarchical structures. It is found that the Pt{100}-specific peptide T7 (Ac-TLTTLTN-CONH2) adopts ST-turn secondary structure, promoting cubic Pt nanocrystal formation at low concentration, and spontaneously transforms into a β-sheet with increased concentration. The β-sheet T7-Pt{100} specificity drives cubic Pt nanocrystals to self-assemble into large-area, long-range, [100] linear assemblies. This study provides a robust demonstration for bio/nonbiogenic material specificity, nanoscale synthesis, and long-range self-organization with biomolecular HOSs and opens vast opportunities for multiscale programmable structures.

Chemistry↗

Impact of Tropical Cyclone Wind Forcing on the Global Climate in a Fully Coupled Climate Model

Tropical cyclones (TCs) alter upper-ocean temperature and influence ocean heat content via enhanced turbulent mixing. A better understanding of the role of TCs within the climate system requires a fully coupled modeling framework, where TC-induced ocean responses feed back to the atmosphere and subsequently to the climate mean state and variability. Here, we investigate the impacts of TC wind forcing on the global ocean and the associated feedbacks within the climate system using the fully coupled Community Earth System Model version 1.3 (CESM1.3). Using the low-resolution version of CESM1.3 (1° atmosphere and ocean grid spacing) with no intrinsic TCs, we conduct a suite of sensitivity experiments by inserting TC winds extracted from a high-resolution (0.25° atmosphere grid spacing) TC-permitting simulation into the low-resolution model. Results from the low-resolution TC experiment are compared to a low-resolution control simulation to diagnose TCs’ impact. We found that the added TC winds can increase ocean heat content by affecting ocean vertical mixing, air–sea enthalpy fluxes, and cloud amount. Here, the added TCs can influence mean SST, precipitation, ocean subsurface temperature, and ocean mixed layer depth. We found a strengthening of the wind-driven subtropical cells and a weakening of the Atlantic meridional overturning circulation due to the changes of surface buoyancy fluxes. TCs in the model cause anomalous equatorward ocean heat convergence in the deep tropics and an increase of poleward ocean heat transport out of the subtropics. Our modeling results provide new insights into the multiscale interactions between TCs and the coupled climate system.

54 ENVIRONMENTAL SCIENCES↗

Micro-Mechanically Guided High-Throughput Alloy Design Exploration Towards Metastability-Induced H Embrittlement Resistance

We develop a high-throughput approach for studying H embrittlement (HE)-resistance in alloys, which is based on combinatorial compositional screening of metastability effects by in situ scanning electron microscopy H-analyses. The project objective included: (i) technique development of high-throughput screening (HTS) for HE-resistance, (ii) discovery of new metallic materials with superior HE-resistance, (iii) Multiscale verification of HE-resistance of the new alloys and H-barrier layers, from atomic scale to an engineering scale. The investigation focused on a model system that enables exploration of different metastable states in Fe-based complex-concentrated alloys. Composition spread islands with hundreds of varying compositions are fabricated on a single substrate using a combinatorial co-sputtering technique. To screen these alloys, we developed and employed a home-built SEM-based integrated analysis system capable of characterizing H-permeability, H-trapping, H-influence on mechanical properties, and other material properties, as well as atomistic to continuum simulations to study the underlying physics.

08 HYDROGEN↗

New Developments in the Embedded Statistical Coupling Method: Atomistic/Continuum Crack Propagation

A concurrent multiscale modeling methodology that embeds a molecular dynamics (MD) region within a finite element (FEM) domain has been enhanced. The concurrent MD-FEM coupling methodology uses statistical averaging of the deformation of the atomistic MD domain to provide interface displacement boundary conditions to the surrounding continuum FEM region, which, in turn, generates interface reaction forces that are applied as piecewise constant traction boundary conditions to the MD domain. The enhancement is based on the addition of molecular dynamics-based cohesive zone model (CZM) elements near the MD-FEM interface. The CZM elements are a continuum interpretation of the traction-displacement relationships taken from MD simulations using Cohesive Zone Volume Elements (CZVE). The addition of CZM elements to the concurrent MD-FEM analysis provides a consistent set of atomistically-based cohesive properties within the finite element region near the growing crack. Another set of CZVEs are then used to extract revised CZM relationships from the enhanced embedded statistical coupling method (ESCM) simulation of an edge crack under uniaxial loading.

Saether, E.↗

Multiscale Electricity Modeling for Evaluating Carbon Capture and Sequestration Technologies (Final Report)

Carbon capture and sequestration (CCS) technologies that can operate with a high degree of operating flexibility could provide necessary electric grid flexibility in a system with high shares of variable renewables. This project examines the deployment and dispatch potential of twelve unique flexible CCS (FLECCS) technologies that encompass post-combustion carbon dioxide (CO 2 ) capture designs, concepts using a storage media to enable energy arbitrage, and hybrid processes that integrate CCS with direct air capture (DAC) for flexibility with net zero or negative CO 2 emissions. FLECCS technology potential is explored with a multi-model, multi-scale framework including the Regional Energy Deployment System (ReEDS) electric sector capacity expansion model (CEM) and the PLEXOS production cost model (PCM). Innovative methods were developed to represent FLECCS technology operating modes, performance, and cost in the two models. ReEDS was then used to simulate nine scenarios for each FLECCS technology, three CO 2 emissions price futures reaching $\$$150, $\$$225, and $\$$300/tCO 2 in 2050; and three scenarios for FLECCS technology deployment favorability relative to competing technologies. For each CO 2 price and reference FLECCS favorability, the 2050 infrastructures from ReEDS model are downscaled and implemented in PLEXOS to examine hourly dispatch under detailed operational constraints that are not included in ReEDS. FLECCS technologies exhibited a wide range of deployment potential ranging from none to several hundred gigawatts of capacity, with outcomes highly sensitive to input cost and performance parameters that are inherently highly uncertain. When deployed, FLECCS tended to displace a combination of wind, solar, and natural gas-based technologies rather than supporting increased renewable deployment. As a result, CO 2 emissions reductions facilitated by FLECCS deployment tended to come with higher overall system costs and electricity prices. When economically competitive, FLECCS technologies can contribute significant flexible generation and firm capacity to the grid, but continued technology development and an expanded analytical scope are necessary to fully understand FLECCS deployment potential its impact on the electric power sector. Follow-on analysis incorporating captured CO 2 tax credit value from the Inflation Reduction Act (IRA) and other potential policy scenarios could be particularly valuable, as this policy can substantially change the relative competitiveness of FLECCS technologies.

03 NATURAL GAS↗

Designing Physicochemically‐Ordered Interphases for High‐Performance Composites

To enhance the mechanical properties of carbon fiber‐reinforced polymer composites, a physicochemical scaffold is designed incorporating microscopically architected chemically reactive nanofibers that act as a multiscale bridge between the carbon fibers and the matrix. Thermally activated nanofibers leverage their morphologically driven mechanochemical properties to form covalent bonds with adjacent polymer molecules, creating a co‐continuous network that dramatically enhances fiber‐matrix load transfer. By meticulously controlling the nanofiber architecture through variable surface area, functional group availability, and polymer chain alignment effects, the extent of covalent bonding between nanofibers and the matrix is manipulated ultimately resulting in improved carbon fiber‐matrix adhesion. Further, the concept was validated using polyacrylonitrile nanofibers within an acrylonitrile butadiene styrene matrix in a discontinuous carbon fiber‐reinforced composite system. Nanomechanical studies using atomic force microscopy and low‐field nuclear magnetic resonance spectroscopy confirmed immobilized, chemically transferred, and ordered nanostructures at the interphase. The resulting composites demonstrated ≈56% and ≈175% improvements in tensile strength and toughness, respectively, compared to composites without nanofiber. Comprehensive thermal, rheological, and X‐ray scattering analyses, along side all‐atomic molecular dynamics simulations, revealed the fundamental mechanisms behind these improvements in mechanical behavior. The versatility and efficacy of the approach have the potential to address longstanding interphase challenges in the composite industry.

36 MATERIALS SCIENCE↗

Denoising diffusion algorithm for inverse design of microstructures with fine-tuned nonlinear material properties

Here we introduce a denoising diffusion algorithm to discover microstructures with nonlinear fine-tuned properties. Denoising diffusion probabilistic models are generative models that use diffusion-based dynamics to gradually denoise images and generate realistic synthetic samples. By learning the reverse of a Markov diffusion process, we design an artificial intelligence to efficiently manipulate the topology of microstructures to generate a massive number of prototypes that exhibit constitutive responses sufficiently close to designated nonlinear constitutive behaviors. To identify the subset of microcstructures with sufficiently precise fine-tuned properties, a convolutional neural network surrogate is trained to replace high-fidelity finite element simulations to filter out prototypes outside the admissible range. Results of this study indicate that the denoising diffusion process is capable of creating microstructures of fine-tuned nonlinear material properties within the latent space of the training data. More importantly, this denoising diffusion algorithm can be easily extended to incorporate additional topological and geometric modifications by introducing high-dimensional structures embedded in the latent space. Numerical experiments are conducted on the open-source mechanical MNIST data set (Lejeune, 2020). Consequently, this algorithm is not only capable of performing inverse design of nonlinear effective media, but also learns the nonlinear structure–property map to quantitatively understand the multiscale interplay among the geometry, topology, and their effective macroscopic properties.

42 ENGINEERING↗

Unveiling the catalytic mechanism of GTP hydrolysis in microtubules

Microtubules (MTs) are large cytoskeletal polymers, composed of αβ-tubulin heterodimers, capable of stochastically converting from polymerizing to depolymerizing states and vice versa. Depolymerization is coupled with hydrolysis of guanosine triphosphate (GTP) within β-tubulin. Hydrolysis is favored in the MT lattice compared to a free heterodimer with an experimentally observed rate increase of 500- to 700-fold, corresponding to an energetic barrier lowering of 3.8 to 4.0 kcal/mol. Mutagenesis studies have implicated α-tubulin residues, α:E254 and α:D251, as catalytic residues completing the β-tubulin active site of the lower heterodimer in the MT lattice. The mechanism for GTP hydrolysis in the free heterodimer, however, is not understood. Additionally, there has been debate concerning whether the GTP-state lattice is expanded or compacted relative to the GDP state and whether a “compacted” GDP-state lattice is required for hydrolysis. In this work, extensive quantum mechanics/molecular mechanics simulations with transition-tempered metadynamics free-energy sampling of compacted and expanded interdimer complexes, as well as a free heterodimer, have been carried out to provide clear insight into the GTP hydrolysis mechanism. α:E254 was found to be the catalytic residue in a compacted lattice, while in the expanded lattice, disruption of a key salt bridge interaction renders α:E254 less effective. The simulations reveal a barrier decrease of 3.8 ± 0.5 kcal/mol for the compacted lattice compared to a free heterodimer, in good agreement with experimental kinetic measurements. Additionally, the expanded lattice barrier was found to be 6.3 ± 0.5 kcal/mol higher than compacted, demonstrating that GTP hydrolysis is variable with lattice state and slower at the MT tip.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Machine-learning modeling of magnetization dynamics in quasi-equilibrium and driven metallic spin systems

Here, we present a perspective on recent progress in machine-learning (ML) force-field approaches for large-scale Landau–Lifshitz–Gilbert (LLG) simulations of metallic spin systems. Building on a generalization of the Behler–Parrinello (BP) architecture originally developed for quantum molecular dynamics, we develop scalable and transferable ML models that faithfully capture the complex, environment-dependent electron-mediated exchange fields characteristic of itinerant magnets. A central ingredient of this framework is the implementation of symmetry-aware magnetic descriptors based on group-theoretical bispectrum formalisms. Leveraging these ML force fields, LLG simulations faithfully reproduce hallmark non-collinear magnetic orders—such as the 120° and tetrahedral states—on the triangular lattice, and successfully capture the complex spin textures emerging in the mixed-phase states of a square-lattice double-exchange model under thermal quench. We further discuss a generalized potential theory that extends the BP formalism to incorporate both conservative and nonconservative electronic torques, thereby enabling ML models to learn nonequilibrium exchange fields from computationally demanding microscopic approaches such as nonequilibrium Green’s-function techniques. This extension yields quantitatively accurate predictions of voltage-driven domain-wall motion and establishes a foundation for quantum-accurate, multiscale modeling of nonequilibrium spin dynamics and spintronic functionalities.

Descriptors↗

Hierarchical Multiscale Process Modeling of A Textile Composite Y-Joint for the Aurora D8 Aircraft

Woven polymer matrix composites (PMCs) are widely used in aerospace applications for their favorable specific properties. However, the behavior of the inter- and intra-tow matrix during manufacturing of woven PMCs is difficult to isolate experimentally due to geometric and processing complexities. Classical closed-form homogenization approaches allow for composite thermos-chemo-mechanical properties to be characterized at any fiber volume fraction. In this work, virtual curing is applied to the inter- and intra-tow matrix materials of a plain weave repeating unit cell (RUC). Thermo-chemo-mechanical property evolution within the tows is calculated using the Concentric Cylinder Model (CCM) and the Rule of Mixtures (ROM). Curing is implemented through user-written subroutines within the commercial finite element method software Abaqus. Virtual curing results are compared with those from a linear-elastic thermal cool down simulation in Abaqus. Results indicate that tow property evolution during curing has a non-negligible effect on the final stress state.

curing↗

Influences of shock imprinting on mix in a 3D-printed porous media

Mixing of materials in porous media can cause a significant impact on fusion yield as previously demonstrated by the National Ignition Facility (NIF) MARBLE Campaign. Initially, the reactants are separated, with deuterium in the lattice struts and a tritium gas fill in the voids. Lattice parameters such as the strut thickness and relative pitch, provide a control for the mix parameters in the experiment. Los Alamos National Laboratory’s (LANL) BOSQUE project looks to better understand how the mix of the reactants and shell materials impact the fusion burn and resultant yield on various laser platforms. xRAGE’s Eulerian hydrodynamics and adaptive mesh refinement (AMR) provide the unique ability to study the impacts of multiscale features of complex lattice structures. This modeling provides the ability to measure shock front variations as the wave progress’ through a given media. Initial conditions of the lattice are essential to accurately model mix and burn measured by experiment. By varying the initial orientation and densities of these lattice regions the early time dynamics of how the shock is launched into the system is changed and advocates for the study of resulting effects. In this work, we will study the sensitivities of shock effects in varying 3D printed geometric systems and how these shocks alter the structure and mix in the lattice. We will discuss both preliminary experimental results and simulations to help plan and constrain future experiments where we will study the impact of different lattice geometries and lattice bulk densities. This work concludes with the relative impacts of lattice geometries on shock speeds at different bulk densities and the resultant mix widths due to those shock interactions. We see agreement with theory at the higher end of our bulk density study, however, as we approach lower bulk density systems the dynamics of these interactions begin to change.

2 Photon Polymerization↗

Driving macro-scale transformations in three-dimensional-printed biopolymers through controlled induction of molecular anisotropy at the nanoscale

Motivated by the need to harness the properties of renewable and biodegradable polymers for the design and manufacturing of multi-scale structures with complex geometries, we have employed our additive manufacturing platform that leverages molecular self-assembly for the production of metre-scale structures characterized by complex geometries and heterogeneous material composition. As a precursor material, we used chitosan, a chemically modified form of chitin, an abundant and sustainable structural polysaccharide. We demonstrate the ability to control concentration-dependent crystallization as well as the induction of the preferred orientation of the polymer chains through the combination of extrusion-based robotic fabrication and directional toolpathing. Anisotropy is demonstrated and assessed through high-resolution micro-X-ray diffraction in conjunction with finite element simulations. Using this approach, we can leverage controlled and user-defined small-scale propagation of residual stresses to induce large-scale folding of the resulting structures.

36 MATERIALS SCIENCE↗