Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical relativity”

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 523 records · Page 29

Simulating a numerical UV completion of quartic Galileons

The Galileon theory is a prototypical effective field theory that incorporates the Vainshtein screening mechanism—a feature that arises in some extensions of general relativity, such as massive gravity. The Vainshtein effect requires that the theory contain higher order derivative interactions, which results in Galileons, and theories like them, failing to be technically well posed. While this is not a fundamental issue when the theory is correctly treated as an effective field theory, it nevertheless poses significant practical problems when numerically simulating this model. These problems can be tamed using a number of different approaches: introducing an active low-pass filter and/or constructing a UV completion at the level of the equations of motion, which controls the high momentum modes. These methods have been tested on cubic Galileon interactions, and have been shown to reproduce the correct low-energy behavior. Here we show how the numerical UV-completion method can be applied to quartic Galileon interactions, and present the first simulations of the quartic Galileon model using this technique. We demonstrate that our approach can probe physics in the regime of the effective field theory in which the quartic term dominates, while successfully reproducing the known results for cubic interactions. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

An unconditionally stable, time-implicit algorithm for solving the one-dimensional Vlasov–Poisson system

The development of an implicit, unconditionally stable, numerical method for solving the Vlasov–Poisson system in one dimension using a phase-space grid is presented. The algorithm uses the Crank–Nicolson discretization scheme and operator splitting allowing for direct solution of the finite difference equations. This method exactly conserves particle number, enstrophy and momentum. A variant of the algorithm which does not use splitting also exactly conserves energy but requires the use of iterative solvers. This algorithm has no dissipation and thus fine-scale variations can lead to oscillations and the production of negative values of the distribution function. We find that overall, the effects of negative values of the distribution function are relatively benign. We consider a variety of test cases that have been used extensively in the literature where numerical results can be compared with analytical solutions or growth rates. We examine higher-order differencing and construct higher-order temporal updates using standard composition methods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Ensemble Kalman filter for data assimilation coupled with low-resolution computations techniques applied in fluid dynamics

This paper presents an innovative Reduced-order model (ROM) for merging experimental and simulation data using data assimilation (DA) to estimate the "True" state of a fluid dynamics system, leading to more accurate predictions. Our methodology introduces a novel approach by implementing the ensemble Kalman filter (EnKF) within a reduced-dimensional framework, grounded in a robust theoretical foundation and applied to fluid dynamics. To address the substantial computational demands of DA, the proposed ROM employs low-resolution (LR) techniques to drastically reduce computational costs. This innovative approach involves downsampling datasets for DA computations, followed by an advanced reconstruction technique based on low-cost singular value decomposition (lcSVD). The lcSVD method, a key innovation in this paper, has never been applied to DA before and offers a highly efficient way to enhance resolution with minimal computational resources. Our results demonstrate significant reductions in both computation time and RAM usage through these LR techniques without compromising the accuracy of the estimations. For instance, in a turbulent test case, for a data compression rate of 15.9, the LR approach can achieve a speed-up of 13.7 and a RAM compression of 90.9% while maintaining a low relative root mean square error (RRMSE) of 2.6%, compared to 0.8% in the high-resolution (HR) reference. Furthermore, we highlight the effectiveness of the EnKF in estimating and predicting the state of fluid flow systems based on limited observations and given low-fidelity numerical data. This paper highlights the potential of the proposed DA method in fluid dynamics applications, particularly for improving computational efficiency in CFD and related fields. Its ability to balance accuracy with low computational and memory costs makes it especially suitable for large-scale and real-time applications, such as environmental monitoring or engineering design. This method will be incorporated into ModelFLOWs-app.

Data Assimilation↗

A self-floating oscillating surge wave energy converter

Here this paper proposes a new conceptual configuration for oscillating surge wave energy converters. The concept is a self-floating device consisting of two vertical flaps with a common hinge and an averaged density equal to that of water. The two flaps will oscillate in opposite directions when driven by incident waves. Kinetic energy will be extracted from the relative speed between the two flaps by a power take-off embedded at the hinge. Numerical studies reveal that the self-floating device has a smaller resonant period and optimal power take-off damping compared with a bottom-hinged floating device of the same overall dimensions. As long as the power take-off damping is appropriately tuned, the opposite-rotation mechanism is valid in both regular and irregular waves, and no matter whether the hinge is fixed or connected with mooring lines. With the self-floating feature and opposite-rotation mechanism, the proposed concept can potentially benefit deep ocean applications by exerting a significantly less reaction loading on mooring lines and eliminating the need for a supporting frame to house the power take-off. The influence of various design parameters such as flap dimensions, submerged depth, incident wave direction, and water depth is also analyzed for a comprehensive evaluation of the proposed concept.

16 TIDAL AND WAVE POWER↗

A Feasibility Study on Use of an Accelerometer to Measure the Dynamic Forces in Turning

Metal cutting is a highly dynamic process that generates continuously varying forces. Measurement of these forces is essential to characterizing a cutting process and to fully utilize the machine tool. A force dynamometer is often used to measure these varying forces; however, it is limited by the natural frequency of the sensor and can sometimes be hard to set up on a machine. Therefore, this study aims to estimate the dynamic component of the cutting force using an accelerometer placed directly below the cutting edge. The placement of the sensor helps achieving better signal to noise ratio (S/N). The cutting tool is modelled as a single degree of freedom (SDOF) system relative to the workpiece and dynamic loads are experimentally simulated on it. The measured acceleration is numerically integrated to obtain velocity and deflection, which along with the natural frequency (ω n ) and damping ratio (ζ), is used to predict the dynamic load using the system equation. The dynamic forces tested over a wide range of frequencies, show good agreement with the data from simultaneous dynamometer measurement and the force estimated using the proposed method. The study shows the feasibility of this method in a real cutting scenario and the capability to apply it to a multi DOF system. One major limitation is the inability to capture the static or quasi-static component of the cutting force.

Cutting vibrations↗

Limitations of traditional tools for beyond design basis external hazard PRA

Probabilistic risk assessment (PRA) is being used increasingly by the nuclear industry for safety during normal operations as well as for the protection against external hazards. Computation of total risk in an external hazard PRA is dependent on hazard assessment, fragility assessment, and systems analysis. A systems analysis for propagation of component fragilities is conducted using event and fault trees. The event and fault trees for an actual power plant can be fairly large in size, which imposes computational challenges. Hence, certain assumptions are employed for computational efficiency. These assumptions typically represent the conditions imposed during the design basis (DB) scenario. The traditional PRA tools based on these assumptions are also widely applied to perform risk assessment in the context of beyond design basis (BDB) scenarios. However, some of these assumptions may not be valid for certain BDB scenarios. In addition, the probability of dependent failures also increases in BDB scenarios due to common cause failures (CCF) which usually results from design modifications, human errors, etc. In this manuscript, a simple and a relatively more complex illustrative examples are used to show the limitation of these assumptions in the numerical quantification of risk for the case of BDB conditions. Case studies with CCF events across multiple fault trees are also presented to illustrate the effect of these assumptions when traditional approach is used in BDB risk assessment. It is shown that the assumptions are valid for the case of DB conditions but may lead to excessively conservative risk estimates in the case of BDB conditions. Finally, a Bayesian network based top-down algorithm is proposed as an alternative tool for accurate numerical quantification of total risk in systems analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Defect-Enabled Phase Programming of Transition Metal Dichalcogenide Monolayers

The ability to tune the local electronic transport properties of group VI transition metal dichalcogenide (TMD) monolayers by strain-induced structural phase transformations (“phase programming”) has stimulated much interest in the potential applications of such layers as ultrathin programmable and dynamically switchable nanoelectronics components. In this manuscript, we propose a new approach toward controlling TMD monolayer phases by employing macroscopic in-plane strains to amplify heterogeneous strains arising from tailored, spatially extended defects within the monolayer. Furthermore, the efficacy of our proposed approach is demonstrated via numerical simulations of emerging domains localized around arrays of holes, grain boundaries, and compositional heterointerfaces. Quantitative relations between the macroscopic strains required, spatial resolution of domain patterns, and defect configurations are developed. In particular, the introduction of arrays of holes is identified as the most feasible phase programming route.

2D materials↗

A New Crater Near InSight: Implications for Seismic Impact Detectability on Mars

A new 1.5 m diameter impact crater was discovered on Mars only ~40 km from the InSight lander. Context camera images constrained its formation between 21 February and 6 April 2019; follow-up High Resolution Imaging Science Experiment images resolved the crater. During this time period, three seismic events were identified in InSight data. In this paper, we derive expected seismic signal characteristics and use them to evaluate each of the seismic events. However, none of them can definitively be associated with this source. Atmospheric perturbations are generally expected to be generated during impacts; however, in this case, no signal could be identified as related to the known impact. Using scaling relationships based on the terrestrial and lunar analogs and numerical modeling, we predict the amplitude, peak frequency, and duration of the seismic signal that would have emanated from this impact. The predicted amplitude falls near the lowest levels of the measured seismometer noise for the predicted frequency. Hence, it is not surprising this impact event was not positively identified in the seismic data. Finding this crater was a lucky event as its formation this close to InSight has a probability of only ~0.2, and the odds of capturing it in before and after images are extremely low. We revisit impact-seismic discriminators in light of real experience with a seismometer on the Martian surface. Using measured noise of the instrument, we revise our previous prediction of seismic impact detections downward, from ~a few to tens, to just ~2 per Earth year, still with an order of magnitude uncertainty.

58 GEOSCIENCES↗

Two-Step Procedure to Detect Cosmological Gravitational Wave Backgrounds with Next-Generation Terrestrial Gravitational-Wave Detectors

Cosmological gravitational-wave backgrounds are an exciting science target for next-generation ground-based detectors, as they encode invaluable information about the primordial Universe. However, any such background is expected to be obscured by the astrophysical foreground from compact-binary coalescences. We propose a novel framework to detect a cosmological gravitational-wave background in the presence of binary black holes and binary neutron star signals with next-generation ground-based detectors, including Cosmic Explorer and the Einstein Telescope. Our procedure involves first removing all the individually resolved binary black hole signals by notching them out in the time-frequency domain. Then, we perform joint Bayesian inference on the individually resolved binary neutron star signals, the unresolved binary neutron star foreground, and the cosmological background. For a flat cosmological background, we find that we can claim detection at 5⁢𝜎 level when Ω ref ≥ 2.7 × 10 −12 /$\sqrt{𝑇_{obs}/yr}$, where 𝑇 obs is the observation time (in years), which is within a factor of ≲ 2 from the sensitivity reached in the absence of these astrophysical foregrounds.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Spin-1 Kitaev-Heisenberg model on a honeycomb lattice

We study the Kitaev-Heisenberg model with spin-1 local degree of freedom on a honeycomb lattice numerically by the infinite density matrix renormalization group method on a cylinder geometry. By tuning the relative value of the Kitaev and Heisenberg exchange couplings, we obtain the phase diagram with two spin liquid phases and four symmetry-broken phases. Here, we identify that the spin liquid phases are gapless by calculating the central charge at the pure Kitaev points without the Heisenberg interactions. Comparing to its spin-1/2 counterpart, the position and number of gapless modes of the spin-1 case are quite different. Due to the approximate Z 2 local conservations, the expectation value of the Wilson loop operator measuring the flux of each plaquette stays near to 1, and the static spin-spin correlations remain short range in the entire spin liquid phases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Average-atom model with Siegert states

In plasmas, electronic states can be well-localized bound states or itinerant free states, or something in between. In self-consistent treatments of plasma electronic structure such as the average-atom model, all states must be accurately resolved in order to achieve a converged numerical solution. Furthermore, this is a challenging numerical and algorithmic problem in large part due to the continuum of free states which is relatively expensive and difficult to resolve accurately. Siegert states are an appealing alternative. They form a complete eigenbasis with a purely discrete spectrum while still being equivalent to a representation in terms of the usual bound states and free states. However, many of their properties are unintuitive, and it is not obvious that they are suitable for self-consistent plasma electronic structure calculations. Here it is demonstrated that Siegert states can be used to accurately solve an average-atom model and offer advantages over the traditional finite-difference approach, including a concrete physical picture of pressure ionization and continuum resonances.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Extracting the Quantum Hall Conductance from a Single Bulk Wave Function

In this work, we propose a new formula that extracts the quantum Hall conductance from a single (2+1)⁢ D gapped wave function. The formula applies to general many-body systems that conserve particle number, and is based on the concept of modular flow, i.e., unitary dynamics generated from the entanglement structure of the wave function. The formula is shown to satisfy all formal properties of the Hall conductance: it is odd under time reversal and reflection, even under charge conjugation, universal and topologically rigid in the thermodynamic limit. Further evidence for relating the formula to the Hall conductance is obtained from conformal field theory arguments. Finally, we numerically check the formula by applying it to a noninteracting Chern band where excellent agreement is obtained.

2-dimensional systems↗

Data-Driven Learning for the Mori--Zwanzig Formalism: A Generalization of the Koopman Learning Framework

A theoretical framework which unifies the conventional Mori--Zwanzig formalism and the approximate Koopman learning of deterministic dynamical systems from noiseless observation is presented. In this framework, the Mori--Zwanzig formalism, developed in statistical mechanics to tackle the hard problem of construction of reduced-order dynamics for high-dimensional dynamical systems, can be considered as a natural generalization of the Koopman description of the dynamical system. We next show that, similar to the approximate Koopman learning methods, data-driven methods can be developed for the Mori--Zwanzig formalism with Mori's linear projection operator. We have developed two algorithms to extract the key operators, the Markov and the memory kernel, using time series of a reduced set of observables in a dynamical system. We have adopted the Lorenz `96 system as a test problem and solved for the above operators. These operators exhibit complex behaviors, which are unlikely to be captured by traditional modeling approaches in Mori--Zwanzig analysis. The nontrivial generalized fluctuation-dissipation relationship, which relates the memory kernel with the two-time correlation statistics of the orthogonal dynamics, was numerically verified as a validation of the solved operators. Here we present numerical evidence that the generalized Langevin equation, a key construct in the Mori--Zwanzig formalism, is more advantageous in predicting the evolution of the reduced set of observables than the conventional approximate Koopman operators.

97 MATHEMATICS AND COMPUTING↗

Fundamental absorption bandwidth to thickness limit for transparent homogeneous layers

Abstract Past work has considered the analytic properties of the reflection coefficient for a metal-backed slab. The primary result established a fundamental relationship for the minimal layer thickness to bandwidth ratio achievable for an absorber. There has yet to be establishment of a similar relationship for non-metal-backed layers, and here we present the universal result based on the Kramers–Kronig relations. Our theory is validated with transfer matrix calculations of homogeneous materials, and full-wave numerical simulations of electromagnetic metamaterials. Our results place more general fundamental limits on absorbers and thus will be important for both fundamental and applied studies.

Padilla, Willie J.↗

Leveraging machine learning to improve understanding and predictability of weather/climate extremes and the resilience of human systems

Focal Area(s): Developing hybrid models for predicting weather/climate extremes at the intersection of the urban environment (focal area #2). Science Challenge: Limited by many sources of uncertainty, numerical models are unable to accurately simulate weather and climate processes and this problem is particularly acute for modeling extreme events, including hurricanes. Moreover, computational limits within these numerical models bound our ability to develop large ensembles to explore this uncertainty and the implications relative to characterization of risk and resilience within the human systems that are impacted by changes in extremes.

54 ENVIRONMENTAL SCIENCES↗

On the Estimation of Circumbinary Orbital Properties

We describe a fast, approximate method to characterize the orbits of satellites around a central binary in numerical simulations. A goal is to distinguish the free eccentricity—the random motion of a satellite relative to a dynamically cool orbit—from oscillatory modes driven by the central binary’s time-varying gravitational potential. We assess the performance of the method using the Kepler-16, Kepler-47, and Pluto–Charon systems. We then apply the method to a simulation of orbital damping in a circumbinary environment, resolving relative speeds between small bodies that are slow enough to promote mergers and growth. These results illustrate how dynamical cooling can set the stage for the formation of Tatooine-like planets around stellar binaries and the small moons around the Pluto–Charon binary planet.

79 ASTRONOMY AND ASTROPHYSICS↗

Modeling and simulations of hydrodynamic shocks in a plasma flowing across randomized ICF scale laser beams

High-energy laser beams interacting with flowing plasmas can produce a plasma response that leads to deflection of the beam, beam bending. Such beams have usually a speckle structure generated by optical smoothing techniques that reduce the spatial and temporal coherence in the laser field pattern. The cumulative plasma response from laser speckles slows down the velocity of the incoming flow by momentum conservation. For slightly super-sonic flow the cumulative plasma response to the ponderomotive force exerted by the beam speckle ensemble is the strongest, such that slowing down the flow to subsonic velocities leads eventually to the generation of a shock around the cross section of the beam. This scenario has been predicted theoretically and is confirmed here by our hydrodynamic simulations in two dimensions with speckled beams and in one dimension with a reduced model. The conditions of shock generation are given in terms of the ponderomotive pressure, speckle size and the flow velocity. The nonlinear properties of the shocks are analyzed using Rankine–Hugoniot relations. According to linear theory, temporally smoothed laser beams exhibit a higher threshold for shock generation. Numerical simulations with beams that are smoothed by spectral dispersion compare well with the linear theory results, diverging from those produced by beams with only a random phase plates in the nonlinear regime. The conditions necessary for shock generation and their effects on the laser plasma coupling in inertial confinement fusion (ICF) experiments are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Decoding the formation of hammerhead ion populations observed by Parker Solar Probe

Context. In situ observations by the Parker Solar Probe (PSP) have revealed new properties of the proton velocity distributions (VDs), including hammerhead features that suggest a non-isotropic broadening of the beams. Aims. The present work proposes a very plausible explanation for the formation of hammerhead proton populations through the action of a proton firehose-like instability triggered by the proton beam. Methods. We investigated a self-generated firehose-like instability driven by the relative drift of ion populations using a simplified moment-based quasi-linear (QL) theory. While simpler and faster than advanced numerical simulations, this toy model provided rapid insights and concisely highlighted the role of plasma micro-instabilities in relaxing the observed anisotropies of particle VDs in the solar wind and space plasmas. Results. The QL theory proposed here shows that the resulting transverse waves are right-hand polarized and have two consequences on the protons: (i) They reduce the relative drift between the beam and the core, but above all, (ii) they induce a strong perpendicular temperature anisotropy specific to the observed hammerhead ion beam. Moreover, the long-run QL results suggest that these hammerhead distributions are rather transitory states that are still subject to relaxation mechanisms, in which instabilities such as the one discussed here are very likely involved. This foundational work motivates future detailed studies using advanced methods.

Shaaban, Shaaban M. (ORCID:000000030465598X)↗