Search NASA⌕ Search

SEARCH · Search NASA

Results for “COUPLED STRESS”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Simulation of Multiphase Flow and Poromechanical Effects Around Injection Wells in CO 2 Storage Sites

In geological CO 2 storage operations, wellbore deformations and leakage pathways formations can occur around injection and abandoned wells subjected to high rates and long-term CO 2 injection. To guide engineering design and prevent CO 2 leakage risks, a full understanding of the underlying physics and robust numerical models is necessary to evaluate the response of underground formations in the near wellbore region and in the reservoir. In this study, a multi-scale and multi-physics open-source simulator (GEOS) is used to simulate multiphase flow and poromechanical deformations over time in three dimensions. The governing equations for mechanical deformations of the rock body and multiphase compositional fluid flow within the rock matrix are solved with a fully coupled finite element and finite volume approach. The Drucker–Prager model with friction hardening is applied to simulate elastoplastic deformation and a multiphase fluid model with power-law correlations for relative permeability is used to model the migration of CO 2 plume, which are coupled with numerical implicit scheme. Simulation results are verified against multiple analytical solutions for multiphase flow and wellbore problems, thus demonstrating the accuracy of this advanced simulator. In two engineering applications, here we highlight the impact of elastoplastic deformation and coupled modeling for assessing induced displacements and stress perturbations, which are more pronounced in the near wellbore regions. This work focuses on short-term processes in the vicinity of injection wells where stress evolutions, rock deformations and multiphase compositional flow and transport are simulated jointly to ensure wellbore stability and prevent damage. This fully coupled geomechanical model can simulate multiphase flow and any associated poromechanical effects within the CO 2 storage site and in the surrounding formations. Such a large-scale, long-term, multi-physics simulation model is useful in many ways: it can guide operational decisions for CO 2 injection, assess the containment potential and risks of a site, and analyze the wellbore stability and integrity during and after CO 2 injection.

58 GEOSCIENCES↗
NASA NTRSDate not supplied

Structural Design and Aeromechanical Analysis of Unconventional Blades for Future Mars Rotorcraft

The structural design for rotor blades with thin, unconventional airfoils for Earth-based testing is obtained using three-Dimensional (3D) aeromechanical analysis. This analysis is in support of the NASA ROAMX project. The outer mold was provided by NASA, but the internal structural design was developed at the University of Maryland and is presented here, along with a thorough aeromechanical analysis. The main objectives are to verify the structural integrity of the design and understand the unique aeroelastic behavior of the non-conventional airfoils designed for low-Reynolds number and high-Mach number flow. Six different blade models are considered, with the pitch axis varied from quarter-chord to mid-chord to determine the effect of C.G. offset on natural frequencies, blade deformations, root loads, and 3D stresses. All six blade designs are first studied under rotation in pure vacuum, and significant reductions in root loads and 3D stresses are achieved by moving the pitch axis closer to mid-chord to reduce the C.G. offset. Based on the vacuum analysis, the blade design with the pitch axis at 40% chord is selected for aerodynamic analysis. The blade control load, airloads, deformations, and 3D stresses are studied for steady hover. Dynamic control load and dynamic 3D stresses are studied for unsteady hover achieved using cyclic. Significant elastic twist is observed due to the trapeze effect and propeller moment, affecting the spanwise distribution of aerodynamic loads on the blades. However, the effect of aerodynamic loads on 3D stresses and the blade control load is found to be minimum. The dynamic control load is found to increase significantly due to inertial coupling from the C.G. offset. The dynamic stresses also increase, although still have factors of safety greater than two for both tensile and compressive stress.

Structural Design↗
NASA NTRS1991

Models of recurrent strike-slip earthquake cycles and the state of crustal stress

Numerical models of the strike-slip earthquake cycle, assuming a viscoelastic asthenosphere coupling model, are examined. The time-dependent simulations incorporate a stress-driven fault, which leads to tectonic stress fields and earthquake recurrence histories that are mutually consistent. Single-fault simulations with constant far-field plate motion lead to a nearly periodic earthquake cycle and a distinctive spatial distribution of crustal shear stress. The predicted stress distribution includes a local minimum in stress at depths less than typical seismogenic depths. The width of this stress 'trough' depends on the magnitude of crustal stress relative to asthenospheric drag stresses. The models further predict a local near-fault stress maximum at greater depths, sustained by the cyclic transfer of strain from the elastic crust to the ductile asthenosphere. Models incorporating both low-stress and high-stress fault strength assumptions are examined, under Newtonian and non-Newtonian rheology assumptions. Model results suggest a preference for low-stress (a shear stress level of about 10 MPa) fault models, in agreement with previous estimates based on heat flow measurements and other stress indicators.

Lyzenga, Gregory A.↗
NASA NTRS2015

Multiscale and Multiphysics Modeling of Additive Manufacturing of Advanced Materials

The objective of this proposed project is to research and develop a prediction tool for advanced additive manufacturing (AAM) processes for advanced materials and develop experimental methods to provide fundamental properties and establish validation data. Aircraft structures and engines demand materials that are stronger, useable at much higher temperatures, provide less acoustic transmission, and enable more aeroelastic tailoring than those currently used. Significant improvements in properties can only be achieved by processing the materials under nonequilibrium conditions, such as AAM processes. AAM processes encompass a class of processes that use a focused heat source to create a melt pool on a substrate. Examples include Electron Beam Freeform Fabrication and Direct Metal Deposition. These types of additive processes enable fabrication of parts directly from CAD drawings. To achieve the desired material properties and geometries of the final structure, assessing the impact of process parameters and predicting optimized conditions with numerical modeling as an effective prediction tool is necessary. The targets for the processing are multiple and at different spatial scales, and the physical phenomena associated occur in multiphysics and multiscale. In this project, the research work has been developed to model AAM processes in a multiscale and multiphysics approach. A macroscale model was developed to investigate the residual stresses and distortion in AAM processes. A sequentially coupled, thermomechanical, finite element model was developed and validated experimentally. The results showed the temperature distribution, residual stress, and deformation within the formed deposits and substrates. A mesoscale model was developed to include heat transfer, phase change with mushy zone, incompressible free surface flow, solute redistribution, and surface tension. Because of excessive computing time needed, a parallel computing approach was also tested. In addition, after investigating various methods, a Smoothed Particle Hydrodynamics Model (SPH Model) was developed to model wire feeding process. Its computational efficiency and simple architecture makes it more robust and flexible than other models. More research on material properties may be needed to realistically model the AAM processes. A microscale model was developed to investigate heterogeneous nucleation, dendritic grain growth, epitaxial growth of columnar grains, columnar-to-equiaxed transition, grain transport in melt, and other properties. The orientations of the columnar grains were almost perpendicular to the laser motion's direction. Compared to the similar studies in the literature, the multiple grain morphology modeling result is in the same order of magnitude as optical morphologies in the experiment. Experimental work was conducted to validate different models. An infrared camera was incorporated as a process monitoring and validating tool to identify the solidus and mushy zones during deposition. The images were successfully processed to identify these regions. This research project has investigated multiscale and multiphysics of the complex AAM processes thus leading to advanced understanding of these processes. The project has also developed several modeling tools and experimental validation tools that will be very critical in the future of AAM process qualification and certification.

Liou, Frank↗
NASA NTRS1997

Initialization and Predictability of a Coupled ENSO Forecast Model

The skill of a coupled ocean-atmosphere model in predicting ENSO has recently been improved using a new initialization procedure in which initial conditions are obtained from the coupled model, nudged toward observations of wind stress. The previous procedure involved direct insertion of wind stress observations, ignoring model feedback from ocean to atmosphere. The success of the new scheme is attributed to its explicit consideration of ocean-atmosphere coupling and the associated reduction of "initialization shock" and random noise. The so-called spring predictability barrier is eliminated, suggesting that such a barrier is not intrinsic to the real climate system. Initial attempts to generalize the nudging procedure to include SST were not successful; possible explanations are offered. In all experiments forecast skill is found to be much higher for the 1980s than for the 1970s and 1990s, suggesting decadal variations in predictability.

Chen, Dake↗
DOE OSTI2021

Stepwise Dynamic Calibration of a Hydromechanical Simulation Using Time-Lapse Vertical Seismic Profile

This study aims to develop a methodology for calibrating subsurface stress changes through time-lapse Vertical Seismic Profiling (VSP) integration. The selected study site is the 13-10A injector well within the ongoing CO2-EOR operation of the Farnsworth Field Unit. The Time-lapse VSP dataset carries the combined effects of fluid substitution and mean effective stress changes, thereby providing a dataset amenable for the calibration of production and injection-induced stress changes. The concept is similar to calibrating a reservoir simulation model in that the process honor real field data to set up an inverse problem. The solution optimizes the independent and impactful geomechanical parameters that replicate the observed time-lapse seismic velocity changes. This stress calibration is enabled by 4D geomechanical modeling and the VSP Integration workflow. This calibration benefits from extensive geological, geophysical and geomechanical characterization through 3D seismic data, geophysical well logs, and core assessed as part of the 1D MEM conducted on the 13-10A subject well. These data are used to develop a site-specific rock physics model. The Biot Gassmann workflow combines rock physics and reservoir simulation outputs to determine the fluid substitution contribution to seismic velocity change. Additionally, modeled seismic velocity attributed to mean effective stress are determined from the geomechanical simulation outputs, and the stress-velocity relationship developed from the ultrasonic seismic velocity measurements on the extracted Morrow B core. A penalty function is then formed between the modeled seismic velocities and the observed time-lapse VSP dataset. Four independent and impactful geomechanical parameters have been determined. These are the bulk modulus and shear modulus for zero porosity and the shear and compressional seismic velocity to mean effective stress derivatives. The dataset of numerous coupled hydromechanical- geomechanical simulation realizations is built by combining variations of the four stated geomechanical parameters. A machine learning-assisted workflow comprised of an artificial neural network and a particle swarm optimizer are used to converge on the optimal geomechanical parameters. The successful execution of this workflow has affirmed the suitability of acoustic time-lapse measurements for 4D-VSP geomechanical stress calibration pending measurable stress sensitivities within the anticipated effective stress changes and the availability of suitable and reliable datasets for petroelastic modeling.

02 PETROLEUM↗
DOE OSTI2020

The effects of stress, temperature and facet structure on growth of { 10 1 ¯ 2 } twins in Mg: a molecular dynamics and phase field study

The present study focuses on the roles of structure, stress and temperature on the mobility of facets bounding twin domains and on the overall impact on the growth kinetics of { 10 1 ¯ 2 } twins in Mg. We first use molecular dynamics simulations (MD) to study the relationship between interface structures, stress, temperature and mobility. As facets can exhibit distinct atomic scale structures, we also quantify the effect of stress relaxation mediated by misfit dislocations on the mobility of the basal-prismatic interface. The information gathered for ten different facets is used to calibrate an anisotropic phase field model predicting the effects of temperature and stress on the kinetics and shape of a growing { 10 1 ¯ 2 } 3D twin, by comparison with MD data. The comparative analysis suggests that, as compared to simulations of individual facets, the effective mobility of facets can be significantly different during growth of a 3D twin. This can be attributed to the complex three-dimensional internal stress state of the twin and to coupling in facet motion triggered by facet junctions. As such, by directly fitting both the twin growth kinetics and morphology as predicted by MD for different stresses and temperatures, we can identify a limited set of rate-limiting facets that control twin growth and morphology, and measure their effective mobility.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Predicting Microstructural Void Nucleation in Discontinuous Fiber Composites through Coupled in-situ X-ray Tomography Experiments and Simulations

Composite materials have become widely used in engineering applications, in order to reduce the overall weight of structures while retaining their required strength. In this work, a composite material consisting of discontinuous glass fibers in a polypropylene matrix is studied at the microstructural level through coupled experiments and simulations, in order to uncover the mechanisms that cause damage to initiate in the microstructure under macroscopic tension. Specifically, we show how hydrostatic stresses in the matrix can be used as a metric to explain and predict the exact location of microvoid nucleation that occurs during damage initiation within the composite’s microstructure. Furthermore, this work provides evidence that hydrostatic stresses in the matrix can lead to coupled microvoid nucleation and early fiber breakage, and that small fragments of fibers can play an important role in the process of microvoid nucleation. These results significantly improve our understanding of the mechanics that drive the initiation of damage in the complex microstructures of discontinuous fiber reinforced thermoplastics, while also allowing scientists and engineers to predict the microstructural damage behavior of these composites at sub-fiber resolution and with high accuracy.

36 MATERIALS SCIENCE↗
DOE OSTI2023

TReactMech v4.217

TReactMech couples geomechanical processes (poroelasticity, failure, and inelastic strain) with multiphase nonisothermal flow (derived from TOUGH2) and reactive geochemical transport. At its core is the reactive-transport code TOUGHREACT v4.13. TReactMech is an efficient hybrid parallel simulator, solving the geomechanics using finite elements and MPI/PetSc, the multiphase flow using integrated finite difference and MPI/PETSc, and the reactive chemistry using OpenMP. The advantages of TReactMech are in its multiphase flow capabilities (e.g., supercritical CO2, supercritical water, air) and parallel geomechanics including full 3-D stress tensor, shear and tensile failure, coupled to porosity and permeability changes. It is backwardly compatible with TOUGH2 and TOUGHREACT v4.13, allowing for easier transitions between the codes. TReactMech can be used to simulate many natural and engineered subsurface systems, including geothermal reservoirs, borehole heat exchangers, geologic carbon sequestration, geologic storage of nuclear waste, groundwater resources, weathering, sediment diagenesis, seafloor hydrothermal circulation, hydrofracturing in unconventional reservoirs, and injection/production-induced surface deformation.

Sonnenthal, Eric↗
NASA NTRS2000

Carotid baroreflex responsiveness in heat-stressed humans

The effects of whole body heating on human baroreflex function are relatively unknown. The purpose of this project was to identify whether whole body heating reduces the maximal slope of the carotid baroreflex. In 12 subjects, carotid-vasomotor and carotid-cardiac baroreflex responsiveness were assessed in normothermia and during whole body heating. Whole body heating increased sublingual temperature (from 36.4 +/- 0.1 to 37.4 +/- 0.1 degrees C, P < 0.01) and increased heart rate (from 59 +/- 3 to 83 +/- 3 beats/min, P < 0. 01), whereas mean arterial blood pressure (MAP) was slightly decreased (from 88 +/- 2 to 83 +/- 2 mmHg, P < 0.01). Carotid-vasomotor and carotid-cardiac responsiveness were assessed by identifying the maximal gain of MAP and heart rate to R wave-triggered changes in carotid sinus transmural pressure. Whole body heating significantly decreased the responsiveness of the carotid-vasomotor baroreflex (from -0.20 +/- 0.02 to -0.13 +/- 0.02 mmHg/mmHg, P < 0.01) without altering the responsiveness of the carotid-cardiac baroreflex (from -0.40 +/- 0.05 to -0.36 +/- 0.02 beats x min(-1) x mmHg(-1), P = 0.21). Carotid-vasomotor and carotid-cardiac baroreflex curves were shifted downward and upward, respectively, to accommodate the decrease in blood pressure and increase in heart rate that accompanied the heat stress. Moreover, the operating point of the carotid-cardiac baroreflex was shifted closer to threshold (P = 0.02) by the heat stress. Reduced carotid-vasomotor baroreflex responsiveness, coupled with a reduction in the functional reserve for the carotid baroreflex to increase heart rate during a hypotensive challenge, may contribute to increased susceptibility to orthostatic intolerance during a heat stress.

NASA Discipline Cardiopulmonary↗
NASA NTRS2010

Experimental Validation of FE/BEM Dynamic Strain Model Under Diffuse Acoustic Field Loading

Structural finite element (FE) models naturally output displacement or acceleration response data. However, they can also be used to compute stress, internal forces, and strain response. When coupled with a boundary element model (BEM) of the fluid surrounding the structure, a fully coupled analysis can be performed. Modeling a diffuse acoustic field in the BEM fluid provides an excitation like that found when the structure is placed in a reverberation chamber. Fully coupling the structural FE model to the acoustic BEM model provides a means to predict not only the acceleration response of the panel to diffuse field loading, but also the ability to predict the dynamic stress and strain response. This type of model has been available with current predictive tools, but experimental validation of the prediction of dynamic stress or strain is difficult to find. An aluminum panel was instrumented with accelerometers and strain gages and hung in a reverberation room and subjected to a diffuse acoustic field. This paper presents the comparison of the experimental and predicted results.

Tsoi, W. Ben↗
NASA NTRS1999

Tensegrity and mechanoregulation: from skeleton to cytoskeleton

OBJECTIVE: To elucidate how mechanical stresses that are applied to the whole organism are transmitted to individual cells and transduced into a biochemical response. DESIGN: In this article, we describe fundamental design principles that are used to stabilize the musculoskeletal system at many different size scales and show that these design features are embodied in one particular form of architecture that is known as tensegrity. RESULTS: Tensegrity structures are characterized by use of continuous tension and local compression; architecture, prestress (internal stress prior to application of external force), and triangulation play the most critical roles in terms of determining their mechanical stability. In living organisms, use of a hierarchy of tensegrity networks both optimizes structural efficiency and provides a mechanism to mechanically couple the parts with the whole: mechanical stresses applied at the macroscale result in structural rearrangements at the cell and molecular level. CONCLUSION: Due to use of tensegrity architecture, mechanical stress is concentrated and focused on signal transducing molecules that physically associate with cell surface molecules that anchor cells to extracellular matrix, such as integrins, and with load-bearing elements within the internal cytoskeleton and nucleus. Mechanochemical transduction may then proceed through local stress-dependent changes in molecular mechanics, thermodynamics, and kinetics within the cell. In this manner, the entire cellular response to stress may be orchestrated and tuned by altering the prestress in the cell, just as changing muscular tone can alter mechanical stability and structural coordination throughout the whole musculoskeletal system.

NASA Discipline Cell Biology↗
DOE OSTI2023

Graded microstructure and mechanical properties of spark plasma sintered Fe-Cr alloys

Developing graded microstructure and mechanical properties is critical for accelerating the design and optimization of structural materials for a wide range of applications. Here, in this study, Fe-Cr alloys were fabricated by spark plasma sintering (SPS) technique. The microstructure, microhardness, and tensile properties of the as-fabricated and thermally annealed variants were investigated. Graded grain structures were created in the as-fabricated cylinders along both axial and radial directions. Grain size was gradually reduced from the sample periphery to the center. Microhardness measurements reveal a gradual decrease of hardness towards the periphery of the samples. The as-fabricated Fe-Cr alloys show a desired combination of tensile strength and elongation, primarily due to the formation of a high density of oxides, voids, dislocations, and grain boundaries. After in-situ thermal annealing at 600 °C under 60 MPa for 2 h using SPS, the Fe-Cr alloys underwent minimal grain growth, and the graded grain structures were retained. The hardness was more uniformly distributed in the annealed variants, and the tensile strength was reduced with an increase in the total elongation, which is attributed to the dissolution of nano-sized oxide particles and the relief of the residual stress. This study demonstrates that SPS coupled with subsequent heat treatment can tailor the graded microstructure and control the mechanical properties of Fe-Cr alloys, showing potential applications in other alloy systems.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Osmotic control of the spacing of parallel shear cracks in shale growing subcritically in geologic past

The geological genesis of natural cracks in sedimentary rocks such as shale is a problem that needs to be understood to improve the technology of hydraulic fracturing as well as deep sequestration of harmful fluids. Why are the vertical natural cracks roughly parallel and equidistant, and why is the spacing roughly 10 cm rather than 1 cm or 100 cm? Fracture mechanics of critical cracks cannot answer this question. Neither can the material heterogeneity. The growth of critical parallel cracks is impossible because the relative crack face displacements would immediately localize into one crack, leading to an earthquake. The cracks must have formed, on the tectonic time scale, by a slow growth of subcritical shear cracks governed by the Charles-Evans law. The idea advanced here is that what controls the crack spacing is the balance between the reduction, due to shear dilatancy, of the concentration of ions such as Na + and Cl - in each fracture process zone (PFZ), which decelerates the cracks, and the restoration of ion concentration by diffusion of ions from the space between the cracks into the FPZ. This diffusion of water is driven mainly by the osmotic pressure gradient, which offsets the deceleration and depends strongly on the crack spacing. A simple analytical solution of the steady state is rendered possible by approximating the ion concentration profiles between adjacent cracks by parabolic arcs. Applying this theory to Woodford shale yields the approximate crack spacing of 10 cm, which is realistic. Furthermore, the stability of unlimited parallel mode II frictional crack growth is proven by examining the second variation of the free energy. Water concentration drop in the FPZ due to shear dilatancy and its restoration by water diffusion from the inter-crack space have similar effect, although probably much weaker.

42 ENGINEERING↗
DOE OSTI2021

Updated Thermofluid Performance of the Simplified Flat Variant of the HEMJ

Our group has recently developed and studied “finger”-type divertors that are a simplified version of the helium-cooled modular divertor with multiple jets (HEMJ) using coupled computational fluid dynamics and thermal stress simulations. Such a simplified geometry could reduce complexity and cost given the large number of fingers required to cover the total divertor target area. Previous experimental studies for this simplified flat design reported lower heat transfer coefficients and higher pressure drops than the HEMJ, contrary to numerical predictions. Subsequent measurements determined that the original test section had significant dimensional variations in the jet exit holes. A new test section was therefore manufactured and tested in the Georgia Tech (GT) helium loop. The experimental results presented here for this test section at maximum heat flux of 7.1 MW/m 2 are in good agreement with numerical predictions. Correlations developed from these experimental data are extrapolated to predict the maximum heat flux that can be accommodated by the flat design and the coolant pumping power requirements under prototypical conditions. Lastly, numerical simulations are used to estimate the sensitivity of the flat design to geometric variations typical of manufacturing tolerances and variations in the gap width.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Thermo-Mechanical Distortion of Tungsten-Coated Steel During High Heat Flux Testing Using Plasma Arc Lamps

An experimental setup and a test section were designed and fabricated for high heat flux testing (HHFT) of neutron-irradiated specimens using water-wall plasma arc lamps. Because of the radiological considerations and limitations of reactor irradiation, the size of the test articles was limited to disks less than 10 mm in diameter. The specimen was clamped onto an actively cooled block, and clamping allowed the insertion of several thermocouples on the back surface of the specimen through a copper (Cu) block. Five vacuum plasma sprayed tungsten (W)–coated F82H steel specimens were subjected to HHFT. Surface profilometry measurements, which were conducted after HHFT, revealed central bowing of the top W surface. This type of residual distortion occurred for all of the specimens, and the larger the specimens were, the larger was the distortion.In an attempt to understand specimen distortion and address the science questions related to the testing of subsize specimens during HHFT, a simplified thermo-mechanical model was developed. By using a measured temperature in the Cu as an isothermal boundary condition, the model eliminated the need for coupling cooling fluid flow models with stress models, greatly simplifying the analysis. The main variable in the proposed model is hC, i.e., the thermal contact conductance between the F82H and the Cu washer. Inelastic properties, including hardening properties, were considered for F82H steel and Cu. Numerical simulation results demonstrated a buildup of residual deformation during HHFT and a very complex state of stress and deformation during typical heat flux (HF) cycling. Additionally, hoop stress evolution during a high heat flux cycle reveals that F82H at an interface with W would be mainly in compression during HF application and experienced a transition to a tension state during cooldown. Also, specimen distortion evolves during each HF cycle, as the specimen bows downward during HF application and upward during the cooldown period between HF cycles. The final specimen distortion, i.e., upward bowing of the specimen center, was qualitatively predicted for hC values of 4000 to 5000 W/(m 2 ·K). This hC range of values, for which bulging is obtained, is at the lower spectrum of the range of values for hC, consistent with the low thermal contact conductance expected from the unpolished F82H surface.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Tardigrade, Version 1.1x

Tardigrade is a collection of functions, libraries, and stand-alone code which enables the simulation of materials using the higher-order micromorphic framework. The tools present are the interface to INL's MOOSE code which is used as the underlying FEA solver, several constitutive equations (linear elasticity, a pressure-sensitive elasto-plastic model, crystal plasticity, and additional models for geo-materials) useful for demonstration purposes and the simulation of some materials, the micromorphic filter which converts DNS information to micromorphic stress-deformation quantities, and the overlap coupling framework which enables multi-scale simulations of a coupled higher-order macroscale and a classically defined microscale.

Miller, Nathan↗