SEARCH · Engineering Papers

Results for “continuum 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 613 records · Page 34

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Sierra/SolidMechanics 5.4 User's Guide

Sierra/SolidMechanics (Sierra/SM) is a Lagrangian, three-dimensional code for finite element analysis of solids and structures. It provides capabilities for explicit dynamic, implicit quasistatic and dynamic analyses. The explicit dynamics capabilities allow for the efficient and robust solution of models with extensive contact subjected to large, suddenly applied loads. For implicit problems, Sierra/SM uses a multi-level iterative solver, which enables it to effectively solve problems with large deformations, nonlinear material behavior, and contact. Sierra/SM has a versatile library of continuum and structural elements, and a large library of material models. The code is written for parallel computing environments enabling scalable solutions of extremely large problems for both implicit and explicit analyses. It is built on the SIERRA Framework, which facilitates coupling with other SIERRA mechanics codes. This document describes the functionality and input syntax for Sierra/SM.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2022

Sierra/SolidMechanics 5.8 User's Manual

42 ENGINEERING↗

DOE OSTI2022

Sierra/SolidMechanics 5.10 User's Guide

42 ENGINEERING↗

DOE OSTI2023

Sierra/Solid Mechanics 5.16 User's Guide

36 MATERIALS SCIENCE↗

DOE OSTI2020

Sierra/SolidMechanics 4.56 User's Guide

Siera/SolidMechanics (Sierra / SM) is a Lagrangian, three-dimensional code for finite element analysis of solids and structures. It provides capabilities for explicit dynamic, implicit quasistatic and dynamic analyses. The explicit dynamics capabilities allow for the efficient and robust solution of models with extensive contact subjected to large, suddenly applied loads. For implicit problems, Sierra / SM uses a multi-level iterative solver, which enables it to effectively solve problems with large deformations, nonlinear material behavior, and contact. Sierra / SM has a versatile library of continuum and structural elements, and a large library of material models. The code is written for parallel computing environments enabling scalable solutions of extremely large problems for both implicit and explicit analyses. It is built on the SIERRA Framework, which facilitates coupling with other SIERRA mechanics codes . This document describes the functionality and input syntax for Sierra/SM.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

A diffusion–deformation model with damage for polymer undergoing rapid decompression failure

Green hydrogen is emerging as one of the most promising alternatives to replace fossil fuels. While hydrogen gas has a good energy density by weight, its poor energy density by volume requires it to be stored under high pressure for commercial use. The hydrogen infrastructure developed to handle this high pressure hydrogen contains multiple components made of polymers as sealing agents. Although polymers do not react chemically with hydrogen gas, they undergo a mechanical failure when high pressure hydrogen gas is suddenly depressurized. This phenomenon known as 'rapid decompression failure' occurs due to the diffusion of hydrogen through polymer and getting trapped inside pre-existing cavities or voids. Here, in this paper, a continuum mechanics-based coupled diffusion-deformation-damage model was developed to predict the hydrogen distribution, stress distribution, and damage propagation inside the polymer while it undergoes rapid decompression failure. The polymer was modeled as a hyperelastic material because it represents the nonlinear material response observed in uniaxial tensile tests perfectly. The effects of hydrogen diffusivity, pre-existing cavity size, cavity location, applied hydrogen pressure, and depressurization rate on damage initiation were studied. It was found that the coefficient of diffusion plays an important role in damage initiation and damage was mostly concentrated in the inside areas rather than near the surface. Experiments were conducted with EPDM polymer which agreed well with the predicted trends using the given model. The effect of adding carbon black and silica filler particles and plasticizer to the pure EPDM polymer was also studied. It was found that damage during RDF decreases with the addition of fillers, but increases with the addition of the plasticizer. Finally, the damage evolution in the presence of two cavities was also studied, and was found that the interaction of stress fields around the cavities alters the damage occurring during RDF.

36 MATERIALS SCIENCE↗

DOE OSTI2026

From Reproducible Edge–Cloud Experimentation to Real-World Practice: The E2Clab Experience

Reproducibility is already difficult in distributed systems; on the computing continuum, it becomes substantially harder. Applications that span sensing devices, edge and fog resources, and cloud platforms must be evaluated across heterogeneous hardware, variable network conditions, cross-layer orchestration decisions, and long-running workflow lifecycles. We use E2Clab as a case study to examine these challenges and their implications for experimental methodology. We explain why reproducible experimentation is harder on the continuum, then revisit E2Clab as an initial response based on explicit modeling of infrastructure, workflow lifecycle, and artifacts. Lastly, we discuss how its evolution toward more realistic application settings can be understood through the lens of Translational Computer Science. We argue that reproducible continuum experimentation requires methods that are rigorous enough for research while remaining adaptable to real-world practice.

42 ENGINEERING↗

DOE OSTI2023

Distance-dependent dielectric constant at the calcite/electrolyte interface: Implication for surface complexation modeling

Hypothesis: The electrical double layer formed at the mineral/electrolyte interface is often modeled using mean-field approaches based on a continuum description of the solvent whose dielectric constant is assumed to decrease monotonically with decreasing distance to the surface. In contrast, molecular simulations show that the solvent polarizability oscillates near the surface similar to the water density profile - as shown previously, for example, by Bonthuis et al. (D.J. Bonthuis, S. Gekle, R.R. Netz, Dielectric Profile of Interfacial Water and its Effect on Double-Layer Capacitance, Phys Rev Lett 107(16) (2011) 166102). Here, we showed that molecular and mesoscale pictures agree by spatially averaging the dielectric constant obtained from molecular dynamics simulations over the distances relevant to the mean-field representation. In addition, the values of capacitances used to describe the electrical double layer in Surface Complexation Models (SCMs) of the mineral/electrolyte interface can be estimated using molecularly informed spatially averaged dielectric constants and positions of hydration layers. Experiments: First, we used molecular dynamics simulations to model the calcite 101¯4/electrolyte interface. Next, by using atomistic trajectories, we calculated the distance-dependent static dielectric constant and water density in the direction normal to the. Finally, we applied spatial compartmentalization consistent with the model of parallel-plate capacitors connected in series to estimate SCM capacitances. Findings: Computationally expensive simulations are required to determine the dielectric constant profile of interfacial water near the mineral surface. On the other hand, water density profiles are readily assessable from much shorter simulation trajectories. Our simulations confirmed that dielectric and water density oscillations at the interface are correlated. Here, we parametrized linear regression models to estimate the dielectric constant directly from the local water density. This is a significant computational shortcut compared to slowly converging calculations relying on total dipole moment fluctuations. The amplitude of the interfacial dielectric constant oscillation can exceed the dielectric constant of the bulk water, suggesting an ice-like frozen state, but only if there are no electrolyte ions. The interfacial accumulation of electrolyte ions causes a decrease in the dielectric constant due to the reduction of water density and re-orientation of water dipoles in ion hydration shells. Finally, we show how to use the computed dielectric properties to estimate SCM's capacitances.

58 GEOSCIENCES↗

DOE OSTI2021

Nuclear Safeguards: Feature Extraction for Machine Learning Enrichment Analysis of UF 6 Cylinders

With the increasing international interest in using nuclear material as sources of energy, comes the growing concern that we may see an increase in proliferation threats. Uranium hexafluoride (UF 6 ) is used in the nuclear fuel cycle for uranium enrichment. Inspectors from the International Atomic Energy Agency (IAEA) monitor the enrichment levels of UF 6 stored in transportation cylinders, to ensure the enrichments match that of a facility operator’s declarations. However, only a characteristic subset of UF 6 cylinders can be measured by inspectors from the IAEA during these inspections. In turn, the inspector may not identify a cylinder whose enrichment levels do not match the facility declarations. Therefore, our team attempts to develop a machine learning network that could determine the enrichment percentage of cylinders as they enter and exit facilities. The machine learning model must be robust against spectral variations due to factors that are internal and external to the UF 6 cylinder. Such factors include but are not limited to the speed and distance of the moving vehicle, cylinder type, cylinder orientation, and fill level. Many features in the spectral continuum can be used to identify and correct for some of these variations. As a consequence, the model must extract numerous features in the continuum. We investigate new approaches in continuum subtraction to verify enrichment percentage: interpolation and extrapolation of lines at the notable characteristic spectral peaks. It was determined that there was no clear answer on which method proved superior, therefore the decision was made to implement both new methods into the current codebase. We currently rely on synthetic training and testing data to analyze the results and fine tune our models but hope to test the system on measured data soon.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

DOE OSTI2023

Stable accretion in young stars: the cases of EX Lupi and TW Hya

ABSTRACT We examine the long-term spectroscopic and photometric variability of EX Lupi and TW Hya, studying the presence of stable accretion and the role it plays in the observed variability. Analysing the velocity modulations of the emission lines with star-melt, we obtain information on the structure of the accretion columns and the disc–star connection. The emission-line radial velocities reveal that TW Hya, like EX Lupi, has a remarkably stable or slow-varying accretion column footprint, locked to the star for several years. The line-emitting regions are non-polar for both EX Lupi and TW Hya, and species with different energies differ in position. In contrast, the continuum emission as observed in the photometry is very variable and can be modelled by hotspot(s) that change over time in phase, shape, temperature, size, and location with respect to the emission-line region. The continuum emission region may not be limited to the stellar surface, especially during episodes of high accretion. The broad-line emission observed in EX Lupi during episodes of increased accretion reveals a further structure, which can be fitted by non-axisymmetric disc in Keplerian rotation inwards of the corotation radius. Since the radial velocity modulation due to accretion footprints is so stable, we used it to search for further velocity modulations. While no residual modulation (other than caused by stellar rotation) is found in these objects, a similar analysis could help to identify young planets/companions. Therefore, determining whether stable accretion footprints are common among young stars is a key to detect young planets.

Sicilia-Aguilar, A.↗

DOE OSTI2023

Scientific machine learning for modeling and simulating complex fluids

The formulation of rheological constitutive equations—models that relate internal stresses and deformations in complex fluids—is a critical step in the engineering of systems involving soft materials. While data-driven models provide accessible alternatives to expensive first-principles models and less accurate empirical models in many engineering disciplines, the development of similar models for complex fluids has lagged. The diversity of techniques for characterizing non-Newtonian fluid dynamics creates a challenge for classical machine learning approaches, which require uniformly structured training data. Consequently, early machine-learning based constitutive equations have not been portable between different deformation protocols or mechanical observables. Here, we present a data-driven framework that resolves such issues, allowing rheologists to construct learnable models that incorporate essential physical information, while remaining agnostic to details regarding particular experimental protocols or flow kinematics. These scientific machine learning models incorporate a universal approximator within a materially objective tensorial constitutive framework. By construction, these models respect physical constraints, such as frame-invariance and tensor symmetry, required by continuum mechanics. We demonstrate that this framework facilitates the rapid discovery of accurate constitutive equations from limited data and that the learned models may be used to describe more kinematically complex flows. This inherent flexibility admits the application of these “digital fluid twins” to a range of material systems and engineering problems. We illustrate this flexibility by deploying a trained model within a multidimensional computational fluid dynamics simulation—a task that is not achievable using any previously developed data-driven rheological equation of state.

Science & Technology - Other Topics↗

DOE OSTI2023

The outburst of the changing-look AGN IRAS 23226-3843 in 2019

IRAS 23226-3843 has previously been classified as a changing-look active galactic nucleus (AGN) based on observations taken in the 1990s in comparison to X-ray data (Swift, XMM-Newton, and NuSTAR) and optical spectra taken after a very strong X-ray decline in 2017. In 2019, Swift observations revealed a strong rebrightening in X-ray and UV fluxes. We aimed to study this outburst in greater detail. We took follow-up Swift, XMM-Newton, and NuSTAR observations of IRAS 23226-3843 together with optical spectra (SALT and SAAO 1.9 m telescope) from 2019 until 2021. IRAS 23226-3843 showed a strong X-ray and optical outburst in 2019. It varied in the X-ray continuum by a factor of 5 and in the optical continuum by a factor of 1.6 within two months. This corresponds to a factor of 3 after correction for the host galaxy contribution. The Balmer and Fe II emission-line intensities showed comparable variability amplitudes during the outburst in 2019. The Hα emission-line profiles of IRAS 23226-3843 changed from a blue-peaked profile in the years 1997 and 1999 to a broad double-peaked profile in 2017 and 2019. However, there were no major profile variations in the extremely broad double-peaked profiles despite the strong intensity variations in 2019. One year after the outburst, IRAS 23226-3843 changed its optical spectral type and became a Seyfert type 2 object in 2020. Blue outflow components are present in the optical Balmer lines and in the Fe band in the X-rays. A deep broadband XMM-Newton/NuSTAR spectrum was taken during IRAS 23226-3843’s maximum state in 2019. This spectrum is qualitatively very similar to a spectrum taken in 2017, but by a factor of 10 higher. The soft X-ray band appears featureless. The soft excess is well modeled with a Comptonization model. A broadband fit with a power-law continuum, Comptonized soft excess, and Galactic absorption gives a good fit to the combined EPIC-pn and NuSTAR spectrum. In addition, we see a complex and broadened Fe K emission-line profile in the X-rays. The changing-look character in IRAS 23226-3843 is most probably caused by changes in the accretion rate – based on the short-term variations on timescales of weeks to months.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Sierra/SolidMechanics 4.58 User's Guide

Sierra / SolidMechanics (Sierra / SM) is a Lagrangian, three-dimensional code for finite element analysis of solids and structures. It provides capabilities for explicit dynamic, implicit quasistatic and dynamic analyses. The explicit dynamics capabilities allow for the efficient and robust solution of models with extensive contact subjected to large, suddenly applied loads. For implicit problems, Sierra / SM uses a multi-level iterative solver, which enables it to effectively solve problems with large deformations, nonlinear material behavior, and contact. Sierra / SM has a versatile library of continuum and structural elements, an d a large library of material models. The code is written for parallel computing environments enabling scalable solutions of extremely large problems for both implicit and explicit analyses. It is built on the SIERRA Framework, which facilitates coupling with other SIERRA mechanics codes . This document describes the functionality and input syntax for Sierra / SM.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Multiphysics Modeling Framework for Soft PVC Gel Sensors with Experimental Comparisons

Polyvinyl chloride (PVC) gels have recently been found to exhibit mechanoelectrical transduction or sensing capabilities under compressive loading applications. This phenomenon is not wholly understood but has been characterized as an adsorption-like phenomena under varying amounts and types of plasticizers. A different polymer lattice structure has also been tested, thermoplastic polyurethane, which showed similar sensing characteristics. This study examines mechanical and electrical properties of these gel sensors and proposes a mathematical framework of the underlying mechanisms of mechanoelectrical transduction. COMSOL Multiphysics is used to show solid mechanics characteristics, electrostatic properties, and transport of interstitial plasticizer under compressive loading applications. The solid mechanics takes a continuum mechanics approach and includes a highly compressive Storakers material model for compressive loading applications. The electrostatics and transport properties include charge conservation and a Langmuir adsorption migration model with variable diffusion properties based on plasticizer properties. Results show both plasticizer concentration gradient as well as expected voltage response under varying amounts and types of plasticizers. Experimental work is also completed to show agreeance with the modeling results.

42 ENGINEERING↗

DOE OSTI2022

A Brief Overview of Reactive Transport Codes Used in CO 2 Applications

Reactive transport modelling consists of computational and numerical models that describe the coupled physical, chemical, mechanical, and biological processes interacting with each other over a broad range of spatial and temporal scales. This report reviews some of the reactive transport codes available in the literature for carbon dioxide (CO 2 ) applications and can assist the scientific community with the applicability of each code to solve a specific problem. Approaches implemented in these codes can be categorized in three groups: (i) continuum-scale, (ii) pore-scale, and (iii) hybrid-scale approaches. Theoretical foundations, numerical implementations, and application examples using the models are described. In the continuum approach, flow and transport are formulated in terms of a representative control volume of the medium and coupled with geochemical reactions. The medium is characterized by bulk parameters such as porosity, permeability, or reactive surface area. In the pore-scale approach, each point of space in the pore network is occupied by either a fluid or solid phase. The pore-scale approach requires an exact knowledge of the spatial and temporal phase distribution. In the hybrid-scale approach, media are described by multiple characteristic length-scales, with some regions using a pore-scale approach while others used a continuum approach.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2026

Plasma screening in mid-charged ions observed by K-shell line emission

Dense plasma environment affects the electronic structure of ions via variations of the microscopic electrical fields, also known as plasma screening. This effect can be either estimated by simplified analytical models, or by computationally expensive and to date unverified numerical calculations. We have experimentally quantified plasma screening from the energy shifts of the bound-bound transitions in matter driven by the x-ray free electron laser (XFEL). This was enabled by identification of detailed electronic configurations of the observed Kα, Kβ and Kγ lines. This work paves the way for improving plasma screening models including connected effects like ionization potential depression and continuum lowering, which will advance the understanding of atomic physics in the Warm Dense Matter regime.

Physics - Plasma physics↗

DOE OSTI2025

Forecasting the Ly-𝛼 forest × CMB lensing bispectrum from ACT, SO, S4, and DESI

We forecast the sensitivity of future correlations between cosmic microwave background (CMB) lensing and the Lyman-𝛼 forest power spectrum. This squeezed-limit bispectrum probes the connection between the Lyman-𝛼 transmission and the underlying large-scale matter density field. We recover the measured signal-to-noise (SNR) ratio obtained with Planck × BOSS, and forecast a SNR of 10.3, 14.8 and 20.2 for DESI combined with ACT, SO, and CMB-S4, respectively. For DESI and SO/CMB-S4, the correlation should be detectable at SNR ≥ 5 in 5 redshift bins, useful for distinguishing our signal from several contaminants. Indeed, our forecast is affected by large theoretical uncertainties in the current modeling of the signal and its contaminants, such as continuum misestimation bias or damped-Lyman-𝛼 absorbers. Quantifying these more accurately will be crucial to enable a reliable cosmological interpretation of this observable in future measurements. We attempt to enumerate the features required in future Lyman-𝛼 forest simulations required to do so.

Cosmic microwave background↗

DOE OSTI2024

Phylogeny, morphology, virulence, ecology, and host range of Ordospora pajunii (Ordosporidae), a microsporidian symbiont of Daphnia spp.

ABSTRACT The impacts of microsporidia on host individuals are frequently subtle and can be context dependent. A key example of the latter comes from a recently discovered microsporidian symbiont of Daphnia , the net impact of which was found to shift from negative to positive based on environmental context. Given this, we hypothesized low baseline virulence of the microsporidian; here, we investigated the impact of infection on hosts in controlled conditions and the absence of other stressors. We also investigated its phylogenetic position, ecology, and host range. The genetic data indicate that the symbiont is Ordospora pajunii , a newly described microsporidian parasite of Daphnia . We show that O. pajunii infection damages the gut, causing infected epithelial cells to lose microvilli and then rupture. The prevalence of this microsporidian could be high (up to 100% in the lab and 77% of adults in the field). Its overall virulence was low in most cases, but some genotypes suffered reduced survival and/or reproduction. Susceptibility and virulence were strongly host-genotype dependent. We found that North American O. pajunii were able to infect multiple Daphnia species, including the European species Daphnia longispina , as well as Ceriodaphnia spp. Given the low, often undetectable virulence of this microsporidian and potentially far-reaching consequences of infections for the host when interacting with other pathogens or food, this Daphnia–O. pajunii symbiosis emerges as a valuable system for studying the mechanisms of context-dependent shifts between mutualism and parasitism, as well as for understanding how symbionts might alter host interactions with resources. IMPORTANCE The net outcome of symbiosis depends on the costs and benefits to each partner. Those can be context dependent, driving the potential for an interaction to change between parasitism and mutualism. Understanding the baseline fitness impact in an interaction can help us understand those shifts; for an organism that is generally parasitic, it should be easier for it to become a mutualist if its baseline virulence is relatively low. Recently, a microsporidian was found to become beneficial to its Daphnia hosts in certain ecological contexts, but little was known about the symbiont (including its species identity). Here, we identify it as the microsporidium Ordospora pajunii . Despite the parasitic nature of microsporidia, we found O. pajunii to be, at most, mildly virulent; this helps explain why it can shift toward mutualism in certain ecological contexts and helps establish O. pajunii is a valuable model for investigating shifts along the mutualism-parasitism continuum.

59 BASIC BIOLOGICAL SCIENCES↗