Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computing continuum”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Impact of artificial topological changes on flow and transport through fractured media due to mesh resolution

Abstract We performed a set of numerical simulations to characterize the interplay of fracture network topology, upscaling, and mesh refinement on flow and transport properties in fractured porous media. We generated a set of generic three-dimensional discrete fracture networks at various densities, where the radii of the fractures were sampled from a truncated power-law distribution, and whose parameters were loosely based on field site characterizations. We also considered five network densities, which were defined using a dimensionless version of density based on percolation theory. Once the networks were generated, we upscaled them into a single continuum model using the upscaled discrete fracture matrix model presented by Sweeney et al. (2019). We considered steady, isothermal pressure-driven flow through each domain and then simulated conservative, decaying, and adsorbing tracers using a pulse injection into the domain. For each simulation, we calculated the effective permeability and solute breakthrough curves as quantities of interest to compare between network realizations. We found that selecting a mesh resolution such that the global topology of the upscaled mesh matches the fracture network is essential. If the upscaled mesh has a connected pathway of fracture (higher permeability) cells but the fracture network does not, then the estimates for effective permeability and solute breakthrough will be incorrect. False connections cannot be eliminated entirely, but they can be managed by choosing appropriate mesh resolution and refinement for a given network. Adopting octree meshing to obtain sufficient levels of refinement leads to fewer computational cells (up to a 90% reduction in overall cell count) when compared to using a uniform resolution grid and can result in a more accurate continuum representation of the true fracture network.

58 GEOSCIENCES↗
DOE OSTI2021

Data driven learning of robust nonlocal models

Nonlocal models use integral operators that embed length-scales 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 nonlocal model as a surrogate for MD displacements. Our 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. The efficacy of this approach is demonstrated with several numerical tests for single layer graphene both in the case of perfect crystal and in the presence of thermal noise.

36 MATERIALS SCIENCE↗
DOE OSTI2022

On the Viability of Quantum Annealers to Solve Fluid Flows

This paper explores the suitability of upcoming novel computing technologies, particularly adiabatic annealing based quantum computers, to solve fluid dynamics problems that form a critical component of several science and engineering applications. For our experiments, we start with a well-studied one-dimensional simple flow problem, and provide a framework to convert such problems in continuum to a form amenable for deployment on such quantum annealers. Since the DWave annealer returns multiple states sampling the energy landscape of the problem, we explore multiple solution selection strategies to approximate the solution of the problem. We analyze the continuum solutions obtained both qualitatively and quantitatively as well as their sensitivities to the particular solution selection scheme.

42 ENGINEERING↗
NASA NTRS2006

Blunt Body Aerodynamics for Hypersonic Low Density Flows

Numerical simulations are performed for the Apollo capsule from the hypersonic rarefied to the continuum regimes. The focus is on flow conditions similar to those experienced by the Apollo 6 Command Module during the high altitude portion of its reentry. The present focus is to highlight some of the current activities that serve as a precursor for computational tool assessments that will be used to support the development of aerodynamic data bases for future capsule flight environments, particularly those for the Crew Exploration Vehicle (CEV). Results for aerodynamic forces and moments are presented that demonstrate their sensitivity to rarefaction; that is, free molecular to continuum conditions. Also, aerodynamic data are presented that shows their sensitivity to a range of reentry velocities, encompassing conditions that include reentry from low Earth orbit, lunar return, and Mars return velocities (7.7 to 15 km/s). The rarefied results obtained with direct simulation Monte Carlo (DSMC) codes are anchored in the continuum regime with data from Navier-Stokes simulations.

Moss, James N.↗
NASA NTRS1995

Some Developments of the Equilibrium Particle Simulation Method for the Direct Simulation of Compressible Flows

The direct simulation Monte Carlo (DSMC) method is the established technique for the simulation of rarefied gas flows. In some flows of engineering interest, such as occur for aero-braking spacecraft in the upper atmosphere, DSMC can become prohibitively expensive in CPU time because some regions of the flow, particularly on the windward side of blunt bodies, become collision dominated. As an alternative to using a hybrid DSMC and continuum gas solver (Euler or Navier-Stokes solver) this work is aimed at making the particle simulation method efficient in the high density regions of the flow. A high density, infinite collision rate limit of DSMC, the Equilibrium Particle Simulation method (EPSM) was proposed some 15 years ago. EPSM is developed here for the flow of a gas consisting of many different species of molecules and is shown to be computationally efficient (compared to DSMC) for high collision rate flows. It thus offers great potential as part of a hybrid DSMC/EPSM code which could handle flows in the transition regime between rarefied gas flows and fully continuum flows. As a first step towards this goal a pure EPSM code is described. The next step of combining DSMC and EPSM is not attempted here but should be straightforward. EPSM and DSMC are applied to Taylor-Couette flow with Kn = 0.02 and 0.0133 and S(omega) = 3). Toroidal vortices develop for both methods but some differences are found, as might be expected for the given flow conditions. EPSM appears to be less sensitive to the sequence of random numbers used in the simulation than is DSMC and may also be more dissipative. The question of the origin and the magnitude of the dissipation in EPSM is addressed. It is suggested that this analysis is also relevant to DSMC when the usual accuracy requirements on the cell size and decoupling time step are relaxed in the interests of computational efficiency.

Macrossan, M. N.↗
DOE OSTI2021

Origin of the background absorption in carbon nanotubes: Phonon-assisted excitonic continuum

Excitonic effects in 1D semiconductors can be qualitatively different from those in higher dimensions. In particular, the Sommerfeld factor, the ratio of the above-band-edge excitonic continuum absorption to free electron-hole pair generation, has been shown to be less than 1 (suppressed) in 1D systems while it is larger than 1 (enhanced) in 2D and 3D systems. Strong continuum suppression indeed exists in semiconducting single-wall carbon nanotubes, a prototypical 1D semiconductor. However, absorption spectra for carbon nanotubes are typically fit with a combination of Lorentzians and a polynomial background baseline with little physical meaning. We performed absorption measurements in aligned single-chirality (6,5) carbon nanotube films. The obtained spectra were fit with our theoretical model obtained by solving the Boltzmann scattering equation involving fifty-nine different types of transitions among three different types of quasiparticles. Specifically, we took into account microscopic interactions between photons, phonons, and excitons, including their dispersions, which unambiguously demonstrated that the background absorption is due to phonon-assisted transitions from the semiconductor vacuum to finite-momentum continuum states of excitons. The excellent agreement we obtained between experiment and theory suggests that our numerical technique can be seamlessly extended to compute strongly out-of-equilibrium many-body dynamics and time-resolved spectra in low-dimensional materials.

36 MATERIALS SCIENCE↗
NASA NTRS1994

Absolute calibration and atmospheric versus mineralogic origin of absorption features in 2.0 to 2.5 micron Mars spectra obtained during 1993

We obtained new high resolution reflectance spectra of Mars during the 1993 opposition from Mauna Kea Observatory using the UKIRT CGS4 spectrometer. Fifty spectra of 1600-2000 km surface regions and a number of standard star spectra were obtained in the 2.04 to 2.44 micron wavelength region on 4 February 1993 UT. Near-simultaneous observations of bright standard stars were used to perform terrestrial atmospheric corrections and an absolute flux calibration. Using the known magnitude of the stars and assuming blackbody continuum behavior, the flux from Mars could be derived. A radiative transfer model and the HITRAN spectral line data base were used to compute atmospheric transmission spectra for Mars and the Earth in order to simulate the contributions of these atmospheres to our observed data. Also, we examined the ATMOS solar spectrum in the near-IR to try to identify absorption features in the spectrum of the Sun that could be misinterpreted as Mars features. Eleven absorption features were detected in our Mars spectra. Our data provide no conclusive identification of the mineralogy responsible for the absorption features we detected. However, examination of terrestrial spectral libraries and previous high spectral resolution mineral studies indicates that the most likely origin of these features is either CO3(sup 2-), HCO3(-), or HSO4(-) anions in framework silicates or possibly (Fe, Mg)-OH bonds in sheet silicates.

Bell, James F., III↗
NASA NTRS2009

Simulation of RCC Crack Growth Due to Carbon Oxidation in High-Temperature Gas Environments

The carbon wall oxidation technique coupled with a CFD technique was employed to study the flow in the expanding crack channel caused by the oxidation of the channel carbon walls. The recessing 3D surface morphing procedure was developed and tested in comparison with the arcjet experimental results. The multi-block structured adaptive meshing was used to model the computational domain changes due to the wall recession. Wall regression rates for a reinforced carbon-carbon (RCC) samples, that were tested in a high enthalpy arcjet environment, were computationally obtained and used to assess the channel expansion. The test geometry and flow conditions render the flow regime as the transitional to continuum, therefore Navier-Stokes gas dynamic approach with the temperature jump and velocity slip correction to the boundary conditions was used. The modeled mechanism for wall material loss was atomic oxygen reaction with bare carbon. The predicted channel growth was found to agree with arcjet observations. Local gas flow field results were found to affect the oxidation rate in a manner that cannot be predicted by previous mass loss correlations. The method holds promise for future modeling of materials gas-dynamic interactions for hypersonic flight.

Titov, E. V.↗
DOE OSTI2023

COWALKER:EFFECTIVE TRANSPORT PROPERTIES OF COMPOSITE MATERIALS

SF-23-026 This software computes effective transport properties of composite materials involving fibers and nanoparticles using a random-walk algorithm that efficiently scales to an arbitrary number of processes and cores. Effective transport properties (thermal, electrical) are key to bridge the microstructure of complex materials with its macroscopic behavior. Traditional approaches either use effective medium approximations (closed mathematical expressions that are approximation for certain conditions) or continuum simulation models such as finite element or finite volume, which require the generation of a mesh for each configuration explored. cowalker leverages the equivalence between laplacian or heat equation-based models and random walks to compute the asymptotic transport properties from an ensemble of first sojourn times of a random walker moving through the composite material. This allows us to directly define a composite material as a collection of particles and use algorithms developed for molecular dynamics to quickly compute the intersection of the walker with the different interfaces in the material. cowalker is developed in C++, and it relies on the GNU Scientific Library for random generation. cowalker is currently delivered as source code, so the GSL library is not included in cowalker's distribution. A more userfriendly version, cowalker.jl is currently in development and will be released as part of cowalker.

YANGUAS-GIL, ANGEL↗
NASA NTRS1999

Some Continuum Aspects of Data Assimilation

A long-sought goal in data assimilation is to build more fully the time dimension into the data assimilation process. Truly four-dimensional data assimilation requires evolving second-moment information, namely the estimation error covariance, along with the state estimate itself. The continuum evolution equations for the estimation error covariance constitute a system of partial differential equations in six space dimensions, forced by a model error covariance and by the observations themselves. The high dimensionality of this system poses a difficult computational problem. Numerous methods have been proposed for approximate, discrete solution of this system of equations. Four-dimensional variational schemes solve these equations indirectly, while schemes based on Kalman filtering solve the equations more directly, usually on a small subspace of the full six-dimensional space. In both cases, most proposed solution methods are not derived from the continuum covariance evolution equations themselves. Instead, most methods simply inherit the discretization used for evolving the state estimate, either in the form of a discrete tangent linear model or, in the case of ensemble schemes, in the form of a fully nonlinear discrete model. In this lecture we show that solutions of the continuum covariance evolution equations possess simple properties that are not readily satisfied upon discretization of the equations. we give numerical examples illustrating that "inherited" discretizations sometimes fail to produce meaningfully accurate solutions. Finally, we suggest alternative discretization methods that may yield more faithful approximate solutions of the continuum problem.

Cohn, Stephen E.↗
DOE OSTI2022

Crossing Versus Locking: Bit Threads and Continuum Multiflows

Bit threads are curves in holographic spacetimes that manifest boundary entanglement, and are represented mathematically by continuum analogues of network flows or multiflows. Subject to a density bound, the maximum number of threads connecting a boundary region to its complement computes the Ryu–Takayanagi entropy. When considering several regions at the same time, for example in proving entropy inequalities, there are various inequivalent density bounds that can be imposed. Here, we investigate for which choices of bound a given set of boundary regions can be “locked”, in other words can have their entropies computed by a single thread configuration. We show that under the most stringent bound, which requires the threads to be locally parallel, non-crossing regions can in general be locked, but crossing regions cannot, where two regions are said to cross if they partially overlap and do not cover the entire boundary. We also show that, under a certain less stringent density bound, a crossing pair can be locked, and conjecture that any set of regions not containing a pairwise crossing triple can be locked, analogously to the situation for networks.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Estimation of micro-Hall-Petch coefficients for prismatic slip system in Mg-4Al as a function of grain boundary parameters

Grain size strengthening, referred to as the Hall-Petch effect, is a common strategy to improve the yield strength of magnesium (Mg) alloys. Several experimental studies have reported that the Hall-Petch slope strongly depends on the texture of the alloy. This effect arises from altering grain boundaries (GBs) resistance to different slip systems to transfer across adjacent grains. The grain boundary barrier strength of certain grain boundaries to basal slip, referred to as basal micro-Hall-Petch, was investigated in the previous work. Here in this study, the micro-Hall-Petch coefficient values for the prismatic slip $(k^{prismatic}_{μ})$ in Mg-4Al and their correlation with the grain boundary parameters were investigated. An experimental method was developed to initiate the prismatic slip band at low-stress levels. High-resolution electron backscatter diffraction (HR-EBSD) was used to measure the residual stress tensor, from which the resolved shear stress ahead of blocked prismatic slip bands was computed for seven different grain boundaries. $(k^{prismatic}_{μ})$ values for each individual GB were calculated by coupling the stress profile information with a continuum dislocation pile-up model. The $(k^{prismatic}_{μ})$ values vary from 0.138 MPa.m 1/2 to 0.685 MPa. m 1/2 which are almost three times larger than the calculated values for the basal micro-Hall-Petch. The $(k^{prismatic}_{μ})$ values were correlated with the GB parameters, and a functional relationship depending on the two most effective angles, the angle between the traces of the slip planes on the GB plane ($\phi$) and the angle between incoming and outgoing slip directions (κ), was proposed to estimate the Hall-Petch barrier for prismatic slip system. The work provides coefficients that can be supplied as input to crystal plasticity models to couple the effect of texture and grain size effectively.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Modeling electrokinetic flows with the discrete ion stochastic continuum overdamped solvent algorithm

In this article we develop an algorithm for the efficient simulation of electrolytes in the presence of physical boundaries. In previous work the discrete ion stochastic continuum overdamped solvent (DISCOS) algorithm was derived for triply periodic domains, and was validated through ion-ion pair correlation functions and Debye-Hückel-Onsager theory for conductivity, including the Wien effect for strong electric fields. In extending this approach to include an accurate treatment of physical boundaries we must address several important issues. First, the modifications to the spreading and interpolation operators necessary to incorporate interactions of the ions with the boundary are described. Next we discuss the modifications to the electrostatic solver to handle the influence of charges near either a fixed potential or dielectric boundary. An additional short-ranged potential is also introduced to represent interaction of the ions with a solid wall. Finally, the dry diffusion term is modified to account for the reduced mobility of ions near a boundary, which introduces an additional stochastic drift correction. Herein, several validation tests are presented confirming the correct equilibrium distribution of ions in a channel. Additionally, the methodology is demonstrated using electro-osmosis and induced-charge electro-osmosis, with comparison made to theory and other numerical methods. Notably, the DISCOS approach achieves greater accuracy than a continuum electrostatic simulation method. We also examine the effect of under-resolving hydrodynamic effects using a “dry diffusion” approach, and find that considerable computational speedup can be achieved with a negligible impact on accuracy.

97 MATHEMATICS AND COMPUTING↗
NASA NTRSDate not supplied

TPSAS-NF1676L-13792-DND

Extreme structural requirements of future aerospace vehicles motivate the development of new, ultra-durable materials and game-changing methodologies for material certification and sustainment. Computational simulations spanning many orders of magnitude in length and time scales—from the nanoscale of underlying damage processes to the larger scales of continuum cracks—are being developed to support these requirements. Results of these simulations are used to deduce key aspects of material response to loads and environments, including: - The microstructural mechanics that govern metal fatigue crack initiation and growth - The variety and complexity of the dislocation-precipitate interac- tions that underpin plastic behavior and damage evolution - The energetic principles that govern the interaction of water with crack surfaces within a structural component - The mechanics of grain boundary separation Modeling the physics behind these metal fatigue behaviors is extremely complex, and requires intensive, high-fidelity simulations. The knowledge gained through these simulations, however, is forming the keystone for advanced material design and will enable development of more capable, reliable structures for aerospace vehicles.

Jacob Hochhalter↗
DOE OSTI2021

SU(n) and Quantum SU(n) Symmetries in Physical Systems [Slides]

Presence of SU(n) or other Lie group symmetry in a physical system is its powerful, usually underutilized property. In many cases it allows for finding analytical solutions to nonlinear differential equations describing this system. Power of the method is presented on diversified examples from mathematical physics: Lie-group symmetries in finding solutions of generalized, multidimensional theory of gravity; analytical Dirac–equation solutions for description of conducting polymers; stability of qubit states in quantum computers; spatial defects in condensed matter; reconstruction of 3D object from its 2D tomographic image; significant improvement of numerical solutions stability for Euler equations. The next question after obtaining such Lie group symmetric solution is: does a generalized solution with appropriate quantum group symmetry exists for the given physical system, and if yes what is the physical meaning of the deformation parameter q introduced by such solution. In many cases it can be identified. Any SU(n) solution is by its nature singular, assuming a perfect symmetry of the physical system discussed. Such solution gives a powerful insight to theoretical physics, yet the assumption may be too demanding for experimental applications. Deformation parameter q from a quantum group symmetry allows for a continuum of solutions, more applicable to experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Strongly nonlinear wave propagation in elasto-plastic metamaterials: Low-order dynamic modeling

Nonlinear elastic metamaterials are known to support a variety of dynamic phenomena that enhance our capacity to manipulate elastic waves. Since these properties stem from complex, subwavelength geometry, full-scale dynamic simulations are often prohibitively expensive at scales of interest. Prior studies have therefore utilized low-order effective medium models, such as discrete mass-spring lattices, to capture essential properties in the long-wavelength limit. While models of this type have been successfully implemented for a wide variety of nonlinear elastic systems, they have predominantly considered dynamics depending only on the instantaneous kinematics of the lattice, neglecting history-dependent effects, such as wear and plasticity. Here, to address this limitation, the present study develops a lattice-based modeling framework for nonlinear elastic metamaterials undergoing plastic deformation. Due to the history- and rate-dependent nature of plasticity, the framework generally yields a system of differential-algebraic equations whose computational cost is significantly greater than an elastic system of comparable size. We demonstrate the method using several models inspired by classical lattice dynamics and continuum plasticity theory and explore means to obtain empirical plasticity models for general geometries, thereby gaining insight into the influence of microstructural plasticity on effective material performance, which can be used to improve the design of nonlinear mechanical metamaterials.

Dynamic simulation↗
NASA NTRS1975

Fundamental analysis of the failure of polymer-based fiber reinforced composites

A mathematical model predicting the strength of unidirectional fiber reinforced composites containing known flaws and with linear elastic-brittle material behavior was developed. The approach was to imbed a local heterogeneous region surrounding the crack tip into an anisotropic elastic continuum. This (1) permits an explicit analysis of the micromechanical processes involved in the fracture, and (2) remains simple enough to be useful in practical computations. Computations for arbitrary flaw size and orientation under arbitrary applied loads were performed. The mechanical properties were those of graphite epoxy. With the rupture properties arbitrarily varied to test the capabilities of the model to reflect real fracture modes, it was shown that fiber breakage, matrix crazing, crack bridging, matrix-fiber debonding, and axial splitting can all occur during a period of (gradually) increasing load prior to catastrophic failure. The calculations also reveal the sequential nature of the stable crack growth process proceding fracture.

Kanninen, M. F.↗
NASA NTRS1977

Atmospheric interaction plume

Computer program can be used to predict gas dynamic and chemical properties of underexpanded rocket plumes from sea level to altitude above which viscous continuum-flow assumption with distinct shocks is no longer valid.

Kelly, J. T.↗