Search NASA⌕ Search

SEARCH · Search NASA

Results for “continuum mechanics of solids”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

The effects of strain heating in lithospheric stretching models

The deformation by stretching of a continental type lithosphere has been formulated so that the problem can be solved by a continuum mechanical approach. The deformation, stress state, and temperature distribution are constrained to satisfy the physical laws of conservation of mass, energy, momentum, and an experimentally defined rheological response. The conservation of energy equation including a term of strain energy dissipation is given. The continental lithosphere is assumed to have the rheology of an isotropic, incompressible, nonlinear viscous, two layered solid.

Stanton, M.↗
DOE OSTI2020

Emergence of p+ip superconductivity in two-dimensional doped Dirac systems

Searching for the $p + ip$ superconducting (SC) state has become a fascinating subject in condensed matter physics recently, as a dream application awaiting in topological quantum computation. Unfortunately, so far there is no universal principle for realizing $p + ip$ in generic solid-state systems. Here we report a theoretical discovery of a $p + ip$ SC ground state (coexisting with ferromagnetic order) in the honeycomb lattice Hubbard model in the extremely strong-coupling limit (e.g., infinite U) at low doping $(δ < 0.2)$, by using both the state-of-art Grassmann tensor product state approach and a continuum quantum field theory approach. Our discovery suggests a mechanism for the $p + ip$ SC state in generic strongly correlated systems based on spin-charge separation and the charge current–current coupling scenario, which opens a door towards experimental realization. The $p + ip$ SC state has an instability towards a potential non-Fermi liquid with a large but finite U . Nevertheless, by applying an in-plane Zeeman field, such a $p + ip$ SC state can be stabilized with finite U in a very wide range of doping. Relevant realistic materials are also proposed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Young but fading radio sources: searching for remnants among compact steep-spectrum radio sources

The incidence of young but fading radio sources provides important information on the life cycle of radio emission in radio-loud active galactic nuclei. Despite its importance for constraining the models of radio source evolution, there are no systematic studies of remnants in complete samples of young radio sources. We report results of the study of 18 compact steep-spectrum (CSS) radio sources, selected from the statistically complete B3-VLA (Very Large Array) CSS sample, characterized by a steep optically thin spectrum (α ≥ 1.0) and no core detection in earlier studies. Our deep multifrequency VLA, pc-scale Very Long Baseline Array, and enhanced Multi Element Remotely Linked Interferometer Network (e-MERLIN) observations allowed us to locate the core component in 10 objects. In three CSS sources, there is no clear evidence of present-time active regions, suggesting that they are likely in a remnant phase. Among sources with core detection, we find three objects that have no clear active regions (hotspots) at the edges of the radio structure, suggesting that the radio emission may have just restarted. Our results support a power-law distribution of the source ages, although the poor statistics prevents us from setting solid constraints on the percentage of remnants and restarted sources in subpopulations of radio sources.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

A data-driven peridynamic continuum model for upscaling molecular dynamics

Nonlocal models, including peridynamics, often use integral operators that embed lengthscales in their definition. However, the integrands in these operators are difficult to define from the data that are typically available for a given physical system, such as laboratory mechanical property tests. In contrast, molecular dynamics (MD) does not require these integrands, but it suffers from computational limitations in the length and time scales it can address. To combine the strengths of both methods and to obtain a coarse-grained, homogenized continuum model that efficiently and accurately captures materials’ behavior, we propose a learning framework to extract, from MD data, an optimal Linear Peridynamic Solid (LPS) model as a surrogate for MD displacements. To maximize the accuracy of the learnt model we allow the peridynamic influence function to be partially negative, while preserving the well-posedness of the resulting model. To achieve this, we provide sufficient well-posedness conditions for discretized LPS models with sign-changing influence functions and develop a constrained optimization algorithm that minimizes the equation residual while enforcing such solvability conditions. This framework guarantees that the resulting model is mathematically well-posed, physically consistent, and that it generalizes well to settings that are different from the ones used during training. We illustrate the efficacy of the proposed approach with several numerical tests for single layer graphene. Our two-dimensional tests show the robustness of the proposed algorithm on validation data sets that include thermal noise, different domain shapes and external loadings, and discretizations substantially different from the ones used for training.

homogenization↗
DOE OSTI2023

Fierro

FIERRO is a parallel C++ code designed to simulate fluid mechanics, heat transfer, and solid mechanics in two- and three dimensional space. FIERRO is written to run on homogeneous (CPU) and heterogeneous (CPU+GPU) high performance computing machines. Fierro can aid a) modeling and design efforts that have historically relied on commercial implicit and explicit finite element codes, b) numerical methods research, c) manufacturing research, and d) computer science research. The code contains diverse numerical methods to solve the governing physics equations for both quasi-static and dynamic problems. Mathematical optimization solvers are coupled to the numerical methods to research topology and shape optimization that has application to additive manufacturing, and to create novel numerical approaches. Phase-field methods with micromechanical solvers are provided to simulate microstructure formation and evolution in manufacturing processes. The micromechanical solvers can also help research efforts create continuum-scale constitutive models for solids, as a function of the microstructure, in situ in a calculation or in a stand-alone manner. No physical data exists within the code.

Morgan, Nathaniel↗
NASA NTRS2002

A Thermodynamic Theory Of Solid Viscoelasticity. Part 1: Linear Viscoelasticity.

The present series of three consecutive papers develops a general theory for linear and finite solid viscoelasticity. Because the most important object for nonlinear studies are rubber-like materials, the general approach is specified in a form convenient for solving problems important for many industries that involve rubber-like materials. General linear and nonlinear theories for non-isothermal deformations of viscoelastic solids are developed based on the quasi-linear approach of non-equilibrium thermodynamics. In this, the first paper of the series, we analyze non-isothermal linear viscoelasticity, which is applicable in a range of small strains not only to all synthetic polymers and bio-polymers but also to some non-polymeric materials. Although the linear case seems to be well developed, there still are some reasons to implement a thermodynamic derivation of constitutive equations for solid-like, non-isothermal, linear viscoelasticity. The most important is the thermodynamic modeling of thermo-rheological complexity , i.e. different temperature dependences of relaxation parameters in various parts of relaxation spectrum. A special structure of interaction matrices is established for different physical mechanisms contributed to the normal relaxation modes. This structure seems to be in accord with observations, and creates a simple mathematical framework for both continuum and molecular theories of the thermo-rheological complex relaxation phenomena. Finally, a unified approach is briefly discussed that, in principle, allows combining both the long time (discrete) and short time (continuous) descriptions of relaxation behaviors for polymers in the rubbery and glassy regions.

Freed, Alan D.↗
DOE OSTI2023

Fierro Version 2.x

FIERRO is a parallel C++ code designed to simulate fluid mechanics, heat transfer, and solid mechanics in two- and three-dimensional space. FIERRO is written to run on homogeneous (CPU) and heterogeneous (CPU+GPU) high performance computing machines. Fierro can aid a) modeling and design efforts that have historically relied on commercial implicit and explicit finite element codes, b) numerical methods research, c) manufacturing research, and d) computer science research. The code contains diverse numerical methods to solve the governing physics equations for both quasi-static and dynamic problems. Mathematical optimization solvers are coupled to the numerical methods to research topology and shape optimization that has application to additive manufacturing, and to create novel numerical approaches. Phase-field methods with micromechanical solvers are provided to simulate microstructure formation and evolution in manufacturing processes. The micromechanical solvers can also help research efforts create continuum-scale constitutive models for solids, as a function of the microstructure, in situ in a calculation or in a stand-alone manner. No physical data exists within the code.

Morgan, Nathaniel↗
NASA NTRS1999

Dynamics of the Molten Contact Line

The purpose of this program is to develop a basic understanding of how a molten material front spreads over a solid that is below its melting point, arrests, and freezes. Our hope is that the work will contribute toward a scientific knowledge base for certain new applications involving molten droplet deposition, including the "printing" of arbitrary three-dimensional objects by precise deposition of individual molten microdrops that solidify after impact. Little information is available at this time on the capillarity-driven motion and arrest of molten contact line regions. Schiaffino and Sonin investigated the arrest of the contact line of a molten microcrystalline wax spreading over a subcooled solid "target" of the same material. They found that contact line arrest takes place at an apparent liquid contact angle that depends primarily on the Stefan number S=c(T(sub f) -T(sub t)/L based on the temperature difference between the fusion point and the target temperature, and proposed that contact line arrest occurs when the liquid's dynamic contact angle approaches the angle of attack of the solidification front just behind the contact line. They also showed, however, that the conventional continuum equations and boundary conditions have no meaningful solution for this angle. The solidification front angle is determined by the heat flux just behind the contact line, and the heat flux is singular at that point. By comparing experiments with numerical computations, Schiaffino and Sonin estimated that the conventional solidification model must break down within a distance of order 0.1 - 1 microns of the contact line. The physical mechanism for this breakdown is as yet undetermined, and no first-principles theory exists for the contact angle at arrest. Schiaffino and Sonin also presented a framework for understanding how to moderate Weber number molten droplet deposition in terms of similarity laws and experimentation. The study is based on experiments with three molten materials- molten wax on solid wax, water on ice, and mercury on frozen mercury- which between them span a considerable range of the deposition/solidification similarity parameters. Correlations are obtained for the spreading velocity, spreading time scales, the spreading factor (i.e. ratio of deposited drop's final footprint radius and the drop's initial radius), post-spreading liquid oscillation amplitudes and time scales, and bulk solidification time scales. Duthaler carried out an experimental and theoretical investigation of the relationship between the liquid's apparent contact angle and the Capillary number Ca=mu U/sigma based on contact line speed, for molten materials spreading over subcooled solids. This relationship is required for modeling of melt spreading. We have adapted Voinov's methodology to the molten contact line and formulated a theoretical model for the Ca vs. contact angle relationship, based Schiaffino and Sonin#s (1997a,b) wedge-like solidification front model. With the solidification front angle taken from Schiaffino and Sonin, the model is in good agreement with the experimental results for Ca vs. contact angle. Duthaler also extended the experimental investigation of droplet deposition and contact line freezing to more materials, including solder on glass, solder on solder, water on ice, and molten microcrystalline wax on wax. The latter also included tests on inclined targets. Deposition tests have also been done with molten octacosane (C28H58) on various targets. An important objective of our program has been the development of micron-scale sensors for measuring the transient temperature at a point on the substrate surface as a molten contact line moves over it. The expectation is that this temperature history will yield a better understanding of the thermal process in the contact line region. The sensors are of the thermistor type, either 2.5 microns or 1.5 microns square, microfabricated with silicon-based technology on either pure silicon or amorphous silicon dioxide chips. Each chip has 32 sensors distributed on its surface in arrays. The time response is better than 10 ms. At the time of writing, sensor calibration is in progress. Results on thermal transients during contact line passage will be discussed at the conference. While we expect that the data will provide information on the near-contact-line heat transfer process, we also foresee possible problems. First, the spatial resolution of the sensors may be insufficient to resolve the near-contact-line region. Second, the sensors protrude about 0.5 microns above the substrate surface, and may affect the contact line motion. Third, a sensor's temperature history depends on both the heat flux distribution into it from the fusion front and the thermal properties of the substrate below it and the solidified melt between it and the fusion front. The heat flux distribution in the contact line region must therefore be unfolded from computations of the overall system's transient thermal response.

Sonin, Ain A.↗
DOE OSTI2026

Electron-impact ionization from excited states method

We present a distorted-wave with exchange (DWE) formulation for computing electron-impact ionization (EII) cross sections within the recently developed excited states method (ESM) for plasmas. The ESM provides a self-consistent quantum-mechanical description of both bound and continuum electronic states in dense plasmas, incorporating finite-temperature and screening effects absent in isolated-atom approaches. Using this framework, we calculate EII cross sections for lithium across a wide range of temperatures and densities. We show that the ESM reproduces the isolated-atom limit at low densities while capturing strong plasma effects, such as pressure ionization and shape resonances, at solid density. Comparisons with isolated-atom DWE and average-atom (AA) calculations reveal that, although AA cross sections can approximate configuration average behavior, the ESM provides more accurate, state-resolved answers. In conclusion, these results demonstrate that the ESM allows one to obtain rate coefficients for collisional-radiative modeling of plasmas out of local thermodynamical equilibrium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

A novel framework for elastoplastic behaviour of anisotropic solids

A general framework for developing nonlinear hyperelastic/plastic constitutive laws for anisotropic solids experiencing large strains and strain rates has been developed. The proposed framework does not rely on the “a priori” known strain energy function, but instead introduces a physical decomposition of the material element into seven physically independent stress bearing mechanisms, each of which has a constitutive law in terms of internal moments described by a scalar function of a single variable. Furthermore, the model has been encoded into a combined finite-discrete element method and tested against static geomechanical test data. The numerical validation experiments show the model can reproduce plastic anisotropic behavior in both biaxial and uniaxial loading of a geomaterial.

58 GEOSCIENCES↗
DOE OSTI2021

Situating the Vector Density Approach Among Contemporary Continuum Theories of Dislocation Dynamics

Abstract For the past century, dislocations have been understood to be the carriers of plastic deformation in crystalline solids. However, their collective behavior is still poorly understood. Progress in understanding the collective behavior of dislocations has primarily come in one of two modes: the simulation of systems of interacting discrete dislocations and the treatment of density measures of varying complexity that are considered as continuum fields. A summary of contemporary models of continuum dislocation dynamics is presented. These include, in order of complexity, the two-dimensional statistical theory of dislocations, the field dislocation mechanics treating the total Kröner–Nye tensor, vector density approaches that treat geometrically necessary dislocations on each slip system of a crystal, and high-order theories that examine the effect of dislocation curvature and distribution over orientation. Each of theories contain common themes, including statistical closure of the kinetic dislocation transport equations and treatment of dislocation reactions such as junction formation. An emphasis is placed on how these common themes rely on closure relations obtained by analysis of discrete dislocation dynamics experiments. The outlook of these various continuum theories of dislocation motion is then discussed.

Engineering↗
NASA NTRS1973

Theory of rheology

The structure of the modern theory of rheology is discussed to show the assumptions and limitations. Rheology is discussed as a branch of continuum mechanics to determine the relationships between stress, strain, and strain rate which will give a closer representation of lubricant properties than the Newtonian flow equation. Rheology is also investigated as a branch of chemical physics. Consideration is limited to those theories of nonpolymeric and polymeric fluids which can represent viscoelasticity in terms of identifiable and measureable molecular characteristics. The possibility that elastic liquids may rupture in shear and linear tension analogous to the failure of solids is proposed.

Hutton, J. F.↗
DOE OSTI2021

A High Resolution Simulation of a Single Shock-Accelerated Particle

We report particle drag models, which capture macroviscous and pressure effects, have been developed over the years for various flow regimes to enable cost effective simulations of particle-laden flows. The relatively recent derivation by Maxey and Riley has provided an exact equation of motion for spherical particles in a flow field based on the continuum assumption. Many models that have been simplified from these equations have provided reasonable approximations; however, the sensitivity of particle-laden flows to particle drag requires a very accurate model to simulate. To develop such a model, a two-dimensional axisymmetric Navier–Stokes direct numerical simulation of a single particle in a transient, shock-driven flow field was conducted using the hydrocode FLAG. FLAGs capability to run arbitrary Lagrangian-Eulerian hydrodynamics coupled with solid mechanic models makes it an ideal code to capture the physics of the flow field around and in the particle as it is shock-accelerated—a challenging regime to study. The goal of this work is twofold: to provide a validation for FLAGs Navier–Stokes and heat diffusion solutions and to provide a rationale for recent experimental particle drag measurements.

42 ENGINEERING↗
DOE OSTI2020

A Numerical Simulation of a Single Shock-Accelerated Particle

Particle drag models, which capture macro viscous and pressure effects, have been developed over the years for various flow regimes to enable cost effective simulations of particle-laden flows. The relatively recent derivation by Maxey and Riley has provided an exact equation of motion for spherical particles in a flow field based on the continuum assumption. Many models that have been simplified from these equations have provided reasonable approximations; however, the sensitivity of particle-laden flows to particle drag requires a very accurate model to simulate. To develop such a model, a 2D axisymmetric Navier-Stokes direct numerical simulation of a single particle in a transient, shock-driven flow field was conducted in the hydrocode FLAG. FLAG’s capability to run arbitrary Lagrangian-Eulerian hydrodynamics coupled with solid mechanic models makes it an ideal code to capture the physics of the flow field around and in the particle as it is shock-accelerated – a challenging regime to study. The goal of this work is twofold: to provide a validation for FLAG’s Navier-Stokes and heat diffusion solutions, and to provide a rationale for recent experimental particle drag measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Data-driven nonlocal model for fragmentation in the crushing of solids

A technique is proposed for reproducing particle size distributions in three-dimensional simulations of the crushing and comminution of solid materials. The method is designed to produce realistic distributions over a wide range of loading conditions, especially for small fragments. In contrast to most existing methods, the new model does not explicitly treat the small-scale process of fracture. Instead, it uses measured fragment distributions from laboratory tests as the basic material property that is incorporated into the algorithm, providing a data-driven approach. The algorithm is implemented within a nonlocal peridynamic solver, which simulates the underlying continuum mechanics and contact interactions between fragments after they are formed. Finally, the technique is illustrated in reproducing fragmentation data from drop weight testing on sandstone samples.

58 GEOSCIENCES↗
DOE OSTI2020

Soft X-ray and EUV emission spectra of beryllium plasma produced by neodymium-glass laser radiation with broad frequency and angular spectra

We present the results of an experimental study of soft X-ray (SXR) and extreme ultraviolet emission spectra of the plasma produced by exposing a plane solid beryllium target to laser radiation with broad frequency and angular spectra. SXR lines up to 1s – 9p of Be IV as well as the plasma continuum are recorded for a laser focal-spot intensity of 5.3 × 10{sup 13} W cm{sup −2}. To model the SXR beryllium plasma spectra, simulations are carried out using the INDHAUS programme and the FLYCHK code in the framework of local thermodynamic equilibrium model, which agree nicely with experimentally obtained results. (laser plasma)

36 MATERIALS SCIENCE↗
DOE OSTI2020

A mathematical model of solid-state dewetting of barium thin films on W (112)

Solid state dewetting occurs when a thin solid film is heated. The temperature of dewetting depends on the thickness of the film; dewetting can be observed in the range of 1/3 to 1/2 of the bulk melting temperature. While remaining solid, the film behaves in a manner similar to liquids dewetting and agglomerating to forming islands or droplets. One of the possible mechanisms is the conversion of a metastable thin film geometry into a more stable form. Heating the metastable film gives the film atoms higher mobility, and the films retract, dewetting the surface. This atomic motion can be restricted due to surface anisotropy. Herein we present in situ emission microscopy observations of barium thin films deposited onto W (112) by thermal evaporation. From the modeling viewpoint, the evolution of the film in this system could be divided in four stages: (i) the nucleation and growth of the thin film as a simply connected region; (ii) formation of droplets/islands/hillocks; (iii) nucleation of holes; (iv) evolution of the components of the disconnected film to their equilibrium state. We present a continuum model that is qualitatively consistent with the evolution of the film observed at the initial stage of the experiment and discuss the later stages of the evolution of surface structures.

36 MATERIALS SCIENCE↗