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 289 records · Page 16

Linear Velocity Problems Part 2: Thermodynamic Considerations

A compressible polytropic gas, in one dimensional planar, cylindrical, or spherical coordinate systems, is examined through four thermodynamic properties: mass density, pressure, specific internal energy (SIE), and entropy, under the assumptions that the fluid velocity is linearly proportional to the radial coordinate and that the Euler gas dynamics equations assert a self similar solution class. The behavior of the thermodynamics is almost entirely determined by an arbitrary function that results from the derivation of the density distribution. Through physical and other arguments, the arbitrary function can be plotted spatially with time profiles, providing a deeper understanding of the system and how the arbitrary function affects the behavior of the compressible gas.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Deformations of Flat Unsymmetric Laminates Subjected to Inplane Loads

The geometrically nonlinear deformation response of initially flat unsymmetric cross-ply laminates subjected to an inplane compressive load and two sets of boundary conditions is studied. Stability of the deformations is considered. At issue is whether or not the plate remains flat with increased compressive loading, and whether it buckles. A semi-infinite unsymmetric cross-ply laminate is used to show the combined effects of geometric nonlinearities and bending-stretch coupling. Finite element results for finite laminates are then presented, and it is shown that to a large degree the boundary conditions control the character of the deformation response. It appears that clamped boundary conditions support buckling behavior, in the classic sense of bifurcation, whereas simply-supported conditions do not.

Hyer, Michael W.↗

Accurate numerical simulations of open quantum systems using spectral tensor trains

Decoherence between qubits is a major bottleneck in quantum computations. Decoherence results from intrinsic quantum and thermal fluctuations as well as noise in the external fields that perform the measurement and preparation processes. With prescribed colored noise spectra for intrinsic and extrinsic noise, we present a numerical method, Quantum Accelerated Stochastic Propagator Evaluation (Q-ASPEN), to solve the time-dependent noise-averaged reduced density matrix in the presence of intrinsic and extrinsic noise. Q-ASPEN is arbitrarily accurate and can be applied to provide estimates for the resources needed to error-correct quantum computations. We employ spectral tensor trains, which combine the advantages of tensor networks and pseudospectral methods, as a variational ansatz to the quantum relaxation problem and optimize the ansatz using methods typically used to train neural networks. Here, the spectral tensor trains in Q-ASPEN make accurate calculations with tens of quantum levels feasible. We present benchmarks for Q-ASPEN on the spin-boson model in the presence of intrinsic noise and on a quantum chain of up to 32 sites in the presence of extrinsic noise. In our benchmark, the memory cost of Q-ASPEN scales as a low-order polynomial in the size of the system once the number of system states surpasses the number of basis functions used in the spectral expansion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Velocity-coupled flow oscillations in a simulated solid-propellant rocket environment

A comprehensive numerical analysis has been carried out to study the unsteady flowfields in a simulated rocket-motor environment. The model is based on the time-dependent compressible Navier-Stokes equations with a two-equation turbulence closure scheme. Various important aspects of the coupling between acoustic oscillations and mean flowfields, including flow reversal, modification of transport properties, etc., are addressed. Results indicate that multi-dimensional effects play important roles in determining local flow structures and wave characteristics. In much of the domain, acoustic velocity nodal points are observed in the near-wall region. The classical one-dimensional theory fails to describe several important mechanisms associated with velocity-induced flow instabilities.

Yang, Vigor↗

A new re-redistribution scheme for weighted state redistribution with adaptive mesh refinement

State redistribution (SRD) is a recently developed technique for stabilizing cut cells that result from finite-volume embedded boundary methods. SRD has been successfully applied to a variety of compressible and incompressible flow problems. When used in conjunction with adaptive mesh refinement (AMR), additional steps are needed to preserve the accuracy and conservation properties of the solution if the embedded boundary is not restricted to a single level of the mesh hierarchy. In this work, we extend the weighted state redistribution algorithm to cases where cut cells live at or near a coarse-fine interface within the domain. Here, we present numerical results that demonstrate that the algorithm is conservative when the coarse-fine interface intersects the embedded boundary. Additionally we compare the numerical solution of the Sod shock tube problem in an inclined cylinder with the analytic solution, and we compare the simulation of a shock hitting a cylindrical obstacle with experimental data. Finally we demonstrate the methodology for simulation of the multicomponent compressible Navier-Stokes equations in a piston-bowl geometry, and discuss the computational efficiency gained by not requiring the entire embedded boundary to be defined at the finest level.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Postbuckling response of long thick plates loaded in compression including higher order transverse shearing effects

This paper presents buckling and postbuckling results for compression-loaded simply-supported aluminum plates and composite plates with a symmetric lay-up of thin +/-45-deg plies composed of many layers. Buckling results for aluminum plates of finite length are given for various length-to-width ratios. Asymptotes to the curves based on the buckling results give N(XCr) for plates of infinite length. Postbuckling results for plates with transverse shearing flexibility are compared to results from classical theory for various width-to-thickness ratios. Characteristic curves indicating the average longitudinal direct stress resultant as a function of the applied displacements are calculated based on four different theories: classical von Karman theory using the Kirchhoff assumptions, first-order shear deformation theory, higher-order shear deformation theory, and three-dimensional flexibility theory. Present results indicate that the three-dimensional flexibility theory gives the lowest buckling loads. The higher-order shear deformation theory has fewer unknowns than the three-dimensional flexibility theory but does not take into account through-the-thickness effects. The figures presented show that small differences occur in the average longitudinal direct stress resultants from the four theories that are functions of applied end-shortening displacement.

Stein, Manuel↗

Variational Asymptotic Homogenization of Finitely Deformed Viscoelastic-Viscoplastic Composites

The objective of this paper is to develop a constitutive model for finitely deformed viscoelastic-viscoplastic materials and a micromechanics approach to homogenizing composites consisting of such materials. The development of the constitutive model involves establishing a thermodynamic framework based on finite strain theory, developing a viscoelasticity and a viscoplasticity model based on the thermodynamic framework, developing a radial return algorithm based on a classic framework, and deriving a closed-from incremental constitutive relation in the spatial configuration. The development of the micromechanics approach involves pulling-back the above constitutive relation to the material configuration, formulating a variational statement with the resulting constitutive relation, discretizing variational statement in a finite-dimensional space, and solving the discretized variational statement using an Euler–Newton predictor–corrector method. The constitutive model is calibrated via monotonic uniaxial compression tests on a polymer, and the calibrated model is validated by comparing its predictions with the cyclic test data. It is shown capable of characterizing viscoelasticity, viscoplasticity, and complex loading paths. The micromechanics approach’s capabilities are demonstrated through homogenizing a unidirectional fiber-reinforced composite, subjected to uniaxial, biaxial, and shear loading, at different strain rates. It is demonstrated to be capable of handling rate dependence and complex loading paths. The present framework can be further improved by implementing more sophisticated viscoelasticity and viscoplasticity models in future work.

Finite element analysis↗

High-temperature superconductivity in LaH 10

The recent discovery of a high critical temperature $T_c$ in compressed H 3 S has been followed by the prediction of Liu et al. of $T_c$ ≈ 250 K in the clathrate LaH 10 structure. This report has been confirmed experimentally by Somayazulu et al. and Drozdov et al. Additional theoretical work by Wang et al. and Quan et al. further established the mechanism of electron-phonon interaction and the dominant role of hydrogen. In the present Rapid Communication we follow the classic McMillan paper, which separates the electron and phonon contributions to the electron-phonon coupling λ. We first compute the numerator of McMillan's expression, the Hopfield parameter η, using the theory of Gaspari and Gyorffy (GG), and obtain the force constants in the denominator from Wang et al. and Quan et al. The resulting λ is used in the Allen-Dynes equation to calculate $T_c$. The value of $T_c$ reaches a maximum in the range of 236–263 K at pressures of 255 GPa and decreases for smaller or larger pressures. Further, we provide a thorough analysis of the different terms of the GG equation and draw the conclusion that the $\textit{sp}$ channel of the hydrogen is the most important contribution to obtain high values of $T_c$ in this material. Consistent with Wang et al., we find large values of λ that decrease with increasing pressure. In addition, we find that the hydrogen sites are the largest contributors to the total value of the coupling constant λ. That is, the acoustic mode associated with La contributes only 2% to the total λ, while the optic modes associated with H contribute 18% for the H1 site and 80% for the H2 site. These relative contributions to λ are consistent with those given by Wang et al. and by Quan et al. Thus, our results strongly support the view that LaH 10 is a metallic hydrogen superconductor.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Postbuckling response of long thick plates loaded in compression including higher order transverse shearing effects

Buckling and postbuckling results are presented for compression-loaded simply-supported aluminum plates and composite plates with a symmetric lay-up of thin + or - 45 deg plies composed of many layers. Buckling results for aluminum plates of finite length are given for various length-to-width ratios. Asymptotes to the curves based on buckling results give N(sub xcr) for plates of infinite length. Postbuckling results for plates with transverse shearing flexibility are compared to results from classical theory for various width-to-thickness ratios. Characteristic curves indicating the average longitudinal direct stress resultant as a function of the applied displacements are calculated based on four different theories: Classical von Karman theory using the Kirchoff assumptions, first-order shear deformation theory, higher-order shear deformation theory, and 3-D flexibility theory. Present results indicate that the 3-D flexibility theory gives the lowest buckling loads. The higher-order shear deformation theory has fewer unknowns than the 3-D flexibility theory but does not take into account through-the-thickness effects. The figures presented show that small differences occur in the average longitudinal direct stress resultants from the four theories that are functions of applied end-shortening displacement.

Stein, Manuel↗

Postbuckling analysis of shear deformable composite flat panels taking into account geometrical imperfections

The effects of initial geometrical imperfections on the postbuckling response of flat laminated composite panels to uniaxial and biaxial compressive loading are investigated analytically. The derivation of the mathematical model on the basis of first-order transverse shear deformation theory is outlined, and numerical results for perfect and imperfect, single-layer and three-layer square plates with free-free, clamped-clamped, or free-clamped edges are presented in graphs and briefly characterized. The present approach is shown to be more accurate than analyses based on the classical Kirchhoff plate model.

Librescu, L.↗

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 that include 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. Using several simplifying assumptions, linearized, stability equations are derived that describe the response of the plate just after bifurcation buckling occurs. These equations are then modified to allow the plate reference surface to be located a distance z(c), from the centroid surface which is convenient for modeling stiffened-plate assemblies. The implementation of the new theory into the VICONOPT buckling and vibration analysis and optimum design program code is described. Either classical plate theory (CPT) or first-order shear-deformation plate theory (SDPT) may be selected in VICONOPT. Comparisons of numerical results for several example problems with different loading states are made. Results from the new curved-plate analysis compare well with closed-form solution results and with results from known example problems in the literature. Finally, a design-optimization study of two different cylindrical shells subject to uniform axial compression is presented.

McGowan, David M.↗

Some Aeronautical Communications Experiments

Classically there has existed an asymmetry between the computing and communicating sides of aerospace systems. Over the past few decades, this asymmetry has shifted to favoring communication link technologies, meaning that advancements in available central processing units (CPUs), storage devices, and internal data buses have stagnated. Indeed, the increased emphasis placed on refining subsystem performance such as with antenna bandwidth in phased arrays, amplifier power efficiency, software defined radio (SDR) flexibility and encoding for data compression and error correction has given rise to successful debuts of multi-gigabit-per-second data return from long space-link distances. These accomplishments are easily quantifiable through link budgets and illustrate what is possible, but also reveal the deficiencies in overall communications capabilities. In particular, the ever-accelerating presence of aerospace vehicles gives rise to newer and larger classes of challenges to address the needs of 21st century systems. Furthermore remote sensing and imaging capabilities have far outpaced our ability to transmit their products to the ground, so we are increasingly dependent on pre-processing and downselection to contend with the communications bottleneck. No longer may we depend upon the constrained logistics in delivering end-to-end data delivery through manual reconfigurations, static event scheduling and execution on a per-vehicle basis, for these methods do not scale and therefore must give way to dynamic, networked approaches with an overall systems view in mind. Emerging mission requirements exhibit a trend toward multiple smaller-scale vehicles working together to perform dissimilar observations. Such operations necessitate sensor fusion across a constellation, and where data processing may be distributed throughout a fairly disconnected network whose topology changes over time in non-deterministic manners. Individual communications link performance is still very relevant to deploying an effective communications system, but now must be embedded within a greater architecture of capability to optimally utilize the bandwidth available from each link to generate an ultimate end-to-end quality of service. The deleterious effects of timing uncertainty across the arrangement presents a challenge to measurement synchronization and delivery, so a successful deployed system needs to be tolerant to the delays inherent in time-of-light between elements and digital processing latencies existing at each node. In this presentation we share the flight test results from a high performance Gbps laser communications terminal evaluated with a suite of store and forward capabilities called High-rate Delay Tolerant Networking (HDTN). The communications payload is operated over Lake Erie across a range of configurations including several convergence layers, and is evaluated to determine recovery time after link disruptions, information loss, efficiency and speed. The effectiveness of utilizing a flying laboratory to increase the Technology Readiness Level (TRL) of an integrated system in relevant environments is discussed, as well as the value of conducting aeronautics experiments to retire risk for technology infusion into space missions. Upcoming flight campaigns will be presented, including opportunities to demonstrate secure command and control, data intensive hyperspectral imaging, quantum link characterization, 4k High Definition (HD) video streaming and internetworked space-ground-aero relay operations. These experiments will pave the way for future missions which will depend upon interoperability across disparate government and privately owned networks, involve contention with uncertain and dynamic timing, and require agility to autonomously configure optimal parameters across networks of ever-increasing size and complexity to ensure data delivery. https://www1.grc.nasa.gov/space/scan/acs/tech-studies/dtn/

Daniel Raible↗

Thermophysical Properties of Liquid Tritium: A Path Integral Monte Carlo Study

Here, we present worm-algorithm, path integral Monte Carlo simulations of bulk liquid tritium. The simulations are benchmarked against empirically known thermophysical properties of liquid deuterium and liquid tritium. Results for the pair correlation function, chemical potential, isothermal compressibility, isochoric heat capacity, and single-particle momentum distributions are reported. Given the benchmark comparisons, our predictions of liquid tritium properties are expected to be accurate to within a few percent. Our simulations unambiguously demonstrate the significance of nuclear quantum effects to the properties of liquid tritium. In particular, under saturated vapor pressure, the average molecular kinetic energy of the liquid is found to be more than 60% higher than the value expected from the classical equipartition theorem.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Colloids and Nucleation

The objectives of the work funded under this grant were to develop a microphotographic technique and use it to monitor the nucleation and growth of crystals of hard colloidal spheres. Special attention is given to the possible need for microgravity studies in future experiments. A number of persons have been involved in this work. A masters student, Keith Davis, began the project and developed a sheet illumination apparatus and an image processing system for detection and analysis. His work on a segmentation program for image processing was sufficient for his master's research and has been published. A post doctoral student Bernie Olivier and a graduate student Yueming He, who originally suggested the sheet illumination, were funded by another source but along with Keith made photographic series of several samples (that had been made by Keith Davis). Data extraction has been done by Keith, Bernie, Yueming and two undergraduates employed on the grant. Results are published in Langmuir. These results describe the sheet lighting technique as one which illuminates not only the Bragg scattering crystal, but all the crystals. Thus, accurate crystal counts can be made for nucleation rate measurements. The strange crystal length scale reduction, observed in small angle light scattering (SALS) studies, following the initial nucleation and growth period, has been observed directly. The Bragg scattering (and dark) crystal size decreases in the crossover region. This could be an effect due to gravitational forces or due to over- compression of the crystal during growth. Direct observations indicate a complex morphology for the resulting hard sphere crystals. The crystal edges are fairly sharp but the crystals have a large degree of internal structure. This structure is a result of (unstable) growth and not aggregation. As yet unpublished work compares growth exponents data with data obtained by SALS. The nucleation rate density is determined over a broad volume fraction range (0.505 to 0.575) and compared to classical nucleation rate calculations. The experimental nucleation rate density varies with volume fraction more slowly than expected theoretically. Gravity may influence nucleation because crystal nuclei are more dense than the surrounding fluid and may be gravitationally convected. The largest observed crystal size was near volume fraction 0.52. Near the freezing point, crystals should grow to large size because nuclei are sparse and there is plenty of room to grow without contacting other crystals. Near melting the opposite is true, with crystals growing into one another shortly after nucleation. The maximum crystal size seen near 0.52 volume fraction suggests that gravitational stresses are limiting crystal growth near freezing.

Ackerson, Bruce↗

Experimental quantum learning of a spectral decomposition

Currently available quantum hardware allows for small-scale implementations of quantum machine learning algorithms. Such experiments aid the search for applications of quantum computers by benchmarking the near-term feasibility of candidate algorithms. Here we demonstrate the quantum learning of a two-qubit unitary by a sequence of three parameterized quantum circuits containing a total of 21 variational parameters. Moreover, we variationally diagonalize the unitary to learn its spectral decomposition, i.e., its eigenvalues and eigenvectors. We illustrate how this can be used as a subroutine to compress the depth of dynamical quantum simulations. One can view our implementation as a demonstration of entanglement-enhanced machine learning, as only a single (entangled) training data pair is required to learn a 4 × 4 unitary matrix.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Viscous aerodynamic design using the adjoint variable approach

The use of classical optimal control methods, in particular variational methods, to solve the airfoil optimization problem, by deriving a set of adjoint (costate) equations and boundary conditions has already been done for inviscid (potential and Euler flows) and two dimensional, steady state, incompressible flow governed by the Navier-Stokes equations. The interior and boundary terms of the volume integral have been derived (in this work) for the steady Navier-Stokes equations in three dimensions for a viscous, compressible heat conducting fluid. This can be used to derive the adjoint equations and numerical boundary conditions for general classes of problems and hence paves the way for a solution to the aerodynamic optimization problem for compressible viscous flows. The next steps to the realization of that goal are projected as below. The usual square integral pressure functional as an objective function is being replaced by a more realistic drag functional subject to a lift constraint. The feasibility of attempting the more difficult time dependent problem is being investigated. It remains to get the full system of adjoint equations and boundary conditions with the new functional. The state and adjoint equations must be discretized and coded. An appropriate optimization program must be used (steepest descents seems inadequate) and various known airfoil shapes should be recovered in test cases of the computer program.

DeRise, George↗

Quantum-Inspired Bayesian Sampling for Uncertainty Quantification and Machine Learning (Final Technical Report)

With increasing simulation and measurement data, machine learning and artificial intelligence have been widely used in computational decision-making of complex engineering systems. The resulting tools, such as uncertainty quantification solvers, reinforcement learning, and physics-informed machine learning, have achieved great success in critical DOE tasks such as material discovery and design, energy system modeling and control, and numerical weather and climate prediction. A core topic in scientific machine learning and artificial intelligence is Bayesian inference: given an observed data set, people want to estimate the posterior distribution of a (possibly large) number of hidden parameters. Due to the flexibility and weak assumptions, Bayesian sampling has been the mainstream Bayesian inference solvers despite the rapid progress of approximate Bayesian inference. Classical Bayesian sampling methods such as Markov-chain Monte Carlo suffer from a low-acceptance rate due to the random walk nature, therefore state-of-the-art techniques use Hamiltonian Monte Carlo and its variants to efficiently draw posterior samples in a high dimension. The key idea of Hamiltonian Monte Carlo and its variants is to simulate the Hamiltonian dynamics of a classical particle with a fixed mass, and their performance significantly degrades when the posterior distribution is highly spiky or has multiple modes. Leveraging the idea of quantum physics, this project has investigated new theory, algorithms and applications of Bayesian inference (especially Bayesian sampling). The main results include: (1) novel quantum-inspired Bayesian sampling methods that can lead to better accuracy for challenging multi-modal or spiky distributions, (2) more scalable machine learning framework leveraging tensor-compressed Bayesian inference, and (3) Bayesian and sampling approaches for verifying the robustness of continuous and binary neural networks.

97 MATHEMATICS AND COMPUTING↗

Building a new multiphysics workflow in MOOSE: application to tritium migration, trapping and advection in TMAP8

Fusion devices are anticipated to produce and consume several kilograms of tritium per year. This rare fuel resource is both highly mobile and radioactive, making tracking inventories a priority for operation and safety. The fusion safety program at the Idaho National Laboratory has been developing the Tritium Migration and Analysis Program (TMAP), of which the latest version is a MOOSE-based application. TMAP8 is verified against its predecessors and possesses additional multi-dimensional tritium migration modeling capabilities. As we extend its capabilities towards both whole device (in multiple dimensions) and whole plant (with multiple components) simulations, the syntax of inputs must become compact, descriptive, compatible with quality assurance processes, and as error-proof as achievable. The new Physics system developed MOOSE can set up equations and instantiating them on plant components. The system permits the automatic definition of complex discretization with a consistency between object parameters achieved programmatically. The Physics system can currently instantiate the equations for heat conduction and Navier Stokes weakly compressible flow. In MOOSE-terms, it automates the definition of kernels, boundary conditions, and several core and helper materials and fields. As part of this effort, Physics classes were developed for tritium migration, trapping and advection within either a multi-dimensional Navier Stokes fluid dynamics simulation, or a 1D thermal hydraulics piping system. In this presentation, we will showcase the new syntax, its application to several verification and validation cases which were already studied using the classical TMAP8 syntax, and a demonstration of the new coupling capabilities for the migration of tritium into blanket coolant channels and the subsequent advection into the coolant loop.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗