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.

SourceResearch connectionFrom yearThrough year

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

At least 487 records · Page 27

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Consistent linearization of micromorphic continuum theories

We outline the procedure of consistent linearization and apply it to the micromorphic, microstretch, and micropolar theories of continua. This yields tractable linear theories for nonlinearly elastic microstructured materials undergoing finite deformations. The results may be readily utilized in computational mechanics, stability and bifurcation analyses, and small-deformation problems in the context of these types of continua. Our results generalize those existing in current literature and facilitate their recovery upon incorporating appropriate kinematic and constitutive assumptions.

McAvoy, Ryan C.↗

DOE OSTI2025

A Geometric Volume of Fluid-Based Multiphase Flow Solver Extension to the Reacting Flow Solver, PeleLM

A new algorithm is presented to simulate multiphase flows with surface tension in a pathway for spray combustion simulation. The algorithm combines capabilities from two open-source packages, including the interface reconstruction library (IRL), a library of computational geometry routines to enable the volume of fluid (VOF) method, and PeleLM, a solver for the reacting Navier-Stokes equations. Additionally, surface tension is implemented using the continuum surface force (CSF) model with an improved height function technique in the volume fraction field. Spurious errors in volume fraction arising from our combined strategy are corrected through a topology-based method that improves both numerical stability and accuracy. Multiple validation simulations are conducted, including (i) translations and rotations of Zalesak's disk, (ii) a stationary circular droplet with surface tension, (iii) an oscillating elliptical droplet, and (iv) three-dimensional deformation of a spherical droplet. Results indicate that the combined scheme retains the favorable properties of each of the component algorithms.

42 ENGINEERING↗

DOE OSTI2023

The evolution of pit morphology and growth kinetics in aluminum during atmospheric corrosion

Abstract Assessing the lifetimes of alloys in humid, corrosive environments requires growth kinetic information regarding individual instances of damage, e.g. pit growth rates. Corrosion rates measured at the continuum scale using mass change convolute the rate of pit nucleation and growth, providing limited information on local kinetics. The current study used in-situ X-ray computed tomography to measure growth rates of individual pits in aluminum over 100 h of exposure in a humid, chloride environment. While pits grew at relatively constant rates over the first hours after nucleation, significant growth-rate nonlinearities subsequently occurred. These were linked to both droplet spreading, which altered the cathode size, and changes in the mode of pit growth. Pit morphology appeared to influence the dominant growth mode and the duration of pit growth. Post-mortem serial sectioning revealed pits preferentially attacked grain-boundary triple junctions and dislocation boundaries.

Noell, Philip J. (ORCID:0000000346998999)↗

NASA NTRS1986

A review of modelling techniques for the open and closed-loop dynamics of large space systems

This paper reviews the steps in the development of mathematical models that can be used to simulate the in-orbit dynamic behavior of large flexible systems. A general continuum formulation is compared with the hybrid coordinate formulation and also a finite element representation of the total system. A review of structural analysis routines emphasizes the use of computer generated graphics to help understand the different modal elastic shape functions of complex systems. Numerical techniques employed to synthesize shape and attitude control laws are summarized. Finally, the modeling of environmental disturbance torques due to the interaction of solar radiation pressure on vibrating and thermally deflected systems is discussed.

Bainum, Peter M.↗

NASA NTRS1994

Computing Incompressible Flows With Free Surfaces

RIPPLE computer program models transient, two-dimensional flows of incompressible fluids with surface tension on free surfaces of general shape. Surface tension modeled as volume force derived from continuum-surface-force model, giving RIPPLE both robustness and accuracy in modeling surface-tension effects at free surface. Also models wall adhesion effects. Written in FORTRAN 77.

Kothe, D.↗

DOE OSTI2022

The performance of ensemble-based free energy protocols in computing binding affinities to ROS1 kinase

Optimization of binding affinities for compounds to their target protein is a primary objective in drug discovery. Herein we report on a collaborative study that evaluates a set of compounds binding to ROS1 kinase. We use ESMACS (enhanced sampling of molecular dynamics with approximation of continuum solvent) and TIES (thermodynamic integration with enhanced sampling) protocols to rank the binding free energies. The predicted binding free energies from ESMACS simulations show good correlations with experimental data for subsets of the compounds. Consistent binding free energy differences are generated for TIES and ESMACS. Although an unexplained overestimation exists, we obtain excellent statistical rankings across the set of compounds from the TIES protocol, with a Pearson correlation coefficient of 0.90 between calculated and experimental activities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1993

Computational strategies in the dynamic simulation of constrained flexible MBS

This research focuses on the computational dynamics of flexible constrained multibody systems. At first a recursive mapping formulation of the kinematical expressions in a minimum dimension as well as the matrix representation of the equations of motion are presented. The method employs Kane's equation, FEM, and concepts of continuum mechanics. The generalized active forces are extended to include the effects of high temperature conditions, such as creep, thermal stress, and elastic-plastic deformation. The time variant constraint relations for rolling/contact conditions between two flexible bodies are also studied. The constraints for validation of MBS simulation of gear meshing contact using a modified Timoshenko beam theory are also presented. The last part deals with minimization of vibration/deformation of the elastic beam in multibody systems making use of time variant boundary conditions. The above methodologies and computational procedures developed are being implemented in a program called DYAMUS.

Amirouche, F. M. L.↗

DOE OSTI2021

Reusability First: Toward FAIR Workflows

The FAIR principles of open science (Findable, Accessible, Interoperable, and Reusable) have had transformative effects on modern large-scale computational science. In particular, they have encouraged more open access to and use of data, an important consideration as collaboration among teams of researchers accelerates and the use of workflows by those teams to solve problems increases. How best to apply the FAIR principles to workflows themselves, and software more generally, is not yet well understood. We argue that the software engineering concept of technical debt management provides a useful guide for application of those principles to workflows, and in particular that it implies reusability should be considered as ‘first among equals’. Moreover, our approach recognizes a continuum of reusability where we can make explicit and selectable the tradeoffs required in workflows for both their users and developers.To this end, we propose a new abstraction approach for reusable workflows, with demonstrations for both synthetic workloads and real-world computational biology workflows. Through application of novel systems and tools that are based on this abstraction, these experimental workflows are refactored to rightsize the granularity of workflow components to efficiently fill the gap between end-user simplicity and general customizability. Our work makes it easier to selectively reason about and automate the connections between trade-offs across user and developer concerns when exposing degrees of freedom for reuse. Additionally, by exposing fine-grained reusability abstractions we enable performance optimizations, as we demonstrate on both institutional-scale and leadership-class HPC resources.

Wolf, Matthew↗

DOE OSTI2026

When Do Band Gap Calculations Agree with Experiments in Monolayer-Protected Cu 14 and Au 20 Atomically Precise Nanoclusters? A (TD)-DFT Comparison of HOMO–LUMO, Fundamental, Optical, and Electrochemical Energy Gaps

In view of the tremendous progress in atomically precise metal nanoclusters where electrochemical and optical energetics are routinely supported by computations to establish structure−function correlations, we explore the relationship between different protocols for measuring and computing band gaps of two distinct organic ligand-protected nanoclusters: [Cu 14 H 10 (MBN) 3 (PPh 3 ) 8 ] + and Au 20 (TBBT) 16 . Through UV/visible spectroscopy and differential pulse voltammetry, we measure optical and electrochemical band gaps in those systems. We then compare these experimentally determined gaps to HOMO−LUMO gaps, fundamental gaps, vertical excitation energies, and E o ox − E o red potentials computed using different density functional theory (DFT) or time-dependent DFT (TDDFT) methods. Specifically, in both copper and gold nanoclusters, we test the effect of truncating inert ligands from the model and compare density functionals with varying degrees of Hartree−Fock (HF) exchange from 0 to 50%, range-separated hybrids with a varying long-range tuning parameter, different correlation functionals, basis sets, and (equilibrium and nonequilibrium) continuum solvation models. Despite having different frontier orbital characters (the copper nanocluster has a metal-to-ligand charge transfer character while the gold nanocluster has metal-centered frontier orbitals), both nanoclusters display a similar sensitivity of the HOMO−LUMO gap to the HF exchange that is partially mitigated when computing the fundamental, optical, and electrochemical gaps. Other factors, such as the nature of the correlation functional, basis set, and geometry relaxation, have a considerably smaller effect on computed band gaps in these systems. Overall, this work provides guidelines for factors of varied importance for correlating computed and experimental band gap values.

Chemical calculations↗

DOE OSTI2025

Comparison of interlaminar damage modeling strategies for hybrid composite/aluminum laminates subjected to low-velocity impact

Low-velocity impact of hybrid metal-composite structures was investigated experimentally and computationally. Composite laminates consisting of 2D woven glass fiber reinforced polymer (GFRP) and carbon fiber reinforced polymer (CFRP) were joined with a 6061-T6 aluminum plate using an epoxy adhesive. Two variations of the structure were studied; one consisting of all plies oriented at 0° and one consisting of all plies oriented at 45°. A drop tower was used to impact structures at a range of energies, including energies above and below the threshold at which the aluminum layer was perforated. Numerical simulations were implemented using Sierra/SM, an in-house transient dynamics finite element code developed at Sandia National Laboratories. A Hosford plasticity model was used to describe the response of the aluminum layer. A newly implemented orthotropic continuum damage mechanics (CDM) constitutive model was used to represent the composite laminate. This 3D-CDM model was compared to a cohesive zone model (2D-CDM/CZM) to investigate efficacy of aluminum perforation energy prediction, delamination prediction, and computational cost. Accuracy of each model was evaluated using the experimental results. Each showed good agreement with the tests for both the force and velocity histories, as well as the observed damage mechanisms. The 2D-CDM/CZM model was marginally more accurate in capturing both the composite and aluminum behavior — this model averaged error percentages of -11.2% and 10.8% for residual velocity and peak force, respectively. Meanwhile, the 3D-CDM model predictions yielded average error percentages of -35.5% (velocity) and 22.6% (force). However, the 3D-CDM model generally resulted in a decreased computational cost; the average run time was 14% shorter than the 2D-CDM/CZM model and 3x as many timesteps per hour were computed using the same computational resources. In conclusion, new experimental data on the impact and perforation resistance of metal-composite laminates is presented in addition to numerical predictions of the impact behavior.

Carbon fiber↗

DOE OSTI2025

Competition between roughness and strength for scale-dependent surfaces

Rocks famously have scale-dependent strength, yet the actual dependence is notoriously hard to measure or incorporate into any theoretical framework. Natural rough surfaces present an opportunity to solve the problem. Surfaces sliding in shear evolve as protrusions collide. These asperities can deform or break, thus creating a new surface shape. In particular, natural surfaces have roughness at all scales as well as scale-dependent strength. Based on a scaling analysis, we have previously suggested that the scale-dependent aspect ratio of steady-state surfaces should be proportional to the scale-dependent shear strain at yield. If true, scale-dependent strength could easily be inferred from natural surfaces. Thus, moving beyond the scaling argument to a rigorous treatment of scale-dependent strength for multiscale rough surfaces in shear is important. However, analytic frameworks for analyzing multiscale problems are challenging, as conventional continuum mechanics typically involves a single value for a material property across scales. Here, in this work, we build on the formalism of Persson (2001) that presents a method to compute contact area for rough surfaces with a prescribed topographic spectrum using a stochastic differential equation. The Persson formalism allows for plastic yield under normal loading of otherwise elastic materials and leaves open the possibility of scale-dependent yield stress. In this study, we pursue this route to develop a theory and numerical results for the yielding of a rough, elastoplastic surface with scale-dependent yield stress. Here, we examine surfaces for which the power spectrum of the topography 𝐶 and yield stress 𝑌 follow power laws as a function of scale 𝜆, such that 𝐶∼𝜆 −𝑚 and 𝑌∼𝜆 −𝑛 , respectively. In this formal treatment of the problem, we focus on surfaces in contact and the resulting yield and do not impose shear. Numerical solutions show that the deviation from the elastic scaling solution is bounded as expected by the prior 1D heuristic scaling argument that anticipates the Hurst exponent as 1−𝑛. We also show that the plasticity is expected to erode the contacts if 𝑚 is lower than 𝑛−3, which corresponds to a Hurst exponent lower than 1−𝑛/2. This result is rigorously sound for 2D, i.e., realistic surfaces, and quantitatively different than the prior scaling argument. The theory now permits a correspondingly quantitative approach to interpreting natural surfaces.

elasticity↗

DOE OSTI2023

Electroweak box diagrams on the lattice for pion and neutron decay

CKM matrix is unitary by construction in the standard model(SM). The recent analyses on the first row of CKM matrix show ≈3σ tension with unitarity. Nonperturbative calculations of the radiative corrections can reduce the theory uncertainty in CKM matrix elements. Here we compute the electroweak box contribution to the pion and kaon β decays using seven Nf = 2 + 1 + 1 HISQ-Clover lattice with various pion mass and lattice spacing. The continuum and chiral limit is taken using the leading dependence on M π and a, where M π extrapolation is taken to the physical pion mass and SU(3) symmetric mass for pion and kaon box contribution, respectively. Our results are $\square^{VA}_{γW}|_{π}$ = 2.820(28) × 10 –3 and $\square^{VA}_{γW}|_{K}$ = 2.384(17) × 10 –3 .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Surrogate Constructed Scalable Circuits ADAPT-VQE in the Schwinger model

Inspired by recent advancements of simulating periodic systems on quantum computers, we develop a new approach, (SC)$^2$-ADAPT-VQE, to further advance the simulation of these systems. Our approach extends the scalable circuits ADAPT-VQE framework, which builds an ansatz from a pool of coordinate-invariant operators defined for arbitrarily large, though not arbitrarily small, volumes. Our method uses a classically tractable ``Surrogate Constructed'' method to remove irrelevant operators from the pool, reducing the minimum size for which the scalable circuits are defined. Bringing together the scalable circuits and the surrogate constructed approaches forms the core of the (SC)$^2$ methodology. Our approach allows for a wider set of classical computations, on small volumes, which can be used for a more robust extrapolation protocol. While developed in the context of lattice models, the surrogate construction portion is applicable to a wide variety of problems where information about the relative importance of operators in the pool is available. As an example, we use it to compute properties of the Schwinger model - quantum electrodynamics for a single, massive fermion in $1+1$ dimensions - and show that our method can be used to accurately extrapolate to the continuum limit.

Gustafson, Erik [RIACS, Mtn. View] (ORCID:00000001↗

NASA NTRS1987

Positron annihilation at the Galactic center

The efficiency of production of positrons by gamma-gamma interaction under conditions which may occur in the Galactic center. Their annihilation line and three-photon annihilation continuum contribute to the emitted radiation. If a geometrical model is assumed, it will be possible to determine the parameters of the gamma-ray source region. By comparing the computed spectra to observations, a rough estimate is made using extant data.

Carrigan, B. J.↗

NASA NTRS1992

Role of shock dynamics in transonic flutter

A computational study of the influence of shock motion on flutter and divergence in transonic flow is presented. The numerical scheme models the entire fluid-structure system as a single continuum dynamics problem, by using a mixed Eulerian-Lagrangian formulation. No assumptions of small displacements are made, but the effect of viscosity is neglected. The results from this study indicate that the shock dynamics gives rise to limit cycles and highly nonlinear aeroelastic phenomena, such as weak divergence and flutter-divergence interactions. Although the shocks typically are destabilizing at the linear flutter boundary, they often have a strongly stabilizing effect for moderate-amplitude motions. The shocks are thus capable of quenching an emerging bending-torsion flutter motion and turning it into limit cycle flutter. The usefulness of classical flutter and divergence boundary diagrams is severely limited in transonic flow, because much of the global dynamic stability information is lost in such a presentation.

Bendiksen, Oddvar O.↗

NASA NTRS1993

One-parameter scaling and exponential-sum fitting for water vapor and CO2 infrared transmission functions

A medium-sized band model for water vapor and CO2 absorption is developed using the one-parameter scaling approximation. The infrared spectrum is divided into 10 bands. The Planck-weighted diffuse transmittance is reduced to a function dependent only upon the scaled absorber amount and fit by an exponential sum. By selecting specific sets of absorption coefficients for exponential-sum fitting, computations of fluxes and cooling rate are made very fast. Compared to a broadband model, the accuracy, speed, and versatility are all enhanced. With absorption due to water vapor line, continuum, CO2 as well as O3 included, the parameterization introduces an error of less than 1.5 W/sq m in fluxes and less than 0.15 C/day in the tropospheric and lower stratospheric cooling rates.

Chou, Ming-Dah↗

NASA NTRS1993

Two color satellite laser ranging upgrades at Goddard's 1.2m telescope facility

The ranging laboratory at Goddard's 1.2 m telescope tracking facility has recently been upgraded to include a single photoelectron sensitive Hamamatsu streak camera-based range receiver which uses doubled and tripled Nd:YAG frequencies for satellite laser ranging. Other ranging system upgrades include a new continuum laser, which will deliver up to 30 millijoules (mJ) at both 532 and 355 nm at a pulsewidth of 30 picoseconds (FWHM), and replacement of both ranging and tracking computers with COMPAQ 386 based systems. Preliminary results using a photomultiplier-tube based receiver and waveform digitizer indicate agreement within the accuracy of the measurement with the theoretical Marini and Murray model for atmospheric refraction. Two color streak camera measurements are used to further analyze the accuracy of these and other atmospheric refraction models.

Zagwodzki, Thomas W.↗

NASA NTRS1994

A quantitative analysis of IRAS maps of molecular clouds

We present an analysis of IRAS maps of five molecular clouds: Orion, Ophiuchus, Perseus, Taurus, and Lupus. For the classification and description of these astrophysical maps, we use a newly developed technique which considers all maps of a given type to be elements of a pseudometric space. For each physical characteristic of interest, this formal system assigns a distance function (a pseudometric) to the space of all maps: this procedure allows us to measure quantitatively the difference between any two maps and to order the space of all maps. We thus obtain a quantitative classification scheme for molecular clouds. In this present study we use the IRAS continuum maps at 100 and 60 micrometer(s) to produce column density (or optical depth) maps for the five molecular cloud regions given above. For this sample of clouds, we compute the 'output' functions which measure the distribution of density, the distribution of topological components, the self-gravity, and the filamentary nature of the clouds. The results of this work provide a quantitative description of the structure in these molecular cloud regions. We then order the clouds according to the overall environmental 'complexity' of these star-forming regions. Finally, we compare our results with the observed populations of young stellar objects in these clouds and discuss the possible environmental effects on the star-formation process. Our results are consistent with the recently stated conjecture that more massive stars tend to form in more 'complex' environments.

Wiseman, Jennifer J.↗