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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Predicting Fluid Flow Regime, Permeability, and Diffusivity in Mudrocks from Multiscale Pore Characterisation

In geoenergy applications, mudrocks prevent fluids to leak from temporary (H 2 , CH 4 ) or permanent (CO 2 , radioactive waste) storage/disposal sites and serve as a source and reservoir for unconventional oil and gas. Understanding transport properties integrated with dominant fluid flow mechanisms in mudrocks is essential to better predict the performance of mudrocks within these applications. In this study, small-angle neutron scattering (SANS) experiments were conducted on 71 samples from 13 different sets of mudrocks across the globe to capture the pore structure of nearly the full pore size spectrum (2 nm–5 μm). We develop fractal models to predict transport properties (permeability and diffusivity) based on the SANS-derived pore size distributions. The results indicate that transport phenomena in mudrocks are intrinsically pore size-dependent. Depending on hydrostatic pore pressures, transition flow develops in micropores, slip flow in meso- and macropores, and continuum flow in larger macropores. Fluid flow regimes progress towards larger pore sizes during reservoir depletion or smaller pore sizes during fluid storage, so when pressure is decreased or increased, respectively. Capturing the heterogeneity of mudrocks by considering fractal dimension and tortuosity fractal dimension for defined pore size ranges, fractal models integrate apparent permeability with slip flow, Darcy permeability with continuum flow, and gas diffusivity with diffusion flow in the matrix. This new model of pore size-dependent transport and integrated transport properties using fractal models yields a systematic approach that can also inform multiscale multi-physics models to better understand fluid flow and transport phenomena in mudrocks on the reservoir and basin scale.

36 MATERIALS SCIENCE↗
DOE OSTI2023

ProvLight: Efficient Workflow Provenance Capture on the Edge-to-Cloud Continuum

Modern scientific workflows require hybrid infrastructures combining numerous decentralized resources on the IoT/Edge interconnected to Cloud/HPC systems (aka the Computing Continuum) to enable their optimized execution. Understanding and optimizing the performance of such complex Edge-to-Cloud workflows is challenging. Capturing the provenance of key performance indicators, with their related data and processes, may assist in understanding and optimizing workflow executions. However, the capture overhead can be prohibitive, particularly in resource-constrained devices, such as the ones on the IoT/Edge.To address this challenge, based on a performance analysis of existing systems, we propose ProvLight, a tool to enable efficient provenance capture on the IoT/Edge. We leverage simplified data models, data compression and grouping, and lightweight transmission protocols to reduce overheads. We further integrate ProvLight into the E2Clab framework to enable workflow provenance capture across the Edge-to-Cloud Continuum. This integration makes E2Clab a promising platform for the performance optimization of applications through reproducible experiments.We validate ProvLight at a large scale with synthetic workloads on 64 real-life IoT/Edge devices in the FIT IoT LAB testbed. Evaluations show that ProvLight outperforms state-of-the-art systems like ProvLake and DfAnalyzer in resource-constrained devices. ProvLight is 26—37x faster to capture and transmit provenance data; uses 5—7x less CPU; 2x less memory; transmits 2x less data; and consumes 2—2.5x less energy. ProvLight [1] and E2Clab [2] are available as open-source tools.

Rosendo, Daniel↗
NASA NTRS1976

Solar limb brightening in submillimeter wavelengths

Differential two-beam scans of the sun in submillimeter wavelengths (350 microns to 1 millimeter) indicate limb brightening approaching 1 percent when the cosine of the angle from the normal equals 0.60. The observations also show considerable chromospheric structure, both in active and quiet regions, but with less relative amplitude than at millimeter and centimeter wavelengths. The limited angular resolution of the observing system, together with photometric errors due to fluctuating atmospheric transparency, make the brightness profile of the extreme limb uncertain. The observed degree of limb brightening is considerably less than that consistent with spherically symmetric model atmospheres based on continuum brightness-temperature measurements. The suppression of limb brightening suggests the existence of irregular granular structure with both horizontal and vertical characteristic sizes of the order of 1500 km. High-resolution images in the wings of the K-line show granular structure of about this horizontal scale.

Lindsey, C.↗
NASA NTRS1980

Discovery of intrinsic linear polarization in SS433

Observations of variable linear polarization in the optical range in the unique radio and X-ray source SS 433, which is characterized by moving emission lines in its visible and near-IR spectra, are reported. The observations were made with an unfiltered GaAs photocathode covering the spectral region 3000-9000 A, with a passband characteristic wavelength near 7000 A. Measurements reveal the presence of both a long-term and a very short time-scale variation in linear polarization in the continuum and in H-alpha, which indicate that a substantial fraction of the polarization is intrinsic to SS 433 and not produced in the interstellar medium. It is suggested that the most likely source of polarization in the optical continuum is electron scattering, and constraints on models of SS 433 implied by the discovery of intrinsic polarization are considered.

Mclean, I. S.↗
NASA NTRS1984

Detection of the J-10 Manifold of the Pure Rotational Band of Phosphine on Saturn

The detection of the J = 10 manifold of the pure rotational band of PH3 on Saturn is reported. The observations were made from the far-infrared cooled grating spectrometer. The wavelengths and observed brightness temperatures for the full disk plus rings are 89 + or - 3 K at 97.04 micrometer, 77 + or - 3 K at 102.72 micrometer, 77 + or - 3 K at 102.94 micrometer, and 83 + or - 3 K at 105.12 micrometers. The points of 97.04 and 105.12 micrometers establish the continuum level and the two points near 103 micrometers measure the depth of the PH3 manifold. After the flux due to the rings is subtracted, the depth of the feature is 16 + or - 6 K relative to the nearby 102 K continuum. These results are compared to theoretical models which parameterize the PH3 mixing ratio as x = x sub zero (P/P sub zero)(alpha) for P P sub zero and as x = x sub zero for P or = P, where P is the total pressure and alpha = H/h is the ratio of the dynamical scale height (H) and the scale height for decreasing the PH3 mixing ratio (h). The parameters x sub zero, P sub zero, and h were varied, as well as the H/He mixing ratio and the pressure-temperature profile. The data are well fitted using pressure-temperature profiles. The preferred values of h, P sub zero, and x sub zero imply that there is little or no PH3 above the thermal inversion and that the mixing ratio below the inversion is consistent with PH3 being 1 to 4 times overabundant relative to the solar P/H ratio.

Haas, M. R.↗
NASA NTRS1985

On the accuracy of modelling the dynamics of large space structures

Proposed space missions will require large scale, light weight, space based structural systems. Large space structure technology (LSST) systems will have to accommodate (among others): ocean data systems; electronic mail systems; large multibeam antenna systems; and, space based solar power systems. The structures are to be delivered into orbit by the space shuttle. Because of their inherent size, modelling techniques and scaling algorithms must be developed so that system performance can be predicted accurately prior to launch and assembly. When the size and weight-to-area ratio of proposed LSST systems dictate that the entire system be considered flexible, there are two basic modeling methods which can be used. The first is a continuum approach, a mathematical formulation for predicting the motion of a general orbiting flexible body, in which elastic deformations are considered small compared with characteristic body dimensions. This approach is based on an a priori knowledge of the frequencies and shape functions of all modes included within the system model. Alternatively, finite element techniques can be used to model the entire structure as a system of lumped masses connected by a series of (restoring) springs and possibly dampers. In addition, a computational algorithm was developed to evaluate the coefficients of the various coupling terms in the equations of motion as applied to the finite element model of the Hoop/Column.

Diarra, C. M.↗
NASA NTRS1985

On the accuracy of modelling the dynamics of large space structures

Proposed space missions will require large scale, light weight, space based structural systems. Large space structure technology (LSST) systems will have to accommodate (among others): ocean data systems; electronic mail systems; large multibeam antenna systems; and, space based solar power systems. The structures are to be delivered into orbit by the Space Shuttle. Because of their inherent size, modelling techniques and scaling algorithms must be developed so that system performance can be predicted accurately prior to launch and assembly. When the size and weight-to-area ratio of proposed LSST systems dictate that the entire system be considered flexible, there are two basic modelling methods which can be used. The first is a continuum approach, a mathematical formulation for predicting the motion of a general orbiting flexible body, in which elastic deformations are considered small compared with characteristic body dimensions. This approach is based on an a priori knowledge of the frequencies and shape functions of all modes included within the system model. Alternatively, finite element techniques can be used to model the entire structure as a system of lumped masses connected by a series of (restoring) springs and possibly dampers. In addition, a computational algorithm was developed to evaluate the coefficients of the various coupling terms in the equations of motion as applied to the finite element model of the Hoop/Column.

Diarra, C. M.↗
DOE OSTI2024

Computational Methods for Modeling Electrospray Microdroplet Chemistry for Improved Quantitative Mass Spectrometry

This project aimed at enhancing the quantitative analysis capabilities of electrospray ionization mass spectrometry (ESI-MS) by developing advanced computational methods. The primary focus was to integrate continuum and molecular dynamics simulations to study the behavior of microdroplets in the ESI process, from formation to evaporation. Through this research, we sought to bridge significant length and time scales to provide a comprehensive understanding of how analyte concentrations evolve from bulk solutions into gas-phase ions. This understanding is crucial for addressing challenges such as ionization efficiency, solvent effects, and ion suppression, which currently limit the accuracy of quantitative ESI-MS.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1987

Large-scale Galactic dust morphology and physical conditions from IRAS observations

In this paper, the zodiacal component is subtracted from the 60 and 100 micron Galactic plane emission by applying an empirical model derived from IRAS data in regions of the sky not dominated by the Galaxy. The corrected observations are used to derive the large-scale physical conditions such as temperature, optical depth, and total FIR brightness of the dust residing in the Galactic disk. (C-12)O, H I, and 5 GHz radio continuum observations are also used to compare the large-scale properties of the gas and dust distributions in the Galaxy. Possible scenarios to explain the findings are suggested.

Sodroski, T. J.↗
DOE OSTI2021

Discrete ion stochastic continuum overdamped solvent algorithm for modeling electrolytes

In this paper we develop a methodology for the mesoscale simulation of strong electrolytes. The methodology is an extension of the fluctuating immersed-boundary approach that treats a solute as discrete Lagrangian particles that interact with Eulerian hydrodynamic and electrostatic fields. In both algorithms the immersed-boundary method of Peskin is used for particle-field coupling. Hydrodynamic interactions are taken to be overdamped, with thermal noise incorporated using the fluctuating Stokes equation, including a "dry diffusion" Brownian motion to account for scales not resolved by the coarse-grained model of the solvent. Long-range electrostatic interactions are computed by solving the Poisson equation, with short-range corrections included using an immersed-boundary variant of the classical particle-particle particle-mesh technique. Also included is a short-range repulsive force based on the Weeks-Chandler-Andersen potential. This methodology is validated by comparison to Debye-Hückel theory for ion-ion pair correlation functions, and Debye-Hückel-Onsager theory for conductivity, including the Wien effect for strong electric fields. In each case, good agreement is observed, provided that hydrodynamic interactions at the typical ion-ion separation are resolved by the fluid grid.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1998

Hubble Space Telescope Eclipse Observations of the Nova Like Cataclysmic Variable UX Ursae Majoris

We present and analyze Hubble Space Telescope observations of the eclipsing nova-like cataclysmic variable UX UMa obtained with the Faint Object Spectrograph. Two eclipses each were observed with the G160L grating (covering the ultraviolet waveband) in 1994 August and with the PRISM (covering the near-ultraviolet to near-infrared) in November of the same year. The system was about 50% brighter in November than in August, which, if due to a change in the accretion rate, indicates a fairly substantial increase in Mass accretion by about 50%. The eclipse light curves are qualitatively consistent with the gradual occultation of an accretion disk with a radially decreasing temperature distribution. The light curves also exhibit asymmetries about mideclipse that are likely due to a bright spot at the disk edge. Bright-spot spectra have been constructed by differencing the mean spectra observed at pre- and posteclipse orbital phases. These difference spectra contain ultraviolet absorption lines and show the Balmer jump in emission. This suggests that part of the bright spot may be optically thin in the continuum and vertically extended enough to veil the inner disk and/or the outflow from UX UMa in some spectral lines. Model disk spectra constructed as ensembles of stellar atmospheres provide poor descriptions of the observed posteclipse spectra, despite the fact that UX UMa's light should be dominated by the disk at this time. Suitably scaled single temperature model stellar atmospheres with T(sub eff) approximately equals 12,500-14,500 K actually provide a better match to both the ultraviolet and optical posteclipse spectra. Evidently, great care must be taken in attempts to derive accretion rates from comparisons of disk models to observations. One way to reconcile disk models with the observed posteclipse spectra is to postulate the presence of a significant amount of optically thin material in the system. Such an optically thin component might be associated with the transition region ("chromosphere") between the disk photosphere and the fast wind from the system whose presence has been suggested by Knigge and Drew. In any event, the wind/ chromosphere is likely to be the region in which many, if not most, of the UV lines are formed. This is clear from the plethora of emission lines that appear in the mideclipse spectra, some of which appear as absorption features in spectra taken at out-of-eclipse orbital phases.

Knigge, Christian↗
NASA NTRS1994

(abstract) The Circumstellar Environment of the Extremely Young Protostellar Source L1448IRS3

The class 0 sources form an interesting new category of protostellar objects. Many have strong millimeter continuum emission and exhibit jetlike outflows. There are suggestions that these objects are systematically younger than typical embedded (class I) sources. We are investigating the properties of class 0 sources to determine whether they are indeed very young or perhaps represent extreme physical conditions, such as rapid rotation. We present millimeter interferometric data for the class 0 object known as L1448 IRS3. This young low-mass star has extremely strong millimeter continuum emission. The interferometer data show the emission is resolved on a scale of a few arcseconds (1000 AU) This suggests the bulk of the dust continuum emission originates in an 'infall' envelope rather than a protostar disk. In addition, the C180 line data display a strong velocity gradiant which indicates the dense core is rapidly rotating. We compare our data with the predictions of protostellar collapse models.

protostellar objects young protostellar sources mi↗
NASA NTRS1994

HST ultraviolet spectropolarimetry of NGC 1068

We present multiaperture spectropolarimetric observations of the Seyfert 2 galaxy NGC 1068. The nuclear region of NGC 1068 was observed with three apertures (4.3 sec x 1.4 sec., 1.0 sec., and 0.3 sec.) by the Hubble Space Telescope (HST). These ultraviolet data allow us to study the scattered, nuclear light without the complication of polarization dilution caused by strong unpolarized starlight. In particular, we can determine the mechanism of the scattering and the size scale of the scattering region. Our data confirm the findings of Antonicci & Miller (1985) and Code et al. (1993). Shortward of 2700 A, the continuum polarization is constant, P is approximately 16%, at a position angle of 97 deg, indicating electron scattering as the mechanism for the polarization. In addition, the narrow-line emission is much less polarized than the continuum, and two broad lines, Mg II lambda lambda 2796, 2804 and C III lambda 1909, are seen in polarized flux Ly-alpha and C IV are not included in our spectral region). The total flux data also indicate that from approximately 2200 A to approximately 3000 A, the spectrum of NGC 1068 is strongly affected by blended Fe II emission. These polarimetry data are explained by the occulation/reflection model proposed by Antonucci & Miller (1985). According to this model, NGC 1068 harbors a Seyfert 1 nucleus which is obscured along the line of sight by an opaque torus. The symmetry axis of the torus is alinged with the radio axis. Radiation from the broad-line region and feature continuum source can escape along the poles of the torus where it is scattered to us. This scattered light receives a partial linear polarization, perpendicular to the symmetry axis of the torus. Our measurements of the relative fluxes and polarized fluxes in the three apertures indicated that the scattering region is spatially extended on a scale of approximatley 1 sec. Also, in the 0.3 sec aperture, and to a lesser extent in the 1.0 sec aperture, the polarization is is somewhat higher, and the position angle somewhat lower, than in largest aperture. Presumably, isolating the central part of the relecting region results in less geometrical cancellation of the polarization and an angle indicative of the location of the innermost scattering paricles.

Antonucci, Robert↗
DOE OSTI2023

Near-continuum, hypersonic oxygen flow over a double cone simulated by direct simulation Monte Carlo informed from quantum chemistry

A large-scale, fully resolved direct simulation Monte Carlo (DSMC) computation of a non-equilibrium, reactive flow of pure oxygen over a double cone is presented. Under the simulated near-continuum conditions, the computational demands are shown to be significant because of the wide range of length scales that must be resolved. Therefore, robust grid adaption capabilities and efficient parallelization of the Stochastic PArallel Rarefied-gas Time-accurate Analyzer (SPARTA) code that is utilized in this work are essential. The thermochemical and transport collision models were selected for efficiency and simplicity. First-principles data, obtained from the highly accurate direct molecular simulation method, were used to inform the collision models’ parameters. Importantly, because SPARTA implements molecular collision models using collision-specific energies, the resulting macroscopic relaxation rates were evaluated a posteriori via zero-dimensional heat bath simulations. The comparisons of surface properties, namely heat flux and pressure, show very close agreement with previous computational fluid dynamics (CFD) results. Differences with the measurements were found to be similar to the CFD simulations. The unresolved discrepancy with the measurements could be due to inconsistent free stream conditions with the actual experimental data or missing physical phenomena altogether, for example atomic and molecular oxygen electronically excited states, three-dimensional effects, or more complex gas–surface interactions. As shown in this work, the advantages of obtaining a DSMC particle solution for these flows reside in the method's ability to be directly informed from first principles and to seamlessly describe internal energy non-equilibrium for all modes. With the advent of exascale computing and beyond, particle methods will be an increasingly important tool to verify the validity of physical assumptions in reduced-order models via fully resolved, experimental-scale simulations, down to the level of molecular-level distributions.

Mechanics↗
DOE OSTI2024

Fracture Dynamics in Silicon Anode Solid-State Batteries

Solid-state batteries (SSBs) with silicon anodes could enable improved safety and energy density compared to lithium-ion batteries. However, degradation arising from the massive volumetric changes of silicon anodes during cycling is not well understood in solid-state systems. Here, we use operando X-ray computed microtomography to reveal micro- to macro-scale chemo-mechanical degradation processes of silicon anodes in SSBs. Mud-type channel cracks driven by biaxial tensile stress form across the electrode during delithiation. We also find detrimental cracks at the silicon/ solid electrolyte interface that form due to local reaction competition between neighboring domains of different sizes. Continuum phase-field damage modeling quantifies stress-driven channel cracking and shows that the lithiated silicon stress state is critical for determining the extent of interfacial fracture. This work reveals the mechanisms that govern SSBs compared to conventional lithium-ion batteries and provides guidelines for engineering chemomechanically resilient electrodes for high-energy batteries.

25 ENERGY STORAGE↗
NASA NTRS1987

Hydrogen emission from Jupiter: Hydrogen emission from sunlit atmosphere of Saturn

Successful IUE observations of the equatorial sunlit atmosphere of Jupiter and Saturn have been obtained. Spectra containing atomic and molecular hydrogen and solar reflection continuum emissions have been analyzed, with the purpose of determining the long term temporal behavior of the electroglow process. Quantitative estimates have been established for the first time using a model analysis of the short wavelength region of the spectrum. Both systems show varying degrees of long term variability in hydrogen emission rate, but the time scale is too short to determine whether there is a dependence on solar cycle activity. As part of the emission modeling program, a preliminary point source spreading function for the IUE SWP instrument has been established, suggesting a wavelength dependence in spectral line width different from previous analyses. Further IUE observations are planned for both Jupiter and Saturn.

Shemansky, D. E.↗
NASA NTRS1986

Interpretation of the prominence differential emissions measure for 3 geometries

Researchers have used prominence extreme ultraviolet line intensities observed from Skylab to derive the differential emission measure Q(T) in the prominence-corona (PC) interface from 3 x 10,000 to 3 times 1 million K, including the effects of Lyman Continuum absorption. Using lines both shortward and longward of the Lyman limit, researchers have estimated the importance of absorption as function of temperature. The magnitude of the absorption, as well as its rate of increase as a function of temperature, place limits on the thread scales and the character of the interfilar medium. Researchers have calculated models based on three assumed geometries: (1) threads with hot sheaths and cool cores; (2) isothermal threads; and (3) threads with longitudinal temperature gradients along the magnetic field. Comparison of the absorption computed from these models with the observed absorption in prominences shows that none of the geometries is totally satisfactory.

Schmahl, E. J.↗
NASA NTRS1992

SUSIM/UARS observations of the 120 to 300 nm flux variations during the maximum of the solar cycle - Inferences for the 11-year cycle

Since October 1991, the Solar Ultraviolet Spectral Irradiance Monitor on board the Upper Atmosphere Research Satellite has been measuring the solar spectral irradiance from 120 to 400 nm with low (5 nm) and moderate (1.1 nm) resolution on a daily basis. By scaling the modulation measured over four solar rotations to proxy indicators of solar activity, we estimate the amplitude of the solar ultraviolet irradiance variation during the 11-year activity cycle, extending the wavelength coverage of existing empirical variability models to 120 nm and distinguishing the variability of the emission lines from that of the underlying solar continuum. These estimates are compared with results from another empirical variability model and from direct measurements.

Lean, Judith↗