Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic models”

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 1,585 records · Page 88

Atomistic Investigation of Plastic Deformation and Dislocation Motion in Uranium Mononitride

Uranium mononitride (UN) is a promising advanced nuclear fuel due to its high thermal conductivity and high fissile density. However, many aspects of its mechanical behavior, particularly at reactor-relevant conditions, remain unclear. In this study, molecular dynamics (MD) simulations were employed to investigate the deformation behavior and dislocation motion in UN. We found that the Kocevski potential predicts the principal slip system as $\frac{1}{2}$ $\langle110\rangle${110}, aligning with experimental data. On the other hand, the Tseplyaev potential predicts slip to primarily occur on $\frac{1}{2}$ $\langle110\rangle${111}. MD simulations of stress–strain behavior were used to estimate the nanoindentation hardness, revealing that the Kocevski potential accurately predicts hardness even though it fails to model dynamic plasticity. Complete dislocation mobility functions have been fitted for the edge and screw dislocations in both the thermally activated and phonon-drag regimes. The 300 K linear mobility of the edge dislocation using the Tseplyaev potential was found to be 817 Pa -1 ·s -1 , whereas that of the screw dislocation using the Kocevski potential was found to be 4546 Pa -1 ·s -1 . At intermediate stresses, we observed that the subsonic steady-state motion of the edge dislocation in UN is intermittently interrupted by velocity jumps, reaching the average sound velocity. Finally, the threshold Schmid stress is calculated as 179–197 MPa, which gives an upper-limit estimate of the uniaxial yield stress of polycrystalline UN of 548–603 MPa. These findings, including the fitted dislocation mobility function, provide essential input for future plasticity and dislocation dynamics models of nuclear fuels.

36 MATERIALS SCIENCE↗

Causal hydrodynamic fluctuations in non-static and inhomogeneous backgrounds

Highlights: •The fluctuation–dissipation relation is modified by inhomogeneous backgrounds. •The constitutive equations are restricted by causality and fluctuation properties. •The structure of constitutive equation gives a natural cutoff scale of the theory. •The memory function with tensor structures in inhomogeneous backgrounds is obtained. -- Abstract: To integrate hydrodynamic fluctuations, namely thermal fluctuations of hydrodynamics, into dynamical models of high-energy nuclear collisions based on relativistic hydrodynamics, the property of the hydrodynamic fluctuations given by the fluctuation–dissipation relation should be carefully investigated. The fluctuation–dissipation relation for causal dissipative hydrodynamics with the finite relaxation time is naturally given in the integral form of the constitutive equation by the linear-response theory. While, the differential form of the constitutive equation is commonly used in analytic investigations and dynamical calculations for practical reasons. We give the fluctuation–dissipation relation for the general linear-response differential form and discuss the restrictions to the structure of the differential form, which comes from the causality and the positive semi-definiteness of the noise autocorrelation, and also the relation of those restrictions to the cutoff scale of the hydrodynamic fluctuations. We also give the fluctuation–dissipation relation for the integral form in non-static and inhomogeneous background by introducing new tensors, the pathline projectors. We find new modification terms to the fluctuation–dissipation relation for the differential form in non-static and inhomogeneous background which are particularly important in dynamical models to describe rapidly expanding systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Time Dependence of Collision Probabilities During Satellite Conjunctions

The NASA Conjunction Assessment Risk Analysis (CARA) team has recently implemented updated software to calculate the probability of collision (P (sub c)) for Earth-orbiting satellites. The algorithm can employ complex dynamical models for orbital motion, and account for the effects of non-linear trajectories as well as both position and velocity uncertainties. This “3D P (sub c)” method entails computing a 3-dimensional numerical integral for each estimated probability. Our analysis indicates that the 3D method provides several new insights over the traditional “2D P (sub c)” method, even when approximating the orbital motion using the relatively simple Keplerian two-body dynamical model. First, the formulation provides the means to estimate variations in the time derivative of the collision probability, or the probability rate, R (sub c). For close-proximity satellites, such as those orbiting in formations or clusters, R (sub c) variations can show multiple peaks that repeat or blend with one another, providing insight into the ongoing temporal distribution of risk. For single, isolated conjunctions, R (sub c) analysis provides the means to identify and bound the times of peak collision risk. Additionally, analysis of multiple actual archived conjunctions demonstrates that the commonly used “2D P (sub c)” approximation can occasionally provide inaccurate estimates. These include cases in which the 2D method yields negligibly small probabilities (e.g., P (sub c)) is greater than 10 (sup -10)), but the 3D estimates are sufficiently large to prompt increased monitoring or collision mitigation (e.g., P (sub c) is greater than or equal to 10 (sup -5)). Finally, the archive analysis indicates that a relatively efficient calculation can be used to identify which conjunctions will have negligibly small probabilities. This small-P (sub c) screening test can significantly speed the overall risk analysis computation for large numbers of conjunctions.

Hall, Doyle T.↗

Hierarchical optimization for neutron scattering problems

We present a scalable optimization method for neutron scattering problems that determines confidence regions of simulation parameters in lattice dynamics models used to fit neutron scattering data for crystalline solids. The method uses physics-based hierarchical dimension reduction in both the computational simulation domain and the parameter space. We demonstrate for silicon that after a few iterations the method converges to parameters values (interatomic force-constants) computed with density functional theory simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Track-train dynamics

The purpose of this paper is to outline one of the serious dynamic problems involving railroad cars and to propose an orderly approach to the definition of a dynamic model through analysis and testing that can be used to investigate various fixes. The problem in question is the hunting phenomenon involving lateral instability of the car suspension (trucks). This phenomenon has been recognized as a major contributing factor in a significant number of derailments causing millions of dollars in losses. This paper expounds on the mechanics and characteristics of hunting. The planned solution is outlined as a building block approach directed at formulating a test verified dynamic model of an open hopper freight car-truck assembly. Results obtained to date are presented together with conclusions related to the program.

Abbott, P. W.↗

Comparison of Surrogate Models with Physical Models for Dynamic Probabilistic Risk Analysis Using the RAVEN Code

Dynamic probabilistic risk analysis (DPRA) has emerged in the nuclear industry as an effective tool to identify points of weakness in nuclear reactors and create a comprehensive approach to quantifying risk. While DPRA allows for detailed analysis, it usually requires thousands of runs due to the necessity of exploring low probability, high consequence events, which can be computationally prohibitively difficult to assess. Risk Analysis Virtual Environment (RAVEN) is a software package being developed under the Nuclear Energy Advanced Modeling and Simulation program at the Idaho National Laboratory. It is able to manage complex control logic to drive a RELAP-7 simulation and to handle sampling strategies including Monte Carlo, Latin Hypercube, Grid and Adaptive samplers. By performing a thorough investigation of the possibility space of an accident scenario with RELAP-7, RAVEN is able to identify consequences of accident scenarios that have not been encountered. The most computational resource intensive step of a RAVEN analysis is running the RELAP-7 code modeling the physical system behavior. Using surrogate models to represent the outcomes of more complex simulations rather than the full model can greatly reduce the computation effort required. In situations where the uncertainty distributions of a model are unknown, or when there are multiple sets of uncertainty distributions, DPRA generally requires rerunning the physical model for each set of uncertainty distributions used. By taking an unbiased sample of the physical model and creating a surrogate model from the results, the surrogate model can be sampled under a different set of uncertainty distributions to provide equivalent results. This paper presents an example of this approach. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Control and dynamics of a flexible spacecraft during stationkeeping maneuvers

A case study of a spacecraft having flexible solar arrays is presented. A stationkeeping attitude control mode using both earth and rate gyro reference signals and a flexible vehicle dynamics modeling and implementation is discussed. The control system is designed to achieve both pointing accuracy and structural mode stability during stationkeeping maneuvers. Reduction of structural mode interactions over the entire mode duration is presented. The control mode using a discrete time observer structure is described to show the convergence of the spacecraft attitude transients during Delta-V thrusting maneuvers without preloading thrusting bias to the onboard control processor. The simulation performance using the three axis, body stabilized nonlinear dynamics is provided. The details of a five body dynamics model are discussed. The spacecraft is modeled as a central rigid body having cantilevered flexible antennas, a pair of flexible articulated solar arrays, and to gimballed momentum wheels. The vehicle is free to undergo unrestricted rotations and translations relative to inertial space. A direct implementation of the equations of motion is compared to an indirect implementation that uses a symbolic manipulation software to generate rigid body equations.

Liu, D.↗

A Method to Capture Macroslip at Bolted Interfaces [PowerPoint]

Relative motion at bolted connections can occur for large shock loads as the internal shear force in the bolted connection overcomes the frictional resistive force. This macroslip in a structure dissipates energy and reduces the response of the components above the bolted connection. There is a need to be able to capture macroslip behavior in a structural dynamics model. A linear model and many nonlinear models are not able to predict marcoslip effectively. The proposed method to capture macroslip is to use the multi-body dynamics code ADAMS to model joints with 3-D contact at the bolted interfaces. This model includes both static and dynamic friction. The joints are preloaded and the pinning effect when a bolt shank impacts a through hole inside diameter is captured. Substructure representations of the components are included to account for component flexibility and dynamics. This method was applied to a simplified model of an aerospace structure and validation experiments were performed to test the adequacy of the method.

42 ENGINEERING↗

System dynamics for predicting the impact of traffic noise on cardiovascular mortality in Madrid

Highlights: • Short-term impact of traffic noise has changing patterns over time. • System dynamics models help understand the evolution of noise impact. • In Madrid, noise impact on mortality exhibits an epidemic behaviour since 1980. • Growing use of electric vehicles would stop the epidemic no sooner than 2040. • Major changes in urban life are needed to reduce noise impact in large cities. The impact of road traffic noise on health is an important issue worldwide. In Madrid (Spain), the existing high daily noise levels all year round have been associated with short-term cardiovascular mortality. A system dynamics model is constructed to estimate the evolution of the traffic noise impact on the population older than 65 in Madrid over the upcoming years. This strategy allows to make predictions based on yearly variations in traffic intensity, population at risk, motor vehicle-related technological advances and lifestyles, and urban planning policies. Madrid harbours a persisting noise epidemic since 1980, which would be eradicated in 2040 at the earliest if significant action for noise abatement were taken.

54 ENVIRONMENTAL SCIENCES↗

COSMIC REIONIZATION ON COMPUTERS. ULTRAVIOLET CONTINUUM SLOPES AND DUST OPACITIES IN HIGH REDSHIFT GALAXIES

We compare the properties of stellar populations of model galaxies from the Cosmic Reionization On Computers (CROC) project with the exiting ultraviolet (UV) and IR data. Since CROC simulations do not follow cosmic dust directly, we adopt two variants of the dust-follows-metals ansatz to populate model galaxies with dust. Using the dust radiative transfer code Hyperion, we compute synthetic stellar spectra, UV continuum slopes, and IR fluxes for simulated galaxies. We find that the simulation results generally match observational measurements, but, perhaps, not in full detail. The differences seem to indicate that our adopted dust-follows-metals ansatzes are not fully sufficient. While the discrepancies with the exiting data are marginal, the future James Webb Space Telescope (JWST) data will be of much higher precision, rendering highly significant any tentative difference between theory and observations. It is, therefore, likely, that in order to fully utilize the precision of JWST observations, fully dynamical modeling of dust formation, evolution, and destruction may be required.

79 ASTRONOMY AND ASTROPHYSICS↗

Dynamic properties of micro-particles in ultrasonic transportation using phase-controllable standing waves

Ultrasonic manipulation has become an attractive method for surface-sensitive objects in micro-technology. Related phenomena, such as radiation force, multiple scattering, and acoustic streaming, have been widely studied. However, in current studies, the behavior of micro-particles in potential force fields is always analyzed in a quasi-static manner. We developed a dynamic model of a dilute micro-particle in the commonly used two-dimensional ultrasonic manipulation system to provide a systemic and quantitative analysis of the transient properties of particle movement. In this model, the acoustic streaming and hydrodynamic forces, omitted in previous work, were both considered. The trajectory of a spherical silica particle with different initial conditions was derived by numerically solving the established nonlinear differential integral equation system, which was then validated experimentally. The envelope of the experimental data on the x-axis showed good agreement with the theoretical calculation, and the greater influence on the y-axis of the deviation between the actual sound field and the ideal distribution employed in our dynamic model could account for the differences in displacement in that direction. Finally, the influence of particle size on its movement and the effect of acoustic streaming on calculating the hydrodynamic forces for an isolated particle with motion relative to the fluid were analyzed theoretically. It was found that the ultrasonic manipulation system will translate from an under-damped system to an over-damped system with a decrease in particle size and the micro-scale acoustic streaming velocity was negligible when calculating the hydrodynamic forces on the particle in the ultrasonic manipulation system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Numerical investigation of electromagnetic pulse welded interfaces between dissimilar metals

The quality of an electromagnetic pulse weld is strongly affected by a variety of interacting physical phenomena including large plastic deformation, materials mixing, localised heating and rapid cooling, and possible localised melting and interfacial cracks and void. In the present study, a thermomechanically coupled dynamic model has been developed to quantitatively resolve the high speed impact joining interface characteristics as well as the process induced interface temperature evolution, defect formation and possible microstructural composition variations. Reasonably good agreement has been obtained between the predicted results and the experimental measurements in terms of interfacial morphology characteristics. The modelling framework is expected to provide further understanding of the hierarchical interfacial features of the non-equilibrium material joining process and the weld formation mechanisms.

Xu, W.↗

Modeling Longitudinal Dynamics in the Fermilab Booster Synchrotron

The PIP-II project will replace the existing 400 MeV linac with a new, CW-capable, 800 MeV superconducting one. With respect to current operations, a 50% increase in beam intensity in the rapid cycling Booster synchrotron is expected. Booster batches are combined in the Recycler ring; this process limits the allowed longitudinal emittance of the extracted Booster beam. To suppress eddy currents, the Booster has no beam pipe; magnets are evacuated, exposing the beam to core laminations and this has a substantial impact on the longitudinal impedance. Noticeable longitudinal emittance growth is already observed at transition crossing. Operation at higher intensity will likely necessitate mitigation measures. We describe systematic efforts to construct a predictive model for current operating conditions. A longitudinal only code including a laminated wall impedance model, space charge effects, and feedback loops is developed. Parameter validation is performed using detailed measurements of relevant beam, rf and control parameters. An attempt is made to benchmark the code at operationally favorable machine settings.

43 PARTICLE ACCELERATORS↗

Library of Advanced Models of large-scale PV (LAMP) (Final Technical Report)

A suite of dynamic models of photovoltaic (PV) systems is developed in this project to accurately capture PV system dynamics during different contingencies. The following models are included: (i) quasi-dynamic, (ii) high-fidelity electromagnetic transient (EMT), (iii) advanced, and (iv) baseline. The models have the flexibility to consider different control functionalities at different controller levels. The study’s requirements determine the type of model needed. For example, high-fidelity EMT models are recommended for contingency studies. Advanced models are recommended for online or real-time implementations with resource constraints, whereas quasi-dynamic models are useful for studying the impact of weather variations.

14 SOLAR ENERGY↗

Anomalous thermal conductivity of monolayer boron nitride

In this paper, we use nonequilibrium molecular dynamics modeling to investigate the thermal properties of monolayer hexagonal boron nitride nanoribbons under uniaxial strain along their longitudinal axis. Our simulations predict that hexagonal boron nitride shows an anomalous thermal response to the applied uniaxial strain. Contrary to three dimensional materials, under uniaxial stretching, the thermal conductivity of boron nitride nanoribbons first increases rather than decreasing until it reaches its peak value and then starts decreasing. Under compressive strain, the thermal conductivity of monolayer boron nitride ribbons monolithically reduces rather than increasing. We use phonon spectrum and dispersion curves to investigate the mechanism responsible for the unexpected behavior. Our molecular dynamics modeling and density functional theory results show that application of longitudinal tensile strain leads to the reduction of the group velocities of longitudinal and transverse acoustic modes. Such a phonon softening mechanism acts to reduce the thermal conductivity of the nanoribbons. On the other hand, a significant increase in the group velocity (stiffening) of the flexural acoustic modes is observed, which counteracts the phonon softening effects of the longitudinal and transverse modes. The total thermal conductivity of the ribbons is a result of competition between these two mechanisms. At low tensile strain, the stiffening mechanism overcomes the softening mechanism which leads to an increase in the thermal conductivity. At higher tensile strain, the softening mechanism supersedes the stiffening and the thermal conductivity slightly reduces. Our simulations show that the decrease in the thermal conductivity under compressive strain is attributed to the formation of buckling defects which reduces the phonon mean free path.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Lattice Anharmonicity and Thermal Conductivity from Compressive Sensing of First-Principles Calculations

First-principles prediction of lattice thermal conductivity K L of strongly anharmonic crystals is a long-standing challenge in solid state physics. Using recent advances in information science, we propose a systematic and rigorous approach to this problem, compressive sensing lattice dynamics (CSLD). Compressive sensing is used to select the physically important terms in the lattice dynamics model and determine their values in one shot. Non-intuitively, high accuracy is achieved when the model is trained on first-principles forces in quasi-random atomic configurations. The method is demonstrated for Si, NaCl, and Cu 12 Sb 4 S 13 , an earth-abundant thermoelectric with strong phononphonon interactions that limit the room-temperature K L to values near the amorphous limit.

36 MATERIALS SCIENCE↗