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 109 records · Page 6

Post-Newtonian Initial Data with Waves: Progress in Evolution

"In Kelly \et al [Phys. Rev. D, 76:024008, (2007)], we presented new binary black-hole initial data adapted to puncture evolutions in numerical relativity. This data satisfies the constraint equations to 2.5 post-Newtonian order, and contains a transverse-traceless "wavy" metric contribution, violating the standard assumption of conformal flatness. We report on progress in evolving this data with a modem moving-puncture implementation of the BSSN equations in several numerical codes. We discuss the effect of the new metric terms on junk radiation and continuity of physical radiation extracted."

Kelly, Bernard J.↗

Numerical aperture limits on efficient ball lens coupling of laser diodes to single-mode fibers with defocus to balance spherical aberration

The potential capabilities and limitations of single ball lenses for coupling laser diode radiation to single-mode optical fibers have been analyzed; parameters important to optical communications were specifically considered. These parameters included coupling efficiency, effective numerical apertures, lens radius, lens refractive index, wavelength, magnification in imaging the laser diode on the fiber, and defocus to counterbalance spherical aberration of the lens. Limiting numerical apertures in object and image space were determined under the constraint that the lens perform to the Rayleigh criterion of 0.25-wavelength (Strehl ratio = 0.80). The spherical aberration-defocus balance to provide an optical path difference of 0.25 wavelength units was shown to define a constant coupling efficiency (i.e., 0.56). The relative numerical aperture capabilities of the ball lens were determined for a set of wavelengths and associated fiber-core diameters of particular interest for single-mode fiber-optic communication. The results support general continuing efforts in the optical fiber communications industry to improve coupling links within such systems with emphasis on manufacturing simplicity, system packaging flexibility, relaxation of assembly alignment tolerances, cost reduction of opto-electronic components and long term reliability and stability.

Wilson, R. Gale↗

On the Importance of Small Ice Crystals in Tropical Anvil Cirrus

In situ measurements of ice crystal concentrations and sizes made with aircraft instrumentation over the past two decades have often indicated the presence of numerous relatively small (< 50 m diameter) crystals in cirrus clouds. Further, these measurements frequently indicate that small crystals account for a large fraction of the extinction in cirrus clouds. The fact that the instruments used to make these measurements, such as the Forward Scattering Spectrometer Probe (FSSP) and the Cloud Aerosol Spectrometer (CAS), ingest ice crystals into the sample volume through inlets has led to suspicion that the indications of numerous small ]crystals could be artifacts of large ]crystal shattering on the instrument inlets. We present new aircraft measurements in anvil cirrus sampled during the Tropical Composition, Cloud, and Climate Coupling (TC4) campaign with the 2 ] Dimensional Stereo (2D ]S) probe, which detects particles as small as 10 m. The 2D ]S has detector "arms" instead of an inlet tube. Since the 2D ]S probe surfaces are much further from the sample volume than is the case for the instruments with inlets, it is expected that 2D ]S will be less susceptible to shattering artifacts. In addition, particle inter ]arrival times are used to identify and remove shattering artifacts that occur even with the 2D ]S probe. The number of shattering artifacts identified by the 2D ]S interarrival time analysis ranges from a negligible contribution to an order of magnitude or more enhancement in apparent ice concentration over the natural ice concentration, depending on the abundance of large crystals and the natural small ]crystal concentration. The 2D ]S measurements in tropical anvil cirrus suggest that natural small ]crystal concentrations are typically one to two orders of magnitude lower than those inferred from CAS. The strong correlation between the CAS/2D ]S ratio of small ]crystal concentrations and large ]crystal concentration suggests that the discrepancy is likely caused by shattering of large crystals on the CAS inlet. We argue that past measurements with CAS in cirrus with large crystals present may contain errors due to crystal shattering, and past conclusions derived from these measurements may need to be revisited. Further, we present correlations between CAS spurious concentration and 2D ]S large ]crystal mass from spatially uniform anvil cirrus sampling periods as an approximate guide for estimating quantitative impact of large ]crystal shattering on CAS concentrations in previous datasets. We use radiative transfer calculations to demonstrate that in the maritime anvil cirrus sampled during TC4, small crystals indicated by 2D ]S contribute relatively little cloud extinction, radiative forcing, or radiative heating in the anvils, regardless of anvil age or vertical location in the clouds. While 2D ]S ice concentrations in fresh anvil cirrus may often exceed 1 cm.3, and are observed to exceed 10 cm.3 in turrets, they are typically ~0.1 cm.3 and rarely exceed 1 cm.3 (<1.4% of the time) in aged anvil cirrus. We hypothesize that isolated occurrences of higher ice concentrations in aged anvil cirrus may be caused by ice nucleation driven by either small ]scale convection or gravity waves. It appears that the numerous small crystals detrained from convective updrafts do not persist in the anvil cirrus sampled during TC ]4.

Jensen, E. J.↗

An implementation of a high-order generalized finite difference method for solving the time-harmonic cold plasma wave equation in toroidal geometry

A high-order physics-informed meshless finite difference numerical technique is introduced for solving the time-harmonic cold plasma wave equation in toroidal geometries, presenting a novel application of the generalized finite difference (GFD) method to plasma wave simulations. The algorithm employs an irregular distribution of computational points, with local point density informed by the shortest wavelength derived from the cold plasma dispersion relation. Numerical stability and robustness are addressed using regularization techniques. The algorithm, implemented for two spatial dimensions, solves for the wave electric field and is demonstrated to achieve convergence rates of $\mathcal{O}$($\mathcal{h}$ $\mathcal{P}$ )⁠. Verification tests reproduce plane wave solutions, and example simulations of ion cyclotron resonance heating and electron cyclotron resonance heating demonstrate its capability, approaching realistic tokamak plasma scenarios. This work contributes to laying a foundation for the GFD method to be used in more sophisticated, optimized, and physically realistic full-wave simulations in time-harmonic plasma wave research.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Class of High-Resolution Explicit and Implicit Shock-Capturing Methods

The development of shock-capturing finite difference methods for hyperbolic conservation laws has been a rapidly growing area for the last decade. Many of the fundamental concepts, state-of-the-art developments and applications to fluid dynamics problems can only be found in meeting proceedings, scientific journals and internal reports. This paper attempts to give a unified and generalized formulation of a class of high-resolution, explicit and implicit shock capturing methods, and to illustrate their versatility in various steady and unsteady complex shock waves, perfect gases, equilibrium real gases and nonequilibrium flow computations. These numerical methods are formulated for the purpose of ease and efficient implementation into a practical computer code. The various constructions of high-resolution shock-capturing methods fall nicely into the present framework and a computer code can be implemented with the various methods as separate modules. Included is a systematic overview of the basic design principle of the various related numerical methods. Special emphasis will be on the construction of the basic nonlinear, spatially second and third-order schemes for nonlinear scalar hyperbolic conservation laws and the methods of extending these nonlinear scalar schemes to nonlinear systems via the approximate Riemann solvers and flux-vector splitting approaches. Generalization of these methods to efficiently include real gases and large systems of nonequilibrium flows will be discussed. Some perbolic conservation laws to problems containing stiff source terms and terms and shock waves are also included. The performance of some of these schemes is illustrated by numerical examples for one-, two- and three-dimensional gas-dynamics problems. The use of the Lax-Friedrichs numerical flux to obtain high-resolution shock-capturing schemes is generalized. This method can be extended to nonlinear systems of equations without the use of Riemann solvers or flux-vector splitting approaches and thus provides a large savings for multidimensional, equilibrium real gases and nonequilibrium flow computations.

Yee, H. C.↗

On the Susceptibility of Cold Tropical Cirrus to Ice Nuclei Abundance

Numerical simulations of cirrus formation in the tropical tropopause layer (TTL) are used to evaluate the impact heterogeneous ice nuclei (IN) abundance on cold cloud microphysical properties. The model includes homogeneous and heterogeneous ice nucleation, deposition growthsublimation, and sedimentation. Reanalysis temperature and wind fields are used to force the simulations, with addition of high-frequency wave temperature variability based on in situ measurements. The model results are constrained by recent in situ observations of TTL cirrus and relative humidity, as well as satellite measurements of cirrus bulk properties and occurrence frequencies. Temperature variability driven by high-frequency waves has a dominant influence on TTL cirrus microphysical properties and occurrence frequencies, and inclusion of these waves is required to produce agreement between the simulated and observed cloud properties as well as the abundance of TTL cirrus. With a composition-independent supersaturation threshold for homogeneous freezing of aqueous aerosols, the model produces excessive ice concentrations compared with in situ observations. Inclusion of relatively numerous heterogeneous nuclei ( 100L) in the simulations improves the agreement with observed ice concentrations. However, when IN contribute significantly to TTL cirrus ice nucleation, the occurrence frequency of large supersaturations with respect to ice is less than indicated by in situ measurements. The simulated TTL cirrus extinction statistics agree with observations (within uncertainties) from both in situ and remote-sensing measurements. We find that the sensitivity of TTL cirrus extinction and ice water content statistics to heterogeneous ice nuclei abundance is relatively weak. Likewise, the simulated occurrence frequencies of TTL cirrus are insensitive to ice nuclei abundance, both in terms of cloud frequency height distribution and regional distribution throughout the tropics.

occurences↗

Observing Mergers of Non-Spinning Black-Hole Binaries

Advances in the field of numerical relativity now make it possible to calculate the final, most powerful merger phase of binary black-hole coalescence for generic binaries. The state of the art has advanced well beyond the equal-mass case into the unequal-mass and spinning regions of parameter space. We present a study of the nonspinning portion of parameter space, primarily using an analytic waveform model tuned to available numerical data, with an emphasis on observational implications. We investigate the impact of varied m8BS ratio on merger signal-to-noise ratios (SNR~) for several detectors, and compare our results with expectations from the test-mass limit. We note a striking similarity of the waveform phasing of the merger waveform across the available mass ratios. Motivated by this, we calculate the match between our equal-mass and 4:1 mass-ratio waveforms during the merger as a function of location on the source sky, using a new formalism for the match that accounts for higher harmonics. This is an indicator of the amount of degeneracy in mass ratio for mergers of moderate mass ratio systems.

McWilliams, Sean T.↗

Modeling and Analysis of Wrinkled Membranes: An Overview

Thin-film membranes are basic elements of a variety of space inflatable/deployable structures. Wrinkling degrades the performance and reliability of these membrane structures, and hence has been a topic of continued interest. Wrinkling analysis of membranes for general geometry and arbitrary boundary conditions is quite challenging. The objective of this presentation is two-fold. Firstly, the existing models of wrinkled membranes and related numerical solution methods are reviewed. The important issues to be discussed are the capability of a membrane model to characterize taut, wrinkled and slack states of membranes in a consistent and physically reasonable manner; the ability of a wrinkling analysis method to predict the formation and growth of wrinkled regions, and to determine out-of-plane deformation and wrinkled waves; the convergence of a numerical solution method for wrinkling analysis; and the compatibility of a wrinkling analysis with general-purpose finite element codes. According to this review, several opening issues in modeling and analysis of wrinkled membranes that are to be addressed in future research are summarized, The second objective of this presentation is to discuss a newly developed membrane model of two viable parameters (2-VP model) and associated parametric finite element method (PFEM) for wrinkling analysis are introduced. The innovations and advantages of the proposed membrane model and PFEM-based wrinkling analysis are: (1) Via a unified stress-strain relation; the 2-VP model treat the taut, wrinkled, and slack states of membranes consistently; (2) The PFEM-based wrinkling analysis has guaranteed convergence; (3) The 2-VP model along with PFEM is capable of predicting membrane out-of-plane deformations; and (4) The PFEM can be integrated into any existing finite element code. Preliminary numerical examples are also included in this presentation to demonstrate the 2-VP model and PFEM-based wrinkling analysis approach.

Yang, B.↗

The Stability of Radiatively Cooling Jets I. Linear Analysis

The results of a spatial stability analysis of a two-dimensional slab jet, in which optically thin radiative cooling is dynamically important, are presented. We study both magnetized and unmagnetized jets at external Mach numbers of 5 and 20. We model the cooling rate by using two different cooling curves: one appropriate to interstellar gas, and the other to photoionized gas of reduced metallicity. Thus, our results will be applicable to both protostellar (Herbig-Haro) jets and optical jets from active galactic nuclei. We present analytical solutions to the dispersion relations in useful limits and solve the dispersion relations numerically over a broad range of perturbation frequencies. We find that the growth rates and wavelengths of the unstable Kelvin-Helmholtz (K-H) modes are significantly different from the adiabatic limit, and that the form of the cooling function strongly affects the results. In particular, if the cooling curve is a steep function of temperature in the neighborhood of the equilibrium state, then the growth of K-H modes is reduced relative to the adiabatic jet. On the other hand, if the cooling curve is a shallow function of temperature, then the growth of K-H modes can be enhanced relative to the adiabatic jet by the increase in cooling relative to heating in overdense regions. Inclusion of a dynamically important magnetic field does not strongly modify the important differences between an adiabatic jet and a cooling jet, provided the jet is highly supermagnetosonic and not magnetic pressure-dominated. In the latter case, the unstable modes behave more like the transmagnetosonic magnetic pressure-dominated adiabatic limit. We also plot fluid displacement surfaces associated with the various waves in a cooling jet in order to predict the structures that might arise in the nonlinear regime. This analysis predicts that low-frequency surface waves and the lowest order body modes will be the most effective at producing observable features in the jet.

Hardee, Philip E.↗

Performance-portable Binary Neutron Star Mergers with AthenaK

We introduce an extension to the AthenaK code for general-relativistic magnetohydrodynamics (GRMHD) in dynamical spacetimes using a 3+1 conservative Eulerian formulation. Like the fixed-spacetime GRMHD solver, we use standard finite-volume methods to evolve the fluid and a constrained-transport scheme to preserve the divergence-free constraint for the magnetic field. We also utilize a first-order flux correction (FOFC) scheme to reduce the need for an artificial atmosphere and optionally enforce a maximum principle to improve robustness. We demonstrate the accuracy of AthenaK using a set of standard tests in flat and curved spacetimes. Using a SANE accretion disk around a Kerr black hole, we compare the new solver to the existing solver for stationary spacetimes using the so-called "HARM-like" formulation. We find that both formulations converge to similar results. We also include the first published binary neutron star (BNS) mergers performed on graphical processing units (GPUs). Thanks to the FOFC scheme, our BNS mergers maintain a relative error of $\mathcal{O}$(10 –11 ) or better in baryon mass conservation up to collapse. Finally, we perform scaling tests of AthenaK on OLCF Frontier, where we show excellent weak scaling of ≥80% efficiency up to 32,768 GPUs and 74% up to 65,536 GPUs for a GRMHD problem in dynamical spacetimes with six levels of mesh refinement. AthenaK achieves an order-of-magnitude speedup using GPUs compared to CPUs, demonstrating that it is suitable for performing numerical relativity problems on modern exascale resources.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Smoothing and flattening the universe through slow contraction versus inflation

In a systematic study, we use an equivalent pair of improved numerical relativity codes based on a tetrad-formulation of the classical Einstein-scalar field equations to examine whether slow contraction or inflation (or both) can resolve the homogeneity, isotropy and flatness problems. Our finding, based on a set of gauge/frame invariant diagnostics and the models considered, is that slow contraction robustly and rapidly smooths and flattens spacetime beginning from initial conditions that are outside the perturbative regime of the flat Friedmann-Robertson-Walker metric, whereas inflation fails these tests. We present new numerical evidence supporting the conjecture that the combination of ultralocal evolution and an effective equation-of-state with pressure much greater than energy density is the key to having robust and rapid smoothing. The opposite of ultralocality occurs in expanding spacetimes, which is the leading obstruction to smoothing following a big bang.

79 ASTRONOMY AND ASTROPHYSICS↗

Binary Black Hole Mergers and Recoil Kicks

Recent developments in numerical relativity have made it possible to follow reliably the coalescence of two black holes from near the innermost stable circular orbit to final ringdown. This opens up a wide variety of exciting astrophysical applications of these simulations. Chief among these is the net kick received when two unequal mass or spinning black holes merge. The magnitude of this kick has bearing on the production and growth of supermassive black holes during the epoch of structure formation, and on the retention of black holes in stellar clusters. Here we report the first accurate numerical calculation of this kick, for two nonspinning black holes in a 1.5:1 mass ratio, which is expected based on analytic considerations to give a significant fraction of the maximum possible recoil. We have performed multiple runs with different initial separations, orbital angular momenta, resolutions, extraction radii, and gauges. The full range of our kick speeds is 86-116 kilometers per second, and the most reliable runs give kicks between 86 and 97 kilometers per second. This is intermediate between the estimates from two recent post-Newtonian analyses and suggests that at redshifts z greater than 10, halos with masses less than 10(exp 9) M(sub SUN) will have difficulty retaining coalesced black holes after major mergers.

Centrella, Joan↗

Comment on "An Analysis of VLF Electric Field Spectra Measured in Titan's Atmosphere by The Huygens Probe" By J. A. Morente et al.

Morente et al. have recently revisited the VLF electric field measurements made with the Permittivity, Wave and Altimetry (PWA) instrument during the descent of the Huygens Probe through the atmosphere of Titan. They assert that they have identified several harmonics of the transverse resonance mode of the surface?]ionosphere cavity, which would prove the existence of an electrical activity in the atmosphere of the largest satellite of Saturn. We refute this finding on the basis that it results from an artifact due to an improper analysis of the data set. [2] The investigators of the Permittivity, Wave and Altimetry (PWA) experiment on the Huygens Probe have reported the extremely low frequency (ELF) and very low frequency (VLF) electric signals recorded during the descent through the atmosphere of Titan. The PWA data are archived in the Planetary Science Archive (PSA) of ESA, and an extensive description of the instrument is at the disposal of the scientific community. Morente and his coworkers have revisited this data set and reported the results of their investigations in two papers. In a first paper, they claim that they have detected in the ELF range (0.100 Hz) several harmonics of a global resonance allegedly generated by lightning activity in the spherical cavity guide formed by the surface of Titan and the inner boundary of the ionosphere, a phenomenon similar to the Schumann resonance observed at EartH In the second paper dedicated to the VLF electric signal recorded by PWA, in the range 0.10 kHz, they argue that they can also bring out the transverse resonance and its harmonics, a more local phenomenon that develops around the excitation source and whose frequency is controlled by the separation between Titan?fs surface and the inner ionospheric boundary. [3] The PWA investigators have analyzed the narrowband ELF signal at about 36 Hz effectively observed during the entire descent. They have not endorsed, however, the alternative approach of Morente et al. that discloses additional ELF narrow-band signals and numerous related harmonics. Hamelin et al. and Beghin et al. have questioned the work of Morente and his coworkers and demonstrated unambiguously that their findings are mere artifacts entirely due to a faulty procedure. Morente et al. refuted this critical analysis of their ELF paper in a rebuttal. The present comment draws a parallelism between the ELF and VLF papers and demonstrates that the VLF signal carries no evidence whatever of any transverse resonance. The analyses of the ELF and VLF signals by Morente and his coworkers are flawed, and their conclusions are shown here to be invalid. A more comprehensive and thorough assessment of the numerical approaches proposed by Morente et al. is given by Berthelin et al..

Grard, Rejean↗

Modelling populations of kilonovae

Abstract The 2017 detection of a kilonova coincident with gravitational-wave emission has identified neutron star mergers as the major source of the heaviest elements and dramatically constrained alternative theories of gravity. Observing a population of such sources has the potential to transform cosmology, nuclear physics, and astrophysics. However, with only one confident multi-messenger detection currently available, modelling the diversity of signals expected from such a population requires improved theoretical understanding. In particular, models that are quick to evaluate and are calibrated with more detailed multi-physics simulations are needed to design observational strategies for kilonovae detection and to obtain rapid-response interpretations of new observations. We use grey-opacity models to construct populations of kilonovae, spanning ejecta parameters predicted by numerical simulations. Our modelling focuses on wavelengths relevant for upcoming optical surveys, such as the Rubin Observatory Legacy Survey of Space and Time (LSST). In these simulations, we implement heating rates that are based on nuclear reaction network calculations. We create a Gaussian-process emulator for kilonova grey opacities, calibrated with detailed radiative transfer simulations. Using recent fits to numerical relativity simulations, we predict how the ejecta parameters from binary neutron star (BNS) mergers shape the population of kilonovae, accounting for the viewing-angle dependence. Our simulated population of BNS mergers produce peak i-band absolute magnitudes of −20 ≤ Mi ≤ −11. A comparison with detailed radiative transfer calculations indicates that further improvements are needed to accurately reproduce spectral shapes over the full light curve evolution.

79 ASTRONOMY AND ASTROPHYSICS↗

Radiation hydrodynamics modelling of kilonovae with SNEC

ABSTRACT We develop a method to compute synthetic kilonova light curves that combine numerical relativity simulations of neutron star mergers and the SNEC radiation–hydrodynamics code. We describe our implementation of initial and boundary conditions, r-process heating, and opacities for kilonova simulations. We validate our approach by carefully checking that energy conservation is satisfied and by comparing the SNEC results with those of two semi-analytic light-curve models. We apply our code to the calculation of colour light curves for three binaries having different mass ratios (equal and unequal mass) and different merger outcome (short-lived and long-lived remnants). We study the sensitivity of our results to hydrodynamic effects, nuclear physics uncertainties in the heating rates, and duration of the merger simulations. We find that hydrodynamics effects are typically negligible and that homologous expansion is a good approximation in most cases. However, pressure forces can amplify the impact of uncertainties in the radioactive heating rates. We also study the impact of shocks possibly launched into the outflows by a relativistic jet. None of our models match AT2017gfo, the kilonova in GW170817. This points to possible deficiencies in our merger simulations and kilonova models that neglect non-LTE effects and possible additional energy injection from the merger remnant and to the need to go beyond the assumption of spherical symmetry adopted in this work.

79 ASTRONOMY AND ASTROPHYSICS↗

Impact of in situ nuclear networks and atomic opacities on neutron star merger ejecta dynamics, nucleosynthesis, and kilonovae

Context. Binary neutron star merger (BNSM) ejecta are key sites of rapid neutron capture (r-process) nucleosynthesis and they produce kilonovae powered by the radioactive decay of freshly synthesized nuclei. Modeling their evolution requires multi-physics simulations involving hydrodynamics, nuclear reactions, and radiative processes. The impact of nuclear burning and atomic opacity is poorly understood and often treated with simplified prescriptions. Aims. We systematically investigate different treatments of nuclear heating, particle thermalization, and atomic opacities in radiation-hydrodynamics simulations of BNSM ejecta and kilonova light curves. Methods. Ejecta profiles from long-term numerical-relativity simulations of asymmetric neutron star binaries with a massive neutron star remnant were evolved to ∼30 days using a 2D ray-by-ray approach. We compared simplified heating-rate and thermalization prescriptions with in situ Nuclear reaction Network (NN) calculations that track nuclear energy deposition and include a composition-dependent thermalization scheme. We also contrasted various gray opacity models with a frequency-dependent treatment based on atomic calculations. Results. Coupling NN and hydrodynamics significantly affects nucleosynthesis and kilonova emission. Assuming homologous expansion alters abundance evolution and produces a narrower, less populated second r-process peak and a third peak shifted to higher mass numbers. The back-reaction of nuclear heating affects the temperature evolution enough to delay and redden the early (t∼ hours) kilonova peaks. A constant thermalization efficiency underestimates and reddens the early emission while overestimating the late-time luminosity compared to the composition-dependent treatment. Analytical opacity prescriptions yield a more extended, colder photosphere, resulting in dimmer, redder kilonovae at early times (t≲ hour), while the delayed recession of the photosphere prolongs the red emission at t ≳ 5 days. Conclusions. Coupling hydrodynamics to an in situ NN is crucial for reliable nucleosynthesis and kilonova predictions. Resolving the first several hundred milliseconds of the hydrodynamics is essential for robust nucleosynthesis calculations. Composition-dependent thermalization and frequency-dependent, atomic-physics-based opacities are needed to accurately capture the temperature evolution of the ejecta and the brightness and color evolution of the kilonova. Calibrated analytic nuclear-power fits with simplified thermalization and opacity prescriptions can still reproduce the density and temperature evolution of the ejecta.

74 ATOMIC AND MOLECULAR PHYSICS↗

Building bulk geometry from the tensor Radon transform

Using the tensor Radon transform and related numerical methods, we study how bulk geometries can be explicitly reconstructed from boundary entanglement entropies in the specific case of AdS 3 /CFT 2 . We find that, given the boundary entanglement entropies of a 2d CFT, this framework provides a quantitative measure that detects whether the bulk dual is geometric in the perturbative (near AdS) limit. In the case where a well-defined bulk geometry exists, we explicitly reconstruct the unique bulk metric tensor once a gauge choice is made. We then examine the emergent bulk geometries for static and dynamical scenarios in holography and in many-body systems. Apart from the physics results, our work demonstrates that numerical methods are feasible and effective in the study of bulk reconstruction in AdS/CFT.

79 ASTRONOMY AND ASTROPHYSICS↗

An investigation of the multi-mode Richtmyer-Meshkov instability at a gas/HE interface using Pagosa

Here, in this work, we present a hydrocode Pagosa and explore the Richtmyer-Meshkov Instability (RMI) at an air/high explosive (HE) interface for the first time that is important but has not received much attention yet in the high explosive safety field. Thus, the presented Pagosa can be expected to predict the whole deflagration-to-detonation transition (DDT) process in future. In Pagosa, spatial discretization is implemented on cubic staggered grids by computing different variables at the vertex and the cell center, respectively, a special operator-splitting technique is employed to reduce the computational cost, and an artificial viscosity is added to handle the discontinuous shock waves in our simulations. To quantitatively evaluate the capability of Pagosa to solve these kinds of instabilities, the single mode Rayleigh-Taylor instability (RTI) and the multimode RMI at an air/SF 6 interface are first performed, respectively. The Pagosa results are compared with the related numerical solutions in the existing references and the experimental result. Moreover, a theoretical derivation of growth of RTI is also provided based on our numerical method. Subsequently, we explore the multi-mode RMI at an air/HE interface as well as the effects of several factors using Pagosa. Numerical results show that Pagosa is a powerful toolset to generate the right structures and the amplitude of RTI and RMI at an air/SF 6 interface. The solid HE can be penetrated by a strong shock wave and forms RMI deformations. The RMI at an air/HE interface behaves very different than at an air/SF 6 interface, periodic, decreased oscillation is observed due to material character, and is very sensitive to the initial simulation settings, that is, a tiny change in physical quantities will lead to a remarkable RMI structure, which is also observed in a shock bubble interaction. The findings in this work are significant, and will present a new insight for the high explosive field.

97 MATHEMATICS AND COMPUTING↗