Natural and anthropogenic infrasound may travel vast distances, making it an invaluable resource for monitoring phenomena such as nuclear explosions, volcanic eruptions, severe storms, and many others. Typically, these waves are captured using pressure sensors, which cannot encode the direction of arrival—critical information when the source location is not known beforehand. Obtaining this information therefore requires arrays of sensors with apertures ranging from tens of meters to kilometers depending on the wavelengths of interest. This is often impractical in locations that lack the necessary real estate (urban areas, rugged regions, or remote islands); in any case, it requires multiple power, digitizer, and telemetry deployments. Here, the theoretical basis behind a compact infrasound direction of arrival sensor based on the acoustic metamaterials is presented. This sensor occupies a footprint that is orders of magnitude smaller than the span of a typical infrasound array. The diminutive size of the unit greatly expands the locations where it can be deployed. The sensor design is described, its ability to determine the direction of arrival is evaluated, and further avenues of study are suggested.
Sierra/SD provides a massively parallel implementation of structural dynamics finite element analysis, required for high-fidelity, validated models used in modal, vibration, static and shock analysis of weapons systems. This document provides a user’s guide to the input for Sierra/SD. Details of input specifications for the different solution types, output options, element types and parameters are included. The appendices contain detailed examples, and instructions for running the software on parallel platforms.
Sierra/SD provides a massively parallel implementation of structural dynamics finite element analysis, required for high fidelity, validated models used in modal, vibration, static and shock analysis of structural systems. This manual describes the theory behind many of the constructs in Sierra/SD. For a more detailed description of how to use Sierra/SD, we refer the reader to User’s Manual. Many of the constructs in Sierra/SD are pulled directly from published material. Where possible, these materials are referenced herein. However, certain functions in Sierra/SD are specific to our implementation. We try to be far more complete in those areas. The theory manual was developed from several sources including general notes, a programmer_notes manual, the user’s notes and of course the material in the open literature.
This document presents tests from the Sierra Structural Mechanics verification test suite. Each of these tests is run nightly with the Sierra/SD code suite and the results of the test checked versus the correct analytic result. For each of the tests presented in this document the test setup, derivation of the analytic solution, and comparison of the Sierra/SD code results to the analytic solution is provided. This document can be used to confirm that a given code capability is verified or referenced as a compilation of example problems.
Sierra/SD is an engineering structural dynamics code that provides Sandia and other customers a tool to model structural and acoustic physics on large complex physical systems using massively parallel processing. This report provides a detailed overview on Sierra/SD’s most recent physics package: coupled electro-mechanical physics. This capability uses the finite element method to model coupled electro-mechanical physics exhibited by piezoelectric materials. This report provides an applications overview, theory overview, and verification examples demonstrating the electro-mechanical physics modeling capabilities of Sierra/SD.
The “how to” document guides the user through complicated aspects of software usage. It should supplement both the User’s manual and the Theory document, by providing examples and detailed discussion that reduce learning time for complex set ups. These documents are intended to be used together. We will not formally list all parameters for an input here – see the User’s manual for this. All the examples in the “How To” document are part of the Sierra/SD test suite, and each will run with no modification. The nature of this document casts together a number of rather unrelated procedures. Grouping them is difficult. Please try to use the table of contents and the index as a guide in finding the analyses of interest.
This paper presents a topology optimization formulation for frequency-domain dynamics to reduce solution dependence upon initial guess and considered loading conditions. Due to resonance phenomena in undamped steady-state dynamics, objectives measuring dynamic response possess many local minima that may represent poor solutions to a design problem, an issue exacerbated for design with respect to multiple frequencies. We propose an extension of the modified error-in-constitutive-equations (MECE) method, used previously in material identification inverse problems, as a new approach for frequency-domain dynamics topology optimization to mitigate these issues. The main idea of the proposed framework is to incorporate an additional penalty-like term in the objective function that measures the discrepancy in the constitutive relations between stresses and strains and between inertial forces and displacements. Then, the design problem is cast within a PDE-constrained optimization formulation in which we seek displacements, stresses, inertial forces, and a density-field solution that minimize our new objective subject to conservation of linear momentum plus some additional constraints. We show that this approach yields superior designs to conventional gradient-based optimization approaches that solely use a functional of displacements as the objective, while strictly enforcing the constitutive equations. The MECE strategy integrates into a density-based topology optimization scheme for void–solid or two-phase material structural design. We highlight the merits of our approach in a variety of scenarios for direct frequency response design, considering multiple frequency load cases and structural objectives.
We develop a generalized stress inversion technique (or the generalized inversion method) capable of recovering stresses in linear elastic bodies subjected to arbitrary cuts. Specifically, given a set of displacement measurements found experimentally from digital image correlation (DIC), we formulate a stress estimation inverse problem as a partial differential equation-constrained optimization problem. We use gradient-based optimization methods, and we accordingly derive the necessary gradient and Hessian information in a matrix-free form to allow for parallel, large-scale operations. By using a combination of finite elements, DIC, and a matrix-free optimization framework, the generalized inversion method can be used on any arbitrary geometry, provided that the DIC camera can view a sufficient part of the surface. We present numerical simulations and experiments, and we demonstrate that the generalized inversion method can be applied to estimate residual stress.
While elastic metasurfaces offer a remarkable and very effective approach to the subwavelength control of stress waves, their use in practical applications is severely hindered by intrinsically narrow band performance. In applications to electromagnetic and photonic metamaterials, some success in extending the operating dynamic range was obtained by using nonlocality. However, while electronic properties in natural materials can show significant nonlocal effects, even at the macroscales, in mechanics, nonlocality is a higher-order effect that becomes appreciable only at the microscales. This study introduces the concept of intentional nonlocality as a fundamental mechanism to design passive elastic metasurfaces capable of an exceptionally broadband operating range. The nonlocal behavior is achieved by exploiting nonlocal forces, conceptually akin to long-range interactions in nonlocal material microstructures, between subsets of resonant unit cells forming the metasurface. These long-range forces are obtained via carefully crafted flexible elements, whose specific geometry and local dynamics are designed to create remarkably complex transfer functions between multiple units. The resulting nonlocal coupling forces enable achieving phase-gradient profiles that are functions of the wavenumber of the incident wave. The identification of relevant design parameters and the assessment of their impact on performance are explored via a combination of semianalytical and numerical models. The nonlocal metasurface concept is tested, both numerically and experimentally, by embedding a total-internal-reflection design in a thin-plate waveguide. Results confirm the feasibility of the intentionally nonlocal design concept and its ability to achieve a fully passive and broadband wave control.
The “how to” document is designed to help walk the analyst through difficult aspects of software usage. It should supplement both the User’s manual and the Theory document, by providing examples and detailed discussion that reduce learning time for complex set ups. These documents are intended to be used together. We will not formally list all parameters for an input here – see the User’s manual for this. All the examples in the “How To” document are part of the Sierra/SD test suite, and each will run with no modification. The nature of this document casts together a number of rather unrelated procedures. Grouping them is difficult. Please try to use the table of contents and the index as a guide in finding the analyses of interest.
This document presents tests from the Sierra Structural Mechanics verification test suite. Each of these tests is run nightly with the Sierra/SD code suite and the results of the test checked versus the correct analytic result. For each of the tests presented in this document the test setup, derivation of the analytic solution, and comparison of the Sierra/SD code results to the analytic solution is provided. This document can be used to confirm that a given code capability is verified or referenced as a compilation of example problems.
This work explores how High Perforrnance Computing is enabling acoustic solutions across a wide-range ofscience and engineering applications that were historically intractable.