Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuum scale modeling”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS1977

Diffuse gamma radiation

An examination of the intensity, energy spectrum, and spatial distribution of the diffuse gamma-radiation observed by SAS-2 satellite away from the galactic plane in the energy range above 35 MeV has shown that it consists of two components. One component is generally correlated with galactic latitudes, the atomic hydrogen column density was deduced from 21 cm measurements, and the continuum radio emission, believed to be synchrotron emission. It has an energy spectrum similar to that in the plane and joins smoothly to the intense radiation from the plane. It is therefore presumed to be of galactic origin. The other component is apparently isotropic, at least on a coarse scale, and has a steep energy spectrum. No evidence is found for a cosmic ray halo surrounding the galaxy in the shape of a sphere or oblate spheroid with galactic dimensions. Constraints for a halo model with significantly larger dimensions are set on the basis of an upper limit to the gamma-ray anisotropy.

Fichtel, C. E.↗
NASA NTRS1988

Multifrequency spectra of BL Lac objects

A program to obtain simultaneous multifrequency spectra of BL Lacertae objects that are known X-ray sources is discussed. The IUE spectra are generally featureless and well-fitted by power law models. For the faintest exposures, Gaussian extraction of the spectrum can greatly impove the signal-to-noise. Most program objects vary in the ultraviolet, although the time scales are not known because of limited observing time. The broadband spectra of BL Lacs exhibit a range of characteristics but the curvature is always downward and the shape is generally smooth. This can be interpreted as synchrotron emission from a relativistic jet; different jet models are possible, and each allows a range of values for the bulk velocity, magnetic field strength, and electron density. Synchrotron models are not required, however an accretion disk model also gives a good fit to the ultraviolet-through-X-ray continuum.

Urry, C. M.↗
DOE OSTI2022

GLUE Code: A framework handling communication and interfaces between scales

Many scientific applications are inherently multiscale in nature. Such complex physical phenomena often require simultaneous execution and coordination of simulations spanning multiple time and length scales. This is possible by combining expensive small-scale simulations (such as molecular dynamics simulations) with larger scale simulations (such continuum limit/hydro solvers) to allow for considerably larger systems using task and data parallelism. However, the granularity of the tasks can be very large and often leads to load imbalance. Traditionally, we use approximations to streamline the computation of the more costly interactions and this introduces trade-offs between simulation cost and accuracy. In recent years, the available computational power and the advances in machine learning have made computing these scale-bridging interactions and multiscale simulations more feasible. One driving application has been plasma modeling in inertial confinement fusion (ICF), which is fundamentally multiscale in nature. This requires deep understanding of how to extrapolate microscopic information into macroscopically relevant scales. For example, in ICF one needs an accurate understanding of the connection between experimental observables and the underlying microphysics. The properties of the larger scales are often affected by the microscale behavior incorporated usually into the equations of state and ionic and electronic transport coefficients (Liboff, 1959; Rinderknecht et al., 2014; Rosenberg et al., 2015; Ross et al., 2017). Instead of incorporating this information using reliable molecular dynamics (MD) simulations, one often needs to use theoretical models, due to the inability of MD to reach engineering scales (Glosli et al., 2007; Marinak et al., 1998). One approach to resolve this issue is by coupling two MD simulations of different scales via force interpolation, e.g., the AdResS method (Krekeler et al., 2018; Nagarajan et al., 2013). Another approach, which we will pursue in the scope of this work, is by enabling scale bridging between MD simulations and meso/macro-scale models through the development and support of application programming interfaces that these different applications can interact through.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Micro-Mechanically Guided High-Throughput Alloy Design Exploration Towards Metastability-Induced H Embrittlement Resistance

We develop a high-throughput approach for studying H embrittlement (HE)-resistance in alloys, which is based on combinatorial compositional screening of metastability effects by in situ scanning electron microscopy H-analyses. The project objective included: (i) technique development of high-throughput screening (HTS) for HE-resistance, (ii) discovery of new metallic materials with superior HE-resistance, (iii) Multiscale verification of HE-resistance of the new alloys and H-barrier layers, from atomic scale to an engineering scale. The investigation focused on a model system that enables exploration of different metastable states in Fe-based complex-concentrated alloys. Composition spread islands with hundreds of varying compositions are fabricated on a single substrate using a combinatorial co-sputtering technique. To screen these alloys, we developed and employed a home-built SEM-based integrated analysis system capable of characterizing H-permeability, H-trapping, H-influence on mechanical properties, and other material properties, as well as atomistic to continuum simulations to study the underlying physics.

08 HYDROGEN↗
NASA NTRS1992

X-ray observations of sunspot penumbral structure

High-resolution coronal observations with the Normal Incidence X-ray Telescope (NIXT) reveal previously unobserved structure in the magnetic fields above a sunspot. The X-ray images are precisely aligned with a continuum photospheric image taken at the same time. The X-ray brightness traces magnetic field lines and shows: (1) that none of the bright loops originate in the spot umbra and (2) that some field lines emerging from the inner penumbra connect to regions far away from the spot. Such large-scale structures must remain distinct from the shallowly inclined fields in the outer penumbra. In particular, they cannot be involved in any interchange between the bright and dark filaments. This imposes constraints on models of penumbral convection.

Sams, B. J., III↗
NASA NTRS2010

Study of Plume Impingement Effects in the Lunar Lander Environment

Plume impingement effects from the descent and ascent engine firings of the Lunar Lander were analyzed in support of the Lunar Architecture Team under the Constellation Program. The descent stage analysis was performed to obtain shear and pressure forces on the lunar surface as well as velocity and density profiles in the flow field in an effort to understand lunar soil erosion and ejected soil impact damage which was analyzed as part of a separate study. A CFD/DSMC decoupled methodology was used with the Bird continuum breakdown parameter to distinguish the continuum flow from the rarefied flow. The ascent stage analysis was performed to ascertain the forces and moments acting on the Lunar Lander Ascent Module due to the firing of the main engine on take-off. The Reacting and Multiphase Program (RAMP) method of characteristics (MOC) code was used to model the continuum region of the nozzle plume, and the Direct Simulation Monte Carlo (DSMC) Analysis Code (DAC) was used to model the impingement results in the rarefied region. The ascent module (AM) was analyzed for various pitch and yaw rotations and for various heights in relation to the descent module (DM). For the ascent stage analysis, the plume inflow boundary was located near the nozzle exit plane in a region where the flow number density was large enough to make the DSMC solution computationally expensive. Therefore, a scaling coefficient was used to make the DSMC solution more computationally manageable. An analysis of the effectiveness of this scaling technique was performed by investigating various scaling parameters for a single height and rotation of the AM. Because the inflow boundary was near the nozzle exit plane, another analysis was performed investigating three different inflow contours to determine the effects of the flow expansion around the nozzle lip on the final plume impingement results.

Marichalar, Jeremiah↗
NASA NTRS2017

Retrieval of Water Vapor Column Abundance and Aerosol Properties from ChemCam Passive Sky Spectroscopy

We derive water vapor column abundances and aerosol properties from Mars Science Laboratory (MSL) ChemCam passive mode observations of scattered sky light. This paper covers the methodology and initial results for water vapor and also provides preliminary results for aerosols. The data set presented here includes the results of 113 observations spanning from Mars Year 31 L(s) = 291 deg (March 30, 2013) to Mars Year 33 L(s) = 127 deg (March 24, 2016). Each ChemCam passive sky observation acquires spectra at two different elevation angles. We fit these spectra with a discrete-ordinates multiple scattering radiative transfer model, using the correlated-k approximation for gas absorption bands. The retrieval proceeds by first fitting the continuum of the ratio of the two elevation angles to solve for aerosol properties, and then fitting the continuum-removed ratio to solve for gas abundances. The final step of the retrieval makes use of the observed CO 2 absorptions and the known CO 2 abundance to correct the retrieved water vapor abundance for the effects of the vertical distribution of scattering aerosols and to derive an aerosol scale height parameter. Our water vapor results give water vapor column abundance with a precision of +/-0.6 precipitable microns and systematic errors no larger than +/-0.3 precipitable microns, assuming uniform vertical mixing. The ChemCam-retrieved water abundances show, with only a few exceptions, the same seasonal behavior and the same timing of seasonal minima and maxima as the TES, CRISM, and REMS-H data sets that we compare them to. However ChemCam-retrieved water abundances are generally lower than zonal and regional scale from-orbit water vapor data, while at the same time being significantly larger than pre-dawn REMS-H abundances. Pending further analysis of REMS-H volume mixing ratio uncertainties, the differences between ChemCam and REMS-H pre-dawn mixing ratios appear to be much too large to be explained by large scale circulations and thus they tend to support the hypothesis of substantial diurnal interactions of water vapor with the surface. Our preliminary aerosol results, meanwhile, show the expected seasonal pattern in dust particle size but also indicate a surprising interannual increase in water-ice cloud opacities.

Mcconnochie, Timothy H.↗
NASA NTRS2001

Motions and Initial Conditions in Star-Forming Dense Cores

Under this grant in the past year we have pursued spectral-line observations of star-forming regions over size scales from 0.01 pc to 0.5 pc. Our main goal has been to measure the systematic and turbulent motions of condensing and collapsing gas. In this area, our results include (1) in 67 starless dense cores, some 19 show clear evidence of spatially extended inward motions, with typical line-of-sight inward speed 0.05-0.09 km s(sup -1) and with typical plane-of-the-sky extent 0.1-0.3 pc, (2) In some 40 nearby regions with embedded groups and clusters, we see extended infall asymmetry in lines of CS and HCO(+) clearly in 4 regions and less clearly in 4 others, (3) Using finer resolution (15 arcsec or 0.01-0.02 pc) and lines tracing higher density, we see spatial concentration of infall asymmetry near the protostars in NGC 1333 IRS 4A and B, L483, and L1251B, and with still finer resolution (2 arcsec or 0.003 pc or 600 AU) we detect inverse P Cyg profiles, indicating absorption of continuum emission from the protostellar envelope by infalling gas in NGC 1333 IRS 4A and 4B. Further, at high resolution we identify regions of stellar mass and low turbulence ("kernels") which are good candidates to become the next generation of stars in embedded clusters. In addition we have completed a survey for the OH Zeeman effect in absorption against nearby H II regions, indicating that the large-scale magnetic field may be nearly critical if it typically threads a flattened structure. We have also developed a model of spatially extended infall motions based on dissipation of turbulence in a magnetized, selfgravitating layer. In the following we describe some of these results in more detail.

Myers, Philip C.↗
DOE OSTI2023

A subdivision-stabilized B-spline mixed material point method

Subjected to external loadings, polymeric materials, e.g., biological tissues, hydrogels, and elastomers, may undergo extreme, nearly incompressible, (self-)contact deformations. For numerical modeling employing mesh-based techniques such as the finite element method (FEM), these deformations pose significant challenges due to large distortions in the deformed geometry, accuracy issues stemming from volumetric locking effects, and increased computational cost from complex contact searches. As an alternative to mesh-based methods, the material point method (MPM), a continuum-based particle technique, is gaining attention for its ability to handle extreme distortions and capture no-slip contact without added cost. For nearly incompressible material behaviors, while mixed formulations can address locking effects by treating displacements and pressure as independent fields, they can suffer from numerical instabilities close to the incompressibility limit due to the violation of the inf-sup condition, leading to inaccurate nodal pressure solutions. Here we propose an efficient and stable mixed B-spline material point method with highest achievable regularity for quasi-compressible polymeric materials. Using the two-scale relation of B-splines, we introduce a subdivision-stabilization for the two-field mixed MPM and obtain numerically stable, oscillation-free nodal solutions with equal-order interpolations with optimal regularity. Building on the Eulerian-Lagrangian nature of MPM, a previously-converged solution framework is adopted to mitigate issues related to cell-crossing and numerical fracture artifact present in standard MPM. We assess the stability and accuracy of the developed mixed MPM at large deformations for soft materials through the benchmark Cook’s membrane problem. Additionally, we test the robustness of the proposed MPM by modeling several examples, including the compression and indentation of a circular block into a quasi-compressible substrate and the twisting deformation of a rectangular block. The findings demonstrate the MPM’s capabilities for modeling practical soft material applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Chemo‐Mechanical Coupling in Hydrogels: Dynamics in the Diffusion‐Limited Regime

Hydrogels are characterized by substantial volume changes in response to external stimuli, making them promising candidates for developing smart materials with enhanced adaptability and responsiveness. By integrating chemical reactions, hydrogels acquire dynamic and tunable responsiveness to external stimuli through chemo‐mechanical coupling, expanding their potential in emerging applications. However, capturing their transient behavior remains challenging due to the complex interplay of chemical reactions, solvent transport, and polymer network deformation. Classical theories capture equilibrium swelling but fail to describe time‐dependent phenomena. To address this, a time‐dependent continuum model is developed that explicitly couples these processes. Volume phase transition in hydrogels with homogeneous chemical reactions is investigated, then the effects of reaction kinetics on these transitions are analyzed. The coupling of these mechanisms is further explored through a study of transient mechanical instabilities. To illustrate the impact of distinct reaction and diffusion timescales, a photo‐active gripper is studied for robotic applications. Finally, a photo‐active microswimmer that exhibits non‐reciprocal motion is proposed to highlight the solvent diffusion for locomotion at the micro‐ and nano‐ scales. The work establishes that transient dynamics of chemo‐mechanical hydrogels generate functions not accessible by steady state models and provides a predictive platform for designing adaptive materials in emerging applications.

chemo-mechanical coupling↗
NASA NTRS1987

Explosive evaporation in solar flares

This paper develops a simple analytical model for the phenomenon of 'explosive evaporation' driven by nonthermal electron heating in solar flares. The model relates the electron energy flux and spectrum, plus details of the preflare atmosphere, to the time scale for explosive evaporation to occur, the maximum pressure and temperature to be reached, rough estimates for the UV pulse emission flux and duration, and the evolution of the blueshifted component of the soft X-ray lines. An expression is given for the time scale for buildup to maximum pressures and the onset of rapid motion of the explosively evaporating plasma. This evaporation can excite a rapid response of UV line and continuum emission. The emission lines formed in the plasma approach a given emissivity-weighted blueshift speed.

Fisher, George H.↗
NASA NTRS1993

Asteroid collisions: Target size effects and resultant velocity distributions

To study the dynamic fragmentation of rock to simulate asteroid collisions, we use a 2-D, continuum damage numerical hydrocode which models two-body impacts. This hydrocode monitors stress wave propagation and interaction within the target body, and includes a physical model for the formation and growth of cracks in rock. With this algorithm we have successfully reproduced fragment size distributions and mean ejecta speeds from laboratory impact experiments using basalt, and weak and strong mortar as target materials. Using the hydrocode, we have determined that the energy needed to fracture a body has a much stronger dependence on target size than predicted from most scaling theories. In addition, velocity distributions obtained indicate that mean ejecta speeds resulting from large-body collisions do not exceed escape velocities.

Ryan, Eileen V.↗
NASA NTRS1993

Plasma-screening effects on the electron-impact excitation of hydrogenic ions in dense plasmas

Plasma-screening effects are investigated on electron-impact excitation of hydrogenic ions in dense plasmas. Scaled cross sections Z(exp 4) sigma for 1s yields 2s and 1s yields 2p are obtained for a Debye-Hueckel model of the screened Coulomb interaction. Ground and excited bound wave functions are modified in the screened Coulomb potential (Debye-Hueckel model) using the Ritz variation method. The resulting atomic wave functions and their eigenenergies agree well with the numerical and high-order perturbation theory calculations for the interesting domain of the Debye length not less than 10. The Born approximation is used to describe the continuum states of the projectile electron. Plasma screening effects on the atomic electrons cannot be neglected in the high-density cases. Including these effects, the cross sections are appreciably increased for 1s yields 2s transitions and decreased for 1s yields 2p transitions.

Jung, Young-Dae↗
NASA NTRS1984

Voyager 1 imaging and IRIS observations of Jovian methane absorption and thermal emission: Implications for cloud structure

Images from three filters of the Voyager 1 wide angle camera are used to measure the continuum reflectivity and spectral gradient near 6000 A and the 6190 A band methane/continuum ratio for a variety of cloud features in Jupiter's atmosphere. The dark barge features in the North Equatorial Belt have anomalously strong positive continuum spectral gradients suggesting unique composition. Methane absorption is shown at unprecedented spatial scales for the Great Red Spot and its immediate environment, for a dark barge feature in the North Equatorial Belt, and for two hot spot and plume regions in the North Equatorial Belt. Methane absorption and five micrometer emission are correlated in the vicinity of the Great Red Spot but are anticorrelated in one of the plume hot spot regions. Methane absorption and simultaneous maps of five micrometer brightness temperature is quantitatively compared to realistic cloud structure models which include multiple scattering at five micrometer as well as in the visible. Variability in H2 quadrupole lines are also investigated.

West, R. A.↗
NASA NTRS1985

Voyager 1 imaging and IRIS observations of Jovian methane absorption and thermal emission - Implications for cloud structure

Images from three filters of the Voyager 1 wide angle camera are used to measure the continuum reflectivity and spectral gradient near 6000 A and the 6190 A band methane/continuum ratio for a variety of cloud features in Jupiter's atmosphere. The dark barge features in the North Equatorial Belt have anomalously strong positive continuum spectral gradients suggesting unique composition. Methane absorption is shown at unprecedented spatial scales for the Great Red Spot and its immediate environment, for a dark barge feature in the North Equatorial Belt, and for two hot spot and plume regions in the North Equatorial Belt. Methane absorption and five micrometer emission are correlated in the vicinity of the Great Red Spot but are anticorrelated in one of the plume hot spot regions. Methane absorption and simultaneous maps of five micrometer brightness temperature are quantitatively compared to realistic cloud structure models which include multiple scattering at five micrometer as well as in the visible. Variability in H2 quadrupole lines are also investigated.

West, R. A.↗
NASA NTRS1993

Boundary formulations for sensitivity analysis without matrix derivatives

A new hybrid approach to continuum structural shape sensitivity analysis employing boundary element analysis (BEA) is presented. The approach uses iterative reanalysis to obviate the need to factor perturbed matrices in the determination of surface displacement and traction sensitivities via a univariate perturbation/finite difference (UPFD) step. The UPFD approach makes it possible to immediately reuse existing subroutines for computation of BEA matrix coefficients in the design sensitivity analysis process. The reanalysis technique computes economical response of univariately perturbed models without factoring perturbed matrices. The approach provides substantial computational economy without the burden of a large-scale reprogramming effort.

Kane, J. H.↗
NASA NTRS2004

A Hierarchical Approach to Fracture Mechanics

Recent research conducted under NASA LaRC's Creativity and Innovation Program has led to the development of an initial approach for a hierarchical fracture mechanics. This methodology unites failure mechanisms occurring at different length scales and provides a framework for a physics-based theory of fracture. At the nanoscale, parametric molecular dynamic simulations are used to compute the energy associated with atomic level failure mechanisms. This information is used in a mesoscale percolation model of defect coalescence to obtain statistics of fracture paths and energies through Monte Carlo simulations. The mathematical structure of predicted crack paths is described using concepts of fractal geometry. The non-integer fractal dimension relates geometric and energy measures between meso- and macroscales. For illustration, a fractal-based continuum strain energy release rate is derived for inter- and transgranular fracture in polycrystalline metals.

Saether, Erik↗
NASA NTRS1993

BBXRT and Ginga observations of the Seyfert 1 galaxy Markarian 335

We present new X-ray observations of Mrk 335 in the 0.4-10 keV band using BBXRT and in the 2-20 keV band with Ginga. The BBXRT data show a break to a softer spectrum below about 1 keV. The Ginga data show evidence for variability in absorption, uncorrelated with flux changes of the continuum source. The inconsistency between the opacity observed at the iron K edge and at softer X-ray energies implies that the absorber has an ionization state about FeXV. We suggest that the nuclear emission can be explained by a self-consistent pair cascade model as developed by Zdziarski et al. (1990) and Zdziarski and Coppi (1991), where the compactness parameter inferred from the spectral fits to both data sets is in good agreement with the value calculated from the ratio of observed X-ray luminosity and variability time scale.

Turner, T. J.↗