Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuum mechanics”

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 487 records · Page 27

Vibrational Stabilization of Flexible Structures

It has been demonstrated that a high frequency excitation (HFE) field significantly changes the fundamental properties of mechanical systems. The most important contribution of HFE is the stiffening effect of an elastic continuum in the direction of the wave vector. This effect allows control of stiffness in any selected direction by the corresponding changes in the intensity of HFE. This new approach can be effective for large flexible space structures. Such an approach may prove to be very practical in the sense that large structures need to be made as flimsy as possible for low cost under ordinary situations. However, for certain operations such as development, orbital transfer, docking, and other circumstances, it would be vital to have a means of temporarily stiffening certain structural members. If the structure was designed to meet these occasional loads without temporary stiffening, it would be considerably more massive and more expensive.

Zak, M.↗

Comparison of solar hard X-ray and UV line and continuum bursts with high time resolution

A comparison of data sets from the UV Spectrometer and Polarimeter and Hard X-ray Burst Spectrometer instruments on SMM has established the close relationship of the impulsive phase hard X-ray and UV continuum and OV line emissions, lending support to the notion that they have a similar origin low in the solar atmosphere. These results severely constrain models that attempt to explain impulsive phase hard X-rays and UV emission; alternative processes of impulsive-phase UV continuum production should accordingly be considered. Attention is given to an electron beam 'hole boring' mechanism and a photoionization radiation transport mechanism.

Orwig, L. E.↗

Continuum radiation in planetary magnetospheres

With the completion of the Voyager tour of the outer planets, radio and plasma wave instruments have executed the first survey of the wave spectra of Earth, Jupiter, Saturn, Uranus, and Neptune. One of the most notable conclusions of this survey is that there is a great deal of qualitative similarity in both the plasma wave and radio wave spectra from one magnetosphere to the next. In particular, in spite of detailed differences, most of the radio emissions at each of the planets have been tentatively classified into two primary categories. First, the most intense emissions are generally associated with the cyclotron maser instability. Second, a class of weaker emissions can be found at each of the magnetospheres which appears to be the result of conversion from intense electrostatic emissions at the upper hybrid resonance frequency into (primarily) ordinary mode radio emission. It is this second category, often referred to as nonthermal continuum radiation, which we will discuss in this review. We review the characteristics of the continuum spectrum at each of the planets, discuss the source region and direct observations of the generation of the emissions where available, and briefly describe the theories for the generation of the emissions. Over the past few years evidence has increased that the linear mode conversion of electrostatic waves into the ordinary mode can account for at least some of the continuum radiation observed. There is no definitive evidence which precludes the possibility that a nonlinear mechanism may also be important.

Kurth, W. S.↗

A Viscoplastic Constitutive Theory for Monolithic Ceramic Materials

With increasing use of ceramic materials in high temperature structural applications such as advanced heat engine components, the need arises to accurately predict thermomechanical behavior. This paper, which is the first of two in a series, will focus on inelastic deformation behavior associated with these service conditions by providing an overview of a viscoplastic constitutive model that accounts for time-dependent hereditary material deformation (e.g., creep, stress relaxation, etc.) in monolithic structural ceramics. Early work in the field of metal plasticity indicated that inelastic deformations are essentially unaffected by hydrostatic stress. This is not the case, however, for ceramic-based material systems, unless the ceramic is fully dense. The theory presented here allows for fully dense material behavior as a limiting case. In addition, ceramic materials exhibit different time-dependent behavior in tension and compression. Thus, inelastic deformation models for ceramics must be constructed in a fashion that admits both sensitivity to hydrostatic stress and differing behavior in tension and compression. A number of constitutive theories for materials that exhibit sensitivity to the hydrostatic component of stress have been proposed that characterize deformation using time-independent classical plasticity as a foundation. However, none of these theories allow different behavior in tension and compression. In addition, these theories are somewhat lacking in that they are unable to capture creep, relaxation, and rate-sensitive phenomena exhibited by ceramic materials at high temperature. When subjected to elevated service temperatures, ceramic materials exhibit complex thermomechanical behavior that is inherently time-dependent, and hereditary in the sense that current behavior depends not only on current conditions, but also on thermo-mechanical history. The objective of this work is to present the formulation of a macroscopic continuum theory that captures these time-dependent phenomena. Specifically, the overview contained in this paper focuses on the multiaxial derivation of the constitutive model, and examines the scalar threshold function and its attending geometrical implications.

Janosik, Lesley A.↗

Updates on the Predictive Materials Modeling Software Tools

Updates on NASA‘s efforts to build a Predictive Material Modeling (PMM) framework from the micro-scale to the macro-scale are presented in this abstract. The PMM effort is part of the Entry Systems Modeling (ESM) project under NASA’s Game Changing Development (GCD) program. To reduce the need for extensive testing and accelerate the design cycle process, ESM is developing simulation and modeling tools that enable the characterization of the properties of thermal protection materials and their response to extremely hot plasma. The Porous Microstructure Analysis (PuMA) software has been developed to compute effective material properties and perform material response simulations on digitized microstructures of porous media. PuMA is able to import three-dimensional digital images obtained from X-ray microtomography or to generate artificial microstructures that mimic real materials. PuMA also provides a module for interactive 3D visualizations. Version 3, which was recently released as open-source, includes modules to compute simple morphological properties such as porosity, volume fractions, pore diameter, and specific surface area. Additional capabilities include the determination of effective thermal and electrical conductivity (both radiative and solid conduction - including the ability to simulate local anisotropy for the latter); effective diffusivity and tortuosity from the continuum to the rarefied regime; techniques to determine the local material orientation, as well as mechanical properties (elasticity coefficients), and permeability. Computed properties are then used to inform a macro-scale material response model, such as those implemented in the Porous material Analysis Toolbox based on OpenFOAM (PATO) software developed within ESM. The computational model in PATO is a generic heat and mass transfer model for porous reactive materials containing several solid phases and a single gas phase. The detailed chemical interactions occurring between the solid phases and the gas phase are modeled at the pore scale, assuming Local Thermal Equilibrium. Recent efforts include the development of a mechanical erosion model as well as a unified model allowing an intrinsic coupling between fluid and material. Comparison to flight data (Mars Science Laboratory [MSL] Entry Descent and Landing Instrument [MEDLI] and Mars 2020 MEDLI2) is critical in order to validate these computational tools. Examples of ablative material response using the code will be presented, including 3D simulations of the full-scale heatshield of the MSL capsule. The simulations demonstrated the ability of the modern material response code, PATO, to handle the material response of geometrically complex and large domains through the use of massively parallel computations.

material modeling↗

Coarse-grained methods for heterogeneous vesicles with phase-separated domains: Elastic mechanics of shape fluctuations, plate compression, and channel insertion

We develop coarse-grained particle approaches for studying the elastic mechanics of vesicles with heterogeneous membranes having phase-separated domains. We perform simulations both of passive shape fluctuations and of active systems where vesicles are subjected to compression between two plates or subjected to insertion into narrow channels. Analysis methods are developed for mapping particle configurations to continuum fields with spherical harmonics representations. Heterogeneous vesicles are found to exhibit rich behaviors where the heterogeneity can amplify surface two-point correlations, reduce resistance during compression, and augment vesicle transport times in channels. The developed methods provide general approaches for characterizing the mechanics of coarse-grained heterogeneous systems taking into account the roles of thermal fluctuations, geometry, and phase separation.

97 MATHEMATICS AND COMPUTING↗

Determinants of the mechanical properties of bones

The mechanical properties of bones are governed by the same principles as those of man-made load-bearing structures, but the organism is able to adapt its bone structure to changes in skeletal loading. In this overview of the determinants of the strength and stiffness of bone, a continuum approach has been taken, in which the behavior of a macroscopic structure depends on its shape and size, and on the mechanical properties of the material within. The latter are assumed to depend on the composition (porosity and mineralization) and organization (trabecular or cortical bone architecture, collagen fiber orientation, fatigue damage) of the bone. The effects of each of these factors are reviewed. Also, the possible means of non-invasively estimating the strength or other mechanical properties of a bone are reviewed, including quantitative computed tomography, photon absorptiometry, and ultrasonic measurements. The best estimates of strength have been obtained with photon absorptiometry and computed tomography, which at best are capable of accounting for 90% of the strength variability in a simple in vitro test, but results from different laboratories have been highly variable.

Non-NASA Center↗

Complications in the interpretation of the charge-asymmetry-dependent π flow for the chiral magnetic wave

The charge asymmetry ($A_ch$) dependence of the $π^₋$ and $π^+$ elliptic flow difference, $Δv_2(A_{\text{ch}})$, has been regarded as a sensitive observable for the possible chiral magnetic wave (CMW) in relativistic heavy ion collisions. In this work, we first demonstrate that, due to nonflow backgrounds, the flow measurements by the Q-cumulant method using all charged particles as reference introduce a trivial linear term to $Δv_2(A_{\text{ch}})$. The trivial contribution to the slope of triangle flow difference $Δv_3( A_{\text{ch}})$ can be negative if the nonflow is dominated by back-to-back pairs. After eliminating the trivial term, we find that the nonflow between like-sign pairs gives rise to an additional positive slope to $Δv_2(A_{\text{ch}})$ because of the larger multiplicity dilution to $π^+$ ($π^₋$) at positive (negative) $A_{\text{ch}}$. We further find that the competition between different π sources can introduce another nontrivial linear-$A_{\text{ch}}$ term due to their different multiplicity fluctuations and anisotropic flows. We then study the effect of neutral cluster (resonance) decays as a mechanism for local charge conservation on the slope parameter of $Δv_2(A_{\text{ch}})$. We find that the slope parameter is sensitive to the kinematics of those neutral clusters. Light resonances give positive slopes while heavy resonances give negative slopes. Local charge conservation from continuum cluster mass distribution can give a positive slope parameter comparable to experimental data. Our studies indicate that many non-CMW physics mechanisms can give rise to a $A_{\text{ch}}$-dependent $Δv_2(A_{\text{ch}})$ and the interpretation of $Δv_2(A_{\text{ch}})$ in terms of the CMW is delicate.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Modeling the formation of Sedan Crater using the FLAG and HOSS codes

Numerical modeling of explosion crater formation requires accounting for complex physical processes. Numerical validation of explosion cratering is an important step in modeling and requires experimental data for comparison. Models using discrete elements and continuum models have both benefits and drawbacks to their approaches. In this work, we consider both an arbitrary Lagrangian–Eulerian (ALE) hydrocode and a finite discrete element method (FDEM) approach to modeling the formation of the Sedan crater, the largest human-made crater in the United States. The Sedan crater formed from an underground nuclear detonation in the Nevada desert as part of Project Plowshare. Our models show that the continuum approach of the hydrocode matched well compared to early test time prior to the mound rupture and subsequent fireball venting, when most of the alluvium exhibited fluid behavior. Our FDEM approach matched the final crater dimensions well, after material had settled back into the crater, when material strength and solid mechanics play key roles. Our work shows how leveraging the benefits of multiple numerical approaches can lead to better understanding of complex physical problems, especially problems with limited experimental data. By using a continuum approach to early-time hydrodynamics and an FDEM approach to later-time solid mechanics, we can better understand the different physical regimes of explosion crater formation.

36 MATERIALS SCIENCE↗

Micro-Mechanical Modeling of Ductile Fracture in Welded Aluminum-Lithium Alloys

This study describes a 3-D computational framework to model stable extension of a macroscopic crack under mode I conditions in ductile metals. Material separation occurs through a local fracture mechanism based on the growth and coalescence of micro voids dispersed in the material. The Gurson-Tvergaard plasticity model for void materials describes the eventual loss of material stress capacity under sustained loading. Fixed-sized computational cell elements defined on a thin layer adjacent to the crack plane provide an implicit length scale for continuum damage model. Outside this layer, the material remains undamaged by void growth, consistent with metallurgical observations. The key micro-mechanics parameters are D, the thickness of cell layer, and f 0 , the initial cell porosity.

Ahmed Ibrahim↗

Multiscale Modeling of Ultra High Temperature Ceramics (UHTC) ZrB2 and HfB2: Application to Lattice Thermal Conductivity

We are developing a multiscale framework in computational modeling for the ultra high temperature ceramics (UHTC) ZrB2 and HfB2. These materials are characterized by high melting point, good strength, and reasonable oxidation resistance. They are candidate materials for a number of applications in extreme environments including sharp leading edges of hypersonic aircraft. In particular, we used a combination of ab initio methods, atomistic simulations and continuum computations to obtain insights into fundamental properties of these materials. Ab initio methods were used to compute basic structural, mechanical and thermal properties. From these results, a database was constructed to fit a Tersoff style interatomic potential suitable for atomistic simulations. These potentials were used to evaluate the lattice thermal conductivity of single crystals and the thermal resistance of simple grain boundaries. Finite element method (FEM) computations using atomistic results as inputs were performed with meshes constructed on SEM images thereby modeling the realistic microstructure. These continuum computations showed the reduction in thermal conductivity due to the grain boundary network.

Lawson, John W.↗

Multiscale Modeling of UHTC: Thermal Conductivity

We are developing a multiscale framework in computational modeling for the ultra high temperature ceramics (UHTC) ZrB2 and HfB2. These materials are characterized by high melting point, good strength, and reasonable oxidation resistance. They are candidate materials for a number of applications in extreme environments including sharp leading edges of hypersonic aircraft. In particular, we used a combination of ab initio methods, atomistic simulations and continuum computations to obtain insights into fundamental properties of these materials. Ab initio methods were used to compute basic structural, mechanical and thermal properties. From these results, a database was constructed to fit a Tersoff style interatomic potential suitable for atomistic simulations. These potentials were used to evaluate the lattice thermal conductivity of single crystals and the thermal resistance of simple grain boundaries. Finite element method (FEM) computations using atomistic results as inputs were performed with meshes constructed on SEM images thereby modeling the realistic microstructure. These continuum computations showed the reduction in thermal conductivity due to the grain boundary network.

Lawson, John W.↗

Measurements of the structure of an ionizing shock wave in a hydrogen-helium mixture

Shock structure during ionization of a hydrogen-helium mixture was studied using hydrogen line and continuum emission measurements. A reaction scheme is proposed which includes hydrogen dissociation and a two-step excitation-ionization mechanism for hydrogen ionization by atom-atom and atom-electron collisions. Agreement was achieved between numerical calculations and measurements of emission intensity as a function of time for shock velocities from 13 to 20 km/sec in a 0.208 H2 - 0.792 He mixture. The electron temperature was found to be significantly different from the heavy particle temperature during much of the ionization process. Similar time histories for H beta and continuum emission indicate upper level populations of hydrogen in equilibrium with the electron concentration during the relaxation process.

Leibowitz, L. P.↗

A data-driven multiscale model for reactive wetting simulations

Here, we describe a data-driven, multiscale technique to model reactive wetting of a silver–aluminum alloy on a Kovar™ (Fe-Ni-Co alloy) surface. We employ molecular dynamics simulations to elucidate the dependence of surface tension and wetting angle on the drop’s composition and temperature. A design of computational experiments is used to efficiently generate training data of surface tension and wetting angle from a limited number of molecular dynamics simulations. The simulation results are used to parameterize models of the material’s wetting properties and compute the uncertainty in the models due to limited data. The data-driven models are incorporated into an engineering-scale (continuum) model of a silver–aluminum sessile drop on a Kovar™ substrate. Model predictions of the wetting angle are compared with experiments of pure silver spreading on Kovar™ to quantify the model-form errors introduced by the limited training data versus the simplifications inherent in the molecular dynamics simulations. The paper presents innovations in the determination of “convergence” of noisy MD simulations before they are used to extract the wetting angle and surface tension, and the construction of their models which approximate physio-chemical processes that are left unresolved by the engineering-scale model. Together, these constitute a multiscale approach that integrates molecular-scale information into continuum scale models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Measurements of the structure of an ionizing shock wave in a hydrogen-helium mixture.

Shock structure during ionization of a hydrogen-helium mixture has been followed using hydrogen line and continuum emission measurements. A reaction scheme is proposed which includes hydrogen dissociation and a two-step excitation-ionization mechanism for hydrogen ionization by atom-atom and atom-electron collisions. Agreement has been achieved between numerical calculations and measurements of emission intensity as a function of time for shock velocities from 13 to 20 km/sec in a 0.208 H2-0.792 He mixture. The electron temperature was found to be significantly different from the heavy particle temperature during much of the ionization process. Similar time histories for H beta and continuum emission indicate upper level populations of hydrogen in equilibrium with the electron concentration during the relaxation process.

Leibowitz, L. P.↗

PowderSim: Lagrangian Discrete and Mesh-Free Continuum Simulation Code for Cohesive Soils

PowderSim is a calculation tool that combines a discrete-element method (DEM) module, including calibrated interparticle-interaction relationships, with a mesh-free, continuum, SPH (smoothed-particle hydrodynamics) based module that utilizes enhanced, calibrated, constitutive models capable of mimicking both large deformations and the flow behavior of regolith simulants and lunar regolith under conditions anticipated during in situ resource utilization (ISRU) operations. The major innovation introduced in PowderSim is to use a mesh-free method (SPH-based) with a calibrated and slightly modified critical-state soil mechanics constitutive model to extend the ability of the simulation tool to also address full-scale engineering systems in the continuum sense. The PowderSim software maintains the ability to address particle-scale problems, like size segregation, in selected regions with a traditional DEM module, which has improved contact physics and electrostatic interaction models.

Johnson, Scott↗

Γ valley transition metal dichalcogenide moiré bands

Significance At small twist angles, moiré patterns form in transition metal dichalcogenide bilayers. Using sophisticated ab initio calculations, we develop a model able to describe the band structure of these systems. We found that the moiré pattern acts as a long-wavelength modulating potential, with a sixfold rotational symmetry not present in the underlying lattice. The modulation potential reduces hole band widths and localizes electrons or holes around honeycomb or kagome lattice sites.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Investigation of mode activity in NBI-heated experiments of Wendelstein 7-X

The 2018 operation phase (OP 1.2b) of the stellarator Wendelstein 7-X (W7-X) included, for the first time, neutral beam injection (NBI) to heat the plasma. Since the injection geometry at W7-X is not parallel, this generates both passing and trapped fast particles. During longer phases of NBI injection, with the primary purpose to study the heating efficiency of this system, Alfvén eigenmodes (AEs) were observed by a number of diagnostics, including the phase contrast imaging (PCI) system, the magnetic pick-up coils (Mirnov coils), and the soft x-ray multi-camera tomography system (XMCTS). Alfvén eigenmodes are of great interest for future fusion reactors as it has been shown that the resonant interaction of fast ions with self-excited AEs can lead to enhanced transport of fast ions and potentially to energy losses. This is especially true for so-called gap-modes, Alfvén eigenmodes with frequencies in gaps of the continuous spectrum, since they lack continuum damping. These modes are commonly known to be excited by fast ions, but other destabilizing mechanisms, e.g. the electron-pressure gradient are also possible. In this article we present a first analysis of the experimentally observed frequencies from the theoretical side. The calculation of shear Alfvén wave continua for selected cases and the assignment of observed frequencies to the gaps of the continuous spectra are presented. Using the ideal-MHD code CKA (Könies A. 2007 10th IAEA TM on Energetic Particles in Magnetic Confinement System), we find gap modes that match the experimental measurements in terms of the observed frequencies. We emphasize the crucial roles played by the coupling of sound and Alfvén waves as well as of the Doppler shift arising as a consequence of the radial electric field in W7-X. We employ the perturbative gyrokinetic code CKA-EUTERPE (Feh´er 2013 Simulation of the interaction between Alfv´en waves and fast particles), using a slowing-down distribution function for the fast ions as calculated by the Monte-Carlo particle following code ASCOT (Hirvijoki et al 2014 Comput. Phys. Commun. 185 1310–21) to assess the fast-ion drive. Here we find that the fast-ion drive is insufficient to overcome the background-plasma damping. The fact that unstable modes were observed experimentally may point to problems with the modelling or indicate the existence of other destabilizing mechanisms, e.g. associated with the electron-pressure gradient (Windisch et al 2017 Plasma Phys. Control. Fusion 59 105002) that sensitively depend on the profiles of the background plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗