Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “classical compression”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

An engineering model of woven composites based on micromechanics

Composites with three-dimensional woven architectures exhibit large strains to failure when compared to composites made up of the same materials but not with three-dimensional interlocking tows. The fracture mechanics of such three-dimensional architectures is a subject requiring substantial investigation and experimental testing. Classical fracture mechanics concepts (for instance, an isolated defect in a homogeneous body) will not be applicable to the woven fracture test specimen. The use of an isolated singularity to characterize an entire specimen is inadequate when the density of defects is considerable and the material is heterogeneous. Modelling of such a complex system requires a great deal of insight and consideration as well as prudent choices of model sizes to make numerical schemes feasible. The purpose of this manuscript is to review our recently acquired knowledge of damage accumulation in woven composites and to describe a practicable model of the macroscopic behavior in these and other complex composite architectures based on such knowledge. In this manuscript, discussion will be limited to uniaxial compressive loading; considerations of general loading (monotonic and cyclic) will appear in a subsequent manuscript. Our modelling efforts may be briefly described as follows: the composite is subdivided into microstructural elements (microelements) in which the micromechanical modelling is either understood rigorously or can be represented adequately by statistical parameters. There can be microstructural elements for many different types of composite components, such as the various types of warp and weft and matrix for three-dimensional woven composites. The physical dimensions of microelements are made as large as possible while the response within the element can still be represented by a single micromechatlical calculation. The various elements are linked together(sometimes by associating distinct corners and edges, sometimes by superposition) in a pattern which resembles a particular weave architecture. The model can then be loaded in any manner and the linear and nonlinear elastic responses of representative weaves can be calculated. After the elastic regime, the fracture response is determined by monitoring the damage accumulation.

Carter, W. C.↗
DOE OSTI2020

Impact of Laboratory-Accelerated Aging Methods to Study Alkali–Silica Reaction and Reinforcement Corrosion on the Properties of Concrete

This study focuses on two separate investigations of the main aging mechanisms: alkali–silica reactivity (ASR) and the corrosion of reinforcing steel (rebar) concrete, both of which may result in a premature failure to meet the serviceability or strength requirements of a concrete structure. However, these processes occur very slowly, spanning decades. The impact of direct chemical additives to fresh concrete to accelerate ASR and the corrosion of reinforcing steel on the fresh and hardened properties of the ensuing material are investigated to inform the potential use of chemicals in large-scale studies. The deterioration of reinforced concrete (RC) is determined by means of expansion, cracking, bulk diffusivity and surface resistivity measurements, and compressive, split tensile and flexural strength tests. The results indicate that the addition of sodium hydroxide and calcium chloride can effectively accelerate the crack formation and propagation in concrete due to ASR and the corrosion of rebar, respectively. The ASR-induced cracks maintained a constant crack width from 0.05 mm to 0.1 mm over the measurement period regardless of the intensity of aging acceleration. Adding 4% chloride by weight of cement for accelerating rebar corrosion resulted in an average crack that was 82% larger than in the case of ASR accelerated with the addition of sodium hydroxide. The addition of alkali resulted in an increase in early-age (7-day) strength. At a total alkali loading of 2.98 kg/m3, 3.84 kg/m3 and 5.57 kg/m3, the 28-day compressive strength of concrete decreased by 3%, 10% and 24%, respectively. Similarly, a higher early-age strength and a lower later-age strength was observed for the concrete in the presence of corrosive calcium chloride. The results from this research are expected to inform future studies on the long-term performance of RC structures under accelerated ASR and corrosion.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A Tutorial on Bayesian analysis of linear shock compression data

Gas gun and other shock compression experiments often produce shock wave velocity measurements that are linearly associated with particle velocity. Traditionally, this empirical relationship is quantified with a single Hugoniot curve that is estimated using least squares regression. However, for downstream modeling and simulation tasks, it is often more useful to have multiple Hugoniot curves in the pressure–volume plane that are consistent with the data. We employ Bayesian uncertainty quantification methods as a framework for propagating measurement uncertainty through to model parameters and predictions. Specifically, this Tutorial shows how to sample multiple Hugoniot curves in the pressure–volume plane that are consistent with the shock wave-particle velocity measurements in a two-step Bayesian approach. First, we obtain an analytical expression for the posterior distribution of the linear model parameters using Bayesian linear regression. Second, we propagate samples from the posterior distribution through the Rankine–Hugoniot equations to yield Hugoniot curves in the pressure–volume plane. The procedure is demonstrated with publicly available data on argon, copper, and nickel, and compared against bootstrapping and linear regression. The Bayesian procedure is shown to be interpretable, computationally inexpensive, and less sensitive than an alternative bootstrapping approach to the removal of the point in the copper dataset that has the largest particle velocity. As a Tutorial on Bayesian methodology for the shock compression community, we provide several derivations and explanations that make this paper self-contained, and make all code and data available at github.com/llnl/BALSCD.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2010

Loop Shaping Control Design for a Supersonic Propulsion System Model Using Quantitative Feedback Theory (QFT) Specifications and Bounds

This paper covers the propulsion system component modeling and controls development of an integrated mixed compression inlet and turbojet engine that will be used for an overall vehicle Aero-Propulso-Servo-Elastic (APSE) model. Using previously created nonlinear component-level propulsion system models, a linear integrated propulsion system model and loop shaping control design have been developed. The design includes both inlet normal shock position control and jet engine rotor speed control for a potential supersonic commercial transport. A preliminary investigation of the impacts of the aero-elastic effects on the incoming flow field to the propulsion system are discussed, however, the focus here is on developing a methodology for the propulsion controls design that prevents unstart in the inlet and minimizes the thrust oscillation experienced by the vehicle. Quantitative Feedback Theory (QFT) specifications and bounds, and aspects of classical loop shaping are used in the control design process. Model uncertainty is incorporated in the design to address possible error in the system identification mapping of the nonlinear component models into the integrated linear model.

Connolly, Joseph W.↗
DOE OSTI2023

Supplementary Material for ORNL_AISD-Ex

This dataset provides supplementary material for the previously published dataset ORNL_AISD-Ex (1), which is available at the following website: https://www.osti.gov/biblio/1907919 The dates comprises one compressed folder called ornl_aisd_ex.zip. The compressed folder ornl_aisd_ex.zip contains 1,000 CSV files, each of them titled ornl_aisd_ex_ID.csv, where ID is a number that ranges between 1 and 1,000. The information contained in ornl_aisd_ex_ID.csv corresponds to the information of the molecules compresses inside the ornl_aisd_ex_ID.tar.gz of the dataset ORNL_AISD-Ex. Each row in each .CSV file is associated with a molecule, and the columns contain the following information: 1) molecules ID 2) SMILES string representation 3) DFTB-PE (eV): formation energy 4) first 50 electronic excitation modes 5) oscillators strengths of the first 50 electronic excitation modes The compression of this information into CSV files will allow a more agile extraction and management of information to the users that do not have access to large scale HPC platforms. REFERENCES (1) Lupo Pasini, Massimiliano, Mehta, Kshitij, Yoo, Pilsun, and Irle, Stephan. ORNL_AISD-Ex: Quantum chemical prediction of UV/Vis absorption spectra for over 10 million organic molecules. United States: N. p., 2023. Web. doi:10.13139/OLCF/1907919.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

A Tensor Network-Based Quantum Algorithm for the Nonlinear 1D Burgers' Equation

In this work, we implement a tensor network-based quantum algorithm to solve unsteady, nonlinear partial differential equations (PDEs). The challenge lies in how to effectively represent, encode, process, and evolve the nonlinear system of PDEs on quantum computers. We will discuss the new techniques using the compressible 1-dimensional (1D) Burgers' equation as an example, because it represents the fundamental nonlinear feature and yet removes certain complexity in physics, allowing us to focus on the design of quantum algorithms. Previous attempts to solve nonlinear PDEs in quantum computation have often involved storing multiple copies of solutions or employing linearizations. Neither is practical due to exponential scaling with evolution time or insufficient solution accuracy. Our framework is based on matrix product states (MPSs) and matrix product operators (MPOs). For example, the velocity field is represented by MPS, whereas the linear and nonlinear spatial differential terms of the velocity field are processed by MPOs. Our primary focus herein is to verify and validate the various tensor network components of the algorithm using solutions obtained by the classical algorithms on high performance computing (HPC) architectures. We use a classical time marching method to demonstrate the functionality of the tensor network operations to model the PDE and their robustness with the time evolution of the system. Our classical simulation results demonstrate the utility of tensor network-based operations in modeling nonlinear PDEs and highlight the necessity as well as potential advantages of using quantum simulations for these techniques.

Gopalakrishnan Meena, Murali [ORNL] (ORCID:0000000↗
DOE OSTI2020

Thermal conductivity and its relation to atomic structure for symmetrical tilt grain boundaries in silicon

We perform a systematic study of thermal resistance and conductance of tilt grain boundaries (GBs) in Si using classical molecular dynamics. The GBs studied are naturally divided into three groups according to the structural units forming the GB core. In this work we find that, within each group, the GB thermal conductivity strongly correlates with the excess GB energy. All three groups predict nearly the same GB conductivity extrapolated to the high-energy limit. This limiting value is close to the thermal conductivity of amorphous Si, suggesting similar heat transport mechanisms. While the lattice thermal conductivity decreases with temperature, the GB conductivity slightly increases. However, at high temperatures it turns over and starts decreasing if the GB structure undergoes a premelting transformation. Analysis of vibrational spectra of GBs resolved along different directions sheds light on the mechanisms of their thermal resistance. The existence of alternating tensile and compressive atomic environments in the GB core gives rise to localized vibrational modes, frequency gaps creating acoustic mismatch with lattice phonons, and anharmonic vibrations of loosely bound atoms residing in open atomic environments.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Modeling of shock-induced force on an isolated particle in water and air

The prediction of force on an isolated particle, while a shock is passing over it, is an important problem in many natural and industrial applications. Although the flow monotonically changes from the pre-shock to the post-shock state, the particle's force has been observed to behave nonmonotonically with a sharp peak when the shock is located halfway across the particle. This nonmonotonic behavior is due to the unsteady nature of the compression and rarefaction waves that radiate as the shock diffracts around the particle and, therefore, cannot be predicted by a quasi-steady model. An accurate force model must account for the unsteady nature of the flow and the sharp discontinues in the flow properties across the shock. In this work, we test four different inviscid models and observe that the compressible Maxey–Riley–Gatignol (C-MRG) model is the most accurate based on comparison with results from particle-resolved inviscid simulations at two different Mach numbers for both water and air as the medium. The C-MRG model is first demonstrated to predict the force on a stationary particle accurately and then extended to capture the force on a moving particle. Numerical complexities regarding the implementation of the C-MRG model are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1999

Development of Curved-Plate Elements for the Exact Buckling Analysis of Composite Plate Assemblies Including Transverse-Shear Effects

The analytical formulation of curved-plate non-linear equilibrium equations including transverse-shear-deformation effects is presented. A unified set of non-linear strains that contains terms from both physical and tensorial strain measures is used. Linearized, perturbed equilibrium equations (stability equations) that describe the response of the plate just after buckling occurs are derived. These equations are then modified to allow the plate reference surface to be located a distance z(sub c) from the centroidal surface. The implementation of the new theory into the VICONOPT exact buckling and vibration analysis and optimum design computer program is described. The terms of the plate stiffness matrix using both classical plate theory (CPT) and first-order shear-deformation plate theory (SDPT) are presented. The effects of in-plane transverse and in-plane shear loads are included in the in-plane stability equations. Numerical results for several example problems with different loading states are presented. Comparisons of analyses using both physical and tensorial strain measures as well as CPT and SDPT are made. The computational effort required by the new analysis is compared to that of the analysis currently in the VICONOPT program. The effects of including terms related to in-plane transverse and in-plane shear loadings in the in-plane stability equations are also examined. Finally, results of a design-optimization study of two different cylindrical shells subject to uniform axial compression are presented.

McGowan, David M.↗
DOE OSTI2024

Nonequilibrium formulation of varying-temperature bit erasure

Landauer's principle states that erasing a bit of information at fixed temperature T costs at least $k$ B $T$ ln 2 units of work. Here we investigate erasure at varying temperature, to which Landauer's result does not apply. Here we formulate bit erasure as a stochastic nonequilibrium process involving a compression of configuration space, with physical and logical states associated in a symmetric way. Erasure starts and ends at temperature T, but temperature can otherwise vary with time in an arbitrary way. Defined in this way, erasure is governed by a set of nonequilibrium fluctuation relations that show that varying-temperature erasure can done with less work than $k$ B $T$ ln 2. As a result, erasure and the complementary process of bit randomization can be combined to form a work-producing engine cycle.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1999

Nanomechanics of Carbon and CxByNz Nanotubes: Via a Quantum Molecular Dynamics Method

Nanomechanics of single-wall C, BN and BC$_3$ and B doped C nanotubes under axial compression and tension are investigated through a generalized tight-binding molecular dynamics (GTBMD) and {\it ab-initio} electronic structure methods. The dynamic strength of BN, BC$_3$ and B doped C nanotubes for small axial strain are comparable to each other. The main difference is in the critical strain at which structural collapse occurs. For example, even a shallow doping with B lowers the value of critical strain for C nanotubes. The critical strain for BN nanotube is found to be more than that for the similar C nanotube. Once the structural collapse starts to occur we find that carbon nanotubes irreversibly go into plastic deformation regime via the formation of tetrahedral (four-fold coordinated) bonds at the location of sharp pinches or kinks. This finding is considerably different from the classical MD (molecular dynamics) simulation results known so far. The energetics and electronic densities of states of the collapsed structures, investigated with {\it ab-initio) methods, will also be discussed.

Srivastava, Deepak↗
DOE OSTI2023

Towards critical and supercritical electromagnetic fields

Abstract The availability of ever stronger, laser-generated electromagnetic fields underpins continuing progress in the study and application of nonlinear phenomena in basic physical systems, ranging from molecules and atoms to relativistic plasmas and quantum electrodynamics. This raises the question: how far will we be able to go with future lasers? One exciting prospect is the attainment of field strengths approaching the Schwinger critical field ${E}_{\mathrm{cr}}$ in the laboratory frame, such that the field invariant ${E}^2-{c}^2{B}^2>{E}_{\mathrm{cr}}^2$ is reached. The feasibility of doing so has been questioned, on the basis that cascade generation of dense electron–positron plasma would inevitably lead to absorption or screening of the incident light. Here we discuss the potential for future lasers to overcome such obstacles, by combining the concept of multiple colliding laser pulses with that of frequency upshifting via a tailored laser–plasma interaction. This compresses the electromagnetic field energy into a region of nanometre size and attosecond duration, which increases the field magnitude at fixed power but also suppresses pair cascades. Our results indicate that laser facilities with peak power of tens of PW could be capable of reaching ${E}_{\mathrm{cr}}$ . Such a scenario opens up prospects for the experimental investigation of phenomena previously considered to occur only in the most extreme environments in the universe.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Computational investigation of water glasses using machine-learning potentials

The molecular origins of water’s anomalous properties have long been a subject of scientific inquiry. The liquid–liquid phase transition hypothesis, which posits the existence of distinct low-density and high-density liquid states separated by a first-order phase transition terminating at a critical point, has gained increasing experimental and computational support and offers a thermodynamically consistent framework for many of water’s anomalies. However, experimental challenges in avoiding crystallization near the postulated liquid–liquid critical point have focused attention to water’s canonical glassy states: low-density and high-density amorphous ice. Here, we use two Deep Potential machine-learning models, trained on the Strongly Constrained and Appropriately Normed density functional and the highly accurate Many-Body Polarizable potential, to conduct an investigation of water’s glassy phenomenology based on quantum mechanical calculations. Despite not being explicitly trained on amorphous ices, both models accurately capture the structure and transformation of the water glasses, including their interconversion along different thermodynamic paths. Isobaric quenching of liquid water at various pressures generates a continuum of intermediate amorphous ices and density fluctuations increase near the liquid–liquid critical pressure. The glass transition temperatures of the amorphous ices produced at different pressures exhibit two distinct branches, corresponding to low-density and high-density amorphous ice behaviors, consistent with experiment and the liquid–liquid transition hypothesis. Extrapolating transformation pressures from isothermal compressions to experimental compression rates brings our simulations into excellent agreement with data. Our findings demonstrate that machine-learning potentials trained on equilibrium phases can effectively model nonequilibrium glassy behavior and pave the way for studying long-timescale, out-of-equilibrium processes with quantum mechanical accuracy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

FSISPH: An SPH formulation for impacts between dissimilar materials

Here, we present an SPH formulation with several new features designed to better model the fully-compressible interaction of dissimilar materials. We developed the new method to simulate the atmospheric entry and break-up of small celestial bodies in planetary atmospheres. The formulation uses a unity-based, density-energy discretization of the hydrodynamic conservation laws with linear-corrected kernel gradients. To account for variations in compressibility, we use an HLLC approximate Riemann solver to adjust the velocity gradient at material interfaces. To handle large transverse velocity discontinuities, we introduce a simple slip interface model that limits the artificial viscosity at material interfaces. Diffusion is optionally applied through the velocity gradient and this allows the density and specific thermal energy to evolve in a manner more consistent with the first law of thermodynamics in comparison to other more direct diffusion schemes. We also introduce a material-local second-order artificial conduction scheme used to smooth the specific thermal energy field. Material damage fits neatly under this framework by treating the damage front as a material interface. The method has been implemented as a solver, FSISPH, within the code, Spheral++, and is publicly available on github. We test our new solver on a number of classic shock, mixing, and multi-material problem. The components we outline can significantly improve accuracy of SPH for problems with sharp contact discontinuities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Dimensionality Reduction with Variational Encoders Based on Subsystem Purification

Efficient methods for encoding and compression are likely to pave the way toward the problem of efficient trainability on higher-dimensional Hilbert spaces, overcoming issues of barren plateaus. Here, we propose an alternative approach to variational autoencoders to reduce the dimensionality of states represented in higher dimensional Hilbert spaces. To this end, we build a variational algorithm-based autoencoder circuit that takes as input a dataset and optimizes the parameters of a Parameterized Quantum Circuit (PQC) ansatz to produce an output state that can be represented as a tensor product of two subsystems by minimizing $Tr(ρ^2)$. The output of this circuit is passed through a series of controlled swap gates and measurements to output a state with half the number of qubits while retaining the features of the starting state in the same spirit as any dimension-reduction technique used in classical algorithms. The output obtained is used for supervised learning to guarantee the working of the encoding procedure thus developed. We make use of the Bars and Stripes (BAS) dataset for an 8 × 8 grid to create efficient encoding states and report a classification accuracy of 95% on the same. Thus, the demonstrated example provides proof for the working of the method in reducing states represented in large Hilbert spaces while maintaining the features required for any further machine learning algorithm that follows.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Variational quantum reinforcement learning via evolutionary optimization

Abstract Recent advances in classical reinforcement learning (RL) and quantum computation point to a promising direction for performing RL on a quantum computer. However, potential applications in quantum RL are limited by the number of qubits available in modern quantum devices. Here, we present two frameworks for deep quantum RL tasks using gradient-free evolutionary optimization. First, we apply the amplitude encoding scheme to the Cart-Pole problem, where we demonstrate the quantum advantage of parameter saving using amplitude encoding. Second, we propose a hybrid framework where the quantum RL agents are equipped with a hybrid tensor network-variational quantum circuit (TN-VQC) architecture to handle inputs of dimensions exceeding the number of qubits. This allows us to perform quantum RL in the MiniGrid environment with 147-dimensional inputs. The hybrid TN-VQC architecture provides a natural way to perform efficient compression of the input dimension, enabling further quantum RL applications on noisy intermediate-scale quantum devices.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

MechElastic: A Python library for analysis of mechanical and elastic properties of bulk and 2D materials

We report the MechElastic Python package evaluates the mechanical and elastic properties of bulk and 2D materials using the elastic coefficient matrix ( C ij ) obtained from any ab-initio density-functional theory (DFT) code. The current version of this package reads the output of VASP, ABINIT, and Quantum Espresso codes (but it can be easily generalized to any other DFT code) and performs the appropriate post-processing of elastic constants as per the requirement of the user. This program can also detect the input structure's crystal symmetry and test the mechanical stability of all crystal classes using the Born-Huang criteria. Various useful material-specific properties such as elastic moduli, longitudinal and transverse elastic wave velocities, Debye temperature, elastic anisotropy, 2D layer modulus, hardness, Pugh's ratio, Cauchy's pressure, Kleinman's parameter, and Lame's coefficients, can be estimated using this program. Another existing feature of this program is to employ the ELATE package (2016) [29] and plot the spatial variation of several elastic properties such as Poisson's ratio, linear compressibility, shear modulus, and Young's modulus in three dimensions. Further, the MechElastic package can plot the equation of state (EOS) curves for energy and pressure for a variety of EOS models such as Murnaghan, Birch, Birch-Murnaghan, and Vinet, by reading the inputted energy/pressure versus volume data obtained via numerical calculations or experiments. This package is particularly useful for the high-throughput analysis of elastic and mechanical properties of materials.

2D materials↗
NASA NTRS1990

Instabilities and subharmonic resonances of subsonic heated round jets, volume 2

When a jet is perturbed by a periodic excitation of suitable frequency, a large-scale coherent structure develops and grows in amplitude as it propagates downstream. The structure eventually rolls up into vortices at some downstream location. The wavy flow associated with the roll-up of a coherent structure is approximated by a parallel mean flow and a small, spatially periodic, axisymmetric wave whose phase velocity and mode shape are given by classical (primary) stability theory. The periodic wave acts as a parametric excitation in the differential equations governing the secondary instability of a subharmonic disturbance. The (resonant) conditions for which the periodic flow can strongly destabilize a subharmonic disturbance are derived. When the resonant conditions are met, the periodic wave plays a catalytic role to enhance the growth rate of the subharmonic. The stability characteristics of the subharmonic disturbance, as a function of jet Mach number, jet heating, mode number and the amplitude of the periodic wave, are studied via a secondary instability analysis using two independent but complementary methods: (1) method of multiple scales, and (2) normal mode analysis. It is found that the growth rates of the subharmonic waves with azimuthal numbers beta = 0 and beta = 1 are enhanced strongly, but comparably, when the amplitude of the periodic wave is increased. Furthermore, compressibility at subsonic Mach numbers has a moderate stabilizing influence on the subharmonic instability modes. Heating suppresses moderately the subharmonic growth rate of an axisymmetric mode, and it reduces more significantly the corresponding growth rate for the first spinning mode. Calculations also indicate that while the presence of a finite-amplitude periodic wave enhances the growth rates of subharmonic instability modes, it minimally distorts the mode shapes of the subharmonic waves.

Ng, Lian Lai↗