Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “long-time integration”

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 19 records

Finite difference schemes for long-time integration

Finite difference schemes for the evaluation of first and second derivatives are presented. These second order compact schemes were designed for long-time integration of evolution equations by solving a quadratic constrained minimization problem. The quadratic cost function measures the global truncation error while taking into account the initial data. The resulting schemes are applicable for integration times fourfold, or more, longer than similar previously studied schemes. A similar approach was used to obtain improved integration schemes.

Haras, Zigo↗

Environmental effects on the torsion and flexure fatigue properties of graphite-epoxy composites

Effects of torsional and flexural fatigue on the long-time integrity of advanced graphite-epoxy structural composites have been investigated. Reversed (R = - 1) and unidirectional (R = 0) torsional fatigue tests were made on unidirectional, cross-plied and woven graphite fiber materials in air and water at room temperature and at 74 C. Flexural fatigue tests (4-point bending) with R = - 1 were run in air and water at room temperature, and at R = 0 in air. Results show that in torsional cycling, both water environment and higher test temperature contribute to significant degradation of torsional stiffness. The degradation of stiffness from torsional stress cycling was observed to be much greater with R = - 1 than with simple R = 0 cycling. The effect of environment is also greater in the fully reversed cycling.

Sumsion, H. T.↗

The effects of rotation on initially anisotropic homogeneous flows

Rotation of initially anisotropic homogeneous flows is studied using a model spectral tensor. It is shown that the anisotropy changes because of the influence of rotation through phase scrambling. Phase scrambling causes the Reynolds stresses to develop with damped oscillations. The final Reynolds stress anisotropy is found to be proportional to the initial structural tensor anisotropy. Closure models for the rapid pressure strain terms should reflect this change in anisotropy, and should drive the anisotropy to reach its final predicted state. Finally, it is shown that long-time integration using direct numerical simulations should be treated with care because phase scrambling effects on a discrete wave space can cause loss of resolution when time becomes large.

Mansour, Nagi N.↗

Environmental Effects on Graphite-Epoxy Fatigue Properties

Effects of torsional and flexural fatigue on the long-time Integrity of advanced graphite-epoxy structural composites have been investigated. Torsional fatigue tests were run at stress ratios of R = 0 (zero to maximum, repeated) and R = -1 (zero mean stress) on unidirectional, angleply, and woven graphite fiber materials in air and water at room temperature and at 74 C. Flexural fatigue tests (four-point bending) with R = -1 were run in air and water at room temperature, and with R = 0 in air. Results show that, in torsional cycling, both water environment and higher test temperature contribute to significant degradation of torsional stiffness. The degradation of stiffness from torsional stress cycling was observed to be much greater with R = -1 than with simple R = 0 cycling. The effect of environment also is greater in the fully reversed cycling. Flexural fatigue results on +/- 30 deg material show a large fatigue effect, with fatigue limits of less than 50% and 30% of the static failure strength for specimens tested under stress ratios of R = 0 and R = -1, respectively. Compliance measurements indicate that the final failures are preceded by damage initiation and accumulation, which begins at about 1% of the specimen life.

Sumsion, H. T.↗

Data Assimilation with Machine Learning Surrogate Models: A Case Study with FourCastNet

Modern data-driven surrogate models for weather forecasting provide accurate short-term predictions but inaccurate and nonphysical long-term forecasts. This paper investigates online weather prediction using machine learning surrogates supplemented with partial and noisy observations. We empirically demonstrate and theoretically justify that, despite the long-time instability of the surrogates and the sparsity of the observations, filtering estimates can remain accurate in the long-time horizon. As a case study, we integrate the Fourier Forecasting Neural Network (FourCastNet), a weather surrogate model, within a variational data assimilation framework using partial, noisy ERA5 global reanalysis data from the European Centre for Medium-Range Weather Forecasts (ECMWF). Here, our results show that filtering estimates remain accurate over a year-long assimilation window and provide effective initial conditions for forecasting tasks, including extreme event prediction.

Data assimilation↗

Data Assimilation with Machine Learning Surrogate Models: A Case Study with FourCastNet

Modern data-driven surrogate models for weather forecasting provide accurate short-term predictions but inaccurate and nonphysical long-term forecasts. This paper investigates online weather prediction using machine learning surrogates supplemented with partial and noisy observations. We empirically demonstrate and theoretically justify that, despite the long-time instability of the surrogates and the sparsity of the observations, filtering estimates can remain accurate in the long-time horizon. As a case study, we integrate FourCastNet, a weather surrogate model, within a variational data assimilation framework using partial, noisy ERA5 data. Our results show that filtering estimates remain accurate over a year-long assimilation window and provide effective initial conditions for forecasting tasks, including extreme event prediction.

Adrian, Melissa [Univ. of Chicago, IL (United Stat↗

Quantum chaos, integrability, and late times in the Krylov basis

Quantum chaotic systems are conjectured to display a spectrum whose fine-grained features (gaps and correlations) are well described by random matrix theory (RMT). We propose and develop a complementary version of this conjecture: quantum chaotic systems display a Lanczos spectrum whose local means and covariances are well described by RMT. To support this proposal, we first demonstrate its validity in examples of chaotic and integrable systems. We then show that for Haar-random initial states in RMTs the mean and covariance of the Lanczos spectrum suffice to produce the full long-time behavior of general survival probabilities including the spectral form factor, as well as the spread complexity. In addition, for initial states with continuous overlap with energy eigenstates, we analytically find the long-time averages of the probabilities of Krylov basis elements in terms of the mean Lanczos spectrum. This analysis suggests a notion of eigenstate complexity, the statistics of which differentiate integrable systems and classes of quantum chaos. Lastly, we clarify the relation between spread complexity and the universality classes of RMT by exploring various values of the Dyson index and Poisson distributed spectra.

Combinatorics↗

Long-time predictions in nonlinear dynamics

It is known that nonintegrable dynamical systems do not allow precise predictions concerning their behavior for arbitrary long times. The available series solutions are not uniformly convergent according to Poincare's theorem and numerical integrations lose their meaningfulness after the elapse of arbitrary long times. Two approaches are the use of existing global integrals and statistical methods. This paper presents a generalized method along the first approach. As examples long-time predictions in the classical gravitational satellite and planetary problems are treated.

Szebehely, V.↗

Quasiperiodic behavior in beam-driven strong Langmuir turbulence

The evolution of unmagnetized beam-driven strong Langmuir turbulence is studied in two dimensions by numerically integrating the Zakharov equations for systems pumped by monochromatic and broadband negative-damping drivers with nonzero central wavenumber. Long-time statistically steady states are reached for which the dependence of the evolution on the driver wavenumber, growth rate, and bandwidth is examined in detail. For monochromatic drivers, a quasiperiodic cycle is found to develop if the driver wavenumber is sufficiently large. The characteristic frequency of the quasiperiodic cycle and the average system energy are both approximately proportional to the growth rate. Broadening of the driver in wavenumber tends to degrade the system-wide coherence of the cycle, but its main features appear to survive on the scale of the coherence length of the driver.

Robinson, P. A.↗

Replica symmetry breaking in spin glasses in the replica-free Keldysh formalism

We show that the algebra of Parisi ultrametric matrices is recovered by the real-time, replica-free, Dyson-Keldysh equations of infinite-range quantum spin glasses in the late time glassy limit. This connects to earlier results on classical and quantum systems showing how ultrametricity emerges from the persistent slow aging dynamics of the glass phase. The stationary spin glass state thereby spontaneously breaks thermal symmetry, or the Kubo-Martin-Schwinger relation of a state in global thermal equilibrium. We describe the Keldysh path integral of the infinite-range Ising model in transverse and longitudinal fields, and in the context of the Landau expansion of the action functional, show how the long-time limit connects to the full replica symmetry breaking obtained in the equilibrium formalism. We also illustrate our formalism by applying it to the spherical quantum p p -spin model, which only exhibits one-step replica symmetry breaking.

Lang, Johannes (ORCID:0000000283533392)↗

On the buoyancy-driven motion of a drop towards a rigid surface or a deformable interface

The deformation of a viscous drop, driven by buoyancy toward a solid surface or a deformable interface, is analyzed in the asymptotic limit of small Bond number, for which the deformation becomes important only when the drop is close to the solid surface or interface. Lubrication theory is used to describe the flow in the thin gap between the drop and the solid surface or interface, and boundary-integral theory is used in the fluid phases on either side of the gap. The evolution of the drop shape is traced from a relatively undeformed state until a dimple is formed and a long-time quasi-steady-state pattern is established. A wide range of drop to suspending phase viscosity ratios is examined. It is shown that a dimple is always formed, independently of the viscosity ratio, and that the long-time thinning rates take simple forms as inverse fractional powers of time.

Yiantsios, Stergios G.↗

Predicting nonequilibrium Green’s function dynamics and photoemission spectra via nonlinear integral operator learning

Understanding the dynamics of nonequilibrium quantum many-body systems is an important research topic in a wide range of fields across condensed matter physics, quantum optics, and high-energy physics. However, numerical studies of large-scale nonequilibrium phenomena in realistic materials face serious challenges due to intrinsic high-dimensionality of quantum many-body problems and the absence of time-invariance. The nonequilibrium properties of many-body systems can be described by the dynamics of the correlator, or the Green's function of the system, whose time evolution is given by a high-dimensional system of integro-differential equations, known as the Kadanoff–Baym equations (KBEs). The time-convolution term in KBEs, which needs to be recalculated at each time step, makes it difficult to perform long-time numerical simulation. In this paper, we develop an operator-learning framework based on recurrent neural networks (RNNs) to address this challenge. We utilize RNNs to learn the nonlinear mapping between Green's functions and convolution integrals in KBEs. By using the learned operators as a surrogate model in the KBE solver, we obtain a general machine-learning scheme for predicting the dynamics of nonequilibrium Green's functions. Besides significant savings per each time step, the new methodology reduces the temporal computational complexity from $O(N_t^3)$ to $O(N_t)$ where N t is the number of steps taken in a simulation, thereby making it possible to study large many-body problems which are currently infeasible with conventional KBE solvers. Through various numerical examples, we demonstrate the effectiveness of the operator-learning based approach in providing accurate predictions of physical observables such as the reduced density matrix and time-resolved photoemission spectra. Moreover, our framework exhibits clear numerical convergence and can be easily parallelized, thereby facilitating many possible further developments and applications.

97 MATHEMATICS AND COMPUTING↗

Mitigation for roof alterations to building 06-cp-65 at the area 6 control point, nevada national security site, nye county, nevada

Building 06-CP-65 has been determined to be a contributing element to the Area 6 Control Point Historic District (O’Neill et al. 2021; Reed 2022). It contributes to the significance of the historic district under Criterion A as one of the principal buildings that supported timing and firing operations for nuclear testing on the NNSS from 1966 to 1992. As such, the building served as a major warehouse with office space within the district. It was used by both REECo, a long-time general contractor at the NNSS, as well as EG&G, which provided technical support to the national laboratories and the DOD. The building served as an important staging area and electrical power supply point for the NNSS diagnostic trailer fleet and its unique location immediately along Mercury Highway allowed easy accessibility to the forward areas of the NNSS. The building also contributes to the historic district under Criterion C as it is one of the large, unadorned, precast concrete buildings at the Control Point. These types of buildings were the prominent elements of the compound that convey the district’s overall utilitarian, military-industrial character. Building 06-CP-65 retains all aspects of integrity to a high degree and easily conveys its significance as a warehouse that supported Control Point operations. Building 06-CP-65 is not recommended eligible for listing in the NRHP as an individual resource. While it served an important support function as part of the Control Point Historic District, an archival and literature review did not reveal information linking it to any specific test, series of tests, programs, or for any specific role on the NNSS other than as a warehouse (Criterion A). It has no direct association with any important individual (Criterion B). It also is not architecturally significant in its own right beyond reflecting the overall aesthetic of the Control Point Historic District (Criterion C), and it does not have potential to yield information important to the history of nuclear testing beyond what can be learned from historic texts, drawings, and other documents (Criterion D).

54 ENVIRONMENTAL SCIENCES↗

Some considerations on velocity vector accuracy in dust trajectory analysis

The relative contributions of comets and asteroids to the reservoir of dust in the interplanetary medium is not known. There are direct observations of dust released from comets and there is evidence to associate the IRAS dust bands with possible collisions of asteroids in the main belt. A means towards sorting out the parent sources has been proposed in the establishment of a dust collector in orbit about the Earth. The purpose of such a facility would be to collect not only cosmic dust particles intact but also the state vectors, as they arrive at the detector, the idea being that one may combine analytical laboratory analysis of the physics and chemistry of the captured particles with orbital data in order to help distinguish between bodies and identify parent bodies. The theoretical study of dust particle orbits in the solar system takes on greatly more importance if we use collected trajectory data. The orbital motion of dust when radiation and forces alone are acting is well understood. When gravitational forces due to the planets are included, the motion can become quite complex. In order to characterize the orbits of particles as they crossed the Earth's orbits, a study of the long-time dust orbital evolution was undertaken. We have considered various parameters associated with these dust orbits to see if one may in a general way discriminate between particles evolved from comets and asteroids. We proceed in this study as we have done previously. That is, we considered the dust particles as ideal black bodies, of density 1 gm/cc, spherical, with radii 10-100 microns. Particles of this size are affected by radiation forces, photon pressure, and Poynting-Robertson drag. Account was also taken of solar wind drag, which amounts to about 30 percent of the Poynting-Robertson drag negligible. The gravitational forces due to the planets are included, unlike in our previous study; the planetary orbits are those of true n-body interaction so that the possibility of secular resonance is included. Our method was to calculate explicitly by a numerical procedure the orbits of dust particles after they left their parent bodies. The motion is determined numerically with the implicit Runge-Kutta integrator using Gauss-Radau spacings.

Jackson, A. A.↗

Turbulent transport measurements with a laser Doppler velocimeter.

The power spectrum of phototube current from a laser Doppler velocimeter operating in the heterodyne mode has been computed. The spectral width and shape predicted by the theory are in agreement with experiment. For normal operating parameters the time-average spectrum contains information only for times shorter than the Lagrangian-integral time scale of the turbulence. To examine the long-time behavior, one must use either extremely small scattering angles, much-longer-wavelength radiation, or a different mode of signal analysis, e.g., FM detection.

Edwards, R. V.↗

Life prediction modeling based on strainrange partitioning

Strainrange partitioning (SRP) is an integrated low-cycle-fatigue life predicting system. It was created specifically for calculating cyclic crack initiation life under severe high-temperature fatigue conditions. The key feature of the SRP system is its recognition of the interacting mechanisms of cyclic inelastic deformation that govern cyclic life at high temperatures. The SRP system bridges the gap between the mechanistic level of understanding that breeds new and better materials and the phenomenological level wherein workable engineering life prediction methods are in great demand. The system was recently expanded to address engineering fatigue problems in the low-strain, long-life, nominally elastic regime. This breakthrough, along with other advances in material behavior and testing technology, has permitted the system to also encompass low-strain thermomechanical loading conditions. Other important refinements of the originally proposed method include procedures for dealing with life-reducing effects of multiaxial loading, ratcheting, mean stresses, nonrepetitive (cumulative loading) loading, and environmental and long-time exposure. Procedure were also developed for partitioning creep and plastic strain and for estimating strainrange versus life relations from tensile and creep rupture properties. Each of the important engineering features of the SRP system are discussed and examples shown of how they help toward predicting high-temperature fatigue life under practical, although complex, loading conditions.

Halford, Gary R.↗

Computation of vortex sheet roll-up in the Trefftz plane

Two vortex-sheet evolution problems arising in aerodynamics are studied numerically. The approach is based on desingularizing the Cauchy principal value integral which defines the sheet's velocity. Numerical evidence is presented which indicates that the approach converges with respect to refinement in the mesh-size and the smoothing parameter. For elliptic loading, the computed roll-up is in good agreement with Kaden's asymptoic spiral at early times. Some aspects of the solution's instability to short-wavelength perturbations, for a small value of the smoothing parameter, are inferred by comparing calculations performed with different levels of computer round-off error. The tip vortices' deformation, due to their mutual interaction, is shown in a long-time calculation. Computations for a simulated fuselage-flap configuration show a complicated process of roll-up, deformation and interaction involving the tip vortex and the inboard neighboring vortices.

Krasny, Robert↗

The asymptotic state of rotating homogeneous turbulence at high Reynolds numbers

The long-time, asymptotic state of rotating homogeneous turbulence at high Reynolds numbers has been examined using large-eddy simulation of the incompressible Navier-Stokes equations. The simulations were carried out using 128 x 128 x 512 collocation points in a computational domain that is four times longer along the rotation axis than in the other directions. Subgrid-scale motions in the simulations were parameterized using a spectral eddy viscosity modified for system rotation. Simulation results show that in the asymptotic state the turbulence kinetic energy undergoes a power-law decay with an exponent which is independent of rotation rate, depending only on the low-wavenumber form of the initial energy spectrum. Integral lengthscale growth in the simulations is also characterized by power-law growth; the correlation length of transverse velocities exhibiting much more rapid growth than observed in non-rotating turbulence.

Squires, Kyle D.↗