SEARCH · Engineering Papers
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.
Here, this article introduces ways to identify dynamic system models using measurement data. In power system analysis, a static model represents the time-invariant input-output relationship of a system, while a dynamic model describes the behavior of the system over time. For example, how will a system transit from one steady-state operation point to another?
Inspired by the self-organization of unicellular species into multicellular organisms, we use theory and simulation to design a system of mobile, active microcapsules that produce chemicals according to a catalytic reaction network (CRN). In solution, the catalytic reactions generate a force that drives the flow of the surrounding fluid. The convective flow, in turn, transports the capsules to new chemical surroundings with each successive translation. Consequently, the chemical signal produced by a cluster of capsules is critically dependent on the proximity and spatial configuration of the neighboring capsules. If all of the capsules and chemical products involved in the CRN lie within sufficient proximity, then the system develops oscillatory behavior that drives the dynamic self-assembly of capsules into larger clusters that are capable of collective action. Finally, these simulations indicate potential chemo-mechanic mechanisms that enable single cellular units, such as ameba, to organize into multicellular life forms.
Elucidating the nature of spin excitations is important to understanding the mechanism of superconductivity in the iron pnictides. Motivated by recent inelastic neutron scattering measurements in the nearly 100% detwinned BaFe 2 As 2 , we study the spin dynamics of an $\textit{S}$ = 1 frustrated bilinear-biquadratic Heisenberg model in the antiferromagnetic phase with wave vector (π, 0). The biquadratic interactions are treated in a dynamical way using a flavor-wave theory in an $\textit{SU}$(3) representation. Besides the dipolar spin wave (magnon) excitations, the biquadratic interactions give rise to quadrupolar excitations at high energies. In this work, we find that the quadrupolar wave significantly influences, in an energy dependent way, the anisotropy between the spin excitation spectra along the (π, 0) and (0, π) directions in the wave vector space. Our theoretical results capture the essential behavior of the spin dynamics measured in the antiferromagnetic phase of the detwinned BaFe 2 As 2 . More generally, our results underscore the importance of electron correlation effects for the microscopic physics of the iron pnictides.
Two methods for determining the location of and load level to produce instability of compressed cylindrical shells are presented. The first relates the variation in the wall normal stiffness as a function of applied compressive force to the critical load. It uses the distribution of stiffness over the surface of the shell as a guide to buckle location. The second method associates the local dynamic mass with instability behavior. The test data presented show that either method will give excellent prediction capability from low-load-level data for shells of orthodox form. Neither method appears to apply to spirally stiffened shells. This is thought to be due to the fact that there is a substantial difference between the buckle pattern under axial compression and the imperfection shape induced by the normal displacement which is used to ascertain the wall stiffness and the dynamic mass.
The present study employs molecular dynamics simulations of Ni-based superalloy to investigate the creep behavior under uniaxial compression test. Dislocation dynamics is analyzed for the nickel-based single crystal superalloy with the presence of void and with varying the distribution of gamma-prime phase The results show that multiple-void systems are more prone to yield than single-void systems and single-void systems are more prone to yield than the system withoutvoid. From the simulations, it has been determined that the creep mechanism in Ni/Ni3Al is subject to change on the applied stress depending on the distribution of gamma-prime phases change. Dislocation behavior is also observed in addition to the time evolution of dislocation nucleation and growth to quantify the propensity and mechanism of creep in nickel-based single crystal superalloy. The studies confirmed that dislocation length prior to plastic deformation decreased as the number of voids increased. In addition, the deformation behavior for finegrained nickel-based polycrystalline superalloy upon grain size variation is also investigated. Our findings show that both the movement mechanisms of the grains and grain boundaries and the relationship between grain size and yield strength are highly affected by the grain size of the materials. All the databases have then been classified to accurately predict the target class vi machine learning using the data mining tool WEKA.
The application of nonlinear dynamics to improve the understanding of numerical uncertainties in computational fluid dynamics (CFD) is reviewed. Elementary examples in the use of dynamics to explain the nonlinear phenomena and spurious behavior that occur in numerics are given. The role of dynamics in the understanding of long time behavior of numerical integrations and the nonlinear stability, convergence, and reliability of using time-marching, approaches for obtaining steady-state numerical solutions in CFD is explained. The study is complemented with spurious behavior observed in CFD computations.
The application of nonlinear dynamics to improve the understanding of numerical uncertainties in computational fluid dynamics (CFD) is reviewed. Elementary examples in the use of dynamics to explain the nonlinear phenomena and spurious behavior that occur in numerics are given. The role of dynamics in the understanding of long time behavior of numerical integrations and the nonlinear stability, convergence, and reliability of using time-marching approaches for obtaining steady-state numerical solutions in CFD is explained. The study is complemented with examples of spurious behavior observed in CFD computations.
The accumulation of gas atoms in tungsten is a topic of long-standing interest to the plasma-facing materials community due the metal's use as a divertor material in some tokamak fusion reactors. The nucleation and growth of He/H gas bubbles (along with their isotopes) can result from impinging fluxes of these gases which give rise to damage at the W divertor surface. The inclusion of He or H in W has been studied extensively by the community, finding that He bubbles modify the surface through periodic dislocation punching and bursting mechanisms while H bubbles impact the metal through plastic-strain induced material failure. However, the mechanisms which are present during the combined flux of both He and H is not well-studied atomistically. Motivated by this, an atomistic modeling study is conducted using molecular dynamics to assess the behavior of mixed concentration He:H bubbles in W. Here we find that the introduction of H into a growing He bubble results in a dramatic change in the nature and presence of dislocation loops which are typically generated via dislocation punching in over-pressurized He bubbles. Most notably, at high H concentrations, there is a switchover in energetic favorability from glissile 1/2<111> dislocations to sessile <100> dislocations. This thermodynamic crossover could imply a significant reduction in W surface morphology changes than with pure He bubbles and, additionally, we show this to have implications on the trapping of H in the bubbles and their associated dislocations.
Per-/polyfluoroalkyl substances (PFAS) have been commonly used over several decades for a variety of products and are very persistent in the environment. However, not much is known about their direct effects on aquatic invertebrates and their ecosystems. We examined the survival, behavior, development, and predation susceptibility of mosquito larvae Culex quinquefasciatus exposed to concentrations of perfluorooctanesulfonic acid (PFOS) ranging from 0.02 to 453.7 µg/L. PFOS exposure resulted in reduced larval survival, with a 48 hour LC50 (concentration with 50 % mortality) of 255.99 µg/L. PFOS exposure also resulted in reduced developmental success and slower maturation to adulthood (and thus slower emergence from the water) compared to control larvae. PFOS also resulted in delays in reaction to prodding stimuli, which were meant to simulate a predator attack, and longer reactions to prodding. Larvae exposed to PFOS also spent more time at the bottom of the water column, rather than at the surface where respiration takes place. Effects of larval mosquito PFOS exposure on predation by unexposed damselfly naiads (genus Ischnura) were not detected. Overall, this research suggests that PFOS impacts mosquito larvae survival, behavior, development, and adult emergence, which could have important implications for food webs or public health given the role of mosquitos as disease vectors.
Nuclei are used for high-precision tests of the Standard Model and for studies of physics beyond the Standard Model. Without a thorough understanding of nuclei, we will not be able to meaningfully interpret the growing body of experimental data nor will we be able to disentangle new physics signals from underlying nuclear effects. This calls for accurate calculations of nuclear structure and reactions. In this work, we focus on electroweak decays in nuclei with mass number A ≤ 10 and report on ab initio Quantum Monte Carlo calculations of reduced matrix elements entering beta decays and electron captures in nuclei with mass number A ≤ 10. The many-body wave functions are calculated using selected Norfolk two- and three-nucleon potential models and associated one- and two-body axial currents at tree-level obtained from a chiral effective field theory with pions, nucleons, and Δ. The agreement with the experimental data is satisfactory except for transitions in A = 8 nuclei. In this specific case, the theory significantly underpredicts the experimental data, which indicates the need of further improvements in the corresponding nuclear wave functions. In this study, emphasis is placed on the contributions of two-body axial currents that are carefully analyzed using two-body transition densities. This allow us to study the spatial distribution and short-range behavior of two-body dynamics. In particular, the transition densities when scaled to peak at 1.0 exhibit universal short-range behavior across the considered nuclei, while they differ in the long-range tails.
Some grouping is necessary when constructing a Leslie matrix model because it involves discretizing a continuous process of births and deaths. Here, the level of grouping is determined by the number of age classes and frequency of sampling. It is largely unknown what is lost or gained by using fewer age classes, and I address this question using aggregation theory. I derive an aggregator for a Leslie matrix model using weighted least squares, determine what properties an aggregated matrix inherits from the original matrix, evaluate aggregation error, and measure the influence of aggregation on asymptotic and transient behaviors. To gauge transient dynamics, I employ reactivity of the standardized Leslie matrix. I apply the aggregator to 10 Leslie models developed for animal populations drawn from a diverse set of species. Several properties are inherited by the aggregated matrix: (a) it is a Leslie matrix; (b) it is irreducible whenever the original matrix is irreducible; (c) it is primitive whenever the original matrix is primitive; and (d) its stable population growth rate and stable age distribution are consistent with those of the original matrix if the least squares weights are equal to the original stable age distribution. In the application, depending on the population modeled, when the least squares weights do not follow the stable age distribution, the stable population growth rate of the aggregated matrix may or may not be approximately consistent with that of the original matrix. Transient behavior is lost with high aggregation.
The STS Centaur was designed to be a high energy upper stage for use with the Space Shuttle. Two versions were designed under development when the program was cancelled. The first version, designated G-prime, was for planetary missions. The second version, designated G, was to place spacecraft in geosynchronous orbit. As a part of the STS Centaur finite-element model verification effort, test articles of both versions were subjected to a series of static tests. In addition the Centaur G-prime test article was subjected to a series of dynamic tests including a modal survey. Both the static and dynamic tests showed that nonlinearities existed in the Centaur and its support system. The support system included flight-like latches. The nonlinearities were particularly apparent in tests that loaded the forward support structure of the Centaur. These test results were used to aid in the development of two improved finite-element models. The first was a linear model, while the second contained nonlinear elements at the boundaries. Results from both models were compared with the transient response obtained from a step-relaxation or twang test. The linear model was able to accurately match the low frequency response found in the test data. However, only the nonlinear model was able to match higher frequency response that was present in some of the test data. In addition the nonlinear model was able to predict other nonlinear behavior such as the dynamic jump that occurs in systems with nonlinear stiffness.
The STS Centaur was designed to be a high energy upper stage for use with the Space Shuttle. Two versions were designed and under development when the program was cancelled. The first version, designated G-prime, was designed for planetary missions. The second version, designated G, was designed to place spacecraft in geosynchronous orbit. As a part of the STS Centaur finite-element model verification effort, test articles of both versions were subjected to a series of static tests. In addition the Centaur G-prime test article was subjected to a series of dynamic tests including a modal survey. Both the static and dynamic tests showed that nonlinearities existed in the Centaur and its support system. The support system included flight-like latches. The nonlinearities were particularly apparent in tests that loaded the forward support structure of the Centaur. These test results were used to aid in the development of two improved finite-element models. The first was a linear model, while the second contained nonlinear elements at the boundaries. Results from both models were compared with the transient response obtained from a step-relaxation or twang test. The linear model was able to accurately match the low frequency response found in the test data. However, only the nonlinear model was able to match higher frequency response that was present in some of the test data. In addition the nonlinear model was able to predict other nonlinear behavior such as the dynamic 'jump' that occurs in systems with nonlinear stiffness.
Understanding the dynamics of circulating tumor cell (CTC) behavior within the vasculature has remained an elusive goal in cancer biology. To elucidate the contribution of hydrodynamics in determining sites of CTC vascular colonization, the physical forces affecting these cells must be evaluated in a highly controlled manner. To this end, we have bioprinted endothelialized vascular beds and perfused these constructs with metastatic mammary gland cells under physiological flow rates. By pairing these in vitro devices with an advanced computational flow model, we found that the bioprinted analog was readily capable of evaluating the accuracy and integrated complexity of a computational flow model, while also highlighting the discrete contribution of hydrodynamics in vascular colonization. This intersection of these two technologies, bioprinting and computational simulation, is a key demonstration in the establishment of an experimentation pipeline for the understanding of complex biophysical events.
A dynamic Variational Multiscale Method (Hughes et al. 1998) is developed by leveraging the Germano procedure from classical Large-eddy Simulations (LES). The similarity between the classical and variational approaches is analyzed in the context of incompressible flow. This analysis leads to a consistent modeling approach for both incompressible and compressible flows, the latter being demonstrated in a priori testing for low-speed attached and separated boundary layers. Similar to the classical LES procedure from which it is derived, the variational dynamic procedure does not guarantee a positive semi-definite coefficient in the general case. However, reproducing the behavior of the classical LES dynamic approach is seen as a necessary first step to develop a VMM that automatically adjusts to the local resolution and flow physics.
Superradiance, in which the collective behavior of emitters can generate enhanced radiative decay, was first predicted by a model, now known as the Dicke model, that contains a collection of two-level systems (the emitters) all interacting with the same photonic mode. In this article, we extend the original Dicke model to elucidate the influence of nuclear motion on superradiant emission. Our dynamical simulations of the combined electronic, nuclear, and photonic system reveal a new time scale attributed to the population leakage of the dark, subradiant states. Furthermore, this dark-state emission pathway can be controlled by tuning the nuclear potential energy landscape. These findings impact how superradiant states and molecular degrees of freedom can be leveraged and utilized in quantum optical systems. Published by the American Physical Society 2025
Here, in this work, expanding cylinder experiments are performed using a gas gun on the aluminum 6061 alloy to understand the dynamic deformation and failure behavior of this material in the presence of surface defects in the form of scratches. The experiments are then modeled with the modified Tepla (for tensile plasticity) model, which is equipped with the damage yield function of Gurson [“Continuum theory of ductile rupture by void nucleation and growth: Part I—Yield criteria and flow rules for porous ductile media,” (1977)] and Tvergaard and Needleman [Acta Metall. 32(1), 157–169 (1984)], over-stress viscosity [F. Addessio and J. Johnson, J. Appl. Phys. 74(3), 1640–1648 (1993)], and the underlying viscoplastic model of Preston et al. [J. Appl. Phys. 93(1), 211–220 (2003)] (Preston–Tonks–Wallace model) for the dense material. The Tepla model parameters were calibrated using the Bayesian approach against experimentally measured free-surface velocity data obtained from plate impact experiments on the aluminum 6061 alloy. The evolution of radial expansion and shapes of the deformed cylinder from the simulations on a pristine cylinder wall and a cylinder wall with a longitudinal surface defect shows excellent agreement with the corresponding experimental results. Finally, the role of surface defect size and cylinder wall geometry on the dynamic strength of the cylinder wall was also investigated through a series of simulations.
This work presents an approach to modeling the dynamic spreading and curing behavior of thermosets in nozzle-based extrusions. Thermosets cover a wide range of materials, some of which permit low-temperature processing with subsequent high-temperature and high-strength working properties. Extruding thermosets may overcome the limited working temperatures and strengths of conventional thermoplastic materials used in additive manufacturing. This project aims to produce technology for the fabrication of thermoset-based structures leveraging advances made in nozzle-based extrusion, such as fused deposition modeling (FDM), material jetting, and direct writing. Understanding the synergistic interactions between spreading and fast curing of extruded thermosetting materials will provide essential insights for applications that require accurate dimensional controls, such as additive manufacturing [1], [2] and centrifugal coating/forming [3]. Two types of thermally curing thermosets -- one being a soft silicone (Ecoflex 0050) and the other being a toughened epoxy (G/Flex) -- served as the test materials in this work to obtain models for cure kinetics and viscosity. The developed models align with extensive measurements made with differential scanning calorimetry (DSC) and rheology. DSC monitors the change in the heat of reaction, which reflects the rate and degree of cure at different crosslinking stages. Rheology measures the change in complex viscosity, shear moduli, yield stress, and other properties dictated by chemical composition. By combining DSC and rheological measurements, it is possible to establish a set of models profiling the cure kinetics and chemorheology without prior knowledge of chemical composition, which is usually necessary for sophisticated mechanistic modeling. In this work, we conducted both isothermal and dynamic measurements with both DSC and rheology. With the developed models, numerical simulations yielded predictions of diameter and height of droplets, along with width and height of extruded lines cured at varied temperatures. Experimental results carried out on a goniometric platform and a nozzle-based 3D printer showed agreement with the numerical simulations. Finally, this presentation will show how the models are adaptable to the planning of tool paths and designs in additive manufacturing.