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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Recent Developments to the Porous Microstructure Analysis (PuMA) Software

The Porous Microstructure Analysis (PuMA) software is an open source framework for image-based simulation, primarily used to determine effective properties based on material microstructure. PuMA was originally developed for the study of NASA thermal protection materials; however, many of the solvers in PuMA have applicability to a broad range of materials science applications. PuMA version 3.2 computes material surface area, pore diameters, effective thermal conductivity, continuum and rarefied tortuosity, and permeability. For anisotropic materials, PuMA can estimate material orientation and compute anisotropic thermal conductivity and elasticity. In this talk, a brief overview of the PuMA software and underlying methods will be presented, as well as some recent and ongoing developments, including the use of immersed boundary methods for image-based simulation and the development of a new weave segmentation tool, called TomoSAM. Cut-Cell method for heat and mass transfer For simulations on complex microstructures, traditional unstructured meshing techniques often prove to be difficult and time-intensive. Voxel-based solvers, which represent the surface as a staircase structure, are relatively simple to implement but can lose accuracy when feature resolution is poor. In this work, we present a novel 3D cut-cell method for solving the variable coefficient Poisson equation on complex microstructures, suitable for the determination of effective thermal conductivity or tortuosity of a material. The method uses a Marching Cubes/Marching Squares surface reconstruction to create cut-cells and determine geometric quantities. A flux-correction method is extended to 3D, with least squares gradient reconstruction, to solve for the boundary fluxes in the cut-cells. Verification cases show the solver achieves globally 2nd order accuracy on complex microstructures. TomoSAM TomoSAM, a module of the PuMA software, has been developed as a plugin for 3D Slicer, a software platform used for 3D image processing and visualization. It utilizes the Segment Anything Model (SAM), a deep learning model capable of identifying objects and generating image masks based on minimal user input. This feature enables efficient segmentation of complex 3D datasets, particularly of woven materials, from tomography or similar imaging methods, reducing the need for manual segmentation.

Tomography↗
NASA NTRS1980

Comments on stellar boundary cooling and the reality of supermetallicity

The paper discusses the 'super-metal-rich' (SMR) stars and reexamines Peterson's analysis of the SMR prototype mu Leo (1978) with regard to a postulated error in continuum error. Model atmospheres are used to compute theoretical equivalent widths and to explore the sensitivity of these widths to metallicity, temperature, surface gravity, and microturbulence. It is shown that Peterson's results are sensitive to continuum placement, and that her data does not indicate that the temperature gradient is steeper in mu Leo than in normal giants. It is concluded that the SMR stars are very metal rich and are also somewhat boundary cooled, possibly due to high metallicity.

Deming, D.↗
NASA NTRS1986

Preliminary development of a global 3-D magnetohydrodynamic computational model for solar wind-cometary and planetary interactions

This is the final summary report by Resource Management Associates, Inc., of the first year's work under Contract No. NASW-4011 to the National Aeronautics and Space Administration. The work under this initial phase of the contract relates to the preliminary development of a global, 3-D magnetohydrodynamic computational model to quantitatively describe the detailed continuum field and plasma interaction process of the solar wind with cometary and planetary bodies throughout the solar system. The work extends a highly-successful, observationally-verified computational model previously developed by the author, and is appropriate for the global determination of supersonic, super-Alfvenic solar wind flows past planetary obstacles. This report provides a concise description of the problems studied, a summary of all the important research results, and copies of the publications.

Stahara, S. S.↗
DOE OSTI2022

Entanglement entropy across the lattice-continuum correspondence

This paper revisits the standard calculations of free field entanglement entropy in light of the newly developed lattice-continuum correspondence. This correspondence prescribes an explicit method to extract an approximately continuum quantum field theory out of a fully regularized lattice theory. This prescription will here be extended to subregion algebras, and it will be shown how entropies of continuum boson and fermion theories can be computed by working purely with lattice quantities. This gives a clear picture of the origin of divergences in entanglement entropy while also presenting a concise and detailed recipe for calculating this important quantity in continuum theories.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2023

Statistical Study of the Non-Thermal Continuum Radiation Beaming Angle Measured By the High Frequency Receiver on Van Allen Probes-A

The nonthermal continuum (NTC) radiation beaming angle is computed over the entire Van Allen Probes-A mission when the spacecraft was in the dawn sector. The conditions in the dawn sector are favorable for the wave vector to lie near/in the spacecraft's spin plane allowing a favorable estimate of the beaming angle, and the dawn sector is also advantageous in that previous studies show NTC occurrence to peak in this sector. We found that scatter plots, over the entire mission, of beaming angle versus magnetic latitude form a distinct inverted-V pattern, with the apex at/near the magnetic equator. This pattern was sharpest for frequencies (f) ≲ 100 kHz. Using the NTC beaming formula from the linear mode conversion theory (LMCT), we show that such an inverted-V pattern is expected due to the large variation in the plasmapause location over the entire mission. The theoretical derived pattern qualitatively reproduces the observed pattern but not quantitatively. The lack of quantitative agreement is discussed and is attributed to several factors, one factor is off-centered emissions from the radio window. The qualitative agreement strongly supports LMCT as being the dominant mechanism generating NTC for f ≲ 100 kHz. For f ≳ 100 kHz, the inverted-V pattern becomes less distinct, and strong near-equatorial beaming is observed. After considering contamination of our selections by left-handed polarized AKR, our study suggests that besides LMCT another unidentified NTC generation mechanism becomes important for f ≳ 100 kHz.

Linear Mode Conversion Theory↗
DOE OSTI2023

Statistical Study of the Nonthermal Continuum Radiation Beaming Angle Measured by the High Frequency Receiver on Van Allen Probes-A

The nonthermal continuum (NTC) radiation beaming angle is computed over the entire Van Allen Probes-A mission when the spacecraft was in the dawn sector. The conditions in the dawn sector are favorable for the wave vector to lie near/in the spacecraft's spin plane allowing a favorable estimate of the beaming angle, and the dawn sector is also advantageous in that previous studies show NTC occurrence to peak in this sector. We found that scatter plots, over the entire mission, of beaming angle versus magnetic latitude form a distinct inverted-V pattern, with the apex at/near the magnetic equator. This pattern was sharpest for frequencies ( f ) ≲ 100 kHz. Using the NTC beaming formula from the linear mode conversion theory (LMCT), we show that such an inverted-V pattern is expected due to the large variation in the plasmapause location over the entire mission. The theoretical derived pattern qualitatively reproduces the observed pattern but not quantitatively. The lack of quantitative agreement is discussed and is attributed to several factors, one factor is off-centered emissions from the radio window. The qualitative agreement strongly supports LMCT as being the dominant mechanism generating NTC for f ≲ 100 kHz. For f ≳ 100 kHz, the inverted-V pattern becomes less distinct, and strong near-equatorial beaming is observed. Finally, after considering contamination of our selections by left-handed polarized AKR, our study suggests that besides LMCT another unidentified NTC generation mechanism becomes important for f ≳ 100 kHz.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Multiscale modeling-enabled design of multifunctional composites

This study aims to create a comprehensive model that considers multiple scales and physics for predicting the electromechanical behavior of fiber-reinforced composites enhanced with barium titanate (BaTiO3). In our earlier work, we have demonstrated that depositing BaTiO3 microparticles of 200-nm-diameter, on fiber surfaces during fiber-reinforced composite fabrication enhances mechanical strength, passive self-sensing, and energy harvesting properties. The key is to carefully control the microparticle concentration to prevent agglomeration. Since the particles are micron-sized, understanding how agglomeration affects the composites' electromechanical properties is crucial for guiding such multifunctional materials’ design. This study introduces a micromechanics-based approach to explore the impact of microparticle dispersion on the bulk composites' electromechanical properties. Insights gained from this investigation are applied in experiments, enabling accurate predictions of mechanical and self-sensing responses in BaTiO3-enhanced fiber-reinforced composites. Micro-level findings from this computational approach can be integrated into larger continuum models to comprehensively capture the electromechanical behavior of the composite structures at bulk scale. The proposed model is validated by comparing predictions with experimental results, accounting for the nonlinear mechanical and electromechanical behaviors of constituent materials. Consequently, this computational model serves as a digital platform for efficiently designing multifunctional composites.

Gupta, Sumit↗
NASA NTRSDate not supplied

Recommended DSMC Collision Model Parameters for Planetary Entry

Hybrid simulation frameworks combining Computational Fluid Dynamics (CFD) and Direct Simulation Monte Carlo (DSMC) are frequently employed to efficiently perform high-fidelity simulations of environments consisting of both continuum and rarified flow. DSMC is a stochastic, particle-based method which solves the fundamental Boltzmann equation and is therefore necessary for high-Knudsen flow where continuum-based assumptions governing CFD break down. However, the DSMC methodology is generally computationally inefficient to model the continuum regime. In a CFD/DSMC hybrid approach, obtaining an accurate, high-fidelity solution hinges on the consistent treatment of transport properties and the thermo-chemical models employed within the two solvers. In principle, in regions where CFD and DSMC are both employed, the same gas mixture under the same conditions should have the same properties, regardless of simulation type. Observed differences should be due to non-equilibrium processes, rather than differences in physical models. The goal of this work is to provide a comprehensive DSMC transport database encompassing the vast majority of known gas species encountered during Earth or other planetary atmospheric entry. This goal is largely possible due to recently performed ab initio quantum chemistry calculations. Combined with other high-fidelity data, the planned database will consist of collision integral data for over 200 neutral and ionized species and over 20000 binary collisions. From these collision integrals, Nelder-Mead optimization is used to compute collision-specific Variable Soft Sphere (VSS) collision model parameters, fitted from 300 K to 20000 K. Initial comparisons of transport properties of relevant equilibrium gas mixtures show great agreement between CFD and DSMC-derived results. The completed database can be readily applied to model binary collisions of any gas mixture containing the included species over the specified temperature range, making it a valuable tool for future planetary probe modeling efforts. An example is shown below in Fig. 1. Equilibrium mixture transport properties for a 35-species mixture composed originally of 10% air and 90% pyrolysis species of a carbon-phenolic ablator material [4] are computed using both fitted VSS Parameters and the original CFD collision integral values. Deviations in computed properties between the two approaches are less than 5% for the entire temperature range.

M. R. Gosma↗
DOE OSTI2020

Softwarized Federations of Science Instruments with Edge-Continuum Containers

Significant expansion of capabilities of DOE science complex of supercomputers, instruments and networks is expected as powerful experimental facilities, exascale computers and terabit networks are added. Combined with the advances in edge and cloud computing technologies, DOE science users now have the promise of unprecedented execution of complex, continuum workflows, namely, small and latency-sensitive computations at the edge and site, massive computations at remote HPC systems, and everything in between on the cloud. But, bringing this capability to the science user requires overcoming the overwhelming complexity of forming the federations of systems, and efficiently and effectively orchestrating the workflows while ensuring high utilization of the expensive facilities. Current manual configuration of the federated systems simply will not scale, since the coordination across sites may take weeks to months, often leading to under-utilized and hard-to-diagnose compositions. A powerful, composable software stack will be developed to (i) wrap the systems so that federations can be composed fast in software, and (ii) containerize computations to be orchestrated across the edge, site, cloud and HPC resources.

Rao, Nageswara S.↗
DOE OSTI2020

Deformation accommodating periodic computational domain for a uniform velocity gradient

Many multiscale methods for granular materials use periodic computational domains to consider particle scale interactions and then calculate the stress to drive the continuum scale calculations. For problems involving large material deformations, the computation domain often needs to be reinitialized because of the distortion, causing the loss of the history information of the system. This work introduces an algorithm to accommodate a large deformation of the material while maintaining the computational domain cuboid to avoid the domain reinitialization during the computation. The algorithm uses a rotating frame of reference, in which the velocity gradient can be represented by an upper triangular matrix. The deformation caused by the upper triangular matrix is treated by the image system implied in the periodicity to maintain the computational domain cuboid. The effect from the rotation of the reference frame is considered using the inertial forces. Finally, simulations of simple and pure shear motions are carried out to illustrate the algorithm.

42 ENGINEERING↗
DOE OSTI2023

Multiphase Modeling in a Parallel Plate Fluidized Bed Receiver for Concentrating Solar Power

A novel high temperature particle solar receiver is developed by using a light trapping planar cavity configuration. As particles fall through the cavity, the concentrated solar radiation warms the boundaries of the receiver and in turn heats the particles. Particles flow through the system, forming a packed bed at the lower end, leaving the system from the bottom at a constant flow rate. Air is introduced to the system as the fluidizing medium to improve particle heat transfer and mixing. A laboratory scale cavity receiver is built and a near IR quartz lamp is used to provide flux to the vertical wall of the heat exchanger. The system is modeled using a continuum two-fluid method. The computational model matches the experimental system size and the particle size distribution is assumed monodisperse. A conduction model that accounts for the effects of solid concentration is implemented, and the heat flux boundary condition matches the experimental setup. Radiative heat transfer is estimated using a widely used correlation during the post-processing step to determine an overall heat transfer coefficient. The model is validated against testing data and achieves less than 30% discrepancy and a heat transfer coefficient greater than 1000 W/m2K.

CSP↗
NASA NTRS1987

Spectral synthesis in the ultraviolet. I - Far-ultraviolet stellar library

A library of mean stellar energy distributions for use in population synthesis of the ultraviolet spectra of active star-forming galaxies is derived from the ultraviolet spectrophotometry in the IUE Spectral Atlas. The spectra extend from 1230 to 1930 A with a resolution of 6 A. The library contains eight main-sequence groups from O3 to A7 V, four giant groups from O5 to B9 III, and three supergiant groups from O9 to A0 I. Several continuum and spectral line indices are computed, and their usefulness as temperature and luminosity discriminants is discussed. The advantages of analyzing composite stellar populations in the far-ultraviolet, where cool starlight is strongly suppressed, are examined.

Fanelli, Michael N.↗
NASA NTRS1994

A coupled/uncoupled deformation and fatigue damage algorithm utilizing the finite element method

A fatigue damage computational algorithm utilizing a multiaxial, isothermal, continuum based fatigue damage model for unidirectional metal matrix composites has been implemented into the commercial finite element code MARC using MARC user subroutines. Damage is introduced into the finite element solution through the concept of effective stress which fully couples the fatigue damage calculations with the finite element deformation solution. An axisymmetric stress analysis was performed on a circumferentially reinforced ring, wherein both the matrix cladding and the composite core were assumed to behave elastic-perfectly plastic. The composite core behavior was represented using Hill's anisotropic continuum based plasticity model, and similarly, the matrix cladding was represented by an isotropic plasticity model. Results are presented in the form of S-N curves and damage distribution plots.

Wilt, Thomas E.↗
NASA NTRS2001

Ranges of applicability for the continuum beam model in the mechanics of carbon nanotubes and nanorods

Limitations in the validity of the continuum beam model for carbon nanotubes (NTs) and nanorods are examined. Applicability of all assumptions used in the model is restricted by the two criteria for geometric parameters that characterize the structure of NTs. The key non-dimensional parameters that control the NT buckling behavior are derived via dimensional analysis of the nanomechanical problem. A mechanical law of geometric similitude for NT buckling is extended from continuum mechanics for different molecular structures. A model applicability map, where two classes of beam-like NTs are identified, is constructed for distinct ranges of non-dimensional parameters. Expressions for the critical buckling loads and strains are tailored for two classes of NTs and compared with the data provided by the molecular dynamics simulations. copyright 2001 Elsevier Science Ltd. All rights reserved.

Continuum beam model↗
NASA NTRS2008

Multiscale Modeling of Damage Processes in fcc Aluminum: From Atoms to Grains

Molecular dynamics (MD) methods are opening new opportunities for simulating the fundamental processes of material behavior at the atomistic level. However, current analysis is limited to small domains and increasing the size of the MD domain quickly presents intractable computational demands. A preferred approach to surmount this computational limitation has been to combine continuum mechanics-based modeling procedures, such as the finite element method (FEM), with MD analyses thereby reducing the region of atomic scale refinement. Such multiscale modeling strategies can be divided into two broad classifications: concurrent multiscale methods that directly incorporate an atomistic domain within a continuum domain and sequential multiscale methods that extract an averaged response from the atomistic simulation for later use as a constitutive model in a continuum analysis.

Glaessgen, E. H.↗
NASA NTRS1987

Inverse Compton scattering of ambient radiation by a cold relativistic jet - A source of beamed, polarized continuum in blazars?

A general formalism is developed for computing the intensity and polarization of unpolarized radiation scattered by a group of electrons with a specified distribution of momenta. The case of isotropic incident radiation with a power-law spectrum is addressed, and the 'head-on' approximation is developed to describe the nearly unidirectional character of the incident radiation in the electron rest frame. The formalism is implemented numerically, verifying the validity of the 'head-on' approximation and illustrating the polarization and intensity properties of radiation scattered by jets with different electron momentum distribution. The results are placed in an astrophysical setting, determining the optical depths and kinetic energies required for the scattered radiation to dominate over the isotropic component in blazars. It is argued that Comptonization of unbeamed AGN radiation by a relativistic jet can account for the polarization and the overall spectral shape of the IR-optical continuum in blazars.

Begelman, Mitchell C.↗
NASA NTRS1973

Path-length distributions of photons diffusely reflected from a semi-infinite atmosphere.

Knowledge of the probability distribution of photon path lengths for diffuse reflection in the continuum from a planetary atmosphere greatly simplifies computation of the shape and strength of absorption lines formed under conditions of multiple scattering. Graphs of this probability distribution are presented for diffuse reflection from a semiinfinite conservative atmosphere for both isotropic and forward-directed scattering, and for several angles of incidence and reflection. Graphs showing distributions for the case of absorption in the continuum are also discussed. The distributions are obtained from an analytic expression for the path-length distribution for nth-order scattering and a tabulation of nth-order intensities. It is found that te shape of the probability distribution depends markedly on the degree of forward scattering due to the phase function. The results are compared with those obtained by a Monte Carlo computation of the sort which might be used if tabulated nth-order intensities were not available.

Appleby, J. F.↗